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

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

Plate Motion of India and Interseismic Strain in the Nepal Himalaya from GPS and DORIS Measurements

We analyse geodetically estimated deformation across the Nepal Himalaya in order to determine the geodetic rate of shortening between Southern Tibet and India, previously proposed to range from 12 to 21 mm yr^?1. The dataset includes spirit-levelling data along a road going from the Indian to the Tibetan border across Central Nepal, data from the DORIS station on Everest, which has been analysed since 1993, GPS campaign measurements from surveys carried on between 1995 and 2001, as well as data from continuous GPS stations along a transect at the logitude of Kathmandu operated continuously since 1997. The GPS data were processed in International Terrestrial Reference Frame 2000 (ITRF2000), together with the data from 20 International GNSS Service (IGS) stations and then combined using quasi- observation combination analysis (QOCA). Finally, spatially complementary velocities at stations in Southern Tibet, initially determined in ITRF97, were expressed in ITRF2000. After analysing previous studies by different authors, we determined the pole of rotation of the Indian tectonic plate to be located in ITRF2000 at 51.409±1.560° N and ?10.915±5.556°E, with an angular velocity of 0.483±0.015°. Myr^?1. Internal deformation of India is found to be small, corresponding to less than about 2 mm yr^?1 of baseline change between Southern India and the Himalayan piedmont. Based on an elastic dislocation model of interseismic strain and taking into account the uncertainty on India plate motion, the mean convergence rate across Central and Eastern Nepal is estimated to 19±2.5 mm yr^?1, (at the 67% confidence level). The main himalayan thrust (MHT) fault was found to be locked from the surface to a depth of about 20 km over a width of about 115 km. In these regions, the model parameters are well constrained, thanks to the long and continuous time-series from the permanent GPS as well as DORIS data. Further west, a convergence rate of 13.4±5 mm yr^?1, as well as a fault zone, locked over 150 km, are proposed. The slight discrepancy between the geologically estimated deformation rate of 21±1.5 mm yr^?1 and the 19±2.5 mm yr^?1 geodetic rate in Central and Eastern Nepal, as well as the lower geodetic rate in Western Nepal compared to Eastern Nepal, places bounds on possible temporal variations of the pattern and rate of strain in the period between large earthquakes in this region.     

Estimation and evaluation of real-time precipitable water vapor from GLONASS and GPS

The revitalized Russian GLONASS system provides new potential for real-time retrieval of zenith tropospheric delays (ZTD) and precipitable water vapor (PWV) in order to support time-critical meteorological applications such as nowcasting or severe weather event monitoring. In this study, we develop a method of real-time ZTD/PWV retrieval based on GLONASS and/or GPS observations. The performance of ZTD and PWV derived from GLONASS data using real-time precise point positioning (PPP) technique is carefully investigated and evaluated. The potential of combining GLONASS and GPS data for ZTD/PWV retrieving is assessed as well. The GLONASS and GPS observations of about half a year for 80 globally distributed stations from the IGS (International GNSS Service) network are processed. The results show that the real-time GLONASS ZTD series agree quite well with the GPS ZTD series in general: the RMS of ZTD differences is about 8 mm (about 1.2 mm in PWV). Furthermore, for an inter-technique validation, the real-time ZTD estimated from GLONASS-only, GPS-only, and the GPS/GLONASS combined solutions are compared with those derived from very long baseline interferometry (VLBI) at colocated GNSS/VLBI stations. The comparison shows that GLONASS can contribute to real-time meteorological applications, with almost the same accuracy as GPS. More accurate and reliable water vapor values, about 1.5–2.3 mm in PWV, can be achieved when GLONASS observations are combined with the GPS ones in the real-time PPP data processing. The comparison with radiosonde data further confirms the performance of GLONASS-derived real-time PWV and the benefit of adding GLONASS to stand-alone GPS processing.     

Homogeneous reprocessing of GPS,GLONASS and SLR observations

The International GNSS Service (IGS) provides operational products for the GPS and GLONASS constellation. Homogeneously processed time series of parameters from the IGS are only available for GPS. Reprocessed GLONASS series are provided only by individual Analysis Centers (i. e. CODE and ESA), making it difficult to fully include the GLONASS system into a rigorous GNSS analysis. In view of the increasing number of active GLONASS satellites and a steadily growing number of GPS+GLONASS-tracking stations available over the past few years, Technische Universität Dresden, Technische Universität München, Universität Bern and Eidgenössische Technische Hochschule Zürich performed a combined reprocessing of GPS and GLONASS observations. Also, SLR observations to GPS and GLONASS are included in this reprocessing effort. Here, we show only SLR results from a GNSS orbit validation. In total, 18 years of data (1994–2011) have been processed from altogether 340 GNSS and 70 SLR stations. The use of GLONASS observations in addition to GPS has no impact on the estimated linear terrestrial reference frame parameters. However, daily station positions show an RMS reduction of 0.3 mm on average for the height component when additional GLONASS observations can be used for the time series determination. Analyzing satellite orbit overlaps, the rigorous combination of GPS and GLONASS neither improves nor degrades the GPS orbit precision. For GLONASS, however, the quality of the microwave-derived GLONASS orbits improves due to the combination. These findings are confirmed using independent SLR observations for a GNSS orbit validation. In comparison to previous studies, mean SLR biases for satellites GPS-35 and GPS-36 could be reduced in magnitude from \(-35\) and \(-38\)  mm to \(-12\) and \(-13\)  mm, respectively. Our results show that remaining SLR biases depend on the satellite type and the use of coated or uncoated retro-reflectors. For Earth rotation parameters, the increasing number of GLONASS satellites and tracking stations over the past few years leads to differences between GPS-only and GPS+GLONASS combined solutions which are most pronounced in the pole rate estimates with maximum 0.2 mas/day in magnitude. At the same time, the difference between GLONASS-only and combined solutions decreases. Derived GNSS orbits are used to estimate combined GPS+GLONASS satellite clocks, with first results presented in this paper. Phase observation residuals from a precise point positioning are at the level of 2 mm and particularly reveal poorly modeled yaw maneuver periods.     

Joint estimation of GPS/BDS real-time clocks and initial results

We present the joint estimation model for Global Positioning System/BeiDou Navigation Satellite System (GPS/BDS) real-time clocks and present the initial satellite clock solutions determined from 106 stations of the international GNSS service multi-GNSS experiment and the BeiDou experimental tracking stations networks for 1 month in December, 2012. The model is shown to be efficient enough to have no practical computational limit for producing 1-Hz clock updates for real-time applications. The estimated clocks were assessed through the comparison with final clock products and the analysis of post-fit residuals. Using the estimated clocks and corresponding orbit products (GPS ultra-rapid-predicted and BDS final orbits), the root-mean-square (RMS) values of coordinate differences from ground truth values are around 1 and 2–3 cm for GPS-only and BDS-only daily mean static precise point positioning (PPP) solutions, respectively. Accuracy of GPS/BDS combined static PPP solutions falls in between that of GPS-only and BDS-only PPP results, with RMS values approximately 1–2 cm in all three components. For static sites, processed in the kinematic PPP mode, the daily RMS values are normally within 4 and 6 cm after convergence for GPS-only and BDS-only results, respectively. In contrast, the combined GPS/BDS kinematic PPP solutions show higher accuracy and shorter convergence time. Additionally, the BDS-only kinematic PPP solutions using clock products derived from the proposed joint estimation model were superior compared to those computed using the single-system estimation model.     

基于GPS的四川地区现时地壳形变分析

对四川区域网多年数据进行处理,得到欧亚框架下的测站运动速度,采用刚性地块假设下的最小二乘拟合方法,得到川滇、川青、扬子地块的运动速度分别为(17.6±1.0)mm/a、(11.4±0.7)mm/a、(9.6±1.1)mm/a,地块运动方向逐渐由SE向SEE变化;作为川滇地块东边界的鲜水河断裂、安宁河-则木河断裂运动速度远大于龙门山断裂,显示了其较强的活动性。     

Estimation of aboveground biomass and carbon in a tropical rain forest in Gabon using remote sensing and GPS data

The knowledge of biomass stocks in tropical forests is critical for climate change and ecosystem services studies. This research was conducted in a tropical rain forest located near the city of Libreville (the capital of Gabon), in the Akanda Peninsula. The forest cover was stratified in terms of mature, secondary and mangrove forests using Landsat-ETM data. A field inventory was conducted to measure the required basic forest parameters and estimate the aboveground biomass (AGB) and carbon over the different forest classes. The Shuttle Radar Topography Mission (SRTM) data were used in combination with ground-based GPS measurements to derive forest heights. Finally, the relationships between the estimated heights and AGB were established and validated. Highest biomass stocks were found in the mature stands (223 ± 37 MgC/ha), followed by the secondary forests (116 ± 17 MgC/ha) and finally the mangrove forests (36 ± 19 MgC/ha). Strong relationships were found between AGB and forest heights (R² > 0.85).     

顾及共模误差的大区域GPS网坐标时间序列噪声分析

噪声分析对GPS时间序列分析有着重要影响,然而针对时间跨度较长的大尺度GPS网的共模误差相关研究较少。本文选取了平均基线长度大于2000 km的欧洲地区9个GPS台站2006-2014年的数据,使用主成分分析法剔除坐标时间序列的共模误差,同时利用极大似然估计的方法对滤波前后的时间序列进行了噪声分析。结果表明,欧洲地区广域GPS网的噪声模型存在多样性,各个分量具有不同的噪声特性,主要表现为白噪声+闪烁噪声、白噪声+幂率噪声,少部分台站N、E两个方向含有随机漫步噪声。经过空间滤波后,部分台站最优噪声模型发生改变,但仍以白噪声+闪烁噪声、白噪声+幂率噪声为主。滤波对N、E方向速度场影响为0.2 mm/a,U方向速度场影响为0.5 mm/a。     

Assessment of PPP integer ambiguity resolution using GPS,GLONASS and BeiDou (IGSO,MEO) constellations

Although integer ambiguity resolution (IAR) can improve positioning accuracy considerably and shorten the convergence time of precise point positioning (PPP), it requires an initialization time of over 30 min. With the full operation of GLONASS globally and BDS in the Asia–Pacific region, it is necessary to assess the PPP–IAR performance by simultaneous fixing of GPS, GLONASS, and BDS ambiguities. This study proposed a GPS + GLONASS + BDS combined PPP–IAR strategy and processed PPP–IAR kinematically and statically using one week of data collected at 20 static stations. The undifferenced wide- and narrow-lane fractional cycle biases for GPS, GLONASS, and BDS were estimated using a regional network, and undifferenced PPP ambiguity resolution was performed to assess the contribution of multi-GNSSs. Generally, over 99% of a posteriori residuals of wide-lane ambiguities were within ±0.25 cycles for both GPS and BDS, while the value was 91.5% for GLONASS. Over 96% of narrow-lane residuals were within ±0.15 cycles for GPS, GLONASS, and BDS. For kinematic PPP with a 10-min observation time, only 16.2% of all cases could be fixed with GPS alone. However, adding GLONASS improved the percentage considerably to 75.9%, and it reached 90.0% when using GPS + GLONASS + BDS. Not all epochs could be fixed with a correct set of ambiguities; therefore, we defined the ratio of the number of epochs with correctly fixed ambiguities to the number of all fixed epochs as the correct fixing rate (CFR). Because partial ambiguity fixing was used, when more than five ambiguities were fixed correctly, we considered the epoch correctly fixed. For the small ratio criteria of 2.0, the CFR improved considerably from 51.7% for GPS alone, to 98.3% when using GPS + GLONASS + BDS combined solutions.     

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

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

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.     

Impact of seasonal station motions on VLBI UT1 intensives results

UT1 estimates obtained from the very long baseline interferometry (VLBI) Intensives data depend on the station displacement model used during processing. In particular, because of seasonal variations, the instantaneous station position during the specific intensive session differs from the position predicted by the linear model generally used. This can cause systematic errors in UT1 Intensives results. In this paper, we first investigated the seasonal signal in the station displacements for the 5 VLBI antennas participating in UT1 Intensives observing programs, along with the 8 collocated GPS stations. It was found that a significant annual term is present in the time series for most stations, and its amplitude can reach 8 mm in the height component, and 2 mm in horizontal components. However, the annual signals found in the displacements of the collocated VLBI and GPS stations at some sites differ substantially in amplitude and phase. The semiannual harmonics are relatively small and unstable, and for most stations no prevailing signal was found in the corresponding frequency band. Then two UT1 Intensives series were computed with and without including the seasonal term found in the previous step in the station movement model. Comparison of these series has shown that neglecting the seasonal station position variations can cause a systematic error in UT1 estimates, which can exceed 1  $\upmu $ s, depending on the observing program.     

Precise geoid determination over Sweden using the Stokes–Helmert method and improved topographic corrections

 Four different implementations of Stokes' formula are employed for the estimation of geoid heights over Sweden: the Vincent and Marsh (1974) model with the high-degree reference gravity field but no kernel modifications; modified Wong and Gore (1969) and Molodenskii et al. (1962) models, which use a high-degree reference gravity field and modification of Stokes' kernel; and a least-squares (LS) spectral weighting proposed by Sj?berg (1991). Classical topographic correction formulae are improved to consider long-wavelength contributions. The effect of a Bouguer shell is also included in the formulae, which is neglected in classical formulae due to planar approximation. The gravimetric geoid is compared with global positioning system (GPS)-levelling-derived geoid heights at 23 Swedish Permanent GPS Network SWEPOS stations distributed over Sweden. The LS method is in best agreement, with a 10.1-cm mean and ±5.5-cm standard deviation in the differences between gravimetric and GPS geoid heights. The gravimetric geoid was also fitted to the GPS-levelling-derived geoid using a four-parameter transformation model. The results after fitting also show the best consistency for the LS method, with the standard deviation of differences reduced to ±1.1 cm. For comparison, the NKG96 geoid yields a 17-cm mean and ±8-cm standard deviation of agreement with the same SWEPOS stations. After four-parameter fitting to the GPS stations, the standard deviation reduces to ±6.1 cm for the NKG96 geoid. It is concluded that the new corrections in this study improve the accuracy of the geoid. The final geoid heights range from 17.22 to 43.62 m with a mean value of 29.01 m. The standard errors of the computed geoid heights, through a simple error propagation of standard errors of mean anomalies, are also computed. They range from ±7.02 to ±13.05 cm. The global root-mean-square error of the LS model is the other estimation of the accuracy of the final geoid, and is computed to be ±28.6 cm. Received: 15 September 1999 / Accepted: 6 November 2000     

中国大陆GPS速度场的球面小波模型及多尺度特征分析

利用1999—2009年间中国大陆共1068个GPS站点在东方向、北方向的速度值,采用DOG球面小波多尺度分析方法,建立了中国大陆东方向、北方向多尺度速度场。球面小波模型的尺度主要根据观测站点的密度来确定,利用检核点上的已知速度与模型速度之间的均方差来评定模型的精度。利用球面小波模型可以更加清晰地表示速度场的大尺度特征和复杂的局部变化特征。站点稠密区域,模型在东方向、北方向上的精度分别为±0.95mm/a、±0.97mm/a,稀疏区域对应的精度分别为±1.32mm/a和±1.30mm/a。     

Integrating GPS and GLONASS to accelerate convergence and initialization times of precise point positioning

The main challenge of dual-frequency precise point positioning (PPP) is that it requires about 30 min to obtain centimeter-level accuracy or to succeed in the first ambiguity-fixing. Currently, PPP is generally conducted with GPS only using the ionosphere-free combination. We adopt a single-differenced (SD) between-satellite PPP model to combine the GPS and GLONASS raw dual-frequency carrier phase measurements, in which the GPS satellite with the highest elevation is selected as the reference satellite to form the SD between-satellite measurements. We use a 7-day data set from 178 IGS stations to investigate the contribution of GLONASS observations to both ambiguity-float and ambiguity-fixed SD PPP solutions, in both kinematic and static modes. In ambiguity-fixed PPP, we only attempt to fix GPS integer ambiguities, leaving GLONASS ambiguities as float values. Numerous experimental results show that PPP with GLONASS and GPS requires much less convergence time than that of PPP with GPS alone. For ambiguity-float PPP, the average convergence time can be reduced by 45.9 % from 22.9 to 12.4 min in static mode and by 57.9 % from 40.6 to 17.7 min in kinematic mode, respectively. For ambiguity-fixed PPP, the average time to the first-fixed solution can be reduced by 27.4 % from 21.6 to 15.7 min in static mode and by 42.0 % from 34.4 to 20.0 min in kinematic mode, respectively. Experimental results also show that the less the GPS satellites are used in float PPP, the more significant is the reduction in convergence time when adding GLONASS observations. In addition, on average, more than 4 GLONASS satellites can be observed for most 2-h observation sessions. Nearly, the same improvement in convergence time reduction is achieved for those observations.     

GPS-PWV estimation and validation with radiosonde data and numerical weather prediction model in Antarctica

Three permanent GPS tracking stations in the trans Antarctic mountain deformation (TAMDEF) network were used to estimate precipitable water vapor (PWV) using measurement series covering the period of 2002–2005. TAMDEF is a National Science Foundation funded joint project between The Ohio State University and the United States Geological Survey. The TAMDEF sites with the longest GPS data spans considered in this research are Franklin Island East (FIE0), the International GNSS Service site McMurdo (MCM4), and Cape Roberts (ROB1). For the experiment, PWV was extracted from the ionosphere-free double-difference carrier phase observations, processed using the adjustment of GPS ephemerides (PAGES) software. The GPS data were processed with a 30 s sampling rate, 15-degree cutoff angle, and precise GPS orbits disseminated by IGS. The time-varying part of the zenith wet delay is estimated using the Marini mapping function, while the constant part is evaluated using the corresponding Marini tropospheric model. Previous studies using TAMDEF data for PWV estimation show that the Marini mapping function performs the best among the models offered by PAGES. The data reduction to compute the zenith wet delay follows the step piecewise linear strategy, which is subsequently transformed to PWV. The resulting GPS-based PWV is compared to the radiosonde observations and to values obtained from the Antarctic mesoscale prediction system (AMPS). This comparison revealed a consistent bias of 1.7 mm between the GPS solution and the radiosonde and AMPS reference values.     

Influence of clock jump on the velocity and acceleration estimation with a single GPS receiver based on carrier-phase-derived Doppler

Since the Selective Availability was turned off, the velocity and acceleration can be determined accurately with a single GPS receiver using raw Doppler measurements. The carrier-phase-derived Doppler measurements are normally used to determine velocity and acceleration when there is no direct output of the raw Doppler observations in GPS receivers. Due to GPS receiver clock drifts, however, a GPS receiver clock jump occurs when the GPS receiver clock resets itself (typically with 1 ms increment/decrement) to synchronize with the GPS time. The clock jump affects the corresponding relationship between measurements and their time tags, which results in non-equidistant measurement sampling in time or incorrect time tags. This in turn affects velocity and acceleration determined for a GPS receiver by the conventional method which needs equidistant carrier phases to construct the derived Doppler measurements. To overcome this problem, an improved method that takes into account, GPS receiver clock jumps are devised to generate non-equidistant-derived Doppler observations based on non-equidistant carrier phases. Test results for static and kinematic receivers, which are obtained by using the conventional method without reconstructing the equidistant continuous carrier phases, show that receiver velocity and acceleration suffered significantly from clock jumps. An airborne kinematic experiment shows that the greatest impact on velocity and acceleration reaches up to 0.2 m/s, 0.1 m/s² for the horizontal component and 0.5 m/s, 0.25 m/s² for the vertical component. Therefore, it can be demonstrated that velocity and acceleration measurements by using a standalone GPS receiver can be immune to the influence of GPS receiver clock jumps with the proposed method.     

Plate kinematics of Nubia–Somalia using a combined DORIS and GPS solution

We have used up to 12 years of data to assess DORIS performance for geodynamics applications. We first examine the noise characteristics of the DORIS time-series of weekly station coordinates to derive realistic estimates of velocity uncertainties. We find that a combination of white and flicker noise best explains the DORIS time-series noise characteristics. Second, weekly solutions produced by the Institut Géographique National/Jet Propulsion Laboratory (IGN/JPL) DORIS Analysis Centre are combined to derive a global velocity field. This solution is combined with two independent GPS solutions, including 11 sites on Nubia and 5 on the Somalia plate. The combination indicates that DORIS horizontal velocities have an average accuracy of 3 mm/year, with best-determined sites having velocity accuracy better than 1 mm/year (one-sigma levels). Using our combined velocity field, we derive an updated plate kinematics model with a focus on the Nubia–Somalia area. Including DORIS data improves the precision of the angular velocity vector for Nubia by 15%. Our proposed model provides robust bounds on the maximum opening rates along the East African Rift (4.7–6.7 mm/year). It indicates opening rates 15 and 7% slower than values predicted by NUVEL-1A for the southern Atlantic Ocean and Indian Ocean, respectively. These differences are likely to arise from the fact that NUVEL-1A considered Africa as a single non-deforming plate, while here we use a more refined approach.     

基于GPS时间序列分析的地震前后噪声特性差异性研究

针对现有地震预测研究中噪声特性分析研究较少,就较大地震对GPS时间序列噪声特性的影响进行分析.选取较大地震及相应的临近测站,将GPS时间序列按照时间跨度分为震前4～3 a、前3～2 a、前1～0 a、震后0～1 a、后1～2 a五个时段,分别求取各时段的谱指数,采用不同的噪声模型组合对它们进行噪声特性分析,然后对结果进...