Unique Pre-Earthquake Deformation Patterns in the Spatial Domains from GPS in Taiwan

Most earthquakes are considered to be caused by stress accumulating in and subsequently releasing from the crust. To extract non-linear and non-stationary earthquake-induced signals associated with stress accumulation, the Hilbert–Huang transform was utilized to filter long-term movements, short-term noise, and frequency-dependent (annual and semi-annual) variations from surface displacements measured by the global positioning system (GPS) in Taiwan. Earthquakerelated surface displacements were expressed as horizontal directions (i.e., GPS azimuths) using the north–south and east–west components of residual GPS data to bypass influences resulted from the inhomogeneous nature of the crust. Analytical results showed that the relationships between earthquake occurrence and the aligned GPS azimuth passed the statistical test of the Molchan’s error diagram. Aligned GPS azimuths were in agreement with direction of earthquake-related P axes for 81% (26/32) studied events. Areas with the highest paralleling orientations of GPS azimuths appeared around epicenters several days to weeks before earthquake occurrence. Durations from aligned GPS azimuths to earthquake occurrence are roughly proportional to earthquake magnitude. Similar variations of the GPS azimuths were observed in GPS data containing or excluding coseismic dislocation (i.e., one day before) in the temporal and spatial domain. These suggest that the aligned GPS azimuth could be a promising anomalous phenomenon for studying crustal deformation before earthquakes.

[1]  Mark A. Meadows,et al.  Shear-wave polarizations and subsurface stress directions at Lost Hills field , 1991 .

[2]  Peng Han,et al.  Evaluation of ULF seismo-magnetic phenomena in Kakioka, Japan by using Molchan's error diagram , 2017 .

[3]  J. Saastamoinen,et al.  Contributions to the theory of atmospheric refraction , 1972 .

[4]  Katsumi Hattori,et al.  Surface displacements in Japan before the 11 March 2011 M9.0 Tohoku-Oki earthquake , 2014 .

[5]  G. Molchan,et al.  Structure of optimal strategies in earthquake prediction , 1991 .

[6]  Demitris Paradissis,et al.  GPS constraints on continental deformation in the Africa‐Arabia‐Eurasia continental collision zone and implications for the dynamics of plate interactions , 2005 .

[7]  R. Robinson,et al.  The enigma of the Arthur's Pass, New Zealand, earthquake: 2. The aftershock distribution and its relation to regional and induced stress fields , 2000 .

[8]  R. Bürgmann,et al.  Coseismic slip distribution of the 2003 Mw 6.6 San Simeon earthquake, California, determined from GPS measurements and seismic waveform data , 2006 .

[9]  S. S. Shen,et al.  A confidence limit for the empirical mode decomposition and Hilbert spectral analysis , 2003, Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[10]  J. Chéry,et al.  Active surface deformation and sub-lithospheric processes in the western Mediterranean constrained by numerical models , 2010 .

[11]  Jing-Shan Hong,et al.  Determining the precipitable water vapor thresholds under different rainfall strengths in Taiwan , 2017 .

[12]  Stuart Crampin,et al.  The fracture criticality of crustal rocks , 1994 .

[13]  Guojie Meng,et al.  Groundwater–strain coupling before the 1999 Mw 7.6 Taiwan Chi-Chi earthquake , 2015 .

[14]  K. Feigl,et al.  Coseismic and Postseismic Fault Slip for the 17 August 1999, M = 7.5, Izmit, Turkey Earthquake. , 2000, Science.

[15]  Jer-Ming Chiu,et al.  Temporal changes in shear wave splitting during an earthquake swarm in Arkansas , 1990 .

[16]  J. Chéry.  Geodetic strain across the San Andreas fault reflects elastic plate thickness variations (rather than fault slip rate) , 2008 .

[17]  Shu-Hao Chang,et al.  Surface Deformation and Seismic Rebound: Implications and Applications , 2011 .

[18]  Yehuda Bock,et al.  Instantaneous global plate motion model from 12 years of continuous GPS observations , 2004 .

[19]  N. Huang,et al.  The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis , 1998, Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[20]  V. Gahalaut,et al.  Constraints on 2004 Sumatra Andaman earthquake rupture from GPS measurements in Andaman Nicobar Islands , 2006 .

[21]  P. Vernant,et al.  GPS evidence for northward motion of the Sinai Block: Implications for E. Mediterranean tectonics , 2004 .

[22]  Theodora Volti,et al.  A successfully stress-forecast earthquake , 1999 .

[23]  K. W. Cattermole.  The Fourier Transform and its Applications , 1965 .

[24]  T. Yeh,et al.  Ground-based GPS remote sensing for precipitable water vapor: A case study of the heat-island effect in Taipei , 2019, Terrestrial, Atmospheric and Oceanic Sciences.

[25]  Guochang Xu,et al.  GPS Height and Gravity Variations Due to Ocean Tidal Loading Around Taiwan , 2008 .

[26]  J. Ray,et al.  Anomalous harmonics in the spectra of GPS position estimates , 2008 .

[27]  Y. Hobara,et al.  Observation of surface displacements from GPS analyses before and after the Jiashian earthquake (M = 6.4) in Taiwan , 2013 .

[28]  Geoffrey Blewitt,et al.  Effect of annual signals on geodetic velocity , 2002 .

[29]  M. Farge.  Wavelet Transforms and their Applications to Turbulence , 1992 .

[30]  Tzay-Chyn Shin,et al.  NOTES AND CORRESPONDENCE Continuous CWB GPS Array in Taiwan and Applications to Monitoring Seismic Activity , 2011 .

[31]  Frederic Masson,et al.  Present‐day crustal deformation and plate kinematics in the Middle East constrained by GPS measurements in Iran and northern Oman , 2004 .

[32]  Mustafa Aktar,et al.  Space and Time Evolution of Rupture and Faulting during the 1999 İzmit (Turkey) Earthquake , 2002 .

[33]  D. Winterstein,et al.  Shear-wave polarizations and subsurface stress directions at Lost Hills field , 1991 .

[34]  I. P. Dobrovolsky,et al.  Estimation of the size of earthquake preparation zones , 1979 .

[35]  Jessica R. Murray,et al.  Coseismic and initial postseismic deformation from the 2004 Parkfield, California, earthquake, observed by global positioning system, electronic distance meter, creepmeters, and borehole strainmeters , 2006 .

[36]  Paul Lundgren,et al.  Geodetic and Seismic Constraints on some Seismogenic Zone Processes in Costa Rica , 2003 .

[37]  T. Dixon,et al.  17. Elastic and Viscoelastic Models of Crustal Deformation in Subduction Earthquake Cycles , 2007 .

[38]  Geoffrey Blewitt,et al.  Crustal displacements due to continental water loading , 2001 .

[39]  Simon McClusky,et al.  GPS constraints on Africa (Nubia) and Arabia plate motions , 2003 .

[40]  J. Normandeau,et al.  Tectonic implications of a dense continuous GPS velocity field at Yucca Mountain, Nevada , 2004 .

[41]  Chuan-Sheng Wang,et al.  Determination of global positioning system (GPS) receiver clock errors: impact on positioning accuracy , 2009 .

[42]  R. Rau,et al.  Anomalous frequency characteristics of groundwater level before major earthquakes in Taiwan , 2012 .

[43]  Norden E. Huang,et al.  A review on Hilbert‐Huang transform: Method and its applications to geophysical studies , 2008 .

[44]  Harry Fielding Reid,et al.  The mechanics of the earthquake , 1910 .

[45]  Freysteinn Sigmundsson,et al.  Current plate movements across the Mid-atlantic ridge determined from 5 years of continuous GPS measurements in Iceland , 2006 .