DORIS and GPS monitoring of the Gavdos calibration site in Crete

Due to its specific geographical location as well as its geodetic equipment (DORIS, GNSS, microwave transponder and tide gauges), the Gavdos station in Crete, Greece is one of the very few sites around the world used for satellite altimetry calibration. To investigate the quality of the Gavdos geodetic coordinates and velocities, we analyzed and compared here DORIS and GPS-derived results obtained during several years of observations. The DORIS solution is the latest ignwd11 solution at IGN, expressed in ITRF2008, while the GPS solution was obtained using the GAMIT software package. Current results show that 1–2 mm/yr agreement can be obtained for 3-D velocity, showing a good agreement with current geophysical models. In particular, the agreement obtained for the vertical velocity is around 0.3–0.4 mm/yr, depending on the terrestrial reference frame. As a by-product of these geodetic GPS and DORIS results, Zenith Tropospheric Delays (ZTDs) estimations were also compared in 2010 between these two techniques, and compared to ECMWF values, showing a 6.6 mm agreement in dispersion without any significant difference between GPS and DORIS (with a 97.6% correlation), but with a 13–14 mm agreement in dispersion when comparing to ECMWF model (with only about 90% correlation for both techniques). These tropospheric delay estimations could also provide an external calibration of the tropospheric correction used for the geophysical data of satellite altimetry missions.     

Temporal and spatial stochastic behaviour of high-frequency slant tropospheric delays from simulations and real GPS data

In this paper, a turbulence theory-based simulation procedure for slant tropospheric delay variations is presented. Based on this procedure tropospheric delay variations are simulated for three different geometric scenarios. The stochastic behaviour of the generated time series is assessed in terms of temporal structure functions. It is shown that the temporal structure functions – in general – follow a 5/3 to 2/3 power-law behaviour. Deviations from this behaviour due to the complex interaction between varying observation geometry and atmospheric/turbulent conditions are discussed.     

An Integer Precise Point Positioning technique for sea surface observations using a GPS buoy

GPS data dedicated to sea surface observation are usually processed using differential techniques. Unfortunately, the precision of resulting kinematic positions is baseline-length dependent. So, high precision sea surface observations using differential GPS techniques are limited to coasts, lakes, and rivers. Recent improvements in GPS satellite products (orbits, clocks, and phase biases) make phase ambiguity fixing at the zero difference level achievable and opens up the observation of the sea surface without geographical constraints. This paper recalls the concept of the Integer Precise Point Positioning technique and discusses the precision of GPS buoy positioning. A sequential version of the GINS software has been implemented to achieve single epoch GPS positioning. We used 1 Hz data from a two week GPS campaign conducted in the Kerguelen Islands. A GPS buoy has been moored close to a radar gauge and 90 m away from a permanent GPS station. This infrastructure offers the opportunity to compare both kinematic Integer Precise Point Positioning and classical differential GPS positioning techniques to in situ radar gauge data. We found that Precise Point Positioning results are not significantly biased with respect to radar gauge data and that horizontal time series are consistent with differential processing at the sub-centimetre precision level. Nevertheless, standard deviations of height time series with respect to radar gauge data are typically [4–5] cm. The dominant driver for noise at this level is attributed to errors in tropospheric estimates which propagate into position solutions.     

Effects of physical correlations on long-distance GPS positioning and zenith tropospheric delay estimates

The global positioning system (GPS) has become an essential tool for the high precision navigation and positioning. The quality of GPS positioning results mainly depends on the model’s formulations regarding GPS observations, including both a functional model, which describes the mathematical relationships between the GPS measurements and unknown parameters, and a stochastic model, which reflects the physical properties of the measurements. Over the past two decades, the functional models for GPS measurements have been investigated in considerable detail. However, the stochastic models of GPS observation data are simplified, assuming that all the GPS measurements have the same variance and are statistically independent. Such assumptions are unrealistic. Although a few studies of GPS stochastic models were performed, they are restricted to short baselines and short time session lengths. In this paper, the stochastic modeling for GPS long-baseline and zenith tropospheric delay (ZTD) estimates with a 24-h session is investigated using the residual-based and standard stochastic models. Results show that using the different stochastic modelling methods, the total differences can reach as much as 3–6 mm in the baseline component, especially in the height component, and 10 mm in the ZTD estimation. Any misspecification in the stochastic models will result in unreliable GPS baseline and ZTD estimations. Using the residual-based stochastic model, not only the precision of GPS baseline and ZTD estimation is obviously improved, but also the baseline and ZTD estimations are closer to the reference value.     

A DORIS determination of the absolute velocities of the Sorsdal and Mellor glaciers in Antarctica

The Lambert–Amery System is the largest glacier–ice shelf system in East Antarctica, draining a significant portion of the ice sheet. Variation in ice sheet discharge from Antarctica or Greenland has an impact on the rate of change in global mean sea level; which is a manifestation of climate change. In conjunction with a measure of ice thickness change, ice sheet discharge can be monitored by determining the absolute velocities of these glaciers.     

Improving the estimation of zenith dry tropospheric delays using regional surface meteorological data

Global Navigation Satellite Systems (GNSS) are emerging as possible tools for remote sensing high-resolution atmospheric water vapour that improves weather forecasting through numerical weather prediction models. Nowadays, the GNSS-derived tropospheric zenith total delay (ZTD), comprising zenith dry delay (ZDD) and zenith wet delay (ZWD), is achievable with sub-centimetre accuracy. However, if no representative near-site meteorological information is available, the quality of the ZDD derived from tropospheric models is degraded, leading to inaccurate estimation of the water vapour component ZWD as difference between ZTD and ZDD. On the basis of freely accessible regional surface meteorological data, this paper proposes a height-dependent linear correction model for a priori ZDD. By applying the ordinary least-squares estimation (OLSE), bootstrapping (BOOT), and leave-one-out cross-validation (CROS) methods, the model parameters are estimated and analysed with respect to outlier detection. The model validation is carried out using GNSS stations with near-site meteorological measurements. The results verify the efficiency of the proposed ZDD correction model, showing a significant reduction in the mean bias from several centimetres to about 5 mm. The OLSE method enables a fast computation, while the CROS procedure allows for outlier detection. All the three methods produce consistent results after outlier elimination, which improves the regression quality by about 20% and the model accuracy by up to 30%.     

Demonstration of the use of the Doppler Orbitography and Radio positioning Integrated by Satellite (DORIS) measurements to validate GPS ionospheric imaging

There is a lack of independent ionospheric data that can be used to validate GPS imaging results at mid latitudes over severe storm times. Doppler Orbitography and Radio positioning Integrated by Satellite (DORIS), a global network of dual-frequency ground to satellite observations, provides this missing data and here is employed as verification to show the accuracy of the ionospheric GPS images in terms of the total electron content (TEC). In this paper, the large-scale ionospheric structures that appeared during the strong geomagnetic storm of 20 November 2003 are reconstructed with a GPS tomographic algorithm, known as MIDAS, and validated with DORIS TEC measurements. The main trough shown in an extreme equatorward position in the ionospheric imaging over mainland Europe is confirmed by DORIS satellite measurements. Throughout the disturbed day, the variations of relative slant TECs between DORIS data and MIDAS results agree quite well, with the average of the mean differences about 2 TECu. We conclude that as a valuable supplement to GPS data, DORIS ionospheric measurements can be used to analyse TEC variations with a relatively high resolution, ∼10 s in time and tens of kilometres in space. This will be very helpful for identification of some highly dynamic structures in the ionosphere found at mid-latitudes, such as the main trough, TID (Travelling Ionospheric Disturbances) and SED (Storm Enhanced Density), and could be used as a valuable auxiliary data source in ionospheric imaging.     

Vertical crustal motion derived from satellite altimetry and tide gauges,and comparisons with DORIS measurements

A somewhat unorthodox method for determining vertical crustal motion at a tide-gauge location is to difference the sea level time series with an equivalent time series determined from satellite altimetry. To the extent that both instruments measure an identical ocean signal, the difference will be dominated by vertical land motion at the gauge. We revisit this technique by analyzing sea level signals at 28 tide gauges that are colocated with DORIS geodetic stations. Comparisons of altimeter-gauge vertical rates with DORIS rates yield a median difference of 1.8 mm yr⁻¹ and a weighted root-mean-square difference of 2.7 mm yr⁻¹. The latter suggests that our uncertainty estimates, which are primarily based on an assumed AR(1) noise process in all time series, underestimates the true errors. Several sources of additional error are discussed, including possible scale errors in the terrestrial reference frame to which altimeter-gauge rates are mostly insensitive. One of our stations, Malè, Maldives, which has been the subject of some uninformed arguments about sea-level rise, is found to have almost no vertical motion, and thus is vulnerable to rising sea levels.     

An estimation of tropospheric corrections using GPS and synoptic data: Improving Urmia Lake water level time series from Jason-2 and SARAL/AltiKa satellite altimetry

Tropospheric correction is one of the most important corrections in satellite altimetry measurements. Tropospheric wet and dry path delays have strong dependence on temperature, pressure and humidity. Tropospheric layer has particularly high variability over coastal regions due to humidity, wind and temperature gradients. Depending on the extent of water body and wind conditions over an inland water, Wet Tropospheric Correction (WTC) is within the ranges from a few centimeters to tens of centimeters. Therefore, an extra care is needed to estimate tropospheric corrections on the altimetric measurements over inland waters. This study assesses the role of tropospheric correction on the altimetric measurements over the Urmia Lake in Iran. For this purpose, four types of tropospheric corrections have been used: (i) microwave radiometer (MWR) observations, (ii) tropospheric corrections computed from meteorological models, (iii) GPS observations and (iv) synoptic station data. They have been applied to Jason-2 track no. 133 and SARAL/AltiKa track no. 741 and 356 corresponding to 117–153 and the 23–34 cycles, respectively. In addition, the corresponding measurements of PISTACH and PEACHI, include new retracking method and an innovative wet tropospheric correction, have also been used. Our results show that GPS observation leads to the most accurate tropospheric correction. The results obtained from the PISTACH and PEACHI projects confirm those obtained with the standard SGDR, i.e., the role of GPS in improving the tropospheric corrections. It is inferred that the MWR data from Jason-2 mission is appropriate for the tropospheric corrections, however the SARAL/AltiKa one is not proper because Jason-2 possesses an enhanced WTC near the coast. Furthermore, virtual stations are defined for assessment of the results in terms of time series of Water Level Height (WLH). The results show that GPS tropospheric corrections lead to the most accurate WLH estimation for the selected virtual stations, which improves the accuracy of the obtained WLH time series by about 5%.     

Near real-time estimation of tropospheric water vapour content from ground based GNSS data and its potential contribution to weather now-casting in Austria

The importance of high resolution meteorological analysis of the atmosphere increased over the past years. A detailed analysis of the humidity field is an important precondition for a better monitoring of local and regional extreme precipitation events and for forecasts with improved spatial resolution. For this reason, the Austrian Meteorological Agency (ZAMG) is operating the spatial and temporal high resolution INCA system (Integrated Now-casting through Comprehensive Analysis) since begin of 2005. Errors in this analysis occur mainly in the areas of rapidly changing and hard to predict weather conditions or rugged topography with extreme differences in height such as the alpine area of Austria. The aim of this work is to provide GNSS based measurements of the tropospheric water vapour content with a temporal resolution of 1 h and a temporal delay of less than 1 h to assimilate these estimates into the INCA system. Additional requirement is an accuracy of better than 1 mm of the precipitable water (PW) estimates.     

Dynamics of global electron content in 1998–2005 derived from global GPS data and IRI modeling

We analyzed the dynamics of global electron content (GEC) for the period 1998–2005 and compared the estimated GEC with variations of the 10.7-cm solar radio emission and with and with GEC values obtained with IRI-2001. We found a strong resemblance between the curves’ shapes for the experimental and modeled GEC: strong semiannual variations are discernible in these series and both curves tend to increase the absolute GEC value during the period of maximum of solar activity. However, there are some significant distinctions, such as absence of 27-day fluctuations in the series of GEC computed by the IRI-2001. On the contrary, observational GEC reflects well dynamics of solar activity: 27-day variations of GEC are very similar to the ones of the index F10.7, but GEC undergoes a lagging of about of 30–60 h as compared to value of the F10.7 index. The relative amplitude of 27-day variations decreases from 8% at the rising and falling solar activity to 2% at the period of its maximum.     

Propagation of plasma bubbles observed in Brazil from GPS and airglow data

Equatorial spread-F is a common occurrence in the equatorial ionosphere that is associated with large variations in plasma density that often cause scintillation and interference in communication signals. These events are known to result from Rayleigh–Taylor instability, but the day-to-day variability of their occurrence is not well understood. The triggering mechanism of plasma depletions is still a matter of debate, but may be linked to gravity waves that under favorable conditions propagate to the middle atmosphere. Understanding the triggering of ESF was the focus of the SpreadFEx campaign near Brasilia, Brazil in 2005. The campaign provided co-located airglow and GPS observations to study the onset of plasma depletions and their evolution as they traversed the region. Comparisons between the 630.0 nm airglow data and GPS data demonstrate the ability of the compact dual frequency GPS array to detect the plasma bubbles and retrieve reliable propagation characteristics of the depletions. In this case study, a plasma depletion was detected and moved over the array at velocities of 85–110 m/s, slowing as it moved towards the east. Correlation of consecutive airglow images gives consistent estimates of the eastward drift over the same time period. Mapping the airglow data to the GPS line-of-sight geometry allows direct comparison and reveals a resolvable westward tilt of the plasma depletion that may be due to vertical shear. The uniqueness of this study is the ability to resolve locally the characteristics of the plasma depletion without relying on assumptions about the mapping of the depletion along magnetic field lines to large latitudinal distances. It presents new information for understanding ESF development and the development of depletions strong enough to produce scintillation.     

TEC derived from some GPS stations in Nigeria and comparison with the IRI and NeQuick models

Total electron content (TEC) measured simultaneously using Global Positioning System (GPS) ionospheric monitors installed at some locations in Nigeria during the year 2011 (Rz = 55.7) was used to study the diurnal, seasonal, and annual TEC variations. The TEC exhibits daytime maximum, seasonal variation and semiannual variations. Measured TEC were compared with those predicted by the improved versions of the International Reference Ionosphere (IRI) and NeQuick models. The models followed the diurnal and seasonal variation patterns of the observed values of TEC. However, IRI model produced better estimates of TEC than NeQuick at all locations.     

Use of varying shell heights derived from ionosonde data in calculating vertical total electron content (TEC) using GPS – New method

The dispersive nature of the ionosphere makes it possible to measure its total electron content (TEC). Thus Global Positioning System, which uses dual-frequency radio signals, is an ideal system to measure TEC. When data from an ionosonde situated in polar region was observed, the height of an approximated thin shell of electrons (shell height) used in GPS studies was seen not to be fixed but rather changing with time. Here we introduce a new method in which we included the varying shell heights derived from the ionosonde to map the slant total electron content from GPS to obtain a more precise vertical total electron content of the ionosphere contrary to some previous methods which used fixed shell heights. In this paper we also compared the ionosonde derived TEC with the GPS derived vertical TEC (vTEC) values. These GPS vTEC values were obtained from GPS slant TEC (sTEC) measurements using both fixed shell height and varying shell heights (from ionosonde measurements). For the polar regions, the varying shell height approach produced better results than the fixed shell height and compared to exponential function, Chapman function seems to be a better function to model the topside ionosphere.     

Ionospheric anomalies of local earthquakes detected by GPS TEC measurements using data from Tashkent and Kitab stations

This paper reports the ionospheric anomalies observed during strong local earthquakes (M?5.0) which occurred mostly in and around Uzbekistan in seismically active zones, during years 2006 to 2009 within approximately 1000 km distance from the observing GPS stations located in Tashkent and Kitab, Uzbekistan. The solar and geomagnetic conditions were quiet during occurrence of the selected strong earthquakes. We produce Total Electron Content (TEC) time series over both sites and apply them to detect anomalous TEC signals preceding or accompanying the local earthquakes. The results show anomalous increase or decrease of TEC before or during the earthquakes. In general the anomalies occurred 1–7 days before the earthquakes as ionospheric electromagnetic precursors. To identify the anomalous values of TEC we calculated differential TEC (dTEC). dTEC is obtained by subtracting monthly averaged diurnal vTEC from the values of observed vTEC at each epoch. This procedure removes normal diurnal variations of vTEC. The present results are in good agreement with the previous observations on ionospheric earthquake precursors reported by various researchers.     

Plasmaspheric electron content derived from GPS TEC and FORMOSAT-3/COSMIC measurements: Solar minimum condition

The plasmaspheric electron content (PEC) was estimated by comparison of GPS TEC observations and FORMOSAT-3/COSMIC radio occultation measurements at the extended solar minimum of cycle 23/24. Results are retrieved for different seasons (equinoxes and solstices) of the year 2009. COSMIC-derived electron density profiles were integrated up to the height of 700 km in order to retrieve estimates of ionospheric electron content (IEC). Global maps of monthly median values of COSMIC IEC were constructed by use of spherical harmonics expansion. The comparison between two independent measurements was performed by analysis of the global distribution of electron content estimates, as well as by selection specific points corresponded to mid-latitudes of Northern America, Europe, Asia and the Southern Hemisphere. The analysis found that both kinds of observations show rather similar diurnal behavior during all seasons, certainly with GPS TEC estimates larger than corresponded COSMIC IEC values. It was shown that during daytime both GPS TEC and COSMIC IEC values were generally lower at winter than in summer solstice practically over all specific points. The estimates of PEC (h > 700 km) were obtained as a difference between GPS TEC and COSMIC IEC values. Results of comparative study revealed that for mid-latitudinal points PEC estimates varied weakly with the time of a day and reached the value of several TECU for the condition of solar minimum. Percentage contribution of PEC to GPS TEC indicated the clear dependence from the time with maximal values (more than 50–60%) during night-time and lesser values (25–45%) during day-time.     

GPS data processing of networks with mixed single- and dual-frequency receivers for deformation monitoring

In order to obtain crustal deformations of higher spatial resolution, existing GPS networks must be densified. This densification can be carried out using single-frequency receivers at moderate costs. However, ionospheric delay handling is required in the data processing. We adapt the Satellite-specific Epoch-differenced Ionospheric Delay model (SEID) for GPS networks with mixed single- and dual-frequency receivers. The SEID model is modified to utilize the observations from the three nearest dual-frequency reference stations in order to avoid contaminations from more remote stations. As data of only three stations are used, an efficient missing data constructing approach with polynomial fitting is implemented to minimize data losses. Data from large scale reference networks extended with single-frequency receivers can now be processed, based on the adapted SEID model. A new data processing scheme is developed in order to make use of existing GPS data processing software packages without any modifications. This processing scheme is evaluated using a sub-network of the German SAPOS network. The results verify that the new scheme provides an efficient way to densify existing GPS networks with single-frequency receivers.     

A gravity wave analysis near to the Andes Range from GPS radio occultation data and mesoscale numerical simulations: Two case studies

Global maps of potential wave energy per unit mass, recently performed with the Global Positioning System (GPS) Radio Occultation (RO) technique and different satellite missions (CHAMP and SAC-C since 2001, GRACE and COSMIC since 2006) revealed in Argentina, at the eastern side of the highest Andes Mountains, a considerable wave activity (WA) in comparison with other extra-tropical regions. The main gravity wave (GW) sources in this natural laboratory are deep convection (mainly during late Spring and Summer), topographic forcing and geostrophic adjustment.     

An analysis on combined GPS/COMPASS data quality and its effect on single point positioning accuracy under different observing conditions

With the rapid development of the COMPASS system, it is currently capable of providing regional navigation services. In order to test its data quality and performance for single point positioning (SPP), experiments have been conducted under different observing conditions including open sky, under trees, nearby a glass wall, nearby a large area of water, under high-voltage lines and under a signal transmitting tower. To assess the COMPASS data quality, the code multipath, cycle slip occurrence rate and data availability were analyzed and compared to GPS data. The datasets obtained from the experiments have also been utilized to perform combined GPS/COMPASS SPP on an epoch-by-epoch basis using unsmoothed single-frequency code observations. The investigation on the regional navigation performance aims at low-accuracy applications and all tests are made in Changsha, China, using the “SOUTH S82-C” GPS/COMPASS receivers. The results show that adding COMPASS observations can significantly improve the positioning accuracy of single-frequency GPS-only SPP in environments with limited satellite visibility. Since the COMPASS system is still in an initial operational stage, all results are obtained based on a fairly limited amount of data.     

Location of layered structures in the ionosphere and atmosphere by use of GPS occultation data

A method is introduced to locate the layered structures in the atmosphere and ionosphere based on simultaneous observations of radio wave intensity and phase variations in trans-ionospheric satellite-to-satellite links. The method determines location of a tangent point on the trans-ionospheric ray trajectory where gradient of refractivity is perpendicular to the ray trajectory and influence of a layered structure on radio wave parameters is maximal. An estimate of the location of a layer can be obtained from a combination of the phase and intensity variations. This new technique was applied to measurements provided during FORMOSAT-3 and CHAMP radio occultation (RO) missions. For the considered RO events the location of the inclined plasma layer in the lower ionosphere is found and the electron density distribution is retrieved. The method is checked by measuring the location of the tangent point on the ray trajectory in the neutral gas in the atmosphere. The results showed a fairly good agreement.