Geology topics

F. Pollitz

Publications and source records attributed to F. Pollitz.

12 recordsLinked to original sources

USGS2011

High-frequency Born synthetic seismograms based on coupled normal modes

High-frequency and full waveform synthetic seismograms on a 3-D laterally heterogeneous earth model are simulated using the theory of coupled normal modes. The set of coupled integral equations that describe the 3-D response are simplified into a set of uncoupled integral equations by using the Born approximation to calculate scattered wavefields and the pure-path approximation to modulate the phase of incident and scattered wavefields. This depends upon a decomposition of the aspherical structure into smooth and rough components. The uncoupled integral equations are discretized and solved in the frequency domain, and time domain results are obtained by inverse Fourier transform. Examples show the utility of the normal mode approach to synthesize the seismic wavefields resulting from interaction with a combination of rough and smooth structural heterogeneities. This approach is applied to an ∼4 Hz shallow crustal wave propagation around the site of the San Andreas Fault Observatory at Depth (SAFOD).

Geophysical Journal International↗

USGS2011

Lower crustal relaxation beneath the Tibetan Plateau and Qaidam Basin following the 2001 Kokoxili earthquake

In 2001 November a magnitude 7.8 earthquake ruptured a 400 km long portion of the Kunlun fault, northeastern Tibet. In this study, we analyse over five years of post-seismic geodetic data and interpret the observed surface deformation in terms of stress relaxation in the thick Tibetan lower crust. We model GPS time-series (first year) and InSAR line of sight measurements (years two to five) and infer that the most likely mechanism of post-seismic stress relaxation is time-dependent distributed creep of viscoelastic material in the lower crust. Since a single relaxation time is not sufficient to model the observed deformation, viscous flow is modelled by a lower crustal Burgers rheology, which has two material relaxation times. The optimum model has a transient viscosity 9 × 10 17 Pa s, steady-state viscosity 1 × 10 19 Pa s and a ratio of long term to Maxwell shear modulus of 2:3. This model gives a good fit to GPS stations south of the Kunlun Fault, while displacements at stations north of the fault are over-predicted. We attribute this asymmetry in the GPS residual to lateral heterogeneity in rheological structure across the southern margin of the Qaidam Basin, with thinner crust/higher viscosities beneath the basin than beneath the Tibetan Plateau. Deep afterslip localized in a shear zone beneath the fault rupture gives a reasonable match to the observed InSAR data, but the slip model does not fit the earlier GPS data well. We conclude that while some localized afterslip likely occurred during the early post-seismic phase, the bulk of the observed deformation signal is due to viscous flow in the lower crust. To investigate regional variability in rheological structure, we also analyse post-seismic displacements following the 1997 Manyi earthquake that occurred 250 km west of the Kokoxili rupture. We find that viscoelastic properties are the same as for the Kokoxili area except for the transient viscosity, which is 5 × 10 17 Pa s. The viscosities estimated for the Manyi and Kokoxili areas are consistent with constraints obtained from other earthquakes in the northwest and south central parts of the Tibetan Plateau.

Tibetan Plateau and Qaidam Basin↗

USGS2011

Lithosphere-asthenosphere interaction beneath the western United States from the joint inversion of body-wave traveltimes and surface-wave phase velocities

The relation between the complex geological history of the western margin of the North American plate and the processes in the mantle is still not fully documented and understood. Several pre-USArray local seismic studies showed how the characteristics of key geological features such as the Colorado Plateau and the Yellowstone Snake River Plains are linked to their deep mantle structure. Recent body-wave models based on the deployment of the high density, large aperture USArray have provided far more details on the mantle structure while surface-wave tomography (ballistic waves and noise correlations) informs us on the shallow structure. Here we combine constraints from these two data sets to image and study the link between the geology of the western United States, the shallow structure of the Earth and the convective processes in mantle. Our multiphase DNA10-S model provides new constraints on the extent of the Archean lithosphere imaged as a large, deeply rooted fast body that encompasses the stable Great Plains and a large portion of the Northern and Central Rocky Mountains. Widespread slow anomalies are found in the lower crust and upper mantle, suggesting that low-density rocks isostatically sustain part of the high topography of the western United States. The Yellowstone anomaly is imaged as a large slow body rising from the lower mantle, intruding the overlying lithosphere and controlling locally the seismicity and the topography. The large E-W extent of the USArray used in this study allows imaging the 'slab graveyard', a sequence of Farallon fragments aligned with the currently subducting Juan de Fuca Slab, north of the Mendocino Triple Junction. The lithospheric root of the Colorado Plateau has apparently been weakened and partly removed through dripping. The distribution of the slower regions around the Colorado Plateau and other rigid blocks follows closely the trend of Cenozoic volcanic fields and ancient lithospheric sutures, suggesting that the later exert a control on the locus of magmato-tectonic activity today. The DNA velocity models are available for download and slicing at http://dna.berkeley.edu. ?? 2011 The Authors Geophysical Journal International ?? 2011 RAS.

Geophysical Journal International↗

USGS2008

Effect of 3-D viscoelastic structure on post-seismic relaxation from the 2004 M = 9.2 Sumatra earthquake

The 2004 M=9.2 Sumatra-Andaman earthquake profoundly altered the state of stress in a large volume surrounding the ???1400 km long rupture. Induced mantle flow fields and coupled surface deformation are sensitive to the 3-D rheology structure. To predict the post-seismic motions from this earthquake, relaxation of a 3-D spherical viscoelastic earth model is simulated using the theory of coupled normal modes. The quasi-static deformation basis set and solution on the 3-D model is constructed using: a spherically stratified viscoelastic earth model with a linear stress-strain relation; an aspherical perturbation in viscoelastic structure; a 'static'mode basis set consisting of Earth's spheroidal and toroidal free oscillations; a "viscoelastic" mode basis set; and interaction kernels that describe the coupling among viscoelastic and static modes. Application to the 2004 Sumatra-Andaman earthquake illustrates the profound modification of the post-seismic flow field at depth by a slab structure and similarly large effects on the near-field post-seismic deformation field at Earth's surface. Comparison with post-seismic GPS observations illustrates the extent to which viscoelastic relaxation contributes to the regional post-seismic deformation. ?? Journal compilation ?? 2008 RAS.

Geophysical Journal International↗

USGS2007

Coseismic slip distributions of the 26 December 2004 Sumatra-Andaman and 28 March 2005 Nias earthquakes from GPS static offsets

Static offsets produced by the 26 December 2004 M ???9 Sumatra-Andaman earthquake as measured by Global Positioning System (GPS) reveal a large amount of slip along the entire ???1300 km-long rupture. Most seismic slip inversions place little slip on the Andaman segment. whereas both near-field and far-field GPS offsets demand large slip on the Andaman segment. We compile available datasets of the static offset to render a more detailed picture of the static-slip distribution. We construct geodetic offsets such that postearthquake positions of continuous GPS sites are reckoned to a time 1 day after the earthquake and campaign GPS sites are similarly corrected for postseismic motions. The newly revised slip distribution (Mw 9.22) reveals substantial segmentation of slip along the Andaman Islands, with the southern quarter slipping ???15 m in unison with the adjacent Nicobar and northern Sumatran segments of length ???700 km. We infer a small excess of geodetic moment relative to the seismic moment. A similar compilation of GPS offsets from the 28 March 2005 Nias earthquake is well explained with dip slip averaging several meters (Mw = 8.66) distributed primarily at depths greater than 20 km.

Bulletin of the Seismological Society of America↗

USGS2007

Coseismic and post-seismic signatures of the Sumatra 2004 December and 2005 March earthquakes in GRACE satellite gravity

The GRACE satellite mission has been measuring the Earth's gravity field and its temporal variations since 2002 April. Although these variations are mainly due to mass transfer within the geofluid envelops, they also result from mass displacements associated with phenomena including glacial isostatic adjustment and earthquakes. However, these last contributions are difficult to isolate because of the presence of noise and of geofluid signals, and because of GRACE's coarse spatial resolution (>400 km half-wavelength). In this paper, we show that a wavelet analysis on the sphere helps to retrieve earthquake signatures from GRACE geoid products. Using a wavelet analysis of GRACE geoids products, we show that the geoid variations caused by the 2004 December ( M w = 9.2) and 2005 March ( M w = 8.7) Sumatra earthquakes can be detected. At GRACE resolution, the 2004 December earthquake produced a strong coseismic decrease of the gravity field in the Andaman Sea, followed by relaxation in the area affected by both the Andaman 2004 and the Nias 2005 earthquakes. We find two characteristic timescales for the relaxation, with a fast variation occurring in the vicinity of the Central Andaman ridge. We discuss our coseismic observations in terms of density changes of crustal and upper-mantle rocks, and of the vertical displacements in the Andaman Sea. We interpret the post-seismic signal in terms of the viscoelastic response of the Earth's mantle. The transient component of the relaxation may indicate the presence of hot, viscous material beneath the active Central Andaman Basin.

Geophysical Journal International↗

USGS2006

Body-wave traveltime and amplitude shifts from asymptotic travelling wave coupling

We explore the sensitivity of finite-frequency body-wave traveltimes and amplitudes to perturbations in 3-D seismic velocity structure relative to a spherically symmetric model. Using the approach of coupled travelling wave theory, we consider the effect of a structural perturbation on an isolated portion of the seismogram. By convolving the spectrum of the differential seismogram with the spectrum of a narrow window taper, and using a Taylor's series expansion for wavenumber as a function of frequency on a mode dispersion branch, we derive semi-analytic expressions for the sensitivity kernels. Far-field effects of wave interactions with the free surface or internal discontinuities are implicitly included, as are wave conversions upon scattering. The kernels may be computed rapidly for the purpose of structural inversions. We give examples of traveltime sensitivity kernels for regional wave propagation at 1 Hz. For the direct SV wave in a simple crustal velocity model, they are generally complicated because of interfering waves generated by interactions with the free surface and the Mohorovic??ic?? discontinuity. A large part of the interference effects may be eliminated by restricting the travelling wave basis set to those waves within a certain range of horizontal phase velocity. ?? Journal compilation ?? 2006 RAS.

Geophysical Journal International↗

USGS2004

A physical model for strain accumulation in the San Francisco Bay region: Stress evolution since 1838

Understanding of the behavior of plate boundary zones has progressed to the point where reasonably comprehensive physical models can predict their evolution. The San Andreas fault system in the San Francisco Bay region (SFBR) is dominated by a few major faults whose behavior over about one earthquake cycle is fairly well understood. By combining the past history of large ruptures on SFBR faults with a recently proposed physical model of strain accumulation in the SFBR, we derive the evolution of regional stress from 1838 until the present. This effort depends on (1) an existing compilation of the source properties of historic and contemporary SFBR earthquakes based on documented shaking, geodetic data, and seismic data (Bakun, 1999) and (2) a few key parameters of a simple regional viscoelastic coupling model constrained by recent GPS data (Pollitz and Nyst, 2004). Although uncertainties abound in the location, magnitude, and fault geometries of historic ruptures and the physical model relies on gross simplifications, the resulting stress evolution model is sufficiently detailed to provide a useful window into the past stress history. In the framework of Coulomb failure stress, we find that virtually all M ??? 5.8 earthquakes prior to 1906 and M ??? 5.5 earthquakes after 1906 are consistent with stress triggering from previous earthquakes. These events systematically lie in zones of predicted stress concentration elevated 5-10 bars above the regional average. The SFBR is predicted to have emerged from the 1906 "shadow" in about 1980, consistent with the acceleration in regional seismicity at that time. The stress evolution model may be a reliable indicator of the most likely areas to experience M ??? 5.5 shocks in the future.

Journal of Geophysical Research B: Solid Earth↗

USGS2004

Postseismic deformation and stress changes following the 1819 Rann of Kachchh, India earthquake: Was the 2001 Bhuj earthquake a triggered event?

The 2001 Mw 7.6 Bhuj earthquake occurred in an intraplate region with rather unusual active seismicity, including an earlier major earthquake, the 1819 Rann of Kachchh earthquake (M7.7). We examine if static coseismic and transient postseismic deformation following the 1819 earthquake contributed to the enhanced seismicity in the region and the occurrence of the 2001 Bhuj earthquake, ???100 km away and almost two centuries later. Based on the Indian shield setting, great rupture depth of the 2001 event and lack of significant early postseismic deformation measured following the 2001 event, we infer that little viscous relaxation occurs in the lower crust and choose an upper mantle effective viscosity of 1019 Pas. The predicted Coulomb failure stress (DCFS) on the rupture plane of the 2001 event increased by more than 0.1 bar at 20 km depth, which is a small but possibly significant amount. Stress change from the 1819 event may have also affected the occurrence of other historic earthquakes in this region. We also evaluate the postseismic deformation and ??CFS in this region due to the 2001 event. Positive ??CFS from the 2001 event occur to the NW and SE of the Bhuj earthquake rupture. Copyright 2004 by the American Geophysical Union.

Geophysical Research Letters↗

USGS2003

Motion of the Scotia sea plates

Earthquake data from the Scotia Arc to early 2002 are reviewed in the light of satellite gravity and other data in order to derive a model for the motion of plates in the Scotia Sea region. Events with magnitude ???5, which occurred on or near the boundaries of the Scotia and Sandwich plates, and for which Centroid Moment Tensor (CMT) solutions are available, are examined. The newer data fill some of the previous sampling gaps along the boundaries of the Scotia and Sandwich plates, and provide tighter constraints on relative motions. Variations in the width of the Brunhes anomaly on evenly spaced marine magnetic profiles over the East Scotia Ridge provide new estimates of Scotia-Sandwich plate spreading rates. Since there are no stable fracture zones in the east Scotia Sea, the mean azimuth of sea floor fabric mapped by sidescan is used to constrain the direction of spreading. 18 new rate estimates and four azimuths from the East Scotia Ridge are combined with 68 selected earthquake slip vectors from the boundaries of the Scotia Sea in a least-squares inversion for the best-fitting set of Euler poles and angular rotation rates describing the 'present-day' motions of the Scotia and Sandwich plates relative to South America and Antarctica. Our preferred model (TLP2003) gives poles that are similar to previous estimates, except for Scotia Plate motion with respect to South America, which is significantly different from earlier estimates; predicted rates of motion also differ slightly. Our results are much more robust than earlier work. We examine the implications of the model for motion and deformation along the various plate boundaries, with particular reference to the North and South Scotia Ridges, where rates are obtained by closure.

Geophysical Journal International↗

USGS2003

Constraints on the viscosity of the continental crust and mantle from GPS measurements and postseismic deformation models in western Mongolia

We use GPS measurements and models of postseismic deformation caused by seven M6.8 to 8.4 earthquakes that occurred in the past 100 years in Mongolia to assess the viscosity of the lower crust and upper mantle. We find an upper mantle viscosity between 1 ?? 1018 and 4 ?? 1018 Pa s. The presence of such a weak mantle is consistent with results from independent seismological and petrological studies that show an abnormally hot upper mantle beneath Mongolia. The viscosity of the lower crust is less well constrained, but a weak lower crust (3 ?? 1016 to 2 ?? 1017 Pa s) is preferred by the data. Using our best fit upper mantle and lower crust viscosities, we find that the postseismic effects of viscoelastic relaxation on present-day horizontal GPS velocities are small (<2 mm yr-1) but still persist 100 years after the 1905, M8.4, Bolnay earthquake. This study shows that the GPS velocity field in the Baikal-Mongolia area can be modeled as the sum of (1) a rigid translation and rotation of the whole network, (2) a 3-5 mm yr-1 simple shear velocity gradient between the Siberian platform to the north and northern China to the south, and (3) the contribution of postseismic deformation, mostly caused by the 1905 Bolnay-Tsetserleg sequence and by the smaller, but more recent, 1957 Bogd earthquake. Copyright 2003 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

USGS2003

Fault interaction and stress triggering of twentieth century earthquakes in Mongolia

A cluster of exceptionally large earthquakes in the interior of Asia occurred from 1905 to 1967: the 1905 M7.9 Tsetserleg and M8.4 Bolnai earthquakes, the 1931 M8.0 Fu Yun earthquake, the 1957 M8.1 Gobi-Altai earthquake, and the 1967 M7.1 Mogod earthquake (sequence). Each of the larger (M ??? 8) earthquakes involved strike-slip faulting averaging more than 5 m and rupture lengths of several hundred kilometers. Available geologic data indicate that recurrence intervals on the major source faults are several thousands of years and distances of about 400 km separate the respective rupture areas. We propose that the occurrences of these and many smaller earthquakes are related and controlled to a large extent by stress changes generated by the compounded static deformation of the preceding earthquakes and subsequent viscoelastic relaxation of the lower crust and upper mantle beneath Mongolia. We employ a spherically layered viscoelastic model constrained by the 1994-2002 GPS velocity field in western Mongolia [Vergnolle et al., 2003]. Using the succession of twentieth century earthquakes as sources of deformation, we then analyze the time-dependent change in Coulomb failure stress (????f). At remote interaction distances, static ????f values are small. However, modeled postseismic stress changes typically accumulate to several tenths of a bar over time intervals of decades. Almost all significant twentieth century regional earthquakes (M ??? 6) with well-constrained fault geometry lie in positive ????f lobes of magnitude about +0.5 bar. Our results suggest that significant stress transfer is possible among continental faults separated by hundreds of kilometers and on timescales of decades. Copyright 2003 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗