Domain Editor-in-Chief: Jody W. Deming; Department of Biological Oceanography, University of Washington, US
Associate Editor: Mary-Louise Timmermans; Yale University, US

1. Introduction

As incident solar radiation increases and the air warms during the onset of summer, the Arctic sea ice cover undergoes a spectacular transformation from a high albedo (>0.8), snow-covered surface during the winter months to a lower albedo surface dominated by melt ponds and open leads. Perovich and Polashenski (2012) divided the surface changes into distinct stages beginning with the wetting of the winter snow layer and the appearance of ponds as meltwater accumulates on the surface. Several weeks after the onset of melt and rapid increase in pond coverage, melt pond drainage occurs when the higher ice permeability and number of open macroscopic flaws during this stage allow water to be lost to the ocean. As a result of drainage, pond coverage drops significantly on first-year ice, where melt ponds are more extensive and shallower, but only slightly on multiyear ice, where thicker and more deformed ice can retain more surface meltwater. Following this stage, pond coverage continues to increase until the seasonal maximum is reached. The variability in the rate of progression through these stages of melt, however, has received less attention likely due to the lack of direct measurements of pond coverage.

A recent study by Schröder et al. (2014), using a melt pond model with more realistic process representations that include surface topography development (Scott and Feltham, 2010), showed that spring melt pond coverage is a skillful predictor of summer minimum ice extent (SMIE). This skill is thought to be derived from improved simulation of the positive albedo feedback in response to pond formation: reduced surface albedo leads to increased absorption of solar radiation, increased ice melt, and enhanced melt pond formation. These results suggest that variability in the spatial and temporal development of surface albedo, an expression of the initial state of the ice cover and the atmosphere/ocean forcing during the summer, could be useful indicators for SMIE forecasts. Despite this potential, and increased interest in more accurate seasonal sea ice forecasts (Eicken, 2013), we are not aware of any published investigation that utilizes retrieved surface albedo in a statistical forecast of SMIE. This may be due to the infrequent satellite retrievals of surface albedo due to clouds during the melt season, solar elevations at polar latitudes, and the latency of these retrievals for use in forecasts. Instead, passive microwave brightness temperatures have been used to provide signals of melt (e.g., onset dates) for the forecast of SMIE (e.g., Petty et al., 2017) rather than measures from albedo changes.

Specular returns from the sea ice cover in CryoSat-2 (CS-2) waveforms are robust indicators of smooth open water or thin ice in leads and melt ponds within the radar-altimeter footprint. This type of distinctive return, easily identified in altimetric waveforms, has been used effectively in ice/water discrimination (e.g., Peacock and Laxon, 2004; Laxon et al., 2013; Armitage and Davidson, 2014) and the identification of local sea levels for the calculation of sea ice freeboard (e.g., Kwok and Cunningham, 2015). Hence, we expect that changes in the density of specular returns in CS-2 data can be compared to surface albedo, as both are dominated by development of melt pond coverage during the melt season.

In this paper, we examine the fractional coverage of specular returns (FSR) as an indicator of open water coverage on the ice cover. The utility of FSR is demonstrated by comparing the relative effectiveness of using FSR and surface albedo for statistical forecasting of summer minimum ice extent (SMIE). After describing the datasets used (Section 2), we describe the computation of FSR, the comparison of the time-varying spatial patterns of FSR with the development of area-averaged surface albedo during summer, and the skill of using FSR and surface albedo to the forecast of the SMIE (Section 3), ending with discussion and conclusions (Section 4).

2. Data Description

2.1. Specular returns in CryoSat-2 waveforms (2011–2017)

We used radar data acquired by the Synthetic Aperture Radar (SAR) Interferometric Radar ALtimeter (SIRAL) instrument on CryoSat-2 (CS-2) (Wingham et al., 2006). The altimeter waveforms (Release C) are available through ESA’s data portal (URL: The synthesized footprint of SIRAL in SAR operation is nominally 0.31 km by 1.67 km in the along- and across-track directions. Specular and near-specular returns within the radar footprint are identified using the combined peak power and width at the half-power point of the first return; the details of this approach are described by Kwok and Morison (2016), where they used the procedure to retrieve the heights of the sea surface over ice-covered oceans. Specular and near-specular returns from areas of open water much smaller than the radar footprint typically dominate and mask out diffuse returns from the surrounding ice cover (Drinkwater, 1991), making SIRAL returns highly sensitive to the presence of surface water. Section 3.1 provides a more detailed discussion of the sources of specular returns during early melt.

2.2. Surface albedo fields (2011–2017)

Weekly-mean surface albedo fields of Arctic sea ice (on a 15-km map grid), between 2011 and 2017, are from EUMETSAT’s Satellite Application Facility on Climate Monitoring (CM SAF CLARA-SAL data set). The retrievals are based on the aggregation of radiance data acquired by the Advanced Very High Resolution Radiometer (AVHRR) instrument on board polar-orbiting NOAA and Metop satellites (Karlsson et al., 2012). Raw radiances are corrected for atmospheric, radiometric, and topographic effects and are expanded into hemispherical spectral albedos before shortwave albedo is obtained using a narrow-to-broadband conversion. Expected uncertainty in the albedo estimates is around 5% (or 0.05), and retrievals of temporal mean area-averaged albedo derived from only a few observations are vulnerable to errors in cloud masking. The reader is referred to Karlsson et al. (2012) for a more detailed description of this product.

2.3. Sea ice extent and high resolution imagery

The sea ice extent dataset was taken from the National Snow and Ice Data Center (Fetterer et al., 2016). Declassified high-resolution visible imagery of the Arctic Ocean, described by Kwok (2014), is from the Global Fiducial Library – GFL (URL:; accessed in January 2017, website currently under revision). These high spatial resolution (~1 m) pan-chromatic images easily resolve melt ponds on the surface and cover the period between 2011 and 2014. The GFL is managed by the U.S. Geological Survey under the National Civil Applications Program.

3. Methods and Results

In this section, we discuss the sources of specular returns in the CS-2 data during early melt, followed by a description of the computation and interpretation of FSR. The time-varying spatial patterns of FSR are then compared to the development of surface albedo during summer. Last, we assess the skill of FSR and surface albedo for forecasts of the summer minimum ice extent (SMIE). The objective is to show the relationships between FSR, albedo, and SMIE.

3.1. Sources of specular returns

During early melt, the period of interest in this work, the surface is composed of open leads and a melting snow cover with increasing coverage of ponds. Specular returns from open leads have been used as a source of sea surface height within the sea ice cover during winter and summer (e.g., Kwok and Morison, 2016), and also during early melt. The smooth open water surfaces of melt ponds are another source of specular returns as melt progresses and likely the main contributor to the large increases in specular returns seen in early melt. The melting snow surfaces tend to remain rough compared to water surfaces or the wavelength of the radar: ~2 cm, and thus their returns are expected to remain diffused and hence not specular. Away from the marginal ice zone at this stage of melt, the contribution of open leads to the increases in specular returns is not expected to be significant as the break-up of the ice cover does not typically occur until later in the melt season (late July to August). However, in our analysis, specular returns from open leads and melt ponds are indistinguishable from each other; we are not aware of any investigation in the published literature that has characterized the evolution of the relative coverage of open water in leads and ponds after the onset of melt.

In CS-2 data, we used specular and near-specular returns that span a range of 25 dB (~300 times in return power) above the surrounding ice cover (Kwok and Morison, 2016), which translates roughly into a potential sensitivity to specular areas covering less than ~40 m2 (i.e., 300 times smaller than the first Fresnel Zone area of ~12,468 m2 calculated using the radar altitude and frequency). Even though the specular returns are sensitive to small areas of open water, these returns could be from individual melt ponds (with scales of meters) as well as from the surfaces that contain a collection of ponds and leads.

3.2. Fractional coverage of specular returns (FSR)

We define FSR as the fractional coverage of specular returns within 25-km segments of CS-2 ground tracks, which are sampled by ~80 contiguous synthetic aperture footprints at a spacing of ~300 m. FSR varies between 0 and 1. Even though the radar returns are highly sensitive to the presence of open water or thin ice, FSR is not an accurate indicator of ice concentration because, as mentioned above, specular returns from areas of open water much smaller than the radar footprint tend to dominate and mask diffuse scattering from the surrounding ice cover. Hence, FSR tends to underestimate ice concentration or overestimate open water fraction.

The maps in Figure 1a show gridded (50 km) composites of FSR constructed from all 25-km along-track segments for consecutive 15-day periods (~210 CS-2 orbits) between June and July for seven years (2011–2017). North-south striping in the composites (see, e.g., 15–30 June 2014) are due to combining data that sampled the surface at different times, as melt progresses rapidly poleward within the 15-day intervals. For each year, the space-time progression of melt is evident in the FSR composites. In early June, FSR is low (~0.1) over most of the central Arctic but slightly higher (~0.2) over the thinner, more deformable ice, with lower ice concentration occupying the peripheral seas (Beaufort, Chukchi, East Siberian and Laptev Seas). As the snow cover melts and ponds develop when the air temperature climbs with the onset of summer, the number of specular returns increases dramatically in the lower latitudes (FSR > 0.8). Given the temporal coarseness of the composites, these fields are not able to capture the finer details of the initial stages of surface melt (e.g., decreases in pond coverage associated with drainage) described by Perovich and Polashenski (2012). Increasing FSR progresses from south to north as solar illumination increases and the Arctic warms, peaking after 15 July with high FSR observed over large regions of the Arctic. FSR begins to return to winter values at the end of the summer, progressing southwards as the temperature falls.

Figure 1 

Gridded 15-day composites of specular returns fraction and albedo between June and July (2011–2017). Correlation between the two composites of (a) specular returns fraction (FSR) and (b) albedo (α) are shown in the bottom left corner of the images in (b). After mid-July (below the dashed line), FSR saturates and the correspondence between the two composites is reduced. DOI:

This seasonal evolution is also illustrated in the seasonal cycle of area-averaged FSR (FSR) within two latitudinal bands: 70–80°N and 80–90°N (Figure 2c) to contrast the progression of melt between the lower and higher latitude Arctic Ocean. As pond coverage increases steadily, higher FSR (>0.8) is seen over most of the Arctic Ocean by the latter half of July, with the higher latitudes lagging the lower latitudes by ~16 days (average over 2011–2017). As the ice cover remains relatively compact during early melt (i.e., low lead fraction), increases in surface water coverage due to melt ponds dominate (see also examples of imagery in Figure 3) while the ice concentration remains high (average ice concentration within the ice in June for seven years is ~96%). By the early half of September, FSR begins its downward trend in the central Arctic and expands southwards as surface meltwater disappears as the air temperature falls.

Figure 2 

Concurrent changes in specular returns fraction and albedo during the melt season (2011–2017). Specular returns fraction (FSR) and albedo (α) for (a) multiyear ice fraction (MYf) < 0.7 and (b) MYf ≥ 0.7; black horizontal lines indicate FSR threshold of 0.5. (c) Area-averaged FSR and α over two latitudinal bands: 70–80°N (black lines) and 80–90°N (red lines); dashed horizontal line indicates FSR threshold of 0.5. (d) Dates (day of year, DOY) when FSR (colored vertical lines) crosses 0.5 and the albedo decreases to 0.5 (colored dots) in two latitudinal bands (70–80°N in black and 80–90°N in red), and the length of the melt season defined by the two crossings FSR. DOI:

3.3. Comparison of time-varying fields of FSR and surface albedo

To assess the similarity between FSR and albedo fields, we compared the time-varying spatial patterns and the interannual variability of the FSR composites with fields of surface albedo. For this, we constructed gridded albedo composites (Figure 1b) with the same temporal sampling as the FSR composites, i.e., consecutive 15-day periods for seven years (2011–2017). Only the June and July composites are shown in Figure 1. The albedo composites do not cover the month of September due to reduced retrievals as the sun begins to set in high Arctic latitudes. In contrast, the FSR composites span the entire melt period because the CS-2 radar acquisitions are not sensitive to clouds or solar illumination of the surface, though they are less useful in winter.

Although there are differences, the large-scale correspondence between the spatial patterns of the FSR and albedo composites is clear during June and early July (spatial correlations of the 15-day fields are shown in Figure 1b). During early melt, these fields are negatively correlated where increases in FSR correspond to decreases in albedo due to higher pond coverage associated with melt. It is not expected that these fields should have exact correspondence. Differences between the 15-day FSR and albedo composites can be attributable to: 1) the sparsity in altimeter footprints and albedo retrievals during summer; 2) the lack of spatial and temporal coincidence between the two measurements; and 3) the saturation of FSR at low albedos (at ~0.6). Past mid-July, spatial correlations with the albedo composites are reduced and become spurious (see correlation values in Figure 1) as the FSR composites become largely saturated as FSR becomes insensitive to further increase in melt pond coverage. More quantitatively, Figure 2a and 2b compares the evolution of FSR and albedo over seasonal (i.e., multiyear ice fraction [MYf] < 0.7) and multiyear ice (MYf ≥ 0.7), and Figure 2c shows their seasonal cycle in two latitudinal bands: 70–80°N and 80–90°N. FSR climbs to >0.9 in seasonal ice but only to ~0.8 in multiyear ice. Similarly, area-averaged albedo (α¯) falls to 0.6 in seasonal ice, and only to ~0.5 in multiyear ice. These behaviors are consistent with lower melt pond coverage over older ice in higher latitudes (i.e., lower FSR and higher albedo).

FSR and albedo can also be used to identify dates delineating the progress of the melt season. Figure 2d shows the dates at which FSR climbs/falls to 0.5, and the dates at which α¯ falls to 0.5 in the two latitudinal bands. The selection of 0.5 as a threshold is somewhat arbitrary but is discussed further below. These two dates from FSR can be used to calculate the duration over which FSR remains above 0.5 (we designate as the duration as D0.5). Figure 2d shows that 2013 has both the shortest D0.5 for both latitudinal bands and the lowest excursions of FSR and albedo during the summer in the 7-year record.

A closer examination of 2013 (composites in Figure 1a and b) shows that there were large areas of lower FSR and higher albedo over the central Arctic ice cover throughout much of the 2013 melt season, indicative of lower melt pond coverage compared to other years. The notably higher albedo in 2013 in late June through August (compared to other years) can be seen in the high-resolution visible-band imagery in Figure 3: there was a complete absence of melt ponds on 20 and 21 August (in 2013) at a location just north of Ellesmere Island, whereas imagery from 22 August 2011, 16 August 2012 and 18 August 2014 all show at least partial melt pond coverage at the same location (map in Figure 3 shows location of the imagery).

Figure 3 

Available high-resolution optical imagery of sea ice north of the Canadian Arctic Archipelago and Greenland. The map (a) shows the location of the imagery in (b). Each high-resolution (1 m) image is 2 km on a side. The imagery used here is described by Kwok (2014) and available through the Global Fiducials Library at USGS (URL: DOI:

The impacts of the anomalous surface melt conditions on the ice cover in summer 2013 were discussed by Kwok and Cunningham (2015), Tilling et al. (2015), and Kwok (2015). In addition to the cooler summer air temperatures in 2013 (compared to 2011–2014 averages), Kwok (2015) also noted anomalously strong mechanical convergence of the ice cover against the Canadian Arctic Archipelago (CAA) due to persistent on-shore winds, which resulted in a net area decrease of ~23% between June and August. Cooler air temperatures and closing of the ice cover reduced ice melt and the impact of openings on air temperatures (i.e., lower albedo and reduced ice-albedo feedback) in the ice cover. The estimates of FSR and albedo in this region (Figure 1) are consistent with the reported behavior. Also of note is the joint behavior of FSR and α¯ in 2017, which reveals a cooler eastern Arctic during most of the summer. The above comparison with albedo suggests that FSR is sensitive to surface melt and is a useful geophysical parameter.

3.4. Summer minimum sea-ice extent (SMIE): Regression analysis

Here, we compare the relative utility of FSR and surface albedo by examining the relationship between FSR, albedo and the Arctic SMIE, then assess the skill of using FSR and α to forecast SMIE. We considered two models wherein the SMIE is dependent on threshold-crossing dates derived from the composites of FSR and surface albedo (α):




The terms ΦSR.L(x) and ΦSR.H(x) are the dates at which FSR climbs above the threshold x in the two latitudinal bands, designated here by L (lower, corresponding to the 70–80°N band) and H (higher, corresponding to the 80–90°N band). Similarly, we define ΛL(y) and ΛH(y) as the dates at which α¯ falls below the threshold y in the same two latitudinal bands L and H; a, b, c, d, e, and f are linear coefficients. These two linear models are used to examine the relationship between SMIE and the progress of surface melt as informed by the crossing dates of the FSR and α¯ thresholds (i.e., x and y) for the two latitudinal bands during early melt.

Figure 4a and b summarizes results from the regression of SMIE against the independent variables in Equations (1) and (2) for the years 2011–2017. To assess the sensitivity of the regression model, we varied the FSR and α¯ thresholds (i.e., x and y) for determining the crossing dates ΦSR.L(x), ΦSR.H(x), ΛL(y), and ΛH(y). The regression analysis was carried out for x = 0.4, 0.5, and 0.6 and for y = 0.5 and 0.6. As a measure of the fraction of explained variance in the regression analysis, we show both the squared correlation ρ2 as well as the adjusted squared correlation ρ2¯, where ρ2¯ is always less than ρ2. The adjusted measure takes into account automatic increases in ρ2 when more than one explanatory variable is added (Draper and Smith, 1998), and is based on the number of explanatory terms relative to the number of data samples. For the above FSR and α¯ thresholds, the best regression results (i.e., estimated SMIE; Figure 4) were obtained with x = 0.5 and y = 0.6, where the models explain ~72/57% and 93/89% (ρ2/ρ¯2) of the variance in SMIE, with standard errors of 0.37 × 106 km2 and 0.19 × 106 km2, respectively (the crossing dates for these two cases can be found in the Appendix). For the seven years (2011–2017), the mean crossing dates are 10 June (day of year = 161) and 26 June (177) for ΦSR.L(0.5) and ΦSR.H(0.5), and 5 June (156) and 19 June (170) for ΛL(0.4) and ΛH(0.4). Not unexpectedly, in all cases examined, the linear model using the albedo crossing dates was able to explain higher percentages of the SMIE variance than the one using the FSR crossing dates, as the albedo fields provide denser coverage of the surface and are a more direct representation of the pond coverage. Figure 4e shows the crossing dates for these two cases and compares the regression results with the observed SMIE.

Figure 4 

Regression of Arctic minimum sea ice extent against observed crossing dates over two latitudinal bands. Shown are correlation results (a–d) and comparisons (e) of observed summer minimum sea ice extent with those forecast using different thresholds for fractional coverage of specular returns (FSR, threshold x) and surface albedo (α¯, threshold y). In (a), ΦSR.L(x) and ΦSR.H(x) refer to the dates when the area-averaged FSR (FSR) for 2011–2017 climbs above threshold x (see Equation 1) in the lower (L, 70–80°N) and higher (H, 80–90°N) latitudinal bands. Squared correlation/adjusted squared correlation (first row) and standard error (second row) are given for three thresholds; underlined values are the statistical significance (p values) of selected regression analyses. In b), same, except that ΛL(y) and ΛH(y) indicate the dates when surface albedo (α¯) falls below threshold y (see Equation 2) in latitudinal bands L and H, for two thresholds. In c) and d), same as a) and b), respectively, except using only crossings from 2011–2016. Additional row (in italics) shows the difference between the forecast and observed ice extent (IE, September minimum) in 2017 (of 4.636 × 106 km2). In upper e), day of year for the different threshold crossings in 2011–2017 is shown. Blue lines/triangles show ΛL(0.6) and ΛH(0.6), red lines/circles show ΦSR.L(0.5) and ΦSR.H(0.5), and thin and thick lines are the crossing dates at the lower (L, 70–80°N) and higher (H, 80–90°N) latitudinal bands, respectively. In lower e), observed ice extent (black line) is compared with those estimated using regression coefficients from crossing dates of 2011–2017 (colored lines, corresponding to information in a) and b) in the colored line boxes) and from only 2011–2016 (colored dashed lines, corresponding to information in c) and d) in colored dashed-line boxes). DOI:

3.5. SMIE forecasts with regression models

To examine the effectiveness of the observed crossing dates in forecasting SMIE, we forecast the 2017 SMIE with linear coefficients from a regression analysis using only the crossing dates from 2011–2016 (i.e., excluding 2017). Similar to above, the best regression results for the six years were obtained with x = 0.5 and y = 0.6 (Figure 4c and 4d). The two models explain ~79/65% and 93/88% (ρ2/ρ¯2) of the variance in SMIE, with standard errors of 0.37 × 106 km2 and 0.22 × 106 km2, comparable to the results using all seven years of crossing dates. With the observed crossing dates from 2017 and Equation 1, the regression coefficients from 2011–2016 were able to forecast the 2017 SMIE with residuals of –0.06 × 106 km2 (2% of the 2017 SMIE) and –0.17 × 106 km2 (4%), respectively.

We repeated this analysis using only crossing dates from 2011 to 2015 (i.e., excluding those from 2016 and 2017). Again, the best fit for the five years (not shown here) was obtained with x = 0.5 and y = 0.6, where the regression models explain ~88/56% and 98/93% (ρ2/ρ¯2) of the SMIE variance, and were able to forecast the 2016 SMIE with residuals of 0.06/0.31 × 106 km2, and the 2017 SMIE with residuals of 0.5/–0.07 × 106 km2, for FSR and α¯, respectively. These comparisons suggest that during early melt (up to ~mid-July), using FSR from CS-2 altimeter returns as indicators of open water coverage is comparable to using surface albedo to the extent measured by the explained variance. Further, the low SMIE residuals obtained from the forecasts demonstrate that threshold crossing dates from FSR or α¯ have potential for SMIE forecast.

4. Discussion and Conclusions

In this paper, we examined the relationship between the fractional coverage of specular returns (FSR) derived from the CryoSat-2 radar altimeter, surface albedo of the Arctic ice cover, and the summer minimum ice extent (SMIE). Specular reflections in altimeter returns are sensitive indicators of open water in sea ice leads and melt ponds within the Arctic sea ice cover. This is especially evident in the large concurrent increases in specular returns in CS-2 waveforms and decreases in surface albedo as the sea ice cover transitions from a high albedo snow-covered surface to a lower albedo surface dominated by ponds during early melt (Figure 1). Our analysis addressed the geophysical utility of FSR by comparing the relative effectiveness of using FSR with surface albedo for forecasting SMIE; results show that both FSR and surface albedo can be used to forecast SMIE, as summarized below.

Regression analysis of the area-averaged FSR and albedo time series, between 2011 and 2017, shows that ~72% and ~93% of the variance in SMIE can be explained by the dates when area-averaged FSR climbs to 0.5 and the area-averaged albedo falls to 0.6 in two latitudinal bands (70–80°N and 80–90°N). The results suggest that these threshold-crossing dates are effective in capturing the temporal and poleward progress of surface melt, and are useful forecasting SMIE. Excluding 2017 from the analysis, we were able to forecast the 2017 SMIE with residuals of 0.06 (~2% of the total extent) and –0.17 × 106 km2 (~4%) (see Figure 4). These residuals are comparable to the top two statistical forecasts of SMIE in June (see SIPN report in Further, excluding both 2016 and 2017 in our analysis yielded forecast residuals that are less than 0.5 × 106 km2 for both the 2016 and 2017 SMIE. Even though the CS-2 record is only seven years, the analyses here show the relative robustness of the regression models, whether using the full record or only a truncated portion of it.

In an examination of the summer ice cover of 2013, both the FSR and albedo fields (Figure 1) captured the lower pond coverage in the central Arctic that summer, as also evident in the later FSR and α¯ threshold crossing dates. The surface conditions were associated with cooler air temperatures and a convergent ice cover (noted by Kwok, 2015; Kwok and Cunningham, 2015; Tilling et al., 2015). The reduced summer melt left a thicker ice cover at the end of summer that was significant in CS-2 thickness estimates. Our analysis demonstrates that both the FSR and albedo fields can be used to inform understanding of SMIE and subsequent winter ice thickness variability.

While the satellite retrievals of surface albedo may be prone to uncertainties due to clouds in the melt season and their dependence on solar illumination at polar latitudes, the FSR composites from CS-2 radar acquisitions are insensitive to clouds and solar illumination. Given this insensitivity and the ~2 day latency of Near Real-Time CS-2 products, both FSR from CS-2 and surface albedo (with latency of ~1 week) during early melt may be useful parameters to include in forecasts of SMIE, when they are available around the end of June. We also note that even though large-scale change in albedo during summer melt is a characteristic feature of the sea ice surface, the available satellite-derived albedo fields have not been used directly in statistical forecasts of SMIE to date. Here, we used simple regression models to forecast end-of-summer ice extent and demonstrate the utility of FSR and albedo; the results also suggest that a more complex model, exploiting the spatial nature of the satellite data, could be explored to investigate the spatial variability of end-of-summer ice coverage.

Data Accessibility Statement

The CryoSat-2 data set is from the European Space Agency data portal (URL:, the albedo dataset is from the EUMETSAT Satellite Application Facility on Climate Monitoring (CM SAF, URL:, and the high-resolution visible imagery is from the Global Fiducial Library (URL:, accessed in January 2017, website currently under revision).