Combined processing of the IGS and the IGEX network

 Since the beginning of the International Global Navigation Satellite System (GLONASS) Experiment, IGEX, in October 1998, the Center for Orbit Determination in Europe (CODE) has acted as an analysis center providing precise GLONASS orbits on a regular basis. In CODE's IGEX routine analysis the Global Positioning System (GPS) orbits and Earth rotation parameters are introduced as known quantities into the GLONASS processing. A new approach is studied, where data from the IGEX network are combined with GPS observations from the International GPS Service (IGS) network and all parameters (GPS and GLONASS orbits, Earth rotation parameters, and site coordinates) are estimated in one processing step. The influence of different solar radiation pressure parameterizations on the GLONASS orbits is studied using different parameter subsets of the extended CODE orbit model. Parameterization with three constant terms in the three orthogonal directions, D, Y, and X (D = direction satellite–Sun, Y = direction of the satellite's solar panel axis), and two periodic terms in the X-direction, proves to be adequate for GLONASS satellites. As a result of the processing it is found that the solar radiation pressure effect for the GLONASS satellites is significantly different in the Y-direction from that for the GPS satellites, and an extensive analysis is carried out to investigate the effect in detail. SLR observations from the ILRS network are used as an independent check on the quality of the GLONASS orbital solutions. Both processing aspects, combining the two networks and changing the orbit parameterization, significantly improve the quality of the determined GLONASS orbits compared to the orbits stemming from CODE's IGEX routine processing. Received: 10 May 2000 / Accepted: 9 October 2000     

CODE’s new ultra-rapid orbit and ERP products for the IGS

The International GNSS Service (IGS) issues four sets of so-called ultra-rapid products per day, which are based on the contributions of the IGS Analysis Centers. The traditional (“old”) ultra-rapid orbit and earth rotation parameters (ERP) solution of the Center for Orbit Determination in Europe (CODE) was based on the output of three consecutive 3-day long-arc rapid solutions. Information from the IERS Bulletin A was required to generate the predicted part of the old CODE ultra-rapid product. The current (“new”) product, activated in November 2013, is based on the output of exactly one multi-day solution. A priori information from the IERS Bulletin A is no longer required for generating and predicting the orbits and ERPs. This article discusses the transition from the old to the new CODE ultra-rapid orbit and ERP products and the associated improvement in reliability and performance. All solutions used in this article were generated with the development version of the Bernese GNSS Software. The package was slightly extended to meet the needs of the new CODE ultra-rapid generation.     

GPS measurements of ocean loading and its impact on zenith tropospheric delay estimates: a case study in Brittany, France

 The results from a global positioning system (GPS) experiment carried out in Brittany, France, in October 1999, aimed at measuring crustal displacements caused by ocean loading and quantifying their effects on GPS-derived tropospheric delay estimates, are presented. The loading effect in the vertical and horizontal position time series is identified, however with significant disagreement in amplitude compared to ocean loading model predictions. It is shown that these amplitude misfits result from spatial tropospheric heterogeneities not accounted for in the data processing. The effect of ocean loading on GPS-derived zenith total delay (ZTD) estimates is investigated and a scaling factor of 4.4 between ZTD and station height for a 10° elevation cut-off angle is found (i.e. a 4.4-cm station height error would map into a 1-cm ZTD error). Consequently, unmodeled ocean loading effects map into significant errors in ZTD estimates and ocean loading modeling must be properly implemented when estimating ZTD parameters from GPS data for meteorological applications. Ocean loading effects must be known with an accuracy of better than 3 cm in order to meet the accuracy requirements of meteorological and climatological applications of GPS-derived precipitable water vapor. Received: 16 July 2001 / Accepted: 25 April 2002 Acknowledgments. The authors are grateful to H.G. Scherneck for fruitful discussions and for his help with the ocean loading calculations. They thank H. Vedel for making the HIRLAM data available; D. Jerett for helpful discussions; and the city of Rostrenen, the Laboratoire d'Océanographie of Concarneau, and the Institut de Protection et de S?reté Nucléaire (BERSSIN) for their support during the GPS measurement campaign. Reviews by C.K. Shum and two anonymous referees significantly improved this paper. This work was carried out in the framework of the MAGIC project (http://www.acri.fr/magic), funded by the European Commission, Environment and Climate Program (EC Contract ENV4-CT98–0745). Correspondence to: E. Calais, Department of Earth and Atmospheric Sciences, Purdue University, West Lafayette, IN 47907-1397, USA. e-mail: ecalais@purdue.edu Tel. : +1-765-496-2915; Fax:+1-765-496-1210     

The international GLONASS experiment: products, progress and prospects

 In October 1998 the IGEX field campaign, the first coordinated international effort to monitor GLONASS satellites on global basis, was started. Currently about 40 institutions worldwide support this effort either by providing GLONASS tracking data or in operating related data and analysis centers. The increasing quality and consistency of the calculated GLONASS orbits (about 25 cm early in 2000), even after the end of the official IGEX field campaign, are shown. Particular attention is drawn to the combination of precise ephemerides in order to generate a robust, reliable and complete IGEX orbits product. Some problems in modeling the effect of solar radiation pressure on GLONASS satellites are demonstrated. Finally, the expected benefits and prospects of the upcoming International GLOnass Service-Pilot Project (IGLOS-PP) of the International GPS Service (IGS) are discussed in more detail. Received: 17 August 2000 / Accepted: 12 April 2001     

Prediction of Earth orientation parameters by artificial neural networks

 Earth orientation parameters (EOPs) [polar motion and length of day (LOD), or UT1–UTC] were predicted by artificial neural networks. EOP series from various sources, e.g. the C04 series from the International Earth Rotation Service and the re-analysis optical astrometry series based on the HIPPARCOS frame, served for training the neural network for both short-term and long-term predictions. At first, all effects which can be described by functional models, e.g. effects of the solid Earth tides and the ocean tides or seasonal atmospheric variations of the EOPs, were removed. Only the differences between the modeled and the observed EOPs, i.e. the quasi-periodic and irregular variations, were used for training and prediction. The Stuttgart neural network simulator, which is a very powerful software tool developed at the University of Stuttgart, was applied to construct and to validate different types of neural networks in order to find the optimal topology of the net, the most economical learning algorithm and the best procedure to feed the net with data patterns. The results of the prediction were analyzed and compared with those obtained by other methods. The accuracy of the prediction is equal to or even better than that by other prediction methods. Received: 6 February 2001 / Accepted: 23 October 2001     

Gravitational perturbation theory for intersatellite tracking

 An analytical gravitational perturbation theory for the intersatellite tracking range and range-rate measurement between two satellites is developed. The satellite-to-satellite tracking (SST) range data measure the difference between the position perturbations of two satellites along the direction of the intersatellite range. The SST range-rate data measure the difference between the velocity perturbations along the direction of the intersatellite range, and the difference of the position perturbation along the direction perpendicular to the intersatellite range (cross-range). The SST range and range rate depend on different orbital excitations for mapping the gravity field. For the Gravity Recovery and Climate Experiment (GRACE), approximately 97% of the geopotential coefficient pairs produce perturbations with a root-mean-square larger than 1 m on the range and 0.1 m/sec on the range rate based on the EGM96 gravity field truncated at degree and order 140. Results in this study showed that ocean tides produce significant perturbations in the range and range-rate measurements. An ocean tide field with a higher degree and order (>70) is required to model the ocean tide perturbations on the intersatellite range and range-rate measurement. Received: 17 May 2000 / Accepted: 3 September 2001     

Auto-covariance estimation of variable samples (ACEVS) and its application for monitoring random ionospheric disturbances using GPS

 Ionospheric variation may be considered as a stationary time series under quiet conditions. However, the disturbance of a stationary random process from stationarity results in the bias of corresponding samples from the stationary observations, and in the change of statistical model parameters of the process. From a general mathematical aspect, a new method is presented for monitoring ionospheric variations, based on the characteristic of time-series observation of GPS, and an investigation of the statistical properties of the estimated auto-covariance of the random ionospheric delay when changing the number of samples in the time series is carried out. A preliminary scheme for monitoring ionospheric delays is proposed. Received: 18 August 2000 / Accepted: 12 April 2001     

Time transfer using GPS carrier phase: error propagation and results

 A joint time-transfer project between the Astronomical Institute of the University of Berne (AIUB) and the Swiss Federal Office of Metrology and Accreditation (METAS) was initiated to investigate the power of the time transfer using GPS carrier phase observations. Studies carried out in the context of this project are presented. The error propagation for the time-transfer solution using GPS carrier phase observations was investigated. To this purpose a simulation study was performed. Special interest was focussed on errors in the vertical component of the station position, antenna phase-center variations and orbit errors. A constant error in the vertical component introduces a drift in the time-transfer results for long baselines in east–west directions. The simulation study was completed by investigating the profit for time transfer when introducing the integer carrier phase ambiguities from a double-difference solution. This may reduce the drift in the time-transfer results caused by constant vertical error sources. The results from the present time-transfer solution are shown in comparison to results obtained with independent time-transfer techniques. The interpretation of the comparison benefits from the investigations of the error propagation study. Two types of solutions are produced on a regular basis at AIUB: one based on the rapid orbits from CODE, the other on the CODE final orbits. The rapid solution is available the day after the observations and has nearly the same quality as the final solution, which has a latency of about one week. The differences between these two solutions are below the nanosecond level. The differences from independent time-transfer techniques such as TWSTFT (two-way satellite time and frequency transfer) are a few nanoseconds for both products. Received: 15 November 2001 / Accepted: 6 September 2002 Correspondence to:R. Dach     

IGS Earth Rotation Parameters

Since its official start in January 1994, the International GPS Service (IGS) has been distributing, as part of its product combination, two distinct Earth rotation parameter (ERP) series: the IGS Rapid series and the IGS Final series. Initially, the IGS Rapid ERP values were interpolations of the International Earth Rotation Service (IERS) Bulletin A, whereas the IGS Final ERP series was based on the IERS Bulletin B. Since June 1996, the IGS has been generating its own Final ERP series consistent with the IGS combined orbit products and based on weighted means of individual IGS analysis center (AC) solutions. At first, only the polar motion (PM) coordinates and their rates were combined. Length of Day (LOD) and Universal Time (UT) solutions, also based on separate weighted mean combinations, followed in March 1997. Currently, the IGS Rapid and Final combinations are produced and made available within 17 hours and 11 days, respectively, after the last observation. Both IGS and the best AC series are consistent and precise at the 0.1-milliarcsecond (mas) level for PM and at about 30 μs for LOD. Biases in some AC solutions may exceed these consistency levels. Comparisons of both IGS ERP series with external standards, such as the IERS multitechnique Bulletins and atmospheric angular momentum series, confirm the estimated precisions. ? 1999 John Wiley & Sons, Inc.     

Combination procedure for length-of-day time series according to the noise frequency behavior

 Length-of-day (LOD) estimates from seven Global Positioning System (GPS) and three satellite laser ranging (SLR) analysis centers were combined into an even-spaced time series for a 27-month period (1996–1998). This time series was compared to the multi-technique Earth-orientation-parameter (EOP) combined solution (C04) derived at the Central Bureau of the International Earth Rotation Service (IERS/CB). Due to inhomogeneities in the different series derived from the various techniques (time, length, quality, and spatial resolution), the concept of a combined solution is justified. The noise behavior in LOD for different techniques varies with frequency; the data series were divided into frequency windows after removing both biases and trends. Different weight factors were assigned in each window, discriminating by technique, and produced one-technique combined solutions. Finally, these one-technique combined solutions were combined to obtain the final multi-technique solution. The LOD combined time series obtained by the present method based on the frequency windows combined series (FWCS) is very close to the IERS C04 solution. It could be useful to generate a new LOD reference time series to be used in the study of high-frequency variations of Earth rotation. Received: 28 March 2000 / Accepted: 15 February 2001     

海潮改正模型对GPS基线的影响

利用GAMIT软件解算了2011年前200天国内9个IGS站的数据,解算时考虑采用未加入和分别加入4种不同的海潮改正模型等情况,对4种不同海潮改正模型在基线分量上的影响进行了比较,说明了不同模型对GPS基线的影响差异较小;并对含沿海地区GPS测站和只含内陆地区GPS测站的基线进行了分析,结果表明含沿海地区GPS测站的基线比只含内陆地区GPS测站的基线受海潮改正模型的影响大;最后重点采用功率谱分析了海潮改正对GPS基线分量影响的周期,得到其3个方向的变化周期。     

Impact of Earth radiation pressure on GPS position estimates

GPS satellite orbits available from the International GNSS Service (IGS) show a consistent radial bias of up to several cm and a particular pattern in the Satellite Laser Ranging (SLR) residuals, which are suggested to be related to radiation pressure mismodeling. In addition, orbit-related frequencies were identified in geodetic time series such as apparent geocenter motion and station displacements derived from GPS tracking data. A potential solution to these discrepancies is the inclusion of Earth radiation pressure (visible and infrared) modeling in the orbit determination process. This is currently not yet considered by all analysis centers contributing to the IGS final orbits. The acceleration, accounting for Earth radiation and satellite models, is introduced in this paper in the computation of a global GPS network (around 200 IGS sites) adopting the analysis strategies from the Center for Orbit Determination in Europe (CODE). Two solutions covering 9 years (2000–2008) with and without Earth radiation pressure were computed and form the basis for this study. In previous studies, it has been shown that Earth radiation pressure has a non-negligible effect on the GPS orbits, mainly in the radial component. In this paper, the effect on the along-track and cross-track components is studied in more detail. Also in this paper, it is shown that Earth radiation pressure leads to a change in the estimates of GPS ground station positions, which is systematic over large regions of the Earth. This observed “deformation” of the Earth is towards North–South and with large scale patterns that repeat six times per GPS draconitic year (350 days), reaching a magnitude of up to 1 mm. The impact of Earth radiation pressure on the geocenter and length of day estimates was also investigated, but the effect is found to be less significant as compared to the orbits and position estimates.     

Local geoid determination combining gravity disturbances and GPS/levelling: a case study in the Lake Nasser area, Aswan, Egypt

 The use of GPS for height control in an area with existing levelling data requires the determination of a local geoid and the bias between the local levelling datum and the one implicitly defined when computing the local geoid. If only scarse gravity data are available, the heights of new data may be collected rapidly by determining the ellipsoidal height by GPS and not using orthometric heights. Hence the geoid determination has to be based on gravity disturbances contingently combined with gravity anomalies. Furthermore, existing GPS/levelling data may also be used in the geoid determination if a suitable general gravity field modelling method (such as least-squares collocation, LSC) is applied. A comparison has been made in the Aswan Dam area between geoids determined using fast Fourier transform (FFT) with gravity disturbances exclusively and LSC using only the gravity disturbances and the disturbances combined with GPS/levelling data. The EGM96 spherical harmonic model was in all cases used in a remove–restore mode. A total of 198 gravity disturbances spaced approximately 3 km apart were used, as well as 35 GPS/levelling points in the vicinity and on the Aswan Dam. No data on the Nasser Lake were available. This gave difficulties when using FFT, which requires the use of gridded data. When using exclusively the gravity disturbances, the agreement between the GPS/levelling data were 0.71 ± 0.17 m for FFT and 0.63 ± 0.15 for LSC. When combining gravity disturbances and GPS/levelling, the LSC error estimate was ±0.10 m. In the latter case two bias parameters had to be introduced to account for a possible levelling datum difference between the levelling on the dam and that on the adjacent roads. Received: 14 August 2000 / Accepted: 28 February 2001     

Apropos laser tracking to GPS satellites

. Laser tracking to GPS satellites (PRN5 and 6) provides an opportunity to compare GPS and laser systems directly and to combine data of both in a single solution. A few examples of this are given in this study. The most important results of the analysis are that (1) daily SLR station coordinate solutions could be generated with a few cm accuracy; (2) coordinates of nine stations were determined in a 2.3-year-long arc solution; (3) the contribution of laser data on the `SLR-GPS' combined orbit, resulting from the simultaneous processing of SLR and GPS data, is significant and (4) laser-only orbits have an accuracy of 10–20 cm, 1-day predictions of SLR orbits differ from IGS orbits by about 20–40 cm, 2-day predictions by 50–60 cm. Received: 1 October 1996 / Accepted: 14 February 1997     

Tectonic plate motion and co-seismic steps surveyed by DORIS field beacons: the New Hebrides experiment

 The New Hebrides experiment consisted of setting up a pair of DORIS beacons in remote tropical islands in the southwestern Pacific, between 1993 and 1997. Because of orbitography requirements on TOPEX/Poséidon, the beacons were only transmitting to SPOT satellites. Root-mean-square (RMS) scatters at the centimeter level on the latitude and vertical components were achieved, but 2-cm RMS scatters affected the longitude component. Nevertheless, results of relative velocity (123 mm/year N250°) are very consistent with those obtained using the global positioning system (GPS) (126 mm/yr N246°). The co-seismic step (12 mm N60°) related to the Walpole event (M _W = 7.7) is consistent with that derived from GPS (10 mm N30°) or from the centroid moment tensor (CMT) of the quake (12 mm N000°). Received: 19 November 1999 / Accepted: 17 May 2000     

GPS network monitors the Western Alps' deformation over a five-year period: 1993–1998

  The Western Alps are among the best studied collisional belts with both detailed structural mapping and also crustal geophysical investigations such as the ECORS and EGT seismic profile. By contrast, the present-day kinematics of the belt is still largely unknown due to small relative motions and the insufficient accuracy of the triangulation data. As a consequence, several tectonic problems still remain to be solved, such as the amount of N–S convergence in the Occidental Alps, the repartition of the deformation between the Alpine tectonic units, and the relation between deformation and rotation across the Alpine arc. In order to address these problems, the GPS ALPES group, made up of French, Swiss and Italian research organizations, has achieved the first large-scale GPS surveys of the Western Alps. More than 60 sites were surveyed in 1993 and 1998 with a minimum observation of 3 days at each site. GPS data processing has been done by three independent teams using different software. The different solutions have horizontal repeatabilities (N–E) of 4–7 mm in 1993 and 2–3 mm in 1998 and compare at the 3–5-mm level in position and 2-mm/yr level in velocity. A comparison of 1993 and 1998 coordinates shows that residual velocities of the GPS marks are generally smaller than 2 mm/yr, precluding a detailed tectonic interpretation of the differential motions. However, these data seem to suggest that the N–S compression of the Western Alps is quite mild (less than 2 mm/yr) compared to the global convergence between the African and Eurasian plate (6 mm/yr). This implies that the shortening must be accomodated elsewhere by the deformation of the Maghrebids and/or by rotations of Mediterranean microplates. Also, E–W velocity components analysis supports the idea that E–W extension exists, as already suggested by recent structural and seismotectonic data interpretation. Received: 27 November 2000 / Accepted: 17 September 2001     

Random simulation and GPS decorrelation

 (i) A random simulation approach is proposed, which is at the centre of a numerical comparison of the performances of different GPS decorrelation methods. The most significant advantage of the approach is that it does not depend on nor favour any particular satellite–receiver geometry and weighting system. (ii) An inverse integer Cholesky decorrelation method is proposed, which will be shown to out-perform the integer Gaussian decorrelation and the Lenstra, Lenstra and Lovász (LLL) algorithm, and thus indicates that the integer Gaussian decorrelation is not the best decorrelation technique and that further improvement is possible. (iii) The performance study of the LLL algorithm is the first of its kind and the results have shown that the algorithm can indeed be used for decorrelation, but that it performs worse than the integer Gaussian decorrelation and the inverse integer Cholesky decorrelation. (iv) Simulations have also shown that no decorrelation techniques available to date can guarantee a smaller condition number, especially in the case of high dimension, although reducing the condition number is the goal of decorrelation. Received: 26 April 2000 / Accepted: 5 March 2001     

The activities of the Ionosphere Working Group of the International GPS Service (IGS)

This article is based on a position paper presented at the IGS Network, Data and Analysis Center Workshop 2002 in Ottawa, Canada, 8–11 April 2002, and introduces the IGS Ionosphere Working Group (Iono_WG). Detailed information about the IGS in general can be found on the IGS Central Bureau Web page: http://igscb.jpl.nasa.gov. The Iono_WG commenced working in June 1998. The working group's main activity currently is the routine production of ionosphere Total Electron Content (TEC) maps with a 2-h time resolution and daily sets of GPS satellite and receiver hardware differential code bias (DCB) values. The TEC maps and DCB sets are derived from GPS dual-frequency tracking data recorded with the global IGS tracking network. In the medium- and long-term, the working group intends to refine algorithms for the mapping of ionospheric parameters from GPS measurements and to realize near–real–time availability of IGS ionosphere products. The paper will give an overview of the Iono_WG activities that include a summary of activities since its establishment, achievements and future plans. Electronic Publication     

Stochastic assessment of GPS carrier phase measurements for precise static relative positioning

 Global positioning system (GPS) carrier phase measurements are used in all precise static relative positioning applications. The GPS carrier phase measurements are generally processed using the least-squares method, for which both functional and stochastic models need to be carefully defined. Whilst the functional model for precise GPS positioning is well documented in the literature, realistic stochastic modelling for the GPS carrier phase measurements is still both a controversial topic and a difficult task to accomplish in practice. The common practice of assuming that the raw GPS measurements are statistically independent in space and time, and have the same accuracy, is certainly not realistic. Any mis-specification in the stochastic model will inevitably lead to unreliable positioning results. A stochastic assessment procedure has been developed to take into account the heteroscedastic, space- and time-correlated error structure of the GPS measurements. Test results indicate that the reliability of the estimated positioning results is improved by applying the developed stochastic assessment procedure. In addition, the quality of ambiguity resolution can be more realistically evaluated. Received: 13 February 2001 / Accepted: 3 September 2001     

Laser-based validation of GLONASS orbits by short-arc technique

 The International GLONASS Experiment (IGEX-98) was carried out between 19 October 1998 and 19 April 1999. Among several objectives was the precise orbit determination of GPS and GLONASS satellites and its validation by laser ranging observations. Local laser-based orbit corrections (radial, tangential and normal components in a rotating orbital local reference frame) are computed using a geometrical short-arc technique. The order of magnitude of these corrections is at the level of few decimeters, depending on the considered components. The orbit corrections are analyzed as a function of several parameters (date, orbital plane, geographical area). The mean corrections are at the level of several centimeters. However, when averaging over the entire campaign and for all the satellites, no mean radial, tangential and normal orbit corrections are found. The origin of the observed corrections is considered (errors due to the geocentric gravitational constant, the non-gravitational forces, the thermal equilibrium of on-board equipment, the reference systems, the location and the signature of the retroreflector array, and the precision of the satellite laser ranges). Some features are also due to errors in the radio-tracking GLONASS orbits. Further investigations will be needed to better understand the origin of various biases. Received: 17 February 2000 / Accepted: 31 January 2001