Relationship between RADARSAT-2 Derived Snow Thickness on Winter First Year Sea-Ice and Aerial Melt-Pond Distribution using Geostatistics and GLCM Texture

A contemporary shift from primarily thicker, older multi-year sea ice (MYI) to thinner, smoother first-year sea ice (FYI) has been attributed to increased atmospheric and oceanic warming in the Arctic, with a steady diminishing of Arctic sea ice thickness due to a reduction of thick MYI compared to FYI. With an increase in FYI fraction, increased melting takes place during the summer months, exposing the sea ice to additional incoming solar radiation. With this change, an increase in melt pond fraction has been observed during the summer melt season. Prior research advocated that thin/thick snow leads to dominant surface flooding/snow patches during summer because of an enhanced ice-albedo feedback. For instance, thin snow cover areas form melt ponds first. Therefore, aerial measurements of melt pond fraction provide a proxy for relative snow thickness. RADARSAT-2 polarimetric SAR data can provide enhanced information about both surface scattering and volume scattering mechanisms, as well as recording the phase difference between polarizations. These polarimetric parameters can be computed that have a useful physical interpretation. The principle research focus is to establish a methodology to determine the relationship between selected geostatistics and image texture measures of pre-melt RADARSAT-2 parameters and aerially-measured melt pond fraction. Overall, the notion of this study is to develop an algorithm to estimate relative snow thickness variability in winter through an integrated approach utilizing SAR polarimetric parameters, geostatistical analysis and texture measures. Results are validated with test sets of melt pond fractions, and in situ snow thickness measurements. Preliminary findings show significant correlations with pond fraction for the standard deviation of HH and HV parameters at small incidence angles, and for the mean of the co-pol phase difference parameter at large incidence angles.