Geosynchronous satellites expanding a future GNSS satellite constellation: A precise orbit determination study

The integration of geosynchronous orbit (GSO) satellites in Global Navigation Satellite Systems (GNSS) is mostly discussed to enable a regional enhancement for tracking. But how do GSO satellites affect the orbit determination of the rest of the constellation? How accurately can these orbits be determined in a future GNSS tracking scenario with optical links? In this simulation study we analyze the benefit of GSO satellites as an expansion of a MEO (Medium Earth Orbit) satellite constellation – we selected the Galileo satellite constellation – for MEO Precise Orbit Determination (POD). We address not only the impact on POD of MEO satellites but also the possibility to precisely determine the GSO satellites – geostationary orbits (GEO) and inclined geosynchronous orbits (IGSO) – in such an expanded MEO constellation. In addition to GNSS microwave observations, we analyze the influence of different optical links between the participating entities: Optical two-way Inter-Satellite Links (OISL) and ground-space oriented Optical Two-Way Links (OTWL). These optical measurements together with the GNSS microwave observations give a remarkable benefit for the POD capability. In the case of GNSS and OTWL, we simulate the measurements with regard to a network of 16 ground stations. We pay great attention to the simulation of systematic effects of all measurement techniques. We discuss the influence on the systematic errors as well as the formal orbit uncertainties. A MEO constellation expanded with GSO satellites as well as the use of optical links together with GNSS observations not only improves the MEO satellite orbits but also the GSOs to a great extent.     

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

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

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.     

全球个人卫星移动通信系统

全球个人通信系统（ＰＣＳ）是一个发展中的新概念。１９９０年初开始提出利用静止轨道以外的其他轨道通信卫星的设想,最好是利用低轨道（ＬＥＯ）卫星。使用低轨道卫星就像把地面蜂窝网搬到天上,并能提供与之相仿的业务。文章评价了利用移动的和个人终端的通信卫星网;给出了未来低轨道卫星网的前景;对其中一些有关问题作了归纳。同时,描述了个人通信系统的功能性构成,并对一些网络的主要技术参数作表比较。     

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.     

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.     

Differential code bias estimation based on uncombined PPP with LEO onboard GPS observations

Differential Code Bias (DCB) is an essential correction that must be provided to the Global Navigation Satellite System (GNSS) users for precise position determination. With the continuous deployment of Low Earth Orbit (LEO) satellites, DCB estimation using observations from GNSS receivers onboard the LEO satellites is drawing increasing interests in order to meet the growing demands on high-quality DCB products from LEO-based applications, such as LEO-based GNSS signal augmentation and space weather research. Previous studies on LEO-based DCB estimation are usually using the geometry-free combination of GNSS observations, and it may suffer from significant leveling errors due to non-zero mean of multipath errors and short-term variations of receiver code and phase biases. In this study, we utilize the uncombined Precise Point Positioning (PPP) model for LEO DCB estimation. The models for uncombined PPP-based LEO DCB estimation are presented and GPS observations acquired from receivers onboard three identical Swarm satellites from February 1 to 28, 2019 are used for the validation. The results show that the average Root Mean Square errors (RMS) of the GPS satellite DCBs estimated with onboard data from each of the three Swarm satellites using the uncombined PPP model are less than 0.18 ns when compared to the GPS satellite DCBs obtained from IGS final daily Global Ionospheric Map (GIM) products. Meanwhile, the corresponding average RMS of GPS satellite DCBs estimated with the conventional geometry-free model are 0.290, 0.210, 0.281 ns, respectively, which are significantly larger than those obtained with the uncombined PPP model. It is also noted that the estimated GPS satellite DCBs by Swarm A and C satellites are highly correlated, likely attributed to their similar orbit type and space environment. On the other hand, the Swarm receiver DCBs estimated with uncombined PPP model, with Standard Deviation (STD) of 0.065, 0.037 and 0.071 ns, are more stable than those obtained from the official Swarm Level 2 products with corresponding STD values of 0.115, 0.101, and 0.109 ns, respectively. The above indicates that high-quality DCB products can be estimated based on uncombined PPP with LEO onboard observations.     

4D Lattice Flower Constellations

4D Lattice Flower Constellations is a new constellation design framework, based on the previous 2D and 3D Lattice theories of Flower Constellations, that focus on the generation of constellations whose satellites can have different semi-major axis and still present a constellation structure that is maintained during the dynamic of the system. This situation can arise when dealing with satellites with very different instruments, or when it is of interest to coordinate two different constellations. In that sense, 4D Lattice Flower Constellations constitutes the most general representation of the Flower Constellation formulation. In addition, the effects of the $J_2$ perturbation are taken into account in order to generate distributions that maintain their initial design configuration under this perturbation for longer periods of time with a low fuel budget. Finally, examples of application are presented, showing the possibilities in satellite constellation design of this new approach.     

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.     

VLBI observations to the APOD satellite

The APOD (Atmospheric density detection and Precise Orbit Determination) is the first LEO (Low Earth Orbit) satellite in orbit co-located with a dual-frequency GNSS (GPS/BD) receiver, an SLR reflector, and a VLBI X/S dual band beacon. From the overlap statistics between consecutive solution arcs and the independent validation by SLR measurements, the orbit position deviation was below 10?cm before the on-board GNSS receiver got partially operational. In this paper, the focus is on the VLBI observations to the LEO satellite from multiple geodetic VLBI radio telescopes, since this is the first implementation of a dedicated VLBI transmitter in low Earth orbit. The practical problems of tracking a fast moving spacecraft with current VLBI ground infrastructure were solved and strong interferometric fringes were obtained by cross-correlation of APOD carrier and DOR (Differential One-way Ranging) signals. The precision in X-band time delay derived from 0.1?s integration time of the correlator output is on the level of 0.1?ns. The APOD observations demonstrate encouraging prospects of co-location of multiple space geodetic techniques in space, as a first prototype.     

附加先验轨道约束的LEO几何法实时精密定轨验证分析

低轨卫星的实时精密定轨能够极大拓展其完成复杂科学任务的能力,例如实时环境监测、机动控制和卫星自主导航等.本文根据几何法实时精密定轨模型,提出了附加LEO先验轨道约束从而改善实时定轨的精度、收敛速度和稳定性的构想.分别采用广播星历、超快速星历预报部分和实时精密星历,设计了6种实时定轨方案,并利用Swarm-A,B,C星7天的观测数据进行方案验证与分析.结果表明,使用广播星历、IGU和IGC星历的方案精度递增,附加先验轨道约束能够进一步提升精度.使用IGC星历并附加标准差为1m的先验轨道约束后,在径向、切向和法向的定轨精度分别达到6.12cm,5.55cm和4.98cm.此外,附加先验轨道约束能够显著提升收敛速度,使用IGC星历平均收敛时间约为31min,附加标准差为1m的先验轨道约束后收敛仅需约4min.      

Medium Earth Orbit dynamical survey and its use in passive debris removal

The Medium Earth Orbit (MEO) region hosts satellites for navigation, communication, and geodetic/space environmental science, among which are the Global Navigation Satellites Systems (GNSS). Safe and efficient removal of debris from MEO is problematic due to the high cost for maneuvers needed to directly reach the Earth (reentry orbits) and the relatively crowded GNSS neighborhood (graveyard orbits). Recent studies have highlighted the complicated secular dynamics in the MEO region, but also the possibility of exploiting these dynamics, for designing removal strategies. In this paper, we present our numerical exploration of the long-term dynamics in MEO, performed with the purpose of unveiling the set of reentry and graveyard solutions that could be reached with maneuvers of reasonable $\mathrm{\Delta }V$ cost. We simulated the dynamics over 120–200?years for an extended grid of millions of fictitious MEO satellites that covered all inclinations from 0 to 90°, using non-averaged equations of motion and a suitable dynamical model that accounted for the principal geopotential terms, 3rd-body perturbations and solar radiation pressure (SRP). We found a sizeable set of usable solutions with reentry times that exceed $～40$ years, mainly around three specific inclination values: 46°, 56°, and 68°; a result compatible with our understanding of MEO secular dynamics. For $\mathrm{\Delta }V?300$ m/s (i.e., achieved if you start from a typical GNSS orbit and target a disposal orbit with $e<0.3$), reentry times from GNSS altitudes exceed $～70$ years, while low-cost ($\mathrm{\Delta }V?5–35$ m/s) graveyard orbits, stable for at lest 200?years, are found for eccentricities up to $e\approx 0.018$. This investigation was carried out in the framework of the EC-funded “ReDSHIFT” project.     

Track detection of high-velocity resident space objects in Low Earth Orbit

The world’s economy has become heavily dependent on the services provided by satellites. With the exponential increase in satellite launches, the population of defunct or inactive hardware in space has grown substantially. This is especially true in sensitive orbits such as the Low Earth Orbit (LEO) and Geostationary Earth Orbit (GEO) regimes. These objects, collectively known as orbital debris, can reach speeds of up to 28 000km h^?1 in LEO. At these orbital speeds, even the smallest of objects can pose a considerable threat to operational satellites or astronauts. This makes the monitoring, and detection, of these objects of the utmost importance. This work describes the latest detection strategy used in one of Europe’s largest Space Situational Awareness (SSA) installation; the BIstatic RAdar for LEo Survey (BIRALES) space debris radar. We present a novel bottom-up approach that makes use of single-linkage clustering to identify faint radar streaks in spectrogram data. Tests on synthetic data have shown that the detection strategy presented in this study obtains a higher detection rate when it is compared against existing methods. Unlike other approaches, this detection strategy, using the Multi-beam streak detection strategy (MSDS) algorithm, was still able to recall 90% of the track information at an Signal-to-Noise Ratio (SNR) of 2dB.