Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “seismic data”

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

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

At least 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A Fusion of Geothermal and InSAR Data with Machine Learning for Enhanced Deformation Forecasting at the Geysers

The Geysers geothermal field in California is experiencing land subsidence due to the seismic and geothermal activities taking place. This poses a risk not only to the underlying infrastructure but also to the groundwater level which would reduce the water availability for the local community. Because of this, it is crucial to monitor and assess the surface deformation occurring and adjust geothermal operations accordingly. In this study, we examine the correlation between the geothermal injection and production rates as well as the seismic activity in the area, and we show the high correlation between the injection rate and the number of earthquakes. This motivates the use of this data in a machine learning model that would predict future deformation maps. First, we build a model that uses interferometric synthetic aperture radar (InSAR) images that have been processed and turned into a deformation time series using LiCSBAS, an open-source InSAR time series package, and evaluate the performance against a linear baseline model. The model includes both convolutional neural network (CNN) layers as well as long short-term memory (LSTM) layers and is able to improve upon the baseline model based on a mean squared error metric. Then, after getting preprocessed, we incorporate the geothermal data by adding them as additional inputs to the model. This new model was able to outperform both the baseline and the previous version of the model that uses only InSAR data, motivating the use of machine learning models as well as geothermal data in assessing and predicting future deformation at The Geysers as part of hazard mitigation models which would then be used as fundamental tools for informed decision making when it comes to adjusting geothermal operations.

Yazbeck, Joe (ORCID:0000000302235260)↗
DOE OSTI2025

Improving the Transportability of a Deep Learning Denoising Model Using Transfer Learning Techniques

The adoption of machine learning techniques in the seismology community has led to great performance improvements in several areas, including signal processing. Specifically, the development of deep learning–based seismic waveform denoising models has the potential to yield improvements in signal detection capabilities for networks operating in particularly noisy environments. Recent advancements in the design of these deep learning denoising models have included the incorporation of continuous and discrete wavelet transform functions into the network architecture to improve the learning capabilities and efficiency of said models. These wavelet transform–based seismic denoising models have shown improved denoising capabilities in regions where there is good agreement between the data features present in the training and evaluation datasets. However, questions remain about the overall transportability of these models to other monitoring regions. Here, in this study, we will determine the baseline transportability of a newly developed multilevel wavelet‐transform convolutional neural network (MWCNN) seismic denoising model. We accomplish this by taking a version of the MWCNN denoising model trained on data collected from the Utah region and evaluating its denoising performance on datasets collected from the neighboring Nevada region, which differ with regard to monitoring sensor types and event histories. We find that there is a notable variability in denoising performance related to the degree of similarity between the initial and new target datasets. The most notable difference in denoising performance is the ability of the denoising model to preserve accurate amplitude information associated with the signal energy present in the waveform data. Finally, we evaluate the ability of transfer learning techniques to improve the transportability of the MWCNN denoising model. We find that although there is still a performance gap present in the denoising results of the MWCNN model, transfer learning did yield improved results.

Quinones, Louis [Sandia National Laboratories (SNL↗
DOE OSTI2020

Play Fairway Analysis: Structurally Controlled Geothermal Systems in the Eastern Great Basin Extensional Regime, Utah

A research team with membership from the University of Utah/Energy & Geoscience Institute, the University of Utah/Dept. of Geology & Geophysics, and the Utah Geological Survey, undertook a play fairway analysis (PFA) for geothermal resources in the Eastern Great Basin (EGB) extensional tectonic regime of western Utah. This is a high-priority region for geothermal exploration because active Basin and Range (B&R) extension with volcanism having a N-S strike is superimposed upon pre-existing E-W belts of plutonic rocks and large-scale structural lineaments. Cumulative heat flow along the N-S strike of the state totals approximately 5 GWt above background stable interior. Three electricity producing power plants currently exist with substantial potential for increase. Succinctly, our PFA approach aims to resolve potential sources of heat and permeability in the region, which are the two principal criteria for establishing a geothermal resource. An initial Phase 1 was carried out using only existing geoscientific data in the area. Criteria selected for focusing heat potential include direct heat flow measurements in boreholes, magnetotelluric (MT) low resistivity anomalies, fluid/gas geochemistry, and proximity to recent volcanic eruptions. Permeability is established through geological structures (fault density, critically stressed areas, seismicity, gravity), and MT low resistivity anomalies. In Phase II of this PFA project, additional geological, geophysical and geochemical data were acquired and analysis carried out primarily over promising composite common risk (CCR) areas initially identified in Phase I in order to focus prospectivity and prepare for drilling recommendations. These prospects are near the Twin Peaks rhyolite field, high heat flow areas north of the producing Cove Fort system, and geophysical structure beneath the Crater Knoll area off the northeast flank of the Mineral Mountains. Additional data included MT site fill-in, structural mapping and analysis using high-resolution imagery, gravity and on-ground mapping, Nodal 3C passive seismic collection, and passive 3He surveying. Heat source and permeability potential are again expressed in terms of their individual common risk segment (CRS) maps, with a color scheme using green for most favorable (low risk) and red for least favorable (high risk). Diverse data types are united through the technique of probability kriging, which establishes prospectivity thresholds for each data type and then computes probability of exceeding that threshold over the PFA area. Modified CRS and CCRS maps are compared to those of Phase I to highlight tighter prospectivity focus. In doing so, the prospectivity threshold for heat was increased significantly to narrow the targeting. In the final Phase III of this project where a recommended deep thermal gradient hole was sited, additional geophysical, geochemical and geological field collection and analysis was carried out to refine drill hole targeting. This includes prospect-scale MT, gravity, structure, passive seismic deployment (Cove Fort area), and detailed 3He isotope profiling. It was the recommendation of the DOE Technical Monitoring Team (TMT) that one or more holes be sited in the north-ern Cove Fort area where legacy TG gradient holes showed high cumulative heat flow. These were to be of moderate depth, 2000-3000 feet, to reach the geothermal fluid table expected to start in excess of 1000 feet depth. The drilling organization stipulated by the DOE/GTO was that of the USGS Research Drilling Program (RDP) centered in Las Vegas, NV. A detailed well plan, appended to this report, was developed principally by Dr. Ben Barker consulting to University of Utah, Dr. Steve Pye on the DOE TMT, Mr. Steven Crawford of the USGS-RDP, and the project PI Phil Wannamaker. However, temperature and possible H 2 S at the systems lead to cancellation of the drilling last-minute as this appeared outside the experience base of the USGS-RDP. We hope to have the opportunity to revisit the test drilling and expand the Play Fairway Analysis of this region at some point in the future.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Aquifer Monitoring Using Ambient Seismic Noise Recorded With Distributed Acoustic Sensing (DAS) Deployed on Dark Fiber

Abstract Groundwater is a critical resource for human activities worldwide, and a vital component of many natural ecosystems. However, the state and dynamics of water‐bearing aquifers remain uncertain, mostly due to the paucity of subsurface data at high spatial and temporal resolution. Here, we show that analysis of infrastructure‐generated ambient seismic noise acquired on distributed acoustic sensing (DAS) arrays has potential as a tool to track variations in seismic velocities (d v / v ) caused by groundwater level fluctuations. We analyze 5 months of ambient noise acquired along an unused, 23 km‐long telecommunication fiber‐optic cable in the Sacramento Valley, CA, a so‐called “dark fiber." Three array subsections, ∼6 km apart, are processed and the stretching technique is applied to retrieve daily d v / v beneath each location. Near the Sacramento river, d v / v variations in the order of 2%–3% correlate with precipitation events and fluctuations in river stage of ∼1.5 m. In contrast, regions away (2.5 km) from the river do not experience large d v / v variations. These observations reveal short‐scale spatial variability in aquifer dynamics captured by this approach. Dispersion analysis and surface wave inversion of noise gathers reveal that seismic velocity perturbations occur at depths of 10–30 m. Rock physics modeling confirms that observed d v / v are linked to pore pressure changes at these depths, caused by groundwater table fluctuations. Our results suggest that DAS combined with ambient noise interferometry provides a means of tracking aquifer dynamics at high spatial and temporal resolutions at local to regional scales, relevant for effective groundwater resource management.

58 GEOSCIENCES↗
DOE OSTI2024

Gulf of Mexico Miocene CO2 Site Characterization Mega Transect

The project Gulf of Mexico Miocene CO2 Site Characterization Mega Transect is one of 9 site characterization projects that were implemented as part of ARRA (American Recovery and Reinvestment Act). Data from this project was used to improve resolution of data in NATCARB in the area of study. Data related to this study has already been incorporated in NATCARB Atlas. The University of Texas at Austin study focused on (1) collecting high-resolution, 3-D (HR3D) seismic datasets of the shallow geologic section to evaluate the capacity of potential reservoirs and the sealing capacity of the confining system, and (2) producing a regional CO2 atlas of Miocene age units of the upper Texas coast submerged lands (southeast Texas). The goal of the study was to evaluate the capability of the Miocene age geologic section of Texas submerged lands to permanently store large volumes of anthropogenic CO2. The geographic locations of three HR3D, “P-Cable” seismic surveys collected in the near-shore waters of southeast Texas offers a vertical and horizontal view (data “cube”) example from the second survey. The dataset's high quality allowed for defining recognizable geologic morphologies (shapes) and possible fluid content, which are important for determining an area’s prospects for geologic storage. The regional CO2 atlas of Miocene age units highlights several topics, including: (1) a regional analyses of petroleum systems as analogs for CO2 storage; (2) petrography, petrology, and extent of confining systems of the Miocene section; (3) static capacity estimate; and (4) examples of characterization methodologies of prospect areas. The HR3D datasets and regional analyses within the CO2 atlas of the offshore Texas Miocene provide a sound basis for a future generation of specific CO2 prospects in the study area.

.las↗
DOE OSTI2025

MTUQ: a framework for estimating moment tensors, point forces, and their uncertainties

SUMMARY We introduce MTUQ, an open-source Python package for seismic source estimation and uncertainty quantification, emphasizing flexibility and operational scalability. MTUQ provides MPI-parallelized grid search and global optimization capabilities, compatibility with 1-D and 3-D Green’s function database formats, customizable data processing, C-accelerated waveform and first-motion polarity misfit functions, and utilities for plotting seismic waveforms and visualizing misfit and likelihood surfaces. Applicability to a range of full- and constrained-moment tensor, point force, and centroid inversion problems is possible via a documented application programming interface, accompanied by example scripts and integration tests. We demonstrate the software using three different types of seismic events: (1) a 2009 intraslab earthquake near Anchorage, Alaska; (2) an episode of the 2021 Barry Arm landslide in Alaska; and (3) the 2017 Democratic People’s Republic of Korea underground nuclear test. With these events, we illustrate the well-known complementary character of body waves, surface waves, and polarities for constraining source parameters. We also convey the distinct misfit patterns that arise from each individual data type, the importance of uncertainty quantification for detecting multimodal or otherwise poorly constrained solutions, and the software’s flexible, modular design.

58 GEOSCIENCES↗
DOE OSTI2022

Optimizing Earthquake Nowcasting With Machine Learning: The Role of Strain Hardening in the Earthquake Cycle

Abstract Nowcasting is a term originating from economics, finance, and meteorology. It refers to the process of determining the uncertain state of the economy, markets or the weather at the current time by indirect means. In this paper, we describe a simple two‐parameter data analysis that reveals hidden order in otherwise seemingly chaotic earthquake seismicity. One of these parameters relates to a mechanism of seismic quiescence arising from the physics of strain‐hardening of the crust prior to major events. We observe an earthquake cycle associated with major earthquakes in California, similar to what has long been postulated. An estimate of the earthquake hazard revealed by this state variable time series can be optimized by the use of machine learning in the form of the Receiver Operating Characteristic skill score. The ROC skill is used here as a loss function in a supervised learning mode. Our analysis is conducted in the region of 5° × 5° in latitude‐longitude centered on Los Angeles, a region which we used in previous papers to build similar time series using more involved methods (Rundle & Donnellan, 2020, https://doi.org/10.1029/2020EA001097 ; Rundle, Donnellan et al., 2021, https://doi.org/10.1029/2021EA001757 ; Rundle, Stein et al., 2021, https://doi.org/10.1088/1361-6633/abf893 ). Here we show that not only does the state variable time series have forecast skill, the associated spatial probability densities have skill as well. In addition, use of the standard ROC and Precision (PPV) metrics allow probabilities of current earthquake hazard to be defined in a simple, straightforward, and rigorous way.

58 GEOSCIENCES↗
DOE OSTI2024

Enhancing Regional Seismic Velocity Models With Higher‐Resolution Local Results Using Sparse Dictionary Learning

Abstract We use sparse dictionary learning to develop transformations between seismic velocity models of different resolution and spatial extent. Starting with data in the common region of both models, the method can enhance a regional lower‐resolution model to match the style and resolution of local higher‐resolution results while preserving its regional coverage. The method is demonstrated by applying it to two‐dimensional V S and three‐dimensional V P and V S regional and local velocity models in southern California. The enhanced reconstructed regional results exhibit clear visual improvements, especially in the reconstructed V P / V S ratios, and better correlations with geological features. Moreover, the reconstructed regional V P , V S models outperform the original ones in comparison of simulated earthquake waveforms to observations. The improved fitting to observed waveforms extends beyond the domain of the overlapping region. The developed dictionary learning approach provides physically interpretable results and offers a powerful tool for additional applications of data enhancement in earth sciences.

Geochemistry & Geophysics↗
DOE OSTI2023

PubDAS: A PUBlic Distributed Acoustic Sensing Datasets Repository for Geosciences

During the past few years, distributed acoustic sensing (DAS) has become an invaluable tool for recording high-fidelity seismic wavefields with great spatiotemporal resolutions. However, the considerable amount of data generated during DAS experiments limits their distribution with the broader scientific community. Such a bottleneck inherently slows down the pursuit of new scientific discoveries in geosciences. Here, we introduce PubDAS—the first large-scale open-source repository where several DAS datasets from multiple experiments are publicly shared. PubDAS currently hosts eight datasets covering a variety of geological settings (e.g., urban centers, underground mines, and seafloor), spanning from several days to several years, offering both continuous and triggered active source recordings, and totaling up to ~90 TB of data. Here this article describes these datasets, their metadata, and how to access and download them. Some of these datasets have only been shallowly explored, leaving the door open for new discoveries in Earth sciences and beyond.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Mapping Glacier Basal Sliding Applying Machine Learning

During the RESOLVE project (“High-resolution imaging in subsurface geophysics: development of a multi-instrument platform for interdisciplinary research”), continuous surface displacement and seismic array observations were obtained on Glacier d’Argentière in the French Alps for 35 days in May 2018. The data set is used to perform a detailed study of targeted processes within the highly dynamic cryospheric environment. In particular, the physical processes controlling glacial basal motion are poorly understood and remain challenging to observe directly. Especially in the Alpine region for temperate based glaciers where the ice rapidly responds to changing climatic conditions and thus, processes are strongly intermittent in time and heterogeneous in space. Spatially dense seismic and Global Positioning System (GPS) measurements are analyzed applying machine learning to gain insight into the processes controlling glacial motions of Glacier d’Argentière. Using multiple bandpass-filtered copies of the continuous seismic waveforms, we compute energy-based features, develop a matched field beamforming catalog and include meteorological observations. Features describing the data are analyzed with a gradient boosting decision tree model to directly estimate the GPS displacements from the seismic noise. We posit that features of the seismic noise provide direct access to the dominant parameters that drive displacement on the highly variable and unsteady surface of the glacier. The machine learning model infers daily fluctuations and longer term trends. The results show on-ice displacement rates are strongly modulated by activity at the base of the glacier. The techniques presented provide a new approach to study glacial basal sliding and discover its full complexity.

58 GEOSCIENCES↗
DOE OSTI2025

Evaluation of Station Performance of the Idaho National Laboratory Seismic Monitoring Network Using Network Detection Thresholds

The Idaho National Laboratory (INL) Seismic Monitoring Network is located in eastern Idaho and monitors a portion of the intermountain seismic belt. It has been in place for 50 yr and has undergone several major changes, the most recent of which has been the transition to the Antelope real‐time acquisition system and the implementation of automatic phase picking algorithms to aid in analysis. This study discusses the efforts to evaluate the performance of the INL seismic monitoring network (and other surrounding stations) using the new real‐time acquisition system. The method outlined by Wilson et al. (2021) is used to develop an empirical relationship between the observability of local earthquakes as a function of magnitude and distance. This relationship is used to produce detection thresholds for Pwaves for all stations of interest. The INL seismic network has two main goals: monitor tectonic‐and volcanic‐related events and measure ground motions for input into seismic hazard analysis. Because of these two overall objectives, several seismic stations have been installed near critical facilities and, therefore, are not as quiet as stations that are used primarily for earthquake detection. This is reflected in their detection thresholds, which are much smaller for stations away from facilities. This study shows that the INL Seismic Monitoring Network is able to detect earthquakes near INL facilities with M L > 1.2, with redundancies built in to ensure this sensitivity even if data became unavailable from some stations. This study also shows “holes” in the monitoring network where the detection of smaller earthquakes is highly dependent on sparsely placed seismic stations. In conclusion, the results of this study will be used to govern plans for expansion of earthquake monitoring in Idaho and the surrounding region and to fine‐tune the detection thresholds for individual stations.

58 - GEOSCIENCES↗
DOE OSTI2022

The Edge of Exploration: An Edge Storage and Computing Framework for Ambient Noise Seismic Interferometry Using Internet of Things Based Sensor Networks

Recent technological advances have reduced the complexity and cost of developing sensor networks for remote environmental monitoring. However, the challenges of acquiring, transmitting, storing, and processing remote environmental data remain significant. The transmission of large volumes of sensor data to a centralized location (i.e., the cloud) burdens network resources, introduces latency and jitter, and can ultimately impact user experience. Edge computing has emerged as a paradigm in which substantial storage and computing resources are located at the “edge” of the network. In this paper, we present an edge storage and computing framework leveraging commercially available components organized in a tiered architecture and arranged in a hub-and-spoke topology. The framework includes a popular distributed database to support the acquisition, transmission, storage, and processing of Internet-of-Things-based sensor network data in a field setting. We present details regarding the architecture, distributed database, embedded systems, and topology used to implement an edge-based solution. Lastly, a real-world case study (i.e., seismic) is presented that leverages the edge storage and computing framework to acquire, transmit, store, and process millions of samples of data per hour.

58 GEOSCIENCES↗
DOE OSTI2022

Unsupervised learning from three-component accelerometer data to monitor the spatiotemporal evolution of meso-scale hydraulic fractures

Enhanced geothermal systems can provide a substantial share of the global energy demand. There exist several hurdles in the engineering implementations of such geothermal systems. One such hurdle is the accurate monitoring of the fracture networks created in subsurface through hydraulic stimulation of these systems. Micro seismicity associated with the stimulation is the primary means to locate the event hypocenters for estimating the stimulated rock volume. Existing methods for location the hypocenters are restricted to only the highest amplitude impulsive signals that are simultaneously detected on several sensors. Consequently, a large portion (usually ~99%) of the measurements are left unused. In this paper, an unsupervised manifold-approximation followed by clustering of 3-component accelerometer data is used to analyze the seismicity recorded on a monitoring well. With this method, a larger portion of the measured signal is used for the monitoring of the hydraulic fracture network. We analyze the EGS Collab experiment 1 microseismic data, recorded at the Sanford Underground Research Facility, South Dakota. Using the data from a single three-component accelerometer, the polarization features viz. Azimuth, incidence, rectilinearity, and planarity are used as inputs for the unsupervised manifold approximation followed by clustering. Our study shows that density-based clusters in the projected 3D space correspond to distinct types of hydraulically fractured zones around the injection point. Finally, we show that the temporal evolution of these clusters can be used to track fracture creation and propagation.

58 GEOSCIENCES↗
DOE OSTI2023

The effect of shear strain and shear localization on fault healing

The seismic cycle of repeated earthquake failure requires that faults regain frictional strength during the interseismic phase, when the fault is locked or undergoing quasi-static creep. Fault healing plays a central role in determining earthquake stress drop, recurrence interval, elastic radiation frequency and other source parameters. In particular, the longer a fault remains quasi-stationary, the stronger it becomes and the larger the potential stress drop can be for the next event. Here, we address the role of shear strain and strain localization on fault healing and healing rate. We performed slide-hold-slide friction experiments on quartz gouge in the double-direct shear configuration for shear strain up to 25 and hold times from 10 to 1000 s. The results show that both healing and healing rate increase nonlinearly with increasing shear strain. Frictional healing scales with volumetric strain within the laboratory fault zone. Using the volumetric strain upon reshear as a proxy for strain localization, we demonstrate that the capacity of a fault to heal is directly proportional to shear bandwidth and degree of strain localization. The more the deformation is localized, the higher are the healing and healing rate, and thus, the fault strength. Our data provide a framework for understanding variations in fault strength over the seismic cycle and the role of brecciation and strain localization on spatiotemporal variations in fault strength.

58 GEOSCIENCES↗
DOE OSTI2023

Seismoacoustic Explosion Yield and Depth Estimation: Insights from the Large Surface Explosion Coupling Experiment

The Large Surface Explosion Coupling Experiment (LSECE) is a chemical explosion experiment conducted in Yucca Flat at the Nevada National Security Site in 2020. The experiment included two surface detonations of ~1000 kg trinitrotoluene equivalent. The main goal of this experiment was to provide the ground-truth data for seismoacoustic wave excitation by large chemical explosions near the ground surface. The seismic and acoustic energy partitioning between the surface is strongly governed by the depth or height of explosions, and either seismic or acoustic-only analysis may have inherent ambiguity in determining explosion yield and depth simultaneously. Previous studies suggested that joint seismoacoustic analysis can resolve the trade-off and reduce the uncertainty of yield and depth estimation dramatically. Here we demonstrate the capability of seismoacoustic analysis to improve the accuracy of explosion yield and depth estimation with the LSECE data. Local acoustic wave propagation in the atmosphere can be substantially affected by constantly varying weather conditions. Consisting of two detonations before dawn and in the afternoon, LSECE provides unique data to evaluate the model accuracy of acoustic wave propagation and seismoacoustic energy partitioning depending on local atmospheric conditions. We quantitatively evaluate the accuracy of yield and depth estimation depending on atmospheric variability and the improvement achieved by the joint seismoacoustic approach.

58 GEOSCIENCES↗
DOE OSTI2022

Post-hazard condition assessment of nuclear piping-equipment systems: Novel approach to feature extraction and deep learning

Over the past decade, the use of artificial intelligence techniques in the field of health-monitoring has gained significant interest, especially for structures such as building and bridges. However, applications to industrial systems such as equipment-piping systems in nuclear plants have not been explored. Here, in this paper, it is shown that the existing techniques developed for buildings and bridges cannot be extended directly to equipment-piping systems as the response of such systems is governed by multiple localized modes unlike that in buildings and bridges. This paper proposes a new approach that consists of three key aspects: (i) a novel vector of degradation-sensitive features extracted from measured data, (ii) using a deep Artificial Neural Network (ANN) for diagnosis of degradation location and degradation severity, and (iii) consideration of uncertainty in degradation severity when training the ANN. Degradation in piping-equipment systems can occur due to flow-accelerated erosion and corrosion. These locations can potentially exhibit damage such as localized yielding or initiation of cracking due to an external event such as an earthquake. Moreover, such locations can at times go undetected by current inspection techniques. Therefore, a robust framework is needed for detection of degradation after a seismic event. This manuscript proposes a proof-of-concept framework, which utilizes data collected from sensors to generate a deep ANN database for predicting degraded locations and severity in a piping-equipment system. Degradation severity is classified as minor, moderate, and severe. In the suggested methodology, a novel vector of degradation-sensitive features is extracted from the sensor data to train the ANN. A simple piping-equipment system is selected to demonstrate feature extraction as a means to simplify pattern recognition, explore the design and parameters of an ANN, and develop a sensor placement strategy. The effectiveness of the proposed framework is demonstrated on a realistic primary safety system of a two-loop nuclear reactor. It is shown that the proposed post-hazard condition assessment framework is able to detect degraded locations along with the severity levels, including minor degradation, with considerably higher accuracy.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Physical properties, internal structure, and the three‐dimensional petrography of CI chondrites

physical properties and the nature of their breccation, we investigated nine samples of the Ivuna and Orgueil CI chondrites ranging in size from 1 mm to 4 cm in approximate diameter. The combined mass of unique material investigated in this work is 113 g. For our investigations, we use ideal gas pycnometry, 3-D laser scanning, x-ray computed microtomography (μCT), and accompanying digital data extraction techniques. We found that the bulk density of the samples ranged from 1.61 to 2.10 g cm −3 . Larger samples tend to have a lower bulk density. Grain density (ranging from 2.44 to 2.55 g cm −3 ) is significantly less variable than the bulk density in our samples and the quantity of porosity (ranging from 14.6% to 33.8%) is the dominant factor in determining the bulk density of CI chondrite material. Our μCT results show that the visible porosity across all sizes of our CI chondrite samples is in the form of cracks, but these cracks can account for less than two-thirds of the porosity in the CI chondrites. Other porosity is not visible, even at μCT resolutions of 2.7 μm voxel edge −1 and we conclude that it is sub-micron in nature. It is not clear if the cracks seen in our samples are indigenous to the chondrites or are a result of terrestrial processes. We also find that the CI chondrites are excellent examples of the fractal-like nature of brecciation, where clasts can be observed at all scales we imaged. The breccias are composed of sub-equant-shaped and sub-rounded-textured clasts like melt-free impact breccias on other solar system bodies. From our μCT volume and digital data extraction, we determine that the Ivuna CI chondrite breccia is organized: the mostly sub-equant clasts within our ~2 cm chunk of Ivuna have a mean diameter of 1.33 mm and their aligned longest axes define a lineation structure. We speculate that the lineation was imparted after fragmentation of the clasts by slight shear on the parent asteroid which could be the result of seismic-related granular flow or mild non-axial impact-related compaction. These data will help to place returned asteroidal material from asteroids 162173 Ryugu and 101955 Bennu and the CI chondrites into a mutual geological context.

CI chondrite↗
DOE OSTI2021

Evaluating Probability of Containment Effectiveness at a GCS Sites using integrated assessment modeling approach with Bayesian decision Networks

Improved scientific and engineering understanding of the behavior of geologic CO2 storage together with established regulatory framework and incentive structures raise the prospects for accelerated, large-scale deployment of this greenhouse gas emissions reduction approach. Incentive structures call for the establishment of appropriate verification and accounting approaches to support claims of the integrity of a geologic storage complex and to justify taking credit for long-term storage. In this study, we present a framework for assessing the probability of containment effectiveness over the lifetime of a geologic carbon storage site (e.g., after 70 years of injection and post-injection site performance) using forward stochastic model realizations based on site characterization data and using a monitoring-informed Bayesian network based on hypothetical detectability from surface seismic surveys over the site injection and post-injection phases. The National Risk Assessment Partnership’s open-source Integrated Assessment Model (NRAP-Open-IAM) was utilized to develop an ensemble of 10,000 a priori stochastic forecasts of CO2 containment. Those simulations were used to train the Bayesian network model to estimate the prior probabilities of the CO2 leakage mass into overlying, monitorable aquifers considering the uncertainties in the reservoir properties, permeability of potentially leaky wells and the overlying aquifers. The conditional probabilities in the Bayesian network were either learned from the NRAP-Open-IAM simulations or derived from the predefined detection thresholds for the monitoring method. Observations obtained from monitoring, over time during the site operation phases were then used to generate updated posterior probabilities of containment (and any loss from containment) in the Bayesian network by propagating the prior probabilities through the conditional probabilities. We demonstrate how to construct and use the Bayesian network for verifying the long-term storage complex effectiveness informed by monitoring based on the NRAP-Open-IAM simulations previously developed for the FutureGen 2.0 site. This approach may have relevance for stake holders to demonstrate secure geologic storage, provide a defensible, probabilistic approach to claim credit for geologic storage, and to estimate the likelihood that any fraction of the claimed credit may need to be refunded to the creditor based on available monitoring information.

Bayesian network, Risk assessment, Monitoring, car↗