Analysis of coseismic deformation using EOF method on dense, continuous GPS data in Taiwan

Abstract Taking advantage of the dense continuous GPS network data, we demonstrate the effective utility of the method of empirical orthogonal function (EOF) in obtaining the coseismic deformation (vector) field on the crustal surface caused by significant earthquakes. We conduct case studies w.r.t. two earthquakes that have occurred in Taiwan: one off the east coast of Taiwan in 2002 (Event I, Mw 7.1) and one inland in southern Taiwan in 2010 (Event II, Mw 5.7 or 6.3). We demonstrate that, superior to the simple least-squares estimates obtained directly from individual GPS position time series, the EOF analysis is capable of extracting the coseismic deformation in the form of coherent spatial pattern and time evolution. The prevalent common-mode errors in the GPS data are also well resolved as distinctive and separate EOF modes. The extracted coseismic spatial patterns provide evidences linking the regional tectonics with the orogenic process in Taiwan under the plate convergence. We further contrast our EOF results against those calculated via idealized elastic dislocation models, and find notable discrepancies especially for the far field, which can be attributed to the complexity of the tectonic environment of Taiwan.

[1]  Chia-Yu Lu,et al.  Oblique convergence, indentation and rotation tectonics in the Taiwan Mountain Belt: Insights from experimental modelling , 1994 .

[2]  Y. Bock,et al.  Anatomy of apparent seasonal variations from GPS‐derived site position time series , 2001 .

[3]  J. Angelier,et al.  Structural, geodetic and seismological evidence for tectonic escape in SW Taiwan , 2001 .

[4]  Robert W. King,et al.  Measurement of Crustal Deformation Using the Global Positioning System , 1991 .

[5]  J. Avouac,et al.  Inverting geodetic time series with a principal component analysis-based inversion method , 2010 .

[6]  Yehuda Bock,et al.  Error analysis of continuous GPS position time series , 2004 .

[7]  Thomas H. Heaton,et al.  The slip history of the 1994 Northridge, California, earthquake determined from strong-motion, teleseismic, GPS, and leveling data , 1996, Bulletin of the Seismological Society of America.

[8]  Y. Hsu,et al.  Source complexity of the 4 March 2010 Jiashian, Taiwan, Earthquake determined by joint inversion of teleseismic and near field data , 2013 .

[9]  Chin-Hsiung Loh,et al.  Reconnaissance Report on the 31 March 2002 Earthquake on the East Coast of Taiwan , 2003 .

[10]  J. Brian Gray,et al.  Leverage, residual, and interaction diagnostics for subsets of cases in least squares regression , 1997 .

[11]  Hiromichi Tsuji,et al.  Coseismic crustal deformation from the 1994 Hokkaido‐Toho‐Oki Earthquake Monitored by a nationwide continuous GPS array in Japan , 1995 .

[12]  Shui-Beih Yu,et al.  GPS measurement of postseismic deformation following the 1999 Chi‐Chi, Taiwan, earthquake , 2003 .

[13]  J. B. Rundle,et al.  Using Eigenpattern Analysis to Constrain Seasonal Signals in Southern California , 2004 .

[14]  H. Kao,et al.  Source Parameters of Regional Earthquakes in Taiwan: 1999-2000 Including the Chi-Chi Earthquake Sequence , 2002 .

[15]  Yehuda Bock,et al.  Detection of arbitrarily large dynamic ground motions with a dense high‐rate GPS network , 2004 .

[16]  Y. Okada Surface deformation due to shear and tensile faults in a half-space , 1985 .

[17]  R. Preisendorfer,et al.  Principal Component Analysis in Meteorology and Oceanography , 1988 .

[18]  William J. Emery,et al.  Data Analysis Methods in Physical Oceanography , 1998 .

[19]  C. S. Ho A synthesis of the geologic evolution of Taiwan , 1986 .

[20]  A. Deschamps,et al.  Arc-continent collision in Taiwan: New marine observations and tectonic evolution , 2002 .

[21]  J. Angelier,et al.  Shear concentration in a collision zone: kinematics of the Chihshang Fault as revealed by outcrop-scale quantification of active faulting, Longitudinal Valley, eastern Taiwan , 1997 .

[22]  Yehuda Bock,et al.  Spatiotemporal filtering using principal component analysis and Karhunen-Loeve expansion approaches for regional GPS network analysis , 2006 .

[23]  R. Rau,et al.  Contemporary deformation of tectonic escape in SW Taiwan from GPS observations, 1995–2005 , 2007 .

[24]  G. Petit,et al.  IERS Conventions (2010) , 2010 .

[25]  Rémi Michel,et al.  Horizontal coseismic deformation of the 1999 Chi-Chi earthquake measured from SPOT satellite images: Implications for the seismic cycle along the western foothills of central Taiwan , 2003 .

[26]  Co-seismic surface deformation of the 2011 off the Pacific coast of Tohoku Earthquake: Spatio-temporal EOF analysis of GPS data , 2011 .

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

[28]  Kenneth W. Hudnut,et al.  Detection of crustal deformation from the Landers earthquake sequence using continuous geodetic measurements , 1993, Nature.

[29]  A. R. H. Swan,et al.  Introduction to Geological Data Analysis , 1995 .

[30]  Timothy H. Dixon,et al.  Constraints on present‐day Basin and Range deformation from space geodesy , 1995 .

[31]  G. Grenerczy,et al.  Tectonic Processes in the Eurasian‐African Plate Boundary Zone Revealed by Space Geodesy , 2013 .

[32]  S. Yoo,et al.  Joint inversion of seismic and geodetic data for the source of the 2010 March 4, Mw 6.3 Jia-Shian, SW Taiwan, earthquake , 2013 .

[33]  K. Hudnut Geodesy tracks plate motion , 1992, Nature.

[34]  R. Rau,et al.  Inferred fault geometry and slip distribution of the 2010 Jiashian, Taiwan, earthquake is consistent with a thick-skinned deformation model , 2011 .

[35]  J. Angelier,et al.  Active deformation of Taiwan from GPS measurements and numerical simulations , 2001 .

[36]  C. Demets,et al.  Crustal velocity field of Mexico from continuous GPS measurements, 1993 to June 2001: Implications for the neotectonics of Mexico , 2003 .

[37]  S. Lallemand,et al.  AN INTRODUCTION TO ACTIVE COLLISION IN TAIWAN , 1997 .

[38]  John Langbein,et al.  Postearthquake relaxation after the 2004 M6 Parkfield, California, earthquake and rate‐and‐state friction , 2008 .

[39]  Wenke Sun,et al.  General formulations of global co-seismic deformations caused by an arbitrary dislocation in a spherically symmetric earth model-applicable to deformed earth surface and space-fixed point , 2009 .

[40]  Chung-Pai Chang,et al.  Strain and stress field in Taiwan oblique convergent system: constraints from GPS observation and tectonic data , 2003 .

[41]  W. Menke Geophysical data analysis : discrete inverse theory , 1984 .

[42]  Shui-Beih Yu,et al.  Velocity field of GPS stations in the Taiwan area , 1997 .

[43]  Yehuda Bock,et al.  Southern California permanent GPS geodetic array: Spatial filtering of daily positions for estimating coseismic and postseismic displacements induced by the 1992 Landers earthquake , 1997 .

[44]  Tsai Yi-Ben,et al.  Seismotectonics of Taiwan , 1986 .

[45]  Cheinway Hwang,et al.  Vertical Displacement due to Ocean Tidal Loading Around Taiwan Based on GPS Observations , 2011 .

[46]  Walter H. F. Smith,et al.  Gridding with continuous curvature splines in tension , 1990 .

[47]  T. Dixon,et al.  Noise in GPS coordinate time series , 1999 .

[48]  W. Hammond The Ghost of an Earthquake , 2005, Science.

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

[50]  Tao Yuan,et al.  Concurrent OH imager and sodium temperature//wind lidar observation of a mesopause region undular bore event over Fort Collins//Platteville, Colorado , 2004 .

[51]  Shui-Beih Yu,et al.  Coseismic Movement and Seismic Ground Motion Associated with the 31 March 2002 off Hualien, Taiwan, Earthquake , 2004 .

[52]  Michael B. Heflin,et al.  The effect of the second order GPS ionospheric correction on receiver positions , 2003 .