Geology topics

Research about San Andreas Fault

Source-linked reports with geographic coverage including San Andreas Fault.

At least 19 recordsLinked to original sources

USGS2022

A unified perspective of seismicity and fault coupling along the San Andreas Fault

The San Andreas Fault (SAF) showcases the breadth of possible earthquake sizes and occurrence behavior; in particular, the central SAF is a microcosm of such diversity. This section also exhibits the spectrum of fault coupling from locked to creeping. Here, we show that the observations of aseismic slip, temporal clustering of seismicity, and spatial variations in earthquake size distributions are tightly connected. Specifically, the creep rate along the central SAF is shown to be directly proportional to the fraction of nonclustered earthquakes for the period 1984–2020. This relationship provides a unified perspective of earthquake phenomenology along the SAF, where lower coupling manifests in weaker temporal clustering, with repeating earthquakes as an end-member. This new paradigm provides additional justification for characterizing the northwest ∼75 kilometers of the creeping segment as a transition zone, with potential implications for seismic hazard.

California↗

USGS2018

The San Andreas Fault System--Complexities along a major transform fault system and relation to earthquake hazards

The San Andreas Fault System is a 1300-km-long transform boundary that accommodates motion between the North American and Pacific Plates. New technologies and data reveal rich details about the present configuration of faults, distribution of strain and associated seismic hazard on this complex network of faults. This contribution provides a brief summary of the geologic history of the San Andreas Fault System, followed by an introduction to recent research that has changed understanding of the hazards along the main faults. Organized by region, we highlight a selection of recent research using new geodetic techniques, improved topographic data, advanced geochronologic methods, and high-resolution geophysics. The contribution ends with a review of the historic earthquakes on the San Andreas and San Jacinto Faults, comparing these to past ruptures interpreted from paleoseismic studies.

California↗

USGS2017

Evaluation of hypotheses for right-lateral displacement of Neogene strata along the San Andreas Fault between Parkfield and Maricopa, California

We used geological field studies and diatom biostratigraphy to test a published hypothesis that Neogene marine siliceous strata in the Maricopa and Parkfield areas, located on opposite sides of the San Andreas Fault, were formerly contiguous and then were displaced by about 80–130 kilometers (km) of right-lateral slip along the fault. In the Maricopa area on the northeast side of the San Andreas Fault, the upper Miocene Bitterwater Creek Shale consists of hard, siliceous shale with dolomitic concretions and turbidite sandstone interbeds. Diatom assemblages indicate that the Bitterwater Creek Shale was deposited about 8.0–6.7 million years before present (Ma) at the same time as the uppermost part of the Monterey Formation in parts of coastal California. In the Parkfield area on the southwest side of the San Andreas Fault, the upper Miocene Pancho Rico Formation consists of soft to indurated mudstone and siltstone and fossiliferous, bioturbated sandstone. Diatom assemblages from the Pancho Rico indicate deposition about 6.7–5.7 Ma (latest Miocene), younger than the Bitterwater Creek Shale and at about the same time as parts of the Sisquoc Formation and Purisima Formation in coastal California. Our results show that the Bitterwater Creek Shale and Pancho Rico Formation are lithologically unlike and of different ages and therefore do not constitute a cross-fault tie that can be used to estimate rightlateral displacement along the San Andreas Fault. In the Maricopa area northeast of the San Andreas Fault, the Bitterwater Creek Shale overlies conglomeratic fan-delta deposits of the upper Miocene Santa Margarita Formation, which in turn overlie siliceous shale of the Miocene Monterey Formation from which we obtained a diatom assemblage dated at about 10.0–9.3 Ma. Previous investigations noted that the Santa Margarita Formation in the Maricopa area contains granitic and metamorphic clasts derived from sources in the northern Gabilan Range, on the opposite side of the San Andreas Fault, that have moved relatively northwestward by 254 ± 5 km of right-lateral displacement along the fault. Our new diatom ages suggest that Santa Margarita deposition and fault displacement began about 10–8 Ma and imply long-term average slip rates along the San Andreas Fault of about 25–32 millimeters per year (mm/yr), Evaluation of Hypotheses for Right-Lateral Displacement of Neogene Strata Along the San Andreas Fault Between Parkfield and Maricopa, California By Richard G. Stanley, John A. Barron, and Charles L. Powell, II about the same as published estimates of Quaternary average slip rates based on geologic and geodetic studies.

California↗

USGS2017

A 15-year catalog of more than 1 million low-frequency earthquakes: Tracking tremor and slip along the deep San Andreas Fault

Low-frequency earthquakes (LFEs) are small, rapidly recurring slip events that occur on the deep extensions of some major faults. Their collective activation is often observed as a semi-continuous signal known as tectonic (or non-volcanic) tremor. This manuscript presents a catalog of more than 1 million LFEs detected along the central San Andreas Fault from 2001-2016. These events have been detected via a multi-channel matched-filter search, cross-correlating waveform templates representing 88 different LFE families with continuous seismic data. Together, these source locations span nearly 150 km along the central San Andreas Fault, ranging in depth from ~16-30 km. This accumulating catalog has been the source of numerous studies examining the behavior of these LFE sources and the inferred slip behavior of the deep fault. The relatively high temporal and spatial resolution of the catalog has provided new insights into properties such as tremor migration, recurrence, and triggering by static and dynamic stress perturbations. Collectively, these characteristics are inferred to reflect a very weak fault likely under near-lithostatic fluid pressure, yet the physical processes controlling the stuttering rupture observed as tremor and LFE signals remain poorly understood. This paper aims to document the LFE catalog assembly process and associated caveats, while also updating earlier observations and inferred physical constraints. The catalog itself accompanies this manuscript as part of the electronic supplement, with the goal of providing a useful resource for continued future investigations.

California↗

USGS2017

Ground-rupturing earthquakes on the northern Big Bend of the San Andreas Fault, California, 800 A.D. to Present

Paleoseismic data on the timing of ground-rupturing earthquakes constrain the recurrence behavior of active faults and can provide insight on the rupture history of a fault if earthquakes dated at neighboring sites overlap in age and are considered correlative. This study presents the evidence and ages for 11 earthquakes that occurred along the Big Bend section of the southern San Andreas Fault at the Frazier Mountain paleoseismic site. The most recent earthquake to rupture the site was the M w 7.7–7.9 Fort Tejon earthquake of 1857. We use over 30 trench excavations to document the structural and sedimentological evolution of a small pull-apart basin that has been repeatedly faulted and folded by ground-rupturing earthquakes. A sedimentation rate of 0.4 cm/yr and abundant organic material for radiocarbon dating contribute to a record that is considered complete since 800 A.D. and includes 10 paleoearthquakes. Earthquakes have ruptured this location on average every ~100 years over the last 1200 years, but individual intervals range from ~22 to 186 years. The coefficient of variation of the length of time between earthquakes (0.7) indicates quasiperiodic behavior, similar to other sites along the southern San Andreas Fault. Comparison with the earthquake chronology at neighboring sites along the fault indicates that only one other 1857-size earthquake could have occurred since 1350 A.D., and since 800 A.D., the Big Bend and Mojave sections have ruptured together at most 50% of the time in M w ≥ 7.3 earthquakes.

California↗

USGS2016

Fortnightly modulation of San Andreas tremor and low-frequency earthquakes

Earth tides modulate tremor and low-frequency earthquakes (LFEs) on faults in the vicinity of the brittle−ductile (seismic−aseismic) transition. The response to the tidal stress carries otherwise inaccessible information about fault strength and rheology. Here, we analyze the LFE response to the fortnightly tide, which modulates the amplitude of the daily tidal stress over a 14-d cycle. LFE rate is highest during the waxing fortnightly tide, with LFEs most strongly promoted when the daily stress exceeds the previous peak stress by the widest margin. This pattern implies a threshold failure process, with slip initiated when stress exceeds the local fault strength. Variations in sensitivity to the fortnightly modulation may reflect the degree of stress concentration on LFE-producing brittle asperities embedded within an otherwise aseismic fault.

California↗

USGS2016

Constraints on the source parameters of low-frequency earthquakes on the San Andreas Fault

Low-frequency earthquakes (LFEs) are small repeating earthquakes that occur in conjunction with deep slow slip. Like typical earthquakes, LFEs are thought to represent shear slip on crustal faults, but when compared to earthquakes of the same magnitude, LFEs are depleted in high-frequency content and have lower corner frequencies, implying longer duration. Here we exploit this difference to estimate the duration of LFEs on the deep San Andreas Fault (SAF). We find that the M ~ 1 LFEs have typical durations of ~0.2 s. Using the annual slip rate of the deep SAF and the average number of LFEs per year, we estimate average LFE slip rates of ~0.24 mm/s. When combined with the LFE magnitude, this number implies a stress drop of ~10 4 Pa, 2 to 3 orders of magnitude lower than ordinary earthquakes, and a rupture velocity of 0.7 km/s, 20% of the shear wave speed. Typical earthquakes are thought to have rupture velocities of ~80–90% of the shear wave speed. Together, the slow rupture velocity, low stress drops, and slow slip velocity explain why LFEs are depleted in high-frequency content relative to ordinary earthquakes and suggest that LFE sources represent areas capable of relatively higher slip speed in deep fault zones. Additionally, changes in rheology may not be required to explain both LFEs and slow slip; the same process that governs the slip speed during slow earthquakes may also limit the rupture velocity of LFEs.

California↗

USGS2015

Using a modified time-reverse imaging technique to locate low-frequency earthquakes on the San Andreas Fault near Cholame, California

We present a new method to locate low-frequency earthquakes (LFEs) within tectonic tremor episodes based on time-reverse imaging techniques. The modified time-reverse imaging technique presented here is the first method that locates individual LFEs within tremor episodes within 5 km uncertainty without relying on high-amplitude P-wave arrivals and that produces similar hypocentral locations to methods that locate events by stacking hundreds of LFEs without having to assume event co-location. In contrast to classic time-reverse imaging algorithms, we implement a modification to the method that searches for phase coherence over a short time period rather than identifying the maximum amplitude of a superpositioned wavefield. The method is independent of amplitude and can help constrain event origin time. The method uses individual LFE origin times, but does not rely on a priori information on LFE templates and families.We apply the method to locate 34 individual LFEs within tremor episodes that occur between 2010 and 2011 on the San Andreas Fault, near Cholame, California. Individual LFE location accuracies range from 2.6 to 5 km horizontally and 4.8 km vertically. Other methods that have been able to locate individual LFEs with accuracy of less than 5 km have mainly used large-amplitude events where a P-phase arrival can be identified. The method described here has the potential to locate a larger number of individual low-amplitude events with only the S-phase arrival. Location accuracy is controlled by the velocity model resolution and the wavelength of the dominant energy of the signal. Location results are also dependent on the number of stations used and are negligibly correlated with other factors such as the maximum gap in azimuthal coverage, source–station distance and signal-to-noise ratio.

California↗

USGS2015

Delayed dynamic triggering of deep tremor along the Parkfield-Cholame section of the San Andreas Fault following the 2014 M6.0 South Napa earthquake

Large, distant earthquakes are known to trigger deep tectonic tremor along the San Andreas Fault and in subduction zones. However, there are relatively few observations of triggering from regional distance earthquakes. Here we show that a small tremor episode about 12–18 km NW of Parkfield was triggered during and immediately following the passage of surface waves from the 2014 M w 6.0 South Napa main shock. More notably, a major tremor episode followed, beginning about 12 h later, and centered SE of Parkfield near Cholame. This major episode is one of the largest seen over the past several years, containing intense activity for ~3 days and taking more than 3 weeks to return to background levels. This episode showed systematic along-strike migration at ~5 km/d, suggesting that it was driven by a slow-slip event. Our results suggest that moderate-size earthquakes are capable of triggering major tremor and deep slow slip at regional distances.

California↗

USGS2015

San Andreas tremor cascades define deep fault zone complexity

Weak seismic vibrations - tectonic tremor - can be used to delineate some plate boundary faults. Tremor on the deep San Andreas Fault, located at the boundary between the Pacific and North American plates, is thought to be a passive indicator of slow fault slip. San Andreas Fault tremor migrates at up to 30 m s - 1 , but the processes regulating tremor migration are unclear. Here I use a 12-year catalogue of more than 850,000 low-frequency earthquakes to systematically analyse the high-speed migration of tremor along the San Andreas Fault. I find that tremor migrates most effectively through regions of greatest tremor production and does not propagate through regions with gaps in tremor production. I interpret the rapid tremor migration as a self-regulating cascade of seismic ruptures along the fault, which implies that tremor may be an active, rather than passive participant in the slip propagation. I also identify an isolated group of tremor sources that are offset eastwards beneath the San Andreas Fault, possibly indicative of the interface between the Monterey Microplate, a hypothesized remnant of the subducted Farallon Plate, and the North American Plate. These observations illustrate a possible link between the central San Andreas Fault and tremor-producing subduction zones.

California↗

USGS2014

A large mantle water source for the northern San Andreas Fault System: A ghost of subduction past

Recent research indicates that the shallow mantle of the Cascadia subduction margin under near-coastal Pacific Northwest U.S. is cold and partially serpentinized, storing large quantities of water in this wedge-shaped region. Such a wedge probably formed to the south in California during an earlier period of subduction. We show by numerical modeling that after subduction ceased with the creation of the San Andreas Fault System (SAFS), the mantle wedge warmed, slowly releasing its water over a period of more than 25 Ma by serpentine dehydration into the crust above. This deep, long-term water source could facilitate fault slip in San Andreas System at low shear stresses by raising pore pressures in a broad region above the wedge. Moreover, the location and breadth of the water release from this model gives insights into the position and breadth of the SAFS. Such a mantle source of water also likely plays a role in the occurrence of Non-Volcanic Tremor (NVT) that has been reported along the SAFS in central California. This process of water release from mantle depths could also mobilize mantle serpentinite from the wedge above the dehydration front, permitting upward emplacement of serpentinite bodies by faulting or by diapiric ascent. Specimens of serpentinite collected from tectonically emplaced serpentinite blocks along the SAFS show mineralogical and structural evidence of high fluid pressures during ascent from depth. Serpentinite dehydration may also lead to tectonic mobility along other plate boundaries that succeed subduction, such as other continental transforms, collision zones, or along present-day subduction zones where spreading centers are subducting.

California↗

USGS2014

Southern San Andreas Fault evaluation field activity: approaches to measuring small geomorphic offsets--challenges and recommendations for active fault studies

In southern California, where fast slip rates and sparse vegetation contribute to crisp expression of faults and microtopography, field and high‐resolution topographic data (<1 m/pixel) increasingly are used to investigate the mark left by large earthquakes on the landscape (e.g., Zielke et al., 2010; Zielke et al., 2012; Salisbury, Rockwell, et al., 2012, Madden et al., 2013). These studies measure offset streams or other geomorphic features along a stretch of a fault, analyze the offset values for concentrations or trends along strike, and infer that the common magnitudes reflect successive surface‐rupturing earthquakes along that fault section. Wallace (1968) introduced the use of such offsets, and the challenges in interpreting their “unique complex history” with offsets on the Carrizo section of the San Andreas fault; these were more fully mapped by Sieh (1978) and followed by similar field studies along other faults (e.g., Lindvall et al., 1989; McGill and Sieh, 1991). Results from such compilations spurred the development of classic fault behavior models, notably the characteristic earthquake and slip‐patch models, and thus constitute an important component of the long‐standing contrast between magnitude–frequency models (Schwartz and Coppersmith, 1984; Sieh, 1996; Hecker et al., 2013). The proliferation of offset datasets has led earthquake geologists to examine the methods and approaches for measuring these offsets, uncertainties associated with measurement of such features, and quality ranking schemes (Arrowsmith and Rockwell, 2012; Salisbury, Arrowsmith, et al., 2012; Gold et al., 2013; Madden et al., 2013). In light of this, the Southern San Andreas Fault Evaluation (SoSAFE) project at the Southern California Earthquake Center (SCEC) organized a combined field activity and workshop (the “Fieldshop”) to measure offsets, compare techniques, and explore differences in interpretation. A thorough analysis of the measurements from the field activity will be provided separately; this paper discusses the complications presented by such offset measurements using two channels from the San Andreas fault as illustrative cases. We conclude with best approaches for future data collection efforts based on input from the Fieldshop.

California↗

USGS2014

Quaternary landscape development, alluvial fan chronology and erosion of the Mecca Hills at the southern end of the San Andreas Fault zone

Quantitative geomorphic analysis combined with cosmogenic nuclide 10 Be-based geochronology and denudation rates have been used to further the understanding of the Quaternary landscape development of the Mecca Hills, a zone of transpressional uplift along the southern end of the San Andreas Fault, in southern California. The similar timing of convergent uplifts along the San Andreas Fault with the initiation of the sub-parallel San Jacinto Fault suggest a possible link between the two tectonic events. The ages of alluvial fans and the rates of catchment-wide denudation have been integrated to assess the relative influence of climate and tectonic uplift on the development of catchments within the Mecca Hills. Ages for major geomorphic surfaces based on 10 Be surface exposure dating of boulders and 10 Be depth profiles define the timing of surface stabilization to 2.6 +5.6/–1.3 ka (Qyf1 surface), 67.2 ± 5.3 ka (Qvof2 surface), and 280 ± 24 ka (Qvof1 surface). Comparison of 10 Be measurements from active channel deposits (Qac) and fluvial terraces (Qt) illustrate a complex history of erosion, sediment storage, and sediment transport in this environment. Beryllium-10 catchment-wide denudation rates range from 19.9 ± 3.2 to 149 ± 22.5 m/Ma and demonstrate strong correlations with mean catchment slope and with total active fault length normalized by catchment area. The lack of strong correlation with other geomorphic variables suggests that tectonic uplift and rock weakening have the greatest control. The currently measured topography and denudation rates across the Mecca Hills may be most consistent with a model of radial topographic growth in contrast to a model based on the rapid uplift and advection of crust.

California↗

USGS2013

S -wave triggering of tremor beneath the Parkfield, California, section of the San Andreas fault by the 2011 Tohoku, Japan earthquake: observations and theory

The dynamic stresses that are associated with the energetic seismic waves generated by the M w 9.0 Tohoku earthquake off the northeast coast of Japan triggered bursts of tectonic tremor beneath the Parkfield section of the San Andreas fault (SAF) at an epicentral distance of ∼8200 km. The onset of tremor begins midway through the ∼100‐s‐period S‐wave arrival, with a minor burst coinciding with the SHSH arrival, as recorded on the nearby broadband seismic station PKD. A more pronounced burst coincides with the Love arrival, followed by a series of impulsive tremor bursts apparently modulated by the 20‐ to 30‐s‐period Rayleigh wave. The triggered tremor was located at depths between 20 and 30 km beneath the surface trace of the fault, with the burst coincident with the S wave centered beneath the fault 30 km northwest of Parkfield. Most of the subsequent activity, including the tremor coincident with the SHSH arrival, was concentrated beneath a stretch of the fault extending from 10 to 40 km southeast of Parkfield. The seismic waves from the Tohoku epicenter form a horizontal incidence angle of ∼14°, with respect to the local strike of the SAF. Computed peak dynamic Coulomb stresses on the fault at tremor depths are in the 0.7–10 kPa range. The apparent modulation of tremor bursts by the small, strike‐parallel Rayleigh‐wave stresses (∼0.7 kPa) is likely enabled by pore pressure variations driven by the Rayleigh‐wave dilatational stress. These results are consistent with the strike‐parallel dynamic stresses (δτ s ) associated with the S, SHSH, and surface‐wave phases triggering small increments of dextral slip on the fault with a low friction (μ∼0.2). The vertical dynamic stresses δτ d do not trigger tremor with vertical or oblique slip under this simple Coulomb failure model.

California↗

USGS2013

Simulations of tremor-related creep reveal a weak crustal root of the San Andreas Fault

Deep aseismic roots of faults play a critical role in transferring tectonic loads to shallower, brittle crustal faults that rupture in large earthquakes. Yet, until the recent discovery of deep tremor and creep, direct inference of the physical properties of lower-crustal fault roots has remained elusive. Observations of tremor near Parkfield, CA provide the first evidence for present-day localized slip on the deep extension of the San Andreas Fault and triggered transient creep events. We develop numerical simulations of fault slip to show that the spatiotemporal evolution of triggered tremor near Parkfield is consistent with triggered fault creep governed by laboratory-derived friction laws between depths of 20–35 km on the fault. Simulated creep and observed tremor northwest of Parkfield nearly ceased for 20–30 days in response to small coseismic stress changes of order 10 4 Pa from the 2003 M6.5 San Simeon Earthquake. Simulated afterslip and observed tremor following the 2004 M6.0 Parkfield earthquake show a coseismically induced pulse of rapid creep and tremor lasting for 1 day followed by a longer 30 day period of sustained accelerated rates due to propagation of shallow afterslip into the lower crust. These creep responses require very low effective normal stress of ~1 MPa on the deep San Andreas Fault and near-neutral-stability frictional properties expected for gabbroic lower-crustal rock.

San Andreas Fault↗

USGS2013

Volatile ﬂuxes through the Big Bend section of the San Andreas Fault, California: helium and carbon-dioxide systematics

To investigate the source of volatiles and their relationship to the San Andreas Fault System (SAFS), 18 groundwater samples were collected from wells near the Big Bend section of the SAFS in southern California and analyzed for helium and carbon abundance and isotopes. Concentrations of 4 He, corrected for air-bubble entrainment, vary from 4.15 to 62.7 (× 10 − 8 ) cm 3 STP g − 1 H 2 O. 3 He/ 4 He ratios vary from 0.09 to 3.52 R A (where R A = air 3 He/ 4 He), consistent with up to 44% mantle helium in samples. A subset of 10 samples was analyzed for the major volatile phase (CO 2 ) — the hypothesized carrier phase of the helium in the mantle–crust system: CO 2 / 3 He ratios vary from 0.614 to 142 (× 10 11 ), and δ 13 C (CO 2 ) values vary from − 21.5 to − 11.9‰ (vs. PDB). 3 He/ 4 He ratios and CO 2 concentrations are highest in the wells located in the Mil Potrero and Cuddy valleys adjacent to the SAFS. The elevated 3 He/ 4 He ratios are interpreted to be a consequence of a mantle volatile flux though the SAFS diluted by radiogenic He produced in the crust. Samples with the highest 3 He/ 4 He ratios also had the lowest CO 2 / 3 He ratios. The combined helium isotope, He–CO 2 elemental relationships, and δ 13 C (CO 2 ) values of the groundwater volatiles reveal a mixture of mantle and deep crustal (metamorphic) fluid origins. The flux of fluids into the seismogenic zone at high hydrostatic pressure may cause fault rupture, and transfer volatiles into the shallow crust. We calculate an upward fluid flow rate of 147 mm a − 1 along the SAFS, up to 37 times higher than previous estimates (Kennedy et al., 1997). However, using newly identified characteristics of the SAFS, we calculate a total flux of 3 He along the SAFS of 7.4 × 103 cm 3 STP a − 1 (0.33 mol 3 He a − 1 ), and a CO 2 flux of 1.5 × 10 13 cm 3 STP a − 1 (6.6 × 10 8 mol a − 1 ), ~ 1% of previous estimates. Lower fluxes along the Big Bend section of the SAFS suggest that the flux of mantle volatiles alone is insufficient to cause the super hydrostatic pressure in the seismogenic zone; however, results identify crustal (metamorphic) fluids as a major component of the CO 2 volatile budget, which may represent the additional flux necessary for fault weakening pressure in the SAFS.

California↗

USGS2012

Strain‐rate profile across the Elsinore, San Jacinto, and San Andreas Faults near Palm Springs, California, 1973‐81

A profile of the strain accumulation rate along a line trending N50°E across the subparallel Elsinore, San Jacinto, and San Andreas faults near Palm Springs, California, has been constructed from trilateration surveys in the 1973‐81 interval. The strain accumulation is principally right‐lateral shear across a vertical plane parallel to fault strike (N40°W). The strain rate profile for that component exhibits two clearly resolved maxima, one centered on the San Jacinto fault (γ max = 0.35 ± 0.02 µrad/a) and the other on the San Andreas fault (γ max = 0.40 ± 0.02 µrad/a); no maximum is associated with the Elsinore fault. This result clearly implies that slip at depth on both the San Andreas and San Jacinto faults is loading the shallower sections of these faults, and that eventually rupture can be expected on both.

California↗

USGS2012

Determination of three‐dimensional velocity anomalies under a seismic array using first P arrival times from local earthquakes: 1. A homogeneous initial model

Geiger's method of locating local earthquakes has been extended to include the effect of P velocity variation along the ray paths in three dimensions. The crustal structure was modeled by rectangular blocks, and a parameter was assigned to each block describing the perturbation of P wave slowness in the block. On the basis of an initial model, a set of linear equations for the observed first P arrival times was formulated in terms of the source and medium parameters. The source parameters for all the earthquakes in the data set and medium parameters for all the blocks penetrated by the seismic rays were then determined simultaneously by the damped least squares method. A computer program has been written for the case of a homogeneous initial medium model with constant P velocity. It has been tested with artificial data and applied to a set of local earthquakes recorded by a dense seismic array in Bear Valley, California. The resultant velocity distribution in the top 5 km shows a narrow low‐velocity zone of about 5 km/s in the San Andreas fault zone sandwiched between high‐velocity regions of about 6 km/s.

California↗