Search

The potential of capillary forces to retain water in pores is an important property of snow and firn at glaciers. Meltwater suspended in pores does not contribute to runoff and may refreeze during winter, which can affect the climatic mass balance and the subsurface density and temperature. However, measurement of firn water content is challenging and few values have been reported in the literature. Here, we use subsurface temperature and density measured at the accumulation zone of Lomonosovfonna (1200 m a.s.l.), Svalbard, to derive water content of the firn profiles after the 2014 and 2015 melt seasons. We do this by comparing measured and simulated rates of freezing front propagation. The calculated volumetric water content of firn is ~1.0–2.5 vol.% above the depth of 5 m and <0.5 vol.% below. Results derived using different thermistor strings suggest a prominent lateral variability in firn water content. Reported values are considerably lower than those commonly used in snow/firn models. This is interpreted as a result of preferential water flow in firn leaving dry volumes within wetted firn. This suggests that the implementation of irreducible water content values below 0.5 vol.% within snow/firn models should be considered at the initial phase of water infiltration., Authors in thesis list of papers: Marchenko S., van Pelt W, Pohjola V., Pettersson R., Lötstedt, P. and Reijmer C.

Projected climate warming and wettening will have a major impact on the state of glaciers and seasonal snow in High Arctic regions. Following up on a historical simulation (1957–2018) for Svalbard, we make future projections of glacier climatic mass balance (CMB), snow conditions on glaciers and land, and runoff, under Representative Concentration Pathways (RCP) 4.5 and 8.5 emission scenarios for 2019–60. We find that the average CMB for Svalbard glaciers, which was weakly positive during 1957–2018, becomes negative at an accelerating rate during 2019–60 for both RCP scenarios. Modelled mass loss is most pronounced in southern Svalbard, where the equilibrium line altitude is predicted to rise well above the hypsometry peak, leading to the first occurrences of zero accumulation-area ratio already by the 2030s. In parallel with firn line retreat, the total pore volume in snow and firn drops by as much as 70–80% in 2060, compared to 2018. Total refreezing remains largely unchanged, despite a marked change in the seasonal pattern towards increased refreezing in winter. Finally, we find pronounced shortening of the snow season, while combined runoff from glaciers and land more than doubles from 1957–2018 to 2019–60, for both scenarios.

Since the first estimates of Svalbard-wide glacier mass balance were made in the early 2000s, there has been great progress in remote sensing and modeling of mass balance, existing field records have been extended, field records at new locations have been added, and there has been considerable environmental change. There is a wide spread in the available estimates of both total mass balance and surface or climatic mass balance, but there is overall agreement that the glaciers on Svalbard have been losing mass since the 1960s, with a tendency toward more negative mass balance since 2000. We define criteria to select data that are representative and of high credibility; this subset shows a more coherent evolution and reduced spread. In addition, we combine individual field mass balance records collected by different groups into a single dataset that samples glaciers across Svalbard and a range of different size classes. We find a close relationship between measured specific surface mass balance and size of the glacier, in such a way that smaller glaciers experience more negative surface mass balances. A qualitatively similar relationship between the accumulation area ratio and glacier area is found for all glaciers in the Svalbard, suggesting that the relation derived from glaciological records is not only an artifact caused by the limited number of samples (n = 12). We apply this relation to upscale measured surface mass balance for a new estimate for all glaciers of Svalbard. Our reconciled estimates are −7 ± 4 Gt a–1 (2000–2019) for the climatic mass balance, and −8 ± 6 Gt a–1 for the total mass balance. The difference between the two represents the sum of frontal ablation and the combined uncertainty, which together amount to ca. −2 ± 7 Gt a–1. While this is consistent with a previous estimate of Svalbard-wide frontal ablation, the uncertainties are large. Furthermore, several large and long-lasting surges have had considerable and multi-year impact on the total m

In this paper, we study the effects of basal friction, sub-aqueous undercutting and glacier geometry on the calving process by combining six different models in an offline-coupled workflow: a continuum-mechanical ice flow model (Elmer/Ice), a climatic mass balance model, a simple sub-glacial hydrology model, a plume model, an undercutting model and a discrete particle model to investigate fracture dynamics (Helsinki Discrete Element Model, HiDEM). We demonstrate the feasibility of reproducing the observed calving retreat at the front of Kronebreen, a tidewater glacier in Svalbard, during a melt season by using the output from the first five models as input to HiDEM. Basal sliding and glacier motion are addressed using Elmer/Ice, while calving is modelled by HiDEM. A hydrology model calculates subglacial drainage paths and indicates two main outlets with different discharges. Depending on the discharge, the plume model computes frontal melt rates, which are iteratively projected to the actual front of the glacier at subglacial discharge locations. This produces undercutting of different sizes, as melt is concentrated close to the surface for high discharge and is more diffuse for low discharge. By testing different configurations, we show that undercutting plays a key role in glacier retreat and is necessary to reproduce observed retreat in the vicinity of the discharge locations during the melting season. Calving rates are also influenced by basal friction, through its effects on near-terminus strain rates and ice velocity.

Knowledge of the ice thickness distribution of glaciers and ice caps is an important prerequisite for many glaciological and hydrological investigations. A wealth of approaches has recently been presented for inferring ice thickness from characteristics of the surface. With the Ice Thickness Models Intercomparison eXperiment (ITMIX) we performed the first coordinated assessment quantifying individual model performance. A set of 17 different models showed that individual ice thickness estimates can differ considerably – locally by a spread comparable to the observed thickness. Averaging the results of multiple models, however, significantly improved the results: on average over the 21 considered test cases, comparison against direct ice thickness measurements revealed deviations on the order of 10 ± 24 % of the mean ice thickness (1σ estimate). Models relying on multiple data sets – such as surface ice velocity fields, surface mass balance, or rates of ice thickness change – showed high sensitivity to input data quality. Together with the requirement of being able to handle large regions in an automated fashion, the capacity of better accounting for uncertainties in the input data will be a key for an improved next generation of ice thickness estimation approaches.

Spatial heterogeneity of snow and firn properties on glaciers introduces uncertainty in interpretation of point and profile observations and complicates modelling of meltwater percolation and runoff. Here we present a study of the temporal and spatial dynamics of firn density and stratigraphy at the plot-scale (approximate to 10 m x 10 m x 10 m) repeated annually during 2012-14 at the Lomonosovfonna ice-field, Svalbard. Results from cores, video inspections in boreholes and radar grid surveys are compared. Ice layers 0.1-50 cm thick comprised approximate to 8% of the borehole length. Most of them are 1-3 cm thick and could not be traced between boreholes separated by 3 m. Large lateral variability of firn structure affects representativeness of observations in single holes and calls for repeated studies in multiple points to derive a representative stratigraphy signal. Radar reflections are poorly correlated with ice layers in individual boreholes. However, the match between the high amplitude peaks in the grid-averaged radar signal and horizons of preferential ice layer formation revealed by averaging the video surveys over multiple boreholes is higher. These horizons are interpreted as buried firn layers previously exposed to melt-freeze or wind-driven densification and several of them are consistently recovered throughout three field campaigns.

Svalbard climate is undergoing amplified change with respect to the global mean. Changing climate conditions directly affect the evolution of the seasonal snowpack, through its impact on accumulation, melt, and moisture exchange. We analyze long-term trends and spatial patterns of seasonal snow conditions in Svalbard between 1961 and 2012. Downscaled regional climate model output is used to drive a snow modeling system (SnowModel), with coupled modules simulating the surface energy balance and snowpack evolution. The precipitation forcing is calibrated and validated against snow depth data on a set of glaciers around Svalbard. Climate trends reveal seasonally inhomogeneous warming and a weakly positive precipitation trend, with strongest changes in the north. In response to autumn warming the date of snow onset increased (2days decade(-1)), whereas in spring/summer opposing effects cause a nonsignificant trend in the snow disappearance date. Maximum snow water equivalent (SWE) in winter/spring shows a modest increase (+0.01 meters water equivalent (mwe)decade(-1)), while the end-of-summer minimum snow area fraction declined strongly (from 48% to 36%). The equilibrium line altitude is highest in relatively dry inland regions, and time series show a clear positive trend (25mdecade(-1)) as a result of summer warming. Finally, rain-on-snow in the core winter season, affecting ground ice formation and limiting access of grazing animals to food supplies, peaks during specific years (1994, 1996, 2000, and 2012) and is found to be concentrated in the lower lying coastal regions in southwestern Svalbard.

Physical and chemical properties of four different ice cores (LF-97, LF-08, LF-09 and LF-11) drilled at Lomonosovfonna, Svalbard, were compared to investigate the effects of meltwater percolation on the chemical and physical stratigraphy of these records. A synthetic ice core approach was employed as reference record to estimate the ionic relocation and meltwater percolation length at this site during the period 2007-2010. Using this method, a partial ion elution sequence obtained for Lomonosovfonna was NO(3)(-)aEuro-> aEuro-SO42-, Mg2+, Cl-, K+, Na+ with nitrate being the most mobile within the snowpack. The relocation length of most of the ions was on the order of 1aEuro-m during this period. In addition, by using both a positive degree day (PDD) and a snow-energy model approaches to estimate the percentage of melt at Lomonosovfonna, we have calculated a melt percentage (MP) of the total annual accumulation within the range between 48 and 70aEuro-%, for the period between 2007 and 2010, which is above the MP range suggested by the ion relocation evidenced in the LF-syn core (i.e., MPaEuro-aEuro parts per thousand= aEuro-30aEuro-%). Using a firn-densification model to constrain the melt range, a MP of 30aEuro-% was found over the same period, which is consistent with the results of the synthetic ice core approach, and a 45aEuro-% of melt for the last 60 years. Considering the ionic relocation lengths and annual melt percentages, we estimate that the atmospheric ionic signal remains preserved in recently drilled Lomonosovfonna ice cores at an annual or bi-annual resolution when weather conditions were similar to those during the 2007-2010 period.