[1] The underestimation of the size of recent megathrust earthquakes illustrates our limited understanding of their spatiotemporal occurrence and governing physics. To unravel their relation to associated subduction dynamics and long-term deformation, we developed a 2-D continuum viscoelastoplastic model that uses an Eulerian-Lagrangian finite difference framework with similar on-and off-fault physics. We extend the validation of this numerical tool to a realistic subduction zone setting that resembles Southern Chile. The resulting quasi-periodic pattern of quasi-characteristic M8-M9 megathrust events compares quantitatively with observed recurrence and earthquake source parameters, albeit at very slow coseismic speeds. Without any data fitting, surface displacements agree with GPS data recorded before and during the 2010 M8.8 Maule earthquake, including the presence of a second-order flexural bulge. These surface displacements show cycle-to-cycle variations of slip deficits, which overall accommodate 5% of permanent internal shortening. We find that thermally (and stress) driven creep governs a spontaneous conditionally stable downdip transition zone between temperatures of 350 ı C and 450 ı C. Ruptures initiate above it (and below the forearc Moho), propagate within it, interspersed by small intermittent events, and arrest below it as ductile shearing relaxes stresses. Ruptures typically propagate upward along lithological boundaries and widen as pressures drop. The main thrust is constrained to be weak due to fluid-induced weakening required to sustain regular subduction and to generate events with natural characteristics (fluid pressures of 75-99% of solid pressures). The agreement with a range of seismological, geodetic, and geological observations demonstrates the validity and strength of this physically consistent seismo-thermo-mechanical approach.

The September 2018, M w 7.5 Sulawesi earthquake occurring on the Palu-Koro strike-slip fault system was followed by an unexpected localized tsunami. We show that direct earthquakeinduced uplift and subsidence could have sourced the observed tsunami within Palu Bay. To this end, we use a physics-based, coupled earthquake-tsunami modeling framework tightly constrained by observations. The model combines rupture dynamics, seismic wave propagation, tsunami propagation and inundation. The earthquake scenario, featuring sustained supershear rupture propagation, matches key observed earthquake characteristics, including the moment magnitude, rupture duration, fault plane solution, teleseismic waveforms and inferred horizontal ground displacements. The remote stress regime reflecting regional transtension applied in the model produces a combination of up to 6 m left-lateral slip and up to 2 m normal slip on the straight fault segment dipping 65 East beneath Palu Bay. The time-dependent, 3D seafloor displacements are translated into bathymetry perturbations with a mean vertical offset of 1.5 m across the submarine fault segment. This sources a tsunami with wave amplitudes and periods that match those measured at the Pantoloan wave gauge and inundation that reproduces observations from field surveys. We conclude that a source related to earthquake displacements is probable and that landsliding may not have been the primary source of the tsunami. These results have important implications for submarine strike-slip fault systems worldwide. Physics-based modeling offers rapid response specifically in tectonic settings that are currently underrepresented in operational tsunami hazard assessment.

We present a 2‐D numerical modeling approach for simulating a wide slip spectrum in a viscoelastoplastic continuum. The key new model component is an invariant reformulation of the classical rate‐ and state‐dependent friction equations, which is designed for earthquake simulations along spontaneously evolving faults. Here we describe the methodology and demonstrate that it is accurate and stable in a setup consisting of a mature strike‐slip fault zone. We show that the nucleation and propagation of an earthquake are well resolved, as supported by a good agreement with various analytical approximations, including those of the nucleation and cohesive zone lengths. Results generally converge with respect to grid size, time step, and other numerical parameters. The convergence rate with respect to grid size depends on the internodal averaging scheme, is influenced by wave reflections, and deteriorates for inclined faults. The simulated slip spectrum, ranging from stable sliding at the loading rate to periodic aseismic slip to periodic seismic slip as a function of nucleation size, is in general agreement with the literature. In this simple setup, dynamic pressure does not play a significant role. By analyzing the role of viscous deformation, we identify and confirm by our simulations a theoretical viscosity threshold below which earthquakes cannot nucleate. This threshold is shown to depend on the reference strength of rate‐ and state‐dependent friction and the loading strain rate, which is in agreement with previous work on the brittle‐ductile transition.

[1] The physics governing the seismic cycle at seismically active subduction zones remains poorly understood due to restricted direct observations in time and space. To investigate subduction zone dynamics and associated interplate seismicity, we validate a continuum, visco-elasto-plastic numerical model with a new laboratory approach (Paper 1). The analogous laboratory setup includes a visco-elastic gelatin wedge underthrusted by a rigid plate with defined velocity-weakening and -strengthening regions. Our geodynamic simulation approach includes velocity-weakening friction to spontaneously generate a series of fast frictional instabilities that correspond to analog earthquakes. A match between numerical and laboratory source parameters is obtained when velocity-strengthening is applied in the aseismic regions to stabilize the rupture. Spontaneous evolution of absolute stresses leads to nucleation by coalescence of neighboring patches, mainly occurring at evolving asperities near the seismogenic zone limits. Consequently, a crack-, or occasionally even pulse-like, rupture propagates toward the opposite side of the seismogenic zone by increasing stresses ahead of its rupture front, until it arrests on a barrier. The resulting surface displacements qualitatively agree with geodetic observations and show landward and, from near the downdip limit, upward interseismic motions. These are rebound and reversed coseismically. This slip increases adjacent stresses, which are relaxed postseismically by afterslip and thereby produce persistent seaward motions. The wide range of observed physical phenomena, including back-propagation and repeated slip, and the agreement with laboratory results demonstrate that visco-elasto-plastic geodynamic models with rate-dependent friction form a new tool that can greatly contribute to our understanding of the seismic cycle at subduction zones.Citation: van Dinther, Y., T. V. Gerya, L. A. Dalguer, F. Corbi, F. Funiciello, and P. M. Mai (2012), The seismic cycle at subduction thrusts: 2. Dynamic implications of geodynamic simulations validated with laboratory models, J. Geophys.

[1] Subduction megathrust earthquakes occur at the interface between the subducting and overriding plates. These hazardous phenomena are only partially understood because of the absence of direct observations, the restriction of the instrumental seismic record to the past century, and the limited resolution/completeness of historical to geological archives. To overcome these restrictions, modeling has become a key-tool to study megathrust earthquakes. We present a novel model to investigate the seismic cycle at subduction thrusts using complementary analog (paper 1) and numerical (paper 2) approaches. Here we introduce a simple scaled gelatin-on-sandpaper setup including realistic tectonic loading, spontaneous rupture nucleation, and viscoelastic response of the lithosphere. Particle image velocimetry allows to derive model deformation and earthquake source parameters. Analog earthquakes are characterized by "quasi-periodic" recurrence. Consistent with elastic theory, the interseismic stage shows rearward motion, subsidence in the outer wedge and uplift of the "coastal area" as a response of locked plate interface at shallow depth. The coseismic stage exhibits order of magnitude higher velocities and reversal of the interseismic deformation pattern in the seaward direction, subsidence of the coastal area, and uplift in the outer wedge. Like natural earthquakes, analog earthquakes generally nucleate in the deeper portion of the rupture area and preferentially propagate upward in a crack-like fashion. Scaled rupture width-slip proportionality and seismic moment-duration scaling verifies dynamic similarities with earthquakes. Experimental repeatability is statistically verified. Comparing analog results with natural observations, we conclude that this technique is suitable for investigating the parameter space influencing the subduction interplate seismic cycle.