Articles | Volume 21, issue 11
Measurement report
07 Jun 2021
Measurement report |  | 07 Jun 2021

Measurement report: Strong light absorption induced by aged biomass burning black carbon over the southeastern Tibetan Plateau in pre-monsoon season

Tianyi Tan, Min Hu, Zhuofei Du, Gang Zhao, Dongjie Shang, Jing Zheng, Yanhong Qin, Mengren Li, Yusheng Wu, Limin Zeng, Song Guo, and Zhijun Wu

During the pre-monsoon season, biomass burning (BB) activities are intensive in southern Asia. Facilitated by westerly circulation, those BB plumes can be transported to the Tibetan Plateau (TP). Black carbon (BC), the main aerosol species in BB emissions, is an important climate warming agent, and its absorbing property strongly depends on its size distribution and mixing state. To elucidate the influence of those transported BB plumes on the TP, a field campaign was conducted on the southeast edge of the TP during the pre-monsoon season. It was found that the transported BB plumes substantially increased the number concentration of the atmospheric BC particles by a factor of 13 and greatly elevated the number fraction of thickly coated BC from 52 % up to 91 %. Those transported BC particles had slightly larger core size and much thicker coatings than the background BC particles. However, the coating mass was not evenly distributed on BC particles with different sizes. The smaller BC cores were found to have larger shell / core ratios than the larger cores. Besides, the transported BB plumes strongly affected the vertical variation in the BC's abundance and mixing state, resulting in a higher concentration, larger number fraction, and higher aging degree of BC particles in the upper atmosphere. Resulting from both increase in BC loading and aging degree, the transported BB plumes eventually enhanced the total light absorption by a factor of 15, of which 21 % was contributed by the BC aging, and 79 % was contributed from the increase in BC mass. Particularly, the light absorption enhancement induced by the aging process during long-range transport has far exceeded the background aerosol light absorption, which implicates a significant influence of BC aging on climate warming over the TP region.

1 Introduction

Biomass burning (BB) is an important source of atmospheric aerosols, exerting a significant impact on the regional and global climate system (von Schneidemesser et al., 2015; Jacobson, 2014). BB is widely occurring around the world as the means of heating, cooking, farming, etc. In some regions, BB has been becoming a serious environmental issue, which is harmful for both the public health and air quality, such as the raging wildfires in Australia and California, USA, this year. Due to the climate warming, this kind of natural BB event may increase in prevalence and severity in the future (Westerling et al., 2006; Abatzoglou and Williams, 2016). Black carbon (BC) is the main aerosol species in BB emissions. It can strongly absorb light across the whole solar spectrum. BC is an important global warming agent, with an estimated radiative forcing of +1.1 W m−2 (only secondary to CO2) (Bond et al., 2013). During the transport in the atmosphere, BC aerosols will experience a complex aging process and gradually become internally mixed with other aerosol species (non-BC coatings), which will lead to a great change in BC's properties (Jacobson, 2001). Those coatings on BC will enhance the light absorption of BC through the so-called lensing effect (Lack and Cappa, 2010), which will further enhance the refractive heating of BC (Jacobson, 2001; Peng et al., 2016). Quantifying the light absorption enhancement of BC requires the information of its size distribution, mixing state, morphology, and chemical compositions (Fierce et al., 2016, 2020).

The Tibetan Plateau (TP), known as the third pole and the roof of the world, is the highest plateau on the Earth and has the largest store of ice outside of the polar regions (Qiu, 2008; Ramanathan and Carmichael, 2008). The TP plays a vital role in the large-scale atmospheric circulation and the hydrological cycle of the entire Asian continent (Lau et al., 2010; Yao et al., 2019). The TP has long been considered to be one of the remote regions that are relatively less influenced by the pollutants emitted from human activities. However, growing evidence proves that anthropogenic emissions from the surrounding settlements have influenced this pristine place (Li et al., 2016; Cong et al., 2013). The warming rate of the TP has been found higher than the rate of global warming (Xu et al., 2009; Chen et al., 2015). Recent studies show that light-absorbing aerosols (e.g. black carbon and mineral dust) are partially responsible for the rapid warming and the accelerated glacier retreat over the TP (Qiu, 2008; Yao et al., 2019). The TP is surrounded by several global hotspots of BC emissions (Bond et al., 2013). Every year open biomass burning (BB) is prevailing in South and Southeast Asia during the pre-monsoon season, resulting in large amounts of BC emissions (Streets et al., 2003; Duncan et al., 2003; Z. Y. Cong et al., 2015; Pani et al., 2020). With the facilitation of the meteorological condition during this period, the significant amounts of BC can be transported into the TP and become a major source of BC in the Tibetan region (Bonasoni et al., 2010; Marinoni et al., 2010; Zhao et al., 2013). Previous studies reveal that southern Asia contributes more than half the amount of BC to the Tibetan region on an annual basis and becomes the dominant source during pre-monsoon seasons (Kopacz et al., 2011; Li et al., 2016; Lu et al., 2012; Zhang et al., 2015).

Nevertheless, existing research on BC aerosols over the TP is limited to the mass concentration measurement, and only a few studies pay attention to BC's mixing state. By using the electron microscopy techniques, researchers observed the morphology of single BC particles in the TP and found that most sampled BC particles were internally mixed with sulfate and organic matter (Cong et al., 2009; Li et al., 2015; Yuan et al., 2019; Zhang et al., 2001). In addition, there are two studies quantifying the number fraction of thickly coated BC in northeastern TP, which found a percentage of around 50 % (Wang et al., 2015, 2014). Therefore, more accurate information on the BC mixing state, especially the size-resolved mixing state, is still in need when estimating the BC's warming effects on the TP.

In this study, an intensive field campaign was conducted during the pre-monsoon season at a remote background site of the Tibetan Plateau. To characterize the mixing state of BC, a combined system of a differential mobility analyzer and single-particle soot photometer (DMA–SP2) was applied to obtain a robust size-resolved coating thickness of BC-containing particles (BC particles for short). By probing into the microphysical properties of BC particles, this study revealed that the recurring pre-monsoonal BC concentration peaks in the southeastern TP (Z. Cong et al., 2015; Zhao et al., 2013) could have resulted from regional transport of aged BC particles, which led to a strong light absorption over the region.

2 Methodology

2.1 Site description

The measurement station is located on Mt. Yulong (27.2 N, 100.2 E), with an altitude of 3410 m a.s.l. (Fig. 1). Mt. Yulong is situated at the southeastern edge of the Tibetan Plateau (TP), a transition zone that connects the vast low-altitude areas in southern Asia to the high-altitude TP. Different from the southern edge of the TP, the gentle slope of the southeastern TP provides a geographical convenience for the transport of air pollutants from southern Asia to the TP. The measurement period lasted from 22 March to 4 April 2015, corresponding to the pre-monsoon season in the region (Bonasoni et al., 2010; Z. Cong et al., 2015). The meteorological condition is shown in Fig. S1 in the Supplement. More detailed information about the measurement site can also be found in Zheng et al. (2017).

Figure 1Location of measurement site at Mt. Yulong. Elevation data from the Shuttle Radar Topography Mission (SRTM) 90 m digital elevation model (DEM) Version 4; base map from the Environmental Systems Research Institute (ESRI) World Ocean Base; manipulating by ArcMap Version 10.2, Redlands, CA.

2.2 Experimental setup

During the field measurement, BC's properties were analyzed by a combined DMA and SP2 system. After being dried by a silica-gel diffusion dryer, the ambient particles were introduced into a differential mobility analyzer (DMA;, Model 3081, TSI Inc., USA), in which particles were separated according to their electrical mobility (mainly determined by the particle size). Afterwards, size-resolved particles were driven into a single-particle soot photometer (SP2; Droplet Measurement Technologies, USA) and a paralleled condensation particle counter (CPC; Model 3022, TSI Inc., USA). In the SP2, the particles passed through an intracavity laser beam (Nd:YAG, λ=1064 nm) separately. During this process, the scattering signals of each particle were detected and recorded. If the particle contains black carbon (BC), it would be heated to its incandescence point, and the incandescence signals would also be recorded. The configurations of the DMA–SP2 coupled system were set as follows: the sheath flow rate of the DMA was 3.0 L min−1, and the aerosol flow rate was 0.45 L min−1; the CPC consumed 0.3 L min−1, and the remaining 0.15 L min−1 air was introduced into the SP2.

2.3 Data analysis

2.3.1 SP2 data

For non-absorbing particles, the peak intensity of the scattering signal is proportional to the scattering cross-section of the particle, which mainly depends on the particle size (Schwarz et al., 2006; Stephens et al., 2003). However, for absorbing particles (mainly BC-containing particles), the scattering signal will be lost due to the particle evaporation (Gao et al., 2007; Moteki and Kondo, 2008). To obtain the full scattering signal of these BC-containing particles, the LEO-fit method (Gao et al., 2007) was applied in this study, and the leading-edge signals were determined by using the method proposed by Liu et al. (2014) in their paper. Then, based on the calibration of the scattering signal, the scattering cross-section of both absorbing and non-absorbing particles could be determined. For BC-containing particles, the peak intensity of the incandescence signal is proportional to the BC mass in each particle, which is unbiased by the mixing state of BC. The calibration of the incandescence signal was also conducted by using the combined DMA–SP2 system, and Aquadag was used as the reference material. The formula provided by Gysel et al. (2011) was applied to convert the mobility diameter into mass for Aquadag. Considering the different sensitivity of the SP2 response to Aquadag and ambient BC, a correction of 0.75 was applied to derive the BC mass in the ambient particles (Moteki and Kondo, 2010; Laborde et al., 2012).

2.3.2 Mixing state of black carbon

BC's mixing state is referred to the state of how BC is mixed with other particle-phase species. As the aging process proceeds, the freshly emitted BC particles will be gradually coated by other chemical species and convert from bare (or thinly coated) BC into thickly coated BC. To quantify BC's mixing state, there are three parameters that were used in this study: the number fraction of thickly coated BC, the shell / core ratio, and the coating thickness. The separation of thinly and thickly coated BC is based on the time delay between the occurrence of the peaks of the scattering signal and the incandescence signal (Moteki and Kondo, 2007; Shiraiwa et al., 2007). Since this method does not require any assumptions or complex calculation, the parameter of the thickly coated fraction can be compared between different studies. The shell / core ratio denotes the diameter ratio of the whole particle diameter (Dp) to the BC core diameter (Dc), following Eq. (1). The coating thickness can be seen as the thickness of the shell, which equals half of the difference between Dp and Dc (Eq. 2).

(1)Shell/core ratio=DpDc(2)Coating thickness=Dp-Dc2

In previous SP2 studies, the Dp usually refers to the optically equivalent diameter of the whole particle diameter, which is derived from the Mie calculation with several presumed input parameters (Taylor et al., 2015). Meanwhile, the Dc refers to the volume-equivalent diameter (or mass-equivalent diameter) of the BC core by assuming a fixed density for a void-free spherical BC core (1.8 g cm−3 is widely used). However, the BC cores in ambient BC particles always contain some inside voids and tend to contain more voids in thinly coated BC particles. Thus, using a density of 1.8 g cm−3 will underestimate the core size and overestimate the coating thickness (and shell / core ratio) to some extent. In this study, with the benefit of the DMA–SP2 coupled system, the Dp can be directly measured by the DMA, which is the mobility diameter of the whole particle. To derive the Dc, instead of using the fixed density, we applied the closure study suggested by Zhang et al. (2018) to derive the size-dependently effective density for thinly coated BC particles (Fig. S2) and adopted a density of 1.2 g cm−3 for thickly coated BC particles (more details on the calculation method can be found in Zhang et al., 2018). The multiply charged particles induced by the DMA can be removed by comparing the mobility diameter with the optical diameter derived from the SP2 data. Figure S3 shows the comparison between the mobility diameter and optical diameter of thinly and thickly coated BC for single-charged particles. It should be noted here that, when it comes to the BC size distribution, the mass-equivalent diameter of BC cores (assuming a density of 1.8 g cm−3) was adopted in this study for direct comparison with previous studies.

2.3.3 Light absorption of black carbon

Light absorption of coated BC particles was calculated based on Mie theory (Bohren and Huffman, 1983). The refractive index of BC coatings is the volume average of different chemical species in BC's coating, with the assumption that the coatings share the same chemical composition with the whole submicron aerosols. The chemical composition was measured by a high-resolution time-of-flight aerosol mass spectrometer (AMS; Aerodyne Research Inc., Billerica, MA, USA). The density and refractive index of different chemical species are from Table 1 of Barnard et al. (2010) and references therein. For the simplification of optical calculation, all the BC particles were sized into different bins according to their core diameters, and all the BC particles in each bin were assumed to have the same shell / core ratio (median), which also means a united absorption cross-section (Cabs,I, calculated from the Mie theory) value in each bin. To calculate the bulk mass absorption cross-section (MAC) for the particle population, the absorption cross-section of BC particles in each bin (Cabs,i) was summed and divided by the sum of BC mass in each bin (Eq. 3).

(3) MAC = i = 1 n C abs , i × N i i = 1 n M bc , i × N i

In Eq. (3), n stands for the total number of size bins. Cabs,i and Mbc,i are the absorption cross-section and BC mass of a particle in the ith bin, respectively. Ni is the total particle number in the ith bin.

Figure 2(a) Coating thickness and (b) number concentration of BC particles. (c) Number fraction of thickly coated BC particles; 72 h back trajectories and active fire spots during (d) Period I and (f) Period II. The back trajectories were calculated by using the Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model (, last access: 19 May 2021). The active fire points were obtained from the Fire Information for Resource Management System (FIRMS), provided by the Moderate Resolution Imaging Spectroradiometer (MODIS) satellite (, last access: 19 May 2021).

3 Results

3.1 Observational evidence for the transport of black carbon aerosols

Figure S4 shows the time series of aerosol mass concentration during this campaign. A pollution episode is marked out for the first few days (23 March to 29 March, marked as Period I) because they are characterized by high concentrations of all aerosol species. After that, the aerosol mass concentration gradually decreases and remains at an extremely low level in the last few days (1 April to 4 April, marked as Period II). In the following discussion, Period II is regarded as the background condition in this region. Then, it can be noticed from Fig. 2a that the coating thickness of BC particles shows a distinct difference between Period I and II: during Period I, the coating thickness is mainly distributed within 30–60 nm, but this mode disappears in Period II. The coating thickness in Period II is distributed around zero. It should be mentioned that the coating thickness here is only calculated for the particles with BC core diameter between 100–120 nm (mass-equivalent diameter). Aside from the coating thickness, the number fraction of thickly coated BC is also higher in Period I than in Period II, with an average of 81.8 % in Period I compared with an average of 51.7 % in Period II (Fig. 1c). The discrepancy in the mixing state between the two periods adequately addresses that the BC particles observed during Period I are highly aged, while the BC particles observed during Period II are relatively fresh. Corresponding to this variation, the abundance of BC-containing particles also exhibits a large difference between Period I and II (Fig. 2b). The average BC number concentration for Period I is 165 cm−3, with a maximum value of 460 cm−3, whilst the average BC number concentration for Period II is only 13 cm−3.

To further explore the origin of those abundant aged BC particles, the air mass back trajectories were investigated. During Period I, the air mass is mainly from the southwest and passes through the areas covered by dense fire spots (Fig. 2d). Differently, the trajectories during Period II are mainly from the west to northwest, with much longer trajectory length and in the absence of fire spots along the pathway. Hence, there is a strong possibility that those abundant aged BC particles observed in Period I are originated from the fire smoke in South and Southeast Asia. This deduction is supported by the chemistry-related evidence provided by Zheng et al. (2017) and Wang et al. (2019). The concentrations of biomass burning tracers, including the acetonitrile, levoglucosan, and potassium ion, are much higher in Period I than those in Period II (Fig. S5). Hence, those abundant aged BC particles observed in Period I are believed to have a highly possible origin of biomass burning.

The high concentration of BC particles in the pre-monsoon season over the TP is not an accidental phenomenon discovered only by our observation. Zhao et al. (2013), Z. Cong et al. (2015), and Sarkar et al. (2015) also found that the BC (or elemental carbon, EC) concentration reaches its maximum of a year in the pre-monsoon season in the Tibetan region, and they also correlated this phenomenon to the transport of biomass burning (BB) aerosols in South Asia. However, our observation provides direct evidence that those BC particles are highly aged, indicating that those BC particles have experienced long-range transport to arrive on the TP. The occurrence of this kind of transport event in the pre-monsoon season can result from the prevalence of BB activities in South and Southeast Asia (Streets et al., 2003; Duncan et al., 2003) and the facilitation of the prevailing westerly wind during this period (Zhao et al., 2013; Cao et al., 2010). Since this kind of BB transport event commonly occurs in a specific season every year, the research on its impact on the TP area is of great importance.

Figure 3Diurnal variation in BC number concentration (a), BC number fraction (b), and shell / core ratio (c) during Period I and Period II. The solid dots are the mean values. The rectangles cover the data from the 25 % to 75 % quantile. The short lines in the middle of the rectangles represent the median values.


3.2 Changes in the diurnal cycle of BC abundance and mixing state under the impact of transported BB plumes (implication for vertical distribution)

Vertical transport of the air within the planet boundary layer (PBL) is an important source of the particles in the free troposphere. It can be reflected by the high particle concentration during the daytime due to the strong convective mixing (Lugauer et al., 2000; Nyeki et al., 1998). In a previous publication based on the same measurement, Shang et al. (2018) found an increasing trend of Aitken-mode particle number concentration during the daytime. By examining the available evidence, they believed that those Aitken-mode particles came from the anthropogenic particles which are uplifted by the elevated PBL. Here, our results show that the BC number concentration substantially increases (by almost a factor of 3) during the period from 08:00 to 12:00 (local time) in Period II (Fig. 3a2), which coincides with the time when the water vapor concentration increases (Fig. S6a2). This indicates that the observation site was probably influenced by the vertical transport during this time. Because the water vapor in the atmosphere is supposed to mostly come from the evaporation of the surface (Kivekäs et al., 2009). Here, the water vapor concentration is calculated based on the temperature and relative humidity. After 13:00, both the BC and water vapor concentration show a decreasing trend, which could have resulted from the scavenging effect of the high-speed wind (Fig. S6b2). However, under the influence of transported BB plumes, the daytime increase in BC number concentration becomes less pronounced Fig. 3a1, and the concentration during 01:00 to 06:00 becomes higher. If the data during 01:00 to 06:00 reflect the condition of the free troposphere as suggested by Shang et al. (2018), our results show that the regional transport of BB plumes considerably increases the BC aerosols in the free troposphere.

The diurnal cycle of BC relative abundance is also changed due to the transported BB plumes. In the background case of Period II (Fig. 3b2), the number fraction of BC-containing particles shows an increasing trend in the morning as the BC concentration does. A similar pattern was also found by Liu et al. (2020) at a mountain site located in the North China Plain, where they used the mass fraction instead of the number fraction. Liu et al. (2020) believed that the increase in the relative abundance of BC during the daytime was owing to more efficient vertical transport of BC than other particle species. Besides, they also suggested that higher BC fraction led to a lower single-scattering albedo, which resulted in an increase in positive radiative forcing. However, the diurnal pattern observed during Period II is almost reversed during Period I. The number fraction of BC particles is higher in the night than that in the day during Period I (Fig. 3b1). This indicates that the transport of BB plumes may lead to a higher BC fraction in the upper atmosphere. Therefore, when the vertical convection lifts the surface air to a higher altitude and mixes with the transported BB plumes, it will reduce the proportion of BC particles in the upper atmosphere.

Beside the absolute and relative abundance, the mixing state of BC is also essential to evaluate the radiative effects of BC. Here, we use the mean shell / core ratio of BC particles with 100–120 nm core diameter (mass-equivalent diameter) to indicate the overall aging degree. A noticeable difference in the diurnal pattern of BC mixing state is also found between Period I and II. For the background condition (Period II), the shell / core ratios fluctuate at a lower level and do not show an obvious pattern (Fig. 3c2). This implicates that, in the background air, there is no significant difference in the BC aging degree within and above the PBL. However, a clear drop in the shell / core ratios is observed during the daytime of Period I (Fig. 3c1). Based on the above discussion, this drop can be caused by the mixing with the less aged BC within the PBL.

Briefly, by comparing the different diurnal patterns during Period I and II, we find that the regional transport of BB plumes may change the vertical variation in BC abundance and mixing state in background condition, which may lead to a profound impact on the radiative balance in the region.

Figure 4(a) Number and (b) mass size distribution of BC cores. (c) Number size distribution of whole BC-containing particles. The left axis corresponds to the values in Period I, and the right axis corresponds to the values in Period II. The mass-equivalent diameters of BC cores are obtained by assuming a density of 1.8 g cm−3.


3.3 Variation in BC size distributions under the impact of transported BB plumes

The aging of BC particles is often caused by the secondary formation of coating materials or by the coagulation with co-existing particles, both resulting in an increase in the particle size. Meanwhile, the BC mass in each BC-containing particle is generally considered to remain unchanged during the aging process. Hence the size distribution of BC cores may give some hints of BC's origin. Figure 4a shows the number size distribution of BC cores in Period I and II. A slightly larger size is found in Period I. The fitted mean diameters for the two periods are 109 and 97 nm, respectively. For the mass distribution, the fitted mean diameters of each period are 181 and 177 nm, respectively (Fig. 4b). Additionally, the mean diameter of mass distribution observed in Period I is quite close to the values (181–182 nm) measured from an aircraft flying over northern India also in the pre-moon season (Brooks et al., 2019). Figure 4c shows the number size distribution of whole BC particles (including the BC core and non-BC coatings) during Period I and II. The fitted mean diameters of each period are 257 and 186 nm, respectively. Since there is only a minor increase on the BC core size from the background air (Period II) to the BB-plume-influenced air (Period I), the significant increase in the whole BC particle diameter mainly results from the atmospheric aging, which is consistent with the discussion above (Sect. 3.1).

Figure 5Variations in the number fraction of thickly coated BC and the shell / core ratio among the particles containing varying BC mass during Period I and Period II.


3.4 Size-resolved mixing state of black carbon

In previous sections, we demonstrate that BC particles of Period I have a higher aging degree than the background BC particles. However, the mixing state between single particles could be various, namely, different BC particles possibly have a variety of aging degree. As illustrated in Fig. 5, the percentage of thickly coated BC (light-blue bars) decreases as the BC core size increases, from 94 % at ∼100 nm core diameter to 30 % at ∼300 nm core diameter in Period I. The median shell / core ratios in each core size bin also show a significant declining trend in Period I, from 1.95 dropping to 1.00. During Period II, although both the thickly coated fractions and the shell / core ratios are much lower than those in Period I, they also show a decreasing trend as the BC core size increases.

Nevertheless, this kind of heterogeneity in mixing state of BC is always neglected in most models, which can lead to an overestimate of the BC absorption enhancement if a uniformed shell / core ratio is assumed. Fierce et al. (2016) used a particle-resolved model to examine the influence of the heterogeneity in mixing state on the absorption enhancement and found that the absorption is overestimated by as much as a factor of 2 if the diversity is neglected. In addition, their model-based analysis also revealed that the coatings tend to be less on particles that contain large BC amounts, which is in line with our discovery above. Thus, it is essential to consider the heterogeneity of mixing state when modeling the optical properties of black carbon.

Figure 6Comparison of BC mass, mass absorption cross-section, absorption enhancement, and absorption coefficient in Period I and II (at the wavelength of 550 nm).


3.5 Strong enhancement of the light absorption due to the transported BB plumes

Given that the transported BB plumes significantly increase the BC loading and BC aging degree in the atmosphere, they will exert a large impact on the light absorption level of the atmosphere. Here, this impact is accordingly considered from two different aspects: the impact on the abundance and the impact on the mixing state. As shown in Fig. 6, the BC mass concentration is elevated by a factor of 13 relative to the background air, indicating that the number of absorbing aerosols is substantially increased. Furthermore, due to the lensing effect, the thick coatings on those transported BC particles enhance the absorption ability of single BC particles compared with the thinly coated BC particles in the background air. Also as illustrated in Fig. 6, the mass absorption cross-section (MAC) is amplified by a factor of 1.2 relative to the background air, from 10.9 to 13.1 m2 g−1. If compared to the uncoated BC particles, the enhancement factor (Eabs) will be up to 1.25 during Period I, while the enhancement factor of background BC particles is only 1.03. It also should be noted that the MAC and Eabs values here reflect the absorption capacity of dry particles. If considering the uptake of water, the enhancement of light absorption will be more significant.

Combining these two aspects of the impact, the transported BB plumes eventually increase the total light absorption by a factor of nearly 15, from 0.6 Mm−1 in Period II to 8.8 Mm−1 in Period I (Fig. 6). In the total light absorption of Period I, 21 % is contributed from the lensing effect, which is related to the atmospheric aging during the long-range transport, and 79 % is contributed from the increase in BC abundance. Particularly, the light absorption enhancement contributed by the aging has far exceeded the background light absorption in Period II, emphasizing an important role of BC aging in the light absorption of the atmosphere in remote and clean areas. Meanwhile, the light absorption of the background air (Period II) is almost all contributed from the uncoated BC particles, and the contribution from the lensing effect is neglectable.

In our previous publication, Wang et al. (2019) have shown that the light absorption at short wavelengths by brown carbon is enhanced because of the transport of BB plumes during Period I. In addition to that research, this study reveals that the light absorption by black carbon, which can strongly absorb light at visible and infrared wavelengths, is also strongly enhanced by the transported BB plumes. Overall, it can be concluded that the transported BB plumes occurring during the pre-monsoon season strongly enhance the light absorption of the atmosphere across the whole solar spectrum, which is bound to have a profound radiative effect on the TP climate system.

4 Conclusions

Every year in the pre-monsoon season, the biomass burning plumes which originated in South and Southeast Asia can be transported to the Tibetan Plateau by the convenience of westerly wind. Based on field measurement, this study reveals the impact of the transported BB plumes on the abundance and mixing state of BC particles over the southeastern TP. As a result of the transport of BB plumes, the number concentration of BC-containing particles greatly increased, with the highest concentration exceeding the background air by more than a factor of 35. The percentage of thickly coated BC also increased significantly, from 52 % to as high as 91 %. Those transported BC particles were found to have much thicker coatings and slightly larger core size than the background BC particles. The coating thickness varied among BC particles of different core sizes, showing a trend of decreasing shell / core ratios with increasing core size. Due to the thicker coatings, the mass absorption cross-section of BC particles was enhanced by a factor of 1.2 compared to the background BC particles. Furthermore, the vertical variations in BC's abundance and mixing state substantially changed because of the transport of BB plumes, resulting in a higher concentration, larger number fraction, and higher aging degree of BC in the upper atmosphere compared to the background condition. The increase in both BC's abundance and aging degree eventually resulted in a 15-fold enhancement of the total light absorption on average. In the total light absorption, 79 % was contributed by the increase in BC's loading, and 21 % was contributed by the aging process during the long-range transport. Particularly, the light absorption enhancement induced by the BC aging has far exceeded the background light absorption. This study highlights the strong impact of the BC particles in the aged BB plumes on the TP climate system. The results in this study may be applicable to other regions beyond the TP, which is also influenced by aged BB plumes brought by long-range transport.

Code availability

The code used in the data processing is available from the authors upon request (

Code and data availability

The data presented in this article are available from the authors upon request (


The supplement related to this article is available online at:

Author contributions

TT wrote the manuscript and conducted the analysis of the datasets. MH, GZ, and ZW contributed important scientific thoughts and result analysis. ZD, DS, JZ, and YQ conducted the online measurement during the field campaign. ML, YW, LZ, and SG contributed to the field campaign.

Competing interests

The authors declare that they have no conflict of interest.


We sincerely thank the China National Environmental Monitoring Center and Lijiang Monitoring Station for support with the field campaign.

Financial support

This work was supported by the National Natural Science Foundation of China (grant no. 91844301) and the China Ministry of Environmental Protection Special Funds for Scientific Research on Public Welfare (grant no. 201309016).

Review statement

This paper was edited by Alex Huffman and reviewed by two anonymous referees.


Abatzoglou, J. T. and Williams, A. P.: Impact of anthropogenic climate change on wildfire across western US forests, P. Natl. Acad. Sci. USA, 113, 11770–11775,, 2016. 

Barnard, J. C., Fast, J. D., Paredes-Miranda, G., Arnott, W. P., and Laskin, A.: Technical Note: Evaluation of the WRF-Chem ”Aerosol Chemical to Aerosol Optical Properties” Module using data from the MILAGRO campaign, Atmos. Chem. Phys., 10, 7325–7340,, 2010. 

Bohren, C. F. and Huffman, D. R.: Absorption and Scattering of Light by Small Particles, John Wiley & Sons, Inc., New York, 1983. 

Bonasoni, P., Laj, P., Marinoni, A., Sprenger, M., Angelini, F., Arduini, J., Bonafè, U., Calzolari, F., Colombo, T., Decesari, S., Di Biagio, C., di Sarra, A. G., Evangelisti, F., Duchi, R., Facchini, MC., Fuzzi, S., Gobbi, G. P., Maione, M., Panday, A., Roccato, F., Sellegri, K., Venzac, H., Verza, GP., Villani, P., Vuillermoz, E., and Cristofanelli, P.: Atmospheric Brown Clouds in the Himalayas: first two years of continuous observations at the Nepal Climate Observatory-Pyramid (5079 m), Atmos. Chem. Phys., 10, 7515–7531,, 2010. 

Bond, T. C., Doherty, S. J., Fahey, D. W., Forster, P. M., Berntsen, T., DeAngelo, B. J., Flanner, M. G., Ghan, S., Karcher, B., Koch, D., Kinne, S., Kondo, Y., Quinn, P. K., Sarofim, M. C., Schultz, M. G., Schulz, M., Venkataraman, C., Zhang, H., Zhang, S., Bellouin, N., Guttikunda, S. K., Hopke, P. K., Jacobson, M. Z., Kaiser, J. W., Klimont, Z., Lohmann, U., Schwarz, J. P., Shindell, D., Storelvmo, T., Warren, S. G., and Zender, C. S.: Bounding the role of black carbon in the climate system: A scientific assessment, J. Geophys. Res.-Atmos., 118, 5380–5552,, 2013. 

Brooks, J., Liu, D., Allan, J. D., Williams, P. I., Haywood, J., Highwood, E. J., Kompalli, S. K., Babu, S. S., Satheesh, S. K., Turner, A. G., and Coe, H.: Black carbon physical and optical properties across northern India during pre-monsoon and monsoon seasons, Atmos. Chem. Phys., 19, 13079–13096,, 2019. 

Cao, J. J., Tie, X. X., Xu, B. Q., Zhao, Z. Z., Zhu, C. S., Li, G. H., and Liu, S. X.: Measuring and modeling black carbon (BC) contamination in the SE Tibetan Plateau, J. Atmos. Chem., 67, 45–60,, 2010. 

Chen, D. L., Xu, B. Q., Yao, T. D., Guo, Z. T., Cui, P., Chen, F. H., Zhang, R. H., Zhang, X. Z., Zhang, Y. L., Fan, J., Hou, Z. Q., and Zhang, T. H.: Assessment of past, present and future environmental changes on the Tibetan Plateau, Chinese Sci. Bull., 60, 3025–3035,, 2015. 

Cong, Z., Kang, S., Kawamura, K., Liu, B., Wan, X., Wang, Z., Gao, S., and Fu, P.: Carbonaceous aerosols on the south edge of the Tibetan Plateau: concentrations, seasonality and sources, Atmos. Chem. Phys., 15, 1573–1584,, 2015. 

Cong, Z. Y., Kang, S. C., Dong, S. P., and Zhang, Y. L.: Individual Particle Analysis of Atmospheric Aerosols at Nam Co, Tibetan Plateau, Aerosol Air Qual. Res., 9, 323–331, 2009. 

Cong, Z. Y., Kang, S. C., Gao, S. P., Zhang, Y. L., Li, Q., and Kawamura, K.: Historical Trends of Atmospheric Black Carbon on Tibetan Plateau As Reconstructed from a 150-Year Lake Sediment Record, Environ. Sci. Technol., 47, 2579–2586,, 2013. 

Cong, Z. Y., Kawamura, K., Kang, S. C., and Fu, P. Q.: Penetration of biomass-burning emissions from South Asia through the Himalayas: new insights from atmospheric organic acids, Sci. Rep.-UK, 5, 9580,, 2015. 

Duncan, B. N., Martin, R. V., Staudt, A. C., Yevich, R., and Logan, J. A.: Interannual and seasonal variability of biomass burning emissions constrained by satellite observations, J. Geophys. Res.-Atmos., 108, D2,, 2003. 

Fierce, L., Bond, T. C., Bauer, S. E., Mena, F., and Riemer, N.: Black carbon absorption at the global scale is affected by particle-scale diversity in composition, Nat. Commun., 7, 12361,, 2016. 

Fierce, L., Onasch, T. B., Cappa, C. D., Mazzoleni, C., China, S., Bhandari, J., Davidovits, P., Al Fischer, D., Helgestad, T., Lambe, A. T., Sedlacek, A. J., Smith, G. D., and Wolff, L.: Radiative absorption enhancements by black carbon controlled by particle-to-particle heterogeneity in composition, P. Natl. Acad. Sci. USA, 117, 5196–5203,, 2020. 

Gao, R. S., Schwarz, J. P., Kelly, K. K., Fahey, D. W., Watts, L. A., Thompson, T. L., Spackman, J. R., Slowik, J. G., Cross, E. S., Han, J. H., Davidovits, P., Onasch, T. B., and Worsnop, D. R.: A novel method for estimating light-scattering properties of soot aerosols using a modified single-particle soot photometer, Aerosol Sci. Tech., 41, 125–135,, 2007. 

Gysel, M., Laborde, M., Olfert, J. S., Subramanian, R., and Gröhn, A. J.: Effective density of Aquadag and fullerene soot black carbon reference materials used for SP2 calibration, Atmos. Meas. Tech., 4, 2851–2858,, 2011. 

Jacobson, M. Z.: Strong radiative heating due to the mixing state of black carbon in atmospheric aerosols, Nature, 409, 695–697,, 2001. 

Jacobson, M. Z.: Effects of biomass burning on climate, accounting for heat and moisture fluxes, black and brown carbon, and cloud absorption effects, J. Geophys. Res.-Atmos., 119, 8980–9002,, 2014. 

Kivekäs, N., Sun, J., Zhan, M., Kerminen, V.-M., Hyvärinen, A., Komppula, M., Viisanen, Y., Hong, N., Zhang, Y., Kulmala, M., Zhang, X.-C., Deli-Geer, and Lihavainen, H.: Long term particle size distribution measurements at Mount Waliguan, a high-altitude site in inland China, Atmos. Chem. Phys., 9, 5461–5474,, 2009. 

Kopacz, M., Mauzerall, D. L., Wang, J., Leibensperger, E. M., Henze, D. K., and Singh, K.: Origin and radiative forcing of black carbon transported to the Himalayas and Tibetan Plateau, Atmos. Chem. Phys., 11, 2837–2852,, 2011. 

Laborde, M., Mertes, P., Zieger, P., Dommen, J., Baltensperger, U., and Gysel, M.: Sensitivity of the Single Particle Soot Photometer to different black carbon types, Atmos. Meas. Tech., 5, 1031–1043,, 2012. 

Lack, D. A. and Cappa, C. D.: Impact of brown and clear carbon on light absorption enhancement, single scatter albedo and absorption wavelength dependence of black carbon, Atmos. Chem. Phys., 10, 4207–4220,, 2010. 

Lau, W. K. M., Kim, M. K., Kim, K. M., and Lee, W. S.: Enhanced surface warming and accelerated snow melt in the Himalayas and Tibetan Plateau induced by absorbing aerosols, Environ. Res. Lett., 5, 025204,, 2010. 

Li, C. L., Bosch, C., Kang, S. C., Andersson, A., Chen, P. F., Zhang, Q. G., Cong, Z. Y., Chen, B., Qin, D. H., and Gustafsson, O.: Sources of black carbon to the Himalayan-Tibetan Plateau glaciers, Nat. Commun., 7, 12574,, 2016. 

Li, W. J., Chen, S. R., Xu, Y. S., Guo, X. C., Sun, Y. L., Yang, X. Y., Wang, Z. F., Zhao, X. D., Chen, J. M., and Wang, W. X.: Mixing state and sources of submicron regional background aerosols in the northern Qinghai–Tibet Plateau and the influence of biomass burning, Atmos. Chem. Phys., 15, 13365–13376,, 2015. 

Liu, D., Allan, J. D., Young, D. E., Coe, H., Beddows, D., Fleming, Z. L., Flynn, M. J., Gallagher, M. W., Harrison, R. M., Lee, J., Prevot, A. S. H., Taylor, J. W., Yin, J., Williams, P. I., and Zotter, P.: Size distribution, mixing state and source apportionment of black carbon aerosol in London during wintertime, Atmos. Chem. Phys., 14, 10061–10084,, 2014. 

Liu, D. T., Hu, K., Zhao, D. L., Ding, S., Wu, Y. F., Zhou, C., Yu, C. J., Tian, P., Liu, Q., Bi, K., Wu, Y. Z., Hu, B., Ji, D. S., Kong, S. F., Ouyang, B., He, H., Huang, M. Y., and Ding, D. P.: Efficient Vertical Transport of Black Carbon in the Planetary Boundary Layer, Geophys. Res. Lett., 47, e2020GL088858,, 2020. 

Lu, Z. F., Streets, D. G., Zhang, Q., and Wang, S. W.: A novel back-trajectory analysis of the origin of black carbon transported to the Himalayas and Tibetan Plateau during 1996–2010, Geophys. Res. Lett., 39, L01809,, 2012. 

Lugauer, M., Baltensperger, U., Furger, M., Gaggeler, H. W., Jost, D. T., Nyeki, S., and Schwikowski, M.: Influences of vertical transport and scavenging on aerosol particle surface area and radon decay product concentrations at the Jungfraujoch (3454 m above sea level), J. Geophys. Res.-Atmos., 105, 19869–19879,, 2000. 

Marinoni, A., Cristofanelli, P., Laj, P., Duchi, R., Calzolari, F., Decesari, S., Sellegri, K., Vuillermoz, E., Verza, G. P., Villani, P., and Bonasoni, P.: Aerosol mass and black carbon concentrations, a two year record at NCO-P (5079 m, Southern Himalayas), Atmos. Chem. Phys., 10, 8551–8562,, 2010. 

Moteki, N. and Kondo, Y.: Effects of mixing state on black carbon measurements by laser-induced incandescence, Aerosol Sci. Tech., 41, 398–417,, 2007. 

Moteki, N. and Kondo, Y.: Method to measure time-dependent scattering cross sections of particles evaporating in a laser beam, J. Aerosol Sci., 39, 348–364,, 2008. 

Moteki, N. and Kondo, Y.: Dependence of Laser-Induced Incandescence on Physical Properties of Black Carbon Aerosols: Measurements and Theoretical Interpretation, Aerosol Sci. Tech., 44, 663–675,, 2010. 

Nyeki, S., Li, F., Weingartner, E., Streit, N., Colbeck, I., Gaggeler, H. W., and Baltensperger, U.: The background aerosol size distribution in the free troposphere: An analysis of the annual cycle at a high-alpine site, J. Geophys. Res.-Atmos., 103, 31749–31761,, 1998. 

Pani, S. K., Wang, S. H., Lin, N. H., Chantara, S., Lee, C. T., and Thepnuan, D.: Black carbon over an urban atmosphere in northern peninsular Southeast Asia: Characteristics, source apportionment, and associated health risks, Environ. Pollut., 259, 113871,, 2020. 

Peng, J. F., Hu, M., Guo, S., Du, Z. F., Zheng, J., Shang, D. J., Zamora, M. L., Zeng, L. M., Shao, M., Wu, Y. S., Zheng, J., Wang, Y., Glen, C. R., Collins, D. R., Molina, M. J., and Zhang, R. Y.: Markedly enhanced absorption and direct radiative forcing of black carbon under polluted urban environments, P. Natl. Acad. Sci. USA, 113, 4266–4271,, 2016. 

Qiu, J.: The third pole, Nature, 454, 393–396,, 2008. 

Ramanathan, V. and Carmichael, G.: Global and regional climate changes due to black carbon, Nat. Geosci., 1, 221–227,, 2008. 

Sarkar, C., Chatterjee, A., Singh, A. K., Ghosh, S. K., and Raha, S.: Characterization of Black Carbon Aerosols over Darjeeling – A High Altitude Himalayan Station in Eastern India, Aerosol Air Qual. Res., 15, 465–478,, 2015. 

Schwarz, J. P., Gao, R. S., Fahey, D. W., Thomson, D. S., Watts, L. A., Wilson, J. C., Reeves, J. M., Darbeheshti, M., Baumgardner, D. G., Kok, G. L., Chung, S. H., Schulz, M., Hendricks, J., Lauer, A., Karcher, B., Slowik, J. G., Rosenlof, K. H., Thompson, T. L., Langford, A. O., Loewenstein, M., and Aikin, K. C.: Single-particle measurements of midlatitude black carbon and light-scattering aerosols from the boundary layer to the lower stratosphere, J. Geophys. Res.-Atmos., 111, D16207,, 2006. 

Shang, D., Hu, M., Zheng, J., Qin, Y., Du, Z., Li, M., Fang, J., Peng, J., Wu, Y., Lu, S., and Guo, S.: Particle number size distribution and new particle formation under the influence of biomass burning at a high altitude background site at Mt. Yulong (3410 m), China, Atmos. Chem. Phys., 18, 15687–15703,, 2018. 

Shiraiwa, M., Kondo, Y., Moteki, N., Takegawa, N., Miyazaki, Y., and Blake, D. R.: Evolution of mixing state of black carbon in polluted air from Tokyo, Geophys. Res. Lett., 34, L16803,, 2007. 

Stephens, M., Turner, N., and Sandberg, J.: Particle identification by laser-induced incandescence in a solid-state laser cavity, Appl. Optics, 42, 3726–3736,, 2003. 

Streets, D. G., Yarber, K. F., Woo, J. H., and Carmichael, G. R.: Biomass burning in Asia: Annual and seasonal estimates and atmospheric emissions, Global Biogeochem. Cy., 17, 4,, 2003. 

Taylor, J. W., Allan, J. D., Liu, D., Flynn, M., Weber, R., Zhang, X., Lefer, B. L., Grossberg, N., Flynn, J., and Coe, H.: Assessment of the sensitivity of core / shell parameters derived using the single-particle soot photometer to density and refractive index, Atmos. Meas. Tech., 8, 1701–1718,, 2015. 

von Schneidemesser, E., Monks, P. S., Allan, J. D., Bruhwiler, L., Forster, P., Fowler, D., Lauer, A., Morgan, W. T., Paasonen, P., Righi, M., Sindelarova, K., and Sutton, M. A.: Chemistry and the Linkages between Air Quality and Climate Change, Chem. Rev., 115, 3856–3897,, 2015. 

Wang, Q. Y., Schwarz, J. P., Cao, J. J., Gao, R. S., Fahey, D. W., Hu, T. F., Huang, R. J., Han, Y. M., and Shen, Z. X.: Black carbon aerosol characterization in a remote area of Qinghai-Tibetan Plateau, western China, Sci. Total Environ., 479, 151–158,, 2014. 

Wang, Q. Y., Huang, R.-J., Cao, J. J., Tie, X. X., Ni, H. Y., Zhou, Y. Q., Han, Y. M., Hu, T. F., Zhu, C. S., Feng, T., Li, N., and Li, J. D.: Black carbon aerosol in winter northeastern Qinghai–Tibetan Plateau, China: the source, mixing state and optical property, Atmos. Chem. Phys., 15, 13059–13069,, 2015. 

Wang, Y. J., Hu, M., Lin, P., Tan, T. Y., Li, M. R., Xu, N., Zheng, J., Du, Z. F., Qin, Y. H., Wu, Y. S., Lu, S. H., Song, Y., Wu, Z. J., Guo, S., Zeng, L. W., Huang, X. F., and He, L. Y.: Enhancement in Particulate Organic Nitrogen and Light Absorption of Humic-Like Substances over Tibetan Plateau Due to Long-Range Transported Biomass Burning Emissions, Environ. Sci. Technol., 53, 14222–14232,, 2019. 

Westerling, A. L., Hidalgo, H. G., Cayan, D. R., and Swetnam, T. W.: Warming and earlier spring increase western US forest wildfire activity, Science, 313, 940–943,, 2006. 

Xu, B. Q., Cao, J. J., Hansen, J., Yao, T. D., Joswia, D. R., Wang, N. L., Wu, G. J., Wang, M., Zhao, H. B., Yang, W., Liu, X. Q., and He, J. Q.: Black soot and the survival of Tibetan glaciers, P. Natl. Acad. Sci. USA, 106, 22114–22118,, 2009. 

Yao, T. D., Xue, Y. K., Chen, D. L., Chen, F. H., Thompson, L., Cui, P., Koike, T., Lau, W. K. M., Lettenmaier, D., Mosbrugger, V., Zhang, R. H., Xu, B. Q., Dozier, J., Gillespie, T., Gu, Y., Kang, S. C., Piao, S. L., Sugimoto, S., Ueno, K., Wang, L., Wang, W. C., Zhang, F., Sheng, Y. W., Guo, W. D., Ailikun, Yang, X. X., Ma, Y. M., Shen, S. S. P., Su, Z. B., Chen, F., Liang, S. L., Liu, Y. M., Singh, V. P., Yang, K., Yang, D. Q., Zhao, X. Q., Qian, Y., Zhang, Y., and Li, Q.: Recent Third Pole's Rapid Warming Accompanies Cryospheric Melt and Water Cycle Intensification and Interactions between Monsoon and Environment: Multidisciplinary Approach with Observations, Modeling, and Analysis, B. Am. Meteorol. Soc., 100, 423–444,, 2019. 

Yuan, Q., Xu, J., Wang, Y., Zhang, X., Pang, Y., Liu, L., Bi, L., Kang, S., and Li, W.: Mixing State and Fractal Dimension of Soot Particles at a Remote Site in the Southeastern Tibetan Plateau, Environ. Sci. Technol., 8227–8234,, 2019. 

Zhang, D. Z., Iwasaka, Y., and Shi, G. Y.: Soot particles and their impacts on the mass cycle in the Tibetan atmosphere, Atmos. Environ., 35, 5883–5894,, 2001. 

Zhang, R., Wang, H., Qian, Y., Rasch, P. J., Easter, R. C., Ma, P.-L., Singh, B., Huang, J., and Fu, Q.: Quantifying sources, transport, deposition, and radiative forcing of black carbon over the Himalayas and Tibetan Plateau, Atmos. Chem. Phys., 15, 6205–6223,, 2015.  

Zhang, Y. X., Su, H., Ma, N., Li, G., Kecorius, S., Wang, Z. B., Hu, M., Zhu, T., He, K. B., Wiedensohler, A., Zhang, Q., and Cheng, Y. F.: Sizing of Ambient Particles From a Single-Particle Soot Photometer Measurement to Retrieve Mixing State of Black Carbon at a Regional Site of the North China Plain, J. Geophys. Res.-Atmos., 123, 12778–12795,, 2018. 

Zhao, Z. Z., Cao, J. J., Shen, Z. X., Xu, B. Q., Zhu, C. S., Chen, L. W. A., Su, X. L., Liu, S. X., Han, Y. M., Wang, G. H., and Ho, K. F.: Aerosol particles at a high-altitude site on the Southeast Tibetan Plateau, China: Implications for pollution transport from South Asia, J. Geophys. Res.-Atmos., 118, 11360–11375,, 2013. 

Zheng, J., Hu, M., Du, Z., Shang, D., Gong, Z., Qin, Y., Fang, J., Gu, F., Li, M., Peng, J., Li, J., Zhang, Y., Huang, X., He, L., Wu, Y., and Guo, S.: Influence of biomass burning from South Asia at a high-altitude mountain receptor site in China, Atmos. Chem. Phys., 17, 6853–6864,, 2017. 

Short summary
Every year in the pre-monsoon season, the black carbon (BC) aerosols originated from biomass burning in southern Asia are easily transported to the Tibetan Plateau (TP) by the convenience of westerly wind. This study reveals that the BC aerosols in the aged biomass burning plumes strongly enhance the total light absorption over the TP, and the aging process during the long-range transport will further strengthen the radiative heating of those BC aerosols.
Final-revised paper