Noise analysis of continuous GPS coordinate time series for CMONOC

The Crustal Movement Observation Network of China (CMONOC) is one of the major scientific infrastructures, mainly using Global Positioning System (GPS) measurements, to monitor crustal deformation in the Chinese mainland. In this paper, decade-long coordinate time series of 26 continuous GPS sites of CMONOC are analyzed for their noise content using maximum likelihood estimation (MLE). We study the noise properties of continuous GPS time series of CMONOC for the unfiltered, filtered solutions and also the common mode signals in terms of power law plus white noise model. In the spatial filtering, we remove for every time series a common mode error that was estimated from a modified stacking of position residuals from other sites within ∼1000 km of the selected site. We find that the common mode signal in our network has a combination of spatially correlated flicker noise and a common white noise with large spatial extent. We demonstrate that for the unfiltered solutions of CMONOC continuous GPS sites the main colored noise is a flicker process, with a mean spectral index of ∼1. For the filtered solutions, the main colored noise is a general power law process, indicating that a major number of the filtered regional solutions have a combination of noise sources or local effects. The velocity uncertainties from CMONOC continuous GPS coordinate time series may be underestimated by factors of 8–16 if a pure white noise model is assumed. In addition, using a white plus flicker noise model, the median values of velocity errors for the unfiltered solutions are 0.16 (north), 0.17 (east) and 0.58 (vertical) mm/yr, and the median values for the filtered solutions are 0.09 (north), 0.10 (east) and 0.40 (vertical) mm/yr.     

Ionospheric scintillations at Guilin detected by GPS ground-based and radio occultation observations

The occurrence of ionospheric scintillations with S4 ? 0.2 was studied using GPS measurements at Guilin, China (25.29°N, 110.33°E; geomagnetic: 15.04°N, 181.98°E), a station located near the northern crest of the equatorial anomaly. The results are presented for data collected from January 2009 to March 2010. The results show that nighttime amplitude scintillations only took place in February and March of the considered years, while daytime amplitude scintillations occurred in August and December of 2009. Nighttime amplitude scintillations, observed in the south of Guilin, always occurred with phase scintillations, TEC (Total Electron Content) depletions, and ROT (Rate Of change of TEC) fluctuations. However, TEC depletions and ROT fluctuations were weak during daytime amplitude scintillations, and daytime amplitude scintillations always took place simultaneously for most of the GPS satellites which appeared over Guilin in different azimuth directions. Ground-based GPS scintillation/TEC observations recorded at Guilin and signal-to-noise-ratio (SNR) measurements obtained from GPS-COSMIC radio occultation indicate that nighttime and daytime scintillations are very likely caused by ionospheric F region irregularities and sporadic E, respectively. Moreover, strong daytime amplitude scintillations may be associated with the plasma density enhancements in ionospheric E region caused by the Perseid and Geminid meteor shower activities.     

Evaluating the feasibility of GPS measurements of SSH on board a ship along the Portuguese West Coast

Aiming to validate coastal altimetry data, an exploratory experiment was carried out, using a survey ship with onboard GPS and motion compensation systems. The ship navigation plan followed a track as coincident as possible with the passage of two altimetry satellites, Envisat and Jason-1, along the Portuguese coast, in June 2007.     

Stochastic modelling considering ionospheric scintillation effects on GNSS relative and point positioning

Global Navigation Satellite Systems (GNSS), in particular the Global Positioning System (GPS), have been widely used for high accuracy geodetic positioning. The Least Squares functional models related to the GNSS observables have been more extensively studied than the corresponding stochastic models, given that the development of the latter is significantly more complex. As a result, a simplified stochastic model is often used in GNSS positioning, which assumes that all the GNSS observables are statistically independent and of the same quality, i.e. a similar variance is assigned indiscriminately to all of the measurements. However, the definition of the stochastic model may be approached from a more detailed perspective, considering specific effects affecting each observable individually, as for example the effects of ionospheric scintillation. These effects relate to phase and amplitude fluctuations in the satellites signals that occur due to diffraction on electron density irregularities in the ionosphere and are particularly relevant at equatorial and high latitude regions, especially during periods of high solar activity. As a consequence, degraded measurement quality and poorer positioning accuracy may result.     

An inter-comparison of zenith tropospheric delays derived from DORIS and GPS data

Doppler Orbitography Radiopositioning Integrated by Satellite (DORIS) and Global Positioning System (GPS) techniques are similarly affected by propagation delays in the neutral atmosphere (troposphere) and hence make use of similar data processing strategies for reducing this effect. We compare Zenith Tropospheric Delays (ZTDs) estimated from 52 DORIS and GPS station pairs co-located at 35 sites over the 2005–2008 period. We find an overall systematic negative mean bias of −4 mm and a median bias of −2 mm, with a large site-to-site scatter and especially stronger biases over South America, potentially linked to remaining problems related to the South Atlantic Anomaly (SAA) in the current DORIS data processing. The standard deviation of ZTD differences is in the range 4–12 mm over the globe (8 mm on average), with larger values located in the southern hemisphere. The spatial variability of differences is consistent with previous work but remains largely unexplained. DORIS is shown to be much less sensitive to instrumental changes than GPS (only the switch from Alcatel to Starec antenna at Toulouse is detected as an offset of −4 mm in the ZTD time series). On the opposite, discontinuities and spurious annual signals are found in the GPS ZTD solutions. A discontinuity of +5 mm is found on 5 November 2006, linked to the switch from relative to absolute GPS antenna models used in the data processing. The use of modified GPS antennas (e.g. at GODE) or improved antenna models is shown to reduce the spurious annual signal (e.g. from 5 mm to 2 mm at METS). Overall, the agreement between both techniques is good, though DORIS shows a significantly larger random scatter. The high stability and good spatial and temporal coverage make DORIS a potential candidate technique for meteorology and climate studies as long as reasonable time averaging can be applied (e.g. differences are reduced from 8.6 to 2.4 mm with 5-day averages) and no real-time application is considered. This technique could be considered as a potential contributor to Global Geodetic Observing System (GGOS) for climatology.     

Eliminating negative VTEC in global ionosphere maps using inequality-constrained least squares

Currently, ground-based Global Navigation Satellite System (GNSS) stations of the International GNSS Service (IGS) are distributed unevenly around the world. Most of them are located on the mainland, while only a small part of them are scattered on some islands in the oceans. As a consequence, many unreasonable zero values (in fact negative values) appear in Vertical Total Electron Content (VTEC) of European Space Agency (ESA) and Center for Orbit Determination in Europe (CODE) IONEX products, especially in 2008 and 2009 when the solar activities were rather quiet. To improve this situation, we directly implement non-negative physical constraints of ionosphere for global ionosphere maps (GIM) with spherical harmonic functions. Mathematically, we propose an inequality-constrained least squares method by imposing non-negative inequality constraints in the areas where negative VTEC values may occur to reconstruct GIM models. We then apply the new method to process the IGS data in 2008. The results have shown that the new algorithm efficiently eliminates the unwanted behavior of negative VTEC values, which could otherwise often be seen in the current CODE and ESA GIM products in both middle and high latitude areas of the Southern Hemisphere (45°S∼90°S) and the Northern Hemisphere (50°N∼90°N). About 64% of GPS receivers’ DCBs have been significantly improved. Finally, we compare the GIM results between with and without the inequality constraints, which has clearly shown that the GIM result with inequality constraints is significantly better than that without the inequality constraints. The inequality-constrained GIM result is also highly consistent with the final IGS products in terms of root mean squared (RMS) and mean VTEC.     

GPS phase difference variation statistics: A comparison between phase scintillation index and proxy indices

Scintillated GPS phase observations are traditionally characterized by the phase scintillation index, derived from specialized GPS receivers usually tracking at 50 Hz. Geodetic quality GPS receivers, on the other hand, are normally tracking at frequencies up to 1 Hz. However, availability of continuously operating geodetic receivers both in time and geographical location are superior to scintillation receiver’s coverage in many parts of the world. This motivates scintillation studies using regional and global geodetic GPS networks. Previous studies have shown the usefulness of GPS estimated total electron content variations for detecting ionospheric irregularities. In this paper, collocated geodetic and scintillation receivers are employed to compare proxy indices derived from geodetic receivers with the phase scintillation index during quiet and moderately disturbed ionospheric conditions. Sensitivity of the phase scintillation indices at high latitude stations to geomagnetic activity is discussed. Global mapping of ionospheric disturbances using proxy indices from real-time 1 Hz GPS stations are also presented.     

A global survey of COSMIC ionospheric peak electron density and its height: A comparison with ground-based ionosonde measurements

With a network of ground-based ionosondes distributed around the world, the ionospheric peak electron density and its height measured by FORMOSAT-3/COSMIC satellites in terms of GPS radio occultation technique are extensively examined in this article. It is found that, in spite of the latitude, the mean values of the peak electron density measured by COSMIC satellites are systematically smaller than those observed by ground-based ionosondes. The discrepancy between them is dependent on the latitude, namely, it is small in low and mid-latitudes and large in high-latitude region. Moreover, statistical analysis shows that the slopes of the regression line that is best fitted to the scatter diagram of occultation-retrieved peak electron density (ordinate axis) versus ionosonde-observed peak density (abscissa axis) are universally less than one. This feature is believed to be the result of path average effect of non-uniform distribution of the electron density along the GSP ray during the occultation. A comparison between COSMIC-measured peak height and ionosonde-derived peak height hmF2 indicates that the former is systematically higher than the latter. The difference in the two can be as large as 20% or more in equatorial and low-latitude regions. This result implies that the peak height hmF2 derived from the virtual height through true height analysis based on Titheridge method seems to underestimate the true peak height. The correlation between COSMIC and ionosonde peak electron densities is analyzed and the result reveals that correlation coefficient seems to be dependent on the fluctuation of the occultation-retrieved electron density profile. The correlation will be higher (lower) for the electron density profiles with smaller (larger) fluctuations. This feature suggests that the inhomogeneous distribution of the electron density along the GPS ray path during the occultation plays an important role affecting the correlation between COSMIC and ionosonde measurements.     

Response of the EIA ionosphere to the 7-8 May 2005 geomagnetic storm

In this paper, response of low latitude ionosphere to a moderate geomagnetic storm of 7–8 May 2005 (SSC: 1920 UT on 7 May with Sym-H minimum, ∼−112 nT around 1600 UT on 8 May) has been investigated using the GPS measurements from a near EIA crest region, Rajkot (Geog. 22.29°N, 70.74°E, Geomag.14°), India. We found a decrease in total electron content (TEC) in 12 h after the onset of the storm, an increase during and after 6 h of Sym-H deep minimum with a decrease below its usual-day level on the second day during the recovery phase of the storm. On 8 May, an increase of TEC is observed after sunset and during post-midnight hours (maximum up to 170%) with the formation of ionospheric plasma bubbles followed by a nearly simultaneous onset of scintillations at L-band frequencies following the time of rapid decrease in Sym-H index (−30 nT/h around 1300 UT).     

Dynamical evolution of high area-to-mass ratio debris released into GPS orbits

A large set of simulations, including all the relevant perturbations, was carried out to investigate the long-term dynamical evolution of fictitious high area-to-mass ratio (A/M) objects released, with a negligible velocity variation, in each of the six orbital planes used by Global Positioning System (GPS) satellites. As with similar objects discovered in near synchronous trajectories, long lifetime orbits, with mean motions of about 2 revolutions per day, were found possible for debris characterized by extremely high area-to-mass ratios. Often the lifetime exceeds 100 years up to A/M ∼ 45 m²/kg, decreasing rapidly to a few months above such a threshold. However, the details of the evolution, which are conditioned by the complex interplay of solar radiation pressure and geopotential plus luni-solar resonances, depend on the initial conditions. Different behaviors are thus possible. In any case, objects like those discovered in synchronous orbits, with A/M as high as 20–40 m²/kg, could also survive in this orbital regime, with semi-major axes close to the semi-synchronous values, with maximum eccentricities between 0.3 and 0.7, and with significant orbit pole precessions (faster and wider for increasing values of A/M), leading to inclinations between 30° and more than 90°.