Digital Commons @ Michigan Tech Digital Commons @ Michigan Tech Evolution of multispectral aerosol absorption properties in a Evolution of multispectral aerosol absorption properties in a biogenically-influenced urban environment during the CARES biogenically-influenced urban environment during the CARES campaign campaign

: We present the evolution of multispectral optical properties through urban aerosols that have aged and interacted with biogenic emissions, resulting in stronger short wavelength absorption and the formation of moderately brown secondary organic aerosols. Ground-based aerosol measurements were made in June 2010 within the Sacramento urban area (site T0) and at a 40-km downwind location (site T1) in the forested Sierra Nevada foothills area. Data on black carbon (BC) and non-refractory aerosol mass and composition were collected at both sites. In addition, photoacoustic (PA) instruments with integrating nephelometers were used to measure spectral absorption and scattering coefﬁcients for wavelengths ranging from 355 to 870 nm. The daytime absorption Ångström exponent (AAE) indicated a modest wavelength-dependent enhancement of absorption at both sites throughout the study. From 22 to 28 June 2010, secondary organic aerosol mass increased signiﬁcantly at both sites, which was due to increased biogenic emissions coupled with intense photochemical activity and air mass recirculation in the area. During this period, the median BC mass-normalized absorption cross-section (MAC) values for 405 nm and 532 nm at T1 increased by ~23% and ~35%, respectively, compared with the relatively less aged urban emissions at the T0 site. In contrast, the average MAC values for the 870 nm wavelength were similar for both sites. These results suggest the formation of moderately brown secondary organic aerosols in biogenically-inﬂuenced urban air.


Introduction
Atmospheric aerosols impact the Earth's radiation budget by scattering and absorbing solar radiation, which is referred to as the aerosol direct effect [1,2]. This can lead to heating or cooling, depending on the ratio of aerosol scattering and extinction coefficients (single scattering albedo, or SSA), the particle scattering asymmetry parameter, and the albedo of the underlying surface [3][4][5][6]. In addition to these direct effects, aerosols can also have indirect effects on the radiation balance by changing the properties of clouds, thus influencing their optical properties and lifetime [7,8]. Inorganic salts such as sulfates are an important aerosol species that strongly scatter and backscatter solar radiation, usually cooling the climate [1,9]. Black carbon (BC) aerosol [10] is the dominant absorber over the solar spectrum [11], and is associated with a warming effect [9,12,13]. Furthermore, a number of studies [14][15][16][17][18][19][20][21] have identified a fraction of organic aerosols (OA) as major absorbers at shorter wavelengths, which are often referred to as "brown carbon" (BrC). OA constitutes a large fraction of the atmospheric aerosol mass; typical measurements range from~20% to 90% of the submicron-diameter aerosol mass [22]. These aerosols scatter light efficiently because their size is comparable to the wavelength of the visible solar radiation, and typical values of the real part of their refractive index are fairly high. Additionally, they can be transported over long distances; hence, their climatic impact is not limited to their source or formation region. However, the role of OA on the Earth's radiation budget is still poorly understood [23].
The Department of Energy's (DOE) Carbonaceous Aerosols and Radiative Effects Study (CARES) field campaign was carried out during June 2010 in the Sacramento (CA, USA) area to better understand the evolution of aerosol properties resulting from the various physical and chemical interactions between anthropogenic and biogenic emissions. A wide range of ground-based and airborne instruments was deployed to monitor comprehensive aerosol properties, meteorological variables, and gaseous species. Two ground sites, one within the Sacramento urban area (T0) and another in Cool, CA,~40 km downwind in the forested Sierra Nevada foothills area (T1), were set up to study the evolution of the aerosols as they aged during transport. Zaveri et al. [24] presented detailed descriptions of the overall research objectives of the CARES campaign, the various sampling platforms and instrumentation, and an overview of several key observations. Fast et al. [25] described the meteorological conditions during CARES.
Previously published aerosol observations during CARES are summarized in the following paragraph. Shilling et al. [26] analyzed aircraft measurements obtained during CARES, and found that the production of secondary organic aerosol (SOA) was dramatically enhanced when anthropogenic emissions from the Sacramento urban area mixed with isoprene-rich air from the Sierra Nevada foothills. A similar conclusion was reached by Setyan et al. [27] and Setyan et al. [28], based on measurements made at the T1 site. Cappa et al. [29] reported that BC absorption enhancement due to coating with non-refractory material during photochemical ageing at T0 was only~6% at 532 nm, and attributed the slightly higher absorption enhancement at 405 nm to the influence of BrC on absorption in this wavelength region. Kassianov et al. [30] demonstrated an unexpectedly large contribution of coarse mode aerosols to aerosol radiative forcing, with coarse mode aerosols sometimes contributing between 50 and 85% of the total aerosol volume observed at T0 and T1 during CARES.
Previous laboratory and field studies employing direct measurements and inference from the use of Mie theory have reported absorption enhancement. As the day progresses, the fractal-like BC particles begin to increase in mass and transform into a more compact structure upon coating by organics, sulfate, ammonium, nitrate, and water [31]. Such coatings on BC particles can lead to absorption enhancement of up to~1.6 [31][32][33], while the AAE can be as high as~1.6 due to non-absorbing coating [18,34]. Laboratory studies have also shown that organic dioctyl sebacate coatings enhanced BC absorption [35].
In this paper, we focus on a comparative analysis of aerosol absorption, scattering, and single scattering albedo over wavelengths ranging from 355 nm to 870 nm at T0 and T1. The objective of this study was to characterize the evolution of these multispectral optical properties as urban aerosols and precursor gases mixed and interacted with biogenic emissions as they were routinely transported to the forested Sierra Nevada foothills area. We begin with a brief description of the measurement sites and instruments, followed by a discussion of the error analysis of aerosol measurements. The BC mass concentration measurements from a single particle soot photometer (SP2) were analyzed together with the aerosol light absorption measurements to obtain BC mass normalized light absorption cross-sections (MAC).

Ground Sites
The CARES campaign took place in the Sacramento valley from 2 to 28 June 2010. Figure 1 shows the location of the ground sampling sites overlaid with the dominant summer daytime wind patterns. Two ground sites were strategically chosen to study the evolution and transport of urban emissions and their interaction with biogenic emissions in the forested area. The urban site, T0 (altitude~30 m above sea level (ASL)), was located on the American River College campus,~14 km northeast of the downtown Sacramento area. The downwind rural site, T1 (altitude~450 m ASL), was located~40 km northeast of the T0 site on the Northside School campus in Cool, CA. In addition to local urban and agricultural sources, emissions from the Bay Area and the petrochemical refineries located in the Carquinez Strait are major sources of air pollution in Sacramento and the Central Valley. During strong southwesterly flow events, coarse-mode sea salt aerosols from the Pacific Ocean can also be advected to the Sacramento area via the Carquinez Strait [24]. non-absorbing coating [18,34]. Laboratory studies have also shown that organic dioctyl sebacate coatings enhanced BC absorption [35].
In this paper, we focus on a comparative analysis of aerosol absorption, scattering, and single scattering albedo over wavelengths ranging from 355 nm to 870 nm at T0 and T1. The objective of this study was to characterize the evolution of these multispectral optical properties as urban aerosols and precursor gases mixed and interacted with biogenic emissions as they were routinely transported to the forested Sierra Nevada foothills area. We begin with a brief description of the measurement sites and instruments, followed by a discussion of the error analysis of aerosol measurements. The BC mass concentration measurements from a single particle soot photometer (SP2) were analyzed together with the aerosol light absorption measurements to obtain BC mass normalized light absorption cross-sections (MAC).

Ground Sites
The CARES campaign took place in the Sacramento valley from 2 to 28 June 2010. Figure 1 shows the location of the ground sampling sites overlaid with the dominant summer daytime wind patterns. Two ground sites were strategically chosen to study the evolution and transport of urban emissions and their interaction with biogenic emissions in the forested area. The urban site, T0 (altitude ~30 m above sea level (ASL)), was located on the American River College campus, ~14 km northeast of the downtown Sacramento area. The downwind rural site, T1 (altitude ~450 m ASL), was located ~40 km northeast of the T0 site on the Northside School campus in Cool, CA. In addition to local urban and agricultural sources, emissions from the Bay Area and the petrochemical refineries located in the Carquinez Strait are major sources of air pollution in Sacramento and the Central Valley. During strong southwesterly flow events, coarse-mode sea salt aerosols from the Pacific Ocean can also be advected to the Sacramento area via the Carquinez Strait [24].  All measurements are reported at ambient temperature and pressure, and in Pacific Daylight Time (PDT), the local time in California during summer.

Measurements of Aerosol Optical Properties
During the CARES campaign, a number of photoacoustic (PA) instruments equipped with integrating reciprocal nephelometers were deployed at the T0 and T1 sites to study the magnitude and spectral dependence of aerosol light absorption (β abs ) and scattering (β sca ) coefficients. To accomplish this, two groups from the University of Nevada, Reno (UNR) and from the Los Alamos National Laboratory (LANL) deployed several PA instruments operating at multiple wavelengths. At T0, wavelengths used included 375, 405, 532, and 870 nm, while at T1, they were 355, 405, 532, and 870 nm.
The PA technique provides real time, in situ measurements of ambient aerosol light absorption [36,37]. In this method, aerosol samples are illuminated in an acoustic resonator by a power-modulated laser. The absorbed energy heats the particles, and the heat is quickly released to the surrounding air, producing a periodic variation in air pressure (i.e., a sound wave). The pressure of the sound wave is measured with a calibrated microphone located in the resonator, thereby producing a first-principle measurement of the β abs . These PA instruments simultaneously measure β sca using a reciprocal nephelometer [38,39].
All PA instruments were calibrated before the CARES campaign using kerosene-flame soot and table salt particles. The absorption and scattering calibration procedure shown and discussed in and around Figure 1 of Gyawali et al. [40] was performed for all of the PA instruments.
Since PA instruments are designed for the concurrent measurement of β abs and β sca , we can obtain the extinction coefficient (β ext = β abs + β sca ), as well as other optical parameters such as the single scattering albedo (SSA) and the absorption Ångström exponent (AAE), through using multiple wavelength measurements. The SSA is the ratio of the β sca and β ext , and can be expressed as The spectral dependence of SSA is related to those of β abs and β sca [41]. At T0, SSA was studied at λ = 405, 532, and 870 nm, and at T1, it was investigated at λ = 355, 532, and 870 nm. The SSA and AAE calculations were reported for the most abundant PA measurements at each site. The AAE expressed by characterizes the spectral dependence of aerosol light absorption assuming a simple power-law dependence between the wavelengths λ 1 and λ 2 [42]. The AAE at T0 was derived from the PA measurements of β abs at 405 and 870 nm, while at T1, it was derived from similar measurements at 355 and 870 nm. These wavelengths were chosen to emphasize light absorption by organic aerosol coatings and intrinsic absorption [18].

Collocated Measurements Used in This Study
BC mass concentrations were measured at T0 and T1 with a single particle soot photometer (SP2) [43] calibrated with Acheson Aquadag [44]. While the SP2 was calibrated with Aquadag, recent studies have shown that diesel soot and wood smoke BC are better represented by fullerene soot [45,46]. Therefore, the SP2-measured BC mass was increased by 43% as a first-order approximation toward a fullerene soot-equivalent value, based on the results of Moteki and Kondo [46].
Within a particle mobility diameter range from 8 to 858 nm, particle size distributions were measured at T1 in 70 logarithmically-spaced size bins by a scanning mobility particle sizer (SMPS) system consisting of a differential mobility analyzer (DMA) and a condensation particle counter (CPC, model 3772; TSI Inc., Shoreview, MN, USA) [27]. Similar measurements were conducted at T0 for particle mobility diameters from 12 to 736 nm in 115 size bins. Measurements of the size-resolved chemical composition of non-refractory species in particles of an aerodynamic diameter less than 1 µm (NR-PM 1 ) were simultaneously acquired with a high-resolution aerosol mass spectrometer (HR-AMS) [27]. The HR-AMS provides mass concentrations of inorganic and organic aerosols.
Aerosol samples for electron microscopy analysis were collected on 13-mm diameter nuclepore polycarbonate membranes with 100-um pore size (Whatman Inc., Chicago, IL, USA) using a custom made sequential sampler, as described in China et al. [47]. The samples were sputter coated with a thin layer of platinum, and were imaged using a Hitachi S-4700 field emission scanning electron microscope (FE-SEM) [24,47].

Error Analysis
We have used basic uncertainty propagation theory to compute the error bars for derived aerosol optical properties, f (β abs , β sca , BC mass ), such as SSA, AAE, and MAC. Assuming that β abs , β sca , and BC mass are uncorrelated, the error δf in f follows from the method of propagation of errors [48] as The three partial derivatives appearing on the right hand side of this expression can be evaluated using Equations (1) and (2), and the well-known expression defining the MAC as MAC = β abs /BC mass . In this study, both instrument calibration uncertainties and statistical variations are considered to estimate the error bars for calculated aerosol optical properties. The terms ∆β abs , ∆β sca , and ∆BC mass correspond to the root square sum of the squares of the uncertainties due to instrument calibrations and the statistical errors. The calibration uncertainties associated with PA absorption coefficient measurements and reciprocal nephelometer (RN) scattering coefficient measurements are 5% and 15%, respectively. Gyawali et al. [40], Lewis [49], Lewis et al. [21], and Arnott et al. [50] have discussed these uncertainties. The uncertainties in BC mass concentrations are about 22% [44]. Additionally, the following equations were used to conservatively estimate the electronic uncertainties associated with measurement precisions of the PA absorption measurements (∆β abs ) precision and SP2 black carbon measurements (∆BC mass ) precision as and In Equations (4) and (5), P L and τ are the laser power and the measurement time, respectively [36]. The uncertainty associated with the electronic precision of SP2 measurements of BC were not considered in the error analysis below due to their negligible contribution to the total MAC error, as the large laser irradiance used by SP2 provides much more precise measurements than those obtained with the PA. The maximum value of τ is 600 s, as the measurement begins to drift after this time. Our PA instruments all use the same type of microphone, acoustical resonator, and detection electronics. The MAC error associated with PA measurement precisions is estimated by using the following equation The statistical errors in regression-derived MAC were obtained from the slope error. Table 1 summarizes the mean, median, and regression values of MAC for 22-28 June 2010 at T0 and T1 including the standard deviation (SD), instrument calibration, and precision error. The total error in the regression MAC values is obtained by using the error propagation technique to combine the slope errors (SE) and calibration uncertainties of the PA and SP2 instruments. Table 1

Results and Discussion
A study of the meteorological parameters characterizing the transport of aerosols during the CARES campaign was published by Fast et al. [25]. In brief, near-surface thermally-driven upslope flows, including synoptic southwesterly (SW) winds, were favorable for transport of the urban plume from Sacramento, CA to the T1 site near Cool, CA. These patterns were occasionally interrupted by a northwesterly (NW) flow, which carried the urban plume southeast to the San Joaquin Valley (Figure 1), thereby reducing the interaction of the urban plume with biogenic emissions at T1 and providing recirculation to T0.

Daily Average of Aerosol Chemical, Physical, and Optical Properties
The left panels of Figure 2 show the daily daytime (06:00 to 18:00 PDT) average of PM 1 organic mass concentration, the PA measurements of β abs and β sca , AAE, SSA, and the particle mobility diameter (particle mean number-weighted diameter) measured by the SMPS at T0 from 12 to 28 June 2010. Likewise, the right panels of Figure 2 present similar measurements for the T1 site, from 2 to 28 June 2010. The vertical bars in PM 1 organic mass concentration, β abs , β sca , and PM 1 mean diameter represent one standard deviation of the averaged measurements during daylight hours. The uncertainties of SSA and AAE were determined by the error propagation technique, which is described in Section 2.4. PM 1 organic mass concentration measurements at T0 were available only from 17 to 25 June 2010. The gray bar at the top of Figure 2 indicates periods during which the synoptic wind was either southwesterly (SW) or northwesterly (NW). Semi-transparent gray shading for the NW flow periods is also shown over all the plots.
Site T0: The daily daytime BC mass concentrations varied between 0.05 ± 0.02 (average ± 1σ) and 0.25 ± 0.09 µg m -3 , with an overall average of 0.14 ± 0.05 µg m -3 (Figure 2a). The daily daytime average β abs plotted for 405 and 870 nm (Figure 2c) exhibited similar variations with BC mass concentrations and PM 1 organic mass concentrations (Figure 2a). The peaks of β abs~5 .5 ± 2.8 and 1.6 ± 0.6, and 6.0 ± 2.4 and 1.9 ± 0.5 Mm -1 at 405 and 870 nm, respectively were observed on 14 and 17 June 2010, and found to be associated with the corresponding peak values of BC mass concentrations.
Scattering at 405 nm ( Figure 2g) correlated especially well with the mean diameter of PM 1 from after 22 June 2010 to the end of the campaign, when PM 1 particle growth contributed to enhanced scattering at 405 nm relative to 870 nm. The daily daytime maximum PM 1 mean mobility diameter (Figure 2k) was 68.5 ± 11.5 nm with β sca of 57 ± 20 Mm −1 for 28 June, while the daily daytime β sca at 870 nm remained constant at~11 Mm −1 . The SSA (Figure 2i) was lower for both 405 and 870 nm during NW flows, due to the increased influence of fresh traffic emissions. The daily daytime SSA varied between 0.72 ± 0.31 (average ± uncertainty) and 0.95 ± 0.03 at 405 nm, and 0.81 ± 0.04 and 0.97 ± 0.01 at 870 nm. The difference between the SSA at 405 and 870 nm decreased continuously as the campaign progressed due to the enhanced formation of SOA and recirculation of aged aerosols. The SSA exhibited a value of~0.90 for both 405 and 870 nm after 25 June. The daily daytime average AAE (Figure 2e) at 405 and 870 nm remained fairly constant, as it varied only between 1.1 ± 1.6 (±uncertainty) and 1.7 ± 1.0 with an average value of~1.6 ± 0.3, indicating enhanced absorption at shorter wavelengths at T0, even in background air.   (g-h) aerosol scattering coefficient (βsca), (i-j) single scattering albedo (SSA), and (k-l) PM1 particle mean mobility diameter. Vertical bars (Figure 2a, b, c, d, g, h, k, and l) represent one standard deviation (1σ) around the mean calculated for each vertical level. The error bars in AAE and SSA (Figure 2e, f, i, and j) correspond to the total uncertainty, the square root of the sum of squares of the 1σ around the mean, and the propagated error from the systematic errors (5% and 15% relative uncertainties in photoacoustic (PA) βabs measurements and reciprocal nephelometer (RN) βsca measurements, respectively). Each data point is averaged from 06:00 to 18:00 Pacific Daylight Time (PDT). The gray bar at the top indicates periods during which the synoptic wind was either southwesterly (SW) or northwesterly (NW). Semi-transparent gray shading for the NW flow periods is also shown over all the plots. ; (e,f) absorption Ångström exponent (AAE); (g,h) aerosol scattering coefficient (β sca ); (i,j) single scattering albedo (SSA); and (k,l) PM 1 particle mean mobility diameter. Vertical bars (Figure 2a-d,g,h,k,l) represent one standard deviation (1σ) around the mean calculated for each vertical level. The error bars in AAE and SSA (Figure 2e,f,i,j) correspond to the total uncertainty, the square root of the sum of squares of the 1σ around the mean, and the propagated error from the systematic errors (5% and 15% relative uncertainties in photoacoustic (PA) β abs measurements and reciprocal nephelometer (RN) β sca measurements, respectively). Each data point is averaged from 06:00 to 18:00 Pacific Daylight Time (PDT). The gray bar at the top indicates periods during which the synoptic wind was either southwesterly (SW) or northwesterly (NW). Semi-transparent gray shading for the NW flow periods is also shown over all the plots.  (Figure 3i). Freshly emitted BC particles can be mixed with partly coated BC particles (Figure 3ii), or coated by other material (Figure 3iii) (darker region), such as unburned fuel, lubricating oil, or organic matter [47,51]. Figure 3 (left image) shows a field emission scanning electron microscope (FE-SEM) image of particles collected at the T0 site on a nuclepore membrane. The image shows a broad size distribution of BC aggregates as large as 4.5 µm in the longest dimension (Figure 3i). Freshly emitted BC particles can be mixed with partly coated BC particles (Figure 3ii), or coated by other material (Figure 3iii) (darker region), such as unburned fuel, lubricating oil, or organic matter [47,51]. Site T1: The daily daytime average BC mass concentration varied from 0.02 ± 0.01 to 0.14 ± 0.08 µg m −3 (Figure 2b), with an overall average of 0.06 ± 0.02 µg m −3 . The lower BC mass concentrations were primarily observed during the NW flows. As mentioned previously, the NW flows disrupted the transport of the urban plume from Sacramento to the T1 site. Between 8-11 June 2010, the average mass concentration of PM1 organics decreased from 3.80 ± 0.40 (±1σ) to 0.80 ± 0.13 µg m −3 , while the mean mobility diameter of these particles decreased from 65 ± 20 to 39 ± 12 nm. In Figure  2b, the time series of total PM1 organic mass concentration has been shown in more detail, with hydrocarbon-like organic aerosol (HOA) and the sum of oxygenated organic aerosol (OOA) mass concentrations plotted individually; these data show that most of the organic aerosol mass fraction observed corresponds to SOA. While this does not give any indication on the ageing of particles, we also note that according to Setyan et al. [27,28], (1) the toluene/benzene mass ratio was between 1.5 and 2.0, suggesting moderate ageing, (2) the O/C ratios of organic aerosol were between 0.35 and 0.45, which also suggest moderate ageing, and (3) the change in organic aerosol (OA) mass concentration relative to that of CO (ΔOA/ΔCO) was the highest when air masses were dominated by anthropogenic emissions in the presence of a high concentration of biogenic volatile organic compounds, indicating enhanced SOA formation in mixed anthropogenic-biogenic plumes. Additional information relevant to aerosol ageing includes the meteorological modeling described by Fast et al. [25] and the temperature time series given by Zaveri et al. [24].
The daily daytime average AAE (Figure 2f) varied from 1.3 ± 0.6 to 2.3 ± 0.8, and remained ~1.8 after 22 June, suggesting enhanced absorption at the shorter wavelengths, which is likely due to coating by mildly absorbing organic species on BC. Higher temperatures result in higher emission Site T1: The daily daytime average BC mass concentration varied from 0.02 ± 0.01 to 0.14 ± 0.08 µg m −3 (Figure 2b), with an overall average of 0.06 ± 0.02 µg m −3 . The lower BC mass concentrations were primarily observed during the NW flows. As mentioned previously, the NW flows disrupted the transport of the urban plume from Sacramento to the T1 site. Between 8 and 11 June 2010, the average mass concentration of PM 1 organics decreased from 3.80 ± 0.40 (±1σ) to 0.80 ± 0.13 µg m −3 , while the mean mobility diameter of these particles decreased from 65 ± 20 to 39 ± 12 nm. In Figure 2b, the time series of total PM 1 organic mass concentration has been shown in more detail, with hydrocarbon-like organic aerosol (HOA) and the sum of oxygenated organic aerosol (OOA) mass concentrations plotted individually; these data show that most of the organic aerosol mass fraction observed corresponds to SOA. While this does not give any indication on the ageing of particles, we also note that according to Setyan et al. [27,28], (1) the toluene/benzene mass ratio was between 1.5 and 2.0, suggesting moderate ageing; (2) the O/C ratios of organic aerosol were between 0.35 and 0.45, which also suggest moderate ageing; and (3) the change in organic aerosol (OA) mass concentration relative to that of CO (∆OA/∆CO) was the highest when air masses were dominated by anthropogenic emissions in the presence of a high concentration of biogenic volatile organic compounds, indicating enhanced SOA formation in mixed anthropogenic-biogenic plumes. Additional information relevant to aerosol ageing includes the meteorological modeling described by Fast et al. [25] and the temperature time series given by Zaveri et al. [24].
The daily daytime average AAE (Figure 2f) varied from 1.3 ± 0.6 to 2.3 ± 0.8, and remained 1.8 after 22 June, suggesting enhanced absorption at the shorter wavelengths, which is likely due to coating by mildly absorbing organic species on BC. Higher temperatures result in higher emission rates for biogenic VOCs (e.g., isoprene, monoterpenes) [53]. Enhanced SOA formation was observed when the urban plume interacted with biogenic emissions during the campaign [27]. The steadily increasing organic mass from 22 to 28 June at T1 (Figure 2b) produced a strong impact on aerosol light scattering at 355 nm (Figure 2h) with~84% increase in β sca , while β abs at 355 nm increased bỹ 35%. For comparison, during the same period, β sca and β abs at 870 nm increased by~5% and~11%, respectively. These results suggest that the SOA formed from the oxidation of mixed biogenic and anthropogenic emissions in the forested region (near T1) developed a modest apparent or intrinsic absorption at 355 nm, which is much less than the strong intrinsic SOA ultraviolet (UV) absorption observed for biomass burning aerosols [18]. Furthermore, at the T1 site, OOA comprised 90% of the total mass of organics, while HOA comprised only 10% [27]. The OOA fraction corresponds to SOA, and the HOA fraction corresponds to traffic-related primary emissions. From this PM 1 classification, it is apparent that the enhanced scattering of 355 nm radiation at T1 was dominated by the OOA components of the SOA.
Emissions of VOC precursors for the observed SOA have been discussed by Zaveri et al. [24]; maps of anthropogenic VOC and biogenic isoprene emission rates along with locations of T0 and T1 sites are shown in Figure 2 [24]. In addition, Shilling et al. [26] analyzed airborne mass spectrometer measurements made during mornings and afternoons, and observed that isoprene concentrations of up to 12.9 ppbv dominated the biogenic VOC concentrations, with monoterpene concentrations generally below the detection limit of~300 pptv. Setyan et al. [27] noted that at T1, monoterpene concentrations have a very strong diurnal pattern; they are close to the detection limit during the day, and increase up to 1 ppb during the night, which is still five times lower than the isoprene concentration. A possible explanation for this result is that the main sources of monoterpenes are quite far from T1, and given that monoterpenes are very reactive with OH radicals, they may be already almost completely oxidized before reaching T1. Another possible explanation is that the main sources of monoterpenes are located northeast of T1. In this case, monoterpenes would be observed at T1 only during the night, when we have the recirculation of air from the foothills of the Sierra Nevada back to T1.
Days with NW flows exhibited lower SSA (Figure 2j) in comparison with days with SW flows both at 355 and 870 nm, with values 0.72 ± 0.13 (±uncertainty) and 0.80 ± 0.20, respectively. As observed at the T0 site, the SSA at both wavelengths became increasingly similar after 22 June, with values~0.90. This was consistent with the concurrently observed increasing mass of PM 1 organics and the diameter of the particles. The aerosol size increases during atmospheric processing as a result of the coagulation and condensation of both organic and inorganic species, which affect aerosol optical properties [54]. After 22 June, the growing mean diameter of PM 1 (Figure 2l) resulted in enhanced β sca at 355 nm at T1, and also at 405 nm at T0, though the effect was more pronounced at T1. The maximum PM 1 mean diameter at T1 reached~98 ± 8 nm, while β sca at 355 nm peaked at 54.0 ± 8.8 Mm −1 on 28 June. The β sca at 870 nm was slightly higher after 22 June, but remained nearly constant at~10 Mm −1 . Figure 3 (right image) shows a FE-SEM image of particles collected at the T1 site, and displays: (a) fractal-like BC aggregates with open ( Figure 3i) and compact structures (Figure 3ii), which are consistent with the presence of fresh and aged particles; (b) partially-encapsulated BC with a spherical particle (Figure 3iii), which are representative of the complex mixing state of BC particles; (c) spherical and near-spherical particles (Figure 3iv), which are consistent with the presence of organic aerosols. For comparison, Chakrabarty et al. [55] reported the morphology of incense smoke aerosols composed of mostly BrC to be non-coalescing and weakly-bound aggregates.

Aerosol Light Absorption Efficiency
Aerosol light absorption efficiency is obtained from the ratio of aerosol light absorption coefficients measured with the PA instruments and BC mass concentrations measured with the SP2, and is referred to as mass-normalized absorption cross-section (MAC). MAC is of great importance in understanding the role of BC in climate forcing. Estimations of MAC and the associated uncertainties are presented in this section.

MAC from Linear Regression, Mean, and Median Analyses
Results are presented for the period from 22 to 28 June, as this was the longest time period during which all SP2 and PA instruments were operating continuously at both sites. As mentioned previously, these days were favorable for the transport of the urban plume from Sacramento to the T1 site and its vicinity. Additionally, a steady buildup of aged pollutants was noted.
Regression analyses were carried out for PA measurements of absorption coefficients and BC mass concentrations from the SP2. The slope of the linear regression yields an estimate of BC MAC (m 2 g −1 ). Better correlation between the absorption coefficients and BC mass concentration is generally expected for the longer wavelengths (e.g., 532 and 870 nm) than for the shorter wavelengths (e.g., 355, 375, and 405 nm), as light absorption by the latter is potentially influenced by the presence of BrC. Perfect correlation is not necessarily expected for any wavelength due to the increase in MAC associated with particle coatings and BrC, and the decrease in MAC associated with the shielding effect on light absorption of the aerosol interior for larger aerosol particles. Measurement precision also affects the correlation.
Error bars for MAC were calculated from the propagation of the standard error of the slope, instrument precision, and calibration uncertainty for both PA and SP2 instruments. MAC values for T0 and T1, from near UV to near IR wavelengths, are shown in Figure 4. Mean and median MAC values are also presented in the same Figure. These mean and median MAC values were derived by dividing continuous measurements of absorption coefficients by continuous measurements of BC mass concentrations. Mean and median values are quite similar for all wavelengths except 355 nm, suggesting that MAC measurements are symmetrically distributed around the mean values for these wavelengths, as discussed further in the following section. The linear regression model assumes that the MAC is the same for all points, and is only a coarse statistical measure of the relationship between BC and light absorption; this is especially true for wavelengths where BrC absorbs strongly as well. We prefer to use the median MAC values, as in Doran et al. [56,57], when reporting a single value for MAC at a specific wavelength. MAC values obtained from the slope of the regression analysis are substantially lower for the shorter wavelengths compared with the mean and median values, because offsets (due to lack of correlation BC and BrC) will positively affect the mean and median values. Table 1 summarizes the linear regression, mean, and median MAC values with the corresponding slope errors, standard deviations, and instrumental errors. This table represents our current state of knowledge for MAC estimations for places such as Sacramento in early summer, where signal levels are modest and measurement uncertainty is similar to the range of MAC values caused by various physical and chemical processes. One use of this table by modelers could be to assume that the MAC values are normally distributed (see next section) about the mean values with a width (due to chemical and physical processes) that is equal to the square root of the standard deviation squared (last column), minus the precision error squared (column 6). In doing so, note that the width of the MAC distribution functions are generally larger (except for 870 nm) at T1 compared with T0, perhaps suggesting that BrC contributes more to light absorption at T1. Note that the slope error from linear regression is proximally close in values to the precision error calculated from the use of Equation (6). The large standard deviation at 355 nm really stands out compared with other values, perhaps indicating that both dust and brown carbon light absorption contribute to the width of the MAC probability distribution, substantially more so than at 375 nm. The median MAC values for T1 ranged from 32.2 ± 18.9 (±standard deviation) at 355 nm to 7.2 ± 0.5 m 2 g −1 at 870 nm, while for T0, the values ranged from 15.3 ± 7.6 m 2 g −1 at 375 nm to 8.2 ± 1.6 m 2 g −1 at 870 nm. Even though the vertical bars (standard deviation) in the average MAC values at T0 and T1 for similar wavelengths were touching or overlapping, the median MAC values at T0 were lower than those at T1. For example, the median MAC values at T1 were about 23% and 35% higher than those at T0 for 405, and 532 nm, respectively. Brown carbon could account for ~10-30% of total absorption at mid-visible to UV wavelengths [58]. Our MAC values at 870 nm for T0 and T1 are larger than Bond and Bergstrom's [11] suggested value of 7.5 ± 1.2 m 2 g −1 for freshly emitted BC at 550 nm (or 4.7 m 2 g −1 at 870 nm, assuming a λ −1 dependence). Using a PA and an organic carbon/elemental carbon OC/EC analyzer instrument, Doran [56,57] found similar MAC values of 5.6 and 5.5 m 2 g − 1 at 870 nm for partially coated BC at sites T1 (north of Mexico City) and T2 (~35 km farther to the northeast) downwind of Mexico City. It is possible that the differences in the determination of BC from the SP2 and EC from the Sunset Labs OC/EC analyzer causes large MAC value differences between Sacramento and Mexico City, as the PA 870 nm instrument was the same in both studies.
The rather small difference in absorption at 870 nm at T1 and T0 suggests that the additional coating on the aged BC particles at T1 may not have produced an appreciable absorption enhancement. However, the observed increases in absorption at 532 nm, 405 nm, and UV wavelengths at T1 may be mainly due to the influence of brown SOA and perhaps mineral dust [59,60]. While these results are generally consistent with the findings of Cappa et al. [29], a more detailed analysis of the absorption data with the core-shell Mie theory, constrained with the observed BC coating thickness (from SP2 measurements) and morphological information (based on SEM images), is needed to estimate the relative contributions of the absorption enhancement and The median MAC values for T1 ranged from 32.2 ± 18.9 (±standard deviation) at 355 nm to 7.2 ± 0.5 m 2 g −1 at 870 nm, while for T0, the values ranged from 15.3 ± 7.6 m 2 g −1 at 375 nm to 8.2 ± 1.6 m 2 g −1 at 870 nm. Even though the vertical bars (standard deviation) in the average MAC values at T0 and T1 for similar wavelengths were touching or overlapping, the median MAC values at T0 were lower than those at T1. For example, the median MAC values at T1 were about 23% and 35% higher than those at T0 for 405, and 532 nm, respectively. Brown carbon could account for 10-30% of total absorption at mid-visible to UV wavelengths [58]. Our MAC values at 870 nm for T0 and T1 are larger than Bond and Bergstrom's [11] suggested value of 7.5 ± 1.2 m 2 g −1 for freshly emitted BC at 550 nm (or 4.7 m 2 g −1 at 870 nm, assuming a λ −1 dependence). Using a PA and an organic carbon/elemental carbon OC/EC analyzer instrument, Doran [56,57] found similar MAC values of 5.6 and 5.5 m 2 g −1 at 870 nm for partially coated BC at sites T1 (north of Mexico City) and T2 (~35 km farther to the northeast) downwind of Mexico City. It is possible that the differences in the determination of BC from the SP2 and EC from the Sunset Labs OC/EC analyzer causes large MAC value differences between Sacramento and Mexico City, as the PA 870 nm instrument was the same in both studies.
The rather small difference in absorption at 870 nm at T1 and T0 suggests that the additional coating on the aged BC particles at T1 may not have produced an appreciable absorption enhancement. However, the observed increases in absorption at 532 nm, 405 nm, and UV wavelengths at T1 may be mainly due to the influence of brown SOA and perhaps mineral dust [59,60]. While these results are generally consistent with the findings of Cappa et al. [29], a more detailed analysis of the absorption data with the core-shell Mie theory, constrained with the observed BC coating thickness (from SP2 measurements) and morphological information (based on SEM images), is needed to estimate the relative contributions of the absorption enhancement and brown SOA to the total multispectral absorption observations at both T0 and T1. Such a study is planned, and results will be published elsewhere. brown SOA to the total multispectral absorption observations at both T0 and T1. Such a study is planned, and results will be published elsewhere.  The Monte Carlo method was used to simulate the MAC values obtained from the regression analysis presented in section 3.2.1. The simulation depicts the spread of the MAC values with instruments that suffer from error due to measurement precision alone. The algorithm for simulating the experiment is as follows: first, randomly choose 〈BC〉 values in the range from BCmin to BCmax, then calculate the corresponding 〈βabs〉 value from 〈βabs〉 = slope〈BC〉, where the slope is the MAC value from the regression analysis. Second, choose a trial measurement of βabs by picking a random number from the range = −4 ∆ + , 4 ∆ + . The probability distribution in 〈βabs〉 is normalized, and the peak value of the probability distribution P max is calculated by using the equation = √2 ∆ . A random number is chosen; if this random number is ≤P max , then accept as βabs. This process is repeated for finding the same number of pairs of acceptable values of βabs and BC, as reported in the measured values in the histograms. The Monte Carlo method was used to simulate the MAC values obtained from the regression analysis presented in Section 3.2.1. The simulation depicts the spread of the MAC values with instruments that suffer from error due to measurement precision alone. The algorithm for simulating the experiment is as follows: first, randomly choose BC values in the range from BC min to BC max , then calculate the corresponding β abs value from β abs = slope BC , where the slope is the MAC value from the regression analysis. Second, choose a trial measurement of β abs by picking a random number from the range β trial abs = [−4 ∆β abs + β abs , 4 ∆β abs + β abs ]. The probability distribution in β abs is normalized, and the peak value of the probability distribution P max is calculated by using the equation P max = √ 2π ∆β abs −1 . A random number is chosen; if this random number is ≤P max , then accept β trial abs as β abs . This process is repeated for finding the same number of pairs of acceptable values of β abs and BC, as reported in the measured values in the histograms.

MAC Frequency Distribution and Monte Carlo Simulation
The modeled MAC distributions are very narrow in comparison to the corresponding observed distributions. This was expected, since the modeled MAC values only consider uncertainty in measurement precision. Additionally, the shortwave lengths (355 and 375 nm) show a very broad histogram for observed MAC, while for longer wavelengths, the histograms are narrower, and the data move to lower MAC values, as do the frequency peaks. Furthermore, for similar wavelengths, the MAC distribution shifted to the higher values for T1 than for T0, implying a prevalence of BrC at T1.

Conclusions
As part of the CARES field campaign carried out in June 2010 in Sacramento Valley, California, the evolution of multispectral aerosol optical properties was characterized for the Sacramento urban area (site T0) and a 40-km downwind location (site T1) in the forested Sierra Nevada foothills area, as the urban aerosols aged and interacted with biogenic emissions. The weather conditions at the beginning of the campaign were spring-like (high wind speeds and mild temperatures), in which aerosol optical properties were mostly driven by primary emissions, and later transitioned to summer-like conditions (low wind speeds and higher temperatures), during which optical properties were increasingly influenced by secondary aerosol formation. Enhanced production of SOA was observed during the southwesterly flow events, during which the average SSA remained high (>0.90) at both sites. In contrast, the SSA values were appreciably lower during the occasional northwesterly flow events (as seen in Figure 3) when both sites were impacted by fresh local emissions of BC.
During the last week of the campaign, under sustained southwesterly flow, the SOA mass increased significantly at both sites, which was due to increased biogenic emissions coupled with intense photochemical activity and air mass recirculation in the area. This resulted in a strong impact on 355 nm light scattering due to larger submicron particles. The month-long daytime average AAE was 1.6 for the wavelength pair 405 and 870 nm at T0, while it was~1.8 for the wavelength pair 355 and 870 nm at T1, indicating a modest wavelength-dependence of absorption at both sites throughout the study. During the last week of the campaign (i.e., the period when SOA mass increased significantly), the BC mass-normalized absorption cross-section (MAC) for 405 nm at T1 increased by~23% compared with that of the relatively less aged urban emissions at the T0 site. In contrast, the MAC values for the 870 nm wavelength were nearly identical at T0 and T1. These results suggest the formation of moderately brown SOA in biogenically-influenced urban air, as opposed to absorption enhancement caused by non-absorbing SOA coating. These results are generally consistent with the conclusions of Cappa et al. [29]. However, a detailed model analysis of the absorption data, constrained with the observed BC coating thickness and morphological information, is needed to estimate the relative contributions of brown SOA and the absorption enhancement caused by non-absorbing SOA coating to the total multispectral absorption observations at both sites. and manuscript preparation. All authors were given the opportunity to write sections of the manuscript concerning their measurements and interpretations, and to comment on the overall manuscript.

Conflicts of Interest:
The authors declare no conflict of interest.