Geology Reports⌕ Search

Geology topics

Alison D. Gordon

Publications and source records attributed to Alison D. Gordon.

14 recordsLinked to original sources
USGS2026

Computation of regional groundwater budgets for the Virginia Coastal Plain aquifer system

Computation of detailed groundwater flow budgets for subdivisions of the Virginia Coastal Plain aquifer system has enabled quantification and more thorough understanding of groundwater flow within this important water resource. A zone budget analysis based on previously published groundwater models of the Virginia Coastal Plain and Virginia Eastern Shore indicates that groundwater conditions vary substantially throughout the Coastal Plain aquifer system because of local variations in hydrogeology and historical and ongoing variations in groundwater use and management. Decades of substantial groundwater withdrawal from the Coastal Plain aquifer system have altered groundwater flow from predevelopment conditions. Rates of sustainable withdrawal are limited because the downward groundwater flow rate into confined aquifers is a relatively small part of the total groundwater budget for the aquifer system compared to the rate of recharge at the land surface. Analyses of groundwater budgets from the Virginia Coastal Plain model indicate that groundwater flow is generally outward from the surficial aquifer to rivers and coastal waterbodies and downward through a series of underlying aquifers and confining units to the Potomac aquifer, which is the deepest aquifer and the source of most groundwater withdrawals. Downward flow into the Potomac aquifer is estimated to be only 7 percent of total net precipitation-derived net recharge at the land surface but makes up about 66 percent of inflow to the aquifer in Virginia, with much of the remaining inflow occurring laterally from outside of defined groundwater budget regions in Virginia. For several decades prior to 2010, high rates of withdrawal from the Potomac aquifer resulted in substantial decline in groundwater storage in the aquifer and in most overlying aquifers and confining units. From 2010 to 2023, rates of withdrawal substantially lower than the historical maximum resulted in small net increases in groundwater storage in the confined aquifer system for most regions of the Virginia Coastal Plain. Nevertheless, for the same period, groundwater storage for the entire model domain continues to incrementally decline, indicating that storage recovery in Virginia is offset by a continued decrease in storage in areas beneath the Chesapeake Bay or adjacent areas of Maryland and North Carolina. Withdrawals from the Potomac aquifer have induced substantial downward flow which is a large part of groundwater budgets for confined aquifers such as the Potomac. For the most recent simulated conditions (2023) downward groundwater flow continues, but because vertical flow rates are a function of the difference between water pressure in the upper surficial systems and lower confined units, rates of downward flow are lower than those in earlier decades as the confined water levels partially recover from larger groundwater withdrawals in the past. Geographically, groundwater flow is generally inward from perimeter regions of the Virginia Coastal Plain toward central regions with the largest withdrawal rates. Groundwater inflow from coastal regions could be contributing to saltwater intrusion, even though that was not measured in this study. Analyses of groundwater budgets from the Virginia Eastern Shore peninsula, a geographic region of the Virginia Coastal Plain, indicate that groundwater flow for that isolated aquifer system is generally outward from the surficial aquifer to coastal water bodies and downward into the confined Yorktown-Eastover aquifer system, which is the source of most withdrawals. Downward groundwater flow into the confined Yorktown-Eastover aquifer system is estimated to be less than 2 percent of total recharge and less than 9 percent of net recharge at the water table but makes up more than 93 percent of all inflow to the confined aquifer system. Decades of substantial but relatively consistent groundwater withdrawals have induced greater downward flow rates into the confined aquifer system but also have resulted in loss of groundwater from storage. For the most recent simulated period (2023), estimated storage loss accounts for slightly under 7 percent of withdrawals from the confined aquifer system. The reported withdrawal rate for this period from the confined Yorktown-Eastover system is near the highest reported rate for the Virginia Eastern Shore, which means that the storage depletion is expected to continue, even though groundwater levels appear to be relatively stable. Estimated groundwater flow rates upward from the confining unit underlying the Yorktown-Eastover system and low rates of inflow from coastal water bodies underscore ongoing concerns about up-coning and lateral intrusion of salty groundwater.

Virginia↗
USGS2026

Computation of regional groundwater budgets for the Virginia Coastal Plain aquifer system

Computation of detailed groundwater flow budgets for subdivisions of Virginia’s Coastal Plain aquifer system has enabled quantification and more thorough understanding of groundwater flow within this important water resource. A zone budget analysis conducted on previously published groundwater models of the Virginia Coastal Plain and Virginia Eastern Shore shows that groundwater conditions vary substantially throughout the Coastal Plain aquifer system due to local variations in hydrogeology and historical and ongoing variations in groundwater use and management. Decades of substantial groundwater withdrawal from the Coastal Plain aquifer system have fundamentally altered groundwater flow from pre-development conditions. Rates of sustainable withdrawal are limited because the downward groundwater flow rate into confined aquifers supplying groundwater is a relatively small portion of the total groundwater water budget for the aquifer system. Analyses of groundwater budgets from the Virginia Coastal Plain model show that groundwater flow is generally outward from the surficial aquifer to rivers and coastal water bodies and downward through a series of underlying aquifers and confining units to the Potomac aquifer, which is the deepest aquifer and the source of most groundwater withdrawals. Downward flow into the Potomac aquifer currently is estimated to be only 7 percent of total net precipitation-derived net recharge at the land surface but makes up about 66 percent of inflow to the aquifer in Virginia, with much of the remaining inflow occurring laterally from areas outside of defined groundwater budget regions in Virginia. For several decades prior to 2010, high rates of withdrawal from the Potomac aquifer resulted in substantial decline in groundwater storage in the aquifer and in most overlying aquifers and confining units. From 2010 to 2025, rates of withdrawal substantially lower than the historical maximum have resulted in small net increases in groundwater storage in the confined aquifer system for most regions of the Virginia Coastal Plain. Nevertheless, for the same period, groundwater storage for the entire model domain continues to incrementally decline, indicating that storage recovery in Virginia is offset by a continued decrease in storage in areas beneath the Chesapeake Bay or in adjacent areas of Maryland and North Carolina. Withdrawals from the Potomac aquifer have induced substantial downward flow which is a large part of groundwater budgets for confined aquifers such as the Potomac. Downward groundwater flow continues under current conditions, but because vertical flow rates are a function of the difference between water pressure in the upper surficial systems and lower confined units, those rates are lower than those in earlier decades as the confined water levels partially recover from larger groundwater withdrawals in the past. Geographically, groundwater flow is generally inward from perimeter regions of the Virginia Coastal Plain toward central regions with the largest withdrawal rates. Estimated groundwater inflow from coastal regions could be contributing to saltwater intrusion, though that was not measured directly in this study. Analyses of groundwater budgets from the Virginia Eastern Shore peninsula, a geographic region of the Virginia Coastal Plain, show that groundwater flow for that isolated aquifer system is generally outward from the surficial aquifer to coastal water bodies and downward into the confined Yorktown-Eastover aquifer system, which is the source of most withdrawals. Downward groundwater flow into the confined Yorktown-Eastover aquifer system is estimated to be less than 2 percent of total recharge and less than 9 percent of net recharge at the water table but makes up over 93 percent of all inflow to the confined aquifer system. Decades of substantial but relatively consistent groundwater withdrawals have induced greater downward flow rates into the confined aquifer system but also have resulted in loss of groundwater from storage. Currently, estimated storage loss accounts for slightly under 7 percent of withdrawals from the confined aquifer system. The current withdrawal rate from the confined Yorktown-Eastover system is near the highest reported rate for the Eastern Shore, which means that the storage depletion is expected to continue, even though groundwater levels appear to be relatively stable. Estimated groundwater flow rates upward from the confining unit underlying the Yorktown-Eastover system and small rates of inflow from coastal water bodies underscore ongoing concerns about up-coning and lateral intrusion of salty groundwater.

EarthArXiv↗
USGS2023

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↗
USGS2021

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↗
USGS2011

Simulated effects of allocated and projected 2025 withdrawals from the Potomac-Raritan-Magothy aquifer system, Gloucester and Northeastern Salem Counties, New Jersey

Withdrawals from the Potomac-Raritan-Magothy aquifer system in New Jersey, which includes the Upper, Middle, and Lower Potomac-Raritan-Magothy aquifers, are the principal source of groundwater supply in northern Gloucester and northeastern Salem Counties in the New Jersey Coastal Plain. Water levels in these aquifers have declined in response to pumping. With increased population growth and development expected in Gloucester County and parts of Salem County over the next 2 decades (2005-2025), withdrawals from these aquifers also are expected to increase. A steady-state groundwater-flow model, developed to simulate flow in the Potomac-Raritan-Magothy aquifer system in northern Gloucester and northeastern Salem Counties, was calibrated to withdrawal conditions in 1998, when groundwater withdrawals from the Potomac-Raritan-Magothy aquifer system in the model area were more than 10,100 Mgal/yr (million gallons per year). Withdrawals from water-purveyor wells accounted for about 63 percent of these withdrawals, and withdrawals from industrial self-supply wells accounted for about 32 percent. Withdrawals from agricultural-irrigation, commercial self-supply, and domestic self-supply wells accounted for the remaining 5 percent. Results of the 2000 baseline groundwater-flow simulation, incorporating average annual 1999-2001 groundwater withdrawals, indicate that the average simulated water levels in the Upper, Middle, and Lower Potomac-Raritan-Magothy aquifers are 31, 27, and 30 feet below the National Geodetic Vertical Datum of 1929 (NGVD 29), respectively, and the lowest simulated water levels are 77, 65, and 59 feet below NGVD 29, respectively. In the full-allocation scenario, the maximum State-permitted (allocated) groundwater withdrawals totaled 16,567 Mgal/yr, an increase of 72 percent from the 2000 baseline simulation. Results of the full-allocation simulation indicate that the average simulated water levels in the Upper, Middle, and Lower Potomac-Raritan-Magothy aquifers are 49, 43, and 48 feet below NGVD 29, respectively, which are 18, 16, and 18 feet lower, respectively, than in the 2000 baseline simulation. The lowest simulated water levels are 156, 95, and 69 feet below NGVD 29, respectively, which are 79, 30, and 10 feet lower, respectively, than in the 2000 baseline simulation. Simulated net flow from the Potomac-Raritan-Magothy aquifer system to streams is 8,441 Mgal/yr in the 2000 baseline simulation but is 6,018 Mgal/yr in the full-allocation scenario, a decrease of 29 percent from the 2000 baseline simulation. Simulated net flow in the 2000 baseline simulation is 1,183 Mgal/yr from the aquifer system to the Delaware River but in the full-allocation scenario is 1,816 Mgal/yr from the river to the aquifer system. Four other simulations were conducted that incorporated full-allocation conditions at water-purveyor wells in Critical Area 2 but increased or decreased withdrawals at selected water-purveyor wells outside Critical Area 2 and agricultural-irrigation and industrial-self-supply wells in the study area. The results of the four simulations also indicate net flow from the Delaware River to the Potomac-Raritan-Magothy aquifer system. A growth scenario was developed to simulate future withdrawals in 2025 estimated from population projections for municipalities in the Salem-Gloucester study area. Simulated withdrawals for this scenario totaled 10,261 Mgal/yr, an increase of 6 percent from the 2000 baseline simulation. This total includes about 25 Mgal/yr withdrawn from the Englishtown aquifer system for domestic self-supply. This scenario incorporated full-allocation withdrawals at water-purveyor wells in Critical Area 2, and increased withdrawals at water-purveyor wells outside Critical Area 2. Results of this simulation indicate that the average simulated water levels in the Upper, Middle, and Lower Potomac-Raritan-Magothy aquifers are 32, 29, and 32 feet below NGVD 29, respectively, which are 1, 2, and

Scientific Investigations Report↗
USGS2007

Hydrogeology of, and simulation of ground-water flow In, the Pohatcong Valley, Warren County, New Jersey

A numerical ground-water-flow model was constructed to simulate ground-water flow in the Pohatcong Valley, including the area within the U.S. Environmental Protection Agency Pohatcong Valley Ground Water Contamination Site. The area is underlain by glacial till, alluvial sediments, and weathered and competent carbonate bedrock. The northwestern and southeastern valley boundaries are regional-scale thrust faults and ridges underlain by crystalline rocks. The unconsolidated sediments and weathered bedrock form a minor surficial aquifer and the carbonate rocks form a highly transmissive fractured-rock aquifer. Ground-water flow in the carbonate rocks is primarily downvalley towards the Delaware River, but the water discharges through the surficial aquifer to Pohatcong Creek under typical conditions. The hydraulic characteristics of the carbonate-rock aquifer are highly heterogeneous. Horizontal hydraulic conductivities span nearly five orders of magnitude, from 0.5 feet per day (ft/d) to 1,800 ft/d. The maximum transmissivity calculated is 37,000 feet squared per day. The horizontal hydraulic conductivities calculated from aquifer tests using public supply wells open to the Leithsville Formation and Allentown Dolomite are 34 ft/d (effective hydraulic conductivity) and 85 to 190 ft/d (minimum and maximum hydraulic conductivity, respectively, yielding a horizontal anisotropy ratio of 0.46). Stream base-flow data were used to estimate the net gain (or loss) for selected reaches on Brass Castle Creek, Shabbecong Creek, three smaller tributaries to Pohatcong Creek, and for five reaches on Pohatcong Creek. Estimated mean annual base flows for Brass Castle Creek, Pohatcong Creek at New Village, and Pohatcong Creek at Carpentersville (from correlations of partial- and continuous-record stations) are 2.4, 25, and 45 cubic feet per second (ft3/s) (10, 10, and 11 inches per year (in/yr)), respectively. Ground-water ages estimated using sulfur hexafluoride (SF6), chlorofluorocarbon (CFC), and tritium-helium age-dating techniques range from 0 to 27 years, with a median age of 6 years. Land-surface and ground-water water budgets were calculated, yielding an estimated rate of direct recharge tothe surficial aquifer of about 23 in/yr, and an estimated net recharge to the ground-water system within the area underlain by carbonate rock (11.4 mi2) of 29 in/yr (10 in/yr over the entire 33.3 mi2 basin). A finite-difference, numerical model was developed to simulate ground-water flow in the Pohatcong Valley. The four-layer model encompasses the entire carbonate-rock part of the valley. The carbonate-rock aquifer was modeled as horizontally anisotropic, with the direction of maximum transmissivity aligned with the longitudinal axis of the valley. All lateral boundaries are no-flow boundaries. Recharge was applied uniformly to the topmost active layer with additional recharge added near the lateral boundaries to represent infiltration of runoff from adjacent crystalline-rock areas. The model was calibrated to June 2001 water levels in wells completed in the carbonate-rock aquifer, August 2000 stream base-flow measurements, and the approximate ground-water age. The ground-water-flow model was constructed in part to test possible site contamination remediation alternatives. Four previously determined ground-water remediation alternatives (GW1, GW2, GW3, and GW4) were simulated. For GW1, the no-action alternative, simulated pathlines originating in the tetrachloroethene (PCE) and trichloroethene (TCE) source areas within the Ground-Water Contamination Site end at Pohatcong Creek near the confluence with Shabbecong Creek, although some particles went deeper in the aquifer system and ultimately discharge to Pohatcong Creek about 10 miles downvalley in Pohatcong Township. Remediation alternatives GW2, GW3, and GW4 include ground-water withdrawal, treatment, and reinjection. The design for GW2 includes wells in the TCE and PCE source areas that wit

New Jersey↗
USGS2007

Simulated effects of projected 2010 withdrawals on ground-water flow and water levels in the New Jersey coastal plain – A task of the New Jersey Water Supply Plan, 2006 revision

A ground-water flow model previously developed as part of a Regional Aquifer System Analysis (RASA) of the New Jersey Coastal Plain was used to simulate ground-water flow in eight major confined aquifers to help evaluate ground-water resources in support of the New Jersey Department of Environmental Protection's revision of the New Jersey State Water Supply Plan. This model was calibrated to 1998 steady-state and transient conditions. Withdrawals at wells in operation in 1998 were varied in three scenarios to evaluate their effects on flow directions, water levels, and water budgets in the confined aquifers. The scenarios used to predict changes in pumpage from 1998 to 2010 were based on (1) a continuation of 1990-99 trends in water use, (2) public-supply withdrawals estimated from county population projections, and (3) restricted withdrawals in Water-Supply Critical Areas. Total withdrawals in these three scenarios were approximately 366, 362, and 355 million gallons per day, respectively. The results of these simulations are used by New Jersey water-management officials to help address water-supply concerns for the State. In the revision of the New Jersey State Water Supply Plan, the eight major confined aquifers of the New Jersey Coastal Plain and their outcrop areas are divided into 41 hydrologic budget areas (HBAs). Simulation results were used to assess the effects of changing ground-water withdrawals on water levels and the flow budgets in each budget area. Simulation results for each scenario were compared with 1998 (baseline) simulated water levels and flow budgets. The 41 hydrologic budget areas are in areas of large ground-water withdrawals, water-level declines, and (or) saltwater-intrusion potential. Their boundaries are based on various hydrologic, geohydrologic, and withdrawal conditions, such as aquifer extent, location of the 250-milligram-per-liter isochlor, aquifer outcrop area, and ground-water divides. The budget areas include primarily the onshore, freshwater portions of the aquifers. A budget analysis was done for each of the hydrologic budget areas for each scenario. Ground-water withdrawals, leakage to streams, net leakage to overlying and underlying aquifers, lateral flow to adjacent budget areas, and the flow direction at the 250-milligram-per-liter isochlor were evaluated. Although three different methods were applied to predict future pumping rates, the simulated water levels for scenarios 1 and 2 were generally within 2 feet of each other in most areas in the confined aquifers, but differences of more than 2 feet occurred locally. Differences in values of flow-budget components between scenarios 1 and 2 as a percentage change from 1998 values were generally within 2 percent in most hydrologic budget areas, but values of some budget components in some hydrologic budget areas differed by more than 2 percent. Simulated water levels recovered as much as 4 feet more in northeastern Camden and northwestern Burlington Counties in the Lower Potomac-Raritan-Magothy aquifer, and as much as 3 feet more in the same area in the Upper and Middle Potomac-Raritan-Magothy aquifers when pumpage restrictions were imposed in Critical Area 2 (scenario 3). In the Wenonah-Mount-Laurel aquifer, water levels declined continually in Monmouth County (HBA 8) downdip from the outcrop (in Critical Area 1) from 1988 to 2010 in all three scenarios, although most of the water levels farther downdip from this area in Critical Area 1 are still recovering because of mandated reductions in pumpage in the 1990s. In the Englishtown aquifer system, water levels declined continually in small areas in HBA 13 in central Monmouth County (in Critical Area 1) and in western Monmouth County downdip from the outcrop from 1988 to 2010 in all three scenarios, although most of the water levels farther downdip from this area are still recovering because of the mandated reductions in pumpage. In the Upper Potomac-Raritan-Magothy aquif

New Jersey↗
USGS2003

Simulation of the Ground-Water Flow System in 1992, and Simulated Effects of Projected Ground-Water Withdrawals in 2020 in the New Jersey Coastal Plain

In 1992, ground-water withdrawals from the unconfined and confined aquifers in the New Jersey Coastal Plain totaled about 300 million gallons per day, and about 70 percent (200 million galllons per day) of this water was pumped from confined aquifers. The withdrawals have created large cones of depression in several Coastal Plain aquifers near populated areas, particularly in Camden and Ocean Counties. The continued decline of water levels in confined aquifers could cause saltwater intrusion, reduction of stream discharge near the outcrop areas of these aquifers, and depletion of the ground-water supply. Because of this, withdrawals from wells located within these critical areas have been reduced in the Potomac-Raritan-Magothy aquifer system, the Englishtown aquifer system, and the Wenonah-Mount Laurel aquifer. A computer-based model that simulates freshwater and saltwater flow was used to simulate transient ground-water flow conditions and the location of the freshwater-saltwater interface during 1989-92 in the New Jersey Coastal Plain. This simulation was used as the baseline for comparison of water levels and flow budgets. Four hypothetical withdrawal scenarios were simulated in which ground-water withdrawals were either increased or decreased. In scenario 1, withdrawals from wells located within critical area 2 in the Potomac-Raritan-Magothy aquifer system were reduced by amounts ranging from 0 to 35 percent of withdrawals prior to 1992. Critical area 2 is mainly located in Camden County, and most of Burlington and Gloucester Counties. With the reductions, water levels recovered about 30 feet in the regional cone of depression centered in Camden County in the Upper Potomac-Raritan-Magothy aquifer and by 20 ft in the Lower and Middle Potomac-Raritan-Magothy aquifers. In scenarios 2 to 4, withdrawals projected for 2020 were input to the model. In scenario 2, withdrawal restrictions within the critical areas were imposed in the Potomac-Raritan-Magothy aquifer system, the Englishtown aquifer system, and the Wenonah-Mount Laurel aquifer, but withdrawals were increased outside the critical areas to the projected 2020 demand. With withdrawals restrictions in critical areas, water levels recovered about 20 feet at the center of the regional cone of depression in the Upper Potomac-Raritan Magothy aquifer. Water levels recovered by about 20 feet at the center of a regional cone of depression in the Englishtown aquifer system in Ocean County, and by about 20 feet in the Wenonah-Mount Laurel aquifer in the same area. In scenario 3, withdrawals were increased to the projected 2020 demand inside and outside the critical areas. As a result, water levels declined as much as 20 feet at the center of a regional cone of depression in the Englishtown aquifer system in Ocean County, and as much as 10 feet in the Wenonah-Mounty Laurel aquifer near this area. The Englishtown aquifer system and the Wenonah-Mount Laurel aquifer are particularly sensitive to increases and decreases in withdrawals because in certain areas the transmissivities of these aquifers are lower than the transmissivities of other confined aquifers of the New Jersey Coastal Plain, and because these aquifers are hydraulically connected. Simulated water levels declined by as much as 10 ft at the center of the regional cone of depression in Atlantic County. In scenario 4, withdrawal amounts were equal to that in scenario 2, except an additional 13.2 million gallons per day was withdrawn from hypothetical wells located outside the critical areas in the Upper Potomac-Raritan-Magothy aquifer, Englishtown aquifer system, and the Wenonah-Mount Laurel aquifer. The additional withdrawals resulted in increased leakage from overlying aquifers to the Wenonah-Mount Laurel aquifer and subsequently to the Englishtown aquifer system.

Water-Resources Investigations Report↗
USGS2002

Simulation of Transient Ground-Water Flow in the Valley-Fill Aquifers of the Upper Rockaway River Basin, Morris County, New Jersey

More than 90 percent of the public water supply in the upper Rockaway River Valley in Morris County, New Jersey, is obtained from ground-water withdrawals from the valley-fill aquifers. During 1997, an average of 9.6 million gallons per day of ground water was withdrawn from these aquifers. The aquifer system consists of an unconfined aquifer (upper aquifer) and a locally confined aquifer (lower aquifer), which are composed of sands and gravels. These aquifers are separated by a discontinuous confining unit that consists mostly of silt and clay. Increases in ground-water withdrawals can induce movement of water from streams to wells, increase flow from the upper aquifer to the lower aquifer, and reduce base flow in the Rockaway River downstream. A ground-water-flow model was used to simulate and quantify the effects of current withdrawals on the valley-fill aquifer system under transient monthly conditions. Recharge over the model area varies both spatially and temporally. Part of model calibration consisted of adjusting percentages of monthly precipitation that recharges the valley-fill aquifer system. More recharge occurs during winter and spring than during summer and fall. This seasonal variation affects ground-water discharge to the Rockaway River. Ground-water withdrawals from the valleyfill aquifers also affect ground-water discharge to the Rockaway River. Three scenarios were simulated to observe the effects of ground-water withdrawals on ground-water discharge to the Rockaway River and to determine the extent to which variations in rates of withdrawals correspond to variations in rates of streamflow depletion. Streamflow depletion was estimated by comparing model-computed ground-water discharge for the three scenarios with the modelcomputed ground-water discharge under transient conditions. In scenario 1, all pumpage was removed from the model. In scenarios 2 and 3, 1 million gallons per day of ground-water withdrawals in excess of the current pumpage was withdrawn from the valley-fill aquifers. In scenario 2, the additional 1 million gallons per day of withdrawals were made from a hypothetical well located in the upper aquifer about 250 feet from the river. In scenario 3, the additional withdrawals were made from a hypothetical well located in the lower aquifer about 1,750 feet from the river. Results of scenario 1 indicated that the difference between the streamflow depletion and withdrawals is small; increases in ground-water withdrawals from the valley-fill aquifers correspond to decreases in ground-water discharge to the Rockaway River of approximately the same amount. Results of scenario 2 and 3 indicated that a lag time could occur between the introduction of withdrawals and the full magnitude of the effects of the withdrawals on streamflow depletion. A lag time of about seven months occurred for scenario 2 with the well placed in the upper aquifer. A longer lag time of more than 1.5 years occurred with the well placed in the lower aquifer and separated from the upper aquifer by a confining unit (scenario 3). Extreme low flow in the Rockaway River is mostly base flow. A flow-duration analysis of the Rockaway River at the surface-water gaging station upstream from the Boonton Reservoir during the drought of 1961-66 indicated that streamflow from the upper Rockaway River Basin alone might not be sufficient to meet the minimum passing flow of 7 million gallons per day during a drought. Under similar drought conditions today, during 3.2 percent of the drought time, streamflow at this station upstream from the reservoir would be less than the minimum passing flow requirement downstream from the reservoir.

Water-Resources Investigations Report↗
USGS1999

Simulation of ground-water flow and movement of the freshwater-saltwater interface in the New Jersey coastal plain

The confined aquifers of the New Jersey Coastal Plain are sands that range in thickness from 50 to 600 feet and are separated by confining units. The confining units are composed of silts and clays that range in thickness from 500 to 1,000 feet. The aquifers are recharged by precipitation on their outcrop areas. This water then flows laterally downdip and vertically to the deeper confined aquifers. The confined aquifers ultimately discharge to the Raritan and Delaware Bays and to the Atlantic Ocean. In 1988, ground-water withdrawals from confined and unconfined New Jersey Coastal Plain aquifers were approximately 345 million gallons per day, more than 75 percent of which was pumped from the confined aquifers. These withdrawals have created large cones of depression in several Coastal Plain aquifers near populated areas, particularly in Camden and Monmouth Counties. The continued decline of water levels in confined aquifers can cause saltwater intrusion, reduce stream discharge near the outcrop areas, and threaten the quality of the ground-water supply. SHARP, a quasi-three-dimensional finite-difference computer model that can simulate freshwater and saltwater flow, was used to simulate the ground-water flow system in the New Jersey Coastal Plain, including the location and movement of the freshwater-saltwater interface in nine aquifers and eight intervening confining units. The freshwater-saltwater interface is defined as the hypothetical line seaward of which the chloride concentration is equal to or greater than 10,000 milligrams per liter. Model simulations were used to estimate the location and movement of the freshwater-saltwater interface resulting from (1) eustatic sea-level changes over the past 84,000 years, (2) ground-water withdrawals from 1896 through 1988, (3) and future ground-water withdrawals from 1988 to 2040 from Coastal Plain aquifers. Simultion results showed that the location and movement of the freshwater-saltwater interface are more dependent on the historical sea level than on the stresses imposed on the flow system by ground-water withdrawals from the Coastal Plain aquifers from 1896 to 1988. Results of a predictive simulation in which pumpage from existing wells was increased by 30 percent indicate that additional withdrawals from each of the eight confined aquifers in the Coastal Plain would broaden and deepen the existing cones of depression and result in significant drawdowns from the 1988 potentiometric surfaces. Drawdowns of 30 feet were simulated at the center of the cone of depression in the Upper, Middle, and Lower Potomac-Raritan-Magothy aquifers in Camden and Ocean Counties. Simulated drawdowns exceeded 80 feet at the center of the cone of depression in the Wenonah-Mount Laurel and Englishtown aquifers in Monmouth County. Drawdowns of 30 feet were simulated in the lower Kirkwood-Cohansey and confined Kirkwood aquifers in Cape May County. Simulation results showed that the increase in ground-water withdrawals would result in only minimal movement of the freshwater-saltwater interface by 2040, despite large drawdowns.

New Jersey↗
USGS1993

Hydrogeology of, and simulated ground-water flow in, the valley-fill aquifers of the upper Rockaway River basin, Morris County, New Jersey

Public water supply in the Rockaway River valley depends almost entirely on ground water from wells in the valley-fill deposits. Ground-water withdrawals from these deposits in 1986 were about 9.1 million gallons per day. A steady-state ground-water flow model was developed to quantify the effects of ground-water withdrawals on water levels in the valley-fill aquifers and on ground-water discharge to the Rockaway River. The ground-water-flow model, which represents an aquifer system consisting of an unconfined and a confined aquifer separated by a discontinuous confining unit, was implemented to examine aquifer resonse to current and predicted ground-water withdrawals in areas of proposed well sites and the effect of increased ground-water withdrawals on ground-water discharge to the river. Ground-water flow to wells in the valley-fill aquifers is sustained by increased vertical flow between the two aquifers, the diversion of ground water that had discharged to the Rockaway River, and induced seepage resulting from pumping near the Rickaway River. If the rate of ground-water recharge decreases or if the rate of ground-water withdrawals from the valley-fill aquifers increases, ground-water discharge to the Rockaway River above the Boonton Reservoir will decrease by an equivalent amount. The average annual base flow of the Rockaway River above the Boonton Reservoir will meet the minimum passing flow requirement of 7 million gallons per day under conditions of average annual ground-water recharge, and increased ground-water withdrawals anticipated by the years 2000 and 2040. For anticipated increases in withdrawals to 11.5 million gallons per day by the the year 2000, and 14.6 million gallons per day by the year 2040, base flow to the Rockaway River above the Boonton Reservoir may not be sufficient to meet the minimum required reservoir outflow during extended periods of decreased recharge, such as drought, as much as 5 and 11 percent of the time, respectively.

Water-Resources Investigations Report↗
USGS1989

Water-quality data for the Potomac-Raritan-Magothy aquifer system in the northern coast plain of New Jersey, 1923-86

Ground-water-quality data for the upper and middle aquifers of the Potomac-Raritan-Magothy aquifer system in Middlesex and Monmouth Counties are compiled for the period 1923-86. A total of 330 wells were sampled: 192 wells in the upper aquifer and 138 wells in the middle aquifer. Most of the complete water-quality analyses were collected after September 1984, as part of a regional ground-water assessment. Well-construction data for the sampled wells also are presented. Public-supply, domestic-supply, industrial, commercial, irrigation, and observation wells were sampled for the study. Field measurements made at the time of sample collection include water temperature, specific conductance , dissolved oxygen, pH, alkalinity, and bicarbonate concentration. Laboratory determinations include common ions, silica, dissolved solids, trace metals, volatile organic compounds, and pesticides. A quality-assurance program was followed to evaluate and assure the quality of the data. The report also contains a table of lithologic and hydrologic characteristics of the geologic units in the study area, a table of chloride concentrations and field measurements from 1923-86, and statistical summaries of selected water-quality data for the upper and middle aquifers. Many constituents were found in a wide range of concentrations. Water from more than 25 percent of the wells sampled contained lead concentrations above the detection limit of 10 ug/L (micrograms per liter). Included in this number are some wells that had lead concentrations greater than the U.S. Environmental Protection Agency (USEPA) primary drinking-water regulation of 50 ug/L. Cadmium concentrations, although lower than lead concentrations, followed a similar pattern. Water from approximately 25 percent of the wells in the upper aquifer, contain cadmium concentrations equal to or greater than the detection limit of 1 ug/L. Dissolved iron concentrations ranged from 5 ug/L to 480,000 ug/L. Water from more than 50 percent of the wells sampled contained iron concentrations in excess of the USEPA secondary drinking-water recommended limit of 300 ug/L. Chloride concentrations greater than the USEPA secondary drinking-water recommended limit of 250 milligrams per liter were found in samples from wells located in the cities of Perth Amboy and South Amboy; in the boroughs of Keansburg, Sayreville, Keyport, and Union Beach; and in the townships of Old Bridge and Woodbridge. Of 21 samples collected from wells screened in the upper aquifer and analyzed for 30 volatile organic compounds (VOCs), 5 samples contained at least 1 VOC at or above the detection limit. In the middle aquifer, 12 of the 21 samples collected and analyzed for VOCs contained at least 1 VOC greater than the detection limit. Concentrations of pesticides generally were low. Of the 43 samples collected from wells screened in the upper aquifer and analyzed for pesticides, 4 samples contained concentrations of pesticides at or greater than the detection limit. In the middle aquifer, 6 of 38 samples collected and analyzed for 32 pesticides had at least 1 pesticide with a concentration greater than the detection limit.

New Jersey↗