Bridging the gap between airborne and spaceborne imaging spectroscopy for mountain glacier surface property retrievals

Observations of glacier albedo are sparse, despite being a first-order control on net solar radiation and melt rates. Recent and forthcoming airborne and satellite imaging spectrometer missions increase our ability to record changes in albedo over time and space but will require an understanding of uncertainties, especially in areas of rugged terrain. Some of these uncertainties arise from the timing of acquisition and resolution of digital elevation models (DEM) that are used for topographic correction. We use the fine-resolution data (2 m) from the Airborne Coastal Observatory (ACO), an imaging spectroscopy and lidar remote sensing platform, to quantify the expected error from other airborne and satellite scale platforms. First, we describe our retrieval workflow - the Imaging Spectroscopy Snow and Ice Algorithm (ISSIA) that yields spatial retrievals of broadband albedo (BBA), optical grain size, and radiative forcing by light absorbing particles (RFLAP) using the coincident 2 m lidar DEM over mountain glaciers. We use two acquisitions over Place Glacier, British Columbia, Canada that represent typical spring and summer conditions as a baseline from which synthetic datasets were generated to represent data from coarser-scale imaging spectroscopy platforms. Each synthetic dataset was motivated by specific airborne and satellite platform scenarios and quantifies four sources of error: (1) DEM spatial resolution and elevation accuracy, (2) coincident vs non-coincident DEM collection, (3) imagery spatial resolution, and (4) spectral resolution. Use of a 30 m DEM introduced the most error in steep mountain terrain, which primarily increased error in the BBA and RFLAP retrievals. The timing of DEM acquisition mattered most for high-resolution imagery and in areas surrounding the glacier toe where mass loss occurred more rapidly. Coarse-resolution imagery did not fully capture spatial variability in surface properties and topography, which lowered variance and increased overall errors. Coarsening of spectral resolution from 5 nm to 10 nm had minimal impact on surface property retrieval and associated errors were independent of DEM or imagery resolution. A simple energy analysis showed that surface albedo error, due to topographic correction with coarse and non-coincident DEMs, resulted in melt rate errors ranging from 6 to 20% for airborne and 11-20% for satellite scenarios.