Engineering PapersSearch

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 55 records · Page 3

Learning earthquake ground motions via conditional generative modeling

Predicting high-fidelity ground motions for future earthquakes is crucial for seismic hazard assessment and infrastructure resilience. Conventional empirical simulations suffer from sparse sensor distribution and geographically localized earthquake locations, while physics-based methods are computationally intensive and require accurate representations of Earth structures and earthquake sources. We propose an artificial intelligence (AI) spectrogram generator, Conditional Generative Modeling for Ground Motion (CGM-GM). CGM-GM leverages earthquake magnitudes and geographic coordinates of earthquakes and sensors as inputs, when postprocessed with phase information, capturing spatially continuous Fourier amplitude spectra (FAS) as well as properties such as P and S arrivals, and waveform durations, without explicit physics constraints. This is achieved through a probabilistic autoencoder that extracts latent distributions in the time-frequency domain and variational sequential models for prior and posterior distributions. We evaluate the performance of CGM-GM using small-magnitude earthquake records from the San Francisco Bay Area, a region with high seismic risks. Here, we report that CGM-GM demonstrates potential for complementing physics-based simulations and non-ergodic empirical ground motion models, as well as shows promise in seismology and beyond.

geophysics

Increased flood exposure in the Pacific Northwest following earthquake-driven subsidence and sea-level rise

Climate-driven sea-level rise is increasing the frequency of coastal flooding worldwide, exacerbated locally by factors like land subsidence from groundwater and resource extraction. However, a process rarely considered in future sea-level rise scenarios is sudden (over minutes) land subsidence associated with great (>M8) earthquakes, which can exceed 1 m. Along the Washington, Oregon, and northern California coasts, the next great Cascadia subduction zone earthquake could cause up to 2 m of sudden coastal subsidence, dramatically raising sea level, expanding floodplains, and increasing the flood risk to local communities. Here, we quantify the potential expansion of the 1% floodplain (i.e., the area with an annual flood risk of 1%) under low (~0.5 m), medium (~1 m), and high (~2 m) earthquake-driven subsidence scenarios at 24 Cascadia estuaries. If a great earthquake occurred today, floodplains could expand by 90 km 2 (low), 160 km 2 (medium), or 300 km 2 (high subsidence), more than doubling the flooding exposure of residents, structures, and roads under the high subsidence scenario. By 2100, when climate-driven sea-level rise will compound the hazard, a great earthquake could expand floodplains by 170 km 2 (low), 240 km 2 (medium), or 370 km 2 (high subsidence), more than tripling the flooding exposure of residents, structures, and roads under the high subsidence scenario compared to the 2023 floodplain. Our findings can support decision-makers and coastal communities along the Cascadia subduction zone as they prepare for compound hazards from the earthquake cycle and climate-driven sea-level rise and provide critical insights for tectonically active coastlines globally.

Science & Technology - Other Topics

Transform fault earthquakes in the North Atlantic: Source mechanisms and depth of faulting

The centroid depths and source mechanisms of 12 large earthquakes on transform faults of the northern Mid-Atlantic Ridge were determined from an inversion of long-period body waveforms. The earthquakes occurred on the Gibbs, Oceanographer, Hayes, Kane, 15 deg 20 min, and Vema transforms. The depth extent of faulting during each earthquake was estimated from the centroid depth and the fault width. The source mechanisms for all events in this study display the strike slip motion expected for transform fault earthquakes; slip vector azimuths agree to 2 to 3 deg of the local strike of the zone of active faulting. The only anomalies in mechanism were for two earthquakes near the western end of the Vema transform which occurred on significantly nonvertical fault planes. Secondary faulting, occurring either precursory to or near the end of the main episode of strike-slip rupture, was observed for 5 of the 12 earthquakes. For three events the secondary faulting was characterized by reverse motion on fault planes striking oblique to the trend of the transform. In all three cases, the site of secondary reverse faulting is near a compression jog in the current trace of the active transform fault zone. No evidence was found to support the conclusions of Engeln, Wiens, and Stein that oceanic transform faults in general are either hotter than expected from current thermal models or weaker than normal oceanic lithosphere.

Bergman, Eric A.

Transform fault earthquakes in the North Atlantic - Source mechanisms and depth of faulting

The centroid depths and source mechanisms of 12 large earthquakes on transform faults of the northern Mid-Atlantic Ridge were determined from an inversion of long-period body waveforms. The earthquakes occurred on the Gibbs, Oceanographer, Hayes, Kane, 15 deg 20 min, and Vema transforms. The depth extent of faulting during each earthquake was estimated from the centroid depth and the fault width. The source mechanisms for all events in this study display the strike slip motion expected for transform fault earthquakes; slip vector azimuths agree to 2 to 3 deg of the local strike of the zone of active faulting. The only anomalies in mechanism were for two earthquakes near the western end of the Vema transform which occurred on significantly nonvertical fault planes. Secondary faulting, occurring either precursory to or near the end of the main episode of strike-slip rupture, was observed for 5 of the 12 earthquakes. For three events the secondary faulting was characterized by reverse motion on fault planes striking oblique to the trend of the transform. In all three cases, the site of secondary reverse faulting is near a compression jog in the current trace of the active transform fault zone. No evidence was found to support the conclusions of Engeln, Wiens, and Stein that oceanic transform faults in general are either hotter than expected from current thermal models or weaker than normal oceanic lithosphere.

Bergman, Eric A.

Earthquake swarms on the Mid-Atlantic Ridge - Products of magmatism or extensional tectonics?

The spatial and temporal patterns and other characteristics of earthquakes in 34 earthquake swarms on the Mid-Atlantic Ridge were compared with those of well-studied earthquake swarms which accompany terrestrial volcanic eruptions, to test the assumption that the teleseismically observed earthquake swarms along mid-ocean ridges are indicators of volcanism. Improved resolution of these patterns for the mid-ocean ridge events was achieved by a multiple-event relocation technique. It was found that the teleseismically located earthquake swarms on the mid-ocean ridge system have few features in common with swarms directly associated with active magmatism in terrestrial volcanic rift zones such as Hawaii and Iceland. While the possibility that some of the mid-ocean earthquake swarms might be directly associated with a current episode of eruptive activity on the Mid-Atlantic Ridge cannot be excluded, none of the 34 swarms studied in this work was found to be a conspicuously attractive candidate for such a role.

Bergman, Eric A.

Earthquake recurrence on the southern San Andreas modulated by fault-normal stress

Earthquake recurrence data from the Pallett Creek and Wrightwood paleoseismic sites on the San Andreas fault appear to show temporal variations in repeat interval. We investigate the interaction between strike-slip faults and auxiliary reverse and normal faults as a physical mechanism capable of producing such variations. Under the assumption that fault strength is a function of fault-normal stress (e.g. Byerlee's Law), failure of an auxiliary fault modifies the strength of the strike-slip fault, thereby modulating the recurrence interval for earthquakes. In our finite element model, auxiliary faults are driven by stress accumulation near restraining and releasing bends of a strike-slip fault. Earthquakes occur when fault strength is exceeded and are incorporated as a stress drop which is dependent on fault-normal stress. The model is driven by a velocity boundary condition over many earthquake cycles. Resulting synthetic strike-slip earthquake recurrence data display temporal variations similar to observed paleoseismic data within time windows surrounding auxiliary fault failures. Our simple model supports the idea that interaction between a strike-slip fault and auxiliary reverse or normal faults can modulate the recurrence interval of events on the strike-slip fault, possibly producing short term variations in earthquake recurrence interval.

Palmer, Randy

Co-Seismic Energy Changes Induced by Earthquakes on a Rotating, Gravitating Earth

Besides operating its own energy budget, an earthquake acts as an agent transferring a much greater amount of energy among the Earth's rotation, elastic field, gravitational field and internal heat. We compute the co-seismic, globally integrated gravitational and rotation changes induced by some 20,000 large earthquakes that occurred in the last quarter century, according to Chao et al. (1995, GJI, 122,776- 783,784-789) and using the Harvard CMT catalog. The result confirms an extremely strong tendency for the earthquakes to decrease the global gravitational energy and to increase the spin energy. It is found that energy is being extracted from the Earth's gravitational field by the action of earthquakes at an average rate of about approx. 2 TeraW during the studied period, larger by far than the approx. 7 GigaW for the average rate of the earthquake-induced rotational energy increase and the approx. 5 GigaW for the seismic energy release. Based on energetics considerations and assuming the inability of the Earth to build up elastic energy continuously over time, it is argued that earthquakes, by converting gravitational energy, may make a significant contribution to the global hedflow.

Chao, Benjamin F.

Atmosphere-Ionosphere Response to the M9 Tohoku Earthquake Revealed by Joined Satellite and Ground Observations. Preliminary Results

The recent M9 Tohoku Japan earthquake of March 11, 2011 was the largest recorded earthquake ever to hit this nation. We retrospectively analyzed the temporal and spatial variations of four different physical parameters - outgoing long wave radiation (OLR), GPS/TEC, Low-Earth orbit tomography and critical frequency foF2. These changes characterize the state of the atmosphere and ionosphere several days before the onset of this earthquake. Our first results show that on March 8th a rapid increase of emitted infrared radiation was observed from the satellite data and an anomaly developed near the epicenter. The GPS/TEC data indicate an increase and variation in electron density reaching a maximum value on March 8. Starting on this day in the lower ionospheric there was also confirmed an abnormal TEC variation over the epicenter. From March 3-11 a large increase in electron concentration was recorded at all four Japanese ground based ionosondes, which return to normal after the main earthquake. We found a positive correlation between the atmospheric and ionospheric anomalies and the Tohoku earthquake. This study may lead to a better understanding of the response of the atmosphere/ionosphere to the Great Tohoku earthquake.

Ouzounov, Dimitar

The Effect of Sonic Booms on Earthquake Warning Systems

Several aerospace companies are designing quiet supersonic business jets for service over the United States. These aircraft have the potential to increase the occurrence of mild sonic booms across the country. This leads to interest among earthquake warning (EQW) developers and the general seismological community in characterizing the effect of sonic booms on seismic sensors in the field, their potential impact on EQW systems, and means of discriminating their signatures from those of earthquakes. The SonicBREWS project (Sonic Boom Resistant Earthquake Warning Systems) is a collaborative effort between Seismic Warning Systems, Inc. (SWS) and NASA Dryden Flight Research Center. This project aims to evaluate the effects of sonic booms on EQW sensors. The study consists of exposing high-sample-rate (1000 sps) triaxial accelerometers to sonic booms with overpressures ranging from 10 to 600 Pa in the free field and the built environment. The accelerometers record the coupling of the sonic boom to the ground and surrounding structures, while microphones record the acoustic wave above ground near the sensor. Sonic booms are broadband signals with more high-frequency content than earthquakes. Even a 1000 sps accelerometer will produce a significantly aliased record. Thus the observed peak ground velocity is strongly dependent on the sampling rate, and increases as the sampling rate is reduced. At 1000 sps we observe ground velocities that exceed those of P-waves from M L 3 earthquakes at local distances, suggesting that sonic booms are not negligible for EQW applications. We present the results of several experiments conducted under SonicBREWS showing the effects of typical-case low amplitude sonic booms and worst-case high amplitude booms. We show the effects of various sensor placements and sensor array geometries. Finally, we suggest possible avenues for discriminating sonic booms from earthquakes for the purposes of EQW.

Wurman, Gilead

Multi-Sensor Observations of Earthquake Related Atmospheric Signals over Major Geohazard Validation Sites

We are conducting a scientific validation study involving multi-sensor observations in our investigation of phenomena preceding major earthquakes. Our approach is based on a systematic analysis of several atmospheric and environmental parameters, which we found, are associated with the earthquakes, namely: thermal infrared radiation, outgoing long-wavelength radiation, ionospheric electron density, and atmospheric temperature and humidity. For first time we applied this approach to selected GEOSS sites prone to earthquakes or volcanoes. This provides a new opportunity to cross validate our results with the dense networks of in-situ and space measurements. We investigated two different seismic aspects, first the sites with recent large earthquakes, viz.- Tohoku-oki (M9, 2011, Japan) and Emilia region (M5.9, 2012,N. Italy). Our retrospective analysis of satellite data has shown the presence of anomalies in the atmosphere. Second, we did a retrospective analysis to check the re-occurrence of similar anomalous behavior in atmosphere/ionosphere over three regions with distinct geological settings and high seismicity: Taiwan, Japan and Kamchatka, which include 40 major earthquakes (M>5.9) for the period of 2005-2009. We found anomalous behavior before all of these events with no false negatives; false positives were less then 10%. Our initial results suggest that multi-instrument space-borne and ground observations show a systematic appearance of atmospheric anomalies near the epicentral area that could be explained by a coupling between the observed physical parameters and earthquake preparation processes.

Ouzounov, D.

Transient Effects in Atmosphere and Ionosphere Preceding the 2015 M7.8 and M7.3 Gorkha–Nepal Earthquakes

We analyze retrospectively/prospectively the transient variations of six different physical parameters in the atmosphere/ionosphere during the M7.8 and M7.3 earthquakes in Nepal, namely: 1) outgoing longwave radiation (OLR) at the top of the atmosphere (TOA); 2) GPS/TEC; 3) the very-low-frequency (VLF/LF) signals at the receiving stations in Bishkek (Kyrgyzstan) and Varanasi (India); 4) Radon observations; 5) Atmospheric chemical potential from assimilation models; and; 6) Air Temperature from NOAA ground stations. We found that in mid-March 2015, there was a rapid increase in the radiation from the atmosphere observed by satellites. This anomaly was located close to the future M7.8 epicenter and reached a maximum on April 21–22. The GPS/TEC data analysis indicated an increase and variation in electron density, reaching a maximum value during April 22–24. A strong negative TEC anomaly in the crest of EIA (Equatorial Ionospheric Anomaly) occurred on April 21, and a strong positive anomaly was recorded on April 24, 2015. The behavior of VLF-LF waves along NWC-Bishkek and JJY-Varanasi paths has shown abnormal behavior during April 21–23, several days before the first, stronger earthquake. Our continuous satellite OLR analysis revealed this new strong anomaly on May 3, which was why we anticipated another major event in the area. On May 12, 2015, an M7.3 earthquake occurred. Our results show coherence between the appearance of these pre-earthquake transient’s effects in the atmosphere and ionosphere (with a short time-lag, from hours up to a few days) and the occurrence of the 2015 M7.8 and M7.3 events. The spatial characteristics of the pre-earthquake anomalies were associated with a large area but inside the preparation region estimated by Dobrovolsky-Bowman. The preearthquake nature of the signals in the atmosphere and ionosphere was revealed by simultaneous analysis of satellite, GPS/TEC, and VLF/LF and suggest that they follow a general temporal-spatial evolution pattern that has been seen in other large earthquakes worldwide

Dimitar Ouzounov

Time Series Foundation Models and Deep Learning Architectures for Earthquake Temporal and Spatial Nowcasting

Advancing the capabilities of earthquake nowcasting, the real-time forecasting of seismic activities, remains crucial for reducing casualties. This multifaceted challenge has recently gained attention within the deep learning domain, facilitated by the availability of extensive earthquake datasets. Despite significant advancements, the existing literature on earthquake nowcasting lacks comprehensive evaluations of pre-trained foundation models and modern deep learning architectures; each focuses on a different aspect of data, such as spatial relationships, temporal patterns, and multi-scale dependencies. This paper addresses the mentioned gap by analyzing different architectures and introducing two innovative approaches called Multi Foundation Quake and GNNCoder. We formulate earthquake nowcasting as a time series forecasting problem for the next 14 days within 0.1-degree spatial bins in Southern California. Earthquake time series are generated using the logarithm energy released by quakes, spanning 1986 to 2024. Our comprehensive evaluations demonstrate that our introduced models outperform other custom architectures by effectively capturing temporal-spatial relationships inherent in seismic data. The performance of existing foundation models varies significantly based on the pre-training datasets, emphasizing the need for careful dataset selection. However, we introduce a novel method, Multi Foundation Quake, that achieves the best overall performance by combining a bespoke pattern with Foundation model results handled as auxiliary streams.

97 MATHEMATICS AND COMPUTING

How earthquakes organize stress

Stress is not uniform in the Earth. Therefore, we must use natural experiments to measure the distribution of stresses and related quantities, rather than single values. For instance, dynamic triggering shows that faults are uniformly distributed over their loading cycles in Southern California. The probability that a fault ruptures across a barrier measures the in situ energy distribution. Fault roughness reflects the distribution of strength. These natural experiments produce observable distributions that are surprisingly consistent and suggest some degree of self-organization in the Earth’s crust. Once established, the functional form of the distributions can be used to track changes in response to earthquakes as well as to distinguish fundamentally different fault systems. Transient fault locking before stress release in laboratory experiments can be interpreted as a consequence of self-organization of fault stress. The robust self-organization of multiple variables in earthquake systems suggests that the most consequential mechanical outcome of earthquakes may be the redistribution of stress and the strain energy associated with it. The low friction on a fault during seismic slip as inferred by temperature measurements of the Tohoku earthquake is consistent with dissipation playing a secondary role to this redistribution process. Through stress redistribution and interaction, subduction zone faults tend to synchronize, perhaps due to their geometric simplicity, while the continental system of Southern California cannot synchronize, perhaps due to the complexity of the fault network. Earthquakes organize stress in the crust and produce a suite of well-defined, consistent distributions.

earthquakes

Earthquake detection in a simulated lunar regolith using distributed acoustic sensing

Current models of inner lunar geology have largely been inferred from the seismic experiments and observations performed during the Apollo missions that comprised a relatively small number of seismic instruments. Refining constraints on fundamental lunar relationships such as crust-mantle and mantle-core boundaries in the future will require seismic arrays spanning larger epicentral distances. A promising technology for installing dense seismic arrays with minimal human effort is distributed acoustic sensing (DAS), an approach that allows a single length of fiber optic cable to act as hundreds or thousands of sensors when coupled with a DAS interrogator. While terrestrial uses of DAS technology for seismic monitoring rely on burying the cable to maximize fidelity of seismic signal transmission to the fiber, digging meters of trench to bury optical fiber on lunar or planetary surfaces is logistically infeasible. To evaluate DAS signal attenuation due to surface deployment of cable in lunar regolith, we completed earthquake detection analyses that evaluated the sensitivity of an optic-fiber DAS system to seismic signals at different burial depths. We deployed a single-mode fiber in a 10-m open-bottom wooden box filled with a lunar regolith simulant (LRS) with fiber buried at different depths within the LRS and recorded signals for four regional and local earthquakes. The results were used to identify and evaluate signal attenuation in surface-deployed fiber compared to buried fiber in the LRS. Burial depth responses to active-source signals were also evaluated similar to previous studies characterizing DAS sensitivity of surface-deployed fiber. Atmospheric noise was minimal as the cable was deployed in an indoor environment; however, where observed, atmospheric and anthropogenic noise was filtered out using the same bandpass filtering used to identify earthquake events. We found that signal attenuation of the surface-deployed fiber compared to buried fiber was relatively high in active-source experiments but was not consistently observed in earthquake signals. That burial depth is not highly correlated to attenuation of the observed earthquake signals indicates that in a noise-limited environment, placing DAS-interrogated fiber directly at the regolith surface may be a promising deployment strategy to consider for sensing remote seismic signals during lunar exploration.

58 GEOSCIENCES

Constructing a High‐Resolution Aftershock Catalog for the 2017 Mw 8.2 Tehuantepec Earthquake Sequence Using a Machine Learning–Based Workflow

The 8 September 2017 Mw 8.2 Tehuantepec earthquake was the largest instrumentally recorded normal‐faulting earthquake in Mexico. The mainshock occurred offshore within the Tehuantepec seismic gap, generating >30,000 aftershocks in the following year. We applied an open‐source, machine learning (ML)–assisted workflow to construct a high‐resolution aftershock catalog using data from temporary and permanent seismic networks in southern Mexico. The workflow integrates PhaseNet for phase detection; GaMMA for phase association; and VELEST, HypoInverse, and HypoDD for velocity modeling and relocation. We processed seven months of continuous waveform data from 29 broadband stations, including a temporary rapid‐response deployment that improved station coverage of the offshore rupture zone. To evaluate performance, we compared our results against analyst‐reviewed picks and event locations from the Servicio Sismológico Nacional catalog. The resulting catalog contains 11,374 relocated earthquakes and represents the most comprehensive published dataset for this sequence, incorporating the first full use of the temporary network. Relocated hypocenters show improved depth control and align well with the Slab2.0 subduction geometry, revealing clearer separation between offshore slab events and onshore crustal seismicity. This study demonstrates that combining ML‐based detection with established methods provides a scalable and reproducible approach for constructing high‐quality earthquake catalogs in tectonically complex environments and offers practical guidance for adapting similar workflows to other earthquake sequences.

Garcia, Marc [The University of Texas at El Paso,

The rupture process and asperity distribution of three great earthquakes from long-period diffracted P-waves

The variation of maximum earthquake size along the subduction zones has been interpreted as a variation in the seismic coupling ostensibly related to the mechanical conditions of the fault zone. Great differences are noted between the seismographs of the three great earthquakes whose rupture processes are presently considered: in the Kurile Islands (1963), The Rat Islands (1965) and Alaska (1964). On-scale long period P waves were recorded in all cases. Source time functions are deconvolved from the observed periods. It is concluded that maximum earthquake size is related to the asperity distribution on the fault. The subduction zones with the largest earthquakes have very large asperities, as in the Alaskan case, while the zones with the smaller great earthquakes, such as the Kurile Islands, have smaller scattered asperities.

Ruff, L.

Scaling differences between large interplate and intraplate earthquakes

A study of large intraplate earthquakes with well determined source parameters shows that these earthquakes obey a scaling law similar to large interplate earthquakes, in which M sub o varies as L sup 2 or u = alpha L where L is rupture length and u is slip. In contrast to interplate earthquakes, for which alpha approximately equals 1 x .00001, for the intraplate events alpha approximately equals 6 x .0001, which implies that these earthquakes have stress-drops about 6 times higher than interplate events. This result is independent of focal mechanism type. This implies that intraplate faults have a higher frictional strength than plate boundaries, and hence, that faults are velocity or slip weakening in their behavior. This factor may be important in producing the concentrated deformation that creates and maintains plate boundaries.

Scholz, C. H.

Intraplate earthquakes and the state of stress in oceanic lithosphere

The dominant sources of stress relieved in oceanic intraplate earthquakes are investigated to examine the usefulness of earthquakes as indicators of stress orientation. The primary data for this investigation are the detailed source studies of 58 of the largest of these events, performed with a body-waveform inversion technique of Nabelek (1984). The relationship between the earthquakes and the intraplate stress fields was investigated by studying, the rate of seismic moment release as a function of age, the source mechanisms and tectonic associations of larger events, and the depth-dependence of various source parameters. The results indicate that the earthquake focal mechanisms are empirically reliable indicators of stress, probably reflecting the fact that an earthquake will occur most readily on a fault plane oriented in such a way that the resolved shear stress is maximized while the normal stress across the fault, is minimized.

Bergman, Eric A.