Geology Reports⌕ Search

USGS · 70185321

Combined use of groundwater dating, chemical, and isotopic analyses to resolve the history and fate of nitrate contamination in two agricultural watersheds, Atlantic coastal plain, Maryland

Abstract

The history and fate of groundwater nitrate (NO 3 − ) contamination were compared in 2 small adjacent agricultural watersheds in the Atlantic coastal plain by combined use of chronologic (CCl 2 F 2 , 3 H), chemical (dissolved solids, gases), and isotopic (δ 15 N,δ 13 C, δ 34 S) analyses of recharging groundwaters, discharging groundwaters, and surface waters. The results demonstrate the interactive effects of changing agricultural practices, groundwater residence times, and local geologic features on the transfer of NO 3 − through local flow systems. Recharge dates of groundwaters taken in 1990–1992 from the surficial aquifer in the Chesterville Branch and Morgan Creek watersheds near Locust Grove, Maryland, ranged from pre-1940 to the late 1980’s. When corrected for localized denitrification by use of dissolved gas concentrations, the dated waters provide a 40-year record of the recharge rate of NO 3 − , which increased in both watersheds by a factor of 3–6, most rapidly in the 1970's. The increase in groundwater NO 3 − over time was approximately proportional to the documented increase in regional N fertilizer use, and could be accounted for by oxidation and leaching of about 20–35% of the fertilizer N. Groundwaters discharging upward beneath streams in both watersheds had measured recharge dates from pre-1940 to 1975, while chemical data for second-order reaches of the streams indicated average groundwater residence times in the order of 20+ years. At the time of the study, NO 3 − discharge rates were less than NO 3 − recharge rates for at least two reasons: (1) discharge of relatively old waters with low initial NO 3 − concentrations, and (2) local denitrification. In the Chesterville Branch watershed, groundwaters remained oxic throughout much of the surficial aquifer and discharged relatively unaltered to the stream, which had a relatively high NO 3 − concentration (9–10 mg/L as N). In the Morgan Creek watershed, groundwaters were largely reduced and denitrified before discharging to the stream, which had a relatively low NO 3 − concentration (2–3 mg/L as N). Chemical and isotopic data indicate that quantitative denitrification occurred within buried calcareous glauconitic marine sediments that are present at relatively shallow depths beneath the Morgan Creek watershed. NO 3 − removal by forests, wetlands, and shallow organic-rich soils in near-stream environments was largely avoided by groundwaters that followed relatively deep flow paths before converging and discharging rapidly upward to the streams.

Explore related subjects

90° N90° S · 180° W ← longitude → 180° E
Source-reported bounding extent: 39.11727568585598° to 39.46588451142044° latitude; -76.2506103515625° to -75.76171875° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

J.K. Böhlke, J. M. Denver. 2010-07-09. Combined use of groundwater dating, chemical, and isotopic analyses to resolve the history and fate of nitrate contamination in two agricultural watersheds, Atlantic coastal plain, Maryland. https://doi.org/10.1029/95wr01584

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

A demonstrated workflow to efficiently identify better performing process-based model configurations

There are often competing options for process representation when designing process-based hydrologic models. Comparing candidate model configurations during model development can be computationally expensive and existing methods for evaluating prior-data conflict often cannot distinguish between multiple models with low prior-data conflict. Here, we demonstrate a workflow to efficiently identify better performing models by comparing prior model output distributions to observed data distributions using a metric referred to as Mahalanobis distance deviation (𝑀 D ). The 𝑀 D metric was developed to incorporate the notion of balancing model fit with model complexity into prior-data conflict approaches. Ninety-six different configurations of the Precipitation-Runoff Modeling System were designed for the East River watershed, Colorado by considering uncertainty in representing spatial discretization, forcing data, potential evapotranspiration, soil zone flow, snow depletion curve, and transpiration period processes. The prior-to-calibration adequacy of each model was estimated by calculating 𝑀 D for observed streamflow, areal runoff, actual evapotranspiration, groundwater recharge, soil moisture, and snow-covered area. All 96 models were then calibrated using the PEST++ iterative ensemble smoother to quantify their performance and validate the 𝑀 D results. The 𝑀 D metrics calculated from prior parameter ensembles as small as 25 exhibited strong correlation with calibrated objective function values highlighting the substantial information on model performance that can be gained prior to calibration through the workflow. Spatial discretization and forcing data were most influential on model performance and their better performing representations were most easily identified. This work demonstrates an efficient approach to identify better performing model configurations and the processes driving those improvements.

Water Resources Research↗

Evaluating SWOT satellite data with gage information in the United States and Canada

The Surface Water and Ocean Topography (SWOT) satellite mission represents a transformative approach to water resources management, offering high-resolution data on water-surface elevations (WSE) and extents. This study evaluates SWOT WSE accuracy in the United States and Canada with in situ streamgage and lake and reservoir gage data. We compare SWOT Level 2 Hydrology data products (Version D) to in situ timeseries of WSE on 335,531 paired observations across 2,629 gages. SWOT data accuracy and temporal coverage are explored as functions of various SWOT quality filters. After applying a standard set of defined filters and thresholds to compare paired measurements of SWOT and in situ observations (both referenced to the Earth Gravitational Model 2008), we find 1-σ (68th percentile) absolute differences of 20 and 11 cm, and median biases of 2.8 and 0.1 cm, for filtered river node and lake data, respectively. River reach and lake comparisons may be biased due to scale mismatches between in situ point-based data and spatially averaged SWOT data. Applying SWOT data quality filters improved median paired elevation R2 between unfiltered and filtered data from 0.12 to 0.92 across all reach, node, and lake observations. Overall, we find that 60.0% and 73.8% of all paired SWOT node and lake WSE observations that satisfy SWOT-observable criteria (ice-free observations from sufficiently large waterbodies within the SWOT cross track range) have elevations within 25 cm of reported gage observations. These comparisons also provide a baseline benchmark that can help inform continued SWOT algorithm refinement.

Water Resources Research↗

Big data for small streams: Leveraging a coarse national model to enhance fine-resolution stream temperature predictions

Predictive stream temperature models are commonly used to address spatial and temporal gaps in observational data because stream temperature information is essential for environmental management decisions. National models capture broad patterns driven by climate, topography, and seasonality, but often lack the spatial resolution to simulate fine-scale variation in groundwater discharge, riparian shading, or tributary inflows. In contrast, watershed models may simulate finer-resolution patterns but miss regional drivers or be limited by scarcity of observations. To generate watershed-extent, fine-resolution stream temperature predictions, we developed 18 multiscale machine-learning models that leverage (a) a broad-extent, coarse-resolution stream temperature model and (b) finer-resolution observations and predictor data. The coarse-resolution model simulates daily stream temperature across the conterminous United States on the National Hydrologic Geospatial Fabric (13.7 km median reach length; 865,364 km total). The fine-resolution models simulate daily stream temperature on the National Hydrography Dataset Plus version 2 (1.5 km median; 5,192,808 km total). We compared multiscale models to seven purely broad- (national) or watershed-extent models, testing them in eight watersheds (436–985 km 2 ) with differing hydrologic settings and under realistic data availability scenarios. Especially when observations within a watershed were scarce, multiscale models generated more accurate predictions than national or watershed-specific models. A recurrent graph convolution neural network pretrained on coarse predictions and fine-tuned on watershed-specific observations was most accurate across watersheds and data scenarios. These results highlight that leveraging both coarse, national results and watershed-specific data can increase model accuracies, particularly when observations within a watershed are spatially or temporally limited.

Water Resources Research↗