利用GPS与环境负荷形变数据研究台风引起的垂向地表位移

台风造成的强降雨、低气压、海面高度变化均会引起地表的形变.本文利用中国大陆构造环境监测网(陆态网)7个GPS台站每日的垂向位移和环境负荷形变模型分析2018年9月10—26日台风"山竹"期间不同负荷引起的区域垂向地表形变.结果表明,台风期间大气负荷和非潮汐海洋负荷垂向形变最大分别达到5.1 mm和-9.2 mm.模型能较好地反映河流区域地表水文负荷变化造成的垂向形变,但不同模型之间存在系统偏差.由于缺少地下水等信息,模型反映负荷长期形变效应的效果不佳,且形变的量级明显小于GPS观测的结果.迅速增加的水文负荷使北海GPS站从开始下沉到最低点(-15.6 mm)5天的下沉量达到25.7 mm;珠海、广州GPS站均观测到河流汇水作用造成地表的二次下沉,且珠海站一周后才抬升到正常位置;湛江和北海GPS站能较好地反映河流水位变化,相关系数分别为-0.66和-0.50.研究结果表明,相比于形变模型,GPS能更有效地监测台风短期水文负荷形变,可为台风洪水等灾害监测与预报提供有用的信息.     

全球水质量迁移对海平面空间模式周年变化的影响

大气、陆地水和海洋之间的水质量迁移对海平面的影响一般假定为均匀薄层分布.但实际上,水质量负荷重新分布一方面会使地壳产生形变,另一方面会引起重力位势场变化(引力位和离心力位),这都会对海平面时空变化特征产生影响,两方面之和称为负荷自吸引效应(SAL).海洋模式模拟的时变洋底压力结果一般符合Boussinesq假设即体积守恒,忽略了大气、陆地水和海洋之间水质量交换的影响.本文基于海平面变化方程,联合陆地水模型、大气地表气压模型、海洋洋底压力模型和GRACE反演的冰川质量变化,详细讨论了2003-2010年SAL对海平面周年变化的影响.主要结论有:(1)SAL对全球海平面周年变化有显著影响,振幅在1.3~19 mm.其中近海岸和低纬度区域受影响较大.(2)在SAL引起的海平面周年振幅变化的因素中,陆地水储量变化因素最大,大气因素次之,非潮汐海洋影响最小.但非潮汐海洋对海平面周年相位空间变化的影响最为复杂.(3)通过与国际长期验潮站观测数据结果比较,在ECCO海洋模式估计的洋底压力结果中引入SAL,能多解释约5.3%观测信号方差.     

海潮负荷对沿海地区宽幅InSAR形变监测的影响

海岸带地区是全球自然生态环境最为复杂和脆弱的地域之一,合成孔径雷达干涉测量（InSAR）技术可以为全球人类活动、气候变暖和俯冲带剧烈构造运动等背景下的大范围海岸带地理环境变化研究提供重要观测资料.海洋潮汐导致固体地球长周期形变,波长尺度为10²~10³ km的海潮负荷引入mm级至cm级的形变梯度,此类非构造信号对海岸带InSAR精密形变分析（如：大范围、微小、缓慢且非稳态构造过程等）造成显著影响.本文以宽幅模式SAR数据为例,基于多种海潮模型研究了全球典型海岸带地区（福建、智利和阿拉斯加湾）海潮负荷效应对宽幅InSAR形变监测的影响,给出了宽幅InSAR海潮负荷三维分量估计与差分相位提取方法,并进一步讨论了基于不同海潮模型估计海潮负荷位移的差异.海潮负荷影响不仅与研究范围大小有关,其形变梯度变化与研究区域地形特征存在强相关,对于长波长形变分析而言,传统平面或者曲面拟合方法难以有效分离海潮负荷位移.     

近海潮汐效应对测站位移的负荷影响

利用TPXO6全球海潮模型和中国东海和南海潮汐资料研究了近海潮汐效应对我国沿海部分台站负荷位移的影响.结果表明,近海效应对位移的水平和垂直分量的影响均达到了mm的量级,其中在厦门站的影响最大,水平和垂直分量的振幅分别接近3 mm和15 mm,因此近海潮汐效应是沿海地区GPS动力学应用中应该考虑的负荷效应中的一个因素.另外,各台站负荷位移水平分量随时间变化的规律是不一样的,最大值对于不同的台站出现在不同的方位,并且对于垂直分量,在上海站和在其他沿海站的影响存在着明显的反相关系.本文结果可为GPS资料处理提供有益参考.     

非构造形变对GPS连续站位置时间序列的影响和修正

GPS观测得到的地壳形变场通常包含有构造形变与非构造形变二类信息, 去除其中的非构造形变信息对于有效运用GPS数据研究构造形变场至关重要. 本文运用国际卫星对地观测资料及各类地球物理模型, 定量计算海潮、大气、积雪和土壤水、海洋非潮汐4项负荷效应造成的地壳非构造形变, 并以此研究和修正这些非构造形变对中国地壳运动观测网络GPS基准站位置时间序列的影响. 研究发现此4项负荷效应, 特别是大气、积雪和土壤水, 对于测站垂向位置的影响显著. 通过模型改正可以使测站垂向位置的RMS降低～1 mm, 占其总量的～11％. 对于垂向时间序列的周年项部分, 这一改正可降低其振幅的37％. 研究还表明经过地球物理模型改正和周年、半周年谐波拟合改正的时间序列比起仅经过周年、半周年谐波拟合改正的时间序列更为平滑, 表明地球物理模型改正对于消除非构造形变场的作用不是周年、半周年谐波拟合改正所能替代的.     

海南井水位对热带气旋响应特征的分析

大量观测资料及实验研究证实,地下水对地壳应力-应变的响应频率较宽,可记录到地震波中的瑞利波(10～20 s)、潮汐波动(1/3～1/2 d)、气压波动(儿小时至几十天)等.在文献中,已见热带气旋引起地震仪的震颤波,倾斜仪和重力仪的“脉动加粗”等方面研究成果,但井水位对大气环流、雷暴高压、热带气旋等的响应,则未见报导.     

香港GPS基准站坐标序列特征分析

利用香港GPS连续运行参考站网络2001年1月至2007年8月的观测资料,全面深入地分析了12个基准站坐标序列特征.本文采用主成分空间滤波算法去除公共误差,来提高坐标序列的信噪比,并采用最大似然估计准则定量估计滤波后坐标序列的噪声特性,计算了地球表面质量负荷（包括大气、非潮汐海洋、积雪和土壤水）对香港GPS基准站坐标序列的影响.研究结果表明：香港GPS基准站坐标序列具有高度的空间相关性,其公共误差具有较强的季节性变化特征;地表质量负荷变化引起的香港地壳形变可以解释公共误差序列中约为3mm的垂向周年变化,经过质量负荷改正后的公共误差序列与高阶电离层误差高度相关;滤波后坐标序列的噪声特性可以用可变白噪声加闪烁噪声模型来描述,顾及闪烁噪声所计算的速度误差要比只考虑可变白噪声计算的速度误差大2~6倍;基准站间存在达1.5 mm/yr的相对水平运动,揭示香港地区存在活动断层;部分基准站坐标具有明显的振幅为1～2 mm本地季节性变化,所有测站的残差序列也表现出强烈的季节性变化.     

大气负载引起地表位移:对重力和基线测量的影响

大气质量对地壳加载并使之形变.通过将逐日的全球气压数据和质量负载的格林函数进行摺积求和,我们怡算大气负载(包括那些与短期大气扰动有关的大气负载)对地表刚点定位测量和地表重力观测的时间依赖关系.我们茸出了无海洋地球和具有一负压海洋(inverted barometer ocean)的地球的响应.负压海洋显著地影响近海岸台站的负载响应.峰峰垂直位移常常是15一20mm,伴随有3一6卜Gal的扰动.基线变化可能大到20mm或更大.这种扰动在高纬度区和冬季达到最大.这些振幅的变化与Rabbel和Zs。hau(1985)的结果一致.他们用距某一点半径有关的高斯函数模拟了气压扰动.使用台距约I00okm以内测点的现有气压数据就足以计算地壳的形变.单是局部气压数据是不够的.Rabbel和Zsc五-au为使用局部气压和区域平均气压算得的这些位移位所假设的改正值,对正好在内地的测点证明是精确的,但一般地说在距海岸几百公里之内,改正值是不适当的.     

SNREI地球对表面负荷和引潮力的形变响应

基于PREM模型，利用非自转、球型分层、各向同性、理想弹性（SNREI）地球的形变理论，讨论了地球在不同驱动力作用下的形变特征．采用地球位移场方程的4阶Runge Kutta数值积分方法，解算了在表面负荷和日月引潮力作用下地球表面和内部形变和扰动位，并给出了地球表面的负荷Love数和体潮Love数．结果表明在固体内核中的形变很小，液核中低阶(n＜10)负荷位移随半径的变化非常复杂．当负荷阶数超过10时，地核中的形变和扰动位都很小，地球的响应主要表现为弹性地幔中的径向位移，且随深度增加急剧减弱，负荷阶数越高这种衰减的速度越快．SNREI地球的地表负荷Love数和体潮Love数与信号频率的依赖关系很弱．在计算体潮Love数的过程中，采用了SNREI地球的运动方程，同时考虑了由于地球自转和椭率引起的核幔边界附加压力，这一近似处理方法获得的结果能很好地符合地球表面重力潮汐实际观测结果．     

环境负载对云南地区基准站时间序列的噪声模型特性研究

以云南地区27个GPS基准站坐标时间序列为研究对象,使用赤池信息量和贝叶斯信息量估计准则（AIC/BIC）对其进行噪声分析,计算并扣除时间序列中大气负载、非潮汐海洋负载和水文负载3种环境负载引起的位移量,得到各基准站分量在环境负载改正前后的最优噪声模型。结果表明,部分基准站分量经过负载改正后最优噪声模型会发生变化,改正前后的大部分基准站噪声特性均体现为闪烁+白噪声和幂律噪声。环境负载对基准站的垂向位移影响比水平向更为显著,水文负载成为影响基准站的最大因素,最大位移量达到厘米级。分析环境负载改正前后噪声特性的变化表明,环境负载改正在U分量上的影响最大,N分量次之,E分量最小,噪声模型的变化在地域上并未呈现明显规律。     

An estimate of the influence of loading effects on tectonic velocities in the Pyrenees

Surface displacements due to temporal changes in environmental mass redistributions are observable in the coordinate time series of many Global Navigation Satellite System (GNSS) sites. In this study, we investigated the effect of loading on estimates of tectonic velocity computed from campaign-style GNSS observations. The study region is in the Pyrenees mountain range between France and Spain (ResPyr campaigns). In this area, seismic activity is continuous and moderate and the expected amplitude of the horizontal tectonic velocity is less than 0.5 mm/yr. In order to determine the velocity, 4 sparse GNSS campaigns were carried out from 1995 to 2010. Considering this small rate of deformation, loading phenomena can contribute a non-negligible artifact to the velocity computation that could affect our geodynamical interpretation. In this investigation, we specifically considered the atmospheric, hydrological, and non-tidal ocean loading phenomena. The computed loading deformations for this region show the horizontal displacements are dominated by the non-tidal ocean loading (maximum 4 mm for the North and 3.1 mm for the East components); the main vertical contributions are due to the atmospheric and continental water storage loading (maximum 14.3 for the atmosphere and 8.1 mm for the hydrology, respectively). We have found that the dominant loading effect on the horizontal velocity is the non-tidal ocean loading (mean of 0.11 mm/yr), whereas the vertical component is dominated by the hydrological loading (mean of 0.21 mm/yr). Since the study area is in a mountainous region, we also analyzed the difference between the atmospheric and the topography dependent atmospheric loading models at our GNSS campaign sites. We did not find any significant difference between the two atmospheric loading models in terms of horizontal velocity. Finally, we performed simulations to identify the optimum timing and frequency of future GNSS campaigns in this area that would minimize the loading effects on tectonic velocity estimates.     

Determination of vertical displacements over the coastal area of Korea due to the ocean tide loading using GPS observations

This paper describes the GPS applicability for detecting the vertical displacements of ground stations caused by ocean tide loading effects. An experiment was carried out using 12 permanent GPS stations located in the coastal area of Korea using data in the period 1 July until 26 August 2003. The relative height differences were calculated from hourly DGPS data processing based on the carrier-phase observation. The power spectra of the M₂ and N₂ constituents of ocean tide loading were derived using the CLEAN algorithm. The differential vertical displacements generated by the ocean tide loading effect are typically 3–25 mm in coastal area of the Korea. We compared the results from GPS with those of the ocean tide models, NAO.99Jb regional model and GOT00.2, FES99 global models. The M₂ (N₂) amplitude differences of vertical displacements between GPS and GOT00.2 is 1.22 ± 3.61 mm (1.01 ± 1.48 mm), and that of the M₂ (N₂) amplitude difference between GPS and FES99 is 0.04 ± 4.64 mm (0.64 ± 1.75 mm), whereas the M₂ (N₂) amplitude difference between GPS and NAO.99Jb are 0.05 ± 1.03 mm (0.86 ± 1.18 mm). The highest vertical displacements at the PALM station are found for 24.5 ± 0.7 mm from GPS observation, and 22.9 mm from the regional model NAO.99Jb and 13.17 and 10.00 mm from the global models GOT00.2 and FES99, respectively. These values show that the vertical displacements derived from GPS are in good agreement with those of the regional model NAO.99Jb around Korea, more than with the global models. This result indicated that GPS is an effective tool to measure the vertical displacement caused by the ocean tide loading effect in the coastal area, and we need to use the NAO.99Jb ocean tide model rather than the global ocean tide models in and around the Korean peninsula for position determination with permanent GPS installations. This work demonstrates that vertical displacement caused by the M₂ and N₂ constituents of ocean tide loading can be measured by carrier-phase DGPS.     

Annual variation detected by GPS,GRACE and loading models

Most GPS coordinate time series, surface displacements derived from the Gravity Recovery and Climate Experiment (GRACE), and loading models display significant annual signals at many regions. This paper compares the annual signals of the GPS position time series from the Crustal Dynamics Data Information System (CDDIS), estimates of loading from GRACE monthly gravity field models calculated by three processing centers (Center of Spatial Research, CSR; Jet Propulsion Laboratory, JPL; GeoForschungsZentrum, GFZ) and three geophysical fluids models (National Center for Environmental Prediction, NCEP; Estimating the Circulation and Climate of the Ocean, ECCO; Global Land Data Assimilation System, GLDAS) for 270 globally distributed stations for the period 2003-2011. The results show that annual variations derived from the level-2 products from the three GRACE product centers are very similar. The absolute difference in annual amplitude between any two centers is never larger than 1.25 mm in the vertical and 0.11 mm in horizontal displacement. The mean phase differences of the GRACE results are less than ten days for all three components. When we correct the GPS vertical coordinate time series using the GRACE annual amplitudes using the products from three GRACE analysis centers, we find that we are able to reduce the GPS annual signal in the vertical at about 80% stations and the average reduction is about 47%. In the north and the east, the annual amplitude is reduced on 77% and 72% of the stations with the average reduction 32% and 33%. We also compare the annual surface displacement signal derived from two environmental models; the two models use the same atmospheric and non-tidal ocean loading and differ only in the continental water storage model that we use, either NCEP or GLDAS. We find that the model containing the GLDAS continental water storage is able to better reduce the annual signal in the GPS coordinate time series.     

Crustal loading in vertical GPS time series in Fennoscandia

We compare time series of vertical position from GPS with modelled vertical deformation caused by variation in continental water storage, variation in the level of the Baltic Sea, and variation in atmospheric pressure. Monthly time series are used. The effect of continental water storage was calculated from three different global models. The effect of non-tidal variation in Baltic Sea level was calculated using tide gauge observations along the coasts. Atmospheric loading was computed from a numerical weather model. The loading time series are then compared with three different GPS time series at seven stations in Fennoscandia. A more detailed analysis is computed at three coastal stations. When the monthly GPS time series are corrected using the load models, their root-mean-square scatter shows an improvement between 40 and 0%, depending on the site and on the GPS solution. The modelled load effect shows a markedly seasonal pattern of 15 mm peak-to-peak, of which the uncorrected GPS time series reproduce between 60 and 0%.     

中国大陆精密重力潮汐改正模型

利用理论和实验重力固体潮模型,充分考虑全球海潮和中国近海潮汐的负荷效应,建立了中国大陆的精密重力潮汐改正模型.结果表明,采用不同的固体潮模型会对重力潮汐结果产生相对变化幅度小于0.06％的差异;在沿海地区海潮负荷的影响约为整个潮汐的4％,而中部地区约为1％,其中中国近海潮汐模型的影响约占整个海潮负荷的10%,内插或外推潮波的负荷约占海潮负荷的3％.通过比较实测的重力数据表明,本文给出的重力潮汐改正模型的精度远远优于0.5×10^-8 m·s^-2,说明了本文构建的模型的实用性,可为中国大陆高精度重力测量提供有效参考和精密的改正模型.     

武汉台重力潮汐长期观测结果

采用武汉台超导重力仪(SG C032)14年多的长期连续观测资料,研究了固体地球对二阶和三阶引潮力的响应特征,精密测定了重力潮汐参数,系统研究了最新的固体潮模型和海潮模型在中国大陆的有效性.采用最新的8个全球海潮模型计算了海潮负荷效应,从武汉台SG C032的观测中成功分离出63个2阶潮汐波群和15个3阶潮汐波群信号,3阶潮波涵盖了周日、半日和1/3日三个频段.重力潮汐观测的精度非常高,标准偏差达到1.116 nm·s^-2,系统反映了非流体静力平衡、非弹性地球对2阶和3阶引潮力的响应特征.结果表明,现有的武汉国际重力潮汐基准在半日频段非常精确,但在周日频段存在比较明显的偏差,需要进一步精化.对于中国大陆的大地测量观测,固体潮可以采用Dehant等考虑地球内部介质非弹性和非流体静力平衡建立的固体潮理论模型或Xu 等基于全球SG观测建立的重力潮汐全球实验模型作为参考和改正模型,海潮负荷效应应该采用Nao99作为改正模型.     

Modelling atmospheric and induced non-tidal oceanic loading contributions to surface gravity and tilt measurements

We investigate the contribution of atmospheric and its induced non-tidal oceanic loading effects on surface time-varying gravity and tilt measurements for several stations in Western Europe. The ocean response to pressure forcing can be modelled accordingly to the inverted barometer, i.e. assuming that air pressure variations are fully compensated by static sea height changes, or using ocean general circulation models. We validate two runs of the HUGO-m barotropic ocean model by comparing predicted sea surface height variations with hundred tide-gauge measurements along the European coasts. We then show that global surface pressure field, as well as a barotropic high-resolution ocean model forced by air pressure and winds allow in most cases a significant reduction of the variance of gravity residuals and, to a smaller extends tilt residuals.We finally show that precise gravity measurements with superconducting gravimeters allow the observation of large storm surges, occurring in the North Sea, even for inland stations. However, we also confirm that the continental hydrology contribution cannot be neglected. Thanks to their specific sensitivity feature, only tiltmeters closest to the coast can clearly detect the loading due to these storm surges.