This paper evaluates orbit accuracy and systematic error for altimeter satellite precise orbit determination on TOPEX, Jason-1, Jason-2 and Jason-3 by comparing the use of four SLR/DORIS station complements from the International Terrestrial Reference System (ITRS) 2014 realizations with those based on ITRF2008. The new Terrestrial Reference Frame 2014 (TRF2014) station complements include ITRS realizations from the Institut National de l'Information Géographique et Forestière (IGN) ITRF2014, the Jet Propulsion Laboratory (JPL) JTRF2014, the Deutsche Geodätisches Forschungsinstitut (DGFI) DTRF2014, and the DORIS extension to ITRF2014 for Precise Orbit Determination, DPOD2014. The largest source of error stems from ITRF2008 station position extrapolation past the 2009 solution end time. The TRF2014 SLR/DORIS complement impact on the ITRF2008 orbit is only 1–2 mm RMS radial difference between 1992–2009, and increases after 2009, up to 5 mm RMS radial difference in 2016. Residual analysis shows that station position extrapolation error past the solution span becomes evident even after two years, and will contribute to about 3–4 mm radial orbit error after seven years. Crossover data show the DTRF2014 orbits are the most accurate for the TOPEX and Jason-2 test periods, and the JTRF2014 orbits for the Jason-1 period. However for the 2016 Jason-3 test period only the DPOD2014-based orbits show a strong and statistically significant margin of improvement. The positive results with DTRF2014 suggest the new approach to correct station positions or normal equations for non-tidal loading before combination is beneficial. We did not find any compelling POD advantage in using non-linear over linear station velocity models in our SLR & DORIS orbit tests on the Jason satellites. The JTRF2014 proof-of-concept ITRS realization demonstrates the need for improved SLR+DORIS orbit centering when compared to the Ries (2013) CM annual model. Orbit centering error is seen as an annual radial signal of 0.4 mm amplitude with the CM model. The unmodeled CM signals show roughly a 1.8 mm peak-to-peak annual variation in the orbit radial component. We find the TRF network stability pertinent to POD can be defined only by examination of the orbit-specific tracking network time series. Drift stability between the ITRF2008 and the other TRF2014-based orbits is very high, the relative mean radial drift error over water is no larger than 0.04 mm/year over 1993–2015. Analyses also show TRF induced orbit error meets current altimeter rate accuracy goals for global and regional sea level estimation.

We have investigated the quality of precise orbits for the SARAL altimeter satellite using Satellite Laser Ranging (SLR) and Doppler Orbitography and Radiopositioning Integrated by Satellite (DORIS) data from March 14, 2013 to August 10, 2014. We have identified a 4.31±0.14cm error in the Z (cross-track) direction that defines the center-of-mass of the SARAL satellite in the spacecraft coordinate system, and we have tuned the SLR and DORIS tracking point offsets. After these changes, we reduce the average RMS of the SLR residuals for seven-day arcs from 1.85 to 1.38cm. We tuned the non-conservative force model for SARAL, reducing the amplitude of the daily adjusted empirical accelerations by eight percent. We find that the best dynamic orbits show altimeter crossover residuals of 5.524cm over cycles 7–15. Our analysis offers a unique illustration that high-elevation SLR residuals will not necessarily provide an accurate estimate of radial error at the 1-cm level, and that other supporting orbit tests are necessary for a better estimate. Through the application of improved models for handling time-variable gravity, the use of reduced-dynamic orbits, and through an arc-by-arc estimation of the C22 and S22 coefficients, we find from analysis of independent SLR residuals and other tests that we achieve 1.1–1.2cm radial orbit accuracies for SARAL. The limiting errors stem from the inadequacy of the DPOD2008 and SLRF2008 station complements, and inadequacies in radiation force modeling, especially with respect to spacecraft self-shadowing and modeling of thermal variations due to eclipses.

Forward modeling of non-tidal atmospheric loading displacements at geodetic tracking stations have not routinely been included in Doppler Orbitography and Radiopositionning Integrated by Satellite (DORIS) or Satellite Laser Ranging (SLR) station analyses for either POD applications or reference frame determination. The displacements which are computed from 6-hourly models such as the ECMWF and can amount to 3-10 mm in the east, north and up components depending on the tracking station locations. We evaluate the application of atmospheric loading in a number ways using the NASA GSFC GEODYN software: First we assess the impact on SLR & DORIS-determined orbits such as Jason-2, where we evaluate the impact on the tracking data RMS of fit and how the total orbits are changed with the application of this correction. Preliminary results show an RMS radial change of 0.5 mm for Jason-2 over 54 cycles and a total change in the Z-centering of the orbit of 3 mm peak-to-peak over one year. We also evaluate the effects on other DORIS-satellites such as Cryosat-2, Envisat and the SPOT satellites. In the second step, we produce two SINEX time series based on data from available DORIS satellites and assess the differences in WRMS, scale and Helmert translation parameters. Troposphere refraction is obviously an important correction for radiometric data types such as DORIS. We evaluate recent improvements in DORIS processing at GSFC including the application of the Vienna Mapping Function (VMF1) grids with a-priori hydrostatic (VZHDs) and wet (VZWDs) zenith delays. We reduce the gridded VZHD at the stations height using pressure and temperature derived from GPT (strategy 1) and Saastamoinen. We discuss the validation of the VMF1 implementation and its application to the Jason-2 POD processing, compared to corrections using the Niell mapping function and the GMF. Using one year of data, we also assess the impact of the new troposphere corrections on the DORIS-only solutions, most importantly on the scale of the weekly solutions.

The JASON-2 satellite, launched in June 2008, is the latest follow-on to the successful TOPEX/Poseidon (T/P) and JASON-I altimetry missions. JASON-2 is equipped with a TRSR Blackjack GPS dual-frequency receiver, a laser retroreflector array, and a DORIS receiver for precise orbit determination (POD). The most recent time series of orbits computed at NASA GSFC, based on SLR/DORIS data have been completed using both ITRF2005 and ITRF2008. These orbits have been shown to agree radially at 1 cm RMS for dynamic vs SLRlDORIS reduced-dynamic orbits and in comparison with orbits produced by other analysis centers (Lemoine et al., 2010; Zelensky et al., 2010; Cerri et al., 2010). We have recently upgraded the GEODYN software to implement model improvements for GPS processing. We describe the implementation of IGS standards to the Jason2 GEODYN GPS processing, and other dynamical and measurement model improvements. Our GPS-only JASON-2 orbit accuracy is assessed using a number of tests including analysis of independent SLR and altimeter crossover residuals, orbit overlap differences, and direct comparison to orbits generated at GSFC using SLR and DORIS tracking, and to orbits generated externally at other centers. Tests based on SLR and the altimeter crossover residuals provide the best performance indicator for independent validation of the NASAlGSFC GPS-only reduced dynamic orbits. For the ITRF2005 and ITRF2008 implementation of our GPS-only obits we are using the IGS05 and IGS08 standards. Reduced dynamic versus dynamic orbit differences are used to characterize the remaining force model error and TRF instability. We evaluate the GPS vs SLR & DORIS orbits produced using the GEODYN software and assess in particular their consistency radially and the stability of the altimeter satellite reference frame in the Z direction for both ITRF2005 and ITRF2008 as a proxy to assess the consistency of the reference frame for altimeter satellite POD.

The long-term stability and the precision of the satellite orbit is a critical component of the Jason-1 and Jason-2 (OSTM) Missions, providing the reference frame for ocean mapping using altimeter data. DORIS tracking in combination with SLR has provided orbits, which are both highly accurate and consistent across missions using the latest and most accurate POD models. These models include GRACE-derived static and time varying gravity fields and a refined Terrestrial Reference Frame based on SLR and DORIS data yielding a uniform station complement. Additional improvements have been achieved based on advances in modeling the satellite surface forces and the tropospheric path delay for DORIS measurements. This paper presents these model improvements for Jason-1 and Jason-2, including a description of DORIS sensitivity to error in tropospheric path delay. We show that the detailed University College London (UCL) radiation pressure model for Jason-1, which includes self-shadowing and thermal re-radiation, is superior to the use of a macromodel for radiation pressure surface force modeling. Improvements in SLR residuals are seen over all Beta-prime angles for both Jason-1 and Jason-2 using the UCL model, with the greatest improvement found over regimes of low Beta-prime where orbit Earth shadowing is maximum. The overall radial orbit improvement for Jason-1 using the UCL model is 3 mm RMS, as corroborated by the improvement in the independent altimeter crossover data. Special attention is paid to Jason-2 POD to assess improvements gained with the latest advances in DORIS receiver technology. Tests using SLR and altimeter crossover residuals suggest the Jason-2 reduced-dynamic DORIS-only, SLR/DORIS, and GPS orbits have all achieved 1-cm radial accuracy. Tests using independent SLR data acquired at high elevation show an average fit value of 1.02 cm for the DORIS-only and 0.94 cm for the GPS reduced-dynamic orbits. Orbit differences suggest that the largest remaining errors in the Jason-2 dynamic orbit solutions are due to radiation pressure mis-modeling and variations in the geopotential not captured in the GRACE-derived annual terms.

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 17years 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.79cm compared to 2.21cm for the MGDR-B orbits. These new TOPEX orbits agree radially with independent SLR/crossover orbits at 0.70cm RMS, and the orbit accuracy is estimated at 1.5–2.0cm RMS over the entire TOPEX time series. The recomputed Jason-1 orbits agree radially with the Jason-1 GDR-C orbits at 1.08cm 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.06cm 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.4mm/yr.

Lunar laser ranging (LLR) measurements are crucial for advanced exploration of the laws of fundamental gravitational physics and geophysics as well as for future human and robotic missions to the Moon. The corner-cube reflectors (CCR) currently on the Moon require no power and still work perfectly since their installation during the project Apollo era. Current LLR technology allows its to measure distances to the Moon with a precision approaching I mm. As NASA pursues the vision of taking humans back to the Moon, new, more precise laser ranging applications will be demanded, including continuous tracking from more sites on Earth, placing, new CCR arrays on the Moon, and possibly installing other devices such as transponders, etc. for multiple scientific and technical purposes. Since this effort involves humans in space, then in all situations the accuracy, fidelity, and robustness of the measurements, their adequate interpretation, and any products based on them, are of utmost importance. Successful achievement of this goal strongly demands further significant improvement of the theoretical model of the orbital and rotational dynamics of the Earth-Moon system. This model should inevitably be based on the theory of general relativity, fully incorporate the relevant geophysical processes, lunar librations, tides, and Should rely upon the most recent standards and recommendations of the IAU for data analysis. This paper discusses methods and problems in developing such a mathematical model. The model will take into account all the classical and relativistic effects in the orbital and rotational motion of the Moon and Earth at the sub-centimeter level. The model is supposed to be implemented as a part of the computer code underlying NASA Goddard's orbital analysis and geophysical parameter estimation package GEODYN and the ephemeris package PMOE 2003 of the Purple Mountain Observatory. The new model will allow us to navigate a spacecraft precisely to a location on the Moon. It will also greatly improve our understanding of the structure of the lunar interior and the nature of the physical interaction at the core-mantle interface layer. The new theory and upcoming millimeter LLR will give us the means to perform one of the most precise fundamental tests of general relativity in the solar system. (C) 2008 COSPAR. Published by Elsevier Ltd. All rights reserved.

The Terrestrial Reference Frame is the fundamental means by which we relate observations in space and time. For example, in order to generate a homogeneous and consistent time series of geo-referenced altimeter measurements over the span of the Topex/Poseidon and Jason-1 missions, we must examine carefully the role of improvements in measurement modelling, force modelling, and improved reference frame realizations. In this paper, we quantify the effects of improvements in force modelling, for example the use of new GRACE-derived gravity models, the effect of time-variable gravity derived from GRACE on altimeter satellite orbits. In addition, we examine the effects of modelling geocenter in altimeteric satellite POD, and look at how the application of atmospheric loading might affect the time-series of precise orbits for Topex/Poseidon and Jason-1.

We have analyzed GRACE Level 1-B data in 2003 and assessed a new approach for extracting time variable gravity that isolates the gravity signal in both time and space. The Level-1B satellite-to-satellite range rate (KBRR) data and accelerometry are processed in daily arcs using the precise orbit products produced by the GRACE team from GPS to calibrate both the accelerometer and KBRR data. We then adjusted select components of the intersatellite baseline vector for each data segment isolated to the region of interest. Herein, we solved for mass anomalies in 45 deg x 45 deg blocks over the Amazon and the nearby Atlantic Ocean and estimate mass flux in units of cm of water over each block. We show with this approach that we can recover mass anomalies on a submonthly basis with 10 to 15 day temporal resolution. We discuss the important issues related to this solution, including the size of the mascon blocks, the weight given to the temporal and spatial constraint used to stabalize the solutions, as well as the optimal correlation in time and distance. We compare the the mascon results with solutions obtained from the more standard approach using spherical harmonics and with independent hydrology models and lake data. This technique demonstrates that sub-monthly medium wavelength mass flux phenomena are well sensed by the hyper-precise line of sight velocity data produced from GRACE.