Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “geologic characterization”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

ED-cPSD: Fast Phase-Size Distribution via Sequential Erosion-Dilation

The Erosion-Dilation continuous Phase-Size Distribution, ED-cPSD, is an application for calculating continuous pore and particle-size distribution from digital reconstructions and/or image-based structural data. It is based on the erosion-dilation continuous phase-size distribution method. A continuous size distribution is a measure of the probability density of finding a particle or pore of a certain size. These distributions are of interest in any field of study involving porous media, including but not limited to electrochemistry, petroleum engineering, geology, and food science. The algorithm behind the software provides a computationally efficient way to calculate phase-size distributions for large domains. For a 3D battery electrode reconstruction with 1.3 x 10 8 voxels, the particle size distribution is derived in under 2 min on a desktop, while also retaining flexibility and computational efficiency for HPC-scale multi-threading. The software can handle structures with over 10 9 voxels. The algorithm is roughly 280 times faster than a previous version on the same task.

Characterization↗
DOE OSTI2024

Idaho National Laboratory Integrated Multisite SSHAC Level 3: Probabilistic Volcanic Hazards Assessment

The Idaho National Laboratory (INL) resides on the eastern Snake River Plain (ESRP), part of the Snake River Plain (SRP) with a complex origin and geologic history of volcanism. Much of the Quaternary (last 2.58 million years) volcanism within 400 km of INL is genetically associated with a major thermal anomaly referred to as the Yellowstone hotspot, which is currently located more than 180 km northeast of INL. Near INL, local volcanic sources include silicic domes near its southern border, and numerous dike-fed basaltic vents of the ESRP, some of which are near or within the INL boundaries. Quantitative probabilistic assessments of screened-in volcanic hazardous phenomena for nine different facility complexes at the INL are presented for a Senior Seismic Hazard Analysis Committee (SSHAC) Level 3 (SL3) study. The comprehensive, integrated multisite SSHAC study consists of a single regional Probabilistic Volcanic Hazards Assessment (PVHA) that pertains to all nine INL facility complexes, with site-specific information developed at each respective area of interest (AOI), referred to as the "INL facility AOI" (Figure ES-1). Hazard products are generated for specified INL facility AOIs for use by multiple stakeholders from different agencies in risk-informed decision-making regarding site selection, operations, and design of nuclear facilities at INL, consistent with U.S. Department of Energy (DOE) and U.S. Nuclear Regulatory Commission (NRC) regulatory guidance. The study also serves as the basis for future periodic safety assessments required for existing DOE facilities, such as 10-year evaluations of natural-phenomena hazards. Elements of the INL PVHA including initial characterization, screening, quantitative assessments of eruption and hazard potential, and consideration of facility needs are conducted using the three phases of the SSHAC process: evaluation, integration, and documentation. As per regulatory guidance, the SSHAC framework provided the necessary processes and procedures for the PVHA Technical Integration (TI) team, at three workshops, five formal working meetings and many TI team remote meetings, to conduct initial characterization, screen the volcanic hazards, create PVHA model inputs, exercise those models in the PVHA, consider facility-specific volcanic hazard needs, and perform final hazard calculations. The Participatory Peer Review Panel (PPRP) provided independent oversight and performed process and technical reviews of the PVHA throughout its duration. The study included an extensive New Data Collection and Analyses (NDCA) program developed by the PVHA TI team to reduce uncertainties in hazard-significant elements in the PVHA model. NDCA activities generated 20 reports providing important contributory datasets and results to the project database for characterizing the ESRP. For example, a report compiling the dimensions of ESRP shield volcanoes and lava fields was used to construct volcanic footprints (areas of impact), was compared with data from INL subsurface cores, and was used to validate the results of lava-flow inundation modeling on the contemporary terrain. Another example is the acquisition of aeromagnetic data over INL and its surrounding area, with maps of buried magmatic features (e.g., subsurface volcanoes and swarms of feeder dikes) that informed the PVHA conceptual model of volcanism. The SSHAC evaluation process was used to conduct all elements of the PVHA including initial characterization and screening. Existing data and NDCA activities provided the foundation to develop the tectonomagmatic conceptual model of volcanism for the region of geographical interest in the SRP and Yellowstone hotspot volcanic system, and for considering volcanoes in the western US. Considering Quaternary volcanic sources active during this period, the screening approach identified and evaluated magma compositions, types of eruptive and intrusive phenomena, types of hazardous phenomena, and proximity of sources to INL. The PVHA TI team evaluated 18 types of magmatic sources in terms of 20 potentially hazardous phenomena, resulting in 360 screening decisions. The screening process resulted in 140 screened-in hazardous volcanic phenomena for Quaternary volcanic sources 1) proximal to INL facility complexes in the ESRP (64 basaltic and 54 silicic), 2) regional sources associated with the Yellowstone caldera system and Blackfoot Reservoir volcanic field (7), and 3) more distal sources from thirteen Cascade volcanoes and two volcanoes at Long Valley caldera (CA).

58 GEOSCIENCES↗
DOE OSTI2022

Validation of earthquake ground-motion models in southern California, USA, using precariously balanced rocks

Accurate estimates of earthquake ground shaking rely on uncertain ground-motion models derived from limited instrumental recordings of historical earthquakes. A critical issue is that there is currently no method to empirically validate the resultant ground-motion estimates of these models at the timescale of rare, large earthquakes; this lack of validation causes great uncertainty in ground-motion estimates. Here, we address this issue and validate ground-motion estimates for southern California utilizing the unexceeded ground motions recorded by 20 precariously balanced rocks. We used cosmogenic 10Be exposure dating to model the age of the precariously balanced rocks, which ranged from ca. 1 ka to ca. 50 ka, and calculated their probability of toppling at different ground-motion levels. With this rock data, we then validated the earthquake ground motions estimated by the Uniform California Earthquake Rupture Forecast, Version 3 (UCERF3) seismic-source characterization and the Next Generation Attenuation (NGA)-West2 ground-motion models. We found that no ground-motion model estimated levels of earthquake ground shaking consistent with the observed continued existence of all 20 precariously balanced rocks. The ground-motion model I14 estimated ground-motion levels that were inconsistent with the most rocks; therefore, I14 was invalidated and removed. At a 2475 year mean return period, the removal of this invalid ground-motion model resulted in a 2−7% reduction in the mean and a 10−36% reduction in the 5th−95th fractile uncertainty of the ground-motion estimates. Our findings demonstrate the value of empirical data from precariously balanced rocks as a validation tool for removing invalid ground-motion models and, in turn, reducing the uncertainty in earthquake ground-motion estimates.

Geology↗
DOE OSTI2023

Connecting Large‐Scale Meteorological Patterns to Extratropical Cyclones in CMIP6 Climate Models Using Self‐Organizing Maps

Abstract Extratropical cyclones (ETCs) are responsible for the majority of cool‐season extreme events in the northeastern United States (NEUS), often leading to high‐impact weather conditions that can have wide‐ranging socioeconomic impacts. Evaluating the ability of climate models to adequately simulate ETC dynamics is essential for improving model performance and increasing confidence in future projections used by stakeholders and policymakers. ETCs are traditionally studied using techniques such as case studies and synoptic typing, however, these approaches can be time‐consuming, require subjective analysis, and do not necessarily identify the coincident large‐scale meteorological patterns (LSMPs). Here, we apply self‐organizing maps (SOMs) as an automated machine‐learning approach to characterize the LSMPs and associated frequency and intensity of discrete ETC events over NEUS. The dominant patterns of geopotential height variability are identified through SOM analysis of five reanalysis products during the last four decades. ETC events are tracked using TempestExtremes and are integrated with SOMs to classify the accumulated cyclone activity (ACA) associated with each pattern. We then evaluate the skill of CMIP6 historical experiments in simulating the LSMPs and ETC events identified in the SOM. Our results identify a robust bias toward more zonal patterns, with models struggling to reproduce the more amplified patterns typically associated with the highest cyclone activity. While model resolution has some impact on simulation credibility, model configuration appears to be more important in LSMP representation. The vast majority of CMIP6 models produce too few ETCs, although model errors are distributed around historical reanalyses when ACA is normalized by storm frequency.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Contribution of El Niño Southern Oscillation (ENSO) Diversity to Low‐Frequency Changes in ENSO Variance

Abstract El Niño Southern Oscillation (ENSO) diversity is characterized based on the longitudinal location of maximum sea surface temperature anomalies (SSTA) and amplitude in the tropical Pacific, as Central Pacific events are typically weaker than Eastern Pacific events. SSTA pattern and intensity undergo low‐frequency modulations, affecting ENSO prediction skill and remote impacts, and resulting in low‐frequency changes in ENSO variance. Yet, how different ENSO types contribute to these decadal variance changes remains unclear. Here, we decompose the low‐frequency changes of ENSO variance into contributions from ENSO diversity categories. We propose a fuzzy clustering of monthly SSTA to allow for non‐binary event category memberships, where each event can belong to different clusters. Our approach identifies two La Niña and three El Niño categories and shows that the major shift of ENSO variance in the mid‐1970s was associated with an increasing likelihood of strong La Niña and extreme El Niño events.

Geology↗
DOE OSTI2020

Zero-Valent Iron Permeable Reactive Barrier to Remediate Volatile Organic Compounds in Groundwater - 20209

The US DOE-SRS, the US EPA, and the South Carolina Department of Health and Environmental Control determined it was appropriate to perform a non-time critical removal action at the P-Area Groundwater Operable Unit at the SRS to reduce the mass and downgradient transport of trichloroethylene in the P-Area groundwater plume. Contaminated groundwater discharges to a nearby stream, Steel Creek, within the SRS boundaries, resulting in trichloroethylene concentrations above the maximum contaminant level. Impact to surface water is limited in areal extent and supported by recently collected characterization data. The P-Area Groundwater Operable Unit encompasses the groundwater beneath an industrial area within SRS, P Area, where the P-Reactor once operated. The boundaries of the P-Area Groundwater Operable Unit extend northwest to Steel Creek, northeast toward PAR Pond, and southeast to Meyers Branch. Groundwater in the Upper Three Runs Aquifer of the P-Area Groundwater Operable Unit has been impacted by reactor and facility operations between 1954 and 1991, including tritium and volatile organic compounds. The P-Area surface units contributing to groundwater contamination were remediated as part of the P-Area Operable Unit in 2011. The P-Reactor closure is one of the first of its kind in the DOE Complex and is one of only a few full-sized production reactors in the US to undergo completion of final closure activities. The nature and extent of groundwater contamination was determined using a variety of investigative approaches such as groundwater monitoring wells, direct-push technology, and surface water samples. Groundwater contamination associated with trichloroethylene is primarily exhibited in a narrow plume that extends from the source area at P-Reactor and west to Steel Creek. Maximum contaminant level exceedances in groundwater occur over an area of ∼6.9 hectares for trichloroethylene with concentrations as high as 7.7 milligrams per liter. To the west of the P-Area facility area, the trichloroethylene groundwater plume is controlled by a buried geologic feature, assumed to be an old stream bed, that further narrows the groundwater plume in what has been designated as the 'neck area.' This narrowing of the groundwater plume provides an ideal location for a treatment barrier. The non-time critical removal action alternative chosen is to install a zero-valent iron permeable reactive barrier within the neck area of the trichloroethylene groundwater plume, perpendicular to groundwater flow direction. This technology will provide a treatment barrier that will reduce trichloroethylene groundwater concentrations by 90% and has an anticipated useful life of at least 25 years. A pre-design investigation was performed in the neck area to confirm site lithology, hydrogeology, geochemistry, and extent of trichloroethylene contamination prior to a final design. A treatability study, conducted as part of the pre-design investigation, indicated that the subsurface and groundwater in the PArea Groundwater Operable Unit is compatible with the zero-valent iron and will not lead to excessive buildup from mineralization/precipitation or biofouling. Probabilistic modeling was conducted using field and laboratory data to determine the expected performance of the zero-valent iron permeable reactive barrier. The model simulations indicated that a 3.81-centimeter thick barrier would provide greater than 90% reduction of trichloroethylene groundwater concentration. The final design of the zero-valent iron permeable reactive barrier is a barrier that will extend 80.5 linear meters in a 'zigzag' orientation to best transect the trichloroethylene plume and account for varying groundwater flow. The barrier will be installed from 13.7 meters below ground surface to 41.1 meters below ground surface for 65.8 linear meters and from 13.7 meters below ground surface to 36.6 meters below ground surface for 14.6 linear meters, the base of which is 'keyed' into a low permeability zone. The barrier is designed to a thickness of 10.2 centimeters, which was determined to reduce the trichloroethylene groundwater concentrations by greater than 90% with a safety factor of 2.67. A total of approximately 689 metric tons of zero-valent iron will be injected through 22 injection wells spaced 3.66 meters apart, using guar to suspend the zero-valent iron. Zero-valent iron permeable reactive barrier construction will be monitored through 23 installed resistivity receivers offset 7.32 meters from the zerovalent iron permeable reactive barrier. The zero-valent iron will be energized with a low-voltage 100 Hertz signal during injection and will be monitored using the resistivity receivers to ensure complete coalescence of the zero-valent iron permeable reactive barrier. The zero-valent iron was sized to have a hydraulic conductivity greater than the natural subsurface, thus promoting groundwater flow through the barrier. As contaminated groundwater contacts the zero-valent iron, volatile organic compounds, including trichloroethylene, are immediately degraded to harmless compounds such as ethylene. The performance of the zero-valent iron permeable reactive barrier will be monitored using three upgradient monitoring well clusters, six downgradient monitoring well clusters, and four in-wall monitoring wells. The in-wall monitoring wells will indicate immediate reduction of trichloroethylene mass in the groundwater and allow analyses of the zero-valent iron permeable reactive barrier health. (authors)

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

Interpretation of Hydraulic Fractures Based on Microseismic Response in the Marcellus Shale, Monongalia County, West Virginia, USA: Implications for Shale Gas Production

Summary Hydraulic fracturing is critical for extracting shale gas in the subsurface. The treatment technique of multistage hydraulic fracturing is widely used to maximize production. In multistage hydraulic fracturing, not all pumping stages make the same contribution for production, although the designed stimulation process is almost same in every pumping stage. In this study, we characterize the microseismic responses of hydraulic fracturing. Significant variations of microseismic characteristics are observed among the different pumping stages. Pre-existing natural fractures are examined along the horizontal well in selected stages. Combined with the Mohr’s circle analysis, results of the fracture study show that induced hydraulic fractures can be captured by pre-existing natural fractures. Induced hydraulic fractures are simulated by the unconventional fracture model (UFM), and the result reveals that the stress-shadow effect diverts the direction of hydraulic fractures. The diverted hydraulic fracture affects the development of the hydraulic fracture network, which has influence on production. The treatment of two-step pumping is investigated by application of hydraulic diffusivity. Hydraulic fracturing performance and production could be optimized by the treatment of two-step pumping in a single stage. The second pumping creates new fractures and fills the fractures with additional proppants to maintain production for a long duration.

Energy & Fuels↗
DOE OSTI2024

Growth and Stability of Stratiform Carrollite (CuCo2S4) in the Tenke-Fungurume Ore District, Central African Copperbelt

Abstract Carrollite (CuCo2S4) is the main ore mineral of the critical battery metal cobalt, yet, surprisingly, detailed characterization of its thermodynamic properties and of its trace element contents remain sparse in the scientific literature. To fill these knowledge gaps, we generated the first thermodynamic data set for carrollite using well-characterized samples obtained from the stratiform sulfide mineralization style at Fungurume 88 deposit in the Democratic Republic of Congo. Detailed mineralogical investigation of these same samples suggests that carrollite has limited ability to incorporate other ‘sweetener’ elements (except Cu and Ni) by substitution mechanisms. Modeling using CHNOSZ 2.0.0. software and an integrated Cu-Co-S-O thermodynamic data set provides insights into the stability of carrollite relative to other Co-bearing minerals under physicochemical conditions relevant to the Central African Copperbelt. Further exploration of these models indicates that dropping redox potential in the mineralizing system is insufficient to explain the mineral-scale and deposit-scale zonation patterns commonly observed in Cu-Co deposits of the Central African Copperbelt. Instead, factors such as increasing pH, decreasing T, decreasing sulfur activity, or some combination of changing physicochemical parameters may need to be invoked to explain these zonation trends. At the Fungurume 88 deposit, the absence of linnaeite (Co3S4) suggests that ore precipitation occurred at temperatures cooler than 211 °C. From our pH-logfO2 diagrams, and relative to higher temperatures, such a lower T regime is marked by a greater degree of offset between the modeled Cu and Co solubility contours, thus facilitating a more differentiated and zoned distribution of these two metals at the deposit scale. Our models are offered to the geoscience community as exceptionally useful tools for modeling and visualizing the Cu-Co-S-O systematics and may be applied to a variety of geological and geometallurgical questions in the Central African Copperbelt and other deposits elsewhere.

Mineralogy↗
DOE OSTI2025

Geochemical Impact of Acid Spearhead and Slickwater Stimulation on Wolfcamp Shale from the Hydraulic Fracturing Test Site

The Hydraulic Fracturing Test Site 1 (HFTS-1) was a field study performed in the Wolfcamp Formation in the West Texas Permian (Midland) Basin, USA, with a focus on improving the efficiency of hydraulic fracturing. Investigating site-specific rock-fluid geochemical interactions during hydraulic fracturing is an important step to understanding the impact on formation shale porosity, permeability, and long-term shale gas production. During field operations in this region, hydraulic fracturing fluid (HFF) injection usually starts with a concentrated acid spearhead for rapid rock dissolution, followed by the injection of near-neutral pH slickwater containing chemicals and proppants. A multistep sequential injection approach was used to investigate different stages of rock-fluid interactions. The carbonate content in the host rock is important when acid spearhead is considered, as carbonate mineral dissolution is rapid and can result in porosity and permeability changes in the shale matrix. Here, in this study, we designed flow-through experiments using fractured carbonate-rich and clay-rich Wolfcamp shale cores with (1) a short-time acid soaking step and (2) a long-term slickwater flow-through step to simulate the injection method used at HFTS-1. The fluid chemistry was analyzed. A thorough mineralogical progression [e.g., Calcium (Ca) dissolution and iron (Fe) redox progression] in the cores during HFF injection was also characterized and imaged by synchrotron microprobe. Reactive transport modeling was performed based on the experimental setup. The results showed that the acid spearhead is a crucial step in creating a reaction front by mineral dissolution, especially in carbonate-rich shales. A slight layer of ferrihydrite precipitated during the slickwater flow-through period. This study provides insights into potential geochemical impact due to hydraulic fracturing operations in the Permian Basin.

clastic rock↗
DOE OSTI2021

Newly acquired core enhances geologic understanding of the northern Paradox Basin Cane Creek play, southeastern Utah

The Cane Creek interval is an emerging unconventional heterolithic play within the Pennsylvanian Paradox Formation of the northern Paradox Basin, southeastern Utah. Structural controls (basement uplifts, salt movement, faults, fractures) likely control production from the Cane Creek, where total oil production of ~10 MMBO is just a fraction of current estimates of the undiscovered resource (up to 1.2 BBOE). Oil production has been most successful from unstimulated horizontal wells in the Cane Creek Unit in the central play area, whereas areas to the north (Greentown and Gunnison Valley) have to date proven to be largely unproductive from vertical and horizontal wells drilled without the benefit of 3D seismic or the deployment of modern hydraulic fracturing techniques. However, the limited historical drilling has shown significant promise (e.g., significant shows). The lack of core data in the northern part of the play has impeded thorough assessment of reservoir quality, character, and lithological facies heterogeneity. The Cane Creek is informally divided into three distinct zones. The upper A and basal C zones are composed of fabric destructive and wavy bedded anhydrite, dolomitic mudstone, and organic-rich algal laminated mudstones (source rock). The middle B zone is generally wave rippled and burrowed, low-permeability sandstone-siltstone (reservoir) with subordinate displacive anhydrite. Recent potash cores collected from the Salt Valley anticline near Crescent Junction, Utah (but east of the main play area), show thicker reservoir sandstone/siltstone packages (~40 ft) in the B zone, and parts of the C zone, in the northern part of the Paradox Basin in comparison to the central (~35 ft) and southern (~30 ft) areas. Source rock analyses from northern core/cuttings also indicate that it is well positioned within the oil window (VRo ~0.8), with thin organic-rich mudstone beds containing up to 20 wt% TOC. Because the siliciclastic reservoir rocks have low permeabilities (0.009–0.202 mD) and variable porosities (6–17%) due to dolomite-anhydrite cement, naturally occurring and possibly stimulated fractures may be essential for hydrocarbon recovery. To better understand the northern part of the Cane Creek play, a new core was obtained from the research stratigraphic well State 16-2 (January 2021). Details from this new core will help characterize reservoir quality via core analyses, petrography, fracture orientations, geomechanics, and correlation of facies and depositional trends. Additionally, an evaluation of matrix porosity/permeability will be conducted to determine its contribution to productivity.

02 PETROLEUM↗
DOE OSTI2025

UV-Vis spectrophotometric determination of rare earth elements (REE) speciation at near-neutral to alkaline pH. Part II: hydrolysis of Er from 25 to 75 °C

Aqueous speciation of rare earth elements (REE) controls their mobilization, fractionation, and enrichment in the natural waters. Geochemical modeling of their speciation is key to improve our understanding of the formation of economic mineral deposits, for developing mineral separation and mine tailing recovery technologies, and for characterizing the geochemistry of thermal water. However, our ability to predict the fate of REE in a wide pH and temperature range is limited by the scarcity of thermodynamic data for the REE hydroxyl complexes. In part I of this study (H. J. Han and A. P. Gysi, Dalton Trans., 2024, 53, 13129–13141), the optical properties of m-cresol purple (mCP) were determined using UV-Vis spectrophotometry between 25 and 75 °C in order to develop a method for deriving the hydrolysis constants of erbium (Er). Here, UV-Vis spectrophotometry experiments were conducted as a function of temperature between 35 and 75 °C to determine the hydrolysis of Er in near-neutral to alkaline solutions using mCP as an in situ pH color indicator. Here, the experiments were conducted with Er concentrations from 0 to ~0.253 mmol kg –1 in low ionic strength solutions (≤0.001 mol kg –1 ).

58 GEOSCIENCES↗
DOE OSTI2023

Microbial Models for Simulating Soil Carbon Dynamics: A Review

Abstract Soils store the largest amount of carbon (C) in the biosphere, and the C pool in soil is critical to the global C balance. Numerous microbial models have been developed over the last few decades to represent microbial processes that regulate the responses of soil organic carbon (SOC) to climate change. However, the representation of microbial processes varies, and how microbial processes are incorporated into SOC models has not been well explored. Here, we reviewed 71 microbial models to characterize the microbial processes incorporated into SOC models and analyzed variations in mechanistic complexity. We revealed that (a) four processes (microbial‐mediated decomposition, mineral interaction, microbial necromass recycling, and active and dormant microbial dynamics) are commonly incorporated in microbial models, (b) ∼48% of models simulate only one microbial process (i.e., microbial‐mediated decomposition) and 35% of models simulate two microbial processes: for example, microbial‐mediated decomposition and mineral interaction, (c) more than 80% microbial models use nonlinear equations, such as forward Michaelis‐Menten kinetics, to represent SOC decomposition, (d) the concept of persistence of SOC due to its intrinsic properties has been replaced by organo‐mineral interaction (∼39% of microbial models) that protects SOC from decomposition, and (e) various temperature and moisture modifiers and pH effects have been used to explain the environmental effect on microbial processes. In the future, to realistically incorporate microbial processes into Earth System Models, it is imperative to identify experimental evidence on rate limitation processes and firmly ground model structure on the field and laboratory data.

Environmental Sciences & Ecology↗
DOE OSTI2020

Sedimentary Geothermal Resources in Nevada, Utah, Colorado, and Texas

The objectives of this project were to (1) perform a literature review of sedimentary geothermal resources, (2) identify data sources and develop data-collection methodologies that characterize selected resources, (3) screen sedimentary basins and formations for sedimentary geothermal potential, and (4) evaluate the technical feasibility of one or more selected locations. Numerous publications have characterized geothermal resources within sedimentary basins. A literature search reviewed publications describing resources located in Colorado, Louisiana, Nevada, Texas, Utah, and Wyoming. The most attractive resources have high temperature gradients, low drilling costs, and reservoir permeabilities greater than 10 millidarcies (mD). Prospects in Colorado, Nevada, Texas, and Utah exhibit attractive characteristics and were chosen for further analysis. Sedimentary resources in Nevada and Utah are most attractive, followed by tested resources in Texas and untested resources in Colorado. The identified resources in Wyoming and Louisiana had lower geothermal gradients and were not evaluated. Reservoir modeling and techno-economic analysis were performed at Marys River Basin–North in Nevada. Geothermal energy production at this location is expected to have a levelized cost of energy (LCOE) ranging between 10 and 20 cents/kWh. Additional work may result in lower LCOE estimates at this location and at other attractive prospects in these three regions. The Great Basin carbonate and alluvial aquifer system of eastern Nevada and western Utah includes a lower carbonate aquifer unit, which has the potential for hosting both conduction- and convection-dominated geothermal systems. Mapping has identified lateral thickness variability expressed as a roughly 120–150-km-wide central corridor, which hosts the thickest and most continuous formations and extends from near Las Vegas north to the Idaho border. Purchase and analysis of privately held legacy seismic data could potentially compensate for the lack of sufficient data documenting measured depth to the lower carbonate aquifer unit. Multiple orogenies, extension episodes, and intrusive events have deformed and displaced the target formations of the lower carbonate aquifer unit. This structural complexity and potential for dipping reservoirs emphasizes the need for detailed geological, geophysical, and reservoir modeling. Heat flow within three Colorado sedimentary basins reviewed as part of this study was calculated in targeted studies by the Colorado Geologic Survey and Colorado School of Mines. These calculations are based on bottom-hole temperature data sets with significant limitations and some variability but produce values consistently higher than the global continental average of 65 mW/m2 for all three basins. Heat flow in the Raton Basin is the highest; however, permeability measurements from specific sedimentary formations with high heat flow have not been obtained. Promising formations for sedimentary geothermal systems were found in all three regions studied—Nevada-Utah, Colorado, and Texas. The next steps for developing sedimentary geothermal resources vary due to differences in available data and resource uncertainties. Additional scopes of work are recommended for identified basins in these three regions.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Bentonite buffer under high temperature: laboratory experiments and coupled process modeling

Abstract. Bentonite buffer in a geological repository will be simultaneously heated from decaying radioactive waste and hydrated from the surrounding host rock, triggering complex and coupled THMC (thermal–hydrological–mechanical–chemical) processes. Understanding the THMC behavior of bentonite-based engineered barrier system (EBS) is key to the evaluation and prediction of its long-term performance. Studies on the THMC process have been focused on conditions under 100 ∘C, as most design concepts impose a thermal limit of 100 ∘C in bentonite. Recently, studies under high-temperature conditions have been conducted to evaluate the possibility of raising the thermal limit and expanding the data/knowledge base to increase the confidence level. In this abstract, we present a series of bench-scale laboratory experiments at high temperatures (up to 200 ∘C) and the corresponding modeling work. Two sets of column tests were conducted, and each set consisted of two test columns: a control column undergoing only hydration (non-heated) and an experiment column experiencing both heating and hydration (heated). During the experiment, frequent X-ray computed tomography (CT) images were collected to provide a 3D visualization of the density distribution and present the spatiotemporal evolution of (1) hydration/dehydration, (2) clay swelling/shrinkage, (3) displacement, and (4) mineral precipitation. The two sets of tests differ with respect to several experimental conditions, such as bentonite type, compacted density and water content, water chemistry, and hydration pressure, but the important difference is that the first set used bentonite powder with a dry density of 1.28 g cm−3, whereas the second set used granulated bentonite (mixture of pellets and powder) with a dry density of 1.45–1.5 g cm−3. In both sets of experiments, a comprehensive post-dismantling characterization of bentonite samples was carried out after the column tests had been running for 1.5 years. Comparing non-heated and heated columns, the temperature gradient led to lower degree of homogenization of bentonite after bentonite became fully saturated; comparing the first and second sets, granulated and powdered bentonite exhibited drastically different hydration behavior. A THM model with a 2D axisymmetric grid system was used to interpret the data from the first set of tests. The model considers the combined impact of saturation, fluid pressure, and porosity change due to swelling/compression on the spatiotemporal distribution of bulk density and movement of the thermocouple modules. Observations from the tests help us understand the early perturbation of bentonite buffer under high temperature, and data from these tests improve the calibration of key constitutive hydrological and mechanical models and, therefore, enhance the modeling capability with respect to calculating the long-term evolution of bentonite buffer.

Zheng, Liange↗
DOE OSTI2024

Spatial variability in Arctic–boreal fire regimes influenced by environmental and human factors

Abstract Wildfire activity in Arctic and boreal regions is rapidly increasing, with severe consequences for climate and human health. Regional long-term variations in fire frequency and intensity characterize fire regimes. The spatial variability in Arctic–boreal fire regimes and their environmental and anthropogenic drivers, however, remain poorly understood. Here we present a fire tracking system to map the sub-daily evolution of all circumpolar Arctic–boreal fires between 2012 and 2023 using 375 m Visible Infrared Imaging Radiometer Suite active fire detections and the resulting dataset of the ignition time, location, size, duration, spread and intensity of individual fires. We use this dataset to classify the Arctic–boreal biomes into seven distinct ‘pyroregions’ with unique climatic and geographic environments. We find that these pyroregions exhibit varying responses to environmental drivers, with boreal North America, eastern Siberia and northern tundra regions showing the highest sensitivity to climate and lightning density. In addition, anthropogenic factors play an important role in influencing fire number and size, interacting with other factors. Understanding the spatial variability of fire regimes and its interconnected drivers in the Arctic–boreal domain is important for improving future predictions of fire activity and identifying areas at risk for extreme events.

Geology↗
DOE OSTI2022

Sea ice breakup and freeze-up indicators for users of the Arctic coastal environment

Abstract. The timing of sea ice retreat and advance in Arctic coastal waters varies substantially from year to year. Various activities, ranging from marine transport to the use of sea ice as a platform for industrial activity or winter travel, are affected by variations in the timing of breakup and freeze-up, resulting in a need for indicators to document the regional and temporal variations in coastal areas. The primary objective of this study is to use locally based metrics to construct indicators of breakup and freeze-up in the Arctic and subarctic coastal environment. The indicators developed here are based on daily sea ice concentrations derived from satellite passive-microwave measurements. The “day of year” indicators are designed to optimize value for users while building on past studies characterizing breakup and freeze-up dates in the open pack ice. Relative to indicators for broader adjacent seas, the coastal indicators generally show later breakup at sites known to have landfast ice. The coastal indicators also show earlier freeze-up at some sites in comparison with freeze-up for broader offshore regions, likely tied to earlier freezing of shallow-water regions and areas affected by freshwater input from nearby streams and rivers. A factor analysis performed to synthesize the local indicator variations shows that the local breakup and freeze-up indicators have greater spatial variability than corresponding metrics based on regional ice coverage. However, the trends towards earlier breakup and later freeze-up are unmistakable over the post-1979 period in the synthesized metrics of coastal breakup and freeze-up and the corresponding regional ice coverage. The findings imply that locally defined indicators can serve as key links between pan-Arctic or global indicators such as sea ice extent or volume and local uses of sea ice, with the potential to inform community-scale adaptation and response.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Predicting Short-Term Deformation in the Central Valley Using Machine Learning

Land subsidence caused by excessive groundwater pumping in Central Valley, California, is a major issue that has several negative impacts such as reduced aquifer storage and damaged infrastructures which, in turn, produce an economic loss due to the high reliance on crop production. This is why it is of utmost importance to routinely monitor and assess the surface deformation occurring. Two main goals that this paper attempts to accomplish are deformation characterization and deformation prediction. The first goal is realized through the use of Principal Component Analysis (PCA) applied to a series of Interferomtric Synthetic Aperture Radar (InSAR) images that produces eigenimages displaying the key characteristics of the subsidence. Water storage changes are also directly analyzed by the use of data from the Gravity Recovery and Climate Experiment (GRACE) twin satellites and the Global Land Data Assimilation System (GLDAS). The second goal is accomplished by building a Long Short-Term Memory (LSTM) model to predict short-term deformation after developing an InSAR time series using LiCSBAS, an open-source InSAR time series package. The model is applied to the city of Madera and produces better results than a baseline averaging model and a one dimensional convolutional neural network (CNN) based on a mean squared error metric showing the effectiveness of machine learning in deformation prediction as well as the potential for incorporation in hazard mitigation models. The model results can directly aid policy makers in determining the appropriate rate of groundwater withdrawal while maintaining the safety and well-being of the population as well as the aquifers’ integrity.

58 GEOSCIENCES↗
DOE OSTI2025

Stratigraphy‐Induced Localization of Microseismicity During CO 2 Injection in Illinois Basin

Abstract Subsurface fluid injection stimulates complex hydromechanical interaction, necessitating the integration of geomechanical data across spatial and temporal scales to consider the sophisticated behavior. Induced seismic response is usually associated with the complex reservoir architecture and pre‐existing features that are three‐dimensional, such as local stratigraphy, fractures, faults, and other discontinuities. This study encompasses laboratory characterization of the coupled hydromechanical response of cores extracted from rock formations in Illinois Basin: reservoir ‐ Mt. Simon sandstone, basal seal ‐ Argenta sandstone, and crystalline basement ‐ Precambrian rhyolite. High‐resolution numerical modeling allows considering the three‐dimensional complexity of the Illinois Basin Decatur Project with spatial resolution comparable to one of the active seismic surveys. A detailed reconstruction of the evolving state of stress in formations lacking direct stress measurements is achieved by numerical modeling that integrated laboratory‐derived hydromechanical properties, a porosity‐permeability relationship, active seismic data, and an inverted three‐dimensional porosity distribution. It appears that the microseismic clusters, mainly observed in the crystalline basement during the injection, are linked to zones experiencing more critically stressed conditions prior to injection. These zones have a potential for reactivation during the injection and are attributed to the specific local stratigraphy of the injection site, as well as transfer of triggering perturbations during the injection.

Bondarenko, N. [University of Illinois Urbana‐Cham↗