The Bayesian detection of discontinuities in a polynomial regression and its application to the cycle-slip problem

This paper deals with the problem of detecting and correcting cycle-slips in Global Navigation Satellite System (GNSS) phase data by exploiting the Bayesian theory. The method is here applied to undifferenced observations, because repairing cycle-slips already at this stage could be a useful pre-processing tool, especially for a network of permanent GNSS stations. If a dual frequency receiver is available, the cycle-slips can be easily detected by combining two phase observations or phase and range observations from a single satellite to a single receiver. These combinations, expressed in a distance unit form, are completely free from the geometry and depend only on the ionospheric effect, on the electronic biases and on the initial integer ambiguities; since these terms are expected to be smooth in time, at least in a short period, a cycle-slip in one or both the two carriers can be modelled as a discontinuity in a polynomial regression. The proposed method consists in applying the Bayesian theory to compute the marginal posterior distribution of the discontinuity epoch and to detect it as a maximum a posteriori (MAP) in a very accurate way. Concerning the cycle-slip correction, a couple of simultaneous integer slips in the two carriers is chosen by maximazing the conditional posterior distribution of the discontinuity amplitude given the detected epoch. Numerical experiments on simulated and real data show that the discontinuities with an amplitude 2 or 3 times larger than the noise standard deviation are successfully identified. This means that the Bayesian approach is able to detect and correct cycle-slips using undifferenced GNSS observations even if the slip occurs by one cycle. A comparison with the scientific software BERNESE 5.0 confirms the good performance of the proposed method, especially when data sampled at high frequency (e.g. every 1 s or every 5 s) are available.     

GNSS three carrier ambiguity resolution using ionosphere-reduced virtual signals

This paper presents a general modeling strategy for ambiguity resolution (AR) and position estimation (PE) using three or more phase-based ranging signals from a global navigation satellite system (GNSS). The proposed strategy will identify three best “virtual” signals to allow for more reliable AR under certain observational conditions characterized by ionospheric and tropospheric delay variability, level of phase noise and orbit accuracy. The selected virtual signals suffer from minimal or relatively low ionospheric effects, and thus are known as ionosphere-reduced virtual signals. As a result, the ionospheric parameters in the geometry-based observational models can be eliminated for long baselines, typically those of length tens to hundreds of kilometres. The proposed modeling comprises three major steps. Step 1 is the geometry-free determination of the extra-widelane (EWL) formed between the two closest L-band carrier measurements, directly from the two corresponding code measurements. Step 2 forms the second EWL signal and resolves the integer ambiguity with a geometry-based estimator alone or together with the first EWL. This is followed by a procedure to correct for the first-order ionospheric delay using the two ambiguity-fixed widelane (WL) signals derived from the integer-fixed EWL signals. Step 3 finds an independent narrow-lane (NL) signal, which is used together with a refined WL to resolve NL ambiguity with geometry-based integer estimation and search algorithms. As a result, the above two AR processes performed with WL/NL and EWL/WL signals respectively, either in sequence or in parallel, can support real time kinematic (RTK) positioning over baselines of tens to hundreds of kilometres, thus enabling centimetre-to-decimentre positioning at the local, regional and even global scales in the future.     

Integer aperture bootstrapping: a new GNSS ambiguity estimator with controllable fail-rate

In this contribution, we introduce a new bootstrap-based method for Global Navigation Satellite System (GNSS) carrier-phase ambiguity resolution. Integer bootstrapping is known to be one of the simplest methods for integer ambiguity estimation with close-to-optimal performance. Its outcome is easy to compute due to the absence of an integer search, and its performance is close to optimal if the decorrelating Z-transformation of the LAMBDA method is used. Moreover, the bootstrapped estimator is presently the only integer estimator for which an exact and easy-to-compute expression of its fail-rate can be given. A possible disadvantage is, however, that the user has only a limited control over the fail-rate. Once the underlying mathematical model is given, the user has no freedom left in changing the value of the fail-rate. Here, we present an ambiguity estimator for which the user is given additional freedom. For this purpose, use is made of the class of integer aperture estimators as introduced in Teunissen (2003). This class is larger than the class of integer estimators. Integer aperture estimators are of a hybrid nature and can have integer outcomes as well as non-integer outcomes. The new estimator is referred to as integer aperture bootstrapping. This new estimator has all the advantages known from integer bootstrapping with the additional advantage that its fail-rate can be controlled by the user. This is made possible by giving the user the freedom over the aperture of the pull-in region. We also give an exact and easy-to-compute expression for its controllable fail-rate.     

Preliminary results on kinematic model of tectonic blocks derived from high precision GPS observations in Southwest China

ＰｒｅｌｉｍｉｎａｒｙｒｅｓｕｌｔｓｏｎｋｉｎｅｍａｔｉｃｍｏｄｅｌｏｆｔｅｃｔｏｎｉｃｂｌｏｃｋｓｄｅｒｉｖｅｄｆｒｏｍｈｉｇｈｐｒｅｃｉｓｉｏｎＧＰＳｏｂｓｅｒｖａｔｉｏｎｓｉｎＳｏｕｔｈｗｅｓｔＣｈｉｎａＬＩＲＥＮＨＵＡＮＧ１）（黄立人）...     

丰城洛市地区推,滑覆构造特征及其成因探讨

介绍了洛市地区推、滑覆构造特征、运动学标志及其形成时期、成因模式。     

Recent Advances in Prediction and Processing of Strong Ground Motions

We present an overview of our recent results on utilizing small earthquakes in the earthquake engineering practice. Site-specific ground motion time-histories of large earthquakes can be successfully simulated using recordings of small earthquakes which are often referred to as 'empirical Green's functions' in seismology. Another important practical problem is whether and how these observations can be used in seismic risk studies which are based on empirical attenuation relations for ground motion parameters. We study a possibility of extrapolating attenuation relations for small earthquakes, to larger magnitudes using the data from the Garner Valley downhole array in Southern California. Finally we introduce efficient ground motion processing techniques in frequency- and time-domains and apply them to site response estimation.     

改进的能量守恒方法及其在CHAMP重力场恢复中的应用

An efficient method for gravity field determination from CHAMP orbits and accelerometer data is referred to as the energy balance approach. A new CHAMP gravity field recovery strategy based on the improved energy balance approach IS developed in this paper. The method simultaneously solves the spherical harmonic coefficients, daily Integration constant, scale and bias parameters. Two 60 degree and order gravitational potential models, XISM-CHAMPO1S from the classical energy balance approach, and XISM-CHAMPO2S from the improved energy balance, are determined using about one year＇s worth of CHAMP kinematic orbits from TUM and accelerometer data from GFZ. Comparisons among XISM-CHAMPO1S, XISM-CHAMPO2S, EIGEN-CGO3C, EIGEN-CHAMPO3S, EIGEN2, ENIGNIS and EGM96 are made. The results show that the XISM-CHAMPO2S model is more accurate than EGM96, EIGENIS, EIGEN2 and XISM-CHAMPO1S at the same degree and order, and has almost the same accuracy as EIGEN-CHAMPO3S.     

基于FFT与小波变换的SOI与GNSS ZTD的周期变化影响研究

为探究ENSO事件对GNSSZTD(水汽)周期变化的影响及其相互关系,以河北省为例开展ENSO事件对GNSSZTD及其周期变化的影响研究.首先利用快速傅里叶变换方法筛选出南方涛动指数(SOI)与GNSS水汽的共同周期,再利用小波变换提取GNSS水汽与SOI共同周期所在的高频项,并将重构的高频项与SOI进行相关性分析.结...     

中国沿海HY-2A校正微波辐射计水汽含量数据精度检验方法

为评价“海洋二号”卫星（HaiYang-2A, HY-2A）校正微波辐射计（Calibration Microwave Radiometer, CMR）近海水汽产品精度,以中国沿海全球导航卫星系统（Global Navigation Satellite System, GNSS）业务观测站数据和欧洲中期天气预报中心（European Centre for Medium-Range Weather Forecasts, ECMWF)发布的第五代大气再分析资料（ECMWF Reanalysis 5, ERA5）作为验证数据。首先对选取的GNSS业务观测站数据和CMR水汽含量数据进行时空匹配,两者的观测时间一致、空间范围取为100 km;然后利用精密单点定位方法反演GNSS业务观测站上空的大气可降水量（Precipitable Water Vapor, PWV）,同时对1 h分辨率的ERA5再分析资料内插计算,得到CMR水汽数据点处的ERA5 PWV;最后以GNSS PWV和ERA5 PWV为参考,分析2015年CMR水汽产品精度和偏差时空分布。结果表明,CMR水汽含量和GNSS PWV、E...