Design and experiment of onboard laser time transfer in Chinese Beidou navigation satellites

High-precision time synchronization between satellites and ground stations plays the vital role in satellite navigation system. Laser time transfer (LTT) technology is widely recognized as the highest accuracy way to achieve time synchronization derived from satellite laser ranging (SLR) technology. Onboard LTT payload has been designed and developed by Shanghai Astronomical Observatory, and successfully applied to Chinese Beidou navigation satellites. By using the SLR system, with strictly controlling laser firing time and developing LTT data processing system on ground, the high precise onboard laser time transfer experiment has been first performed for satellite navigation system in the world. The clock difference and relative frequency difference between the ground hydrogen maser and space rubidium clocks have been obtained, with the precision of approximately 300 ps and relative frequency stability of 10E−14. This article describes the development of onboard LTT payload, introduces the principle, system composition, applications and LTT measuring results for Chinese satellite navigation system.     

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.     

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.     

Spin axis orientation of Ajisai determined from Graz 2 kHz SLR data

The Graz 2 kHz Satellite Laser Ranging (SLR) measurements allow determination of the spin axis orientation of the geodetic satellite Ajisai. The high repetition rate of the laser makes it possible to determine the epoch time when the laser is pointing directly between two corner cube reflector (CCR) rings of the satellite. Identification of many such events during a few (up to 3) consecutive passes allows to state the satellite orientation in the celestial coordinate system. Six years of 2 kHz SLR data (October 2003–October 2009) delivered 331 orientation values which clearly show precession of the axis along a cone centered at 14^h56^m2.8^s in right ascension and 88.512° in declination (J2000.0 celestial reference frame) and with an half-aperture angle θ of 1.405°. The spin axis precesses with a period of 117 days, which is equal to the period of the right ascension of the ascending node of Ajisai’s orbit. We present a model of the axis precession which allows prediction of the satellite orientation – necessary for the envisaged laser time transfer via Ajisai mirrors.     

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

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

Estimation of the Love and Shida numbers: h2, l2 using SLR data for the low satellites

In this paper we present results for the global elastic parameters: Love number h₂ and Shida number l₂ derived from the analysis of Satellite Laser Ranging (SLR) data. SLR data for the two low satellites STELLA (H = 800 km) and STARLETTE (H = 810 km) observed during 2.5 years from January 3, 2005 until July 1, 2007 with 18 globally distributed ground stations were analyzed. The analysis was done separately for the two satellites. We do a sequential analysis and study the stability and convergence of the estimates as a function of length of the data set used.     

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.     

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.     

22 Years of AJISAI spin period determination from standard SLR and kHz SLR data

Satellite Laser Ranging (SLR) is a powerful and efficient technique to measure spin parameters of satellites equipped with corner cube reflectors. We obtained spin period determination of the satellite AJISAI from SLR data only: 17246 pass-by-pass estimates from standard 1–15 Hz SLR data (14/Aug/1986–30/Dec/2008) and 1444 pass-by-pass estimates (9/Oct/2003–30/Dec/2008) from data of the first 2 kHz SLR system from Graz, Austria. A continuous history of the slowing down of AJISAI spin is derived from frequency analysis, and corrected for the apparent effects. The apparent corrections, elaborated here, allowed very accurate determination of AJISAI initial spin period: 1.4855 ± 0.0007 [s]. The paper identifies also non-gravitational effects as a source of the periodical changes in the rate of slowing down of the satellite.     

Photon counting detector for space debris laser tracking and lunar laser ranging

We are reporting on a design, construction and performance of solid state photon counting detector package which has been designed for laser tracking of space debris. The detector has been optimized for top photon detection efficiency and detection delay stability. The active area of the commercially available avalanche photodiode manufactured on Si (SAP500 supplied by Laser Components, Inc.) is circular with a diameter of 500 μm. The newly designed control circuit enables to operate the detection sensor at a broad range of biases 5–50 V above its breakdown voltage of 125 V. This permits to select a right trade-off between photon detection efficiency, timing resolution and dark count rate. The photon detection efficiency exceeds 70% at the wavelength of 532 nm. This is the highest photon detection efficiency ever reported for such a device. The timing properties of the detector have been investigated in detail. The timing resolution is better than 80 ps r.m.s, the detection delay is stable within units of picoseconds over several hours of operation. The detection delay stability in a sense of time deviation of 800 fs has been achieved. The temperature change of the detection delay is 0.5 ps/K. The detector has been tested as an echo signal detector in laser tracking of space debris at the satellite laser station in Graz, Austria. Its application in lunar laser ranging is under consideration by several laser stations.     

Monitoring atomic clocks on board GNSS satellites

A method for monitoring atomic clocks on board Global Navigation Satellites System (GNSS) satellites is described to address the issue of clock related signal integrity in safety–critical applications of GNSS. The carrier-phase time transfer is employed in the clock monitoring method which enables tight tracking of the satellite onboard clocks and thus improves detectability of clock anomalies. Detecting onboard clock anomalies requires the ability to monitor clocks in real time, and a Kalman filter can then be utilized to estimate the phase offsets between the satellite clocks and ground clocks. This study, using the difference between the measured and predicted phase offset as a test statistic, sets a threshold for clock anomalies based on the prediction interval approach. Finally the validity of the monitoring method is examined by processing a set of real GNSS data that includes two recent incidents of clock anomalies in GNSS satellites.     

Achievable debris orbit prediction accuracy using laser ranging data from a single station

Earlier studies have shown that an orbit prediction accuracy of 20 arc sec ground station pointing error for 1–2 day predictions was achievable for low Earth orbit (LEO) debris using two passes of debris laser ranging (DLR) data from a single station, separated by about 24 h. The accuracy was determined by comparing the predicted orbits with subsequent tracking data from the same station. This accuracy statement might be over-optimistic for other parts of orbit far away from the station. This paper presents the achievable orbit prediction accuracy using satellite laser ranging (SLR) data of Starlette and Larets under a similar data scenario as that of DLR. The SLR data is corrupted with random errors of 1 m standard deviation so that its accuracy is similar to that of DLR data. The accurate ILRS Consolidated Prediction Format orbits are used as reference to compute the orbit prediction errors. The study demonstrates that accuracy of 20 arc sec for 1–2 day predictions is achievable.     

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

Design and performances of laser retro-reflector arrays for Beidou navigation satellites and SLR observations

Beidou is the regional satellite navigation system in China, consisting of three kinds of orbiting satellites, MEO, GEO and IGSO, with the orbital altitudes of 21500–36000 km. For improving the accuracy of satellites orbit determination, calibrating microwave measuring techniques and providing better navigation service, all Beidou satellites are equipped with laser retro-reflector arrays (LRAs) to implement high precision laser ranging. The paper presents the design of LRAs for Beidou navigation satellites and the method of inclined installation of LRAs for GEO satellites to increase the effective reflective areas for the regional ground stations. By using the SLR system, the observations for Beidou satellites demonstrated a precision of centimeters. The performances of these LRAs on Beidou satellites are very excellent.     

Estimation of the elastic Earth parameters (h2, l2) using SLR data

We present results for the global elastic parameters h₂ and l₂ derived from the analysis of Satellite Laser Ranging (SLR) data. SLR data for the two satellites LAGEOS 1 and LAGEOS 2 observed during 2.5 years from January 3, 2005 until July 1, 2007 with 18 globally distributed ground stations were analysed using different approaches. The analysis was done separately for the two satellites and approaches to estimate the two elastic parameters independently and together were performed. We do a sequential analysis and study the stability of the estimates as a function of length of the data set used. The adjusted final values for h₂ equal to 0.6151 ± 0.0008 and 0.6152 ± 0.0008, and those for l₂ equal to 0.0886 ± 0.0003 and 0.0881 ± 0.0003 for LAGEOS 1 and LAGEOS 2 tracking data are compared to other independently derived estimates. These parameters and their errors achieve stability at about the 24 and 27 month time interval for h₂ and l₂, respectively.     

The impact of orbital errors on the estimation of satellite clock errors and PPP

Precise satellite orbit and clocks are essential for providing high accuracy real-time PPP (Precise Point Positioning) service. However, by treating the predicted orbits as fixed, the orbital errors may be partially assimilated by the estimated satellite clock and hence impact the positioning solutions. This paper presents the impact analysis of errors in radial and tangential orbital components on the estimation of satellite clocks and PPP through theoretical study and experimental evaluation. The relationship between the compensation of the orbital errors by the satellite clocks and the satellite-station geometry is discussed in details. Based on the satellite clocks estimated with regional station networks of different sizes (∼100, ∼300, ∼500 and ∼700 km in radius), results indicated that the orbital errors compensated by the satellite clock estimates reduce as the size of the network increases. An interesting regional PPP mode based on the broadcast ephemeris and the corresponding estimated satellite clocks is proposed and evaluated through the numerical study. The impact of orbital errors in the broadcast ephemeris has shown to be negligible for PPP users in a regional network of a radius of ∼300 km, with positioning RMS of about 1.4, 1.4 and 3.7 cm for east, north and up component in the post-mission kinematic mode, comparable with 1.3, 1.3 and 3.6 cm using the precise orbits and the corresponding estimated clocks. Compared with the DGPS and RTK positioning, only the estimated satellite clocks are needed to be disseminated to PPP users for this approach. It can significantly alleviate the communication burdens and therefore can be beneficial to the real time applications.     

New results on spin determination of nanosatellite BLITS from High Repetition Rate SLR data

The nanosatellite BLITS (Ball Lens In The Space) demonstrates a successful design of the new spherical lens type satellite for Satellite Laser Ranging (SLR). The spin parameters of the satellite were calculated from more than 1000 days of SLR data collected from 6 High Repetition Rate (HRR) systems: Beijing, Changchun, Graz, Herstmonceux, Potsdam, Shanghai.     

Laser measurements to space debris from Graz SLR station

In order to test laser ranging possibilities to space debris objects, the Satellite Laser Ranging (SLR) Station Graz installed a frequency doubled Nd:YAG pulse laser with a 1 kHz repetition rate, a pulse width of 10 ns, and a pulse energy of 25 mJ at 532 nm (on loan from German Aerospace Center Stuttgart – DLR). We developed and built low-noise single-photon detection units to enable laser ranging to targets with inaccurate orbit predictions, and adapted our standard SLR software to include a few hundred space debris targets. With this configuration, we successfully tracked – within 13 early-evening sessions of each about 1.5 h – 85 passes of 43 different space debris targets, in distances between 600 km and up to more than 2500 km, with radar cross sections from >15 m² down to <0.3 m², and measured their distances with an average precision of about 0.7 m RMS.     

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.