Hydrogeologic Framework of Fractured Sedimentary Rock,Newark Basin,New Jersey

The hydrogeologic framework of fractured sedimentary bedrock at the former Naval Air Warfare Center (NAWC), ¹ 1 For more information on the USGS Toxics Substances Hydrology Program at the Naval Air Warfare Center visit the NAWC website at http://nj.usgs.gov/nawc/
Trenton, New Jersey, a trichloroethylene (TCE)-contaminated site in the Newark Basin, is developed using an understanding of the geologic history of the strata, gamma-ray logs, and rock cores. NAWC is the newest field research site established as part of the U.S. Geological Survey Toxic Substances Hydrology Program, Department of Defense (DoD) Strategic Environmental Research and Development Program, and DoD Environmental Security Technology Certification Program to investigate contaminant remediation in fractured rock. Sedimentary bedrock at the NAWC research site comprises the Skunk Hollow, Byram, and Ewing Creek Members of the Lockatong Formation and Raven Rock Member of the Stockton Formation. Muds of the Lockatong Formation that were deposited in Van Houten cycles during the Triassic have lithified to form the bedrock that is typical of much of the Newark Basin. Four lithotypes formed from the sediments include black, carbon-rich laminated mudstone, dark-gray laminated mudstone, light-gray massive mudstone, and red massive mudstone. Diagenesis, tectonic compression, off-loading, and weathering have altered the rocks to give some strata greater hydraulic conductivity than other strata. Each stratum in the Lockatong Formation is 0.3 to 8 m thick, strikes N65°E, and dips 25° to 70°NW. The black, carbon-rich laminated mudstone tends to fracture easily, has a relatively high hydraulic conductivity and is associated with high natural gamma-ray count rates. The dark-gray laminated mudstone is less fractured and has a lower hydraulic conductivity than the black carbon-rich laminated mudstone. The light-gray and the red massive mudstones are highly indurated and tend to have the least fractures and a low hydraulic conductivity. The differences in gamma-ray count rates for different mudstones allow gamma-ray logs to be used to correlate and delineate the lithostratigraphy from multiple wells. Gamma-ray logs and rock cores were correlated to develop a 13-layer gamma-ray stratigraphy and 41-layer lithostratigraphy throughout the fractured sedimentary rock research site. Detailed hydrogeologic framework shows that black carbon-rich laminated mudstones are the most hydraulically conductive. Water-quality and aquifer-test data indicate that groundwater flow is greatest and TCE contamination is highest in the black, carbon- and clay-rich laminated mudstones. Large-scale groundwater flow at the NAWC research site can be modeled as highly anisotropic with the highest component of permeability occurring along bedding planes.     

Magnetic Anisotropy as an Aid to Identifying CRM and DRM in Red Sedimentary Rocks

To further evaluate the potential of magnetic anisotropy techniques for determining the origin of the natural remanent magnetization (NRM) in sedimentary rocks, several new remanence anisotropy measurement techniques were explored. An accurate separation of the remanence anisotropy of magnetite and hematite in the same sedimentary rock sample was the goal.In one technique, Tertiary red and grey sedimentary rock samples from the Orera section (Spain) were exposed to 13 T fields in 9 different orientations. In each orientation, alternating field (af) demagnetization was used to separate the magnetite and hematite contributions of the high field isothermal remanent magnetization (IRM). Tensor subtraction was used to calculate the magnetite and hematite anisotropy tensors. Geologically interpretable fabrics did not result, probably because of the presence of goethite which contributes to the IRM. In the second technique, also applied to samples from Orera, an anisotropy of anhysteretic remanence (AAR) was applied in af fields up to 240 mT to directly measure the fabric of the magnetite in the sample. IRMs applied in 2 T fields followed by 240 mT af demagnetization, and thermal demagnetization at 90°C to remove the goethite contribution, were used to independently measure the hematite fabric in the same samples. This approach gave geologically interpretable results with minimum principal axes perpendicular to bedding, suggesting that the hematite and magnetite grains in the Orera samples both carry a depositional remanent magnetization (DRM). In a third experiment, IRMs applied in 13 T fields were used to measure the magnetic fabric of samples from the Dome de Barrot area (France). These samples had been demonstrated to have hematite as their only magnetic mineral. The fabrics that resulted were geologically interpretable, showing a strong NW-SE horizontal lineation consistent with AMS fabrics measured in the same samples. These fabrics suggest that the rock's remanence may have been affected by strain and could have originated as a DRM or a CRM.Our work shows that it is important to account for the presence of goethite when using high field IRMs to measure the remanence anisotropy of hematite-bearing sedimentary rocks. It also shows that very high magnetic fields (>10 T) may be used to measure the magnetic fabric of sedimentary rocks with highly coercive magnetic minerals without complete demagnetization between each position, provided that the field magnetically saturates the rock.     

Water Well Hydrographs: An Underutilized Resource for Characterizing Subsurface Conditions

Many of the world's major aquifers are under severe stress as a result of intensive pumping to support irrigated agriculture and provide drinking water supplies for millions. The question of what the future holds for these aquifers is one of global importance. Without better information about subsurface conditions, it will be difficult to reliably assess an aquifer's response to management actions and climatic stresses. One important but underutilized source of information is the data from monitoring well networks that provide near-continuous records of water levels through time. Most organizations running these networks are, by necessity, primarily focused on network maintenance. The result is that relatively little attention is given to interpretation of the acquired hydrographs. However, embedded in those hydrographs is valuable information about subsurface conditions and aquifer responses to natural and anthropogenic stresses. We demonstrate the range of insights that can be gleaned from such hydrographs using data from the High Plains aquifer index well network of the Kansas Geological Survey. We show how information about an aquifer's hydraulic state and lateral extent, the nature of recharge, the hydraulic connection to the aquifer and nearby pumping wells, and the expected response to conservation-based pumping reductions can be extracted from these hydrographs. The value of this information is dependent on accurate water-level measurements; errors in those measurements can make it difficult to fully exploit the insights that water-well hydrographs can provide. We therefore conclude by presenting measures that can help reduce the potential for such errors.     

An Evaluation of Soil-Gas Surveying for H2S for Locating Subsurface Hydrocarbon Contamination

A soil-gas survey was conducted at a gasoline service station and a former fire training facility to determine if surveying for hydrogen sulfide could be useful in locating nonaqueous phase hydrocarbon fuel in the subsurface. Relative; to total organic vapor, oxygen, and carbon dioxide distributions, detectable hydrogen sulfide concentrations were much more restricted to the suspected source vicinity at both sites. Appreciable levels of soil-gas hydrogen sulfide. up to 600 Vppb. were observed in areas characterized by anaerobic or microaerophilic conditions having bulk oxygen levels below 4 percent. Based on the hydrogen sulfide distribution, nonaqueous phase hydrocarbon fuel was located at each site. These results suggest that soil-gas surveying for hydrogen sulfide may help locale mobile or residual gasoline and other nonaqueous phase hydrocarbons in the subsurface.     

An Evaluation of Interpolation Methodologies for Generating Three-Dimensional Hydraulic Property Distributions from Measured Data

The process of attempting to model ground-water systems requires a good understanding of the spatial variation of aquifer hydraulic properties. The capabilities of the more recent innovative flowmeters such as the electromagnetic and heat pulse flowmeters provide the sensitivity to measure ambient flows and pump-induced flows. These flowmeters provide the measurements of pump-induced vertical flows which are analyzed to obtain vertical variations in horizontal hydraulic conductivity, K(z). With discrete areal K-values, K(x, y), and vertical profiles of K, provided by multiwell testing, the essential elements are present to produce a three-dimensional hydraulic conductivity field. The advent of these new flow measuring devices has contributed much to the motivation behind this paper. This paper presents the results of applying deterministic and stochastic methodology to the three-dimensional interpolation of hydraulic properties, specifically, hydraulic conductivity, K. Three of the approaches applied in this paper are deterministic in nature, inverse-distance weighting, inverse-distance-squared weighting, and ordinary kriging, while the fourth is a stochastic approach based on self-affine fractals. All of the methods are applied to measured data collected from 14 wells at a site in the United States near Mobile, Alabama. The three-dimensional K-distributions generated by each of the methods are used as inputs to an advective based transport model with the resulting model output compared to a two-well tracer study run previously at the same site.     

Anisotropy rock physics model for the Longmaxi shale gas reservoir,Sichuan Basin,China

The preferred orientation of clay minerals dominates the intrinsic anisotropy of shale. We introduce the clay lamination (CL) parameter to the Backus averaging method to describe the intrinsic shale anisotropy induced by the alignment of clay minerals. Then, we perform the inversion of CL and the Thomsen anisotropy parameters. The direct measurement of anisotropy is difficult because of the inability to measure the acoustic velocity in the vertical direction in boreholes and instrument limitations. By introducing the parameter CL, the inversion method provides reasonable estimates of the elastic anisotropy in the Longmaxi shale. The clay content is weakly correlated with the CL parameter. Moreover, the parameter CL is abnormally high at the bottom of the Longmaxi and Wufeng Formations, which are the target reservoirs. Finally, we construct rock physics templates to interpret well logging and reservoir properties.     

Procedures for Sampling Deep Subsurface Microbial Communities in Unconsolidated Sediments

A corehole sampling project utilizing a wireline coring system provided sediment samples for microbiological characterization from deep unconsolidated sediments. Sampling tools were developed or modified to minimize contamination during sample acquisition and to facilitate stringent decontamination requirements. Quality assurance procedures, including the use of tracers, were implemented to minimize and quantify contamination from drilling hardware, drilling fluids and sample processing. Tracers included microspheres, potassium bromide, rhodamine dye, and perfluorocarbons, which enabled the detection and measurement of 1mg of drilling fluid per kg of sediment. In addition, sample processing was performed on-site in an anaerobic chamber to prevent exposure of the subsurface materials to atmospheric oxygen concentrations. Sediment samples were then disbursed to investigators at National Laboratories and universities funded through the Department of Energy Subsurface Science Microbiology Program for microbiological characterization. Results of these efforts demonstrated that representative subsurface samples were collected and disbursed.     

A Genetic Algorithm for Stochastic Inversion in Contaminant Subsurface Hydrology

Identifying the spatial distribution of hydrological properties of aquifers is a key problem in subsurface hydrology. The aquifer structure plays an important role in contaminant transport. Identifying the properties (primarily the hydraulic conductivity) is essentially an inversion problem that is ill-posed, non-unique and computationally intensive by definition. In this work, the non-uniqueness of the inverse problem is tackled via a novel Genetic Algorithm approach combined with a geostatistical method (Sequential Indicator Simulations) for construction of realizations of properties spatial distributions, which are modeled as random. The Genetic Algorithm cross-over operator is based on a novel concept of pilot-planes: daughter realizations adopt pilot-planes from one of their parents. In addition, each aquifer realization is conditioned on the geological hard data and is constructed by sampling the facies distribution, evaluated by indicator variograms. The approach is illustrated in two test cases: a synthetic two-dimensional (2D) case and an actual three-dimensional (3D) case. The results have shown the ability of the proposed approach to generate a set of realizations, where each individual exhibits minor deviations from the measurements. Further, a comparison between the proposed approach and direct (Monte Carlo) approach shows that the Genetic Algorithm was able to generate an ensemble of solutions with a better fitting of the measurements than the direct approach by a significantly reduced computational effort.     

A New Automatic Subsurface Gas Monitoring System for Seismogeochemical Studies, Installed in Haruno Borehole, Shizuoka Prefecture, Central Japan

The results of subsurface gas monitoring by application of gas chromatography (GC) to the gas composition of bubbles associated with groundwater for seismogeochemical studies are reported. An automated gas monitoring system was used to determine gas compositions in a 500-m borehole at the Haruno Crustal Movement Observation Site (HOS), central Japan during period 1, from December 1999 to December 2000. The average ± two standard deviation (2SD) compositions of gases in this period were He = 82 ± 29 ppmV, H₂ = 170 ± 62 ppmV, Ar = 0.05 ± 0.07%, N₂ = 50 ± 8%, and CH₄ = 45 ± 6%. A new automated gas monitoring system equipped with a micro-GC was installed in the borehole at the HOS, and gas bubbles from the borehole were monitored during period 2, from December 2006 to March 2007. The average ± two standard deviation (2SD) compositions of gases in this period were He = 8 ± 7 ppmV, H₂ = 13 ± 15 ppmV, Ar = 0.6 ± 0.3%, N₂ = 66 ± 7%, and CH₄ = 14 ± 14%. The gas concentration ratios (He/Ar, H₂/Ar, N₂/Ar, and CH₄/Ar) fluctuated significantly over time and repeatedly showed abrupt spike-like increases during period 2. The gas compositions obtained in period 1 and 2 were markedly different. Over the period from 2006 to 2007, the gas bubbles were depleted in He, H₂, and CH₄ of deep origin, but enriched in Ar and N₂ of atmospheric origin. This difference can be interpreted as being due to an irreversible change of the aquifer/gas system. The present deep component in the HOS gas is estimated to have composition He = 63 ppmV, H₂ = 37 ppmV, Ar = 0.17%, N₂ = 63%, and CH₄ = 37%. The new monitoring system is able to analyze the gas composition using a smaller volume of sample gas and with greater precision than the previous system. During the 3-month monitoring period 2, the separation capacity of the capillary column of the micro-GC was sufficiently maintained to determine gas-chromatographic peak areas for the five gaseous species examined. This study confirms that the new monitoring system with micro-GC is promising for continuous subsurface gas monitoring for earthquake prediction studies.     

Residual dune ridges: Sedimentary architecture,genesis, and implications for palaeo-climate reconstructions

Sedimentary architecture and genesis of residual dune ridges in a temperate climate are presented and implications for their use as archive of changes in long-term precipitation and wind climate are discussed. Residual dunes are common features of wet aeolian systems, where they form sets of shallow ridges, oriented perpendicular to the prevailing wind direction. Residual dune ridges of the study area are vegetated and typically elevate 0.6 to 2.5 m above the surrounding interdune flats. They develop on the lower stoss side of active transgressive dunes, triggered by periods of elevated groundwater table and hence colonization of the foot of the dune by rapid growing pioneer vegetation. Stabilized by plants, the growing ridge detaches from the active transgressive dune and gets abandoned within years in the course of the downwind-migration of the transgressive dune. Grain-size data suggest a main sediment supply from the transgressive dune and only minor input from other sources. Ground-penetrating radar reveals that the residual dune ridges are composed of windward-dipping as well as leeward-dipping sedimentary beds. Leeward-dipping strata reflect sediment supply from the parental dune, whereas windward-dipping beds are seen to result from sediment redistribution along the ridge and sediment supply from the adjacent swales during the ridge growth period. Multi-annual to multi-decadal variability in precipitation leads to the development of sequences composed of tens of ridges, spanning time periods of several centuries. Spacing of individual ridges in these sequences is controlled not by long-term variability in precipitation alone, but probably also reflects variable wind intensity which affects the migration rate of the parental dune. The important role of vegetation in ridge construction makes these landforms a demonstrative example of landscape development by geo-biosphere interacting processes.