Tile Drainage Density Reduces Groundwater Travel Times and Compromises Riparian Buffer Effectiveness

Improving model capability in simulating spatiotemporal variations and flow contributions of nitrate export in tile-drained catchments

It is essential to identify the dominant flow paths, hot spots and hot periods of hydrological nitrate-nitrogen (NO3-N) losses for developing nitrogen loads reduction strategies in agricultural watersheds. Coupled biogeochemical transformations and hydrological connectivity regulate the spatiotemporal dynamics of water and NO3-N export along surface and subsurface flows. However, modeling performance is usually limited by the oversimplification of natural and human-managed processes and insufficient representation of spatiotemporally varied hydrological and biogeochemical cycles in agricultural watersheds. In this study, we improved a spatially distributed process-based hydro-ecological model (DLEM-catchment) and applied the model to four tile-drained catchments with mixed agricultural management and diverse landscape in Iowa, Midwestern US. The quantitative statistics show that the improved model well reproduced the daily and monthly water discharge, NO3-N concentration and loading measured from 2015 to 2019 in all four catchments. The model estimation shows that subsurface flow (tile flow + lateral flow) dominates the discharge (70-75%) and NO3-N loading (77-82%) over the years. However, the contributions of tile drainage and lateral flow vary remarkably among catchments due to different tile-drained area percentages and the presence of farmed potholes (former depressional wetlands that have been drained for agricultural production). Furthermore, we found that agricultural management (e.g. tillage and fertilizer management) and catchment characteristics (e.g. soil properties, farmed potholes, and tile drainage) play important roles in predicting the spatial distributions of NO3-N leaching and loading. The simulated results reveal that the model improvements in representing water retention capacity (snow processes, soil roughness, and farmed potholes) and tile drainage improved model performance in estimating discharge and NO3-N export at a daily time step, while improvement of agricultural management mainly impacts NO3-N export prediction. This study underlines the necessity of characterizing catchment properties, agricultural management practices, flow-specific NO3-N movement, and spatial heterogeneity of NO3-N fluxes for accurately simulating water quality dynamics and predicting the impacts of agricultural conservation nutrient reduction strategies.

Characterizing Differences in Sources of and Contributions to Fecal Contamination of Sediment and Surface Water with the Microbial FIT Framework

Surface water monitoring and microbial source tracking (MST) are used to identify host sources of fecal pollution and protect public health. However, knowledge of the locations of spatial sources and their relative impacts on the environment is needed to effectively mitigate health risks. Additionally, sediment samples may offer time-integrated information compared to transient surface water. Thus, we implemented the newly developed microbial find, inform, and test framework to identify spatial sources and their impacts on human (HuBac) and bovine (BoBac) MST markers, quantified from both riverbed sediment and surface water in a bovine-dense region. Dairy feeding operations and low-intensity developed land-cover were associated with 99% (p-value < 0.05) and 108% (p-value < 0.05) increases, respectively, in the relative abundance of BoBac in sediment, and with 79% (p-value < 0.05) and 39% increases in surface water. Septic systems were associated with a 48% increase in the relative abundance of HuBac in sediment and a 56% increase in surface water. Stronger source signals were observed for sediment responses compared to water. By defining source locations, predicting river impacts, and estimating source influence ranges in a Great Lakes region, this work informs pollution mitigation strategies of local and global significance.

Quantifying the effectiveness of a saturated buffer to reduce tile NO3-N concentrations in eastern Iowa

Agricultural drainage tiles are primary contributors to NO₃-N export from Iowa croplands. Saturated buffers are a relatively new conservation practice that diverts tile water into a distribution tile installed in a riparian buffer parallel to a stream with the intent of enhancing NO₃-N processing within the buffer. In this study, tile NO₃-N concentration reductions were characterized through two different saturated buffers at a working farm site in eastern Iowa. Study objectives were to (1) evaluate the hydrogeology and water quality patterns in the saturated buffer and (2) quantify the reduction in tile NO₃-N concentration from the saturated buffer installation. Results showed that the two saturated buffers are reducing NO₃-N concentrations in tile drainage water from input concentrations of approximately 15 mg/l to levels < 1.5 mg/l at the streamside well locations. The reduction occurs rapidly in the fine-textured and organic-rich alluvial soils with most of the reduction occurring within 1.5 m of the distribution line. Denitrification is hypothesized as being primarily responsible for the concentration reductions based on soil and water chemistry conditions, completion of a geophysical survey (quantifying low potential for N loss to deeper aquifers), and comparisons to other similar Iowa sites. The study provides more assurance to new adopters that this practice can be installed in many areas throughout the Midwestern Cornbelt region.

Controls on subsurface nitrate and dissolved reactive phosphorus losses from agricultural fields during precipitation-driven events

The magnitude of nitrogen (N) and phosphorus (P) exported from agricultural fields via subsurface tile drainage systems is determined by site-specific interactions between weather, soil, field, and management characteristics. Here, we used multiple regression analyses to evaluate the influence of 29 controls of precipitation event-driven discharge, nitrate (NO₃^--N) load, and dissolved reactive P (DRP) load from subsurface tile drains, leveraging a unique dataset of ~7000 precipitation events observed across 40 agricultural fields (n = 190 site years) instrumented to collect continuous water quality samples. We calculated marginal effects of significant controls and assessed the modifying influence of event rainfall, duration, and intensity, and antecedent precipitation. Tile discharge was strongly and positively influenced by previous 7-day precipitation and total rainfall and negatively influenced by daily temperature and tile spacing. Both tile NO₃^--N and DRP loads were positively influenced by transport and source variables, including event discharge and total fertilizer applied as well as soil test P (STP) in the case of tile DRP load; factors with the strongest negative influence on tile NO₃^--N and DRP loads were related to time of year. The strength and direction of both positive and negative controls also varied with precipitation characteristics. For example, the positive influence of event discharge on nutrient loads lessened as event duration, event intensity, and previous 7-day precipitation increased, while the positive influence of N and P sources strengthened, particularly in response to extreme (or maximum) events. Results here demonstrate the predominant role of transport and source controls while accounting for interactive effects among site-specific characteristics and underscore the importance of storm dynamics when managing N and P loss from agricultural fields.

Drivers of excess phosphorus and stream sediments in a nested agricultural catchment during base and stormflow conditions

A variety of landscape and hydrological characteristics influence nutrient concentrations and suspended sediments in freshwater systems, yet the combined influence of these characteristics within nested agricultural catchments is still poorly understood, particularly across varying flow states. To tease apart potential drivers at within-catchment scales, it is necessary to sample at a spatiotemporal resolution that captures how landscape drivers change with time. The overall objective of this study was to evaluate the relative influence of landscape and hydrological characteristics at sub-catchment scales in relation to total P (TP), soluble reactive P (SRP), the ratio of SRP and TP (SRP/TP), and total suspended solids (TSS) across varying flow conditions. Synoptic surveys were conducted at 13 longitudinal sampling sites under a variety of flow conditions (n = 14) between 2016 and 2017 in the Innisfil Creek watershed, southern Ontario. The surveys were grouped into baseflow and stormflow conditions, and partial least squares regression (PLSR) was used to characterize the relationships between catchment characteristics, median concentrations of P, and TSS. Soil texture (i.e., clay dominated), winter wheat (Triticum aestivum L.), and constructed drain density had the largest influences on stormflow SRP and SRP/TP ratios, but measures of soil erosion, like the Bank Erosion Hazard Index and sinuosity, had the largest influence on stormflow TSS. During baseflow periods, these landscape characteristics were not informative, and they were difficult to tie to in-stream conditions. Overall, our PLSR models indicated that buried tile drainage was a major source of SRP in Innisfil Creek, whereas bank erosion was a dominant source of TSS.

Quantifying the contribution of tile drainage to basin-scale water yield using analytical and numerical models

The Des Moines Lobe (DML) of north-central Iowa has been artificially drained by subsurface drains and surface ditches to provide some of the most productive agricultural land in the world. Herein we report on the use of end-member mixing analysis (EMMA) models and the numerical model Soil and Water Assessment Tool (SWAT) to quantify the contribution of tile drainage to basin-scale water yields at various scales within the 2370 km² Boone River watershed (BRW), a subbasin within the Des Moines River watershed. EMMA and SWAT methods suggested that tile drainage provided approximately 46 to 54% of annual discharge in the Boone River and during the March to June period, accounted for a majority of flow in the river. In the BRW subbasin of Lyons Creek, approximately 66% of the annual flow was sourced from tile drainage. Within the DML region, tile drainage contributes to basin-scale water yields at scales ranging from 40 to 16,000 km², with downstream effects diminishing with increasing watershed size. Developing a better understanding of water sources contributing to river discharge is needed if mitigation and control strategies are going to be successfully targeted to reduce downstream nutrient export.

In Situ Denitrification in Saturated Riparian Buffers

Excess NO leaching from the agricultural Midwest via tile drainage water has contributed to both local drinking water and national Gulf of Mexico benthic hypoxia concerns. Both in-field and edge-of-field practices have been designed to help mitigate NO flux to surface waters. Edge-of-field practices focus on maximizing microbial denitrification, the conversion of NO to N gas. This study assessed denitrification rates from two saturated riparian buffers (SRBs) for 2 yr and a third SRB for 1 yr, for a total of five sample years. These SRBs were created by diverting NO-rich tile drainage water into riparian buffers soils. The SRBs in this study removed between 27 and 96% of the total diverted NO load. Measured cumulative average denitrification rate for each SRB sample year accounted for between 3.7 and 77.3% of the total NO removed. Both the cumulative maximum and 90% confidence interval denitrification rates accounted for all of the NO removed by the SRBs in three of the five sample years, indicating that denitrification can be a dominant NO removal mechanism in this edge-of-field practice. When adding the top 20 cm of each core to the cumulative denitrification rates for each SRB, denitrification accounted for between 33 and over 100% of the total NO removed. Buffer age (time since establishment) was speculated to enhance denitrification rates, and there was a trend of the soil closer to the surface making up the majority of the total denitrification rate. Finally, both NO and C could limit denitrification in these SRBs.

Performance of Saturated Riparian Buffers in Iowa, USA

Nitrate from artificial drainage pipes (tiles) underlying agricultural fields is a major source of reactive N, especially NO, in surface waters. A novel approach for reducing NO loss is to intercept a field tile where it crosses a riparian buffer and divert a fraction of the flow as shallow groundwater within the buffer. This practice is called a saturated riparian buffer (SRB), and although it is promising, little data on the performance of the practice is available. This research investigated the effectiveness of SRBs in removing NO at six sites installed across Iowa, resulting in a total of 17 site-years. Water flow and NO in the tile outlets, diverted into the buffers, and NO concentration changes within the buffers were monitored throughout the year at each site. Results showed that all the SRBs were effective in removing NO from the tile outlet, with the average annual NO load removal ranging from 13 to 179 kg N for drainage areas ranging from 3.4 to 40.5 ha. This is NO that would have otherwise discharged directly into the adjoining streams. The annual removal effectiveness, which is the total NO removed in the SRB divided by the total NO draining from the field, ranged from 8 to 84%. This corresponds to an average removal rate of 0.040 g N m d with a range of 0.004 to 0.164 g N m d. Assuming a 40-yr life expectancy for the structure and a 4% discount rate, we computed a mean equal annual cost for SRBs of US$213.83. Given the average annual removal of 73 kg for all site-years, this cost equates to $2.94 kg N removed, which is very competitive with other field-edge practices such as denitrification bioreactors and constructed wetlands. Thus, SRBs continue to be a promising practice for NO removal in tile-drained landscapes.

Estimation of tile drainage contribution to streamflow and nutrient loads at the watershed scale based on continuously monitored data

Nitrogen losses from artificially drained watersheds degrade water quality at local and regional scales. In this study, we used an end-member mixing analysis (EMMA) together with high temporal resolution water quality and streamflow data collected in the 122 km² Otter Creek watershed located in northeast Iowa. We estimated the contribution of three end-members (groundwater, tile drainage, and quick flow) to streamflow and nitrogen loads and tested several combinations of possible nitrate concentrations for the end-members. Results indicated that subsurface tile drainage is responsible for at least 50% of the watershed nitrogen load between April 15 and November 1, 2015. Tiles delivered up to 80% of the stream N load while providing only 15-43% of the streamflow, whereas quick flows only marginally contributed to N loading. Data collected offer guidance about areas of the watershed that should be targeted for nitrogen export mitigation strategies.