Modelling the Earth's Geomagnetic Field to High Degree and Order

作者：Dave Schmitz, J. Meyer, Joseph C. Cain · 发表于：Geophysical Journal International · 年份：1989 · DOI：10.1111/j.1365-246x.1989.tb00512.x · 被引用次数：22 · 研究领域：Geomagnetism and Paleomagnetism Studies、Geophysics and Gravity Measurements、Solar and Space Plasma Dynamics

SUMMARY We present a method for modelling the Earth’s magnetic field to very high degree and order in terms of spherical harmonics. The method exploits the orthogonality of the spherical functions, using, in part, the method of Gauss-Legendre quadrature. This method is compared to a simpler quadrature method (Newton-Cotes). We show that the GaussLegendre technique is more accurate in most cases than Newton-Cotes quadrature, and in all cases, even where the two give about the same results, that the Gauss-Legendre method is more efficient in that it requires less data and hence less computation. The two quadrature methods are applied to sets of radial field data computed from an n = 29 model which simulate Mugsut observations. The results are that direct integration of a complete global coverage of observations using Newton-Cotes produces errors in the spatial spectrum comparable to that of the geomagnetic field at n = 9, whereas Gauss-Legendre gives exact recovery. When 7” of polar region data are removed to simulate the Mugsat orbit, both methods fail, although Gauss-Legendre gives somewhat less noisy results. However, when the analysis is performed on residuals to a field truncated at n = 15, both methods give comparable levels of noise. Simple interpolation of data over the pole is seen to reduce the errors significantly beyond n = 50. Addition of synthetic noise is seen to provide a means of evaluating the accuracy of coefficients derived from actual data. A theoretical re...