Impacts of coagulation on the appearance time method for new particle growth rate evaluation and their corrections
 ^{1}Aerosol and Haze Laboratory, Beijing Advanced Innovation Center for Soft Matter Science and Engineering, Beijing University of Chemical Technology, Beijing, 100029, China
 ^{2}Institute for Atmospheric and Earth System Research/Physics, Faculty of Science, University of Helsinki, 00014 Helsinki, Finland
 ^{3}School of Environmental Science and Engineering, Shanghai Jiao Tong University, Shanghai, 200240, China
 ^{4}State Key Joint Laboratory of Environmental Simulation and Pollution Control, School of Environment, Tsinghua University, Beijing, 100084, China
 ^{5}Shanghai Key Laboratory of Atmospheric Particle Pollution and Prevention (LAP3), Department of Environmental Science and Engineering, Fudan University, Shanghai, 200433, China
 ^{1}Aerosol and Haze Laboratory, Beijing Advanced Innovation Center for Soft Matter Science and Engineering, Beijing University of Chemical Technology, Beijing, 100029, China
 ^{2}Institute for Atmospheric and Earth System Research/Physics, Faculty of Science, University of Helsinki, 00014 Helsinki, Finland
 ^{3}School of Environmental Science and Engineering, Shanghai Jiao Tong University, Shanghai, 200240, China
 ^{4}State Key Joint Laboratory of Environmental Simulation and Pollution Control, School of Environment, Tsinghua University, Beijing, 100084, China
 ^{5}Shanghai Key Laboratory of Atmospheric Particle Pollution and Prevention (LAP3), Department of Environmental Science and Engineering, Fudan University, Shanghai, 200433, China
Correspondence: Runlong Cai (runlong.cai@helsinki.fi)
Hide author detailsCorrespondence: Runlong Cai (runlong.cai@helsinki.fi)
The growth rate of atmospheric new particles is a key parameter that determines their survival probability of becoming cloud condensation nuclei and hence their impact on the climate. There have been several methods to estimate the new particle growth rate. However, due to the impact of coagulation and measurement uncertainties, it is still challenging to estimate the initial growth rate of new particles, especially in polluted environments with high background aerosol concentrations. In this study, we explore the influences of coagulation on the appearance time method to estimate the growth rate of sub3 nm particles. The principle of the appearance time method and the impacts of coagulation on the retrieved growth rate are clarified via derivations. New formulae in both discrete and continuous spaces are proposed to correct for the impacts of coagulation. Aerosol dynamic models are used to test the new formulae. New particle formation in urban Beijing is used to illustrate the importance of considering the impacts of coagulation on the sub3 nm particle growth rate and its calculation. We show that the conventional appearance time method needs to be corrected when the impacts of coagulation sink, coagulation source, and particle coagulation growth are nonnegligible compared to the condensation growth. Under the simulation conditions with a constant concentration of nonvolatile vapors, the corrected growth rate agrees with the theoretical growth rates. However, the uncorrected parameters, e.g., vapor evaporation and the variation in vapor concentration, may impact the growth rate obtained with the appearance time method. Under the simulation conditions with a varying vapor concentration, the average bias in the corrected 1.5–3 nm particle growth rate ranges from 6 %–44 %, and the maximum bias in the sizedependent growth rate is 150 %. During the test new particle formation event in urban Beijing, the corrected condensation growth rate of sub3 nm particles was in accordance with the growth rate contributed by sulfuric acid condensation, whereas the conventional appearance time method overestimated the condensation growth rate of 1.5 nm particles by 80 %.
 Article
(4177 KB) 
Supplement
(1246 KB)  BibTeX
 EndNote
New particle formation (NPF) is frequently observed in various atmospheric environments (Kulmala et al., 2004; Kerminen et al., 2018; Nieminen et al., 2018; Lee et al., 2019). It contributes significantly to the number concentrations of aerosols and cloud condensation nuclei (CCN) and hence impacts the global climate (Kuang et al., 2009; Kerminen et al., 2012). The new particle growth rate is one of the key parameters to characterize NPF events. On the one hand, the newly formed particles (∼ 1 nm) have to survive coagulation scavenging before they grow to the CCN size (∼ 100 nm). Given the same background aerosol concentration, i.e., the same coagulation loss rate, it is the growth rate that determines the survival probability of new particles (Weber et al., 1997; Lehtinen et al., 2007). Therefore, measuring the new particle growth rate accurately contributes to understanding the impact of NPF on the climate. On the other hand, the particle growth rate is a key to investigating growth mechanisms. Theoretical particle growth rates contributed by condensing vapors are usually compared to measured growth rates to reveal the possible particle growth mechanisms (Ehn et al., 2014; Yao et al., 2018; Mohr et al., 2019). A nonbiased and accurate determination of measured growth rates is an important fundament of these comparisons.
Although new particle growth rates are frequently reported in various environments around the world, it remains difficult to retrieve accurate particle growth rates from an ambient dataset. Due to the varying atmospheric conditions, significant Kelvin effect, and sizedependent particle compositions, the particle growth rate is a function of both time and particle size. The measured evolution of aerosol size distribution does not directly indicate the sizeandtimeresolved growth rate of single particles because one cannot directly track single particles from the size distributions. There are several methods to obtain the sizeandtimeresolved growth rate by solving aerosol general dynamic equations (GDE; Kuang et al., 2012; Pichelstorfer et al., 2018). However, few applications of these GDE methods have been reported for particle growth analysis in the real atmosphere (e.g., Kuang et al., 2012). The most likely reason for this is that these GDE methods are sensitive to measurement uncertainties caused by atmospheric instability and instruments, which needs to be solved in future studies.
Apart from solving the GDEs, the widely used methods to estimate the particle growth rate are based on finding the representative particle diameter or time. The representativediameter method usually uses the peak diameter of the size distribution of new particles and estimates its increase rate from its temporal evolution. The rate of increase in peak diameter is then taken as the particle growth rate (Kulmala et al., 2012) after correcting for (or sometimes neglecting) the influence of coagulation on the peak shifting (Stolzenburg et al., 2005). During the correction, coagulation is often classified into innermodal coagulation (selfcoagulation) and intermodal coagulation (Anttila et al., 2010; Kerminen et al., 2018). The peak diameter is usually obtained by fitting a lognormal function to the measured aerosol size distribution of new particles. With a distinct peak diameter in the growing particle population, this method is theoretically feasible for estimating new particle formation rates. However, the mode fitting is usually tricky, especially when there is no welldefined mode in the growing distribution, due to either the aerosol distribution itself or the measurement uncertainties.
The representativetime method estimates the corresponding time for a series of diameters according to a certain criterion and then calculates the growth rate according to the relationship between the diameters and their corresponding time (Dada et al., 2020). The corresponding time is determined as the time to reach either the maximum concentration (maximum concentration method; Lehtinen and Kulmala, 2003) or a certain proportion of the maximum concentration (appearance time method; Lehtipalo et al., 2014) of a given particle size bin. Previous studies have tested the appearance time under various modeling conditions. Their results indicate that some appearance time methods can reproduce the theoretical growth rate within acceptable uncertainties under certain test conditions (Lehtipalo et al., 2014) but not under other test conditions (Olenius et al., 2014; Kontkanen et al., 2016; Li and McMurry, 2018). As shown in the “Theory” section below, the discrepancy is because the slope of particle size against appearance time usually convolves other information (e.g., coagulation) in addition to particle growth.
Determining the growth rate of sub3 nm particles is more challenging than determining that of larger particles. Firstly, there are considerable uncertainties in the measured sub3 nm aerosol size distributions (Kangasluoma et al., 2020) compared to largersized particles (e.g., > 10 nm; Wiedensohler et al., 2012). These uncertainties pose a great challenge to the methods based on solving aerosol general dynamic equations. Secondly, during a typical atmospheric NPF event, the sub3 nm particle size distribution function usually decreases monotonically with the increasing diameter (Jiang et al., 2011b). As a result, the representativediameter method usually faces difficulties in covering the sub3 nm size range. In contrast, despite lacking a clear mathematical understanding of the information convolved in the slope of appearance time against particle diameter, the appearance time method is usually favored for sub3 nm particles and clusters because of the existence of the concentration peak of new particles during an atmospheric NPF event. In addition, the appearance time method is not significantly affected by the systematic instrumental uncertainties because the appearance time of each size bin is only determined by the relative signal rather than by the absolute particle concentration.
Coagulation impacts both particle growth and the growth rate calculation, especially for polluted environments and some chamber studies with high aerosol concentrations. The impact of coagulation on aerosol dynamics has been known for decades (e.g., McMurry, 1983). Recent studies have discussed the importance of considering coagulation when estimating the new particle growth rate (Cai and Jiang, 2017), the influence of transport on measured size distributions (Cai et al., 2018), and primary particle emissions (Kontkanen et al., 2020) under a high aerosol concentration. Similarly, neglecting particle coagulation may cause a bias in the retrieved particle growth rate. Therefore, the coagulation growth has to be considered before investigating the contributions of various condensing vapors to particle growth.
In this study, the feasibility and limitations of the appearance time method are investigated based on theoretical derivations. The impact of coagulation on the retrieved growth rate using the appearance time method is explored and then corrected for. Aerosol dynamic models are used to test the conventional and corrected methods. After that, the corrected appearance time method is applied in a typical NPF event in urban Beijing to show the impact of coagulation on growth rate evaluation in the real atmosphere. In addition, the potential uncertainties in the corrected appearance time method caused by vapor evaporation and the variation in vapor concentration are discussed.
2.1 Particle growth rate
Before deriving the formulae for the appearance time method, the definitions of particle growth and coagulation loss have to be clarified to avoid potential misunderstanding. Although widely used in NPF analyses, the exact meanings of these two concepts vary with their applied conditions.
The particle growth rate, by definition, is the rate of increase in particle diameter as a function of time for a given particle. Assuming that there is a sufficient number of particles of the same size and compositions, it is reasonable to neglect the influence of the stochastic effect due to a low particle number and use the expectation of the singleparticle growth rate to characterize the growth of the aerosol population with the same size. When there is only one nonvolatile condensing vapor, the formula for the expectation of the singleparticle condensation growth rate (referred to as the condensation growth rate below for simplicity) is shown in Eq. (1):
where GR_{cond} is the condensation growth rate (nm s^{−1}) that neglects evaporation, d_{p} is the particle diameter (nm), t is time (s), d_{1} is the diameter of the condensing vapor (nm), β_{1,p} is the coagulation coefficient between d_{1} and d_{p} (cm^{3} s^{−1}), and N_{1} is the vapor concentration (cm^{−3}). Particle evaporation is assumed to be negligible, and the particle shape is assumed to be spherical both before and after the growth. Note that Eq. (1) is expressed in the discrete form; i.e., it does not assume a continuum particle size (d_{1}→0). When multiple vapors contribute to particle growth simultaneously, the total condensation growth rate is the sum of the condensation growth rates contributed by every single vapor. When considering particle evaporation, i.e., monomer dissociation, particle growth due to the net effect of vapor association and dissociation will be explicitly referred to as net condensation growth.
In addition to the condensation of vapors, coagulation also contributes to particle growth. For a given particle with the size of d_{p}, the coagulation with a particle much smaller than d_{p} is usually considered as a contribution to its growth. In contrast, the coagulation with a particle much larger than d_{p} is usually considered the coagulation loss of particle d_{p}. We follow this convention to distinguish between coagulation growth and loss; i.e., particle coagulation with another particle no larger than itself is taken as coagulation growth, and otherwise, it is taken as coagulation loss. Hence, the formula for the expectation of the singleparticle coagulation growth rate (referred to as the coagulation growth rate for simplicity) in the discrete form is
where GR_{coag} is the coagulation growth rate (nm s^{−1}), d_{min} is the minimum particle size (nm), β_{p,i} is the coagulation coefficient (cm^{3} s^{−1}) between d_{p} and d_{i}, and N_{i} is the concentration (cm^{−3}) of particles with the size d_{i}. Since both condensation and coagulation contribute to particle growth, the total singleparticle growth rate is equal to the sum of GR_{cond} and GR_{coag}.
When retrieving the particle growth rate from the measured aerosol size distributions, the retrieved value is named the apparent growth rate. “Apparent” emphasizes that the method does not necessarily guarantee that the retrieved growth rate is equal to the condensation or total growth rate of a single particle or the investigated aerosol population. When using the representativediameter method, the retrieved apparent growth rate is the increasing rate of the peak diameter and does not directly characterize the growth of any particle(s). For instance, the coagulation loss rate is a function of particle diameter; as a result, the peak diameter shifts towards larger sizes with time because smaller particles are scavenged faster by coagulation than larger particles. Similarly, other sizedependent processes such as condensation and coagulation growth also cause the shift of peak diameter. As a result, the apparent growth rate sometimes needs to be corrected before it is taken as the total growth rate or the condensation growth rate (Stolzenburg et al., 2005). When using the representativetime method, although the retrieved apparent growth rate is close to the condensation growth rate under some modeling conditions (Lehtipalo et al., 2014), their deviation can be significant under other conditions, on which we elaborate in Sect. 4.
2.2 Coagulation sink and source
For a given particle, its coagulation with another particle can be classified into coagulation growth and coagulation loss as aforementioned. This classification is based on the Lagrangian specification that tracks the growth of a single particle. In contrast, according to the Eulerian specification that focuses on given particle diameters, each coagulation causes a sink of two particles and a source of one new particle with a larger diameter regardless of the particle sizes. Herein, we define the coagulation sink and source as the loss and production rate for particle size bins in the Eulerian specification. According to these definitions, the coagulation of a large particle with another smaller particle is counted as the coagulation sink (in the Eulerian specification) but not as the coagulation loss (in the Lagrangian specification) of the large particle. Following previous studies, we use CoagS (s^{−1}) to represent the sink coefficient (Kulmala et al., 2001) and CoagSrc (cm^{−3} s^{−1}) to represent the production rate due to coagulation (Kuang et al., 2012). Their formulae in the discrete form are given below:
where d_{max} is the maximum particle diameter (nm); d_{j} is defined by ${d}_{j}^{\mathrm{3}}={d}_{\text{p}}^{\mathrm{3}}{d}_{i}^{\mathrm{3}}$; N_{i} and N_{j} are the concentrations of d_{i} and d_{j}, respectively; and other variables have been introduced above. Note that CoagS and CoagSrc are defined differently and their units are also different. Table 1 summarizes the differences between the definitions in the Lagrangian and Eulerian specifications.
2.3 Formulae for the new appearance time method
The conventional appearance time method has been discussed in detail in Lehtipalo et al. (2014), and its derivation under ideal conditions has been reported in He et al. (2020). Here we briefly describe the procedure to retrieve the particle growth rate from the temporal evolution of a measured aerosol size distribution using the appearance time method. For each aerosol size bin, its corresponding appearance time is determined as the moment that the measured aerosol concentration in this size bin reached 50 % or any other given proportion of its maximum concentration during the event. The maximum and halfmaximum concentration of this size bin can be either taken from the smoothed temporal evolution of the measured concentration or determined by fitting a sigmoid function to the measured data. The growth rate is then estimated as the slope of the diameter of aerosol size bins versus their corresponding appearance time; i.e.,
where GR_{conv} is the total growth rate (nm s^{−1}) retrieved by the conventional appearance time method, Δd_{p} (nm) is the size difference between two adjacent measured size bins, and Δt (s) is the time difference in the appearance time of these two size bins. Note that since the appearance time is a function of d_{p} and GR_{conv} is estimated for different appearance times, the apparent relationship between GR_{conv} and d_{p} does not necessarily indicate the size dependency of GR_{conv} at any given moment.
The correction formulae for influences of coagulation on the appearance time method in the discrete space are
where GR_{corr,tot} is the total growth rate (nm s^{−1}) after correcting for the impact of both coagulation sink and source; GR_{corr,cond} is the condensation growth rate (nm s^{−1}) after correction; GR_{coag} is the coagulation growth rate (nm s^{−1}); CoagS (s^{−1}) and CoagSrc (cm^{−3} s^{−1}) are the coagulation sink and coagulation source term for d_{p}, respectively; and N_{p} is the number concentration of particles with the size d_{p} at their appearance time. Note that coagulation growth is corrected in Eq. (7) but not Eq. (6).
When measuring aerosol size distribution using size spectrometers, the measured distributions are usually reported in a certain number of sectional bins. Therefore, in addition to the formula in the discrete form (Eq. 6), the correction formulae for the appearance time in the sectional form are given below:
where d_{p,u} (nm) and d_{p,l} (nm) are the upper and lower size limits of a given size bin, d_{p} (nm) is the representative diameter (usually the geometric mean diameter) of this size bin, CoagS is the coagulation sink (s^{−1}) for d_{p}, CoagSrc is the coagulation source term (cm^{−3} s^{−1}) for the given size bin, ${N}_{[{d}_{\text{p,l}},{d}_{\text{p,u}}]}$ is the measured concentration (cm^{−3}) of the size bin [d_{p,l}, d_{p,u}] at t_{p}, d_{1} is the diameter (nm) of the condensing vapor, and n_{i} is the aerosol size distribution function (dN_{i} $/$ dlogd_{i}, cm^{−3}, where N_{i} is the cumulative distribution) for the given size d_{i}.
The derivation for Eq. (6) is detailed in Sect. 4.1 and 4.2. The new and conventional appearance time methods are tested using a discretesectional model in Sect. 4.3 and a measured atmospheric NPF event in Sect. 4.4.
The proposed correction formulae for influences of coagulation on the appearance time method (Eqs. 6 and 7) are derived in Sect. 4.1 and 4.2. To test the validity of these corrections, numerical models are used to simulate an evolving aerosol size distribution and provide the theoretical growth rate according to the input monomer concentration. The growth rate is retrieved from the simulated aerosol size distribution using the conventional and corrected appearance time methods and then compared to the theoretical growth rate. Since the correction formulae are proposed only for the influence of coagulation, the influences of evaporation and varying vapor concentration on the retrieved growth rate with the appearance time method are discussed in Sect. 4.3. After that, the conventional and corrected appearance time methods are applied to a new particle formation event measured in urban Beijing. The impact of coagulation on the appearance time of new particles and hence growth rates are indicated by the differences between the growth rates retrieved from the measured aerosol size distributions using the conventional and corrected appearance time methods.
3.1 Numerical models
A discrete aerosol model and a discretesectional aerosol model based on aerosol dynamics were used to provide an evolving aerosol size distribution and hence to test the conventional and corrected appearance time methods. The discrete model assumed that new particle formation is driven by the nucleation and condensation of a certain singlecomponent condensing vapor (H_{2}SO_{4}). The vapor concentration was set as a constant. Condensation, coagulation, external loss, and evaporation were considered in this discrete model. The concentrations of particles up to the size of 100 vapor molecules (∼ 3.5 nm) were numerically solved using the Julia programming language. The theoretical condensation and coagulation growth rates were calculated using Eqs. (1) and (2).
The discretesectional model was composed of 30 discrete bins (up to 2.2 nm) and 400 sectional bins (up to 230 nm). The vapor concentration was assumed to follow a normal distribution to simulate its diurnal variation in the real atmosphere. A growth enhancement factor as a function of particle size (Kuang et al., 2010) was used to account for the condensation of multiple vapors. A certain concentration of 100 nm particles was used as background particles, and their concentration and size were kept constant during each simulation. The discretesection model was coded in MATLAB, and it is detailed in Li and Cai (2020).
The simulation conditions for Figs. 1 and 3–7 are summarized in Table 2. The simulations with varying vapor concentration are summarized in Table A1. Note that the aerosol dynamic models are only used to provide a benchmark to compare the conventional and corrected appearance time methods in this study.
3.2 Measurements
The NPF event measured on 24 February 2018, in urban Beijing, was used to test the influences of coagulation on the appearance time method. During 08:00–16:00 (local time, UTC+8) on this NPF day, the mean temperature, relative humidity, and wind speed were 1.3 ^{∘}C, 22 %, and 1.4 m s^{−1}, respectively.
The measurement site is located on the west campus of the Beijing University of Chemical Technology, which is close to the west 3rd Ring Road of Beijing. The aerosol size distributions were measured using a homemade particle size distribution system (PSD; Liu et al., 2016) and a homemade diethylene glycol scanning mobility particle spectrometer (Jiang et al., 2011a; DEGSMPS; Cai et al., 2017) equipped with a coresampling apparatus (Fu et al., 2019). The sulfuric acid monomer and dimer concentrations were measured using a long chemical ionization timeofflight mass spectrometer (ToFCIMS; Aerodyne Research, Inc.; Jokinen et al., 2012). The meteorological data were measured using a local weather station (Vaisala, AWS310). More details on this measurement site and the instruments have been introduced elsewhere (Deng et al., 2020b).
4.1 Appearance time method under ideal conditions
Prior to investigating the impacts of coagulation on the appearance time method, we briefly illustrate the principle of the appearance time method. He et al. (2020) recently demonstrated that the appearance time method can retrieve the condensation growth rate under ideal conditions. The ideal conditions are as follows:

The vapor concentration is constant.

The initial concentrations of new particles are equal to zero.

Condensation is the only cause of the change in particle concentrations; i.e., there is no coagulation, evaporation, external loss, etc.

The condensation rate (i.e., coagulation rate between vapor and particles) is independent of the particle diameter.
It can be seen that none of these ideal conditions is consistent with real atmospheric environments. Since the conventional appearance time method is derived based on these conditions, violating them may cause biases in the appearance time method. We will first show brief derivations of the conventional appearance time and then discuss the correction for the influences of coagulation and other remaining potential uncertainties in the corrected appearance time method.
Under the above ideal conditions, the population balance equation for a particle containing i molecules is
where N_{1} is the concentration (cm^{−3})) of the condensing vapor (monomer); N_{i−1}, and N_{i} are the concentrations (cm^{−3}) of particles containing i and i−1 monomer molecules, respectively; t is time (s); and β is the coagulation coefficient (cm^{3} s^{−1}) between a vapor molecule and any particle, which is assumed to be independent of the particle size in Eq. (10). For the case i=2, the last term in Eq. (10) should be modified as 0.5$\mathit{\beta}{N}_{\mathrm{1}}^{\mathrm{2}}$.
Solving the differential equations in Eq. (10) yields the analytical solution for N_{i}:
where N_{i,∞} is the concentration limit (cm^{−3}) of N_{i} when t approaches infinity (dN_{i} $/$ dt=0), and it is equal to 0.5$\mathit{\beta}{N}_{\mathrm{1}}^{\mathrm{2}}$ (for i>1) under these ideal conditions.
Figure 1a shows the concentrations of N_{i} normalized by dividing by their corresponding N_{i,∞}. It can be seen that the distance between two adjacent concentration curves is approximately a constant though these curves are not parallel. Hence, the appearance time method takes the moment that N_{i} reaches a certain percentage of its maximum value (N_{i,∞}) as its representative time. Previous studies indicate that the 50 % sizeresolved appearance time method which chooses the certain percent as 50 % is more robust against nonideal conditions compared to using the criteria of other percent values (Lehtipalo et al., 2014). As shown in Fig. 2, an approximate solution of t for ${N}_{i}\left(t\right)=\mathrm{0.5}{N}_{i,\mathrm{\infty}}$ (referred to as t_{i}) is
Equation (12) indicates that the slope of particle size in terms of its molecule number versus its appearance time is approximately equal to its condensation growth rate; i.e.,
where GR_{cond,n} is the condensation growth rate in terms of the molecule number (s^{−1}). The relationship between GR_{cond,n} herein and GR_{cond} in Eq. (1) (which is defined with respect to particle diameter) is GR${}_{\text{cond}}={\text{GR}}_{\text{cond},n}\times \mathrm{\Delta}{d}_{i}$, where Δd_{i} (nm) is the increase in d_{i} due to the condensation of one vapor molecule.
According to the derivations above, the 50 % sizeresolved appearance time (referred to as 50 % appearance time for short) method can retrieve the particle growth rate under the given ideal conditions. The slope of particle size versus the appearance time is approximately equal to the condensation growth rate. That is, this slope is mainly determined by condensation growth under these ideal conditions. Note that Eq. (13) is only valid for the 50 % appearance time, whereas other thresholds to determine the appearance time may cause systematic bias. This bias comes from the nonparallelism of particle concentration curves (see Fig. 1). As shown in Fig. 2, in this test case, the 5 % appearance time method overestimates the growth rate by 15 % and the 95 % appearance time method underestimates the growth rate by 12 %. It should be clarified that since these biases are not huge, it is acceptable to use other thresholds instead of 50 % to reduce the impact of measurement uncertainties in practical applications.
Particle evaporation is assumed to be negligible in the above derivations, yet Fig. 2 indicates that sizeindependent evaporation does not significantly impact the validity of the 50 % appearance time method. Assuming a sizeindependent evaporation rate, E (s^{−1}), t_{i} is approximately equal to (ln2+i) $/$ (βN_{1}−E), which is bigger than that without evaporation. Meanwhile, considering evaporation, the net condensation growth rate (in the Lagrangian specification) is equal to the vapor condensation rate subtracted by the particle evaporation rate, i.e., βN_{1}−E. That is, the increase in appearance time agrees with the decrease in the net condensation growth rate. In practice, the particle evaporation rate is usually size dependent due to the significant Kelvin effect and it impacts N_{i,∞} in addition to the net condensation growth rate. It will be shown below in Sect. 4.3 that for particles close to the critical size (corresponding to $\mathit{\beta}{N}_{\mathrm{1}}E=\mathrm{0}$), evaporation may cause a substantial bias in the net condensation growth rate retrieved using the appearance time method.
Note that the equality between the slope of particle size versus the 50 % appearance time and the net condensation growth rate holds only under the above ideal conditions. The following derivations and results in Sect. 4.2 will show how the slope is affected by coagulation while maintaining the same condensation growth rate.
4.2 The impacts of coagulation and their corrections
We first show the impact of an external sink on the appearance time method. The external sink is herein referred to as the sink due to coagulation with background particles, wall loss, dilution, transport, etc. For the convenience of comparison with Fig. 1a, the external sink is assumed to be temporally independent of the particle diameter. The impact of its size dependency will be discussed later. Considering the constant external sink, the population balance equation for N_{i} is
where ES is the external sink (s^{−1}) and other variables have been introduced in Eq. (10).
Similarly to Eqs. (11) and (12), the approximate solutions for N_{i} and its corresponding appearance time (t_{i}) are
Equations (15) and (16) indicate that the impact of ES on t_{i} is mathematically equivalent to vapor condensation (βN_{1}). In the presence of a nonnegligible external sink, the particle concentration will approach its limit faster than the scenario without an external sink (Fig. 1b). As a result, the slope of particle diameter versus appearance time is affected by both condensation growth and the external sink.
Combining Eqs. (1), (13), and (16), the impact of the external sink can be readily corrected for. The correction formula is
where GR_{EScorr} is the growth rate (nm s^{−1}) after correcting for the external sink, GR_{conv} is the growth rate (nm s^{−1}) retrieved by the conventional appearance time method, and d_{1} is the diameter (nm) of the condensing vapor.
As shown in Fig. 3a, with a vapor concentration of 5×10^{6} cm^{−3} and an external sink ranging from $\mathrm{1}\times {\mathrm{10}}^{\mathrm{3}}$ to $\mathrm{5}\times {\mathrm{10}}^{\mathrm{3}}$ s^{−1}, the conventional appearance time method overestimates the condensation growth rate substantially. Such an overestimation caused by mistaking the external sink for condensation growth has also been reported in previous studies (Olenius et al., 2014; Li and McMurry, 2018). In contrast, the corrected growth rate agrees well with the theoretical condensation growth rate.
Practically, the coagulation coefficient (β) and ES are functions of the particle diameter. For the convenience of illustration, we use the sizedependent coagulation sink, CoagS, as an example to represent the total particle sink due to coagulation, wall loss, dilution, and transport. Similarly to Eq. (15), the approximate analytical solution for N_{i} with the sizedependent β and CoagS is
where β_{1,i} (or β_{1,g}) is the coagulation coefficient between a vapor molecule and a particle containing i (or g) molecules (cm^{3} s^{−1}), CoagS_{i} (or CoagS_{g}) is the coagulation sink of particles containing i (or g) molecules (s^{−1}), and other variables have been introduced above. Correspondingly, the ES term in Eq. (17) should be replaced with CoagS_{i} to correct for the impact of the sizedependent coagulation sink. When deriving Eq. (18), it is assumed that ${\mathit{\beta}}_{\mathrm{1},i}{N}_{\mathrm{1}}+{\text{CoagS}}_{i}$ is close to ${\mathit{\beta}}_{\mathrm{1},i\mathrm{1}}{N}_{\mathrm{1}}+{\text{CoagS}}_{i\mathrm{1}}$. This approximation is reasonable because both β_{1,i} and CoagS_{i} change gradually with the particle size, yet it introduces minor systematic biases in N_{i} and its corresponding appearance time.
As shown in Fig. 3b, when β_{1,i} and CoagS_{i} are size dependent, the corrected appearance time method is still able to reproduce the condensation growth rate. It is assumed that the CoagS_{i} in Fig. 3b is contributed by only the large background particles. Hence, the CoagS_{i} in Fig. 3b is estimated from the condensation sink (CS) using an empirical formula (Eq. 8 in Lehtinen et al., 2007), where CS indicates the condensation loss rate of the vapor.
In addition to the coagulation with another background particle, the coagulation between two new particles also contributes to the CoagS of both these two particles. As explained in Sect. 2.2, no matter how small the coagulating particle is, the coagulation between a given particle and any other particle should be accounted for in CoagS. This is because the appearance time method is derived in the Eulerian specification and the CoagS is defined for a certain particle diameter rather than for a certain particle. In contrast, when focusing on the survival probability of new particles (Weber et al., 1997; Lehtinen et al., 2007), CoagS should be calculated in the Lagrangian specification; i.e., only the coagulation with a larger particle that causes particle loss should be accounted for. To emphasize the difference between the two definitions of CoagS, the corrected growth rates using the total (Eulerian) CoagS and the background (Lagrangian) CoagS are compared in Fig. 4. The particle source due to coagulation is not considered in this comparison. A constant concentration of 100 nm particles is used as the background particles. The background CoagS refers to the sink due to coagulation with all larger particles, including both the background particles and new particles. The total CoagS is calculated using Eq. (3). Note that due to the contribution of new particles, the total CoagS does not follow a simple decreasing trend with the increasing particle diameter (see Fig. B1 in Cai and Jiang, 2017). Hence, the empirical formula (Lehtinen et al., 2007) to generate a sizedependent CoagS in Fig. 4b should not be used for the total CoagS.
As shown in Fig. 4, the condensation growth rate after correcting for the background CoagS is still overestimated. In contrast, the growth rate after correcting for the total CoagS agrees with the theoretical growth rate for particles larger than 1.3 nm. For sub1.3 nm particles, the systematic bias in the growth rate after correcting for the total CoagS is mainly caused by the violation of the assumption that ${\mathit{\beta}}_{\mathrm{1},i}{N}_{\mathrm{1}}+{\text{CoagS}}_{i}$ is close to ${\mathit{\beta}}_{\mathrm{1},i\mathrm{1}}{N}_{\mathrm{1}}+{\text{CoagS}}_{i\mathrm{1}}$. The size dependency of the particle coagulation coefficient under the influence of new particle coagulation increases with decreasing particle size. For instance, under the test conditions, (${\mathit{\beta}}_{\mathrm{1},\mathrm{4}}{N}_{\mathrm{1}}+{\text{CoagS}}_{\mathrm{4}})/({\mathit{\beta}}_{\mathrm{1},\mathrm{3}}{N}_{\mathrm{1}}+{\text{CoagS}}_{\mathrm{3}})=\mathrm{1.12}$, while (${\mathit{\beta}}_{\mathrm{1},\mathrm{50}}{N}_{\mathrm{1}}+{\text{CoagS}}_{\mathrm{50}})/({\mathit{\beta}}_{\mathrm{1},\mathrm{49}}{N}_{\mathrm{1}}+{\text{CoagS}}_{\mathrm{49}})=\mathrm{1.01}$. As a result, the corrected appearance time method still overestimates the growth rate for sub1.3 nm particles.
In addition to CoagS, the impact of the coagulation source on the appearance time also needs correcting for. Adding the coagulation source term to the population balance equation of N_{i} yields
where CoagS_{i} is the coagulation sink (s^{−1}) corresponding to N_{i}, CoagSrc_{i} is the coagulation source term (cm^{−3} s^{−1}) corresponding to N_{i}, and other variables have been introduced above.
Since CoagSrc_{i} is determined by the concentrations of all particles containing 2 to i−2 molecules, it is difficult to obtain an analytical solution of Eq. (19) without approximation. Here we provide an approximation method to correct for the impact of the coagulation source on the appearance time. Compared to the scenario where there is no coagulation source term in the balance equation (Eq. 19), the coagulation source has three impacts on particle size distribution and particle growth: (1) the coagulation with a smaller particle contributes to particle growth; (2) the coagulation source increases the maximum particle concentrations; (3) the coagulation source shortens the time for particles to reach their maximum concentration.
These three impacts of the coagulation source are accounted for in the corrected appearance time method (Eqs. 6 and 7). Impact (1) is corrected for using Eq. (7). To correct for impacts (2) and (3), we simply assume that CoagSrc_{i} is a constant during the increasing period of N_{i} and use Eq. (6) (see the Supplement for the derivation of Eq. 6) to estimate the growth rate.
The corrected appearance time method was tested under various modeling conditions. As the example in Fig. 5 indicates, the growth rates estimated using the corrected appearance time method agree with the theoretical growth rates. Equation (5) is able to retrieve the growth rate of sub3 nm particles unless coagulation source is a governing reason for the change in particle concentration; i.e., CoagSrc $/$ 2N_{p} is comparable or larger than β_{1,p}N_{1}. Under these conditions, CoagSrc may not be a constant, and, hence, the approximation of CoagSrc $/$ 2N_{p} may cause bias. Fortunately, CoagSrc usually decreases with increasing particle size due to the decreasing particle concentration. Furthermore, it will be shown in Sect. 4.4 that CoagSrc does not have a major impact on the apparent growth rate of sub10 nm particles even during an intensive atmospheric NPF event in urban Beijing. Hence, we consider Eq. (6) as a rough but sufficient formula to correct for the impact of coagulation on the appearance time for most atmospheric NPF events.
4.3 Uncertainties in the appearance time method
In the above analysis, the correction for the influences of coagulation on the appearance time method was validated by derivations and simulations. However, there are still potential uncertainties in the corrected appearance time method because of other uncorrected influences. This section will discuss the uncertainties caused by vapor evaporation and a varying vapor concentration.
Due to the uncorrected vapor evaporation, the appearance time method may overestimate the growth rate for when the net condensation rate in the Lagrangian specification (βN_{1}−E) is close to zero or negative. As shown in Fig. 6a, a sizedependent evaporation rate is assumed for the clusters and particles. The evaporation rate of a dimer (i.e., cluster containing two vapor molecules) is assumed to be ∼ 0.2 s^{−1} according to the stability of H_{2}SO_{4}–NH_{3} clusters reported in Myllys et al. (2019). The evaporation rate of particles containing more than three vapor molecules is estimated using the Kelvin equation with arbitrarily assumed vapor saturation pressure. With these assumptions of vapor evaporation, the evaporation rate (E) of a single cluster exceeds its monomer association rate (βN_{1}) in the sub1.5 nm size range (Fig. 6a). The influence of the coagulation source is not accounted for in this test, yet its influences can be corrected for according to the above discussions in Sect. 4.2. As shown in Fig. 6b, the conventional and corrected appearance time method overestimates the net condensation growth rate of sub2 nm particles. This bias is mainly caused by the influence of the sizedependent evaporation on N_{i,∞}. Assuming a prior knowledge of the sizedependent evaporation rate, the net condensation growth rate retrieved using the appearance time method can be corrected using Eq. (21):
where the subscripts i and i+1 indicate the number of molecules contained in the clusters; N_{i,50} and ${N}_{i+\mathrm{1},\mathrm{50}}$ are the concentrations of N_{i} and N_{i+1}, respectively, at the 50 % appearance time of N_{i}; and other variables have been declared above. N_{i,50} and ${N}_{i+\mathrm{1},\mathrm{50}}$ can be approximated by N_{i,∞} and ${N}_{i+\mathrm{1},\mathrm{\infty}}$, respectively. The last term in Eq. (21) corrects for the influence of evaporation on N_{i,∞} and the source of N_{i.}
As shown in Fig. 6b, the appearance time method corrected using Eq. (21) follows the theoretical net condensation growth rate. However, since Eq. (21) requires prior knowledge of the sizedependent evaporation rate, it is difficult to use this a correction for an atmospheric NPF event. Hence, one should note the uncertainties in the appearance time method for particles in the neighborhood of the critical size.
Note that the above discussions are for the case that nucleation and growth are driven by a single volatile vapor. The condensation of a volatile vapor onto a nonvolatile particle may not introduce a significant bias to the growth rate estimated using the appearance time method, as shown in Fig. S1 in the Supplement. This is because with an existing nonvolatile vapor, N_{i,∞} may not be significantly affected by the evaporation of other volatile vapors.
The constant vapor concentration assumption may be valid for some chamber experiments; however, the vapor concentration usually follows a diurnal pattern in the real atmosphere. The varying vapor concentration may impact the appearance time and hence the retrieved apparent growth rate. As reported in previous studies (Lehtipalo et al., 2014; Olenius et al., 2014), the retrieved appearance time is sensitive to the variation in vapor concentration. In the presence of coagulation, it is difficult to correct for the impact of the varying vapor concentration. Herein, we use the discretesectional model to test the uncertainties in the corrected appearance time method under a varying vapor concentration. The vapor concentration is assumed to follow a normal distribution. The condensation sink (10^{−3}–10^{−2} s^{−1}) is contributed simultaneously by a certain number concentration (7.2×10^{2}–7.2×10^{3} cm^{−3}) of 100 nm background particles and the new particles. The growth rate is firstly estimated using the 50 % sizeresolved appearance time method and then corrected using Eq. (8). The 50 % appearance time is herein calculated using the maximum sizeresolved particle concentration because N_{i,∞} is not available, and this approximation introduces biases to the retrieved growth rate. Since the vapor concentration varies with time, the retrieved growth rate characterizes particle growth at both different diameters and different times instead of the sizedependent growth at a certain moment. To keep in accordance with the appearance time method, the theoretical condensation and coagulation growth rates of each d_{p} value are calculated at their corresponding t_{p} values.
In general, neglecting the variation in the vapor concentration introduces biases to the appearance time method. As the example shown in Fig. 7 (test no. 8 in Table A1) indicates, the deviation between the corrected particle growth rate and the theoretical growth rate is smaller than the deviation between the conventional growth rate and the theoretical growth rate. However, for particles larger than ∼ 5 nm and smaller than ∼ 2 nm, the appearance time method overestimates the particle growth rate even after correcting for the impact of coagulation in the test case. These overestimations are caused by approximating the influence of the coagulation source on particle growth with the coagulation source term in Eq. (6). As summarized in Table A1, the relative discrepancy depends on the exact conditions. The average discrepancy of the corrected appearance time method for 1.5–3 nm particles ranges from 6 % to 44 % in the test conditions, which is smaller than that of the conventional method. The maximum sizedependent discrepancy of the corrected appearance time method reaches 150 % under the test conditions.
Although it is difficult to correct for the bias due to the varying vapor concentration, one can try to avoid large uncertainties because the bias seems to follow a certain pattern. Compared to the scenario of an increasing vapor concentration, it is found that the discrepancy between the real and retrieved growth rates are usually larger after the peak time of the vapor concentration. As shown in Fig. 7, the appearance time of ∼ 4.9 nm particles is 12:00 (local time) and a substantial discrepancy between the theoretical and retrieved growth rate is observed for particles larger than 4.9 nm. Fortunately, during a typical atmospheric NPF event, new particles usually grow large before the vapor concentration starts to decrease. To reduce this systematic bias due to a decreasing vapor concentration, we suggest using the other methods, e.g., the representativediameter method to estimate particle growth rate when the vapor concentration decreases.
For particles close to the size of vapor molecules (sub2 nm in these tests), the appearance time usually convolves other information (e.g., the varying vapor concentration and the sizedependent coagulation coefficient) in addition to particle growth. Figure S2 shows that with larger vapor molecules, the size range for the discrepancy between the theoretical and retrieved growth rate shifts towards larger diameters. Considering the influences of vapor evaporation and varying vapor concentration on the appearance time method, one should be cautious about the sizeresolved growth rate for particles close to the size of vapor molecules.
CoagS is assumed to be independent of time in the above discussions, whereas it may vary significantly during an NPF event in the atmosphere. The varying CoagS influences N_{i,∞} and hence the appearance time. Figure S3 shows a bias in the growth rate retrieved using the appearance time method caused by a varying CoagS.
In addition to vapor evaporation and the variation in vapor concentration, there may be other limitations for determining the appearance time in the atmosphere. Differently from controlled chamber studies (Dada et al., 2020), the uncertainties in atmospheric measurements pose challenges to growth rate estimation. These uncertainties come from instrumental biases; atmospheric turbulence; and the omitted contributions from transport, mixing, and emissions to the measured aerosol size distribution. Since the growth rate is calculated using a differential formula (Eq. 5), it is usually more sensitive to uncertainties than a physical quantity calculated using an integral formula (e.g., CoagS). For instance, the appearance time as a function of particle diameter in Fig. 7 had to be smoothed before calculating the growth rate using Eq. (5); otherwise, the calculated growth rate at some certain size bins would be negative. The applications of other methods to estimate the particle growth rate face the same challenge. As discussed in the Introduction, the appearance time method is used to estimate the new aerosol growth rate because other methods sometimes cannot report a growth rate in the concerned size range. Hence, further investigations concerning the uncertainties are needed for a better estimation of the growth rate in the atmosphere.
4.4 Application in atmospheric measurements
A typical intense NPF event measured in urban Beijing is used to test the correction for the influences of coagulation on the appearance time method. During the event, the peak sulfuric acid concentration is $\sim \mathrm{6}\times {\mathrm{10}}^{\mathrm{6}}$ cm^{−3} and the average CS for sulfuric acid is 0.024 s^{−1}. The theoretical condensation and coagulation growth rates are calculated using the measured sulfuric acid concentration and aerosol size distribution, respectively. Particle growth due to the uptake of sulfuric acid dimers is herein accounted for as the condensation growth. Note that the sum of the theoretical condensation and coagulation growth rate is not necessarily equal to the theoretical total growth rate for the measured NPF event. This is because only the condensation of sulfuric acid is considered although other vapors may also contribute to new particle growth. The enhancement due to Van der Waals forces is considered when calculating the coagulation coefficient (Alam, 1987; Chan and Mozurkewich, 2001; Stolzenburg et al., 2020). The appearance time retrieved from the measured aerosol size distributions was smoothed before estimating the particle growth rate. In addition to the example given in Fig. 8, the average growth rates of 1.5–3 nm particles measured in urban Beijing retrieved using the conventional and corrected appearance time methods were reported in Deng et al. (2020a).
As shown in Fig. 8, the measured H_{2}SO_{4} concentration started to increase around 07:00 (local time, UTC+8), and its 50 % appearance time was 07:30. Because SO_{2} concentration declined between 08:00 and 14:00, the peak concentration of H_{2}SO_{4} was observed at 08:00 rather than at noon. Shortly after the increase in H_{2}SO_{4} concentration, new particles down to the cluster size (∼ 1 nm in geometric diameter) were observed. The high concentration of sub2 nm aerosol during 07:00–15:00 and the continuous growth pattern of new particles from 1.5 to 10 nm indicates that the observed new particle formation event was a typical regional event (Kulmala et al., 2012), though transport, mixing, and emissions might have influenced the observed aerosol size distributions.
The impacts of particle coagulation are nonnegligible compared to particle growth and the growth rate calculation in urban Beijing. On one hand, the conventional appearance time method overestimates the particle growth rate for sub3 nm particles in urban Beijing due to the impact of CoagS. The deviation between the conventional and corrected growth rate decreases with the increasing diameter because CoagS decreases with particle diameter. As illustrated above, the correction for CoagSrc is only an approximation rather than obtained based on solid derivations. However, the negligible impact of CoagSrc on the measured growth rate in urban Beijing indicates that this approximation does not cause significant bias. Differently from coagulation growth which is weighted by particle size, the CoagSrc of d_{p} is only determined by the number concentrations of particles smaller than d_{p} (and their coagulation coefficient). Even under such an intense NPF event (with the maximum formation rate exceeding 200 cm^{−3} s^{−1}), the new particle concentration is usually much smaller than the vapor concentration due to the high CoagS and possibly cluster evaporation. Hence, it is sometimes acceptable to neglect the CoagSrc $/$ N_{p} term in Eqs. (6) and (8) to facilitate calculation. On the other hand, the coagulation with smaller particles enhances particle growth and this enhancement increases with the increasing particle size. This emphasizes that during an intensive NPF event with a high new particle concentration, the condensation growth rate contributed by condensing vapors cannot be taken as the total growth rate that determines the survival probability of new particles.
The difference between the measured and theoretical growth rates in Fig. 8 also indicates the growth mechanism of new particles. Considering the uncertainties in the appearance time, the sum of condensation and coagulation flux of sulfuric acid molecules and clusters is approximately equal to the measured particle growth rate for ∼ 3 nm particles, which indicates that sulfuric acid is a governing species that contributes to the initial growth of sub3 nm particles during the test event. The deviation between the measured growth and theoretical growth for particles larger than ∼ 3 nm indicates that there are other chemical species in addition to sulfuric acid (and the bases to neutralize it) contributing to particle growth. Note that the above discussion is only based on a single case study. Hence, further investigations based on longterm measurements are needed to reveal the growth mechanism in the polluted atmospheric environment.
Summarizing all the analysis above, the growth rate retrieved using the conventional appearance time method may be systematically overestimated due to the impact of coagulation, especially for intensive NPF events in polluted environments. Such an overestimation may be significant for sub3 nm particles because CoagS increases with decreasing particle size. In addition, the coagulation growth rate also needs to be corrected before investigating the condensation growth mechanism. For example, in the test case shown in Fig. 8, the retrieved condensation growth rate of 1.5 nm particles using the conventional appearance time method without correcting for the impact of CoagS and the coagulation growth rate is overestimated by 80 %. Figure 8 also indicates that the impact of CoagS may be negligible for larger particles and clean environments (see also Fig. S4). However, external sinks (e.g., dilution) may also cause an overestimation of the growth rate retrieved using the appearance time method if they are not properly corrected for.
The impact of coagulation on the particle growth rate retrieved using the appearance time method was investigated based on theoretical derivations and aerosol dynamics modeling. It was found that the oftenused 50 % sizeresolved appearance time method can reproduce the condensation growth rate only under the idealized condition without particle coagulation. When using the appearance time method in the real world, the coagulation sink, coagulation source, and coagulation growth need to be considered. Equations (5)–(9) provide a method in both discrete and sectional forms to correct for the impacts of the coagulation sink and coagulation source on the appearance time method. The feasibility of the correction for the influences of coagulation was verified using discrete and discretesectional aerosol models. In addition, vapor evaporation and the variation in vapor concentration was found to impact the appearance time method. The average uncertainties in the corrected 1.5–3 nm particle growth rate for each NPF event were 6 %–44 % in the test cases, and the maximum sizedependent uncertainty was 150 %. These uncertainties indicate that even after the correction for coagulation, one should be cautious about the appearance time method for particles close to the size of vapor molecules. Further, the growth rate of vapors and clusters is recommended to be estimated based on cluster dynamics instead of their representative time.
A typical NPF event measured in urban Beijing was used to show the quantitative impacts of coagulation on the retrieved growth rate. The systematic bias in the conventional appearance time method was observed for sub3 nm particles due to the uncorrected impact of the coagulation sink. Besides, coagulation growth was nonnegligible compared to the growth due to sulfuric acid condensation, which emphasizes the importance of distinguishing between the condensation and total growth rates. During the test event, the apparent growth rate of 1.5 nm particles retrieved using the conventional method was 80 % higher than the corrected condensation growth rate, whereas the corrected condensation growth rate was approximately equal to the theoretical growth rate contributed by sulfuric acid condensation.
Consider a particle population with a uniform diameter of d_{p} (nm) and a concentration of N_{0}. N_{0} is assumed to be sufficiently large so that the stochastics in particle growth are negligible. At the initial moment t_{0} (s), the mean particle diameter is d_{p}. During a short time interval dt (s), β_{1,p}N_{1}N_{0}dt particles collide with the condensing vapor with a diameter of d_{1}, where β_{1,p} is the coagulation coefficient (cm^{3} s^{−1}) and N_{1} is the vapor concentration. Hence, the mean diameter ($\stackrel{\mathrm{\u203e}}{{d}_{\text{p}}})$ weighted by the particle number concentration at the moment t_{0}+dt is
Comparing $\stackrel{\mathrm{\u203e}}{{d}_{\text{p}}}\left({t}_{\mathrm{0}}\right)$ and $\stackrel{\mathrm{\u203e}}{{d}_{\text{p}}}\left({t}_{\mathrm{0}}+\mathrm{d}t\right)$ yields the condensation growth rate:
The Taylor series of Eq. (A2) is
where $o\left[{\mathit{\beta}}_{\mathrm{1},\mathrm{p}}{N}_{\mathrm{1}}{d}_{\text{p}}{\left({d}_{\mathrm{1}}/{d}_{\mathrm{p}}\right)}^{\mathrm{6}}\right]$ is the Peano form of the remainder. The first term on the righthand side of Eq. (A3) is the formula for the particle growth rate in the continuous form, and the second term (the remainder) is the difference between the growth rate formula in the continuous and discrete forms (Olenius et al., 2018). When d_{p} is sufficiently larger than d_{1}, Eq. (A2) is reduced to ${\mathit{\beta}}_{\mathrm{1},\mathrm{p}}{N}_{\mathrm{1}}{d}_{\mathrm{1}}^{\mathrm{3}}/\left(\mathrm{3}{d}_{\mathrm{p}}^{\mathrm{2}}\right)$.
The Julia code for the discrete model is available upon request. The MATLAB code for the discretesectional model can be found via the link in Li and Cai (2020).
The measured data for Fig. 8 are available at https://doi.org/10.5281/zenodo.4506345 (Cai, 2021). The simulation results shown in Figs. 1–7 can be reproduced using the codes.
The supplement related to this article is available online at: https://doi.org/10.5194/acp2122872021supplement.
RC, JK, JJ, and MK initialized the study. RC and CL developed the models. RC and XCH derived the formulae. CD, YL, RY, CY, LW, and JJ performed the measurements and analyzed the data. RC performed the simulation and wrote the manuscript with the help of the other coauthors.
The authors declare that they have no conflict of interest.
This research has been supported by the Academy of Finland project (grant nos. 332547 and 1325656), a UHEL 3year grant (grant no. 75284132), the National Key R&D Program of China (grant no. 2017YFC0209503), and the Samsung PM_{2.5} SRP.
Openaccess funding was provided by the Helsinki University Library.
This paper was edited by Radovan Krejci and reviewed by Wolfgang Junkermann, Nikolaos Kalivitis, and one anonymous referee.
Alam, M. K.: The Effect of van der Waals and Viscous Forces on Aerosol Coagulation, Aerosol Sci. Technol., 6, 41–52, https://doi.org/10.1080/02786828708959118, 1987.
Anttila, T., Kerminen, V.M., and Lehtinen, K. E. J.: Parameterizing the formation rate of new particles: The effect of nuclei selfcoagulation, J. Aerosol Sci., 41, 621–636, https://doi.org/10.1016/j.jaerosci.2010.04.008, 2010.
Cai, R.: Impacts of coagulation on the appearance time method for new particle growth rate evaluation and their corrections [Data set], Zenodo, https://doi.org/10.5281/zenodo.4506345, 2021.
Cai, R. and Jiang, J.: A new balance formula to estimate new particle formation rate: reevaluating the effect of coagulation scavenging, Atmos. Chem. Phys., 17, 12659–12675, https://doi.org/10.5194/acp17126592017, 2017.
Cai, R., Chen, D.R., Hao, J., and Jiang, J.: A miniature cylindrical dierential mobility analyzer for sub3nm particle sizing, J. Aerosol Sci., 106, 111–119, https://doi.org/10.1016/j.jaerosci.2017.01.004, 2017.
Cai, R., Chandra, I., Yang, D., Yao, L., Fu, Y., Li, X., Lu, Y., Luo, L., Hao, J., Ma, Y., Wang, L., Zheng, J., Seto, T., and Jiang, J.: Estimating the influence of transport on aerosol size distributions during new particle formation events, Atmos. Chem. Phys., 18, 16587–16599, https://doi.org/10.5194/acp18165872018, 2018.
Chan, T. W. and Mozurkewich, M.: Measurement of the coagulation rate constant for sulfuric acid particles as a function of particle size using tandem differential mobility analysis, J. Aerosol Sci., 32, 321–339, https://doi.org/10.1016/S00218502(00)000811, 2001.
Dada, L., Lehtipalo, K., Kontkanen, J., Nieminen, T., Baalbaki, R., Ahonen, L., Duplissy, J., Yan, C., Chu, B., Petaja, T., Lehtinen, K., Kerminen, V. M., Kulmala, M., and Kangasluoma, J.: Formation and growth of sub3nm aerosol particles in experimental chambers, Nat. Protoc., 15, 1013–1040, https://doi.org/10.1038/s415960190274z, 2020.
Deng, C., Cai, R., Yan, C., Zheng, J., and Jiang, J.: Formation and growth of sub3nm particles in megacities: impacts of background aerosols, Faraday Discuss., Advance Article, https://doi.org/10.1039/d0fd00083c, 2020a.
Deng, C., Fu, Y., Dada, L., Yan, C., Cai, R., Yang, D., Zhou, Y., Yin, R., Lu, Y., Li, X., Qiao, X., Fan, X., Nie, W., Kontkanen, J., Kangasluoma, J., Chu, B., Ding, A., Kerminen, V., Paasonen, P., Worsnop, R. D., Bianchi, F., Liu, Y., Zheng, J., Wang, L., Kulmala, M., and Jiang, J.: Seasonal characteristics of new particle formation and growth in urban Beijing, Environ. Sci. Technol., 54, 8547–8557, https://doi.org/10.1021/acs.est.0c00808, 2020b.
Ehn, M., Thornton, J. A., Kleist, E., Sipila, M., Junninen, H., Pullinen, I., Springer, M., Rubach, F., Tillmann, R., Lee, B., LopezHilfiker, F., Andres, S., Acir, I. H., Rissanen, M., Jokinen, T., Schobesberger, S., Kangasluoma, J., Kontkanen, J., Nieminen, T., Kurten, T., Nielsen, L. B., Jorgensen, S., Kjaergaard, H. G., Canagaratna, M., Maso, M. D., Berndt, T., Petaja, T., Wahner, A., Kerminen, V. M., Kulmala, M., Worsnop, D. R., Wildt, J., and Mentel, T. F.: A large source of lowvolatility secondary organic aerosol, Nature, 506, 476–479, https://doi.org/10.1038/nature13032, 2014.
Fu, Y., Xue, M., Cai, R., Kangasluoma, J., and Jiang, J.: Theoretical and experimental analysis of the core sampling method: Reducing diffusional losses in aerosol sampling line, Aerosol Sci. Technol., 53, 793–801, https://doi.org/10.1080/02786826.2019.1608354, 2019.
He, X.C., Iyer, S., Sipilä, M., Ylisirniö, A., Peltola, M., Kontkanen, J., Baalbaki, R., Simon, M., Kürten, A., Tham, Y. J., Pesonen, J., Ahonen, L. R., Amanatidis, S., Amorim, A., Baccarini, A., Beck, L., Bianchi, F., Brilke, S., Chen, D., Chiu, R., Curtius, J., Dada, L., Dias, A., Dommen, J., Donahue, N. M., Duplissy, J., El Haddad, I., Finkenzeller, H., Fischer, L., Heinritzi, M., Hofbauer, V., Kangasluoma, J., Kim, C., Koenig, T. K., Kubečka, J., Kvashnin, A., Lamkaddam, H., Lee, C. P., Leiminger, M., Li, Z., Makhmutov, V., Xiao, M., Marten, R., Nie, W., Onnela, A., Partoll, E., Petäjä, T., Salo, V.T., Schuchmann, S., Steiner, G., Stolzenburg, D., Stozhkov, Y., Tauber, C., Tomé, A., Väisänen, O., VazquezPufleau, M., Volkamer, R., Wagner, A. C., Wang, M., Wang, Y., Wimmer, D., Winkler, P. M., Worsnop, D. R., Wu, Y., Yan, C., Ye, Q., Lehtinen, K., Nieminen, T., Manninen, H. E., Rissanen, M., Schobesberger, S., Lehtipalo, K., Baltensperger, U., Hansel, A., Kerminen, V.M., Flagan, R. C., Kirkby, J., Kurtén, T., and Kulmala, M.: Determination of the Collision Rate Coefficient between Charged Iodic Acid Clusters and Iodic Acid using the Appearance Time Method, Aerosol Sci. Technol., 1–17, https://doi.org/10.1080/02786826.2020.1839013, 2020.
Jiang, J., Chen, M., Kuang, C., Attoui, M., and McMurry, P. H.: Electrical Mobility Spectrometer Using a Diethylene Glycol Condensation Particle Counter for Measurement of Aerosol Size Distributions Down to 1 nm, Aerosol Sci. Technol., 45, 510–521, https://doi.org/10.1080/02786826.2010.547538, 2011a.
Jiang, J., Zhao, J., Chen, M., Eisele, F. L., Scheckman, J., Williams, B. J., Kuang, C., and McMurry, P. H.: First Measurements of Neutral Atmospheric Cluster and 1–2nm Particle Number Size Distributions During Nucleation Events, Aerosol Sci. Technol., 45, ii–v, https://doi.org/10.1080/02786826.2010.546817, 2011b.
Jokinen, T., Sipilä, M., Junninen, H., Ehn, M., Lönn, G., Hakala, J., Petäjä, T., Mauldin III, R. L., Kulmala, M., and Worsnop, D. R.: Atmospheric sulphuric acid and neutral cluster measurements using CIAPiTOF, Atmos. Chem. Phys., 12, 4117–4125, https://doi.org/10.5194/acp1241172012, 2012.
Kangasluoma, J., Cai, R., Jiang, J., Deng, C., Stolzenburg, D., Ahonen, L. R., Chan, T., Fu, Y., Kim, C., Laurila, T. M., Zhou, Y., Dada, L., Sulo, J., Flagan, R. C., Kulmala, M., Petäjä, T., and Lehtinen, K.: Overview of measurements and current instrumentation for 110nm aerosol particle number size distributions, J. Aerosol Sci., 148, 105584, https://doi.org/10.1016/j.jaerosci.2020.105584, 2020.
Kerminen, V.M., Paramonov, M., Anttila, T., Riipinen, I., Fountoukis, C., Korhonen, H., Asmi, E., Laakso, L., Lihavainen, H., Swietlicki, E., Svenningsson, B., Asmi, A., Pandis, S. N., Kulmala, M., and Petäjä, T.: Cloud condensation nuclei production associated with atmospheric nucleation: a synthesis based on existing literature and new results, Atmospheric Chemistry and Physics, 12, 1203712059, 10.5194/acp12120372012, 2012.
Kerminen, V.M., Chen, X., Vakkari, V., Petäjä, T., Kulmala, M., and Bianchi, F.: Atmospheric new particle formation and growth: review of field observations, Environ. Res. Lett., 13, https://doi.org/10.1088/17489326/aadf3c, 2018.
Kontkanen, J., Olenius, T., Lehtipalo, K., Vehkamäki, H., Kulmala, M., and Lehtinen, K. E. J.: Growth of atmospheric clusters involving cluster–cluster collisions: comparison of different growth rate methods, Atmos. Chem. Phys., 16, 5545–5560, https://doi.org/10.5194/acp1655452016, 2016.
Kontkanen, J., Deng, C., Fu, Y., Dada, L., Zhou, Y., Cai, J., Daellenbach, K. R., Hakala, S., Kokkonen, T. V., Lin, Z., Liu, Y., Wang, Y., Yan, C., Petäjä, T., Jiang, J., Kulmala, M., and Paasonen, P.: Sizeresolved particle number emissions in Beijing determined from measured particle size distributions, Atmos. Chem. Phys., 20, 11329–11348, https://doi.org/10.5194/acp20113292020, 2020.
Kuang, C., McMurry, P. H., and McCormick, A. V.: Determination of cloud condensation nuclei production from measured new particle formation events, Geophys. Res. Lett., 36, L09822, https://doi.org/10.1029/2009gl037584, 2009.
Kuang, C., Riipinen, I., Sihto, S.L., Kulmala, M., McCormick, A. V., and McMurry, P. H.: An improved criterion for new particle formation in diverse atmospheric environments, Atmos. Chem. Phys., 10, 8469–8480, https://doi.org/10.5194/acp1084692010, 2010.
Kuang, C., Chen, M., Zhao, J., Smith, J., McMurry, P. H., and Wang, J.: Size and timeresolved growth rate measurements of 1 to 5 nm freshly formed atmospheric nuclei, Atmos. Chem. Phys., 12, 3573–3589, https://doi.org/10.5194/acp1235732012, 2012.
Kulmala, M., Dal Maso, M., Mäkelä, J. M., Pirjola, L., Väkevä, M., Aalto, P., Miikkulainen, P., Hämeri, K., and O'dowd, C. D.: On the formation, growth and composition of nucleation mode particles, Tellus, 53, 479–490, https://doi.org/10.3402/tellusb.v53i4.16622, 2001.
Kulmala, M., Vehkamäki, H., Petäjä, T., Dal Maso, M., Lauri, A., Kerminen, V.M., Birmili, W., and McMurry, P. H.: Formation and growth rates of ultrafine atmospheric particles: a review of observations, J. Aerosol Sci., 35, 143–176, https://doi.org/10.1016/j.jaerosci.2003.10.003, 2004.
Kulmala, M., Petäjä, T., Nieminen, T., Sipilä, M., Manninen, H. E., Lehtipalo, K., Dal Maso, M., Aalto, P. P., Junninen, H., Paasonen, P., Riipinen, I., Lehtinen, K. E., Laaksonen, A., and Kerminen, V.M.: Measurement of the nucleation of atmospheric aerosol particles, Nat. Protoc., 7, 1651–1667, https://doi.org/10.1038/nprot.2012.091, 2012.
Lee, S.H., Gordon, H., Yu, H., Lehtipalo, K., Haley, R., Li, Y., and Zhang, R.: New Particle Formation in the Atmosphere: From Molecular Clusters to Global Climate, J. Geophys. Res.Atmos., 124, 7098–7146, https://doi.org/10.1029/2018jd029356, 2019.
Lehtinen, K. E. J. and Kulmala, M.: A model for particle formation and growth in the atmosphere with molecular resolution in size, Atmos. Chem. Phys., 3, 251–257, https://doi.org/10.5194/acp32512003, 2003.
Lehtinen, K. E. J., Dal Maso, M., Kulmala, M., and Kerminen, V.M.: Estimating nucleation rates from apparent particle formation rates and vice versa: Revised formulation of the KerminenKulmala equation, J. Aerosol Sci., 38, 988–994, https://doi.org/10.1016/j.jaerosci.2007.06.009, 2007.
Lehtipalo, K., Leppä, J., Kontkanen, J., Kangasluoma, J., Franchin, A., Wimmer, D., Schobesberger, S., Junninen, H., Petäjä, T., Sipilä, M., Mikkilä, J., Vanhanen, J., Worsnop, D. R., and Kulmala, M.: Methods for determining particle size distribution and growth rates between 1 and 3 nm using the Particle Size Magnifier, Boreal Environ. Res., 19, 215–236, 2014.
Li, C. and Cai, R.: Tutorial: The discretesectional method to simulate an evolving aerosol, J. Aerosol Sci., 150, 105615, https://doi.org/10.1016/j.jaerosci.2020.105615, 2020.
Li, C. and McMurry, P. H.: Errors in nanoparticle growth rates inferred from measurements in chemically reacting aerosol systems, Atmos. Chem. Phys., 18, 8979–8993, https://doi.org/10.5194/acp1889792018, 2018.
Liu, J., Jiang, J., Zhang, Q., Deng, J., and Hao, J.: A spectrometer for measuring particle size distributions in the range of 3 nm to 10 μm, Front. Env. Sci. Eng., 10, 63–72, https://doi.org/10.1007/s117830140754x, 2016.
McMurry, P. H.: New particle formation in the presence of an aerosol: rates, time scales, and sub0.01 μm size distributions, J. Colloid Interf. Sci., 95, 72–80, https://doi.org/10.1016/00219797(83)900735, 1983.
Mohr, C., Thornton, J. A., Heitto, A., LopezHilfiker, F. D., Lutz, A., Riipinen, I., Hong, J., Donahue, N. M., Hallquist, M., Petaja, T., Kulmala, M., and YliJuuti, T.: Molecular identification of organic vapors driving atmospheric nanoparticle growth, Nat. Commun., 10, 4442, https://doi.org/10.1038/s41467019124732, 2019.
Myllys, N., Kubečka, J., Besel, V., Alfaouri, D., Olenius, T., Smith, J. N., and Passananti, M.: Role of base strength, cluster structure and charge in sulfuricaciddriven particle formation, Atmos. Chem. Phys., 19, 9753–9768, https://doi.org/10.5194/acp1997532019, 2019.
Nieminen, T., Kerminen, V.M., Petäjä, T., Aalto, P. P., Arshinov, M., Asmi, E., Baltensperger, U., Beddows, D. C. S., Beukes, J. P., Collins, D., Ding, A., Harrison, R. M., Henzing, B., Hooda, R., Hu, M., Hõrrak, U., Kivekäs, N., Komsaare, K., Krejci, R., Kristensson, A., Laakso, L., Laaksonen, A., Leaitch, W. R., Lihavainen, H., Mihalopoulos, N., Németh, Z., Nie, W., O'Dowd, C., Salma, I., Sellegri, K., Svenningsson, B., Swietlicki, E., Tunved, P., Ulevicius, V., Vakkari, V., Vana, M., Wiedensohler, A., Wu, Z., Virtanen, A., and Kulmala, M.: Global analysis of continental boundary layer new particle formation based on longterm measurements, Atmos. Chem. Phys., 18, 14737–14756, https://doi.org/10.5194/acp18147372018, 2018.
Olenius, T., Riipinen, I., Lehtipalo, K., and Vehkamäki, H.: Growth rates of atmospheric molecular clusters based on appearance times and collisionevaporation fluxes: Growth by monomers, J. Aerosol Sci., 78, 55–70, https://doi.org/10.1016/j.jaerosci.2014.08.008, 2014.
Olenius, T., Pichelstorfer, L., Stolzenburg, D., Winkler, P. M., Lehtinen, K. E. J., and Riipinen, I.: Robust metric for quantifying the importance of stochastic effects on nanoparticle growth, Sci. Rep., 8, 14160, https://doi.org/10.1038/s4159801832610z, 2018.
Pichelstorfer, L., Stolzenburg, D., Ortega, J., Karl, T., Kokkola, H., Laakso, A., Lehtinen, K. E. J., Smith, J. N., McMurry, P. H., and Winkler, P. M.: Resolving nanoparticle growth mechanisms from size and timedependent growth rate analysis, Atmos. Chem. Phys., 18, 1307–1323, https://doi.org/10.5194/acp1813072018, 2018.
Stolzenburg, D., Simon, M., Ranjithkumar, A., Kürten, A., Lehtipalo, K., Gordon, H., Ehrhart, S., Finkenzeller, H., Pichelstorfer, L., Nieminen, T., He, X.C., Brilke, S., Xiao, M., Amorim, A., Baalbaki, R., Baccarini, A., Beck, L., Bräkling, S., Caudillo Murillo, L., Chen, D., Chu, B., Dada, L., Dias, A., Dommen, J., Duplissy, J., El Haddad, I., Fischer, L., Gonzalez Carracedo, L., Heinritzi, M., Kim, C., Koenig, T. K., Kong, W., Lamkaddam, H., Lee, C. P., Leiminger, M., Li, Z., Makhmutov, V., Manninen, H. E., Marie, G., Marten, R., Müller, T., Nie, W., Partoll, E., Petäjä, T., Pfeifer, J., Philippov, M., Rissanen, M. P., Rörup, B., Schobesberger, S., Schuchmann, S., Shen, J., Sipilä, M., Steiner, G., Stozhkov, Y., Tauber, C., Tham, Y. J., Tomé, A., VazquezPufleau, M., Wagner, A. C., Wang, M., Wang, Y., Weber, S. K., Wimmer, D., Wlasits, P. J., Wu, Y., Ye, Q., ZaunerWieczorek, M., Baltensperger, U., Carslaw, K. S., Curtius, J., Donahue, N. M., Flagan, R. C., Hansel, A., Kulmala, M., Lelieveld, J., Volkamer, R., Kirkby, J., and Winkler, P. M.: Enhanced growth rate of atmospheric particles from sulfuric acid, Atmos. Chem. Phys., 20, 7359–7372, https://doi.org/10.5194/acp2073592020, 2020.
Stolzenburg, M. R., McMurry, P. H., Sakurai, H., Smith, J. N., Mauldin, R. L., Eisele, F. L., and Clement, C. F.: Growth rates of freshly nucleated atmospheric particles in Atlanta, J. Geophys. Res., 110, D22S05, https://doi.org/10.1029/2005jd005935, 2005.
Weber, R. J., Marti, J. J., McMurry, P. H., Eisele, F. L., Tanner, D. J., and Jefferson, A.: Measurements of new particle formation and ultrafine particle growth rates at a clean continental site, J. Geophys. Res.Atmos., 102, 4375–4385, https://doi.org/10.1029/96jd03656, 1997.
Wiedensohler, A., Birmili, W., Nowak, A., Sonntag, A., Weinhold, K., Merkel, M., Wehner, B., Tuch, T., Pfeifer, S., Fiebig, M., Fjäraa, A. M., Asmi, E., Sellegri, K., Depuy, R., Venzac, H., Villani, P., Laj, P., Aalto, P., Ogren, J. A., Swietlicki, E., Williams, P., Roldin, P., Quincey, P., Hüglin, C., FierzSchmidhauser, R., Gysel, M., Weingartner, E., Riccobono, F., Santos, S., Grüning, C., Faloon, K., Beddows, D., Harrison, R., Monahan, C., Jennings, S. G., O'Dowd, C. D., Marinoni, A., Horn, H.G., Keck, L., Jiang, J., Scheckman, J., McMurry, P. H., Deng, Z., Zhao, C. S., Moerman, M., Henzing, B., de Leeuw, G., Löschau, G., and Bastian, S.: Mobility particle size spectrometers: harmonization of technical standards and data structure to facilitate high quality longterm observations of atmospheric particle number size distributions, Atmos. Meas. Tech., 5, 657–685, https://doi.org/10.5194/amt56572012, 2012.
Yao, L., Garmash, O., Bianchi, F., Zheng, J., Yan, C., Kontkanen, J., Junninen, H., Mazon, S. B., Ehn, M., Paasonen, P., Sipilä, M., Wang, M., Wang, X., Xiao, S., Chen, H., Lu, Y., Zhang, B., Wang, D., Fu, Q., Geng, F., Li, L., Wang, H., Qiao, L., Yang, X., Chen, J., Kerminen, V.M., Petäjä, T., Worsnop, D. R., Kulmala, M., and Wang, L.: Atmospheric new particle formation from sulfuric acid and amines in a Chinese megacity, Science, 361, 278–281, https://doi.org/10.1126/science.aao4839, 2018.