Search USGS⌕ Search

Geology topics

C. Thurber

Publications and source records attributed to C. Thurber.

10 recordsLinked to original sources
USGS2011

Imaging the San Andreas fault with explosion and earthquake sources

Mounting evidence suggests that fault zone heterogeneity may play a crucial role in the localization of rupture in earthquakes [ Aki , 1995]. The heterogeneity can take several forms: spatial variations in physical properties (elastic properties, pore fluid pressure, etc.) or complexity in the fault surface (bends, offsets, etc.). High-resolution, three-dimensional models of the P and S wave velocity ( V p and V s ) structure and accurate hypocenters provide tools for studying the effects of fault zone heterogeneity on rupture localization. A dense, passive array of 48 seismic instruments was deployed south of Hollister, California, from mid-November 1994 to late May 1995. The array obtained local-earthquake seismograms for high-resolution, three-dimensional imaging of the ( V p and V s ) structure in the San Andreas fault zone. In mid-May, anactive seismic experiment was also carried out. Figure 1 shows a 7-km wide, 3- to 4-km deep zone of very low V p imaged during the experiment. The zone is bounded on the southwest by the San Andreas fault and the adjacent Northern Gabilan Range. The San Andreas fault appears to be approximately vertical to at least 5 km depth (the limit of good model resolution), contrary to an existing model based on gravity data.

California↗
USGS2009

Joint inversion for Vp, Vs, and Vp/Vs at SAFOD, Parkfield, California

We refined the three-dimensional (3-D) Vp, Vs and Vp/Vs models around the San Andreas Fault Observatory at Depth (SAFOD) site using a new double-difference (DD) seismic tomography code (tomoDDPS) that simultaneously solves for earthquake locations and all three velocity models using both absolute and differential P, S, and S-P times. This new method is able to provide a more robust Vp/Vs model than that from the original DD tomography code (tomoDD), obtained simply by dividing Vp by Vs. For the new inversion, waveform cross-correlation times for earthquakes from 2001 to 2002 were also used, in addition to arrival times from earthquakes and explosions in the region. The Vp values extracted from the model along the SAFOD trajectory match well with the borehole log data, providing in situ confirmation of our results. Similar to previous tomographic studies, the 3-D structure around Parkfield is dominated by the velocity contrast across the San Andreas Fault (SAF). In both the Vp and Vs models, there is a clear low-velocity zone as deep as 7 km along the SAF trace, compatible with the findings from fault zone guided waves. There is a high Vp/Vs anomaly zone on the southwest side of the SAF trace that is about 1-2 km wide and extends as deep as 4 km, which is interpreted to be due to fluids and fractures in the package of sedimentary rocks abutting the Salinian basement rock to the southwest. The relocated earthquakes align beneath the northeast edge of this high Vp/Vs zone. We carried out a 2-D correlation analysis for an existing resistivity model and the corresponding profiles through our model, yielding a classification that distinguishes several major lithologies. ?? 2009 by the American Geophysical Union.

Geochemistry, Geophysics, Geosystems↗
USGS2009

Regional three-dimensional seismic velocity model of the crust and uppermost mantle of northern California

We present a three-dimensional (3D) tomographic model of the P wave velocity (Vp) structure of northern California. We employed a regional-scale double-difference tomography algorithm that incorporates a finite-difference travel time calculator and spatial smoothing constraints. Arrival times from earthquakes and travel times from controlled-source explosions, recorded at network and/or temporary stations, were inverted for Vp on a 3D grid with horizontal node spacing of 10 to 20 km and vertical node spacing of 3 to 8 km. Our model provides an unprecedented, comprehensive view of the regional-scale structure of northern California, putting many previously identified features into a broader regional context and improving the resolution of a number of them and revealing a number of new features, especially in the middle and lower crust, that have never before been reported. Examples of the former include the complex subducting Gorda slab, a steep, deeply penetrating fault beneath the Sacramento River Delta, crustal low-velocity zones beneath Geysers-Clear Lake and Long Valley, and the high-velocity ophiolite body underlying the Great Valley. Examples of the latter include mid-crustal low-velocity zones beneath Mount Shasta and north of Lake Tahoe. Copyright 2009 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2009

High precision relocation of earthquakes at Iliamna Volcano, Alaska

In August 1996, a period of elevated seismicity commenced beneath Iliamna Volcano, Alaska. This activity lasted until early 1997, consisted of over 3000 earthquakes, and was accompanied by elevated emissions of volcanic gases. No eruption occurred and seismicity returned to background levels where it has remained since. We use waveform alignment with bispectrum-verified cross-correlation and double-difference methods to relocate over 2000 earthquakes from 1996 to 2005 with high precision (~ 100 m). The results of this analysis greatly clarify the distribution of seismic activity, revealing distinct features previously hidden by location scatter. A set of linear earthquake clusters diverges upward and southward from the main group of earthquakes. The events in these linear clusters show a clear southward migration with time. We suggest that these earthquakes represent either a response to degassing of the magma body, circulation of fluids due to exsolution from magma or heating of ground water, or possibly the intrusion of new dikes beneath Iliamna's southern flank. In addition, we speculate that the deeper, somewhat diffuse cluster of seismicity near and south of Iliamna's summit indicates the presence of an underlying magma body between about 2 and 4 km depth below sea level, based on similar features found previously at several other Alaskan volcanoes.

Alaska↗
USGS2006

Three-dimensional compressional wavespeed model, earthquake relocations, and focal mechanisms for the Parkfield, California, region

We present a new three-dimensional (3D) compressional vvavespeed (V p) model for the Parkfield region, taking advantage of the recent seismicity associated with the 2003 San Simeon and 2004 Parkfield earthquake sequences to provide increased model resolution compared to the work of Eberhart-Phillips and Michael (1993) (EPM93). Taking the EPM93 3D model as our starting model, we invert the arrival-time data from about 2100 earthquakes and 250 shots recorded on both permanent network and temporary stations in a region 130 km northeast-southwest by 120 km northwest-southeast. We include catalog picks and cross-correlation and catalog differential times in the inversion, using the double-difference tomography method of Zhang and Thurber (2003). The principal Vp features reported by EPM93 and Michelini and McEvilly (1991) are recovered, but with locally improved resolution along the San Andreas Fault (SAF) and near the active-source profiles. We image the previously identified strong wavespeed contrast (faster on the southwest side) across most of the length of the SAF, and we also improve the image of a high Vp body on the northeast side of the fault reported by EPM93. This narrow body is at about 5- to 12-km depth and extends approximately from the locked section of the SAP to the town of Parkfield. The footwall of the thrust fault responsible for the 1983 Coalinga earthquake is imaged as a northeast-dipping high wavespeed body. In between, relatively low wavespeeds (<5 km/sec) extend to as much as 10-km depth. We use this model to derive absolute locations for about 16,000 earthquakes from 1966 to 2005 and high-precision double-difference locations for 9,000 earthquakes from 1984 to 2005, and also to determine focal mechanisms for 446 earthquakes. These earthquake locations and mechanisms show that the seismogenic fault is a simple planar structure. The aftershock sequence of the 2004 mainshock concentrates into the same structures defined by the pre-2004 seismicity, confirming earlier observations (Waldhauser et al., 2004) that the seismicity pattern at Parkfield is long lived and persists through multiple cycles of mainshocks.

Bulletin of the Seismological Society of America↗
USGS2004

Joint inversion of gravity and arrival time data from Parkfield: New constraints on structure and hypocenter locations near the SAFOD drill site

Taking advantage of large datasets of both gravity and elastic wave arrival time observations available for the Parkfield, California region, we generated an image consistent with both types of data. Among a variety of strategies, the best result was obtained from a simultaneous inversion with a stability requirement that encouraged the perturbed model to remain close to a starting model consisting of a best fit to the arrival time data. The preferred model looks essentially the same as the best-fit arrival time model in areas where ray coverage is dense, with differences being greatest at shallow depths and near the edges of the model where ray paths are few. Earthquake locations change by no more than about 100 m, the general effect being migration of the seismic zone to the northeast, closer to the surface trace of the San Andreas Fault. Copyright 2004 by the American Geophysical Union.

Geophysical Research Letters↗
USGS2004

Fine-scale structure of the San Andreas fault zone and location of the SAFOD target earthquakes

We present results from the tomographic analysis of seismic data from the Parkfield area using three different inversion codes. The models provide a consistent view of the complex velocity structure in the vicinity of the San Andreas, including a sharp velocity contrast across the fault. We use the inversion results to assess our confidence in the absolute location accuracy of a potential target earthquake. We derive two types of accuracy estimates, one based on a consideration of the location differences from the three inversion methods, and the other based on the absolute location accuracy of “virtual earthquakes.” Location differences are on the order of 100–200 m horizontally and up to 500 m vertically. Bounds on the absolute location errors based on the “virtual earthquake” relocations are ∼50 m horizontally and vertically. The average of our locations places the target event epicenter within about 100 m of the SAF surface trace.

California↗
USGS2002

Ground truth seismic events and location capability at Degelen mountain, Kazakhstan

We utilized nuclear explosions from the Degelen Mountain sub-region of the Semipalatinsk Test Site (STS), Kazakhstan, to assess seismic location capability directly. Excellent ground truth information for these events was either known or was estimated from maps of the Degelen Mountain adit complex. Origin times were refined for events for which absolute origin time information was unknown using catalog arrival times, our ground truth location estimates, and a time baseline provided by fixing known origin times during a joint hypocenter determination (JHD). Precise arrival time picks were determined using a waveform cross-correlation process applied to the available digital data. These data were used in a JHD analysis. We found that very accurate locations were possible when high precision, waveform cross-correlation arrival times were combined with JHD. Relocation with our full digital data set resulted in a mean mislocation of 2 km and a mean 95% confidence ellipse (CE) area of 6.6 km2 (90% CE: 5.1 km2), however, only 5 of the 18 computed error ellipses actually covered the associated ground truth location estimate. To test a more realistic nuclear test monitoring scenario, we applied our JHD analysis to a set of seven events (one fixed) using data only from seismic stations within 40?? epicentral distance. Relocation with these data resulted in a mean mislocation of 7.4 km, with four of the 95% error ellipses covering less than 570 km2 (90% CE: 438 km2), and the other two covering 1730 and 8869 km2 (90% CE: 1331 and 6822 km2). Location uncertainties calculated using JHD often underestimated the true error, but a circular region with a radius equal to the mislocation covered less than 1000 km2 for all events having more than three observations. ?? 2002 Elsevier Science B.V. All rights reserved.

Physics of the Earth and Planetary Interiors↗
USGS1997

Two-dimensional seismic image of the San Andreas Fault in the Northern Gabilan Range, central California: Evidence for fluids in the fault zone

A joint inversion for two-dimensional P-wave velocity (Vp), P-to-S velocity ratio (Vp/Vs), and earthquake locations along the San Andreas fault (SAF) in central California reveals a complex relationship among seismicity, fault zone structure, and the surface fault trace. A zone of low Vp and high Vp/Vs lies beneath the SAF surface trace (SAFST), extending to a depth of about 6 km. Most of the seismic activity along the SAF occurs at depths of 3 to 7 km in a southwest-dipping zone that roughly intersects the SAFST, and lies near the southwest edge of the low Vp and high Vp/Vs zones. Tests indicate that models in which this seismic zone is significantly closer to vertical can be confidently rejected. A second high Vp/Vs zone extends to the northeast, apparently dipping beneath the Diablo Range. Another zone of seismicity underlies the northeast portion of this Vp/Vs high. The high Vp/Vs zones cut across areas of very different Vp values, indicating that the high Vp/Vs values are due to the presence of fluids, not just lithology. The close association between the zones of high Vp/Vs and seismicity suggests a direct involvement of fluids in the faulting process.

Geophysical Research Letters↗