SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
In this talk, five specific major GCE improvements: (1) ice microphysics, (2) longwave and shortwave radiative transfer processes, (3) land surface processes, (4) ocean surface fluxes and (5) ocean mixed layer processes are presented. The performance of these new GCE improvements will be examined. Observations are used for model validation.
This study presents results of comparisons between instrumental radiation data in the elevated Tibetan Plateau and two global satellite products: the Global Energy and Water Cycle Experiment - Surface Radiation Budget (GEWEX-SRB) and International Satellite Cloud Climatology Project - Flux Data (ISCCP-FD). In general, shortwave radiation (SW) is estimated better by ISCCP-FD while longwave radiation (LW) is estimated better by GEWEX-SRB, but all the radiation components in both products are under-estimated. Severe and systematic errors were found in monthly-mean SRB SW (on plateau-average, -48 W/sq m for downward SW and -18 W/sq m for upward SW) and FD LW (on plateau-average, -37 W/sq m for downward LW and -62 W/sq m for upward LW) for radiation. Errors in monthly-mean diurnal variations are even larger than the monthly mean errors. Though the LW errors can be reduced about 10 W/sq m after a correction for altitude difference between the site and SRB and FD grids, these errors are still higher than that for other regions. The large errors in SRB SW was mainly due to a processing mistake for elevation effect, but the errors in SRB LW was mainly due to significant errors in input data. We suggest reprocessing satellite surface radiation budget data, at least for highland areas like Tibet.
We demonstrate the practicality and effectiveness of using a Green's functions estimation approach for adjusting uncertain parameters in an Earth system model (ESM). This estimation approach has previously been applied to an intermediate-complexity climate model and to individual ESM components, e.g., ocean, sea ice, or carbon cycle components. Here, the Green's functions approach is applied to a state-of-the-art ESM that comprises a global atmosphere/land configuration of the Goddard Earth Observing System (GEOS) coupled to an ocean and sea ice configuration of the Massachusetts Institute of Technology general circulation model (MITgcm). Horizontal grid spacing is approximately 110 km for GEOS and 37–110 km for MITgcm. In addition to the reference GEOS-MITgcm simulation, we carried out a series of model sensitivity experiments, in which 20 uncertain parameters are perturbed. These “control” parameters can be used to adjust sea ice, microphysics, turbulence, radiation, and surface schemes in the coupled simulation. We defined eight observational targets: sea ice fraction, net surface shortwave radiation, downward longwave radiation, near-surface temperature, sea surface temperature, sea surface salinity, and ocean temperature and salinity at 300 m. We applied the Green's functions approach to optimize the values of the 20 control parameters so as to minimize a weighted least-squares distance between the model and the eight observational targets. The new experiment with the optimized parameters resulted in a total cost reduction of 9 % relative to a simulation that had already been adjusted using other methods. The optimized experiment attained a balanced cost reduction over most of the observational targets. We also report on results from a set of sensitivity experiments that are not used in the final optimized simulation but helped explore options and guided the optimization process. These experiments include an assessment of sensitivity to the number of control parameters and to the selection of observational targets and weights in the cost function. Based on these sensitivity experiments, we selected a specific definition for the cost function. The sensitivity experiments also revealed a decreasing overall cost as the number of control variables was increased. In summary, we recommend using the Green's functions estimation approach as an additional fine-tuning step in the model development process. The method is not a replacement for modelers' experience in choosing and adjusting sensitive model parameters. Instead, it is an additional practical and effective tool for carrying out final adjustments of uncertain ESM parameters.
Analysis of four years of earth radiation budget, cloud, and sea surface temperature data confirms that cloud parameters change dramatically when and where sea surface temperatures increase above approximately 300K. These results are based upon monthly mean values within 2.5 deg x 2.5 deg grid points over the 'warm pool' region of the western tropical Pacific. The question of whether sea surface temperatures are influenced, in turn, by the radiative effects of these clouds (Ramanathan and Collins) is less clear. Such a feedback, if it exists, is weak. The reason why clouds might have so little influence, despite large changes in their longwave and shortwave radiative effects, might be that the sea surface responds to both the longwave heating and the shortwave cooling effects of clouds, and the two effects nearly cancel. There are strong correlations between the rate of change of sea surface temperature and any of the radiation budget parameters that are highly correlated with the incident solar flux-implying that season and latitude are the critical factors determining sea surface temperatures. With the seasonal or both seasonal and latitudinal variations removed, the rate of change of sea surface temperature shows no correlation with cloud-related parameters in the western tropical Pacific.
Theoretical and satellite-based assessments of the sensitivity of broadband shortwave radiative fluxes in cloudy atmospheres to small perturbations in the cloud droplet number concentration (N) of liquid water clouds under constant water conditions are performed. Two approaches to study this sensitivity are adopted: absolute increases in N, for which the radiative response is referred to as absolute cloud susceptibility, and relative increases in N or relative cloud susceptibility. Estimating the former is more challenging as it requires an assumed value for either cloud liquid water content or geometrical thickness; both susceptibilities require an assumed relationship between the droplet volume and effective radius. Expanding upon previous susceptibility studies, present radiative calculations include the effect of AN perturbations on droplet asymmetry parameter and single-scattering albedo, in addition to extinction. Absolute cloud susceptibility has a strong nonlinear dependence on the droplet effective radius as expected, while relative cloud susceptibility is primarily dependent on optical thickness. Molecular absorption and reflecting surfaces both reduce the relative contribution of the cloud to the top-of-atmosphere (TOA) flux and therefore also reduce the TOA albedo susceptibility. Transmittance susceptibilities are negative with absolute values similar to albedo susceptibility, while atmospheric absorptance susceptibilities are about an order of magnitude smaller than albedo susceptibilities and can be either positive or negative. Observation-based susceptibility calculations are derived from MODIS pixel-level retrievals of liquid water cloud optical thickness, effective radius, and cloud top temperature; two data granule examples are shown. Susceptibility quantifies the aerosol indirect effect sensitivity in a way that can be easily computed from model fields. As such, susceptibilities derived from MODIS observations provide a higher-order test of model cloud properties used for indirect effect studies. MODIS-derived global distributions of cloud susceptibility and radiative forcing calculations are presented in a companion paper.
The first paper was Multidecadal changes in near-global cloud cover and estimated cloud cover radiative forcing, by J. R. Norris (2005, J. Geophys. Res. - Atmos., 110, D08206, doi: lO.l029/2004JD005600). This study examined variability in zonal mean surface-observed upper-level (combined midlevel and high-level) and low-level cloud cover over land during 1971-1 996 and over ocean during 1952-1997. These data were averaged from individual synoptic reports in the Extended Edited Cloud Report Archive (EECRA). Although substantial interdecadal variability is present in the time series, long-term decreases in upper-level cloud cover occur over land and ocean at low and middle latitudes in both hemispheres. Near-global upper-level cloud cover declined by 1.5%-sky-cover over land between 1971 and 1996 and by 1.3%-sky-cover over ocean between 1952 and 1997. Consistency between EECRA upper-level cloud cover anomalies and those from the International Satellite Cloud Climatology Project (ISCCP) during 1984-1 997 suggests the surface-observed trends are real. The reduction in surface-observed upper-level cloud cover between the 1980s and 1990s is also consistent with the decadal increase in all-sky outgoing longwave radiation reported by the Earth Radiation Budget Satellite (EMS). Discrepancies occur between time series of EECRA and ISCCP low-level cloud cover due to identified and probable artifacts in satellite and surface cloud data. Radiative effects of surface-observed cloud cover anomalies, called "cloud cover radiative forcing (CCRF) anomalies," are estimated based on a linear relationship to climatological cloud radiative forcing per unit cloud cover. Zonal mean estimated longwave CCRF has decreased over most of the globe. Estimated shortwave CCRF has become slightly stronger over northern midlatitude oceans and slightly weaker over northern midlatitude land areas. A long-term decline in the magnitude of estimated shortwave CCRF occurs over low-latitude land and ocean, but comparison with EMS all-sky reflected shortwave radiation during 1985-1997 suggests this decrease may be underestimated.
Because of the importance and complexity of dynamical, microphysical, and radiative processes taking place near cloud edges, the transition zone between clouds and cloud free air has been the subject of intense research both in the ASR program and in the wider community. One challenge in this research is that the one-dimensional (1D) radiative models widely used in both remote sensing and dynamical simulations become less accurate near cloud edges: The large horizontal gradients in particle concentrations imply that accurate radiative calculations need to consider multi-dimensional radiative interactions among areas that have widely different optical properties. This study examines the way the importance of multidimensional shortwave radiative interactions changes as we approach cloud edges. For this, the study relies on radiative simulations performed for a multiyear dataset of clouds observed over the NSA, SGP, and TWP sites. This dataset is based on Microbase cloud profiles as well as wind measurements and ARM cloud classification products. The study analyzes the way the difference between 1D and 2D simulation results increases near cloud edges. It considers both monochromatic radiances and broadband radiative heating, and it also examines the influence of factors such as cloud type and height, and solar elevation. The results provide insights into the workings of radiative processes and may help better interpret radiance measurements and better estimate the radiative impacts of this critical region.
Marine boundary layer clouds, including the transition from stratocumulus to cumulus, are poorly represented in numerical weather prediction and general circulation models. Further uncertainties in the cloud structure arise in the presence of biomass burning carbonaceous aerosol, as is the case over the southeast Atlantic Ocean where biomass burning aerosol is transported from the African continent. As the aerosol plume progresses across the southeast Atlantic Ocean, radiative heating within the aerosol layer has the potential to alter the thermodynamic environment and therefore the cloud structure; however, limited work has been done to quantify this along the trajectory of the aerosol plume in the region. The deployment of the First Atmospheric Radiation Measurement Mobile Facility in support of the Layered Atlantic Smoke Interactions with Clouds field campaign provided a unique opportunity to collect observations of cloud and aerosol properties during two consecutive biomass burning seasons during July through October of 2016 and 2017 over Ascension Island (7.96 °S, 14.35 °W). Using observed profiles of temperature, humidity, and clouds from the field campaign, alongside aerosol optical properties from MERRA-2 as input for the Rapid Radiation Transfer Model, profiles of the radiative heating rate due to aerosols and clouds were computed. Radiative heating is also assessed across the southeast Atlantic Ocean using an ensemble of back trajectories from the Hybrid Single Particle Lagrangian Integrated Trajectory Model. Idealized experiments using the Rapid Radiation Transfer Model with and without aerosols and a range of values for the single scattering albedo demonstrate that shortwave heating within the aerosol layer above Ascension Island can locally range between 2 and 8 K per day depending on the aerosol optical properties, though impacts of the aerosol can be felt elsewhere in the atmospheric column. When considered under clear conditions, the aerosol has a cooling effect at the TOA and based on the observed cloud properties at Ascension Island, the cloud albedo is not large enough to overcome this. Shortwave radiative heating due to biomass burning aerosol is not balanced by additional longwave cooling, and the net radiative impact results in a stabilization of the lower troposphere. However, these results are extremely sensitive to the single scattering albedo assumptions in models.
Abstract Multiple Earth system models (ESMs) discretize surface albedo into two semi‐broadbands comprising the UV/visible and near‐infrared (NIR) wavelengths. Here, we use an offline single‐column radiative transfer model to investigate the radiative effects of spectrally resolving the surface albedo. We use the Snow, Ice, and Aerosol Radiative model, extended to simulate liquid water, to calculate snow, ice, and liquid water albedo. We flux‐weight the hyperspectral albedo into the coarser spectral bands used by the atmospheric shortwave radiative transfer model. We establish representative atmospheric profiles for the three surface types and compare their shortwave fluxes and atmospheric warming rates with the spectrally resolved albedo to those calculated with the semi‐broadband approximation. Spectrally resolved surface albedo over snow and ice reduces atmospheric warming by darkening the albedo of NIR bands, correcting the too‐strong surface absorption in visible bands, and too‐weak surface absorption in shortwave infrared bands caused by the semi‐broadband approximation. We explore the effects on surface and atmospheric warming rates of varying solar zenith angle, cloud cover, relative humidity, and snow grain/air bubble radii. The semi‐broadband albedo biases can exceed 10% and 2% for the surface and atmospheric net flux respectively, being particularly strong under conditions which alter the distributions of surface insolation (i.e., cloud cover or increased atmospheric water vapor). These results show that transmitting a higher resolution spectral radiation field between the atmosphere and surface reduces biases in surface absorption and atmospheric heating present in ESMs that currently use the semi‐broadband approximation.
The Clouds and Earth’s Radiant Energy System (CERES) projects provides satellite-based observations of the radiative fluxes and clouds systems. CERES climate quality data products typically take several months of calibration and validation before release to the public. An alternative data product, Fast Longwave and Shortwave radiative Flux (FLASHFlux), was created to provide data to the applied sciences and educational users. FLASHFlux provides Top-of-Atmosphere radiative fluxes, Clouds properties, and parameterized surface radiative fluxes within four days for footprint (Level 2) data. We investigate the use of Artificial Neural Network (ANN) using MODerate resolution Imaging Spectroradiometer (MODIS) derived clouds properties and meteorology from the Global Assimilation and Meteorology Office (GMAO) scaled to the CERES footprint from the CERES Clouds Radiative Swath (CRS) data product to compute surface radiative fluxes. We test ANN produce fluxes against surface fluxes produced from the Fu-Liou model used in CRS and the Langley Parameterized Shortwave Algorithm (LPSA) and Langley Parameterized Longwave Algorithm (LPLA) used in FLASHFlux. We also validated each model with ground-based observations. Furthermore, we investigate Leave-One-Feature-Out Importance (LOFO) to evaluate the significance of each feature in our training and provide insight for future models. Advances in machine learning, along with increases in computational capabilities and available data allow us to estimate effects of unresolved processes in our climate without direct modeling. This work evaluates the ability to create accurate data-driven models to supplement or replace current models that estimate surface radiative fluxes.
To study Arctic stratus cloud properties and their effect on the surface radiation balance during the spring transition season, analyses are performed using data taken during three cloudy and two clear days in May 1998 as part of the First ISCCP Regional Experiment (FIRE) Arctic Cloud Experiment (ACE). Radiative transfer models are used in conjunction with surface- and satellite-based measurements to retrieve the layer-averaged microphysical and shortwave radiative properties. The surface-retrieved cloud properties in Cases 1 and 2 agree well with the in situ and satellite retrievals. Discrepancies in Case 3 are due to spatial mismatches between the aircraft and the surface measurements in a highly variable cloud field. Also, the vertical structure in the cloud layer is not fully characterized by the aircraft measurements. Satellite data are critical for understanding some of the observed discrepancies. The satellite-derived particle sizes agree well with the coincident surface retrievals and with the aircraft data when they were collocated. Optical depths derived from visible-channel data over snow backgrounds were overestimated in all three cases, suggesting that methods currently used in satellite cloud climatologies derive optical depths that are too large. Use of a near-infrared channel with a solar infrared channel to simultaneously derive optical depth and particle size appears to alleviate this overestimation problem. Further study of the optical depth retrieval is needed. The surface-based radiometer data reveal that the Arctic stratus clouds produce a net warming of 20 W m(exp -2) in the surface layer during the transition season suggesting that these clouds may accelerate the spring time melting of the ice pack. This surface warming contrasts with the net cooling at the top of the atmosphere (TOA) during the same period. All analysis of the complete FIRE ACE data sets will be valuable for understanding the role of clouds during the entire melting and refreezing process that occurs annually in the Arctic.
Shallow cloud decks residing in or near the boundary layer cover a large fraction of the Southern Ocean (SO) and play a major role in determining the amount of shortwave radiation reflected back to space from this region. In this article, we examine the macrophysical characteristics and thermodynamic phase of low clouds (tops <3 km) and precipitation using ground-based ceilometer, depolarization lidar and vertically-pointing W-band radar measurements collected during the Macquarie Island Cloud and Radiation Experiment (MICRE) from April 2016 to March 2017. During MICRE, low clouds occurred ~65% of the time on average (slightly more often in austral winter than summer). About 2/3 of low clouds were cold-topped (temperatures ≤0°C). These were thicker and had higher bases on average than warm-topped clouds. 83%–88% of cold-topped low clouds were liquid phase at cloud base (depending on the season). The majority of low clouds had precipitation in the vertical range 150–250 m below cloud base, a significant fraction of which did not reach the surface. Phase characterization is limited to the period between April 2016 and November 2016. Small-particle (low-radar-reflectivity) precipitation (which dominates precipitation occurrence) was mostly liquid below-cloud, while large-particle precipitation (which dominates total accumulation) was predominantly mixed/ambiguous or ice phase. Approximately 40% of cold-topped clouds had mixed/ambiguous or ice phase precipitation below (with predominantly liquid phase cloud droplets at cloud base). Below-cloud precipitation with radar reflectivity factors below about -10 dBZ were predominantly liquid, while reflectivity factors above about 0 dBZ were predominantly ice.
Modeled differences in predicted effects are increasingly used to help quantify the uncertainty of these effects. Here, we examine modeled differences in the aerosol indirect effect in a series of experiments that help to quantify how and why model-predicted aerosol indirect forcing varies between models. The experiments start with an experiment in which aerosol concentrations, the parameterization of droplet concentrations and the autoconversion scheme are all specified and end with an experiment that examines the predicted aerosol indirect forcing when only aerosol sources are specified. Although there are large differences in the predicted liquid water path among the models, the predicted aerosol first indirect effect for the first experiment is rather similar, about -0.6 W/sq m to -0.7 W/sq m. Changes to the autoconversion scheme can lead to large changes in the liquid water path of the models and to the response of the liquid water path to changes in aerosols. Adding an autoconversion scheme that depends on the droplet concentration caused a larger (negative) change in net outgoing shortwave radiation compared to the 1st indirect effect, and the increase varied from only 22% to more than a factor of three. The change in net shortwave forcing in the models due to varying the autoconversion scheme depends on the liquid water content of the clouds as well as their predicted droplet concentrations, and both increases and decreases in the net shortwave forcing can occur when autoconversion schemes are changed. The parameterization of cloud fraction within models is not sensitive to the aerosol concentration, and, therefore, the response of the modeled cloud fraction within the present models appears to be smaller than that which would be associated with model "noise". The prediction of aerosol concentrations, given a fixed set of sources, leads to some of the largest differences in the predicted aerosol indirect radiative forcing among the models, with values of cloud forcing ranging from -0.3 W/sq m to -1.4W/sq m. Thus, this aspect of modeling requires significant improvement in order to improve the prediction of aerosol indirect effects.
Here we improved the treatments of convective cloud microphysics in the NCAR Community Atmosphere Model version 5.3 (CAM5.3) by 1) implementing new terminal velocity parameterizations for convective ice and snow particles, 2) adding graupel microphysics, 3) considering convective snow detrainment, and 4) enhancing rain initiation and generation rate in warm clouds. Furthermore, we evaluated the impacts of improved microphysics on simulated global climate, focusing on simulated cloud radiative forcing, graupel microphysics, convective cloud ice amount, and tropical precipitation. Compared to CAM5.3 with the default convective microphysics, the too-strong cloud shortwave radiative forcing due primarily to excessive convective cloud liquid is largely alleviated over the tropics and midlatitudes after rain initiation and generation rate is enhanced, in better agreement with the CERES-EBAF estimates. Geographic distributions of graupel occurrence are reasonably simulated over continents; whereas the graupel occurrence remains highly uncertain over the oceanic storm-track regions. When evaluated against the CloudSat–CALIPSO estimates, the overestimation of convective ice mass is alleviated with the improved convective ice microphysics, among which adding graupel microphysics and the accompanying increase in hydrometeor fall speed play the most important role. The probability distribution function (PDF) of rainfall intensity is sensitive to warm rain processes in convective clouds, and enhancement in warm rain production shifts the PDF toward heavier precipitation, which agrees better with the TRMM observations. Common biases of overestimating the light rain frequency and underestimating the heavy rain frequency in GCMs are mitigated.
Aerosol optical properties inferred from the Sea-viewing Wide Field-of-view Sensor (SeaWiFS) radiance measurements are used to compute the aerosol shortwave radiative forcing using a radiative transfer model. The aerosol optical thickness at the wavelength of 865-nm is taken from the SeaWIFS archive. It is found that the nominal optical thickness over oceans ranges from 0.1 to 0.2. Using a maritime aerosol model and the radiances measured at the various SeaWiFS channels, the Angstrom exponent is determined to be 0.2174, the single-scattering albedo to be 0.995, and the asymmetry factor to be 0.786. The radiative transfer model has eight bands in the visible and ultraviolet spectral regions and three bands in the near infrared. It includes the absorption due to aerosols, water vapor, carbon dioxide, and oxygen, and the scattering due to aerosols and gases (Rayleigh scattering). The radiative forcing is computed over global oceans for four months (January, April, July, and October, 1998) to represent four seasons. It is found that the aerosol radiative forcing is large and changes significantly with seasons near the continents with large-scale forest fires and desert dust. Averaged over oceans and the four months, the aerosol radiative forcing is approximately 7 W/sq m at the top of the atmosphere. This large radiative forcing is expected to have a significant cooling effect on the Earth's climate as implied from simulations of a number of general circulation models.
Abstract. Mountain hydrology is controlled by interacting processes extending from the atmosphere through the bedrock. Integrated process models (IPMs), one of the main tools needed to interpret observations and refine conceptual models of the mountainous water cycle, require meteorological forcing that simulates the atmospheric process to predict hydroclimate then subsequently impacts surface–subsurface hydrology. Complex terrain and extreme spatial heterogeneity in mountainous environments drive uncertainty in several key considerations in IPM configurations and require further quantification and sensitivity analyses. Here, we present an IPM using the Weather Research and Forecasting (WRF) model which forces an integrated hydrologic model, ParFlow-CLM, implemented over a domain centered over the East River watershed (ERW), located in the Upper Colorado River basin (UCRB). The ERW is a heavily instrumented 300 km2 region in the headwaters of the UCRB near Crested Butte, CO, with a growing atmosphere-through-bedrock observation network. Through a series of experiments in the water year 2019 (WY19), we use four meteorological forcings derived from commonly used reanalysis datasets, three subgrid-scale physics scheme configurations in WRF, and two terrain shading options within WRF to test the relative importance of these experimental design choices for key hydrometeorological metrics including precipitation and snowpack, as well as evapotranspiration, groundwater storage, and discharge simulated by the ParFlow-CLM. Our hypothesis is that uncertainty from synoptic-scale forcings produces a much larger spread in surface–subsurface hydrologic fields than subgrid-scale physics scheme choice. Results reveal that the WRF subgrid-scale physics configuration leads to larger spatiotemporal variance in simulated hydrometeorological conditions, whereas variance across meteorological forcing with common subgrid-scale physics configurations is more spatiotemporally constrained. Despite reasonably simulating precipitation, a delay in simulated discharge peak is due to a systematic cold bias across WRF simulations, suggesting the need for bias correction. Discharge shows greater variance in response to the WRF simulations across subgrid-scale physics schemes (26 %) rather than meteorological forcing (6 %). The topographic radiation option has minor effects on the watershed-average hydrometeorological processes but adds profound spatial heterogeneity to local energy budgets (±30 W m−2 in shortwave radiation and 1 K air temperature differences in late summer). This is the first presentation of sensitivity analyses that provide support to help guide the scientific community to develop observational constraints on atmosphere-through-bedrock processes and their interactions.
The main objective of the International Satellite Land Surface Climatology Project (ISLSCP) has been stated as 'the development of techniques that may be applied to satellite observations of the radiation reflected and emitted from the Earth to yield quantitative information concerning land surface climatological conditions'. The major field study, FIFE (the First ISLSCP Field Experiment), was conducted in 1987-89 to accomplish this objective. Four intensive field campaigns (IFC's) were carried out in 1987 and one in 1989. Factors contributing to observed reflected radiation from the FIFE site must be understood before the radiation observed by satellites can be used to quantify surface processes. Our last report (Walter-Shea et al., 1992b) focused on slope effects on incoming and outgoing shortwave radiation and net radiation from data collected in 1989. We report here on the final analysis of the slope data as well as results from thermal radiation studies conducted during the FIFE experiment. The specific areas reported are the following: (1) analysis of slope effects on measured reflectance values and estimates of surface albedo; (2) using remotely-measured surface temperatures as a means of estimating sensible heat flux from the Konza Prairie; (3) extracting canopy temperatures from remotely-measured composite surface temperatures; (4) modeling the measured composite temperature of partially vegetated surfaces; and (5) estimating gap distribution in partially vegetated surfaces from reflectance measurements.
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.