The Influence of Depth-Varying Elastic Properties of the Upper Plate on Megathrust Earthquake Rupture Dynamics and Tsunamigenesis
18 pages, 10 figures, supporting information https://doi.org/10.1029/2021JB022328.-- Data Availability Statement: The results of the simulations through each fault domain, scripts needed to build the unstructured mesh, and the files containing the frictional setup can be found here: https://figshare.com/projects/Dynamic_rupture_simulation_results_from_Prada_et_al_preprint_https_doi_org_10_1002_essoar_10507023_1_/120540.-- This is a contribution of the Barcelona Center for Subsurface Imaging that is a Grup de Recerca de la Generalitat de Catalunya (2017 SGR 1662) ; Megathrust earthquakes are strongly influenced by the elastic properties of rocks surrounding the fault. In contrast to friction, these properties can be derived in situ from geophysical measurements along the seismogenic zone. However, they are often overestimated in numerical simulations, particularly in the shallow megathrust. Here we explore the influence that realistic depth-varying upper-plate elastic properties along the megathrust have on earthquake rupture dynamics and tsunamigenesis using 3D dynamic rupture and tsunami simulations. We compare results from three subduction zone scenarios with homogeneous and heterogeneous elastic media, and bimaterial fault. Elastic properties in the heterogeneous model follow a realistic depth-distribution derived from controlled-source tomography models of subduction zones. We assume the same friction properties for all scenarios. Simulations in the heterogeneous and homogeneous models show that rigidity depth variations explain the depth-varying behavior of slip, slip rate, frequency content, and rupture time. The depth-varying behavior of slip, frequency content, and rupture duration quantitatively agree with previous predictions based on worldwide data compilations, explaining the main depth-dependent traits of tsunami earthquakes and large shallow megathrust earthquakes. Large slip, slow rupture and slip rate amplification in bimaterial simulations are largely controlled by elastic rock properties of the most compliant side of the fault, which in subduction zones is the upper plate. Large shallow slip and trenchward increasing upper-plate compliance of the heterogeneous model lead to the largest co-seismic seafloor deformation and tsunami amplitude. This highlights the importance of considering realistic upper-plate rigidity variations to properly assess the tsunamigenic potential of megathrust earthquakes ; M. Prada is funded by the Beatriu de Pinós postdoctoral programme of the Government of Catalonia's Secretariat for Universities and Research of the Ministry of Economy and Knowledge (Ref # 2017BP00170). Jean-Paul Ampuero is supported by the Investments in the Future project UCAJEDI (ANR-15-IDEX-01) managed by the French National Research Agency (ANR). First author and ICM have had funding support of the "Severo Ochoa Centre of Excellence" accreditation (CEX2019-000928-S) of the Spanish Research Agency (AEI). [.] Tsunami-HySEA code development is supported by the Spanish Government-FEDER funded project MEGAFLOW (RTI2018-096064-B-C21) ; Peer reviewed