Microscale Rearrangement of Ammonium Induced by Snow Metamorphism

Introduction

Ice cores recovered from cold glaciers are invaluable archives of past climate and atmospheric composition, covering up to 800,000 years (Members Epica Community, 2004). Among the various types of impurities trapped in the ice, major ions provide information about past changes in anthropogenic pollution, atmospheric transport, forest fires, and atmospheric temperature (Legrand and Mayewski, 1997; Kaspari et al., 2007; Kellerhals et al., 2010; Eichler et al., 2011; Thompson et al., 2013). During snowfall, atmospheric constituents such as aerosol particles and gases are scavenged and deposited by wet and/or dry deposition. After deposition in polar or high-altitude regions, snow is successively transformed into firn and glacier ice, leading to an embedding of chemical impurities including major ions. However, the ice-core concentration records of major ions are not only determined by their initial atmospheric concentrations. Post-depositional processes can strongly alter the originally deposited signal. Such processes include wind erosion (Lorius et al., 1969; Dansgaard et al., 1973), re-emission of volatile compounds from the snow surface (Lalonde et al., 2002; Röthlisberger et al., 2002; Bartels-Rausch et al., 2008), sublimation (Stichler et al., 2001), migration within the snow and firn layer (Saltzman, 1995; Wolff, 1996), and relocation during melt-water percolation (Eichler et al., 2001; Li et al., 2006; Grannas et al., 2013). With increasing global temperatures, the fate of impurities during melt events is of increasing concern (Eichler et al., 2001; Meyer and Wania, 2008; Grannas et al., 2013). Species embedded in the interior of the snow matrix are less prone to being washed away by melt or rainwater than those species located at the air-ice interface. Eichler et al. (2001) have shown that in sections of ice core records that had experienced melt processes traces of impurities with high solubility in ice can still reliably be used as archives to reconstruct past atmospheric composition.

Further interest in the location of impurities comes from their role in chemistry. Snow is a fully connected porous ice structure and readily exchanges trace gases and aerosols with the overlying air (Dominé and Shepson, 2002; Grannas et al., 2007; Bartels-Rausch et al., 2014). These fluxes are an interplay of physical exchange processes and chemistry. The chemical complexity of snow with its multiphase and multicomponent character has recently raised much interest (Abbatt et al., 2012; Bartels-Rausch et al., 2014; Kahan et al., 2014). Laboratory and field studies have shown that fluxes originate from impurities hosted in various compartments of snow. Volatile impurities forming a solid solution may portioning to the gas phase, as shown for formaldehyde (Barret et al., 2011a,b). Surface adsorption and related fluxes from and to the air-ice interface is of high relevance for halogens, whereas contributions from their solid solution are thought to be negligible (Abbatt et al., 2012). Similar, the rate of higher-order chemical reactions critical depends on the location with significant differences for reactions occurring at the air-ice interface compared to those in liquid fractions of the snow, such as brine (Bartels-Rausch et al., 2014) or aerosol deposits (Koop et al., 2000; Bartels-Rausch et al., 2014). Chemistry in solid solutions is generally negligible. That the chemistry critically depends on the origin of impurities and the chemical properties of these compartments rather than simply the total impurity concentrations was recently impressively illustrated by a field study on halogen chemistry in arctic snow (Pratt et al., 2013).

While it is clear that the fate of impurities and thus their preservation in ice cores and their role in atmospheric chemistry, critical depends on their location - in particular, the distribution between solid solution, liquid fraction, and at the air-ice interface- the role of snow dynamics on this distribution is essentially open. Dry metamorphism is the most common dynamic process of snow “aging.” Snow on Earth’s ground is almost always exposed to a temperature gradient, for example between bedrock and the overlying air. Consequently, the relatively warmer ice surface sublimates and the water vapor is deposited on the cold side of the structure. Studies of snow exposed to a temperature gradient have shown that dry metamorphism leads to complete restructuring of the snowpack, i.e., to a mass turnover of the entire ice mass of up to 60% per day at a gradient of ∼50 K/m (Pinzer and Schneebeli, 2009; Pinzer et al., 2012). Although physical changes during snow metamorphism are well studied, very little is known about concurrent rearrangement of chemical species. Hewitt et al. (1991) investigated the effect of snow metamorphism on the redistribution of Cl^–, NO₃⁻, and SO₄^-. This study suggested incorporation of Cl^– and NO₃⁻ into the ice interior and exclusion of SO₄^2- to the air-ice interface with progressing snow metamorphism. Motivation to study the incorporation of ions during metamorphism came from its influence on ionic concentrations in meltwater, where excluded ions lead to ionic pulses early in the melt with potentially adverse impacts on the environmental systems.

Of the major ions in ice cores, NH₄⁺ is an essential proxy for biogenic, forest fire, and anthropogenic agricultural emissions (Döscher et al., 1996; Fuhrer et al., 1996; Eichler et al., 2009, 2011; Kellerhals et al., 2010). NH₄⁺ is widely used to date ice cores by counting annual layers due to the strong seasonal variability (Fuhrer et al., 1993). The majority of previous studies point to good preservation of NH₄⁺ concentration records even in areas strongly affected by melting (Alps (Eichler et al., 2001), South America (Ginot et al., 2010), Tian Shan (Li et al., 2006), and Svalbard (Virkkunen et al., 2007). The good preservation has been suggested to be due to rearrangement processes during snow metamorphism, including a burial of NH₄⁺ into the ice matrix, not accessible by meltwater. To the best of our knowledge, no studies on the relocation of ammonium (NH₄⁺) during snow metamorphism are available, so far. Further interest is motivated by NH₃ being the most abundant alkaline gas in the atmosphere. Consequently, NH₃/NH₄⁺ concentrations in the air and snow strongly impact the uptake of acidic species into ice (Hoog et al., 2007). Therefore, relocation of NH₄⁺ during snow metamorphism might significantly change acid-base equilibria and buffer capacities within the snow microstructure, with direct consequences for air–ice chemical exchanges and chemical reactivity (Bartels-Rausch et al., 2014).

In this work, we monitor the redistribution of NH₄⁺ with respect to five other major ions (Ca²⁺, Cl^–, F^–, Na⁺, and SO₄^2-) during temperature gradient snow metamorphism with a series of elution experiments. Elution was chosen as method because it reveals the accessibility of different ions in complex multiphase settings and thus gives an indication of their distribution in snow. Due to the complexity of the underlying mechanisms, this study focuses on the overall change in accessibility of NH₄⁺ rather than disentangling the underlying processes in detail. To address the diversity of snow in the environment, we investigate changes to the accessibility of ions in natural and artificial snow. Great care was taken to work with snow samples with concentration levels of solutes less than 1000 μg/l to ensure direct relevance and comparability of the results to the Alpine snowpack. The concentrations are comparable to the study of Hewitt et al. (1991) but significantly reduced compared to most laboratory studies (Grannas et al., 2007; Hullar and Anastasio, 2016) and correspond to the range measured in Alpine snow (Wagenbach et al., 1988; Baltensperger et al., 1993; Schwikowski et al., 1999). One key aspect of the study was to expose the snow samples to well defined temperature settings. For comparison, snow samples were stored at isothermal and temperature gradient conditions to clarify the impact of the vastly different metamorphism conditions. Metamorphism was performed for artificial and natural snow samples in a well-controlled laboratory setting for up to 3 months. The structure of the samples was monitored using computer tomography. Additionally, we investigated structural differences and ion distribution in samples taken from different depths of a natural snowpack that had undergone natural snow metamorphism to discuss the relevance of our findings for environmental settings. This approach allowed, for the first time, investigation of the microscale redistribution of NH₄⁺ with respect to other major ions during snow metamorphism.

Materials and Methods

This study combines the preparation of doped artificial snow or sampling of natural snow, metamorphism experiments, and recording elution series of the ions Ca²⁺, Cl^–, F^–, Na⁺, NH₄⁺, and SO₄^2-. All tools and containers used during snow production, sampling, and elution were carefully pre-cleaned 5 times with ultrapure water (18 MΩcm quality, arium^® pro, Sartorius, Göttingen, Germany).

Preparation of Artificial and Natural Snow Samples

Preparation of Artificial Snow

Artificial snow was produced by shock freezing sprayed droplets of a sample solution in liquid nitrogen. This method of snow production has already been successfully applied in other studies (Bartels-Rausch et al., 2004, 2013; Matykiewiczová et al., 2007; Kerbrat et al., 2010). The sample solutions with defined concentrations of the major ions Ca²⁺, Cl^–, F^–, Na⁺, NH₄⁺, and SO₄^2- were prepared by dissolving specific amounts of salts (p.a. quality, CaCl₂, NaF, and (NH₄)₂SO₄) in ultrapure water (18 MΩcm quality, arium^® pro, Sartorius, Göttingen, Germany). In the cold room of the Paul Scherrer Institute at −20°C the ice droplets were sieved to a size fraction of 300–600 μm in diameter. About 60 g of the snow samples were filled in pre-cleaned 160 ml polypropylene (PP) containers (Faust, Schaffhausen, Switzerland). The bottom of each container was already pre-filled with 4 × 10 ml layers of frozen ultrapure water, to prevent the formation of an air gap at the bottom of the snow sample caused by the vapor flux during the aging. The containers were stored isothermally at −20°C until the beginning of the metamorphism experiment. The bulk concentrations, c_Bulk, of the different major ions in the sieved snow samples were determined in triplicate analyses using ion chromatography (IC). Concentrations varied between ∼380 and 1000 μg/l (Table 1).

Sampling of Natural Snow at WFJ

After a major snowfall event on April 4th, 2017, 20 kg of fresh snow were collected from the Weissfluhjoch site (WFJ) in the Swiss Alps above Davos (2536 m a.s.ls, 46°49′47^″ N 9°48′33^″ E). The site is a reference CryoNet station (Marty and Meister, 2012). The snow was sampled into a plastic box using a pre-cleaned Teflon shovel. The plastic box was sealed and placed in an insulated container during the transport. In a −20°C laboratory at the WSL-Institute for Snow and Avalanche Research SLF, the snow was stirred using a Teflon rod to obtain a homogeneous snow sample. Bulk concentrations, c_Bulk, of the major ions in the WFJ snow samples were determined in triplicate analyses by IC and are presented in Table 1. The natural snow samples were subsampled in 160 ml PP containers (as described in the section “Preparation of Artificial Snow”) and isothermally stored at −20°C until the beginning of the metamorphism experiment.

Controlled and Natural Metamorphism

In total, we examined three different combinations of snow type and metamorphism:

(a) artificial snow exposed to a laboratory-controlled temperature gradient.

(b) natural snow exposed to a laboratory-controlled temperature gradient.

(c) natural snow from different depths of a snow pit at WFJ, that was exposed to natural metamorphism.

The combinations and the number of samples are listed in Table 2.

Laboratory-Controlled Metamorphism Experiment

The 160 ml PP containers holding the snow samples were mounted onto a heating plate set to −4°C at the base of a metamorphism box (Figure 1A). The metamorphism box (outside dimensions: 140 cm × 70 cm × 60 cm, located at the SLF) offered space for up to 28 containers. The box was filled with a 4 cm ice layer, that is, to the same height as the ice layer in the sample containers. The remaining space around the sample holders was filled with sieved snow (about 6 cm depth to the cap of the sample holders). The metamorphism box was covered with a thin aluminum plate to inhibit sublimation of the filled-in snow. The laboratory temperature was held at −8°C. Effective temperatures at the top and bottom of the 6 cm thick snow samples were −6.5 and −4.1°C, respectively, corresponding to a temperature gradient of 40 K/m. The snow samples were stored in the metamorphism box for 0, 3, 6, 12, 30, 60, and 90 days, respectively. At the end of each storage period, the corresponding containers were removed and prepared for the elution experiments. In order to keep the conditions steady throughout the experiment, empty spaces from removed containers were filled up with snow-filled dummy containers.

In addition to the snow samples aged at the defined temperature gradient in the metamorphism box, we stored samples of artificial and natural snow at isothermal conditions (−20°C) for a 90-day period. Samples were stored inside a multilayer box out of steel plates (1 cm thick) and styrofoam (2 and 5 cm thick) according to Löwe et al. (2011) to maintain a homogeneous temperature of −20°C. In Table 2, an overview of all used containers is given. For the artificial snow samples, we used two containers for the 0-, 6-, 12-, and 30-day batches, respectively: One for the elution experiment and one for monitoring structural changes using X-ray micro-computed tomography (microCT). For the 3-day batch, we had only one container (elution), for the 60-day batch 4 containers (triplicate elution/microCT) and for the 90-day isothermal sample only 1 container (microCT) as well. For the natural snow, we used 4 containers for each 0-, 3-, 6-, 12-, 30-, 60-, and 90-day batch (triplicate elution/microCT), and two containers (elution/microCT) for the 90-day isothermal batch. As initial state “0 day” we define all states of metamorphism up to and including a duration of 24 h.

Natural Metamorphism: Samples of Natural Snow From a WFJ Snow Pit

In addition to the collected surface snow (see section “Sampling of Natural Snow at WFJ”), we sampled 1 kg of snow at 3 different depths of a snow pit at the field site, WFJ. Each sample corresponded to a snowfall in January, February, March, and April 2017 and an age of 90, 60, 30, and 0 days, respectively. This snow was filled into pre-cleaned plastic boxes and transported frozen to the SLF cold lab. In the laboratory at −20°C, the snow was stirred with a Teflon rod and four 60 g subsamples of each sample were filled into 160 ml PP containers (Table 2, triplicate elution plus microCT sample). Bulk concentrations, c_Bulk, of the different major ions in the snow pit samples are presented in Table 1. The containers were isothermally stored at −20°C until the beginning of the elution. These samples were not exposed to artificial temperature gradient metamorphism in the metamorphism box.

The age of the snow in the snow profile was determined based on the density evolution. Daily characterization of the snowpack at the WFJ field is conducted by the SLF every winter season; this includes density measurements with the Snow Micro Penetrometer SMP (Calonne et al., 2019). These measurements allowed us to monitor the stratigraphic evolution within the snowpack and to determine the age and depth of particular snow layers.

Elution Experiments

The elution experiments were performed in the cold room of the SLF at a temperature of 0°C. All snow samples were tempered at −0.2°C for 2 h prior to the elution experiments. Then the snow samples were extracted by removing the ice bottom using a ceramic knife and transferred into the elution column. The elution was performed in a 10 cm long plastic column with a diameter of 6 cm (Figure 1B). The samples were rinsed with an eluent of ultrapure water (0°C). The eluent flow was set to 140 ml/h with the help of a drop regulator. After a specific saturation time of between 30 min and 1:45 h (Table 3) the eluate started to drip out at the outflow of the column and was collected in 15 ml PP vials (Sarstedt, Nümbrecht, Germany). We collected 10 eluates for each elution experiment: the first 4 with a volume of 3 ml, the second 4 with 6 ml, and finally 2 with 12 ml volume. During the elution experiment, the surface temperature of the elution station was monitored by an infrared thermometer to stay within −0.5 ± 0.5°C. We assume that melting of the samples at this temperature and applying a 0°C eluent is negligible but cannot be excluded entirely. After the collection of the 10 eluates, the snow sample in the elution column, including the retained ultrapure water, was allowed to melt at room temperature. Two aliquots were collected in 50 ml PP vials (Sarstedt, Nümbrecht, Germany). Procedural blanks taken during the elution experiments were always close to the detection limit for all investigated ions.

Vertical Distribution of Ions After Temperature Gradient Storage

To investigate the vertical distribution of the ion concentration within one sample, we split the microCT samples of the natural snow batch (0, 3, 6, 12, 60, 90, and 90 days isothermal, Table 2) into two parts: One part was used to determine the microstructure in the microCT (see section “X-Ray Microtomography (microCT)”). The second part was used to analyze the vertical distribution of the ions. The 6 cm high snow samples were vertically cut into five ∼1.2 cm slices. The potentially contaminated surface of these slices was removed with a pre-cleaned ceramic knife and the decontaminated samples stored in 50 ml PP containers at −20°C until IC analysis.

Ion Chromatography (IC)

Samples were kept at −20°C until analysis at the Paul Scherrer Institute. Concentrations of the major cations Ca²⁺, Na⁺, NH₄⁺, and anions Cl^–, F^–, SO₄^2- in the different elution and snow samples were determined using ion chromatography (IC). A Metrohm (Herisau, Switzerland) 850 Professional IC combined with an 872 Extension Module and an 858 Professional Sample Processor autosampler was used for the analyses. Cations were separated using a Metrosep C4 column (Metrohm) and 2.8 mM HNO₃ as eluent at a flow rate of 1 ml min^–1. Anions were separated using a Metrosep A Supp 10 column (Metrohm) and were eluted stepwise using first, a 1.5 mM Na₂CO₃/0.3 mM NaHCO₃ (1:1 mixture) eluent, then an 8 mM Na₂CO₃/1.7 mM NaHCO₃ (1:1 mixture) eluent at a flow rate of 0.9 ml min^–1. Possible instrumental drifts were monitored by measuring an in-house standard after every 20th sample. Detection limits were 10, 1, 0.5, 0.5, 0.5, and 5 μg/l for Ca²⁺, Cl^–, F^–, Na⁺, NH₄⁺, and SO₄^2-, respectively.

X-Ray Microtomography (microCT)

One sample of each batch (Table 2) was scanned for structural analysis using the microCT. The snow samples were extracted from the snow sample holder by separating the snow from the ice base layer. The snow samples were cut into a 1.5 cm × 1.5 cm × 6 cm column and mounted into a microCT sample holder (diameter 2 cm). The X-ray computed tomography scanner (Scanco microCT 40, SCANCO Medical AG, Switzerland) at the SLF was operated at a temperature of −15°C and had a resolution of 10 μm. The reconstructed microCT images were filtered with a Gaussian filter (support 2 voxels, standard deviation 1 voxel) and the threshold for segmentation was applied according to Hagenmuller et al. (2014). Structural parameters of the segmented ice structure were extracted with the software tools of the microCT device (Image Processing Language, Scanco Medical).

Results

Structural Development

Temperature Gradient Metamorphism

The microCT images in Figure 2 show the structural characteristics of the snow samples at different stages in metamorphism. These images display the drastic structural changes in the snow samples during metamorphism.

The initial structure of the two types of snow, natural and artificial, differed strongly. In artificial snow, the individual spheres of homogeneous diameter are partially recognizable while the natural snow shows the typical characteristics of fresh snow with branching (see Figures 2, 3, first column). In both types, the bonds between the individual particles are still little developed, so that the shape of the original deposited crystals can be presumed. Already after 12 days of storage time at a gradient of 40 K/m, the apparent differences between the two types of snow were much smaller. A significant coarsening of the structure was observed and a fully connected porous ice structure had grown in both the artificial and the natural snow. The reduction of the specific surface area (SSA) during this re-crystallization is qualitatively visible: With increasing exposure time to the temperature gradient, the coarsening of the snow particles is further enhanced (Figure 2, 3rd to 5th column). The natural snowpack samples as third snow type showed a similar structural development. In fact, its 90-day sample had undergone the largest transformation compared to the initial state at day 0 (Figure 2, last row). The evolution of the SSA of the individual snow samples during the metamorphism experiment is shown in Figure 4. SSA of natural snow decreased by a factor of 2.2 from 27.0 to 12.2 m²/kg within 30 days. The initial decrease was less pronounced for artificial snow, showing a decrease from 15.9 to 12.2 m²/kg during the 1st month. The development of the SSA in our experiments concur with results from experimental as well as modeling studies on dry metamorphism (Flanner and Zender, 2006; Taillandier et al., 2007; Pinzer et al., 2012). The snow samples taken from the WFJ depth profile showed similar SSA decrease during the first 30 days compared to the natural snow aged in the laboratory at 40 K/m (Figure 4). Afterward, the SSA of the WFJ snow pit samples decreased below that of the natural and artificial samples exposed to a controlled temperature gradient.

The natural snow samples had a lower density (∼400 kg/m³) than the artificial snow samples (∼600 kg/m³). As the containers were sealed, bulk density remained almost constant during the course of the experiment.

Isothermal Metamorphism

The two samples isothermally stored at −20°C for 90 days behaved differently than the ones exposed to the temperature gradient: Both artificial and natural snow samples experienced much less change without a temperature gradient, as shown in Figure 4. This is consistent with previous studies (Kaempfer et al., 2005; Taillandier et al., 2007). For artificial snow, the decrease in SSA from 15.9 to 12.4 m²/kg after 90 days was ∼22%, corresponding to an aging time of ∼6–30 days at a gradient of 40 K/m. The SSA of natural snow declined from 27 to 19.7 m²/kg after 90 days of isothermal storage; this decrease of ∼27% corresponds to an aging time of only ∼3–6 days at a gradient of 40 K/m.

Redistribution of NH₄⁺ and Other Ions

Elution Experiments

Figure 5 presents the release profile of Ca²⁺ and NH₄⁺ during the elution of WFJ snow samples exposed to controlled metamorphism for 0 and 30 days. Concentrations of all six investigated major ions decreased between the 1st and 10th eluate, as exemplarily shown for NH₄⁺ and Ca²⁺ in Figure 5. The concentration in the first eluate, c_A, reflects the concentration of impurities that are best accessible by the eluent (Figures 5A,B, left gray bars). Their concentration in the eluate is given by the amount of eluted ions and the volume of eluent. In this work, concentrations were normalized to the bulk concentrations, c_Bulk, of the original snow sample to allow for direct comparison between the different ions. Since c_Bulk of NH₄⁺ is about 10 times larger than that of Ca²⁺ in the natural snow sample (see Table 1), Figure 5A reveals that the amount of NH₄⁺ in the first eluate exceeds that of Ca²⁺ by this same factor of 10 at the beginning of the metamorphism experiment (0 days). The concentration in the residual elution sample c_R reflects the concentration of impurities that are least accessible by the eluent (Figures 5A,B, right gray bars). Generally, the ions with the highest concentration ratio c_A/c_Bulk in the first eluate show the lowest concentration ratio c_R/c_Bulk in the residual rest and vice versa (Figures 5A,B, right gray bars). A mass balance could be achieved within an uncertainty of <40%. The shape of the release profiles in Figure 5, where the majority of the ion load is found in the first few eluates is typical for contaminants that are readily eluted (Grannas et al., 2013). For example, laboratory elution experiments of water-soluble organics with melt water have shown that ¾ of each species is found in the first quarter of total eluates (Meyer et al., 2009a). This shape of the elution curve and in particular the absence of increasing concentration ratio toward the later part of the elution further indicates that release from insoluble aerosol deposits or other processes acting on longer time scales, such as chromatographic effects, can be neglected (Meyer et al., 2009a,b, 2011; Plassmann et al., 2010; Mann et al., 2011; Grannas et al., 2013). Support comes from a study by Hewitt et al. (1989) showing that ice surfaces do not lead to a chromatographic separation of Cl^–, NO₃⁻, and SO₄^2- in elution experiments. In this work, we use the first eluate and the residual elution sample for further analysis.

The relative concentration in the residual elution sample c_R/c_Bulk are smaller than 1 for NH₄⁺ and Ca²⁺ in both sampling sets of 0 and 30 days (Figure 5, right gray bars). Obviously, c_R of each, which is given by the amount of ions not accessible to the eluent and the dilution of the residual snow with eluent, is smaller than the respective bulk concentration. To assess the dilution, the mass of residual water and of the snow sample was noted for each experiment. Assuming that the water filling the void space in the residual sample still contains ions (0.2 concentration ratio is a typical value in the 10th eluate, Figure 5), a dilution to a concentration ratio of ∼0.7–0.8 can be estimated. Clearly, the c_R/c_Bulk ratio of 0.3 and 0.5 for NH₄⁺ and 0.4 and 0.3 for Ca²⁺ at 0 and 30 days (Figure 5), respectively, are significantly lower. Despite the uncertainties in potential melting and refreezing of ice by the eluent (see below), this back-of-the-envelope estimate carefully suggests that the remaining snow after the elution is depleted in both ion concentration relative to their bulk concentrations. Since the accessible ion content was removed during the elution, these depleted ratios confirm qualitatively the ion’s and ice’s mass balance. As the dilution factors despite their large uncertainties are the same for all ions in a specific sample, we will focus on discussing their relative behavior in this work.

The relative increase of the NH₄⁺ c_R/c_Bulk concentration ratio after 30 days compared to that of Ca²⁺ shows that NH₄⁺ is preferable found in the least accessible fraction after metamorphism. The comparison of the relative concentrations c_A/c_Bulk in the first eluate for the samples aged for 0 and 30 days (Figure 5A vs. Figure 5B, left gray bars), show a strong trend as well: Ca²⁺ became even better accessible after 30 days. On the contrary, the NH₄⁺ concentration ratio in the first eluate did not significantly change between day 0 and 30.

Temporal Changes of Ion Concentrations in First Eluate and Residual Snow

Figure 6 shows the temporal trend of the normalized ion concentrations in the first eluate, c_A/c_Bulk (Figures 6A–C), and c_R/c_Bulk in the molten residual sample (Figures 6D–F). The 1σ uncertainty for samples, where a 3-fold determination was performed, is also shown. The partially large uncertainties result from potential sampling inhomogeneities, alterations in snow metamorphism (between different sites in the metamorphism box or in the field), variations during the elution process and analytical errors.

Generally, the NH₄⁺ concentration ratio showed a pronounced increase in the residual sample with time by about a factor of 2–3 (Figures 6D,F) and ∼8 (Figure 6E), whereas SO₄^2- and Ca²⁺ ratios showed an increase in the first eluate. Depending on the different snow types, we observed varying trends with time for Cl^–, F^–, and Na⁺, the reasons for which will be discussed later:

Controlled metamorphism, artificial snow (Figures 6A,D)

There was a continuous rise of the NH₄⁺, Cl^–, F^–, and Na⁺ concentration ratios in the residual snow during the 2 months of aging; NH₄⁺ and Cl^– showed the strongest increase. In the first eluate, concentrations ratios of all ions decreased during the first 12 days and afterward there was a strong enhancement in the SO₄^2- and Ca²⁺ ratios.

Controlled metamorphism, natural snow (Figures 6B,E)

There was a continuous increase of the F^–, NH₄⁺, Cl^–, and SO₄^2- concentration ratios in the residual snow during the 1st month; NH₄⁺ and F^– showed the strongest enrichment. Between 2 and 3 months, the concentration ratios of the four ions remained rather stable. The normalized concentrations of Na⁺ and Ca²⁺ stayed almost constant with time. In the first eluate, we observed declining concentration ratios of all major ions during the first 12 days, followed by increasing values, most pronounced for Na⁺ and Ca²⁺, during the 1st month. Concentration ratios of all ions remained relatively constant after the first 30 days.

Natural metamorphism, natural snowpit (Figures 6C,F)

NH₄⁺ concentrations in the residual sample increased with time, whereas the concentrations of Cl^–, SOSO₄^2-, Na⁺, and Ca²⁺ did not significantly change. The concentration of F^– was below the detection limit. In the first eluate, concentrations of all ions peaked after 30 days. After 90 days, SOSO₄^2- and Ca²⁺ had the highest enrichment in the first eluate compared to other ions. Note: the 4 different snow samples (0, 30, 60, and 90 days) resulted from 4 different precipitation events and the metamorphism proceeded under unknown conditions.

Temporal Changes in Ion Concentrations Under Isothermal Conditions

Ion concentration ratios in the first eluate and the residual sample from the isothermally stored reference samples showed that there was little to no redistribution of the ions without temperature gradient compared to the samples of the same age exposed to temperature gradient (Figures 6B,E, gray inset). Chemical redistribution after 90 days of isothermal storage corresponds to the change after ∼6–12 days at the temperature gradient.

Enrichment of NH₄⁺ in the Residual Sample

NH₄⁺ shows an overall very high ratio in the residual samples, indicating its low tendency to be expelled from the ice bulk (Figures 6D–F). After 30 days this ratio was significantly higher than that of any other ion except F- in the natural snow.

To qualitatively compare the enrichment of NH₄⁺ with that of the other ions in the residual sample and to correct for potential bias related to changes in the elution behavior and dilution caused by structural changes of the samples with time, we investigated the ratio of the ion concentration c_R/c_Bulk (ion) to the NH₄⁺ concentration ratio c_R/c_Bulk (NH₄⁺), as shown in Figure 7. The comparison clearly shows that in artificial snow NH₄⁺ had the highest ratio in the residual sample compared to all other ions during the whole experiment, whereas Ca²⁺ and SO₄^2- revealed always the lowest (Figure 7A). Furthermore, records of the Na⁺, F^– and Ca²⁺, SO₄^2- ratios agreed well for the whole 60 days period (Figure 7A). This is different for natural snow and the natural snowpack, where this “grouping” was not observed. Here, at the beginning of the experiment, the Ca²⁺ ratio was more enhanced in the residual snow compared to NH₄⁺ (Figures 7B,C). This changed after ∼12 days in natural snow and ∼30 days in the natural snowpack. In natural snow F^– in the residual snow sample was higher in relation to NH₄⁺ right from the begin - this was even more pronounced with increasing storage time. We discuss this finding in the section “Accessibility of the Ions”.

Changes of the ions’ distribution in the samples stored isothermally for 90 days was insignificant compared to the same age with gradient (Figure 7B).

Vertical Distribution of Ions After Temperature Gradient Storage

The comparison of the concentration ratio of the ions within the 5 parts (∼1.2 cm) of the vertically split snow samples did not show any significant differences. No vertical gradient was formed during the different stages of snow metamorphism.

Discussion

Snow metamorphism alters the ice matrix in which impurities are embedded, and in particular, continues creating new ice surfaces. We debate the processes involved in snow metamorphism in the following two sections. This forms the basis for discussing the driving processes for the relocation of chemical species, especially why the structural and chemical changes were more pronounced during temperature gradient metamorphism compared to isothermal metamorphism.

Structural Changes

Temperature Gradient Metamorphism

A temperature gradient in snow creates a water vapor pressure gradient. This gradient induces a diffusive flux from high-pressure zones (warmer side of ice particles) to low pressure zones (colder side of ice particles) (Sokratov and Maeno, 1997; Pinzer et al., 2012). The net vapor flux is not a continuous flow-through from the bottom to the top of a snowpack. Rather it takes place between adjacent crystals. Yosida et al. (1955) describes this process as “hand to hand transport.” Further process details are described in Pinzer et al. (2012).

According to Pinzer et al. (2012) numerical simulations and experiments show that the vapor flux is almost structure-independent. Therefore, a constant temperature gradient leads to a constant water vapor flux. Some earlier studies based on simplified snow geometries concluded that the vapor diffusion coefficient is enhanced (Sommerfeld et al., 1987; Colbeck, 1993), while others (Voitkovskii et al., 1988; Sokratov and Maeno, 2000) concluded that there is no enhancement. A detailed discussion of this issue is given in Pinzer et al. (2012, p. 1042). Over time, the local rearrangements completely renew and thereby transform the structure of the original ice crystals. The lifetime of a particular ice volume is defined by the complete mass turnover and allows calculating its re-crystallization rate. Pinzer et al. (2012) showed a characteristic lifetime of 2–3 days for an ice volume at a gradient of 50 K/m. Accordingly, during our storage period of 90 days at a gradient of 40 K/m, the entire ice structure was continuously sublimated and rebuilt up to 40 times. This continuous rebuild is accompanied by a decrease in SSA (Taillandier et al., 2007), as observed in all experiments (Figure 4). However, even though the SSA decrease decelerated during the artificially induced metamorphism experiments (Figure 4), the re-crystallization rate did not change due to the constant temperature gradient.

Figure 4 shows that the initial SSA decrease in the artificial snow was much less pronounced compared to natural snow. This difference is due to the absence of fine crystal structures that would induce rapid SSA reduction during metamorphism. The absence of small structures in artificial snow resulted, on the one hand, from the round shape of the snow particles due to the production method and, on the other hand, from the subsequent sieving (see section “Preparation of Artificial Snow”). After 12 days, the artificial and the natural snow showed a comparable pattern: An anisotropic structure was formed from the isotropic starting material as shown in Figure 3. Moreover, the SSA values converged to a very similar value within an uncertainty of ∼10%; in agreement to earlier work concluding that the SSA of artificial snow agrees well with that of metamorphosed natural snow (Bartels-Rausch et al., 2004).

The samples taken directly from the different depths of the WFJ snow pit had been exposed to varying temperature gradients in the field. Within the temperature profile of a natural snowpack often a diurnal cycle is visible: The lower temperatures at night cool the snow surface, leading to stronger gradients at night than during the day, while warming by the sun during the day leads to an increase in the snow surface temperature which may reverse the gradient in the snowpack. The SSA after 60 and 90 days of natural metamorphism was lower compared to that of the laboratory study. This is caused by the lower density and the higher average temperature of the WFJ site (Kaempfer et al., 2005).

Isothermal Conditions

The development of the SSA under isothermal condition showed only a small change (Figure 4, dotted lines). This can be explained by a much lower re-crystallization rate. In absence of a temperature gradient, the re-crystallization is driven by the differences in surface energy (Löwe et al., 2011). There is a natural tendency to minimize those differences by reducing the local curvature of particular ice structures (Kerbrat et al., 2008). Dendritic snow crystals are thereby progressively rounded (Flin et al., 2004). The flux due to the curvature reduction is significantly smaller than the one due to the temperature gradient and therefore, this process is much slower (Kaempfer et al., 2005).

Figure 2A (last column) clearly shows that the resulting structural changes are relatively small. The images confirm that the structure is still isotropic: After 90 days of isothermal storage, the original structure of the artificial snow was still visible. Only the bonding between the original particles increased. In natural snow, the single crystals were generally preserved and visible. However, these crystals were visibly rounded, as the edged contours had disappeared due to the reduction of the surface energy (Löwe et al., 2011).

In contrast to the temperature gradient metamorphism, changes in SSA during isothermal storage can be used as direct measure of the re-crystallization rate: A reduction in SSA of 22% for artificial snow represents roughly a mass turnover of around 80% during 90 days. The estimation is based on the conversion of the SSA into equivalent sphere diameters (D = 6/SSA) and the resulting change in volume of the spheres. This compares to a turnover of up to 4000% during temperature gradient metamorphism, corresponding to the 40 re-crystallization cycles during the same time interval. This strong transformation, in which the entire ice structure is altered, forms the basis of our hypothesis that temperature gradient metamorphism has a strong effect on the (re-)distribution of the chemical species contained in snow, whereas less chemical changes are expected under isothermal conditions.

Rearrangement of NH₄⁺ and Other Major Ions During Snow Metamorphism

Progressing snow metamorphism in the laboratory and the field led to an enhanced incorporation of NH₄⁺ into less accessible ice interior regions. This occurred for all snow types. In the following we will relate structural changes of the snow during metamorphism to the observed locality of NH₄⁺ in the snow.

Initial Chemical Characterization of the Snow

Artificial snow

The artificial snow was produced by shock freezing a stock solution of Ca²⁺, Cl^–, F^–, Na⁺, NH₄⁺, and SO₄^2- and stored at −20°C. The initial location of impurities in ice is a strong function of concentration, species, and freezing conditions (Bartels-Rausch et al., 2014). Generally, since the solubility of impurities is much lower in ice than in water, solutes are largely expelled from the forming ice when freezing from a solution. A second phase at air-ice interfaces, ice-ice interfaces of grain boundaries, or in internal micropockets forms (see section “Accessibility of the Ions”). Interestingly, the analysis of the artificial snow revealed prominent ion-pairs that persist during the entire metamorphism period. Here, Ca²⁺ and SO₄^2-, Na⁺ and F^–, NH₄⁺, and Cl^– are grouped (Figures 6A,D, 7A). It appears that we have the following salts precipitating during the shock-freezing of the artificial snow: NaF at −3°C, CaSO₄ and NH₄Cl at −15°C as predicted by salt-water binary phase diagrams (Purdon and Slater, 1946; Melnikov, 1997; Lewis et al., 2010; Thomas and Dieckmann, 2010; Stefan-Kharicha et al., 2018). These ions-pairs may stay close also during the metamorphism experiments. The location of the second phase and other reservoirs in the snow will be discussed below.

Natural snow

The natural snow used for the 90-day metamorphism experiment in the laboratory contained mainly anthropogenically derived NH₄NO₃ and (NH₄)₂SO₄, explaining the majority of the ion balance (see Table 1). The origin of the other ions is less certain: Ca²⁺ (dust, CaCO₃, CaSO₄), Na⁺ (sea salt NaCl, NaCO₃), Cl^– (NaCl, NH₄Cl), and F^– (NaF, NH₄F). In natural snow, the initial location of an ion in the snow structure is determined by scavenging processes in the atmosphere. Aerosol particles containing major ions can be incorporated in the snow crystal matrix, by acting as ice nuclei. On the other hand, direct impaction of aerosol particles on snow crystals or accumulation of supercooled water droplets onto ice crystal surfaces (riming) may lead to the enrichment of these species at the edges of the crystal (Mosimann et al., 2002; Pruppacher and Klett, 2004). NH₄⁺ is less likely to be within the ice nuclei (contrary to e.g., Ca²⁺), but rather scavenged by the latter two processes.

With progressing snow metamorphism, there is a strong relocation of the chemical compounds. Throughout our elution experiment, we monitored this change of ion location by studying their accessibility during elution with water.

Accessibility of the Ions

Here, we differentiate between accessibility and preservation based on the elution behavior. The accessible fraction of the water-soluble ions are those that are eluted from the snow sample in the first 3 ml eluate after a saturation time in the order of 1 to 2 h (Table 3). Preserved are those ions that are not washed out even after 2–3 h of complete elution.

Ions that are present at the ice surface are directly dissolved with the first eluent. Further, ions in grain boundaries are well known for their relatively fast diffusion and thus, effective transport of impurities from the interior of ice to the surface (Huthwelker et al., 2006). A back-of-the-envelope calculation based on the water self-diffusion in grain boundaries of 7 × 10^–13 m²/s (Lu et al., 2009) suggests that within 1 h of saturation time during our elution experiments (Table 3), impurities could migrate 100 μm through grain boundaries. To the best of our knowledge, the diffusion of ions from grain boundaries into interstitial water has not been investigated or compared to the release from the ice interface. Concerning gas - ice exchange processes, Bartels-Rausch et al. (2004, 2013) has shown that grain boundaries have no significant impact on the adsorption of volatile trace gasses. Extrapolation to the elution, i.e., the transfer from ice to water, is hampered by the very different time scales of the investigations and by the very distinct behavior of strong acids at the ice-air interface (Zimmermann et al., 2016; Kong et al., 2017). Given these uncertainties, we refrain from attributing the source of the accessible fraction of ions. Despite this uncertainty, the important point is that ions from either reservoir are accessible to percolating water in natural and artificial snow.

The fact that Ca²⁺, SO₄^2-, and Na⁺ showed higher relative concentrations in the first eluate compared to other ions reflects the ease of their accessibility. A large fraction of these ions is either located at the air-ice interface or the ice-ice interface (within 100 μm of the surface). This finding is in line with their low solubility in ice. In an analysis of ice core data, solubilities of Na⁺, Ca²⁺, and SO₄^2- in the order of 2–10 μg/l have been suggested (Eichler et al., 2001). Thus, the low solubility of the latter 3 ions in ice and taken that we worked with concentrations of 500–1000 μg/l in artificial snow and 20–1000 μg/l in natural snow favors their presence at interfaces. In natural snow, an accumulation of Ca²⁺ in the residual part is visible on day 0 (Figures 6, 7). Ca²⁺ seems to have been exceptionally well protected from the eluent at the beginning of the experiment. This may be an indication that Ca²⁺ acted as ice nuclei (see section “Initial Chemical Characterization of the Snow”). This finding is in good agreement with observations that mineral dust particles preferably act as ice nuclei (Szyrmer and Zawadzki, 1997).

In general, NH₄⁺, F^–, and Cl^– were found at very low concentration ratios in the first eluate, but with the highest concentration ratios in the residual sample. This general finding is consistent with their high solubility in ice. Phase diagrams of solid solutions are sparse, and we currently lack predictive capabilities of solute concentration in ice. Based on the atomic radii of N and F being comparable to that of an O atom, NH₄⁺ and F^– are generally thought to substituting water molecules, in the ice lattice (Hobbs, 1974; Pruppacher and Klett, 2010). Besides substitution, solutes may also fill interstitial space in crystals. Hobbs (1974) gives a solubility of about 2 g/l for NH₄⁺ and F^–. As detailed by Hobbs and others (Huthwelker et al., 2006; Dominé et al., 2008) such solubility measurements are delicate and the apparent solubility can be significantly enhanced in presence of grain boundaries or if the thermodynamic equilibrium is not established, or if micropockets forms. We therefore caution when referring to this value. For HCl a solid solution solubility of 3 mg/l was derived by Thibert and Dominé (1997) in thermodynamically carefully controlled experiments on single crystalline ice. Whether this solubility reflects that of Cl^– or of HCl is open to debate, as is the mechanism. However, concentrations of NH₄⁺, F^–, and Cl^– in our snow samples are well below the solubility limits reported in the literature and suggests that they generally could be quantitatively incorporated into the ice structure.

Furthermore, micropockets can host impurities in ice. They are potentially present in shock-frozen ice (Hullar and Anastasio, 2016), such as our artificial snow, as well as in ice cores (Eichler et al., 2017, 2019). In principle, micropockets may be composed of salt deposits or liquid brine. In ice cores, for example, sulfate deposits appear to be common (Eichler et al., 2017, 2019). This might explain why we find sulfate at similar concentration ratio in the residual samples of natural snow as chloride, for which a very high solubility in ice crystals was found. Detailed laboratory studies have revealed that the number of micropockets in ice frozen from the liquid critically depends on the freezing rate and conditions as well as on the concentration of salt (Hullar and Anastasio, 2016; McFall et al., 2018). In shock-frozen ice, evidence was found for the presence of solutes in micropockets with a radius below 1 μm. In our shock frozen samples, we used significantly lower concentrations of ions by about a factor of 1000, also to minimize the occurrence of micropockets. However, with a total ion load of 4 10^–6 mol, one can estimate a total liquid volume of 5 μl in the artificial snow. This would correspond to thousands of 1 μm (radius) micropockets in each shock-frozen ice sphere of the artificial snow. For this back-of-the-envelope estimate, a typical concentration of brine that is in equilibrium with ice at −5°C of 1.4 M was used.

In summary, we have discussed solid solution, interfaces, micropockets and condensation nuclei as potential location of impurities in our snow samples. Despite the inability to clearly identify the location in our samples, we argued that impurities at the air-ice interface and to some extent at the ice-ice interface would represent easily accessible compartments. Whether the impurities form patches, micropockets or are present as induvial molecules there, is beyond the scope of this work. Solid solutions and internal grain micropockets represent reservoirs where the solutes are rather inaccessible.

Impact of Structural Changes During Snow Metamorphism on the Redistribution of Ions

During the first ∼12 days of the snow metamorphism, the concentration ratio of all accessible ions in the first eluate decreased. After ∼12 days, NH₄⁺, Cl^–, and F^– concentrations were gradually enhanced in the less accessible interior of the ice matrix, whereas SO₄^2- and Ca²⁺ accumulated at the ice surface. Here we investigate how those findings relate to structural changes.

The general decrease in the concentration ratios of all ions during the first 12 days was concurrent with a reduction in the SSA (Figure 4). However, the SSA decline during the first 12 days was much more pronounced in the natural snow (drop from 27.0 to 15.5 m²/kg) than in the artificial snow (drop from 15.9 to 13.5 m²/kg; see section “Temperature Gradient Metamorphism”). Changes in the accessible fraction c_A/c_Bulk of the ions in artificial and natural snow were rather similar in the meantime (decrease by a factor of ∼1.2–2 between day 0 and 12, see Figures 6A,B). Therefore, we conclude that changes in the SSA did not directly influence the accessibility of ions during the elution. This is not surprising. The surface area does indeed play a crucial role for exchange processes such as adsorption from the gas phase (Bartels et al., 2002; Kerbrat et al., 2008). In this study, however, the amount of accessible ions is not the result of a dynamic partitioning equilibrium. Rather, the absence of the SSA impact again confirms that the elution is not limited by the interface-water flux during elution.

In the first few days at best, the grain boundaries could play a role. Riche et al. (2012) demonstrated a loss of internal grain boundaries in artificial snow of up to 40% at isothermal storage at a temperature of −5°C. In our experiments, we expect even faster changes due to the high water vapor turnaround in temperature gradient metamorphism. With a gradient, grain boundaries move with the growing ice, eventually forming ice-ice interfaces when two growing areas meet (see Figure 8). How impurities respond to this is essentially unknown. However, the loss of internal grain boundaries during the initial stage of metamorphism would lead to a less effective transport of impurities from the interior of the ice to accessible surfaces and thus explain the decreasing concentrations of the ions.

Another hypothesis, that the initial decrease of c_A could be the result of changes in the accessibility due to changes in residence time of the eluent, was not confirmed. In artificial snow, the shortest saturation time is observed for day 0 and 3 (Table 3). This can be explained by the poorly connected structure at the beginning of the experiment (see section “Structural Development”) that led to a higher permeability. Thus, we expected a lower uptake of ions from the ice surface, as well as less exchange with potential pores and consequently a lower c_A. For day 6 and 12 the saturation time increased significantly without increasing c_A. Similarly, in the case of natural snow there was no correlation between residence time and c_A; during the first 12 days, the longest saturation time was contrasted with the smallest c_A. We conclude that the elution times have no influence on the results of the elution experiments. This is an important finding with regard to a further process that could influence the elution; wet snow metamorphism. As soon as there is liquid water present in the snowpack, new processes of transformation and particle growth occur (Brun, 1989). From the lack of correlation, we conclude that possible wet snow metamorphism processes play no significant role during the elution procedure. This finding is confirmed by the elution results of the isothermal samples (Figures 6B,E). Although these samples occasionally have the longest elution times (Table 3), there is barely any change in the ion distribution.

To conclude the most likely explanation for the observed decrease in the concentration ratios of all accessible ions in the first eluate during the initial ∼12 days of the snow metamorphism is its relation to changes in the grain boundary content - and the re-location of impurities initially hosted there - driven by the continual complete re-structuring of the ice matrix during the temperature gradient metamorphism.

In the second phase of the performed elution experiment, between day 12 and 30, despite negligible changes in the SSA, there was a significant enhancement of SO₄^2-, Ca²⁺ concentration ratios compared to the other ions at the accessible ice surface in artificial and natural snow samples. At the same time, NH₄⁺, Cl^–, and F^– ratios increased in the ice interior.

Diffusion might act on these timescales to redistribute impurities particularly along grain-boundaries. Yet, diffusion is driven by concentration gradients (for example, between the interfaces and the interior). Due to the lack of knowledge about diffusion constants, we could only speculate why this should lead to a separation of NH₄⁺ and Ca²⁺. More importantly, diffusion as the main process can be excluded by comparing these results to those of the isothermal experiments, showing no such separation. The clear difference to the isothermal data strongly indicates that the observed distribution manifests themselves during the temperature gradient metamorphism.

Furthermore Hörhold et al. (2012) suggest that the densification of firn is enhanced by impurities. The preferred occurrence of Ca²⁺ at the accessible ice surface is in support of this hypothesis, as Ca²⁺ could enhance grain-boundary mobility during the densification process.

Volatile solutes can migrate with the water vapor from warm to cold areas of snow. Wagnon et al. (1999) have detected a displacement of Cl^–, F^–, and NO₃⁻ in firn layers at high sulfuric acid concentrations. Indeed, the sorption to growing ice has also been found to be enhanced significantly for strong, volatile acids (De Angelis and Legrand, 1994; Dominé et al., 1995; Domine and Rauzy, 2004; Kärcher and Basko, 2004; Ullerstam and Abbatt, 2005; Kärcher et al., 2009). To test the relevance of this process, the concentration of ions was determined vertically along the sample. These vertical concentration profiles (see section “Vertical Distribution of Ions After Temperature Gradient Storage”) give no indication of an enrichment along the vapor flux. This is in line with the initial presence of non-volatile compounds in our work (such as NH₄NO₃ and (NH₄)₂SO₄ in the natural snow (see section “Initial Chemical Characterization of the Snow”), and a low overall acidity, preventing the formation of volatile species. If therefore only pure water vapor sublimates during the metamorphism, the ions must remain at the interfaces of the remaining ice. Cragin et al. (1996) sketches a similar mechanism. However, this was based on the false assumption that large snow crystals continue to grow at the expense of small ones; although older and younger particles and ice surfaces exist in parallel, the maximum lifetime of an ice volume is limited to about 3 days (see section “Structural Change Under Temperature Gradient Metamorphism”). The noteworthy result is thus that the rearrangement continues over much longer timescales and is not related to changes in the SSA. Indeed, Pinzer et al. (2012) have shown that SSA is only a poor indicator of the dynamics in the snow during metamorphism and therefore not well suited to characterize the redistribution of embedded impurities. The change in SSA is only in the first phase of metamorphism related to vapor flux (and therefore, re-crystallization) because a quasi-steady state is reached with respect to the SSA, but not for the re-crystallization during a constant temperature gradient driven water vapor flux.

At first view, the tendency of the ions to be hosted in the interior of ice - either as solid solution or in micropockets - goes along with their solubility in ice. For species where the solubility is very likely exceeded such as Ca²⁺ and SO₄^2-, we find a larger part at the interfaces (Figure 8, green crosses). The explanation by solubility, however, does not go far enough to explain the observed relocation of ions at lower concentration than their solubility. If solubility would be the only argument, there should be no difference and no change with time between the ions NH₄⁺, F^–, Cl^–, since solubility limit in ice was not reached for any of them. We therefore propose a more complex picture describing the fate of chemicals during temperature gradient metamorphism. It is based on the distribution equilibria between solid solution, partitioning to interfaces, and into micropockets. First it is important to note that ion redistribution in the isothermally stored samples was insignificant compared to the temperature gradient samples (Figure 7B), giving strong experimental evidence for the pivotal role of the re-crystallization of the ice matrix. Therefore, exchange between the solute’s reservoirs is limited at isothermal conditions. Secondly, change in the distribution seems to continue for more than 60 days. This timescale is much longer than the lifetime of the ice matrix and the question arises why the preferential distribution of solutes is not reached earlier, when the ice matrix is completely reformed for the first time: Clearly, the evaporating ice leaves the impurities in whatever compartment behind irrespective of the specific compartment (see above), remaining at the newly formed interfaces (Figure 8). Ion-free ice grows at colder neighboring areas on interfaces that might hold ions with an ice growth rate of roughly 2 nm/s. For this estimation we assumed an optically equivalent diameter D = 0.4 mm (D = 6/SSA) and a lifetime of an imaginary ice particle of about 60 h (see section “Structural Change Under Temperature Gradient Metamorphism”). The ice growth rate has thus the same order of magnitude as the migration distance based on typical diffusivities in ice that are in the range of 7 × 10¹⁶ to 2 × 10¹⁴ m²/s for HCl and HNO₃ (Thibert and Dominé, 1997, 1998). Thus NH₄⁺ and any impurity with both a diffusivity in that range and beneficial energetics of stabilization in the ice, can – driven by the strong concentration gradient – migrate into the newly formed ice forming a solid solution there. Taken the uncertainties in diffusivity, this estimate servers as rough upper limit. For ions that have a lower tendency to enter and be stabilized in the ice a smaller fraction enters the ice and a larger part remains at their initial location. With proceeding metamorphism, this location will end as air-ice interface, i.e., the location of these ions will more and more shift toward the accessible interfaces. Our data indicate that ions with a high solubility in ice also have a high tendency to be efficiently incorporated into growing ice (Figure 8, red dots), presumably because the same physical properties are at play in both processes.

For the qualitative comparison of the NH₄⁺ incorporation into the ice interior regions in relation to that of the other 5 ions during snow metamorphism, we formed enrichment sequences based on the ratio of the normalized ion concentration to the normalized NH₄⁺ concentration in the residual snow after 90 days (Figure 7). We obtained the following enrichment sequences for the 3 different experiments:

(a) laboratory-controlled metamorphism, artificial snow: NH₄⁺ Cl^– > F^– ∼ Na⁺ > SO₄^2-∼ Ca²⁺

(b) laboratory-controlled metamorphism, WFJ snow: F^– > NH₄⁺ Cl^– ∼ SO₄^2- Ca²⁺ > Na⁺

(c) uncontrolled natural metamorphism, WFJ snow: NH₄⁺ Cl^– ∼ Na⁺ > SO₄^2-∼ Ca²⁺ (F^– was below detection limit)

After 90 days of laboratory-controlled and natural snow metamorphism, NH₄⁺ was always one of the two ions most enriched in the less accessible ice interior.

The differences in the obtained enrichment sequences and thus the incorporation of the ions into the growing ice can be explained by different factors. The absolute concentration of ions plays a major role in the distribution of impurities in ice, as has been observed for trace elements in natural snow (Avak et al., 2019) and in earlier laboratory studies (Workman and Reynolds, 1950; Bartels-Rausch et al., 2014). For example, F^– concentrations in the artificial snow are two orders of magnitude higher than that of the WFJ snow samples (see Table 1). At higher concentrations, solubility might have been exceeded, explaining the less favorable F^– incorporation in the less accessible ice interior in artificial compared to natural snow. Furthermore, chemical composition has a strong influence on the solubility of certain ions in ice. As an example, NH₄⁺ enhances the solubility of F^– and Cl^– in ice (Gross and Svec, 1997).

Environmental Implications

Our findings can directly be applied to interpret results of elution studies in the field. Climate warming induces melt processes significantly influencing the preservation of atmospheric pollutants in high-altitude snowpacks and glacier ice. Understanding the post-depositional behavior of different compounds is crucial for future reconstruction of past environmental conditions from these natural archives. A variety of studies has investigated the elution behavior of NH₄⁺ and other major ions from melting snowpacks and glacier ice. Table 4 presents a compilation of published elution sequences. This compilation contains only sites where NH₄⁺ was studied together with other major ions. Furthermore, for better clarity, the listed sequences contain only the 6 major ions investigated in our study, even though further ions may be available in some of the works.

The compilation shows that in all elution studies (except one) NH₄⁺ is together with Cl^– the least mobile ion. Our study allows us, for the first time, to explicitly relate the preservation of NH₄⁺ at various sites affected by melting around the world to the burial of NH₄⁺ into the less accessible ice interior during snow metamorphism. Furthermore, the observed higher mobility of SO₄^2- and Ca²⁺ compared to NH₄⁺ in the elution studies (Table 4) can be clearly related to a stronger gradual enrichment of these two ions at the accessible ice surface during snow metamorphism, favoring a meltwater relocation.

The obvious differences between the published elution sequences might be explained by the different metamorphism stages and temperature gradients of the investigated snowpack, as well as varying chemical composition, different sources and concentration levels, and acid input (Brimblecombe et al., 1987, 1988; Cragin et al., 1996; Gross and Svec, 1997; Wu et al., 2018). Another factor is the partial occurrence of wet snow metamorphism in some study areas compared to the assumed pure dry metamorphism investigated in our work (see section “Impact of Structural Changes During Snow Metamorphism on the Redistribution of Ions”). However, despite the dissimilarities between the different study sites, the general low mobility of NH₄⁺ due to the gradual burial in the less accessible ice interior indicates that NH₄⁺ is one of the best-suited environmental proxies even in snowpacks and ice cores affected by melting. The establishment of such meltwater-resistant proxies is particularly relevant as many high-mountain glaciers worldwide are retreating or are in danger of melting.

Conclusion

In this study, we investigated the redistribution of NH₄⁺ and five other major ions (Ca²⁺, Na⁺, Cl^–, F^–, and SO₄^2-) during dry snow metamorphism of different snow types, i.e., shock-frozen artificial snow and natural snow.

Our results show that temperature gradient snow metamorphism in the laboratory and in the field leads to the incorporation of NH₄⁺ into less accessible ice interior regions. Exposed to a gradient of 40 K/m, the complete ice matrix of every sample was continuously rebuilt within 2 to 3 days, induced by the water vapor flux. We propose that the continuous reconstruction drives the selective incorporation of ions with a high solubility in ice (NH₄⁺, F^–, and Cl^–) into interior regions and the rejection of others (including Ca²⁺, SO₄^2-) to the exterior (air-ice or ice-ice interface).

This distribution of the ions results from the specific balance of ice growth rate, diffusivity, chemical composition, and the tendency to enter the ice matrix. The absence of a relation to changes in the specific surface area confirms that this parameter is not a valid indicator of water vapor fluxes in temperature gradient metamorphism. During isothermal snow metamorphism, the chemical rearrangement was not significant. This can be explained by the low re-crystallization rate, which was two orders of magnitude lower than in the experiments with a temperature gradient of 40 K/m.

With progressing snow metamorphism, NH₄⁺ was always one of the two ions most enriched in the interior ice regions relative to the total amount. The comparability of our results for all studied snow types suggests that the redistribution of NH₄⁺ in metamorphosing snow is independent on the snow type and thus, the findings can be widely applied to different snowpacks worldwide. This work explicitly links the observed preservation of NH₄⁺ at a number of melt-affected study sites to the incorporation into the ice interior during snow metamorphism. Accordingly, NH₄⁺ is one of the best-suited environmental proxies in snowpacks, firn and ice cores affected by melting.

Author Contributions

AE, MS, and TB-R designed the study. AE supervised the experiments and data analysis. JT, SA, JE, MS, and AE designed, planned, and performed the experiments. JT, SA, and JE produced the artificial snow samples. JT, SA, JE, AE, and SB collected natural snow samples in the field. JT, SA, JE, and AE performed the elution experiments. JT prepared and analyzed the microCT samples, wrote the manuscript with input from all authors, and designed the figures. SB analyzed the chemical samples. JT and AE analyzed the data. AE, TB-R, and MS edited the manuscript. All authors approved the submitted version of the manuscript.

Funding

Funding for this study was provided by the Swiss National Science Foundation (SNSF) under Grant No. 155999.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We are grateful to Matthias Jaggi (SLF) and Mario Birrer (PSI) for their technical assistance. We would like to thank Margret Matzl for her help in evaluating the microCT data and Bettina Richter for the help with Python. We greatly appreciate the comments and suggestions of the two reviewers FD and TK, which helped to improve the clarity of the manuscript.

References