Interactive comment on “ Simulating secondary organic aerosol from anthropogenic and biogenic precursors : comparison to outdoor chamber experiments , effect of oligomerization on SOA formation and reactive uptake of aldehydes

lines 6-7: replace "were" with "where" p.3, line 20: It is recommended that the authors be more specific about the ranges of conditions covered by the experiments, particularly for atmospherically-relevant conditions that are outside the available datasets. p.5, line 3-4: Chamber data studies are referenced twice. p. 7, lines 27-30: How are "short", "big", and "bigger" oligomers defined?

Abstract.New parameterizations for the formation of organic aerosols were developed.These parameterizations cover secondary organic aerosol (SOA) formation from biogenic and anthropogenic precursors, NO x dependency, oligomerization and the reactive uptake of pinonaldehyde.These parameterizations were implemented in a box model in which the condensation and/or evaporation of semi-volatile organic compounds was simulated by the Secondary Organic Aerosol Processor (SOAP) model to take the dynamic evolution of concentrations into account.
The parameterizations were tested against several experiments carried out in previous studies in the EUPHORE outdoor chamber.Two datasets of experiments were used: the anthropogenic experiments (in which SOA is formed mainly from a mixture of toluene, 1,3,5-trimethylbenzene and o-xylene) and the biogenic experiments (in which SOA is formed mainly from α-pinene and limonene).
When assuming no wall deposition of organic vapors, satisfactory results (bias lower than 20 %) were obtained for the biogenic experiments and for most of the anthropogenic experiments.However, a decrease of SOA concentrations (up to 30 %) was found when taking wall deposition of organic vapors into account (with the parameters of Zhang et al., 2014).The anthropogenic experiments seem to indicate a complex NO x dependency that could not be reproduced by the model.Oligomerization was found to have a strong effect on SOA composition (oligomers were estimated to account for up to 78 % of the SOA mass) and could therefore have a strong effect on the formation of SOA.The uptake of pinonaldehyde (which is a high-volatility semi-volatile organic compound, SVOC) onto acidic aerosol was found to be too slow to be significant under atmospheric conditions (no significant amount of SOA formed after 3 days of evolution), indicating that the parameterization of Pun and Seigneur (2007) used in some air quality models may lead to an overestimation of SOA concentrations.The uptake of aldehydes could nevertheless be an important SOA formation pathway for less volatile or more reactive aldehydes than pinonaldehyde.
Regarding viscosity, a low effect of viscosity on SOA concentrations was estimated by the model, although a decrease of SVOC evaporation was found when taking it into account, as well as a lower sensitivity of concentrations to changes of temperature during the experiments.

Introduction
Because of the effect of fine particles on human health (WHO, 2003) and ecosystems (Kanakidou et al., 2005), the use of models has become a common practice to evaluate impacts and mitigation strategies.Particulate organic matter (OM) represents a large fraction of the total fine particulate mass, typically between 20 % and 60 % (Kanakidou et al., 2005;Yu et al., 2007;Zhang et al., 2007), with the secondary fraction (secondary organic aerosol, SOA) representing most of it (90 % according to the best estimate of Kanakidou et al., 2005).Therefore, efforts have to be made to represent OM as accurately as possible in models.
Numerous models have been developed to simulate OM in 3-D air quality models (Schell et al., 2001;Donahue et al., 2006Donahue et al., , 2011;;Pun et al., 2002;Couvidat et al., 2012;Griffin et al., 2003;Jathar et al., 2015;Tulet et al., 2006;Carlton et al., 2010;Menut et al., 2013).Most of these models use simple parameterizations based on SOA yields estimated from smog chamber experiments conducted under specific conditions, which can be different from atmospheric conditions (low humidity, specific NO x conditions).Nevertheless, SOA yields from chamber experiments have traditionally been estimated considering particle wall losses, but not gas wall losses.Several studies (Zhang et al., 2014(Zhang et al., , 2015;;Bian et al., 2015;Cappa et al., 2016) indicate that this fact could have led to a strong underestimation of SOA yields.On the other side, none of these models take the whole complexity of the processes involved in organic aerosol formation (nonideality, multiphase partitioning, viscosity of the aerosol, phase separation, aging, oligomerization and organosulfate formation, effects of NO x concentrations, etc.) into account, which can highly affect the level of SOA predicted (e.g., Ng et al., 2007a, b;Pun and Seigneur, 2007;Couvidat and Seigneur, 2011;Hall IV and Johnston, 2011).Therefore, the development of parameterizations addressing these aspects can provide insights about the SOA formation processes and improve current model estimates.
Oligomer formation has been addressed in some modeling studies.This process may be important for organic aerosol as it can transform semi-volatile and volatile organic compounds into less volatile compounds.Trump and Donahue (2014) studied the effect of oligomerization on the dynamic of organic aerosol formation.Pun and Seigneur (2007) developed a parameterization for the oligomerization of aldehydes by increasing their partitioning toward the particle.This parameterization (treating oligomerization with an equilibrium constant) has been used in Couvidat et al. (2012) to increase SOA formation from pinonaldehyde.High aerosol concentrations were simulated with this parameterization.However, Liggio and Li (2006b) showed that the uptake of pinonaldehyde onto acidic aerosol is a slow process and that this process could be due to oligomerization but also to organosulfate formation.Carlton et al. (2010) used a simple first-order rate constant of oligomerization for all organic compounds based on the results of Kalberer et al. (2004).Jathar et al. (2016) used the parameterization of Carlton et al. (2010) and showed that, depending on the mechanism, it could either have a strong effect on SOA formation or not.In all the cases, the authors found a strong effect on composition (or volatility distribution).Lemaire et al. (2016) compared these different approaches for oligomerization and emphasized the need to simulate oligomerization properly in air quality models.
Several approaches have been used to represent SOA formation in 3-D air quality models.Most of these approaches have in common that they are based on results from smog chamber experiments but differ in the methodology for lumping organic compounds and the treatment of the processes involved in SOA formation.Among these approaches, the volatility basis set (VBS) approach (Donahue et al., 2006;Robinson et al., 2007) has been widely used (Lane et al., 2008;Shrivastava et al., 2008;Murphy and Pandis, 2009;Shrivastava et al., 2011;Han et al., 2016).It uses a basis set in which semi-volatile organic compounds (SVOCs) are distributed on an equally spaced logarithmic scale for volatility (volatility bins).Aging was treated in the VBS by assuming a shift of volatility per aging step (gas-phase oxidation creates a change of volatility bins).A 1.5-D VBS (Koo et al., 2014) and a 2-D VBS (Donahue et al., 2011) were also developed to take changes in the oxidation degree into account.
Another approach used in 3-D air quality models is the molecular surrogate approach (Pun et al., 2002(Pun et al., , 2006;;Griffin et al., 2003;Couvidat et al., 2012).In this approach, surrogate compounds (that gather a large number of SVOCs with similar thermodynamic properties) are associated with molecular structures to extrapolate SOA formation from smog chambers to the atmosphere.In the molecular surrogate approach, several processes, which are often not taken into account in the other approaches, can be readily estimated (e.g., absorption into an aqueous phase, hygroscopicity, nonideality) and can be treated explicitly in the model.Based on this methodology, the Secondary Organic Aerosol Processor (SOAP) (Couvidat and Sartelet, 2015) was developed to simulate the gas-particle partitioning of SVOCs by taking nonideality (by considering the interactions between organic and inorganic compounds), multiphase partitioning, phase separation, hygroscopicity and the viscous state of organic aerosols into account.It was recently implemented in the CHIMERE air quality model (Couvidat et al., 2018) to simulate the gas-particle partitioning of SVOC formed with the H 2 O mechanism (Couvidat et al., 2012).However, in H 2 O, a simplified aging mechanism was only applied for the first step of oxidation of anthropogenic primary SVOC due to lack of information, and SVOCs were supposed to be only hydrophilic (condense only the aqueous phase of particles) or only hydrophobic (condense only the organic phase of particles).Moreover, for anthropogenic SOA, comparison to measurements carried out inside EUPHORE has shown that for anthropogenic SOA the H 2 O mechanism may overestimate SOA in the conditions of the experiments.
In order to create a SOA mechanism that could be implemented in 3-D air quality models and that considers the main phenomena involved in SOA formation (nonideality, multiphase partitioning, viscosity of the aerosol, phase separation, aging, oligomerization, effects of NO x concentrations), the aim of this study is to improve and update the H 2 O mechanism for the formation of SOA from the oxidation of monoterpenes and aromatics by reconsidering some of the assumptions made in the development of this mechanism (no aging, compounds only hydrophilic or only hydrophobic, use of the parameterization of Pun and Seigneur, 2007) and have a more realistic representation of SOA formation.The update of the mechanism was done based on SOA yields available in the literature with the methodology of Odum et al. (1996); an aging mechanism was added and SVOC surrogates are assumed to be able to condense on both the aqueous phase and the organic phase.In order to evaluate this SOA mechanism, it was tested against some experiments carried out in the EUPHORE outdoor chamber under conditions different from the conditions under which it was developed (mixture of different compounds, noncontrolled temperature) to ensure the good performances of the model.
Two new parameterizations on particle-phase reactions were also developed covering oligomerization and the uptake of pinonaldehyde.For oligomerization, instead of representing this process with a first-order complete reaction as done by Carlton et al. (2010), it was chosen to reproduce this phenomenon (as an effort to provide a more realistic representation) with a second-order reversible reaction unfavored by humid conditions by analogy to oligomerization reactions expected to occur in the particle like esterification, hemiacetalization, aldolization and peroxyhemiacetalization (Roldin et al., 2014).For the uptake of pinonaldehyde, instead of assuming equilibrium (as done by Pun and Seigneur, 2007), a parameterization for the dynamic of the uptake, which was shown to limit strongly the uptake of pinonaldehyde (Liggio and Li, 2006b) as a function of the particle acidity, was taken into account.
In addition, the effect of particle-phase reactions was investigated.Oligomerization was simulated with the secondorder reversible reaction to evaluate the impact of this process and how it may impact SOA formation.For the uptake of pinonaldehyde onto acidic aerosols, the results obtained with the equilibrium-based parameterization of Pun and Seigneur (2007) were also compared to the results given by the dynamic parameterization developed in this study.SOA concentrations simulated with these parameterizations were compared to measured concentrations from experiments in the presence of SO 2 (that leads to the formation of acidic aerosols).Finally, the effects of low diffusion inside the particle due to high viscosity and of wall losses of organic vapors were also studied.

Method
The new parameterizations described hereafter were developed and tested against experiments that were carried out in previous studies inside the outdoor chamber EUPHORE in Valencia, Spain.

Experimental datasets
EUPHORE is a 204 m 3 hemispherical Teflon outdoor chamber located in Valencia (CEAM), Spain.The chamber is surrounded by retractable steel housing, which can be opened or closed to control the time of exposure to sunlight.The housing also serves to protect the chamber from rain and strong winds.The floor of the reactor consists of aluminum panels covered with Teflon, with a cooling system designed to compensate the heating of the chamber caused by solar irradiation.Two high-powered mixing fans (with an air flow of 4000 m 3 h −1 ), located on the floor of the chamber, are used to ensure the homogeneous mixing of reactants and products.The EUPHORE facility is described in detail in Volkamer et al. (2001).
The experiments used for the comparison between the model and experiments have been published in previous studies (Vivanco et al., 2011(Vivanco et al., , 2013)).The experiments used in this study were initially planned to increase the experimental information regarding SOA formation and to be mainly used in model performance evaluation processes.As most of the studies at that time were focused on the oxidation of individual organic gases or simple mixtures of them, those experiments constitute a contribution to the experimental database on SOA formation, by considering the oxidation of different mixtures of organic precursors.Although the experimental dataset was designed for model evaluation purposes, it was not designed to evaluate the model developed in this study.Vivanco et al. (2016) and Santiago et al. (2012b) used these experiments to evaluate simple parameterizations existing in several air quality models.The experiments may not cover the full range of values covered by the model for some parameters, especially for relative humidity (RH), which was lower than 40 % for all experiments.Experimental conditions are described in Tables 1 and 2.
OH was formed by HONO photolysis.Ozone was formed throughout the experiment and often reached concentrations around 200 ppb.However, at the start of the experiments, in which concentrations of ozone are lower, monoterpenes will mainly react with OH radicals.As almost all the monoterpenes reacted during the first hour of the reaction, monoterpenes can be considered to react both with OH and O 3 .Calculations carried with RACM2 showed that 10 % to 70 % of the monoterpenes should have reacted with O 3 .For experiments B4 and B5 with very low relative humidities (around 10 %), most of the monoterpenes reacted with OH (around 90 %).OH photo-oxidation time (determined for a concen-F.Couvidat et al.: Simulating SOA from anthropogenic and biogenic precursors tration of 1.5×10 6 molecules cm −3 ) was determined to range from a few hours to more than a day.Two datasets of experiments are used: (i) the anthropogenic experiments where SOA is formed from the oxidation of toluene (TOL), 1,3,5-trimethylbenzene (TMB) and to a lesser extent o-xylene and octane and (ii) the biogenic experiments where SOA is formed from α-pinene (API) and limonene (LIM) and to a lesser extent isoprene.Experimental conditions for the anthropogenic and biogenic experiments are described, respectively, in Tables 1 and 2. For the biogenic dataset, two experiments were carried out in the presence of SO 2 to evaluate the parameterization of Pun and Seigneur (2007).Due to its low SOA yield, SOA from isoprene oxidation should not represent a significant amount of total SOA.
Although octane is a precursor of SOA and is present in the anthropogenic experiments, SOA formation from octane oxidation was not taken into account in this study as Vivanco et al. (2016) showed with an experiment in the same chamber that octane leads to an insignificant amount of SOA.Due to its low yield (0.5 % according to Lim and Ziemann, 2005), octane SOA should not represent more than a few percent of total SOA.PM volume concentrations were measured with a Scanning Mobility Particle Sizer Spectrometer (SMPS).Data are not corrected for wall losses as the deposition of particles is taken into account in the model.

Model development
Secondary organic aerosols inside the chamber were simulated by coupling the gas-phase mechanism RACM2 (Goliff et al., 2013) with the SOAP model (Couvidat and Sartelet, 2015) to compute the dynamic formation of SOA.RACM2 was used because it has been shown to perform well for oxidant formation (Kim et al., 2009).To represent the chemical evolution of SVOC, RACM2 was modified to take the formation of the surrogate species into account according to the mechanisms described hereafter.The ROS2 algorithm (Verwer et al., 1999) was used to solve the chemical kinetic equations.
SOAP is a model designed to be modular with user options depending on the computation time.SOAP uses the molecular surrogate approach to estimate several properties and parameters (hygroscopicity, absorption into the aqueous phase of particles, activity coefficients and phase separation) and to evaluate the partitioning of organic compounds between one or several organic phases (the number of organic phases is determined by the Gibbs energy minimization) and the inorganic phase.It accounts for the influence of interactions between organic and inorganic compounds by using the AIOMFAC algorithm (Zuend et al., 2008(Zuend et al., , 2011;;Zuend and Seinfeld, 2012).Secondary inorganic aerosol formation was added to the SOAP model by using the equilibrium parameters of ISORROPIA v2.1 (Fountoukis and Nenes, 2007).
Currently, SOAP assumes that inorganic aerosols are metastable liquids and therefore does not take efflorescence or deliquescence processes into account.This assumption could be wrong in the presence of ammonia due to the low humidity inside the chamber, as ammonium sulfate would probably be solid as the humidity in the chamber (below 40 %) is always far below the deliquescence relative humidity (80 % at 298 K).However, no ammonia was present in these experiments and SO 2 was only introduced for two experiments.For those two experiments, SO 2 oxidation will lead to sulfuric acid formation which remains liquid over the full RH range (Seinfeld and Pandis, 1998).
Several experimental data were used to constrain the model.Temperature and relative humidity measurements inside the chamber were used as inputs for the model.SMPS measurements were used to compute the mean diameter of particles.As modeling properly nucleation and coagulation of particles would be needed to simulate the size distribution of particles adequately, particles were gathered inside a single diameter bin.The mean diameter of particles was constrained to provide the model with a realistic estimation of the diameter without modeling the nucleation in SOAP.The diameter of the particle is used to compute the kinetic rate of condensation and/or evaporation and/or diffusion and the Kelvin effect by the dynamic approach of SOAP.
The model results were compared to non-corrected (according to wall losses) results.Wall deposition of particles was simulated via the first-order loss rate parameter constrained to reproduce the loss of particles during the last hours of the experiment when the chamber was enclosed by the retractable steel housing.
Wall losses of vapors were not taken into account due to the lack of information on deposition onto EUPHORE walls.However, in the last section, the effect of vapor wall losses on SOA formation in the chamber was investigated by taking vapor wall losses into account according to parameters estimated by Zhang et al. (2014).

SOA mechanism
A new mechanism was developed for SOA formation from toluene (TOL), o-xylene (XYL) and trimethylbenzene (TMB).Parameters were fitted on data issued from several studies under low and high NO x conditions: Ng et al. (2007b) for toluene SOA, Cocker III et al. (2001) for TMB SOA and the results from both studies for xylene SOA.No information was found on SOA formation from TMB oxidation under low NO x conditions; therefore, the stoichiometric coefficient from the low NO x conditions for xylene was used.Similarly to Vivanco et al. (2016), the low NO x condition yields are only used if radicals formed from the oxidation of the precursors react at least twice with the HO 2 radical, to prevent high formation of low NO x SOA under intermediate NO x conditions.Molecular structures (used by SOAP to estimate several properties and to compute activity coefficients) were selected based on the results of Im et al. (2014) by selecting the compound with a similar volatility, formed in the largest quantity and ensuring the best reproduction of the O / C and H / C ratios.Properties of surrogate species and reactions are shown in Table 3.Reactions leading to SOA formation are shown in Table 4 for TOL SOA and in Table 5 for XYL and TMB SOA.
For monoterpenes, the mechanism of Pun et al. (2006) was updated using results from more recent studies.For SOA formation from α-pinene, β-pinene and limonene ozonolysis, the mechanism is based on the parameterizations described in Lee et al. ( 2011) for high and low NO x conditions.The API + OH reaction is based on the results of Svendby et al. (2008), whereas for the LIM + OH reaction, the yield of BiA2D was optimized to give the best results.Reactions and properties of surrogate species are shown in Tables 6 and 3.The terpene + NO 3 yields are based on Fry et al. (2014).
In the surrogate SOA approach, results of the model may depend on the choice of the surrogate structure to represent SOA formation.Although different VOCs form different compounds (or the same compounds in different quantities), this approach lumps together species that have similar properties (like saturation vapor pressure) and that are expected to have similar chemical structure.For example, firstgeneration aldehyde products from terpene oxidation were represented by the surrogate species BiA0D.Moreover, the compounds were selected (based on the best information available on SOA products) to reproduce the mean properties of SOA and of the O / C and H / C ratios.Kim et al. (2014) showed that the mean O / C and H / C ratios are around 0.34-0.36and 1.4-1.5 for SOA from α-pinene and limonene, respectively.The selected surrogates give a O / C ratio between 0.3 and 0.44 and a H / C ratio between 1.55 and 1.6.For TOL SOA, the chosen surrogates give a O / C ratio (between 0.71 and 0.8) and H / C ratio (between 1.14 and 1.6) similar to the reported values by Schilling (2015) (between 0.7 and 0.8 for O / C and between 1.2 and 1.6 for H / C).For TMB SOA, Sato et al. (2012) indicate a O / C ratio around 0.4-0.6 and a H / C ratio around 1.5-1.7.The chosen surrogates reproduce the O / C ratio (between 0.38 and 0.55) but seem to underestimate the high H / C ratio (between 1.13 and 1.44 due to the low H / C ratio of the AnIP2 compound).However, for AnIP2, a better molecular structure with a similar O / C ratio and saturation vapor pressure could not be found.
The mechanism of Couvidat and Seigneur (2011) was used to simulate SOA formation from isoprene.).b Saturation vapor pressure (torr) at 298 K. c Enthalpy of vaporization (kJ mol −1 ).

Oligomerization
Carlton et al. ( 2010) developed a parameterization for oligomerization based on a simple first-order rate constant of oligomerization.Based on the results of Kalberer et al. (2004) indicating 50 % of polymers after 20 h in TMB SOA, Carlton et al. (2010) determined a kinetic constant of 9.6 × 10 −6 s −1 .However, this number was obtained with laser desorption ionization-mass spectrometry (LDI-MS) by taking all compounds with m/z higher than 400 (with m the ion mass and z the ion charge) into account and not always taking small oligomers (dimers or trimers that may be formed more rapidly) into account.Kalberer et al. (2004) and Kalberer et al. (2006) also studied the oligomer fraction based on a volatility tandem differential mobility analyzer (VTDMA) giving the remaining volume fraction of particles at different temperatures.They found that for TMB SOA, the remaining fraction at 100 • C after 5 h ranges from 50 % to 62 % and is composed of small oligomers or very low-volatility  See TOLP + A reaction * Oxidants may be present as both reactants and products, so a reaction added to RACM2 will not affect the original photochemical oxidant concentrations.MO 2 and ACO 3 are the methyl peroxy radical and the peroxyacetyl radical, respectively.A is either NO, NO 3 , MO 2 or ACO 3 .organic compounds.The remaining fraction at 100 • C can even reach 80 % to 90 % after 25 h.Based on these results, and assuming that the remaining fraction is mainly composed of small oligomers and that oligomerization is irreversible, the first-order constant of oligomerization should be around 3.85 × 10 −5 s −1 to take small oligomer formation into account.
However, a first-order complete reaction may not be appropriate to represent oligomerization.Indeed, Roldin et al. (2014) showed that oligomerization should involve secondorder reversible reactions like esterification, hemiacetalization, aldolization and peroxyhemiacetalization.The equilibrium of these reactions is unfavored by humid conditions and the reaction is catalyzed under acidic conditions.Indeed, oligomer formation by esterification was reported in the case of isoprene SOA (Surratt et al., 2006).
In this study, oligomerization is represented by a reversible process which is mainly due to mechanisms like esterifica- tion, unfavored by humid conditions.It is represented by either a simple "reduced" reaction (or "bulk oligomerization") or an extended representation.For the extended representation, each of the oligomerization steps (each combination up to the formation of tetramer) is represented: where A is a monomer compound, A dimer is the dimer, A trimer is the trimer of compound A and A tretramer the tetramer of compound A, k oligo, extended is the forward kinetic of oligomerization and k reverse, extended is the reverse kinetic of oligomerization.In the case of a mixture of several compounds, each combination and each product would need to be represented.The same values of k oligo, extended and k reverse, extended are used for each steps.If more data were available, the parameterization could be improved to take different kinetic rate parameters of oligomerization per combination into account.The forward kinetic of oligomerization can be linked to the reverse kinetic of oligomerization via the oligomerization equilibrium constant: For the reduced (or bulk oligomerization) parameterization, all the steps are described by a single reaction: where A is a monomer compound, A oligo is the monomer blocks of compound A inside oligomers, m oligo is the number of monomer blocks inside oligomers, k oligo is the kinetic rate parameter of oligomerization and k reverse is the kinetic rate parameter of the reverse reaction.In this parameterization, oligomers are represented by a single model species.A represents the monomer and A oligo any monomer blocks present inside the oligomers.It can be use to represent oligomerization in a mixture of compounds (for example a mixture of compounds A and B).If a monomer block from a compound A reacts with another monomer block (the same compound A or another compound B), the compound will be converted from monomer A to oligomer A oligo (independently of the reaction with A or B).The parameterization currently does not take different kinetic rate parameters between each combination of compounds into account due to lack of data.It assumes that a compound A reacting with itself and a compound B reacting with itself will have the same kinetic parameter and will have the same value when compounds A and B react together.If more data were available, the param-eterization could be improved to take different kinetic rate parameters of oligomerization per combination into account.
In this study, the net flux of oligomerization J oligo is computed using activities.Activity is often seen as the "apparent concentration" of a compound in thermodynamics.It is linked to the chemical potential (molar Gibbs free energy of a particular component) by the following equation: where a i is the activity of compound i, µ i is the chemical potential of compound i, µ 0 i is the chemical potential under standard conditions, R is the gas constant and T is the thermodynamic temperature.
Activities (calculated here on the mole fraction basis) are used instead of concentrations for two main reasons.First, chemical rates are more consistent with thermodynamic equilibrium by computing rates using activities.For example, in the case of a simple one product (A) giving one product (B) equilibrium reaction, if chemical reactions are written using concentrations, the net flux of reaction J would be computed with the following equation: where k 1 is the forward kinetic parameter, k −1 is the reverse kinetic parameter, C A is the concentration of compound A and C B is the concentration of compound B. At equilibrium, J would be equal to zero and the equilibrium constant would then correspond to the ratio of concentrations instead of a ratio of activities.This paradox can be lifted by using activities instead of concentrations.Second, some studies (Madon and Iglesia, 2000;Rahimpour, 2004) expressed the need to compute chemical rates using activities and showed that better results are obtained for nonideal systems.
The net flux of oligomerization J oligo for the reduced parameterization (Eq. 3) is therefore computed with the following equations: where X a, monomer is the molar fraction of compound a, a a, monomer is the activity of compound a on a molar fraction basis and a a, oligomer is the activity of the oligomer formed from compound a on a molar fraction basis.Activities are computed with the AIOMFAC model (Zuend et al., 2008(Zuend et al., , 2011;;Zuend and Seinfeld, 2012;Ganbavale et al., 2015).The kinetic rate of oligomerization k oligo is computed as follows: where k max oligo is the maximum kinetic rate parameter for oligomerization and a monomer is the sum of monomer activities.
To calculate the reverse kinetic rate parameter, the computation is based on the equilibrium oligomerization constant K eq oligo .The equilibrium constant for oligomerization (due to esterification or a similar oligomerization mechanism) is computed with where a H 2 O is the activity of water on a molar basis.At equilibrium, the rate of oligomerization is zero.Therefore, To represent oligomerization, k max oligo , m oligo and K eq oligo were estimated based on the results of Kalberer et al. (2006).In their study, the molar masses of heavy molecules in SOA from several precursors were measured using matrix-assisted laser desorption-ionization mass spectrometry (MALDI-MS).To represent oligomerization, we represented the formation of oligomers up to tetramers from a single monomer explicitly using the extended parameterization (Eq. 1) to simulate the evolution of mass averaged molar mass of oligomers (which cannot be simulated with our reduced parameterization (Eq. 3) for oligomerization because the molar mass of oligomers does not vary with this parameterization).The equilibrium constant of oligomerization K eq oligo and the oligomerization constant k oligo, extended of the extended parameterization (Eqs. 1 and 2) were fitted so that the evolution as a function of time of the weight average molar mass of oligomers during the first hours of the experiment is reproduced by the model as shown by Fig. 1.A molar mass of the TMB SOA monomer of 155 g mol −1 was chosen so that the molar mass of TMB SOA is between 140 and 170 g mol −1 based on Im et al. (2014).m oligo and k max oligo are chosen so that the molar averaged mean molar mass and the total mass of oligomers (with oligomers and monomers) of the reduced (Eq.3) and the extended (Eq. 1) parameterizations were similar.The reduced parameterization values (Eq. 3) found for k max oligo , m oligo and K eq oligo were respectively 2.2 × 10 −4 s −1 , 3.35 and 2.94 for TMB SOA.k max oligo is in the same order of magnitude as the kinetic rate parameter reported for the reaction acetaldehyde and methanol by Roldin et al. (2014), 4.9 × 10 6 a[H + ] M h −1 (a[H + ] being the activity of H + ), which should be around 6 × 10 −4 s −1 on an activity basis for a pH of 4.6 (order of magnitude of pKa for carboxylic acids).Equation (3) could also be used to take the oligomerization inside an aqueous acidic phase into account if the effect of acidity is taken into account.Based on a pH of 4.6, a kinetic rate of 8.76 a[H + ] (instead of k max oligo in Eq. 3) should be used to take the effect of acidity on oligomerization into account.
This reduced parameterization (Eq. 3) only takes the formation of short oligomers (oligomers of two to four monomer blocks that can be formed quickly during the first hours) into account but should give a good insight into the impact of oligomerization on SOA formation.Formation of short oligomers can impact SOA formation by transforming monomers and increasing the condensation of semi-volatile compounds onto the particle.However, the formation of big oligomers (more than four blocks of monomers) should affect SOA formation less as it should mainly lead to the transformation of oligomers into bigger oligomers (oligomers with higher molar masses).Big oligomers could affect the partitioning of monomers indirectly by increasing the mean molar mass of the organic phase, the partitioning constant of monomers being inversely proportional to the mean molar mass.
This parameterization also gives good results for the isoprene SOA oligomerization using a molar mass of 120 g mol −1 .This molar mass corresponds to the molar mass of methyl glyceric acid, which was shown to undergo oligomerization by esterification.For α-pinene, Kalberer et al. (2006) did not find any significant temporal evolution of oligomer molar masses, which is consistent with the formation of dimers that cannot react further and therefore parameters for oligomerization of α-pinene SOA cannot be evaluated.
The results for SOA oligomerization from TMB and isoprene for the extended (Eq. 1) and reduced parameterizations (Eq. 3) are shown in Fig. 1.For the formation of dimers from α-pinene, we used the same parameters except that m oligo is set to two to limit the formation of oligomers to dimers.De-Palma et al. ( 2013) confirmed that particle-phase dimer formation is possible.However, Kristensen et al. (2014) studied the formation of four dimers and determined that those four dimers are not formed from particle-phase reaction but through gas-phase reactions of the stabilized Criegee inter-mediate formed from the ozonolysis of α-pinene.It could then be possible that in the case of α-pinene not all the oligomers are formed via particle-phase reactions.
In this study, the second-order reduced parameterization (Eq. 3) was used for simulations.

Uptake of pinonaldehyde onto acidic aerosols
Several studies (Liggio and Li, 2006a, b) reported an uptake of pinonaldehyde onto acidic aerosols higher than what could be predicted by assuming equilibrium between the gas and particle phases and no chemical reaction inside the particles.This phenomenon is attributed to oligomer and/or organosulfate formation.Gao et al. (2004b) reported a similar phenomenon for various aldehydes.To represent this phenomenon, Pun and Seigneur (2007) developed a parameterization by computing an effective Henry's law constant H eff for aldehyde compounds: where H is the monomer Henry's law constant of the compound and a(H + ) is the activity of protons in the aqueous phase.As fine particles are generally very acidic (Ludwig and Klemm, 1990;Keene et al., 2004), uptake of pinonaldehyde will in fact appear to be an irreversible process, even though the parameterization is formulated as a reversible process.Using this parameterization for pinonaldehyde in 3-D air quality models leads to very high concentrations of SOA from monoterpenes (Couvidat et al., 2012).However, this parameterization does not take the uptake rate of aldehydes into account.Liggio and Li (2006b) measured the uptake rate coefficients for various acidities of the aerosols.The authors showed that the uptake of pinonaldehyde onto aerosols can only be significant for very high acidities (which can be reached with low ammonia concentration and at low humidities).Based on these results, a parameterization for the chemical evolution of pinonaldehyde has been developed in this study: where BiA0D is the surrogate species used by Couvidat et al. (2012) to represent pinonaldehyde, BiA0D part is the amount of pinonaldehyde inside the particle and k trans (in s −1 ) is the kinetic rate of BiA0D transformation inside the particle into a product assumed non-volatile.Liggio and Li (2006b) did not evaluate a kinetic rate of pinonaldehyde transformation inside the particle but measured a kinetic rate of uptake.However, the kinetic parameter of uptake k uptake can be linked to the kinetic parameter of transformation k trans by assuming equilibrium between the gas and particle phases: where K aq is the partition coefficient of pinonaldehyde between the gas phase and the particle and AQ is the particle mass.
As in Liggio and Li (2006b), the pH of particles and activities of compounds was computed with AIOMFAC (Zuend et al., 2008(Zuend et al., , 2011;;Zuend and Seinfeld, 2012;Ganbavale et al., 2015) depending on the conditions (humidity, temperature, concentrations, etc.) of the experiments of Liggio and Li (2006b).The estimated flux of transformation based on the results of Liggio and Li (2006b) was plotted against several variables to find the variables exhibiting the best correlations.m H + , the molality of ion H + , and a Two possible parameterizations can therefore be used to take the possible transformation of pinonaldehyde in the particle phase into account: the pH-dependent parameterization and the HSO − 4 -dependent parameterization.For the pHdependent parameterization, the kinetic rate is computed as whereas the HSO − 4 -dependent parameterization is computed with In that case, we assumed that the reaction leads to the formation of the organosulfate formed from pinonaldehyde.The reaction is assumed to be complete as according to the parameterization of Pun and Seigneur (2007); the uptake of pinonaldehyde onto an acidic aerosol should be an irreversible process.
The pH-dependent parameterization (Eq.13) could be representative of an oligomerization mechanism catalyzed by a H + ion.However, in that case, the uptake of pinonaldehyde should be seen as a reversible pathway, which also depends on the humidity of the experiments.The HSO − 4dependent parameterization (Eq.14) could be representative of organosulfate formation, and therefore pinonaldehyde should in that case act as a sink for sulfates due to the reaction between sulfates and pinonaldehyde.

Aging mechanism
To test the influence of aging on SOA formation, an aging mechanism was developed.The aging of BiA0D is based on the results of Chacon-Madrid and Donahue (2013) who studied SOA formation from pinonaldehyde oxidation.SOA formation from BiA1D is based on the SOA yield from the oxidation of pinonic acid as measured by Müller et al. (2012).2015) * Oxidants may be present as both reactants and products so that a reaction added to RACM2 will not affect the original photochemical oxidant concentrations.MO 2 and ACO 3 are the methyl peroxy radical and the peroxyacetyl radical, respectively.
For the aging of BiA2D, the parameterizations of Jathar et al. (2015) were used to determine the amount of functionalization and the decrease of volatility due to aging.Nopinone, which is formed from the oxidation of β-pinene, was shown to form a significant amount of SOA and low-volatility products (Sato et al., 2016) and was included in the mechanism.SOA yields from nopinone oxidation were based on Mutzel et al. (2016).The yields of formation of nopinone were based on Hakola et al. (1994) for the reaction of β-pinene with OH and O 3 and on Hallquist et al. (1999) for the reaction of β-pinene with NO 3 .Kinetics for aging were taken from the Master Chemical Mechanism v3.3.1 (Jenkin et al., 1997;Saunders et al., 2003).
For the aging of aromatic SOA compounds, the parameterizations of Jathar et al. (2015) were used to determine the amount of functionalization and the decrease of volatility due to aging.According to this study, oxidation of aromatics should mainly lead to the addition of one oxygen atom.For simplification purposes, we assumed that aging of aromatic SOA leads to the addition of a single hydroxy group.
Table 7 shows the aging mechanism used in this study.This simple mechanism is used here in order to be able to assess the capacity to form oligomers and the reactive uptake of pinonaldehyde in the long-term formation of SOA.

Long-term SOA formation simulations
The experimental conditions (low relative humidity, oxidation over a few hours, high organic mass loading) are probably very different from typical ambient conditions.To eval- uate the potential effect of oligomerization more precisely, SOA formation from the oxidation of α-pinene, toluene and TMB over 3 days of evolution with an organic mass loading of 5 µg m −3 was simulated for various humidities with or without oligomerization.The simulated effect of oligomerization is compared to the simulated effect of aging on SOA formation.Octane and NO x concentrations are set in order to reproduce a level of OH concentrations similar to summer conditions (around 0.001 ppb during daytime).The diurnal cycle was simulated by computing the evolution of the zenith angle at Valencia, Spain, as a function of local time.This diurnal profile was simulated to take the fact into account that oxidation and aging slow down during night (due to the low concentrations of radicals), whereas oligomerization continues.On the other hand, the relative humidity was assumed constant (whereas under ambient conditions relative humidity probably has a diurnal profile).The structures of pinonic acid for simulations with α-pinene and of AnRP2 for simulations with toluene and TMB oxidation are used for the structure of the organic loading.Those simulations F. Couvidat et al.: Simulating SOA from anthropogenic and biogenic precursors are not representative of atmospheric conditions but can be used to illustrate the effect of long-term oligomerization on SOA formation.For simulations with oligomerization, the fraction of monomers inside the aerosol is assumed to be 1; i.e., the aerosol is assumed to be mainly constituted of monomers that can react with absorbed compounds.These simulations therefore provide information on the maximal effect that oligomerization can have on SOA yields.

Results
The simulated degradation of precursors is illustrated in Figs.S1 and S2 in the Supplement.Whereas the mechanism was able to reproduce the degradation of VOCs during the beginning of the experiments, medium-term oxidation of toluene seems to be underestimated for some experiments (A1, A3, A4 and A7).Due to this, medium-term formation of toluene SOA may be underestimated.

Comparison with measurements
To ensure the satisfactory performances of the developed SOA mechanism under uncontrolled conditions, the results of the SOA mechanism are compared to the measurements carried out inside EUPHORE.Figures 3 and 4 show the results of the model for the biogenic experiments without SO 2 and the anthropogenic experiments, respectively.Simulated SOA composition for simulations considering oligomerization is illustrated in Figs.S3 (biogenic experiments) and S4 (anthropogenic experiments).
Usually, in environmental chamber studies, the wall deposition rate is based on the evolution of concentrations during the last hours of experiments.However, in the last stage of the experiments, oligomerization could still occur and affect the evolution of concentrations and therefore affect the estimation of the wall deposition rate.The estimated deposition rate could be biased if the evolution of concentrations due to oligomerization at the end of the experiment is significant compared to the deposition rate.In these simulations, the wall deposition rate was constrained to reproduce with the model the decrease of SOA volume concentrations (measured with the SMPS) during the last hours of the experiments.Because the computed evolution of SOA concentrations during the last hours can be slightly different with or without oligomerization, the wall deposition rates used with and without oligomerization are different.To examine the effect of the wall deposition rate, the formation of SOA with oligomerization was also simulated with the deposition rate as computed by the simulation without oligomerization to ensure that the effects of changes on deposition remain low (changes of concentrations of a few percent) and that the hypothesis on oligomerization does not affect the estimated deposition rate.
For the biogenic experiments, the model gives good results (bias lower than 20 %) with or without oligomerization for all experiments, with slightly better results without oligomerization for experiment B5 and slightly better results with oligomerization for experiment B1.The experiments B2, B3, B4 and B5 have previously been modeled by Santiago et al. ( 2012a) using a box model version of two air quality models: CMAQ (Carlton et al., 2010) and Chimere (Menut et al., 2013).The authors found a significant overestimation of modeled SOA which was not found with the model used in this study.
For the anthropogenic experiments, the model gives satisfactory results (bias lower than 20 %) with or without oligomerization for experiment A1, A3, A6 and A7 but overestimates concentrations by 30-40 % for A2 and A5 and underestimates concentrations by 25 % for A8 and A9.
Concentrations for experiment A4 are overestimated by a factor of 2. However, experimental conditions for experiment A4 are close to those of experiment A1, except that HONO concentrations are 2 times higher than for experiment A4.The model gives similar results for A1 (for which the model gives satisfactory results) and A4.It could not explain why concentrations were so different between these two experiments.One possibility is that anthropogenic SOA formation is inhibited at high HONO concentrations.
The alternation of overestimation and underestimation events could be related to differences in chemical regimes.To compare the chemical regimes of the experiments, a ratio representative of the chemical regime (called hereafter "chemical regime ratio") was computed with the sum of VOC concentrations multiplied by their reactivity with OH divided by NO x concentrations: where C r is the chemical regime ratio, C i is the concentration of each VOC in ppb, C NO x is the concentration of NO x in ppb and k OH, i is the reactivity with OH in molecule −1 cm 3 s −1 ; 10 10 is a factor set so that the values are close to unity.The chemical regime ratio was used instead of the VOC / NO x ratio because in this study a mixture of VOCs (and not a single VOC) was present in the chamber.The chemical regime ratio takes the reactivity of the compounds into account and can therefore be used to compare different experiments with different mixtures of VOCs.
It appears that depending on the chemical regime ratio, the model may overestimate or underestimate the concentrations.Indeed, for the experiments with a low chemical regime ratio (experiments A8 and A9, which have ratios of 0.13 and 0.7, respectively) SOA concentrations are underestimated.On the contrary, SOA concentration experiments with a high chemical regime ratio (experiments A2 and A5, with ratios equal to 2.2 and 1.64, respectively) are overestimated.Except for A4, the model gives satisfactory results for the other experiments (having chemical regime ratios between 0.77 and 1.05).It may indicate a more complex NO x dependency than what is represented in the model.

Effect of oligomerization on SOA formation
For the biogenic and anthropogenic experiments, a low effect of oligomerization on SOA formation was simulated by the reduced parameterization (Eq. 3) as shown by Figs. 3  and 4. Concentrations are indeed similar with and without oligomerization with low differences compared to the amount of SOA.Although the amount of SOA does not significantly change, the composition of SOA is strongly affected by oligomerization.Indeed, the oligomer content (in mass) at the end of the experiment varies from 68 % to 78 % for the biogenic experiments and from 38 % to 58 % for the anthropogenic experiments with a similar range of humidity (between 0.4 % and 37 %).For a 55 % relative humidity, Gao et al. (2004a) determined for α-pinene that the oligomer content of SOA could be well over 50 %.This result is consistent with the results of our parameterization as an oligomer content around 60 % was simulated for such an humidity.On the contrary, for the biogenic experiments, the first-order complete oligomerization reaction of Carlton et al. (2010) with a kinetic rate parameter of 9.6 × 10 −6 s −1 (halftime of 20 h) simulates a low number of oligomers in SOA (below 15 %).Even with the kinetic rate parameter of 3.85 × 10 −5 s −1 (kinetic estimated from the remaining fraction at 100 • C measured by Kalberer et al., 2006), the number of oligomers estimated (between 29 % and 40 %) is below the number estimated by the reduced parameterization (Eq.3).These results indicate that even though the second-order reduced parameterization (Eq. 3) is not complete, it leads to a faster formation of oligomers and therefore, the parameterization of Carlton et al. (2010) may underestimate the short-term formation of oligomers.On the other hand, as the parameterization of F. Couvidat et al.: Simulating SOA from anthropogenic and biogenic precursors Carlton et al. (2010) is complete, it may also lead to an overestimation of the long-term formation of oligomers.
The results of long-term SOA formation simulations are shown in Figs. 5,6 and 7. Without taking aging into account, oligomerization leads to a significant increase of concentrations after 3 days of evolution, even at high humidity.Increase factors due to oligomerization are 2.5 for α-pinene, 2.0 for toluene and 6.2 for TMB for a relative humidity of 30 % and are 1.8 for α-pinene, 1.4 for toluene and 3.3 for TMB for a relative humidity of 70 %.Oligomerization is therefore a process that could significantly affect long-term SOA formation.
Assuming aging leads to a slight decrease of SOA mass for toluene SOA due to fragmentation, while it leads to an increase of concentrations for TMB and α-pinene SOA due to functionalization.However, the effect of aging on SOA formation simulated here seems less important than the effect of oligomerization.
The different parameterizations of oligomerization are also compared in Fig. 5: the equilibrium second-order reaction parameterization developed in this study, the first-order complete reaction of Carlton et al. (2010) with a kinetic rate parameter of 9.6 × 10 −6 s −1 (halftime of 20 h) and a kinetic rate parameter of 3.85×10 −5 s −1 (halftime of 5 h).The three different parameterizations have different impacts on SOA formation.The first-order parameterization with a kinetic rate parameter of 9.6 × 10 −6 s −1 has a low impact on SOA yields except after a few days of oligomerization.With the kinetic rate parameter of 3.85 × 10 −5 s −1 , the increase of SOA yields can be significant after 10 h.For RH = 70 % the SOA yield increases by 72 % after 1 day.The second-order parameterization gives a faster SOA production.Therefore, oligomerization may be a faster process than simulated in previous modeling studies.However, this figure shows the maximal effect of the second-order parameterization (molar fraction of oligomers of 1).Depending on conditions, the second-order parameterization may only have a low impact on SOA concentrations and the formation of oligomers could rapidly reach an equilibrium.
Depending on the conditions, oligomerization could also lead to a decrease of SOA concentrations as the partitioning is inversely proportional to the mean molar mass of organic aerosols (Pankow, 1994), which will increase with the oligomer formation.Increasing the mean molar mass by a factor of 2 leads to a decrease around 40 % of the SOA formed without oligomerization, indicating that the partitioning of monomers is sensitive to the value of the mean molar mass.

Uptake of pinonaldehyde onto acidic aerosols
The pH (Eq.13) and HSO − 4 (Eq.14) parameterizations and the parameterization of Pun and Seigneur (2007) (Eq.10) assuming equilibrium are compared and evaluated against experiments B6 and B7, in which α-pinene and limonene are oxidized in the presence of SO 2 .The oxidation of SO 2 leads to the formation of sulfuric acid and therefore to acidic aerosols.The amount of SOA formed with each parameterization is compared to the results of the experiments (Fig. 8).
The parameterization of Pun and Seigneur (2007) (Eq.10) leads to a significant overestimation of SOA concentrations.With this parameterization, BiA0D is entirely absorbed by the acidic aerosol.On the contrary, using the pH (Eq.13) and HSO − 4 (Eq.14) parameterizations which take the dynamic of the uptake into account, no significant formation of SOA is formed by this pathway.The two parameterizations give almost the same results as the assumption of no uptake.These results indicate that the dynamics of the uptake of pinonaldehyde need to be taken into account.
Pinonaldehyde was found to be too volatile to form significant SOA by this pathway.The long-term formation of SOA by this pathway is tested with the simulation described in the previous section but in the presence of 2 µg m −3 of sulfuric acid.Fig. 5 shows that even at low humidity (RH = 30 %), the amount of SOA formed is not significant after 3 days of evolution.Therefore, SOA formation from the reactive uptake by acidic aerosol of pinonaldehyde probably does not contribute significantly to SOA formation.
However, it could be possible that aldehyde compounds less volatile than pinonaldehyde react and form a significant amount of SOA by this pathway.To test this hypothesis, we assumed that organosulfate could be formed from the compounds AnRP1 and AnRP2 (which have an aldehyde group in their molecular surrogate structure) using the HSO 4 parameterization (Eq.14).The kinetic rate parameter should depend on the compound, but the HSO 4 parameterization (Eq.13) probably provides a good estimation of this phenomenon.Figure 9 shows the amount of organosulfate that would be formed by this pathway for a humidity of 70 % from the oxidation of toluene.Long-term organosulfate formation seems possible by this pathway (even at high humidity) as a significant amount of organosulfate (13 % to 18 % of SOA) is formed with this assumption and a significant mass of aldehydes has been converted into organosulfates.However, if high concentrations of ammonia are present in the atmosphere, pH will increase and HSO − 4 will decrease.In the case of high concentrations of ammonia, the kinetic rate could be too low for this process to be significant.This process would need for sulfate to not be fully neutralized by ammonia.

Investigation of the effect of particle-phase diffusion and wall losses of organic vapors
To evaluate the effect of the particle-phase viscosity on SOA formation, the SOAP model was run for a particle-phase diffusion coefficient of 2 × 10 −16 molecules cm −2 s −1 , which corresponds to the order of magnitude of values determined for toluene SOA at low humidities (Song et al., 2016).To investigate wall losses of organic vapors, wall losses were treated as in Zhang et al. (2014) by using a first-order wall loss rate and the organic aerosol equivalent wall (set to 10 mg m −3 following Zhang et al., 2014).According to Mc-Murry and Grosjean (1985), the wall vapor loss rate k w can be computed with the following formula: where A / V is the surface on volume ratio of the chamber (equal to 1 m −1 for EUPHORE), a w is the mass accommodation coefficient of organic vapors onto the wall, c is the mean thermal speed of the molecules, k e is the coefficient of eddy diffusion and D gas is the molecular diffusivity.
As the mixing is carried out by two fans with a flux of 4000 m 3 h −1 , the characteristic time for the mixing inside the chamber is estimated to be around 90 s for EUPHORE and k w is estimated to be around 10 −4 s −1 .Following Zhang et al. (2014), the kinetic of evaporation of compounds from the wall k w, off is computed with where K p, w is the gas-wall partition coefficient (chosen equal to the gas-particle partition coefficient) and C w is the equivalent wall organic aerosol (OA) concentration (set to 10 000 µg m −3 ).The effects of the particle-phase viscosity and of wall losses of vapor are illustrated in Fig. 10 for experiments A5, A7 and B1.Results for other experiments (biogenic and anthropogenic) are similar to these experiments and are illustrated in Figs.S5 and S6.
Generally, similar results are obtained between nonviscous and viscous aerosols.Assuming viscous aerosols can lead to an increase of SOA concentrations due to the limitation of evaporation with low diffusion.As found by Couvidat and Sartelet (2015), condensation of low-volatility com-   pounds is possible without diffusing into the particle by condensing at the gas-particle interface.The condensation of low-volatility compounds can create a layer onto which more volatile compounds can condense to respect Raoult's law at the interface.However, viscosity can prevent the evaporation of those more volatile compounds by preventing their diffusion from the core of the particle to the interface.For non-viscous aerosols, the deposition of particles to the walls leads to a decrease of the absorbing mass (mass of the organic aerosol).As the gas-particle partitioning depends on the absorbing mass, SVOCs will evaporate to maintain the gasparticle partitioning, whereas this evaporation will be limited for a viscous aerosol.
It can be assumed that viscous aerosols can also lead to a decrease of SOA sensitivity to changes of conditions during the experiments.For example, experiment A5 is characterized by a decrease of the temperature during the experiment leading to a decrease of volatility.Whereas concentrations increase when assuming non-viscous aerosols, concentrations seem unaffected by the change of temperature when assuming viscous aerosols (as compounds are absorbed by the core of particles with a slow kinetic rate).The shape of the modeled SOA concentration curve is closer to measurements when assuming viscous aerosols.
Taking wall losses of vapors into account leads to a decrease of SOA concentrations (as illustrated in Fig. 10), particularly for the anthropogenic experiments (decrease between 8 % and 30 %) due to the lower organic aerosol loading compared to the biogenic experiments.As the developed mechanism is based on the methodology of Odum et al. (1996) by using experimental results from Teflon chambers, it could be possible that the parameters of the mechanism (saturation vapor pressures and stoichiometric coefficients for the formation of SVOCs) are biased.As stronger organic aerosol loading favors condensation onto organic aerosol over condensation onto walls, wall losses of vapor may produce a shift in apparent volatilities; the organic aerosol formed in chambers may appear more volatile than it should be.For both the biogenic and anthropogenic experiments, results can be corrected by decreasing the saturation vapor pressures by 20 % as better or similar results could be obtained than without taking wall losses of organic vapors into account.These results indicate that volatility determined by the use of Odum's curves may be slightly overestimated.However, wall deposition of organic vapors could also lead to an underestimation in the SOA mechanism of stoichiometric coefficients.

Conclusions
Several parameterizations were developed in this study.First, in the spirit of updating the H 2 O mechanism (Couvidat et al., 2012), a new mechanism was developed to take hydrophilic and hydrophobic properties of monoterpene and aromatic SOA into account and introduce aging in the mechanism.The performance of the new mechanism was evaluated by comparison to experimental results from previous studies carried out in the EUPHORE chamber in Valencia, Spain.Second, parameterizations to take oligomerization and the uptake of pinonaldehyde onto acidic aerosols into account were devel-oped.Finally, the effects of particle viscosity and wall deposition of vapors were investigated.
Without taking wall losses of vapors into account, satisfactory results are obtained for the biogenic experiments.However, the experiments used for the comparison present a high organic aerosol loading (between 50 and 150 µg m −3 ), much higher than the values normally observed under atmospheric conditions.More experiments at low organic aerosol loading should be carried out to provide a wider observational dataset that could be used to evaluate model developments.A good performance is also found for most of the anthropogenic experiments.However, for the experiments with the lower chemical regime ratio (Eq.15; product of VOC concentrations and their reactivity divided by NO x concentrations), the model underestimates concentrations by 30 %, whereas for the experiments with the higher chemical regime ratio, the model overestimates concentrations by 30 %.This could indicate a more complex NO x chemistry than the one taken into account in the model.However, this discrepancy could be due to the difficulty of simulating radicals inside the chamber under some conditions properly.More experiments should be carried out to confirm these results and improve the NO x dependency inside the model.
Contrary to Carlton et al. (2010), who represented oligomerization as a first-order irreversible process, oligomerization was represented in the present study (as an effort to provide a more realistic representation) as a second-order reversible reaction unfavored by humid conditions by analogy to oligomerization reactions expected to occur in the particles, such as esterification, hemiacetalization, aldolization and peroxyhemiacetalization.With this new parameterization, oligomerization was shown to have little impact on SOA mass during the experiments.However, oligomerization was found to significantly affect SOA composition as a large part of SOA is constituted (more than 50 % for the biogenic experiments) by oligomers according to the simulations.These results seem to be consistent with the high number of oligomers reported by Gao et al. (2004a).Conversely, a low number of oligomers was simulated with the parameterization of Carlton et al. (2010) (below 15 % at the end of the experiments).The results of the present study with the new parameterization indicate that oligomerization may be faster than simulated by the parameterization of Carlton et al. (2010) but at the same time may reach an equilibrium.However, more effort should be given to improve this parameterization.Indeed, this parameterization represents a bulk oligomerization and does not account for differences in the reactivity of monomers.This parameterization should nonetheless be employed and tested in 3-D air quality models to diagnose the effect of oligomerization on SOA formation.
Instead of representing the uptake of pinonaldehyde by an equilibrium phenomenon (like the parameterization of Pun and Seigneur, 2007), we chose to take the kinetic of this uptake into account, which was shown to limit the uptake of pinonaldehyde strongly (Liggio and Li, 2006b).The uptake of pinonaldehyde was found to be too slow to contribute significantly to the formation of organic aerosols, suggesting that the approach of Pun and Seigneur (2007) (Eq.10) greatly overestimates the effect of the particle acidity on SOA formation.However, it could be possible that the uptake of aldehydes less volatile than pinonaldehyde (or more reactive) is strongly influenced by the particle acidity, leading to the formation of oligomers or organosulfate.
Taking into account wall losses of organic vapors with the parameters of Zhang et al. (2014) led to a decrease of SOA concentrations of up to 30 %.As the mechanism developed in the present study is based on the methodology of Odum et al. (1996) by using experimental results from Teflon chambers, the parameters of the mechanism (saturation vapor pressures and stoichiometric coefficients for the formation of SVOC) could be biased.The simulated decrease of SOA concentrations could however be compensated by decreasing the volatility of SVOC by 20 % (decrease of the value of saturation vapor pressures by 20 %).
Regarding RH, the model indicates that RH may strongly influence SOA yields, which is consistent with the results of Healy et al. (2009).Unfortunately, the conditions of the experiments used in this study do not cover the full range of RH.More experiments carried out at high humidity could be used to evaluate the model performance for other humidities.
Finally, these parameterizations should be implemented in a 3-D air quality models to evaluate their impact on SOA formation in the atmosphere.Oligomerization should especially be studied in detail as it could lead either to an increase or a decrease (by increasing the mean molar mass that would lead to a decrease of the partitioning constant of monomers) of SOA concentrations.
Data availability.The data are accessible by contacting the corresponding author (florian.couvidat@ineris.fr).
The Supplement related to this article is available online at: https://doi.org/10.5194/acp-18-15743-2018supplementAuthor contributions.FC wrote the article, developed the model and simulations and analysed results.MV co-developed the model, provided the observational dataset and analysed results.BB was the initial investigator of the study and analysed results.
Competing interests.The authors declare that they have no conflict of interest.

Figure 1 .
Figure 1.Temporal evolution of molar masses (a) and of the mass fraction of oligomer (b) for the extended (Eq. 1) and reduced (Eq.3) parameterizations.
HSO − 4 on a molality basis, were found to be the best variables with high correlation coefficients (0.98 for m H + and 0.92 for a the kinetic of transformation pinonaldehyde k trans inside the particle.Figure2shows k trans (in s −1 ) as a function of m H + (in mol kg −1 ) and a

Figure 2 .
Figure 2. Kinetic rate of transformation of pinonaldehyde as a function of m H + , the molality of ion H + (a), and a (m) HSO − 4

Figure 3 .
Figure 3. Aerosol concentration formation for the biogenic experiments without SO 2 .Black lines correspond to SMPS measurements.

Figure 4 .
Figure 4. Aerosol concentration formation for the anthropogenic experiments.Black lines correspond to SMPS measurements.

Figure 5 .
Figure 5. Evolution of the SOA yield from α-pinene oxidation as a function of time for an organic mass loading of 5 µg m −3 .Solid lines correspond to SOA formation with aging.Dashed-dotted lines correspond to SOA formation without aging.

Figure 6 .
Figure 6.Evolution of the SOA yield from toluene oxidation as a function of time for an organic mass loading of 5 µg m −3 .Solid lines correspond to SOA formation with aging.Dashed-dotted lines correspond to SOA formation without aging.

Figure 7 .
Figure 7. Evolution of the SOA yield from trimethylbenzene oxidation as a function of time for an organic mass loading of 5 µg m −3 .Solid lines correspond to SOA formation with aging.Dashed-dotted lines correspond to SOA formation without aging.

Figure 8 .
Figure 8. Aerosol concentration formation for the biogenic experiments with SO 2 .Black lines correspond to SMPS measurements.

Figure 9 .
Figure 9. Yield of formation of SOA, aldehydes and organosulfates with (solid lines) or without (dashed lines) the conversion of aldehydes into organosulfates as a function of time for an organic mass loading of 5 µg m −3 with 2 µg m −3 of sulfates.

Figure 10 .
Figure10.Effect of viscosity and gas wall losses on SOA concentrations for several experiments.The black line corresponds to SMPS measurements, the red lines correspond to modeled SOA concentrations for the non-viscous aerosol assumption, the blue lines correspond to modeled SOA concentrations for the viscous aerosol assumption and the green line corresponds to modeled SOA concentrations for the viscous aerosol assumption with a decrease by 20 % of volatilities.Solid lines correspond to simulations assuming no gas wall losses, and dotted lines correspond to those assuming gas wall losses.

Table 1 .
Initial concentrations in ppb, temperature (T ) in Kelvin and relative humidity in % for biogenic experiments.

Table 3 .
Properties of the different surrogate SOA species.
a Molecular weight (g mol −1