ERS-1 and ERS-2 operational and precise orbit determination

This paper presents the European Space Operations Centre's orbit determination and prediction systems for the ERS-1 mission. The routine operational orbit determination and prediction subsystem is discussed briefly, and statistics of the accuracy compared to the requirements are given. The precise orbit determination subsystem is then described, and the accuracy of its results are compared to those of the operational orbit system and to the D-PAF preliminary orbit solutions. Some geophysical results from the altimeter data, processed in these orbit determinations, are also presented. The ESOC/OAD ‘ERS-1 Orbit Report’ is introduced as a document providing this information on a monthly basis. Finally, this paper describes how the experience gained with the precise orbit determination will be exploited to further improve the accuracy of the routine system that will be used for ERS-2, and provides an estimate of this accuracy.     

ERS-1 radial position refinement by dual crossover analysis with TOPEX/Poseidon

Achieving orbital accuracies in the radial direction for ERS-1 commensurate with those for TOPEX/Poseidon is of utmost importance for the integration of the two altimeter data sets. This paper outlines a procedure whereby the radial orbit error for ERS-1 is recovered as a time series expansion in the form of a finite Fourier series with additional terms for atmospheric drag, solar radiation pressure, and initial state vector mismodelling. Using a least squares collocation method with constraints derived from the JGM2 gravity field co-variance matrix, the radial error is recovered using both dual crossovers and ERS-1 single satellite crossovers. Aggregate arcs are then used to derive the ERS-1 orbit error over the repeat period of 35 days. The results are presented in the improvement of fit in the dual crossover, ERS-1 crossover and altimetry data sets as well as the recovery of an altimeter bias for the two satellites.     

Dynamic and non-dynamic ERS-1 radial orbit improvement from ERS-1/TOPEX dual-satellite altimetry

Dual-satellite altimeter crossover differences between ERS-1 and TOPEX/Poseidon have been included as supplementary tracking data in ERS-1 orbit computations from SLR and single-satellite crossover differences. It was found that including the dual-satellite crossover differences slightly improves the ERS-1 radial orbit accuracy of about 12 cm for orbits computed with the JGM-2 gravity field and also leads to a better ‘centering’ of the ERS-1 orbit in the terrestrial reference frame defined for TOPEX/Poseidon. In addition to this dynamic orbit improvement technique, a non-dynamic technique has been investigated that removes the larger part of the ERS-1 radial orbit error from the dual-satellite crossover difference residuals. For ERS-1 orbits computed with the GEM-T2 gravity field, it was found that the non-dynamic technique could improve the radial orbit accuracy from 140 cm to the same level of accuracy as the ERS-1 JGM-2 orbits.     

New improved orbit solutions for the ERS-1 and ERS-2 satellites

Improved orbit solutions of the European Remote Sensing Satellites ERS-1 and ERS-2 have been computed in the ITRF2005 terrestrial reference frame using the recent models based mainly on IERS Conventions 2003. These solutions cover the periods 3 August 1991 to 8 July 1996 for ERS-1, and 3 May 1995 to 4 July 2003 for ERS-2. For each satellite, the final orbit solution is based on a combination of three separate orbit solutions independently computed at the Delft Institute of Earth Observation and Space Systems (DEOS) of the Delft University of Technology (The Netherlands), the Navigation Support Office of the European Space Operations Centre (ESOC, Germany) and the Helmholtz Centre Potsdam GFZ German Research Centre for Geosciences (Germany) using three different software packages for precise orbit determination, but using the same models in the same terrestrial reference frame within the European Space Agency (ESA) project ‘Reprocessing of Altimeter Products for ERS (REAPER)’. Validation using radar altimeter data indicates that the new combined orbits of ERS-1 and ERS-2 computed by us are significantly more accurate, approaching the 2–3 cm level in radial direction, than previously available orbit solutions.     

Gravity field model adjustment from ERS-1 and TOPEX altimetry and SLR tracking using an analytical orbit perturbation theory

The JGM-2 gravity field model has been adjusted using 70 days of ERS-1 and 19 10-day repeat cycles of TOPEX/Poseidon SLR and single satellite altimeter crossover differences. In addition, dual satellite altimeter crossover differences between ERS-1 and TOPEX for the selected 70-day period of ERS-1 and TOPEX repeat cycle 18 have been used in the JGM-2 model adjustment. In the computation of the normal equations, use was made of the analytical Lagrange linear perturbation theory. The single satellite altimeter crossover difference rms was brought down from 15.2 to 14.1 cm for ERS-1 and from 10.6 to 10.5 cm for TOPEX. The dual satellite altimeter crossover difference rms was reduced from 17.7 to 16.9 cm. Furthermore, the weighted rms of fit of SLR measurements was brought down from 16.0 to 14.4 cm for ERS-1, and from 5.4 to 5.0 cm for TOPEX/Poseidon.     

Precise orbit determination and sea-surface topography from ERS-1

Results of ERS - 1 orbit refinement procedures for long-arcs are summarised with particular reference to gravity field modelling and sea-surface topography. Long arcs from the three day and thirty - five day repeat cycles have been used to refine the earth's gravity field constrained by the JGM- 2 covariance matrix and simultaneously to solve for the quasi - stationary sea - surface topography. Orbital computations are presented for arcs not included within the refinement process with improvements identified through rms of fit to laser, altimeter and single satellite crossover data. Additional results utilising dual crossovers with TOPEX / POSEIDON in a full dynamic process are also presented. In particular, reference is made to the relative altimeter biases between the respective satellites and the characteristics of the orbital adjustment through inclusion of dual crossover data.     

High-fidelity geometry models for improving the consistency of CHAMP,GRACE, GOCE and Swarm thermospheric density data sets

During the last two decades, accelerometers on board of the CHAMP, GRACE, GOCE and Swarm satellites have provided high-resolution thermosphere density data to improve our knowledge on atmospheric dynamics and coupling processes in the thermosphere-ionosphere region. Most users of the data have focused on relative density variations. Scale differences between datasets and models have been largely neglected or removed using ad hoc scale factors. The origin of these scale differences arises from errors in the aerodynamic modelling, specifically in the modelling of the satellite outer surface geometry and of the gas-surface interactions. Therefore, the first step to remove the scale differences is to enhance the geometry modelling. This work forms the foundation for the future improvement of characterization of satellite aerodynamics and gas-surface interactions models at TU Delft, as well as for extending the use of sideways and angular accelerations in the aerodynamic analysis of accelerations and derivation of thermosphere datasets. Although work to improve geometry and aerodynamic force models by other authors has focused on CHAMP and GRACE, this paper includes the GOCE and Swarm satellites as well. In addition, it uses a density determination algorithm that is valid for arbitrary attitude orientations, enabling a validation making use of attitude manoeuvres. The results show an improvement in the consistency of density data between these four missions, and of data obtained before, during and after attitude manoeuvres of CHAMP and Swarm. The new models result in larger densities, compared to the previously used panel method. The largest average rescaling of density, by switching to the new geometry models is reached for Swarm at 32%, the smallest for GRACE at 5%. For CHAMP and GOCE, mean differences of 11% and 9% are obtained respectively. In this paper, an overview of the improvements and comparisons of data sets is provided together with an introduction to the next research phase on the gas-surface interactions.     

陆地卫星地面接收站微波系统扩充方案

讨论了为接收Landsat-6而扩充的陆地卫星地面接收站微波系统方案。该工作完成后,不仅能实现对Landsat-5、Landsat-6、SPOT 卫星的兼容,而且配上相应的解调器也能接收ERS-1、JERS-1及中国即将发射的资源卫星及其它工作在8025～8400MHz 的任何卫星的发射数据。     

Optical response of nanosatellite BLITS measured by the Graz 2 kHz SLR system

The nanosatellite BLITS (Ball Lens In The Space) is the first object designed as a passive, spherical retroreflector of the Luneburg type, dedicated for Satellite Laser Ranging (SLR). The optical response of BLITS has been measured by the Graz 2 kHz SLR station and compared with the response of the classical retroreflector arrays (RRA) of the Low Earth Orbiting satellites such as ERS-2 and Stella. This work demonstrates that the optical response of BLITS is flat and featureless, comparable with the signature of a point-source or a flat target, and suggests that this innovative design will deliver a higher normal point (NP) accuracy (2.55 mm) than any other SLR target currently in orbit. The high reflectivity of the glassy BLITS (about 60% of the return rate from the multi-reflector Stella) is found to be decreasing by about 30% per year, probably due to the solar irradiation. Detailed analysis of the reflective half-shell demonstrates that a high return rate of SLR measurements can be achieved regardless of the incident angle of the laser beam, thus making the spherical lens a perfect successor of the classical RRA panels mounted on active satellites such as CHAMP, GOCE and GRACE.     

不同后掠角三角翼的静态地面效应数值模拟

采用数值模拟的方法研究了不同后掠角三角翼的静态地面效应,通过对气动力和流场特性的分析发现,随着后掠角的减小,地面对迎风面下流动的阻滞作用增强,地效导致的迎风面气动力增量也随之增大。地效导致的背风面气动力增量同样随着后掠角的减小而增大,但在不同的后掠角范围内,地效诱导背风面气动力增量的机理不同:中大后掠角下,其主要通过增强前缘涡强度诱导更大的吸力,而小后掠角下,其主要通过促进前缘涡向内扩散增大吸力范围。      

ESA’s process for the identification and assessment of high-risk conjunction events

ESA’s Space Debris Office provides an operational service for the assessment of collision risks of ESA satellites. Currently, the ENVISAT and ERS-2 missions in low Earth orbits are covered by this service. If an upcoming high-risk conjunction event is predicted based on analysis of Two-Line Element (TLE) data from the US Space Surveillance Network, then independent tracking data of the potential high-risk conjunction object are acquired to improve the knowledge of its orbit. This improved knowledge and the associated small error covariances derived from the orbit determination process scale down the position error ellipsoid at the conjunction epoch. Hence, for the same miss-distance, in most cases an avoidance manoeuvre can be suppressed with an acceptable residual risk.     

超低轨卫星气动舵机辅助姿态控制方法设计

针对超低轨卫星姿态控制差异化需求,开展了基于气动舵机辅助的姿态控制策略研究。完成了超低轨道稀薄大气下卫星气动舵机布局设计与气动特性研究,理论气动力可达10^-1 N量级,气动力矩可达10^-1 N·m量级。在此基础上,完成了基于气动舵机辅助的姿态控制策略研究。通过仿真验证,在x轴采用动量轮控制、y轴和z轴采用气动舵机辅助控制情况下,可实现优于0.004°的三轴指向精度和优于0.000 7(°)/s的三轴姿态稳定度。所设计气动舵机辅助姿态控制策略对超低轨卫星技术应用与发展具有重要技术价值和工程意义。     

Evaluation of JGM 2 geopotential errors from Geosat, TOPEX/Poseidon and ERS-1 crossover altimetry

World-ocean distribution of the crossover altimetry data from Geosat, TOPEX/Poseidon (T/P) and the ERS 1 missions have provided strong independent evidence that NASA's/CSR's JGM 2 geopotential model (70 × 70 in spherical harmonics) yields accurate radial ephemerides for these satellites. In testing the sea height crossover differences found from altimetry and JGM 2 orbits for these satellites, we have used the sea height differences themselves (of ascending minus descending passes averaged at each location over many exact repeat cycles) and the Lumped Latitude Coefficients (LLC) derived from them. For Geosat we find the geopotential-induced LLC errors (exclusive of non-gravitational and initial state discrepancies) mostly below 6 cm, for TOPEX the corresponding errors are usually below 2 cm, and for ERS 1 (35-day cycle) they are generally below 5 cm. In addition, we have found that these observations agree well overall with predictions of accuracy derived from the JGM 2 variance-covariance matrix; the corresponding projected LLC errors for Geosat, T/P, and ERS 1 are usually between 1 and 4 cm, 1 – 2 cm, and 1 – 4 cm, respectively (they depend on the filtering of long-periodic perturbations and on the order of the LLC). This agreement is especially impressive for ERS 1 since no data of any kind from this mission was used in forming JGM 2.

The observed crossover differences for Geosat, T/P and ERS 1 are 8, 3, and 11 cm (rms), respectively. These observations also agree well with prediction of accuracy derived from the JGM 2 variance-covariance matrix; the corresponding projected crossover errors for Geosat and T/P are 8 cm and 2.3 cm, respectively. The precision of our mean difference observations is about 3 cm for Geosat (approx. 24,000 observations), 1.5 cm for T/P (approx. 6,000 observations) and 5 cm for ERS 1 (approx. 44,000 observations). Thus, these “global” independent data should provide a valuable new source for improving geopotential models. Our results show the need for further correction of the low order JGM 2 geopotential as well as certain resonant orders for all 3 satellites.     

翼型近波浪水面气动特性研究

采用求解Navier-Stokes方程的数值方法研究了翼型NACA4412近距离经过波浪水面时的气动特性。对数值方法的准确性进行了验证。计算了翼型经过波浪水面和固壁波浪地面2种边界条件下的气动力系数,并进行了对比。研究结果表明:翼型在经过波浪水面时气动力系数会发生周期性的变化,与固壁波浪的情况相比,气动力的变化曲线存在显著差异,波动幅度更大。通过对流场结构的分析,发现了翼型和波浪水面之间的作用机理。波浪水面的质点存在竖直方向上的运动,在小地面间隙时,水面质点向上运动会挤压翼型和水面之间的空气,从而造成翼型气动力大幅波动。同时解释了来流速度越大,气动力系数波动幅度减小的原因。      

Radar altimetry measurements over antarctic ice sheet: A focus on antenna polarization and change in backscatter problems

In this paper, we investigate the impact of the error due to the penetration of the altimetric wave within the snowpack. The phenomenon has two different impacts. The first one, due to temporal change in snow characteristics, affects the ice sheet volume trend as derived from altimetric series. The second one, because of both the anisotropy of the ice sheet surface properties and of the linear antenna polarization, introduces a difference in measurements at crossover points. These two phenomena are the cause of what are probably the most critical limitations to the interpretation of long-term altimetric series in term of mass balance and to the comparison between or data fusion of different missions. Moreover, they will lead to the largest error when comparing data from EnviSat with data from CryoSat, because of the different orbits, or with data from AltiKa, because of the different radar frequencies.     

无人直升机系留气动载荷CFD计算分析

系留气动载荷作为无人直升机系留装置的设计输入,通常以机身大风侧向角气动特性风洞试验数据为基础进行计算。采用CFD计算方法对某无人直升机算例样机的机身大风侧向角气动特性进行了计算,包括自由来流、停放在开阔地面和船艉甲板3个状态,以机身气动特性CFD计算结果为基础计算了其系留气动载荷。结果表明:无人直升机在开阔地面停放时的系留气动载荷与自由来流时基本一致。而受船体上层建筑的影响,停放在船艉甲板时的系留气动载荷与自由来流时有较大的差别,除部分风侧向角状态的偏航力矩之外,力和部分力矩的绝对值相对较小,部分风侧向角状态的力和力矩方向相反。研究结果可为选取无人直升机系留气动载荷计算方法和不同停放环境下的机身气动特性的CFD计算及风洞试验状态提供一定的参考。      

Decentralized differential drag based control of nanosatellites swarm spatial distribution using magnetorquers

Nanosatellites in the swarm initially move along arbitrary unbounded relative trajectories according to the launch initial conditions. Control algorithms developed in the paper are aimed to achieve the required spatial distribution of satellites in the along-track direction. The paper considers a swarm of 3U CubeSats in LEO, their form-factor is suitable for the aerodynamic control since the ratio of the satellite maximum to minimum cross-section areas is 3. Each satellite is provided with the information about the relative motion of neighboring satellites inside a specified communication area. The paper develops the corresponding decentralized control algorithms using the differential drag force. The required attitude control for each satellite is implemented by the active magnetic attitude control system. A set of decentralized control strategies is proposed taking into account the communicational constraints. The performance of these strategies is studied numerically. The swarm separation effect is demonstrated and investigated.     

Long-term photometric signature study of two GEO satellites

Geostationary earth orbit satellites have been extensively used for unique high-orbit stationary characteristics. Long-term precise investigation is an important issue in the observation of GEO satellites, since it can provide valuable information on the satellites’ operation state, discrimination and early warning analysis. Ground-based optical-electronic devices play a significant role in the observation. 4-month photometric signature variation of two satellites is presented based on the successive observations using the 1.56-meter telescope of Shanghai Astronomical Observatory (SHAO). It can be concluded that the long-term brightness change mainly results from the sun declination angle and regular orbit maneuver. Moreover, the solar panel offsets of the two satellites are analyzed and found to be approximately 4 degrees. Estimation of photometric accuracy reaches 0.15 mag for the application of CCD drift-scan optical-electronic technique.     

GNSS global real-time augmentation positioning: Real-time precise satellite clock estimation,prototype system construction and performance analysis

Lots of ambiguities in un-differenced (UD) model lead to lower calculation efficiency, which isn’t appropriate for the high-frequency real-time GNSS clock estimation, like 1 Hz. Mixed differenced model fusing UD pseudo-range and epoch-differenced (ED) phase observations has been introduced into real-time clock estimation. In this contribution, we extend the mixed differenced model for realizing multi-GNSS real-time clock high-frequency updating and a rigorous comparison and analysis on same conditions are performed to achieve the best real-time clock estimation performance taking the efficiency, accuracy, consistency and reliability into consideration. Based on the multi-GNSS real-time data streams provided by multi-GNSS Experiment (MGEX) and Wuhan University, GPS + BeiDou + Galileo global real-time augmentation positioning prototype system is designed and constructed, including real-time precise orbit determination, real-time precise clock estimation, real-time Precise Point Positioning (RT-PPP) and real-time Standard Point Positioning (RT-SPP). The statistical analysis of the 6 h-predicted real-time orbits shows that the root mean square (RMS) in radial direction is about 1–5 cm for GPS, Beidou MEO and Galileo satellites and about 10 cm for Beidou GEO and IGSO satellites. Using the mixed differenced estimation model, the prototype system can realize high-efficient real-time satellite absolute clock estimation with no constant clock-bias and can be used for high-frequency augmentation message updating (such as 1 Hz). The real-time augmentation message signal-in-space ranging error (SISRE), a comprehensive accuracy of orbit and clock and effecting the users’ actual positioning performance, is introduced to evaluate and analyze the performance of GPS + BeiDou + Galileo global real-time augmentation positioning system. The statistical analysis of real-time augmentation message SISRE is about 4–7 cm for GPS, whlile 10 cm for Beidou IGSO/MEO, Galileo and about 30 cm for BeiDou GEO satellites. The real-time positioning results prove that the GPS + BeiDou + Galileo RT-PPP comparing to GPS-only can effectively accelerate convergence time by about 60%, improve the positioning accuracy by about 30% and obtain averaged RMS 4 cm in horizontal and 6 cm in vertical; additionally RT-SPP accuracy in the prototype system can realize positioning accuracy with about averaged RMS 1 m in horizontal and 1.5–2 m in vertical, which are improved by 60% and 70% to SPP based on broadcast ephemeris, respectively.     

飞行载荷外部气动力的二次规划等效映射方法

传统的气动力映射方法中，守恒型方法能保证总积分载荷的守恒，但计算复杂；非守恒型方法实现简单，能保证气动压力的分布，但不能保证总积分载荷的守恒。利用等式约束保证积分载荷相等，并基于非守恒型加入线性基的径向基函数（RPIM）构造最小二乘保持原气动压力的分布特征，将气动力等效传递问题转化为一个带等式约束的二次规划问题。通过在非守恒方法的基础上构造出新的守恒型方法，兼顾了2种方法的优点。对某一大展弦比机翼上表面的三维气动力等效映射至二维平板面元气动力的结果表明，所提方法在牺牲了一定气动压力分布精度的基础上，严格保证了总积分载荷的相等，同时也继承了非守恒插值方法的高效性。所提方法可用于各类复杂外部气动力数据的快速等效映射，以及其他学科类似的数据映射问题。