Search

The 1989 Mw 6.9 Loma Prieta earthquake is the first major event to occur along the San Andreas fault (SAF) zone in central California since the 1906 M 7.9 San Francisco earthquake. Given the complexity of this event, uncertainty has persisted as to whether this earthquake ruptured the SAF itself or a secondary fault. Recent work on the SAF in the Coachella Valley in southern California has revealed similar complexity, arising from a nonplanar, nonvertical fault geometry, and has led us to reexamine the Loma Prieta event. We have compiled data sets and data analyses in the vicinity of the Loma Prieta earthquake, including the 3D seismic velocity model and aftershock relocations of Lin and Thurber (2012), potential field data collected by the U.S. Geological Survey following the earthquake, and seismic refraction and reflection data from the 1991 profile of Catchings et al. (2004). The velocity model and aftershock relocations of Lin and Thurber (2012) reveal a geometry for the SAF that appears similar to that in the Coachella Valley (although rotated 180°): at Loma Prieta the fault dips steeply near the surface and curves with depth to join the moderately southwest-dipping main rupture below 6 km depth, itself also nonplanar. The SAF is a clear velocity boundary, with higher velocities on the northeast, attributable to Mesozoic accretionary and other rocks, and lower velocities on the southwest, attributable to Cenozoic sedimentary and volcanic rocks of the La Honda block. Rocks of the La Honda block have been offset right-laterally hundreds of kilometers from similar rocks in the southern San Joaquin Valley and vicinity, providing evidence that the curved northeast fault boundary of this block is the plate boundary. Thus, we interpret that the Loma Prieta earthquake occurred on the SAF and not on a secondary fault.

SUMMARY We derive a detailed earthquake catalogue and Vp, Vs and Vp/Vs models for the region around the 2019 Mw 6.4 and Mw7.1 Ridgecrest, California, earthquake sequence using data recorded by rapid-response, densely deployed sensors following the Ridgecrest main shock and the regional network. The new catalogue spans a 4-month period, starting on 1 June 2019, and it includes nearly 95 000 events detected and located with iterative updates to our velocity models. The final Vp and Vs models correlate well with surface geology in the top 4 km of the crust and spatial seismicity patterns at depth. Joint interpretation of the derived catalogue, velocity models, and surface geology suggests that (i) a compliant low-velocity zone near the Garlock Fault arrested the Mw 7.1 rupture at the southeast end; (ii) a stiff high-velocity zone beneath the Coso Mountains acted as a strong barrier that arrested the rupture at the northwest end and (iii) isolated seismicity on the Garlock Fault accommodated transtensional-stepover strain triggered by the main events. The derived catalogue and velocity models can be useful for multiple future studies, including further analysis of seismicity patterns, derivations of accurate source properties (e.g. focal mechanisms) and simulations of earthquake processes and radiated seismic wavefields.

The 2019 Ridgecrest, California, earthquake sequence included Mw 6.4 and 7.1 earthquakes that occurred on successive days beginning on 4 July 2019. These two largest earthquakes of the sequence occurred on orthogonal faults that ruptured the Earth’s surface. To better evaluate the 3D subsurface fault structure, (P- and S-wave) velocity, 3D and temporal variations in seismicity, and other important aspects of the earthquake sequence, we recorded aftershocks and ambient noise using up to 461 three-component nodal seismographs for about two months, beginning about one day after the Mw 7.1 mainshock. The ∼30,000Mw≥1 earthquakes that were recorded on the dense arrays provide an unusually large volume of data with which to evaluate the earthquake sequence. This report describes the recording arrays and is intended to provide metadata for researchers interested in evaluating various aspects of the 2019 Ridgecrest earthquake sequence using the nodal data set.

Surface ruptures from the 18 April 1906 M∼7.9 San Francisco earthquake were distributed over an ∼35-meter-wide zone at San Andreas Lake on the San Francisco Peninsula in California (Schussler, 1906). Since ∼1906, the surface ruptures have been largely covered by water, but with water levels at near-historic low levels in 2008–2011, we observed that the 1906 surface ruptures were no longer visible. As a fault imaging test, we acquired refraction tomography and guided-wave data across the 1906 surface ruptures in 2011. We found that individual fault traces, as mapped by Schussler (1906), can be identified on the basis of discrete low-velocity zones (VS and VP, reduced ∼40% and ∼34%, respectively) and high-amplitude guided waves. Guided waves have traditionally been observed as large-amplitude waveforms over wide (hundreds of meters to kilometers) zones of faulting, but we demonstrate that by evaluating guided waves (including Rayleigh/Love- and P/SV-types) in terms of peak ground velocity (PGV), individual near-surface fault traces within a fault zone can be precisely located, even more than 100 yr after the surface ruptures. Such precise exploration can be used to focus paleoseismic trenching efforts and to identify or exclude faulting at specific sites. We evaluated PGV of both S-wave-type and Fϕ-mode-type guided waves and found that both wave types can be used to identify subsurface fault traces. At San Andreas Lake (main fault), S-wave-type guided waves travel up to 18% slower than S body waves, and Fϕ-mode guided waves travel ∼60% slower than P body waves but ∼15% faster than S body waves. We found that guided-wave amplitudes vary with frequency but are up to five times higher than those of body waves, including the S wave. Our data are consistent with the concept that guided waves can be a strong-shaking hazard during large-magnitude earthquakes.

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

Poor knowledge of how faults slip and distribute deformation in the shallow crust hinders efforts to mitigate hazards where faults increasingly intersect with the expanding global population at Earth’s surface. Here we analyze two study sites along the 2014 M 6.0 South Napa, California, earthquake rupture, each dominated by either co- or post-seismic shallow fault slip. We combine mobile laser scanning (MLS), active-source seismic tomography, and finite element modeling to investigate how deformation rate and mechanical properties of the shallow crust affect fault behavior. Despite four orders-of-magnitude difference in the rupture velocities, MLS-derived shear strain fields are remarkably similar at the two sites and suggest deceleration of the co-seismic rupture near Earth’s surface. Constrained by the MLS and seismic data, finite element models indicate shallow faulting is more sensitive to lithologic layering and plastic yielding than to the presence of fault compliant zones (i.e., regions surrounding faults with reduced stiffness). Although both elastic and elastoplastic models can reproduce the observed surface displacement fields within the uncertainty of MLS data, elastoplastic models likely provide the most reliable representations of subsurface fault behavior, as they produce geologically reasonable stress states and are consistent with field, geodetic, and seismological observations.

We measure peak ground velocities from fault‐zone guided waves (FZGWs), generated by on‐fault earthquakes associated with the 24 August 2014 M w 6.0 South Napa earthquake. The data were recorded on three arrays deployed across north and south of the 2014 surface rupture. The observed FZGWs indicate that the West Napa fault zone (WNFZ) and the Franklin fault (FF) are continuous in the subsurface for at least 75 km. Previously published potential‐field data indicate that the WNFZ extends northward to the Maacama fault (MF), and previous geologic mapping indicates that the FF extends southward to the Calaveras fault (CF); this suggests a total length of at least 110 km for the WNFZ–FF. Because the WNFZ–FF appears contiguous with the MF and CF, these faults apparently form a continuous Calaveras–Franklin–WNFZ–Maacama (CFWM) fault that is second only in length (∼300 km) to the San Andreas fault in the San Francisco Bay area. The long distances over which we observe FZGWs, coupled with their high amplitudes (2–10 times the S waves) suggest that strong shaking from large earthquakes on any part of the CFWM fault may cause far‐field amplified fault‐zone shaking. We interpret guided waves and seismicity cross sections to indicate multiple upper crustal splays of the WNFZ–FF, including a northward extension of the Southhampton fault, which may cause strong shaking in the Napa Valley and the Vallejo area. Based on travel times from each earthquake to each recording array, we estimate average P ‐, S ‐, and guided‐wave velocities within the WNFZ–FF (4.8–5.7, 2.2–3.2, and 1.1–2.8 km/s, respectively), with FZGW velocities ranging from 58% to 93% of the average S ‐wave velocities.

The aftershocks of the 24 August 2014 M w 6.0 South Napa earthquake generated prominent fault‐zone trapped waves (FZTWs) that were recorded on two 1.9‐km‐long seismic arrays deployed across the northern projection (array 1, A1) and the southern part (A2) of the surface rupture of the West Napa fault zone (WNFZ). We also observed FZTWs on an array (A3) deployed across the intersection of the Franklin and Southampton faults, which appear to be the southward continuations of the WNFZ. A1, A2, and A3 consisted of 20, 20, and 10 L28 (4.5 Hz) three‐component seismographs. We analyzed waveforms of FZTWs from 55 aftershocks in both time and frequency to characterize the fault damage zone associated with this M w 6.0 earthquake. Post‐ S coda durations of FZTWs increase with epicentral distances and focal depths from the recording arrays, suggesting a low‐velocity waveguide along the WNFZ to depths in excess of 5–7 km. Locations of the aftershocks showing FZTWs, combined with 3D finite‐difference simulations, suggest the subsurface rupture zone having an S ‐wave speed reduction of ∼40%–50% between A1 and A2, coincident with the ∼14‐km‐long mapped surface rupture zone and at least an ∼500‐m‐wide deformation zone. The low‐velocity waveguide along the WNFZ extends further southward to at least A3, but with a more moderate‐velocity reduction of 30%–35% at ray depth. This last FZTW observation suggests continuity between the WNFZ and Franklin fault. The waveguide effect may have localized and amplified ground shaking along the WNFZ and the faults farther to the south (see a companion paper by Catchings et al. , 2016).

The Koyna-Warna area of India is one of the best worldwide examples of reservoir-induced seismicity, with the distinction of having generated the largest known induced earthquake (M6.3 on 10 December 1967) and persistent moderate-magnitude (>M5) events for nearly 50 years. Yet, the fault structure and tectonic setting that has accommodated the induced seismicity is poorly known, in part because the seismic events occur beneath a thick sequence of basalt layers. On the basis of the alignment of earthquake epicenters over an ~50 year period, lateral variations in focal mechanisms, upper-crustal tomographic velocity images, geophysical data (aeromagnetic, gravity, and magnetotelluric), geomorphic data, and correlation with similar structures elsewhere, we suggest that the Koyna-Warna area lies within a right step between northwest trending, right-lateral faults. The sub-basalt basement may form a local structural depression (pull-apart basin) caused by extension within the step-over zone between the right-lateral faults. Our postulated model accounts for the observed pattern of normal faulting in a region that is dominated by north-south directed compression. The right-lateral faults extend well beyond the immediate Koyna-Warna area, possibly suggesting a more extensive zone of seismic hazards for the central India area. Induced seismic events have been observed many places worldwide, but relatively large-magnitude induced events are less common because critically stressed, preexisting structures are a necessary component. We suggest that releasing bends and fault step-overs like those we postulate for the Koyna-Warna area may serve as an ideal tectonic environment for generating moderate- to large- magnitude induced (reservoir, injection, etc.) earthquakes.

The Salton Sea Geothermal Field is one of the most geothermally and seismically active areas in California and presents an opportunity to study the effect of high-temperature metamorphism on the properties of seismogenic faults. The area includes numerous active tectonic faults that have recently been imaged with active source seismic reflection and refraction. We utilize the active source surveys, along with the abundant microseismicity data from a dense borehole seismic network, to image the 3-D variations in seismic velocity in the upper 5 km of the crust. There are strong velocity variations, up to ~30%, that correlate spatially with the distribution of shallow heat flow patterns. The combination of hydrothermal circulation and high-temperature contact metamorphism has significantly altered the shallow sandstone sedimentary layers within the geothermal field to denser, more feldspathic, rock with higher P wave velocity, as is seen in the numerous exploration wells within the field. This alteration appears to have a first-order effect on the frictional stability of shallow faults. In 2005, a large earthquake swarm and deformation event occurred. Analysis of interferometric synthetic aperture radar data and earthquake relocations indicates that the shallow aseismic fault creep that occurred in 2005 was localized on the Kalin fault system that lies just outside the region of high-temperature metamorphism. In contrast, the earthquake swarm, which includes all of the M > 4 earthquakes to have occurred within the Salton Sea Geothermal Field in the last 15 years, ruptured the Main Central Fault (MCF) system that is localized in the heart of the geothermal anomaly. The background microseismicity induced by the geothermal operations is also concentrated in the high-temperature regions in the vicinity of operational wells. However, while this microseismicity occurs over a few kilometer scale region, much of it is clustered in earthquake swarms that last from hours to a few days and are localized near the MCF system.

The San Andreas fault (SAF) is one of the most studied strike‐slip faults in the world; yet its subsurface geometry is still uncertain in most locations. The Salton Seismic Imaging Project (SSIP) was undertaken to image the structure surrounding the SAF and also its subsurface geometry. We present SSIP studies at two locations in the Coachella Valley of the northern Salton trough. On our line 4, a fault‐crossing profile just north of the Salton Sea, sedimentary basin depth reaches 4 km southwest of the SAF. On our line 6, a fault‐crossing profile at the north end of the Coachella Valley, sedimentary basin depth is ∼2–3 km and centered on the central, most active trace of the SAF. Subsurface geometry of the SAF and nearby faults along these two lines is determined using a new method of seismic‐reflection imaging, combined with potential‐field studies and earthquakes. Below a 6–9 km depth range, the SAF dips ∼50°–60° NE, and above this depth range it dips more steeply. Nearby faults are also imaged in the upper 10 km, many of which dip steeply and project to mapped surface fault traces. These secondary faults may join the SAF at depths below about 10 km to form a flower‐like structure. In Appendix D, we show that rupture on a northeast‐dipping SAF, using a single plane that approximates the two dips seen in our study, produces shaking that differs from shaking calculated for the Great California ShakeOut, for which the southern SAF was modeled as vertical in most places: shorter‐period (T

Although seismicity near Koyna Reservoir (India) has persisted for ~50 years and includes the largest induced earthquake (M 6.3) reported worldwide, the seismotectonic framework of the area is not well understood. We recorded ~1800 earthquakes from 6 January 2010 to 28 May 2010 and located a subset of 343 of the highest-quality earthquakes using the tomoDD code of Zhang and Thurber (2003) to better understand the framework. We also inverted first arrivals for 3-D Vp, Vs, and Vp/Vs and Poisson's ratio tomography models of the upper 12 km of the crust. Epicenters for the recorded earthquakes are located south of the Koyna River, including a high-density cluster that coincides with a shallow depth (

The determination of near‐surface (vadose zone and slightly below) fault locations and geometries is important because assessment of ground rupture, strong shaking, geologic slip rates, and rupture histories occurs at shallow depths. However, seismic imaging of fault zones at shallow depths can be difficult due to near‐surface complexities, such as weathering, groundwater saturation, massive (nonlayered) rocks, and vertically layered strata. Combined P ‐ and S ‐wave seismic‐refraction tomography data can overcome many of the near‐surface, fault‐zone seismic‐imaging problems because of differences in the responses of elastic (bulk and shear) moduli of P and S waves to shallow‐depth, fault‐zone properties. We show that high‐resolution refraction tomography images of P ‐ to S ‐wave velocity ratios ( V P / V S ) can reliably identify near‐surface faults. We demonstrate this method using tomography images of the San Andreas fault (SAF) surface‐rupture zone associated with the 18 April 1906 ∼ M 7.9 San Francisco earthquake on the San Francisco peninsula in California. There, the SAF cuts through Franciscan melange, which consists of an incoherent assemblage of greywacke, chert, greenstone, and serpentinite. A near‐vertical zone (∼75° northeast dip) of high P ‐wave velocities (up to 3000 m/s), low S ‐wave velocities (∼150–600 m/s), high V P / V S ratios (4–8.8), and high Poisson’s ratios (0.44–0.49) characterizes the main surface‐rupture zone to a depth of about 20 m and is consistent with nearby trench observations. We suggest that the combined V P / V S imaging approach can reliably identify most near‐surface fault zones in locations where many other seismic methods cannot be applied.

In a comment on our 2008 paper (Catchings, Gandhok, et al. , 2008) on the Santa Monica fault in Los Angeles, California, Pratt and Dolan (2010) (herein referred to as P&D) cite numerous objections to our work, inferring that our study is flawed. However, as shown in our reply, their objections contradict their own published works, published works of others, and proven seismic methodologies. Rather than responding to each repeated invalid objection, we address their objections by topic in the subsequent sections. In Catchings, Gandhok, et al. (2008), we presented high-resolution seismic-reflection images that showed two near-surface faults in the upper 50 m beneath the grounds of the Wadsworth Veterans Administration Hospital (WVAH). Although P&D suggest we effectively duplicated their seismic acquisition, our survey was not a duplication of their efforts. Rather, we conducted a seismic-imaging survey over a similar profile as Pratt et al. (1998) but used a different data acquisition system and different data processing methods to evaluate methods of seismically imaging blind faults in the wake of the 17 January 1994 M 6.7 Northridge earthquake. We used an acquisition method that provides both tomographic seismic velocities and reflection images. Our combined-data approach allowed for shallower imaging (∼2.5 m minimum) than the ∼20-m minimum of Pratt et al. (1998), clearer images of the fault zone, and more accurate depth determinations (rather than time images). In processing the reflection images, we used prestack depth migration, which is generally accepted as the only proper imaging method for imaging subsurface structures with strong lateral velocity variations (Versteeg, 1993), a condition shown to exist at the WVAH site. We correlated our reflection images with refraction tomography images, borehole lithology, and velocity data, Interferometric Synthetic Aperture Radar images, and changes in groundwater depths. Except for some minor differences, our seismic-reflection images coincide with previously published seismic-reflection …

The Mission Creek and Banning faults are two of the principal strands of the San Andreas fault zone in the northern Coachella Valley of southern California. Structural characteristics of the faults affect both regional earthquake hazards and local groundwater resources. We use seismic, gravity, and geological data to characterize the San Andreas fault zone in the vicinity of Desert Hot Springs. Seismic images of the upper 500 m of the Mission Creek fault at Desert Hot Springs show multiple fault strands distributed over a 500 m wide zone, with concentrated faulting within a central 200 m wide area of the fault zone. High-velocity (up to 5000 m/sec) rocks on the northeast side of the fault are juxtaposed against a low-velocity ( M >6.0) earthquakes in the area (in 1948 and 1986) occurred at or near the depths (∼10 to 12 km) of the merged (San Andreas) fault. Large-magnitude earthquakes that nucleate at or below the merged fault will likely generate strong shaking from guided waves along both fault zones and from amplified seismic waves in the low-velocity basin between the two fault zones. The Mission Creek fault zone is a groundwater barrier with the top of the water table varying by 60 m in depth and the aquifer varying by about 50 m in thickness across a 200 m wide zone of concentrated faulting.

To better understand the structure of the San Andreas fault (SAF) at Burro Flats in southern California, we acquired a three-dimensional combined set of seismic reflection and refraction profiles centered on the main active trace at Burro Flats. In this article, we discuss the variation in shallow-depth velocities along each seismic profile, with special emphasis on the 1500 m/sec P -wave velocity contour, which can be an indicator of shallow-depth water-saturated unconsolidated sediments. Along the four seismic profiles, minimum depths of the groundwater table, as inferred from 1500 m/sec velocity contour, range from 10 to about 20 m. The largest variations in depth to the top of the groundwater table occur in areas near mapped faults, suggesting that the groundwater flow in Burro Flats is strongly affected by the locations of fault traces. We also used the seismic data to develop seismic reflection images that show multiple strands of the SAF in the upper 60 m. Reflectors above the 10 m depth probably correspond to Holocene alluvial deposits; reflectors below the 15 m depth probably arise from velocity or density variations within the Precambrian gneiss complex, likely due to weathering. Apparent vertical offsets of reflectors are observed along profiles (lines 1 and 2) that are normal to the SAF, indicating minor apparent vertical offsets on the SAF at shallow depths. Along line 2, the apparently vertically offset reflectors correlate with zones of relatively low P -wave velocity. Along the central part of lines 1 and 2, the faults form a flower structure, which is typical of strike-slip faults such as the SAF.

The seismic quality factor Q for P waves in Atlantic Coastal Plain sediments is estimated using data from the 2004 U.S. Geological Survey seismic survey in eastern Virginia. The estimates are based on spectral ratios derived from reflections and sediment-guided P waves in Late Cretaceous and Tertiary sediments within the annular trough of the Late Eocene Chesapeake Bay impact structure. The estimates of Q for the frequency range of 10–150 Hz are from 75 to 100, with the best estimate of 80 based on multichannel stacking of spectral ratios from receivers in the offset range of 200–2000 m. This result is approximately a factor of 2 larger than the results previously reported for the Charleston, South Carolina, area, and it is approximately one-half of that recently reported for the Mississippi Embayment.

High-resolution seismic-reflection and seismic-refraction imaging, combined with existing borehole, earthquake, and paleoseismic trenching data, suggest that the Santa Monica fault zone in Los Angeles consists of multiple strands from several kilometers depth to the near surface. We interpret our seismic data as showing two shallow-depth low-angle fault strands and multiple near-vertical (∼85°) faults in the upper 100 m. One of the low-angle faults dips northward at about 28° and approaches the surface at the base of a topographic scarp on the grounds of the Wadsworth VA Hospital (WVAH). The other principal low-angle fault dips northward at about 20° and projects toward the surface about 200 m south of the topographic scarp, near the northernmost areas of the Los Angeles Basin that experienced strong shaking during the 1994 Northridge earthquake. The 20° north-dipping low-angle fault is also apparent on a previously published seismic-reflection image by Pratt et al. (1998) and appears to extend northward to at least Wilshire Boulevard, where the fault may be about 450 m below the surface. Slip rates determined at the WVAH site could be significantly underestimated if it is assumed that slip occurs only on a single strand of the Santa Monica fault or if it is assumed that the near-surface faults dip at angles greater than 20–28°. At the WVAH, tomographic velocity modeling shows a significant decrease in velocity across near-surface strands of the Santa Monica fault. P -wave velocities range from about 500 m/sec at the surface to about 4500 m/sec within the upper 50 m on the north side of the fault zone at WVAH, but maximum measured velocities on the south side of the low-angle fault zone at WVAH are about 3500 m/sec. These refraction velocities compare favorably with velocities measured in nearby boreholes by Gibbs et al. (2000). This study illustrates the utility of combined seismic-reflection and seismic-refraction methods, which allow more accurate reflection imaging and compositional estimations across areas with highly variable velocities, a property that is characteristic of most fault zones.

A 5-km-long, high-resolution seismic imaging survey across the San Andreas fault (SAF) zone and the proposed San Andreas Fault Observatory at Depth (SAFOD) drill site near Parkfield, California, shows that velocities vary both laterally and vertically. Velocities range from 4.0 km/sec) probably correspond to granitic rock of the Salinian block, which is exposed a few kilometers southwest of the SAF. The depth to the top of probable granitic rock varies laterally along the seismic profile but is about 600 m below the surface at the proposed SAFOD site. We observe a prominent, lateral low-velocity zone (LVZ) beneath and southwest of the surface trace of the SAF. The LVZ is about 1.5 km wide at 300-m depth but tapers to about 600 m wide at 750-m depth. At the maximum depth of the velocity model (750 m), the LVZ is centered approximately 400 m southwest of the surface trace of the SAF. Similar velocities and velocity gradients are observed at comparable depths on both sides of the LVZ, suggesting that the LVZ is anomalous relative to rocks on either side of it. Velocities within the LVZ are lower than those of San Andreas fault gouge, and the LVZ is also anomalous with respect to gravity, magnetic, and resistivity measurements. Because of its proximity to the surface trace of the SAF, it is tempting to suggest that the LVZ represents a zone of fractured crystalline rocks at depth. However, the LVZ instead probably represents a tectonic sliver of sedimentary rock that now rests adjacent to or encompasses the SAF. Such a sliver of sedimentary rock implies fault strands on both sides and possibly within the sliver, suggesting a zone of fault strands at least 1.5 km wide at a depth of 300 m, tapering to about 600 m wide at 750-m depth. Fluids within the sedimentary sliver are probably responsible for observed low-resistivity values.

Swanson et al. (2002), herein referred to as the “comment authors,” take issue with numerous aspects of our 1998 article (Catchings et al. , 1998) on Potrero Canyon, Los Angeles County, California. The comment authors take considerable liberty in paraphrasing our conclusions and discounting our interpretations, principally because they believe that our seismic images do not represent the geologic structure beneath Potrero Canyon. In particular, the comment authors contend that our subsurface images are inconsistent with surface geologic mapping. The comment authors attempt to support their objections by suggesting a combination of reflection artifacts and errors in our velocity model result in an “inaccurate seismic profile.” In response to their contentions, we show a published map, two published cross sections, and two sets of published borehole data that are consistent with our seismic image and the known geologic structure of the Potrero Canyon area. We present velocity measurements from a borehole along our seismic profile that are consistent with our velocity model and inconsistent with the comment authors’ assumed velocities. Furthermore, we show that the comment authors’ assumptions about our acquisition and processing techniques are based on inadequate research and lack of familiarity with seismological techniques. Below, we address 16 objections by the comment authors on a topical basis. Comment: “Catchings et al. (… the authors) stated that Potrero Canyon is believed to be within or overlie a broad syncline that is faulted along its axis. They interpreted this syncline to be present in their seismic profiles based on the shallow-dipping (0°–20°) reflectors shown in their figures 5 and 6, but they did not provide a geologic cross section at the seismic line to verify this. In contrast to the author’s statements and seismic profile, however, Winterer and Durham (1962) did not show a syncline below Potrero Canyon …

Models of P- and S-wave velocity, Vp/Vs ratios, Poisson's ratios, and density for the crust and upper mantle are presented along a 400-km-long profile trending from Memphis, Tennessee, to St. Louis, Missouri. The profile crosses the New Madrid seismic zone and reveals distinct regional variations in the crustal velocity structure north and south of the latitude of New Madrid. In the south near Memphis, the upper few kilometers of the crust are dominated by upper crustal sedimentary basins or graben with P-wave velocities less than 5 km/sec and S-wave velocities of about 2 km/sec. P-wave velocities of the upper and middle crust range from 6.0 to 6.5 km/sec at depths above 25 km, and corresponding S-wave velocities range from 3.5 to 3.7 km/sec. The lower crust consists of a high-velocity layer (Vp = 7.4 km/sec; Vs ∼4.2 km/sec) that is up to 20-km thick at the latitude of New Madrid but thins to about 15 km near Memphis. To the north, beneath the western-most Illinois basin, low-velocity (Vp < 5 km/sec; Vs < 2.3 km/sec) sedimentary basins are less than 1-km deep. The average velocities (Vp = 6.0 km/sec; Vs = 3.5 km/sec) of the underlying, near-surface rocks argue against large thicknesses of unconsolidated or noncarbonate sediments within 50 km of the western edge of the Illinois basin. Most of the crust beneath the Illinois basin is modeled as one layer, with velocities up to 6.8 km/sec (Vs = 3.7 km/sec) at 37-km depth. The thick, high-velocity (Vp = 7.4 km/sec; Vs ∼4.2 km/sec) lower crustal layer thins from about 20 km near New Madrid to about 6 km beneath the western Illinois basin. Refractions from the Moho and upper mantle occur as first arrivals over distances as a great as 160 km and reveal upper mantle layering to 60 km depth. Upper mantle layers with P-wave velocities of 8.2 km/sec (Vs = 4.5 km/sec) and 8.4 km/sec (Vs = 4.7 km/sec) are modeled at 43 and 60 km depth, respectively. Crustal Vp/Vs ratios range between 1.74 and 1.83, and upper mantle Vp/V s ratios range from 1.78 to 1.84. Poisson's ratios range from about 0.26 to 0.33 in the crust and from about 0.27 to 0.29 in the upper mantle. Modeled average densities range from about 2.55 in the sedimentary basins to 3.43 in the upper mantle. Geophysical characteristics of the crust and upper mantle within the New Madrid seismic zone are consistent with other continental rifts, but the crustal structure of the Illinois basin is not characteristics of most continental rift settings. Seismic and gravity data suggest a buried horst near the middle of Reelfoot rift, beneath which is a vertical zone of seismicity and velocity anomalies. The relative depth of the Reelfoot rift north and south of the Reelfoot graben suggests that the rift and its bounding faults may extend eastward beneath the city of Memphis.

Apparent southward-dipping, reverse-fault zones are imaged to depths of about 1.5 km beneath Potrero Canyon, Los Angeles County, California. Based on their orientation and projection to the surface, we suggest that the imaged fault zones are extensions of the Oak Ridge fault. Geologic mapping by others and correlations with seismicity studies suggest that the Oak Ridge fault is the causative fault of the 17 January 1994 Northridge earthquake (Northridge fault). Our seismically imaged faults may be among several faults that collectively comprise the Northridge thrust fault system. Unusually strong shaking in Potrero Canyon during the Northridge earthquake may have resulted from focusing of seismic energy or co-seismic movement along existing, related shallow-depth faults. The strong shaking produced ground-surface cracks and sand blows distributed along the length of the canyon. Seismic reflection and refraction images show that shallow-depth faults may underlie some of the observed surface cracks. The relationship between observed surface cracks and imaged faults indicates that some of the surface cracks may have developed from nontectonic alluvial movement, but others may be fault related. Immediately beneath the surface cracks, P-wave velocities are unusually low (

The 17 January 1994, Northridge, California, earthquake produced strong ground shaking at the Cedar Hills Nursery (referred to here as the Tarzana site) within the city of Tarzana, California, approximately 6 km from the epicenter of the mainshock. Although the Tarzana site is on a hill and is a rock site, accelerations of approximately 1.78 g horizontally and 1.2 g vertically at the Tarzana site are among the highest ever instrumentally recorded for an earthquake. To investigate possible site effects at the Tarzana site, we used explosive-source seismic refraction data to determine the shallow (

Our data indicate that critical and postcritical reflections from crustal layers and the Moho produced increased shaking at discrete distances along the San Francisco Peninsula during the 1989 Loma Prieta earthquake. These reflections may have produced an increase in amplitude that is as much as 10 times greater than that of the direct arrival. Peak amplitude-distance patterns measured from explosive sources, synthetic seismograms, aftershocks, and the mainshock of the 1989 Loma Prieta earthquake indicate that (1) point sources, such as explosions, produce similar peak amplitude-distance relations as distributed, double-couple sources when the crustal structure is approximately a Poisson solid; (2) peak amplitudes from smaller point sources (explosions) may be scaled to those of larger sources (earthquakes); and (3) reflections caused a pattern of high amplitudes at specific distances along the San Francisco Peninsula that geographically correlates with areas of reported damage following from the Loma Prieta mainshock. Our study indicates that critical and postcritical reflections were stronger influences on the locations of strong shaking than local geology because a number of sites around the San Francisco Bay that are underlain by unconsolidated sediments experienced much less shaking than sites that were underlain by hard rock. Furthermore, some severely shaken hard-rock sites were farther from the epicenter than many of the less severely shaken “soft-sediment” sites. Models of the regional crustal structure and seismic-wave propagation paths may be useful in identifying sites around the San Francisco Bay region that are at risk due to reflected seismic energy, and this type of analysis may be useful in other seismically active regions.

The New Madrid seismic zone (NMSZ) is the most active seismic zone in the United States east of the Rocky Mountains. Frequency-magnitude relationships show that an earthquake magnitude of 6.3 or greater can be expected in this area about once every 100 years, on average [Johnston and Nava, 1985]. During the winter of 1811–1812, three earthquakes of moment magnitude near 8 devastated New Madrid, Mo., and the surrounding region. Today, if an earthquake in the magnitude 6–7 range were to strike, in addition to loss of life and injuries, residents would suffer about $3.6 billion in immediate property losses [Algermissen, 1991] and feel the pinch of a long-term economic impact.

The San Andreas Fault Observatory at Depth pilot hole is located on the southwestern side of the Parkfield San Andreas fault. This observatory includes a vertical seismic profiling (VSP) array. VSP seismograms from nearby microearthquakes contain signals between thePandSwaves. These signals may bePandSwaves scattered by the local geologic structure. The collected scattering points form planar surfaces that we interpret as the San Andreas fault and four other secondary faults. The scattering process includes conversions betweenPandSwaves, the strengths of which suggest large contrasts in material properties, possibly indicating the presence of cracks or fluids.

The authors use a method of traveltime inversion of high-resolution seismic data to provide the first reliable images of internal details of the Columbia River Basalt Group (CRBG), the subsurface basalt/sediment interface, and the deeper sediment/basement interface. Velocity structure within the basalts, delineated on the order of 1 km horizontally and 0.2 km vertically, is constrained to within [plus minus]0.1 km/s for most of the seismic profile. Over 5,000 observed traveltimes fit their model with an rms error of 0.018 s. The maximum depth of penetration of the basalt diving waves (truncated by underlying low-velocity sediments) provides a reliable estimate of the depth to the base of the basalt, which agrees with well-log measurements to within 0.05 km (165 ft). The authors use image blurring, calculated from the resolution matrix, to estimate the aspect ratio of images velocity anomaly widths to true widths for velocity features within the basalt. From their calculations of image blurring, they interpret low velocity zones (LVZ) within the basalts at Boylston Mountain and the Whiskey Dick anticline to have widths of 4.5 and 3 km, respectively, within the upper 1.5 km of the model. At greater depth, the widths of these imaged LVZs thin tomore » approximately 2 km or less. They interpret these linear, subparallel, low-velocity zones imaged adjacent to anticlines of the Yakima Fold Belt to be brecciated fault zones. These fault zones dip to the south at angles between 15 to 45 degrees.« less

Clear subsurface seismic images have been obtained at low cost on the Columbia Plateau, Washington. The Columbia Plateau is perhaps the most notorious of all “bad‐data” areas because large impedance contrasts in surface flood basalts severely degrade the seismic wavefield. This degradation was mitigated in this study via a large‐explosive source, wide‐recording aperture shooting method. The shooting method emphasizes the wide‐angle portion of the wavefield, where Fermat’s principle guarantees reverberation will not interfere with the seismic manifestations of crucial geologic interfaces. The basalt diving wave, normally discarded in standard common midpoint (CMP) seismic profiling, can be used to image basalt velocity structure via traveltime inversion. Maximum depth‐penetration of the diving wave tightly constrains basalt‐sediment interface depth. An arrival observed only at shot‐receiver offsets greater than 15 km can be used to determine the velocity and geometry of basement via simultaneous inversion. The results from this study suggest that previous geologic hypotheses and hydrocarbon play concepts for the Columbia Plateau may have been in error.

Near-vertical and wide-angle seismic reflection data provide evidence for the presence of a magma body at the base of the crust beneath Buena Vista Valley in northwestern Nevada. The seismic response of this hypothesized magma body is characterized by high-amplitude, near-vertical P wave reflections and a comparably strong P-to-S converted phase. The magma body, referred to here as the Buena Vista Magma Body, is probably a single sill with thickness no greater than 200 m and length no greater than 1.8 km. The melt fraction in the sill is undoubtedly greater than 20–30%, and probably exceeds 50%. Melt composition is unconstrained. Although the age of the Buena Vista Magma Body is difficult to determine precisely, it is probably no older than 500,000 years. This suggests that magmatism in the Basin and Range Province is an ongoing process, despite the relative paucity of volcanic rocks erupted at the surface during the last 6 m.y.

Seismic reflection and refraction images illuminate the San Andreas Fault to a depth of 1 kilometer. The prestack depth-migrated reflection image contains near-vertical reflections aligned with the active fault trace. The fault is vertical in the upper 0.5 kilometer, then dips about 70° to the southwest to at least 1 kilometer subsurface. This dip reconciles the difference between the computed locations of earthquakes and the surface fault trace. The seismic velocity cross section shows strong lateral variations. Relatively low velocity (10 to 30%), high electrical conductivity, and low density indicate a 1-kilometer-wide vertical wedge of porous sediment or fractured rock immediately southwest of the active fault trace.

Structural details of the crust determined from P-wave velocity models can be improved with S-wave velocity models, and S-wave velocities are needed for model-based predictions of strong ground motion in southern California. We picked P- and S-wave travel times for refracted phases from explosive-source shots gathers of the Los Angeles Region Seismic Experiment, Phase II (LARSE II), and we developed refraction velocity models from these picks using two different inversion algorithms. Vp/Vs ratios were calculated from the resulting P- and S-wave models where both models are constrained by ray coverage. The two P-wave velocity models are compared to each other and to results from forward modeling. Generally, the P-wave inverse and forward models agree well for velocities lower than 5.0 km/s but only broadly agree with each other for velocities above 5.0 km/s. Similarly, the S-wave inverse models agree well with each other for velocities lower than 2.5 km/s but only broadly agree for velocities higher than 2.5 km/s. The most prominent structures in our S-wave models are two north-dipping low-velocity zones in the Central Transverse Ranges that we interpret as faults. These low-velocity zones differ somewhat between the two inversion models, but the Vp/Vs models (one model for each technique) show these features to be remarkably similar. Interestingly, both Vp/Vs models have several features that are not visible in either the P- or S-wave models alone. Two of these features (relatively high Vp/Vs ratios) occur in the vicinity of wells that bottom in "granitic" rocks and we interpret these high Vp/Vs ratios to indicate that the granitic rocks are highly fractured or even brecciated. Finally, to evaluate the Southern California Earthquake Center (SCEC) Community Velocity Model (CVM), which predicts Vs based on the Vp model, we compare data from our Vp and Vp/Vs models to empirical formulas that relate P- to S-wave velocities (see Brocher, 2005). These empirical curves provide an adequate average relationship between Vp and Vs, but our model Vs varies as much as ± 20 % for Vp ≤ 5.0 km/s and as much as ± 35 % for Vp ≥ 5.0 km/s. This large variation in the predicted S- wave velocity demonstrates the value of determining Vs independently from Vp.

In the northwestern to central Nevada Basin and Range, there are correlations between velocity and specific geologic structures of the crust. Mapped range-bounding faults at the surface can be traced to appreciable (10 km) depths based on velocity variations and are consistent with subsurface projections of the faults based on seismic reflection images. The limiting depth of the faults, as indicated by the velocity variations, corresponds to the maximum depth of earthquakes along the seismic profile. Correlations between velocity and the surface geology show that in the upper crust the pre-Cenozoic rocks are underlain by high-velocity (6.0 km/s) rocks, whereas the Tertiary ranges are underlain by lower-velocity (4.0-5.7 km/s) rocks to depths as great as 10 km. Although the Tertiary rocks differ in composition from the Mesozoic rocks, the lower-velocity Tertiary rocks may also be attributed to rock masses which are broken (4.0-5.7 km/s), and the higher-velocity Mesozoic rocks (6.0 km/s) may be attributed to largely unbroken rock masses. The regional seismicity pattern is consistent with this interpretation, as earthquakes are largely confined within or near the base of the low-velocity rocks. These low-velocity, highly fractured rocks are laterally distributed in discrete zones, suggesting that extension is not uniformly distributed but occurs in discrete, highly extended zones. Beneath these highly extended zones, the lower-crustal layers show structural evidence of extension, and velocity measurements suggest that the lowermost crust has been magmatically underplated. The superposition of Tertiary volcanic rocks, highly fractured upper crust, and lower-crustal magmatic underplating suggests that the Tertiary volcanic rocks originated from lower-crustal magmas that migrated to the surface via the highly extended zones. The velocity structure of one of the highly extended zones and the Lahonton Basin resembles that of many continental rifts. The velocity structure beneath central Nevada, however, is much more like normal continental crust. On the basis of isotopic studies, it is concluded that the transition between highly extended crust and more normal crust occurs in the area inferred to be the edge of the North American craton.

We present an interpretation of the crustal and uppermost mantle structure of the Basin and Range of northwestern Nevada based on seismic refraction/wide-angle reflection, near-vertical reflection, and gravity data. In comparison to most previous estimates, we find that the crust is somewhat thicker (32–36 km versus 22–30 km), and the uppermost mantle velocity is somewhat higher (8.0 km/s versus 7.3–7.9 km/s). Along our transects, the crust is thinnest (32 km) in the Carson Sink-Buena Vista Valley region and increases by 2–4 km to the west and east, respectively. There is considerable complexity throughout the crust where velocities range from of 2.5 km/s at the surface to 7.4 km/s in the lowermost crust. Variations in velocity and structure of the upper crustal layers reveal apparent basement velocity depressions (areas of lower velocities extending up to 10 km in depth) that underlie some surface ranges as well as the basins. The middle crust rises from about 20 km beneath central Nevada to within 12 km of the surface beneath the area of thinnest crust and is characterized by a modest (∼0.1 km/s) change in velocity and low-velocity gradients. These midcrustal layers mark the onset of high crustal reflectivity and the apparent limiting depth to which Basin and Range faults can be traced in near-vertical reflection profiles, suggesting that these midcrustal layers represent the transition between the brittle and ductile zones of the crust. The lower crust is more structurally complex, with layers thickening and thinning in a systematic manner with the upper crustal layers; generally, where there are velocity depressions in the upper crust, the lower crust is thickest and shallowest. The geometry of these lower crustal layers (derived from refraction modeling) coincides with changes in the crustal reflectivity, determined from the Consortium of Continental Reflection Profiling reflection data. The lower crustal layer is unusually high in velocity (7.4 km/s) and is likely the layer identified as mantle in some previous studies. We do not identify the 7.4 km/s layer as mantle because (1) there is an underlying layer with a velocity (8.0 km/s) that is more consistent with the worldwide average velocity for the upper mantle, and (2) the 7.4 km/s layer does not correspond to the “reflection” Moho. Gravity modeling and comparison to existing seismic models show a general consensus in many aspects with respect to crustal structure. This new model forms the basis for speculation on some of the processes associated with rifting of the Basin and Range Province. One such process, lithospheric magmatism, is inferred from the strong attenuation of transmitted seismic waves, which occurs at the same interface at which high-amplitude, bright spot reflections originate. Unlike previous models, the overall structure and velocity of the crust and uppermost mantle of our new model are similar to other regions worldwide which have undergone high degrees of extension.

[1] A 30-km-long, radial seismic reflection and refraction survey completed across the northern part of the late Eocene Chesapeake Bay impact structure (CBIS) on the Delmarva Peninsula, Virginia, USA, confirms that the CBIS is a complex central-peak crater. We used a tomographic P wave velocity model and low-fold reflection images, constrained by data from two deep boreholes located on the profile, to interpret the structure and composition of the upper 5 km of crust. The seismic images exhibit well-defined structural features, including (with increasing radial distance) a collapsed central uplift, a breccia-filled moat, and a collapsed transient-crater margin (which collectively constitute a ∼40-km-wide collapsed transient crater), and a shallowly deformed annular trough. These seismic images are the first to resolve the deep structure of the crater (>1 km) and the boundaries between the central uplift, moat, and annular trough. Several distinct seismic signatures distinguish breccia units from each other and from more coherent crystalline rocks below the central uplift, moat, and annular trough. Within the moat, breccia extends to a minimum depth of 1.5 km or a maximum of 3.5 km, depending upon the interpretation of the deepest layered materials. The images show ∼350 to 500 m of postimpact sediments above the impactites. The imaged structure of the CBIS indicates a complex sequence of events during the cratering process that will provide new constraints for numerical modeling.