-
2020年中国进行了第3次国家珠穆朗玛峰(以下简称“珠峰”)高程测量[1],珠峰是否增高,2015年尼泊尔地震对珠峰有多大影响,这些问题受到社会和地学学者的广泛关注,珠峰高海拔地区地壳形变机制与特征成为地学研究的热点[2-6]。
珠峰位于中国青藏高原南端的喜马拉雅山脉,是印度板块与欧亚板块的交界处,是中国和世界最高峰。喜马拉雅山脉位于青藏高原最南端,是印度板块和欧亚板块汇聚的主要形变带,是研究大陆构造形变机制的最具代表性的区域,是地学研究的宝库[7-8]。喜马拉雅山脉特别是在珠峰南北两侧区域地壳形变具有较为明显的差异,呈现阶梯状分布特征。在珠峰南麓的尼泊尔南端属于恒河平原,向北进入海拔2 000 m以下的次喜马拉雅地区,位于加德满都以南的中坡度平缓地区的低喜马拉雅地区、加德满都以北的高喜马拉雅地区,再向北则进入喜马拉雅山脉以北的青藏高原。短短200多公里的距离,出现从平原到世界屋脊的巨大变化,跨越印度和欧亚两大板块[9-10]。由于板块的快速汇聚作用,使得该区域一直是一个强震活动带,有历史记录以来,珠峰周边地区多次发生震级Ms 8.0以上强震,最近一次是2015-04-25发生的Ms 8.1级地震,研究表明,这次地震对珠峰地区地壳形变趋势及珠峰高程都带来较大影响[11-13]。
研究跨喜马拉雅地区,特别是在珠峰周边地壳三维形变特征具有重要的科学价值和现实意义。地学学者针对跨喜马拉雅山脉进行了长期的研究,20世纪90年代以来,随着全球导航卫星系统(global navigation satellite system,GNSS)技术的发展,利用GNSS长期、连续、高精度的技术优势,地学学者在横跨喜马拉雅山脉布设了连续运行参考站(continuously operating reference stations,CORS)网,获得了大量的GNSS观测数据。陆态网络在中国青藏高原地区先后建立了数十个CORS,为研究跨喜马拉雅山脉地壳形变研究提供了丰富的数据资源[14]。
研究表明,喜马拉雅山脉不同区域地壳形变特征差异巨大,低喜马拉雅地区地势平缓,地壳垂直和水平运动与印度板块一致性较好,而到高喜马拉雅地区后,地形存在角度较大的陡坡带,其垂直形变变化较大,青藏高原由南到北则有明显的南北压缩、东西拉伸的区域形变特征[15-17]。为精细化研究该区域地壳形变特征,需要根据不同地质构造分段分区域研究地壳运动与形变。
本文利用珠峰南北两侧跨喜马拉雅山脉分布的CORS网多年GNSS的观测数据,以及2005年和2020年珠峰高程测量布设的形变监测网精确连续监测珠峰周边地区地壳形变动态变化,研究了珠峰地区地壳垂直形变趋势,定量分析珠峰周边区域2015年尼泊尔地震同震位移,以及地震对区域地壳垂直形变长期趋势的影响,其中,青藏高原研究区域主要集中在西藏块体[18-22]。
-
本文收集了2005年和2020年珠峰高程测量形变监测网GNSS数据,以及1999—2020年分布在跨喜马拉雅山脉的32个CORS连续观测数据,研究了该区域地壳长期形变趋势。CORS分布如图 1所示,其中,14个CORS位于尼泊尔,18个CORS中国青藏高原区域。
喜马拉雅山脉南北两侧CORS高程差异很大,如图 2所示,最低位于尼泊尔的BRN2站,大地高约20 m,最高位于中国青藏高原的XZSH站,大地高近5 000 m。这些CORS站点高程呈明显的阶梯状分布,可以精细化监测珠峰周边地壳分段三维形变特征。
本文利用GAMIT(GNSS at MIT)软件对GNSS数据进行单天解基线计算,解算方案主要参数的设置如表 1所示。
表 1 单天解基线解算设置
Table 1. Settings of Baseline Resolution
参数 处理方式 基线处理模式 RELAX 观测值 LC+PC组合 坐标框架 ITRF2014 对流层延迟 GMF 海潮改正 otl_FES2004 固体潮模型 IERS2010 光压模型 BERNE 解算获取各站单日解时间序列如图 3所示,单天解精度水平方向优于3 mm,高程方向优于5 mm。利用GLOBK(global Kalman filter)软件进行整体平差,获取区域CORS三维线性速率,平均精度优于0.1 mm/a,其中,NAST、CHLM、KKN4、SNDL等4个CORS因为受2015年尼泊尔地震影响,出现较大跳跃,在水平和高程方向上最大位移超过1 m,无法计算线性速率,故在1999—2020年速率解算中删除这4个站,见图 4。
-
1999—2020年CORS三维运动速率如图 5所示。站点速率水平方向均优于1 mm/a,垂直方向均优于1.5 mm/a。由图 5显示,喜马拉雅山脉南侧CORS水平方向东西速率变化不大,北向速率从南向北有比较明显的递减趋势。恒河平原的BRN2站北向运动速率约35 mm/a至高喜马拉雅地区的SYBC站约17 mm/a;喜马拉雅山脉北侧由南到北则有明显的东向运动速率逐次递增、北向运动速率逐次递增趋势,东向速率由XZZF站约37 mm/a至位于羌塘块体的QHTT站约48 mm/a,北向速率由XZZF站约25 mm/a至QHTT站约6.5 mm/a,整体呈右旋趋势。
垂直方向有明显的阶梯型特征。喜马拉雅山脉南侧,恒河平原和次喜马拉雅地区两个CORS(BRN2和GRIH站),高程方向呈略微下降趋势,其他CORS均呈上升趋势,总体上升为2~3 mm/a;位于喜马拉雅山脉北侧青藏高原各站点均呈上升趋势,总体隆升为3~5 mm/a。研究结果表明,青藏高原垂直方向长期保持整体隆升趋势。
为分析区域垂直形变特征,分别利用位于印度板块国际GNSS服务组织(International Global Navigation Satellite System Service,IGS)站(HYDE和IISC站)和欧亚板块稳定IGS站(DEAJ和ULAB站),分析珠峰地区相对印度板块和欧亚板块的高程变化如图 6所示。
图 6 1999—2020年相对于位于印度板块和欧亚板块上的IGS站的垂直形变
Figure 6. Vertical Deformation Relative to Indian Plate and Eurasian Plate from 1999 to 2020
由图 6(a)可知,由于1999—2020年HYDE和IISC站高程方向变化相对稳定,高程方向速率很小,平均约0.2 mm/a,因此,相对这两个IGS站,研究区域各CORS高程方向速率基本不变,区域整体垂直形变与在国际地球参考框架(international terrestrial reference frame,ITRF)下的一致。相对位于欧亚板块上的IGS站(DEAJ和ULAB站)站高程方向速率平均约2 mm/a,由图6(b)可知,在喜马拉雅山脉北侧的由低喜马拉雅地区向高喜马拉雅地区过渡地带的CORS相对DEAJ和ULAB站呈下降趋势,而在高喜马拉雅地区和青藏高原则仍在整体上呈上升趋势。研究表明,地壳整体隆升是珠峰地区乃至青藏高原总体运动趋势,其隆升速率远大于印度板块和欧亚板块。
-
2015-04-25,尼泊尔发生了Ms 8.1地震。本文利用CORS连续观测数据,分析了地震影响的同震位移;对比地震前后CORS速率,分析了地震对区域地壳形变趋势影响。
-
本文避开了主震及后Ms 7.0级余震后0.5 h及2015-04-26 Ms 5.0和Ms 7.1级余震。选用UTC时间2015-04-25T03:00:00—06:00:00和07:30:00—10:30:00、2015-04-24T18:00:00—24:00:00和2015-04-25T07:30:00—19:30:00、2015-04-24和2015-04-27、2015-04 22日—24日和2015-04 28日—30日的GNSS观测数据,地震前后3 h、12 h、24 h、72 h的CORS坐标变化见图 7。
由图 7结果显示,较大地表位移变化主要集中在地震发生的短时间内。地震对距离震中较近(约100 km以内)站点影响较大,其中,CHLM、KKN4、NAST等站同震位移超过1 m。而随着距离震中越远,则影响快速衰减,对于超过200 km站点,同震位移普遍为cm级。中国青藏高原CORS垂直方向普遍出现下降,地震后3 h普遍下降1 cm左右,最大约2 cm;地震前后3 d平差对比,最小下降约2.5 mm,最大下降约15 mm。
-
尼泊尔地震导致珠峰及周边地区地壳高程下降。长期是否会导致该地区隆升趋势发生变化,本文对地震前后时段(1999-01-01—2015-03-30、2015-07-01—2020-12-31)CORS三维速率变化进行了对比。为避免地震及其余震引起的跳变对长期趋势产生影响,地震前后两个时段中不包括2015年4月—6月数据,以分析地震对珠峰及周边地区形变趋势的影响。2015年地震前后CORS站水平和垂直方向速率差异见图 8,N、E和高程方向运动速率差异见图 9。图 8中,蓝色箭头为地震前站点速率,红色箭头为地震后站点速率。
图 8 2015年尼泊尔地震前后CORS水平和垂直方向速率差异
Figure 8. Velocities Difference of Horizontal and Vertical Directions of CORS Before and After 2015 Nepal Earthquake
图 9 2015年尼泊尔地震前后CORS N、E和高程方向运动速率差异
Figure 9. Velocities Difference of N, E and Height Directions of CORS Before and After the 2015 Nepal Earthquake
珠峰作为印度板块与欧亚板块碰撞汇聚的结果,其垂直形变与水平形变有一定的相关性。采用Delaunay三角形法分别计算了地震前后应变率场[23-24]如图 10和图 11所示。利用误差传播定律得到应变率精度[25],两期应变率精度在同一量级,东西向应变率中误差平均为1.5×10-8/a,南北向应变率中误差平均为1.8×10-8/a,剪应变率中误差平均为2.0×10-8/a。对比地震前后垂直形变特征,为分析地震前后垂直形变趋势变化动力学机制提供佐证资料。
图 10 1999—2015年应变率场(箭头表示主应变率)
Figure 10. Strain Rate Field of 1999 to 2015 (Arrow Indicates the Principal Strain Rate)
从图 10和图 11发现,在2015年尼泊尔地震前,该区域就有明显的挤压收缩形变,地震后印度板块北向运动趋势基本未变,而青藏高原特别是西藏块体北向运动减速,导致该区域挤压形变更为明显;在2015年尼泊尔地震震中至珠峰一线,有明显的应力集中现象,震后应变更为明显,这说明地震有可能导致了印度板块与欧亚板块相对汇聚加速。
由图 8和图 9可知,喜马拉雅山脉以南地区,震后高程方向速率上升趋势非常明显,平均隆升速率由平均约1 mm/a上升至约4 mm/a。其中,KKN4站除外,该站在地震同震中高程方向发生了127 mm的同震位移,震后出现明显的下降趋势,可能存在震后回弹现象。在青藏高原,除个别站点外,大部分站点震后速率上升,特别是西藏块体,上升趋势较为明显,由平均约2.3 mm/a上升至约3.6 mm/a。分析可能的原因为印度板块与欧亚板块加速汇聚导致该区域隆升趋势加速。
-
中国在2005年[26-27]和2020年[1]分别对珠峰高程进行了测量,重合点高程差异如图 12所示。测量期间在周边设立多个GNSS形变监测网,分别在2005年1月—5月和2020年1月—5月,在测绘部门设立的测量标志点,按照国家GPS B级网要求进行的数据采集,每个点连续观测3 d,每天观测8 h。
图 12 2020年与2005年珠峰高程测量重合点高程差异
Figure 12. Elevation Difference of Coincidence Points of Everest Elevation Survey in 2020 and 2005
选择距离珠峰较近、大地高4 000 m以上的18个重合站点,分别解算得到2005年点位和2020年点位在ITRF2014框架下的三维坐标。各站大地高中误差最大为±3.92 mm,平均为±2.94 mm。对点位大地高进行对比,以此分析在较长时间尺度里,珠峰周边区域垂直方向变化。研究发现,18个重合点大地高均出现不同程度上升,最大达到约3.9 cm,18个重合点平均上升约2.2 cm。
对比2005年和2020年珠峰高程测量结果,18个重合点均出现不同程度隆升,再次证明,长期珠峰地区整体保持隆升趋势。
-
本文利用珠峰形变监测网及阶梯型分布的CORS 1999—2020年GNSS连续观测数据,研究分析珠峰地区地壳长期形变特征,特别是垂直形变特征;并研究了2015年尼泊尔地震对该区域CORS短期位移影响及地壳长期垂直形变趋势的影响。
GNSS监测结果表明,珠峰及周边地区地壳垂直方向具有明显的阶梯型特征。喜马拉雅山脉南侧,恒河平原和次喜马拉雅地区地壳垂直形变很小,低喜马拉雅和高喜马拉雅地区均呈上升趋势,总体上升趋势约为2~3 mm/a。喜马拉雅山脉北侧青藏高原各站点全部呈上升趋势。从长期来看,青藏高原保持着整体隆升的趋势,总体隆升趋势为3~5 mm/a。
2015年地震前后地壳形变趋势变化显示,震后印度板块和欧亚板块交界处应力集中和挤压收缩形变更加显著,说明地震导致了印度板块与欧亚板块有加速汇聚趋势,该趋势有导致震后珠峰区域地壳垂直方向出现加速隆升趋势,该趋势对珠峰南侧影响更为显著,珠峰南侧震后地壳垂直形变由平均约1 mm/a上升至约4 mm/a,珠峰北侧则由平均约2.3 mm/a上升至约3.6 mm/a。
对比2020珠峰高程测量与2005年珠峰高程测量形变监测网GNSS重合点解算结果,再次证明珠峰地区地壳整体隆升是该区域地壳长期垂直形变趋势。
Vertical Deformation Characteristics Affected by Strong Earthquakes in Mount Qomolangma and Surrounding Areas
-
摘要: 2015年尼泊尔地震对珠穆朗玛峰高程的影响,近年一直受到全世界关注。2020年珠穆朗玛峰高程测量在珠穆朗玛峰及周边地区布设了高精度的全球导航卫星系统(global navigation satellite system, GNSS)形变监测网,收集了1999—2020年跨喜马拉雅山脉的32个连续运行参考站(continuously operating reference stations, CORS)的GNSS连续观测数据。利用GNSS数据监测了珠穆朗玛峰周边地区地壳三维形变特征,定量获取了2015年尼泊尔强震对珠穆朗玛峰周边CORS同震位移,以及地震对区域地壳三维形变长期趋势的影响,特别是对该地区垂直形变的影响。研究结果表明,该区域地壳垂直形变由南至北跨喜马拉雅山脉呈明显的阶梯型分布特征;震后印度板块与欧亚板块存在加速汇聚趋势,导致震后地壳隆升速率同步增大。Abstract:
Objectives The impact of the 2015 Nepal earthquake on the height of Mount Qomolangma has attracted worldwide attention in recent years. 2020 height measurement of Mount Qomolangma has set up a high-precision GNSS(global navigation satellite system) deformation monitoring network in Mount Qomolangma and its surrounding areas. In addition, we have collected the GNSS observation datum of 32 CORS(continuously operating reference stations) across the Himalayas during 1999 to 2020. Methods We monitor the three-dimensional crustal deformation characteristics of Mount Qomolangma and surrounding areas, and quantitatively obtained the co-seismic displacement of the CORS around the Mount Qomolangma of the 2015 Nepal earthquake, as well as the impact of the earthquake on the long-term trend of regional three-dimensional crustal deformation, especially on the vertical deformation in this region. Results The results show that the vertical crustal deformation in this area has obvious stepped distribution characteristics from south to north across the Himalayas. After the earthquake, there was an accelerated convergence trend between the Indian plate and the Eurasian plate, which led to the synchronous increase of the crustal uplift rate after the earthquake. Conclusions The crustal rise in Mount Qomolangma and surrounding areas is the main trend. The 2015 Nepal earthquake caused the crustal decline in a short time, but did not change the long-term uplift trend. -
表 1 单天解基线解算设置
Table 1. Settings of Baseline Resolution
参数 处理方式 基线处理模式 RELAX 观测值 LC+PC组合 坐标框架 ITRF2014 对流层延迟 GMF 海潮改正 otl_FES2004 固体潮模型 IERS2010 光压模型 BERNE -
[1] 党亚民, 郭春喜, 蒋涛, 等. 2020珠峰测量与高程确定[J]. 测绘学报, 2021, 50(4): 556-561 https://www.cnki.com.cn/Article/CJFDTOTAL-CHXB202104014.htm Dang Yamin, Guo Chunxi, Jiang Tao, et al. 2020 Height Measurement and Determination of Mount Qomolangma[J]. Acta Geodaetica et Cartographica Sinica, 2021, 50(4): 556-561 https://www.cnki.com.cn/Article/CJFDTOTAL-CHXB202104014.htm [2] Bilham R, Gaur V K, Molnar P. Himalayan Seismic Hazard[J]. Science, 2001, 293 (5534): 1442-1444 doi: 10.1126/science.1062584 [3] Tapponnier P, Xu Z Q, Roger F, et al. Oblique Stepwise Rise and Growth of the Tibet Plateau[J]. Science, 2001, 294(5547): 1671-1677 doi: 10.1126/science.105978 [4] Bilham R, Larson K, Freymueller J. GPS Measurements of Present-Day Convergence Across the Nepal Himalaya[J]. Nature, 1997, 386(6620): 61-64 doi: 10.1038/386061a0 [5] 许志琴, 杨经绥, 侯增谦, 等. 青藏高原大陆动力学研究若干进展[J]. 中国地质, 2016, 43(1): 1-42 https://www.cnki.com.cn/Article/CJFDTOTAL-DIZI201601001.htm Xu Zhiqin, Yang Jingsui, Hou Zengqian, et al. The Progress in the Study of Continental Dynamics of the Tibetan Plateau[J]. Geology in China, 2016, 43 (1): 1-42 https://www.cnki.com.cn/Article/CJFDTOTAL-DIZI201601001.htm [6] 杨强, 党亚民. 利用GPS速度场估算青藏高原地壳韧性层等效粘滞系数分布的研究[J]. 测绘学报, 2010, 39(5): 497-502 https://www.cnki.com.cn/Article/CJFDTOTAL-CHXB201005012.htm Yang Qiang, Dang Yamin. A Research About Effective Viscosity of Tibetan Plateau Lithosphere Viscoelastic Ductile Layer Using GPS Velocity Fields [J]. Acta Geodaetica et Cartographica Sinica, 2010, 39(5): 497-502 https://www.cnki.com.cn/Article/CJFDTOTAL-CHXB201005012.htm [7] 许志琴, 杨经绥, 戚学祥, 等. 印度/亚洲碰撞: 南北向和东西向拆离构造与现代喜马拉雅造山机制再讨论[J]. 地质通报, 2006, 25(Z1): 1-14 https://www.cnki.com.cn/Article/CJFDTOTAL-ZQYD2006Z1003.htm Xu Zhiqin, Yang Jingsui, Qi Xuexiang, et al. India-Asia Collision: A Further Discussion of N-S- and E-W-Trending Detachments and the Orogenic Mechanism of the Modern Himalayas[J]. Geological Bulletin of China, 2006, 25(Z1): 1-14 https://www.cnki.com.cn/Article/CJFDTOTAL-ZQYD2006Z1003.htm [8] 许志琴, 杨经绥, 李海兵, 等. 青藏高原与大陆动力学: 地体拼合、碰撞造山及高原隆升的深部驱动力[J]. 中国地质, 2006, 3(2): 221-238 https://www.cnki.com.cn/Article/CJFDTOTAL-DIZI200602001.htm Xu Zhiqin, Yang Jingsui, Li Haibing, et al. The Qinghai-Tibet Plateau and Continental Dynamics: A Review on Terrain Tectonics, Collisional Orogenesis, and Processes and Mechanisms for the Rise of the Plateau[J]. Geology in China, 2006, 33(2): 221-238 https://www.cnki.com.cn/Article/CJFDTOTAL-DIZI200602001.htm [9] 姚华建, 尹九洵. 喜马拉雅造山带构造特征与2015年尼泊尔Mw 7.8级大地震[J]. 科学通报, 2015, 60(27): 2656-2658 https://www.cnki.com.cn/Article/CJFDTOTAL-KXTB201527011.htm Yao Huajian, Yin Jiuxun. Structural Characteristics of Himalayan Orogenic Belt and the 2015 Nepal Mw 7.8 Earthquake[J]. Chinese Science Bulletin, 2015, 60(27): 2656-2658 https://www.cnki.com.cn/Article/CJFDTOTAL-KXTB201527011.htm [10] 刘静, 纪晨, 张金玉, 等. 2015年4月25日尼泊尔Mw 7.8级地震的孕震构造背景和特征[J]. 科学通报, 2015, 60(27): 2640-2655 https://www.cnki.com.cn/Article/CJFDTOTAL-KXTB201527010.htm Liu Jing, Ji Chen, Zhang Jinyu, et al. Tectonic Setting and General Features of Coseismic Rupture of the 25 April, 2015 Mw 7.8 Gorkha, Nepal Earthquake[J]. Chinese Science Bulletin, 2015, 60(27): 2640-2655 https://www.cnki.com.cn/Article/CJFDTOTAL-KXTB201527010.htm [11] 苏小宁, 王振, 孟国杰, 等. GPS观测的2015年尼泊尔Ms 8.1级地震震前应变积累及同震变形特征[J]. 科学通报, 2015, 60 (22): 2115-2123 https://www.cnki.com.cn/Article/CJFDTOTAL-KXTB201522007.htm Su Xiaoning, Wang Zhen, Meng Guojie, et al. Pre-seismic Strain Accumulation and Co-seismic Deformation of the 2015 Nepal Ms 8.1 Earthquake Observed by GPS[J]. Chinese Science Bulletin, 2015, 60(22): 2115-2123 https://www.cnki.com.cn/Article/CJFDTOTAL-KXTB201522007.htm [12] 王小瑞, 党引群, 高新妍, 等. 利用GPS监测网分析尼泊尔地震对珠峰地区及周边影响[J]. 全球定位系统, 2016, 41(3): 73-77 https://www.cnki.com.cn/Article/CJFDTOTAL-QUDW201603015.htm Wang Xiaorui, Dang Yinqun, Gao Xinyan, et al. Analysis Influence Everest Area Surrounding Nepal Earthquake with GPS Monitoring Network[J]. GNSS World of China, 2016, 41(3): 73-77 https://www.cnki.com.cn/Article/CJFDTOTAL-QUDW201603015.htm [13] 王虎, 李建成, 党亚民, 等. PPP网解UPD模糊度固定技术监测尼泊尔Ms 8.1级地震对中国珠峰地区及周边地震同震位移[J]. 测绘学报, 2016, 45(S2): 147-155 https://www.cnki.com.cn/Article/CJFDTOTAL-CHXB2016S2018.htm Wang Hu, Li Jiancheng, Dang Yamin, et al. Co-seismic Displacement of the 25 April 2015 Nepal Ms 8.1 Earthquake Effects on the China's Mount Everest Area Derived from GNSS Data Using the PPP Network Solution by UPD Ambiguity Fixed Technology[J]. Acta Geodaetica et Cartographica Sinica, 2016, 45(S2): 147-155 https://www.cnki.com.cn/Article/CJFDTOTAL-CHXB2016S2018.htm [14] 梁诗明. 基于GPS观测的青藏高原现今三维地壳运动研究[D]. 北京: 中国地震局地质研究所, 2014 Liang Shiming. Three-Dimensional Velocity Field of Present-Day Crustal Motion of the Tibetan Plateau Inferred from GPS Measurements[D]. Beijing: Institute of Geology, China Earthquake Administration, 2014 [15] Freymueller J, Bilham R, Bürgmann R, et al. Global Positioning System Measurements of Indian Plate Motion and Convergence Across the Lesser Himalaya[J]. Geophysical Research Letters, 1996, 23 (22): 3107-3110 doi: 10.1029/96GL02518 [16] Bilham R, Larson K, Freymueller J. GPS Measurements of Present-Day Convergence Across the Nepal Himalaya[J]. Nature, 1997, 386(6620): 61-64 doi: 10.1038/386061a0 [17] Bettinelli P, Avouac J P, Flouzat M, et al. Plate Motion of India and Interseismic Strain in the Nepal Himalaya from GPS and DORIS Measurements[J]. Journal of Geodesy, 2006, 80 (8/9/10/11): 567-589 [18] 张国民, 马宏生, 王辉, 等. 中国大陆活动地块边界带与强震活动[J]. 地球物理学报, 2005, 48(3): 602-610 https://www.cnki.com.cn/Article/CJFDTOTAL-DQWX200503017.htm Zhang Guomin, Ma Hongsheng, Wang Hui, et al. Boundaries Between Active-Tectonic Blocks and Strong Earthquakes in the China Mainland[J]. Chinese Journal of Geophysics, 2005, 48(3): 602-610 https://www.cnki.com.cn/Article/CJFDTOTAL-DQWX200503017.htm [19] 张培震, 邓起东, 张国民, 等. 中国大陆的强震活动与活动地块[J]. 中国科学: 地球科学, 2003, 33(S1): 12-20 https://www.cnki.com.cn/Article/CJFDTOTAL-JDXK2003S1001.htm Zhang Peizhen, Deng Qidong, Zhang Guomin, et al. Active Tectonic Blocks and Strong Earthquakes in Continent of China[J]. Science in China: Terrae, 2003, 33(S1): 12-20 https://www.cnki.com.cn/Article/CJFDTOTAL-JDXK2003S1001.htm [20] 邓起东, 张培震, 冉勇康, 等. 中国活动构造基本特征[J]. 中国科学: 地球科学, 2002, 32(12): 1020-1030 https://www.cnki.com.cn/Article/CJFDTOTAL-JDXK200212006.htm Deng Qidong, Zhang Peizhen, Ran Yongkang, et al. Basic Characteristics of Active Tectonics of China [J]. Science in China: Terrae, 2002, 32(12): 1020-1030 https://www.cnki.com.cn/Article/CJFDTOTAL-JDXK200212006.htm [21] 王敏, 沈正康, 牛之俊, 等. 现今中国大陆地壳运动与活动块体模型[J]. 中国科学: 地球科学, 2003, 33(S1): 21-32 https://www.cnki.com.cn/Article/CJFDTOTAL-JDXK2003S1002.htm Wang Min, Shen Zhengkang, Niu Zhijun, et al. The Present Crustal Movement of Chinese Mainland and Activities Block Model[J]. Science in China: Terrae, 2003, 33(S1): 21-32 https://www.cnki.com.cn/Article/CJFDTOTAL-JDXK2003S1002.htm [22] 党亚民, 杨强, 梁诗明, 等. 川滇区域活动块体运动与应变特征地震影响分析[J]. 测绘学报, 2018, 47(5): 559-566 https://www.cnki.com.cn/Article/CJFDTOTAL-CHXB201805002.htm Dang Yamin, Yang Qiang, Liang Shiming, et al. Block Movement and Strain Characteristics Effected by Earthquake in Sichuan-Yunnan Region[J]. Acta Geodaetica et Cartographica Sinica, 2018, 47(5): 559-566 https://www.cnki.com.cn/Article/CJFDTOTAL-CHXB201805002.htm [23] 江在森, 刘经南. 应用最小二乘配置建立地壳运动速度场与应变场的方法[J]. 地球物理学报, 2010, 53(5): 1109-1117 https://www.cnki.com.cn/Article/CJFDTOTAL-DQWX201005012.htm Jiang Zaisen, Liu Jingnan. The Method in Establishing Strain Field and Velocity Field of Crustal Movement Using Least Squares Collocation[J]. Chinese Journal of Geophysics, 2010, 53(5): 1109-1117 https://www.cnki.com.cn/Article/CJFDTOTAL-DQWX201005012.htm [24] 顾国华. 大地形变测量中的应变计算[J]. 地壳形变与地震, 1989, 9(1): 38-45 https://www.cnki.com.cn/Article/CJFDTOTAL-DKXB198901004.htm Gu Guohua. Some Problems in Calculating Strains from Geodetic Data[J]. Crustal Deformation and Earthquake, 1989, 9(1): 38-45 https://www.cnki.com.cn/Article/CJFDTOTAL-DKXB198901004.htm [25] 姚宜斌, 刘强, 江国焰, 等. 华北地区应变率及其精度评定[J]. 武汉大学学报·信息科学版, 2015, 40(10): 1317-1323 doi: 10.13203/j.whugis20140001 Yao Yibin, Liu Qiang, Jiang Guoyan, et al. Strain Rate and Its Accuracy Assessment in North China [J]. Geomatics and Information Science of Wuhan University, 2015, 40(10): 1317-1323 doi: 10.13203/j.whugis20140001 [26] 党亚民, 程传录, 陈俊勇, 等. 2005珠峰测高GPS测量及其数据处理[J]. 武汉大学学报·信息科学版, 2006, 31(4): 297-300 http://ch.whu.edu.cn/article/id/2422 Dang Yamin, Cheng Chuanlu, Chen Junyong, et al. GPS Data Processing of the 2005 Qomolangma Height Surveying[J]. Geomatics and Information Science of Wuhan University, 2006, 31 (4): 297-300 http://ch.whu.edu.cn/article/id/2422 [27] 党亚民. 珠峰高程复测有关问题的探讨[J]. 测绘科学, 2005, 30(3): 101-103 https://www.cnki.com.cn/Article/CJFDTOTAL-CHKD200503031.htm Dang Yamin. Investigation on the Height Repetition Determination of Qomolangma Peak[J]. Science of Surveying and Mapping, 2005, 30(3): 101-103 https://www.cnki.com.cn/Article/CJFDTOTAL-CHKD200503031.htm -