2007/05/22 by F. A. Dahlen, Frederik J. Simons, Frederik J Simons · 2 citations
Environmental Science · Physics and Astronomy · #Cosmology and Gravitation Theories #Radio Astronomy Observations and Technology #Soil Geostatistics and Mapping #astro-ph
paper · pdf · doi:10.1111/j.1365-246x.2008.03854.x
Submitted to the Geophysical Journal International
arxiv created 2007/05/22 · openalex publication_date 2008/08/22 · arxiv updated 2009/12/01 · openalex created_date 2016/06/24 · openalex updated_date 2026/07/28
We address the problem of estimating the spherical-harmonic power spectrum of a statistically isotropic scalar signal from noise-contaminated data on a region of the unit sphere. Three different methods of spectral estimation are considered: (i) the spherical analogue of the one-dimensional (1-D) periodogram, (ii) the maximum likelihood method, and (iii) a spherical analogue of the 1-D multitaper method. The periodogram exhibits strong spectral leakage, especially for small regions of area A≪ 4π, and is generally unsuitable for spherical spectral analysis applications, just as it is in 1-D. The maximum likelihood method is particularly useful in the case of nearly-whole-sphere coverage, A≈ 4π, and has been widely used in cosmology to estimate the spectrum of the cosmic microwave background radiation from spacecraft observations. The spherical multitaper method affords easy control over the fundamental trade-off between spectral resolution and variance, and is easily implemented regardless of the region size, requiring neither non-linear iteration nor large-scale matrix inversion. As a result, the method is ideally suited for most applications in geophysics, geodesy or planetary science, where the objective is to obtain a spatially localized estimate of the spectrum of a signal from noisy data within a pre-selected and typically small region.