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

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

Monument-antenna effects on GPS coordinate time series with application to vertical rates in Antarctica

We examine the electromagnetic coupling of a GPS antenna–monument pair in terms of its simulated affect on long GPS coordinate time series. We focus on the Earth and Polar Observing System (POLENET) monument design widely deployed in Antarctica and Greenland in projects interested particularly in vertical velocities. We base our tests on an absolute robot calibration that included the top ~0.15 m of the monument and use simulations to assess its effect on site coordinate time series at eight representative POLENET sites in Antarctica over the period 2000.0–2011.0. We show that the neglect of this calibration would introduce mean coordinate bias, and most importantly for velocity estimation, coordinate noise which is highly sensitive to observation geometry and hence site location and observation period. Considering only sub-periods longer than 2.5 years, we show vertical site velocities may be biased by up to ±0.4 mm/year, and biases up to 0.2 mm/year may persist for observation spans of 8 years. Changing between uniform and elevation-dependent observation weighting alters the time series but does not remove the velocity biases, nor does ambiguity fixing. The effect on the horizontal coordinates is negligible. The ambiguities fixed series spectra show noise between flicker and random walk with near-white noise at the highest frequencies, with mean spectral indices (frequencies <20 cycles per year) of approximately −1.3 (uniform weighting) and −1.4 (elevation-dependent weighting). While the results are likely highly monument specific, they highlight the importance of accounting for monument effects when analysing vertical coordinate time series and velocities for the highest precision and accuracy geophysical studies.     

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基准站时间序列的影响分析

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

华北地区季节性垂直位移信号提取及物理机制分析

华北地区地表垂直位移呈现显著的季节性波动。为探究其驱动源信号,利用独立分量分析方法对全球导航定位系统20个基准站的垂直位移时间序列进行时空分解,获取前6个独立分量,并与环境负载造成的垂直位移进行比较分析。结果表明,第一个独立分量呈现显著周年变化和相对均匀的空间响应,与大气负载位移时间序列的平均相关系数为0.60,两者吻合良好。第二和第三个独立分量均是周年和年内信号的组合,其空间响应反映出局部集聚特征,与华北地区水资源的分布状况相类似,合并这两个独立分量后与水文负载位移时间序列的平均相关系数为0.50。由此表明大气和水文负载是驱动华北地区季节性垂直位移的主要源信号。     

Tight integration of GPS observations and persistent scatterer InSAR for detecting vertical ground motion in Hong Kong

A method for tightly integrating GPS observations and the persistent scatterer (PS) interferometric synthetic aperture radar (InSAR) is proposed to detect vertical ground motion in Hong Kong, China. The tropospheric zenith wet delays (ZWD) derived from GPS observations at sites in the SAR scenes are used first to model and correct for the tropospheric effects in the interferograms generated from the SAR images. The vertical motion rates (VMR) of the GPS sites determined based on the continuous GPS observations obtained at the sites are then used as constraints in the PS InSAR solutions to further enhance the quality of the solutions. Data from 12 continuous tracking GPS sites in Hong Kong and 8 ENVISAT ASAR images acquired during 2006–2007 are used to demonstrate the application of the proposed method. The results from the example show that the method is capable of detecting accurately the vertical ground motion.     

垂直形变监测中引入重力观测值的可行性研究

全面分析重力测量的各种误差和目前能达到的精度后，对垂直形变监测中引入重力观测值的可行性进行了具体研究。结果表明，若仅是为了测定点的垂直位移的高程速率，增加重力测量数据后对形变分析的精度增益十分有限。但是对于大跨度垂直形变监测和局部地区突发性垂直形变测定，重力测量方法可以及时掌握垂直形变的趋势     

热膨胀效应对GNSS短基线坐标时间序列影响研究

天线观测墩及基岩的热膨胀效应会造成全球卫星导航系统(GNSS)基准站坐标时间序列非线性变化,本文利用长期观测数据评估了热膨胀效应对GNSS短基线时间序列周期特性的影响. 结果表明:热膨胀效应在高程上可以引起约1 mm的季节性周年信号,可以很好地解释GNSS短基线时间序列中的周期信号,而对平面方向的贡献不大,但是观测墩因日光照射受热不均产生的热弯曲对平面方向上的影响应予以考虑.      

Deriving time-series three-dimensional displacements of mining areas from a single-geometry InSAR dataset

This paper presents a method for deriving time-series three-dimensional (3-D) displacements of mining areas from a single-geometry interferometric synthetic aperture radar (InSAR) dataset (hereafter referred to as the SGI-based method). This is mainly aimed at overcoming the limitation of the traditional multi-temporal InSAR techniques that require SAR data from at least three significantly different imaging geometries to fully retrieve time-series 3-D displacements of mining areas. The SGI-based method first generates the multi-temporal observations of the mining-induced vertical subsidence from the single-geometry InSAR data, using a previously developed method for retrieving 3-D mining-related displacements from a single InSAR pair (SIP). The weighted least-squares solutions of the time series of vertical subsidence are estimated from these generated multi-temporal observations of vertical subsidence. Finally, the time series of horizontal motions in the east and north directions are estimated using the proportional relationship between the horizontal motion and the subsidence gradient of the mining area, on the basis of the SGI-derived time series of vertical subsidence. Seven ascending ALOS PALSAR images from the Datong mining area of China were used to test the proposed SGI-based method. The results suggest that the SGI-based method is effective. The SGI-based method not only extends the SIP-based method to time-series 3-D displacement retrieval from a single-geometry InSAR dataset, but also limits the uncertainty propagation from InSAR-derived deformation to the estimated 3-D displacements.     

Vertical displacement rates in the Upper Rhine Graben area derived from precise leveling

The recent vertical displacement field of the Upper Rhine Graben (URG) located in the tri-national region between Germany, France and Switzerland is investigated using repeatedly measured leveling data. We estimate vertical displacement rates at leveling benchmarks by applying a kinematic network adjustment on more than 40,000 height differences measured by German, French and Swiss surveying agencies. Focusing on an optimal solution for the adjusted rates in the URG area also historical data (measured before 1900) are used, significantly increasing the time span of available measurements and the number of transnational connections between the three countries. To account for inhomogeneities apparent in the database, we apply an iterative variance component estimation within the adjustment procedure, particularly revealing more realistic information on the accuracy of the estimated rates. A special focus within our analysis is put on the statistical testing of gross errors in the observations and model-related errors at benchmarks with non-linear movement. As some of the estimated vertical rates behave significantly different compared to the vertical rates of adjacent benchmarks, a filtering of outliers is applied after the adjustment procedure. The resulting map of linear height changes in an area of 280 km in N–S and 230 km in E–W direction provides detailed insight into the recent vertical displacements of the URG and neighboring regions. In the German part of the study area, it was possible for the first time to consistently constrain an average subsidence rate of 0.5 mm/a ( \(\pm \) 0.2 mm/a) of the Graben interior w.r.t. the Black Forest. In addition to the tectonic displacements, some man-induced surface movements, e.g., caused by oil and groundwater extraction, are observed and discussed.     

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与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%。     

太阳角的计算方法及其在遥感中的应用

大多数遥感探测器所接收的电磁辐射能量主要源自地球表面对太阳光的反时,而反射能量的大小和当时太阳与卫星星下点的相对位置有密切的关系,因此,反映太阳位置的量──太阳赤经、太阳赤纬、太阳时角、太阳高度角和太阳天顶角等在遥感工作中极其有用。本文全面论述了这些有关太阳角度的计算方法,并给出其在遥感工作中的几个应用实例。     

Real time estimation of position and the gravity vector with an inertial survey system

Based on sine functions for fitting horizontal velocity errors in an inertial survey system and a step or slope function approximation of the vertical deflection between stops, a positioning accuracy of 10 m (CEP) can be achieved with a car-borne system – a modified -21 aircraft navigation system, with a 10 min travel period between Zero Velocity Update Times (ZUPTs) and a total 1.5 hour survey time. A rough estimation of the vertical deflection can also be expected in real time out of the approach. Also a forward-backward-forward run is presented which demonstrates determination of the gravity anomaly. Received 10 March 1995; Accepted 8 July 1996     

The spherical horizontal and spherical vertical boundary value problem – vertical deflections and geoidal undulations – the completed Meissl diagram

 In a comparison of the solution of the spherical horizontal and vertical boundary value problems of physical geodesy it is aimed to construct downward continuation operators for vertical deflections (surface gradient of the incremental gravitational potential) and for gravity disturbances (vertical derivative of the incremental gravitational potential) from points on the Earth's topographic surface or of the three-dimensional (3-D) Euclidean space nearby down to the international reference sphere (IRS). First the horizontal and vertical components of the gravity vector, namely spherical vertical deflections and spherical gravity disturbances, are set up. Second, the horizontal and vertical boundary value problem in spherical gravity and geometry space is considered. The incremental gravity vector is represented in terms of vector spherical harmonics. The solution of horizontal spherical boundary problem in terms of the horizontal vector-valued Green function converts vertical deflections given on the IRS to the incremental gravitational potential external in the 3-D Euclidean space. The horizontal Green functions specialized to evaluation and source points on the IRS coincide with the Stokes kernel for vertical deflections. Third, the vertical spherical boundary value problem is solved in terms of the vertical scalar-valued Green function. Fourth, the operators for upward continuation of vertical deflections given on the IRS to vertical deflections in its external 3-D Euclidean space are constructed. Fifth, the operators for upward continuation of incremental gravity given on the IRS to incremental gravity to the external 3-D Euclidean space are generated. Finally, Meissl-type diagrams for upward continuation and regularized downward continuation of horizontal and vertical gravity data, namely vertical deflection and incremental gravity, are produced. Received: 10 May 2000 / Accepted: 26 February 2001     

降雨对Ku波段旋转扇形波束散射计影响

星载散射计测量的归一化后向散射截面是有关海面毛细重力波的大小和方向的函数:散射计的回波强度与海面毛细重力波的振幅成正比;风向对后向散射系数具有调制作用。因此,可以利用散射计的数据,根据地球物理模型反演出具有高精确度、无雨和中低风速条件下的海面风场矢量。然而高达10%的散射计测量数据会受到降雨影响(Nie,etal.,2008),尤其工作在Ku波段的散射计。降雨对海面散射幅度的影响主要包括:(1)降雨对雷达波的衰减和散射;(2)降雨改变海洋水面形态和海面粗糙度。     

Time transfer using GPS carrier phase: error propagation and results

 A joint time-transfer project between the Astronomical Institute of the University of Berne (AIUB) and the Swiss Federal Office of Metrology and Accreditation (METAS) was initiated to investigate the power of the time transfer using GPS carrier phase observations. Studies carried out in the context of this project are presented. The error propagation for the time-transfer solution using GPS carrier phase observations was investigated. To this purpose a simulation study was performed. Special interest was focussed on errors in the vertical component of the station position, antenna phase-center variations and orbit errors. A constant error in the vertical component introduces a drift in the time-transfer results for long baselines in east–west directions. The simulation study was completed by investigating the profit for time transfer when introducing the integer carrier phase ambiguities from a double-difference solution. This may reduce the drift in the time-transfer results caused by constant vertical error sources. The results from the present time-transfer solution are shown in comparison to results obtained with independent time-transfer techniques. The interpretation of the comparison benefits from the investigations of the error propagation study. Two types of solutions are produced on a regular basis at AIUB: one based on the rapid orbits from CODE, the other on the CODE final orbits. The rapid solution is available the day after the observations and has nearly the same quality as the final solution, which has a latency of about one week. The differences between these two solutions are below the nanosecond level. The differences from independent time-transfer techniques such as TWSTFT (two-way satellite time and frequency transfer) are a few nanoseconds for both products. Received: 15 November 2001 / Accepted: 6 September 2002 Correspondence to:R. Dach     

联合小波和方差分量估计方法分析中国IGS测站时间序列变化特征

利用小波谱分析了中国区域IGS（International GNSS Service）站坐标时间序列中周期信号成分及其随时间的变化,根据小波谱中小波系数的能量在时域上的变化对时间序列进行分段,然后采用方差-协方差验后估计法分析了白噪声+闪烁噪声组合模型下测站在不同时间段的运动特征及变化规律。结果表明,中国区域的IGS站运动特征在不同时间段内是变化的,其中URUM站垂直速度变化大于4 mm,BJFS以及CHAN站水平速度变化大于2 mm,其余站速度变化均在2 mm以内;周年、半周年信号的振幅也均有差异,其中BJFS站垂直分量周年振幅变化最大,超过5 mm;此外,强度较弱的其他周期项存在于分段时间序列中,但其对速度、周期振幅以及噪声分量估计结果影响较小;对时间序列进行分段估计,所得的参数不同,其相差程度可反映出不同时间段内的板块运动等地球物理因素变化的差异。     

CORRECTION OF PLANE-TABLE HEIGHTS FOR CURVATURE OF THE EARTH AND ATMOSPHERIC REFRACTION

Abstract

Most text-books on surveying limit their discussion of the correction of vertical angles for curvature of the earth and atmospheric refraction to the correction of angles taken with a theodolite during triangulation and omit any reference to those taken with a clinometer. This is rather illogical, as in well-observed triangulation, with all vertical angles measured in both directions, no correction for these effects is necessary, whilst in plane-tabling on small scales where sketching at considerable distances is frequently employed the application of corrections for these effects is essential.     

联合EMD-HD和小波分解的GNSS坐标时间序列降噪分析EI北大核心CSCD

针对经验模态分解(empirical mode decomposition,EMD)在GNSS坐标时间序列的降噪过程中存在筛选准则的选取和模态混叠效应等问题,本文引入Hausdorff距离(Hausdorff distance,HD)筛选准则并结合小波分解(wavelet decomposition,WD),提出EMD-HD&WD算法。通过对我国大陆构造环境监测网络149个GNSS测站的垂向坐标时间序列降噪处理,分别利用复合指标T值、测站的速度不确定度和闪烁噪声振幅验证算法的可靠性和普适性。结果显示:HD优于现有的筛选准则;EMD-HD&WD算法对测站的速度不确定度和闪烁噪声振幅的平均改正率均为88.4%。分析表明,本文算法能够有效识别和剔除噪声并且改善EMD的模态混叠效应,提高GNSS垂向坐标时间序列的模型精度。