USGS Science⌕ Search

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.

At least 343 records · Page 19Linked to original sources

Earthquakes and faults in the San Francisco Bay area (1970-2003)

The map depicts both active and inactive faults and earthquakes magnitude 1.5 to 7.0 in the greater San Francisco Bay area. Twenty-two earthquakes magnitude 5.0 and greater are indicated on the map and listed chronologically in an accompanying table. The data are compiled from records from 1970-2003. The bathymetry was generated from a digital version of NOAA maps and hydrogeographic data for San Francisco Bay. Elevation data are from the USGS National Elevation Database. Landsat satellite image is from seven Landsat 7 Enhanced Thematic Mapper Plus scenes. Fault data are reproduced with permission from the California Geological Survey. The earthquake data are from the Northern California Earthquake Catalog.

California↗

History of late Holocene earthquakes at the Willow Creek site on the Nephi segment, Wasatch fault zone, Utah

This 43-page report presents new data from the Willow Creek site that provides well-defined and narrow bounds on the times of the three youngest earthquakes on the southern strand of the Nephi segment, Wasatch Fault zone, and refines the time of the youngest earthquake to about 200 years ago. This is the youngest surface rupture on the entire Wasatch fault zone, which occurred about a century or less before European settles arrived in Utah. Two trenches at the Willow Creek site exposed three scarp-derived colluvial wedges that are evidence of three paleoearthquakes. OxCal modeling of ages from Willow Creek indicate that paleoearthquake WC1 occurred at 0.2 ± 0.1 ka, WC2 occurred at 1.2 ± 0.1 ka, and WC3 occurred at 1.9 ± 0.6 ka. Stratigraphic constraints on the time of paleoearthquake WC4 are extremely poor, so OxCal modeling only yields a broadly constrained age of 4.7 ± 1.8 ka. Results from the Willow Creek site significantly refine the times of late Holocene earthquakes on the Southern strand of the Nephi segment, and this result, when combined with a reanalysis of the stratigraphic and chronologic information from previous investigations at North Creek and Red Canyon, yield a stronger basis of correlating individual earthquakes between all three sites.

Paleoseismology of Utah↗

Does paleoseismology forecast the historic rates of large earthquakes on the San Andreas fault system?

The 98-year open interval since the most recent ground-rupturing earthquake in the greater San Andreas boundary fault system would not be predicted by the quasi-periodic recurrence statistics from paleoseismic data. We examine whether the current hiatus could be explained by uncertainties in earthquake dating. Using seven independent paleoseismic records, 100 year intervals may have occurred circa 1150, 1400, and 1700 AD, but they occur in a third or less of sample records drawn at random. A second method sampling from dates conditioned on the existence of a gap of varying length suggests century-long gaps occur 3-10% of the time. A combined record with more sites would lead to lower probabilities. Systematic data over-interpretation is considered an unlikely explanation. Instead some form of non-stationary behaviour seems required, perhaps through long-range fault interaction. Earthquake occurrence since 1000 AD is not inconsistent with long-term cyclicity suggested from long runs of earthquake simulators.

Conference Paper↗

Geotechnical reconnaissance of the 2002 Denali fault, Alaska, earthquake

The 2002 M7.9 Denali fault earthquake resulted in 340 km of ruptures along three separate faults, causing widespread liquefaction in the fluvial deposits of the alpine valleys of the Alaska Range and eastern lowlands of the Tanana River. Areas affected by liquefaction are largely confined to Holocene alluvial deposits, man-made embankments, and backfills. Liquefaction damage, sparse surrounding the fault rupture in the western region, was abundant and severe on the eastern rivers: the Robertson, Slana, Tok, Chisana, Nabesna and Tanana Rivers. Synthetic seismograms from a kinematic source model suggest that the eastern region of the rupture zone had elevated strong-motion levels due to rupture directivity, supporting observations of elevated geotechnical damage. We use augered soil samples and shear-wave velocity profiles made with a portable apparatus for the spectral analysis of surface waves (SASW) to characterize soil properties and stiffness at liquefaction sites and three trans-Alaska pipeline pump station accelerometer locations. ?? 2004, Earthquake Engineering Research Institute.

Earthquake Spectra↗

Surface rupture and slip distribution of the Denali and Totschunda faults in the 3 November 2002 M 7.9 earthquake, Alaska

The 3 November 2002 Denali fault, Alaska, earthquake resulted in 341 km of surface rupture on the Susitna Glacier, Denali, and Totschunda faults. The rupture proceeded from west to east and began with a 48-km-long break on the previously unknown Susitna Glacier thrust fault. Slip on this thrust averaged about 4 m (Crone et al. , 2004) . Next came the principal surface break, along 226 km of the Denali fault, with average right-lateral offsets of 4.5–5.1 m and a maximum offset of 8.8 m near its eastern end. The Denali fault trace is commonly left stepping and north side up. About 99 km of the fault ruptured through glacier ice, where the trace orientation was commonly influenced by local ice fabric. Finally, slip transferred southeastward onto the Totschunda fault and continued for another 66 km where dextral offsets average 1.6–1.8 m. The transition from the Denali fault to the Totschunda fault occurs over a complex 25-km-long transfer zone of right-slip and normal fault traces. Three methods of calculating average surface slip all yield a moment magnitude of M w 7.8, in very good agreement with the seismologically determined magnitude of M 7.9. A comparison of strong-motion inversions for moment release with our slip distribution shows they have a similar pattern. The locations of the two largest pulses of moment release correlate with the locations of increasing steps in the average values of observed slip. This suggests that slip-distribution data can be used to infer moment release along other active fault traces.

Alaska↗

Paleoearthquakes on the southern San Andreas Fault, Wrightwood, California, 3000 to 1500 B.C.: A new method for evaluating paleoseismic evidence and earthquake horizons

We present evidence of 11–14 earthquakes that occurred between 3000 and 1500 b.c. on the San Andreas fault at the Wrightwood paleoseismic site. Earthquake evidence is presented in a novel form in which we rank (high, moderate, poor, or low) the quality of all evidence of ground deformation, which are called “event indicators.” Event indicator quality reflects our confidence that the morphologic and sedimentologic evidence can be attributable to a ground-deforming earthquake and that the earthquake horizon is accurately identified by the morphology of the feature. In four vertical meters of section exposed in ten trenches, we document 316 event indicators attributable to 32 separate stratigraphic horizons. Each stratigraphic horizon is evaluated based on the sum of rank (Rs), maximum rank (Rm), average rank (Ra), number of observations (Obs), and sum of higher-quality event indicators (Rs >1 ). Of the 32 stratigraphic horizons, 14 contain 83% of the event indicators and are qualified based on the number and quality of event indicators; the remaining 18 do not have satisfactory evidence for further consideration. Eleven of the 14 stratigraphic horizons have sufficient number and quality of event indicators to be qualified as “probable” to “very likely” earthquakes; the remaining three stratigraphic horizons are associated with somewhat ambiguous features and are qualified as “possible” earthquakes. Although no single measurement defines an obvious threshold for designation as an earthquake horizon, Rs, Rm, and Rs >1 correlate best with the interpreted earthquake quality. Earthquake age distributions are determined from radiocarbon ages of peat samples using a Bayesian approach to layer dating. The average recurrence interval for the 10 consecutive and highest-quality earthquakes is 111 (93–131) years and individual intervals are ±50% of the average. With comparison with the previously published 14–15 earthquake record between a.d. 500 and present, we find no evidence to suggest significant variations in the average recurrence rate at Wrightwood during the past 5000 years.

California↗

A grid-doubling finite-element technique for calculating dynamic three-dimensional spontaneous rupture on an earthquake fault

We present a new finite-element technique for calculating dynamic 3-D spontaneous rupture on an earthquake fault, which can reduce the required computational resources by a factor of six or more, without loss of accuracy. The grid-doubling technique employs small cells in a thin layer surrounding the fault. The remainder of the modelling volume is filled with larger cells, typically two or four times as large as the small cells. In the resulting non-conforming mesh, an interpolation method is used to join the thin layer of smaller cells to the volume of larger cells. Grid-doubling is effective because spontaneous rupture calculations typically require higher spatial resolution on and near the fault than elsewhere in the model volume. The technique can be applied to non-planar faults by morphing, or smoothly distorting, the entire mesh to produce the desired 3-D fault geometry. Using our FaultMod finite-element software, we have tested grid-doubling with both slip-weakening and rate-and-state friction laws, by running the SCEC/USGS 3-D dynamic rupture benchmark problems. We have also applied it to a model of the Hayward fault, Northern California, which uses realistic fault geometry and rock properties. FaultMod implements fault slip using common nodes, which represent motion common to both sides of the fault, and differential nodes, which represent motion of one side of the fault relative to the other side. We describe how to modify the traction-at-split-nodes method to work with common and differential nodes, using an implicit time stepping algorithm.

Geophysical Journal International↗

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↗

Large-magnitude, late Holocene earthquakes on the Genoa fault, West-Central Nevada and Eastern California

The Genoa fault, a principal normal fault of the transition zone between the Basin and Range Province and the northern Sierra Nevada, displays a large and conspicuous prehistoric scarp. Three trenches excavated across this scarp exposed two large-displacement, late Holocene events. Two of the trenches contained multiple layers of stratified charcoal, yielding radiocarbon ages suggesting the most recent and penultimate events on the main part of the fault occurred 500-600 cal B.P. , and 2000-2200 cal B.P. , respectively. Normal-slip offsets of 3-5.5 m per event along much of the rupture length are comparable to the largest historical Basin and Range Province earthquakes, suggesting these paleoearthquakes were on the order of magnitude 7.2-7.5. The apparent late Holocene slip rate (2-3 mm/yr) is one of the highest in the Basin and Range Province. Based on structural and behavioral differences, the Genoa fault is here divided into four principal sections (the Sierra, Diamond Valley, Carson Valley, and Jacks Valley sections) and is distinguished from three northeast-striking faults in the Carson City area (the Kings Canyon, Carson City, and Indian Hill faults). The conspicuous scarp extends for nearly 25 km, the combined length of the Carson Valley and Jacks Valley sections. The Diamond Valley section lacks the conspicuous scarp, and older alluvial fans and bedrock outcrops on the downthrown side of the fault indicate a lower activity rate. Activity further decreases to the south along the Sierra section, which consists of numerous distributed faults. All three northeast-striking faults in the Carson City area ruptured within the past few thousand years, and one or more may have ruptured during recent events on the Genoa fault.

California, Nevada↗

Inferring fault rheology from low-frequency earthquakes on the San Andreas

Families of recurring low-frequency earthquakes (LFEs) within nonvolcanic tremor (NVT) on the San Andreas fault in central California show strong sensitivity to shear stress induced by the daily tidal cycle. LFEs occur at all levels of the tidal shear stress and are in phase with the very small, ~400 Pa, stress amplitude. To quantitatively explain the correlation, we use a model from the existing literature that assumes the LFE sources are small, persistent regions that repeatedly fail during shear of a much larger scale, otherwise aseismically creeping fault zone. The LFE source patches see tectonic loading, creep of the surrounding fault which may be modulated by the tidal stress, and direct tidal loading. If the patches are small relative to the surrounding creeping fault then the stressing is dominated by fault creep, and if patch failure occurs at a threshold stress, then the resulting seismicity rate is proportional to the fault creep rate or fault zone strain rate. Using the seismicity rate as a proxy for strain rate and the tidal shear stress, we fit the data with possible fault rheologies that produce creep in laboratory experiments at temperatures of 400 to 600°C appropriate for the LFE source depth. The rheological properties of rock-forming minerals for dislocation creep and dislocation glide are not consistent with the observed fault creep because strong correlation between small stress perturbations and strain rate requires perturbation on the order of the ambient stress. The observed tidal modulation restricts ambient stress to be at most a few kilopascal, much lower than rock strength. A purely rate dependent friction is consistent with the observations only if the product of the friction rate dependence and effective normal stress is ~ 0.5 kPa. Extrapolating the friction rate strengthening dependence of phyllosilicates (talc) to depth would require the effective normal stress to be ~50 kPa, implying pore pressure is lithostatic. If the LFE source is on the order of tens of meters, as required by the model, rate-weakening friction rate dependence (e.g., olivine) at 400 to 600°C requires that the minimum effective pressure at the LFE source is ~ 2.5 MPa.

California↗

Holocene behavior of the Brigham City segment: implications for forecasting the next large-magnitude earthquake on the Wasatch fault zone, Utah

The Brigham City segment (BCS), the northernmost Holocene‐active segment of the Wasatch fault zone (WFZ), is considered a likely location for the next big earthquake in northern Utah. We refine the timing of the last four surface‐rupturing (~Mw 7) earthquakes at several sites near Brigham City (BE1, 2430±250; BE2, 3490±180; BE3, 4510±530; and BE4, 5610±650 cal yr B.P.) and calculate mean recurrence intervals (1060–1500 yr) that are greatly exceeded by the elapsed time (~2500 yr) since the most recent surface‐rupturing earthquake (MRE). An additional rupture observed at the Pearsons Canyon site (PC1, 1240±50 cal yr B.P.) near the southern segment boundary is probably spillover rupture from a large earthquake on the adjacent Weber segment. Our seismic moment calculations show that the PC1 rupture reduced accumulated moment on the BCS about 22%, a value that may have been enough to postpone the next large earthquake. However, our calculations suggest that the segment currently has accumulated more than twice the moment accumulated in the three previous earthquake cycles, so we suspect that additional interactions with the adjacent Weber segment contributed to the long elapse time since the MRE on the BCS. Our moment calculations indicate that the next earthquake is not only overdue, but could be larger than the previous four earthquakes. Displacement data show higher rates of latest Quaternary slip (~1.3 mm/yr) along the southern two‐thirds of the segment. The northern third likely has experienced fewer or smaller ruptures, which suggests to us that most earthquakes initiate at the southern segment boundary.

Utah↗

Catalog of earthquakes along the San Andreas fault system in central California for the year 1971

A network of seismograph stations was established by the National Center for Earthquake Research (NCER) to study in detail the earthquakes along the San Andreas fault system in central California (Eaton, Lee, and Pakiser, 1970). This report, in the form of a catalog, summarizes the results of routine earthquake locations from that network for the year 1971, and provides the basic data for further studies to be published elsewhere. Similar catalogs for 1969 and 1970 have been prepared by Lee, Roller, Bauer, and Johnson, (1972), and Lee, Roller, Meagher, and Bennett, (1972). We located a total of 2,429 events that occurred in 1971 from primarily 34,111 first P-wave arrival times recorded at 121 seismograph stations. Among these events, 339 were blasts (mostly in quarries), 8 were suspected blasts, and 41 were either too far away from available stations or too small in magnitude to be adequately located. The remaining 2,041 events which were identified as earthquakes are presented in the catalog. The 121 seismograph stations consisted of 93 telemetered stations operated by NCER, and 28 stations operated by several other organizations. Readings from the latter stations were obtained through the courtesy of the Seismographic Stations, University of California, Berkeley, the Earthquake Mechanism Laboratory, National Oceanic and Atmospheric Administration, San Francisco, and the California Department of Water Resources, Sacramento. Data on local earthquakes in northern California, Nevada, and Oregon have been published in the Bulletins of the Seismographic Stations of the University of California, Berkeley (UCB) for many decades. The present catalog covers earthquakes which occurred in a smaller area but in greater detail than those covered by the UCB Bulletins. For example, the UCB Bulletins for 1971 (Cloud and Qamar, 1972; Cloud, Tomblin, and Litehiser, 1972) list a total of 116 earthquakes (most with magnitude greater than 2.5) in the region covered by this report (a coastal strip from Santa Rosa to Parkfield). In addition, we attempted to use all available data in locating earthquakes within and near our seismic network.

California↗

An earthquake mechanism based on rapid sealing of faults

Recent seismological, heat flow and stress measurements in active fault zones such as the San Andreas have led to the suggestion 1,2 that such zones can be relatively weak. One explanation for this may be the presence of overpressured fluids along the fault 3–5 , which would reduce the shear stress required for sliding by partially 'floating' the rock. Although several mechanisms have been proposed for overpressurizing fault fluids 3,4,6,7 , we recall that 'pressure seals' are known to form in both sedimentary 8 and igneous 9 rocks by the redistribution of materials in solution; the formation of such a seal along the boundaries of a fault will prevent the communication of fluids between the porous, deforming fault zone and the surrounding country rock. Compaction of fault gouge, under hydrostatic loading and/or during shear, elevates pore pressure in the sealed fault and allows sliding at low shear stress. We report the results of laboratory sliding experiments on granite, which demonstrate that the sliding resistance of faults can be significantly decreased by sealing and compaction. The weakening that results from shear-induced compaction can be rapid, and may provide an instability mechanism for earthquakes.

Nature↗

Tilt precursors before earthquakes on the San Andreas fault, California

An array of 14 biaxial shallow-borehole tiltmeters (at 10-7 radian sensitivity) has been installed along 85 kilometers of the San Andreas fault during the past year. Earthquake-related changes in tilt have been simultaneously observed on up to four independent instruments. At earthquake distances greater than 10 earthquake source dimensions, there are few clear indications of tilt change. For the four instruments with the longest records (>10 months), 26 earthquakes have occurred since July 1973 with at least one instrument closer than 10 source dimensions and 8 earthquakes with more than one instrument within that distance. Precursors in tilt direction have been observed before more than 10 earthquakes or groups of earthquakes, and no similar effect has yet been seen without the occurrence of an earthquake.

California↗

Rupture branching structure of the 2014 Mw 6.0 South Napa, California earthquake inferred from explosion-generated fault-zone trapped waves

We present evidence for multiple fault branches of the West Napa fault zone (WNFZ) based on fault‐zone trapped waves (FZTWs) generated by two explosions that were detonated within the main surface rupture zone produced by the 24 August 2014 Mw 6.0 South Napa earthquake. The FZTWs were recorded by a 15‐kilometer‐long dense (100 m spacing) linear seismic array consisting of 155 4.5‐hertz three‐component seismometers that were deployed across the surface ruptures and adjacent faults in Napa Valley in the summer of 2016. The two explosions were located ∼3.5 km north and ∼5 km south of the 2016 recording array. Prominent FZTWs, with large amplitudes and long wavetrains following the P and S waves, are observed on the seismograms. We analyzed FZTW waveforms in both time and frequency domains to characterize the branching structure of subsurface rupture zones along the WNFZ. The 2014 surface rupture zone was ∼12 km in length along the main trace of the WNFZ, which appears to form an ∼400–600‐meter‐wide low‐velocity waveguide to depths in excess of 5–7 km. Seismic velocities within the main rupture are reduced by 40%–50% relative to the surrounding‐rock velocities. Within 1.5 km of the main trace of the WNFZ, there are at least two subordinate fault traces that formed 3‐ to 6‐kilometer‐long surface breaks during the 2014 mainshock. Our modeling suggests that these subordinate fault traces are also low‐velocity waveguides that connect with the main rupture at depths of ∼2–3 km⁠ , forming a flower structure. FZTWs were also recorded at seismic stations across the Carneros fault (CF), which is ∼1 km west of the WNFZ; this suggests that the CF connects with the WNFZ at shallow depths, even though the CF did not experience surface rupture during the 2014 Mw 6.0 mainshock. 3D finite‐difference simulations of recorded FZTWs imply a branching structure along multiple fault strands associated with the WNFZ.

California↗

A method for mapping apparent stress and energy radiation applied to the 1994 Northridge earthquake fault zone-revisited

McGarr and Fletcher (2000) introduced a technique for estimating apparent stress and seismic energy radiation associated with small patches of a larger fault plane and then applied this method to the slip model of the Northridge earthquake (Wald et al., 1996). These results must be revised because we did not take account of the difference between the seismic energy near the fault and that in the farfield. The fraction f(VR) of the near-field energy that propagates into the far-field is a monotonic function that ranges from 0.11 to 0.40 as rupture velocity VR increases from 0.6?? to 0.95??, where ?? is the shear wave speed. The revised equation for apparent stress for subfault ij is taij = f(VR) ????/ 2 Dij??? D(t)ij2dt, where ?? is density, D(t)ij is the time-dependent slip, and Dij is the final slip. The corresponding seismic energy is Eaij = ADijtaij, where A is the subfault area. Our corrected distributions of apparent stress and radiated energy over the Northridge earthquake fault zone are about 35% of those published before.

Geophysical Research Letters↗

Soufrière Hills eruption, Montserrat, 1995 - 1997: volcanic earthquake locations and fault plane solutions

A total of 9242 seismic events, recorded since the start of the eruption on Montserrat in July 1995, have been uniformly relocated with station travel-time corrections. Early seismicity was generally diffuse under southern Montserrat, and mostly restricted to depths less than 7 km. However, a NE-SW alignment of epicentres beneath the NE flank of the volcano emerged in one swarm of volcano-tectonic earthquakes (VTs) and later nests of VT hypocentres developed beneath the volcano and at a separated location, under St. George's Hill. The overall spatial distribution of hypocentres suggests a minimum depth of about 5 km for any substantial magma body. Activity associated with the opening of a conduit to the surface became increasingly shallow, with foci concentrated below the crater and, after dome building started in Fall 1995, VTs diminished and repetitive swarms of ‘hybrid’ seismic events became predominant. By late-1996, as magma effusion rates escalated, most seismic events were originating within a volume about 2 km diameter which extended up to the surface from only about 3 km depth - the diminution of shear failure earthquakes suggests the pathway for magma discharge had become effectively unconstricted. Individual and composite fault plane solutions have been determined for a few larger earthquakes. We postulate that localised extensional stress conditions near the linear VT activity, due to interaction with stresses in the overriding lithospheric plate, may encourage normal fault growth and promote sector weaknesses in the volcano.

Geophysical Research Letters↗

Slip partitioning of the Calaveras Fault, California, and prospects for future earthquakes

Examination of main shock and microearthquake data from the Calaveras Fault during the last 20 years reveals that main shock hypocenters occur at depths of 8-9 km near the base of the zone of microearthquakes. Microseismicity extends between depths of 4 and 10 km and defines zones of concentrated microseismicity and aseismic zones. Estimates of the fault regions which slipped during the Coyote Lake and Morgan Hill earthquakes as derived from seismic radiation coincide with zones which are otherwise aseismic. We propose that these persistent aseismic zones represent stuck patches which slip only during moderate earthquakes. From the pattern of microearthquake locations we recognize six aseismic zones where we expect future main shocks will rupture the Calaveras Fault. -from Authors

Journal of Geophysical Research Solid Earth↗