Geology ReportsSearch

SEARCH · Geology Reports

Results for “Geologic report (New Jersey)”

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.

376 records · Page 21Linked to original sources

Relations of Water Quality to Streamflow, Season, and Land Use for Four Tributaries to the Toms River, Ocean County, New Jersey, 1994-99

The effects of nonpoint-source contamination on the water quality of four tributaries to the Toms River in Ocean County, New Jersey, have been investigated in a 5-year study by the U.S. Geological Survey (USGS), in cooperation with the New Jersey Department of Environmental Protection (NJDEP). The purpose of the study was to relate the extent of land development to loads of nutrients and other contaminants to these streams, and ultimately to Barnegat Bay. Volumetric streamflow (discharge) was measured at 6 monitoring sites during 37 stormflow and base-flow sampling events over a 5-year period (May 1994-September 1999). Concentrations and yields (area-normalized instantaneous load values) of nitrogen and phosphorus species, total suspended solids, and fecal coliform bacteria were quantified, and pH, dissolved oxygen, and stream stage were monitored during base-flow conditions and storms. Sufficient data were collected to allow for a statistical evaluation of differences in water quality among streams in subbasins with high, medium, and low levels of land development. Long Swamp Creek, in a highly developed subbasin (64.2 percent developed); Wrangle Brook, in a moderately developed subbasin (34.5 percent); Davenport Branch, in a slightly developed subbasin (22.8 percent); and Jakes Branch, in an undeveloped subbasin (0 percent) are the subbasins selected for this study. No point-source discharges are known to be present on these streams. Water samples were collected and analyzed by the NJDEP, and discharge measurements and data analysis were conducted by the USGS. Total nitrogen concentrations were lower in Davenport Branch than in Long Swamp Creek and Wrangle Brook during base flow and stormflow. Concentrations of total nitrogen and nitrate were highest in Wrangle Brook (as high as 3.0 mg/L and 1.6 mg/L, respectively) as a result of high concentrations of nitrate in samples collected during base flow; nitrate loading from ground-water discharge is much higher in Wrangle Brook than in any of the other streams, possibly as a result of an experimental wastewater-(secondary effluent) disposal site that was in operation during the 1980's. Ammonia concentrations were higher in samples from Long Swamp Creek than in those from the other two monitoring sites under all flow conditions, and ammonia yields were higher during stormflow than base flow at all monitoring sites. Concentrations and yields of fecal coliform bacteria and total suspended solids were higher during stormflow than during base flow at all monitoring sites. Concentrations and yields were significantly higher in Long Swamp Creek, a highly developed subbasin and Wrangle Brook, a moderately developed subbasin than in Davenport Branch, a slightly developed subbasin. Concentrations and yields of phosphate species, which also are strongly related to stormflow, were higher during stormflow in Long Swamp Creek than in the other subbasins. Base-flow separation techniques were used on hydrographs generated for storms to distinguish the fraction of discharge and constituent loading attributable to storm runoff (overland flow) from the fraction contributed by ground-water discharge. Precipitation records were used to determine the total annual volumes of ground-water discharge and runoff at each monitoring site. These volumes were used in conjunction with water-quality data to calculate total annual loads of each constituent at each monitoring site, separated into ground-water discharge and runoff fractions. It was determined that loads of ammonia, nitrate, organic nitrogen, total nitrogen, and orthophosphate in ground-water discharge were significantly higher in the moderately developed Wrangle Brook subbasin than in the highly developed Long Swamp Creek subbasin, and that no relation was apparent between the percent of land development and constituent loads from ground-water discharge. The loading of each constituent contributed by ground-water discharge is specific

Scientific Investigations Report

Water-level conditions in the confined aquifers of the New Jersey Coastal Plain, 2013

The Coastal Plain aquifers of New Jersey provide an important source of water for more than 3.5 million people. In 2013, groundwater withdrawals from 10 confined aquifers of the New Jersey Coastal Plain totaled about 190 million gallons per day. 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 many confined Coastal Plain aquifers. Starting in 1978, the U.S. Geological Survey (USGS) began mapping the potentiometric surfaces of the major confined Coastal Plain aquifers every 5 years to provide a regional assessment of groundwater conditions. In a study conducted by the USGS, in cooperation with the New Jersey Department of Environmental Protection, water levels in 10 confined aquifers of the New Jersey Coastal Plain were measured and evaluated to provide a regional overview of groundwater conditions during fall 2013. Water levels were measured in 987 wells in New Jersey, and parts of Pennsylvania and Delaware. Potentiometric-surface maps were prepared for, in ascending order of age, the confined Cohansey aquifer of Cape May County, Rio Grande water-bearing zone, Atlantic City 800-foot sand, Piney Point aquifer, Vincentown aquifer, Wenonah-Mount Laurel aquifer, Englishtown aquifer system, and the Upper, Middle, and Lower aquifers of the Potomac-Raritan-Magothy (PRM) aquifer system. Persistent, regionally extensive cones of depression were present in the potentiometric surfaces of 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. Changes in water levels from 2008 to 2013 were measured in many Coastal Plain aquifers in New Jersey. In some areas, water levels continued to decline as a result of pumping, but in other areas water levels continued to recover as a result of regulated decreases in groundwater withdrawals. Since 2008, in the confined Cohansey aquifer in Cape May County, water levels generally did not change; however, cones of depression in the potentiometric surface of the Piney Point aquifer in some areas of Cumberland County deepened by more than 20 feet (ft). In Critical Area 1, an area of restricted withdrawals, measured water levels in the Wenonah-Mount Laurel aquifer declined in parts of southern Monmouth County by more than 10 ft; however, rises in water levels of more than 10 ft were measured in parts of northern Ocean and Monmouth Counties. Since 2008, in Critical Area 2, also an area of restricted withdrawals, measured water levels in the Wenonah-Mount Laurel aquifer rose more than 20 ft in parts of western Burlington County and more than 20 ft in parts of western Camden County. Since 2008, in Critical Area 1, measured water levels in the Englishtown aquifer system declined in parts of eastern Ocean County by more than 10 ft and in southeastern Monmouth County by more than 20 ft; however, rises in water levels of more than 10 ft were measured in other parts of Ocean and Monmouth Counties. In general, since 2008 in Critical Area 2, in the Upper PRM aquifer, measured water levels continued to rise by 10 ft or more in central and western Burlington and central Camden Counties. In the Middle PRM aquifer in Critical Area 2, measured water levels rose in parts of central Camden County by 10 ft or more. However, measured water levels in the Lower PRM aquifer in Critical Area 2 were more than 10 ft lower in the center of the cone of depression in central Camden County, but measured water levels continued to rise updip from this area in Critical Area 2. Seasonal water-level fluctuations are presented in time-series hydrographs for 77 wells during 1978–2013. Analyses of long-term water-level changes for the period 2008–13 indicate downward water-level trends at 14 wells (18 percent), upward trends at 34 wells (44 percent), and no substantial change at 29 wells (38 percent). Downward trends were most often observed for wells screened in the Piney Point aquifer and the Atlantic City 800-foot sand. Upward water-level trends were most often measured for wells screened in the PRM aquifer system. Upward water-level trends also were measured for wells in the Englishtown aquifer system and the Wenonah-Mount Laurel aquifer in Critical Area 1 in some areas; however, downward trends and no substantial changes were measured in other areas.

New Jersey

Results from laboratory and field testing of nitrate measuring spectrophotometers

Five ultraviolet (UV) spectrophotometer nitrate analyzers were evaluated by the U.S. Geological Survey (USGS) Hydrologic Instrumentation Facility (HIF) during a two-phase evaluation. In Phase I, the TriOS ProPs (10-millimeter (mm) path length), Hach NITRATAX plus sc (5-mm path length), Satlantic Submersible UV Nitrate Analyzer (SUNA, 10-mm path length), and S::CAN Spectro::lyser (5-mm path length) were evaluated in the HIF Water-Quality Servicing Laboratory to determine the validity of the manufacturer's technical specifications for accuracy, limit of linearity (LOL), drift, and range of operating temperature. Accuracy specifications were met in the TriOS, Hach, and SUNA. The stock calibration of the S::CAN required two offset adjustments before the analyzer met the manufacturer's accuracy specification. Instrument drift was observed only in the S::CAN and was the result of leaching from the optical path insert seals. All tested models, except for the Hach, met their specified LOL in the laboratory testing. The Hach's range was found to be approximately 18 milligrams nitrogen per liter (mg-N/L) and not the manufacturer-specified 25 mg-N/L. Measurements by all of the tested analyzers showed signs of hysteresis in the operating temperature tests. Only the SUNA measurements demonstrated excessive noise and instability in temperatures above 20 degrees Celsius (°C). The SUNA analyzer was returned to the manufacturer at the completion of the Phase II field deployment evaluation for repair and recalibration, and the performance of the sensor improved significantly. In Phase II, the analyzers were deployed in field conditions at three diferent USGS sites. The measured nitrate concentrations were compared to discrete (reference) samples analyzed by the Direct UV method on a Shimadzu UV1800 bench top spectrophotometer, and by the National Environmental Methods Index (NEMI) method I-2548-11 at the USGS National Water Quality Laboratory. The first deployment at USGS site 0249620 on the East Pearl River in Hancock County, Mississippi, tested the ability of the TriOs ProPs (10-mm path length), Hach NITRATAX (5 mm), Satlantic SUNA (10 mm), and the S::CAN Spectro::lyser (5 mm) to accurately measure low-level (less than 2 mg-N/L) nitrate concentrations while observing the effect turbidity and colored dissolved organic matter (CDOM) would have on the analyzers' measurements. The second deployment at USGS site 01389005 Passaic River below Pompton River at Two Bridges, New Jersey, tested the analyzer's accuracy in mid-level (2-8 mg-N/L) nitrate concentrations. This site provided the means to test the analyzers' performance in two distinct matrices—the Passaic and the Pompton Rivers. In this deployment, three instruments tested in Phase I (TriOS, Hach, and SUNA) were deployed with the S::CAN Spectro::lyser (35 mm) already placed by the New Jersey Water Science Center (WSC). The third deployment at USGS site 05579610 Kickapoo Creek at 2100E Road near Bloomington, Illinois, tested the ability of the analyzers to measure high nitrate concentrations (greater than 8 mg-N/L) in turbid waters. For Kickapoo Creek, the HIF provided the TriOS (10 mm) and S::CAN (5 mm) from Phase I, and a SUNA V2 (5 mm) to be deployed adjacent to the Illinois WSC-owned Hach (2 mm). A total of 40 discrete samples were collected from the three deployment sites and analyzed. The nitrate concentration of the samples ranged from 0.3–22.2 mg-N/L. The average absolute difference between the TriOS measurements and discrete samples was 0.46 mg-N/L. For the combined data from the Hach 5-mm and 2-mm analyzers, the average absolute difference between the Hach samples and the discrete samples was 0.13 mg-N/L. For the SUNA and SUNA V2 combined data, the average absolute difference between the SUNA samples and the discrete samples was 0.66 mg-N/L. The average absolute difference between the S::CAN samples and the discrete samples was 0.63 mg-N/L.

Open-File Report

Geology and water resources of the Wharton Tract and the Mullica River basin in southern New Jersey

The Wharton Tract is an area of 150 square miles located in the Mullica River basin in southern New Jersey's Pine Barrens region. The tract is a relatively flat, low-lying, generally sandy area containing shallowly incised streams. The larger streams are commonly bordered by swamps. The tract was purchased by the State primarily as a water-supply preserve, but also for conservation and recreational purposes. Mean streamflow at three continuous record gaging stations in the Wharton Tract is: Mullica River near Batsto (46.1-square mile drainage area), 6 mgd (million gallons per day); Batsto River at Batsto (70.5-square mile drainage area), 81 mgd; Oswego River at Harrisville (64.0-square mile drainage area), 55 mgd. Thus mean discharge from this 180.6-sq mi (square mile) area is 202 mgd or 1.12 mgd per sq mi. Principal aquifers in the Wharton Tract and Mullic River basin are in the Kirkwood Formation of middle Miocene age, Cohansey Sand of Miocene(?) and Pliocene(?) age, and in overlying hydraulically connected deposits of Quaternary age. The Kirkwood Formation is composed of sand, silt, and clay. Diverse lithologies represent deposition in different environments such as nearshore marine, barrier bar, lagoonal, estuarine, and tidal marsh. Hydraulic characteristics of the Kirkwood in the Mullica River basin are virtually unknown. Most Kirkwood aquifers in the basin are believed to be hydraulically connected with the overlying Cohansey Sand. The Cohansey Sand is dominantly a quartz sand containing minor amounts of pebbly sand, silty sand, and interbedded clay. Almost all of the Wharton Tract and most of the Mullica River basin lie within the sandier area of the Cohansey, which contains approximately 75 percent sand beds and 25 percent silt and clay beds. Data from test drilling show that the upper 100 feet of sediments in the Wharton Tract contain about 93 percent sand beds, 3.5 percent clay beds, and 3.5 percent silt beds. The Cohansey in the Tract ranges in thickness from less than 50 feet to about 180 feet; the average thickness is about 125 feet. The Cohansey Sand is believed to be, in overall aspect, a deltaic deposit. It contains materials that were deposited locally in nearshore-marine, fluvial, estuarine, lagoonal, and beach environments. The hydraulic conductivity of Cohansey aquifer material ranges from about 90 to 250 feet per day (660 to 1,885 gallons per day per square foot) in southern New Jersey. One aquifer test in the Wharton Tract gives an average value of 130 feet per day (1,000 gallons per day per square foot). The transmissivity of the Cohansey Sand aquifer through most of the Wharton Tract is typically between 10,000 and 20,000 square feet per day (75,000 and 150,000 gallons per day per foot). Deposits of Quaternary age form a discontinuous veneer lying unconformably above the Cohansey Sand. The most important hydrologic function of most of these deposits is to absorb precipitation and transmit the water to the underlying Cohansey Sand. Thicker deposits of estuarine sand and clay of the Cape May Formation fill a channel in the underlying Cohansey Sand along the lower reaches of the Mullica River. This channel deposit is 85 feet thick near Batsto. Ground water and surface water in the Mullica River basin are low in dissolved solids, generally less than 50 mg/1 (milligrams per liter). Iron concentrations are generally high, up to 49,000 micrograms per liter (49 mg/1) in ground water and up to 7,100 micrograms per liter (7.1 mg/1) in the streams. The water is acidic as indicated by typical pH values of from 4.5 to 6.5. Color of the surface water is commonly high, ranging from 3 - 150 platinum-cobalt units. After appropriate treatment these waters are acceptable for most uses. The Wharton Tract is well situated to support the growing water needs of nearby New Jersey communities. Maximum development of water can be achieved by conjunctive use of ground and surface water. During most of the year, some water would be withdrawn either directly from streams or from adjacent wells. During periods of low flow during summer or fall water would be pumped from wells farther from the stream~. From analysis of flow-duration curves it is estimated that 70 mgd of water could be developed with minimal effect upon low flows in the half of the tract above the gaging stations on the Mullica and Batsto Rivers. The quantity available in the entire tract is greater, possibly in the order of 150 mgd. With augmentation of streamflow by pumping from ground water, it is likely that considerably more water could be safely used on a perennial basis. A possibility exists for multiple use of the water resources of the Mullica River through construction of an inexpensive tide barrier at the Garden State Parkway. This would create a fresh water lake in a State forest, park, and recreational area, which would also provide a flexible and economical water supply for much of the Atlantic coastal resort development.

New Jersey

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

Naturally occurring contaminants in the Piedmont and Blue Ridge crystalline-rock aquifers and Piedmont Early Mesozoic basin siliciclastic-rock aquifers, eastern United States, 1994–2008

Groundwater quality and aquifer lithologies in the Piedmont and Blue Ridge Physiographic Provinces in the eastern United States vary widely as a result of complex geologic history. Bedrock composition (mineralogy) and geochemical conditions in the aquifer directly affect the occurrence (presence in rock and groundwater) and distribution (concentration and mobility) of potential naturally occurring contaminants, such as arsenic and radionuclides, in drinking water. To evaluate potential relations between aquifer lithology and the spatial distribution of naturally occurring contaminants, the crystalline-rock aquifers of the Piedmont and Blue Ridge Physiographic Provinces and the siliciclastic-rock aquifers of the Early Mesozoic basin of the Piedmont Physiographic Province were divided into 14 lithologic groups, each having from 1 to 16 lithochemical subgroups, based on primary rock type, mineralogy, and weathering potential. Groundwater-quality data collected by the U.S. Geological Survey (USGS) National Water-Quality Assessment (NAWQA) Program from 1994 through 2008 from 346 wells and springs in various hydrogeologic and land-use settings from Georgia through New Jersey were compiled and analyzed for this study. Analyses for most constituents were for filtered samples, and, thus, the compiled data consist largely of dissolved concentrations. Concentrations were compared to criteria for protection of human health, such as U.S. Environmental Protection Agency (USEPA) drinking water maximum contaminant levels and secondary maximum contaminant levels or health-based screening levels developed by the USGS NAWQA Program in cooperation with the USEPA, the New Jersey Department of Environmental Protection, and Oregon Health & Science University. Correlations among constituent concentrations, pH, and oxidation-reduction (redox) conditions were used to infer geochemical controls on constituent mobility within the aquifers. Of the 23 trace-element constituents evaluated, arsenic, manganese, and zinc were detected in one or more water samples at concentrations greater than established human health-based criteria. Arsenic concentrations typically were less than 1 microgram per liter (µg/L) in most groundwater samples; however, concentrations of arsenic greater than 1 µg/L frequently were detected in groundwater from clastic lacustrine sedimentary rocks of the Early Mesozoic basin aquifers and from metamorphosed clastic sedimentary rocks of the Piedmont and Blue Ridge crystalline rock aquifers. Groundwater from these rock units had elevated pH compared to other rock units evaluated in this study. Of the nine samples for which arsenic concentration was greater than 10 µg/L, six were classified as oxic and three as anoxic, and seven had pH of 7.2 or greater. Manganese concentrations typically were less than 10 µg/L in most samples; however, 8.3 percent of samples from the Piedmont and Blue Ridge crystalline-rock aquifers and 3.0 percent of samples from the Early Mesozoic basin siliciclastic rock aquifers had manganese concentrations greater than the 300-µg/L health-based screening level. The positive correlation of manganese with iron and ammonia and the negative correlation of manganese with dissolved oxygen and nitrate are consistent with the reductive dissolution of manganese oxides in the aquifer. Zinc concentrations typically were less than 10 µg/L in the groundwater samples considered in the study, but 0.4 percent and 5.5 percent of the samples had concentrations greater than the health-based screening level of 2,000 µg/L and one-tenth of the health-based screening level, respectively. The mean rank concentration of zinc in groundwater from the quartz-rich sedimentary rock lithologic group was greater than that for other lithologic groups even after eliminating samples collected from wells constructed with galvanized casing. Approximately 90 percent of 275 groundwater samples had radon-222 concentrations that were greater than the proposed alternative maximum contaminant level of 300 picocuries per liter. In contrast, only 2.0 percent of 98 samples had combined radium (radium-226 plus radium-228) concentrations greater than the maximum contaminant level of 5.0 picocuries per liter, and 0.6 percent of 310 samples had uranium concentrations greater than the maximum contaminant level of 30 µg/L. Radon concentrations were highest in the Piedmont and Blue Ridge crystalline-rock aquifers, especially in granite, and elevated median concentrations were noted in the Piedmont Early Mesozoic basin aquifers, but without the extreme maximum concentrations found in the crystalline rocks (granites). Although the siliciclastic lithologies had a greater frequency of elevated uranium concentrations, radon and radium were commonly detected in water from both siliciclastic and crystalline lithologies. Uranium concentrations in groundwater from clastic sedimentary and clastic lacustrine/evaporite sedimentary lithologic groups within the Early Mesozoic basin aquifers, which had median concentrations of 3.6 and 3.1 µg/L, respectively, generally were higher than concentrations for other siliciclastic lithologic groups, which had median concentrations less than 1 µg/L. Although 89 percent of the 260 samples from crystalline-rock aquifers had uranium concentrations less than 1 µg/L, 0.8 percent had uranium concentrations greater than the 30-µg/L maximum contaminant level, and 6.5 percent had concentrations greater than 3 µg/L.

Alabama;Delaware;Georgia;Maryl;New Jersey;North Ca

Hydrogeologic framework, ground-water quality, and simulation of ground-water flow at the Fair Lawn Well Field Superfund site, Bergen County, New Jersey

Production wells in the Westmoreland well field, Fair Lawn, Bergen County, New Jersey (the 'Fair Lawn well field Superfund site'), are contaminated with volatile organic compounds, particularly trichloroethylene, tetrachloroethylene, and 1,1,1-trichloroethane. In 1983, the U.S. Environmental Protection Agency (USEPA) placed the Westmoreland well field on its National Priority List of Superfund sites. In an effort to determine ground-water flow directions, contaminant-plume boundaries, and contributing areas to production wells in Fair Lawn, and to evaluate the effect of present pump-and-treat systems on flowpaths of contaminated ground water, the U.S. Geological Survey (USGS), in cooperation with the USEPA, developed a conceptual hydrogeologic framework and ground-water flow model of the study area. MODFLOW-2000, the USGS three-dimensional finite-difference model, was used to delineate contributing areas to production wells in Fair Lawn and to compute flowpaths of contaminated ground water from three potential contaminant sources to the Westmoreland well field. Straddle-packer tests were used to determine the hydrologic framework of, distribution of contaminants in, and hydrologic properties of water-bearing and confining units that make up the fractured-rock aquifer underlying the study area. The study area consists of about 15 square miles in and near Fair Lawn. The area is underlain by 6 to 100 feet of glacial deposits and alluvium that, in turn, are underlain by the Passaic Formation. In the study area, the Passaic Formation consists of brownish-red pebble conglomerate, medium- to coarse-grained feldspathic sandstone, and micaceous siltstone. The bedrock strata strike N. 9o E. and dip 6.5o to the northwest. The bedrock consists of alternating layers of densely fractured rocks and sparsely fractured rocks, forming a fractured-rock aquifer. Ground-water flow in the fractured-rock aquifer is anisotropic as a result of the interlayering of dipping water-bearing and confining units. Wells of similar depth aligned along the strike of the bedding intersect the same water-bearing units, but wells aligned along the dip of the bedding may intersect different water-bearing units. Consequently, wells aligned along strike are in greater hydraulic connection than wells aligned along dip. The Borough of Fair Lawn pumps approximately 770 million gallons per year from 13 production wells. Hydrographs from six observation wells ranging in depth from 162 to 505 feet in Fair Lawn show that water levels in much of the study area are affected by pumping. Straddle packers were used to isolate discrete intervals within six open-hole observation wells owned by the Fair Lawn Water Department. Transmissivity, water-quality, and static-water-level data were obtained from the isolated intervals. Measured transmissivity ranged from near 0 to 8,900 feet squared per day. The broad range in measured transmissivity is a result of the heterogeneity of the fractured-rock aquifer. Eight water-bearing units and eight confining units were identified in the study area on the basis of transmissivity. The water-bearing units range in thickness from 21 to 95 feet; the mean thickness is 50 feet. The confining units range in thickness from 22 to 248 feet; the mean thickness is 83 feet. Water-level and water-quality data indicate effective separation of water-bearing units by the confining units. Water-quality samples were collected from the six observation wells at 16 depth intervals isolated by the straddle packers in 2000 and 2001. Concentrations of volatile organic compounds generally were low in samples from four of the wells, but were higher in samples from a well in Fair Lawn Industrial Park and in a well in the Westmoreland well field. The digital ground-water flow model was used to simulate steady-state scenarios representing conditions in the study area in 1991 and 2000. These years were chosen because during the intervening period,

Scientific Investigations Report

Hydrologic assessment of the Edwin B. Forsythe National Wildlife Refuge

The Edwin B. Forsythe National Wildlife Refuge (hereafter Forsythe refuge or the refuge) is situated along the central New Jersey coast and provides a mixture of freshwater and saltwater habitats for numerous bird, wildlife, and plant species. Little data and information were previously available regarding the freshwater dynamics that support the refuge’s ecosystems. In cooperation with the U.S. Fish and Wildlife Service, the U.S. Geological Survey conducted an assessment of the hydrologic resources and processes in the refuge and surrounding areas to provide baseline information for evaluating restoration projects and future changes in the hydrologic system associated with climate change and other anthropogenic stressors. During spring 2015, water levels were measured at groundwater and surface-water sites in and near the Forsythe refuge. These water-level measurements, along with surface-water elevations obtained from digital elevation models, were used to construct water-table-elevation and depth-to-water maps of the refuge and surrounding areas. Water-table elevations in the refuge ranged from sea level to approximately 65 feet above sea level; in most of the refuge, the water-table elevation was within 3 feet of sea level. The water-table-elevation map indicates that the direction of shallow groundwater flow at the regional scale is generally from west to east (much of it from the northwest to the southeast), and groundwater moves downgradient from the uplands toward major groundwater discharge areas consisting of coastal streams and wetlands. The depth to water is estimated to be less than 2 feet for approximately 86 percent of the refuge, which coincides closely with the percentage of wetland area in the refuge. Depth to water in excess of 20 feet below land surface is limited to higher elevation areas of the refuge. Streamflow data collected at continuous-record streamgages and partial-record stations within the Mullica-Toms Basin were summarized. Hydrograph separation of streamflow data for eight streamgages (2004–13) reveals that base flow accounts for 68–94 percent of streamflow in basins upstream from the refuge. The high base-flow inputs underscore the importance of groundwater as a source of freshwater that supports both the streams that flow into the refuge and the hydroecology of the contributing basins. Mean annual flow typically ranged from 1.7 to 2.1 cubic feet per second per square mile at the streamgages (2004–13) and between 1.2 and 2.3 cubic feet per second per square mile at the partial-record stations (1965–2015) but was notably greater or lower than these ranges at several stations. Mean annual water budgets were estimated for multiple regions of the refuge for 2004–13 using data compiled from nearby meteorological stations and groundwater flows derived from previously calibrated groundwater-flow models. Precipitation, groundwater recharge, and evapotranspiration were estimated from available data; direct runoff was calculated as the residual component of the water balance. Groundwater recharge rates were greatest in the upland-dominated areas of the refuge with estimates of 14.4 to 18.9 inches per year, which are equivalent to 30 to 40 percent of precipitation. Groundwater recharge rates were nearly zero in the central coastal areas because these areas are major groundwater discharge zones, the water table is near land surface, the subsurface is close to saturation and cannot accept much recharge, and much of the area is underlain by thick marsh deposits likely with low permeability. Estimates of evapotranspiration varied from about 26 inches per year in the upland-dominated areas to more than 35 inches per year in the coastal wetlands, equivalent to 55–79 percent of mean annual precipitation, indicating that it is a major component of the hydrodynamics of the Forsythe refuge. On the basis of output from previously calibrated groundwater-flow models, nearly all of the groundwater exiting the surficial aquifer system in the central coastal areas of the refuge is discharged to wetlands, which highlights the importance of groundwater discharge in supporting the ecosystems of the Forsythe refuge. In the central coastal areas, horizontal flow contributes more than 90 percent of the groundwater flow to the surficial system, indicating that the upbasin areas are a substantial source of water that ultimately discharges to streams and wetlands in the refuge.

New Jersey

High-water marks from Hurricane Sandy for coastal areas of Connecticut, Rhode Island, and Massachusetts, October 2012

Because coastal areas in Connecticut, Rhode Island, and Massachusetts were heavily affected by Hurricane Sandy in October 2012, the U.S. Geological Survey (USGS), under a mission agreement with the Federal Emergency Management Agency, collected storm tide high-water marks in those coastal areas. This effort was undertaken to better understand the areal extent and impact of storm tides resulting from strong storms. On October 27–29, 2012, Hurricane Sandy moved up the eastern coast of the United States after passing over the Bahamas. On October 29 at about 7:30 p.m. eastern daylight time, Hurricane Sandy made landfall its final time near Brigantine, New Jersey, with recorded wind speeds of about 80 miles per hour. The damages from Hurricane Sandy exceeded $50 billion in total, making it the second most costly Atlantic hurricane at that time, second only to Hurricane Katrina in 2005. Hurricane Sandy also resulted in 147 deaths, and about 650,000 homes and many businesses being damaged along the eastern coast of the United States. The severity of Hurricane Sandy’s effects resulted in presidential disaster declarations being declared in 10 States from Virginia to Massachusetts and the District of Columbia in the months following Hurricane Sandy; the list of States affected included Connecticut, Rhode Island, and Massachusetts. In response to the approach of Hurricane Sandy, the USGS deployed 60 temporary storm tide sensors and 2 temporary real-time rapid deployment gages to collect tide elevation data during the storm along the coastal areas of Connecticut, Rhode Island, and Massachusetts. This activity was done from Virginia to Maine before the storm. Following Hurricane Sandy, in October and November 2012, 371 storm tide high-water marks were identified and flagged in the coastal areas of Connecticut, Rhode Island, and Massachusetts. High-water marks near USGS temporary storm tide sensors, real-time rapid deployment gages, and streamgages affected by the tides as well as high-water marks on Block Island, R.I., and Martha’s Vineyard and Nantucket, Mass., were surveyed at the same time the high-water marks were identified and flagged in October and November 2012. The remaining high-water marks flagged during October and November 2012 were surveyed from December 2013 through June 2014 and in December 2016. Elevations of all high-water marks were referenced to the North American Vertical Datum of 1988 and horizontal coordinates to the North American Datum of 1983 using the Global Navigation Satellite System, survey-grade Digital Global Positioning System receivers, and total station surveying equipment. Of the 371 storm tide high-water marks flagged following Hurricane Sandy, only 364 high-water marks were surveyed; the remaining 7 could not be found or had been destroyed when locations were revisited to conduct surveys. The 157 high-water marks surveyed in Connecticut had elevations that ranged from 2.5 to 12.2 feet (ft) with an average elevation of 8.1 ft and a median elevation of 8.3 ft. The 76 high-water marks in Rhode Island had elevations that ranged from 3.6 to 16.2 ft and averaged 7.1 ft with a median of 6.6 ft. The 131 high-water marks in Massachusetts had elevations that ranged from 2.8 to 22.7 ft and averaged 7.3 ft with a median of 6.6 ft. Individual information on the location, type, accuracy, and elevation of the 371 high-water marks can be found in an accompanying USGS data release and at the USGS Flood Event Viewer website for Hurricane Sandy ( https://stn.wim.usgs.gov/fev/#Sandy ). The high-water marks along the coast line of Connecticut and eastern Massachusetts, including Nantucket, generally had higher storm tide elevations than the coast line of Rhode Island including Block Island and southern Massachusetts, including Martha’s Vineyard. The high-water mark elevations compare well with recorded peak-storm tide data at USGS temporary storm tide sensors and real-time rapid deployment gages deployed for Hurricane Sandy in Connecticut, Rhode Island, and Massachusetts. High-water mark data collected following Hurricane Sandy will be used by Federal, State, and local government agencies, nongovernmental organizations, universities, and the public for better understanding the areal extent and impact of the storm tides. Additionally, these data can be used for such activities as land-use planning, flood risk studies, flood resiliency studies, and coastal modeling. These data from this historic storm can be compared with other regional hurricanes and tropical storms for planning into the future.

Connecticut, Massachusetts, Rhode Island

Trends in chemical concentration in sediment cores from three lakes in New Jersey and one lake on Long Island, New York

Sediment cores were extracted from three lakes in northeastern New Jersey and one lake on western Long Island, New York, as part of the U.S. Geological Survey National Water-Quality Assessment Program. Sediment layers were dated by use of cesium-137 (137Cs), copper, lead, or dichlorodiphenyl-trichloroethane (DDT) profiles. Sediment layers were analyzed for seven selected trace elements, including arsenic, cadmium, chromium, lead, mercury, nickel, and zinc, and five hydrophobic organochlorine compounds, including chlordane, dieldrin, total DDT, total polychlorinated biphenyls (PCBs), and total polycyclic aromatic hydrocarbons (PAHs). All seven trace elements were detected throughout the cores from all four lakes. Concentrations of all elements, except arsenic, were elevated in the three cores from lakes within urbanized watersheds (Packanack Lake, Orange Reservoir, and Newbridge Pond) relative to the concentrations in the lake core collected below the largely forested, reference watershed (Clyde Potts Reservoir). Results of trend analyses indicate that concentrations of all trace elements, with the exception of arsenic and lead, were relatively constant throughout the core from the minimally urbanized Clyde Potts Reservoir. In urban lakes, significant upward trends in concentrations from deeper to shallower sediments were observed either to peak concentrations or throughout the core for all elements, with the exception of chromium at all lakes and arsenic and nickel at Orange Reservoir. This finding indicates that changes in population and land use in the urbanized watersheds over the period of sedimentary record have contributed to upward trends in trace-element concentrations. Although downward trends in concentrations were observed for some trace elements in the years after their concentrations peaked, concentrations of all trace elements in urban lake cores were higher in the most recently deposited sediments than at the base of each respective core. Lead concentrations over time were highly correlated with the population in the vicinity of the lake until the concentration peak in sediment deposited in the mid-1970’s. Concentrations of lead in lake sediment appear to be closely related to the use of leaded gasoline because lead concentrations generally decreased after the use of leaded gasoline was phased-out in the mid-1970’s. Zinc concentrations were highly correlated with population over the entire length of the core. In general, zinc concentrations increased in the three urbanized watersheds, probably in response to increasing population and vehicular use. This trend was not evident at Clyde Potts Reservoir, however, where vehicular traffic in the watershed is minimal. Detectable concentrations of chlordane, total DDT, and total PCBs were present in cores from all lakes; however, dieldrin was detected only in the Newbridge Pond and Packanack Lake cores. Concentrations generally were higher in cores from the urbanized Newbridge Pond and Orange Reservoir watersheds than in those from the minimally urbanized Clyde Potts Reservoir watershed. With the exception of chlordane in the Clyde Potts and Orange Reservoir cores, concentrations of the four organochlorine compounds had significant downward trends from peak concentrations to recently deposited sediment or non-significant trends throughout the core. On the basis of these findings and as a result of regulatory actions prohibiting the production and use of these compounds, downward trends in sedimentary concentrations are expected to continue; however, the persistence of these compounds indicates that a substantial amount of time may be required to purge them from the watersheds. Concentrations of PAHs in sediment generally increased with population growth and urbanization, probably as a result of increased fossil-fuel combustion (gasoline and home-heating fuels and other uses (roads and parking lots paved with asphalt) associated with increased urban development and vehicular traffic. This finding is supported by low concentrations of PAHs in Packanack Lake sediments in the 1930’s, before the watershed was urbanized and when automobiles were comparatively rare. As vehicular use and urbanization increase in these watersheds, the general increase of PAH concentrations in lake sediments can be expected to continue. Data from this study indicate that changes in population, land use, and chemical use in the urbanized watersheds over the period of sedimentary record have contributed to upward trends in concentrations of trace elements and hydrophobic organic compounds. Although downward trends were observed for some constituents in the years after their concentrations peaked, concentrations of most constituents in urban lake cores were higher in the most recently deposited sediments than at the base of each respective core and in the reference lake cores. Similar trends in concentrations of these constituents have been observed in sediment cores from other urban lakes across the United States.

New Jersey, New York

A comparison of synthetic flowpaths derived from light detection and ranging topobathymetric data and National Hydrography Dataset High Resolution Flowlines

Bathymetric and topobathymetric light detection and ranging (lidar) digital elevation models created for the Delaware River were provided to the National Geospatial Program and used to evaluate synthetic flowpath extraction from bathymetric/topobathymetric lidar survey data as a data source for improving the density, distribution, and connectivity of the National Hydrography Dataset High Resolution Flowline Network. As the surface-water component of The National Map, the National Hydrography Dataset maintains the Nation’s drainage network flow information and geometries for surface-water features used in hydrologic, hydraulic, and other science and engineering disciplines. The regional lidar survey for the Delaware River between Hancock, New York, and Trenton, New Jersey, was collected for the U.S. Geological Survey using the Experimental Advanced Airborne Research Lidar sensor system and processed by the Coastal National Elevation Database Applications Program. Using 1 percent of the maximum flow accumulation value for the surveyed Delaware River corridor as the flow accumulation threshold for grid cells at 1-, 5-, and 10-meter resolution created 223 to 283 kilometers of synthetic flowpaths potentially representing the river channel thalweg, which is the deepest point in a riverbed cross-section. There was potential for improving the High Resolution National Hydrography Dataset (HR NHD) Flowline network in places where the Delaware River channel, depicted as an Artificial Path in the HR NHD, is offset from the extracted synthetic river flowpath which sometimes appeared better positioned than the Artificial Path to represent the river thalweg. For the same area, using 0.05 percent of the maximum flow accumulation at the 1-, 5-, and 10-meter resolutions extracted 744 to 1,317 kilometers of synthetic flowpaths, with extracted synthetic flowpaths representing the main river channel and additional synthetic flowpaths representing tributaries or streams adjacent to the main channel. Overlaying these results with the HR NHDFlowline Network indicates that some of the additional synthetic flowpaths are connected to or extend HR NHD stream/river feature types. Some disconnected or isolated synthetic flowpaths not included in stream/river feature types were validated in orthoimagery and U.S. Topo Maps and provide examples of how extracted synthetic flowpaths could be used to delineate new stream/river features. Other additional extracted synthetic flowpaths depict linear features such as canals, tree lines, roads, or linear topographic depressions. For some river reaches where obstructions to flow or where low-relief topographic or bathymetric surfaces alter the flow direction, the software tool used to develop the flow direction grid did not calculate a primary flowpath for the river channel. Based on the results of this analysis, site conditions for the Delaware River corridor did not affect the quality of lidar bathymetric survey data. However, depending on the resolution of the lidar bathymetric digital elevation models (BDEMs), site conditions do have different effects on results for extracted synthetic flowpaths. We found that synthetic flowpaths extracted from 1-meter resolution lidar DEMs had more varied flow directions around in-channel landforms that obstructed flow than synthetic flowpaths extracted from 5- or 10-meter resolution lidar DEMs. As a result the 1-meter resolution DEM created some isolated or discontinuous synthetic flowpath segments where the 5- and 10-meter DEMs developed more continuous flowpaths. In this case the river bed upstream from the in-channel obstruction is shallower than the river bed downstream. Under these conditions the 1-meter resolution DEM provided synthetic flowpaths delineating a potential river thalweg. In this same area, the software solution modified (virtually raised) the river bed in the 5- and 10-meter resolution DEMs and flattened the bathymetric surface to create a continuous downstream flow direction, which caused trellis-patterned synthetic flowpaths to form. Under different site conditions and converse to the above development of synthetic flowpaths at different resolutions, at an abandoned river flood plain (terrace) with low relief that is adjacent to the river channel, the flow direction grid for the 1-meter resolution DEM developed continuous synthetic flowpath corresponding to a HR NHD Flowline network stream/river feature that connected to the main river channel but the larger resolution DEMs created isolated or disconnected synthetic flowpaths. A project to continue an evaluation of benefits of or issues caused by extracting synthetic flowpaths to enhance the HR NHD could include a study to assess the potential for merging surface-water flowpaths extracted from lidar topobathymetry and 3D Elevation Program digital elevation models. The merged DEM approach to synthetic flowpath extraction could extend the HR NHDFlowline network and enhance flow accumulations that might develop better flow direction grids in low-relief areas. Because of the confined lateral extent of the Delaware River, the lidar DEMs were not used to create catchments or watersheds; however, the merged DEM approach could also be tested as a resource for enhancing HR NHD catchments and watersheds. This lidar DEM synthetic flowpath extraction project supports the National Geospatial Program efforts to collect and produce high-quality lidar data to provide 3-dimensional representations of natural feature and aligns with the National Spatial Data Infrastructure to improve utilization of geospatial data. The results also can be useful for understanding strategies that can help maintain quality data in the HR NHD programs. KEYWORDS: bathymetric, digital elevation model, extracted synthetic flowpath, lidar, High Resolution National Hydrography Dataset, topobathymetric

New Jersey

Simulated effects of alternative withdrawal strategies on groundwater flow in the unconfined Kirkwood-Cohansey aquifer system, the Rio Grande water-bearing zone, and the Atlantic City 800-foot sand in the Great Egg Harbor and Mullica River Basins, New Jersey

Groundwater is essential for water supply and plays a critical role in maintaining the environmental health of freshwater and estuarine ecosystems in the Atlantic Coastal basins of New Jersey. The unconfined Kirkwood-Cohansey aquifer system and the confined Atlantic City 800-foot sand are major sources of groundwater in the area, and each faces different water-supply concerns. The U.S. Geological Survey (USGS), in cooperation with the New Jersey Department of Environmental Protection (NJDEP), conducted a study to simulate the effects of withdrawals in the Kirkwood-Cohansey aquifer system, the Atlantic City 800-foot sand, and the Rio Grande water-bearing zone and to evaluate potential scenarios. The study area encompasses Atlantic County and parts of Burlington, Camden, Gloucester, Ocean, Cape May, and Cumberland Counties. The major hydrogeologic units affecting water supply in the study area are the surficial Kirkwood-Cohansey aquifer system, a thick diatomaceous clay confining unit in the upper part of Kirkwood Formation; the Rio Grande water-bearing zone; and the Atlantic City 800-foot sand of the Kirkwood Formation. Hydrogeologic data from 18 aquifer tests and specific capacity data from 230 wells were analyzed to provide horizontal hydraulic conductivity of the aquifers. Groundwater withdrawals are greatest from the Kirkwood-Cohansey aquifer system, and 65 percent of the water is used for public supply. Groundwater withdrawals from the Atlantic City 800-foot sand are about half those from the Kirkwood-Cohansey aquifer system. Ninety-five percent of the withdrawals from the Atlantic City 800-foot sand is used for public supply. Data from six streamgaging stations and 51 low-flow partial record sites were used to estimate base flow in the area. Base flow ranges from 60 to 92 percent of streamflow. A groundwater flow model of the Kirkwood-Cohansey aquifer system, the Rio Grande water-bearing zone, and the Atlantic City 800-foot sand was developed and calibrated using water-level data from 148 wells and base-flow data from 22 gaging or low-flow partial record stations. The Kirkwood-Cohansey aquifer system within the Great Egg Harbor River and the Mullica River Basins was simulated on a monthly basis from 1998 through 2006. An existing regional model of the New Jersey Coastal Plain was revised to provide boundary conditions for the Great Egg Harbor and Mullica River Basin model (referred to as the Great Egg-Mullica model). In the Great Egg-Mullica model, monthly groundwater recharge rates used in the model ranged from 10-15 inches per year in 2001 to 20-25 inches per year in 2005. The mean-absolute error for 10 of the 14 long-term hydrographs used in model calibration was less than 5 ft. Groundwater flow budgets for the Great Egg-Mullica model calibration periods, May 2005 and September 2006, and for the entire model calibration period 1998 to 2006, showed that nearly 70 percent of the water entering the Atlantic City 800-foot sand came from the horizontal connection with the Kirkwood-Cohansey aquifer system in updip areas. The groundwater flow model was used to simulate scenarios under three possible conditions: average 1998 to 2006 withdrawals (Average scenario), full-allocation withdrawals (Full Allocation scenario), and projected 2050-demand withdrawals (2050 Demand scenario). Withdrawals in the Full Allocation scenario are nearly twice the withdrawals from the Average scenario, primarily because of the potential for large agricultural withdrawals if all allocations are used. Withdrawals for the 2050 Demand scenario are about 50 percent greater than those for the Average scenario, primarily due to expected increases in withdrawals for public supply. Monthly base-flow depletion criteria were determined using the Low-Flow Margin method, currently under consideration by NJDEP, to estimate available water on an annual basis at the Hydrologic Unit Code 11 (HUC11) level and to determine whether a water-supply deficit exists. Simulations of various groundwater-withdrawal scenarios were made using the calibrated model, and results were compared with baseline conditions (no withdrawals) to determine where and when base-flow deficits may be occurring and may be expected to occur in the future. Scenarios were simulated to assess base-flow depletion that could occur from different groundwater-withdrawal situations. In the Average scenario, deficits occurred in 7 of the 14 subbasins. In the Full Allocation scenario, deficits occurred in 11 of the subbasins. In the 2050 Demand scenario, deficits occurred in 9 of the 14 subbasins. The largest deficits occurred in the Absecon Creek subbasin because the base-flow depletion criteria for this subbasin is small due to the surface-water diversions that are already occurring there and because existing groundwater withdrawals in the subbasin have resulted in base-flow depletion under current (1998-2006) conditions. Three adjusted scenarios, variations of the Average, Full Allocation, and 2050 Demand scenarios, were simulated; for the adjusted scenarios, the withdrawals were modified in stages with the intent to successively eliminate or minimize the base-flow deficits. Modifications included shifting withdrawals to a deeper part of the Kirkwood-Cohansey aquifer system, implementing seasonal conjunctive use of shallow and deep aquifers, and specifying reductions in withdrawals within a HUC11 subbasin in deficit. The adjusted scenarios are intended to show the relative effectiveness of each of the three approaches in reducing the deficits. Most of the deficits under the Average, Full Allocation, and 2050 Demand scenarios were eliminated by reductions in withdrawals or allocations. Shifting withdrawals to a deeper part of the Kirkwood-Cohansey aquifer system or seasonal conjunctive use did not eliminate deficits for any subbasin. Reductions in withdrawals accounted for more than 95 percent of the total reduction of deficits in all but one subbasin.

New Jersey

Quality of water from crystalline rock aquifers in New England, New Jersey, and New York, 1995-2007

Crystalline bedrock aquifers in New England and parts of New Jersey and New York (NECR aquifers) are a major source of drinking water. Because the quality of water in these aquifers is highly variable, the U.S. Geological Survey (USGS) statistically analyzed chemical data on samples of untreated groundwater collected from 117 domestic bedrock wells in New England, New York, and New Jersey, and from 4,775 public-supply bedrock wells in New England to characterize the quality of the groundwater. The domestic-well data were from samples collected by the USGS National Water-Quality Assessment (NAWQA) Program from 1995 through 2007. The public-supply-well data were from samples collected for the U.S. Environmental Protection Agency (USEPA) Safe Drinking Water Act (SDWA) Program from 1997 through 2007. Chemical data compiled from the domestic wells include pH, specific conductance, dissolved oxygen, alkalinity, and turbidity; 6 nitrogen and phosphorus compounds, 14 major ions, 23 trace elements, 222 radon gas (radon), 48 pesticide compounds, and 82 volatile organic compounds (VOCs). Additional samples were collected from the domestic wells for the analysis of gross alpha- and gross beta-particle radioactivity, radium isotopes, chlorofluorocarbon isotopes, and the dissolved gases methane, carbon dioxide, nitrogen, and argon. Chemical data compiled from the public-supply wells include pH, specific conductance, nitrate, iron, manganese, sodium, chloride, fluoride, arsenic, uranium, radon, combined radium ( 226 radium plus 228 radium), gross alpha-particle radioactivity, and methyl tert -butyl ether (M t BE). Patterns in fluoride, arsenic, uranium, and radon distributions were discernable when the data were compared to lithology groupings of the bedrock, indicating that the type of bedrock has an effect on the quality of groundwater from NECR aquifers. Fluoride concentrations were significantly higher in groundwater samples from the alkali granite, peraluminous granite, and metaluminous granite lithology groups than from samples in the other lithology groups. Water samples from 1.4 percent of 2,167 studied wells had fluoride concentrations that were equal to or greater than the maximum contaminant level (MCL) of 4 milligrams per liter (mg/L) and 7.5 percent of the wells had fluoride concentrations that were equal to or greater than the secondary MCL of 2 mg/L. For arsenic, groundwater samples from the calcareous metasedimentary rocks in the New Hampshire-Maine geologic province, peraluminous granite, and pelitic rocks lithology groups had higher concentrations than did samples from the other lithology groups. Water samples from 13.3 percent of 2,054 studied wells had arsenic concentrations that were equal to or greater than the MCL of 10 micrograms per liter (μg/L), about double the national rate of occurrence in community-supply systems and in domestic wells of the United States. Uranium concentrations were significantly higher in groundwater samples from the peraluminous granite, alkali granite, and calcareous metasedimentary rocks in the New Hampshire-Maine geologic province lithology groups than from samples in the other lithology groups. Water samples from 14.2 percent of 556 studied wells had uranium concentrations equal to or greater than the MCL of 30 μg/L. Radon activities were equal to or greater than the proposed MCL of 300 picocuries per liter (pCi/L) in 95 percent of 943 studied wells, and 33 percent of the wells had radon activities were equal to or greater than the proposed alternative maximum contaminant level (AMCL) of 4,000 pCi/L. Radon activities exceeded the proposed AMCL in 20 percent or more of groundwater samples in each of the studied lithology groups with a minimum of 9 samples, but radon activities were significantly higher in groundwater samples from the alkali granite, peraluminous granite, and Narragansett basin metasedimentary rocks lithology groups. Water samples from 3.2 percent of 564 studied wells had combined radium activities equal to or greater than the MCL of 5 pCi/L; however, combined radium activities were not significantly different among the studied lithology groups. Land use and population density also were evaluated to explain patterns in water quality. Concentrations of nitrate, sodium, chloride, and MtBE from the studied wells were significantly greater in areas of high population density (≥50 persons per square kilometer) than in areas of low population density (<50 persons per square kilometer). Concentrations of sodium, chloride, and M t BE from the studied wells were significantly greater in areas classified as developed (urban lands) than in areas classified as undeveloped (forested), agricultural, or mixed (no dominant land use). Nitrate concentrations from the public-supply wells were not significantly different among the four land use categories, but nitrate concentrations from the domestic wells were significantly greater in areas classified as developed than in areas classified as undeveloped, agricultural, or mixed. Chloride to bromide mass ratios in the domestic well samples indicate that the groundwater was probably affected by at least three halogen sources: local precipitation and recharge waters, remnant seawater and connate waters evolved from seawater, and recharge waters affected by road salt. The groundwater in the NECR aquifers generally contained low concentrations of nitrate, VOCs, and pesticides. Less than 1 percent of water samples from 4,781 studied wells had concentrations of nitrate greater than the MCL of 10 mg/L. Less than 1 percent of water samples from 1,299 studied wells exceeded the USEPA advisory level of 20 to 40 μg/L for M t BE. None of the other studied VOCs exceeded a human health benchmark. M t BE (36 percent frequency detection) and chloroform (32.9 percent frequency detection) were the most frequently detected (>0.02 μg/L) VOCs in the domestic wells. M t BE was detected more often in water samples with apparent ages of less than 25 years than in water samples with apparent ages greater than 25 years. This finding is consistent with the time period of high M t BE use in areas in the United States where reformulated gasoline was mandated. The largest pesticide concentration was an estimated concentration of 0.06 μg/L for the herbicide metolachlor. Deethylatrazine, a degradate of atrazine, (18 percent frequency detection) and atrazine (8 percent frequency detection) were the only pesticide compounds detected (>0.001 μg/L) in more than 3 percent of the domestic wells. None of the detected pesticide compounds exceeded human health benchmarks. Concentrations of nitrate and gross alpha-particle activities were significantly greater in the water samples from the domestic wells than in samples from the public-supply wells. Concentrations of sodium, chloride, iron, manganese, and uranium were significantly greater in the water samples from the public-supply wells than in the samples from the domestic wells. One possible explanation may be related to differences in field processing (filtered samples from the domestic wells compared to unfiltered samples from the public-supply wells). The high frequency of detections for a wide variety of manmade and naturally occurring contaminants in both domestic and public-supply wells shows the vulnerability of NECR aquifers to contamination. The highly variable water quality and the association with highly variable lithology of crystalline bedrock underscores the importance of testing individual wells to determine if concentrations for the most commonly detected contaminants exceed human health benchmarks.

Connecticut, Maine, Massachusetts, New Hampshire,

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

Occurrence of organic compounds and trace elements in the upper Passaic and Elizabeth Rivers and their tributaries in New Jersey, July 2003 to February 2004: Phase II of the New Jersey toxics reduction workplan for New York-New Jersey Harbor

Samples of surface water and suspended sediment were collected from the Passaic and Elizabeth Rivers and their tributaries in New Jersey from July 2003 to February 2004 to determine the concentrations of selected chlorinated organic and inorganic constituents. This sampling and analysis was conducted as Phase II of the New York-New Jersey Harbor Estuary Workplan—Contaminant Assessment and Reduction Program (CARP), which is overseen by the New Jersey Department of Environmental Protection. Phase II of the New Jersey Workplan was conducted to define upstream tributary and point sources of contaminants in those rivers sampled during Phase I work, with special emphasis on the Passaic and Elizabeth Rivers. Samples were collected from three groups of tributaries: (1) the Second, Third, and Saddle Rivers; (2) the Pompton and upper Passaic Rivers; and (3) the West Branch and main stem of the Elizabeth River. The Second, Third, and Saddle Rivers were sampled near their confluence with the tidal Passaic River, but at locations not affected by tidal flooding. The Pompton and upper Passaic Rivers were sampled immediately upstream from their confluence at Two Bridges, N.J. The West Branch and the main stem of the Elizabeth River were sampled just upstream from their confluence at Hillside, N.J. All tributaries were sampled during low-flow discharge conditions using the protocols and analytical methods for organic constituents used in low-flow sampling in Phase I. Grab samples of streamflow also were collected at each site and were analyzed for trace elements (mercury, methylmercury, cadmium, and lead) and for suspended sediment, particulate organic carbon, and dissolved organic carbon. The measured concentrations and available historical suspended-sediment and stream-discharge data (where available) were used to estimate average annual loads of suspended sediment and organic compounds in these rivers. Total suspended-sediment loads for 1975–2000 were estimated using rating curves developed from historical U.S. Geological Survey (USGS) suspended-sediment and discharge data, where available. Average annual loads of suspended sediment, in millions of kilograms per year (Mkg/yr), were estimated to be 0.190 for the Second River, 0.23 for the Third River, 1.00 for the Saddle River, 1.76 for the Pompton River, and 7.40 for the upper Passaic River. On the basis of the available discharge records, the upper Passaic River was estimated to provide approximately 60 percent of the water and 80 percent of the total suspended-sediment load at the Passaic River head-of-tide, whereas the Pompton River provided roughly 20 percent of the total suspended-sediment load estimated at the head-of-tide. The combined suspended-sediment loads of the upper Passaic and Pompton Rivers (9.2 Mkg/yr), however, represent only 40 percent of the average annual suspended-sediment load estimated for the head-of-tide (23 Mkg/yr) at Little Falls, N.J. The difference between the combined suspended-sediment loads of the tributaries and the estimated load at Little Falls represents either sediment trapped upriver from the dam at Little Falls, additional inputs of suspended sediment downstream from the tributary confluence, or uncertainty in the suspended-sediment and discharge data that were used. The concentrations of total suspended sediment-bound polychlorinated biphenyls (PCBs) in the tributaries to the Passaic River were 194 ng/g (nanograms per gram) in the Second River, 575 ng/g in the Third River, 2,320 ng/g in the Saddle River, 200 ng/g in the Pompton River, and 87 ng/g in the upper Passic River. The dissolved PCB concentrations in the tributaries were 563 pg/L (picograms per liter) in the Second River, 2,510 pg/L in the Third River, 2,270 pg/L in the Saddle River, 887 pg/L in the Pompton River, and 1,000 pg/L in the upper Passaic River. Combined with the sediment loads and discharge, these concentrations resulted in annual loads of suspended sediment-bound PCBs, in grams per year (g/yr), of 37 in the Second River; 132 in the Third River; 2,320 in the Saddle River; 352 in the Pompton River; and 644 in the upper Passaic River. Annual loads of dissolved PCBs, in grams per year, are 9.2 in the Second River; 47 in the Third River; 212 in the Saddle River; 349 in the Pompton River; and 549 in the upper Passaic River. Concentrations of total suspended sediment-bound polychlorinated dibenzo-p-dioxins and polychlorinated dibenzo-p-difurans (PCDD/PCDFs) were 6,000 pg/g (picograms per gram) in the Second River; 11,300 pg/g in the Third River; 37,700 pg/g in the Saddle River; 7,140 pg/g in the Pompton River; and 9,640 pg/g in the upper Passaic River. Total toxic equivalence quotients (TEQs), which included PCDD/PCDFs and coplanar PCBs, ranged from 2.7 pg/g in the Second River to 132 pg/g (as 2,3,7,8-TCDD) in the Saddle River. Average annual loads of PCDD/PCDFs were from 1.1 g/yr in the Second River to 71 g/yr in the upper Passaic River. The load of TEQs (as 2,3,7,8-TCDD) from PCDD/PCDFs and coplanar PCBs in the tributaries were 0.5 mg/yr (milligrams per year) in the Second River, 5.8 mg/yr in the Third River, 130 mg/yr in the Saddle River, 46 mg/yr in the Pompton River, and 100 mg/yr in the upper Passaic River. These loads represent an addition to the TEQ load estimated to cross the head-of-tide of 0.1 percent by the Second River, 0.7 percent by the Third River, and 15 percent by the Saddle River. Loads of sediment-bound trace elements mercury, methylmercury, lead, and cadmium were calculated using concentrations obtained from grab samples, which were assumed to represent average annual concentrations in these rivers. Loads of sediment-bound mercury were estimated to be 1,200 g/yr in the Second River; 130 g/yr in the Third River; 4,200 g/yr in the Saddle River; 3,400 g/yr in the Pompton River; and 6,500 g/yr in the upper Passaic River. Loads of sediment-bound lead were estimated to be 56 kg/yr (kilograms per year) in the Second River; 89 kg/yr in the Third River; 1,140 kg/yr in the Saddle River; 310 kg/yr in the Pompton River; and 1,040 kg/yr in the upper Passaic River. Loads of sediment-bound cadmium were estimated to be 1 kg/yr in the Second River; 0.59 kg/yr in the Third River; 60 kg/yr in the Saddle River; 16 kg/yr in the Pompton River; and 11 kg/yr in the upper Passaic River. These loads indicate the importance of the sediment-bound contributions of organic compounds and trace elements to the upper Passaic and Saddle Rivers. Concentrations of suspended sediment-bound PCBs in the main stem and the West Branch of the Elizabeth River were 806 ng/g and 3,100 ng/g, respectively, representing loads of 40 g/yr and 1,150 g/yr, respectively. These loads were estimated using assumed discharge conditions. Concentrations of suspended sediment-bound PCDD/PCDFs were 7,270 pg/g and 9,980 pg/g in the main stem and West Branch, respectively, representing average annual loads of 0.36 g/yr and 3.7 g/yr, respectively. Total TEQ loads (sum of PCDD/PCDFs and PCBs) were 2.1 mg/yr (as 2,3,7,8-TCDD) in the main stem and 34 mg/yr in the West Branch, respectively. These load estimates, however, were directly related to the assumed annual discharge for the two branches. Long-term measurement of stream discharge and suspended-sediment concentrations would be needed to verify these loads. On the basis of the concentrations measured in this work, it appears that the West Branch is the principal source of PCBs, PCDD/PCDFs, total TEQs, and metals to the main stem of the Elizabeth River. Additional sources of these constituents may exist between the confluence and the head-of-tide.

New Jersey, New York

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