Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “geophysical data”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Multichannel Analysis of Surface Waves Accelerated (MASWAccelerated): Software for efficient surface wave inversion using MPI and GPUs

Multichannel Analysis of Surface Waves (MASW) is a technique frequently used in geotechnical engineering and engineering geophysics to infer 1D layered models of seismic shear wave velocities in the top tens to hundreds of meters of the subsurface. We aim to accelerate MASW calculations by capitalizing on modern computer hardware available in the workstations of most engineers: multiple cores and graphics processing units (GPUs). We propose new parallel and GPU accelerated algorithms for computing 1D MASW inversion, and provide software implementations in C using Message Passing Interface (MPI) and CUDA. These algorithms take advantage of sparsity that arises in the problem, and the work balance between processes considers typical data trends. We compare our methods to an existing open source Matlab MASW tool. Our serial C implementation achieves a 2x speedup over the Matlab software, and we continue to see improvements by parallelizing the problem with MPI. Here we see nearly perfect strong and weak scaling for uniform data, and improve strong scaling for realistic data by repartitioning the problem to process mapping. By utilizing GPUs available on most modern workstations, we observe an additional 1.3x speedup over the serial C implementation on the first use of the method. We typically repeatedly evaluate theoretical dispersion curves as part of an optimization procedure, and on the GPU the kernel can be cached for faster reuse on later runs. We observe a 3.2x speedup on the cached GPU runs compared to the serial C runs. This work is the first open-source parallel or GPU-accelerated software tool for MASW imaging, and should enable geotechnical engineers to fully utilize all computer hardware at their disposal.

58 GEOSCIENCES↗
DOE OSTI2023

Baltimore Ecosystem Study: Soil atmosphere fluxes of carbon dioxide, nitrous oxide and methane, 1998 - ongoing

The Baltimore Ecosystem Study (BES) established a network of long-term permanent biogeochemical study plots in 1998. These plots provide long-term data on vegetation, soil and hydrologic processes in the key ecosystem types within the urban ecosystem. The network of study plots includes forest plots (upland and riparian), chosen to represent the range of forest conditions in the area and grass plots (to represent home lawns). Plots are instrumented with lysimeters (drainage and tension) to sample soil solution chemistry, time domain reflectometry probes to measure soil moisture, dataloggers to measure and record soil temperature, and trace gas flux chambers to measure the flux of carbon dioxide, nitrous oxide and methane from soil to the atmosphere. Measurements of in situ nitrogen mineralization, nitrification and denitrification were made at approximately monthly intervals from Fall 1998 - Fall 2000. Detailed vegetation characterization (all layers) was done in summer 1998 and 2015. Data from these plots has been published in Groffman et al. (2006, 2009), Groffman and Pouyat (2009), Savva et al. (2010), Costa and Groffman (2013), Duncan et al. (2013), Waters et al. (2014), Ni and Groffman (2018), Templeton et al. (2019). Literature Cited Costa, K.H. and P.M. Groffman. 2013. Factors regulating net methane flux in urban forests and grasslands. Soil Science Society of America Journal 77:850 - 855. Duncan, J. M., L. E. Band, and P. M. Groffman. 2013. Towards closing the watershed nitrogen budget: Spatial and temporal scaling of denitrification. Journal of Geophysical Research Biogeosciences 118:1-5; DOI: 10.1002/jgrg.20090 Groffman PM, Pouyat RV, Cadenasso ML, Zipperer WC, Szlavecz K, Yesilonis IC,. Band LE and Brush GS. 2006. Land use context and natural soil controls on plant community composition and soil nitrogen and carbon dynamics in urban and rural forests. Forest Ecology and Management 236:177-192. Groffman, P.M., C.O. Williams, R.V. Pouyat, L.E. Band and I.C. Yesilonis. 2009. Nitrate leaching and nitrous oxide flux in urban forests and grasslands. Journal of Environmental Quality 38:1848-1860. Groffman, P.M. and R.V. Pouyat. 2009. Methane uptake in urban forests and lawns. Environmental Science and Technology 43:5229-5235. DOI: 10.1021/es803720h. Ni, X. and P.M. Groffman. 2018. Declines in methane uptake in forest soils. Proceedings of the National Academies of Science of the United States of America 115:8587-8590. Savva, Y., K. Szlavecz, R. V. Pouyat, P. M. Groffman, and G. Heisler. 2010. Effects of land use and vegetation cover on soil temperature in an urban ecosystem. Soil Science Society of America Journal 74:469-480. Templeton, L., M.L. Cadenasso, J. Sullivan, M. Neel and P.M. Groffman. 2019. Changes in vegetation structure and composition of urban and rural forest patches in Baltimore from 1998 to 2015. Forest Ecology and Management. In press. Waters, E.R., J.L. Morse, N.D. Bettez and P.M. Groffman. 2014. Differential carbon and nitrogen controls of denitrification in riparian zones and streams along an urban to exurban gradient. Journal of Environmental Quality 43:955–963.

Groffman, Peter M↗
DOE OSTI2020

The EGS Collab Project: Learnings from Experiment 1

The primary objective of the EGS Collab Project sponsored by DOE is to increase the understanding needed to efficiently implement enhanced geothermal systems (EGS). One goal of the EGS Collab project is to create a collaborative research environment in which to study stimulation of crystalline rock at the 10 meter scale. Key to this effort is the collection of high quality data to allow comparison to numerical coupled process models in an effort to build confidence in the codes and modeling techniques used. In response to this, the EGS Collab team has created an underground test bed at the Sanford Underground Research Facility (SURF) in Lead SD at a depth of approximately 1.5 km to examine hydraulic fracturing (Experiment 1). We are currently designing a second test bed aimed at investigating shear stimulation (Experiment 2). At the Experiment 1 location, we have characterized our host rock using laboratory testing and numerous field-based geophysical and geological techniques, and created a well-instrumented test bed to allow us to carefully monitor stimulation events and flow tests. In addition to the installed geophysical sensors, we have used tracer tests, differences in the ambient microbial communities at flow collection locations, and cold water injection to inform us about dynamic flow pathways. In Experiment 1, we have hydraulically stimulated the host rock in a number times at several locations in one well, creating new fractures that connect to existing fractures between the injection and production boreholes. We have performed long-term ambient and chilled water injection tests as an analog to EGS, and have monitored system changes resulting from these water injections through geophysical monitoring, flow and pressure measurements, tracer tests, and microbiology. Here, we summarize the tests performed, issues identified including poroelastic and thermoelastic effects, Joule-Thomson effects, restarting effects, indications of flow channeling, and the primary learnings from Experiment 1.

Enhanced Geothermal Systems, EGS Collab, stimulati↗
DOE OSTI2023

Geothermal Play Fairway Analysis, Part 2: GIS methodology

Play Fairway Analysis (PFA) in geothermal exploration originates from a systematic methodology developed within the petroleum industry and is based on a geologic, geophysical, and hydrologic framework of identified geothermal systems. We tailored this methodology to study the geothermal resource potential of the Snake River Plain and surrounding region, but it can be adapted to other geothermal resource settings. We adapted the PFA approach to geothermal resource exploration by cataloging the critical elements controlling exploitable hydrothermal systems, establishing risk matrices that evaluate these elements in terms of both probability of success and level of knowledge, and building a code-based ‘processing model’ to process results. A geographic information system was used to compile a range of different data types, which we refer to as elements (e.g., faults, vents, heat flow, etc.), with distinct characteristics and measures of confidence. Discontinuous discrete data (points, lines, or polygons) for each element were transformed into continuous interpretive 2D grid surfaces called evidence layers. Because different data types have varying uncertainties, most evidence layers have an accompanying confidence layer which reflects spatial variations in these uncertainties. Risk layers, as defined here, are the product of evidence and confidence layers, and are the building blocks used to construct Common Risk Segment (CRS) maps for heat, permeability, and seal, using a weighted sum for permeability and heat, but a different approach with seal. CRS maps quantify the variable risk associated with each of these critical components. In a final step, the three CRS maps were combined into a Composite Common Risk Segment (CCRS) map, using a modified weighted sum, for results that reveal favorable areas for geothermal exploration. Additional maps are also presented that do not mix contributions from evidence and confidence (to allow an isolated view of evidence and confidence), as well as maps that calculate favorability using the product of components instead of a weighted sum (to highlight where all components are present). Our approach helped to identify areas of high geothermal favorability in the western and central Snake River Plain during the first phase of study and helped identify more precise local drilling targets during the second phase of work. By identifying favorable areas, this methodology can help to reduce uncertainty in geothermal energy exploration and development.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

Compositional evidence for chondrule origins of low‐Ca pyroxenes in comet Wild 2 and a giant cluster IDP

Abstract A literature compilation of 1136 low‐Ca pyroxene compositions from chondrules from 12 primitive type 2–3 carbonaceous, ordinary and enstatite chondrite groups define unique regions on an Al 2 O 3 and Cr 2 O 3 diagram when compared to low‐Ca pyroxenes from equilibrated type 4‐6 chondrites. Measured compositions of 100 low‐Ca pyroxenes from comet Wild 2 and a giant cluster IDP of probable cometary origin are similar to each other and fall in the type 2–3 chondrite chondrule region suggesting that most of the pyroxenes likely formed in the solar nebula like conventional chondrules. The data imply that most low Ca‐pyroxenes from comet Wild 2 and the giant cluster IDP formed from igneous crystallization processes and did not experience significant thermal metamorphism, indicating that the low‐Ca pyroxenes were unlikely incorporated into large parent bodies prior to accretion in their respective comet bodies. An intriguing group of nine low‐Ca pyroxenes from comet Wild 2 with low Cr and Al that fall where type 4–6 chondrites are located are interpreted as products of condensation. The compositional data combined with previously measured oxygen isotopes on 17 low‐Ca pyroxenes support earlier conclusions that comet samples have links with carbonaceous, ordinary, and possibly enstatite chondrite groups. Our results provide additional evidence that comets accreted materials from multiple chondrule reservoirs throughout the solar nebula.

Geochemistry & Geophysics↗
DOE OSTI2024

3D modeling of deep borehole electromagnetic measurements with energized casing source for fracture mapping at the Utah Frontier Observatory for Research in Geothermal Energy

Here, we present a 3D numerical modelling analysis evaluating the deployment of a borehole electromagnetic measurement tool to detect and image a stimulated zone at the Utah Frontier Observatory for Research in Geothermal Energy geothermal site. As the depth to the geothermal reservoir is several kilometres and the size of the stimulated zone is limited to several 100 m, surface-based controlled-source electromagnetic measurements lack the sensitivity for detecting changes in electrical resistivity caused by the stimulation. To overcome the limitation, the study evaluates the feasibility of using a three-component borehole magnetic receiver system at the Frontier Observatory for Research in Geothermal Energy site. To provide sufficient currents inside and around the enhanced geothermal reservoir, we use an injection well as an energized casing source. To efficiently simulate energizing the injection well in a realistic 3D resistivity model, we introduce a novel modelling workflow that leverages the strengths of both 3D cylindrical-mesh-based electromagnetic modelling code and 3D tetrahedral-mesh-based electromagnetic modelling code. The former is particularly well-suited for modelling hollow cylindrical objects like casings, whereas the latter excels at representing more complex 3D geological structures. In this workflow, our initial step involves computing current densities along a vertical steel-cased well using a 3D cylindrical electromagnetic modelling code. Subsequently, we distribute a series of equivalent current sources along the well's trajectory within a complex 3D resistivity model. We then discretize this model using a tetrahedral mesh and simulate the borehole electromagnetic responses excited by the casing source using a 3D finite-element electromagnetic code. This multi-step approach enables us to simulate 3D casing source electromagnetic responses within a complex 3D resistivity model, without the need for explicit discretization of the well using an excessive number of fine cells. We discuss the applicability and limitations of this proposed workflow within an electromagnetic modelling scenario where an energized well is deviated, such as at the Frontier Observatory for Research in Geothermal Energy site. Using the workflow, we demonstrate that the combined use of the energized casing source and the borehole electromagnetic receiver system offer measurable magnetic field amplitudes and sensitivity to the deep localized stimulated zone. The measurements can also distinguish between parallel-fracture anisotropic reservoirs and isotropic cases, providing valuable insights into the fracture system of the stimulated zone. Besides the magnetic field measurements, vertical electric field measurements in the open well sections are also highly sensitive to the stimulated zone and can be used as additional data for detecting and imaging the target. We can also acquire additional multiple-source data by grounding the surface electrode at various locations and repeating borehole electromagnetic measurements. This approach can increase the number of monitoring data by several factors, providing a more comprehensive dataset for analysing the deep-localized stimulated zone. The numerical analysis indicates that it is feasible to use the combination of the energized casing and downhole electromagnetic measurements in monitoring localized stimulated zone at large depths.

58 GEOSCIENCES↗
DOE OSTI2025

Impact of Salinity on Ground Ice Distribution Across an Arctic Coastal Polygonal Tundra Environment

The heterogeneous distribution of ground ice in the Arctic is a key driver of uneven ground subsidence as permafrost thaws, significantly impacting infrastructure and surface/subsurface hydrology. These topographic and hydrological changes contribute to major uncertainties in energy and carbon fluxes and storage in a warming Arctic. This study aims to improve our understanding of the controls on ground ice and organic matter distribution within the top 3 m of permafrost in coastal polygonal tundra near Utqiagvik, Alaska. To this end, we apply a neural network approach to bulk density distributions derived from nondestructive X-ray tomography of soil cores, trained with laboratory analyses, to improve the resolution and spatial coverage of estimates of dry bulk density, ice content, and organic matter content. In addition, we use capacitively coupled geophysical imaging to map soil electrical conductivity and salinity variations. The results show that sedimentary deposits from ocean transgressions, along with subsequent ice wedge polygon geomorphological processes, jointly influence the distribution of ice content at various scales. The impact of the latter decreases with depth, whereas the influence of salinity and sedimentary history increases. Although the controls on the distribution of soil organic matter content (g/cm 3 ) remain unclear, the pronounced heterogeneity in bulk density strongly influences its calculation from laboratory mass fraction measurements (g/g). From a methodological perspective, the interdependencies among soil components and the need for increased data coverage underscore the value of high-resolution density measurements, such as using X-ray tomography. Overall, this study emphasizes the importance of considering salinity constraints on ice content distribution in coastal permafrost regions. The results are expected to aid in the development of data products and process representations in geomorphological and ecosystem models.

Dafflon, Baptiste [Lawrence Berkeley National Labo↗
DOE OSTI2023

Using a physics-informed neural network and fault zone acoustic monitoring to predict lab earthquakes

Abstract Predicting failure in solids has broad applications including earthquake prediction which remains an unattainable goal. However, recent machine learning work shows that laboratory earthquakes can be predicted using micro-failure events and temporal evolution of fault zone elastic properties. Remarkably, these results come from purely data-driven models trained with large datasets. Such data are equivalent to centuries of fault motion rendering application to tectonic faulting unclear. In addition, the underlying physics of such predictions is poorly understood. Here, we address scalability using a novel Physics-Informed Neural Network (PINN). Our model encodes fault physics in the deep learning loss function using time-lapse ultrasonic data. PINN models outperform data-driven models and significantly improve transfer learning for small training datasets and conditions outside those used in training. Our work suggests that PINN offers a promising path for machine learning-based failure prediction and, ultimately for improving our understanding of earthquake physics and prediction.

42 ENGINEERING↗
DOE OSTI2020

Physical Properties and Gas Hydrate at a Near-Seafloor Thrust Fault, Hikurangi Margin, New Zealand

The Papaku Fault Zone, drilled at International Ocean Discovery Program (IODP) Site U1518, is an active splay fault in the frontal accretionary wedge of the Hikurangi Margin. In logging-while-drilling data, the 33-m-thick fault zone exhibits mixed modes of deformation associated with a trend of downward decreasing density, P-wave velocity, and resistivity. Methane hydrate is observed from ~30 to 585 m below seafloor (mbsf), including within and surrounding the fault zone. Hydrate accumulations are vertically discontinuous and occur throughout the entire logged section at low to moderate saturation in silty and sandy centimeter-thick layers. In this paper, we argue that the hydrate distribution implies that the methane is not sourced from fluid flow along the fault but instead by local diffusion. This, combined with geophysical observations and geochemical measurements from Site U1518, suggests that the fault is not a focused migration pathway for deeply sourced fluids and that the near-seafloor Papaku Fault Zone has little to no active fluid flow.

58 GEOSCIENCES↗
DOE OSTI2023

Behavior of Hydrogarnet‐Type Defects in Hydrous Stishovite at Various Temperatures and Pressures

Abstract Dense polymorphs of silica have been demonstrated experimentally to incorporate from 1.5 wt% to as much as 11.6 wt% H 2 O as OH groups, with implications for the hydrogen budgets of Earth and other planets. This OH is thought to enter the SiO 2 structure via a charge‐balanced substitution in which silicon vacancies (V Si ) are compensated by protonating four of the surrounding six oxygen atoms, often referred to as a hydrogarnet‐type defect. There are many possible configurations for this defect structure in dense silica, but the nature of these configurations and whether they can be distinguished experimentally is unknown. We present here density functional theory calculations that systematically assess the possible configurations of a hydrogarnet‐type defect in stishovite (rutile‐type SiO 2 ), with direct comparisons to experimental vibrational spectroscopy data. We predict that stishovite synthesized at 450 K and 10 GPa quenched to room temperature is dominated by a single defect type with tetrahedral geometry. This leads to OH stretching modes (2,500–3,000 cm −1 ) and SiOH bending modes (∼1,400–1,450 cm −1 ) largely consistent with experimentally observed modes. One remaining issue is that our calculations produce results compatible with experimental data on H to D exchange, but do not explain why a considerable fraction of the 1,420 cm −1 mode shifts by only 40 cm −1 in deuterated samples. At elevated pressures and temperatures, we find that a second square planar defect configuration also becomes favorable, leading to modes that should allow differentiation from the tetrahedral configuration.

Geochemistry & Geophysics↗
DOE OSTI2025

Multi-frequency progressive refinement for learned inverse scattering

Interpreting scattered acoustic and electromagnetic wave patterns is a computational task that enables remote imaging in a number of important applications, including medical imaging, geophysical exploration, sonar and radar detection, and nondestructive testing of materials. However, accurately and stably recovering an inhomogeneous medium from far-field scattered wave measurements is a computationally difficult problem, due to the nonlinear and non-local nature of the forward scattering process. We design a neural network, called Multi-Frequency Inverse Scattering Network (MFISNet), and a training method to approximate the inverse map from far-field scattered wave measurements at multiple frequencies. We consider three variants of MFISNet, with the strongest performing variant inspired by the recursive linearization method — a commonly used technique for stably inverting scattered wavefield data — that progressively refines the estimate with higher frequency content. MFISNet outperforms past methods in regimes with high-contrast, heterogeneous large objects, and inhomogeneous unknown backgrounds.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Basal Melting and Oceanic Observations Beneath Central Fimbulisen, East Antarctica

Abstract Basal melting of ice shelves is fundamental to Antarctic ice sheet mass loss, yet direct observations remain sparse. We present the first year‐round melt record (2017–2021) from a phase‐sensitive radar on Fimbulisen, one of the fastest flowing ice shelves in Dronning Maud Land, East Antarctica. The observed long‐term mean ablation rate at 350 m depth below the central ice shelf was 1.0 ± 0.5 m yr −1 , marked by substantial sub‐weekly variability ranging from 0.4 to 3.5 m yr −1 . 36‐h filtered basal melt rate fluctuations closely align with ocean velocity. On seasonal time scales, melt rates peak during austral spring to autumn (September–March), driven by both elevated ocean velocities and thermal driving near the base. The combined effect of thermal driving and current speed explains the majority of the melt rate variability ( r = 0.84), highlighting the dominant role of shear‐driven turbulence. This relationship enables parameterization of melt rates for the decade‐long ocean record (2010–2021), although deviations appear under low and high forcing conditions. Both observed and parameterized melt rates show similar yearly mean magnitudes compared to satellite‐derived melt rates but with a tenfold lower seasonal amplitude and a 3‐month delay in seasonality. These detailed concurrent ice–ocean observations provide essential validation data for remote sensing and numerical models that aim to quantify and project ice‐shelf response to a change in ocean forcing. In situ measurements and continued monitoring are crucial for accurately assessing and modeling future basal melt rates, and for understanding the complex dynamics driving ice‐shelf stability and sea‐level change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Heat capacity and thermodynamic functions of partially dehydrated sodium and zinc zeolite A (LTA)

Zeolite A (LTA) is an industrially important zeolite that exhibits sorption-induced framework flexibility, the thermodynamics of which are poorly understood. In this work, we report heat capacity measurements on zinc and sodium zeolite A from 1.8 to 300 K and compare the heat capacity of water in sodium zeolite A with that of water in other zeolites. The heat capacity of zeolitic water varies significantly depending on the hydration level and identity of the host zeolite, and more tightly bound water exhibits strong inflections in its heat capacity curve. This suggests a combination of effects, including differences in water-framework binding strength and hydration-dependent flexibility transitions. We also report fits of the heat capacity data using theoretical functions, and we report values for $C_{P,m}^°$,$Δ_0^TS_m^°$,$Δ_0^TH_m^°$, and $Φ_m^°$ from 0 to 300 K. These results contribute to a systematic thermodynamic understanding of the effects of cation exchange, guest molecule confinement, and sorbate-dependent flexibility transitions in zeolites.

Geochemistry & Geophysics↗
DOE OSTI2024

Global Geo-processed Data of Aquifer Properties by 0.5° Grid, Country and Water Basins

This repository of global hydrogeologic datasets contains aquifer properties on 0.5° scale, including depth to groundwater (Fan et al., 2013), aquifer thickness (de Graaf et al., 2015), WHYMap aquifer classes (Richts et al., 2011), recharge (Döll and Fiedler, 2008; Gleeson et al., 2016), lakes (Messager et al., 2016), porosity and permeability (Gleeson et al., 2014), digitized and geo-processed from their respective sources. Globally gridded aquifer properties could be used independently to estimate global groundwater availability or used as critical inputs to the superwell model to simulate groundwater extraction and provide estimates of pumped volumes and unit costs under user-specific scenarios. Key resources related to this data are: Niazi, H., Ferencz, S. B., Graham, N. T., Yoon, J., Wild, T. B., Hejazi, M., Watson, D. J., & Vernon, C. R. (2025). Long-term hydro-economic analysis tool for evaluating global groundwater cost and supply: Superwell v1.1. Geoscientific Model Development, 18(5), 1737-1767. https://doi.org/10.5194/gmd-18-1737-2025 superwell model repository which uses this data to simulate groundwater extraction and provides estimates of the global extractable volumes and unit-costs ($/km3) of accessible groundwater production under user-specified extraction scenarios. Repository Overview Main output: aquifer_properties_rec.csv contains all processed outputs, including aquifer properties like porosity, permeability, recharge, lake areas, aquifer thickness, and depth to groundwater. shapefiles.zip: contains all digitized GIS databases and shapefile for all aquifer properties prep_inputs.R and prep_inputs_recharge_lakes.R: R scripts that process the shapefiles to produce the aquifer_properties_rec.csv file plot_inputs.R: R script for plotting the maps and conducting preliminary analysis on the available groundwater volume basin_to_country_mapping.csv, basin_country_region_mapping.csv and continent_county_mapping.csv provide the mapping between continents, 32 energy-economic macro regions, countries, and water basins for post-processing aquifer_properties_rec.csv Maps: Each map visualizes the spatial distribution of one of the aquifer properties across the globe map_in_Porosity.png map_in_Permeability.png map_in_Aquifer_thickness.png map_in_Depth_to_water.png map_in_Recharge.png map_in_Grid_area_km.png map_in_Lake_area_km.png map_in_WHYClass.png Sample inputs sample_inputs.py: this script samples inputs from the aquifer_properties_rec dataset, ensuring the sampled and original inputs maintain the same distributions sampled_data_100.csv contains 100 sampled data points and sampled_data_100.png compares their distributions Dataset Overview The main outputs are consolidated in a comprehensive aquifer_properties_rec.csv file and include the following fields: GridCellID: Unique identifier for each (roughly 0.5°) grid cell Continent: Continent name Country: Country name GCAM_basin_ID: Identifier for GCAM hydrologic basin Basin_long_name: Full name of the basin WHYClass: Hydrogeologic classification based on WHYMap aquifer classes (Richts et al., 2011) Porosity: Soil porosity (%) (Gleeson et al., 2014) Permeability: Soil permeability (in square meters; Gleeson et al., 2014) Aquifer_thickness: Thickness of the aquifer (in meters; de Graaf et al., 2015) Depth_to_water: Depth to groundwater (in meters; Fan et al., 2013) Recharge: long-term annual averaged recharge rates (in m/yr; Döll and Fiedler, 2008; Gleeson et al., 2016) Grid_area: Area of the grid cell (in square meters) Lakes_area: Area of inland lakes (in square meters; Messager et al., 2016) Key References The datasets are digitized versions of global hydrogeologic properties from the following key literature sources: Depth to Groundwater: Fan, Y., Li, H., & Miguez-Macho, G. (2013). Global Patterns of Groundwater Table Depth. Science, 339(6122), 940-943. https://doi.org/10.1126/science.1229881 Aquifer Thickness: de Graaf, I. E. M., Sutanudjaja, E. H., van Beek, L. P. H., & Bierkens, M. F. P. (2015). A high-resolution global-scale groundwater model. Hydrol. Earth Syst. Sci., 19(2), 823-837. https://doi.org/10.5194/hess-19-823-2015 Porosity and Permeability: Gleeson, T., Moosdorf, N., Hartmann, J., & van Beek, L. P. H. (2014). A glimpse beneath earth's surface: GLobal HYdrogeology MaPS (GLHYMPS) of permeability and porosity. Geophysical Research Letters, 41(11), 3891-3898. https://doi.org/10.1002/2014GL059856 Aquifer classes: Richts, A., Struckmeier, W. F., & Zaepke, M. (2011). WHYMAP and the Groundwater Resources Map of the World 1:25,000,000. In J. A. A. Jones (Ed.), Sustaining Groundwater Resources: A Critical Element in the Global Water Crisis (pp. 159-173). Springer Netherlands. https://doi.org/10.1007/978-90-481-3426-7_10 Recharge: Döll, P., & Fiedler, K. (2008). Global-scale modeling of groundwater recharge. Hydrol. Earth Syst. Sci., 12(3), 863-885. https://doi.org/10.5194/hess-12-863-2008; Gleeson, T., Befus, K. M., Jasechko, S., Luijendijk, E., & Cardenas, M. B. (2016). The global volume and distribution of modern groundwater. Nature Geoscience, 9(2), 161-167. https://doi.org/10.1038/ngeo2590 Inland Lakes: Messager, M. L., Lehner, B., Grill, G., Nedeva, I., & Schmitt, O. (2016). Estimating the volume and age of water stored in global lakes using a geo-statistical approach. Nature Communications, 7(1), 13603. https://doi.org/10.1038/ncomms13603 Cite as Niazi, H., Watson, D., Hejazi, M., Yonkofski, C., Ferencz, S., Vernon, C., Graham, N., Wild, T., & Yoon, J. (2024). Global Geo-processed Data of Aquifer Properties by 0.5° Grid, Country and Water Basins. MultiSector Dynamics-Living, Intuitive, Value-adding, Environment. https://doi.org/10.57931/2484226 Contact Reach out to Hassan Niazi or open an issue in the superwell repository for questions or suggestions.

aquifer thickness↗
DOE OSTI2020

Heat Flow on the U.S. Beaufort Margin, Arctic Ocean: Implications for Ocean Warming, Methane Hydrate Stability, and Regional Tectonics

The Arctic Ocean represents one of the last geological frontiers on Earth. The formation of the Arctic Ocean remains unclear, and although it is well understood that the Arctic surface temperature is warming at a rate approximately twice as fast as the rest of the Earth, it is unclear how deep Arctic Ocean temperatures are changing. Understanding whether deep Arctic Ocean water is warming is important, since it can lead to the breakdown and release of huge quantities of frozen methane molecules trapped below the Arctic seafloor, which can increase ocean acidity and destabilize the seafloor. To understand both Arctic Ocean temperature change and geologic evolution, we collected temperature measurements at 103 sites. These measurements tell us how temperature increases with depth below the seafloor and can be used to understand both ocean temperature change and regional geology. Analysis of our data shows that at depths of 300–900 m below sea level, the Arctic Ocean has been warming steadily for perhaps several decades—nearly twice as deep as previous studies suggest. At ocean depths greater than 1,000 m, our analysis also reveals surprisingly high temperature increases with depth in the seafloor. The cause of these significant increases is unclear.

58 GEOSCIENCES↗
DOE OSTI2022

Trans-crustal structural control of CO 2 -rich extensional magmatic systems revealed at Mount Erebus Antarctica

Erebus volcano, Antarctica, with its persistent phonolite lava lake, is a classic example of an evolved, CO 2 -rich rift volcano. Seismic studies provide limited images of the magmatic system. Here we show using magnetotelluric data that a steep, melt-related conduit of low electrical resistivity originating in the upper mantle undergoes pronounced lateral re-orientation in the deep crust before reaching shallower magmatic storage and the summit lava lake. The lateral turn represents a structural fault-valve controlling episodic flow of magma and CO 2 vapour, which replenish and heat the high level phonolite differentiation zone. This magmatic valve lies within an inferred, east-west structural trend forming part of an accommodation zone across the southern termination of the Terror Rift, providing a dilatant magma pathway. Unlike H 2 O-rich subduction arc volcanoes, CO 2 -dominated Erebus geophysically shows continuous magmatic structure to shallow crustal depths of < 1 km, as the melt does not experience decompression-related volatile supersaturation and viscous stalling.

58 GEOSCIENCES↗
DOE OSTI2023

Computed Tomography Scanning and Geophysical Measurements of the CarbonSAFE Seal Integrity Wells in the Illinois Basin

The U.S. Department of Energy’s (DOE) National Energy Technology Laboratory (NETL) researchers in Morgantown, West Virginia, utilized computed tomography (CT) facilities and the Multi-Sensor Core Logger (MSCL) to evaluate the integrity of three cores (Geophysical Monitor #2 (GM2), Verification Well #1 (VW #1), and Wabash #1 for the Illinois Basin Carbon Storage Assurance Facility Enterprise (CarbonSAFE) efforts. This work was carried out in collaboration with the Illinois State Geological Survey (ISGS) as part of their efforts to analyze cores from two field locations, including the Illinois Basin-Decatur Project (IBDP) and second stage Illinois Industrial Carbon Capture and Storage (IL-ICCS) sites at Decatur, Illinois, and the Wabash CarbonSAFE project at the Wabash Valley Resources Integrated Gasification Combined Cycle plant in Terre Haut, Indiana. The results of this study are presented in several formats and are available online on the Energy Data eXchange (EDX). The rock characterization was conducted using non-destructive techniques, allowing for future analysis of the cores. While the equipment used did not provide direct visualization of shale pore space, it allowed for detection of fractures and discontinuities. Low resolution CT imagery with the NETL medical CT scanner was performed on the entire core. Qualitative analysis of the medical CT images, coupled with x-ray fluorescence (XRF), P-wave, and magnetic susceptibility measurements from the MSCL were utilized to identify areas of interest for further study as well as fractured zones. The ability to quickly identify key areas for more detailed study with higher resolution will save time and resources in future studies. The combination of methods used provided a multi-scale analysis of this core and provided both a macro and micro description of the core that is relevant for many subsurface energy-related examinations traditionally performed at NETL.

58 GEOSCIENCES↗