Search NASASearch

SEARCH · Search NASA

Results for “microearthquakes”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

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

Deep learning forecasts the spatiotemporal evolution of fluid-induced microearthquakes

Microearthquakes generated by subsurface fluid injection record the evolving stress state and permeability of reservoirs. Forecasting their spatiotemporal evolution is therefore critical for applications such as enhanced geothermal systems, carbon dioxide sequestration and other geoengineering applications. Here we propose a transformer neural network model that ingests hydraulic stimulation history and prior microearthquake observations to forecast four key quantities: cumulative microearthquake count, cumulative logarithmic seismic moment, and the 50th- and 95th-percentile extents of the microearthquake cloud. Applied to the EGS Collab Experiment 1 dataset, the model achieves R2 > 0.98 for the 1-s forecast horizon and R2 > 0.88 for the 15-s forecast horizon across all targets, and supplies uncertainty estimates through a learned standard deviation term. These accurate, uncertainty-quantified forecasts enable real-time inference of fracture propagation and permeability evolution, demonstrating the strong potential of deep-learning approaches to improve seismic-risk assessment and guide mitigation strategies in future fluid-injection operations.

Chung, Jaehong

Tectonic implications of the microearthquake seismicity and fault plane solutions in southern Peru

Because the contortion in the seismic zone in southern Peru is aligned approximately parallel to the direction of relative plate motion, rather than perpendicular to the coast of Peru, the position of the contortion need not migrate with respect to the overriding South American plate as the Nazca plate subducts beneath it, and the flow in the surrounding asthenosphere could be in a steady state. In addition, the position of the contortion defines the northern boundary of the volcanic arc in southern Peru. The inference that a wedge of asthenospheric material must overlie the downgoing slab for subduction-related volcanism to occur is thereby strengthened.

Grange, F.

Tectonic deformation of the Andes and the configuration of the subducted slab in central Peru: Results from a micro-seismic experiment

The vast majority of the microearthquakes recorded occurred to the east: on the Huaytapallana fault in the Eastern Cordillera or in the western margin of the sub-Andes. The sub-Andes appear to be the physiographic province subjected to the most intense seismic deformation. Focal depths for the crustal events here are as deep as 50 km, and the fault plane solutions, show thrust faulting on steep planes oriented roughly north-south. The Huaytapallana fault in the Cordillera Oriental also shows relatively high seismicity along a northeast-southwest trend that agrees with the fault scarp and the east dipping nodal plane of two large earthquakes that occurred on this fault in 1969. The recorded microearthquakes of intermediate depth show a flat seismic zone about 25 km thick at a depth of about 100 km. This agrees with the suggestion that beneath Peru the slab first dips at an angle of 30 deg to a depth of 100 km and then flattens following a quasi-horizontal trajectory. Fault plane solutions of intermediate depth microearthquakes have horizontal T axes oriented east-west.

Suarez, G.

Machine Learning Approaches to Predicting Induced Seismicity and Imaging Geothermal Reservoir Properties

This project developed machine learning (ML) methods, lab data sets, and field data to advance geothermal exploration and geothermal energy production. The work had three focus areas. One involved the development of ML methods to use microearthquakes (MEQs) for imaging geothermal reservoir properties and improving subsurface characterization – most importantly the evolution of permeability within the evolving reservoir. This part of the work included development of ML approaches for automated MEQ location, focal mechanism determination and identification of earthquake precursors. The second area focused on using MEQ signals generated by geothermal exploration and production to predict the relationship between fluid injection and seismicity. Here, we extended to reservoir scale our success in using ML to predict laboratory earthquakes and fault zone stress state. The third focus area was on lab experiments. Here, we developed new ML models for lab earthquake prediction and identification of precursors to failure to improve earthquake forecasting and early warning in geothermal settings. Major outcomes of our work include ML models that learn from MEQ signals during geothermal exploration and production to predict induced seismicity. MEQs occur naturally in connection with drilling and energy production. We developed ML methods to use the seismic waves from these events to characterize the elastic, hydraulic and poromechanical properties of reservoirs. Our work illuminated fracture geometry and the evolution of fracture permeability by incorporating seismic coda wave analysis and ML methods to relate fluid injection and seismicity. We significantly expanded laboratory earthquake prediction to include methods that use both passive measurements of microearthquakes within the lab fault zones and also active source acoustic measurements of fault zone elastic properties. These methods can now predict fault zone stress state, time to failure and the magnitude of lab earthquakes. Our work showed that repetitive stick- slip failure events during frictional sliding (the lab equivalent of earthquakes) are preceded by a cascade of micro-failure events that radiate energy in a manner that foretells unstable failure – manifest as laboratory MEQs. We documented a mapping between fracture properties and statistical attributes of elastic radiation. We extended existing works to geothermal reservoir scale and developed ML methods to determine reservoir permeability, fracture properties, and their evolution during geothermal energy production. An attractive feature of ML algorithms is their ability to handle big datasets and reveal patterns and correlations that may remain invisible to conventional analyses. Our work connected data from field, laboratory and intermediate scales to study permeability, stress, strength, fracture stiffness and geometry. At the field scale we used data from the Newberry Volcano field site, UtahFORGE, EGS Collab, and also the Bedretto underground research lab in Switzerland. These data sets are bridging the gap between the lab scale, theory, and reservoir scale. Our work produced plain language summaries to improve public understanding of DOE research. We also developed openly distributed ML and seismicity datasets for use by all researchers and we published connections between induced seismicity in geothermal areas and reservoir properties including permeability, fracture properties, and stress state. Our models are designed for the large data sets of induced seismicity typically associated with geothermal sites. We produced labeled event catalogs and used them on geothermal data to assess how ML can facilitate geothermal production and exploration. All datasets are available on the GDR Productivity: The project produced 32 publications in peer reviewed journals (two are in review). It supported the work of 6 PhD students, 40 conference presentations, 6 keynote talks at national meetings, and mentoring and professional development for 4 postdoctoral fellows.

15 GEOTHERMAL ENERGY

Strain buildup and release, earthquake prediction and selection of VLB sites for the margins of the North Pacific

Microearthquake studies of the southern section of the San Andreas fault system were conducted in order to elucidate the role of known but little studied complications in the fault system which could affect the SAFE measurement. Data reduction from microearthquake studies in Baja California, Mexico, was completed. In addition, styles of deformation on the Alpine fault in New Zealand and on the San Andreas fault were compared. Sections of these faults with comparable physical characteristics seem to deform in a similar manner, indicating that deformation style is controlled by certain fundamental relations and will continue in the future. The geology of eastern New Guinea is also discussed. Deformation in the area is proposed to be related to collision of a series of island arcs with central New Guinea, which is part of the Australian continent. Among other unusual properties, eastern New Guinea contains the highest anomaly of the gravimetrically determined geoid based on GEM 6 and surface gravity data.

Scholz, C. H.

Complementary Subsurface Characterization Methods to Develop a Geologic Model for the EGS Collab Experiment, Sanford Underground Research Facility

The EGS (Enhanced Geothermal Systems) Collab project was performed within the Sanford Underground Research Facility (SURF) with a goal of understanding processes and evaluation of models related to hydraulic stimulation of rock at depth. The present work deals with the development of Testbed 2 where experiments were conducted at a depth of 1.25 km and were located within a well-characterized testbed in a metamorphic, amphibolite host rock. A total of eleven boreholes varying in length between 10.6 m and 81.2 m were continuously cored to develop the testbed. In addition to the continuous coring of the amphibolite host rock, geophysical instrumentation supporting electrical resistivity tomography (ERT), microearthquake (MEQ) detection, and optical fiber providing distributed temperature sensing (DTS) and distributed acoustic sensing (DAS) were installed in the monitoring boreholes, all of which produced a comprehensive complementary suite of characterization and monitoring technologies. Pre-stimulation characterization of the groundwater conditions identified only two hydraulically significant fractures, neither of which transects the central portion of the testbed. Flow and pressure monitoring indicated low preexisting pore pressure conditions likely affected by the mine openings.

15 GEOTHERMAL ENERGY

Complementary Subsurface Characterization Methods to Develop a Geologic Model for the EGS Collab Experiment, Sanford Underground Research Facility

The EGS (Enhanced Geothermal Systems) Collab project was performed within the Sanford Underground Research Facility (SURF) with a goal of understanding processes and evaluation of models related to hydraulic simulation of rock at depth. The present work deals with the development of Testbed 2 where experiments were conducted at a depth of 1.25 km and were located within a well-characterized testbed in a metamorphic, amphibolite host rock. A total of eleven boreholes varying in length between 10.6 m and 81.2 m were continuously cored to develop the testbed. In addition to the continuous coring of the amphibolite host rock, geophysical instrumentation supporing electrical resisivity tomography (ERT), microearthquake (MEQ) detection, and optical fiber providing distributed temperature sensing (DTS) and distributed acoustic sensing (DAS) were installed in the monitoring boreholes, all of which produced a comprehensive complementary suite of characterization and monitoring technologies.

15 GEOTHERMAL ENERGY

DASEventNet: AI‐Based Microseismic Detection on Distributed Acoustic Sensing Data From the Utah FORGE Well 16A (78)‐32 Hydraulic Stimulation

Abstract Distributed acoustic sensing (DAS) has emerged as a promising seismic technology for monitoring microearthquakes (MEQs) with high spatial resolution. Efficient algorithms are needed for processing large DAS data volumes. This study introduces a deep learning (DL) model based on a Residual Convolutional Neural Network (ResNet) for detecting MEQs using DAS data, named as DASEventNet. The test data were collected from the Utah FORGE 16A (78)‐32 hydraulic stimulation experiments conducted in April 2022. The DASEventNet model achieves a remarkable accuracy of 100% when discriminating MEQs from noise in the raw test set of 260 examples. Surprisingly, the model identified weak MEQ signatures that have been manually categorized as noise. The decision‐making process with the model is decoded by the classic activation map, which illuminates learning features of the DASEventNet model. These features provide clear illustrations of weak MEQs and varied noise types. Finally, we apply the trained model to the entire period (∼7 days) of continuous DAS recordings and find that it discovers >5,700 new MEQs, previously unregistered in the public Silixa DAS catalog. The DASEventNet model significantly outperforms the traditional seismic method Short‐Term Average/Long‐Term Average (STA/LTA), which detected only 1,307 MEQs. The DASEventNet detection threshold is M w −1.80 compared to the minimum magnitude of M w −1.14 detected by STA/LTA. The spatiotemporal distribution of the newly identified MEQs defines an extensive stimulation zone and more accurately characterizes fracture geometry. Our results highlight the potential of DL for long‐term, real‐time microseismic monitoring that can improve enhanced geothermal systems and other activities that include subsurface hydraulic fracturing.

15 GEOTHERMAL ENERGY

Co‐Evolution of Specific Stiffness and Permeability of Rock Fractures Offset in Shear

Abstract Fractures and faults represent planes of weakness and compliance in rock masses that serve as focal points for both microearthquakes and fluid transport, with seismicity and permeability evolution closely linked. Contact stiffness is highly stress‐sensitive and directly influences permeability. We explore the co‐evolution of specific stiffness and permeability of rough fractures under normal stress and shear offset using numerical simulations. Individual rough fractures are represented by variable amplitude (Root mean square) and wavelength ( λ ) using a granular mechanics model. Contacting rough surfaces are mated, offset in shear, and then compacted in displacement mode. The compacting fractures generate stress‐dependent changes in contact porosity, which govern both permeability and stiffness evolution. We establish a universal dimensionless relationship linking specific stiffness and permeability that inherently incorporates the effects of surface roughness, shear offset, and microcracking. The observed cracking effect—where local stress redistribution and pressure‐driven microcrack propagation dynamically alter the aperture field—introduces a nonlinear permeability response at high stress. Increased roughness amplitude and larger shear offsets reduce stiffness while dampening permeability sensitivity to stress, demonstrating a strong interplay between surface texture and hydro‐mechanical behavior. While the model captures this behavior effectively, deviations emerge at very low porosities due to extreme aperture sensitivity in this limit.

He, Xinxin [Energy and Mineral Engineering Pennsyl

Fault Slip and Fluid Flow: Seismic Source Analysis to Assess Role of Multiple Slip Patches in Fault Permeability

The relationship between fault reactivation, microearthquakes (MEQs), and permeability evolution during fluid injection plays a critical role in energy harvesting and waste disposal. Recent studies have demonstrated the possibility of predicting fault permeability using cumulative seismic moments of MEQs quantitatively. To understand the underlying physical processes, we conduct fault reactivation experiments using Utah FORGE granitoid and analyze acoustic emission (AE) signals generated during stepwise increases in fluid injection pressure. Frequency analysis of thousands of calibrated AE signals reveals that fault reactivation produces multiple AE source patches with millimeter-scale radii—smaller than the sample fault radius. The cumulative area of the reactivated patches covers the fault multiple times over (∼10x–50x area) for each pressure step. These findings provide mechanistic insight that measured permeability enhancement is not driven by a single large slip event, but by the sequential and interacting activation of multiple slip patches that create a continuous flow pathway.

Nurshal, M. E. M. [Pennsylvania State University,

Importance of Considering Near-Surface Attenuation in Earthquake Source Parameter Estimation: Insights from Kappa at a Dense Array in Oklahoma

ABSTRACT Separating earthquake source spectra from propagation effects is challenging. The propagation effect contains a site-dependent term related to the high attenuation of shallow sediments. Neglecting the site-dependent attenuation can cause large biases and scattering in the corner-frequency (fc) estimates, resulting in significant stress-drop deviations. In this study, we investigate shallow attenuation at the LArge-n Seismic Survey in Oklahoma (LASSO) and site-related biases and scattering in source parameter measurements due to simplified attenuation models. We measure the high-frequency spectral decay parameter kappa on the vertical acceleration spectra of regional earthquakes (125 km away). The site-dependent kappa (κ0,acc) suggests that attenuation increases rapidly at shallow depth and is highly site-dependent. 10%–75% of the attenuation is site-dependent for S waves and even larger for P waves. The quality factor for S waves (QS) ranges from 10 to 100 in the upper 400 m. QP for P waves is mostly below 10 within the same depth. The Quaternary sediments tend to be more attenuating (QS<30), but the Permian rocks also can have high attenuation. We demonstrate that using a non-site-dependent attenuation model in single-spectra fitting leads to large scattering in fc estimates among stations with apparent good fits. The apparent fc can significantly deviate when the range of site-dependent kappa is large or with a higher assumed source spectral fall-off rate n. The biases in apparent fc depend on site condition and distance; however, the correlation between fc and these factors might not be obvious, depending on model assumptions. An apparent increase of stress drop with magnitude in a previous study for local microearthquakes (1.3

Chang, Hilary

Satellite relay telemetry in the surveillance of active volcanoes and major fault zones

A review was made of efforts to develop a dense telemetered microearthquake network to study earthquake mechanics along the San Andreas fault and the strain mechanics of the Kilauea Volcano. The principle elements and objectives of the ERTS-A proposal are outlined. Some of the aspects of the earthquake network and the results obtained from it as well as some promising experiments in computerized record processing are discussed.

Eaton, J. P.

Heat flow and geothermal potential of the East Mesa KGRA, Imperial Valley, California

The East Mesa KGRA (Known Geothermal Resource Area) is located in the southeast part of the Imperial Valley, California, and is roughly 150 kilometers square in areal extent. A new heat flow technique which utilizes temperature gradient measurements across best clays is presented and shown to be as accurate as conventional methods for the present study area. Utilizing the best clay gradient technique, over 70 heat flow determinations have been completed within and around the East Mesa KGRA. Background heat flow values range from 1.4 to 2.4 hfu (1 hfu = .000001 cal. per square centimeter-second) and are typical of those throughout the Basin and Range province. Heat flow values for the northwest lobe of the KGRA (Mesa anomaly) are as high as 7.9 hfu, with the highest values located near gravity and seismic noise maxima and electrical resistivity minima. An excellent correlation exists between heat flow contours and faults defined by remote sensing and microearthquake monitoring.

Swanberg, C. A.

A study of microseismicity in northern Baja California, Mexico

Five microearthquake instruments were operated for 2 months in 1974 in a small mobile array deployed at various sites near the Agua Blanca and San Miguel faults. An 80-km-long section of the San Miguel fault zone is presently active seismically, producing the vast majority of recorded earthquakes. Very low activity was recorded on the Agua Blanca fault. Events were also located near normal faults forming the eastern edge of the Sierra Juarez suggesting that these faults are active. Hypocenters on the San Miguel fault range in depth from 0 to 20 km although two-thirds are in the upper 10 km. A composite focal mechanism showing a mixture of right-lateral and dip slip, east side up, is similar to a solution obtained for the 1956 San Miguel earthquake which proved consistent with observed surface deformation.

Johnson, T. L.

Three-dimensional P and S wave velocity structures of southern Peru and their tectonic implications

Arrival times of compressional and shear (S) waves from microearthquakes recorded in 1981 by an 18-station regional array are used to study the three-dimensional velocity structure of the crust and upper mantle of the central Andes. The data suggest a crustal thickness of about 40 km beneath the coast, increasing to about 70 km beneath the Cordillera Occidental. The inverse correlation between the dip in the Moho and the dip of the slab may indicate a broad-scale causal relation between the two. S wave velocities in the mantle between 70 and 130 km depth above the 30-degree dipping slab are low, possibly indicating the presence of a partially melted asthenosphere that may be responsible for the magmatic activity recorded in southern Peru.

Cunningham, Paul S.

Terrestrial geophysics in the SeaRISE project

Some areas of research in the SeaRISE project are briefly discussed. They are as follows: (1) Radar Sounding serves multiple purposes. The most general and obvious is mapping ice thickness and the surface and bedrock topography of the ice sheet. (2) The purpose of Seismic Shooting, in addition to water depth measurements on floating ice, is to provide information about the internal physical characteristics of the ice sheet, the rock beneath it, and the interface between the two. (3) Passive Seismic monitoring of microearthquakes can be used to study brittle fracture within the ice or the rock beneath it. Common parameters available from these studies are fault location, orientation, and displacement, as well as the size of the rupture area, stress drop, and energy released. (4) There is a large contrast in Electrical Resistivity between ice or permafrost on the one hand and liquid water or wet rock on the other hand. Thus, electrical resistivity profiles have the ability of revealing the depth to the melting point, whether it is found at the base of the ice or in the subglacial rock. (5) Gravity anomalies, especially combined with seismic measurements, are an effective tool for determining deeper crustal structure. Anomalies averaged over extensive areas are useful also for their potential to reveal isostatic imbalance, which is a measure of average glacial change over the last several hundred years.

Bentley, C. R.