Nitrate stable isotopes and major ions in snow and ice samples from four Svalbard sites

Increasing reactive nitrogen (N-r) deposition in the Arctic may adversely impact N-limited ecosystems. To investigate atmospheric transport of N-r to Svalbard, Norwegian Arctic, snow and firn samples were collected from glaciers and analysed to define spatial and temporal variations (1 10 years) in major ion concentrations and the stable isotope composition (delta N-15 and delta O-18) of nitrate (NO3-) across the archipelago. The delta N-15(NO3-) and delta O-18(NO3-) averaged -4 parts per thousand and 67 parts per thousand in seasonal snow (2010-11) and -9 parts per thousand and 74 parts per thousand in firn accumulated over the decade 2001-2011. East-west zonal gradients were observed across the archipelago for some major ions (non-sea salt sulphate and magnesium) and also for delta N-15(NO3-) and delta O-18(NO3-) in snow, which suggests a different origin for air masses arriving in different sectors of Svalbard. We propose that snowfall associated with long-distance air mass transport over the Arctic Ocean inherits relatively low delta N-15(NO3-) due to in-transport N isotope fractionation. In contrast, faster air mass transport from the north-west Atlantic or northern Europe results in snowfall with higher delta N-15(NO3-) because in-transport fractionation of N is then time-limited

Controlled experiments on the diffusion rate of stable isotopes of water in artificial firn

We studied the diffusion rate of stable isotopes of water ice in a controlled laboratory experiment, where isotopically different layers of fabricated firn were stacked. This experiment was done in order to study the diffusion rate between the isotopically different layers and to compare the measurements with forward-calculated diffusion rates of the fabricated firn stack, using three different analytical and numerical methods. The layers varied in thickness within the stack, and the studied period includes variations in firn temperature. We sampled the firn stack on five different dates along the 144 day long experiment. We obtained a quantitatively good match between the calculated and measured diffusion rates, with a root-mean-square similarity of 67-80% in the depth-averaged diffusion rates. Despite the good average match, we found some qualitative discrepancy between our calculations and our measurements. The calculated diffusion rates predict >5% faster isotopic decay of thinner layers and 13% slower decay of thicker layers than our measurements show. We speculate that grain size and grain fabric may influence the tortuosity of the firn pack, and that the current generalization of tortuosity, being a function of firn density as the only free parameter, is an oversimplification

The changing impact of snow conditions and refreezing on the mass balance of an idealized Svalbard glacier.

Glacier surface melt and runoff depend strongly on seasonal and perennial snow (firn) conditions. Not only does the presence of snow and firn directly affect melt rates by reflecting solar radiation, it may also act as a buffer against mass loss by storing melt water in refrozen or liquid form. In Svalbard, ongoing and projected amplified climate change with respect to the global mean change has severe implications for the state of snow and firn and its impact on glacier mass loss. Model experiments with a coupled surface energy balance—firn model were done to investigate the climatic mass balance and the changing role of snow and firn conditions for an idealized Svalbard glacier. A climate forcing for the past, present and future (1984–2104) is constructed, based on observational data from Svalbard Airport and a seasonally dependent projection scenario. With this forcing we mimic conditions for a typical inland Svalbard glacier. Results illustrate ongoing and future firn degradation in response to an elevational retreat of the equilibrium line altitude (ELA) of 31 m decade−1. The temperate firn zone is found to retreat and expand, while cold ice in the ablation zone warms considerably. In response to pronounced winter warming and an associated increase in winter rainfall, the current prevalence of refreezing during the melt season gradually shifts to the winter season in a future climate. Sensitivity tests reveal that in a present and future climate the density and thermodynamic structure of Svalbard glaciers are heavily influenced by refreezing. Refreezing acts as a net buffer against mass loss. However, the net mass balance change after refreezing is substantially smaller than the amount of refreezing itself, which can be ascribed to melt-enhancing effects after refreezing, which partly offset the primary mass-retaining effect of refreezing

Simulating melt, runoff and refreezing on Nordenskiöldbreen, Svalbard, using a coupled snow and energy balance model

A distributed energy balance model is coupled to a multi-layer snow model in order to study the mass balance evolution and the impact of refreezing on the mass budget of Nordenski¨oldbreen, Svalbard. The model is forced with output from the regional climate model RACMO and meteorological data from Svalbard Airport. Extensive calibration and initialisation are performed to increase the model accuracy. For the period 1989–2010, we find a mean net mass balance of −0.39mw.e. a−1. Refreezing contributes on average 0.27mw.e. a−1 to the mass budget and is most pronounced in the accumulation zone. The simulated mass balance, radiative fluxes and subsurface profiles are validated against observations and are generally in good agreement. Climate sensitivity experiments reveal a non-linear, seasonally dependent response of the mass balance, refreezing and runoff to changes in temperature and precipitation. It is shown that including seasonality in climate change, with less pronounced summer warming, reduces the sensitivity of the mass balance and equilibrium line altitude (ELA) estimates in a future climate. The amount of refreezing is shown to be rather insensitive to changes in climate

Parameterizing deep water percolation improves subsurface temperature simulations by a multilayer firn model.

Deep preferential percolation of melt water in snow and firn brings water lower along the vertical profile than a laterally homogeneous wetting front. This widely recognized process is an important source of uncertainty in simulations of subsurface temperature, density, and water content in seasonal snow and in firn packs on glaciers and ice sheets. However, observation and quantification of preferential flow is challenging and therefore it is not accounted for by most of the contemporary snow/firn models. Here we use temperature measurements in the accumulation zone of Lomonosovfonna, Svalbard, done in April 2012–2015 using multiple thermistor strings to describe the process of water percolation in snow and firn. Effects of water flow through the snow and firn profile are further explored using a coupled surface energy balance - firn model forced by the output of the regional climate model WRF. In situ air temperature, radiation, and surface height change measurements are used to constrain the surface energy and mass fluxes. To account for the effects of preferential water flow in snow and firn we test a set of depth-dependent functions allocating a certain fraction of the melt water available at the surface to each snow/firn layer. Experiments are performed for a range of characteristic percolation depths and results indicate a reduction in root mean square difference between the modeled and measured temperature by up to a factor of two compared to the results from the default water infiltration scheme. This illustrates the significance of accounting for preferential water percolation to simulate subsurface conditions. The suggested approach to parameterization of the preferential water flow requires low additional computational cost and can be implemented in layered snow/firn models applied both at local and regional scales, for distributed domains with multiple mesh points

An iterative inverse method to estimate basal topography and initialize ice flow models

We evaluate an inverse approach to reconstruct distributed bedrock topography and simultaneously initialize an ice flow model. The inverse method involves an iterative procedure in which an ice dynamical model (PISM) is run multiple times over a prescribed period, while being forced with space- and time-dependent climate input. After every iteration bed heights are adjusted using information of the remaining misfit between observed and modeled surface topography. The inverse method is first applied in synthetic experiments with a constant climate forcing to verify convergence and robustness of the approach in three dimensions. In a next step, the inverse approach is applied to Nordenski¨oldbreen, Svalbard, forced with height- and timedependent climate input since 1300 AD. An L-curve stopping criterion is used to prevent overfitting. Validation against radar data reveals a high correlation (up to R = 0.89) between modeled and observed thicknesses. Remaining uncertainties can mainly be ascribed to inaccurate model physics, in particular, uncertainty in the description of sliding. Results demonstrate the applicability of this inverse method to reconstruct the ice thickness distribution of glaciers and ice caps. In addition to reconstructing bedrock topography, the method provides a direct tool to initialize ice flow models for forecasting experiments

Computing the volume response of the Antarctic Peninsula ice sheet to warming scenarios to 2200

An ice core extracted from Holtedahlfonna ice cap, western Spitsbergen, record spanning the period 1700–2005, was analyzed for major ions. The leading empirical orthogonal function (EOF) component is correlated with an index of summer melt (log([Na+]/[Mg2+]) from 1850 and shows that almost 50% of the variance can be attributed to seasonal melting since the beginning of the industrial revolution. The Holtedahlfonna d18O value is less negative than in the more easterly Lomonosovfonna ice core, suggesting that moist air masses originate from a closer source, most likely the Greenland Sea. During the Little Ice Age the lower methanesulfonic acid (MSA) concentration and MSA non-sea-salt sulfate fraction are consistent with the Greenland Sea as the main source for biogenic ions in the ice cores. Both the melt index and the MSA fraction suggest that the early decades of the 18th century may have exhibited the coldest summers of the last 300 years in Svalbard. Ammonium concentrations rise from 1880, which may result from the warming of the Greenland Sea or from zonal differences in atmospheric pollution transport over Svalbard. During winter, neutralized aerosols are trapped within the tropospheric inversion layer, which is usually weaker over open seas than over sea ice, placing Holtedahlfonna within the inversion more frequently than Lomonosovfonna

Statistical extraction of volcanic sulphate from non-Polar ice cores

Ice cores from outside the Greenland and Antarctic ice sheets are difficult to date because of seasonal melting and multiple sources (terrestrial, marine, biogenic and anthropogenic) of sulfates deposited onto the ice. Here we present a method of volcanic sulfate extraction that relies on fitting sulfate profiles to other ion species measured along the cores in moving windows in log space. We verify the method with a well dated section of the Belukha ice core from central Eurasia. There are excellent matches to volcanoes in the preindustrial, and clear extraction of volcanic peaks in the post-1940 period when a simple method based on calcium as a proxy for terrestrial sulfate fails due to anthropogenic sulfate deposition. We then attempt to use the same statistical scheme to locate volcanic sulfate horizons within three ice cores from Svalbard and a core from Mount Everest. Volcanic sulfate i