Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “earthquake”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

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↗
DOE OSTI2026

A Nonergodic Ground-Motion Model for the San Francisco Bay Area for Small-Magnitude Earthquakes

ABSTRACT Recently, generative models have become a computationally efficient alternative to physics-based numerical simulations of ground motions. Neural networks can learn from existing ground-motion data to generate unobserved ground-motion data at new source and site locations. A key challenge with generative models is ensuring that predicted ground motions remain within a physically realistic range. For this purpose, we developed an empirical, nonergodic ground-motion model (GMM) for small-magnitude earthquakes in the San Francisco Bay area based on about 5000 recordings per component for Mw ≤ 4 earthquakes. The nonergodic GMM predicts spatially varying median source, site, and path effects for both the Fourier amplitude spectrum (FAS) and the Fourier phase derivative (a proxy for duration), as well as the corresponding epistemic uncertainty for each term. For FAS, our model shows above-average source and site effects in the western part of the region and below-average effects in the eastern part, with regional effects exhibiting larger spatial correlation lengths with increasing frequency. For duration, the source term is negligible for small-magnitude earthquakes, and the site term leads to site-specific variations up to 5 s. Path effects for FAS and duration depend on the source–site pair and are extrapolated spatially using recent methods for path-effect modeling. The aleatory variability of the within-site within-path residuals is similar to the variability found in previous studies for other regions. The nonergodic model provides two key contributions: first, median adjustment terms that are transferable to larger magnitude earthquakes, further reducing aleatory variability in probabilistic seismic hazard analysis; second, region-specific criteria for validating machine learning-based ground-motion generators to evaluate whether synthetic ground motions exhibit physically realistic source, site, and path effects.

Lacour, Maxime↗
DOE OSTI2025

Probabilistic Seismic Hazard Analysis for Iraq Based on the Updated Earthquake Catalog (1900-2021) and Ground Motion Characteristics

Onur et al. (2017) compiled the first comprehensive earthquake catalog for Iraq, covering 1900 to 2009 within 26°–40°N latitude and 36°–51°E longitude. This catalog was utilized in a probabilistic seismic hazard assessment (PSHA) by Abdulnaby et al. (2020) to aid in updating Iraq’s building code seismic provisions. Recently, we have updated the earthquake catalog for Iraq by adding earthquakes recorded from 2010 to 2021 and directly calculating moment magnitude (Mw) for about 2,800 earthquakes using the coda envelope methodology and waveform data from the Mesopotamian Seismological Network (MPSN) in Iraq.

58 GEOSCIENCES↗
NASA NTRS1974

ERTS-1, earthquakes, and tectonic evolution in Alaska

In comparing seismicity patterns in Alaska with ERTS-1 imagery, it is striking to see the frequency with which earthquake epicenters fall on, or near, lineaments visible on the imagery. Often these lineaments prove to be tectonics faults which have been mapped in the field. But equally as often, existing geologic and tectonic maps show no evidence of these features. The remoteness and inaccessibility of most of Alaska is responsible, in large part, for the inadequacy of the mapping. ERTS-1 imagery is filling a vital need in providing much of the missing information, and is pointing out many areas of potential earthquake hazard. Earthquakes in central and south-central Alaska result when the northeastern corner of the north Pacific lithospheric plate underthrusts the continent. North of Mt. McKinley, the seismicity is continental in nature and of shallow origin, with earthquakes occurring on lineaments, and frequently at intersections of lineaments. The shallower events tend to align themselves with lineaments visible on the imagery.

Gedney, L.↗
NASA NTRS1980

Shallow moonquakes - How they compare with earthquakes

Of three types of moonquakes strong enough to be detectable at large distances - deep moonquakes, meteoroid impacts and shallow moonquakes - only shallow moonquakes are similar in nature to earthquakes. A comparison of various characteristics of moonquakes with those of earthquakes indeed shows a remarkable similarity between shallow moonquakes and intraplate earthquakes: (1) their occurrences are not controlled by tides; (2) they appear to occur in locations where there is evidence of structural weaknesses; (3) the relative abundances of small and large quakes (b-values) are similar, suggesting similar mechanisms; and (4) even the levels of activity may be close. The shallow moonquakes may be quite comparable in nature to intraplate earthquakes, and they may be of similar origin.

Nakamura, Y.↗
NASA NTRS1979

Earthquake Testing

During NASA's Apollo program, it was necessary to subject the mammoth Saturn V launch vehicle to extremely forceful vibrations to assure the moonbooster's structural integrity in flight. Marshall Space Flight Center assigned vibration testing to a contractor, the Scientific Services and Systems Group of Wyle Laboratories, Norco, California. Wyle-3S, as the group is known, built a large facility at Huntsville, Alabama, and equipped it with an enormously forceful shock and vibration system to simulate the liftoff stresses the Saturn V would encounter. Saturn V is no longer in service, but Wyle-3S has found spinoff utility for its vibration facility. It is now being used to simulate earthquake effects on various kinds of equipment, principally equipment intended for use in nuclear power generation. Government regulations require that such equipment demonstrate its ability to survive earthquake conditions. In upper left photo, Wyle3S is preparing to conduct an earthquake test on a 25ton diesel generator built by Atlas Polar Company, Ltd., Toronto, Canada, for emergency use in a Canadian nuclear power plant. Being readied for test in the lower left photo is a large circuit breaker to be used by Duke Power Company, Charlotte, North Carolina. Electro-hydraulic and electro-dynamic shakers in and around the pit simulate earthquake forces.

Source record↗
DOE OSTI2025

Automatic speech recognition predicts contemporaneous earthquake fault displacement

Abstract Significant progress has been made in probing the state of an earthquake fault by applying machine learning to continuous seismic waveforms. The breakthroughs were originally obtained from laboratory shear experiments and numerical simulations of fault shear, then successfully extended to slow-slipping faults. Here we apply the Wav2Vec-2.0 self-supervised framework for automatic speech recognition to continuous seismic signals emanating from a sequence of moderate magnitude earthquakes during the 2018 caldera collapse at the Kīlauea volcano on the island of Hawai’i. We pre-train the Wav2Vec-2.0 model using caldera seismic waveforms and augment the model architecture to predict contemporaneous surface displacement during the caldera collapse sequence, a proxy for fault displacement. We find the model displacement predictions to be excellent. The model is adapted for near-future prediction information and found hints of prediction capability, but the results are not robust. The results demonstrate that earthquake faults emit seismic signatures in a similar manner to laboratory and numerical simulation faults, and artificial intelligence models developed for encoding audio of speech may have important applications in studying active fault zones.

58 GEOSCIENCES↗
DOE OSTI2025

Rapid wavefield forecasting for earthquake early warning via deep sequence to sequence learning

We propose a deep learning model, WaveCastNet, to forecast high-dimensional wavefields. WaveCastNet integrates a convolutional long expressive memory architecture into a sequence-to-sequence forecasting framework, enabling it to model long-term dependencies and multiscale patterns in both space and time. By sharing weights across spatial and temporal dimensions, WaveCastNet requires significantly fewer parameters than more resource-intensive models such as transformers, resulting in faster inference times. Crucially, WaveCastNet also generalizes better than transformers to rare and critical seismic scenarios, such as high-magnitude earthquakes. Here, we show the ability of the model to predict the intensity and timing of destructive ground motions in real time, using simulated data from the San Francisco Bay Area. Furthermore, we demonstrate its zero-shot capabilities by evaluating WaveCastNet on real earthquake data. Our approach does not require estimating earthquake magnitudes and epicenters, steps that are prone to error in conventional methods, nor does it rely on empirical ground-motion models, which often fail to capture strongly heterogeneous wave propagation effects.

Geophysics↗
DOE OSTI2025

Factors That Influence Variability in Stress-Drop Measurements Using Spectral Decomposition and Spectral-Ratio Methods for the 2019 Ridgecrest Earthquake Sequence

Stress drop is a fundamental parameter related to earthquake source physics, but is hard to measure accurately. To better understand how different factors influence stress-drop measurements, we compare two different methods using the Ridgecrest stress-drop validation data set: spectral decomposition (SD) and spectral ratio (SR), each with different processing options. Here, we also examine the influence of spectral complexity on source parameter measurement. Applying the SD method, we find that frequency bandwidth and time-window length could influence spectral magnitude calibration, while depth-dependent attenuation is important to correctly map stress-drop variations. For the SR method, we find that the selected source model has limited influence on the measurements; however, the Boatwright model tends to produce smaller standard deviation and larger magnitude dependence than the Brune model. Variance reduction threshold, frequency bandwidth, and time-window length, if chosen within an appropriate parameter range, have limited influence on source parameter measurement. For both methods, wave type, attenuation correction, and spectral complexity strongly influence the result. The scale factor that quantifies the magnitude dependence of stress drop show large variations with different processing options, and earthquakes with complex source spectra deviating from the Brune-type source models tend to have larger scale factor than earthquakes without complexity. Based on these detailed comparisons, we make a few specific suggestions for data processing workflows that could help future studies of source parameters and interpretations.

58 GEOSCIENCES↗
DOE OSTI2026

Reply to “Comment on ‘The Reduced Detection Rate of Signals That Are Hidden by Earthquakes: Case Studies with Spotlight Detectors That Operate at Seismic Arrays,’ by Joshua D. Carmichael, Brent G. Delbridge, and Richard Alfaro-Diaz” by Paul G. Richards

We reply to a comment that Paul G. Richards (Richards, 2026) directed to an article by Carmichael et al. (2025). Richards expressed concern over two issues. The first issue relates to the size of the explosions discussed in the article. The second relates to the probability that earthquakes and explosions sourced near the same location coincide in time, by chance. We do not address the first issue because the original article does not assign significance to the size of the explosions from the point of view of an experimental party. We do respond to Richards’ comment on the second issue and claim that it can be explained through a clarification about conditional versus joint probabilities. Our present discussion, therefore, describes the difference between the conditional probability that a correlator fails, given that an explosion and earthquake coincidently occur in the same place, and the joint probability that a correlator fails while an earthquake and explosion also coincidentally occur in the same place. This latter, joint probability appears to be rare. We clarify that our original article did not treat the joint probability, only the conditional probability.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2025

Shear-Wave Radiation Patterns from Explosive and Earthquake Sources in Scattering, Heterogeneous Media

Distinguishing whether a seismic event is an earthquake or an explosion is a core problem in explosion monitoring. For simplistic models, earthquakes produce a predictable S-wavefield based on their radiation pattern, while an explosion produces no S-wavefield. However, observations from nuclear tests show that explosions can produce significant S-wave energy on both horizontal components. We perform numerical experiments using SPECFEM to constrain when an S-wavefield generated from an explosion might differ from one produced by an earthquake. We generate the S-wavefield for our explosions by placing the source location within a small region consisting of large velocity heterogeneities. For the 2D cases investigated here, a crucial condition is that the size of the region of heterogeneities must be much smaller (approximately 1/12th) than the minimum wavelength of the wavefield for our simulation setups. We investigate both isotropic heterogeneities and anisotropic heterogeneities. Initial results demonstrate that the heterogeneous region lowers the P/S amplitude ratios of an explosion. The presence of strong anisotropic heterogeneities near the source can also produce an S-wavefield with similar amplitudes to the P-wavefield. Lastly, we show that S-wave polarization angles vary as a function of azimuth and distance for explosive sources in the presence of small-scale heterogeneities.

58 GEOSCIENCES↗
NASA NTRS1971

Satellite Relay Telemetry of Seismic Data in Earthquake Prediction and Control

The Satellite Telemetry Earthquake Monitoring Program was started to evaluate the applicability of satellite relay telemetry in the collection of seismic data from a large number of dense seismograph clusters laid out along the major fault systems of western North America. Prototype clusters utilizing phone-line telemetry were then being installed by the National Center for Earthquake Research in 3 regions along the San Andreas fault in central California; and the experience of installing and operating the clusters and in reducing and analyzing the seismic data from them was to provide the raw materials for evaluation in the satellite relay telemetry project. The principal advantages of the satellite relay system over commercial telephone or microwave systems were: (1) it could be made less prone to massive failure during a major earthquake; (2) it could be extended readily into undeveloped regions; and (3) it could provide flexible, uniform communications over large sections of major global tectonic zones. Fundamental characteristics of a communications system to cope with the large volume of raw data collected by a short-period seismograph network are discussed.

Jackson, W. H.↗
NASA NTRS1975

A surface-wave investigation of the rupture mechanism of the Gobi-Altai (4 December 1957) earthquake

Long period records of multiple Love waves from the 1957 earthquake in Mongolia at Pasadena are analyzed and compared to synthetic seismograms, generated by the method of Kanamori. A fit in the time domain shows that the records are not consistent with the previous solution, achieved through a frequency domain analysis of directivity. The solution asks for a shorter rupture of 270 km at a velocity of 3.5 km/s. The focal parameters are constrained by updating all the reported first motion and are found to be: Strike = 103 deg, Dip = 53 deg, Slip = 32 deg. A seismic moment of 1.8 10 to the 28th power dynes-cm is obtained. These figures are also consistent with a time domain analysis of Love waves at Palisades and Strasbourg, and of Rayleigh waves at Pasadena, with a directivity study of Love waves at Pasadena, and with static deformation and isoseismal data. A discussion is given of the relation between moment, magnitude and rupture area, and a comparison is made with other events in the same region: It is concluded that this earthquake does not exhibit an intra-plate behavior, but rather compares better with inter-plate events, such as the great Assam earthquake.

Okal, E. A.↗
NASA NTRS1977

A local earthquake coda magnitude and its relation to duration, moment M sub O, and local Richter magnitude M sub L

A relationship was found between the seismic moment, M sub O, of shallow local earthquakes and the total duration of the signal, t, in seconds, measured from the earthquakes origin time, assuming that the end of the coda is composed of backscattering surface waves due to lateral heterogenity in the shallow crust following Aki. Using the linear relationship between the logarithm of M sub O and the local Richter magnitude M sub L, a relationship between M sub L and t, was found. This relationship was used to calculate a coda magnitude M sub C which was compared to M sub L for Southern California earthquakes which occurred during the period from 1972 to 1975.

Suteau, A. M.↗
NASA NTRS1978

Finite element modeling of stress in the Nazca plate - Driving forces and plate boundary earthquakes

The state of stress within the Nazca plate due to plate driving forces and large plate boundary earthquakes has been analyzed by applying a finite element method using the wave front solution technique to models of the intraplate stress field in a single plate using a refined grid. Although only static elastic models have been explicitly calculated, certain limiting cases of an elastic plate over a viscous asthenosphere were also treated. A state of nearly east-west compression inferred from the source mechanism of thrust earthquakes in the interior of the plate requires ridge pushing forces. The net pulling force on the oceanic plate by the subducted slab has a maximum value comparable to pushing forces. The estimated horizontal deviatoric stress in intraplate regions, based on potential forces associated with the ridge, is on the order of a few hundred bars. The intraplate stress field in the region of the 1960 earthquake may change by a few tens of bars at most once the asthenosphere has relaxed, with changes on the order of one bar occurring at greater distances into the plate. The changes in the intraplate stress field are probably not noticeable unless the lithosphere is near failure.

Richardson, R. M.↗
NASA NTRS1979

Potential utilization of the NASA/George C. Marshall Space Flight Center in earthquake engineering research

Earthquake engineering research capabilities of the National Aeronautics and Space Administration (NASA) facilities at George C. Marshall Space Flight Center (MSFC), Alabama, were evaluated. The results indicate that the NASA/MSFC facilities and supporting capabilities offer unique opportunities for conducting earthquake engineering research. Specific features that are particularly attractive for large scale static and dynamic testing of natural and man-made structures include the following: large physical dimensions of buildings and test bays; high loading capacity; wide range and large number of test equipment and instrumentation devices; multichannel data acquisition and processing systems; technical expertise for conducting large-scale static and dynamic testing; sophisticated techniques for systems dynamics analysis, simulation, and control; and capability for managing large-size and technologically complex programs. Potential uses of the facilities for near and long term test programs to supplement current earthquake research activities are suggested.

Scholl, R. E.↗
NASA NTRS1980

Regional analysis of earthquake occurrence and seismic energy release

The historic temporal variation in earthquake occurrence and seismic energy release on a regional basis throughtout the world were studied. The regionalization scheme employed divided the world into large areas based either on seismic and tectonic considerations (Flinn-Engdahl Scheme) or geographic (longitude and latitude) criteria. The data set is the wide earthquake catalog of the National Geophysical Solar-Terrestrial Data Center. An apparent relationship exists between the maximum energy released in a limited time within a seismic region and the average or background energy per year averaged over a long time period. In terms of average or peak energy release, the most seismic regions of the world during the 50 to 81 year period ending in 1977 were Japanese, Andean South American, and the Alaska-Aleutian Arc regions. The year to year fluctuations in regional seismic energy release are greater, by orders of magnitude, than the corresponding variations in the world-wide seismic energy release. The b values of seismic regions range from 0.7 to 1.4 where earthquake magnitude is in the range 6.0 to 7.5.

Cohen, S. C.↗
NASA NTRS1981

Strain buildup and release, earthquake prediction and selection of VBL sites for margins of the north Pacific

During the past year, the grant supported research on several aspects of crustal deformation. The relation between earthquake displacements and fault dimensions was studied in an effort to find scaling laws that relate static parameters such as slip and stress drop to the dimensions of the rupture. Several implications of the static relations for the dynamic properties of earthquakes such as rupture velocity and dynamic stress drop were proposed. A theoretical basis for earthquake related phenomena associated with slow rupture growth or propagation, such as delayed multiple events, was developed using the stress intensity factor defined in fracture mechanics and experimental evidence from studies of crack growth by stress corrosion. Finally, extensive studies by Japanese geologists have established the offset across numerous faults in Japan over the last one hundred thousand years. These observations of intraplate faulting are being used to establish the spatial variations of the average strain rate of subregions in southern Japan.

Scholz, C. H.↗