USGS ScienceSearch

SEARCH · USGS Science

Results for “Maryland Geological Survey Report of Investigations”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

108 records · Page 6Linked to original sources

Development of relations of stream stage to channel geometry and discharge for stream segments simulated with Hydrologic Simulation Program-Fortran (HSPF), Chesapeake Bay Watershed and adjacent parts of Virginia, Maryland, and Delaware

The U.S. Geological Survey (USGS), U.S. Environmental Protection Agency (USEPA), Chesapeake Bay Program (CBP), Interstate Commission for the Potomac River Basin (ICPRB), Maryland Department of the Environment (MDE), Virginia Department of Conservation and Recreation (VADCR), and University of Maryland (UMD) are collaborating to improve the resolution of the Chesapeake Bay Regional Watershed Model (CBRWM). This watershed model uses the Hydrologic Simulation Program-Fortran (HSPF) to simulate the fate and transport of nutrients and sediment throughout the Chesapeake Bay watershed and extended areas of Virginia, Maryland, and Delaware. Information from the CBRWM is used by the CBP and other watershed managers to assess the effectiveness of water-quality improvement efforts as well as guide future management activities. A critical step in the improvement of the CBRWM framework was the development of an HSPF function table (FTABLE) for each represented stream channel. The FTABLE is used to relate stage (water depth) in a particular stream channel to associated channel surface area, channel volume, and discharge (streamflow). The primary tool used to generate an FTABLE for each stream channel is the XSECT program, a computer program that requires nine input variables used to represent channel morphology. These input variables are reach length, upstream and downstream elevation, channel bottom width, channel bankfull width, channel bankfull stage, slope of the floodplain, and Manning's roughness coefficient for the channel and floodplain. For the purpose of this study, the nine input variables were grouped into three categories: channel geometry, Manning's roughness coefficient, and channel and floodplain slope. Values of channel geometry for every stream segment represented in CBRWM were obtained by first developing regional regression models that relate basin drainage area to observed values of bankfull width, bankfull depth, and bottom width at each of the 290 USGS streamflow-gaging stations included in the areal extent of the model. These regression models were developed on the basis of data from stations in four physiographic provinces (Appalachian Plateaus, Valley and Ridge, Piedmont, and Coastal Plain) and were used to predict channel geometry for all 738 stream segments in the modeled area from associated basin drainage area. Manning's roughness coefficient for the channel and floodplain was represented in the XSECT program in two forms. First, all available field-estimated values of roughness were compiled for gaging stations in each physiographic province. The median of field-estimated values of channel and floodplain roughness for each physiographic province was applied to all respective stream segments. The second representation of Manning's roughness coefficient was to allow roughness to vary with channel depth. Roughness was estimated at each gaging station for each 1-foot depth interval. Median values of roughness were calculated for each 1-foot depth interval for all stations in each physiographic province. Channel and floodplain slope were determined for every stream segment in CBRWM using the USGS National Elevation Dataset. Function tables were generated by the XSECT program using values of channel geometry, channel and floodplain roughness, and channel and floodplain slope. The FTABLEs for each of the 290 USGS streamflow-gaging stations were evaluated by comparing observed discharge to the XSECT-derived discharge. Function table stream discharge derived using depth-varying roughness was found to be more representative of and statistically indistinguishable from values of observed stream discharge. Additionally, results of regression analysis showed that XSECT-derived discharge accounted for approximately 90 percent of the variability associated with observed discharge in each of the four physiographic provinces. The results of this study indicate that the methodology developed to generate FTABLEs for every s

Chesapeake Bay Watershed

Hydrogeologic controls on ground-water discharge to the Washington METRO subway tunnel near the Medical Center station and Crossover, Montgomery County, Maryland

Excessive water intrusion has been observed inside several of the Washington Metropolitan Area Transit Authority subway tunnels, with the worst leakage occurring along the Red Line tunnels and stations north of Dupont Circle in Washington, D.C. These tunnels were constructed in bedrock that contains permeable (water-bearing) joints and fractures. Excessive water leakage through the walls and water inside the underground facilities has damaged mechanical and electrical components in the tunnel, and has escalated the deterioration rate of the rail system. The U.S. Geological Survey and the Washington Metropolitan Area Transit Authority have worked cooperatively on a study from 200003 to describe and quantify the factors controlling ground-water flow into the Red Line subway tunnel near the Medical Center Station and Crossover in Montgomery County, Maryland. The Red Line near the Medical Center Station and Crossover passes through or beneath the gneissic Sykesville Formation and the biotite-hornblende tonalite member of the Georgetown Intrusive Suite, both of which contain numerous fractures. The mapped foliation and joints of the Sykesville Formation in the vicinity of the Medical Center Station and Crossover are generally orientated north-south. Fractures in the Sykesville Formation in outcrops appear to be poorly connected. In the biotite-hornblende tonalite member of the Georgetown Intrusive Suite, the general orientation of the mapped foliation and joints is east-west. In contrast to the fractures in the Sykesville Formation, the fractures in the Georgetown Intrusive Suite in outcrops appear to be more numerous and have a greater degree of connectivity. Fractures intersecting four bedrock wells near the Medical Center Station and Crossover that were drilled into the biotite-hornblende tonalite member of the Georgetown Intrusive Suite show an east-west orientation matching the foliation and joints shown on geologic maps. The excessive water intrusion at the Medical Center Station and Crossover could be the result of its location within the Georgetown Intrusive Suite. The abrupt changes in the mapped directions of ground-water flow based on the hydraulic heads at the contact between the Sykesville Formation and biotite-hornblende tonalite member of the Georgetown Intrusive Suite could also be the result of the change in fracturing between these two lithologies. Saprolite, a residual of soft, red/brown to gray clay from decomposed crystalline rock, overlies the bedrock and varies from about 20 to 55 feet thick, depending on location. On the basis of a slug test conducted in the lower part of the saprolite near the Medical Center Station and Crossover, transmissivity and storativity of the saprolite were estimated to be 10 feet squared per day and 10-6 , respectively. The transmissivity of fractures intersecting bedrock boreholes drilled in the biotite-hornblende tonalite member of the Georgetown Intrusive Suite varies over five orders of magnitude, from a maximum of approximately 10 feet squared per day to the detection limit of the in situ testing apparatus, which is approximately 10-4 feet squared per day. In general, the transmissivity of fractures intersecting the boreholes increases with depth. The low transmissivity of bedrock fractures in close proximity to the saprolite is likely to be caused by the fractures being filled with byproducts of rock weathering, resulting in reduced permeability. 2 Hydrogeologic Controls on Ground-Water Discharge to the Washington METRO Subway Tunnel The bulk transmissivity of the bedrock aquifer is approximately 3.7 feet squared per day, as determined from an aquifer test conducted by pumping a 240-foot-deep borehole and monitoring the drawdown over 3 days in the pumped borehole and several observation boreholes. In general, the hydraulic head decreases with depth in bedrock boreholes, indicating the potential for downward ground-water flow. Based on hydraulic head values mea

Water-Resources Investigations Report

Microbial consortia development and microcosm and column experiments for enhanced bioremediation of chlorinated volatile organic compounds, West Branch Canal Creek wetland area, Aberdeen Proving Ground, Maryland

Chlorinated solvents, including 1,1,2,2-tetrachloroethane, tetrachloroethene, trichloroethene, carbon tetrachloride, and chloroform, are reaching land surface in localized areas of focused ground-water discharge (seeps) in a wetland and tidal creek in the West Branch Canal Creek area, Aberdeen Proving Ground, Maryland. In cooperation with the U.S. Army Garrison, Aberdeen Proving Ground, Maryland, the U.S. Geological Survey is developing enhanced bioremediation methods that simulate the natural anaerobic degradation that occurs without intervention in non-seep areas of the wetland. A combination of natural attenuation and enhanced bioremediation could provide a remedy for the discharging ground-water plumes that would minimize disturbance to the sensitive wetland ecosystem. Biostimulation (addition of organic substrate or nutrients) and bioaugmentation (addition of microbial consortium), applied either by direct injection at depth in the wetland sediments or by construction of a permeable reactive mat at the seep surface, were tested as possible methods to enhance anaerobic degradation in the seep areas. For the first phase of developing enhanced bioremediation methods for the contaminant mixtures in the seeps, laboratory studies were conducted to develop a microbial consortium to degrade 1,1,2,2-tetrachloroethane and its chlorinated daughter products under anaerobic conditions, and to test biostimulation and bioaugmentation of wetland sediment and reactive mat matrices in microcosms. The individual components required for the direct injection and reactive mat methods were then combined in column experiments to test them under groundwater- flow rates and contaminant concentrations observed in the field. Results showed that both direct injection and the reactive mat are promising remediation methods, although the success of direct injection likely would depend on adequately distributing and maintaining organic substrate throughout the wetland sediment in the seep area. For bioaugmentation, two mixed anaerobic cultures, named the 'West Branch Consortia' (WBC-1 and WBC-2), were developed by enrichment of wetland sediment collected from two contaminated sites in the study area where rapid and complete reductive dechlorination naturally occurs. WBC are capable of degrading 1,1,2,2-tetrachloroethane, 1,1,2-trichloroethane, 1,2-dichloroethane, tetrachloroethene, trichloroethene, cis- and trans-1,2-dichloroethene, and vinyl chloride to the non-chlorinated end-products ethene and ethane. In addition, the column experiments showed that the consortia could completely degrade carbon tetrachloride and chloroform, although they were not grown on these contaminants. No other cultures are known that can degrade the broad mixture of chlorinated alkanes, alkenes, and methanes as shown for WBC. WBC-2 (suspended in the culture media) is capable of complete dechlorination of 50 micromolar 1,1,2,2-tetrachloroethane to ethene in 1 to 2 days with little transient accumulation of chlorinated daughter products. Only about 5 percent of the clones sequenced from WBC-1 and WBC-2 were related to dechlorinating bacteria that have been studied previously in culture, indicating the presence of unknown dechlorinators. Dehalococcoides spp. comprised about 1 percent of WBC-1 and WBC-2, which is minor compared to the population size of about 30 percent in other dechlorinating consortia for chlorinated alkenes. Although both WBC-1 and WBC-2 showed efficient degradation in laboratory tests in this study, long-term cultivation of WBC-1, which was developed using hydrogen as the organic substrate, was determined to be infeasible. Thus, WBC-2, cultivated with lactate as the organic substrate, would be used in future tests. Nutrient (ammonia and phosphate mixture) addition to anaerobic microcosms constructed with wetland sediment and ground water collected from the study area showed some enhancement in the degradation rate of 1,1,2,2-tetrachloroethane, but degrada

Maryland

Water quality of groundwater and stream base flow in the Marcellus Shale Gas Field of the Monongahela River Basin, West Virginia, 2011-12

The Marcellus Shale gas field underlies portions of New York, Pennsylvania, Ohio, Virginia, Maryland, Tennessee, and West Virginia. Development of hydraulic fracturing and horizontal drilling technology led to extensive development of gas from the Marcellus Shale beginning about 2007. The need to identify and monitor changes in water-quality conditions related to development of the Marcellus Shale gas field prompted the U.S. Geological Survey, in cooperation with the West Virginia Department of Environmental Protection, Division of Water and Waste Management, to document water quality for comparison with water quality in samples collected at a future date. The identification of change in water-quality conditions over time is more difficult if baseline water-quality conditions have not been documented. U.S. Geological Survey personnel sampled groundwater and surface water in West Virginia’s Monongahela River Basin during 2011–12. A groundwater survey, in which 39 wells and 2 springs were sampled, was conducted during June through September 2011. A base-flow survey was conducted during July through October 2012; 50 stream sites were sampled under base-flow conditions in this survey. Because additives to hydraulic fracturing fluids are variable and decrease in flowback water over a relatively short time, water-quality analyses for this study focused on documenting the water-quality characteristics typical of water from shallow aquifers; water derived from contact with the Marcellus Shale (flowback from hydraulic fracturing or formation water); and water with constituents from conventional oil and gas development, sewage effluent, and coal-mine drainage. All samples were analyzed for field properties (water temperature, pH, specific conductance, dissolved oxygen, and turbidity), major ions, trace elements, naturally occurring radioactive materials, and stable isotopes. In addition to documenting baseline water-quality conditions for an area of shale-gas development, these data were examined for patterns in water quality. Groundwater and base-flow survey data were compared to historical data from the Monongahela River Basin in West Virginia. Additionally, groundwater- and base-flow survey samples were grouped by Marcellus Shale gas production in the subbasin in which that sampling site was located. The comparisons of data collected as part of this study with historical data identified few differences. No significant difference was found in a comparison of groundwater survey data and historical data. Base-flow survey samples differed significantly from historical data for pH, chloride, and strontium, all of which had higher concentrations in the base-flow survey samples. Differences in pH are likely related to changes in mining regulation beginning in 1977. Concentrations of chloride and strontium elevated above background concentrations may be related to saline groundwater; saline water is within 300 feet of the land surface in parts of the study area. In the comparison of base-flow survey samples grouped by shale-gas-production setting, significant differences were found for fluoride and barium. Concentrations of fluoride and barium were higher in stream subbasins with active Marcellus Shale production than in subbasins not near active Marcellus Shale production. Elevated fluoride and barium are associated with deep brines. Generally, naturally occurring radioactive materials were not found in elevated concentrations in either groundwater or base-flow samples. Only 3 samples, 2 from the groundwater survey and one from the base-flow survey, exceeded the U.S. Environmental Protection Agency maximum contaminant level for radium isotopes of 5.0 picocurie per liter for either a single isotope or a combined value of radium-226 and radium-228. Stable isotope composition indicates broad similarity among surface water, shallow groundwater, and precipitation in the region. Neither shallow groundwater nor surface water showed a marked similarity with the deep brines associated with shale gas. In most of the groundwater survey samples, 38 of 41 samples, dissolved gas profiles were similar to those previously found in samples from shallow, domestic wells in the region. This study provides a baseline of water-quality conditions in the Monongahela River Basin in West Virginia during the early phases of development of the Marcellus Shale gas field. Although not all inclusive, the results of this study provide a set of reliable water-quality data against which future data sets can be compared and the effects of shale-gas development may be determined.

West Virginia

Hydrogeologic framework and hydrologic conditions of the Piney Point aquifer in Virginia

The Piney Point aquifer in Virginia is newly described and delineated as being composed of six geologic units, in a study conducted by the U.S. Geological Survey in cooperation with the Virginia Department of Environmental Quality (VA DEQ). The eastward-dipping geologic units include, in stratigraphically ascending order, the sand of the Nanjemoy Formation Woodstock Member, interbedded limestone and sand of the Piney Point Formation, silty and clayey sand of the Gosport Formation equivalent sediments, silty sand of the Oligocene-age sediments, silty fine-grained sand of the Old Church Formation, and silty sand of the Calvert Formation, Newport News unit and basal Plum Point Member. Identification of geologic units is based on typical sediment lithologies of geologic formations. Fine-grained sediments that compose confining units positioned immediately above and below the Piney Point aquifer are also described. The Piney Point aquifer is one of several confined aquifers within the Virginia Coastal Plain and includes a highly porous and solution-channeled indurated limestone within the Piney Point Formation from which withdrawals are made. The limestone is relatively continuous laterally across central parts of the Northern Neck, Middle Peninsula, and York-James Peninsula. Other geologic units are of variable extent. The configurations of most of the geologic units are further affected by newly identified faults that are aligned radially from the Chesapeake Bay impact crater and create constrictions or barriers to groundwater flow. Some geologic units are also truncated beneath the lower Rappahannock River by a resurge channel associated with the impact crater. Groundwater withdrawals from the Piney Point aquifer increased from approximately 1 million gallons per day (Mgal/d) during 1900 to 7.35 Mgal/d during 2004. As a result, a water-level cone of depression in James City and northern York Counties was estimated to be as low as 70 feet (ft) below the National Geodetic Vertical Datum of 1929 (NGVD 29) by 2005. Withdrawals decreased to 5.01 Mgal/d by 2009 as withdrawals were shifted toward other sources, and by 2015 water levels had recovered to approximately 50 ft below NGVD 29. The mean estimated transmissivity of the Piney Point aquifer in York and James City Counties is 16,300 feet squared per day (ft 2 /d), but farther north it is only 925 ft 2 /d. The mean well specific capacity in York and James City Counties is 11.4 gallons per minute per foot (gal/min/ft). Farther north in Virginia, mean specific capacity is only 2.26 gal/min/ft, and in Maryland it is 0.99 gal/min/ft. The northward decrease in specific capacity probably reflects the northward decrease in transmissivity, which results from poor development of the solution-channeled limestone. An aquifer test in northern York County induced vertical leakage to the solution-channeled limestone from overlying silty sand and a change in response of the aquifer to pumping from a single layer to two layers. Transmissivity of the limestone of approximately 19,800 ft 2 /d was distinguished from the silty sand of approximately 2,500 ft 2 /d. Most of the water in the Piney Point aquifer is slightly alkaline with moderate concentrations primarily of sodium and bicarbonate that are slightly undersaturated with respect to calcite. Iron concentrations are generally less than 0.3 milligrams per liter (mg/L). Mixing of freshwater with seawater has elevated chloride concentrations to the southeast to as much as 7,120 mg/L. Information on the Piney Point aquifer can benefit water-resource management in siting production wells, predicting likely well yield, and anticipating water-level response to withdrawals. Models that vertically discretize individual geologic units can potentially be used to evaluate groundwater flow in greater detail by representing lateral flow and vertical leakage among the geologic units. Because groundwater withdrawals are made primarily from the limestone and sand of the Piney Point Formation, the VA DEQ has considered regarding the limestone and sand singly as a regulated aquifer apart from the other geologic units. Under current policy in Virginia, if only the limestone and sand were regarded as a regulated aquifer, a greater amount of drawdown would be allowed than is allowed for the Piney Point aquifer consisting of six geologic units. Some production wells intercept multiple geologic units, and the units can undergo water-level decline and vertical leakage induced by pumping from the limestone and sand. Whether the other geologic units are to be regarded as regulated aquifers is an additional consideration for the VA DEQ.

Virginia

Geomorphic characteristics of Tenmile Creek, Montgomery County, Maryland, 2014–16

Data collected from April 2014 through September 2016 were used to assess geomorphic characteristics and geomorphic changes over time in a selected reach of Tenmile Creek, a small rural watershed near Clarksburg, Maryland. Longitudinal profiles of the channel bed, water surface, and bank features were developed from field surveys. Changes in cross-section geometry between field surveys were documented. Grain-size distributions for the channel bed were developed from pebble counts. Continuous-record streamflow and precipitation data were also collected in the Tenmile Creek watershed and used to supplement the geomorphic analyses. The Rosgen system of stream classification was used to classify the stream channel according to morphological measurements of slope, entrenchment ratio, width-to-depth ratio, sinuosity, and median particle diameter of the channel materials. Boundary shear stress near the U.S. Geological Survey (USGS) streamflow-gaging station was assessed by using hydraulic variables computed from the cross-section surveys and slope measurements derived from crest-stage gages and temporary data loggers installed along the study reach. Analysis of the longitudinal profiles indicated relatively small changes in the percentage and distribution of riffles, pools, and runs in the study reach between April 2014 and March 2015. More noticeable changes were observed during surveys conducted in March 2016 and September 2016. The channel-bed slope showed a net reduction over time from 0.0072 to 0.0040 feet per foot (ft/ft). The low-flow water-surface slope also showed a net reduction over time from 0.0065 to 0.0045 ft/ft. Net aggradation in the lower section of the study reach combined with net degradation in the upper section of the study reach contributed to the net reduction in channel-bed and water-surface slope. The large storm and resulting flood on July 30, 2016 was a major factor in observed changes in the longitudinal profiles between the March 2016 and September 2016 surveys. Comparison of data from the cross-sectional surveys indicated vertical changes in all cross sections, with more extreme changes observed between surveys in the lower section of the study reach due in part to alternating periods of net storage and transport of sand. Lateral erosion was not a major factor in the study reach, with the exception of cross section Dd, where considerable lateral erosion was documented during the study period. The flood that resulted from the large storm on July 30, 2016 was a major factor in some of the vertical changes observed in the channel bed of the study reach cross sections. Particle-size analyses of the channel bed from pebble counts indicated median particle diameters ranging from 15.5 millimeters (mm) to 23.1 mm, which is characterized as medium to coarse gravel. Sand percentages ranging from 3.4 percent to 16.4 percent of the total counts were observed over time. Net increases in storage of fine sediment in the reach were observed between April 2014 and March 2016, and a considerable reduction in storage was observed between March 2016 and September 2016. The Tenmile Creek stream channel was classified as a C4 channel, based on morphological descriptions from the Rosgen system of stream classification. The C4 classification describes a single-thread channel with a slight entrenchment ratio; a moderate to high width-to-depth ratio; moderate to high sinuosity; a water-surface slope of less than 2 percent; and a median particle diameter in the gravel range of 2 to 64 mm. The analysis of boundary shear stress indicated a range of 0.35 to 1.18 pounds per square foot for instantaneous streamflow ranging from 79 to 2,860 cubic feet per second during the study period. The relation between discharge and boundary shear stress for Tenmile Creek was compared to similar relations that were previously developed for Minebank Run, a small, urban watershed in the eastern section of the Piedmont Physiographic Province in Baltimore County, Md. that was physically restored during 2004–05. The comparison indicated a much flatter slope in the relation for Minebank Run in both its unrestored and restored conditions. This difference in the relations indicates that the erosive power in the urban watershed of Minebank Run is much more sensitive to increases in discharge magnitude than in the non-urban watershed of Tenmile Creek.

Maryland

The water-quality monitoring program for the Baltimore reservoir system, 1981-2007—Description, review and evaluation, and framework integration for enhanced monitoring

The City of Baltimore, Maryland, and parts of five surrounding counties obtain their water from Loch Raven and Liberty Reservoirs. A third reservoir, Prettyboy, is used to resupply Loch Raven Reservoir. Management of the watershed conditions for each reservoir is a shared responsibility by agreement among City, County, and State jurisdictions. The most recent (2005) Baltimore Reservoir Watershed Management Agreement (RWMA) called for continued and improved water-quality monitoring in the reservoirs and selected watershed tributaries. The U.S. Geological Survey (USGS) conducted a retrospective review of the effectiveness of monitoring data obtained and analyzed by the RWMA jurisdictions from 1981 through 2007 to help identify possible improvements in the monitoring program to address RWMA water-quality concerns. Long-term water-quality concerns include eutrophication and sedimentation in the reservoirs, and elevated concentrations of (a) nutrients (nitrogen and phosphorus) being transported from the major tributaries to the reservoirs, (b) iron and manganese released from reservoir bed sediments during periods of deep-water anoxia, (c) mercury in higher trophic order game fish in the reservoirs, and (d) bacteria in selected reservoir watershed tributaries. Emerging concerns include elevated concentrations of sodium, chloride, and disinfection by-products (DBPs) in the drinking water from both supply reservoirs. Climate change and variability also could be emerging concerns, affecting seasonal patterns, annual trends, and drought occurrence, which historically have led to declines in reservoir water quality. Monitoring data increasingly have been used to support the development of water-quality models. The most recent (2006) modeling helped establish an annual sediment Total Maximum Daily Load to Loch Raven Reservoir, and instantaneous and 30-day moving average water-quality endpoints for chlorophyll-a (chl-a) and dissolved oxygen (DO) in Loch Raven and Prettyboy Reservoirs. Modelers cited limitations in data, including too few years with sufficient stormflow data, and (or) a lack of (readily available) data, for selected tributary and reservoir hydrodynamic, water-quality, and biotic conditions. Reservoir monitoring also is too infrequent to adequately address the above water-quality endpoints. Monitoring data also have been effectively used to generally describe trophic states, changes in trophic state or conditions related to trophic state, and in selected cases, trends in water-quality or biotic parameters that reflect RWMA water-quality concerns. Limitations occur in the collection, aggregation, analyses, and (or) archival of monitoring data in relation to most RWMA water-quality concerns. Trophic, including eutrophic, conditions have been broadly described for each reservoir in terms of phytoplankton production, and variations in production related to typical seasonal patterns in the concentration of DO, and hypoxic to anoxic conditions, where the latter have led to elevated concentrations of iron and manganese in reservoir and supply waters. Trend analyses for the period 1981-2004 have shown apparent declines in production (algal counts and possibly chl-a). The low frequency of phytoplankton data collection (monthly or bimonthly, depending on the reservoir), however, limits the development of a model to quantitatively describe and relate temporal variations in phytoplankton production including seasonal succession to changes in trophic states or other reservoir water-quality or biotic conditions. Extensive monitoring for nutrients, which, in excessive amounts, cause eutrophic conditions, has been conducted in the watershed tributaries and reservoirs. Data analyses (1980-90s) have (a) identified seasonal patterns in concentrations, (b) characterized loads from (non)point sources, and (c) shown that different seasonal patterns and trends in nutrient concentrations occur between watershed tributaries and downstream reservoir.

Maryl;Pennsylvania

Natural attenuation of chlorinated volatile organic compounds in a freshwater tidal wetland, Aberdeen Proving Ground, Maryland

Ground-water contaminant plumes that are flowing toward or currently discharging to wetland areas present unique remediation problems because of the hydrologic connections between ground water and surface water and the sensitive habitats in wetlands. Because wetlands typically have a large diversity of microorganisms and redox conditions that could enhance biodegradation, they are ideal environments for natural attenuation of organic contaminants, which is a treatment method that would leave the ecosystem largely undisturbed and be cost effective. During 1992-97, the U.S. Geological Survey investigated the natural attenuation of chlorinated volatile organic compounds (VOC's) in a contaminant plume that discharges from a sand aquifer to a freshwater tidal wetland along the West Branch Canal Creek at Aberdeen Proving Ground, Maryland. Characterization of the hydrogeology and geochemistry along flowpaths in the wetland area and determination of the occurrence and rates of biodegradation and sorption show that natural attenuation could be a feasible remediation method for the contaminant plume that extends along the West Branch Canal Creek. The aquifer, which received contaminants in the past from sources that were located upgradient of the current eastern boundary of the wetland, is about 40 to 45 feet thick in the study area. The overlying wetland sediments consist of two distinct layers that have a combined thickness of about 6 to 12 feet--an upper unit of peat and a lower unit of silty to sandy clay or clayey sand. Head distributions show strong vertical gradients, and flow directions are predominantly upward through the wetland sediments. Tidal fluctuations, however, cause some reversals in ground-water-flow directions and changes in discharge locations, which affect the distribution and transport of contaminants. The average linear velocity of ground water is estimated to be about 2 to 3 feet per year along vertical flow lines in the wetland area, and the total ground-water discharge along a 1-foot-wide strip of the wetland extending from the eastern wetland boundary to the creek is in the range of 0.13 to 0.25 feet squared per day. The major parent contaminants in the aquifer were trichloroethylene (TCE); 1,1,2,2-tetrachloroethane (PCA); carbon tetrachloride (CT); and chloroform (CF). The aquifer was typically aerobic, but ground water in the wetland sediments became increasingly anaerobic along the upward flow direction. Iron-reducing conditions were predominant in the lower wetland sediment unit composed of clayey sand and silt, and methanogenesis was predominant in an upper unit composed of peat. Total concentrations of VOC's and the relative proportions of parent compounds to anaerobic daughter products changed substantially as the contaminants were transported upward through these changing redox environments. Concentrations of the parent compounds TCE and PCA ranged from about 100 to 2,000 micrograms per liter in the aquifer beneath the wetland, whereas concentrations of daughter products were low or undetectable. In contrast, the parent compounds commonly were not detected in ground water in the wetland sediments, but the daughter compounds 1,2- dichloro-ethylene (12DCE), vinyl chloride (VC), 1,1,2-trichloroethane (112TCA), and 1,2-dichloroethane (12DCA) were observed. 12DCE and VC were the dominant daughter products in the wetland sediments, even where PCA was the primary parent contaminant discharging to the wetland. The presence of these daughter products in water in the wetland sediments indicates that TCE and PCA are degraded by reductive dechlorination reactions in the naturally anaerobic wetland sediments. Although production of daughter products was observed, total concentrations of the VOC's decreased upward through the 6- to 12-feet-thick wetland sediments and were less than 5 micrograms per liter near the surface. The daughter products are apparently further degraded by these anaerobic processes to non-chlorinated, non-toxic endproducts, or are removed by other attenuation processes such as aerobic degradation or volatilization. Several sets of anaerobic microcosms, some amended with TCE concentrations of 400 or 990 micrograms per liter (3.0 or 7.5 micromoles per liter) and one amended with PCA concentration of 480 micrograms per liter (2.9 micromoles per liter), were constructed in the laboratory using wetland sediment and ground water from the study area. These microcosm experiments confirmed field observations that 12DCE and VC are the dominant daughter products from anaerobic biodegradation of both TCE and PCA. Production of 112TCA and 12DCA also was observed in the PCA-amended microcosms, but concentrations were lower than those of 12DCE and VC. These results indicate that degradation of PCA occurs through both hydrogenolysis and dichloroelimination pathways under anaerobic conditions. Daughter products were not observed in sterile controls, except for relatively low concentrations of 12DCE in the PCA-amended controls, suggesting that the reactions were predominantly microbially mediated. Parent and daughter concentrations in the microcosms decreased to less than 5 micrograms per liter in less than 34 days under methanogenic conditions, thus showing extremely rapid biodegradation rates in these organic-rich wetland sediments. Maximum potential rate constants for biodegradation of TCE and PCA, which were calculated from the microcosm exper-iments assuming first-order degradation rates, ranged from 0.10 to 0.31 per day under methanogenic conditions, corresponding to half-lives of 2 to 7 days. The rate constant for TCE biodegradation under sulfate-reducing conditions was 0.045 West Branch, Aberdeen Proving Ground, Maryland 3 per day (half-life of 15 days), which is lower than under methanogenic conditions. These estimated rate constants for the wetland sediments are two to three orders of magnitude higher than those reported in the literature for TCE biodegradation in microcosms constructed with sandy aquifer sediments. Aerobic biodegradation rates for cis-12DCE, trans-12DCE, and VC were in the same range as those measured for TCE and PCA under anaerobic conditions in microcosm experiments. Thus, production of these daughter products by anaerobic biodegradation of TCE and PCA could be balanced by their consumption where oxygen is available in the wetland sediment, such as near roots and land surface. In the aerobic microcosm experiments, biodegradation of cis-12DCE, trans-12DCE, and VC occurred only if methane consumption occurred, indicating that methanotrophs were involved in the process. Aerobic biodegradation was fastest for vinyl chloride and slowest for TCE in the microcosm experiments. These results agree with other laboratory and field studies that have shown faster degradation by methane-utilizing cultures when the compounds are less halogenated. Equilibrium sorption isotherms were measured in 24-hour batch tests and used to estimate distribution coefficients ( K d 's) to describe the ratios of sorbed to aqueous concentrations of PCA and the daughter products cis-12DCE, trans-12DCE, and VC. The estimated K d's for PCA, cis-12DCE, trans-12DCE, and VC were about 2.3, 1.8, 2.4, and 1.3 liters per kilogram of sediment, respectively. Sorbed concentrations of PCA, cis-12DCE, and trans-12DCE in the wetland sediments, therefore, would be expected to be about twice the concentration measured in the water, whereas sorbed concentrations of VC would not be much greater than the aqueous concentrations. Coefficients of retardation, which were calculated using the K d 's and an advective flow velocity of about 2 feet per year, indicate that sorption alone would cause the movement of the contaminants in the wetland sediments to be 6 to 10 times slower than the advective ground-water flow. Biodegradation and sorption, therefore, are significant natural attenuation mechanisms for chlorinated hydrocarbons in these wetland sediments. The relatively thin layers of wetland sediments are critical in reducing contaminant concentrations and toxicity of the ground water before it discharges to the wetland surface and the creek, and natural attenuation in the wetland sediments could be an effective remediation method for the ground-water contaminants.

Maryland

Water-level conditions in selected confined aquifers of the New Jersey and Delaware coastal plain, 2003

The Coastal Plain aquifers of New Jersey provide an important source of water for more than 2 million people. Steadily increasing withdrawals from the late 1800s to the early 1990s resulted in declining water levels and the formation of regional cones of depression. In addition to decreasing water supplies, declining water levels in the confined aquifers have led to reversals in natural hydraulic gradients that have, in some areas, induced the flow of saline water from surface-water bodies and adjacent aquifers to freshwater aquifers. In 1978, the U.S. Geological Survey began mapping the potentiometric surfaces of the major confined aquifers of New Jersey every 5 years in order to provide a regional assessment of ground-water conditions in multiple Coastal Plain aquifers concurrently. In 1988, mapping of selected potentiometric surfaces was extended into Delaware. During the fall of 2003, water levels measured in 967 wells in New Jersey, Pennsylvania, northeastern Delaware, and northwestern Maryland were used estimate the potentiometric surface of the principal confined aquifers in the Coastal Plain of New Jersey and five equivalent aquifers in Delaware. Potentiometric-surface maps and hydrogeologic sections were prepared for the confined Cohansey aquifer of Cape May County, the Rio Grande water-bearing zone, the Atlantic City 800-foot sand, the Vincentown aquifer, and the Englishtown aquifer system in New Jersey, as well as for the Piney Point aquifer, the Wenonah-Mount Laurel aquifer, and the Upper Potomac-Raritan-Magothy, the Middle and undifferentiated Potomac-Raritan-Magothy, and the Lower Potomac-Raritan-Magothy aquifers in New Jersey and their equivalents in Delaware. From 1998 to 2003, water levels in many Coastal Plain aquifers in New Jersey remained stable or had recovered, but in some areas, water levels continued to decline as a result of pumping. In the Cohansey aquifer in Cape May County, water levels near the center of the cone of depression underlying the southern part of the peninsula remained about the same as in 1998. To the south, recoveries up to 8 feet were observed in southern Lower Township as withdrawals had decreased since 1998. In the northern part of Cape May County, water levels had not changed substantially from historic conditions. In the Rio Grande water-bearing zone, water levels rose by as much as 13 ft at the Rio Grande well field; elsewhere across the aquifer, little change had occurred. In the Atlantic City 800-foot sand, water-level changes were greatest in southern Cape May County; at the Cape May desalination wells, water levels were as much as 32 ft lower in 2003 than in 1998. In contrast, water levels at the center of a regional cone of depression near Atlantic City rose by as much as 10 ft. Within the Piney Point aquifer water levels rose by 46 ft near Seaside Park. Similarly, water levels increased by more than 30 ft in and around the major cone of depression underlying Dover, Delaware. In the Vincentown aquifer, water levels stabilized or recovered by 2 ft to 6 ft from 1998 to 2003 in most of the wells measured; the exception is near Adelphia in Monmouth County, where water levels rose by as much as 18 ft. From 1998 to 2003, water levels near the center of a large cone of depression that extends from Monmouth to Ocean County recovered by as much as 20 ft in the Wenonah-Mount Laurel aquifer. Concurrently, ground-water levels within the Englishtown aquifer system declined by as much as 13 ft in the same area. Water levels across much of the Upper Potomac-Raritan-Magothy aquifer in the northern Coastal Plain remained about the same as 5 years previous, except in northern Ocean County where ground-water levels declined 10 ft to 33 ft. Water levels in the Middle Potomac-Raritan-Magothy aquifer declined from 5 to 9 ft along the border between Monmouth and Middlesex County. Elsewhere, across the northern part of the Coastal Plain, water levels stabilized within the Cretaceous-a

Delaware, New Jersey

Estimation of selected streamflow statistics for a network of low-flow partial-record stations in areas affected by Base Realignment and Closure (BRAC) in Maryland

The U.S. Geological Survey, in cooperation with the Maryland Department of the Environment, operated a network of 20 low-flow partial-record stations during 2008 in a region that extends from southwest of Baltimore to the northeastern corner of Maryland to obtain estimates of selected streamflow statistics at the station locations. The study area is expected to face a substantial influx of new residents and businesses as a result of military and civilian personnel transfers associated with the Federal Base Realignment and Closure Act of 2005. The estimated streamflow statistics, which include monthly 85-percent duration flows, the 10-year recurrence-interval minimum base flow, and the 7-day, 10-year low flow, are needed to provide a better understanding of the availability of water resources in the area to be affected by base-realignment activities. Streamflow measurements collected for this study at the low-flow partial-record stations and measurements collected previously for 8 of the 20 stations were related to concurrent daily flows at nearby index streamgages to estimate the streamflow statistics. Three methods were used to estimate the streamflow statistics and two methods were used to select the index streamgages. Of the three methods used to estimate the streamflow statistics, two of them--the Moments and MOVE1 methods--rely on correlating the streamflow measurements at the low-flow partial-record stations with concurrent streamflows at nearby, hydrologically similar index streamgages to determine the estimates. These methods, recommended for use by the U.S. Geological Survey, generally require about 10 streamflow measurements at the low-flow partial-record station. The third method transfers the streamflow statistics from the index streamgage to the partial-record station based on the average of the ratios of the measured streamflows at the partial-record station to the concurrent streamflows at the index streamgage. This method can be used with as few as one pair of streamflow measurements made on a single streamflow recession at the low-flow partial-record station, although additional pairs of measurements will increase the accuracy of the estimates. Errors associated with the two correlation methods generally were lower than the errors associated with the flow-ratio method, but the advantages of the flow-ratio method are that it can produce reasonably accurate estimates from streamflow measurements much faster and at lower cost than estimates obtained using the correlation methods. The two index-streamgage selection methods were (1) selection based on the highest correlation coefficient between the low-flow partial-record station and the index streamgages, and (2) selection based on Euclidean distance, where the Euclidean distance was computed as a function of geographic proximity and the basin characteristics: drainage area, percentage of forested area, percentage of impervious area, and the base-flow recession time constant, t. Method 1 generally selected index streamgages that were significantly closer to the low-flow partial-record stations than method 2. The errors associated with the estimated streamflow statistics generally were lower for method 1 than for method 2, but the differences were not statistically significant. The flow-ratio method for estimating streamflow statistics at low-flow partial-record stations was shown to be independent from the two correlation-based estimation methods. As a result, final estimates were determined for eight low-flow partial-record stations by weighting estimates from the flow-ratio method with estimates from one of the two correlation methods according to the respective variances of the estimates. Average standard errors of estimate for the final estimates ranged from 90.0 to 7.0 percent, with an average value of 26.5 percent. Average standard errors of estimate for the weighted estimates were, on average, 4.3 percent less than the best average standard errors of estima

Scientific Investigations Report

Water-quality assessment of the Lower Susquehanna River Basin, Pennsylvania and Maryland: sources, characteristics, analysis and limitations of nutrient and suspended-sediment data, 1975-90

This report describes analyses of available information on nutrients and suspended sediment collected in the Lower Susquehanna River Basin during water years 1975-90. Most of the analyses were applied to data collected during water years 1980-89. The report describes the spatial and temporal availability of nutrient and suspended-sediment data and presents a preliminary concept of the spatial and temporal patterns of concentrations and loads within the basin. Where data were available, total and dissolved forms of nitrogen and phosphorus species from precipitation, surface water, ground water, and springwater, and bottom material from streams and reservoirs were evaluated. Suspended-sediment data from streams also were evaluated. The U.S. Geological Survey National Water Information System (NWIS) database was selected as the primary database for the analyses. Precipitation-quality data from the National Atmospheric Deposition Program (NADP) and bottom-material-quality data from the National Uranium Resource Evaluation (NURE) were used to supplement the water-quality data from NWIS. Concentrations of nutrients were available from 3 precipitation sites established for longterm monitoring purposes, 883 wells (854 synoptic areal survey sites and 29 project and research sites), 23 springs (17 synoptic areal survey sites and 6 project and research sites), and 894 bottom-material sites (840 synoptic areal survey sites and 54 project and research sites). Concentrations of nutrients and (or) suspended sediment were available from 128 streams (36 long-term monitoring sites, 51 synoptic areal survey sites, and 41 project and research sites). Concentrations of nutrients and suspended sediment in streams varied temporally and spatially and were related to land use, agricultural practices, and streamflow. A general north-to-south pattern of increasing median nitrate concentrations, from 2 to 5 mg/L, was detected in samples collected in study unit streams. In streams that drain areas dominated by agriculture, concentrations of nutrients and suspended sediment tend to be elevated with respect to those found in areas of other land-use types and are related to the amount of commercial fertilizer and animal manure applied to the area drained by the streams. Animal manure is the dominant source of nitrogen for the streams in the lower, agricultural part of the basin. Concentrations of nutrients in samples from wells varied with season and well depth and were related to hydrogeologic setting. Median concentrations of nitrate were 2.5 and 3.5 mg/L for wells drawing water at depths of 0 to 100 ft and 101 to 200 ft, respectively. The lowest median concentrations for nitrate in ground water from wells were generally found in siliciclastic-bedrock, forested settings of the Ridge and Valley Physiographic Province, and the highest were found in carbonate-bedrock agricultural settings of the Piedmont Physiographic Province. Twenty-five percent of the measurements from wells in carbonate rocks in the Piedmont Physiographic Province exceeded the Pennsylvania drinking-water standard. An estimate of mass balance of nutrient loads within the Lower Susquehanna River Basin was produced by combining the available information on stream loads, atmosphericdeposition loads, commercial-fertilizer applications, animal-manure production, privateseptic-system nonpoint-source loads, and municipal and industrial point-source loads. The percentage of the average annual nitrate load carried in base flow of streams in the study unit ranged from 45 to 76 percent, and the average annual phosphorus load carried in base flow ranged from 20 to 33 percent. Average annual yields of nutrients and suspended sediment from tributary basins are directly related to percentage of drainage area in agriculture and inversely to drainage area. Information required to compute loads of nitrogen and phosphorus were available for all sources except atmospheric deposition, for which only nitrogen data were available. Atmospheric deposition is the dominant source of nitrogen for the mostly forested basins draining the upper half of the study unit. The estimate of total annual nitrogen load to the study unit from precipitation is 98.8 million pounds. Nonpoint and point sources of nutrients were estimated. Nonpoint and point sources combined, including atmospheric deposition, provide a potential annual load of 390 million pounds of nitrogen and 79.5 million pounds of phosphorus. The range of percentages of the estimated nonpoint and point sources that were measured in the stream was 20 to 47 percent for nitrogen and 6 to 14 percent for phosphorus. On the average, the Susquehanna River discharges 141,000 pounds of nitrogen and 7,920 pounds of phosphorus to the Lower Susquehanna River reservoir system each year. About 98 percent of the nitrogen and 60 percent of the phosphorus passes through the reservoir system. Interpretations of available water-quality data and conclusions about the water quality of the Lower Susquehanna River Basin were limited by the scarcity of certain types of water-quality data and current ancillary data. A more complete assessment of the water quality of the basin with respect to nutrients and suspended sediment would be enhanced by the availability of additional data for multiple samples over time from all water environments; samples from streams in the northern and western part of the basin; samples from streams and springs throughout the basin during high base-flow or stormflow conditions; and information on current land-use, and nutrient loading from all types of land-use settings.

Maryland, Pennsylvania

Refining sources of polychlorinated biphenyls in the Back River watershed, Baltimore, Maryland, 2018–2020

Older urban landscapes present unique and complex stressors to urban streams and their habitats through the introduction of legacy and emerging toxic contaminants. Contaminant sources are often associated with various developed land uses such as older residential areas, active and former industrial sites, contaminated sites, and effluents from municipal wastewater treatment plant discharges. These landscapes have a history of legacy contaminant use such as polychlorinated biphenyls (PCBs) resulting in impacts to sediment and water in these complex environments. Despite the ban of PCBs in new commercial use in 1979, PCB contamination is still widespread in the environment, with many fish consumption advisories throughout the Chesapeake Bay region based on elevated PCBs. Several watersheds in the Baltimore region have mandated reductions in PCBs per total maximum daily loads in tidal waters of the watersheds in order to promote compliance with water quality standards. Some of these mandated reductions (for example, regulated watershed runoff) specified in the total maximum daily loads are the responsibility of the local jurisdictions as part of their phase 1 National Pollutant Discharge Elimination System municipal separate storm sewer system permit. In cooperation with the Baltimore City Department of Public Works and Maryland Department of the Environment, the U.S. Geological Survey and University of Maryland, Baltimore County conducted a study from 2018 to 2020 to refine the sources of PCBs from the City of Baltimore into Back River and to use the results to improve the conceptual site model of PCBs in the Back River watershed. PCB concentrations in the water column of the nontidal streams in Back River watershed are relatively consistent throughout both tributaries, with greater concentrations detected in samples collected from Moores Run but greater loads estimated in samples collected from Herring Run. PCB concentrations measured in the bed sediments and analysis of the flux between sediment porewater (hereafter porewater) and surface water within the tributaries suggest that there are no stationary legacy sources within the stream channels. The bulk of PCB mass entering the system from these nontidal tributaries appears to be introduced primarily during storm events. While only one storm event was sampled and concentrations were quantified only in Herring Run, solids captured during the storm were characterized by increases in PCB mass and overall suspended solids concentrations. Although the bioavailability of the PCB-associated sediment is unknown, this mechanism appears to warrant additional attention to better understand how concentrations vary under different storm conditions and temporally. The importance of contaminated stormwater in loading to Herring Run is further supported by the PCB concentrations in storm drain sediments collected near the tributary, which were present in higher concentrations and were characterized by different homolog signatures compared to that in bed sediments. The observations in the tributaries differed from PCB concentrations and sediment characteristics downstream from the City of Baltimore boundary, in the upper tidal area of the main stem of Back River, particularly at the passive sampler locations BRT–1 and BRT–3. This depositional environment is characterized by higher organic content in sediments and higher concentrations of PCBs in porewater, which result in the possible flux of contaminants from sediment to the water column. This flux is generally opposite of that observed in the nontidal tributaries and the farthest upstream tidal site (BRT–2) and may be a result of the possible settling of sediment particles introduced via suspended solids in stormwater. Despite an observed considerable reduction in overall PCB mass loading to and from the Back River Wastewater Treatment Plant (BRWWTP) (and similar reductions observed in biosolids) compared to the estimates previously reported from 2015, effluent from the BRWWTP continues to be a primary source of PCBs to Back River. The current study confirmed the likeliness of fat, oil, and grease deposits within the miles of sewer pipe as a source of PCBs to the BRWWTP influent. The differences between PCB concentrations in fat, oil, and grease deposits found in pipes (during replacement) compared to that of the BRWWTP suggest that legacy deposits may contain higher PCB concentrations and may act as a source of PCBs to passing sewage, eventually entering the BRWWTP. Variation in freely dissolved concentrations in the sewer system was apparent through the analysis of PCBs in the primary pump stations using passive sampling, with the largest contribution to the influent attributed to a single pump station and associated piping. The contribution of PCBs to Herring Run and Moores Run via sanitary sewer overflows compared to the BRWWTP effluent is negligible, similar to reports from another large urban wastewater treatment plant. Therefore, decreased occurrence of sanitary sewer overflows is not expected to largely decrease PCB loads. Results of this study suggest that targeted, sediment-capture best management practices in Back River watershed could be an effective way to reduce PCB mass loading assuming that deposited contaminated sediments are effectively isolated. Recent studies of some common urban best management practices such as bioretention have shown removal of PCBs within the stormwater control structures. In addition, appropriately timed street sweeping practices with appropriate collection equipment may be an effective way to reduce contaminants such as PCBs from road runoff sources. Reductions in concentrations and mass loading within the sewer system measured in this study compared to that estimated 5 years prior reflect the possible success of ongoing gray infrastructure management actions. Reductions may be attributable to enhanced nutrient reduction upgrades to the BRWWTP and extensive capital improvements and maintenance to the sewer system. This study employed a combined sampling approach and a variety of sampling methods to include low-density polyethylene passive samplers, high-volume water samples, and grab samples of both water and sediment to characterize the PCB inputs to Herring Run, Moores Run, and Back River. Incorporating the passive samplers provided a time-weighted average of the freely dissolved concentration in the surface water, porewater, WWTP influent and effluent, and pump station influent over the deployment period with picogram per liter detection limits. A similar monitoring approach from this study could be implemented within other subwatersheds or municipal separate storm sewer system jurisdictions to assist in refining primary sources of PCBs in order to inform appropriate mitigation approaches.

Maryland

Updates to the regional groundwater-flow model of the New Jersey Coastal Plain, 1980–2013

A 21-layer three-dimensional transient groundwater-flow model of the New Jersey Coastal Plain was developed and calibrated by the U.S. Geological Survey (USGS) in cooperation with the New Jersey Department of Environmental Protection to simulate groundwater-flow conditions during 1980–2013, incorporating average annual groundwater withdrawals and average annual groundwater recharge. This model is the third version of the New Jersey Coastal Plain regional groundwater-flow model that was initially developed as part of the USGS Regional Aquifer System Analysis (RASA) program. The model simulates groundwater flow in 11 aquifers and 10 intervening confining units of the New Jersey Coastal Plain to provide a regional overview of groundwater conditions. Averaged groundwater withdrawal data for 1980 to 2013 were used in the model. The 11 aquifers in New Jersey are, from shallowest to deepest, the Holly Beach water-bearing zone and the confined Cohansey aquifer in Cape May County; the Rio Grande water-bearing zone; the Atlantic City 800-foot sand; the Piney Point, Vincentown, and Wenonah-Mount Laurel aquifers; the Englishtown aquifer system; and the upper, middle, and lower aquifers of the Potomac-Raritan-Magothy (PRM) aquifer system. The model was developed with the MODFLOW–2005 numerical code and the UCODE parameter estimation technique and calibrated using water-level and base-flow observations. A total of 3,453 water-level observations from 392 wells in New Jersey and 48 wells in Delaware from 1983 to 2013 were used in model calibration, which includes historical water-level trends for 29 wells in New Jersey during 1980–2013 presented in time-series hydrographs. In addition, derived observations also were included by calculating the vertical gradient at 33 pairs of nested observation wells in New Jersey, for a total of 210 observations. Changes in water levels over time were calculated for 134 wells in New Jersey and four wells in Delaware where water levels had varied substantially (approximately 10 ft) over the 30-year span of synoptic water-level measurements, for a total of 767 observations. A total of 1,485 base-flow observations in 47 surface-water basins in New Jersey from 1980 to 2013 were used in model calibration. Updates to the groundwater-flow model include the conversion to a fully three-dimensional model from the previous quasi-three-dimensional model. The new model will allow for potential future uses such as particle tracking or simulation of variable-density groundwater flow that could not be accomplished with earlier versions of the model. Spatially and temporally variable recharge estimated by using a soil-water balance model resulted in a spatially and temporally finer discretization. The Rio Grande water-bearing zone was added to the model as an aquifer layer to refine estimates of simulated flow in Atlantic and Cape May Counties, New Jersey. Hydrogeologic parameters were updated to include the confining units in New Jersey and corresponding hydrogeologic units in Delaware and eastern Maryland. The simulated water levels for the New Jersey Coastal Plain aquifers were compared to water-level measurements made during 1980–2013. The average residual for 4,243 water-level observations for New Jersey (simulated water levels minus measured water levels) is 1.5 feet. The simulated water-level contours for the confined aquifers for 2013 were compared to potentiometric surfaces produced from water levels measured during 2013. Simulated water levels generally matched the 2013 potentiometric surfaces of the confined aquifers in the areas of large withdrawals. Hydrographs of wells in the confined Coastal Plain aquifers of New Jersey show that simulated water levels generally match the magnitude and seasonal variation of the observed water levels. Hydrographs of base flow for the 47 streamgaging stations in New Jersey indicate that most of the simulated and estimated data match reasonably well. Groundwater withdrawals are an important resource for water supply, agricultural, industrial, and commercial needs in the New Jersey Coastal Plain. Groundwater withdrawals from the New Jersey Coastal Plain aquifers have resulted in persistent, regionally extensive cones of depression in the Englishtown aquifer system and Wenonah-Mount Laurel aquifer in Ocean and Monmouth Counties; Wenonah-Mount Laurel and upper, middle, and lower PRM aquifers in Camden County; and Atlantic City 800-foot sand in Atlantic County. Because hydrologic stresses and water-management needs change with time, periodic updates to the groundwater-flow model are required to provide current information about hydrologic conditions in the New Jersey Coastal Plain and to maintain its usefulness as a tool to manage water resources and develop water-resource strategies. The current updates will support the continued application of this model as a tool for evaluating the regional effects of changes in groundwater withdrawals and of current and potential future water-management strategies on groundwater levels in the New Jersey Coastal Plain.

New Jersey Coastal Plain

Simulation of streamflow and water quality in the Christina River subbasin and overview of simulations in other subbasins of the Christina River Basin, Pennsylvania, Maryland, and Delaware, 1994-98

The Christina River Basin drains 565 square miles (mi 2 ) in Pennsylvania and Delaware and includes the major subbasins of Brandywine Creek, Red Clay Creek, White Clay Creek, and Christina River. The Christina River subbasin (exclusive of the Brandywine, Red Clay, and White Clay Creek subbasins) drains an area of 76 mi 2 . Streams in the Christina River Basin are used for recreation, drinking water supply, and support of aquatic life. Water quality in some parts of the Christina River Basin is impaired and does not support designated uses of the stream. A multi-agency water-quality management strategy included a modeling component to evaluate the effects of point- and nonpoint-source contributions of nutrients and suspended sediment on stream water quality. To assist in nonpoint-source evaluation, four independent models, one for each of the four main subbasins of the Christina River Basin, were developed and calibrated using the model code Hydrological Simulation Program–Fortran (HSPF). Water-quality data for model calibration were collected in each of the four main subbasins and in small subbasins predominantly covered by one land use following a nonpoint- source monitoring plan. Under this plan, stormflow and base-flow samples were collected during 1998 at two sites in the Christina River subbasin and nine sites elsewhere in the Christina River Basin. The HSPF model for the Christina River subbasin simulates streamflow, suspended sediment, and the nutrients, nitrogen and phosphorus. In addition, the model simulates water temperature, dissolved oxygen, biochemical oxygen demand, and plankton as secondary objectives needed to support the sediment and nutrient simulations. For the model, the basin was subdivided into nine reaches draining areas that ranged from 3.8 to 21.9 mi 2 . Ten different pervious land uses and two impervious land uses were selected for simulation. Land-use areas were determined from 1995 land-use data. The predominant land uses in the Christina River subbasin are residential, urban, forested, agricultural, and open. The hydrologic component of the model was run at an hourly time step and calibrated using streamflow data from two U.S. Geological Survey (USGS) streamflow-measurement stations for the period of October 1, 1994, through October 29, 1998. Daily precipitation data from one National Oceanic and Atmospheric Administration (NOAA) meteorologic station and hourly data from one NOAA meteorologic station were used for model input. The difference between observed and simulated streamflow volume ranged from -2.3 to 5.3 percent for a 10-month portion of the calibration period at the two calibration sites. Annual differences between observed and simulated streamflow generally were greater than the overall error for the 4-year period. For example, at Christina River at Coochs Bridge, near the bottom of the free-flowing part of the subbasin (drainage area of 21 mi 2 ), annual differences between observed and simulated streamflow ranged from -6.9 to 6.5 percent and the overall error for the 4-year period was -1.1 percent. Calibration errors for 36 storm periods at the three calibration sites for total volume, low-flow recession rate, 50-percent lowest flows, 10-percent highest flows, and storm peaks were within the recommended criteria of 20 percent or less. Much of the error in simulating storm events on an hourly time step can be attributed to uncertainty in the rainfall data. The water-quality component of the model was calibrated using nonpoint-source monitoring data collected at two USGS streamflow-measurement stations and other water-quality monitoring data. The period of record for water-quality monitoring was variable at the stations, with a start date ranging from October 1994 to January 1998 and an end date of October 1998. Because of availability, monitoring data for suspended-solids concentrations were used as surrogates for suspended-sediment concentrations, although suspended-solids data may underestimate suspended sediment and affect apparent accuracy of the suspended-sediment simulaion. Comparison of observed to simulated loads for up to six storms in 1998 at the two nonpoint-source monitoring sites (Little Mill Creek near Newport and Christina River at Coochs Bridge, Del.) indicate that simulation error is commonly as large as an order of magnitude for suspended sediment and nutrients. The simulation error tends to be smaller for dissolved nutrients than for particulate nutrients. Errors of 40 percent or less for monthly or annual values indicate a fair to good water-quality calibration according to recommended criteria; much larger errors are possible for individual events. Assessment of the water-quality calibration under stormflow conditions is limited by the relatively small amount of available water-quality data in the subbasin. Users of the Christina River subbasin HSPF model and HSPF models for other subbasins in the Christina River Basin should be aware of model limitations and consider the following if the model is used for predictive purposes: streamflow-duration curves suggest the model simulates streamflow reasonably well when measured over a broad range of conditions and time although streamflow and the corresponding water quality for individual storm events may not be well simulated; streamflow-duration curves for the simulation period compare well with duration curves for the 8-year period ending in 2001 at Christina River at Coochs Bridge, Del., and include all but the extreme high-flow and low-flow events; and calibration for water quality was based on limited data, with the result of increasing uncertainty in the water-quality simulation.

Water-Resources Investigations Report

Simulation of streamflow and water quality in the White Clay Creek subbasin of the Christina River Basin, Pennsylvania and Delaware, 1994-98

The Christina River Basin drains 565 square miles (mi 2 ) in Pennsylvania, Maryland, and Delaware. Water from the basin is used for recreation, drinking water supply, and to support aquatic life. The Christina River Basin includes the major subbasins of Brandywine Creek, White Clay Creek, and Red Clay Creek. The White Clay Creek is the second largest of the subbasins and drains an area of 108 mi 2 . Water quality in some parts of the Christina River Basin is impaired and does not support designated uses of the streams. A multi-agency water-quality management strategy included a modeling component to evaluate the effects of point and nonpoint-source contributions of nutrients and suspended sediment on stream water quality. To assist in non point-source evaluation, four independent models, one for each of the three major subbasins and for the Christina River, were developed and calibrated using the model code Hydrological Simulation Program—Fortran (HSPF). Water-quality data for model calibration were collected in each of the four main subbasins and in smaller subbasins predominantly covered by one land use following a nonpoint-source monitoring plan. Under this plan, stormflow and base- flow samples were collected during 1998 at two sites in the White Clay Creek subbasin and at nine sites in the other subbasins. The HSPF model for the White Clay Creek Basin simulates streamflow, suspended sediment, and the nutrients, nitrogen and phosphorus. In addition, the model simulates water temperature, dissolved oxygen, biochemical oxygen demand, and plankton as secondary objectives needed to support the sediment and nutrient simulations. For the model, the basin was subdivided into 17 reaches draining areas that ranged from 1.37 to 13 mi 2 . Ten different pervious land uses and two impervious land uses were selected for simulation. Land-use areas were determined from 1995 land-use data. The predominant land uses in the White Clay Creek Basin are agricultural, forested, residential, and urban. The hydrologic component of the model was run at an hourly time step and primarily calibrated using streamflow data from two U.S. Geological Survey (USGS) streamflow-measurement stations for the period of October 1, 1994, through October 29, 1998. Additional calibration was done using data from two other USGS streamflow-measurement stations with periods of record shorter than the calibration period. Daily precipitation data from two National Oceanic and Atmospheric Administration (NOAA) gages and hourly precipitation and other meteorological data for one NOAA gage were used for model input. The difference between simulated and observed streamflow volume ranged from -0.9 to 1.8 percent for the 4-year period at the two calibration sites with 4-year records. Annual differences between observed and simulated streamflow generally were greater than the overall error. For example, at a site near the bottom of the basin (drainage area of 89.1 mi 2 ), annual differences between observed and simulated streamflow ranged from -5.8 to 14.4 percent and the overall error for the 4-year period was -0.9 percent. Calibration errors for 36 storm periods at the two calibration sites for total volume, low-flowrecession rate, 50-percent lowest flows, 10-percent highest flows, and storm peaks were within the recommended criteria of 20 percent or less. Much of the error in simulating storm events on an hourly time step can be attributed to uncertainty in the hourly rainfall data. The water-quality component of the model was calibrated using data collected by the USGS and state agencies at three USGS streamflow-measurement stations with variable water-quality monitoring periods ending October 1998. Because of availability, monitoring data for suspended-solids concentrations were used as surrogates for suspended-sediment concentrations, although suspended solids may underestimate suspended sediment and affect apparent accuracy of the suspended-sediment simulation. Comparison of observed to simulated loads for up to five storms in 1998 at each of the two nonpoint-source monitoring sites in the White Clay Creek Basin indicate that simulation error is commonly as large as an order of magnitude for suspended sediment and nutrients. The simulation error tends to be smaller for dissolved nutrients than for particulate nutrients. Errors of 40 percent or less for monthly or annual values indicate a fair to good water-quality calibration according to recommended criteria, with much larger errors possible for individual events. The accuracy of the water-quality calibration under stormflow conditions is limited by the relatively small amount of water-quality data available for the White Clay Creek Basin. Users of the White Clay Creek HSPF model should be aware of model limitations and consider the following if the model is used for predictive purposes: streamflow and water quality for individual storm events may not be well simulated, but the model performance is reasonable when evaluated over longer periods of time; the observed flow-duration curve for the simulation period is similar to the long-term flow-duration curve at White Clay Creek near Newark, Del., indicating that the calibration period is representative of all but highest 0.1 percent and lowest 0.1 percent of flows at that site; relative errors in streamflow and water-quality simulations are greater for smaller drainage areas than for larger areas; and calibration for water-quality was based on sparse data.

Water-Resources Investigations Report

Nutrient concentrations and loads in the northeastern United States - Status and trends, 1975-2003

The U.S. Geological Survey (USGS) National Water-Quality Assessment Program (NAWQA) began regional studies in 2003 to synthesize information on nutrient concentrations, trends, stream loads, and sources. In the northeastern United States, a study area that extends from Maine to central Virginia, nutrient data were evaluated for 130 USGS water-quality monitoring stations. Nutrient data were analyzed for trends in flow-adjusted concentrations, modeled instream (non-flow-adjusted) concentrations, and stream loads for 32 stations with 22 to 29 years of water-quality and daily mean streamflow record during 1975-2003 (termed the long-term period), and for 46 stations during 1993-2003 (termed the recent period), by using a coupled statistical model of streamflow and water quality developed by the USGS. Recent trends in flow-adjusted concentrations of one or more nutrients also were analyzed for 90 stations by using Tobit regression. Annual stream nutrient loads were estimated, and annual nutrient yields were calculated, for 47 stations for the long-term and recent periods, and for 37 additional stations that did not have a complete streamflow and water-quality record for 1993-2003. Nutrient yield information was incorporated for 9 drainage basins evaluated in a national NAWQA study, for a total of 93 stations evaluated for nutrient yields. Long-term downward trends in flow-adjusted concentrations of total nitrogen and total phosphorus (18 and 19 of 32 stations, respectively) indicate regional improvements in nutrient-related water-quality conditions. Most of the recent trends detected for total phosphorus were upward (17 of 83 stations), indicating possible reversals to the long-term improvements. Concentrations of nutrients in many streams persist at levels that are likely to affect aquatic habitat adversely and promote freshwater or coastal eutrophication. Recent trends for modeled instream concentrations, and modeled reference concentrations, were evaluated relative to ecoregion-based nutrient criteria proposed by the U.S. Environmental Protection Agency. Instream concentrations of total nitrogen and total phosphorus persist at levels higher than proposed criteria at more than one-third and about one-half, respectively, of the 46 stations analyzed. Long-term trends in nutrient loads were primarily downward, with downward trends in total nitrogen and total phosphorus loads detected at 12 and 17 of 32 stations, respectively. Upward trends were rare, with one upward trend for total nitrogen loads and none for total phosphorus. Trends in loads of nitrite-plus-nitrate nitrogen included 7 upward and 8 downward trends among 32 stations. Downward trends in loads of ammonia nitrogen and total Kjeldahl nitrogen were detected at all six stations evaluated. Long-term downward trends detected in four of the five largest drainage basins evaluated include: total nitrogen loads for the Connecticut, Delaware, and James Rivers; total Kjeldahl nitrogen and ammonia nitrogen loads for the Susquehanna River; ammonia nitrogen and nitrite-plus-nitrate nitrogen loads for the James River; and total phosphorus loads for the Connecticut and Delaware Rivers. No trends in load were detected for the Potomac River. Nutrient yields were evaluated relative to the extent of land development in 93 drainage basins. The undeveloped land-use category included forested drainage basins with undeveloped land ranging from 75 to 100 percent of basin area. Median total nitrogen yields for the 27 undeveloped drainage basins evaluated, including 9 basins evaluated in a national NAWQA study, ranged from 290 to 4,800 pounds per square mile per year (lb/mi 2 /yr). Total nitrogen yields even in the most pristine drainage basins may be elevated relative to natural conditions, because of high rates of atmospheric deposition of nitrogen in parts of the northeastern United States. Median total phosphorus yields ranged from 12 to 330 lb/mi 2 /yr for the 26 undeveloped basins evaluated. The undeveloped category includes some large drainage basins with point-source discharges and small percentages of developed land; in these basins, streamflow from undeveloped headwater areas dilutes streamflow in more urbanized reaches, and dampens but does not eliminate the point-source "signal" of higher nutrient loads. Median total nitrogen yields generally do not exceed 1,700 lb/mi 2 /yr, and median total phosphorus yields generally do not exceed 100 lb/mi 2 /yr, in the drainage basins that are least affected by human land-use and waste-disposal practices. Agricultural and urban land use has increased nutrient yields substantially relative to undeveloped drainage basins. Median total nitrogen yields for 24 agricultural basins ranged from 1,700 to 26,000 lb/mi 2 /yr, and median total phosphorus yields ranged from 94 to 1,000 lb/mi 2 /yr. The maximum estimated total nitrogen and total phosphorus yields, 32,000 and 16,000 lb/mi 2 /yr, respectively, for all stations in the region were in small (less than 50 square miles (mi 2 )) agricultural drainage basins. Median total nitrogen yields ranged from 1,400 to 17,000 lb/mi 2 /yr in 26 urbanized drainage basins, and median total phosphorus yields ranged from 43 to 1,900 lb/mi 2 /yr. Urbanized drainage basins with the highest nutrient yields are generally small (less than 300 mi2) and are drained by streams that receive major point-source discharges. Instream nutrient loads were evaluated relative to loads from point-source discharges in four drainage basins: the Quinebaug River Basin in Connecticut, Massachusetts, and Rhode Island; the Raritan River Basin in New Jersey; the Patuxent River Basin in Maryland; and the James River Basin in Virginia. Long-term downward trends in nutrient loads, coupled with similar trends in flow-adjusted nutrient concentrations, indicate long-term reductions in the delivery of most nutrients to these streams. However, the absence of recent downward trends in load for most nutrients, coupled with instream concentrations that exceed proposed nutrient criteria in several of these waste-receiving streams, indicates that challenges remain in reducing delivery of nutrients to streams from point sources. During dry years, the total nutrient load from point sources in some of the drainage basins approached or equaled the nutrient load transported by the stream.

Connecticut;Delaware;Maine;Maryl;Massachusetts;New

Pesticides in groundwater in the Anacostia River and Rock Creek watersheds in Washington, D.C., 2005 and 2008

The U.S. Geological Survey (USGS), in cooperation with the District Department of the Environment, conducted a groundwater-quality investigation to (a) determine the presence, concentrations, and distribution of selected pesticides in groundwater, and (b) assess the presence of pesticides in groundwater in relation to selected landscape, hydrogeologic, and groundwater-quality characteristics in the shallow groundwater underlying the Anacostia River and Rock Creek watersheds in Washington, D.C. With one exception, well depths were 100 feet or less below land surface. The USGS obtained or compiled ancillary data and information on land use (2001), subsurface sediments, and groundwater samples from 17 wells in the lower Anacostia River watershed from September through December 2005, and from 14 wells in the lower Anacostia River and lower Rock Creek watersheds from August through September 2008. Twenty-seven pesticide compounds, reflecting at least 19 different types of pesticides, were detected in the groundwater samples obtained in 2005 and 2008. No fungicides were detected. In relation to the pesticides detected, degradate compounds were as or more likely to be detected than applied (parent) compounds. The detected pesticides chiefly reflected herbicides commonly used in urban settings for non-specific weed control or insecticides used for nonspecific haustellate insects (insects with specialized mouthparts for sucking liquid) or termite-specific control. Detected pesticides included a combination of pesticides currently (2008) in use, banned or under highly restricted use, and some that had replaced the banned or restricted-use pesticides. The presence of banned and restricted-use pesticides illustrates their continued persistence and resistance to complete degradation in the environment. The presence of the replacement pesticides indicates the susceptibility of the surficial aquifer to contamination irrespective of the changes in the pesticides used. A preliminary review of the data collected in 2005 and 2008 indicated that differences in the surficial geology, land use (as a surrogate for pesticide use), and above-average precipitation for most of 2004 through 2008, as well as differences in the number and performance of USGS laboratory methods used, could have led to more pesticides detected in groundwater samples collected in 2008 than in groundwater samples collected in 2005. Thus, although data from both years of collection were used for interpretive analysis, emphasis was placed on the analysis of the data obtained in 2008. The presence of pesticides in shallow groundwater (less than approximately 100 ft (feet), or 30 m (meters), below land surface) indicated at least the upper surficial aquifer in Washington, D.C. was susceptible to contamination. One or more herbicides or insecticides were detected in groundwater samples collected from 50 percent of the shallow wells sampled in 2005, and from 62 percent of the shallow wells sampled in 2008. Differences among types of pesticides in shallow groundwater were apparent. The most frequently detected class of herbicides was the s-triazine compounds—atrazine, simazine, or prometon, or the atrazine-degradate compounds—2-chloro-4-ethylamino-6-amino-s-triazine (desethylatrazine or CIAT) and 2-chloro-4-isopropylamino-6-amino-s-triazine (hydroxyatrazine or OIET). The next most frequently detected classes of herbicides were the chloroacetanilides, including metolachlor and acetochlor, and the ureic herbicides, including diuron (and degradate, 3,4-dichloroaniline), fluometuron, metsulfuron methyl, sulfameturon, bromacil, and tebuthiuron. Insecticides also were detected, but less frequently than herbicides, with one or more insecticides present in groundwater samples from 38 percent of shallow wells sampled in 2008. Detected insecticides included parent or degradate compounds commonly used for either nonspecific or haustellate (sucking) insects, including chlorpyrifos and dichlorodiphenyldichloroethane (p,p’-DDD; a degradate of dichlorodiphenyltrichloroethane, DDT), and for termite control, including dieldrin, chlordane, heptachlor epoxide, (a degradate of heptachlor), fipronil, and the sulfone and sulfide degradates of fipronil. The concentrations of individual pesticides in shallow groundwater in both years were low. Maximum concentrations were no greater than a few tenths of a microgram per liter (μg/L); typical concentrations often were less than 0.1 μg/L. Multiple pesticides, however, commonly were present in groundwater. For example, in 2008, approximately 88 percent (7 of 8) of the wells that yielded a sample with at least one detectable pesticide contained five or more pesticides. The highest number of detections occurred in a groundwater sample from well WE Ca 32, which is located in a highly developed urban area; this sample contained 15 different pesticide residues. In relation to human and aquatic health, no pesticide concentration in either 2005 or 2008 exceeded Federal drinking-water standards. Groundwater samples from a few sites, however, contained levels of chiefly banned or restricted-use pesticides that exceeded other human-health and (or) aquatic-health guidelines. For example, concentrations of dieldrin in 2008 groundwater samples from three wells—WE Ca 32 (0.028 μg/L), WE Ba 11 (0.016 μg/L), and WW Ac 8 (0.014 μg/L)—fell within the range of concern for 2004 Federally approved non-regulatory USGS Health-Based Assessment benchmarks (0.002 to 0.2 μg/L), and exceeded earlier (1999) Federal criteria for drinking water (0.000052 μg/L). Other individual compounds whose concentrations exceeded 1999 Federal guidelines for samples from one or more of these three sites, or another site, included p,p’-DDD, dichlorodiphenyldichloroethylene (p,p’-DDE; another degradate of DDT), chlordane, and heptachlor epoxide. Pesticide concentrations in groundwater also were compared to three aquatic-health guidelines for freshwater (United States, Great Lakes, or Canada). One or more of these guidelines were exceeded in groundwater samples obtained in 2005 or 2008 for one or more of the compounds chlordane, dieldrin, heptachlor epoxide, p,p’-DDE, p,p’-DDD, and chlorpyrifos. The spatial distribution of pesticides in the shallow groundwater appeared to be related, in part, to land use, a surrogate for pesticide use. Although most of the wells sampled in this study are in parklands or other relatively open and accessible space, multiple pesticides most often were detected in 2008 groundwater samples collected from wells where a considerable percentage (in excess of 60 percent) of the land within a 500-m radius is developed space (residential, commercial, or other urban infrastructure). Insecticides were detected in wells surrounded by at least 50 percent, and most commonly by more than 80 percent, development. Well WE Ca 32, the site associated with the highest number of pesticide residues in groundwater (8 herbicides and 7 insecticides), is in a small residential park, where 99 percent of the surrounding land is well-maintained residential and commercial development. The vertical distribution of detected pesticides in shallow groundwater appeared to be related, in part, to depth below land surface, surficial-bedrock type, and differences in the chemistry of shallow groundwater. Pesticides were detected at relatively shallow depths in wells that may not have fully penetrated the shallow aquifer. For wells in which at least one pesticide was detected, the median depth below land surface to the top of the well screen was 5.8 m, and the maximum depth was 8.5 m. Among the types of surficial materials in which wells were completed—alluvium, terrace deposits, or Potomac Formation sub- or outcrops in the Coastal Plain Province, and saprolite or fractured bedrock (Laurel and Sykesville Formations) underlying saprolite in the Piedmont Province—no pesticides were detected in groundwater associated with wells completed in the alluvium or fractured bedrock. Detections occurred in some but not all wells completed in the other surficial materials. Overall, the pattern in occurrence appeared related to the local permeability of these sediments and groundwater chemistry. Groundwater with multiple pesticide detections tended to occur in permeable sediments (absent any appreciable overlying clay, silt, or clay-silt layers), in conjunction with other common urban contaminants (elevated chloride in excess of tens to hundreds of milligrams per liter (mg/L), and oxic, rather than reduced, groundwater as evidenced by elevated (in excess of 5 mg/L) concentrations of nitrate). The results of this investigation were compared to results from two other similar and recent studies on pesticide occurrence in the shallow aquifer. These included a study in the nearby Maryland and Delaware Coastal Plain Physiographic Province and one in the Maryland and Virginia Piedmont Physiographic Province. Results from these studies were similar to the current study in relation to (a) the types, frequencies, concentrations, and mixtures of pesticides detected; (b) compounds that exceeded human-and aquatic-health criteria; and (c) the occurrence and distribution of pesticides within the surficial aquifer in relation to depths, sediment types, and groundwater chemistries.

Washington D.C.

Assessment of groundwater availability in the Northern Atlantic Coastal Plain aquifer system From Long Island, New York, to North Carolina

Executive Summary The U.S. Geological Survey began a multiyear regional assessment of groundwater availability in the Northern Atlantic Coastal Plain (NACP) aquifer system in 2010 as part of its ongoing regional assessments of groundwater availability of the principal aquifers of the Nation. The goals of this national assessment are to document effects of human activities on water levels and groundwater storage, explore climate variability effects on the regional water budget, and provide consistent and integrated information that is useful to those who use and manage the groundwater resource. As part of this nationwide assessment, the USGS evaluated available groundwater resources within the NACP aquifer system from Long Island, New York, to northeastern North Carolina. The northern Atlantic Coastal Plain physiographic province depends heavily on groundwater to meet agricultural, industrial, and municipal needs. The groundwater assessment of the NACP aquifer system included an evaluation of how water use has changed over time; this evaluation primarily used groundwater budgets and development of a numerical modeling tool to assess system responses to stresses from future human uses and climate trends. This assessment focused on multiple spatial and temporal scales to examine changes in groundwater pumping, storage, and water levels. The regional scale provides a broad view of the sources and demands on the system with time. The sub-regional scale provides an evaluation of the differing response of the aquifer system across geographic areas allowing for closer examination of the interaction between different aquifers and confining units and the changes in these interactions under pumping and recharge conditions in 2013 and hydrologic stresses as much as 45 years in the future. By focusing on multiple scales, water-resource managers may utilize this study to understand system response to changes as they affect the system as a whole. The NACP aquifer system extends from Long Island to northeastern North Carolina, and includes aquifers primarily within New York, New Jersey, Delaware, Maryland, Virginia, and North Carolina. The seaward-dipping sedimentary wedge that underlies the northern Atlantic Coastal Plain physiographic province forms a complex groundwater system. Although the NACP aquifer system is recognized by the U.S. Geological Survey as one of the smallest of the 66 principal aquifer systems in the Nation, it ranks 13th overall in terms of total groundwater withdrawals and is 7th in population served. Despite abundant precipitation [about 45 inches per year (in/yr)], the supply of fresh surface water in this region is limited because many of the surface waters in this area are brackish estuaries, contributing to why many communities in the northern Atlantic Coastal Plain physiographic province rely heavily on groundwater to meet their water needs. Increases in population and changes in land use during the past 100 years have resulted in diverse increased demands for freshwater throughout the northern Atlantic Coastal Plain physiographic province with groundwater serving as a vital source of drinking water for the nearly 20 million people who live in the region. Total groundwater withdrawal in 2013 was estimated to be about 1,300 million gallons per day (Mgal/d) and accounts for about 40 percent of the drinking water supply with the densely populated areas tending to have the highest rates of withdrawals and, therefore, being most susceptible to effects from these withdrawals over time. Water levels in many of the confined aquifers are decreasing by as much as 2 feet per year (ft/yr) in response to extensive development and subsequent increased withdrawals throughout the region. Total water-level decreases (drawdowns) are more than 100 feet (ft) in some aquifers from their predevelopment (before 1900) levels. These drawdowns extend across state lines and under the Chesapeake and Delaware Bays, creating the potential for interstate aquifer management issues. Regional water-resources managers in the northern Atlantic Coastal Plain physiographic province face challenges beyond competing local domestic, industrial, agricultural, and environmental demands for water. Large changes in regional water use have made the State-level management of aquifer resources more difficult because of hydrologic effects that extend beyond State boundaries. The northern Atlantic Coastal Plain physiographic province is underlain by a wedge of unconsolidated to partially consolidated sediments that are typically thousands of feet thick along the coastline with a maximum thickness of about 10,000 ft near the edge of the continental shelf. The NACP aquifer system consists of nine confined aquifers and nine confining units capped by an unconfined surficial aquifer that is bounded laterally from the west by the contact between Coastal Plain sediments and the upland Piedmont bedrock. This aquifer system extends to the east to the limit of the Continental Shelf, however, the boundary between fresh and saline groundwater is considered to be much closer to the shoreline and varies vertically by aquifer. Precipitation over the region for average conditions from 2005 to 2009 is about 61,800 Mgal/d, but about 70 percent of it is lost to evapotranspiration resulting in an inflow of about 19,600 Mgal/d entering the groundwater system as aquifer recharge. Most of this recharge enters the aquifer system and flows through the shallow unconfined aquifer and either discharges to streams or directly to coastal waters without reaching the deep, confined aquifer system. In addition to recharge from precipitation, other sources of water include the return of wastewater from domestic septic systems of about 240 Mgal/d, about 60 Mgal/d of water released from storage in the confined system, and about 30 Mgal/d of lateral inflow at the boundary between freshwater and saltwater in response to pumping for conditions in 2013. The outflow needed to balance the inflows was subdivided between streamflow, discharge to tidal portions of streams, and coastal discharge. The hydrologic budget developed for current [2013] conditions determined that 93 percent of the total outflow was to surface waters with about 70 percent divided evenly between streamflow and shallow coastal discharge and 23 percent as discharge to tidal waters. The remaining 7 percent of the total outflow components include withdrawals from both the surficial and confined aquifers of the groundwater system. The groundwater availability assessment of the NACP aquifer system highlights the importance of analyses at both the regional and local scales to understand how changes in land use, water use, and climate have affected groundwater resources and how these resources may change in the future. The investigation included assessments of the regional changes in water levels and budgets across State lines, the importance of considering storage change in the confining units, the response of the aquifer system to a continuation of current [2013] hydrologic stresses into the future, and the potential effects of climate change and sea-level rise on the aquifer system. The Potomac aquifer group includes two of the most widely used aquifers in the NACP aquifer system, the Potomac-Patapsco and Potomac-Patuxent regional aquifers, providing about 24 percent of the total groundwater used in the region. Withdrawals from large pumping centers in this deep, confined aquifer group have resulted in substantial decreases in water-levels across state lines, particularly between southern Virginia and northeastern North Carolina as well as between southern New Jersey and northern Delaware where water levels in the Potomac-Patapsco aquifer have decreased by as much as 200 ft and 50 ft, respectively from predevelopment to current [2013] conditions. This response in water levels also is reflected in changes in water budgets where, for example, about 20 percent of the total response to pumping in Virginia is met by inducing flow from adjacent States. Understanding and quantifying these hydrologic effects that extend beyond State boundaries is critical for the State- and regional-level management of aquifer resources. The cumulative storage loss from the intervening confining units throughout the entire NACP aquifer system was about 35 percent of the total storage loss from predevelopment to current [2013] conditions. In geographic areas such as Delmarva Peninsula, Maryland, and New Jersey, the water released from storage in the confining units makes up the majority of the total storage release from the groundwater system and is becoming proportionally more important over time as the surficial aquifer approaches equilibrium with respect to pumping and recharge stresses as of 2013. Storage loss from the confining units is of particular concern because, unlike in the sands that comprise the confined aquifers, water removed from the clayey confining unit sediments cannot be replenished as these units gradually compress. This non-recoverable storage loss, if great enough, can result in land subsidence where these units are thick and the release from storage is relatively large and contributes to increased concerns for sea-level rise in areas such as the lower portion of the Chesapeake Bay. Groundwater usage increased dramatically in the NACP aquifer system during post-World War II era from the mid-1940s to early the 1980s, with withdrawals increasing from about 400 Mgal/d to more than 1,300 Mgal/d. Although groundwater withdrawals have been relatively constant since the early 1980s, about half of the total groundwater withdrawn from the NACP aquifer system since 1900 was withdrawn in the past 30 years. An analysis of the response of the groundwater system to a continuation of the current [2013] pumping for an additional 30 years into the future shows that the flow system continues to adjust in terms of changes in water budget components, water levels, and the boundary between freshwater and saltwater as it approaches equilibrium. The largest change in water budget components is the reduction in the amount of water released from storage. Across the entire NACP aquifer system, the reduction of storage release from 7 to 4 percent of the total water budget change is accounted for by reductions in groundwater discharge to streams and coastal waters. Locally, a similar response is calculated for each of the geographic areas except for Virginia where the amount of water released from storage accounts for about 25 percent of the total change in water budget. This finding suggests that the groundwater flow system in Virginia is not approaching equilibrium under the current [2013] stresses and, therefore, water levels will continue to decrease even if the pumping remains constant. An analysis of the change in water levels in the Potomac-Patapsco aquifer as pumping is continued 30 years into the future reveals that the largest decreases in water levels throughout the NACP aquifer system will occur in the southern Virginia and northeastern North Carolina parts of the study area. It is these areas that also see the greatest potential for increased lateral movement of saline groundwater in the deep, confined portion of the groundwater flow system in response to a continuation of the current [2013] pumping rates. The potential effects of long-term climate change and variability on the hydrologic system and availability of water resources in the NACP aquifer system continue to be of serious societal concern. These concerns include the effects of changes in aquifer recharge and in sea-level rise on the groundwater flow system. An assessment of the potential effects of a prolonged drought during current [2013] pumping conditions indicated that the reductions in recharge associated with droughts, including additional irrigation withdrawals required to meet increased crop water demand, have the greatest effects on water levels and streamflows in the surficial aquifer, and changes in water levels in the confined aquifers primarily resulted from the increased withdrawals associated with increased irrigation pumping; this response was most apparent in the Delmarva Peninsula. These results suggest that water levels may not be susceptible to the effects of droughts in the confined aquifers of the NACP aquifer system not used for irrigation, unlike in the unconfined surficial aquifer. A second analysis also was conducted to assess the effects of sea-level rise on the groundwater system throughout the northern Atlantic Coastal Plain physiographic province because recent analyses of the relative rates of sea-level rise along the Atlantic coast indicate that the Mid-Atlantic region represents a hot spot with anomalously higher rates of sea-level rise than observed elsewhere in the United States. Groundwater levels rose from 0 to 3 ft in response to a 3-ft simulated change in sea-level position, with the largest response occurring along the shoreline and away from non-tidal streams. About 37 percent (or 10,000 square miles) of the area of the northern Atlantic Coastal Plain physiographic province may experience about a 0.5-ft or more increase in water levels with the 3-ft increase in sea-level position, whereas about 18 percent (almost 5,000 square miles) of land of the northern Atlantic Coastal Plain physiographic province may experience a 2-ft or more increase in water levels with the 3-ft increase in sea-level position. These increases in the water table are of particular concern in low-lying areas where the unsaturated (vadose) zone is already thin, thus creating concerns for groundwater inundation of subsurface infrastructure, such as basements, septic systems, and subway systems. This increase in the water table also will likely alter the distribution of groundwater discharge to surface-water bodies thus increasing groundwater flow to streams that would have otherwise discharged directly to coastal waters. Throughout the NACP aquifer system, this redistribution of groundwater discharge results in an additional 2 percent of base flow in streams. Although the increases in groundwater discharge to streams (and corresponding decreases in discharge to coastal waters) calculated for the entire NACP aquifer system and its geographic areas represent only a small increase compared with current [2013] conditions, this redistribution of groundwater discharge from the coast to streams locally can alter the delivery of freshwater input to coastal receiving waters and have ecohydrological implications on the sensitive ecosystems which rely on a balance of groundwater discharge and surface-water flow.

Delaware, Maryland, New Jersey, New York, North Ca