Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Earthquakes”

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

Fault Network Geometry Modulates Earthquake Source Spectra Across Scales

Earthquake source spectra provide unique insights into the earthquake rupture process. Motivated by previous research suggesting that complex fault geometries enhance high‐frequency seismic radiation, we study the influence of fault network geometry on earthquake source spectra using multiple independent observations. At regional scales, we examine correlations of stress drop measurements with surface fault trace misalignment in Southern California, Japan, and Central Italy. At a global scale, we examine correlations of moment‐rate function complexity of large earthquakes with focal mechanism variability, a proxy for local fault complexity. Despite significant scatter in the observations, we find overall consistent positive correlations. The concept that elastic interactions of discrete fault structures during the earthquake rupture process generates high‐frequency ground motions offers a coherent framework for interpreting our observations. These findings suggest that variations in fault complexity explain why some earthquakes produce stronger high‐frequency ground motions than others.

Lee, Jaeseok [Brown Univ., Providence, RI (United ↗
DOE OSTI2026

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

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

Incorporating Full Elastodynamic Effects and Dipping Fault Geometries in Community Code Verification Exercises for Simulations of Earthquake Sequences and Aseismic Slip (SEAS)

ABSTRACT Numerical modeling of earthquake dynamics and derived insight for seismic hazard relies on credible, reproducible model results. The sequences of earthquakes and aseismic slip (SEAS) initiative has set out to facilitate community code comparisons, and verify and advance the next generation of physics-based earthquake models that reproduce all phases of the seismic cycle. With the goal of advancing SEAS models to robustly incorporate physical and geometrical complexities, here we present code comparison results from two new benchmark problems: BP1-FD considers full elastodynamic effects, and BP3-QD considers dipping fault geometries. Seven and eight modeling groups participated in BP1-FD and BP3-QD, respectively, allowing us to explore these physical ingredients across multiple codes and better understand associated numerical considerations. With new comparison metrics, we find that numerical resolution and computational domain size are critical parameters to obtain matching results. Codes for BP1-FD implement different criteria for switching between quasi-static and dynamic solvers, which require tuning to obtain matching results. In BP3-QD, proper remote boundary conditions consistent with specified rigid body translation are required to obtain matching surface displacements. With these numerical and mathematical issues resolved, we obtain excellent quantitative agreements among codes in earthquake interevent times, event moments, and coseismic slip, with reasonable agreements made in peak slip rates and rupture arrival time. We find that including full inertial effects generates events with larger slip rates and rupture speeds compared to the quasi-dynamic counterpart. For BP3-QD, both dip angle and sense of motion (thrust versus normal faulting) alter ground motion on the hanging and foot walls, and influence event patterns, with some sequences exhibiting similar-size characteristic earthquakes, and others exhibiting different-size events. These findings underscore the importance of considering full elastodynamics and nonvertical dip angles in SEAS models, as both influence short- and long-term earthquake behavior and are relevant to seismic hazard.

Geochemistry & Geophysics↗
DOE OSTI2024

Inferring the Focal Depths of Small Earthquakes in Southern California Using Physics-Based Waveform Features

Determining the depths of small crustal earthquakes is challenging in many regions of the world, because most seismic networks are too sparse to resolve trade-offs between depth and origin time with conventional arrival-time methods. Precise and accurate depth estimation is important, because it can help seismologists discriminate between earthquakes and explosions, which is relevant to monitoring nuclear test ban treaties and producing earthquake catalogs that are uncontaminated by mining blasts. Here, we examine the depth sensitivity of several physics-based waveform features for ~8000 earthquakes in southern California that have well-resolved depths from arrival-time inversion. We focus on small earthquakes (2 < M L < 4) recorded at local distances (<150 km), for which depth estimation is especially challenging. We find that differential magnitudes (M w /M L –M c ) are positively correlated with focal depth, implying that coda wave excitation decreases with focal depth. We analyze a simple proxy for relative frequency content, Φ≡log 10 (M 0 )+3log 10 (f c ), and find that source spectra are preferentially enriched in high frequencies, or “blue-shifted,” as focal depth increases. Here, we also find that two spectral amplitude ratios Rg 0.5–2 Hz/Sg 0.5–8 Hz and Pg/Sg at 3–8 Hz decrease as focal depth increases. Using multilinear regression with these features as predictor variables, we develop models that can explain 11%–59% of the variance in depths within 10 subregions and 25% of the depth variance across southern California as a whole. We suggest that incorporating these features into a machine learning workflow could help resolve focal depths in regions that are poorly instrumented and lack large databases of well-located events. Some of the waveform features we evaluate in this study have previously been used as source discriminants, and our results imply that their effectiveness in discrimination is partially because explosions generally occur at shallower depths than earthquakes.

58 GEOSCIENCES↗
DOE OSTI2024

Proximal Observations of Epicentral Infrasound Generated by Shallow Low-Magnitude Earthquakes in the Permian Basin, West Texas

Infrasound generated by earthquakes and explosions is generally detected at receivers at epicentral distances of hundreds to thousands of kilometers. However, proximal (<50 km) observations are especially important for low-magnitude earthquakes and low-yield explosions that may not generate signals capable of being detected at great ranges. Here, in this study, we present on the signals detected on an infrasound array 3 km away from two M L 2.9 earthquakes in the Permian Basin of west Texas. Local infrasound (LIS), generated at receivers during the passage of seismic waves, was detected following each earthquake. Epicentral infrasound (EIS), created at or near the epicenter and propagating away as a sound wave, was also detected. Array processing methods show that the EIS signals arrive from the same direction as the earthquake epicenters and at acoustic speeds. To our knowledge, these are the first observations of laterally propagating EIS at proximal ranges following an earthquake of any magnitude.

58 GEOSCIENCES↗
DOE OSTI2023

Ground Motion Models (GMMs) Improvements Using Earthquake Simulations on High Performance Computers

A computationally efficient simulation platform was developed that can provide representative synthetic ground motions from crustal earthquakes in the Stable Continental Regions of Central and Eastern US (CEUS), using 3D modeling and high-performance computing. The main objective was to use synthetic ground motion to provide constrains to refinements of exiting ergodic Ground Motion Models (GMMs), for large magnitude earthquakes and near-fault distances, for which these models are less reliable. Physics-based broadband (0-5Hz) ground motion simulations were used to estimate the near-fault ground motion amplitudes and within event and between-event variabilities associated with fault rupture characteristics. As part of a strategy for selecting a reginal velocity model and validation of developed rupture modeling technique, ground motions from the moment magnitude Mw5.0 November 7, 2016, Cushing Oklahoma, and Mw5.8 September 3, 2016, Pawnee Oklahoma earthquakes were simulated. In our simulations we used a 3D regional velocity model that was based on Saikia’s 1D velocity model. Saikia’s model demonstrated better performance in modelling high frequency regional wave propagation for CEUS region. The proposed 3D model includes lateral variations added to the 1D background model using the stochastic scheme of Pitarka and Mellors. Comparisons of the simulations with recordings of both earthquakes demonstrated the reliability of our deterministic simulation approach while emphasizing the importance of including small-scale variability in the regional velocity model needed to reproduce the observed high-frequency wave scattering effects. As part of validation analysis, comparisons with different GMMs for a Mw6.5 earthquake in the CESUS region resulted in a very good match between the simulated and empirical ground motion models. Initial investigations of within-event and between-event ground motion variabilities for Mw6.5 scenario earthquakes on a strike-slip fault, suggest that they are strongly related to spatial slip and slip rate variations, average rupture velocity, rupture area and rupture initiation location. For certain scenarios we found that the ground motion variability observed at near-fault distances (< 5 km) also persists at longer distances. Regardless of the rupture scenario, the simulated ground motion tends to fully saturate at short distances and for all periods. The near-fault saturation has to do with the attenuation of waves propagating along the fault and local rupture radiation pattern that also contribute to stronger ground motion variation at such distances. Analysis of effects of rupture initiation location suggest that the peak ground motion (PGV) and spectral acceleration (SA) can be quite variable due to rupture directivity effects. Such effects are stronger at periods longer than 1s. The effect of the 1D velocity models and surface topography on simulated ground motion were investigated by comparing three component synthetic seismograms computed at selected sites. Effect of surface topography was considered using the ratio between spectral accelerations simulated for two 1D models with flat surface topography and realistic model with surface topography. Overall, the topography slightly amplifies (by ~30%) the ground motion amplitude in the frequency range 1-3Hz. The effect of topography is more visible in the surface and coda waves portion of the seismograms.

58 GEOSCIENCES↗
NASA NTRS1987

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.↗
NASA NTRS1988

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.↗
NASA NTRS1990

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.↗
NASA NTRS1995

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↗
NASA NTRS2003

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.↗
NASA NTRS2011

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↗
NASA NTRS2011

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↗
NASA NTRS2012

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.↗
NASA NTRS2021

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

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

Periodicity Analysis of Earthquake Occurrence and Hypocenter Depth Near Parkfield, California, 1994–2002 Versus 2006–2014

In central California, periodic earthquake occurrence suggests a relationship with annual hydrological, atmospheric, thermal, and tidal loadings. Here, in catalogs of declustered earthquakes within 100 km from Parkfield, CA, we study a semiannual periodicity for the monthly number of ≤7.2-km-deep earthquakes and the monthly median hypocenter depth over 1994–2002 and 2006–2014. Peak-trough months in fitted periodic components differ between time spans and event populations. Deeper earthquakes present no semiannual and a weak annual periodicity. Although modeled pore pressure shows a Spring peak, when added to the elastic Coulomb stress from surface hydrospheric loads, it fails to predict a 6-month periodicity for 2006–2014. In 1994–2002, the pore-pressure amplitude appears to be of same order as the elastic stress and may have had a stronger effect. In 2006–2014, we did not find load model parameters explaining the observed pattern or apparent changes following the 2003 San Simeon and 2004 Parkfield earthquakes.

58 GEOSCIENCES↗