Annual variability of ice-nucleating particle concentrations at different Arctic locations

Abstract

Number concentrations of ice-nucleating particles ($N_{\mathrm{INP}}$) in the Arctic were derived from ground-based filter samples. Examined samples had been collected in Alert (Nunavut, northern Canadian archipelago on Ellesmere Island), Utqiaġvik, formerly known as Barrow (Alaska), Ny-Ålesund (Svalbard), and at the Villum Research Station (VRS; northern Greenland). For the former two stations, examined filters span a full yearly cycle. For VRS, 10 weekly samples, mostly from different months of one year, were included. Samples from Ny-Ålesund were collected during the months from March until September of one year. At all four stations, highest concentrations were found in the summer months from roughly June to September. For those stations with sufficient data coverage, an annual cycle can be seen. The spectra of$N_{\mathrm{INP}}$ observed at the highest temperatures, i.e., those obtained for summer months, showed the presence of INPs that nucleate ice up to$-\mathrm{5}$ ${}^{\circ }$C. Although the nature of these highly ice-active INPs could not be determined in this study, it often has been described in the literature that ice activity observed at such high temperatures originates from the presence of ice-active material of biogenic origin. Spectra observed at the lowest temperatures, i.e., those derived for winter months, were on the lower end of the respective values from the literature on Arctic INPs or INPs from midlatitude continental sites, to which a comparison is presented herein. An analysis concerning the origin of INPs that were ice active at high temperatures was carried out using back trajectories and satellite information. Both terrestrial locations in the Arctic and the adjacent sea were found to be possible source areas for highly active INPs.

1 Introduction

The Arctic warms faster than any other region on Earth, a phenomenon known as Arctic amplification . Many different processes, some of which are heavily interconnected, contribute to this . However, not all of these processes and feedbacks are fully understood, and some might even still be unknown. Clouds in the Arctic are special in that they often form extended, persistent, low-level stratiform cloud layers, which are kept stable for days by different feedback processes (; ). These clouds influence the energy budget and generally warm the surface compared to clear skies . They often contain supercooled liquid water. In the range of temperatures ($T$) down to $-\mathrm{20}$ ${}^{\circ }$C, fractions of supercooled liquid clouds were reported to be above 50 % based on annual mean data for Europe and North America (both including the Arctic) from satellite remote sensing . For a multiyear analysis of all clouds, based on ground-based remote sensing at two western Arctic locations (Eureka and Utqiaġvik), clouds containing only liquid water occurred at least 20 % of the time in all months with a maximum of 56 % in September . Also, during two Arctic aircraft campaigns operating out of Inuvik, each in April and May of two different years, based on in situ measurements, at least 60 % of the clouds observed down to $-\mathrm{18}$ ${}^{\circ }$C were characterized as mostly liquid .

Ice nucleation forms primary ice in clouds, and for $T$ from 0 to roughly $-\mathrm{38}$ ${}^{\circ }$C, ice-nucleating particles (INPs) are needed to induce this nucleation process. Not many measurements of number concentrations of INPs ($N_{\mathrm{INP}}$) in the Arctic have been done up to now; however, these particles play an important role in the lifetime and radiative effects of Arctic stratiform clouds. A number of effects of ice in clouds are described in , including the fact that ice clouds are optically thinner than supercooled liquid water clouds so that the former emit less longwave radiation towards the surface. A modeling study showed that an increase in $N_{\mathrm{INP}}$ may cause a faster dissipation of these stratiform clouds , which in turn will influence the surface energy budget. But it should also be mentioned that $N_{\mathrm{INP}}$ values of $>\mathrm{1}$ L${}^{-\mathrm{1}}$ were needed to obtain an effect in this modeling study, and such concentrations were observed for midlatitude regions only for $T$ below $\approx -\mathrm{15}$ ${}^{\circ }$C . Recycling of INPs was assumed to be possible in Arctic clouds, again based on a modeling study ; i.e., ice crystals falling from a cloud could sublimate and re-entrain into clouds from below, which might potentially enhance the effect of changes in $N_{\mathrm{INP}}$. It also has been shown with large-eddy simulations that ice crystal number concentrations significantly influence cloud structure and the evolution of Arctic mixed-phase clouds . Overall, the cloud phase (i.e., supercooled water versus ice) is important for the radiation budget and hence the effect of Arctic stratiform clouds on climate. examined in detail a mixed-phase stratiform Arctic cloud and its phase transitions for roughly 1.5 d. Observed changes in the cloud were related to changes in air mass, but it was also said explicitly that for a better understanding of cloud phase transitions, among other observations, measurements of $N_{\mathrm{INP}}$ are also needed. All of this highlights the importance of insight on the abundance of INPs in the Arctic and on their sources and sinks.

In the past, some studies on Arctic INPs were done. However, most studies only included samples collected during short-term deployments. A comparison of results of the present study with some of those from the literature will be made further down in Sect. , while the main outcomes of these previous studies are already described in the following. In general, it can still be said that data on Arctic $N_{\mathrm{INP}}$ are scarce, which is particularly true for data at high $T$. derived $N_{\mathrm{INP}}$ based on ground-based and aircraft measurements, respectively. It was suggested that midlatitude pollution did not contribute INPs to Arctic aerosol, as $N_{\mathrm{INP}}$ values were found to be lowest in winter when Arctic haze, originating from anthropogenic pollution, was present. measured INPs during a ship cruise in the Arctic in the months from August to October. It was concluded that INPs were oceanic in origin, while land was only a weak source, and that the upper troposphere was deficient in INPs. derived $N_{\mathrm{INP}}$ also during an Arctic ship cruise in July to September and found a decline in $N_{\mathrm{INP}}$ during that phase. At least for some of the detected INPs, the most likely sources were assumed to be bacteria and fragments from marine biota emitted via bubble bursting from the open sea. detected INPs during aircraft measurements in the Arctic during the month of May. They reported strongly varying concentrations and found some INPs that contained Si that were likely mineral dust particles, while other INPs seemed to consist of low-molecular-weight components. also took aircraft measurements of $N_{\mathrm{INP}}$ in the vicinity of Arctic clouds during fall, but they reported lower values than those obtained by in spring. It was concluded that typical Arctic values for $N_{\mathrm{INP}}$ might be overestimated by current parameterizations. derived $N_{\mathrm{INP}}$ for size-segregated aerosol samples collected between the end of March and July in Alert. $N_{\mathrm{INP}}$ values derived from these Alert samples were slightly below values reported for other more southerly stations (mostly in Canada) in that study. They also found that at all stations large fractions of INPs were contributed by supermicron particles. These fractions were generally largest at the highest $T$ at which measurements were made, i.e., at $-\mathrm{15}$ ${}^{\circ }$C; for the Alert samples, $>\mathrm{90}$ % and 70 % of all INPs were $>\mathrm{1}$ and $>\mathrm{2.5}$ $\mathrm{\mu }$m, respectively. Similarly, found a size-dependent ability of particles nucleating ice for samples collected mostly in coastal areas in southern Canada and one sample collected in Lancaster Sound in the Canadian Arctic, with larger particles being more ice active. They also concluded that sea spray aerosol was not a major contributor to INPs for the samples taken in southern Canada. Based on concentrations of K-feldspar taken from a global model, $N_{\mathrm{INP}}$ values measured at $-\mathrm{25}$ ${}^{\circ }$C were modeled well, while INPs ice active at $-\mathrm{15}$ ${}^{\circ }$C were missing in this model. Also, reported a strong size dependence of the ice activity in samples collected on land in the northern Alaskan Arctic, where again the largest particles in the supermicron size range were the most efficient INPs. During sampling phases from March until mid-May, when grounds were covered in snow and ice, number concentrations of the supermicron INPs were lower than those in late May by up to 2 orders of magnitude. The increase in $N_{\mathrm{INP}}$ was suggested to originate from open leads in the sea ice and from open tundra. Similarly, for a coastal mountain station in northern Norway (at 70${}^{\circ }$ N), based on four filters sampled during July, observed that air masses were enriched in INPs ice active at $-\mathrm{15}$ ${}^{\circ }$C when they had passed over land. An origin of these INPs from decaying leaves was suggested. derived $N_{\mathrm{INP}}$ during a ship cruise in the Canadian Arctic marine boundary layer in summer. They suggest that mineral dust contributed more strongly to the observed INPs than sea spray, with mineral dust particles likely originating in the Arctic (Hudson Bay, eastern Greenland, northwest continental Canada).

Different substances are known to contribute to atmospheric INPs, as outlined in a number of review articles . In general, it is known that $N_{\mathrm{INP}}$ increases roughly exponentially with decreasing $T$, although at higher $T$ steep increases may be observed, followed by a weaker increase or even a plateau region down to roughly $-\mathrm{20}$ ${}^{\circ }$C as seen in, e.g.,. Ice nucleation at higher $T$ is typically related to macromolecules from biogenic entities as bacteria, fungal spores, lichen, pollen, and marine biota. These ice-active macromolecules nucleate ice from just below 0 ${}^{\circ }$C down to roughly $-\mathrm{20}$ ${}^{\circ }$C . Biogenic INPs typically occur in low concentrations in the atmosphere, but nevertheless, at remote marine locations such as the Southern Ocean, where $N_{\mathrm{INP}}$ is generally low, marine biogenic INPs might make up a large fraction or even the entire INP population . At less remote locations, the majority of atmospheric INPs consist of mineral dust particles originating from deserts or soils. Pure mineral dust particles of atmospherically relevant sizes typically are ice active below $-\mathrm{15}$ ${}^{\circ }$C or even below $-\mathrm{20}$ ${}^{\circ }$C and, with the abovementioned exception of remote marine locations, typically occur at much higher concentrations than biogenic INPs . However, mineral dust particles might also occur together with biogenic ice-active material , and such a mixed particle acts like a biogenic INP and should be attributed to the aforementioned group of biogenic INPs.

The existence of particularly high fractions of supercooled water observed in Arctic stratiform clouds, as, e.g., observed in in the temperature range above $-\mathrm{20}$ ${}^{\circ }$C in comparison to more convective clouds in midlatitudes and the tropics, could be expected to be linked to a lack of biogenic INPs in the Arctic due to sparse biological activity. However, it is known that biogenic INPs are contained in seawater and the oceanic surface microlayer (SML) and are emitted to the atmosphere by sea spray production . An increase in INP concentrations in the SML and for the biosphere in general from equatorial regions towards the poles has been observed. Also present in the Arctic are fungi , lichen, and bacteria , which could potentially contribute biogenic INPs.

In the present study, we aimed at increasing the knowledge of Arctic surface concentrations of INPs active in the immersion freezing mode, as described in the following. The immersion freezing mode was examined as it has been described as the most important heterogenous ice nucleation mode in mixed-phase clouds . No specific measurement campaign was organized for the examinations described herein. Instead, use was made of already existing filter samples. Besides for deriving temperature spectra of $N_{\mathrm{INP}}$ for 104 filter samples (Sect. ), we also determined possible sources for INPs that are ice active at high $T$ for selected samples. For that, correlations between $N_{\mathrm{INP}}$ and some chemical compounds were made (Sect. ), and an analysis was done concerning possible regions of origin of INPs that are ice active at high $T$ (Sect. ) for a selection of samples. The results will also be discussed in light of literature data (Sect. ).

2 Measurements

Quartz-fiber filters were sampled regularly at the four Arctic measurement stations of Alert, Ny-Ålesund, Utqiaġvik, and Villum Research Station (VRS) during the past years. Figure  shows the locations of these four stations, which are all in close proximity to the ocean ($<\mathrm{3}$ km). A portion of the filters was provided for the analysis presented herein. In the following, some detail will be given on these four different stations, including the filter handling, and on the measurement and evaluation method used to obtain INP number concentrations. We also describe in detail the temperature history of the filters, although it is not yet known with certainty how the temperature during storage will affect INP concentrations. Generally, the filters were kept frozen whenever possible. Transport from the four institutes where the samples had been kept to TROPOS was done in insulated boxes, together with cooling elements. The shipment was organized such that transport was fast (1–3 d) and that upon arrival at TROPOS the temperature in the boxes was still below 0 ${}^{\circ }$C. At TROPOS, samples were again stored at $-\mathrm{18}$ ${}^{\circ }$C until the measurements were done. These measures during storage and transport are precautions, as for biogenic INPs, storage at temperatures above 0 ${}^{\circ }$C or even storage under freezing conditions has been found to reduce their ice activity respectively. In this study, unless mentioned otherwise, all samples from all stations (also including field blanks) were treated similarly during all procedures. In the next sections, peculiarities of the separate four stations are described, followed by details of the measurements and their evaluation.

Figure 1

The location of the four stations from which filter samples were included herein: Utqiaġvik (red diamond), Alert (yellow diamond), VRS (green diamond), and Ny-Ålesund (blue diamond).

[Figure omitted. See PDF]

2.1 Alert

A custom-built high-volume aerosol sampler was used at the Dr. Neil Trivett Global Atmosphere Watch Observatory in Alert, Canada (82${}^{\circ }$30${}^{\prime }$ N, 62${}^{\circ }$22${}^{\prime }$ W; 210 m above sea level, a.s.l.), to collect 38 samples between April 2015 and April 2016. The sampler is installed at a walk-up deck about 4 m above the ground. The flow rate is approximately 1.4 m${}^{\mathrm{3}}$ min${}^{-\mathrm{1}}$ at standard temperature and pressure (STP) conditions. Quartz filters ($\mathrm{8}\times \mathrm{10}$ in; Pall Life Sciences, Pallflex filters, USA) were pre-fired at 900 ${}^{\circ }$C overnight and then shipped to Alert while already loaded on cartridges. During transport and storage, the filter-containing cartridges were wrapped with aluminum foil, and they were inside sealed plastic bags. Sampling time for those filters was either 1 week or 2 weeks, the latter being used from August until October (due to operational issues, no filter was sampled in July). A total of nine field blanks (roughly one every month) were collected. These field blanks were treated similarly to the other filters, i.e., inserted into the sampler for 2 min, but without an airflow through them. They were also stored similarly to the sampled filters at all times. After sampling, filters were stored (at room temperature $\approx \mathrm{20}$ ${}^{\circ }$C) in their sampling cartridges (wrapped in aluminum foil inside sealed plastic bags) at the Alert station and shipped in cardboard boxes (containing five sampling cartridges each) to the Toronto lab at Environment Climate Change Canada where they were stored frozen at $<-\mathrm{30}$ ${}^{\circ }$C. From these filters, a circular piece 47 mm in diameter was shipped to Leipzig for this study.

The total sampling area on the filters was 17.8 cm $\times$ 22.8 cm. For the measurements at TROPOS, described in detail in Sect.  below, circles with 1 mm diameter were punched out from the samples using sterile biopsy punches and immersed in ultrapure water separately. The volume of air sampled per 1 mm piece of filter differed for the different samples and varied from roughly 270 to 540 L.

2.2

Ny-Ålesund

Filter sampling in Ny-Ålesund on Svalbard (at 78${}^{\circ }$55${}^{\prime }$ N, 11${}^{\circ }$55${}^{\prime }$ E; 11 m a.s.l.) is done by the University of Florence, Italy. Quartz-fiber filters have been sampled regularly since 2010 using a high-volume sampler with quartz microfiber filters (CHMLAB Group QF1 grade, Barcelona, Spain). The filters were pretreated at 400 ${}^{\circ }$C prior to sampling. The filters had a diameter of 47 mm, one-quarter of which was provided for the present study. Sampling duration was 4 d in 2012, and of these filters, 13 sampled from late March until the beginning of September were examined in the present study, together with two field blanks. The total air volume collected on each filter was roughly 200 m${}^{\mathrm{3}}$. Each circular 1 mm filter piece used for the analysis sampled particles from roughly 130 L. Once sampled, filters were stored in a freezer at the Italian base in Ny-Ålesund and then shipped to Italy via cargo. At the home university, they were then stored in a cold room at $-\mathrm{20}$ ${}^{\circ }$C.

2.3

Utqiaġvik (formerly known as Barrow)

Filter sampling in Utqiaġvik, Alaska (at 71${}^{\circ }$18${}^{\prime }$ N, 156${}^{\circ }$46${}^{\prime }$ W; 11 m a.s.l.), is done by Baylor University, US, as described in . For the present study, quartz-fiber filters were used that had been sampled regularly in an annual campaign from June 2012 to June 2013 using a high-volume sampler (Tisch Environmental, Cleves, OH, USA). Filters were stored frozen prior to and immediately following all sampling. Two rectangular filter pieces (with a lateral length of 1.5 cm) from each of 41 different filters and two field blanks were provided for the present study. The sampled area of each filter was 399 cm${}^{\mathrm{2}}$. Sampling on each filter was done for 4 up to 13 d (7 d on average), collecting particles from a total air volume of roughly 14 000 to 51 000 m${}^{\mathrm{3}}$. This yields an air volume of roughly 270 to 1000 L collected on each circular 1 mm filter piece used for the analysis. Prior to sampling, filters were pretreated at 500 ${}^{\circ }$C. After sampling, the filters were stored in a freezer on-site and transported to the home university in coolers with cooling elements, where they were then stored at $-\mathrm{18}$ ${}^{\circ }$C.

2.4 Villum Research Station

Villum Research Station (VRS) at Station Nord in northern Greenland (at 81${}^{\circ }$36${}^{\prime }$ N, 16${}^{\circ }$40${}^{\prime }$ W; 24 m a.s.l.) is operated by Aarhus University, Denmark (in cooperation with the Danish Defense, the Arctic Command). Quartz-fiber filters have been sampled regularly since 2008 using a high-volume sampler (DIGITEL Hegnau, Switzerland) and employing weekly sampling . The filters had an exposed area of 154 cm${}^{\mathrm{2}}$ and sampled a total air volume of roughly 5000 m${}^{\mathrm{3}}$. From filters sampled in 2015, a 2 cm diameter piece was cut from each of the filters and provided for this study. Due to the large interest in shares of the filters, only samples from 11 different filters, all from different months in 2015 and 1 from December 2013, could be used herein. As for all samples used in this study, from the 2 cm pieces punches of 1 mm in diameter were cut at TROPOS directly prior to the measurements. The resulting small pieces were then used for INP analysis. The area of these 1 mm pieces corresponds to a sample volume of 255 L of air. Prior to sampling in the field, filters were pretreated at 450 ${}^{\circ }$C. Storage of the filters at VRS was done in freezers. Filters are transported from Greenland to Denmark around three times per year by the Danish Royal Air Force and then shipped to Roskilde, where they were then stored at $-\mathrm{18}$ ${}^{\circ }$C.

2.5 Freezing device INDA, the Ice Nucleation Droplet Array

For freezing experiments examining immersion freezing, a device comparable to one introduced in was used, but deploying PCR trays instead of separate tubes. The same device had been used in . From each filter piece that had been shipped to TROPOS, circles with a diameter of 1 mm were punched out directly before measurements were done, and each of the 96 wells of a PCR tray was filled with such a filter piece together with 50 $\mathrm{\mu }$L of ultrapure water (background measurements of ultrapure water are given in the Supplement). After sealing the PCR tray with a transparent foil, it was immersed into a bath thermostat such that the water table in the wells was below the surface of the liquid in the thermostat. The bath of the thermostat was then cooled with a cooling rate of 1 K min${}^{-\mathrm{1}}$, and the freezing process was monitored by a camera, taking a picture every 6 s. An LED light source installed below the PCR tray ensured that wells in which the water was still liquid could be easily distinguished from frozen ones. Frozen fractions ($f_{\mathrm{ice}}$) were then determined as the number of frozen tubes divided by the total number of tubes. Typically, fresh water to be used in the experiments was taken once a day and stored in a glass bottle. Whenever fresh water was taken, an experiment was run with this water in the tubes only to ensure that the water was satisfyingly clean. Similarly, experiments were run with field blank filters that had the same history as the samples but without sampling (see the Supplement), and signals from the field blanks were well below those of the sampled filters. A subtraction of the signals of the field blank from those of the measurements was not done. This is justified in a detailed discussion in the Supplement. The interpretation of the results from the filters presented in this study is the same for both uncorrected and background-corrected samples.

2.6

Deriving $N_{\mathrm{INP}}$

Equation () was used to derive $N_{\mathrm{INP}}$ from the measured $f_{\mathrm{ice}}$ . This equation accounts for the possibility of the presence of multiple INPs in one vial by assuming that the INPs are Poisson distributed. Additionally it normalizes the values resulting from the measurement with the air volume sampled on each 1 mm filter piece. This yields concentrations of ice-nucleating particles per volume of sampled air. 1$$N_{\mathrm{INP}}=-\left(\ln (\mathrm{1}-f_{\mathrm{ice}})\right)/\left(F\cdot A_{\mathrm{p}}/A_{\mathrm{f}}\right)$$ $F$ is the total volume of air drawn through the filter, and $A_{\mathrm{p}}$ and $A_{\mathrm{f}}$ are the surface area of a single 1 mm filter piece and the whole sampled area of the filter, respectively.

The temperature and concentration regions for which data were obtained for the different samples depend on a number of factors. The measured value $f_{\mathrm{ice}}$ is a fraction ranging from 0 to 1. Therefore, $N_{\mathrm{INP}}$, as derived using Eq. (), can only take on a limited range of values. This range is based on the negative natural logarithm of 0.01 and 0.99 (4.6 to 0.01). The absolute values of $N_{\mathrm{INP}}$ then also depend on the volume of air sampled onto each 1 mm filter piece, i.e., on the volume of air drawn through the filter during the sampling period and on the relation of the surface area of one filter piece to the total sampled surface area. The volume collected per 1 mm filter piece was within a factor of 4.5 for all filters (120 to 540 L). Altogether, the range of $N_{\mathrm{INP}}$ that can be obtained herein is roughly from $\mathrm{2}\times {\mathrm{10}}^{-\mathrm{5}}$ to 0.04 L${}^{-\mathrm{1}}$. From this limitation, it also follows that the range of $T$ for which $N_{\mathrm{INP}}$ could be obtained is limited, as it is tied to the concentrations that can be measured. Times with more ice-active INPs show up as $N_{\mathrm{INP}}$ at higher $T$. The highest $T$ at which ice activity was observed was close to $-\mathrm{5}$ ${}^{\circ }$C, as will be shown in the next section.

It should also be noted that samples that had less than 60 L of air volume collected on each 1 mm filter piece were also examined, but $f_{\mathrm{ice}}$ was close to the background and therefore these samples were not considered in this study. Results from background measurements are given in the Supplement. Measurement uncertainty as shown in this work was derived based on , i.e., following the assumption that the INPs are Poisson distributed between the different examined droplets. The first few droplets that freeze in each experiment therefore show the highest uncertainties.

2.7 Using back trajectories and satellite maps

A more in-depth analysis concerning possible INP source regions was done for a selection of filter samples from each measurement station. For that, 5 d back trajectories were calculated with HYSPLIT to determine the origin of sampled air masses. These calculations were based on GDAS (Global Data Assimilation System) meteorological data using an hourly time resolution. A new trajectory was started every 6 h during the whole time for which sampling was done on the respective filter. Back trajectories were initiated at an altitude of 100 m above the sampling locations, as this altitude still has a high likelihood of being connected to the ground and as lower elevations are more prone to uncertainties. In the Supplement, these back trajectories are shown separately for the selected examined filter samples.

Using these back trajectories, we examined over which ground the air masses collected on the filters had passed in the 5 d prior to arrival at the measurement station. The aim was to see contributions from open land or open water for the different selected filter samples. Therefore, a distinction was made between snow, open land, sea ice, and open water. To do so, maps from the Interactive Multisensor Snow and Ice Mapping System (IMS) 2008 were used. IMS maps are a composite product produced by NOAA/NEDIS (National Oceanic and Atmospheric Administration's National Environmental Satellite Data and Information Service) combining information on both sea ice and snow cover. Information from 15 different sources of input is included in the production of these maps . These maps have been provided for 20 years. We used the daily Northern Hemisphere maps with a resolution of 4 km 2008. For each time step we applied nearest-neighbor interpolation in space and time to find the corresponding satellite coordinate along the back trajectory. With that, for each back trajectory, we determined the conditions on the ground during the passage of the air mass, i.e., if the ground was covered by snow or ice or if open water or open land was present. The resulting information is shown exemplarily in Fig.  for the filter samples collected in Alert starting 10 June 2015. Using these maps, we counted how often air masses that were collected on one filter were above open land or open water while the air mass was below 100 m. (As will be discussed in detail in Sect. , for Utqiaġvik, results reported here always refer to an upper altitude of 500 m and 10 d back trajectories.) An altitude restriction was used as we were trying to geographically locate INP sources on the surface. Additionally, the back trajectories were only considered back in time until an integral amount of 2 mm of precipitation (taken from the information included in the back trajectories) was reached. This was done as precipitation formation occurs via the ice phase so that precipitation is assumed to lead to a washout of INPs.

3 Results

In the following, $N_{\mathrm{INP}}$ derived from filter samples will shortly be introduced. A correlation with some available chemical composition data is made. Finally, an analysis of air mass origins is introduced to analyze possible source regions for INPs that are ice active at high $T$. A comparison to literature data and further discussion of the results are then presented in Sects.  and , respectively.

3.1 Arctic atmospheric INP concentrations

Quartz-fiber filters from the four different Arctic stations shown in Fig.  were analyzed to derive $N_{\mathrm{INP}}$. Figure  shows $N_{\mathrm{INP}}$ for all different samples separately for the four sampling locations. Due to the comparably large number of filters analyzed for Alert and Utqiaġvik, separate curves cannot be seen easily in Fig. . Therefore, Fig.  shows time series of $N_{\mathrm{INP}}$ for $T$ at $-\mathrm{7}$, $-\mathrm{10}$, $-\mathrm{13}$, and $-\mathrm{9.5}$ ${}^{\circ }$C for Alert, Utqiaġvik, Ny-Ålesund, and VRS, respectively. $T$ was chosen such that $N_{\mathrm{INP}}$ could be obtained for the largest possible number of all curves (data for the three samples with the lowest $N_{\mathrm{INP}}$ are missing for Alert (29 April and 27 May 2015 and 4 April 2016), as is the one with the highest $N_{\mathrm{INP}}$ for Utqiaġvik (3 May 2013), indicated by arrows in Fig. ). The yellow background shows for which samples an more in-depth analysis is presented in Sect. . Error bars in Fig.  show the 95 % confidence interval.

Figure 2

$N_{\mathrm{INP}}$ for all samples. Symbol color distinguishes between data obtained for the different months, while symbol type indicates the day of the month when sampling started. Curves with particularly high and low $N_{\mathrm{INP}}$ are explicitly listed in the legend. (The $T$ axis is the same for panels a–c and different for panel d; see values for $T$ given on the top and bottom of the figure, respectively.)

[Figure omitted. See PDF]

Figure 3

Time series of $N_{\mathrm{INP}}$ at $T$ of $-\mathrm{7}$, $-\mathrm{10}$, $-\mathrm{13}$, and $-\mathrm{9.5}$ ${}^{\circ }$C for Alert, Utqiaġvik, Ny-Ålesund, and VRS, respectively. Small arrows in the panels indicate times when values were either below (for Alert) or above (for Utqiaġvik) the detection limit. The yellow background shows for which samples a more in-depth analysis is presented in Sect. . Error bars show the 95 % confidence interval.

[Figure omitted. See PDF]

It is worth noting that once a sample has a comparably high concentration at one $T$ this is generally observed at all $T$ at which measurements are available and vice versa; i.e., the curves do not intersect much (see Fig. ). Therefore, the curves shown in Fig.  can be used to discuss observed trends for INPs that are ice active at high $T$. It should also be noted that for Ny-Ålesund, data only exist for March until September and that for VRS there is mostly only one data point per month, if any (no data exist for February, March, and May).

Figure 4

Exemplary results from the analysis of the IMS maps for the filter collected in Alert starting 10 June 2015. In both panels, each row (indicated as track number) represents one trajectory going back in time for 5 d, starting from 0 (when sampling took place), and displaying 120 separate time steps. The colors indicate the altitude of the air mass for the different time steps in panel (a) and the respective nature of the ground at the location of the air mass in panel (b).

[Figure omitted. See PDF]

In Figs.  and it can be seen that in general $N_{\mathrm{INP}}$ obtained for the summer months is higher than for the winter months. A decrease in $N_{\mathrm{INP}}$ is observed starting in fall (October or November). For months with the lowest observed ice activity, which are generally winter and early spring months, values for $N_{\mathrm{INP}}$ were measured down to below $-\mathrm{20}$ ${}^{\circ }$C. From June until September mostly INPs that were ice active between $-\mathrm{5}$ and $-\mathrm{15}$ ${}^{\circ }$C were detected (see Fig. ), and for Utqiaġvik and VRS such highly ice-active INPs were observed as early as April. These highly ice-active INPs will be the focus of the next two sections (Sect.  and ). Similarly highly ice-active INPs have been suggested to be biogenic in origin based on tests such as heat treatment , which due to the limited available amount of filter material could not be done in the present study.

3.2 Sources of INPs

3.2.1 Correlation to chemical composition

Typical $N_{\mathrm{INP}}$ values measured in the atmosphere are several orders of magnitude below total particle number concentrations, and therefore mass concentrations of INPs are so small that a correlation between bulk chemical composition and $N_{\mathrm{INP}}$ might not be expected, particularly not for the very rare INPs that are ice active at high $T$. This is in line with recent findings for a long-term study of INPs at Cape Verde by . There, no correlation between $N_{\mathrm{INP}}$ and bulk chemical composition was found for $T$ down to $-\mathrm{16}$ ${}^{\circ }$C for a number of different compounds, which included ${\mathrm{Ca}}^{\mathrm{2}+}$, ${\mathrm{Na}}^{+}$, and elemental carbon as tracers of continental, marine, and combustion sources, respectively. At a lower $T$ of $-\mathrm{25}$ ${}^{\circ }$C, recently reported that mineral dust tracers correlated with INPs, which suggests that mineral dust was a major contributor to the INP population at that $T$. found that for three coastal sites in Canada a model based on K-feldspar as the only INP calculated $N_{\mathrm{INP}}$ that fit measurements well at $-\mathrm{25}$ ${}^{\circ }$C, while at $-\mathrm{15}$ ${}^{\circ }$C measurements were underpredicted, suggesting a missing source of INPs that are active at higher $T$. In the following we examine whether there is a correlation between the INPs detected in the present study that are ice active at high $T$ and chemical composition.

The examined filters pieces were not particularly sampled for this study and were entirely needed for the above-described INP analysis. No dedicated chemical analysis could be done additionally. But as other parts of most of the filters were also used in other studies, some information on chemical composition was available. This was used to derive the correlations with $N_{\mathrm{INP}}$ shown in Fig. , i.e., with those INPs that are ice active at high $T$. Table  shows values for $R$, $R^{\mathrm{2}}$, and $p$ for linear correlations between $N_{\mathrm{INP}}$ and different bulk chemical properties.

Table 1

$R$, $R^{\mathrm{2}}$, and $p$ values for linear correlations between $N_{\mathrm{INP}}$ (as shown in Fig. 2) and different bulk chemical properties.

Location	Species	$R$	$R^{\mathrm{2}}$	$p$
Utqiaġvik	EC	$-\mathrm{0.36}$	0.13	0.02
	OC	$-\mathrm{0.04}$	$<\mathrm{0.01}$	0.82
	fluoride	0.14	0.02	0.41
	chloride	0.13	0.02	0.43
	nitrite	0.22	0.05	0.19
	bromide	$-\mathrm{0.01}$	$<\mathrm{0.01}$	0.98
	sulfate	$-\mathrm{0.12}$	0.01	0.47
	nitrate	0.03	$<\mathrm{0.01}$	0.86
Ny-Ålesund	PM10	$-\mathrm{0.36}$	0.13	0.22
	OC	$-\mathrm{0.39}$	0.15	0.18
	ammonium	$-\mathrm{0.44}$	0.20	0.13
	K	$-\mathrm{0.57}$	0.32	0.04
	Mg	$-\mathrm{0.36}$	0.13	0.22
	nitrite	0.15	0.02	0.62
	nitrate	$-\mathrm{0.16}$	0.03	0.59
	sulphate	$-\mathrm{0.60}$	0.36	0.03
	Na	$-\mathrm{0.30}$	0.09	0.32
	Ca	$-\mathrm{0.50}$	0.25	0.08
	Cl	$-\mathrm{0.13}$	0.02	0.68
	MSA	0.18	0.03	0.55
Alert	POC+CC	0.59	0.35	0.01
	OC	$-\mathrm{0.12}$	0.01	0.62
	EC	0.05	$<\mathrm{0.01}$	0.84

In general, no correlations were found. The only case with a positive value for $R$ (0.59) and a low value for $p$ (0.01) was found for POC+CC (pyrolyzed organic carbon and carbonate carbon) in Alert. Both POC and CC contain carbon that pyrolyzes at 870 ${}^{\circ }$C in a pure He stream . CC might indicate the presence of soil dust . POC includes some charred carbon formed at 550 ${}^{\circ }$C, which is the lower temperature step of the applied analysis EnCan-total-900 method;, and highly oxidized organic compounds and/or high-molecular-weight refractory carbon. Based on previous studies, the POC mass is proportional to the oxygen mass in organic aerosols , releasing as carbon monoxide at 870 ${}^{\circ }$C. POC was observed to form from sucrose and glucose and therefore might be indicative of biogenic material. This points towards a direction in which more detailed studies should be undertaken in the future. In Sect. , we will discuss a range of possible sources for the observed INPs.

3.2.2 Determination of possible source regions

Results from the more in-depth analysis concerning possible INP source regions, based on back trajectories and satellite maps as described in Sect. , are presented in the following. For this analysis, samples were chosen that had been collected in spring, directly before and after the transition from typical winter to typical summer conditions (see yellow background in Fig. ). A total of 17 separate filter samples were included: 2 for VRS, 4 each for Alert and Ny-Ålesund, and 7 for Utqiaġvik. The aim was to see if contributions from open land or open water were potentially more pronounced during times when INPs active at high $T$ were observed.

Figure  shows the number of time steps when air masses were over open land or open water for the separate filter samples. Additionally, gray bars in the background indicate the percentage of time the air masses collected on one filter were below 100 m for Alert, VRS, and Ny-Ålesund. It can already be seen that the presence of highly ice-active INPs on a filter is related to air masses that fulfill the above criteria, i.e., that traveled over open land or open water at a low altitude. It also can be seen that this was not found for Utqiaġvik. Initially, no open land and hardly any open water had been found for this site when 5 d back trajectories were used, together with an altitude restriction of 100 m, which means that air masses did not travel over open land or open water at altitudes below 100 m. To check if the length of the back trajectory or the chosen maximum altitude influenced our results for Utqiaġvik, an analysis was also done using 10 d back trajectories and 500 m as the altitude limit, which is presented in Fig. . This extension only resulted in larger percentages of time for which the air masses were below this altitude limit. However, there were still not a large number of time steps found for which air masses traveled over open land or open water for Utqiaġvik. We will get back to this again below.

Figure 5

The number of time steps when air masses were at low altitudes and over open land or open water are shown in green and blue, respectively, for different filter samples. The analysis shown was done for altitudes up to 100 m for Alert, Ny-Ålesund, and VRS and up to 500 m for Utqiaġvik. Gray bars in the background indicate the percentage of time the air masses collected on each filter were below that altitude. Triangles indicate samples for which highly ice-active INPs were detected, i.e., for which INP spectra were measured at high $T$ (single-colored triangles) and at medium $T$ (triangles with yellow interior). The respective INP spectra are shown in Fig. .

[Figure omitted. See PDF]

Figure 6

Panels (a), (b), (c), and (d) each show a map on the left side in which locations are indicated where air masses collected on different filters crossed over open land or open water while being at a low altitude (below 100 m for Alert, Ny-Ålesund, and VRS; below 500 m for Utqiaġvik). Black diamonds indicate the location of the measurement site. (Please note: the maps in c and d are rotated by 90${}^{\circ }$ compared to the those in a and b.) The right side shows the INP spectra for the corresponding filters.

[Figure omitted. See PDF]

Figure  shows the spectra of $N_{\mathrm{INP}}$ (called INP spectra for simplicity from now on) for the samples included in this analysis (right side) and the locations where the respective air masses traveled over open land or open water at altitudes below 100 m (or 500 m for Utqiaġvik) (left side). Three different types of INP spectra can be distinguished: first, there are those for which we observed the start of ice activation only at around $-\mathrm{10}$ ${}^{\circ }$C and which went down to well below $-\mathrm{15}$ ${}^{\circ }$C. For these, INP spectra and locations are depicted in magenta or orange. Second, there are INP spectra for which we observed ice nucleation from roughly $-\mathrm{5}$ to above $-\mathrm{15}$ ${}^{\circ }$C. For these, INP spectra and locations are depicted in greenish colors. The third category was used only for Utqiaġvik for INP spectra with medium ice activity depicted in blueish. Error bars shown in Fig.  show the 95 % confidence interval.

For Alert, VRS, and Ny-Ålesund, the absence or scarcity of orange and magenta marks on the maps in Fig.  (maps on the left in a, b, and c) shows that almost no open land or open water contributed to air masses sampled on the respective filters. This is in accordance with the corresponding INP spectra, which showed comparably low ice activity. The magenta locations close to Svalbard for the Alert sample correspond to the sample from 20 May 2015, for which somewhat more ice-active INPs were found than for the subsequent sample from 27 May 2015. For this latter sample, no contributions from open land or open water were observed, and it is, in fact, the sample with the lowest ice activity observed in this study (see Fig. ). Locations depicted in greenish colors potentially contributed INPs that are ice active at high $T$. They can be found on open land as well as on open water. In connection to filters sampled at Alert or at VRS they show up in north Greenland, on Ellesmere Island (on which Alert is located), in Baffin Bay, and along the southern part of the west coast of Greenland. Concerning filters sampled at Ny-Ålesund, greenish marks show up on Svalbard and the adjacent sea.

The above analysis shows that coastal regions may be particularly important as a source for highly ice-active INPs, including open waters close to coasts. Indeed, highly ice-active biogenic INPs were found in Arctic surface waters before e.g.,. For the highly ice-active samples collected on Ny-Ålesund on 12 and 28 June 2012, the surroundings of the measurement station were completely snow free during the times when these samples were collected, whereas for all other cases there was at least partial or total snow cover around the stations. In other words, local terrestrial sources close to the measurement station may also contribute as sources for highly ice-active INPs, as already discussed in . Also, describe Arctic landmasses to be the source for observed Arctic INPs (ice active at $-\mathrm{15}$, $-\mathrm{20}$, and $-\mathrm{25}$ ${}^{\circ }$C), and these INPs were suggested to be mineral dust. On Svalbard, found higher atmospheric $N_{\mathrm{INP}}$ in July than in March, and they additionally described glacial outwash sediments in Svalbard to be highly ice active. This ice activity was assumed to be connected to small amounts of organic (likely biogenic) material. Based on these findings, suggest the higher $N_{\mathrm{INP}}$ in summer to be connected to organic (biogenic) components in glacially sourced dust. Some coastal regions in the Arctic, e.g., the west coast of Greenland together with the region around Baffin Bay and the Canadian Arctic Archipelago as well as the area around the Bering Strait and also Svalbard, are known for their abundance of seabird colonies . These regions partially coincide with regions highlighted as possible INP sources in Fig. . These regions are known to emit ammonia, which plays a role in new particle formation in the Arctic . But clearly, newly formed particles are not expected to contribute to atmospheric INPs at the temperatures examined in this study, and INPs are likely also emitted from regions with high biological activity. In Sect.  we will discuss possible INP sources in more detail.

For Utqiaġvik, data from seven filters were included in the analysis. Two of them showed INP spectra at comparably low $T$, three at medium $T$, and two at high $T$. For all types of INP spectra, no contribution from open land was observed with the back-trajectory analysis. Only minor contributions from open water were found for the latter two types, although the analysis was extended to include 10 d back trajectories, and the maximum altitude up to which air masses were considered was relaxed to 500 m. Air masses did travel below 100 m, and even more often below 500 m (see Fig.  and the back trajectories for Utqiaġvik and their heights profiles in the Supplement). However, the transition to filters on which INPs active at comparably high $T$ was already found to happen earlier at Utqiaġvik than at the other three measurement locations towards the end of March. IMS maps almost exclusively identified the ground as sea ice and snow in the regions that were crossed by the air masses, even until the beginning of May 2013. And in general, air masses spent more time over sea ice than over snow (see back trajectories in the Supplement). There has to be a source for highly ice-active INPs that was not revealed in the analysis done here. Polynyas and open leads may contribute to explaining this inconsistency. The resolution of the IMS maps used may be too coarse so that open water related to polynyas and open leads could have gone unnoticed.

While the analysis introduced here shows regions that may have potentially contributed highly ice-active INPs to the sampled air masses, it does not make any statement about other regions. Other regions could potentially be sources, too, but might have only been crossed by air masses at high altitudes or may not have been crossed at all.

4 Comparison with literature

Figure  shows the ranges of $N_{\mathrm{INP}}$ observed for the four stations as shaded areas. Data from the different stations cover a rather similar range. As explained above (Sect. ), $N_{\mathrm{INP}}$ could only be measured up to some 10${}^{-\mathrm{2}}$ L${}^{-\mathrm{1}}$, depending on the volume of air sampled onto one 1 mm filter piece. Hence, the upper concentration limit of our data is determined by the measurement method. The gray background shows literature data of $N_{\mathrm{INP}}$ determined from precipitation samples collected mostly in North America and Europe . In that data set, samples showing the highest $N_{\mathrm{INP}}$ at $T>-\mathrm{20}$ ${}^{\circ }$C originate from rain and hail samples collected in North Carolina (US) and Alberta, Montreal (CA). Some of the INP spectra we detected at the highest $T$, observed particularly in Alert and Utqiaġvik, show values for $N_{\mathrm{INP}}$ that are similar or only about 1 order of magnitude below data reported in . At the lowest temperatures at which we detected INP spectra, $N_{\mathrm{INP}}$ values are lower than data from , i.e., the lowest Arctic $N_{\mathrm{INP}}$, as those we observed in the winter months might be below the lowest values observed on continents in midlatitudes. It is worth adding that still lower concentrations were observed in marine remote locations in the Southern Ocean and for clean marine air in the northeast Atlantic .

Figure 7

Comparison of $N_{\mathrm{INP}}$ determined in this study for the Arctic with literature data by (gray background), , , , , , , , and . Green and brown symbols represent data from surface-based measurements; black and blue represent airborne measurements. For , brown indicates data they cited from the literature, with the vertical bar indicating the extent of the reported values.

[Figure omitted. See PDF]

Figure  shows additional data on Arctic $N_{\mathrm{INP}}$ from the literature that were already discussed in the Introduction . Not every single data point from these papers is shown, as Fig.  aims to give an overview of the range of data that exists. Data on Arctic $N_{\mathrm{INP}}$ in general are still scarce, which is particularly true for data at high $T$. Also, data scatter over a wide range. The highest values of $N_{\mathrm{INP}}$ shown in Fig.  originate from one set of aircraft measurements made in May , while a second set of aircraft measurements, taken in September and October onboard an aircraft flying out of Alaska, agrees with our data at the highest $T$ . For the data taken from , the reported measurement uncertainty is shown, which is representative for typical uncertainties for the type of instrumentation used in and . Error bars indicate 1 standard deviation at the higher end. The lower end is indicative of the detection limit, and for a substantial fraction of measurements no INPs were detected in .

Going back to Fig. , literature data from ground-based measurements that are in the same range of $N_{\mathrm{INP}}$ in which we measured are also within the same $T$ range. But besides the data from , these data are at the lower end of $T$ that we observed. This also holds for the data from , which were obtained during summer ship cruises in Baffin Bay and in the central Bering Sea.

Data shown for were taken in Ny-Ålesund and Utqiaġvik. A slight tendency towards higher $N_{\mathrm{INP}}$ in summer months, compared to winter months, can be seen. Similarly, as said in the Introduction, found a decreasing trend for $N_{\mathrm{INP}}$ at $-\mathrm{15}$ ${}^{\circ }$C from July to September. The observed decrease was roughly 1 order of magnitude; however, the scatter from sample to sample was almost as large as that trend. Nevertheless, these data sets are early indications of the annual trend that has very clearly been found in the present study.

Concerning possible sources for INPs, assumed that mainly oceanic sources contributed to the observed INPs, with only a weak contribution from land. discussed a marine origin of at least some of the INPs they analyzed. These latter two studies were ship based. The land-based study by showed INPs that were ice active at higher $T$ than those observed in and . traced these INPs back to terrestrial contributions, possibly decaying leaves. During aircraft measurements, identified some INPs as mineral dust particles and others as containing low-molecular-weight components. These latter might have been connected to biogenic INPs.

These different studies report very diverse $N_{\mathrm{INP}}$. In these studies, different instrumentation was used and sometimes different ice nucleation mechanisms were probed. Also, instrumental limitations typically determine the ranges of $T$ and $N_{\mathrm{INP}}$ that can be probed. All of this might add to the diversity in the data. But, as can be seen in Fig. , a difference of up to 2 orders of magnitude in $N_{\mathrm{INP}}$ was measured between winter and summer for a single temperature in our data set. Therefore, the diversity in $N_{\mathrm{INP}}$ reported in previous studies will also originate from different times of the year when these studies were conducted.

All of our samples were collected on land, and regions that showed up as possible sources for highly ice-active INPs in Fig.  were on or close to land. While, as said above (Sect. ), other regions cannot be excluded as sources, it will be interesting to see in the future if $N_{\mathrm{INP}}$ values detected further away from terrestrial sources will consistently be lower than those obtained on land.

Concerning an influence of INPs emitted from the ground at higher altitudes, recently compared aerosol particle number size distributions (PNSDs) measured on the ground and during overflights at different heights in May 2014 in Tuktoyaktuk (Northwest Territories of Canada on the Arctic Ocean). PNSDs measured on the ground and at heights up 1200 m were generally similar, while PNSDs obtained at higher altitudes were clearly different. Therefore, atmospheric aerosol, including INPs, can be similar at heights up to levels at which cloud formation is observed. Hence, INPs detected by ground-based measurements may well be able to influence ice formation in clouds, at least during times when the cloud layers are coupled to the surface.

Compared to previous literature data introduced in this study, the new long-term data presented here extend the range of Arctic $N_{\mathrm{INP}}$ towards higher $T$. They also clearly show that Arctic INP concentrations vary throughout the year, with a regular presence of INPs that are ice active at $T$ well above $-\mathrm{10}$ ${}^{\circ }$C throughout the summer months at all four terrestrial measurement stations.

5 Discussion

The annual cycle observed for $N_{\mathrm{INP}}$ in this study is not in tune with what is known for particle number concentrations and size distributions occurring across the Arctic . This is not too surprising: recent studies, including Arctic locations, found that a large fraction of INPs are supermicron in size , while the majority of particles are in the submicron size range.

Concerning the annual Arctic aerosol cycle, there is a maximum in particle number concentrations in early spring, before precipitation sets in, caused by accumulating anthropogenic pollution known as Arctic haze . $N_{\mathrm{INP}}$ values are low during that time of the year, which might indicate that anthropogenic pollution does not contribute to atmospheric INPs, at least in the $T$ range examined in this study. This is in line with the observation that Arctic haze particles are not efficient INPs and also with a recent study showing that anthropogenic pollution did not contribute to INPs in very polluted air in Beijing . It is also in agreement with observations by based on ice cores from Svalbard and Greenland, who found that $N_{\mathrm{INP}}$ in the Arctic did not increase over the past 500 years (from roughly 1480 to 1990), while tracers for anthropogenic pollution did increase markedly.

The formation of Arctic haze is related to the fact that during winter months, air masses from midlatitudes can transport aerosol particles into the Arctic . In contrast, in the summer months the Arctic lower atmosphere is effectively isolated and the transport of atmospheric aerosol particles into the Arctic is low . Hence, the increase in $N_{\mathrm{INP}}$ observed here in late spring and summer has to originate from Arctic sources, including local ones situated close to the measurement site.

Once Arctic haze is removed by precipitation in spring , Arctic particles are mostly newly formed particles in the size range up to 100 nm , originating from gaseous precursors . New particle formation ceases in late summer or early fall, and until Arctic haze sets in again later in winter, the total particle number concentration is at its lowest at well below 100 cm${}^{-\mathrm{3}}$. But newly formed particles do not contribute to INPs unless subject to particular conditions at cirrus level . Nor does their annual cycle follow what we observed for $N_{\mathrm{INP}}$, as concentrations for newly formed particles decrease much earlier in the year than observed for $N_{\mathrm{INP}}$.

Arctic concentrations of mineral dust particles have been measured and modeled and were described to largely depend on long-range transport, with largest concentrations occurring in spring. However, dust from local sources may also contribute . Mineral dust particles themselves likely only become ice active at $T$ at the lower end of or below those examined in the present study, but they may be carriers of biogenic ice-active macromolecules . Continental sources of fungal spores were found to contribute to the organic aerosol observed in the Arctic in summer , and it is known that spores of some fungal species are ice active at comparably high $T$ . The same also holds for lichen and bacteria . Sources and abundance of bacteria were examined for an area in southwest Greenland , and roughly half of all airborne bacterial cells were described to originate from local terrestrial environments such as surface soils, while the other half was said to have been long-range transported, originating from marine, glaciated, and terrestrial surfaces. In summary, microorganisms originating from continental sources in general could contribute to the highly ice-active INPs observed in this study. This is in line with the observation presented in Fig. , in which, particularly for VRS and Ny-Ålesund, highly ice-active INPs were observed for air masses that traveled over open land.

Particle generation at open leads in the Arctic has been observed in relation to sea spray from bubble production mechanisms that exist independently of wind . found particles in Arctic aerosol that could be attributed to algal exopolymer secretions and suggest they became airborne via bubble bursting. This agrees with observations by and , showing that Arctic marine microgels from algae, present both in surface water and SML, contribute to Arctic atmospheric particles. Such microgels were assumed to be related to INPs observed in marine SML . It was recently found that ice activity in Arctic seawater was negatively correlated with salinity , possibly indicating that these aquatic INPs were associated with melting sea ice releasing ice-active biogenic material into the ocean. Marginal ice zones are known to be of importance for phytoplankton blooms due to melting sea ice releasing iron . Overall, although the abovementioned particle production mechanism in Arctic open leads was found to be only of minor importance for the overall atmospheric particle number concentration , this, together with wind-driven sea spray production, could be a source for at least some of the highly ice-active INPs observed in our measurements. This likely holds for Alert, but it may also be the case in the example of Utqiaġvik, for which in April and May very few contributions from open land or open sea, as defined by the IMS maps, were found, while highly ice-active INPs were present. These could potentially originate from open leads or polynyas. It should be added that found a similarly large increase in highly ice-active INPs during May at Oliktok Point in northern Alaska, only roughly 300 km east of Utqiaġvik. This increase was related to INPs from tundra surfaces and open water, particularly the marginal ice zone. Polynyas 700 km away from the measurement site were not found to contribute, likely due to settling of particles during the long traveling times in the slow-moving air masses that were observed.

In summary, while anthropogenic pollution and new particle formation do not explain the highly ice-active INPs observed in this study, these INPs can originate from both terrestrial and marine sources in the Arctic. These sources are strong in summer and weak or absent in winter, depending on the conditions on the ground. Should terrestrial sources be found to contribute more strongly than marine ones, $N_{\mathrm{INP}}$ values further away from land will be lower than those observed here.

Concerning Arctic mixed-phase clouds at temperatures above $-\mathrm{20}$ ${}^{\circ }$C, recently observed a depression of cloud ice in these clouds when they are polluted, as observed during Arctic haze. Under these conditions, they have a lower amount of cloud ice mass for a given amount of condensed liquid mass. Based on the results presented here, this can be explained by the lower concentrations of INPs during times of Arctic haze – a reasoning that sheds light on the importance of the findings reported here, which could not yet have been discussed in . With respect to recent modeling results, did large-eddy simulation (LES) studies of Arctic mixed-phase stratocumulus clouds and found indications that changes in INPs in Arctic aerosol dominate over changes in cloud condensation nuclei (CCN). This shows the potential importance of determining INPs and their changes due to climate change in the Arctic. However, examined the Arctic annual cycle in cloud amounts from 24 different models and found significant disagreements. They conclude that the parameterization of ice microphysics in the models contributes to the observed differences. Overall, the field of Arctic mixed-phase clouds and related ice nucleation is currently one of intense research, and the results presented herein contribute to our understanding and lay a foundation for improved descriptions of INPs in modeling in the future.

6 Summary and conclusions

This study shows a yearly cycle of $N_{\mathrm{INP}}$ at four different Arctic sites, an observation that has not been seen so clearly in previous studies. Maximum values were observed roughly from late spring until well into fall, and minimum values were observed during winter and early spring. INP spectra for the most ice-active INPs are close to some observed at continental locations outside of the Arctic . Potential source regions for INPs that are ice active at the higher examined $T$ were determined based on combining air mass back trajectories with IMS satellite maps. Such source regions were found in the Arctic on open land as well as on open water, particularly in Baffin Bay, along the southern part of the west coast of Greenland, and on Svalbard and the adjacent sea. Contributions from these regions could explain the increase in highly ice-active INPs observed in spring in Alert, Ny-Ålesund, and VRS. However, a possible source region for the highly ice-active INPs observed in Utqiaġvik could not be identified. Here, highly ice-active INPs appeared earlier in the year than at the other stations. Open leads and polynyas present in the Arctic sea ice are not identified as open water in IMS maps and may have contributed to these INPs.

Regions that were not identified as source regions for highly ice-active INPs in this study may still contribute such INPs. It should be mentioned that source regions determined here were mostly on or close to land. We showed that there is a large scatter in the literature values observed for $N_{\mathrm{INP}}$ in the Arctic. As we find a yearly cycle in Arctic $N_{\mathrm{INP}}$, this large scatter may partially originate from the different times of the year when samples were taken. It may also depend on the proximity to land where the sampling was done if terrestrial sources and sources in close proximity to land are shown to be dominant for highly ice-active INPs in the future.

Independent of their origin, these observed INPs might be transported certain distances and might have more than only local influences. This can become important, as biogenic emissions in general can be expected to increase in the Arctic in coming years due to Arctic amplification, which amplifies marine primary production and changes Arctic microbial communities . Feedback mechanisms involving ice formation, along with radiative properties and lifetimes of Arctic stratiform clouds, may exist and be related to biogenic INPs. Therefore, more thorough studies concerning Arctic INPs are needed. Determining the current status of INPs in the Arctic and future changes that are to be expected might help us to understand aerosol–cloud interactions in the Arctic and their significance for the observed strong Arctic warming.

Data availability

Data (frozen fractions and INP concentrations) are available on Pangaea (https://doi.pangaea.de/10.1594/PANGAEA.899701; ).

The supplement related to this article is available online at: https://doi.org/10.5194/acp-19-5293-2019-supplement.

Author contributions

HW and FS designed the research; RJS, TEB, RT, SB, RB, HS, LH, WZ, and HH did field campaign management and/or filter sampling; TEB, CEM, and WZ did chemical analysis; MH calibrated and JL, ML, and OL operated the ice nucleation measurement device; PH calculated trajectories; AH did the work on the satellite data; HW analyzed the data and wrote the paper with contributions from LH, HH, RJS, CEM, RT, FS, MH, PH, ML, RB, HS, and WZ.

Competing interests

The authors declare that they have no conflict of interest.

Special issue statement

This article is part of the special issue “Arctic mixed-phase clouds as studied during the ACLOUD/PASCAL campaigns in the framework of (AC)${}^{\mathrm{3}}$ (ACP/AMT inter-journal SI)”. It is not associated with a conference.

Acknowledgements

We gratefully acknowledge the funding from the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation; project number 268020496, TRR 172) within the Transregional Collaborative Research Center “ArctiC Amplification: Climate Relevant Atmospheric and SurfaCe Processes, and Feedback Mechanisms (AC)${}^{\mathrm{3}}$” for filter analysis and interpretations done in subproject B04. Heike Wex was partially funded by DFG within the Ice Nuclei research UnIT (INUIT, FOR 1525), WE 4722/1-2. We acknowledge Cecilia Shin, Darrell Ernst, and Andrew Platt for logistic support, fieldwork in Alert, and assistance in analysis, as well as the Canadian Forces Station Alert for supporting sample collection. Logistic assistance of the Polar Support Unit of the CNR (Italian National Research Council) Department of Earth and Environment (POLARNET) in coordinating the activities based at the Dirigibile Italia Arctic station at Ny-Ålesund is acknowledged, specifically Mauro Mazzola and Angelo Viola. Financial and technical support for samples from Utqiaġvik was provided by the United States Department of Energy (Atmospheric Radiation Measurement field campaign no. 2010-05876), NOAA (award no. NA14OAR4310150), and the C. Gus Glasscock, Jr. Endowed Fund for Excellence in Environmental Sciences. We also thank Walter Brower and Jimmy Ivanoff of the Ukpeaġvik Iñupiat Corporation for sample collection and field assistance and Fred Helsel, Dan Lucero, Jeffrey Zirzow, and the Sandia National Laboratory for site access and preparation. The work at Villum Research Station at Station Nord was financially supported by the Danish Environmental Protection Agency and Energy Agency via the MIKA/DANCEA funds for Environmental Support to the Arctic Region, which is part of the Danish contribution to the Arctic Monitoring and Assessment Programme (AMAP) and to the Danish research project “Short-Lived Climate Forcers” (SLCF). The Villum Foundation is acknowledged for funding the construction of Villum Research Station, Station Nord. The authors gratefully acknowledge the NOAA Air Resources Laboratory (ARL) for the provision of the HYSPLIT transport and dispersion model and/or READY website (http://www.ready.noaa.gov, last access: 13 April 2019) used in this publication. And furthermore, we thank Richard Leaitch (ECCC) for helpful discussions during the writing process.The publication of this article was funded by the Open Access Fund of the Leibniz Association.

Review statement

This paper was edited by Amy Solomon and reviewed by Paul DeMott and one anonymous referee.