Search USGS⌕ Search

Geology topics

Nicholas J. van der Elst

Publications and source records attributed to Nicholas J. van der Elst.

31 records · Page 2Linked to original sources
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↗
USGS2019

#EarthquakeAdvisory: Exploring discourse between government officials, news media and social media during the Bombay Beach 2016 Swarm

Communicating probabilities of natural hazards to varied audiences is a notoriously difficult task. Many of these challenges were encountered during the 2016 Bombay Beach, California, swarm of ~100 2≤M≤4.3 earthquakes, which began on 26 September 2016 and lasted for several days. The swarm’s proximity to the southern end of the San Andreas fault caused concern that a larger earthquake could be triggered. Within 1–2 days, different forecast models were used to evaluate the likelihood of a larger event with two agencies (the U.S. Geological Survey [USGS] and the California Governor’s Office of Emergency Services) releasing probabilities and forecasts for larger earthquakes. Our research explores communication and news media efforts, as well as how people on a microblogging social media site (Twitter) responded to these forecasts. Our findings suggest that news media used a combination of information sources, basing their articles on what they learned from social media, as well as using information provided by government agencies. As the swarm slowed down, there is evidence of the continued interplay between news media and social media, with the USGS issuing revised probability reports and scientists from the USGS and other institutions participating in media interviews. In reporting on the swarm, news media often used language more generally than the scientists; terms such as probability, likelihood, chance, and possibility were used interchangeably. Knowledge of how news media used scientific information from the 2016 Bombay Beach forecasts can assist local, state, and federal agencies in developing effective communication strategies to respond to future earthquakes.

California↗
USGS2019

Improving earthquake forecasts during swarms with a duration model

Earthquake swarms present a challenge for operational earthquake forecasting because they are driven primarily by transient external processes, such as fluid flow, the behavior and duration of which are difficult to predict. In this study, we develop a swarm duration model to estimate how long a swarm is likely to last based on actuarial statistics of previous swarms in a given region. We demonstrate this approach using swarms that have been identified in the Salton trough in southern California, finding that swarms last an average of ∼ 7 "> ∼ 7 ∼7 days and have a relatively constant 15%–16% chance of terminating each day for the first 14 days of the swarm. Cataloged swarm durations are exponentially distributed, so we use a Poissonian model for swarm termination to encapsulate and extend the actuarial statistics. We then show how using the swarm duration model would have affected the earthquake forecast that was released during the 2016 Bombay Beach swarm. The earthquake forecast is substantially improved by incorporating a probabilistic model for how long the swarm is likely to last.

Bulletin of the Seismological Society of America↗
USGS2018

Updated California aftershock parameters

Reasenberg and Jones (1989) introduced a statistical model for aftershock rate following a mainshock along with estimates of “generic” California parameter values based on past aftershock sequences. The Reasenberg and Jones (1989) model has been used for decades to issue aftershock forecasts following M ≥ 5 "> M ≥ 5 mainshocks in California. Here, we update the “generic” parameters for California through a fit to the aftershock sequences of M ≥ 5 "> M ≥ 5 mainshocks occurring since 1980. We find aftershock productivity values that are lower on average than the generic productivity reported by Reasenberg and Jones (1989) , likely because low‐productivity sequences were omitted from their analysis and possibly because of a trade‐off between productivity and b "> b ‐value. We confirm the observation of Llenos and Michael (2017) that southern California sequences are more productive on average than northern California sequences. The Mendocino area is much less productive; the hydrothermal areas in Long Valley, Coso, and the Salton Sea, in contrast, are much more productive. We also quantify the variability of the Reasenberg and Jones (1989) productivity parameter a "> a a between sequences with a normal distribution. This distribution of a "> a ‐values can be used to compute aftershock forecasts that include epistemic uncertainty and can be used as the prior for Bayesian updating of the a "> a a ‐value as a sequence progresses.

California↗
USGS2018

Turing-style tests for UCERF3 synthetic catalogs

Epidemic-Type Aftershock Sequence (ETAS) catalogs generated from the 3rd Uniform California Earthquake Rupture Forecast (UCERF3) model are unique in that they are the first to combine a complex, fault-based long-term forecast with short-term earthquake clustering statistics. We present Turing-style tests to examine whether these synthetic catalogs can successfully imitate observed earthquake behavior in California. We find that UCERF3-ETAS is more spatially diffuse than the observed historic catalog in California and that it is lacking quiet periods that are present in the real catalog. While mean aftershock productivity of the observed catalog is matched closely by UCERF3-ETAS, the real catalog has more inter-sequence productivity variability and small mainshocks have more foreshocks. In sum, we find that UCERF3-ETAS differs from the observed catalog in ways that are foreseeable from its modeling simplifications. The tests we present here can be used on any model which produces suites of synthetic catalogs; as such, in addition to providing avenues for future improvements to the model, they could also be incorporated into testing platforms such as Collaboratory for the Study of Earthquake Predictability (CSEP).

Bulletin of the Seismological Society of America↗
USGS2017

Accounting for orphaned aftershocks in the earthquake background rate

Aftershocks often occur within cascades of triggered seismicity in which each generation of aftershocks triggers an additional generation, and so on. The rate of earthquakes in any particular generation follows Omori's law, going approximately as 1/ t . This function decays rapidly, but is heavy-tailed, and aftershock sequences may persist for long times at a rate that is difficult to discriminate from background. It is likely that some apparently spontaneous earthquakes in the observational catalogue are orphaned aftershocks of long-past main shocks. To assess the relative proportion of orphaned aftershocks in the apparent background rate, I develop an extension of the ETAS model that explicitly includes the expected contribution of orphaned aftershocks to the apparent background rate. Applying this model to California, I find that the apparent background rate can be almost entirely attributed to orphaned aftershocks, depending on the assumed duration of an aftershock sequence. This implies an earthquake cascade with a branching ratio (the average number of directly triggered aftershocks per main shock) of nearly unity. In physical terms, this implies that very few earthquakes are completely isolated from the perturbing effects of other earthquakes within the fault system. Accounting for orphaned aftershocks in the ETAS model gives more accurate estimates of the true background rate, and more realistic expectations for long-term seismicity patterns.

Geophysical Journal International↗
USGS2017

Tidal triggering of earthquakes suggests poroelastic behavior on the San Andreas Fault

Tidal triggering of earthquakes is hypothesized to provide quantitative information regarding the fault's stress state, poroelastic properties, and may be significant for our understanding of seismic hazard. To date, studies of regional or global earthquake catalogs have had only modest successes in identifying tidal triggering. We posit that the smallest events that may provide additional evidence of triggering go unidentified and thus we developed a technique to improve the identification of very small magnitude events. We identify events applying a method known as inter-station seismic coherence where we prioritize detection and discrimination over characterization. Here we show tidal triggering of earthquakes on the San Andreas Fault. We find the complex interaction of semi-diurnal and fortnightly tidal periods exposes both stress threshold and critical state behavior. Our findings reveal earthquake nucleation processes and pore pressure conditions – properties of faults that are difficult to measure, yet extremely important for characterizing earthquake physics and seismic hazards.

California↗
USGS2016

Three ingredients for Improved global aftershock forecasts: Tectonic region, time-dependent catalog incompleteness, and inter-sequence variability

Following a large earthquake, seismic hazard can be orders of magnitude higher than the long‐term average as a result of aftershock triggering. Because of this heightened hazard, emergency managers and the public demand rapid, authoritative, and reliable aftershock forecasts. In the past, U.S. Geological Survey (USGS) aftershock forecasts following large global earthquakes have been released on an ad hoc basis with inconsistent methods, and in some cases aftershock parameters adapted from California. To remedy this, the USGS is currently developing an automated aftershock product based on the Reasenberg and Jones (1989) method that will generate more accurate forecasts. To better capture spatial variations in aftershock productivity and decay, we estimate regional aftershock parameters for sequences within the García et al. (2012) tectonic regions. We find that regional variations for mean aftershock productivity reach almost a factor of 10. We also develop a method to account for the time‐dependent magnitude of completeness following large events in the catalog. In addition to estimating average sequence parameters within regions, we develop an inverse method to estimate the intersequence parameter variability. This allows for a more complete quantification of the forecast uncertainties and Bayesian updating of the forecast as sequence‐specific information becomes available.

California↗
USGS2016

Fortnightly modulation of San Andreas tremor and low-frequency earthquakes

Earth tides modulate tremor and low-frequency earthquakes (LFEs) on faults in the vicinity of the brittle−ductile (seismic−aseismic) transition. The response to the tidal stress carries otherwise inaccessible information about fault strength and rheology. Here, we analyze the LFE response to the fortnightly tide, which modulates the amplitude of the daily tidal stress over a 14-d cycle. LFE rate is highest during the waxing fortnightly tide, with LFEs most strongly promoted when the daily stress exceeds the previous peak stress by the widest margin. This pattern implies a threshold failure process, with slip initiated when stress exceeds the local fault strength. Variations in sensitivity to the fortnightly modulation may reflect the degree of stress concentration on LFE-producing brittle asperities embedded within an otherwise aseismic fault.

California↗
USGS2016

Induced earthquake magnitudes are as large as (statistically) expected

A major question for the hazard posed by injection-induced seismicity is how large induced earthquakes can be. Are their maximum magnitudes determined by injection parameters or by tectonics? Deterministic limits on induced earthquake magnitudes have been proposed based on the size of the reservoir or the volume of fluid injected. However, if induced earthquakes occur on tectonic faults oriented favorably with respect to the tectonic stress field, then they may be limited only by the regional tectonics and connectivity of the fault network. In this study, we show that the largest magnitudes observed at fluid injection sites are consistent with the sampling statistics of the Gutenberg-Richter distribution for tectonic earthquakes, assuming no upper magnitude bound. The data pass three specific tests: (1) the largest observed earthquake at each site scales with the log of the total number of induced earthquakes, (2) the order of occurrence of the largest event is random within the induced sequence, and (3) the injected volume controls the total number of earthquakes rather than the total seismic moment. All three tests point to an injection control on earthquake nucleation but a tectonic control on earthquake magnitude. Given that the largest observed earthquakes are exactly as large as expected from the sampling statistics, we should not conclude that these are the largest earthquakes possible. Instead, the results imply that induced earthquake magnitudes should be treated with the same maximum magnitude bound that is currently used to treat seismic hazard from tectonic earthquakes.

Journal of Geophysical Research B: Solid Earth↗
USGS2015

Larger aftershocks happen farther away: nonseparability of magnitude and spatial distributions of aftershocks

Aftershocks may be driven by stress concentrations left by the main shock rupture or by elastic stress transfer to adjacent fault sections or strands. Aftershocks that occur within the initial rupture may be limited in size, because the scale of the stress concentrations should be smaller than the primary rupture itself. On the other hand, aftershocks that occur on adjacent fault segments outside the primary rupture may have no such size limitation. Here we use high-precision double-difference relocated earthquake catalogs to demonstrate that larger aftershocks occur farther away than smaller aftershocks, when measured from the centroid of early aftershock activity—a proxy for the initial rupture. Aftershocks as large as or larger than the initiating event nucleate almost exclusively in the outer regions of the aftershock zone. This observation is interpreted as a signature of elastic rebound in the earthquake catalog and can be used to improve forecasting of large aftershocks.

Geophysical Research Letters↗
USGS2008

Complex faulting associated with the 22 December 2003 M w 6.5 San Simeon California, earthquake, aftershocks and postseismic surface deformation

We use data from two seismic networks and satellite interferometric synthetic aperture radar (InSAR) imagery to characterize the 22 December 2003 M w 6.5 San Simeon earthquake sequence. Absolute locations for the mainshock and nearly 10,000 aftershocks were determined using a new three-dimensional (3D) seismic velocity model; relative locations were obtained using double difference. The mainshock location found using the 3D velocity model is 35.704° N, 121.096° W at a depth of 9.7±0.7 km. The aftershocks concentrate at the northwest and southeast parts of the aftershock zone, between the mapped traces of the Oceanic and Nacimiento fault zones. The northwest end of the mainshock rupture, as defined by the aftershocks, projects from the mainshock hypocenter to the surface a few kilometers west of the mapped trace of the Oceanic fault, near the Santa Lucia Range front and the >5 mm postseismic InSAR imagery contour. The Oceanic fault in this area, as mapped by Hall (1991) , is therefore probably a second-order synthetic thrust or reverse fault that splays upward from the main seismogenic fault at depth. The southeast end of the rupture projects closer to the mapped Oceanic fault trace, suggesting much of the slip was along this fault, or at a minimum is accommodating much of the postseismic deformation. InSAR imagery shows ∼72 mm of postseismic uplift in the vicinity of maximum coseismic slip in the central section of the rupture, and ∼48 and ∼45 mm at the northwest and southeast end of the aftershock zone, respectively. From these observations, we model a ∼30-km-long northwest-trending northeast-dipping mainshock rupture surface—called the mainthrust—which is likely the Oceanic fault at depth, a ∼10-km-long southwest-dipping backthrust parallel to the mainthrust near the hypocenter, several smaller southwest-dipping structures in the southeast, and perhaps additional northeast-dipping or subvertical structures southeast of the mainshock plane. Discontinuous backthrust features opposite the mainthrust in the southeast part of the aftershock zone may offset the relic Nacimiento fault zone at depth. The InSAR data image surface deformation associated with both aseismic slip and aftershock production on the mainthrust and the backthrusts at the northwest and southeast ends of the aftershock zone. The well-defined mainthrust at the latitude of the epicenter and antithetic backthrust illuminated by the aftershock zone indicate uplift of the Santa Lucia Range as a popup block; aftershocks in the southeast part of the zone also indicate a popup block, but it is less well defined. The absence of backthrust features in the central part of the zone suggests range-front uplift by fault-propagation folding, or backthrusts in the central part were not activated during the mainshock.

California↗