Identifying the Recurrence Patterns of Nonvolcanic Tremors Using a 2‐D Hidden Markov Model With Extra Zeros

Nonvolcanic tremor activity has been observed in many places worldwide. In some regions, their activity was observed to accompany slow slip events. Before examining whether and how nonvolcanic tremor activity is related to slow slip, it is essential to understand quantitatively the spatiotemporal migration patterns of nonvolcanic tremors. We developed a 2‐D hidden Markov model to automatically analyze and forecast the spatiotemporal behavior of tremor activity in the regions Kii and Shikoku, southwest Japan. This new automated procedure classifies the tremor source regions into distinct segments in 2‐D space and infers a clear hierarchical structure of tremor activity, where each region consists of several subsystems and each subsystem contains several segments. The segments can be quantitatively categorized into three different types according to their occurrence patterns: episodic, weak concentration, and background, extending earlier knowledge gained from handpicked tremor swarms. The Kii region can be categorized into four different subsystems, with two often linked to each other. The Shikoku region can be divided into six subsystems, with two in central Shikoku linked to each other. Moreover, a significant increase in the proportion of tremor occurrence was detected in a segment in southwest Shikoku before the 2003 and 2010 long‐term slow slip events in the Bungo channel. This highlights the possible correlation between nonvolcanic tremor and slow slip events. The model can be used to analyze tremor data from other regions.

[1]  Andrew J. Viterbi,et al.  Error bounds for convolutional codes and an asymptotically optimum decoding algorithm , 1967, IEEE Trans. Inf. Theory.

[2]  L. Baum,et al.  A Maximization Technique Occurring in the Statistical Analysis of Probabilistic Functions of Markov Chains , 1970 .

[3]  Jr. G. Forney,et al.  The viterbi algorithm , 1973 .

[4]  D. Rubin,et al.  Maximum likelihood from incomplete data via the EM - algorithm plus discussions on the paper , 1977 .

[5]  G. Schwarz.  Estimating the Dimension of a Model , 1978 .

[6]  Kazushige Obara,et al.  Nonvolcanic Deep Tremor Associated with Subduction in Southwest Japan , 2002, Science.

[7]  Andrea Donnellan,et al.  A Hidden Markov Model Based Tool for Geophysical Data Exploration , 2002 .

[8]  H. Dragert,et al.  Episodic Tremor and Slip on the Cascadia Subduction Zone: The Chatter of Silent Slip , 2003, Science.

[9]  H. Kao,et al.  The Source‐Scanning Algorithm: mapping the distribution of seismic sources in time and space , 2004 .

[10]  H. Hirose,et al.  Episodic slow slip events accompanied by non‐volcanic tremors in southwest Japan subduction zone , 2004 .

[11]  Makoto Murakami,et al.  Aseismic slip and low‐frequency earthquakes in the Bungo channel, southwestern Japan , 2004 .

[12]  Hitoshi Hirose,et al.  Repeating short- and long-term slow slip events with deep tremor activity around the Bungo channel region, southwest Japan , 2005 .

[13]  G. Beroza,et al.  Low-frequency earthquakes in Shikoku, Japan, and their relationship to episodic tremor and slip , 2006, Nature.

[14]  Richard M. Allen,et al.  Segmentation in episodic tremor and slip all along Cascadia , 2006 .

[15]  Michael R. Brudzinski,et al.  Determination of slow slip episodes and strain accumulation along the Cascadia margin , 2006 .

[16]  M. Bebbington.  Identifying volcanic regimes using Hidden Markov Models , 2007 .

[17]  Kazushige Obara,et al.  Phenomenology of deep slow earthquake family in southwest Japan , 2010 .

[18]  Takuto Maeda,et al.  Spatiotemporal distribution of seismic energy radiation from low‐frequency tremor in western Shikoku, Japan , 2009 .

[19]  Douglas H. Christensen,et al.  Possible relationship between nonvolcanic tremor and the 1998–2001 slow slip event, south central Alaska , 2009 .

[20]  T. Maeda,et al.  Slow Earthquakes Linked Along Dip in the Nankai Subduction Zone , 2010, Science.

[21]  H. Hirose,et al.  Recurrence behavior of short‐term slow slip and correlated nonvolcanic tremor episodes in western Shikoku, southwest Japan , 2010 .

[22]  Gerassimos A. Papadopoulos,et al.  Identifying Seismicity Levels via Poisson Hidden Markov Models , 2010 .

[23]  S. Sekine,et al.  Along‐strike variations in short‐term slow slip events in the southwest Japan subduction zone , 2010 .

[24]  S. Ide.  Striations, duration, migration and tidal response in deep tremor , 2010, Nature.

[25]  Sachiko Tanaka,et al.  Depth‐dependent activity of non‐volcanic tremor in southwest Japan , 2010 .

[26]  Migration properties of non‐volcanic tremor in Shikoku, southwest Japan , 2011 .

[27]  K. Creager,et al.  A continuum of stress, strength and slip in the Cascadia subduction zone , 2011 .

[28]  Kazushige Obara,et al.  Characteristics and interactions between non-volcanic tremor and related slow earthquakes in the Nankai subduction zone, southwest Japan , 2011 .

[29]  T. Maeda,et al.  Depth‐dependent mode of tremor migration beneath Kii Peninsula, Nankai subduction zone , 2012 .

[30]  Ting Wang,et al.  Markov-modulated Hawkes process with stepwise decay , 2012 .

[31]  Mark Bebbington,et al.  Identifying anomalous signals in GPS data using HMMs: An increased likelihood of earthquakes? , 2013, Comput. Stat. Data Anal..

[32]  Richard M. Allen,et al.  Automated detection and location of tectonic tremor along the entire Cascadia margin from 2005 to 2011 , 2015 .

[33]  Total energy of deep low‐frequency tremor in the Nankai subduction zone, southwest Japan , 2016 .

[34]  A. Kato,et al.  Connecting slow earthquakes to huge earthquakes , 2016, Science.

[35]  Ting Wang,et al.  Hidden Markov modelling of sparse time series from non‐volcanic tremor observations , 2017 .