Flow was simulated in 10 aquifers of the New Jersey Coastal Plain using a multilayer finite-difference model for prepumping steady-state conditions and transient conditions from 1896-1981. The highest transmissivity, greater than 10,000 sq ft/day, is in Camden and Gloucester Counties in the Potomac-Raritan-Magothy aquifers; Monmouth and Ocean Counties in the middle aquifer of the Potomac-Raritan Magothy aquifer system; and Ocean, Burlington, Atlantic, and Cape May Counties in the Kirkwood-Cohansey aquifer system. Confining unit leakance is highest, > than 0.001 ft/day/ft in updip areas and lowest, < 0.00001 ft/day/ft, in downdip areas. Areas near the center of the major cones of depression approximate steady-state conditions. However, downdip and offshore areas are under transient conditions. Simulated head changes along the saltwater- freshwater interface boundary indicate that the lower aquifer of the Potomac-Raritan-Magothy aquifer system and the confined Kirkwood aquifer have the greatest potential for updip movement of chlorides. The simulated sources of water to wells in 1978 include: (1) 3% from aquifer storage; (2) 3% from boundary flows; (3) 4% from the ocean and bays; and (4) 90% from streamflow. (USGS)