Abstract:Understanding of postseismic deformation following great subduction zone earthquakes is complicated by the combined effects of viscoelastic relaxation of earthquake-induced stresses in the upper mantle and time-dependent afterslip on the megathrust. We integrate geodetic observations and constraints on afterslip from small repeating earthquakes on the megathrust to better distinguish contributions from these two postseismic processes. We have developed a three-dimensional, spherical viscoelastic finite element model to study the postseismic deformation of the 2011 M(w)9.0 Tohoku earthquake that has been recorded at unprecedented high resolution in space and time. We model stress-driven afterslip in a 2 km thick weak shear zone away from historic rupture zones on the megathrust. We model both the viscoelastic relaxation of the upper mantle and shear zone deformation with a transient Burgers body rheology. The transient Kelvin viscosity is assumed to be one order of magnitude lower than that of the Maxwell viscosity. Viscoelastic relaxation of the mantle wedge alone causes postseismic uplift and seaward motion in the upper plate, opposite to the pattern from relaxation of just the oceanic upper mantle. Afterslip on the fault produces uplift updip of the afterslip zone and subsidence over its downdip edge and mostly seaward motion above the afterslip zone. The best fit Maxwell viscosity of the shear zone at depths <= 50km is 10(17) Pa s, constrained by afterslip estimates from repeating earthquakes. The optimal viscosities of the deep weak shear zone, continental mantle wedge, and oceanic upper mantle are determined to be 5 x 10(17) Pa s, 3 x 10(19) Pa s, and 5 x 10(19) Pa s, respectively. The stress-driven afterslip in the shear zone is up to similar to 3.5 min the first 2 years after the earthquake, equivalent to an M(w)8.4. Our model reproduces the first-order pattern of the GPS observations both in horizontal and in vertical directions. Seafloor geodetic observations of subsidence and landward motions near the high-slip zone of the earthquake provide evidence for a low-viscosity asthenosphere below the oceanic lithosphere.

Weak Anatolian upper mantle inferred from postseismic deformation of the 2023 T&#252;rkiye earthquake doublet

Far‐Field Viscoelastic Post‐Seismic Deformation of the 2011 Mw9 Tohoku Earthquake

Large-scale Variation in Seismic Anisotropy in the Crust and Upper Mantle Beneath Anatolia, Turkey

Stress‐driven Relaxation of Heterogeneous Upper Mantle and Time‐dependent Afterslip Following the 2011 Tohoku Earthquake

Immature characteristics of the East Anatolian Fault Zone from SAR, GNSS and strong motion data of the 2023 Türkiye–Syria earthquake doublet

Isotropic and Anisotropic P Wave Velocity Structures of the Crust and Uppermost Mantle Beneath Turkey

Rheological Structure and Lithospheric Stress Interaction in the Alaska Subduction Zone Gleaned From the 2018 Mw 7.9 Oceanic Crustal Earthquake

Viscoelastic relaxation of the upper mantle and afterslip following the 2014 MW8.1 Iquique earthquake

Crust and upper mantle S wave velocity structure in eastern Turkey based on ambient noise tomography

Decadal Viscoelastic Postseismic Deformation of the 1964 Mw9.2 Alaska Earthquake

Early Postseismic Deformation of the 2010 Mw 6.9 Yushu Earthquake and Its Implication for Lithospheric Rheological Properties

Dominance of transient rheology in far-field postseismic deformation following the 2012 Mw 7.8 Haida Gwaii and 2013 Mw 7.5 Craig earthquakes

Crustal and Uppermost Mantle Isotropic and Anisotropic P-wave Velocity Variations Beneath Turkey&#160;

Coseismic slip distribution of the 2023 earthquake doublet in Turkey and Syria from joint inversion of Sentinel-1 and Sentinel-2 data: An iterative modeling method for mapping large earthquake deformation

Weakness of the Indian Lower Crust Beneath the Himalaya Inferred from Postseismic Deformation of the 2015 Mw 7.8 Gorkha Earthquake

Lithospheric rheology in West Tibet and Pamir constrained from the postseismic deformation of the 2005 Mw7.6 Kashmir (Pakistan) earthquake

High‐Resolution Broadband Lg Attenuation Structure of the Anatolian Crust and Its Implications for Mantle Upwelling and Plateau Uplift

The 2023 M w 7.7 Pazarcik earthquake caused a reversal in vertical motion along the SW branch of the East Anatolian fault

Geodynamics of East Anatolia‐Caucasus Domain: Inferences from 3D Thermo‐mechanical Models, Residual Topography, and Admittance Function Analyses

Integrated Investigation on Heterogeneous Lower Crust Rheology in Kyushu and Afterslip Behavior Following the 2016 Mw7.1 Kumamoto Earthquake

Weak Anatolian upper mantle inferred from postseismic deformation of the 2023 T&amp;#252;rkiye earthquake doublet