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.     

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).     

Impact of Jason-2/T2L2 Ultra-Stable-Oscillator Frequency Model on DORIS stations coordinates and Earth Orientation Parameters

Doppler Orbitography and Radiopositioning Integrated by Satellite (DORIS) is a tracking technique based on a one-way ground to space Doppler link. For Low Earth Orbit (LEO) satellites, DORIS shows a robust capability in terms of data coverage and availability, due to a wide and well-distributed ground network, where data are made available by the International Doris Service (IDS). However, systematic errors remain in the DORIS data, such as instabilities of the on-board clock due to radiation encountered in space, which limit the accurate determination of station positions.The DORIS on-board clock frequency stability is degraded by the increased radiation found in the region of the South Atlantic Anomaly (SAA) and has been shown to degrade station position estimation. This paper introduces a new model correction to the DORIS data for the frequency of the Jason-2 Ultra Stable Oscillator (USO), derived from the Time Transfer by Laser Link (T2L2) experiment (Belli and Exertier, 2018). We show that a multi-satellite DORIS solution including this T2L2-corrected data applied to the frequency modelling for The DORIS data, improves the estimation of station coordinates. We show the tie residuals with respect to collocated GPS stations are improved by several millimeters. We also demonstrate that the 117-day (Jason-2) draconitic signal in the geophysical parameters is reduced, implying that the origin of this signal is not just solar radiation pressure mis-modeling, but also radiation-induced clock perturbations on the Jason-2 DORIS Ultra-Stable-Oscillator (USO). Finally we demonstrate through comparisons with the International Earth Rotations and Reference Systems Service (IERS) C04 series for Earth Orientation Parameters (EOP), that the estimation of EOP is improved in both a Jason-2 DORIS-only and a multi-satellite DORIS solution for EOP.     

In-orbit assessment of laser retro-reflector efficiency onboard high orbiting satellites

The navigation and geodetic satellites that orbit the Earth at altitudes of approximately 20,000 km are tracked routinely by many of the Satellite Laser Ranging (SLR) stations of the International Laser Ranging Service (ILRS). In order to meet increasing demands on SLR stations for daytime and nighttime observations, any new mission needs to ensure a strong return signal so that the target is easily acquirable. The ILRS has therefore set a minimum effective cross-section of 100 million square metres for the on-board laser retro-reflector arrays (LRAs) and further recommends the use of ‘uncoated’ cubes in the arrays. Given the large number of GNSS satellites that are currently supported by SLR, it is informative to make an assessment of the relative efficiencies of the various LRAs employed. This paper uses the laser ranging observations themselves to deduce and then compare the efficiencies of the LRAs on the COMPASS-M1 navigation satellite, two satellites from the GPS and three from the GLONASS constellations, the two GIOVE test satellites from the upcoming Galileo constellation, the two Etalon geodetic spheres and the geosynchronous communications test satellite, ETS-8. All the LRAs on this set of satellites employ back-coated retro-reflector cubes, except those on the COMPASS-M1 and ETS-8 vehicles which are uncoated. A measure of return signal strength, and thus of LRA-efficiency, is calculated using the laser-range full-rate data archive from 2007 to 2010, scaled to remove the effects of variations in satellite range, atmospheric attenuation and retro-reflector target total surface area. Observations from five SLR stations are used in this study; they are Herstmonceux (UK), Yarragadee (Australia), Monument Peak and McDonald (USA) and Wettzell (Germany). Careful consideration is given to the treatment of the observations from each station in order to take account of local working practices and system upgrades. The results show that the uncoated retro-reflector cubes offer significant improvements in efficiency.     

Quasi-biennial oscillation in GPS VTEC measurements

The quasi-biennial oscillation, QBO, a well known periodicity in the equatorial stratospheric zonal winds, is also found in ionospheric parameters and in solar and geomagnetic activity indices. Many authors speculated about the link between the QBO in solar and geomagnetic activity and the QBO in atmospheric parameters. In this work we analyze the presence of the QBO in the ionosphere using the Vertical Total Electron Content (VTEC) values obtained from Global Navigation Satellite System (GNSS) measurements during the period 1999–2012. In particular, we used IONEX files, i.e. the International GNSS Service (IGS) ionospheric products. IONEX provide VTEC values around the world at 2-h intervals. From these data we compute global and zonal averages of VTEC at different local times at mid and equatorial geomagnetic latitudes. VTEC and Extreme Ultra Violet (EUV) solar flux time series are analyzed using a wavelet multi resolution analysis. In all cases the QBO is detected among other expected periodicities.     

Impact of loading effects on determination of the International Terrestrial Reference Frame

The International Terrestrial Reference Frame (ITRF), as a realization of the International Terrestrial Reference System (ITRS), is represented by a set of station positions and linear velocities. They are intended to be used as regularized coordinates to which some corrections should be added to access instantaneous coordinates. The latest ITRS realization is the ITRF2005, which has integrated time series of station positions to form long-term solutions for the four space geodetic techniques. Currently, a purely linear model is used to parameterize station displacements in the estimation process, plus occasional discontinuities in case of earthquakes or equipment changes. However the input data have been derived without applying surface loading models and so surface loading effects are supposed to be embedded in the coordinates as measured quantities. We evaluate the effect of applying a posteriori loading corrections, which include the effect of atmospheric, non-tidal ocean, and continental water loading, to time series of positions estimated from Satellite Laser Ranging (SLR), Very Long Baseline Interferometry (VLBI), and Global Positioning System (GPS) data. We notice that they reduce about 50% or more of the annual signals in the translation and scale parameter time series of the SLR and VLBI techniques, except in SLR Z translation. In general, the estimated secular frame definition is negligibly affected and estimated positions and velocities are not significantly modified for stations that have accumulated a large number of observations. A multi-technique combination of such derived frames allows concluding that, for some cases, loading model corrections might degrade co-located station coordinates almost as much as they benefit them. However, most significant improvement of the estimated secular coordinates is observed for stations with less than 100 estimated positions as demonstrated with a multi-technique combination.     

Contributions of the French Institut Géographique National (IGN) to the International DORIS Service

DORIS is one of the four geodetic space techniques participating to the Global Geodetic Observing System (GGOS). Since the early development of this system, the Institut Géographique National played a specific and active role. Within, the International DORIS Service (IDS), IGN is in a particular position. While it is responsible for the installation and the maintenance of the DORIS ground tracking stations, it also handles one of the two IDS data center delivering DORIS data and products and has been an IDS Analysis Center for years, providing all possible IDS products, in particular the latest ignwd08 time series in preparation of ITRF2008. This paper explains the different aspects of the IGN contribution to IDS from an historical point of view, presents current activities and scientific results and provides a perspective for future activities. Recent DORIS results show a 10 mm precision or better when more than four DORIS satellites are available. Comparisons between recent DORIS solutions (ign07d02 and ign09d02) and past ITRF realizations show that errors are shared between the DORIS and the ITRF realizations. Some problems related to DORIS data processing are also discussed and possible ways to solve them in the future are discussed. In particular, we can now reject the tropospheric origin of the problem detected in the Envisat data after the software upload of October 12, 2004. A few applications in geodesy (terrestrial reference frame, Earth’s polar motion) and geophysics are also discussed as a natural extension of these service-type activities.     

Improving DORIS geocenter time series using an empirical rescaling of solar radiation pressure models

Even if Satellite Laser Ranging (SLR) remains the fundamental technique for geocenter monitoring, DORIS can also determine this geophysical parameter. Gobinddass et al. (2009) found that part of the systematic errors at 118 days and 1 year can be significantly reduced by rescaling the current solar radiation pressure models using satellite-dependent empirical models. Here we extend this study to all DORIS satellites and propose a complete set of empirical solar radiation parameter coefficients. A specific problem related to SPOT-5 solar panel realignment is also detected and explained. New DORIS geocenter solutions now show a much better agreement in amplitude with independent SLR solutions and with recent geophysical models. Finally, the impact of this refined DORIS data strategy is discussed in terms of Z-geocenter monitoring as well as for other geodetic products (altitude of high latitude station such as Thule in Greenland) and Precise Orbit Determination. After reprocessing the full 1993.0-2008.0 DORIS data set, we confirm that the proposed strategy allows a significant reduction of systematic errors at periods of 118 days and 1 year (up to 20 mm), especially for the most recent data after 2002.5, when more DORIS satellites are available for geodetic purposes.     

Precision orbit determination for TOPEX/Poseidon using TDRSS doppler tracking data

Precision orbit determination on the TOPEX/Poseidon (T/P) altimeter satellite is now being routinely achieved with sub-5cm radial and sub-15 cm total positioning accuracy using state-of-the-art modeling with precision tracking provided by a combination of: (a) global Satellite Laser Ranging (SLR) and Doppler Orbitography and Radiopositioning Integrated by Satellite (DORIS), or (b) the Global Positioning System (GPS) Constellation which provides pseudo-range and carrier phase observations. The geostationary Tracking and Data Relay Satellite System (TDRSS) satellites are providing the operational tracking and communication support for this mission. The TDRSS Doppler data are of high precision (0.3 mm/s nominal noise levels). Unlike other satellite missions supported operationally by TDRSS, T/P has high quality independent tracking which enables absolute orbit accuracy assessments. In addition, the T/P satellite provides extensive geometry for positioning a satellite at geostationary altitude, and thus the TDRSS-T/P data provides an excellent means for determining the TDRS orbits. Arc lengths of 7 and 10 days with varying degrees of T/P spacecraft attitude complexity are studied. Sub-meter T/P total positioning error is achieved when using the TDRSS range-rate data, with radial orbit errors of 10.6 cm and 15.5 cm RMS for the two arcs studied. Current limitations in the TDRSS precision orbit determination capability include mismodeling of numerous TDRSS satellite-specific dynamic and electronic effects, and in the inadequate treatment of the propagation delay and bending arising from the wet troposphere and ionosphere.     

LEO enhanced Global Navigation Satellite System (LeGNSS) for real-time precise positioning services

Global Navigation Satellite System (GNSS) has been widely used in many geosciences areas with its Positioning, Navigation and Timing (PNT) service. However, GNSS still has its own bottleneck, such as the long initialization period of Precise Point Positioning (PPP) without dense reference network. Recently, the concept of PNTRC (Positioning, Navigation, Timing, Remote sensing and Communication) has been put forward, where Low Earth Orbit (LEO) satellite constellations are recruited to fulfill diverse missions. In navigation aspect, a number of selected LEO satellites can be equipped with a transmitter to transmit similar navigation signals to ground users, so that they can serve as GNSS satellites but with much faster geometric change to enhance GNSS capability, which is named as LEO constellation enhanced GNSS (LeGNSS). As a result, the initialization time of PPP is expected to be shortened to the level of a few minutes or even seconds depending on the number of the LEO satellites involved. In this article, we simulate all the relevant data from June 8th to 14th, 2014 and investigate the feasibility of LeGNSS with the concentration on the key issues in the whole data processing for providing real-time PPP service based on a system configuration with fourteen satellites of BeiDou Navigation Satellite System (BDS), twenty-four satellites of the Global Positioning System (GPS), and sixty-six satellites of the Iridium satellite constellations. At the server-end, Precise Orbit Determination (POD) and Precise Clock Estimation (PCE) with various operational modes are investigated using simulated observations. It is found out that GNSS POD with partial LEO satellites is the most practical mode of LeGNSS operation. At the user-end, the Geometry Dilution Of Precision (GDOP) and Signal-In-Space Ranging Error (SISRE) are calculated and assessed for different positioning schemes in order to demonstrate the performance of LeGNSS. Centimeter level SISRE can be achieved for LeGNSS.     

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.     

VLBI observations to the APOD satellite

The APOD (Atmospheric density detection and Precise Orbit Determination) is the first LEO (Low Earth Orbit) satellite in orbit co-located with a dual-frequency GNSS (GPS/BD) receiver, an SLR reflector, and a VLBI X/S dual band beacon. From the overlap statistics between consecutive solution arcs and the independent validation by SLR measurements, the orbit position deviation was below 10?cm before the on-board GNSS receiver got partially operational. In this paper, the focus is on the VLBI observations to the LEO satellite from multiple geodetic VLBI radio telescopes, since this is the first implementation of a dedicated VLBI transmitter in low Earth orbit. The practical problems of tracking a fast moving spacecraft with current VLBI ground infrastructure were solved and strong interferometric fringes were obtained by cross-correlation of APOD carrier and DOR (Differential One-way Ranging) signals. The precision in X-band time delay derived from 0.1?s integration time of the correlator output is on the level of 0.1?ns. The APOD observations demonstrate encouraging prospects of co-location of multiple space geodetic techniques in space, as a first prototype.     

Status of the T2L2/Jason2 Experiment

The T2L2 (Time Transfer by Laser Link) project, developed by CNES and OCA will permit the synchronization of remote ultra stable clocks and the determination of their performances over intercontinental distances. The principle of the experiment derives from Satellite Laser Ranging (SLR) technology with dedicated space equipment. T2L2 was accepted in 2005 to be on board the Jason2 altimetry satellite. The payload consists of both event timer and photo detection modules. The system uses the ultra-stable quartz oscillator of DORIS as on-board reference clock on one hand, and the Laser Reflector Array, making T2L2 a real two-way time transfer system on the other hand. The expected time stability of the T2L2 instrument (detection and timing), referenced by the DORIS oscillator and including all internal error sources should be at the level of 10–12 ps at 1 s and <1 ps at 1000 s. The metrological specifications of T2L2 should permit to maintain a precision of 1 to a few ps when measuring the phase of a clock during around 1000 seconds.     

The impact on tsunami detection from space using GNSS-reflectometry when combining GPS with GLONASS and Galileo

The devastating Sumatra tsunami in 2004 demonstrated the need for a tsunami early warning system in the Indian Ocean. Such a system has been installed within the German-Indonesian Tsunami Early Warning System (GITEWS) project. Tsunamis are a global phenomenon and for global observations satellites are predestined. Within the GITEWS project a feasibility study on a future tsunami detection system from space has therefore been carried out. The Global Navigation Satellite System Reflectometry (GNSS-R) is an innovative way of using GNSS signals for remote sensing. It uses ocean reflected GNSS signals for sea surface altimetry. With a dedicated Low Earth Orbit (LEO) constellation of satellites equipped with GNSS-R receivers, densely spaced sea surface height measurements could be established to detect tsunamis. Some general considerations on the geometry between LEO and GNSS are made in this simulation study. It exemplary analyzes the detection performance of a GNSS-R constellation at 900 km altitude and 60° inclination angle when applied to the Sumatra tsunami as it occurred in 2004. GPS is assumed as signal source and the combination with GLONASS and Galileo signals is investigated. It can be demonstrated, that the combination of GPS and Galileo is advantageous for constellations with few satellites while the combination with GLONASS is preferable for constellations with many satellites. If all three GNSS are combined, the best detection performance can be expected for all scenarios considered. In this case an 18 satellite constellation will detect the Sumatra tsunami within 17 min with certainty, while it takes 53 min if only GPS is considered.     

Spin parameters of High Earth Orbiting satellites Etalon-1 and Etalon-2 determined from kHz Satellite Laser Ranging data

The high repetition rate Satellite Laser Ranging (SLR) system developed in Graz, Austria, measures ranges to the High Earth Orbiting satellites Etalon-1 and Etalon-2 with the millimeter accuracy. The 2 kHz repetition rate of the laser and the relatively high return rates allow to use the SLR data to calculate the spin parameters of the Etalon satellites. The analysis of the 10 years (October 2003–September 2013) of the SLR data gives trends of the spin axes orientation (J2000 Inertial Reference Frame):     

Electron density retrieval from occulting GNSS signals using a gradient-aided inversion technique

In the coming years, opportunities for remote sensing of electron density in the Earth’s ionosphere will expand with the advent of Galileo, which will become part of the global navigation satellite system (GNSS). Methods for accurate electron density retrieval from radio occultation data continue to improve. We describe a new method of electron density retrieval using total electron content measurements obtained in low Earth orbit. This method can be applied to data from dual-frequency receivers tracking the GPS or Galileo transmitters. This simulation study demonstrates that the method significantly improves retrieval accuracy compared to the standard Abel inversion approach that assumes a spherically symmetric ionosphere. Our method incorporates horizontal gradient information available from global maps of Total Electron Content (TEC), which are available from the International GNSS Service (IGS) on a routine basis. The combination of ground and space measurements allows us to improve the accuracy of electron density profiles near the occultation tangent point in the E and F regions of the ionosphere.     

Numerical modeling of a Global Navigation Satellite System in a general relativistic framework

In this article we model a Global Navigation Satellite System (GNSS) in a Schwarzschild space–time, as a first approximation of the relativistic geometry around the Earth. The closed time-like and scattering light-like geodesics are obtained analytically, describing respectively trajectories of satellites and electromagnetic signals. We implement an algorithm to calculate Schwarzschild coordinates of a GNSS user who receives proper times sent by four satellites, knowing their orbital parameters; the inverse procedure is implemented to check for consistency. The constellation of satellites therefore realizes a geocentric inertial reference system with no a priori realization of a terrestrial reference frame. We perform a simulation of position determination and show that the determination of the four coordinates with a 25–32 digit accuracy takes only around 60 ms. Effects of non-gravitational perturbations on positioning errors are assessed, and methods to reduce them are sketched. In particular, inter-links between satellites could greatly enhance stability and accuracy of the positioning system. Effects of gravitational perturbations are omitted in this paper in order to make a clearer comparison between the relativistic and non-relativistic scheme, but they will be included in subsequent work. We believe that the final algorithm will be a serious alternative to the usual post-Newtonian scheme.     

Integer phase ambiguity estimation in next-generation geodetic Very Long Baseline Interferometry

Next-generation Very Long Baseline Interferometry (VLBI) system designs are aiming at 1 mm global position accuracy. In order to achieve this, it is not only necessary to deploy improved VLBI systems, but also to develop analysis strategies that take full advantage of the observations taken. Since the new systems are expected to incorporate four independent radio frequency bands, it should be feasible to resolve phase ambiguities directly from post-correlation data, providing roughly an order of magnitude improvement in precision of the delay observable. As the unknown ambiguities are of integer nature, it is discussed here how they the can be resolved analytically using algorithms which have been developed for Global Navigation Satellite System (GNSS) applications. Furthermore, it will be shown that ionosphere contribution and source structure effects, so-called core-shifts, can be solved simultaneously with the delay, which is the main geodetic observable for follow-on analysis. In order to verify the proposed algorithm, simulated observations were created using parameters from actual design studies. It is shown that, even in the case of low signal-to-noise ratio observations, reliable phase ambiguity resolution can be achieved and it is discussed how the integer ambiguity recovery depends on the number of observations and signal-to-noise ratio.     

Seasonal variations of storm-time TEC at European middle latitudes

Global Navigation Satellite System (GNSS) measurements of the Total Electron Content (TEC) from local (Dourbes, 50.1°N, 04.6°E) and European IGS (International GNSS Service) stations were used to obtain the TEC changes during the geomagnetic storms of the latest solar activity cycle. A common epoch analysis, with respect to geomagnetic storm intensity, season, and latitude, was performed on data representing nearly 300 storm events. In general, the storm-time behaviour of TEC shows clear positive and negative phases, relative to the non-storm (median) behaviour, with amplitudes that tend to increase during more intense storms. The most pronounced positive phase is observed during winter, while the strongest and yet shortest negative phase is detected during equinox. Average storm-time patterns in the TEC behaviour are deduced for potential use in ionosphere prediction services.     

基于LabWindows/CVI的IGS实时数据自动下载软件开发

随着全球导航卫星系统（GNSS）的日益发展,国际GNSS服务组织（IGS）所产生的各项数据以及产品日趋精确,应用也非常广泛。本文基于Labwindows/CVI开发平台,设计了从IGS的FTP服务器下载实时数据或产品的自动运行软件。该软件目前运行稳定、可靠,可以为IGS数据或产品的使用者提供自动化服务。