Results of supervised classification of six Landsat 8 images acquired on: 25 May, 3 June, 22 June, 15 July, 4 August and 20 August 2014, covering glaciers in Hornsund fiord. Method of classification: Maximum Likelihood. The results show variability of snow cover areas in melting period of 2014 for glaciers located in Hornsund fiord and larger than 9 km2 (Körberbreen, Samarinbreen, Chomjakovbreen, Mendelejevbreen, Svalisbreen, Hornbreen, Storbreen, Kvalfangarbreen, Mühlbacherbreen, Paierlbreen and Hansbreen). For more information, please check: https://doi.org/10.3390/w9100804 Overview: Results of Maximum Likelihood classification of Landsat 8 images for analysed glaciers. Red - snow cover, yellow - glacier ice, black - debirs, grey - cloud cover.

The monitoring of the mass balance of the Werenskioldbreen (Wedel Jarlsberg Land, Spitsbergen, Svalbard) in the years 1999-2002 and 2009-2018. It is calculated on the base of 4 to 9 ablation stakes (depend on year). The mass balance is determined by conducting field surveys on floating calendar dates (floating-date system). Data have also been submitted to the World Glacier Monitoring Sevice (WGMS, https://wgms.ch)

The positions of the glacier termini in Hornsund are derived with very high frequency in the period 1991–2018. Over 230 multispectral and Synthetic Aperture Radar (SAR) data were used: LANDSAT 5, LANDSAT 7, LANDSAT 8, Terra ASTER, Alos AVNIR, SPOT 5, ERS-1, ERS-2, ENVISAT, Alos PALSAR, TerraSAR-X, TanDEM-X, and Sentinel-1. SAR data were used to detect any variability in the glacier front during the polar night. The satellite data were digitized manually to obtain the ice cliff position. Multispectral images were orthorectified and geocoded in PCI Geomatica and ArcGIS software. SAR data were usually provided at the SLC level, so that both radiometric and geometric corrections could be applied using the same methods, and with the same digital elevation model (2008 DEM SPOT developed by the IPY-SPIRIT Project; Korona et al., 2009). The SAR data were processed in BEAM (http://www.brockmann-consult.de/cms/web/beam). Sentinel data downloaded from the Sentinel’s Data Hub were already processed. Data not published.

Since 2010, continuous monitoring of the front zone of the Hansa Glacier with the use of Canon Eos 1000D photo cameras has been carried out (timelapse). Pictures in different periods of time were taken by 3 different cameras. Two cameras (106 and 107) were located on the Fugleberget slope and one (601) on the Baranowskiodden. The periods for which data are available and the interval of taking pictures are shown in Appendix 1.

Subglacial topography was derived from radio-echo sounding (RES) survey conducted in spring 2008 by the University of Silesia research team (M. Grabiec and J. Jania) in cooperation with the Institute of Geophysics Polish Academy of Sciences (D. Puczko) and the Maria Curie-Sklodowska University (G.Gajek). The profiles were acquired by the radar system equipped with 25 MHz unshielded antenna pulled behind snowmobile. Traces were recorded every 0.5 s, that translates into 1.5-2.0 trace-to-trace distance depending on the vehicle’s velocity. Traces were positioned by GNSS receiver working in differential mode with 3D accuracy ± 1m. In total over 100 km of RES profiles were acquired on Hansbreen, 66 km on Werenskioldbreen and 43 km on Renardbreen. RES data were processed using standard procedure including: DC-offset, time-zero adjustment, 2-D filter, amplitude correction, bandpass filtering and migration. Time-to-depth conversion applied average radio-wave velocity in glacier ice 16.4 cm ns-1 calculated based on CMP analysis performed on Hansbreen in the same season as the GPR profiling. The ice/bed interface was picked up semi-automatically with RMSE 5.3 ns (0.43 m) (Grabiec, 2017). Then the bedrock elevation data were interpolated over studied glaciers taking into account elevation of nonglaciated surroundings (Grabiec 2017) and bathymetry at the front of tidewater Hansbreen (Grabiec et al. 2012). Finally produced 100 m resolution DEMs are in UTM 33X WGS84 reference system. DEM of 300 m resolution is freely available. For 100 m resolution DEM please contact: mariusz.grabiec@us.edu.pl. The data were collected and processed under following projects: • IPY/269/2006 GLACIODYN The dynamic response of Arctic glaciers to global warming • UE FP7-ENV-2008-1 ice2sea Estimating the future contribution of continental ice to sea-level rise • PNRF-22-AI-1/07 AWAKE Arctic Climate and Environment of the Nordic Seas and the Svalbard – Greenland Area • 03/KNOW2/2014 KNOW Leading National Research Centre Reference: Grabiec M., Jania J., Puczko D., Kolondra L., and Budzik T., 2012: Surface and bed morphology of Hansbreen, a tidewater glacier in Spitsbergen. Polish Polar Research 38(2): 111-138. Grabiec M. 2017: Stan i współczesne zmiany systemów lodowcowych południowego Spitsbergenu w świetle badań metodami radarowymi. Wydawnictwo Uniwersytetu Śląskiego, 328 s. Decaux, L., Grabiec, M., Ignatiuk, D., and Jania, J. 2018: Role of discrete recharge from the supraglacial drainage system for modelling of subglacial conduits pattern of Svalbard polythermal glaciers, The Cryosphere Discuss., https://doi.org/10.5194/tc-2017-219, in review.

Snow depth data series contain records obtained by high frequency GPR on selected glaciers of Hornsund area (S Svalbard) since 2008. Currently the largest collection regards Hansbreen. Data for other glaciers are successively appended. The GPR survey on Hansbreen is regularly carried out approximately along the same tracks. Due to dynamically changing glacier surface topography influencing different survey abilities the some parts of profiles are modified in consecutive seasons. The total distance of Hansbreen profiles are as follows (Fig.1): 63.9 km (2008), 117,5 km (2011), 105,1 km (2013), 103,9 km (2014), 98,5 km (2015), 91,1 km (2016), 101,0 km (2017) and 108,4 km (2018).

• Discharge measurements were conducted from May/June/July to September/October in the years 2007–2012. Direct observation periods were 62, 51, 61, 40, 121, and 35 days in the consecutive years. The data set was obtained from CTD-DIVER DI 261 or Mini-Diver (Van Essen Instruments, Delft, The Netherlands) logger with barometric compensation by BaroDiver (Schlumberger, Houston, TX, USA) with 10-min intervals and flow velocities were measured with a SEBA F1 current meter (SEBA Hydrometrie GmbH, Kaufbeuren, Germany). Mean daily discharge and total runoff in the hydrologically active season was calculated on the basis of the 24h running average of the water level and a rating curve (Appendix 1). More details have been reported by Majchrowska E., Ignatiuk D., Jania J., Marszałek H., Wąsik M., 2015: Seasonal and interannual variability in runoff from the Werenskioldbreen catchment, Spitsbergen. Polish Polar Research vol. 36, no. 3, pp. 197–224. doi: 10.1515/popore−2015−0014 • Hydrological data for Werenskioldbreen - archive of the University of Silesia. The base contain mean daily discharge in the years: 1972, 1973, 1974, 1979, 1980, 1983, 1985, 1986, 1988, 1998, 2007, 2008, 2009, 2010, 2011, 2012, 2013, 2017 (18 hydrologically active seasons) – Appendix 2

The base contain: 1. List (*xls) of terrophotogrammetric photographs taken by member of Polish Expedition on Iceland in 1968 (zone of Skeidarárjökull). 2. Skaning above photographs with resolution of 2400 dpi, file format *.tif and *.jpg (400 dpi). All photographs were taken by means of instrument type Phototheodolite Zeiss 19/1318 on glass plates. Authors of photographs: Tadeusz Konysz. Fot better identification of all photographs was prepared presentations (in *.pptx and *.pdf formats), where photographs there are arranged in in common practise, classic form used for terrophotogrammetric stereoscopic photographs: A, AL, AR, B, BL, BR (the meaning of a letter: A - right photogrammetric station and photograph taken in normal case [90°]; B - left photogrammetric station and photograph taken in normal case [90°]; AL - photograph taken from right station in left direction; AR - photograph taken from right station in right direction; BL - photograph taken from left station in left direction; BR - photograph taken from left station in right direction). Other photographs (e.g. panoramic) has additional informations. All names of photogrammetric stations and name of objects (in Polish language) was taken from the author’s collection of data.

Dataset contains detailed measurements of physical features of seasonal snow cover, according to International Classification of Seasonal Snow on the Ground (Fierz et al. 2009). Fieldwork has been done during peak of accumulation on several glaciers in the region of Hornsund, Svalbard: Amundsenisen (2013) Hansbreen (three sites; 2010, 2012, 2013, 2014, 2015, 2017, 2018) Werenskioldbreen (two sites; 2013, 2015) Nannbreen (2013) Ariebreen (2014) Acknowledgements: Research Council of Norway, Arctic Field Grant 2013: Spatial distribution of snow cover and drainage systems on the glaciers on Wedel Jarlsberg Land (RiS ID: 6158); the National Science Centre PRELUDIUM 4: Role of meltwater from snow cover for supplying drainage systems of the Spitsbergen glaciers (2012/07/N/ST10/03784) References: Laska M., Luks B., Budzik T., 2016. Influence of snowpack internal structure on snow metamorphism and melting intensity on Hansbreen, Svalbard. Polish Polar Research, 37(2): 193–218. doi:10.1515/popore-2016-0012 Laska M., Grabiec M., Ignatiuk D., Budzik T., 2017. Snow deposition patterns on southern Spitsbergen glaciers, Svalbard, in relation to recent meteorological conditions and local topography. Geografiska Annaler, Series A: Physical Geography, 99(3): 262–287 doi:10.1080/04353676.2017.1327321 Laska M., Barzycka B., Luks B., 2017. Melting Characteristics of Snow Cover on Tidewater Glaciers in Hornsund Fjord, Svalbard. Water, 9(10), 804. doi:10.3390/w9100804

Radar satellite (SAR) images for Hornsund: ERS-1, ERS-2, ENVISAT, ALOS Palsar, TerraSAR-X, TandemX-1, acquired between 1992 and 2014. 210 archival SAR data were provided at the SLC level, so that both radiometric and geometric corrections were applied using the same methods, and with the same digital elevation model (2008 DEM SPOT developed by the IPY-SPIRIT Project; Korona et al., 2009). The SAR data were processed in BEAM (http://www.brockmann-consult.de/cms/web/beam).