Search NASA⌕ Search

SEARCH · Search NASA

Results for “Seismic Hazard”

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 109 records · Page 6

Compilation of a Comprehensive Earthquake Catalog and Relocations in the Caucasus Region

Instrumental seismic monitoring has a long history in the Caucasus and started in 1899 when the first seismograph was installed in Tbilisi, Georgia. Much of the analog paper records from this time period are preserved in the Tbilisi archives because Georgia served as the regional data center. In the 1990s, due to the collapse of the Soviet Union and the political turmoil in the region, the analog networks and the communication between the newly formed national networks deteriorated. In Georgia, for the next 13 yr, the seismic network coverage was poor until the 2002 Tbilisi earthquake. Following this earthquake, the first permanent digital seismic station in Georgia was established in Tbilisi in 2003. The digital era progressively improved the ability to collect and archive data and today more than a hundred broadband seismic stations (including temporary arrays) are operating in the southern Caucasus. Until recently, the region lacked a coordinated effort to catalog all analog and digital era data collected by different countries into a single repository. As a result of collaboration between Lawrence Livermore National Laboratory, the Ilia State University, and the Republican Seismic Survey Center of Azerbaijan, a comprehensive earthquake catalog was compiled for the Caucasus and neighboring areas as part of a broader probabilistic seismic hazard assessment project. Here this project digitized Soviet-era paper bulletins, compiled a unified earthquake catalog from regional bulletins, developed 1D reference velocity model, and used it to relocate the events. The final catalog contains 16,963 events with magnitudes 3.7 and above, bringing together all the available data sets in the Caucasus region from 1900 to 2015, significantly improving locations, and generating the most complete earthquake catalog in the region, temporally and geographically.

58 GEOSCIENCES↗

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↗

Report on Depth-Dependent Q from Frequency-Dependent Lg Q

As a phase that travels through the crustal waveguide, Lg samples the velocity and attenuation structure of the crust. We have developed methods (Pasyanos et al., 2009a; Pasyanos et al., 2009b) to estimate Qp and Qs of the crust and upper mantle from the amplitudes of regional phases, and have applied it to a number of regions around the world, and over a broad frequency band (0.5 - 10 Hz). In prior attenuation modeling, Lg phase amplitudes are able to effectively distinguish the high Q (low attenuation) of old stable cratonic crust from the low Q (high attenuation) of more recently active tectonic regions. For a variety of purposes (e.g. waveform modeling, seismic hazard, etc.), it is useful to have more information on the depth-dependence of the Q structure in the crust. One obvious candidate for doing this is using the frequency-dependence information of Lg Q, in the manner of surface waves. Unfortunately, however, Lg samples the whole of the crust without much differentiation with respect to frequency that is characteristic of fundamental mode surface waves. Lg is generally thought to be the sum of supercritically reflected S-waves trapped in the crustal waveguide (e.g. Xie and Lay, 1994). The phase can also be modeled as a superposition of surface waves. At lower frequencies, it is primarily composed of fundamental and low-order surface waves, while at higher frequencies, it is comprised purely of many high-order surface waves. The low-order surface waves that comprise Lg in low frequency bands (< 1 Hz) are preferentially sensitive to shallow earth structure, such as sedimentary basins and upper crust, while high-order surface waves have a more homogenized sensitivity to the whole crust. Furthermore, at the lower frequencies, we might expect to see some differentiation in the sensitivity between early Lg (herein referred to as Lg1) and late Lg (referred to as Lg2) where the latter is comprised of lower-order surface waves and, hence, shallower structure. Therefore, we might expect to see the signature of shallow features, primarily the lower Q of sediments as compared to crystalline crust, at lower frequencies and in later arriving Lg amplitudes.

58 GEOSCIENCES↗

Research Introduction [Slides]

In the topic of Reverse Time Imaging, we proposed a new IC to reduce computation cost but reserve image resolution for distributed sensor networks. For Induced Seismicity in Oklahoma, we analyzed fault stress state analysis at state scale, and we applied machine learning techniques to polarity picking and seismicity rate forecasting. The results provide better understanding of fault properties, stress field, and the relationships among fault, stress state, injections, and potential seismic hazards. Lastly, for Microseismic Monitoring, we detected and located 770 low-frequency events (5-50 Hz): (1) Shallow events are highly clustered, consistent with the pathway from injection well 13-10A to monitoring well; moment tensor inversion shows dominant tensile cracking; (2) Deep events are scattered and show migration pattern to the basement; moment tensors show that most events are shear cracks.

58 GEOSCIENCES↗

5-2428: Coupled Investigation of Fracture Permeability Impact on Reservoir Stress and Seismic Slip Behavior [Slides]

Our goal is to develop, apply and validate a holistic thermal, hydrologic, mechanical, and chemical (THMC) workflow that also includes evaluation of induced seismic slip in EGS reservoirs. We will integrate experimental and modelling approaches to reduce parameter uncertainty and better integrate friction constitutive laws for faults to improve our ability to predict and mitigate seismic hazard at Utah FORGE and future EGS sites. We propose a novel approach that incorporates 3D physics-based Earthquake simulations in THMC models, using lab measurements of friction parameters, herein referred to as “THMC+E” models. These simulations will enable exploration of various circumstances that may hinder EGS success and develop mitigation strategies.

58 GEOSCIENCES↗

Simulation of Physics-Based 0-10Hz Strong Motion Using High Performance Computing Supporting Refinements to Regional Ground Motion Models for the Central Eastern US

In collaboration with the U.S. Nuclear Regulatory Commission (NRC) the LLNL has developed a computationally efficient simulation platform designed to perform physics-based ground motion simulations for crustal earthquakes in the Stable Continental Regions of Central and Eastern US (CEUS), using high-performance computing. The main objective of the earthquake simulations was to use synthetic ground motion to provide constrains to refinements of existing ergodic Ground Motion Models (GMMs), for large magnitude earthquakes and near-fault distances, for which these models are less reliable. Physics-based broadband (0-10Hz) 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. In our simulations we used a 3D regional velocity model that was based on Saikia’s 1D velocity model (1994). In simulations performed during the first stage of this project the Saikia’s velocity model demonstrated better performance in modelling high frequency regional wave propagation for the CEUS region recorded during the Mw5.0 November 7, 2016, Cushing Oklahoma (Taylor et al., 2017), and Mw5.8 September 3, 2016, Pawnee Oklahoma earthquakes. The proposed regional 3D model includes random perturbations to the 1D background model using the stochastic scheme of Pitarka and Mellors (2021). In addition, validation analysis of the rupture generator and regional wave propagation models, using comparisons with different GMMs for Mw6.5 and Mw7.0 scenario earthquakes in the CEUS region resulted in a very good match between the simulated and empirical ground motion models. For the purposes of seismic hazard assessment at the existing and planned nuclear power plants, NRC is interested in studies aimed at improving the current ground motion models (GMM) for both Stable Continental Regions (SCR) in the Central and Eastern US and Active Crustal Regions (ACR) in the Western US. Due to lack of recorded data, these improvements require synthetic data for short fault distances and large magnitude earthquakes for which the existing recorded data is not enough to uniquely constrain the GMMs. The need for simulations and strong motion data is especially critical for the CEUS region where we do not have recorded data from potentially large damaging earthquakes with moment magnitudes 6.0 and higher. In this the project, we focused on 10Hz simulations of Mw7.0 scenario earthquakes with strike slip and thrust faulting mechanisms. We used more than 50 Mw7.0 earthquake rupture scenarios to investigate the ground motion uncertainty due to unknown earthquake rupture parameters, in particular, the slip distribution, rupture velocity, and faulting mechanism, and their implication on ground motion amplification due to forward rupture directivity effects.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Activating a Natural Fault Zone in the Swiss Alps

One major hurdle for understanding earthquake mechanics are observational limitations. Important phenomena like strain localisation, fault dilation, and fault healing are readily studied in rock mechanical laboratory experiments and with numerical models. At the scale of natural earthquakes, however, these phenomena are often unresolvable, even by state-of-the-art observatories. To overcome this limitation, we are currently building the Earthquake Physics Testbed at the Bedretto Underground Laboratory for Geosciences and Geoenergies (BedrettoLab), an experimental testbed where we can activate an extensively instrumented natural fault zone via hydraulic stimulation. The goal of the Fault Activation and Earthquake Rupture (FEAR) project is to induce earthquakes of up to Mw~1.0 on this exceptionally well characterised and instrumented fault zone. Here we summarize the main scientific goals and current FEAR project status, and present first results from conducted experiments. We discuss how this large-scale experimental approach may allow us to tackle both fundamental science as well as practical questions on earthquake physics, induced seismicity and seismic hazard.

Meier, Men-Andrin↗

Fracture Network Localization Preceding Catastrophic Failure in Triaxial Compression Experiments on Rocks

We quantify the spatial distribution of fracture networks throughout six in situ X-ray tomography triaxial compression experiments on crystalline rocks at confining stresses of 5–35 MPa in order to quantify how fracture development controls the final macroscopic failure of the rock, a process analogous to those that control geohazards such as earthquakes and landslides. Tracking the proportion of the cumulative volume of fractures with volumes >90th percentile to the total fracture volume, ∑ v 90 / v t o t indicates that the fracture networks tend to increase in localization toward these largest fractures for up to 80% of the applied differential stress. The evolution of this metric also matches the evolution of the Gini coefficient, which measures the deviation of a population from uniformity. These results are consistent with observations of localizing low magnitude seismicity before large earthquakes in southern California. In both this analysis and the present work, phases of delocalization interrupt the general increase in localization preceding catastrophic failure, indicating that delocalization does not necessarily indicate a reduction of seismic hazard. However, the proportion of the maximum fracture volume to the total fracture volume does not increase monotonically. Experiments with higher confining stress tend to experience greater localization. To further quantify localization, we compare the geometry of the largest fractures, with volumes >90th percentile, to the best fit plane through these fractures immediately preceding failure. The r 2 scores and the mean distance of the fractures to the plane indicate greater localization in monzonite than in granite. The smaller mean mineral diameter and lower confining stress in the granite experiments may contribute to this result. Tracking these various metrics of localization reveals a close association between macroscopic yielding and the acceleration of fracture network localization. Near yielding, ∑ v 90 / v t o t and the Gini coefficient increase while the mean distance to the final failure plane decreases. Macroscopic yielding thus occurs when the rate of fracture network localization increases.

58 GEOSCIENCES↗

Quaternary slip rates and most recent surface rupture of the Bitterroot fault, western Montana

The Bitterroot fault is a ~100-km-long Quaternary active normal fault that bounds the eastern margin of the north–south-trending Bitterroot Mountains and accommodates extension near the Intermountain Seismic Belt. Here, new detailed mapping using high-resolution topographic data derived from light detection and ranging (lidar) along the southern Bitterroot Range documents multiple generations of fault scarps in Holocene–Pleistocene deposits with vertical off sets that increase in magnitude with age. Fault mapping indicates a complex fault geometry characterized by an en echelon pattern of discontinuous segments of 41–78° east-dipping normal faults that appear to cut the older Eocene detachment fault, and locally 70–88° west-dipping antithetic normal faults. 10 Be cosmogenic exposure dating provides in situ age control for 32 surface boulders (>1 m) sampled in glacial deposits. Near Como Dam, two Pinedale-age glacial moraine sequences yield peak age distributions of 15.0 ± 0.4 ka and 16.4 ± 0.6 ka as apparent exposure ages (e = 0), and 15.4 ± 0.4 ka and 16.8 ± 0.6 ka based on the maximum allowed boulder surface erosion rate (e = 2 mm/ka). Vertical separation of 3.5 ± 0.2 m across the ~16–17 ka glacial moraine off set by the Bitterroot fault scarp yields a fault slip rate of 0.2–0.3 mm/yr. Glacial Lake Missoula highstand shorelines, inset into the ~15 ka glacial moraine and vertically off set 4.6 ± 1.6 m by an antithetic strand of the Bitterroot fault, yield fault slip rates of 0.2–0.5 mm/yr that overlap with fault slip rates on the main strand near Lake Como (0.2–0.3 mm/yr). At the Ward Creek fan located ~15 km to the north of Lake Como, two glacial debris fan sequences yield peak age distributions of 16.6 ± 0.4 ka and 62.8 ± 1.7 ka (e = 0), and 17.0 ± 0.4 ka and 69.9 ± 2.2 ka (e = 2 mm/ka). Vertical separations of 2.4 ± 0.3 m and 4.5 ± 0.2 m on the ~17 ka and ~63–70 ka fan surfaces off set by the Bitterroot fault yield fault slip rates of 0.2–0.3 mm/yr and 0.1 mm/yr, respectively. Our results indicate broadly consistent fault slip rates for the main fault segments at Lake Como (0.2–0.3 mm/yr) and the Ward Creek fan (0.1–0.3 mm/yr) with an along-strike range of 0.1–0.3 mm/yr for the southern Bitterroot fault. Fault scaling relations and evidence of multiple late Quaternary fault surface ruptures suggest the Bitterroot fault could produce a Mw ~7.2 earthquake. Structural model constraints and our slip rate results indicate both high-angle or low-angle fault geometries are possible at depth. A seismogenic low-angle fault model could generate a larger earthquake of Mw >7.2. Earthquake history is unknown for the Bitterroot fault, but fault scarps in young glacial deposits demonstrate its seismogenic potential. Data from this study suggest seismic hazards from the Bitterroot fault may pose a high level of risk to the Missoula metropolitan area, the State’s second most populous region, and major infrastructures across the Missoula and Bitterroot Valleys.

58 GEOSCIENCES↗

Subsurface stress criticality associated with fluid injection and determined using machine learning

Machine-learning methods and apparatus are disclosed to determine critical state or other parameters related to fluid-driven failure of a terrestrial locale impacted by anthropogenic activities such as hydraulic fracturing, hydrocarbon extraction, wastewater disposal, or geothermal harvesting. Acoustic emission, seismic waves, or other detectable indicators of microscopic processes are sensed. A classifier is trained using time series of microscopic data along with corresponding data of critical state or failure events. In disclosed examples, random forests and artificial neural networks are used, and grid-search or EGO procedures are used for hyperparameter tuning. Once trained, the classifier can be applied to live data from a fluid injection locale in order to assess a frictional state, assess seismic hazard, assess permeability, make predictions regarding a future fluid-driven failure event, or drive engineering solutions for mitigation or remediation. Variations are disclosed.

Johnson, Paul Allan↗

Deep Isolation-Development of the Safety Case for Disposal of Radioactive Wastes in Horizontal Boreholes - 20028

Deep Isolation has developed a safe, secure, and permanent geological disposal method for high-level waste, including spent nuclear fuel as well as sealed sources and other highly radioactive materials. The method leverages well established directional drilling technology to create horizontal repositories deep underground. The combination of great depth (1-3 km) and the ability to precisely position repositories in a horizontal orientation provides access to geologic strata that are inaccessible for typical mined and deep vertical borehole repositories. A number of potential direct and indirect safety benefits accrue including:1) greater depth below regional freshwater aquifers; 2) increased flexibility and more geologic options for siting; 3) access to deep formations with persistent and sustained reducing conditions. In addition there are a number of safety elements related specifically to the horizontal repository geometry including: 1) passive direction of thermally driven fluid and radionuclide movement away from the vertical access hole and toward the 'dead end' portion of the repository; 2) mitigation of seismic hazards by orienting repositories parallel to local and regional fault structures. In this paper we explore and discuss some of the key technologic, geologic and hydrologic elements that support the deep horizontal borehole safety case. The stalemate seen across the globe on the disposal of spent nuclear fuel and high-level waste can be broken. Deep Isolation offers a novel option for safe, secure, and permanent deep geological disposal of nuclear waste that can be developed as centralized repositories or adapted to smaller regional or site-specific repositories located near waste sources. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Advanced High-Performance Computational Modeling of the Seismic Response of High-Hazard and/or Nuclear Facilities and Critical Infrastructure at the NNSS

New methods for predicting the amplitude and variability of ground shaking from earthquakes (and explosions) are needed for seismic hazard analysis for buildings, nuclear power plants, and critical infrastructure at the NNSS. We are comparing existing 1-D and new 3-D geophysical methods for estimating the shear-wave velocity structure in the upper 30 meters of the ground surface (Vs30), which plays a major role in ground motion amplification and seismic response of buildings. We evaluate the performance of these methodologies at the U1a Complex at the NNSS and develop simple 1-D and high-resolution 3-D Vs30 models. We then emplace these high-resolution models into a background seismic velocity model. We will collaborate with Lawrence Livermore National Laboratory (LLNL) to conduct numerical modeling of the ground shaking at the NNSS using their high-performance computing technology and state-of-the-art ground motion simulation methodology. The primary work that was completed in FY 2019 was to acquire the seismic systems and familiarize staff at the NNSS with their use. We also worked on developing a collection plan with the Device Assembly Facility (DAF) at the NNSS, but due to time constraints and other ongoing projects at the DAF, we had to use U1a as a backup. We were able to coordinate the seismic survey, and we will complete the collection of seismic data in FY 2020. Additionally, during FY 2019, we completed the geologic framework model (GFM) for the U1a Complex and modeled the Yucca fault. LLNL worked with us through FY 2019 to prepare the data files for their modeling software and tested the software for reliability. In FY 2020 we will develop the end-to-end capability so that any facility could easily be modeled and the expected shaking from a local earthquake understood. The work in FY 2020 will include building fault models from the GFM and finalizing the velocity model analysis. The final simulations will be run for multiple rupture models, and final assessments will demonstrate the seismic hazard at the U1a Complex.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

2024.03.13 Bushing Subcommittee presentation on Transformer Bushing Seismic Resilience: Design, Hazard, & Prevention Analysis [Slides]

Presents results from power grid simulations on synthetic grids that estimate the impact from power transformer failures due to earthquakes. In the best-case scenario, there is an 82.9% likelihood of load loss greater or equal to 51 MW. Methodologies are described as well as path forward to improving transformer resilience and, consequently, the resilience of the bulk power grid.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Seismic Loading for Short-Term Duration Exposures and Temporary Structures

The purpose of this document is to determine the appropriate seismic demand forces associated with a short-term and temporary seismic exposure for protecting structures, systems, and components (SSCs) from hazards associated with Seismic Design Category 1 (SDC-1), SDC-2, and SDC-3 seismic events based on data provided in the Probabilistic Seismic Hazards Update for LLNL prepared by Fugro Consultants, Inc. in accordance with DOE-STD-1020.

42 ENGINEERING↗

Prestack Waveform Inversion at GC 955: Trials and sensitivity of PWI to high-resolution seismic data (Data Report)

Understanding the spatial extent and variability of hydrate deposits has important implications for potential economic production, climate change, and assessing natural hazards risks. Seismic reflection techniques are capable of determining the extent of gas hydrate deposits and can shed light on hydrate concentrations and properties. Using high resolution prestack time migrated seismic data collected during a 2013 US Geological Survey expedition and prestack waveform inversion (PWI) techniques, we produce highly resolved velocity models and compare them to co-located well logs. Coupling our PWI results with velocity-porosity relationships and nearby well control, we compare hydrate properties at site 955 in the Green Canyon and site 313 at Walker Ridge.

03 NATURAL GAS↗

Stress-Drop and Source Scaling of the 2019 Ridgecrest, California, Earthquake Sequence

Stress drop, while difficult to measure reliably and at scale, is a key source parameter for understanding the earthquake rupture process and its relationship to strong ground motion. Here, we use a P-wave spectral decomposition approach, designed for large and densely sampled datasets, to measure earthquake stress drop in the region surrounding the 2019 Ridgecrest, California, earthquake sequence. With more than 11,000 measurements of earthquake stress drop in the 20-yr time period from 2000 through 2019, this dataset provides an opportunity to understand how coseismic stress changes and how other geophysical factors relate to the distribution of stress drop and its evolution in space and time. We observe a mild but persistent deviation from self-similar scaling, with larger events having systematically higher stress drops, though this trend depends on the assumption of an omega-square source spectral model. Earthquake stress drop increases with hypocentral depth in this study region, and the Ridgecrest aftershocks tend to have higher stress drops than the pre-event seismicity. This is in part due to their deeper hypocenters. Coherent spatial patterns of stress drop in the aftershock sequence correlate with the slip distribution of the M 7.1 mainshock, whose northwest rupture tip terminated in a long-lived zone of enervated stress drop. Although physical interpretation of these results is complicated by the trade-offs between the timing, depth, and location of these earthquakes, the observations provide new insight into the physics of the earthquake source in an area of renewed seismic activity in southern California.

58 GEOSCIENCES↗

High-Quality Revision of the Israeli Seismic Bulletin

Seismic bulletins, with trustworthy phase picks, origin times, and source locations are key for regional seismic studies, such as travel-time (TT) tomography, attenuation tomography, and anisotropy studies. To lay the groundwork for such studies in Israel, we revised the seismic bulletin of Israel and the surrounding area and obtained a trustworthy TT data set. From the earthquake and explosion bulletins of the Geophysical Institute of Israel, we compiled a starting data set of about 123,000 earthquakes and explosions that occurred during the past 40 yr. After screening out the poorly recorded events, we were left with a data set of ~38,000 well-recorded events. We then revised the remaining data set in two consecutive steps. In the first, we reviewed and updated station metadata, including changes in station metadata parameters over time. In the second step, we jointly relocated a list of selected seismic events, using the Bayesian hierarchical location software package (BayesLoc) of Myers et al. (2007) that performs joint relocation of multiple events. We observed striking dissimilarities between the spatial distributions of the newly relocated catalog and the initial locations. Although the depth distribution of the starting catalog is trimodal with peaks at 0, 5, and 10 km, the distribution in this study is unimodal, with a broad peak between 7.5 and 12.5 km. By differencing the observed arrival times and the origin times obtained through relocation with BayesLoc, we obtained a revised TT database that consists of 261,336 Pg, 132,876 Pn, 114,816 Sg, and 60,394 Sn arrivals, from a set of 30,458 jointly relocated seismic sources. In this work, we compared prerevision and postrevision TTs as a function of epicentral distance and concluded that the revised data set contains far fewer outliers and inconsistencies than the original data set. The revised TT data set may be used for seismic studies, such as TT tomography, attenuation tomography, and anisotropy studies.

58 GEOSCIENCES↗

Constraining maximum event magnitude during injection-triggered seismicity

Abstract Understanding mechanisms controlling fluid injection-triggered seismicity is key in defining strategies to ameliorate it. Recent triggered events (e.g. Pohang, Mw 5.5) have exceeded predictions of average energy release by a factor of >1000x, necessitating robust methodologies to both define critical antecedent conditions and to thereby constrain anticipated event size. We define maximum event magnitudes resulting from triggering as a function of pre-existing critical stresses and fluid injection volume. Fluid injection experiments on prestressed laboratory faults confirm these estimates of triggered moment magnitudes for varied boundary conditions and injection rates. In addition, observed ratios of shear slip to dilation rates on individual faults signal triggering and may serve as a measurable proxy for impending rupture. This new framework provides a robust method of constraining maximum event size for preloaded faults and unifies prior laboratory and field observations that span sixteen decades in injection volume and four decades in length scale.

58 GEOSCIENCES↗