Aerosol Optical Properties and Associated Direct Radiative Forcing over the Yangtze River Basin during 2001 – 2015

The spatiotemporal variation of aerosol optical depth at 550 nm (AOD550), Ångström exponent at 470–660 nm (AE470–660), water vapor content (WVC), and shortwave (SW) instantaneous aerosol direct radiative effects (IADRE) at the top-of-atmosphere (TOA) in clear skies obtained from the Moderate Resolution Imaging Spectroradiometer (MODIS) and Clouds and the Earth’s Radiant Energy System (CERES) are quantitatively analyzed over the Yangtze River Basin (YRB) in China during 2001–2015. The annual and seasonal frequency distributions of AE470–660 and AOD550 reveal the dominance of fine aerosol particles over YRB. The regional average AOD550 is 0.49 ± 0.31, with high value in spring (0.58 ± 0.35) and low value in winter (0.42 ± 0.29). The higher AOD550 (≥0.6) is observed in midstream and downstream regions of YRB and Sichuan Basin due to local anthropogenic emissions and long-distance transport of dust particles, while lower AOD550 (≤0.3) is in high mountains of upstream regions. The IADRE is estimated using a linear relationship between SW upward flux and coincident AOD550 from CERES and MODIS at the satellite passing time. The regional average IADRE is −35.60 ± 6.71 Wm−2, with high value (−40.71 ± 6.86 Wm−2) in summer and low value (−29.19 ± 7.04 Wm−2) in winter, suggesting a significant cooling effect at TOA. The IADRE at TOA is lower over Yangtze River Delta (YRD) (≤−30 Wm−2) and higher in midstream region of YRB, Sichuan Basin and the source area of YRB (≥−45 Wm−2). The correlation coefficient between the 15-year monthly IADRE and AOD550 values is 0.63, which confirms the consistent spatiotemporal variation patterns over most of the YRB. However, a good agreement between IADRE and AOD is not observed in YRD and the source area of YRB, which is probably due to the combined effects of aerosol and surface properties.


Introduction
Aerosol is one of the crucial factors contributing to regional and global climate change [1].It can directly absorb and scatter solar radiation at the top-of-atmosphere (TOA) and surface, which is defined as Aerosol Direct Radiative Effect (ADRE).The ADRE can be quantitatively calculated through the difference of upward shortwave (SW) solar fluxes without (Fclr) and with (Faero) aerosols for clear skies [2,3].The aerosol optical properties such as the aerosol optical depth (AOD) and the Ångström exponent (AE) are key parameters for estimating ADRE on the earth-atmosphere system [4].However, due to the diverse spatiotemporal distribution of aerosols, the ADRE remains the major source of uncertainty for assessing the regional climate [5].Therefore, to further explore aerosol effects on regional climate and air quality, long-term systematic estimation of aerosol optical properties and associated direct radiative effects needs to be made.
Ground and satellite remote sensing are common techniques widely used for measuring aerosol optical properties.With the improvement of the spatial and temporal resolution of sensors such as Multi-angle Imaging SpectroRadiometer (MISR), the Moderate Resolution Imaging Spectroradiometer (MODIS), and Clouds and the Earth's Radiant Energy System (CERES), the satellite remote sensing has an advantage over the ground-based approach for measuring aerosol properties at regional and global scales [6][7][8][9][10].Besides, efforts have been also made for estimating the ADRE through satellite remote sensing or model-simulation (i.e., radiative transfer model (RTM)).For the model-simulated method, aerosol optical properties obtained from ground-based or remote-sensing observations have been used as input to the RTM to calculate ADRE at TOA and surface [11][12][13][14][15][16].For the remote sensing method, one is to estimate ADRE by the linear relationship between AOD and concurrent upward SW flux at TOA obtained from satellite observations [17][18][19][20][21][22]; the other is to take Fclr as the minimum upward SW flux at TOA in clear sky during one month, but it is not suitable for the regions with high cloud coverage or aerosol loading such as the study region (YRB) [4].In other words, the remote-sensing technique (linear fitting) is more suitable for estimating ADRE at regional and global scales.
During the past decades, a series of studies on aerosol optical properties and associated direct radiative effects continue to emerge over Amazonia [3,17], Southeast Asia [18,19], Northeast India [11,20,[23][24][25][26][27][28][29], Central India [12], and global land [2,[30][31][32][33] and ocean [4,5,21].The results showed that the averaged ADRE at TOA was negative at global scale [4][5][6], but it varied regionally probably due to the comprehensive effects of aerosol optical properties and underlying surface characteristics.Recently, with the rapid development of urbanization and industrialization, China has become one of the major sources of anthropogenic aerosols.It triggers numerous studies for analyzing the aerosol optical and radiative properties in different regions of China.For instance, Sundström et al. [22] estimated ADRE at TOA using a linear regression of MODIS AOD and CERES SWF in Middle-East China; Xia et al. [14] analyzed the spatiotemporal distributions of aerosol optical and radiative properties at TOA and surface through ground-based data of 21 sites in China; Kang et al. [15] and Yu et al. [16] also reported the aerosol optical properties and ADREs based on RTM for Nanjing and Beijing, respectively.In general, previous studies in China mainly focused on the regions with intensive population and frequent industrial activities (short-term analyses), for example, Pearl River Delta (PRD), Yangtze River Delta (YRD) and Beijing-Tianjin-Hebei region.However, long-term systemic studies on the aerosol optical and radiative characteristics are still insufficient over the whole YRB due to the lack of observation sites, especially over the upstream and midstream regions of YRB.Furthermore, the midstream and downstream regions of YRB and Sichuan Basin have suffered from serious air pollution caused by the industrial production, vehicular emission, straw burning and the long-distance transport of dust particles from North China in spring [8][9][10][34][35][36].The various aerosol sources and diverse surface types make the aerosol optical and radiative characteristics over YRB more complicated [37].Therefore, to improve the knowledge of regional aerosol-climate interactions, it is necessary to estimate the aerosol optical properties and ADREs based on the satellite-retrieved method over the entire YRB.
In this study, comparisons between CARSNET and MODIS AOD 550 are implemented to test the performance of satellite data over the different regions inside YRB during 2007-2011.Integrated analyses of the spatial and temporal (annual, seasonal and monthly) distributions of aerosol optical properties (AOD 550 , AE 470-660 and WVC) retrieved from MOD08 over YRB during 2001-2015 are further performed.The spatiotemporal pattern of SW instantaneous ADRE (IADRE) at TOA for clear skies over YRB is quantitatively assessed by CERES-SFF data.Finally, a systemic uncertainty for IADRE is calculated and compared with previous studies.

Satellite and Ground Data
The China Aerosol Remote Sensing Network (CARSNET) is a ground-based remote sensing aerosol network with about 50 sites all over China, which is established by the China Meteorological Administration and other local meteorological institutes.Just like the Aerosol Robotic Network (AERONET), CARSNET has become a major resource for validating and analyzing aerosol properties over different areas of China [38].There are similar accuracies for retrieving AOD between CARSNET and AERONET; for example, the bias for AOD between Lin'an (CARSNET) and ZFU (AERONET), two ground sites nearby YRD, is less than 0.02 [39].The CARSNET retrieved AOD using the sun-photometer CE318 every 15 min at 4 wavelengths (440, 670, 870 and 1020 nm) for Level-1.0 (unscreened) and Level-1.5 (cloud-screened) data [40].In the present study, Level-1.5 aerosol data at the three CARSNET stations (Figure 1) over different regions of YRB are selected to verify the performance of MODIS aerosol products.The time periods are 2007-2010, 2009-2011 and 2007-2011 for Chengdu (urban site), Dongtan (rural site) and Changde (rural site), respectively [38].More details of CARSNET AOD are shown in

Satellite and Ground Data
The China Aerosol Remote Sensing Network (CARSNET) is a ground-based remote sensing aerosol network with about 50 sites all over China, which is established by the China Meteorological Administration and other local meteorological institutes.Just like the Aerosol Robotic Network (AERONET), CARSNET has become a major resource for validating and analyzing aerosol properties over different areas of China [38].There are similar accuracies for retrieving AOD between CARSNET and AERONET; for example, the bias for AOD between Lin'an (CARSNET) and ZFU (AERONET), two ground sites nearby YRD, is less than 0.02 [39].The CARSNET retrieved AOD using the sun-photometer CE318 every 15 min at 4 wavelengths (440, 670, 870 and 1020 nm) for Level-1.0 (unscreened) and Level-1.5 (cloud-screened) data [40].In the present study, Level-1.5 aerosol data at the three CARSNET stations (Figure 1) over different regions of YRB are selected to verify the performance of MODIS aerosol products.The time periods are 2007-2010, 2009-2011 and 2007-2011 for Chengdu (urban site), Dongtan (rural site) and Changde (rural site), respectively [38].More details of CARSNET AOD are shown in Table 1.The MODIS instrument aboard Terra (1999) and Aqua (2002) satellites provides land, atmosphere and ocean standard products (https://ladsweb.modaps.eosdis.nasa.gov/) in 36 spectral bands with a wide spectral range from 0.4 µm to 14.4 µm and at spatial resolutions of 250 m, 500 m and 1000 m [41].In this study, the daily AOD 550 values obtained from the Level-2.0 aerosol product (MOD04) are used.The Level-3.0 monthly atmosphere product (MOD08) at 1  1) are based on the Deep Blue algorithm (DB), which is suitable for retrieving aerosol optical properties over bright surfaces like urban and desert regions [42,43].
The CERES_SSF_FM1_MODIS_Ed3A product from CERES instrument on Terra (https://ceres.larc.nasa.gov/order_data.php)used in this study is a set of Level-2.0 data (Table 1).It can provide cloud and aerosol properties derived from MODIS, meteorological information offered by the Global Modeling and Assimilation Office (GMAO) and radiances in shortwave (0.3-5.0 µm), window (8.0-12 µm) and total (0.3-200 µm) broadband channels at a spatial resolution of 20 km [44].For this product, the measured radiances are converted to the spatiotemporal concurrent fluxes by the Angular Distribution Model (ADM) [45].In addition, in terms of the different spatial resolutions between MODIS (10 km × 10 km) and CERES (20 km × 20 km), the cloud and aerosol properties such as the clear area percentage and AOD 550 from MODIS are needed to be averaged into CERES footprint through the point spread function [45].Notice that, in the CERES_SSF_FM1_MODIS_Ed3A product, the MODIS aerosol and cloud masking datasets have been corrected into CERES footprint and can be extracted directly.

Validation between CARSNET and MODIS (DB) AODs
CARSNET AODs are measured at four wavelengths (440, 670, 870 and 1020 nm), while MODIS (DB) AODs are retrieved at three wavelengths (470, 550 and 660 nm).To estimate AODs at wavelength λ i , the Ångström equation is used: where AE and β are the Ångström exponent and the turbidity coefficient, respectively.Considering that λ 1 (440 nm) and λ 2 (870 nm) are not affected by water vapor, AE and β can be obtained as follows: Using Equations ( 1) and (2), we can obtain CARSNET AOD at 550 nm.
The MODIS (DB) AOD 550 values are averaged over 50 km × 50 km areas centered at ground-based sites, while the CARSNET AOD 550 values are averaged within ±30 min of the satellite passing time (Terra, 10:30 local time).The correlation coefficients of these two datasets at three stations (Chengdu, Dongtan and Changde) in China are analyzed.

Calculation of Instantaneous ADRE
In clear sky, SW IADRE at TOA can be estimated by where (lat, lon) is the given location coordinate, F aero (lat, lon) and F clr (lat, lon) are upward SW fluxes with and without aerosols, respectively at the satellite's passing time.F aero (lat, lon) can be directly observed from the CERES_SSF_FM1_MODIS_Ed3A product.However, since there is no case where the aerosol loading is zero at TOA in reality [4], F clr (lat, lon) needs to be estimated using the linear regression between the upward SW flux at TOA and the concurrent AOD.The linear regression is implemented as follows: (1) Read data: Read the upward SW flux (20 km × 20 km) at TOA and the concurrent MODIS AOD 550 (20 km × 20 km) from the CERES_SSF_FM1_MODIS_Ed3A product.It should be noted that the MODIS AOD 550 here has been corrected into CERES merged dataset and can be obtained directly from the CERES_SSF_FM1_MODIS_Ed3A product.According to Sena et al. [17], under the condition with excessive aerosol loadings (>2.0), multiple scattering of aerosols was likely to make the regression non-linearity, on the contrary, the relationship of upward SW flux and AOD 550 was usually linear in the low AOD condition (<2.0).Therefore, the corrected MODIS AOD 550 values are required to be smaller than 2.0.In addition, "longitude and colatitude of CERES FOV at surface", "CERES solar zenith at surface" and "Clear area percentage at sub pixel resolution" are also read from the CERES_SSF_FM1_MODIS_Ed3A product for the qualifications in the following steps.
(2) Select valid observations in each 1 • × 1 • grid: Divide the study region (YRB) into 167 grid cells of 1 • × 1 • , and select the valid observations (SWFs and AODs) in each grid cells.The qualifications of valid observations are as follows: (1) The clear area percentage at sub pixel resolution (CERES footprint) is restricted to more than 99.9% for getting rid of cloud contamination.(2) According to Sundström et al. [22], the correction of CERES upward flux due to the normalization of the solar zenith angle (SZA) varied from -15 to 15 Wm -2 in Middle-East China.Therefore, SZA is limited to less than 60 • for avoiding the oversized pixels.After collecting the valid observations for one month, the numbers of the valid values in each grid are counted.To build a successful linear regression, those grid cells containing less than 10 valid observations are eliminated.
(3) Build linear regression: Calculate the correlation coefficients of SWFs and AODs for each selected grid in Step (2).To build a successful linear regression, those grid cells with small absolute values of correlation coefficients (<0.2) are also eliminated [22].After that, the linear regression in each grid is built to estimate F clr , where SWF is the upward SW flux, F clr is referred to the y intercept of the linear regression (i.e., the upward SW flux for AOD = 0), and α is the slope of the linear regression.
In addition, we calculate the IADRE efficiency by regressions from median values of IADRE and AOD.It is conducive to further analyze the effects of other aerosol optical prosperities except AOD on IADRE [19].The detailed flowchart of our satellite-retrieved method for calculating IADRE at TOA over YRB during 2001-2015 is illustrated in Figure 2.

Validation of MODIS (DB) AOD versus CARSNET AOD
The regression analyses between MODIS (DB) and CARSNET AOD550 over Chengdu (2007-2010), Changde (2007)(2008)(2009)(2010)(2011) and Dongtan (2009-2011) are shown in Figure 3.The parameters are slope (M), intercept (C), number of matched data points (N), correlation coefficient (r) and root mean square error (RMSE).The r values are 0.82, 0.70 and 0.82 over Chengdu, Changde and Dongtan, respectively, implying the good agreement between MODIS (DB) and CARSNET AOD550.The high correlations also provide possibility and confidence to analyze aerosol optical and radiative properties over YRB by MODIS DB aerosol retrievals.The MODIS (DB) AOD retrieval at Chengdu site is better than Changde and Dongtan due to the closer to 1.0 slope (M = 0.94), higher correlation (R = 0.82) and smaller error (RMSE = 0.19).The reason may be that the DB algorithm has advantage of retrieving AOD over urban region (e.g., Chengdu) than vegetated rural areas (e.g.,

Temporal Variation of Regional Average
Figure 4 shows multi-year variation of monthly averaged aerosol optical properties (AOD550, AE470-660 and WVC) retrieved from monthly MOD08 during 2001-2015 over YRB.The annual mean AOD550 (Figure 4a) is about 0.49, which is higher than Kunming (0.28) [14], lower than YRD (0.61) and Anhui Province (0.63) [8].There are distinct spatial distributions for aerosol loadings in different regions of YRB, due to the large standard deviation of AOD550 (0.31).In Figure 4b, the multi-year variation of regional average AE470-660 slightly fluctuates, indicating a stable aerosol mode over YRB.Generally, higher (lower) AE470-660 values indicate the dominance of finer (coarser) aerosol particles.Therefore, most AE470-660 values are larger than 1.0 (annual mean value is 1.42), suggesting fine aerosol particles are major aerosol types over most regions of YRB. Figure 4c shows a periodic change in WVC with the annual mean value (2.13 ± 0.94 cm), which is closely related to the rainfall pattern over YRB.

Temporal Variation of Regional Average
Figure 4 shows multi-year variation of monthly averaged aerosol optical properties (AOD 550 , AE 470-660 and WVC) retrieved from monthly MOD08 during 2001-2015 over YRB.The annual mean AOD 550 (Figure 4a) is about 0.49, which is higher than Kunming (0.28) [14], lower than YRD (0.61) and Anhui Province (0.63) [8].There are distinct spatial distributions for aerosol loadings in different regions of YRB, due to the large standard deviation of AOD 550 (0.31).In Figure 4b, the multi-year variation of regional average AE 470-660 slightly fluctuates, indicating a stable aerosol mode over YRB.Generally, higher (lower) AE 470-660 values indicate the dominance of finer (coarser) aerosol particles.Therefore, most AE 470-660 values are larger than 1.0 (annual mean value is 1.42), suggesting fine aerosol particles are major aerosol types over most regions of YRB. Figure 4c shows a periodic change in WVC with the annual mean value (2.13 ± 0.94 cm), which is closely related to the rainfall pattern over YRB. Figure 5 shows the monthly and seasonal variations of regional average aerosol optical properties, including AOD550, AE470-660 and WVC over YRB for the period 2001-2015.The monthly AOD550 (Figure 5a) is higher in April (0.61 ± 0.36) and June (0.58± 0.37) and lower in December (0.34 ± 0.23) and November (0.39 ± 0.26).The seasonal mean AOD550 is 0.58 ± 0.35, 0.52 ± 0.29, 0.44 ± 0.29 and 0.42 ± 0.29 (Figure 5b) for spring (MAM), summer (JJA), autumn (SON) and winter (DJF), respectively.The maximum AOD550 in spring is probably associated with frequent dust events from North China [8,9,49], while the minimum in winter and autumn may be due to the windy and dry weather affected by Mongolian anticyclone [50].In summer, the higher rainfall reduces the presence of aerosol particles in the atmosphere, resulting in a decrease of AOD550 compared to spring, while the hygroscopic growth of fine-mode particles and the continuous biomass burning in June and July account for the increases in AOD550 compared to autumn and winter [34,35].
The monthly variation of regional average AE470-660 (Figure 5c) displays an opposite wave-like mode to AOD550, which increases from May (1.28 ± 0.38) and April (1.31 ± 0.37) to October (1.47 ± 0.24) and September (1.51 ± 0.25), indicating the dominance of fine particles over YRB.In Figure 5d, the seasonal mean AE470-660 in spring (1.33 ± 0.36) is lower than those in summer (1.42 ± 0.37), autumn (1.46 ± 0.27) and winter (1.45 ± 0.38), which is probably due to the long-distance transport of dust particles from North China in spring.Both monthly and seasonal variations of regional average WVC (Figure 5e,f) are unimodal curves, with peak values at July (4.05 ± 1.67 cm) and summer (3.79 ± 1.57 cm), which is associated with the Monsoon climate.
From the above, the monthly and seasonal variations of AOD550 over YRB are associated with the aerosol type (AE470-660), relative humidity (WVC), anthropogenic activity and terrain; these findings are consistent with previous investigations over Sichuan Basin [34] and Nanjing [15].However, a poor correlation coefficient (<0.1) was observed among AOD, AE and WVC at Lanzhou [51]. Figure 5 shows the monthly and seasonal variations of regional average aerosol optical properties, including AOD 550 , AE 470-660 and WVC over YRB for the period 2001-2015.The monthly AOD 550 (Figure 5a) is higher in April (0.61 ± 0.36) and June (0.58± 0.37) and lower in December (0.34 ± 0.23) and November (0.39 ± 0.26).The seasonal mean AOD 550 is 0.58 ± 0.35, 0.52 ± 0.29, 0.44 ± 0.29 and 0.42 ± 0.29 (Figure 5b) for spring (MAM), summer (JJA), autumn (SON) and winter (DJF), respectively.The maximum AOD 550 in spring is probably associated with frequent dust events from North China [8,9,49], while the minimum in winter and autumn may be due to the windy and dry weather affected by Mongolian anticyclone [50].In summer, the higher rainfall reduces the presence of aerosol particles in the atmosphere, resulting in a decrease of AOD 550 compared to spring, while the hygroscopic growth of fine-mode particles and the continuous biomass burning in June and July account for the increases in AOD 550 compared to autumn and winter [34,35].
The monthly variation of regional average AE 470-660 (Figure 5c) displays an opposite wave-like mode to AOD 550 , which increases from May (1.28 ± 0.38) and April (1.31 ± 0.37) to October (1.47 ± 0.24) and September (1.51 ± 0.25), indicating the dominance of fine particles over YRB.In Figure 5d, the seasonal mean AE 470-660 in spring (1.33 ± 0.36) is lower than those in summer (1.42 ± 0.37), autumn (1.46 ± 0.27) and winter (1.45 ± 0.38), which is probably due to the long-distance transport of dust particles from North China in spring.Both monthly and seasonal variations of regional average WVC (Figure 5e,f) are unimodal curves, with peak values at July (4.05 ± 1.67 cm) and summer (3.79 ± 1.57 cm), which is associated with the Monsoon climate.

Spatial Patterns of Temporal Variation of AOD
The spatial patterns of annual and seasonal mean AOD550 derived from MOD08 during 2001-2015 over YRB are illustrated in Figure 6.The high AOD550 (≥0.6) mainly appears over Sichuan Basin and the midstream and downstream regions of YRB.These high-AOD550 regions are usually centered on metropolises such as Shanghai, Wuhan and Chengdu (with intensive anthropogenic activities).Another potential reason for high aerosol loadings over Sichuan Basin is associated with its unique Basin topography, resulting in a low wind speed throughout the year, which is detrimental to the diffusion of aerosol particles.The low AOD550 (≤0.3) is mainly located in high Mountains over the upstream region of YRB, especially in Hengduan Mountains (<0.1) due to less human activities and dense forest vegetation [8].In other words, AOD550 decreases with increasing altitude and increases with rising human activities.This result is consistent with previous studies in Nanjing [15], YRD and Anhui Province [8].Recently, Xia et al. [14] analyzed the spatiotemporal distribution of aerosol optical properties over China using ground-based data at 21 stations; their results showed that the AOD550 over Middle-East China was generally larger than 0.5 with the lowest values (<0.1) appearing in Tibet Plateau.From the above, the monthly and seasonal variations of AOD 550 over YRB are associated with the aerosol type (AE 470-660 ), relative humidity (WVC), anthropogenic activity and terrain; these findings are consistent with previous investigations over Sichuan Basin [34] and Nanjing [15].However, a poor correlation coefficient (<0.1) was observed among AOD, AE and WVC at Lanzhou [51].

Spatial Patterns of Temporal Variation of AOD
The spatial patterns of annual and seasonal mean AOD 550 derived from MOD08 during 2001-2015 over YRB are illustrated in Figure 6.The high AOD 550 (≥0.6) mainly appears over Sichuan Basin and the midstream and downstream regions of YRB.These high-AOD 550 regions are usually centered on metropolises such as Shanghai, Wuhan and Chengdu (with intensive anthropogenic activities).Another potential reason for high aerosol loadings over Sichuan Basin is associated with its unique Basin topography, resulting in a low wind speed throughout the year, which is detrimental to the diffusion of aerosol particles.The low AOD 550 (≤0.3) is mainly located in high Mountains over the upstream region of YRB, especially in Hengduan Mountains (<0.1) due to less human activities and dense forest vegetation [8].In other words, AOD 550 decreases with increasing altitude and increases with rising human activities.This result is consistent with previous studies in Nanjing [15], YRD and Anhui Province [8].Recently, Xia et al. [14] analyzed the spatiotemporal distribution of aerosol optical properties over China using ground-based data at 21 stations; their results showed that the AOD 550 over Middle-East China was generally larger than 0.5 with the lowest values (<0.1) appearing in Tibet Plateau. Figure 7 shows the seasonal distribution of AOD550 over YRB; there is a similar spatial pattern with the annual pattern in Figure 6.However, the magnitudes vary greatly for each season.General speaking, the AODs550 in spring and summer are higher than those in autumn and winter except for the Sichuan Basin.Meanwhile, the seasonal range of AOD varies more significantly over high-value AOD areas.The highest AOD550 appears in spring over most of YRB due to the dust storms from North China.The low AOD550 in autumn and winter for most areas, expect Sichuan Basin, is probably caused by the windy and dry weather.In winter, high mountains around the Sichuan Basin block the invasion of cold air and make the local air flow inadequate, which does not aid the diffusion of the local aerosol particles.Furthermore, the high relative humidity caused by the winter morning fogs also accounts for high AOD550 in the Sichuan Basin [34]. Figure 7 shows the seasonal distribution of AOD 550 over YRB; there is a similar spatial pattern with the annual pattern in Figure 6.However, the magnitudes vary greatly for each season.General speaking, the AODs 550 in spring and summer are higher than those in autumn and winter except for the Sichuan Basin.Meanwhile, the seasonal range of AOD varies more significantly over high-value AOD areas.The highest AOD 550 appears in spring over most of YRB due to the dust storms from North China.The low AOD 550 in autumn and winter for most areas, expect Sichuan Basin, is probably caused by the windy and dry weather.In winter, high mountains around the Sichuan Basin block the invasion of cold air and make the local air flow inadequate, which does not aid the diffusion of the local aerosol particles.Furthermore, the high relative humidity caused by the winter morning fogs also accounts for high AOD 550 in the Sichuan Basin [34].

Frequency Distribution of Temporal Variation
Figure 8 illustrates the frequency distribution of annual mean AOD550, AE470-660 and WVC at 1° × 1° grid cell over YRB for the period 2001-2015.Approximately 92% of AODs are less than 1.0, while only about 1% of AODs are higher than 1.4.These are similar to the observations at Wuhan, where nearly 78% of AODs range from 0.40 to 1.23 [35].In Figure 8b, the frequency distribution of annual mean AE470-660 shows a single peak curve with the maximum frequency (32%) appearing around 1.5; AE470-660 more than 1.0 and less than 0.8 accounts for 90% and 3%, respectively.This indicates that fine aerosol particles are major aerosol types over YRB.In comparison, the annual mean AE470-660 at Wuhan roughly changes from 0.8 to 1.68, and the highest frequency (73.1%) is in the range of 1.05-1.55 [35]; the frequency of AE470-660 less than 0.6 is about 1%, while AE470-660 larger than 1.8 accounts for 0.3% over Shanghai [8].Besides, the annual mean WVC ranges from 0.2 to 7.6 cm and 71% of the WVC values are in the range of 0.2-3 cm (Figure 8c).
Figure 9a shows the seasonal frequency of AOD550; the high AOD550 values are gathered between 0.1-0.9, and the related frequencies are 81%, 88%, 93% and 94% during spring, summer, autumn and winter, respectively.In spring, the maximum frequency (21%) appears around 0.7 and 19% of the AOD550 values are higher than 0.9.There is a wide distribution of AOD550 in the range 0.1-2.3 during summer; 4% of the AOD550 values are less than 0.1 and nearly 6% of them are larger than 1.1, which are due to various aerosol sources (e.g., dust, smoke and pollution) in this season.The highest frequencies are 23% and 28% for AOD550 around 0.4 during autumn and winter, respectively.
The frequency distribution of AE470-660 for each season display a similar unimodal curve (Figure 9b); the maximum frequency is about 26%, 30%, 35% and 32% in the range of 1.5-1.6 during spring, summer, autumn and winter, respectively.Besides, nearly 13% of the AE470-660 values are less than 1.0 in spring, but the AE470-660 values (<1.0)only account for 4%, 3% and 4% in summer, autumn and winter, respectively.This indicates that there are relatively large particles during spring compared to other seasons, mostly due to the frequent dust events in this season.
In Figure 9c, the seasonal mean WVC ranges from 1 to 8 cm in summer, while it is concentrated in a narrow range of 0-3.5 cm in winter.Approximately 64% of the WVC values are around 3-6 cm and 6% are larger than 6 cm during summer, while nearly 90% of the WVC values are less than 1.5 cm in winter.During autumn, almost 90% of WVC values vary from 0 to 4 cm and only 4% of the WVC are more than 5 cm.In spring, the highest frequency value of WVC is 20% around 2 cm and

Frequency Distribution of Temporal Variation
Figure 8 illustrates the frequency distribution of annual mean AOD 550 , AE 470-660 and WVC at 1 • × 1 • grid cell over YRB for the period 2001-2015.Approximately 92% of AODs are less than 1.0, while only about 1% of AODs are higher than 1.4.These are similar to the observations at Wuhan, where nearly 78% of AODs range from 0.40 to 1.23 [35].In Figure 8b, the frequency distribution of annual mean AE 470-660 shows a single peak curve with the maximum frequency (32%) appearing around 1.5; AE 470-660 more than 1.0 and less than 0.8 accounts for 90% and 3%, respectively.This indicates that fine aerosol particles are major aerosol types over YRB.In comparison, the annual mean AE 470-660 at Wuhan roughly changes from 0.8 to 1.68, and the highest frequency (73.1%) is in the range of 1.05-1.55 [35]; the frequency of AE 470-660 less than 0.6 is about 1%, while AE 470-660 larger than 1.8 accounts for 0.3% over Shanghai [8].Besides, the annual mean WVC ranges from 0.2 to 7.6 cm and 71% of the WVC values are in the range of 0.2-3 cm (Figure 8c).
Figure 9a shows the seasonal frequency of AOD550; the high AOD550 values are gathered between 0.1-0.9, and the related frequencies are 81%, 88%, 93% and 94% during spring, summer, autumn and winter, respectively.In spring, the maximum frequency (21%) appears around 0.7 and 19% of the AOD 550 values are higher than 0.9.There is a wide distribution of AOD 550 in the range 0.1-2.3 during summer; 4% of the AOD 550 values are less than 0.1 and nearly 6% of them are larger than 1.1, which are due to various aerosol sources (e.g., dust, smoke and pollution) in this season.The highest frequencies are 23% and 28% for AOD 550 around 0.4 during autumn and winter, respectively.
The frequency distribution of AE 470-660 for each season display a similar unimodal curve (Figure 9b); the maximum frequency is about 26%, 30%, 35% and 32% in the range of 1.5-1.6 during spring, summer, autumn and winter, respectively.Besides, nearly 13% of the AE 470-660 values are less than 1.0 in spring, but the AE 470-660 values (<1.0)only account for 4%, 3% and 4% in summer, autumn and winter, respectively.This indicates that there are relatively large particles during spring compared to other seasons, mostly due to the frequent dust events in this season.
In Figure 9c, the seasonal mean WVC ranges from 1 to 8 cm in summer, while it is concentrated in a narrow range of 0-3.5 cm in winter.Approximately 64% of the WVC values are around 3-6 cm and 6% are larger than 6 cm during summer, while nearly 90% of the WVC values are less than 1.5 cm in winter.During autumn, almost 90% of WVC values vary from 0 to 4 cm and only 4% of the WVC are more than 5 cm.In spring, the highest frequency value of WVC is 20% around 2 cm and nearly 90% of the WVC values are in the range of 0-3 cm.These seasonal frequencies of WVC are consistent with the seasonal variations of monsoon climate over YRB.         2) are all negative, indicating that aerosols scatter solar shortwave radiation back to the space and hence produce a cooling effect at TOA.In this study, the annual mean IADRE over the whole YRB based on the satellite-retrieved method is compared with previous studies (Table 3); for example, the annual mean IADRE was −12.8 ± 2.1 Wm −2 over Southeast Asia using CERES_SSF_Terra datasets from 2001-2010 [19]; similarly, the seasonal IADRE values varied from −4.43 Wm −2 in post-monsoon to −70.66 Wm −2 in monsoon using CERES_SSF_Aqua datasets over Northeast India during 2002-2013 [20].Additionally, the linear regression between the 15-year monthly MODIS AOD 550 and IADRE is also shown in Figure 10b.The IADRE has a good correlation with AOD 550 (r = 0.63), representing that the cooling effect at TOA is partly caused by aerosol loading over YRB.  2) are all negative, indicating that aerosols scatter solar shortwave radiation back to the space and hence produce a cooling effect at TOA.In this study, the annual mean IADRE over the whole YRB based on the satellite-retrieved method is compared with previous studies (Table 3); for example, the annual mean IADRE was −12.8 ± 2.1 Wm −2 over Southeast Asia using CERES_SSF_Terra datasets from 2001-2010 [19]; similarly, the seasonal IADRE values varied from −4.43 Wm −2 in post-monsoon to −70.66 Wm −2 in monsoon using CERES_SSF_Aqua datasets over Northeast India during 2002-2013 [20].Additionally, the linear regression between the 15-year monthly MODIS AOD550 and IADRE is also shown in Figure 10b.The IADRE has a good correlation with AOD550 (r = 0.63), representing that the cooling effect at TOA is partly caused by aerosol loading over YRB.The monthly IADRE (Figure 11a) ranges from a maximum in April (−41.95 ± 6.38 Wm −2 ) to a minimum in December (−24.76 ± 7.33 Wm −2 ).Generally, the seasonal mean IADRE values (Figure 11b) are higher in spring (−40.01 ± 8.30 Wm −2 ) and summer (−40.71 ± 6.86 Wm −2 ), but lower in autumn (−32.52 ± 6.52 Wm −2 ) and winter (−29.19 ± 7.04 Wm −2 ), which are similar to the seasonal variation of MODIS AOD 550 (Figure 5b).The relatively higher IADRE and larger AOD 550 in spring and summer probably illustrate the connection between aerosol loading and cooling effect at TOA.In addition, the monthly and seasonal variations of IADRE over YRB are also associated with the periodic change of WVC (Figure 9).According to previous research [17,22], IADRE is estimated by the linear regression of CERES AOD and SWF in the condition of unrealistically low WVC.This may lead to an overestimation of IADRE due to the neglect of radiation attenuation caused by WVC [17].Besides, the aerosol hygroscopic growth caused by WVC increases the value of AOD in summer (Figures 8 and 9), hence results in the large negative IADRE in this season over YRB.Compared with other regions in China (Table 3), the seasonal variations of IADRE over YRB are consistent with previous research.For example, the 24 h-averaged ADRE values in Nanjing based on the RTM were −21.9 Wm −2 , −30.0 Wm −2 and −31.5 Wm −2 during autumn, summer and spring, respectively [15]; the model-retrieved ADREs in Beijing were in the range of −66 Wm −2 to −116 Wm −2 and −13 Wm −2 to −23 Wm −2 for dust days and non-dust days, respectively [16].
Remote Sens. 2017, 9, x FOR PEER REVIEW 15 of 23 The monthly IADRE (Figure 11a) ranges from a maximum in April (−41.95 ± 6.38 Wm −2 ) to a minimum in December (−24.76 ± 7.33 Wm −2 ).Generally, the seasonal mean IADRE values (Figure 11b) are higher in spring (−40.01 ± 8.30 Wm −2 ) and summer (−40.71 ± 6.86 Wm −2 ), but lower in autumn (−32.52 ± 6.52 Wm −2 ) and winter (−29.19 ± 7.04 Wm −2 ), which are similar to the seasonal variation of MODIS AOD550 (Figure 5b).The relatively higher IADRE and larger AOD550 in spring and summer probably illustrate the connection between aerosol loading and cooling effect at TOA.In addition, the monthly and seasonal variations of IADRE over YRB are also associated with the periodic change of WVC (Figure 9).According to previous research [17,22], IADRE is estimated by the linear regression of CERES AOD and SWF in the condition of unrealistically low WVC.This may lead to an overestimation of IADRE due to the neglect of radiation attenuation caused by WVC [17].Besides, the aerosol hygroscopic growth caused by WVC increases the value of AOD in summer (Figure 8 and 9), hence results in the large negative IADRE in this season over YRB.Compared with other regions in China (Table 3), the seasonal variations of IADRE over YRB are consistent with previous research.For example, the 24 h-averaged ADRE values in Nanjing based on the RTM were −21.9 Wm −2 , −30.0 Wm −2 and −31.5 Wm −2 during autumn, summer and spring, respectively [15]; the model-retrieved ADREs in Beijing were in the range of −66 Wm −2 to −116 Wm −2 and −13 Wm −2 to −23 Wm −2 for dust days and non-dust days, respectively [16].2) are all negative, there are still few positive IADRE values (Figure 12) over 1° × 1° grid cells, probably due to the comprehensive effect of aerosol optical properties and surface reflectivity.Generally, positive IADRE tends to be observed over regions with strong absorbing aerosols and high surface reflectivity [2,3].Assuming that the aerosol optical properties remain constant over bright surfaces, the aerosol layer appears "darker" than its background.This condition may reduce the upward SW flux (i.e., Faero), resulting in less negative or even positive IADRE at TOA.For example, during 2000-2009 in Amazonia, the annual mean ADRE calculated using linear regression between upward CERES SWF and MODIS AOD over cerrado (−4.6 ± 1.6 Wm −2 ) was less negative than the value over forest (−6.2 ± 1.9 Wm −2 ), mainly due to the higher surface reflectivity over cerrado [17].In addition, the absorbing aerosols can cause a warming (positive) effect at TOA by directly absorbing solar SW radiation.Overall, more than 99% of the grid cells in YRB have negative IADRE values (Figure 12), indicating cooling effects at TOA.
The spatial distribution of IADRE in Figure 12 matches well with AOD over Sichuan Basin and the midstream region of YRB, indicating that the high IADRE (>−45 Wm −2 ) is mainly due to the large aerosol loading (AOD550 > 0.5) over these areas.On the other hand, the IADRE values remain largely negative (>−45 Wm −2 ), while the AOD550 values are generally extremely low (<0.1) over the source area of YRB.In contrast, the large aerosol loadings (AOD550 > 0.6) over YRD do not  2) are all negative, there are still few positive IADRE values (Figure 12) over 1 • × 1 • grid cells, probably due to the comprehensive effect of aerosol optical properties and surface reflectivity.Generally, positive IADRE tends to be observed over regions with strong absorbing aerosols and high surface reflectivity [2,3].Assuming that the aerosol optical properties remain constant over bright surfaces, the aerosol layer appears "darker" than its background.This condition may reduce the upward SW flux (i.e., Faero), resulting in less negative or even positive IADRE at TOA.For example, during 2000-2009 in Amazonia, the annual mean ADRE calculated using linear regression between upward CERES SWF and MODIS AOD over cerrado (−4.6 ± 1.6 Wm −2 ) was less negative than the value over forest (−6.2 ± 1.9 Wm −2 ), mainly due to the higher surface reflectivity over cerrado [17].In addition, the absorbing aerosols can cause a warming (positive) effect at TOA by directly absorbing solar SW radiation.Overall, more than 99% of the grid cells in YRB have negative IADRE values (Figure 12), indicating cooling effects at TOA. reflectivity over different ecosystems of YRB.Usually, the IADRE values are more negative over dark surfaces with low reflectivity such as forest areas in the Southeast Tibet Plateau, than those over high reflective surface such as urban areas within YRD.Besides, the absorbing aerosols (e.g., black carbon particles) gradually increase with the rapid development of industrial activities over urban areas within YRD, which are helpful to absorb solar radiation and hence change IADRE from negative to less negative or even positive values.The spatial distribution of IADRE over each 1° × 1° grid cell within YRB (Figure 13) differs greatly for each season, and this is not well consistent with the seasonal spatial pattern of AOD550.In spring, the high IADRE values (>−45 Wm −2 ) appear in the midstream region of YRB, Sichuan Basin and the source area of YRB; the magnitude of negative IADRE in spring is higher than that in other seasons.In summer, the high IADRE values (>−45 Wm −2 ) are gathered from the source region of YRB to the Hengduan Mountain, which is possibly due to the comprehensive effects of surface properties and dust aerosols [52].It is revealed that dust aerosols in the Tibet Plateau were transported from the Taklimakan Desert during summer, resulting in a significantly cooling (negative) effect at TOA [53].However, IADRE over the Sichuan Basin in summer is less negative than that in winter, which may be due to larger aerosol loadings in winter than summer over the Sichuan Basin.There is a decreasing trend for IADRE except the Sichuan Basin from autumn to winter; then the IADRE values over the upstream region of YRB and the northern YRB are The spatial distribution of IADRE in Figure 12 matches well with AOD over Sichuan Basin and the midstream region of YRB, indicating that the high IADRE (>−45 Wm −2 ) is mainly due to the large aerosol loading (AOD 550 > 0.5) over these areas.On the other hand, the IADRE values remain largely negative (>−45 Wm −2 ), while the AOD 550 values are generally extremely low (<0.1) over the source area of YRB.In contrast, the large aerosol loadings (AOD 550 > 0.6) over YRD do not translate into high IADRE values (<−30 Wm −2 ).These interesting phenomena may be due to the reason that IADRE is influenced by the combination of aerosol optical properties and surface reflectivity over different ecosystems of YRB.Usually, the IADRE values are more negative over dark surfaces with low reflectivity such as forest areas in the Southeast Tibet Plateau, than those over high reflective surface such as urban areas within YRD.Besides, the absorbing aerosols (e.g., black carbon particles) gradually increase with the rapid development of industrial activities over urban areas within YRD, which are helpful to absorb solar radiation and hence change IADRE from negative to less negative or even positive values.
The spatial distribution of IADRE over each 1 • × 1 • grid cell within YRB (Figure 13) differs greatly for each season, and this is not well consistent with the seasonal spatial pattern of AOD 550 .In spring, the high IADRE values (>−45 Wm −2 ) appear in the midstream region of YRB, Sichuan Basin and the source area of YRB; the magnitude of negative IADRE in spring is higher than that in other seasons.
In summer, the high IADRE values (>−45 Wm −2 ) are gathered from the source region of YRB to the Hengduan Mountain, which is possibly due to the comprehensive effects of surface properties and dust aerosols [52].It is revealed that dust aerosols in the Tibet Plateau were transported from the Taklimakan Desert during summer, resulting in a significantly cooling (negative) effect at TOA [53].However, IADRE over the Sichuan Basin in summer is less negative than that in winter, which may be due to larger aerosol loadings in winter than summer over the Sichuan Basin.There is a decreasing trend for IADRE except the Sichuan Basin from autumn to winter; then the IADRE values over the upstream region of YRB and the northern YRB are decreasing to the minimum (<−10 Wm −2 ) in winter.One reason for the above may be that the surface reflectivity rises with the decreasing vegetation coverage in autumn and the presence of snow in winter [54].The relatively high concentration of absorbing black carbon aerosols in winter may be another reason [55].decreasing to the minimum (< −10 Wm −2 ) in winter.One reason for the above may be that the surface reflectivity rises with the decreasing vegetation coverage in autumn and the presence of snow in winter [54].The relatively high concentration of absorbing black carbon aerosols in winter may be another reason [55].

Temporal Distribution of IADRE Efficiency
The IADRE efficiency is defined as the value of IADRE per unit AOD (Wm −2 τ −1 ), which is a key parameter to analyze the climatic effect of aerosol.It can be estimated from the slope of linear regression between IADRE and concurrent AOD550 (from CERES_SSF product) in clear skies.However, the slopes are likely to be affected by outliers probably due to SZA, WVC, aerosol type and surface reflectivity [22].In order to minimize these effects, the median values of IADRE and AOD550 over each 1° × 1° grid cell are used for the linear regressions.The procedure for calculating the monthly IADRE efficiency over YRB for the period 2001-2015 is illustrated in Figure 14.The values of IADRE in each month appear more negative with increasing AODs.Generally, the IADRE efficiency becomes more negative with the lower fine mode fraction [19].Meanwhile, the lower surface reflectivity can also increase the magnitude of IADRE efficiency [20].
Figure 15 shows the monthly and seasonal shifts of regional average IADRE efficiency over YRB during 2001-2015.The monthly mean IADRE efficiency (Figure 15a [15].Furthermore, the error bars for the monthly and seasonal mean IADRE efficiencies exhibit large values in Figure 15, revealing large regional differences within YRB.

Temporal Distribution of IADRE Efficiency
The IADRE efficiency is defined as the value of IADRE per unit AOD (Wm −2 τ −1 ), which is a key parameter to analyze the climatic effect of aerosol.It can be estimated from the slope of linear regression between IADRE and concurrent AOD 550 (from CERES_SSF product) in clear skies.However, the slopes are likely to be affected by outliers probably due to SZA, WVC, aerosol type and surface reflectivity [22].In order to minimize these effects, the median values of IADRE and AOD 550 over each 1 • × 1 • grid cell are used for the linear regressions.The procedure for calculating the monthly IADRE efficiency over YRB for the period 2001-2015 is illustrated in Figure 14.The values of IADRE in each month appear more negative with increasing AODs.Generally, the IADRE efficiency becomes more negative with the lower fine mode fraction [19].Meanwhile, the lower surface reflectivity can also increase the magnitude of IADRE efficiency [20].
Figure 15 shows the monthly and seasonal shifts of regional average IADRE efficiency over YRB during 2001-2015.The monthly mean IADRE efficiency (Figure 15a 3) in the range of −48 to −59 Wm −2 τ −1 [14] and −26.05 to −45.46 Wm −2 τ −1 at Nanjing [15].Furthermore, the error bars for the monthly and seasonal mean IADRE efficiencies exhibit large values in Figure 15, revealing large regional differences within YRB.

Uncertainty Analysis
In this section, the uncertainty in IADRE from CERES and MODIS datasets is quantitatively assessed using the satellite-retrieved method.Obviously, it is a challenge due to diverse systematic uncertainties in calibration of CERES radiances (±0.4 Wm −2 ), unfiltering of CERES radiances at TOA (±0.4 Wm −2 ), conversion of CERES TOA radiances to fluxes by ADMs (±0.4 Wm −2 ), cloud contamination (±0.5 Wm −2 ) and estimation of SW fluxes at TOA without aerosols in clear sky (i.e., Fclr) [2,3,18,45].Among them, the calculation of Fclr is the largest error source for uncertainty.Since Fclr is estimated from the linear regression between CERES upward SW flux and concurrent MODIS AOD550, the systematic bias in MODIS AOD retrieval (i.e., the data on the "x-axis") will lead to a systematic uncertainty of the intercept (Fclr).According to previous studies [3,5,18], the systematic MODIS AOD bias is about ±0.05.Besides, the 15-year

Uncertainty Analysis
In this section, the uncertainty in IADRE from CERES and MODIS datasets is quantitatively assessed using the satellite-retrieved method.Obviously, it is a challenge due to diverse systematic uncertainties in calibration of CERES radiances (±0.4 Wm −2 ), unfiltering of CERES radiances at TOA (±0.4 Wm −2 ), conversion of CERES TOA radiances to fluxes by ADMs (±0.4 Wm −2 ), cloud contamination (±0.5 Wm −2 ) and estimation of SW fluxes at TOA without aerosols in clear sky (i.e., Fclr) [2,3,18,45].Among them, the calculation of Fclr is the largest error source for uncertainty.Since Fclr is estimated from the linear regression between CERES upward SW flux and concurrent MODIS AOD550, the systematic bias in MODIS AOD retrieval (i.e., the data on the "x-axis") will lead to a systematic uncertainty of the intercept (Fclr).According to previous studies [3,5,18], the systematic MODIS AOD bias is about ±0.05.Besides, the 15-year

Uncertainty Analysis
In this section, the uncertainty in IADRE from CERES and MODIS datasets is quantitatively assessed using the satellite-retrieved method.Obviously, it is a challenge due to diverse systematic uncertainties in calibration of CERES radiances (±0.4 Wm −2 ), unfiltering of CERES radiances at TOA (±0.4 Wm −2 ), conversion of CERES TOA radiances to fluxes by ADMs (±0.4 Wm −2 ), cloud contamination (±0.5 Wm −2 ) and estimation of SW fluxes at TOA without aerosols in clear sky (i.e., Fclr) [2,3,18,45].Among them, the calculation of Fclr is the largest error source for uncertainty.Since Fclr is estimated from the linear regression between CERES upward SW flux and concurrent MODIS AOD 550 , the systematic bias in MODIS AOD retrieval (i.e., the data on the "x-axis") will lead to a systematic uncertainty of the intercept (Fclr).According to previous studies [3,5,18], the systematic MODIS AOD bias is about ±0.05.Besides, the 15-year mean IADRE efficiency (the median value of slope for the regression) is estimated to be 65.18Wm −2 τ −1 .Therefore, the systemic uncertainty in Fclr caused by the linear regression is approximately translated to ±1.6 Wm −2 .Furthermore, assuming that the above-mentioned uncertainties (Table 4) are uncorrelated, the total uncertainty in IADRE can be estimated to be (0.4 2 + 0.4 2 + 0.4 2 + 0.5 2 + 1.6 2 ) 1/2 = 1.8 Wm −2 in the present study as estimated by [56]: where U t and U i represent the combined uncertainty and the independent uncertainty in IADRE, respectively [3,5,18].Based on the same methodology of estimating systematic uncertainties in IADRE, the Fclr (total) uncertainty was ±1.1 Wm −2 (±1.39 Wm −2 ) in Amazonia [3] and ±1.05 Wm −2 (±1.2 Wm −2 ) in Southeast Asia [18].
Undoubtedly, apart from the above-mentioned sources of uncertainties, the total uncertainties for IADRE are also affected by other complex factors such as WVC, Sun-Earth distance, aerosol type and surface reflectivity [57,58].To reduce uncertainties in the IADRE, instrument calibrations, retrieved algorithms and multivariate data fusion need to be further improved.

Summary and Conclusions
Long-term spatiotemporal distributions of the aerosol optical properties and IADREs at TOA over YRB were investigated using MODIS and CERES datasets.MODIS DB observations were also validated by CARSNET AOD 550 during 2007-2011 to verify the performance of satellite-retrieved aerosol data.The key findings in this study are summarized as follows: 1.

2.
The 15-year averaged AOD 550 , A E470-660 and WVC values were 0.49, 1.42 and 2.13 cm, respectively.High AOD 550 values generally appeared in spring and summer due to long-distance dust transport from North China and the local biomass burning, while low AOD 550 values were observed in winter and autumn as a result of dry and windy weather conditions.Seasonal AE 470-660 values were all larger than 1.0 with the minimum (1.33 ± 0.36) in spring, indicating the dominance of fine-mode particles over YRB.Spatially, high AOD 550 values (>0.6) were observed over the midstream and downstream regions of YRB and Sichuan Basin because of intensive industrial and urban activities, while low AOD 550 values (<0.3) were in high mountains over the upstream region of YRB.

3.
Approximately 92% of the AOD 550 values were less than 1.0; there was a wide distribution of AOD 550 in the range 0.1-2.3 in summer, probably due to the variety of aerosol sources.The frequency analyses indicated 90% of the AE 470-660 values to be higher than 1.0, and about 71% of the WVC values were in the range of 0.2-3 cm, and 64% of the WVC values were around 3-6 cm in summer.
The total uncertainty of IADRE was estimated to be ±1.8Wm −2 , from the calculation of Fclr, CERES radiances unfiltered, conversion of CERES TOA radiances to fluxes and cloud contamination.
This study clearly revealed how to estimate SW aerosol radiative effects at TOA using satellite-retrieved method without any RTM.To explore SW aerosol radiative effects on the earth-atmosphere system more thoroughly, combinations of model-simulated approaches and satellite-retrieved methods should be considered in next study.

Figure 1 .
Figure 1.The location of Yangtze River Basin (YRB) over China.The urban (Chengdu) and rural (Dongtan and Changde) China Aerosol Remote Sensing Network (CARSNET) stations are marked with red point and red square, respectively.

Figure 1 .
Figure 1.The location of Yangtze River Basin (YRB) over China.The urban (Chengdu) and rural (Dongtan and Changde) China Aerosol Remote Sensing Network (CARSNET) stations are marked with red point and red square, respectively.

Figure 2 .
Figure 2. The flow diagram for calculating the Shortwave (SW) instantaneous aerosol direct radiative effect (IADRE) at the top-of-atmosphere (TOA).The IADRE is estimated by linear fitting the concurrent Clouds and the Earth's Radiant Energy System (CERES) SW TOA fluxes and Moderate Resolution Imaging Spectroradiometer (MODIS) aerosol optical depth at 550nm (AOD 550 ) selected from the CERES_SSF_FM1_MODIS_Ed3A product in each 1 • × 1 • grid cell over YRB for the period 2001-2015.

Figure 5 .
Figure 5. Monthly and seasonal variations of regional average AOD550 (a,b), AE470-660 (c, d), and WVC (e,f) over YRB for the period 2001-2015.The error bars represent the standard deviations of AOD550, AE470-660 and WVC in each month or season.

Figure 5 .
Figure 5. Monthly and seasonal variations of regional average AOD 550 (a,b), AE 470-660 (c,d), and WVC (e,f) over YRB for the period 2001-2015.The error bars represent the standard deviations of AOD 550 , AE 470-660 and WVC in each month or season.

Figure 6 .
Figure 6.Spatial distribution of annual mean AOD550 over YRB for the period 2001-2015.

Figure 6 .
Figure 6.Spatial distribution of annual mean AOD 550 over YRB for the period 2001-2015.

Figure 7 .
Figure 7. Spatial distribution of seasonal mean AOD550 values over YRB for the period January 2001 -December 2015.

Figure 7 .
Figure 7. Spatial distribution of seasonal mean AOD 550 values over YRB for the period January 2001-December 2015.
Remote Sens. 2017, 9, x FOR PEER REVIEW 12 of 23 nearly 90% of the WVC values are in the range of 0-3 cm.These seasonal frequencies of WVC are consistent with the seasonal variations of monsoon climate over YRB.

3. 3 .
Spatiotemporal Distribution of IADRE 3.3.1.Temporal Distribution of IADRE The annual variations of SW IADRE at TOA over YRB from 2001 to 2015 are shown in Figure 10a and Table 2.The annual IADRE ranges from the highest value of −46.65 ± 11.32 Wm −2 in 2001 to the lowest value of −27.99 ± 6.12 Wm −2 in 2013, with the multi-year average value (−35.60 ± 6.71 Wm −2 ).The values of annual IADRE at TOA (Table

23 3. 3 .
Remote Sens. 2017, 9, x FOR PEER REVIEW 13 of Spatiotemporal Distribution of IADRE 3.3.1.Temporal Distribution of IADRE The annual variations of SW IADRE at TOA over YRB from 2001 to 2015 are shown in Figure 10a and Table 2.The annual IADRE ranges from the highest value of −46.65 ± 11.32 Wm −2 in 2001 to the lowest value of −27.99 ± 6.12 Wm −2 in 2013, with the multi-year average value (−35.60 ± 6.71 Wm −2 ).The values of annual IADRE at TOA (Table

Figure 10 .
Figure 10.Annual variation (a) and linear regression (b) of AOD550 and SW TOA IADRE over YRB for the period 2001-2015.

Figure 10 .
Figure 10.Annual variation (a) and linear regression (b) of AOD 550 and SW TOA IADRE over YRB for the period 2001-2015.

Figure 11 .
Figure 11.Monthly (a) and seasonal (b) variation of SW TOA IADRE over YRB for the period 2001-2015.The error bars represent the standard deviations of IADRE in each month or season.

Figure 11 .
Figure 11.Monthly (a) and seasonal (b) variation of SW TOA IADRE over YRB for the period 2001-2015.The error bars represent the standard deviations of IADRE in each month or season.

23 Figure 14 .
Figure 14.Calculation of the monthly IADRE efficiency using the median values of IADRE and AOD550 over YRB for the period 2001-2015.

Figure 15 .
Figure 15.Monthly (a) and seasonal (b) variation of IADRE efficiency over YRB during 2001-2015.The bars indicate the relative error in each month or season.

Figure 14 . 23 Figure 14 .
Figure 14.Calculation of the monthly IADRE efficiency using the median values of IADRE and AOD 550 over YRB for the period 2001-2015.

Figure 15 .
Figure 15.Monthly (a) and seasonal (b) variation of IADRE efficiency over YRB during 2001-2015.The bars indicate the relative error in each month or season.

Figure 15 .
Figure 15.Monthly (a) and seasonal (b) variation of IADRE efficiency over YRB during 2001-2015.The bars indicate the relative error in each month or season.

.
Remote Sens. 2017, 9, x FOR PEER REVIEW 3 of 23 analyses of the spatial and temporal (annual, seasonal and monthly) distributions of aerosol optical properties (AOD550, AE470-660 and WVC) retrieved from MOD08 over YRB during 2001-2015 are further performed.The spatiotemporal pattern of SW instantaneous ADRE (IADRE) at TOA for clear skies over YRB is quantitatively assessed by CERES-SFF data.Finally, a systemic uncertainty for IADRE is calculated and compared with previous studies.

Table 1 .
Characteristics of satellite and ground-based data used in this study.

Table 1 .
Characteristics of satellite and ground-based data used in this study.
• × 1 • spatial resolution is used to analyze the aerosol optical properties including AOD 550 , AE 470-660 and WVC during 2001-2015.It should be noted that AOD 550 and AE 470-660 obtained from both MOD04 and MOD08 (Table

Table 2 .
Statistics of multi-year monthly, seasonal and annual mean Shortwave (SW) instantaneous aerosol direct radiative effect (IADRE) at the top-of-atmosphere (TOA).

Table 2 .
Statistics of multi-year monthly, seasonal and annual mean Shortwave (SW) instantaneous aerosol direct radiative effect (IADRE) at the top-of-atmosphere (TOA).

Table 4 .
Statistics of the sources and magnitudes of diverse uncertainties in IADRE.