Towards development of a consistent orbit series for TOPEX,Jason-1, and Jason-2

The TOPEX/Poseidon, Jason-1 and Jason-2 set of altimeter data now provide a time series of synoptic observations of the ocean that span nearly 17 years from the launch of TOPEX in 1992. The analysis of the altimeter data including the use of altimetry to monitor the global change in mean sea level requires a stable, accurate, and consistent orbit reference over the entire time span. In this paper, we describe the recomputation of a time series of orbits that rely on a consistent set of reference frames and geophysical models. The recomputed orbits adhere to the IERS 2003 standards for ocean and earth tides, use updates to the ITRF2005 reference frame for both the SLR and DORIS stations, apply GRACE-derived models for modeling of the static and time-variable gravity, implement the University College London (UCL) radiation pressure model for Jason-1, use improved troposphere modeling for the DORIS data, and apply the GOT4.7 ocean tide model for both dynamical ocean tide modeling and for ocean loading. The new TOPEX orbits have a mean SLR fit of 1.79 cm compared to 2.21 cm for the MGDR-B orbits. These new TOPEX orbits agree radially with independent SLR/crossover orbits at 0.70 cm RMS, and the orbit accuracy is estimated at 1.5–2.0 cm RMS over the entire TOPEX time series. The recomputed Jason-1 orbits agree radially with the Jason-1 GDR-C orbits at 1.08 cm RMS. The GSFC SLR/DORIS dynamic and reduced-dynamic orbits for Jason-2 agree radially with independent orbits from the CNES and JPL at 0.70–1.06 cm RMS. Applying these new orbits, and using the latest altimeter corrections for TOPEX, Jason-1, and Jason-2 from September 1992 to May 2009, we find a global rate in mean sea level of 3.0 ± 0.4 mm/yr.     

Regional in situ validation of satellite altimeters: Calibration and cross-calibration results at the Corsican sites

The in situ validation of the satellite altimeter sea surface heights is generally performed either at a few local points directly flown over by the satellites or using the global tide gauge network. A regional in situ calibration method was developed by NOVELTIS in order to monitor the altimeter data quality in a perimeter of several hundred kilometres around a given in situ calibration site. The primary advantage of this technique is its applicability not only for missions flying over dedicated sites but also for missions on interleaved or non repetitive orbits. This article presents the altimeter bias estimates obtained with this method at the Corsican calibration site, for the Jason-1 mission on its nominal and interleaved orbits as well as for the Jason-2 and Envisat missions. The various regional bias estimates (8.2 cm and 7.4 cm for Jason-1 respectively on the nominal and interleaved orbits in Senetosa, 16.4 cm for Jason-2 in Senetosa and 47.0 cm for Envisat in Ajaccio, with an accuracy between 2.5 cm and 4 cm depending on the mission) are compared with the results obtained by the other in situ calibration teams. This comparison demonstrates the coherency at the centimetre level, the stability and the generic character of the method, which would also be of benefit to the new and future altimeter missions such as Cryosat-2, SARAL/AltiKa, Sentinel-3, Jason-3, Jason-CS.     

Generating precise and homogeneous orbits for Jason-1 and Jason-2

Driven by the GMES (Global Monitoring for Environment and Security) and GGOS (Global Geodetic Observing System) initiatives the user community has a strong demand for high-quality altimetry products. In order to derive such high-quality altimetry products, precise orbits for the altimetry satellites are a necessity. With the launch of the TOPEX/Poseidon mission in 1992 a still on-going time series of high-accuracy altimetry measurements of ocean topography started, continued by the altimetry missions Jason-1 in 2001 and Jason-2/OSTM in 2008. This paper contributes to the on-going orbit reprocessing carried out by several groups and presents the efforts of the Navigation Support Office at ESA/ESOC using its NAPEOS software for the generation of precise and homogeneous orbits referring to the same reference frame for the altimetry satellites Jason-1 and Jason-2. Data of all three tracking instruments on-board the satellites (beside the altimeter), i.e. GPS, DORIS, and SLR measurements, were used in a combined data analysis. About 7 years of Jason-1 data and more than 1 year of Jason-2 data were processed. Our processing strategy is close to the GDR-C standards. However, we estimated slightly different scaling factors for the solar radiation pressure model of 0.96 and 0.98 for Jason-1 and Jason-2, respectively. We used 30 s sampled GPS data and introduced 30 s satellite clocks stemming from ESOC’s reprocessing of the combined GPS/GLONASS IGS solution. We present the orbit determination results, focusing on the benefits of adding GPS data to the solution. The fully combined solution was found to give the best orbit results. We reach a post-fit RMS of the GPS phase observation residuals of 6 mm for Jason-1 and 7 mm for Jason-2. The DORIS post-fit residuals clearly benefit from using GPS data in addition, as the DORIS data editing improves. The DORIS observation RMS for the fully combined solution is with 3.5 mm and 3.4 mm, respectively, 0.3 mm better than for the DORIS-SLR solution. Our orbit solution agrees well with external solutions from other analysis centers, as CNES, LCA, and JPL. The orbit differences between our fully combined orbits and the CNES GDR-C orbits are of about 0.8 cm for Jason-1 and at 0.9 cm for Jason-2 in the radial direction. In the cross-track component we observe a clear improvement when adding GPS data to the POD process. The 3D-RMS of the orbit differences reveals a good orbit consistency at 2.7 cm and 2.9 cm for Jason-1 and Jason-2. Our resulting orbit series for both Jason satellites refer to the ITRF2005 reference frame and are provided in sp3 file format on our ftp server.     

Using altimetry and seafloor pressure data to estimate vertical deformation offshore: Vanuatu case study

Measuring ground deformation underwater is essential for understanding Earth processes at many scales. One important example is subduction zones, which can generate devastating earthquakes and tsunamis, and where the most important deformation signal related to plate locking is usually offshore. We present an improved method for making offshore vertical deformation measurements, that involve combining tide gauge and altimetry data. We present data from two offshore sites located on either side of the plate interface at the New Hebrides subduction zone, where the Australian plate subducts beneath the North Fiji basin. These two sites have been equipped with pressure gauges since 1999, to extend an on-land GPS network across the plate interface. The pressure series measured at both sites show that Wusi Bank, located on the over-riding plate, subsides by 11 ± 4 mm/yr with respect to Sabine Bank, which is located on the down-going plate. By combining water depths derived from the on-bottom pressure data with sea surface heights derived from altimetry data, we determine variations of seafloor heights in a global reference frame. Using altimetry data from TOPEX/Poseidon, Jason-1, Jason-2 and Envisat missions, we find that the vertical motion at Sabine Bank is close to zero and that Wusi Bank subsides by at least 3 mm/yr and probably at most 11 mm/yr.     

Altimetry,bathymetry and geoid variations at the Gavdos permanent Cal/Val facility

The aim of this work has been to examine the relationship of steep bathymetry in the coastal areas around the permanent Cal/Val facility of Gavdos, and their influence on the produced calibration values for the Jason-2 satellite altimeter. The paper describes how changes in seafloor topography (from 200 to 3500 m depth over a distance of 10 km) are reflected on the determined altimeter parameters using different reference surfaces for satellite calibration. Finally, it describes the relation between these parameter trends and the region’s local characteristics.     

Vertical motions in Thailand after the 2004 Sumatra–Andaman Earthquake from GPS observations and its geophysical modelling

Following previous findings from ongoing GPS research in Thailand since 2004 we continue to exploit the GPS technique to monitor and model land motions induced by the Sumatra–Andaman Earthquake. Our latest results show that up to the end of 2010, Thailand has been co-seismically displaced and is subsequently undergoing a post-seismic horizontal deformation with total displacements (co-seismic plus post-seismic) ranging from 10.5 to 74.7 cm. We observed the largest horizontal displacements in the southern part of Thailand and moderate and small displacements in the central and northern parts. In addition to horizontal displacements throughout Thailand, continuous GPS measurements show that large parts of Thailand are subsiding at rates up to 1 cm/yr. It is the first time that such vertical post-seismic deformations at large distances (650–1500 km away from the Earthquake’s epicentre) have been recorded. We have investigated the physical processes leading to the observed subsidence. While after-slip on the subduction interface induces negligible or even slightly positive vertical motions, relaxation in the asthenosphere is associated with a sizable subsidence. Predictions from a 3D finite element model feature an asthenosphere with an effective viscosity of the order of 3 * 10¹⁸ Pas, fit the horizontal post-seismic data and the observed subsidence well. This model is then used to predict the subsidence over the whole seismic cycle. The subsidence should go on with a diminishing rate through the next two decades and its final magnitude should not exceed 10 cm in the Bangkok area.     

GPS-based sea level measurements to help the characterization of land contamination in coastal areas

The Corsica site has been established in 1996 to perform altimeter calibration on TOPEX/Poseidon and then on its successors Jason-1 and Jason-2. The first chosen location was under the #85 ground track that overflight the Senetosa Cape. In 2005, it was decided to develop another location close to Ajaccio, to be able to perform the calibration of Envisat and in a next future of SARAL/AltiKa that will flight over the same ground tracks. Equipped with various instruments (tide gauges, permanent GPS, GPS buoy, weather station…) the Corsica calibration site is able to quantify the altimeter Sea Surface Height bias but also to give an input on the origin of this bias (range, corrections, orbits, …). Due to the size of Corsica (not a tiny island), the altimeter measurement system (range and corrections) can be contaminated by land. The aim of this paper is to evaluate this land contamination by using GPS measurements from a fixed receiver on land and from another receiver onboard a life buoy. Concerning the altimeter land contamination, we have quantify that this effect can reach 8 mm/km and then affects the Sea Surface Height bias values already published in the framework of the Corsica calibration site by 5–8 mm for TOPEX and Jason missions. On the other hand, the radiometer measurements (wet troposphere correction) are also sensitive to land and we have been able to quantify the level of improvement of a dedicated coastal algorithm that reconciles our results with those coming from other calibration sites. Finally, we have also shown that the standard deviation of the GPS buoy sea level measurements is highly correlated (∼87%) with the Significant Wave Height derived from the altimeters and can be used to validate such parameter.     

Calibration of Envisat radar altimeter over Lake Issykkul

This study presents the results of calibration/validation (C/V) of Envisat satellite radar altimeter over Lake Issykkul located in Kyrgyzstan, which was chosen as a dedicated radar altimetry C/V site in 2004. The objectives are to estimate the absolute altimeter bias of Envisat and its orbit based on cross-over analysis with TOPEX/Poseidon (T/P), Jason-1 and Jason-2 over the ocean. We have used a new method of GPS data processing in a kinematic mode, developed at the Groupe de Recherche de Geodesie Spatiale (GRGS), which allows us to calculate the position of the GPS antenna without needing a GPS reference station. The C/V is conducted using various equipments: a local GPS network, a moving GPS antenna along the satellites tracks over Lake Issykkul, In Situ level gauges and weather stations. The absolute bias obtained for Envisat from field campaigns conducted in 2009 and 2010 is between 62.1 and 63.4 ± 3.7 cm, using the Ice-1 retracking algorithm, and between 46.9 and 51.2 cm with the ocean retracking algorithm. These results differ by about 10 cm from previous studies, principally due to improvement of the C/V procedure. Apart from the new algorithm for GPS data processing and the orbit error reduction, more attention has been paid to the GPS antenna height calculation, and we have reduced the errors induced by seiche over Lake Issykkul. This has been assured using cruise data along the Envisat satellite track at the exact date of the pass of the satellite for the two campaigns. The calculation of the Envisat radar altimeter bias with respect to the GPS levelling is essential to allow the continuity of multi-mission data on the same orbit, with the expected launch of SARAL/Altika mission in 2012. Implications for hydrology in particular, will be to produce long term homogeneous and reliable time series of lake levels worldwide.     

The effect of geocenter motion on Jason-2 orbits and the mean sea level

We compute a series of Jason-2 GPS and SLR/DORIS-based orbits using ITRF2005 and the std0905 standards ( Lemoine et al., 2010). Our GPS and SLR/DORIS orbit data sets span a period of 2 years from cycle 3 (July 2008) to cycle 74 (July 2010). We extract the Jason-2 orbit frame translational parameters per cycle by the means of a Helmert transformation between a set of reference orbits and a set of test orbits. We compare the annual terms of these time-series to the annual terms of two different geocenter motion models where biases and trends have been removed. Subsequently, we include the annual terms of the modeled geocenter motion as a degree-1 loading displacement correction to the GPS and SLR/DORIS tracking network of the POD process. Although the annual geocenter motion correction would reflect a stationary signal in time, under ideal conditions, the whole geocenter motion is a non-stationary process that includes secular trends. Our results suggest that our GSFC Jason-2 GPS-based orbits are closely tied to the center of mass (CM) of the Earth consistent with our current force modeling, whereas GSFC’s SLR/DORIS-based orbits are tied to the origin of ITRF2005, which is the center of figure (CF) for sub-secular scales. We quantify the GPS and SLR/DORIS orbit centering and how this impacts the orbit radial error over the globe, which is assimilated into mean sea level (MSL) error, from the omission of the annual term of the geocenter correction. We find that for the SLR/DORIS std0905 orbits, currently used by the oceanographic community, only the negligence of the annual term of the geocenter motion correction results in a – 4.67 ± 3.40 mm error in the Z-component of the orbit frame which creates 1.06 ± 2.66 mm of systematic error in the MSL estimates, mainly due to the uneven distribution of the oceans between the North and South hemisphere.     

Dayside ionospheric uplift during strong geomagnetic storms as detected by the CHAMP,SAC-C,TOPEX and Jason-1 satellites

Ionosphere response to severe geomagnetic storms that occurred in 2001–2003 was analyzed using data of global ionosphere maps (GIM), altimeter data from the Jason-1 and TOPEX satellites, and data of GPS receivers on-board CHAMP and SAC-C satellites. This allowed us to study in detail ionosphere redistribution due to geomagnetic storms, dayside ionospheric uplift and overall dayside TEC increase. It is shown that after the interplanetary magnetic field turns southward and intensifies, the crests of the equatorial ionization anomaly (EIA) travel poleward and the TEC value within the EIA area increases significantly (up to ∼50%). GPS data from the SAC-C satellite show that during the main phase of geomagnetic storms TEC values above the altitude of 715 km are 2–3 times higher than during undisturbed conditions. These effects of dayside ionospheric uplift occur owing to the “super-fountain effect” and last few hours while the enhanced interplanetary electric field impinged on the magnetopause.     

The impact of microwave absorber and radome geometries on GNSS measurements of station coordinates and atmospheric water vapour

We have used microwave absorbing material in different geometries around ground-based Global Navigation Satellite System (GNSS) antennas in order to mitigate multipath effects on the estimates of station coordinates and atmospheric water vapour. The influence of a hemispheric radome – of the same type as in the Swedish GPS network SWEPOS – was also investigated. Two GNSS stations at the Onsala Space Observatory were used forming a 12 m baseline. GPS data from October 2008 to November 2009 were analyzed by the GIPSY/OASIS II software using the Precise Point Positioning (PPP) processing strategy for five different elevation cutoff angles from 5° to 25°. We found that the use of the absorbing material decreases the offset in the estimated vertical component of the baseline from ∼27 mm to ∼4 mm when the elevation cutoff angle varies from 5° to 20°. The horizontal components are much less affected. The corresponding offset in the estimates of the atmospheric Integrated Water Vapour (IWV) decreases from ∼1.6 kg/m² to ∼0.3 kg/m². Changes less than 5 mm in the offsets in the vertical component of the baseline are seen for all five elevation cutoff angle solutions when the antenna was covered by a hemispheric radome. Using the radome affects the IWV estimates less than 0.4 kg/m² for all different solutions. IWV comparisons between a Water Vapour Radiometer (WVR) and the GPS data give consistent results.     

The contribution of local gravimetric geoid models to the calibration of satellite altimetry data and an outlook of the latest GOCE GGM performance in Gavdos

The use of geoid heights has been one of the available methodologies utilized for the independent calibration/validation of altimeters on-board satellites. This methodology has been employed for long in the Gavdos dedicated cal/val facility (Crete, Greece), where calibration results for the Jason satellites have been estimated, both for ascending and descending passes. The present work gives a detailed overview of the methodology followed in order to estimate a high-resolution and accuracy gravimetric geoid model for the wider Gavdos area, in support of the on-going calibration work. To estimate the geoid model, the well-known remove-compute-restore method is used while residual geoid heights are estimated through least-squares collocation so that associated errors are determined as well. It is found that the estimated formal geoid errors from LSC along passes 018 and 109 of Jason satellites, used for the bias estimation, range between ±0.8–1.6 cm. The so-derived geoid heights are employed in the determination of the Jason-2 altimeter bias for all available cycles (cycles 1-114, spanning the period from July 2008 to August 2011) together with the RioMed DOT model. From the results acquired the Jason-2 bias has been estimated to be +196.1 ± 3.2 mm for pass 109 and +161.9 ± 5.1 mm for pass 018. Within the same frame, the GOCE/GRACE-based geopotential model GOCO02s has been used to estimate the mean dynamic ocean topography and the steady-state circulation in the area around Gavdos. The so-derived DOT model was used to estimate the Jason-2 bias in an effort to evaluate the performance of satellite-only geoid models and investigate whether their spatial resolution and accuracy provides some improvement w.r.t. traditional local gravimetric geoids. From the results acquired with geoid heights from GOCO02s, the estimated Jason-2 bias deviates significantly from that of the local gravimetric model, which can be attributed to a possible mean offset and the low resolution of GOCE-based GGMs. On the other hand, when the newly estimated GOCE-based DOT was employed with geoid heights from the local gravimetric geoid model, the Jason-2 bias has been estimated to be +185.1 ± 3.2 mm for pass 109 and +130.2 ± 5.1 mm for pass 018.     

Analysis of electron content variations over Japan during solar minimum: Observations and modeling

Comparative analysis of GPS TEC data and FORMOSAT-3/COSMIC radio occultation measurements was carried out for Japan region during period of the extremely prolonged solar minimum of cycle 23/24. COSMIC data for different seasons corresponded to equinox and solstices of the years 2007–2009 were analyzed. All selected electron density profiles were integrated up to the height of 700 km (altitude of COSMIC satellites), the monthly median estimates of Ionospheric Electron Content (IEC) were retrieved with use of spherical harmonics expansion. Monthly medians of TEC values were calculated from diurnal variations of GPS TEC estimates during considered month. Joint analysis of GPS TEC and COSMIC data allows us to extract and estimate electron content corresponded to the ionosphere (its bottom and topside parts) and the plasmasphere (h > 700 km) for different seasons of 2007–2009. Percentage contribution of ECpl to GPS TEC indicates the clear dependence from the time and varies from a minimum of about 25–50% during day-time to the value of 50–75% at night-time. Contribution of both bottom-side and topside IEC has minimal values during winter season in compare with summer season (for both day- and night-time). On average bottom-side IEC contributes about 5–10% of GPS TEC during night and about 20–27% during day-time. Topside IEC contributes about 15–20% of GPS TEC during night and about 35–40% during day-time. The obtained results were compared with TEC, IEC and ECpl estimates retrieved by Standard Plasmasphere–Ionosphere Model that has the plasmasphere extension up to 20,000 km (GPS orbit).     

Measuring atmospheric density using GPS–LEO tracking data

We present a method to estimate the total neutral atmospheric density from precise orbit determination of Low Earth Orbit (LEO) satellites. We derive the total atmospheric density by determining the drag force acting on the LEOs through centimeter-level reduced-dynamic precise orbit determination (POD) using onboard Global Positioning System (GPS) tracking data. The precision of the estimated drag accelerations is assessed using various metrics, including differences between estimated along-track accelerations from consecutive 30-h POD solutions which overlap by 6 h, comparison of the resulting accelerations with accelerometer measurements, and comparison against an existing atmospheric density model, DTM-2000. We apply the method to GPS tracking data from CHAMP, GRACE, SAC-C, Jason-2, TerraSAR-X and COSMIC satellites, spanning 12 years (2001–2012) and covering orbital heights from 400 km to 1300 km. Errors in the estimates, including those introduced by deficiencies in other modeled forces (such as solar radiation pressure and Earth radiation pressure), are evaluated and the signal and noise levels for each satellite are analyzed. The estimated density data from CHAMP, GRACE, SAC-C and TerraSAR-X are identified as having high signal and low noise levels. These data all have high correlations with anominal atmospheric density model and show common features in relative residuals with respect to the nominal model in related parameter space. On the contrary, the estimated density data from COSMIC and Jason-2 show errors larger than the actual signal at corresponding altitudes thus having little practical value for this study. The results demonstrate that this method is applicable to data from a variety of missions and can provide useful total neutral density measurements for atmospheric study up to altitude as high as 715 km, with precision and resolution between those derived from traditional special orbital perturbation analysis and those obtained from onboard accelerometers.     

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.     

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.     

Detection of Envisat RA2/ICE-1 retracked radar altimetry bias over the Amazon basin rivers using GPS

Altimetry is now routinely used to monitor stage variations over rivers, including in the Amazon basin. It is desirable for hydrologic studies to be able to combine altimetry from different satellite missions with other hydrogeodesy datasets such as leveled gauges and watershed topography. One requirement is to accurately determine altimetry bias, which could be different for river studies from the altimetry calibrated for deep ocean or lake applications. In this study, we estimate the bias in the Envisat ranges derived from the ICE-1 waveform retracking, which are nowadays widely used in hydrologic applications. As a reference, we use an extensive dataset of altitudes of gauge zeros measured by GPS collocated at the gauges. The thirty-nine gauges are spread along the major tributaries of the Amazon basin. The methodology consists in jointly modeling the vertical bias and spatial and temporal slope variations between altimetry series located upstream and downstream of each gauge. The resulting bias of the Envisat ICE-1 retracked altimetry over rivers is 1.044 ± 0.212 m, revealing a significant departure from other Envisat calibrations or from the Jason-2 ICE-1 calibration.     

Ambiguity resolution in precise point positioning with hourly data for global single receiver

Integer ambiguity resolution (IAR) can improve precise point positioning (PPP) performance significantly. IAR for PPP became a highlight topic in global positioning system (GPS) community in recent years. More and more researchers focus on this issue. Progress has been made in the latest years. In this paper, we aim at investigating and demonstrating the performance of a global zero-differenced (ZD) PPP IAR service for GPS users by providing routine ZD uncalibrated fractional offsets (UFOs) for wide-lane and narrow-lane. Data sets from all IGS stations collected on DOY 1, 100, 200 and 300 of 2010 are used to validate and demonstrate this global service. Static experiment results show that an accuracy better than 1 cm in horizontal and 1–2 cm in vertical could be achieved in ambiguity-fixed PPP solution with only hourly data. Compared with PPP float solution, an average improvement reaches 58.2% in east, 28.3% in north and 23.8% in vertical for all tested stations. Results of kinematic experiments show that the RMS of kinematic PPP solutions can be improved from 21.6, 16.6 and 37.7 mm to 12.2, 13.3 and 34.3 mm for the fixed solutions in the east, north and vertical components, respectively. Both static and kinematic experiments show that wide-lane and narrow-lane UFO products of all satellites can be generated and provided in a routine way accompanying satellite orbit and clock products for the PPP user anywhere around the world, to obtain accurate and reliable ambiguity-fixed PPP solutions.     

Cross testing of ionosphere models IRI-2001 and IRI-2007, data from satellite altimeters (Topex/Poseidon and Jason-1) and global ionosphere maps

In this study, we use a great body of statistical data covering the entire 23rd solar cycle to cross test data of satellite altimeters, Global Ionosphere Maps and the International Reference Ionosphere models, IRI-2001 and IRI-2007. It is revealed that experimental TEC values of the satellite altimeters regularly exceed the model ones by ∼3 TECU (1 TECU = 10¹⁶ m⁻²). The best possible value of difference between TECs obtained from altimeter and GIM-map data significantly differs for different laboratories: the maximum for CODG data falls on 2.5 TECU, ESAG – 3 TECU, JPLG – 0 TECU, UPCG – 2 TECU. The dependence of experimental and model data root-mean-square deviation on the F10.7 index is shown to be nearly linear. IRI-2001 and IRI-2007 relative errors are characterized by considerable 11-year and annual variations. Given the geomagnetic planetary index K_p under 7, IRI-2001 and IRI-2007 reproduce TEC in the ionosphere with an accuracy of ∼30% relative to measurement data from satellite altimeters. The amplitude of absolute error variations resulting from the difference in ionization enhancement between the model and the real ionosphere during the morning solar terminator transit is ∼5 TECU.