A new approach to the telescope operation method for satellite tracking using a time synchronization technique

For the development of a telescope that is capable of precisely tracking satellites and high-speed operation such as satellite laser ranging, a special method of telescope operation is required. This study aims to propose a new telescope operation method and system configuration for the independent development of a mount and an operation system which includes the host computer. Considering that the tracking of a satellite is performed in real time, communication and synchronization between the two independent subsystems are important. Therefore, this study applied the concept of time synchronization, which is used in various fields of industry, to the communication between the command computer and the mount. In this case, communication delays do not need to be considered in general, and it is possible to cope with data loss. Above all, when the mount is replaced in the future, only the general communication interface needs to be modified, and thus, it is not limited by replacement in terms of the overall system management. The performance of the telescope operation method developed in this study was verified by applying the method to the first mobile SLR system in Korea. This study is significant in that it proposed a new operation method and system configuration, to which the concept of time synchronization was applied, for the observation system that requires an optical telescope.     

The crustal dynamics data information system: A resource to support scientific analysis using space geodesy

Since 1982, the Crustal Dynamics Data Information System (CDDIS) has supported the archive and distribution of geodetic data products acquired by the National Aeronautics and Space Administration (NASA) as well as national and international programs. The CDDIS provides easy, timely, and reliable access to a variety of data sets, products, and information about these data. These measurements, obtained from a global network of nearly 650 instruments at more than 400 distinct sites, include DORIS (Doppler Orbitography and Radiopositioning Integrated by Satellite), GNSS (Global Navigation Satellite System), SLR and LLR (Satellite and Lunar Laser Ranging), and VLBI (Very Long Baseline Interferometry). The CDDIS data system and its archive have become increasingly important to many national and international science communities, particularly several of the operational services within the International Association of Geodesy (IAG) and its observing system the Global Geodetic Observing System (GGOS), including the International DORIS Service (IDS), the International GNSS Service (IGS), the International Laser Ranging Service (ILRS), the International VLBI Service for Geodesy and Astrometry (IVS), and the International Earth rotation and Reference frame Service (IERS). Investigations resulting from the data and products available through the CDDIS support research in many aspects of Earth system science and global change. Each month, the CDDIS archives more than one million data and derived product files totaling over 90 Gbytes in volume. In turn, the global user community downloads nearly 1.2 Tbytes (over 10.5 million files) of data and products from the CDDIS each month. The requirements of analysts have evolved since the start of the CDDIS; the specialized nature of the system accommodates the enhancements required to support diverse data sets and user needs. This paper discusses the CDDIS, including background information about the system and its user communities, archive contents, available metadata, and future plans.     

Precision orbit determination performance for CryoSat-2

In this paper we discuss our efforts to perform precision orbit determination (POD) of CryoSat-2 which depends on Doppler and satellite laser ranging tracking data. A dynamic orbit model is set-up and the residuals between the model and the tracking data is evaluated. The average r.m.s. of the 10?s averaged Doppler tracking pass residuals is approximately 0.39?mm/s; and the average of the laser tracking pass residuals becomes 1.42?cm. There are a number of other tests to verify the quality of the orbit solution, we compare our computed orbits against three independent external trajectories provided by the CNES. The CNES products are part of the CryoSat-2 products distributed by ESA. The radial differences of our solution relative to the CNES precision orbits shows an average r.m.s. of 1.25?cm between Jun-2010 and Apr-2017. The SIRAL altimeter crossover difference statistics demonstrate that the quality of our orbit solution is comparable to that of the POE solution computed by the CNES. In this paper we will discuss three important changes in our POD activities that have brought the orbit performance to this level. The improvements concern the way we implement temporal gravity accelerations observed by GRACE; the implementation of ITRF2014 coordinates and velocities for the DORIS beacons and the SLR tracking sites. We also discuss an adjustment of the SLR retroreflector position within the satellite reference frame. An unexpected result is that we find a systematic difference between the median of the 10 s Doppler tracking residuals which displays a statistically significant pattern in the South Atlantic Anomaly (SSA) area where the median of the velocity residuals varies in the range of ?0.15 to +0.15?mm/s.     

10 Years of LAGEOS-1 and 15 years of LAGEOS-2 spin period determination from SLR data

Satellite Laser Ranging (SLR) stations measure distance to the satellites equipped with Corner Cube Reflectors (CCRs). These range measurements contain information about spin parameters of the spacecraft. In this paper we present results of spin period determination of two passive satellites from SLR data only: 10 years of LAGEOS-1 (10426 values), and 15 years of LAGEOS-2 (15580 values). The measurements have been made by standard 10 Hz SLR systems and the first 2 kHz SLR system from Graz (Austria). The obtained data allowed calculation of the initial spin period of the satellites: 0.61 s for LAGEOS-1 and 0.906 s for LAGEOS-2. Long time series of the spin period values show that the satellite’s slowing down rate is not constant but is oscillating with a period of 846 days for LAGEOS-1 and 578 days for LAGEOS-2. The results presented here definitely prove that the SLR is a very efficient technique able to measure spin period of the geodetic satellites.     

Time Transfer by Laser Link: Data analysis and validation to the ps level

The Time Transfer by Laser Link (T2L2) is a very high resolution time transfer technique based on the recording of arrival times of laser pulses at the satellite. T2L2 was designed to achieve time stability in the range of 1 ps over 1000 s and an accuracy better than 100 ps. The project is in operation onboard the Jason-2 satellite since June 2008. The principle is based on the Satellite Laser Ranging (SLR) technology; it uses the input of 20–25 SLR stations of the international laser network which participate in the tracking. This paper focuses on the data reduction process which was developed specifically to transform the raw information given by both space instrument and ground network: first to identify the triplets (ground and onboard epochs and time of flight of the laser pulse), second to estimate a usable product in terms of ground-to-space time transfer (including instrumental corrections), and thirdly to produce synchronization between any pair of remote ground clocks. In describing the validation of time synchronizations, the paper opens a way for monitoring the time difference between ultra-stable clocks thanks to a laser link at a few ps level for Common View passes. It highlights however that without accurately characterizing the onboard oscillator of Jason-2 and knowing the unavailability of time calibrations of SLR stations generally, time transfer over intercontinental distances remain difficult to be accurately estimated.     

Precise orbit determination using the batch filter based on particle filtering with genetic resampling approach

In this study, genetic resampling (GRS) approach is utilized for precise orbit determination (POD) using the batch filter based on particle filtering (PF). Two genetic operations, which are arithmetic crossover and residual mutation, are used for GRS of the batch filter based on PF (PF batch filter). For POD, Laser-ranging Precise Orbit Determination System (LPODS) and satellite laser ranging (SLR) observations of the CHAMP satellite are used. Monte Carlo trials for POD are performed by one hundred times. The characteristics of the POD results by PF batch filter with GRS are compared with those of a PF batch filter with minimum residual resampling (MRRS). The post-fit residual, 3D error by external orbit comparison, and POD repeatability are analyzed for orbit quality assessments. The POD results are externally checked by NASA JPL’s orbits using totally different software, measurements, and techniques. For post-fit residuals and 3D errors, both MRRS and GRS give accurate estimation results whose mean root mean square (RMS) values are at a level of 5 cm and 10–13 cm, respectively. The mean radial orbit errors of both methods are at a level of 5 cm. For POD repeatability represented as the standard deviations of post-fit residuals and 3D errors by repetitive PODs, however, GRS yields 25% and 13% more robust estimation results than MRRS for post-fit residual and 3D error, respectively. This study shows that PF batch filter with GRS approach using genetic operations is superior to PF batch filter with MRRS in terms of robustness in POD with SLR observations.     

Design of tracking mount and controller for mobile satellite laser ranging system

In this study, we have proposed and implemented a design for the tracking mount and controller of the ARGO-M (Accurate Ranging system for Geodetic Observation - Mobile) which is a mobile satellite laser ranging (SLR) system developed by the Korea Astronomy and Space Science Institute (KASI) and Korea Institute of Machinery and Materials (KIMM). The tracking mount comprises a few core components such as bearings, driving motors and encoders. These components were selected as per the technical specifications for the tracking mount of the ARGO-M. A three-dimensional model of the tracking mount was designed. The frequency analysis of the model predicted that the first natural frequency of the designed tracking mount was high enough. The tracking controller is simulated using MATLAB/xPC Target to achieve the required pointing and tracking accuracy. In order to evaluate the system repeatability and tracking accuracy of the tracking mount, a prototype of the ARGO-M was fabricated, and repeatability tests were carried out using a laser interferometer. Tracking tests were conducted using the trajectories of low earth orbit (LEO) and high earth orbit (HEO) satellites. Based on the test results, it was confirmed that the prototype of the tracking mount and controller of the ARGO-M could achieve the required repeatability along with a tracking accuracy of less than 1 arcsec.     

Quality assessment of sub-Nyquist recovery from future gravity satellite missions

Drawing on experience from Gravity Recovery and Climate Experiment (GRACE) data analysis, the scientific challenges were already identified in several studies. Any future mission should focus on improvement in both precision and resolution in space and time. For future gravity missions which use high quality sensors, aliasing of high frequency time-variable geophysical signals to the lower frequency signals is one of the most serious problems. The aliasing problem and the spatio-temporal resolution are mainly restricted by two sampling theorems describing the space-time sampling of satellite missions: (i) a Heisenberg-like uncertainty theorem which states that the product of spatial resolution and time resolution is constant, and (ii) the Colombo–Nyquist rule (CNR), which requires the number of satellite revolutions in a repeat period to be at least twice a given maximum spherical harmonic degree. The CNR holds under the assumption of equal ground-track spacing, and limits the spatial resolution of the gravity solution.     

Different space weather effects in anomalies of the high and low orbital satellites

Preliminary results of the EU INTAS Project 00810, which aims to improve the methods of safeguarding satellites in the Earth’s magnetosphere from the negative effects of the space environment, are presented. Anomaly data from the “Kosmos” series satellites in the period 1971–1999 are combined in one database, together with similar information on other spacecraft. This database contains, beyond the anomaly information, various characteristics of the space weather: geomagnetic activity indices (Ap, AE and Dst), fluxes and fluences of electrons and protons at different energies, high energy cosmic ray variations and other solar, interplanetary and solar wind data. A comparative analysis of the distribution of each of these parameters relative to satellite anomalies was carried out for the total number of anomalies (about 6000 events), and separately for high (5000 events) and low (about 800 events) altitude orbit satellites. No relation was found between low and high altitude satellite anomalies. Daily numbers of satellite anomalies, averaged by a superposed epoch method around sudden storm commencements and proton event onsets for high (>1500 km) and low (<1500 km) altitude orbits revealed a big difference in a behavior. Satellites were divided on several groups according to the orbital characteristics (altitude and inclination). The relation of satellite anomalies to the environmental parameters was found to be different for various orbits that should be taken into account under developing of the anomaly frequency models.     

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

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

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.     

星地激光时间比对在卫星导航系统中的应用

星地精密激光时间比对具有精度高、系统误差少等优点,因此利用星地精密激光时间比对,不但可以对无线电时间比对进行外部精度检验,而且可以检验并分离无线电伪距测量的系统误差,分析检测设备时延的不稳定性和卫星钟的中短期性能指标,提高卫星钟差的预报精度。     

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.     

改进S／DB提高FY—2A云图配准精度的方法

针对ＦＹ－２Ａ卫星因扫描辐射计步进引起卫星自旋周期变化从而造成图像精度恶化问题,根据扫描辐射计的运动规律,提出地面同步与数据缓冲系统（Ｓ／ＤＢ）预先补偿系统误差的解决方法;增强了Ｓ／ＤＢ的跟踪水平,使卫星图像的配准精度提高了近三倍,达到了内外先进卫星的水平,实际应用效果良好。     

Contribution of satellite laser ranging to combined gravity field models

In the framework of satellite-only gravity field modeling, satellite laser ranging (SLR) data is typically exploited to recover long-wavelength features. This contribution provides a detailed discussion of the SLR component of GOCO02S, the latest release of combined models within the GOCO series. Over a period of five years (January 2006 to December 2010), observations to LAGEOS-1, LAGEOS-2, Ajisai, Stella, and Starlette were analyzed. We conducted a series of closed-loop simulations and found that estimating monthly sets of spherical harmonic coefficients beyond degree five leads to exceedingly ill-posed normal equation systems. Therefore, we adopted degree five as the spectral resolution for real data analysis. We compared our monthly coefficient estimates of degree two with SLR and Gravity Recovery and Climate Experiment (GRACE) time series provided by the Center for Space Research (CSR) at Austin, Texas. Significant deviations in C₂₀ were noted between SLR and GRACE; the agreement is better for the non-zonal coefficients. Fitting sinusoids together with a linear trend to our C₂₀ time series yielded a rate of (−1.75 ± 0.6) × 10⁻¹¹/yr; this drift is equivalent to a geoid change from pole to equator of 0.35 ± 0.12 mm/yr or an apparent Greenland mass loss of 178.5 ± 61.2 km³/yr. The mean of all monthly solutions, averaged over the five-year period, served as input for the satellite-only model GOCO02S. The contribution of SLR to the combined gravity field model is highest for C₂₀, and hence is essential for the determination of the Earth’s oblateness.