Reducing the draconitic errors in GNSS geodetic products

Systematic errors at harmonics of the GPS draconitic year have been found in diverse GPS-derived geodetic products like the geocenter $Z$ -component, station coordinates, $Y$ -pole rate and orbits (i.e. orbit overlaps). The GPS draconitic year is the repeat period of the GPS constellation w.r.t. the Sun which is about 351 days. Different error sources have been proposed which could generate these spurious signals at the draconitic harmonics. In this study, we focus on one of these error sources, namely the radiation pressure orbit modeling deficiencies. For this purpose, three GPS+GLONASS solutions of 8 years (2004–2011) were computed which differ only in the solar radiation pressure (SRP) and satellite attitude models. The models employed in the solutions are: (1) the CODE (5-parameter) radiation pressure model widely used within the International GNSS Service community, (2) the adjustable box-wing model for SRP impacting GPS (and GLONASS) satellites, and (3) the adjustable box-wing model upgraded to use non-nominal yaw attitude, specially for satellites in eclipse seasons. When comparing the first solution with the third one we achieved the following in the GNSS geodetic products. Orbits: the draconitic errors in the orbit overlaps are reduced for the GPS satellites in all the harmonics on average 46, 38 and 57 % for the radial, along-track and cross-track components, while for GLONASS satellites they are mainly reduced in the cross-track component by 39 %. Geocenter $Z$ -component: all the odd draconitic harmonics found when the CODE model is used show a very important reduction (almost disappearing with a 92 % average reduction) with the new radiation pressure models. Earth orientation parameters: the draconitic errors are reduced for the $X$ -pole rate and especially for the $Y$ -pole rate by 24 and 50 % respectively. Station coordinates: all the draconitic harmonics (except the 2nd harmonic in the North component) are reduced in the North, East and Height components, with average reductions of 41, 39 and 35 % respectively. This shows, that part of the draconitic errors currently found in GNSS geodetic products are definitely induced by the CODE radiation pressure orbit modeling deficiencies.     

Can GPS-Derived Surface Loading Bridge a GRACE Mission Gap?

We investigated two ‘gap-filler’ methods based on GPS-derived low-degree surface loading variations (GPS-I and GPS-C) and a more simple method (REF-S) which extends a seasonal harmonic variation into the expected Gravity Recovery and Climate Experiment (GRACE) mission gap. We simulated two mission gaps in a reference solution (REF), which is derived from a joint inversion of GRACE (RL05) data, GPS-derived surface loading and simulated ocean bottom pressure. The GPS-I and GPS-C methods both have a new type of constraint applied to mitigate the lack of GPS station network coverage over the ocean. To obtain the GPS-C solution, the GPS-I method is adjusted such that it fits the reference solution better in a 1.5 year overlapping period outside of the gap. As can be expected, the GPS-I and GPS-C solutions contain larger errors compared to the reference solution, which is heavily constrained by GRACE. Within the simulated gaps, the GPS-C solution generally fits the reference solution better compared to the GPS-I method, both in terms of spherical harmonic loading coefficients and in terms of selected basin-averaged hydrological mass variations. Depending on the basin, the RMS-error of the water storage variations (scaled for leakage effects) ranges between 1.6 cm (Yukon) and 15.3 cm (Orinoco). In terms of noise level, the seasonal gap-filler method (REF-S) even outperforms the GPS-I and GPS-C methods, which are still affected by spatial aliasing problems. However, it must be noted that the REF-S method cannot be used beyond the study of simple harmonic seasonal variations.     

Rapid crustal uplift in Patagonia due to enhanced ice loss

A vertical crustal uplift rate of 39 mm yr^? 1 is measured between 2003 and 2006 using Global Positioning System (GPS) measurements at the northeastern edge of the Southern Patagonia Icefield (SPI). This is the largest present-day glacial isostatic rate ever recorded. The combination of SPI's rapid melting and the unique regional slab-window tectonics that promotes a relatively low viscosity, is central to our interpretation of the observations. The two effects lead to a strong interaction of short relaxation times with ice loads that change on a comparable time scale. The profile of GPS observations link ice loss to the soft viscoelastic isostatic flow response over the time scale of the Little Ice Age (LIA), including ice loss in the period of observation. The agreement of the results with our model predictions strongly suggests the large crustal uplift in Patagonia is due an accelerated glacier wasting since the termination of the LIA and that the effective regional mantle viscosity is near 4.0–8.0 × 10¹⁸ Pa s. A century-long diminution of the icefields, at rates that are about 1/4 – 1/2 the contemporary loss rates, is consistent with multidecadal-scale temperature trends estimated for the past 50–100 years and is, in fact, a key feature in any model capable of explaining the uplift data.     

Quality of reprocessed GPS satellite orbits

High-precision Global Positioning System (GPS) satellite orbits are one of the core products of the International GNSS Service (IGS). Since the establishment of the IGS in 1994, the quality and consistency of the IGS orbits has steadily been improved by advances in the modeling of GPS observations. However, due to these model improvements and reference frame changes, the time series of operational orbits are inhomogeneous and inconsistent. This problem can only be overcome by a complete reprocessing starting with the raw observation data. The quality of reprocessed GPS satellite orbits for the time period 1994–2005 will be assessed in this paper. Orbit fits show that the internal consistency of the orbits could be improved by a factor of about two in the early years. Comparisons with the operational IGS orbits show clear discontinuities whenever the reference frame was changed by the IGS. The independent validation with Satellite Laser Ranging (SLR) residuals shows an improvement of up to 30% whereas a systematic bias of 5 cm still persists.     

Assessing the Current Evolution of the Greenland Ice Sheet by Means of Satellite and Ground-Based Observations

The present study utilises different satellite and ground-based geodetic observations in order to assess the current evolution of the Greenland Ice Sheet (GIS). Satellite gravimetry data acquired by the Gravity Recovery and Climate Experiment are used to derive ice-mass changes for the period from 2003 to 2012. The inferred time series are investigated regarding long-term, seasonal and interannual variations. Laser altimetry data acquired by the Ice, Cloud, and land Elevation Satellite (ICESat) are utilised to solve for linear and seasonal changes in the ice-surface height and to infer independent mass-change estimates for the entire GIS and its major drainage basins. We demonstrate that common signals can be identified in the results of both sensors. Moreover, the analysis of a Global Positioning System (GPS) campaign network in West Greenland for the period 1995–2007 allows us to derive crustal deformation caused by glacial isostatic adjustment (GIA) and by present-day ice-mass changes. ICESat-derived elastic crustal deformations are evaluated comparing them with GPS-observed uplift rates which were corrected for the GIA effect inferred by model predictions. Existing differences can be related to the limited resolution of ICESat. Such differences are mostly evident in dynamical regions such as the Disko Bay region including the rapidly changing Jakobshavn Isbræ, which is investigated in more detail. Glacier flow velocities are inferred from satellite imagery yielding an accelerated flow from 1999 to 2012. Since our GPS observations cover a period of more than a decade, changes in the vertical uplift rates can also be investigated. It turns out that the increased mass loss of the glacier is also reflected by an accelerated vertical uplift.