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 1,279 records · Page 71Linked to original sources

USGS2022

Very low frequency earthquakes in between the seismogenic and tremor zones in Cascadia?

Megathrust earthquakes and their associated tsunamis cause some of the worst natural disasters. In addition to earthquakes, a wide range of slip behaviors are present at subduction zones, including slow earthquakes that span multiple orders of spatial and temporal scales. Understanding these events may shed light on the stress or strength conditions of the megathrust fault. Out of all types of slow earthquakes, very low frequency earthquakes (VLFEs) are most enigmatic because they are difficult to detect reliably, and the physical nature of VLFEs are poorly understood. Here we show three VLFEs in Cascadia that were dynamically triggered by a 2009 Mw 6.9 Canal de Ballenas earthquake in the Gulf of California. The VLFEs likely locate in between the seismogenic zone and the Cascadia episodic tremor and slip (ETS) zone, including one event with a moment magnitude of 5.7. This is the largest VLFE reported to date, causing clear geodetic signals. Our results show that the Cascadia megathrust fault might slip rapidly at some spots in this gap zone, and such a permissible slip behavior has direct seismic hazard implications for coastal communities and perhaps further inland. Further, the observed seismic sources may represent a new class of slip events, whose characteristics do not fit current understandings of slow or regular earthquakes.

Oregon, Washington↗

USGS1998

Non-double-couple earthquakes 2. Observations

Most studies assume that earthquakes have double-couple (DC) source mechanisms, corresponding to shear motion on planar faults. However, many well-recorded earthquakes have radiation patterns that depart radically from this model, indicating fundamentally different source processes. Seismic waves excited by advective processes, such as landslides and volcanic eruptions, are consistent with net forces rather than DCs. Some volcanic earthquakes also have single-force mechanisms, probably because of advection of magmatic fluids. Other volcanic earthquakes have mechanisms close to compensated linear vector dipoles and may be caused by magmatic intrusions. Shallow earthquakes in volcanic or geothermal areas and mines often have mechanisms with isotropic components, indicating volume changes of either explosive or implosive polarity. Such mechanisms are consistent with failure involving both shear and tensile faulting, which may be facilitated by high-pressure, high-temperature fluids. In mines, tunnels are cavities that may close. Deep-focus earthquakes occur within zones of polymorphic phase transformations in the upper mantle at depths where stick-slip instability cannot occur. Their mechanisms tend to be deviatoric (volume conserving), but non-DC, and their source processes are poorly understood. Automatic global moment tensor services routinely report statistically significant non-DC components for large earthquakes, but detailed reexamination of individual events is required to confirm such results.

Reviews of Geophysics↗

USGS2012

Energy constraints on kinematic models of oblique faulting: Loma Prieta versus Parkfield-Coalinga

The Loma Prieta earthquake was a combination of right lateral strike-slip and reverse faulting. This oblique motion has been explained by a kinematic model that has the fault accommodate both the translation and convergence required by the Pacific-North American plate motions and a local bend in the San Andreas Fault However, in other areas of California, such as Parkfield-Coalinga, this combination of motions is taken up by decoupled strike-slip and thrust faults. In this paper I develop a simple physical model, based on applying Hamilton's principle to the energy released in earthquakes, that allows us to determine if a region should undergo oblique or decoupled faulting. This determination depends on the dip of available faults and the amount of convergence relative to the amount of translation. These energy constraints correctly predict the different types of faulting at Loma Prieta and Parkfield-Coalinga.

California↗

USGS2005

Constraints on fault slip rates of the southern California plate boundary from GPS velocity and stress inversions

We use Global Positioning System (GPS) velocities and stress orientations inferred from seismicity to invert for the distribution of slip on faults in the southern California plate-boundary region. Of particular interest is how long-term slip rates are partitioned between the Indio segment of the San Andreas fault (SAF), the San Jacinto fault (SJF) and the San Bernardino segment of the SAE We use two new sets of constraints to address this problem. The first is geodetic velocities from the Southern California Earthquake Center's (SCEC) Crustal Motion Map (version 3 by Shen et al.), which includes significantly more data than previous models. The second is a regional model of stress-field orientations at seismogenic depths, as determined from earthquake focal mechanisms. While GPS data have been used in similar studies before, this is the first application of stress-field observations to this problem. We construct a simplified model of the southern California fault system, and estimate the interseismic surface velocities using a backslip approach with purely elastic strain accumulation, following Meade et al. In addition, we model the stress orientations at seismogenic depths, assuming that crustal stress results from the loading of active faults. The geodetically derived stressing rates are found to be aligned with the stress orientations from seismicity. We therefore proceed to invert simultaneously GPS and stress observations for slip rates of the faults in our network. We find that the regional patterns of crustal deformation as imaged by both data sets can be explained by our model, and that joint inversions lead to better constrained slip rates. In our preferred model, the SJF accommodates ???15 mm yr-1 and the Indio segment of the SAF ???23 mm yr-1 of right-lateral motion, accompanied by a low slip rate on the San Bernardino segment of the SAF 'Anomalous' fault segments such as around the 1992 Mw = 7.3 Landers surface rupture can be detected. There, observed stresses deviate strongly from the long-term loading as predicted by our simple model. Evaluation of model misfits together with information from palaeoseismology may provide further insights into the time dependence of strain accumulation along the San Andreas system. ?? 2004 RAS.

Geophysical Journal International↗

USGS1976

An experiment in earthquake control at Rangely, Colorado

An experiment in an oil field at Rangely, Colorado, has demonstrated the feasibility of earthquake control. Variations in seismicity were produced by controlled variations in the fluid pressure in a seismically active zone. Precise earthquake locations revealed that the earthquakes clustered about a fault trending through a zone of high pore pressure produced by secondary recovery operations. Laboratory measurements of the frictional properties of the reservoir rocks and an in situ stress measurement made near the earthquake zone were used to predict the fluid pressure required to trigger earthquakes on preexisting fractures. Fluid pressure was controlled by alternately injecting and recovering water from wells that penetrated the seismic zone. Fluid pressure was monitored in observation wells, and a computer model of the reservoir was used to infer the fluid pressure distributions in the vicinity of the injection wells. The results of this experiment confirm the predicted effect of fluid pressure on earthquake activity and indicate that earthquakes can be controlled wherever we can control the fluid pressure in a fault zone.

Colorado↗

USGS2006

San Andreas fault geometry in the Parkfield, California, region

In map view, aftershocks of the 2004 Parkfield earthquake lie along a line that forms a straighter connection between San Andreas fault segments north and south of the Parkfield reach than does the mapped trace of the fault itself. A straightedge laid on a geologic map of Central California reveals a ???50-km-long asymmetric northeastward warp in the Parkfield reach of the fault. The warp tapers gradually as it joins the straight, creeping segment of the San Andreas to the north-west, but bends abruptly across Cholame Valley at its southeast end to join the straight, locked segment that last ruptured in 1857. We speculate that the San Andreas fault surface near Parkfield has been deflected in its upper ???6 km by nonelastic behavior of upper crustal rock units. These units and the fault surface itself are warped during periods between large 1857-type earthquakes by the presence of the 1857-locked segment to the south, which buttresses intermittent coseismic and continuous aseismic slip on the Parkfield reach. Because of nonelastic behavior, the warping is not completely undone when an 1857-type event occurs, and the upper portion of the three-dimensional fault surface is slowly ratcheted into an increasingly prominent bulge. Ultimately, the fault surface probably becomes too deformed for strike-slip motion, and a new, more vertical connection to the Earth's surface takes over, perhaps along the Southwest Fracture Zone. When this happens a wedge of material currently west of the main trace will be stranded on the east side of the new main trace.

Bulletin of the Seismological Society of America↗

USGS2015

GEAR1: A Global Earthquake Activity Rate model constructed from geodetic strain rates and smoothed seismicity

Global earthquake activity rate model 1 (GEAR1) estimates the rate of shallow earthquakes with magnitudes 6–9 everywhere on Earth. It was designed to be reproducible and testable. Our preferred hybrid forecast is a log–linear blend of two parent forecasts based on the Global Centroid Moment Tensor (CMT) catalog (smoothing 4602 m≥5.767 shallow earthquakes, 1977–2004) and the Global Strain Rate Map version 2.1 (smoothing 22,415 Global Positioning System velocities), optimized to best forecast the 2005–2012 Global CMT catalog. Strain rate is a proxy for fault stress accumulation, and earthquakes indicate stress release, so a multiplicative blend is desirable, capturing the strengths of both approaches. This preferred hybrid forecast outperforms its seismicity and strain‐rate parents; the chance that this improvement stems from random seismicity fluctuations is less than 1%. The preferred hybrid is also tested against the independent parts of the International Seismological Centre‐Global Earthquake Model catalog (m≥6.8 during 1918–1976) with similar success. GEAR1 is an update of this preferred hybrid. Comparing GEAR1 to the Uniform California Earthquake Rupture Forecast Version 3 (UCERF3), net earthquake rates agree within 4% at m≥5.8 and at m≥7.0. The spatial distribution of UCERF3 epicentroids most resembles GEAR1 after UCERF3 is smoothed with a 30 km kernel. Because UCERF3 has been constructed to derive useful information from fault geometry, slip rates, paleoseismic data, and enhanced seismic catalogs (not used in our model), this is encouraging. To build parametric catastrophe bonds from GEAR1, one could calculate the magnitude for which there is a 1% (or any) annual probability of occurrence in local regions.

Bulletin of the Seismological Society of America↗

USGS2019

Impact of down-dip rupture limit and high stress drop subevents on coseismic land-level change during Cascadia megathrust earthquakes

Seismic hazard associated with Cascadia megathrust earthquakes is strongly dependent on the landward rupture extent and heterogeneous fault properties. We use 3-D numerical simulations and a seismic velocity model for Cascadia to estimate coseismic deformation due to ~M9 earthquake scenarios. Our earthquake source model is based on observations of the 2010 M8.8 Maule and 2011 M9.0 Tohoku earthquakes, which exhibited distinct strong-motion-generating subevents in the deeper portions of the fault. We compare our estimates for land-level change to paleoseismic estimates for coseismic coastal subsidence during the A.D. 1700 Cascadia earthquake. Results show that megathrust rupture extending to the 1 cm/yr locking contour provides a good match to geologic data, and along-strike variations in coastal subsidence can be produced by including strong-motion-generating subevents in the down-dip regions of the megathrust. This work demonstrates the potential to improve seismic hazard estimates for Cascadia earthquakes by comparing physics-based earthquake simulations with geologic observations.

Cascadian Fault↗

USGS2021

Diffuse deformation and surface faulting distribution from sub-metric image correlation along the 2019 Ridgecrest ruptures (California, USA)

The 2019 M w "> M w 6.4 and 7.1 Ridgecrest, California, earthquake sequence (July 2019) ruptured consecutively a system of high‐angle strike‐slip cross faults (northeast‐ and northwest‐trending) within 34 hr. The complex rupture mechanism was illuminated by seismological and geodetic data, bringing forward the issue of the interdependency of the two fault systems both at depth and at the surface, and of its effect on the final surface displacement pattern. Here, we use high‐resolution (WorldView and Pleiades) optical satellite image correlation to measure the near‐fault horizontal and vertical surface displacement fields at 0.5 m ground resolution for the two earthquakes. We point out significant differences with previous geodetic‐ and geologic‐based measurements, and document the essential role of distributed faulting and diffuse deformation in producing the observed surface displacement patterns. We derive strain fields from the horizontal displacement maps, and highlight the predominant role of rotation and shear strain in the surface rupture process. We discuss the segmentation of the rupture based on the fault geometry and along‐strike slip variations. We also image several northeast‐trending faults with similar orientation to the deeply embedded shear fabric identified in aftershock studies, and show that these cross faults are present all along the rupture, including at a scale <100 m. Finally, we compare our results to kinematic slip inversions, and show that the surface diffuse deformation is primarily associated with areas of shallow slip deficit; however, this diffuse deformation cannot be explained using elastic modeling. We conclude that inelastic processes play an important role in contributing to the total surface deformation associated with the 2019 Ridgecrest sequence.

California↗

USGS2020

Nodal seismograph recordings of the 2019 Ridgecrest Earthquake Sequence

The 2019 Ridgecrest, California earthquake sequence included M w 6.4 and Mw 7.1 earthquakes that occurred on successive days beginning on 4 July 2019. These two largest earthquakes of the sequence occurred on orthogonal faults that ruptured the Earth’s surface. To better evaluate the 3D subsurface fault structure, ( P - and S -wave) velocity, 3D and temporal variations in seismicity, and other important aspects of the earthquake sequence, we recorded aftershocks and ambient noise using up to 461 three-component nodal seismographs for about two months, beginning about one day after the M w 7.1 mainshock. The ~ 30,000 M w ≥1 earthquakes that were recorded on the dense arrays provide an unusually large volume of data with which to evaluate the earthquake sequence. This report describes the recording arrays and is intended to provide metadata for researchers interested in evaluating various aspects of the 2019 Ridgecrest earthquake sequence using the nodal data set.

California↗

USGS1985

Crustal structure of the southern Calaveras fault zone, central California, from seismic refraction investigations

A magnitude 5.7 earthquake on 6 August 1979, within the Calaveras fault zone, near Coyote Lake of west-central California, motivated a seismic-refraction investigation in this area. A northwest-southeast profile along the fault, as well as two fan profiles across the fault were recorded to examine the velocity structure of this region. The analysis of the data reveals a complicated upper crustal velocity structure with strong lateral variations in all directions. The near-surface layers consist of recent sediments with seismic P -wave velocities of 2.6 to 3.2 km/sec. These are underlain by rocks of the Great Valley Sequence which have an average velocity of 4.5 km/sec. The Great Valley Sequence is present along the whole profile; depths range from 4.3 to 4.8 km in the northwest near Anderson Lake and in the southeast of the profile line near Hollister. In the middle of the profile near Coyote Lake, however, this layer thins and we find a laterally limited higher velocity layer (5.1 km/sec) between depths of 2.8 to 4.8 km. The high-velocity zone, which coincides with a gravity high, also correlates spatially with the hypocentral area of the 6 August earthquake and its aftershocks and may therefore represent an asperity on the fault. Velocities within the fault zone were determined from the fan profiles. Near Anderson Lake, a pronounced delay of first arrivals on the fan records indicates a vertical 1- to 2-km-wide near-surface, low-velocity zone along the fault. Near Coyote Lake, the delays observed in the fan records correlate with two subsurface en-echelon fault planes which have been previously identified from lineations in the seismicity pattern. The structure of the lower crust is similar to the neighboring Diablo Range: a 8- to 9-km-thick upper crustal layer with a seismic velocity of 5.7 to 6.3 km/sec is underlain by a 3-km-thick layer with velocity 6.8 km/sec. In accordance with previous studies of the Diablo Range, there are indications of a pronounced lower crustal low-velocity zone between a depth of 17 and 23 km. The presence of this low-velocity zone suggests that higher velocity (igneous?) rocks of the middle crust have been thrust over sedimentary rocks.

California↗

USGS1994

Some observations of landslides triggered by the 29 April 1991 Racha earthquake, Republic of Georgia

On 29 April 1991 an M s 7.0 earthquake occurred in the Racha region of the Great Caucasus Mountains in north-central Republic of Georgia. The earthquake occurred on a thrust fault striking roughly east-west and dipping about 20° to 45° northward; focal depth was 17 ± 2 km. We observed no surface fault rupture, but the earthquake caused extensive structural damage to the many unreinforced stone buildings in the area, and at least 114 people were killed. Many landslides were triggered in a 2500-km 2 epicentral area, and they caused much of the structural damage and at least half the fatalities. We observed the following six types of landslides (in order of decreasing abundance): rock falls, debris slides, slumps, earth slides, rock block slides, and rock avalanches. The types of landslides triggered by the earthquake are controlled primarily by lithology and geologic structure. Enigmatic landslide processes associated with this earthquake include (1) delays of several days between earthquake shaking and significant landslide movement, probably caused by changes in groundwater conditions; (2) small co-seismic displacement of landslides active at the time of the earthquake, a possible result of viscoplastic damping of the seismic shaking; and (3) somewhat unusual failure geometries related to local topography and geologic structure.

Bulletin of the Seismological Society of America↗

USGS2004

A formulation of directivity for earthquake sources using isochrone theory

A functional form for directivity effects can be derived from isochrone theory, in which the measure of the directivity-induced amplification of an S body wave is c, the isochrone velocity. Ground displacement of the near-, intermediate-, and far-field terms of P and S waves is linear in isochrone velocity for a finite source in a whole space. We have developed an approximation c-tilde-prime of isochrone velocity that can easily be implemented as a predictor of directivity effects in empirical ground motion prediction relations. Typically, for a given fault surface, hypocenter, and site geometry, c-tilde-prime is a simple function of the hypocentral distance, the rupture distance, the crustal shear wave speed in the seismogenic zone, and the rupture velocity. c-tilde-prime typically ranges in the interval 0.44, for rupture away from the station, to about 4, for rupture toward the station. In this version of the theory directivity is independent of period. Additionally, we have created another functional form which is c-tilde-prime modified to include the approximate radiation pattern of a finite fault having a given rake. This functional form can be used to model the spatial variations of fault-parallel and fault-normal horizontal ground motions. The strengths of this formulation are 1) the proposed functional form is based on theory, 2) the predictor is unambiguously defined for all possible site locations and source rakes, and 3) it can easily be implemented for well-studied important previous earthquakes. We compare predictions of our functional form with synthetic ground motions calculated for finite strike-slip and dip-slip faults in the magnitude range 6.5 - 7.5. In general our functional form correlates best with computed fault-normal and fault-parallel motions in the synthetic motions calculated for events with M6.5. Correlation degrades but is still useful for larger events and for the geometric average horizontal motions. We have had limited success applying it to geometrically complicated faults.

Open-File Report↗

USGS2022

Earthquake scenario selection for portfolio holders in CEUS: A case study with Oklahoma DOT

Portfolio managers of spatially distributed assets in the central and eastern United States (CEUS) and other low- to moderate seismic hazard regions require scenario-based seismic risk assessment for the purpose of emergency management and planning. Uncertainties regarding the long-term seismicity of the region, unknown faults, and limited historical records complicate the selection of an earthquake scenario. Through a case study with the Oklahoma Department of Transportation (ODOT) and their portfolio of bridges, we look at one such exercise, which consists of two scenario earthquakes: one scenario earthquake selected from the U.S. Geological Survey Building Seismic Safety Commission (BSSC) scenario catalog, a magnitude (M) 7.2 event on the Meers fault, and a second aftershock selected by the consequence-driven earthquake scenario selection (Co-DESS) method. The latter is driven by ODOT’s desired service actions to be included during the earthquake drill; in this case, we identify an earthquake that is likely to trigger inspections for bridges across multiple districts, thereby testing not only inspection protocols but also coordination efforts between district groups. We find that the Co-DESS selected event is smaller in magnitude and offers different geographical options than selection through conventional selection methods, while still meeting necessary consequences for an effective earthquake exercise.

Conference Paper↗

USGS2014

Slip rates and spatially variable creep on faults of the northern San Andreas system inferred through Bayesian inversion of Global Positioning System data

Fault creep, depending on its rate and spatial extent, is thought to reduce earthquake hazard by releasing tectonic strain aseismically. We use Bayesian inversion and a newly expanded GPS data set to infer the deep slip rates below assigned locking depths on the San Andreas, Maacama, and Bartlett Springs Faults of Northern California and, for the latter two, the spatially variable interseismic creep rate above the locking depth. We estimate deep slip rates of 21.5 ± 0.5, 13.1 ± 0.8, and 7.5 ± 0.7 mm/yr below 16 km, 9 km, and 13 km on the San Andreas, Maacama, and Bartlett Springs Faults, respectively. We infer that on average the Bartlett Springs fault creeps from the Earth's surface to 13 km depth, and below 5 km the creep rate approaches the deep slip rate. This implies that microseismicity may extend below the locking depth; however, we cannot rule out the presence of locked patches in the seismogenic zone that could generate moderate earthquakes. Our estimated Maacama creep rate, while comparable to the inferred deep slip rate at the Earth's surface, decreases with depth, implying a slip deficit exists. The Maacama deep slip rate estimate, 13.1 mm/yr, exceeds long-term geologic slip rate estimates, perhaps due to distributed off-fault strain or the presence of multiple active fault strands. While our creep rate estimates are relatively insensitive to choice of model locking depth, insufficient independent information regarding locking depths is a source of epistemic uncertainty that impacts deep slip rate estimates.

California↗

USGS2014

Seismological and geodetic constraints on the 2011 M w 5.3 Trinidad, Colorado earthquake and induced deformation in the Raton Basin

The Raton Basin of southern Colorado and northern New Mexico is an actively produced hydrocarbon basin that has experienced increased seismicity since 2001, including the August 2011 M w 5.3 Trinidad normal faulting event. Following the 2011 earthquake, regional seismic observations were used to relocate 21 events, including the 2011 main shock, two foreshocks, and 13 aftershocks. Additionally, interferometric synthetic aperture radar (InSAR) observations of both the 2011 event and preevent basin deformation place constraint on the spatial kinematics of the 2011 event and localized basin subsidence due to ground water or gas withdrawal. We find that the 2011 earthquake ruptured an 8–10 km long segment of a normal fault at depths of 1.5–6.0 km within the crystalline Precambrian basement underlying the Raton Basin sedimentary rocks. The earthquake also nucleated within the crystalline basement in the vicinity of an active wastewater disposal site. The ensuing aftershock sequence demonstrated statistical properties expected for intraplate earthquakes, though the length of the 2011 earthquake is unexpectedly long for an M w 5.3 event, suggesting that wastewater disposal may have triggered a low stress drop, otherwise natural earthquake. Additionally, preevent and postevent seismicity in the Raton Basin spatially correlates to regions of subsidence observed in InSAR time series analysis. While these observations cannot discern a causal link between hydrocarbon production and seismicity, they constrain spatial relationships between active basin deformation and geological and anthropogenic features. Furthermore, the InSAR observations highlight the utility of space-based geodetic observations for monitoring and assessing anthropogenically induced and triggered deformation.

Journal of Geophysical Research B: Solid Earth↗

USGS2013

Permanently enhanced dynamic triggering probabilities as evidenced by two M ≥ 7.5 earthquakes

The 2012 M 7.7 Haida Gwaii earthquake radiated waves that likely dynamically triggered the 2013 M 7.5 Craig earthquake, setting two precedents. First, the triggered earthquake is the largest dynamically triggered shear failure event documented to date. Second, the events highlight a connection between geologic structure, sedimentary troughs that act as waveguides, and triggering probability. The Haida Gwaii earthquake excited extraordinarily large waves within and beyond the Queen Charlotte Trough, which propagated well into mainland Alaska and likely triggering the Craig earthquake along the way. Previously, focusing and associated dynamic triggering have been attributed to unpredictable source effects. This case suggests that elevated dynamic triggering probabilities may exist along the many structures where sedimentary troughs overlie major faults, such as subduction zones’ accretionary prisms and transform faults’ axial valleys. Although data are sparse, I find no evidence of accelerating seismic activity in the vicinity of the Craig rupture between it and the Haida Gwaii earthquake.

Alaska, British Columbia↗

USGS2016

Lithospheric rheology constrained from twenty-five years of postseismic deformation following the 1989 Mw 6.9 Loma Prieta earthquake

The October 17, 1989 M w 6.9 Loma Prieta earthquake provides the first opportunity of probing the crustal and upper mantle rheology in the San Francisco Bay Area since the 1906 M w 7.9 San Francisco earthquake. Here we use geodetic observations including GPS and InSAR to characterize the Loma Prieta earthquake postseismic displacements from 1989 to 2013. Pre-earthquake deformation rates are constrained by nearly 20 yr of USGS trilateration measurements and removed from the postseismic measurements prior to the analysis. We observe GPS horizontal displacements at mean rates of 1–4 mm/yr toward Loma Prieta Mountain until 2000, and ∼2 mm/yr surface subsidence of the northern Santa Cruz Mountains between 1992 and 2002 shown by InSAR, which is not associated with the seasonal and longer-term hydrological deformation in the adjoining Santa Clara Valley. Previous work indicates afterslip dominated in the early (1989–1994) postseismic period, so we focus on modeling the postseismic viscoelastic relaxation constrained by the geodetic observations after 1994. The best fitting model shows an elastic 19-km-thick upper crust above an 11-km-thick viscoelastic lower crust with viscosity of ∼6 × 10 18 Pas, underlain by a viscous upper mantle with viscosity between 3 × 1018 and 2 × 10 19 Pas. The millimeter-scale postseismic deformation does not resolve the viscosity in the different layers very well, and the lower-crustal relaxation may be localized in a narrow shear zone. However, the inferred lithospheric rheology is consistent with previous estimates based on post-1906 San Francisco earthquake measurements along the San Andreas fault system. The viscoelastic relaxation may also contribute to the enduring increase of aseismic slip and repeating earthquake activity on the San Andreas fault near San Juan Bautista, which continued for at least a decade after the Loma Prieta event.

California↗