HTML
--> --> --> -->2.1. Instrumentation
The characterization of atmospheric aerosols was performed using a sunphotometer (Cimel èlectronique, Paris, France). This instrument measures direct solar radiation every 15 min with a 1.2° full field-of-view at 1020, 870P1, 670, 440, 870P2, 870, 936, and 870P3 nm (P indicates a polarized filter). The full-width at half-maximum of each interference filter is 10 nm. Measurements at 1020, 870, 670, and 440 nm are used to retrieve AOD. The raw AOD is cloud-screened according to the method of (Smirnov et al., 2000), with an uncertainty of < 0.01 for wavelengths >440 nm (Eck et al., 1999).The Cimel sunphotometer also measures sky radiances under almucantar and principal plane scenarios at 440, 670, 870 and 1020 nm following the optical air mass protocol described by (Holben et al., 1998). Aerosol inversion algorithms and software (SKYRAD.pack version 4.2) are used to retrieve aerosol optical and microphysical properties from sky radiances. Other aerosol characteristics, such as SSA, the refractive index, and the particle volume size distribution, are obtained by this technique (Nakajima et al., 1996). Calibration of the sunphotometer was performed by CARSNET following the protocol developed by AERONET. Details about the intercomparison calibration and sphere calibration carried out at the Chinese Academy of Meteorological Sciences are given in (Che et al., 2009a) and (Tao et al., 2014).
Figure1. Location of Panyu——the observation site in the PRD region.
2
2.2. Methodology
2.2.1. Retrievals of AOD, ω, and AEThe retrieval of columnar aerosol optical properties from sky radiances requires an accurate correction for the effects of multiple scattering and for the contribution of light reflected from the earth's surface and scattered downward in the atmosphere (Valenzuela et al., 2012). The SKYRAD.pack version 4.2 software (Nakajima et al., 1996), which is commonly used to retrieve aerosol optical and radiative parameters (e.g., Kim et al., 2005; Che et al., 2008; Liu et al., 2008; Khatri et al., 2014; Wang et al., 2014), was used to calculate the AOD and SSA by using linear and nonlinear inversion schemes developed based on radiative transfer calculations. Additional input parameters for the SKYRAD.pack software include the Moderate Resolution Imaging Spectroradiometer (MODIS) Level 2 Collection 5 total ozone amount product (MOD07, 5 km × 5 km, daily) and the spectral surface reflectance product at seven wavelengths (MOD09, 0.5 km × 0.5 km, every eight days). The SSA is typically used to characterize aerosol absorption and is a key variable in assessing the radiative forcing due to aerosols. The ω describes the loss of photons to absorption, which is useful in identifying particle composition, especially carbonaceous particles (Corrigan et al., 2006). Here, ω is expressed as \begin{equation} \label{eq1} \omega(\lambda)=1-{\rm SSA}(\lambda) , \ \ (1)\end{equation} where Λ is the wavelength.
The AE is a good indicator of particle size and is given by the following equation: \begin{equation} \label{eq2} {\rm AE}=-\dfrac{\ln\left[\frac{{\rm AOD}(\lambda_1)}{{\rm AOD}(\lambda_2)}\right]} {\ln\left(\frac{\lambda_1}{\lambda_2}\right)} , \ \ (2)\end{equation} where Λ1 and Λ2 are the two wavelengths chosen to calculate the AE. The wavelength pairing of 440 nm and 870 nm was used in this study. The AE can range from negative values to greater than 1 (Gobbi et al., 2007), depending on the particle size.
2.2.2. Radiative forcing and radiative forcing efficiency
Radiative fluxes and aerosol direct radiative effects of total, diffuse, and direct shortwave irradiances in the broadband spectral range (0.25-4.00 μm) were calculated using the Santa Barbara Discrete Ordinate Radiative Transfer Model (DISORT) Atmospheric Radiative Transfer (SBDART) model (Ricchiazzi et al., 1998). The SBDART model is based on low-resolution band models developed for LOWTRAN 7 atmospheric transmission and the DISORT radiative transfer model, which has 33 layers and four radiation streams. Various studies have shown that the SBDART model can successfully simulate downwelling broadband fluxes at the surface and upwelling fluxes at the top of the atmosphere (TOA; Xia et al., 2007a; Li et al., 2010). SBDART simulations and measurements of broadband irradiance agree to within 3% (Halthore et al., 2005). This model has been used to estimate aerosol radiation forcing in China (Liu et al., 2007; Xia et al., 2007c, 2007d; Li et al., 2010).
The AOD, SSA, AE, and Asymmetry factor (ASY) at four AERONET wavelengths (i.e., 440, 675, 870, and 1020 nm) have been used to interpolate and extrapolate into the spectral divisions of the SBDART model (Xia et al., 2007a, 2007b). Values for the ASY (ASY440\; nm = 0.73; ASY675\; nm = 0.66, ASY870\; nm = 0.63; ASY1020\; nm = 0.61) from the AERONET site in Hong Kong (about 110 km to the southeast of Panyu) were used in this study and assumed to be constant. Vertical profiles of water vapor and ozone were obtained by partitioning total column water vapor amounts and total ozone amounts from the MOD07 product according to a standard model atmosphere (McClatchey et al., 1972). The surface albedo from the MODIS Albedo product (MCD43B3, https://lpdaac.usgs.gov/dataset_discovery/ modis/modis_products_table/mcd43b3) was derived using the 16-day anisotropy model provided in the MODIS BRDF/Albedo Model Parameters product (MCD43B1) at a resolution of 500 m. The total shortwave broadband albedo was used according to algorithms developed by (Liang, 2001) and set to be constant within the retrieval period.
ADRF denotes the direct effect of aerosols on the atmospheric energy budget. The SBDART model was run twice to simulate shortwave irradiances with and without aerosol particles under cloud-free conditions, and then used to determine the ADRF at the surface (F SFC) and at the TOA (F TOA). The aerosol radiative forcing within the atmosphere (F ATM) was defined as the difference between F TOA and F SFC. Instantaneous values of ADRF at different levels were calculated using the following equations: \begin{eqnarray} \label{eq3} \Delta F&=&F\downarrow -F\uparrow ,\ \ (3)\\ \label{eq4} {\rm F}_{\rm TOA}&=&\Delta F_{\rm TOA,\ with_aerosol}-\Delta {\rm F}_{\rm TOA,\ without_aerosol} ,\ \ (4)\\ \label{eq5} {\rm F}_{\rm SFC}&=&\Delta {\rm F}_{\rm SFC,\ with_aerosol}-\Delta{\rm F}_{\rm SFC,\ without_aerosol} ,\ \ (5)\quad\\ \label{eq6} {\rm F}_{\rm ATM}&=&{\rm F}_{\rm TOA}-{\rm F}_{\rm SFC} ,\ \ (6) \end{eqnarray} where ? F denotes the net downward flux (downward radiation, F↓, minus upward radiation, F↑).
The diurnal mean radiative forcing is given as \begin{equation} \label{eq7} {\rm d}F=\frac{1}{24}\smallint F(t){\rm d}t , \ \ (7)\end{equation} where t is the definite integral of time over the whole day or daytime only, and F(t) represents instantaneous radiative forcing values. On a given day, the number of data points may be limited because of the presence of clouds. To alleviate this limitation, monthly mean values from all instantaneous measurements were determined, from which monthly and annual means of shortwave ADRF were computed (Kim et al., 2005; Li et al., 2010). The uncertainties in the main input parameters of AOD, AE, SSA, ASY, surface reflectance, and ozone amounts caused an error of 8.76 3.44 W m-2 (Li et al., 2010).
We also computed the FE, defined as the change in diurnal mean radiative forcing (F) with respect to the change in AOD at 550 nm: \begin{equation} \label{eq8} {\rm FE}=\frac{F}{{\rm AOD}(\lambda_{550})} .\ \ (8) \end{equation} Λ550 was the wavelength at 550 nm. FE was calculated at the surface (FE SFC), at the TOA (FE TOA), and within the atmosphere (FE ATM) using Eqs. (3)-(8) (Li et al., 2010; Liu et al., 2011). The units are W m-2 for ADRF and W m-2 per AOD for FE. The AOD was calculated by the following equation (?ngstr?m, 1964): \begin{equation} \label{eq9} {\rm AOD}(\lambda)=\beta\lambda^{-{\rm AE}} , \ \ (9)\end{equation} where Λ is the wavelength (here, 550 nm) and β is the turbidity coefficient (AOD at Λ=1 μm).
2.2.3. Data
Data derived from the sky radiation algorithm were used to analyze the AOD, SSA and AE at a temporal resolution of 30 min. The temporal resolution of ADRF was the same as that of the aerosol property parameters. There were 1219 measurements acquired by the sunphotometer from January 2006 to December 2012. However, observations were not continuous because of instrument calibration and maintenance (Zhu et al., 2014). In addition, the PRD region is typically cloudy and rainy during the wet seasons (spring and summer), which greatly influences instrument observations and data inversions. For this reason, this study focused on the dry season (October through February).
-->
3.1. Aerosol optical properties
Figure 2a shows the frequency distribution of AOD at 550 nm for each month. The number of samples in each month (October, November, December, January, and February) was 109, 273, 373, 357 and 107, respectively. In general, the pattern in the frequency of AOD followed a Gaussian distribution, which is in accordance with (Behnert et al., 2007) but different from the log-normal distribution reported by (Liu et al., 2008). The AOD ranged from ~0.1-1.5, with more instances of high AOD in October and February than in other months. The frequencies of occurrence of AOD <0.3 were 7.4% (October), 9.5% (November), 23.9% (December), 5.3% (January), and 0.00% (February). The frequencies of occurrence of AOD >0.7 were 40.74% (October), 20.9% (November), 19.0% (December), 13.2% (January), and 56.1% (February). The frequencies of occurrence of AOD between 0.3 and 0.7 were 51.9%, 69.6%, 57.1%, 81.5% and 43.9%, respectively. This indicates that high AOD prevails during the dry season of the PRD region. About 11.7% of all AOD values were less than 0.2, with the greatest proportion of these values occurring in February. A greater proportion of AOD values less than 0.2 (22%) has been reported at Xianghe, where frequent airflow outbreaks lead to a relatively higher occurrence of background aerosol loading levels (Xia et al., 2005, 2007c). The low AOD situation in the PRD region most likely arises from the wet scavenging of aerosols after precipitation. However, these events are typically short-lived because of the high RH conditions prevalent in the region, in addition to the presence of local aerosol sources (Wang et al., 2014).Figure 2b shows the monthly variability in AOD at different wavelengths. The first and second highest monthly mean values of AOD occurred in February and October (0.9 and 0.7, respectively). The minimum monthly mean value of AOD (0.5) occurred in December. Note that the monthly mean AODs from November to January were not significantly different (p>0.05). The magnitudes of AODs from October to February were similar to those measured in Hefei, where the AOD is influenced by temperature, RH via hygroscopic growth, and stagnant weather conditions (Wang et al., 2014). The air temperature is of significance for the vertical convection of the aerosol loading. Considering that the temperature in February over the PRD region is relatively low compared to the other months, the maximum AOD value occurring in February is most likely due to the presence of frequent stagnant weather systems and the swelling effect caused by high RH.
Figure2. (a) Frequency distributions of AOD at 500 nm and (b) monthly variations in AOD at different wavelengths in the PRD region.
The SSA mainly depends on the chemical composition and size distribution of aerosol particles (Wang et al., 2014). Figure 3a shows the frequency of occurrence of SSA at 440 nm. The histograms of SSA have a Gaussian distribution. The frequencies of occurrence in October and February were higher than in other months. The frequencies of occurrence of SSA values <0.8, >0.92, and from 0.8 to 0.92, were 7.6%, 84.8%, and 14.3%, respectively. This suggests that the aerosols were strongly absorbing.
Relatively high values of SSA were seen in October (0.90 0.04), followed by moderate values in February (0.88 0.05) and lower values in the other months (ranging from 0.86 0.05 to 0.87 0.053) (Fig. 3b). The mean value of SSA was 0.87, which is lower than that reported at Xianghe (39.753°N, 116.961°E), a site located between two megacities (Beijing 70 km to the northwest and Tianjin 70 km to the southeast; Li et al., 2007) and at Xinglong (40.396°N, 117.578°E), which is located at the top of a mountain (970 m above sea level) ~ 100 km away from Beijing to the northwest (Zhu et al., 2014). Local pollutants and the products of their photochemical reactions, which are strongly absorbing, are the primary sources of aerosols in the PRD region. The other sites mentioned are more susceptible to coarse-mode mineral dust and local emissions. The relatively higher values of SSA in October and in February are associated with a stronger swelling effect due to the greater occurrence of a humid atmosphere (Fig. 4).
Figure3. (a) Frequency distributions of SSA at 440 nm and (b) monthly variations in SSA at 440 nm in the PRD region. The notches in the notched box and whisker plots represent the 95% confidence levels of the median values. Non-overlapping notches indicate that the median values are significantly different from each other.
The AE is an indicator of particle size. High values of AE mean that fine particles dominate, while low values mean that coarse particles dominate. The frequency distribution of AE (440-870 nm) in the PRD region is shown in Fig. 5a. The histograms of AE (440-870 nm) follow a Gaussian distribution, with a central value equal to ~1.3. Other studies have shown that if AE (440-870 nm) <0.75, coarse-mode dust particles dominate, and if AE (470-870 nm) >1.0, mainly fine-mode anthropogenic particles are present. Values of AE (470-870 nm) between 0.75 and 1.0 suggest complicated aerosol modes (Eck et al., 2005; Yang et al., 2011; Wang et al., 2014). The frequency of occurrence of AE (440-870 nm) <0.75 ranged from 0 in November and February to 1.83% in October. The proportion of AE (440-870 nm) between 0.75 and 1.2 varied from 6.54% in February to 17.43% in December. The majority of aerosols had an AE (440-870 nm) >1.2, ranging from 80.97% in December to 93.46% in February. On the whole, during the dry season, aerosols in the PRD region are mainly fine-mode pollution particles (~ 84.66%) and floating dust aerosols (~ 14.44%). The presence of heavy dust is negligible (~ 0.90%).
Figure 5b shows notched box whisker plots of the monthly AE (440-870 nm). The mean value of AE (470-870 nm) was 1.33. AE (440-870 nm) values varied from 1.31 to 1.38 and were significantly higher in February and October (p<0.01) and lower in December (p<0.01), suggesting a greater contribution of fine particles to extinction over the PRD region. The PRD region typically experiences high-humidity weather in February, when the background of cold air suddenly becomes warm air during the transition between winter and spring (Zhang et al., 2014). Frequent air mass stagnation episodes in February limit the diffusion of air, which results in higher AE (440-870 nm) due to anthropogenic aerosol pollution.
Figure4. Frequency distributions of RH. The frequency in each month was calculated using 10 RH bins ranging from 5% to 95%.
Figure5. (a) Frequency distributions of AE (440-870 nm) and (b) monthly variations in AE (440-870 nm) in the PRD region. The notches in the notched box and whisker plots represent the 95% confidence levels of the median values. Non-overlapping notches indicate that the median values are significantly different from each other.
2
3.2. Aerosol characterization inferred from the AE (440-870 nm), ω, and volume size distributions
The wavelength dependence of ω provides more information about the physicochemical properties of aerosol types and the AE provides information about particle size (Logan et al., 2013). Therefore, these two parameters were combined to classify aerosol types in the PRD region.As done by (Logan et al., 2013) when examining data from sites in China and Thailand, ω=0.07 was chosen as the demarcation line between strongly and weakly absorbing aerosols, and AE=0.75 was chosen as the demarcation line between fine-mode and coarse-mode particles. These threshold values were deemed reasonable to use in this study because the PRD region is also primarily affected by pollution, biomass-burning aerosols, mineral dust, and complex mixtures of various particle types (Lai et al., 2007; Lee et al., 2007; Deng et al., 2008). Four regions were defined: Region I, where fine-mode, weakly absorbing particles dominate; Region II, where fine-mode, strongly absorbing particles dominate; Region III, where coarse-mode, strongly absorbing mineral dust particles dominate; and Region IV, where coarse-mode, weakly absorbing particles dominate, e.g., desert aerosols.
Figure 6 shows the mean classification results for the months of October to February of the years 2006-12. About 9.5% of all data points fell in Region I, i.e., the region representing fine-mode, weakly absorbing particles. Because few dust events take place in the PRD region, the number of data points in the coarse-mode regions (III and IV) was negligible. Up to 90% of all particles fell in region II, mainly in January (Fig. 6d). The majority of region II data points were centered on ω values between 0.07 and 0.17, and AE values between 1.2 and 1.5.
Figure6. Classification of the mean physicochemical properties of aerosols into four regions using AE (440-870 nm) and ω at 440 nm. The color scale represents the relative density of points. Panels (a-f) correspond to October, November, December, January, February, and all data, respectively, in the years 2006-12. R is the proportion of particle amount.
Overall, fine-mode, strongly absorbing aerosol particles were dominant in the study area. Given the location of Panyu, pollutants at the site are likely a combination of carbon aerosols [e.g., organic and black carbon (BC)] or mixtures of sulfate, nitrite, and carbon aerosols generated by vehicles and stoves (Andreae et al., 2008; Zheng et al., 2011; Wang et al., 2012). Compared to region II, aerosols in region I had a higher AOD (0.65) and larger particle size ( AE=1.30).
The average volume size distributions for different months are shown in Fig. 7. The column-integrated aerosol volume size distribution is bimodal, with a dominating fine mode (radius <1.16 μm) and a coarse mode (radius >1.16 μm). The peak in the fine mode was the same for all months (radius = 0.12 μm). Peak volumes of the fine mode were higher in February and in October, which had almost equal values (0.11 and 0.10 μm3 μm-2, respectively). Peak volumes were lower in the other months and ranged from 0.07-0.08 μm3 μm-2. The coarse mode showed maxima in radii equal to 5.29, 5.29, 11.31, 11.31 and 3.62 μm in October, November, December, January and February, with peak volumes of 0.033, 0.031, 0.038, 0.037 and 0.028 μm3 μm-2, respectively. Because dust events in the PRD region are scarce, the monthly coarse modes were similar, with a mean peak volume 0.43 times that of the peak volume in the fine mode. The lower tail of the coarse mode in February is likely the result of fine-mode contamination.
2
3.3. Diurnal radiative forcing and efficiency
In the PRD region, high aerosol loading prevails all year round, with an average AOD equal to 0.54. As a result, aerosols are expected to seriously impact irradiances reaching the earth's surface, which can significantly influence atmospheric stability and regional climate.Previous studies have revealed that ADRF varies significantly in different regions and at different time scales (Xia et al., 2007a, 2016; Li et al., 2010). Other studies have reported only instantaneous values of aerosol radiative forcing (e.g., Che et al., 2009b, 2014), which should not be compared with daily and annual averages. Annually averaged global values of F SFC, F TOA and F ATM from satellite-model integrated approaches over land are -11.9, -4.9 and 7.0 W\;m-2, respectively, and from model simulations are -7.6, -3.0 and 4.6 W m-2, respectively (Yu et al., 2006). The monthly mean diurnal ADRF in the PRD region is shown in Fig. 8. The largest monthly values were in February, when F SFC and F ATM reached -37.5 and 30.0 W m-2, respectively. This happened because AOD values were high and SSA values were low in February. The average F SFC, F ATM and F TOA was -33.4 7.0, 26.1 5.6 and -7.3 2.7 W m-2, respectively. The large negative F SFC and positive F ATM suggests that more solar radiation is absorbed within the atmosphere, consequently warming the atmosphere, reducing eddy heat convergence, and inducing a reduction in surface temperature, which ultimately affects regional climate (Liu et al., 2008; Ge et al., 2010). Figure 9 shows the diurnal ADRF at the surface, within the atmosphere, and at the TOA, as a function of AOD. The ADRF was almost linearly correlated with AOD. The standard deviations at the surface (3.5 W m-2), at the TOA (1.96 W m-2) and in the atmosphere (3.81 W m-2) were much lower than the values reported by (Li et al., 2010) (14.45, 5.45 and 16.08 W m-2, respectively). This is because the study region here is much smaller and there is little variation in aerosol composition, which is also suggested by the near constant magnitude of the SSA (Fig. 3b).
Figure7. Average volume size distributions of aerosols as a function of particle radius over the PRD region for different months. Error bars represent standard deviations.
Figure8. Monthly mean diurnal ADRF at the surface (F SFC), at the top of the atmosphere (F TOA), and within the atmosphere (F ATM). Error bars represent standard deviations.
Figure 10 shows the interannual trends in diurnal ADRF during the dry season in the PRD region. The relative standard deviations of F SFC, F TOA and F ATM were ~ 20.96%, 36.99% and 21.46%, respectively. These values were less than those reported over China in general [56.96% to 533.33% (Li et al., 2010)] and in the Yangtze Delta region [26.15% to 58.54% (Liu et al., 2012)]. The lower values of the relative standard deviation are mostly due to relatively less variation in aerosol composition and partly due to month-to-month changes in radiative forcing. The magnitudes of F SFC varied from -36.2 W m-2 in 2006 to -32.3 W m-2 in 2012, and the magnitudes of F TOA varied from -7.2 W\;m-2 in 2006 to -10.6 W m-2 in 2012. Both of these trends did not pass the 99% confidence level of the t-test (p<0.01; r SFC=0.53, r TOA=0.24; n=6). By contrast, F ATM varied from -29.1 W m-2 in 2006 to -21.7 W m-2 in 2012, and this trend was statistically significant (p<0.01; r ATM=0.92; n=6). Moreover, the SSA increased from 0.87 in 2006 to 0.91 in 2012, and the fine-mode aerosol concentration such as PM2.5 decreased from 69.5 \(\mu\rm g\;m^-3\) in 2006 to 45.0 \(\mu\rm g\;m^-3\) in 2012. These results suggest a trend of less-absorbing particles being released into the atmosphere during the study period. This phenomenon is most likely due to the pollutant-controlling measures executed by the local government. In recent years, the local government has implemented a series of control measures with respect to emissions from coal burning, motor vehicle exhaust, and industrial sectors, which has resulted in an apparent decrease in the number of hazy days (http://www.gzepb.gov.cn/, in Chinese) and the amount of PM2.5 (Wang et al., 2016). Owing to these measures, concentrations of carbonaceous aerosols, such as organic carbon and elemental carbon, have shown a decreasing trend (Fu et al., 2014; Wang et al., 2016). The weaker radiative forcing in 2010 may be attributable to the stricter emission reduction measures in place during the 2010 Asia Games held in Guanzhou, resulting in AOD reducing to 0.43, and PM2.5 and elemental carbon concentrations decreasing by 29% and 49%, respectively (Hu et al., 2013).
Figure9. Diurnal ADRF at the surface (SFC, black), at the top of the atmosphere (TOA, blue), and within the atmosphere (ATM, red) as a function of AOD at 550 nm over the PRD region. Best-fit lines through the points are shown and linear regression functions and statistics (coefficient of determination, R2, and standard deviation, SD) are given.
Figure10. Interannual trends in diurnal ADRF during the dry season (October to February) in the PRD region. Dashed red, green, and blue lines represent linear-fit lines for ADRF within the atmosphere (ATM), at the top of the atmosphere (TOA), and at the surface (SFC), respectively. Vertical bars represent one standard deviation computed from the monthly mean forcing.
Table 1 summarizes the ADRF over the PRD and other regions. The presence of anthropogenic aerosols over the PRD region can considerably decrease surface and TOA forcings and enhance atmospheric forcing. The magnitudes of the mean F SFC and F ATM in this study region are about two times greater than that reported over the whole of China, but are similar to those reported in the Yangtze Delta region and in northern China where anthropogenic and dust aerosols dominate. Relatively high RH over the PRD region can partly explain the larger aerosol cooling effect at the surface compared to that over the whole China region because RH strongly influences the ADRF estimated within the surface boundary layer (Cheng et al., 2008). Enhanced warming within the atmosphere is not surprising because an abundance of absorbing aerosol particles like BC is predominant in the PRD region (Wu et al., 2009). The F TOA can change from negative to large positive forcing depending on the cloud fraction, surface albedo, and the vertical distributions of aerosols and clouds. However, regional TOA direct forcing is highly uncertain, particularly in regions containing atmospheric pollution like India (Ramanathan et al., 2001) and China (Li et al., 2010). The magnitude of F TOA was higher than values reported for the whole of China (Li et al., 2010), but similar to that in the Yangtze Delta region, which implies that abundant anthropogenic aerosols further cool the atmosphere-surface system over the polluted PRD region. However, the huge amount of solar radiation trapped inside the atmosphere is a significant source of heating for the atmosphere, particularly within the lower atmosphere, which can substantially alter the atmospheric stability and influence the dynamic system (Li et al., 2010).
The aerosol radiative FE is commonly used to quantify the radiative forcing potential of a given type of composite aerosol (Pathak et al., 2010). The diurnal mean FEs over the PRD and other regions are summarized in Table 2. The average FEs in this study were -60.0 7.8 W m-2 per AOD at the surface, -12.8 3.1 W m-2 per AOD at the TOA, and 47.3 8.3 W m-2 per AOD within the atmosphere. The magnitudes of these values are significantly higher than those reported in Taihu and over China as a whole. This is mainly due to the influence of strongly absorbing particles (e.g., carbonaceous aerosols), which reduce the net radiation at the surface and outgoing radiation at the TOA. As a result, atmospheric absorption increases. The magnitudes of the FEs are similar to those reported in Beijing, Xianghe, Yinchuan, East Asia, and at the Kaashidhoo Climate Observatory in the Maldives, but less than the values reported in South Asia and the Shandong Peninsula in China, where anthropogenic aerosols, BC, and organic carbon aerosols dominate (Wang et al., 2007; Xin et al., 2014).