Light absorption of brown carbon in eastern China based on 3-year multi-wavelength aerosol optical property observations and an improved absorption Ångström exponent segregation method

Brown carbon (BrC), a certain group of organic carbon (OC) with strong absorption from the visible (VIS) to ultraviolet (UV) wavelengths, makes a considerable contribution to light absorption on both global and regional scales. A high concentration and proportion of OC has been reported in China, but studies of BrC absorption based on long-term observations are rather limited in this region. In this study, we reported 3-year results of light absorption of BrC based on continuous measurement at the Station for Observing Regional Processes of the Earth System (SORPES) in the Yangtze River Delta, China, combined with Mie theory calculation. Light absorption of BrC was obtained using an improved absorption Ångström exponent (AAE) segregation method. The AAE of non-absorbing coated black carbon (BC) at each time step is calculated based on Mie theory simulation, together with single particle soot photometer (SP2) and aethalometer observations. By using this improved method, the variation of the AAE over time is taken into consideration, making it applicable for long-term analysis. The annual average light absorption coefficient of BrC (babs_BrC) at 370 nm was 6.3 Mm−1 at the SORPES station. The contribution of BrC to total aerosol absorption (PBrC) at 370 nm ranged from 10.4 to 23.9 % (10th and 90th percentiles, respectively), and reached up to∼ 33 % in the openbiomass-burning-dominant season and winter. Both babs_BrC and PBrC exhibited clear seasonal cycles with two peaks in later spring/early summer (May–June, babs_BrC ∼ 6 Mm−1, PBrC ∼ 17 %) and winter (December, babs_BrC ∼ 15 Mm−1, PBrC ∼ 22 %), respectively. Lagrangian modeling and the chemical signature observed at the site suggested that open biomass burning and residential coal/biofuel burning were the dominant sources influencing BrC in the two seasons, respectively.


Introduction
Atmospheric aerosols not only cause adverse impacts on human health, but they also alter earth's radiation balance through their strong light scattering and absorption, substantially influencing regional and even global climate change (IPCC, 2013;Dockery et al., 1993;Wang et al., 2017b).Light absorption of aerosols strongly influences the magnitude and sign of radiative transfer.Black carbon (BC) and dust have been considered to be two dominant contributors to aerosol absorption extinction.However, a certain type of organic aerosol defined as brown carbon (BrC) was revealed to be of strong light-absorbing efficiency (Formenti et al., 2003;Pöschl, 2003;Andreae and Gelencser, 2006;Kirchstetter et al., 2004;Mukai and Ambe, 1986;Patterson and McMahon, 1984), which can cause perturbations on radiation transfer similar to BC.This implies that the aerosol cooling effect could be overestimated by ignoring the light absorption from Published by Copernicus Publications on behalf of the European Geosciences Union.
BrC.Its strong light absorption in the ultraviolet (UV) range can also affect atmospheric oxidizing capacity by restraining the photolysis rates for photochemically active gases (Laskin et al., 2015).It has been reported that the radiative forcing by BrC globally is about one-quarter of that by BC (Feng et al., 2013), while regional radiative forcing by BrC in areas with intensive combustion activities (e.g., South and East Asia, South America and Africa) can be much higher than the global average, indicating the substantial contribution by BrC to aerosol light absorption in these regions and thereby its significant influence on regional climate change.Therefore, investigation on the BrC contribution to light absorption is of great importance to reduce the uncertainties of the aerosol-radiation interaction (ARI) estimation.
BrC is defined as a type of organic carbon which can absorb solar radiation efficiently in the near-UV (300-400 nm) to the visible (VIS) range (Formenti et al., 2003;Pöschl, 2003;Andreae and Gelencser, 2006;Kirchstetter et al., 2004;Mukai and Ambe, 1986;Patterson and McMahon, 1984).It is a group of species with specific physical properties but detailed chemical components that are difficult to characterize.BrC can be produced not only from primary emissions relating to biomass burning (BB) and fuel combustion, but also from secondary organic aerosol (SOA) formation through the oxidation of volatile organic compounds (VOCs) (Andreae and Gelencsér, 2006;Saleh et al., 2014;Chen and Bond, 2010;Lack et al., 2012;Laskin et al., 2015;Healy et al., 2015).Some studies reported that vehicle and ship emissions may also be sources of BrC (Stone et al., 2009;Cavalli et al., 2004).The complication of the emitted mixtures containing BC, BrC, non-absorbing organic aerosol and inorganic materials in different proportions makes it difficult to perform source attribution and to estimate its emission factor.
Light absorption of BrC is usually estimated based on its strong wavelength dependence with higher absorption from the VIS to the UV range.The wavelength dependence of the aerosol absorption coefficient (b abs ) is normally represented by the absorption Ångström exponent (AAE), and its relationship with b abs is b abs (λ) ∝ λ −AAE (Moosmüller et al., 2011;Sun et al., 2007).Based on the difference of wavelength dependence for BC and BrC, previous studies segregated light absorption of BrC from multi-wavelength optical measurements (Lack and Langridge, 2013;Mohr et al., 2013;Shen et al., 2017), which is called the AAE method.Earlier, the similar concept was used to segregate carbonaceous aerosol fractions from different emission sources based on their difference in AAE value (e.g., wood burning and traffic emission) (Sandradewi et al., 2008;Healy et al., 2012).Usually, the AAE of BC (AAE BC ) was set to be 1.0 based on the properties of bulk BC by many previous studies, assuming the unity of AAE BC between any two wavelengths within the UV to the near-infrared (IR) range (Shen et al., 2017;Olson et al., 2015;Lack and Langridge, 2013).However, large uncertainty was mentioned by Lack and Langridge (2013) when using this assumption, and also this method can only be used when the proportion of BrC is high.Moreover, it is noticed that the AAE of BC is not always 1.0; instead, it can be affected by particle size and mixing state (Lack and Cappa, 2010;Liu et al., 2015Liu et al., , 2014)).
China, the largest country in East Asia, with tremendous fossil fuel and biofuel consumption and extensive agricultural burning, is of great concern in terms of its large contribution of carbonaceous aerosols including light-absorbing carbon (BC, BrC) (Ding et al., 2016a).The BC issue has been noticed in recent years in China, including its temporal variations, emission sources and climate effect (Cao et al., 2006(Cao et al., , 2009;;Yang et al., 2009;Wang et al., 2017a;Huang et al., 2014) as well as its "dome effect" in modifying the boundary layer and enhancing haze pollution in megacities (Ding et al., 2016b;Wang et al., 2018).Organic matter (OM) is a large contributor to PM 2.5 in China (15-51 % of PM 2.5 ; Wang et al., 2017c), especially in the Yangtze River Delta (YRD) region, which is one of the most densely populated city clusters in eastern China and also an important agricultural center, with crops planted in both cold and warm seasons (Ding et al., 2013a).In the YRD region, the OM fraction is 20-40 % of PM 2.5 (Wang et al., 2017c, b) due to the influence of complicated combustion sources; therefore, it is a highly possible that BrC is a contributor to aerosol light absorption in this region.However, studies concerning BrC in China, especially in the YRD region, have still been limited up to now.Yuan et al. (2016) reported light absorption contributions of BrC in the Pearl River Delta (PRD) region to be 6.3 to 12.1 % at 405 nm in autumn and winter campaigns by using the AAE method.Shen et al. (2017) conducted light absorption measurements in Xi'an in summer and winter, and reported the average b abs_BrC_370 of 6.4 and 43.0 Mm −1 in the two seasons, respectively, based on the AAE method.Yang et al. (2009) calculated the mass absorption efficiency of BrC at 370 nm to be 2.2 m 2 g −1 and the AAE of BrC to be 3.5, respectively, during a campaign in March 2015 at Xianghe in northern China.It is noticed that long-term investigation of BrC light absorption in the YRD region has not been reported yet.
This study intends to provide a comprehensive analysis of light-absorbing characteristics of BrC and its influencing factors based on field measurements of multi-wavelength aerosol light absorption at the Station for Observing Regional Processes of the Earth System (SORPES) station in Nanjing in the western YRD region and an improved AAE segregation method based on Mie theory simulation and BC single particle measurements.The contribution of BrC to light absorption is estimated based on 3-year observations, and the potential sources of BrC in typical seasons are discussed.

Measurements of optical properties and relevant species
Field observation was conducted at the SORPES station, located on top of a small hill (118 • 57 10 E, 32 • 07 14 N; 40 m a.s.l.) in the Xianlin campus of Nanjing University (Shen et al., 2018;Ding et al., 2013cDing et al., , 2016a)).The SORPES station is a research and experiment platform in the western YRD, which is influenced by intensive human activities (Ding et al., 2016a).This observation site is about 20 km east of the city center.The map and emission characteristics surrounding the site were given in our previous works (e.g., Ding et al., 2013cDing et al., , 2016a)).Since this observation site is generally upwind from the city center and influenced by the East Asian monsoon, this site is generally downwind of a densely populated city cluster including the megacity of Shanghai.Therefore, this station can be considered as a regional background station in the western YRD region and is an ideal site to study the impact of multiple anthropogenic emissions on regional air quality in eastern China (Ding et al., 2016a(Ding et al., , 2013c)).Aerosol optical properties were measured from 1 June 2013 to 31 May 2016.Aerosol light absorption measurements were conducted using a multi-wavelength aethalometer (Model AE-31, Magee Scientific Company Berkeley, California, USA), which performs continuous measurements at seven wavelengths, i.e., 370 nm in the UV wavelength range, 470, 520, 590, 660 and 880 nm in the VIS wavelength range and 950 nm in the IR wavelength range.The time resolution is 5 min.Sample air was obtained through a stainlesssteel inlet with a PM 2.5 cut cyclone (Very Sharp Cut Cyclone, VSCC, BGI Inc.) to avoid the impact of coarse-mode particles (e.g., dust), protected with a rain cap.The sample flow rate of the aethalometer was set to 5.0 liter per minute (L min −1 ).In this study, the measurements at 370, 520 and 880 nm were used for further analysis.The wavelength of 370 nm was used because studies have found that BrC shows strong light absorption in the near-UV wavelength range (Andreae and Gelencsér, 2006;Hansen and Schnell, 2005).Measurement data at the wavelength of 880 nm were chosen because light absorption at this wavelength normally represents the BC absorption (Virkkula et al., 2015).Regarding the wavelength dependence analysis of BC and BrC, the wavelength near 880 nm is better for the calculation of the AAE of BC because BC is the dominant absorption components in that wavelength range.However, the response of the 590 and 660 nm data may be affected by the presence of interfering materials such as hematite mineral dust and tobacco smoke (user manual of aethalometer AE-31, Hansen and Schnell, 2005); hence 520 nm data were used for the following calculations.Light scattering coefficients at three wavelengths (450, 525 and 635 nm) were measured by an integrating nephelometer (Aurora 3000, Ecotech).Sample air passed through a stainless tube with a rain cap and an external heater.In order to maintain sample air at low humidity (RH < 50 %), an internal heater was used.For those data with RH exceeding 50 % due to the malfunction of the heater, scattering coefficients were corrected for hygroscopic growth (Zhang et al., 2015).The light absorption coefficient (b abs ) at each wavelength λ was calculated using the method presented by Collaud Coen et al. (2010) to correct the systematic errors of filter-based absorption measurements.The attenuation coefficient b ATN at each wavelength λ is firstly calculated from where A and Q represent the spot size and flow rate, respectively.ATN(λ) is the attenuation change in a time step t.
The b abs at wavelength λ is then obtained after correction for a filter-loading effect, an embedded aerosol scattering effect and a multiple scattering effect by the filter fiber.The correction is performed using the Collaud Coen correction algorithm, with Schmid scattering correction adopted (Schmid et al., 2006;Collaud Coen et al., 2010).The equation can be presented as where R is the function for the filter-loading correction calculated using the equation from Collaud Coen et al. (2010, Eq. 13).C scat represents the aerosol scattering correction.
To calculate C scat , light scattering coefficients and scattering Ångström exponents measured by the nephelometer are used to obtain scattering at the aethalometer wavelengths, and the constants to calculate C scat are taken from Arnott et al. (2005).Detailed calculation equations of R and C scat can be found in Collaud Coen et al. (2010).C ref is the multiple scattering correction factor, which is set to be 4.26 according to Collaud Coen et al. (2010).A comparison of different aethalometer correction algorithms (Saturno et al., 2017) shows that the AAE derived by the Collaud Coen correction algorithm agrees well with that from multi-wavelength reference measurements, proving the reliable AAE values from this correction.The Collaud Coen correction also shows a good performance in obtaining absorption coefficients at 370 nm (Saturno et al., 2017), which is the critical wavelength in BrC segregation.Absorption coefficients are presented under standard temperature and pressure (STP, i.e., 273.15 K, 1013 hPa).Size distribution and mixing states of refractory BC were measured using the single particle soot photometer (SP2, Droplet Measurement Technologies, USA).The operation mechanism of SP2 has been described well in previous studies (Stephens et al., 2003;Schwarz et al., 2006Schwarz et al., , 2008)).SP2 uses a laser-induced incandescence technique equipped with a Nd:YAG laser (λ = 1064 nm) and optical detectors to quantify the size of a single particle by detecting scattering and a laser-induced incandescence signal.The mass of refractory BC can be determined by its nearly linear relationship with the peak height of the incandescence signal.The incandescence response of SP2 was calibrated using fullerene soot with a known mobility diameter selected by DMA (Gysel et al., 2011).The mass of fullerene soot was calculated using size-resolved effective density presented by Gysel et al. (2011).For pure scattering particles, the peak height of the scattering signal is linearly proportional to the particle scattering cross section.As for BC-containing particles, due to the loss of coating and the vaporization of the BC core when the particle passes through the laser beam, the scattering signal is different from the original particle.To determine the scattering cross section of BC-containing particles and saturated scattering particles, the leading-edgeonly (LEO) fit method developed by Gao et al. (2007) was adopted.D p of BC was then estimated by using the coreshell Mie model, assuming a BC core and shell refractive index of 2.26−1.26i(Moteki et al., 2010) and 1.52−0i (Pitchford et al., 2007), respectively.The scattering signal was calibrated using polystyrene latex (PSL) spheres with known sizes.The detection range of the BC core is 80-600 nm, assuming a density of 1.8 g cm −3 (Bond and Bergstrom, 2006).In this study, the four measurement periods are 20 May-12 June 2016, 8-31 August 2017, 1-30 November 2017 and 1-28 February 2018, representing spring, summer, autumn and winter, respectively.
In the present study, water-soluble ions (K + , Cl − ) were measured by the Monitor for Aerosols and Gases in Air (MARGA, Metrohm Co.); PM 2.5 and meteorological data were used for further supporting discussions.More detailed descriptions of these measurements can be found in Ding et al. (2013bDing et al. ( , c, 2016a)).

Optical calculation
It has been proved that BrC shows strong light absorbance in the UV-VIS wavelength range.To quantify the light absorption of BrC based on optical measurement results, the AAE segregation method is used (Lack and Langridge, 2013;Mohr et al., 2013).Light absorption of BrC is calculated as the result of b abs minus the light absorption coefficient of BC (b abs_BC ) at 370 nm.Here b abs_BC_370 is defined as the absorption coefficient of pure BC or BC with non-absorbing coating at 370 nm, and b abs_BrC is obtained from total absorption at 370 nm deducting absorption of the BC core and lensing effects, as the following equation shows: where b abs_370 and b abs_880 represent the absorption coefficients at 370 and 880 nm, respectively, which are calculated from the light absorption measurement data.AAE BC520-880 and AAE BC370-520 stand for the AAEs of pure BC and BC with non-absorbing coating at long and short wavelength ranges.We calculate AAE BC520-880 and AAE BC370-520 from SP2 data using the core-shell Mie model, which has been widely applied in BC-related studies (Bond and Bergstrom, 2006).It is mentioned that BC morphology can affect the AAE (Liu et al., 2015) and it is possible to overestimate BrC absorption; however, the complex morphology can vary with time and currently it is hard to evaluate its quantitative effect.Also, this site is a regional background station influenced more by aged air plumes (Ding et al., 2016a); therefore here we implement the core-shell model.The AAE at two wavelengths is calculated as the following equation: As mentioned above, an AAE value of 1.0 was adopted for BC by many researchers.However, AAE BC can vary with BC core size, coating thickness, morphology, etc.Evidence has shown that the AAE of pure BC cores can be lower than 1.0 as the diameter is out of the range of Rayleigh theory, and that BC with a clear shell might have an AAE higher than 1.0 (Bond et al., 2013;Lack and Cappa, 2010;Gyawali et al., 2009).It is also observed at the SORPES station that AAE 520-880 , which is expected to be mainly affected by BC absorption, is not always 1.0 and exhibits clear seasonal and diurnal variations (Shen et al., 2018).Hence, assuming an AAE BC of 1.0 in the estimation of BrC may induce large uncertainties or bias (comparison of calculated b abs_BrC assuming AAE BC = 1.0 versus the modified method will be discussed later).Therefore, it is essential to firstly evaluate the quantitative impacts of BC size and coating on the AAE value and determine the proper AAE BC for more accurate b abs_BrC calculation.

Estimation of BC optical properties and BrC segregation
Based on the core-shell Mie theory model, we conducted a series of calculations to discuss the variation pattern of the AAE for BC-containing particles (Bohren and Huffman, 1983).We used Christian Mätzler's code (Mätzler, 2002) for Mie calculations of spherical particles at different wavelengths.
Mie theory simulations were conducted firstly for monodispersed particles.Here particle core diameter (D c ) is defined as the diameter of the core alone, and the shell diameter (D p ) refers to the total particle diameter.Coating thickness was represented using D p /D c .D c increases from 1 to 200 nm with a 1 nm interval, and D p /D c was set to be 1.0-3.0 with 100 bins.The refractive index (RI) of the BC core was set to be 1.56 − 0.47i according to Dalzell and Sarofim (1969), and it was 1.52 − 0i for a clear shell (Pitchford et al., 2007).
Figure 1 shows the variations of AAE BC370-520 and AAE BC520-880 , with D c and D p /D c of monodispersed BCcontaining particles.The black dashed lines in the figure illustrate the D c and D p /D c value range, whereby BC is mostly distributed at the SORPES station (shown in Fig. 2).Within this range, it can be found that AAE BC370-520 and AAE BC520-880 mainly exhibit a deceasing pattern as D c increases but with a different amplitude.Taking pure BC as an example, AAE BC370-520 decreases from 1.1 to 0.6 as D c increases from 80 to 180 nm, while AAE BC520-880 decreases from 1.2 to 1.0 as D c increases.While as D p /D c grows, the AAE of BC shows a non-monotonous variation trend.For D c = 100 nm, as D p /D c grows from 1 to 3, AAE BC370-520 first increases from 1.0 to 1.3, peaks at D p /D c = 2 and then decreases back to 1.0 at D p /D c = 3, while AAE BC520-880 increases with D p /D c at Dc=100 nm.As for the magnitude of AAE, AAE BC370-520 is generally lower than AAE BC520-880 for BC with the same size and coating thickness.The above results suggest that the D c and D p /D c range of BC measured in this study are located in the regime in which AAE BC changes largely and non-monotonously.Moreover, in shortand long-wavelength ranges, AAE BC with the same BC size and coating thickness is also different.Therefore, instead of assuming the AAE to be a constant, real-time AAE BC determination is proposed in this study.
To explore the characteristics of AAE BC at short-and long-wavelength ranges at this observation site, we analyzed the size distribution and mixing state of BC measured by SP2 first.Figure 2 shows the overall D c number size distribution and D p /D c of BC-containing particles in the four seasons.It can be found that the number size distribution of BC cores in spring, summer and autumn shows a similar pattern.In winter, larger BC cores take up a higher proportion than other seasons.However, the coating thickness of BC is relatively lower in winter (peak number at D p /D c ∼ 1.6), possibly due to low photochemical oxidation in this season.The coating thickness of BC is higher in spring than other seasons.
Based on the size and coating thickness of each BCcontaining particle measured by SP2, b abs was calculated using the core-shell Mie model.AAE BC370-520 and AAE BC520-880 were then derived, as Fig. 3a illustrates.It can be observed that the fluctuation of AAE BC at both shortand long-wavelength ranges in different seasons is not significant.Median values of AAE BC520-880 in spring, summer, autumn and winter are 0.80, 0.78, 0.79 and 0.81, respectively, and AAE BC370-520 median values are 0.53, 0.54, 0.51 and 0.50 in the four seasons.Moreover, for BC at this site, AAE BC370-520 is always lower than AAE BC520-880 .Figure 3b shows the calculated AAE of particles measured by the aethalometer.The median values of AAE 520-880 in the four seasons range from 0.83 to 1.03, while AAE 370-520 shows a higher level and larger variation.The seasonal median values range from 1.00 to 1.31.Compared to Fig. 3a, it is noticed that observed AAE 520-880 is comparable to AAE BC520-880 , indicating that BC is the dominant absorbing  Dalzell and Sarofim (1969) and 1.52−0i for a clear shell (Pitchford et al., 2007).component in this wavelength range.Contrarily, the difference between AAE 370-520 and AAE BC370-520 is more obvious.Unlike AAE BC , AAE 370-520 is higher than AAE 520-880 in all seasons, indicating the possibility of the presence of other particle components with strong light absorption at a short wavelength range and a high AAE.
Since normally the SP2 is not used as a long-term continuous observation instrument, an alternative method to derive AAE BC is needed, especially for tracing back the historical level of BrC absorption without real-time SP2 measurement.As shown in Fig. 3, the difference between AAE BC370-520 and AAE BC520-880 at this site is not significant over time.Hence a correction factor R AAE is defined as the ratio of AAE BC370-520 to AAE BC520-880 calculated from SP2 data.Figure S2a (in the Supplement) illustrates the variation of R AAE .During the whole observation period, R AAE ranges between 0.60 and 0.69 (5th and 95th percentile), and the median value is 0.66, 0.69, 0.64 and 0.62 in spring, summer, autumn and winter, respectively.As mentioned above, AAE 520-880 calculated from aethalometer data is approximately equal to AAE BC520-880 in the wavelength range in which the main absorber is BC (Lack and Langridge, 2013).Therefore, AAE 520-880 is used to represent AAE BC520-880 ,    3) and (4).P BrC , defined as the 1.2 0.6 1.9 1.3 1.1 1.0 1.2 AAE 520-880 0.9 0.6 1.2 1.0 0.9 0.8 0.9 K + (µg m −3 ) 0.9 0.1 6.5 1.4 0.7 0.6 0.8 Cl − (µg m −3 ) 2.2 0.0 12.9 4.0 2.0 0.8 1.7  ), was also calculated using the measurement data.A statistical overview is summarized in Table 1.Seasonal cycles of b abs_BrC and P BrC are shown in Fig. 5, together with K + , K + / PM 2.5 , PM 2.5 and AAE at different wavelength ranges.Firstly, it can be found that b abs_BrC exhibits a distinct two-peak seasonal pattern, whereby the peak value occurs in June and December, with a mean b abs_BrC of 5.9 and 15.5 Mm −1 , respectively (Fig. 5a).It is also observed that b abs_BrC during winter, especially December, is much higher than that in the other three seasons (2 to 3 times higher).P BrC also presents a two-peak seasonal trend, with the high P BrC months of May-June and December.The mean P BrC in winter and summer is 19.6 and 14.4 %, respectively, which is lower than that in Xi'an but higher than the PRD region (Shen et al., 2017;Yuan et al., 2016).The 95th percentile of P BrC can reach to 32 % in December, which certainly cannot be ignored in the light absorption estimation in the YRD region.Notably, P BrC has a similar seasonal variation pattern with K + , except in February when intensive fireworks during the Chinese New Year can lead to significantly high values of K + concentrations.Moreover, AAE 370-520 shows distinct seasonal variations, which has a much wider range of change, with 1st and 99th percentiles of 0.6 and 1.9 (Table 1), than that of AAE 520-880 .It is also noticed that the variation pattern of AAE 370-520 is similar with K + .Since K + is mainly emitted from primary combustion processes, the simultaneous variation of P BrC and AAE 370-520 with K + suggests that primary emissions are likely to make a considerable contribution to BrC in this area for the most of time.
Due to the distinct seasonal variation in b abs_BrC observed at the SORPES station and its considerable contribution to the light absorption, it is essential to recognize the potential source areas and types of BrC in the YRD region.The diurnal variation of b abs_BrC was compared to b abs_BC in each season, shown in Fig. 6.Hourly mean values of b abs_BrC and b abs_BC are plotted.The reason to compare b abs_BrC with b abs_BC is that BC is one of the major light absorbers in the atmosphere and it is mostly from primary emission sources.Overall, b abs_BrC shows quite a similar diurnal pattern to b abs_BC in the four seasons, which is high during the night and starts decreasing after sunrise.This indicates that BrC at the site is also dominated by primary emissions.The lowest values occurred in the afternoon due to the development of the planetary boundary layer (PBL).It is observed that b abs_BC exhibits clear morning peaks in all four seasons, suggesting that traffic exhaust is likely to be a considerable emission source of BC at this site.Compared to that, the morning peaks of b abs_BrC in summer and autumn are less obvious than those of BC, while in winter and spring, the morning peaks are not noticeable.Such a difference reveals that during winter and spring, traffic emission is not the main contributor to local BrC.
In order to investigate the potential source region of BrC at the SORPES station, Lagrangian particle dispersion modeling (LPDM) was conducted following the method developed by Ding et al. (2013b) by using the Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model (Draxler and Hess, 1998;Stein et al., 2015).The meteorological data used in this model were GDAS (Global Data Assimilation System) data with a spatial resolution of 0.5 • in both latitude and longitude.In each simulation, 3000 particles were released at an altitude of 100 m above ground level (Wang et al., 2017a) and backwardly run for a 3-day period, and then the retroplume, i.e., the footprint of the surface at 100 m, was obtained following the method of Ding et al. (2013b).
As shown in Fig. 5a, the two distinct peaks of b abs_BrC are in June and December, respectively.It is then necessary to explore the possible emission sources in these 2 months.As mentioned before, since BrC is an operational definition, it is difficult to perform source apportionment for BrC.The feasible way to determine its sources is to compare the relationships between BrC and certain species that are possibly from the same emission sources.Since the emission of BrC is usually related to biomass burning (Laskin et al., 2015;Saleh et al., 2014), maps of fire counts in June and December of 2014 are presented in Fig. 7, together with the monthly averaged 3day backward retroplume, in order to first diagnose whether the majority of air plumes pass through open burning areas in these 2 months.Then, correlations between b abs_BrC and K +  in June and December were compared since K + is normally considered to be a tracer of primary emission from BB (Ding et al., 2013a).The results are shown in Figs.7 and 8.
Figure 7a shows intensive open burning, detected in June from the northwest of the site, but in this month air plumes are mainly transported from the eastern area of Nanjing where fire spots can also be found but are less concentrated than in northwestern regions.Contrarily, very few fire counts can be detected in December, suggesting much fewer open burning events in this month.Therefore, the high level of b abs_BrC in December is not likely to be from open BB emissions.The map of the retroplume reveals that air masses are mainly from the northern area in December.Then, correlations between daily average b abs_BrC and K + mass concentration in June and December are compared, as Fig. 8a shows.It can be observed that the correlations between b abs_BrC and K + in these 2 months exhibit clearly different patterns.The slope of fitted b abs_BrC and K + in June is 4.65 and the correlation coefficient R 2 is 0.92.Knowing the strong correlation between b abs_BrC and K + in June, combined with the observed intense fire counts in this month, it can be presumed that primary open BB emission could be a major contributor to BrC during June.As for December, b abs_BrC and K + present a much higher slope (slope is 10.59), which is approximately twice that in June.The distribution of b abs_BrC vs. K + in these two seasons is shown in Fig. 8b.The result shows that b abs_BrC / K + ratio of June and December shows a significant difference (through a t test, P < 0.05), indicating that the dominant emission source of BrC in December is not open biomass burning (a significant difference test of b abs_BrC / K + ratio is also done for comparison for May and June, which are the main open BB seasons, and the result shows that there is no significant difference in these 2 months; Fig. S5).The above analyses have demonstrated that BrC is not mainly from open biomass burning and vehicle exhausts in December.The increase of P BrC in winter, representing the higher ratio of BrC to BC mass, indicates the change of main emission sources.Nanjing is dominated by northeasterly wind in winter (Fig. S6), with air masses longrange transported from North China by the winter monsoon (Ding et al., 2013c(Ding et al., , 2016a)).Because of cold weather in winter, there is higher residential coal and biomass/biofuel burning emission for household heating (Fu et al., 2018).Therefore, residential coal and household biomass burning might be the main sources of high BrC in December under the influence of the winter monsoon (Zhang et al., 2016).Studies conducted in Beijing have also suggested the important contribution of residential biofuel (Yan et al., 2015;Cheng et al., 2016) and coal combustion (Yan et al., 2017) on BrC in northern China in winter.From the above analyses, it can be concluded that the high BrC absorption contributions in May-June and winter season are caused by different sources.In May-June, strong open biomass burning leads to the high BrC in Nanjing, while in winter, residential coal and household biofuel burning are possibly the major sources of BrC in winter season.Detailed sources of BrC can be further explored, combining field measurements of organic aerosols in the future.

Summary
In this study, light absorption of BrC was quantified using the optical method based on the definition of BrC.Mie theory simulation and observational results were combined to improve this method by calculating AAE BC at each time point instead of assuming a constant.Long-term variation of b abs_BrC and P BrC was then derived.Apparent light absorption contributed by BrC is discovered in the YRD region.b abs_BrC and P BrC both exhibit clear seasonal cycles with two peaks in May-June and December.The light absorption contribution of BrC at 370 nm ranges from 10.4 to 23.9 % (for 10th and 90th percentiles), and can reach 33.3 % in the open-BB-dominant and winter seasons.Comparison between b abs_BrC and b abs_BC suggests that vehicle emission has a negligible impact on the regional BrC level during winter and spring.Source analysis was performed based on temporal variations of BrC and the comparison of possible coemitted pollutants or related parameters.LPDM and MODIS fire data were also used to support analysis.The months of June and December with the peak level of b abs_BrC are chosen to analyze the potential emission sources of BrC, and it is found that the high contributions of BrC in these 2 months are dominated by different emission sources.In June, intensive primary open BB emission is the dominant source of BrC, making b abs_BrC have short-term high values, raising its average level in this month.In December, a high portion of BrC is possibly contributed by residential coal and household biofuel burning.
Overall, this work explores an improved optical method to quantify b abs_BrC from long-term observations.Compared to the conventional one, which may induce large uncertainty due to the assumption of constant AAE BC regardless of its variation with particle size, wavelength and time, this improved method is applicable to those sites at which the BrC proportion is low and makes long-term analysis available.This study also highlights the considerable contribution of BrC to light absorption in the near-UV range in the YRD region.Moreover, different emission sources of BrC are found in different seasons, providing a clearer reference for mitigation measures as well as regional control policies in eastern China.
Data availability.The GDAS data used in the LPDM calculation can be acquired from ftp://gdas-server.iarc.uaf.edu/gdas1/(last access: 15 June 2017).Measurement data at SORPES, including aerosol data and relevant trace gases as well as meteorological data, are available upon request from the corresponding author before the SORPES database is opened publicly.
Author contributions.Aijun Ding designed the study.Jiaping Wang contributed to the measurements, data analysis and the first draft of this paper.Wei Nie, Yicheng Shen, Xuguang Chi, Yuning Xie, Peng Sun, Zheng Xu and Ximeng Qi participated in observation at the SORPES station.Aijun Ding, Wei Nie, Yafang Cheng, Jiandong Wang, Xin Huang and Hang Su contributed to analysis, modeling and editing of the paper.
Competing interests.The authors declare that they have no conflict of interest.Special issue statement.This article is part of the special issue "Pan-Eurasian Experiment (PEEX)".It is not associated with a conference.

Figure 1 .
Figure 1.Variation of absorption Ångström exponents between (a) 370 and 520 nm (AAE BC370-520 ) and between (b) 520 and 880 nm (AAE BC520-880 ) along with core diameter (D c ) and the coating thickness (D p /D c ) of clear (pure scattering) shell simulated with the core-shell Mie model.The refractive index (RI) of the BC core is set to be 1.56 − 0.47i according toDalzell and Sarofim (1969) and 1.52−0i for a clear shell(Pitchford et al., 2007).

Figure 2 .
Figure 2. Normalized number size distributions of (a) the BC core (D c ) and (b) D p /D c of BC-containing particles from SP2 measurement in the four seasons.

Figure 4 .
Figure 4. Comparison between the time series of daily mean b abs_BrC using the modified method and using AAE BC = 1.0.The blue diamonds represent the calculation result using R AAE = 0.65, which is the mean value from SP2 data.Calculated b abs_BrC using AAE BC = 1.0 is plotted by the black line.

Figure 5 .
Figure 5. Seasonal cycle of (a) b abs_BrC , K + , and AAE at different wavelength ranges (AAE 370-520 , AAE 370-880 and AAE 370-880 , shown by the solid line, dashed line and diamonds, respectively) and (b) P BrC , K + / PM 2.5 and PM 2.5 .For b abs_BrC , K + , P BrC , K + / PM 2.5 and PM 2.5 figures, bold solid lines represent median values, diamonds show the monthly averages and thin solid lines forming the shaded area are 25th and 75th percentiles.
contribution of light absorption by BrC at 370 nm (P BrC = b abs_BrC b abs_370

Figure 6 .
Figure 6.Diurnal variations of b abs_BrC and b abs_BC in the four seasons of (a) winter, (b) spring, (c) summer, and (d) autumn, respectively.The plots show the occurrence frequencies of b abs_BrC in each b abs_BrC bin.The dark red and black circle lines represent the hourly mean b abs_BrC and b abs_BC , respectively.

Figure 7 .
Figure 7. Map of averaged 3-day backward retroplume and the fire counts for (a) June and (b) December in 2014.Fire count data are from the MODIS Collection 6 Active Fire Product provided by the NASA fire mapper (https://firms.modaps.eosdis.nasa.gov/map/,last access: 20 October 2017).

Figure 8 .
Figure 8.(a) Correlations between daily average b abs_BrC and K + mass concentration in June (black circles) and December (red triangles) and (b) boxplot of b abs_BrC / K + ratio in June and December (data are from the year 2014), where red lines represent the median value, blue boxes represent 25th and 75th percentile ranges and thin bars are 5th and 95th percentiles.

Table 1 .
Statistical summary of data measured at the SORPES station.