Search NASA⌕ Search

SEARCH · Search NASA

Results for “Geologic Map”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Global ocean dimethyl sulfide climatology estimated from observations and an artificial neural network

Marine dimethyl sulfide (DMS) is important to climate due to the ability of DMS to alter Earth's radiation budget. Knowledge of the global-scale distribution, seasonal variability, and sea-to-air flux of DMS is needed in order to improve understanding of atmospheric sulfur, aerosol/cloud dynamics, and albedo. Here we examine the use of an artificial neural network (ANN) to extrapolate available DMS measurements to the global ocean and produce a global climatology with monthly temporal resolution. A global database of 82 996 ship-based DMS measurements in surface waters was used along with a suite of environmental parameters consisting of latitude–longitude coordinates, time of day, time of year, solar radiation, mixed layer depth, sea surface temperature, salinity, nitrate, phosphate, and silicate. Linear regressions of DMS against the environmental parameters show that on a global-scale mixed layer depth and solar radiation are the strongest predictors of DMS. These parameters capture ~9 % and ~7 % of the raw DMS data variance, respectively. Multilinear regression can capture more of the raw data variance (~39 %) but strongly underestimates DMS in high-concentration regions. In contrast, the artificial neural network captures ~66 % of the raw data variance in our database. Like prior climatologies our results show a strong seasonal cycle in surface ocean DMS with the highest concentrations and sea-to-air fluxes in the high-latitude summertime oceans. We estimate a lower global sea-to-air DMS flux (20.12±0.43 Tg S yr –1 ) than the prior estimate based on a map interpolation method when the same gas transfer velocity parameterization is used. Our sensitivity test results show that DMS concentration does not change unidirectionally with each of the environmental parameters, which emphasizes the interactions among these parameters. The ANN model suggests that the flux of DMS from the ocean to the atmosphere will increase with global warming. Given that larger DMS fluxes induce greater cloud albedo, this corresponds to a negative climate feedback.

54 ENVIRONMENTAL SCIENCES↗
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 OSTI2024

Wellbore cement alteration and roles of CO 2 and shale during underground hydrogen storage

To mitigate climate change and adopt renewable energy, energy storage is crucial and can be done in the form of hydrogen gas (H 2 ). Subsurface geologic reservoirs are positioned to store H 2 on the largest scales for the longest terms of all potential options. However, H 2 injection may boost reactions that consume hydrogen, generate undesired gases, and alter pore structures of geomedia. To explore the extent of H 2 -associated biotic reactions at a near wellbore location, four experiments were conducted under underground storage conditions with wellbore cement cores and, in most instances, shale samples submerged in synthetic formation brine. Post-reaction gas, aqueous, and solid phase samples were analyzed using olfactory screening and, later, gas chromatography (GC-MS), inductively coupled plasma optical emission spectroscopy (ICP-OES), scanning electron microscopy (SEM), and synchrotron micro-scale x-ray fluorescence (μ-XRF). Within a period of 16 weeks, hydrogen sulfide (H 2 S) was generated in systems containing both H 2 and shale. XRF mapping identified a zone enriched in iron(II) and reduced sulfur along the rim of cement cross sections that was largely associated with CO 2 -induced cement carbonation. Shale did not show noticeable alteration, but there is evidence it contributed to the initial inoculation of the system and provided nutrients for microbes via water-rock interactions. Here, this study considers both rock formations and wellbore cement not previously evaluated concurrently. Findings support understanding and modeling of H 2 -associated biogeochemical reactions during underground hydrogen storage.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Small Scale Field Test Demonstrating CO 2 Sequestration In Arbuckle Saline Aquifer And By CO 2 -Eor At Wellington Field, Sumner County, Kansas

The objectives of this project are to understand the processes that occur when a maximum of 70,000 metric tonnes of CO 2 are injected into two different formations to evaluate the response in different lithofacies and depositional environments. The evaluation will be accomplished through the use of both in situ and indirect MVA (monitoring, verification, and accounting) technologies. The project will optimize for carbon storage accounting for 99% of the CO 2 using lab and field testing and comprehensive characterization and modeling techniques. Site characterization and CO 2 injection should demonstrate state-of-the-art MVA tools and techniques to monitor and visualize the injected CO 2 plume and to refine geomodels developed using nearly continuous core, exhaustive wireline logs, and well tests and a multi-component 3-D seismic survey. Reservoir simulation studies will map the injected CO 2 plume and estimate tonnage of CO 2 stored in solution, as residual gas, and by mineralization and integrate MVA results and reservoir models shall be used to evaluate CO 2 leakage. A rapid-response mitigation plan was developed to minimize CO 2 leakage and provide a comprehensive risk management strategy. The CO 2 was intended to be supplied from a reliable facility and have an adequate delivery and quality of CO 2 . However, several unforeseen circumstances complicated this plan: (1) the initially negotiated CO 2 supply facility went offline and contracts associated with CO 2 supply had to be renegotiated, (2) a UIC Class VI permit proved to be difficult to obtain due to the experimental nature of the project. Both subjects are detailed in separate deliverables attached to this report. The CO 2 enhanced oil recovery (EOR) and geologic storage in Mississippian carbonate reservoir was sucessully deployed. Approximately 20,000 metric tons of CO 2 was injected in the upper part of the Mississippian reservoir to verify CO 2 EOR viability in carbonate reservoirs and evaluate a potential of transitioning to geologic CO 2 storage through EOR. A total of 1,101 truckloads, 19,803 metric tons—an average of 120 tonnes per day—were delivered over the course of injection that lasted from January 9 to June 21, 2016. After cessation of CO 2 injection, the KGS 2-32 well was converted to water injector and continues to operate. CO 2 EOR progression in the field was monitored weekly with fluid level, temperature, and production recording and formation fluid composition sampling. It is important to note that normally, CO 2 EOR pilots are less efficient than commercial operations due to lack of directional and precise well control, lack of surface facilities for CO 2 recycling, and other factors. As a result of this pilot CO 2 injection, the observed incremental average oil production increase was ~68% with only ~18% of injected CO 2 produced back. Decline curve analysis forecasts of additional cumulative oil produced were 32.44M STB to the end of 2027. Wellington Mississippian pilot efficiency by the end of forecast calculations is 11 MCF per barrel of produced oil. Using 32M STB oil production and $1,964,063 cost of CO 2 , CO 2 EOR cost per barrel of oil production is ~$60. Simple but robust monitoring technologies proved to be very efficient in detecting and locating CO 2 . High CO 2 reservoir retentions with low yields within an actively producing field could help to estimate real-world risks of CO 2 geological storage for future projects. The Wellington Field CO 2 EOR was executed in a controlled environment with high efficiency. This case study proves that CO 2 EOR could be successfully applied in Kansas carbonate reservoirs if CO 2 sources and associated infrastructure are available. Recent developments in unconventional resources development in Mid-Continent USA and associated large volume disposal of backflow water and the resulting seismic activity have brought more focus and attention to the Arbuckle Group in southern Kansas. Despite the commercial interest, limited essential information about reservoir properties and structural elements has impeded the management and regulation of disposal, an issue brought to the forefront by recent seismicity in and near areas of large volumes and rates of brine disposal. The Kansas Geological Survey (KGS) collected, compiled, and analyzed available data, including well logs, core data, step rate tests, drill stem tests, 2-D and 3-D seismic data, water level measurements, and others types of data. Several exploratory wells were drilled and core was collected and modern suites of logs were analyzed. Reservoir properties were populated into several site-specific geological models. The geological models illustrate the highly heterogeneous nature of the Arbuckle Group. Vertical and horizontal variability results in several distinct hydro-stratigraphic units that are the result of both depositional and diagenetic processes. During the course of this project, it has been demonstrated that advanced seismic interpretation methods can be used successfully for characterization of the Mississippian reservoir and Arbuckle saline aquifer. Analysis of post-stack 3-D seismic data at the Mississippian reservoir showed the response of a gradational velocity transition. Pre-stack gather analysis showed that porosity zones of the Mississippian and Arbuckle reservoirs exhibit characteristic amplitude versus offset (AVO) response. Simultaneous AVO inversion estimated P- and S-impedances. The 3-D survey gather azimuthal anisotropy analysis (AVAZ) provided information about the fault and fracture network and showed good agreement to the regional stress field and well data. Mississippian reservoir porosity and fracture predictions agreed well with the observed mobility of injected CO 2 in KGS well 2-32. Fluid substitution modeling predicted acoustic impedance reduction in the Mississippian carbonate reservoir introduced by the presence of CO 2 . Seismicity in the United States midcontinent has increased by orders of magnitude over the past decade. Spatiotemporal correlations of seismicity to wastewater injection operations have suggested that injection-related pore fluid pressure increases are inducing the earthquakes. In this investigation, we examine earthquake occurrence in southern Kansas and northern Oklahoma and its relation to the change in pore pressure. The main source of data comes from the Wellington Array in the Wellington oil field, in Sumner County, Kansas, which has monitored for earthquakes in central Sumner County, Kansas, since early 2015. The seismometer array was established to monitor CO 2 injection operations at Wellington Field. Although no seismicity was detected in association with the spring 2016 Mississippian CO 2 injection, the array has recorded more than 2,500 earthquakes in the region and is providing valuable understanding to induced seismicity. A catalog of earthquakes was built from this data and was analyzed for spatial and temporal changes, stress information, and anisotropy information. The region of seismic concern has been shown to be expanding through use of the Wellington earthquake catalog, which has revealed a northward progression of earthquake activity reaching the metropolitan area of Wichita. The stress orientation was also calculated from this earthquake catalog through focal mechanism inversion. The calculated stress orientation was confirmed through comparison to other stress measurements from well data and previous earthquake studies in the region. With this knowledge of the stress orientation, the anisotropy in the basement could be understood. This allowed for the anisotropy measurements to be correlated to pore pressure increases. The increase in pore pressure was monitored through time-lapse shear-wave anisotropy analysis. Since the onset of the observation period in 2010, the orientation of the fast shear wave has rotated 90°, indicating a change associated with critical pore pressure build up. The time delay between fast and slow shear wave arrivals has increased, indicating a corresponding increase in anisotropy induced by pore pressure rise. In-situ near-basement fluid pressure measurements corroborate the continuous pore pressure increase revealed by the shear-wave anisotropy analysis over the earthquake monitoring period. This research is the first to identify a change in pore fluid pressure in the basement using seismological data and it was recently published in the AAAS journal Science Advances (Nolte et al., 2017). The shear-wave splitting analysis is a novel application of the technique, which can be used in other regions to identify an increase in pore pressure. This increasing pore fluid pressure has become more regionally extensive as earthquakes are occurring in southern Kansas, where they previously were absent. These monitoring techniques and analyses provide new insight into mitigating induced seismicity’s impact on society.

02 PETROLEUM↗
DOE OSTI2021

Hierarchical Fractional Advection-Dispersion Equation (FADE) to Quantify Anomalous Transport in River Corridor over a Broad Spectrum of Scales: Theory and Applications

Fractional calculus-based differential equations were found by previous studies to be promising tools in simulating local-scale anomalous diffusion for pollutants transport in natural geological media (geomedia), but efficient models are still needed for simulating anomalous transport over a broad spectrum of scales. This study proposed a hierarchical framework of fractional advection-dispersion equations (FADEs) for modeling pollutants moving in the river corridor at a full spectrum of scales. Applications showed that the fixed-index FADE could model bed sediment and manganese transport in streams at the geomorphologic unit scale, whereas the variable-index FADE well fitted bedload snapshots at the reach scale with spatially varying indices. Further analyses revealed that the selection of the FADEs depended on the scale, type of the geomedium (i.e., riverbed, aquifer, or soil), and the type of available observation dataset (i.e., the tracer snapshot or breakthrough curve (BTC)). When the pollutant BTC was used, a single-index FADE with scale-dependent parameters could fit the data by upscaling anomalous transport without mapping the sub-grid, intermediate multi-index anomalous diffusion. Pollutant transport in geomedia, therefore, may exhibit complex anomalous scaling in space (and/or time), and the identification of the FADE’s index for the reach-scale anomalous transport, which links the geomorphologic unit and watershed scales, is the core for reliable applications of fractional calculus in hydrology.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Machine Learning Models for Mapping Groundwater Pollution Risk: Advancing Water Security and Sustainable Development Goals in Georgia, USA

The widespread use of pesticides, such as atrazine and malathion, in agricultural systems raises significant concerns regarding the contamination of groundwater, which serves as a critical resource for drinking water. This study applies machine learning techniques to predict the concentrations of atrazine and malathion in groundwater across Georgia, USA, using 2019 data. A Random Forest classifier was employed to integrate various environmental and demographic factors, including pesticide application rates, precipitation, lithology, and population density, to predict pesticide contamination in groundwater. The models demonstrated high training accuracies of 100% and moderate average testing accuracy of 55% for atrazine and 60% for malathion across five iterations. The low test accuracy of the model, ranging from 50% to 75%, is likely due to overfitting, which can be attributed to the small dataset size and the complex nature of pesticide-contamination patterns, making it challenging for the model to generalize to unseen data. Feature importance analysis revealed that average pesticide usage emerged as the most influential factor for atrazine, while aquifer lithology and precipitation played crucial roles in both models. These results provide valuable insights into the dynamics of pesticide contamination, highlighting areas at greater risk of contamination. The findings underscore the importance of integrating environmental, geological, and agricultural variables for more effective groundwater management and sustainable agricultural practices, contributing to the protection of water resources and public health.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Integration of seismic-pressure-petrophysics inversion of continuous active-seismic monitoring data for monitoring and quantifying CO 2 plume (Final Report)

The overall objective of this project is to develop and validate an integrated package of joint seismic-pressure-petrophysics inversion (jSPPI) of continuous active-source seismic monitoring dataset capable of providing real-time monitoring of CO 2 plume during geologic carbon sequestration (GCS). The three specific developments include: (a) the methodologies for fast seismic full waveform inversion of continuous active source seismic monitoring, (CASSM) datasets for simultaneously estimating velocity and attenuation, and with data assimilation; (b) joint Bayesian petrophysical inversion of seismic models and pressure data for providing and updating CO 2 saturation models; (c) the methods using multiple datasets including (Crainfield and Frio-II borehole) synthetic, laboratory, and field CASSM datasets. The outcomes of jSPPI include (a) a workflow for processing CASSM data, (b) Bayesian inversion algorithms using CASSM data and pressure response data, and (c) integration with data assimilation algorithms for continuously updating site-specific models used for prediction and reservoir management. The validation of joint FWI will be conducted using synthetic models based on the Cranfield and Frio experiments as well as field CASSM datasets collected as part of the Frio-II pilot injection. To quantify and map the mass and distribution of CO 2 (saturation), we will jointly invert velocity and attenuation measurements from the FWI with a Bayesian approach using a rock physics model for attenuation (e.g., White’s attenuation model with two selected patch sizes (White, 1976; Dutta and Seriff, 1979)). The Bayesian inversion will be applied to each time step in the CASSM survey in an updating scheme, which integrates with an ensemble of reservoir simulations at each step. A more complete experimental validation dataset will be collected as part of a mesoscale (2-3 m) gas-CO 2 injection experiment utilizing a higher frequency version of the CASSM system developed for laboratory studies; the integrated inversion will be demonstrated using this dataset which will provide both a dense geometry as well as more precise secondary confirmation measurements (e.g. saturation) typically not available in the field. The resulting real-time map of CO 2 saturation is able to provide a deeper scientific understanding of the complex, time-varying dynamics of subsurface fluid flow migration path as well as the rapid detection of CO 2 leakage hazards.

25 ENERGY STORAGE↗
DOE OSTI2018

Fallon FORGE: Seismic Reflection Profiles

Newly reprocessed Naval Air Station Fallon (1994) seismic lines: pre-stack depth migrations, with interpretations to support the Fallon FORGE (Phase 2B) 3D Geologic model. Data along seven profiles (>100 km of total profile length) through and adjacent to the Fallon site were re-processed. The most up-to-date, industry-tested seismic processing techniques were utilized to improve the signal strength and coherency in the sedimentary, volcanic, and Mesozoic crystalline basement sections, in conjunction with fault diffractions in order to improve the identification and definition of faults within the study area.

15 GEOTHERMAL ENERGY↗
DOE OSTI2026

Prairie State Generating Company Static and Dynamic Modeling

A reservoir modeling study was conducted to assess the feasibility of storing 162.5 million tonnes (8.125 million tonnes annually) of industrially sourced carbon dioxide (CO 2 ) in the St. Peter – Everton and Knox storage complexes at the Prairie State Generating Company’s (PSGC) site near Marissa, Washington County, Illinois (PSGC site). Two separate models were constructed for the St. Peter – Everton and the Knox storage complexes. The St. Peter and Everton sandstones are the target storage units for the St. Peter – Everton storage complex. The Knox Group formations are the target storage units for the Knox storage complex. The Maquoketa Shale is the primary confining unit for both storage complexes. The static reservoir models used as input to the dynamic reservoir models of the storage complexes were developed from interpretations of logs (e.g., gamma ray, resistivity, porosity, photoelectric, and sonic), structure surfaces using data from 36 wells, thickness maps, well test data, seismic data, and permeability data in PetrelTM.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2013

Investigation of uranium attenuation and release at column and pore scales in response to advective geochemical gradients

Column experiments were devised to investigate the role of changing fluid composition on mobility of uranium through a sequence of geologic media. Fluids and media were chosen to be relevant to the ground water plume emanating from the former S-3 ponds at the Oak Ridge Integrated Field Research Challenge (ORIFC) site. Synthetic ground waters were pumped upwards at 0.05 mL/minute for 21 days through layers of quartz sand alternating with layers of uncontaminated soil, quartz sand mixed with illite, quartz sand coated with iron oxides, and another soil layer. Increases in pH or concentration of phosphate, bicarbonate, or acetate were imposed on the influent solutions after each 7 pore volumes while uranium (as uranyl) remained constant at 0.1mM. A control column maintained the original synthetic groundwater composition with 0.1mM U. Pore water solutions were extracted to assess U retention and release in relation to the advective ligand or pH gradients. Following the column experiments, subsamples from each layer were characterized using microbeam X-ray absorption spectroscopy (XANES) in conjunction with X-ray fluorescence mapping and compared to sediment core samples from the ORIFC, at SSRL Beam Line 2-3.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2016

A new algorithm for three‐dimensional joint inversion of body wave and surface wave data and its application to the Southern California plate boundary region

Abstract We introduce a new algorithm for joint inversion of body wave and surface wave data to get better 3‐D P wave ( V p ) and S wave ( V s ) velocity models by taking advantage of the complementary strengths of each data set. Our joint inversion algorithm uses a one‐step inversion of surface wave traveltime measurements at different periods for 3‐D V s and V p models without constructing the intermediate phase or group velocity maps. This allows a more straightforward modeling of surface wave traveltime data with the body wave arrival times. We take into consideration the sensitivity of surface wave data with respect to V p in addition to its large sensitivity to V s , which means both models are constrained by two different data types. The method is applied to determine 3‐D crustal V p and V s models using body wave and Rayleigh wave data in the Southern California plate boundary region, which has previously been studied with both double‐difference tomography method using body wave arrival times and ambient noise tomography method with Rayleigh and Love wave group velocity dispersion measurements. Our approach creates self‐consistent and unique models with no prominent gaps, with Rayleigh wave data resolving shallow and large‐scale features and body wave data constraining relatively deeper structures where their ray coverage is good. The velocity model from the joint inversion is consistent with local geological structures and produces better fits to observed seismic waveforms than the current Southern California Earthquake Center (SCEC) model.

Fang, Hongjian↗
DOE OSTI2023

Toward Improved Regional Hydrological Model Performance Using State-Of-The-Science Data-Informed Soil Parameters

Accurate soil moisture and streamflow data are an aspirational need of many hydrologically relevant fields. Model simulated soil moisture and streamflow hold promise but models require validation prior to application. Calibration methods are commonly used to improve model fidelity but misrepresentation of the true dynamics remains a challenge. In this study, we leverage soil parameter estimates from the Soil Survey Geographic (SSURGO) database and the probability mapping of SSURGO (POLARIS) to improve the representation of hydrologic processes in the Weather Research and Forecasting Hydrological modeling system (WRF-Hydro) over a central California domain. Our results show WRF-Hydro soil moisture exhibits increased correlation coefficients ( r ), reduced biases, and increased Kling-Gupta Efficiencies (KGEs) across seven in situ soil moisture observing stations after updating the model's soil parameters according to POLARIS. Compared to four well-established soil moisture data sets including Soil Moisture Active Passive data and three Phase 2 North American Land Data Assimilation System land surface models, our POLARIS-adjusted WRF-Hydro simulations produce the highest mean KGE (0.69) across the seven stations. More importantly, WRF-Hydro streamflow fidelity also increases, especially in the case where the model domain is set up with SSURGO-informed total soil thickness. The magnitude and timing of peak flow events are better captured, r increases across nine United States Geological Survey stream gages, and the mean KGE across seven of the nine gages increases from 0.12 to 0.66. Our pre-calibration parameter estimate approach, which is transferable to other spatially distributed hydrological models, can substantially improve a model's performance, helping reduce calibration efforts and computational costs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Geothermal Exploration of Mount Baker Hot Springs Through Ground-Based Magnetic and Gravity Surveys

The Washington State Division of Geology and Earth Resources (WDGER), together with the USGS, Altarock Energy Inc., Colorado College, Temple University, Innovate Geothermal Ltd., University of British Columbia, and Western Washington University(WWU), conducted phase two of a geothermal play fairway exploration at three plays (Mount Baker, Mount St. Helens seismic zone, and Wind River valley) along the Washington Cascade Range. The Mount Baker (MB) play is an ideal location for further investigation as a potential geothermal resource due to the presence of thermal features, young volcanic centers, existing geothermal leases, accessibility, and its proximity to existing transmission lines at the Baker Lake Dam. The goal of phase II exploration at all three plays is to collect geologic and geophysical data to help constrain spatial changes in subsurface lithology and structures that may be associated with the geothermal systems. Furthermore, this paper discusses the geothermal potential of the southeast flank of Mount Baker, WA by interpreting newly collected ground-based gravity and magnetic data. The Mount Baker play was surveyed from July-September, 2016 by the USGS, WDGER, and WWU. Results of a detailed magnetic and gravity surveys of the region surrounding the Baker hot springs, the primary surface expression of the hydrothermal system, are presented here. Gravity data were collected at 495 stations in the study area amounting to an area of roughly 150 square km. Approximately 93 km of line magnetometer data were collected within the study area, over a dozen magnetic susceptibility measurements were taken at 50 locations and about 50 hand samples were collected for density measurements. Preliminary data and modeling suggest pronounced magnetic anomalies trending NE/SW, generally in line with previously mapped lineaments identified in LiDAR. We present a series of preliminary two-dimensional forward models of the gravity and magnetic data as well as a 3D inversion model of magnetic susceptibility and interpret these models in terms of the structural controls influencing the Baker hot springs and hydrothermal system.

Geophysics↗
DOE OSTI2026

Qualitative Risk Assessment of Legacy Wells within the Estimated Prairie State Generating Company Area of Review

This report details the digitization of a legacy wellbore database, including data processing assumptions, parameter estimation, and risk assessment methodology. The database, comprising 6,454 documents, was provided by ISGS. It includes valuable data from the Prairie State Generating Company (PSGC) and One Earth Energy (OEE) sites of the CarbonSAFE Phase III – Illinois Storage Corridor project. The report focuses on wells within a 15-mile radius from the Lively Grove #1 (LG#1) well at PSGC site, evaluating subsurface conditions and potential risks. A total of 4,386 wellbores within 15 miles of the LG#1 well were filtered based on depth and formation codes. LG#1 is the stratigraphic well at the PSGC site drilled in 2021. Ninety-four (94) wells penetrating the Maquoketa Shale Group (the primary confining unit) within the estimated area-of-review (AoR) for the PSGC site were evaluated using a qualitative risk assessment (QRA) methodology. The QRA developed by Arbad et al. 2022 focuses on legacy wells within the AoR and categorizes them based on well construction details. The QRA identifies wells that need immediate attention by categorizing them based on penetration depth and protection. Wells within the AoR were categorized into nine groups based on penetrations and protections. These categories range from Type 1 wells, with no documentation, to Type 9 wells, which do not penetrate the primary confining unit or storage reservoir (unit). Well accessibility within the AoR varies based on well status, including Dry & Abandoned (DA), Plugged & Abandoned (PA), Injection (INJ), Oil/Gas Producing (PROD), and Observation (Obs) wells. Accessibility levels were determined by well construction, with DA wells being the least accessible and Observation wells the most accessible, impacting gas leakage detection possibilities. Remedial action priority of wells decreases from Type 1 to Type 9 wells. Type 1 to Type 6 wells with status DA and PA require immediate attention, while Type 7 and Type 8 wells are low priority. A risk matrix used to prioritize corrective actions for legacy wells is proposed to categorize wells within an AoR based on penetrations, protections, and accessibility. The methodology involves data acquisition, well categorization into nine types, and determining CO 2 leakage pathways using well schematics and geospatial mapping. This approach is particularly useful for managing the integrity of legacy wells throughout the lifecycle of a Carbon Capture and Storage (CCS) project. A qualitative risk assessment of 94 wells within the AoR of the PSGC site identified 54 wells with high priority for corrective action due to penetration of the primary containment seal. The assessment utilizes color-coded maps to categorize well types and prioritize corrective actions, providing a comprehensive analysis. Schematics of wells penetrating the primary confining unit were drawn, and leakage pathways were identified. Details of all wells penetrating the confining zone are provided in the appendix, including information on well types, plugging, and casing status.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2021

Carbon content of soil fractions varies with season, rainfall, and soil fertility across a lowland tropical moist forest gradient

Tropical forests contain some of the largest soil carbon (C) stocks on Earth, making them broadly relevant to terrestrial-climate feedbacks, yet our understanding of how their soil organic C (SOC) fractions vary over space and time is limited. We studied effects of season, fertility, and mean annual precipitation (MAP) on the C contents of soil fractions across 14 lowland forests in Panama. We measured free-debris, occluded-debris, and mineral-associated SOC fractions, as well as soluble C associated with each fraction. We hypothesized that mineral-associated SOC would be greatest in infertile, strongly weathered soils with large reactive mineral contents. We also hypothesized that the debris SOC fractions would accumulate during the dry season, reflecting seasonal increases in litterfall. To address this, we compared soil fractions in wet and dry seasons from fertile and infertile soils across a range of 1809 – 2864 mm MAP. The C content (mg C / g soil) of all soil fractions varied with fertility and MAP: specifically, free-debris SOC was greatest in wet, high-fertility soils, and occluded-debris SOC was greater in high-fertility than low-fertility soils. The mineral-associated SOC fraction, which contained the majority of bulk soil C, showed increasing C content with greater MAP in infertile sites, presumably driving similar spatial patterns in the bulk soil. Only the free-debris SOC fraction showed strong seasonal variation, increasing in mass during the dry season. Nitrogen behaved similarly to C. In summary, soil C contents increased with MAP in infertile sites but not fertile sites, driven by the mineral-associated SOC fraction. The dry season had greater free-debris SOC, but this seasonal trend was not apparent in bulk soil C, likely because of the small size of the free-debris SOC fraction. Furthermore, changes in the quantity and seasonality of precipitation, which are projected for tropical forests, might shift spatial and temporal patterns of soil C storage, which would in turn influence forest-climate feedbacks for this C-rich biome.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

Water cycling between ocean and mantle: Super-earths need not be waterworlds

Large terrestrial planets are expected to have muted topography and deep oceans, implying that most super-Earths should be entirely covered in water, so-called waterworlds. This is important because waterworlds lack a silicate weathering thermostat so their climate is predicted to be less stable than that of planets with exposed continents. In other words, the continuously habitable zone for waterworlds is much narrower than for Earth-like planets. A planet's water is partitioned, however, between a surface reservoir, the ocean, and an interior reservoir, the mantle. Plate tectonics transports water between these reservoirs on geological timescales. Degassing of melt at mid-ocean ridges and serpentinization of oceanic crust depend negatively and positively on seafloor pressure, respectively, providing a stabilizing feedback on long-term ocean volume. Motivated by Earth's approximately steady-state deep water cycle, we develop a two-box model of the hydrosphere and derive steady-state solutions to the water partitioning on terrestrial planets. Critically, hydrostatic seafloor pressure is proportional to surface gravity, so super-Earths with a deep water cycle will tend to store more water in the mantle. We conclude that a tectonically active terrestrial planet of any mass can maintain exposed continents if its water mass fraction is less than ∼0.2%, dramatically increasing the odds that super-Earths are habitable. The greatest source of uncertainty in our study is Earth's current mantle water inventory: the greater its value, the more robust planets are to inundation. Lastly, we discuss how future missions can test our hypothesis by mapping the oceans and continents of massive terrestrial planets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Basin-Scale Structural Features Database

The Basin-Scale Structural Features database provides spatial datasets of faults, fractures, folds, and earthquakes compiled from public, authoritative sources (e.g., U.S. Geological Survey and State Geological Surveys) and aggregated into derivative forms to support subsurface assessments. Recognizing that characterizing basin-scale structural features requires interpreting data that are often ambiguous or lack key information, the source data were evaluated using a knowledge-data framework and geospatial fuzzy logic method (Justman et al., 2020) to represent both measured (observed) and predicted (inferred or potential) structural features as derivative datasets. This workflow employs conceptual models for known structural features and predicted structural features, incorporating geospatial data to estimate potential, even with limited data. The aim is to aid and support an understanding of basin-scale features and identify potential gaps in data and knowledge. As of 4/30/2025, the database includes resources for nine sedimentary basins: Appalachian, Denver, U.S. Gulf Coast, Illinois, Michigan, Permian, Sacramento, San Joquin and Williston. The database is organized by basin and then data category: 1) Faults, fractures, folds, 2) Earthquakes, 3) Topographic, 4) Structural contours and isopachs, 5) Geophysical, and 6) Structural feature density assessment maps.

basin scale↗