Search NASA⌕ Search

SEARCH · Search NASA

Results for “seismic data processing”

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

Revealing complex subsurface dynamics with continuous seismic monitoring: Observations using distributed acoustic sensing and surface orbital vibrators during hydraulic fracturing

Understanding hydraulic fracturing is crucial to improving the stimulation of unconventional reservoirs and increasing fluid production. This study develops a novel seismic monitoring technology using distributed acoustic sensing (DAS) and surface orbital vibrators (SOV) to capture fracture seismic response and mechanical properties at high temporal intervals. We analyze continuous time-lapse vertical seismic profiling (VSP) data acquired every hour during the first nine days of treatment of an unconventional reservoir in the Austin Chalk/Eagle Field Laboratory. The VSP data contain clear seismic signals scattered from the activated fractures. The spatiotemporal changes of the fracture reflectivity revealed by the SOV/DAS data correlate well with the observations of fracture locations inferred from low-frequency DAS data. These results capture the fracture opening and closure processes, as well as highlight potential prestage activations of the fractures due to hydraulic connectivity with preexisting fracture systems. Therefore, analysis of the presented data set provides a unique opportunity to understand fracture initiation and subsequent evolution, not only in the context of unconventional resources but also in enhanced geothermal systems.

Correa, Julia↗

Seismic Wave Coherence as a Function of Propagation Azimuth at the Large-N Array, Source Physics Experiment

The Large-N array of the Source Physics Experiment (SPE) consisted, in part, of 496 vertical component geophones that recorded the seismic wave field produced by the SPE-5 buried chemical explosion. Preliminary observations of the data showed a large degree of azimuthally dependent seismic scattering, particularly for post-P wave arrivals, hindering surface wave analysis. We document and quantify the azimuthal dependence of the wave field scattering in order to guide future coherent wave field processing methods. Specifically, we form three linear arrays, with different nominal source-receiver azimuths, by extracting a subset of the Large-N stations. For each linear array, we evaluate wave field coherence as a function of frequency and inter-station distance. For P waves, we observe that there is a strong azimuthal dependence of wave coherence, with the highest degree of scattering occurring in a northwest/southeast propagation direction. This suggests that there are structural elements beneath the Large-N array that affect the direct source to receiver body wave ray path. We also observe that the scattering of the post-P energy displays a coherence that is dependent on both frequency and azimuthal direction. This energy is preferentially coherent in the southwest-to-northeast propagation direction, consistent with the strike of the steeply dipping fault (Boundary fault) adjacent to the northeast side of the Large-N array, but only at low frequencies (<10 Hz). At higher frequencies, the azimuthally dependent wave coherence diminishes, suggesting that the scattering of high frequency portion of the post-P wave field is independent of the large-scale geologic structure at this site.

58 GEOSCIENCES↗

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

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

25 ENERGY STORAGE↗

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↗

Modeling CO 2 Sequestration in Saline Aquifer and Depleted Oil Reservoir To Evaluate Regional CO 2 Sequestration Potential of Ozark Plateau Aquifer System, South-Central Kansas

1. Drilled, cored, and logged three wells to the basement and collecting more than 2,700 ft of conventional core; obtained 20 m i2 of multicomponent 3D seismic imaging and merged and reprocessed more than 125 mi 2 of existing 3D seismic data for use in modeling CO 2 - EOR oil recovery and CO 2 storage in five oil fields in southern Kansas. 2. Determined the technical feasibility of injecting and sequestering CO 2 in a set of four depleted oil reservoirs in the Cutter, Pleasant Prairie South, Eubank, and Shuck fields in southwest Kansas; of concurrently recovering oil from those fields; and of quantifying the volumes of CO2 sequestered and oil recovered during the process. 3. Formed a consortium of six oil operating companies, five of which own and operate the four fields. The consortium became part of the Southwest Kansas CO 2 -EOR Initiative for the purpose of sharing data, knowledge, and interest in understanding the potential for CO 2 -EOR in Kansas. 4. Built a regional well database covering 30,000 mi 2 and containing stratigraphic tops from ~90,000 wells; correlated 30 major stratigraphic horizons; digitized key wells, including wireline logs and sample logs; and analyzed more than 3,000 drill stem tests to establish that fluid levels in deep aquifers below the Permian evaporites are not connected to the surface and therefore pressures are not hydrostatic. Connectivity with the surface aquifers is lacking because shale aquitards and impermeable evaporite layers consist of both halite and anhydrite. 5. Developed extensive web applications and an interactive mapping system that do the following: a. Facilitate access to a wide array of data obtained in the study, including core descriptions and analyses, sample logs, digital (LAS) well logs, seismic data, gravity and magnetics maps, structural and stratigraphic maps, inferred fault traces, earthquakes, Class I and II disposal wells, and surface lineaments. b. Provide real-time analysis of the project dataset, including automated integration and viewing of well logs, core, core analyses, brine chemistry, and stratigraphy using the Java Profile app. A cross-section app allows for the display of log data for up to four wells at a time. 6. Integrated interpretations from the project’s interactive web-based mapping system to gain insights to aid in assessing the efficacy of geologic CO 2 storage in Kansas and insights toward understanding recent seismicity to aid in evaluating induced vs. naturally occurring earthquakes. 7. Developed a digital type-log system, including web-based software to modify and refine stratigraphic nomenclature to provide stakeholders a common means for communication about the subsurface. 8. Contracted use of a nuclear magnetic resonance (NMR) log and ran it slowly to capture response and characterize larger pores common for carbonate reservoirs. Used NMR to extend core analyses to apply permeability, relative permeability to CO 2 , and capillary pressure to the major rock types, each uniquely expressed as a reservoir quality index (RQI), present in the Mississippian and Arbuckle rocks. 9. Characterized and evaluated the possible role of microbes in dense brines. Used microbes to compliment H/O stable isotopes to fingerprint brine systems. Used perforation/swabbing to obtain samples from multiple hydrostratigraphic units and confirmed equivalent results using less expensive drill stem tests (DST). 10. Used an integrated approach from whole core, logs, tests, and seismic to verify and quantify properties of vuggy, brecciated, and fractured carbonate intervals. 11. Used complex geocellular static and dynamic models to evaluate regional storage capacity using large parallel processing. 12. Carbonates are complex reservoirs and CO 2 -EOR needs to move to the next generation to increase effectiveness of CO 2 and efficiency and safety of the injection.

54 ENVIRONMENTAL SCIENCES↗

Modeling CO 2 Sequestration in Saline Aquifer and Depleted Oil Reservoir To Evaluate Regional CO 2 Sequestration Potential of Ozark Plateau Aquifer System, South-Central Kansas

1. Drilled, cored, and logged three wells to the basement and collecting more than 2,700 ft of conventional core; obtained 20 m i2 of multicomponent 3D seismic imaging and merged and reprocessed more than 125 mi 2 of existing 3D seismic data for use in modeling CO 2 - EOR oil recovery and CO 2 storage in five oil fields in southern Kansas. 2. Determined the technical feasibility of injecting and sequestering CO 2 in a set of four depleted oil reservoirs in the Cutter, Pleasant Prairie South, Eubank, and Shuck fields in southwest Kansas; of concurrently recovering oil from those fields; and of quantifying the volumes of CO2 sequestered and oil recovered during the process. 3. Formed a consortium of six oil operating companies, five of which own and operate the four fields. The consortium became part of the Southwest Kansas CO 2 -EOR Initiative for the purpose of sharing data, knowledge, and interest in understanding the potential for CO 2 -EOR in Kansas. 4. Built a regional well database covering 30,000 mi 2 and containing stratigraphic tops from ~90,000 wells; correlated 30 major stratigraphic horizons; digitized key wells, including wireline logs and sample logs; and analyzed more than 3,000 drill stem tests to establish that fluid levels in deep aquifers below the Permian evaporites are not connected to the surface and therefore pressures are not hydrostatic. Connectivity with the surface aquifers is lacking because shale aquitards and impermeable evaporite layers consist of both halite and anhydrite. 5. Developed extensive web applications and an interactive mapping system that do the following: a. Facilitate access to a wide array of data obtained in the study, including core descriptions and analyses, sample logs, digital (LAS) well logs, seismic data, gravity and magnetics maps, structural and stratigraphic maps, inferred fault traces, earthquakes, Class I and II disposal wells, and surface lineaments. b. Provide real-time analysis of the project dataset, including automated integration and viewing of well logs, core, core analyses, brine chemistry, and stratigraphy using the Java Profile app. A cross-section app allows for the display of log data for up to four wells at a time. 6. Integrated interpretations from the project’s interactive web-based mapping system to gain insights to aid in assessing the efficacy of geologic CO 2 storage in Kansas and insights toward understanding recent seismicity to aid in evaluating induced vs. naturally occurring earthquakes. 7. Developed a digital type-log system, including web-based software to modify and refine stratigraphic nomenclature to provide stakeholders a common means for communication about the subsurface. 8. Contracted use of a nuclear magnetic resonance (NMR) log and ran it slowly to capture response and characterize larger pores common for carbonate reservoirs. Used NMR to extend core analyses to apply permeability, relative permeability to CO 2 , and capillary pressure to the major rock types, each uniquely expressed as a reservoir quality index (RQI), present in the Mississippian and Arbuckle rocks. 9. Characterized and evaluated the possible role of microbes in dense brines. Used microbes to compliment H/O stable isotopes to fingerprint brine systems. Used perforation/swabbing to obtain samples from multiple hydrostratigraphic units and confirmed equivalent results using less expensive drill stem tests (DST). 10. Used an integrated approach from whole core, logs, tests, and seismic to verify and quantify properties of vuggy, brecciated, and fractured carbonate intervals. 11. Used complex geocellular static and dynamic models to evaluate regional storage capacity using large parallel processing. 12. Carbonates are complex reservoirs and CO 2 -EOR needs to move to the next generation to increase effectiveness of CO 2 and efficiency and safety of the injection.

54 ENVIRONMENTAL SCIENCES↗

Modeling CO 2 Sequestration in Saline Aquifer and Depleted Oil Reservoir To Evaluate Regional CO 2 Sequestration Potential of Ozark Plateau Aquifer System, South-Central Kansas

1. Drilled, cored, and logged three wells to the basement and collecting more than 2,700 ft of conventional core; obtained 20 m i2 of multicomponent 3D seismic imaging and merged and reprocessed more than 125 mi 2 of existing 3D seismic data for use in modeling CO 2 - EOR oil recovery and CO 2 storage in five oil fields in southern Kansas. 2. Determined the technical feasibility of injecting and sequestering CO 2 in a set of four depleted oil reservoirs in the Cutter, Pleasant Prairie South, Eubank, and Shuck fields in southwest Kansas; of concurrently recovering oil from those fields; and of quantifying the volumes of CO2 sequestered and oil recovered during the process. 3. Formed a consortium of six oil operating companies, five of which own and operate the four fields. The consortium became part of the Southwest Kansas CO 2 -EOR Initiative for the purpose of sharing data, knowledge, and interest in understanding the potential for CO 2 -EOR in Kansas. 4. Built a regional well database covering 30,000 mi 2 and containing stratigraphic tops from ~90,000 wells; correlated 30 major stratigraphic horizons; digitized key wells, including wireline logs and sample logs; and analyzed more than 3,000 drill stem tests to establish that fluid levels in deep aquifers below the Permian evaporites are not connected to the surface and therefore pressures are not hydrostatic. Connectivity with the surface aquifers is lacking because shale aquitards and impermeable evaporite layers consist of both halite and anhydrite. 5. Developed extensive web applications and an interactive mapping system that do the following: a. Facilitate access to a wide array of data obtained in the study, including core descriptions and analyses, sample logs, digital (LAS) well logs, seismic data, gravity and magnetics maps, structural and stratigraphic maps, inferred fault traces, earthquakes, Class I and II disposal wells, and surface lineaments. b. Provide real-time analysis of the project dataset, including automated integration and viewing of well logs, core, core analyses, brine chemistry, and stratigraphy using the Java Profile app. A cross-section app allows for the display of log data for up to four wells at a time. 6. Integrated interpretations from the project’s interactive web-based mapping system to gain insights to aid in assessing the efficacy of geologic CO 2 storage in Kansas and insights toward understanding recent seismicity to aid in evaluating induced vs. naturally occurring earthquakes. 7. Developed a digital type-log system, including web-based software to modify and refine stratigraphic nomenclature to provide stakeholders a common means for communication about the subsurface. 8. Contracted use of a nuclear magnetic resonance (NMR) log and ran it slowly to capture response and characterize larger pores common for carbonate reservoirs. Used NMR to extend core analyses to apply permeability, relative permeability to CO 2 , and capillary pressure to the major rock types, each uniquely expressed as a reservoir quality index (RQI), present in the Mississippian and Arbuckle rocks. 9. Characterized and evaluated the possible role of microbes in dense brines. Used microbes to compliment H/O stable isotopes to fingerprint brine systems. Used perforation/swabbing to obtain samples from multiple hydrostratigraphic units and confirmed equivalent results using less expensive drill stem tests (DST). 10. Used an integrated approach from whole core, logs, tests, and seismic to verify and quantify properties of vuggy, brecciated, and fractured carbonate intervals. 11. Used complex geocellular static and dynamic models to evaluate regional storage capacity using large parallel processing. 12. Carbonates are complex reservoirs and CO 2 -EOR needs to move to the next generation to increase effectiveness of CO 2 and efficiency and safety of the injection.

54 ENVIRONMENTAL SCIENCES↗

Modeling CO 2 Sequestration in Saline Aquifer and Depleted Oil Reservoir To Evaluate Regional CO 2 Sequestration Potential of Ozark Plateau Aquifer System, South-Central Kansas

1. Drilled, cored, and logged three wells to the basement and collecting more than 2,700 ft of conventional core; obtained 20 m i2 of multicomponent 3D seismic imaging and merged and reprocessed more than 125 mi 2 of existing 3D seismic data for use in modeling CO 2 - EOR oil recovery and CO 2 storage in five oil fields in southern Kansas. 2. Determined the technical feasibility of injecting and sequestering CO 2 in a set of four depleted oil reservoirs in the Cutter, Pleasant Prairie South, Eubank, and Shuck fields in southwest Kansas; of concurrently recovering oil from those fields; and of quantifying the volumes of CO2 sequestered and oil recovered during the process. 3. Formed a consortium of six oil operating companies, five of which own and operate the four fields. The consortium became part of the Southwest Kansas CO 2 -EOR Initiative for the purpose of sharing data, knowledge, and interest in understanding the potential for CO 2 -EOR in Kansas. 4. Built a regional well database covering 30,000 mi 2 and containing stratigraphic tops from ~90,000 wells; correlated 30 major stratigraphic horizons; digitized key wells, including wireline logs and sample logs; and analyzed more than 3,000 drill stem tests to establish that fluid levels in deep aquifers below the Permian evaporites are not connected to the surface and therefore pressures are not hydrostatic. Connectivity with the surface aquifers is lacking because shale aquitards and impermeable evaporite layers consist of both halite and anhydrite. 5. Developed extensive web applications and an interactive mapping system that do the following: a. Facilitate access to a wide array of data obtained in the study, including core descriptions and analyses, sample logs, digital (LAS) well logs, seismic data, gravity and magnetics maps, structural and stratigraphic maps, inferred fault traces, earthquakes, Class I and II disposal wells, and surface lineaments. b. Provide real-time analysis of the project dataset, including automated integration and viewing of well logs, core, core analyses, brine chemistry, and stratigraphy using the Java Profile app. A cross-section app allows for the display of log data for up to four wells at a time. 6. Integrated interpretations from the project’s interactive web-based mapping system to gain insights to aid in assessing the efficacy of geologic CO 2 storage in Kansas and insights toward understanding recent seismicity to aid in evaluating induced vs. naturally occurring earthquakes. 7. Developed a digital type-log system, including web-based software to modify and refine stratigraphic nomenclature to provide stakeholders a common means for communication about the subsurface. 8. Contracted use of a nuclear magnetic resonance (NMR) log and ran it slowly to capture response and characterize larger pores common for carbonate reservoirs. Used NMR to extend core analyses to apply permeability, relative permeability to CO 2 , and capillary pressure to the major rock types, each uniquely expressed as a reservoir quality index (RQI), present in the Mississippian and Arbuckle rocks. 9. Characterized and evaluated the possible role of microbes in dense brines. Used microbes to compliment H/O stable isotopes to fingerprint brine systems. Used perforation/swabbing to obtain samples from multiple hydrostratigraphic units and confirmed equivalent results using less expensive drill stem tests (DST). 10. Used an integrated approach from whole core, logs, tests, and seismic to verify and quantify properties of vuggy, brecciated, and fractured carbonate intervals. 11. Used complex geocellular static and dynamic models to evaluate regional storage capacity using large parallel processing. 12. Carbonates are complex reservoirs and CO 2 -EOR needs to move to the next generation to increase effectiveness of CO 2 and efficiency and safety of the injection.

54 ENVIRONMENTAL SCIENCES↗

Direct measurement of building transient and residual drift using an optical sensor system

Interstory drift (ID) is a key response parameter for buildings subjected to lateral loads and is used to define performance-based limit states, allowable deformations and damage states in a number of seismic design codes and standards. Here, an ability to rapidly and accurately measure both transient and residual ID during an earthquake would provide important observables for understanding the seismic demands and post-earthquake condition of a building. Accurate retrieval of ID from accelerometer-based instrumentation systems can be very challenging, if not impossible, as a result of instrumentation limitations and the post-processing associated with strong motion accelerometer data. This is particularly true for the case in which residual drifts occur during inelastic building response. In the study presented herein, a newly developed optical sensor system, designed specifically for directly measuring both transient and residual ID, was experimentally evaluated through shake table testing and computational simulations. The ability of the sensor to accurately measure ID is demonstrated and key operational considerations for sensor system deployment are examined through a model-based investigation.

42 ENGINEERING↗

An efficient method to propagate model uncertainty when inverting seismic data for time domain seismic moment tensors

SUMMARY We present a computationally efficient method to approximately propagate uncertainty when linearly inverting seismic data for point source, time variable moment tensor components. The method is based on the assumption that the data residual, given by the difference between the observed seismic data and the data predicated by a linear inversion, contains the effects of both data and model uncertainty. Our method uses a distribution of data residuals, added directly to the data, in a pseudo-Monte Carlo scheme. Using the assumption that the data residual is a stochastic process, we use the well-known Karhunen–Loève (KL) theorem to construct a distribution of data residuals, where the required basis functions are constructed using Fourier series. The Fourier series are scaled by a product of a random variable and the real-valued spectral amplitudes of the original data residual’s spectrum. Thus, the Fourier series and spectral amplitudes are eigenfunction-eigenvalue pairs used in the KL-based construction of data residual distribution. Using tests with synthetic data, we show that our method compares closely with a Finite Difference Monte Carlo (FDMC) method that we presented previously. More importantly, the method presented here is computationally several orders of magnitude faster than our previous FDMC method, and requires no a priori assumptions of model and/or data uncertainty.

Poppeliers, Christian (ORCID:0000000159526849)↗

Characterization with dense array data of seismic sources in the shallow part of the San Jacinto fault zone

SUMMARY We analyse dominant sources identified in a catalogue of more than 156 000 localizations performed using a 26-d data set recorded by a dense array set on the San Jacinto fault near Anza, in California. Events were localized using an array processing technique called Match Field Processing. As for all array processing techniques, the quality of the event position decrease when the events are outside of the array. We thus separate localizations in and outside the array using simple geometrical conditions. We compare the time distribution of the localization to additional data such as meteorological data, day of human activity as well as existing catalogues to determine the nature of the dominant events located using our method. We find that most of the events located outside of the array could be attributed to a surface structure excited by wind. On the other hand, part of the localizations under the array occur during regional earthquakes and could correspond to diffraction on the fault's heterogeneities. The rest of the localizations inside the array could be generated by the fault itself.

Geochemistry & Geophysics↗

Applicability of Sentinel-1 Terrain Observation by Progressive Scans multitemporal interferometry for monitoring slow ground motions in the San Francisco Bay Area

Data from the Sentinel-1 satellite have already proven useful for investigating seismic and volcanic events since its launch in April 2014. The requirement of ultrahigh coregistration accuracy and the current relatively short time of Sentinel-1 acquisitions make its application challenging for studying slow deformation processes, such as fault creep and land subsidence. Here in this paper we analyze a set of 14 SAR images over the San Francisco Bay Area spanning 1 year from early 2015 to 2016. We show that implementing an existing Enhanced Spectral Diversity algorithm or using precise orbits together with a reference digital elevation model both yield the required coregistration accuracy for making use of the phase measurements in time series analysis of ground deformation. Following a thorough validation test, we update our estimates of Hayward Fault creep rate and confirm uplift due to recharge of the Santa Clara Valley aquifer system during the final summer of 4 year drought.

58 GEOSCIENCES↗

A deep learning-enhanced framework for multiphysics joint inversion

Joint inversion has drawn considerable attention due to the availability of multiple geophysical data sets, ever-increasing computational resources, the development of advanced algorithms, and its ability to reduce inversion uncertainty. A key issue of joint inversion is to develop effective strategies to link different geophysical data in a unified mathematical framework, in which the information obtained from different models can complement each other. We have developed a deep learning-enhanced joint inversion framework to simultaneously reconstruct different physical models by fusing different types of geophysical data. Traditionally, structure similarity constraints are pursued by joint inversion algorithms using manually crafted formulations (e.g., cross gradient). The constraint is constructed by a deep neural network (DNN) during the learning process. The framework is designed to combine the DNN and a traditional independent inversion workflow and improve the joint inversion result iteratively. The network can be easily extended to incorporate multiphysics without structural changes. Numerical experiments on the joint inversion of 2D DC resistivity data and seismic traveltime are used to validate our method. In addition, this learning-based framework demonstrates excellent generalization abilities when tested on data sets using different geologic structures. It also can handle different sensing configurations and nonconforming discretization.

Geochemistry & Geophysics↗

STILGAR End-of-Project Report

The Subsurface Tunnel Imaging LeveraGed by Analysis of Rayleigh wave ellipticity (STILGAR) project demonstrated an integrated geophysical approach for detecting, locating, and characterizing underground structural changes using dense seismic arrays and advanced inversion techniques. Field campaigns were conducted at two operational mines—the Redmond salt mine (Utah) and Graymont Pleasant Gap limestone mine (Pennsylvania)—providing real-world testbeds for monitoring anthropogenic subsurface activity. At the Redmond salt mine, seismic interferometry combined with back-projection inversion successfully identified continuous, low-amplitude signals from mining operations. The approach differentiated stationary from migrating anthropogenic sources, captured daily operational cycles, and validated the potential of passive seismic monitoring for remote detection of underground activity. At the Graymont Pleasant Gap mine, two dense seismic deployments in the spring and fall of 2023 generated over 4 TB of high-resolution data. Key outcomes included the relocation of 199 underground and 8 surface explosions with accuracies within tens of meters and the development of a 3D P-wave velocity model using the triple-difference tomography algorithm (tomoTD) that resolved major structural features such as the mine entrance, low-velocity tunnels, and roof-collapse areas. Ambient noise cross-correlation and back-projection analyses revealed persistent sources linked to ongoing mining activity, whereas horizontal-to-vertical spectral ratio (HVSR) and ellipticity studies confirmed stable site responses across seasons and identified soil thickness trends consistent with regional erosional and depositional processes. Checkerboard and sensitivity tests further validated the robustness of the tomographic results. Overall, the findings emphasize that although significant progress has been made in subsurface imaging, further work is needed to enhance the detection and localization of underground structures. Accurate imaging requires higher frequencies, yet anthropogenic sources tend to dominate the seismic record at those frequencies, and high-frequency surface waves are affected by higher modes that complicate interpretation. The improved detection and localization of human-induced signals enabled detailed temporal and spatial mapping of daily mine operations, demonstrating the feasibility of continuous anthropogenic source monitoring. Sensitivity to signals from nontraditional sources, such as fan operations, highlights the broader applicability of this approach to other industrial environments in which continuous and impulsive signals are present. The field campaigns produced a substantial volume of high-quality seismic data, supporting the development and testing of new methods for seismic source characterization and subsurface imaging. Future deployments should include sensors capable of recording lower frequencies to probe deeper structures, increase bandwidth to enhance resolution and sensitivity to both shallow and deep targets, and collect additional large-scale datasets to refine imaging and source characterization techniques. Moreover, conducting 3D modeling studies of seismic wavefields at higher frequencies will provide a better understanding of wave scattering and cavity–wavefield interactions in complex underground environments. In conclusion, the STILGAR project demonstrated that integrated seismic monitoring can effectively characterize underground operations, capturing both natural and anthropogenic signals. The approaches developed provide a foundation for improved detection, localization, and imaging of subsurface structures and are directly transferable to broader industrial monitoring applications.

58 GEOSCIENCES↗

Center for Frontiers of Subsurface Energy Security (Final Report)

For the past two decades a slowly growing consensus has been that rising temperature is being caused by increases in carbon dioxide (CO 2 ) in the Earth’s atmosphere. What is less consensual is what to do about it. Solutions range from reducing emissions to various ways to capture CO 2 (from effluents of power plants, as plants, etc.) and store it away from the atmosphere. The problem is of such magnitude that it is likely that we will need many methods. Also some methods would cause draconian harm to global economies. This work is about the storage part of the reducing of CO 2 in the atmosphere. There have been many methods proposed for storing CO 2 : in vegetation, at the bottom of oceans, in land farms, and in subsurface formations. The latter is particularly prominent given the large estimated volume in porous formation in the Earth, and the 50-year old experience in injecting CO 2 for enhanced oil recovery. This report is about the science behind subsurface storage, called here geological carbon storage or GCS. Effective GCS requires that billions of metric tons of CO 2 be injected into storage reservoirs annually. GCS feasibility requires balancing the opposing impacts of large emplacement (injection) rates without compromising caprock seal integrity during injection. Once injected, the challenge is to assure the injected CO 2 remains permanently stored. Under most conditions, CO 2 is in a supercritical state, or scCO 2 . scCO 2 is immiscible with the subsurface brine and has a factor of ten smaller viscosity, leading to the development of two-phase fluid flow and scCO 2 bypassing. After the injection of large volumes of scCO 2 , the pore pressure, state of stress, chemical, thermal, and biological steady-state or equilibrium conditions in the subsurface are disrupted, which may lead to unpredictable behavior because the responses to these changes are often non-linear. Current challenges for efficient and reliable scCO 2 storage are: (1) sustaining large injection rates; (2) using available pore space efficiently, and (3) controlling undesired or unexpected behavior, such as scCO 2 leakage. Various aspects of these challenges are addressed by collaborative projects within Center for Frontiers in Subsurface Energy Security (CFSES). Challenge 1 – Sustaining large storage rates. To offset the annual CO 2 emissions of in the United States, ~7 billion metric tons of CO 2 must be injected into the storage reservoirs annually. Injected CO 2 forms a plume that starts spreading, displacing originally resident brine, though inefficiently so. Because of its low viscosity, and buoyant forces, a CO 2 plume rises to the reservoir rock-caprock interface. Injecting large volumes of scCO 2 at sufficiently high rates into subsurface storage formations leads to increases in pore pressures, potential expansion of the reservoir rock, and may result in fracturing of the reservoir rock and/or caprock. CFSES research uses existing and new experimental and modeling approaches to identify chemical-physical controls on the permeability and pore pressure dynamics, anticipating geomechanical (fracturing) events, predicting multi-phase flow patterns and trapping of CO 2 . Challenge 2 – Using pore space with unprecedented efficiency. Since large volumes of scCO 2 are to be stored, the volume available in storage reservoirs must be used efficiently. Storage efficiency is defined as a fraction of pore space occupied by carbon dioxide. Current estimates are that less than 5 percent of available pore volume is available for storage because scCO 2 , having a smaller density and viscosity than reservoir brine, tends to form viscous fingers, and bypass available pore space. This flow regime results in a CO 2 plume spreading over large areas in subsurface. The research goal in Challenge 2 is to advance the fundamental understanding of two-phase flow through porous media, including molecular-to meter- scale experimental and modeling approaches. This work leads to novel strategies to ensure that injected CO 2 occupies more than 5% of the pore space, ideally up to 50%. Challenge 3 – Controlling undesired or unexpected behavior. Following the injection of scCO 2 , the subsurface reservoir re-equilibrates to different physical and chemical states. This re-equilibration in space and time could lead to unexpected and/or undesired behavior: fracturing of reservoir rock and/or caprock, development of preferential flow paths, CO 2 leakage (to the surface or overlaying fresh water reservoirs), induced seismicity, and activation of faults. One of the main research goals of CFSES was to advance the fundamental understanding of how individual chemical and physical processes are coupled. We particularly focused on the chemical-mechanical coupling, and flow-mechanical coupling. We collected experimental and modeling data to determine the temporal and spatial scales of the coupled processes, which may lead to pore collapse caused by creep, activation of fractures and faults, chemically-induced fracture propagation, and self-focusing of scCO 2 flow.

54 ENVIRONMENTAL SCIENCES↗

Numerical modelling of impact seismic sources using the stress glut theory

SUMMARY Meteorite impacts have proved to be a significant source of seismic signal on the Moon, and have now been recorded on Mars by InSight seismometers. Understanding how impacts produce seismic signal is key to the interpretation of this unique data, and to improve their identification in continuous seismic records. Here, we use the seismic Representation Theorem, and particularly the stress glut theory, to model the seismic motion resulting from impact cratering. The source is described by equivalent forces, some resulting from the impactor momentum transfer, and others from the stress glut, which represents the mechanical effect of plasticity and non linear processes in the source region. We condense these equivalent forces into a point-source with a time-varying single force and nine-component moment tensor. This analytical representation bridges the gap between the complex dynamics of crater formation, and the linear point-source representation classically used in seismology. Using the multiphysics modelling software HOSS, we develop a method to compute the stress glut of an impact, and the associated point-source from hypervelocity impact simulations. For a vertical and an oblique impact at 1000 m s−1, we show that the moment tensor presents a significant deviatoric component. Hence, the source is not an ideal isotropic explosion contrary to previous assumptions, and draws closer to a double couple for the oblique impact. The contribution of the point force to the seismic signal appears negligible. We verify this model by comparing two signals: (1) HOSS is coupled to SPECFEM3D to propagate the near-source signal elastically to remote seismic stations; (2) the point-source model derived from the stress-glut theory is used to generate displacements at the same distance. The comparison shows that the point-source model is accurately simulating the low-frequency impact seismic waveform, and its seismic moment is in trend with Lunar and Martian impact data. High-frequencies discrepancies exist, which are partly related to finite-source effects, but might be further explained by the difference in mathematical framework between classical seismology and HOSS’ numerical modelling.

79 ASTRONOMY AND ASTROPHYSICS↗