The adjusted optical properties for Galileo/BeiDou-2/QZS-1 satellites and initial results on BeiDou-3e and QZS-2 satellites

Solar Radiation Pressure (SRP) is the dominant non-gravitational perturbation for GNSS (Global Navigation Satellite System) satellites. In the absence of precise surface models, the Empirical CODE Orbit Models (ECOM, ECOM2) are widely used in GNSS satellite orbit determination. Based on previous studies, the use of an a priori box-wing model enhances the ECOM model, especially if the spacecraft is a stretched body satellite. However, so far not all the GNSS system providers have published their metadata. To ensure a precise use of the a priori box-wing model, we estimate the optical parameters of all the Galileo, BeiDou-2, and QZS-1 (Quasi Zenith Satellite System) satellites based on the physical processes from SRP to acceleration. Validation using orbit prediction proves that the adjusted parameters of Galileo and QZS-1 satellites exhibit almost the same performance as the corresponding published and “best guess” values. Whereas, the estimated parameters of BeiDou-2 satellites demonstrate an improvement of more than 60% over the initial “guess” values. The resulting optical parameters of all the satellites are introduced into an a priori box-wing model, which is jointly used with ECOM and ECOM2 model in the orbit determination. Results show that the pure ECOM2 model exhibits better performance than the pure ECOM model for Galileo, BeiDou-2 GEO and QZS-1 orbits. Combined with the a priori box-wing model the ECOM model (ECOM+BW) results in the best Galileo, BeiDou-2 GEO and QZS-1 orbits. The standard deviation (STD) of satellite laser ranging residuals reduce by about 20% and 5% with respect to the pure ECOM2 model for Galileo and BeiDou-2 GEO orbits, while the reductions are about 40% and 60% for QZS-1 orbits in yaw-steering and orbit-normal mode respectively. BeiDou-2 IGSO and MEO satellite orbits do not benefit much from the a priori box-wing model. In summary, we suggest setting up a unified SRP model of ECOM+BW for Galileo, QZS-1, and BeiDou-2 orbits based on the adjusted metadata. In addition, we estimate the optical parameters of BeiDou-3e and QZS-2 satellites using a limited number of tracking stations. Results regarding the unified SRP model indicate the same advantages, the STD of satellite laser ranging residuals reduces by about 30% and 20% for QZS-2 and BeiDou-3e orbits respectively over orbit products without a priori model. The estimation procedure is effective and easy to apply to the new emerging satellites in the future.     

Analysis on navigation accuracy of high orbit spacecraft based on BDS GEO and IGSO satellites

To make up for the insufficiency of earth-based TT&C systems, the use of GNSS technology for high-orbit spacecraft navigation and orbit determination has become a new technology. It is of great value to applying Geosynchronous Earth Orbit (GEO) and Inclined GeoStationary Orbit (IGSO) navigation satellites for supporting the navigation of high-orbit spacecraft since there are three different types of navigation satellites in BeiDou Navigation Satellite System (BDS): Medium Earth Orbit (MEO), GEO and IGSO. This paper conducts simulation experiments based on Two-Line Orbital Element (TLE) data to analyze and demonstrate the role of these satellites in the navigation of high-orbit spacecraft. Firstly, the spacecraft in GEO was used as the target satellite to conduct navigation experiments. Experiments show that for the spacecraft on the GEO orbit, after adding GEO and IGSO respectively on the basis of receiving MEO navigation satellite signals, the accuracies were improved by 7.22 % and 6.06 % respectively. When adding both GEO and IGSO navigation satellites at the same time, the accuracy can reach 16 m. In the second place, navigation and positioning experiments were carried out on three high elliptical orbit (HEO) satellites with different semimajor axis (32037.2 km, 42385.9 km, 67509.6 km). The experiments show that the number of visible satellites has been improved significantly after adding GEO and IGSO navigation satellites at the same time. The visible satellites in these three orbits were improved by 32.84 %, 41.12 % and 37.68 %, respectively compared with only observing MEO satellites.The RMS values of the navigation positioning errors of these three orbits are 25.59 m, 87.58 m and 712.48 m, respectively.     

Instantaneous BeiDou–GPS attitude determination: A performance analysis

The advent of modernized and new global navigation satellite systems (GNSS) has enhanced the availability of satellite based positioning, navigation, and timing (PNT) solutions. Specifically, it increases redundancy and yields operational back-up or independence in case of failure or unavailability of one system. Among existing GNSS, the Chinese BeiDou system (BDS) is being developed and will consist of geostationary (GEO) satellites, inclined geosynchronous orbit (IGSO) satellites, and medium-Earth-orbit (MEO) satellites. In this contribution, a BeiDou–GPS robustness analysis is carried out for instantaneous, unaided attitude determination.     

遥感卫星传感器的重复观测能力计算

遥感卫星传感器的重复观测能力是体现卫星系统动态监测能力的重要指标之一. 通过分析卫星星下点成像和侧视成像时的重访条件, 给出了卫星可能实现重访的轨道计算方法, 同时以HJ-1A/1B卫星(环境与灾害监测预报小卫星星座的光学星)为例, 给出了特定地区可以实现重访的优选轨道, 并讨论了由两颗卫星组成的卫星星座对重访的影响, 以及纬度对重访的影响.      

Orbit determination of BeiDou navigation satellite system maneuvered satellites based on instantaneous velocity pulses parameters

The BeiDou navigation satellite system (BDS) comprises geostationary earth orbit (GEO) satellites as well as inclined geosynchronous orbit (IGSO) and medium earth orbit (MEO) satellites. Owing to their special orbital characteristics, GEO satellites require frequent orbital maneuvers to ensure that they operate in a specific orbital window. The availability of the entire system is affected during the maneuver period because service cannot be provided before the ephemeris is restored. In this study, based on the conventional dynamic orbit determination method for navigation satellites, multiple sets of instantaneous velocity pulses parameters which belong to one of pseudo-stochastic parameters were used to simulate the orbital maneuver process in the orbital maneuver arc and establish the observed and predicted orbits of the maneuvered and non-maneuvered satellites of BeiDou regional navigation satellite system (BDS-2) and BeiDou global navigation satellite system (BDS-3). Finally, the single point positioning (SPP) technology was used to verify the accuracy of the observed and predicted orbits. The orbit determination accuracy of maneuvered satellites can be greatly improved by using the orbit determination method proposed in this paper. The overlapping orbit determination accuracy of maneuvered GEO satellites of BDS-2 and BDS-3 can improve 2–3 orders of magnitude. Among them, the radial orbit determination accuracy of each maneuvered satellite is basically better than 1 m. simultaneously, the combined orbit determination of the maneuvered and non-maneuvered satellites does not have a great impact on the orbit determination accuracy of the non-maneuvered satellites. Compared with the multi GNSS products (indicated by GBM) from the German Research Centre for Geosciences (GFZ), the impact of adding the maneuvered satellites on the orbit determination accuracy of BDS-2 satellites is less than 9 %. Furthermore, the orbital recovery time and the service availability period are significantly improved. When the node of the predicted orbit is traversed approximately 3 h after the maneuver, the accuracy of the predicted orbit of the maneuvered satellite can reach that of the observed orbit. The SPP results for the BDS reached a normal level when the node of the predicted orbit was 2 h after the maneuver.     

LEO enhanced Global Navigation Satellite System (LeGNSS) for real-time precise positioning services

Global Navigation Satellite System (GNSS) has been widely used in many geosciences areas with its Positioning, Navigation and Timing (PNT) service. However, GNSS still has its own bottleneck, such as the long initialization period of Precise Point Positioning (PPP) without dense reference network. Recently, the concept of PNTRC (Positioning, Navigation, Timing, Remote sensing and Communication) has been put forward, where Low Earth Orbit (LEO) satellite constellations are recruited to fulfill diverse missions. In navigation aspect, a number of selected LEO satellites can be equipped with a transmitter to transmit similar navigation signals to ground users, so that they can serve as GNSS satellites but with much faster geometric change to enhance GNSS capability, which is named as LEO constellation enhanced GNSS (LeGNSS). As a result, the initialization time of PPP is expected to be shortened to the level of a few minutes or even seconds depending on the number of the LEO satellites involved. In this article, we simulate all the relevant data from June 8th to 14th, 2014 and investigate the feasibility of LeGNSS with the concentration on the key issues in the whole data processing for providing real-time PPP service based on a system configuration with fourteen satellites of BeiDou Navigation Satellite System (BDS), twenty-four satellites of the Global Positioning System (GPS), and sixty-six satellites of the Iridium satellite constellations. At the server-end, Precise Orbit Determination (POD) and Precise Clock Estimation (PCE) with various operational modes are investigated using simulated observations. It is found out that GNSS POD with partial LEO satellites is the most practical mode of LeGNSS operation. At the user-end, the Geometry Dilution Of Precision (GDOP) and Signal-In-Space Ranging Error (SISRE) are calculated and assessed for different positioning schemes in order to demonstrate the performance of LeGNSS. Centimeter level SISRE can be achieved for LeGNSS.     

基于参考轨道的Walker星座相对相位保持策略

针对高轨道Walker全球星座相对相位保持问题,在分析星座几何结构演化规律基础上,提出了一种基于动态调整参考轨道的星座相对相位保持策略.该参考轨道的半长轴取星座所有卫星的平均轨道半长轴的平均值,通过计算各卫星相对参考轨道的相对相位偏差以及变化率,而后根据有关维持的约束条件来选择合适的卫星进行站点保持,使得所有卫星的相对参考轨道的相位偏差不超过允许值.当一旦有卫星进行了站点保持,则参考轨道就重新统计确定.文章通过一个12颗MEO卫星构成的Walker-δ星座的分析算例,表明这种相对保持策略可以减少星座的站点维持次数.本文提出的方法可以为我国今后卫星导航星座的维持提供参考.     

An empirical solar radiation pressure model for satellites moving in the orbit-normal mode

We present a family of empirical solar radiation pressure (SRP) models suited for satellites orbiting the Earth in the orbit normal (ON) mode. The proposed ECOM-TB model describes the SRP accelerations in the so-called terminator coordinate system. The choice of the coordinate system and the SRP parametrization is based on theoretical assumptions and on simulation results with a QZS-1-like box-wing model, where the SRP accelerations acting on the solar panels and on the box are assessed separately. The new SRP model takes into account that in ON-mode the incident angle of the solar radiation on the solar panels is not constant like in the yaw-steering (YS) attitude mode. It depends on the elevation angle of the Sun above the satellite’s orbital plane. The resulting SRP vector acts, therefore, not only in the Sun-satellite direction, but has also a component normal to it. Both components are changing as a function of the incident angle. ECOM-TB has been used for precise orbit determination (POD) for QZS-1 and BeiDou2 (BDS2) satellites in medium (MEO) and inclined geosynchronous Earth orbits (IGSO) based on IGS MGEX data from 2014 and 2015. The resulting orbits have been validated with SLR, long-arc orbit fits, orbit misclosures, and by the satellite clock corrections based on the orbits. The validation results confirm that—compared to ECOM2—ECOM-TB significantly (factor 3–4) improves the POD of QZS-1 in ON-mode for orbits with different arc lengths (one, three, and five days). Moderate orbit improvements are achieved for BDS2 MEO satellites—especially if ECOM-TB is supported by pseudo-stochastic pulses (the model is then called ECOM-TBP). For BDS2 IGSOs, ECOM-TB with its 9 SRP parameters appears to be over-parameterized. For use with BDS2 IGSO spacecraft we therefore developed a minimized model version called ECOM-TBMP, which is based on the same axis decomposition as ECOM-TB, but has only 2 SRP parameters and is supported by pseudo-stochastic parameters, as well. This model shows a similar performance as ECOM-TB with short arcs, but an improved performance with (3-day) long-arcs. The new SRP models have been activated in CODE’s IGS MGEX solution in Summer 2018. Like the other ECOM models the ECOM-TB derivatives might be used together with an a priori model.     

Galileo real-time orbit determination with multi-frequency raw observations

Real-time GNSS-based applications require corresponding real-time orbit products. While traditional GNSS orbits are generated with the dual-frequency IF (Ionosphere-Free) model, the increase of multi-frequency signal satellites brings new challenges for the data processing. Therefore, real-time orbit determination with the multi-frequency UC (Uncombined) model is introduced in this study considering its flexibility. With the derived mathematical model conforming to IGS (International GNSS Service) dual-frequency clock definition and one-week triple-frequency Galileo observation data from 90 IGS network stations, the convergence and accuracy of real-time orbits is assessed and the characteristics of satellite IFCB (Inter-Frequency Clock Bias) are analyzed. Results indicate that the model differences, including dual-frequency IF model, dual-frequency UC model and triple-frequency UC model, contribute to only cm-level differences with CODE (Center for Orbit Determination in Europe) final orbits after a convergence time of around 12 h. The constellation-mean RMS (Root Mean Square) differences of the converged real-time orbits with the CODE final orbits reaches about 5.0 cm, 7.0 cm and 5.0 cm for the radial, tangential and normal directions. The convergence of satellite IFCB is much faster than that of satellite orbit, which reflects a loose correlation between these two parameters. While the Galileo satellite IFCB are temporally stable, the modeling of satellite IFCB may be unreliable when over constrained and becomes even more unstable with commonly encountered datum changes. In summary, real-time GNSS orbit determination with multi-frequency raw observations is feasible and extendable with proper treatment of IFCB.     

Near real-time BDS GEO satellite orbit determination and maneuver analysis with reversed point positioning

The Geostationary Earth Orbit (GEO) satellite is a crucial part of the BeiDou Navigation Satellite System (BDS) constellation. However, due to various perturbation forces acting on the GEO satellite, it drifts gradually over time. Thus, frequent orbit maneuvers are required to maintain the satellite at its designed position. During the orbit maneuver and recovery periods, the orbit quality of the maneuvered satellite computed with broadcast navigation ephemeris will be significantly degraded. Furthermore, the conventional dynamic Precise Orbit Determination (POD) approach may not work well, because of a lack of publicly available satellite information for modeling the thrust forces. In this paper, a near real-time approach free of thrust forces modeling is proposed for BDS GEO satellite orbit determination and maneuver analysis based on the Reversed Point Positioning (RPP). First, the station coordinates and receiver clock offsets are estimated by GPS/BDS combined Single Point Positioning (SPP) with single-frequency phase-smoothed pseudorange observations. Then, with the fixed station coordinates and receiver clock offsets, the RPP method can be conducted to determine the GEO satellite orbits. When no orbit maneuvers occur, the proposed method can obtain orbit accuracies of 0.92, 2.74, and 8.30?m in the radial, along-track, and cross-track directions, respectively. The average orbit-only Signal-In-Space Range Error (SISRE) is 1.23?m, which is slightly poorer than that of the broadcast navigation ephemeris. Using four days of GEO maneuvered datasets, it is further demonstrated that the derived orbits can be employed to characterize the behaviors of GEO satellite maneuvers, such as the time span of the maneuver as well as the satellite thrusting accelerations. These results prove the efficiency of the proposed method for near real-time GEO satellite orbit determination during maneuvers.     

Early analysis of precise orbit and clock offset determination for the satellites of the global BeiDou-3 system

Eight new-generation BeiDou satellites (BeiDou-3) have been launched into Medium Earth Orbit (MEO), allowing for global coverage since March 2018, and they are equipped with new hydrogen atomic clocks and updated rubidium clocks. Firstly, we analyzed the signals for the carrier-to-noise-density ratio (C/N0) and pseudorange multipath (MP) by using international GNSS (Global Navigation Satellite System) Monitoring and Assessment System (iGMAS) station data, and found that B1C has a lower C/N0, and B2a has the same level of C/N0 as the B1I and B3I signals. For pseudorange multipath, compared with the BeiDou-2 satellites, the obvious systematic variation of MP scatters related to the elevation angle is greatly improved for the BeiDou-3 and BeiDou-3e satellites signals. For the signals of the BeiDou-3 satellites, the order of the Root Mean Square (RMS) values of multipath and noise is B3I?<?B1I?<?B2a?<?B1C. Then, the comparison of the precise orbit determination and clock offset determination for the BeiDou-2, BeiDou-3, and BeiDou-3 experimental (BeiDou-3e) satellites was done by using 10 stations from iGMAS. The 3D precision of the 24?h orbit overlap is 24.55, 25.61, and 23.35?cm for the BeiDou-3, BeiDou-3e, and BeiDou-2 satellites, respectively. BeiDou-3 satellite has a comparable precision to that of the BeiDou-2 satellite. For the precision of clock offset estimation, the Standard Deviation (STD) of the BeiDou-3 MEO satellite is 0.350?ns, which is an improvement of 0.042?ns over that of the BeiDou-2 MEO satellite. The stabilities of the BeiDou-3 and BeiDou-3e onboard clocks are better than those of BeiDou-2 by factors of 2.84 and 1.61 at an averaging time of 1000 and 10,000?s, respectively.     

基于射频的导航卫星星间链路信息传输与仿真研究

星间链路是导航卫星实现精密定轨和自主导航的关键技术之一。导航卫星通过星间链路完成伪距测量和数据交换,维持系统稳定运行的时空基准,保证系统持续提供精准导航服务。根据全球导航卫星系统的建设情况和发展趋势,首先介绍星间观测和信息传输频段,并从天线特征、多址控制方式和网络拓扑结构等角度分析了射频链路的工作体制。最后,针对实际导航卫星星座,应用OPNET平台建立导航信息传输仿真模型,通过分析信息传输效率,验证了基于射频链路导航信息传输的可行性和有效性,对全球导航卫星系统的星间链路研究具有一定的参考价值。     

风云三号C星GNOS北斗掩星电离层探测初步结果

利用风云三号卫星C星GNOS掩星探测仪电离层数据,分析了2013年10月FY-3C GNOS探测的北斗掩星电离层廓线分布,将2013年10月1日至2015年10月10日期间FY-3C GNOS观测的F₂层峰值电子密度（N_mF₂）与地面电离层测高仪观测结果进行对比,验证了FY-3C GNOS北斗电离层掩星的探测精度.结果表明,FY3-C GNOS北斗电离层掩星与电离层测高仪探测的N_mF₂数据相关系数为0.96,平均偏差为10.21%,标准差为19.61%.在不同情况下其数据精度有如下特征:白天精度高于夜晚;夏季精度高于分季,分季精度高于冬季;中纬地区精度高于低纬地区,低纬地区精度高于高纬地区; BDS倾斜同步轨道（IGSO）卫星精度高于同步轨道（GEO）卫星和中轨道（MEO）卫星.FY-3C GNOS北斗电离层掩星与国际上其他掩星电离层数据精度的一致性对GNSS掩星探测资料的综合利用具有重大意义.      

基于射频的导航卫星星间链路信息传输与仿真研究

星间链路是导航卫星实现精密定轨和自主导航的关键技术之一。导航卫星通过星间链路完成伪距测量和数据交换,维持系统稳定运行的时空基准,保证系统持续提供精准导航服务。根据全球导航卫星系统的建设情况和发展趋势,首先介绍星间观测和信息传输频段,并从天线特征、多址控制方式和网络拓扑结构等角度分析了射频链路的工作体制。最后,针对实际导航卫星星座,应用OPNET平台建立导航信息传输仿真模型,通过分析信息传输效率,验证了基于射频链路导航信息传输的可行性和有效性,对全球导航卫星系统的星间链路研究具有一定的参考价值。     

一种双基站低轨道区域覆盖δ星座设计

在双基站理论研究基础上 ,进行了定倾角区域性覆盖双基站δ星座设计 ,采用倾角为 41°的回归轨道 ,用点覆盖数值仿真的方法设计了对地面站与空间站的一重覆盖与二重覆盖下的星座构型。且星座中各星之间 ,各星与空间站之间均可建立星际链路 ,保证数据在空间站、地面站及星座之间连续传递     

Kalman filter based estimation of differential hardware biases with triangular interpolation technique for IRNSS

Ionosphere delay correction is the main error correction to the computation of single frequency user position using satellite navigation. However ionosphere delay consists of not only delay but also frequency dependent differential hardware biases from satellite and receiver ends. For ionosphere point of view, Indian Regional Navigation Satellite System (IRNSS) service area comes in equatorial anomaly region. It is a unique satellite navigation system which operates at L5 and S frequencies and consists of Geostationary Earth Orbit (GEO) and Geo Synchronous Orbit (GSO) satellite constellation. With IRNSS measurements availability, there is a good opportunity to estimate and analyse differential hardware biases with GEO/GSO combination and with equatorial ionosphere variation. In this paper, Kalman filter based estimation with triangular interpolation technique is used to estimate differential hardware biases for all IRNSS satellites and reference receivers at L5 frequency. The standard deviation of the 15?days of daily estimation of satellite differential hardware biases is in the range of 0.32 to 1.17 TECU for all IRNSS satellites. Similarly, the standard deviation of the 15?days of daily estimation varies up to 2.85 and 6.0 TECU for receiver differential hardware biases during calm and stormy period respectively. The ionosphere delay computed using estimated differential hardware biases is compared with Global Ionosphere Map (GIM) data. A rigorous analysis is carried out to study the error in the estimation in terms of input data noise level, satellite constellation and effect of latitude. Our result reveals that over IRNSS service area, there is an exponential increase in the error in the estimation of receiver differential hardware biases with respect to latitude.     

混合星座导航卫星广播星历相关问题研究

GPS卫星广播星历参数具有参数少、物理意义明确以及精度高等特点,可以考虑将它应用于包含MEO、IGSO和GEO卫星的混合星座卫星导航系统。分析了采用GPs卫星广播星历参数时MEO、IGSO和GEO卫星的广播星历拟合精度,并且比较分析了在一个卫星的轨道周期内,广播星历参数拟合结果的变化规律。仿真结果表明,MEO、IGSO和GEO卫星的广播星历拟合误差最大在分米量级;MEO和IGSO卫星在一个轨道周期内星历参数拟合结果的变化规律相近,但是与GEO卫星的差异较大。     

The impact on tsunami detection from space using GNSS-reflectometry when combining GPS with GLONASS and Galileo

The devastating Sumatra tsunami in 2004 demonstrated the need for a tsunami early warning system in the Indian Ocean. Such a system has been installed within the German-Indonesian Tsunami Early Warning System (GITEWS) project. Tsunamis are a global phenomenon and for global observations satellites are predestined. Within the GITEWS project a feasibility study on a future tsunami detection system from space has therefore been carried out. The Global Navigation Satellite System Reflectometry (GNSS-R) is an innovative way of using GNSS signals for remote sensing. It uses ocean reflected GNSS signals for sea surface altimetry. With a dedicated Low Earth Orbit (LEO) constellation of satellites equipped with GNSS-R receivers, densely spaced sea surface height measurements could be established to detect tsunamis. Some general considerations on the geometry between LEO and GNSS are made in this simulation study. It exemplary analyzes the detection performance of a GNSS-R constellation at 900 km altitude and 60° inclination angle when applied to the Sumatra tsunami as it occurred in 2004. GPS is assumed as signal source and the combination with GLONASS and Galileo signals is investigated. It can be demonstrated, that the combination of GPS and Galileo is advantageous for constellations with few satellites while the combination with GLONASS is preferable for constellations with many satellites. If all three GNSS are combined, the best detection performance can be expected for all scenarios considered. In this case an 18 satellite constellation will detect the Sumatra tsunami within 17 min with certainty, while it takes 53 min if only GPS is considered.     

A scheduling strategy to inter-satellite links assignment in GNSS

By introducing inter-satellite link (ISL), the dependence of the global navigation satellite system (GNSS) on ground infrastructure can be reduced and its performance enhanced via inter-satellite ranging and communication. Owing to platform restrictions, there are usually fewer onboard Ka-band ISL antennas than the number of visible satellites, which poses a problem when optimizing the inter-satellite links assignment of the GNSS. In this study, to optimize inter-satellite ranging and communication, a multi-objective optimization model is built and a scheduling strategy is proposed for the inter-satellite links assignment scheduling problem. The position dilution of precision (PDOP) of links and the transmission time-delay of telemetry data are set as the ranging performance and communication metrics, respectively. We regard the links assignment in each slot as a general graph-matching problem, and apply the Blossom algorithm to obtain the maximum matching. We then generate and optimize the satellite sequences for whole slots using non-dominated sorting genetic algorithm II (NSGA-II). The simulation scenes include 10,080 epochs of GNSS constellation, and the simulation results show that the performance of the proposed strategy is better than that of other methods published recently, and can provide various solutions to meet the different preferences of system managers.