SEARCH · USGS Science

Results for “Geophysics”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,369 records · Page 76Linked to original sources

USGS2008

Hydrogeologic Framework in Three Drainage Basins in the New Jersey Pinelands, 2004-06

The U.S. Geological Survey, in cooperation with the New Jersey Pinelands Commission, began a multi-phase hydrologic investigation in 2004 to characterize the hydrologic system supporting the aquatic and wetland communities of the New Jersey Pinelands area (Pinelands). The Pinelands is an ecologically diverse area in the southern New Jersey Coastal Plain underlain by the Kirkwood-Cohansey aquifer system. The demand for ground water from this aquifer system is increasing as local development increases. To assess the effects of ground-water withdrawals on Pinelands stream and wetland water levels, three drainage basins were selected for detailed hydrologic assessments, including the Albertson Brook, McDonalds Branch and the Morses Mill Stream basins. Study areas were defined surrounding the three drainage basins to provide sub-regional hydrogeologic data for the ground-water flow modeling phase of this study. In the first phase of the hydrologic assessments, a database of hydrogeologic information and a hydrogeologic framework model for each of the three study areas were produced. These framework models, which illustrate typical hydrogeologic variations among different geographic subregions of the Pinelands, are the structural foundation for predictive ground-water flow models to be used in assessing the hydrologic effects of increased ground-water withdrawals. During 2004-05, a hydrogeologic database was compiled using existing and new geophysical and lithologic data including suites of geophysical logs collected at 7 locations during the drilling of 21 wells and one deep boring within the three study areas. In addition, 27 miles of ground-penetrating radar (GPR) surface geophysical data were collected and analyzed to determine the depth and extent of shallow clays in the general vicinity of the streams. On the basis of these data, the Kirkwood-Cohansey aquifer system was divided into 7 layers to construct a hydrogeologic framework model for each study area. These layers are defined by their predominant sediment textures as aquifers and leaky confining layers. The confining layer at the base of the Kirkwood-Cohansey aquifer system, depending on location, is defined as one of two distinct clays of the Kirkwood Formation. The framework models are described using hydrogeologic sections, maps of structure tops of layers, and thickness maps showing variations of sediment textures of the various model layers. The three framework models are similar in structure but unique to their respective study areas. The hydraulic conductivity of the Kirkwood-Cohansey aquifer system in the vicinity of the three study areas was determined from analysis of 16 slug tests and 136 well-performance tests. The mean values for hydraulic conductivity in the three study areas ranged from about 84 feet per day to 130 feet per day. With the exception of the basal confining layers, the variable and discontinuous nature of clay layers within the Kirkwood-Cohansey aquifer system was confirmed by the geophysical and lithologic records. Leaky confining layers and discontinuous clays are generally more common in the upper part of the aquifer system. Although the Kirkwood-Cohansey aquifer system generally has been considered a water-table aquifer in most areas, localized clays in the aquifer layers and the effectiveness of the leaky confining layers may act to impede the flow of ground water in varying amounts depending on the degree of confinement and the location, duration, and magnitude of the hydraulic stresses applied. Considerable variability exists in the different sediment textures. The extent to which this hydrogeologic variability can be characterized is constrained by the extent of the available data. Thus, the hydraulic properties of the modeled layers were estimated on the basis of available horizontal hydraulic conductivity data and the range of sediment textures estimated from geophysical and lithologic data.

Scientific Investigations Report↗

USGS2010

Completion Summary for Well NRF-16 near the Naval Reactors Facility, Idaho National Laboratory, Idaho

In 2009, the U.S. Geological Survey in cooperation with the U.S. Department of Energy's Naval Reactors Laboratory Field Office, Idaho Branch Office cored and completed well NRF-16 for monitoring the eastern Snake River Plain (SRP) aquifer. The borehole was initially cored to a depth of 425 feet below land surface and water samples and geophysical data were collected and analyzed to determine if well NRF-16 would meet criteria requested by Naval Reactors Facility (NRF) for a new upgradient well. Final construction continued after initial water samples and geophysical data indicated that NRF-16 would produce chemical concentrations representative of upgradient aquifer water not influenced by NRF facility disposal, and that the well was capable of producing sustainable discharge for ongoing monitoring. The borehole was reamed and constructed as a Comprehensive Environmental Response Compensation and Liability Act monitoring well complete with screen and dedicated pump. Geophysical and borehole video logs were collected after coring and final completion of the monitoring well. Geophysical logs were examined in conjunction with the borehole core to identify primary flow paths for groundwater, which are believed to occur in the intervals of fractured and vesicular basalt and to describe borehole lithology in detail. Geophysical data also were examined to look for evidence of perched water and the extent of the annular seal after cement grouting the casing in place. Borehole videos were collected to confirm that no perched water was present and to examine the borehole before and after setting the screen in well NRF-16. Two consecutive single-well aquifer tests to define hydraulic characteristics for well NRF-16 were conducted in the eastern SRP aquifer. Transmissivity and hydraulic conductivity averaged from the aquifer tests were 4.8 x 103 ft2/d and 9.9 ft/d, respectively. The transmissivity for well NRF-16 was within the range of values determined from past aquifer tests in other wells near NRF of 4.4 x 102 to 5.1 x 105 ft2/d. Water samples were analyzed for metals, nutrients, total organic carbon, volatile organic compounds, semi-volatile organic compounds, herbicides, pesticides, polychlorinated biphenols, and radionuclides. All chloride, nitrate, and sulfate concentrations were less than background concentrations for the eastern SRP aquifer north of the NRF. Concentrations in water samples for most of the organic compounds and radionuclides were less than the reporting limits and reporting levels.

Scientific Investigations Report↗

USGS2015

Delineation of areas having elevated electrical conductivity, orientation and characterization of bedrock fractures, and occurrence of groundwater discharge to surface water at the U.S. Environmental Protection Agency Barite Hill/Nevada Goldfields Superfund site near McCormick, South Carolina

During October 2012 through March 2013, the U.S. Geological Survey (USGS), in cooperation with the U.S. Environmental Protection Agency (EPA) Region 4, Superfund Section, conducted borehole geophysical logging, surface geophysical surveys, and water-quality profiling in selected wells and areas to characterize or delineate the extent of elevated subsurface electrical conductivity at the EPA Barite Hill/Nevada Goldfields Superfund site near McCormick, South Carolina. Elevated electrical conductivity measured at the site may be related to native rock materials, waste rock disposal areas used in past operations, and (or) groundwater having elevated dissolved solids (primarily metals and major ions) related to waste migration. Five shallow screened wells and four open-borehole bedrock wells were logged by using a suite of borehole tools, and downhole water-quality profiles were recorded in two additional wells. Well depths ranged from about 26 to 300 feet below land surface. Surface geophysical surveys based on frequency-domain electromagnetic and distributed temperature sensing (DTS) techniques were used to identify areas of elevated electrical conductivity (Earth materials and groundwater) and potential high dissolved solids in groundwater and surface water on land and in areas along the northern unnamed tributary at the site. Results from the electromagnetic-induction logging of four selected wells near the Main Pit and one well located about 800 feet southeast of the Main Pit lake indicate that elevated electrical conductivity extends to a depth of about 110 feet below land surface. Groundwater-quality properties recorded in eight selected wells were highly variable, suggesting a broad spectrum of geochemical conditions and contaminant concentrations within the groundwater system. Ranges of field water-quality properties recorded from water-profiling of groundwater in all wells logged were as follows: pH, 3.1 to 9.2; specific conductance, 48 to 5,300 microsiemens per centimeter; dissolved oxygen, 0.2 to 4.4 milligrams per liter; and water temperature, 17.0 to 18.0 degrees Celsius. The highest specific conductance and lowest pH measurements were made in boreholes located between the Main Pit lake and the northern unnamed tributary. Conceptually, these wells may intercept elevated dissolved solids in groundwater leaking from the Main Pit lake along a flow path that discharges into the unnamed tributary to the north. Results from surface geophysical electromagnetic and fiber-optics surveys confirm areas of focused discharge of groundwater near the Main Pit lake along the northern unnamed tributary. The frequency-domain surface electromagnetic surveys also identified an area with higher levels of elevated electrical conductivity located northwest of the former Rainsford Pit area. Bedrock properties were characterized from borehole geophysical logs collected from three open-borehole bedrock wells. The mean strike azimuth of the borehole foliation data measured in bedrock well IR-1 was 221° (N. 41° E.), and the mean dip angle was 78° to the northwest. Dominant strike azimuth orientations of primary fractures measured in three boreholes were from 210° to 250° (N. 30° E. to N. 70° E.) with a mean dip of 68° northwest. Transmissivity estimates interpreted from the heat-pulse flowmeter data from bedrock well IR-1 were about 69 feet squared per day, and the radius of influence was estimated at about 640 feet.

South Carolina↗

USGS2020

Hydrogeologic characterization of the Hualapai Plateau on the western Hualapai Indian Reservation, northwestern Arizona

This study was developed to assess if groundwater from the western Hualapai Plateau could be used to supply developments in the Grand Canyon West area of the Hualapai Indian Reservation and to collect hydrogeologic data for future use in a numerical groundwater model for the reservation. Ground-based geophysical surveys; existing well, spring, and other hydrogeologic information from previous studies; and new well and spring data collected for this study were used to provide a better understanding of the hydrogeology of the western Hualapai Plateau. Surface geophysical data provided information on the depth and geologic structure of lower Paleozoic rock units and Proterozoic crystalline and metamorphic rocks that underlie the western Hualapai Plateau. The surface geophysical data and discharge information from springs were used to select a site to drill and develop the U.S. Geological Survey Hualapai Test Well. The Hualapai Test Well was drilled to understand the geophysical properties of geologic formations at depth. These data were used to verify the results of surface geophysical data and to evaluate if sufficient water was present in the Hualapai Test Well for potential groundwater development. The Hualapai Test Well was drilled to a depth of 2,468 feet and bottomed in Proterozoic granite. Water was expected in the lower part of the Muav Limestone, but water was not observed until the Tapeats Sandstone at a depth of 2,400 feet. The Tapeats Sandstone was determined to be confined with a hydrostatic head of over 900 feet. A 48-hour pumping test was conducted to determine aquifer properties. Low specific capacity indicated that although groundwater is present in the Tapeats Sandstone, well yields are likely to be small. A water-quality sample indicated the sample had a calcium, magnesium-bicarbonate water type with a total dissolved-solids concentration of 371 milligrams per liter. Alpha radioactivity of the sample, 18.3 picocuries per liter, exceeded the U.S. Environmental Protection Agency maximum contaminant level of 15 picocuries per liter for drinking water. Concentrations of iron and manganese in the water sample also exceeded the U.S. Environmental Protection Agency secondary maximum contaminant levels for drinking water. An inventory of wells and springs provided insight into the occurrence of groundwater on the western Hualapai Plateau. Data from 56 springs on and adjacent to the western Hualapai Plateau were compiled for this study, and new data were collected at 31 springs. Discharge from springs visited for this study ranged from dry to about 345 gallons per minute. The temporal data from springs, where repeat measurements were available, indicated that spring flow is highly variable and likely related to seasonal and annual precipitation. Water levels from 36 wells on and adjacent to the western Hualapai Plateau were compiled for this study, and new water levels were collected at 5 wells. The spring and well data in conjunction with the Hualapai Test Well results indicated that on the western Hualapai Plateau, bedrock aquifers have limited discrete flow paths that make extensive groundwater development unlikely.

Arizona↗

USGS2000

Hydrogeology of the regional aquifer near Flagstaff, Arizona, 1994-97

Sandstones, siltstones, and limestones that are Pennsylvanian to Permian in age underlie the southern part of the Colorado Plateau near Flagstaff, Arizona, and contain a complex regional aquifer that has become increasingly important as a source of water for domestic, municipal, and recreational uses. Ground-water flow in the regional aquifer is poorly understood in this area because (1) depth of the aquifer limits exploratory drilling and testing and (2) the geologic structure increases the complexity of the aquifer characteristics and the ground-water flow system. Four methods were used to improve the understanding of the hydrogeology of the regional aquifer near Flagstaff. • Remote-sensing techniques and geologic mapping provided data to identify many structural features that indicate a more complex structural environment and history than previously realized. • Data from surface-geophysical techniques that included ground-penetrating radar, seismic reflection and seismic refraction, and square-array resistivity, verified that some of the geologic structures expressed at land surface propagate deep into the subsurface and through the principal water-bearing zones of the regional aquifer at near-vertical angles. • A well and spring inventory, borehole-geophysical methods, and well and aquifer tests provided additional information relating aquifer and ground-water flow characteristics to geologic structure. • Water-chemistry data, which included major ion, nutrient, trace-element, and radioactive and stable-isotope analyses, provided an independent means of verifying the hydrogeologic characteristics of the aquifer and were used to determine recharge and discharge areas, groundwater movement, and ground-water age. Ground-water recharge occurs throughout the area but is greatest at higher altitudes where precipitation is greater and in areas where heavily fractured rock units of the aquifer are exposed. The estimated annual average recharge to the regional aquifer in the study area is about 290,000 acre-feet. Ground water flows laterally and vertically through pore spaces in the rock and along faults and other fractures from high-altitude areas in the southern part of the study area to regional drains north of the study area along the Little Colorado and Colorado Rivers, and to drains south of the study area along Oak Creek and the Verde Valley. Ground-water discharge in these areas—about 400,000 acre-feet per year—exceeds the annual recharge to the aquifer in the Flagstaff area, but ground water from areas outside the study area contributes to this discharge as well. The saturated thickness of the regional aquifer averages about 1,200 feet, and the amount of water in storage could be as much as 4,800,000 acrefeet, or about 10 percent of the total volume of the aquifer. The quality of water in the regional aquifer in terms of dissolved-solids concentrations is good for most uses throughout the area. Dissolvedsolids concentrations generally are less than 500 milligrams per liter. Water in the regional aquifer is primarily a calcium magnesium bicarbonate type. In some areas near the Rio de Flag, the water has significant nitrate and chloride components, which indicate direct recharge in these areas from the Rio de Flag. Oxygen and deuterium data indicate a common recharge source for water in the aquifer and that some sites receive recharge from surface waters where evaporation has occurred. Estimated carbon-14 ages and tritium activities indicate ground-water ages from less than 200 years in the Lake Mary area to more than 5,000 years in the Wupatki area. The regional aquifer is heterogeneous and anisotrophic and has a complex ground-water flow system. The most productive water-bearing material tends to be fine- to medium-grained sandstones, and ground-water flow and potential well yields are related to geologic structure. Fracturing associated with structural deformation increases recharge locally and also increases the potential for high well yields. Surface-geophysical techniques provided information on the orientation of high-angle, deep-seated structure in the saturated zone. Borehole-geophysical data identified horizontal to near-horizontal fractures as significant components of the fracture-flow system not apparent in the surface-geophysical data. Structural features that strike northwest appear to be areas that have the greatest potential for high well yields. A north-northeastward-striking structure may be just as promising, but additional data are needed to verify this relation.

Arizona↗

USGS2015

Geologic cross sections and preliminary geologic map of the Questa Area, Taos County, New Mexico

In 2011, the senior authors were contacted by Ron Gardiner of Questa, and Village of Questa Mayor Esther Garcia, to discuss the existing and future groundwater supply for the Village of Questa. This meeting led to the development of a plan in 2013 to perform an integrated geologic, geophysical, and hydrogeologic investigation of the Questa area by the New Mexico Bureau of Geology & Mineral Resources (NMBG), the U.S. Geological Survey (USGS), and New Mexico Tech (NMT). The NMBG was responsible for the geologic map and geologic cross sections. The USGS was responsible for a detailed geophysical model to be incorporated into the NMBG products. NMT was responsible for providing a graduate student to develop a geochemical and groundwater flow model. This report represents the final products of the geologic and geophysical investigations conducted by the NMBG and USGS. The USGS final products have been incorporated directly into the geologic cross sections. The objective of the study was to characterize and interpret the shallow (to a depth of approximately 5,000 ft) three-dimensional geology and preliminary hydrogeology of the Questa area. The focus of this report is to compile existing geologic and geophysical data, integrate new geophysical data, and interpret these data to construct three, detailed geologic cross sections across the Questa area. These cross sections can be used by the Village of Questa to make decisions about municipal water-well development, and can be used in the future to help in the development of a conceptual model of groundwater flow for the Questa area. Attached to this report are a location map, a preliminary geologic map and unit descriptions, tables of water wells and springs used in the study, and three detailed hydrogeologic cross sections shown at two different vertical scales. The locations of the cross sections are shown on the index map of the cross section sheet.

Open-File Report↗

USGS2008

Data-resolution matrix and model-resolution matrix for Rayleigh-wave inversion using a damped least-squares method

Inversion of multimode surface-wave data is of increasing interest in the near-surface geophysics community. For a given near-surface geophysical problem, it is essential to understand how well the data, calculated according to a layered-earth model, might match the observed data. A data-resolution matrix is a function of the data kernel (determined by a geophysical model and a priori information applied to the problem), not the data. A data-resolution matrix of high-frequency (>2 Hz) Rayleigh-wave phase velocities, therefore, offers a quantitative tool for designing field surveys and predicting the match between calculated and observed data. We employed a data-resolution matrix to select data that would be well predicted and we find that there are advantages of incorporating higher modes in inversion. The resulting discussion using the data-resolution matrix provides insight into the process of inverting Rayleigh-wave phase velocities with higher-mode data to estimate S-wave velocity structure. Discussion also suggested that each near-surface geophysical target can only be resolved using Rayleigh-wave phase velocities within specific frequency ranges, and higher-mode data are normally more accurately predicted than fundamental-mode data because of restrictions on the data kernel for the inversion system. We used synthetic and real-world examples to demonstrate that selected data with the data-resolution matrix can provide better inversion results and to explain with the data-resolution matrix why incorporating higher-mode data in inversion can provide better results. We also calculated model-resolution matrices in these examples to show the potential of increasing model resolution with selected surface-wave data. ?? Birkhaueser 2008.

Pure and Applied Geophysics↗

USGS2005

Applying petrophysical models to radar travel time and electrical resistivity tomograms: Resolution-dependent limitations

[1] Geophysical imaging has traditionally provided qualitative information about geologic structure; however, there is increasing interest in using petrophysical models to convert tomograms to quantitative estimates of hydrogeologic, mechanical, or geochemical parameters of interest (e.g., permeability, porosity, water content, and salinity). Unfortunately, petrophysical estimation based on tomograms is complicated by limited and variable image resolution, which depends on (1) measurement physics (e.g., electrical conduction or electromagnetic wave propagation), (2) parameterization and regularization, (3) measurement error, and (4) spatial variability. We present a framework to predict how core‐scale relations between geophysical properties and hydrologic parameters are altered by the inversion, which produces smoothly varying pixel‐scale estimates. We refer to this loss of information as “correlation loss.” Our approach upscales the core‐scale relation to the pixel scale using the model resolution matrix from the inversion, random field averaging, and spatial statistics of the geophysical property. Synthetic examples evaluate the utility of radar travel time tomography (RTT) and electrical‐resistivity tomography (ERT) for estimating water content. This work provides (1) a framework to assess tomograms for geologic parameter estimation and (2) insights into the different patterns of correlation loss for ERT and RTT. Whereas ERT generally performs better near boreholes, RTT performs better in the interwell region. Application of petrophysical models to the tomograms in our examples would yield misleading estimates of water content. Although the examples presented illustrate the problem of correlation loss in the context of near‐surface geophysical imaging, our results have clear implications for quantitative analysis of tomograms for diverse geoscience applications.

Journal of Geophysical Research B: Solid Earth↗

USGS1986

Archaeological prospection: An introduction to the Special Issue

Archaeological prospection, as the use of geophysical methods in archaeology is known in Europe, is about four decades old (seven decades, if aerial photography of archaeological sites is included). Virtually the entire range of geophysical methods, perhaps excluding only borehole techniques, has found application in the search for archaeological sites unseen or partially known. Pressures by developers, and the public's growing sensitivity toward the preservation of historic and prehistoric cultural artifacts and sites, has led to an accelerating use of high-resolution geophysical methods in the archaeological sciences. The archaeogeophysical articles in this Special Issue are reasonably representative of the development of this specialty field of geophysics.

Geophysics↗

USGS2020

Passive seismic survey of sediment thickness, Dasht-e-Nawar basin, eastern Afghanistan

Exploration of water resources is needed for public supply, extraction of mineral resources, and economic development in Afghanistan. Remotely-sensed data are useful for identifying the general nature of surface sediments, however, “boots on the ground” geophysical surveys or drilling programs are needed to quantify the thickness of sediments or aquifers. The nature of such investigations presents a risk to field crews that may prohibit exploration of potentially valuable aquifers or mineral resources. The Dasht-e-Nawar basin, in east-central Afghanistan, contains a 400 km2 playa that includes evaporative mineral deposits, particularly lithium, which has been of interest since the 1970s. However, exploration of the basin, as with many areas of Afghanistan, has been hampered by decades of conflict. In 2014, an investigation of the basin was conducted by the U.S. Department of Defense Task Force for Business and Stability Operations (TFBSO), and their contractor, in cooperation with the U.S. Geological Survey (USGS). For this investigation the USGS compared the results of a rapid passive seismic survey of basin sediment thickness to the results of an independently conducted gravity survey of the same area. Each point measurement for the passive seismic method requires less than 30 minutes in the field by one person. The technique utilizes ambient seismic noise without an external sound source such as required by traditional seismic surveys. Additionally, the technique does not require external sensor arrays, which can be kilometers long for some geophysical techniques. The passive seismic equipment used in this study weighs approximately 1 kilogram and is about 10 cm3 in size. Although relatively new for assessment of sediment thickness, several investigations have found this method to be capable of estimating sediment thicknesses, in the 10’s to 1000 meter range, in settings with unconsolidated sediment over bedrock and a contrast in acoustic impedance. In this investigation, the gravity survey was conducted during a period of 3 weeks by an experienced field crew; required a detailed, centimeter-scale land elevation survey; and required laboratory analyses of sediment and rock densities to interpret the gravity data. In contrast, the passive seismic survey was collected by two inexperienced operators over a period of 8 days and required no additional data to interpret. Due to security restrictions, USGS personnel could not visit the site and the seismic operator was trained immediately prior to the field work. Although the quality of the seismic survey was affected by strong afternoon winds, and by the inexperience of the field operator, the results were fairly comparable to the gravity survey. Similar basin sediment thicknesses and patterns in sediment thickness were identified by both surveys in the basin with an estimated maximum thickness of approximately 170 m. The passive seismic technique required substantially less field resources and time than would be required by other geophysical surveys. Although this method will not be effective in all geologic settings, it may be a valuable assessment tool for use before conducting other, more intensive, geophysical efforts or drilling programs, especially in regions with elevated security risks such as Afghanistan.

Dasht-e-Nawar basin↗

USGS2006

Nonlinear inversion of potential-field data using a hybrid-encoding genetic algorithm

Using a genetic algorithm to solve an inverse problem of complex nonlinear geophysical equations is advantageous because it does not require computer gradients of models or "good" initial models. The multi-point search of a genetic algorithm makes it easier to find the globally optimal solution while avoiding falling into a local extremum. As is the case in other optimization approaches, the search efficiency for a genetic algorithm is vital in finding desired solutions successfully in a multi-dimensional model space. A binary-encoding genetic algorithm is hardly ever used to resolve an optimization problem such as a simple geophysical inversion with only three unknowns. The encoding mechanism, genetic operators, and population size of the genetic algorithm greatly affect search processes in the evolution. It is clear that improved operators and proper population size promote the convergence. Nevertheless, not all genetic operations perform perfectly while searching under either a uniform binary or a decimal encoding system. With the binary encoding mechanism, the crossover scheme may produce more new individuals than with the decimal encoding. On the other hand, the mutation scheme in a decimal encoding system will create new genes larger in scope than those in the binary encoding. This paper discusses approaches of exploiting the search potential of genetic operations in the two encoding systems and presents an approach with a hybrid-encoding mechanism, multi-point crossover, and dynamic population size for geophysical inversion. We present a method that is based on the routine in which the mutation operation is conducted in the decimal code and multi-point crossover operation in the binary code. The mix-encoding algorithm is called the hybrid-encoding genetic algorithm (HEGA). HEGA provides better genes with a higher probability by a mutation operator and improves genetic algorithms in resolving complicated geophysical inverse problems. Another significant result is that final solution is determined by the average model derived from multiple trials instead of one computation due to the randomness in a genetic algorithm procedure. These advantages were demonstrated by synthetic and real-world examples of inversion of potential-field data. ?? 2005 Elsevier Ltd. All rights reserved.

Computers & Geosciences↗

USGS2026

Understanding the resource potential of natural hydrogen on Earth: Scientific gaps, uncertainties and recommendations

A comprehensive scientific research roadmap is essential to bridge knowledge gaps and deepen the understanding of key geological, geochemical, and geophysical aspects of natural hydrogen (H 2 ) as a potential new energy resource. This paper reviews major scientific uncertainties on natural H 2 , suggesting research priorities, as a guide for defining exploration strategies, techniques, and data interpretation. The uncertainties concern all phases of the natural H 2 cycle, from generation (source rocks) through migration (advection and diffusion) and accumulation (reservoir and cap rocks) to the application and interpretation of subsurface and surface geochemical and geophysical exploration techniques. Understanding H 2 sources and generation rates (the amount of H 2 generated by a given volume of rock over time) is crucial for determining whether a geological H 2 system operates as a short-term dynamic system with rapid H 2 production and release, or as a conventional gas system with long-term accumulations, analogous to petroleum reservoirs. Preliminary estimates for serpentinisation, radiolysis, and organic matter degradation suggest that H 2 generation is not inherently fast, especially for non-hydrothermal continental systems (crystalline basement of shields, ophiolites, peridotite massifs, sedimentary basins), and long-term accumulations, like those of fossil natural gas systems, represent the most likely scenario. The mechanisms of H 2 migration through geological formations require application of fundamental principles of fluid-flow physics, distinguishing advection and diffusion, as well as their forms (from gas-phase, bubble flows to aqueous solutions). Additional studies of H 2 accumulation and retention in subsurface reservoirs could improve understanding of mechanisms of H 2 migration by focusing on the rock fluid-bearing properties and the factors affecting H 2 preservation, such as the presence of cap rocks impermeable to H 2 , pressure conditions, residence times, and microbial or abiotic consumption. Advanced techniques, including reservoir modelling, flow simulations, 3D imaging (micro-CT) of H 2 -bearing rocks, and extraction and analysis of gas occluded in rocks, can provide insights into the stability and potential recoverability of H 2 accumulations. The interpretation of surface exploration techniques, including gas geochemistry, geophysics, and remote sensing, long employed in mineral and energy resource exploration, is now being adapted for natural H 2 studies, but challenges remain in the data interpretation. Distinguishing H 2 seepage due to geological degassing from H 2 produced near the surface by modern microbial processes or artificial sources, such as hammering or drilling for soil-gas sampling, drilling into aquifers, and corrosion in boreholes, is an essential step in exploration. The simple detection of H 2 in soils, even in morphological structures like sub-circular depressions or “fairy circles”, cannot be cursorily interpreted as a signal of natural H 2 seepage from a deep source. A holistic geochemical approach, including isotopic analyses of gases associated with H 2 , is recommended to distinguish among the variety of possible H 2 origins. Observations of H 2 in wells should be interrogated to rule out possible artifacts such as corrosion and drill bit metamorphism. The integration of multiple geophysical methods, including seismic, gravimetric, magnetic, and electro-magnetic surveys, is recommended to mitigate interpretation ambiguities regarding the structure of a subsurface H 2 system (source and reservoir rocks, including fluid and gas storage), due to the non-uniqueness of rock-specific physical properties.

Earth-Science Reviews↗

USGS2024

Geochemical and geochronologic evidence for a contiguous northeastern Wyoming Province

The extent and nature of the Wyoming Province, an Archean craton in southwestern Laurentia, are poorly understood due to limited exposure between spatially isolated basement-cored uplifts. This lack of exposure has led to debate about whether the northeastern Wyoming Province is underlain by contiguous Archean crust or Proterozoic rocks and suture zone associated with the Trans-Hudson orogeny. To assess these models we analyzed samples recovered from drill cores in the buried northeastern Wyoming Province that straddle the proposed Proterozoic suture. Whole-rock geochemical and Nd-Pb-O isotopic data suggest the rocks formed by partial melting of > 3.5 Ga hydrated mafic to tonalitic sources similar to elsewhere in the northern Wyoming Province. Zircon U-Pb dates record Mesoarchean magmatism and metamorphism (3.0–2.8 Ga). Published whole-rock Rb-Sr, hornblende and biotite K-Ar, and apatite fission track dates suggest these rocks have not been heated above 300 °C since 2.5–2.1 Ga. Geophysical potential field data are consistent across northeastern Wyoming contrasting with major discontinuities associated with documented Proterozoic orogens on all other margins of the Wyoming Province. Hence, geochemical, isotopic, and geochronologic data along with geophysical imaging can be most simply interpreted in terms of continuous Archean crust in the northeastern Wyoming Province. A geophysically-imaged reflector east of the Bighorn Mountains may juxtapose Archean terranes with similar ages and sources, similar to the boundary between the Montana metasedimentary terrane and Beartooth-Bighorn magmatic zone in the northwestern Wyoming Province. This work emphasizes the value of integrating geologic and geophysical constraints to constrain Archean provinces and their tectonic evolution.

Colorado, Idaho, Montana, Nebraska, North Dakota, ↗

USGS2008

Electrical characterization of non‐Fickian transport in groundwater and hyporheic systems

Recent work indicates that processes controlling solute mass transfer between mobile and less mobile domains in porous media may be quantified by combining electrical geophysical methods and electrically conductive tracers. Whereas direct geochemical measurements of solute preferentially sample the mobile domain, electrical geophysical methods are sensitive to changes in bulk electrical conductivity (bulk EC) and therefore sample EC in both the mobile and immobile domains. Consequently, the conductivity difference between direct geochemical samples and remotely sensed electrical geophysical measurements may provide an indication of mass transfer rates and mobile and immobile porosities in situ. Here we present (1) an overview of a theoretical framework for determining parameters controlling mass transfer with electrical resistivity in situ; (2) a review of a case study estimating mass transfer processes in a pilot‐scale aquifer storage recovery test; and (3) an example application of this method for estimating mass transfer in watershed settings between streams and the hyporheic corridor. We demonstrate that numerical simulations of electrical resistivity studies of the stream/hyporheic boundary can help constrain volumes and rates of mobile‐immobile mass transfer. We conclude with directions for future research applying electrical geophysics to understand field‐scale transport in aquifer and fluvial systems subject to rate‐limited mass transfer.

Water Resources Research↗

USGS2011

Critical steps for the continuing advancement of hydrogeophysics

Special hydrogeophysics issues published by hydrology and geophysics journals, special sessions and workshops at conferences, and an increasing number of short courses demonstrate the growing interest in the use of geophysics for hydrologic investigations. The formation of the hydrogeophysics technical subcommittee of AGU's Hydrology section adds further evidence of the recognized significance of this growing interdisciplinary field. Given the clear value of nondestructive and nonintrusive imaging for subsurface investigations, we believe the advances in the adoption of existing geophysical methods, the development of novel methods, and the merging of geophysical and other data made in hydrogeophysics could be applied to a wide range of geological, environmental, and engineering applications.

Eos Science News↗

USGS2018

MoisturEC: a new R program for moisture content estimation from electrical conductivity data

Noninvasive geophysical estimation of soil moisture has potential to improve understanding of flow in the unsaturated zone for problems involving agricultural management, aquifer recharge, and optimization of landfill design and operations. In principle, several geophysical techniques (e.g., electrical resistivity, electromagnetic induction, and nuclear magnetic resonance) offer insight into soil moisture, but data‐analysis tools are needed to “translate” geophysical results into estimates of soil moisture, consistent with (1) the uncertainty of this translation and (2) direct measurements of moisture. Although geostatistical frameworks exist for this purpose, straightforward and user‐friendly tools are required to fully capitalize on the potential of geophysical information for soil‐moisture estimation. Here, we present MoisturEC, a simple R program with a graphical user interface to convert measurements or images of electrical conductivity (EC) to soil moisture. Input includes EC values, point moisture estimates, and definition of either Archie parameters (based on experimental or literature values) or empirical data of moisture vs. EC. The program produces two‐ and three‐dimensional images of moisture based on available EC and direct measurements of moisture, interpolating between measurement locations using a Tikhonov regularization approach.

Groundwater↗

USGS2011

Hydrogeophysical methods for analyzing aquifer storage and recovery systems

Hydrogeophysical methods are presented that support the siting and monitoring of aquifer storage and recovery (ASR) systems. These methods are presented as numerical simulations in the context of a proposed ASR experiment in Kuwait, although the techniques are applicable to numerous ASR projects. Bulk geophysical properties are calculated directly from ASR flow and solute transport simulations using standard petrophysical relationships and are used to simulate the dynamic geophysical response to ASR. This strategy provides a quantitative framework for determining site‐specific geophysical methods and data acquisition geometries that can provide the most useful information about the ASR implementation. An axisymmetric, coupled fluid flow and solute transport model simulates injection, storage, and withdrawal of fresh water (salinity ∼500 ppm) into the Dammam aquifer, a tertiary carbonate formation with native salinity approximately 6000 ppm. Sensitivity of the flow simulations to the correlation length of aquifer heterogeneity, aquifer dispersivity, and hydraulic permeability of the confining layer are investigated. The geophysical response using electrical resistivity, time‐domain electromagnetic (TEM), and seismic methods is computed at regular intervals during the ASR simulation to investigate the sensitivity of these different techniques to changes in subsurface properties. For the electrical and electromagnetic methods, fluid electric conductivity is derived from the modeled salinity and is combined with an assumed porosity model to compute a bulk electrical resistivity structure. The seismic response is computed from the porosity model and changes in effective stress due to fluid pressure variations during injection/recovery, while changes in fluid properties are introduced through Gassmann fluid substitution.

Groundwater↗

USGS1970

Correlation of aeromagnetics and aeroradioactivity with lithology in the Spotsylvania area, Virginia

The U.S. Geological Survey made a detailed aeromagnetic and aeroradioactivity survey of 1050 sq mi of the Spotsylvania area in the Piedmont province of eastern Virginia. The study area consists of sixteen 7 1/2-minute quadrangles in Spotsylvania and adjoining counties, Virginia. East-west traverses were flown at a spacing of 1/2 mi and a flight elevation of 500 ft above ground. Geophysical data were compiled at a scale of 1:24,000, with a contour interval of 25 gammas for the magnetics and 25 counts per sec for the aeroradioactivity. The geophysical data were checked in the field, and a lithologic map was made. Correlation is excellent between geophysical data and lithologic units, and it is evident that good detailed aeromagnetic and aeroradioactivity maps can provide appreciable shortcuts in mapping an area such as the Piedmont, where lithologies are complex and outcrops sparse. Radioactivity highs are generally associated with magnetic lows, and more importantly, for purposes of geologic interpretation where the magnetic expression is flat and featureless, there is considerable radioactivity detail for geologic guidance. By using the two geophysical maps in conjunction, the geology between widely spaced outcrops can be worked out. In the study area, high radioactivity and low magnetic intensity characterize granitic rocks, many of which were not previously mapped. Mafic bodies are characterized by high magnetic values and low radioactivity. Much of the area is underlain by a sequence of metamorphosed sedimentary and volcanic rocks having alternating quartz-muscovite and magnetiferous chlorite-actinolite assemblages. © 1970, The Geological Society of America, Inc.

Virginia↗