Evaluation of satellite remote sensing albedo retrievals over the ablation area of the southwestern Greenland ice sheet

MODerate resolution Imaging Spectroradiometer (MODIS) albedo products have been validated over spatially uniform, snow-covered areas of the Greenland ice sheet (GrIS) using the so-called single 'point-to-pixel' method. This study expands on this methodology by applying a 'multiple-point-to-pixel' method and examination of spatial autocorrelation (here using semivariogram analysis) by using in situ observations, high-resolution World View-2 (WV-2) surface reflectances, and MODIS Collection V006 daily blue-sky albedo over a spatially heterogeneous surfaces in the lower ablation zone in southwest Greenland. Our results using 232 ground based samples within two MODIS pixels, one being more spatial heterogeneous than the other, show little difference in accuracy among narrow and broad band albedos (except for Band 2). Within the more homogenous pixel area, in situ and MODIS albedos were very close (error varied from -4% to +7%) and within the range of ASD standard errors. The semivariogram analysis revealed that the minimum observational footprint needed for a spatially representative sample is 30 m. In contrast, over the more spatially heterogeneous surface pixel, a minimum footprint size was not quantifiable due to spatial autocorrelation, and far exceeds the effective resolution of the MODIS retrievals. Over the high spatial heterogeneity surface pixel, MODIS is lower than ground measurements by 4-7%, partly due to a known in situ undersampling of darker surfaces that often are impassable by foot (e.g., meltwater features and shadowing effects over crevasses). Despite the sampling issue, our analysis errors are very close to the stated general accuracy of the MODIS product of 5%. Thus, our study suggests that the MODIS albedo product performs well in a very heterogeneous, low-albedo, area of the ice sheet ablation zone. Furthermore, we demonstrate that single 'point-to-pixel' methods alone are insufficient in characterizing and validating the variation of surface albedo displayed in the lower ablation area. This is true because the distribution of in situ data deviations from MODIS albedo show a substantial range, with the average values for the 10th and 90th percentiles being -0.30 and 0.43 across all bands. Thus, if only single point is taken for ground validation, and is randomly selected from either distribution tails, the error would appear to be considerable. Given the need for multiple in-situ points, concurrent albedo measurements derived from existing AWSs, (low-flying vehicles (airborne or unmanned) and high-resolution imagery (WV-2)) are needed to resolve high sub-pixel variability in the ablation zone, and thus, further improve our characterization of Greenland's surface albedo. (C) 2017 Elsevier Inc. All rights reserved.