Measurement report: Introduction to the HyICE-2018 campaign for measurements of ice-nucleating particles and instrument inter-comparison in the Hyytiälä boreal forest

Abstract

The formation of ice particles in Earth's atmosphere strongly influences the dynamics and optical properties of clouds and their impacts on the climate system. Ice formation in clouds is often triggered heterogeneously by ice-nucleating particles (INPs) that represent a very low number of particles in the atmosphere. To date, many sources of INPs, such as mineral and soil dust, have been investigated and identified in the low and mid latitudes. Although less is known about the sources of ice nucleation at high latitudes, efforts have been made to identify the sources of INPs in the Arctic and boreal environments. In this study, we investigate the INP emission potential from high-latitude boreal forests in the mixed-phase cloud regime. We introduce the HyICE-2018 measurement campaign conducted in the boreal forest of Hyytiälä, Finland, between February and June 2018. The campaign utilized the infrastructure of the Station for Measuring Ecosystem-Atmosphere Relations (SMEAR) II, with additional INP instruments, including the Portable Ice Nucleation Chamber I and II (PINC and PINCii), the SPectrometer for Ice Nuclei (SPIN), the Portable Ice Nucleation Experiment (PINE), the Ice Nucleation SpEctrometer of the Karlsruhe Institute of Technology (INSEKT) and the Microlitre Nucleation by Immersed Particle Instrument ($\mathrm{\mu }$L-NIPI), used to quantify the INP concentrations and sources in the boreal environment. In this contribution, we describe the measurement infrastructure and operating procedures during HyICE-2018, and we report results from specific time periods where INP instruments were run in parallel for inter-comparison purposes. Our results show that the suite of instruments deployed during HyICE-2018 reports consistent results and therefore lays the foundation for forthcoming results to be considered holistically. In addition, we compare measured INP concentrations to INP parameterizations, and we observe good agreement with the parameterization developed from measurements conducted in a ponderosa pine forest ecosystem in Colorado, USA.

1 Introduction

Atmospheric aerosols are recognized to play an important role in nearly every aspect of the physics and chemistry of the atmosphere . However, the interactions between aerosols and clouds and how these interactions influence Earth's surface energy budget and water cycle represent significant knowledge gaps in climate science. Thus, it is critical to understand fundamental aerosol processes such as aerosol formation, growth and aerosol–cloud interactions to evaluate the impact of aerosols on Earth's radiative balance.

Certain types of atmospheric particles, called ice-nucleating particles (INPs), have the potential to initiate the formation of ice in clouds, thus affecting the properties of the clouds and often the initiation of precipitation. Heterogeneous ice nucleation processes may include (i) deposition nucleation, where bulk liquid water is presumed to be absent and where ice is formed from vapor supersaturated with respect to ice; (ii) immersion freezing, where ice formation is initiated by an INP located within a body of liquid; (iii) condensation freezing, where an INP simultaneously acts as a cloud condensation nucleus (CCN); (iv) contact freezing, where ice formation is triggered at the air–water interface by an INP that comes into contact with a supercooled liquid droplet; and (v) pore condensation and freezing, where ice is formed via liquid water condensation in pores .

Numerous measurements of INPs have been performed in various settings and environments. However, measurements in boreal forests are largely underrepresented, and little is known concerning the INP sources and properties from this environment. A recent study showed that the boreal forest is an important source of biogenic INPs and that the seasonal cycle of INP concentrations is linked to the abundance of biogenic aerosol particles . Boreal forests represent more than one-third of all forests and cover more than 15 $\times {\mathrm{10}}^{\mathrm{6}}$ km${}^{\mathrm{2}}$ of land . The majority of boreal landscape is in the Arctic and sub-Arctic region of the continental Northern Hemisphere, and due to its high-latitude position, the environment experiences strong seasonal changes in meteorological and environmental conditions. Moreover, the boreal forest environment is generally far from large dust sources and strong anthropogenic emissions, which motivates the investigation of biogenic ice nucleation activity in this environment. Finally, INPs from mid- to high-latitude sources may have a disproportionate effect on climate through their influence on shallow clouds , emphasizing the need to know more about INP sources in these regions, including in boreal forests.

To overcome our lack of knowledge, a measurement campaign was organized in 2018 at the Station for Measuring Ecosystem-Atmosphere Relations (SMEAR) II, at the Hyytiälä Forestry Research Station in Juupajoki, Finland. The intensive field campaign, called HyICE-2018, began in February and extended until June, with longer-term INP monitoring efforts continuing for more than 1 year beyond the intensive measurement period .

One aim of the campaign was to utilize the highly instrumented and time-resolved measurements of the SMEAR II station. The station is a well-known atmospheric observatory in a high-latitude location (61${}^{\circ }$51${}^{\prime }$ N) within a boreal forest ecosystem . The SMEAR II site has hosted field campaigns for decades and is well suited for investigating potential links between atmospheric and ecosystem processes. The station is equipped with a suite of advanced aerosol instrumentation that has been utilized for significant advancements in aerosol process studies. Therefore, the SMEAR II station is well suited for an intensive campaign focused on measuring INPs in the boreal forest environment.

Additionally, new particle formation (NPF) events are frequently recorded at the SMEAR II station , and there are relevant scientific questions that remain unanswered concerning the role of NPF in the context of ambient INPs. NPF events occur when volatile organic compounds (VOCs) and biogenic volatile organic compounds (BVOCs) are oxidized and condensed to form very small particulate, secondary organic aerosol (SOA), some of which continue to grow to larger sizes that may form CCN and INPs . The role of NPF in forming particles large enough to participate in cloud activation as CCN and INPs has been investigated in laboratory studies . However, while there is evidence that NPF can play a role in CCN activation , and studies have shown that SOA can nucleate ice under cirrus conditions , the relationship between NPF and INPs has not been explored with field measurements.

In order to augment the standard SMEAR II monitoring instrumentation for the HyICE-2018 campaign, several institutions deployed INP measurement systems, including the University of Helsinki, the University of Gothenburg, the Karlsruhe Institute of Technology, the University of Leeds, the University of Eastern Finland and the Swiss Federal Institute of Technology in Zurich (ETH Zurich). Furthermore, the Finnish Meteorological Institute and Tampere University contributed remote sensing retrievals and bioaerosol monitoring, respectively. In addition to the individual data collection efforts, several days were utilized for instrument inter-comparison analyses to test the reproducibility of results across instruments and scientific teams in the field setting.

The primary objectives of the campaign were to

• •

  quantify and characterize INPs in a boreal environment within different thermodynamic forcing regimes (i.e., different temperatures, $T$, and saturations with respect to ice, $S_{\mathrm{i}}$),

• •

  determine and quantify the existence of INPs seasonal variation,

• •

  assess the vertical distribution of INPs above the boreal forest.

This paper gives an overview of the campaign setting and design as an introduction to the Copernicus special issue “Ice nucleation in the boreal atmosphere”, which is anticipated to include several contributions from HyICE-2018. Data from several days dedicated to evaluate instrument inter-comparison is presented to illustrate instrument-to-instrument agreement and to facilitate future presentation and interpretation of data from subsets of instruments.

2 Methods

2.1 Measurement site – SMEAR II

The HyICE-2018 campaign took place at the SMEAR II station in Hyytiälä, Finland . The station is located in southern Finland (61${}^{\circ }$51${}^{\prime }$ N, 24${}^{\circ }$17${}^{\prime }$ E; 181 m above sea level) and is surrounded by boreal coniferous forest (Fig. ). The conditions at the site are typical for a background location, with the main pollution sources being the city of Tampere (60 km to the southwest with $\approx \mathrm{238}\mathrm{000}$ inhabitants as of 2019) and the activity and buildings at the station . The station has several operational units that span a wide forest area and reach into and above the tree canopy. There are several towers, including a 128 m mast used for atmospheric flux measurements, an 18 m tower for irradiation and flux measurements, a separate 18 m tower for tree physiology measurements, and a 35 m walk-up tower for aerosol measurements. There are also several measurement cottages and containers spread throughout the forest, with a cottage located near the walk-up tower dedicated to measuring the physical properties of aerosols, and a main cottage based near the 128 m mast (Fig. ).

Figure 1

Overview of the SMEAR II station in Hyytiälä, Finland, used for staging the HyICE-2018 measurement campaign, with the available infrastructure and measurement platforms depicted.

[Figure omitted. See PDF]

The SMEAR II station is equipped to monitor the physical and chemical properties of aerosols and gas-phase precursors to aerosol formation with a suite of state-of-the-art monitoring instrumentation. Measurements also cover meteorology, radiation, soil, snow cover and gases. An overview of the instruments in operation at the site during the HyICE-2018 campaign is available in the Appendix (Tables ,  and ). In this study, data from the SMEAR II differential mobility particle sizer (DMPS) and aerosol particle sizer (APS; TSI model 3321) were used. The DMPS measures aerosol particle size distributions from 3 to 1000 nm in mobility diameter, with a 10 min time resolution . During HyICE-2018, the instrument was sampling through a total suspended particulate (TSP) inlet 8 m inside the forest canopy and was operated following the guidelines from Aerosols, Clouds and Trace gases Research InfraStructure ACTRIS;. The APS is used to measure the supermicron aerosol particle size distribution from 0.5 to 20 $\mathrm{\mu }$m in aerodynamic diameter. The instrument was sampling through a TSP inlet (DIGITEL Elektronik GmbH) 6 m above ground level, and a vertical sampling line was used to avoid particle losses. In addition, the APS sampling line was heated to 40 ${}^{\circ }$C to ensure that the relative humidity in the sampling line remained below 40 %, which prevents condensation and dries the aerosol before measurements.

2.2 Additional INP measurements for the HyICE-2018 campaign

One motivation for HyICE-2018 was to compare different INP measurement instruments and techniques in a field setting. Previous studies have performed inter-comparisons of INP instrumentation in a number of intensive measurement campaigns . However, these efforts focused on well-controlled laboratory measurements to assess sampling procedures and to calibrate instruments relative to one another . Only a few INP instruments have been colocated for long field measurements of real atmospheric aerosol , and continued efforts in multi-instrument measurements can be beneficial for the entire community .

During HyICE-2018, several instruments and techniques were used to quantify INPs, as summarized in Table . The working principles of the different INP instruments are depicted in Fig. . Three continuous-flow diffusion chambers (CFDCs, described in Sect. ) and one expansion chamber (described in Sect. ) were used to conduct online measurements of INP concentrations, and two droplet-freezing assays (described in Sect. ) were used for offline measurements of INP concentrations collected on filters. The intensive measurement campaign took place between February and June 2018, and the dates on which the instruments were operating are depicted in Fig.  against a background of the evolving average daytime temperature, which illustrates the changing seasons.

Table 1

Instrumentation installed at the SMEAR II research station specifically for the HyICE-2018 campaign. A list of the instruments permanently running at the SMEAR II research station can be found in Tables , and , including the APS and DMPS which were used in this analysis.

Function	Type	Instrument	Institution	Inlet	Time resolution	Sampling location
INP counters	Continuous flow	PINC	ETH	PM${}_{\mathrm{2.5}}$ + heated inlet	20–30 min NC	Main cottage
	diffusion chamber			+ dryer
	(CFDC)	PINCii	INAR, GU	PM${}_{\mathrm{2.5}}$ + heated inlet	15 min NC	Main cottage
		SPIN (DMT)	UEF	PM${}_{\mathrm{2.5}}$ + dryer	30 min NC	Container
	Expansion chamber	PINE	KIT, Leeds	Heated inlet	6 min continuous	Main cottage
	Droplet freezing assay	INSEKT	KIT	PM${}_{\mathrm{10}}$	10–24 h	Aerosol cottage,
						35 m tower,
		$\mathrm{\mu }$L-NIPI	Leeds	PM${}_{\mathrm{1}}$, PM${}_{\mathrm{2.5}}$, PM${}_{\mathrm{10}}$	4–12 h	Aboard aircraft
Particle	Condensation	CPC 3010	INAR	TSP	Continuous	128 m mast
counters	particle counter	(TSI)
		CPC 3010	INAR	TSP	Continuous	35 m tower
		(TSI)
	Optical particle sizer	OPS 3330	INAR	TSP	Continuous	128 m mast
		(TSI)
	Bioaerosol sensor	WIBS-NEO	TAU	$\ge \mathrm{0.5}\mathrm{\mu }$m	Continuous	Cottages
		(DMT)

NC = not continuous (continuous = measuring without interruption unless maintenance was required); PM = particulate matter; TSP = total suspended particulate.

Figure 2

Timeline of deployment for the instrumentation installed specifically for the HyICE-2018 campaign. Shading depicts mean daytime air temperature ${\stackrel{\mathrm{‾}}{\text{T}}}_{\text{day}}$ (08:00–20:00, UTC+2) measured at 4.2 m above ground level. Instrument inter-comparison days are indicated by the dashed lines.

[Figure omitted. See PDF]

Figure 3

Schematic of the working principles of the INP measurement systems deployed during HyICE-2018. (a) The CFDCs operate with activation and evaporation sections and alternate ambient sampling (pink) and filtered background measurements (green). (b) The PINE expansion chamber samples a volume of ambient air before a short expansion cycle activates INP within the trapped volume. In both (a) and (b), optical particle counting is used at the chamber exits to measure particle number. (c) The droplet-freezing assays rely on filter sampling of ambient aerosols. Collected particles are washed from the filters and analyzed for INP content, generating curves of INP temperature spectra. (d) The thermodynamic space typically accessible to different measurement techniques.

[Figure omitted. See PDF]

2.2.1 Continuous-flow diffusion chambers

Three CFDCs were deployed for online INP measurements (Fig. a). These instruments include the Portable Ice Nucleation Chamber (PINC), the Portable Ice Nucleation Chamber II (PINCii) and the SPectrometer for Ice Nuclei (SPIN). The three models are iterations of a design that consists of two parallel ice-coated walls which are cooled to below freezing temperatures . In the upper part of the chamber, referred to as the main chamber, the walls are held at different temperatures (Fig. a), to produce vapor supersaturation in the region between the walls . The lower chamber, or evaporation section, is held at an isothermal condition of ice saturation. During measurements, a continuous flow of sample air, referred to as the sample lamina, is sandwiched between two particle-free sheath flows and is drawn through the chamber, thus exposing any airborne particles to the supersaturated conditions. At a given lamina temperature, $T_{\mathrm{l}}$, and saturation condition (e.g., saturation with respect to ice, $S_{\mathrm{i}}$), some particles will induce water condensation and/or ice formation by accommodating excess vapor. When the flow leaves the main supersaturated chamber, it immediately enters an evaporation section, which is held at ice saturation, and is thus subsaturated with respect to liquid water. Within the evaporation section, liquid droplets evaporate, creating a size difference between ice crystals and water droplets. When paired with an inlet size cutoff, particles within the exit flow that exceed the cutoff are determined to be ice. They are detected and counted by an optical particle counter (OPC) located at the exit of the chamber. Above a certain $S_{\mathrm{i}}$, CFDCs have an instrument-specific point of “droplet breakthrough” where the residence time of the evaporation section does not enable adequate evaporation for phase differentiation based on cutoff size with standard optical particle counting. During HyICE-2018, the CFDCs were operated at relative humidities below the point of droplet breakthrough in order to prevent such a situation.

CFDCs are online instruments that measure INP concentration in real time with a minimum time resolution determined by the instrument-specific particle counting method. However, since ambient INP concentrations are generally low, measurements typically consist of multi-minute counting averages. For the CFDCs used in this study, sampling intervals varied between 5 and 20 min, depending on the instrument, and were separated by background measurements of clean, filtered air. Ice crystal concentrations were obtained from the OPC counts and averaged over time. Concentrations obtained during background measurements were then subtracted from the concentrations measured during each sampling window to compute the measured INP concentrations, where the lower limit of detection is defined as 1 standard error above the background mean.

Overall, CFDC measurements are time-limited by the quality of the thin ice layer coating the chamber walls, which deteriorates over time and contributes to increasing particle counts as the instruments operate. Single experiments typically last 3 to 5 h, after which the ice coating needs to be regenerated to restore low background conditions. This is done by warming and purging the chamber, before re-cooling and re-coating the walls with ice. This process can last 1 to 3 h, allowing for two to four daily measurement cycles if continuous operation is desired . In the following sections we describe the instrument specifics for each CFDC run during HyICE-2018.

I. Portable Ice Nucleation Chamber – PINC

PINC is a first-generation, field-deployable, parallel-plate CFDC based on the designs by , and it has been operated in many locations around the world for more than a decade . During HyICE-2018, PINC was operated from 19 February to 2 April (Fig. ) in the main aerosol cottage (Figs.  and ), with one to two experimental cycles conducted per day, always during daytime. For the duration of the campaign, PINC was operated at a fixed lamina temperature of $T_{\mathrm{l}}$ = $-$31 ${}^{\circ }$C and a relative humidity with respect to water ${\mathrm{RH}}_{\mathrm{w}}=\mathrm{105}$ %. These conditions were selected to simulate mixed-phase cloud conditions and correspond to the condensation/immersion freezing mode(s) of ice nucleation . Sampling was performed from a total aerosol inlet mounted outside of the building, 6 m above ground level (Fig. ). The inlet was heated to 25–30 ${}^{\circ }$C to evaporate droplets and ice crystals and had a carrier flow rate of 250 L min${}^{-\mathrm{1}}$. Individual instruments then sampled from manifolds on this inlet using their own external pumps. For PINC, a 4 L min${}^{-\mathrm{1}}$ sample flow was extracted, and a cyclone was used to eliminate particles larger than 2.5 $\mathrm{\mu }$m. A molecular sieve dryer was installed to keep the relative humidity of the sample below 30 %. After the dryer, the sample flow was split into four: 1 L min${}^{-\mathrm{1}}$ to a condensation particle counter (CPC; TSI model 3010), 1 L min${}^{-\mathrm{1}}$ to an aerodynamic particle sizer (APS; TSI model 3321), 1 L min${}^{-\mathrm{1}}$ to a scanning mobility particle sizer (SMPS; with one Hauke-type differential mobility analyzer (DMA) and one TSI model 3772 CPC) and 1 L min${}^{-\mathrm{1}}$ to PINC Fig.  and.

Figure 4

Instrumental setup in the main cottage. Note that PINC and its complete setup (PFPC, CPC, APS and SMPS) were used from 19 February to 2 April 2018 and were replaced by PINCii from 22 April to 10 June 2018. Between May 4 and May 23, the SMPS was used in parallel with PINCii, as represented with the dashed lines (see Sect. 2.2.1 for details).

[Figure omitted. See PDF]

On several occasions, a portable fine particle concentrator (PFPC) was used to concentrate aerosol particles in the sample flow upstream of PINC (Fig. ). The PFPC, described by and based on the design by , is a multistage concentrator based on virtual impaction. It concentrates aerosol particles with a certain size-dependent enrichment factor where larger particles are concentrated more efficiently than smaller ones. The size-dependent enrichment factor is determined by measuring the particle size distributions before and after the PFPC. The enrichment factor was estimated as $\mathrm{25}\pm \mathrm{6}$ for ambient particles of diameters between 0.4 and 2.5 $\mathrm{\mu }$m, when the PFPC was operated at sea level in the vertical configuration . During the HyICE-2018 campaign, a second enrichment factor was determined before each ice nucleation experiment by calculating the ratio between a concentrated INP measurement point and an ambient measurement point bypassing the PFPC. Ambient INP concentrations were then back-calculated by multiplying the concentrated INP concentrations by this second enrichment factor .

When the PFPC was used, the 2.5 $\mathrm{\mu }$m cyclone was removed and replaced by the 2.5 $\mathrm{\mu }$m impactor located inside the PFPC (see Fig. ). On 22 March, the cyclone was removed so that PINC directly sampled the inlet air and the operating conditions were changed to $T_{\mathrm{l}}=-\mathrm{29}$ ${}^{\circ }$C, and RH${}_{\mathrm{w}}=\mathrm{105}$ % for inter-comparison purposes.

In this article we summarize the performance of PINC during the targeted days of instrument inter-comparison, while the longer-term results from the PINC measurements during HyICE-2018 are presented and discussed by .

II. Portable Ice Nucleation Chamber II – PINCii

PINCii is a parallel-plate CFDC developed as an upgrade to the PINC instrument. Although many specific engineering details differ, the primary differences between PINC and PINCii are the chamber dimensions and cooling power. While PINC's main chamber and evaporation section are 568 and 230 mm in height, respectively , PINCii is approximately twice as large with a main chamber of 1000 mm and an evaporation section of 440 mm. A manuscript outlining the engineering and operational details of PINCii is in preparation .

During the campaign, PINCii was operating from 22 April to 10 June 2018 and measured INP concentrations at fixed $T_{\mathrm{l}}=-\mathrm{32}$ ${}^{\circ }$C and ${\mathrm{RH}}_{\mathrm{w}}=\mathrm{105}$%. PINCii was located in the main cottage and sampled from the heated total aerosol inlet, essentially acting as a substitute for the earlier PINC measurements as depicted in Fig. . However, during this later sampling interval, the PFPC was not used, and the additional SMPS was used for the first 2 weeks of May only. For the typical setup during the campaign, PINCii sampled downstream of the 2.5 $\mathrm{\mu }$m cyclone.

III. SPectrometer for Ice Nuclei – SPIN

SPIN is a commercially developed CFDC manufactured by Droplet Measurement Technologies in Boulder, CO, USA, and based on the parallel-plate PINC design. The design and use of SPIN have been previously documented , and some SPIN instruments ($\approx$ 10 exist worldwide) have participated in earlier inter-comparison activities . For HyICE-2018, the University of Eastern Finland SPIN instrument (UEF-SPIN) was installed and operated for several days between 20 March and 3 April. SPIN was situated in a measurement container approximately 200 m from the main cottage (see Fig. ). SPIN sampled through a burled PM${}_{\mathrm{2.5}}$ inlet (DIGITEL Elektronik AG, Switzerland) located 110 cm above the roof of the container. The PM${}_{\mathrm{2.5}}$ inlet had a nominal sampling flow rate of 16.67 L min${}^{-\mathrm{1}}$, regulated by an external pump and mass-flow controllers (MFCs) (Fig. ). A virtual impactor (VI)-type concentrator, which was a scaled-up version of the construction described by , was used to enhance SPIN's sampling and detection of INPs (more details about the VI concentrator are provided below). The configuration illustrated in Fig.  provided an amplification of the particle number concentration by approximately 9-fold. A CPC (Airmodus model A20) was operated in parallel with the SPIN to monitor the concentrated particle number concentration during measurements.

Figure 5

SPIN instrumental setup during the HyICE-2018 campaign.

[Figure omitted. See PDF]

The use of the VI concentrator may have generated biases in the SPIN measurements. Indeed, the magnification factor of the VI concentrator was most prominent for particles with diameters from 1.3 to 2.5 $\mathrm{\mu }$m, with an average of $\mathrm{8.45}\pm \mathrm{0.43}$ (Fig. ). For smaller particles, the magnification factor decreased steeply, being 4.2 for 850 nm particles. This resulted in larger supermicron particles being over-represented in the sampled particles in comparison to their ambient number concentrations. Indeed, the APS data showed very low number concentrations for particles with diameters larger than 1.3 $\mathrm{\mu }$m, with a daily average of 0.26 cm${}^{-\mathrm{3}}$. Because of the over-representation of supermicron particles, and since larger particles with larger surface areas are increasingly likely to host active sites, making them better INPs , it can be expected that the role of such particles as INPs is more pronounced in the SPIN observations. Note that, due to the absence of a rigorous correction method for the SPIN measurements during HyICE-2018, the SPIN data presented here are the concentrated INP concentrations per volume of sampled air and not the back-calculated ambient INP concentrations.

During its first 3 measurement days, SPIN was programmed to measure using automated RH scans, where the relative humidity with respect to water, ${\mathrm{RH}}_{\mathrm{w}}$, was increased from 70 % to 110 % while keeping a constant lamina temperature. Such scans were realized for lamina temperatures of $T_{\mathrm{l}}=-\mathrm{40}$, $-$36 and $-$32 ${}^{\circ }$C. On these 3 measurement days, a total of 33 RH scans were performed with 5 min background checks between successive scans. Thereafter, the sampling program was modified to be more suitable for ambient particle measurements, and a modified scanning protocol was used for the remainder of the campaign. The modified sampling protocol consisted of longer scans at $T_{\mathrm{l}}=-$32 and $-$28 ${}^{\circ }$C, where at each $T_{\mathrm{l}}$ measurements were made with two saturation conditions (${\mathrm{RH}}_{\mathrm{w}}=\mathrm{95}$ % and 110 %). At each of the four conditions, 20 min sampling intervals followed by 5 min of background sampling were used. Measurements continued until the background signal exceeded 10–15 particles L${}^{-\mathrm{1}}$, after which the chamber was re-iced to return to a lower background signal .

For inter-comparison purposes, SPIN sampled at $T_{\mathrm{l}}=-\mathrm{29}$ ${}^{\circ }$C and ${\mathrm{RH}}_{\mathrm{w}}=\mathrm{105}$ % on 22 March and at $T_{\mathrm{l}}=-\mathrm{31}$ ${}^{\circ }$C and ${\mathrm{RH}}_{\mathrm{w}}=\mathrm{105}$ % on 28 March.

2.2.2 Expansion chamber – PINE

The Portable Ice Nucleation Experiment (PINE) chamber (Fig. b), as described in , has been developed based on the working principle of the Aerosol Interactions and Dynamics in the Atmosphere (AIDA) cloud chamber , which simulates cloud formation in rising atmospheric air parcels. The PINE chamber is operated in a cycled mode during which it is first flushed with ambient aerosol to renew the sampled volume of air under investigation. During HyICE-2018, this first mode was run for 4 min with a flow rate of 3 L min${}^{-\mathrm{1}}$. Then, an expansion mode is initiated by sealing the main inlet valve and pumping air out of the chamber, leading to a decreasing pressure and a decreasing (but well-mixed) particle number concentration within the chamber. As the pressure decreases to 700 mbar in approximately 40 s, the volume of gas expands and the gas temperature decreases, leading to an increase in the relative humidity within the chamber. Upon reaching water saturation, aerosol particles within the chamber activate to form supercooled liquid droplets. If aerosol particles immersed within the droplets are also active as INPs at the respective droplet temperature, then those droplets will freeze via the immersion-freezing pathway. Frozen droplets are optically detected at the exit of the chamber using an OPC (Welas model 2500p) and are differentiated from liquid droplets based on their size. After the expansion, a final refill procedure is conducted, wherein the chamber continues to be depressurized using dry, filtered air to avoid any icing. The total experimental time of the three-mode cycle is approximately 6 min, and each cycle generates one measurement point. Longer averaging times (e.g., 1, 6 or 24 h) are commonly implemented in postprocessing for purposes of statistical analysis. PINE is operated in a way that ensures frost-free walls such that a subtraction of any background ice counts is not needed.

During HyICE-2018, PINE was located in the main cottage (Figs.  and ) and was operated continuously from 13 March until 11 May 2018 at measurement temperatures between $-$24 and $-$32 ${}^{\circ }$C. During this time, PINE was operating with either a constant temperature or with a stepwise temperature ramping, during which the temperature was lowered three times by 2–3 ${}^{\circ }$C each hour. In the main cottage, PINE sampled from the heated total aerosol inlet without using a size cutoff in order to sample the total aerosol. Particle transmission efficiency as a function of particle size (Fig. ) was investigated by measuring particle concentrations upstream of the PINE chamber inlet and in the ambient air using an OPC (MetOne model GT 526S). The results illustrate that PINE was effectively sampling PM${}_{\mathrm{5}}$ aerosol, although the cutoff was not as sharp as would be expected from a traditional impactor.

During the inter-comparison days, PINE was operating at a constant temperature close to the lamina temperature, $T_{\mathrm{l}}$, selected for the CFDCs. On 22 March, PINE sampled at a temperature of $T=-\mathrm{29}$ ${}^{\circ }$C, while on 28 March the temperature was lowered to $T=-\mathrm{30}$ ${}^{\circ }$C. For the inter-comparison conducted in April, PINE sampled at $T=-\mathrm{29}$ ${}^{\circ }$C.

While this study focuses on PINE's results during the inter-comparison days, the long-term measurements are presented and discussed in .

2.2.3 Filter sampling for droplet-freezing assays

Numerous droplet-freezing assay techniques exist for offline measurements of INP concentrations in aerosol samples. The techniques presented here relied on filter collection and resuspension of collected material into known volumes of liquid (Fig. c). The liquid is generally separated into sample aliquots which are exposed to decreasing temperatures and monitored for freezing using optical techniques. Aliquots are assumed to freeze when the most “active” INP in any given sample volume initiates ice formation. A series of assumptions allows the concentration of INP at a given freezing temperature to be calculated from the number of frozen aliquots . It is important to note that these techniques typically assess the ice nucleation tendencies of sample material that is immersed within the liquid (Fig. d). During HyICE-2018, two droplet-freezing assay techniques using similar sampling protocols but different freezing systems were deployed, and the analyses were performed in the SMEAR II laboratories.

I. Ice Nucleation SpEctrometer of the Karlsruhe Institute of Technology – INSEKT

The INSEKT is a droplet-freezing assay based on the Colorado State University Ice Spectrometer (CSU-IS) design, which was originally developed by . A detailed description of the INSEKT setup and working principle can be found in . For INSEKT, filters were collected beginning 2 March and continuing through the end of the intensive campaign in mid-June 2018. Further INSEKT filters were collected more intermittently at the site for more than 1 more year (until May 2019), and that data yielded detailed information related to the seasonality of observed INPs . Ambient aerosol particles were collected onto 0.2 $\mathrm{\mu }$m filters (Whatman Nuclepore track-etched polycarbonate membranes, 47 mm), which were precleaned with a 10 % H${}_{\mathrm{2}}$O${}_{\mathrm{2}}$ solution and subsequently rinsed with deionized water that was passed through a 0.1 $\mathrm{\mu }$m syringe filter (Whatman). After collection, the filters were stored in sterile petri dishes, wrapped in aluminum foil and frozen until analysis with INSEKT. Sampling was performed at three different SMEAR II locations. The primary sampling unit was installed at the aerosol cottage (Fig. ) at ground level, where daily 24 h filter samples were collected through a PM${}_{\mathrm{10}}$ inlet mounted on a $\approx \mathrm{1.8}$ m vertical sampling line. A secondary filter sampling unit was used in a container on top of the 35 m tower, again using a PM${}_{\mathrm{10}}$ inlet and $\approx \mathrm{1.8}$ m vertical sampling line. Generally, from the 35 m tower, two filters were collected per day ($\approx$10 hours during the day and $\approx \mathrm{14}$ h during the night). However, there were 10 d during the campaign when the tower sampling time was extended to 24 h in order to directly compare with the results from the aerosol cottage. On 28 March, two filter samples were collected at the main cottage for inter-comparison with the other INP instruments. For both filters, the sampling time was approximately 3 h, and sampling was carried out using a total aerosol inlet. Higher-altitude air was also sampled using an aircraft-mounted filter sampling unit (see Sect.  for more details).

After sampling, the filters were suspended in 8 mL of Nanopure water that had been passed through a 0.1 $\mathrm{\mu }$m syringe filter. The sample solution was then spun on a rotator for approximately 20 min in order to wash the collected aerosol particles from the filter. Small 50 $\mathrm{\mu }$L volumes of the resulting liquid were pipetted into the 192 wells of two polymerase chain reaction (PCR) plates. The plates were then inserted into the freezing apparatus and were cooled with a constant 0.25 or 0.33 ${}^{\circ }$C min${}^{-\mathrm{1}}$ cooling rate using an ethanol chiller (LAUDA Dr. R. Wobser GmbH & Co. KG.; Proline RP 890 in combination with a LAUDA command module). Brightness changes of the small sample volumes, which correspond to freezing events, were detected using a camera (EO-23122, Edmund Optics monochrome camera) and a custom-made LabVIEW program for image acquisition and analysis. The number of frozen volumes, which increases as the temperature decreases, is used to determine a temperature spectrum of INP concentration . The INSEKT is able to measure INP concentrations at temperatures between $-$5 and $-$25 ${}^{\circ }$C, which is relevant for heterogeneous freezing conditions within supercooled mixed-phase clouds.

Handling blank filters, which were collected without ambient air flowing through the membranes, were used during the campaign to verify that contamination was not an issue in the absence of ambient particles. Moreover, several control freezing spectra were done using the Nanopure water to ensure that the background freezing due to impurities in the water remained low relative to the number of INPs in the water after aerosol sampling. The INP concentrations reported here were obtained by subtracting the INP concentrations derived from both the handling blank filters and the freezing background measurements.

After quantifying the INP content of the aerosol samples, some heat treatment tests were performed to investigate the heat sensitivity of the sampled aerosol with respect to their freezing ability, which can be used to investigate if biological particles contributed to the INP population . For the heat treatment tests, a polypropylene test tube (CELLSTAR, Greiner Bio-One) filled with 2 mL of the aerosol suspension was placed in a 100 ${}^{\circ }$C water bath for about 20 min. The heat-treated samples were then reanalyzed with the INSEKT to quantify changes in the INP temperature spectra.

II. Microlitre Nucleation by Immersed Particle Instrument – $\mathrm{\mu }$L-NIPI

The $\mathrm{\mu }$L-NIPI, similar to INSEKT, is also a droplet-freezing assay instrument for offline measurements of INP concentrations, requiring the collection of aerosol samples on filters. During HyICE-2018, filters were collected from 7 March until 10 June 2018, using omnidirectional ambient air particulate samplers (BGI PQ100, Mesa Laboratories Inc.). Filters were mainly collected at ground level next to the aerosol cottage using PM${}_{\mathrm{10}}$, PM${}_{\mathrm{2.5}}$ and PM${}_{\mathrm{1}}$ inlets on 0.4 $\mathrm{\mu }$m pore track-etched membrane polycarbonate filters (Nuclepore, Whatman). Some filters were also collected on top of the 35 m tower using a PM${}_{\mathrm{10}}$ inlet and aboard an aircraft. In addition, filters were collected using a prototype version of the Selective Height Aerosol Research Kit SHARK; to quantify more detailed, size-resolved measurements of the INP population, although this was done only sporadically for a few different sampling periods.

After sampling, collected aerosol particles were washed from filters using 5 mL of Nanopure water that had been filtered through a 0.2 $\mathrm{\mu }$m filter (Sartorius, model Minisart). Droplets of the sample solution containing particles were pipetted onto a hydrophobic glass slide that holds approximately 50 droplets of 1 $\mathrm{\mu }$L (Fig. c). The glass slide was placed on the temperature-controlled plate of the $\mathrm{\mu }$L-NIPI, which was cooled to $-$40 ${}^{\circ }$C at 1 ${}^{\circ }$C min${}^{-\mathrm{1}}$. The freezing temperature of each droplet was recorded via a digital camera using changes in contrast to determine when a droplet had frozen . Note that, as for INSEKT, handling blank and water background measurements were done to determine the baseline of the results, and only the data that are statistically above the baseline are reported.

Heat tests were performed using the same technique as described for the INSEKT, except the samples were heated for 30 min instead of 20 min . The heated samples were analyzed using the $\mathrm{\mu }$L-NIPI, and the results were compared to the original unheated sample results in order to quantify any changes in INP activity, which is used to infer information about INPs of biogenic origin.

2.3 Additional aerosol particle characterization

Beyond INP measurements, further enhanced efforts were made to quantify, classify and assess aerosol properties at SMEAR II during the HyICE-2018 campaign. Multiple additional particle counters were installed to provide physical characterization of particles (Table ). In addition to these ground- and tower-based instruments, airplane, drone flights and remote sensing retrievals were used to provide insight into the vertical distribution of aerosols and INPs. Moreover, the SMEAR II station has a unique boreal location and is known for the documented occurrences of NPF events driven by biogenic volatile organic compounds emitted by the vegetation of the boreal forest . Hence, extra effort was made to assess whether any links between biology and INP emerge, with a special focus on the seasonal transition of the forest biome.

2.3.1 Search for bio-ice nucleators

The boreal forest is a diverse ecosystem in which biological drivers of aerosol properties have been previously identified. For example, past research at the SMEAR II research station has shown that gas-phase BVOCs can act as precursors to NPF . In addition, previous studies have identified several biological influences on INPs, including biological detritus , pollen , bacteria such as Pseudomonas syringae , fungi and other microorganisms such as viruses . Therefore, during HyICE-2018, we aimed to assess potential links between biological activity and INPs.

In order to evaluate INPs of biological origin, a wideband integrated bioaerosol sensor (WIBS; Droplet Measurement Technology model WIBS-NEO) was used for online analysis of particles during the campaign. The WIBS measures particles between 0.5 and 30 $\mathrm{\mu }$m using a light-scattering technique. In addition to particle counting, particles trigger two optically filtered (280 and 370 nm) xenon lamps that excite the fluorescence of biological particles. The emission is monitored in two detection bands (310–400 and 420–650 nm), whereafter an additional fluorescence threshold can be used to distinguish biological particles (strong fluorescence response) from other materials (e.g., some types of dust, soot) that are typically more weakly fluorescent. In this study, a threshold of FT $+$ 9$\mathit{\sigma }$ (where FT is the mean value of the forced trigger intensities and where $\mathit{\sigma }$ is their standard deviation) was used to determine if a particle is fluorescent details in. The use of two excitation wavelengths and two detection channels allows for additional resolution in the fluorescence analysis, because different particle types often have different fluorescence intensities and emission bands . An additional benefit of the WIBS instrument is that it allows particle asphericity to be calculated . During HyICE-2018, the WIBS was operated from 11 March to 25 June at a flow rate of 0.3 L min${}^{-\mathrm{1}}$ with an acquisition rate of 13 Hz, and the data were later averaged into user-defined time intervals of 10 min. The WIBS was first installed in the main cottage (11 March–3 April; Fig. ), where it was sampling from the heated total aerosol inlet described previously (Figs.  and ), and was therefore effectively sampling PM${}_{\mathrm{5}}$ aerosol. Later, the instrument was moved into the aerosol cottage (3 April–25 June) and attached to a PM${}_{\mathrm{10}}$ inlet described in . Thus the WIBS data reported here are for particles between 0.5 and 5 or 10 $\mathrm{\mu }$m, depending on the instrument location.

In addition to online measurements, several offline sampling techniques were used to investigate potential biological contributions to INP from primary biological materials. Plant materials were collected and washed to test for ice-nucleation-active biological materials, and additional air samples were collected to cultivate fungal samples for DNA sequencing. For fungal cultivation, filter samples were collected overnight from 21 to 22 May and 22 to 23 May using a ILFH-27 filter holder (HI-Q Environmental Products Company Inc.) with Teflon filters and a flow rate of 20 L min${}^{-\mathrm{1}}$. After sampling, the filters were suspended in a buffer solution (100 mM NaCl, 0.5 mM CaCl${}_{\mathrm{2}}$, 1.0 mM MgCl${}_{\mathrm{2}}$, 20 mM potassium phosphate pH 7.2, filter sterilized through 0.22 $\mathrm{\mu }$m pore size), and 200 $\mathrm{\mu }$L of the sample was spread on Luria-Bertani plates and incubated at 22 ${}^{\circ }$C for $\approx \mathrm{7}$ d or until visible colonies appeared. Single filamentous colonies of different morphology were isolated and inoculated on fresh media and subsequently incubated at 22 ${}^{\circ }$C. The ice nucleation activity of each filamentous colony was then tested using the $\mathrm{\mu }$L-NIPI, and the results are presented in the manuscript in preparation.

2.3.2 Vertical profiling

The vertical distribution of INPs through the atmospheric mixed layer is a key area of investigation within the INP research community. During HyICE-2018, several techniques were used in order to assess if vertical gradients exist in INP concentrations, which might bias ground-level measurements.

Airborne measurements were made between March and May 2018 using a Cessna 172 aircraft equipped with onboard instrumentation for classifying aerosol physical properties and equipment for INP filter sampling for offline characterization using INSEKT and $\mathrm{\mu }$L-NIPI. In-flight monitoring of aerosol physical properties was conducted using a particle size magnifier (PSM; Airmodus model A10) and a CPC (TSI model 3776) to measure particle number concentrations. Furthermore, a SMPS and an optical particle sizer (OPS; TSI model 3330) were used to classify particle number size distributions. A shrouded solid diffuser inlet, designed based on the University of Hawaii inlet , with a 5.0 $\mathrm{\mu }$m aerodynamic diameter cutoff was used. With this setup, particle number size distributions from 1.5 nm to 5.0 $\mathrm{\mu }$m were measured. Additional sensors (Rotronic HygroClip-S, PT1-100 temperature sensor and LI-COR Li-840) measured relative humidity, temperature, and CO${}_{\mathrm{2}}$ and H${}_{\mathrm{2}}$O concentrations during the flights. The aircraft's GPS receiver recorded latitude, longitude and flight altitude. The boundary layer depth was estimated during the flights using the real-time measurements of particle concentration, relative humidity and temperature, and it ranged between 500 and 2500 m (mean = 1300 m, SD = 704 m). Filter samples were collected within the boundary layer and the free troposphere to characterize the vertical distribution of INP concentrations. Measurement flights took off from the nearby ($\approx \mathrm{60}$ km distant) Tampere–Pirkkala airport, with 3 h flight plans that consisted of 20–40 km segments flown over the SMEAR II research station at different altitudes between 100 and 3500 m above ground level.

Drone flights flown by the Finnish Meteorological Institute's Remotely Piloted Aircraft System (FMI-RPAS) rotorcraft hexacopter were conducted on 27 March in an effort to assess the state of lower boundary layer aerosol and for comparisons with the instruments installed at various heights on the measurement towers. The RPAS carried a payload that included two CPCs (TSI model 3007) with different cutoff diameters (7 and 14 nm), an OPC (Alphasense model N2, 16 size bins from 0.38 to 17 $\mathrm{\mu }$m) to classify particles by size, and sensors to measure temperature, relative humidity (Vaisala HMP110 probe), and pressure (Arduino Bosch BME280). More details on the platform and the instrumentation setup are given in . During HyICE-2018, flights were typically conducted below 700 m above ground level.

The SMEAR II station is additionally an ACTRIS cloud profiling station equipped with a cloud radar (94 GHz FMCW

FMCW represents frequency modulated continuous wave.

Doppler cloud radar, RPG-FMCW-94-DP), a Doppler lidar (DL, HALO Photonics) and a profiling microwave radiometer (Humidity And Temperature PROfiler, RPG-HATPRO). The cloud profiling instruments were operational during the HyICE-2018 campaign. The cloud observations are used to quantify cloud properties such as cloud boundaries, phase, ice water content and liquid water content . Moreover, radar observations can be used to identify ice particle growth processes like riming and the onset of ice particle formation . In combination with INP measurements, these observations are useful for identifying cases that may be attributed to secondary ice production .

3 Results and discussion

The wide range of techniques and instrumentation employed during HyICE-2018 has led to many results. In this paper, we summarize the general results for the duration of the intensive campaign and focus on a few instances where defined efforts were made to intercompare different INP measurement techniques and instrumentation. Longer-term and more detailed studies that have emerged from specific instruments and/or activities have been or will be published independently, with many contributions aimed at the Copernicus special issue “Ice nucleation in the boreal atmosphere” .

3.1 Meteorological conditions and seasonal change

The HyICE-2018 campaign aimed to capture the seasonal transition from winter to summer conditions. During winter, the campaign was characterized by deep snow cover (60 cm on average) and cold temperatures (between $-$17 and 0 ${}^{\circ }$C) (Fig. a). In 2018, the transition from winter to summer was rather abrupt at the measurement site, and as depicted in Fig. a, snow cover went from a near maximum to completely melted away within a couple of weeks in April 2018. This transition coincides with an increase in the fraction of fluorescent particles (Fig. c), which are used as a proxy for biological particles . By mid-May, the forest ecosystem fully transitioned into summer and ambient temperatures reached nearly 30 ${}^{\circ }$C. In fact, May 2018 was declared anomalously warm over all of Finland . The seasonal change is also noticeable from particle number size distribution measurements (Fig. b), with an increase in particle concentrations beginning in April. The seasonal change is less evident in the NPF event frequency, which is also plotted in Fig. b (white diamonds; ). In Fig. , aerosol characteristics (concentration and NPF occurrence) are represented as a function of wind direction, air temperature, and time for heights (8.4 and 67.2 m) below and above the forest canopy. Although trends might be expected due to varying source regions , no clear correlation is observed between the aerosol features, the wind direction or the changing seasons. The same conclusion is drawn for the intermediate heights of 16.8 and 33.6 m (see Figs.  and ). Such observations further motivate the search for other seasonally dependent variables, such as aerosol chemical composition and biological activity, that may influence the ice-nucleating potential of the aerosols .

Figure 6

(a) Snow depth (1 min resolution, white trace) as a function of time, with the shading representing the mean daytime air temperature. (b) Particle number size distribution from the SMEAR II DMPS. NPF events are indicated by white diamonds. (c) daily average of total particle concentration from the SMEAR II APS, total and fluorescent particle concentration from the WIBS and fluorescent fraction (f = WIBS fluorescent / WIBS total concentration) as a function of time. See Sect. 2.3.1 for more information concerning the WIBS. In all three panels, the instrument inter-comparison days are indicated by the black dashed lines.

[Figure omitted. See PDF]

Figure 7

Wind roses representing multiple aerosol features, including relative average particle concentration (calculated from the SMEAR II DMPS and APS concentrations) and NPF events, as a function of air temperature, wind direction and time at 8.4 (a) and 67.2 m (b) above ground level on the mast.

[Figure omitted. See PDF]

3.2 INP measurement inter-comparisons

To maximize the number of instruments that were available for the inter-comparison study, 22 March, 28 March, 26 April and 28 April 2018 were chosen for instrument inter-comparison. During these days, the online chambers PINC, PINCii, SPIN and PINE were operated with thermodynamic conditions close to one another, and an effort was made to maximize the temporal overlap of measurements. Filter sampling for droplet-freezing assay inter-comparison measurements were carried out on 28 March 2018.

3.2.1 Online instrument inter-comparison

A summary of the thermodynamic conditions used with the online chambers during the inter-comparison days is shown in Table . Note that, as mentioned earlier, both PINC and SPIN were using concentrators to increase their signal-to-noise ratio, which might lead to potential bias in their results. For PINC, the ambient INP concentration was back-calculated from the concentrated measurements, while the SPIN data presented here are the concentrated INP concentration per volume of sampled air.

Table 2

Thermodynamic conditions and inlet settings used for the online instrument inter-comparison studies conducted on on 22 March, 28 March, 26 April and 28 April 2018.

Inter-comparison date	Instrument	Lamina temperature $T_l$	Relative humidity	Inlet setting	Concentrator
		or temperature $T$	with respect
			to water ${\mathrm{RH}}_{\mathrm{w}}$
22 March 2018	SPIN	$-$29 ${}^{\circ }$C	105 %	PM${}_{\mathrm{2.5}}$ + dryer	Yes
	PINC	$-$29 ${}^{\circ }$C	105 %	Morning: heated	Morning: No
				Afternoon: PM${}_{\mathrm{2.5}}$	Afternoon: Yes
				+ heated + dryer
	PINE	$-$29 ${}^{\circ }$C	–	Heated	No
28 March 2018	SPIN	$-$31 ${}^{\circ }$C	105 %	PM${}_{\mathrm{2.5}}$ + dryer	Yes
	PINC	$-$31 ${}^{\circ }$C	105 %	PM${}_{\mathrm{2.5}}$ + heated + dryer	Yes
	PINE	$-$30 ${}^{\circ }$C	–	Heated	No
26 April 2018	PINCii	$-$32 ${}^{\circ }$C	105 %	PM${}_{\mathrm{2.5}}$ + heated	No
	PINE	$-$29 ${}^{\circ }$C	–	Heated	No
28 April 2018	PINCii	$-$32 ${}^{\circ }$C	105 %	PM${}_{\mathrm{2.5}}$ + heated	No
	PINE	$-$29 ${}^{\circ }$C	–	Heated	No

Figure  presents the time series of each chamber's measurements during the inter-comparison days. The PINE data are presented as a 5-point moving average to reduce the uncertainty associated with poor counting statistics in the periods when INP concentrations were close to the detection limit (below $\approx$ 5 L${}^{-\mathrm{1}}$ for a single expansion). The error bars represent 20 % uncertainty in absolute INP concentrations cf.. The SPIN data are generated from the difference between 15 min sampling averages and 5 min interpolated background concentrations, with error bars that represent $\pm \mathrm{1}$ standard deviation of the processed signal. The PINC and PINCii data are processed in an analogous manner to the SPIN data but with sampling windows of 20 and 15 min, respectively, and background windows of 10 and 15 min, respectively.

Figure 8

INP concentrations measured between $-$29 to $-$32 ${}^{\circ }$C for the online INP chambers SPIN, PINC and PINE on 22 and 28 March, and for PINE and PINCii on 26 and 28 April 2018. Aerosol number concentration, $N_{\mathrm{p}}$ ($\mathit{>}\mathrm{0.5}$ $\mathrm{\mu }$m), is represented by the black solid line, and the shaded areas represent the three INP parameterizations from , and . All parameterization envelopes represent a $-\mathrm{29}\ge T\ge -\mathrm{32}$ ${}^{\circ }$C cloud temperature window. See Table  for more information concerning the thermodynamic conditions at which the chambers were running. Note that the single SPIN data point that includes a black cross on 22 March 2018 is deemed to be below the level of detection, defined as the average background signal plus 1 standard error of the mean.

[Figure omitted. See PDF]

As seen in Fig. , overall there is good agreement in the INP concentrations from the PINC and PINE chambers. The INP concentrations are within the same order of magnitude and generally follow the same trend. On 22 March between 17:00 and 19:00 (UTC+2) and on 28 March around 11:30 (UTC+2), the trends in the INP concentration are however opposite. This might be due to the fact that PINE was measuring with a much higher frequency than PINC, and it was therefore able to capture changes in INP concentrations that were not visible in the 30 min PINC data.

In comparison to the PINC and the PINE chambers, the SPIN chamber tends to measure INP concentrations that are lower by a factor of $\approx \mathrm{10}$, despite the use of a concentrator. This systematic offset is likely partially due to instrumental differences. In the cases presented here, the PINC data are corrected by a uniform scaling factor of 1.14 determined from well-characterized particle losses and lamina spreading measured in the instrument . Currently, no simple scaling factor is available for the SPIN during HyICE-2018; thus, no correction factor is used for the SPIN results. However, the biases in the ice-activated fraction from the SPIN chamber were discussed in , when the chamber was used in a separate laboratory experiment. The study shows that at approximately $-$31 ${}^{\circ }$C, the activated fraction is biased low by a factor of $\approx \mathrm{3}$ due to lamina spreading and particle losses, which may partly explain why SPIN measures systematically lower INP concentrations than PINE and PINC during HyICE-2018.

Figure  also shows good agreement between the PINE and PINCii chambers, with INP concentrations measured within the same order of magnitude and overall similar trends throughout the days. The main deviation is observed on 28 April between 12:00 and 13:00 (UTC+2), where both chambers show opposite trends in the INP concentration. Although this deviation cannot be explained at this time, it is short-lived and only represents a few data points.

The aerosol number concentration for particles with an aerodynamic diameter between 0.5 and 20 $\mathrm{\mu }$m (noted $N_{\mathrm{p}}$ ($\mathit{>}\mathrm{0.5}$ $\mathrm{\mu }$m)) obtained from the SMEAR II APS is plotted for comparison (Fig. , solid black lines) and suggests a reasonably constant activated fraction throughout the inter-comparison days. On 22 March, the aerosol number concentration gradually decreases between 10:00 and 14:00 (UTC+2), and it is interesting to observe that such variability is reflected in both the PINE and PINC data.

In Fig. , we also present three parameterizations that predict INP concentrations for simple comparison purposes and to give context to our ambient measurements. The parameterization by was developed by combining observations from nine different field studies, while the parameterization by proposes a modified version focusing on fluorescent biological aerosol particles. Both parameterizations use the number concentration of aerosol particles with diameters larger than 0.5 $\mathrm{\mu }$m and the cloud temperature. In the figures presented here, the and parameterizations were calculated using the aerosol number concentration obtained from the SMEAR II APS ($N_{\mathrm{p}}$ ($\mathit{>}\mathrm{0.5}$ $\mathrm{\mu }$m)) and a cloud temperature window between $-\mathrm{29}$ and $-\mathrm{32}$ ${}^{\circ }$C. The additional parameterization presented in Fig.  has been developed from long-term observations conducted in Hyytiälä during and after the HyICE-2018 campaign, and it relies on both ambient air and cloud temperatures.

From the figure, the parameterization shows the best agreement with the INP concentrations measured during the inter-comparison study. The parameterization tends to underestimate the INP concentration measured by all the online chambers except for SPIN, whose uncorrected INP concentration is systematically lower than the concentration measured by the other chambers, as discussed previously. Furthermore, the parameterization tends to overestimate the INP concentration measured and does not reflect the variability observed in the INP concentrations. The parameterization relies on ambient air temperature, which was relatively constant over the inter-comparison days, and was developed to capture the seasonal variation of INP concentration with a time resolution of 1 to several days. Thus the longer time resolution inherent in the parameterization might explain why it fails to reflect the daily variability observed in INP concentrations. Moreover, the parameterization was established using INP data with $-\mathrm{23}$ ${}^{\circ }$C minimum temperatures and was not tested for applications at lower temperatures; thus, deviations between the online measured and the calculated INP concentrations were to be expected. On the other hand, the parameterization uses data from nine different campaigns conducted in various environments, sometimes well away from aerosol sources (e.g., in the Arctic), which might explain why it tends to underestimate the INP concentrations measured in the boreal forest of Hyytiälä. The good performance of the parameterization might be related to the fact that the parameterization focuses on biological aerosol particles, which might represent an important proportion of the aerosol population in Hyytiälä (Fig. c).

Despite the reasonable agreement with the parameterization, on 22 March none of the parameterizations successfully represents the measured concentrations. Before 13:00 (UTC+2) on 22 March, the INP concentrations are higher than the concentrations predicted by and are more comparable to the parameterization. In contrast, INP concentrations measured after 20:00 (UTC+2) on 22 March and after 17:00 (UTC+2) on 28 March decrease to below the parameterization and are better represented by the parameterization. Therefore, although the parameterization performs better than the others, the INP concentrations measured during the 4 inter-comparison days seem to be influenced by factors other than those included in the parameterizations. Such an observation highlights the need for new and improved understanding to better represent the INP concentrations in boreal forest environments.

3.2.2 Offline instrument inter-comparison

Filter sampling for the droplet-freezing assay inter-comparison was conducted on 28 March and extended through the morning of 29 March. Figure a depicts the time periods during which sampling took place. Although the filters were not collected at the exact same time, efforts were made to coordinate the measurements, and two additional filters were collected at the main cottage as a complement to the filters collected on a daily basis. The INP temperature spectra of the collected filters are presented in Fig. b. Although the onset freezing temperatures of the $\mathrm{\mu }$L-NIPI samples are 5 ${}^{\circ }$C lower than for the INSEKT samples, the INP temperature spectra show substantial overlap between the two techniques. Moreover, although the inter-comparison period was relatively short and filters were collected over 24 h only, a reasonable temporal agreement is observed, and both techniques show lower INP concentrations for the filters that were collected during the night compared to the filters collected exclusively during the day. As reflected in Fig. , the aerosol number concentration varies slightly during the day, with a minimum of 500 L${}^{-\mathrm{1}}$ around 16:00 (UTC+2). Although the variations remain within 1 order of magnitude, changes in the aerosol number concentration could explain the differences in the INP concentration measured between the day and night filters.

Figure 9

(a) Timeline of the filter sampling carried out on 28 March 2018. The droplet-freezing assays used for the INP analysis (INSEKT; $\mathrm{\mu }$L-NIPI), the sampling locations (MC = main cottage; AC = aerosol cottage; TW = tower) and the inlets used during sampling (Total, PM${}_{\mathrm{10}}$) are indicated. The colored bars indicate the various sampling time windows and correspond to the colors used in the INP temperature spectra presented in the bottom panel and in Fig. . (b) INP temperature spectra measured for each filter collected on 28 March 2018. For the INSEKT data, the error bars represent the statistical as well as the systematic error of the INSEKT assay. More details related to the calculation of these error bars is given in . For the $\mathrm{\mu }$L-NIPI data, the error bars were calculated using the Poisson Monte Carlo procedure as described in .

[Figure omitted. See PDF]

It should be noted that dilutions were used to extend the INSEKT detection range to lower temperatures. For the INSEKT measurements presented in Fig. , the data points between $-$12 and $-$18 ${}^{\circ }$C are derived by investigation of the non-diluted aerosol suspension, while the data points at lower temperatures correspond to a suspension that was diluted by a factor of 10. The discontinuous drop observed between the series of data points of both the “INSEKT Morning” and the “INSEKT Afternoon” filters occurs at the dilution step and is nonphysical. It might be a consequence of the shorter sampling period used for these two filters (Fig. ) or of inhomogeneity in the suspension caused by particle settling . All data points are presented here as measured, and consistent overlap, certainly within the uncertainty envelopes, is observed for the INP concentrations that are uniformly increasing with decreasing temperature.

In Fig. , the INSEKT and $\mathrm{\mu }$L-NIPI methods are directly compared. To quantify the agreement between the two methods, the concordance correlation coefficient (CCC) is calculated following the method developed by . Reasonable agreement is observed for the full measurement set, with a CCC of 0.81 (Fig. a). The data obtained during the day (Fig. b) show the best agreement (CCC $=$ 0.96), while the data obtained during the night (Fig. d) show less agreement (CCC $=$ 0.53). The main deviation in the agreement between the two methods is shown in Fig. c (CCC $=$ 0.26), which is expected due to a shorter temporal overlap in the sample collection for these two filters. Indeed, the filter for $\mathrm{\mu }$L-NIPI was collected from 12:20 to 17:00 (UTC+2), while the filter for INSEKT was collected from 16:30 to 20:00 (UTC+2), representing only 30 min of overlap between the measurements. Such short temporal overlap, together with the variations in the aerosol number concentration, could explain the deviation in the INP concentrations measured by the two methods.

Figure 10

INP concentration measurements from INSEKT and $\mathrm{\mu }$L-NIPI compared to one another for 28 March. Point color represents the ice nucleation temperature. The panels include (a) the full measurement set, (b) INSEKT day and morning (red and orange bars in Fig. ) versus $\mathrm{\mu }$L-NIPI afternoon measurements (dark green in Fig. ), (c) INSEKT afternoon (light green in Fig. ) versus $\mathrm{\mu }$L-NIPI afternoon measurements (dark green in Fig. ) and (d) INSEKT night (dark blue in Fig. ) versus $\mathrm{\mu }$L-NIPI night measurements (light blue in Fig. ). In each case the dashed line represents the $\mathrm{1}:\mathrm{1}$ line, and the concordance correlation coefficient (CCC) is indicated.

[Figure omitted. See PDF]

3.2.3 Inter-comparison summary

In Fig. , the measurements from all INP instruments sampling on 28 March 2018 are presented in the form of INP temperature spectra. The data points from the online INP measurement systems SPIN, PINC and PINE were obtained by averaging the measured INP concentrations over the entire day, and the error bars now represent the standard deviation of the processed data shown in Fig. . The selected INP parameterizations are also depicted in the figure, where the shaded regions represent the average aerosol number concentration $N_{\mathrm{p}}$ ($\mathit{>}\mathrm{0.5}$ $\mathrm{\mu }$m) $\pm \mathrm{1}$ standard deviation for the and parameterizations, and the average air temperature ${\stackrel{\mathrm{‾}}{T}}_{\text{air}}\pm \mathrm{1}$ standard deviation for the parameterization. For each parameterization, the average and standard deviation were calculated between 08:00 (UTC+2) on 28 March and 08:00 (UTC+2) on 29 March 2018.

As expected, in this context the parameterization performs better at warmer temperatures (between $-$12 and $-$25 ${}^{\circ }$C), temperatures for which the parameterization was established. However, it overestimates the INP concentrations measured at colder temperatures by approximately 1 order of magnitude, and only 19 % of the data points fall within its shaded region (Fig. ). The parameterization, on the other hand, is able to reproduce more of the ambient data, with 35 % of the data points within its shaded region. Moreover, although it fails to predict the lowest concentrations obtained from the offline methods, its trend agrees best with both online and offline INP measurements. Conversely, the parameterization only reproduces 3 % of the data points and does not capture the ice nucleation behavior observed by the droplet-freezing assays at warmer temperatures. However, the parameterization is based on various measurements conducted in very different environments, sometimes well away from aerosol sources; thus, bias can be expected. The performance of the various parameterizations highlights the need to critically assess the validity range(s) of such empirical parameterizations, especially when they are suggested for utilization in predictive models.

Figure 11

INP temperature spectra from all the instruments running on 28 March. Error bars assigned to online measurement techniques represent the standard deviation of the processed data shown in Fig. . The , and parameterizations are replotted with the envelopes now indicative of the average $N_{\mathrm{p}}$ ($\mathit{>}\mathrm{0.5}$ $\mathrm{\mu }$m) $\pm \mathrm{1}$ standard deviation for and , as well as the average air temperature ${\stackrel{\mathrm{‾}}{T}}_{\text{air}}\pm \mathrm{1}$ standard deviation for the Schneider parameterization.

[Figure omitted. See PDF]

4 Conclusions

In this measurement report we introduced the HyICE-2018 campaign conducted at the SMEAR II research station in the boreal forest of Hyytiälä, Finland. The campaign took place from February to June 2018 and utilized the infrastructure of the SMEAR II station with additional instrumentation focusing on INP measurements. The main objectives of the campaign were to

• •

  quantify and characterize INPs in a boreal environment within different thermodynamic conditions,

• •

  examine the seasonal variability of the INP concentrations,

• •

  study the vertical distribution of INPs above the boreal forest.

Several days during the campaign were selected to intercompare INP measurement systems in a field setting. The instrumentation included three CFDCs (PINC, PINCii and SPIN) and one expansion chamber (PINE) for online measurements of INPs and two droplet-freezing assays (INSEKT and $\mathrm{\mu }$L-NIPI) for offline analysis of INPs collected on filters. Results from the online measurements show that, despite the differences in setups and measurement settings, there is good agreement between the PINC and PINE chambers and between the PINE and PINCii chambers. The chambers detect INP concentration in the same order of magnitude and show similar trends throughout the inter-comparison days. The SPIN chamber tends to measure lower concentrations than the other chambers, likely in part due to instrumental biases.

Results from the offline measurements show reasonable agreement between the two droplet-freezing assays INSEKT and $\mathrm{\mu }$L-NIPI, with an overall concordance correlation coefficient of 0.81. A 5 ${}^{\circ }$C deviation in the onset freezing temperature is, however, observed between the instruments, and results from the direct comparison of the obtained INP concentrations illustrate the importance of longer sampling durations and longer temporal overlap to improve agreement.

The measurements from the inter-comparison days were compared to three existing parameterizations. Although the comparison is based on limited data and thus might not be representative of the entire HyICE-2018 campaign, results show that the parameterization tends to underestimate the INP concentrations observed by the online INP chambers for temperatures between $-$29 and $-$32 ${}^{\circ }$C, and it does not capture the ice nucleation behavior observed by the droplet-freezing assays at warmer temperatures between $-$20 and $-$12 ${}^{\circ }$C. The parameterization shows better agreement with the INP temperature spectra obtained from the droplet-freezing assays, but it tends to overestimate the INP concentrations measured by the online chambers between $-$29 and $-$32 ${}^{\circ }$C. The parameterization does not capture the daily variability observed in the INP concentration, but such a result is not surprising since the parameterization was developed to predict seasonal variation of INP concentrations with a time resolution of 1 to several days. Finally, the parameterization by shows the best agreement with the measured INP concentrations from both online and offline INP measurement techniques, which is attributed to the fact that it is focused on biological aerosol particles which might represent an important proportion of the aerosol population in the boreal forest of Hyytiälä. The parameterization was developed from measurements conducted in a mid-latitude ponderosa pine forest ecosystem in Colorado, USA, and it is interesting to see that its results are consistent with some of the measurements obtained during the HyICE-2018 campaign. Although more analyses would be required for a systematic comparison of parameterizations with the INP concentrations measured during HyICE-2018, the results suggest that the sources of INPs in the boreal forest might be comparable to those in a North American forest. On the other hand, despite the good agreement with the parameterization, there is 1 measurement day of inter-comparison when none of the parameterizations successfully represent the measured concentrations. This suggests that, on that day, the INP concentration might be influenced by factors that are not included in the parameterizations used here and that continuous efforts are required to improve parameterizations to better represent the INP concentrations from boreal forest environments.

Although not all the highlighted objectives were addressed in this paper, a number of results from the campaign will be or have already been presented in separate contributions to the Copernicus special issue “Ice nucleation in the boreal atmosphere”. The special issue already includes an examination of the PINC data accumulated during winter 2018, which showed that no persistent local sources of INPs could be identified and postulated that the INPs detected were the result of long-range transport and dilution of INPs with sources far from the measurement site . In another study, present a more than 1-year record of INP concentrations measured with INSEKT that shows a clear seasonal cycle of INP concentrations and INP types, most likely driven by the abundance of boreal biogenic aerosol. Forthcoming studies will explore atmospheric vertical profiles of INPs, INP sources and transport modeling, plausible links between INP abundance and NPF, and the ice nucleation activity of boreal biology such as flora and fungi; they will also utilize the full suite of available measurements and applied machine learning techniques in an effort to illuminate non-obvious INP processes and feedbacks .

Appendix A Table A1

Part 1: overview of the SMEAR II instrumentation running throughout the HyICE-2018 campaign.

Quantity/property measured	Method/instrument	Time resolution	Location
Atmospheric aerosol measurements
Aerosol particle and	Particle size magnifier	1 s	Tower (35 m)
ion size distribution	(PSM, Airmodus model A10) (1–3 nm)		and container (2 m)
	Twin differential mobility particle sizer	10 min	Tower (35 m)
	(DMPS) (3–1000 nm)		and aerosol cottage (8 m)
	Aerodynamic particle sizer	10 min	Aerosol cottage (5 m)
	(APS, TSI model 3320) (0.5–20 $\mathrm{\mu }$m)
	Neutral cluster and air ion spectrometer	1 min	Tower (35 m)
	(NAIS, Airel) (0.8–42 nm for ions)		and aerosol cottage (2 m)
Air ion mobility	Balanced scanning mobility analyzer	10 min	Rea hut (2 m)
distribution	(BSMA) (0.4 7.5 nm)
Cluster composition	Chemical ionization atmospheric	1 s	Container (2 m)
and concentration	pressure interface time-of-flight
	mass spectrometer
	(CI-API-TOF, Tofwerk)
Cloud condensation	CCN counter	1 s	Aerosol cottage (8 m)
nuclei (CCN)	(DMT, model CCN-100)
PM${}_{\mathrm{10}}$ mass	Particulate detection monitor	1 min	Aerosol cottage (5 m)
concentration	(Thermo Scientific, model 5030 SHAR Monitor)
Black carbon mass concentration	Multi-angle absorption photometer	1 min	Aerosol cottage (5 m)
	(MAAP, Thermo Scientific)
	Aethalometer	10 min	Aerosol cottage (5 m)
	(Magee Scientific, model AE-33-7)
Aerosol optical properties	Continuous light absorption photometer (CLAP)		Aerosol cottage (5 m)
	Nephelometer (TSI, model 3563)	10 min	Aerosol cottage (5 m)
	Sun photometer (Cimel, model CE-318)
Aerosol chemical composition	Aerosol chemical speciation monitor	30 min
	(ACSM; Aerodyne Research Inc., USA)
	Organic carbon/elemental carbon analyser	4 h
	(OCEC, Sunset Laboratory)
	Time-of-flight aerosol mass spectrometer
	(LToF-AMS)
Size-classified mass collection	PM${}_{\mathrm{10}}$ impactor	3 d	Aerosol cottage (5 m)
of aerosol particles	(DEKATI)
PM${}_{\mathrm{10}}$ heavy metals	Filter sampler (MCZ, model MicroPNS S7)	7 d
PM${}_{\mathrm{10}}$ main ions	Filter sampler (3-stage EMEP)	7 d

Table A2

Part 2: overview of the SMEAR II instrumentation running throughout the HyICE-2018 campaign.

Quantity/property measured	Method/instrument	Time resolution	Location
Meteorology
Ambient air temperature	Ventilated and shielded sensors (Pt-100)	10 min	Rea hut (2 m)
	Rotronic RS12T/RS24T, HygroMet		Mast (125 m)
	MP102H-530300, HygroClip2 (HC2-S3)
Ambient air pressure	Druck DPI 260 barometer	1 min	Mast (0 m)
	Vaisala PTB210 barometer	10 s	Mast (30 m)
	Metallic aneroid capsule capacitance, Vaisala PA-11	1 min	Field (2 m)
Relative humidity and	Vaisala dew point transmitter	10 s	Field (1.5 m)
dew point temperature	Hygroclip2 probe		Tower (35 m)
	Rotronic humidity and temperature sensor MP102H	5 s	Tower (35 m)
	Chilled mirror dew point monitor	1 min	Tower (16 m)
			and mast (23 m)
	EdgeTech Model 200M meteorological system	1 min	Tower (35 m)
Wind speed and direction	Thies ultrasonic anemometer 2D	1 min	Mast (8.4, 16.8, 33.6
			and 67.2 m)
	Metek uSonic-3 anemometer	0.1 s	Mast (67.2 and 125 m)
	Gill 3D sonic anemometer	0.1 s	Field (1.5 m)
	Vaisala 2D sonic anemometer	1 s	Field (3 m)
Precipitation	Vector ARG-100 tipping-bucket rain gauge	1 s	Tower (35 m)
			and catchment (2 m)
	Bucket weighing, OTT Pluvio2, 200 and 400 cm${}^{\mathrm{2}}$	1 min	Field (1.5 and 4 m)
Surface precipitation and	Optical disdrometer, OTT Parsivel2	1 min	Field (4 m)
particle size distribution	2D video disdrometer (particle video imaging)	1 min	Field (4 m)
	Particle video imaging, PIP, NASA	1 min	Field (1 m)
	Particle holographic imaging, SAKU III	1 min	Field (2 m)
Snow depth	Manual measurements at seven locations	7 d	Catchment
	Jenoptik SHM30 optical snow depth sensor	1 min	Field

Table A3

Part 3: overview of the SMEAR II instrumentation running throughout the HyICE-2018 campaign.

Quantity/property measured	Method/instrument	Time resolution	Location
Remote sensing
Vertical profile of hydrometeors and	94 GHz FMCW Doppler cloud radar(RPG-FMCW-94-DP)	3 s	Field
wind profile; cloud and boundary layer height	Microwave radiometer – Humidity And Temperature PROfiler (RPG-HATPRO)	3 s	Field
	HALO Photonics Doppler lidar Ceilometer	16 s 16 s	Roof Aerosol cottage
Solar and terrestrial radiation
Gas measurements
(includes CO${}_{\mathrm{2}}$, H${}_{\mathrm{2}}$O, CH${}_{\mathrm{4}}$, CO, O${}_{\mathrm{3}}$, SO${}_{\mathrm{2}}$, NO${}_x$ and VOC measurements)
Flux measurements
(includes CO${}_{\mathrm{2}}$, H${}_{\mathrm{2}}$O, CH${}_{\mathrm{4}}$, aerosol particles, momentum and heat flux measurements)
Soil and water balance
(includes snow and water collection; soil temperature, water content, matric potential, and heat flux
measurements; measurements of CO${}_{\mathrm{2}}$, CH${}_{\mathrm{4}}$, and VOC concentrations; and fluxes in the forest floor)
Forest ecophysiology and productivity
(includes exchange of gas by shoots and stems measurements, sap flow measurement, litter collection, etc.)

Figure A1

Measured magnification factor for the virtual impactor (VI)-type concentrator used with SPIN during the HyICE-2018 campaign. Flow rates of 20 and 2 L min${}^{-\mathrm{1}}$ were used for the sheath flow and sample flow, respectively.

[Figure omitted. See PDF]

Figure A2

Transmission efficiency as a function of particle size for the sampling system upstream of PINE. The measurements were made with an OPC (MetOne, GT 526S) with an accuracy of $\pm \mathrm{10}\mathit{\%}$.

[Figure omitted. See PDF]

Figure A3

Wind rose illustrating multiple aerosol features, including relative average particle concentration (calculated from the SMEAR II DMPS and APS concentrations) and NPF events, as a function of air temperature, wind direction and time at 16.8 m above ground level on the mast.

[Figure omitted. See PDF]

Figure A4

Wind rose illustrating multiple aerosol features, including relative average particle concentration (calculated from the SMEAR II DMPS and APS concentrations) and NPF events, as a function of air temperature, wind direction and time at 33.6 m above ground level on the mast.

[Figure omitted. See PDF]

Data availability

The aerosol, trace gas and meteorological data are available at the SmartSMEAR data repository https://avaa.tdata.fi/web/smart;. Contact with the original data contributors can be requested from [email protected]. The INP data presented in this study are available at 10.5281/zenodo.5141574 . Other data are available upon request from the corresponding authors. The data discussed and presented from PINC are also included in and are available at 10.3929/ethz-b-000397022.

Author contributions

OM, MK, TP, JD and DM initiated and planned the HyICE-2018 campaign. ZB and JD largely coordinated and oversaw the campaign on the ground with the help of the permanent SMEAR II staff. ZB, DC, MPA, SDvD, PH, KK, JLa, MiP, JuS, FV, YW, NSA, BB, DB, MID, ADH, PH, KH, LL, ML, AM, MaP, GCEP, PP, UP, TS and NSU conducted measurements during the HyICE-2018 campaign. ZB, DC, MPA, PH, KK, MiP, JuS and (JáS). FV and YW performed the data analysis of their respective instruments and took part in interpretation of the results. JPDA and EG provided the PFPC, and EG set up the instrument at the station. EST, JPDA, DHB, ZAK, JK, AV, MK, BJM, TP, OM and JD participated as supervisors of the HyICE-2018 campaign. ZB and EST wrote the article and with the help of DC constructed the figures. MB helped with writing and revising the article. All authors reviewed the article and contributed to the review response.

Competing interests

The contact author has declared that neither they nor their co-authors have any competing interests.

Disclaimer

Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Special issue statement

This article is part of the special issue “Ice nucleation in the boreal atmosphere”. It is not associated with a conference.

Acknowledgements

The authors would like to gratefully acknowledge Janne Levula, Matti Loponen, Heikki Laakso, Turo Salminen and the rest of the technical staff of the Hyytiälä Forestry Field Station for their expertise, hard work and willingness to help throughout the HyICE-2018 campaign. Erkki Järvinen and the pilots at Air Spark Oy are acknowledged for operating the research airplane. Simo Hakala is acknowledged for his help with the new particle formation event classification. Hanna Danielsson is acknowledged for her significant help with creating Fig. 1.

Financial support

This project received funding from the European Union's Horizon 2020 research and innovation program under grand agreement nos. 654109 and 739530 and TransNational Access via ACTRIS-2 HyICE-2018 TNA project. The work of the University of Helsinki was supported by the Academy of Finland Centre of Excellence in Atmospheric Science (grant no. 307331) and NANOBIOMASS (307537), ACTRIS-Finland (328616), ACTRISCF (329274) and Arctic Community Resilience to Boreal Environmental change: Assessing risks from fire and disease (ACRoBEAR, 334792) Belmont Forum project. In addition, the work of the University of Helsinki was financially supported by the European Commission through ACTRIS2 (654109) and ACTRIS-IMP (871115) and ACTRIS2 TransNational Access and through integrative and Comprehensive Understanding on Polar Environments (iCUPE, 689443), ERA-NET-Cofund and by the University of Helsinki (ACTRIS-HY). The work of the KIT Institute for Meteorology and Climate Research (IMK-AAF) was supported through the Research Program “Atmosphere and Climate (ATMO)” of the Helmholtz Association and by the KIT Technology Transfer Project PINE (N059). EST and DC have been supported by the Swedish Research Councils, VR (2013-05153, 2020-03497) and FORMAS (2017-00564), and the Swedish Strategic Research Area MERGE. BJM and MPA acknowledge the European Research Council, ERC, MarineIce 648661, for funding. PH and JK acknowledge the funding from the Arctic Academy program “ARKTIKO” of the Academy of Finland under grant no. 286558, and PH acknowledges support from the Maj and Tor Nessling Foundation. NSU acknowledges the support of the Alexander von Humboldt Foundation, Germany (1188375). MP, SDvD and ZAK acknowledge the funding from the European Union's Horizon 2020 research and innovation program (under the Marie Skłodowska-Curie grant agreement no. 751470, “ATM-METFIN” and grant agreement no. 654109) and the European Union Seventh Framework Programme (FP7/2007–2013; grant agreement no. 262254). EG and JPDA acknowledge the support from the NSERC (grant no. RGPIN-2017-05972), and EG acknowledges support from the University of Toronto Centre for Global Change Science Graduate Student Research Award. NSA and DHB acknowledge the funding from the Academy of Finland under the Postdoctoral Researcher Grant 309570 to NSA.

Open-access funding was provided by the Helsinki University Library.

Review statement

This paper was edited by Yafang Cheng and reviewed by two anonymous referees.