Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “earthquake forecasting”

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

Foreshock properties illuminate nucleation processes of slow and fast laboratory earthquakes

Understanding the connection between seismic activity and the earthquake nucleation process is a fundamental goal in earthquake seismology with important implications for earthquake early warning systems and forecasting. We use high-resolution acoustic emission (AE) waveform measurements from laboratory stick-slip experiments that span a spectrum of slow to fast slip rates to probe spatiotemporal properties of laboratory foreshocks and nucleation processes. We measure waveform similarity and pairwise differential travel-times (DTT) between AEs throughout the seismic cycle. AEs broadcasted prior to slow labquakes have small DTT and high waveform similarity relative to fast labquakes. We show that during slow stick-slip, the fault never fully locks, and waveform similarity and pairwise differential travel times do not evolve throughout the seismic cycle. In contrast, fast laboratory earthquakes are preceded by a rapid increase in waveform similarity late in the seismic cycle and a reduction in differential travel times, indicating that AEs begin to coalesce as the fault slip velocity increases leading up to failure. These observations point to key differences in the nucleation process of slow and fast labquakes and suggest that the spatiotemporal evolution of laboratory foreshocks is linked to fault slip velocity.

58 GEOSCIENCES↗

A proposal to conduct a Caribbean plate project involving the application of space technology to the study of Caribbean geology

The Caribbean plate project is designed to improve current understanding of geological resources and geological hazards within the Caribbean region. Models of mineral occurrence and genesis (including energy resources) on a regional scale, which contribute to nonrenewable resource investigations. Models of lithospheric stress and strain on a regional scale, which contribute to forecasting geological hazards such as earthquakes and major volcanic eruptions are developed. Geological information is synthesize, and research tools provided by space technology the study of the Earth's crust are used. The project was organized in a thematic fashion, to focus on specific geological aspects of the Caribbean plate which are considered to be key factors in developing the types of models described. The project adopts a synoptic perspective in seeking to characterize the three dimensional structure, composition, state of stress, and evolution of the entire Caribbean plate. Geological information derived from analysis of space acquired data is combined with information provided by conventional methods to obtain insight into the structure, composition, and evolution of the Earth's crust. In addition, very long baseline interferometry and laser ranging techniques, which are also based upon the use of space technology, obtain information concerning crustal motion that, in turn, provides insight into the distribution and localization of crustal stress.

Wadge, G.↗

Cyclic Injection Leads to Larger and More Frequent Induced Earthquakes under Volume-Controlled Conditions

As carbon storage technologies advance globally, methods to understand and mitigate induced earthquakes become increasingly important. Although the physical processes that relate increased subsurface pore pressure changes to induced earthquakes have long been known, reliable methods to forecast and control induced seismic sequences remain elusive. Suggested reservoir engineering scenarios for mitigating induced earthquakes typically involve modulation of the injection rate. Some operators have implemented periodic shutdowns (i.e., effective cycling of injection rates) to allow reservoir pressures to equilibrate (e.g., Paradox Valley) or shut-in wells after the occurrence of an event of concern (e.g., Basel, Switzerland). Other proposed scenarios include altering injection rates, actively managing pressures through coproduction of fluids, and preinjection brine extraction. Here, in this work, we use 3D physics-based earthquake simulations to understand the effects of different injection scenarios on induced earthquake rates, maximum event magnitudes, and postinjection seismicity. For comparability, the modeled injection considers the same cumulative volume over the project’s operational life but varies the schedule and rates of fluid injected. Simulation results show that cyclic injection leads to more frequent and larger events than constant injection. Furthermore, with intermittent injection scenario, a significant number of events are shown to occur during pauses in injection, and the seismicity rate remains elevated for longer into the postinjection phase compared to the constant injection scenario.

58 GEOSCIENCES↗

Earthquakes and sea level - Space and terrestrial metrology on a changing planet

A review is presented of the stability and scale of crustal deformation metrology which has particular relevance to monitoring deformation associated with sea level and earthquakes. Developments in space geodesy and crustal deformation metrology in the last two decades have the potential to acquire a homogeneous global data set for monitoring relative horizontal and vertical motions of the earth's surface to within several millimeters. New tools discussed for forecasting sea level rise and damaging earthquakes include: very long baseline interferometry, satellite laser ranging, the principles of GPS geodesy, and new sea level sensors. Space geodesy permits a unified global basis for future metrology of the earth, and the continued availability of the GPS is currently fundamental to this unification.

Bilham, Roger↗

Forecasting induced seismicity in Oklahoma using machine learning methods

Oklahoma earthquakes in the past decade have been mostly associated with wastewater injection. Here we use a machine learning technique—the Random Forest to forecast induced seismicity rate in Oklahoma based on injection-related parameters. We split the data into training (2011.01–2015.05) and test (2015.06–2020.12) periods. The model forecasts seismicity rate during the test period based on input features, including operational parameters (injection rate and pressure), geological information (depth to basement), and modeled pore pressure and poroelastic stress. The results show overall good match with observed seismicity rate (adjusted R 2 of 0.75). The model shows that pore pressure rate and poroelastic stressing rates are the two most important features in forecasting. The absolute values of pore pressure and poroelastic stress, and the injection rate itself, are less important than the stressing rates. These findings further emphasize that temporal changes of stressing rates would lead to significant changes in seismicity rates.

58 GEOSCIENCES↗

Understanding Earthquake Fault Systems Using QuakeSim Analysis and Data Assimilation Tools

We are using the QuakeSim environment to model interacting fault systems. One goal of QuakeSim is to prepare for the large volumes of data that spaceborne missions such as DESDynI will produce. QuakeSim has the ability to ingest distributed heterogenous data in the form of InSAR, GPS, seismicity, and fault data into various earthquake modeling applications, automating the analysis when possible. Virtual California simulates interacting faults in California. We can compare output from long time history Virtual California runs with the current state of strain and the strain history in California. In addition to spaceborne data we will begin assimilating data from UAVSAR airborne flights over the San Francisco Bay Area, the Transverse Ranges, and the Salton Trough. Results of the models are important for understanding future earthquake risk and for providing decision support following earthquakes. Improved models require this sensor web of different data sources, and a modeling environment for understanding the combined data.

earthquake simulations↗

Tsunami Early Warning From Global Navigation Satellite System Data Using Convolutional Neural Networks

Abstract We investigate the potential of using Global Navigation Satellite System (GNSS) observations to directly forecast full tsunami waveforms in real time. We train convolutional neural networks to use less than 9 min of GNSS data to forecast the full tsunami waveforms over 6 hr at select locations, and obtain accurate forecasts on a test data set. Our training and test data consists of synthetic earthquakes and associated GNSS data generated for the Cascadia Subduction Zone using the MudPy software, and corresponding tsunami waveforms in Puget Sound computed using GeoClaw. We use the same suite of synthetic earthquakes and waveforms as in earlier work where tsunami waveforms were used for forecasting, and provide a comparison. We also explore varying the number of GNSS stations, their locations, and their observation durations.

Rim, Donsub↗

Characterizing HVDC Transmission Flexibility under Extreme Operating Conditions

System operators rely on system flexibility, traditionally mainly from generation, to handle unexpected reliability and resilience events, ranging from excessive resource forecast errors to extreme events like heatwaves, earthquakes, and cyberattacks. Flexible transmission, such as controllable high voltage direct current (HVDC) transmission systems present an opportunity to increase overall system flexibility to accommodate operational challenges. This paper provides a methodology to study contributions to system flexibility by controllable, power electronics based transmission. The Western Electricity Coordinating Council (WECC) system is used as an example to study contributions from existing and future HVDC lines. Under extreme system conditions, it is identified that HVDC transmission flexibility can contribute with 24.8 to 28% of avoided unserved energy, and in some areas, the benefits amount to 50 to 70%.

HVDC transmission, PCM, Balancing authorities↗

Neural network wavelet technology: A frontier of automation

Neural networks are an outgrowth of interdisciplinary studies concerning the brain. These studies are guiding the field of Artificial Intelligence towards the, so-called, 6th Generation Computer. Enormous amounts of resources have been poured into R/D. Wavelet Transforms (WT) have replaced Fourier Transforms (FT) in Wideband Transient (WT) cases since the discovery of WT in 1985. The list of successful applications includes the following: earthquake prediction; radar identification; speech recognition; stock market forecasting; FBI finger print image compression; and telecommunication ISDN-data compression.

Szu, Harold↗

ORION: Operational FoRecastIng Of INduced Seismicity

The Operational Forecasting of Induced Seismicity toolkit “ORION” (ORION) is an open-source, observation-based ensemble forecasting toolkit which is geared towards helping operators understand the seismic hazard (i.e., probabilistic assessment of the magnitude and frequency of induced seismic events) at a site. ORION analyzes how the seismic hazard evolves during injection and suggests possible mitigation strategies to employ if an earthquake that exceeds certain threshold is observed. Through its ensemble modeling approach, ORION leverages the benefits of statistical-, physics-, and machine learning-based forecasting methodologies, while reducing the impact of each model’s respective limitations. The ORION toolkit consists of an easy-to-use GUI interface that affords a user as much or as little interaction as desired. Advanced capabilities allow the user to upload local, high-precision earthquake catalogs, projected injection profiles and/or spatiotemporal estimates of pressure/stress, and to tune various model parameters. ORION will then provide a spatial and temporal ensemble forecast of seismicity defined as the probability of exceedance of a given earthquake magnitude over a forecast period. Additionally, ORION will provide probability distribution of the statistically derived maximum possible earthquake magnitude that may be expected. Finally, ORION will provide suggested operational management strategies (e.g. reduce injection volumes at specific wells) based on the level of hazard.

NRAP,NRAP Tools,ORION,Operational Forcasting of In↗

Real-Time Bayesian Inference at Extreme Scale: A Digital Twin for Tsunami Early Warning Applied to the Cascadia Subduction Zone

We present a Bayesian inversion-based digital twin that employs acoustic pressure data from seafloor sensors, along with 3D coupled acoustic–gravity wave equations, to infer earthquake-induced spatiotemporal seafloor motion in real time and forecast tsunami propagation toward coastlines for early warning with quantified uncertainties. Our target is the Cascadia subduction zone, with one billion parameters. Computing the posterior mean alone would require 50 years on a 512 GPU machine. Instead, exploiting the shift invariance of the parameter-to-observable map and devising novel parallel algorithms, we induce a fast offline–online decomposition. The offline component requires just one adjoint wave propagation per sensor; using MFEM, we scale this part of the computation to the full El Capitan system (43,520 GPUs) with 92% weak parallel efficiency. Moreover, given real-time data, the online component exactly solves the Bayesian inverse and forecasting problems in 0.2 seconds on a modest GPU system, a ten-billion-fold speedup.

97 MATHEMATICS AND COMPUTING↗

3-D Simulations of earthquakes rupture jumps: 1. Homogeneous pre-stress conditions

SUMMARY Observational and modelling studies indicate that earthquake ruptures can jump between fault sections as large as ∼3 and ∼5 km for compressional and extensional offsets, respectively. Here, we compare characteristics of the rupture jump process on parallel but offset fault sections from traditional 3-D dynamic rupture simulations governed by slip weakening friction using the finite element code, FaultMod, to those from quasi-dynamic simulations governed by rate- and state-dependent friction (rate-state friction) using the code RSQSim. These simulations use spatially uniform initial stresses. For a variety of measures the rupture renucleation position on the offset fault, the rate-state friction and slip weakening friction models produce very similar results. The principal difference is the additional occurrence of delayed rupture jumps that arise from the time- and stress-dependent nucleation that is characteristic of rate-state friction. For immediate rupture jumps, models with slip weakening friction span greater offsets than those with rate-state friction. However the jump distances are nearly identical when delayed rupture jumps are included in the comparisons. We propose that delayed rupture jumps are the likely mechanism for adjacent large-earthquake pairs and clusters. Based on the similarity of renucleation positions with both dynamic and quasi-dynamic models, we conclude that the renucleation positions for rupture initiation on the receiver fault (separated by less than ∼3 km from the source fault) are primarily controlled by static stress changes induced by slip on the initiating fault. However, in light of the slightly greater maximum jump distances (>3 km) seen with the dynamic slip weakening friction model, dynamic stress changes from seismic waves play an increasingly important role as offset distances increase.

, RSQSim↗

Attention Network Forecasts Time-to-Failure in Laboratory Shear Experiments

Rocks under stress deform by creep mechanisms that include formation and slip on small-scale internal cracks. Intragranular cracks and slip along grain contacts release energy as elastic waves termed acoustic emissions (AE). AEs are thought to contain predictive information that can be used for fault failure forecasting. Here, we present a method using unsupervised classification and an attention network to forecast labquakes using AE waveform features. Our data were generated in a laboratory setting using a biaxial shearing device with granular fault gouge intended to mimic the conditions of tectonic faults. Here, we analyzed the temporal evolution of AEs generated throughout several hundred laboratory earthquake cycles. We used a Conscience Self-Organizing Map (CSOM) to perform topologically ordered vector quantization based on waveform properties. The resulting map was used to interactively cluster AEs. We examined the clusters over time to identify those with predictive ability. Finally, we used a variety of LSTM and attention-based networks to test the predictive power of the AE clusters. By tracking cumulative waveform features over the seismic cycle, the network is able to forecast the time-to-failure (TTF) of lab earthquakes. Our results show that analyzing the data to isolate predictive signals and using a more sophisticated network architecture are key to robustly forecasting labquakes. In the future, this method could be applied on tectonic faults to monitor earthquakes and augment early warning systems.

58 GEOSCIENCES↗

Coastal Land Change Due to Tectonic Processes and Implications for Relative Sea-Level Rise in the Samoan Islands

Of the major coastal land change mechanisms responsible for relative sea-level change, tectonic subsidence is generally quoted as ranging from < mm/yr to 1 cm/yr. However, we documented coseismic and ongoing post-earthquake surface displacements from continuous GPS and tide gauge/altimetry data that indicated rapid subsidence on two of the major Samoan Islands of 12 - 20 cm during and following the 8.1 2009 Tonga-Samoa earthquake. Earlier results and our modeling of GRACE-derived gravimetric data provided a preliminary forecast of future relative sea-level rise through rapid land subsidence [Han et al., 2019]. Of course these numerical forecasts of time-dependent deformation are only as good as our input observations and our assumed rheological models. As part of our current NASA Earth Surface and Interior study, we are obtaining a wider range of data to constrain and test alternate models of ongoing postseismic deformation across American Samoa and Upolu, Samoa: (1) times series of altimetry plus tide gauge data processed to complement the cGPS data available to provide high-temporal resolution, point measurements of uplift/subsidence, (2) InSAR derived observations of surface deformation across the highly vegetated Samoan Islands, (3) evaluating and using NASA satellite lidar data (ICESat-I & ICESat-II, GEDI) for fusion with multi-source topographic data sets and for estimating topographic change on the decadal time scale. We are evaluating and using these new observations to better understand and separate out local, island-wide, and multi-island subsidence patterns and to evaluate the high impact of rising sea-level in a tectonically active region.

Jeanne Sauber-Rosenberg↗

DeepShake: Shaking Intensity Prediction Using Deep Spatiotemporal RNNs for Earthquake Early Warning

Abstract We propose a deep spatiotemporal recurrent neural network, DeepShake, to project future shaking intensity directly from current ground-motion observations. DeepShake is a network-based forecasting model, able to predict future shaking intensity at all stations within a network given previously measured ground shaking. The model is not given any a priori knowledge of station locations; instead, it learns wave propagation amplitudes and delays solely from training data. We developed DeepShake with the 35,679 earthquakes from the 2019 Ridgecrest sequence. Tasked with alerting for modified Mercalli intensity (MMI) IV+ shaking on 3568 validation earthquakes at least 5 s in advance, DeepShake achieves an equal error rate of 11.4%. For the Mw 7.1 earthquake that hit Ridgecrest on 5 July 2019, DeepShake was able to provide targeted alerts to all stations inside the network 5 s prior to the arrival of MMI IV+ waveforms. DeepShake demonstrates that deep spatiotemporal neural networks can effectively provide one-step earthquake early warning with reasonable accuracy and latency.

Geochemistry & Geophysics↗

Volcanism in Eastern Africa

In 1891, the Virunga Mountains of Eastern Zaire were first acknowledged as volcanoes, and since then, the Virunga Mountain chain has demonstrated its potentially violent volcanic nature. The Virunga Mountains lie across the Eastern African Rift in an E-W direction located north of Lake Kivu. Mt. Nyamuragira and Mt. Nyiragongo present the most hazard of the eight mountains making up Virunga volcanic field, with the most recent activity during the 1970-90's. In 1977, after almost eighty years of moderate activity and periods of quiescence, Mt. Nyamuragira became highly active with lava flows that extruded from fissures on flanks circumscribing the volcano. The flows destroyed vast areas of vegetation and Zairian National Park areas, but no casualties were reported. Mt. Nyiragongo exhibited the same type volcanic activity, in association with regional tectonics that effected Mt. Nyamuragira, with variations of lava lake levels, lava fountains, and lava flows that resided in Lake Kivu. Mt. Nyiragongo, recently named a Decade volcano, presents both a direct and an indirect hazard to the inhabitants and properties located near the volcano. The Virunga volcanoes pose four major threats: volcanic eruptions, lava flows, toxic gas emission (CH4 and CO2), and earthquakes. Thus, the volcanoes of the Eastern African volcanic field emanate harm to the surrounding area by the forecast of volcanic eruptions. During the JSC Summer Fellowship program, we will acquire and collate remote sensing, photographic (Space Shuttle images), topographic and field data. In addition, maps of the extent and morphology(ies) of the features will be constructed using digital image information. The database generated will serve to create a Geographic Information System for easy access of information of the Eastem African volcanic field. The analysis of volcanism in Eastern Africa will permit a comparison for those areas from which we have field data. Results from this summer's work will permit further study and monitoring of the volcanic activity in the area. This is of concern due to the large numbers of refugees fleeing into Zaire where they are being positioned at the base of Mt. Nyiragongo. The refugees located at the base of the volcano are in direct hazard of suffocation by gas emission and destruction by lava flow. The results from this summer study will be used to secure future funding to enable continuation of this project.

Cauthen, Clay↗

Annual review of earth observations from space

An overview is given of the present state of satellites making observations of the earth. Satellite systems discussed include the NOAA series of Synchronous Meteorological Satellites (SMS) and sun-synchronous satellites, the two LANDSATS (formerly called ERTS), and the NIMBUS series. Examples are presented of the types of observations being made as well as their purposes. These include observations of synoptic and mesoscale atmospheric processes for daily weather forecasting, global atmospheric processes for long-range weather forecasting, planetary radiation budget and ocean circulation for monitoring climatic trends, earth dynamics and tectonic structure for mineral exploration and assessing earthquake hazards, atmospheric composition and water quality for environment monitoring, and thematic mapping for monitoring land use, managing crops and water resources, and assessing environmental impacts.

Nordberg, W.↗