Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Rayleigh waveforms”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Yield Estimates of three Historical Atmospheric Nuclear Explosions at Lop Nor, from Reduced Order Models of Regionally Recorded Rayleigh Waveforms and One Seismic Station: A Brief Communication

We use a reduced order model (ROM) for Rayleigh waveforms sourced by large atmospheric explosions along with data collected from a single seismic station (TLG) to estimate the yields of three atmospheric nuclear tests with order-megaton (MT) yields. These historical tests were conducted by China between 1973 and 1980 (CHIC 15, CHIC 16, and CHIC 26). We use our ROM to estimate yields of 2.34 MT (CHIC 15), 0.51 MT (CHIC 16), and 1.15 (CHIC 26) in the 20.5 s period band that are consistent with publicly accepted values of 3 MT (CHIC 15), 0.50 MT (CHIC 16), and 0.730 MT (CHIC 26).

58 GEOSCIENCES↗
DOE OSTI2025

Designing Yield Grids for Airborne Nuclear Explosions to Maintain Decrements in Rayleigh Waveform Amplitudes

This document details a method for selecting grid points that sample the yields of airborne, low yield nuclear explosions that constrains the relative size of seismic waveforms sourced by explosions that are parameterized by this grid. In detail, we impose this constraint so that the vertical component of ground displacement that is sourced by explosions with a fixed height of burst (HoB) at grid point k change by a fixed factor of a when compared to the waveform amplitudes sourced by explosions at grid point k + 1, as observed at a fixed standoff distance. This gridding thereby preserves amplitude ratios along the grid so that they are agnostic to k. We assert that this method is useful in tests that apply algorithms against synthetic or augmented datasets for which the performance of the algorithm must be measured against source size, and waveform amplitudes must change uniformly (e.g., Berg et al, 2025).

58 GEOSCIENCES↗
NASA NTRS1985

Fade margin calculation for channels impaired by Rician fading

Excess path loss due to multipath severely restricts the performance of power limited mobile networks such as those using satellite-aided links. To reduce multipath related losses, the higher elevation angle of the spacecraft can be exploited by utilizing mobile antennas which reduce the strength of the multipath reflections in favor of the line-of-sight signal. The presence of a strong and stable path in a fading link will change the envelope statistics of the received waveform from Rayleigh to a more favorable Rician distribution. It is determined that the excess path loss, or fade margin, of a Rician channel when coherent detection of binary phase shift keying (BPSK) or quaternary phase shift keying (QPSK) signals is considered. The results are presented parametrically such that they can be applied to a wide range of propagation characteristics from heavy fading to nonfading situations. Furthermore, similar results are also given for the case where only limited coverage is provided.

Davarian, F.↗
DOE OSTI2021

Rayleigh Wave Propagation in the Bighorn Mountains Region, Wyoming

Short-period Rayleigh waves, Rg , provide strong constraints on the depth of shallow seismic events and are of interest for monitoring small explosions. Characterizing the seismic sources that generate Rg requires an understanding of how shallow crustal structure affects Rayleigh wave propagation. Here, in support of these efforts, this study utilizes observed waveforms from small shallow explosions recorded on temporary seismic network deployments in the Bighorn region, Wyoming. We study regional near-surface geology by measuring changes in surface-wave amplitude and polarization during propagation through basins, foothills, and mountains. We develop additional insight by carrying out surface-wave eigenfunction analyses and numerical-wave simulations, which together reproduce many characteristics seen in the observed waveforms. Our results show how sedimentary basins in the Bighorn region allow for amplified prograde-polarized higher-mode and retrograde-polarized fundamental-mode Rayleigh waves, whereas adjacent mountains only support retrograde motion. These different modes provide distinct constraints on the Earth structure and source characteristics, potentially enabling targeted inversions in future studies. Finally, our findings provide insight into Rg propagation through complex near-surface geology, improving our understanding of shallow propagation and source effects that are relevant to explosion monitoring efforts.

58 GEOSCIENCES↗
DOE OSTI2022

High-Precision Characterization of Seismicity from the 2022 Hunga Tonga-Hunga Ha'apai Volcanic Eruption

The earthquake swarm accompanying the January 2022 Hunga Tonga-Hunga Ha'apai (HTHH) volcanic eruption includes a large number of posteruptive moderate-magnitude seismic events and presents a unique opportunity to use remote monitoring methods to characterize and compare seismic activity with other historical caldera-forming eruptions. We compute improved epicentroid locations, magnitudes, and regional moment tensors of seismic events from this earthquake swarm using regional to teleseismic surface-wave cross correlation and waveform modeling. Precise relative locations of 91 seismic events derived from 59,047 intermediate-period Rayleigh- and Love-wave cross-correlation measurements collapse into a small area surrounding the volcano and exhibit a southeastern time-dependent migration. Regional moment tensors and observed waveforms indicate that these events have a similar mechanism and exhibit a strong positive compensated linear vector dipole component. Precise relative magnitudes agree with regional moment tensor moment magnitude ($M_w$) estimates while also showing that event sizes and frequency increase during the days after the eruption followed by a period of several weeks of less frequent seismicity of a similar size. The combined information from visual observation and early geologic models indicate that the observed seismicity may be the result of a complex series of events that occurred after the explosive eruption on 15 January, possibly involving rapid resupply of the magma chamber shortly after the eruption and additional faulting and instability in the following weeks. In addition, we identify and characterize an $M_w$ 4.5 event five days before the paroxysmal explosion on 15 January, indicating that additional seismic events preceding the main eruption could have been identified with improved local monitoring. As a result, our analysis of the HTHH eruption sequence demonstrates the value of potentially utilizing teleseismic surface-wave cross correlation and waveform modeling methods to assist in the detailed analysis of remote volcanic eruption sequences.

58 GEOSCIENCES↗
NASA NTRS2012

Three-Dimensional Numerical Simulations of Equatorial Spread F: Results and Observations in the Pacific Sector

A three-dimensional numerical simulation of plasma density irregularities in the postsunset equatorial F region ionosphere leading to equatorial spread F (ESF) is described. The simulation evolves under realistic background conditions including bottomside plasma shear flow and vertical current. It also incorporates C/NOFS satellite data which partially specify the forcing. A combination of generalized Rayleigh-Taylor instability (GRT) and collisional shear instability (CSI) produces growing waveforms with key features that agree with C/NOFS satellite and ALTAIR radar observations in the Pacific sector, including features such as gross morphology and rates of development. The transient response of CSI is consistent with the observation of bottomside waves with wavelengths close to 30 km, whereas the steady state behavior of the combined instability can account for the 100+ km wavelength waves that predominate in the F region.

ALTAIR↗
DOE OSTI2021

Report on the LLNL Global Full-waveform Inversion Workflow and Progress

LLNL has developed the SPiRaL global seismic tomography model based upon millions of body wave travel times and surface wave dispersion curves in the period range of 25-200 seconds (Simmons et al., 2021). The motivation to construct the SPiRaL model (which stands for S-, P-, Rayleigh, and Love waves) is based upon the guidance provided by a 3-D model workshop held in Berkeley in 2007 and documented in a paper by Zucca et al. (2009). It was concluded in Zucca et al. (2009) that a “practical and obtainable” model is one that predicts body wave travel times for event location while also predicting surface wave velocities for moment tensor estimation. We have largely followed that guidance which has led to the construction of SPiRaL. Looking forward, we now seek to improve the global model through full waveform inversion (FWI). This necessitates the development/customization of an efficient and flexible modeling workflow to automate the process using LLNL’s supercomputing platforms. This report briefly outlines the work that being performed at LLNL to develop the necessary tools and framework to update the SPiRaL global model with full waveform information.

58 GEOSCIENCES↗
DOE OSTI2020

A multi-terawatt two-color beam for high-power field-controlled nonlinear optics

Two-color laser beams are instrumental in light-field control and enhancement of high-order harmonic, spectral supercontinuum, and terahertz radiation generated in gases, plasmas, and solids. Here we demonstrate a multi-terawatt two-color beam produced using a relativistic plasma mirror, with 110 mJ at 800 nm and 30 mJ at 400 nm. Both color components have high spatial quality and can be simultaneously focused, provided that the plasma mirror lies within a Rayleigh range of the driving fundamental beam. Favorable scaling of second-harmonic generation by plasma mirrors at relativistic intensities suggests them as an excellent tool for multi-color waveform synthesis beyond the petawatt level.

42 ENGINEERING↗
DOE OSTI2025

Simulating Isotropic Rayleigh Wave Horizontal to Vertical Ratios at Rock Valley

Elastic anisotropy in the subsurface can manifest from pre-existing structure or stress-induced cracking – both phenomena are of interest at the Rock Valley site in Nevada. The azimuthal dependence of Rayleigh wave ellipticity, also known as H/V, can reveal elastic anisotropy, but the effect of isotropic variability in H/V must first be accounted for. To that end, we simulate synthetic seismograms from a ring of sources in a fully 3-D isotropic velocity model that captures the geology at Rock Valley. We measure H/V from these waveforms and compare them with H/V predictions obtained using 1D models extracted at specific locations from the 3D model. We find that an average of the azimuthal H/V measurements from the 3-D simulations and 1-D predictions agree, but that the azimuthal variability of H/V is a large fraction of the mean value. These results estimate the noise character in azimuthal H/V so it can be factored out of future anisotropic simulations.

58 GEOSCIENCES↗
NASA NTRS1983

Waveform inversion of mantle Love waves: The born seismogram approach

Normal mode theory, extended to the slightly laterally heterogeneous Earth by the first-order Born approximation, is applied to the waveform inversion of mantle Love waves (200-500 sec) for the Earth's lateral heterogeneity at l=2 and a spherically symmetric anelasticity (Q sub mu) structure. The data are from the Global Digital Seismograph Network (GDSN). The l=2 pattern is very similar to the results of other studies that used either different methods, such as phase velocity measurements and multiplet location measurements, or a different data set, such as mantle Rayleigh waves from different instruments. The results are carefully analyzed for variance reduction and are most naturally explained by heterogeneity in the upper 420 km. Because of the poor resolution of the data set for the deep interior, however, a fairly large heterogeneity in the transition zones, of the order of up to 3.5% in shear wave velocity, is allowed. It is noteworthy that Love waves of this period range can not constrain the structure below 420 km and thus any model presented by similar studies below this depth are likely to be constrained by Rayleigh waves (spheroidal modes) only.

Tanimoto, T.↗
NASA NTRS1984

Waveform inversion of mantle Love waves - The Born seismogram approach

Normal mode theory, extended to the slightly laterally heterogeneous earth by the first-order Born approximation, is applied to the waveform inversion of mantle Love waves (200-500 sec) for the earth's lateral heterogeneity at l = 2 and a spherically symmetric anelasticity (Q sub mu) structure. The data are from the Global Digital Seismograph Network (GDSN). The l = 2 pattern is very similar to the results of other studies that used either different methods, such as phase velocity measurements and multiplet location measurements, or a different data set, such as mantle Rayleigh waves from different instruments. The results are carefully analyzed for variance reduction and are most naturally explained by heterogeneity in the upper 420 km. Because of the poor resolution of the data set for the deep interior, however, a fairly large heterogeneity in the transition zones, of the order of up to 3.5 percent in shear wave velocity, is allowed. It is noteworthy that Love waves of this period range can not constrain the structure below 420 km and thus any model presented by similar studies below this depth are likely to be constrained by Rayleigh waves (spheroidal modes) only.

Tanimoto, T.↗
NASA NTRS1982

A large normal-fault earthquake at the junction of the Tonga trench and the Louisville ridge

Long-period vertical-component Rayleigh waves are inverted in order to determine the source mechanism of the October 10, 1977 earthquake that occurred in the oceanic plate at the junction of the Tonga-Kermadec trench systems with the aseismic Louisville ridge. The cause was predominantly normal faulting on a plane striking roughly parallel to the trench, with a seismic moment of 1.7 x 10 to the 27th dyn cm. A focal depth of 20 km is determined by waveform modeling, but the actual rupture may have extended to 30 or 40 km. Two sources separated by 16 s comprised the event, which experienced an inferred rupture velocity of 3.5 km/sec. The interpretation that the earthquake was caused by gravitational pull due to the sinking slab implies that the Louisville ridge causes some degree of local decoupling between the plates. This event may be associated with the breakup of the Osbourn seamount. Alternatively, the earthquake may have resulted from tensional plate bending stress, as implied by its relatively shallow depth.

Eissler, H.↗
DOE OSTI2023

Accelerating full-waveform inversion using source stacking: synthetic experiments at the global scale in a realistic 3-D earth model

SUMMARY The spectral element method is currently the method of choice for computing accurate synthetic seismic wavefields in realistic 3-D earth models at the global scale. However, it requires significantly more computational time, compared to normal mode-based approximate methods. Source stacking, whereby multiple earthquake sources are aligned on their origin time and simultaneously triggered, can reduce the computational costs by several orders of magnitude. We present the results of synthetic tests performed on a realistic radially anisotropic 3-D model, slightly modified from model SEMUCB-WM1 with three component synthetic waveform ‘data’ for a duration of 10 000 s, and filtered at periods longer than 60 s, for a set of 273 events and 515 stations. We consider two definitions of the misfit function, one based on the stacked records at individual stations and another based on station-pair cross-correlations of the stacked records. The inverse step is performed using a Gauss–Newton approach where the gradient and Hessian are computed using normal mode perturbation theory. We investigate the retrieval of radially anisotropic long wavelength structure in the upper mantle in the depth range 100–800 km, after fixing the crust and uppermost mantle structure constrained by fundamental mode Love and Rayleigh wave dispersion data. The results show good performance using both definitions of the misfit function, even in the presence of realistic noise, with degraded amplitudes of lateral variations in the anisotropic parameter ξ. Interestingly, we show that we can retrieve the long wavelength structure in the upper mantle, when considering one or the other of three portions of the cross-correlation time series, corresponding to where we expect the energy from surface wave overtone, fundamental mode or a mixture of the two to be dominant, respectively. We also considered the issue of missing data, by randomly removing a successively larger proportion of the available synthetic data. We replace the missing data by synthetics computed in the current 3-D model using normal mode perturbation theory. The inversion results degrade with the proportion of missing data, especially for ξ, and we find that a data availability of 45 per cent or more leads to acceptable results. We also present a strategy for grouping events and stations to minimize the number of missing data in each group. This leads to an increased number of computations but can be significantly more efficient than conventional single-event-at-a-time inversion. We apply the grouping strategy to a real picking scenario, and show promising resolution capability despite the use of fewer waveforms and uneven ray path distribution. Source stacking approach can be used to rapidly obtain a starting 3-D model for more conventional full-waveform inversion at higher resolution, and to investigate assumptions made in the inversion, such as trade-offs between isotropic, anisotropic or anelastic structure, different model parametrizations or how crustal structure is accounted for.

Geochemistry & Geophysics↗
NASA NTRS1999

Reflections of AE Waves in Finite Plates: Finite Element Modeling and Experimental Measurements

The capability of a three-dimensional dynamic finite element method for predicting far-field acoustic emission signals in thin plates of finite lateral extent, including their reflections from the plate edges, was investigated. A lead break (Hsu-Neilsen) source to simulate AE was modeled and used in the experimental measurements. For the thin plate studied, the signals were primarily composed of the lowest order symmetric (S0) and antisymmetric (A0) Lamb modes. Experimental waveforms were detected with an absolutely calibrated, wideband, conical element transducer. The conditions of lead fractures both on the surface of the plate as well as on the edge of the plate were investigated. Surface lead breaks preferentially generate the A0 mode while edge lead breaks generate the S0 mode. Reflections of developed plate waves from both normal and oblique incidence angles were evaluated. Particularly interesting for the case of the lead break on the plate edge were S0 waves produced by the interaction of a Rayleigh wave with the plate corner and by a bulk shear wave mode converting at the side edge. The Rayleigh wave, in this case, propagated along the specimen edge. For all cases considered, the experimental measurements were in good agreement with the predictions of the finite element model.

Prosser, W. H.↗
NASA NTRS2011

Multicolor Photometric Observation of Lightning from Space: Comparison with Radio Measurements

This study evaluates the effectiveness of spectrophotometric measurements from space in revealing properties of lightning flash. The multicolor optical waveform data obtained by FORMOSAT-2/Imager of Sprites and Upper Atmospheric Lightning (ISUAL) were analyzed in relation to National Lightning Detection Network (NLDN), North Alabama Lightning Mapping Array (LMA). As of July 2011, we found six lightning events which were observed by ISUAL and North Alabama LMA. In two of these events, NLDN showed clear positive cloud-to-ground (CG) discharges with peak current of +139.9 kA and +41.6 kA and, around that time, LMA showed continuous intra-cloud (IC) leader activities at 4-6 km altitudes. ISUAL also observed consistent optical waveforms of the IC and CG components and, interestingly, it was found that the blue/red spectral ratio clearly decreased by a factor of 1.5-2.5 at the time of CG discharges. Other four lightning events in which NLDN did not detect any CG discharges were also investigated, but such a feature was not found in any of these cases. These results suggest that the optical color of CG component is more reddish than that of IC component and we explain this as a result of more effective Rayleigh scattering in blue light emissions coming from lower-altitude light source. This finding suggests that spectral measurements could be a new useful technique to characterize ICs and CGs from space. In this talk, we will also present a result from lightning statistical analysis of ISUAL spectrophotometric data and ULF magnetic data.

Adachi, Toru↗
NASA NTRS1993

An analysis of nearfield normal mode amplitude anomalies of the Landers earthquake

The 1992 Landers earthquake (M(sub w) = 7.3) occurred in the middle of the TERRAscope network. Long-period Rayleigh waves recorded at the TERRAscope stations (delta less than or = 3 deg) after traveling around the Earth show large amplitude anomalies, one order of magnitude larger than spherical Earth predictions up to a period of about 600 s. The ground motions over the epicentral region at and after the arrival of R4-5 are in phase at all stations. These observations are inconsistent with the nearly vertical strike slip mechanism of the Landers earthquake. Synthetic seismograms for a rotating, elliptic, and laterally heterogeneous Earth model calculated by the variational method agree well with the observed waveforms. Calculations for various 3D Earth models demonstrate that the amplitudes are very sensitive to the large scale aspherical structure in the crust and the mantle. The anomalies for modes shorter than 300 s period can be explained by lateral heterogeneity shallower than the upper mantle. Rotation of the Earth and lower mantle heterogeneity are required to explain mode amplitudes at longer periods. Current whole mantle seismic tomographic models can fully explain the observed amplitudes longer than 300 s. To assess the effect of the high order lateral heterogeneity in the mantle more precise estimate of the crustal correction is required.

Watada, Shingo↗