Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “geologic mapping”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Dating submarine landslides using the transient response of gas hydrate stability

Abstract Submarine landslides are prevalent on the modern-day seafloor, yet an elusive problem is constraining the timing of past slope failure. We present a novel age-dating technique based on perturbations to underlying gas hydrate stability caused by slide-impacted seafloor changes. Using three-dimensional (3-D) seismic data, we mapped an irregular bottom simulating reflection (BSR) underneath a submarine landslide in the Orca Basin, Gulf of Mexico. The irregular BSR mimics the pre-slide seafloor geometry rather than the modern bathymetry. Therefore, we suggest that the gas hydrate stability zone (GHSZ) is still adjusting to the post-slide sediment temperature. We applied transient conductive heat-flow modeling to constrain the response of the GHSZ to the slope failure, which yielded a most likely age of ca. 8 ka, demonstrating that gas hydrate can respond to landslides even on multimillennial time scales. We further provide a generalized analytical solution that can be used to remotely date submarine slides in the absence of traditional dating techniques.

Geology↗
DOE OSTI2024

Performance evaluation of the USGS velocity model for the San Francisco Bay Area

In this study, we evaluated the performance of the United States Geological Survey velocity model developed for the San Francisco Bay Area (SFBA), version 21.1. The evaluation was performed through high-resolution three-dimensional physics-based ground motion simulations of seven small-magnitude earthquakes (ranging from magnitude 3.8 to 4.4) that occurred on the eastern side of the San Francisco Bay. The simulations were performed in the frequency range from 0 to 5 Hz with a minimum shear-wave velocity of 250 m/s, which allowed the capture of wave propagation effects of the near-surface soft materials that characterize local basins. Based on the direct comparison of Fourier amplitude spectra between recorded and simulated ground motions for more than 250 stations, we found that the velocity model generally performs well in the frequency range of 0.2–5 Hz. The median value of the Fourier amplitude residuals was found to be near zero for all seven earthquakes. The slight over-prediction of 0.2 log-natural units at frequencies above 3 Hz in our simulations was attributed to the potentially inaccurate representation of the source radiation pattern by a double-couple point source model, and simple representation of shallow small-scale underground structural complexity in the velocity model. Maps of spectral amplitude differences between the simulated and recorded data were used to identify areas responsible for systematic ground motion over-predictions or under-predictions. For example, while some sub-domains over soft sediments show over-prediction patterns, the block east of the Hayward fault is prone to exhibit patterns of under-prediction. These maps can be used to guide future refinements of the SFBA velocity model. Since our simulation methodology allows for the decoupling of the source and wave propagation effects, the ground motion data generated by our simulations can also be used to quantify the epistemic uncertainty due to the velocity model, in empirically based ground motion estimates for the SFBA.

58 GEOSCIENCES↗
DOE OSTI2021

A multi-dimensional parametric study of variability in multi-phase flow dynamics during geologic CO 2 sequestration accelerated with machine learning

Successful geologic CO 2 storage projects depend on numerical simulations to predict reservoir performance during site selection, injection verification, and post-injection monitoring phases of the project. These numerical simulations solve non-linear sets of coupled partial differential equations, while accounting for multi-phase fluid dynamics on the basis of constitutive equations that are embedded into the solution scheme. As a consequence, individual simulations often require tens to hundreds of hours to complete on high-performance computing clusters. Moreover, laboratory experiments reveal that parametric functions for capillary pressure and relative permeability exhibit substantial variability, even within the same rock type. This combination of computational expense and wide-ranging parametric variability means that there remains substantial uncertainty in the behavior of multi-phase CO 2 -water systems, particularly in the context of feedbacks between relative permeability and capillary pressure. To bridge this knowledge gap, here we develop a novel workflow that utilizes physics-based numerical simulation to train an artificial neural network (ANN) emulator for interrogating the multivariate parameter space that governs both capillary pressure and relative permeability. With this approach, the ANN is trained to emulate both fluid pressure distribution and CO 2 saturation, which are then interrogated quantitatively to generate parametric response surface mappings with high-fidelity resolution. Results from this study initially show that capillary entry pressure is the dominant control on both CO 2 plume geometry and fluid pressure propagation when considering the combined effects of capillary pressure and relative permeability, particularly when phase interference is low and residual CO 2 saturation is high. Moreover, the ANN emulator provides tremendous computational speed-up by computing 2691 individual simulations in several minutes; whereas, the same simulation ensemble would have required ~3 years of simulation time using only physics-based simulation methods (25,000 times speed up).

58 GEOSCIENCES↗
DOE OSTI2026

A Play-Based Exploration of CO 2 Storage in the Illinois Basin (Final Technical Report)

This report documents the results of "A Play-Based Exploration of CO₂ Storage in the Illinois Basin" (DE-FE0032366), a project funded by the U.S. Department of Energy Office of Fossil Energy and Carbon Management and conducted by the Illinois State Geological Survey (ISGS) at the University of Illinois Urbana-Champaign in partnership with Visage Energy. The project adapted play-based exploration (PBE), a systematic basin-scale evaluation methodology from the petroleum industry, to screen areas of Illinois for commercial geologic carbon storage (GCS) in Cambro-Ordovician strata. The traditional play concept was expanded to encompass three play element groups, subsurface geologic factors, surface features, and societal factors, yielding 24 play elements with defensible suitability criteria applied through a five-tier classification scheme. An integrated geospatial database was assembled from ISGS, MGSC, MRCI, NATCARB, and public data sources, supported by significant data-improvement work including correction of legacy well locations, digitization of more than 1,300 well construction records using the DOE CATALOG team's OGRRE tool, compilation of a statewide 2D seismic database, and production of a refined fault and fold geodatabase.

58 GEOSCIENCES↗
DOE OSTI2022

Legacy surface change analyses of the 1993 Rock Valley earthquake sequence for direct comparison to planned NA-22 underground conventional high-explosive experiments (Source Physics Experiment 3 (SPE3) - RV/DC Task 1.3 FY22 Final Report)

Recent work under two previous phases of the NNSA NA-22 Source Physics Experiment (SPE) have shown that underground chemical high-explosive experiments can produce detectable surface changes that differ in spatial extent and vertical magnitude depending on the geologic media at the site (Schultz-Fellenz et al., 2018; 2020; Crawford et al., 2021). Neither of these two prior phases of SPE identified natural earthquake-related surface effects in the same region that occurred at a similar depth as the explosive experiments for direct comparison. While earthquakes can also produce surface changes that are detectable using remote sensing data analyses, it is expected that the pattern and spatial extent of surface changes would vary between earthquakes and explosions. However, no direct-observed surface-change signature comparison between earthquakes and explosions has ever been performed. The SPE Phase 3 Rock Valley/Direct Comparison (RV/DC) program presents a unique opportunity to investigate and characterize co-occurrence of both earthquakes and explosions. In this project, we worked to address five tasks in a workflow, as follows: 1. Identify and obtain existing high-resolution legacy satellite and aerial imagery as close in time before and after the 1993 Rock Valley earthquake sequence to temporally constrain the analyses. 2. Transform these pre-earthquake and post-earthquake datasets into digital elevation models (DEMs) using geospatial analysis software packages (e.g., ArcGIS, Agisoft Metashape, and Google Earth Engine). 3. Perform DEM differencing analyses to assess and quantify earthquake-related changes from the 1993 sequence, and develop map products that visualize these analyses. 4. Use the analyses from (3) to: (a) assess spatial distribution and magnitude of surface changes due to the 1993 Rock Valley earthquake sequence, and (b) determine parameters of forthcoming, planned explosion-related surface change data collection from sensors mounted on unmanned aerial vehicles (UAVs) (e.g., spatial extent of collection, design and density of survey control, sensors to deploy, forward speed and line spacing of UAV flight lines, flight altitude). 5. Develop a summary report on the analyses, including how the analyses define parameters and identify focus areas for any future surface change analytical field campaigns related to the explosive experiment. Analyzing these legacy data and identifying whether they can detect any surface changes related to the earthquake sequence facilitates opportunities for direct signature comparison of surface change from explosions at one location, which has never previously been performed. Comparing the surface change signatures from a co-located and depth-equivalent earthquake and an explosion could help to advance remote sensing event discrimination techniques. This report summarizes the work completed toward this ambitious goal.

42 ENGINEERING↗
DOE OSTI2026

Distributed Mafic Rock Resources for Carbon Mineralization in Arizona

Ex-situ carbon mineralization is a process by which CO2 is reacted with alkaline silicate minerals and rocks to produce stable carbonate materials, which can be used for other industrial processes. Arizona, U.S.A., hosts abundant surficial mafic rocks in three young volcanic fields, Geronimo-San Bernardino, San Francisco, and Springerville, and other distributed locations throughout the state. We created a Mafic Rock Resource Inventory (MRRI) that categorizes geochemical, physical, and textural characteristics of a diverse suite of surficial mafic rock samples and provide a benchmark reaction dataset parameterizing the temperature, pressure, and pH conditions best suited ex-situ mineralization in different rock types. MRRI data is publicly available online via a map-viewer. We establish two reaction condition sets, varied in temperature, pressure, and pH, where crystal-rich and glassy rocks reach maximum reaction extent and different carbonate phases are formed. Systematic ex-situ mineralization experiments on 21 diverse rock types show trends in geochemical, mineralogical, and reactivity behavior and establish maximum effective capture capacity. From this, scoria cones in three Arizona volcanic fields have a ~62 Gt effective CO₂ storage capacity with one of the fields having a ~42 Gt storage capacity in lava flows. Reactivity results have applications to alkaline mafic rock resources exposed globally, including producing additional effective storage capacity estimates and scaled commercialization of mafic rock ex-situ mineralization, should reaction extents be improved through advances in mineralization techniques. MRRI data were used to create a Direct Air Capture to Mineralization (DACM) systems model, technoeconomic analysis, and life-cycle assessment. These documents are presented as three appendices.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

HFTS-1 (Final Report)

Despite the long history of hydraulic fracturing, the optimal number of fracturing stages during multi-stage fracture stimulation in horizontal wells can be a complex subject. In addition to the increased expense of multistage fracturing in horizontal wells, increasing the number of fracturing stages does not always correlate with a commensurate increase in production. The problem is the application of a uniform fracture stimulation design to all stages with no consideration for geological variations along the wellbore. The result is an inefficient use and costly waste of energy and water. Optimization of the fracturing process requires an understanding of the cause-and-effect relationship between fracturing parameters and local geological properties at a given location along the wellbore. Realizing that the generalized rock mechanics theories and hypotheses are not truly applicable to fractured and laminated shales, quantifiable impacts of a shale’s geomechanical and depositional features are a prerequisite for design and implementation of optimized hydraulic fractures. The overarching goal of this project is to understand and define the relationships of shale geology and fracture dynamics using detailed field data that includes coring of the fracture domain. Finally, resource recovery from shale formations is estimated to be less than fifteen percent. Research proposed in this project will establish the foundations for investigating enhanced recovery techniques for increased resource recovery in existing fracture treated wells. Natural gas as an EOR fluid is in the initial stages of being used broadly in the oil and gas industry and as such, many aspects of the process need to be researched and addressed prior to widespread acceptance. In addition, re-fracturing techniques, particularly using cemented liners also need to be further investigated. This project explores several of these EOR techniques in greater detail, not just at the theoretical level, but through actual field test sites. This report delves into the research program, identified tasks, several studies, observations, results, and conclusions under the Hydraulic Fracturing Test Site (HFTS) project. This test project is subdivided into three distinct phases and each phase is discussed in detail under separate sections. The original test site for Phase 1 as well as the Phase 2 site are in Midland Basin in Texas and operated by Laredo Petroleum. In the first two phases, the program targeted the Wolfcamp formation in Midland Basin. Phase 3 of the study targeted the EagleFord play in South Texas with Devon Energy as the operator. Significant conclusions from this study include: 1) lateral and vertical fracture growth mapping as well as proppant transport behavior, 2) Use of field gas and its applicability for huff n puff EOR injection in the Midland Basin, 3) Development and validation of proppant log, both at core scale as well as in-fill well implementation by analyzing drilling mud return samples, 4) Applicability of liner refracturing including best practices for enhanced unit level productivity. This report will detail all these studies as well as technical results and conclusions from said studies. These include general subsurface characterization work, completion designs as implemented, studies looking at fracture geometry, high resolution microseismic study, advanced crosswell strain mapping studies including drainage characterization, core characterization results, proppant analysis (core and observation wells), geochemistry, etc. to name a few.

02 PETROLEUM↗
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 OSTI2021

Displacement Analysis of Geothermal Field Based on PSInSAR And SOM Clustering Algorithms A Case Study of Brady Field, Nevada—USA

The availability of free and high temporal resolution satellite data and advanced SAR techniques allows us to analyze ground displacement cost-effectively. Our aim was to properly define subsidence and uplift areas to delineate a geothermal field and perform time-series analysis to identify temporal trends. A Persistent Scatterer Interferometry (PSI) algorithm was used to estimate vertical displacement in the Brady geothermal field located in Nevada by analyzing 70 Sentinel-1A Synthetic-Aperture Radar (SAR) images, between January 2017 and December 2019. To classify zones affected by displacement, an unsupervised Self-Organizing Map (SOM) algorithm was applied to classify points based on their behavior in time, and those clusters were used to determine subsidence, uplift, and stable regions automatically. Finally, time-series analysis was applied to the clustered data to understand the inflection dates. The maximum subsidence is –19 mm/yr with an average value of –6 mm/yr within the geothermal field. The maximum uplift is 14 mm/yr with an average value of 4 mm/yr within the geothermal field. The uplift occurred on the NE of the field, where the injection wells are located. On the other hand, subsidence is concentrated on the SW of the field where the production wells are located. The coupling of the PSInSAR and the SOM algorithms was shown to be effective in analyzing the direction and pattern of the displacements observed in the field.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Evaluation of global teleconnections in CMIP6 climate projections using complex networks

In climatological research, the evaluation of climate models is one of the central research subjects. As an expression of large-scale dynamical processes, global teleconnections play a major role in interannual to decadal climate variability. Their realistic representation is an indispensable requirement for the simulation of climate change, both natural and anthropogenic. Therefore, the evaluation of global teleconnections is of utmost importance when assessing the physical plausibility of climate projections. We present an application of the graph-theoretical analysis tool δ-MAPS, which constructs complex networks on the basis of spatio-temporal gridded data sets, here sea surface temperature and geopotential height at 500 hPa. Complex networks complement more traditional methods in the analysis of climate variability, like the classification of circulation regimes or empirical orthogonal functions, assuming a new non-linear perspective. While doing so, a number of technical tools and metrics, borrowed from different fields of data science, are implemented into the δ-MAPS framework in order to overcome specific challenges posed by our target problem. Those are trend empirical orthogonal functions (EOFs), distance correlation and distance multicorrelation, and the structural similarity index. δ-MAPS is a two-stage algorithm. In the first place, it assembles grid cells with highly coherent temporal evolution into so-called domains. In a second step, the teleconnections between the domains are inferred by means of the non-linear distance correlation. We construct 2 unipartite and 1 bipartite network for 22 historical CMIP6 climate projections and 2 century-long coupled reanalyses (CERA-20C and 20CRv3). Potential non-stationarity is taken into account by the use of moving time windows. The networks derived from projection data are compared to those from reanalyses. Our results indicate that no single climate projection outperforms all others in every aspect of the evaluation. But there are indeed models which tend to perform better/worse in many aspects. Differences in model performance are generally low within the geopotential height unipartite networks but higher in sea surface temperature and most pronounced in the bipartite network representing the interaction between ocean and atmosphere.

58 GEOSCIENCES↗
DOE OSTI2021

Tree Canopies Reflect Mycorrhizal Composition

Mycorrhizae alter global patterns of CO 2 fertilization, carbon storage, and elemental cycling, yet knowledge of their global distributions is currently limited by the availability of forest inventory data. Here, we show that maps of tree-mycorrhizal associations (hereafter “mycorrhizal maps”) can be improved by the novel technology of imaging spectroscopy because mycorrhizal signatures propagate up from plant roots to impact forest canopy chemistry. We analyzed measurements from 143 airborne imaging spectroscopy surveys over 112,975 individual trees collected across 13 years. Results show remarkable accuracy in capturing ground truth observations of mycorrhizal associations from canopy signals across disparate landscapes (R 2 = 0.92, p < 0.01). Additionally, upcoming imaging spectroscopy satellite missions can reveal new insights into landscape-scale variations in water, nitrogen, phosphorus, carotenoid/anthocyanin, and cellulose/lignin composition. Applied globally, this approach could improve the spatial precision of mycorrhizal distributions by a factor of roughly 104 and facilitate the incorporation of dynamic shifts in forest composition into Earth system models.

58 GEOSCIENCES↗
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 OSTI2025

GIScience in the era of Artificial Intelligence: a research agenda towards Autonomous GIS

The advent of generative AI exemplified by large language models (LLMs) opens new ways to represent and compute geographic information and transcends the process of geographic knowledge production, driving geographic information systems (GIS) towards autonomous GIS. Leveraging LLMs as the decision core, autonomous GIS can independently generate and execute geoprocessing workflows to perform spatial analysis. In this vision paper, we further elaborate on the concept of autonomous GIS and present a conceptual framework that defines its five autonomous goals, five levels of autonomy, five core functions, and three operational scales. We demonstrate how autonomous GIS could perform geospatial data retrieval, spatial analysis, and map making with four proof-of-concept GIS agents. We conclude by identifying critical challenges and future research directions, including fine-tuning and self-growing decision-cores, autonomous modelling, and examining the societal and practical implications of autonomous GIS. By establishing the groundwork for a paradigm shift in GIScience, this paper envisions a future where GIS moves beyond traditional workflows to autonomously reason, derive, innovate, and advance geospatial solutions to pressing global challenges. Meanwhile, we emphasize that as we design and deploy increasingly intelligent geospatial systems, we carry a responsibility to ensure they are developed in socially responsible ways, serve the public good, and support the continued value of human geographic insight in an AI-augmented future.

Autonomous GI↗
DOE OSTI2020

SIMS Analysis of d18O and d13C of Carbonate Cements in Shale & Mudstone. Final report

This project has had two primary goals: 1) improving new technology for in situ microanalysis of stable isotope ratios by SIMS (secondary ion mass spectrometry) and 2) evaluating the importance of in situ analysis for studies of carbonate cements within shales of the Bakken Petroleum System. The systematic relations of stable isotope ratios provide an important tool to help resolve critical questions in geology. Stable isotope analysis of carbonate cements in shale constrain many topics of interest to DOE including: CO2 sequestration, climate change, and the genesis and quality of rock reservoirs for petroleum, natural gas, and water. However, this evidence is often preserved as fine-scale zonation, which is homogenized by conventional bulk analysis of powders. Recently, analysis of ultra-small spots in polished samples at the University of Wisconsin WiscSIMS Lab. has documented stable isotope zoning in carbonates and other diagenetic cements at smaller scales than previously possible. We have made detailed studies of oxygen and carbon isotope ratios in diagenetic cements from shales of the Bakken Formation from the Williston basin and elsewhere. The Bakken contains both reservoir and source rocks for petroleum. Hydrofracturing has greatly enhanced production of oil in these tight reservoirs. The sediments contain zoned carbonate cements with growth textures that record multiple events. The growth of carbonate (dolomite) cements has been proposed to control reservoir quality, but previous stable isotope studies of homogenized powders could not resolve individual generations of cement and this proposal remains untested. Our better characterization of the multiple growth zones in these carbonates will aid in understanding the genesis of these tight reservoirs and the production of unconventional oil and gas. This study has also developed SIMS standards for analysis of all compositions of Ca-Mg-Fe carbonate and selected Mn-carbonates and hydromagnesite. We have refined techniques for improved accuracy of SIMS stable isotope analysis, and for sub-1-micron spot sizes. We have established protocols for mapping and correlating sample analyses at micron-scale using QGIS software.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2023

Fracture-matrix fluid exchange in oil-bearing unconventional mudstones

The poromechanical properties of unconventional reservoir materials are in large part dictated by their mineralogy. Since these properties govern the response to stress experienced during hydraulic fracturing, fluid production, and fluid injection, they play a central role in the formation of microcracks or bedding delaminations which ultimately dominate mass transport. In this work we study access to the porosity of end member unconventional reservoir materials, where the end members are predominantly dictated by carbonate content. Access to the porosity is quantified using state of the art 3D x-ray computed tomography coupled with physics informed data analytics. Xenon gas, which attenuates x-rays, provides a spatiotemporal map of access to the porosity. The accessible porosity is quantified over a range of net confining stress relevant to the manmade disturbances listed above. These experiments demonstrate that heavily carbonated mudstones are nearly impermeable at the core (~ cm) scale, while carbonate free analogues afford better access to the microstructure. Consistent with previous qualitative 2D radiographs, access to the interior of the clastic mudstones is first observed along planar microcracks, followed by slow penetration into the surrounding matrix. Physics informed data analytics of the 3D tomography measurements presented here show that these microcracks do not permit uniform access to the adjacent rock matrix. In addition, variation of the effective pressure elucidates the mechanisms that govern fracture/matrix fluid exchange. Under conditions consistent with hydrocarbon production fluid accumulates in the immediate vicinity of the nearest microcrack. While there is clear evidence that, as intended, part of this accumulation is from the more distant matrix, fluid is also squeezed out of the microcrack. The fluid build-up at the microcrack indicates that migration out of the rock is hindered by the coupled poroelastic response of the microcrack and adjacent rock matrix. We show that these mechanisms ultimately account for the meager oil recovery factors realized in practice. These insights have implications for making reservoir scale predictions based on core scale observations, and provide a basis for devising new asset development techniques to access more porosity, and enhance fluid extraction. Finally, these findings shed light on key features and mechanisms that govern shale storage capacity, with relevance to other important industrial processes, such as geologic CO 2 storage.

58 GEOSCIENCES↗
DOE OSTI2023

Estimating Individual Tree Mortality in the Sierra Nevada Using Lidar and Multispectral Reflectance Data

Widespread tree mortality events occur during periods of severe drought in temperate conifer forests and are expected to become more frequent in many areas due to climate change. Improved mapping of individual tree mortality is needed to identify risk factors and design effective conservation strategies. In this study, we used National Ecological Observatory Network (NEON) lidar and multispectral reflectance airborne observations to map individual tree mortality over a 160 km 2 area during and after the 2012–2016 drought for two sites in California's Sierra National Forest. We used NEON lidar to derive tree locations and crown perimeters and multispectral data to map tree mortality for more than 1 million trees. We found that 25.4% of the trees in our study area died between 2013 and 2017, with considerably higher mortality at the lower-elevation Soaproot Saddle site. Between 2017 and 2019, an additional 2.0%–2.8% of the trees died each year. Two wildfires in 2020 and 2021 increased tree mortality within burned area perimeters by 49%–89% between 2019 and 2021. Consistent with previous work, we found that tree mortality risk increased as a function of tree height. Tree mortality was positively associated with distance from rivers, trees per hectare, and decreasing slope at the lower elevation site. In contrast, increasing slope was positively associated with tree mortality at the higher elevation site. Our approach and dataset provide a means to study the combined effects of drought and wildfire on tree mortality and may improve projections of forest resilience under a changing climate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

University of Hawai‘i, Shallow Geothermal Resources: Energy Technology Innovation Partnership Project (Final Report)

Scientists at Lawrence Berkeley National Laboratory (Berkeley Lab) have teamed up with the University of Hawai‘i at Manoa (UH Manoa) through the U.S. Department of Energy’s Energy Technology Innovation Partnership Project to evaluate the technological and market feasibility of shallow geothermal heat exchanger (GHE) technology. UH requested this analysis to evaluate opportunities in building cooling, energy efficiency, and emissions reduction applications in Hawai‘i. UH has an abundance of geologic and geothermal data and is looking to the national labs’ expertise to execute this analysis. UH is also interested in investigating policy, regulatory, and business conditions advantageous for implementation of a pilot project and more broad deployment of this technology in Hawai‘i. In many locations around the world, the demands for heating and cooling are roughly balanced over the course of the year, so GHEs do not cause significant long-term changes in subsurface temperature. This is not the case in Hawai’i, where the demand for heating is very small, meaning that, over time, GHEs will add heat to the subsurface. If temperatures increase significantly, GHE systems will not work as designed. Regional groundwater flow has the potential to sweep heated water away from boreholes, thereby maintaining the functionality of the GHE system. Significant regional groundwater flow requires two things: a sufficiently large driving hydraulic head gradient (usually closely related to surface topography), and sufficient porosity and permeability to enable groundwater to flow in large enough quantities to enable near-borehole temperatures to be maintained at ambient values. Hawai‘i’s volcanic terrain offers ample surface topographic variation. The lava itself shows an extremely large range of porosity and permeability, so sites with large enough values of these properties must be selected. Numerical modeling of coupled groundwater and heat flow can be used to determine how large is large enough. Primarily, closed-loop systems have been investigated. Other options considered are open-loop systems and using cool seawater as the chilling source. Project work investigated the feasibility of GHE technology at two scales. At the island scale, GIS layers of various attributes relevant for GHE were combined to develop an overall favorability map for employing GHE in Hawai‘i. At the local scale, a hydrogeologic model for the subsurface component of a closed-loop system was developed for the Stan Sheriff Center at the UH Manoa campus. This site is considered promising because the rock below and immediately downgradient of the borefield is highly permeable, consisting of a subsurface karst system (limestone containing high-permeability open channels), which is underlain by a thick, high-permeability fractured basalt. Moreover, the site is near the base of the Ko‘olau Range, providing a large hydraulic head gradient. Thus, groundwater flow through the site is expected to be large, enabling efficient removal of heated groundwater. A full-GHE-system model of the site was also developed, with a simplified representation of the subsurface, in which groundwater flow is not considered and heat transfer is purely by conduction. Using the building cooling load data provided by UH, simulation results show that with groundwater flow present, a GHE can operate successfully for at least 10 years, but with no groundwater flow, the subsurface begins to heat up after only one year of operation, making the GHE unviable within 2-6 years. The team also developed a techno-economic model for this site to compare the cost of cooling using a GHE system with the costs of operating the current air-conditioning system. The GHE system is advantageous economically if favorable tax incentives and interest rates can be obtained.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

Transforming Uinta Basin Earth Materials for Advanced Products (TUBE-MAP)

The objectives of this project were to quantify, assess, and plan to enable the transformation of Uinta Basin earth resources, such as coal, oil shale, resin, rare earth elements, and critical minerals into high value metal, mineral, and carbon-based products. The specific major goals were 1) basinal assessments and initial planning (Task 2), 2) basinal assessment for waste stream reuse with associated plan development (Task 3), 3) basinal strategies development for infrastructure, industries, and business (Task 4), 4) technology assessment, development, and field-testing plan (Task 5), 5) technology innovation center plan (Task 6), and 6) stakeholder outreach and education plan (Task 7).

01 COAL, LIGNITE, AND PEAT↗