Quasi-Dynamic 3D Modeling of the Generation and Afterslip of a Tohoku-oki Earthquake Considering Thermal Pressurization and Frictional Properties of the Shallow Plate Boundary

The generation of the 2011 Tohoku-oki earthquake has been modeled by many authors by considering a dynamic weakening mechanism such as thermal pressurization (TP). Because the effects of TP on afterslip have not been investigated, this study develops a 3D quasi-dynamic model of the earthquake cycle to investigate afterslip of the Tohoku-oki earthquake, considering TP and the geometry of the plate boundary. We employ several velocity-weakening (VW) patches for M-w 7 class events, and two large shallow VW patches. The frictional properties are set as velocity-strengthening (VS) outside the VW patches. The results show that, during megathrust earthquakes, fast slip propagates to the surrounding VS regions near the VW patches owing to weakening by TP. Following M-w 9 events, large afterslips occur in regions below the northern shallow rupture area in the off-Fukushima region close to the Japan Trench, which is consistent with observations. In the VS region near the VW patches, during the early afterslip period, frictional behavior exhibits less VS with increasing slip velocity due to pore pressure reduction. We also consider the frictional properties of the shallow plate boundary fault off Tohoku, which exhibits a transition from VW to VS from low to high slip velocities. The results show the occurrence of slow slip events (SSEs) at intervals of a few decades at the shallow plate boundary. During megathrust events, the VW property at low slip velocity promotes slip along the shallow SSE region more than the case with VS property throughout the entire velocity range.