Fred F. Pollitz

Gravitational viscoelastic postseismic relaxation on a layered spherical Earth

Viscoelastic relaxation of a ductile asthenosphere underlying a purely elastic plate is a strong candidate process for explaining anomalous rates of crustal deformation observed following large earthquakes. The nongravitational treatment of Pollitz [1992], which is valid on a global scale and includes the effects of compressibility, is here extended to permit the calculation of gravitational viscoelastic relaxation for a specified spherically layered viscoelastic rheology following an earthquake in an elastic layer. The simple approximations we adopt make the resulting treatment particularly suitable for near-field calculations. For an asthenosphere with a Maxwell rheology, the effect of gravitational coupling is manifested only many relaxation times after the earthquake, as obtained by previous investigators. Its effect is generally to speed up the long-wavelength component of the relaxation process and attenuate the overall vertical displacement pattern. Several subtle features common to the relaxation behavior from several different fault types (thrust, rift, and strike-slip) are identified. The effect of gravitational coupling on horizontal displacements is consistent with flexure of the upper elastic plate driven by the corresponding effects on the vertical displacement. Stress diffusion away from the source region generally exhibits pulse-like behavior which is dispersive in both space and time. If the asthenosphere is confined to a relatively narrow channel, then the dispersion branches governing relaxation are radically altered, and stress diffusion effects far from the coseismic rupture zone exhibit a complicated time dependence reflecting the competing tendencies of toroidal and spheroidal mode relaxation.

Mobility of continental mantle: Evidence from postseismic geodetic observations following the 1992 Landers earthquake

The crust around the rupture zone of the 1992 Landers earthquake has continued to deform in the years following the earthquake at rates ∼3 times greater than pre-earthquake rates. We use a combination of Global Positioning System (GPS) and synthetic aperture radar (InSAR) data collected during a ∼3-year epoch following the earthquake in order to investigate postseismic mechanisms responsible for the high transient velocities. In order to maximize the potential signal from viscoelastic relaxation we evaluate and model postseismic relaxation following the first few months of documented accelerated deformation. The combination of GPS and InSAR data allows us to establish viscoelastic relaxation of the lower crust and upper mantle as the dominant postseismic process and to discriminate among possible viscoelastic models. The data particularly require the presence of a highly ductile uppermost mantle beneath the central Mojave Domain, with temperature between the wet and dry basalt solidus. This is consistent with independent seismic and geochemical inferences of a regionally warm uppermost mantle. Further consideration of seismic velocity variations in conjunction with faulting patterns within the Mojave Desert suggests that the primary faulting characteristics of the Mojave Desert, namely, the pervasive late Cenozoic deformation within the Eastern California Shear Zone versus the near absence of faults in the Western Mojave Domain, are controlled by the rheology of the uppermost mantle.

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

A paper-machine press section having structure for transferring a web away from the press section toward the drying section of the machine. The press section includes a press roll with which a transfer structure cooperates for transferring a web from this press roll away from the latter and toward the drying section. This transfer structure includes at least a transfer suction roll having an outer shell formed with perforations through which suction is transmitted to a web. A suction chamber communicates with these perforations while between the suction chamber and the outer shell which has the perforations there is a sealing structure which is capable of providing between the perforated shell and the suction chamber a seal which does not require the use of seal water.

Geodetic and seismic constraints on some seismogenic zone processes in Costa Rica

New seismic and geodetic data from Costa Rica provide insight into seismogenic zone processes in Central America, where the Cocos and Caribbean plates converge. Seismic data are from combined land and ocean bottom deployments in the Nicoya peninsula in northern Costa Rica and near the Osa peninsula in southern Costa Rica. In Nicoya, inversion of GPS data suggests two locked patches centered at 14 ± 2 and 39 ± 6 km depth. Interplate microseismicity is concentrated in the more freely slipping intermediate zone, suggesting that small interseismic earthquakes may not accurately outline the updip limit of the seismogenic zone, the rupture zone for future large earthquakes, at least over the short (∼1 year) observation period. We also estimate northwest motion of a coastal “sliver block” at 8 ± 3 mm/yr, probably related to oblique convergence. In the Osa region to the south, convergence is orthogonal to the trench. Cocos-Caribbean relative motion is partitioned here, with ∼8 cm/yr on the Cocos-Panama block boundary (including a component of permanent shortening across the Fila Costena fold and thrust belt) and ∼1 cm/yr on the Panama block–Caribbean boundary. The GPS data suggest that the Cocos plate–Panama block boundary is completely locked from ∼10–50 km depth. This large locked zone, as well as associated forearc and back-arc deformation, may be related to subduction of the shallow Cocos Ridge and/or younger lithosphere compared to Nicoya, with consequent higher coupling and compressive stress in the direction of plate convergence.

The 11 April 2012 east Indian Ocean earthquake triggered large aftershocks worldwide

Large earthquakes trigger very small earthquakes globally during passage of the seismic waves and during the following several hours to days, but so far remote aftershocks of moment magnitude M ≥ 5.5 have not been identified, with the lone exception of an M = 6.9 quake remotely triggered by the surface waves from an M = 6.6 quake 4,800 kilometres away. The 2012 east Indian Ocean earthquake that had a moment magnitude of 8.6 is the largest strike-slip event ever recorded. Here we show that the rate of occurrence of remote M ≥ 5.5 earthquakes (>1,500 kilometres from the epicentre) increased nearly fivefold for six days after the 2012 event, and extended in magnitude to M ≤ 7. These global aftershocks were located along the four lobes of Love-wave radiation; all struck where the dynamic shear strain is calculated to exceed 10−7 for at least 100 seconds during dynamic-wave passage. The other M ≥ 8.5 mainshocks during the past decade are thrusts; after these events, the global rate of occurrence of remote M ≥ 5.5 events increased by about one-third the rate following the 2012 shock and lasted for only two days, a weaker but possibly real increase. We suggest that the unprecedented delayed triggering power of the 2012 earthquake may have arisen because of its strike-slip source geometry or because the event struck at a time of an unusually low global earthquake rate, perhaps increasing the number of nucleation sites that were very close to failure.

Joint estimation of afterslip rate and postseismic relaxation following the 1989 Loma Prieta earthquake

Global Positioning System (GPS) data from campaigns carried out over the 5 years following the 1989 Loma Prieta earthquake and leveling data measured in 1990 and 1992 define the postseismic velocity field around the Loma Prieta rupture zone. Subtraction of a background velocity field yields a residual velocity pattern which we interpret as the product of two physical processes: (1) slow afterslip along distinct planes in the upper cru, st and (2) viscoelastic relaxation of the lower crust and upper mantle. Biirgmann et al. (1997) previously derived an afterslip model involving uniform afterslip on two optimally determined planes, in- cluding oblique reverse slip on the coseismic rupture and reverse slip on a shallow thrust fault to the northeast of the San Andreas fault. We further consider models of distributed slip on these two fault planes plus a viscoelastic relaxation pattern which depends on a suitable cose- ismic rupture model and crust and mantle viscosities. Several fault models from the literature were considered for the 1989 coseismic rupture, with nearly identical impact on the results. Simultaneous maximum likelihood inversion of the GPS and leveling data for afterslip distri- bution and viscosity yields the following results: (1) A good fit to the data is obtained by smooth afterslip distributions without any viscoelastic relaxation being required. (2) Tangible broad-scale viscoelastic relaxation of the lower crust and upper mantle are present in this data set at 97% confidence, and a lower crustal viscosity of- 1029 Pa s is obtained; however, the viscosity of both the lower crust and upper mantle are poorly constrained. (3) For a given misfit, 20% less integrated afterslip and smoother afterslip distributions result when viscoelas- tic relaxation is included. (4) Maximum slip rates on the slip-distributed models are 3-5 crn/yr, the dominant patches estimated on the two planes fill in the entire depth range 4-13 km without significant overlap, and deeper afterslip is not required. The afterslip distribution on the coseismic rupture plane is strongly dominated by reverse slip immediately southeast of the main center of coseismic reverse slip.

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 sug- gests a causal role for the Landers earthquake in triggering the Hector Mine earth- quake. 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, dissi- pation 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 relaxa- tion 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 con- tributed 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.

Temporal evolution of continental lithospheric strength in actively deforming regions

It has been agreed for nearly a century that a strong, loadbearing 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 depth 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 0 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 substrate. 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.

Viscosity structure beneath northeast Iceland

The dynamics of crustal rifting in Iceland depend strongly on the lower crustal rheology, which controls the intensity of upper crustal stress concentration and scale time of heat diffusion from the underlying mantle plume. While magnetotelluric surveys suggest the presence of a pervasive hot and highly ductile lowermost crust with possibly high fraction of partial melt, observations of low seismic attenuation and strong shear wave transmission suggest a much cooler lower crust and upper mantle. Since viscosity is also sensitive to the degree of partial melt present, viscosity estimates for these regions could shed light on the factors responsible for these observations. In this study we utilize horizontal and vertical displacement vectors determined in GPS campaigns in northeast Iceland since 1986. These are modeled in terms of steady state tectonic loading plus postseismic/postdiking relaxation following the 1975–1984 Krafla rifting episode, as first proposed by Foulger and others. With the elastic part of the model fixed by external constraints, these data have a high sensitivity to the viscosity structure beneath Iceland. Lower crust and upper mantle viscosities of about 3 × 1019 Pa s and 3 × 1018 Pa s, respectively, yield the closest agreement with the data. Our lower crustal viscosity estimate is consistent with the low attenuation and low (subsolidus) temperature for the lower crust inferred in recent studies. Inversions for fissure opening during the Krafla rifting episode yield about 7 m of opening centered on the Krafla rift, as is observed. Allowing for contemporaneous deep rifting on vertical faults along the Askja segment partially accounts for the observed increase in separation across the rift during 1987–1992 but does not account for large displacements in the southeastern part of the network or the large relative subsidence around the Askja rift during 1987–1990. Recent deep normal faulting beneath the Askja rift and further south might explain all of these remaining features.

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, 1 infer steady state and transient viscosities of η 1 = 2.8 x 10 18 Pa s and η 2 = 1.0 x 10 17 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.