Distinct surface response to black carbon aerosols

For the radiative impact of individual climate forcings, most previous studies focused on the global mean values at the top of the atmosphere (TOA), and less attention has been paid to surface processes, especially for black carbon (BC) aerosols. In this study, the surface radiative responses to five different forcing agents were analyzed by using idealized model simulations. Our analyses reveal that for greenhouse gases, solar irradiance, and scattering aerosols, the surface temperature changes are mainly dictated by the changes of surface radiative heating, but for BC, surface energy redistribution between different components plays a more crucial role. Globally, when a unit BC forcing is imposed at TOA, the net shortwave radiation at the surface decreases by −5.87± 0.67 W m−2 (W m−2)−1 (averaged over global land without Antarctica), which is partially offset by increased downward longwave radiation (2.32± 0.38 W m−2 (W m−2)−1 from the warmer atmosphere, causing a net decrease in the incoming downward surface radiation of −3.56± 0.60 W m−2 (W m−2)−1. Despite a reduction in the downward radiation energy, the surface air temperature still increases by 0.25± 0.08 K because of less efficient energy dissipation, manifested by reduced surface sensible (−2.88± 0.43 W m−2 (W m−2)−1) and latent heat flux (−1.54± 0.27 W m−2 (W m−2)−1), as well as a decrease in Bowen ratio (−0.20±0.07 (W m−2)−1). Such reductions of turbulent fluxes can be largely explained by enhanced air stability (0.07± 0.02 K (W m−2)−1), measured as the difference of the potential temperature between 925 hPa and surface, and reduced surface wind speed (−0.05± 0.01 m s−1 (W m−2)−1). The enhanced stability is due to the faster atmospheric warming relative to the surface, whereas the reduced wind speed can be partially explained by enhanced stability and reduced Equator-to-pole atmospheric temperature gradient. These rapid adjustments under BC forcing occur in the lower atmosphere and propagate downward to influence the surface energy redistribution and thus surface temperature response, which is not observed under greenhouse gases or scattering aerosols. Our study provides new insights into the impact of absorbing aerosols on surface energy balance and surface temperature response.

Abstract. For the radiative impact of individual climate forcings, most previous studies focused on the global mean values at the top of the atmosphere (TOA), and less attention has been paid to surface processes, especially for black carbon (BC) aerosols. In this study, the surface radiative responses to five different forcing agents were analyzed by using idealized model simulations. Our analyses reveal that for greenhouse gases, solar irradiance, and scattering aerosols, the surface temperature changes are mainly dictated by the changes of surface radiative heating, but for BC, surface energy redistribution between different components plays a more crucial role. Globally, when a unit BC forcing is imposed at TOA, the net shortwave radiation at the surface decreases by −5.87 ± 0.67 W m −2 (W m −2 ) −1 (averaged over global land without Antarctica), which is partially offset by increased downward longwave radiation (2.32 ± 0.38 W m −2 (W m −2 ) −1 from the warmer atmosphere, causing a net decrease in the incoming downward surface radiation of −3.56 ± 0.60 W m −2 (W m −2 ) −1 . Despite a reduction in the downward radiation energy, the surface air temperature still increases by 0.25 ± 0.08 K because of less efficient energy dissipation, manifested by reduced surface sensible (−2.88 ± 0.43 W m −2 (W m −2 ) −1 ) and latent heat flux (−1.54 ± 0.27 W m −2 (W m −2 ) −1 ), as well as a decrease in Bowen ratio (−0.20±0.07 (W m −2 ) −1 ). Such reductions of turbulent fluxes can be largely explained by enhanced air stability (0.07 ± 0.02 K (W m −2 ) −1 ), measured as the difference of the potential temperature between 925 hPa and surface, and reduced surface wind speed (−0.05 ± 0.01 m s −1 (W m −2 ) −1 ). The enhanced stability is due to the faster atmospheric warming relative to the surface, whereas the reduced wind speed can be partially explained by enhanced stability and reduced Equator-to-pole atmospheric temperature gradient. These rapid adjustments under BC forcing occur in the lower atmosphere and propagate downward to influence the surface energy redistribution and thus surface temperature response, which is not observed under greenhouse gases or scattering aerosols. Our study provides new insights into the impact of absorbing aerosols on surface energy balance and surface temperature response. logical cycle (Ramanathan et al., 2001;Menon et al., 2002). The former is realized via absorbing solar radiation, causing positive effective radiative forcing (ERF) at the top of the atmosphere (TOA) and thus warming the climate (Ramanathan and Carmichael, 2008;Bond et al., 2013;Myhre et al., 2013b), while the latter is partly through modifying the microphysical properties of clouds (e.g., albedo and lifetime) (Koch and Del Genio, 2010;Bond et al., 2013;Boucher et al., 2013), which could further impact ERF. For instance, Menon et al. (2002) attributed the cooling and drying trends in north China in the second half of the 20th century to BC aerosols; Meehl et al. (2008) suggested that BC contributed to the precipitation change in India by altering the meridional temperature gradient.
For radiative impacts, however, most previous studies have only focused on TOA forcing. TOA forcing is useful in understanding the climate feedback, climate sensitivity, and future climate change (Andrews et al., 2012), but it is not necessarily predictive of the spatial pattern of surface temperature response, which is more related to surface radiative changes (Wild et al., 2004). One intriguing phenomenon for BC is that, on the global scale, BC could warm the surface even with reduced solar radiation and net radiation at the surface (Ramanathan and Carmichael, 2008), which is somewhat counterintuitive as higher surface temperature generally requires more incoming radiation at the surface. FAQ 7.2 of Boucher et al. (2013) briefly described the heating process induced by BC. Specifically, BC particles firstly heat the atmosphere and cause surface cooling locally, but then they warm both the surface and the atmosphere due to atmospheric circulation and mixing processes. When it comes to surface response, Ramanathan et al. (2001) suggested that the reduced solar radiation at the surface is possibly counteracted by reduced evaporation, which further perturbs the hydrological cycle. Krishnan and Ramanathan (2002) found that the source regions of haze are subject to cooling due to the absorption of solar radiation, whereas regions outside the source can be warming, thus contributing to overall global warming. Liepert et al. (2004) argued that aerosols and clouds could lead to weakened turbulent flux at the surface. Wilcox et al. (2016) reported a reduction of turbulent flux under BC aerosols at the surface and linked such responses to clouds. Based on model simulations, Myhre et al. (2018) concluded that BC aerosols can change the global hydrological cycle by suppressing sensible heat flux at the surface and attributed this suppression to the changes of air stability.
The published studies cited above provide informative insights into the surface radiative responses to BC aerosols, but our understanding is still incomplete, especially from the perspective of the surface energy balance. In this study, we aim to fill this gap and answer the following scientific questions. (i) How does the surface warming under BC aerosols differ from warming due to greenhouse gases and solar forcing? (ii) What are the specific mechanisms that drive such warming responses? (iii) what are the relative contributions to the surface temperature change from each surface energy budget component?
2 Data and methods

Data
This study employs the model output from the Precipitation Driver and Response Model Intercomparison Project (PDRMIP), utilizing simulations examining the climate responses to individual climate drivers (Myhre et al., 2017). The eight models used in this study are CanESM2, GISS-E2R, HadGEM2-ES, HadGEM3, MIROC, CESM-CAM4, CESM-CAM5, and NorESM. The versions of these models are essentially the same as their versions in the 5th Assessment Report of the Intergovernmental Panel on Climate Change (IPCC AR5). The configurations and basic settings are listed in Table 1. In these simulations, five separate perturbations were applied to all the models instantly on a global scale: a doubling of CO 2 concentration (CO 2 × 2), a tripling of CH 4 concentration (CH 4 × 3), a 2 % increase in solar irradiance (solar + 2 %), a 10-fold increase in present-day black carbon concentration/emission (BC × 10), and a 5-fold increase in present-day SO 4 concentration/emission (SO 4 × 5). Each perturbation was run in two parallel configurations, a 15-year fixed sea surface temperature (fsst) simulation and a 100-year coupled simulation. The former is compared with its fsst control simulation to diagnose the ERF at the TOA and fast responses in each model, whereas the latter is used to examine climate responses. One model (CESM-CAM4) used a slab ocean setup for the coupled simulation whereas the others used a full dynamic ocean. For aerosol perturbations, monthly year 2000 concentrations were derived from the AeroCom Phase II initiative (Myhre et al., 2013a) and multiplied by the stated factors in concentration-driven models. Some models were unable to perform simulations with prescribed concentrations. These models multiplied emissions by these factors instead ( Table 1). The aerosol loadings in the NorESM model for the two aerosol perturbations are shown in Fig. 1 for an illustrative purpose; the spatial patterns are similar for other models. In the BC experiment, the concentration is highest in east China (E. China), followed by India, tropical Africa and South America (S. America). In the current study, these four regions are referred to as source regions due to their high emissions while the US and Europe are defined as non-source regions due to their relatively low emissions. For the SO 4 experiment, the aerosols are mainly restricted to the Northern Hemisphere (NH), with the highest loading observed in E. China, followed by India and Europe. The eastern US also has moderately high concentrations. More detailed descriptions of PDRMIP and some PDRMIP findings are given in Samset et al. (2016), Myhre et al. (2017), and Tang et al. (2018).

Methods
In this study, we start from the surface energy balance. We restrict our discussions to land grids only because this is where most people live, and thus the temperature response over land is more important to the well-being of humans. The incoming radiative energy (R in ) includes In Eq. (1), ↓ SW represents downward shortwave radiation and ↑ SW represents reflected SW radiation. ↓ LW denotes the downward LW radiation. The law of energy conservation requires that the R in should be balanced by the outgoing energy (E out ): In Eq.
(2), ↑ LW is the outgoing longwave radiation, which is a function of temperature based on the Stefan-Boltzmann law. H , λE, and G denote sensible heat flux, latent heat flux, and ground heat flux, respectively. For latent heat flux (λE), λ is the specific latent heat of evaporation, and E is the evaporation rate. R in is defined as surface radiative heating, as it is the radiative input provided to the surface to raise the surface temperature (Wild et al., 2004). The surface responds to the imposed energy by redistributing the energy content through each E out component. Since R in is equal to E out , we have (3) The changes of each energy component, denoted by , are obtained by subtracting the control simulations from the perturbations using the data of the years 6-15 in each fsst simulation and the years 71-100 in each coupled simulation. The changes are then normalized by the ERF in the corresponding experiments to obtain the changes per unit global forcing. Negative values for the energy component, denoted by a negative sign, represent decreasing trends. The ERF values for each model are obtained from Tang et al. (2019) and are defined as the combination of net SW radiation plus the net LW radiation at the TOA in the fsst simulations . The multi-model mean (MMM) ERF values are 3.68 ± 0.09 W m −2 (CO 2 ×2), 1.15±0.09 W m −2 (CH 4 ×3), 4.21± 0.05 W m −2 (solar + 2 %), 1.20±0.28 W m −2 (BC ×10), and −3.63 ± 0.71 W m −2 (SO 4 × 5) for indicated experiments (mean ±1 standard error). The MMM changes are estimated by averaging all eight models' results. A two-sided Student t-test is used to examine whether the MMM results are significantly different from zero. The same process was repeated for all variables analyzed in the current study.

Incoming radiation and surface temperature changes under BC forcing
Figure 2a-c show the MMM changes of R in and its components for the fsst simulations of the BC experiment. The fsst simulations are analyzed because we mainly focus on the rapid adjustments when the forcing is instantly imposed. Rapid adjustments are generally referred to as the fast responses that affect the components of the climate system and modify the global energy budget indirectly. Unlike feedbacks, rapid adjustments do not operate through changes in the global mean temperature, and most are thought to occur within a few weeks (Boucher et al., 2013). Specifically, when a unit BC forcing is imposed at the TOA, the net surface SW radiation decreases by −5.87 ± 0.67 W m −2 due to the absorption of solar radiation by BC particles, whereas ↓ LW radiation shows an increase of 2.32 ± 0.38 W m −2 ( Fig. 2a  and b), as a result of the warmer atmosphere. When combined, R in still decreases by −3.56 ± 0.60 W m −2 on the global scale, with some positive changes only in high-latitude regions (Fig. 2c). However, the surface air temperature increases globally by 0.25 ± 0.08 K despite the decreased R in , except in the source regions where some slight cooling trends occurred (Fig. 2d). It is noted that these results are for land grids only. The pattern of cooling in the source regions and warming elsewhere agrees well with the findings reported by Krishnan and Ramanathan (2002). This type of change persists into the near-equilibrium state, where global mean temperature changes and associated feedbacks are included ( Fig. 2e-h). Due to the enhanced warming of the atmosphere and probably water vapor buildup, the ↓ LW radiation shows a stronger increase (Fig. 2f), making the R in mostly positive and therefore the temperature change positive, except for source regions (Fig. 2g and h). An open question is how the temperature increased with a decreasing R in in the rapid adjustment processes. In order to better understand the mechanisms behind this warming phenomenon, we will explore the E out components in the next section. The spatial patterns of R in and ↑ LW for the BC experiment were quite different ( Fig. 3d and i). Another notable feature is the significant reductions of sensible and latent heat flux in the BC experiment ( Fig. 3n and s), which is in agreement with previous studies (Wilcox et al., 2016;Myhre et al., 2018;Suzuki and Takemura, 2019). These changes are obvious when averaged globally ( Fig. 4a and Table 2). The R in decreases by −3.56 ± 0.60 W m −2 , and the H and λE decrease by −2.88 ± 0.43 W m −2 and −1.54 ± 0.27 W m −2 , respectively, making the energy partitioned to ↑ LW positive (0.86 ± 0.36 W m −2 ). In other words, although the radiative heating (R in ) decreases, convective and evaporative cooling decrease by a greater amount (124 % relative to R in ) owing to less efficient energy dissipation, thereby warming the surface and leading to a positive ↑ LW radiation. The reduction of turbulent fluxes (H and λE) is found for both source regions and non-source regions (Fig. 4b). In the source regions, the reductions of turbulent flux are nearly the same as the reduction of R in , making the temperature response negligible or only slightly negative. In the non-source regions, the reduction of turbulent fluxes exceeds the reduction of R in , making the temperature response positive. Another interesting phenomenon is that the reduction of H is greater than the reduction of λE for the BC experiment, both on the global scale and in the source regions (Fig. 4). The greater reduction of H indicates a decrease in Bowen ratio (β), defined as the ratio of sensible heat flux over latent heat flux. Globally, β decreases by −0.20 ± 0.07 in the BC case, and such a drop could reach −0.3 in the source regions. In comparison, under other forcing agents, the changes of β are much smaller (the global MMM changes within ±0.10). The larger change of H is somewhat contradicting to the common sense that λE dominates the turbulent flux on a global mean scale. This is because on a global scale, 85 % of λE is from the ocean (Schmitt, 2008). In this study, we only focus on land grids, in which the λE is largely suppressed.  Under other forcing agents, the spatial patterns of ↑ LW are similar to the patterns of R in . The changes of H and λE are relatively small and sometimes cancel each other out, making little contributions to temperature change compared with BC (Figs. 3 and 4a). Therefore, the temperature change ( ↑ LW) is dominated by the change of radiative heating ( R in ). It is worth noting that the decrease in λE in the CO 2 experiment is due to the physiological effect of vegetation and plantation (Fig. 3p), which is included in the PDRMIP models . It is also noted that the stronger responses in the BC scenario (Figs. 3, 4 and Table   2) could be partially related to its larger changes of surface radiative heating ( R in ) compared with other forcing agents. Taking CO 2 as an example, R in is 1.26 W m −2 , similar to its TOA forcing (1 W m −2 ), whereas for BC, R in is roughly 3 times larger (Table 2). Observations show that the surface forcing of BC could be 10 times larger than TOA forcing on regional scales (Magi et al., 2008), indicating that BC could cause stronger changes of surface forcing than TOA forcing relative to other forcing agents. The contributions of G are negligible (Fig. 3u-y) and will not be further discussed.  1.54 ± 0.14 0.54 ± 0.09 0.32 ± 0.06 0.44 ± 0.07 0.00 ± 0.02 0.24 ± 0.01 0.10 ± 0.02

Attribution of temperature change
In order to quantify the contributions of each component to T , we applied a multi-linear regression model to the MMM values of T and the energy components in each experiment, as T = a × R in + b × H + c × λE. Here R in represents the changes of radiative heating, and H and λE denote changes in the surface energy redistribution. All grid cells were given equal weight. The results are listed in Table 3. A point-wise comparison of original T and fitted T is shown in Fig. S1. These regressions reproduce the T fairly well, since the correlation coefficients between T and fitted T are all above 0.73, and most of the data points align along the one-to-one line. For CO 2 , CH 4 , and solar and sulfate aerosols, the coefficients of R in are 1 order of magnitude larger than the coefficients of the turbulent fluxes, suggesting that R in dominates the temperature change under these forcing agents. When it comes to BC, however, T is more sensitive to H , followed by R in and λE. With the regression coefficients, we estimated the contributions of each energy component to T (Fig. 5). In line with our previous results, T was dominated by surface heating ( R in ) for most forcing agents with a very limited role from turbulent fluxes. BC, nonetheless, is an exception. For BC aerosols, T is influenced by both surface heating and turbulent fluxes, with the cooling from the former being overwhelmed by the warming from the latter (Fig. 5n,  s, and x). The domain-averaged changes for the BC experiment are listed in Table 4. Globally, H produces 0.19 K warming, and λE leads to 0.07 K warming. The combined 0.26 K warming is offset by −0.19 K cooling attributed to reduced R in , producing a net warming of 0.07 K. In terms of percentage, H and λE contribute 73 % and 27 %, respectively, to the total warming. Such patterns are also seen on regional scales. The warming contributions from H were 40 % (US), 47 % (Europe), 76 % (E. China), 87 % (India),  68 % (Africa), and 82 % (S. America), and the remaining part was contributed by λE.

Mechanisms underlying the reduction of turbulent fluxes
The above analyses show that for most of the forcing agents, R in dominates the surface temperature response, while for BC, the surface energy redistribution also comes into play in modifying temperature response as a result of the significant reductions of turbulent fluxes. The next question is why turbulent fluxes decrease substantially in response to BC particles. According to the bulk parameterization of turbulent fluxes, the sensible heat flux is expressed as Q SH = ρC p C H U (T s − T a ) and latent heat flux as Q LH = ρL v C E U (q s − q a ). In these two equations, ρ is air density, C p and L v are air specific heat capacity and latent heat of vaporization, respectively; C H and C E are two exchange coefficients; U denotes surface wind speed; and (T s − T a ) and (q s − q a ) represent temperature gradient and humidity gradient between the surface and air, respectively. In the rapid adjustment stage, wind speed (U ) and temperature gradient are the two possible causes for the changes of sensible heat, while wind speed and humidity gradient (q s − q a ) are likely to drive the change of latent heat flux. Figure 6a-e show the MMM changes of lower tropospheric stability (LTS), defined as the potential temperature difference between 925 hPa and the surface. An enhanced stability is observed for the BC experiment; LTS is 0.07 ± 0.02 K averaged globally, in contrast to near-zero values from other experiments (−0.01 to 0.01 K). The LTS values for the NH are even larger, with 0.09 ± 0.02 K for BC and −0.01 to 0.01 K for other forcing agents. The enhanced LTS, which can significantly impact the sensible heat flux , arises from the fact that the BC layers warm faster relative to the surface due to BC absorption of solar radiation. The changes of LTS patterns are similar for 850 and 700 hPa (Fig. S2). Figure 6f-j portray the MMM changes of surface wind speed. BC causes a much larger decrease in wind speed with respect to other forcing agents, 0.05 ± 0.01 m s −1 globally compared with zero from other forcing agents. The reduction in wind speed explains the weakening in both sensible and latent heat fluxes according to the bulk parameterization. Figure 6k-o show the changes of humidity gradient (q s −q a ), defined as the specific humidity difference between the surface and 850 hPa. For CH 4 , solar radiation, and SO 4 , the gradient increases globally with values of 0.02 g kg −1 (W m −2 ) −1 , 0.01 g kg −1 (W m −2 ) −1 , and 0.01 g kg −1 (W m −2 ) −1 respectively, causing an increase in λE (Figs. 3 and 4). For CO 2 , the gradient shows slightly negative values (−0.002 g kg −1 per W m 2 ), corresponding to reduced λE (Figs. 3 and 4). In terms of BC, the humidity gradient increases by 0.06 ± 0.02 g kg −1 (W m −2 ) −1 , but with reduced λE flux, indicating that humidity gradient is not the primary driver of latent heat change. These analyses illustrate that humidity gradient may also influence latent heat flux for CO 2 , CH 4 , solar radiation, and scattering aerosols. For BC, on the other hand, change of wind speed should be the primary driver of the reduction of λE, and humidity gradient is of less importance. Now the last question is why the surface wind speed decreases under BC forcing. The first potential explanation is the abovementioned enhanced LTS. Jacobson and Kaufman (2006) has clearly demonstrated that the enhanced LTS and reduced turbulent exchange can reduce the turbulent kinetic energy and vertical transport of horizonal momentum, thereby reducing surface wind speed. The second possible explanation is that from the dynamical perspective, wind speed is controlled by the pressure gradient force (PGF), the Coriolis force, the gravitational force, and the frictional force. PGF is the driving force for atmospheric motion and is potentially the main driver for the changes of wind speed in the current idealized experiments. On a global scale, the excessive heating in the tropics with respect to middle and high latitudes causes PGF to point toward polar regions. Here we hypothesize that the decrease in temperature gradient between the Equator and poles under BC forcing weakened the PGF and slowed down the wind speed. Evidence in support of this hypothesis is found in Fig. 7a-e showing the zonal mean atmospheric temperature change. Mechanistically, BC caused a larger atmospheric heating in the middle latitudes of NH (30-60 • N) relative to tropics due to more of the BC forcing being located at middle latitudes of NH (Fig. 7d). The faster warming of middle latitudes weakened the PGF between the Equator and polar regions, as seen from the larger increase in geopotential height of 500 hPa in the middlelatitude regions (Fig. 7i). These patterns are not observed in the other experiments. The changes of geopotential height at other levels show similar results (Fig. S3).
To further understand the relationship between changes in wind speed and temperature, we defined a temperature gradient index (Allen et al., 2012) as 2 × T 30−60 − ( T 0−30 + T 60−90 ), where T is the mass-weighted (300-850 hPa) temperature response, and subscripts 0-30, 30-60, and 60-90 denote low (0-30 • N), middle (30-60 • N), and high (60-90 • N) latitudinal zones, respectively. When the index becomes more positive, middle latitudes warm faster, and a stronger reduction of wind speed is expected. The results for each individual model and experiment are shown in Fig. 8. A reasonably good correlation is seen in the BC scenario (r = −0.59): a larger change in the temperature gradient index corresponds to a stronger decrease in wind speed. The results for other experiments are mostly scattered around zero.
On regional and local scales, several other factors might also contribute to surface wind change (Wu et al., 2018). For instance, the aerosols in Asia have been reported to modify the land-sea temperature contrast and thus modify monsoon circulation (Xu et al., 2006;Meehl et al., 2008). The different phases of internal variability (e.g., El Niño-Southern Oscillation (ENSO) and North Atlantic Oscillation (NAO)) could modulate the circulations on interannual to multi-decadal timescales (Jerez et al., 2013;Hu and Fedorov, 2018). Bichet et al. (2012) suggested that changes in the surface roughness length may also change the wind speed. These factors are not considered in the present study.

Discussion and summary
Our analyses demonstrate that under BC forcing surface energy redistribution plays a vital role in modifying the surface temperature due to the changes in turbulent fluxes. The changes of turbulent fluxes are consequences of a downward influence from the atmosphere. The warming of BC layers in the atmosphere enhances air stability and reduces wind speed. As a result, the surface turbulent fluxes are suppressed. This mechanism is not observed for other forcing agents such as greenhouse gases and scattering aerosols. A similar "top-down" mechanism has been previously observed  in the solar forcing, in which the stratosphere ozone reacts to the UV part of the solar variability and produces additional heating, leading to changes of circulation in the stratosphere. The changes in the stratosphere modify tropical tropospheric circulation that may impact the surface climate (Haigh, 1996;Gray et al., 2010).
As noted in Sect. 3.1, our above analyses mainly focus on the rapid adjustments, which are part of ERF by definition (Boucher et al., 2013). For the BC experiment, these ad-justments drive the surface to respond to the forcing. Most of the changes seen in the rapid adjustment stage extend into the near equilibrium (Figs. S4-S5). For BC forcing, R in in near-equilibrium state is close to zero with large inter-model spread (0.20±1.19 W m −2 ). The equilibrium turbulent flux H and λE are lowered by −2.57 ± 0.39 and −1.54 ± 0.27 W m −2 , respectively, which are comparable in magnitude to changes in the rapid adjustment stage. Such reductions of the equilibrium turbulent fluxes are found in both source regions and non-source regions. Since less energy dissipated away from the surface, more energy (4.34 ± 0.74 W m −2 ) was partitioned into LW, warming the surface by 0.87 ± 0.13 K. Figure S6 shows the slow responses under each forcing, which are obtained by subtracting the rapid adjustments from the coupled simulations. The slow responses are driven by global mean temperature change alone. Interestingly, the spatial patterns are quite similar across different forcing agents. Our finding further confirms that it is the rapid adjustment that led to the different surface responses to BC.
Two limitations exist in our current study. First, the aerosol-cloud interactions could not be fully represented, because for the models with fixed aerosol concentration, the changes of cloud lifetime do not affect aerosols. Second, for the BC simulations, two models (MIROC and NorESM) include aerosol indirect effects while the remaining ones have only aerosol-radiation interactions included (instantaneous and rapid adjustments). The cloud effects in these two models may slightly modify the SW radiation at the surface (Tang et al., 2020), although the results from these two models do not differ qualitatively from the other models without those effects. We suggest that our conclusions are not sensitive to such cloud effects.
In summary, our study shows that for forcing agents such as greenhouse gas (GHG), solar radiation, and scattering aerosol, R in dominates the surface temperature response. For BC forcing, the surface energy redistribution also plays an important role. Under BC forcing, the energy is dissipated less efficiently from the surface to the lower atmosphere, which causes warming at the surface despite the reduced radiative heating. The reductions of sensible heat flux account for 73 % of the surface warming on a global scale and 40 %-80 % of the warming on regional scales, with the remaining part arising from the reductions of latent heat flux. Such reductions of turbulent fluxes can be explained by enhanced lower tropospheric stability and reduced surface wind speed. The former is attributed to a faster atmospheric warming relative to the surface, whereas the latter is associated with enhanced stability and reduced Equator-to-pole atmospheric temperature gradient. These analyses contribute to our understanding of the impact of absorbing aerosols on surface radiation and climate.
Code and data availability. The PDRMIP model output used in this study is available to the public through the Norwegian FEIDE data storage facility. For more information, please see http://cicero. uio.no/en/PDRMIP (Myhre and Forster, 2013). The data shown in the figures are available upon reasonable request.
Author contributions. TT designed this study. TT performed the data analysis and wrote the initial manuscript. All authors contributed to scientific discussion, result framing, and manuscript polishing.
Competing interests. The authors declare that they have no conflict of interest.
Disclaimer. Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. and the Newton Fund through the Met Office Climate Science for Service Partnership Brazil (CSSP Brazil).
Financial support. The research has been partially supported by the NSF of the United States, with the grant no. 14-19398 and grant AGS1933630.
Review statement. This paper was edited by Zhanqing Li and reviewed by two anonymous referees.