THE VERTICAL COSEISMIC DEFORMATION FIELD OF THE WENCHUAN EARTHQUAKE BASED ON THE COMBINATION OF GPS AND INSAR MEASUREMENTS

doi:10.3969/j.issn.0253-4967.2014.03.014

Abstract

Abstract:

Vertical coseisimic deformation near seismogenic fault is one of the most important parameters for understanding the fault behavior, especially for thrust or normal fault, since near field vertical deformation provides meaningful information for understanding the rupture characteristics of the seismogenic fault and focal mechanism. Taking Wenchuan thrust earthquake for an example, we interpolate GPS horizontal observed deformation using Biharmonic spline interpolation and derive them into east-westward or north-southward deformation field. We first use reliable GPS observed value to correct InSAR reference point and unify both GPS and InSAR coordinate frame. We then make a profile using InSAR data and compare it to that from GPS data and we find GPS and InSAR observation reference point has a 9.93cm difference in the hanging wall side, and around -11.49cm in the footwall. After correction, we obtain a continuous vertical deformation field of the Wenchuan earthquake by combined calculation of GPS and InSAR LOS deformation field. The results show that the vertical deformation of both hanging wall and foot wall of the fault decreases rapidly, with deformation greater than 30cm within 50km across the fault zone. The uneven distribution of the vertical deformation has some peak values at near fault, mainly distributed at the southern section(the town of Yingxiu), the middle(Beichuan)and the northern end(Qingchuan)of the seismogenic fault. These three segments have their own characteristics. The southern section of the fault has an obvious asymmetric feature, which exhibits dramatic uplift reaching 550cm on the hanging wall, with the maximum uplift area located in Yingxiu town to Lianshanping. The middle section shows a strong anti-symmetric feature, with one side uplifting and the other subsiding. The largest uplifting of the southern segment reaches around 255cm, located at the east of Chaping, and the largest subsiding is in Yongqing, reaching around -215cm. The vertical deformation of the northern section is relatively small and distributed symmetrically mainly in the north of Qingchuan, with the maximum uplift to be 120cm, locating in the northernmost of the seismogenic fault.

Key words: Wenchuan earthquake, Biharmonic spline interpolation, combination of GPS and InSAR, vertical deformation field

摘要：

发震断层的形变是断层活动的重要参数之一，对认识断层性质、震源机理有重要作用。文中以逆冲性质为主的汶川地震为例，采用符合地表水平形变特征的Biharmonic样条插值对GPS水平形变矢量插值，然后再分解为EW和SN向分量。利用可靠的GPS观测值对InSAR参考点进行校正，统一两者的坐标系。通过对汶川地震视线向形变场剖面与GPS对比分析发现，断层上盘GPS与InSAR观测参考点相差9.93cm，而下盘则为-11.49cm。在此研究基础上，通过GPS水平形变场与InSAR视线向形变场联合解算，获取了汶川地震垂直连续形变场。结果表明，断层两侧垂直形变衰减较快，横跨断裂带形变量>30cm的宽度不超过50km；沿发震断层附近垂直形变高值区分布不均匀，主要集中分布在发震断裂的汶川县城至都江堰段、茶坪—北川—南坝段和青川段。这3段各有特色，南段断层两侧垂直形变极不对称，主要以上盘剧烈抬升为主，最大抬升区域在映秀镇至连山坪一带，抬升量达到5.5m。中段表现为较强的反对称性，断层一侧抬升另一侧沉降。该段上盘最大抬升区域在茶坪东侧，抬升量为255cm，下盘最大沉降量在永庆，沉降量为-215cm。北端垂直形变相对较小，主要分布在青川北侧，呈对称分布，在发震断层最北端，最大抬升量为120cm。

关键词: 汶川地震, Biharmonic样条插值, InSAR和GPS结合, 垂直形变场

CLC Number:

P315.72⁺5

SHAN Xin-jian, QU Chun-yan, GUO Li-min, ZHANG Guo-hong, SONG Xiao-gang, ZHANG Gui-fang, WEN Shao-yan, WANG Chi-sheng, XU Xiao-bo, LIU Yun-hua. THE VERTICAL COSEISMIC DEFORMATION FIELD OF THE WENCHUAN EARTHQUAKE BASED ON THE COMBINATION OF GPS AND INSAR MEASUREMENTS[J]. SEISMOLOGY AND GEOLOGY, 2014, 36(3): 718-730.

单新建, 屈春燕, 郭利民, 张国宏, 宋小刚, 张桂芳, 温少妍, 汪池升, 徐小波, 刘云华. 基于InSAR与GPS观测的汶川同震垂直形变场的获取[J]. 地震地质, 2014, 36(3): 718-730.

References

邓起东, 陈社发, 赵小麟. 1994. 龙门山及其邻区的构造和地震活动及动力学
[J]. 地震地质, 16(4): 389-403.
DENG Qi-dong, CHEN She-fa, ZHAO Xiao-lin. 1994. Tectonics, seismicity and dynamics of Longmenshan Mountains and its adjacent regions
[J]. Seismology and Geology, 16(4): 389-403(in Chinese).
国家重大科学工程 "中国地壳运动观测网络"项目组. 2008. GPS测定的2008年汶川M_S 8.0地震的同震位移场
[J]. 中国科学(D辑), 38(10): 1195-1206.
Working Group of Crustal Movement Observation Network of China(CMONOC). 2008. Coseismic displacement field of the M_S 8.0 Wenchuan earthquake in 2008 measured by GPS surveys
[J]. Science in China(Ser D), 38(10): 1195-1206(in Chinese).
李元金, 罗立民, 张鹏程, 等. 2011. 基于校正靶特征与Biharmonic 样条插值的XRII 图像扭曲校正
[J]. 东南大学学报(自然科学版), 41(6): 1213-1218.
LI Yuan-jin, LUO Li-min, ZHANG Peng-cheng, et al. 2011. Distortion correction of XRII image based on calibration grid characteristic and Biharmonic interpolation
[J]. Journal of Southeast University(Natural Science Edition), 41(6): 1213-1218(in Chinese).
陆鹏. 2011. 基于Biharmonic样条的离散数据三维地形生成方法
[J]. 广西水利水电, 2011(6): 64-73.
LU Peng. 2011. 3D terrain generation out of discrete data by Biharmonic spline method
[J]. Gunagxi Water Resource & Hydropower Engineering, 2011(6): 64-73(in Chinese).
单新建, 屈春燕, 宋小刚, 等. 2009. 汶川M_S 8.0地震InSAR同震形变场观测与研究
[J]. 地球物理学报, 52(2): 496-504.
SHAN Xin-jian, QU Chun-yan, SONG Xiao-gang, et al. 2009. Coseismic surface deformation caused by the Wenchuan M_S 8.0 earthquake from InSAR data analysis
[J]. Chinese J Geophys, 52(2): 496-504(in Chinese).
汪驰升, 单新建, 张国宏, 等. 2009. 基于ASAR升降轨数据解算于田M_S 7.3地震3D同震形变场
[J]. 地震, 29(B10): 105-112.
WANG Chi-sheng, SHAN Xin-jian, ZHANG Guo-hong, et al. 2009. 3D coseismic deformation field of the Yutian M_S7.3 earthquake calculated from ASAR ascending and descending data
[J]. Earthquake, 29(B10): 105-112(in Chinese).
王亚涛, 董兰芳, 倪奎. 2007. 基于Biharmonic样条插值的图像渐变算法及实现
[J]. 中国图像图形学报, 12(12): 2189-2194.
WANG Ya-tao, DONG Lan-fang, NI Kui. 2007. Image morphing algorithm based on Biharmonic spline interpolation and its implementation
[J]. Journal of Image and Graphics, 12(12): 2189-2194(in Chinese).
张培震, 徐锡伟, 闻学泽, 等. 2008. 2008年汶川8.0级地震发震断裂的滑动速率、复发周期和构造成因
[J]. 地球物理学报, 51(4): 1066-1073.
ZHANG Pei-zhen, XU Xi-wei, WEN Xue-ze, et al. 2008. Slip rates and recurrence intervals of the Longmenshan active fault zone and tectonic implications for the mechanism of the May 12 Wenchuan earthquake, 2008, Sichuan, China
[J]. Chinese J Geophys, 51(4): 1066-1073(in Chinese).
Hu Jun, Li Zhi-wei, Ding Xiao-li, et al. 2008. Two-dimensional co-seismic surface displacements field of the Chi-Chi earthquake inferred from SAR image matching
[J]. Sensors, 8: 6484-6495.
Qu Chun-yan, Shan Xin-jian, Zhang Guo-hong, et al. 2010. Coseismic displacement field of the Wenchuan M_S8.0 earthquake in 2008 derived using differential radar interferometry
[J]. Journal of Applied Remote Sensing, 4(1): 043516(19 March 2010). doi: 10.1117/1.3386043.
Samsonov S, Tiampo K. 2006. Analytical optimization of a DInSAR and GPS dataset for derivation of three-dimensional surface motion
[J]. IEEE Geosci Remote Sens Lett, 3(1): 107-111.
Sandwell D T. 1987. Biharmonic spline interpolation of GEOS-3 and SEASAT altimeter data
[J]. Geophysical Research Letters, 14(2): 139-142.
Shen Zheng-kang, Sun Jian-bao, Zhang Pei-zhen, et al. 2009. Slip maxima at fault junctions and rupturing of barriers during the 2008 Wenchuan earthquake
[J]. Nature Geoscience, 2: 718-724, Published online: 27 September 2009. doi: 10.1038/ngeo636.
Sverrir G, Freysteinn S. 2002. Three-dimensional surface motion maps estimated from combined interferometric synthetic aperture radar and GPS data
[J]. Journal of Geophysical Research, B10(107): 1029.
Wang Qi, Qiao Xue-jun, Lan Qi-gui, et al. 2011. Rupture of deep faults in the 2008 Wenchuan earthquake and uplift of the Longmen Shan
[J]. Nature Geoscience, 4: 634-640. doi: 10.1038/nego1210.
Wang Qing-liang, Cui Du-xin, Zhang Xi, et al. 2009. Coseismic vertical deformation of the M_S8.0 Wenchuan earthquake from repeated levelings and its constraint on listric fault geometry
[J]. Earthquake Science, 22: 595-602.
Wright T J, Parsons B E, Lu Z. 2004. Toward mapping surface deformation in three dimensions using InSAR
[J]. Geophys Res Lett, 31: L01607. doi: 10.1029/2003GL018827.
Xu Xi-wei, Wen Xue-ze, Yu Gui-hua, et al. 2009. Coseismic reverse- and oblique-slip surface faulting generated by the 2008 M_W7.9 Wenchuan earthquake, China
[J]. Geology, 37(6): 515-518.
Yuri F, Mark S. 2001. The complete(3-D)surface displacement field in the epicentral area of the 1999 M_W7.1 Hector Mine earthquake, California, from space geodetic observation
[J]. Geophysical Research Letters, 28(16): 3063-3066.

[1]	SONG Dong-mei, WANG Hui, SHAN Xin-jian, WANG Bin, CUI Jian-yong. STRESS-INDUCED HEATING HYPOTHESIS BASED ON CORRELATION ANALYSIS OF GRAVITY AND THERMAL FIELDS BEFORE WENCHUAN EARTHQUAKE [J]. SEISMOLOGY AND GEOLOGY, 2023, 45(5): 1112-1128.
[2]	SONG Dong-mei, WANG Hui, SHAN Xin-jian, WANG Bin, CUI Jian-yong. A NOVEL EXTRACTION METHOD OF PRE-EARTHQUAKE GRACE GRAVITY ANOMALY INFORMATION BASED ON MAXIMUM SHEAR STRAIN [J]. SEISMOLOGY AND GEOLOGY, 2022, 44(6): 1539-1556.
[3]	CHEN Li-juan, CHEN Xue-zhong, LI Yan-e, GONG Li-wen. RELATIONSHIP BETWEEN DECREASING AMPLITUDE OF b-VALUE AND THE SEISMOGENIC ZONE OF THE WENCHUAN M_S8.0 EARTHQUAKE [J]. SEISMOLOGY AND GEOLOGY, 2022, 44(4): 1046-1058.
[4]	WANG Xiao-shan, WAN Yong-ge. CHARACTERISTICS OF THE CRUSTAL STRESS FIELD AND ITS DIRECTION CONVERGENCE BEFORE THE WENCHUAN EARTHQUAKE [J]. SEISMOLOGY AND GEOLOGY, 2022, 44(2): 363-377.
[5]	GUO Shu-song, ZHU Yi-qing, XU Yun-ma, LIU Fang, ZHAO Yun-feng, ZHANG Guo-qing, ZHU Hui. GRAVITY EVIDENCE OF META-INSTABLE STATE BEFORE THE 2008 WENCHUAN EARTHQUAKE [J]. SEISMOLOGY AND EGOLOGY, 2021, 43(6): 1368-1380.
[6]	LAN Jian, CHEN Xiao-li. EVOLUTION CHARACTERISTICS OF LANDSLIDES TRIGGERED BY 2008 M_S8.0 WENCHUAN EARTHQUAKE IN YINGXIU AREA [J]. SEISMOLOGY AND GEOLOGY, 2020, 42(1): 125-146.
[7]	ZHU Chuan-hua, SHAN Xin-jian, ZHANG Guo-hong, JIAO Zhong-hu, ZHANG Ying-feng, LI Yan-chuan, QIAO Xin. TWO DIMENSIONAL MODEL ON RELATION BETWEEN THERMAL ANOMALY BEFORE WENCHUAN EARTHQUAKE AND TECTONIC STRESS [J]. SEISMOLOGY AND GEOLOGY, 2019, 41(6): 1497-1510.
[8]	SONG Dong-mei, XIANG Liang, SHAN Xin-jian, YIN Jing-yuan, WANG Bin, CUI Jian-yong. THE METHOD OF CONSTRUCTING IONOSPHERIC TEC BACKGROUND FIELD BASED ON SVR MODEL [J]. SEISMOLOGY AND GEOLOGY, 2019, 41(6): 1511-1528.
[9]	MA Si-yuan, XU Chong, WANG Tao, LIU Jia-mei. APPLICATION OF TWO SIMPLIFIED NEWMARK MODELS TO THE ASSESSMENT OF LANDSLIDES TRIGGERED BY THE 2008 WENCHUAN EARTHQUAKE [J]. SEISMOLOGY AND GEOLOGY, 2019, 41(3): 774-788.
[10]	XU Zhi-ping, WANG Fu-yun, JIANG Lei, ZHAO Yan-na, YANG Li-pu, TANG Lin. THE THREE DIMENSIONAL DENSITY STRUCTURE OF CRUST AND UPPER MANTLE IN THE CENTRAL-SOUTHERN PART OF LONGMENSHAN [J]. SEISMOLOGY AND GEOLOGY, 2019, 41(1): 84-98.
[11]	YIN De-yu, LIU Qi-fang, LIU Chang, JI Xin-yang. ESTIMATING THE WENCHUAN EARTHQUAKE RUPTURE PROCESS USING NEAR FIELD STRONG MOTION RECORDS AND COSEISMIC DISPLACEMENTS [J]. SEISMOLOGY AND GEOLOGY, 2018, 40(3): 698-717.
[12]	ZHAO You-jia, ZHANG Guo-hong, ZHANG Ying-feng, SHAN Xin-jian, QU Chun-yan. TWO-DIMENSIONAL WHOLE CYCLE SIMULATION OF SPONTANE-OUS RUPTURE OF THE 2008 WENCHUAN EARTHQUAKE USING THE CONTINUOUS-DISCRETE ELEMENT METHOD [J]. SEISMOLOGY AND GEOLOGY, 2018, 40(1): 12-26.
[13]	LIU Yuan-zheng, MA Jin, MA Wen-tao, JIANG Tong. INFLUENCE OF PRESSURE HEAD CHANGE AND ITS CHANGE RATE ON RESERVOIR TRIGGERED SEISMICITY -A CASE STUDY OF ZIPINGPU RESERVOIR [J]. SEISMOLOGY AND GEOLOGY, 2017, 39(3): 437-450.
[14]	YAO Lu, MA Sheng-li, WANG Yu, HE Hong-lin, CHEN Jian-ye, YANG Xiao-song, SHIMAMOTO Toshihiko. THE VITRINITE REFLECTANCE OF FAULT ROCKS FROM THE WENCHUAN EARTHQUAKE FAULT ZONE: CONSTRAINTS ON FRICTIONAL PROPERTIES OF THE FAULT DURING THE EARTHQUAKE [J]. SEISMOLOGY AND GEOLOGY, 2016, 38(4): 817-829.
[15]	WEI Ben-yong, SU Gui-wu. ASSESSMENT ON INDIRECT ECONOMIC LOSS OF WENCHUAN EARTHQUAKE DISASTER BASED ON INPUT-OUTPUT ANALYSIS [J]. SEISMOLOGY AND GEOLOGY, 2016, 38(4): 1082-1094.