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     

First use of IGEX precise ephemerides for intercontinental GLONASS P-code time transfer

 Until recently, the Global Positioning System (GPS) was the only operational means of distributing time to an arbitrary number of users and of synchronizing clocks over large distances with a high degree of precision and accuracy. Over the last few years it has been shown that similar performance can be achieved using the Russian Global Navigation Satellite System (GLONASS). GLONASS time transfer between continents was initially hampered by the lack of post-processed precise ephemerides. Results from the International GLONASS Experiment (IGEX) campaign are now available, however, and this paper reports on the first use of IGEX precise ephemerides for GLONASS P-code intercontinental time links. The results of GLONASS P-code and GPS C/A-code time transfer are compared under similar conditions. Received: 31 January 2000 / Accepted: 10 July 2000     

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     

Integrity Monitoring of IGEX-98 Data,Part I: Availability

The International GLONASS Experiment 1998 (IGEX-98) was the first international tracking campaign of the Russian counterpart to the Global Positioning System (GPS), GLONASS. Started in October 1998, the campaign was originally scheduled to last for three months. However, the launch of additional GLONASS satellites and a widespread enthusiasm among the participants led to an indefinite continuation of the campaign on a “best effort” basis. At the Delft University of Technology, the data of six IGEX-98 stations have been analyzed in detail with integrity monitoring software, developed at the Department of Mathematical Geodesy and Positioning of the University. The software aims to detect outliers and slips in code and phase observations in real time. In addition, the software also allows the validation of the information contained in the broadcast navigation messages. The results of the IGEX-98 data analyses will be presented in a three-part series. In the second part, GLONASS outlier and slips statistics will be discussed, while in the third part the anomaly detection results of the GLONASS and GPS messages will be shown. In this first part of the series, however, the most basic of all statistics will be considered: a simple day-to-day count of the number of GLONASS and GPS observations. Although simple, this statistic yields a surprising amount of information both on the availability of the GLONASS satellites and on the peculiarities of some of the receiver makes participating in the IGEX-98 campaign. ? 2000 John Wiley & Sons, Inc.     

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     

IGEX-98: International GLONASS Experiment

GLONASS is a Russian military navigation satellite system. Even if the constellation is far from being operational, a large community is potentially interested in using its signal in combination with GPS to improve availability, integrity, or precision. The goal of this article is to present the international GLONASS campaign, named IGEX-98, which is really the first attempt to obtain precise GLONASS orbits for geodetic applications. A world-wide network of GLONASS receivers has been deployed and is still operational. Several analysis centers process the GLONASS data on a regular basis in the same way it is done for GPS within the International GPS Service. ? 1999 John Wiley & Sons, Inc.     

Integrity Monitoring of IGEX-98 Data,Part II: Cycle Slip and Outlier Detection

At the Delft University of Technology (DUT), the data of six stations participating in the first international GLONASS (Global Navigation Satellite System, the Russian counterpart of GPS) tracking campaign, IGEX-98, were analyzed with integrity monitoring software. The software was developed at the Department of Mathematical Geodesy and Positioning of DUT. The main function of this software is to detect slips and outliers in phase and code observations in real time. In addition, the software also allowe the validation of the information contained in the broadcast navigation messages. The results of the IGEX-98 data analyses will be presented in a three-part series. In this second part of the series, GLONASS outlier and slip statistics will be discussed in detail. The first part was concerned with the availability of GLONASS observations, while in the third one the broadcast navigation message validation results will be considered. ? 2000 John Wiley & Sons, Inc.     

GPS-assisted GLONASS orbit determination

 Using 1 week of data from a network of GPS/GLONASS dual-tracking receivers, 15-cm accurate GLONASS orbit determination is demonstrated with an approach that combines GPS and GLONASS data. GPS data are used to define the reference frame, synchronize receiver clocks and determine troposphere delay for the GLONASS tracking network. GLONASS tracking data are then processed separately, with the GPS-defined parameters held fixed, to determine the GLONASS orbit. The quality of the GLONASS orbit determination is currently limited by the size and distribution of the tracking network, and by the unavailability of a sufficiently refined solar pressure model. Temporal variations in the differential clock bias of the dual-tracking receivers are found to have secondary impact on the orbit determination accuracy. Received: 5 January 2000 / Accepted: 15 February 2001     

Integrity Monitoring of IGEX-98 Data – Part III: Broadcast Navigation Message Validation

At the Delft University of Technology (DUT), data of six stations participating in the first international GLONASS tracking campaign, IGEX-98, were analyzed with integrity monitoring software. The software was developed at the Department of Mathematical Geodesy and Positioning of DUT. The main function of the software is to detect slips and outliers in phase and code observations in real time. In addition, the software also allows the validation of the information contained in the broadcast navigation messages. The results of the IGEX-98 data analyses are presented in a three part series. In this third and final part of the series, GLONASS and GPS navigation message validation results will be discussed in detail. The first part was concerned with the availability of GLONASS observations (Jonkman & de Jong, 2000a), while in the second part, GLONASS slip and outlier statistics were considered (Jonkman & de Jong, 2000b). The discussion of the navigation message validation results concentrates on anomalies detected in the data of an IGEX-98 station in Switzerland. Data collected at this station from September 1998 to December 1999 was analyized. In all, 19 anomalies were detected by the integrity monitoring software, 8 in GLONASS messages and 11 in GPS messages. The cause of the anomalies was established by evaluating the original data files, studying official bulletins of the GLONASS and GPS ground segments, and in some cases tracing the anomaly through the complete IGEX-98 network. ? 2000 John Wiley & Sons, Inc.     

Analytical solar radiation pressure modelling for GLONASS using a pixel array

 A new method for calculating analytical solar radiation pressure models for GNSS spacecraft has been developed. The method simulates the flux of light from the Sun using a pixel array. The method can cope with a high level of complexity in the spacecraft structure and models effects due to reflected light. Models have been calculated and tested for the Russhar global navigation satellite system GLONASS IIv spacecraft. Results are presented using numerical integration of the force model and long-arc satellite laser ranging (SLR) analysis. The integrated trajectory differs from a precise orbit calculated using a network of global tracking stations by circa 2 m root mean square over a 160 000-km arc. The observed − computed residuals for the 400-day SLR arc are circa 28 mm. Received: 23 December 1999 / Accepted: 28 August 2000     

一个新的GLONASS与GPS的坐标转换关系

介绍了一组新的更精确的PZ-90与WGS-84两种坐标框架之间转换参数。其中详细描述了其转换参数的求解模型和实验方案，对求解的转换参数进行了分析。     

Direct transformation from geocentric coordinates to geodetic coordinates

 The transformation from geocentric coordinates to geodetic coordinates is usually carried out by iteration. A closed-form algebraic method is proposed, valid at any point on the globe and in space, including the poles, regardless of the value of the ellipsoid's eccentricity. Received: 14 August 2000 / Accepted: 26 June 2002     

Analysis of the first year of Earth rotation parameters with a sub-daily resolution gained at the CODE processing center of the IGS

 The solutions of the CODE Analysis Center submitted to the IGS, the International Global Position System (GPS) Service for Geodynamics, are based on three days of observation of about 80–100 stations of the IGS network. The Earth rotation parameters (ERPs) are assumed to vary linearly over the three days with respect to an a priori model. Continuity at the day boundaries as well as the continuity of the first derivatives are enforced by constraints. Since early April 1995 CODE has calculated a new ERP series with an increased time resolution of 2 hours. Again continuity is enforced at the 2-hours-interval boundaries. The analysis method is described, particularly how to deal with retrograde diurnal terms in the ERP series which may not be estimated with satellite geodetic methods. The results obtained from the first year of data covered by the time series (time interval from 4 April 1995 to 30 June 1996) are also discussed. The series is relatively homogeneous in the sense of the used orbit model and the a priori model for the ERPs. The largest source of excitation at daily and sub-daily periods is likely to be the effect of the ocean tides. There is good agreement between the present results and Topex/Poseidon ocean tide models, as well as with models based on Very Long Baseline Interferometry (VLBI) and Satellite Laser Ranging (SLR) data. Non-oceanic periodic variations are also observed in the series. Their origin is most probably a consequence of the GPS solution strategy; other possible sources are the atmospheric tides. Received: 13 July 1999 / Accepted: 21 March 2000     

Somigliana–Pizzetti gravity: the international gravity formula accurate to the sub-nanoGal level

 The Somigliana–Pizzetti gravity field (the International gravity formula), namely the gravity field of the level ellipsoid (the International Reference Ellipsoid), is derived to the sub-nanoGal accuracy level in order to fulfil the demands of modern gravimetry (absolute gravimeters, super conducting gravimeters, atomic gravimeters). Equations (53), (54) and (59) summarise Somigliana–Pizzetti gravity Γ(φ,u) as a function of Jacobi spheroidal latitude φ and height u to the order ?(10⁻¹⁰ Gal), and Γ(B,H) as a function of Gauss (surface normal) ellipsoidal latitude B and height H to the order ?(10⁻¹⁰ Gal) as determined by GPS (`global problem solver'). Within the test area of the state of Baden-Württemberg, Somigliana–Pizzetti gravity disturbances of an average of 25.452 mGal were produced. Computer programs for an operational application of the new international gravity formula with (L,B,H) or (λ,φ,u) coordinate inputs to a sub-nanoGal level of accuracy are available on the Internet. Received: 23 June 2000 / Accepted: 2 January 2001     

PZ-90与WGS-84的坐标转换参数

主要论述了GLONASS坐标系统PZ-90的定义及与GPS坐标系统WGS-84进行坐标转换的基本方法和数学模型。对现有转换参数的获取方法进行论述与评价，并对各种不同的转换参数进行了分析比较。     

现代大地测量参考系统

概述现代大地测量参考系统的定义及不同参考系统之间的关系.主要讨论我国当代大地测量界常使用的3种用以表示几何位置的参考系统1980年国家大地坐标系、全球大地测量系统 1984 (WGS 84)、国际地球参考系统(ITRS),和一种用以表示物理位置,高程的参考系统1985年国家高程基准.并讨论大地测量中框架和基准的概念.     

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     

An approach to GLONASS ambiguity resolution

 When processing global navigation satellite system (GLONASS) carrier phases, the standard double-differencing (DD) procedure cannot cancel receiver clock terms in the DD phase measurement equations due to the multiple frequencies of the carrier phases. Consequently, a receiver clock parameter has to be set up in the measurement equations in addition to baseline components and DD ambiguities. The resulting normal matrix unfortunately becomes singular. Methods to deal with this problem have been proposed in the literature. However, these methods rely on the use of pseudo-ranges. As pseudo-ranges are contaminated by multi-path and hardware delays, biases in these pseudo-ranges are significant, which may result in unreliable ambiguity resolution. A new approach is addressed that is not sensitive to the biases in the pseudo-ranges. The proposed approach includes such steps as converting the carrier phases to their distances to cancel the receiver clock errors, and searching for the most likely single-differenced (SD) ambiguity. Based on the results from the theoretical investigation, a practical procedure for GLONASS ambiguity resolution is presented. The initial experimental results demonstrate that the proposed approach is useable in cases of GLONASS and combined global positioning system (GPS) and GLONASS positioning. Received: 19 August 1998 / Accepted: 12 November 1999     

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/GLONASS定位系统融合的坐标转换研究

GPS和GLONASS卫星定位系统分别采用WGS-84和PZ-90坐标系。为统一两种卫星定位成果,欧、美、俄于20世纪90年代各自求出两坐标系之间的转换参数。目前三种参数尚未统一,对GPS/GLONASS联合定位造成较大影响。本文针对国外介绍PZ-90和WGS-84坐标系相互转换时常用的转换模型及三种不同转换参数进行比较分析。分别选用地面GPS参考站坐标和GLONASS卫星轨道坐标,用三种坐标转换参数进行转换,对转换结果差异及其对单点定位和相对定位精度造成的影响进行全面分析比较,得出一些有益结论。