USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Geophysical Research Earth Surface”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10Linked to original sources
USGS2005

Coseismic slip distribution of the 1923 Kanto earthquake, Japan

The slip distribution associated with the 1923 M = 7.9 Kanto, Japan, earthquake is reexamined in light of new data and modeling. We utilize a combination of first-order triangulation, second-order triangulation, and leveling data in order to constrain the coseismic deformation. The second-order triangulation data, which have not been utilized in previous studies of 1923 coseismic deformation, are associated with only slightly smaller errors than the first-order triangulation data and expand the available triangulation data set by about a factor of 10. Interpretation of these data in terms of uniform-slip models in a companion study by Nyst et al. shows that a model involving uniform coseismic slip on two distinct rupture planes explains the data very well and matches or exceeds the fit obtained by previous studies, even one which involved distributed slip. Using the geometry of the Nyst et al. two-plane slip model, we perform inversions of the same geodetic data set for distributed slip. Our preferred model of distributed slip on the Philippine Sea plate interface has a moment magnitude of 7.86. We find slip maxima of ???8-9 m beneath Odawara and ???7-8 m beneath the Miura peninsula, with a roughly 2:1 ratio of strike-slip to dip-slip motion, in agreement with a previous study. However, the Miura slip maximum is imaged as a more broadly extended feature in our study, with the high-slip region continuing from the Miura peninsula to the southern Boso peninsula region. The second-order triangulation data provide good evidence for ???3 m right-lateral strike slip on a 35-km-long splay structure occupying the volume between the upper surface of the descending Philippine Sea plate and the southern Boso peninsula. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2007

Mechanical properties of sand, silt, and clay containing tetrahydrofuran hydrate

The mechanical behavior of hydrate-bearing sediments subjected to large strains has relevance for the stability of the seafloor and submarine slopes, drilling and coring operations, and the analysis of certain small-strain properties of these sediments (for example, seismic velocities). This study reports on the results of comprehensive axial compression triaxial tests conducted at up to 1 MPa confining pressure on sand, crushed silt, precipitated silt, and clay specimens with closely controlled concentrations of synthetic hydrate. The results show that the stress-strain behavior of hydrate-bearing sediments is a complex function of particle size, confining pressure, and hydrate concentration. The mechanical properties of hydrate-bearing sediments at low hydrate concentration (probably < 40% of pore space) appear to be determined by stress-dependent soil stiffness and strength. At high hydrate concentrations (>50% of pore space), the behavior becomes more independent of stress because the hydrates control both stiffness and strength and possibly the dilative tendency of sediments by effectively increasing interparticle coordination, cementing particles together, and filling the pore space. The cementation contribution to the shear strength of hydrate-bearing sediments decreases with increasing specific surface of soil minerals. The lower the effective confining stress, the greater the impact of hydrate formation on normalized strength.

Journal of Geophysical Research B: Solid Earth↗
USGS2013

InSAR Evidence for an active shallow thrust fault beneath the city of Spokane Washington, USA

In 2001, a nearly five month long sequence of shallow, mostly small magnitude earthquakes occurred beneath the city of Spokane, a city with a population of about 200,000, in the state of Washington. During most of the sequence, the earthquakes were not well located because seismic instrumentation was sparse. Despite poor-quality locations, the earthquake hypocenters were likely very shallow, because residents near the city center both heard and felt many of the earthquakes. The combination of poor earthquake locations and a lack of known surface faults with recent movement make assessing the seismic hazards related to the earthquake swarm difficult. However, the potential for destruction from a shallow moderate-sized earthquake is high, for example Christchurch New Zealand in 2011, so assessing the hazard potential of a seismic structure involved in the Spokane earthquake sequence is important. Using interferometric synthetic aperture radar (InSAR) data from the European Space Agency ERS2 and ENVISAT satellites and the Canadian Space Agency RADARSAT-1, satellite we are able to show that slip on a shallow previously unknown thrust fault, which we name the Spokane Fault, is the source of the earthquake sequence. The part of the Spokane Fault that slipped during the 2001 earthquake sequence underlies the north part of the city, and slip on the fault was concentrated between ~0.3 and 2 km depth. Projecting the buried fault plane to the surface gives a possible surface trace for the Spokane Fault that strikes northeast from the city center into north Spokane.

Washington↗
USGS2002

Prehistoric ruptures of the Gurvan Bulag fault, Gobi Altay, Mongolia

The 1957 Gobi Altay M8.3 earthquake in southern Mongolia was associated with the simultaneous rupture of several faults, including the Gurvan Bulag reverse fault, which is located about 25 km south of the main strike-slip Bogd fault. Our study of paleoseismic excavations across the Gurvan Bulag fault suggests that the penultimate surface rupture occurred after 6.0 ka, most likely between 2.6 and 4.4 ka, and a possible earlier rupture occurred after 7.3 ka. Our interpretation of the stratigraphic relations in one of the exposures suggests that at least five earthquakes have generated surface rupture of the Gurvan Bulag fault since the abandonment of an ancient alluvial fan surface. Luminescence dating of sediment associated with this surface indicates that it formed either 26.6 ± 2.1 ka or 16.1 ± 2.0 ka. These data imply that the recurrence intervals for surface faulting on the Gurvan Bulag and Bogd faults are similar, on the order of several thousands of years, but that the penultimate surface ruptures of the two faults did not occur during the same earthquake.

Gobi Altay↗
USGS1996

Frictional behavior of large displacement experimental faults

The coefficient of friction and velocity dependence of friction of initially bare surfaces and 1-mm-thick simulated fault gouges (<90 μm) of Westerly granite were determined as a function of displacement to >400 mm at 25°C and 25 MPa normal stress. Steady state negative friction velocity dependence and a steady state fault zone microstructure are achieved after ∼18 mm displacement, and an approximately constant strength is reached after a few tens of millimeters of sliding on initially bare surfaces. Simulated fault gouges show a large but systematic variation of friction, velocity dependence of friction, dilatancy, and degree of localization with displacement. At short displacement (<10 mm), simulated gouge is strong, velocity strengthening and changes in sliding velocity are accompanied by relatively large changes in dilatancy rate. With continued displacement, simulated gouges become progressively weaker and less velocity strengthening, the velocity dependence of dilatancy rate decreases, and deformation becomes localized into a narrow basal shear which at its most localized is observed to be velocity weakening. With subsequent displacement, the fault restrengthens, returns to velocity strengthening, or to velocity neutral, the velocity dependence of dilatancy rate becomes larger, and deformation becomes distributed. Correlation of friction, velocity dependence of friction and of dilatancy rate, and degree of localization at all displacements in simulated gouge suggest that all quantities are interrelated. The observations do not distinguish the independent variables but suggest that the degree of localization is controlled by the fault strength, not by the friction velocity dependence. The friction velocity dependence and velocity dependence of dilatancy rate can be used as qualitative measures of the degree of localization in simulated gouge, in agreement with previous studies. Theory equating the friction velocity dependence of simulated gouge to the sum of the friction velocity dependence of bare surfaces and the velocity dependence of dilatancy rate of simulated gouge fails to quantitatively account for the experimental observations.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Strain accumulation across the Wasatch Fault near Ogden, Utah

Deformation of a 70 by 40‐km trilateration network spanning the north trending Wasatch fault near Ogden, Utah, has been monitored from 1972 through 1990. All but nine of the 200 measurements are consistent with deformation that is linear in time. We presume that those nine observations are contaminated by some blunder in making the measurements and that deformation is linear in time. The strain rate over the 1972–1990 interval across a 40‐km‐wide zone west of the Wasatch fault (hanging wall block) is adequately approximated by a 0.04 ± 0.01 μstrain/yr N85°E±5° uniaxial extension. Adequate coverage of the footwall block is available only over the 1981–1990 interval, and the strain rate across a 30‐km‐wide zone there is approximated by a 0.07 ± 0.03 μstrain/yr N20°E±12° uniaxial extension. The observed extension of the hanging wall block in the direction normal to fault strike suggests that the Wasatch fault may be listric; a planar fault would produce a zone of contraction normal to the fault at the surface of the hanging wall block. The component of extension parallel to fault strike observed on the footwall block is unexplained.

Utah↗
USGS1998

The coseismic slip distributions of the 1940 and 1979 Imperial Valley, California, earthquakes and their implications

Geodetic arrays observed by the U.S. Coast and Geodetic Survey span the Imperial fault in southern California. For the 1940 M 7.1 Imperial Valley earthquake, a 1934–1941 triangulation network has sufficient resolution to allow inversion for the coseismic slip distribution on fault segments 5 to 25 km long extending from the surface to a depth of 9 km. The estimated right-lateral slip is 0.8 to 1.7 m on the northern 30 km of the main trace of the Imperial fault, 4.8±0.2 m on a 10-km-long segment straddling the United States - Mexico border, and 1.3±0.4 m on a southern 25-km-long segment in Mexico. Fixing this strike-slip model and inverting 1940 leveling data only for dip slip yields 0.1 m of east-side-down dip slip. The seismic moment for this model is M 0 = (3.2±0.3) ×10 19 N m. The 1979 geodetic data set, mostly elevation changes from leveling routes, has insufficient resolution for inversion. However, it is possible to use this geodetic data set and results published by others to infer that the 1940 and 1979 earthquakes may be similar on the rupture zone common to both events. Our preferred 1940 model is similar to the 1979 geodetic results of Crook [1984] on the segments where both networks have good resolution. Elevation changes from 1940 and 1979 leveling data are very similar. Thus the geodetic data corroborate the surface slip evidence of Sharp [1982b] that the 1940 and 1979 slip distributions are examples of “characteristic slip” on the northern Imperial fault.

California↗
USGS2012

Paleoseismology along the 1980 surface rupture of the Irpinia Fault: Implications for earthquake recurrence in the southern Apennines, Italy

The Irpinia fault was the source of the M s 6.9 1980 Irpinia earthquake and produced the first unequivocal historical surface faulting in Italy. Trenching of the 1980 fault scarp at Piano di Pecore, a flat intermontane basin about 5 km south of the 1980 instrumental epicenter, provides the first data on earthquake recurrence intervals, slip per event, and slip rate on a major normal fault in the Southern Apennines fault zone. The trenches exposed evidence of four pre-1980 paleoearthquakes that occurred during the past 8600 years. A best estimate average recurrence interval is 2150 years, although the time interval between individual events varies by as much as a factor of 2. Each paleo earthquake is similar to the 1980 surface rupture in amount of slip and style of deformation, which suggests that the 1980 event is characteristic for the Irpinia fault. Slip per event values average 61 cm. The net vertical displacement of 2.12–2.36 m since 8600 cal year B.P. observed in the trenches gives a vertical slip rate of 0.25–0.35 mm/yr, a dip slip rate of 0.29–0.40 mm/yr, and an extension rate of 0.14–0.20 mm/yr. Although fault behavior data are only available for the Irpinia fault they provide a starting point for evaluating earthquake recurrence and rates of deformation in southern Apennines. They suggest that (1) fault specific earthquake recurrence intervals based on the historical seismic record overestimates the occurrence of large magnitude ( M 7) earthquakes and (2) the Holocene rate of extension across the Apennines is ≤1 mm/yr. The 1980 earthquake and the paleoseismologic observations show that repeated and localized surface faulting occurs in southern Apennines and leaves subtle but distinct geomorphic evidence that can be detected with detailed and careful investigation.

Apennine Mountains↗
USGS2012

Reflections from midcrustal rocks within the Mesozoic subduction complex near the eastern Aleutian Trench

Seismic reflection data collected in 1973 by Western Geophysical Company show that highly reflective rocks make up the midcrust of the convergent margin adjacent to the eastern Aleutian Trench. These rocks form an arch that strikes obliquely across the strongly expressed northeast-southwest structural grain of exposed Mesozoic rocks. In an earlier report we proposed that the deep events mark the location of either the presently subducting plate or the top of underplated rocks. However, the short radius of curvature of this arch precludes the possibility that the deep events indicate the igneous oceanic crust. Instead, the deep reflections could be from underplated rocks that have been arched by the imbrication or underplating of strata below the reflective rocks. High-amplitude water layer and interbed multiples prevent precise connection of shallow and deep geology. For example, the Border Ranges fault, a suture between major tectonostratigraphic terranes, is not revealed on any of our seismic sections. We speculate, however, that one band of reflections that rises toward but does not reach the surface is from the Eagle River thrust fault, which separates Late Cretaceous melange from deformed turbidite sequences of the same age.

Journal of Geophysical Research Solid Earth↗
USGS2013

Comment on “Historical perspective on seismic hazard to Hispaniola and the northeast Caribbean region” by U. ten Brink et al.

The analysis of historical earthquakes in the northeastern Caribbean by ten Brink et al. [2011, hereafter TB11] addresses the occurrence of large and destructive historical earthquakes associated with the North American-Caribbean plate boundary. One conclusion presented in TB11 is that the recurrence interval for large earthquakes on the left-lateral, strike-slip Septentrional Fault (SF) (Figure 1a) is approximately 300 years. Their Figure 7 shows rupture of the SF across the entire island of Hispaniola in CE 1200, 1542, and 1842. Our comment challenges this model for SF earthquake recurrence because it is inconsistent with our published paleoseismic data that show no large historical earthquake is associated with surface rupture along the SF east of Santiago (Figure 1a)[Prentice et al., 1993; Mann et al., 1998; Prentice et al., 2003].

Hispaniola↗
USGS1997

Resolution of fault slip along the 470-km-long rupture of the great 1906 San Francisco earthquake and its implications

Data from all available triangulation networks affected by the 1906 earthquake have been combined to assess the trade-off between slip resolution and its uncertainty and to construct a conservative image of coseismic slip along the rupture. Because of varying network aperture and station density, slip resolution is very uneven. Although slip is determined within uncertainties of ??1.0 m along 60% of the fault, constraints are poor on the remaining, mostly offshore portions of the rupture. Slip decreases from maxima of 8.6 and 7.5 m at Shelter Cove and Tomales Bay to 4.5 m near Mount Tamalpais and 2.7 m at Loma Prieta. The geodetically derived slip distribution is in poor agreement with estimates based on analysis of S wave seismograms, probably because these waves register only 20-30% of the total seismic moment obtained from longer-period surface waves. Consideration of a range of fault geometries for 1906 slip near Loma Prieta indicates right-lateral motions lie between 2.3 and 3.1 m. These values are considerably greater than the 1.5 m of measured surface slip on which several assessments of high earthquake hazard for this fault segment were based. This factor, along with the absence of 1989 slippage where 1906 surface slip was used to make the forecasts, casts doubt on some claims of success in predicting the 1989 M = 6.9 Loma Prieta earthquake.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Origin of Hawaiian tholeiite: A metasomatic model

Two voluminous magma types generated in the mantle underlying the Pacific plate are mid‐ocean ridge tholeiite (MORB) erupted at the East Pacific Rise spreading center and Hawaiian tholeiite (HT) erupted above the Hawaiian hot spot or melting anomaly. MORB has low initial 87 Sr/ 86 Sr ratios and low amounts of all incompatible trace elements including rare earths; chondrite‐normalized patterns are depleted in light rare earths. HT, by contrast, has higher initial 87 Sr/ 86 Sr and higher amounts of incompatible trace elements; chondrite‐nor‐malized patterns are enriched in the middle and light rare earths. HT is generally poorer in CaO and Al 2 O 3 and much richer in total iron and TiO 2 compared with MORB having the same MgO content. Primary magma compositions for the two volcanic systems are calculated in Fe‐Mg equilibrium with residual olivine (Fo 92 ). MORB is generated by partial melting of a trace element depleted Iherzolite source leaving a residual assemblage dominated by olivine and orthopyroxene. The percentage of partial melting for a primary magma containing 15% MgO is calculated to be 35–42% in a source mantle having a heavy rare earth content of 3×chondrite and 33–35% MgO. HT, represented by Kilauea tholeiite, is generated by partial melting of a mixture of unmelted and residual mantle for MORB which has been modified by metasomatic addition of a nephelinitic fluid, amphibole, and minor amounts of apatite and Fe‐bearing phases such as sulfide and magnetite/ilmenite. This model yields a picritic magma in equilbrium with magnesian dunite at high (>40%) degrees of partial melting. The source also has 35% MgO before partial melting. Melting in both systems in polyvariant and not controlled by lower‐temperature invariant equilibria. The low‐velocity zone is considered to be the source of metasomatic fluids that are driven upward into the lowermost lithosphere in response to a thermal plume. Picritic primary magmas are produced by shear melting, localized in the zone of thinned and metasomatized lithosphere beneath Hawaii. Melt extraction is rapid and episodic at intervals of months to decades; magma is not stored in the mantle but passes upward to a plexus of storage reservoirs located 2–6 km beneath the surface of Kilauea. Kilauea primary magmas fractionate olivine during upward transport to reach bulk compositions of 13–14% MgO in storage. Different magma batches erupted to the surface, distinguished by different major and minor element compositons compared at similar MgO content, represent combinations of differing degrees of metasomatic enrichment, differing degrees of partial melting, and some effects of premelting mantle heterogeneity.

Hawaii↗
USGS2003

Slip history and dynamic implications of the 1999 Chi-Chi, Taiwan, earthquake

[1] We investigate the rupture process of the 1999 Chi‐Chi, Taiwan, earthquake using extensive near‐source observations, including three‐component velocity waveforms at 36 strong motion stations and 119 GPS measurements. A three‐plane fault geometry derived from our previous inversion using only static data [ Ji et al. , 2001 ] is applied. The slip amplitude, rake angle, rupture initiation time, and risetime function are inverted simultaneously with a recently developed finite fault inverse method that combines a wavelet transform approach with a simulated annealing algorithm [ Ji et al. , 2002b ]. The inversion results are validated by the forward prediction of an independent data set, the teleseismic P and SH ground velocities, with notable agreement. The results show that the total seismic moment release of this earthquake is 2.7 × 10 20 N m and that most of the slip occurred in a triangular‐shaped asperity involving two fault segments, which is consistent with our previous static inversion. The rupture front propagates with an average rupture velocity of ∼2.0 km s −1 , and the average slip duration (risetime) is 7.2 s. Several interesting observations related to the temporal evolution of the Chi‐Chi earthquake are also investigated, including (1) the strong effect of the sinuous fault plane of the Chelungpu fault on spatial and temporal variations in slip history, (2) the intersection of fault 1 and fault 2 not being a strong impediment to the rupture propagation, and (3) the observation that the peak slip velocity near the surface is, in general, higher than on the deeper portion of the fault plane, as predicted by dynamic modeling.

Journal of Geophysical Research B: Solid Earth↗
USGS2002

Preeruptive inflation and surface interferometric coherence characteristics revealed by satellite radar interferometry at Makushin Volcano, Alaska: 1993-2000

Pilot reports in January 1995 and geologic field observations from the summer of 1996 indicate that a relatively small explosive eruption of Makushin, one of the more frequently active volcanoes in the Aleutian arc of Alaska, occurred on 30 January 1995. Several independent radar interferograms that each span the time period from October 1993 to September 1995 show evidence of ???7 cm of uplift centered on the volcano's east flank, which we interpret as preeruptive inflation of a ???7-km-deep magma source (??V = 0.022 km3). Subsequent interferograms for 1995-2000, a period that included no reported eruptive activity, show no evidence of additional ground deformation. Interferometric coherence at C band is found to persist for 3 years or more on lava flow and other rocky surfaces covered with short grass and sparsely distributed tall grass and for at least 1 year on most pyroclastic deposits. On lava flow and rocky surfaces with dense tall grass and on alluvium, coherence lasts for a few months. Snow and ice surfaces lose coherence within a few days. This extended timeframe of coherence over a variety of surface materials makes C band radar interferometry an effective tool for studying volcano deformation in Alaska and other similar high-latitude regions.

Journal of Geophysical Research B: Solid Earth↗
USGS2008

Observations and interpretation of fundamental mode Rayleigh wavefields recorded by the Transportable Array (USArray)

Broadband recordings of the dense Transportable Array (TA) in the western United States provide unparalleled detailed images of long-period seismic surface wavefields. With 400 stations spanning most of the western United States, wavefronts of fundamental mode Rayleigh waves may be visualized coherently across the array at periods ???40 s. In order to constrain the Rayleigh wave phase velocity structure in the western United States, I assemble a data set of vertical component seismograms from 53 teleseismic events recorded by the TA from April 2006 to October 2007. Complex amplitude spectra from these recordings at peni ods 27-100 s are interpreted using the multiplane wave tomographic method of Friederich and Wielandt (1995) and Pollitz (1999). This analysis yields detailed surface wave phase velocity and three-dimensional shear wave velocity patterns across the North American plate boundary zone, elucidating the active processes in the highly heterogeneous western U.S. upper mantle.

Journal of Geophysical Research B: Solid Earth↗
USGS2021

Geodetic constraints on a 25-year magmatic inflation episode near Three Sisters, central Oregon

Crustal inflation near the Three Sisters volcanic center documented since the mid-1990s has persisted for more than two decades. We update past analyses of the event through 2020 by simultaneously inverting InSAR interferograms, GPS time series, and leveling data for time-dependent volcanic deformation source parameters. We explore several source models to estimate how the deformation rate varied through time and to identify parameters that can reproduce measured deformation. Our preferred model is a Mogi source 4.1 km below sea level (5.9 km below the surface) about 5 km west of the summit of South Sister. Inflation started in late 1995 or 1996; the rate increased rapidly during 1998–1999, and peaked in late 1999, resulting in maximum surface uplift of about 30 cm by mid-2020. Since 2000, the inflation rate generally declined exponentially with a time constant of about 6 years. Two source inflation scenarios fit the data equally well. In the first, the crust surrounding the source is elastic and the net source-volume increase, which we attribute to persistent magma input, has been about 49 × 10 6 m 3 . The second scenario adds a viscoelastic shell surrounding the Mogi source. In that case, an injection of about 21 × 10 6 m 3 of magma prior to 2000, followed by continuing relaxation of the viscoelastic shell, can account for most of the observed surface deformation. In both scenarios, modeling reveals quasiperiodic increases in the inflation rate (pulses) with a recurrence interval of 3–4 years, both before and after 2000.

Oregon↗
USGS2000

Fracture process zone in granite

In uniaxial compression tests performed on Aue granite cores (diameter 50 mm, length 100 mm), a steel loading plate was used to induce the formation of a discrete shear fracture. A zone of distributed microcracks surrounds the tip of the propagating fracture. This process zone is imaged by locating acoustic emission events using 12 piezoceramic sensors attached to the samples. Propagation velocity of the process zone is varied by using the rate of acoustic emissions to control the applied axial force. The resulting velocities range from 2 mm/s in displacement-controlled tests to 2 μm/s in tests controlled by acoustic emission rate. Wave velocities and amplitudes are monitored during fault formation. P waves transmitted through the approaching process zone show a drop in amplitude of 26 dB, and ultrasonic velocities are reduced by 10%. The width of the process zone is ∼9 times the grain diameter inferred from acoustic data but is only 2 times the grain size from optical crack inspection. The process zone of fast propagating fractures is wider than for slow ones. The density of microcracks and acoustic emissions increases approaching the main fracture. Shear displacement scales linearly with fracture length. Fault plane solutions from acoustic events show similar orientation of nodal planes on both sides of the shear fracture. The ratio of the process zone width to the fault length in Aue granite ranges from 0.01 to 0.1 inferred from crack data and acoustic emissions, respectively. The fracture surface energy is estimated from microstructure analysis to be ∼2 J. A lower bound estimate for the energy dissipated by acoustic events is 0.1 J.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Irregular recurrence of paleoearthquakes along the central Garlock fault near El Paso Peaks, California

New investigations at the El Paso Peaks paleoseismic site have refined and extended the record of paleoearthquakes for the central Garlock fault. Event evidence, in the form of buried fissures, scarps, folds, and upward terminations, has allowed us to identify six well-resolved earthquakes, designated as events W, U, R, Q, K, and F, within the last ???7000 years. The exposures that we excavated did not reveal additional evidence for three poorly resolved events reported by McGill and Rockwell [1998]. Evidence for event Y, identified in the previous study, was reexamined and found not to be a paleoearthquake. Radiocarbon dates of detrital charcoal, combined with a method of interpolating the event ages with revised sedimentation rates, provide our preferred constraints on the timing of faulting events. The most recent surface-rupturing earthquake, event W, occurred between A.D. 1450 and 1640. The penultimate event U occurred between A.D. 675 and 950. Event R, which was unrecognized in the previous study, occurred between A.D. 250 and 475. Event Q occurred relatively shortly before event R, between A.D. 25 and 275. Event K occurred between 3340 and 2930 B.C. The oldest identified paleoearthquake, event F, occurred between 5300 and 4670 B.C. The event ages indicate that earthquake recurrence is highly irregular at the El Paso Peaks site, with individual preferred intervals ranging from as little as 215 years to as many as 3300 years. The irregular recurrence of surface-rupturing earthquakes along the central Garlock fault may be related to temporal and spatial clustering of earthquakes proposed for faults in the Eastern California Shear Zone.

Journal of Geophysical Research B: Solid Earth↗