Impact of ITRS 2014 realizations on altimeter satellite precise orbit determination

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.     

Quality Assessment of the IDS Contribution to ITRF2008

Doppler Orbitography Radiopositionning Integrated by Satellite (DORIS) is one of the four fundamental techniques contributing to the ITRF. The optimal coverage over the globe of the DORIS observing sites and sites co-located with GPS, allow a strong embedding of DORIS within the ITRF network. DORIS contributes to the access to ITRF through precise orbit determination of altimetric satellites with onboard DORIS receivers. The DORIS contribution to the ITRF2008 is enhanced by the fact that the solutions of seven analysis centers were included in the submitted combined time series of weekly station positions and daily polar motion. We evaluate the quality of the DORIS combined solution in terms of its agreement with the other techniques (VLBI, SLR, GPS) contributing to the ITRF2008 combination. We show in particular that the precisions of the current IDS products range between 1.5 to 2.6 mm for station positions (at the epochs of minimum variances); better than 1 mm/yr in velocities and between 170 and 260 micro-arc-seconds for polar motion, a significant improvement by a factor of three to five, compared to past data used in the ITRF2005 combination. This improvement is certainly due to improved analysis strategies employed by the seven IDS analysis centers that contributed to the combined weekly submitted solutions of station positions and polar motion. A spectral analysis of DORIS station height time series indicates that annual and semi-annual signals are dominant. However, TOPEX draconitic period of about 118 days is still detected in about 20% of the station position power spectra. DORIS height annual signals correlate well with GPS annual signal estimated at some co-located stations, which show that DORIS technique is able to detect loading signals.     

Analysis of the DORIS contributions to ITRF2008

In its function as an ITRS Combination Centre, DGFI is in charge with the computation of an ITRF2008 solution. The computation methodology of DGFI is based on the combination of datum-free normal equations (weekly or session data sets, respectively) of station positions and Earth orientation parameters (EOP) from the geodetic space techniques DORIS, GPS, SLR and VLBI. In this paper we focus on the DORIS part within the ITRF2008 computations. We present results obtained from the analysis of the DORIS time series for station positions, network translation and scale parameters, as well as for the terrestrial pole coordinates. The submissions to ITRF2008 benefit from improved analysis strategies of the seven contributing IDS analysis centres and from a combination of the weekly solutions of station positions and polar motion. The results show an improvement by a factor of two compared to past DORIS data submitted to ITRF2005, which has been evaluated by investigating the repeatabilities of position time series. The DORIS position time series were analysed w.r.t. discontinuities and other non-linear effects such as seasonal variations. About 40 discontinuities have been identified which have been compared with the results of an earlier study. Within the inter-technique combination we focus on the DORIS contribution to the integration of the different space geodetic observations and on a comparison of the geodetic local ties with the space geodetic solutions. Results are given for the 41 co-location sites between DORIS and GPS.     

Determination of station positions and velocities from laser ranging observations to Ajisai,Starlette and Stella satellites

The positions and velocities of the four Satellite Laser Ranging (SLR) stations: Yarragadee (7090), Greenbelt (7105), Graz (7839) and Herstmonceux (7840) from 5-year (2001–2005) SLR data of low orbiting satellites (LEO): Ajisai, Starlette and Stella were determined. The orbits of these satellites were computed from the data provided by 20 SLR stations. All orbital computations were performed by means of NASA Goddard’s GEODYN-II program. The geocentric coordinates were transformed to the topocentric North–South, East–West and Vertical components in reference to ITRF2005. The influence of the number of normal points per orbital arc and the empirical acceleration coefficients on the quality of station coordinates was studied. To get standard deviation of the coordinates determination lower than 1 cm, the number of the normal points per site had to be greater than 50. The computed positions and velocities were compared to those derived from LAGEOS-1/LAGEOS-2 data. Three parameters were used for this comparison: station coordinates stability, differences from ITRF2005 positions and velocities. The stability of coordinates of LEO satellites is significantly worse (17.8 mm) than those of LAGEOS (7.6 mm), the better results are for Ajisai (15.4 mm) than for Starlette/Stella (20.4 mm). The difference in positions between the computed values and ITRF2005 were little bit worse for Starlette/Stella (6.6 mm) than for LAGEOS (4.6 mm), the results for Ajisai were five times worse (29.7 mm) probably due to center of mass correction of this satellite. The station velocities with some exceptions were on the same level (≈1 mm/year) for all satellites. The results presented in this work show that results from Starlette/Stella are better than those from Ajisai for station coordinates determination. We can applied the data from LEO satellites, especially Starlette and Stella for determination of the SLR station coordinates but with two times lower accuracy than when using LAGEOS data.     

DPOD2005: An extension of ITRF2005 for Precise Orbit Determination

For Precise Orbit Determination of altimetry missions, we have computed a data set of DORIS station coordinates defined for specific time intervals called DPOD2005. This terrestrial reference set is an extension of ITRF2005. However, it includes all new DORIS stations and is more reliable, as we disregard stations with large velocity formal errors as they could contaminate POD computations in the near future. About 1/4 of the station coordinates need to be defined as they do not appear in the original ITRF2005 realization. These results were verified with available DORIS and GPS results, as the integrity of DPOD2005 is almost as critical as its accuracy. Besides station coordinates and velocities, we also provide additional information such as periods for which DORIS data should be disregarded for specific DORIS stations, and epochs of coordinate and velocity discontinuities (related to either geophysical events, equipment problem or human intervention). The DPOD model was tested for orbit determination for TOPEX/Poseidon (T/P), Jason-1 and Jason-2. Test results show DPOD2005 offers improvement over the original ITRF2005, improvement that rapidly and significantly increases after 2005. Improvement is also significant for the early T/P cycles indicating improved station velocities in the DPOD2005 model and a more complete station set. Following 2005 the radial accuracy and centering of the ITRF2005-original orbits rapidly degrades due to station loss.     

DPOD2014: A new DORIS extension of ITRF2014 for precise orbit determination

To support precise orbit determination of the altimetry missions, the International DORIS Service (IDS) regularly estimates the DPOD (DORIS terrestrial reference frame for Precise Orbit Determination) solution which includes mean positions and velocities of all the DORIS stations. This solution is aligned to the current realization of the International Terrestrial Reference Frame (ITRF) and so, can be seen as a DORIS extension of the ITRF. In 2016, moving to the IDS Combination Center, the DPOD construction scheme changed. The new DPOD solution is produced from a DORIS cumulative position and velocity solution. We present the new methodology used to compute DPOD2014 and its validation procedure. In order to present geophysical applications and interpretations of these results, we show two examples: (1) the Gorkha earthquake (M7.8 – April 2015) generates a 3-D mis-positioning of nearly 55?mm of the EVEB DORIS station at the Everest base camp 90?km from the epicenter. (2) Applying the results the DPOD2014 realization, we show that the most recent vertical velocity of Thule, Greenland is similar to that observed between 2006 and 2010, indicating further ongoing ice mass loss in the Thule region of northwest Greenland.     

A VLBI baseline post-adjustment approach for station velocity estimation in Eurasian continent

Baseline lengths and their time-derivatives among 58 geodetic VLBI stations were fitted by using 4439 observing sessions from the International VLBI Service for Geodesy and Astrometry (IVS). First, the velocities of eight stations in Eurasian continent were set as unknown quantities. Then, two standard global solutions from 3523 IVS sessions and 1110 sessions from database code XA, respectively, were applied prior to all-station coordinates and the non-estimated station velocities. Finally, from the relations among the coordinates, velocities, baseline length and its time-derivative, two types of baseline post-adjustment (BPA) were used to estimate the velocities of the eight stations. We discuss the data processing details, including the effect of different prior values for the stations and the optimal solution.     

The International DORIS Service (IDS): Toward maturity

DORIS is one of the four space-geodetic techniques participating in the Global Geodetic Observing System (GGOS), particularly to maintain and disseminate the Terrestrial Reference Frame as determined by International Earth rotation and Reference frame Service (IERS). A few years ago, under the umbrella of the International Association of Geodesy, a DORIS International Service (IDS) was created in order to foster international cooperation and to provide new scientific products. This paper addresses the organizational aspects of the IDS and presents some recent DORIS scientific results. It is for the first time that, in preparation of the ITRF2008, seven Analysis Centers (AC’s) contributed to derive long-term time series of DORIS stations positions. These solutions were then combined into a homogeneous time series IDS-2 for which a precision of less than 10 mm was obtained. Orbit comparisons between the various AC’s showed an excellent agreement in the radial component, both for the SPOT satellites (e.g. 0.5–2.1 cm RMS for SPOT-2) and Envisat (0.9–2.1 cm RMS), using different software packages, models, corrections and analysis strategies. There is now a wide international participation within IDS that should lead to future improvements in DORIS analysis strategies and DORIS-derived geodetic products.     

Observation level combination of SLR and VLBI with c5++: A case study for TIGO

A multi-technique space geodetic analysis software named c5++ has been developed and allows one to combine data on the observation level. With SLR and VLBI modules being ready and tested, this software has been used to compute coordinate time series of the geodetic fundamental station TIGO, located near Concepción, Chile. It can be shown that the combination of space geodetic data on the observation level leads to a significant improvement of station position repeatability, which is an important measure for the stability of a station in the terrestrial reference frame. Moreover, it could be demonstrated that the geophysical signal of the post-seismic tectonic plate movement is usually more complete than detected by any of the two single-technique solutions. In addition, it has been confirmed that so-called nuisance parameters, which are relying on data from a single technique, are not biased when combing observations from different space geodetic techniques.     

Polar motion and earth rotation series from LAGEOS

As an Associate Analysis Centre for the MERIT campaign ESOC has processed the complete data set of LAGEOS laser ranging normal points for the 14 months period of the campaign (September 1983 to October 1984), using independently developed software. Geocentric coordinates have been obtained for a network of 20 laser stations, and series for the pole, UT1R and length of day have been derived with 5 day resolution. Discrepancies in the node and inclination estimates between successive monthly arcs are attributed to mismodelled tidal effects. Detailed comparisons have been made with solutions obtained at other processing centres based on VLBI and SLR techniques.     

Refining DORIS atmospheric drag estimation in preparation of ITRF2008

In preparation of ITRF2008, all geodetic technique services (VLBI, SLR, GPS and DORIS) are generating new solutions based on combination of individual analysis centers solutions. These data reprocessing are based on a selection of models, parameterization and estimation strategy unique to each analysis center and to each technique. While a good agreement can be found for models between groups, thanks to the existence of the IERS conventions, a great diversity still exist for parameter estimation, allowing possible future improvements in this direction. The goal of this study is to focus on the atmospheric drag estimation used to generate the new DORIS/IGN ignwd08 time series prepared for ITRF2008. We develop here a method to inter-compare different processing strategies. In a first step, by analyzing single-satellite solutions for a few weeks of data but for a large number of possible analysis strategies, we demonstrate that estimating drag coefficient more frequently (typically every 1–2 h instead of previously every 4–8 h) for the lowest DORIS satellites (SPOTs and Envisat) provides better geodetic results for station coordinates and polar motion. This new processing strategy also solved earlier problem found when processing DORIS data during intense geomagnetic events, such as geomagnetic storms. Differences between drag estimation strategies can mostly be found during these few specific periods of extreme geomagnetic activity (few days per year). In such a case, when drag coefficient is only estimated every 6 h or less often for single-satellite solution, a significant degradation in station coordinate accuracy can be observed (120 mm vs. 20 mm) and significant biases arose in polar motion estimation (5 mas vs. 0.3 mas). In a second step, we reprocessed a full year of DORIS data (2003) in a standard multi-satellite mode. We were able to provide statistics on a more reliable data set and to strengthen these conclusions. Our proposed DORIS analysis is easy to implement in all software packages and is now already used by several analysis centers of the International DORIS Service (IDS) when submitting reprocessed solutions for ITRF2008.     

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.     

Creation of the new industry-standard space test of laser retroreflectors for the GNSS and LAGEOS

We built a new experimental apparatus (the “Satellite/lunar laser ranging Characterization Facility”, SCF) and created a new test procedure (the SCF-Test) to characterize and model the detailed thermal behavior and the optical performance of cube corner laser retroreflectors in space for industrial and scientific applications. The primary goal of these innovative tools is to provide critical design and diagnostic capabilities for Satellites Laser Ranging (SLR) to Galileo and other GNSS (Global Navigation Satellite System) constellations. The capability will allow us to optimize the design of GNSS laser retroreflector payloads to maximize ranging efficiency, to improve signal-to-noise conditions in daylight and to provide pre-launch validation of retroreflector performance under laboratory-simulated space conditions. Implementation of new retroreflector designs being studied will help to improve GNSS orbits, which will then increase the accuracy, stability, and distribution of the International Terrestrial Reference Frame (ITRF), to provide better definition of the geocenter (origin) and the scale (length unit).     

A balloon-borne very long baseline interferometry experiment in the stratosphere: Systems design and developments

The balloon-borne very long baseline interferometry (VLBI) experiment is a technical feasibility study for performing radio interferometry in the stratosphere. The flight model has been developed. A balloon-borne VLBI station will be launched to establish interferometric fringes with ground-based VLBI stations distributed over the Japanese islands at an observing frequency of approximately 20?GHz as the first step. This paper describes the system design and development of a series of observing instruments and bus systems. In addition to the advantages of avoiding the atmospheric effects of absorption and fluctuation in high frequency radio observation, the mobility of a station can improve the sampling coverage (“uv-coverage”) by increasing the number of baselines by the number of ground-based counterparts for each observation day. This benefit cannot be obtained with conventional arrays that solely comprise ground-based stations. The balloon-borne VLBI can contribute to a future progress of research fields such as black holes by direct imaging.     

ITRS/GCRS transformation: Uncertainty propagation analysis and short-term modelling of IAU 2006/2000A developments

A study of the uncertainty propagation in ITRS/GCRS transformation is presented in this work. General law of propagation of variances is applied to the ITRS/GCRS transformation matrix, deriving the analytical expressions to compute GCRS position uncertainty. This evaluation is based on EOP uncertainties provided by IERS long-term series and formal uncertainties of ITRS-compatible coordinates. Numerical results for the period 1998–2016 are shown and discussed for ITRS positions in different altitudes and latitudes, providing graphical and numerical insights of the mapping of EOP uncertainties to transformed coordinates.Eventually, an analysis of short-term evolution of the Celestial Intermediate Pole coordinates in the GCRS provided by the IAU2006/2000A precession-nutation model is carried out, in order to assess the feasibility to potentially broadcast these parameters in GNSS navigation message. This approach would facilitate the dissemination of terrestrial-celestial transformation parameters for real time users, given that polar motion and UT1-UTC are already foreseen in GPS interface specification. The results presented in this work confirm the feasibility of this idea.     

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.     

Application of time-variable process noise in terrestrial reference frames determined from VLBI data

In recent years, Kalman filtering has emerged as a suitable technique to determine terrestrial reference frames (TRFs), a prime example being JTRF2014. The time series approach allows variations of station coordinates that are neither reduced by observational corrections nor considered in the functional model to be taken into account. These variations are primarily due to non-tidal geophysical loading effects that are not reduced according to the current IERS Conventions (2010). It is standard practice that the process noise models applied in Kalman filter TRF solutions are derived from time series of loading displacements and account for station dependent differences. So far, it has been assumed that the parameters of these process noise models are constant over time. However, due to the presence of seasonal and irregular variations, this assumption does not truly reflect reality. In this study, we derive a station coordinate process noise model allowing for such temporal variations. This process noise model and one that is a parameterized version of the former are applied in the computation of TRF solutions based on very long baseline interferometry data. In comparison with a solution based on a constant process noise model, we find that the station coordinates are affected at the millimeter level.     

Distance variations between laser ranging stations derived from LAGEOS tracking data

The results obtained by the DGFI analysis center for the MERIT campaign have demonstrated that from LAGEOS laser ranging data the relative positions of the globally distributed tracking stations can be evaluated with an accuracy of a few centimeters. It therefore seemed promising to investigate the stability of tracking station coordinates over a period of several years in order to detect interstation distance variations due to global horizontal deformations of the Earth's surface. The investigation focussed on a selected number of stations which have a good tracking data coverage over the 1980 – 1984 time span. Annual geocentric coordinates of these stations were estimated. From the coordinates spherical distances between stations on five different lithospheric plates were derived and the changes in length from year to year were investigated.     

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.     

Automated and dynamic scheduling for geodetic VLBI – A simulation study for AuScope and global networks

As we move into the next era of geodetic VLBI, the scheduling process is one focus for improvement in terms of increased flexibility and the ability to react with changing conditions. A range of simulations were conducted to ascertain the impact of scheduling on geodetic results such as Earth Orientation Parameters (EOPs) and station coordinates. The potential capabilities of new automated scheduling modes were also simulated, using the so-called ‘dynamic scheduling’ technique. The primary aim was to improve efficiency for both cost and time without losing geodetic precision, particularly to maximise the uses of the Australian AuScope VLBI array.We show that short breaks in observation will not significantly degrade the results of a typical 24?h experiment, whereas simply shortening observing time degrades precision exponentially. We also confirm the new automated, dynamic scheduling mode is capable of producing the same standard of result as a traditional schedule, with close to real-time flexibility. Further, it is possible to use the dynamic scheduler to augment the 3 station Australian AuScope array and thereby attain EOPs of the current global precision with only intermittent contribution from 2 additional stations. We thus confirm automated, dynamic scheduling bears great potential for flexibility and automation in line with aims for future continuous VLBI operations.