Birte
Rörup
a,
Xu-Cheng
He
*ab,
Jiali
Shen
ac,
Rima
Baalbaki
a,
Lubna
Dada
ad,
Mikko
Sipilä
a,
Jasper
Kirkby
ef,
Markku
Kulmala
ag,
Antonio
Amorim
h,
Andrea
Baccarini
i,
David M.
Bell
d,
Lucía
Caudillo-Plath
f,
Jonathan
Duplissy
ac,
Henning
Finkenzeller
aw,
Andreas
Kürten
f,
Houssni
Lamkaddam
d,
Chuan Ping
Lee
d,
Vladimir
Makhmutov
jk,
Hanna E.
Manninen
e,
Guillaume
Marie
f,
Ruby
Marten
d,
Bernhard
Mentler
ab,
Antti
Onnela
e,
Maxim
Philippov
j,
Carolin Wiebke
Scholz
l,
Mario
Simon
f,
Dominik
Stolzenburg
amn,
Yee Jun
Tham
ao,
António
Tomé
p,
Andrea C.
Wagner
fq,
Mingyi
Wang
r,
Dongyu
Wang
d,
Yonghong
Wang
as,
Stefan K.
Weber
ef,
Marcel
Zauner-Wieczorek
f,
Urs
Baltensperger
d,
Joachim
Curtius
f,
Neil M.
Donahue
t,
Imad
El Haddad
d,
Richard C.
Flagan
u,
Armin
Hansel
l,
Ottmar
Möhler
v,
Tuukka
Petäjä
a,
Rainer
Volkamer
w,
Douglas
Worsnop
a and
Katrianne
Lehtipalo
ax
aInstitute for Atmospheric and Earth System Research/Physics, Faculty of Science, University of Helsinki, Helsinki, Finland. E-mail: xh346@cam.ac.uk
bYusuf Hamied Department of Chemistry, University of Cambridge, Cambridge, UK
cHelsinki Institute of Physics, University of Helsinki, Helsinki, Finland
dLaboratory of Atmospheric Chemistry, Paul Scherrer Institute, Villigen, Switzerland
eCERN, European Organisation for Nuclear Research, Geneva, Switzerland
fInstitute for Atmospheric and Environmental Sciences, Goethe University Frankfurt am Main, Frankfurt am Main, Germany
gJoint International Research Laboratory of Atmospheric and Earth System Sciences, School of Atmospheric Sciences, Nanjing University, Nanjing, China
hCENTRA, FCUL, University of Lisbon, Lisbon, Portugal
iLaboratory of Atmospheric Processes and their Impacts, École polytechnique fédérale de Lausanne, Lausanne, Switzerland
jLebedev Physical Institute, Russian Academy of Sciences, Moscow, Russia
kMoscow Institute of Physics and Technology, National Research University, Moscow, Russia
lInstitute for Ion and Applied Physics, University of Innsbruck, Innsbruck, Austria
mInstitute for Materials Chemistry, TU Wien, Vienna, Austria
nFaculty of Physics, University of Vienna, Vienna, Austria
oSchool of Marine Sciences, Sun Yat-sen University, Zhuhai, China
pIDL-UBI, Universidade da Beira Interior, Covilhã, Portugal
qAerosol Physics, Tampere University, Tampere, Finland
rDepartment of the Geophysical Sciences, University of Chicago, Chicago, USA
sState Key Joint Laboratory of Environment Simulation and Pollution Control, Research Center for Eco-Environmental Sciences, Chinese Academy of Sciences, Beijing, China
tCenter for Atmospheric Particle Studies, Carnegie Mellon University, Pittsburgh, USA
uDivision of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, USA
vInstitute of Meteorology and Climate Research, Karlsruhe Institute of Technology, Karlsruhe, Germany
wDepartment of Chemistry & CIRES, University of Colorado Boulder, Boulder, USA
xFinnish Meteorological Institute, Helsinki, Finland
First published on 21st March 2024
Iodine oxoacids are recognised for their significant contribution to the formation of new particles in marine and polar atmospheres. Nevertheless, to incorporate the iodine oxoacid nucleation mechanism into global simulations, it is essential to comprehend how this mechanism varies under various atmospheric conditions. In this study, we combined measurements from the CLOUD (Cosmic Leaving OUtdoor Droplets) chamber at CERN and simulations with a kinetic model to investigate the impact of temperature, ionisation, and humidity on iodine oxoacid nucleation. Our findings reveal that ion-induced particle formation rates remain largely unaffected by changes in temperature. However, neutral particle formation rates experience a significant increase when the temperature drops from +10 °C to −10 °C. Running the kinetic model with varying ionisation rates demonstrates that the particle formation rate only increases with a higher ionisation rate when the iodic acid concentration exceeds 1.5 × 107 cm−3, a concentration rarely reached in pristine marine atmospheres. Consequently, our simulations suggest that, despite higher ionisation rates, the charged cluster nucleation pathway of iodic acid is unlikely to be enhanced in the upper troposphere by higher ionisation rates. Instead, the neutral nucleation channel is likely to be the dominant channel in that region. Notably, the iodine oxoacid nucleation mechanism remains unaffected by changes in relative humidity from 2% to 80%. However, under unrealistically dry conditions (below 0.008% RH at +10 °C), iodine oxides (I2O4 and I2O5) significantly enhance formation rates. Therefore, we conclude that iodine oxoacid nucleation is the dominant nucleation mechanism for iodine nucleation in the marine and polar boundary layer atmosphere.
Environmental significanceUnderstanding atmospheric new particle formation is of paramount importance, since it is impacting the climate. Particle formation through iodine oxoacid nucleation has rarely been incorporated into global simulations, despite observations at various coastal locations and increasing concentration of iodine, primarily driven by climate change. This study employs chamber measurements and a kinetic model to explore the influence of temperature, ionisation rate, and humidity on iodine oxoacid nucleation. Our findings indicate that ion-induced formation rates of iodine oxoacid are not significantly affected by temperature but the neutral formation rates increase with decreasing temperature. Furthermore, increases in ionisation rates decrease the ion-induced nucleation rate under atmospherically relevant iodic acid concentrations, and humidity has no significant effect on nucleation. |
New particle formation has been observed all over the world3 at locations ranging from cities4–7 to pristine regions like the Arctic,8 Antarctica,9 remote marine atmosphere,10 and the boreal forest.11,12 Sulphuric acid and organic vapours are believed to be the major contributors to new particle formation in the continental boundary layer and have been extensively researched.13–15
Iodine species have recently been proposed to be important for particle formation processes as well.16–21 This has been observed in coastal17–20,22 and polar environments.23–25 Iodine precursors, emitted from the sea surface,26 sea ice,27 and macroalgae18 are rapidly photolysed in the atmosphere to produce a series of iodine oxides (e.g., iodine monoxide, IO and iodine dioxide, OIO) and oxoacids (e.g., iodic acid, HIO3 and iodous acid, HIO2). The exact iodine particle formation mechanisms have been under debate for more than 20 years, with iodine oxides such as OIO, iodine tetroxide (I2O4) and iodine pentoxide (I2O5) initially thought to be the nucleating molecules.17,18,28–31
However, recent field observations at the Mace Head observatory have revealed the important role of HIO3 in iodine particle formation processes.19 Further experiments carried out at the Cosmics Leaving OUtdoor Droplets (CLOUD) chamber at CERN found that HIO3 is produced from oxidised iodine in the presence of ozone and water vapour,16 and iodooxy hypoiodite (I2O2) is the critical precursor for producing HIO3.32 This simple and efficient production pathway of HIO3 may help explain the ubiquitous presence of HIO3 globally,16 even including inland polluted cities such as Beijing and Nanjing.33
The CLOUD experiment has also revealed that HIO3 alone is not able to explain iodine particle formation and HIO2 is needed to form the initial clusters together with HIO3 (termed the iodine oxoacid nucleation mechanism).16 Recent theoretical work utilising quantum chemical calculations and kinetic simulations further uncovered the critical role of HIO2 in stabilising HIO3 clusters because HIO2 exhibits strong alkaline-like behaviour.34,35 The fact that iodine oxidation produces both HIO3 and HIO2 makes its nucleation process autonomous; this differentiates iodine species from the sulphuric acid nucleation, which typically requires additional base molecules such as ammonia and amines to form stable clusters.15,36 Importantly, the particle formation rates of iodine oxoacids (HIO3 and HIO2) exceed those of sulphuric acid at comparable acid concentrations, especially in the remote areas under low base concentrations.16 The stabilising effect of HIO2 was also seen with sulphuric acid.21,37,38
Since iodine oxoacid nucleation is a recently discovered mechanism, there exists a critical knowledge gap regarding the factors that influence nucleation rates in this chemical system. These factors include atmospheric ionisation rates, humidity, and temperature. The importance of atmospheric ionisation rates lies in their direct influence on particle formation dynamics, particularly in the upper troposphere where charged clusters play a significant role. Similarly, temperature is a crucial parameter affecting formation rates, as variations in temperature can significantly impact the stability and reactivity of precursor molecules.39 Hence, understanding the precise influence of these factors is paramount for accurately representing them in models.
In this study, we carry out particle formation experiments at the CLOUD chamber at CERN aiming to improve the understanding of iodine oxoacid particle formation and its dependency on temperature, humidity and ionisation rate. As water vapour is a key component in the formation of iodic acid, varying the humidity may alter the ratio between iodine oxides and oxoacids. This may further allow us to disentangle the contribution of iodine oxoacids and oxides to iodine particle formation in the marine boundary layer. Additionally, we use a kinetic model to calculate the formation rate of charged iodic acid clusters at different temperatures, ionisation rates and sink. Recent measurements of iodine in the upper troposphere suggest the presence of HIO3 in this region.40 Combining those parameters allows us to evaluate how significant the charged cluster formation mechanism of iodic acid is in the upper troposphere.
We used three different ionisation settings to investigate the effect of ions on nucleation. First, the chamber was in the so-called neutral mode, in which a high-voltage field was used to remove all ions. Two electrodes created an electric field of up to 30 kV m−1, thereby effectively removing all the ions inside the CLOUD chamber. When the electric field was turned off, ions were created in the chamber by natural galactic cosmic rays (GCRs). The ion production rate under such GCR conditions inside the chamber was determined to be 4.1 ion pairs cm−3 s−1.39 In the last setting, the chamber was exposed to a 3.5 GeV/c pion beam from the CERN Proton Synchrotron. This increased the ion production rate by approximately a factor of twenty compared to the GCR conditions.39
In this work, we focus on the nucleation experiments initiated by molecular iodine (I2). During the experiments, stable concentrations of I2 and ozone were maintained inside the chamber. We initiated each nucleation event by turning on the green light (528 nm), which photolysed I2 to produce iodine atoms. Iodine atoms undergo oxidation processes which eventually form HIO3 and HIO2. Each nucleation experiment was continued until a steady state for each ionisation setting was reached, after which the chamber was cleaned and prepared for the next experiment.16
The experiments relevant to this study were conducted at both −10 °C and +10 °C. An additional series of experiments was carried out at +10 °C, with relative humidity (RH) ranging from below 0.008% to 80%. The cryogenic synthetic air employed in this investigation comprises 21% O2 and 79% N2, with a certified water concentration well below 1 ppmv (0.008% RH). To regulate humidity in the chamber, the primary airflow is humidified using a temperature-controlled Nafion humidifier, which utilises ultrapure water. In order to achieve the driest conditions possible for the experiments without added water vapour, the chamber was purged with cryogenic air for over a week before the dry experiments.
Gas-phase species were measured using a nitrate chemical ionisation mass spectrometer (nitrate-CIMS) and a bromide-CIMS. The nitrate-CIMS is an APi-TOF (Atmospheric Pressure interface Time-Of-Flight) coupled with a chemical ionisation unit using nitric acid as the reagent gas. It is used extensively for detecting sulphuric acid, highly oxygenated organic molecules, and iodic acid (HIO3). In this study, it was used to analyse the concentrations of HIO3, HIO2, I2O4 and I2O5. The details of the nitrate-CIMS instrument in this study can be found in Kürten et al. (2014).45 The nitrate-CIMS has an ion filter integrated into its sampling line to separate ions and charged clusters from the neutral molecules; thus, it measures only neutral molecules and clusters in CLOUD. In this study, a bromide-CIMS was used to measure the water content in the chamber.46 It is an APi-TOF coupled with a chemical ionisation unit using dibromomethane (CH2Br2) as the reagent gas. The CH2Br2 is fed into the sheath flow of the inlet, under the illumination of a soft X-ray source, producing bromide anions (Br−). The Br− ions are directed into the sample flow by a negative electric field, which then cluster with neutral molecules in the sample air.46,47 With the measured water content using the H2OBr−,46 we estimated the relative humidity below 2% in the chamber during the dry experiments. The detection limit in this case is inadequate for measuring the water content in the dry cryogenic air. Nonetheless, it serves the purpose of ensuring that no unintended humidity is introduced into the cryogenic air. Above 2% the relative humidity was measured using the dew point monitor DewMaster from EdgeTech. We chose a time when both the dew point monitor and the bromide-CIMS were working well to calculate a calibration factor, which we then used to estimate the humidity below 2%.
We calculated the formation rate at 2.5 nm (J2.5) for five event days, where the HIO3 concentration was between 1 × 106 and 8 × 106 cm−3 and the temperature between 0 and −10 °C. The sulphuric acid and methanesulphonic acid concentrations were negligible (below 7 × 105 cm−3 and 4 × 105 cm−3 respectively) during the events. The standard deviation was calculated for the J2.5 and the HIO3 concentration to estimate the uncertainty.
(1) |
A more detailed description of the formation rate calculation can be found in Dada et al. (2020).48 The formation rate of each experiment was determined as the median value after reaching a steady state.
For field measurements, a different size range had to be used due to limited instrumentation. The formation rates were calculated for a specific size range (from 2.5 to 25 nm in this study), instead of the total particle population above a certain size, because it is possible that other particle sources could contribute to the particle population, especially at larger sizes. In addition, while wall and dilution losses are not applicable for the atmosphere, freshly formed particles can grow out of the size range (i.e., 25 nm), therefore the loss due to growth out of the size range needs to be accounted for (SGR). Finally, the particle formation from the Arctic field study can be calculated using the following equation:49
(2) |
The formation rates were averaged from the beginning of the new particle formation event until its peak. The growth rates (GR5–20 nm) were calculated by Baccarini et al. (2020)24 for the mean diameter mode of the events. Both for the chamber experiments and the field observations, the standard deviation was calculated to represent the errors.
Jtot = JN + J± + Jrec | (3) |
At CLOUD we can directly measure JN, JGCR, and JBeam by varying the ionisation setting of the chamber between Neutral, GCR and Beam, as explained in Section 2.1. During the Neutral stage, we can directly measure JN and during the GCR and Beam stages, we can measure Jtot at different ionisation rates. Therefore, Jiin can be calculated with
Jiin = JGCR − JN | (4) |
This requires that all the experimental conditions except the ionisation rate stay the same. A time series of the formation rate during a CLOUD experiment is plotted in Fig. 1 to visualise the terms contributing to nucleation at different stages of the experiment. In order to disentangle J± from Jiin, the PANDA520 model16,52 was used to simulate the formation rate of charged clusters J± at 1.7 nm. We note that J± differs from Jiin because the model does not factor in neutralised particles resulting from ion–ion recombination.
Fig. 1 An example of the formation rate time series at different ionisation settings to visualise the terminology used in this study. In the first stage of the experiment, the HV clearing field was turned on to establish ion-free conditions in the chamber and thus we measured the neutral formation rates (JN). At the end of this stage (15:50, indicated with a vertical dashed line), the HV clearing field was switched off. After this point, the ions produced by galactic cosmic rays (GCRs) were no longer removed from the chamber and the particle formation rates increased (JGCR). The difference between JGCR and JN represents the ion-induced formation rate (Jiin).50 |
It was previously shown that ion-induced iodic acid nucleation only occurs in the negative ion channel, while the positive ion channel shows no cluster growth.16 Therefore, positive clusters are treated as one of the sinks for the negative clusters in the model. Additionally, evaporation rates are not taken into account in the model since He et al. (2021)16 has concluded that ion-induced nucleation from iodine oxoacids proceeds at the collision limit with negligible cluster evaporation, specifically at temperatures below 10 °C. Furthermore, the model includes parametrisations for wall and dilution loss, ion–ion recombination and ion–neutral collision.
The wall loss is calculated based on the equation from Crump et al. (1981)53 and takes the temperature, pressure, and particle diameter into account.52 The dilution loss is 1.6 × 10−4 s−1, which is calculated from the volume and the total flow of the CLOUD chamber; 26.1 m3 and 250 litre per minute respectively. In the simulations for atmospheric conditions, the wall loss is used to mimic the sink of charged clusters to large particles in the atmosphere (coagulation loss) and wet and dry deposition.
The ion–ion recombination rate (krec) was calculated according to the parametrisation from Israël (1957)54 which takes the pressure and temperature into account. According to Zauner-Wieczorek et al. (2022),55 this parametrisation has the best agreement with field and model data for an altitude range of 0–22 km. We calculated the ion–ion recombination rate for the upper troposphere by adjusting the temperature and pressure according to the international standard atmosphere,56 which is a standardised model of the vertical temperature and pressure changes.
The ion–neutral collision rates and negative ion transfer rates were calculated with the tool provided by Kummerlöwe et al. (2005).57
The ionisation rate (IR) for most of the modelled results in this simulation was set to 4.1 ion pairs cm−3 s−1, based on steady-state concentrations of ions in the CLOUD chamber. Additionally, we used ionisation rates of 10, 20 and 50 ion pairs cm−3 s−1 to investigate the influence of the ion production rate on the formation rate of charged clusters. We chose 50 ion pairs cm−3 s−1 as the maximum value because it represents the ionisation rate in the atmosphere at around 15 km.58
The formation rate of the charged clusters can be calculated as shown here:
J± = N−crit−1[HIO3]kcol | (5) |
[HIO3] is the iodic acid concentration and kcol is the collision rate. N−crit is the concentration of the negative ion cluster concentration at the selected size, which is 1.7 nm in this study. We need to know how many monomers are needed for a cluster to reach the selected size. Under the assumption that the monomers are spherical, we can calculate this number with eqn (6):
(6) |
When comparing panels A and B, it is apparent that the formation rate increases significantly when the temperature decreases. The average HIO3 and HIO2 concentrations are 2 × 107 cm−3 and 5 × 105 cm−3 at +10 °C (A2) and 8 × 106 cm−3 and 2 × 105 cm−3 at −10 °C (B2). Even with the lower concentration of iodine oxoacids, the neutral particle formation rate at −10 °C exceeds that at +10 °C by approximately two orders of magnitude. This confirms that a lower temperature significantly enhances the neutral iodine oxoacid nucleation, consistent with previous studies.16,62,63 The temperature enhancement can be seen as well in the particle size distributions (A3 and B3); the latter clearly features a stronger formation rate.
Furthermore, the ion enhancement of the formation rate is clearly visible at +10 °C (A1), whereas at −10 °C there is no significant change in the total formation rate when the ionisation rate is changed (B1). This is likely due to the fact that ion-induced nucleation of iodic oxoacids is already proceeding at the kinetic limit at +10 °C and therefore cannot increase further when the temperature decreases.16 On the other hand, the neutral iodine oxoacid nucleation is relatively weak at +10 °C and it is significantly enhanced at −10 °C. Hence, the neutral iodine oxoacid nucleation dominates the total formation rate at −10 °C and the ion-induced channel becomes less significant. In contrast, elevating the ionisation rate from 0 to 4.1 ion pairs cm−3 s−1 at +10 °C results in an increase of the formation rate from 0.04 to 3.5 cm−3 s−1, highlighting the importance of ion-induced nucleation at warmer temperatures. Interestingly, increasing the ionisation rate further from GCR conditions to Beam-enhanced conditions increases the formation rate only by a factor of 2 from 3.5 to 7.2 cm−3 s−1. This differentiates the iodine oxoacid nucleation from, e.g., the sulphuric acid–ammonia system in which a larger beam enhancement was observed.15 As ion-induced nucleation is thought to have potential climate impacts especially in the upper troposphere–lower stratosphere,15,64 resolving the ionisation rate enhancement is crucial. We will investigate this phenomenon in the following section.
The model shows that J± only has a very weak temperature dependence, being slightly lower at lower temperatures (Fig. 3A), which is due to the fact that the recombination loss increases with decreasing temperature. Therefore, more ions will be lost due to ion–ion recombination processes. This confirms the conclusion that the temperature dependence of Jtot, which was seen in Fig. 2, comes from the neutral nucleation pathway, which is strongly enhanced by a lower temperature.
In contrast, the ionisation rate has a clear effect on J± (Fig. 3B). At high iodic acid concentrations, a higher ionisation rate leads to higher formation rates. But interestingly, when the iodic acid concentration is lower than about 1.5 × 107 cm−3, increased ionisation rate suppresses J±. This goes against the general conception that an increased ionisation rate always leads to an enhanced formation rate.15,36 This is because higher ionisation rates result in elevated ion concentrations. However, the higher ion concentrations also lead to increased recombination losses between ions of opposite charge. The higher ionisation rate enhances the total ion concentration essentially linearly. On the other hand, the survival probability of an ion monomer (i.e., IO3−) growing to 1.7 nm relies on the growth by sequential HIO3 additions and the sink by neutralisation with cations. It should be noted that the neutralisation process affects every step of the cluster growth process. Hence, the heightened total ion concentrations and the diminished cluster survival probability exhibit distinct dependencies on the ionisation rate. When a certain level of iodic acid is reached, the gain of 1.7 nm ions due to condensational growth outweighs the loss of ions due to recombination and J± increases with the increasing ionisation rate. However, at low to moderate iodic acid concentrations, this dependence will be reversed. The same principle applies to charged particle coagulation loss to larger charged particles of opposite polarity.65 Our data show that the turning point for HIO3 is at approximately 1.5 × 107 cm−3. This explains, at least partly, why the beam enhancement in the experiment presented in Fig. 2 A1 (HIO3 = 2 × 107 cm−3) is rather weak.
Finally, the sink due to larger particles (coagulation loss) and deposition has the expected effect of decreasing the formation rate (Fig. 3C). The default sink value used in the simulations at +10 °C is 0.002 s−1 which originates from the wall loss rate of CLOUD.52 It is important to note that increasing the sink from its default value has a more significant impact on overall formation rates compared to decreasing the sink by the same factor at HIO3 <2 × 107 cm−3. This is because the charged cluster formation process is controlled by two major sinks in the atmosphere: ion–ion recombination and coagulation with existing neutral particles. Neither of these processes depends on the concentration of HIO3. Conversely, the charged cluster growth process is dominated by HIO3 condensation. Therefore, at lower HIO3 concentrations, increasing the sink would have a more significant impact on reducing the cluster formation rate.
In order to estimate the collective impact of these parameters and to explore the potential significance of J± across different altitudes of the marine atmosphere, we utilised the conditions stipulated by the international standard atmosphere.56 This allowed us to adjust the temperature and pressure in accordance with prescribed values and subsequently compute the formation rate at various altitudes (depicted in Fig. 4). At the sea level (0 km), the pressure and temperature stand at 1013 hPa and 15 °C, respectively. At 5 km, these values change to 540 hPa and −17.5 °C; at 10 km, they are 264 hPa and −50 °C; and at 15 km, they reach 120 hPa and −56.5 °C. As altitude increases, both temperature and pressure decrease, while the ionisation rate due to galactic cosmic rays increases. Additionally, in Fig. 4B, we introduced variation to the sink, simulating its highest magnitude at the ground level and subsequent decrease at higher altitudes.
Fig. 4 The simulated formation rate of charged clusters under conditions representing different altitudes. The pressure and temperature were adjusted according to the international standard atmosphere56 and the ionisation rate increases with increasing altitude (A). The pressure and temperature at the different altitudes are 1013 hPa and 15 °C at 0 km; 540 hPa and −17.5 °C at 5 km; 264 hPa and −50 °C at 10 km and at 15 km it is 120 hPa and −56.5 °C. In figure (B) the sink effect is added. |
The ionisation rate has the greatest influence on J± across the range of altitudes simulated here. At low HIO3 concentrations, the increasing ionisation rate at higher altitudes leads to a decreased J±, whereas the opposite is the case at high HIO3 concentrations (>1.5 × 107 cm−3, Fig. 4A). In Fig. 4B, the sink is also considered, which decreases the formation rate, especially at the lowest altitudes and thus evens out the simulated formation rates at different parts of the atmosphere. Since the HIO3 concentration in marine atmospheres (excluding the coastal regions) is likely close to or below 1 × 107 cm−3,16,40 the low HIO3 concentrations are likely mechanistically more interesting.
All in all, our results indicate that enhancing the ionisation rate may not necessarily accelerate new particle formation since it is also dependent on the chemical system. One unique thing about ion-induced HIO3 nucleation is that it proceeds at the collision limit at temperatures of +10 °C and below, i.e., every collision of HIO3 with charged clusters leads to irreversible growth,16,52 which is not the case for systems such as sulphuric acid and oxidised organic vapours. Hence, in the case of the iodic acid system, a threshold of 1.5 × 107 cm−3 is observed, beyond which increased ion concentrations result in higher rates of charged cluster nucleation. It is important to note that different systems may exhibit varying thresholds.
Our results further indicate that ion-induced nucleation probably does not play a significant role in iodine oxoacid particle formation in the upper troposphere–lower stratosphere, where the low temperatures favour the neutral iodine oxoacid nucleation, and the formation rate of charged clusters (J±) is low according to our simulations. It has to be noted that our model only considers the charged clusters that survive until the critical size, and does not consider the possible enhancement of ion-induced nucleation by ion–ion recombination processes. However, we expect that the combined change of J± and Jrec (i.e., Jiin) is still insignificant compared to the enhancement of neutral formation rate at lower temperatures. This is supported by Fig. 2, as panel B1 clearly shows that the increase in ionisation rate has minimal effect on the total formation rate (Jtot) at −10 °C. This may well extend to lower temperatures as neutral JN will likely be further enhanced at lower temperatures while Jiin changes moderately.
On the other hand, it is also clear from our data that ion-induced nucleation can significantly enhance particle formation rates from HIO3 at temperatures above 0 °C and iodic acid concentrations below about 5 × 107 cm−3,16 which makes the ion-induced nucleation mechanism important in large parts of coastal and marine lower atmospheres.
Fig. 5 Total formation rate, Jtot at 1.7 nm from GCR and Neutral runs at CLOUD (same data as in He et al. (2021)16) compared to the charged cluster formation rates from the simplified PANDA model at +10 °C, as well as Jtot at 2.5 nm from field measurements in the Arctic Ocean.24 The error bars are the standard deviation of the measurement data. Shaded area is the formation rate of charged clusters for HIO3 clusters (lower limit) and HIO3 + 2H2O clusters (upper limit) at 1.7 nm. The model data agree well with the CLOUD GCR points at +10 °C below 5 × 107 cm−3, where the ion-induced nucleation dominates. |
Finally, we compare iodine oxoacid formation rates from CLOUD and model data to field observations. The challenge is that it is difficult to find field data with necessary instrumentation and where we can safely assume that the nucleation and early growth are dominated by iodine oxoacids, since very often other potential precursor vapours such as sulphuric acid, methanesulphonic acid and/or various organic compounds are present simultaneously.25 However, Baccarini et al. (2020)24 reported observations of new particle formation driven by iodic acid in the central Arctic Ocean. We used the data from Baccarini et al. (2020)24 to calculate J2.5 for a subset of reported particle formation events, where particles were present in the smallest size channels of the NAIS so that we could assume that they were formed at or very close to the measurement location (Fig. 5, triangle symbols). It was not possible to calculate J1.7 directly, since no instruments that measured this size range were present. Given the uncertainties in the instrumentation, slightly different size range (J2.5 is commonly smaller than J1.7) and method to calculate the formation rates (see Section 2.4), we consider that the agreement is reasonable. The observed particle formation rates in the Arctic Ocean at temperatures from −1 °C to −9 °C fall mostly between the CLOUD data at +10 °C and −10 °C, and a little bit above the simulated J±. This is in agreement with the interpretation of Baccarini et al. (2020)24 that particle formation during the reported events proceeds as negative ion-induced iodic acid nucleation, with a potential contribution from the neutral channel as well.
For RH >2%, iodine oxoacid particle formation has no obvious dependence on RH. This feature is distinct from the nucleation processes from sulphuric acid and oxidised organic vapours, where a strong RH dependence has been observed.39,66 Also, Jtot can be almost perfectly constrained by iodine oxoacid concentrations (Fig. 6A–C), consistent with the results from our previous studies.16,35 This holds both at low HIO3, where ion-induced nucleation mechanisms dominate, and at high (>5 × 107 cm−3) HIO3 where neutral processes dominate, as shown in Section 3.3.
Interestingly, the formation rates in the dry experiments behave very differently compared to the data points with RH >2%. These represent the driest reachable conditions in the CLOUD chamber with an RH of below 0.008% (or absolute water molecule concentration of 2.52 × 1013 cm−3). Such extremely dry conditions are achieved by using certified cryogenic air which was used to purge the chamber for over a week before the dry experiments.46 The particle formation rates observed in the dry experiments are much higher than in the other experiments (Fig. 6) and these extremely high rates cannot be explained by the iodine oxoacid nucleation mechanism.16,35 In the dry experiments, iodine oxide concentrations (e.g., I2O4 and I2O5) are much higher, which suggests that dry conditions favour the formation of higher iodine oxide concentrations which leads to accelerated particle formation. However, it is important to emphasise that the water concentrations in the dry experiments are not representative of typical tropospheric conditions. Consequently, they are presented here solely as a mechanistic reference.
As soon as we introduced water vapour into the chamber (RH = 2%), I2O4 and I2O5 concentrations dropped sharply, especially the I2O5 concentration decreased by two orders of magnitude. This is expected since I2O5 reacts with water to form iodine oxoacids.32 It should be noted that when the RH is higher than 2%, the I2O5 concentration is too low to be relevant for aerosol nucleation (<2 × 105) and the I2O4 concentration is lower than the HIO2 concentration. The potential role of I2O4 in iodine particle formation through, e.g., HIO3–I2O4 nucleation may be worth considering as I2O4 has been thought to be driving iodine particle formation.28 However, our data indicate a minor contribution of I2O4 to iodine particle formation. First, the correlation between Jtot and I2O4 is clearly impacted by the RH since higher RH lowers the concentration of I2O4 while enhancing Jtot at the same time. This should be the opposite if I2O4 were to contribute to Jtot significantly.
The negative dependence of the I2O4 concentration on RH additionally suggests that its role in the marine boundary layer might be small. For example, at RH = 80%, the average I2O4:HIO3 and I2O4:HIO2 ratios are 0.005 and 0.146, respectively, and one may expect that these ratios may be even lower in the marine boundary layer with an RH frequently higher than 80%.
All in all, our results suggest that the iodine oxoacid nucleation is the dominant iodine particle formation mechanism in the marine and polar boundary layer for all relevant temperatures and humidities and that larger iodine oxides such as I2O4 and I2O5 play minor roles.
A kinetic model (PANDA520) was used to simulate the charged cluster formation rate of iodic acid and its dependency on temperature, ionisation rate, and sinks. The results showed that the formation rate of charged clusters is not significantly dependent on temperature, but it is strongly dependent on the ionisation rate. J± increases with increasing ionisation rate, but only if the iodic acid concentration is above ca. 1.5 × 107 cm−3. This is because with increasing ionisation rate, the ion concentration and therefore the recombination loss also increases. When the HIO3 concentration exceeds 1.5 × 107 cm−3, along with a higher ionisation rate, ion growth due to collisions increases more rapidly than recombination loss. This results in an elevated formation rate. The model results agreed both qualitatively and quantitatively with the chamber experiments at +10 °C. On the other hand, as the neutral nucleation of iodic acid is strongly dependent on temperature,16 the neutral pathway became dominant at −10 °C and the total formation rate was insensitive to the ionisation rate.
This leads to the conclusion that despite higher ionisation rates, charged cluster nucleation of iodic acid is probably negligible at high altitudes with cold temperatures and low iodic acid concentrations. However, in the marine boundary layer, ion-induced iodic acid nucleation is important, as suggested also by the few existing field observations.24
The formation rates showed no dependence on relative humidity, as long as RH was above 2%. Therefore, the formation rates could be constrained by the iodine oxoacid mechanism (HIO3 × HIO2) over a wide range of RH and acid concentrations. I2O4 and I2O5 play a significant role in the nucleation process, only if the atmosphere is extremely dry (below 0.008% at +10 °C). These desiccated conditions do not accurately represent tropospheric conditions. As a result, our findings suggest that the iodine oxoacid nucleation mechanism is the predominant pathway for iodine nucleation within the marine boundary layer atmosphere.
This journal is © The Royal Society of Chemistry 2024 |