2006年12月13日太阳射电暴对GPS观测的影响

日地空间环境不仅影响航天器运行和安全, 也是导航、定位和通信等无线 电应用系统主要的误差源. 其中来自太阳L波段的射电暴被认为是全球导航卫星 系统(GNSS)稳定和性能的潜在威胁因素, 当L波段射电爆发达到一定阈值时, 将给用户带来不同程度的射电噪声干扰, 严重时会引起接收机失锁和定位服务 中断. 本文对2006年12月13日太阳射电暴对GPS造成的影响进行了研究, 利用太阳射电 观测数据、L波段闪烁观测数据和向阳面不同区域的GPS观测网数据, 分析 GPS观测对射电暴的响应. 结果表明, 此次事件对GPS观测产生了明显的影响, 射 电暴期间GPS发生幅度闪烁事件和明显失锁现象, 多个台站上空的多颗GPS 卫星 信号完全中断长达6min左右, 且多个台站上空锁定的卫星数目小于4颗, 使 得GPS定位完全失效. 相对而言, 射电暴期间日下点附近的GPS台站受到的影响 比远离日下点的大.     

航天器应急姿态捕获中GPS信号的应用

星载全球定位系统（GPS）卫星接收机在测量接收各GPS卫星信号时,可同时得到接收信号的信号强度测量辅助数据E。理论分析表明,接收信号的强度E与信号入射天线的法向夹角α强相关。如建立E与α稳定的先验模型,E就可以作为测量值,计算入射天线的角度α。在同一时刻,通过三个以上GPS卫星信号入射天线的角度α,可计算星载GPS卫星接收机接收天线的空间姿态。确定姿态的精度取决于E与α相关先验模型的稳定性。利用CHAMP卫星星载接收机在轨实测数据检验,估算的初始姿态精度为2°~3°。该方法可作为航天器故障状态下应急姿态捕获的一种辅助手段,也可为携带星载GPS而无高精度定姿要求的简易航天器提供一种新的无附加成本的定姿途径。     

A comprehensive study of the 2016 Mw 6.0 Italy earthquake based on high-rate (10 Hz) GPS data

The 2016 M_w 6.0 Italy earthquake is successfully recorded by the near-field 10?Hz GPS and 200?Hz Strong Motion (SM) stations, providing valuable data for this study. A comprehensive study of this earthquake is carried out based on GPS data, which contains coseismic deformations analysis, noise analysis, seismic wave picking, and magnitude determination. The noise of most GPS-derived displacement waveforms can be described as a combination of white noise, flicker noise, and random walk noise after the earthquake occurrence, and the spectral indices vary significantly for most stations, implying that the seismic signals have affected the noise characteristic of GPS-derived displacement waveforms. S-transform is employed to assess the GPS capability to detect the seismic arrival time. The SM station AMT and the GPS station AMAT are in good agreement in seismic wave picking, and the difference is only 1.2?s in the north component, suggesting that the outcome of seismic wave picking using GPS data is reliable. Then, a classic empirical formula is employed to determine the moment magnitude. A robust moment magnitude (Mw 5.90) can be estimated by the nine GPS stations with about 23.9?s. If four GPS stations near the epicenter is chosen to determine the magnitude, it only take 13.0?s to retrieve a reliable preliminary (Mw 5.82) magnitude, which is 5.4?s ahead of nine stations. In addition, Cross Wavelet Transform (XWT) is adopted to measuring the correlation and phase relationship between GPS and SM records. The result of XWT analysis indicates 10?Hz GPS is capable of capturing reliable and accurate coseismic dynamic deformations, as evidenced by the XWT-based semblance being close to 1 between GPS and SM records. The above results confirm the capability of 10?Hz GPS to capture coseismic dynamic deformations, detect seismic arrival time, and determine earthquake magnitude. Moreover, rapid magnitude determination based on 10?Hz GPS data can be regarded as an important supplement to Earthquake Early Warning (EEW).     

精密GPS卫星钟差的改正和应用

分析了GPS卫生钟差的变化特性，探讨了利用GPS地面跟踪站的观测数据估算GPS卫星钟差的可行性，建立了相应的算法和软件系统，并把由地面跟踪站的实测数据估算的卫星钟差用于星载GPS定轨计算，得到优于1m的定轨精度。     

南极GPS定位中的电离层延迟分析

基于南极地区国际GNSS服务组织(IGS,International GNSS Service)跟踪站的全球定位系统(GPS,Global Position System)双频实测数据,分析了南极地区电离层延迟的变化情况及其二阶项延迟对南极GPS定位结果的影响.结果表明:南极地区的总电子含量(TEC,Total Electron Content)日间波动频繁,其日间TEC最大值变化较中纬度地区剧烈;在南极地区夏季,电离层二阶项延迟对GPS定位结果的影响可达cm级.同时,由于欧洲定轨中心(CODE,The European Center for Orbit Determination)提供的全球电离层模型(GIM,Global Ionosphere Maps)在南极区域应用的局限性,通过选取南极地区6个IGS跟踪站作为基准站建立了区域电离层TEC模型(RIM,Regional Ionosphere Model).经实测数据计算证明,对于南极地区,RIM的定位精度在一定程度上优于全球电离层模型GIM.     

GPS data processing of networks with mixed single- and dual-frequency receivers for deformation monitoring

In order to obtain crustal deformations of higher spatial resolution, existing GPS networks must be densified. This densification can be carried out using single-frequency receivers at moderate costs. However, ionospheric delay handling is required in the data processing. We adapt the Satellite-specific Epoch-differenced Ionospheric Delay model (SEID) for GPS networks with mixed single- and dual-frequency receivers. The SEID model is modified to utilize the observations from the three nearest dual-frequency reference stations in order to avoid contaminations from more remote stations. As data of only three stations are used, an efficient missing data constructing approach with polynomial fitting is implemented to minimize data losses. Data from large scale reference networks extended with single-frequency receivers can now be processed, based on the adapted SEID model. A new data processing scheme is developed in order to make use of existing GPS data processing software packages without any modifications. This processing scheme is evaluated using a sub-network of the German SAPOS network. The results verify that the new scheme provides an efficient way to densify existing GPS networks with single-frequency receivers.     

不同太阳活动条件下电离层形态对估算GPS系统硬件延迟的影响

利用两个中纬度台站GPS观测数据提取的GPS卫星硬件延迟,分析了不同太阳活动情况下估算的硬件延迟稳定性和统计特征,结合同期电离层观测数据,研究了电离层状态对硬件延迟估算结果的影响.研究结果表明,基于太阳活动高年(2001年)GPS观测数据估算的硬件延迟稳定性要低于太阳活动低年GPS观测数据的估算结果,利用2001年GPS数据估算的卫星硬件延迟年标准偏差(RMS)平均值约为1TECU,而2009年GPS数据估算的卫星硬件延迟年标准偏差平均值约为0.8TECU.通过对2001年和2009年北京地区电离层F₂层最大电子密度(N_mF₂)变化性分析,结合GPS硬件延迟估算方法对电离层时空变化条件的要求,认为硬件延迟稳定性与太阳活动强度的联系是由不同太阳活动条件下电离层变化的强度差异引起的.     

基于GPS伪距的微小卫星定轨技术研究

高可靠、高精度、低成本的全自主导航是现代微小卫星定轨领域的发展方向,GPS以其在成本、精度、功耗等方面的独特优势,在低轨航天器上得到应用。论文采用GPS伪距作为观测量,在可见星较少的情况下仍可得到观测信息,较GPS直接输出提高了导航精度,并能进一步降低硬件成本。还引入了衰减记忆滤波算法,通过调节协方差矩阵,削弱了历史量测信息对当前数据融合的贡献,从而减小了轨道摄动对定轨精度的影响。仿真结果表明当导航星数目较少时,该算法仍可实现信息融合,且显著提高了定位精度。     

电离层周日变化对解算GPS硬件延迟稳定性的影响

针对电离层周日变化特征分析了其可能对SCORE方法估算的硬件延迟稳定性的影响. 利用BJFS以及XIAM台站的GPS观测数据, 解算了位于太阳活动高年(2001年)和太阳活动低年(2009年)的卫星硬件延迟并分析了估算的硬件延迟的稳定性. 研究发现, 电离层周日变化对估算的硬件延迟稳定性具有一定影响, 但是利用不同台站所得到的卫星硬件延迟稳定性在昼夜不同时间上的解算结果存在一定差异. 电离层周日变化对利用 BJFS台站数据解算的硬件延迟稳定性日夜差异较为明显, 在太阳活动高年利用XIAM 台站数据解算的硬件延迟日夜稳定性差异不很明显, 由于XIAM台站处于电离层赤道异常峰附近, 夜间电离层变化很大, 因此对比中纬度地区, 电离层周日变化对赤道异常峰附近地区硬件延迟稳定性解算结果的影响相对较小, 但在太阳活动低年, 其影响仍较为显著.     

Multivariate bootstrapped relative positioning of spacecraft using GPS L1/Galileo E1 signals

GNSS-based precise relative positioning between spacecraft normally requires dual frequency observations, whereas attitude determination of the spacecraft, mainly due to the stronger model given by the a priori knowledge of the length and geometry of the baselines, can be performed precisely using only single frequency observations. When the Galileo signals will come available, the number of observations at the L1 frequency will increase as we will have a GPS and Galileo multi-constellation. Moreover the L1 observations of the Galileo system and modernized GPS are more precise than legacy GPS and this, combined with the increased number of observations, will result in a stronger model for single frequency relative positioning. In this contribution we will develop an even stronger model by combining the attitude determination problem with relative positioning. The attitude determination problem will be solved by the recently developed Multivariate Constrained (MC-) LAMBDA method. We will do this for each spacecraft and use the outcome for an ambiguity constrained solution on the baseline between the spacecraft. In this way the solution for the unconstrained baseline is bootstrapped from the MC-LAMBDA solutions of each spacecraft in what is called: multivariate bootstrapped relative positioning. The developed approach will be compared in simulations with relative positioning using a single antenna at each spacecraft (standard LAMBDA) and a vectorial bootstrapping approach. In the simulations we will analyze single epoch, single frequency success rates as the most challenging application. The difference in performance for the approaches for single epoch solutions, is a good indication of the strength of the underlying models. As the multivariate bootstrapping approach has a stronger model by applying information on the geometry of the constrained baselines, for applications with large observation noise and limited number of observations this will result in a better performance compared to the vectorial bootstrapping approach. Compared with standard LAMBDA, it can reach a 59% higher success rate for ambiguity resolution. The higher success rate on the unconstrained baseline between the platforms comes without extra computational load as the constrained baseline(s) problem has to be solved for attitude determination and this information can be applied for relative positioning.     

基于单频星载GPS数据的低轨卫星精密定轨

为满足搭载单频GPS接收机低轨卫星的精密定轨需求以及深化单频定轨研究,文中解决了单频星载GPS数据的周跳探测问题,并利用“海洋二号”（HY-2A）卫星及“资源三号”（ZY-3）卫星的单频星载GPS实测数据采用两种方法确定了二者的简化动力学轨道,并通过观测值残差分析、与双频精密轨道比较、激光测卫数据检核等方法对所得轨道精度进行评定。结果表明,在不考虑电离层延迟影响的情况下,HY-2A卫星定轨精度为2~3dm,ZY-3卫星为1m左右;而采用半和改正组合消除电离层延迟一阶项影响后,二者定轨精度均显著提高,HY-2A卫星三维精度提高至1dm左右,ZY-3卫星提高至1~2dm。文章的研究成果表明,搭载单频GPS接收机的低轨卫星也可获得厘米级的定轨精度。     

A regional tropospheric model for airborne GPS applications

A regional tropospheric model can be constructed using surveys from GPS ground networks. Using this model the tropospheric delays of a kinematic station within the region can be interpolated. However, such a model is generally not suitable for an airborne platform high above the ground networks. In this paper, a method of constructing a regional tropospheric model for airborne GPS applications is described. First, the kinematic station in the air is projected onto the ground. Then the tropospheric delays at projected point are interpolated from those of the ground networks. Finally, the tropospheric delays at projected point are extended upward to the airborne platform using pressure and temperature gradients and humidity exponential function. For validation of this method, the data of airborne campaign carried out by BKG in cooperation with GFZ and BGR 2008 in the northern Alps are used. The results show that GPS kinematic positioning precision in height component can be improved using this method.     

Epochwise prediction of GPS single differenced ionospheric delays of formation flying spacecraft

Both single and dual frequency GPS relative navigation filters may benefit from proper predictions of single differenced ionospheric delays. In this article, the single differenced ionospheric delays of GPS observations are predicted for the GRACE formation during the switch manoeuvre.Two prediction methods are considered. The first is based on a Taylor expansion to first order of a mapping function that maps slant total electron content measurements to vertical total electron content estimates. The second method fits a shape profile through undifferenced ionospheric data available. It then raytraces through this profile to estimate the difference in total electron content along the path of the GPS signals.Continuously changing ionospheric conditions hamper the assessment of the quality of the predictions. Comparison of both methods shows that the raytracing method performs better. The difference of predictions and measurements generally shows a smaller RMS than the measurements alone. However, both methods suffer from a number of systematically unpredicted observations, which arise from small unaccounted differential variations in electron densities along the path of the GPS signals. These prediction methods perform better when spacecraft separation is small. Baselines considered here range from tens of kilometres down to several hundred metres. When smallest spacecraft separation occurs (0.4 km), the single differenced ionospheric delays exhibit RMS values of 0.0089 m. The first method shows a difference between measurements and predictions with an RMS of 0.0081 m. For the second method the difference RMS is found to be 0.0067 m.     

GPS sensing of precipitable water vapour during the March 2010 Melbourne storm

The March 2010 Melbourne storm is used as a case study to examine the potential of using Global Positioning System (GPS) observations for studying the precipitable water vapour (PWV) field. The Victorian statewide GPS infrastructure network, i.e. GPSnet, was used in this study. GPSnet is currently the only statewide and densest GPS infrastructure network in Australia, which provides an excellent opportunity to examine the distribution of water vapour as the severe weather system passed over the state. Data from 15 GPSnet stations were processed over a one-week period, i.e. a few days prior to and after the storm passage, during which the course of the storm extended from the west to the southeast corner of the state. In addition, data from two radiosonde sites of the Australian Bureau of Meteorology Upper Air Network were used to compare and validate the GPS derived PWV measurements. The findings demonstrate that there is strong spatial and temporal correlation between variations of the ground-based GPS-PWV estimates and the passage of the storm over the state. This is encouraging as the ground-based GPS water vapour sensing technique can be considered as a supplemental meteorological sensor in studying severe weather events. The advantage of using ground-based GPS-PWV technique is that it is capable of providing continuous observation of the storm passage with high temporal resolution. The spatial resolution of the distribution of water vapour is dependent on the geographical location and density of the GPS stations.     

Determination of crustal deformations based on GPS observing-session duration in Marmara region,Turkey

This paper reports on the estimation of changes of a GPS-based crustal deformation field as function of the observing session duration. This investigation was carried out on the Marmara Continuous GPS Network (Marmara region) by using GPS data collected on a set of 10 stations during the period 2002.50–2005.47 The GPS observations were processed in the ITRF 2005 reference frame using Bernese v5.0 software.     

Enhancing the kinematic precise orbit determination of low earth orbiters using GPS receiver clock modelling

Clock error estimation has been the focus of a great deal of research because of the extensive usage of clocks in GPS positioning applications. The receiver clock error in the spacecraft orbit determination is commonly estimated on an epoch-by-epoch basis, along with the spacecraft’s position. However, due to the high correlation between the spacecraft orbit altitude and the receiver clock parameters, estimates of the radial component are degraded in the kinematic approach. Using clocks with high stability, the predictable behaviour of the receiver oscillator can be exploited to improve the positioning accuracy, especially for the radial component. This paper introduces two GPS receiver clock models to describe the deterministic and stochastic property of the receiver clock, both of which can improve the accuracy of kinematic orbit determination for spacecraft in low earth orbit. In particular, the clock parameters are estimated as time offset and frequency offset in the two-state model. The frequency drift is also estimated as an unknown parameter in the three-state model. Additionally, residual non-deterministic random errors such as frequency white noise, frequency random walk noise and frequency random run noise are modelled. Test results indicate that the positioning accuracy could be improved significantly using one day of GRACE flight data. In particular, the error of the radial component was reduced by over 40.0% in the real-time scenario.     

Application of neural networks to South African GPS TEC modelling

The propagation of radio signals in the Earth’s atmosphere is dominantly affected by the ionosphere due to its dispersive nature. Global Positioning System (GPS) data provides relevant information that leads to the derivation of total electron content (TEC) which can be considered as the ionosphere’s measure of ionisation. This paper presents part of a feasibility study for the development of a Neural Network (NN) based model for the prediction of South African GPS derived TEC. The South African GPS receiver network is operated and maintained by the Chief Directorate Surveys and Mapping (CDSM) in Cape Town, South Africa. Vertical total electron content (VTEC) was calculated for four GPS receiver stations using the Adjusted Spherical Harmonic (ASHA) model. Factors that influence TEC were then identified and used to derive input parameters for the NN. The well established factors used are seasonal variation, diurnal variation, solar activity and magnetic activity. Comparison of diurnal predicted TEC values from both the NN model and the International Reference Ionosphere (IRI-2001) with GPS TEC revealed that the IRI provides more accurate predictions than the NN model during the spring equinoxes. However, on average the NN model predicts GPS TEC more accurately than the IRI model over the GPS locations considered within South Africa.     

Non-parametric regional VTEC modeling with Multivariate Adaptive Regression B-Splines

In this work Multivariate Adaptive Regression B-Splines (BMARS) is applied to regional spatio-temporal mapping of the Vertical Total Electron Content (VTEC) using ground based Global Positioning System (GPS) observations. BMARS is a non-parametric regression technique that utilizes compactly supported tensor product B-splines as basis functions, which are automatically obtained from the observations. The algorithm uses a scale-by-scale model building strategy that searches for B-splines at each scale fitting adequately to the data and provides smoother approximations than the original Multivariate Adaptive Regression Splines (MARS). It is capable to process high dimensional problems with large amounts of data and can easily be parallelized. The real test data is collected from 32 ground based GPS stations located in North America. The results are compared numerically and visually with both the regional VTEC modeling generated via original MARS using piecewise-linear basis functions and another regional VTEC modeling based on B-splines.     

Comparison of GPS derived TEC with the TEC predicted by IRI 2012 model in the southern Equatorial Ionization Anomaly crest within the Eastern Africa region

We have compared the TEC obtained from the IRI-2012 model with the GPS derived TEC data recorded within southern crest of the EIA in the Eastern Africa region using the monthly means of the 5 international quiet days for equinoxes and solstices months for the period of 2012 – 2013. GPS-derived TEC data have been obtained from the Africa array and IGS network of ground based dual-frequency GPS receivers from four stations (Kigali (1.95°S, 30.09°E; Geom. Lat. 11.63°S), Malindi (2.99°S, 40.19°E; Geom. Lat. 12.42°S), Mbarara (0.60°S, 30.74°E; Geom. Lat. 10.22°S) and Nairobi (1.22°S, 36.89°E; Geom. Lat. 10.69°S)) located within the EIA crest in this region. All the three options for topside Ne of IRI-2012 model and ABT-2009 for bottomside thickness have been used to compute the IRI TEC. Also URSI coefficients were considered in this study. These results are compared with the TEC estimated from GPS measurements. Correlation Coefficients between the two sets of data, the Root-Mean Square Errors (RMSE) of the IRI-TEC from the GPS-TEC, and the percentage RMSE of the IRI-TEC from the GPS-TEC have been computed. Our general results show that IRI-2012 model with all three options overestimates the GPS-TEC for all seasons and at all stations, and IRI-2001 overestimates GPS-TEC more compared with other options. IRI-Neq and IRI-01-corr are closely matching in most of the time. The observation also shows that, GPS TEC are underestimated by TEC from IRI model during noon hours, especially during equinoctial months. Further, GPS-TEC values and IRI-TEC values using all the three topside Ne options show very good correlation (above 0.8). On the other hand, the TEC using IRI-Neq and IRI-01- corr had smaller deviations from the GPS-TEC compared to the IRI-2001.     

一种基于星载GPS的月球探测器导航算法

为验证GNSS旁瓣信号对中高轨航天器在轨飞行的导航支持能力,"嫦娥5号"飞行试验器搭载了GPS导航接收机。由于目前使用的星载导航解平稳性不够,为此,根据月球探测器不同飞行段落轨道动力学特性,建立了一种基于星载GPS数据的月球探测器实时导航算法。该算法充分考虑了月球探测器不同段落的动力学特性,基于UKF滤波设计,并引入了自适应策略,以降低滤波参数选取的难度。在分析星载GPS数据特点基础上,月地返回段实测数据的处理结果验证了算法的可行性和有效性。