USGS Science⌕ Search

USGS · 70030759

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

Abstract

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.

Explore related subjects

90° N90° S · 180° W ← longitude → 180° E
Source-reported bounding extent: 32.77276128128182° to 34.473495994688804° latitude; -119.51003975872379° to -115.9080695691091° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

K.M. Scharer, R.J. Weldon, T. E. Fumal, G. P. Biasi. 2007. Paleoearthquakes on the southern San Andreas Fault, Wrightwood, California, 3000 to 1500 B.C.: A new method for evaluating paleoseismic evidence and earthquake horizons. https://doi.org/10.1785/0120060137

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

ETAS-positive: Extracting unbiased epidemic-type aftershock models from incomplete earthquake catalogs

Earthquake catalogs are inherently incomplete, and short-term fluctuations in catalog completeness can strongly bias estimates of seismicity parameters. Here, we extend the positive statistics framework recently introduced for estimating Gutenberg–Richter a and b values to the Epidemic-Type Aftershock Sequence (ETAS) model. The approach focuses on complete magnitude–time intervals defined by pairs of earthquakes in which the second event is larger than the first. Comprehensive synthetic tests demonstrate that the proposed ETAS + method retrieves nearly unbiased parameter estimates even from highly incomplete catalogs, whereas parameters obtained using classical estimation approaches are strongly biased. Application to Southern California reveals a similar trend, with the productivity scaling parameter α approaching the b value, consistent with self-similar earthquake triggering within estimation uncertainties. Owing to its robustness against catalog incompleteness and computational efficiency, the ETAS + method provides a powerful new tool for seismicity analysis and earthquake forecasting.

California↗

Efficient physics‐informed ground‐motion simulations with reduced‐order models: CyberShake implications and high‐resolution site terms for southern San Andreas fault earthquakes

Recent advances in Probabilistic Seismic Hazard Analysis (PSHA) leverage physics‐based ground‐motion simulations to estimate seismic hazard, such as the CyberShake project. However, computational costs quickly escalate when performing PSHA for numerous faults or sites and can become prohibitively expensive. To reduce computational demands, CyberShake uses reciprocity and interpolates physics‐informed corrections from simulations conducted at fewer locations, but the accuracy of these interpolations remains poorly quantified. To quantify the interpolation accuracy, we derive high‐resolution, frequency‐dependent site terms for southern California and compare them with interpolated site terms using the CyberShake approach. We accomplish this by performing a set of earthquake point‐source simulations distributed along the nonplanar fault geometry for the southern San Andreas fault (SSAF) extending from Bombay Beach to Lake Hughes. Using SeisSol, we simulate three minutes of viscoelastic seismic wave propagation for these sources and store the horizontal‐component Green’s functions for 480,000 sites. We then use a scientific machine learning approach based on interpolated proper orthogonal decomposition to construct an accurate reduced‐order model of the Green’s functions to efficiently predict effective amplitude spectra (EAS) for finite‐source rupture models of SSAF earthquakes. Using minimum curvature interpolation with tension, as used in CyberShake, we compare the interpolated site terms against our high‐resolution site terms. We identify local discrepancies with EAS differing by up to a factor of approximately three. Furthermore, we identify locations where unexpectedly high or low ground motions are missed when using the interpolated dataset for these earthquakes. We estimate that our approach may be used within CyberShake to reduce the time‐to‐solution by a factor of 336 for the entire earthquake rupture forecast. Our analysis of physics‐based site terms provides more insight into the seismic hazard due to SSAF ruptures and guides future developments by combining high‐performance computing and reduced‐order modeling techniques for PSHA.

California↗