Mixing times of organic molecules within secondary organic aerosol particles : a global planetary boundary layer perspective

15 When simulating the formation and life cycle of secondary organic aerosol (SOA) with chemical transport models, it is often assumed that organic molecules are well mixed within SOA particles on the time scale of 1 h. While this assumption has been debated vigorously in the literature, the issue remains unresolved in part due to a lack of information on the mixing times within SOA particles as a function of both temperature and relative humidity. Using laboratory data, meteorological fields, and a chemical transport model, we estimated how often mixing times are < 1 20 h within SOA in the planetary boundary layer (PBL), the region of the atmosphere where SOA concentrations are on average the highest. First, a parameterization for viscosity as a function of temperature and RH was developed for pinene SOA using room-temperature and low-temperature viscosity data for -pinene SOA generated in the laboratory using mass concentrations of 1000 g m. Based on this parameterization, the mixing times within -pinene SOA are < 1 h for 98.5 % and 99.9 % of the occurrences in the PBL during January and July, respectively, when 25 concentrations are significant (total organic aerosol concentrations are > 0.5 μg m at the surface). Next, as a starting point to quantify how often mixing times of organic molecules are < 1 h within -pinene SOA generated using low, atmospherically relevant, mass concentrations, we developed a temperature-independent parameterization for viscosity using the room-temperature viscosity data for α-pinene SOA generated in the laboratory using a mass concentration of 70 μg m. Based on this temperature-independent parameterization, mixing times within -pinene 30 SOA are < 1 h for 45 and 38 % of the occurrences in the PBL during January and July, respectively, when concentrations are significant. However, associated with these conclusions are several caveats, and due to these caveats, we are unable to make strong conclusions about how often mixing times of organic molecules are < 1 h within α-pinene SOA generated using low, atmospherically relevant, mass concentrations. Finally, a parameterization for viscosity of anthropogenic SOA as a function of temperature and RH was developed using sucrose-water data. Based 35 on this parameterization, and assuming sucrose is a good proxy for anthropogenic SOA, 70 % and 83 % of the mixing

Abstract. When simulating the formation and life cycle of secondary organic aerosol (SOA) with chemical transport models, it is often assumed that organic molecules are well mixed within SOA particles on the timescale of 1 h. While this assumption has been debated vigorously in the literature, the issue remains unresolved in part due to a lack of information on the mixing times within SOA particles as a function of both temperature and relative humidity. Using laboratory data, meteorological fields, and a chemical transport model, we estimated how often mixing times are < 1 h within SOA in the planetary boundary layer (PBL), the region of the atmosphere where SOA concentrations are on average the highest. First, a parameterization for viscosity as a function of temperature and RH was developed for α-pinene SOA using room-temperature and low-temperature viscosity data for α-pinene SOA generated in the laboratory using mass concentrations of ∼ 1000 µg m −3 . Based on this parameterization, the mixing times within α-pinene SOA are < 1 h for 98.5 % and 99.9 % of the occurrences in the PBL during January and July, respectively, when concentrations are significant (total organic aerosol concentrations are > 0.5 µg m −3 at the surface). Next, as a starting point to quantify how often mixing times of organic molecules are < 1 h within α-pinene SOA generated using low, atmospherically relevant mass concentrations, we developed a temperatureindependent parameterization for viscosity using the roomtemperature viscosity data for α-pinene SOA generated in the laboratory using a mass concentration of ∼ 70 µg m −3 . Based on this temperature-independent parameterization, mixing times within α-pinene SOA are < 1 h for 27 and 19.5 % of the occurrences in the PBL during January and July, respectively, when concentrations are significant. However, associated with these conclusions are several caveats, and due to these caveats we are unable to make strong conclusions about how often mixing times of organic molecules are < 1 h within α-pinene SOA generated using low, atmospherically relevant mass concentrations. Finally, a parameterization for viscosity of anthropogenic SOA as a function of temperature and RH was developed using sucrose-water data. Based on this parameterization, and assuming sucrose is a good proxy for anthropogenic SOA, 70 and 83 % of the mixing times within anthropogenic SOA in the PBL are < 1 h for January and July, respectively, when concentrations are significant. These percentages are likely lower limits due to the assumptions used to calculate mixing times.
reactions of the SVOCs in the particle phase (Ervens et al., 2011;Hallquist et al., 2009). The term "secondary" indicates the aerosol is formed in the atmosphere rather than emitted directly into the atmosphere in the particle phase. Globally, SOA from biogenic sources dominate, with SOA from anthropogenic sources contributing approximately 10 % to the total SOA budget Spracklen et al., 2011). Major contributors to biogenic SOA are oxidation products of α-pinene and isoprene (Hu et al., 2015;Kanakidou et al., 2005;Pathak et al., 2007) and, as a result, SOA derived from α-pinene and isoprene are the most widely used representatives of biogenic SOA in experimental and modelling studies.
The planetary boundary layer (PBL) is the lowest part of the atmosphere, ranging from the Earth's surface to roughly 1 km in altitude, depending on location and time (Wallace and Hobbs, 2006). Within this region vertical mixing of air masses is rapid, on the order of 30 min (Wallace and Hobbs, 2006). In addition, within the PBL the temperature varies from roughly 265 to 305 K and the relative humidity (RH) varies from roughly 20 to 100 % (see below). SOA concentrations are also on average highest in the PBL (Heald et al., 2011;Wagner et al., 2015).
When simulating the formation, growth, and evaporation of SOA particles with chemical transport models, it is often assumed that SVOCs are well mixed within SOA particles on a timescale of 1 h ). If SVOCs are not well mixed within SOA particles on this timescale, then chemical transport models could incorrectly predict SOA mass concentrations by up to an order of magnitude (Shiraiwa and Seinfeld, 2012) and incorrectly predict the size of SOA particles (Zaveri et al., 2014), with implications for air quality and climate predictions (Seinfeld and Pandis, 2006). Recent research has shown that mixing times of organic molecules within SOA particles can be > 1 h at roomtemperature and at low RHs (Abramson et al., 2013;Grayson et al., 2016;Liu et al., 2016;Perraud et al., 2012;Renbaum-Wolff et al., 2013;Song et al., 2016;Ye et al., 2016;Zhang et al., 2015). In addition, studies have shown that proxies of SOA particles can form glass at low RHs and low temperatures Zobrist et al., 2008). Nevertheless, the conditions that lead to slow mixing times in SOA may be infrequent on a global scale in the PBL. If this is the case, then the assumption of well-mixed SOA particles in chemical transport models should be reasonable. How often mixing times are > 1 h under ambient conditions in the PBL is not well constrained, in part due to the lack of information on mixing times of organic molecules in SOA particles as a function of both RH and temperature.
In the following, we (a) developed a parameterization for the viscosity of α-pinene SOA particles as a function of both RH and temperature, (b) determined the distribution of RH and temperature in the PBL from an archive of meteorological fields, (c) determined the conditions in the PBL when SOA concentrations are significant using a chemical trans-port model, and (d) quantified how often mixing times of SVOCs are > 1 h within α-pinene for ambient temperatures and RHs in the PBL. Mixing times within anthropogenic SOA and the effect of SOA mass concentration on mixing times are also discussed.
2 Materials and methods 2.1 Parameterization for the viscosity of α-pinene SOA as a function of temperature and RH The following data were used to develop a parameterization of the viscosity α-pinene SOA as a function of temperature and RH: (a) room-temperature measurements of viscosity of SOA derived from α-pinene ozonolysis by Grayson et al. (2016)  In these experiments, the SOA was generated with mass concentrations of 121 and 520 µg m −3 . We use the viscosity measurements from Grayson et al. (2016) determined with a mass concentration of 520 µg m −3 to be more consistent with the mass concentrations used by Järvinen et al. (2016). Although there are other room-temperature measurements of the viscosity of α-pinene SOA Hosny et al., 2016;Kidd et al., 2014;Pajunoja et al., 2014;Renbaum-Wolff et al., 2013), we used the roomtemperature measurements from Grayson et al. (2016) because (1) viscosity was measured over a range of relative humidities in this study, (2) the mass concentrations used by Grayson et al. (2016) to generate the SOA were similar to the mass concentrations used by Järvinen et al. (2016), and (3) Grayson et al. (2016) measured the viscosity of the total SOA (both the water-soluble component and the waterinsoluble component). Due to the experimental conditions used by Grayson et al. (2016) and Järvinen et al. (2016), the parameterization developed here is applicable to SOA generated using a mass concentration of ∼ 1000 µg m −3 . We focused on ∼ 1000 µg m −3 because both low-temperature and roomtemperature viscosity measurements have been carried out using this mass concentration. The effect of mass concentration on the viscosity α-pinene SOA is discussed in Sect. 3.5.
To develop a parameterization for viscosity as function of temperature and RH, the following equation was fit to the measurements by Grayson et al. (2016), Järvinen et al. (2016), and Crittenden et al. (2012) (Tables S1-S3): where C 1 and C 2 are constants, k GT is the Gordon-Taylor fitting parameter, T gSOA and T gH 2 O are the glass transition temperatures of dry SOA and water, and w SOA and w H 2 O are the weight fractions of the dry SOA and water in the particles. The weight fractions of the dry SOA and water in the particles were determined from the RH using the following equation : where i is the van 't Hoff factor and n is the number of moles of dry SOA and water in the particles. We assumed a value of 1 for the van 't Hoff factor ) and a dry molecular weight for SOA of 175 g mol −1 (Huff Hartz et al., 2005).
Since the glass transition temperature of water is known (135 K) (Corti et al., 2008), the unknowns in Eq. (1) (and hence fitting parameters) were C 1 , C 2 , k GT , and T gSOA . The values for these parameters retrieved by fitting the equation to the viscosity data discussed above (using a nonlinear curve fitting function in MATLAB) are reported in Table S4.
Equation (1) was based on the Williams-Landel-Ferry (WLF) equation and the Gordon-Taylor equation. The WLF equation provides a relationship between viscosity and temperature: where C 1 and C 2 are constants, T is the temperature, T g is the glass transition temperature, η is the viscosity, and η g is the viscosity at the glass transition (10 12 Pa s). The Gordon-Taylor equation provides a relationship between the glass transition temperature of a mixture and the weight fractions of its components: where w 1 and w 2 are the weight fractions of the solute and water, T g1 and T g2 are the glass transition temperatures of the solute and water, and k GT is a fitting parameter that describes the interaction between the two species. Equations (3) and (4) can be combined to give Eq. (1). Equation (3) (and hence Eq. 1) is valid only at or above the glass transition temperature. As a result, we have not used Eq. (1) to predict viscosities > 10 12 Pa s (which corresponds to mixing times > 5×10 5 h). If the temperature and RH in the PBL were such that the viscosity was greater than 10 12 Pa s, we assigned a viscosity of 10 12 Pa s and a mixing time of 5 × 10 5 hours. This assignment does not affect the conclusions in this paper since a mixing time of 5×10 5 hours is already well above the residence time of SOA particles in the atmosphere.

Organic aerosol concentrations in the planetary boundary layer
To determine the conditions in the PBL when SOA concentrations are significant we used the global chemical transport model GEOS-Chem (http://acmg.seas.harvard.edu/ geos/). The version of GEOS-Chem used (v10-01) includes organic aerosol (OA) formation from SVOC and intermediate-volatility organic compounds (IVOC) , plus new aerosol production from nitrate radical oxidation of isoprene and terpenes and NO xdependent aerosol yields from terpenes . In this version IVOC emissions are spatially distributed based on naphthalene. To estimate SVOC emissions we scaled the default GEOS-Chem primary organic aerosol emissions inventory by 1.27 following . GEOS-Chem was run at a horizontal grid resolution of 4 • latitude by 5 • longitude using GEOS-5 meteorology with 47 vertical layers with a 3-year spin-up period. Shown in Fig. 1 are the monthly averaged total OA concentrations at the surface for the months of January and July 2006. These monthly averaged total OA concentrations were used to remove times and locations where SOA concentrations are not expected to be of major importance for climate, health, or visibility.

RH and temperature in the PBL
Information on the RH and temperature distributions in the global PBL in different seasons was also needed to assess mixing times within SOA particles. First, the PBL heights were determined globally using the 6 h averaged GEOS-5 meteorology fields. Then, temperature and RH in each grid cell within the PBL were determined globally using the 6 h averaged GEOS-5 meteorology fields. To determine if a grid cell was within the PBL, the PBL heights mentioned above were used. The GEOS-5 archive provides temperature and RH at a horizontal grid resolution of 4 • latitude by 5 • longitude and 47 vertical layers.
3 Results and discussion 3.1 Parameterization of viscosity and mixing times within α-pinene SOA particles as a function of RH and temperature Shown in Fig. 2a (contours) is the RH and temperaturedependent parameterization for α-pinene SOA viscosities based on viscosities measured at room temperature (Grayson et al., 2016) and low temperature (Järvinen et al., 2016), as well as the viscosity of water as a function of temperature (Crittenden et al., 2012). From the viscosity parameteriza- tion, the diffusion coefficients of organic molecules within α-pinene SOA particles were calculated using the Stokes-Einstein equation: where D is the diffusion coefficient, k is the Boltzmann constant, T is temperature in Kelvin, η is the dynamic viscosity, and R H is the hydrodynamic radius of the diffusing species. For the calculations, a hydrodynamic radius of 0.38 nm was used for the diffusing organic molecules within SOA based on an assumed molecular weight of 175 g mol −1 (Huff Hartz et al., 2005), a density of 1.3 g cm −3 (Chen and Hopke, 2009;Saathoff et al., 2009), and spherical symmetry. The Stokes-Einstein equation should give reasonable values when the radius of the diffusing molecules is roughly greater than or equal to the radius of the matrix molecules and when the viscosity of the matrix is relatively small (400 Pa s) (Chenyakin et al., 2017;Price et al., 2016). When the viscosity of the matrix is large (10 6 Pa s), the Stokes-Einstein equation can underpredict diffusion coefficients of organic molecules in organic matrices (Champion et al., 1997;Chenyakin et al., 2017;Price et al., 2016). Hence, the diffusion coefficients and mixing times estimated here should be considered lower and upper limits, respectively. From the diffusion coefficients, the mixing times of organic molecules within an α-pinene SOA particle were calculated with the following equation : where τ mix is the mixing time, d is the diameter of an SOA particle, and D is the diffusion coefficient estimated from Eq. (5). For these calculations, it was assumed that the α-pinene SOA particles have a diameter of 200 nm, which is roughly the median diameter in the volume distribution of ambient SOA-containing particles Pöschl et al., 2010;Riipinen et al., 2011). Once the mixing time has elapsed, the concentration of the diffusing molecules at the centre of the particle is within 1/e of the equilibrium concentration . The calculated mixing times (Fig. 2b) illustrate that, as expected, inverse relationships exist between both mixing time and RH, as well as mixing time and temperature. Fig. 3a and b are the normalized frequency counts of temperature and RH in the PBL for the months of January and July 2006, respectively, based on the archive of meteorological fields (GEOS-5) used in the global chemical transport model, GEOS-Chem v10-01. We only included grid points in our analysis when the grid points were within the PBL and the monthly average mass concentration of total OA was > 0.5 µg m −3 at the surface, based on GEOS-Chem v10-01 (Fig. 1). In other words, we included all the grid points in a column up to the top of the PBL when determining frequency counts if the monthly averaged total OA concentration was > 0.5 µg m −3 at the surface. This filtering removes cases where SOA concentrations are not expected to be of major importance for climate, health, or visibility. We chose a mass concentration of > 0.5 µg m −3 for filtering because the mass concentration of total OA at the surface was > 0.5 µg m −3 in all but one of the previous field measurements of OA at remote locations (Spracklen et al., 2011). The normalized frequency counts illustrate that the temperature and RH in the PBL are often in the range of 290-300 K and > 50 % RH for the month of January (Fig. 3a) and in the range of 285-300 K and > 30 % RH for the month of July (Fig. 3b). For reference, shown in Figs. S1 and S2 in the Supplement are the average temperature and RH conditions at the Earth's surface and top of the planetary boundary layer, respectively, for January and July, based on the archive of meteorological fields for 2006 (GEOS-5).

Mixing times of organic molecules within α-pinene
SOA particles the PBL Also shown in Fig. 3a and b are the mixing times within 200 nm α-pinene SOA particles predicted with our parameterization (contours). These results, together with the frequency counts of temperature and RH throughout the vertical column of the PBL, indicate that the mixing times of organic molecules within α-pinene SOA are often < 1 × 10 −1 h for conditions in the PBL. Shown in Fig. 4 are the normalized frequency distributions of mixing times within α-pinene SOA for January and July, based on the data in Fig. 3a and b. Figure 4 suggests that the mixing times within α-pinene SOA are < 1 h for 98.5 and shows the conditions for July. Mixing times (contours) are reported in hours. Frequency counts in the PBL were only included for the conditions where the mass concentration of total organic aerosol was > 0.5 µg m −3 at the surface. The viscosity parameterization used to calculate mixing times was based on α-pinene SOA generated using mass concentrations of ∼ 1000 µg m −3 . 99.9 % of the occurrences in the PBL during January and July, respectively, when monthly average total OA concentrations are > 0.5 µg m −3 at the surface.
Within the PBL, RH increases and temperature decreases with altitude, with both changes being substantial and impacting mixing times in opposite directions. Shown in Fig. 5a-c are calculated monthly average afternoon (13:00-15:00, LT) vertical profiles of temperature, RH, and mixing times within α-pinene SOA over Hyytiälä (boreal forest) and the Amazon (rainforest) for the driest month of the year at each location (the method used to calculate vertical profiles is described in the Supplement, Sect. S1). Afternoon vertical profiles were chosen since this is the time of the day when Black symbols correspond to January and red symbols corresponds to July. Frequency counts in the PBL were only included for the conditions where the mass concentration of total organic aerosol was > 0.5 µg m −3 at the surface. The viscosity parameterization used to calculate mixing times was based on α-pinene SOA generated using mass concentrations of ∼ 1000 µg m −3 .
RH is typically lowest and thus mixing times are the longest. Figure 5c shows that mixing times within α-pinene SOA decrease significantly with altitude for these two locations. This is because the plasticizing effect of water on viscosity dominates the temperature effect for these conditions. Shown in Fig. 6 are global maps of the monthly averaged mixing times of organic molecules within α-pinene SOA for conditions at the top of the PBL for the months of January and July. Figure 6 shows that 92 and 98 % of the locations for January and July, respectively, have a mixing time < 0.1 h for conditions at the top of the PBL when monthly averaged total OA surface concentrations are > 0.5 µg m −3 . Within the PBL, vertical mixing of air masses occurs on the order of 30 min. Since the mixing times within α-pinene SOA particles for conditions at the top of the PBL are < 0.1 h for most locations where the SOA concentrations are significant (total OA concentration > 0.5 µg m −3 at the surface), a reasonable upper limit of the mixing time within the α-pinene SOA studied here for most locations in the PBL is 30 min. During this 30 min interval, mixing times within α-pinene SOA particles can cycle between short and long values, though rarely being > 1 h (Figs. 3 and 4).

Sensitivity analysis
To calculate the mixing times discussed above, we assumed that the α-pinene SOA particles have a diameter of 200 nm. We also repeated these calculations assuming a diameter of 500 nm, since aged OA can have larger diameters (Takegawa et al., 2006). Based on the viscosity parameterization shown in Fig. 2a, mixing times within 500 nm α-pinene SOA particles are < 1 h for 95.9 and 99.4 % of the occurrences in the PBL during January and July, respectively (Fig. S3).
The parameterization of viscosity used above was developed using viscosity measurements by Grayson et al. (2016), Järvinen et al. (2016), and Crittenden et al. (2012). As a sensitivity analysis, we developed a second parameterization, using the same procedure as describe above but using the upper limits to the viscosities reported by Grayson et al. (2016) and the upper limits to the RH ranges reported by Järvinen et al. (2016). This should result in an upper limit to the viscosity parameterization discussed above. The uncertainties in the measurements by Crittenden et al. (2012) were not considered since they are small compared to the uncertainties reported by Grayson et al. (2016) and Järvinen et al. (2016). Based on this second parameterization, mixing times are < 1 h for 96.6 and 99.5 % of the occurrences in the PBL during January and July, respectively, when the total OA was > 0.5 µg m −3 at the surface (Fig. S4).

Effect of mass concentration used to generate the SOA
The parameterizations developed above were based on SOA generated using a mass concentration of ∼ 1000 µg m −3 . As mentioned, we focused on ∼ 1000 µg m −3 because lowtemperature and room-temperature viscosity measurements have been carried out using this mass concentration. However, the viscosity of some types of SOA may depend on the mass concentration used to generate the SOA. For example, Grayson et al. (2016) showed that under dry conditions, the viscosity of α-pinene SOA may increase by a factor of 5 as the production mass concentration decreased from 1200 to 120 µg m −3 . In addition, mass concentrations of biogenic SOA are typically ≤ 10 µg m −3 in the atmosphere (Spracklen et al., 2011). As a starting point to quantify how often mixing times of organic molecules are < 1 h within α-pinene SOA generated using low, atmospherically relevant mass concentrations, we developed a temperature-independent parameterization using the room-temperature viscosity data for αpinene SOA from Zhang et al. (2015) (Table S5) and roomtemperature viscosity data for water from (Crittenden et al., 2012) (Table S3). Zhang et al. (2015) measured the viscosity of α-pinene SOA over a range of relative humidities (0-60 %), and the SOA used in these experiments was generated in the laboratory using a mass concentration of ∼ 70 µg m −3 . The median room-temperature viscosities reported by Zhang et al. are higher than the median room-temperature viscosities reported by Grayson et al. (2016) using a mass concentration of 520 µg m −3 (Fig. S5). Although not proven, a reasonable explanation for the difference in median viscosities is the difference in mass concentrations used to generate the SOA. A temperature-independent parameterization was generated by fitting Eq. (1) to the room-temperature viscosity Figure 5. Temperature, RH, and estimated mixing times for α-pinene SOA as a function of altitude for Hyytiälä (boreal forest) and the Amazon (rainforest). The temperature and RH at ground level are the average afternoon values in the driest month of the year for the respective locations. The vertical profiles of temperature and RH are plotted until the RH is 100 % for these locations. The height at which RH reaches 100 % is only slightly lower than the average height of the planetary boundary layer predicted by GEOS-5 meteorology data. For details see the Supplement, Sect. S1. The viscosity parameterization used to calculate mixing times was based on α-pinene SOA generated using mass concentrations of ∼ 1000 µg m −3 . data from Zhang et al. (2015) and Crittenden et al. (2012), but with the temperature (T ) in Eq. (1) replaced by 293 K. The values for the parameters retrieved by fitting the modified Eq. (1) to the viscosity data are reported in Table S6. The temperature-independent parameterization generated using this method is shown in Fig. 7a. Shown in Figs. 7b, 8a, and b (contours) is the parameterization for mixing times within 200 nm α-pinene SOA based on this temperatureindependent viscosity parameterization. Also included in Fig. 8a and b are the normalized frequency counts of temperature and RH in the PBL for the months of January and July 2006, respectively, when the monthly average mass concentration of total OA was > 0.5 µg m −3 at the surface. Shown in Fig. 9 are the normalized frequency distributions of mixing times within α-pinene SOA for January and July, based on the data in Fig. 8a and b. Figure 9 suggests that the mixing times within α-pinene SOA is < 1 h for 27 and 19.5 % of the occurrences in the PBL during January and July, respectively, when monthly average total OA concentrations were > 0.5 µg m −3 at the surface. However, several caveats need to be emphasized: (1) the parameterization was developed based on room-temperature viscosity data only. Viscosities, and hence mixing times, will increase as temperature decreases. As an illustration, the viscosity of sucrose-water mixtures can increase by 2-3 orders of magnitude as the temperature decreases by 10 K close to the glass transition temperature (Champion et al., 1997).
(2) The mixing times were calculated using the Stokes-Einstein relation, which can underpredict diffusion coefficients, and hence overpredict mixing times, when the viscosity of the matrix is high. For example, the Stokes-Einstein equation can underpredict diffusion coefficients of organic molecules in sucrose-water mixtures by at least a factor of 10 to 100 at viscosities ≥ 10 6 Pa s (Chenyakin et al., 2017;Price et al., 2016). (3) The viscosity data from Zhang et al. (2015) have an uncertainty of ±1 order of magnitude, which was not considered in the temperature-independent parameterization. Considering these caveats, we are unable to make strong conclusions about how often mixing times of organic molecules are < 1 h within α-pinene SOA generated at low mass concentrations. To help resolve this issue, temperaturedependent studies of the viscosity of α-pinene SOA generated using low mass concentrations are needed. Also, the accuracy of the Stokes-Einstein equation for predicting diffusion coefficients of organics within α-pinene SOA needs to be determined.

Mixing times of organic molecules within anthropogenic SOA particles in the PBL
Recently it has been shown that the diffusion rates of organics in SOA from toluene photooxidation are slower than the diffusion rates of organics in SOA from α-pinene ozonolysis at room temperature Song et al., 2016;Ye et al., 2016). These results indicate that mixing times are longer in some types of anthropogenic SOA than some types of biogenic SOA, at least at room temperature. SOA derived from anthropogenic sources can be a significant contributor to SOA over polluted regions Spracklen et al., 2011). Viscosities or diffusion rates within toluene SOA or other types of anthropogenic SOA have yet to be measured at temperatures lower than room temperature. As a result, we have used sucrose as a proxy of anthropogenic SOA since the viscosity of sucrose is similar to the viscosity of toluene SOA at room temperature ( Fig. S6) (Power and Reid, 2014;Song et al., 2016), and since a parameterization of the viscosity of sucrose as a function of temperature and RH can be developed using literature data. In the Supplement (Sect. S2, Table S7-S9, and Figs. S7-S10) we carried out a similar analysis for sucrose as for α-pinene SOA above. Assuming sucrose is a good proxy for anthropogenic SOA, the analysis suggests that 70 and 83 % of the mixing times within anthropogenic SOA in the PBL are < 1 h for January and July, respectively, when SOA concentrations are significant (total OA concentration > 0.5 µg m −3 at the surface). In addition, 81 and 87 % of the locations for January and July, respectively, have a mixing time < 0.1 h at the top of the PBL when surface concentrations of total OA are > 0.5 µg m −3 . These percentages for anthropogenic SOA are likely lower limits since, as mentioned earlier, studies have shown that the Stokes-Einstein relation (which is used here to calculate diffusion coefficients of organic molecules from viscosities) can underpredict diffusion coefficients of organic molecules in sucrose-water mixtures by at least a factor of 10 to 100 at viscosities ≥ 10 6 Pa s (Chenyakin et al., 2017;Price et al., 2016). Measurements of diffusion rates of organic molecules within anthropogenic SOA as a function of both tempera- Figure 7. Plot of RH vs. temperature with contour lines representing (a) the viscosity parameterization for α-pinene SOA particles based on the data from Zhang et al. (2015) and (b) mixing times calculated for organic molecules within 200 nm diameter α-pinene SOA particles. The symbols in (a) represent the laboratory data used to develop the parameterization: the square represents the water viscosity at room temperature from Crittenden et al. (2012) and the circles represent the viscosity data from Zhang et al. (2015). The viscosity parameterization is based α-pinene SOA generated using mass concentrations of ∼ 70 µg m −3 .
ture and RH are needed to better constrain how often mixing times are > 1 h within anthropogenic SOA in the PBL. Shiraiwa et al. (2017) recently estimated mixing times of organics within SOA in the troposphere using a global chemistry climate model and a relationship between glass transition temperatures, molar mass, and oxygen-to-carbon elemental ratios. Their results suggest mixing times of organics within SOA are short (< 1 min) over the oceans, tropics, and high latitudes at the surface and 850 hPa. However, their results suggest mixing times are long (> 1 h) over dry regions (i.e. major deserts) at the surface and at 850 hPa and over most continental regions at 850 hPa. The general trends ob- Figure 8. Six-hour normalized frequency counts of temperature and RH in the planetary boundary layer (PBL) (colour scale) together with the mixing times for organic molecules within 200 nm α-pinene SOA particles (contours) calculated based on the parameterization generated using the viscosities from Zhang et al. (2015). Panel (a) shows the conditions for January and panel (b) shows the conditions for July. Mixing times (contours) are reported in hours. Frequency counts in the PBL were only included for the conditions where the mass concentration of total organic aerosol was > 0.5 µg m −3 at the surface. The viscosity parameterization used to calculate mixing times was based on α-pinene SOA generated using mass concentrations of ∼ 70 µg m −3 . served by Shiraiwa et al. (2017) are consistent with the trends observed here. However, the mixing times predicted by Shiraiwa et al. (2017) appear to be longer than the mixing times predicted here using viscosities of α-pinene SOA generated with a mass concentration ∼ 1000 µg m −3 . Quantitative differences between the current work and the work by Shiraiwa et al. (2017) are not surprising since Shiraiwa et al. (2017) considered both anthropogenic SOA and biogenic SOA simultaneously and since they used a very different approach to estimate viscosities of atmospheric SOA. Figure 9. Normalized frequency distributions of mixing times within α-pinene SOA in the planetary boundary layer (PBL) for the parameterization generated using the viscosity data from Zhang et al. (2015). Black symbols correspond to January and red symbols corresponds to July. Frequency counts in the PBL were only included for the conditions where the mass concentration of total organic aerosol was > 0.5 µg m −3 at the surface. The viscosity parameterization used to calculate mixing times was based on α-pinene SOA generated using mass concentrations of ∼ 70 µg m −3 .

Summary and conclusions
A parameterization for viscosity as a function of temperature and RH was developed for α-pinene SOA based on roomtemperature and low-temperature viscosity data of α-pinene SOA generated in the laboratory using mass concentrations of ∼ 1000 µg m −3 . We focused on ∼ 1000 µg m −3 because low-temperature and room-temperature viscosity measurements have been carried out using this mass concentration. Based on this parameterization, as well as RH and temperatures in the PBL, the mixing times within α-pinene SOA are < 1 h for 98.5 and 99.9 % of the occurrences in the PBL during January and July, respectively, when monthly average total OA concentrations are > 0.5 µg m −3 at the surface. Also based on this parameterization, 92 and 98 % of the locations for January and July, respectively, have a mixing time < 0.1 h for conditions at the top of the PBL when monthly averaged total OA surface concentrations are > 0.5 µg m −3 . As a starting point to quantify how often mixing times of organic molecules are < 1 h within α-pinene SOA generated using low mass concentrations, we developed a temperatureindependent parameterization using the room-temperature viscosity data for α-pinene SOA from Zhang et al. (2015). Zhang et al. (2015) measured the viscosity of α-pinene SOA generated using a mass concentration of ∼ 70 µg m −3 . Based on this temperature-independent parameterization, mixing times within α-pinene SOA are < 1 h for 27 and 19.5 % of the occurrences in the PBL during January and July, respectively, when monthly average total OA concentrations are > 0.5 µg m −3 at the surface. However, several caveats need to be emphasized for these results. Most importantly, the results were based on room-temperature viscosity data only and the mixing times were calculated using the Stokes-Einstein relation, which can underpredict diffusion coefficients of organic molecules, and hence overpredict mixing times, when the viscosity of the matrix is high.
As a starting point to quantify how often mixing times of organic molecules are < 1 h within anthropogenic SOA, a parameterization for viscosity as a function of temperature and RH was developed using sucrose-water viscosity data. Based on this parameterization and assuming sucrose is a good proxy for anthropogenic SOA, 70 and 83 % of the mixing times within anthropogenic SOA in the PBL are < 1 h for January and July, respectively, when SOA concentrations are significant (total OA concentration > 0.5 µg m −3 at the surface). These percentages for anthropogenic SOA are likely lower limits since studies have shown that the Stokes-Einstein relation (which is used here to calculate diffusion coefficients of organic molecules from viscosities) can underpredict diffusion coefficients of organic molecules in sucrose-water mixtures by at least a factor of 10 to 100 at viscosities ≥ 10 6 Pa s (Chenyakin et al., 2017;Price et al., 2016).
To improve the predictions presented above the following are needed: (1) viscosities as a function of temperature and RH for α-pinene SOA and anthropogenic SOA generated using low mass concentrations and (2) studies that quantify the accuracy of the Stokes-Einstein equation for predicting diffusion coefficients in SOA. Studies that explore further the effect of oxidation level, oxidation type, and gas-phase precursor on viscosity and diffusion within biogenic and anthropogenic SOA would also be beneficial.
Data availability. Underlying material and related items for this paper are located in the Supplement.