Making the connection between Geosat and TOPEX/Poseidon

We can presently construct two independent time series of sea level, each at a precision of a few centimeters, from Geosat (1985–1988) and TOPEX/Poseidon (1992–1995) collinear altimetry. Both are based on precise satellite orbits computed using a common geopotential model, JGM-2 (Nerem et al. 1994). We have attempted to connect these series using Geosat-T/P crossover differences in order to assess long-term ocean changes between these missions. Unfortunately, the observed result are large-scale sea level differences which appear to be due to a combination of geodetic and geopotential error sources. The most significant geodetic component seems to be a coordinate system bias for Geosat sea level (relative to T/P) of −7 to −9 cm in the y-axis (towards the Eastern Pacific). The Geosat-T/P sea height differences at crossovers (with JGM-2 orbits) probably also contain stationary geopotential-orbit error of about the same magnitude which also distort any oceanographic interpretation of the observed changes. We also found JGM-3 Geosat orbits have not resolved the datum errors evident from the JGM-2 Geosat -T/P results. We conclude that the direct altimetric approach to accurate determination of sea level change using Geosat and T/P data still depends on further improvement in the Geosat orbits, including definition of the geocenter. Received: 11 March 1996; Accepted: 19 September 1996     

Precise ERS-1 orbit calculation

This analysis was performed with the GEOSAT software developed at NDRE for high-precision analysis of satellite tracking and VLBI data for geodetic and geodynamic applications.For applications to ERS-1, a realistic surface force model is used together with the Jacchia 77 atmospheric model, semi-daily drag coefficients, a 1-cpr sinusoidal along-track acceleration, and the GSFC JGM-2 gravity model. ERS-1 orbits have been derived for 5.5-day arcs of laser tracking data between July 6 and August 12, 1992. Results from overlapping orbits and comparison with precise D-PAF orbits indicate an orbital accuracy of 10–15 cm in the radial direction, ~ 60 cm in the along-track direction and ~ 15 cm in the cross-track direction.     

A spectral optimally interpolated mean sea surface from satellite radar altimetry

 A technique is presented for the development of a high-precision and high-resolution mean sea surface model utilising radar altimetric sea surface heights extracted from the geodetic phase of the European Space Agency (ESA) ERS-1 mission. The methodology uses a cubic-spline fit of dual ERS-1 and TOPEX crossovers for the minimisation of radial orbit error. Fourier domain processing techniques are used for spectral optimal interpolation of the mean sea surface in order to reduce residual errors within the initial model. The EGM96 gravity field and sea surface topography models are used as reference fields as part of the determination of spectral components required for the optimal interpolation algorithm. A comparison between the final model and 10 cycles of TOPEX sea surface heights shows differences of between 12.3 and 13.8 cm root mean square (RMS). An un-optimally interpolated surface comparison with TOPEX data gave differences of between 15.7 and 16.2 cm RMS. The methodology results in an approximately 10-cm improvement in accuracy. Further improvement will be attained with the inclusion of stacked altimetry from both current and future missions. Received: 22 December 1999 / Accepted: 6 November 2000     

The ERS-2 altimetric bias and gravity field enhancement using dual crossovers between ERS and TOPEX/Poseidon

 Dual satellite crossovers (DXO) between the two European Remote Sensing satellites ERS-1 and ERS-2 and TOPEX/Poseidon are used to (1) refine the Earth's gravity field and (2) extend the study of the ERS-2 altimetric range stability to cover the first four years of its operation. The enhanced gravity field model, AGM-98, is validated by several methodologies and will be shown to provide, in particular, low geographically correlated orbital error for ERS-2. For the ERS-2 altimetric range study, TOPEX/Poseidon is first calibrated through comparison against in situ tide gauge data. A time series of the ERS-2 altimeter bias has been recovered along with other geophysical correction terms using tables for bias jumps in the range measurements at the single point target response (SPTR) events. On utilising the original version of the SPTR tables the overall bias drift is seen to be 2.6±1.0 mm/yr with an RMS of fit of 12.2 mm but with discontinuities at the centimetre level at the SPTR events. On utilising the recently released revised tables, SPTR2000, the drift is better defined at 2.4±0.6 mm/yr with the RMS of fit reduced to 3.7 mm. Investigations identify the sea-state bias as a source of error with corrections affecting the overall drift by close to 1.2 mm/yr. Received: 25 May 2000 / Accepted: 24 January 2001     

Precise ERS-2 orbit determination using SLR, PRARE, and RA observations

The second European Remote Sensing satellite ERS-2, launched in 1995, is supported by satellite laser ranging (SLR), the German precise range and range-rate equipment (PRARE) microwave tracking system, and radar altimetry (RA). SLR, PRARE and crossover RA observations (XO/RA) from May, June, and July 1996 have been analyzed using the GEOSAT software developed at Forsvarets forskningsinstitutt (FFI, The Norwegian defence research establishment). Orbits computed with the JGM3 gravity model and the more recent EGM96 gravity model have been compared with orbits from other analysis centers. Based on these comparisons in addition to the post-fit observation residuals, and results from internal orbital overlaps, we conclude that the radial component of the ERS-2 orbit can be determined with a precision of 6 cm. Received: 30 June 1997 / Accepted: 4 February 1998     

TOPEX/Poseidon orbit error assessment

This paper discusses the accuracy of TOPEX/Poseidon orbits computed at Delft University, Section Space Research & Technology (DUT/SSR&T), from several types of tracking data,i.e. SLR, DORIS, and GPS. To quantify the orbit error, three schemes are presented. The first scheme relies on the direct altimeter observations and the covariance of the JGM-2 gravity field. The second scheme is based on crossover difference residuals while the third scheme uses the differences of dynamic orbit solutions with the GPS reduced-dynamic orbit. All three schemes give comparable results and indicate that the radial orbit error of TOPEX/Poseidon is 3–4 cm. From the orbit comparisons with GPS reduced dynamic, both the along-track and cross-track errors of the dynamic orbit solutions were found to be within 10–15 cm.     

Global mean sea surface and marine gravity anomaly from multi-satellite altimetry: applications of deflection-geoid and inverse Vening Meinesz formulae

 Global mean sea surface heights (SSHs) and gravity anomalies on a 2^′×2^′ grid were determined from Seasat, Geosat (Exact Repeat Mission and Geodetic Mission), ERS-1 (1.5-year mean of 35-day, and GM), TOPEX/POSEIDON (T/P) (5.6-year mean) and ERS-2 (2-year mean) altimeter data over the region 0^∘–360^∘ longitude and –80^∘–80^∘ latitude. To reduce ocean variabilities and data noises, SSHs from non-repeat missions were filtered by Gaussian filters of various wavelengths. A Levitus oceanic dynamic topography was subtracted from the altimeter-derived SSHs, and the resulting heights were used to compute along-track deflection of the vertical (DOV). Geoidal heights and gravity anomalies were then computed from DOV using the deflection-geoid and inverse Vening Meinesz formulae. The Levitus oceanic dynamic topography was added back to the geoidal heights to obtain a preliminary sea surface grid. The difference between the T/P mean sea surface and the preliminary sea surface was computed on a grid by a minimum curvature method and then was added to the preliminary grid. The comparison of the NCTU01 mean sea surface height (MSSH) with the T/P and the ERS-1 MSSH result in overall root-mean-square (RMS) differences of 5.0 and 3.1 cm in SSH, respectively, and 7.1 and 3.2 μrad in SSH gradient, respectively. The RMS differences between the predicted and shipborne gravity anomalies range from 3.0 to 13.4 mGal in 12 areas of the world's oceans. Received: 26 September 2001 / Accepted: 3 April 2002 Correspondence to: C. Hwang Acknowledgements. This research is partly supported by the National Science Council of ROC, under grants NSC89-2611-M-009-003-OP2 and NSC89-2211-E-009-095. This is a contribution to the IAG Special Study Group 3.186. The Geosat and ERS1/2 data are from NOAA and CERSAT/France, respectively. The T/P data were provided by AVISO. The CLS and GSFC00 MSS models were kindly provided by NASA/GSFC and CLS, respectively. Drs. Levitus, Monterey, and Boyer are thanked for providing the SST model. Dr. T. Gruber and two anonymous reviewers provided very detailed reviews that improved the quality of this paper.     

Ocean tides from T/P, ERS-1, and GEOSAT altimetry

 Aliasing of the diurnal and semi-diurnal tides is a major problem when estimating the ocean tides from satellite altimetry. As a result of aliasing, the tides become correlated and many years of altimeter observations may be needed to seperate them. For the three major satellite altimetry missions to date i.e., GEOSAT, ERS-1, and TOPEX/POSEIDON (T/P), the alias periods as well as the Rayleigh periods over which the tides decorrelate can be identified. Especially in case of GEOSAT and ERS-1, severe correlation problems arise. However, it is shown by means of covariance analyses that the tidal phase advance differences on crossing satellite groundtracks can significantly reduce the correlations among the diurnal and semi-diurnal tides and among these tides and the seasonal cycles of ocean variability. Therefore, it has been attempted to solve a multi-satellite response tidal solution for the diurnal and semi-diurnal bands from a total of 7 years of altimetry. Unfortunately, it could be shown that the GEOSAT and ERS-1 orbit errors are too large to improve a 3-year T/P tidal solution with about 2 years of GEOSAT and 2 years of ERS-1 altimeter observations. However, these results are preliminary and it is expected that more accurate orbits, which have become available recently for ERS-1, and additional altimeter data from ERS-2 and the GEOSAT Follow-On (GFO) should lead to an improved T/P tidal model. Received: 4 May 1999 / Accepted: 24 January 2000     

Spherical panel clustering and its numerical aspects

In modern approximation methods, linear combinations in terms of (space localizing) radial basis functions play an essential role. Areas of application are numerical integration formulae on the unit sphere Ω corresponding to prescribed nodes, spherical spline interpolation and spherical wavelet approximation. The evaluation of such a linear combination is a time-consuming task, since a certain number of summations, multiplications and the calculation of scalar products are required. A generalization of the panel clustering method in a spherical setup is presented. The economy and efficiency of panel clustering are demonstrated for three fields of interest, namely upward continuation of the Earth's gravitational potential, geoid computation by spherical splines and wavelet reconstruction of the gravitational potential. Received: 1 October 1997 / Accepted: 1 April 1998     

Spectral accuracy of Jgm3 from satellite crossover altimetry

A detailed accuracy assessment of the geopotential model Jgm3 is made based on independent single- and dual-satellite sea-height differences at crossovers from altimetry with Jgm3-based orbits. These differences, averaged over long time spans and in latitude bands, are converted to spectra (latitude-lumped coefficients) by least-squares estimation. The observed error spectra so obtained are then compared directly to error projections for them from the Jgm3 variance–covariance matrix. It is found from these comparisons that Jgm3 is generally well calibrated with respect to the crossover altimetry of and between Geosat, TOPEX/Poseidon (T/P), and Ers 1. Some significant discrepancies at a few lower orders (namely m=1 and 3) indicate a need for further improvement of Jgm3. A companion calibration (by order) of the geopotential model Jgm2 shows its variance–covariance matrix also to be generally well calibrated for the same single- and dual-satellite altimeter data sets (but based on Jgm2 orbits), except that the error projections for Geosat are too pessimistic. The analysis of the dual-satellite crossovers reveals possible relative coordinate system offsets (particularly for Geosat with respect to T/P) which have been discussed previously. The long-term detailed seasonally averaged Geosat sea level with respect to T/P (covering 1985–1996) should be useful in gauging the relative change in sea level between different parts of the ocean over the single 4-year gap between these missions (1988–1992). Received: 16 February 1998 / Accepted: 25 November 1998     

A harmonic spline model for local estimation of planetary gravity fields from line-of-sight acceleration data

The spline interpolation technique is applied to estimate locally the radial component of a planetary gravity field from residual acceleration data along a special direction, the direction of observation from the Earth. After the presentation of the theoretical framework, the method is tested on synthetic and real data in the case of Venus. It is shown that this spline technique can be used succesfully to build local models of radial gravity fields at the planet surface. Received: 13 March 1997 / Accepted: 17 November 1998     

The Abel-Poisson kernel and the Abel-Poisson integral in a moving tangent space

The upward-downward continuation of a harmonic function like the gravitational potential is conventionally based on the direct-inverse Abel-Poisson integral with respect to a sphere of reference. Here we aim at an error estimation of the “planar approximation” of the Abel-Poisson kernel, which is often used due to its convolution form. Such a convolution form is a prerequisite to applying fast Fourier transformation techniques. By means of an oblique azimuthal map projection / projection onto the local tangent plane at an evaluation point of the reference sphere of type “equiareal” we arrive at a rigorous transformation of the Abel-Poisson kernel/Abel-Poisson integral in a convolution form. As soon as we expand the “equiareal” Abel-Poisson kernel/Abel-Poisson integral we gain the “planar approximation”. The differences between the exact Abel-Poisson kernel of type “equiareal” and the “planar approximation” are plotted and tabulated. Six configurations are studied in detail in order to document the error budget, which varies from 0.1% for points at a spherical height H=10km above the terrestrial reference sphere up to 98% for points at a spherical height H = 6.3×10⁶km. Received: 18 March 1997 / Accepted: 19 January 1998     

Calibration of satellite gradiometer data aided by ground gravity data

Parametric least squares collocation was used in order to study the detection of systematic errors of satellite gradiometer data. For this purpose, simulated data sets with a priori known systematic errors were produced using ground gravity data in the very smooth gravity field of the Canadian plains. Experiments carried out at different satellite altitudes showed that the recovery of bias parameters from the gradiometer “measurements” is possible with high accuracy, especially in the case of crossing tracks. The mean value of the differences (original minus estimated bias parameters) was relatively large compared to the standard deviation of the corresponding second-order derivative component at the corresponding height. This mean value almost vanished when gravity data at ground level were combined with the second-order derivative data set at satellite altitude. In the case of simultaneous estimation of bias and tilt parameters from ∂² T/∂z ²“measurements”, the recovery of both parameters agreed very well with the collocation error estimation. Received: 10 October 1996 / Accepted 25 May 1998     

The exterior Airy/Heiskanen topographic–isostatic gravity potential, anomaly and the effect of analytical continuation in Stokes' formula

This study deals with the external type of topographic–isostatic potential and gravity anomaly and its vertical derivatives, derived from the Airy/Heiskanen model for isostatic compensation. From the first and the second radial derivatives of the gravity anomaly the effect on the geoid is estimated for the downward continuation of gravity to sea level in the application of Stokes' formula. The major and regional effect is shown to be of order H ³ of the topography, and it is estimated to be negligible at sea level and modest for most mountains, but of the order of several metres for the highest and most extended mountain belts. Another, global, effect is of order H but much less significant Received: 3 October 1997 / Accepted: 30 June 1998     

Role of ocean variability and dynamic ocean topography in the recovery of the mean sea surface and gravity anomalies from satellite altimeter data

Procedures to calculate mean sea surface heights and gravity anomalies from altimeter-derived sea surface heights and along-track sea surface slopes using the least-squares collocation procedure are derived. The slope data is used when repeat track averaging is not possible to reduce ocean variability effects. Tests were carried out using Topex, Geosat, ERS-1 [35-day and 168-day (2 cycle)] data. Calculations of gravity anomalies in the Gulf Stream region were made using the sea surface height and slope data. Tests were also made correcting the sea surface heights for dynamic ocean topography calculated from a degree 360 expansion of data from the POCM-4B global ocean circulation model. Comparisons of the anomaly predictions were carried out with ship data using anomalies calculated for this paper as well as others. Received: 19 August 1996 / Accepted: 14 April 1997     

Unbiased vs biased estimation of GPS phase ambiguities from dual-frequency code and phase observables

With access to dual-frequency pseudorange and phase Global Positioning System (GPS) data, the wide-lane ambiguity can easily be fixed. Advantage is taken of this information in the linear combination of the above four observables for base ambiguity estimation (i.e. of N ₁ and N ₂). Starting points for our analysis are the Best Linear Unbiased Estimators BLUE₁ and BLUE₂. BLUE₁ is the best one (with minimum mean square error, MSE) if the ionosphere effect is negligible. If this is not the case, BLUE₂ has the smallest variance, but not necessarily the least mean square error. Hence, both estimators may suffer from a non-optimal treatment of the ionosphere bias. BLUE₁ ignores possible ionosphere bias, while BLUE₂ compensates for this bias in a less favourable way by eliminating it at the price of increased noise. As an alternative, linear estimators are derived, which make a compromise between the ionosphere bias and the random observation errors. This leads to the derivation of the Best Linear Estimator (BLE) and the Restricted Best Linear Estimator (RBLE) with minimum MSE. The former is generally not very useful, while the RBLE is recommended for practical use. It is shown that the MSE of the RBLE is limited by the variances of BLUE₁ and BLUE₂, i.e.

Theory of satellite ground-track crossovers

The fundamental geometry of satellite ground tracks and their crossover problem are investigated. For idealized nominal ground tracks, the geometry is governed by a few constant parameters whose variations lead to qualitative changes in the crossover solutions. On the basis that the theory to locate crossovers has not been studied in sufficient detail, such changes are described in regard to the number of crossover solutions in conjunction with their bifurcations. Employing the spinor algebra as a tool for establishing the ground-track crossing condition, numerical methodologies to locate crossovers appearing in general dual-satellite ground-track configurations are also presented. The methodologies are applied to precisely determined orbital ephemerides of the GEOSAT, ERS-1, and TOPEX/POSEIDON altimeter satellites. Received: 19 November 1996 / Accepted: 12 May 1997     

On calculating the Earth's C21 and S21 gravity coefficients in the IERS terrestrial reference frame

Modern models of the Earth's gravity field are developed in the IERS (International Earth Rotation Service) terrestrial reference frame. In this frame the mean values for gravity coefficients of the second degree and first order, C ₂₁(IERS) and S ₂₁(IERS), by the current IERS Conventions are recommended to be calculated by using the observed polar motion parameters. Here, it is proved that the formulae presently employed by the IERS Conventions to obtain these coefficients are insufficient to ensure their values as given by the same source. The relevant error of the normalized mean values for C ₂₁(IERS) and S ₂₁(IERS) is 3×10⁻¹², far above the adopted cutoff (10⁻¹³) for variations of these coefficients. Such an error in C ₂₁ and S ₂₁ can produce non-modeled perturbations in motion prediction of certain artificial Earth satellites of a magnitude comparable to the accuracy of current tracking measurements. Received: 14 September 1998 / Accepted: 20 May 1999     

Comment and reply regarding Changyou Zhang (1995) A general formula and its inverse formula for gravimetric transformation by use of convolution and deconvolution techniques.

The incorrectness in the derivations of error formula (52) in Zhang's paper is pointed out and the correction is given. Table 3 [list of error formulae of Eqs (5)–(8)] in Zhang's paper is also corrected. Received: 25 February 1998 / Accepted: 30 June 1998     

Truncation error formulae for the disturbing gravity vector

Recurrence relations have been derived for truncation error coefficients of the extended Stokes' function and its partial derivatives required in the computation of the disturbing gravity vector at any elevation above the earth's surface. The corresponding formulae, the example of values of the truncation error coefficients for H=30.1 km and ψ₀=30^∘ and the estimations of truncation error are given in this article. Received: 26 January 1996 / Accepted: 11 June 1997