Earthquake Cycle Deformation in the Tibetan Plateau with a Weak Mid-Crustal Layer

Abstract

Geodetic observations of interseismic deformation across the Tibetan plateau contain information about both tectonic and earthquake cycle processes. Time-variations in surface velocities between large earthquakes are sensitive to the rheological structure of the subseismogenic crust, and, in particular, the viscosity of the middle and lower crust. Here we develop a semianalytic solution for time-dependent interseismic velocities resulting from viscoelastic stress relaxation in a localized midcrustal layer in response to forcing by a sequence of periodic earthquakes. Earthquake cycle models with a weak midcrustal layer exhibit substantially more near-fault preseismic strain localization than do classic two-layer models at short (<100 yr) Maxwell times. We apply both this three-layer model and the classic two-layer model to geodetic observations before and after the 1997 (M_W = 7.6) Manyi and 2001 (M_W = 7.8) Kokoxili strike-slip earthquakes in Tibet to estimate the viscosity of the crust below a 20 km thick seismogenic layer. For these events, interseismic stress relaxation in a weak ((viscosity \leq10^{18.5} Pa⋅s)) and thin (height ≤20 km) midcrustal layer explains observations of both preseismic near-fault strain localization and rapid (>50 mm/yr) postseismic velocities in the years following the coseismic ruptures. We suggest that earthquake cycle models with a localized midcrustal layer can simultaneously explain both preseismic and postseismic geodetic observations with a single Maxwell viscosity, while the classic two-layer model requires a rheology with multiple relaxation time scales.