Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “earthquake location models”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Geodetic Measurements of Slow-Slip Events Southeast of Parkfield, CA

Tremor and low-frequency earthquakes are presumed to be indicative of surrounding slow, aseismic slip that is often below geodetic detection thresholds. This study uses data from borehole seismometers and long-baseline laser strainmeters to observe both the seismic and geodetic signatures of episodic tremor and slip on the Parkfield region of the San Andreas Fault near Cholame, CA. We note the observed occurrence rates of both the tremors and co-located families of low-frequency earthquakes are not steady but instead exhibit quasiperiodic bursts of increased activity.We show that these periods of elevated seismic activity correlate with statistically significant stacked strain signals consisting of 44 slow-slip events. Modeled individual slow-slip events and their total summed moment, which are constrained by seismic signals and stacked strain, respectively, indicate that the individual moment magnitudes of these events range from $M_w$ 4.6–5.2.We find that the measured geodetic signal likely precedes the seismic signal by several hours, consistent with the aseismic slip preceding and driving the observed seismic tremor activity.We confirm that strike-slip faults, in addition to subduction zones, are capable of producing episodic tremor and slip.

58 GEOSCIENCES↗
NASA NTRS1987

Earthquake source mechanisms and transform fault tectonics in the Gulf of California

The source parameters of 19 large earthquakes in the Gulf of California were determined from inversions of long-period P and SH waveforms. The goal was to understand the recent slip history of this dominantly transform boundary between the Pacific and North American plates as well as the effect on earthquake characteristics of the transition from young oceanic to continental lithosphere. For the better recorded transform events, the fault strike is resolved to + or - 4 deg at 90 percent confidence. The slip vectors thus provide important constraints on the direction of relative plate motion. Most centroid depths are poorly resolved because of tradeoffs between depth and source time function. On the basis of waveform modeling, historical seismicity, and other factors, it is appropriate to divide the Gulf into three distinct zones. The difference in seismic character among the three zones is likely the result of differing levels of maturity of the processes of rifting, generation of oceanic crust, and formation of stable oceanic transform faults. The mechanism of an earthquake on the Tres Marias Escarpment is characterized by thrust faulting and likely indicates the direction of relative motion between the Rivera and North American plates. This mechanism requires revision in plate velocity models which predict strike slip motion at this location.

Goff, John A.↗
NASA NTRS2021

Newmark-Type Pseudo-Three-Dimensional Back-Analysis of Co-Seismic Landslides in Egkremnoi, Lefkada, Greece

A pseudo-three-dimensional (pseudo-3D) methodology to back-analyze coseismic landslides was developed and applied to 68mapped landslides, which occurred over approximately 1 km2area in Egkremnoi, Greece, during the 2015 Mw6.5 Lefkada earthquake. The methodology is based on a one-dimensional (1D) Newmark-type sliding block model to assess instability and a spatial projection in 3D topography to derive landslide geometry. The strength parameters for modeled landslides that best match the landslide location, area, and volume were derived through an iterative scheme that optimizes the match using predefined matching criteria. The range of different-sized landslides produced shear strength estimates from ~10 to 300 kPa and led to the derivation of a regionally averaged strength envelope characterized by a cohesion of 6 kPa and a friction angle of 53°for the highly fractured limestones that are encountered in this area. Compared to previous full 3Dslope stability analyses in this area, the friction angle using this methodology was found to be generally consistent, but the cohesion was lower. The presented methodology can provide a computationally efficient method to estimate the average shear strength of a geologic unit over large areas, especially where extensive field and laboratory tests on the materials are unavailable or difficult to conduct.

Weibing Gong↗
DOE OSTI2025

Prediction and Analysis of Utah FORGE Injection Activities using a Coupled Thermo-hydro-mechanical and Earthquake (THM+E) Modeling Workflow

A coupled thermo-hydro-mechanical (THM) numerical workflow that is capable of modeling seismic slip is critical for the successful development of enhanced geothermal systems (EGS). By integrating key physical processes, this workflow enables accurate simulation of temperature and pressure diffusions, stress changes, and induced seismicity. As a result, it serves as a vital tool for predicting induced seismicity and optimizing reservoir stimulation strategies. The Utah FORGE (Frontier Observatory for Research in Geothermal Energy) project, located near Milford, Utah, is a U.S. Department of Energy initiative aimed at advancing EGS technology. In April 2024, eight new stimulation stages (Stages 3R-10) were conducted in well 16A (injection well) subsequent to the first series of stimulation (Stages 1-3) performed in April, 2022. To monitor the induced seismicity, geophones were deployed in wells 58-32, 56-32, and 78B-32, while fiber optic cables were also installed in wells 16B, 78-32, and 78B-32 to collect microseismic data and detect frac hits Preliminary analyses of microseismic catalogs and fiber optic data suggest that the stimulated fractures in Stages 3R–6 closely align with that generated during Stage 3, indicating that the new stimulations were likely reactivating the previously stimulated fracture. To better understand the underlying process, a comprehensive modeling approach that can accurately capture thermal, hydrological, mechanical, and seismic responses is essential. In this work, we propose and utilize a coupled thermo-hydro-mechanical and earthquake (THM+E) simulation workflow to numerically investigate the stimulation activities on well 16A. The specific objective is to confirm whether the new stimulation stages (Stages 3R–6) reactivated fractures previously stimulated during Stage 3. For this purpose, we perform THM+E simulations individually for Stages 3, 3R, 4, and 5, incorporating the discrete fracture networks (DFNs) created by the plane-fitting technique based on the microseismic catalogs. The simulation workflow consists of two separate models: a THM model and an earthquake model, coupled in a one-way manner. Detailed descriptions of the workflow are provided in Section 3. Simulation results are presented in terms of injection pressure, permeability evolution, and predicted seismic catalogs, which are then compared with field data for further analyses. This report is structured as follows. In Section 2, we present detailed analyses of the field data and propose the hypothesis that the new stimulation stages (Stages 3R–6) were probably reactivating the previously stimulated fractures in Stage 3. In Section 3, we introduce the coupled THM+E workflow and the problem setup to validate our hypothesis, followed by the simulation results for each stage in Section 4. Meanwhile, discussions are included to analyze the model predictions and their comparison with field data. Lastly, we conclude the report and outline future plans in Section 5.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

Applying Machine Learning and Bayesian Inference to Identify and Locate Moving Anthropogenic Sources Using Distributed Acoustic Sensing Data

Distributed acoustic sensing (DAS) systems, which use existing telecommunication fibers, offer high‐resolution capabilities ideal for recording anthropogenic sources. However, the complexity of urban environments and the large amount of data recorded by DAS require automated methods to efficiently detect and categorize anthropogenic sources. Here, we evaluate how well three machine learning models (k‐nearest neighbor [k‐NN], convolutional neural networks, and recurrent‐convolutional neural networks) can identify various anthropogenic sources recorded by DAS. Our findings reveal that both k‐NN and neural network methods perform well in high signal‐to‐noise ratio (SNR) settings. However, their accuracy decreases at SNRs <4. We also use Kalman filtering, a form of Bayesian inference, on backprojected locations of these sources to recover locations that generally fall within standard smartphone Global Positioning System errors. By combining machine learning and Kalman filter results, we calculate a multidimensional model of moving anthropogenic sources. These results demonstrate the potential of DAS data in urban seismology for accurately identifying and locating such sources. Depending on the research objectives, these sources can be further studied or filtered out to improve the quality of seismic data for earthquake studies. Such methods provide a valuable tool for urban seismology and seismic hazard analysis.

Luckie, Thomas William [Sandia National Laboratori↗
NASA NTRS2003

Glacier Ice Mass Fluctuations and Fault Instability in Tectonically Active Southern Alaska

Across southern Alaska the northwest directed subduction of the Pacific plate is accompanied by accretion of the Yakutat terrane to continental Alaska. This has led to high tectonic strain rates and dramatic topographic relief of more than 5000 meters within 15 km of the Gulf of Alaska coast. The glaciers of this area are extensive and include large glaciers undergoing wastage (glacier retreat and thinning) and surges. The large glacier ice mass changes perturb the tectonic rate of deformation at a variety of temporal and spatial scales. We estimated surface displacements and stresses associated with ice mass fluctuations and tectonic loading by examining GPS geodetic observations and numerical model predictions. Although the glacial fluctuations perturb the tectonic stress field, especially at shallow depths, the largest contribution to ongoing crustal deformation is horizontal tectonic strain due to plate convergence. Tectonic forces are thus the primary force responsible for major earthquakes. However, for geodetic sites located < 10-20 km from major ice mass fluctuations, the changes of the solid Earth due to ice loading and unloading are an important aspect of interpreting geodetic results. The ice changes associated with Bering Glacier s most recent surge cycle are large enough to cause discernible surface displacements. Additionally, ice mass fluctuations associated with the surge cycle can modify the short-term seismicity rates in a local region. For the thrust faulting environment of the study region a large decrease in ice load may cause an increase in seismic rate in a region close to failure whereas ice loading may inhibit thrust faulting.

SauberRosenberg, Jeanne M.↗
DOE OSTI2024

Heterogeneous mantle effects on the behavior of SmKS waves and outermost core imaging

Seismic traveltime anomalies of waves that traverse the uppermost 100–200 km of the outer core have been interpreted as evidence of reduced seismic velocities (relative to radial reference models) just below the core–mantle boundary (CMB). These studies typically investigate differential traveltimes of SmKS waves, which propagate as P waves through the shallowest outer core and reflect from the underside of the CMB m times. The use of SmKS and S(m-1)KS differential traveltimes for core imaging are often assumed to suppress contributions from earthquake location errors and unknown and unmodelled seismic velocity heterogeneity in the mantle. The goal of this study is to understand the extent to which differential SmKS traveltimes are, in fact, affected by anomalous mantle structure, potentially including both velocity heterogeneity and anisotropy. Velocity variations affect not only a wave's traveltime, but also the path of a wave, which can be observed in deviations of the wave's incoming direction. Since radial velocity variations in the outer core will only minimally affect the wave path, in contrast to other potential effects, measuring the incoming direction of SmKS waves provides an additional diagnostic as to the origin of traveltime anomalies. Here we use arrays of seismometers to measure traveltime and direction anomalies of SmKS waves that sample the uppermost outer core. We form subarrays of EarthScope's regional Transportable Array stations, thus measuring local variations in traveltime and direction. We observe systematic lateral variations in both traveltime and incoming wave direction, which cannot be explained by changes to the radial seismic velocity profile of the outer core. Moreover, we find a correlation between incoming wave direction and traveltime anomaly, suggesting that observed traveltime anomalies may be caused, at least in part, by changes to the wave path and not solely by perturbations in outer core velocity. Modelling of 1-D ray and 3-D wave propagation in global 3-D tomographic models of mantle velocity anomalies match the trend of the observed traveltime anomalies. Overall, we demonstrate that observed SmKS traveltime anomalies may have a significant contribution from 3-D mantle structure, and not solely from outer core structure.

58 GEOSCIENCES↗
DOE OSTI2024

Heterogeneous mantle effects on the behavior of SmKS waves and outermost core imaging

Seismic traveltime anomalies of waves that traverse the uppermost 100–200 km of the outer core have been interpreted as evidence of reduced seismic velocities (relative to radial reference models) just below the core–mantle boundary (CMB). These studies typically investigate differential traveltimes of SmKS waves, which propagate as P waves through the shallowest outer core and reflect from the underside of the CMB m times. The use of SmKS and S(m- 1)KS differential traveltimes for core imaging are often assumed to suppress contributions from earthquake location errors and unknown and unmodelled seismic velocity heterogeneity in the mantle. The goal of this study is to understand the extent to which differential SmKS traveltimes are, in fact, affected by anomalous mantle structure, potentially including both velocity heterogeneity and anisotropy. Velocity variations affect not only a wave’s traveltime, but also the path of a wave, which can be observed in deviations of the wave’s incoming direction. Since radial velocity variations in the outer core will only minimally affect the wave path, in contrast to other potential effects, measuring the incoming direction of SmKS waves provides an additional diagnostic as to the origin of traveltime anomalies. Here we use arrays of seismometers to measure traveltime and direction anomalies of SmKS waves that sample the uppermost outer core. We form subarrays of EarthScope’s regional Transportable Array stations, thus measuring local variations in traveltime and direction. We observe systematic lateral variations in both traveltime and incoming wave direction, which cannot be explained by changes to the radial seismic velocity profile of the outer core. Moreover, we find a correlation between incoming wave direction and traveltime anomaly, suggesting that observed traveltime anomalies may be caused, at least in part, by changes to the wave path and not solely by perturbations in outer core velocity. Modelling of 1-D ray and 3-D wave propagation in global 3-D tomographic models of mantle velocity anomalies match the trend of the observed traveltime anomalies. Overall, we demonstrate that observed SmKS traveltime anomalies may have a significant contribution from 3-D mantle structure, and not solely from outer core structure.

58 GEOSCIENCES↗
DOE OSTI2024

Machine Learning‐Assisted Microearthquake Location Workflow for Monitoring the Newberry Enhanced Geothermal System

Abstract Enhanced geothermal systems (EGS) offer a sustainable energy source but face challenges in accurately locating microearthquakes induced during reservoir stimulation. Locating these microearthquakes provides reliable feedback on the stimulation progress. Current deep learning methods for locating earthquakes require extensive data sets for training, which is problematic as detected microearthquakes are often limited. To address the scarcity of training data, we propose a practical workflow using probabilistic multilayer perceptron (PMLP) which predicts microearthquake locations from cross‐correlation time lags in waveforms. Utilizing a 3D velocity model of Newberry site derived from ambient noise interferometry, we generate numerous synthetic microearthquakes and 3D acoustic waveforms for PMLP training. Accurate synthetic tests prompt us to apply the trained network to the 2012 and 2014 stimulation field waveforms. To enhance the accuracy of source localization, we carefully handpick the P‐arrival times. Predictions on the 2012 stimulation data set show major microseismic activity at depths of 0.5–1.2 km, correlating with a known casing leakage scenario. In the 2014 data set, the majority of predictions concentrate at 2.0–2.9 km depths, consistent with results obtained from conventional physics‐based inversion, and align with the presence of natural fractures from 2.0 to 2.7 km. We validate our findings by comparing the synthetic and field picks, demonstrating a satisfactory match for the first arrivals. By combining the benefits of quick inference speeds and accurate location predictions, we demonstrate the feasibility of using realistic synthetic data set to locate microseismicity for EGS monitoring.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Bayesian Optimal Sensor Augmentation Via Estimated Mutual Information

We consider the problem of designing a sensor network to most efficiently locate the source of a seismic event. The left panel of Figure 1 shows 982 possible sensor locations in a region of interest. The objective is to infer the unknown location of an earthquake, indicated by the gold dot for illustration purposes. Assume that three sensors have already been placed for data collection. Travel times of the P-wave from the source to the three sensors are observed and used to infer the unknown source location via a Bayesian statistical model described subsequently. The right panel of Figure 1 shows samples from the resulting posterior distribution for source location along with a posterior estimate of this location.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Sequences of seismic and aseismic slip on bimaterial faults show dominant rupture asymmetry and potential for elevated seismic hazard

Here we perform numerical simulations of sequences of earthquake and aseismic slip on planar rate and state faults separating dissimilar material within the 2-D plane strain approximation. We resolve all stages of the earthquake cycle from aseismic slip to fast ruptures while incorporating full inertia effects during seismic event propagation. We show that bimaterial coupling results in favorable nucleation site and subsequent asymmetric rupture propagation. We demonstrate that increasing the material contrast enhances this asymmetry leading to higher slip rates and normal stress drops in the preferred rupture propagation direction. The normal stress drop, induced by the bimaterial effect, leads to strong dynamic weakening of the fault and may destabilize the creeping region on a heterogeneous rate and state fault, resulting in extended rupture propagation. Such rupture penetration into creeping patches may lead to more frequent opening of earthquake gates, causing increased seismic hazard. Furthermore, bimaterial coupling may lead to irregular seismicity pattern in terms of event length, peak slip rates, and hypocenter location, depending on the properties of the creeping patches bordering the seismogenically active part of the fault. Our results highlight robust characteristics of bimaterial interfaces that persist over long sequence of events and suggest the need for further exploration of the role of material contrast in earthquake physics and models of seismic hazard.

58 GEOSCIENCES↗
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↗
NASA NTRS1998

GPS Monitoring of Surface Change During and Following the Fortuitous Occurrence of the M(sub w) = 7.3 Landers Earthquake in our Network

Accomplishments: (1) Continues GPS monitoring of surface change during and following the fortuitous occurrence of the M(sub w) = 7.3 Landers earthquake in our network, in order to characterize earthquake dynamics and accelerated activity of related faults as far as 100's of kilometers along strike. (2) Integrates the geodetic constraints into consistent kinematic descriptions of the deformation field that can in turn be used to characterize the processes that drive geodynamics, including seismic cycle dynamics. In 1991, we installed and occupied a high precision GPS geodetic network to measure transform-related deformation that is partitioned from the Pacific - North America plate boundary northeastward through the Mojave Desert, via the Eastern California shear zone to the Walker Lane. The onset of the M(sub w) = 7.3 June 28, 1992, Landers, California, earthquake sequence within this network poses unique opportunities for continued monitoring of regional surface deformation related to the culmination of a major seismic cycle, characterization of the dynamic behavior of continental lithosphere during the seismic sequence, and post-seismic transient deformation. During the last year, we have reprocessed all three previous epochs for which JPL fiducial free point positioning products available and are queued for the remaining needed products, completed two field campaigns monitoring approx. 20 sites (October 1995 and September 1996), begun modeling by development of a finite element mesh based on network station locations, and developed manuscripts dealing with both the Landers-related transient deformation at the latitude of Lone Pine and the velocity field of the whole experiment. We are currently deploying a 1997 observation campaign (June 1997). We use GPS geodetic studies to characterize deformation in the Mojave Desert region and related structural domains to the north, and geophysical modeling of lithospheric behavior. The modeling is constrained by our existing and continued GPS measurements, which will provide much needed data on far-field strain accumulation across the region and on the deformational response of continental lithosphere during and following a large earthquake, forming the basis for kinematic and dynamic modeling of secular and seismic-cycle deformation. GPS geodesy affords both regional coverage and high precision that uniquely bear on these problems.

Miller, M. Meghan↗
NASA NTRS2012

Simplex GPS and InSAR Inversion Software

Changes in the shape of the Earth's surface can be routinely measured with precisions better than centimeters. Processes below the surface often drive these changes and as a result, investigators require models with inversion methods to characterize the sources. Simplex inverts any combination of GPS (global positioning system), UAVSAR (uninhabited aerial vehicle synthetic aperture radar), and InSAR (interferometric synthetic aperture radar) data simultaneously for elastic response from fault and fluid motions. It can be used to solve for multiple faults and parameters, all of which can be specified or allowed to vary. The software can be used to study long-term tectonic motions and the faults responsible for those motions, or can be used to invert for co-seismic slip from earthquakes. Solutions involving estimation of fault motion and changes in fluid reservoirs such as magma or water are possible. Any arbitrary number of faults or parameters can be considered. Simplex specifically solves for any of location, geometry, fault slip, and expansion/contraction of a single or multiple faults. It inverts GPS and InSAR data for elastic dislocations in a half-space. Slip parameters include strike slip, dip slip, and tensile dislocations. It includes a map interface for both setting up the models and viewing the results. Results, including faults, and observed, computed, and residual displacements, are output in text format, a map interface, and can be exported to KML. The software interfaces with the QuakeTables database allowing a user to select existing fault parameters or data. Simplex can be accessed through the QuakeSim portal graphical user interface or run from a UNIX command line.

Donnellan, Andrea↗
DOE OSTI2023

Pore Pressure Diffusion and Onset of Induced Seismicity

The Raton Basin of Southern Colorado and Northern New Mexico experienced a significant increase in seismicity over the last several decades, including two M5+ earthquakes. The increase in seismicity started soon after an increase in wastewater disposal associated with coal-bed methane production. Twenty-nine wastewater injection wells have operated in the Basin, with several high-rate injectors located within a kilometer of the most seismically active regions. We developed a pore pressure model to investigate the relationship between pore pressure diffusion and onset of induced seismicity, with an emphasis on the Trinidad fault zones. Our model indicates that sufficient pore pressure increase had accumulated on the northern extent of the Trinidad fault zone to induce seismicity at the time of early-stage seismic activity. The hydrogeologic architecture of the Basin played a key role in controlling pore pressure perturbation. The basal sedimentary interval with moderately low diffusivity restricted vertical pore pressure diffusion, creating a substantial lag time between wells becoming active and the expression of perturbation in the fault zones. The spatiotemporal relationship between early-stage seismicity and wastewater injection in conjunction with modeled pore pressures suggests that the early-stage seismicity in the Raton Basin was likely triggered by pore pressure diffusion.

58 GEOSCIENCES↗
DOE OSTI2021

Widespread deep seismicity in the Delaware Basin, Texas, is mainly driven by shallow wastewater injection

Significance The Delaware Basin, western Texas, has experienced a surge in the number of earthquakes. We investigate basin-wide earthquake, hydrogeologic, industrial, and geodetic data spanning 1993–2020 to identify the reasons for the increase. We use these data and models to quantify the stress changes from wastewater injection and hydrocarbon extraction. We demonstrate that the basin-wide seismicity is dominated by the stresses caused by shallow injection that are transmitted vertically by deforming the underlying crust. We analyze both injection and production data to understand how subsurface properties control earthquake locations and timings. We suggest that injecting wastewater into porous formations will minimize seismic hazard. We establish a framework explaining why earthquakes sometimes occur and sometimes do not occur within human activity-impacted areas.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Potential Link Between 2020 Mentone, West Texas M5 Earthquake and Nearby Wastewater Injection: Implications for Aquifer Mechanical Properties

Abstract The M5 Mentone earthquake that occurred on March 26, 2020, was the largest event recorded over the last 2 decades in West Texas within the Delaware Basin, a U.S. major petroleum‐producing area. Also, numerous hydrofracturing and wastewater disposal wells are spread across this region. Within a 30 km distance to mainshock, eight class‐II injection wells for industrial wastewater disposal target the deep porous Ellenburger aquifer at an average rate of 1.36 × 10 6 barrel (BBL) per month during 2012–2020. Poroelastic models of fluid diffusion show these nearby injectors collectively imparted up to 80.5 kPa of Coulomb stress at the mainshock location, capable of triggering this M5 event. Assuming the Mentone event occurs when pore‐pressure increase is maximum, the time delay between peak injection and the M5 occurrence corresponds with an optimal permeability of 6.76 × 10 ‐14 m 2 for the Ellenburger aquifer layer, in agreement with independent estimates.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Crust and Upper Mantle Structure Beneath the Eastern United States

Abstract The Eastern United States (EUS) has a complex geological history and hosts several seismic active regions. We investigate the subsurface structure beneath the broader EUS. To produce reliable images of the subsurface, we simultaneously invert smoothed P‐wave receiver functions, Rayleigh‐wave phase and group velocity measurements, and Bouguer gravity observations for the 3D shear‐wave speed. Using surface‐wave observations (3–250 s) and spatially smoothed receiver functions, our velocity models are robust, reliable, and rich in detail. The shear‐wave velocity models fit all three types of observations well. The resulting velocity model for the eastern U.S. shows thinner crust beneath New England, the east coast, and the Mississippi Embayment (ME). A relatively thicker crust was found beneath the stable North America craton. A relatively slower upper mantle was imaged beneath New England, the east coast, and western ME. A comparison of crust thickness derived from our model against four recent published models shows first‐order consistency. A relatively small upper mantle low‐speed region correlates with a published P‐wave analysis that has associated the anomaly with a 75 Ma kimberlite volcanic site in Kentucky. We also explored the relationship between the subsurface structure and seismicity in the eastern U.S. We found that earthquakes often locate near regions with seismic velocity variations, but not universally. Not all regions of significant subsurface wave speed changes are loci of seismicity. A weak correlation between upper mantle shear velocity and earthquake focal mechanism has been observed.

58 GEOSCIENCES↗