SEARCH · Engineering Papers

Results for “seismic”

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.

SourceResearch connectionFrom yearThrough year

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

At least 289 records · Page 16

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Urban Seismic Site Characterization by Fiber-Optic Seismology

ccurate ground motion prediction requires detailed site effect assessment, but in urban areas where such assessments are most important, geotechnical surveys are difficult to perform, limiting their availability. Distributed acoustic sensing (DAS) offers an appealing alternative by repurposing existing fiber-optic cables, normally employed for telecommunication, as an array of seismic sensors. We present a proof-of-concept demonstration by using DAS to produce high-resolution maps of the shallow subsurface with the Stanford DAS array, California. We describe new methods and their assumptions to assess H/V spectral ratio—a technique widely used to estimate the natural frequency of the soil—and to extract Rayleigh wave dispersion curves from ambient seismic field. These measurements are jointly inverted to provide models of shallow seismic velocities and sediment thicknesses above bedrock in central campus. The good agreement with an independent survey validates the methodology and demonstrates the power of DAS for microzonation.

Geochemistry & Geophysics↗

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 OSTI2021

Mitigating Injection‐Induced Seismicity Along Basement Faults by Extraction: Application to 2016–2018 Pohang Earthquakes

Abstract The risk of induced earthquake is considered as the major obstacle for sustainable subsurface energy activities. As a mitigation approach, mechanistic risk assessment can be essential to minimize the potential of seismic hazard. Coupled hydro‐mechanical modeling approach is integrated with geomechanical analysis to better assess the occurrence of seismicity induced by injection and/or extraction, which will provide a basis to update operational conditions (e.g., number and location of mitigating wells, rate and duration of operations) for mitigating earthquake nucleation. The Coulomb stability analysis and seismicity rate estimates from this generic study indicate that simultaneous extraction can relieve injection‐induced pore‐pressure buildup, but also transfer and accumulate elastic energy along the faults. Field‐scale simulation with variation in well design for the Pohang EGS site, South Korea, show that multiple wells aligned across a preexisting fault form gradients of pressure and stress fields normal to the fault, which can enhance the direct impact of diffusion process and poroelastic stressing on the fault stability perturbation. Our mechanistic study suggests multiphysics‐oriented mitigation strategies that consider (1) possible mechanical deformation, caused by hydro‐mechanical coupling effects, and (2) proper configuration of well operations based on site‐specific geological characteristics.

Chang, K. W.↗

DOE OSTI2021

Spectral Element Modeling of Acoustic to Seismic Coupling Over Topography

Acoustic waves in the atmosphere are commonly recorded on seismometers as they couple into the ground. These signals, here called ground coupled airwaves, are not commonly considered in numerical modeling of infrasound propagation, which often assumes a rigid unmeshed boundary. Starting from an analytically-tractable spherical wave model, we analyze how the coupling of an acoustic wave into a planar elastic halfspace changes over a wide range of scenarios. We use energy admittance to quantify acoustic to seismic coupling over both a planar elastic halfspace and meshed topography. Our spectral element and analytic calculations have different maxima as a function of incidence angle, with very high admittance values for near-vertical incidence (maximum ≈78%). Energy admittance calculations at shallow incidence angles are much smaller (less than 1%). In simulations over the complex topography of Sakurajima Volcano, we attribute the variable spatial pattern of energy admittance to changes in earth parameters between each model. The observed pressure difference over the simulated 15 km region appears to be < 2%. While this value is relatively small, the cumulative addition over 100s of km and multiple acoustic bounce points may be significant. Acoustic to seismic coupling along the propagation path may bias long distance yield estimates, particularly when infrasound propagates over regions with steep topography or particularly slow seismic velocities, such as alluvial planes.

58 GEOSCIENCES↗

DOE OSTI2021

Seasonal Seismicity in the Lake Biwa Region of Central Japan Moderately Modulated by Lake Water Storage Changes

Seasonal seismicity in Japan has been observed and attributed to variations induced by hydrological loading, as well as atmospheric and pore-fluid pressure. We examine the seasonal variations of crustal (<15 km depth) and low magnitude (M < 4.5) earthquakes in central Japan near the Biwako-Seigan Fault Zone (BSFZ), and analyze their correlation with the annual hydrological cycle of nearby Lake Biwa, the largest freshwater lake in Japan. The 2002–2018 lake water storage is estimated using water-level gauge and satellite altimetry data. Regional surface mass loading from soil moisture and snow predicted by the Global Land Data Assimilation System is used to augment the lake mass variation. We compute the stress and stressing-rate changes induced by the surface mass loading from lake storage, snow, and soil moisture at seismogenic depth (10 km) along the BSFZ using a finite element model incorporating a 3D Earth structure. The results show that seasonal surface mass loading generates 1–5 kPa peak-to-peak Coulomb stress change and 2.5–15 kPa/year stressing-rate variations on the fault segments along the BSFZ. Two of the four fault segments present significant correlation between local seismicity rate and loading-derived stress and/or stressing rate changes. Among the hydrological load components, lake water, snow, and soil moisture exhibit different annual phases, and the lake water storage dominates the total annual stress changes. Our analysis indicates that seasonal surface hydrological loads, especially lake water variations, modulate the seismicity rate along the BSFZ in central Japan.

58 GEOSCIENCES↗

DOE OSTI2023

The Variable Continuous Bimaterial Interface in the San Jacinto Fault Zone Revealed by Dense Seismic Array Analysis of Fault Zone Head Waves

Key factors controlling earthquake ruptures include fault geometry, continuity, and seismic velocity structure around the fault. We present a novel tool that better informs deep bimaterial fault geometry embedded in distributed damage and seismicity, associated velocity contrasts across the fault, and their correlations with surface complexities. The method employs fault zone head and direct body waves and is applied to recordings from five spatiotemporally different seismic arrays along the complex San Jacinto fault zone (SJFZ) in southern California. We detect and distinguish these signals based on instantaneous phase coherence and relative energy in a cascading manner from one scale array to another. The analysis reveals a >70-km long continuous bimaterial interface within the SJFZ with several deep northeast dipping fault segments. The northern SJFZ, for instance, locates ~7 km northeast of its surface expression at 18-km depth. P-wave velocity contrasts range from near 0% to >15%, consistent with other bimaterial faults, and differ by a few % depending on fault-array azimuth, implying directional-dependent velocity contrasts. S-wave head waves and velocity contrasts are also imaged for the first time at the southern SJFZ, averaging to 2.9% in agreement with tomography results. The imaged geometry and continuity suggest the SJFZ initiated along remnant tectonic structures and translates to a rupture potential of M > 7.2, i.e., the sizes of its largest paleo-earthquakes. The P and S contrasts, and their ratios, have important implications for earthquake rupture speed, mode, directivity, and frictional heating along the SJFZ and other major faults globally.

58 GEOSCIENCES↗

DOE OSTI2023

Examining Infrasound Propagation at High Spatial Resolution Using a Nodal Seismic Array

Infrasound—acoustic waves in the atmosphere below 20 Hz—is a useful monitoring tool. Topography and atmospheric structure strongly control infrasound propagation, and at common source–receiver distances neither of these effects can be ignored when quantitative source constraints are sought. Detailed spatial measurements of the infrasound wavefield would inform propagation models and improve source estimates. However, the “large-N” deployment strategy now well-known in seismology has not yet been realized for infrasound studies. Here, we use the 900-node seismic array from the 2014 Imaging Magma Under St. Helens (iMUSH) experiment as a proxy for a large-N infrasound network, by leveraging acoustic–seismic coupled arrivals. The active-source component of iMUSH consisted of 23 shallowly buried explosions around Mount Saint Helens volcano; these explosions produced epicentral infrasound recorded on the nodes. We find that the bulk presence of ground-coupled infrasound on the nodes is controlled by wind noise and source–receiver distance, with observed arrivals for eight explosions. Explosions with the most extensive coupling produce complex spatial waveform patterns across the array. These patterns are related to both topographic and atmospheric propagation effects, as well as spatially variable site (coupling) effects. We compare our observations to simple topographic diffraction and high-resolution wind advection models, and full-wave numerical simulations. We find strong spatial correlations between (a) coupled arrival strength and modeled topographic obstruction and (b) coupled arrival time and along-path winds. Our seismoacoustic analyses and results are applicable to other existing and future nodal seismic data sets and can expand the utility of such deployments.

58 GEOSCIENCES↗

DOE OSTI2025

Analysis of Seismic Wave Amplification in Sedimentary Basins Using 3D Wavefield Simulations: Nenana Basin, Central Alaska

We investigate seismic wave amplification in Nenana basin, central Alaska, using 3D seismic wavefield simulations. We quantify frequency-dependent amplification by comparing synthetic seismograms among four different seismic velocity models: (a) a tomographic model, (b) a tomographic model with the uppermost 6.5 km replaced by a 1D basin profile, (c) a tomographic model with an embedded ellipsoidal basin, and (d) a tomographic model with an embedded realistic basin. For each model we perform wavefield simulations accurate up to 2 Hz for a set of 10 different earthquakes, which provide a range of frequencies and incident angles for waves interacting with the basin. Average amplification ratios are four on the horizontal components and seven on the vertical component. The dominant influence of the amplification is the 3D geometry of the basin, not the slow-velocity profile near the surface. Our synthetic amplification ratios are in general agreement with amplification estimates obtained from 14 stations that recorded the same set of earthquakes. Our approach offers a general strategy for documenting frequency-dependent basin amplification for a region with realistic basin structures and local earthquakes.

Geosciences↗

DOE OSTI2020

An exponential build-up in seismic energy suggests a months-long nucleation of slow slip in Cascadia

Slow slip events result from the spontaneous weakening of the subduction megathrust and bear strong resemblance to earthquakes, only slower. This resemblance allows us to study fundamental aspects of nucleation that remain elusive for classic, fast earthquakes. We rely on machine learning algorithms to infer slow slip timing from statistics of seismic waveforms. We find that patterns in seismic power follow the 14-month slow slip cycle in Cascadia, arguing in favor of the predictability of slow slip rupture. Here, we show that seismic power exponentially increases as the slowly slipping portion of the subduction zone approaches failure, a behavior that shares a striking similarity with the increase in acoustic power observed prior to laboratory slow slip events. Our results suggest that the nucleation phase of Cascadia slow slip events may last from several weeks up to several months.

58 GEOSCIENCES↗

DOE OSTI2024

Crustal permeability generated through microearthquakes is constrained by seismic moment

We link changes in crustal permeability to informative features of microearthquakes (MEQs) using two field hydraulic stimulation experiments where both MEQs and permeability evolution are recorded simultaneously. The Bidirectional Long Short-Term Memory (Bi-LSTM) model effectively predicts permeability evolution and ultimate permeability increase. Our findings confirm the form of key features linking the MEQs to permeability, offering mechanistically consistent interpretations of this association. Transfer learning correctly predicts permeability evolution of one experiment from a model trained on an alternate dataset and locale, which further reinforces the innate interdependency of permeability-to-seismicity. Models representing permeability evolution on reactivated fractures in both shear and tension suggest scaling relationships in which changes in permeability (Δk) are linearly related to the seismic moment (M) of individual MEQs as Δk ∝ M. This scaling relation rationalizes our observation of the permeability-to-seismicity linkage, contributes to its predictive robustness and accentuates its potential in characterizing crustal permeability evolution using MEQs.

15 GEOTHERMAL ENERGY↗

DOE OSTI2022

Cubic-meter scale laboratory fault re-activation experiments to improve the understanding of induced seismicity risks

Abstract To understand fluid induced seismicity, we have designed a large-scale laboratory experiment consisting of a one-cubic-meter sandstone with an artificial fault cut and fluid-injection boreholes. The sandstone block is assembled in a true triaxial loading frame and equipped with 38 piezoelectric sensors to locate and characterise acoustic emission events. The differential stress on the artificial fault is increased in stages to bring it towards a critically stressed state. After each stage of differential stress increase, fluids are injected at low pressures through boreholes to test the potential of fault re-activation. In addition, a high-pressure injection was conducted that created a hydraulic fracture from the injection borehole towards the artificial fault. The newly generated fluid pathway resulted in an activation of the complete block through a stick–slip movement. We compare acoustic emission measurements from the laboratory experiment with seismicity observations from the field-scale CO 2 injection at Decatur, Illinois, U.S., and conclude that the existence of fluid pathways plays a decisive role for the potential of induced seismicity.

47 OTHER INSTRUMENTATION↗

DOE OSTI2021

Network analysis of earthquake ground motion spatial correlation: a case study with the San Jacinto seismic nodal array

SUMMARY The spatial correlation of earthquake ground motion intensity can be measured from strong motion data; however, the data used in past studies is sparsely sampled in space, and only the interstation distance was considered as a correlation variable. These limitations mean that we have only weak constraints on the true correlation structure of ground motion and that potentially important aspects of spatial correlation are unconstrained. In this study, we combine a large-N seismic array and graph analytics to explore this issue at a local scale using small local and regional earthquakes. Our result suggests site conditions, and how they interact with the incident seismic wavefield, strongly condition the spatial correlation of ground motion. Future progress in characterizing ground motion spatial variability will require dense wavefield measurements, either through nodal deployments, or perhaps distributed acoustic sensing measurements, of seismic wavefields. Aftershock sequences of major earthquakes would provide particularly data-rich targets of opportunity.

58 GEOSCIENCES↗

DOE OSTI2022

Imaging the crust and uppermost mantle structure of Portugal (West Iberia) with seismic ambient noise

SUMMARY We present a new high-resolution 3-D shear wave velocity (Vs) model of the crust and uppermost mantle beneath Portugal, inferred from ambient seismic noise tomography. We use broad-band seismic data from a dense temporary deployment covering the entire Portuguese mainland between 2010 and 2012 in the scope of the WILAS project. Vertical component data are processed using phase correlation and phase weighted stack to obtain empirical Green functions (EGFs) for 2016 station pairs. Further, we use a random sampling and subset stacking strategy to measure robust Rayleigh-wave group velocities in the period range 7–30 s and associated uncertainties. The tomographic inversion is performed in two steps: First, we determine group-velocity lateral variations for each period. Next, we invert them at each grid point using a new trans-dimensional inversion scheme to obtain the 3-D shear wave velocity model. The final 3-D model extends from the upper crust (5 km) down to the uppermost mantle (60 km) and has a lateral resolution of ∼50 km. In the upper and middle crusts, the Vs anomaly pattern matches the tectonic units of the Variscan Massif and Alpine basins. The transition between the Lusitanian Basin and the Ossa Morena Zone is marked by a contrast between moderate- and high-velocity anomalies, in addition to two arched earthquake lineations. Some faults, namely, the Manteigas–Vilariça–Bragança fault and the Porto–Tomar–Ferreira do Alentejo fault, have a clear signature from the upper crust down to the uppermost mantle (60 km). Our 3-D shear wave velocity model offers new insights into the continuation of the main tectonic units at depth and contributes to better understanding the seismicity of Portugal.

Silveira, Graça (ORCID:0000000221102554)↗

DOE OSTI2023

Emulation of seismic-phase traveltimes with machine learning

SUMMARY We present a machine learning (ML) method for emulating seismic-phase traveltimes that are computed using a global-scale 3-D earth model and physics-based ray tracing. Accurate traveltime predictions based on 3-D earth models are known to reduce the bias of event location estimates, increase our ability to assign phase labels to seismic detections and associate detections to events. However, practical use of 3-D models is challenged by slow computational speed and the unwieldiness of pre-computed lookup tables that are often large and have prescribed computational grids. In this work, we train a ML emulator using pre-computed traveltimes, resulting in a compact and computationally fast way to approximate traveltimes that are based on a 3-D earth model. Our model is trained using approximately 850 million P-wave traveltimes that are based on the global LLNL-G3D-JPS model, which was developed for more accurate event location. The training-set consists of traveltimes between 10 393 global seismic stations and randomly sampled event locations that provide a prescribed, distance-dependent geographic sample density for each station. Prediction accuracy is dependent on event-station distance and whether the station was included in the training set. For stations included in the training set the mean absolute deviation (MAD) of the difference between traveltimes computed using ray tracing through the 3-D model and the ML emulator for local, regional, and teleseismic distances are 0.090, 0.125 and 0.121 s, respectively. For tested station locations not included in the training set, MAD values for the three distance ranges increase to 0.173, 0.219 and 0.210 s, respectively. Empirical traveltime residuals for a global reference data are indistinguishable when ML emulation or the 3-D model is used to compute traveltimes. This result holds regardless of whether the recording station is used in ML training or not.

58 GEOSCIENCES↗

DOE OSTI2025

Divide and conquer: separating the two probabilities in seismic phase picking

There are two fundamental probabilities in the seismic phase picking process—the probability of the existence of a seismic phase (detection probability) and the probability associated with the phase arrival time estimation (timing probability). The nearly ubiquitous approach in developing deep learning phase picking models is to use a kernel, such as a truncated Gaussian, to mask the labelled phase arrival time and train a segmentation model. Once a model is trained, the times of the peaks in the output are taken as phase arrival times (picks), and the height of the peaks are taken as ‘probability’ of the picks. Here, we show that this ‘probability’ represents neither the detection nor the timing probability because this approach forces the output to follow the shape of the kernel. We introduce an approach using two models to estimate these two distinct probabilities. We use a binary classifier with a calibrated confidence to address the detection probability and a multiclass classifier to obtain a probability mass function to address the timing probability. This new approach can make the deep learning-based phase picking process more interpretable and provide options to logically control seismic monitoring workflows.

58 GEOSCIENCES↗

DOE OSTI2026

Geophysical Monitoring Program Core Seismic Libraries

The Geophysical Monitoring Program Core Seismic Libraries contain a broad array of seismic and waveform analysis libraries useful for a wide range of seismic analysis processes. These include signal processing, travel time calculation, waveform correlation, array processing, geometry and geodesy, and plotting libraries.

Barno, JustinG [Lawrence Livermore National Labora↗

DOE OSTI2022

Host-region parameters for an adjustable model for crustal earthquakes to facilitate the implementation of the backbone approach to building ground-motion logic trees in probabilistic seismic hazard analysis

The backbone approach to constructing a ground-motion logic tree for probabilistic seismic hazard analysis (PSHA) can address shortcomings in the traditional approach of populating the branches with multiple existing, or potentially modified, ground-motion models (GMMs) by rendering more transparent the relationship between branch weights and the resulting distribution of predicted accelerations. To capture epistemic uncertainty in a tractable manner, there are benefits in building the logic tree through the application of successive adjustments for differences in source, path, and site characteristics between the host region of the selected backbone GMM and the target region for which the PSHA is being conducted. The implementation of this approach is facilitated by selecting a backbone GMM that is amenable to such host-to-target adjustments for individual source, path, and site characteristics. The NGA-West2 GMM of Chiou and Youngs (CY14) has been identified as a highly adaptable model for crustal seismicity that is well suited to such adjustments. Rather than using generic source, path, and site characteristics assumed appropriate for the host region, the final suite of adjusted GMMs for the target region will be better constrained if the host-region parameters are defined specifically on the basis of their compatibility with the CY14 backbone GMM. To this end, making use of a recently developed crustal shear-wave velocity profile consistent with CY14, we present an inversion of the model to estimate the key source and path parameters, namely the stress parameter and the anelastic attenuation. With these outputs, the effort in constructing a ground-motion logic tree for any PSHA dealing with crustal seismicity can be focused primarily on the estimation of the target-region characteristics and their associated uncertainties. The inversion procedure can also be adapted for any application in which different constraints might be relevant.

Stafford, Peter J.↗

DOE OSTI2021

Seismically Detecting Nuclear Reactor Operations Using a Power Spectral Density (PSD) Misfit Detector

To explore the ability to indirectly detect and attribute various operations conducted at a nuclear reactor using waveform data, we investigated the seismic signals recorded near the High Flux Isotope Reactor (HFIR) located at Oak Ridge National Laboratory in Oak Ridge, Tennessee. Specifically, we processed seismic data collected from a single seismoacoustic station, WACO, near the HFIR facility, and employed a power spectral density misfit detector to identify signals of interest and associate the detections with operational events. Initial results suggest that this method provides a promising means of regularly detecting at least 19 unique operations. Furthermore, with additional station deployment and more comprehensive data logs, we anticipate that future analysis will offer an additional means to seismically monitor nuclear reactors (such as HFIR) health and performance more accurately.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗