Articles | Volume 14, issue 5
Earth Syst. Sci. Data, 14, 2487–2500, 2022
Earth Syst. Sci. Data, 14, 2487–2500, 2022
Data description paper
01 Jun 2022
Data description paper | 01 Jun 2022

A decade of glaciological and meteorological observations in the Arctic (Werenskioldbreen, Svalbard)

A decade of glaciological and meteorological observations in the Arctic (Werenskioldbreen, Svalbard)
Dariusz Ignatiuk1, Małgorzata Błaszczyk1, Tomasz Budzik1, Mariusz Grabiec1, Jacek A. Jania1, Marta Kondracka1, Michał Laska1, Łukasz Małarzewski1, and Łukasz Stachnik2 Dariusz Ignatiuk et al.
  • 1Faculty of Natural Sciences, University of Silesia in Katowice, Katowice, 40-007, Poland
  • 2Faculty of Earth Sciences and Environmental Management, University of Wrocław, Wrocław, 50-137, Poland

Correspondence: Dariusz Ignatiuk ( and Marta Kondracka (


The warming of the Arctic climate is well documented, but the mechanisms of Arctic amplification are still not fully understood. Thus, monitoring of glaciological and meteorological variables and the environmental response to accelerated climate warming must be continued and developed in Svalbard. Long-term meteorological observations carried out in situ on glaciers in conjunction with glaciological monitoring are rare in the Arctic and significantly expand our knowledge about processes in the polar environment. This study presents glaciological and meteorological data collected for 2009–2020 in southern Spitsbergen (Werenskioldbreen). The meteorological data are composed of air temperature, relative humidity, wind speed, short-wave and long-wave upwelling and downwelling radiation on 10 min, hourly and daily resolution (2009–2020). The snow dataset includes 49 data records from 2009 to 2019 with the snow depth, snow bulk density and snow water equivalent data. The glaciological data consist of seasonal and annual surface mass balance measurements (point and glacier-wide) for 2009–2020. The paper also includes modelling of the daily glacier surface ablation (2009–2020) based on the presented data. The datasets are expected to serve as local forcing data in hydrological and glaciological models as well as validation of calibration of remote sensing products. The datasets are available from the Polish Polar Database (, last access: 24 May 2022) and Zenodo (, Ignatiuk, 2021a;, Ignatiuk, 2021b).

1 Introduction

Long-term meteorological observations carried out in situ on glaciers in conjunction with glaciological monitoring are rare in the Arctic and may be used to expand knowledge about processes in the polar environment. Terrestrial meteorological monitoring alone does not always adequately address the needs of numerical modelling as well as validation and calibration of satellite products regarding glaciers (Pellicciotti et al., 2014; Gabbi et al., 2014). The warming of the Arctic climate is well documented, but the mechanisms of Arctic amplification are still not fully understood (IPCC, 2019). Both, climate and ocean variables have fluctuated in Svalbard in recent decades (Nordli et al., 2014; Schuler et al., 2014; Isaksen et al., 2016; Vikhamar-Schuler et al., 2016; Walczowski et al., 2017; Osuch and Wawrzyniak, 2017; Førland et al., 2020; Wawrzyniak and Osuch, 2020), which causes progressive and ongoing changes in the cryosphere (Błaszczyk et al., 2013; Wawrzyniak et al., 2016; Box et al., 2019; Grabiec et al., 2018; Nuth et al., 2019; van Pelt et al., 2019; Schuler et al., 2020; Błaszczyk et al., 2021). According to the data published in SIOS data access portal (, last access: 24 May 2022) and the meteorological bulletin “Spitsbergen-Hornsund” (, last access: 24 May 2022), 2020 was the year with the warmest summer in the history of instrumental observations in Svalbard. (The mean June–July–August temperature was 7.2 C, about 3 C above the climatological normal at the Svalbard Airport meteorological station.) In Hornsund the same summer months mean was 4.8 C (only 1.2 C higher than the local normal). The highest air temperature since the beginning of measurements was recorded on 25 July 2020: 21.7 and 16.5 C at Svalbard Airport and the Polish Polar Station in Hornsund respectively. Moreover, in 2019 the sea ice area on the Arctic Ocean reached the second minimum extent in the history of satellite measurements since 1979 (Yadav et al., 2020). While the summer of 2021 was colder and the minimal Arctic sea ice extent significantly larger, acceleration of the climate warming trend is proved despite interannual variations (Hanssen-Bauer et al., 2019). Such acceleration causes significant changes in the cryosphere of Svalbard and is particularly reflected in the faster melting of glaciers and thawing of the permafrost (Schuler et al., 2020; Christiansen et al., 2021). It also stimulates faster energy and mass exchange between the atmosphere, cryosphere and ocean. The above examples of transition in air temperature, sea ice extent or glacier and permafrost melting demonstrate regional differences in climate warming and subsequently response of other environmental components. Therefore, monitoring of such parameters and the environmental response to climate change is recommended to be carried out in Svalbard, where climate warming is one of the most dynamic (Nordli et al., 2014; Isaksen et al., 2016). Long-term observations allow for better quantification of observed changes and facilitate more profound understanding of these changes. This study presents the unique Arctic glaciological and meteorological data collected for 2009–2020 in southern Spitsbergen.

2 Study area

Werenskioldbreen is a well-studied, polythermal glacier located in South Spitsbergen (Fig. 1) (Baranowski, 1982; Pälli et al., 2003; Grabiec et al., 2012; Ignatiuk et al., 2015; Stachnik et al., 2016a, b; Sułowicz et al., 2020). This valley-type glacier covered an area of 27.1 km2 in 2008 (Ignatiuk et al., 2015) and 25.7 km2 in 2017 (current study) in a catchment area of 44 km2. Werenskioldbreen is divided by a medial moraine into the south-eastern part Slyngfjellbreen and the northern part Skilryggbreen accumulation area (Fig. 1). The glacier's forefield is closed by a distinct arc of the ice-cored terminal moraine with one river gorge. Such a hydrological system allows the glacier basin to be treated as a well-defined research laboratory for many hydrological and interdisciplinary studies (Majchrowska et al., 2015; Stachnik et al., 2016b; Łepkowska and Stachnik, 2018; Gwizdała et al., 2018; Stachnik et al., 2019; Osuch et al., 2022). The glacier is situated 15 km to the north from the Polish Polar Station Hornsund. The Stanisław Baranowski Spitsbergen Polar Station (University of Wrocław), a small field station, is located at the southern edge of the Werenskioldbreen terminal moraine. Both facilities greatly simplify the accessibility and logistics of research and monitoring projects.

Figure 1Location of mass-balance stakes (1–9) for 2009–2020 and the automatic weather station (AWS) on Werenskioldbreen (background: GeoEye, 10 August 2010).

3 Instruments and methodology

3.1 Meteorological monitoring

The automatic weather station (AWS) is located at an altitude of 380 m above sea level (Fig. 1), close to the average equilibrium line altitude (ELA) for the years 1959–2008 (Noël et al., 2020). The station was installed on the glacier on 15 April 2010. The AWS was mounted on a long steel mast placed in the ice drilling hole (ca. 6 m deep). In the following years, as ablation progressed, the sensors were lowered or the mast was replaced with a new one in close proximity to the original location. The recording of variables (air temperature, humidity, wind speed, short-wave and long-wave radiation) started on 17 April 2010 (Table 1). The Kipp & Zonnen CNR4 consists of two CM3 pyranometers, two CG3 pyrgeometers and temperature sensors (PT100). Pyranometers (180 solid angle) have a glass dome and measure radiation in the range from 300 to 2800 nm. One of the pyranometers directed upwards measures downwelling radiation, and the second one directed downwards measures solar radiation reflected from the earth's surface (upwelling radiations). Pyrgeometers (180–150 solid angle) have silicone windows, which allow radiation measurements in the range from 4500 to 42 000 nm. Like the pyranometers, the pyrgeometers point in opposite directions (upwards and downwards). One of the pyrgeometers measures the long-wave radiation coming from the atmosphere, the second one the long-wave radiation from the ground surface. In 2016, the A100R cup anemometer was replaced with the Gill WindSonic sensor, which allowed for measuring wind speed. Damaged sensors were replaced during the spring or autumn service to maintain data continuity.

Table 1Automatic weather station (AWS) sensors specification (Werenskioldbreen, Svalbard). All sensors were set in the spring so that they were about 1.5 m above the snow surface, and height has been varied during the season.

Download Print Version | Download XLSX

In the years 2010–2011, the sampling time was 1 min. Due to high energy demand during the polar night, the sampling time was changed to an instantaneous measurement every 10 min. Calibration and testing of sensors were performed regularly during spring expeditions based on the infrastructure of the Polish Polar Station Hornsund.

3.2 Glaciological monitoring

In 2009–2010, nine mass-balance stakes were installed on the Werenskioldbreen. The locations have been chosen to cover the elevation range from 117 to 515 m a.s.l. to create altitudinal profiles along the northern and southern tributaries of the glacier (Fig. 1).

The stakes, 6–8 m long, were embedded in the glacier by a steam drilling rig or by Kovacs ice coring system (ICS). The mass-balance stakes were measured twice a year (spring–autumn, 2009–2013) during the winter maximum accumulation (April–May) and at the end of the ablation season (September–October) or once a year (at spring, since 2014). The measure of winter accumulation was determined during the spring campaigns. The properties of snow cover (bulk snow density, snow depth and snow water equivalent (SWE)) were measured in snow pits (a 100 cm3 snow gauge by Winter Engineering was used to determine the snow density of subsequent layers) or shallow core boreholes (ICS). During the measurements, repeated soundings of the snow depth were also performed with avalanche probes. In the absence of the autumn campaign, boreholes have been drilled near each stake in order to accurately determine the amount of summer ablation and possible summer accumulation. Measurements during the autumn campaign did not always take place after or at the end of the ablation season. This was due to the logistics of the expeditions and the extension of the ablation season. In the case of availability of data from the SR50A sensor, ablation or accumulation corrections were also made if the winter or summer season ended later than the date of field observations. Some of the ablation stakes have been damaged every few years. They have been broken by wind, polar bears, melt out from the ice or been buried by snow. The network of ablation stakes was supplemented and renovated during maintenance visits. Unfortunately, recent years have resulted in large gaps in measurements due to the pandemic travel restrictions (years 2020 and 2021). Detailed information on the temporal availability of glaciological data is presented in Table 2.

Table 2Overview of mass balance and snow cover measurements on ablation stakes and infrastructure maintenance for the years 2009–2020 on Werenskioldbreen (Svalbard). S: spring campaign (winter balance, April–May); A: autumn campaign (summer balance, August–September, can be performed next year spring); X: lack of stake; SP: snow pit (SWE data); KD: ICS drilling (SWE data).

Download Print Version | Download XLSX

Based on the data collected, the following glaciological variables are available: seasonal and annual point and glacier-wide surface mass balance, snow depth, bulk snow density and SWE at the measuring points.

The analyses of the glacier's surface were based on altitude zones determined from digital elevation models (DEM). Two DEMs with geoidal height (EGM2008) were used, one generated from SPOT image acquired on 1 September 2008 (Ignatiuk et al., 2015) for the years 2008–2019 and Pleiades high-resolution images taken on 20 August 2017 (Błaszczyk et al., 2019) for year 2020.

4 Meteorological observations

4.1 Air temperature and radiation

The air temperature data forms the most homogeneous series for 2010–2016 (Fig. 2a). For 2017–2020 the data gaps were already significant due to a series of failures of the instruments. Also, for 2010–2016, net radiation balance data (short- and long-wave radiation) are available (Fig. 2b and c). For 2017–2020, radiation data included only downwelling short-wave radiation.

Figure 2Time series of meteorological variables from 2010 to 2016 on Werenskioldbreen, including daily average air temperature (a), hourly average short-wave radiation (b) and hourly average long-wave radiation (c) for downwelling (black line) and upwelling (red line) radiation.


For years with full data continuity (Fig. 2a), air temperature monthly and yearly averages were calculated and then compared with the data from the Polish Polar Station Hornsund (Wawrzyniak and Osuch, 2020). The average difference in the annual temperature (2011, 2012, 2014 and 2015) between the glacier (380 m a.s.l.) and the Polish Polar Station Hornsund (8 m a.s.l.) was −2.7C, which gives an average temperature lapse rate of 0.72 per 100 m (annual values vary from 0.55 to 0.80 per 100 m). We have calculated the significance of the trend presented in this study using the non-parametric modified Mann–Kendall test (Hamed and Rao, 1998) while considering the effect of autocorrelation of time series. The slope of the trend was calculated using Sen's method (Sen, 1968). The test indicated the statistically significant increasing temperature trend in the period 2010–2016 (with the significance level α= 0.05 and p value = 0.036), and taking into account the 12-month seasonality the Sen's slope was 0.02. Gaps in the data were filled based on the relations between air temperature measured on the PPS in Hornsund and air temperature on the WRN glacier (R2=0.96). According to the Wawrzyniak and Osuch (2020), the estimated slope of trend for air temperature between 1979 and 2018 at the Hornsund station was estimated as 1.14 C per decade. Increasing distance between the glacier surface and the sensors during the season is affecting the air temperature measurements. Periodic measurements of the vertical temperature gradient between 0.5 and 4 m carried out at the AWS indicate that the air temperature in the atmospheric boundary layer changes by about 0.2 C m−1 during the ablation season. Downwelling short-wave radiation reaches its maximum during the middle of the polar day (June). Its annual course is governed by the occurrence of polar day and night, while its daily course is governed the height of the sun above the horizon. The reflected short-wave radiation (upwelling) is a function of the surface albedo. In the spring and early summer, we observed the highest values of reflected radiation due to the presence of snow cover on the glacier. In the further part of the ablation season (July–August), we noticed a sudden decrease in reflected radiation (Fig. 2b) as a result of the disappearance of snow cover at the measurement site (i.e. the AWS) and the appearance of glacial ice on the surface. The decrease in upwelling short-wave radiation can be slight (e.g. 2012) when the melting of snow cover occurs mainly as a result of surface ablation or abrupt (e.g. 2015) when significant rainfall led to a sudden change in the albedo on the glacier. The maximum values of downwelling or upwelling long-wave radiation (Fig. 2c) usually occurred in summer and autumn. The values in winter and spring are lower, which in general shows similar patterns with the seasonal variations in air temperature. The values above 316 W m−2 of outgoing long-wave radiation may be caused by the presence of the water under the station or the presence of sediment or cryoconite. Both of these situations were observed at the AWS. In the view of the CNR4 sensor (180), there is also a mast with a logger, sensors and a solar panel, which cause distortions in the observations. These problems are difficult to eliminate. Assuming the homogeneity of the surface around the AWS, increasing the distance of the CNR4 sensor from the glacier surface should not affects its measurements.

4.2 Other variables

The AWS measured relative humidity, wind speed, ablation and accumulation of snow (for the time span see Table 1). These sensors were installed at the station depending on the needs of the ongoing projects. Not all of them could be connected to the data logger at the same time. Therefore, the data obtained for these variables are not continuous and not homogenous for the entire observation period. Nevertheless, these data are of great value for solving specific scientific problems like rain on snow events (Łupikasza et al., 2019) or supplying data to other models (Decaux et al., 2019).

5 Glaciological observations

5.1 Point ablation and accumulation

Measurements on mass-balance stakes (Fig. 1; Table 2) were performed in accordance with the recommendations and guidelines contained in the Glossary of Glacier Mass Balance and Related Terms (Cogley et al., 2011). After Cogley et al. (2011) it was assumed that accumulation is always positive, while ablation is negative. Therefore, the calculation of the point mass balance is Eq. (1):

(1) b a = c a + a a = b w + b s = c w + a w + c s + a s ,

where ba is annual balance at a point, ca is annual accumulation, aa is annual ablation, bw is winter balance, bs is summer balance, cw is winter accumulation, aw is winter ablation, cs is summer accumulation and as is summer ablation.

A method of determining point mass balance on the glacier surface includes measurements at stakes and in snow pits or boreholes. The high of each stake above snow or ice is measured twice in the maximum of winter accumulation and at the end of the summer ablation. The measurements also include depth probing and density sampling of the snow and firn (see Sect. 5.2). They are made at single points, the results from a number of points being extrapolated and integrated to yield the surface mass balance over the entire glacier (see Sect. 5.3; Cogley et al., 2011). The error of point mass balance was estimated using the total differential function using the general equation:

(2) Δ f x 1 , x 2 , x n = f x 1 2 Δ x 1 2 + f x 2 2 Δ x 2 2 + + f x n 2 Δ x n 2 ,

where Δf, Δx represent errors of the variables, f/xn is the partial derivative and x1x2,…xn are the variables.

Based on Eq. (2) we have calculated the error of winter, summer and annual point surface mass balance.

The dataset includes point winter and summer mass balance measurements on mass-balance stakes for 2009–2020 and the calculated point annual mass balance. The data allow the analysis of the spatial and seasonal variability of accumulation and ablation at points on the glaciers at different altitudes. The analysis of the winter balance (Fig. 3) shows the interannual fluctuations in snow accumulation in the entire altitude profile of the glacier. The analysis of point winter balance shows the smallest interannual fluctuations on the glacier snout (stake 9) and in the sheltered upper glacier cirque (stake 8). In the case of the point summer balance, the greatest interannual changes are observed in the middle zone of the glacier (200–400 m a.s.l.). This is due to a longer ablation period and higher temperatures not previously recorded at these altitudes. The variability of the annual surface mass balance is dominated by the summer surface mass balance (Østby et al., 2017; Grabiec et al., 2012; Van Pelt et al., 2019).

Each of the balance years can be considered separately (Fig. 4). Winter accumulation in the analysed period was generally low. The last significant accumulation on Werenskioldbreen took place in 2011. A slight accumulation in the highest parts of the glacier was also observed in 2013, 2015 and 2020. Observations from 2020, however, may be biased by unqualified substitutive observers due to the pandemic situation. In all other years, the ELA was above the highest monitored stake, no. 1 (Table 2). The data on point mass balance components are crucial for calculations of the glacier-wide surface mass balance. These data have very high importance for different modelling purposes, e.g. hydrology and glacier drainage modelling, total water discharge from the glacier, sea level rise models and validation of remote sensing products.

Figure 3Winter and summer balance at the point on 1–9 mass-balance stakes for the years 2009–2020 For the location and elevation of the stake see Fig. 1 and Table 2. Stakes have been lined up according to a height above sea level.


Figure 4Examples of winter, summer and annual point mass balance on Werenskioldbreen (season 2010–2011). Whiskers show an error (total differential function).


5.2 Snow water equivalent (SWE)

For the years 2009–2019, 49 samplings (shallow drilling or snow pits) were made on the glacier during the spring measurement campaigns in order to determine the bulk snow density, and thus SWE. In the case of snow pit measurements, the density was measured for each homogeneous layer. The bulk snow density for the snow profile was then averaged weighted by layer thickness.

Whereas the bulk snow density during the drilling of snow cores was calculated based on the length and weight of each core in the profile, SWE was calculated based on the following equation (Sturm et al., 2010):


where hs is snow depth (m), ρs is bulk snow density (kg m−3), ρw is density of water (1 kg m−3), Δhs is error of the snow depth (0.01 m) and Δρs is error of the bulk snow density (10 kg m−3).

The average density of snow cover ranged from 386 to 447 kg m−3 (Table 3). The highest snow density values were noted in 2012. They are related to the extremely warm conditions in the winter season 2011–2012 with the heavy rainfall (Łupikasza et al., 2019) during the winter and probably caused by the inflow of warm Atlantic water. (The fjords of south-west Spitsbergen did not freeze.)

Table 3Average snow depth and bulk snow density based on data from sampling points (snow pits and drilling cores) on Werenskioldbreen for the years 2009–2019.

Download Print Version | Download XLSX

Figure 5Relationship between SWE and snow depth (a) and SWE and elevation (m a.s.l.) (b). Whiskers show an error (total differential function; Eq. 4).


The SWE values show a very high correlation with the snow depth (R2=0.96; Fig. 5a) and lower correlation with the altitude above sea level (R2=0.62; Fig. 5b). Uszczyk et al. (2019) found the relationship between the bulk snow density and the snow depth on Hansbreen located next to Werenskioldbreen. It was observed that the bulk snow density increases with snow depth. The long-term data collected on Werenskioldbreen have not confirmed this correlation. In fact, in some seasons it is the opposite, i.e. thinner snow cover in the lower zones of the glacier has the highest bulk snow density. Seasonal variability can be explained by various meteorological conditions during the accumulation season. The differences between Werenskioldbreen and Hansbreen can most likely be explained by different orographic conditions and exposure, which affects snow blowing and snow deposition.

5.3 Glacier-wide surface mass balance

While the mass balance is measured on many glaciers, the data series rarely exceeds 10 years (Schuler et al., 2020). Multi-year data series, such as those from Werenskioldbreen, represent a unique value for tracking long-term changes in the Arctic environment. Calculation of the mass balance was based on point winter and summer balance analyses and digital elevation models. The point measurements are extrapolated over the glacier surface determining the balance as a function of altitude and averaging them, using the weights determined from the distribution of the glacier surface as a function of altitude (Cogley et al., 2011). The error was estimated using the total differential function.

Figure 6Annual surface mass balance and its components of Werenskioldbreen during 2009–2020. Blue bars: winter mass balance; red bars: summer mass balance; black bars: net mass balance. The results for 2019 may be understated (i.e. field measurements were performed by a non-expert crew).


The largest fluctuations are observed in the summer balance, which depends on the interannual changes in the duration of the positive air temperatures and thus the length of the ablation season. The winter balance shows greater stability; however, over the decade, the amount of snow accumulation is downward. This entails a negative trend in the surface mass balance of Werenskioldbreen (Fig. 6). Based on the trend lines, it can be concluded that acceleration of mass loss decreases by 0.09 m water equivalent (w.e.) per decade, while the summer balance decreases by −0.14 m w.e. per decade. This gives us an acceleration of the ablation by 0.23 m w.e. per decade on Werenskioldbreen.

The significance of the trends using the non-parametric modified Mann–Kendall test (Hamed and Rao, 1998) showed there is no statistically significant trend for 2009–2020 (α=0.05). The Sen's slope was −0.01, −0.12 and −0.02 for winter, summer and annual glacier-wide mass balance respectively. Hagen et al. (2012) have shown that it is impossible to give any trend for the glacier mass balance data for such a short time period.

Grabiec et al. (2012) have used monthly values of air temperature and precipitation from the meteorological station at Hornsund and the reanalysis ERA-40 data to hindcast the mass balance of Werenskioldbreen for the years 1959–2002. The average correlation coefficient of the modelled and observed mass balance in the five seasons (1994, 1999–2002) was 0.67 (meteorological data) and 0.70 (ERA-40). The average glacier-wide winter surface mass balance for 1959–2002, according to the model, was 0.81 m w.e. (ERA-40 data) and 0.87 m w.e. (meteorological data) which, compared with the years 2009-2020 (average winter glacier-wide surface mass balance 0.72 m w.e.), decreased by 7 % and 13 % respectively. The glacier-wide summer surface mass balance decreased from −1.23 m w.e. during 1959–2002 to −1.68 m w.e. during 2009–2020 (37 %) in comparison with the meteorological model, and from −1.14 m w.e. during 1959–2002 (47 %) for the ERA-40 data model. Detailed analysis of glacier-wide summer surface mass balance data from modelling (1979–2005) and observations (2009–2020) shows an increase in the average 10-year glacier-wide summer surface mass balance from −1.16 m w.e. for 1979–1988, through −1.35 and −1.55 for 1989–1998 and 1999–2005, respectively, to −1.68 for 2009–2020. A natural consequence of increasing the glacier-wide summer surface mass balance is also a much more negative average annual glacier-wide surface mass balance in the years 2009-2020 (−0.97 m w.e.) compared with the years 1959–2002 (−0.35 m w.e. for the meteorological data model and −0.34 m w.e. for ERA-40 data model).

5.4 Daily surface ablation

The influence of air temperature on the glacier surface ablation has been the subject of numerous studies. The coefficient of determination between the annual ablation and the sum of positive daily air temperature was calculated as 0.96 (R2) by Braithwaite and Olsen (1989). High correlation is caused by the strong dependence between the air temperature and the components of the energy balance (Hock, 2003). Ohmura (2001) presented the physical basis for the application of temperature ablation models, the relationship between air temperature and long-wave radiation of the atmosphere, sensible heat and incident short-wave radiation. The basic temperature ablation model is given by the equations (Braithwaite, 1995):


where T+ is sum of positive air temperatures (K) during the same period of n time steps Δt (h), DDF is the degree-day factor (mm d−1 K−1), Mm is measured ablation (m), M is melting (m w.e.) and ϱ is density (kg m−3). Melting is assumed to be zero when the air temperature is ≤0C.

On the basis of glaciological and meteorological data collected on Werenskioldbreen, daily surface ablation for May–November 2009–2020 was calculated (Fig. 7). In the case of gaps in meteorological data collected by the AWS on Werenskioldbreen, data from the Polish Polar Station located 16 km south-east were used (Wawrzyniak and Osuch, 2020). Linear regression was used to fill the gaps (R2=0.96).

Figure 7Cumulative ablation [106 m3] in May–November (2009–2020) for Werenskioldbreen.


Seasonal sums of surface ablation oscillate between about 23.7±1.7 (2019) and 64.2±4.5×106 m3 (2013), with an average of 44.7±3.1×106 m3 for 2009–2020. The value in 2019 may be underestimated due to problems with field measurements caused by the pandemic travel restrictions. The length of the ablation season determines meltwater runoff volumes. It varied in the analysed period from 134 d in 2014 to 203 d in 2016 (the average for 2009–2020 was 163 d). The surface ablation is affected by the decrease in the number of sunny days and the increase of days with precipitation and cloud cover (Wawrzyniak and Osuch, 2020). The amount of water produced by surface ablation is the largest component of the total runoff from the catchment, but precipitation can also be an important element of the water balance (Majchrowska et al., 2015).

6 Quality control and data processing

Data quality assurance includes additional measurements and calibration of equipment performed during the observation period and post-processing of the collected data. The analysis differed for the meteorological data constituting the time series and for the glaciological data.

The first stage of quality control for meteorological data consisted of visualizing each of the measurement series and reviewing the disrupted data caused by interruptions in the operation of sensors. Due to its location, the automatic measurement station (AWS) operating on Werenskioldbreen could not be maintained with high frequency. As a result, there were periodic problems with the power supply as well as with the freezing of some sensors. Power shortages manifested in the disappearance of measurements and the occurrence of isolated measurements, the correctness of which could not be confirmed, and therefore they were removed. Similarly, malfunctioning of sensors manifested in “blocking” the measurement at one value for a longer time. It mainly concerned wind speed measurements. As such values are unnatural, they were identified as erroneous and removed from the set during visual inspection. The next stage of the control was the identification of individual measurements where the values were too different compared with the previous and following measurements and that did not fit in the short-term trend. These data spikes were averaged with respect to adjacent measurements. It mainly concerned air temperature and humidity records, where such spikes are believed to be artefacts. Similarly, the analysis of the measurement series was performed in terms of unnatural values, i.e. values exceeding the permissible variability of the relative humidity or air temperature; these were a few cases. In these situations, such values were eliminated or averaged over adjacent measurements. In the last step, the same variables were compared with records from other weather stations in Svalbard. Air temperature time series have been tested with observations at the Polish Polar Station Hornsund (Wawrzyniak and Osuch, 2020). Mainly, the correlations of the variability of parameters were checked in comparison with the stations accepted as reference. Nevertheless, it should be remembered that even in the case of close points, this correlation does not have to be high or consistent due to the specificity of these stations, i.e. different shading conditions, ground, topography or exposure.

Analysis of collected data revealed some imperfections in 10 min measurements of the air humidity. There have been some measurements that slightly exceeded the value of 100 %. In this case, one of two procedures was undertaken: if the neighbouring measurements to the questionable record showed high air humidity, the exceeding value was reduced to 100 %; or if the neighbouring measurements to the questionable measurement showed low air humidity, the value was averaged from these neighbouring measurements.

A separate analysed issue is the variability or fluctuation in measurements in this 10 min series. A series of unjustified peaks in the measured values were identified, with a relatively small variability of the parameters recorded by the sensors earlier and later. The sites for this potential correction were searched for when reviewing the series on the chart. At the time of identifying such a value, the variability of consecutive measurements was analysed. Twelve measurements were analysed before the questioned measurement (±2 h). In this situation, when these fluctuations exceeded 2 standard deviations of the variability, they were averaged with the direct measurements before and after the questioned measurement. When more than 3 standard deviations of the measurements assessed according to this criterion were registered, directly after each other, they were completely removed and marked in the published files by missing values (described in the attributes of NetCDF files). This procedure was used for both air humidity and temperature. The tested measurements were compared, if possible, with their counterparts measured at the station in Hornsund, for similar dynamics of variation, which could justify similar dynamics of variability in the published measurement series. When large dynamics of the variability of the measured parameters were identified in both locations, the criterion of 3 standard deviations was used instead of the previously described 2 standard deviations. Unfortunately, the conditions of the measuring station in Hornsund are very different from the location on the glacier. For this reason, the direct possibility of comparing measurements between these locations may be limited only to the analysis of short-term trends and dynamics of the variability of both air temperature and humidity. The stations differ in height above sea level, distance from the sea and the ground on which they are located.

In the case of wind measurements, the most common problem was the one resulting from the sensor icing, which manifested in recording the same value over a longer period. The only possible correction here was to remove erroneous values throughout the occurrence.

In the case of radiation measurements, the fewest corrections were made as these sensors proved to be the most reliable. The introduced corrections concerned only sporadic jumps in the measured parameters. However, in this case, due to the large impact of cloudiness on the measured parameters, which may be marked in the measurements, only the evident cases were removed or averaged with such fluctuations. Jumps of single measurements against a background of relatively low long-term variability were identified as such cases.

Measurement series prepared and tested in such a way were used to calculate series with an hourly and daily resolution (24 h). The series was created as a result of averaging or summing up depending on the parameter under development.

Glaciological data are not collected automatically in large amounts but are based on single, unique observations that must be made with great care as they are not possible to repeat or relate to observations from other areas.

Each measurement of the ablation stake was performed twice. If the funnel melted in ice or snow around the stake, the measurement was made to the theoretical flat surface joining the edges of the funnel. In the event of a stake skewing, its total length was measured and then, if possible, the stake was replaced with a new one. Measurements of the snow depth, apart from making snow pits or shallow drilling, were always verified by taking two to three measurements with an avalanche probe. In order to obtain comparable measurements of bulk snow density (and SWE), these measurements were performed with two different methods (snow pit and shallow drilling), and a series of parallel measurements were performed showing that the difference in the calculated SWE did not exceed 5 %. In order to obtain the most accurate data from the ICS drill, the quality of the obtained ice and snow cores was checked in order to determine the precise diameter of the obtained cores.

The obtained point and glacier-wide surface mass balance calculations were compared with the data published by the World Glacier Monitoring Service (, last access: 24 May 2022) for other glaciers on Svalbard in order to verify the consistency of trends (Schuler et al., 2020). Data on surface ablation in seasons where it was possible were controlled by comparison with the data collected by the SR50 (sonic ranger) sensor, which was also used to verify the duration of the ablation season. The glaciological data were saved in the CSV files.

The quality of DEM generated from the SPOT images in 2008 was validated with the height of stakes on Werenskioldbreen (Ignatiuk et al., 2015). The median value and standard deviation of the accuracy of the DEM were −0.85 and 2.2 m respectively. Validation of the DEM generated from Pleiades images taken in 2017 was based on stake positions over neighbouring Hansbreen (Błaszczyk et al., 2019). The median value and standard deviation of DEM accuracy were −0.36 and 0.24 m respectively.

7 Dataset structure

Prepared measurement series were saved in the NetCDF (Network Common Data Form) format and placed on the server supporting OPeNDAP (, last access: 24 May 2022). The choice of this type of file is due to its universal nature. NetCDF files are in line with the modern trend of storing and publishing measurement series meeting the FAIR data principles. The collections are compliant with Unidata's Attribute Convention for Dataset Discovery (ACDD-1.3) and Climate and Forecast (CF) conventions (CF-1.8). The Attribute Convention for Dataset Discovery identify and define a list of NetCDF global attributes recommended for describing a NetCDF dataset to discovery systems such as Digital Libraries. Software tools can use these attributes for extracting metadata from datasets, and exporting to Dublin Core, DIF, ADN, FGDC, ISO 19115 etc. metadata formats. The CF metadata conventions are designed to promote the processing and sharing of files created with the NetCDF API. The conventions define metadata that provide a definitive description of what the data in each variable represents and the spatial and temporal properties of the data. This enables users of data from different sources to decide which quantities are comparable and facilitates building applications with powerful extraction, regridding and display capabilities. The CF convention includes a standard name table, which defines strings that identify physical quantities. Global Attributes of prepared NetCDF files comply with the recommendations of The Arctic Data Center (ADC) which is a service provided by the Norwegian Meteorological Institute (MET) (, last access: 24 May 2022).

All ACDD 1.3 variable attributes recommended were used. They were supplemented with the so-called _FillValue =−999.9 indicating data gaps and valid_max and valid_min describing the natural and allowed variability of these parameters in the measurement area. All measurement parameter names follow CF Standard Name Table version 77 which was available on the day when the dataset was published.

The keywords vocabulary used is consistent with the Global Change Master Directory (GCMD) keywords (, last access: 24 May 2022) developed for 20 years by The National Aeronautics and Space Administration (NASA)/gcmd-keywords) which are a hierarchical set of controlled earth science vocabularies that help ensure earth science data, services and variables are described in a consistent and comprehensive manner, and allow for the precise searching of metadata and subsequent retrieval of data, services and variables.

8 Data availability

The data are stored in two repositories that provide long-term availability, open access, DOI and license according to the FAIR principles.

Zenodo (, last access: 24 May 2022):
– meteorological data
(; Ignatiuk, 2021a),
– glaciological data
(; Ignatiuk, 2021b).
During the review of the article, successive versions of the data were corrected and updated (versions 1–4). The final version of the datasets for meteorological data is version 4 and for glaciological data is version 1.

Polish Polar Database (, last access: 24 May 2022):
– air temperature
(;jsessionid=7A0C3C8EAEA1B8F61D8F0B57177B7098#/metadata/abc6becf-97f0-4dca-b597-2fa3438f43ab; Ignatiuk and Małarzewski, 2022a),
– relative humidity
(;jsessionid=7A0C3C8EAEA1B8F61D8F0B57177B7098#/metadata/bdd6b724-d75c-49a1-83c6-eb2007107cde; Ignatiuk and Małarzewski, 2022b),
– wind speed
(;jsessionid=7A0C3C8EAEA1B8F61D8F0B57177B7098#/metadata/d0ad64ab-ad70-43d7-9383-8a9213e6c40f; Ignatiuk and Małarzewski, 2022c),
– short-wave flux
(;jsessionid=7A0C3C8EAEA1B8F61D8F0B57177B7098#/metadata/12ed9717-8cd7-4583-b2c6-089d50e6ad61; Ignatiuk and Małarzewski, 2022d,
and;jsessionid=7A0C3C8EAEA1B8F61D8F0B57177B7098#/metadata/fa3bd41b-dfbb-49e8-bdf6-7c56e9bb902f; Ignatiuk and Małarzewski, 2022e),
– long-wave flux
(;jsessionid=7A0C3C8EAEA1B8F61D8F0B57177B7098#/metadata/9309a6b1-663c-4227-9eb6-39761c1d868d; Ignatiuk and Małarzewski, 2022f,
and;jsessionid=7A0C3C8EAEA1B8F61D8F0B57177B7098#/metadata/5aa3b739-af33-4e57-bf68-7a8757985b2d; Ignatiuk and Małarzewski, 2022g).

In addition, the glacier mass balance data are stored in the World Glacier Monitoring Service database (, WGMS, 2021) and INTAROS Data Catalogue (; Ignatiuk, 2022).

All the data are also available through the Svalbard Integrated Arctic Earth Observing System (SIOS) data access portal (, last access: 24 May 2022).

9 Summary

This paper has presented details of the glaciological and meteorological dataset (2009–2020) from the Werenskioldbreen (Svalbard). The meteorological dataset includes 10 min, hourly and daily air temperature, relative humidity, short- and long-wave radiation as well as wind speed. The glaciological dataset includes point surface mass balance (winter, summer and net), snow depth, bulk density and snow water equivalent (SWE) for the mass-balance stakes, annual glacier-wide surface mass balance and modelled daily surface ablation. These data allow observations of the rapid changes taking place in the Arctic. In particular, they allow determination of the rate of climate change directly on glaciers. Werenskioldbreen mass loss is accelerating at a rate of −0.23 m w.e. per decade. These observation data already have been used to assess the hydrological models and glaciological studies. The objective of releasing these data is to improve the usage of them for calibration and validation of the remote sensing products and models, as well as to increase data reuse (Moholdt et al., 2010; Möller et al., 2011; Claremar et al., 2012; Østby et al., 2014; Błaszczyk et al., 2019).

Author contributions

Conceptualization and formal analysis were performed by DI. The original draft was written by DI, MK, and ŁM. Conceptualization, funding acquisition, and supervision were performed by JAJ. Investigation and Methodology were performed by DI, MB, TB, MG, MK, ML, and ŁS. Data curation and data preparation were provided by TB and ŁM. All authors wrote, reviewed, and edited the paper.

Competing interests

The contact author has declared that neither they nor their co-authors have any competing interests.


Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


This study presents part of the results from the project “Hindcasting and projections of hydro-climatic conditions of Southern Spitsbergen” and the “Arctic climate system study of ocean, sea ice, and glaciers interactions in Svalbard area” (AWAKE2), supported by the National Centre for Research and Development within the Polish–Norwegian Research Cooperation Programme and the SvalGlac – Sensitivity of Svalbard glaciers to climate change. Meteorological data have been processed under assessment of the University of Silesia in Katowice data repository within the Svalbard Integrated Arctic Earth Observing System (SIOS). Glaciological data have been processed under assessment of the University of Silesia in Katowice data repository within the project Integrated Arctic Observing System (INTAROS). A legacy from the ice2sea 7th FP UE was used. The studies were carried out as part of the scientific activity of the Centre for Polar Studies (University of Silesia in Katowice) with the use of research and logistic equipment (monitoring and measuring equipment, sensors, multiple AWS, GNSS receivers, snowmobiles and other supporting equipment) of the Polar Laboratory of the University of Silesia in Katowice. The authors thank the employees, doctoral students and students of the University of Silesia in Katowice for their help in the fieldwork. We would also like to thank our colleagues from the University of Wrocław for their hospitality at the Stanisław Baranowski Polar Station. Also thanks to the members of the expedition from the Polish Polar Station Hornsund for their cooperation. We would like to acknowledge Inger Jennings for the linguistic proofreading and editor Kenneth Mankoff for the comments and support.

Financial support

This research has been supported by the Polish National Science Centre (grant no. 2017/27/B/ST10/01269), the National Centre for Research and Development (AWAKE2, grant no. Pol-Nor/198675/17/2013), the European Science Foundation (SvalGlac), the EU Horizon 2020 project (INTAROS, grant no. 727890), and the Research Council of Norway (project no. 291644, Svalbard Integrated Arctic Earth Observing System – Knowledge Centre, operational phase).

Review statement

This paper was edited by Kenneth Mankoff and reviewed by two anonymous referees.


Baranowski, S.: Naled ice in front of some Spitsbergen glaciers, J. Glaciol., 28, 211–214,, 1982. 

Błaszczyk, M., Jania, J. A., and Kolondra, L.: Fluctuations of tidewater glaciers in Hornsund Fjord (Southern Svalbard) since the beginning of the 20th century, Polish Polar Res., 34, 327–352,, 2013. 

Błaszczyk, M., Ignatiuk, D., Grabiec, M., Kolondra, L., Laska, M., Decaux, L., Jania, J., Berthier, E., Luks, B., Barzycka, B., and Czapla, M.: Quality assessment and glaciological applications of digital elevation models derived from space-borne and aerial images over two tidewater glaciers of Southern Spitsbergen, Remote Sens., 11, 1121,, 2019. 

Błaszczyk, M., Jania, J. A., Ciepły, M., Grabiec, M., Ignatiuk, D., Kolondra, L., Kruss, A., Luks, B., Moskalik, M., Pastusiak, T., Strzelewicz, A., Walczowski, W., and Wawrzyniak, T.: Factors controlling terminus position of Hansbreen, a tidewater glacier in Svalbard, J. Geophys. Res.-Earth Surf., 126, e2020JF005763,, 2021. 

Box, J. E., Colgan, W. T., Wouters, B., Burgess, D. O., O'Neel, S., Thomson, L. I., and Mernild, S. H.: Global sea-level contribution from Arctic land ice: 1971 to 2017, Environ. Res. Lett., 13, 1–11,, 2019. 

Braithwaite, R. J.: Positive degree-day factors for ablation on the Greenland ice sheet studied by energy-balance modelling, J. Glaciol., 41, 153–160,, 1995. 

Braithwaite, R. J. and Olesen, O. B.: Calculation of glacier ablation from air temperature, West Greenland, in: Glacier fluctuations and climatic change, edited by: Oerlemans, J., Kluwer Academic Publishers, Dordrecht, 219–233, calculation-of-glacier-ablation-from-air-temperature-west-greenland(5823a1a2-551e-4270-8d60-b5d23422ebdb)/export.html (last access: 23 December 2021), 1989. 

Christiansen, H. H., Gilbert, G. L., Neumann, U., Demidov, N., Guglielmin, M., Isaksen, K., Osuch, M., and Boike, J.: Ground ice content, drilling methods and equipment and permafrost dynamics in Svalbard 2016–2019 (PermaSval), Zenodo,, 2021. 

Claremar, B., Obleitner, F., Reijmer, C., Pohjola, V., Waxegård, A., Karner, F., and Rutgersson, A.: Applying a mesoscale atmospheric model to Svalbard glaciers, Adv. Meteorol., 2012, 321649,, 2012. 

Cogley, J. G., Hock, R., Rasmussen, L. A., Arendt, A. A., Bauder, A., Jansson, P., Braithwaite, R. J., Kaser, G., Möller, M., Nicholson, L. and Zemp, M.: Glossary of glacier mass balance and related terms, IHP-VII Te., UNESCO-IHP, Paris, (last access: 23 December 2021), 2011. 

Decaux, L., Grabiec, M., Ignatiuk, D., and Jania, J.: Role of discrete water recharge from supraglacial drainage systems in modeling patterns of subglacial conduits in Svalbard glaciers, The Cryosphere, 13, 735–752,, 2019. 

Førland, E. J., Isaksen, K., Lutz, J., Hanssen-Bauer, I., Schuler, T. V., Dobler, A., Gjelten, H. M., and Vikhamar-Schuler, D.: Measured and modeled historical precipitation trends for Svalbard, J. Hydrometeorol., 21, 1279–1296,, 2020. 

Gabbi, J., Carenzo, M., Pellicciotti, F., Bauder, A., and Funk, M.: A comparison of empirical and physically based glacier surface melt models for long-term simulations of glacier response, J. Glaciol., 60, 1140–1154,, 2014. 

Grabiec, M., Budzik, T., and Głowacki, P.: Modeling and hindcasting of the mass balance of Werenskioldbreen (Southern Svalbard), Arct. Antarct. Alp. Res., 44, 164–179,, 2012. 

Grabiec, M., Ignatiuk, D., Jania, J. A., Moskalik, M., Głowacki, P., Błaszczyk, M., Budzik, T., and Walczowski, W.: Coast formation in an Arctic area due to glacier surge and retreat: The Hornbreen–Hambergbreen case from Spistbergen, Earth Surf. Process. Landf., 43, 387–400,, 2018. 

Gwizdała, M., Jeleńska, M., and Łęczyński, L.: The magnetic method as a tool to investigate the Werenskioldbreen environment (south-west Spitsbergen, Arctic Norway), Polar Res., 37, 1–9,, 2018. 

Hagen, J., Dunse, T., Eiken, T., Kohler, J., Moholdt, G., Nuth, C., Schuler, T., and Østby, T.: The mass balance of the Austfonna Ice Cap, Svalbard, 2004–2010, Geophys. Res. Abstr. 14, EGU2012-6085-1, 2012. 

Hamed, K. H. and Rao, R.: A modified Mann-Kendall trend test for autocorrelated data, J. Hydrol., 204, 182–196, 1998. 

Hanssen-Bauer, I., Førland, E.J., Hisdal, H., Mayer, S., Sandø, A. B., and Sorteberg, A. (Eds.): Climate in Svalbard 2100 – a knowledge base for climate adaptation, Norway, Norwegian Centre of Climate Services (NCCS) for Norwegian Environment Agency (Miljødirektoratet), 208 pp., (NCCS report 1/2019),, 2019. 

Hock, R.: Temperature index melt modelling in mountain areas, J. Hydrol., 282, 104–115,, 2003. 

Ignatiuk, D.: Meteorological data from the Werenskioldbreen (Svalbard) 2009–2020, Zenodo [data set],, 2021a. 

Ignatiuk, D.: Glaciological data (point mass balance, SWE, snow depth, bulk snow density, modelled runoff) from Werenskioldbreen (Svabard) 2009–2020, Zenodo [data set],, 2021b. Ignatiuk, D.: Glacier mass balance: Werenskioldbreen,, INTAROS Data Catalogue [data set], last access: 24 May 2022. 

Ignatiuk, D. and Małarzewski, Ł.:;jsessionid=7A0C3C8EAEA1B8F61D8F0B57177B7098#/metadata/abc6becf-97f0-4dca-b597-2fa3438f43ab, Geonetwork [data set], last access: 24 May 2022a. 

Ignatiuk, D. and Małarzewski, Ł.:;jsessionid=7A0C3C8EAEA1B8F61D8F0B57177B7098#/metadata/bdd6b724-d75c-49a1-83c6-eb2007107cde, Geonetwork [data set], last access: 24 May 2022b. 

Ignatiuk, D. and Małarzewski, Ł.:;jsessionid=7A0C3C8EAEA1B8F61D8F0B57177B7098#/metadata/d0ad64ab-ad70-43d7-9383-8a9213e6c40f, Geonetwork [data set], last access: 24 May 2022c. 

Ignatiuk, D. and Małarzewski, Ł.:;jsessionid=7A0C3C8EAEA1B8F61D8F0B57177B7098#/metadata/12ed9717-8cd7-4583-b2c6-089d50e6ad61, Geonetwork [data set], last access: 24 May 2022d. 

Ignatiuk, D. and Małarzewski, Ł.:;jsessionid=7A0C3C8EAEA1B8F61D8F0B57177B7098#/metadata/fa3bd41b-dfbb-49e8-bdf6-7c56e9bb902f, Geonetwork [data set], last access: 24 May 2022e. 

Ignatiuk, D. and Małarzewski, Ł.:;jsessionid=7A0C3C8EAEA1B8F61D8F0B57177B7098#/metadata/9309a6b1-663c-4227-9eb6-39761c1d868d, Geonetwork [data set], last access: 24 May 2022f. 

Ignatiuk, D. and Małarzewski, Ł.:;jsessionid=7A0C3C8EAEA1B8F61D8F0B57177B7098#/metadata/5aa3b739-af33-4e57-bf68-7a8757985b2d, Geonetwork [data set], last access: 24 May 2022g. 

Ignatiuk, D., Piechota, A., Ciepły, M., and Luks, B.: Changes of altitudinal zones of Werenskioldbreen and Hansbreen in period 1990–2008, Svalbard, AIP Conf. Proc., 1618, 280,, 2015. 

IPCC: The ocean and cryosphere in a changing climate, (last access: 1 April 2020), 2019. 

Isaksen, K., Nordli, Førland, E. J., Łupikasza, E., Eastwood, S., and Niedźwiedź, T.: Recent warming on Spitsbergen – Influence of atmospheric circulation and sea ice cover, J. Geophys. Res.-Atmos., 121, 11913–11931,, 2016. 

Łepkowska, E. and Stachnik, L.: which drivers control the suspended sediment flux in a High Arctic glacierized basin (Werenskioldbreen, Spitsbergen)?, Water, 10, 1408,, 2018. 

Łupikasza, E. B., Ignatiuk, D., Grabiec, M., Cielecka-Nowak, K., Laska, M., Jania, J., Luks, B., Uszczyk, A., and Budzik, T.: The role of winter rain in the glacial system on Svalbard, Water, 11, 334,, 2019. 

Majchrowska, E., Ignatiuk, D., Jania, J., Marszałek, H., and Wasik, M.: Seasonal and interannual variability in runoff from the Werenskioldbreen catchment, Spitsbergen, Polish Polar Res., 36, 197–224,, 2015. 

Moholdt, G., Nuth, C., Hagen, J. O., and Kohler, J.: Recent elevation changes of Svalbard glaciers derived from ICESat laser altimetry, Remote Sens. Environ., 114, 2756–2767,, 2010. 

Möller, M., Finkelnburg, R., Braun, M., Hock, R., Jonsell, U., Pohjola, V. A., Scherer, D., and Schneider, C.: Climatic mass balance of the ice cap Vestfonna, Svalbard: A spatially distributed assessment using ERA-Interim and MODIS data, J. Geophys. Res.-Earth Surf., 116, 3009,, 2011. 

Noël, B., Jakobs, C. L., van Pelt, W. J. J., Lhermitte, S., Wouters, B., Kohler, J., Hagen, J. O., Luks, B., Reijmer, C. H., van de Berg, W. J., and van den Broeke, M. R.: Low elevation of Svalbard glaciers drives high mass loss variability, Nat. Commun., 11, 4597,, 2020. 

Nordli, Ø., Przybylak, R., Ogilvie, A. E. J., and Isaksen, K.: Long-term temperature trends and variability on Spitsbergen: The extended Svalbard airport temperature series, 1898–2012, Polar Res., 33, 1898–2012,, 2014. 

Nuth, C., Gilbert, A., Köhler, A., McNabb, R., Schellenberger, T., Sevestre, H., Weidle, C., Girod, L., Luckman, A., and Kääb, A.: Dynamic vulnerability revealed in the collapse of an Arctic tidewater glacier, Sci. Rep.-UK, 9, 5541,, 2019. 

Ohmura, A.: Physical basis for the temperature-based melt-index method, J. Appl. Meteorol., 40, 753–761,<0753:PBFTTB>2.0.CO;2, 2001. 

Osuch, M. and Wawrzyniak, T.: Inter- and intra-annual changes in air temperature and precipitation in western Spitsbergen, Int. J. Climatol., 37, 3082–3097,, 2017. 

Osuch, M., Wawrzyniak, T., and Łepkowska, E.: Changes in the flow regime of High Arctic catchments with different stages of glaciation, SW Spitsbergen, Sci. Total Environ., 817, 152924,, 2022. 

Østby, T. I., Schuler, T. V., and Westermann, S.: Severe cloud contamination of MODIS Land Surface Temperatures over an Arctic ice cap, Svalbard, Remote Sens. Environ., 142, 95–102,, 2014. 

Østby, T. I., Schuler, T. V., Hagen, J. O., Hock, R., Kohler, J., and Reijmer, C. H.: Diagnosing the decline in climatic mass balance of glaciers in Svalbard over 1957–2014, The Cryosphere, 11, 191–215,, 2017. 

Pälli, A., Moore, J. C., Jania, J., Kolondra, L., and Głowacki, P.: The drainage pattern of Hansbreen and Werenskioldbreen, two polythermal glaciers in Svalbard, Polar Res., 22, 355–371,, 2003. 

Pellicciotti, F., Ragettli, S., Carenzo, M., and McPhee, J.: Changes of glaciers in the Andes of Chile and priorities for future work, Sci. Total Environ., 493, 1197–1210,, 2014. 

Schuler, T. V., Dunse, T., Østby, T. I., and Hagen, J. O.: Meteorological conditions on an Arctic ice cap – 8 years of automatic weather station data from Austfonna, Svalbard, Int. J. Climatol., 34, 2047–2058,, 2014. 

Schuler, T. V., Kohler, J., Elagina, N., Hagen, J. O. M., Hodson, A. J., Jania, J. A., Kääb, A. M., Luks, B., Małecki, J., Moholdt, G., Pohjola, V. A., Sobota, I., and Van Pelt, W. J. J.: Reconciling Svalbard Glacier Mass Balance, Front. Earth Sci., 8, 156,, 2020. 

Sen, P. K.: Estimates of the regression coefficient based on Kendall's tau, J. Am. Stat. Assoc., 63, 1379–1389,, 1968. 

Stachnik, L., Majchrowska, E., Yde, J. C., Nawrot, A. P., Cichała-Kamrowska, K., Ignatiuk, D., and Piechota, A.: Chemical denudation and the role of sulfide oxidation at Werenskioldbreen, Svalbard, J. Hydrol., 538, 177–193,, 2016a. 

Stachnik, Ł., Yde, J. C., Kondracka, M., Ignatiuk, D., and Grzesik, M.: Glacier naled evolution and relation to the subglacial drainage system based on water chemistry and GPR surveys (Werenskioldbreen, SW Svalbard), Ann. Glaciol., 57, 19–30,, 2016b. 

Stachnik, Ł., Yde, J. C., Nawrot, A., Uzarowicz, Ł., Łepkowska, E., and Kozak, K.: Aluminium in glacial meltwater demonstrates an association with nutrient export (Werenskiöldbreen, Svalbard), Hydrol. Process., 33, 1638–1657,, 2019. 

Sturm, M., Taras, B., Liston, G. E., Derksen, C., Jonas, T., and Lea, J.: Estimating snow water equivalent using snow depth data and climate classes, J. Hydrometeorol., 11, 1380–1394,, 2010. 

Sułowicz, S., Bondarczuk, K., Ignatiuk, D., Jania, J. A., and Piotrowska-Seget, Z.: Microbial communities from subglacial water of naled ice bodies in the forefield of Werenskioldbreen, Svalbard, Sci. Total Environ., 723, 138025,, 2020. 

Uszczyk, A., Grabiec, M., Laska, M., Kuhn, M., and Ignatiuk, D.: Importance of snow as component of surface mass balance of Arctic glacier (Hansbreen, southern Spitsbergen), Polish Polar Res., 40, 311–338,, 2019.  

van Pelt, W., Pohjola, V., Pettersson, R., Marchenko, S., Kohler, J., Luks, B., Hagen, J. O., Schuler, T. V., Dunse, T., Noël, B., and Reijmer, C.: A long-term dataset of climatic mass balance, snow conditions, and runoff in Svalbard (1957–2018), The Cryosphere, 13, 2259–2280,, 2019. 

Vikhamar-Schuler, D., Isaksen, K., Haugen, J. E., Tømmervik, H., Luks, B., Schuler, T. V., and Bjerke, J. W.: Changes in winter warming events in the Nordic Arctic Region, J. Climate, 29, 6223–6244,, 2016. 

Walczowski, W., Beszczynska-Möller, A., Wieczorek, P., Merchel, M., and Grynczel, A.: Oceanographic observations in the Nordic Sea and Fram Strait in 2016 under the IO PAN long-term monitoring program AREX, Oceanologia, 59, 187–194,, 2017. 

Wawrzyniak, T. and Osuch, M.: A 40-year High Arctic climatological dataset of the Polish Polar Station Hornsund (SW Spitsbergen, Svalbard), Earth Syst. Sci. Data, 12, 805–815,, 2020. 

Wawrzyniak, T., Osuch, M., Napiórkowski, J., and Westermann, S.: Modelling of the thermal regime of permafrost during 1990–2014 in Hornsund, Svalbard, Polish Polar Res., 37, 219–242,, 2016. 

WGMS: Fluctuations of Glaciers Database, World Glacier Monitoring Service [data set], Zurich, Switzerland,, 2021. 

Yadav, J., Kumar, A., and Mohan, R.: Dramatic decline of Arctic sea ice linked to global warming, Nat. Hazards, 103, 2617–2621,, 2020. 

Short summary
This paper presents details of the glaciological and meteorological dataset (2009–2020) from the Werenskioldbreen (Svalbard). These high-quality and long-term observational data already have been used to assess hydrological models and glaciological studies. The objective of releasing these data is to improve their usage for calibration and validation of the remote sensing products and models, as well as to increase data reuse.