Geology topics

Faults & earthquakes

Explore USGS earthquake and fault publications. Search fault names, locations and seismic hazards; follow the original report for maps, methods and publication context.

Follow this topic with RSS

At least 937 records · Page 52Linked to original sources

USGS2026

Wavelet Inversion for SliP (WISP): Open-source earthquake slip modeling software

Models of the spatiotemporal evolution of earthquake slip, termed finite-fault models, are a critical component of rapid earthquake and tsunami response, earthquake forecasting, seismic ground-motion estimates, and studies of earthquake kinematics. Here, we detail a newly released finite-fault modeling software, Wavelet Inversion for SliP (WISP), in use at the U.S. Geological Survey’s National Earthquake Information Center (NEIC) and available to the public. WISP version 1.1.0 allows inversion of teleseismic body and surface waves, as well as local strong-motion, static and dynamic Global Navigation Satellite System, and satellite imagery (e.g., Interferometric Synthetic Aperture Radar) observations on single or multiple planar fault segments. The software is used in NEIC rapid response of earthquakes M w ≥ 7, generally resulting in a published model within the first few hours after the event origin time. The rupture location and dimensions are then used as inputs to downstream products to estimate earthquake shaking, predict loss, and model the likelihood of secondary hazards, namely landslides and liquefaction. WISP is also used in research studies to evaluate the characteristics of complex ruptures including multifault ruptures and earthquake doublets, among others. The WISP version 1.1.0 software release is composed of Python-wrapped FORTRAN code to accomplish the inversion procedure. A simple command line interface facilitates ease of use even for those with only a cursory knowledge of Python scripting. WISP version 1.1.0 includes a Jupyter Notebook tutorial demonstrating use of the software for modeling the 2015 M w 8.3 Illapel, Chile, earthquake. In parallel with the tutorial, we demonstrate the typical usage of the WISP software using the M w 8.3 Illapel earthquake example here.

Seismological Research Letters↗

USGS2010

The MW 7.0 Haiti Earthquake of January 12, 2010: USGS/EERI Advance Reconnaissance Team Report

Executive Summary A field reconnaissance in Haiti by a five-member team with expertise in seismology and earthquake engineering has revealed a number of factors that led to catastrophic losses of life and property during the January 12, 2010, Mw 7.0 earthquake. The field study was conducted from January 26 to February 3, 2010, and included investigations in Port-au-Prince and the heavily damaged communities to the west, including Leogane, Grand Goave, Petite Goave, and Oliver. Seismology Despite recent seismic quiescence, Haiti has suffered similar devastating earthquakes in the historical past (1701, 1751, 1770 and 1860). Despite this knowledge of historical seismicity, Haiti had no seismograph stations during the main earthquake, so it is impossible to estimate accurately the intensity of ground motions. Nonetheless, the wide range of buildings damaged by the January 12, 2010, earthquake suggests that the ground motions contained seismic energy over a wide range of frequencies. Another earthquake of similar magnitude could strike at any time on the eastern end of the Enriquillo Fault, directly to the south of Port-au-Prince. Reconstruction must take this hazard into account. The four portable seismographs installed by the team recorded a series of small aftershocks. As expected, the ground motions recorded at a hard-rock site contained a greater proportion of high frequencies than the motions recorded at a soil site. Two of the stations continue to monitor seismic activity. A thorough field investigation of the mapped Enriquillo Fault south of the city of Leogane failed to find any evidence of surface faulting. This led the team to conclude that the earthquake was unlikely to have produced any surface rupture in the study area. Geotechnical Aspects Soil liquefaction, landslides and rockslides in cut slopes, and road embankment failures contributed to extensive damage in Port-au-Prince and elsewhere. A lack of detailed knowledge of the physical conditions of the soils (for example, lithology, stiffness, density, and thickness) made it difficult for us to quantitatively assess the role of ground-motion amplification in the widespread damage. Buildings The Haitian Ministry of Statistics and Informatics reported that one-story buildings represent 73 percent of the building inventory. Most ordinary, one-story houses have roofs made of sheet metal (82 percent), whereas most multistory houses and apartments have roofs made of concrete (71 percent). Walls made of concrete/block/stone predominate both in ordinary houses and apartments. It appears that the widespread damage to residences and commercial and government buildings was attributable to a great extent to the lack of earthquake-resistant design. In many cases, the structural types, member dimensions, and detailing practices were inadequate to resist strong ground motions. These vulnerabilities may have been exacerbated by poor construction practices. Reinforced concrete frames with concrete block masonry infill appeared to perform particularly poorly. Structures with light (timber or sheet metal) roofs performed better compared to structures with concrete roofs and slabs. The seismic performance of some buildings was adequate, and some of the damaged buildings appeared to have had low deformation demands. These observations suggest that structures designed and constructed with adequate stiffness and reinforcing details would have resisted the earthquake without being damaged severely. A damage survey of 107 buildings in downtown Port-au-Prince indicated that 28 percent had collapsed and another 33 percent were damaged enough to require repairs. A similar survey of 52 buildings in Leogane found that 62 percent had collapsed and another 31 percent required repairs. Bridges There was no evidence of bridge collapses attributable to the earthquake. Most bridges in Port-au-Prince are simple box culverts consisting of box girders 2.0 to 2.

Open-File Report↗

USGS2003

Imaging the complexity of an active normal fault system: The 1997 Colfiorito (central Italy) case study

Six moderate magnitude earthquakes (5 < Mw < 6) ruptured normal fault segments of the southern sector of the North Apennine belt (central Italy) in the 1997 Colfiorito earthquake sequence. We study the progressive activation of adjacent and nearby parallel faults of this complex normal fault system using ???1650 earthquake locations obtained by applying a double-difference location method, using travel time picks and waveform cross-correlation measurements. The lateral extent of the fault segments range from 5 to 10 km and make up a broad, ???45 km long, NW trending fault system. The geometry of each segment is quite simple and consists of planar faults gently dipping toward SW with an average dip of 40??-45??. The fault planes are not listric but maintain a constant dip through the entire seismogenic volume, down to 8 km depth. We observe the activation of faults on the hanging wall and the absence of seismicity in the footwall of the structure. The observed fault segmentation appears to be due to the lateral heterogeneity of the upper crust: preexisting thrusts inherited from Neogene's compressional tectonic intersect the active normal faults and control their maximum length. The stress tensor obtained by inverting the six main shock focal mechanisms of the sequence is in agreement with the tectonic stress active in the inner chain of the Apennine, revealing a clear NE trending extension direction. Aftershock focal mechanisms show a consistent extensional kinematics, 70% of which are mechanically consistent with the main shock stress field.

Journal of Geophysical Research B: Solid Earth↗

USGS2015

Structure of the Koyna-Warna Seismic Zone, Maharashtra, India: A possible model for large induced earthquakes elsewhere

The Koyna-Warna area of India is one of the best worldwide examples of reservoir-induced seismicity, with the distinction of having generated the largest known induced earthquake (M6.3 on 10 December 1967) and persistent moderate-magnitude (>M5) events for nearly 50 years. Yet, the fault structure and tectonic setting that has accommodated the induced seismicity is poorly known, in part because the seismic events occur beneath a thick sequence of basalt layers. On the basis of the alignment of earthquake epicenters over an ~50 year period, lateral variations in focal mechanisms, upper-crustal tomographic velocity images, geophysical data (aeromagnetic, gravity, and magnetotelluric), geomorphic data, and correlation with similar structures elsewhere, we suggest that the Koyna-Warna area lies within a right step between northwest trending, right-lateral faults. The sub-basalt basement may form a local structural depression (pull-apart basin) caused by extension within the step-over zone between the right-lateral faults. Our postulated model accounts for the observed pattern of normal faulting in a region that is dominated by north-south directed compression. The right-lateral faults extend well beyond the immediate Koyna-Warna area, possibly suggesting a more extensive zone of seismic hazards for the central India area. Induced seismic events have been observed many places worldwide, but relatively large-magnitude induced events are less common because critically stressed, preexisting structures are a necessary component. We suggest that releasing bends and fault step-overs like those we postulate for the Koyna-Warna area may serve as an ideal tectonic environment for generating moderate- to large- magnitude induced (reservoir, injection, etc.) earthquakes.

Journal of Geophysical Research↗

USGS1975

The nature of surface tilt along 85 km of the San Andreas fault-preliminary results form a 14-instrument array

The continuous monitoring of surface deformation near active faults is clearly necessary for an understanding of elastic strain accumulation and elastic and anelastic strain release associated with earthquakes. Fourteen 2-component tiltmeters have been installed in shallow boreholes along 85 km of the currently most active section of the San Andreas fault in the western United States. These instruments operate at a sensitivity of 10 −8 radians. Five of these tiltmeters, extending along one 35 km section of the fault, have been in operation since June 1973. The results indicate that regional tectonic tilting has occurred before more than ten individual earthquakes or groups of earthquakes with epicenters within ten earthquake source dimensions of one or more instruments. This tilting has a time scale of up to a month depending on earthquake magnitude. The amplitude of these tilts exceeds by almost an order of magnitude that expected from a dislocation model of the source using seismically determined parameters. No indication of rapid or accelerated tilt just prior to these earthquakes has been seen.

California↗

USGS1972

Prospects for earthquake prediction and control

The San Andreas fault is viewed, according to the concepts of seafloor spreading and plate tectonics, as a transform fault that separates the Pacific and North American plates and along which relative movements of 2 to 6 cm/year have been taking place. The resulting strain can be released by creep, by earthquakes of moderate size, or (as near San Francisco and Los Angeles) by great earthquakes. Microearthquakes, as mapped by a dense seismograph network in central California, generally coincide with zones of the San Andreas fault system that are creeping. Microearthquakes are few and scattered in zones where elastic energy is being stored. Changes in the rate of strain, as recorded by tiltmeter arrays, have been observed before several earthquakes of about magnitude 4. Changes in fluid pressure may control timing of seismic activity and make it possible to control natural earthquakes by controlling variations in fluid pressure in fault zones. An experiment in earthquake control is underway at the Rangely oil field in Colorado, where the rates of fluid injection and withdrawal in experimental wells are being controlled.

California↗

USGS2025

Probing lower-crustal fault properties with frequency-dependent tidal tremor triggering

The way seismicity responds to periodic stress perturbations offers crucial insights into the processes that can trigger an earthquake. Laboratory and theoretical analyses have shown that the period of imposed forcing and source properties affect the sensitivity to triggering, but frequency-dependent triggering of tectonic faults is poorly understood. The rate of low-frequency earthquakes (LFEs) near Parkfield, California has been found to be strongly correlated with solid earth tides. Tidal forcing acts over multiple frequencies, and the sensitivity to tidal triggering of LFEs therefore provides a unique opportunity to probe the physics of earthquake triggering and underlying fault properties. Here, we find that the response of LFEs to solid earth tides at diurnal and semi-diurnal frequencies is highly variable but spatially coherent along the San Andreas Fault. Using rate-state friction modeling, we find that the variation of the amplitude of tidal modulation is mainly affected by the spatial variation of the background effective normal stress, whereas the spatially varied tidal modulation at diurnal and semi-diurnal frequencies is mainly affected by the fault frictional property. The spatial variations of the LFE response reveal a heterogeneous lower-crustal geologic structure and complex physical faulting processes below the rupture zone of eventual great earthquakes.

California↗

USGS2004

Glacier ice mass fluctuations and fault instability in tectonically active Southern Alaska

Across the plate boundary zone in south central Alaska, tectonic strain rates are high in a region that includes large glaciers undergoing wastage (glacier retreat and thinning) and surges. For the coastal region between the Bering and Malaspina Glaciers, the average ice mass thickness changes between 1995 and 2000 range from 1 to 5 m/year. These ice changes caused solid Earth displacements in our study region with predicted values of -10 to 50 mm in the vertical and predicted horizontal displacements of 0-10 mm at variable orientations. Relative to stable North America, observed horizontal rates of tectonic deformation range from 10 to 40 mm/year to the north-northwest and the predicted tectonic uplift rates range from approximately 0 mm/year near the Gulf of Alaska coast to 12 mm/year further inland. The ice mass changes between 1995 and 2000 resulted in discernible changes in the Global Positioning System (GPS) measured station positions of one site (ISLE) located adjacent to the Bagley Ice Valley and at one site, DON, located south of the Bering Glacier terminus. In addition to modifying the surface displacements rates, we evaluated the influence ice changes during the Bering glacier surge cycle had on the background seismic rate. We found an increase in the number of earthquakes (ML???2.5) and seismic rate associated with ice thinning and a decrease in the number of earthquakes and seismic rate associated with ice thickening. These results support the hypothesis that ice mass changes can modulate the background seismic rate. During the last century, wastage of the coastal glaciers in the Icy Bay and Malaspina region indicates thinning of hundreds of meters and in areas of major retreat, maximum losses of ice thickness approaching 1 km. Between the 1899 Yakataga and Yakutat earthquakes (Mw=8.1, 8.1) and prior to the 1979 St. Elias earthquake (M s=7.2), the plate interface below Icy Bay was locked and tectonic strain accumulated. We used estimated ice mass change during the 1899-1979 time period to calculate the change in the fault stability margin (FSM) prior to the 1979 St. Elias earthquake. Our results suggest that a cumulative decrease in the fault stability margin at seismogenic depths, due to ice wastage over 80 years, was large, up to ???2 MPa. Ice wastage would promote thrust faulting in events such as the 1979 earthquake and subsequent aftershocks.

Conference Paper↗

USGS1998

Review of magnetic field monitoring near active faults and volcanic calderas in California: 1974-1995

Differential magnetic fields have been monitored along the San Andreas fault and the Long Valley caldera since 1974. At each monitoring location, proton precession magnetometers sample total magnetic field intensity at a resolution of 0.1 nT or 0.25 nT. Every 10 min, data samples are transmitted via satellite telemetry to Menlo Park, CA for processing and analysis. The number of active magnetometer sites has varied during the past 21 years from 6 to 25, with 12 sites currently operational. We use this network to identify magnetic field changes generated by earthquake and volcanic processes. During the two decades of monitoring, five moderate earthquakes (M5.9 to M7.3) have occurred within 20 km of magnetometer sites located along the San Andreas fault and only one preseismic signal of 1.5 nT has been observed. During moderate earthquakes, coseismic magnetic signals, with amplitudes from 0.7 nT to 1.3 nT, have been identified for 3 of the 5 events. These observations are generally consistent with those calculated from simple seismomagnetic models of these earthquakes and near-fault coseismic magnetic field disturbances rarely exceed one nanotesla. These data are consistent with the concept of low shear stress and relatively uniform displacement of the San Andreas fault system as expected due to high pore fluid pressure on the fault. A systematic decrease of 0.8-1 nT/year in magnetic field has occurred in the Long Valley caldera since 1989. These magnetic field data are similar in form to observed geodetically measured displacements from inflation of the resurgent dome. A simple volcanomagnetic model involving pressure increase of 50 MPa/a at a depth of 7 km under the resurgent dome can replicate these magnetic field observations. This model is derived from the intrusion model that best fits the surface deformation data. ?? 1998 Elsevier Science B.V.

Physics of the Earth and Planetary Interiors↗

USGS2009

Constraining earthquake source inversions with GPS data: 1. Resolution-based removal of artifacts

We present a resolution analysis of an inversion of GPS data from the 2004 M w 6.0 Parkfield earthquake. This earthquake was recorded at thirteen 1-Hz GPS receivers, which provides for a truly coseismic data set that can be used to infer the static slip field. We find that the resolution of our inverted slip model is poor at depth and near the edges of the modeled fault plane that are far from GPS receivers. The spatial heterogeneity of the model resolution in the static field inversion leads to artifacts in poorly resolved areas of the fault plane. These artifacts look qualitatively similar to asperities commonly seen in the final slip models of earthquake source inversions, but in this inversion they are caused by a surplus of free parameters. The location of the artifacts depends on the station geometry and the assumed velocity structure. We demonstrate that a nonuniform gridding of model parameters on the fault can remove these artifacts from the inversion. We generate a nonuniform grid with a grid spacing that matches the local resolution length on the fault and show that it outperforms uniform grids, which either generate spurious structure in poorly resolved regions or lose recoverable information in well-resolved areas of the fault. In a synthetic test, the nonuniform grid correctly averages slip in poorly resolved areas of the fault while recovering small-scale structure near the surface. Finally, we present an inversion of the Parkfield GPS data set on the nonuniform grid and analyze the errors in the final model. Copyright 2009 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

USGS2019

Rupture model of the M5.8 Pawnee, Oklahoma earthquake from regional and teleseismic waveforms

The 2016 M5.8 Pawnee, Oklahoma earthquake is the largest earthquake to have been induced by wastewater disposal. We infer the coseismic slip history from analysis of apparent source time functions and inversion of regional and teleseismic P‐waveforms, using aftershocks as empirical Green's functions. The earthquake nucleated on the shallow part of the fault, initially rupturing towards the surface, followed shortly thereafter by slip deeper on the fault. Deeper slip occurred below the aftershocks and at greater depths than most induced seismicity in the region, suggesting that small‐ to moderate‐sized earthquakes may not occur on deeper parts of faults in Oklahoma because they are further from failure than shallower fault sections. Comparisons with models of pore pressure perturbations further suggest that the earthquake may have initiated within a region of higher pore pressure perturbation but was not confined to this zone. These observations inform source physics and understanding of maximum magnitudes.

Oklahoma↗

USGS2025

Timescales of surface faulting preservation in low-strain intraplate regions from landscape evolution modeling and the geomorphic and historical record

Large surface-rupturing intraplate earthquakes in stable continental regions (SCRs) are uncommon globally and have recurrence intervals of thousands to hundreds of thousands of years based on the paleoseismic and geomorphic record, challenging accurate active fault identification in these regions. To constrain the timescales of preservation for scarps created by surface ruptures from dip-slip earthquakes, we use a two-dimensional scarp diffusion model for typical intraplate settings and explore which parameters influence fault scarp preservation. These parameters include the coseismic vertical surface offset, the recurrence interval of similar magnitude earthquakes, diffusivity (as a proxy for mean annual precipitation rate), and the erodibility of the surficial material. We constrain parameter ranges from a compilation of historical surface ruptures in intraplate settings in a variety of climates, including the Central and Eastern United States, Australia, Europe, Central Asia (Mongolia, China), India, and West Africa. The timescales of scarp preservation from landscape evolution modeling agree well with observations of scarp preservation in low-strain SCR and intraplate tectonic settings, with some notable exceptions for Australian scarps. We find that the erodibility of the surficial material and earthquake recurrence interval have a stronger effect on the timescales of scarp preservation than diffusivity or coseismic vertical surface offset. Our model results may aid in identifying and characterizing subtle, slow-moving active faults in low-strain SCR and intraplate tectonic settings for different tectonic, geomorphic, and climatic characteristics. Accurate fault locations and characterization from the landscape record has implications for both probabilistic seismic and fault displacement hazard analyses.

Journal of Geophysical Research Solid Earth↗

USGS2025

Induced earthquakes are generally not tidally triggered in Oklahoma and Kansas

Human-induced earthquakes occur along critically stressed faults as injected wastewater simultaneously heightens fluid pressure and pushes faults to failure. We investigate the possibility that small stresses imposed by Earth tides could trigger earthquakes in the induced seismicity region of Oklahoma and Kansas from 2011 to 2018. We decluster a catalog consisting of ∼110,000 earthquakes using three methods (Reasenberg, nearest-neighbor distance, and phase-bin). We find no significant tidal earthquake triggering using Schuster's p -value test for the declustered catalogs as a whole. We search for localized triggering using discretized space-time cells and find ∼0–6% of cells have significant tidal triggering which is close to what is randomly expected (5%) and indicates there is an insignificant amount of tidal triggering for the full study region. One area that has significant p -values across multiple time windows, ∼2014–2016 is ∼15 km from a region of large wastewater injection volume. It is possible that localized tidal triggering occurs for this time and area because faults remain critically stressed and are particularly susceptible to slip under the small stress load from the semidiurnal tide. Possible explanations for the lack of tidal triggering in our broader study are that the pre-seismic stressing rate in the earthquake nucleation area is faster than the tidal stressing rate (∼3 kPa/day), faults are not close enough to critically stressed to be affected by tidal forcing, and that nucleation occurs over longer periods than the tides considered in this study (∼1, ∼14 days). Fluid injection could be the source of a higher pre-seismic stress rate.

Kansas, Oklahoma↗

USGS1997

Stressing of the New Madrid seismic zone by a lower crust detachment fault

A new mechanical model for the cause of the New Madrid seismic zone in the central United States is analyzed. The model contains a subhorizontal detachment fault which is assumed to be near the domed top surface of locally thickened anomalous lower crust (“rift pillow”). Regional horizontal compression induces slip on the fault, and the slip creates a stress concentration in the upper crust above the rift pillow dome. In the coseismic stage of the model earthquake cycle, where the three largest magnitude 7–8 earthquakes in 1811–1812 are represented by a single model mainshock on a vertical northeast trending fault, the model mainshock has a moment equivalent to a magnitude 8 event. During the interseismic stage, corresponding to the present time, slip on the detachment fault exerts a right-lateral shear stress on the locked vertical fault whose failure produces the model mainshock. The sense of shear is generally consistent with the overall sense of slip of 1811–1812 and later earthquakes. Predicted rates of horizontal strain at the ground surface are about 10 −7 year −1 and are comparable to some observed rates. The model implies that rift pillow geometry is a significant influence on the maximum possible earthquake magnitude.

Journal of Geophysical Research B: Solid Earth↗

USGS2012

Earthquakes in the Orozco transform zone: seismicity, source mechanisms, and tectonics

As part of the Rivera Ocean Seismic Experiment, a network of ocean bottom seismometers and hydrophones was deployed in order to determine the seismic characteristics of the Orozco transform fault in the central eastern Pacific. We present hypocentral locations and source mechanisms for 70 earthquakes recorded by this network. All epicenters are within the transform region of the Orozco Fracture Zone and clearly delineate the active plate boundary. About half of the epicenters define a narrow line of activity parallel to the spreading direction and situated along a deep topographic trough that forms the northern boundary of the transform zone (region 1). Most focal depths for these events are very shallow, within 4 km of the seafloor; several well-determined focal depths, however, are as great as 7 km. No shallowing of seismic activity is observed as the rise-transform intersection is approached; to the contrary, the deepest events are within 10 km of the intersection. First motion polarities for most of the earthquakes in region 1 are compatible with right-lateral strike slip faulting along a nearly vertical plane, striking parallel to the spreading direction. Another zone of activity is observed in the central part of the transform (region 2). The apparent horizontal and vertical distribution of activity in this region is more scattered than in the first, and the first motion radiation patterns of these events do not appear to be compatible with any known fault mechanism. Pronounced lateral variations in crustal velocity structure are indicated for the transform region from refraction data and measurements of wave propagation directions. The effect of this lateral heterogeneity on hypocenters and fault plane solutions is evaluated by tracing rays through a three-dimensional velocity grid. While findings for events in region 1 are not significantly affected, in region 2, epicentral mislocations of up to 10 km and azimuthal deflections of up to 45° may result from assuming a laterally homogeneous velocity structure. When corrected for the effects of lateral heterogeneity, the epicenters and fault plane solutions for earthquakes in region 2 are compatible with predominantly normal faulting along a topographic trough trending NW–SE; the focal depths, however, are poorly constrained. These results suggest an en echelon spreading center or leaky transform regime in the central transform region.

Journal of Geophysical Research B: Solid Earth↗

USGS1999

Frequency-magnitude statistics and spatial correlation dimensions of earthquakes at Long Valley caldera, California

Intense earthquake swarms at Long Valley caldera in late 1997 and early 1998 occurred on two contrasting structures. The first is defined by the intersection of a north-northwesterly array of faults with the southern margin of the resurgent dome, and is a zone of hydrothermal upwelling. Seismic activity there was characterized by high b-values and relatively low values of D, the spatial fractal dimension of hypocentres. The second structure is the pre-existing South Moat fault, which has generated large-magnitude seismic activity in the past. Seismicity on this structure was characterized by low b-values and relatively high D. These observations are consistent with low-magnitude, clustered earthquakes on the first structure, and higher-magnitude, diffuse earthquakes on the second structure. The first structure is probably an immature fault zone, fractured on a small scale and lacking a well-developed fault plane. The second zone represents a mature fault with an extensive, coherent fault plane.

Geophysical Journal International↗

USGS2018

Maturity of nearby faults influences seismic hazard from hydraulic fracturing

Understanding the causes of human-induced earthquakes is paramount to reducing societal risk. We investigated five cases of seismicity associated with hydraulic fracturing (HF) in Ohio since 2013 that, because of their isolation from other injection activities, provide an ideal setting for studying the relations between high-pressure injection and earthquakes. Our analysis revealed two distinct groups: ( i ) deeper earthquakes in the Precambrian basement, with larger magnitudes (M > 2), b-values < 1, and many post–shut-in earthquakes, versus ( ii ) shallower earthquakes in Paleozoic rocks ∼400 m below HF, with smaller magnitudes (M < 1), b-values > 1.5, and few post–shut-in earthquakes. Based on geologic history, laboratory experiments, and fault modeling, we interpret the deep seismicity as slip on more mature faults in older crystalline rocks and the shallow seismicity as slip on immature faults in younger sedimentary rocks. This suggests that HF inducing deeper seismicity may pose higher seismic hazards. Wells inducing deeper seismicity produced more water than wells with shallow seismicity, indicating more extensive hydrologic connections outside the target formation, consistent with pore pressure diffusion influencing seismicity. However, for both groups, the 2 to 3 h between onset of HF and seismicity is too short for typical fluid pressure diffusion rates across distances of ∼1 km and argues for poroelastic stress transfer also having a primary influence on seismicity.

Ohio↗

USGS2014

CyberShake-derived ground-motion prediction models for the Los Angeles region with application to earthquake early warning

Real-time applications such as earthquake early warning (EEW) typically use empirical ground-motion prediction equations (GMPEs) along with event magnitude and source-to-site distances to estimate expected shaking levels. In this simplified approach, effects due to finite-fault geometry, directivity and site and basin response are often generalized, which may lead to a significant under- or overestimation of shaking from large earthquakes ( M > 6.5) in some locations. For enhanced site-specific ground-motion predictions considering 3-D wave-propagation effects, we develop support vector regression (SVR) models from the SCEC CyberShake low-frequency (<0.5 Hz) and broad-band (0–10 Hz) data sets. CyberShake encompasses 3-D wave-propagation simulations of >415 000 finite-fault rupture scenarios (6.5 ≤ M ≤ 8.5) for southern California defined in UCERF 2.0. We use CyberShake to demonstrate the application of synthetic waveform data to EEW as a ‘proof of concept’, being aware that these simulations are not yet fully validated and might not appropriately sample the range of rupture uncertainty. Our regression models predict the maximum and the temporal evolution of instrumental intensity (MMI) at 71 selected test sites using only the hypocentre, magnitude and rupture ratio, which characterizes uni- and bilateral rupture propagation. Our regression approach is completely data-driven (where here the CyberShake simulations are considered data) and does not enforce pre-defined functional forms or dependencies among input parameters. The models were established from a subset (∼20 per cent) of CyberShake simulations, but can explain MMI values of all >400 k rupture scenarios with a standard deviation of about 0.4 intensity units. We apply our models to determine threshold magnitudes (and warning times) for various active faults in southern California that earthquakes need to exceed to cause at least ‘moderate’, ‘strong’ or ‘very strong’ shaking in the Los Angeles (LA) basin. These thresholds are used to construct a simple and robust EEW algorithm: to declare a warning, the algorithm only needs to locate the earthquake and to verify that the corresponding magnitude threshold is exceeded. The models predict that a relatively moderate M 6.5–7 earthquake along the Palos Verdes, Newport-Inglewood/Rose Canyon, Elsinore or San Jacinto faults with a rupture propagating towards LA could cause ‘very strong’ to ‘severe’ shaking in the LA basin; however, warning times for these events could exceed 30 s.

California↗