CHAPTER 5 GRAVITY ANOMALY OVER EAST CHINA SEA AND TAIWAN
5.5 Case Studies
5.5.2 The Taiwan Strait…
Next we carried out experiments over TS using the same four cases as in the ECS. We used shipborne gravity data along Tracks 1 and 2 (Figure 5.1) to evaluate the altimeter-derived gravity anomalies. These shipborne gravity data were compiled
by Hsu et al. (1998), who has crossover adjusted the shipborne gravity data and removed bad values. Table 5.5 shows the results of the comparisons between altimeter-derived and shipborne gravity anomalies in the four cases. The conclusion from Table 5.5 is similar to what has been drawn from Table 5.4, that is, the least square collocation method with DOV produces the worst result. Again, use of differenced heights produces the best results provided that the same altimeter-gravity conversion method is used.
Since land gravity data are available along the coast of TS, we also assessed the impact of land gravity data on the accuracy of altimeter-derived gravity anomaly.
Figure 5.11 shows the distribution of land gravity and altimeter data around Taiwan.
Note that there is no altimeter at the immediate coastal waters off the west coast of Taiwan. We experimented with the method of LSC using differenced heights, and with and without land gravity data (two cases). Table 5.6 shows the statistics of the differences between shipborne and altimeter-derived gravity anomalies. Figure 5.12 shows the shipborne and altimeter-derived gravity anomalies along Tracks 1 and 2.
The patterns of the difference along Tracks 1 and 2 are similar to those for Tracks c1217 and dmm07 (Figure 5.8). That is, the agreement between shipborne and altimeter-derived gravity anomalies at large spatial scales is better than that at small spatial scales. Furthermore, Tracks 1 and 2 are only tens of km off the west coast of Taiwan and the density of altimeter data is low along theses two tracks (Figure 5.11).
In general, difference of gravity anomaly increases with decreasing density of altimeter data. The density of land gravity data is relatively high over a zone from
5° .
22 N to 24.5°N, and here the agreement between altimeter-derived (with land gravity) and shipborne gravity anomalies is better than that in other parts of the ship tracks. In the areas south of 22.5°N and north of 24.5°N, only few land gravity data
are used, so there is virtually no difference between the altimeter-only gravity anomalies and the altimeter-land gravity anomalies.
Figure 5.11: Distribution of land gravity data and altimeter data around Taiwan.
Figure 5.12: Gravity anomalies along Track 1 and Track 2 in the Taiwan Strait.
Again, we investigated the possible causes of the difference between altimeter-derived and shipborne gravity anomalies in the TS. Figure 5.13 and 5.14 shows the same quantities as in Figure 5.9 and 5.10, but in the TS. Again, along the two ship tracks, the patterns of gravity anomaly difference are quite random and are
not correlated with depth and tidal height difference. However, there is a strong correlation (about 0.9) between standard error of ERS-1 SSH and tidal height difference along both Tracks 1 and 2. Again, tidal height difference increases as decreasing depth. The conclusions from the analyses associated with Figures 5.9 &
5.10 and Figures 5.13 & 5.14 are in agreement with the conclusion drawn in Section 5.2: the major source of large standard error of SSH and large error of altimeter-derived gravity anomaly is tide model error. Another cause of degraded altimeter-derived gravity anomaly near coasts, which is not investigated in this paper, is low altimeter data density due to data editing. Even the data editing criterion near coasts is relaxed to increase data density, the additional altimeter data may not be of good quality for gravity derivation. One way to improve altimeter data quality near coasts is to retrack waveforms of altimeter ranging. For example, Deng et al. (2003) have obtained improved T/P SSHs by waveform retracking over the Australian coasts.
Currently, globally retracked ERS-1 and Geosat waveforms are available (Lillibridge et al., 2004; Smith et al., 2004), and have been shown to produce improved marine gravity fields.
Figure 5.13: Time series of normalized difference of gravity anomaly, standard error of ERS-1 SSH, tide model difference and depth, along Track 1
Figure 5.14: Time series of normalized difference of gravity anomaly, standard error of ERS-1 SSH, tide model difference and depth, along Track 2.
CHAPTER 6
CONCLUSIONS AND RECOMMENDATIONS
6.1 Conclusions
In this thesis, Chapter 2 first introduces 4 methods for gravity anomaly derivations from altimetry, which include a completely different kind of case, LSC with differenced height. Chapter 2 and 3 describe the work of an improved computation of global MSSH and gravity anomaly grids using improved multi-satellite altimeter data sets and a new procedure. Various tests have been made in order to find an optimal set of parameters for the computation. The procedure is based on the deflection-geoid and inverse Vening Meinesz formulas, as well as the remove-restore concept. The crossover adjustment of SSHs is not used because of the use of DOV as the altimeter data type in the computation. Using DOV is particularly advantageous over areas with sparse altimeter data and in the case of long wavelength error contained in altimeter data. The comparisons of the NCTU01 MSSH with the T/P and the ERS-1 MSSH result in overall RMS differences of 5.0 and 3.1 cm in SSH, respectively, and 7.1 and 3.2 µrad in SSH gradient, respectively.
The agreements between the predicted and shipborne gravity anomalies range from 3.0 to 13.4 mgals, depending on the gravity signatures and the altimeter data noise, the later being affected by instrument noise, sea state and accuracies of geophysical corrections. The NCTU01 MSSH model outperforms the NASA/GSFC model (Wang, 2001) and the CLS model, and the NCTU01 gravity anomaly model has a better accuracy than those of the models of Sandwell and Smith (1997) and Hwang et al.
(1998). The gravity anomalies described in Chapter 3 can be used to derive mean gravity anomalies for modeling high-degree geopotential coefficients, and to study the spectral properties of the Earth’s gravity field.
Chapter 4 compares two methods with three kinds of altimeter data types of gravity anomaly derivation from altimetry data. For the four computations of two methods, the RMS differences between altimetry-derived gravity anomalies and shipborne gravity anomalies are 9.06 (differenced height) and 9.59, 10.26 (height slope) and 9.77, 10.44 and 13.10, 10.73 and 11.86 mgals, over TS and ECS, respectively. Comparing to NCTU01 gravity anomaly model, the RMS differences are further reduced by 1.67 and 2.27 mgal in these two areas. The method of LSC with differenced heights yields promising results and produced the best accuracy over TS and ECS. An iterative method to remove outliers in along-track altimeter data improves altimeter data quality and improves the accuracy of predicted gravity anomaly. Preprocessing of altimeter data is important in obtaining good estimation of gravity. Good estimation of gravity depends on both good altimeter data and good method. Use of along-track DOV with LSC yields a better result over TS, but gives a worse result over ECS compared to the result from the inverse Vening Meinesz formula. However, the inverse Vening-Meinesz formula with 1D FFT is the fastest methods among the three methods.
Chapter 5 studies the sources of the differences among two global gravity anomaly grids over ECS and TS. We conclude that tide model error is the biggest contributor to the differences. It is found that tide model error, standard deviation of SSH and ocean depth are highly correlated. Also, the complicated sea states over these two areas increase the roughness of the sea surface and hence the noise level of altimeter ranging. As case studies, we experimented with the method of LSC using
differenced heights, and with and without land gravity data. Including land gravity data in the vicinity of coasts enhances the accuracy of altimeter-derived gravity anomalies. Future studies of gravity anomalies recovery will focus on shallow waters over the Yellow Sea, South China Sea and Southeast Asia.
6.2 Recommendations
For future work, we recommend a procedure to improve the accuracy of altimeter-derived gravity anomalies: (1) retrack near-shore waveforms of altimeter to produce waveform-corrected SSHs, (2) use the corrected SSHs to improve tide model, (3) use improved tide model to correct for the ocean tide effect in SSH. Finally, the improved tide-corrected and waveform-corrected altimeter data will lead to an improved gravity field over shallow waters. In addition, the altimeter derived gravity could be improved by merging the altimeter data available today with those derived from future and more accurate missions, and by combining with other sources of data, such as airborne observations and shipborne gravity data to do the prediction procedure.
REFERENCE
Andersen, O.B., Knudsen, P., and Tscherning, C.C. (1996) Investigation of methods for global gravity field recovery from the dense ERS-1 geodetic mission altimetry. In: Rapp R, Nerem RS, Cazenave A (eds), Global gravity field and temporal variations. IAG Symposium 116. Springer, Berlin Heidelberg New York, pp. 218-226.
Andersen, O. and Knudsen, P. (1998) Global marine gravity from the ERS-1 and Geosat geodetic mission altimetry. J. Geophys. Res., 103, pp. 8129-8137.
Andersen, O.B. and Kundsen, P. (2000) The Role of satellite altimetry in gravity field modeling in coastal areas, Phys. Chem. Earth (A), 25, 1, pp. 17-24.
Andersen, O.B., Knudsen, P. and Trimmer, R. (2001) The KMS2001 Global Mean Sea Surface and Gravity Field, IAG2001 Scientific Assembly, Budapest, Hungary, September.
Andersen, O.B., Knudsen, P. and Trimmer, R. (2005) Improved high resolution altimetric gravity field mapping (KMS2002 global marine gravity field). Proc., IAG Symposia, 127, IUGG Sapporo, 2003, Springer, Berlin.
Anzenhofer, M. and Shum, C.K. (2001) Coastal Altimetry and Applications. Rep of Dept of Geodetic Science and Surveying, in press, Ohio State University, Columbus.
Archiving, Validation, and Interpretation of Satellite Oceanographic Data (1996) AVISO User Handbook for Merged TOPEX/Poseidon Products, 3rd Ed.
Toulouse, France.
Balmino Moynot, G.B., Sarrailh, M. and Vales, N. (1987) Free air gravity anomalies over the oceans from SEASAT and GEOS-3 altimetry data, EOS Trans. AGU, 68, pp. 17-19.
Basic, T. and Rapp, R.H. (1992) Ocean wide prediction of gravity anomalies and sea surface heights using GEOS-3, SEASAT and GEOSAT altimeter data, and ETOPO5U bathematric data, Report No. 416, Department of Geodetic Science and Surveying, The Ohio State University, Columbus, Ohio.
Bilitza, D. (1997) International reference ionosphere - status 1995/96.Adv. Space Res., 20, (9), pp. 1751-1754.
Cazenave, A. and Royer, J.Y. (2001) Applications to marine geophysics, in: Satellite Altimetry and Earth Sciences: a Handbook of Techniques and Applications, Fu,
L.L., and Cazenave, A. (ed.), Academic Press, San Diego.
Chen, S.A. (2001) Determinations of tide, circulation and eddied over the South China Sea using TOPEX/POSEIDON altimeter data. PhD thesis, National Chiao Tung University, Hsinchu, Taiwan.
De Boor, C. (1978) A practical guide to splines. Springer, Berlin.
Deng, A., Featherstone, W., Hwang, C. and Berry, P. (2002) Waveform retracking of ERS-1, Mar. Geod., in press.
Deng, X., Featherstone, W.E., Hwang, C. and Berry, P.A.M. (2001) Estimation of Contamination of ERS-2 and POSEIDON Satellite Radar Altimetry Close to the Coasts of Australia, paper submitted to Mar. Geod..
Eanes, R. (1999) Improved ocean tide model from satellite altimetry. Fall Meeting 1999, Am. Geophys Un, San Francisco.
Eren, K. (1980) Spectral Analysis of Geos-3 Altimeter Data and Frequency Domain Collocation, OSU, Rep. 297.
Fairhead, J.D. and Green, C.M. (2001) Pushing the resolution boundaries of marine satellite gravity, IAG 2001 Scientific Assembly, Budapest.
Fu, L.-L., and Cazenave, A. (Eds), (2001) Satellite Altimetry and Earth Sciences: A Handbook of Techniques and Applications. Academic Press, San Diego.
Fu, L.-L. (editor) (2003) Wide-Swath Altimetric Measurement of Ocean Surface Topography, JPL Publication 03-002, Jet Propulsion Laboratory, Pasadena, CA, pp. 67. (ftp://ftp-oceans.jpl.nasa.gov/pub/llf/WSOAreportFinal2.pdf)
Fu, L.L. and Rodriguez, E. (2004) High-Resolution Measurement of Ocean Surface Topography by Radar Interferometry for Oceanographic and Geophysical Applications, American Geophysical Union, San Francisco.
Gerald, C.F., and Wheatley, P.O. (1994) Applied Numerical Analysis, 5th Ed, Addison-Wesley, New York.
Gomez, V., Maravall, A. and Pena, D. (1999) Missing observations in ARIMA models:
skipping approach versus additive outlier approach, J. Econometrics, 88, (2), pp.
341-363.
Haxby, W.F., Karner, G.D., La Brecque, J.L., and Weissel, J.K. (1983) Digital images of combined oceanic and continental data sets and their use in tectonic studies,
Haxby, W.F. (1987) Gravity field of the world’s oceans. A portrayal of gridded geophysical data derived from Seasat radar altimeter measurement of the shape of the ocean surface, National Geophysical Data Center, Boulder, CO.
Heiskanen, W.A. and Moritz, H. (1985) Physical Geodesy. Reprint, Institute of Physical Geodesy, Tech. University of Graz, Austria.
Hernandez, F. and Schaeffer, P. (2000) Altimetric Mean Sea Surfaces and Gravity Anomaly maps inter-comparisons AVI-NT-011-5242-CLS, CLS Ramonville St Agne.
Hipkin, R., (2000) Modeling the geoid and sea-surface topography in coastal areas.
Phys. Chem. Earth (A), 25, (1), pp. 9-16.
Hsu, H.Y. (1997) Global Marine Geoid from Multi-Satellite Altimetry and It’s application, MSD thesis of department of civil engineering, National Chiao Tung University, Hsinchu, Taiwan. (in Chinese)
Hsu, S., Liu, C., Shyu, C., Liu, S., Sibue, J., Lallemand, S., Wang, C. and Reed, D.
(1998) New gravity and Magnetic Anomaly Maps in the Taiwan-Luzon region and their preliminary interpretation, TAO, 9, 3, pp. 509-532.
Hwang, C. (1989) High precision gravity anomaly and sea surface height estimation from Geo-3/Seasat altimetry date, Report No. 399, Department of Geodetic Science and Surveying, The Ohio State University, Columbus, Ohio.
Hwang, C. and Parsons, B. (1995) Gravity anomalies derived from Seasat, Geosat, ERS-1 and TOPEX/POSEIDON altimetry and ship gravity : A case study over the Reykjanes Ridge, Geophys. J. Int.,122, pp. 551-568.
Hwang, C. and Parsons, B. (1996) An optimal procedure for deriving marine gravity from multi-satellite altimetry, Geophys. J. Int., 125, pp. 705-719.
Hwang, C. (1997) Analysis of some systematic errors affecting altimeter-derived sea surface gradient with application to geoid determination over Taiwan. J. Geod., 71, pp. 113-130.
Hwang, C. (1998) Inverse Vening Meinesz formula and deflection-geoid formula:
application to prediction of gravity and geoid over the South China Sea. J.
Geod., 72, pp. 304-312.
Hwang, C., Kao, E.C. and Parsons, B. (1998) Global Derivation of marine gravity anomalies from Seasat, Geosat, ERS-1 and TOPEX/POSEIDON altimeter data. Geophys. J. Int., 134, pp. 449-459.
Hwang, C., Hsu, H. and Jang, R. (2002) Global mean sea surface and marine gravity anomaly from multi-satellite altimetry: applications of deflection-geoid and inverse Vening Meinesz formulae, J. Geod., 76, pp. 407-418.
Hwang, C., Hsu, H.Y. and Deng, X. (2004) Marine gravity anomaly from satellite altimetry: a comparison of methods over Shallow Waters, International Association of Geodesy Symposia, 126, pp. 59-66, Springer, Berlin.
Hwang, C., Shum, C.K. and Li, J.C. (Eds.) (2004a) Satellite Altimetry for Geodesy, Geophysics and Oceanography. International Association of Geodesy Symposia, 126, Springer, Berlin.
Jacobs, G.A., Hur, H.B. and Riedlinger, S.K. (2000) Yellow and East China Seas response to winds and currents. J. Geophys. Res., 105, (C9), pp. 21947-21968.
Jan, S., Chern, C.S., Wang, J. and Chao, S.Y. (2004) The anomalous amplification of M-2 tide in the TS. Geophys. Res. Lett., 31, (7), Art. No. L07308.
Kaiser, R. (1999) Detection and estimation of structural changes and outliers in unobserved components, Comp. Stat., 14, (4), pp. 533-558 .
Kalnay, E., Kanamitsu, M., Kistler, R., Collins, W., Deaven, D., Gandin, L., Iredell, M., Saha, S., White, G., Woollen, J., Zhu, Y., Chelliah, M., Ebisuzaki, W., Higgins, W., Janowiak, J., Mo, K.C., Ropelewski, C., Wang, J., Leetmaa, A., Reynolds, R., Jenne, R. and Joseph, D. (1996) The NCEP/NCAR 40-year reanalysis project. Bull. Am. Meteol. Soc., 77, pp.
437-471.
Knudsen, P. (1987) Estimation and modeling of the local empirical covariance function using gravity and satellite altimeter data, Bull. Geodesique., 61, pp. 145-160.
Koch, K.R. (1970) Gravity anomalies for ocean areas from satellite altimetry, Proceedings of the Second Marine Geodesy Symposium, Marine Technology Society, Washington, D.C.
Koch, K.R. (1987) Parameter Estimation and Hypothesis Testing in Linear Models, Springer, Berlin.
Lefevre, F., Le Provost, C. and Lyard, F.H. (2000) How can we improve a global ocean tide model at a regional scale? A test on the Yellow Sea and the East China Sea. J. Geophys. Res., 105, (C4), pp. 8707-8725.
Cox C.M., Klosko, S.M., Luthcke, S.B., Torrence, M.H., Wang, Y.M., Williamson, R.G., Pavlis, E.C., Rapp, R.H. and Olson, T.R. (1998) The Development of Joint NASA GSFC and the National Imagery and Mapping Agency (NIMA) Geopotential Model EGM96. Rep NASA/TP-1998-206861, National Aeronautics and Space Administration, Greenbelt, MD.
Le Traon, P.Y. and Ogor, F. (1998) ERS-1/2 orbit improvement using TOPEX/POSEIDON: The 2 cm challenge. J. Geophys. Res., 103, pp.
8045-8057.
Lemoine, F.G, Kenyon, S.C., Factor, J.K., Trimmer, R.G., Pavlis, N.K., Chinn, D.S., Cox, C.M., Klosko, S.M, Luthcke, S.B., Torrence, M.H., Wang, Y.M., Williamson, R.G., Pavlis, E.C., Rapp, R.H. and Olson, T.R. (1998) The Development of Joint NASA GSFC and the National Imagery and Mapping Agency (NIMA) Geopotential Model EGM96. Rep NASA/TP-1998-206861, National Aeronautics and Space Administration, Greenbelt, MD.
Levitus, S., Monterey, G.I. and Boyer, T. (1997) Seasonal Variability of Dynamic Height and its Fourier Analysis. NOAA NESDIS Atlas 15, U.S. Gov. Printing Office. Wash., D.C.
Li, J. and Sideris, M. (1997) Marine gravity and geoid determination by optimal combination of satellite altimetry and shipborne gravimetry data. J. Geod., 71, pp. 209-216.
Liang, C.K. (1983) The Adjustment and Combination of Geos-3 and Seasat Altimeter Data. Rep 346, Dept of Geodetic Science and Surveying, Ohio State University, Columbus.
Lillibridge, J.L., Smith, W.H.F, Scharroo, R. and Sandwell, D.T. (2004) The Geosat geodetic mission twentieth anniversary edition data product. AGU 2004 Fall meeting, San Francisco.
Matsumoto, K., Takanezawa, T. and Ooe, M. (2000) Ocean Tide Models Developed by Assimilating TOPEX/POSIDON Altimeter Data into Hydrodynamical Model:
A Global Model and a Regional Model around Japan. J. Oceanogr., 56, pp.
567-581.
Matsumoto, K., Sato, T., Takanezawa, T., and Ooe, M. (2001) GOTIC2: A Program for Computation of Oceanic Tidal Loading Effect. J. Geod. Soc., Japan, 47, pp.
243-248.
Moritz, H. (1980) Advanced Physical Geodesy, Herbert Wichmann, Karlsruhe.
Niwa, Y. and Hibiya, T. (2004) Three-dimensional numerical simulation of M-2 internal tides in the East China Sea. J. Geophys. Res., 109, (C4), Art. No.
C04027.
NOAA (1997) The Complete Geosat Altimeter GDR Handbook, NODC Laboratory for Satellite Altimetry. Silver Spring, Maryland.
Olgiati, A., Balmino, G., Sarrailh, M., and Green, C.M. (1995) Gravity anomalies from satellite altimetry: comparison between computation via geoid heights and via deflections of the vertical. Bull. Geod., 69, pp. 252-260.
Pearson, R.K. (2002) Outliers in process modelling and identification, IEE Trans on Control Tech, 10, pp. 55-63.
Pope, A.J. (1976) The statistics of residual and the detection of outliers, Tech. Rep.
TR-NOS-65-NGS-1, National Ocean Survey, Rockville.
Randel, D.L., VonderHaar, T.H., Ringerud, M.A., Stephens, G.L., Greenwald, T.J. and Combs, C.L. (1996) A new global water vapor dataset. Bull. Am. Meteol. Soc., 77, pp. 1233-1246.
Rapp, R.H. (1983) The determination of geoid undulations and gravity anomalies from Seasat altimeter data. J. Geophys. Res., 88, pp. 1552-1562.
Rapp, R.H. (1985) Detailed gravity anomalies and sea surface height derived from Geo-3/Seasat altimeter data, Report No. 365, Department of Geodetic Science and Surverying, The Ohio State University, Columbus, Ohio.
Rapp, R. and Balasubramania, N. (1992) A Conceptual Formulation of a World Height System. Rep. No. 421, Dept. of Geod. Sci. and Surveying, Ohio State University, Columbus, Ohio.
Rapp, R.H. and Basic, T. (1992) Oceanwide gravity anomalies from Geos-3, Seasat and Geosat altimeter data, Geophys. Res. Lett., 19, pp. 1979-1982.
Rapp, R.H. and Yi, Y. (1997) Role of ocean variability and dynamic ocean topography in the recovery of the mean sea surface and gravity anomalies from satellite altimeter data. J. Geod., 71, pp. 617-629.
Sandwell, D.T. (1992) Antarctic marine gravity field from high-density satellite altimetry, Geophys. J. Int., 109, pp. 437-448.
Sandwell, D.T. and Smith, W.H.F. (1997) Marine gravity anomaly from Geosat and ERS 1 satellite altimetry. J. Geophys. Res., 102, pp. 10039-10054.
Sandwell, D.T. and Smith, W.H.F. (2001) Bathymetric estimation. In: Fu. L,-L. and Canezave, A. (eds), Satellite Altimetry and Earth Sciences, pp 442-457.
Sandwell, D.T. and Smith, W.H.F. (2001) Bathymetric estimation. In Fu, L.-L. and Cazenave, A. (Eds) Satellite Altimetry and Earth Sciences: A Handbook of Techniques and Applications, Academic Press, New York.
Schwarz, K.P., Sideris, M.G. and Forsberg, R. (1990) The use of FFT techniques in physical geodesy. Geophys. J. Int., 100, pp. 485-514.
Seeber, G. (1993) Satellite Geodesy. de Gruyter, Berlin.
Shum, C.K., Woodworth, P.L., Andersen, O.B., Egbert, G.D., Francies, O., King, C., Klosko, S.M., Le Provost, C., Li, X., Molines, J.M., Parke, M.E., Ray, R.D., Schlax, M.G., Stammer, D., Tierney, C.C., Vincent, P. and Wunsch, C.
(1997) Accuracy assessment of recent ocean tide models. J. Geophys. Res., 102, pp. 25173-25194.
Sideris, M.G. (1996) On the use of heterogeneous noisy data in spectral gravity field modelling methods, J. Geod., 70, pp. 470-479.
Smith, W.H.F. and Sandwell, D.T. (2004) Improved global marine gravity field from reprocessing of Geosat and ERS-1 radar altimeter waveforms. AGU 2004 Fall meeting, San Francisco.
Smith, W.H.F. and Wessel, P. (1990) Gridding with continuous curvature splines in tension, Geophys., 55, pp. 293-305.
Torge, W. (1989) Gravimetry, de Gruyter, Berlin.
Trimmer, R.G., Andersen, O. and Driscoll, M.L. (2001) Estimating altimetry 5´×5´
mean gravity anomaly regional accuracies. IAG 2001 Scientific Assembly, Budapest.
Tscherning, C.C. and Rapp, R.H. (1974) Closed covariance expressions for gravity anomalies geoid undulations and deflections of the vertical implied by anomaly degree variance models, Rep 208, Dept of Geod Sci, Ohio State University, Columbus.
Tscherning, C. C. (1997) Geoid determination by Least Squares Collocation, International school for the determination and use of the Geoid, International
Tscherning, C. C. (1997) Geoid determination by Least Squares Collocation, International school for the determination and use of the Geoid, International