Cascading Rupture of Patches of High Seismic Energy Release Controls The
Total Page:16
File Type:pdf, Size:1020Kb
Nakamoto et al. Earth, Planets and Space (2021) 73:59 https://doi.org/10.1186/s40623-021-01384-6 FULL PAPER Open Access Cascading rupture of patches of high seismic energy release controls the growth process of episodic tremor and slip events Keita Nakamoto1, Yoshihiro Hiramatsu2* , Takahiko Uchide3 and Kazutoshi Imanishi3 Abstract Slip phenomena on plate interfaces refect the heterogeneous physical properties of the slip plane and, thus, exhibit a wide variety of slip velocities and rupture propagation behaviors. Recent fndings on slow earthquakes reveal simi- larities and diferences between slow and regular earthquakes. Episodic tremor and slip (ETS) events, a type of slow earthquake widely observed in subduction zones, likewise show diverse activity. We investigated the growth of 17 ETS events beneath the Kii Peninsula in the Nankai subduction zone, Japan. Analyses of waveform data recorded by a seismic array enabled us to locate tremor hypocenters and estimate the migration patterns and spatial distribu- tion of the energy release of tremor events. Here, we describe three major features in the growth of ETS events. First, independent of their start point and migration pattern, ETS events exhibit patches of high seismic energy release on the up-dip part of the ETS zone, suggesting that the location of these patches is controlled by inherent physical or frictional properties of the plate interface. Second, ETS events usually start outside the high-energy patches, and their fnal extent depends on whether the patches participate in the rupture. Third, we recognize no size dependence in the initiation phase of ETS events of diferent sizes with comparable start points. These features demonstrate that the cascading rupture of high-energy patches governs the growth of ETS events, just as the cascading rupture of asperi- ties governs the growth of regular earthquakes. Keywords: Slow earthquakes, Tectonic tremor, Subduction dynamics, Rapid tremor reversal Introduction Recent fndings about slow earthquakes have broad- Earthquakes feature heterogeneous slip distributions on ened our understanding of subduction dynamics and fault planes that are characterized by one to several areas the generation of megathrust earthquakes (e.g., Obara of large slip, or asperities. In general, earthquake hypo- and Kato 2016). In subduction zones, strain on the plate centers are located outside of asperities (e.g., Yamanaka interface is released by slow slips in slow earthquakes and and Kikuchi 2004; Mai et al. 2005). Te complex growth high-speed ruptures in regular earthquakes. Interestingly, and rupture processes of earthquakes are fundamental the spatial distributions of slow and regular earthquakes features that have fueled a long-running debate on how are complementary: slow earthquakes are distributed and when the eventual size of an earthquake is deter- around and outside the locked zones that are the source mined, specifcally the existence of a nucleation phase areas of megathrust earthquakes, both along dip (Obara (e.g., Ellsworth and Beroza 1995; Ide 2019). and Kato 2016) and along strike (Nishikawa et al. 2019). *Correspondence: [email protected] Furthermore, a temporal relationship between slow and 2 School of Geosciences and Civil Engineering, College of Science regular earthquakes has been proposed. Kato et al. (2012) and Engineering, Kanazawa University, Kakuma, Kanazawa 920-1192, reported the occurrence of two slow slip events (SSEs) Japan Full list of author information is available at the end of the article from earthquake activities that migrated toward the future hypocenter of the 2011 Tohoku earthquake before © The Author(s) 2021, corrected publication 2021. This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/. Nakamoto et al. Earth, Planets and Space (2021) 73:59 Page 2 of 17 the mainshock. Te migration speed of those hypocent- Te moment rate of slow slips is often estimated from ers, 2–5 km/day, is close to the typical migration speed the number, amplitude, or radiated energy of tremors of episodic tremor and slip (ETS) events, 10 km/day or LFEs. Ide et al. (2008) demonstrated that the radi- (Dragert et al. 2001; Obara 2002). Slow slips preceding a ated seismic energy rate of tremors in the 2–8 Hz band dynamic rupture have also been reported in laboratory is proportional to the seismic moment rate of very-low- experiments and numerical simulations (e.g., Dieterich frequency earthquakes in the Kii Peninsula. Ide and Yabe 1979; Ohnaka 1992; Shibazaki and Matsu’ura 1992). Mat- (2014) confrmed this proportionality between tremor suzawa et al. (2010) showed that in a numerical simula- energy rate and seismic moment rate of very-low-fre- tion, the recurrence interval of SSEs grew shorter in the quency earthquakes in the entire Nankai subduction time interval preceding a megathrust earthquake. Tere- zone. Hirose and Obara (2010) showed that the moment fore, investigating the features of slow earthquakes yields rate function of an SSE estimated from tiltmeter records clues to the generation of megathrust earthquakes in sub- was strongly correlated with the temporal change in the duction zones. number of tremors in Shikoku. Maeda and Obara (2009) In some subduction zones, slow earthquakes occur as found a clear correspondence between the temporal short-term SSEs in the form of ETS events, in which tec- variation in tilt records and the accumulated seismic tonic (nonvolcanic) tremor and low-frequency earthquake energy from tremor sources for ETS events in Shikoku. (LFE) are synchronized with slip events in time and space Hawthorne and Rubin (2013) found that the strain rate (Rogers and Dragert 2003). Te slip areas and slip rates of is proportional to the amplitude of tremor on time scales ETS events are strongly correlated with tremor epicent- shorter than one day. Tey interpreted this relation to ers (Hirose and Obara 2010; Bartlow et al. 2011; Michel mean that the moment rate of slow slip is correlated with et al. 2019a). Tus, tremor and LFE can be used to moni- tremor amplitude. tor ETS events (Shelly et al. 2007a; Obara 2010; Houston Seismic arrays are a powerful tool for detecting, locat- et al. 2011; Ghosh et al. 2012; Rubin and Armbruster ing, and characterizing tremor (Ghosh et al. 2009; Imani- 2013; Peng et al. 2015), although short-term SSEs are not shi et al. 2011). Te Geological Survey of Japan, National always accompanied by tremors (Obara et al. 2011; Wech Institute of Advanced Industrial Science and Technology and Bartlow 2014; Michel et al. 2019a). Tremor and LFE (GSJ), installed a seismic array of 39 sensors on the Kii amplitudes are also have been used to detect small epi- Peninsula from February 2011 to November 2016 (Fig. 1), sodes of transient aseismic slip hidden in the geodetic producing a set of well-recorded waveforms that has noise (Frank 2016; Frank et al. 2018). Te growth process been valuable for examining the detailed growth process of ETS events is not simple and is typifed by diverse pat- of ETS events from tectonic tremor. Sagae et al. (2021) terns of migration and eventual sizes (Obara et al. 2010), applied the MUSIC high-resolution frequency–wave- heterogeneous slip distributions (Hirose and Obara number (f–k) method (Schmidt 1986) to a dataset from 2010), and brief episodes of secondary slip front in which July 2012 to July 2014, revealing a more detailed view of tremor migrates faster than the main front: sometimes tremor migration than the conventional envelope cross- in the direction opposite to its overall migration (Hou- correlation method. We report here the spatio-temporal ston et al. 2011). Ghosh et al. (2012) reported the exist- distribution of the radiated seismic energy of tectonic ence of tremor asperities, areas of high radiated seismic tremor during 17 ETS events between April 2011 and energy, on the plate interface analogous to the asperities December 2014. We show that persistent patches of rela- of regular earthquakes and suggested that tremor asperi- tively high seismic energy release are located on the plate ties in the transition zone might control the evolution of interface in the ETS zone, and that the growth process of slow earthquakes in space and time. Savard and Bostock ETS events is controlled by cascading ruptures of these (2015) identifed the location of fve asperities from LFE patches. epicenters in Cascadia. Hirose and Obara (2010) con- ducted geodetic inversion and confrmed the existence of Data a persistent asperity of slow slip approximately 30 km in We analyzed waveform data recorded by the GSJ array diameter in Shikoku, Japan, that was active during several network in the Kii Peninsula from April 2011 to Decem- diferent ETS events. In the Nankai subduction zone of ber 2014 (Fig. 1). Te array consisted of 39 three-compo- southwestern Japan, the ETS zone can be divided into sev- nent velocity seismographs with a natural frequency of eral segments (e.g., Obara 2002) within which ETS events 2 Hz, and the sampling frequency of waveform data was have similar recurrence intervals, such as 3 or 6 months. 200 Hz.