Dual-satellite crossover latitude-lumped coefficients, their use in geodesy and oceanography

Latitude-lumped coefficients (LLC) are defined, representing geopotential-orbit variations for dual-satellite crossovers (DSC). Formulae are derived for their standard errors from the covariances of geopotential field models. Numerical examples are presented for pairs of the altimeter-bearing satellites TOPEX/Poseidon, ERS 1, and Geosat, using the error matrices of recent gravity models. The DSC, connecting separate missions, will play an increasingly important role in oceanography spanning decades only when its nonoceanographic signals are thoroughly understood. In general, the content of even the long-term averaged DSC is more complex then their single satellite crossover (SSC) counterpart. The LLC, as the spatial spectra for the geopotential-caused crossover effects, discriminate these source-differences sharply. Thus, the zero-order LLC in DSC data contains zonal gravity information not present in SSC data. In addition, zero- and first-order LLC of DSC data can reveal a geocenter discrepancy between the orbit tracking of the separate satellite missions. For example, DSC analysis from orbits computed with JGM 2 show that the y-axis of the geocenter for Geosat in 1986–1988 is shifted with respect to T/P by 6–9 cm towards the eastern Pacific. Also, where the time-gap is necessarily large (as between, say, Geosat and T/P missions) oceanographic (sea-level) differences in DSC may corrupt the geopotential interpretation of the data. Most importantly, as we illustrate, media delays for the altimeter (from the ionosphere, wet troposphere and sea-state bias) are more likely sources of contamination across two missions than in SSC analyses. Again, the LLC of zero order best shows this contrast. Using the higher-order LLC of DSC for both Geosat-T/P and ERS 1-T/P as likely representation of geopotential-only error, we show by comparison with the predicted standard errors of JGM 2 that the latter's previously calibrated covariance matrix is generally valid. Received: 14 February 1996 / Accepted: 27 March 1997     

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     

Seasonal sea level change from TOPEX/Poseidon observation and thermal contribution

Seasonal steric sea-level change due to temperature variation in the mixing layer is assessed using space-measured sea-surface temperature data and historical in situ temperature measurements. The results are compared with TOPEX/Poseidon satellite altimeter measurement at different large spatial scales. It is indicated that thermal effect accounts for much of the observed seasonal variability, especially when averaging over zonal regions. Some regional seasonal patterns of sea-level anomalies in the tropical oceans are well represented by the thermal model prediction. Systematic differences are shown between TOPEX/Poseidon observation and thermal contribution at a 1–2 cm level. The potential causes for these differences are discussed, including water mass exchanges among the atmosphere, land, and oceans, and error sources in the steric result and geophysical corrections applied in TOPEX/Poseidon data. Received: 25 September 1998 / Accepted: 13 July 1999     

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     

Evaluation of pre-CHAMP gravity models `GRIM5-S1' and `GRIM5-C1' with satellite crossover altimetry

 The new GFZ/GRGS gravity field models GRIM5-S1 and GRIM5-C1, currently used as initial models for the CHAMP mission, have been compared with other recent models (JGM 3, EGM 96) for radial orbit accuracy (by means of latitude lumped coefficients) in computations on altimetry satellite orbits. The bases for accuracy judgements are multi-year averages of crossover sea height differences from Geosat and ERS 1/2 missions. This radially sensitive data is fully independent of the data used to develop these gravity models. There is good agreement between the observed differences in all of the world's oceans and projections of the same errors from the scaled covariance matrix of their harmonic geopotential coefficients. It was found that the tentative scale factor of five for the formal standard deviations of the harmonic coefficients of the new GRIM fields is justified, i.e. the accuracy estimates, provided together with the GRIM geopotential coefficients, are realistic. Received: 20 February 2001 / Accepted: 24 October 2001     

利用卫星测高数据反演海洋重力异常研究

全面研究了利用卫得测高数据反演海洋重力异常3种主要方法（即Stokes数据解析反解以及逆Vening-Meinesz公式）的技术特点，建立了3种算法的数学模型及其谱计算式，在以1440阶次位模型定义的标准场中完成了3种算法的数值比较和内部检核，通过仿真试验实现了3种算法的可靠性和稳定性检验，最后，本文利用卫得测高实测对南中国海地区的海洋重力异常进行了实际反演，并将反演结果同船测数据进行了比较。     

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     

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.     

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     

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     

The filtering effect of orbit correction on geopotential errors

Summary The geophysical interpretation of satellite tracking residuals generally ignores the filtering effect of initial orbit correction on the true errors of the model. While the filtered information is usually regarded as lost, knowing the spectral characteristics of the filter is a great aid in the detailed interpretation of residuals, especially of global data sets. In this regard, we derive the filter characteristics (admittances) of orbit correction in the presence of geopotential-caused trajectory errors. We then apply the filter to determine the likely power of the lost radial information in crossover differences of sea heights determined from satellite altimetry or in the latitude lumped coefficients derived from them. For example, we find that resonant geopotential information with periods longer than the corrected orbit's arc length is largely lost in residual crossover data. Results are given for GEOSAT, ERS-1 and TOPEX/Poseidon in their Exact Repeat Missions, using calibrated variancecovariance matrices of the harmonic geopotential coefficients of several recent Earth gravity models. To prove that filtering is important, we first employed a simple cut of all perturbing terms with periods longer than the general tracking period (4 days for GEOSAT and ERS-1, and 10 days for TOPEX). But the cut is too crude a method from a theoretical viewpoint, and thus, we developed two new filters. A comparison of their admittances explains the differences (and sometimes anomalous behaviour) between them and the cut. Many numerical examples (single-satellite crossover errors and latitude lumped coefficient errors, as projected from the variance-covariance matrices) are presented.This paper has been presented during the Panel on Satellite Dynamics, at COSPAR 1994, in Hamburg, Germany.     

Accuracy assessment and refinement of the JGM-2 and JGM-3 gravity fields for radial positioning of ERS-1

ERS-1 radial positioning using the JGM-2 and JGM-3 gravity fields is assessed by analysing dual crossovers with TOPEX/Poseidon, neither field containing ERS-1 data. This method allows a more complete recovery of ERS-1 radial orbit error, specifically of the previously unattainable mean geographical error. The global analysis shows that the theoretical error derived from the JGM-2 covariance matrix is realistic and that JGM-3 represents a slight improvement, at least at the inclination of ERS-1. A latitudinal-based study in the southern ocean indicates possible weaknesses in both fields, notably for low and resonant geopotential orders m. A refinement of JGM-2, RGM-2, is undertaken through inclusion of ERS-1 and STELLA laser tracking and ERS-1 altimetry, reducing several of its deficiencies. Received: 14 May 1996 / Accepted: 17 February 1997     

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.     

卫星测高在陆地湖泊水位变化监测中的应用

本文以呼伦湖为例,利用TOPEX/Poseidon测高卫星(1993～2002年)的GDR和WDR数据以及Jason-1(2002～2005年)的GDR数据对湖泊水位变化进行研究,得到湖泊水位异常的变化趋势。通过利用小波分析的方法对卫星测高获得的湖水位异常时间序列分析,得出了呼伦湖的水位年际变化周期,证明了利用卫星测高进行卫星测高监测湖泊水位变化的可行性。     

Quasi-stationary sea surface topography estimation by the multiple input/output method

 Multiple input/multiple output system theory (MIMOST) is briefly presented, and the application of the method to the quasi-stationary sea surface topography (QSST) estimation and the filtering of the input observations are discussed. The repeat character of satellite altimetry missions provides more than one sample of the measured sea surface height (SSH) field, and an approximation of the input signal and error power spectral densities can be determined using this successive information. A case study in the Labrador Sea is considered using SSHs from ERS1 phases C and G, ERS1-GM, ERS2 phase A and TOPEX/POSEIDON altimetric missions in combination with shipborne gravity anomalies. The time period of the observations in this study is from 1993 to 1998. Some comparisons between the techniques used for the power spectral density approximation are carried out and some remarks on the properties of the estimated QSST are presented. Received: 19 October 1999 / Accepted: 23 October 2000     

Analysis of some systematic errors affecting altimeter-derived sea surface gradient with application to geoid determination over Taiwan

This paper analyzes several systematic errors affecting sea surface gradients derived from Seasat, Geosat/ERM, Geosat/GM, ERS-1/35d, ERS-1/GM and TOPEX/POSEIDON altimetry. Considering the data noises, the conclusion is: (1) only Seasat needs to correct for the non-geocentricity induced error, (2) only Seasat and Geosat/GM need to correct for the one cycle per revolution error, (3) only Seasat, ERS-1/GM and Geosat/GM need to correct for the tide model error; over shallow waters it is suggested to use a local tide model not solely from altimetry. The effects of the sea surface topography on gravity and geoid computations from altimetry are significant over areas with major oceanographic phenomena. In conclusion, sea surface gradient is a better data type than sea surface height. Sea surface gradients from altimetry, land gravity anomalies, ship gravity anomalies and elevation data were then used to calculate the geoid over Taiwan by least-squares collocation. The inclusion of sea surface gradients improves the geoid prediction by 27% when comparing the GPS-derived and the predicted geoidal heights, and by 30% when comparing the observed and the geoid-derived deflections of the vertical. The predicted geoid along coastal areas is accurate to 2 cm and can help GPS to do the third-order leveling. Received 22 January 1996; Accepted 13 September 1996     

利用TOPEX/Poseidon卫星测高资料监测全球海平面变化

本文利用１９９３年１月至１９９９年５月的ＴＯＰＥＸ／Ｐｏｓｅｉｄｏｎ卫星测高数据计算了全球海平变化。海平面模型误差和不恰当的加权方法都会影响全球海平面序列的周年变化,潮汐模型误差对于季节性变化有明显的影响。利用两个高度计算的海平面变化存在周年变化存在明显的差异,利用ＴＯＰＥＸ高度计和Ｐｏｓｅｉｄｏｎ高度计得到的该时段全球平均海平面的变化率分别为２．０ｍｍ／ａ和－０．５ｍｍ／ａ。与南方涛动指数的比较     

基于Laplace方程的垂线偏差法反演全球海域重力异常

海洋重力场模型反演的质量主要依赖于采用测高数据的精度、空间分辨率和数据分布密集程度。本文联合Geosat GM/ERM、ERS-1 GM/ERM、TOPEX/Poseidon、Envisat、Cryosat-2、Jason-1 ERM/GM和SARAL/AltiKa等多种测高观测数据集,深入比较了多种波形重跟踪算法的效果,回波数据重跟踪处理后的沿轨海面高标准差。统计表明,Sandwell算法优于MLE-4算法、Davis阈值法、改进阈值法和β参数拟合法;基于不同测高数据波形重采样的结果给出了沿轨海面梯度计算中低通滤波的参数选择方法,并采用Sandwell提出的垂线偏差法,反演了全球海域1′×1′的重力场模型。检核表明,反演结果与DTU13和SIO V23.1模型检核的差值均方根分别为3.4、1.8 mGal,与NGDC船测数据的检核精度为4~8 mGal,且本文模型在部分典型海区内精度更优。     

Gravity recovery using COSMIC GPS data: application of orbital perturbation theory

 COSMIC is a joint Taiwan–US mission to study the atmosphere using the Global Positioning System (GPS) occultation technique. Improved formulas are developed for the radial, along-track, and cross-track perturbations, which are more accurate than the commonly used order-zero formulas. The formulas are used to simulate gravity recovery using the geodetic GPS data of COSMIC in the operational phase. Results show that the EGM96 model can be improved up to degree 26 using 1 year of COSMIC data. TOPEX/POSEIDON altimeter data are used to derive a temporal gravity variation. COSMIC cannot reproduce this gravity variation perfectly because of data noise and orbital configuration, but the recovered field clearly shows the gravity signature due to mass movement in an El Ni?o. Received: 3 March 2000 / Accepted: 10 November 2000     

Prediction of surface horizontal displacements, and gravity and tilt changes caused by filling the Three Gorges Reservoir

 Horizontal displacements, and gravity and tilt changes induced by filling the Three Gorges Reservoir are modeled using elastic loading Green functions. When the water surface reaches its highest level, the effects become maximum on the reservoir banks. The longitudinal and latitudinal components of the horizontal displacements reach −8.2 and 7.7 mm respectively, gravity is increased by up to 3.4 mGal, and the prime vertical and meridian components of the tilt changes are −7.8 and −17.5 arcseconds respectively. Accordingly, the filling of the reservoir will influence values observed from global positioning system (GPS), gravimetry and tilt measurements in the area. The results given can be used to provide important corrections for extracting earthquake-related signals from observed data. Received: 19 January 2001 / Accepted: 3 September 2001