基于CORS站网的三峡库区负荷形变与重力场时空变化监测汇编.docx

基于CORS站网的三峡库区环境负荷形变与重力场时空变化监测章传银①王伟①甘卫军②李辉③张庆涛④①中国测绘科学研究院，100830，北京②中国地震局地质研究所，北京③中国地震局地震研究所，武汉④国家测绘地理信息局第一大地测量队，西安摘要：受地球动力学因素特别是地球表层大气、地表水及地下水动力环境影响，地面站点位置、地球重力场及大地水准面随时间变化。本文以三峡地区CORS站网为主，少量重力台站为辅，采用负荷形变与地球重力场严密组合方法，综合确定了2011年1月至2015年6月三峡地区环境负荷驱动的地壳形变与重力场月变化，结果显示：⑴CORS站网具备地壳垂直形变、大地水准面及地面重力变化监测能力；⑵三峡地区地壳垂直形变年变化幅度36.1mm，大地水准面年变化幅度28.2mm，地面重力年变化幅度117.4μGal；⑶CORS站网地面重力变化监测精度水平不低于流动重力场重复测量；CORS站网地壳垂直形变与地面重力变化监测具有一定的外推预报能力。关键词：GNSS CORS站；环境负荷形变；区域重力场变化一、引言时变地球重力场包含地球系统密度分布及物质运动的丰富信息，直接反映了地球各圈层最基本的物质及其变化特性，是地球环境变化导致地球形变最基本、最直接和最重要的物理量之一（许厚泽，2010）。受引潮天体、海平面及大气变化、地表及地下水变化（包括冰川雪山、江河湖库、土壤水和地下水等变化，也称大陆水变化）、核幔耦合及内核动力作用等地球动力学因素影响，地面站点位置、地球重力场及大地水准面随时间变化。为方便高精度大地测量数据处理及定量分析，通常将地球重力场随时间变化分成两个部分：一部分是按大地测量参考系统约定的动力学模型，计算的固体潮、海潮负荷及大气潮负荷的影响，称为潮汐影响；另一部分是扣除潮汐影响后的残差重力场时变，称为地球重力场形变。由于潮汐影响能随时随地从大地测量观测量和参数中予以移去或恢复，因此，时变地球重力场监测和研究的重点通常是地球重力场形变（党亚民等，2015）。在大地测量基准维持涉及的数年到十年时间内，由非潮汐海平面变化、大陆水变化和非潮汐大气负荷变化等构成地球表层环境负荷变化，对重力场变化及地壳非构造垂直形变的贡献占其时变量的97%以上。正因为如此，卫星重力场时变监测也极大地推动了地球水循环与海洋、大气、水资源环境及相关地球科学的发展。本文以三峡地区CORS站网非潮汐垂直形变为主，少量重力台站非潮汐重力变化为辅，综合确定三峡地区环境负荷驱动的地壳形变与重力场时空变化。二、监测原理与方法地球表层环境负荷变化可用等效水高变化（或面密度变化）表示。这里的等效水高变化通常是指某一时刻等效水高相对于参考时刻等效水高的差值，或相对于某一段时间等效水高平均值之差。等效水高变化有时也称为等效水高异常。地面等效水高变化的规格化负荷球谐展开式为：⑴由地球负荷形变和地球重力场理论，地面等效水高变化引起的大地水准面变化为：⑵引起的地面重力变化为：⑶引起的地面站点大地高变化（负荷垂直形变）为：⑷式⑴~⑷中，为阶次规格化负荷球谐系数；为阶次规格化连带勒让德函数；为地面点的地心纬度和经度；为阶位负荷（勒夫）数；为阶径向负荷数；为水的密度；为固体地球平均密度；为万有引力常数；为地球总质量；为地球平均半径；为地面平均重力。地面等效水高变化引起的地面站点位置或地面重力场参数变化还可采用负荷格林函数积分方法表示为如下通用形式：⑸式中，为地面积分流动点；为地面计算点；为流动点到计算点的空间距离；为球面角距；为地球表面；为通用的负荷格林函数，当为地面重力变化时，取重力格林函数：⑹式中，为阶勒让德函数。当为地面大地高变化（即径向变化）时，取径向格林函数：⑺同理，可列出环境负荷变化引起的垂线偏差、负荷水平形变、地面扰动重力变化球谐展开式和负荷格林函数式。由⑵~⑸式不难看出，通过构建以卫星重力、卫星测高或各种地面大地测量变化量为观测量，以负荷等效水高或其球谐系数为待估参数的观测方程，采用适当的参数估计算法，可反演地壳负荷形变及地面重力场变化。将⑺代入⑸式，可将CORS站大地高变化表示为：⑻式中：为构造垂直形变与地下水均衡垂直形变。三峡地区构造垂直形变很小，一年的地下水变化基本处于平衡状态，地下水均衡效应很不明显，本文将作为随机噪声处理（这种处理方式不适合活动构造边缘、地下水趋势性变化明显的CORS站）。同理，可将重力台站的重力变化表示为⑼三、数据处理与综合反演⑴采用的数据资源及预处理①三峡地区26座CORS站非潮汐大地高变化。CORS网解算时的地球动力学改正采用IERS2010协议，移去了固体潮、海潮负荷、周日和半日大气潮负荷影响。以2010年11月平均值为基准，生成2011年1月至2015年6月CORS站大地高月变化②三峡地区8座重力台站非潮汐重力变化。重力台站都配置便携式潮