Exploratory Mapping of Blue Ice Regions in Antarctica Using Very High-Resolution Satellite Remote Sensing Data

Mapping spatiotemporal changes in the distribution of blue ice regions (BIRs) in Antarctica requires repeated, precise, and high-resolution baseline maps of the blue ice extent. This study demonstrated the design and application of a newly-developed semi-automatic method to map BIRs in the Antarctic environment using very high-resolution (VHR) WorldView-2 (WV-2) satellite images. We discussed the potential of VHR satellite data for the mapping of BIRs in the Antarctic environment using a customized normalized-difference blue-ice index (NDBI) method devised using yellow, green, and near-infrared spectral bands of WV-2 data. We compared the viability of the newly developed customized NDBI approach against state-of-the-art target detection (TD), spectral processing (SP) and pixel-wise supervised (PSC) feature extraction (FE) approaches. Four semi-automatic FE approaches (three existing plus one newly developed) consisting of 16 standalone FE methods (12 existing + four customized) were evaluated using an extensive quantitative and comparative assessment for mapping BIRs in the vicinity of Schirmacher Oasis, on the continental Antarctic coastline. The results suggested that the customized NDBI approach gave a superior performance and the highest statistical stability when compared with existing FE techniques. The customized NDBI generally rendered the lowest level of misclassification (average RMSE = 654.48 ± 58.26 m2), followed by TD (average RMSE = 987.81 ± 55.05 m2), SP (average RMSE = 1327.09 ± 127.83 m2) and PSC (average RMSE = 2259.43 ± 115.36 m2) for mapping BIRs. Our results indicated that the use of the customized NDBI approach can greatly improve the semi-automatic mapping of BIRs in the Antarctic environment. This study presents the first refined map of distribution of BIRs around the Schirmacher Oasis. The total area of blue ice in the study area was estimated to be 106.875 km2, approximately 61% of the study area. The WV-2 derived BIR map area presented in this study locally refined the existing BIR map derived using Landsat Enhanced Thematic Mapper Plus (ETM+) and the Moderate Resolution Imaging Spectroradiometer (MODIS)-based mosaic of Antarctica (MOA) dataset by ~31% (~33.40 km2). Finally, we discussed the practical challenges and future directions in mapping BIRs across Antarctica.