Geology topics

Elizabeth S. Cochran

Publications and source records attributed to Elizabeth S. Cochran.

At least 55 records · Page 3Linked to original sources

USGS2021

Evidence for latent crustal fluid injection transients in southern California from long-duration earthquake swarms

Earthquake swarms are manifestations of aseismic driving processes deep in the crust. We examine the spatiotemporal distribution of aseismic processes in Southern California using a 12-years catalog of swarms derived with deep learning algorithms. In a core portion of the plate boundary region, which is not associated with elevated heat flow, we identify 92 long-duration swarms ranging from 6 months to 7 years that constitute 26.4% of the total seismicity. We find that 53% of the swarms exhibit ultra-slow diffusive patterns with propagating backfronts, consistent with expectations for natural fluid injection processes. The chronology of the swarms indicates that the aseismic driving processes were active at all times during 2008–2020. The observations challenge common views about the nature of swarms, which would characterize any one of these sequences as anomalous. The regional prevalence of these sequences suggests that transient fluid injection processes play a key role in crustal fluid transport.

California↗

USGS2021

A framework for evaluating earthquake early warning for an infrastructure network: An idealized case study of a northern California rail system

Earthquake early warning (EEW) systems provide a few to tens of seconds of warning before shaking hits a site. Despite the recent rapid developments of EEW systems around the world, the optimal alert response strategy and the practical benefit of using EEW are still open-ended questions, especially in areas where EEW systems are new or have not yet been deployed. Here, we use a case study of a rail system in California’s San Francisco Bay Area to explore potential uses of EEW for rail systems. Rail systems are of particular interest not only because they are important lifeline infrastructure and a common application for EEW around the world, but also because their geographically broad yet networked infrastructure makes them almost uniquely well suited for utilizing EEW. While the most obvious potential benefit of EEW to the railway is to prevent derailments by stopping trains before the arrival of shaking, the lead time for warnings is usually not long enough to significantly reduce a train’s speed. In reality, EEW’s greatest impact is preventing derailment by alerting trains to slow down or stop before they encounter damaged track. We perform cost-benefit analyses of different decision-making strategies for several EEW system designs to find an optimal alerting strategy. On-site EEW provides better outcomes than source-parameter-based EEW when warning at a threshold of 120 gal (the level of shaking at which damage might occur) regardless of false alarm tolerance. A source-parameter-based EEW system with a lower alerting threshold (e.g., 40 gal) can reduce the exposure to potentially damaged track compared to an on-site system alerting at 120 gal, but a lower alerting threshold comes at the cost of additional precautionary system stops. The optimal EEW approach for rail systems depends strongly on the ratio of the cost of stopping the system unnecessarily to the potential loss from traversing damaged tracks.

California↗

USGS2021

Fluid-earthquake and earthquake-earthquake interactions in southern Kansas, USA

An increase in injection activity associated with energy production in southern Kansas starting in 2013 has been linked to the occurrence of more than 130,000 earthquakes (M −1.5 to 4.9) between 2014 and 2017. Studies suggest that the dramatic increase in seismicity rate is related to wastewater injection into the highly permeable Arbuckle formation. Most of the seismicity is located in the underlying crystalline basement, for which hydrological properties and specific fault geometries are unknown. Additionally, some earthquake clusters occurred relatively far (tens of kilometers) from the main injection wells. Therefore, the effect of pore pressure diffusion may be insufficient to explain the relation between the volume of injected fluids and the spatiotemporal evolution of seismicity. Combining physical models (static stress and poroelasticity) and a statistical cluster analysis applied to a high-resolution relocated catalog, we analyze the evolution of seismicity in southern Kansas. We find that pore pressure changes (Δ p ) and Coulomb stress changes (ΔCFS) due to fluid diffusion smaller than 0.1 MPa are enough to initiate seismic sequences, which then evolve depending on their distance from the major injection wells. However, we find that earthquake sequences have different seismogenic responses to Δ p and ΔCFS in terms of triggering threshold. In regions located close to disposal wells (Harper area) our cluster analysis suggests that both earthquake interactions and fluid diffusion control the evolution of seismicity. On the other hand, at greater distances (Milan area), where clustering behavior suggests greater earthquake interactions, we find that coseismic ΔCFS are larger than Δ p .

Kansas, Oklahoma↗

USGS2021

Characterizing stress orientations in southern Kansas

Induced seismicity predominantly occurs along faults that are optimally oriented to the local principal compressive stress direction, and the characterization of these stress orientations is an important component of understanding seismic hazards. The seismicity rate in southern Kansas rapidly increased in 2013 primarily due to the disposal of large volumes of wastewater into the Arbuckle Group. Previously, local stress orientations in this area were poorly constrained, which limited our understanding of the complex faulting and diverse earthquake mechanisms in this region. We use shear‐wave splitting and focal mechanism inversion techniques to create multiple, independent estimates of maximum horizontal stress directions ( ⁠ S H max ">

Kansas↗

USGS2021

Characteristics of frequent dynamic triggering of microearthquakes in Southern California

Dynamic triggering of earthquakes has been reported at various fault systems. The triggered earthquakes are thought to be caused either directly by dynamic stress changes due to the passing seismic waves, or indirectly by other nonlinear processes that are initiated by the passing waves. Distinguishing these physical mechanisms is difficult because of the general lack of high‐resolution earthquake catalogs and robust means to quantitatively evaluate triggering responses, particularly, delayed responses. Here we use the high‐resolution Quake Template Matching catalog in Southern California to systematically evaluate teleseismic dynamic triggering patterns in the San Jacinto Fault Zone and the Salton Sea Geothermal Field from 2008 to 2017. We develop a new statistical approach to identify triggered cases, finding that approximately 1 out of every 5 global M w ≥ 6 earthquakes dynamically trigger microearthquakes in Southern California. The triggering responses include both instantaneous and delayed triggering, showing a highly heterogeneous pattern and indicating possible evolving triggering thresholds. We do not observe a clear peak ground velocity triggering threshold that can differentiate triggering earthquakes from non‐triggering events, but there are subtle differences in the frequency content that may possibly differentiate the earthquakes. In contrast to the depth distribution of background seismicity, the identified triggered earthquakes tend to concentrate at the edges of the seismogenic zones. Although instantaneously triggered earthquakes are likely a result of dynamic Coulomb stress changes, the cases of delayed dynamic triggering are best explained by nonlinear triggering processes, including cyclic material fatigue, accelerated transient creep, and stochastic frictional heterogeneities.

California↗

USGS2021

Stress controls rupture extent and maximum magnitude of induced earthquakes

Seismic hazard forecasts of induced seismicity often require estimates of the maximum possible magnitude ( M max ). Empirical models suggest that maximum magnitudes, or expected number of earthquakes, are related to the volume of injected fluid. We perform a suite of 3D physics-based earthquake simulations with rate- and state-dependent friction, systematically varying the area of the pressurized region and the amplitude of the initial homogeneous or heterogeneous shear stress. Using the resulting catalog we explore the conditions that result in pressure-controlled versus runaway ruptures that extend outside the pressurized zone. We find that proposed empirical scaling laws correctly predict M max when shear stresses are further from failure (≤90% of maximum shear stress) and for high amplitude stress fields. Runaway ruptures are observed for higher initial shear stresses and smoother stress fields. In these cases, runaway ruptures occur early after the onset of injection and rarely preceded by foreshock activity.

Geophysical Research Letters↗

USGS2020

The induced Mw 5.0 March 2020 west Texas seismic sequence

On March 26, 2020, a M 5.0 earthquake occurred in the Delaware Basin, Texas, near the border between Reeves and Culberson Counties. This was the third largest earthquake recorded in Texas and the largest earthquake in the Central and Eastern United States since the three M 5.0–5.8 induced events in Oklahoma during 2016. Using multistation waveform template matching, we detect 3,940 earthquakes in the sequence with the first event in the area occurring in May 2018. The M 5.0 earthquake sequence occurred on a ENE (∼082°) normal fault dipping ∼37° toward the south. The earthquake caused 6 mm of oblique surface deformation, and geodetic slip inversion suggests slip was isolated above 6 km depth. We find that the sequence was most likely induced by nearby wastewater disposal operations, and seismicity rates in the region surrounding the M 5.0 will likely continue to increase in the future if disposal operations continue unaltered.

Texas↗

USGS2020

Apparent earthquake rupture predictability

To what extent can the future evolution of an ongoing earthquake rupture be predicted? This question of fundamental scientific and practical importance has recently been addressed by studies of teleseismic source time functions (STFs) but reaching contrasting conclusions. One study concludes that the initial portion of STFs is the same regardless of magnitude. Another study concludes that the rate at which earthquakes grow increases systematically and strongly with final event magnitudes. Here we show that the latter reported trend is caused by a selection bias towards events with unusually long durations, and by estimates of STF growth made when the STF is already decaying. If these invalid estimates are left out, the trend is no longer present, except during the first few seconds of the smallest events in the dataset, Mw5–6.5, for which the reliability of the STF amplitudes is questionable. Simple synthetic tests show that the observations are consistent with statistically indistinguishable growth of smaller and larger earthquakes. A much weaker trend is apparent among events of comparable duration, but we argue that its significance is not resolvable by the current data. Finally, we propose a nomenclature to facilitate further discussions of earthquake rupture predictability and determinism.

Geophysical Journal International↗

USGS2020

Shaking is almost always a surprise: The earthquakes that produce significant ground motion

Although small earthquakes are expected to produce weak shaking, ground motion is highly variable and there are outlier earthquakes that generate more shaking than expected—sometimes significantly more. We explore datasets of M 0.5–8.3 earthquakes to determine the relative impact of frequent, smaller-magnitude earthquakes that rarely produce strong ground motion, to rare, large earthquakes that always cause strong shaking. We find that the natural variability of ground motion, combined with the Gutenberg–Richter magnitude–frequency relationship, ensures that most occurrences of any ground motion come from earthquakes of smaller magnitude than expected, often > 2 magnitude units smaller. This holds even for very strong shaking ( > 20% g ), suggesting that M < 7 earthquakes could be a significant source of damage.

Seismological Research Letters↗

USGS2020

Using a large-n seismic array to explore the robustness of spectral estimations

Spectral analysis is widely used to estimate and refine earthquake source parameters such as source radius, seismic moment, and stress drop. This study aims to quantify the precision of the single spectra and empirical Green's function spectral ratio approach using the Large‐ n Seismic Survey in Oklahoma (LASSO) array. The dense station coverage in an area of local saltwater disposal offers a unique opportunity to observe and quantify radiation pattern effects and subsequent precision of spectral estimates of small earthquakes ( M < 3). The results suggest that the precision of source properties estimated from direct phase arrivals for arrays with less than 20 stations should be assumed to be not less than 30% and could be as high as 150% if less than five stations are used. Furthermore, we do not see clear evidence for, or against, a scaling of stress drop with magnitude of small earthquakes ( M < 3) as observed by other studies.

Oklahoma↗

USGS2020

VS30 and Dominant Site Frequency (⁠fd⁠) as Provisional Station ML Corrections (⁠dML⁠) in California

New seismic stations added to a regional seismic network cannot be used to calculate local magnitude ( ⁠ M L "> M L ⁠ ) until a revised regionwide amplitude decay function is developed. Each station must record a minimum number of local and regional earthquakes that meet specific amplitude requirements prior to recalibration of the amplitude decay function. Station component adjustments ( ⁠ d M L "> d M L ⁠ ; Uhrhammer et al. , 2011 ) are then calculated after inverting for a new regional amplitude decay function, constrained by the sum of d M L "> d M L ⁠ for long‐running stations. Therefore, there can be significant delay between when a new station starts contributing real‐time waveform packets and when data can be included in magnitude determinations. We propose the use of known estimates of seismic site conditions such as the time‐averaged shear‐wave velocity ( ⁠ V S "> V S ⁠ ) of the upper 30 m ( ⁠ V S 30 "> V S 30 ⁠ ) and the site dominant frequency ( ⁠ f d "> f d ⁠ ) to calculate d M L "> d M L ⁠ ⁠ . Previously established d M L "> d M L ⁠ ⁠ , measured V S 30 "> V S 30 ⁠ ⁠ , and f d "> f d data are available for between 126 and 458 horizontal components (east–west and north–south) at 81 seismic stations in the California Integrated Seismic Network; d M L "> d M L ⁠ d M L "> data range from −1.10 to 0.39, V S 30 "> V S 30 ⁠ values range from 202 to 1464    m / s "> 1464 m / s ⁠ , and 440 f d "> f d values are compiled from earthquake and microseismic records that range from 0.13 to 21 Hz. We find V S 30 "> V S 30 ⁠ and d M L "> d M L ⁠ exhibit a positive coefficient of determination ( ⁠ R = 0.59 "> R = 0.59 ⁠ ), indicating that as V S 30 "> V S 30 ⁠ increases, d M L "> d M L ⁠ increases. This implies that greater site amplification (lower V S 30 "> V S 30 ⁠ ⁠ ) results in smaller d M L "> d M L ⁠ ⁠ . f d "> f d and d M L "> d M L ⁠ also generally exhibit a positive correlation ( ⁠ R 2 < 0.56 "> R 2 < 0.56 ⁠ ), which implies lower d M L "> d M L values are related to site resonance at depth‐dependent frequencies. Using the developed relationships, V S 30 "> V S 30 ⁠ or f d "> f d measurements can be used to establish a provisional d M L "> d M L ⁠ for newly established stations. This procedure allows new stations to contribute to regional network M L "> M L determinations immediately without the need to wait until a minimum set of earthquake data has been recorded.

California↗

USGS2020

Real-time performance of the PLUM earthquake early warning method during the 2019 M6.4 and M7.1 Ridgecrest, California, Earthquakes

We evaluate the timeliness and accuracy of ground‐motion‐based earthquake early warning (EEW) during the July 2019 M 6.4 and 7.1 Ridgecrest earthquakes. In 2018, we began retrospective and internal real‐time testing of the propagation of local undamped motion (PLUM) method for earthquake warning in California, Oregon, and Washington, with the potential that PLUM might one day be included in the ShakeAlert EEW system. A real‐time version of PLUM was running on one of the ShakeAlert EEW system’s development servers at the time of the 2019 Ridgecrest sequence, allowing us to evaluate the timeliness and accuracy of PLUM’s warnings for the M 6.4 and 7.1 mainshocks in real time with the actual data availability and latencies of the operational ShakeAlert EEW system. The latter is especially important because high‐data latencies during the M 7.1 earthquake degraded ShakeAlert’s performance. PLUM proved to be largely immune to these latencies. In this article, we present a retrospective analysis of PLUM performance and explore three potential regional alerting strategies ranging from spatially large regions (counties), to moderate‐size regions (National Weather Service public forecast zones), to high‐spatial specificity (50 km regular geographic grid). PLUM generated initial shaking forecasts for the two mainshocks 5 and 6 s after their respective origin times, and faster than the ShakeAlert system’s first alerts. PLUM was also able to accurately forecast shaking across southern California for all three alerting strategies studied. As would be expected, a cost‐benefit analysis of each approach illustrates trade‐offs between increasing warning time and minimizing the area receiving unneeded alerts. Choosing an optimal alerting strategy requires knowledge of users’ false alarm tolerance and minimum required warning time for taking protective action, as well as the time required to distribute alerts to users.

California↗

USGS2020

Minimal clustering of injection-induced earthquakes observed with a large-n seismic array

The clustering behavior of injection‐induced earthquakes is examined using one month of data recorded by the LArge‐ n Seismic Survey in Oklahoma (LASSO) array. The 1829‐node seismic array was deployed in a 25 km×32 km area of active saltwater disposal in northern Oklahoma between 14 April and 10 May 2016. Injection rates in the study area are nearly constant around the time of the deployment. We develop a local magnitude ( ⁠ M L ⁠ ) equation for the region and estimate magnitudes for 1104 earthquakes recorded by the deployment. The determined earthquake magnitudes range from M L 0.01 to 3.0. The majority of earthquakes occurred between 1.5 and 5.5 km depth, and the shallowest earthquake depths overlap with the base of injection wells at depths between 1.5 and 2.5 km. We compute focal mechanisms of the largest events ( ⁠ M L > 2.0 ⁠ ), and find a mix of normal‐ and strike‐slip‐faulting types. Earthquakes occur regularly in time during the deployment, but are not evenly distributed in space across the study area, that is, they are spatially clustered. Analysis of the nearest‐neighbor distances in the space–time–magnitude domain shows the seismicity is dominated by single‐event clusters (i.e., independent events). This high proportion of single‐event clusters compared with multievent clusters has been previously noted for induced events at geothermal sites. When clustering occurs, the number of events in a cluster is typically small. We observe only four clusters with 10 or more events. For these larger clusters, we find equivalent numbers of foreshocks and aftershocks; however, the foreshock sequences are significantly longer in duration lasting days to tens of days, while aftershock sequences are observed only on the order of one day. The minimal clustering observed for events in the LASSO array suggests that the majority of events are being directly driven by stress changes due to local saltwater disposal.

Oklahoma↗

USGS2020

Activation of optimally and unfavourably oriented faults in a uniform local stress field during the 2011 Prague, Oklahoma, sequence

The orientations of faults activated relative to the local principal stress directions can provide insights into the role of pore pressure changes in induced earthquake sequences. Here, we examine the 2011 M 5.7 Prague earthquake sequence that was induced by nearby wastewater disposal. We estimate the local principal compressive stress direction near the rupture as inferred from shear wave splitting measurements at spatial resolutions as small as 750 m. We find that the dominant azimuth observed is parallel to previous estimates of the regional compressive stress with some secondary azimuths oriented subparallel to the strike of the major fault structures. From an extended catalogue, we map ten distinct fault segments activated during the sequence that exhibit a wide array of orientations. We assess whether the five near-vertical fault planes are optimally oriented to fail in the determined stress field. We find that only two of the fault planes, including the M 5.7 main shock fault, are optimally oriented. Both the M 4.8 foreshock and M 4.8 aftershock occur on fault planes that deviate 20–29° from the optimal orientation for slip. Our results confirm that induced event sequences can occur on faults not optimally oriented for failure in the local stress field. The results suggest elevated pore fluid pressures likely induced failure along several of the faults activated in the 2011 Prague sequence.

Oklahoma↗

USGS2020

Near-field ground motions from the July, 2019 Ridgecrest, California, earthquake sequence

The 2019 Ridgecrest, California, earthquake sequence, including an Mw 6.4 event on 4 July and an Mw 7.1 approximately 34 hr later, was recorded by 15 instruments within 55 km nearest‐fault distance. To characterize and explore near‐field ground motions from the Mw 6.4 foreshock and Mw 7.1 mainshock, we augment these records with available macroseismic information, including conventional intensities and displaced rocks. We conclude that near‐field shaking intensities were generally below modified Mercalli intensity 9, with concentrations of locally high values toward the northern and southern termini of the mainshock rupture. We further show that, relative to near‐field ground motions at hard‐rock sites, instrumental ground motions at alluvial near‐field sites for both the Mw 6.4 foreshock and Mw 7.1 mainshock were depleted in energy at frequencies higher than 2–3 Hz, as expected from ground‐motion models. Both the macroseismic and instrumental observations suggest that sediments in the Indian Wells Valley experienced a pervasively nonlinear response, which helps explain why shaking intensities and damage in the closest population center, Ridgecrest, were relatively modest given its proximity to the earthquakes.

California↗

USGS2020

How often can Earthquake Early Warning systems alert sites with high intensity ground motion?

Although numerous Earthquake Early Warning (EEW) algorithms have been developed we still lack a detailed understanding of how often and under what circumstances useful ground motion alerts can be provided to end-users. Here we analyze the alerting performance of the PLUM, EPIC and FinDer algorithms by running them retrospectively on the seismic strong motion data of the 219 earthquakes in Japan since 1996 that exceeded Modified Mercalli Intensity (MMI) of 4.5 on at least 10 sites (Mw 4.5-9.1). Our analysis suggests that, irrespective of the algorithm, EEW end-users should be prepared that EEW can often but not always provide useful ground motion alerts. A majority of sites with moderate-strong ground motion (MMI 5-6) can generally get at least a few seconds of warning time from all algorithms. If such shaking is caused by a shallow crustal event, around 50% of such sites receive alerts with warning times >5 s. Many sites with severe-extreme ground motion (MMI >=8) can be alerted successfully in the case of very large offshore earthquakes, but less than 20% can be alerted ahead of time if such shaking is caused by a shallow crustal event. Our results provide detailed quantitative insight into the expected alerting performance for EEW algorithms under realistic conditions. The main caveat is that the largest shallow crustal event in our data set has Mw7.0, i.e. the data set does not contain very large strike slip events.

Journal of Geophysical Research↗

USGS2020

The U.S. Geological Survey’s Rapid Seismic Array Deployment for the 2019 Ridgecrest Earthquake Sequence

Rapid seismic deployments following large earthquakes capture ephemeral near‐field recordings of aftershocks and ambient noise that can provide valuable data for seismological studies. The U.S. Geological Survey installed 19 temporary seismic stations following the 4 July 2019 M w 6.4 and 6 July 2019 (UTC) M w 7.1 earthquakes near the city of Ridgecrest, California. The stations record the aftershock sequence beginning two days after the mainshock and are expected to remain in the field through approximately January 2020. The deployment augments the permanent seismic network in the area to improve azimuthal coverage and provide additional near‐field observations. This article summarizes the motivation and goals of the deployment; details of station installation, instrumentation, and configurations; and initial data quality and observations from the network. We expect these data to be useful for a range of studies including detailing near‐field variability in strong ground motions, determining stress drops and rupture directivity of small events, imaging the fault zone, documenting the evolution of crustal properties within and outside of the fault zone, and others.

California↗

USGS2020

3D fault architecture controls the dynamism of earthquake swarm

The vibrant evolutionary patterns made by earthquake swarms are incompatible with standard, effectively two-dimensional (2D) models for general fault architecture. We leverage advances in earthquake monitoring with a deep-learning algorithm to image a fault zone hosting a 4-year-long swarm in southern California. We infer that fluids are naturally injected into the fault zone from below and diffuse through strike-parallel channels while triggering earthquakes. A permeability barrier initially limits up-dip swarm migration but ultimately is circumvented. This enables fluid migration within a shallower section of the fault with fundamentally different mechanical properties. Our observations provide high-resolution constraints on the processes by which swarms initiate, grow, and arrest. These findings illustrate how swarm evolution is strongly controlled by 3D variations in fault architecture.

Science↗