Faults In Southern California Research Articles

Refined Earthquake Focal Mechanism Catalog for Southern California Derived With Deep Learning Algorithms

AbstractEarthquake focal mechanisms, determined with P‐wave polarities and S/P amplitude ratios, are primary data for analyzing fault zone geometry, sense of slip, and the crustal stress field. Solving for the focal mechanisms of small earthquakes is often challenging because phase arrivals and first‐motion polarities are hard to be separated from noise. To overcome this challenge, we implement convolutional‐neural‐network algorithms (Ross, Meier, & Hauksson, 2018, Ross, Meier, Hauksson, & Heaton, 2018, https://doi.org/10.1029/2017jb015251, https://doi.org/10.1785/0120180080) to detect additional phases and polarities. Using both existing and these new data, we build a high‐quality focal mechanism catalog of 297,478 events that occurred from 1981 to 2021 in southern California with the HASH method of Hardebeck and Shearer (2002), https://doi.org/10.1785/0120010200, Hardebeck and Shearer (2003), https://doi.org/10.1785/0120020236. The new focal mechanism catalog is overall consistent with the standard catalog (Yang et al., 2012, https://doi.org/10.1785/0120110311) but includes 40% more focal mechanisms, and is more consistent with moment tensor solutions derived using waveform‐fitting methods. We apply the new catalog to identify changes in focal mechanism properties caused by the occurrences of large mainshocks such as the 2010 Mw7.2 El Mayor‐Cucapah and 2019 Mw7.1 Ridgecrest earthquakes. Such changes may be associated with co‐seismic stress drops, post‐seismic deformation processes, and static stress changes on a regional scale. The new high‐resolution catalog will contribute to improved understanding of the crustal stress field, earthquake triggering mechanisms, fault zone geometry, and sense of slip on the faults in southern California.

Spatial Variations of Upper Crustal Anisotropy Along the San Jacinto Fault Zone in Southern California: Constraints From Shear Wave Splitting Analysis

AbstractTo discern spatial and explore possible existence of temporal variations of upper crustal anisotropy in an ∼15 km section of the San Jacinto Fault Zone (SJFZ) that is composed of the Buck Ridge and Clark faults in southern California, we conduct a systematic shear wave splitting investigation using local S‐wave data recorded by three broadband seismic stations located near the surface expression of the SJFZ. An automatic data selection and splitting measurement procedure is first applied, and the resulting splitting measurements are then manually screened to ensure reliability of the results. Strong spatial variations in crustal anisotropy are revealed by 1,694 pairs of splitting parameters (fast polarization orientation and splitting delay time), as reflected by the dependence of the resulting splitting parameters on the location and geometry of the raypaths. For raypaths traveling through the fault zones, the fast orientations are dominantly WNW‐ESE which is parallel to the faults and may be attributed to fluid‐filled fractures in the fault zones. For non‐fault‐zone crossing raypaths, the fast orientations are dominantly N–S which are consistent with the orientation of the regional maximum compressive stress. A three‐dimensional model of upper crustal anisotropy is constructed based on the observations. An increase in the raypath length normalized splitting times is observed after the 03/11/2013 M4.7 earthquake, which is probably attributable to changes in the spatial distribution of earthquakes before and after the M4.7 earthquake rather than reflecting temporal changes of upper crustal anisotropy.

On the Rupture Propagation of the 2019 M6.4 Searles Valley, California, Earthquake, and the Lack of Immediate Triggering of the M7.1 Ridgecrest Earthquake

AbstractThe 2019 M6.4 Searles Valley, Southern California, earthquake nucleated on a buried right‐lateral fault segment and propagated around a perpendicular fault intersection to a surface‐outcropping left‐lateral segment, but it did not propagate coseismically to the intersecting fault of the subsequent M7.1 Ridgecrest mainshock. We use the three dimensional finite element method to explore the physical reasons for this curious rupture path. Rather than model the details of the measured and inferred slip distribution, we use simple constant‐traction assumptions to explore how initial stress, hypocenter location, and the depth of burial of the initial right‐lateral segment may have influenced rupture propagation. The results suggest that only a narrow range of fault stresses, fault burial depths, and hypocenter locations would result in the observed rupture path in this earthquake. The results may have implications for interactions between faults in Southern California and beyond.

Humming Trains in Seismology: An Opportune Source for Probing the Shallow Crust

Abstract Seismologists are eagerly seeking new and preferably low-cost ways to map and track changes in the complex structure of the top few kilometers of the crust. By understanding it better, they can build on what is known regarding important, practical issues. These include telling us whether imminent earthquakes and volcanic eruptions are generating telltale underground signs of hazard, about mitigation of induced seismicity such as from deep injection of wastewater, how the Earth and its atmosphere couple, and where accessible natural resources are. Passive seismic imaging usually relies on blind correlations within extended recordings of Earth’s ceaseless “hum” or coda of well-mixed, small vibrations. In this article, we propose a complementary approach. It is seismic interferometry using opportune sources—specifically ones not stationary in time and moving in a well-understood configuration. Its interpretation relies on an accurate understanding of how these sources radiate seismic waves, precise timing, careful placement of pairs of listening stations, and seismic phase differentiation (surface and body waves). Massive freight trains were only recently recognized as such a persistent, powerful cultural (human activity-caused) seismic source. One train passage may generate a tremor with an energy output of a magnitude 1 earthquake and be detectable for up to 100 km from the track. We discuss the source mechanisms of train tremors and review the basic theory on sources. Finally, we present case studies of body- and surface-wave retrieval as an aid to mineral exploration in Canada and to monitoring of a southern California fault zone. We believe noise recovery from this new signal source, together with dense data acquisition technologies such as nodes or distributed acoustic sensing, will deeply transform our ability to monitor activity in the shallow crust at sharpened resolution in time and space.

New geodetic constraints on southern San Andreas fault-slip rates, San Gorgonio Pass, California

AbstractAssessing fault-slip rates in diffuse plate boundary systems such as the San Andreas fault in southern California is critical both to characterize seismic hazards and to understand how different fault strands work together to accommodate plate boundary motion. In places such as San Gorgonio Pass, the geometric complexity of numerous fault strands interacting in a small area adds an extra obstacle to understanding the rupture potential and behavior of each individual fault. To better understand partitioning of fault-slip rates in this region, we build a new set of elastic fault-block models that test 16 different model fault geometries for the area. These models build on previous studies by incorporating updated campaign GPS measurements from the San Bernardino Mountains and Eastern Transverse Ranges into a newly calculated GPS velocity field that has been removed of long- and short-term postseismic displacements from 12 past large-magnitude earthquakes to estimate model fault-slip rates. Using this postseismic-reduced GPS velocity field produces a best-fitting model geometry that resolves the long-standing geologic-geodetic slip-rate discrepancy in the Eastern California shear zone when off-fault deformation is taken into account, yielding a summed slip rate of 7.2 ± 2.8 mm/yr. Our models indicate that two active strands of the San Andreas system in San Gorgonio Pass are needed to produce sufficiently low geodetic dextral slip rates to match geologic observations. Lastly, results suggest that postseismic deformation may have more of a role to play in affecting the loading of faults in southern California than previously thought.

A Statistical Method for Associating Earthquakes with Their Source Faults in Southern California

ABSTRACT We present a new statistical method for associating earthquakes with their source faults in the Southern California Earthquake Center’s 3D Community Fault Models (CFMs; Plesch et al., 2007) in near-real time and for historical earthquakes. The method uses the hypocenter location, focal mechanism orientation, and earthquake sequencing to produce the probabilities of association between a given earthquake and each fault in the CFM as well as the probability that the event occurred on a fault not represented in the CFM. We used a set of known likely associations (the Known Likely Sets) as training or testing data and demonstrated that our models perform effectively on these examples and should be expected to perform well on other earthquakes with similar characteristics including the full catalog of southern California earthquakes (Hauksson et al., 2012). To produce near-real-time associations for future earthquakes, the models have been implemented as an R script and connected to the Southern California Seismic Network data processing system operated by the California Institute of Technology and the U.S. Geological Survey to automatically produce fault associations for earthquakes of M≥3.0 as they occur. To produce historical associations, we apply the method to the most recent CFM version (v.5.2), yielding modeled historical associations for all events of M≥3.0 in the catalog of southern California earthquakes from 1981 to 2016. More than 80% of these events and 99% of moment within the geography covered by the CFM had a primary association with a CFM fault. The models can help identify clusters of small earthquakes that indicate the onset of activity associated with major faults. The method will also assist in communicating objective information about the faults that source earthquakes to the scientific community and general public. In the event of a damaging southern California earthquake, the near-real-time association will provide valuable information regarding the similarity of the current event to forecast scenarios, potentially aiding in earthquake response.

Three‐Dimensional Basin and Fault Structure From a Detailed Seismic Velocity Model of Coachella Valley, Southern California

AbstractThe Coachella Valley in the northern Salton Trough is known to produce destructive earthquakes, making it a high seismic hazard area. Knowledge of the seismic velocity structure and geometry of the sedimentary basins and fault zones is required to improve earthquake hazard estimates in this region. We simultaneously inverted first P wave travel times from the Southern California Seismic Network (39,998 local earthquakes) and explosions (251 land/sea shots) from the 2011 Salton Seismic Imaging Project to obtain a 3‐D seismic velocity model. Earthquakes with focal depths ≤10 km were selected to focus on the upper crustal structure. Strong lateral velocity contrasts in the top ~3 km correlate well with the surface geology, including the low‐velocity (<5 km/s) sedimentary basin and the high‐velocity crystalline basement rocks outside the valley. Sediment thickness is ~4 km in the southeastern valley near the Salton Sea and decreases to <2 km at the northwestern end of the valley. Eastward thickening of sediments toward the San Andreas fault within the valley defines Coachella Valley basin asymmetry. In the Peninsular Ranges, zones of relatively high seismic velocities (~6.4 km/s) between 2‐ and 4‐km depth may be related to Late Cretaceous mylonite rocks or older inherited basement structures. Other high‐velocity domains exist in the model down to 9‐km depth and help define crustal heterogeneity. We identify a potential fault zone in Lost Horse Valley unassociated with mapped faults in Southern California from the combined interpretation of surface geology, seismicity, and lateral velocity changes in the model.

Absolute Stress Fields in the Source Region of the 1992 Landers Earthquake

AbstractEarthquake focal mechanisms are often inverted to obtain the deviatoric stress field. Because shear stress is equal to the frictional strength of the fault at the time of an earthquake, six components of the absolute stress tensor at the hypocenter can be obtained from a focal mechanism by combining deviatoric stress fields with the Coulomb failure criterion. For a data set of focal mechanisms determined for southern California earthquakes, including the 1992 Landers earthquake sequence, we calculated the absolute stress tensors at their hypocenters using a standard intrinsic friction coefficient under three pore pressure conditions, parameterized by the reference pore pressure at the optimally oriented faults to the stress field. Three absolute stress fields were obtained for southern California immediately before the Landers main shock by applying each data set of the stress tensors to an inversion scheme based on Bayesian statistical inference and Akaike's Bayesian Information Criterion. The coseismic stress field was calculated to obtain the absolute stress fields immediately after the main shock. The variations of the coseismic stress rotation were related to the reference pore pressure. Comparing this relation with that obtained through stress inversion, we determined the absolute stress field and the most plausible reference pore pressure to be hydrostatic. On average, the maximum shear stresses immediately before the main shock were 44 ± 15 and 79 ± 24 MPa at depths of 5 and 10 km, respectively. Earthquakes on the off‐plate boundary faults in southern California occur on faults that are loaded by Anderson‐Byerlee stress conditions.

The Palos Verdes Fault offshore Southern California: Late Pleistocene to present tectonic geomorphology, seascape evolution, and slip rate estimate based on AUV and ROV surveys

AbstractThe Palos Verdes Fault (PVF) is one of few active faults in Southern California that crosses the shoreline and can be studied using both terrestrial and subaqueous methodologies. To characterize the near‐seafloor fault morphology, tectonic influences on continental slope sedimentary processes and late Pleistocene to present slip rate, a grid of high‐resolution multibeam bathymetric data, and chirp subbottom profiles were acquired with an autonomous underwater vehicle (AUV) along the main trace of PVF in water depths between 250 and 600 m. Radiocarbon dates were obtained from vibracores collected using a remotely operated vehicle (ROV) and ship‐based gravity cores. The PVF is expressed as a well‐defined seafloor lineation marked by subtle along‐strike bends. Right‐stepping transtensional bends exert first‐order control on sediment flow dynamics and the spatial distribution of Holocene depocenters; deformed strata within a small pull‐apart basin record punctuated growth faulting associated with at least three Holocene surface ruptures. An upper (shallower) landslide scarp, a buried sedimentary mound, and a deeper scarp have been right‐laterally offset across the PVF by 55 ± 5, 52 ± 4 , and 39 ± 8 m, respectively. The ages of the upper scarp and buried mound are approximately 31 ka; the age of the deeper scarp is bracketed to 17–24 ka. These three piercing points bracket the late Pleistocene to present slip rate to 1.3–2.8 mm/yr and provide a best estimate of 1.6–1.9 mm/yr. The deformation observed along the PVF is characteristic of strike‐slip faulting and accounts for 20–30% of the total right‐lateral slip budget accommodated offshore Southern California.

Rapid mapping of ultrafine fault zone topography with structure from motion

Structure from Motion (SfM) generates high-resolution topography and coregistered texture (color) from an unstructured set of overlapping photographs taken from varying viewpoints, overcoming many of the cost, time, and logistical limitations of Light Detection and Ranging (LiDAR) and other topographic surveying methods. This paper provides the first investigation of SfM as a tool for mapping fault zone topography in areas of sparse or low-lying vegetation. First, we present a simple, affordable SfM workflow, based on an unmanned helium balloon or motorized glider, an inexpensive camera, and semiautomated software. Second, we illustrate the system at two sites on southern California faults covered by existing airborne or terrestrial LiDAR, enabling a comparative assessment of SfM topography resolution and precision. At the first site, an ∼0.1 km2 alluvial fan on the San Andreas fault, a colored point cloud of density mostly >700 points/m2 and a 3 cm digital elevation model (DEM) and orthophoto were produced from 233 photos collected ∼50 m above ground level. When a few global positioning system ground control points are incorporated, closest point vertical distances to the much sparser (∼4 points/m2) airborne LiDAR point cloud are mostly <3 cm. The second site spans an ∼1 km section of the 1992 Landers earthquake scarp. A colored point cloud of density mostly >530 points/m2 and a 2 cm DEM and orthophoto were produced from 450 photos taken from ∼60 m above ground level. Closest point vertical distances to existing terrestrial LiDAR data of comparable density are mostly <6 cm. Each SfM survey took ∼2 h to complete and several hours to generate the scene topography and texture. SfM greatly facilitates the imaging of subtle geomorphic offsets related to past earthquakes as well as rapid response mapping or long-term monitoring of faulted landscapes.

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

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

Slip rates and off‐fault deformation in Southern California inferred from GPS data and models

In Southern California, fault slip rate estimates along the San Andreas fault (SAF) and Garlock fault from geodetically constrained kinematic models are systematically at the low end or lower than geologic slip rate estimates. The sum of geodetic model slip rates across the eastern California shear zone is higher than the geologic sum. However, the ranges of reported model and geologic slip rate estimates in the literature are sufficiently large that it remains unclear whether these apparent discrepancies are real or attributable to epistemic uncertainties in the two types of estimates. We further examine uncertainties in geodetically derived slip rate estimates on major faults in Southern California by conducting a suite of inversions with four kinematic models. Long‐term rigid elastic block models constrained by the geologic slip rates cannot fit the present‐day GPS‐derived velocity field. Deforming (permanent off‐fault strain) elastic block models and viscoelastic earthquake cycle block models constrained by geologic slip rates can fit the present‐day GPS‐derived velocity field with 28–33% of the total geodetic moment rate occurring as distributed deformation off of the major faults. Models incorporating viscoelastic mantle flow predict systematically higher slip rates than purely elastic models on many of the major Southern California faults with ranges of (elastic/viscoelastic) 29–34/30–37 mm/yr for the Carrizo SAF segment, 20–24/20–32 mm/yr for the Mojave SAF segment, 14–17/18–22 mm/yr for the Coachella SAF segment, 13–19/14–22 mm/yr for the San Jacinto fault, and 5–11/5–11 mm/yr for the western Garlock fault.

Particle-size distributions of low-angle normal fault breccias: Implications for slip mechanisms on weak faults

Slip on low-angle normal faults is not well understood because they slip at high angles to the maximum principal stress directions. These faults are considered weak and their motion cannot be explained using standard Byerlee friction and Andersonian fault mechanics. One proposed mechanism for weak fault slip is reduction of effective normal stress induced by high pore-fluid pressure. This mechanism is likely to allow dilation of the fault zone and, therefore, affect the particle-size distribution of fault breccia, which has been shown to differ for unconstrained versus constrained comminution. High pore-fluid pressure can cause dilation which leads to unconstrained comminution. We analyze samples from the footwalls of two low-angle normal faults in southern California (West Salton and Whipple detachment faults) to determine the fault-rock textures and grain-size distributions (GSDs). The GSDs are fractal with fractal dimensions ranging from ∼2.6 to 3.4. The lower end of this range is thought to reflect constrained comminution and only occurs in samples from the footwall of a small-offset “minidetachment” fault about 100 m below the Whipple detachment. The higher fractal dimensions are common in cataclasites related to the main faults and also reflect constrained comminution but are overprinted by shear localization. Our GSDs are similar to those from natural and laboratory-deformed fault rocks from strong faults. We conclude that if high pore-fluid pressure aided slip on these faults, it did not strongly affect mechanisms by which brecciation occurs, implying that fluid pressure generally was sublithostatic. Independent evidence exists for lithostatic fluid pressure that having dropped or cycled to hydrostatic levelsin the minidetachment, but our GSD results suggest that periods of high fluid pressure were too short or infrequent for unconstrained comminution to have been the dominant cataclastic mechanism. Fractal dimensions of ∼2.6 for these samples suggest that little subsequent abrasion occurred due to shear localization, consistent with minor offset on the minidetachment. Main detachment footwall samples with fractal dimensions ≥3 reflect constrained comminution followed by shear-related abrasion, and suggest that seismic cycling was important in formation of main detachment cataclasites.

How do “ghost transients” from past earthquakes affect GPS slip rate estimates on southern California faults?

In this study, we investigate the extent to which viscoelastic velocity perturbations (or “ghost transients”) from individual fault segments can affect elastic block model‐based inferences of fault slip rates from GPS velocity fields. We focus on the southern California GPS velocity field, exploring the effects of known, large earthquakes for two end‐member rheological structures. Our approach is to compute, at each GPS site, the velocity perturbation relative to a cycle average for earthquake cycles on particular fault segments. We then correct the SCEC CMM4.0 velocity field for this perturbation and invert the corrected field for fault slip rates. We find that if asthenosphere viscosities are low (3 × 1018 Pa s), the current GPS velocity field is significantly perturbed by viscoelastic earthquake cycle effects associated with the San Andreas Fault segment that last ruptured in 1857 (Mw = 7.9). Correcting the GPS velocity field for this perturbation (or “ghost transient”) adds about 5 mm/a to the SAF slip rate along the Mojave and San Bernardino segments. The GPS velocity perturbations due to large earthquakes on the Garlock Fault (most recently, events in the early 1600s) and the White Wolf Fault (most recently, the Mw = 7.3 1952 Kern County earthquake) are smaller and do not influence block‐model inverted fault slip rates. This suggests that either the large discrepancy between geodetic and geologic slip rates for the Garlock Fault is not due to a ghost transient or that un‐modeled transients from recent Mojave earthquakes may influence the GPS velocity field.

Virtual California Earthquake Simulator

The model that would eventually become Virtual California started as a limited simulation model developed by Rundle (1988) for the distributed seismicity on the San Andreas and adjacent faults in southern California. This model included stress accumulation and release as well as stress interactions between faults. An updated version of this model was developed (Rundle et al. 2001, 2002, 2004), including the major strike‐slip faults in California, and was called Virtual California. Details of the model were given in Rundle et al. (2006a,b). Yakovlev et al. (2006) utilized Virtual California simulations to examine the recurrence time statistics on faults in California. They concluded that the distribution of return times on a fault is well approximated by a Weibull distribution. Yikilmaz et al. (2010) used Virtual California to simulate earthquakes on the Nankai trough, Japan, and found an excellent agreement with the historical sequence of 13 great earthquakes. Yikilmaz et al. (2011) gave a composite simulation of seismicity in northern California utilizing Virtual California for earthquakes on mapped faults, a random background of smaller earthquakes, and a branching aftershock sequence (BASS) model simulation of aftershocks. In its current incarnation, Virtual California is a sophisticated tool for simulating earthquakes on a wide variety of fault geometries in a high‐performance computing environment and is part of a larger effort by the Southern California Earthquake Center to unify the results from several different earthquake simulators. A description of this effort and a comparison of the results can be found in Tullis et al. (2012a) and Tullis et al. (2012b), respectively. There are three major components that make up Virtual California: a fault model, a set of quasi‐static interactions (Green’s …

The magnitude distribution of earthquakes near Southern California faults

We investigate seismicity near faults in the Southern California Earthquake Center Community Fault Model. We search for anomalously large events that might be signs of a characteristic earthquake distribution. We find that seismicity near major fault zones in Southern California is well modeled by a Gutenberg-Richter distribution, with no evidence of characteristic earthquakes within the resolution limits of the modern instrumental catalog. However, the b value of the locally observed magnitude distribution is found to depend on distance to the nearest mapped fault segment, which suggests that earthquakes nucleating near major faults are likely to have larger magnitudes relative to earthquakes nucleating far from major faults.

Origin of pulverized rocks during earthquake fault rupture

The origin of pulverized rocks (PR) in surface outcrops adjacent to the fault cores of the San Andreas and other major faults in Southern California is not clear, but their structural context indicates that they are clearly associated with faulting. An understanding of their origin might allow inferences to be drawn about the nature of dynamic slip on faults, including rupture mechanisms and their speed during earthquakes. In the present study, we use split Hopkinson bar recovery experiments to investigate whether PR can be produced under dynamic stress wave loading conditions in the laboratory and whether PR is diagnostic of any particular process of formation. The results of the study indicate that in Westerly granite for transition from sparse fracture to pervasive pulverization requires high strain rates in excess of 250/s and that the formation of PR may be inhibited at the larger burial depths. The constraint imposed by field observations of the relatively low strains (1-3%) in PR recovered from the field and the laboratory derived threshold for the critical strain rate (similar to 250/s and higher) together indicate that a dynamic supershear-type rupture may be necessary for the origin of pulverized rocks at distances of tens of meters away from the fault plane as observed in the field for both large strike-slip-type and the relatively small dip-slip-type fault ruptures in nature.

High-resolution imaging of the deep anisotropic structure of the San Andreas Fault system beneath southern California

SUMMARY Geological and geodetic records show that motion is distributed over a broad and complex system of faults in southern California, with many large earthquakes occurring off the main San Andreas fault. This unusually complex plate boundary geometry is imposed by the big bend of the San Andreas Fault and the dynamics of the underlying lithosphere, but lithosphere rheology and distribution of applied forces are still poorly constrained. Geodynamic modelling is thus very uncertain, and distribution of lithospheric deformations highly controversial. Here, exploiting the property that deformation orients olivine crystals, we show that deep lithospheric deformations can be mapped directly with a new imaging approach of seismic anisotropy, with a resolution fine enough for detailed geodynamic interpretation. This method relies entirely on finite-frequency effects in the splitting of SKS waves. We build an extensive data set of more than 3400 SKS splitting measurements from the analysis of SKS waves recorded by all the broad-band stations in southern California. Resolution tests demonstrate that, using this data set, we are able to resolve anisotropic structures smaller than the size of the first Fresnel zone of SKS waves. The good vertical resolution in our tomographic images results from the short inter-station spacing in southern California. In our 3-D anisotropy model, we find no evidence for localized lithospheric shear deformation beneath the San Andreas Fault. Instead, the lithospheric plate boundary is localized along a broad shear zone beneath the East California Shear Zone. Therefore, surface and deep deformation patterns are poorly correlated and most likely decoupled. Active N–S convergence in the Transverse Ranges results in strongly coherent E–W alignment of olivine fast axes in the shallow lithosphere. At the top of the asthenosphere, the observed low level of anisotropy suggests weak lithospheric basal tractions and thus a strong decoupling between the lithosphere and the asthenosphere.

Tectonic evolution of the Salton Sea inferred from seismic reflection data

The Salton Sea is an evolving pull-apart basin located between the San Andreas and Imperial faults in Southern California. Seismic and geological data reveal a rapidly subsiding southern sub-basin that is bounded by a hinge zone to the north, and northwest-dipping normal faults to the south. Oblique extension across strike-slip faults causes subsidence and leads to the formation of pull-apart basins such as the Salton Sea in southern California. The formation of these basins has generally been studied using laboratory experiments or numerical models1,2,3,4. Here we combine seismic reflection data and geological observations from the Salton Sea to understand the evolution of this nascent pull-apart basin. Our data reveal the presence of a northeast-trending hinge zone that separates the sea into northern and southern sub-basins. Differential subsidence (>10 mm yr−1) in the southern sub-basin suggests the existence of northwest-dipping basin-bounding faults near the southern shoreline, which may control the spatial distribution of young volcanism. Rotated and truncated strata north of the hinge zone suggest that the onset of extension associated with this pull-apart basin began after ∼0.5 million years ago. We suggest that slip is partitioned spatially and temporally into vertical and horizontal domains in the Salton Sea. In contrast to previous models based on historical seismicity patterns5, the rapid subsidence and fault architecture that we document in the southern part of the sea are consistent with experimental models for pull-apart basins1.

Strain rate patterns from dense GPS networks

Abstract. The knowledge of the crustal strain rate tensor provides a description of geodynamic processes such as fault strain accumulation, which is an important parameter for seismic hazard assessment, as well as anthropogenic deformation. In the past two decades, the number of observations and the accuracy of satellite based geodetic measurements like GPS greatly increased, providing measured values of displacements and velocities of points. Here we present a method to obtain the full continuous strain rate tensor from dense GPS networks. The tensorial analysis provides different aspects of deformation, such as the maximum shear strain rate, including its direction, and the dilatation strain rate. These parameters are suitable to characterize the mechanism of the current deformation. Using the velocity fields provided by SCEC and UNAVCO, we were able to localize major active faults in Southern California and to characterize them in terms of faulting mechanism. We also show that the large seismic events that occurred recently in the study region highly contaminate the measured velocity field that appears to be strongly affected by transient postseismic deformation. Finally, we applied this method to coseismic displacement data of two earthquakes in Iceland, showing that the strain fields derived by these data provide important information on the location and the focal mechanism of the ruptures.