Nontectonic signals in reprocessed continuous GPS position time series of CMONOC

The nearly nine-year continuous GPS data collected since 1 March 1999 from the Crustal Motion Observation Network of China(CMONOC) were consistently analyzed.Most of the nonlinear movements in the cumulative position time series pro-duced by CMONOC data center disappeared;and more accurate vertical terms and tectonic signals were extracted.Displacements caused by atmospheric pressure loading,nontidal ocean loading,soil moisture mass loading,and snow cover mass loading using the National Centers for Environmental Prediction(NCEP) Reanalysis I/II models and Estimation of the Circulation and Climate of the Ocean(ECCO) data can explain most of the vertical annual terms at many stations,while only parts can be explained at Lhasa and southern coastal sites,indicating that there are some deformation mechanisms that are still unknown or not modeled accurately.The remarkable differences in vertical position time series for short-baseline sites reveal that GPS stations can be greatly affected by lo-cal factors;and attention should be paid when explaining observed GPS velocity vectors.     

Assessment of second- and third-order ionospheric effects on regional networks: case study in China with longer CMONOC GPS coordinate time series

Higher-order ionospheric (HOI) delays are one of the principal technique-specific error sources in precise global positioning system analysis and have been proposed to become a standard part of precise GPS data processing. In this research, we apply HOI delay corrections to the Crustal Movement Observation Network of China’s (CMONOC) data processing (from January 2000 to December 2013) and furnish quantitative results for the effects of HOI on CMONOC coordinate time series. The results for both a regional reference frame and global reference frame are analyzed and compared to clarify the HOI effects on the CMONOC network. We find that HOI corrections can effectively reduce the semi-annual signals in the northern and vertical components. For sites with lower semi-annual amplitudes, the average decrease in magnitude can reach 30 and 10 % for the northern and vertical components, respectively. The noise amplitudes with HOI corrections and those without HOI corrections are not significantly different. Generally, the HOI effects on CMONOC networks in a global reference frame are less obvious than the results in the regional reference frame, probably because the HOI-induced errors are smaller in comparison to the higher noise levels seen when using a global reference frame. Furthermore, we investigate the combined contributions of environmental loading and HOI effects on the CMONOC stations. The largest loading effects on the vertical displacement are found in the mid- to high-latitude areas. The weighted root mean square differences between the corrected and original weekly GPS height time series of the loading model indicate that the mass loading adequately reduced the scatter on the CMONOC height time series, whereas the results in the global reference frame showed better agreements between the GPS coordinate time series and the environmental loading. When combining the effects of environmental loading and HOI corrections, the results with the HOI corrections reduced the scatter on the observed GPS height coordinates better than the height when estimated without HOI corrections, and the combined solutions in the regional reference frame indicate more preferred improvements. Therefore, regional reference frames are recommended to investigate the HOI effects on regional networks.     

陆态网基准站高程时间序列的环境负载效应分析

GPS坐标时间序列呈现显著的季节性变化,以中国大陆构造环境监测网络(CMONOC) 233个GPS测站的高程方向坐标时间序列为研究对象,通过计算环境负载效应改正前后的均方根(RMS)减少率,对比了德国波茨坦地学研究中心(GFZ)与法国斯特拉斯堡大学提供的EOST不同环境负载产品对CMONOC坐标时间序列数据的改善程度.统计结果表明：经过GFZ环境负载产品和EOST环境负载产品改正后分别有90.2%和75.54%的测站RMS减小,说明GFZ的环境负载产品对CMONOC基准站坐标时间序列数据的改善效果更好.同时发现环境负载效应不仅对GPS坐标时间序列中的季节项有明显的改正效果,对非季节项的改善效果也非常明显,经过GFZ提供的环境负载产品改正后有87.8%的GPS测站周年项振幅削弱,有82.4%的测站非季节项的RMS减小.     

环境负载对区域GPS基准站时间序列的影响分析

基于1999-2011年中国地壳运动观测网络(CMONOC)及2007-2012年武汉市连续运行卫星定位服务系统(WHCORS)的观测数据,联合中国及周边35个IGS基准站数据,采用GAMIT软件解算,获得了ITRF2008框架下的GPS基准站坐标时间序列.然后采用QOCA计算了环境负载位移,并利用其对GPS时间序列进行改正.研究结果表明环境负载造成的中国区域基准站位移呈现显著的区域性特征,东北、华北、华中区域变化较为一致且较大,西南区域基准站(KUNM)垂直方向的负载位移均方根(RMS)最大值达6.09 mm.通过对改正前后的数据进行时间序列分析,认为环境负载改正能够削弱中国区域大多数GPS基准站(约70%)垂直及水平东方向位移时间序列的非线性变化,其加权均方根(WRMS)减小量最大达1.5 mm,但是对于水平北方向坐标时间序列的改善作用不明显.     

中国大陆GPS基准站的时间序列特征

利用1999年3月至2002年3月期间“中国地壳运动观测网络”25个GPS基准站与周边国家IGS站的观测资料，进行了时间序列分析研究，并对提高GPS基准站的观测精度进行了探讨。结果表明，GPS基准站的时间序列具有一定的周期性，高程分量的周期性最为明显，GPS点位高程时间序列拟舍曲线的波峰和波谷出现的时间表现为区域性的不同。     

“陆态网”GPS与GRACE的中国大陆地表垂直形变对比分析

基于中国大陆构造环境监测网络的连续GPS观测数据,比较分析了中国大陆234个GPS台站和GRACE得到的地表垂直形变。GPS和GRACE垂直形变具有较好的一致性,反映了地表质量变化是引起GPS垂直形变非线性变化的重要因素之一,但二者也存在着一定的差异。为定量分析GPS和GRACE垂直形变的差异,探讨了热膨胀效应对GPS垂直位移的影响及区域地壳结构对GRACE估算地表垂直负荷形变的影响。结果表明,中国大陆50%以上的GPS台站热膨胀垂直形变周年振幅不小于1 mm;对GPS进行热膨胀效应改正后,中国大陆GPS与GRACE垂直形变具有更好的一致性;GPS与GRACE垂直形变周年振幅比值由1.07±0.06变为1.01±0.05;热膨胀效应可以解释6.2%的GPS与GRACE垂直形变的差异,热膨胀效应改正可使GPS和GRACE垂直形变的一致性相对增加11.2%。是否顾及区域地壳结构引起的GRACE估算中国大陆垂直负荷形变的相对差异为2.5%。     

利用GPS垂向位移监测西南地区干旱事件

GPS垂向位移包含了陆地水储量变化引起的地球弹性形变信息。本文利用中国大陆构造环境监测网(陆态网)31个GPS台站2010年8月—2016年12月的垂向位移数据,分析我国西南地区的干旱事件。结果表明,GPS垂向位移的季节性变化与降雨、GRACE反演的水储量、河流水位变化具有较强的负相关关系,其异常变化(非季节性变化)与我国使用的综合气象干旱指数(CI)、GRACE干旱强度指数(GRACE-DSI)相关性较强(相关系数约为-0.70),而与标准化降水蒸散指数(SPEI)的相关性较弱。相比于CI干旱指数,GPS垂向位移能更好地探测短期降雨急剧减少引起的干旱事件。数据缺失影响了GRACE-DSI探测干旱的准确性,而SPEI的波动较大,且存在高估某些短期降雨严重性的情况。研究结果展示了GPS在干旱监测方面的可靠性和优势。     

北斗二代导航卫星系统地壳运动监测能力

为了评估北斗卫星导航系统（BeiDou navigation satellite system, BDS）监测中国大陆地区地壳变形的技术能力,利用GAMIT/GLOBK软件处理了2017—2019年中国大陆构造环境监测网络23个基准站的全球定位系统（global positioning system, GPS）与BDS-2双模观测数据。结果显示,北斗二代的水平和垂向单日测站定位精度分别约为5~7 mm和13 mm,基线相对定位精度水平分量达到3~4 mm+（1~2）×10-8,水平位移速度测定精度约为0.6 mm/a。北斗二代的精密定位水平大致与20世纪90年代初GPS相当,可用于测定大尺度的板块运动及板内变形,但受卫星星座和定轨精度限制,不能准确反映季节性变动状态。作为对现有GPS监测的补充,可将基准站3年尺度的地壳运动监测精度最多提高20%。     

利用GPS垂直位移反演云南省陆地水储量变化

地表质量的重分布会引起固体地球的弹性形变,GPS连续运行观测站能够精确测定地表负荷引起的地壳形变。本文通过模拟数据对利用云南省及其周边47个中国大陆构造环境监测网（陆态网）台站反演云南地区陆地水储量的可行性进行分析：以水文模型周年振幅为真值,计算47个台站点的负荷形变,同时加入随机误差构成模拟观测数据,最后采用模型反演陆地水储量变化;1000次的随机模拟试验表明利用当前GPS台站数据可有效地反演云南地区陆地水储量变化。基于上述结论,笔者反演了云南省2010—2014年陆地水储量变化,GPS反演结果表明：云南省陆地水变化呈现明显的地域分布特征,西南部高山地区的水储量周年变化高于东部平原地区;在时间尺度上,云南省大部分地区水储量在10月（夏季末）达到最大值,在4月（冬季末）达到最小值;云南省2010—2014年陆地水呈缓慢增长趋势,约为20 mm/a。通过GPS陆地水储量反演结果与GRACE、GLDAS以及TRMM数据综合对比分析,表明利用云南地区当前GPS台站可以作为独立观测量用于GRACE与GRACE Follow-on衔接期间的陆地水储量变化监测。     

Spatiotemporal filtering for regional GPS network in China using independent component analysis

Removal of the common mode error (CME) is a routine procedure in postprocessing regional GPS network observations, which is commonly performed using principal component analysis (PCA). PCA decomposes a network time series into a group of modes, where each mode comprises a common temporal function and corresponding spatial response based on second-order statistics (variance and covariance). However, the probability distribution function of a GPS time series is non-Gaussian; therefore, the largest variances do not correspond to the meaningful axes, and the PCA-derived components may not have an obvious physical meaning. In this study, the CME was assumed statistically independent of other errors, and it was extracted using independent component analysis (ICA), which involves higher-order statistics. First, the ICA performance was tested using a simulated example and compared with PCA and stacking methods. The existence of strong local effects on some stations causes significant large spatial responses and, therefore, a strategy based on median and interquartile range statistics was proposed to identify abnormal sites. After discarding abnormal sites, two indices based on the analysis of the spatial responses of all sites in each independent component (east, north, and vertical) were used to define the CME quantitatively. Continuous GPS coordinate time series spanning \(\sim \)4.5 years obtained from 259 stations of the Tectonic and Environmental Observation Network of Mainland China (CMONOC II) were analyzed using both PCA and ICA methods and their results compared. The results suggest that PCA is susceptible to deriving an artificial spatial structure, whereas ICA separates the CME from other errors reliably. Our results demonstrate that the spatial characteristics of the CME for CMONOC II are not uniform for the east, north, and vertical components, but have an obvious north–south or east–west distribution. After discarding 84 abnormal sites and performing spatiotemporal filtering using ICA, an average reduction in scatter of 6.3% was achieved for all three components.     

地表质量负载对格陵兰岛区域GPS站点坐标时序特征的影响

选取ITRF2008框架下格陵兰岛区域12个GPS站2013年1月-2016年12月期间的日解坐标时间序列作为研究对象,并利用极大似然估计分析地表质量负载改正前后各站点的噪声特性、速度场及周期项振幅。结果表明：站点最优噪声模型主要为白噪声+幂律噪声与白噪声+闪烁噪声,地表质量负载形变修正GPS坐标时序后,明显增加U方向闪烁噪声的成分,平均降低其速度约0.36 mm/a,对水平方向影响较小;同时分别降低高程方向44.1%、14.2%的1 a项、0.5 a项振幅,相反,却增加了水平方向的周期项振幅。     

GPS坐标时间序列中的异常高频周期性噪声

本文分析了中国地壳运动观测网络GPS基准站及中国周边IGS连续站坐标中的周期性噪声,在北、东和垂向分量坐标时间序列中均发现了周期约351/n(n=1,…,6)天的正弦项,而地表质量负荷造成的位移序列中并没有对应的成分,不能解释异常周期项的来源。在估计台站的运动速率时,如考虑此类"异常"周期项,速率误差略微减小,而周年项振幅的估计误差显著增大,且残差离散度未见明显改善。因此,不推荐在GPS速率矢量反演中加入异常周期项参数。     

Vertical deformations from homogeneously processed GRACE and global GPS long-term series

Temporal variations in the geographic distribution of surface mass cause surface displacements. Surface displacements derived from GRACE gravity field coefficient time series also should be observed in GPS coordinate time series, if both time series are sufficiently free of systematic errors. A successful validation can be an important contribution to climate change research, as the biggest contributors to mass variability in the system Earth include the movement of oceanic, atmospheric, and continental water and ice. In our analysis, we find that if the signals are larger than their precision, both geodetic sensor systems see common signals for almost all the 115 stations surveyed. Almost 80% of the stations have their signal WRMS decreased, when we subtract monthly GRACE surface displacements from those observed by GPS data. Almost all other stations are on ocean islands or small peninsulas, where the physically expected loading signals are very small. For a fair comparison, the data (79 months from September 2002 to April 2009) had to be treated appropriately: the GPS data were completely reprocessed with state-of-the-art models. We used an objective cluster analysis to identify and eliminate stations, where local effects or technical artifacts dominated the signals. In addition, it was necessary for both sets of results to be expressed in equivalent reference frames, meaning that net translations between the GPS and GRACE data sets had to be treated adequately. These data sets are then compared and statistically analyzed: we determine the stability (precision) of GRACE-derived, monthly vertical deformation data to be ~1.2 mm, using the data from three GRACE processing centers. We statistically analyze the mean annual signals, computed from the GPS and GRACE series. There is a detailed discussion of the results for five overall representative stations, in order to help the reader to link the displayed criteria of similarity to real data. A series of tests were performed with the goal of explaining the remaining GPS–GRACE residuals.     

利用不同气候类型的GPS可降水量比较研究

利用高精度定位定轨软件GAMIT对2000～2004年的中国地壳运动监测网络数据进行了处理,获得各测站的可降水量时间序列,并以不同气候类型将各测站的可降水量进行了分类比较。在不同气候类型的可降水量的比较中,青藏高原高寒地区的可降水量最低,温带大陆性气候的可降水量次之,温带季风气候的可降水量居中,热带季风气候的可降水量最高,亚热带季风气候的可降水量为次高。在可降水量的图形变化比较中,不同气候类型的可降水量序列的波形变化也不相同,对可降水量峰值时间进行了量化比较。     

Estimates of ocean tide loading displacements and its impact on position time series in Hong Kong using a dense continuous GPS network

Three-dimensional ocean tide loading (OTL) displacements of eight diurnal and semidiurnal constituents at 12 sites in Hong Kong were estimated using 3–7 years of continuous global positioning system (GPS) observations. OTL displacements were estimated using the precise point positioning (PPP) technique on a daily basis and then combined. The OTL displacements obtained by GPS were compared with predictions using seven recent global ocean tide models. The effect of OTL displacements on GPS position time series was also investigated. The study shows that the GPS-derived OTL displacements (excluding K1 and K2 constituents) agree best with those predicted by the GOT4.7 and NAO99b models. The GPS/model agreement is generally at the sub-millimeter level, except for S2, K1, and K2 constituents with relatively large errors. After systematic biases between the GPS and model values are removed, the misfits of all sites for M2, S2, N2, O1, P1, and Q1 are less than 0.5 and 1.0 mm in the horizontal and vertical components, respectively, while larger misfits (within 2.5 mm) are observed for K1 and K2. Integer ambiguity fixing slightly improves the east component of OTL displacement estimates. The study also finds that GPS-derived OTL corrections, instead of model predicts, can be used in daily data processing with the exception of K1 and K2. Including K2 corrections, a secular vertical rate of up to 1 mm/year in position time series can be induced, which needs to be confirmed by further studies.     

利用GPS和GRACE研究澳大利亚地壳垂向季节性变化

为探究重力场恢复与气候实验(gravity recovery and climate experiment,GRACE)卫星与全球定位系统(global positioning system,GPS)两种独立技术获取的因陆地水储量变化引起的地壳垂向季节性位移的一致性,选取澳大利亚27个GPS站点5～10 a的高程时间序...     

Spatiotemporal filtering of regional GNSS network’s position time series with missing data using principle component analysis

The existing spatiotemporal analysis methods suppose that the involved time series are complete and have the same data interval. However missing data inevitably occur in the position time series of Global Navigation Satellite Systems networks for many reasons. In this paper, we develop a modified principal component analysis to extract the Common Mode Error (CME) from the incomplete position time series. The principle of the proposed method is that a time series can be reproduced from its principle components. The method is equivalent to the method of Dong et al. (J Geophys Res 111:3405–3421, 2006) in case of no missing data in the time series and to the extended ‘stacking’ approach under the assumption of a uniformly spatial response. The new method is first applied to extract the CME from the position time series of the Crustal Movement Observation Network of China (CMONOC) over the period of 1999–2009 where the missing data occur in all stations with the different gaps. The results show that the CMEs are significant in CMONOC. The size of the first principle components for the North, East and Up coordinates are as large as 40, 41 and 37 % of total principle components and their spatial responses are not uniform. The minimum amplitudes of the first eigenvectors are only 41, 15 and 29 % for the North, East and Up coordinate components, respectively. The extracted CMEs of our method are close to the data filling method, and the Root Mean Squared error (RMS) values computed from the differences of maximum CMEs between two methods are only 0.31, 0.52 and 1.55 mm for North, East and Up coordinates, respectively. The RMS of the position time series is greatly reduced after filtering out the CMEs. The accuracies of the reconstructed missing data using the two methods are also comparable. To further comprehensively test the efficiency of our method, the repeated experiments are then carried out by randomly deleting different percentages of data at some stations. The results show that the CMEs can be extracted with high accuracy at the non missing-data epochs. And at the missing-data epochs, the accuracy of extracted CMEs has a strong dependence on the number of stations with missing data.     

质量负荷引起地表形变的格林函数和球谐函数方法对比研究

固体地球对地表质量负荷响应,一般采用格林函数或球谐函数方法来计算。两者在数学上是等价的,在实际计算中存在差异。应用地表流体质量负荷变化数据,定量分析了二者在计算地表位移时的精度。结果表明,同一负荷作用下对单点地表位移的计算,两者的精度在误差水平上是一致的,其计算效率也是一致的;但在计算全球1°规则网格站点时,球谐函数方法要比格林函数方法在计算效率上快近100倍。对地表流体质量变化而言,2°空间网格的分辨率一般可满足改正GPS位移数据的精度,且地表流体引起的站点垂直位移变化可解释超过50%的GPS观测方差。     

汶川地震前后四川盆地CORS站运动特性分析

为探讨四川盆地在汶川地震前后地壳运动的规律,本文对2006年3月至2012年9月四川12个连续运行参考站（Continuously Operating Reference Stations, CORS）6年多的观测资料进行了详细的处理和分析。通过比较震前震后CORS站运动的相似性和差异性发现：由于震后形变的影响,震后第一主成分(Principle Component First, PC1)的贡献率下降20%~40%;震后多站的速度场发生了明显的变化,突出表现在PIXI,CHDU,MYAN等站,YAAN,QLAI站处于闭锁状态;盆地垂直分量的周年振幅在震前一年偏大,震后逐年减小。以上发现表明汶川地震可能改变了四川盆地原有的运动特性。     

基于 EEMD 和 Husrt 指数的 GNSS 基准站的垂向速率估计

连续的全球卫星导航系统(GNSS)基准站的坐标时间序列中包含了复杂的噪声信号、非构造形变及其它因素的影响,尤其在垂直方向,对GNSS基准站在国际地球参考框架(ITRF)下的运动速率估计产生了较大的干扰。为进一步提高速率精度,文中采用整体模态分解(EEMD)方法对GNSS基准站的垂向观测时间序列进行分解,并根据各种信号的Hurst值进行分类及重构为噪声信号、季节性信号和长期趋势信号,采用最小二乘方法拟合长期趋势信号得到垂向速率。通过对中国大陆构造环境监测网络(CMONOC)的GNSS台站从2001—2013年近13a的垂向坐标时间序列的实例分析,采用基于EEMD和Husrt指数的最小二乘法能够准确地估计GNSS基准站的垂向速率。