Coupled afterslip and viscoelastic flow following the 2002 Denali Fault, Alaska earthquake

Abstract

SUMMARY Models of postseismic deformation following the 2002 M 7.9 Denali Fault, Alaska earthquake provide insight into the rheologic structure of the Alaskan lithosphere and the physical processes activated following a large earthquake. We model coseismic GPS displacements and 4 yr of postseismic GPS position time‐series with a coupled model of afterslip on the fault in the lithosphere and distributed viscous flow in the asthenosphere. Afterslip is assumed to be governed by a simplified version of a laboratory‐derived rate‐strengthening friction law that is characterized with a single parameter, σ(a−b), where σ is the effective normal stress on the fault and a−b is a dimensionless empirical parameter. Afterslip is driven by coseismic stress changes on the fault generated by the main shock. The lithosphere is modelled as an elastic plate overlying a linear, Maxwell, viscoelastic asthenosphere. We devise a scheme to simultaneously estimate the distributions of coseismic slip and afterslip, friction parameters, lithosphere thickness and asthenosphere viscosity. The postseismic GPS time‐series are best reproduced with a 45–85 km thick elastic lithosphere overlying an asthenosphere of viscosity 0.6–1.5 × 1019 Pa s. The 45–85 km elastic lithosphere thickness is greater than or equal to the average crustal thickness in the region suggesting that distributed postseismic flow occurs largely within the mantle while postseismic deformation in the crust is confined largely to slip on a discrete fault penetrating the entire crust and perhaps into the uppermost mantle. For typical laboratory values of a−b of order 10−2 at temperatures corresponding to the inferred depth of afterslip, the estimated effective normal stress on the fault is ∼50 MPa, which is about an order of magnitude lower than effective normal stresses at hydrostatic pore pressure. Previous studies showed that models with linear Newtonian rheology cannot reproduce the observed GPS time‐series but that models incorporating non‐linear or biviscous flow of the mantle do fit the data. We show that a model with afterslip governed by a non‐linear fault zone rheology coupled to Newtonian mantle flow is sufficient to reproduce the GPS time‐series.

Glacial Isostatic Adjustment modelling: historical perspectives, recent advances, and future directions

“…The low viscosity values inferred from GIA for regions such as Iceland and Alaska are similar to the values determined in 815 studies of post-seismic deformation (Arnadottir et al, 2005;Johnson et al, 2009), and it has been suggested that power-law Earth Surf. Dynam.…”

Using a spatially realistic load model to assess impacts of Alaskan glacier ice loss on sea level

“…It does not include tectonic effects. It is therefore possible that mascons 4–6 and 8–9 may be affected by postseismic deformation following the M w 9.2 1964 Alaska earthquake [ Zweck et al , 2002; Suito and Freymueller , 2009], while estimates for mascons 2 and 3 may have been affected by the 2002 M 7.9 Denali Fault earthquake [ Freed et al , 2006; Johnson et al , 2009]. …”

Simulations of tremor‐related creep reveal a weak crustal root of the San Andreas Fault

“…The effective normal stress on the fault is low, not exceeding 10 MPa for the range of parameters explored here. Assuming typical laboratory values of order a = 0.01, s(a À b) is at least 1-2 orders of magnitude lower than values of 0.1-0.7 MPa inferred from numerous afterslip studies [e.g., Hearn et al, 2002;Perfettini and Avouac, 2007;Johnson et al, 2008]. Unless effective normal stress is exceedingly small, a-b must be small, of order 10 À 6 to 10 À 3 (Figure 2d).…”