Crustal Deformation Prior to the 2017 Jiuzhaigou, Northeastern Tibetan Plateau (China), Ms 7.0 Earthquake Derived from GPS Observations

The 2017 Jiuzhaigou Ms 7.0 earthquake occurred on the northeastern margin of the Tibetan Plateau, with no noticeable rupture surface recognized. We characterized the pre-seismic deformation of the earthquake from GPS (Global Positioning System) data at eight continuous and 73 campaign sites acquired over the 2009–2017 period. With respect to the Eurasian plate, the velocity field showed a noticeable decrease, from west of the epicenter of the Jiuzhaigou earthquake to the western edge of the Longmenshan fault, in the southeast direction. The total northwest west–southeast east shortening rate in the vicinity of the epicentral area was in the range of 1.5 mm/y to 3.1 mm/y. With a GPS velocity transect across the Huya fault (HYF), where the epicenter was located, we estimated the activity of the HYF, showing a dominant left-lateral slip rate of 3.3 ± 0.2 mm/y. We calculated strain rates using a spherical wavelet-based multiscale approach that solved for the surface GPS velocity according to multiscale wavelet basis functions while accounting for spatially variable spacing of observations. Multiscale components of the two-dimensional strain rate tensor showed a complex crustal deformation pattern. Our estimates of strain rate components at the scale of seven and eight revealed extensional strain rate on the northern extension of the HYF. The Jiuzhaigou earthquake occurred at the buffer zone between extensional and compressional deformation, and with significant maximum shear rates being 100–140 nanostrain/y. In addition, a maximum shear strain rate of 60–120 nanostrain/y appeared around the epicenter of the 2013 Ms 6.6 Minxian–Zhangxian earthquake. These findings imply that inherent multiscale strain rates could be separated to identify strain accumulation related to mediumand large-sized earthquakes.

[1]  Z. Altamimi,et al.  ITRF2014: A new release of the International Terrestrial Reference Frame modeling nonlinear station motions , 2016 .

[2]  O. Francis,et al.  Modelling the global ocean tides: modern insights from FES2004 , 2006 .

[3]  Mark Simons,et al.  Interseismic crustal deformation in the Taiwan plate boundary zone revealed by GPS observations, seismicity, and earthquake focal mechanisms , 2008 .

[4]  Gene H. Golub,et al.  Tikhonov Regularization and Total Least Squares , 1999, SIAM J. Matrix Anal. Appl..

[5]  Z. Altamimi,et al.  ITRF2014 plate motion model , 2017 .

[6]  Leigh H. Royden,et al.  Topographic ooze: Building the eastern margin of Tibet by lower crustal flow , 2000 .

[7]  Mariusz Figurski.  Investigation of the impact of ITRF2014/IGS14 on the positions of the reference stations in Europe , 2017 .

[8]  S. Ward.  On the consistency of earthquake moment rates, geological fault data, and space geodetic strain: the United States , 1998 .

[9]  T. Herring,et al.  Introduction to GAMIT/GLOBK , 2006 .

[10]  Yan Zhan,et al.  The 2017 Jiuzhaigou Earthquake: A Complicated Event Occurred in a Young Fault System , 2018 .

[11]  Zhen Li,et al.  Effect of Surface Mass Loading on Geodetic GPS Observations , 2018 .

[12]  Simon D. P. Williams,et al.  Fast error analysis of continuous GNSS observations with missing data , 2013, Journal of Geodesy.

[13]  David D. Jackson,et al.  Crustal deformation across and beyond the Los Angeles basin from geodetic measurements , 1996 .

[14]  T. Nilsson,et al.  GPT2: Empirical slant delay model for radio space geodetic techniques , 2013, Geophysical research letters.

[15]  J. C. Savage,et al.  Strain accumulation and rotation in the Eastern California Shear Zone , 2001 .

[16]  Chieh‐Hung Chen,et al.  Heterogeneous strain regime in the eastern margin of Tibetan Plateau and its tectonic implications , 2015 .

[17]  Gene H. Golub,et al.  Generalized cross-validation as a method for choosing a good ridge parameter , 1979, Milestones in Matrix Computation.

[18]  I. Bogdanovaa,et al.  Stereographic wavelet frames on the sphere , 2005 .

[19]  David LaVallee,et al.  Higher‐order ionospheric effects on the GPS reference frame and velocities , 2009 .

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

[21]  Robert W. King,et al.  Fault locking, block rotation and crustal deformation in the Pacific Northwest , 2007 .

[22]  Xiaoping Yang,et al.  Basic characteristics of active tectonics of China , 2003, Science in China Series D Earth Sciences.

[23]  Jianping Wu,et al.  Relocation of mainshock and aftershock sequence of the M s 7.0 Sichuan Jiuzhaigou earthquake , 2018 .

[24]  Richard W. Allmendinger,et al.  Strain and rotation rate from GPS in Tibet, Anatolia, and the Altiplano , 2007 .

[25]  Roland Bürgmann,et al.  Contemporary crustal deformation around the southeast borderland of the Tibetan Plateau: TIBET SOUTHWEST BORDERLAND DEFORMATION , 2005 .

[26]  Wen Xue,et al.  Correlation of major earthquake sequences on the northern and eastern boundaries of the Bayan Har block, and its relation to the 2008 Wenchuan earthquake , 2011 .

[27]  Peter Molnar,et al.  Active faulting and tectonics in China , 1977 .

[28]  K. Whipple,et al.  Neotectonics of the Min Shan, China: Implications for mechanisms driving Quaternary deformation along the eastern margin of the Tibetan Plateau , 2000 .

[29]  Alain Geiger,et al.  GPS‐derived strain rate field within the boundary zones of the Eurasian, African, and Arabian Plates , 2000 .

[30]  Xi-wei Xu,et al.  Latest Quaternary paleoseismology and slip rates of the Longriba fault zone, eastern Tibet: Implications for fault behavior and strain partitioning , 2013 .

[31]  Xiaohong Zhang,et al.  The Effects of Higher-Order Ionospheric Terms on GPS Tropospheric Delay and Gradient Estimates , 2018, Remote. Sens..

[32]  Galina Dick,et al.  Systematic errors of mapping functions which are based on the VMF1 concept , 2014, GPS Solutions.

[33]  R. L. Hardy,et al.  Least squares prediction of gravity anomalies, geoidal undulations, and deflections of the vertical with multiquadric harmonic functions , 1975 .

[34]  Bertrand Meyer,et al.  Oblique Stepwise Rise and Growth of the Tibet Plateau , 2001, Science.

[35]  Brendan J. Meade,et al.  Present-day kinematics at the India-Asia collision zone , 2007 .

[36]  Yanqiang Wu,et al.  Comparison of GPS strain rate computing methods and their reliability , 2011 .

[37]  Chen Chang.  Division,deformation and tectonic implication of active blocks in the eastern segment of Bayan Har block , 2013 .

[38]  张会平,et al.  2017年8月8日九寨沟 M 7.0地震及余震震源机制解与发震构造分析 , 2017 .

[39]  Guojie Meng,et al.  Positioning Performance of BDS Observation of the Crustal Movement Observation Network of China and Its Potential Application on Crustal Deformation , 2018, Sensors.

[40]  John Haines,et al.  Contemporary horizontal velocity and strain rate fields of the Pacific‐Australian plate boundary zone through New Zealand , 2001 .

[41]  B. Burchfiel,et al.  The Geological Evolution of the Tibetan Plateau , 2008, Science.

[42]  P. Molnar,et al.  Cenozoic Tectonics of Asia: Effects of a Continental Collision: Features of recent continental tectonics in Asia can be interpreted as results of the India-Eurasia collision. , 1975, Science.

[43]  Michael A. Ellis,et al.  Active tectonics of the Beichuan and Pengguan faults at the eastern margin of the Tibetan Plateau , 2007 .

[44]  J. Ray,et al.  The IGS contribution to ITRF2014 , 2016, Journal of Geodesy.

[45]  Brendan J. Meade,et al.  Block models of crustal motion in southern California constrained by GPS measurements , 2005 .

[46]  J. Shaw,et al.  Uplift of the Longmen Shan and Tibetan plateau, and the 2008 Wenchuan (M = 7.9) earthquake , 2009, Nature.

[47]  Carl Tape,et al.  Multiscale estimation of GPS velocity fields , 2008 .