低轨道卫星GPS精密定轨动力学方法的研究

结合低轨卫星简化动力学定轨算法,以及不同几何信息精度条件下的纯几何定轨和动力定轨精度比较,定量分析星载双频GPS实现精密定轨过程中的主要因素,得到星载GPS接收机性能设计所需的关键技术指标,为卫星精密定轨系统的顶层设计提供了科学合理的参考依据。     

Stochastic modelling considering ionospheric scintillation effects on GNSS relative and point positioning

Global Navigation Satellite Systems (GNSS), in particular the Global Positioning System (GPS), have been widely used for high accuracy geodetic positioning. The Least Squares functional models related to the GNSS observables have been more extensively studied than the corresponding stochastic models, given that the development of the latter is significantly more complex. As a result, a simplified stochastic model is often used in GNSS positioning, which assumes that all the GNSS observables are statistically independent and of the same quality, i.e. a similar variance is assigned indiscriminately to all of the measurements. However, the definition of the stochastic model may be approached from a more detailed perspective, considering specific effects affecting each observable individually, as for example the effects of ionospheric scintillation. These effects relate to phase and amplitude fluctuations in the satellites signals that occur due to diffraction on electron density irregularities in the ionosphere and are particularly relevant at equatorial and high latitude regions, especially during periods of high solar activity. As a consequence, degraded measurement quality and poorer positioning accuracy may result.     

Eliminating negative VTEC in global ionosphere maps using inequality-constrained least squares

Currently, ground-based Global Navigation Satellite System (GNSS) stations of the International GNSS Service (IGS) are distributed unevenly around the world. Most of them are located on the mainland, while only a small part of them are scattered on some islands in the oceans. As a consequence, many unreasonable zero values (in fact negative values) appear in Vertical Total Electron Content (VTEC) of European Space Agency (ESA) and Center for Orbit Determination in Europe (CODE) IONEX products, especially in 2008 and 2009 when the solar activities were rather quiet. To improve this situation, we directly implement non-negative physical constraints of ionosphere for global ionosphere maps (GIM) with spherical harmonic functions. Mathematically, we propose an inequality-constrained least squares method by imposing non-negative inequality constraints in the areas where negative VTEC values may occur to reconstruct GIM models. We then apply the new method to process the IGS data in 2008. The results have shown that the new algorithm efficiently eliminates the unwanted behavior of negative VTEC values, which could otherwise often be seen in the current CODE and ESA GIM products in both middle and high latitude areas of the Southern Hemisphere (45°S∼90°S) and the Northern Hemisphere (50°N∼90°N). About 64% of GPS receivers’ DCBs have been significantly improved. Finally, we compare the GIM results between with and without the inequality constraints, which has clearly shown that the GIM result with inequality constraints is significantly better than that without the inequality constraints. The inequality-constrained GIM result is also highly consistent with the final IGS products in terms of root mean squared (RMS) and mean VTEC.     

An inter-comparison of zenith tropospheric delays derived from DORIS and GPS data

Doppler Orbitography Radiopositioning Integrated by Satellite (DORIS) and Global Positioning System (GPS) techniques are similarly affected by propagation delays in the neutral atmosphere (troposphere) and hence make use of similar data processing strategies for reducing this effect. We compare Zenith Tropospheric Delays (ZTDs) estimated from 52 DORIS and GPS station pairs co-located at 35 sites over the 2005–2008 period. We find an overall systematic negative mean bias of −4 mm and a median bias of −2 mm, with a large site-to-site scatter and especially stronger biases over South America, potentially linked to remaining problems related to the South Atlantic Anomaly (SAA) in the current DORIS data processing. The standard deviation of ZTD differences is in the range 4–12 mm over the globe (8 mm on average), with larger values located in the southern hemisphere. The spatial variability of differences is consistent with previous work but remains largely unexplained. DORIS is shown to be much less sensitive to instrumental changes than GPS (only the switch from Alcatel to Starec antenna at Toulouse is detected as an offset of −4 mm in the ZTD time series). On the opposite, discontinuities and spurious annual signals are found in the GPS ZTD solutions. A discontinuity of +5 mm is found on 5 November 2006, linked to the switch from relative to absolute GPS antenna models used in the data processing. The use of modified GPS antennas (e.g. at GODE) or improved antenna models is shown to reduce the spurious annual signal (e.g. from 5 mm to 2 mm at METS). Overall, the agreement between both techniques is good, though DORIS shows a significantly larger random scatter. The high stability and good spatial and temporal coverage make DORIS a potential candidate technique for meteorology and climate studies as long as reasonable time averaging can be applied (e.g. differences are reduced from 8.6 to 2.4 mm with 5-day averages) and no real-time application is considered. This technique could be considered as a potential contributor to Global Geodetic Observing System (GGOS) for climatology.     

Ionospheric scintillations at Guilin detected by GPS ground-based and radio occultation observations

The occurrence of ionospheric scintillations with S4 ? 0.2 was studied using GPS measurements at Guilin, China (25.29°N, 110.33°E; geomagnetic: 15.04°N, 181.98°E), a station located near the northern crest of the equatorial anomaly. The results are presented for data collected from January 2009 to March 2010. The results show that nighttime amplitude scintillations only took place in February and March of the considered years, while daytime amplitude scintillations occurred in August and December of 2009. Nighttime amplitude scintillations, observed in the south of Guilin, always occurred with phase scintillations, TEC (Total Electron Content) depletions, and ROT (Rate Of change of TEC) fluctuations. However, TEC depletions and ROT fluctuations were weak during daytime amplitude scintillations, and daytime amplitude scintillations always took place simultaneously for most of the GPS satellites which appeared over Guilin in different azimuth directions. Ground-based GPS scintillation/TEC observations recorded at Guilin and signal-to-noise-ratio (SNR) measurements obtained from GPS-COSMIC radio occultation indicate that nighttime and daytime scintillations are very likely caused by ionospheric F region irregularities and sporadic E, respectively. Moreover, strong daytime amplitude scintillations may be associated with the plasma density enhancements in ionospheric E region caused by the Perseid and Geminid meteor shower activities.     

A global survey of COSMIC ionospheric peak electron density and its height: A comparison with ground-based ionosonde measurements

With a network of ground-based ionosondes distributed around the world, the ionospheric peak electron density and its height measured by FORMOSAT-3/COSMIC satellites in terms of GPS radio occultation technique are extensively examined in this article. It is found that, in spite of the latitude, the mean values of the peak electron density measured by COSMIC satellites are systematically smaller than those observed by ground-based ionosondes. The discrepancy between them is dependent on the latitude, namely, it is small in low and mid-latitudes and large in high-latitude region. Moreover, statistical analysis shows that the slopes of the regression line that is best fitted to the scatter diagram of occultation-retrieved peak electron density (ordinate axis) versus ionosonde-observed peak density (abscissa axis) are universally less than one. This feature is believed to be the result of path average effect of non-uniform distribution of the electron density along the GSP ray during the occultation. A comparison between COSMIC-measured peak height and ionosonde-derived peak height hmF2 indicates that the former is systematically higher than the latter. The difference in the two can be as large as 20% or more in equatorial and low-latitude regions. This result implies that the peak height hmF2 derived from the virtual height through true height analysis based on Titheridge method seems to underestimate the true peak height. The correlation between COSMIC and ionosonde peak electron densities is analyzed and the result reveals that correlation coefficient seems to be dependent on the fluctuation of the occultation-retrieved electron density profile. The correlation will be higher (lower) for the electron density profiles with smaller (larger) fluctuations. This feature suggests that the inhomogeneous distribution of the electron density along the GPS ray path during the occultation plays an important role affecting the correlation between COSMIC and ionosonde measurements.     

The geomagnetic solar flare effect identified by SIIG as an indicator of a solar flare observed by GOES satellites

This paper studies the efficiency of geomagnetic solar flare effects (gsfe) in X solar flare detection; so during the period 1999–2007 a comparison between solar flare (sf) observed by satellites of the Geostationary Operational Environmental Satellite (GOES) programme and gsfe published by the Service International des Indices Geomagnetiques (SIIG) is made.     

GPS姿态测量的载波相位整周模糊度快速解算

 根据 ARCE方法可以将模糊度搜索空间从所观测到的 m颗卫星的 m -1维降为 3维独立整周模糊度搜索空间,结合整数高斯变换及基线长度约束减小模糊度搜索空间,利用 Cholesky分解提高模糊度搜索效率,实验结果证明该方法能够快速解算整周模糊度适于实时姿态确定等应用。     

一种消除GPS模糊度相关性的新算法

 通过分析几种典型的整周模糊度去相关处理的方法 ,提出了一种简单快速消除GPS模糊度相关性的算法 ,介绍了算法中排序、整数LDLT 分解和模糊度筛选的处理过程 ,比较了该算法与其它算法性能上的差异 ,最后通过计算说明该算法的高效、可靠。     

NONLINEAR FILTERING ALGORITHM FOR INS INITIAL ALIGNMENT

The initial alignment error equation of an INS (Inertial Navigation System) with large initial azimuth error has been derived and nonlinear characteristics are included. When azimuth error is fairly small, the nonlinear equation can be reduced to a linear one. Extended Kalman filter, iterated filter and second order filter formulas are derived for the nonlinear state equation with linear measurement equation. Simulations results show that the accuracy of azimuth error estimation using extended Kalman filter is better than that of using standard Kalman filter while the iterated filter and second order filter can give even better estimation accuracy.