Ultra-high-Degree Surface Spherical Harmonic Analysis Using the Gauss–Legendre and the Driscoll/Healy Quadrature Theorem and Application to Planetary Topography Models of Earth, Mars and Moon

详细信息    查看全文

• 作者：Moritz Rexer ; Christian Hirt
• 关键词：Spherical harmonic analysis ; Quadrature ; Gauss–Legendre ; Driscoll/Healy ; Topography ; Digital elevation model ; Earth ; Mars ; Moon
• 刊名：Surveys in Geophysics
• 出版年：2015
• 出版时间：November 2015
• 年：2015
• 卷：36
• 期：6
• 页码：803-830
• 全文大小：10,169 KB
• 参考文献：Abrykosov O, F?rste C, Gruber C, Shako R, Barthelmes F (2012) Harmonic analysis of the DTU10 global gravity anomalies. In: Abbasi A, Giesen N (eds) EGU General Assembly conference abstracts, vol 14, p 4945
  Andersen O, Knudsen P, Kenyon S, Factor J, Holmes S (2013) The dtu13 global marine gravity field—first evaluation. Technical report, DTU Space - National Space Institute
  Arabelos D, Tscherning C (1998) The use of least squares collocation method in global gravity field modeling. Phys. Chem. Earth 23(1):1-2. doi:10.-016/?S0079-1946(97)00234-6 CrossRef
  Balmino G, Vales N, Bonvalot S, Briais A (2012) Spherical harmonic modelling to ultra-high degree of Bouguer and isostatic anomalies. J Geodesy 86(7):499-20. doi:10.-007/?s00190-011-0533-4 CrossRef
  Bartusch M, Berg H, Siebertz O (2008) The TanDEM-X Mission. In: 7th European conference on synthetic aperture radar (EUSAR), pp 1-
  Becker J, Sandwell D, Smith W, Braud J, Binder B, Depner J, Fabre D, Factor J, Ingalls S, Kim SH, Ladner R, Marks K, Nelson S, Pharaoh A, Trimmer R, Von Rosenberg J, Wallace G, Weatherall P (2009) Global bathymetry and elevation data at 30 arc seconds resolution: Srtm30\_plus. Mar Geodesy 32(4):355-71. doi:10.-080/-149041090329776- CrossRef
  Bucha B, Janák J (2013) A MATLAB-based graphical user interface program for computing functionals of the geopotential up to ultra-high degrees and orders. Comput Geosci 56:186-96CrossRef
  Claessens S (2006) Solutions to ellipsoidal boundary value problems for gravity field modelling. PhD thesis, Curtin University of Technology
  Claessens S, Hirt C (2013) Ellipsoidal topographic potential—new solutions for spectral forward gravity modelling of topography with respect to a reference ellipsoid. J Geophys Res 118(11):5991-002. doi:10.-002/-013JB010457 CrossRef
  Colombo O (1981) Numerical methods for harmonic analysis on the sphere. Technical report. Report no. 310, The Ohio State University
  Dassios G (2012) Ellipsoidal harmonics. Cambridge University Press, Cambridge. doi:10.-017/?CBO9781139017749-/span> CrossRef
  Driscoll J, Healy D (1994) Computing Fourier transforms and convolutions on the 2-sphere. Adv. Appl. Math. 15(2):202-50. doi:10.-006/?aama.-994.-008 CrossRef
  ESA (1999) Gravity field and steady-state ocean circulation mission. Report for the mission selection of the four candidate earth explorer missions (ESA SP-1233(1)), European Space Agency
  Farr T, Rosen P, Caro E, Crippen R, Duren R, Hensley S, Kobrick M, Paller M, Rodriguez E, Roth L, Seal D, Shaffer S, Shimada K, Umland J, Werner M, Oskin M, Burbank D, Alsdorf D (2007) The shuttle radar topography mission. Rev Geophys 45:RG2004. doi:10.-029/-005RG000183
  Fecher T, Pail R, Gruber T (2013) Global gravity field modeling based on GOCE and complementary gravity data. Int J Appl Earth Obs Geoinf 35:120-27. doi:10.-016/?j.?jag.-013.-0.-05 CrossRef
  Fukushima T (2012) Numerical computation of spherical harmonics of arbitrary degree and order by extending exponent of floating point numbers. J Geodesy 86(4):271-85. doi:10.-007/?s00190-011-0519-2 CrossRef
  Fukushima T (2015) Personal communication, IUGG Prague
  Gruber C (2011) A study on the Fourier composition of the associated Legendre functions; suitable for applications in ultra-high resolution. Scientific technical report 11/04, German Research Centre for Geosciences (GFZ), Potsdam. doi:10.-312/?GFZ.?b103-11041
  Gruber C, Novak P, Sebera J (2011) FFT-based high-performance spherical harmonic transformation. Stud Geophys Geod 55:489-00CrossRef
  Gruber C, Barthelmes F, Flechtner F, Novak P (2014) Derivation of topographic potential from global DEM models. In: Rizos C, Willis P (eds) Earth on the edge: science for a sustainable planet: proceedings of the IAG General Assembly, Melbourne, Australia, 28 June- July 2011, vol 139. Springer, Berlin, pp 535-42
  Hirt C, Kuhn M (2014) A band-limited topographic mass distribution generates a full-spectrum gravity field—gravity forward modelling in the spectral and spatial domain revisited. J Geophys Res Solid Earth 119(4):3646-661. doi:10.-002/-013JB010900 CrossRef
  Hirt C, Rexer M (2015) Earth 2014: 1-shape, topography, bedrock and ice-sheet models—available as gridded data and degree 10,800 spherical harmonics. Int J Appl Earth Obs Geoinf 39:103-12. doi:10.-016/?j.?jag.-015.-3.-01 CrossRef
  Hofsommer D (1957) On the expansion of a function in a series of spherical harmonics. Technical report. Report no. R344A, Computation Department of the Mathematical Centre, Amsterdam
  Hofsommer D, Potters M (1960) Table of Fourier coefficients of associated Legendre functions. Report r 478. knaw, Computational Department of the Mathematical Centre, Amsterdam
  Holmes S, Featherstone W (2002) A unified approach to the Clenshaw summation and the recursive computation of very high degree and order normalised associated Legendre functions. J Geodesy 76:279-99. do
• 作者单位：Moritz Rexer (1) (2)
  Christian Hirt (1) (2) (3)
  
  1. Institute for Astronomical and Physical Geodesy, Technische Universit?t München, Arcisstrasse 21, 80333, Munich, Germany
  2. Institute for Advanced Study, Technische Universit?t München, Lichtenbergstr. 2 a, 85748, Garching, Germany
  3. Western Australian Geodesy Group, Department of Spatial Sciences, The Institute for Geophysical Research, Curtin University of Technology, GPO Box U1987, Perth, WA, 6845, Australia
  
• 刊物类别：Earth and Environmental Science
• 刊物主题：Earth sciences
  Geophysics and Geodesy
  Geosciences
  Astronomy
  
• 出版者：Springer Netherlands
• ISSN：1573-0956

文摘

In geodesy and geophysics, spherical harmonic techniques are popular for modelling topography and potential fields with ever-increasing spatial resolution. For ultra-high-degree spherical harmonic modelling, i.e. degree 10,000 or more, classical algorithms need to be extended to avoid under- or overflow problems associated with the computation of associated Legendre functions (ALFs). In this work, two quadrature algorithms—the Gauss–Legendre (GL) quadrature and the quadrature following Driscoll/Healy (DH)—and their implementation for the purpose of ultra-high (surface) spherical harmonic analysis of spheroid functions are reviewed and modified for application to ultra-high degree. We extend the implementation of the algorithms in the SHTOOLS software package (v2.8) by (1) the X-number (or Extended Range Arithmetic) method for accurate computation of ALFs and (2) OpenMP directives enabling parallel processing within the analysis. Our modifications are shown to achieve feasible computation times and a very high precision: a degree-21,600 band-limited (=frequency limited) spheroid topographic function may be harmonically analysed with a maximum space-domain error of \(3 \times 10^{-5}\) and \(5 \times 10^{-5}\) m in 6 and 17 h using 14 CPUs for the GL and for the DH quadrature, respectively. While not being inferior in terms of precision, the GL quadrature outperforms the DH algorithm in terms of computation time. In the second part of the paper, we apply the modified quadrature algorithm to represent for—the first time—gridded topography models for Earth, Moon and Mars as ultra-high-degree series expansions comprising more than 2 billion coefficients. For the Earth’s topography, we achieve a resolution of harmonic degree 43,200 (equivalent to ~500 m in the space domain), for the Moon of degree 46,080 (equivalent to ~120 m) and Mars to degree 23,040 (equivalent to ~460 m). For the quality of the representation of the topographic functions in spherical harmonics, we use the residual space-domain error as an indicator, reaching a standard deviation of 3.1 m for Earth, 1.9 m for Mars and 0.9 m for Moon. Analysing the precision of the quadrature for the chosen expansion degrees, we demonstrate limitations in the implementation of the algorithms related to the determination of the zonal coefficients, which, however, do not exceed 3, 0.03 and 1 mm in case of Earth, Mars and Moon, respectively. We investigate and interpret the planetary topography spectra in a comparative manner. Our analysis reveals a disparity between the topographic power of Earth’s bathymetry and continental topography, shows the limited resolution of altimetry-derived depth (Earth) and topography (Moon, Mars) data and detects artefacts in the SRTM15 PLUS data set. As such, ultra-high-degree spherical harmonic modelling is directly beneficial for global inspection of topography and other functions given on a sphere. As a general conclusion, our study shows that ultra-high-degree spherical harmonic modelling to degree ~46,000 has become possible with adequate accuracy and acceptable computation time. Our software modifications will be freely distributed to fill a current availability gap in ultra-high-degree analysis software. Keywords Spherical harmonic analysis Quadrature Gauss–Legendre Driscoll/Healy Topography Digital elevation model Earth Mars Moon