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.     

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.     

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.     

Gravity model development for precise orbit computations for satellite altimetry

Satellite altimeters have been developed which have the capability of measuring the surface topography of the ocean over entire basins for time periods of several years with a precision of a few cm. The combination of these data with accurate satellite position information, subsurface measurements and models of the density field can be used to determine the general circulation of the ocean and its variability. The Topex altimeter oceanographic experiment proposed for the early 1990's requires a radial orbit accuracy approaching 10 cm. Considerable progress has been made at the Goddard Space Flight Center in developing precision orbit determination to support satellite altimetry analysis. This progress is attributed to the development of more accurate models of the earth gravity field, earth and ocean tidal effects, a reference coordinate system as well as improved drag and solar radiation pressure computations. This computational system has been used for the Seasat data acquired during the 1978 time period and has produced an rms radial accuracy of about 50 cm. A first step towards the development of a new Interim gravity model is presented. Global orbit determination accuracies have been improved using these new PGS-T1 and T2 fields.     

Selenodetic experiments of SELENE: Relay subsatellite,differential VLBI,and laser altimeter

Since 1960s, the gravitational potential of the Moon has been extensively studied from Doppler tracking data between a ground station and spacecraft orbiting in front of the Moon (e. g., Lorell and Sjogren, 1968; Bills and Ferrari, 1980; Konopliv et al., 1993; Lemoine et al., 1997). Because direct radio communication is interrupted while spacecraft is orbiting behind the Moon, however, the coverage of tracking data has been limited mostly to the nearside of the Moon so far. In order to compensate for such lack of tracking data, we propose satellite-to-satellite Doppler measurement by using a relay subsatellite in Japanese mission to the Moon in 2003. A complete coverage of Doppler tracking from an orbiter at sufficiently low altitude will significantly improve lunar gravity model and will contribute for future geophysical study of interior and tectonics on the Moon. Further, we propose differential VLBI experiment between the subsatellite and a propulsion module landed on the surface of the Moon. The differential VLBI is about 10 times more accurate than conventional Doppler measurement for long-wavelength gravity field. Besides, differential VLBI is sensitive to the displacement perpendicular to the line of sight. Thus the VLBI experiment provides precise estimates of the lunar gravity potential at low degree. The last proposal for selenodetic experiments is a laser altimeter. Global topography model has been already developed from the analysis of Clementine LIDAR data (Zuber et al., 1994), but it is suggested that the model includes appreciable anisotropy between NS and E-W directions due to highly eccentric orbit of Clementine spacecraft (Bills and Lemoine, 1995). The laser altimeter experiment from an orbiter in nearly circular orbit will provide a new reference for the isotropic lunar topography model.     

Analysis of radial orbit errors of ERS-1, and the development of super-tailored gravity models

The altimetry mission of the future ESA remote sensing satellite ERS-1 requires very accurate orbit solutions, of which in particular the radial position component should have an accuracy of approximately 10 cm. This paper presents some investigations into the possibility of reducing the radial position error due to the earth's gravity field, which is by far the largest contributing error source.With a detailed harmonic analysis of the ERS-1 orbit a number of gravity field model terms are identified which produce the major radial orbit perturbations. These dominant terms are adjusted in a least-squares orbit determination and parameter estimation procedure using actual SEASAT laser tracking observations and altimeter height measurements. The initial gravity model is the NASA GEM-L2 model derived from satellite tracking data only, with an emphasis on LAGEOS data. The resulting super-tailored model yields a significantly improved radial accuracy relative to GEM-L2, but fails to reach the accuracy of the SEASAT-tailored model PGS-S4.Finally, the SEASAT altimeter residuals and the residuals of the cross-over differences are analyzed in the frequency domain by applying a special filtering technique which separates the major radial orbit error and geoid error contributions.     

Co-registration of Chang’E-1 stereo images and laser altimeter data with crossover adjustment and image sensor model refinement

Topographic information is of fundamental importance for various scientific investigations in lunar and planetary exploration. To provide high-precision, seamless mapping capability, it is critical to co-register image and altimetry data, the two major data sources for topographic modeling. This paper presents a method for co-registration of Chang’E-1 (CE-1) stereo images and laser altimeter (LAM) data with crossover adjustment and refinement of the image sensor model. The crossover adjustment is tested in a larger area (0–60°N, 50–0°W); the image refinement and co-registration with LAM data are tested in an area (46.2–50.0°N, 31.8–28.8°W) within the larger area. Experimental results demonstrate that this co-registration reduces the mean differences of inconsistency from more than 200 m to just 3.21 m in the Z direction of object space. In image space, the mean errors of homologous points both in the column and row directions are reduced to below 0.1 pixel. This indicates that the proposed crossover adjustment of LAM data and refinement of the CE-1 stereo image model can effectively improve co-registration of the two data sets.     

The challenges in long-term altimetry calibration for addressing the problem of global sea level change

Long-term change of the global sea level resulting from climate change has become an issue of great societal interest. The advent of the technology of satellite altimetry has modernized the study of sea level on both global and regional scales. In combination with in situ observations of the ocean density and space observations of Earth’s gravity variations, satellite altimetry has become an essential component of a global observing system for monitoring and understanding sea level change. The challenge of making sea level measurements with sufficient accuracy to discern long-term trends and allow the patterns of natural variability to be distinguished from those linked to anthropogenic forcing rests largely on the long-term efforts of altimeter calibration and validation. The issues of long-term calibration for the various components of the altimeter measurement system are reviewed in the paper. The topics include radar altimetry, the effects of tropospheric water vapor, orbit determination, gravity field, tide gauges, and the terrestrial reference frame. The necessity for maintaining a complete calibration effort and the challenges of sustaining it into the future are discussed.     

Towards development of a consistent orbit series for TOPEX,Jason-1, and Jason-2

The TOPEX/Poseidon, Jason-1 and Jason-2 set of altimeter data now provide a time series of synoptic observations of the ocean that span nearly 17 years from the launch of TOPEX in 1992. The analysis of the altimeter data including the use of altimetry to monitor the global change in mean sea level requires a stable, accurate, and consistent orbit reference over the entire time span. In this paper, we describe the recomputation of a time series of orbits that rely on a consistent set of reference frames and geophysical models. The recomputed orbits adhere to the IERS 2003 standards for ocean and earth tides, use updates to the ITRF2005 reference frame for both the SLR and DORIS stations, apply GRACE-derived models for modeling of the static and time-variable gravity, implement the University College London (UCL) radiation pressure model for Jason-1, use improved troposphere modeling for the DORIS data, and apply the GOT4.7 ocean tide model for both dynamical ocean tide modeling and for ocean loading. The new TOPEX orbits have a mean SLR fit of 1.79 cm compared to 2.21 cm for the MGDR-B orbits. These new TOPEX orbits agree radially with independent SLR/crossover orbits at 0.70 cm RMS, and the orbit accuracy is estimated at 1.5–2.0 cm RMS over the entire TOPEX time series. The recomputed Jason-1 orbits agree radially with the Jason-1 GDR-C orbits at 1.08 cm RMS. The GSFC SLR/DORIS dynamic and reduced-dynamic orbits for Jason-2 agree radially with independent orbits from the CNES and JPL at 0.70–1.06 cm RMS. Applying these new orbits, and using the latest altimeter corrections for TOPEX, Jason-1, and Jason-2 from September 1992 to May 2009, we find a global rate in mean sea level of 3.0 ± 0.4 mm/yr.     

Chang’E-1 precision orbit determination and lunar gravity field solution

In this paper we present results assessing the role of Very Long Baseline Interferometry (VLBI) tracking data through precision orbit determination (POD) during the check-out phase for Chang’E-1, and the lunar gravity field solution CEGM-01 based on the orbital tracking data acquired during the nominal phase of the mission. The POD of Chang’E-1 is performed using S-band two-way Range and Range Rate (R&RR) data, together with VLBI delay and delay rate observations. The role of the VLBI data in the POD of Chang’E-1 is analyzed, and the resulting orbital accuracies are estimated for different solution strategies. The final orbital accuracies proved that the VLBI tracking data can improve the Chang’E-1 POD significantly. Consequently, CEGM-01 based on six-month tracking data during Chang’E-1 nominal mission phase is presented, and the accuracy of the model is assessed by means of the gravity field power spectrum, admittance and coherence between gravity and topography, lunar surface gravity anomaly and POD for both Chang’E-1 and Lunar Prospector (LP). Our analysis indicates that CEGM-01 has significant improvements over a prior model (i.e. GLGM-2), and shows the potential of Chang’E-1 tracking data in high resolution lunar gravity field model solution by combining with SELENE and LP tracking data.     

DORIS-based point mascons for the long term stability of precise orbit solutions

In recent years non-tidal Time Varying Gravity (TVG) has emerged as the most important contributor in the error budget of Precision Orbit Determination (POD) solutions for altimeter satellites’ orbits. The Gravity Recovery And Climate Experiment (GRACE) mission has provided POD analysts with static and time-varying gravity models that are very accurate over the 2002–2012 time interval, but whose linear rates cannot be safely extrapolated before and after the GRACE lifespan. One such model based on a combination of data from GRACE and Lageos from 2002–2010, is used in the dynamic POD solutions developed for the Geophysical Data Records (GDRs) of the Jason series of altimeter missions and the equivalent products from lower altitude missions such as Envisat, Cryosat-2, and HY-2A. In order to accommodate long-term time-variable gravity variations not included in the background geopotential model, we assess the feasibility of using DORIS data to observe local mass variations using point mascons. In particular, we show that the point-mascon approach can stabilize the geographically correlated orbit errors which are of fundamental interest for the analysis of regional Mean Sea Level trends based on altimeter data, and can therefore provide an interim solution in the event of GRACE data loss. The time series of point-mass solutions for Greenland and Antarctica show good agreement with independent series derived from GRACE data, indicating a mass loss at rate of 210 Gt/year and 110 Gt/year respectively.     

Determination of the mean dynamic ocean topography model through combining multi-source gravity data and DTU15 MSS around China's coast

Mean dynamic ocean topography (or MDT) is closely related to ocean circulation and global climate change. It has important scientific significance and application value for the development and utilization of marine resources in China's coastal areas. Based on the terrain gravity, marine gravity, and SRTM 3?s data, an algorithm to reduce the problem of gravity data gaps between land and sea is proposed. A consistent land-sea gravity model is established based on point-mass fusion method. Then geoid model, which accuracy was estimated to be 8.5?cm through the verification of 348 GNSS/level data from the coastal provinces, of China's coastal areas was calculated through remove-restore technique. Connecting the above geoid model with DTU15 MSS model to establish a MDT model in China's coastal areas using the direct method in space domain. The effect of gravity field model, dominant factors of sea surface topography, and low pass filter are analyzed. Taking Bohai Sea and Yellow Sea as an example, and comparing MDT with the two international models CNES_2013_MDT and DTU15_MDT. The results show that the MDT has the potential to construct a vertical datum of the ocean and carry out related scientific research and application.     

Altimetry,bathymetry and geoid variations at the Gavdos permanent Cal/Val facility

The aim of this work has been to examine the relationship of steep bathymetry in the coastal areas around the permanent Cal/Val facility of Gavdos, and their influence on the produced calibration values for the Jason-2 satellite altimeter. The paper describes how changes in seafloor topography (from 200 to 3500 m depth over a distance of 10 km) are reflected on the determined altimeter parameters using different reference surfaces for satellite calibration. Finally, it describes the relation between these parameter trends and the region’s local characteristics.     

Necessary conditions for the optimality of singular arcs of spacecraft trajectories subject to multiple gravitational bodies

This document analyzes the optimality of intermediate thrust arcs (singular arcs) of spacecraft trajectories subject to multiple gravitational bodies. A series of necessary conditions for optimality are formally derived, including the generalized Legendre–Clebsch condition. As the order of singular optimality turns out to be two, an explicit formula for the singular optimal control is also presented. These analytical outcomes are validated by showing that they are identical to Lawden’s classical result if the equations of motion are reduced for a central gravity field. Practical utility is demonstrated by applying these analytical derivations to a candidate optimal trajectory near the Moon subject to solar and Earth perturbation. While the candidate optimal trajectory turns out to be bang-singular-bang, the intermediate thrust arc satisfies all the necessary conditions for optimality.     

基于UKF的雷达高度计自主定轨

探讨了利用推广卡尔曼滤波估计非线性系统状态时存在的问题,从而介绍了目前广泛使用的分步逼近的卡尔曼滤波(UKF,Unscented Kalman Filter).为了提高导航的可靠性和准确性,在星敏感器导航系统中引入雷达高度计作为一个新的测量设备,提出了一种基于星上雷达测高仪及星敏感器联合进行卫星自主定轨的算法.建立了比较复杂的地球海平面模型,并考虑了其中风生重力波的影响. 利用雷达测高仪的测量结果和地球形状模型,计算地心矢量在卫星本体中坐标系的方向.利用UKF滤波定轨算法,明显提高了自主定轨的精度.数值仿真结果表明,UKF定轨精度要远优于推广卡尔曼滤波.      

An improved altimeter-derived gravity anomaly from shipborne gravity based on the mean sea surface height constraint factors method

Coastal marine gravity modeling faces challenges due to the degradation of the quality and poor coverage of altimeter data in coastal regions. The effective fusion of shipborne gravity data and altimeter-derived marine gravity data can make shipborne gravity data more useful for the accurate estimation of altimeter-derived coastal marine gravity. A mean sea surface height constraint factor (MSSHCF) method based on the ordinary kriging method and the remove-restore technique is proposed to fuse altimeter-derived gravity model with shipborne gravity data. In this method, all data are standardized during the interpolation process to reduce the error and mean sea surface as a vertical variable is added to the semi-variance function in ordinary kriging to obtain the residual shipborne gravity as corrected data source. The coastal marine gravity models V2.1 and V3.1 which fused altimeter-derived gravity data with shipborne gravity data and V1.1 without shipborne gravity data at a spatial resolution of 1′×1′ can be obtained. Validation experiments show that the accuracy of the gravity model V3.1 obtained by the MSSHCF method more closely agrees with the validated gravity model DTU17 and SS V31 than the model V2.1 obtained by the ordinary kriging interpolation method and the V1.1 model. Our results were validated against shipborne gravity data; the accuracy of model V3.1 was 4.95 % higher than the model V1.1 in South China Sea area A and 2.48 % higher in South China Sea area B. Meanwhile, the accuracy of model V3.1 was 2.07 % higher than model V2.1 in South China Sea area A and 2.42 % higher in South China Sea area B. The effects of distance from the coast and sea depth on the marine gravity model were also evaluated. The results show that the gravity model V3.1 has higher accuracy with the change in ocean distance and depth than the V2.1 and V1.1 gravity models. Thus, our study shows that the MSSHCF method effectively refines coastal altimeter-derived gravity using shipborne gravity data.     

Calibration of Envisat radar altimeter over Lake Issykkul

This study presents the results of calibration/validation (C/V) of Envisat satellite radar altimeter over Lake Issykkul located in Kyrgyzstan, which was chosen as a dedicated radar altimetry C/V site in 2004. The objectives are to estimate the absolute altimeter bias of Envisat and its orbit based on cross-over analysis with TOPEX/Poseidon (T/P), Jason-1 and Jason-2 over the ocean. We have used a new method of GPS data processing in a kinematic mode, developed at the Groupe de Recherche de Geodesie Spatiale (GRGS), which allows us to calculate the position of the GPS antenna without needing a GPS reference station. The C/V is conducted using various equipments: a local GPS network, a moving GPS antenna along the satellites tracks over Lake Issykkul, In Situ level gauges and weather stations. The absolute bias obtained for Envisat from field campaigns conducted in 2009 and 2010 is between 62.1 and 63.4 ± 3.7 cm, using the Ice-1 retracking algorithm, and between 46.9 and 51.2 cm with the ocean retracking algorithm. These results differ by about 10 cm from previous studies, principally due to improvement of the C/V procedure. Apart from the new algorithm for GPS data processing and the orbit error reduction, more attention has been paid to the GPS antenna height calculation, and we have reduced the errors induced by seiche over Lake Issykkul. This has been assured using cruise data along the Envisat satellite track at the exact date of the pass of the satellite for the two campaigns. The calculation of the Envisat radar altimeter bias with respect to the GPS levelling is essential to allow the continuity of multi-mission data on the same orbit, with the expected launch of SARAL/Altika mission in 2012. Implications for hydrology in particular, will be to produce long term homogeneous and reliable time series of lake levels worldwide.     

On the information contents and regularisation of lunar gravity field solutions

Till the present day the recovery of the lunar gravity field from satellite tracking data depends in a crucial way on the level and method of regularisation. With Earth-based tracking only, the spatial data coverage is limited to only slightly more than 50% and the inverse problem remains severely ill-posed. The development of global gravity models suitable for precise orbit modeling as well as geophysical studies therefore requires a significant level of regularisation, limiting the solution power over the far-side where no gravity information is available. Unconstrained solutions, within the framework of global harmonic base functions, are only possible for very low degrees (< 10). Any significant change to this situation is only to be expected when global satellite-to-satellite tracking data of high quality becomes available early in the next decade. Yet, a rigorous analysis of the impact of the chosen method and level of regularisation is lacking. Most gravity models employ a Kaula-type signal smoothness constraint of 15 × 10⁻⁵ /l², which allows a good overall data fit as well as a smooth field over the far-side. Furthermore, a geographical type of constraint has been suggested, where surface accelerations have been introduced in areas of no data coverage. Modern numerical methods, on the other hand, offer direct tools and search mechanisms for the optimal level of regularisation. This paper presents a study of Tikhonov-type regularisation of lunar gravity solutions, with emphasis on the so-called L-curve and quasi-optimality methods for regularisation parameter estimation. Furthermore, new quality measures of lunar gravity solutions are presented, which account for the bias introduced by the regularisation.     

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.