USGS ScienceSearch

USGS · 70022871

Fractured-aquifer hydrogeology from geophysical logs: Brunswick group and Lockatong Formation, Pennsylvania

Abstract

The Brunswick Group and the underlying Lockatong Formation are composed of lithified Mesozoic sediments that constitute part of the Newark Basin in southeastern Pennsylvania. These fractured rocks form an important regional aquifer that consists of gradational sequences of shale, siltstone, and sandstone, with fluid transport occurring primarily in fractures. An extensive suite of geophysical logs was obtained in seven wells located at the borough of Lansdale, Pennsylvania, in order to better characterize the areal hydrogeologic system and provide guidelines for the refinement of numerical ground water models. Six of the seven wells are approximately 120 m deep and the seventh extends to a depth of 335 m. Temperature, fluid conductivity, and flowmeter logs are used to locate zones of fluid exchange and to quantify transmissivities. Electrical resistivity and natural gamma logs together yield detailed stratigraphic information, and digital acoustic televiewer data provide magnetically oriented images of the borehole wall from which almost 900 fractures are identified. Analyses of the geophysical data indicate that the aquifer penetrated by the deep well can be separated into two distinct structural domains, which may, in turn, reflect different mechanical responses to basin extension by different sedimentary units: 1. In the shallow zone (above 125 m), the dominant fracture population consists of gently dipping bedding plane partings that strike N46°E and dip to the northwest at about 11 degrees. Fluid flow is concentrated in the upper 80 m along these subhorizontal fractures, with transmissivities rapidly diminishing in magnitude with depth. 2. The zone below 125 m marks the appearance of numerous high-angle fractures that are orthogonal to the bedding planes, striking parallel but dipping steeply southeast at 77 degrees. This secondary set of fractures is associated with a fairly thick (approximately 60 m) high-resistivity, low-transmissivity sandstone unit that is abruptly terminated by a thin shale bed at a depth of 190 m. This lower contact effectively delineates the aquifer's vertical extent at this location because no detectable evidence of ground water movement is found below it. Thus, fluid flow is controlled by fractures, but fracture type and orientation are related to lithology. Finally, a transient thermal-conduction model is successfully applied to simulate observed temperature logs, thereby confirming the effects of ground-surface warming that occurred in the area as a result of urbanization at the turn of the century. The systematic warming of the upper 120 m has increased the transmissivity of this aquifer by almost 10%, simply due to changes in fluid viscosity and density.

Explore related subjects

90° N90° S · 180° W ← longitude → 180° E
Source-reported bounding extent: 40.20833° to 40.25833° latitude; -75.31667° to -75.25° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Roger H. Morin, Lisa A. Senior, Edward R. Decker. 2005-08-04. Fractured-aquifer hydrogeology from geophysical logs: Brunswick group and Lockatong Formation, Pennsylvania. https://doi.org/10.1111/j.1745-6584.2000.tb00329.x

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

Evaluating micrometeorological estimates of groundwater discharge from Great Basin desert playas

Groundwater availability studies in the arid southwestern United States traditionally have assumed that groundwater discharge by evapotranspiration (ETg) from desert playas is a significant component of the groundwater budget. However, desert playa ETg rates are poorly constrained by Bowen Ratio energy budget (BREB) and eddy-covariance (EC) micrometeorological measurement approaches. Best attempts by previous studies to constrain ETg from desert playas have resulted in ETg rates that are within the measurement error of micrometeorological approaches. This study uses numerical models to further constrain desert playa ETg rates that are within the measurement error of BREB and EC approaches, and to evaluate the effect of hydraulic properties and salinity-based groundwater-density contrasts on desert playa ETg rates. Numerical models simulated ETg rates from desert playas in Death Valley, California and Dixie Valley, Nevada. Results indicate that actual ETg rates from desert playas are significantly below the uncertainty thresholds of BREB- and EC-based micrometeorological measurements. Discharge from desert playas likely contributes less than 2 percent of total groundwater discharge from Dixie and Death Valleys, which suggests discharge from desert playas also is negligible in other basins. Simulation results also show that ETg from desert playas primarily is limited by differences in hydraulic properties between alluvial fan and playa sediments and, to a lesser extent, by salinity-based groundwater density contrasts.

Ground Water

A semi-structured MODFLOW-USG model to evaluate local water sources to wells for decision support

In order to better represent the configuration of the stream network and simulate local groundwater-surface water interactions, a version of MODFLOW with refined spacing in the topmost layer was applied to a Lake Michigan Basin (LMB) regional groundwater-flow model developed by the U.S. Geological. Regional MODFLOW models commonly use coarse grids over large areas; this coarse spacing precludes model application to local management issues (e.g., surface-water depletion by wells) without recourse to labor-intensive inset models. Implementation of an unstructured formulation within the MODFLOW framework (MODFLOW-USG) allows application of regional models to address local problems. A “semi-structured” approach (uniform lateral spacing within layers, different lateral spacing among layers) was tested using the LMB regional model. The parent 20-layer model with uniform 5000-foot (1524-m) lateral spacing was converted to 4 layers with 500-foot (152-m) spacing in the top glacial (Quaternary) layer, where surface water features are located, overlying coarser resolution layers representing deeper deposits. This semi-structured version of the LMB model reproduces regional flow conditions, whereas the finer resolution in the top layer improves the accuracy of the simulated response of surface water to shallow wells. One application of the semi-structured LMB model is to provide statistical measures of the correlation between modeled inputs and the simulated amount of water that wells derive from local surface water. The relations identified in this paper serve as the basis for metamodels to predict (with uncertainty) surface-water depletion in response to shallow pumping within and potentially beyond the modeled area, see Fienen et al. (2015a).

Illinois, Indiana, Michigan, Wisconsin

Extending the MODPATH algorithm to rectangular unstructured grids

The recent release of MODFLOW-USG, which allows model grids to have irregular, unstructured connections, requires a modification of the particle-tracking algorithm used by MODPATH. This paper describes a modification of the semi-analytical particle-tracking algorithm used by MODPATH that allows it to be extended to rectangular-based unstructured grids by dividing grid cells with multi-cell face connections into sub-cells. The new method will be incorporated in the next version of MODPATH which is currently under development.

Ground Water