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 451 records · Page 25

Compressional residual stress in Bastogne boudins revealed by synchrotron X-ray microdiffraction

Lattice distortions in crystals can be mapped at the micron scale using synchrotron X-ray Laue microdiffraction (μXRD). From lattice distortions the shape and orientation of the elastic strain tensor can be derived and interpreted in terms of residual stress. We apply the new method to vein quartz from the original boudinage locality at Bastogne, Belgium. Furthermore, a long-standing debate surrounds the kinematics of the Bastogne boudins. The μXRD measurements reveal a shortening residual elastic strain, perpendicular to the vein wall, corroborating the model that the Bastogne boudins formed by layer-parallel shortening and not by layer-parallel extension, as is in the geological community generally inferred by the process of boudinage.

58 GEOSCIENCES↗

Sampling Density in Regional Exploration and Environmental Geochemical Studies: A Review

A variety of factors contribute to the proper and efficient evaluation of geochemical programs including sampling medium, sample size and sample fraction, sampling density, digestion procedure, and analytical techniques. Among the aforementioned factors, sampling density is of great importance, affecting project ability to detect anomalous sites, which are either economically viable or environmentally hazardous. It is quite critical to have a sufficient number of samples to delineate geochemical anomalies or baselines and to reveal the relative variations in elemental concentrations in a study area. The resolution of the subsequent maps obtained by sampling greatly depends on sampling density. While a large number of samples may generate redundant data, a small number of samples smooth out the patterns and may completely miss anomalies. This review provides a platform for applied geochemists who desire to design a sampling scheme for a regional geochemical program and do not have the financial resources to perform an orientation survey, which will give them the necessary information for designing, according to project objectives, a cost-effective regional geochemical survey for their study area. It is advisable to integrate during the desktop study all available information, such as geological, structural, geophysical, remote sensing data and topographical data by using different data fusion algorithms to arrive at a sampling density fit-for-purpose of the target area. This in turn facilitates an effective spatial positioning of sample locations. Quantitative analyses of the efficiency of different sampling densities for delineating geochemical patterns have not been well documented, warranting more attention and research.

58 GEOSCIENCES↗

A configurational force for adaptive re-meshing of gradient-enhanced poromechanics problems with history-dependent variables

Here, we introduce a mesh-adaption framework that employs a multi-physical configurational force and Lie algebra to capture multiphysical responses of fluid-infiltrating geological materials while maintaining the efficiency of the computational models. To resolve sharp gradients of both displacement and pore pressure, we introduce an energy-estimate-free re-meshing criterion by extending the configurational force theory to consider the energy dissipation due to the fluid diffusion and the gradient-dependent plastic flow. To establish new equilibria after remeshing, the local tensorial history-dependent variables at the integration points are first decomposed into spectral forms. Then, the principal values and directions are projected onto smooth fields interpolated by the basis function of the finite element space via the Lie-algebra mapping. Our numerical results indicate that this Lie algebra operator in general leads to a new trial state closer to the equilibrium than the ones obtained from the tensor component mapping approach. A new configurational force for dissipative fluid-infiltrating porous materials that exhibit gradient-dependent plastic flow is introduced such that the remeshing may accommodate the need to resolve the sharp pressure gradient as well as the strain localization. The predicted responses are found to be not influenced by the mesh size due to the micromorphic regularization, while the adaptive meshing enables us to capture the width of deformation bands without the necessity of employing fine mesh everywhere in the domain.

42 ENGINEERING↗

The formation mechanisms of sedimentary pyrite nodules determined by trace element and sulfur isotope microanalysis

Redox sensitive trace elements in pyrite, including nodules, are increasingly used to infer the chemical conditions of ancient oceans—but considerable uncertainty remains regarding the mechanism and timing of nodule formation. Resolving these uncertainties is important because pyrite nodules must form in connection with the overlying water column, rather than during late diagenesis, to reflect the composition of the global ocean. Existing models for pyrite nodule formation have been specific to pyrite textures from individual sites, and we lack a unified model that can explain the compositional and textural diversity observed in nodules from different localities. In this study we examine ten pyrite nodules from several geological periods (Neoarchean to Carboniferous) using in situ LA-ICP-MS and SHRIMP-SI analyses. We present transects of spot analyses of trace elements (As, Ag, Cu, Co, Ni, Sb and Se) and S isotope ratios for each nodule. The pyrite nodules can be classified according to three main categories: those with (1) little to no trace element or isotopic zonation of the nodule from core to margin, (2) strong zonation from core to margin, and (3) minor zonation near the core but more significant zonation near the margin of the nodule. We further illustrate this zonation with a NanoSIMS element map from an additional pyrite nodule. These results are interpreted to indicate nodule formation along a spectrum between two end-member mechanisms. We suggest that the absence of trace element or isotopic zonation reflects nodule growth by a pathway that is analogous to the pervasive growth mechanism for carbonate nodules. This model involves the production of many nucleation sites that are evenly distributed within the volume that the nodule eventually occupies. Consequently, this mechanism results in a chemically homogenous nodule. Pyrites formed this way are suitable for paleoceanographic reconstruction. The other end-member mechanism is analogous to the concentric growth model for carbonate concretions. In this scenario, the core of the nodule forms first and is followed by the addition of concentric layers—each with a progressively different trace element content and δ 34 S signature as diagenesis progresses. Despite having limited utility for econstructing ancient seawater, these late forming nodules may track the evolving availability of bioessential trace elements for the subsurface biosphere with important implications for global biogeochemical cycles. Finally, spatial trends in trace elements and S isotopes thus speak to the mechanisms of pyrite nodule formation and provide a framework for evaluating nodule suitability for a range of paleoenvironmental studies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dynamic phases of active matter systems with quenched disorder

Depinning and nonequilibrium transitions within sliding states in systems driven over quenched disorder arise across a wide spectrum of size scales ranging from atomic friction at the nanoscale, flux motion in type II superconductors at the mesoscale, colloidal motion in disordered media at the microscale, and plate tectonics at geological length scales. Here we show that active matter or self-propelled particles interacting with quenched disorder under an external drive represents a class of system that can also exhibit pinning-depinning phenomena, plastic flow phases, and nonequilibrium sliding transitions that are correlated with distinct morphologies and velocity-force curve signatures. When interactions with the substrate are strong, a homogeneous pinned liquid phase forms that depins plastically into a uniform disordered phase and then dynamically transitions first into a moving stripe coexisting with a pinned liquid and then into a moving phase-separated state at higher drives. We numerically map the resulting dynamical phase diagrams as a function of external drive, substrate interaction strength, and self-propulsion correlation length. These phases can be observed for active matter moving through random disorder. Lastly, our results indicate that intrinsically nonequilibrium systems can exhibit additional nonequilibrium transitions when subjected to an external drive.

74 ATOMIC AND MOLECULAR PHYSICS↗

Plutonium Migration during the Leaching of Cemented Radioactive Waste Sludges

One of the most challenging components of the UK nuclear legacy is Magnox sludge, arising from the corrosion of Mg alloy-clad irradiated metallic U fuel that has been stored in high pH ponds. The sludges mainly comprise Mg hydroxide and carbonate phases, contaminated with fission products and actinides, including Pu. Cementation and deep geological disposal is one option for the long-term management of this material, but there is a need to understand how Pu may be leached from the waste, if it is exposed to groundwater. Here, we show that cemented Mg(OH) 2 powder prepared with Pu(IV) aq is altered on contact with water to produce a visibly altered ‘leached zone’, which penetrates several hundred microns into the sample. In turn, this zone shows slow leaching of Pu, with long-term leaching rates between 1.8–4.4 × 10 -5 % of total Pu per day. Synchrotron micro-focus X-ray fluorescence mapping identified decreased Pu concentration within the ‘leached zone’. A comparison of micro-focus X-ray absorption spectroscopy (µ-XAS) spectra collected across both leached and unleached samples showed little variation, and indicated that Pu was present in a similar oxidation state and coordination environment. Fitting of the XANES spectra between single oxidation state standards and EXAFS modeling showed that Pu was present as a mixture of Pu(IV) and Pu(V). The change in Pu oxidation from the stock solution suggests that partial Pu oxidation occurred during sample ageing. Similarity in the XAS spectra from all samples, with different local chemistries, indicated that the Pu oxidation state was not perturbed by macro-scale variations in cement chemistry, surface oxidation, sample aging, or the leaching treatment. These experiments have demonstrated the potential for leaching of Pu from cementitious waste forms, and its underlying significance requires further investigation.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Evaluation of Geologic Carbon Storage Resource Estimates (SREs) of Cambrian Ordovician Units within the CarbonSAFE Prefeasibility Study Region Subtasks 4.2 and 4.3 – Topical Report

This report by the Indiana Geological & Water Survey (IGWS) assessed the storage resource estimate (SRE) of Cambrian-Ordovician stratigraphic units located within the CarbonSAFE prefeasibility area within the Illinois Basin. A comprehensive data set of wireline logs and petrophysical information was used. These data include core analysis for porosity and permeability, and existing well data including location and stratigraphic information. This report includes storage resource estimates (SREs) for four potential storage reservoirs (Upper Ordovician Trenton Limestone/Black River Group, Middle Ordovician St. Peter Sandstone, Lower Ordovician and Upper Cambrian Knox Supergroup, including a separate assessment for the lower Knox Potosi Dolomite. The resultant SREs volumes are displayed using georeferenced maps and boxplots. The latter allow for comparing data statistics (mean values and variability) between five methods employed for the SRE calculations. Differences observed in SRE results from the five methods are mainly attributable to differences in the data used for porosity in each method. ArcMap 10.5.1 software was used to portray SREs to help visualize spatial variance of estimates for each methodology, and more importantly, to highlight those areas having the greatest total storage potential estimates.

01 COAL, LIGNITE, AND PEAT↗

Geological Carbon Sequestration Storage Resource Estimates for the Ordovician St. Peter Sandstone, Illinois and Michigan Basins, USA

The Cambro-Ordovician strata of the Midwest of the United States is a primary target for potential geological storage of CO2 in deep saline formations. The objective of this project is to develop a comprehensive evaluation of the Cambro-Ordovician strata in the Illinois and Michigan Basins above the basal Mount Simon Sandstone since the Mount Simon is the subject of other investigations including a demonstration-scale injection at the Illinois Basin Decatur Project. The primary reservoir targets investigated in this study are the middle Ordovician St Peter Sandstone and the late Cambrian to early Ordovician Knox Group carbonates. The topic of this report is a regional-scale evaluation of the geologic storage resource potential of the St Peter Sandstone in both the Illinois and Michigan Basins. Multiple deterministic-based approaches were used in conjunction with the probabilistic-based storage efficiency factors published in the DOE methodology to estimate the carbon storage resource of the formation. Extensive data sets of core analyses and wireline logs were compiled to develop the necessary inputs for volumetric calculations. Results demonstrate how the range in uncertainty of storage resource estimates varies as a function of data availability and quality, and the underlying assumptions used in the different approaches. In the simplest approach, storage resource estimates were calculated from mapping the gross thickness of the formation and applying a single estimate of the effective mean porosity of the formation. Results from this approach led to storage resource estimates ranging from 3.3 to 35.1 Gt in the Michigan Basin, and 1.0 to 11.0 Gt in the Illinois Basin at the P10 and P90 probability level, respectively. The second approach involved consideration of the diagenetic history of the formation throughout the two basins and used depth-dependent functions of porosity to derive a more realistic spatially variable model of porosity rather than applying a single estimate of porosity throughout the entire potential reservoir domains. The second approach resulted in storage resource estimates of 3.0 to 31.6 Gt in the Michigan Basin, and 0.6 to 6.1 Gt in the Illinois Basin. The third approach attempted to account for the local-scale variability in reservoir quality as a function of both porosity and permeability by using core and log analyses to calculate explicitly the net effective porosity at multiple well locations, and interpolate those results throughout the two basins. This approach resulted in storage resource estimates of 10.7 to 34.7 Gt in the Michigan Basin, and 11.2 to 36.4 Gt in the Illinois Basin. A final approach used advanced reservoir characterization as the most sophisticated means to estimating storage resource by defining reservoir properties for multiple facies within the St Peter formation. This approach was limited to the Michigan Basin since the Illinois Basin data set did not have the requisite level of data quality and sampling density to support such an analysis. Results from this approach led to storage resource estimates of 15.4 Gt to 50.1 Gt for the Michigan Basin. The observed variability in results from the four different approaches is evaluated in the context of data and methodological constraints, leading to the conclusion that the storage resource estimates from the first two approaches may be conservative, whereas the net porosity based approaches may over-estimate the resource.

Barnes, David↗

Integrating Magnetotellurics, Soil Gas Geochemistry and Structural Analysis to Identify Hidden, High Enthalpy, Extensional Geothermal Systems

Our research group comprising personnel from University of Utah/EGI (UU/EGI), University of Nevada at Reno (UNR) and Lawrence Berkeley National Laboratory (LBNL) has been applying magnetotellurics (MT), diagnostic structural affiliations, and soil gas flux and geochemistry to assist in identifying hidden, high-enthalpy geothermal systems in extensional regimes of the U.S. Great Basin. We are specifically looking for upwelling low-resistivity and geologic structures that can carry fluids from magmatic or high-grade metamorphic conditions in the deep crust upward to exploitable depths, and to verify the nature of the deep sources through soil gas geochemical compositions. The project was motivated by prior MT transect coverage of western and central Nevada centered upon the Dixie Valley geothermal system where such favorable indicators were first recognized. The high-angle structures are taken to be fluidized fault zones connecting deep magmatic/metamorphic activity with the geothermal system but the concept requires verification by testing at other systems. The project was set up with a two-phased organization. Phase I was carried out at the McGinness Hills system, central Nevada, where Ormat’s flagship power facility is located and a considerable amount of pre-existing data was available. Transect MT data also showed a strong low-resistivity upwelling originating from interpreted deep crustal magmatic underplating. Controlling structures on production as indicated by Ormat data and new mapping by UNR were favorable to dilatency being an accommodation zone between major normal faults of opposing dip. A 3D MT survey and inversion confirmed the existence of the steep low-resistivity zone dipping ESE toward the deep crust and placed N-S bounds upon the feature. Ormat personnel cooperated with PI Wannamaker in the field in providing fluid samples from the production intervals for He isotope composition. Elevated 3He was verified through analysis at LBNL confirming that there was a magmatic connection into the producing system. High CO2 soil gas flux including possibly metamorphic 13C and 14C component was measured over the area of dilatent structures. Hence, the triad of indicators posed above was confirmed in Phase I. Subsequently, Phase II of the project proceeded in the greenfield Kumiva-Blackrock Desert district of northwestern Nevada to see if a new system could be identified. Transect MT data also showed a strong low-resistivity upwelling originating from interpreted deep crustal magmatic underplating. However, given that there was no exploitable geothermal system recognized at the outset in the region, the collection of new geoscientific data was required. An areal MT survey of 131 sites over the Kumiva-Blackrock Desert district was completed, and a 3D inversion produced using the in-house, DOE-supported finite element inversion algorithm. Particular MT low resistivity upwellings that received followup study occurred under the west flank of the Seven Troughs Range, under the playa immediately west of the Blue Wing Mountains, and under northern Granite Springs Valley. Structural assessment of the project area by Co-I J. Faulds at UNR provided numerous favorable Quaternary fault settings, which were correlated to the MT upwelling structures. Extensive profiling for CO2 soil gas flux over favorable geophysical and geological structure was performed by LBNL personnel. Soil gas flux anomalies generally were not large but did show correlation with MT upwelling structure and favorably dilatent geological structures. Isotope analyses showed presence of possible inorganic/metamorphic 13C but 14C concentrations did not exceed background values. We view the initial concept of a confluence of MT low-resistivity upwelling, favorably dilatent 3D geological structure, and elevated soil gas flux including 13C component to be supported by the further evidence of this project although the indicators in the Phase II study were more diffuse. The northern Granite Springs Valley structure is receiving followup in the UNR Play Fairway Analysis, being in the southern termination of a major normal fault with silicified soil and warm temperature probe results

15 GEOTHERMAL ENERGY↗

Jurassic deserts and CO2 sequestration: the Glen Canyon Group in Utah, western USA

The Early Jurassic in the Western Interior of USA consisted of an extensive desert environment, including one of the largest ergs in geologic history. The Navajo erg has been estimated to extend as much as 2.2 million square kilometers, though the preserved extent is somewhat smaller. The Mesozoic was a global greenhouse phase, and during the Jurassic this region experienced fluctuations in climate aridity, reflected in the depositional environment. In addition to extensive dunefields, interdune deposits (lakes and oases) and fluvial systems are also documented within formations of the Glen Canyon Group. The Navajo Sandstone has received much attention as a potential CO2 injection target in recent years. It consists of thick, aeolian sandstones with high porosity and permeability, occurs in both outcrop and subcrop, and has industry data. Mapping of stratigraphic and sedimentologic changes within the Navajo Sandstone has been undertaken in localized areas, however piecing together these studies and developing regional models for CO2 potential has not received as much attention. There is significant industry data across Utah, and utilizing this data to pivot from a focus on hydrocarbon extraction to CO2 injection is an effective way to move forward with green energy, and to meet carbon neutral emission goals. Using legacy well data, we are developing a more comprehensive study of the Navajo Sandstone as a potential CO2 reservoir, in addition to identifying other zones of interest within the Glen Canyon Group, and improving understanding of stratigraphic complexity within one of the most significant aeolian systems in the world.

Mahon, Elizabeth↗

The Colorado East River Community Observatory Data Collection

Abstract The U.S. Department of Energy's (DOE) Colorado East River Community Observatory (ER) in the Upper Colorado River Basin was established in 2015 as a representative mountainous, snow‐dominated watershed to study hydrobiogeochemical responses to hydrological perturbations in headwater systems. The ER is characterized by steep elevation, geologic, hydrologic and vegetation gradients along floodplain, montane, subalpine, and alpine life zones, which makes it an ideal location for researchers to understand how different mountain subsystems contribute to overall watershed behaviour. The ER has both long‐term and spatially‐extensive observations and experimental campaigns carried out by the Watershed Function Scientific Focus Area (SFA), led by Lawrence Berkeley National Laboratory, and researchers from over 30 organizations who conduct cross‐disciplinary process‐based investigations and modelling of watershed behaviour. The heterogeneous data generated at the ER include hydrological, genomic, biogeochemical, climate, vegetation, geological, and remote sensing data, which combined with model inputs and outputs comprise a collection of datasets and value‐added products within a mountainous watershed that span multiple spatiotemporal scales, compartments, and life zones. Within 5 years of collection, these datasets have revealed insights into numerous aspects of watershed function such as factors influencing snow accumulation and melt timing, water balance partitioning, and impacts of floodplain biogeochemistry and hillslope ecohydrology on riverine geochemical exports. Data generated by the SFA are managed and curated through its Data Management Framework. The SFA has an open data policy, and over 70 ER datasets are publicly available through relevant data repositories. A public interactive map of data collection sites run by the SFA is available to inform the broader community about SFA field activities. Here, we describe the ER and the SFA measurement network, present the public data collection generated by the SFA and partner institutions, and highlight the value of collecting multidisciplinary multiscale measurements in representative catchment observatories.

54 ENVIRONMENTAL SCIENCES↗

Model Data Archive for Manuscript Titled "Evaluation of a Coupled Surface–Subsurface Hydrologic Model Using Dense Water‑Level Sensors in a Mixed Urban–Rural Watershed"

This archive provides scripts, input files, and datasets used for the implementation and evaluation of a fully coupled surface–subsurface hydrologic model in the Neches River Basin, southeast Texas. The study uses the Advanced Terrestrial Simulator (ATS) to simulate coupled surface–subsurface hydrologic processes over a mixed urban–rural watershed and evaluates model performance using a dense network of 136 in situ water-level sensors, nine U.S. Geological Survey (USGS) stream gauges, and SSEBop-derived evapotranspiration estimates during the period October 2014–June 2024. The workflow is implemented primarily in Python 3 using the Watershed Workflow package. The Jupyter notebooks can be executed using open-source software such as Anaconda JupyterLab or Visual Studio Code. Other data files include TXT, CSV, XML, SHP, TIF, NetCDF, HDF5, and ExodusII files, which can be processed using the provided Python scripts. ATS input files are provided in XML format and can be edited using any commonly used text editor. This archive contains: *Scripts and input files used to generate the ATS model setup, including watershed discretization, mesh generation, parameter mapping, and model configuration. *Jupyter notebooks used for preprocessing observational data, evaluating streamflow, water levels, and evapotranspiration, computing performance metrics, and generating the figures presented in the manuscript. *ATS simulation outputs and processed observational datasets, including OneRain and DD6 water-level sensors, USGS streamflow observations, GIS data, and supporting spatial datasets used throughout the study.

Dense water-level sensor network↗

DEEPEN 3D PFA Favorability Models and 2D Favorability Maps at Newberry Volcano

DEEPEN stands for DE-risking Exploration of geothermal Plays in magmatic ENvironments. Part of the DEEPEN project involved developing and testing a methodology for a 3D play fairway analysis (PFA) for multiple play types (conventional hydrothermal, superhot EGS, and supercritical). This was tested using new and existing geoscientific exploration datasets at Newberry Volcano. This GDR submission includes images, data, and models related to the 3D favorability and uncertainty models and the 2D favorability and uncertainty maps. The DEEPEN PFA Methodology is based on the method proposed by Poux et al. (2020), which uses the Leapfrog Geothermal software with the Edge extension to conduct PFA in 3D. This method uses all available data to build a 3D geodata model which can be broken down into smaller blocks and analyzed with advanced geostatistical methods. Each data set is imported into a 3D model in Leapfrog and divided into smaller blocks. Conditional queries can then be used to assign each block an index value which conditionally ranks each block's favorability, from 0-5 with 5 being most favorable, for each model (e.g., lithologic, seismic, magnetic, structural). The values between 0-5 assigned to each block are referred to as index values. The final step of the process is to combine all the index models to create a favorability index. This involves multiplying each index model by a given weight and then summing the resulting values. The DEEPEN PFA Methodology follows this approach, but split up by the specific geologic components of each play type. These components are defined as follows for each magmatic play type: 1. Conventional hydrothermal plays in magmatic environments: Heat, fluid, and permeability 2. Superhot EGS plays: Heat, thermal insulation, and producibility (the ability to create and sustain fractures suitable for and EGS reservoir) 3. Supercritical plays: Heat, supercritical fluid, pressure seal, and producibility (the proper permeability and pressure conditions to allow production of supercritical fluid) More information on these components and their development can be found in Kolker et al., 2022. For the purposes of subsurface imaging, it is easier to detect a permeable fluid-filled reservoir than it is to detect separate fluid and permeability components. Therefore, in this analysis, we combine fluid and permeability for conventional hydrothermal plays, and supercritical fluid and producibility for supercritical plays. More information on this process is described in the following sections. We also project the 3D favorability volumes onto 2D surfaces for simplified joint interpretation, and we incorporate an uncertainty component. Uncertainty was modeled using the best approach for the dataset in question, for the datasets where we had enough information to do so. Identifying which subsurface parameters are the least resolved can help qualify current PFA results and focus future efforts in data collection. Where possible, the resulting uncertainty models/indices were weighted using the same weights applied to the respective datasets, and summed, following the PFA methodology above, but for uncertainty. There are two different versions of the Leapfrog model and associated favorability models: - v1.0: The first release in June 2023 - v2.1: The second release, with improvements made to the earthquake catalog (included additional identified events, removed duplicate events), to the temperature model (fixed a deep BHT), and to the index models (updated the seismicity-heat source index models for supercritical and EGS, and the resistivity-insulation index models for all three play types). Also uses the jet color map rather than the magma color map for improved interpretability. - v2.1.1: Updated to include v2.0 uncertainty results (see below for uncertainty model versions) There are two different versions of the associated uncertainty models: - v1.0: The first release in June 2023 - v2.0: The se...

15 GEOTHERMAL ENERGY↗

Newberry Volcano EGS Demonstration Stimulation Modeling

As a part of Phase I of the Newberry Volcano EGS Demonstration project, several data sets were collected to characterize the rock volume around the well. Fracture, fault, stress, and seismicity data has been collected by borehole televiewer, LiDAR elevation maps, and microseismic monitoring. Well logs and cuttings from the target well (NWG 55-29) and core from a nearby core hole (USGS N-2) have been analyzed to develop geothermal, geochemical, mineralogical and strength models of the rock matrix, altered zones, and fracture fillings (see Osborn et al., this volume). These characterization data sets provide inputs to models used to plan and predict EGS reservoir creation and productivity. One model used is AltaStim, a stochastic fracture and flow software model developed by AltaRock. The software's purpose is to model and visualize EGS stimulation scenarios and provide guidance for final planning. The process of creating an AltaStim model requires synthesis of geologic observations at the well, the modeled stress conditions, and the stimulation plan. Any geomechanical model of an EGS stimulation will require many assumptions and unknowns; thus, the model developed here should not be considered a definitive prediction, but a plausible outcome given reasonable assumptions. AltaStim is a tool for understanding the effect of known constraints, assumptions, and conceptual models on plausible outcomes.

15 GEOTHERMAL ENERGY↗

Joint impedance and facies inversion of time-lapse seismic data for improving monitoring of CO 2 incidentally stored from CO 2 EOR

Time-lapse seismic monitoring is an effective and proven technology for mapping the distribution of CO 2 in a subsurface reservoir. When injected CO 2 displaces other reservoir fluids, porous-medium properties are changed and thus the seismic impedance changes, causing time-lapse seismic amplitude differences in the injection zones. The analysis and interpretation of images created from these amplitude differences can provide information about reservoir architecture and the CO 2 migration within the reservoir. Incorporating seismic inversion and rock physics into the interpretation of time-lapse seismic data can considerably improve the modeling and monitoring to detect and assess the location of CO 2 over time. The joint inversion method presented in this paper has an integral representation of the geology in the inversion algorithm using elastic facies, which provides information about the spatial distribution of the geologic heterogeneities controlling the movement of fluids in the reservoir. The method was successfully applied to time-lapse seismic data from a mature oil field undergoing CO 2 enhanced oil recovery. The estimated seismic acoustic impedances and facies reflect the characteristics of individual geologic facies and fluid conditions of the reservoir subject to CO 2 injection. Finally, the probabilities estimated by the joint impedance and facies inversion for the reservoir's litho-fluid facies can be used for forecasting CO 2 saturation and pressure changes within the target reservoir.

4D seismic↗

Three‐Phase CO 2 Flow in a Basalt Fracture Network

Abstract Geologic CO 2 sequestration in basalt reservoirs is predicated on permanent CO 2 isolation via rapid mineralization reactions. This process is supported by a substantial body of evidence, including laboratory experiments documenting rapid mineralization rates, regional storage estimates indicating large, accessible storage reservoirs, and two successful pilot‐scale studies. Nevertheless, there remains significant uncertainty in the behavior of CO 2 flow within basalt fracture networks, particularly in the context estimating physical trapping potential in early time and as CO 2 undergoes phase change. In this study, a Monte Carlo numerical model is designed to simulate a supercritical CO 2 plume infiltrating a low‐permeability flood basalt entablature. The fracture network model is based on outcrop‐scale LiDAR mapping of Columbia River Basalt, and CO 2 flow is simulated within fifty equally probable realizations of the fracture network. The spatial distribution of fracture permeability for each realization is randomly drawn from a basalt aperture distribution, and ensemble results are analyzed with e‐type estimates to compute mean and standard deviation of fluid pressure and CO 2 saturation. Results of this model after 10 years of simulation suggest that (1) CO 2 flow converges on a single dominant flow path, (2) CO 2 accumulates at fracture intersections, and (3) variability in permeability can account for a 1.6 m depth interval within which free CO 2 may change phase from supercritical fluid to subcritical liquid or gas. In the context of CO 2 sequestration in basalt, these results suggest that physical CO 2 trapping may be substantially enhanced as carbonate minerals precipitate within the basalt fracture network.

Gierzynski, Alec O.↗

Overview of Observed and Measured Physical and Geophysical Effects of Historic Underground Nuclear Tests at the Nevada National Security Site

As detailed in a variety of sources (e.g., Allen et al. [1997]; Grasso [2001]), underground explosions can leave visible evidence of their occurrence in the geologic medium. The appearance and extent of the effects visible at the surface depend on a variety of factors, including scaled depth of burial (SDOB; actual depth of burial divided by the cube root of the yield of the explosion) and geologic medium characteristics. Surface effects were routinely documented for most of the underground nuclear explosions (UNEs) conducted at the former Nevada Test Site (now Nevada National Nuclear Site) and spatial data for them are available to the public (Grasso, 2001). Other effects, including collapse chimneys, glass-lined surfaces indicative of the initial explosion cavity, and shock damage to the medium, have also been characterized in various ways. These studies were done for a variety of purposes such as determining how closely UNEs could safely be positioned to conserve testing “real estate,” as input for containment design, to obtain survivability data for various types of equipment and instrumentation, and to understand device performance. Here we present examples of various mapped and measured explosion effects for the two main types of historical UNEs conducted at the NNSS, those in vertical shafts and those in horizontal tunnels. We do not include here cratering tests, which were shots conducted at shallow SDOBs for the purpose of creating a surface excavation.

explosion effects, Nevada National Security Site↗

Early Holocene cold snaps and their expression in the moraine record of the eastern European Alps

Glaciers preserve climate variations in their geological and geomorphological records, which makes them prime candidates for climate reconstructions. Investigating the glacier–climate system over the past millennia is particularly relevant first because the amplitude and frequency of natural climate variability during the Holocene provides the climatic context against which modern, human-induced climate change must be assessed. Second, the transition from the last glacial to the current interglacial promises important insights into the climate system during warming, which is of particular interest with respect to ongoing climate change. Evidence of stable ice margin positions that record cooling during the past 12 kyr are preserved in two glaciated valleys of the Silvretta Massif in the eastern European Alps, the Jamtal (JAM) and the Laraintal (LAR). We mapped and dated moraines in these catchments including historical ridges using beryllium-10 surface exposure dating (10Be SED) techniques and correlate resulting moraine formation intervals with climate proxy records to evaluate the spatial and temporal scale of these cold phases. The new geochronologies indicate the formation of moraines during the early Holocene (EH), ca. 11.0 ± 0.7 ka (n = 19). Boulder ages along historical moraines (n = 6) suggest at least two glacier advances during the Little Ice Age (LIA; ca. 1250–1850 CE) around 1300 CE and in the second half of the 18th century. An earlier advance to the same position may have occurred around 500 CE. The Jamtal and Laraintal moraine chronologies provide evidence that millennial-scale EH warming was superimposed by centennial-scale cooling. The timing of EH moraine formation coincides with brief temperature drops identified in local and regional paleoproxy records, most prominently with the Preboreal Oscillation (PBO) and is consistent with moraine deposition in other catchments in the European Alps and in the Arctic region. This consistency points to cooling beyond the local scale and therefore a regional or even hemispheric climate driver. Freshwater input sourced from the Laurentide Ice Sheet (LIS), which changed circulation patterns in the North Atlantic, is a plausible explanation for EH cooling and moraine formation in the Nordic region and in Europe.

54 ENVIRONMENTAL SCIENCES↗