Threshold changes in rainfall‐runoff generation commonly represent shifts in runoff mechanisms and hydrologic connectivity controlling water and solute transport and transformation. In watersheds with limited human influence, threshold runoff responses reflect interaction between precipitation event and antecedent soil moisture. Similar analyses are lacking in intensively managed landscapes where installation of subsurface drainage tiles has altered connectivity between the land surface, groundwater, and streams, and where application of fertilizer has created significant stores of subsurface nitrogen. In this study, we identify threshold patterns of tile‐runoff generation for a drained agricultural field in Illinois and evaluate how antecedent conditions—including shallow soil moisture, groundwater table depth, and the presence or absence of crops—control tile response. We relate tile‐runoff thresholds to patterns of event nitrate load observed across multiple storm events and evaluate how antecedent conditions control within‐event nitrate concentration‐discharge relationships. Our results demonstrate that an event tile‐runoff threshold emerges relative to the sum of gross precipitation and indices of antecedent shallow soil moisture and antecedent below‐tile groundwater moisture deficit, indicating that both shallow soil and below‐tile storages must be filled to generate significant runoff. In turn, event nitrate load shows a linear dependence on runoff for most time periods, suggesting that subsurface nitrate export and storage can be estimated using runoff threshold relationships and long‐term average nitrate concentrations. Finally, within‐event nitrate concentration‐discharge relationships are controlled by event size and the antecedent tile flow state because these factors dictate the sequence of flow path activation and tile connectivity over a storm event.

Core Ideas: IML‐CZO is structured to study system responses through event‐based monitoring. Management legacy has shaped critical zone processes. Management and weather affect landscape heterogeneity and surface–subsurface pathways. In intensively managed landscapes, interactions between surface (tillage) and subsurface (tile drainage) management with prevailing climate/weather alter landscape characteristics, transport pathways, and transformation rates of surface/subsurface water, soil/sediment, and particulate/dissolved nutrients. To capture the high spatial and temporal variability of constituent transport and residence times in the critical zone (between the bedrock and canopy) of these altered landscapes, both storm event and continuous measurements are needed. The Intensively Managed Landscapes Critical Zone Observatory (IML‐CZO) is comprised of three highly characterized, well instrumented, and representative watersheds (i.e., Clear Creek, Iowa; Upper Sangamon River, Illinois; and Minnesota River, Minnesota). It is organized to quantify the heterogeneity in structure and dynamic response of critical zone processes to human activities in the context of the glacial and management (anthropogenic) legacies. Observations of water, sediment, and nutrients are made at nested points of the landscape in the vertical and lateral directions during and between storm events (i.e., continuously). The measurements and corresponding observational strategy are organized as follows. First, reference measurements from surface soil and deep core extractions, geophysical surveys, lidar, and hyperspectral data, which are common across all Critical Zone Observatories, are available. The reference measurements include continuous quantification of energy, water, solutes, and sediment fluxes. The reference measurements are complemented with event‐based measurements unique to IML‐CZO. These measurements include water table fluctuations, enrichment ratios, and roughness as well as bank erosion, hysteresis, sediment sources, and lake/floodplain sedimentation. The coupling of reference and event‐based measurements support testing of the central hypothesis (i.e., system shifts from transformer to transporter in IML‐CZO due to the interplay between management and weather/climate). Data collected since 2014 are available through a data repository and through the Geodashboard interface, which can be used for process‐based model simulations.

The multifaceted interface of plant roots, microbes, water, and soil can be considered a critical zone within the Critical Zone as it is host to many important dynamically linked processes, including the promotion of nutrient cycling through absorption and rhizodeposition, interaction and feedbacks with microorganisms and fungi, root‐facilitated hydraulic redistribution, and soil carbon dynamics. Such important processes in the Critical Zone have not been fully characterized and modeled in an ecohydrologic framework linking above‐ground natural and/or anthropogenic processes to below‐ground biogeochemical cycling. Specifically, the relation between root exudates and nutrient cycling remains an open challenge. Here we present the model REWT (Root Exudation in Watershed‐scale Transport) to demonstrate the systematic modeling of root exudation in an interconnected ecohydrologic framework. REWT incorporates an explicit dynamic root exudation transport model, nutrient absorption, and coupled microbial processes within the framework of a validated ecohydrologic model. Model simulations demonstrate the influence of root exudation of glucose, a polysaccharide that serves as fuel for microbes, and flavonoids, which can act as a biological nitrification inhibitor on microbial processes linked to soil carbon and nitrogen cycling. To demonstrate the capabilities of this theoretical framework, we parameterize REWT for corn and soybean crops in the Midwestern United States, and simulations indicate that rates of nitrification and respiration were substantially altered compared to model simulations in which root exudation was not explicitly included. This work demonstrates the importance of systematically incorporating root exudates into hydrobiogeochemical models and can serve to inform experimental design for active root zone processes.

We explore the impacts of tile drains in agricultural fields on the coupled age and concentration dynamics of nitrate, immobile ammonium, mobile ammonia and ammonium, and nonreactive tracers such as chloride. We implement two mobile interacting pore domains to capture matrix and preferential flow paths in a coupled ecohydrology and biogeochemistry model, Dhara. We apply this model to an agricultural farm that utilizes a corn‐soybean rotation in the Midwestern United States located in the Intensively Managed Landscapes Critical Zone Observatory. In general, we observe both low concentration and age of nitrate in the areas that are classified as topographic depressions even with the presence of tile drains. Also, an increase in the age of mobile ammonia/ammonium is observed after installing tile drains. This is in contrast to the cases for nitrate, immobile ammonium, and nonreactive tracer. These results arise because the depletion of mobile ammonia/ammonium due to tile drainage causes a high mobility flux from immobile ammonium to mobile ammonia/ammonium, which also carries a considerable amount of relatively old age of nitrogen from immobile ammonium to mobile ammonia/ammonium. These results illustrate how storm event scale dynamics impact spatial heterogeneity and temporal variability of the efflux, which helps in disentangling the complexity of nitrogen dynamics in the soil. This understanding can contribute to precision agriculture for nitrogen applications to reduce environmental impacts.