skip to main content

Title: Spatially similar surface energy flux perturbations due to greenhouse gases and aerosols

Despite distinct geographic distributions of top-of-the-atmosphere radiative forcing, anthropogenic greenhouse gases and aerosols have been found to produce similar patterns of climate response in atmosphere-and-ocean coupled climate model simulations. Understanding surface energy flux changes, a crucial pathway by which atmospheric forcing is communicated to the ocean, is a vital bridge to explaining the similar full atmosphere-and-ocean responses to these disparate forcings. Here we analyze the fast, atmosphere-driven change in surface energy flux caused by present-day greenhouse gases vs aerosols to elucidate its role in shaping the subsequent slow, coupled response. We find that the surface energy flux response patterns achieve roughly two-thirds of the anti-correlation seen in the fully coupled response, driven by Rossby waves excited by symmetric changes to the land–sea contrast. Our results suggest that atmosphere and land surface processes are capable of achieving substantial within-hemisphere homogenization in the climate response to disparate forcers on fast, societally-relevant timescales.

; ; ;
Publication Date:
Journal Name:
Nature Communications
Nature Publishing Group
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract. For the radiative impact of individual climate forcings,most previous studies focused on the global mean values at the top of theatmosphere (TOA), and less attention has been paid to surface processes,especially for black carbon (BC) aerosols. In this study, the surface radiativeresponses to five different forcing agents were analyzed by using idealizedmodel simulations. Our analyses reveal that for greenhouse gases, solarirradiance, and scattering aerosols, the surface temperature changes aremainly dictated by the changes of surface radiative heating, but for BC,surface energy redistribution between different components plays a morecrucial role. Globally, when a unit BC forcing is imposed at TOA, the netshortwave radiation at the surface decreases by -5.87±0.67 W m−2 (W m−2)−1 (averaged over global land without Antarctica), which ispartially offset by increased downward longwave radiation (2.32±0.38 W m−2 (W m−2)−1 from the warmer atmosphere, causing a netdecrease in the incoming downward surface radiation of -3.56±0.60 W m−2 (W m−2)−1. Despite a reduction in the downward radiationenergy, the surface air temperature still increases by 0.25±0.08 Kbecause of less efficient energy dissipation, manifested by reduced surfacesensible (-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 inBowen ratio (-0.20±0.07 (W m−2)−1). Such reductions of turbulentfluxes can be largely explained by enhanced air stability (0.07±0.02 K (W m−2)−1), measured as the difference of the potential temperaturebetween 925 hPa and surface,more »and reduced surface wind speed (-0.05±0.01 m s−1 (W m−2)−1). The enhanced stability is due to the fasteratmospheric warming relative to the surface, whereas the reduced wind speedcan be partially explained by enhanced stability and reduced Equator-to-poleatmospheric temperature gradient. These rapid adjustments under BC forcingoccur in the lower atmosphere and propagate downward to influence thesurface energy redistribution and thus surface temperature response, whichis not observed under greenhouse gases or scattering aerosols. Our studyprovides new insights into the impact of absorbing aerosols on surfaceenergy balance and surface temperature response.« less
  2. Abstract. Equilibrium climate sensitivity (ECS) has been directly estimated using reconstructions of past climates that are different than today's. A challenge to this approach is that temperature proxies integrate over the timescales of the fast feedback processes (e.g., changes in water vapor, snow, and clouds) that are captured in ECS as well as the slower feedback processes (e.g., changes in ice sheets and ocean circulation) that are not. A way around this issue is to treat the slow feedbacks as climate forcings and independently account for their impact on global temperature. Here we conduct a suite of Last Glacial Maximum (LGM) simulations using the Community Earth System Model version 1.2 (CESM1.2) to quantify the forcingand efficacy of land ice sheets (LISs) and greenhouse gases (GHGs) in order to estimate ECS. Our forcing and efficacy quantification adopts the effective radiative forcing (ERF) and adjustment framework and provides a complete accounting for the radiative, topographic, and dynamical impacts of LIS on surface temperatures. ERF and efficacy of LGM LIS are −3.2 W m−2 and 1.1, respectively. The larger-than-unity efficacy is caused by the temperature changes over land and the Northern Hemisphere subtropical oceans which are relatively larger than those in response to a doubling of atmospheric CO2. The subtropical sea-surface temperature (SST) response ismore »linked to LIS-induced wind changes and feedbacks in ocean–atmosphere coupling and clouds. ERF and efficacy of LGM GHG are −2.8 W m−2 and 0.9, respectively. The lower efficacy is primarily attributed to a smaller cloud feedback at colder temperatures. Our simulations further demonstrate that the direct ECS calculation using the forcing, efficacy, and temperature response in CESM1.2 overestimates the true value in the model by approximately 25 % due to the neglect of slow ocean dynamical feedback. This is supported by the greater cooling (6.8 ∘C) in a fully coupled LGM simulation than that (5.3 ∘C) in a slab ocean model simulation with ocean dynamics disabled. The majority (67 %) of the ocean dynamical feedback is attributed to dynamical changes in the Southern Ocean, where interactions between upper-ocean stratification, heat transport, and sea-ice cover are found to amplify the LGM cooling. Our study demonstrates the value of climate models in the quantification of climate forcings and the ocean dynamical feedback, which is necessary for an accurate direct ECS estimation.« less
  3. Abstract. The Arctic is warming 2 to 3 times faster than the global average, partly due to changes in short-lived climate forcers (SLCFs) including aerosols. In order to study the effects of atmospheric aerosols in this warming, recent past (1990–2014) and future (2015–2050) simulations have been carried out using the GISS-E2.1 Earth system model to study the aerosol burdens and their radiative and climate impacts over the Arctic (>60∘ N), using anthropogenic emissions from the Eclipse V6b and the Coupled Model Intercomparison Project Phase 6 (CMIP6) databases, while global annual mean greenhouse gas concentrations were prescribed and kept fixed in all simulations. Results showed that the simulations have underestimated observed surface aerosol levels, in particular black carbon (BC) and sulfate (SO42-), by more than 50 %, with the smallest biases calculated for the atmosphere-only simulations, where winds are nudged to reanalysis data. CMIP6 simulations performed slightly better in reproducing the observed surface aerosol concentrations and climate parameters, compared to the Eclipse simulations. In addition, simulations where atmosphere and ocean are fully coupled had slightly smaller biases in aerosol levels compared to atmosphere-only simulations without nudging. Arctic BC, organic aerosol (OA), and SO42- burdens decrease significantly in all simulations by 10 %–60 % following the reductionsmore »of 7 %–78 % in emission projections, with the Eclipse ensemble showing larger reductions in Arctic aerosol burdens compared to the CMIP6 ensemble. For the 2030–2050 period, the Eclipse ensemble simulated a radiative forcing due to aerosol–radiation interactions (RFARI) of -0.39±0.01 W m−2, which is −0.08 W m−2 larger than the 1990–2010 mean forcing (−0.32 W m−2), of which -0.24±0.01 W m−2 was attributed to the anthropogenic aerosols. The CMIP6 ensemble simulated a RFARI of −0.35 to −0.40 W m−2 for the same period, which is −0.01 to −0.06 W m−2 larger than the 1990–2010 mean forcing of −0.35 W m−2. The scenarios with little to no mitigation (worst-case scenarios) led to very small changes in the RFARI, while scenarios with medium to large emission mitigations led to increases in the negative RFARI, mainly due to the decrease in the positive BC forcing and the decrease in the negative SO42- forcing. The anthropogenic aerosols accounted for −0.24 to −0.26 W m−2 of the net RFARI in 2030–2050 period, in Eclipse and CMIP6 ensembles, respectively. Finally, all simulations showed an increase in the Arctic surface air temperatures throughout the simulation period. By 2050, surface air temperatures are projected to increase by 2.4 to 2.6 ∘C in the Eclipse ensemble and 1.9 to 2.6 ∘C in the CMIP6 ensemble, compared to the 1990–2010 mean. Overall, results show that even the scenarios with largest emission reductions leads to similar impact on the future Arctic surface air temperatures and sea-ice extent compared to scenarios with smaller emission reductions, implying reductions of greenhouse emissions are still necessary to mitigate climate change.« less
  4. Abstract Under increasing greenhouse gases, the Arctic warms about twice as fast as elsewhere, known as Arctic amplification (AA). AA weakens meridional temperature gradients and is hypothesized to weaken zonal wind and cause wavier circulation with stronger meridional wind ( υ ) over northern mid-to-high latitudes. Here model simulations are analyzed to examine the υ response to increased CO 2 and AA alone. Total υ changes are found to be dominated by the effect of increased CO 2 without AA, with a zonal wavenumber-4 (wavenumber-3) change pattern over the northern (southern) extratropics that generally enhances current υ and results partly from changes in zonal temperature gradients. The extratropical υ change patterns are quasi-barotropic and are more pronounced during boreal winter. The CO 2 forcing also causes baroclinic υ changes over the tropics tied to convection changes. The impact of AA on υ is mainly over the northern extratropics and is opposite to the effect of increased CO 2 but with smaller magnitude. An eastward shift (∼5° longitude) and an amplitude increase (∼1 m s −1 ) in the climatology of the northerlies over Europe caused mainly by CO 2 forcing contribute to the drying in southern Europe, while both AAmore »and CO 2 forcing enhance the climatology of the northerlies over East Asia. Over the northern mid-to-high latitudes, Arctic sea ice loss and AA enhance the land–ocean thermal contrast in winter, while increased CO 2 alone weakens it, resulting in opposite changes in zonal temperature gradients and thus υ . Different warming rates over land and ocean also contribute to the intermodel spread in υ response patterns among climate models. Significance Statement Meridional wind ( υ ) greatly contributes to thermal and moisture advection due to large meridional gradients in these fields. It is hypothesized that the enhanced Arctic warming under anthropogenic global warming could weaken meridional temperature gradients, decelerate westerly jets, and cause wavier circulation with stronger υ over northern extratropics. Using novel climate model simulations, we found that the effect of increased CO 2 without AA determines the total υ changes. AA generally weakens the climatological υ , contrary to the direct effect of increased CO 2 . The υ changes are small relative to its climatology but may have large impacts on regional climate over central Europe, East Asia, and interior North America. More research is needed to examine the mechanisms causing regional υ changes.« less
  5. Abstract. Since the middle Miocene (15 Ma, million years ago), the Earth's climate has undergone a long-term cooling trend, characterised by a reduction in ocean temperatures of up to 7–8 ∘C. The causes of this cooling are primarily thought to be due to tectonic plate movements driving changes in large-scale ocean circulation patterns, and hence heat redistribution, in conjunction with a drop in atmospheric greenhouse gas forcing (and attendant ice-sheet growth and feedback). In this study, we assess the potential to constrain the evolving patterns of global ocean circulation and cooling over the last 15 Ma by assimilating a variety of marine sediment proxy data in an Earth system model. We do this by first compiling surface and benthic ocean temperature and benthic carbon-13 (δ13C) data in a series of seven time slices spaced at approximately 2.5 Myr intervals. We then pair this with a corresponding series of tectonic and climate boundary condition reconstructions in the cGENIE (“muffin” release) Earth system model, including alternative possibilities for an open vs. closed Central American Seaway (CAS) from 10 Ma onwards. In the cGENIE model, we explore uncertainty in greenhouse gas forcing and the magnitude of North Pacific to North Atlantic salinity flux adjustment required in the modelmore »to create an Atlantic Meridional Overturning Circulation (AMOC) of a specific strength, via a series of 12 (one for each tectonic reconstruction) 2D parameter ensembles. Each ensemble member is then tested against the observed global temperature and benthic δ13C patterns. We identify that a relatively high CO2 equivalent forcing of 1120 ppm is required at 15 Ma in cGENIE to reproduce proxy temperature estimates in the model, noting that this CO2 forcing is dependent on the cGENIE model's climate sensitivity and that it incorporates the effects of all greenhouse gases. We find that reproducing the observed long-term cooling trend requires a progressively declining greenhouse gas forcing in the model. In parallel to this, the strength of the AMOC increases with time despite a reduction in the salinity of the surface North Atlantic over the cooling period, attributable to falling intensity of the hydrological cycle and to lowering polar temperatures, both caused by CO2-driven global cooling. We also find that a closed CAS from 10 Ma to present shows better agreement between benthic δ13C patterns and our particular series of model configurations and data. A final outcome of our analysis is a pronounced ca. 1.5 ‰ decline occurring in atmospheric (and ca. 1 ‰ ocean surface) δ13C that could be used to inform future δ13C-based proxy reconstructions.« less