A velocity field estimation of the Brazilian portion of the SOAM plate

With the proposition for the adoption of Geocentric Reference System for the Americas (SIRGAS) as a terrestrial reference frame for South America, the need for temporal monitoring of station coordinates used in its materialization has become apparent. This would provide a dynamic characterization of the frame. The Brazilian Network for Continuous Monitoring of GPS (RBMC) has collected high accuracy GPS measurements since 1996. The Brazilian Institute of Geography and Statistics (IBGE) maintains this network in collaboration with several universities and organizations. Most of the stations are also part of the SIRGAS network. The RBMC also contributes data to the International Terrestrial Reference System (ITRS) to densify the global frame. Two of the RBMC stations are also part of the International GPS Service (IGS). This paper reports initial results from these stations. To estimate the velocity field defined by these stations, ten IGS stations located on the border of the South American plate and in adjacent plates, along with nine RBMC stations, were used. Observations covering five groups of 15 days each were used. These groups of observations were at epochs 1997.3, 1997.9, 1998.3, 1998.9 and 1999.2. Seven IGS stations were chosen to have their coordinates constrained to those epochs. IGS products (precise ephemeris and clocks) were used to process the daily solutions, which were carried out with Bernese software. Carrier phase double differences were formed using the ionospheric-delay free observable. The troposphere was modeled using a combination of the Saastamoinen model and the Niell mapping function. A tropospheric parameter was estimated every two hours. The results of the daily baseline solutions were combined using the summation of normal equations technique, in which the final coordinates and velocities were estimated. The results were compared with various models, such as the NNR-NUVEL1 and the APKIM8.80. Velocity vectors estimated for the RBMC stations show good agreement with those two models, with rates approximately equal to 2 cm/year.     

IGS reference frames: status and future improvements

The hierarchy of reference frames used in the International GPS Service (IGS) and the procedures and rationale for realizing them are reviewed. The Conventions of the International Earth Rotation and Reference Systems Service (IERS) lag developments in the IGS in a number of important respects. Recommendations are offered for changes in the IERS Conventions to recognize geocenter motion (as already implemented by the IGS) and to enforce greater model consistency in order to achieve higher precision for combined reference frame products. Despite large improvements in the internal consistency of IGS product sets, defects remain which should be addressed in future developments. If the IGS is to remain a leader in this area, then a comprehensive, long-range strategy should be formulated and pursued to maintain and enhance the IGS reference frame, as well as to improve its delivery to users. Actions should include the official designation of a high-performance reference tracking network whose stations are expected to meet the highest standards possible.Also published in the proceedings of the workshop and symposium Celebrating a Decade of the International GPS Service, Astronomical Institute, University of Bern, Switzerland.     

IGS08: the IGS realization of ITRF2008

On April 17, 2011, the International GNSS Service (IGS) stopped using the IGS05 reference frame and adopted a new one, called IGS08, as the basis of its products. The latter was derived from the latest release of the International Terrestrial Reference Frame (ITRF2008). However, the simultaneous adoption of a new set of antenna phase center calibrations by the IGS required slight adaptations of ITRF2008 positions for 65 of the 232 IGS08 stations. The impact of the switch from IGS05 to IGS08 on GNSS station coordinates was twofold: in addition to a global transformation due to the frame change from ITRF2005 to ITRF2008, many station coordinates underwent small shifts due to antenna calibration updates, which need to be accounted for in any comparison or alignment of an IGS05-consistent solution to IGS08. Because the heterogeneous distribution of the IGS08 network makes it sub-optimal for the alignment of global frames, a smaller well-distributed sub-network was additionally designed and designated as the IGS08 core network. Only 2?months after their implementation, both the full IGS08 network and the IGS08 core network already strongly suffer from the loss of many reference stations. To avoid a future crisis situation, updates of IGS08 will certainly have to be considered before the next ITRF release.     

Geocenter variations derived from a combined processing of LEO- and ground-based GPS observations

GNSS observations provided by the global tracking network of the International GNSS Service (IGS, Dow et al. in J Geod 83(3):191–198, 2009) play an important role in the realization of a unique terrestrial reference frame that is accurate enough to allow a detailed monitoring of the Earth’s system. Combining these ground-based data with GPS observations tracked by high-quality dual-frequency receivers on-board low earth orbiters (LEOs) is a promising way to further improve the realization of the terrestrial reference frame and the estimation of geocenter coordinates, GPS satellite orbits and Earth rotation parameters. To assess the scope of the improvement on the geocenter coordinates, we processed a network of 53 globally distributed and stable IGS stations together with four LEOs (GRACE-A, GRACE-B, OSTM/Jason-2 and GOCE) over a time interval of 3 years (2010–2012). To ensure fully consistent solutions, the zero-difference phase observations of the ground stations and LEOs were processed in a common least-squares adjustment, estimating all the relevant parameters such as GPS and LEO orbits, station coordinates, Earth rotation parameters and geocenter motion. We present the significant impact of the individual LEO and a combination of all four LEOs on the geocenter coordinates. The formal errors are reduced by around 20% due to the inclusion of one LEO into the ground-only solution, while in a solution with four LEOs LEO-specific characteristics are significantly reduced. We compare the derived geocenter coordinates w.r.t. LAGEOS results and external solutions based on GPS and SLR data. We found good agreement in the amplitudes of all components; however, the phases in x- and z-direction do not agree well.     

The international global navigation satellite systems service (IGS): development and achievements

Since 21 June 1992 the International GPS Service (IGS), renamed International GNSS Service in 2005, produces and makes available uninterrupted time series of its products, in particular GPS observations from the IGS Global Network, GPS orbits, Earth orientation parameters (components x and y of polar motion, length of day) with daily time resolution, satellite and receiver clock information for each day with different latencies and accuracies, and station coordinates and velocities in weekly batches for further analysis by the IERS (International Earth Rotation and Reference Systems Service). At a later stage the IGS started exploiting its network for atmosphere monitoring, in particular for ionosphere mapping, for troposphere monitoring, and time and frequency transfer. This is why new IGS products encompass ionosphere maps and tropospheric zenith delays. This development became even more important when more and more space-missions carrying space-borne GPS for various purposes were launched. This article offers an overview for the broader scientific community of the development of the IGS and of the spectrum of topics addressed today with IGS data and products.     

Single epoch estimation of the Galileo integrity chain sensor station clock offsets

The Galileo integrity chain depends on a number of key factors, one of which is contamination of the signal-in-space errors with residual errors other than imperfect modelling of satellite orbits and clocks. A potential consequence of this is that the user protection limit is driven not by the errors associated with the imperfect orbit and clock modelling, but by the distortions induced by noise and bias in the integrity chain. These distortions increase the minimum bias the integrity chain can guarantee to detect, which is reflected in the user protection limit. A contributor to this distortion is the inaccuracy associated with the estimation of the offset between the Galileo sensor station (GSS) receiver clocks and the Galileo system time (GST). This offset is termed the receiver clock synchronization error (CSE). This paper describes the research carried out to determine both the CSE and its associated error using GPS data as captured with the Galileo System Test Bed Version 1 (GSTB-V1). In the study we simulate open access to a time datum using IGS data. Two methods are compared for determining CSE and the corresponding uncertainty (noise) across a global network of tracking stations. The single-epoch single-station method is an ‘averaging’ technique that uses a single epoch of data, and is carried out at individual sensor stations, without recourse to the data from other stations. The global network solution method is also single epoch based, but uses the inversion of a linearised model of the global system to solve for the CSE simultaneously at all GSS along with a number of other parameters that would otherwise be absorbed into the CSE estimate in the averaging technique. To test the effectiveness of various configurations in the two methods the estimated synchronisation errors across the GSS network (comprising 25 stations) are compared to the same values as estimated by the International GPS Service (IGS) using a global tracking network of around 150 stations, as well as precise orbit and satellite clock models determined by a combination of global analysis centres. The results show that the averaging technique is vulnerable to unmodelled errors in the satellite clock offsets from system time, leading to receiver CSE errors in the region of 12 ns (3.7 m), this value being largely driven by the satellite CSE errors. The global network approach is capable of delivering CSE errors at the level of 1.5 ns (46 cm) depending on the number of parameters in the linearised model. The International GNSS Service (IGS) receiver clock estimates were used as a truth model for comparative assessment.     

Implementation and testing of the gridded Vienna Mapping Function 1 (VMF1)

The new gridded Vienna Mapping Function (VMF1) was implemented and compared to the well-established site-dependent VMF1, directly and by using precise point positioning (PPP) with International GNSS Service (IGS) Final orbits/clocks for a 1.5-year GPS data set of 11 globally distributed IGS stations. The gridded VMF1 data can be interpolated for any location and for any time after 1994, whereas the site-dependent VMF1 data are only available at selected IGS stations and only after 2004. Both gridded and site-dependent VMF1 PPP solutions agree within 1 and 2 mm for the horizontal and vertical position components, respectively, provided that respective VMF1 hydrostatic zenith path delays (ZPD) are used for hydrostatic ZPD mapping to slant delays. The total ZPD of the gridded and site-dependent VMF1 data agree with PPP ZPD solutions with RMS of 1.5 and 1.8 cm, respectively. Such precise total ZPDs could provide useful initial a priori ZPD estimates for kinematic PPP and regional static GPS solutions. The hydrostatic ZPDs of the gridded VMF1 compare with the site-dependent VMF1 ZPDs with RMS of 0.3 cm, subject to some biases and discontinuities of up to 4 cm, which are likely due to different strategies used in the generation of the site-dependent VMF1 data. The precision of gridded hydrostatic ZPD should be sufficient for accurate a priori hydrostatic ZPD mapping in all precise GPS and very long baseline interferometry (VLBI) solutions. Conversely, precise and globally distributed geodetic solutions of total ZPDs, which need to be linked to VLBI to control biases and stability, should also provide a consistent and stable reference frame for long-term and state-of-the-art numerical weather modeling.     

The International GNSS Service in a changing landscape of Global Navigation Satellite Systems

The International GNSS Service (IGS) is an international activity involving more than 200 participating organisations in over 80 countries with a track record of one and a half decades of successful operations. The IGS is a service of the International Association of Geodesy (IAG). It primarily supports scientific research based on highly precise and accurate Earth observations using the technologies of Global Navigation Satellite Systems (GNSS), primarily the US Global Positioning System (GPS). The mission of the IGS is “to provide the highest-quality GNSS data and products in support of the terrestrial reference frame, Earth rotation, Earth observation and research, positioning, navigation and timing and other applications that benefit society”. The IGS will continue to support the IAG’s initiative to coordinate cross-technique global geodesy for the next decade, via the development of the Global Geodetic Observing System (GGOS), which focuses on the needs of global geodesy at the mm-level. IGS activities are fundamental to scientific disciplines related to climate, weather, sea level change, and space weather. The IGS also supports many other applications, including precise navigation, machine automation, and surveying and mapping. This article discusses the IGS Strategic Plan and future directions of the globally-coordinated ~400 station IGS network, tracking data and information products, and outlines the scope of a few of its numerous working groups and pilot projects as the world anticipates a truly multi-system GNSS in the coming decade.     

基于GAMIT/GLOBK和GPS观测确定加拿大地区的冰川均衡调整导致的垂向形变

本文采用加拿大地区31个IGS台站的观测数据,利用GAMIT/GLOBK软件处理了时间跨度为2000—2018年的GPS原始观测资料,得到了ITRF2014参考框架下的台站位置垂向运动速率的时间序列。对GPS时间序列进行阶跃探测及修复、异常值探测及剔除、趋势项估计、去除近期冰川融化导致的地表弹性变形后,得到了由冰川均衡调整(GIA)导致的GPS台站抬升速率。本文的结果与前人基于GPS观测得到的结果(在ITRF2008框架下)相差在2 mm/a以内,与ICE6G系列GIA模型预测值相差在3 mm/a以内,因此验证了本文结果的正确性和可靠性,为进一步利用全球GPS台站观测数据研究GIA垂直形变速率,进而约束和改进GIA模型打下了坚实基础。     

Plate Motion of India and Interseismic Strain in the Nepal Himalaya from GPS and DORIS Measurements

We analyse geodetically estimated deformation across the Nepal Himalaya in order to determine the geodetic rate of shortening between Southern Tibet and India, previously proposed to range from 12 to 21 mm yr^?1. The dataset includes spirit-levelling data along a road going from the Indian to the Tibetan border across Central Nepal, data from the DORIS station on Everest, which has been analysed since 1993, GPS campaign measurements from surveys carried on between 1995 and 2001, as well as data from continuous GPS stations along a transect at the logitude of Kathmandu operated continuously since 1997. The GPS data were processed in International Terrestrial Reference Frame 2000 (ITRF2000), together with the data from 20 International GNSS Service (IGS) stations and then combined using quasi- observation combination analysis (QOCA). Finally, spatially complementary velocities at stations in Southern Tibet, initially determined in ITRF97, were expressed in ITRF2000. After analysing previous studies by different authors, we determined the pole of rotation of the Indian tectonic plate to be located in ITRF2000 at 51.409±1.560° N and ?10.915±5.556°E, with an angular velocity of 0.483±0.015°. Myr^?1. Internal deformation of India is found to be small, corresponding to less than about 2 mm yr^?1 of baseline change between Southern India and the Himalayan piedmont. Based on an elastic dislocation model of interseismic strain and taking into account the uncertainty on India plate motion, the mean convergence rate across Central and Eastern Nepal is estimated to 19±2.5 mm yr^?1, (at the 67% confidence level). The main himalayan thrust (MHT) fault was found to be locked from the surface to a depth of about 20 km over a width of about 115 km. In these regions, the model parameters are well constrained, thanks to the long and continuous time-series from the permanent GPS as well as DORIS data. Further west, a convergence rate of 13.4±5 mm yr^?1, as well as a fault zone, locked over 150 km, are proposed. The slight discrepancy between the geologically estimated deformation rate of 21±1.5 mm yr^?1 and the 19±2.5 mm yr^?1 geodetic rate in Central and Eastern Nepal, as well as the lower geodetic rate in Western Nepal compared to Eastern Nepal, places bounds on possible temporal variations of the pattern and rate of strain in the period between large earthquakes in this region.     

Time Transfer for TAI Using a Geodetic Receiver: An Example with the Ashtech ZXII-T

The International Atomic Time scale (TAI) is computed by the Bureau International des Poids et Mesures (BIPM) from a set of atomic clocks distributed in about 40 time laboratories around the world. The time transfer between these remote clocks is mostly performed by the so-called GPS common view method: The clocks are connected to a GPS time receiver whose internal software computes the offsets between the remote clocks and GPS time. These data are collected in a standard formal called CCTF. In the present study we develop both the procedure and the software tool that allows us to generate the CCTF files needed for time transfer to TAI, using RINEX files produced by geodetic receivers driven by an external frequency. The CCTF files are then generated from the RINEX observation files. The software is freely available at ftp://omaftp.oma.be/dist/astro/time/RINEX_CCTF. Applied to IGS (International GPS Service) receivers, this procedure will provide a direct link between TAI and the IGS clock combination. We demonstrate here the procedure using the RINEX files from the Ashtech Metronome (ZXII-T) GPS receiver, to which we apply the conventional analysis to compute the CCTF data. We compared these results with the CCTF files produced by a time receiver R100-30T from 3S-Navigation. We also used this comparison with the results of a calibrated time receiver to determine the hardware delay of the geodetic receiver. ? 2001 John Wiley & Sons, Inc.     

Time and Frequency Transfer Using GPS Codes and Carrier Phases: Onsite Experiments

Recent studies have shown the capabilities of Global Positioning System (GPS) carrier phases for frequency transfer based on the observations from geodetic GPS receivers driven by stable atomic clocks. This kind of receiver configuration is the kind primarily used within the framework of the International GPS Service (IGS). The International GPS Service/Bureau International des Poids et Mesures (IGS/BIPM) pilot project aims at taking advantage of these GPS receivers to enlarge the network of Time Laboratories contributing to the realization of the International Atomic Time (TAI). In this article, we outline the theory necessary to describe the abilities and limitations of time and frequency transfer using the GPS code and carrier phase observations. We report on several onsite tests and evaluate the present setup of our 12-channel IGS receiver (BRUS), which uses a hydrogen maser as an external frequency reference, to contribute to the IGS/BIPM pilot project. In the initial experimental setup, the receivers had a common external frequency reference; in the second setup, separate external frequency references were used. Independent external clock monitoring provided the necessary information to validate the results. Using two receivers with a common frequency reference and connected to the same antenna, a zero baseline, we were able to use the carrier phase data to derive a frequency stability of 6 × 10⁻¹⁶ for averaging times of one day. The main limitation in the technique originates from small ambient temperature variations of a few degrees Celsius. While these temperature variations have no effect on the functioning of the GPS receiver within the IGS network, they reduce the capacities of the frequency transfer results based on the carrier phase data. We demonstrate that the synchronization offset at the initial measurement epoch can be estimated from a combined use of the code and carrier phase observations. In our test, the discontinuity between two consecutive days was about 140 ps. ? 1999 John Wiley & Sons, Inc.     

Application of geodetic receivers to timing and time transfer

Two methods for smoothing pseudorange observable by Carrier and Doppler are discussed. Then the procedure based on the RINEX observation files is tested using the Ashtech Z-XII3T geodetic receivers driven by a stable external frequency at UNSO. This paper proposes to adapt this procedure for the links between geodetic receivers, in order to take advantage of theP codes available onL1 andL2. This new procedure uses the 30-second RINEX observations files, the standard of the International GPS Service (IGS), and processes the ionosphere-free combination of the codesP1 andP2; the satellite positions are deduced from the IGS rapid orbits, available after two days.     

Application of geodetic receivers to timing and time transfer

Two methods for smoothing pseudorange observable by Carrier and Doppler are discussed. Then the procedure based on the RINEX observation files is tested using the Ashtech Z-XII3T geodetic receivers driven by a stable external frequency at UNSO. This paper proposes to adapt this procedure for the links between geodetic receivers, in order to take advantage of the P codes available on L1 and L2, This new procedure uses the 30-second RINEX observations files, the standard of the International GPS Service (IGS), and processes the ionosphere-free combination of the codes P1 and P2 ; the satellite positions are deduced from the IGS rapid orbits, available after two days.     

Improving the estimation of fractional-cycle biases for ambiguity resolution in precise point positioning

Ambiguity resolution dedicated to a single global positioning system (GPS) station can improve the accuracy of precise point positioning. In this process, the estimation accuracy of the narrow-lane fractional-cycle biases (FCBs), which destroy the integer nature of undifferenced ambiguities, is crucial to the ambiguity-fixed positioning accuracy. In this study, we hence propose the improved narrow-lane FCBs derived from an ambiguity-fixed GPS network solution, rather than the original (i.e. previously proposed) FCBs derived from an ambiguity-float network solution. The improved FCBs outperform the original FCBs by ensuring that the resulting ambiguity-fixed daily positions coincide in nature with the state-of-the-art positions generated by the International GNSS Service (IGS). To verify this improvement, 1?year of GPS measurements from about 350 globally distributed stations were processed. We find that the original FCBs differ more from the improved FCBs when fewer stations are involved in the FCB estimation, especially when the number of stations is less than 20. Moreover, when comparing the ambiguity-fixed daily positions with the IGS weekly positions for 248 stations through a Helmert transformation, for the East component, we find that on 359 days of the year the daily RMS of the transformed residuals based on the improved FCBs is smaller by up to 0.8?mm than those based on the original FCBs, and the mean RMS over the year falls evidently from 2.6 to 2.2?mm. Meanwhile, when using the improved rather than the original FCBs, the RMS of the transformed residuals for the East component of 239 stations (i.e. 96.4% of all 248 stations) is clearly reduced by up to 1.6?mm, especially for stations located within a sparse GPS network. Therefore, we suggest that narrow-lane FCBs should be determined with ambiguity-fixed, rather than ambiguity-float, GPS network solutions.     

Robust combination of IGS analysis center GLONASS clocks

GPS Solutions - The International GNSS Service (IGS) Analysis Centers (ACs) generate precise GNSS products by integrating tracking data from globally distributed IGS stations. The ACs’...     

Quality assessment of GPS reprocessed terrestrial reference frame

The International GNSS Service (IGS) contributes to the construction of the International Terrestrial Reference Frame (ITRF) by submitting time series of station positions and Earth Rotation Parameters (ERP). For the first time, its submission to the ITRF2008 construction is based on a combination of entirely reprocessed GPS solutions delivered by 11 Analysis Centers (ACs). We analyze the IGS submission and four of the individual AC contributions in terms of the GNSS frame origin and scale, station position repeatability and time series seasonal variations. We show here that the GPS Terrestrial Reference Frame (TRF) origin is consistent with Satellite laser Ranging (SLR) at the centimeter level with a drift lower than 1 mm/year. Although the scale drift compared to Very Long baseline Interferometry (VLBI) and SLR mean scale is smaller than 0.4 mm/year, we think that it would be premature to use that information in the ITRF scale definition due to its strong dependence on the GPS satellite and ground antenna phase center variations. The new position time series also show a better repeatability compared to past IGS combined products and their annual variations are shown to be more consistent with loading models. The comparison of GPS station positions and velocities to those of VLBI via local ties in co-located sites demonstrates that the IGS reprocessed solution submitted to the ITRF2008 is more reliable and precise than any of the past submissions. However, we show that some of the remaining inconsistencies between GPS and VLBI positioning may be caused by uncalibrated GNSS radomes.     

Estimation of phase center corrections for GLONASS-M satellite antennas

Driven by the comprehensive modernization of the GLONASS space segment and the increased global availability of GLONASS-capable ground stations, an updated set of satellite-specific antenna phase center corrections for the current GLONASS-M constellation is determined by processing 84 weeks of dual-frequency data collected between January 2008 and August 2009 by a worldwide network of 227 GPS-only and 115 combined GPS/GLONASS tracking stations. The analysis is performed according to a rigorous combined multi-system processing scheme providing full consistency between the GPS and the GLONASS system. The solution is aligned to a realization of the International Terrestrial Reference Frame 2005. The estimated antenna parameters are compared with the model values currently used within the International GNSS Service (IGS). It is shown that the z-offset estimates are on average 7 cm smaller than the corresponding IGS model values and that the block-specific mean value perfectly agrees with the nominal GLONASS-M z-offset provided by the satellite manufacturer. The existence of azimuth-dependent phase center variations is investigated and uncertainties in the horizontal offset estimates due to mathematical correlations and yaw-attitude modeling problems during eclipse seasons are addressed. Finally, it is demonstrated that the orbit quality benefits from the updated GLONASS-M antenna phase center model and that a consistent set of satellite antenna z-offsets for GPS and GLONASS is imperative to obtain consistent GPS- and GLONASS-derived station heights.     

基于我国全球连续监测与评估系统的ERP测定的精度分析

利用GPS观测资料解算地球自转参数,用全球均匀分布的22个IGS跟踪站（IGS05）的连续观测资料估计地球自转参数（ERP）,并与IERSC04（UTC0时）的结果相比较,二者相差很小,均在IERS的ERP估计精度范围之内。基于即将建成的COMPASS全球连续监测与评估系统跟踪站,选择其网的8个IGS跟踪站的资料进行了解算并进行了分析和比对。     

A model of present-day tectonic plate motions from 12 years of DORIS measurements

In the frame of the International DORIS Service (IDS), the Laboratoire d’Etudes en Géophysique et Océanographie Spatiales (LEGOS)/Collecte Localisation Satellites (CLS) Analysis Center (LCA) processes DORIS measurements from the SPOT, TOPEX/Poseidon and Envisat satellites and provides weekly station coordinates of the whole network to the IDS. Based on DORIS measurements, the horizontal and vertical velocities of 57 DORIS sites are computed. The 3D positions and velocities of the stations with linear motion are estimated simultaneously from the 12-year (1993–2004) combined normal equation matrix. We include 35 DORIS sites assumed to be located in the stable zones of 9 tectonic plates. For the motion of these plates, we propose a model (LCAVEL-1) of angular velocities in the ITRF2000 reference frame. Based on external comparison with the most recent global plate models (PB2002, REVEL, GSRM-1 and APKIM2000) and on internal analysis, we estimate an average velocity error of the DORIS solution of less than 3 mm/year. The LCAVEL-1 model presents new insights of the Somalia/Nubia pair of plates, as the DORIS technique has the advantage of having a few stations located on those two plates. We also computed (and provide in this article) the horizontal motion of the sites located close to plate boundaries or in the deformation zones defined in contemporary models. These computations could be used in further analysis for these particular regions of the Earth not moving as rigid plates.