Dynamics of critical source areas: does connectivity explain chemistry?

Phosphorus and iron-oxide transport from a hydrologically isolated grassland hillslope

Dissolved reactive phosphorus (DRP) loss from agricultural soils can negatively affect water quality. Shallow subsurface pathways can dominate P losses in grassland soils, especially in wetter months when waterlogging is common. This study investigated the processes controlling intra- and inter-event and seasonal DRP losses from poorly drained permanent grassland hillslope plots. Temporal flow related water samples were taken from surface runoff and subsurface (in-field pipe) discharge, analysed, and related to the likelihood of anaerobic conditions and redoximorphic species including nitrate (NO₃^-) over time. Subsurface drainage accounted for 89% of total losses. Simple linear regression and correlation matrices showed positive relationships between DRP and iron and soil moisture deficit; and negative relationships between these three factors and NO₃^- concentrations in drainage. These data indicate that waterlogging and low NO₃^- concentrations control the release of P in drainage, potentially via reductive dissolution. The relationship between DRP and metal release was less obvious in surface runoff, as nutrients gathered from P-rich topsoil camoflaged redox reactions. The data suggest a threshold in NO₃^- concentrations that could exacerbate P losses, even in low P soils. Knowledge of how nutrients interact with soil drainage throughout the year can be used to better time soil N and P inputs via, for example, fertiliser or grazing to avoid to excessive P loss that could harm water quality.

Upscaling of Denitrification Rates from Point to Catchment Scales for Modeling of Nitrate Transport and Retention

The spatial and temporal variability of denitrification makes it challenging to integrate conceptual, process-based understandings of nitrate transport and retention into numerical modeling at the catchment scale, although it is critical for the realism and predictive power of the model. In this study, we propose a novel approach where the conceptual understandings of the spatial structure of denitrification zones and the corresponding representative denitrification rates are transformed into a form that can be integrated into a multi-point statistical simulation framework. This is done by constructing a denitrification training image (TI) coupled to a geophysically based TI of the hydrogeological structure. The field observations and laboratory analyses of denitrification rates and the chemistry of water and sediment revealed that the study catchment's subsurface can be characterized by three zones: (1) the oxic zone with no nitrate reduction; (2) the slow-denitrification zone (mean of ln-transformed rate = -1.19 ± 0.52 mg N L^-1 yr^-1); and (3) the high-denitrification zone (mean of ln-transformed rate = 3.86 ± 1.96 mg N L^-1 yr^-1). The underlying controls on the spatial distribution of these zones and the representativeness of denitrification rates were investigated. Then, a TI illustrating the subsurface structure of the denitrification zone was constructed by synthesizing the results of these geochemical interpretations and the hydrogeology TI.

Livestock as vectors of organic matter and nutrient loading in aquatic ecosystems in African savannas

Populations of large wildlife have declined in many landscapes around the world, and have been replaced or displaced by livestock. The consequences of these changes on the transfer of organic matter (OM) and nutrients from terrestrial to aquatic ecosystems are not well understood. We used behavioural data, excretion and egestion rates and C: N: P stoichiometry of dung and urine of zebu cattle, to develop a metabolism-based estimate of loading rates of OM (dung), C, N and P into the Mara River, Kenya. We also directly measured the deposition of OM and urine by cattle into the river during watering. Per head, zebu cattle excrete and/or egest 25.6 g dry matter (DM, 99.6 g wet mass; metabolism) - 27.7 g DM (direct input) of OM, 16.0-21.8 g C, 5.9-9.6 g N, and 0.3-0.5 g P per day into the river. To replace loading rates OM of an individual hippopotamus by cattle, around 100 individuals will be needed, but much less for different elements. In parts of the investigated sub-catchments loading rates by cattle were equivalent to or higher than that of the hippopotamus. The patterns of increased suspended materials and nutrients as a result of livestock activity fit into historical findings on nutrients concentrations, dissolved organic carbon and other variables in agricultural and livestock areas in the Mara River basin. Changing these patterns of carbon and nutrient transport and cycling are having significant effects on the structure and functioning of both terrestrial and aquatic ecosystems.

Global mapping of freshwater nutrient enrichment and periphyton growth potential

Periphyton (viz. algal) growth in many freshwater systems is associated with severe eutrophication that can impair productive and recreational use of water by billions of people. However, there has been limited analysis of periphyton growth at a global level. To predict where nutrient over-enrichment and undesirable periphyton growth occurs, we combined several databases to model and map global dissolved and total nitrogen (N) and phosphorus (P) concentrations, climatic and catchment characteristics for up to 1406 larger rivers that were analysed between 1990 and 2016. We predict that 31% of the global landmass contained catchments may exhibit undesirable levels of periphyton growth. Almost three-quarters (76%) of undesirable periphyton growth was caused by P-enrichment and mapped to catchments dominated by agricultural land in North and South America and Europe containing 1.7B people. In contrast, undesirable periphyton growth due to N-enrichment was mapped to parts of North Africa and parts of the Middle East and India affecting 280 M people. The findings of this global modelling approach can be used by landowners and policy makers to better target investment and actions at finer spatial scales to remediate poor water quality owing to periphyton growth.

Synthesis of nutrient and sediment export patterns in the Chesapeake Bay watershed: Complex and non-stationary concentration-discharge relationships

Derived from river monitoring data, concentration-discharge (C-Q) relationships are useful indicators of riverine export dynamics. A top-down synthesis of C-Q patterns was conducted for suspended sediment (SS), total phosphorus (TP), and total nitrogen (TN) for nine major tributaries (15 monitoring sites) to Chesapeake Bay, which represent diverse characteristics in terms of land use, physiography, and hydrological settings. Model coefficients from the recently-developed Weighted Regressions on Time, Discharge, and Season (WRTDS) method were used to make informative interpretation of C-Q relationships. Unlike many previous C-Q studies that focused on stormflow conditions, this approach allows simultaneous examination of various discharge conditions within an uncertainty framework. This synthesis on WRTDS coefficients (i.e., the sensitivity of concentration to discharge) has offered new insights on the complexity of watershed function. Results show that watershed export has been dominated by mobilization patterns for SS and TP (particulate-dominated species) and chemostasis patterns for TN (dissolved-dominated species) under many river discharge conditions. Among nine possible modalities of low-flow vs. high-flow patterns, the three most frequent modalities are mobilization vs. mobilization (17 cases), chemostasis vs. mobilization (13 cases), and chemostasis vs. chemostasis (7 cases), representing 82% of all 45 watershed-constituent pairs. The general lack of dilution patterns may suggest that none of these constituents has been supply-limited in these watersheds. For many watershed-constituent combinations, results show clear temporal non-stationarity in C-Q relationships under selected time-invariant discharges, reflecting major changes in dominant watershed sources due to anthropogenic actions. These results highlight the potential pitfalls of assuming fixed C-Q relationships in the record. Overall, this work demonstrates the utility of WRTDS model coefficients for interpretation of river water-quality data and for generation of sensible hypotheses on dominant processes in different watersheds. The approach is readily adaptable to other river systems, where long-term discretely-sampled data are available, to decipher complex interactions between hydrological and biogeochemical processes.

High frequency monitoring of pesticides in runoff water to improve understanding of their transport and environmental impacts

Rainfall-induced peaks in pesticide concentrations can occur rapidly. Low frequency sampling may therefore largely underestimate maximum pesticide concentrations and fluxes. Detailed storm-based sampling of pesticide concentrations in runoff water to better predict pesticide sources, transport pathways and toxicity within the headwater catchments is lacking. High frequency monitoring (2min) of seven pesticides (Dimetomorph, Fluopicolide, Glyphosate, Iprovalicarb, Tebuconazole, Tetraconazole and Triadimenol) and one degradation product (AMPA) were assessed for 20 runoff events from 2009 to 2012 at the outlet of a vineyard catchment in the Layon catchment in France. The maximum pesticide concentrations were 387μgL^-1. Samples from all of the runoff events exceeded the legal limit of 0.1μgL^-1 for at least one pesticide (European directive 2013/39/EC). High resolution sampling used to detect the peak pesticide levels revealed that Toxic Units (TU) for algae, invertebrates and fish often exceeded the European Uniform principles (25%). The point and average (time or discharge-weighted) concentrations indicated up to a 30- or 4-fold underestimation of the TU obtained when measuring the maximum concentrations, respectively. This highlights the important role of sampling methods for assessing peak exposure. High resolution sampling combined with concentration-discharge hysteresis analyses revealed that clockwise responses were predominant (52%), indicating that Hortonian runoff is the prevailing surface runoff trigger mechanism in the study catchment. The hysteresis patterns for suspended solids and pesticides were highly dynamic and storm- and chemical-dependent. Intense rainfall events induced stronger C-Q hysteresis (magnitude). This study provides new insights into the complexity of pesticide dynamics in runoff water and highlights the ability of hysteresis analysis to improve understanding of pesticide supply and transport.

Impact of legacy soil phosphorus on losses in drainage and overland flow from grazed grassland soils

Rates and quantities of legacy soil phosphorus (P) lost from agricultural soils, and the timescales for positive change to water quality, remain unclear. From 2000 to 2004 five 0.2ha grazed grassland plots located on a drumlin hillslope in Northern Ireland, received chemical fertiliser applications of 0, 10, 20, 40, 80kgPha^-1yr^-1 resulting in soil Olsen P concentrations of 19, 24, 28, 38 and 67mgPL^-1, respectively, after which applications ceased. Soil Olsen P and losses to overland flow and drainage were monitored from 2005 to 2011 on an event and weekly flow proportional basis, respectively. Soluble reactive P and total P time series were synchronised with daily rainfall and modelled soil moisture deficits. From 2005 to 2011 soil Olsen P decline was proportional to soil P status with a 43% reduction in the plot at 67mgPL^-1 in 2004 and a corresponding 12% reduction in the plot with lowest soil P. However, there was no significant difference in the flow-weighted mean concentration for overland flow among plots, all of which exceeded 0.035mgL^-1 in >98% of events. Strong interannual and event variations in losses were observed with up to 65% of P being lost during a single rainfall event. P concentrations in drainage flow were independent of Olsen P and drain efficiency was potentially the primary control on concentrations, with the highest concentrations recorded in the plot at 38mgL^-1 Olsen P in 2004 (up to 2.72mgL^-1). Hydrological drivers, particularly antecedent soil moisture, had a strong influence on P loss in both overland and drainage flow, with higher concentrations recorded above a soil moisture deficit threshold of 7mm. This study demonstrates that on some soil types, legacy P poses a significant long term threat to water quality, even at agronomically optimum soil P levels.

Improving the identification of hydrologically sensitive areas using LiDAR DEMs for the delineation and mitigation of critical source areas of diffuse pollution

Identifying critical source areas (CSAs) of diffuse pollution in agricultural catchments requires the accurate identification of hydrologically sensitive areas (HSAs) at highest propensity for generating surface runoff and transporting pollutants. A new GIS-based HSA Index is presented that improves the identification of HSAs at the sub-field scale by accounting for microtopographic controls. The Index is based on high resolution LiDAR data and a soil topographic index (STI) and also considers the hydrological disconnection of overland flow via topographic impediment from flow sinks. The HSA Index was applied to four intensive agricultural catchments (~7.5-12km(2)) with contrasting topography and soil types, and validated using rainfall-quickflow measurements during saturated winter storm events in 2009-2014. Total flow sink volume capacities ranged from 8298 to 59,584m(3) and caused 8.5-24.2% of overland-flow-generating-areas and 16.8-33.4% of catchment areas to become hydrologically disconnected from the open drainage channel network. HSA maps identified 'breakthrough points' and 'delivery points' along surface runoff pathways as vulnerable points where diffuse pollutants could be transported between fields or delivered to the open drainage network, respectively. Using these as proposed locations for targeting mitigation measures such as riparian buffer strips reduced potential costs compared to blanket implementation within an example agri-environment scheme by 66% and 91% over 1 and 5years respectively, which included LiDAR DEM acquisition costs. The HSA Index can be used as a hydrologically realistic transport component within a fully evolved sub-field scale CSA model, and can also be used to guide the implementation of 'treatment-train' mitigation strategies concurrent with sustainable agricultural intensification.

Physico-chemical characteristics affect the spatial distribution of pesticide and transformation product loss to an agricultural brook

Diffuse entry of pesticide residues from agriculture into rivers is spatially unevenly distributed. Therefore, the identification of critical source areas (CSAs) may support water quality management in agricultural catchments. In contrast to former studies, we followed the hypothesis that not only hydrological and topographical characteristics but also physico-chemical properties of pesticide residues have a major influence on their loss to rivers and on corresponding formation of CSAs. We designed a virtual experiment, i.e. a numerical experiment as close as possible to environmental conditions, in a headwater catchment where pronounced spatial differences in hydrological transport processes were identified in the past. 144 scenarios with different combinations of adsorption coefficients (KOC = 10-1000 ml/g) and transformation half-lives (DT50 = 3-60 days) for pesticide parent compounds (PCs) and their transformation products (TPs) were simulated using the catchment-scale spatially distributed reactive transport model ZIN-AgriTra. Export fractions of substances in the virtual experiment ranged from 0.001-15% for pesticides and 0.001-1.8% for TPs. The results of the scenario investigations suggest that more of the calculated export mass variability could be attributed to KOC than to DT50 for both PCs and TPs. CSAs for TPs were spatially more equally distributed in the catchment than for PC export which was likely an effect of changing physico-chemical properties during transformation. The ranking of highest export fields was different between PCs and TPs for most of the investigated scenarios but six fields appeared among the top ten export fields in 95% of the scenarios, which shows the influence of site characteristics such as tile drains or soil properties in the catchment. Thus, the highest export fields were determined by a combination of site characteristics and substance characteristics. Therefore, despite the challenge of widely differing physico-chemical characteristics of pesticides on the market, these characteristics are an important consideration when delineating pesticide residue CSAs.

Implications of hydrologic connectivity between hillslopes and riparian zones on streamflow composition

Hydrological responses in mountainous headwater catchments are often highly non-linear with a distinct threshold-related behavior, which is associated to steep hillslopes, shallow soils and strong climatic variability. A holistic understanding of the dominant physical processes that control streamflow generation and non-linearity is required in order to assess potential negative effects of agricultural land use and water management in those areas. Therefore, streamflow generation in a small pre-Alpine headwater catchment (Upper Rietholzbach (URHB), ~1km(2)) was analyzed over a 2-year period by means of rainfall-response analysis and water quality data under explicit consideration of the joint behaviors of climate forcing and shallow groundwater dynamics. The runoff coefficients indicate that only a small fraction of the total catchment area (1-26%) generates streamflow during rainfall events. Hereby, the valley bottom areas (riparian zones) were the most important event-water source whereas only the lower parts of the hillslopes became hydrologically connected to the river network with higher antecedent moisture conditions. However, a distinct threshold-like behavior could not be observed, suggesting a more continuous shift from a riparian-zone to a more hillslope-dominated streamflow hydrograph. Regular manure application on the hillslopes in combinations with lateral hillslope groundwater flux and long groundwater residence times in the riparian zones resulted in a higher mineralization (e.g., total phosphorous) and significant denitrification in the valley bottom area. Despite the important role of the riparian zones for event-flow generation in the URHB, their nutrient buffer capacity is expected to be small due to the low permeability of the local subsurface material. The findings of this integrated analysis are summarized in a conceptual framework describing the hydrological functioning of hillslopes and riparian zones in the URHB.

Evaluating the critical source area concept of phosphorus loss from soils to water-bodies in agricultural catchments

Using data collected from six basins located across two hydrologically contrasting agricultural catchments, this study investigated whether transport metrics alone provide better estimates of storm phosphorus (P) loss from basins than critical source area (CSA) metrics which combine source factors as well. Concentrations and loads of P in quickflow (QF) were measured at basin outlets during four storm events and were compared with dynamic (QF magnitude) and static (extent of highly-connected, poorly-drained soils) transport metrics and a CSA metric (extent of highly-connected, poorly-drained soils with excess plant-available P). Pairwise comparisons between basins with similar CSA risks but contrasting QF magnitudes showed that QF flow-weighted mean TRP (total molybdate-reactive P) concentrations and loads were frequently (at least 11 of 14 comparisons) more than 40% higher in basins with the highest QF magnitudes. Furthermore, static transport metrics reliably discerned relative QF magnitudes between these basins. However, particulate P (PP) concentrations were often (6 of 14 comparisons) higher in basins with the lowest QF magnitudes, most likely due to soil-management activities (e.g. ploughing), in these predominantly arable basins at these times. Pairwise comparisons between basins with contrasting CSA risks and similar QF magnitudes showed that TRP and PP concentrations and loads did not reflect trends in CSA risk or QF magnitude. Static transport metrics did not discern relative QF magnitudes between these basins. In basins with contrasting transport risks, storm TRP concentrations and loads were well differentiated by dynamic or static transport metrics alone, regardless of differences in soil P. In basins with similar transport risks, dynamic transport metrics and P source information additional to soil P may be required to predict relative storm TRP concentrations and loads. Regardless of differences in transport risk, information on land use and management, may be required to predict relative differences in storm PP concentrations between these agricultural basins.

Fungicides transport in runoff from vineyard plot and catchment: contribution of non-target areas

Surface runoff and erosion during the course of rainfall events are major processes of pesticides transport from agricultural land to aquatic ecosystem. These processes are generally evaluated either at the plot or the catchment scale. Here, we compared at both scales the transport and partitioning in runoff water of two widely used fungicides, i.e., kresoxim-methyl (KM) and cyazofamid (CY). The objective was to evaluate the relationship between fungicides runoff from the plot and from the vineyard catchment. The results show that seasonal exports for KM and CY at the catchment were larger than those obtained at the plot. This underlines that non-target areas within the catchment largely contribute to the overall load of runoff-associated fungicides. Estimations show that 85 and 62 % of the loads observed for KM and CY at the catchment outlet cannot be explained by the vineyard plots. However, the partitioning of KM and CY between three fractions, i.e., the suspended solids (>0.7 μm) and two dissolved fractions (i.e., between 0.22 and 0.7 µm and <0.22 µm) in runoff water was similar at both scales. KM was predominantly detected below 0.22 μm, whereas CY was mainly detected in the fraction between 0.22 and 0.7 μm. Although KM and CY have similar physicochemical properties and are expected to behave similarly, our results show that their partitioning between two fractions of the dissolved phase differs largely. It is concluded that combined observations of pesticide runoff at both the catchment and the plot scales enable to evaluate the sources areas of pesticide off-site transport.