Felix Norman Teferle

Using continuous GPS and absolute gravity to separate vertical land movements and changes in sea-level at tide-gauges in the UK

Researchers investigating climate change have used historical tide-gauge measurements from all over the world to investigate the changes in sea-level that have occurred over the last century or so. However, such estimates are a combination of any true sea-level variations and any vertical movements of the land at the specific tide-gauge. For a tide- gauge record to be used to determine the climate related component of changes in sea-level, it is therefore necessary to correct for the vertical land movement component of the observed change in sea-level. In 1990, the Institute of Engineering Surveying and Space Geodesy and Proudman Oceanographic Laboratory started developing techniques based on the Global Positioning System (GPS) for measuring vertical land movements (VLM) at tide-gauges in the UK. This paper provides brief details of these early developments and shows how they led to the establishment of continuous GPS (CGPS) stations at a number of tide-gauges. The paper then goes on to discuss the use of absolute gravity (AG), as an independent technique for measuring VLM at tide-gauges. The most recent results, from CGPS time-series dating back to 1997 and AG time-series dating back to 1995/1996, are then used to demonstrate the complementarity of these two techniques and their potential for providing site-specific estimates of VLM at tide-gauges in the UK.

Modelling the glacial isostatic adjustment of the UK region

The glacial isostatic adjustment of the UK region has been considered in a number of recent studies. We have revisited this problem in order to: (i) highlight some key issues with regard to limitations in the ice modelling approach adopted in these studies and (ii) consider the constraints provided from observations of crustal motion available via continuous global positioning system monitoring. With regard to the first aim, we have found that: (i) previous studies have significantly overestimated ice thicknesses in regions where trim line field constraints were adopted and (ii) the duration of the glaciation phase of the UK ice sheet is a critical aspect of the model and that discrepancies in this model component have led to inconsistent inferences of Earth model parameters. With regard to the second aim, we have found that predictions of horizontal velocities (relative to a chosen site) based on a UK ice model calibrated to fit the regional sea-level database capture the geometry of the signal well but only account for 10% of the magnitude (for a range of Earth models).

Performance of Precise Point Positioning with Ambiguity Resolution for 1- to 4-Hour Observation Periods

Abstract Recent progress in integer ambiguity resolution at a single station has made it possible to achieve high positioning accuracy in static precise point positioning (PPP) using a short period of observations. In this paper, 12 stations across Europe are used to conduct short-period (i.e. one, two, three and four hours) static PPP with ambiguity resolution from Day 245 to 251 in 2007. It is demonstrated that, when over three hours of observations are used, PPP can achieve a success rate of 100% for ambiguity resolution, a 3D positioning accuracy of about 1.0 cm and an occurrence of less than 1.0% for degraded solutions. Moreover, for the fixed solutions, increasing the observation period hardly improves the horizontal positioning accuracy while still improving the vertical one. Therefore, it is proposed that at least three hours of observations should be used in the ambiguity-fixed static PPP if a reliable millimetre positioning accuracy is required in the engineering applications.

Using GPS to Separate Crustal Movements and Sea Level Changes at Tide Gauges in the UK

Changes in mean sea level recorded by tide gauges are corrupted by crustal movements, which can be of a similar order of magnitude. A network of permanent, continuous GPS (CGPS) stations has been established in the UK with five stations being sited at tide gauges. Data from these and four other CGPS stations have been analysed. A common mode filtering technique was successfully applied in order to reduce the effect of annually repeating signals on station velocity estimates. The effect of time-dependent correlations in the coordinate time series were accounted for when computing station velocity uncertainties.

A Bayesian Monte Carlo Markov Chain Method for Parameter Estimation of Fractional Differenced Gaussian Processes

We present a Bayesian Monte Carlo Markov Chain method to simultaneously estimate the spectral index and power amplitude of a fractional differenced Gaussian process at low frequency, in presence of white noise, and a linear trend and periodic signals. This method provides a sample of the likelihood function and thereby, using Monte Carlo integration, all parameters and their uncertainties are estimated simultaneously. We test this method with simulated and real Global Positioning System height time series and propose it as an alternative to optimization methods currently in use. Furthermore, without any mathematical proof, the results from the simulations suggest that this method is unaffected by the stationary regime and hence, can be used to check whether or not a time series is stationary.

Application of the dual-CGPS concept to monitoring vertical land movements at tide gauges ☆

Abstract Recently the concept of using dual-continuous GPS (dual-CGPS) stations for monitoring vertical land movements at tide gauge sites was proposed by a working group of the European Sea Level Observing System (EOSS). In this concept, one CGPS station is established at the tide gauge, in order to monitor the local vertical land movements, and a second CGPS station is established on stable rock, within a few kilometers of the tide gauge further inland. Since 1997, a network of CGPS stations has been established to monitor vertical land movements in the UK. For the dual-CGPS station concept, pairs formed by the CGPS stations at Newlyn tide gauge and Camborne and at Lowestoft tide gauge and Hemsby, have been analyzed over a time span of approximately 2.7 and 2.0 years respectively. It has been shown by numerous authors that temporal and spatial correlations in CGPS coordinate time series can introduce biases in the estimated station velocities and their uncertainties. By analyzing the coordinate time series of two CGPS stations close to each other, spatial correlations can be removed successfully by differencing and a cleaner, difference time series for each coordinate component can be obtained. In this paper, the observed coordinate time series of four CGPS stations have been used to produce synthetic coordinate time series of 6 years in length. Station velocities and their uncertainties have been derived for both observed and synthetic coordinate time series and compared in order to assess the capabilities of the dual-CGPS concept to better describe local and geophysical vertical land movements.

Land motion in the urban area of Nottingham observed by ENVISAT-1

Persistent scatterer interferometry (PSI) analysis of land motion in the City of Nottingham has been undertaken using advanced synthetic aperture radar (ASAR) data from the ENVISAT-1 satellite covering the period between November 2002 and February 2009, and some interesting areas of motion are evident from the results. The aims of this article are twofold: first to relate these areas to other sources of information such as geology, GPS, and precise levelling, and second to discuss a novel method of atmospheric correction for PSI using precise point positioning (PPP) derived GPS zenith wet delay (ZWD) estimates. PSI validation using the various sources of data is discussed with contextual information regarding historical activity. This article shows that land motion in an urban environment can be widespread and anisotropic, with many different causes, with the consequences that sparse networks of ground-based instruments may be incapable of monitoring land motion in such a complex environment. A ZWD correction method is also presented and results are shown, which improve upon those without the correction. The ZWD estimates are produced using an undifferenced global navigation satellite system (GNSS) processing technique known as PPP. This, coupled with the chosen PSI method, provides a low noise solution when compared to similar experiments involving the integration of GNSS and interferometric synthetic aperture radar (InSAR) data. The results show that there is good temporal and spatial correlation between the PPP ZWD estimates and the PSI-derived atmospheric phase contributions. Noticeable improvements in the PSI solution when comparing before and after processing runs to precise levelling observations are shown.

Optimum stochastic modeling for GNSS tropospheric delay estimation in real-time

In GNSS data processing, the station height, receiver clock and tropospheric delay (ZTD) are highly correlated to each other. Although the zenith hydrostatic delay of the troposphere can be provided with sufficient accuracy, zenith wet delay (ZWD) has to be estimated, which is usually done in a random walk process. Since ZWD temporal variation depends on the water vapor content in the atmosphere, it seems to be reasonable that ZWD constraints in GNSS processing should be geographically and/or time dependent. We propose to take benefit from numerical weather prediction models to define optimum random walk process noise. In the first approach, we used archived VMF1-G data to calculate a grid of yearly and monthly means of the difference of ZWD between two consecutive epochs divided by the root square of the time lapsed, which can be considered as a random walk process noise. Alternatively, we used the Global Forecast System model from National Centres for Environmental Prediction to calculate random walk process noise dynamically in real-time. We performed two representative experimental campaigns with 20 globally distributed International GNSS Service (IGS) stations and compared real-time ZTD estimates with the official ZTD product from the IGS. With both our approaches, we obtained an improvement of up to 10% in accuracy of the ZTD estimates compared to any uniformly fixed random walk process noise applied for all stations.

Impact of Antenna Phase Centre Calibrations on Position Time Series: Preliminary Results

Advances in GPS error modelling and the continued effort of re-processing have considerably decreased the scatter in position estimates over the last decade. The associated reduction of noise in derived position time series has revealed the presence of previously undetected periodic signals. It has been shown that these signals have frequencies related to the orbits of the GPS satellites. A number of potential sources for these periodicities at the draconitic frequency and its harmonics have already been suggested in the literature and include, e.g., errors in the sub-daily tidal models, multipath and unresolved integer ambiguities. Due to the geometrical relationship between the observing site and the orbiting satellite, deficiencies in the modelling of electromagnetic phase centres of receiving antennas have the potential to also contribute to the discovered periodic signals. The change from relative to absolute type mean antenna/radome calibrations within the International GNSS Service (IGS) led to a significant improvement, but the use of individual calibrations could possibly add further refinements to computed solutions. However, at this stage providing individual calibrations for all IGS stations is not feasible. Furthermore, antenna near-field electromagnetic effects might outweigh the benefits of individual calibrations once an antenna is permanently installed. In this study, we investigate the differences between position estimates obtained using individual and type mean antenna/radome calibrations as used by the IGS community. We employ position time series derived from precise point positioning (PPP) as implemented in two scientific GNSS software packages. Our results suggest that the calibration differences propagate directly into the position estimates, affecting both sub-daily and daily results and yielding periodic variations. The sub-daily variations have periods close to half a sidereal day and one sidereal day with peak-to-peak amplitudes of up to 10 mm in all position components. The stacked power spectra of the daily difference time series reveal peaks at the GPS draconitic frequency and its harmonics with peak-to-peak amplitudes of up to 1 mm. Although these results are still preliminary, they confirm that small differences between individual and type mean antenna/radome calibrations propagate into position time series and may be partly responsible for the spurious signals with draconitic frequency and its harmonics.

New Estimates of Present-Day Crustal/Land Motions in the British Isles Based on the BIGF Network

In this study we present results from a recent re-processing effort that included data from more than 120 continuous Global Positioning System (CGPS) stations in the British Isles for the period from 1997 to 2008. Not only was the CGPS network dramatically densified from previous investigations by the authors, it now also includes, for the first time, stations in Northern Ireland, providing new constraints on glacio-isostatic processes active in the region. In our processing strategy we apply a combination of re-analysed satellite orbit and Earth rotation products together with updated models for absolute satellite and receiver antenna phase centers, and for the computation of atmospheric delays. Our reference frame implementation uses a semi-global network of 37 stations, to align our daily position estimates, using a minimal constraints approach, to ITRF2005. This network uses a combination of current IGS reference frame stations plus additional IGS stations in order to provide similar network geometries throughout the complete time span. The derived horizontal and vertical station velocities are used to investigate present-day crustal/land motions in the British Isles. This first solution provides the basis for our contribution to the Working Group on Regional Dense Velocity Fields, 2007–2011 of the International Association of Geodesy Subcommission 1.3 on Regional Reference Frames.