Understanding long-term baseflow water quality trends using a synoptic survey of the ground water-surface water interface, central Wisconsin

Forecasting and Hindcasting PFAS Concentrations in Groundwater Discharging to Streams near a PFAS Production Facility

Per- and polyfluoroalkyl substances (PFAS) are known to be highly persistent in groundwater, making it vital to develop new approaches to important practical questions such as the time scale for future persistence of PFAS in contaminated groundwater. In the approach presented here, groundwater from beneath streambeds was analyzed for PFAS and age-dated using SF6 and 3H/3He. The results were coupled with groundwater flux measurements in a convolution approach to estimate past and future PFAS concentrations in groundwater discharge to the streams. At our test site near the Cape Fear River (CFR) of North Carolina, PFAS were detected in groundwater up to 43 years old, suggesting that some PFAS entered groundwater immediately or shortly after fluorochemical production began at the nearby Fayetteville Works. Results are consistent with little to no retardation in groundwater for perfluoroethers such as hexafluoropropylene oxide-dimer acid (HFPO-DA) and perfluoro-2-methoxypropanoic acid (PMPA), the two most abundant PFAS, with mean concentrations of 229 and 498 ng/L, respectively. Future PFAS concentrations in groundwater discharge to streams were estimated to remain above current MCL or health advisory levels through at least 2050 or 2060 (using 3H/3He and SF6, respectively). Recent atmospheric deposition data suggest lower but non-negligible amounts of PFAS may continue to enter groundwater, likely further extending PFAS persistence in groundwater and the adjacent CFR. This approach shows promise for giving an overall perspective on persistence of PFAS in groundwater discharge from a broad contaminated area.

Predicting Redox Conditions in Groundwater at a National Scale Using Random Forest Classification

Redox conditions in groundwater may markedly affect the fate and transport of nutrients, volatile organic compounds, and trace metals, with significant implications for human health. While many local assessments of redox conditions have been made, the spatial variability of redox reaction rates makes the determination of redox conditions at regional or national scales problematic. In this study, redox conditions in groundwater were predicted for the contiguous United States using random forest classification by relating measured water quality data from over 30,000 wells to natural and anthropogenic factors. The model correctly predicted the oxic/suboxic classification for 78 and 79% of the samples in the out-of-bag and hold-out data sets, respectively. Variables describing geology, hydrology, soil properties, and hydrologic position were among the most important factors affecting the likelihood of oxic conditions in groundwater. Important model variables tended to relate to aquifer recharge, groundwater travel time, or prevalence of electron donors, which are key drivers of redox conditions in groundwater. Partial dependence plots suggested that the likelihood of oxic conditions in groundwater decreased sharply as streams were approached and gradually as the depth below the water table increased. The probability of oxic groundwater increased as base flow index values increased, likely due to the prevalence of well-drained soils and geologic materials in high base flow index areas. The likelihood of oxic conditions increased as topographic wetness index (TWI) values decreased. High topographic wetness index values occur in areas with a propensity for standing water and overland flow, conditions that limit the delivery of dissolved oxygen to groundwater by recharge; higher TWI values also tend to occur in discharge areas, which may contain groundwater with long travel times. A second model was developed to predict the probability of elevated manganese (Mn) concentrations in groundwater (i.e., ≥50 μg/L). The Mn model relied on many of the same variables as the oxic/suboxic model and may be used to identify areas where Mn-reducing conditions occur and where there is an increased risk to domestic water supplies due to high Mn concentrations. Model predictions of redox conditions in groundwater produced in this study may help identify regions of the country with elevated groundwater vulnerability and stream vulnerability to groundwater-derived contaminants.

Assessment of River-Aquifer Interaction and Nitrogen Contamination in an Agricultural Zone Located in the Guarani Aquifer System Outcrop Area

The excessive use of nitrogen fertilizers is responsible for an increase in nitrate concentrations in water bodies, which in the future could led to an irreversible contamination compromising the water resource quality. In this way, understand the water movement within a watershed and evaluate the impacts related to agricultural practices is relevant for water management, especially in an environmental fragile region, such as the outcrop area of the Guarani Aquifer. Water samples from a small watershed located at the Guarani Aquifer region in São Paulo state, representing surface water and groundwater discharge in riverbeds from two creeks, as well as groundwater (springs and wells) were collected for isotopic ratios (δ18O e δ2H) and nitrate determination. The results indicated that the river flow is mostly supplied by groundwater discharge, and despite the observed concentrations of nitrate in groundwater reaching the creeks, the current scenario indicate contamination in the surface water, above the regulatory levels. Therefore, the expansion in sugarcane production increases the possibility that the released nitrate reaches high levels in the future in this watershed.

Influence of redox gradients on nitrate transport from the landscape to groundwater and streams

Increases in nitrogen applications to the land surface since the 1950s have led to a cascade of negative environmental impacts, including degradation of drinking water supplies, nutrient enrichment of aquatic ecosystems and contributions to global climate change. In this study, groundwater, streambed porewater, and stream sampling were used to establish trends in nitrate concentrations and how redox gradients influence nitrate transport across diverse glacial terranes. Decadal sampling has found that elevated nitrate concentrations in shallow groundwater beneath cropland have been sustained for decades. Redox gradients established in the saturated zone using dissolved O₂, iron, nitrate and excess N₂ from denitrification suggest that nitrate-bearing zones are thin in glacial terranes dominated by fine materials. These thin nitrate-bearing zones lead to suboxic, low nitrate streambed porewater and limit the contributions of nitrate to streams from slow-flow groundwater. In contrast, thick oxic zones in more coarse-grained glacial terranes allow nitrate to reach deeper groundwater, resulting in streambed porewater with elevated nitrate concentrations and causing a large portion of stream nitrate to be derived from slow-flow groundwater. Groundwater age tracer data indicate that denitrification occurs more quickly in the terrane dominated by fine material than in the more coarse-grained terrane. The quicker depletion of nitrate in the more fine-grained terrane suggests that the thinner oxic zone in this terrane is due, in part, to the greater availability and reactivity of electron donors in this terrane than in the more coarse-grained terrane. Groundwater age tracer data and hydrograph separation analysis suggest that saturated zone lag times between when changes in land use practices occur and when changes in stream water are fully observed may vary widely across hydrogeologic settings.

Re-Purposing Groundwater Flow Models for Age Assessments: Important Characteristics

Groundwater flow model construction is often time-consuming and costly, with development ideally focused on a specific purpose, such as quantifying well capture from water bodies or providing flow fields for simulating advective transport. As environmental challenges evolve, the incentive to re-purpose existing groundwater flow models may increase. However, few studies have evaluated which characteristics of groundwater flow models deserve greatest consideration when re-purposing models for groundwater age and advective transport simulations. In this paper, we compare simulated age metrics produced by three MODFLOW-MODPATH models of the same area but with differing levels of complexity (layering and heterogeneity). Comparisons are made at three watershed scales (HUC 8 to HUC 12). Groundwater age metrics, specifically the young fraction and median age of the young and old fractions, are used for evaluation because they relate to intrinsic susceptibility of aquifers and are simpler to interpret than full age distributions used for advective transport. Results indicate that: (1) the young fraction is less sensitive to model layering than the median age of young and old fractions, suggesting that simple models may suffice for basic intrinsic susceptibility assessments; (2) water table mounding and associated discharge into partially penetrating boundaries, such as head-water streams, is important for simulating both the young fraction and the median age of the young fraction; and (3) the influence of partially penetrating head-water streams is maintained regardless of the porosity distribution. Results of this work should aid modelers with evaluating the appropriateness of re-purposing existing groundwater flow models for age simulations.

Streambed Flux Measurement Informed by Distributed Temperature Sensing Leads to a Significantly Different Characterization of Groundwater Discharge

Groundwater discharge though streambeds is often focused toward discrete zones, indicating that preliminary reconnaissance may be useful for capturing the full spectrum of groundwater discharge rates using point-scale quantitative methods. However, many direct-contact reconnaissance techniques can be time-consuming, and remote sensing (e.g., thermal infrared) typically does not penetrate the water column to locate submerged seepages. In this study, we tested whether dozens of groundwater discharge measurements made at “uninformed” (i.e., selected without knowledge on high-resolution temperature variations at the streambed) point locations along a reach would yield significantly different Darcy-based groundwater discharge rates when compared with “informed” measurements, focused at streambed thermal anomalies that were identified a priori using fiber-optic distributed temperature sensing (FO-DTS). A non-parametric U-test showed a significant difference between median discharge rates for uninformed (0.05 m·day−1; n = 30) and informed (0.17 m·day−1; n = 20) measurement locations. Mean values followed a similar pattern (0.12 versus 0.27 m·day−1), and frequency distributions for uninformed and informed measurements were also significantly different based on a Kolmogorov–Smirnov test. Results suggest that even using a quick “snapshot-in-time” field analysis of FO-DTS data can be useful in streambeds with groundwater discharge rates <0.2 m·day−1, a lower threshold than proposed in a previous study. Collectively, study results highlight that FO-DTS is a powerful technique for identifying higher-discharge zones in streambeds, but the pros and cons of informed and uninformed sampling depend in part on groundwater/surface water exchange study goals. For example, studies focused on measuring representative groundwater and solute fluxes may be biased if high-discharge locations are preferentially sampled. However, identification of high-discharge locations may complement more randomized sampling plans and lead to improvements in interpolating streambed fluxes and upscaling point measurements to the stream reach scale.

Vulnerability of streams to legacy nitrate sources

The influence of hydrogeologic setting on the susceptibility of streams to legacy nitrate was examined at seven study sites having a wide range of base flow index (BFI) values. BFI is the ratio of base flow to total streamflow volume. The portion of annual stream nitrate loads from base flow was strongly correlated with BFI. Furthermore, dissolved oxygen concentrations in streambed pore water were significantly higher in high BFI watersheds than in low BFI watersheds suggesting that geochemical conditions favor nitrate transport through the bed when BFI is high. Results from a groundwater-surface water interaction study at a high BFI watershed indicate that decades old nitrate-laden water is discharging to this stream. These findings indicate that high nitrate levels in this stream may be sustained for decades to come regardless of current practices. It is hypothesized that a first approximation of stream vulnerability to legacy nutrients may be made by geospatial analysis of watersheds with high nitrogen inputs and a strong connection to groundwater (e.g., high BFI).

Agriculture-related trends in groundwater quality of the glacial deposits aquifer, central Wisconsin

Measuring and understanding trends in groundwater quality is necessary for determining whether changes in land-management practices have an effect on groundwater quality. This paper describes an approach that was used to measure and understand trends using data from two groundwater studies conducted in central Wisconsin as part of the USGS NAWQA program. One of the key components of this approach, determining the age of sampled groundwater, gave a temporal component to the snapshots of water quality that were obtained through synoptic-sampling efforts. This approach can be used at other locations where groundwater quality data are collected, groundwater age can be determined, and associated temporal data are available. Results of these studies indicate measured concentrations of nitrate and atrazine plus deethylatrazine were correlated to historical patterns of fertilizer and atrazine use. Concentrations of nitrate in groundwater have increased over time; concentrations of atrazine plus deethylatrazine increased and then decreased. Concentrations of nitrate also were correlated to screen depth below the water level and concentrations of dissolved O(2); concentrations of atrazine plus deethylatrazine were correlated to dissolved O(2) and annual precipitation. To measure trends in concentrations of atrazine plus deethylatrazine, the data, collected over a near-decadal period, were adjusted to account for changes in laboratory-reporting levels and analytical recoveries. Only after accounting for these changes was it apparent that the median concentrations of atrazine plus deethylatrazine decreased over the near-decadal interval between sampling efforts.

Collateral geochemical impacts of agricultural nitrogen enrichment from 1963 to 1985: a southern Wisconsin ground water depth profile

In this study, we used chlorofluorocarbon (CFC) age-dating to investigate the geochemistry of N enrichment within a bedrock aquifer depth profile beneath a south central Wisconsin agricultural landscape. Measurement of N(2)O and excess N(2) allowed us to reconstruct the total NO(3)(-) and total nitrogen (TN) leached to ground water and was essential for tracing the separate influences of soil nitrification and ground water denitrification in the collateral geochemical chronology. We identify four geochemical impacts due to a steady ground water N enrichment trajectory (39 +/- 2.2 micromol L(-1) yr(-1), r(2) = 0.96) over two decades (1963-1985) of rapidly escalating N use. First, as a by-product of soil nitrification, N(2)O entered ground water at a stable (r(2) = 0.99) mole ratio of 0.24 +/- 0.007 mole% (N(2)O-N/NO(3)-N). The gathering of excess N(2)O in ground water is a potential concern relative to greenhouse gas emissions and stratospheric ozone depletion after it discharges to surface water. Second, excess N(2) measurements revealed that NO(3)(-) was a prominent, mobile, labile electron acceptor comparable in importance to O(2.) Denitrification transformed 36 +/- 15 mole% (mol mol(-1) x 100) of the total N within the profile to N(2) gas, delaying exceedance of the NO(3)(-) drinking water standard by approximately 6 yr. Third, soil acids produced from nitrification substantially increased the concentrations of major, dolomitic ions (Ca, Mg, HCO(3)(-)) in ground water relative to pre-enrichment conditions. By 1985, concentrations approximately doubled; by 2006, CFC age-date projections suggest concentrations may have tripled. Finally, the nitrification induced mobilization of Ca may have caused a co-release of P from Ca-rich soil surfaces. Dissolved P increased from an approximate background value of 0.02 mg L(-1) in 1963 to 0.07 mg L(-1) in 1985. The CFC age-date projections suggest the concentration could have reached 0.11 mg L(-1) in ground water recharge by 2006. These results highlight an intersection of the N and P cycles potentially important for managing the quality of ground water discharged to surface water.

Linking ground-water age and chemistry data along flow paths: implications for trends and transformations of nitrate and pesticides

Tracer-based ground-water ages, along with the concentrations of pesticides, nitrogen species, and other redox-active constituents, were used to evaluate the trends and transformations of agricultural chemicals along flow paths in diverse hydrogeologic settings. A range of conditions affecting the transformation of nitrate and pesticides (e.g., thickness of unsaturated zone, redox conditions) was examined at study sites in Georgia, North Carolina, Wisconsin, and California. Deethylatrazine (DEA), a transformation product of atrazine, was typically present at concentrations higher than those of atrazine at study sites with thick unsaturated zones but not at sites with thin unsaturated zones. Furthermore, the fraction of atrazine plus DEA that was present as DEA did not increase as a function of ground-water age. These findings suggest that atrazine degradation occurs primarily in the unsaturated zone with little or no degradation in the saturated zone. Similar observations were also made for metolachlor and alachlor. The fraction of the initial nitrate concentration found as excess N2 (N2 derived from denitrification) increased with ground-water age only at the North Carolina site, where oxic conditions were generally limited to the top 5 m of saturated thickness. Historical trends in fluxes to ground water were evaluated by relating the times of recharge of ground-water samples, estimated using chlorofluorocarbon concentrations, with concentrations of the parent compound at the time of recharge, estimated by summing the molar concentrations of the parent compound and its transformation products in the age-dated sample. Using this approach, nitrate concentrations were estimated to have increased markedly from 1960 to the present at all study sites. Trends in concentrations of atrazine, metolachlor, alachlor, and their degradates were related to the timing of introduction and use of these compounds. Degradates, and to a lesser extent parent compounds, were detected in ground water dating back to the time these compounds were introduced.