Impacts of water partitioning and polarity of organic compounds on secondary organic aerosol over eastern China

Secondary organic aerosol (SOA) is an important component of fine particular matter (PM2.5). Most air quality models use an equilibrium partitioning method along with the saturation vapor pressure (SVP) of semivolatile organic compounds (SVOCs) to predict SOA formation. However, the models typically assume that the organic particulate matter (OPM) is an ideal mixture and ignore the partitioning of water vapor to OPM. In this study, the Community Multiscale Air Quality model (CMAQ) is updated to investigate the impacts of water vapor partitioning and nonideality of the organic–water mixture on SOA formation during winter (January) and summer (July) of 2013 over eastern China. The updated model treats the partitioning of water vapor molecules into OPM and uses the universal functional activity coefficient (UNIFAC) model to estimate the activity coefficients of species in the organic–water mixture. The modified model can generally capture the observed surface organic carbon (OC) with a correlation coefficient R of 0.7 and the surface organic aerosol (OA) with the mean fractional bias (MFB) and mean fractional error (MFE) of−0.28 and 0.54, respectively. SOA concentration shows significant seasonal and spatial variations, with high concentrations in the North China Plain (NCP), central China, and the Sichuan Basin (SCB) regions during winter (up to 25 μg m−3) and in the Yangtze River Delta (YRD) during summer (up to 16 μg m−3). In winter, SOA decreases slightly in the updated model, with a monthly averaged relative change of 10 %– 20 % in the highly concentrated areas, mainly due to organic– water interactions. The monthly averaged concentration of SOA increases greatly in summer, by 20 %–50 % at the surface and 30 %–60 % in the whole column. The increase in SOA is mainly due to the increase in biogenic SOA in inland areas and anthropogenic SOA in coastal areas. As a result, the averaged aerosol optical depth (AOD) is increased Published by Copernicus Publications on behalf of the European Geosciences Union. 7292 J. Li et al.: Impacts of water partitioning and polarity of organic compounds on secondary organic aerosol by up to 10 %, and the cooling effect of aerosol radiative forcing (ARF) is enhanced by up to 15 % over the YRD in summer. The aerosol liquid water content associated with OPM (ALWorg) at the surface is relatively high in inland areas in winter and over the ocean in summer, with a monthly averaged concentration of 0.5–3.0 and 5–7 μg m−3, respectively. The hygroscopicity parameter κ of OA based on the κ–Köhler theory is determined using the modeled ALWorg. The correlation of κ with the O : C ratio varies significantly across different cities and seasons. Analysis of two representative cities, Jinan (in the NCP) and Nanjing (in the YRD), shows that the impacts of water partitioning and nonideality of the organic–water mixture on SOA are sensitive to temperature, relative humidity (RH), and the SVP of SVOCs. The two processes exhibit opposite impacts on SOA in eastern China. Water uptake increases SOA by up to 80 % in the organic phase, while including nonunity activity coefficients decreases SOA by up to 50 %. Our results indicate that both water partitioning into OPM and the activity coefficients of the condensed organics should be considered in simulating SOA formation from gas–particle partitioning, especially in hot and humid environments.

by up to 10 %, and the cooling effect of aerosol radiative forcing (ARF) is enhanced by up to 15 % over the YRD in summer. The aerosol liquid water content associated with OPM (ALW org ) at the surface is relatively high in inland areas in winter and over the ocean in summer, with a monthly averaged concentration of 0.5-3.0 and 5-7 µg m −3 , respectively. The hygroscopicity parameter κ of OA based on the κ-Köhler theory is determined using the modeled ALW org . The correlation of κ with the O : C ratio varies significantly across different cities and seasons. Analysis of two representative cities, Jinan (in the NCP) and Nanjing (in the YRD), shows that the impacts of water partitioning and nonideality of the organic-water mixture on SOA are sensitive to temperature, relative humidity (RH), and the SVP of SVOCs. The two processes exhibit opposite impacts on SOA in eastern China. Water uptake increases SOA by up to 80 % in the organic phase, while including nonunity activity coefficients decreases SOA by up to 50 %. Our results indicate that both water partitioning into OPM and the activity coefficients of the condensed organics should be considered in simulating SOA formation from gas-particle partitioning, especially in hot and humid environments.

Introduction
Secondary organic aerosol (SOA) is formed via a complex interaction of volatile organic compounds (VOCs) with oxidants and primary particles emitted from anthropogenic and biogenic sources in the atmosphere. As an important component of fine particular matter (PM 2.5 ), SOA can cause severe air pollution in urban and suburban areas (Huang et al., 2014) and exhibits adverse health effects (Atkinson et al., 2014). SOA also plays an important role in new particle formation and particle growth that further contribute to the enhancement of cloud condensation nuclei (CCNs; Wiedensohler et al., 2009;Ehn et al., 2014). This will, in turn, impact the atmospheric aerosol burden, precipitation and water circulation, solar radiation budget, and climate (Ramanathan et al., 2001). However, the extents of those influences are not well understood so far due to the high uncertainties associated with the formation as well as the physical and chemical properties of SOA (Shrivastava et al., 2017). Large gaps still exist in SOA mass loading and properties between models and observations (Gentner et al., 2017;Ervens et al., 2011;Hayes et al., 2015). Therefore, it is crucial to explore and resolve this issue to improve our knowledge of the roles of SOA in the environment, human health, and climate.
Gas-particle partitioning of semivolatile and low-volatile organic compounds (SVOCs and LVOCs) generated from VOC oxidation is an important pathway of SOA formation. In most current chemical transport models (CTMs), this process is treated as an equilibrium partitioning process that depends on the mass concentration of the organic particulate matter (OPM), ambient temperature (T ), the mean molecular weight of the OPM, and the volatility of condensed organics (Pankow, 1994). The formation of condensed organic products is commonly represented by lumped surrogate SVOCs in a two-product model with volatilities and SVOC yields fitted to chamber experiments (Odum et al., 1996). To better represent the volatility of primary organic aerosol (POA) and the multigeneration oxidation of SVOCs to a wider range, Donahue et al. (2006) proposed the volatility basis set (VBS) model, in which the mass yields of SVOCs are fitted to a fixed number of volatility bins (usually 0.01-10 5 µg m −3 ). The VBS model has been adopted by several CTMs (WRF-Chem, GEOS-Chem, etc.).
Although the models can capture the general trend of SOA evolution and mass concentration to some extent (J. Li et al., 2017;Bergström et al., 2012;Woody et al., 2016), two key factors currently neglected in models may lead to biases: (1) the molecular structures and interactions of functional groups (−OH, −C = O, −COOH, etc.) of condensed organics (nonideality) and (2) partitioning of water vapor, the most abundant atmospheric constituent besides O 2 and N 2 , to OPM. Nonideality alters the volatility of condensed organics and thus their contributions to the total SOA mass loading (Cappa et al., 2008). Water partitioning into OPM can reduce the partial pressure of organics due to the effect of Raoult's law (Prisle et al., 2010) and lead to increases in SOA mass. The amount of aerosol liquid water associated with organics (ALW org ) may vary for different precursors (Healy et al., 2009;Prisle et al., 2010). The above two aspects will not only affect the chemical composition of SOA but also the inorganic portion (Ansari and Pandis, 2000) and optical properties (Denjean et al., 2015) of aerosols.
Laboratory and field studies have observed water absorbed by SOA from a variety of precursor VOCs (Lambe et al., 2011;D. F. Zhao et al., 2016;Asa-Awuku et al., 2010;Varutbangkul et al., 2006). The hygroscopicity of SOA, quantitatively described by the hygroscopicity parameter κ is correlated with the oxygen-to-carbon ratio (O : C) and increases with more oxidized SOA during photochemical aging (Lambe et al., 2011;D. F. Zhao et al., 2016). The OPMassociated water can be estimated using the κ-Köhler theory under the Zdanovskii-Stokes-Robinson (ZSR) assumption of no interactions between any constituents in aerosols (Petters and Kreidenweis, 2007). The total water content is the summation of water associated with each solute at the same water activity. Guo et al. (2015) found that this simplified method, along with the ISORROPIA model, which is used to predict aerosol liquid water (ALW) associated with the inorganic portion of aerosols, could reproduce the observed total ALW in the ambient environment. Pye et al. (2017) found that the modeled organic aerosol (OA) improved significantly but biased high at nighttime when ALW org is included in the calculation. However, as the interaction among organic species and between organics and water in the organic-water mixture has been shown to play an important role in SOA formation and water partitioning to OPM (Kim et al., 2019), the ALW org estimated by the κ-Köhler theory and its impact on SOA might not be accurate. Using the universal functional activity coefficient (UNIFAC) method (Fredenslund et al., 1975) for calculating activity coefficients of the organicwater mixture, it was found that in the eastern US, where biogenic SOA dominated the OA, considering ALW org leads to a significant increase in predicted SOA Jathar et al., 2016).
China has been suffering from severe PM 2.5 pollution especially in the eastern region, with fast urbanization and economic development (Fu and Chen, 2017). SOA is a very significant component of PM 2.5 in China, contributing about 20 %-50 % (Y. J. . The fraction of SOA in OA increases during haze events (Huang et al., 2014;. Previous modeling studies in China indicate that SOA was underpredicted (Lin et al., 2016;Jiang et al., 2012), and the impacts of nonideality and water-OPM partitioning on modeled SOA have not been evaluated.
In this study, regional simulations of SOA during January and July of 2013 over eastern China were conducted to investigate the seasonal variation in SOA due to water partitioning into OPM. The model performance was evaluated against observed meteorological parameters (temperature and relative humidity, RH) as well as PM 2.5 , organic carbon (OC), and OA at ground monitoring sites. The regional and seasonal impacts on SOA, ALW org , and properties of aerosols were quantified. Lastly, influences of the results by several factors including meteorological parameters, estimations of saturation vapor pressures (SVP) of condensed organics, and the individual impacts of ALW org and nonideality of the organicwater mixture on SOA prediction were analyzed.

Model description
The Community Multiscale Air Quality model (CMAQ v5.0.1), coupled with a modified SAPRC-11, was used in this study. Model configurations were largely based on that used by Hu et al. (2016) as summarized below. Firstly, SAPRC-11 was expanded for a more detailed treatment of isoprene oxidation and for tracking dicarbonyl products (glyoxal and methylglyoxal) from different groups of major precursors . Secondly, SOA from isoprene epoxydiols (IEPOX), methacrylic acid epoxide (MAE), and dicarbonyls through surface-controlled irreversible reactive uptake was added Li et al., 2015;Liu et al., 2020;Ying et al., 2015). Thirdly, the heterogeneous formation of secondary nitrate and sulfate from NO 2 and SO 2 reactions on the particle surfaces  was added, which is an important source of secondary inorganic aerosol (Zheng et al., 2015) and improves model estimates of nitrate and sulfate (Qiao et al., 2018;Shi et al., 2017). Fourthly, SOA yields were corrected for vapor wall loss . Impacts of the above updates on model performances have been extensively discussed in the cited work and are not further investigated in the current study.
The SOA module mostly follows Pankow et al. (2015). Two types of SOA as traditionally treated in CMAQ were considered: the "semivolatile" (SV) portion that formed via equilibrium absorption partitioning of SVOCs and the "nonvolatile" (NV) portion that includes the oligomers and SOA formed via direct oxidation of aromatics at low NO x . SOA from dicarbonyls, IEPOX, and MAE was formed by irreversible reactive uptake and categorized as NV SOA in the current model as well. Some studies investigated SOA from glyoxal, methylglyoxal, and IEPOX using detailed reactions and reversible pathways in models or observed as reversible processes in chamber experiments, leading to a relatively lower SOA yield compared to the surface-controlled irreversible uptake (Lim et al., 2013;Knote et al., 2014;Galloway et al., 2009;El-Sayed et al., 2018;. The nonvolatile assumption used in this paper allows an upper-limit estimation of the importance of these additional SOA formation pathways. POA was treated as nonvolatile and nonreactive. The mass distribution of SVOCs between the gas phase and particle phase follows the equation where K p,i (m 3 µg −1 ) is the gas-particle partitioning constant for compound i, F i (µg m −3 ) is the concentration of species i in the particle phase, A i (µg m −3 ) is the concentration of species i in the gas phase, and M (µg m −3 ) is the total mass concentration of the absorbing organic phase (i.e., OPM). The gas-particle partitioning constant K p,i is dependent on the chemical composition of the OPM. Pankow (1994) derived K p,i for SVOCs partitioning into an absorbing organic phase as where p o L,i (atm) is the SVP of the pure compound i at temperature T (K), ξ i is the activity coefficient of species i in the absorbing organic phase, MW (g mol −1 ) is the average molecular weight of the OPM, R (8.314 J mol −1 K −1 ) is the gas constant, and 10 6 is used to convert the units to cubic meters per microgram.
There are 12 lumped SVOCs generated by oxidation of alkanes, alkenes, and aromatics under different NO x conditions and 8 NV organic products as listed in Tables S1 and S2 in the Supplement. More details about the lumped precursors such as formation conditions ("high" and "low" NO x ), lumping species and method, and yields from parent VOCs can be found in Carlton et al. (2010) and are summarized in the Supplement. Activity coefficients of SVOCs were calculated based on the composition of the OPM using the UNI-FAC method, with carbon number (n c ), functional groups, and energy interaction parameters assigned to both SV and NV compounds . The UNIFAC model is one of the most commonly used models with which activity coefficients of condensed organics and their interactions with water can be estimated. This method has been adopted to investigate the impacts of nonideality and water partitioning into OPM on SOA for different precursors in box models (Seinfeld et al., 2001;Bowman and Melton, 2004) and CTMs Kim et al., 2019). In the current model, POA was assumed to have a bulk composition of 10 categories of surrogate species (Table S3), as used by . POA is also involved in the calculation of activity coefficients for the organic-water mixture. Detailed information about the surrogate species, including their structures and properties, can be found in Li et al. (2015) and references therein.
In addition to organic compounds, water partitioning into OPM is enabled according to Eqs. (1) and (2). In such a case, the absorbing phase in Eq. (1) includes both organic aerosol and water associated with OPM. As water condenses in the absorbing organic phase, it will further alter the molar fraction of each composition, the activity coefficient of SVOCs, and the SV SOA mass concentrations as a result. In the current model, we assumed no interactions between the inorganic and organic phases.

Estimating the hygroscopicity of OA
Based on the κ-Köhler theory with linearly additive hygroscopic behavior of each component of the mixed particle, ALW org is related to the hygroscopicity parameter for the organic mixture (κ org ) by Eq. (3) (Petters and Kreidenweis, 2007): where ρ w is the density of water (assumed to be 1 g cm −1 ), V org is the volume concentration of organics, and a w is the water activity (assumed to be the same as RH). Since ALW org in this study is calculated mechanistically using the partitioning theory, κ org can be estimated by rearranging Eq. (3): V org can be estimated from the modeled mass concentration of OA, assuming the density of OA to be 1.2 g cm −3 . In many studies, κ org is assumed to increase linearly with the oxidation state of OA, expressed as the O : C ratio (Massoli et al., 2010;Duplissy et al., 2011;Lambe et al., 2011). The correlation of κ org and the O : C ratio in nine representative cities was evaluated during January and July of 2013 with the reduced major axis regression method (Ayers, 2001). The O : C ratio was calculated using Eq. (5): where f i and (O : C) i are the molar fraction and O : C ratio of organic aerosol component i, respectively. For POA, a fixed molar fraction and composition were assumed following Li et al. (2015). For SOA, the O : C ratio was calculated by using the ratio of organic matter to organic carbon (OM : OC) following Simon and Bhave (2012): The OM : OC ratio of each SOA component follows Pankow et al. (2015) as shown in Tables S1 and S2.

Model application
The simulation domain has a horizontal resolution of 36 km × 36 km (100 × 100 grids) and a vertical structure of 18 layers up to 21 km, which covers eastern China as shown in Fig Table S4 and Fig. S2. Meteorological fields were generated using the Weather Research and Forecasting (WRF) model v3.6.1 with initial and boundary conditions from the National Centers for Environmental Prediction Final (NCEP FNL) Operational Model Global Tropospheric Analyses dataset. More details about the model application can be found in Hu et al. (2016).
Four scenarios are investigated in this study. The base case (BC) applies the default secondary organic aerosol module of CMAQ v5.0.1. In this case, no water partitioning into OPM is considered. Lumped semivolatile products from the oxidation of various precursors partition into a single organic phase, which is considered as an ideal mixture of POA and SOA with γ org = 1. The water case (C1) includes water partitioning into OPM, which is again considered as an ideal solution (γ org =1 and γ H 2 O = 1). The UNIFAC case (C2) considers the interaction between organic constituents with UNIFAC-calculated activity coefficients (γ org = 1) but does not allow water partitioning into OPM. The combined case (C3) allows both water partitioning and interactions between all constituents (including water and organics) using UNIFAC-calculated activity coefficients (γ org = 1 and γ H 2 O = 1). The results of the BC and C3 are used to examine the overall impacts of water partitioning into OPM and polarity of organics on SOA and ALW org , as shown in Sect. 3.1-3.4. The separate influences of those two processes on SOA from C1 and C2 are discussed in Sect. 3.5.
3 Results and discussion

Model evaluation
The meteorological inputs and emissions have been used in several previous publications. Model performance on meteorological parameters (temperature and RH), gaseous species, and gas and aerosol concentrations have been extensively evaluated (Hu et al., 2016Qiao et al., 2018;Shi et al., 2017). A summary of the model performance related to this study is provided below. Observed meteorological data were obtained from the National Climatic Data Center (ftp://ftp.ncdc.noaa.gov/pub/data/noaa, last access: 15 February 2019). Observations of OC at two urban locations, Beijing Wang et al., 2015) and Guangzhou (Lai et al., 2016), and of OA in Beijing (Sun et al., 2013)  Temperature and RH are the two meteorological factors that affect SOA formation. Table 1 Table S5. Overall, WRF tends to underestimate both temperature and RH. The model shows better agreement with observed temperature as R is higher than that of RH. Both temperature and RH are well captured by the model in the YRD, the Pearl River Delta (PRD), and the central regions of China (the major regions of eastern China). In these re- gions, MB and GE of temperature are −1.2-0.7 and 1.8-2.6 K, respectively, which are −11.8 %-5.6 % and 9.2 %-16.8 % for RH, respectively.
Model estimates of daily organic carbon (OC) from case C3 are compared with measurements at monitoring sites in Beijing and Guangzhou in January of 2013 (Fig. 1a). The factors used to convert SOA to OC (SOC) are listed in Tables S1 and S2. OC from POA (POC) is directly predicted by the model. Generally, the ratio between modeled and observed OC concentration falls in the range of 1 : 2 to 2 : 1, with an R value of 0.7. The model tends to underestimate OC on high-concentration days. Overall, the mean fractional bias (MFB) and mean fractional error (MFE) of OC are −0.20 and 0.27, which are within the criteria (MFB ≤ ± 0.6; MFE ≤ 0.75) suggested by the EPA (2007). The bias in OC might be due to underestimated POA emissions and underpredicted SOA in CMAQ from missing precursors . No significant differences in OC are observed in C3 compared to BC (not shown), likely due to the low-biased SOA predicted in the current model that limits the impact of ALW org on SOA formation.
The underestimation of SOA can be seen from Fig. 1b as well. CMAQ can well capture the observed diurnal variation in OA in Beijing during wintertime except for the underestimates of peak values. The correlation coefficient of modeled to observed OA is 0.55. We find a 25 % underestimate of OA on average. Better agreement between the model and observations is shown on nonpolluted days (daily averaged concentration less than 75 µg m −3 ). The MFB and MFE of Table 1. Statistical analysis of modeled temperature (K) and relative humidity (%) of January and July of 2013 at the monitoring sites in different geographical areas as shown in Fig. S1. The cities are listed in Table S5. polluted days are −0.38 and 0.64, respectively, which are worse than those of nonpolluted days (−0.26 for MFB and 0.52 for MFE). The overall MFB and MFE of OA during January are −0.28 and 0.54, respectively, which are within the criteria (MFB ≤ ± 0.6; MFE ≤ 0.75) suggested by the EPA (2007). Again, no apparent changes in SOA or OA are observed between case C3 and BC (not shown) since POA is predicted to be the primary contributor to OA in Beijing in winter in the current model, with an averaged SOA/POA ratio of 0.12. This ratio is much lower than the field observation of about 0.45-1.94 (Zhao et al., 2019;Sun et al., 2013Sun et al., , 2016. The bias might be due to the missing SOA converted by partitioning and aging of semivolatile POA as well as oxidation from intermediate volatile organic compounds (IVOCs) and VOC oxidation products. Those pathways are shown to be important for SOA formation by modeling, field, and chamber studies (Hodzic et al., 2010;Jimenez et al., 2009;Murphy et al., 2017;Robinson et al., 2007;Shrivastava et al., 2008;Tkacik et al., 2012;Zhao et al., 2014;.
A sensitivity test was performed by using the newest CMAQ model version 5.3.1, which includes all the above processes in the aerosol module. The SOA/POA ratio in Beijing is improved greatly in winter, at 0.83. However, high uncertainties still exist in the emissions of the involved precursors and characterization of SOA formation through these processes, requiring further constraints by observations. Their influences on SOA due to water partitioning into OPM and nonideality of the organic-water mixture will be evaluated in a future study.
Due to the lack of observed OC and OA in July of 2013, model performances are evaluated by comparing predicted and observed PM 2.5 at ground sites (Fig. S1) as an alternative, as shown in Fig. S3. Generally, the model can well reproduce the diurnal variation in PM 2.5 in most regions. Predicted PM 2.5 on highly concentrated days is biased low, especially in the North China Plain (NCP). The NCP has the highest PM 2.5 , ranging from 60 to 300 µg m −3 . The bias in modeled PM 2.5 is significant in cities in the northwest. This might be due to missing dust emissions in the current inventory (Hu et al., 2016). To further evaluate the model performance, the averaged MFB and MFE of modeled PM 2.5 are plotted against observations of each site as shown in Fig. S4. The criteria of MFB and MFE followed recommendations by Boylan and Russell (2006). Our model performs well since most of the predictions meet the criteria, and a large fraction (> 58 %) meets the goal. The averaged MFB and MFE across all the sites are −0.28 and 0.39, respectively, indicating slight underestimation of PM 2.5 by the model.

Impacts on SOA and ALW org
The spatial distribution of SOA varies greatly in the two seasons. In winter, SOA are relatively high in the eastern SCB and the central and eastern provinces of Shandong, Henan, Anhui, and Hubei (Figs. 2 and S5). Monthly averaged Figure 2. Monthly averaged SOA from C3 and changes in SOA due to water partitioning into OPM and nonideality of the organic-water mixture. "Abs. Diff." represents absolute differences (C3-BC); "Rel. Diff." represents relative differences ((C3-BC)/BC; %). SOA concentrations in the above areas are up to 25 and 15-20 µg m −3 , respectively. Anthropogenic emissions such as dicarbonyl products from the oxidation of xylene and toluene  are the major sources of SOA (Fig. S6). In summer, surface SOA is high in the northeast, the NCP, and the YRD. The highest SOA occurs in Shanghai and Jiangsu provinces as well as the coastal area of the Yellow Sea, with values of ∼ 9-16 µg m −3 at the surface and ∼ 20-25 mg m −2 in the column (col-SOA) of the atmosphere below 21 km (Fig. S5). In contrast to winter SOA, a significant fraction of summer SOA originates from biogenic emissions (Fig. S7). Anthropogenic SOA is highly concentrated in the coastal areas of the Yellow Sea and Bohai Bay.
The combined effects of water partitioning into OPM and nonideality on SOA formation also exhibit strong seasonal variation. In winter, SOA is slightly decreased by 1.5 µg m −3 (10 %-20 %) on average at the surface (Fig. 2) and less than ∼ 1 mg m −2 (20 %) in the column (Fig. S5) over high-SOA regions where anthropogenic sources dominate. We show later that the decrease in SOA is mainly due to the large activity coefficients of SVOCs, which decrease K p,i . In summer, higher temperature and RH promote water partitioning and SOA formation so that SOA apparently increases over the entire domain, with the highest enhancement of 2-4 µg m −3 (20 %-50 %) at the surface (Fig. 2) and 4-6 mg m −2 (30 %-60 %) in the column (Fig. S5) over the YRD and the coastal area of the Yellow Sea. Anthropogenic SOA dominates the total change in winter as shown in Fig. S6. In summer, the  The relationship between κ org and O : C is fitted by a linear function with reduced major-axis regression (blue lines) and an exponential function (red lines), respectively. κ 01 and κ 07 represent the fitted correlation for January and July, respectively. increase in SOA is attributed to biogenic sources in inland areas and anthropogenic sources over the ocean (Fig. S7).
Regional distribution of ALW org is similar to the change in SOA as shown in Fig. 3. In winter, a maximum averaged concentration of 3.0 µg m −3 for ALW org occurs in the high-SOA region, where significant changes in SOA also occurs. In other areas, the averaged concentration of ALW org is about 0.5-1.5 µg m −3 . Overall, the average ratio of ALW org to SOA is about 0.1-0.3 in winter. In summer, water partitioning into OPM mostly occurs in the eastern coastal area at the surface, where a significant increase in anthropogenic SOA (such as those from toluene and xylenes) is observed. This might be due to the high polarity of anthropogenic SVOCs (having more −COOH groups) that absorb more water. In the coastal area, the averaged concentration of ALW org is about 5-7 µg m −3 , with an ALW org /SOA ratio of 0.5-1.0. Over land, the averaged concentration of ALW org is about 1-3 µg m −3 (ALW org /SOA ratio of 0.2-0.5) in northeastern and eastern China. Water partitioning is mostly associated with biogenic SOA originated from isoprene and monoterpene oxidation that produces SVOCs with abundant OH groups.
Based on the column concentrations of ALW org and the ALW org /SOA ratio (Fig. S8), more ALW org must have occurred in winter in the southern and southwestern regions at higher levels, where a significant increase in col-SOA occurs (Fig. S5). The averaged col-ALW org /col-SOA ratio in the high-SOA area is 0.1-0.3. In summer, the ALW org must be high at higher altitudes over the central regions in China. The maximum col-ALW org is about 7 mg m −2 over the YRD, with a col-ALW org /col-SOA ratio of about 0.3.

Impacts on aerosol properties
Since ALW org is determined in C3, the values of κ org can be estimated from the modeled ALW org , OA, and RH using Eq. (4). Nine representative cities were selected to investigate the relationship of κ org and O : C and its seasonal variation as shown in Figs. S9 and S10. The results of all the cities in winter and summer are merged and analyzed as shown in Fig. 4. Pairs of κ org and O : C data are grouped into 10 O : C bins, and the averaged κ org in each bin is then calculated. Overall, the estimated O : C ratio is within the range of 0.2-0.8. The averaged κ org in each O : C bin is less than 0.1 in winter, with the highest value in Guangzhou. As more ALW org is formed in summer, the averaged κ org also increases greatly, with the highest value of 0.35 in Beijing. The linear correlation between κ org and O : C shows significant spatial and seasonal variations. For example, the slope of κ org -O : C is 70 %-90 % smaller in winter than in summer in northern cities such as Shenyang, Beijing, Zhengzhou, and Xi'an. However, in Guangzhou, the slope is 83% higher in winter than in summer. In Chengdu, the slope is quite similar in both seasons. Overall, the slope of κ org vs. O : C in the nine cities is 0.16 in winter and 0.40 in summer. Most of the fitted linear correlations of the individual cities fall outside of the range of 0.18-0.37 suggested in previous studies (Duplissy et al., 2011;Lambe et al., 2011;Massoli et al., 2010;Chang et al., 2010), indicating that the hygroscopicity of organic aerosol cannot be simply represented by a single parameter such as the O : C ratio (Rickards et al., 2013). In both seasons, κ org approaches zero and negative values as O : C decreases, which might be due to the linear regression of κ org and O : C. To avoid this, an exponential fitting of the two variables is performed so that κ org falls in the range of (0, 1) and is positively correlated with O : C. In this case, The impacts on aerosol optical depth (AOD) and aerosol radiative forcing (ARF) are further investigated. Figure 5 shows the monthly averaged AOD at 550 nm in January and July of 2013. It is calculated by summarizing the product of the model-estimated extinction coefficient of fine particles (b ext,i ) multiplied by the thickness (HL i ) in each layer: where N is the number of layers. There are two methods to estimate the aerosol extinction coefficient in CMAQ v5.0.1. One is based on the Mie theory and the predicted aerosol component concentrations (b ext,m ), and the other is based on correlation with the IMPROVE monitoring network, which considers the impacts of hygroscopicity of different aerosol components (b ext,r ; Malm et al., 1994). AOD calculated with the two types of extinction coefficient are denoted as AOD m and AOD r , respectively. In Fig. 5, a clear pattern of high AOD r in the SCB and the NCP and low AOD r in western China is observed in both winter and summer, which is consistent with previous studies (He et al., , 2019Luo et al., 2014). An identical trend in AOD m is shown in Fig. S11. The monthly averaged AOD r ranges from 1.0 to 3.2 in January and from 0.3 to 0.9 in July. AOD m is lower than AOD r , falling within the range of 0.7-2.2 in January and 0.3-0.6 in July. The model significantly overestimates AOD in January but agrees bet- Figure 6. Monthly averaged shortwave direct aerosol radiative forcing at the top of the atmosphere from C3 and the relative changes due to water partitioning into OPM and nonideality of the organicwater mixture during January and July of 2013. "Rel. Diff." represents relative differences ((C3-BC)/BC; %). ter with observations from MODIS, where AOD is high in July (Fig. S12). The bias in the predicted AOD might be partially due to the empirical equation applied in the calculation of AOD in CMAQ Liu et al., 2010) and partially due to the uncertainties of fine AOD over land from the MODIS data Levy et al., 2010). The increase in AOD due to ALW org shows a strong spatial and seasonal pattern. In winter, there are no significant changes in AOD r across the whole domain due to insignificant changes in SOA. In summer, AOD r increases significantly in the YRD and the adjacent areas by up to 10 %.
ARF represents the change in the radiative flux at the top of the atmosphere due to aerosols. An offline version of the Shortwave Radiative Transfer Model for GCMs (RRTMG_SW) was used to calculate the direct radiative effect of aerosols on shortwave radiation (Iacono et al., 2008). Generally, fine aerosols exhibit cooling effects on the shortwave radiation in both winter and summer over the entire domain as shown in Fig. 6. This impact is much stronger in the areas where AOD is high (Fig. 5). ARF is highest in Shandong in winter and in the coastal area near Jiangsu province in summer, amounting to about −5 and −6 W m −2 , respectively. In winter, no significant changes in ARF are observed in eastern China (Fig. 6b). This is likely attributed to an insignificant contribution of SOA to PM 2.5 in winter compared to other components with cooling effects, such as sulfate. In summer, SOA is an important component of PM 2.5 (20 %-60 %), and the effects of water partitioning on shortwave radiation are relatively stronger. An enhancement in the cooling effects of ARF of up to 15 % occurs near the YRD region, where AOD significantly changes as well.

Sensitivity to T , RH, and SVP
Meteorological conditions and SOA precursors affect the impacts of water partitioning on SOA. To better illustrate the dependency of SOA on temperature, RH, and SVP of SVOCs, an offline calculation of SOA formation was performed in two representative cities (Jinan, in the NCP, during winter and Nanjing, in the YRD, during summer) when the daily maximum SOA increases occurred. We assume temperature (T ) and water vapor mixing ratio (QV) to be within the range of X±σ , where X and σ are the mean and standard deviation, respectively, calculated based on WRF predictions at each location. We choose 10 evenly distributed values for T and QV within the range of X ± σ . The temperature dependence parameter of SVP ( H /R) is also scaled separately by 0.2, 0.8, 1.4, and 2.0 for all the SVOCs. As shown in Fig. 7, SOA exhibits a negative correlation with T and a positive correlation with QV in both cities. SOA is more sensitive to QV under cold conditions (Jinan) and to temperature under hot conditions (Nanjing). When the temperature is fixed, Figure 9. Monthly averaged impacts of nonideality of the organicwater mixture on SOA. "Abs. Diff." represents absolute differences (C3-C1); "Rel. Diff." represents relative differences ((C3-C1)/C1; %). the sensitivity of SOA to H /R is different in the two cities. We find more changes in SOA across H /R in Jinan. This is attributed to the temperature correction factor (ζ corr ) of K p in CMAQ as defined below: where T ref is the reference temperature (298 K), and T is temperature. According to Fig. 7, the range of T is 265-274 K in Jinan and 300-307 K in Nanjing. The deviation of temperature from the reference value (298 K) is greater in Jinan than in Nanjing. Therefore, the unit change in H /R causes greater variations of ζ corr and thus K p in Jinan. As a result, SOA is more sensitive to H /R. The impacts of SVP estimation on SOA are more significant in winter than in summer.

Separate impacts of water partitioning and polarity of condensed organics
The impacts of water partitioning into OPM and nonideality of the organic-water mixture on SOA are in opposite directions. Water partitioning increases SOA by ∼ 10 %-20 % in winter and ∼ 30 %-80 % in summer in most areas of the domain (Fig. 8). This is because the molecular weight of water is quite small and will reduce the mole-averaged molecular weight of the OPM (MW) in Eq. (2)  . The reduced MW further increases K p and promotes mass transfer of SVOCs from the gas phase to OPM. On the other hand, by considering nonideality of the organic-water mixture, activity coefficients of SVOCs are usually greater than 1.0 in this study, leading to a decrease in K p . As a result, the total SOA concentration is reduced by ∼ 10 %-20 % in winter and ∼ 10 %-50 % in summer (Fig. 9). Overall, the final impacts are the combined consequences of the two processes. In winter, the increase in SOA caused by water partitioning is offset by the decrease in SOA due to the polarity of SVOCs in most areas of the domain, resulting in slight decreases in SOA. In summer, the effect of water partitioning overcomes that of SVOC polarity, so the total SOA loading increases. This further leads to an enhanced attenuation of shortwave solar radiation and cooling of the atmosphere.

Conclusions
The WRF-CMAQ model was used to investigate the impacts of water partitioning into OPM and nonideality of the organic-water mixture on SOA formation over eastern China during January and July of 2013. SOA is greatly enhanced in summer, especially in the YRD and over the Yellow Sea, by up to 50 % and 60 % at the surface and in the whole column, respectively. No significant impacts on SOA are observed in winter. This might be due to the underestimation of SOA in the current model. ALW org is highly correlated with the change in SOA, with a ratio of ALW org to SOA of 0.1-0.3 and 0.2-1.0 at the surface, where significant changes in SOA occur, in winter and summer, respectively. By using the modeled ALW org , correlations between κ org and O : C are examined in nine representative cities, showing significant spatial and seasonal variations. The increases in SOA lead to an enhancement in the averaged AOD and the cooling effects of aerosols by up to 10 % and 15 %, respectively, in summer. The model-predicted SOA is sensitive to temperature and QV in both seasons, with higher sensitivity to QV during winter and to temperature during summer. Estimation of SVP also affects modeled SOA, especially in a cold environment. The effects of water partitioning into OPM and nonideality of the organic-water mixture on SOA are the opposite.
Since the activity coefficients of SVOCs are mostly greater than 1.0 during the simulated episode, SOA concentration decreases when the nonideality effect is considered. The averaged SOA concentration decreases by up to 20 % in winter and 50 % in summer. Water partitioning alone increases SOA by 10 %-20 % in winter and 30 %-80 % in summer. It should be noticed that the results shown in this study are the lower limit as the current model tends to underestimate SOA. It is crucial to consider both effects in simulating SOA formation under hot and humid conditions in CTMs.
Data availability. The data used in this paper can be provided upon request from the corresponding authors: Qi Ying (qy-ing@civil.tamu.edu) and Jianlin Hu (jianlinhu@nuist.edu.cn).
Author contributions. JL, QY, and JH designed the research and contributed to model development and configurations. YLZ, XW, XL, and YS provided the observations of OC and OA. JL, HZ, QY, and JH analyzed the data. JL prepared the manuscript and all coauthors helped improve the manuscript.
Competing interests. The authors declare that they have no conflict of interest.