Geology topics

F. F. Pollitz

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

30 records · Page 2Linked to original sources

USGS2005

Transient rheology of the upper mantle beneath central Alaska inferred from the crustal velocity field following the 2002 Denali earthquake

The M7.9 2002 Denali earthquake, Alaska, is one of the largest strike-slip earthquakes ever recorded. The postseismic GPS velocity field around the 300-km-long rupture is characterized by very rapid horizontal velocity up to ???300 mm/yr for the first 0.1 years and slower but still elevated horizontal velocity up to ???100 mm/yr for the succeeding 1.5 years. I find that the spatial and temporal pattern of the displacement field may be explained by a transient mantle rheology. Representing the regional upper mantle as a Burghers body, I infer steady state and transient viscosities of ??1 = 2.8 ?? 1018 Pa s and ??2 = 1.0 ?? 1017 Pa s, respectively, corresponding to material relaxation times of 1.3 and 0.05 years. The lower crustal viscosity is poorly constrained by the considered horizontal velocity field, and the quoted mantle viscosities assume a steady state lower crust viscosity that is 7??1. Systematic bias in predicted versus observed velocity vectors with respect to a fixed North America during the first 3-6 months following the earthquake is reduced when all velocity vectors are referred to a fixed site. This suggests that the post-Denali GPS time series for the first 1.63 years are shaped by a combination of a common mode noise source during the first 3-6 months plus viscoelastic relaxation controlled by a transient mantle rheology.

Journal of Geophysical Research B: Solid Earth↗

USGS2005

Geophysics: The size and duration of the Sumatra-Andaman earthquake from far-field static offsets

The 26 December 2004 Sumatra earthquake produced static offsets at continuously operating GPS stations at distances of up to 4500 kilometers from the epicenter. We used these displacements to model the earthquake and include consideration of the Earth's shape and depth-varying rigidity. The results imply that the average slip was >5 meters along the full length of the rupture, including the ???650-kilometer-long Andaman segment. Comparison of the source derived from the far-field static offsets with seismically derived estimates suggests that 25 to 35% of the total moment release occurred at periods greater than 1 hour. Taking into consideration the strong dip dependence of moment estimates, the magnitude of the earthquake did not exceed Mw = 9.2.

Science↗

USGS2004

A physical model for strain accumulation in the San Francisco Bay Region

Strain accumulation in tectonically active regions is generally a superposition of the effects of background tectonic loading, steady-state dislocation processes, such as creep, and transient deformation. In the San Francisco Bay region (SFBR), the most uncertain of these processes is transient deformation, which arises primarily in association with large earthquakes. As such, it depends upon the history of faulting and the rheology of the crust and mantle, which together determine the pattern of longer term (decade-scale) post-seismic response to earthquakes. We utilize a set of 102 GPS velocity vectors in the SFBR in order to characterize the strain rate field and construct a physical model of its present deformation. We first perform an inversion for the continuous velocity gradient field from the discrete GPS velocity field, from which both tensor strain rate and rotation rate may be extracted. The present strain rate pattern is well described as a nearly uniform shear strain rate oriented approximately N34??W (140 nanostrain yr-1) plus a N56??E uniaxial compression rate averaging 20 nanostrain yr-1 across the shear zone. We fit the velocity and strain rate fields to a model of time-dependent deformation within a 135-kin-wide, arcuate shear zone bounded by strong Pacific Plate and Sierra Nevada block lithosphere to the SW and NE, respectively. Driving forces are purely lateral, consisting of shear zone deformation imposed by the relative motions between the thick Pacific Plate and Sierra Nevada block lithospheres. Assuming a depth-dependent viscoelastic structure within the shear zone, we account for the effects of steady creep on faults and viscoelastic relaxation following the 1906 San Francisco and 1989 Loma Prieta earthquakes, subject to constant velocity boundary conditions on the edges of the shear zone. Fault creep is realized by evaluating dislocations on the creeping portions of faults in the fluid limit of the viscoelastic model. A priori plate-boundary(PB)-parallel motion is set to 38 mm yr -1. A grid search based on fitting the observed strain rate pattern yields a mantle viscosity of 1.2 ?? 1019 Pa s and a PB-perpendicular convergence rate of ???3 mm yr-1. Most of this convergence appears to be uniformly distributed in the Pacific-Sierra Nevada plate boundary zone. ?? 2005 RAS.

Geophysical Journal International↗

USGS2003

Transient rheology of the uppermost mantle beneath the Mojave Desert, California

Geodetic data indicate that the M7.1 Hector Mine, California, earthquake was followed by a brief period (a few weeks) of rapid deformation preceding a prolonged phase of slower deformation. We find that the signal contained in continuous and campaign global positioning system data for 2.5 years after the earthquake may be explained with a transient rheology. Quantitative modeling of these data with allowance for transient (linear biviscous) rheology in the lower crust and upper mantle demonstrates that transient rheology in the upper mantle is dominant, its material properties being typified by two characteristic relaxation times ???0.07 and ???2 years. The inferred mantle rheology is a Jeffreys solid in which the transient and steady-state shear moduli are equal. Consideration of a simpler viscoelastic model with a linear univiscous rheology (2 fewer parameters than a biviscous model) shows that it consistently underpredicts the amplitude of the first ???3 months signal, and allowance for a biviscous rheology is significant at the 99.0% confidence level. Another alternative model - deep postseismic afterslip beneath the coseismic rupture - predicts a vertical velocity pattern opposite to the observed pattern at all time periods considered. Despite its plausibility, the advocated biviscous rheology model is non-unique and should be regarded as a viable alternative to the non-linear mantle rheology model for governing postseismic flow beneath the Mojave Desert. Published by Elsevier B.V.

Earth and Planetary Science Letters↗

USGS2003

The relationship between the instantaneous velocity field and the rate of moment release in the lithosphere

Instantaneous velocity gradients within the continental lithosphere are often related to the tectonic driving forces. This relationship is direct if the forces are secular, as for the case of loading of a locked section of a subduction interface by the downgoing plate. If the forces are static, as for the case of lateral variations in gravitational potential energy, then velocity gradients can be produced only if the lithosphere has, on average, zero strength. The static force model may be related to the long-term velocity field but not the instantaneous velocity field (typically measured geodetically over a period of several years) because over short time intervals the upper lithosphere behaves elastically. In order to describe both the short- and long-term behaviour of an (elastic) lithosphere-(viscoelastic) asthenosphere system in a self-consistent manner, I construct a deformation model termed the expected interseismic velocity (EIV) model. Assuming that the lithosphere is populated with faults that rupture continually, each with a definite mean recurrence time, and that the Earth is well approximated as a linear elastic-viscoelastic coupled system, I derive a simple relationship between the instantaneous velocity field and the average rate of moment release in the lithosphere. Examples with synthetic fault networks demonstrate that velocity gradients in actively deforming regions may to a large extent be the product of compounded viscoelastic relaxation from past earthquakes on hundreds of faults distributed over large ( ≥106 km2) areas.

Geophysical Journal International↗

USGS2003

Post-seismic relaxation theory on laterally heterogeneous viscoelastic model

Investigation was carried out into the problem of relaxation of a laterally heterogeneous viscoelastic Earth following an impulsive moment release event. The formal solution utilizes a semi-analytic solution for post-seismic deformation on a laterally homogeneous Earth constructed from viscoelastic normal modes, followed by application of mode coupling theory to derive the response on the aspherical Earth. The solution is constructed in the Laplace transform domain using the correspondence principle and is valid for any linear constitutive relationship between stress and strain. The specific implementation described in this paper is a semi-analytic discretization method which assumes isotropic elastic structure and a Maxwell constitutive relation. It accounts for viscoelastic-gravitational coupling under lateral variations in elastic parameters and viscosity. For a given viscoelastic structure and minimum wavelength scale, the computational effort involved with the numerical algorithm is proportional to the volume of the laterally heterogeneous region. Examples are presented of the calculation of post-seismic relaxation with a shallow, laterally heterogeneous volume following synthetic impulsive seismic events, and they illustrate the potentially large effect of regional 3-D heterogeneities on regional deformation patterns.

Geophysical Journal International↗

USGS2002

Regional seismic wavefield computation on a 3-D heterogeneous Earth model by means of coupled traveling wave synthesis

I present a new algorithm for calculating seismic wave propagation through a three-dimensional heterogeneous medium using the framework of mode coupling theory originally developed to perform very low frequency (f < ???0.01-0.05 Hz) seismic wavefield computation. It is a Greens function approach for multiple scattering within a defined volume and employs a truncated traveling wave basis set using the locked mode approximation. Interactions between incident and scattered wavefields are prescribed by mode coupling theory and account for the coupling among surface waves, body waves, and evanescent waves. The described algorithm is, in principle, applicable to global and regional wave propagation problems, but I focus on higher frequency (typically f ??????0.25 Hz) applications at regional and local distances where the locked mode approximation is best utilized and which involve wavefields strongly shaped by propagation through a highly heterogeneous crust. Synthetic examples are shown for P-SV-wave propagation through a semi-ellipsoidal basin and SH-wave propagation through a fault zone.

Pure and Applied Geophysics↗

USGS2002

Stress triggering of the 1999 Hector Mine earthquake by transient deformation following the 1992 Landers earthquake

The M 7.3 June 28, 1992 Landers and M 7.1 October 16, 1999 Hector Mine earthquakes, California, both right lateral strike-slip events on NNW-trending subvertical faults, occurred in close proximity in space and time in a region where recurrence times for surface-rupturing earthquakes are thousands of years. This suggests a causal role for the Landers earthquake in triggering the Hector Mine earthquake. Previous modeling of the static stress change associated with the Landers earthquake shows that the area of peak Hector Mine slip lies where the Coulomb failure stress promoting right-lateral strike-slip failure was high, but the nucleation point of the Hector Mine rupture was neutrally to weakly promoted, depending on the assumed coefficient of friction. Possible explanations that could account for the 7-year delay between the two ruptures include background tectonic stressing, dissipation of fluid pressure gradients, rate- and state-dependent friction effects, and post-Landers viscoelastic relaxation of the lower crust and upper mantle. By employing a viscoelastic model calibrated by geodetic data collected during the time period between the Landers and Hector Mine events, we calculate that postseismic relaxation produced a transient increase in Coulomb failure stress of about 0.7 bars on the impending Hector Mine rupture surface. The increase is greatest over the broad surface that includes the 1999 nucleation point and the site of peak slip further north. Since stress changes of magnitude greater than or equal to 0.1 bar are associated with documented causal fault interactions elsewhere, viscoelastic relaxation likely contributed to the triggering of the Hector Mine earthquake. This interpretation relies on the assumption that the faults occupying the central Mojave Desert (i.e., both the Landers and Hector Mine rupturing faults) were critically stressed just prior to the Landers earthquake.

California↗

USGS2001

Viscoelastic shear zone model of a strike-slip earthquake cycle

I examine the behavior of a two-dimensional (2-D) strike-slip fault system embedded in a 1-D elastic layer (schizosphere) overlying a uniform viscoelastic half-space (plastosphere) and within the boundaries of a finite width shear zone. The viscoelastic coupling model of Savage and Prescott [1978] considers the viscoelastic response of this system, in the absence of the shear zone boundaries, to an earthquake occurring within the upper elastic layer, steady slip beneath a prescribed depth, and the superposition of the responses of multiple earthquakes with characteristic slip occurring at regular intervals. So formulated, the viscoelastic coupling model predicts that sufficiently long after initiation of the system, (1) average fault-parallel velocity at any point is the average slip rate of that side of the fault and (2) far-field velocities equal the same constant rate. Because of the sensitivity to the mechanical properties of the schizosphere-plastosphere system (i.e., elastic layer thickness, plastosphere viscosity), this model has been used to infer such properties from measurements of interseismic velocity. Such inferences exploit the predicted behavior at a known time within the earthquake cycle. By modifying the viscoelastic coupling model to satisfy the additional constraint that the absolute velocity at prescribed shear zone boundaries is constant, I find that even though the time-averaged behavior remains the same, the spatiotemporal pattern of surface deformation (particularly its temporal variation within an earthquake cycle) is markedly different from that predicted by the conventional viscoelastic coupling model. These differences are magnified as plastosphere viscosity is reduced or as the recurrence interval of periodic earthquakes is lengthened. Application to the interseismic velocity field along the Mojave section of the San Andreas fault suggests that the region behaves mechanically like a ???600-km-wide shear zone accommodating 50 mm/yr fault-parallel motion distributed between the San Andreas fault system and Eastern California Shear Zone. Copyright 2001 by the American Geophysical Union.

California↗

USGS2001

Remarks on the travelling wave decomposition

In elastic wave propagation on a spherically symmetric earth model, a normal mode sum is converted into a sum of equivalent travelling waves by means of a travelling wave decomposition (TWD). For two decades, seismologists have assumed that each travelling wave in the TWD is associated with only real phase velocities, that is, no evanescent waves travel on a spherically symmetric earth model. In this paper, this assumption is proven false. By including a countably infinite set of waves travelling as evanescent waves, several conceptual difficulties confronting the TWD are resolved.

Geophysical Journal International↗

USGS2001

Mantle flow beneath a continental strike-slip fault: Postseismic deformation after the 1999 Hector Mine earthquake

Two recent large earthquakes in the Mojave Desert, California - the magnitude 7.3 1992 Landers and magnitude 7.1 1999 Hector Mine earthquakes - have each been followed by elevated crustal strain rates over periods of months and years. Geodetic data collected after the Hector Mine earthquake exhibit a temporally decaying horizontal velocity field and a quadrant uplift pattern opposite to that expected for localized shear beneath the earthquake rupture. We interpret the origin of this accelerated crustal deformation to be vigorous flow in the upper mantle in response to the stress changes generated by the earthquake. Our results suggest that transient flow in the upper mantle is a fundamental component of the earthquake cycle and that the lower crust is a coherent stress guide coupling the upper crust with the upper mantle.

Science↗

USGS2001

Sinking mafic body in a reactivated lower crust: A mechanism for stress concentration at the New Madrid seismic zone

We propose a geodynamic model for stress concentration in the New Madrid seismic zone (NMSZ). The model postulates that a high-density (mafic) body situated in the deep crust directly beneath the most seismically active part of the NMSZ began sinking several thousands of years ago when the lower crust was suddenly weakened. Based on the fact that deformation rates in the NMSZ have accelerated over the past 9 k.y., we envision the source of this perturbation to be related to the last North American deglaciation. Excess mass of the mafic body exerts a downward pull on the elastic upper crust, leading to a cycle of primary thrust faulting with secondary strike-slip faulting, after which continued sinking of the mafic body reloads the upper crust and renews the process. This model is consistent with the youth of activity, the generation of a sequence of earthquakes, and the velocity evolution during interseismic periods, which depend upon the density contrast of the mafic body with respect to the surrounding crust, its volume, and the viscosity of the lower crust.

Bulletin of the Seismological Society of America↗