Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “seismic sources”

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 181 records · Page 10

Challenges and Opportunities in Turning Large U.S. Archives of Analog Seismograms into a Modern Usable Resource

Abstract Quantitative seismology is based firmly on the analysis of actual ground motions, and the transition to digital recording in the 1980s enabled sophisticated new capabilities to extract useful results from waveforms. With some effort, these tools can also be applied to analog records. Focusing on assets available within U.S. institutions, we review the necessary steps and the challenges in enabling “data rescue”—that is, preserving the scientific information latent in large analog seismogram archives and making it usable. They include: determining what assets are available (the analog seismogram archives held by various institutions, with associated metadata on instrument responses, station locations, and timing information); developing a consensus on the top level of a triage process (which analog records most definitely should be rescued?); deciding the level of quality needed in copying original seismograms to media suitable for digitizing; assessing the relative merits of scanning and digitizing; and, the need for a community service in distributing scans and digital records, as they accumulate. The necessary level of effort can benefit from practical experience. For example, specific studies have used digitized versions of analog recordings to model earthquake sources and assess seismic hazard. Other studies have used them to gain experience with nuclear explosion signals recorded at regional distances, noting that regional signals enable explosions to be monitored down to levels much lower than those attainable teleseismically. The opportunities presented by large archives of analog seismograms include the insights they present to current and future seismologists studying earthquakes and explosions, into the practical areas of assessing seismic hazard, monitoring for test ban compliance down to low explosion yields—and prompt characterization of actual explosions should they occur, as well the traditional academic pursuit of a better understanding of earthquake physics.

Geochemistry & Geophysics↗

Origin and evolution of the lunar surface - The major questions remaining

The major factors in the evolution of the lunar surface have not been determined yet. Huge lava flows and lunar differentiation, though commonly assumed, is in discord with much of the evidence. The alternative is for most of the surface to represent the last stages of accretion of the moon only, with the chemical differentiation having taken place previously in the source material. Radar, seismic, surface exposure, and mascon evidence can then be accounted for. A large-scale surface transport mechanism of soil must then have been present.

Gold, T.↗

Statistical description of tectonic motions

The behavior of stochastic processes was studied whose power spectra are described by power-law behavior. The details of the analysis and the conclusions that were reached are presented. This analysis was extended to compare detection capabilities of different measurement techniques (e.g., gravimetry and GPS for the vertical, and seismometers and GPS for horizontal), both in general and for the specific case of the deformations produced by a dislocation in a half-space (which applies to seismic of preseismic sources). The time-domain behavior of power-law noises is also investigated.

Agnew, Duncan Carr↗

Coseismic fault slip associated with the 1992 M(sub w) 6.1 Joshua Tree, California, earthquake: Implications for the Joshua Tree-Landers earthquake sequence

Coseismic surface deformation associated with the M(sub w) 6.1, April 23, 1992, Joshua Tree earthquake is well represented by estimates of geodetic monument displacements at 20 locations independently derived from Global Positioning System and trilateration measurements. The rms signal to noise ratio for these inferred displacements is 1.8 with near-fault displacement estimates exceeding 40 mm. In order to determine the long-wavelength distribution of slip over the plane of rupture, a Tikhonov regularization operator is applied to these estimates which minimizes stress variability subject to purely right-lateral slip and zero surface slip constraints. The resulting slip distribution yields a geodetic moment estimate of 1.7 x 10(exp 18) N m with corresponding maximum slip around 0.8 m and compares well with independent and complementary information including seismic moment and source time function estimates and main shock and aftershock locations. From empirical Green's functions analyses, a rupture duration of 5 s is obtained which implies a rupture radius of 6-8 km. Most of the inferred slip lies to the north of the hypocenter, consistent with northward rupture propagation. Stress drop estimates are in the range of 2-4 MPa. In addition, predicted Coulomb stress increases correlate remarkably well with the distribution of aftershock hypocenters; most of the aftershocks occur in areas for which the mainshock rupture produced stress increases larger than about 0.1 MPa. In contrast, predicted stress changes are near zero at the hypocenter of the M(sub w) 7.3, June 28, 1992, Landers earthquake which nucleated about 20 km beyond the northernmost edge of the Joshua Tree rupture. Based on aftershock migrations and the predicted static stress field, we speculate that redistribution of Joshua Tree-induced stress perturbations played a role in the spatio-temporal development of the earth sequence culminating in the Landers event.

Bennett, Richard A.↗

Assessing Potential Geophysical and Environmental Impacts from Frequent Rocket Launch Missions at Kennedy Space Center

Kennedy Space Center (KSC) in Florida has been utilized for space missions over many years with a gradually increasing number of rocket launches. There have been multiple studies where air-coupled acoustic waves and infrasound originating from launched rockets were used for operational purposes, such as locating booster trajectories as a function of changing atmospheric conditions. However, no study has utilized the acoustoelastic waves as a signal source for subsurface seismic investigations. We conducted a dispersion analysis using the seismic energy recorded from the Artemis I rocket launch at KSC in November 2022, and from these results we generated depth-sensitivity kernels at different wave frequencies. The kernels were compared with sedimentary core data to verify boundaries of carbonate layers above the Floridan Aquifer System. Accumulative information of bedrock-sediment boundary across the sedimentary platform could especially provide geo-structural evidence that manifests the configuration of coastal features. Continuous dispersion analysis of seismic recordings from consecutive rocket launches also has potential to identify non-stationary environmental effects, such as reorientation of sedimentary structures and fluctuation of the saltwater/groundwater lens from gravitational tides, which may affect erosional susceptibility of coastal features.

Han Byul Woo↗

Assessing Potential Geophysical and Environmental Impacts from Frequent Rocket Launch Missions at Kennedy Space Center

Kennedy Space Center (KSC) in Florida has been utilized for space missions over many years with a gradually increasing number of rocket launches. There have been multiple studies where air-coupled acoustic waves and infrasound originating from launched rockets were used for operational purposes, such as locating booster trajectories as a function of changing atmospheric conditions. However, no study has utilized the acoustoelastic waves as a signal source for subsurface seismic investigations. We conducted a dispersion analysis using the seismic energy recorded from the Artemis I rocket launch at KSC in November 2022, and from these results we generated depth-sensitivity kernels at different wave frequencies. The kernels were compared with sedimentary core data to verify boundaries of carbonate layers above the Floridan Aquifer System. Accumulative information of bedrock-sediment boundary across the sedimentary platform could especially provide geo-structural evidence that manifests the configuration of coastal features. Continuous dispersion analysis of seismic recordings from consecutive rocket launches also has potential to identify non-stationary environmental effects, such as reorientation of sedimentary structures and fluctuation of the saltwater/groundwater lens from gravitational tides, which may affect erosional susceptibility of coastal features.

Rocket Launch Seismicity↗

The Seismic Signature of a High-Energy Density Physics Laboratory and Its Potential for Measuring Time-Dependent Velocity Structure

The Z Machine at Sandia National Laboratories is a pulsed power facility for high-energy density physics experiments that can shock materials to extreme temperatures and pressures through a focused energy release of up to ~ 25 MJ in < 100 nanoseconds. It has been in operation for more than two decades and conducts up to ~ 100 experiments, or “shots,” per year. Based on a set of 74 known shot times from 2018, we determined that Z Machine shots produce detectable ~ 3–17 Hz ground motion 12 km away at the Albuquerque Seismological Laboratory, New Mexico (ANMO), borehole seismograph, with peak signal at ~ 7 Hz. The known shot waveforms were used to create a three-component template, leading to the detection of 2339 Z Machine shots since 1998 through single-station cross-correlation. Local seismic magnitude estimates range from local magnitude (M L ) —2 to —1.3 and indicate that only a small fraction of the shot energy is transmitted by seismic phases observable at 12 km distance. The most recent major facility renovation, which was intended to decrease mechanical dissipation, is associated with an abrupt decrease in observed seismic amplitudes at ANMO despite stable maximum shot energy. The highly repetitive impulsive sources are well suited to coda-wave interferometry to investigate time-dependent velocity structures. Relative velocity variations (dv/v) show an annual cycle with amplitude of ~ 0.2%. Local minima are observed in the late spring, and dv/v increases through the summer monsoon rainfall, possibly reflecting patchy saturation as rainfall infiltrates near the eastern edge of the Albuquerque basin. Importantly, the cumulative results demonstrate that forensic seismology can provide insight into long-term operation of facilities such as pulsed-power laboratories, and that their recurring signals may be valuable for studies of time-dependent structure

54 ENVIRONMENTAL SCIENCES↗

Economic assessment of seismic monitoring for underground hydrogen storage

Underground hydrogen storage (UHS) plays a key role in the energy landscape. However, like other subsurface engineering technologies, UHS may cause leakage into the groundwater or atmosphere and possibly induce local seismicity. To reduce these risks, seismic monitoring could be a viable technique to track the UHS plume, detect leakages, and locate induced seismicity events. Seismic monitoring has been proposed to safely monitor UHS, but research in this area is still new and requires field studies. Lab and theoretical studies have demonstrated the validity of seismic monitoring for UHS. Therefore, it is imperative to analyze the economic feasibility of seismic monitoring for UHS. Hence, we develop a cost model and open-source Python code for seismic monitoring that considers types of seismometers, comprehensive operational scenarios, detection thresholds, and long-term leakage monitoring. A case study is further provided to validate the cost model on reservoir simulations of UHS. We find that the levelized cost for a 10-year operating UHS site will range on the order of ∼0.003 $\$$/kg. The methods developed in this study could also be applied to the monitoring of groundwater, gas, and/or wastewater injection.

08 HYDROGEN↗

Seismic moment tensor classification using elliptical distribution functions on the hypersphere

Discrimination of underground explosions from naturally occurring earthquakes and other anthropogenic sources is one of the fundamental challenges of nuclear explosion monitoring. In an operational setting, the number of events that can be thoroughly investigated by analysts is limited by available resources. The capability to rapidly screen out events that can be robustly identified as not being explosions is, therefore, of great potential benefit. Nevertheless, possible mis-classification of explosions as earthquakes currently limits the use of screening methods for verification of test-ban treaties. Moment tensors provide a physics-based classification tool for the characterization of different seismic sources and have enabled the advent of new techniques for discriminating between earthquakes and explosions. Following normalization and projection of their six-degree vectors onto the hypersphere, existing screening approaches use spherically symmetric metrics to determine whether any new moment tensor may have been an explosion. Here, we show that populations of moment tensors for both earthquakes and explosions are anisotropically distributed on the hypersphere. Distributions possessing elliptical symmetry, such as the scaled von Mises–Fisher distribution, therefore provide a better description of these populations than the existing spherically symmetric models. We describe a method that uses these elliptical distributions in combination with a Bayesian classifier to achieve successful classification rates of 99 per cent for explosions and 98 per cent for earthquakes using existing catalogues of events from the western United States. The 1983 May 5 Crowdie underground nuclear test and 2018 July 20 DAG-1 deep-borehole chemical explosion are the only two explosions out of 140 that are incorrectly classified. Application of the method to the 2006–2017 nuclear tests in the Democratic People’s Republic of Korea yields 100 per cent identification rates and we provide a simple routine MTid for general usage. The approach provides a means to rapidly assess the likelihood of an event being an explosion and can be built into monitoring workflows that rely on simultaneously assessing multiple different discrimination metrics.

58 GEOSCIENCES↗

Code for Calculating Regional Seismic Travel Time

The RSTT software computes predictions of the travel time of seismic energy traveling from a source to a receiver through 2.5D models of the seismic velocity distribution within the Earth. The two primary applications for the RSTT library are tomographic inversion studies and seismic event location calculations. In tomographic inversions studies, a seismologist begins with number of source-receiver travel time observations and an initial starting model of the velocity distribution within the Earth. A forward travel time calculator, such as the RSTT library, is used to compute predictions of each observed travel time and all of the residuals (observed minus predicted travel time) are calculated. The Earth model is then modified in some systematic way with the goal of minimizing the residuals. The Earth model obtained in this way is assumed to be a better model than the starting model if it has lower residuals. The other major application for the RSTT library is seismic event location. Given an Earth model, an initial estimate of the location of a seismic event, and some number of observations of seismic travel time thought to have originated from that event, location codes systematically modify the estimate of the location of the event with the goal of minimizing the difference between the observed and predicted travel times. The second application, seismic event location, is routinely implemented by the military as part of its effort to monitor the Earth for nuclear tests conducted by foreign countries.

BALLARD, SANFORD↗

Regional Seismic Travel Time v.3.2.0

The RSTT software computes predictions of the travel time of seismic energy traveling from a source to a receiver through 2.5D models of the seismic velocity distribution within the Earth. The two primary applications for the RSTT library are tomographic inversion studies and seismic event location calculations. In tomographic inversions studies, a seismologist begins with number of source-receiver travel time observations and an initial starting model of the velocity distribution within the Earth. A forward travel time calculator, such as the RSTT library, is used to compute predictions of each observed travel time and all of the residuals (observed minus predicted travel time) are calculated. The Earth model is then modified in some systematic way with the goal of minimizing the residuals. The Earth model obtained in this way is assumed to be a better model than the starting model if it has lower residuals. The other major application for the RSTT library is seismic event location. Given an Earth model, an initial estimate of the location of a seismic event, and some number of observations of seismic travel time thought to have originated from that event, location codes systematically modify the estimate of the location of the event with the goal of minimizing the difference between the observed and predicted travel times. The second application, seismic event location, is routinely implemented by the military as part of its effort to monitor the Earth for nuclear tests conducted by foreign countries.

BALLARD, SANFORD↗

Seismic monitoring of a small CO 2 injection using a multi-well DAS array: Operations and initial results of Stage 3 of the CO 2 CRC Otway project

Active time-lapse seismic is widely employed for monitoring CO 2 geosequestration due to its ability to track the distribution of fluids in space and time. However, standard 4D seismic monitoring suffers from several challenges, including high cost, disruption to other land uses, and, consequently, relatively large intervals between monitor surveys. Some of these challenges can be mitigated using permanently installed sources and receivers. Such an approach was tested at the CO 2 CRC Otway site by continuous offset VSP monitoring of 15,000 t of supercritical CO 2 injected into an aquifer 1,500 m deep with nine permanent seismic sources (surface orbital vibrators or SOVs) and five downhole fibre-optic receivers. This continuous monitoring is complemented by multi-well 4D VSP using a mobile vibroseis source and the same DAS receivers, which included one baseline and two monitor surveys after injection of 4,000 and 12,000 t of CO 2 . The continuous DAS-SOV monitoring detected an abrupt increase of travel times below the injection interval on the second day of injection (after injection of 300 t of CO 2 ) and tracked the growth of the areal CO 2 plume by mapping changes of reflection amplitudes. The plume is also detected by time-lapse changes of reflection amplitudes in multi-well 4D VSPs. The plume images obtained from continuous offset VSP and 4D VSP are broadly consistent with each other but with some differences due to differences in illumination, lateral variations of velocities and seismic anisotropy. Furthermore, these differences also serve as a measure of uncertainty of 4D VSP images.

4D VSP↗

Spatially Dependent Seismic Wavefield Scattering from an Underground Chemical Explosion: Analysis of the Source Physics Experiment Dry Alluvium Geology Large- N Array

Explosion sources have been observed to generate significant shear‐wave energy despite their isotropic nature. To investigate this phenomenon, we conduct an analysis of the seismic data collected as part of the Source Physics Experiment (SPE): Dry Alluvium Geology (DAG) and investigate the generation of shear‐wave energy via scattering. The data were produced by three underground chemical explosions and consist of three‐component seismograms, which were recorded by the DAG Large‐N array. Synthetic tests suggest that for the DAG experiments, small‐scale stochastic heterogeneities, defined as features with correlation lengths of 10–100s of meters, are more effective than large‐scale geologic structure (scales >1–10 km) at reproducing the scattering of explosion generated wavefields observed at DAG. We analyze the seismic data for spatially variable ratios between transversely and radially polarized seismic energy, and then estimate the mean free path of P and S waves. All analyses are conducted within a frequency band of 5–50 Hz. The ratio of transversely to radially polarized energy is the highest in the east and west portion of the Large‐N array. In addition, the magnitude of the estimated S‐wave mean free path is shorter in the eastern portion of the Large‐N array. This variation indicates that the eastern area of the DAG array is where more scattering is occurring, suggesting azimuthal dependence of P‐to‐P and P‐to‐S scattering. Furthermore, this azimuthal dependence of P‐to‐S scattering can have implications for explosion discrimination based on spectral ratios of seismic wave types, because the general assumption is that explosions do not generate shear‐wave energy. Synthetic tests modeling only larger‐scale geologic structure had lower transversely polarized energy (only four stations showing a transversely to radially polarized energy ratio greater than 1) and fewer stations (<10) displaying shorter (<300 m) mean free paths than what was observed in the DAG data results.

Darrh, Andréa Nicole↗

Fourier-DeepONet: Fourier-enhanced deep operator networks for full waveform inversion with improved accuracy, generalizability, and robustness

In this article, full waveform inversion (FWI) infers the subsurface structure information from seismic waveform data by solving a non-convex optimization problem. Data-driven FWI has been increasingly studied with various neural network architectures to improve accuracy and computational efficiency. Nevertheless, the applicability of pre-trained neural networks is severely restricted by potential discrepancies between the source function used in the field survey and the one utilized during training. Here, we develop a Fourier-enhanced deep operator network (Fourier-DeepONet) for FWI with the generalization of seismic sources, including the frequencies and locations of sources. Specifically, we employ the Fourier neural operator as the decoder of DeepONet, and we utilize source parameters as one input of Fourier-DeepONet, facilitating the resolution of FWI with variable sources. To test Fourier-DeepONet, we develop three new and realistic FWI benchmark datasets (FWI-F, FWI-L, and FWI-FL) with varying source frequencies, locations, or both. Our experiments demonstrate that compared with existing data-driven FWI methods, Fourier-DeepONet obtains more accurate predictions of subsurface structures in a wide range of source parameters. Moreover, the proposed Fourier-DeepONet exhibits superior robustness when handling data with Gaussian noise or missing traces and sources with Gaussian noise, paving the way for more reliable and accurate subsurface imaging across diverse real conditions.

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)↗

The feasibility of MT tipper data to monitor CO2 storage sites

Monitoring carbon storage sites using geophysical techniques is a critical component to the success and safety of storage programs. Currently, the primary methods of monitoring such sites are seismic and, to a much lesser extent, electromagnetics using active sources. The cost of such methods, especially seismic, can be prohibitively expensive. Natural source lectromagnetics, or magnetotellurics (MT), represents a low-cost, and underutilized method with the potential to aid CO2 monitoring efforts. Specifically, the tipper of MT data gives insight to the dimensionality of the subsurface and is able to detect the expansion front of a CO2 plume in a saline reservoir. We analyze the feasibility of using the tipper to monitor two shallow CO2 plumes and conclude that the tipper may be a suitable method for long-term monitoring.<br>

Kohnke, Colton↗

Ambient field seismology in critical zone hydrological sciences

Passive ambient noise monitoring is an emerging tool in environmental seismology, leveraging the ambient seismic field to assess temporal variations in shallow subsurface properties. This review focuses on the potential and challenges of using scattered coda waves from noise correlation functions to monitor critical zone dynamics. The sensitivity of seismic velocities to various environmental factors, including precipitation, snowmelt, atmospheric pressure, and groundwater fluctuations, underscores the method’s versatility. While coda waves excel in detecting subtle changes due to their scattered nature, ballistic waves provide higher spatial resolution, albeit with challenges in source stability. Advances in seismic sensing, including distributed acoustic sensing and low-cost geophone networks, have enabled high-resolution monitoring of hydrological processes, subsurface deformation, and seismic hazards. Integrating seismic data with hydrological models provides insights into water storage, pore pressure changes, and soil moisture dynamics. However, limitations in spatial resolution, calibration with ground truth data, and coupled effects between environmental factors remain key challenges. This review emphasizes the importance of interdisciplinary approaches in refining methodologies, enhancing sensor deployments, and addressing data gaps. Passive seismic monitoring offers opportunities to understand critical zone processes and their broader impacts on seismic hazards and environmental sustainability.

58 GEOSCIENCES↗