Determining the groundwater contribution of nonpoint source pollution at a watershed scale is a challenging issue. In this study, we utilized a top‐down approach to characterize representative groundwater response units (GRUs) based on land use and landscape position (e.g., upland, sideslope, or floodplain) in the 275‐km2Clear Creek Watershed, Iowa. Groundwater monitoring wells were then established along downslope transects in representative GRUs. This unique combination of top‐down/bottom‐up approaches allowed us to estimate groundwater pollutant loads at the watershed scale with minimal monitoring. For the 2015 study period, results indicated that more groundwater recharge occurred in the floodplain (404 mm) compared to the uplands or sideslopes (281 and 165 mm, respectively), irrespective of land use. Recharge in the floodplains consisted of 37% of the annual precipitation, whereas upland wells averaged 26% and sideslopes averaged 15% of the annual precipitation. Less recharge was found to occur beneath perennial grass compared to row crop and urbanized areas. Baseflow discharge accounted for 69% of the total NO3‐N exported from the Clear Creek Watershed, with row crop areas contributing approximately 95% of the annual load. Orthophosphorus (OP) yields were approximately 0.72 kg/ha beneath urban and suburban areas, three times higher than those in row crop or perennial areas. Urban and suburban areas accounted for 21.4% of groundwater orthophosphorus and chloride loads in the watershed compared to only 8.5% of the land area. Overall, the groundwater load allocation model for baseflow nutrient discharge to Clear Creek can be used to target future nonpoint source load reduction strategies at the watershed scale. The use of GRUs can pinpoint better areas of concern for controlling nutrient loads.
Groundwater transit time distributions (TTDs) describe the spectrum of flow‐weighted apparent ages of groundwater from aquifer recharge to discharge. Regional‐scale TTDs in stream baseflow are often estimated from numerical models with limited calibration from groundwater sampling and suggest much younger groundwater discharge than has been observed by discrete age‐dating techniques. We investigate both local and regional‐scale groundwater TTDs in the Upper Middle Loup watershed (5,440 km2) overlying the High Plains Aquifer in the Nebraska Sand Hills, USA. We determined flow‐weighted apparent ages of groundwater discharging through the streambed at 88 discrete points along a 99 km groundwater‐dominated stream segment using3H, noble gases,14C, and groundwater flux measurements at the point‐scale (<7.6 cm diameter). Points were organized in transects across the stream width (3–10 points per transect) and transects were clustered in five sampling areas (10–610 m in stream length) located at increasing distances along the stream. Groundwater apparent ages ranged from 0 to 8,200 years and the mean groundwater transit time along the 99 km stream is >3,000 years. TTDs from upstream sampling areas were best fit by distributions with a narrow range of apparent ages, but when older groundwater from downstream sampling areas is included, the regional TTD is scale dependent and the distribution is better described by a gamma model (
- NSF-PAR ID:
- 10494512
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Water Resources Research
- Volume:
- 60
- Issue:
- 2
- ISSN:
- 0043-1397
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
more » « less
Abstract We utilized 251 measurements from a recently developed automated seepage meter (ASM) in streambeds in the Nebraska Sand Hills, USA to investigate the small‐scale spatial variability of groundwater seepage flux (q) and the ability of the ASM to estimate mean q at larger scales. Small‐scale spatial variability of q was analyzed in five dense arrays, each covering an area of 13.5–28.0 m2(169 total point measurements). Streambed vertical hydraulic conductivity (K) was also measured. Results provided: (a) high‐resolution contour plots of q and K, (b) anisotropic semi‐variograms demonstrating greater correlation scales of q and K along the stream length than across the stream width, and (c) the number of rows of points (perpendicular to streamflow) needed to represent the groundwater flux of areas up to 28.0 m2. The findings suggest that representative streambed measurements are best conducted perpendicular to streamflow to accommodate larger seepage flux heterogeneity in this direction and minimize sampling redundancy. To investigate the ASM's ability to produce accurate mean q at larger scales, seepage meters were deployed in four stream reaches (170–890 m), arranged in three to six transects (three to eight points each) per reach across the channel. In each reach, the mean seepage flux from ASMs was compared to the seepage flux from bromide tracer dilution. Agreement between the two methods indicates the viability of a modest number of seepage meter measurements to determine the overall groundwater flux to the stream and can guide sampling for solutes and environmental tracers.
-
more » « less
Abstract Groundwater discharge to streams is a nonpoint source of nitrogen (N) that confounds N mitigation efforts and represents a significant portion of the annual N loading to watersheds. However, we lack an understanding of where and how much groundwater N enters streams and watersheds. Nitrogen concentrations at the end of groundwater flowpaths are the culmination of biogeochemical and physical processes from the contributing land area where groundwater recharges, within the aquifer system, and in the near-stream riparian area where groundwater discharges to streams. Our research objectives were to quantify the spatial distribution of N concentrations at groundwater discharges throughout a mixed land-use watershed and to evaluate how relationships among contributing and riparian land cover, modeled aquifer characteristics, and groundwater discharge biogeochemistry explain the spatial variation in groundwater discharge N concentrations. We accomplished this by integrating high-resolution thermal infrared surveys to locate groundwater discharge, biogeochemical sampling of groundwater, and a particle tracking model that links groundwater discharge locations to their contributing area land cover. Groundwater N loading from groundwater discharges within the watershed varied substantially between and within streambank groundwater discharge features. Groundwater nitrate concentrations were spatially heterogeneous ranging from below 0.03–11.45 mg-N/L, varying up to 20-fold within meters. When combined with the particle tracking model results and land cover metrics, we found that groundwater discharge nitrate concentrations were best predicted by a linear mixed-effect model that explained over 60% of the variation in nitrate concentrations, including aquifer chemistry (dissolved oxygen, Cl−, SO42−), riparian area forested land cover, and modeled physical aquifer characteristics (discharge, Euclidean distance). Our work highlights the significant spatial variability in groundwater discharge nitrate concentrations within mixed land-use watersheds and the need to understand groundwater N processing across the many spatiotemporal scales within groundwater cycling.
-
more » « less
Abstract Most terrestrial allochthonous organic matter enters river networks through headwater streams during high flow events. In headwaters, allochthonous inputs are substantial and variable, but become less important in streams and rivers with larger watersheds. As allochthonous dissolved organic matter (DOM) moves downstream, the proportion of less aromatic organic matter with autochthonous characteristics increases. How environmental factors converge to control this transformation of DOM at a continental scale is less certain. We hypothesized that the amount of time water has spent travelling through surface waters of inland systems (streams, rivers, lakes, and reservoirs) is correlated to DOM composition. To test this hypothesis, we used established river network scaling relationships to predict relative river network flow‐weighted travel time (FWTT) of water for 60 stream and river sites across the contiguous United States (3090 discrete samples over 10 water years). We estimated lentic contribution to travel times with upstream in‐network lake and reservoir volume. DOM composition was quantified using ultraviolet and visible absorption and fluorescence spectroscopy. A combination of FWTT and lake and reservoir volume was the best overall predictor of DOM composition among models that also incorporated discharge, specific discharge, watershed area, and upstream channel length. DOM spectral slope ratio (R2 = 0.77) and Freshness Index (R2 = 0.78) increased and specific ultraviolet absorbance at 254 nm (R2 = 0.68) and Humification Index (R2 = 0.44) decreased across sites as a function of FWTT and upstream lake volume. This indicates autochthonous‐like DOM becomes continually more dominant in waters with greater FWTT. We assert that river FWTT can be used as a metric of the continuum of DOM composition from headwaters to rivers. The nature of the changes to DOM composition detected suggest this continuum is driven by a combination of photo‐oxidation, biological processes, hydrologically varying terrestrial subsidies, and aged groundwater inputs.
-
Abstract Non-perennial rivers and streams make up over half the global river network and are becoming more widespread. Transitions from perennial to non-perennial flow are a threshold-type change that can lead to alternative stable states in aquatic ecosystems, but it is unknown whether streamflow itself is stable in either wet (flowing) or dry (no-flow) conditions. Here, we investigated drivers and feedbacks associated with regime shifts between wet and dry conditions in an intermittent reach of the Arkansas River (USA) over the past 23 years. Multiple lines of evidence suggested that these regimes represent alternative stable states, including (a) significant jumps in discharge time series that were not accompanied by jumps in flow drivers such as precipitation and groundwater pumping; (b) a multi-modal state distribution with 92% of months experiencing no-flow conditions for <10% or >90% of days, despite unimodal distributions of precipitation and pumping; and (c) a hysteretic relationship between climate and flow state. Groundwater levels appear to be the primary control over the hydrological regime, as groundwater levels in the alluvial aquifer were higher than the stream stage during wet regimes and lower than the streambed during dry regimes. Groundwater level variation, in turn, was driven by processes occurring at both the regional scale (surface water inflows from upstream, groundwater pumping) and the reach scale (stream–aquifer exchange, diffuse recharge through the soil column). Historical regime shifts were associated with diverse pressures including network disconnection caused by upstream water use, increased flow stability potentially associated with reservoir operations, and anomalous wet and dry climate conditions. In sum, stabilizing feedbacks among upstream inflows, stream–aquifer interactions, climate, vegetation, and pumping appear to create alternative wet and dry stable states at this site. These stabilizing feedbacks suggest that widespread observed shifts from perennial to non-perennial flow will be difficult to reverse.more » « less
