USGS Science⌕ Search

Geology topics

F. F. Pollitz

Publications and source records attributed to F. F. Pollitz.

At least 19 recordsLinked to original sources
USGS2011

Epistemic uncertainty in California-wide synthetic seismicity simulations

The generation of seismicity catalogs on synthetic fault networks holds the promise of providing key inputs into probabilistic seismic-hazard analysis, for example, the coefficient of variation, mean recurrence time as a function of magnitude, the probability of fault-to-fault ruptures, and conditional probabilities for foreshock–mainshock triggering. I employ a seismicity simulator that includes the following ingredients: static stress transfer, viscoelastic relaxation of the lower crust and mantle, and vertical stratification of elastic and viscoelastic material properties. A cascade mechanism combined with a simple Coulomb failure criterion is used to determine the initiation, propagation, and termination of synthetic ruptures. It is employed on a 3D fault network provided by Steve Ward (unpublished data, 2009) for the Southern California Earthquake Center (SCEC) Earthquake Simulators Group. This all-California fault network, initially consisting of 8000 patches, each of ∼12 square kilometers in size, has been rediscretized into Graphic patches, each of ∼1 square kilometer in size, in order to simulate the evolution of California seismicity and crustal stress at magnitude M∼5–8. Resulting synthetic seismicity catalogs spanning 30,000 yr and about one-half million events are evaluated with magnitude-frequency and magnitude-area statistics. For a priori choices of fault-slip rates and mean stress drops, I explore the sensitivity of various constructs on input parameters, particularly mantle viscosity. Slip maps obtained for the southern San Andreas fault show that the ability of segment boundaries to inhibit slip across the boundaries (e.g., to prevent multisegment ruptures) is systematically affected by mantle viscosity.

California↗
USGS2010

Rayleigh-wave phase-velocity maps and three-dimensional shear velocity structure of the western US from local non-plane surface wave tomography

We utilize two-and-three-quarter years of vertical-component recordings made by the Transportable Array (TA) component of Earthscope to constrain three-dimensional (3-D) seismic shear wave velocity structure in the upper 200 km of the western United States. Single-taper spectral estimation is used to compile measurements of complex spectral amplitudes from 44 317 seismograms generated by 123 teleseismic events. In the first step employed to determine the Rayleigh-wave phase-velocity structure, we implement a new tomographic method, which is simpler and more robust than scattering-based methods (e.g. multi-plane surface wave tomography). The TA is effectively implemented as a large number of local arrays by defining a horizontal Gaussian smoothing distance that weights observations near a given target point. The complex spectral-amplitude measurements are interpreted with the spherical Helmholtz equation using local observations about a succession of target points, resulting in Rayleigh-wave phase-velocity maps at periods over the range of 18–125 s. The derived maps depend on the form of local fits to the Helmholtz equation, which generally involve the nonplane-wave solutions of Friederich et al. In a second step, the phase-velocity maps are used to derive 3-D shear velocity structure. The 3-D velocity images confirm details witnessed in prior body-wave and surface-wave studies and reveal new structures, including a deep (>100 km deep) high-velocity lineament, of width ∼200 km, stretching from the southern Great Valley to northern Utah that may be a relic of plate subduction or, alternatively, either a remnant of the Mojave Precambrian Province or a mantle downwelling. Mantle seismic velocity is highly correlated with heat flow, Holocene volcanism, elastic plate thickness and seismicity. This suggests that shallow mantle structure provides the heat source for associated magmatism, as well as thinning of the thermal lithosphere, leading to relatively high stress concentration. Our images also confirm the presence of high-velocity mantle at 100 km depth beneath areas of suspected mantle delamination (southern Sierra Nevada; Grande Ronde uplift), low velocity mantle underlying active rift zones, and high velocity mantle associated with the subducting Juan de Fuca plate. Structure established during the Proterozoic appears to exert a lasting influence on subsequent volcanism and tectonism up to the Present.

Geophysical Journal International↗
USGS2010

Viscoelastic-cycle model of interseismic deformation in the northwestern United States

We apply a viscoelastic cycle model to a compilation of GPS velocity fields in order to address the kinematics of deformation in the northwestern United States. A viscoelastic cycle model accounts for time-dependent deformation following large crustal earthquakes and is an alternative to block models for explaining the interseismic crustal velocity field. Building on the approach taken in Pollitz et al., we construct a deformation model for the entire western United States-based on combined fault slip and distributed deformation-and focus on the implications for the Mendocino triple junction (MTJ), Cascadia megathrust, and western Washington. We find significant partitioning between strike-slip and dip-slip motion near the MTJ as the tectonic environment shifts from northwest-directed shear along the San Andreas fault system to east-west convergence along the Juan de Fuca Plate. By better accounting for the budget of aseismic and seismic slip along the Cascadia subduction interface in conjunction with an assumed rheology, we revise a previous model of slip for the M~ 9 1700 Cascadia earthquake. In western Washington, we infer slip rates on a number of strike-slip and dip-slip faults that accommodate northward convergence of the Oregon Coast block and northwestward convergence of the Juan de Fuca Plate. Lateral variations in first order mechanical properties (e.g. mantle viscosity, vertically averaged rigidity) explain, to a large extent, crustal strain that cannot be rationalized with cyclic deformation on a laterally homogeneous viscoelastic structure. Our analysis also shows that present crustal deformation measurements, particularly with the addition of the Plate Boundary Observatory, can constrain such lateral variations.

Geophysical Journal International↗
USGS2009

The 1911 M ~6.6 Calaveras earthquake: Source parameters and the role of static, viscoelastic, and dynamic coulomb stress changes imparted by the 1906 San Francisco earthquake

The occurrence of a right-lateral strike-slip earthquake in 1911 is inconsistent with the calculated 0.2-2.5 bar static stress decrease imparted by the 1906 rupture at that location on the Calaveras fault, and 5 yr of calculated post-1906 viscoelastic rebound does little to reload the fault. We have used all available first-motion, body-wave, and surface-wave data to explore possible focal mechanisms for the 1911 earthquake. We find that the event was most likely a right-lateral strikeslip event on the Calaveras fault, larger than, but otherwise resembling, the 1984 M w 6.1 Morgan Hill earthquake in roughly the same location. Unfortunately, we could recover no unambiguous surface fault offset or geodetic strain data to corroborate the seismic analysis despite an exhaustive archival search. We calculated the static and dynamic Coulomb stress changes for three 1906 source models to understand stress transfer to the 1911 site. In contrast to the static stress shadow, the peak dynamic Coulomb stress imparted by the 1906 rupture promoted failure at the site of the 1911 earthquake by 1.4-5.8 bar. Perhaps because the sample is small and the aftershocks are poorly located, we find no correlation of 1906 aftershock frequency or magnitude with the peak dynamic stress, although all aftershocks sustained a calculated dynamic stress of ???3 bar. Just 20 km to the south of the 1911 epicenter, we find that surface creep of the Calaveras fault at Hollister paused for ~17 yr after 1906, about the expected delay for the calculated static stress drop imparted by the 1906 earthquake when San Andreas fault postseismic creep and viscoelastic relaxation are included. Thus, the 1911 earthquake may have been promoted by the transient dynamic stresses, while Calaveras fault creep 20 km to the south appears to have been inhibited by the static stress changes.

Bulletin of the Seismological Society of America↗
USGS2008

Observations and interpretation of fundamental mode Rayleigh wavefields recorded by the Transportable Array (USArray)

Broadband recordings of the dense Transportable Array (TA) in the western United States provide unparalleled detailed images of long-period seismic surface wavefields. With 400 stations spanning most of the western United States, wavefronts of fundamental mode Rayleigh waves may be visualized coherently across the array at periods ???40 s. In order to constrain the Rayleigh wave phase velocity structure in the western United States, I assemble a data set of vertical component seismograms from 53 teleseismic events recorded by the TA from April 2006 to October 2007. Complex amplitude spectra from these recordings at peni ods 27-100 s are interpreted using the multiplane wave tomographic method of Friederich and Wielandt (1995) and Pollitz (1999). This analysis yields detailed surface wave phase velocity and three-dimensional shear wave velocity patterns across the North American plate boundary zone, elucidating the active processes in the highly heterogeneous western U.S. upper mantle.

Journal of Geophysical Research B: Solid Earth↗
USGS2008

Probabilistic seismic hazard in the San Francisco Bay area based on a simplified viscoelastic cycle model of fault interactions

We construct a viscoelastic cycle model of plate boundary deformation that includes the effect of time-dependent interseismic strain accumulation, coseismic strain release, and viscoelastic relaxation of the substrate beneath the seismogenic crust. For a given fault system, time-averaged stress changes at any point (not on a fault) are constrained to zero; that is, kinematic consistency is enforced for the fault system. The dates of last rupture, mean recurrence times, and the slip distributions of the (assumed) repeating ruptures are key inputs into the viscoelastic cycle model. This simple formulation allows construction of stress evolution at all points in the plate boundary zone for purposes of probabilistic seismic hazard analysis (PSHA). Stress evolution is combined with a Coulomb failure stress threshold at representative points on the fault segments to estimate the times of their respective future ruptures. In our PSHA we consider uncertainties in a four-dimensional parameter space: the rupture peridocities, slip distributions, time of last earthquake (for prehistoric ruptures) and Coulomb failure stress thresholds. We apply this methodology to the San Francisco Bay region using a recently determined fault chronology of area faults. Assuming single-segment rupture scenarios, we find that fature rupture probabilities of area faults in the coming decades are the highest for the southern Hayward, Rodgers Creek, and northern Calaveras faults. This conclusion is qualitatively similar to that of Working Group on California Earthquake Probabilities, but the probabilities derived here are significantly higher. Given that fault rupture probabilities are highly model-dependent, no single model should be used to assess to time-dependent rupture probabilities. We suggest that several models, including the present one, be used in a comprehensive PSHA methodology, as was done by Working Group on California Earthquake Probabilities.

Journal of Geophysical Research B: Solid Earth↗
USGS2008

Dislocation models of interseismic deformation in the western United States

The GPS-derived crustal velocity field of the western United States is used to construct dislocation models in a viscoelastic medium of interseismic crustal deformation. The interseismic velocity field is constrained by 1052 GPS velocity vectors spanning the ???2500-km-long plate boundary zone adjacent to the San Andreas fault and Cascadia subduction zone and extending ???1000 km into the plate interior. The GPS data set is compiled from U.S. Geological Survey campaign data, Plate Boundary Observatory data, and the Western U.S. Cordillera velocity field of Bennett et al. (1999). In the context of viscoelastic cycle models of postearthquake deformation, the interseismic velocity field is modeled with a combination of earthquake sources on ???100 known faults plus broadly distributed sources. Models that best explain the observed interseismic velocity field include the contributions of viscoelastic relaxation from faulting near the major plate margins, viscoelastic relaxation from distributed faulting in the plate interior, as well as lateral variations in depth-averaged rigidity in the elastic lithosphere. Resulting rigidity variations are consistent with reduced effective elastic plate thickness in a zone a few tens of kilometers wide surrounding the San Andreas fault (SAF) system. Primary deformation characteristics are captured along the entire SAF system, Eastern California Shear Zone, Walker Lane, the Mendocino triple junction, the Cascadia margin, and the plate interior up to ???1000 km from the major plate boundaries.

Journal of Geophysical Research B: Solid Earth↗
USGS2008

Temporal evolution of continental lithospheric strength in actively deforming regions

It has been agreed for nearly a century that a strong, load-bearing outer layer of earth is required to support mountain ranges, transmit stresses to deform active regions and store elastic strain to generate earthquakes. However the dept and extent of this strong layer remain controversial. Here we use a variety of observations to infer the distribution of lithospheric strength in the active western United States from seismic to steady-state time scales. We use evidence from post-seismic transient and earthquake cycle deformation reservoir loading glacio-isostatic adjustment, and lithosphere isostatic adjustment to large surface and subsurface loads. The nearly perfectly elastic behavior of Earth's crust and mantle at the time scale of seismic wave propagation evolves to that of a strong, elastic crust and weak, ductile upper mantle lithosphere at both earthquake cycle (EC, ???10?? to 103 yr) and glacio-isostatic adjustment (GIA, ???103 to 104 yr) time scales. Topography and gravity field correlations indicate that lithosphere isostatic adjustment (LIA) on ???106-107 yr time scales occurs with most lithospheric stress supported by an upper crust overlying a much weaker ductile subtrate. These comparisons suggest that the upper mantle lithosphere is weaker than the crust at all time scales longer than seismic. In contrast, the lower crust has a chameleon-like behavior, strong at EC and GIA time scales and weak for LIA and steady-state deformation processes. The lower crust might even take on a third identity in regions of rapid crustal extension or continental collision, where anomalously high temperatures may lead to large-scale ductile flow in a lower crustal layer that is locally weaker than the upper mantle. Modeling of lithospheric processes in active regions thus cannot use a one-size-fits-all prescription of rheological layering (relation between applied stress and deformation as a function of depth) but must be tailored to the time scale and tectonic setting of the process being investigated.

GSA Today↗
USGS2007

Finite-frequency traveltime tomography of San Francisco Bay region crustal velocity structure

Seismic velocity structure of the San Francisco Bay region crust is derived using measurements of finite-frequency traveltimes. A total of 57 801 relative traveltimes are measured by cross-correlation over the frequency range 0.5–1.5 Hz. From these are derived 4862 ‘summary’ traveltimes, which are used to derive 3-D P -wave velocity structure over a 341 × 140 km 2 area from the surface to 25 km depth. The seismic tomography is based on sensitivity kernels calculated on a spherically symmetric reference model. Robust elements of the derived P -wave velocity structure are: a pronounced velocity contrast across the San Andreas fault in the south Bay region (west side faster); a moderate velocity contrast across the Hayward fault (west side faster); moderately low velocity crust around the Quien Sabe volcanic field and the Sacramento River delta; very low velocity crust around Lake Berryessa. These features are generally explicable with surface rock types being extrapolated to depth ∼10 km in the upper crust. Generally high mid-lower crust velocity and high inferred Poisson's ratio suggest a mafic lower crust.

California↗
USGS2006

The 1923 Kanto earthquake reevaluated using a newly augmented geodetic data set

This study revisits the mechanism of the 1923 Ms = 7.9 Kanto earthquake in Japan. We derive a new source model and use it to assess quantitative and qualitative aspects of the accommodation of plate motion in the Kanto region. We use a new geodetic data set that consists of displacements from leveling and angle changes from triangulation measurements obtained in surveys between 1883 and 1927. Two unique aspects of our analysis are the inclusion of a large number of second-order triangulation measurements and the application of a correction to remove interseismic deformation. The geometry of the fault planes is adopted from a recent seismic reflection study of the Kanto region. We evaluate the minimum complexity necessary in the model to fit the data optimally. Our final uniform-slip elastic dislocation model consists of two adjacent ???20?? dipping low-angle planes accommodating reverse dextral slip of 6.0 in on the larger, eastern plane and 9.5 m on the smaller, western plane with azimuths of 163?? and 121??, respectively. The earthquake was located in the Sagami trough, where the Philippine Sea plate subducts under Honshu. Compared to the highly oblique angle of plate convergence, the coseismic slip on the large fault plane has a more orthogonal orientation to the strike of the plate boundary, suggesting that slip partitioning plays a role in accommodation of plate motion. What other structure is involved in the partitioning is unclear. Uplift records of marine coastal terraces in Sagami Bay document 7500 years of earthquake activity and predict average recurrence intervals of 400 years for events with vertical displacement profiles similar to those of the 1923 earthquake. This means that the average slip deficit per recurrence interval is ???50% of the relative plate convergence. These findings of plate motion partitioning and slip deficit lead us to suggest that instead of a simple recurrence model with characteristic earthquakes, additional mechanisms are necessary to describe the accommodation of deformation in the Kanto region. So far, obvious candidates for these alternative mechanisms have not been discovered. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2006

Mechanical deformation model of the western United States instantaneous strain-rate field

We present a relationship between the long-term fault slip rates and instantaneous velocities as measured by Global Positioning System (GPS) or other geodetic measurements over a short time span. The main elements are the secularly increasing forces imposed by the bounding Pacific and Juan de Fuca (JdF) plates on the North American plate, viscoelastic relaxation following selected large earthquakes occurring on faults that are locked during their respective interseismic periods, and steady slip along creeping portions of faults in the context of a thin-plate system. In detail, the physical model allows separate treatments of faults with known geometry and slip history, faults with incomplete characterization (i.e. fault geometry but not necessarily slip history is available), creeping faults, and dislocation sources distributed between the faults. We model the western United States strain-rate field, derived from 746 GPS velocity vectors, in order to test the importance of the relaxation from historic events and characterize the tectonic forces imposed by the bounding Pacific and JdF plates. Relaxation following major earthquakes (M ??? 8.0) strongly shapes the present strain-rate field over most of the plate boundary zone. Equally important are lateral shear transmitted across the Pacific-North America plate boundary along ???1000 km of the continental shelf, downdip forces distributed along the Cascadia subduction interface, and distributed slip in the lower lithosphere. Post-earthquake relaxation and tectonic forcing, combined with distributed deep slip, constructively interfere near the western margin of the plate boundary zone, producing locally large strain accumulation along the San Andreas fault (SAF) system. However, they destructively interfere further into the plate interior, resulting in smaller and more variable strain accumulation patterns in the eastern part of the plate boundary zone. Much of the right-lateral strain accumulation along the SAF system is systematically underpredicted by models which account only for relaxation from known large earthquakes. This strongly suggests that in addition to viscoelastic-cycle effects, steady deep slip in the lower lithosphere is needed to explain the observed strain-rate field. ?? 2006 The Authors Journal compilation ?? 2006 RAS.

Geophysical Journal International↗
USGS2006

Direct test of static stress versus dynamic stress triggering of aftershocks

Aftershocks observed over time scales of minutes to months following a main shock are plausibly triggered by the static stress change imparted by the main shock, dynamic shaking effects associated with passage of seismic waves from the main shock, or a combination of the two. We design a direct test of static versus dynamic triggering of aftershocks by comparing the near-field temporal aftershock patterns generated by aseismic and impulsive events occurring in the same source area. The San Juan Bautista, California, area is ideally suited for this purpose because several events of both types of M???5 have occurred since 1974. We find that aftershock rates observed after impulsive events are much higher than those observed after aseismic events, and this pattern persists for several weeks after the event. This suggests that, at least in the near field, dynamic triggering is the dominant cause of aftershocks, and that it generates both immediate and delayed aftershock activity.

Geophysical Research Letters↗
USGS2006

Inference of postseismic deformation mechanisms of the 1923 Kanto earthquake

Coseismic slip associated with the M7.9, 1923 Kanto earthquake is fairly well understood, involving slip of up to 8 m along the Philippine Sea-Honshu interplate boundary under Sagami Bay and its onland extension. Postseismic deformation after the 1923 earthquake, however, is relatively poorly understood. We revisit the available deformation data in order to constrain possible mechanisms of postseismic deformation and to examine the consequences for associated stress changes in the surrounding crust. Data from two leveling lines and one tide gage station over the first 7-8 years postseismic period are of much greater amplitude than the corresponding expected interseismic deformation during the same period, making these data suitable for isolating the signal from postseismic deformation. We consider both viscoelastic models of asthenosphere relaxation and afterslip models. A distributed coseismic slip model presented by Pollitz et al. (2005), combined with prescribed parameters of a viscoelastic Earth model, yields predicted postseismic deformation that agrees with observed deformation on mainland Honshu from Tokyo to the Izu peninsula. Elsewhere (southern Miura peninsula; Boso peninsula), the considered viscoelastic models fail to predict observed deformation, and a model of ???1 in shallow afterslip in the offshore region south of the Boso peninsula, with equivalent moment magnitude Mw = 7.0, adequately accounts for the observed deformation. Using the distributed coseismic slip model, layered viscoelastic structure, and a model of interseismic strain accumulation, we evaluate the post-1923 stress evolution, including both the coseismic and accumulated postseismic stress changes and those stresses contributed by interseismic loading. We find that if account is made for the varying tectonic regime in the region, the occurrence of both immediate (first month) post-1923 crustal aftershocks as well as recent regional crustal seismicity is consistent with the predicted stress pattern. This suggests that the influence of the 1923 earthquake on regional seismicity is fairly predictable and has persisted for at least seven decades following the earthquake.

Journal of Geophysical Research B: Solid Earth↗
USGS2006

Stress changes along the Sunda trench following the 26 December 2004 Sumatra-Andaman and 28 March 2005 Nias earthquakes

The 26 December 2004 Mw = 9.2 and 28 March 2005 Mw = 8.7 earthquakes on the Sumatra megathrust altered the state of stress over a large region surrounding the earthquakes. We evaluate the stress changes associated with coseismic and postseismic deformation following these two large events, focusing on postseismic deformation that is driven by viscoelastic relaxation of a low-viscosity asthenosphere. Under Coulomb failure stress (CFS) theory, the December 2004 event increased CFS on the future hypocentral zone of the March 2005 event by about 0.25 bar, with little or no contribution from viscous relaxation. Coseismic stresses around the rupture zones of the 1797 and 1833 Sunda trench events are negligible, but postseismic stress perturbations since December 2004 are predicted to result in CFS increases of 0.1 to 0.2 bar around these rupture zones between 2 and 8 years after the December 2004 event. These are considerable stress perturbations given that the 1797 and 1833 rupture zones are likely approaching the end of a complete seismic cycle. Copyright 2006 by the American Geophysical Union.

Geophysical Research Letters↗
USGS2006

Post-seismic relaxation following the great 2004 Sumatra-Andaman earthquake on a compressible self-gravitating Earth

he Mw ??? 9.0 2004 December 26 Sumatra-Andaman and Mw =8.7 2005 March 28 Nias earthquakes, which collectively ruptured approximately 1800 km of the Andaman and Sunda subduction zones, are expected to be followed by vigorous viscoelastic relaxation involving both the upper and lower mantle. Because of these large spatial dimensions it is desirable to fully account for gravitational coupling effects in the relaxation process. We present a stable method of computing relaxation of a spherically-stratified, compressible and self-gravitating viscoelastic Earth following an impulsive moment release event. The solution is cast in terms of a spherical harmonic expansion of viscoelastic normal modes. For simple layered viscoelastic models, which include a low-viscosity oceanic asthenosphere, we predict substantial post-seismic effects over a region several 100s of km wide surrounding the eastern Indian Ocean. We compare observed GPS time-series from ten regional sites (mostly in Thailand and Indonesia), beginning in 2004 December, with synthetic time-series that include the coseismic and post-seismic effects of the 2004 December 26 and 2005 March 28 earthquakes. A viscosity structure involving a biviscous (Burgers body) rheology in the asthenosphere explains the pattern and amplitude of post-seismic offsets remarkably well. ?? 2006 The Authors Journal compilation ?? 2006 RAS.

Geophysical Journal International↗
USGS2005

Postseismic deformation following the June 2000 earthquake sequence in the south Iceland seismic zone

We observe postseismic deformation on two spatiotemporal scales following Mw = 6.5 earthquakes in the south Iceland seismic zone on 17 and 21 June 2000. We see a rapidly decaying deformation transient lasting no more than 2 months and extending about 5 km away from the two main shock ruptures. This local, month-scale transient is captured by several radar interferograms and is also observed at a few campaign GPS sites located near the faults. A slower transient with a characteristic timescale of about a year is detected only by GPS measurements. The month-scale deformation pattern has been explained by poroelastic rebound due to postearthquake pore pressure changes. In contrast, the year-scale deformation can be explained by either afterslip at 8-14 km depth or viscoelastic relaxation of the lower crust and upper mantle in response to the coseismic stress changes. The optimal viscoelastic models have lower crustal viscosities of 0.5-1 ?? 1019 Pa s and upper mantle viscosity of ???3 ?? 1018 Pa s. Because of the limitations of our GPS campaign data, we consider both afterslip and viscoelastic relaxation as plausible mechanisms explaining the deformation field. Both types of postseismic deformation models suggest that the areas of large coseismic stress increase east of the 17 June and west of the 21 June ruptures continue to be loaded by the postseismic deformation. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2005

Coseismic slip distribution of the 1923 Kanto earthquake, Japan

The slip distribution associated with the 1923 M = 7.9 Kanto, Japan, earthquake is reexamined in light of new data and modeling. We utilize a combination of first-order triangulation, second-order triangulation, and leveling data in order to constrain the coseismic deformation. The second-order triangulation data, which have not been utilized in previous studies of 1923 coseismic deformation, are associated with only slightly smaller errors than the first-order triangulation data and expand the available triangulation data set by about a factor of 10. Interpretation of these data in terms of uniform-slip models in a companion study by Nyst et al. shows that a model involving uniform coseismic slip on two distinct rupture planes explains the data very well and matches or exceeds the fit obtained by previous studies, even one which involved distributed slip. Using the geometry of the Nyst et al. two-plane slip model, we perform inversions of the same geodetic data set for distributed slip. Our preferred model of distributed slip on the Philippine Sea plate interface has a moment magnitude of 7.86. We find slip maxima of ???8-9 m beneath Odawara and ???7-8 m beneath the Miura peninsula, with a roughly 2:1 ratio of strike-slip to dip-slip motion, in agreement with a previous study. However, the Miura slip maximum is imaged as a more broadly extended feature in our study, with the high-slip region continuing from the Miura peninsula to the southern Boso peninsula region. The second-order triangulation data provide good evidence for ???3 m right-lateral strike slip on a 35-km-long splay structure occupying the volume between the upper surface of the descending Philippine Sea plate and the southern Boso peninsula. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2005

Waveform tomography of crustal structure in the south San Francisco Bay region

We utilize a scattering-based seismic tomography technique to constrain crustal tructure around the southern San Francisco Bay region (SFBR). This technique is based on coupled traveling wave scattering theory, which has usually been applied to the interpretation of surface waves in large regional-scale studies. Using fully three-dimensional kernels, this technique is here applied to observed P, S, and surface waves of intermediate period (3-4 s dominant period) observed following eight selected regional events. We use a total of 73 seismograms recorded by a U.S. Geological Survey short-period seismic array in the western Santa Clara Valley, the Berkeley Digital Seismic Network, and the Northern California Seismic Network. Modifications of observed waveforms due to scattering from crustal structure include (positive or negative) amplification, delay, and generation of coda waves. The derived crustal structure explains many of the observed signals which cannot be explained with a simple layered structure. There is sufficient sensitivity to both deep and shallow crustal structure that even with the few sources employed in the present study, we obtain shallow velocity structure which is reasonably consistent with previous P wave tomography results. We find a depth-dependent lateral velocity contrast across the San Andreas fault (SAF), with higher velocities southwest of the SAF in the shallow crust and higher velocities northeast of the SAF in the midcrust. The method does not have the resolution to identify very slow sediment velocities in the upper approximately 3 km since the tomographic models are smooth at a vertical scale of about 5 km. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗