Precision orbit determination performance for CryoSat-2

In this paper we discuss our efforts to perform precision orbit determination (POD) of CryoSat-2 which depends on Doppler and satellite laser ranging tracking data. A dynamic orbit model is set-up and the residuals between the model and the tracking data is evaluated. The average r.m.s. of the 10?s averaged Doppler tracking pass residuals is approximately 0.39?mm/s; and the average of the laser tracking pass residuals becomes 1.42?cm. There are a number of other tests to verify the quality of the orbit solution, we compare our computed orbits against three independent external trajectories provided by the CNES. The CNES products are part of the CryoSat-2 products distributed by ESA. The radial differences of our solution relative to the CNES precision orbits shows an average r.m.s. of 1.25?cm between Jun-2010 and Apr-2017. The SIRAL altimeter crossover difference statistics demonstrate that the quality of our orbit solution is comparable to that of the POE solution computed by the CNES. In this paper we will discuss three important changes in our POD activities that have brought the orbit performance to this level. The improvements concern the way we implement temporal gravity accelerations observed by GRACE; the implementation of ITRF2014 coordinates and velocities for the DORIS beacons and the SLR tracking sites. We also discuss an adjustment of the SLR retroreflector position within the satellite reference frame. An unexpected result is that we find a systematic difference between the median of the 10 s Doppler tracking residuals which displays a statistically significant pattern in the South Atlantic Anomaly (SSA) area where the median of the velocity residuals varies in the range of ?0.15 to +0.15?mm/s.     

Standing window of oblique detonation with pathological behaviour

Standing of an Oblique Detonation Wave (ODW) on a wedge within combustor is the prerequisite of thrust generation for ODW engine which is regarded as a novel and conceptual propulsion device with hypersonic flight Mach number. Usually a standing window of ODW is defined as the wedge angle ranged from the ODW detached angle from wedge (upper limit) to the angle that a Chapman-Jouguet (CJ) detonation occurs (lower limit). For pathological detonation cases, however, the CJ detonation cannot be achieved, and thus the lower limit of the standing window of ODW should be revisited. In present study, two types of reactions in hypersonic incoming flow that include the behavior of pathological detonation, that is, the single-step irreversible reaction with mole variation and the two-step irreversible reactions with exothermic process followed by endothermic process, have been used for studying standing behavior of ODW. The steady detonation polar analysis of ODW is carried out for both reaction systems. The results reveal that the reaction with more mole decrement and the reactions with stronger endothermic process show the pathological detonation feature and therefore modify the lower limit of standing window of ODW. Three equivalent parameters are proposed to quantitatively measure the standing window range of ODW from points of view of thermodynamics, Mach number of incoming flow and heat effect of reactions. It is found that the standing window of ODW is determined by the specific heat ratio, the overdrive degree of detonation and the endothermic level of the hypersonic incoming flow, regardless of whether the detonation is pathological or not.     

Inter-agency comparison of TanDEM-X baseline solutions

TanDEM-X (TerraSAR-X add-on for Digital Elevation Measurement) is the first Synthetic Aperture Radar (SAR) mission using close formation flying for bistatic SAR interferometry. The primary goal of the mission is to generate a global digital elevation model (DEM) with 2 m height precision and 10 m ground resolution from the configurable SAR interferometer with space baselines of a few hundred meters. As a key mission requirement for the interferometric SAR processing, the relative position, or baseline vector, of the two satellites must be determined with an accuracy of 1 mm (1D RMS) from GPS measurements collected by the onboard receivers. The operational baseline products for the TanDEM-X mission are routinely generated by the German Research Center for Geosciences (GFZ) and the German Space Operations Center (DLR/GSOC) using different software packages (EPOS/BSW, GHOST) and analysis strategies. For a further independent performance assessment, TanDEM-X baseline solutions are generated at the Astronomical Institute of the University of Bern (AIUB) on a best effort basis using the Bernese Software (BSW).     

Analysis of BDS-3 distributed autonomous navigation based on BeiDou system simulation platform

Satellite autonomous navigation is an important function of the BeiDou-3 navigation System (BDS-3). Satellite autonomous navigation means that the navigation satellite uses long-term forecast ephemeris and Inter-Satellite Link (ISL) measurements to determinate its own spatial position and time reference without the support of the ground Operation and Control System (OCS) for a long time to ensure that the navigation system can normally maintain the time and space reference. This paper aims to analyze the feasibility of distributed autonomous navigation algorithms. For the first time, a ground parallel autonomous navigation test system (GPANTS) is built. The performance of distributed autonomous navigation is then analyzed using the two-way ISL ranging of BDS-3 satellites. First, the BDS simulation platform and the GPANTS are introduced. Then, the basic principles of distributed satellite autonomous orbit determination and time synchronization based on ISL measurements are summarized. Preliminary evaluation of the performance of the BDS-3 constellation autonomous navigation service under ideal conditions through simulation data. Then the performance of autonomous navigation for 22 BeiDou-3 satellites using ISL measurements is evaluated. The results show that when satellites operate autonomously for 50 days without the support of any ground station, the User Range Error (URE) of autonomous orbit determination is better than 3 m, and the time synchronization accuracy is better than 4 ns.     

摆动式扫描红外地球敏感器在轨探头级备份工作研究

针对摆动式扫描红外地球敏感器两个探测器工作时输出弦宽信息的特点以及现役通信卫星的一般配置,提出了使用两个探测器信息进行姿态确定及利用双份地球敏感器进行探头级备份工作的方法,该方法在工程实现上简单易行,实现精度较高,且提高了在轨卫星地球敏感器的使用可靠性及冗余度.     

北斗星座仿真的单天线姿态确定技术研究

载体姿态测量的实际需要和多天线姿态测量的局限性决定了单天线姿态测量系统拥有广泛的发展前景。分析了单天线测姿的应用价值和基本原理,从我国的具体国情出发,提出了基于北斗卫星星座仿真的单天线姿态确定技术研究的方法。利用STK平台对目前组网运行的北斗卫星星座和理想的北斗满卫星星座进行了仿真,利用Matlab平台进行了实验,仿真实验结果和误差分析表明,北斗卫星星座仿真的单天线姿态确定技术在适用条件和测姿准确性等方面效果理想,在一定精度要求下满足应用需要,对建立实用可靠的单天线测姿具有积极的推动作用。     

Characteristics of GOCE orbits based on Satellite Laser Ranging

The Gravity field and steady-state Ocean Circulation Explorer (GOCE) was the first European Space Agency’s (ESA) Earth Explorer core mission. Through its extremely low, about 260?km above the Earth, circular, sun-synchronous orbit, the satellite gained high spatial resolution and accuracy gravity gradient, and ocean circulation data. Global Positioning System (GPS) receivers, mounted on the spacecraft, allowed the determination of reduced-dynamic and kinematic GOCE orbits, whereas Laser Retroreflector Array (LRA) dedicated to Satellite Laser Ranging (SLR) allowed an independent validation of GPS-derived orbits. In this paper, residuals between different GPS-based orbit types and SLR observations are used to investigate the sensitivity and the influence of solar, geomagnetic, and ionospheric activities on the quality of kinematic and reduced-dynamic GOCE orbits. We also analyze the quality of data provided by individual SLR sites, by detecting time biases using ascending and descending sun-synchronous GOCE orbit passes, and the residual analysis of the measurement characteristics, i.e., the dependency of SLR residuals as a function of nadir and horizontal angles. Results show a substantial vulnerability of kinematic orbit solutions to the solar F10.7 index and the ionospheric activity measured by the variations of the Total Electron Content (TEC) values. The sensitivity of kinematic orbits to the three-hour-range KP index is rather minor. The reduced-dynamic orbits are almost insensitive to indices describing ionospheric, solar, and geomagnetic activities. The investigation of individual SLR sites shows that some of them are affected by time bias errors, whereas other demonstrate systematics, such as a dependency between observation residuals and the satellite nadir angle or the horizontal azimuth angle from the SLR station to the direction of the satellite.     

Analysis of laser ranges and angular measurements data fusion for space debris orbit determination

In the framework of space debris, the orbit determination process is a fundamental step, both, for researchers and for satellite operators. The accurate knowledge of the orbit of space debris objects is needed to allow space debris characterization studies and to avoid unnecessary collision avoidance maneuvers.The accuracy of the results of an orbit determination process depends on several factors as the number, the accuracy, the kind of processed measurements, their distribution along the orbit, and the object-observer relative geometry. When the observation coverage of the target orbit is not homogeneous, the accuracy of the orbit determination can be improved processing different kind of observables. Recent studies showed that the satellite laser ranging technique can be successfully applied to space debris.In this paper, we will investigate the benefits of using laser ranges and angular measurements for the orbit determination process. We will analyze the influence of the number of used observations, of the covered arc of orbit, of each observable, and of the observation geometry on the estimated parameters. Finally, using data acquired on short observation arcs, we analyze the achievable accuracies for the orbital regimes with the highest space debris density, and to the consequences of the data fusion on catalog maintenance operations. The results shown are obtained using only real data (both angular and laser measurements) provided by sensors of the Swiss Optical Ground Station and Geodynamics Observatory Zimmerwald owned by the Astronomical Institute of the University of Bern (AIUB) and for some studies also using ranges provided from other stations of the International Laser Ranging Service (ILRS).     

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.     

Proposal of application of same beam VLBI measurements in precision orbit determination of lunar orbiter and return capsule and lunar gravity field simulation in Chang’E-3 mission

High accuracy differenced phase delay can be obtained by observing multiple point frequencies of two spacecraft using the same beam Very Long Baseline Interferometry (VLBI) technology. Its contribution in lunar spacecraft precision orbit determination has been performed during the Japanese lunar exploration mission SELENE. In consideration that there will be an orbiter and a return capsule flying around the moon during the Chinese lunar exploration future mission Chang’E-3, the contributions of the same beam VLBI in spacecraft precision orbit determination and lunar gravity field solution have been investigated. Our results show that the accuracy of precision orbit determination can be improved more than one order of magnitude after including the same beam VLBI measurements. There are significant improvements in accuracy of low and medium degree coefficients of lunar gravity field model obtained from combination of two way range and Doppler and the same beam VLBI measurements than the one that only uses two way range and Doppler data, and the accuracy of precision orbit determination can reach meter level.