Measurement report: Spatial variations in ionic chemistry and water-stable isotopes in the snowpack on glaciers across Svalbard during the 2015–2016 snow accumulation season

The Svalbard archipelago, located at the Arctic sea-ice edge between 74 and 81 N, is ∼ 60 % covered by glaciers. The region experiences rapid variations in atmospheric flow during the snow season (from late September to May) and can be affected by air advected from both lower and higher latitudes, which likely impact the chemical composition of snowfall. While long-term changes in Svalbard snow chemistry have been documented in ice cores drilled from two high-elevation glaciers, the spatial variability of the snowpack composition across Svalbard is comparatively poorly understood. Here, we report the results of the most comprehensive seasonal snow chemistry survey to date, carried out in April 2016 across 22 sites on seven glaciers across the archipelago. At each glacier, three snowpits were sampled along the altitudinal profiles and the collected samples were analysed for major ions (Ca2+, K, Na, Mg2+, NH+4 , SO2− 4 , Br , Cl, and NO−3 ) and stable water isotopes (δ 18O, δ2H). The main aims were to investigate the natural and anthropogenic processes influencing the snowpack and to better understand the influence of atmospheric aerosol transport and deposition patterns on the snow chemical composition. The snow deposited in the southern region of Svalbard is characterized by the highest total ionic loads, mainly attributed to sea-salt particles. Both NO−3 and NH + 4 in the seasonal snowpack reflect secondary aerosol formation and post-depositional changes, resulting in very different spatial deposition patterns: NO−3 has its highest loading in northwestern Spitsbergen and NH+4 in the south-west. The Br − Published by Copernicus Publications on behalf of the European Geosciences Union. 3164 E. Barbaro et al.: Ionic chemistry in the snowpack across Svalbard enrichment in snow is highest in north-eastern glacier sites closest to areas of extensive sea-ice coverage. Spatial correlation patterns between Na and δ18O suggest that the influence of long-range transport of aerosols on snow chemistry is proportionally greater above 600–700 m a.s.l.


Introduction
Svalbard is a region of the Arctic experiencing rapid climate change. The mean warming rate is +1.35 K per decade, much faster than the global average (Isaksen et al., 2016;Maturilli et al., 2013;Nordli et al., 2014). This archipelago is located at the southern edge of the perennial Arctic sea ice in the North Atlantic Ocean and is characterized by a maritime climate with large, rapid temperature variations during winter (Brage et al., 2014). South-westerly inflow of mild oceanic air, associated with a low-pressure system east of Iceland, often brings relatively warm and moist air in the winter months, while Arctic air intrusions from the northeast, driven by a high-pressure system over Greenland, result in much colder temperatures (Rinke et al., 2017). In addition to these synoptic fluctuations, intense autumn or winter cyclonic storms associated with anomalous warming events sometimes occur, transporting both heat and moisture from lower latitudes to Svalbard (Rinke et al., 2017).
The aforementioned meteorological conditions also favour long-range transport of aerosols to the archipelago, including pollutants from continental sources. Depending on the predominant air flow pattern at the time of snowfall, the archipelago may experience regionally different amounts of both snow accumulation (Eneroth et al., 2003;Forland et al., 2011) and chemical loads, the latter reflecting contrasting mixtures of aerosols, varying by source area (Aas et al., 2016;Forsström et al., 2009;Möller and Kohler, 2018). These regional differences are also associated with contrasts in seaice cover. While all Svalbard coasts are usually ice-free in summer, sea ice can form and cover large parts of the ocean surface in the eastern and northern parts of the archipelago. Contrastingly, the southern and western parts often remain ice-free (Dahlke et al., 2020) and therefore tend to experience greater snowfall owing to the proximity of open water. In addition, the West Spitsbergen Current, a branch of the Atlantic Meridional Overturning Circulation (AMOC) that flows to the west of the archipelago, causes markedly different regional climatic conditions between its eastern and western parts (van Pelt et al., 2019): the west exhibits higher temperatures and precipitation, while the east is less humid and cooler and has also experienced a stronger warming trend since 1957.
The seasonal snowpack contains a complex mixture of impurities that are either scavenged from the atmosphere during snowfall or directly received through dry deposition (Kuhn, 2001). On land, the majority of impurities found in seasonal snow are usually eluted during summer melt, influencing terrestrial and aquatic systems Brimblecombe et al., 1987). However, in the accumulation area of Arctic glaciers and ice caps, impurities can be retained within or below the seasonal snow layer Pohjola et al., 2002;Vega et al., 2016). For this reason, chemical impurities such as major ions (Ca 2+ , K + , Na + , Mg 2+ , NH + 4 , SO 2− 4 , Br − , Cl − , and NO − 3 ) in ice cores have been widely used to study the past trends of atmospheric and climatic conditions Isaksson et al., 2003;Thompson et al., 2002;Wolff et al., 2010). Previous studies in Svalbard (Isaksson et al., 2001;Matoba et al., 2002;Nawrot et al., 2016;Semb et al., 1984) have shown that the chemistry of the seasonal snowpack is dominated by sea-salt ions (Hodgkins and Tranter, 1998). However, the region is also a sink for atmospheric contaminants brought in by long-range transport (Vecchiato et al., 2018). Investigations of precipitation and snow cover chemistry have predominantly focused on the central and western parts of the archipelago (Kühnel et al., 2011;Nawrot et al., 2016;Vega et al., 2015;Virkkunen et al., 2007) due to the accessibility of research facilities in these sectors.
In general, stable water isotope measurements in different components of the water cycle are available in isotope databases, maintained and updated by the International Atomic Energy Agency (IAEA) but also by national or international organizations (West et al., 2010). Moreover, in Svalbard, stable water isotope investigations are performed in ice cores and surface snowpit samples because δ 18 O and δ 2 H are still the most common tools for finding the depthtime relationship in ice cores Punning et al., 1986). The preservation of uninterrupted annual isotope cycles varies depending on the site: in sites such as central Greenland annual isotope cycles are well preserved, while in sites with high intra-seasonal variations or with different pre-and post-depositional processes the annual layers can be difficult to distinguish (Pohjola et al., 2002). To investigate to what degree δ 18 O in snow changes after accumulation, Igarashi et al. (2001) carried out the observation of the precipitation at Ny-Ålesund. These authors concluded that the fluctuation of δ 18 O could not be explained by changes in surface air temperature only but that the characteristics of the air masses also influenced the isotope signature of the precipitation (Igarashi et al., 2001). At the moment, there is a lack of data regarding the stable water isotope composition of surface snow from Svalbard, and this survey is a substantial contribution to fill that gap.
In the present study, the concentration, mass loading, spatial and altitudinal distribution of major ion species (Ca 2+ , K + , Na + , Mg 2+ , NH + 4 , SO 2− 4 , Br − , Cl − , and NO − 3 ) in snow, together with its stable oxygen and hydrogen isotope composition (δ 18 O and δ 2 H), were evaluated in the late winter snowpack at 22 glacier sites across Svalbard. Stable isotope ratios were used as supporting data to define the accumula- tion seasonality in the snowpack and to identify the moisture sources that feed snowfall, thereby providing clues as to the predominant air transport pathways to the snowpit sites (Gat et al., 2001). This study was part of the larger Community Coordinated Snow Study in Svalbard (C2S3) project and is the most comprehensive survey of seasonal snow chemistry in Svalbard to date. The snowpack survey, which was carried out by coordinated teams using a common sampling protocol (Gallet et al., 2018), aimed to map and characterize regional differences in the chemical composition and impurity load of the winter snowpack. We further interpret the observed differences in chemical loading in relation to meteorological and other environmental factors.
Thereby, we aim to identify the conditions controlling the chemistry of Svalbard snow that are susceptible to the variable climate warming impact across the region.

Sampling location and strategy
During April 2016, the seasonal snowpack was sampled at 22 sites on seven glaciers across Svalbard (Table 1 and Fig. 1). The glaciers are of different sizes and hypsometry. Wind fields for each glacier are not available. Indeed, the wind direction can change in concomitance with cyclonic events that could occur during the season and can act differently for each glacier.
The first glacier considered in this study is Austfonna (AF), located on Nordaustlandet, the second largest island in Svalbard, with approximately 80 % of its area covered by ice. AF is the largest ice cap in Svalbard with a geographic area of 8357 km 2 and has one main central dome of up to 600 m ice thickness feeding several drainage basins (Dallmann et al., 2015;Schuler et al., 2020).
The other six glaciers investigated here are located on Spitsbergen. In north-western Spitsbergen, we studied three glaciers near Ny-Ålesund: Austre Lovénbreen (ALB), Kongsvegen (KVG), and Holtedahlfonna (HDF). ALB is a small land-based valley glacier, 4 km long from south to north along the Brøgger Peninsula. The glacier area was 4.48 km 2 in 2013 and its elevation ranges from 50 to 550 m a.s.l. The total catchment area spreads over 10 577 km 2 , taking into account an outlet where the main stream crosses a compact calcareous outcrop 400 m upstream from the coastline (Marlin et al., 2017). KVG is a north-westflowing grounded glacier located about 20 km east of Ny-Ålesund (Melvold and Hagen, 1998), with an average ice thickness of 190 m and a maximum of 450 m (Lindbäck et al., 2018). It has a total length of ca. 24 km with an average 3.5 km width, has a maximum elevation of 800 m a.s.l., and flows from south-east to north-west . HDF glacier is the largest ice field (ca. 300 km 2 ) on north-western Spitsbergen, about 40 km from the Ny-Ålesund station. It is distributed over an elevation range of 0-1241 m a.s.l. (Beaudon et al., 2017;Nuth et al., 2017). Lomonosovfonna (LF) is one of the highest ice fields on Spitsbergen and is located in the central part of the island. The summit lies at 1250 m a.s.l. and has a pronounced cupola shape with an approximate radius of 500 m. The total accumulation area of the entire LF ice system was about 600 km 2 at the beginning of the 21st century (Isaksson et al., 2001). Even though this is the highest point in our survey, the air temperature can pass above zero during the summer, resulting, although not significantly, in the relocation of ions (Pohjola et al., 2002;Vega et al., 2016). In southern Spitsbergen, two different glaciers were investigated, Hansbreen (HB) and Werenskiöldbreen (WB), close to the Hornsund station. The HB is a medium-sized (56 km 2 ) tidewater glacier located in the southern part of Wedel Jarlsberg Land. The glacier is ∼ 16 km long, and its elevation extends up to 550 m a.s.l. The WSB glacier has an area of 27 km 2 , is a land-terminating valley glacier to the west of HB, and ranges in elevation from 50 to 600 m a.s.l. (Schuler et al., 2020).
Each glacier was sampled in the ablation zone, close to the equilibrium line altitude (ELA), and in the accumulation zone ( Table 1). The ELA is the elevation at which the surface mass balance is zero, i.e. where the accumulation of snow is exactly balanced by ablation over a period of a year (Cogley et al., 2011). Although the exact elevation ranges of these zones (accumulation, ablation, and ELA) differ for each glacier, they share enough glaciological similarities to support intersite comparisons. All snowpits have been collected from the glacier central line in order to minimize the side accumulation effect due to orography.

Sampling procedure
Snowpit sampling was performed using a common protocol (Gallet et al., 2018) with pre-cleaned equipment (i.e. tubes, plastic scrapers, and plastic shovels cleaned with ultrapure water) and protective clothing (powder-free plastic gloves, clean suits, and face masks). This protocol allowed sampling and field data collection in a consistent manner, obtaining comparable datasets from different research sites.
Samples for ionic chemistry were taken from each discrete snowpit layer, according to the visible stratigraphy, and were directly filled into pre-cleaned 50 mL polypropylene "Falcon" centrifuge tubes. This type of sampling facilitates linking a snow layer (and its properties) to a specific weather event (i.e. precipitation or surface melt). Moreover, sampling by discrete layers makes it possible to correlate the intervals of snow accumulation between separate snowpits at different altitudes, as reported in this paper when we compare three different areas of the same glacier (ablation, ELA, and accumulation). It is also more accurate for chemical load calculations where ice layers occur in snowpits.
Samples for the isotopic composition of water were collected at a 5 cm resolution for sites in the Ny-Ålesund area and at a 10 cm or stratigraphic layer resolution for other sites. All sampling was conducted at a safe distance and upwind from potential local pollution sources, such as the snowmobiles used for transport by the sampling team.

Major ion analyses
Samples from glaciers in the Hornsund area (HB, WB) were analysed at the Polish Polar Station Hornsund (Institute of Geophysics, Polish Academy of Sciences), while samples from glaciers near Ny-Ålesund (KVG, ALB, HDF) were shipped frozen to the Institute of Polar Sciences (ISP-CNR) in Venice (Italy). Snow sampled in central Spitsbergen (AF, LF) was shipped frozen to the Department of Earth Sciences at Uppsala University (Sweden). Due to a temporary equipment malfunction in Uppsala, only cations could be analysed there, and the refrozen samples were forwarded to ISP-CNR for anion analysis. All samples and standards were handled and prepared under clean room conditions, wearing powderfree gloves. In all labs except at the Polish Polar Station Hornsund, laminar flow hoods (class 100) were used. Samples were melted immediately before analysis.

Hornsund
Samples were filtered through 0.45 µm mixed cellulose ester membranes (Merck Millipore S-pak ) prior to analysis. Ion concentrations were determined on a Metrohm 761 Compact IC ion chromatograph equipped with an autosampler (Metrohm, Herisau, Switzerland), with isocratic flow of 0.69 mL min −1 and chemical suppression for anions (column Metrosep A Supp S + Metrosep A Supp 4/5 Guard 4.0, eluent: NaHCO 3 1.0 mmol L −1 + Na 2 CO 3 3.2 mmol L −1 ). Cations were determined without suppression (column Metrosep C4 + Metrosep C4 Guard; mobile phase: HNO 3 1.7 mmol L −1 + 2,6-pyridinecarboxylic acid (dipicolinic acid, DPA) 0.7 mmol L −1 at a flow rate of 0.9 mL min −1 ). Cation samples were acidified with 2 µL of 2 mM HNO 3 per 10 mL sample prior to analysis. The injection volume was 20 µL in the anion system and 100 µL in the cation system. Nitric acid solutions were prepared from POCH S.A. (Poland) concentrated weighed amounts, while sodium carbonate and hydrogen carbonate as well as DPA were dissolved from the solid phase (Merck Millipore).

Uppsala
Samples were filtered using 0.22 µm polyethersulfone membranes (Minisart , Sartorius) and cation determination was performed using a Metrohm ProfIC850 ion chromatograph (Metrohm, Herisau, Switzerland), equipped with an autosampler and a Metrosep C4 column. The mobile phase of 0.02 M DPA and 0.1 M HNO 3 was run in an isocratic flow of 0.7 mL min −1 . Very low detection limits (≤ 0.006 mg L −1 ) were achieved thanks to the sample injection volume of 500 µL.

Instrumental performance of each laboratory
For all laboratories, calibrations for ions were evaluated using analytical standards (Merck/Sigma Aldrich). The calibrations in each lab delivered different linear ranges for each ion due to the different methods used (Table S1 in the Supplement). Good linearity was demonstrated in each lab and all calibration curves had R 2 > 0.99. Samples that had ion concentrations beyond the calibration range were diluted with ultrapure water before re-analysis. Analytical blanks of ultrapure water (> 18 M cm) were included in the analysis at all three labs. The method detection limit (MDL) was set to 3 times the standard deviation of the blank values (Table S1). For Na + , Mg 2+ , Cl − , and SO 2− 4 , values < MDL occurred in less than 10 % of cases, and for Ca 2+ and NO − 3 the < MDL concentrations were noted in 12 % and 17 % of all cases, respectively. However, K + and Br − were detected only in 53 % and 46 % of all samples, respectively, while NH + 4 concentration exceeded the MDL only in 36 % of all measurements. For the calculation of the bulk ionic loading in snowpits, values < MDL were assumed to be equal to half the MDL.
Accuracy and precision are important parameters to be evaluated during method validation. Checks for accuracy were made using certified multi-element standard solutions for anions (F − , Cl − , Br − , NO − 3 , SO 2− 4 , no. 89886-50ML-F, Sigma Aldrich) and cations (Na + , K + , Mg 2+ , Ca 2+ , no. 89316-50ML-F, Sigma Aldrich) at a concentration of 10 mg L −1 ± 0.2 %. Accuracy is expressed as a relative error calculated as (Q − T )/T × 100, where Q is the determined value and T is the "true" value. The accuracy for each ion in all labs was always < ±10 %, except for Mg 2+ measurements at the Hornsund laboratory. The analytical precision was quantified as the relative standard deviation (RSD) for replicates (n > 3) of standard solutions and was always < 10 % for each ion (Table S1).

Stable water isotopes
The determination of stable isotope ratios of O and H was performed at Tallinn University of Technology (Estonia). The isotopic ratios were determined by laser spectroscopy, using a Picarro model L2120-i water isotope analyser (Picarro Inc., Sunnyvale, USA), which allows for the simultaneous determination of 18 O/ 16 O and 2 H/ 1 H in H 2 O with a high-precision AO211 vaporizer. Results are reported in the standard delta notation as δ 18 O and δ 2 H relative to Vi-enna Standard Mean Ocean Water (VSMOW). Reproducibility was ±0.1 ‰ for δ 18 O and ±1 ‰ for δ 2 H, respectively. Seven injections were carried out for each sample, but only the last four injections (4 to 7) were used for calculations to minimize the memory effect. Laboratory standards TLN-A2 (−10.15; −77.5) and TLN-B2 (−21.95; −162.5) were regularly calibrated against international V-SMOW, GNIP, and V-SLAP standards. Standards (TLN-A2, TLN-B, and TLN-D4) were measured at the beginning, in the middle, and at the end of each set of measurements (54 bottles). Additionally, every seven samples, the laboratory standard TLN-D4 (−17.5; −133.0) was measured and used for drift correction if needed.

Spatial distribution of ionic species
To investigate differences in snowpack composition across all glaciers, we compared the total mass of ions that accumulated in snow at the different sampling sites. On average, the snow cover season on Svalbard lasts from early September to early May, but snow may also fall in summer months at high elevations. The snowpits in this study were sampled in early to late April 2016 and therefore might not contain the full annual ionic burden, since deposition can still occur before the beginning of the snowmelt season. Therefore, we report these data as ionic loads (mg m −2 ) rather than annual fluxes. In each snowpit, the total ionic load was calculated as the cumulative sum of the ionic loads in each discrete layer, i.e. ionic concentrations multiplied by the snow water equivalent of the layer (Table 2). On the other hand, to evaluate the transport processes of chemical species from other regions to Svalbard, we evaluate the volume-weighted mean concentrations of major ions in each snowpit. These values are calculated as the total ionic load of each snowpit divided by its total SWE (snow water equivalent) ( Table 3). The snowpack chemical characteristics were then compared between glacier zones (ablation zone, ELA, and accumulation zone; Fig. 1).
Snowpit samples collected in the Hornsund area (HB and WB, southern Spitsbergen) show a markedly higher total load for all major ions (Figs. 1 and 2) than at all other sites. The samples collected in the accumulation zones of WB and HB have total ionic loads of 8161 and 8023 mg m −2 , respectively, 4 times higher than those collected in the accumulation zone of KVG (2861), AF (2607), and ALB (1934) and 16 times higher than those sampled at LF (639) and HDF (583 mg m −2 ). Similar differences are observed for the snowpits collected at lower altitudes (Fig. 2).
In the accumulation zone of all glaciers (Fig. 3), Na + and Cl − are generally the most abundant ionic species, with percentages ranging from 29 % (HDF) to 36 % (AF) for Na + and from 34 % (LF) to 48 % (HB and WB) for Cl − , respectively. The snowpack on Hornsund glaciers (HB, WB) has higher Cl − percentages (48 %-49 %) compared to that of other glaciers (34 %-39 %), while conversely the SO 2− 4 percentage is lower there (9 %) than on other glaciers (11 %-23 %). The ionic loads are generally highest in the accumulation zone of glaciers and lowest in the ablation zone (Fig. 2), mostly due to the lower snow accumulation. This pattern holds true for Na + , Cl − , NH + 4 , K + , Ca 2+ , and Mg 2+ at most glacier sites, except in the Hornsund area. The load of Br − is similar on glaciers of the Ny-Ålesund sector (ALB, HDF, KVG) and on LF but is higher in AF and Hornsund glaciers (HB, WB; Fig. 2). The load of NO − 3 is similar for all glaciers, except for LF, where very low loads are found. Unlike total SO 2− 4 , the non-sea-salt fraction of sulfate (nss-SO 2− 4 ), calculated using a seawater SO 2− 4 : Na + mass ratio of 0.252 (Millero et al., 2008), shows lower loads on Hornsund glaciers (15-107 mg m −2 ) when compared to glaciers in other parts of the archipelago (Fig. 1, Table 2). The nss-SO 2− 4 loads vary between 22 and 131 at HDF and LF, between 75 and 266 at KVG and ALB, and between 153 and 206 mg m −2 at AF.

Stable water isotopes (δ 18 O and δ 2 H)
Our results provide the first picture of spatial variations in the mean stable water isotope composition of the seasonal snowpack across Svalbard (Table 3, Fig. 4). The SWEweighted mean δ 18 O and δ 2 H are calculated using the formula SWE − δ = (δ i × SWE i )/SWE t , where δ i are the δ values of each layer, SWE i are the SWE of each layer, and SWE t is the SWE of the entire snowpit. These SWEweighted mean values decrease significantly from south to north (Spearman rank correlation ρ with latitude is −0.69 and −0.65 for δ 18 O and δ 2 H, with p < 0.001 and p < 0.01, respectively). The isotopically heaviest snow (least negative δ values) occurs on glaciers of the Hornsund area (δ 18 O:  Glacier sites in north-western Spitsbergen (KVG, ALB, and HDF) and on LF have mean δ 18 O and δ 2 H values that fall within these ranges. On KVG, ALB, HDF, and LF, δ 18 O and δ 2 H values in snow decrease monotonically (becoming gradually more negative) with increasing elevation. On the other hand, on AF, WB, and HB there is no statistical difference between the mean δ 18 O and δ 2 H values of all snowpits (Fig. 4). A general significant anticorrelation with altitude is found for SWE-weighted mean δ 2 H (ρ = −0.63, p < 0.01) and δ 18 O (ρ = −0.65, p < 0.01).

Discussion
There have been few published studies on recent seasonal snow or firn chemistry in Svalbard, hence comparisons of our data with these earlier results are limited to a few sites. Virkkunen et al. (2007) and Vega et al. (2015) quantified the annual chemical loads of Na + , Ca 2+ , NO − 3 , and nss-SO 2− 4 at the Lomonosovfonna summit (LF3) from 2002 to 2011 using snow and firn cores. In Table 4, we also report the unpublished data of samples collected in 2009-2011 by C. Vega, obtained using the methods outlined in Sect. 2.2. Our study extends these data to 2016. The range of annual ionic loads at LF3 over the 15-year period is wide, and no clear temporal trend can be identified (Table 4). At the Holtedahlfonna sum-mit (HDF3), firn core measurements by Spolaor et al. (2013) found a mean Na + concentration of 110 ± 73 ng g −1 over the period 2003-2012, while the mean concentration in the April 2016 snowpack (this study) was 191 ng g −1 , hence within the range reported in earlier years.

The main ion sources in the seasonal snow of Svalbard
The composition of the Svalbard seasonal snowpack sampled during the C2S3 project clearly indicates that the ocean is the main source of ions in snow, as was shown by Hodgkins and Tranter (1998). At all sites, the dominant ions are Na + , Cl − , and SO 2− 4 , with comparatively minor amounts of K + , Ca 2+ , and Mg 2+ (Fig. 3). To help clarify the possible sources and modes of deposition of ions in snow, we computed Spearman rank correlations between total ionic loads (ρ load ), as well as between volume-weighted mean ionic concentrations (ρ conc ), across all snowpits (n = 22; Table 5). The chemical species that are predominantly wet-deposited, sharing common sources and not undergoing significant composition changes in transport should exhibit similar concentration patterns (high ρ conc ) (Schüpbach et al., 2018). The concentrations of Mg 2+ , K + , and Ca 2+ are all positively correlated with those of Na + and Cl − , indicating a common sea-spray source. Moreover, this input is the single significant source of K + and Mg 2+ , as indicated by near-zero calculated values of nss-K + and nss-Mg 2+ in the sampled snowpits (Table 3). The ρ load correlations are very similar for these ionic species, which points to both wet and dry deposition being a significant mechanism in their accumulation in the snowpack. The concentrations of Mg 2+ are positively and significantly correlated with both Ca 2+ and nss-Ca 2+ (ρ conc = 0.70 and 0.47, respectively; the latter coefficient is higher for loads at 0.56; Table 5), suggesting that they share some non-marine source(s). Furthermore, all glaciers have greater Ca 2+ : Mg 2+ ratios than seawater (0.32; Fig. 5; Millero et al., 2008). It is likely that the excess Ca 2+ and Mg 2+ come from mineral particles, i.e. CaCO 3 (calcite) and CaMg(CO 3 ) 2 (dolomite), derived from local rock (or soil) dust (Kekonen et al., 2005), especially limestone, dolostone, and marble, which are abundant in Svalbard (Dallmann, 1999). The presence of carbonate ions in the collected snow samples would explain the missing negative charge in the ionic balance (anion X − ; Fig. S1 in the Supplement).
Sulfate (SO 2− 4 ) is highly and significantly correlated (p < 0.05) with both Na + (ρ load = 0.92; ρ conc = 0.80) and Cl − (ρ load = 0.93; ρ conc = 0.75), indicating that sea spray is its main source (Table 5) deposition in Svalbard is long-range atmospheric transport of secondary aerosols containing SO 2− 4 , such as ammonium sulfate. This sulfate can be formed by SO x emitted from coal combustion throughout the winter and biomass burning in the spring (Barrie, 1986;Law and Stohl, 2007;Nawrot et al., 2016). The nss-SO 2− 4 does not correlate significantly with other ionic species, suggesting a separate origin. However, we need to caution that in the southern region of Svalbard, the estimation of ss-SO 2− 4 is subject to higher uncertainty because of the higher amount of Na + in the atmospheric deposition there.
Bulk ionic loads of SO 2− 4 in the snowpits are significantly and positively correlated with those of NO − 3 (ρ load = 0.55) and NH + 4 (ρ load = 0.68), but the correlations between weighted mean ionic concentrations are not significant, hinting at co-deposition (wet) rather than shared sources (Table 5). These species are known to form secondary aerosols (Karl et al., 2019;Schaap et al., 2004), and thus their proportions in aerosols may differ significantly from those in their source emissions. It is also possible that nitrogen species underwent further post-depositional photochemical reduction and evasion, thereby reducing their concentrations in snow Table 3. Volume-weighted mean concentrations of major ions in each snowpit (calculated as the sum of loads in all layers divided by the total SWE of the snowpit): the nss (non-sea-salt) fractions have been calculated in each layer before the volume-weighting procedure. Average SWE-weighted stable water isotope ratios (δ 18 O and δ 2 H expressed as ‰) and average deuterium excess (d) are also reported. Site  Table 4. Loads (mg m −2 ) of selected major ions from the 2016 sampling and from earlier studies at the Lomonosovfonna summit (LF3), corresponding to the concentrations given in Fig. 4 Table 5. Spearman rank order correlations of (a) ionic loads (mg m −2 ) and (b) SWE-weighted mean concentrations of major ions across all seven glaciers (n = 22 locations). ns: non-significant correlations (p value > 0.05). Ionic loads were calculated from all snowpit layers, while SWE-weighted mean concentrations were calculated by dividing the total ionic loads in each snowpit by its total SWE. Non-sea-salt (nss) components were estimated based on seawater ratios: Ca 2+ /Na + is 0.038, while SO 2− 4 /Na + is 0.252 (w w −1 ; Millero et al., 2008).  (Curtis et al., 2018). Finally, we remark here that the snowpit sampling was done in April, earlier than the beginning of the oceanic algal bloom in the surrounding Svalbard basin, which could have led to an underrepresentation of biological emissions from late spring in our samples. Spatial variations of ammonium load (NH + 4 ) across Svalbard glaciers mirror the pattern shown by sea-salt ions, with higher loads in the Hornsund area and lower loads in other areas. This is also reflected by significant correlations (Table 5) of the bulk loads of NH + 4 with those of Na + and Cl − (ρ load = 0.64 and 0.73, respectively) and with Na + , K + , and Mg 2+ by concentration (ρ conc = 0.47, 0.62, and 0.47, respectively). Ammonium has been linked to biogenic, forest fire, and anthropogenic agricultural emissions (Trachsel et al., 2019). The higher annual snowpack load of NH + 4 determined in the Hornsund area is more likely connected with biological sources than anthropogenic activities, although some contribution from biomass burning events cannot be excluded. The marine primary productivity in spring 2016 (April and May) was higher in the south-eastern ocean sector of the Svalbard archipelago (Fig. S2), which could partially explain the higher NH + 4 load. This would also explain the correlation between ammonium and sea-salt ions (Table 5). Locally, especially for HB, there may be extra NH + 4 emissions from bird colonies (Keslinka et al., 2019;Wojczulanis et al., 2008).
Unlike NH + 4 , the bulk loading of NO − 3 in snow is highest in north-western Spitsbergen (Ny-Ålesund area) when compared to other parts of Svalbard. Deposition of NO − 3 in Arctic snow is often associated with the long-range atmospheric transport of NO x and related N species from anthropogenic source regions at lower latitudes Fibiger et al., 2016;Vega et al., 2015). Differences in NO − 3 loads in snow in various parts of Svalbard might therefore reflect differences in the transport pathways of precipitating air masses, including the formation of secondary aerosols, or post-depositional processes, rather than local emissions. While local shipping routes and the settlement of Ny-Ålesund itself may contribute NO − 3 emissions (Winther et al., 2014), the highest share of the total ionic load of NO − 3 was found in the accumulation zone of HDF (9 % of the total ionic load; Fig. 3). Given that HDF is the most remote site from Ny-Ålesund relative to KVG or ALB, it should not capture a high share of local pollution. The highest correlation coefficient for NO − 3 , in terms of both concentrations and loads, was found with nss-Ca 2+ . This would support both the formation of calcium nitrate in the atmosphere (Gibson et al., 2006) or post-depositional processes removing the NO − 3 from layers poor in Ca 2+ , since calcium has been hypothesized to stabilize the nitrate in the snowpack against post-depositional losses (Kekonen et al., 2017).

Chlorine depletion
Although Na + and Cl − , the main species of sea salt, are significantly correlated (ρ conc = 0.95, Table 5), the values of the Cl − /Na + ratio (1.8 w w −1 ) in snow are lower than that in seawater on most studied glaciers, except those near Hornsund (Fig. 5), suggesting a Cl − deficit at the more northerly sites. Whillow et al. (1992) found an opposite situation in the snowpack of Greenland, indicating Cl − /Na + values higher than the ratio of seawater. This Cl − enrichment relative to the Cl − /Na + ratio in seawater may reflect Cl derived from anthropogenic sources as well as from gas-phase chlorine transportation and deposition in central Greenland.
Contrastingly, a possible explanation of Cl − deficit in the Svalbard snowpack might be de-chlorination of the seaspray aerosol during transport or, less likely, at the snowatmosphere interface. This reaction occurs between sea-salt particles, containing NaCl, and HNO 3 , H 2 SO 4 , or organic acids to release gaseous HCl (Zhuang et al., 1999). We calculated the percentage of Cl − depletion (Cl − dep ) as Cl − dep = (Cl − ss −Cl − meas )/Cl − ss ×100 %, where Cl − ss = 1.174Na + meas , and Cl − meas and Na + meas are the measured equivalent concentrations (Yao et al., 2003). Except for site HDF2 (Cl − dep = 2 %), the lowest mean Cl − dep values were obtained for Hornsund glaciers (WB,, while values at other glacier sites ranged between 21 % and 75 % ( Table 2). This suggests that sea-salt aerosols travel along a route from southern to northern Svalbard, which gives more time for Cl − depletion in the ionic mixtures reaching more northerly locations.

Bromine enrichment
In addition to Cl − , snowfall can scavenge Br − (Peterson et al., 2019;Spolaor et al., 2019). Br − loads on Svalbard glaciers surveyed in April 2016 are positively and significantly correlated with those of primary sea-salt ions Na + (ρ load = 0.48), Cl − (ρ load = 0.53), and Mg 2+ (ρ load = 0.51) ( Table 5). Correlations between weighted mean concentrations are not significant, however, suggesting departures of the Br − concentrations in snow from typical seawater ionic ratios at some glacier sites. A Br − enrichment factor (Br enr ) can be calculated as Br enr = Br − /(0.0065Na + ), where 0.0065 is the Br − : Na + seawater mass ratio (Maffezzoli et al., 2017). The Br enr reflects specific processes (in particular sea-ice Br emission) that affect the Br − concentration and load in the snowpack (Spolaor et al., 2014). Results of our calculations (Table 2, Fig. S3) show that for glaciers of the Hornsund area (HB and WB) and north-western Spitsbergen (KVG, ALB and HDF), the mean Br enr values are often < 1, indicating some Br − depletion, in agreement with the findings of Jacobi et al. (2019) for glaciers in the Ny-Ålesund area. The depletion could be a result of snowpack Br re-emission, but this seems unlikely since field measurements near Ny-Ålesund found no evidence of such volatilization of snow-bound Br (Spolaor et al., 2019).
Alternatively, Br − depletion could occur through BrO loss from marine aerosols and subsequent deposition of these Brdepleted aerosols in snow. In contrast to southern and northwestern Spitsbergen, glaciers in central Spitsbergen (LF) and in Nordaustlandet (AF) show Br enr values > 1. These glaciers lie relatively close to areas to the east of the archipelago that are often covered by first-year sea ice. Newly formed sea ice has been shown to release gas-phase Br into the polar atmosphere, thus supplying an extra Br source in addition to sea spray (Spolaor et al., 2016). The spatial distribution of the Brenriched snowpit sites supports this; i.e. the sites closest to the areas covered by first-year sea ice have the largest Br en-richments and the latter decrease with greater distance from the eastern shores of Svalbard (Fig. S3). A survey of the average sea-ice coverage in the period March-May 2016, which is relevant to the Br enrichment phenomenon, confirms that the north-eastern and eastern shores of Svalbard were indeed covered much more frequently by close and open drift ice than the south or north-west (https://cryo.met.no/en/sea-ice, last access: 1 March 2021).

Distribution pattern of δ 18 O and δ 2 H
As described earlier, the SWE-weighted mean δ 18 O and δ 2 H values in glacier snowpits decrease significantly with increasing latitude across Svalbard, the least negative values occurring on glaciers of the Hornsund area and the most negative in Austfonna (Table 3). This pattern follows the climate gradient across the archipelago, milder in the south, colder in the north. Part of the south-north contrast in δ values can be explained by the lower mean altitude of glacier sites in the Hornsund area compared to some of the higher-elevation sites further north on Spitsbergen or on Austfonna. The relationship with elevation is similar for both isotopic ratios in the collected dataset, except that with AF the isotopic signals might be influenced by additional processes since it is an isolated ice cap mainly surrounded by ocean or sea ice in winter.
Deuterium excess (d = δ 2 H − (8 × δ 18 O)) is mainly influenced by the source region of the precipitating moisture and in particular by the sea surface temperature, but also relative humidity and wind speed (Gat, 1996;Uemura et al., 2008). In addition, d is also influenced by the temperature gradient between the moisture source and precipitation area (Johnsen et al., 1989). The SWE-weighted mean d values in Svalbard snowpits vary within a relatively narrow range of 6.74 ‰ (from 10.10 ‰ to 16.84 ‰), and, similarly to δ 18 O, show no clear gradient with elevation or longitude (Table 3). Deuterium excess shows a significant correlation with latitude, at ρ = 0.60 (p < 0.01). A more detailed analysis of d by latitude shows that significantly different values are only obtained in snowpits sampled beyond 79.2 • N, i.e. in Austfonna snowpits. This is confirmed by the Kruskal-Wallis test, i.e. rank-based ANOVA, calculated with two groups of d values divided by the latitude threshold 79.2 • N (z = 4.23, p < 0.04). In fact, drawing the latitude threshold anywhere between 78.7 and 79.7 • N, a statistically significant difference with p < 0.05 is obtained. This is consistent with lower temperatures and evaporation rates in the more northern waters around Svalbard and suggests that snowfall on AF is at least partly affected by a different, more northerly moisture source than the rest of the archipelago.

Effect of elevation: a case study of Na
The glacier survey carried out during the C2S3 project afforded the opportunity to investigate the possible effect of Figure 6. Sodium load in snowpits ordered by increasing elevation in m a.s.l., indicated by the red lines. The colours identify areas where the snowpits have been excavated: each colour represents a separate glacier (HB -blue; WB -purple; ALB -red; LF -green; KVG -pink; AF -light blue; HDF -orange). IQR = inter-quartile range, i.e. the difference between the values of quartiles 3 and 1. elevation on the ionic composition of the snowpack. To do this, we compared the bulk load and SWE-weighted mean concentration of Na + across all studied snowpits, ordered by elevation (Fig. 6). Overall, both Na + load and concentration decrease with increasing altitude (ρ load = −0.24, not significant; ρ conc = −0.72, p < 0.05). This likely reflects greater local sea-spray aerosol deposition at lower compared to higher glacier sites. We then computed linear (Pearson) correlation coefficients (R, with associated p values) between log-transformed Na + loading (log(Na load )) and δ 18 O for all snowpits in the accumulation zones of all seven glaciers (Fig. 7). The calculation was performed with all snow layers. The Na + load was used as a sea-spray tracer, while the δ 18 O was assumed to vary with moisture source between discrete snowfall events. We find that the positive correlation between log(Na load ) and δ 18 O increases with elevation from R = 0.1 (HB3; 396 m a.s.l.) to R = 0.65 (LF3; 1193 m a.s.l.) and reaches a 95 % threshold of significance (R > 0.3) for glacier sites above 600 m a.s.l. (KVG, AF, LF, and HDF; Fig. 7). The average distance from the sea is a comparatively negligible factor in explaining the correlation between log(Na load ) and δ 18 O. The increase in strength and significance of the log(Na load )-δ 18 O correlation with altitude might be explained by different contributions of locally emitted ss-Na + relative to those of Na + from more distant sources. Sites located at lower altitudes are proportionally more affected by local sea-spray deposition, with or without snowfall. Conversely, sites at higher elevations likely receive a larger share of their ionic load from more distant sources and by wet deposition through snowfall. At the four sites (KVG, AF, LF, and HDF) where the log(Na load )-δ 18 O correlation is significant, increases in δ 18 O in snow layers are often associated with higher Na + concentrations. It is rather difficult to propose a precise explanation for this association. However, we would indicate that the isotopically heavier (less negative) δ 18 O values suggest that the co-registered Na + enhancements are associated with precipitation of a relatively warm air event, probably advected from lower latitudes. The snowfall associated with a warm event is able to wet scavenge the sea-spray aerosol present in the atmosphere. By contrast, when cold air masses (Arctic type) dominate, snowfall events are relatively limited due to the low air humidity causing a lower efficiency of wet scavenging. This results in lower δ 18 O and (likely) Na sodium loads, suggesting that wet deposition dominates the chemical load of the snowpack. Although this process should occur also at lower-elevation sites, the local emission and associated dry deposition are likely more important than wet deposition. More frequent meltrefreeze episodes at lower elevations would also mask the proposed relationship (as suggested by the vertical profiles of stratigraphy reported in Fig. S4).
Another possible explanation is that in the Arctic, air masses are transported from low-to high-elevation sites without any strong disturbance of the atmospheric conditions. In this case, isotopically heavier molecules and seaspray particles are gradually scavenged from the air masses. If this was the main process, we should find the correlation across all studied sites, assuming that Na + scavenged at a similar rate to that of isotopic fractionation. Since this has not been found, we propose that the correlation at higher elevation cannot be explained by atmospheric distillation alone. The possibility that the correlation is due to different sources of air masses seems unsupported due to the absence of a correlation between d-excess and sodium.

Summary and conclusion
We have quantified and described, for the first time, the spatial distribution of major ion loads (Ca 2+ , K + , Na + , Mg 2+ , NH + 4 , SO 2− 4 , Br − , Cl − and NO − 3 ) and variations of δ 18 O and δ 2 H in the snowpack on glaciers across Svalbard for a single accumulation season (2015)(2016). The highest total ionic loads are found in the southern region of Spitsbergen (Hornsund area) and exceed 8 g m −2 . Conversely, the lowest total ionic loads (≤ 0.6 g m −2 ) are found at sites in central or north-western Spitsbergen (LF and HDF). Sea-salt ions (Cl − , Na + , and SO 2− 4 ) dominate the ionic loads at all sites, but their share is highest at sites near Hornsund for e.g. 48 % Cl − compared to only 29 % on Holtedahlfonna. Relatively elevated Ca 2+ /Mg 2+ ratios in snow at all sites indicate nonsea-salt Ca 2+ inputs, most likely in the form of carbonate dust. Unlike other ions, NO − 3 has the highest loads in glaciers of north-western Spitsbergen and the lowest at LF. The nitrogen species, NO − 3 and NH + 4 , show distinct spatial distribution patterns. The highest NO − 3 loads are found in the northwestern part of Svalbard, while the highest NH + 4 loads are in the south-west. Bromide (Br − ) is most enriched in snow relative to seawater at AF and LF, the glacier sites located closest to areas with first-year sea-ice cover. This supports first-year sea ice being an important source of non-sea-salt Br − in the polar atmosphere.
An increasing positive correlation between log(Na load ) and δ 18 O as a function of elevation sites suggests that locations above 600-700 m a.s.l. are influenced by a proportionally higher share of ions from distant sources, while the lower sites are exposed to more local sources, especially sea spray. These findings confirm that the optimal sites to study the effects of long-range pollution deposition in Svalbard are those at higher elevations, such as the accumulation zones of HDF or LF, because they are the sites least impacted by local aerosol emissions. The current study gives the first picture of the ionic composition in the Svalbard snowpack in different regions across the archipelago, in the context of which processes are relevant in controlling the annual snowpack chemical composition, especially the influence of local and longrange transport.
Data availability. The data used in this manuscript are available from Table 2 and upon request to the authors.
Supplement. The supplement related to this article is available online at: https://doi.org/10.5194/acp-21-3163-2021-supplement. Competing interests. The authors declare that they have no conflict of interest.