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 127 records · Page 7

Real-time deep-learning inversion of seismic full waveform data for CO 2 saturation and uncertainty in geological carbon storage monitoring

Deep-learning inversion has recently drawn attention in geological carbon storage research due to its potential of imaging and monitoring carbon storage in real time, significantly improving efficiency and safety of carbon storage operations. We present a deep-learning full waveform inversion method that after the neural network has been trained can image CO 2 saturation and its uncertainty in real time. Our deep-learning inversion method is based on the U-Net architecture with the neural network trained on pairs of synthetic seismic data and CO 2 saturation models. Accordingly, our training establishes a mapping relationship between seismic data and CO 2 saturation models and once fully trained directly estimates CO 2 saturation as a function of subsurface location. We further quantify uncertainties of CO 2 saturation estimates using the Monte Carlo dropout method and a bootstrap aggregating method. For this proof-of-concept study, the CO 2 training models and data are derived from the Kimberlina 1.2 model, a hypothetical 3D geological carbon storage model that is constructed based on various geological and hydrological data from the Southern San Joaquin Basin, California. We perform deep-learning inversion experiments using noise-free and noisy training and test data sets and compare the results. Our modelling experiments show that (1) the deep-learning inversion can estimate 2D distributions of CO 2 fairly well even in the presence of Gaussian random noise and (2) both CO 2 saturation imaging and uncertainty quantification can be done in real time. Our results suggest that the deep-learning inversion method can serve as a robust real-time monitoring tool for geological carbon storage and/or other time-varying reservoir/aquifer properties that result from injection, extraction, and/or other subsurface transport phenomena.

58 GEOSCIENCES↗

Three-dimensional seismic characterization and imaging of the Soda Lake geothermal field

Accurate characterization of subsurface geophysical properties and detection of the fault system are essential for geothermal energy exploration and production. The Soda Lake geothermal field is in western Nevada with a complex fault system. Previous seismic characterization only produced a low-resolution, smooth velocity model along with a simple, conceptual fault model. Using optimized correlation-based full-waveform inversion, wavefield-separation-based reverse-time migration, and automatic fault detection techniques, we present 3D seismic characterization for the Soda Lake geothermal field using 3D surface seismic data acquired with Vibroseis sources. Here, we obtain 3D high-resolution velocity, density, and acoustic impedance models, 3D seismic images with different grid spacings, and a high-resolution fault system. Consistency check between the constructed faults and currently active injection and production geothermal wells verifies that our seismic inversion and imaging results and detected faults are reliable. These results can provide valuable information for optimizing well placement and geothermal energy production at the Soda Lake geothermal field.

15 GEOTHERMAL ENERGY↗

Physics-based Induced Earthquake Forecasting: Process Understanding, and Hazards Mitigation

Disposal of saltwater co-produced with oil and gas is linked to elevated seismicity in the Central and Midwest US. There is a concern that these events may lead to widespread damage and an overall increase in seismicity. Thus an improved understanding of the spatially and temporally variable deformation and stress field associated with fluid injection operation is critically necessary for evaluating time-varying seismic hazards. Despite the improvements in seismic monitoring capacity and the resulting decrease in the magnitude detection threshold, estimates of induced earthquake probability remain elusive due to insufficient models incapable of accounting for the complex physics governing the process of induced seismicity. The proposed research effort will comprehensively analyze, integrate, and interpret geodetic, injection, and seismic data in the vicinity of the injection sites in Oklahoma to resolve the 4-dimensional distribution of pore pressure and stress in the shallow crust. This project, in particular, is focused on exploring the statistical relation between injection operation and increased earthquake hazard. The amplitude of and the extent to which pore pressure changes are determined by some factors, in particular, the hydrogeological properties of the rocks, such as diffusivity. Thus the available deformation data will be used to constrain hydrogeological properties of the medium, to accurately resolve the evolution of crustal stresses due to fluid injection. Having the time-varying models of stress changes, a statistical framework will be implemented to estimate the time-dependent probability of large earthquakes on the nearby fault systems. These data and models help to improve seismic hazard estimates and aid in constructing operational-induced earthquake forecast models. This information can also be integrated into the updated U.S. National Seismic Hazard Map, which local communities and authorities use in their earthquake risk estimates and mitigation efforts.

58 GEOSCIENCES↗

A Methodological Overview of Seismic Analysis for Nuclear Event Detection

Underground explosions generate potentially detectable signatures, including energy waves that travel through the Earth’s subsurface (i.e., seismic waves), low-frequency sound waves (i.e., infrasound and hydroacoustic waves), and radioactive gases and/or particles that might leak from the test cavity (if the event was nuclear). There can also be intelligence indicators of a test, such as observations of modified patterns of life and activity at a suspected test site. If all of these detectable signatures and intelligence indicators are present and self-consistent, then analysts have high confidence in classifying a signature generating event as an explosion. However, because only partial information about an event is likely to be available, determining whether an event was natural (e.g., an earthquake or landslide) or manmade (e.g., a chemical or nuclear explosion) is much more challenging. This primer describes how one category of event signatures—seismic signatures—can augment event analyses. While universities and government organizations have generated detailed technical descriptions of seismic analytic techniques, we seek to translate seismic event analysis for a broad, non-technical audience. When the geologic conditions near an event are well-characterized, seismic data can be used to calculate critical information, such as event location and depth, with relatively high accuracy. Moreover, specific features within seismic datasets can help determine whether an event was an explosion. However, a key challenge in seismic analysis is that geologic site conditions are often poorly characterized, complicating the ability to discern the true nature of the event. To overcome this challenge, geologists answer a series of questions (discussed in section 1) to guide seismic event analysis and determine the most probable nature of an event. As more information is gathered during each analytic step, confidence grows regarding the nature of the event. Section 2 addresses uncertainties in seismic analysis and the vital nature of high-fidelity geologic data for accurate seismic event analysis.

58 GEOSCIENCES↗

A Multimodal Event Catalog and Waveform Data Set That Supports Explosion Monitoring from Nevada, U.S.A.

Multimodal, curated data sets and nuisance event catalogs remain rare in the explosion monitoring community relative to curated seismic data sets. The source of this relative absence is the difficultly in deploying multimodal receivers that sense the seismic, acoustic, and other modalities from multiphysics sources. We provide such a data set in this study that delivers seismic, infrasound, and electromagnetic (magnetometer) sensor records collected over a two–week period, within 255 km of a 10 ton buried chemical explosion called DAG–4 that was located at 37.1146°, –116.0693° on 22 June 2019 21:06:19.88 UTC. This catalog includes 485 seismic, seismoacoustic, and infrasound–only events that an expert analyst manually built by reviewing waveforms from 29 seismic and infrasound sensors. Our data release includes waveforms from these 29 seismic, infrasound, and seismoacoustic stations and two magnetometer stations and their station metadata. We deliver these waveforms in NNSA KB Core CSS.w format (i4) with a corresponding wfdisc table that provides the header information. Here, we expect that this data set will provide a valuable, benchmark resource to develop signal processing algorithms and explosion monitoring methods against manual, human observations.

58 GEOSCIENCES↗

Induced microearthquakes predict permeability creation in the brittle crust

Predicting the evolution of permeability accurately during stimulation at the reservoir scale and at the resolution of individual fractures is essential to characterize the fluid transport and the reactive/heat-transfer characteristics of reservoirs where stress exerts significant control. Here, we develop a hybrid machine learning (ML) model to visualize in situ permeability evolution for an intermediate-scale (~10 m) hydraulic stimulation experiment. This model includes an ML model that was trained using the well history of flow rate and wellhead pressure and MEQ data from the first three stimulation episodes to predict average permeability from the statistical features of the MEQs alone for later episodes. Moreover, a physics-inspired model is integrated to estimate in situ fracture permeability spatially. This method relates fracture permeability to fracture dilation and scales dilation to the equivalent MEQ magnitude, according to laboratory observations. The seismic data are then applied to define incremental changes in permeability in both space and time. Our results confirm the excellent agreement between the ground truth and model-predicted permeability evolution. The resulting permeability map defines and quantifies flow paths in the reservoir with the averaged permeability comparing favorably with the ground truth of permeability.

induced seismicity↗

TCCSP Site Characterization and Feasibility Report

The Tulare County Carbon Storage Project (TCCSP) focused on establishing the geological and commercial viability of a 50 million metric ton (MMT) carbon storage complex located in Tulare County, California. The geologic scope of work included collecting new and existing subsurface data, such as two-dimensional (2D) seismic data, advanced geophysical logs, core and fluid samples, and specialized core-analysis datasets, to support the integration and interpretation of findings from earlier tasks. These efforts contributed to the comprehensive feasibility study for commercial carbon capture and storage (CCS) development at the TCCSP site and in the surrounding region. The TCCSP project team evaluated the feasibility of the Lower Monterey Group (Santa Margarita Sand Member or its equivalent), Temblor Formation sand members (Olcese, Jewett, Vedder), and Lodo-Martinez Zone (Domengine Sandstone, Walker Formation, Lodo Formation, and the Martinez Sand) to serve as commercial CCS reservoirs along with potential for the Upper Monterey Group (Reef Ridge, Antelope Members) and the Kreyenhagen Formation to vertically seal underlying reservoirs from shallow above-zones.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Dense Seismic Array Study of a Legacy Underground Nuclear Test at the Nevada National Security Site

The complex postdetonation geologic structures that form after an underground nuclear explosion are difficult to constrain because increased heterogeneity around the damage zone affects seismic waves that propagate through the explosion site. Generally, a vertical rubble-filled structure known as a chimney is formed after an underground nuclear explosion that is composed of debris that falls into the subsurface cavity generated by the explosion. Compared with chimneys that collapse fully, leaving a surface crater, partially collapsed chimneys can have remnant subsurface cavities left in place above collapsed rubble. The 1964 nuclear test HADDOCK, conducted at the Nevada test site (now the Nevada National Security Site), formed a partially collapsed chimney with no surface crater. Understanding the subsurface structure of these features has significant national security applications, such as aiding the study of suspected underground nuclear explosions under a treaty verification. In this study, we investigated the subsurface architecture of the HADDOCK legacy nuclear test using hybrid 2D–3D active source seismic reflection and refraction data. The seismic data were acquired using 275 survey shots from the Seismic Hammer (a 13,000 kg weight drop) and 65 survey shots from a smaller accelerated weight drop, both recorded by ~1000 three-component 5 Hz geophones. First-arrival, P-wave tomographic modeling shows a low-velocity anomaly at ~200 m depth, likely an air-filled cavity caused by partial collapse of the rock column into the temporary postdetonation cavity. A high-velocity anomaly between 20 and 60 m depth represents spall-related compaction of the shallow alluvium. Hints of low velocities are also present near the burial depth (~364 m). The reflection seismic data show a prominent subhorizontal reflector at ~300 m depth, a short-curved reflector at ~200 m, and a high-amplitude reflector at ~50 m depth. Comparisons of the reflection sections to synthetic data and borehole stratigraphy suggest that these features correspond to the alluvium–tuff contact, the partial collapse cavity, and the spalled layer, respectively.

58 GEOSCIENCES↗

Improving Subsurface Stress Characterization for Carbon Dioxide Storage Projects by Incorporating Machine Learning Techniques

The overall objective of this project is to develop a framework for reliable characterization and prediction of the state of stress in the overburden and underburden (including the basement) in CO 2 storage reservoirs using machine learning and integrated geomechanics and geophysical methods. Specifically, we propose to develop workflow encompassing of technologies and/or methods to predict stress and pressure changes due to CO 2 injection in an active tertiary recovery site and their impacts on subtle fault activation, fractures and occurrence of microseismic events and compare responses to field observations. In this project, we anticipate using dataset from the Farnsworth field Unit (FWU) which is operated by Purdure Petroleum. A novel elastic-waveform VSP inversion technique will be used to estimate high-resolution spatial and temporal changes of elastic moduli in CO 2 storage reservoirs, which will be combined with velocity-stress relationship derived from laboratory tests to obtain subsurface pressure and stress. Clustered microseismic data will be jointly inverted for improved focal mechanisms. Least-squares reverse-time migration of microseismic waveform data will be performed to directly image fracture/fault zones. Additionally, a deep neural network machine learning technique with convolutional and recurrent layers will be used for learning the spectro-temporal structures in microseismic waveforms. The results of this geotechnical data analysis will be integrated to develop a high-resolution 3D mechanical earth model extending from the overburden sealing formations to the underburden including the basement. Mechanical properties will be derived through integration of mechanical logs, tests, available results from chemo-mechanical laboratory tests, and elastic inversion of seismic data using a combination of Bayesian and stochastic methods as well as machine learning technique. Failure features (faults/fractures) will be represented and/or modeled based on seismic and core data analysis. A transient hydrodynamic-geomechanical model will be developed through coupling with the calibrated FWU reservoir simulation model. The full physics coupled model will be used to train a reduced order proxy model using machine learning algorithm for estimating stress which will then be used with appropriate constitutive relationships and forward seismological models to simulate pressure changes and induced microseismicity. An advanced optimization framework will be developed to perform a history match to minimize error between field observations and simulated. The history matched proxy model will be verified against the full-physics equivalent. The field observations that will be used in the coupled model calibration process include pressure/stress inverted from VSP, moment magnitude from microseismic analysis, real time downhole pressure measurements, production and injection data. Parameter sensitivity and uncertainty analysis will be performed to characterize the impact of model parameter uncertainty on stress estimates. The proposed project will have significant impact on future field implementation of the proposed technology. Because the project field site is an ongoing CO 2 EOR development, the value of the new technology will be demonstrated in an operational context and evaluated as a viable risk mitigation strategy. Cost/benefit will be evaluated together with the various commercial incentives for CO 2 sequestration available to oil and gas operators. The extensive available dataset and ongoing data acquisition under the SWP Phase III work plan provides flexibility for investigation of multiple approaches and reduces technical risk.

58 GEOSCIENCES↗

Real-time reconstruction of ground motion during small magnitude earthquakes: A pilot study

This study presents a pilot investigation into a novel method for reconstructing real-time ground motion during small magnitude earthquakes (M < 4.5), removing the need for computationally expensive source characterization and simulation processes to assess ground shaking. Small magnitude earthquakes, which occur frequently and can be modeled as point sources, provide ideal conditions for evaluating real-time reconstruction methods. Utilizing sparse observation data, the method applies the Gappy Auto-Encoder (Gappy AE) algorithm for efficient field data reconstruction. This is the first study to apply the Gappy AE algorithm to earthquake ground motion reconstruction. Numerical experiments conducted with SW4 simulations demonstrate the method’s accuracy and speed across varying seismic scenarios. The reconstruction performance is further validated using real seismic data from the Berkeley area in California, USA, demonstrating the potential for practical application of real-time earthquake data reconstruction using Gappy AE. As a pilot investigation, it lays the groundwork for future applications to larger and more complex seismic events.

58 GEOSCIENCES↗

Seismic Survey 2016 Data at Crescent Valley, Nevada

In September 2016, 989 vertical-component seismic instruments were deployed for 75 hours at the Crescent Valley greenfield geothermal play area in Nevada. Data were recorded 12/9/16 - 12/15/16. Data are stored in individual files in one-minute increments in SEGD formats. See the metadata in Crescent Valley Station Locations.sp1, Crescent Valley Observer Logs.xls, and the attached Metadata link for details about the seismic station locations, seismic data logger specifications, instrumentation specifications, descriptions of data, a fracture finding summary, and the final report for the 2016 seismic survey done in Crescent Valley, Nevada. Note that no discrete events were identified, and a 1750m x 2250m x 300m PSET volume of acoustic energy was delivered.

15 GEOTHERMAL ENERGY↗

Modulation of seismic noise near the San Jacinto fault in southern California: origin and observations of the cyclical time dependence and associated crustal properties

SUMMARY We examine the cyclic amplitude variation of seismic noise recorded by continuous three-component broad-band seismic data with durations spanning 91–713 d (2008–2011) from three different networks: Anza seismic network, IDA network and the transportable seismic array. These stations surround the San Jacinto Fault Zone (SJFZ) in southern California. We find the seismic noise amplitudes exhibit a cyclical variation between 0.3 and 7.2 Hz. The high-frequency (≥0.9 Hz) noise variations can be linked to human activity and are not a concern. Our primary interest is signals in the low frequencies (0.3–0.9 Hz), where the seismic noise is modulated by semi-diurnal tidal mode M2. These long-period (low-frequency) variations of seismic noise can be attributed to a temporal change of the ocean waves breaking at the shoreline, driven by ocean tidal loading. We focus on the M2 variation of seismic noise at f = 0.6 Hz, travelling distances of ∼92 km through the crust from offshore California to the inland Anza, California, region. Relative to the shoreline station, data from the inland stations show a phase lag of ∼ –12°, which we attribute to the cyclic change in M2 that can alter crustal seismic attenuation. We also find that for mode M2 at 0.6 Hz, the amplitude variations of the seismic quality factor (Q) depend on azimuth and varies from 0.22 per cent (southeast to northwest) to 1.28 per cent (northeast to southwest) with Q = 25 for Rayleigh waves. We propose the direction dependence of the Q variation at 0.6 Hz reflects the preferred orientation of subfaults parallel to the main faulting defined by the primarily N45°W strike of the SJFZ.

Geochemistry & Geophysics↗

Development of Time Lapse VSP Integration Workflow: A Case Study at Farnsworth CO2-EOR Project

Abstract This study aims to develop a 4D Vertical Seismic Profile (VSP) integration workflow to improve the prediction of subsurface stress changes. The selected study site is a 5-spot pattern within the ongoing CO2-EOR operations at the Farnsworth Field Unit FWU in Ochiltree County, Texas. The specific pattern has undergone extensive geological and geomechanical characterization through the acquisition of 3D seismic data, geophysical well logs, and core. This workflow constrains a numerical hydromechanical model by applying a penalty function formed between "modeled" versus "observed" time-lapse compressional and shear seismic velocity changes. Analyses of geophysical logs and ultra-sonic measurements on core exhibit measurable sensitivities to changes in both fluid saturation and mean effective stress. These data are used to develop a site-specific rock physics model and stress-velocity relationship, which inform the numerical models used to generate the "modeled" portion of the penalty function. The "observed" portion of the penalty function is provided by a novel elastic full-waveform inversion of the available 3D baseline and three monitor surveys to produce high-quality estimates of time-lapse compressional and shear seismic velocity changes. The modeling workflow accounts sequentially for fluid substitution and stress impacts. Hydrodynamic and geomechanical properties of the 3D coupled numerical model are estimated through geostatistical integration of well log and core data with 3D seismic inversion products. Changes in seismic velocities due to fluid substitution are computed using the Biot-Gassmann workflow and site-specific rock physics. Stress impacts on time-lapse seismic velocity changes are modeled from the effective stress output of the hydromechanical model and are initially based on the velocity versus effective stress relationship extracted from core mechanical testing. Based on the principle of superposition of seismic wavefields, seismic velocity changes attributed to fluid substitution and that due to changes in mean effective stress are treated as linearly additive. The modeled results are upscaled using Backus averaging to reconcile scale discrepancies between the modeled and measured datasets to formulate the penalty function. This manuscript presents the forward modeling process and concludes that for the base case, the seismic velocity changes due to mean effective stress dominates over the seismic velocity changes attributed to fluid substitution because of the extensive range of the pressure perturbations. Successful minimization of this penalty function calibrates the coupled hydrodynamic geomechanical numerical model and affirms the suitability of acoustic time-lapse measurements such as 4D-VSP for geomechanical calibration.

02 PETROLEUM↗

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

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

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

Using Multiple Geophysical Methods to Refine a Stratigraphic Conceptual Site Model at a Nuclear Waste Site

At contaminated sites, there is a critical reliance on conceptual site models of which a key component is the geologic framework model (GFM), which is used to build predictive fate and transport modeling in support of remediation goals. Here, we detail an ongoing study at a nuclear waste site (Hanford Site, Washington, USA) where surface geophysical methods were used to complement existing borehole information and to help site new boreholes to update and refine the GFM. Ground-based electrical resistivity tomography (ERT), time-domain electromagnetics (TEM), and seismic data was acquired on the Hanford Site’s Central Plateau, where the presence of a hydraulically transmissive subsurface feature, or paleochannel, was suspected but uncertain based on contaminant concentrations at sparsely located boreholes. To assess and compare ERT, TEM, and seismic geophysical datasets, co-located data was acquired in an area where the existing GFM was more certain. Surrounding the suspected paleochannel, ERT and seismic refraction tomography showed similar subsurface structures consistent with a channelized feature. While TEM interpretation was limited in the suspected paleochannel area due to a thick resistive top layer and high electromagnetic noise, TEM showed greater utility in a different area where these factors were less prevalent. Based on the results of quasi-3D ERT inversions, we propose updates to the GFM, including identifying a paleochannel. Further, we evaluate the utility of geophysical methods and examine lessons learned that will be used as a basis for continued characterization efforts to support site management decisions and implementation efforts.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Central Asia Seismic Hazard Assessment (CASHA): A Probabilistic Seismic Hazard Assessment for Kazakhstan, Kyrgyzstan and Tajikistan

Probabilistic seismic hazard assessments (PSHA) underpin the calculation of earthquake loads in most building codes around the world. In Central Asia, the building codes are slowly being updated to incorporate some of the contemporary concepts of seismic hazard representation. There is also a regional desire to coordinate hazard assessments and building code modernization. However, some challenges remain. Expertise in the region related to seismic hazard assessments is still largely compartmentalised, requiring a significant amount of training and capacity building in seismic hazard assessment related topics. In addition, there are vast amounts of seismic data (bulletin and waveforms), both from analogue and digital eras, that the region’s countries stored but until recently did not use or share among themselves or with the broader seismological community around the world. Finally, after the collapse of the Soviet Union in the 1990s, many of the countries’ seismic networks suffered a major setback with the lack of attention and budget to update existing equipment and installation of new instruments. In order to address these issues, the United States Department of Energy through Lawrence Livermore National Laboratory (LLNL) initiated a project in 2016 to engage and train local scientists in Central Asia to install new equipment, to enhance the quality of seismic monitoring and reporting, to improve and harmonise the regional earthquake catalogue, and to conduct national probabilistic seismic hazard assessments using the new and improved datasets. To achieve the seismic hazard assessment related goals, a series of workshops were held in Almaty, Kazakhstan; Bishkek, Kyrgyzstan; and Dushanbe, Tajikistan from 2016 until 2020. During the time that the COVID-19 pandemic restricted travel, workshops continued online (22 online workshops were hosted in two years). Finally, in May 2022, an in-person workshop in Istanbul, Turkey brought together all project participants along with civil engineers engaged with building code activities in their respective countries, providing a platform to discuss the implementation of the hazard models into updates of building codes in each country, as well as to discuss model parameters, sensitivity analyses and model results in terms of hazard maps, uniform hazard spectra and hazard deaggregations. The workshops were a combination of lectures and hands-on exercises, and included international participation as well as local scientists and engineers. The workshops served several purposes, including training, coordination of data collection, interactions between local earth scientists and engineers, and brainstorming and knowledge exchange among local and international experts. This report outlines the new earthquake catalogue compilation effort and the PSHA project undertaken in Kyrgyzstan, Tajikistan, and Kazakhstan as part of this initiative. The southern part of this region is tectonically active with moderate to high levels of both shallow crustal seismic activity and occurrence of deeper earthquakes under the Hindu Kush and Pamir mountain ranges. Deeper earthquakes also occur near southwestern Kazakhstan, under the eastern Greater Caucasus and Caspian Sea. Large portions of central and northern Kazakhstan, on the other hand, are in stable continental regions with low levels of seismic activity. This study systematically compiles and improves all available data on local seismicity, active faults, and ground motion attenuation characteristics of the region; and builds a framework to enable a contemporary PSHA to be carried out with the engagement of local scientists. While the project was regional, the seismic hazard assessments are primarily driven by the countries’ own national preferences and understanding of data collection, interpretation, and validation of results.

58 GEOSCIENCES↗

Hydrologic Impacts of a Strike-Slip Fault Zone: Insights from Joint 3D Body-Wave Tomography of Rock Valley

The Rock Valley fault zone (RVFZ), an intraplate strike-slip fault zone in the southern Nevada National Security Site (NNSS), hosted a series of very shallow (<3 km) earthquakes in 1993. The RVFZ may also have hydrological significance within the NNSS, potentially playing a role in regional groundwater flow, but there is a lack of local hydrological data. In the Spring of 2021, we collected active-source accelerated weight drop seismic data over part of the RVFZ to better characterize the shallow subsurface. Here, we manually picked ~17,000 P-wave travel times and over 14,000 S-wave travel times, which were inverted for P-wave velocity (V P ), S-wave velocity (V S ), and V P /V S ratio in a 3D joint tomographic inversion scheme. Seismic velocities are imaged as deep as ~700 m in areas and generally align with geologic and structural expectations. V P and V S are relatively reduced near mapped and inferred faults, with the most prominent lower V P and V S zone around the densest collection of faults. We image V P /V S ratios ranging from ~1.5 to ~2.4, the extremes of which occur at a depth of ~100 m and are juxtaposed across a fault. One possible interpretation of the imaged seismic velocities is enhanced fault damage near the densest collection of faults with relatively higher porosity and/or crack density at ~100 m depth, with patches of semiperched groundwater present in the sedimentary rock in higher V P /V S areas and drier rock in lower V P /V S areas. A relatively higher V P /V S area beneath the densest faults persists at depth, which suggests percolation of groundwater via the fault damage zone to the regionally connected lower carbonate aquifer. Potentially, the presence and movement of groundwater may have played a role in the 1993 earthquake aftershocks.

58 GEOSCIENCES↗