Spatiospectral concentration in the Cartesian plane

Abstract

We pose and solve the analogue of Slepian's time-frequency concentration problem in the two-dimensional plane, for applications in the natural sciences. We determine an orthogonal family of strictly bandlimited functions that are optimally concentrated within a closed region of the plane, or, alternatively, of strictly spacelimited functions that are optimally concentrated in the Fourier domain. The Cartesian Slepian functions can be found by solving a Fredholm integral equation whose associated eigenvalues are a measure of the spatiospectral concentration. Both the spatial and spectral regions of concentration can, in principle, have arbitrary geometry. However, for practical applications of signal representation or spectral analysis such as exist in geophysics or astronomy, in physical space irregular shapes, and in spectral space symmetric domains will usually be preferred. When the concentration domains are circularly symmetric in both spaces, the Slepian functions are also eigenfunctions of a Sturm-Liouville operator, leading to special algorithms for this case, as is well known. Much like their one-dimensional and spherical counterparts with which we discuss them in a common framework, a basis of functions that are simultaneously spatially and spectrally localized on arbitrary Cartesian domains will be of great utility in many scientific disciplines, but especially in the geosciences.Comment: 34 pages, 7 figures. In the press, International Journal on Geomathematics, April 14th, 201

Accelerated West Antarctic ice mass loss continues to outpace East Antarctic gains

“…It is difficult to compare spatial resolution between modeling approaches without detailed considerations to quantify bias, variance, signal-to-noise ratios, and the effect of bandwidth selection Slobbe et al, 2012). Nevertheless, a detailed comparative study by Longuevergne et al (2010), carried out on the scale of hydrological basins, independently of our own efforts, was indicative of favorable behavior of the approaches based on Slepian functions.…”

“…The effort in constructing the Slepian functions is minimal, as they are efficiently computed via variety of numerical methods, both on the sphere (Simons and Dahlen, 2007) and in the plane (Simons and Wang, 2011). Compared to other gravitybased approaches, we do not a priori remove correlations between spherical-harmonic coefficients (Swenson and Wahr, 2006), and neither smooth nor project on a predefined basin structure (Wouters et al, 2008;King et al, 2012;Sasgen et al, 2012Sasgen et al, , 2013.…”

The spherical Slepian basis as a means to obtain spectral consistency between mean sea level and the geoid

“…Since the efficient generation of Slepian functions on domains of arbitrary geometry, whether in spherical or Cartesian (Simons and Wang 2011) coordinates, presents no intrinsic difficulties, we are here able to consider their use in the context of the work by , on which we build. Our main goal is to design a bandlimited Slepian basis for the ocean basins in spherical geometry and to evaluate the utility of this basis in solving the problem stated above, which is to derive suitable representations of altimetric MSL while maintaining spectral consistency with the gravimetric geoid.…”

On the robustness of estimates of mechanical anisotropy in the continental lithosphere: A North American case study and global reanalysis

“…For simplicity, we use square data patches or windows. Selective methods for regions of arbitrary description, such as irregular tectonic provinces, have also been developed, and may help avoid blending contrasting T e from different features into a single estimate (Simons and Wang, 2011). However, most workers to date have used square or circular windows, making a simple geometric patch a better test domain for the reliability of existing analyses.…”