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.     

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.     

Influence of time variable geopotential models on precise orbits of altimetry satellites,global and regional mean sea level trends

During the last decade a significant progress has been reached in the investigation of the gravity field of the Earth. Besides static, also time variable geopotential models have been recently created. In this paper we investigate the impact of the recent time variable geopotential models on altimetry satellite orbits and such altimetry products based on these orbits, as global and regional mean sea level trends. We show that the modeling of time variable gravity improves the orbit solutions, at least for the GRACE period where time variable gravity is sufficiently accurately observed by this mission. Our analysis includes six geopotential models jointly developed by GFZ German Research Centre for Geosciences and Space Geodesy Research Group (CNES/GRGS) Toulouse: the stationary model EIGEN-GL04S, a stationary version of EIGEN-6S (EIGEN-6S_stat), a corrected version of EIGEN-6S and three enhanced versions of EIGEN-6S called EIGEN-6S2, EIGEN-6S2A and EIGEN-6S2B. By “stationary” we mean “containing periodic parameters such as annual and semi-annual variations, but no secular (drift) terms”. We computed precise orbits for the radar altimetry satellites ERS-1, ERS-2, TOPEX/Poseidon, and Envisat over 20 years between 1991 and 2011. The orbit, single-mission and multi-mission altimetry crossover analyses show that the time variable models EIGEN-6S_corrected, EIGEN-6S2 and its two precursors EIGEN-6S2A/B perform notably better than the stationary models for the GRACE period from 2003 onwards. Thus, using EIGEN-6S2 and EIGEN-6S2A/B we have got 3.6% smaller root mean square fits of satellite laser ranging observations for Envisat, as when using EIGEN-GL04S. However, for the pre-GRACE period 1991–2003, the stationary geopotential models EIGEN-GL04S and EIGEN-6S_stat as well as EIGEN-6S2 having no drift terms for degree 3–50 at this time interval perform superior compared to EIGEN-6S_correct and EIGEN-6S2A/B which contain drifts for this period. We found, that the time variable geopotential models have a low (0.1–0.2 mm/yr) impact on our results for the global mean sea level trend. However, we found strong East/West differences up to 3 mm/yr in the regional mean sea level trends when using orbits of all four satellites based on time variable and stationary geopotential models. We show that these differences are related to the relative drifts of the centers-of-origin between the orbit solutions based on the time variable and stationary geopotential models. From the results of our detailed study, we conclude that the final version of the time variable gravity field model EIGEN-6S2 performs best for the four satellites tested. This model provides the most reliable and mission-consistent sea level estimates for the whole time period from 1992 to 2010. This model is of maximum spherical harmonic degree and order 260 and contains time series for drifts as well as annual and semiannual variations of the spherical harmonic coefficients for degree 2–50.     

Improvement of the SLR Borowiec station position in the global network ITRF91

This paper describes briefly two station coordinates solutions, the first one computed at Space Research Centre (SRC) using tracking data from Lageos satellite, and the second one computed at European Space Operations Centre (ESOC) using tracking data from Lageos 1, ERS-1 and TOPEX/Poseidon in a multi-arc solution. In particular the solution computed for the Borowiec station in ITRF91 system is described extensively. The Borowiec station position was estimated simultaneousely, considering or not the existence of range biases, with other geophysical parameters such as: daily polar motion parameters x_p, y_p, ocean tide coefficients, earth gravitational constant GM, etc.     

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.     

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.     

GRACE gravity field modeling with an investigation on correlation between nuisance parameters and gravity field coefficients

The GRACE (Gravity Recovery And Climate Experiment) monthly gravity models have been independently produced and published by several research institutions, such as Center for Space Research (CSR), GeoForschungsZentrum (GFZ), Jet Propulsion Laboratory (JPL), Centre National d’Etudes Spatiales (CNES) and Delft Institute of Earth Observation and Space Systems (DEOS). According to their processing standards, above institutions use the traditional variational approach except that the DEOS exploits the acceleration approach. The background force models employed are rather similar. The produced gravity field models generally agree with one another in the spatial pattern. However, there are some discrepancies in the gravity signal amplitude between solutions produced by different institutions. In particular, 10%–30% signal amplitude differences in some river basins can be observed. In this paper, we implemented a variant of the traditional variational approach and computed two sets of monthly gravity field solutions using the data from January 2005 to December 2006. The input data are K-band range-rates (KBRR) and kinematic orbits of GRACE satellites. The main difference in the production of our two types of models is how to deal with nuisance parameters. This type of parameters is necessary to absorb low-frequency errors in the data, which are mainly the aliasing and instrument errors. One way is to remove the nuisance parameters before estimating the geopotential coefficients, called NPARB approach in the paper. The other way is to estimate the nuisance parameters and geopotential coefficients simultaneously, called NPESS approach. These two types of solutions mainly differ in geopotential coefficients from degree 2 to 5. This can be explained by the fact that the nuisance parameters and the gravity field coefficients are highly correlated, particularly at low degrees. We compare these solutions with the official and published ones by means of spectral analysis. It is found that our solutions are, in general, consistent with others in the spatial pattern. The water storage variations of the Amazon, Chari and Ganges river basins have also been computed. The variations computed with the NPARB approach are closer to those produced by JPL and DEOS solutions, while the variations produced with the NPESS approach are in good agreement with those produced by the CSR and GFZ solutions. A simulation study is implemented with considering realistic noise and low-frequency error. The two approaches are used to recover the true model. The NPESS solution appears closer to the true one. Therefore we are inclined to estimate the nuisance parameters simultaneously with the geopential coefficients.     

The effect of geocenter motion on Jason-2 orbits and the mean sea level

We compute a series of Jason-2 GPS and SLR/DORIS-based orbits using ITRF2005 and the std0905 standards ( Lemoine et al., 2010). Our GPS and SLR/DORIS orbit data sets span a period of 2 years from cycle 3 (July 2008) to cycle 74 (July 2010). We extract the Jason-2 orbit frame translational parameters per cycle by the means of a Helmert transformation between a set of reference orbits and a set of test orbits. We compare the annual terms of these time-series to the annual terms of two different geocenter motion models where biases and trends have been removed. Subsequently, we include the annual terms of the modeled geocenter motion as a degree-1 loading displacement correction to the GPS and SLR/DORIS tracking network of the POD process. Although the annual geocenter motion correction would reflect a stationary signal in time, under ideal conditions, the whole geocenter motion is a non-stationary process that includes secular trends. Our results suggest that our GSFC Jason-2 GPS-based orbits are closely tied to the center of mass (CM) of the Earth consistent with our current force modeling, whereas GSFC’s SLR/DORIS-based orbits are tied to the origin of ITRF2005, which is the center of figure (CF) for sub-secular scales. We quantify the GPS and SLR/DORIS orbit centering and how this impacts the orbit radial error over the globe, which is assimilated into mean sea level (MSL) error, from the omission of the annual term of the geocenter correction. We find that for the SLR/DORIS std0905 orbits, currently used by the oceanographic community, only the negligence of the annual term of the geocenter motion correction results in a – 4.67 ± 3.40 mm error in the Z-component of the orbit frame which creates 1.06 ± 2.66 mm of systematic error in the MSL estimates, mainly due to the uneven distribution of the oceans between the North and South hemisphere.     

Impact of GPS antenna phase center variations on precise orbits of the GOCE satellite

The first European Space Agency Earth explorer core mission GOCE (Gravity field and steady-state Ocean Circulation Explorer) has been launched on March 17, 2009. The 12-channel dual-frequency Global Positioning System receiver delivers 1 Hz data and provides the basis for precise orbit determination (POD) on the few cm-level for such a very low orbiting satellite (254.9 km). As a member of the European GOCE Gravity Consortium, which is responsible for the GOCE High-level Processing Facility (HPF), the Astronomical Institute of the University of Bern (AIUB) provides the Precise Science Orbit (PSO) product for the GOCE satellite. The mission requirement for 1-dimensional POD accuracy is 2 cm. The use of in-flight determined antenna phase center variations (PCVs) is necessary to meet this requirement. The PCVs are determined from 154 days of data and the magnitude is up to 3-4 cm. The impact of the PCVs on the orbit determination is significant. The cross-track direction benefits most of the PCVs. The improvement is clearly seen in the orbit overlap analysis and in the validation with independent Satellite Laser Ranging (SLR) measurements. It is the first time that SLR could validate the cross-track component of a LEO orbit.     

Generating precise and homogeneous orbits for Jason-1 and Jason-2

Driven by the GMES (Global Monitoring for Environment and Security) and GGOS (Global Geodetic Observing System) initiatives the user community has a strong demand for high-quality altimetry products. In order to derive such high-quality altimetry products, precise orbits for the altimetry satellites are a necessity. With the launch of the TOPEX/Poseidon mission in 1992 a still on-going time series of high-accuracy altimetry measurements of ocean topography started, continued by the altimetry missions Jason-1 in 2001 and Jason-2/OSTM in 2008. This paper contributes to the on-going orbit reprocessing carried out by several groups and presents the efforts of the Navigation Support Office at ESA/ESOC using its NAPEOS software for the generation of precise and homogeneous orbits referring to the same reference frame for the altimetry satellites Jason-1 and Jason-2. Data of all three tracking instruments on-board the satellites (beside the altimeter), i.e. GPS, DORIS, and SLR measurements, were used in a combined data analysis. About 7 years of Jason-1 data and more than 1 year of Jason-2 data were processed. Our processing strategy is close to the GDR-C standards. However, we estimated slightly different scaling factors for the solar radiation pressure model of 0.96 and 0.98 for Jason-1 and Jason-2, respectively. We used 30 s sampled GPS data and introduced 30 s satellite clocks stemming from ESOC’s reprocessing of the combined GPS/GLONASS IGS solution. We present the orbit determination results, focusing on the benefits of adding GPS data to the solution. The fully combined solution was found to give the best orbit results. We reach a post-fit RMS of the GPS phase observation residuals of 6 mm for Jason-1 and 7 mm for Jason-2. The DORIS post-fit residuals clearly benefit from using GPS data in addition, as the DORIS data editing improves. The DORIS observation RMS for the fully combined solution is with 3.5 mm and 3.4 mm, respectively, 0.3 mm better than for the DORIS-SLR solution. Our orbit solution agrees well with external solutions from other analysis centers, as CNES, LCA, and JPL. The orbit differences between our fully combined orbits and the CNES GDR-C orbits are of about 0.8 cm for Jason-1 and at 0.9 cm for Jason-2 in the radial direction. In the cross-track component we observe a clear improvement when adding GPS data to the POD process. The 3D-RMS of the orbit differences reveals a good orbit consistency at 2.7 cm and 2.9 cm for Jason-1 and Jason-2. Our resulting orbit series for both Jason satellites refer to the ITRF2005 reference frame and are provided in sp3 file format on our ftp server.     

CubeSail: A low cost CubeSat based solar sail demonstration mission

CubeSail is a nano-solar sail mission based on the 3U CubeSat standard, which is currently being designed and built at the Surrey Space Centre, University of Surrey. CubeSail will have a total mass of around 3 kg and will deploy a 5 × 5 m sail in low Earth orbit. The primary aim of the mission is to demonstrate the concept of solar sailing and end-of-life de-orbiting using the sail membrane as a drag-sail. The spacecraft will have a compact 3-axis stabilised attitude control system, which uses three magnetic torquers aligned with the spacecraft principle axis as well as a novel two-dimensional translation stage separating the spacecraft bus from the sail. CubeSail’s deployment mechanism consists of four novel booms and four-quadrant sail membranes. The proposed booms are made from tape-spring blades and will deploy the sail membrane from a 2U CubeSat standard structure. This paper presents a systems level overview of the CubeSat mission, focusing on the mission orbit and de-orbiting, in addition to the deployment, attitude control and the satellite bus.     

On prediction of re-entry time of an upper stage from GTO

The evolution of objects in geostationary transfer orbit (GTO) is determined by a complex interplay of atmospheric drag and luni-solar gravity. These orbits are highly eccentric (eccentricity >0.7) and have large variations in velocity and perturbations during a revolution. The periodic changes in the perigee altitudes of these orbits are mainly due to the gravitational perturbations of the Sun and the Moon. The re-entry time of the objects in such orbits is sensitive to the initial conditions. The aim of this paper is to study the re-entry time of the cryogenic stage of the Indian geo-synchronous launch vehicle, GSLV-F04/CS, which has been decaying since 2 September 2007 from initial orbit with eccentricity equal to 0.706. Two parameters, initial eccentricity and ballistic coefficient, are chosen for optimal estimation. It is known that the errors are more in eccentricity for the observations based on two line elements (TLEs). These two parameters are computed with response surface method using a genetic algorithm for the selected eight different zones, based on rough linear variation of the mean apogee altitude during 200 days orbit evolution. The study shows that the GSLV-F04/CS will re-enter between 5 December 2010 and 7 January 2011. The methodology is also applied to study the re-entry of six decayed objects (cryogenic stages of GSLV and Molniya satellites). Good agreement is noticed between the actual and the predicted re-entry times. The absolute percentage error in re-entry prediction time for all the six objects is found to be less than 7%. The present methodology is being adopted at Vikram Sarabhai Space Centre (VSSC) to predict the re-entry time of GSLV-F04/CS.     

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.     

HUMEX, a study on the survivability and adaptation of humans to long-duration exploratory missions, part I: lunar missions.

The European Space Agency has recently initiated a study of the human responses, limits and needs with regard to the stress environments of interplanetary and planetary missions. Emphasis has been laid on human health and performance care as well as advanced life support developments including bioregenerative life support systems and environmental monitoring. The overall study goals were as follows: (i) to define reference scenarios for a European participation in human exploration and to estimate their influence on the life sciences and life support requirements; (ii) for selected mission scenarios, to critically assess the limiting factors for human health, wellbeing, and performance and to recommend relevant countermeasures; (iii) for selected mission scenarios, to critically assess the potential of advanced life support developments and to propose a European strategy including terrestrial applications; (iv) to critically assess the feasibility of existing facilities and technologies on ground and in space as testbeds in preparation for human exploratory missions and to develop a test plan for ground and space campaigns; (v) to develop a roadmap for a future European strategy towards human exploratory missions, including preparatory activities and terrestrial applications and benefits. This paper covers the part of the HUMEX study dealing with lunar missions. A lunar base at the south pole where long-time sunlight and potential water ice deposits could be assumed was selected as the Moon reference scenario. The impact on human health, performance and well being has been investigated from the view point of the effects of microgravity (during space travel), reduced gravity (on the Moon) and abrupt gravity changes (during launch and landing), of the effects of cosmic radiation including solar particle events, of psychological issues as well as general health care. Countermeasures as well as necessary research using ground-based test beds and/or the International Space Station have been defined. Likewise advanced life support systems with a high degree of autonomy and regenerative capacity and synergy effects were considered where bioregenerative life support systems and biodiagnostic systems become essential. Finally, a European strategy leading to a potential European participation in future human exploratory missions has been recommended.     

Survey-only optical strategies for cataloguing space debris objects in the future European space surveillance system

In the framework of a potential European Space Situational Awareness System (ESSAS), we propose some optical strategies such that try to minimize the requirement of tracking measurements for the orbit determination computation when the catalogue is under construction. We will analyse them in terms of coverage, timeliness and orbit determination accuracy by means of the AS4 simulator (developed by Deimos Space S.L.U.). Moreover, observation campaigns have been performed from La Sagra Observatory in order to check the applicability of those strategies. These strategies are used for defining different choices for the future European Optical Space Surveillance System in the framework of the ESA contract no. 22738/09/D/HK.     

Accuracy assessment of TDRSS-based TOPEX/Poseidon orbit determination

Orbit determination results are obtained for the Ocean Topography Experiment (TOPEX)/Poseidon spacecraft by the Goddard Space Flight Center (GSFC) Flight Dynamics Division (FDD) using a batch-least-squares estimator available in the Goddard Trajector Determination System (GTDS) to process Tracking and Data Relay Satellite (TDRS) System (TDRSS) measurements. The GTDS orbit solutions are compared with the definitive Precision Orbit Determination (POD) orbit solutions. The root-mean-square (RMS) solution difference in the radial component is 28 centimeters.     

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.     

Assessing attitude error of FORMOSAT-3/COSMIC satellites and its impact on orbit determination

An attitude determination and control system (ADCS) is critical to satellite attitude maneuvers and to the coordinate transformation from the inertial frame to the spacecraft frame. This paper shows specific sensors in the ADCS of the satellite mission FORMOSAT-3/COSMIC (F3/C) and the impact of the ADCS quality on orbit accuracy. The selection of main POD antenna depends on the beta angles of the different F3/C satellites (for FM2 and FM4) during the inflight phase. In particular, under the eclipse, alternative attitude sensors are activated to replace the Sun sensors, and such a sensor change leads to anomalous GPS phase residuals and a degraded orbit accuracy. Since the nominal attitude serves as a reference for ADCS, the 3-dimensional attitude-induced errors in reduced dynamic orbits over selected days in 2010 show 9.35, 10.78, 4.97, 5.48, 7.18, and 6.89 cm for FM1–FM6. Besides, the 3-dimensional velocity errors induced by the attitude effect are 0.10, 0.10, 0.07, 0.08, 0.09, and 0.10 for FM1–FM6. We analyze the quality of the observed attitude transformation matrix of F3/C and its impact on kinematic orbit determination. With 249 days of GPS in 2008, the analysis leads to the following averaged 3-dimensional attitude-induced orbit errors: 2.72, 2.62, 2.37, 1.90, 1.70, and 1.99 cm for satellites FM1–FM6. Critical suggestions of geodetic payloads for the follow-on mission of F3/C are presented based on the current result.     

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.