USGS ScienceSearch

USGS · wri034035

Residence times and nitrate transport in ground water discharging to streams in the Chesapeake Bay Watershed

Abstract

One of the major water-quality problems in the Chesapeake Bay is an overabundance of nutrients from the streams and rivers that discharge to the Bay. Some of these nutrients are from nonpoint sources such as atmospheric deposition, agricultural manure and fertilizer, and septic systems. The effects of efforts to control nonpoint sources, however, can be difficult to quantify because of the lag time between changes at the land surface and the response in the base-flow (ground water) component of streams. To help resource managers understand the lag time between implementation of management practices and subsequent response in the nutrient concentrations in the base-flow component of streamflow, a study of ground-water discharge, residence time, and nitrate transport in springs throughout the Chesapeake Bay Watershed and in four smaller watersheds in selected hydrogeomorphic regions (HGMRs) was conducted. The four watersheds were in the Coastal Plain Uplands, Piedmont crystalline, Valley and Ridge carbonate, and Valley and Ridge siliciclastic HGMRs. A study of springs to estimate an apparent age of the ground water was based on analyses for concentrations of chlorofluorocarbons in water samples collected from 48 springs in the Chesapeake Bay Watershed. Results of the analysis indicate that median age for all the samples was 10 years, with the 25th percentile having an age of 7 years and the 75th percentile having an age of 13 years. Although the number of samples collected in each HGMR was limited, there did not appear to be distinct differences in the ages between the HGMRs. The ranges were similar between the major HGMRs above the Fall Line (modern to about 50 years), with only two HGMRs of small geographic extent (Piedmont carbonate and Mesozoic Lowland) having ranges of modern to about 10 years. The median values of all the HGMRs ranged from 7 to 11 years. Not enough samples were collected in the Coastal Plain for comparison. Spring samples showed slightly younger water under wet conditions than under dry conditions. The apparent age of water from wells, springs, and other ground-water discharge points in the four targeted watersheds was modern to 60 years, which was similar to the apparent ages from the spring study. In the Pocomoke River Watershed in the Coastal Plain Uplands HGMR, the apparent age of ground-water samples ranged from 0 to 60 years; the ages in the vicinity of the streams ranged from 0 to 23 years. The apparent ages of ground water in the Polecat Creek Watershed in the Piedmont crystalline HGMR ranged from 2 to 30 years. The apparent ages of water from wells in the Muddy Creek Watershed in the Valley and Ridge carbonate HGMR ranged from 10 to 20 years (except for a single sample that was 45 years). The ages in the East Mahantango Creek Watershed in the Valley and Ridge siliciclastic HGMR ranged from 0 to 50 years. The distribution in apparent age of water from wells in the targeted watersheds, however, generally is older than that for water from the springs. The median age of water from wells in the Muddy Creek Watershed, for example, was 15 years, compared to 11 years for the water from the springs in that watershed, and less than 10 years for water from all springs in the spring study. The similarity in the ranges in apparent age of water from the wells and from the springs shows that the samples from the targeted watersheds and springs have bracketed the range of apparent ages that would be expected in the shallow ground-water-flow systems throughout the Chesapeake Bay Watershed. The apparent age of water from individual wells does not necessarily represent the entire distribution of ages of the discharging ground water, and it is this distribution of ages that affects the response of nutrient concentrations in stream base flow. Nutrient-reduction scenarios were modeled for two watersheds for which the distribution of apparent ground-water ages was available, the East Mahantango Creek Watershed in the Valley and Ridge siliciclastic HGMR and the Locust Grove Watershed in the Coastal Plain Uplands HGMR. A nutrient-reduction scenario was created for East Mahantango Creek, where the average residence time was determined to be approximately 10 years on the basis of the output of particle tracking from a ground-water-flow model. This scenario showed decreases of nearly 50 percent in base-flow concentrations of nitrate in streams within the first year after the reduction in nitrogen input; smaller reductions in nitrate concentration occurred in each subsequent year. A second scenario for that same watershed, in which the same 10-year average residence time was assumed and an exponential model was used for analysis, showed that a 50-percent reduction in base-flow concentrations of nitrate could take up to 5 years. For the Locust Grove Watershed, in which an average residence time of 32 years was assumed, simulation with the exponential model showed that it may take more than 20 years to achieve a 50-percent reduction in base-flow concentra-tions of nitrate. Although it was not possible to construct such scenarios for all watersheds, these examples show the range of possible responses to changes in nutrient inputs in two very different types of watersheds. Findings from this study include information on factors that affect ground-water age, spatial distribution of ages, and nitrogen transport. In the East Mahantango Creek Watershed and the Polecat Creek Watershed, the residence time varied spatially depending on the position of the flow path, and temporally depending on the recharge conditions. Generally, ground water in areas near the stream had short residence times and the water in upland areas had longer residence times. Water traveling through deep layers had longer residence times than water traveling through shallow layers, and residence times were faster under high recharge conditions than low recharge conditions. Ground water in the Pocomoke Watershed exhibits a similar pattern: younger water discharges to small order streams in headwater basins and older water discharges to larger streams near the basin outlet. Factors affecting nitrogen transport in ground water include spatial and temporal variation in input sources, ground-water age, and aquifer processes that lead to denitrification. Spatial and temporal variations in nitrogen sources affect all the watersheds. Tributaries with higher inputs of nitrogen have higher concentrations in stream base flow. Areas where nitrogen application rates have increased over time show an age-nitrate relation in ground-water samples. The age-nitrate relation can be affected by denitrification, which occurs in Pocomoke and East Mahantango Creeks but is not evident in Polecat and Muddy Creeks. In East Mahantango Creek, the level of denitrification is significant in water with residence times greater than 20 years, but because this is a small component of overall ground-water discharge to a stream, it may not remove a significant quantity of nitrogen from the system. Denitrification in Pocomoke Creek is significant and appears to affect mostly older water discharging to streams. Therefore, if most of the nitrogen entering these two streams is associated with the discharge of younger ground water, denitrification may not greatly affect the overall nitrogen delivery to these streams. Other findings of this study show that nitrate in ground water discharging along preferential flow paths may not be affected by natural processes, such as denitrification or uptake by riparian vegetation. Seeps to swales and ditches beneath the north uplands at Polecat Creek indicate a shallow water table and discharge of young ground water whereas the absence of such seeps on the south side indicates a deep water table and a lack of young ground water. Similarly, discharge at the base of the slope and to the valley wetland south of the creek but not north of the creek indicates a different role for the riparian forest on the two sides of the creek. In many of the systems where water discharges at the base of slopes to wetlands, ditches have been dug to drain the valley. Such drainage circumvents possible removal of nitrate by riparian vegetation. Because ground-water residence times do not appear directly related to the HGMRs, the targeting of management practices will achieve the most rapid response in water quality if directed at 1) watersheds with large agricultural sources of nitrate, 2) areas with the shortest ground-water-flow paths and 3) areas not affected by significant denitrification. The fastest response in stream base-flow concentrations of nitrogen to implementation of management practices would be to implement practices in those areas with the highest loads rather than attempt to target practices on the basis of HGMR stratification. Overall findings of the study indicate that 1) ground-water contributions to nitrogen in streamflow are significant, 2) some response to management practices should be evident in base-flow concentrations of nitrogen and loads within 1 to 5 years in watersheds with the shortest average residence times, but response time may be closer to 20 years in watersheds with longer average ground-water residence times, 3) the majority of the response in ground-water discharge to any changes in management practices will be distributed over a 10-year time period even in the watersheds with the fastest response times, and 4) given that half the streamflow is from ground-water discharge and the other half is runoff or soil water, about 90 percent of total water being discharged to a stream will be less than about a decade old; therefore, full implementation of nutrient reductions may result in improved streamwater quality in about a decade. In the more-likely scenario of gradual source reduction, the reduction in concentrations of nitrate in streams and aquifers would take longer than the examples shown here.

Explore related subjects

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

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Bruce D. Lindsey, Scott Phillips, Colleen A. Donnelly, Gary K. Speiran, Niel Plummer, John Karl Bohlke, Michael J. Focazio, William C. Burton, Eurybiades Busenberg. 2003. Residence times and nitrate transport in ground water discharging to streams in the Chesapeake Bay Watershed. https://doi.org/10.3133/wri034035

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

KEEP EXPLORING

Related USGS reports

Hydrogeology of the Tully Valley and characterization of mudboil activity, Onondaga County, New York

Mudboil activity in the Tully Valley, in central New York, is causing turbidity in nearby Onondaga Creek, where it has caused a bridge to collapse; it also has threatened or damaged other structures and has caused extensive land subsidence. Mudboil activity was intermittent from its first reported appearance in the 1890's until the 1970's, when the rates of mudboil discharge and land subsidence began to increase. Historically, the water discharged from mudboils was reported as fresh, but chemical analyses in the late 1970's indicated an increase in specific conductance and chloride concentration. Mudboil discharge is driven by artesian pressure in unconsolidated sediments that are confined by a 60-foot layer of silt and red clay. This process, once begun, has been self-propagating. Artesian pressures are about 20 feet above land surface over most of the valley floor but exceed 30 feet above land surface along Onondaga Creek where Rattlesnake Gulf and Rainbow Creek enter the Tully Valley. The source of artesian pressure is recharge from the Tully (Valley Heads) Moraine at the south end of the valley, and the alluvial fans of Rattlesnake Gulf and Rainbow Creek. The mudboils are found within a 300-foot-wide by 1,500-foot-long corridor along Onondaga Creek just upstream from the two alluvial fans, and in a 5-acre subsided area just west of that corridor. Remediation efforts have entailed (1) diversion of flow from the tributary that feeds the subsided area, (2) installation of depressurizing wells at several locations, and (3) construction of a dam and settling impoundment to detain mudboil sediment that would normally discharge to Onondaga Creek. These efforts have been partly successful, but further work is needed to slow the mudboil activity, which is expected to persist in both areas. Mudboil activity is normally greatest during the early spring and late fall, when artesian pressures increase in response to seasonal ground-water recharge. Suspended-sediment concentrations at the out-flow of the subsidence area ranged from 31,210 mg/L (milligrams per liter) in October 1991 to 17 mg/L after remediation efforts in the summer of 1993. Yearly average suspended-sediment loads to Onondaga Creek from the subsidence area for water years 1992, 1993, 1994, and 1995 were 29.8, 9.75, 1.41, and 1.80 tons per day, respectively. Sediment discharged from the mudboils initially was 30 to 60 percent clay and 80 to 100 percent silt-sized or smaller sediment, and the sand fraction never exceeded 20 percent. After the remediation projects, 50 to 80 percent was clay, and nearly all sediment was silt size or smaller. Analyses of water from upstream and downstream of the subsidence area, as well as from mudboil vents within that area, indicate that the source of water for some mudboils is a confined freshwater aquifer, whereas for others it is an underlying, brackish-water aquifer. Water from the freshwater aquifer has specific conductance values ranging from about 400 ( μS /cm (microsiemens per centimeter at 25° Celsius) to almost 900 ( μS /cm, dissolved chloride concentrations range from 37 to 430 mg/L, and dissolved-solids concentrations range from 215 to 463 mg/L. Specific conductance of water from the brackish-water aquifer ranges from 17,000 to 28,000 ( μS /cm, chloride concentrations range from 2,000 to 7,100 mg/L, and dissolved-solids concentrations range from 4,200 to 12,800 mg/L. The largest landslide in New York State in the last 75 years occurred at the foot of Bare Mountain, 1 mile downstream from the mudboil area, in April 1993 and was the fourth in a series of slides that have occurred at the base of this hill. Slope instability was reported as early as May 1990. After the slide, intermittent mudboil-like activity was observed at several springs within the backscarp of the slide; water from these springs ranged from fresh to brackish. The chemical similarity between water from some springs in the backscarp area and water in the lower (brackish) aquifer beneath the mudboil area may indicate a hydraulic connection between this aquifer and the surficial deposits. Hydrologic changes in the valley during the last 100 years have been attributed to salt-solution mining in the upstream (southern) end of the valley. The removal of nearly 150 feet of salt from four evaporite beds in the Syracuse Shale of the Salina Group has caused the collapse of bedrock and unconsolidated deposits in and near the brine field, 3 miles south of the mudboil area. These collapses have created a hydraulic connection among bedding plane aquifers in the bedrock and increased the hydraulic connection with unconsolidated aquifers. The ground-water flow system after brine field closure in 1988 may have reached a new semiequilibrium, but mudboil activity will likely continue because artesian pressures remain. Whether mudboils were present before salt solution-mining began is unknown.

New York

Evaluation of the Source and Transport of High Nitrate Concentrations in Ground Water, Warren Subbasin, California

Ground water historically has been the sole source of water supply for the Town of Yucca Valley in the Warren subbasin of the Morongo ground-water basin, California. An imbalance between ground-water recharge and pumpage caused ground-water levels in the subbasin to decline by as much as 300 feet from the late 1940s through 1994. In response, the local water district, Hi-Desert Water District, instituted an artificial recharge program in February 1995 using imported surface water to replenish the ground water. The artificial recharge program resulted in water-level recoveries of as much as 250 feet in the vicinity of the recharge ponds between February 1995 and December 2001; however, nitrate concentrations in some wells also increased from a background concentration of 10 milligrams per liter to more than the U.S. Environmental Protection Agency (USEPA) maximum contaminant level (MCL) of 44 milligrams per liter (10 milligrams per liter as nitrogen). The objectives of this study were to: (1) evaluate the sources of the high-nitrate concentrations that occurred after the start of the artificial-recharge program, (2) develop a ground-water flow and solute-transport model to better understand the source and transport of nitrates in the aquifer system, and (3) utilize the calibrated models to evaluate the possible effect of a proposed conjunctive-use project. These objectives were accomplished by collecting water-level and water-quality data for the subbasin and assessing changes that have occurred since artificial recharge began. Collected data were used to calibrate the ground-water flow and solute-transport models. Data collected for this study indicate that the areal extent of the water-bearing deposits is much smaller (about 5.5 square miles versus 19 square miles) than that of the subbasin. These water-bearing deposits are referred to in this report as the Warren ground-water basin. Faults separate the ground-water basin into five hydrogeologic units: the west, the midwest, the mideast, the east and the northeast hydrogeologic units. Water-quality analyses indicate that septage from septic tanks is the primary source of the high-nitrate concentrations measured in the Warren ground-water basin. Water-quality and stable-isotope data, collected after the start of the artificial recharge program, indicate that mixing occurs between imported water and native ground water, with the highest recorded nitrate concentrations in the midwest and the mideast hydrogeologic units. In general, the timing of the increase in measured nitrate concentrations in the midwest hydrogeologic unit is directly related to the distance of the monitoring well from a recharge site, indicating that the increase in nitrate concentrations is related to the artificial recharge program. Nitrate-to-chloride and nitrogen-isotope data indicate that septage is the source of the measured increase in nitrate concentrations in the midwest and the mideast hydrogeologic units. Samples from four wells in the Warren ground-water basin were analyzed for caffeine and selected human pharmaceutical products; these analyses suggest that septage is reaching the water table. There are two possible conceptual models that explain how high-nitrate septage reaches the water table: (1) the continued downward migration of septage through the unsaturated zone to the water table and (2) rising water levels, a result of the artificial recharge program, entraining septage in the unsaturated zone. The observations that nitrate concentrations increase in ground-water samples from wells soon after the start of the artificial recharge program in 1995 and that the largest increase in nitrate concentrations occur in the midwest and mideast hydrogeologic units where the largest increase in water levels occur indicate the validity of the second conceptual model (rising water levels). The potential nitrate concentration resulting from a water-level rise in the midwest and mideast hydrogeologic units was estimated using a simple mixing-cell model. The estimated value is within the range of concentrations measured in samples from wells, further indicating the validity of the second conceptual model. A ground-water flow model and a solute-transport model were developed for the Warren ground-water basin for the period 1956-2001. MODFLOW-96 was used for the ground-water flow model and MOC3D was used for the solute-transport model. The model cell size is about 500 feet by 500 feet and the models were discretized vertically into three layers. The models were calibrated using a trial-and-error approach using water-level and nitrate-concentration data collected between 1956-2001. In order to better match the measured data, low fault hydraulic characteristic values were required, thereby compartmentalizing the ground-water basin. In addition, it was necessary to parameterize the specific yield distribution for the top model layer where unconfined ground-water conditions occur into three homogeneous zones. Separate sets of specific- yield values were needed to simulate the drawdown and subsequent water-level recovery. In addition, the calibrated natural recharge was about 83 acre-feet per year. The entrainment of unsaturated-zone septage was simulated as recharge having an associated nitrate concentration. The volume of recharge was a function of the measured water-level rise between 1994-98 and the moisture content of the unsaturated zone. The nitrate concentration of the recharge water was a weighted function of the assumed nitrate concentration in the infiltrating water associated with the overlying land use. The simulated hydraulic head and nitrate concentration results were in good agreement with the measured data indicating that the mechanism for the increase in nitrate concentrations was rising water levels entraining high-nitrate septage in the unsaturated-zone. The calibrated models were used to simulate the possible effects of a planned conjunctive-use project in the western part of the ground-water basin. The simulated project included the addition of a new recharge pond and a new extraction well. In addition, recharge at two existing recharge ponds was increased and three existing production wells were pumped, treated in a nitrate-removal facility, and used for water supply. The simulated hydraulic heads increased in the west, the mideast, and parts of the east hydrogeologic units; however, the simulated hydraulic heads decreased in the midwest and northeast hydrogeologic units. The simulated nitrate concentrations increased to above the MCL of 44 milligrams per liter (10 milligrams per liter as nitrogen) in parts of the west as a result of the increase in simulated hydraulic head. The simulated nitrate concentrations decreased in part of the midwest hydrogeologic unit as a result of the artificial recharge and pumping from the nitrate-removal wells. The simulated nitrate concentrations increased to above the MCL of 44 milligrams per liter in part of the mideast and parts of the east hydrogeologic units beneath commercial land-use areas.

California

Clear-water abutment and contraction scour in the Coastal Plain and Piedmont Provinces of South Carolina, 1996-99

The U.S. Geological Survey, in cooperation with the South Carolina Department of Transportation, collected observations of clear-water aburment and contraction scour at 146 bridges in the Coastal Plain and Piedmont of South Carolina. Scour depths ranged from 0 to 23.6 feet. Theoretical scour depths were computed at each bridge and compared with observed scour. This comparison showed that theoretical scour depths, in general, exceeded the observed scour depths and often were excessive. A comparison of field data with dimensionless relations for laboratory data showed that the range of dimensionless variables used in laboratory investigations was outside of the range for field data in South Carolina, suggesting laboratory relations may not be applicable to field conditions in South Carolina. Variables determined to be important in developing scour within laboratory studies were investigated to understand their influence within the South Carolina field data, and many of these variables appeared to be insignificant under field conditions found in South Carolina. The strongest explanatory variables were embankment length, geometric-contraction ratio, approach velocity, and soil cohesion. Envelope curves developed with the field data are useful tools for assessing reasonable ranges of scour depth in South Carolina. These tools are simple to apply and are an improvement over the current methods for predicting theoretical scour. Data from this study have been compiled into a database that includes photographs, figures, observed scour depths, theoretical scour depths, limited basin characteristics, limited soil data, and theoretical hydraulic data. The database can be used to compare studied sites with unstudied sites to assess the potential for scour at the unstudied sites. In addition, the database can be used to assess the performance of various theoretical methods for predicting clear-water abutment and contraction scour.

South Carolina