Earth Gravitational Model

Summary

The Earth Gravitational Models (EGM) are a series of geopotential models of the Earth published by the National Geospatial-Intelligence Agency (NGA). They are used as the geoid reference in the World Geodetic System.

A view of the Earth's geoid, as provided by EGM96

The NGA provides the models in two formats: as the series of numerical coefficients to the spherical harmonics which define the model, or a dataset giving the geoid height at each coordinate at a given resolution.[1]

Three model versions have been published: EGM84 with n=m=180, EGM96 with n=m=360, and EGM2008 with n=m=2160. n and m are the degree and orders of harmonic coefficients; the higher they are, the more parameters the models have, and the more precise they are. EGM2008 also contains expansions to n=2190.[1] Developmental versions of the EGM are referred to as Preliminary Gravitational Models (PGMs).[2]

Each version of EGM has its own EPSG code as a vertical datum.

History edit

EGM84 edit

The first EGM, EGM84, was defined as a part of WGS84 along with its reference ellipsoid. WGS84 combines the old GRS 80 with the then-latest data, namely available Doppler, satellite laser ranging, and Very Long Baseline Interferometry (VLBI) observations, and a new least squares method called collocation.[3] It allowed for a model with n=m=180 to be defined, providing a raster for every half degree (30', 30 minute) of latitude and longitude of the world.[4] NIMA also computed and made available 30′×30′ mean altimeter derived gravity anomalies from the GEOSAT Geodetic Mission. 15′×15′ is also available.[5]

EGM96 edit

EGM96 from 1996 is the result of a collaboration between the National Imagery and Mapping Agency (NIMA), the NASA Goddard Space Flight Center (GSFC), and the Ohio State University. It took advantage of new surface gravity data from many different regions of the globe, including data newly released from the NIMA archives. Major terrestrial gravity acquisitions by NIMA since 1990 include airborne gravity surveys over Greenland and parts of the Arctic and the Antarctic, surveyed by the Naval Research Lab (NRL) and cooperative gravity collection projects, several of which were undertaken with the University of Leeds. These collection efforts have improved the data holdings over many of the world's land areas, including Africa, Canada, parts of South America and Africa, Southeast Asia, Eastern Europe, and the former Soviet Union. In addition, there have been major efforts to improve NIMA's existing 30' mean anomaly database through contributions over various countries in Asia. EGM96 also included altimeter derived anomalies derived from ERS-1 by Kort & Matrikelstyrelsen (KMS), (National Survey and Cadastre, Denmark) over portions of the Arctic, and the Antarctic, as well as the altimeter derived anomalies of Schoene [1996] over the Weddell Sea. The raster from EGM96 is provided at 15'x15' resolution.[1]

EGM96 is a composite solution, consisting of:[6]

  1. a combination solution to degree and order 70,
  2. a block diagonal solution from degree 71 to 359,
  3. and the quadrature solution at degree 360.

PGM2000A is an EGM96 derivative model that incorporates normal equations for the dynamic ocean topography implied by the POCM4B ocean general circulation model.

EGM2008 edit

The official Earth Gravitational Model EGM2008 has been publicly released by the National Geospatial-Intelligence Agency (NGA) EGM Development Team. Among other new data sources, the GRACE satellite mission provided a very high resolution model of the global gravity. This gravitational model is complete to spherical harmonic degree and order 2159 (block diagonal), and contains additional coefficients extending to degree 2190 and order 2159. It provides a raster of 2.5′×2.5′ and an accuracy approaching 10 cm. 1'×1' is also available[7] in non-float but lossless PGM,[5][8] but original .gsb files are better.[9] Indeed, some libraries like GeographicLib use uncompressed PGM, but it is not original float data as was present in .gsb format. That introduces an error of up to 0.3 mm because of 16 bit quantisation, using lossless float GeoTIFF or original .gsb files is a good idea.[5][10] The two grids can be recreated by using program in Fortran and source data from NGA.[11] "Test versions" of EGM2008 includes PGM2004, 2006, and 2007.[2]

As with all spherical harmonic models, EGM2008 can be truncated to have fewer coefficients with lower resolution.

EGM2020 edit

EGM2020 is to be a new release (still not released as of January 2024) with the same structure as EGM2008, but with improved accuracy by incorporating newer data.[12] It was originally planned to be released in April 2020.[13] The precursor version XGM2016 (X stands for experimental) was released in 2016 up to degree and order (d/o) 719.[14] XGM2019e was released in 2020 up to spheroidal d/o 5399 (that corresponds to a spatial resolution of 2′ which is ~4 km) and spherical d/o 5540 with a different spheroidal harmonic construction followed by conversion back into spherical harmonics.[15][16] XGM2020 was also released recently.[17]

See also edit

References edit

  1. ^ a b c "WGS 84 Earth Gravitational Model". earth-info.nga.mil. Archived from the original on 27 March 2013. Retrieved 30 July 2019.
  2. ^ a b Pavlis, Nikolaos K.; Holmes, Simon A.; Kenyon, Steve C.; Factor, John K. (April 2012). "The development and evaluation of the Earth Gravitational Model 2008 (EGM2008)". Journal of Geophysical Research: Solid Earth. 117 (B4). Bibcode:2012JGRB..117.4406P. doi:10.1029/2011JB008916.
  3. ^ Ruffhead, A. (April 1987). "An introduction to least-squares collocation". Survey Review. 29 (224): 85–94. doi:10.1179/003962687791512662.
  4. ^ "WGS 84, N=M=180 Earth Gravitational Model". earth-info.nga.mil.
  5. ^ a b c "GeographicLib: Geoid height". geographiclib.sourceforge.io. Retrieved 2022-01-31.
  6. ^ Lemoine, F. G., S. C. Kenyon, J. K. Factor, R.G. Trimmer, N. K. Pavlis, D. S. Chinn, C. M. Cox, S. M. Klosko, S. B. Luthcke, M. H. Torrence, Y. M. Wang, R. G. Williamson, E. C. Pavlis, R. H. Rapp and T. R. Olson (1998). The Development of the Joint NASA GSFC and the National Imagery and Mapping Agency (NIMA) Geopotential Model EGM96. NASA/TP-1998-206861, July 1998. Partly available online.
  7. ^ "EPSG:3859". Retrieved 2022-02-04.
  8. ^ "EGM2008 height". Retrieved 2022-02-04.
  9. ^ "EGM96 and EGM2008 Geoids". www.usna.edu. Retrieved 2022-01-31.
  10. ^ Jochem, Lesparre (8 December 2019). "[PROJ] Is 16-bit quantization of values / (sub-)millimetric error in grids (sometimes) acceptable ?". Retrieved 2022-01-31.
  11. ^ "EGM2008 - WGS 84 Version". Archived from the original on 2021-02-18. Retrieved 2022-02-04.
  12. ^ Barnes, D.; Factor, J. K.; Holmes, S. A.; Ingalls, S.; Presicci, M. R.; Beale, J.; Fecher, T. (1 December 2015). Earth Gravitational Model 2020. AGU Fall Meeting. pp. G34A–03. Bibcode:2015AGUFM.G34A..03B.
  13. ^ Daniel Barnes; Jim Beale; Sarah Ingalls; Howard Small; Rose Ganley; Cliff Minter; Manny Presicci (September 18, 2019). "EGM2020: Updates" (PDF).
  14. ^ Pail, R.; Fecher, T.; Barnes, D.; Factor, J. F.; Holmes, S. A.; Gruber, T.; Zingerle, P. (April 2018). "Short note: the experimental geopotential model XGM2016". Journal of Geodesy. 92 (4): 443–451. Bibcode:2018JGeod..92..443P. doi:10.1007/s00190-017-1070-6. S2CID 126360228.
  15. ^ Zingerle, P.; Pail, R.; Gruber, T.; Oikonomidou, X. (July 2020). "The combined global gravity field model XGM2019e" (PDF). Journal of Geodesy. 94 (7): 66. Bibcode:2020JGeod..94...66Z. doi:10.1007/s00190-020-01398-0.
  16. ^ "The experimental gravity field model XGM2019e". dataservices.gfz-potsdam.de. Retrieved 2022-01-30.
  17. ^ Zingerle, Philipp (2020-05-02). "High-resolution combined global gravity field modelling: the d/o 5,400 XGM2020 mode" (PDF).

External links edit

  • EGM96: The NASA GSFC and NIMA Joint Geopotential Model
  • Earth Gravitational Model 2008 (EGM2008)
  • GeographicLib provides a utility GeoidEval (with source code) to evaluate the geoid height for the EGM84, EGM96, and EGM2008 Earth gravity models. Here is an online version of GeoidEval.
  • The Tracker Component Library from the United States Naval Research Laboratory is a free Matlab library with a number of gravitational synthesis routines. The function getEGMGeoidHeight can be used to evaluate the geoid height under the EGM96 and EGM2008 models. Additionally, the gravitational potential, acceleration, and gravity gradient (second spatial derivatives of the potential) can be evaluated using the spherHarmonicEval function, as demonstrated in DemoGravCode.