Premium
Efficient modeling of the Earth's magnetic field with harmonic splines
Author(s) -
Parker Robert L.,
Shure Loren
Publication year - 1982
Publication title -
geophysical research letters
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 2.007
H-Index - 273
eISSN - 1944-8007
pISSN - 0094-8276
DOI - 10.1029/gl009i008p00812
Subject(s) - spline (mechanical) , mathematics , basis function , harmonic , spherical harmonics , series (stratigraphy) , mathematical analysis , algorithm , computer science , physics , geology , quantum mechanics , paleontology , thermodynamics
The construction of smooth potential field models has many geophysical applications. The recently‐developed method of harmonic splines produces magnetic field models at the core surface which are maximally smooth in the sense of minimization of certain special norms. They do not exhibit the highly oscillatory fields produced by models derived from a least‐squares analysis with a truncated spherical harmonic series. Modeling the data by harmonic splines requires solving a square system of equations with dimension equal to the number of data. Too many data have been collected since the 1960s for this method to be practical. We produce almost optimally smooth models by the following method. Since each spline function for the optimal model corresponds to an observation location (called a knot), we select a subset of these splines with knots well‐distributed around the Earth’s surface. In this depleted basis we then find the smoothest model subject to an appropriate fit to all of the data. This reduces the computational problem to one comparable to least‐squares analysis while nearly preserving the optimality inherent in the original harmonic spline models.