留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

应用GPS连续观测修正流动观测垂向速率的分析与实验

占伟 李斐 朱爽 张杰

占伟, 李斐, 朱爽, 张杰. 应用GPS连续观测修正流动观测垂向速率的分析与实验[J]. 武汉大学学报 ● 信息科学版, 2016, 41(7): 963-968. doi: 10.13203/j.whugis20140251
引用本文: 占伟, 李斐, 朱爽, 张杰. 应用GPS连续观测修正流动观测垂向速率的分析与实验[J]. 武汉大学学报 ● 信息科学版, 2016, 41(7): 963-968. doi: 10.13203/j.whugis20140251
ZHAN Wei, LI Fei, ZHU Shuang, ZHANG Jie. Analysis and Test of Correction to Vertical Velocity Measured by GPS Mobile Observation with Continuous Observation[J]. Geomatics and Information Science of Wuhan University, 2016, 41(7): 963-968. doi: 10.13203/j.whugis20140251
Citation: ZHAN Wei, LI Fei, ZHU Shuang, ZHANG Jie. Analysis and Test of Correction to Vertical Velocity Measured by GPS Mobile Observation with Continuous Observation[J]. Geomatics and Information Science of Wuhan University, 2016, 41(7): 963-968. doi: 10.13203/j.whugis20140251

应用GPS连续观测修正流动观测垂向速率的分析与实验

doi: 10.13203/j.whugis20140251
基金项目: 

中国地震局地震行业科研专项基金 No.201308009

国家自然科学基金 No.41474002

详细信息
    作者简介:

    占伟,博士生,研究方向为GPS数据处理与应用。zw000373@163.com

  • 中图分类号: P228.41

Analysis and Test of Correction to Vertical Velocity Measured by GPS Mobile Observation with Continuous Observation

Funds: 

The Special Fund for China Earthquake Research Project No.201308009

the National Natural Science Foundation of China No.41474002

  • 摘要: 目前,基于GPS流动观测确定垂向变化速率受年周期运动影响较大,但因观测数据较少,不能有效消除周期性运动的影响,使得GPS流动观测得到的垂向速率精度较低。提出利用同一区域GPS连续观测修正流动观测垂向速率的方法,首先依据GPS连续站垂向时间序列计算得到连续站的年周期运动振幅,然后进行空间插值得到各流动站的年周期运动振幅,进而对流动站垂向时间序列进行周期性运动修正,最后依据修正后的流动站垂向时间序列求取垂向线性速率。分别选用云南地区陆态网络GPS连续(2010~2013年)和流动观测数据(2011~2013年)进行了验证。结果表明,在年周期运动空间一致性较好的区域,该方法能够有效修正GPS流动观测垂向速率,提高垂向速率的解算精度。
  • 图  1  云南地区GPS连续站年周期运动 (2010~2013年,参考点2010.0)

    Figure  1.  Annual Terms of Vertical Position Variations for GPS Fiducial Stations(from 2010 to 2013, the Reference Time as 2010.0) in Yunnan Area

    图  2  云南地区GPS垂向速度场(流动观测,2011~2013年)

    Figure  2.  GPS Vertical Velocity Fields in Yunnan Area (Mobile Observations,from 2011 to 2013)

    图  3  云南地区GPS垂向速度场(连续观测,2010~2013年)

    Figure  3.  GPS Vertical Velocity Fields in Yunnan Area (Continuous Observations,from 2010 to 2013)

    表  1  云南地区GPS连续站垂向运动信息(2010~2013年)

    Table  1.   Vertical Motion for GPS Fiducial Station in Yunnan Area from 2010 to 2013

    点名速率 /(mm\5a-1)年周期观测天数
    振幅/mm初相/(°)
    KMIN2.939.4517.611 305
    XIAG0.179.5357.31 365
    YNCX*1.528.554.111 108
    YNDC*2.219.449.891 050
    YNGM-0.246.35350.08988
    YNHZ1.218.298.56973
    YNJD1.0410.978.311126
    YNJP-1.98.8817.381 134
    YNLA-0.6814.319.051 147
    YNLC*0.611.495.491 008
    YNLJ*1.3910.5618.841 115
    YNMH-1.0211.67.861 150
    YNMJ-0.39.890.151 129
    YNML1.07712.251 149
    YNMZ*-3.1710.6524.81 163
    YNRL-0.6411.810.571 087
    YNSD*2.1812.2711.65975
    YNSM*-1.5410.945.441 162
    YNTC0.8412.736.871 149
    YNTH*1.168.2717.311 164
    YNWS3.178.4721.88968
    YNXP2.157.956.631 215
    YNYA2.389.6814.571 006
    YNYL1.9512.8516.351 163
    YNYM*1.4510.2217.141 121
    YNYS2.1111.1112.931 160
    YNZD0.868.9722.98980
    注: 标*的点为检核站,未标注的点为基准站。
    下载: 导出CSV

    表  2  检核站解算结果

    Table  2.   Results of Check Sites

    检核站年周期运动振幅/mm速率结果/(mm\5a-1)
    ABA0B0V0V1V2V1-V0V2-V0
    YNCX8.530.619.551.661.52±0.122.04±0.181.55±0.120.510.03
    YNDC9.301.628.251.392.21±0.191.21±0.252.09±0.19-1.01-0.12
    YNLC11.441.108.950.230.6±0.151.21±0.260.71±0.160.610.11
    YNLJ9.993.4110.432.381.39±0.111.29±0.21.24±0.1-0.10-0.15
    YNMZ9.664.478.172.36-3.17±0.16-3.01±0.23-3.22±0.160.16-0.06
    YNSD12.022.4811.001.282.18±0.17-0.22±0.291.86±0.17-2.40-0.32
    YNSM10.891.0410.811.01-1.54±0.12-1.14±0.21-1.54±0.120.400.00
    YNTH7.902.468.091.631.16±0.111.34±0.171.1±0.120.18-0.05
    YNYM9.763.019.362.101.45±0.131.31±0.211.34±0.13-0.14-0.11
    下载: 导出CSV
  • [1] Wang Qi, Zhang Peizhen, Freymueller J T, et al. Present-day Crustal Deformation in China Constrained by Global Positioning System Measurements[J]. Science , 2001, 294(5 542):574-577
    [2] 梁洪宝,占伟,杨博,等.新疆地区现今水平形变及其动态变化[J]. 大地测量与地球动力学, 2012, 32(6):24-28

    Liang Hongbao, Zhan Wei, Yang Bo, et al. Present-day Horizontal Deformation and Its Dynamic Character in Xinjiang Area[J]. Journal of Geodesy and Geodynamics , 2012, 32(6):24-3
    [3] 武艳强,江在森,杨国华,等.南北地震带北段近期地壳变形特征研究[J]. 武汉大学学报·信息科学版, 2012, 37(9):1 045-1 049

    Wu Yanqiang, Jiang Zaisen,Yang Guohua,et al. Deformation Characteristics of North Section of North-south Seimic Zone in Recent Period[J]. Geomatics and Information Science of Wuhan University , 2012, 37(9):1 045-1 049
    [4] 刘志广,占伟,杨博,等.青藏高原中南部近期地壳水平形变[J]. 大地测量与地球动力学, 2013, 33(3):1-4

    Liu Zhiguang, Zhan Wei, Yang Bo, et al. Recent Horizontal Crust Deformation in Central and Southern Qinghai-Tibet Plateau[J]. Journal of Geodesy and Geodynamics , 2013, 33(3):1-4
    [5] Wu Yanqiang, Jiang Zaiseng, Wang Ming, et al. Preliminary Results Pertaining to Coseismic Displacement and Preseismic Strain Accumulation of the Lushan MS7.0 Earthquake, as Reflected by GPS Surveying[J]. Chinese Science Bulletin , 2013, 58: 3 460-3 466
    [6] Liang Shiming, Gan Weijun, Shen Chunzheng, et al. Three-dimensional Velocity Field of Present-Day Crustal Motion of the Tibetan Plateau Derived from GPS Measurements[J]. Journal of Geophysical Research , 2013, 118(10):5 722-5 732
    [7] Meng Xiangang, Liu Zhiguang. The Ms7.0 Lushan Earthquake and the Activity of the Longmenshan Fault Zone[J]. Geodesy and Geodynamics , 2013, 4(3):40-47
    [8] Dong D N, Fang P, Bock Y, et al. Anatomy of Apparent Seasonal Variations from GPS-derived Site Position Time Series[J]. Journal of Geophysical Research , 2002, doi: 10.1029/2001JB000573
    [9] 王敏,沈正康,董大南.非构造形变对GPS连续站位置时间序列的影响和修正[J].地球物理学报, 2005, 48(5):1 045-1 052

    Wang Min, Shen Zhengkang, Dong Danan. Effects of Non-Tectonic Crustal Deformation on Continuous GPS Position Time Series and Correction to Them[J]. Chinese J Geophys , 2005, 48 (5):1 045-1 052
    [10] Wdowinski S, Bock Y, Zhang J, et al.Southern California Permanent GPS Geodetic Array: Spatial Filtering of Daily Positions for Estimating Coseismic and Post Seismic Displacement Induced by the Landers Earthquake[J]. Journal of Geophysical Research , 1997, 102(B8):18 057-18 070
    [11] 盛传贞,甘卫军,梁诗明,等. 滇西地区GPS时间序列中陆地水载荷形变干扰的GRACE分辨与剔除[J]. 地球物理学报, 2014, 57(1):42-52

    Sheng Chuanzhen, Gan Weijun, Liang Shiming, et al. Identification and Eliminationof Non-Tectonic Crustal Deformation Caused by Land Water from GPS Time Series in the Western Yunnan Province Based on GRACE Observations[J]. Chinese J Geophys , 2014, 57(1):42-52
    [12] Herring T A, King R W, McClusky S C.GAMIT Reference Manual: GPS Analysis at MIT, Version 10.4[R]. Massachusetts Institute of Technology, Cambridge, 2010
    [13] 郝明.基于精密水准数据的青藏高原东缘现今地壳垂直运动与典型地震同震及震后垂直形变研究[D]. 北京:中国地震局地质研究所,2012

    Hao Ming. Present Crustal Vertical Movement of Eastern Tibetan Plateau and Coseismic and Postseimic Vertical Deformation of Two Typical Earthquakes[D]. Beijing: Institute of Geology, Chinese Earthquake Administration, 2012
    [14] 郝明,王庆良,崔笃信,等.滇西南龙陵-澜沧断裂带现今地壳垂直运动研究[J]. 地震研究, 2013, 36(3):281-285

    Hao Ming, Wang Qingliang, Cui Duxin, et al. Present Crustal Vertical Movement of Longling-Lancang Fault in Southeast of Yunnan Province[J]. Journal of Seismological Research , 2013, 36(3):281-285
    [15] 刘任莉,李建成,姜卫平,等. 联合GRACE与GPS比较山西省垂向地表形变[J]. 武汉大学学报·信息科学版,2013,38(4):426-430

    Liu Renli, Li Jiancheng, Jiang Weiping, et al. Comparing Vertical Surface Displacements Using GRACE and GPS over Shanxi Province[J]. Geomatics and Information Science of Wuhan University, 2013,38(4):426-430
    [16] Tesmer V, Steigenberger P, Dam T V, et al. Vertical Deformations from Homogeneously Processed GRACE and Global GPS Long-Term Series[J]. Journal of Geodesy, 2011, 85(5):291-310
    [17] Fu Y, Freymueller J T, Jens en T. Seasonal Hydrological Loading in Southern Alaska Observed by GPS and GRACE[J]. Geophysical Research Letters, 2012, 39(15):151-155
    [18] Zou Rong, Freymueller J T, Ding Kaihua, et al. Evaluating Seasonal Loading Models and Their Impact on Global and Regional Reference Frame Alignment[J]. Journal of Geophysical Research Solid Earth, 2014, 119(2):1 337-1 358
  • [1] 顾国华, 王武星.  GPS测得的2018年夏威夷6.9级地震与火山喷发地壳运动 . 武汉大学学报 ● 信息科学版, 2019, 44(8): 1191-1197, 1204. doi: 10.13203/j.whugis20180463
    [2] 顾国华, 王武星.  2016年新西兰7.8级大地震GPS观测结果与弹性回跳模型 . 武汉大学学报 ● 信息科学版, 2017, 42(11): 1673-1680. doi: 10.13203/j.whugis20170286
    [3] 周东旭, 周兴华, 张化疑, 王朝阳, 唐秋华.  利用GPS连续观测进行中国沿海验潮站地壳垂直形变分析 . 武汉大学学报 ● 信息科学版, 2016, 41(4): 516-522. doi: 10.13203/j.whugis20140714
    [4] 魏二虎, 刘文杰, WEI Jianan, 金双根, 刘经南.  VLBI和GPS观测联合解算地球自转参数和日长变化 . 武汉大学学报 ● 信息科学版, 2016, 41(1): 66-71,92. doi: 10.13203/j.whugis20130435
    [5] 陈德忠, 叶世榕, 刘炎炎, 刘郑.  基于观测值域的GPS多路径误差应用分析 . 武汉大学学报 ● 信息科学版, 2014, 39(2): 147-151. doi: 10.13203/j.whugis20120719
    [6] 王振杰, 聂志喜, 欧吉坤.  一种基于TurboEdit改进的GPS双频观测值周跳探测方法 . 武汉大学学报 ● 信息科学版, 2014, 39(9): 1017-1021. doi: 10.13203/j.whugis20130021
    [7] 魏二虎, 李广文, 畅柳, 曹起.  利用GPS观测数据研究高频地球自转参数 . 武汉大学学报 ● 信息科学版, 2013, 38(7): 818-821.
    [8] 魏文薪, 江在森, 武艳强, 赵静.  利用GPS数据研究川滇块体东边界主要断裂带运动特性 . 武汉大学学报 ● 信息科学版, 2012, 37(9): 1041-1044.
    [9] 瞿伟, 张勤, 王庆良, 李振洪.  利用GPS观测资料分析山西地堑系现今地壳构造活动特征 . 武汉大学学报 ● 信息科学版, 2012, 37(2): 195-198.
    [10] 艾松涛, 王泽民, 鄂栋臣, 闫明.  利用GPS的北极冰川运动监测与分析 . 武汉大学学报 ● 信息科学版, 2012, 37(11): 1337-1340.
    [11] 李一鹤, 沈云中.  GPS观测值的时间相关性对基线解算影响 . 武汉大学学报 ● 信息科学版, 2011, 36(4): 427-430.
    [12] 吴继忠.  基于GPS观测的Baja California地震地壳变形分析 . 武汉大学学报 ● 信息科学版, 2011, 36(4): 437-440.
    [13] 李博峰, 沈云中, 楼立志.  GPS中长基线观测值随机特性分析 . 武汉大学学报 ● 信息科学版, 2010, 35(2): 176-180.
    [14] 乔学军, 陈顒, 王琪, 杨少敏.  首都圈地区现今地壳运动的GPS观测与构造活动模拟 . 武汉大学学报 ● 信息科学版, 2008, 33(7): 692-696.
    [15] 李军, 王继业, 熊熊, 许厚泽.  东北亚地区GPS观测数据质量检测和分析 . 武汉大学学报 ● 信息科学版, 2006, 31(3): 209-212.
    [16] 朱新慧, 孙付平, 秦勇.  组合GPS和VLBI数据建立板块运动模型 . 武汉大学学报 ● 信息科学版, 2005, 30(7): 604-608.
    [17] 刘根友, 朱耀仲, 许厚泽, 张为民.  GPS监测中国沿海验潮站垂直运动观测研究 . 武汉大学学报 ● 信息科学版, 2005, 30(12): 1044-1047.
    [18] 吴云, 孙建中, 乔学军, 王辉.  GPS在现今地壳运动与地震监测中的初步应用 . 武汉大学学报 ● 信息科学版, 2003, 28(S1): 79-82,136.
    [19] 许才军, 刘经南, 晁定波, 李延兴.  利用 GPS复测资料研究华北地块旋转运动 . 武汉大学学报 ● 信息科学版, 2000, 25(1): 74-78.
    [20] 刘希康, 丁志峰, 李媛, 刘志广.  EMD在GNSS时间序列周期项处理中的应用 . 武汉大学学报 ● 信息科学版, 0, 0(0): 0-0. doi: 10.13203/j.whugis20210029
  • 加载中
图(3) / 表(2)
计量
  • 文章访问数:  960
  • HTML全文浏览量:  28
  • PDF下载量:  213
  • 被引次数: 0
出版历程
  • 收稿日期:  2014-05-29
  • 刊出日期:  2016-07-05

应用GPS连续观测修正流动观测垂向速率的分析与实验

doi: 10.13203/j.whugis20140251
    基金项目:

    中国地震局地震行业科研专项基金 No.201308009

    国家自然科学基金 No.41474002

    作者简介:

    占伟,博士生,研究方向为GPS数据处理与应用。zw000373@163.com

  • 中图分类号: P228.41

摘要: 目前,基于GPS流动观测确定垂向变化速率受年周期运动影响较大,但因观测数据较少,不能有效消除周期性运动的影响,使得GPS流动观测得到的垂向速率精度较低。提出利用同一区域GPS连续观测修正流动观测垂向速率的方法,首先依据GPS连续站垂向时间序列计算得到连续站的年周期运动振幅,然后进行空间插值得到各流动站的年周期运动振幅,进而对流动站垂向时间序列进行周期性运动修正,最后依据修正后的流动站垂向时间序列求取垂向线性速率。分别选用云南地区陆态网络GPS连续(2010~2013年)和流动观测数据(2011~2013年)进行了验证。结果表明,在年周期运动空间一致性较好的区域,该方法能够有效修正GPS流动观测垂向速率,提高垂向速率的解算精度。

English Abstract

占伟, 李斐, 朱爽, 张杰. 应用GPS连续观测修正流动观测垂向速率的分析与实验[J]. 武汉大学学报 ● 信息科学版, 2016, 41(7): 963-968. doi: 10.13203/j.whugis20140251
引用本文: 占伟, 李斐, 朱爽, 张杰. 应用GPS连续观测修正流动观测垂向速率的分析与实验[J]. 武汉大学学报 ● 信息科学版, 2016, 41(7): 963-968. doi: 10.13203/j.whugis20140251
ZHAN Wei, LI Fei, ZHU Shuang, ZHANG Jie. Analysis and Test of Correction to Vertical Velocity Measured by GPS Mobile Observation with Continuous Observation[J]. Geomatics and Information Science of Wuhan University, 2016, 41(7): 963-968. doi: 10.13203/j.whugis20140251
Citation: ZHAN Wei, LI Fei, ZHU Shuang, ZHANG Jie. Analysis and Test of Correction to Vertical Velocity Measured by GPS Mobile Observation with Continuous Observation[J]. Geomatics and Information Science of Wuhan University, 2016, 41(7): 963-968. doi: 10.13203/j.whugis20140251
  • 近年来,随着各类GPS监测网络的建立,尤其是中国地壳运动观测网络(简称“网络工程”)和中国大陆构造环境监测网络(简称“陆态网络”)的投入运行以来,产生了大量的GPS连续和流动观测数据。依据这些数据得到的GPS速度场及其衍生产品在地壳运动监测、地球动力学研究等领域发挥了重要作用[1-7]。在上述研究中,GPS流动水平速度场应用较多,而GPS流动垂向速度场受非构造运动影响较大且缺乏足够的观测数据来消除这些影响,使得GPS流动观测解算的垂向速度场精度不高,应用较少。

    目前,我国部分区域(如京津唐、川滇地区)GPS连续观测网络已初具规模,如何利用这些较高密度的连续观测网络提高同一区域的流动观测垂向速率精度,是挖掘GPS垂向速度场价值的当务之急,具有重要的意义。鉴于此,本文提出了利用GPS连续观测修正流动观测垂向速率的方法,并选用云南地区陆态网络GPS连续和流动观测数据进行了验证。

    • GPS垂向时间序列存在较为显著的季节性运动,且年周期运动最为突出[8-9],因此,本文只考虑年周期运动,则有:

      (1)

      式中,Y为垂向时间序列;a为常数项;V0为速率;AB为年周期运动正弦波和余弦波振幅;C为年周期运动总振幅;φ为年周期运动初相;ε为残差;t为时间。

      由于网络工程和陆态网络采用的流动观测模式为每隔2~3 a观测一期,每期连续观测4 d,这种观测模式难以获取准确的年周期运动信息,因此,求解时只能顾及线性运动:

      (2)

      在用最小二乘理论求解参数时,如果忽略部分参数,则估计参数的平差值会包含未估计参数的影响[8],即按式(2)求得的速率V1包含了年周期运动的影响,与V0存在一定偏差。如果能够获取每个流动站的年周期运动振幅并在时间序列Y中扣除年周期运动得到Y′,再根据Y′只顾及线性运动求解速率V2,则V2应该能够扣除年周期运动的影响,等于(或者接近于)V0

      (3)

      由于一些区域内GPS时间序列中的非构造噪声存在区域相关性[10-11],因此,可以利用同一个区域的连续站年周期运动信息通过空间加权内插得到流动站的年周期运动信息。对于每个流动站,年周期运动振幅A0B0分别取其周边200 km内所有连续站的年周期运动振幅AB的加权平均(权取距离d的平方倒数):

      (4)
    • 为了检核修正方法,选取陆态网络GPS连续站密度较高的云南地区进行试验。陆态网络在该地区布设了27个连续站,其中KMIN和XIAG站从1999年开始观测,其他25个测站是从2010年陆续开始观测。本文选取这27个测站2010-01-01以后的观测数据,至2013年底这些测站观测时间均超过960 d,能够较为准确地获取该地区垂向速率和年周期运动信息。利用GAMIT/GLOBK[12]和QOCA[8]软件解算各连续站观测数据,得到各连续站2010~2013年的时间序列,然后按式(1)分别求解各连续站垂向运动速率和年周期运动振幅(表 1图 1)。

      图  1  云南地区GPS连续站年周期运动 (2010~2013年,参考点2010.0)

      Figure 1.  Annual Terms of Vertical Position Variations for GPS Fiducial Stations(from 2010 to 2013, the Reference Time as 2010.0) in Yunnan Area

      表 1  云南地区GPS连续站垂向运动信息(2010~2013年)

      Table 1.  Vertical Motion for GPS Fiducial Station in Yunnan Area from 2010 to 2013

      点名速率 /(mm\5a-1)年周期观测天数
      振幅/mm初相/(°)
      KMIN2.939.4517.611 305
      XIAG0.179.5357.31 365
      YNCX*1.528.554.111 108
      YNDC*2.219.449.891 050
      YNGM-0.246.35350.08988
      YNHZ1.218.298.56973
      YNJD1.0410.978.311126
      YNJP-1.98.8817.381 134
      YNLA-0.6814.319.051 147
      YNLC*0.611.495.491 008
      YNLJ*1.3910.5618.841 115
      YNMH-1.0211.67.861 150
      YNMJ-0.39.890.151 129
      YNML1.07712.251 149
      YNMZ*-3.1710.6524.81 163
      YNRL-0.6411.810.571 087
      YNSD*2.1812.2711.65975
      YNSM*-1.5410.945.441 162
      YNTC0.8412.736.871 149
      YNTH*1.168.2717.311 164
      YNWS3.178.4721.88968
      YNXP2.157.956.631 215
      YNYA2.389.6814.571 006
      YNYL1.9512.8516.351 163
      YNYM*1.4510.2217.141 121
      YNYS2.1111.1112.931 160
      YNZD0.868.9722.98980
      注: 标*的点为检核站,未标注的点为基准站。

      表 1是云南地区GPS连续站2010~2013年垂向运动信息,从中可以看出该地区各GPS连续站年周期运动较为显著,年周期运动振幅较大(7~14 mm),垂向运动速率相对较小(4 mm/a以内),因此,测站年周期运动会对垂向运动速率的求解产生较大的影响,在该地区有必要对流动观测垂向运动进行修正。

      图 1是云南地区GPS连续站2010~2013年年周期运动,时间参考点为2010.0,箭头长度表示振幅,从正东逆时针旋转的方位角表示初相位方向,箭头指向东、南、西、北方向依次表示极值发生在0.25、0.50、0.75、0.0 年。结合表 1图 1可以看出云南地区GPS连续站的年周期运动空间一致性较好,各连续站年周期运动振幅和初相均比较一致,满足本方法的使用条件。

    • 在云南地区27个连续站中选取18个测站(图 1中蓝色点)记为基准站,9个测站(表 1标*点和图 1中黑色点)记为检核站,用以检核修正方法效果。首先根据每个连续站垂向时间序列Y按式(1)用最小二乘方法求解速率V0和年周期运动振幅A和B;然后根据每个检核站垂向时间序列Y按式(2)用最小二乘方法求解速率V1;再次,每个检核站按式(4)根据基准站的年周期运动振幅插值求得年周期运动振幅A0B0,并在时间序列Y中扣除依据A0B0计算得到的年周期运动,得到时间序列Y′,然后按式(3)用最小二乘方法求解速率V2。比较V1V2V0的差异(见表 2),如果V2V0的差异比V1V0的差异小,则说明修正方法有效。

      表 2  检核站解算结果

      Table 2.  Results of Check Sites

      检核站年周期运动振幅/mm速率结果/(mm\5a-1)
      ABA0B0V0V1V2V1-V0V2-V0
      YNCX8.530.619.551.661.52±0.122.04±0.181.55±0.120.510.03
      YNDC9.301.628.251.392.21±0.191.21±0.252.09±0.19-1.01-0.12
      YNLC11.441.108.950.230.6±0.151.21±0.260.71±0.160.610.11
      YNLJ9.993.4110.432.381.39±0.111.29±0.21.24±0.1-0.10-0.15
      YNMZ9.664.478.172.36-3.17±0.16-3.01±0.23-3.22±0.160.16-0.06
      YNSD12.022.4811.001.282.18±0.17-0.22±0.291.86±0.17-2.40-0.32
      YNSM10.891.0410.811.01-1.54±0.12-1.14±0.21-1.54±0.120.400.00
      YNTH7.902.468.091.631.16±0.111.34±0.171.1±0.120.18-0.05
      YNYM9.763.019.362.101.45±0.131.31±0.211.34±0.13-0.14-0.11

      表 2可以看出,不顾及年周期运动求得的速率(V1)与顾及年周期运动求得的速率(V0)差别较大,最大值达到-2.4 mm/a,说明在年周期运动较为显著的地区,在求解线性速率时必须考虑年周期运动。由于该区域年周期运动空间一致性较好,9个检核站内插得到的年周期运动振幅比较接近于实际值(差异在3 mm以内)。由于内插得到的检核站年周期运动振幅精度较高,使得本方法对垂向速率的修正效果较好,各检核站经过年周期运动修正后得到的速率(V2)比未经过修正的速率(V1)更为接近V0,9个检核站V2V0的差异量值在0.32 mm/a以内,且V2的误差比V1小,接近V0的误差。以上结果说明,在垂向时间序列年周期运动显著且空间一致性较好的区域,通过本方法能够较好地消除年周期运动对垂向速率的影响,提高垂向速率的解算精度。

    • 陆态网络在云南地区布设了157个GPS流动站,分别于2011年3~8月和2013年4~9月进行了一期观测,每期观测时段长度为4 d。利用GAMIT/GLOBK和QOCA软件解算了这两期流动观测数据(解算模式与§2中GPS连续站时间序列的解算模式保持一致,即GAMIT解算参数、QOCA解算中的框架点的选取等设置一致),得到了各流动站2011~2013年的时间序列,按式(2)求得了各流动站未扣除年周期运动的2011~2013年垂向速率(图 2(a))。利用§2计算得到的云南地区27个连续站年周期运动振幅,按式(4)内插得到各流动站年周期运动振幅,然后对各流动站的时间序列扣除年周期运动,再按式(3)求得扣除年周期运动的2011~2013年垂向速率(图 2(b))。

      图  2  云南地区GPS垂向速度场(流动观测,2011~2013年)

      Figure 2.  GPS Vertical Velocity Fields in Yunnan Area (Mobile Observations,from 2011 to 2013)

      图 2可知,未进行年周期改正的垂向速度场(图 2(a))和改正后的垂向速度场(图 2(b))整体趋势类似,即大体以北纬24°为界,北纬24°以北的区域以隆升为主,北纬24°以南以沉降为主,大体与由GPS连续站测定的云南地区垂向速度场(图 3)特征一致。改正后的垂向速度场比未改正的速度场在空间一致性上有所改善,且更加接近GPS连续站测定的垂向速度场和水准观测结果得到的该区地壳垂直运动特征[13-14],如图 2中两个红色椭圆框内区域。未改正的速度场显得较为混乱,方向上有升也有降,量值有大有小,改正后的速度场基本上都是上升,且量值空间分布更为均匀。图 2中绿色椭圆框区域的速度场经过改正后更接近郝明等的文章中图 7结果[14]。总体而言,本方法对GPS流动垂向速度场的改进取得了较好的效果。

      图  3  云南地区GPS垂向速度场(连续观测,2010~2013年)

      Figure 3.  GPS Vertical Velocity Fields in Yunnan Area (Continuous Observations,from 2010 to 2013)

    • 本文提出了利用连续观测站年周期运动信息通过空间插值获取流动观测站年周期运动信息、并在流动观测站时间序列中加以扣除、进而求解垂向速率的方法,并利用云南地区陆态网络2010~2013年连续观测数据和2011~2013年流动观测数据进行了验证,结果表明在年周期运动空间一致性较好的区域(云南地区)能够取得较好的效果。

      本方法的前提是区域内年周期运动空间一致性较好,只有满足这个条件,才能利用连续站年周期运动信息内插得到较为准确的流动站年周期运动振幅,才能得到较为准确的垂向速率。由于本文的主要目的是验证利用连续站信息修正流动观测垂向速率的方法有效性,因此,只考虑了季节性运动中最主要的年周期运动,但本方法也适用于其他季节性运动(如半年周期运动),只需在式(1)和式(3)增加相关参数求取并加以修正。

      我国国土辽阔,各区域情况不尽相同,可先处理区域内连续观测数据,分析区域垂向季节性运动机制。如果与趋势(线性)运动相比,季节性运动振幅较小,对垂向速率的影响可忽略不计,则无需进行季节性改正;如果区域内季节性运动影响较大,且有多种较为显著的周期运动(年周期、半年周期),则可以根据分析结果选取对垂向速率影响较大的季节项加以修正。为了进一步提高方法的准确性,还将进行空间插值方法的择优选取、在GPS连续站垂向季节性运动信息求取前扣除一些能够定量计算的非构造形变(如大气、积雪等负荷效应造成的非构造形变)[15-18])等分析研究。

参考文献 (18)

目录

    /

    返回文章
    返回