A method for constructing directional surface wave spectra from ICESat-2 altimetry

Sea ice is important for Earth's energy budget as it influences surface albedo and air-sea fluxes in polar regions. On its margins, waves heavily impact sea ice. Routine and repeat observations of waves in sea ice are currently lacking, and therefore a comprehensive understanding of how waves interact with sea ice and are attenuated by it is elusive.In this paper, we develop methods to separate the two-dimensional (2D) surface wave spectra from sea-ice height observations made by the ICESat-2 (IS2) laser altimeter, a polar-orbiting satellite. A combination of a linear inverse method, called generalized Fourier transform (GFT), to estimate the wave spectra along each beam and a Metropolis-Hastings (MH) algorithm to estimate the dominant wave's incident angle was developed. It allows us to estimate the 2D wave signal and its uncertainty from the high-density, unstructured ATL03 ICESat-2 photon retrievals. The GFT is applied to re-binned photon retrievals on 25 km segments for all six beams and outperforms a discrete Fourier transform (DFT) in accuracy while having fewer constraints on the data structure.The MH algorithm infers wave direction from beam pairs every 25 km using coherent crests of the most energetic waves. Assuming a dominant incident angle, both methods together allow a decomposition into 2D surface wave spectra with the advantage that the residual surface heights can potentially be attributed to other sea-ice properties. The combined GFT-MH method shows promise in routinely isolating waves propagating through sea ice in ICESat-2 data. We demonstrate its ability on a set of example ICESat-2 tracks, suggesting a detailed comparison against in situ data is necessary to understand the quality of retrieved spectra.