Temporal evolution of continental lithospheric strength in actively deforming regions,GSA

ABSTRACT 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 10 3 yr) and glacio-isostatic adjustment (GIA, ~10 3 to 10 4 yr) time scales. Topography and gravity field correlations indicate that lithosphere isostatic adjustment (LIA) on ~10 6 -10 7 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

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

International audienceWe 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

GPS Measurements Athe Northern Caribbean Plate Boundary Zone: Impact of Postseismic Relaxation Following Historic Earthquakes

GPS measurements in the northern Caribbean suggest that the rate of Caribbean plate motion relative to North America is about 10 mm/yr faster than predicted by global plate motion model NUVEL-1A. Several of the key sites used in the GPS study are located in the Dominican Republic, near the rupture zones of large earthquakes in 1946 and in the previous two centuries. Postseismic relaxation of the crust and upper mantle is a possible explanation for the plate velocity discrepancy. We explore a range of fault mechanisms and crustal and mantle rheology to place an upper bound on postseismic relaxation effects. The upper bound velocity contribution in the southern Dominican Republic is 5–6 mm/yr, and the most plausible contribution is 1–2 mm/yr, suggesting that postseismic effects cannot account for the discrepancy. This implies that the NUVEL-1A model underestimates the rate of motion of the Caribbean plate

Geodetic slip model of the 2011 M9.0 Tohoku earthquake

The three-dimensional crustal displacement field as sampled by GPS is used to determine the coseismic slip of the 2011 M9.0 Tohoku Earthquake. We employ a spherically layered Earth structure and use a combination of onland GPS, out to ∼4000 km from the rupture, and offshore GPS, which samples the high-slip region on the interplate boundary along the Japan trench. Inversion of the displacement field for dip slip, assuming an interplate boundary of variable dip and striking 195°, yields a compact slip maximum of about 33 m located 200 km east of Sendai. The geodetic moment is 4.06 × 10^22 N m, corresponding to Mw = 9.0. The area of maximum slip is concentrated at a depth of about 10 km, is updip of the rupture areas of the M≳7 Miyagi-oki earthquakes of 1933, 1936, 1937, and 1978, and roughly coincides with the rupture area of the M7.1 1981 Miyagi-oki earthquake. The overlap of the 2011 slip area with several preceding ruptures suggests that the same asperities may rupture repeatedly with M≳7 events within several decades of one another.Published versio

Coseismic slip distribution of the February 27, 2010 Mw 8.8 Maule, Chile earthquake

International audienceStatic offsets produced by the February 27, 2010 Mw = 8.8 Maule, Chile earthquake as measured by GPS and InSAR constrain coseismic slip along a section of the Andean megathrust of dimensions 650 km (in length) × 180 km (in width). GPS data have been collected from both campaign and continuous sites sampling both the near-field and far field. ALOS/PALSAR data from several ascending and descending tracks constrain the near-field crustal deformation. Inversions of the geodetic data for distributed slip on the megathrust reveal a pronounced slip maximum of order 15 m at ∼15-25 km depth on the megathrust offshore Lloca, indicating that seismic slip was greatest north of the epicenter of the bilaterally propagating rupture. A secondary slip maximum appears at depth ∼25 km on the megathrust just west of Concepción. Coseismic slip is negligible below 35 km depth. Estimates of the seismic moment based on different datasets and modeling approaches vary from 1.8 to 2.6 × 1022 N m. Our study is the first to model the static displacement field using a layered spherical Earth model, allowing us to incorporate both near-field and far-field static displacements in a consistent manner. The obtained seismic moment of 1.97 × 1022 N m, corresponding to a moment magnitude of 8.8, is similar to that obtained by previous seismic and geodetic inversions