Modeling organic aerosols in a megacity : potential contribution of semi-volatile and intermediate volatility primary organic compounds to secondary organic aerosol formation

It has been established that observed local and regional levels of secondary organic aerosols (SOA) in polluted areas cannot be explained by the oxidation and partitioning of anthropogenic and biogenic VOC precursors, at least using current mechanisms and parameterizations. In this study, the 3-D regional air quality model CHIMERE is applied to estimate the potential contribution to SOA formation of recently identified semi-volatile and intermediate volatility organic precursors (S/IVOC) in and around Mexico City for the MILAGRO field experiment during March 2006. The model has been updated to include explicitly the volatility distribution of primary organic aerosols (POA), their gas-particle partitioning and the gas-phase oxidation of the vapors. Two recently proposed parameterizations, those of Robinson et al. (2007) (“ROB”) and Grieshop et al. (2009) (“GRI”) are compared and evaluated against surface and aircraft measurements. The 3-D model results are assessed by comparing with the concentrations of OA components from Positive Matrix Factorization of Aerosol Mass Spectrometer (AMS) data, and for the first time also with oxygen-to-carbon ratios derived from high-resolution AMS measurements. The results show a substantial enhancement in predicted SOA concentrations (2–4 times) with respect to the previously published base case without S/IVOCs (Hodzic et al., 2009), both within and downwind of the city leading to much reduced Correspondence to: A. Hodzic (alma@ucar.edu) discrepancies with the total OA measurements. Model improvements in OA predictions are associated with the bettercaptured SOA magnitude and diurnal variability. The predicted production from anthropogenic and biomass burning S/IVOC represents 40–60% of the total measured SOA at the surface during the day and is somewhat larger than that from commonly measured aromatic VOCs, especially at the T1 site at the edge of the city. The SOA production from the continued multi-generation S/IVOC oxidation products continues actively downwind. Similar to aircraft observations, the predicted OA/1CO ratio for the ROB case increases from 20–30 μg sm−3 ppm−1 up to 60–70 μg sm−3 ppm−1 between a fresh and 1-day aged air mass, while the GRI case produces a 30% higher OA growth than observed. The predicted average O/C ratio of total OA for the ROB case is 0.16 at T0, substantially below observed value of 0.5. A much better agreement for O/C ratios and temporal variability (R 2=0.63) is achieved with the updated GRI treatment. Both treatments show a deficiency in regard to POA ageing with a tendency to over-evaporate POA upon dilution of the urban plume suggesting that atmospheric HOA may be less volatile than assumed in these parameterizations. This study highlights the important potential role of S/IVOC chemistry in the SOA budget in this region, and highlights the need for further improvements in available parameterizations. The agreement observed in this study is not sufficient evidence to conclude that S/IVOC are the major missing SOA source in megacity environments. The model is still very underconstrained, and Published by Copernicus Publications on behalf of the European Geosciences Union. 5492 A. Hodzic et al.: Modeling organic aerosols in a megacity other possible pathways such as formation from very volatile species like glyoxal may explain some of the mass and especially increase the O/C ratio.

discrepancies with the total OA measurements.Model improvements in OA predictions are associated with the bettercaptured SOA magnitude and diurnal variability.The predicted production from anthropogenic and biomass burning S/IVOC represents 40-60% of the total measured SOA at the surface during the day and is somewhat larger than that from commonly measured aromatic VOCs, especially at the T1 site at the edge of the city.The SOA production from the continued multi-generation S/IVOC oxidation products continues actively downwind.Similar to aircraft observations, the predicted OA/ CO ratio for the ROB case increases from 20-30 µg sm −3 ppm −1 up to 60-70 µg sm −3 ppm −1 between a fresh and 1-day aged air mass, while the GRI case produces a 30% higher OA growth than observed.The predicted average O/C ratio of total OA for the ROB case is 0.16 at T0, substantially below observed value of 0.5.A much better agreement for O/C ratios and temporal variability (R 2 =0.63) is achieved with the updated GRI treatment.Both treatments show a deficiency in regard to POA ageing with a tendency to over-evaporate POA upon dilution of the urban plume suggesting that atmospheric HOA may be less volatile than assumed in these parameterizations.This study highlights the important potential role of S/IVOC chemistry in the SOA budget in this region, and highlights the need for further improvements in available parameterizations.The agreement observed in this study is not sufficient evidence to conclude that S/IVOC are the major missing SOA source in megacity environments.The model is still very underconstrained, and 1 Introduction Organic aerosols (OA) are currently one of the most active areas of research in atmospheric chemistry and air quality.Composed of both directly emitted primary organic particles (POA) and photochemically produced secondary organic matter (SOA), they have been identified as one of the two largest constituents of submicron aerosols throughout the world, with similar concentrations as sulfate (e.g.Murphy et al., 2006;Zhang et al., 2007).Despite their abundance in the troposphere and recognized effects on climate and human health, OA emissions, ageing, and properties remain poorly understood, and as a result, poorly represented in chemical-transport models (Kanakidou et al., 2005;Hallquist et al., 2009).Biogenic SOA levels formed in clean regions appear to be well represented overall by "traditional" SOA models (e.g.Tunved et al., 2006;Hodzic et al., 2009;Chen et al., 2009;Slowik et al., 2010).In contrast, many studies have reported that observed levels of secondary organic species in polluted regions at both local and regional scales cannot be explained by current modeling assumptions and from the oxidation of known VOC gaseous precursors ("V-SOA", following the terminology introduced by Tsimpidi et al., 2010) (e.g.Heald et al., 2005;de Gouw et al., 2005;Volkamer et al., 2006;Pun and Seigneur, 2007;Simpson et al., 2007;Hodzic et al., 2009;Tsimpidi et al., 2010).Model discrepancies can be substantial for short oxidation times but have been found to increase downwind of source regions (e.g.Heald et al., 2005;Volkamer et al., 2006), which consequently affect the ability of models to predict aerosol regional concentrations and radiative forcing.The most commonly used traditional semi-empirical V-SOA formation approach (Odum et al., 1996) considers only the early generation oxidation of gaseous VOC precursors into two nonreactive semi-volatile products that are partitioned between gas and aerosol phases depending on temperature and OA mass concentration.This approach does not reflect the full spectrum of potential SOA precursors, or the complexity of the OA and VOC chemical aging (Robinson et al., 2007;Jimenez et al., 2010).
In the past few years, evidence for new precursors and pathways of SOA formation has been identified.In particular, the semi-volatile character of primary organic aerosols (POA) suggested by Robinson et al. (2007) has been the focus of much recent attention.Traditionally, POA that are directly emitted into the atmosphere from various combustion sources (e.g.vehicle exhaust, biomass burning, etc.) have been modeled as hydrophobic non-volatile compounds that remain in the particulate phase throughout their lifetime.Recent studies (Lipsky and Robinson, 2006;Robinson et al., 2007;Huffman et al., 2008Huffman et al., , 2009a, b;, b;Sage et al., 2008;Grieshop et al., 2009) have shown that POA species should be treated as thermodynamically active compounds that can evaporate from the particulate phase, react in the gas-phase, and repartition as SOA.In the conceptual model of Robinson et al. (2007), POA is co-emitted with two types of gases: (i) semi-volatile compounds (SVOC) which are characterized by a relatively low volatility and are in the substantial partitioning with POA, and (ii) intermediate volatility vapors (IVOC) which are highly volatile and remain preferentially in the gas-phase.S/IVOC species have traditionally been ignored in SOA modeling and are also typically neglected in VOC emission inventories.They are composed of species such as long-chain alkanes, and can undergo multiple oxidation steps leading to a large amount of lower volatility precursors that are likely to partition to the aerosol phase downwind of source regions.
Following Tsimpidi et al. (2010), we will refer to SOA formed from SVOC as "S-SOA, to SOA formed from IVOC as "I-SOA," and to the combination of S-SOA plus I-SOA as "SI-SOA."Although IVOC and SVOC species have been known to be present in urban atmospheres for at least a decade (Fraser et al., 1997), their role as potentially important SOA precursors was only identified in the recent work of Robinson et al. (2007) and related publications.SOA formation from some IVOCs has been established experimentally, e.g.naphthalene and heptadecane are two IVOCs whose SOA formation has been investigated in some detail in the Caltech and Carnegie Mellon chambers (Chan et al., 2009;Presto et al., 2009).Unfortunately many IVOCs and SVOCs are not separable by GC-MS and appear in the "unresolved complex mixture" (UCM), which led to the need to parameterize their amounts and SOA formation (Robinson et al., 2007).
One of the key uncertainties, in addition to the chemical properties of S/IVOC vapors, concerns the treatment of the POA emissions, i.e., whether the POA in the emissions inventory is assumed to be the aerosol fraction before or after the evaporation of SVOC has occurred (Shrivastava et al., 2008;Murphy and Pandis, 2009).Air quality modeling studies based on traditional emission inventories (before evaporation of POA) have shown a modest increase in total predicted OA (from SI-SOA) in comparison to the traditional non-volatile POA approach, when the POA in the inventory is assumed to be the amount before evaporation of SVOC has taken place (Robinson et al., 2007;Shrivastava et al., 2008).Shrivastava et al. (2008) even reported a slightly decreased agreement between modeled and observed total OA concentrations throughout the continental US because of the evaporation of a large fraction of the POA mass.Although modest changes were obtained in predicted total OA concentrations, the fraction of oxidized OA (SOA) was substantially increased and the regional distribution of the SOA/OA ratio was more consistent with observations.More recently, Tsimpidi et al. (2010) have suggested that in the center of Mexico City SI-SOA contribute several µg/m 3 during the afternoon.Unlike Shrivastava et al. (2008), these authors assumed, based on comparisons with observations when using non-volatile POA, that the POA reported in the Mexico City Metropolitan Area (MCMA) emission inventory was the remaining aerosol fraction after evaporation of associated semi-volatile vapors.They determined the SVOC emissions by multiplying POA emissions by a factor of 3, and IVOC emissions by multiplying POA+SVOC by a factor of 1.5, as specified by Robinson et al. (2007).A box model study using data from the same campaign (MCMA-2003) (Dzepina et al., 2009) has reported a larger fraction of total SOA from SI-SOA during a stagnant day in which the local emissions underwent photochemistry for several hours with limited transport.By applying thermodynamic equilibrium with the parameters specified by Robinson et al. (2007) and the measured HOA concentrations, Dzepina et al. (2009) derived an amount of S/IVOC very similar to that of Tsimpidi et al. (2010), and achieved a much better agreement with observations during the morning and early afternoon hours in Mexico City and suggested that primary SI-SOA accounts for ∼55% of the observed SOA mass (∼ twice higher SI-SOA than V-SOA).The improved agreement between observed and predicted SOA when including SI-SOA in both studies also suggested that the POA emissions inventory in Mexico City does not include SVOC that evaporate after POA emission.
In addition to anthropogenic emissions (e.g. from diesel exhaust), a large fraction of POA emitted from biomass burning (BB) should also be considered as semi-volatile, as demonstrated by recent thermodenuder experiments (e.g.Grieshop et al., 2009;Huffman et al., 2009a, b).These SVOC should be oxidized quickly in the gas-phase, which is consistent with field studies that have observed production and/or oxidation of organic particles emitted from biomass burning (Reid et al., 2005;Capes et al., 2008;Yokelson et al., 2009).Grieshop et al. (2009) also showed that the SOA formed from BB woodstove emissions in their chamber could not be explained by V-SOA from the measured precursors (monoterpenes and aromatics) and proposed that SI-SOA represents most of the SOA formation potential from BB smoke.The initially proposed volatility distribution and aging mechanism of semi-volatile vapors (Robinson et al., 2007) has been recently updated based on comparisons of simulations with chamber experiments (Grieshop et al., 2009).Fire emissions are no longer assumed to contribute substantial IVOC in this updated approach.The lack of strong measurement-based constraints in the current treatment of POA emissions and chemistry, especially at long aging times beyond those that can be captured in smog chambers, could lead to substantial errors in predicted concentrations and spatial distributions of both POA and SOA, and therefore deserve further evaluation.
The overall objective of the present study is to provide an estimate of the potential importance of SI-SOA in the highly polluted environment of Mexico City.Specific objectives include: 1) quantifying the contribution of SI-SOA both within and downwind from Mexico City on much larger temporal and spatial scales than attempted by Dzepina et al. (2009) and Tsimpidi et al. (2010); 2) quantifying SOA production associated specifically with biomass burning activity, and 3) evaluating the model sensitivity to the S/IVOC parameterization by comparing two approaches, that of Robinson et al. (2007) and the recently revised parameterization of Grieshop et al. (2009).
We use SOA predictions from the meso-scale chemistry transport model CHIMERE in combination with measurements by the Aerosol Mass Spectrometers (AMS) from the MILAGRO 2006 field project, which were analyzed by the Positive Matrix Factorization evaluation toolkit (PMF, Ulbrich et al., 2009) to determine the contribution of different OA components.In addition to OA component concentrations, OA/ CO ratios (where CO is the CO concentration above the regional background), and O/C ratios will be used to evaluate the model.Comparisons are performed at the surface at two locations (urban T0 and suburban T1 sites) as well as over the city and its outflow along the G-1 and C-130 aircraft flight tracks (Kleinman et al., 2008;DeCarlo et al., 2008).Mexico City is a particularly suitable region for investigating SOA formation processes because of its highly active photochemistry and large amount of available SOA precursors.POA average levels are close to 7 µg/m 3 (Aiken et al., 2009a) and are expected to be associated with a large amount of primary S/IVOC vapors.Several aspects of OA formation in Mexico City have already been addressed in recent studies which have shown that the emissions of POA from both anthropogenic activities and biomass burning are reasonably represented in the emission inventories and models, with perhaps some underprediction of the urban POA, and mixed results for the BB POA with overprediction downwind of some large fires but underprediction in the urban area in the early morning (e.g.Fast et al., 2009;Hodzic et al., 2009;Zavala et al., 2009).A better prediction of the temporal trend of BBOA was achieved when the emissions profile took into account the emissions during the late evening into the shallow nighttime boundary layer (Aiken et al., 2009b).The fair agreement between HOA derived from measurements and predicted POA concentrations reported by Fast et al. (2009), Hodzic et al. (2009), and Tsimpidi et al. (2010) confirms that SVOC are unaccounted for in the POA emissions in Mexico City (i.e., that the inventory is consistent with the POA present after SVOC evaporation).Chemical analyses of submicron aerosols have revealed an active production of oxygenated organic species (OOA) within the city and their continuous production in the outflow region (e.g.Volkamer et al., 2006;DeCarlo et al., 2008;Kleinman et al., 2008;Aiken et al., 2008Aiken et al., , 2009a, b;, b;Herndon et al., 2008).V-SOA modeling showed once again its limitations in the context of Mexico City (Hodzic et al., 2009), where consistent with previous studies (Volkamer et al., 2006;Dzepina et al., 2009) it failed to reproduce the amount of OOA, including their production rate in the late morning, aging during the afternoon hours, and evolving background concentrations.Hodzic et al. (2009) estimated the model underprediction of the SOA mass to exceed 20 Tons during the afternoon peak over the Mexico City Metropolitan Area (MCMA).
Previous studies of SI-SOA in Mexico City have focused on SOA production within the Mexico City basin during a short period of time (Dzepina et al., 2009;Tsimpidi et al., 2010) (15 h and 4 days, respectively).Dzepina et al. ( 2009) also performed a 3-day box-model aging simulation that suggested that SI-SOA could increase substantially in the outflow from Mexico City.However the spatially and temporally limited framework of the previous studies did not allow investigating the significance of SI-SOA over regional scales and the multi-day processing of S/IVOC in the Mexico City plume in a 3-D context.Modeling results suggest strongly that the transport of OA from the regional scale is a substantial contributor to the observed SOA levels within Mexico City.Tsimpidi et al. (2010) assumed this SOA background to be between 3 and 7 µg/m 3 at the boundaries of their ∼150×150 km 2 model domain to achieve a better agreement with observations.Hodzic et al. (2009) showed that the regional background is partially due to biogenic SOA (∼1.5 µg/m 3 ) advected from the coastal mountain ranges.Regional POA emissions associated with surrounding populated areas (e.g.Cuernavaca, Puebla, Toluca) and biomass burning are also important and can contribute substantial amounts of S/IVOC vapors that, as they age, partition to the aerosol.Hodzic et al. (2009) reported that the V-SOA formed from biomass burning emissions with a traditional SOA module was very small.However it is likely that S-SOA from biomass burning may be the dominant source of SOA from these emissions (Grieshop et al., 2009) and may also contribute to the regional background.These unanswered questions provide motivation for reevaluating the significance of this formation process in a regional context over a two week time period.
The paper is organized as follows: in Sect.2, the CHIMERE model configuration and the updates to include S/IVOCs are presented; in Section 3, model predictions of secondary OA presented are assessed against OA measurements both within the city and in the outflow region.The contribution of anthropogenic and biomass burning SI-SOA is quantified and compared to the traditional V-SOA formation pathway.Section 4 summarizes the conclusions of this study.The terminology used for the various fractions and sources of organic compounds is given for reference in Table 1.

Model formulation and experiment design
For this study, the CHIMERE model is run from 11 to 31 March 2006 over the Mexico City region at both regional (35×35 km 2 ) and urban (5×5 km 2 ) scales using the same configuration and the same forcing (i.e.meteorology, emissions, boundary and initial conditions) as for our previous study presented by Hodzic et al. (2009).A detailed description of the model formulation and the aerosol treatment was given by Hodzic et al. (2009) and Bessagnet et al. (2008).Here, we only report modifications that were specifically implemented to the organic aerosol module (i.e.dynamic POA partitioning and oxidation) for the purpose of the present study.

Modeling the SOA formation from primary S/IVOCs
To examine the importance of SI-SOA formation, we implemented within the CHIMERE model an explicit treatment of emissions, gas-particle partitioning and gas-phase oxidation of primary organic aerosols as suggested by Robinson et al. (2007).The previous approach in which POA was considered inert and non-volatile has been replaced by a dynamic treatment in which POA material can evaporate, oxidize and re-condense over time to form S-SOA, and where IVOC are emitted in an amount proportional to the POA.The first step of this approach was to estimate the amount of semi-volatile and intermediate volatility vapors that are co-emitted with the modeled particulate POA.Following Robinson et al. (2007), this was achieved by assigning a volatility distribution of POA concentrations.S/IVOCs are distributed according to their volatility into nine surrogate species (with the lower volatility species referred to as SVOC and the three of higher volatility as IVOC; see Table 2) that partition between the gas and particle phases according to absorption into the total OA (with activity coefficient of 1 in all cases) under thermodynamic equilibrium (Pankow et al., 1994;Donahue et al., 2006).
Per the discussion above and similar to Dzepina et al. ( 2009) and Tsimpidi et al. (2010), we assume here that POA reported in the emission inventories, and therefore the modeled POA concentrations, correspond to the aerosol mass after evaporation of semi-volatile compounds.Therefore, this evaporated mass was added to the model in order to achieve the equilibrium with the emitted POA.According to Donahue et al. (2006) the fraction (X aero,i ) of particulate organic mass in equilibrium with S/IVOC is given by: where f i is the mass fraction of species i in the traditionally defined POA emissions, C i * is the effective saturation concentration of species i, C OM is the mass concentration of  a Molecular weight for all species for the ROB case is assumed to be 250 g/mol.
b The majority of intermediate volatile species (IVOC) is not accounted for in the traditional POA emission inventories.This additional contribution has been added here to SVOC 7 , SVOC 8 and SVOC 9 and represent 1.5 times the mass reported in POA emissions as in Robinson et al. (2007).c The gas-phase oxidation by OH of all of the nine S/IVOC is assigned a rate constant of 4×10 −11 cm 3 molec −1 s −1 (ROB simulation), and of 2×10 −11 cm 3 molec −1 s −1 (GRI simulation).
Following Tsimpidi et al. (2010), we assumed that approximately one third of the total emitted urban POA mass is in the aerosol phase under average ambient conditions found in Mexico-City (average OA concentration of 20 µg/m 3 ).Based on this assumption, total semivolatile plus aerosol phase emissions of organic species (POA+SVOC) were determined by multiplying the POA emission fluxes by a factor of 3. In addition to SVOC, an additional mass (also not present in the emissions inventory) has been added to account for emissions of higher-volatility vapors i.e.IVOC 7 , IVOC 8 and IVOC 9 .Following Robinson et al. (2007) this mass was assumed to be 1.5 times the mass of total gas and aerosol organic carbon emissions (POA+SVOC).This ratio of total IVOC to POA was experimentally determined, and it ranges from 1.5 to 3.0 for diesel engine emissions.Here we consider the conservative value of 1.5 for all emissions.Emission factors applied to each S/IVOC species are summarized in Table 2.Note that in presence of 10 µg/m 3 of ambient POA and very low SOA concentrations, the total amount of primary organic material (POA+SVOC+IVOC) derived from the Robinson parameters is ∼75 µg/m 3 at 20 • C, of which 87% remains in the gas-phase (SVOC + IVOC = 65 µg/m 3 ) (Table SI -3 of Dzepina et al., 2009).The extra amount of organic material introduced here (as S/IVOC) is small compared to the total organic gas + particle material in Mexico City, and is consistent with overall constraints from OH reactivity and integral FTIR stretches measured during MCMA-2003(Heald et al., 2008;Sheehy et al., 2008;Dzepina et al., 2009), but is poorly constrained.
All POA surface anthropogenic sources (e.g.traffic) as well as biomass burning emissions were considered as semivolatile following the S/IVOC distribution determined by Robinson et al. (2007) for diesel exhaust and wood burning.
POA emissions associated with elevated industrial sources (i.e. point sources) were not treated as semi-volatile (i.e. they were considered non-volatile with the same emission flux as in the inventory) in this study because their volatility distribution has not been reported.Contribution of point sources is relatively small in the MCMA region except for the Tula refinery complex located 45 km north of the Mexico City and the flights performed over the site have not detected significant organic aerosol emissions (in contrast to SO 2 for example, DeCarlo et al., 2008).Primary emissions and secondary trace gases and particulates were occasionally transported over the T0 and T1 sites during periods of northwesterly near-surface winds (de Foy et al., 2009).
The second modeling step consists in adding the photochemistry of S/IVOC species and their partitioning into the aerosol phase.
Gas-phase oxidation by OH of all of the nine S/IVOC is assigned a rate constant of 4×10 −11 cm 3 molec −1 s −1 (value used in the Robinson et al.Each reaction is assumed to reduce the volatility of the oxidized S/IVOCs by one order of magnitude (i.e.C* is divided by a factor of 10), with a small net increase in mass (7.5%) to account for added oxygen.Oxidized semi-volatile products from this reaction can partition to the OA phase to form SOA. The more the S/IVOC molecules become oxidized, the more likely they are to be in the aerosol phase.A continuous drop in volatility of the gas-phase material upon oxidation will result in an efficient partitioning to the aerosol phase of all the material, given sufficient time for oxidation.This may eventually lead to much higher net production of SOA downwind of the urban center (Dzepina et al., 2009).
In addition to the original parameterization of Robinson et al. (2007), here we also consider the updated version as described in Grieshop et al. (2009).Several updates have been considered compared to the original Robinson's approach, including: i) no biomass burning emissions of IVOC are considered; ii) each oxidation reaction is assumed to reduce the volatility of the oxidized S/IVOC vapors by two orders of magnitude (i.e.C* drops by 100), and a larger net increase in mass of 40% is assumed to account for the addition of oxygen.The addition of oxygen was parameterized based on comparison with laboratory measurements, although it is not clear mechanistically how so much oxygen could be added in a single reaction step; iii) the OH reaction rate constants, enthalpies of vaporization, and molecular weights of the surrogate species have also been updated (Table 2).
Similar to Hodzic et al. (2009), in this work dry deposition for the organic vapors is considered (Wesely, 1989).As the deposition velocities for these species have not yet been experimentally determined, deposition velocities have been calculated similar to that of NO 2 .The actual deposition velocities of the semivolatile and at least somewhat polar species that partition to SOA are likely to be larger than that of NO 2 , however this assumption provides a conservative upper limit to the amount of organic vapors.

Model application
For this study, two model simulations were carried out in order to calculate SI-SOA.The "ROB" simulation implements the parameterization of Robinson et al. (2007) in addition to V-SOA formation from anthropogenic, biomass burning and biogenic sources.The "GRI" simulation investigates the sensitivity of the SOA formation to changes in the assumptions following the updated parameterization of Grieshop et al. (2009) as described in Table 2.These two simulations are also compared to the OA simulation from Hodzic et al. (2009) that is based on the traditional SOA approach in which POA are considered as inert species.In this traditional simulation (referred as "REF" here, and as "BIO-T" by Hodzic et al., 2009) SOA is formed from anthropogenic, biomass burning and biogenic VOC precursors.Accounting for S/IVOC emissions and chemistry increased the computational cost of the model by about a factor of two compared to the REF simulation.
For comparison with measurements, the simulated parameters are spatially and temporally interpolated to the location of the measurement sites.The ground measurements are compared with the model's lowest level ( z ∼ 25 m).It should be noted that measured and predicted concentrations of all species are reported under ambient conditions of pressure and temperature unless indicated otherwise.

Results and discussion
In this section we examine the performance of the model in simulating POA and SOA in and near Mexico City.Predictions of SI-SOA (ROB and GRI simulations) are presented and evaluated against the highly time-resolved AMS measurements obtained at the surface within the city (T0 and T1 sites, Sect.3.2), as well as over the city and in the outflow region aboard the G-1 and C-130 aircraft (Sect.3.4).The model evaluation is limited to changes in organic species because CHIMERE predictions of meteorological variables, major gaseous pollutants, and primary and secondary organic aerosols using the traditional SOA treatment (i.e.REF simulation) has already been evaluated for March 2006.Simulated submicron OA fraction including the first four and half model size bins (from 0.04 to 1.0 µm) is used for comparison with the AMS/PMF-derived components: hydrocarbon-like organic aerosol (HOA), biomass burning organic aerosols (BBOA), and total oxygenated organic aerosol (OOA) (e.g.Aiken et al., 2009a).The modeled POA from anthropogenic (biomass burning) sources is compared with HOA (BBOA) respectively, while predicted SOA is assessed against OOA.A more detailed description on AMS apportionment using PMF analysis can be found in Ulbrich et al. (2009).Uncertainties associated with rotationality within sets of potential solutions of PMF analysis for the present dataset are documented in Aiken et al. (2009) and(DeCarlo et al., 2010).The implications of those uncertainties are much smaller than differences between the REF, ROB and GRI simulations.

Predicted OA concentrations in and around Mexico City
Figure 1 shows the spatial distribution of average TOA and SOA surface concentrations during 15-31 March 2006 as predicted by two S/IVOC (ROB and GRI) simulations and the traditional simulation.Predicted TOA mean concentrations feature similar spatial distribution in all three simulations with highest values associated with the Mexico City basin and surrounding smaller cities.Substantial differences are however obtained in terms of absolute TOA values.The GRI simulation predicts 60% higher peak TOA concentrations (∼22 µg/m 3 ) over Mexico City in comparison to the traditional simulation (<14 µg/m 3 ), while ROB predicts a more moderate increase of about 20% (∼18 µg/m 3 ).Background concentrations are also increased by a factor of 2-3 in the ROB and GRI simulations with respect to the REF simulation.The increase in predicted TOA levels by non-traditional simulations can be explained by a strong enhancement in the SOA production both near the emissions centers and in the outflow regions to the north and south of Mexico City basin as illustrated in Fig. 1b.Indeed, SVOC and IVOC vapors are co-emitted in equilibrium with POA mainly in the city center and are progressively transported from the city to areas located downwind.During their transport, semi-volatile vapors undergo several oxidation steps, which tend to increase their ability to partition to the aerosol phase as they age.Regional emissions of POA associated with smaller cities and the occurrence of fires also contribute to emissions of fresh S/IVOCs.This combination of local and regional emissions of semi-volatile vapors and their chemical processing leads to an increase in average SOA levels throughout the domain.
The GRI simulation features the strongest SOA production with average concentrations ranging from 12 to 15 µg/m 3 inside Mexico City, which is 3-5 times higher than those predicted with the REF model.Enhancement in the regional SOA concentrations is also seen and it is a result of the photochemical aging of S/IVOC vapors that have the potential to continuously produce SOA in the outflow region, and as some of the air advected away from Mexico City can be recirculated to the region (Emmons et al., 2010).Predicted SOA levels from the GRI simulation vary from 8 to 10 µg/m 3 in the outflow region northeast of the basin, and 6 µg/m 3 in the background areas, which is substantially higher than those predicted by the REF simulation.The increase in SOA production appears to be more moderate in the ROB simulation with SOA levels close to 8 µg/m 3 in the near field, 6 µg/m 3 in the outflow region and ∼3 µg/m 3 at remote locations.This difference of 30-50% in the predicted SOA levels inside the city is related to the differences in the treatment of the chemistry of S/IVOC vapors and will be further discussed in the next section.
Including the chemistry and partitioning of semi-volatile primary organic vapors leads to substantially higher and likely more realistic predictions of SOA concentrations in comparison to the traditional simulation (Hodzic et al., 2009).Predicted SOA now represents a much larger fraction (40-60%) of the total OA, which is more consistent with the observed values reported during the campaign (e.g.Aiken et al., 2009a, DeCarlo et al., 2009;Herndon et al., 2008) that are described in more detail next.1240 Fig. 2. Time series of modeled and observed surface concentrations of major carbonaceous aerosol compounds including primary organic aerosol from anthropogenic (HOA) and biomass burning emissions (BBOA), secondary organic matter (OOA), and total organic aerosols (TOA).Time series and average diurnal profile plots compare the ROB (red line) and GRI (green line) simulations, with the AMS measurements at T0 (given at ambient atmospheric conditions).The variability associated with average observations and model predictions is given in shaded area and red vertical bars, respectively, in the right panels.Diurnal profiles for the traditional simulation are also indicated for POA, SOA and TOA plots (blue line), where POA includes both anthropogenic and biomass burning contributions.

Primary Organic Aerosols
In both non-traditional model simulations (ROB and GRI) the emitted POA is considered as a semi-volatile species.Therefore, POA can evaporate from the particle and undergo further chemical processing in the gas phase leading to production of lower volatility organic oxygenated vapors that are likely to re-condense and form SI-SOA.This dynamic treatment of POA can have major consequences on the predicted amount of POA and its effects need to be examined first.
Figures 2 and 3 compare the time series of observed primary organic aerosols from anthropogenic sources (HOA) and biomass burning (BBOA) with the corresponding mod-eled primary species at the T0 and T1 sites, respectively.At T0, the model reproduces fairly well the magnitude and the diurnal variations of the observed HOA for both ROB and GRI simulations (Fig. 2a) as suggested by a relatively low bias (<20%) and correlation coefficients greater than 0.55. the temporal agreement with observations and the correlation coefficients.Small differences exist between the ROB and GRI predictions, with a slightly higher modeled HOA for the GRI simulation (<2% for the average concentration), which is consistent with the fact that a larger total organic mass is available for their partitioning (Dzepina et al., 2009).Overall, the model skill in reproducing HOA has not been deteriorated or considerably changed by the introduction of the dynamic partitioning of primary organic species as shown in Fig. 2, indicating that the treatment of anthropogenic emissions in the model is reasonable in the near field.POA concentrations predicted in the ROB and GRI simulations stay close to the REF model, with however a 1 to 2 µg/m 3 decrease during afternoon hours (12:00-18:00 p.m. LT) which leads to a better agreement with observations.
At the T1 site at the edge of the city (Fig. 3a), discrepancies between observed and predicted HOA concentrations are larger both in terms of magnitude and the temporal variability as already found for the traditional simulation (Hodzic et al., 2009).The model simulations show an overprediction (∼30-50%) which is mainly caused by a factor of 2-3 overprediction of nighttime values as suggested by the comparison of diurnal profiles.Although the occurrence and the magnitude of the morning HOA peak associated with the traffic emissions is well captured, the correlation coefficient stays low for both non-traditional simulations (∼0.3) due to nighttime inconsistencies.POA concentrations predicted in the ROB and GRI simulations are lower by ∼1 µg/m 3 then the REF model suggesting that the volatility may be a bit too high in the model causing excessive POA evaporation downwind of the city.As for T0, the dynamic treatment of primary organic species and assumptions made for their emissions allow reasonable prediction of the anthropogenic HOA.
The predicted BBOA shows larger discrepancies with the observations.At T0 (Fig. 2b), BBOA is present throughout the period of interest (average concentration of 1.54 µg/m 3 ), with however higher concentrations from 15-24 March and occurrence of three large fire events on 18, 20 and 21 March characterized by peak values in the range 10-25 µg/m 3 .Both simulations generally underpredict the amount of smoke  Hodzic et al., 2009); and for two S/IVOC model simulations that in addition of traditional precursors includes also the semi-volatile primary organics following the approach of Robinson et al. (2007), or Grieshop et al. (2009) POA reaching the T0 site.While the models capture the afternoon increase in BBOA concentrations, they fail to reproduce the early morning (∼06:00 LT) peaks that occurred on the 18 and 21 March in a very shallow boundary layer.As already discussed in Hodzic et al. (2009), reproducing the nighttime concentrations of pollutants is very challenging as pollutant dispersion is highly sensitive to small errors in boundary layer mixing and near-surface winds.Although the gap between the model and observations seems large (here ∼3 µg/m 3 ), it corresponds to relatively small errors in terms of emitted mass of smoke.The smoldering regime in the fire lifecycle that usually occurs in the early morning and can produce large quantities of smoke is not well depicted in the emissions.In our simulations, diurnal profiles with maximum emissions in the early afternoon (12:00-16:00 LT) are applied.Aiken et al. (2009b) achieved a better correlation between the BBOA at T0 and model prediction by using a modified emissions profile which includes smoldering emissions into the late evening and night.It is also likely that some fire sources are missing from the biomass burning emission inventory as already discussed by Fast et al. (2009).Fast et al. (2009) presented a detailed comparison of simulated non-volatile BBOA (generated with the same inventory used here) with the same set of AMS observations, and have found similar BBOA underprediction at T0.At the T1 site (Fig. 3b), similar difficulties are encountered.Predictions of both model simulations are substantially below the observed levels with an underprediction of ∼70% of observed BBOA, largely associated with the early morning discrepancies.A somewhat better agreement is reached in the late afternoon, although the temporal variability is not captured (r 2 ∼ 0.1).

Secondary Organic Aerosols
We next examine the model skill in predicting secondary organic aerosols at the surface and the role of S/IVOC vapors in the formation of SOA.According to our earlier study (Hodzic et al., 2009) SOA production from anthropogenic and biogenic precursors based on the traditional modeling approach could explain less than 35% (see Table 3) of OOA within the city during the MILAGRO experiment (15% attributed to anthropogenic precursors).The comparison of observed and predicted SOA time series shown in Figs.2c and 3c suggest a substantial enhancement in SOA concentrations at both urban and suburban locations when SI-SOA is taken into account.At T0 (Fig. 2c), the ROB simulation captures the overall magnitude of observations reasonably well with an increase in average modeled SOA from 2.6 µg/m 3 (REF simulation) to 6.3 µg/m 3 for ROB which is in a better agreement with the observed value of 7.9 µg/m 3 (Table 3).Besides the bias reduction, the appreciable decrease of the RMS error (from ∼7.4 to 5.25 µg/m 3 ) confirms the improved agreement with observations.The increase in correlation coefficient from 0.38 to 0.44 also denotes a more consistent temporal agreement between model and observations.The ROB model captures within a factor of two the photochemical production of SOA in the late morning (∼09:00-11:00 a.m.LT) and the production is maintained despite the growth of the convective boundary layer that was causing SOA to evaporate in the traditional model (Hodzic et al., 2009).If SI-SOA evaporation occurs in the ROB model during this photochemically active period, the vapors can be quickly oxidized and re-partition to the particle phase (Dzepina et al., 2009).The predicted nighttime values are about 4 µg/m 3 and are consistent with OOA.The ROB model predicts somewhat lower production of SOA during the day near the source www.atmos-chem-phys.net/10/5491/2010/Atmos.Chem.Phys., 10, 5491-5514, 2010 region leading to underpredicted SOA peak values, in particular after 24 March by a factor of 2. The passage of a cold surge led to precipitation on the 24 March and suppression of biomass burning, so that the underestimation by the ROB model stems directly from the representation of urban SOA production.At the T1 site, strong enhancement in predicted SOA concentrations is obtained with an increase in average levels from 1.1 µg/m 3 (REF simulation) to 4.2 µg/m 3 for the ROB simulation, leading to a fairly good agreement with observations particularly in terms of absolute values (the model bias ∼12%).This improvement can also be seen in terms of enhanced temporal agreement with the increase of the correlation coefficient from 0.24 to 0.28.The GRI parameterization tends to produce substantially higher SOA mass than ROB both in the immediate vicinity of sources and at the edge of the city.At the T0 site (Fig. 2c), 80% higher average SOA levels are computed (∼11.4 µg/m 3 ).The late morning SOA enhancement (09:00-11:00 a.m.LT) reaches 12 µg/m 3 in the model which is 30% higher than observed (∼9 µg/m 3 ).The model has the tendency to overpredict the observed levels on 16, 17, and 21 March, while a much better correspondence with OOA levels is achieved after 24 March.This overprediction coincides with periods of high atmospheric stability during which the wind speed was underpredicted in the model (and close to zero) leading to an unrealistic accumulation of pollutants during the day, which was also visible in predicted levels of CO (Hodzic et al., 2009).
A larger overprediction by the GRI model can be noted at the T1 site (Fig. 3c) almost everyday.The model predicts 70% higher SOA levels.The bias and RMS error exceed 3.4 and 5.7 µg/m 3 , respectively.In particular, the photochemical production from 09:00 a.m.-01:00 p.m. LT is far larger than observed (6 µg/m 3 for predicted SOA increase vs. 1 µg/m 3 increase for OOA).The correlation coefficient is also decreased in comparison to the ROB simulation.The results suggest that the updated treatment of S/IVOC aging used in the GRI model leads to unrealistically high SOA production downwind of the source region, consistent with the results of Dzepina et al. (2009).The likely reason for this enhancement in GRI SOA production in comparison to ROB is the addition of 40% oxygen in mass and the drop of two orders of magnitude in volatility upon each oxidation step (instead of 7.5%, and 1 order of magnitude for ROB) that facilitate the transfer to the aerosol phase.Our findings are consistent with results of a box model study of Dzepina et al. (2009) in which it was shown that the GRI approach leads to 62% higher SOA production compared to ROB.
It should also be noted that ROB and GRI treatments involve several additional OH reactions and are likely to consume some OH in comparison to the REF simulation (Hodzic et al., 2009).The OH mid-day average levels were found to decrease from 0.42 ppt (REF) down to 0.30 ppt (GRI) and 0.25 ppt (ROB), which led to a better agreement with the observed mid-day values of ∼0.30 ppt.Nighttime levels were still underpredicted.OH is the only reactive gas directly susceptible to modification by the addition of SIVOC gas-phase reactions, since these species are assumed to not react with O 3 or NO 3 .

Total Organic Aerosols
Comparisons of predicted total OA suggest that the implementation of the SI-SOA mechanisms improves model ability to predict SOA formation.At the urban T0 site (Fig. 2d), predicted TOA lies in the variability range associated with the observations: both ROB and GRI simulations show lower absolute bias (<20%) than the traditional simulation (∼45%).No substantial improvement is however seen in simulated TOA temporal variability, which correlates slightly better for ROB (R 2 =0.48) than for GRI (R 2 =0.46) and REF (R 2 =0.45) simulations.As suggested by diurnal profiles (Fig. 2), this temporal mismatch is likely caused by the dominance of the nighttime errors featuring too high POA concentrations , due to a too low model mixing layer (Hodzic et al., 2009;Fast et al., 2009).Model predictions correlate somewhat better with daytime (06:00-18:00 LT) observations of TOA for all three simulations i.e.ROB (R 2 =0.63),GRI (R 2 =0.52) and REF (R 2 =0.58).
At the T1 site (Fig. 3d), the overprediction by GRI is obvious, with the bias in simulated TOA approaching 41%, and a drop in correlation coefficients down to 0.11.The excessive photochemical SI-SOA production in the early afternoon (12:00-15:00 p.m. LT) leads to almost a factor of 1.7 higher TOA levels than observed.The ROB predictions match the observed levels within 20%, although the temporal variability is not well captured (R 2 =0.23).Similar to Hodzic et al. (2009), nighttime comparisons appear more challenging at both sites as they are influenced by nighttime boundary layer effects, errors in the emissions or unaccounted oxidation processes.

Relative contributions to predicted SOA
The relative contribution of various sources to simulated SOA concentrations is shown in Fig. 4, including V-SOA from anthropogenic and biogenic precursors, and SI-SOA from anthropogenic and from biomass burning sources.For the ROB simulation, traditional and non-traditional precursors seem to equally contribute to the SOA budget in the source region (T0) during the 15-31 March period (Fig. 4a).During the day, about 2 µg/m 3 can be attributed to each anthropogenic and biogenic V-SOA, while about 3.5 µg/m 3 is due to anthropogenic SI-SOA and 1 µg/m 3 to biomass burning SI-SOA.The ratio of Ant.SI-SOA/V-SOA is similar to that found by Dzepina et al. (2009) and substantially higher than of Tsimpidi et al. (2010).The differences are probably explained by the fact that Tsimpidi et al. (2010) used very aggressive aging of V-SOA for which there is limited experimental evidence, but which greatly increases V-SOA formation.Such aging of V-SOA was not implemented here or by Dzepina et al. (2009).The relative contribution of SI-SOA increases substantially at the city edge (Fig. 4b) as vapors are continuously being oxidized and driven into the particle phase by the multi-generation oxidation mechanism.The anthropogenic SI-SOA stays elevated (i.e.close to 2.5 µg/m 3 ), and exceeds the anthropogenic V-SOA (<0.5 µg/m 3 ), or biogenic V-SOA (1.2 µg/m 3 ).SI-SOA from biomass burning contributes about 0.3 µg/m 3 throughout the day, which is of comparable magnitude to the contribution of anthropogenic V-SOA.A very different behavior is observed for the GRI simulation shown in Fig. 5, in which anthropogenic SI-SOA dominates the total SOA concentrations at both sites.Anthropogenic SI-SOA increase by 3-5 µg/m 3 at T0 (Fig. 5a) and 2-3 µg/m 3 at T1 (Fig. 5b) in comparison to the ROB simulation, while contributions of the other sources remain at similar levels.

O/C ratios
In addition to concentrations of organic aerosol species, AMS data provide the opportunity for evaluating the predicted oxygen to carbon (O/C) atomic ratio for the total OA.The measured O/C ratio of total OA is a good indicator of the level of processing and therefore the age of the organic aerosol (Aiken et al., 2008;DeCarlo et al., 2008;Jimenez et al., 2009;Ng et al., 2010;Heald et al., 2010).Typically, the values range from very low (<0.1) for e.g.vehicular emissions, to values close to 0.35-0.4 in freshly formed SOA, and up to 0.8-1 for aged air masses that contain more processed organic material.Therefore the O/C ratio provides an additional constraint for evaluating SOA formation mechanisms.Model O/C ratios were estimated using measured chamber values for V-SOA (Ant V-SOA=0.37;Bio V-SOA=0.40),hydrocarbon-like organic aerosol (Ant POA=0.06), and biomass burning organic aerosol (BB POA=0.30) reported by Aiken et al. (2008) and Shilling et al. (2009).The organic carbon content in the secondary organic matter formed from the S/IVOC vapors was directly determined from modeled species as both carbon and oxygen fractions are explicitly modeled.This difference can also be seen in Fig. 6c that compares the simulated O/C values with estimated O/C values derived from AMS measurements at T0.The observed ratios range from 0.2 to 0.7 and display a strong diurnal amplitude with lower values in the morning (∼06:00 a.m.) during the rush hour and in absence of photochemistry, and higher values during the day (∼02-04:00 p.m.) when more secondary oxidized matter is being formed by photochemistry.The O/C ratios calculated from the ROB treatment are substantially lower than the observed levels, with values ranging from 0.08 in the morning to 0.2 in the early afternoon.The predicted SOA seems to be much less oxidized than the measured OOA, which is indicative of a too low addition of oxygen during oxidation in the ROB parameterization.A much better agreement with the observations is achieved for the GRI treatment in terms of magnitude and the temporal vari-ability, although the predictions were consistently somewhat lower than observed after 24 March.The reason for this large difference between ROB and GRI models arises from the assumptions concerning the mass of oxygen that is added at each oxidation stage.Indeed, in ROB it is assumed that each generation of oxidation of S/IVOC vapors adds 7.5% to the mass of products due to the addition of oxygen atoms.This percentage is much larger for GRI, i.e. 40%.However the large addition of oxygen upon each oxidation step in the GRI parameterization has not been mechanistically explained.It is also possible that other mechanisms that are not considered in the present simulations, such as aerosol reactive uptake of very oxygenated such as glyoxal (Volkamer et al., 2007(Volkamer et al., , 2009) ) is causing some of the rapid increase of oxygen in reality.Although the ROB treatment seems to perform better in terms of predicted surface concentrations, it results in far too low predicted O/C in OOA.This comparison show that the amount of oxygen in OA is an important constraint for determining the nature and the age of the aerosol, parameters that will play a key role in predicting the radiative forcing of organics in regional and climate models.

Model evaluation along G-1 and C-130 flight paths
S/IVOC vapors can undergo several generations of oxidation and have the potential to continuously contribute to the production of SI-SOA further downwind of the city.AMS measurements from aboard the G-1 aircraft (Kleinman et al., 2008) and the C-130 aircraft (DeCarlo et al., 2008) during MILAGRO are available study the evolution and growth of organic aerosols aloft and downwind of Mexico City.Data from 8 G-1 flights are compared with the model outputs (including flights from 15a, b, 18a, 19a, b, 20a, b, 22a, where a and b refers to morning and afternoon flights).For the C-130, PMF-analyzed AMS data for only one flight (29 March) are available for comparison with the model.Figures 7 and 8 compare CO mixing ratios and organic aerosol species concentrations as predicted by ROB and GRI treatments and observed during two particular flights.It should be noted that G1 AMS data have been recently updated which resulted in a slightly different repartition between HOA and OOA aerosol components from the one reported in Fast et al. (2009) (based on 3-component PMF analysis).The updated data are used here.We also note that PMF performs better when using high-resolution AMS data instead of unit-resolution data, as the spectra of the different components are more distinct and easier to separate under high-resolution (DeCarlo et al., 2006;Aiken et al., 2009a).
On the morning flight of 19 March (Fig. 7), the aircraft arrived from the eastern part of the domain, crossed over the T1 site and performed a loop over Mexico City passing over T0 at 17:00-17:20 UTC, before sampling the outflow region north of the city around 18:00 UTC.As suggested by model OA predictions, the aircraft crossed several fire plumes that originated from the hills surrounding the city and were transported northward in the southerly flow.The comparison of observed and predicted CO mixing ratios shows the model skill in reproducing the magnitude and spatial variability observed in CO fields as indicated by a low overall bias (−50 ppbv) and the correlation coefficient of 0.6.The highest concentrations are observed at 17:00 UTC over the city with the predicted values close to 450 ppbv which is somewhat below the observed peak of 500 ppbv.This agreement for a chemically inert species such as CO indicates that the plume location and the transport of the pollutants are reasonably predicted on this particular day.
In contrast, modeling of OA appears to be more challenging and substantial differences can be seen between the ROB and GRI treatments.Although the observed total OA is reasonably well reproduced, with an overall underprediction of about 30% for the ROB simulation and almost no bias for the GRI simulation (∼7%), the spatial variability is not well captured as indicated by low correlation coefficients (∼0.16).The highest discrepancies are seen in the outflow region north of T1 (∼17:50 UTC) and coincide with the un-derpredicted CO mixing ratios.This model error has likely a dynamic origin as explained by Fast et al. (2009).Relatively small errors in the wind speed and wind direction are likely to occur in this complex terrain as already reported in Fast et al. (2009) and Hodzic et al. (2009), and can cause plume misplacement of several km or the too rapid dispersion of the city plume.In terms of OA constituents, major inconsistencies are found for the BBOA fraction.Two out of four smoke plumes are either mislocated in the model or the amount of emitted smoke particles is highly overpredicted in the model estimates.The predicted smoke concentrations west of the T0 site and east of the T1 site both reach 4-6 µg/m 3 while these plumes are not reported in the observations.The comparison suggests that the plume along 98 • 40 N is slightly shifted towards east and that concentrations in the plume center are lower.This inconsistency in BBOA predictions was also seen by Fast et al. (2009) that reported smoke concentrations above 15 µg/m 3 in the WRF/Chem model.It is likely the result of an inaccurate representation of the fire emissions (e.g.fuel types and consumption rates, flaming vs. smoldering regimes, etc.) which is not well represented in current fire emission inventories.Also, at this stage, rotational ambiguity in the unit mass resolution PMF results for BBOA and HOA cannot be discounted.
In addition, the modeled anthropogenic POA fraction for both the ROB and GRI simulations appears to be substantially lower than the observed HOA.The predicted average POA values (∼0.35 µg/m 3 ) only account for about 30% of the observed HOA (∼0.9 µg/m 3 ).This model underprediction seems equally important within the city boundary layer and in the downwind region.In comparison with the REF case, the model bias is increased from −0.1 to −0.55 µg/m 3 when the dynamic POA partitioning is considered.As already discussed in Sect.3.2, the comparison with observations suggests that atmospheric HOA is less volatile than what is being assumed in the ROB and GRI approaches.It is likely that the simulated POA evaporates too steadily upon dilution of the urban plume in the growing PBL.Once evaporated, the primary organic vapors are oxidized into less volatile species that are likely to re-condense and form SOA. It should also be mentioned that our results for the REF case (POA considered as inert) show a better agreement with observations than the ones reported by Fast et al. (2009) due to changes in HOA estimates in the more recent PMF results.According to their comparison results, a large model POA underprediction was a common feature to all G1 flights, and the underprediction could reach up to a factor of 4 (i.e. the morning flight of 18 March).In our study the POA underprediction occurs but is generally more moderate.
Interestingly the best agreement is found for the predicted secondary organic fraction, with a good correspondence with OOA over the city as well as north of T1.Both ROB and GRI treatments are somewhat lower in terms of average predicted levels with a negative bias of 0.9 and 0.5 µg/m 3 respectively for the observed average OOA of ∼2.8 µg/m 3 .Similar to  POA particles, the dilution of the plume also has a tendency to evaporate the freshly formed SOA.However this evaporation competes with the continuous formation of the new SOA downwind leading to a net production of the SOA mass as reported by Kleinman et al. (2008).The SOA concentrations predicted by the traditional SOA treatment are also displayed in Fig. 7 as a blue line.It can be clearly seen that the SOA production aloft over the city and downwind is insufficient when only traditional mostly aromatic precursors are considered, and is much improved when SI-SOA is added.
Another flight of interest to analyze SOA formation took place in the afternoon on 20 March.The airplane sampled exclusively the city outflow region as shown in Fig. 8, with three consecutive passes over the T1 site between 21:00 and 21:30 UTC and three passes north of T1 (the T2 site not shown here).Model simulations suggest the presence of a long narrow plume transported north of the city in the southwesterly flow.The outflow appears to be a mixture of the anthropogenic pollution originated from the city and the biomass burning emissions that took place on the hills southwest of T0.The six consecutive airplane encounters with the plume can clearly be identified in observed and predicted CO mixing ratios.The model was able to approximately predict the locations and absolute enhancements of the CO levels associated with the main plume most of the time, as well as the background levels of CO.The first and the last airplane overpasses were however performed at higher altitudes (between 4.5 and 5 km above ground level vs. 3-4 km for other periods) and suggest the model underprediction of the vertical extent of the plume.Concentrations at these higher altitudes are more challenging to simulate accurately due to the coarser model resolution at higher altitudes.In addition, the spatial structure of the plume is not well captured north of T1.A relatively small misplacement of the plume toward west can be seen, and is responsible for a reduced correlation with the observations (∼0.42) during this flight.Similar to CO mixing ratios, the predicted OA concentrations match reasonably well the observed levels within the plume, except for the elevated air mass at 21:00 UTC.The transition between elevated OA levels of 6-8 µg/m 3 found within the plume and the background concentrations of a few µg/m 3 seems well reproduced in the model.Up to 1.5 µg/m 3 higher OA levels are predicted with the GRI model compared to the ROB one.
The comparison of individual OA components suggests the occurrence of error compensation between a high SOA production, and low levels of anthropogenic POA.In particular, the SOA production from the GRI treatment is above (∼0.5 µg/m 3 ) OOA concentrations within the plume.However, the overall average OOA concentrations are captured well as suggested by relatively low bias (<50%) for both simulations.Comparing these SOA predictions from the traditional treatment shows that the S/IVOC chemistry is responsible for this increased agreement with the observations.On the other hand, the magnitude of the predicted anthropogenic POA is substantially lower than the observed HOA within the plume for both simulations, although somewhat higher concentrations are simulated inside the plume.Discrepancies can also be seen between simulated and observed BBOA, with more pronounced gradients in the model between the plume area (0.4 to 1.6 µg/m 3 ) and the remote background area that is not influenced by smoke.The observed BBOA levels show more uniform spatial structures, with higher background levels i.e. 0.4-0.6,indicating that greater dispersion of the smoke plume during this day than simulated by the model.It should be noted that our model results for primary organic aerosols generated from both anthropogenic and biomass burning sources agree well with WRF/Chem analysis of primary aerosol compounds reported by Fast et al. (2009).
A relatively similar model behavior with respect to OA formation is found for the C130 flight of 29 March (Fig. 9).This day is characterized by more stable atmospheric conditions leading to a greater accumulation of pollutants within the city basin than observed for the previous two flights (19)(20).The pollution plume is located slightly west and northwest of T0.Modeled concentrations display a fairly good spatial agreement with measurements along the aircraft track for all species (e.g.CO: R 2 =0.86), with exception of BBOA.The observed magnitude and temporal variability in  CO concentrations is captured remarkably well.The predicted SOA levels for both S/IVOC runs match within 50% the observed concentrations (5-12 µg/m 3 ) along the first aircraft leg (∼18:30-19:30 UTC) north of the city.High OOA levels found within the urban plume (∼15-20 µg/m 3 ) are reasonably reproduced by the GRI simulation, and within 30% for the ROB one.As for previous aircraft comparisons, the HOA levels are underpredicted by a factor of 2 by the model over the urban area, leading to an overall underprediction of observed OA concentrations reaching 50% for ROB and 40% for GRI respectively.
Finally, scatter plots shown in Fig. 10 summarize the comparisons of observed and simulated organic aerosol concentrations for the ensemble of G-1 flights considered in this study.The distinction is made between points measured aloft the MCMA area and in the remote regions.Most of the observed and simulated OA concentrations stay below 12 µg/m 3 , and no systematic bias can be seen in model simulations.The higher OA levels (>12 µg/m 3 ) are observed over    < 25 ppb the city, but no clear difference in model performance can be noted between the city and the remote regions.The strong scatter around the 1:1 slope is likely associated with inaccuracies found in the predicted spatial features of the plume (i.e.plume extent and location), and therefore are related to dynamic processes rather than aerosol chemistry.To examine independently the SOA production in the plume from the transport (eliminate the sensitivity to the advection and mixing processes), in the next section we will consider the OA concentrations normalized to CO.

Continuous SOA growth in the aged plume
The comparison with aircraft data allows quantifying the changes in the OA concentrations as a function of photochemical age of the urban plume.Similar to Kleinman et al. (2008), the SOA production downwind of the city was estimated from the OA/ CO ratio (where CO is CO minus its assumed 100 ppb background concentration) and the age of the plume was determined by −Log 10 (NO x /NO y ).CO and organic aerosols concentrations display very different age dependences.CO levels decrease as the urban plume is transported downwind of the emission sources upon dilution, tending towards the background levels.Chemically long-lived, CO can be considered as a reasonably conserved tracer for the urban emissions and their dispersion on shorter timescales.Conversely, OA concentrations decrease upon dilution at the same rate as CO, but also increase due to the continuous photochemical production of SOA.Therefore, concentrations of OA normalized to CO mixing ratios are a good indicator of the SOA growth as a function of the chemical age of the plume.
Figure 11a shows the OA/ CO ratio for Standard Pressure and Temperature (STP) conditions (273 K and 1 atm).STP conditions were used here to account for changes in OA concentrations due to changes in altitude of the airplane (i.e.changes in volume).Similar to Kleinman et al. (2008), comparison was performed only for aircraft points that satisfy the conditions summarized in Table 4.It can be seen that the observed OA/ CO ratio increases substantially as the plume ages as already observed in the outflow of several   urban areas (de Gouw and Jimenez, 2009).The values increase from 20-30 µg sm −3 ppm −1 in the fresh air mass up to 60-70 µg sm −3 ppm −1 in a 1-day aged plume.Ratios predicted by the ROB treatment are in a reasonably good agreement with the observations, except for the very last point that seems overpredicted by nearly 20 µg sm −3 ppm −1 .The GRI treatment, however, overpredicts the OA growth rate (i.e.∼30%) with values ranging from 20-40 µg sm −3 ppm −1 to 70-95 µg sm −3 ppm −1 .This overprediction for the GRI treatment is consistent with the results already discussed at the surface.(MILAGRO 2006) and reported in recent studies.MCMA 2003 summarizes results reported by (Volkamer et al., 2006;Dzepina et al., 2009;and Tsimpidi et al., 2010); CARBOSOL 2006 refers to (Simpson et al., 2007);TORCH 2003refers to (Johnson et al., 2006);NEAQS 2003(de Gouw et al., 2005); ACD-Asia 2001 (Heald et al., 2005).In blue are represented estimates of the traditional SOA approach, while red refers to the volatility basis set.
OA/ CO ratios are, however, influenced by model excessive loss of POA by evaporation during the transport of the plume, and therefore it is difficult to appreciate the aerosol growth associated with the photochemistry.The comparison of OOA/ CO ratios reveals an excessive production of SOA downwind of the city mainly for the GRI parameterization as shown in Fig. 11b.The ROB simulation captures relatively well the SOA growth downwind of the city.The SOA/ CO ratios are in range 30-60 µg sm −3 ppm −1 and 40-80 µg sm −3 ppm −1 for ROB and GRI simulations respectively, while the observed values suggest a more limited increase of 30-50 µg sm −3 in SOA per ppm of CO in the 1/2day to 1-day aged plume.The SOA growth from the updated GRI treatment is 30-40% higher on average mainly due to larger decrease of volatility and higher oxygen gain at each oxidation step of S/IVOC vapors.Similar values (∼62%) are reported by a box model study of Dzepina et al. (2009).This active production of SOA downwind of Mexico City contrast greatly with the severe underprediction of the SOA production simulated by the traditional SOA model (Hodzic et al., 2009) in which SOA was below 20 µg sm −3 ppm −1 and was mainly associated with the presence of biogenic background SOA concentrations.

Conclusions
In this study a meso-scale chemistry-transport model has been applied to investigate whether or not semi-volatile and intermediate volatility primary organic vapors contribute to the formation of SOA in the vicinity of Mexico City.The model has been updated to account for the dynamic treatment of POA emissions which includes both the gas/particle partitioning of primary species and the photochemical processing of associated vapors.Two approaches based on Robinson et al. (2007) and Grieshop et al. (2009) have been applied and inter-compared.The following are the principal findings of this study: i) Model results suggest that organic aerosol formation from primary organic vapors has the potential to increase SOA levels substantially both within Mexico City where POA emissions are the highest, and in the surrounding regions that are influenced by the transport and ageing of S/IVOC vapors.The ROB and GRI parameterizations produce somewhat different SOA concentrations and properties due to differences in the assumptions associated with volatility, reaction rates and oxygen gain at each oxidation step.
ii) In both treatments an enhancement in average SOA levels is predicted within the city in presence of S/IVOC relative to the traditional mechanism (3.8 µg/m 3 or a factor of 2.5 for ROB, and 8.7 µg/m 3 or a factor of 4 for GRI), which substantially reduces the gap with OOA observations.The diurnal variability characterized by the rapid SOA increase after sunrise is reproduced leading to higher temporal correlations for both simulations.
For the ROB case, the predicted SOA lies slightly below OOA values but within 20% and 13% on average at the T0 and T1 sites.Although the ROB case seems to capture relatively well the amount and the evolution of SOA in Mexico City, the predicted average O/C ratio of total OA is 3 times lower than the observed AMS value.The updated GRI treatment produces a much better agreement for O/C ratios at the surface in terms of their absolute values and temporal variability, but the SOA production is generally overestimated.The predicted SOA mass for the GRI case is higher compared to the OOA measurements, as indicated by an overprediction of 45% at T0 and 70% at T1.
iii) Downwind SOA production persists actively due to continued photochemistry of multi-generation S/IVOC oxidation products.Similar to aircraft observations, the predicted OA/ CO ratio for the ROB case increases from 20-30 µg sm −3 ppm −1 up to 60-70 µg sm −3 ppm −1 between a fresh and 1-day aged air mass, while the GRI case produces a 20-30% higher OA growth than observed.Comparison with individual flights of 19, 20 and 29 March indicates a much improved prediction of SOA in comparison to the traditional approach (which was a factor of 5-10 too low) and an agreement generally within 30% with observations (see Fig. 12).These results indicate that this new framework proposed by Robinson et al. (2007) based on the volatility of primary organic species constitutes a genuine revision in our understanding of SOA formation and is likely to be one of the important SOA formation pathways, although details of the parameterization need further improvement.
iv) The dynamic treatment of POA does not have major effects on its surface concentrations in the source region, and the model results remain close to ones obtained with the non-volatile POA assumption.However, both non-traditional treatments show a deficiency in regard to POA ageing with a tendency to over-evaporate POA upon dilution of the urban plume during the regional transport.Comparisons with measurements suggest that atmospheric HOA may be less volatile than assumed in these parameterizations and the sensitivity of SOA production to this parameter should be investigated in future studies.
v) This study also highlights the need for improvements in current parameterizations.Experimentally constrained amounts and partitioning/ageing parameterizations of S/IVOC are urgently needed for further progress in this area.In particular, the effect of fragmentation of oxygenated species leading to higher volatility species (Kroll et al., 2009) needs to be taken into account.
Parameterizations which explicitly track both volatility and oxidation state (e.g.Jimenez et al., 2009) may also allow improved predictions of OA evolution and aging.Opportunities for improvement also lie in the better constrained deposition velocities of the S/IVOC species.Their values are likely higher than the ones considered in this study (based on NO 2 , Hodzic et al., 2009) and should be experimentally and/or theoretically constrained.Finally, we also note that the modeling of SOA is still underconstrained.The improved agreement with observations achieved here does not prove that S/IVOC are the major missing SOA source in megacity environments, and other possible pathways such as formation from very volatile species like glyoxal need to be explored in future studies.

Fig. 1 .
Fig. 1.Spatial distribution of secondary organic (a) and total organic (b) aerosols concentrations (µg/m 3 ) in PM 2.5 particles over the Mexico City basin as predicted by ROB (left panels), GRI (middle), and REF (right) simulations at the surface averaged between 15 and 31 March 2006.

45 1231Figure 2 :
Figure 2: Time series of modeled and observed surface concentrations of major 1232

Figure 4 :Fig. 4 .
Figure 4: Average diurnal profiles of oxygenated organic aerosol (OOA or SOA, µg m -3 ) as observed and modeled (ROB run) during the 15-30 March comparison period (left panels) and the high biomass burning period (18-22 March, right panels) at (a) T0 and (b) T1.Black dots represent PMF values derived from AMS observations.Color polygons represent the contribution of each modeled SOA compound including SOA from biogenics (Bio V-SOA), traditional anthropogenics (Ant V-SOA), and non-traditional anthropogenics (Ant SI-SOA) and biomass burning (BB SI-SOA).

Figure 5 :Fig. 5 .
Figure 5: Same as Figure 4 but for GRI simulation.12561257 Figure 6a and b shows the spatial distribution of O/C ratio over the Mexico City region as predicted by the two SOA treatments.The lowest values are found in Mexico City and other smaller urban areas which have a stronger influence of primary unprocessed organic material.The O/C increases progressively in the city outflow region as the POA gets diluted and as more SOA is formed with a higher contribution www.atmos-chem-phys.net/10/5491/2010/Atmos.Chem.Phys., 10, 5491-5514, 2010 b) 1255

Figure 6 :
Figure 6: Average mass O/C ratios over the Mexico City region as predicted by the ROB (a) 1261

Figure 6 :Fig. 6 .
Figure 6: Average mass O/C ratios over the Mexico City region as predicted by the ROB (a) and GRI (b) 1261 1265

Figure 7 :Fig. 7 .
Figure 7: Observed and predicted concentrations of CO, OOA, HOA, BBOA and total organic 1266 matter (µg m -3 ) along the G1 flight path during March 19, 2006.The blue line refers to the REF 1267

Figure 9 :Fig. 9 .
Figure 9: Same as Figure 7 except for the C130 flight of March 29.Scatter plot comparing observed and 1274 1280

Figure 10 :
Figure 10: Observed and predicted OA concentrations during eight MILAGRO G-1 flights in the vicinity 1281

Figure 11 :Fig. 10 .
Figure 11: OA (a) and OOA (b) concentration normalized to CO mixing ratios as a function of 1287

Figure 11 :
Figure 11: OA (a) and OOA (b) concentration normalized to CO mixing ratios as a function of

Figure 12 :Fig. 12 .
Figure 12: Comparison of observed to predicted SOA ratios as determined in the present study 1293

Table 1 .
The terminology used for the various fractions and sources of organic compounds.
Robinson et al. (2007)phase, and i is the number of components in the emissions (the number of lumped species).The Clausius-Clapeyron equation is used to account for changes in C i * with temperature with the enthalpies of vaporiza-tion estimated for each surrogate species byRobinson et al. (2007).www.atmos-chem-phys.net/10/5491/2010/Atmos.Chem.Phys., 10, 5491-5514, 2010

Table 2 .
Parameters used to compute the partitioning of semi-volatile primary organics.

Table 3 .
Comparison of observed and simulated OOA concentrations during the MILAGRO experiment for the (REF) reference model simulation with the SOA formation from traditional anthropogenic and biogenic VOC precursors (from . To measure the linear dependence between modeled and observed values, Pearson correlation coefficients are used throughout the paper.− O i ) 2 ; where N is the number of samples, O i are observations and M i are model predictions.These indicators were computed on all OOA available data from the MILAGRO ground sites during March 2006.
i (M i

Table 4 .
Kleinman et al. (2008)be satisfied for the comparisons of OA/CO ratios.This choice follows the recommendations ofKleinman et al. (2008).