Science Center Objects

Simulation of Groundwater Flow in the Potomac-Raritan-Magothy Aquifer System, Pennsauken Township and Vicinity

By Daryll A. Pope and Martha K. Watt


MODEL VERSION/TYPE: MODFLOW-2000, steady state and transient

AREA STUDIED: Northwest Camden County

AQUIFERS SIMULATED: Potomac-Raritan-Magothy aquifer system


MODEL SIZE: 11 layers, 85 rows, 108 columns

MINIMUM GRID SPACING:206 feet x 219 feet

MODEL ARCHIVE is available by email request at


The Potomac-Raritan-Magothy aquifer system is one of the primary sources of potable water in the Coastal Plain of New Jersey, particularly in heavily developed areas along the Delaware River. In Pennsauken Township, Camden County, local drinking-water supplies from this aquifer system have been contaminated by hexavalent chromium at concentrations that exceed the New Jersey maximum contaminant level. In particular, groundwater at the Puchack well field has been adversely affected to the point where, since 1984, water is no longer withdrawn from this well field for public supply. The area that contains the Puchack well field was added to the National Priorities List in 1998 as a Superfund site.

The U.S. Geological Survey (USGS) conducted a reconnaissance study from 1996 to 1998 during which hydrogeologic and water-quality data were collected and a groundwater-flow model was developed to describe the conditions in the aquifer system in the Pennsauken Township area. The current investigation by the USGS, in cooperation with the U.S. Environmental Protection Agency (USEPA), is an extension of the previous study. Results of the current study can be applied to a Remedial Investigation and Feasibility Study conducted at the Puchack well field Superfund site.

The USGS study collected additional data on the hydrogeology and water-quality in the area. These data were incorporated into a refined model of the groundwater-flow system in the Potomac-Raritan-Magothy aquifer system. A finite-difference model was developed to simulate groundwater flow and the advective transport of chromium-contaminated groundwater in the aquifers of the Potomac-Raritan-Magothy aquifer system in the Pennsauken Township area. An 11-layer model was used to represent the complex hydrogeologic framework. The model was calibrated using steady-state water-level data from March 1998, April 1998, and April 2001. Water-level recovery during the shutdown of Puchack 1 during March to April 1998 was simulated to evaluate model performance in relation to changing stresses. The Delaware River contributes appreciable flow to the groundwater system from areas where the Middle and Lower aquifers crop out beneath the river. A transient simulation of an aquifer test near the Delaware River was run to help characterize the hydraulic conductivity of the riverbed sediments represented in the model. Vertical flow across confining units between the aquifers is highly variable and is important in the movement of water and associated contaminants through the flow system. The model was embedded within a regional model of the Potomac-Raritan-Magothy aquifer system in Camden County.

In general, a simulation of baseline conditions, which can provide a representation on which simulations of various alternatives can be based for the feasibility study, incorporated average conditions from 1998 to 2000. Groundwater withdrawals within the model area during this period averaged about 14 Mgal/d. Regional groundwater flow is from recharge areas and from the Delaware River to downgradient pumped wells located just east of the model area in central Camden County. Simulation results show an important connection between the Intermediate sand and the Lower aquifer of the Potomac-Raritan-Magothy aquifer system in the vicinity of the chromium-contaminated area. The Delaware River contributes nearly 10 Mgal/d to the flow system, whereas recharge contributes about 6 Mgal/d. Groundwater withdrawals within the model area account for nearly 14 Mgal/d (mostly from the Lower aquifer of the Potomac-Raritan-Magothy aquifer system).