Can a Coupled Meteorology-chemistry Model Reproduce the Historical Trend in Aerosol Direct Radiative Effects over the Northern Hemisphere?

The copyright of individual parts of the supplement might differ from the CC-BY 3.0 licence.


Introduction
Solar radiation, the Earth's primary energy source, plays a crucial role in the climate system.The decadal variations in surface shortwave radiation (SWR) starting from the 1950s and related climate impacts have been well documented (Ohmura and Wild, 2002;Mercado et al., 2009;Ohmura, 2009;Wild et al., 2005Wild et al., , 2007;;Wild, 2009).Many studies suggest that such variations were caused by the changes in anthropogenic aerosol loading over this time period (Streets et al., 2006;Ruckstuhl and Norris, 2009;Ohmura, 2009).Atmospheric aerosols not only directly scatter or absorb solar radiation (McCormick and Ludwig, 1967), but also affect surface solar radiation indirectly by altering cloud optical properties and lifetime (Twomey, 1977;Albrecht, 1989).The global aerosol effective radiative forcing at the top of atmosphere (TOA) is estimated to have cooling effects from both aerosol-radiation interactions (−0.95 to +0.05 W m −2 ) and aerosol-cloud interactions (−1.2 to 0.0 W m −2 ), which are comparable in magnitude to the warming effects by anthropogenic greenhouse gases (2.83 W m −2 ) (IPCC, 2014).However, large discrepancies in aerosol radiative forcing exist amongst the results estimated by different approaches, particularly at a regional scale.Aerosol direct radiative ef-Published by Copernicus Publications on behalf of the European Geosciences Union.fects (DREs) simulated by global chemical transport models (CTMs) are 30-40 % smaller than those derived from measurements (Yu et al., 2006).Additionally, decadal changes in surface SWR are considerably underestimated by global climate models (GCMs) (Wild, 2009) for both clear-and all-sky conditions.For example, Ruckstuhl and Norris (2009) stated that the surface SWR trends in Europe under cloud-free conditions simulated by the IPCC-AR4 models are smaller on average than that indicated by the observational evidence, and such large discrepancies in sign and magnitude between modeled and observed trends may be caused by large uncertainties in historic emission inventories and associated aerosol burdens.Thus, to draw robust interpretations of longterm records and assessments of aerosol impacts on climate forcing, more accurate descriptions of atmospheric aerosol loading, optical properties, and spatial-temporal distribution are necessary.
A two-way coupled meteorology and atmospheric chemistry model, i.e., the Weather Research and Forecast (WRF) model coupled with the Community Multiscale Air Quality (CMAQ) model, has been developed by the US Environmental Protection Agency (EPA) (Pleim et al., 2008;Mathur et al., 2010Mathur et al., , 2014;;Wong et al., 2012;Yu et al., 2014;Wang et al., 2014).This model system can be applied as an integrated regional climate and chemistry model, serving as an important tool for downscaling future projections of global climate to higher resolution as well as assessing the interactions between atmospheric chemistry, radiation, and meteorology.A preliminary analysis on a 10-day WRF-CMAQ simulation of a wildfire event in California suggests that including the radiative effects of aerosols improves the accuracy of both the meteorology and air quality simulations (Wong et al., 2012).However, an assessment of the performance of such coupled models in reproducing the aerosol radiative effects is needed to build confidence in their use as regional climatechemistry models over decadal time periods.Decadal hemispheric WRF-CMAQ simulations from 1990 to 2010 were conducted and evaluated through comparison with long-term surface observations of gaseous and particle species in our recent study (Xing et al., 2015).The current study focuses on aerosol direct radiative effects (DREs) and aims to answer the following questions: (1) how well does the new model represent the regional and temporal variability of aerosol burden and DRE?And (2) is the model able to capture past trends in aerosol loading and the associated radiative effects?
A brief description of the model configuration and observations is given in Sect. 2. The evaluation of the model simulated historical AOD and clear-sky SWR is presented in Sect.3.1 and 3.2.The estimates of DRE are provided in Sect.3.3.Aerosol radiative efficiency is further discussed in Sect. 4.

Model configuration
This study focuses on the summer months (June, July, and August) over a 21-year period (1990-2010) using WRF-CMAQ (WRFv3.4coupled with CMAQv5.0)driven by internally consistent historical emission inventories obtained from EDGAR.Details about emission processing and WRF configuration are described in Xing et al. (2015).The strength of nudging coefficients for four-dimensional data assimilation and indirect soil temperature nudging employed in WRF have been tested and chosen to improve model performance for meteorological variables without dampening the effects of radiative feedbacks.The nudging coefficient for both u/v wind and potential temperature is set to 0.00005 s −1 , while 0.00001 s −1 is used for nudging of the water vapor mixing ratio.The simulation domain is shown in Fig. 1 and covers most of the Northern Hemisphere, discretized with a grid of 108 km × 108 km resolution and 44 vertical layers of variable thickness between the surface and 50 mb (Xing et al., 2015).For further analysis and comparison with measurements, we selected four sub-regions over land (including adjacent ocean areas that might be impacted by transport from land) and four sub-regions mostly over the ocean.Three of these land regions are mostly impacted by anthropogenic emissions, i.e., eastern China (ECH, 20-40 • N, 100-125 • E), the eastern US (EUS, 28-50 • N, 100-70 • W), and Europe (EUR, 35-65 • N, 10 • W-30 • E), while the fourth land region is largely impacted by dust emissions, i.e., the Sahara Desert and the Arabian Desert (SHR, 10-25 • N, 10 • W-50 • E).Out of the four ocean regions, the North Pacific (NPA, 30-50 • N, 150 • E-130 • W) is in between ECH and EUS, the North Atlantic (NAT, 35-50 • N, 60-15 • W) is in between EUR and EUS, and the remaining two regions are downwind of SHR, i.e., the central Atlantic (CAT, 10-25 • N, 60-15 • W) and the northern Indian Ocean (NIN, 10-25 • N, 55-75 • E).
Both the feedback and non-feedback cases were simulated using the same model configuration and initial condition obtained from our previous continuous 21-year WRF-CMAQ simulations (Xing et al., 2015), except that aerosol direct radiative effects updated in the rapid radiative transfer model for general circulation models (RRTMG) (Clough et al., 2005) were considered in the feedback case.Similar to Hogrefe et al. (2015), in the non-feedback case, no default or climatological aerosol profiles were provided to the RRTMG model employed in WRF.Therefore, there are no aerosol effects on the radiation calculations in the nonfeedback case.DRE is thus estimated as the difference between feedback and non-feedback simulations.In addition, the coupled model used in this study allows the DRE to affect the dynamical fields and also represent the subsequent modulation of aerosol quantities associated with the "updated" dy- namical fields, such as soil dust emission flux and photolysis rates that are calculated online with the dynamical model.
Aerosol indirect effects by altering cloud optical properties and lifetime were not considered in the current study.Efficiencies for extinction, total scattering, backscattering, and asymmetry factor for a single particle were calculated using the BHCOAT coated-sphere module approach (Bohren and Huffman, 1983).Aerosol effects are treated dynamically in the coupled system, where the CMAQ chemistry and radiation-feedback modules are called every 5 and 20 time steps of WRF, respectively.The time step of WRF was set to be 60 s in simulation; thus, the meteorology fields will be updated from the feedback module every 20 min in simulation.Thus, space-and time-varying aerosol optical properties estimated from simulated aerosol composition and size distribution are used in the RRTM-based radiation calculations and impact the simulated dynamical and chemical state of the atmosphere.

Long-term observations
Table 1 summarizes the observation data used in this study, including satellite-retrieved and surface-based measured AOD and clear-sky radiation.

Aerosol optical depth (AOD)
Historical satellite remote sensing of AOD has been recorded going back to the 1980s.Despite limitations in accuracy (Chu et al., 2002), the satellite-measured AOD has global distributions with relative long-time coverage that makes it suitable for evaluating model simulation performance (Roy et al., 2007;Liu et al., 2010;Wang et al., 2010Wang et al., , 2011) ) and for understanding the evolution of atmospheric aerosols (Chin et al., 2014).To minimize the influences from uncertainties of a single satellite product, we conducted our analysis based on all available satellite retrievals (or averages) with more robust estimates (e.g., the level-3 data set).The Advanced Very High Resolution Radiometer (AVHRR) provides the longestrunning global satellite observations, starting in 1981.Longterm AVHRR-retrieved AOD over the oceans (AVHRR retrievals are unavailable over land due to relatively large uncertainties associated with surface reflectance over the land surface) has been produced by the NOAA Climate Data Record (CDR) project (Chan et al., 2013;Zhao et al., 2013).Five satellite retrieval products that cover both land and ocean are also used in this study (Table 1).The Total Ozone Mapping Spectrometer (TOMS) is one of the earliest satellites providing AOD measurement products with data from 1979 to 2001 (Torres et al., 1998(Torres et al., , 2002)).However, it is not specifically designed to measure aerosols and thus has limited accuracy (Chin et al., 2014).Both the TOMS product and AVHRR product are only used here for trend analysis.The other four data sources are more recent satellite products, i.e., the Seaviewing Wide Field-of-view Sensor (SeaWiFS) (McClain et al., 1998;Hsu et al., 2012;Sayer et al., 2012), the Multiangle Imaging Spectroradiometer (MISR) (Kahn et al., 2005;2010) and the Moderate Resolution Imaging and Spectroradiometer (MODIS) on both NASA's Earth Observing System (EOS)-Terra and EOS-Aqua satellites (Kaufman et al., 1997;Remer et al., 2005Remer et al., , 2008;;Levy et al., 2010) that have improved accuracy mostly covering the period of the 2000s.The bright desert surfaces that are missing in the "dark target" retrieval (ver.5.1) in the MODIS AOD product are filled by the "deep blue" retrieval that uses the 412 nm channel of MODIS to enable the retrieval of AOD over bright surfaces over land (Hsu et al., 2004).An exception to this is for MODIS-Terra data after 2007 when there is no deep-blue retrieval available.
In conducting the modeled-observed comparisons, the satellite data sets are interpolated to the 108 km × 108 km grid of the CMAQ northern hemispheric domain as shown in Fig. 1.The grid matching was conducted by using the simple inverse distance weighting method.To be consistent with the specific time that the satellite crosses the Equator, especially for the four more recent satellites, i.e., MISR (10:30 local time), SeaWiFS (noon local time), MODIS-Terra (10:30 local time) and MODIS-Aqua (13:30 local time), we chose to extract the AOD from model outputs at 11:00 (local time).Model-simulated AOD is also compared with measurements from the worldwide ground-based Aerosol Robotic Network (AERONET) that has been widely used for satellite product validation and model evaluations because of its direct measurements and unified high standard for instrument calibration (Holben et al., 2001;Chin et al., 2014).We select eleven sites from the AERONET with relatively long-term historical records to enable comparison of AOD trends, as additional evidence demonstrating the decadal changes in the tropospheric aerosol burden.For the purpose of comparison with the model, all the AERONET AOD observations were converted to the 533 nm band based on the equation of the Ångström exponent.

Clear-sky radiation
Since this study focuses on the aerosol direct radiative effect, the clear-sky SWR product is examined in this analysis.The Clouds and the Earth's Radiant Energy System (CERES) instruments were developed for NASA's EOS report observations starting from 1997 (Wielicki et al., 1996(Wielicki et al., , 1998)).The CERES Energy Balanced and Filled (EBAF) data set provides satellite-derived clear-sky shortwave radiation at TOA over the globe on a monthly averaged basis start-ing from 2000 (Loeb et al., 2009(Loeb et al., , 2012)).Such observationderived clear-sky SWR at TOA with global coverage is suitable for comparison with model simulations (Satheesh and Ramanathan, 2000;Rajeev and Ramanathan, 2001;Anantharaj et al., 2010) and estimation of the DRE (Yu et al., 2006;Patadia et al., 2008) at a regional scale.
Efforts have also been made to establish records of observation-derived clear-sky shortwave radiation at the surface, for the purpose of reducing the uncertainties in the fluxes at the surface, which are significantly larger than those at the TOA (Wild et al., 2006).For surface observations with high frequency such as the World Climate Research Programs Baseline Surface Radiation Network (WCRP-BSRN), it is possible to process a stratification of records into cloudy and clear-sky periods on the basis of an advanced clear-sky detection algorithm (Long and Ackerman, 2000;Wild et al., 2005) and to validate the model (Wild et al., 2006;Freidenreich and Ramaswamy, 2011).Unlike those at the TOA, the surface fluxes cannot be directly measured by satellites.However, the CERES mission estimates a global clear-sky surface SWR through radiative transfer calculations using satellite-retrieved surface, cloud, and aerosol properties as input (Kato et al., 2013), which agrees with the surface ob- servations on a global mean level (Wild et al., 2013).The satellite-derived surface product has also been used in recent analyses (Hakuba et al., 2014).
For comparison with model simulations on a regional scale, we used the CERES-derived clear-sky shortwave radiations at both TOA and the surface.In addition, we used clear-sky SWR at seven sites of the Surface Radiation Budget Network (SURFRAD, a component of the BSRN network) processed by Gan et al. (2014).These seven sites are grouped into the eastern US (including Bondville, IL; Goodwin Creek, MS; Penn State, PA; and Southern Great Plains, OK) and the western US (including Table Mountain, CO; Desert Rock, NV; and Fort Peck, MT) for comparison with simulations.
In consideration of the limited length of the record, as in our previous study (Xing et al., 2015), this study only focuses on linear trends.The linear least square fit method was employed, and significance of trends was examined with a Student's t test at the 80, 90 and 95 % confidence levels (p = 0.2, 0.1 and 0.05).To be consistent with the available period of satellite products mostly covering the period of the 2000s, the two decades were separated into two time periods spanning, i.e., the 1990s (1990-2000) and 2000s (2000-2010) for analysis.

Comparison with satellite-retrieved AOD
Due to discrepancies in the sensitivities and retrieval algorithms of the various sensors, satellite-retrieved AODs from the different platforms are not necessarily consistent with each other.Figure 2 displays the spatial distribution of JJAaveraged AOD from multiple-platform satellite retrieval data sets as well as from the model simulation.The "polluted half-circle" that starts from the Sahara Desert to East Asia is evident in MISR, MODIS-Terra and MODIS-Aqua.Missing data are found in large areas of China and India in both TOMS-and SeaWiFS-retrieved AOD.Inconsistencies among satellite retrievals are also noticeable in the North Pacific Ocean as well as in the southwestern United States, where higher AOD is shown in MODIS-Terra and MODIS-Aqua but low or no value in other satellite retrievals.In general, the model captures the spatial pattern of AOD but underestimates the AOD levels in dust areas including the Sahara Desert, the Arabian Desert, the Indo-Gangetic basin in northern India that is affected by dust storms coming from the west during summer (Prasad and Singh, 2007) 5 1990 1992 1994 1996 1998 2000 2002 2004 2006 2008 2010 WRF-CMAQ(608) Satellite_avg MODIS_TERRA( 608) MODIS_AQUA( 608) SeaWiFS( 371) MISR( 608) TOMS( 478 3 1990 1992 1994 1996 1998 2000 2002 2004 2006 2008 2010 WRF-CMAQ(533) Satellite_avg MODIS_TERRA( 533) MODIS_AQUA( 533) SeaWiFS( 531) MISR( 533) TOMS( 517) 5 1990 1992 1994 1996 1998 2000 2002 2004 2006 2008 2010 WRF-CMAQ(787) Satellite_avg MODIS_TERRA( 787) MODIS_AQUA( 787) SeaWiFS( 555) MISR( 787) TOMS( 486 in which parameters related to the calculation of "threshold friction velocity" need future investigation (Fu et al., 2014).Ridley et al. (2014) found similar underestimation of summer AOD over the Sahel in their modeling study, suggesting that such biases might be associated with the Haboobs-driven dust emission, which is not well represented in the current dust module (Marsham et al., 2011).However, the model captures the moderate AOD level in the eastern US (EUS) and Europe (EUR) in early years and the high AOD level in eastern China (ECH) where anthropogenic aerosols are the dominant contributors.Slight underestimation in ECH and EUS might be associated with underestimation of fine particles (Xing et al., 2015), and poor representation of secondary organic aerosol formation is also more evident in summer (see Fig. 3a, b, c).
The model successfully captured AOD trends after 2000 in three non-dust land regions (i.e., ECH, EUS and EUR).The spatial distribution of the AOD trend (yr −1 ) and time series of region-averaged AOD are provided in Figs. 2 and  3, respectively.In ECH, a significant increase in the AOD level is suggested from satellite retrievals, with growth rates of +0.001 to +0.004 (+0.3 to +1.7 %) yr −1 .Such an increase in AOD has also been simulated by the model, though with a faster growth of +0.014 (+5.3 %) yr −1 (see Table 2).The overestimation of the AOD trend in ECH is caused by the discrepancy between the simulated and observed AOD trend over the southern part of eastern China, which can be strongly influenced by biomass burning in Southeast Asia (Deng et al., 2008;Fu et al., 2012).As displayed in Fig. 2, the model failed to capture the negative trend of the satelliteretrieved AODs in Southeast Asia.Variations of biomass burning activity in Southeast Asia are difficult to capture in the model without an accurate temporally resolved biomass emission inventory, currently not available.The declining trends shown in observed AOD might also be associated with the recession since late 2008 (Lin et al., 2010) that may be not well represented in the emission inventory.The SeaWiFS retrieved AOD presents a more significant (at the p = 0.1 level) increasing trend compared to other satellite products, though the growth rate (+0.004 yr −1 ) is still lower than that of the simulated AOD trend (+0.014 yr −1 ).The simulated declining trends of −0.004 (−2.7 %) yr −1 in EUS and −0.003 (−2.1 %) yr −1 in EUR are very close to the observed trends of −0.003 to −0.008 (−1.5 to −4.4 %) yr −1 and −0.002 to −0.007 (−1.0 to −3.8 %) yr −1 , respectively.Additionally, declining trends in both observed and simulated AODs in EUS and EUR are significant (at least at the p = 0.2 level).
For the period before 2000, the simulated AOD in ECH also shows an increasing trend of +0.002 (+1.0 %) yr −1 , which is smaller in magnitude than that after 2000.The opposite trend (though not significant at the p = 0.2 level) estimated from both AVHRR and TOMS might be explained by the limited number of grid values available for calculation (refer to Fig. 3).Both simulation and satellite retrievals (except TOMS in EUR, which has no trend) show declining trends in EUS and EUR before 2000, with magnitudes comparable to those for the period post 2000, and the declining trend in EUR is significant (at the p = 0.2 level) in both simulated and AVHRR-retrieved AOD.
The model captures the decreasing trends in the North Atlantic (NAT) for both the before-2000 and after-2000 periods (as seen in Table 2 and Fig. 3f).Such declining trends are attributable to decreases in both EUS and EUR.The AOD trends in the North Pacific (NPA) before 2000 could possibly be associated with a declining trend in global background AOD from biomass burning in Southeast Asia, dust variations in the Sahara Desert and anthropogenic emission reduction in Japan, Europe and North America.For the period before 2000, both simulated and satellite-retrieved AOD in NPA exhibit declining trends of −0.003 (−2.08 %) yr −1 and −0.001 to −0.002 (−0.5 to −0.6 %) yr −1 , respectively.However, after 2000, the increase in transport from ECH results in an increasing trend shown in both simulated and most satellite-retrieved AOD (significant at the p = 0.2 level for SeaWiFS and MISR) in NPA of +0.003 (+2 %) yr −1 and +0.001 to +0.002 (+0.3 to +1.7 %) yr −1 , respectively, except in AVHRR-retrieved AOD, which exhibits the opposite trend.
The model also captures declining trends in regions influenced mostly by natural mineral dust aerosol, i.e., SHR, CAT and NIN for the period before 2000.Such declining trends might result from a reduction in surface winds over dust source regions in Africa (Ridley et al., 2014).The AOD in CAT also shows a declining trend in most satellite retrievals for the period after 2000; however, the magnitude, i.e., up to −0.006 (−1.6 %) yr −1 , is much smaller than that for the period before 2000, which is up to −0.011 (−2.8 %) yr −1 .Simulations show a similar transition from stronger declining trends of −0.005 (−1.9 %) yr −1 before 2000 to smaller declining trends of 0.000 (−0.02 %) yr −1 after 2000.However, the model fails to capture the increasing trends in the Arabian Desert as well as its downwind ocean area NIN that can be observed from satellite-retrieved AOD (see Figs. 2,3h).Additionally, the simulated increasing trend in SHR has a smaller magnitude of +0.001 (+0.33 %) than satelliteretrieved AOD trends of +0.002 to +0.15 (+0.5 to +2.3 %), except for MODIS-Terra.This discrepancy might be caused by the poor representation of dust emissions over the Arabian Desert.

Comparison with AERONET measurement
The spatial location of the 693 AERONET sites can be seen in Fig. 4a, which shows the averaged summer AOD of all available records for each site.The corresponding simulated AOD averaged from the same selected periods that were used to average the observations at each site is given in Fig. 4b for comparison.Generally, the spatial gradients of AOD shown in the model are consistent with that from observations.However, the model tends to underestimate AOD both over regions where anthropogenic aerosols are dominant (ECH, EUS and EUR) as well as over regions where natural aerosols are dominant (SHR, CAT and NIN).Such a discrepancy is consistent with the comparisons against satellite retrievals in the previous section.
Eleven AERONET sites that have relatively long-term records were selected for trend analysis and are marked by red boxes in Fig. 4a.The simulated AOD as well as the satellite-retrieved AOD is chosen from the corresponding grid cell for each site based on its spatial location.A siteby-site comparison is presented in Fig. 5.
At the Beijing site (Fig. 5a), which is in northeastern China, the model agrees with satellite retrievals that show a continual increasing trend from 1990 to 2010.However, the AOD level from AERONET is much higher than both simulated and satellite-retrieved AOD, and shows an opposite trend.Such a discrepancy was also found in our previous analysis when comparing concentrations of precursors against surface observations (Xing et al., 2015).Extremely high AOD levels from AERONET also indicate that the coarse spatial resolution might limit the model's ability to represent the pollution distribution at a finer scale.Another site within the ECH area is Chen-Kung University (Fig. 5b) in Taiwan.Simulated and satellite-retrieved AOD trends show opposite directions, possibly because the model fails to capture the variations in biomass burning activities in Southeast Asia as discussed in the previous section.How- ever, the correlation coefficient is relatively high (R = 0.44) and the bias is relatively small (NMB = −5 %).The AOD measured by AERONET shows a similar increasing trend as the simulation, but the correlation is poor.Simulated, surface-observed and satellite-retrieved AODs are found to be consistent with each other in the other two land areas where anthropogenic aerosols are dominant, i.e., EUS and EUR.At the GSFC site (Fig. 5c) in EUS, all AODs show declining trends of −0.007 to −0.013 yr −1 , which is slightly higher than that of the region-averaged trend of −0.003 to −0.008 yr −1 .The correlation coefficient between the simulated and observed AOD ranges between 0.42 and 0.57, but the underestimation is > 40 % due to the underestimation of fine particles (Xing et al., 2015).At the Lille site (Fig. 5d) in EUR, simulated, surface-observed and satellite-retrieved AODs show consistent decreasing trends of −0.003 to −0.007 yr −1 , which are comparable with the region-averaged level.The correlation between the simulated and observed AOD ranges between 0.17 and 0.19 and the smaller NMB of < ±15 % compared to the EUS region likely is attributable to a better simulation of fine particle concentrations (Xing et al., 2015).At the HJAndrews site (Fig. 5e) on the western coast of the US where fewer emission reductions occurred than in the EUS, the simulated AOD shows a slightly increasing trend of +0.001 yr −1 that is also found in surface-observed and satellite-retrieved AOD.Such increasing trends might be associated with the increased AOD level in NPA after 2000.At Bonanza Creek (Fig. 5f), the model fails to capture the abnormal high AOD level observed by both AERONET and satellites for the years of 2004 and 2009, when large forest fires occurred in Alaska (National Interagency Fire Center, http://www.nifc.gov/fireInfo/fireInfo_stats_lgFires.html),re-sulting in a significantly enhanced PM burden.At Mauna Loa (Fig. 5g) located in the Hawaiian Islands, the model overestimates the AOD level compared to both AERONET and satellite observations, resulting from the uncertainties of windblown dust emissions.The AERONET-monitored AOD is even lower than the satellite retrievals, which might be explained by the high elevation (3400 m) of this AERONET site.
However, the model tends to underestimate AOD by 30-60 % in areas where mineral dust aerosols are dominant.At the Banizoumbou site (Fig. 5h) in the Sahel, south of the Sahara Desert, increasing trends of +0.004 to +0.017 yr −1 shown in both surface-observed and satellite-retrieved AOD are successfully captured by the model, which shows a similar increasing trend of +0.010 yr −1 .Also, the correlation between the simulated and observed AOD is fairly good (R > 0.5).At the Cape Verde site (Fig. 5i) in CAT where most satellites show continual decreasing trends, AOD observed by AERONET shows a slight increasing trend of +0.005 yr −1 .A similar increasing trend is also evident in the simulated AOD, i.e., +0.002 yr −1 .Another site where the model agrees with AERONET better than satellites is Sede Boker (Fig. 5j) in Israel, between the Sahara Desert and the Arabian Desert.The small declining trend of AOD observed by AERONET is captured by the model, with better correla-tion and smaller bias than compared with satellite-retrieved AOD.However, at Solar Village (Fig. 5k) in the Arabian Desert, the model fails to capture increasing trends of AOD observed in both AERONET and satellites, which is consistent with the finding we discussed in the previous section.

Comparison with CERES
The upwelling clear-sky SWR at TOA represents the solar radiation reflected back to space by surface albedo as well  as atmospheric aerosols.If surface albedo changes are negligible, its variation is indicative of the changes in top of the atmosphere cooling effects due to aerosols.The spatial distributions of observed and simulated clear-sky SWR at TOA are presented in Fig. 6 for 2000 and 2010.Higher values are shown in the areas with highly reflective surfaces, i.e., icecovered and desert land areas where ground-surface albedo is higher.Because the aerosol radiative effect has been considered in the feedback case, the simulated TOA SWR is higher than that in the non-feedback case, particularly in regions with high AOD levels such as ECH and SHR, and also compares better with the CERES data relative to the simulation without aerosol feedback effects.However, simulated clearsky TOA SWRs over desert areas, EUS, EUR and southern Asia are still lower than those derived by CERES due to the underestimation of natural and anthropogenic aerosols.
The spatial distribution of land-use types in the model is kept unchanged for the 21-year period.Therefore, there is no clear trend of SWR in the non-feedback case over the entire domain, except over the Arctic Ocean where simulated ice cover is melting, leading to a significant declining trend of TOA SWR.This is also evident in the CERES data.The TOA SWR trend in the feedback case that also represents the trend of DRE shows a similar spatial pattern as the AOD trend shown in Fig. 2. A comparison of the time series of the TOA SWR anomaly (anomaly in the average of 2000-2010, as seen in Fig. 7) also suggests better agreement with observations for the feedback simulation in which the correlation with observations is higher than that in the no-feedback case.
For the period of 2000-2010, the simulated TOA SWR from the feedback case exhibits a significant increasing trend of +0.301 W m −2 (+42 %) yr −1 in ECH, a decreasing trend of −0.095 W m −2 (−0.14 %) yr −1 in EUS and a decreasing trend of −0.106 W m −2 (−0.15 %) yr −1 in EUR (as seen in Fig. 7a-c and Table 3).Similar trends are evident in the CERES-derived TOA SWR, which shows a relatively smaller increasing trend of +0.103 W m −2 (+14 %) in ECH, but fairly comparable decreasing trends of −0.099 W m −2 (−0.14 %) and −0.116 W m −2 (−0.17 %) yr −1 in EUS and EUR, respectively.No observed satellite-derived TOA SWR is available back to the 1990s; however, the signs of the simulated TOA SWR trends agree well with the corresponding sign of the simulated and observed AOD trends in each region (Table 3).
The simulated TOA SWR in NPA (Fig. 7e) shows a declining trend of −0.092 W m −2 (−0.17 %) yr −1 before 2000 and an increasing trend of +0.107 W m −2 (+0.20 %) yr −1 after 2000, which are consistent with the AOD trends for Note: formatted entries are significant at the p = 0.05 level: italic denotes a significant decrease; bold denotes a significant increase.
Simulation results are all from WRF-CMAQ with feedback; simulated SWR is the monthly 24 h average.each period.The increasing trend in TOA SWR after 2000 is noticeable in CERES data as well, with a relatively smaller increasing trend of +0.057 W m −2 (+0.13 %) yr −1 .
In NAT (Fig. 7f), the simulated TOA SWR is decreasing for the whole period, with a stronger declining trend of −0.136 W m −2 (−0.25 %) yr −1 from 1990 to 2000 compared to the trend of −0.039 W m −2 (−0.07 %) yr −1 for the period after 2000.The decrease in TOA SWR in this region is associated with the corresponding reduction in AOD that also shows more significant decreasing for the period before 2000.However, the declining trend from the model for 2000-2010 is too weak to be observed in CERES data, which show a small increasing trend of +0.023 W m −2 (+0.05 %) yr −1 .The decreasing AOD in the Sahara Desert leads to a decreasing TOA SWR trend in the downwind ocean region CAT that is noticeable in both simulation and CERES data.More evidence of a decreasing trend in TOA SWR is shown in the simulation for the period before 2000 than that for the period after 2000 because the reduction of AOD in CAT is greater before 2000 than after.The model also captured the increasing TOA SWR trend in the Sahel that is also evident in CERES.
The significant declining trend in SHR for the period of 1990-2000 simulated by the model is associated with the corresponding decreases in AOD.However, the model fails to capture the slightly increasing trend observed by CERES over the Arabian Desert (as seen in Fig. 6a), also resulting in a failure to capture the increasing trend in NIN for the period of 2000-2010 (Figs. 6a, 7h) (Figs. 6a, 7h).This discrepancy is consistent with those noted in the AOD comparison, i.e., that simulated AOD in NIN shows a decreasing trend but observed AOD is increasing.
Higher surface SWR is shown in the mid-latitude regions where the solar zenith angle is high during June-August.The surface SWR is also influenced by the water vapor in the atmosphere; therefore, the desert areas (and downwind ocean regions) with relative dry air show higher surface SWR compared to other areas located at the same latitude (see Fig. 6b).Due to aerosol direct radiative effects, the simulated surface SWR decreases in the feedback case compared to the nofeedback case and is more comparable to the surface SWR derived from CERES.The surface SWR in ECH and EUS shows a noticeable reduction from the no-feedback to the feedback case, but is still overestimated compared to that derived by CERES.A more pronounced reduction in surface SWR caused by aerosol direct radiative effects is shown in regions dominated by natural dust aerosols, because the reduction of surface SWR in dust areas (and downwind ocean regions) is also caused by the decrease in near-ground albedo due to decreases in wind-blown dust emissions stemming from the lower wind speed, which is one of the climate responses to the aerosol radiative effects.
For the period of 2000-2010, the surface SWR simulated in the feedback case exhibits a decreasing trend of −0.56 W m −2 (−0.17 %) yr −1 in ECH (Fig. 7a) and an increasing trend of +0.097 W m −2 (+0.03 %) yr −1 and +0.141 W m −2 (+0.04 %) yr −1 in EUS (Fig. 7b) and EUR (Fig. 7c), respectively, with fairly good agreement with those derived from CERES (R > 0.5, which is much better than that calculated from the no-feedback case).Note that the simulated surface SWR trends in both EUS and EUR are considerably lower than CERES.
Similar to TOA SWR, the surface SWR trend in NPA also shows opposite directions for the periods before and after 2000 (Fig. 7e): an increasing surface SWR trend of +0.129 W m −2 (+0.04 %) yr −1 before 2000 but decreasing by −0.106 W m −2 (−0.03 %) yr −1 after 2000.Good agreement (R = 0.45) is found with CERES-derived surface SWR, which also exhibits a significantly declining trend of −0.334 W m −2 (−0.10 %) yr −1 from 2000 to 2010.The simulated increasing trend in NAT (Fig. 7f) is stronger before 2000 than for the period after 2000, whereas CERES-derived surface SWR shows a small trend in the opposite direction, which is consistent with the discrepancy found in TOA SWR.
The model fails to capture the increasing trend in surface SWR shown in most parts of SHR (Fig. 6b).However, the simulation has a good agreement with that derived from CERES (R = 0.49) from the series comparison of the surface SWR anomaly, suggesting that the variation of dust is too large to present a clear trend within a relatively short time period.Such a discrepancy might contribute to the bias in reproducing surface SWR in the downwind CAT region where the simulated and observed surface SWR trends present opposite directions for the period of 2000-2010.However, for the period of 1990-2000 when the reduction of dust aerosol in SHR is significant, the surface SWR in both SHR and CAT as expected presents increasing trends.In NIN (Fig. 7h), the decreasing trend from 2000 to 2010 shown in CERES-derived surface SWR is associated with the increase in dust aerosol from the Arabian Desert where the model fails to reproduce the observed AOD trend.The simulated surface SWR has a good correlation (R = 0.57) with CERES data, though shows a slight trend in the opposite direction.The increasing trend of surface SWR as well as the decreasing trend of TOA SWR in NIN for the period of 1990-2010 is associated with the decreasing trend of AOD for the same period.

Comparison with SURFRAD
To further investigate the model's ability to reproduce the historical clear-sky SWR at the surface, we compared results from simulations of feedback and no-feedback cases with observations at SURFRAD sites.The simulation data are selected from grid cells corresponding to the spatial locations of each SURFRAD site at the time of the measurement, and then grouped to a regional level, i.e., the eastern and western US.The CERES data are also selected for each corresponding site for comparison, but not necessarily for the same time periods, since CERES data are monthly mean including all hours.
In general, the simulated clear-sky surface SWR agrees better with SURFRAD than with CERES, with higher correlation coefficients (Fig. 8).The increasing trend in the   eastern US simulated by the model in the feedback case of +0.31 W m −2 yr −1 is more comparable with that observed by SURFRAD (+0.52 W m −2 yr −1 ) than CERES (+1.37 W m −2 yr −1 ), even though the CERES data are a monthly mean of 24 h, which is lower than daytime averages from both SURFRAD and the simulations.However, the model fails to reproduce some of the yearly variations (e.g., the sharp decrease during [1998][1999][2000] that are evident in SURFRAD, suggesting the need for simulations conducted on a finer spatial scale with more accurate spatially resolved emissions. The clear increasing trend of +0.62 W m −2 yr −1 shown in SURFRAD in the western US is associated with the increasing clear-sky diffuse radiation that might be influenced by factors other than aerosols (e.g., air traffic activities) (Gan et al., 2014).The model shows no trends because most reductions in aerosols occurred in EUS.CERES even shows a declining trend of −0.93 W m −2 yr −1 in the western US.The discrepancy in the two observation-derived trends suggests that there is a need to further improve the accuracy of observed surface SWR data as well.Note: formatted entries are significant at the p = 0.05 level: italic = significant decrease; bold = significant increase.The regions defined in Yu et al. (2006) are slightly different from this study, with ECH compared to land in zones 4 and 8 of Yu et al. (2006), EUS to land in zone 2, EUR to land in zone 3, SHR to land in zone 7, NPA to ocean in zones 1 and 4, NAT to ocean in zone 2, CAT to ocean in zone 6, and NIN to ocean in zone 7.

Aerosol direct radiative effect (DRE)
We estimated the DRE as the difference in clear-sky SWR with and without aerosols.It is easy to estimate the DRE from the model simulations by taking the difference of clearsky SWR between the feedback case and no-feedback case.
The simulated DRE at both TOA and surface is shown in Table 4 and Fig. 9. Also, we compared the modeled results with the values of measurement-based assessments summarized in Yu et al. (2006).All the values of "measured" DREs mentioned in subsequent discussions are obtained from Yu et al. (2006).Since it is impossible to remove all aerosols from the real atmosphere, the DRE cannot be directly observed.The measurement-based method needs the adoption of a radiative transfer model or other assumptions to estimate DRE.For example, Remer and Kaufman (2006) put the results of the MODIS aerosol retrieval as an internally consistent set of aerosol optical properties into a column radiative transfer climate model (i.e., CLIRAD-SW) to calculate the upwelling hemispheric broadband fluxes and the aerosol effects at TOA. Zhang et al. (2005a, b) calculated the shortwave aerosol radiative forcing over the global oceans by using the 20 km resolution CERES measurements as well as the aerosol-dependent angular distribution models.
The TOA-DRE is evident (| DRE | > 4 W m −2 ) in all eight regions except SHR, where the surface-reflected upwelling radiation can be reduced by the aerosol absorption, leading to a significant reduction of TOA-DRE (Kim et al., 2005).From the 1990s to the 2000s on a decadal average basis, the simulated TOA-DRE in ECH increased by 30 % due to the increased aerosol burden.A more significant increasing trend is shown in the 2000s compared to the 1990s.The simulated TOA-DRE in ECH shows comparable values with the measured one over the land of eastern Asia (Yu et al., 2006).The reduction of aerosols in EUS and EUR mitigated the TOA-DRE by 19 % and 20 % from the 1990s to the 2000s.The simulated TOA-DRE agrees with the measured one in EUR, but has lower values in EUS.Also, the simulated TOA-DRE in SHR is significantly lower than measurements.The underestimation of TOA-DRE in EUS and SHR might be associated with the underestimation of the aerosol burden (as indicated by the comparison of simulated and observed AOD).In NPA, the TOA-DRE shows similar opposite trends (before and after 2000), which is consistent with the trends of AOD and clear-sky TOA SWR, and is slightly lower than measurements but comparable with those simulated by GCMs as −3.6 to −11.7 W m −2 (Yu et al., 2006).The decrease in the AOD level in NAT also reduced the TOA-DRE by 9 %, with a more significant trend in the 1990s.The simu- The surface-DRE shows a similar spatial distribution as TOA-DRE, except that a stronger DRE is shown in SHR.Trends in surface-DRE at each site are consistent with those in TOA-DRE that have the same direction during the corresponding periods.Over land regions the simulated surface-DRE is about 2 times as much as the TOA-DRE for ECH, EUS and EUR, and about 6 times for SHR.On a global scale, the surface-DRE (−13.5 to −17.4 W m −2 ) derived by the satellite-model integrated approaches (Yu et al., 2006) is about 3.4 times as large as the TOA-DRE (−5.3 to −6.6 W m −2 ).The surface-DRE (−14.4 to −30.4 W m −2 ) derived by AERONET is about 2.2-3.8 times as large as the derived TOA-DRE (−5.2 to −11.1 W m −2 ) (Yu et al., 2006).Like many GCMs, the simulated surface-DRE over land is very likely to be underestimated in this study, because the simulated TOA-DRE by the model is comparable (in ECH and EUR) or lower (in EUS and SHR) than measurementderived values, but the ratio of surface-DRE/TOA-DRE estimates in this study is lower than measurement-derived values.
Over ocean regions, the simulated surface DRE (−7.3 to −18.9 W m −2 ) is about 1.4-1.8times as much as the TOA DRE (−5.1 to −11.1 W m −2 ).The ratio is comparable with that found from measurement-based estimates that show the global ocean averaged surface DRE (−9.3 to −11.9 W m −2 ) about 1.6 times as much as the measured TOA-DRE (−3.5 to −7.0 W m −2 ) (Yu et al., 2006).

Discussion
Successful estimates of DRE trends depend on accurate estimates of AOD as well as the aerosol direct radiative efficiency (E τ ), which is defined as the DRE per unit aerosol optical depth and has been used for comparisons between different methods (Yu et al., 2006).The E τ in this study can be estimated simply by using DRE divided by AOD.The relationship between DRE and AOD is close to linear under low-AOD conditions (Chung, 2012).Thus, the slope of the linear regression between SWR and AOD can provide a close estimate of E τ , denoted as E * τ (e.g., Satheesh and Ramanathan, 2000).
However, the E * τ (estimated as the slope of the linear regression between DRE and AOD) is smaller in magnitude (i.e., absolute value) than E τ at high AOD levels.This is illustrated in Fig. 10, which displays the relationship between average daytime clear-sky DRE and AOD (note that each point in these figures represents a daily average value).It is noticeable that the response of DRE to AOD becomes nonlinear at high AOD, more apparent in the relationships for the SHR region where the AOD is higher.Consequently, in this regime of high AOD, the estimated E * τ < E τ .Additionally, Yu et al. (2006) noted smaller radiative efficiencies for anthropogenic aerosols with higher absorbing components.This may also result in E * τ < E τ in regions of anthropogenic aerosols (such as ECH, EUS and EUR).To further explore the response of clear-sky SWR to the AOD level, we further analyzed data over an 11-year period (2000)(2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010).To minimize the influence of month-tomonth variability (e.g., associated with the solar zenith angle) between June, July and August, monthly averaged SWR and AOD were deseasonalized by subtracting the average of 11-year data for the corresponding month (see Fig. S1 in the Supplement).We chose 24 h-averaged SWR but AOD at noon (local time) to be consistent with the observationderived data from CERES-EBAF and satellite retrievals.The estimated radiative efficiency based on 24 h average SWR to noon AOD (denoted as E * τ 2 ) is about half of the value of E τ .Table 5 presents the estimates of E τ , E * τ and observed and simulated E * τ 2 at both TOA and the surface for eight regions.The measurement-based E τ summarized in Yu et al. (2006) (denoted as E τ -yu ) is used to compare the estimates in this study.
Over land regions where anthropogenic aerosols are dominant, the simulated E τ at TOA in ECH, EUS and EUR is about −45.4,−49.6 and −57.2 W m −2 τ −1 , respectively, which is slightly larger in magnitude than the measurementderived TOA-E τ that is about −9 to −33, −24 to −37 and −11 to −34 W m −2 τ −1 , respectively.Furthermore, the simulated TOA-E * τ 2 is also much higher than that derived from observations.One possible reason for this discrepancy is the moderate underestimation of AOD, since in most cases TOA-E τ is larger when AOD is low.Such overestimation of TOA-E τ may compensate for the underestimation of AOD, resulting in a comparable TOA-DRE in these regions.
The simulated E τ at the surface is about −90.8, −101.5 and −114.3W m −2 τ −1 for ECH, EUS and EUR, respectively, which is higher than the measurement-derived surface-E τ of −51 to −106, −65 to −84 and −57 to −98 W m −2 τ −1 .However, surface-E * τ has much lower values than surface-E τ ; thus, it is comparable with measurement-derived surface-E τ .The simulated surface-E * τ 2 also shows good agreement with observed estimates (−40 to −60 W m −2 τ −1 ).Discrepancies are found in the EUS where simulated surface-E * τ 2 is the lowest but observed surface-E * τ 2 is the highest of all regions.Lower surface-E * τ 2 in the US in the feedback case might be associated with the unexpected positive correlation (with slope = +14.97W m −2 τ −1 ) between surface SWR and AOD in the non-feedback case (see Table 5).With the agreement in surface-E τ , the underestimation of surface-DRE in EUS, EUS and EUR is primarily associated with the underestimation of the AOD level.
The feedback from reduction in TOA-DRE by aerosol absorption results in a relative lower TOA-E * τ 2 over SHR than those in other regions for both observation and simulation.However, the simulated TOA-E * τ 2 is twice as much as observed, suggesting overestimation of TOA-E τ in SHR as well.
Over the ocean regions, both simulated E τ and E * τ at TOA in NPA (−66.9 and −60.9 W m −2 τ −1 ) and NAT (−65.3 and −59.2 W m −2 τ −1 ) are higher than those measurementderived TOA-E τ that are about −31 to −52 and −32 to −41 W m −2 τ −1 , respectively.Also, the simulated TOA-E *  Eτ is calculated by using DRE divided by AOD.E * τ is the slope of DRE against AOD.E * τ 2 is calculated by deseasonalized SWR and AOD, observed SWR is from CERES, and observed AOD is the weighted average of four EOS satellite retrievals by the available number of grids for calculation; see Fig. 3. Correlation coefficients (R) between SWR and AOD are calculated and shown in parentheses for the Obs., Sim.-feedback and Sim.-no feedback cases.Eτ -yu is the data reported in Yu et al. (2006) that summarized both measurement-based and GCM estimations.The regions defined in Yu et al. (2006) are slightly different from this study, with ECH compared to land in zones 4 and 8 of Yu et al. (2006), EUS to land in zone 2, EUR to land in zone 3, SHR to land in zone 7, NPA to ocean in zones 1 and 4, NAT to ocean in zone 2, CAT to ocean in zone 6, and NIN to ocean in zone 7. is higher than that derived from observations.The overestimation (by 60-70 %) of TOA-E * τ 2 in NPA and NAT is not as significant as that in land regions (110-190 %); therefore, the simulated TOA-DRE over NPA and NAT is still within the range of measurement-based estimates.The simulated surface-E τ and E * τ for NPA (−96.9 and −68.8 W m −2 τ −1 ) and NAT (−98.7 and −59.2 W m −2 τ −1 ) are consistent with those measurement-derived values, which are −61 to −90 and −60 to −90 W m −2 τ −1 .Also, the simulated surface-E * τ 2 agrees with the observation-derived values.Therefore, the model successfully reproduces the observed surface-DRE over NAP and NAT.
In the other two ocean regions where natural dust is dominant (i.e., CAT and NIN), the simulated TOA-E τ and E * τ are overestimated compared to measurement-derived data.The simulated TOA-E * τ 2 is also higher than that derived from observations.The surface-E τ and E * τ tend to be overestimated as well in these regions, and the observed surface-E * τ 2 is much lower than that in other regions, whereas the simulated surface-E * τ 2 in these regions are comparable with others.A possible reason is that surface SWR shows a nonlinear response with AOD (showing a smaller slope, i.e., E τ ) when AOD levels are high.The model fails to capture the extremely high AOD observed by satellites, and consequently overestimates the E * τ 2 in these regions.Such overestimation of both TOA-and surface-E τ may compensate for the underestimation of AOD, resulting in a comparable TOA-and surface-DRE in these regions.

Conclusion
A 21-year simulation from 1990 to 2010 over the Northern Hemisphere was conducted with a coupled meteorologychemistry model.In general, the model captured historical AOD trends of 21 years in most regions, including the continual increasing trend in ECH, a decreasing trend for EUS, EUR and NAT, and the decreasing trend in the 1990s but increasing trend in the 2000s in SHR and NPA.The model also captured the decreasing trend in NIN before 2000, but failed to capture its increasing trend after 2000.That discrepancy as well as the underestimation of AOD over regions that have substantial natural dust aerosol contributions might be associated with the uncertainty in estimation of wind-blown dust in the model.Slight underestimations in ECH and EUS might be associated with underestimation of fine particles as well as the poor representation of secondary organic aerosol formation that is more evident in summer.The simulation with aerosol radiative effects successfully captured the historical clear-sky SWR at both TOA and surface in ECH, EUS and EUR.The model also captured the enhanced DRE in NPA at both TOA and surface after the 2000s.However, discrepancies are found over dust regions due to the uncertainties from AOD estimations in these regions.
The DRE estimation requires not only accurate estimates of AOD, but also of correct aerosol radiative efficiency.Unfortunately, neither has a certain value from measurementbased studies that still present a wide range of values (Yu et al., 2006;Chin et al., 2014).Uncertainties therefore limit the accuracy of DRE estimation.The same issue applies to the model as well.Estimates of TOA-DRE are comparable with those derived from measurements for most regions; however, overestimates of TOA-E τ may be compensating for the underestimations of AOD.Such overestimates of TOA-E τ might be associated with the uncertainties of the ratio of scatting aerosols to total aerosols (e.g., a larger TOA-E τ is expected when the change in absorption aerosols is larger than scattering aerosols).Compared to previous estimates from GCMs, the simulated surface-E τ in this study agrees better with those derived by the measurements.However, surface-DRE trends are underestimated due to underestimated AOD in land regions.Further improvement of simulated TOA-E τ as well as both anthropogenic and natural AOD is important.An accurate temporally resolved biomass emission inventory (van der Werf et al., 2006;Shi and Yamaguchi, 2014), an improved dust emission model (Kok et al., 2014) and an advance scheme to model atmospheric organic aerosol (Koo et al., 2014;Zhao et al., 2015) are suggested for future model investigations.We are currently conducting a similar study with a finer-scale simulation and relatively better spatially and temporally resolved emission inventories over the continental US domain.Further analysis of those model calculations and assessment of the impacts of the higher resolution emissions can be found in Gan et al. (2015).
The Supplement related to this article is available online at doi:10.5194/acp-15-9997-2015-supplement.

Figure 2 .
Figure 2. Spatial distribution of summertime AOD and its trend from satellite retrievals and the WRF-CMAQ model.
as well as the Gobi Desert in western China.Such underestimations as well as unexpected high AOD estimates in Hawaii might be associated with the uncertainty of the wind-blown dust module www.atmos-chem-phys.net/15

Figure 3 .
Figure3.Temporal series of regional JJA-averaged AOD from satellite retrievals and WRF-CMAQ.The number of grids (> 80 % coverage) involved in calculation is shown in the parentheses; only four EOS satellite-retrieved AODs, i.e., MISR, SeaWiFS, MODIS-Terra and MODIS-Aqua, were averaged into "Satellite-avg".

Figure 4 .
Figure 4. Comparison of spatial distributions of summer AOD between AERONET and WRF-CMAQ (dot size indicates the extent of data coverage for the period 1990-2010; all converted into 533 nm; sites marked by red boxes are selected for trend analysis).

Figure 5 .
Figure 5.Comparison of JJA-averaged AOD of selected AERONET sites.

Figure 6 .
Figure 6.Spatial distribution of summertime clear-sky shortwave radiation (SWR) from the CERES satellite and the WRF-CMAQ model: (a) at TOA (upwelling), and (b) at the surface (downwelling) (unit: W m −2 ; trends are computed on the basis of the JJA average over the 2000-2010 period with a linear least square fit method).

Figure 7 .
Figure 7. Observed and simulated time series of summertime clear-sky SWR at TOA and the surface (regional JJA average, anomaly in the average of 2000-2010; R is calculated from Sim.-feedback and Sim.-no feedback with CERES).

Figure 8 .
Figure 8.Comparison of SURFRAD/Ceres observed and WRF-CMAQ clear-sky SWR at the surface: (a) the eastern US and (b) the western US (site-averaged, anomaly in the average of 1995-2010).

Table 1 .
Summary of long-term observations used in this study.

Table 2 .
Trends in regional AOD from satellite retrievals and the WRF-CMAQ model.
P values: the probabilities that the trends are not statistically significant.No values when P > 0.2.Simulation results are all from WRF-CMAQ with feedback; simulated AOD is reported at 11:00 LT for each grid cell.

Table 3 .
Trends in regional SWR from satellite retrievals and WRF-CMAQ model.

Table 4 .
Estimates of 24 h-mean JJA-averaged aerosol direct radiative effect (DRE) and its trends.