Spectral and spatial decomposition of lithospheric magnetic field models using spherical Slepian functions

Abstract

Global magnetic field models are typically expressed as spherical-harmonic expansion coefficients. Slepian functions are linear combinations of spherical harmonics that produce new basis functions, which vanish approximately outside chosen geographical boundaries but also remain orthogonal within the spatial region of interest. Hence, they are suitable for decomposing spherical-harmonic models into portions that have significant magnetic field strength only in selected areas. Slepian functions are spatio-spectrally concentrated, balancing spatial bias and spectral leakage. Here, we employ them as a basis to decompose the global lithospheric magnetic field model MF7 up to degree and order 72, into two distinct regions. One of the resultant fields is concentrated within the ensemble of continental domains, and the other is localized over its complement, the oceans. Our procedure neatly divides the spectral power at each harmonic degree into two parts. The field over the continents dominates the overall crustal magnetic field, and each region has a distinct power-spectral signature. The oceanic power spectrum is approximately flat, while that of the continental region shows increasing power as the spherical-harmonic degree increases. We provide a further breakdown of the field into smaller, non-overlapping continental and oceanic regions, and speculate on the source of the variability in their spectral signatures.