# Quantifying climate feedbacks in polar regions

Hugues Goosse, Jennifer E Kay, Kyle C Armour, Alejandro Bodas-Salcedo, Hélène Chepfer, David Docquier, Alexandra Jonko, Paul J Kushner, Olivier Lecomte, François Massonnet, et al.

<figure><img src="image_1.png" /></figure>

## To cite this version:

Hugues Goosse, Jennifer E Kay, Kyle C Armour, Alejandro Bodas-Salcedo, Hélène Chepfer, et al. Quantifying climate feedbacks in polar regions. Nature Communications, 2018, 9, pp.1919. 10.1038/s41467-018-04173-0. hal-01801558

HAL Id: hal-01801558

https://hal.sorbonne-universite.fr/hal-01801558

Submitted on 28 May 2018

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
PERSPECTIVE

DOI: 10.1038/s41467-018-04173-0

OPEN

Quantifying climate feedbacks in polar regions

Hugues Goosse ID<fnref n="1" />, Jennifer E. Kay², Kyle C. Armour³,
Alejandro Bodas-Salcedo⁴, Helene Chepfer⁵, David Docquier ID<fnref n="1" />,
Alexandra Jonko ID⁶, Paul J. Kushner ID<fnref n="7" />, Olivier Lecomte<fnref n="1" />,
François Massonnet ID<fnref n="1" />,⁸, Hyo-Seok Park⁹, Felix Pithan ID¹⁰,
Gunilla Svensson ID¹¹ & Martin Vancoppenolle ID¹²

The concept of feedback is key in assessing whether a perturbation to a system is amplified or damped by mechanisms internal to the system. In polar regions, climate dynamics are controlled by both radiative and non-radiative interactions between the atmosphere, ocean, sea ice, ice sheets and land surfaces. Precisely quantifying polar feedbacks is required for a process-oriented evaluation of climate models, a clear understanding of the processes responsible for polar climate changes, and a reduction in uncertainty associated with model projections. This quantification can be performed using a simple and consistent approach that is valid for a wide range of feedbacks, offering the opportunity for more systematic feedback analyses and a better understanding of polar climate changes.

The climate of polar regions is highly sensitive to changes in climate forcing, but also displays large internal variability. Over recent decades, northern polar regions have warmed more than twice the global average with sea-ice decreasing trends for all months of the year, especially in late summer<fnref n="1" />⁻³. In contrast, the southern polar regions have warmed less rapidly with some regions experiencing cooling and sea ice advance and others experiencing warming and sea ice loss⁴⁻⁶.

Observed changes in polar regions result from numerous interactions involving the atmo-
sphere, land surfaces, ocean and sea ice. Due to the complexity of the underlying processes, we
do not fully understand them. Advancing scientific understanding in polar regions is particularly
challenging due to a short and incomplete observational record⁵, ⁶, large internal climate
variability⁶⁻⁹ and the large biases of climate models in these regions¹⁰.

The feedback framework¹¹, ¹² offers a standard method to analyze such complex dynamics (Box 1). The first step is to define a simple reference system and to estimate the response of this reference system to a perturbation. In a second step, the internal dynamics processes are

<footnote n="1">¹Earth and Life Institute, Université catholique de Louvain, Louvain-la-Neuve B 1348, Belgium. ²Department of Atmospheric and Oceanic Sciences, and Cooperative Institute for Research in Environmental Science, University of Colorado, Boulder, CO 80309, USA. ³School of Oceanography and Department of Atmospheric Sciences, University of Washington, Seattle, WA 98105, USA. ⁴Met Office Hadley Centre, Exeter EX1 3PB, UK. ⁵Sorbonne Université, UPMC Paris 6, LMD-IPSL, CNRS, Paris 75005, France. ⁶Earth and Environmental Sciences Division, Los Alamos National Laboratory, New Mexico, NM 87545, USA.</footnote>

<footnote n="7">⁷ Department of Physics, University of Toronto, Toronto M5S 1A7, Canada. ⁸ Earth Sciences Department, Barcelona Supercomputing Center, Barcelona 08034, Spain. ⁹ Korea Institute of Geoscience and Mineral Resources, Daejeon 34132, South Korea. ¹⁰ Alfred Wegener Institute, Helmholtz Centre for Polar and Marine Research, Bremerhaven D-27570, Germany. ¹¹ Department of Meteorology and Bolin Center for Climate Research, Stockholm University, Stockholm 10691, Sweden. ¹² Sorbonne Université, CNRS, IRD, MNHN, LOCEAN-IPSL, Paris 75252, France. Correspondence and requests for materials should be addressed to H.G. (email: hugues.goosse@uclouvain.be)</footnote>
Box 1 | the standard radiative feedback framework

The radiative feedback framework is based on the analysis of changes to the energy balance at the top of the atmosphere (TOA) caused by a perturbation. An initial perturbation to TOA radiation, F (in W m⁻²), is termed the 'radiative forcing' and is due, for instance, to a change in the atmospheric concentration of carbon dioxide (CO₂) or in solar irradiance. Consider as an example the response to a positive radiative forcing resulting from an increase in greenhouse gas concentrations. This will initially lead to a decrease in outgoing longwave radiation, resulting in a TOA radiative imbalance and accumulation of energy within the climate system¹¹⁷. This in turn will trigger changes in the climate, in particular a temperature increase that leads to larger emissions of infrared radiation by the Earth. Ultimately, those larger emissions will compensate for the additional energy input due to the forcing. After some time, the climate system will come into a new equilibrium characterized by higher temperatures than before the perturbation was applied.

When studying the energy budget of the whole Earth, it is convenient to assume that the modification of the radiative fluxes emitted by the Earth is proportional to changes in global mean surface temperature T_s (in K). The imbalance of the energy budget at the TOA averaged over the whole Earth at any time (ΔR, in W m⁻²) is then expressed as

$$ \Delta R = F + \lambda \Delta T_s \quad (1) $$

where λ is the net climate feedback parameter (W m⁻² K⁻¹), which is a key characteristic of the climate system response, and ΔT_s is the global surface air temperature change following the perturbation. λ is negative for a stable climate and a larger absolute value corresponds to a less sensitive climate characterized by a smaller temperature change for a specific forcing. At equilibrium, when by definition the heat budget is balanced at the TOA (i.e., ΔR = 0), the surface temperature change in response to the perturbation is simply ΔT_s = -F/λ. The equilibrium climate sensitivity, estimated as the global mean temperature change in response to a doubling of the CO₂ concentration in the atmosphere, which corresponds to a radiative forcing F of roughly 3.7 W m⁻², is thus equal to -3.7/λ. The transient imbalance at the TOA corresponds to heat storage, which is mainly accounted for by the ocean, so the term ΔR is often approximated by the ocean heat uptake⁶⁶.

The net climate feedback parameter λ can be separated into contributions from changes in surface albedo, clouds, water vapor and temperature, referred to as the feedback variables. The feedback related to temperature is itself the sum of a contribution from vertically homogenous warming or cooling (black-body response or Planck feedback, denoted by λ₀) and one from changes in vertical temperature gradient (the lapse rate feedback). For each process, a specific feedback parameter λᵢ can be computed. Their sum approximatively gives back the net climate feedback parameter λ = λ₀ + Σλᵢ + ε, where ε is a small residual accounting for non-linearities. A positive value of the feedback parameter λᵢ corresponds to a positive feedback, a negative one to a negative feedback

<figure><img src="image_3.png" /></figure>

Schematic illustration of the radiative feedback framework (based on ref. ¹²). With no feedback, the reference response of the system to a radiative forcing F is considered to be ΔT₀ = -F/λ₀. A feedback, with a feedback parameter λᵢ, will induce a change in the radiative balance λᵢΔT_s that will reinforce or dampen the effect of the radiative forcing, leading to a response of the system of ΔT_s = -(F + λᵢΔT_s)/λ₀.

represented as feedbacks that are triggered by the initial response and amplify (positive feedback) or dampen (negative feedback) it.

In climate dynamics, the classical radiative feedback framework links global surface temperature changes to perturbations of Earth's top-of-atmosphere energy budget¹²⁻¹⁶ and serves as a critical tool for quantifying climate response to greenhouse gas forcing. For instance, the magnitude of radiative feedbacks can be directly related to equilibrium climate sensitivity, commonly defined as the equilibrium global mean temperature change in response to a doubling of the CO₂ concentration in the atmosphere¹²⁻¹⁶.

As well as radiative feedbacks, other types of feedbacks affect polar regions (Fig. 1). While analyses of radiative feedbacks in polar regions have provided clear insights into processes controlling high latitude climate change, there is much less agreement on the relative importance of non-radiative feedbacks and on the way to quantify them.

Here, we provide an overview of key radiative and non-radiative feedbacks in polar regions, how they are currently evaluated and discuss why they are important for our understanding of polar climate change. We also propose an inclusive methodology that can be applied to quantify the influence of all those feedbacks, and eventually stimulate more systematic

analyses in observational and model ensembles. Estimating the magnitude of feedbacks is essential for improving our understanding of the dynamics of polar climate and to identify the relative contribution of various processes to observed high-latitude changes. In addition, it is a powerful tool to identify the origin of model biases and to reduce the uncertainty in the response to anthropogenic forcing which is directly linked to feedbacks.

## Feedbacks in polar regions

**Radiative feedbacks.** The temperature feedback represents the changes in infrared (longwave) radiative fluxes due to changes in surface and tropospheric temperatures (Table 1). It can be decomposed into a Planck feedback due to radiation changes caused by vertically uniform warming of the surface and troposphere and a lapse rate feedback due to vertically non-uniform warming¹⁷. The negative Planck feedback is the climate system's basic response to forcing that drives the system to a new equilibrium temperature. Due to the dependence of blackbody radiation on temperature, the Planck feedback, or in other words the increase in outgoing longwave radiation per unit of local warming, is less negative in polar regions than at lower latitudes¹⁸.
<figure><img src="image_4.png" /><figcaption>**Fig. 1** A schematic of some important radiative and non-radiative feedbacks in polar regions involving the atmosphere, the ocean, sea ice and ice sheets. TOA refers to the top of the atmosphere. Solar radiation (in yellow) and Infrared Radiation (in red) represent the shortwave (solar) and longwave (infrared) radiation exchanges. A red plus sign means that the feedback is positive, a negative blue sign corresponds to a negative feedback. Both signs are present for cloud feedbacks as both positive and negative feedbacks are occurring simultaneously and the net effect is not known. The gray line on the right represents a simplified temperature profile in polar regions for the atmosphere and the ocean, the dashed line corresponding to a strong surface inversion. Oceanic and atmospheric heat transport are mentioned but without signs as the processes involved are not restricted to polar regions and it is not clear if they could be formally expressed using a closed feedback loop</figcaption></figure>

While the lapse rate feedback is negative in the tropics, it is often positive in polar regions because stable stratification, especially in non-summer months, suppresses vertical mixing and warming remains largely confined to a thin near-surface layer¹⁹, ²⁰.

As the surface warms, additional water vapor amplifies the greenhouse effect and induces further warming²¹, ²². This water vapor feedback is largest in the tropics where the climatological temperature is higher and the increase in water vapor is at its maximum. In polar regions, the positive water vapor feedback is weaker than in the tropics but it still plays a relevant role in the polar response to the forcing¹⁹, 23, 24.

The surface albedo feedback is a first-order visible (shortwave) positive radiative climate feedback mechanism in polar regions²⁵⁻²⁸. As the climate warms, snow and ice cover melt, exposing underlying surfaces that typically have much lower albedos. This leads to an increased absorption of shortwave radiation by the surface, and as a result amplifies the initial warming. When melting, the snow covering Arctic sea ice contributes to forming melt ponds, increasing the absorption of solar radiation and amplifying the surface albedo feedback²⁹. Melt ponds do not form in the Southern Ocean as surface melting is very limited there, providing an illustration of different ways snow and ice interactions affect the surface albedo feedback²⁹, ³⁰.

Clouds influence the heat balance of the Earth by affecting the radiative fluxes in both visible and infrared bands and are involved in a variety of feedbacks¹⁴, ³¹, ³². The sign of any cloud feedback depends on the balance of shortwave cooling and longwave heating by the clouds. Cloud feedbacks are the most uncertain of all the radiative feedbacks as the cloud radiative effect depends on several factors that can be modified by the initial response to the perturbation¹⁴, ³³⁻³⁵. Among all mechanisms involved, two polar-specific cloud feedback examples are listed in Table 1: the cloud sea-ice feedback³⁶⁻³⁹ and the cloud optical depth feedback³², ³⁴, ⁴⁰. When sea ice melts and new open water is exposed, surface turbulent heat fluxes can increase humidity in the lower atmosphere and increase low-level clouds. During polar night, increasing low cloud increases downwelling longwave radiation, leading to further sea ice loss and thus to a positive feedback. Observational evidence shows that this cloud-sea ice feedback operates in non-summer months in both the Arctic³⁷, ³⁹ and the Antarctic⁴¹. The cloud optical depth feedback operates both at mid- and high-latitudes. Cloud liquid particles are smaller than cloud ice particles, and are therefore more efficient at reflecting solar radiation back to space. As the climate warms, the total amount of cloud water in mixed phase clouds increases, which increases the amount of reflected solar radiation
<table><caption>Table 1 Key radiative and non-radiative feedbacks in polar regions that are related to the atmosphere, ocean, sea ice, ice sheets and land surfaces and can be measured using a feedback factor</caption><thead><tr><th></th><th>Name</th><th>Description</th><th>Measure</th><th>Reference(s)</th></tr></thead><tbody><tr><td rowspan="10">Radiative feedbacks</td><td>Planck (-)</td><td>Higher surface and atmospheric temperatures increase outgoing longwave radiation, avoiding runaway warming</td><td>Change of TOA flux due to temperature change at constant lapse rate divided by surface temperature change</td><td>12-14,18</td></tr><tr><td>Lapse rate (+ in Arctic, close to 0 in Antarctic)</td><td>In a warmer world and at high latitudes, stable stratification conditions in the lower troposphere result in a larger warming of the lower than of the upper troposphere, leading to a smaller increase in outgoing longwave radiation compared to vertically uniform warming, and thus to further warming</td><td>Change of TOA flux due to lapse rate changes divided by surface temperature change (normalized by Planck feedback)</td><td>19,20</td></tr><tr><td>Surface albedo (+)</td><td>Melting ice and snow lowers surface albedo, leading to increased absorption of shortwave radiation and amplified warming</td><td>Change of TOA flux due to surface albedo change divided by surface temperature change (normalized by Planck feedback)</td><td>19,27,28,101</td></tr><tr><td>Water vapor (+)</td><td>In a warming climate, the amount of water vapor in the atmosphere increases, which amplifies the greenhouse effect and leads to further warming</td><td>Change of TOA flux due to water vapor change divided by surface temperature change (normalized by Planck feedback)</td><td>22-24</td></tr><tr><td>Cloud (+/-) Two examples are provided below</td><td>Warming of the atmosphere leads to changes in the amount and characteristics of clouds, modifying the radiative balance. The cloud contribution can be decomposed in several ways, two examples being given below</td><td>Change of TOA flux due to changes in cloud properties divided by surface temperature change (normalized by Planck feedback)</td><td>14,31-41</td></tr><tr><td>Example 1: Cloud-sea ice (+ in non-summer months, close to 0 in summer)</td><td>Decreased sea ice extent in non-summer months results in greater cloud cover and increased downwelling longwave radiation, leading to further sea ice loss</td><td>Change of TOA flux due to changes in cloud amount and opacity resulting from varying sea ice concentration divided by surface temperature change</td><td>36-39</td></tr><tr><td>Example 2: Cloud optical depth (-)</td><td>As the climate warms, the fraction of liquid water in mixed-phase clouds increases, resulting in higher cloud albedo, more reflection of shortwave radiation and reduced warming</td><td>Change of TOA flux due to changes in cloud optical depth divided by surface temperature change</td><td>32,34,40</td></tr><tr><td>Ice production-entrainment (-) (mostly active in Southern Ocean)</td><td>Brine rejection during sea ice formation induces an ocean mixed layer deepening that brings to the surface warmer water from deeper levels, melting a part of the ice initially formed and inhibiting further ice production.</td><td>Ratio of the sea ice melt due to the entrainment of warmer water in the mixed layer to the initial ice formation</td><td>50,51</td></tr><tr><td>Ice production-ocean heat storage (+) (mostly active in Southern Ocean)</td><td>Anomalous sea ice production induces vertical exchanges of salt, a higher stratification, storage of heat at depth and finally lower oceanic heat fluxes that favor further ice production.</td><td>Ratio of the latent heat associated to ice production to the heat content change of the ocean subsurface layer</td><td>52,53</td></tr><tr><td>Ice growth-thickness (-)</td><td>Thin sea ice grows more rapidly than thick sea ice due to its higher heat conduction, dampening the response to an initial decrease imposed by a perturbation.</td><td>Normalized difference in the thickness response to an energetic perturbation with and without thickness dependence of the ice growth rate</td><td>48,49</td></tr><tr><td>Surface mass balance-elevation (+) (mostly active in Greenland Ice Sheet)</td><td>Increased air temperature leads to ice melting, which lowers the surface elevation of the ice sheet, hence leading to ice exposure to warmer air temperatures and further ice melting.</td><td>Ratio of the additional sea level contribution due to this feedback to the sea level contribution without feedback</td><td>56,57</td></tr><tr><td>Ice shelf melting sea ice (-) (mostly active in Southern Ocean)</td><td>Ocean warming leads to ice shelf melting, which releases freshwater into the ocean and reduces vertical mixing. This results in sea ice expansion and reduced ocean warming.</td><td>Ratio of the additional change in sea ice extent caused by this feedback to the total change in extent without feedback</td><td>63,64</td></tr><tr><td>Marine ice sheet instability (+) (mostly active in West Antarctic Ice Sheet)</td><td>An initial retreat in the grounding line position of a marine ice sheet on an upward-sloping bed towards the ocean leads to increased ice discharge, ice thinning and further retreat.</td><td>Ratio of the additional sea level contribution due to this feedback to the sea level contribution without feedback</td><td>58-60</td></tr></tbody></table>

The proposed selection is illustrative rather than exhaustive. The sign in the first column indicates whether the feedback is positive or negative in polar regions
(i.e., increase the planetary albedo), acting as a negative feedback³². At the same time, the fraction of cloud water that is liquid also increases, enhancing the effect of this cloud optical depth feedback (Supplementary Fig. 1). Climate models robustly show this feedback³⁴ but in a manner that is stronger than implied by observations⁴², ⁴³ —in many cases due to excessive cloud ice in the present-day modeled climate⁴⁴⁻⁴⁶.

**Feedbacks related to sea ice and the ocean.** As the magnitude of some radiative feedbacks is modulated by processes that may appear hidden in their overall evaluation, focused analyses have been proposed to reveal elements specifically related to sea ice and ocean. For instance, the surface albedo feedback is more efficient for thin than for thick sea ice as a similar change in thickness induced by a given perturbation will lead to a larger increase in the open water area and thus a larger change in albedo. This has led to the definition of the open water formation efficiency as the percent open water formation per cm of ice melt over the melt season⁴⁷.

Some other feedbacks are not directly related to radiative processes. For example, basal sea ice growth rate is largely driven by heat conduction, which varies as the inverse of sea ice thickness: thin ice grows much faster than thick ice⁴⁸. At the same time, sea ice melt rate is nearly independent of ice thickness. This leads to the negative ice growth-thickness feedback. When a positive radiative perturbation is applied to the sea ice surface energy balance leading to an initial ice thinning, ice formation in winter is enhanced so that the ice adjusts its thickness to match the new growth rates to the new melting rates⁴⁹, resulting in a new equilibrium for the sea ice thickness. As the thermal insulation power of snow is even more efficient than that of sea ice, its influence on the ice growth-thickness feedback is important, but not fully quantified due to the complex relationship between snow depth and sea-ice thickness³⁰.

Because of the large heat fluxes at the ocean-ice interface in the Southern Ocean, feedback analyses there have often focused on the interactions between ocean and sea ice. In the ice-covered parts of the Southern Ocean, the stability of the upper water column is controlled by salinity, with the ocean temperature and salinity increasing with depth below the surface mixed layer (the ocean layer that has a nearly homogeneous density). In winter, when brine is released by sea ice formation, surface mixed layer density increases, inducing a mixed layer deepening and the entrainment of relatively warm water into the surface oceanic layer. This warm water reduces ice formation and can even melt ice, which partly compensates for the initial ice formation, leading to a negative ice production-entrainment feedback⁵⁰, ⁵¹

The ice production-entrainment feedback is acting at the seasonal scale but the magnitude of the ice-ocean flux can also be modulated at inter-annual to decadal timescales, leading to the positive ice production heat-storage feedback⁵², ⁵³. If ice production is very large during a particular year, the mixed layer will be deeper and the salt released by ice formation will be distributed over a larger depth range. In summer, the mixing is weaker and the freshwater input due to ice melting will be spread over a shallower layer, leading to a net downward vertical salt transport. This can restrain the vertical exchanges in the water column the subsequent winter, leading to less heat transfer to the surface and more heat storage at depth. Finally, the weaker heat flux at the ice-ocean interface would favor additional ice formation, leading to a positive feedback amplifying the initial perturbation. The heat storage at depth can also be reinforced by a net freshwater input at the surface (due for instance to a net transport of sea ice to the region) that further stabilizes the water column.

**Feedbacks related to land surfaces and ice sheets.** At low and mid-latitudes, a drying of the soils in response to an initial temperature rise can amplify the warming as the evaporation transpiration that normally cools the surface is reduced⁵⁴. This positive soil moisture-temperature feedback is less active at high latitudes but, despite the low temperatures, evapotranspiration strongly contributes to moderate the summer warming over land⁵⁵.

Ice sheets, glaciers and snow cover over land provide key components to the surface albedo feedback²⁵⁻²⁸. They also give rise to a number of specific feedback mechanisms. Three important ones are discussed here. Compared to the feedbacks mainly involving the atmosphere and sea ice, ice sheets generally, but not exclusively, play a role on longer time scales. In the positive surface mass balance-elevation feedback, increased air temperature leads to ice melting, which lowers the surface elevation of the ice sheet, exposing the ice to warmer air temperatures and thus further melting⁵⁶, ⁵⁷. This positive feedback is mostly relevant for the Greenland ice sheet where surface melting is substantial, while currently the Antarctic ice sheet hardly experiences it due to very low surface temperature.

The marine ice sheet instability has the potential to destabilize large ice sheet regions⁵⁸⁻⁶⁰. The stability of a marine ice sheet such as the West Antarctic Ice Sheet is determined by the position of the grounding line, i.e., the boundary between the grounded ice sheet and the floating ice shelf. If it is located on a bedrock sloping downward toward the interior of the ice sheet, an initial retreat of the grounding line, for instance due to basal ice melting, leads to an increase in ice discharge, which results in a further retreat of the grounding line inland until a new stable position is reached. Rapid changes in ice sheets may also be linked to the acceleration of the ice transport due to basal lubrication caused by meltwater penetration at the bed or to breakup of the ice shelves because of a weakening of the ice due to surface melting⁶¹, ⁶².

Another feedback mainly acting in the Southern Ocean is related to the interactions between floating ice shelves, sea ice and the surrounding ocean. A subsurface Southern Ocean warming leads to increased basal ice shelf melting, and the upper ocean layers get fresher due to the resulting cold freshwater input. This results in lower heat flux from the ocean interior to the surface, sea ice expansion and reduced ocean surface warming, providing a negative ice shelf melting sea ice feedback⁶³, ⁶⁴.

**Non-local feedbacks and feedbacks involving other components of the climate system.** While this Perspective focuses on feedbacks that act through physical processes in polar regions (Table 1, Fig. 1), we should mention that many other feedback processes exist, some of which involve biological processes and biogeochemical cycles⁶⁵⁻⁶⁸. One example is the bio-optical feedback, which occurs when climate warming and sea ice retreat in the Arctic Ocean lead to intense phytoplankton blooms⁶⁹. These blooms trap the penetrating solar heat flux at the ocean surface, which increases sea surface temperature. As a result, sea ice concentration decreases, which leads to enhanced absorption of solar energy into the ocean and further warming of the Arctic⁷⁰.

The response to a perturbation also implies a redistribution of the energy between different latitudes. First, the warming of the tropics under greenhouse gas forcing leads to enhanced poleward energy transport by the atmospheric circulation to higher latitudes, contributing to warming there⁷¹⁻⁷⁵. This indicates a coupling between radiative feedbacks and atmospheric heat transport⁷⁴, ⁷⁶, ⁷⁷. Moreover, radiative feedbacks at low latitudes may influence polar warming through their effect on poleward energy transport, while changes in polar regions may affect dynamics in the lower latitudes⁷⁸. This is an area of ongoing
<figure><img src="image_5.png" /><figcaption>**Fig. 2** Nonlinearity in the surface albedo feedback factor for three consecutive doublings of CO₂. The feedback factor, defined as the ratio of the magnitude of the albedo feedback on the Planck feedback, is calculated using the radiative kernel technique⁸⁵ and zonal averages are plotted for three consecutive doublings of CO₂ concentrations in CCSM3. The global average feedback factor decreases⁸⁷ from 0.097 for 2xCO₂-CNTL to 0.053 for 8xCO₂-4xCO₂.</figcaption></figure>

research, and it is not yet clear if the response of the system has a direct impact on the original perturbation itself such that a closed feedback loop can be identified.

Second, the ocean heat transport has been found to strongly shape polar climate change, with increased poleward heat transport into the Arctic⁷⁹⁻⁸¹ and decreased poleward heat transport into the Southern Ocean⁴, ⁷⁹ under global warming. Here too, it is unclear whether these changes can be represented in terms of a closed feedback loop (e.g., sea ice thinning enhancing ocean heat transport into the Arctic⁸¹) or should be classified as important drivers of polar climate change that cannot be expressed within a feedback framework.

## Quantitative evaluation of feedbacks in polar regions

**Radiative feedbacks.** The feedback parameters provide classical measures of the magnitude of the radiative feedbacks (Box 1). They are defined as the change in radiative fluxes due to the impact of a change in surface temperature upon the variable of interest (e.g., surface albedo, water vapor amount, cloud cover also referred to as climate variable) and are quantified in W m⁻² K⁻¹. The net climate feedback parameter λ, which is equal to the sum of all the parameters for the individual feedbacks, can be estimated by measuring all the terms of the equation describing the global mean radiative balance (Eq. (1) in Box 1) or by regressing the change in radiative flux at the top of the atmosphere (TOA) against the global mean surface temperature change⁸². It is somewhat more complex to evaluate specific feedback parameters λᵢ, as this requires isolating the impact of each feedback variable on the Earth's energy balance³¹, ³³, ⁸³⁻⁸⁷.

Since TOA fluxes determine the total energy budget of the Earth's climate system, they are a natural reference point for computing climate feedbacks at a global scale. They are also closely connected to surface temperature change in the Tropics, where deep convection leads to a vertically well-mixed atmosphere. In the Arctic, where deep vertical mixing is suppressed by strong static stability in the troposphere, computing feedback parameters based on surface fluxes can lead to important additional insights¹⁹, ²³, ⁵⁵, ⁷⁴. For example, a change in clouds that raises atmospheric emissivity in the Arctic inversion layer can lead to increases in both upwelling and downwelling longwave radiation, and thus lead to energy loss and a negative cloud feedback at TOA but energy gain and a positive cloud feedback at the surface¹⁹, ⁸⁸.

Individual feedback parameters defined at the surface or TOA can be diagnosed using several different methods, including partial radiative perturbations³¹, the less computationally

expensive approximate partial radiative perturbations⁸³, and the even more idealized radiative kernel technique⁸⁵, ⁸⁶. Using this now widely used method, changes in TOA radiative fluxes due to a uniform, idealized perturbation in the feedback variable are first computed using a radiative transfer model to obtain the so-called kernel. The kernel thus only depends on the radiative transfer algorithm and the mean state of the system. λᵢ can then be derived by multiplying the kernel by the response of the feedback variable to changes in global mean surface temperature.

In parallel to feedback parameters, other expressions can sometimes be easier to interpret or be more convenient. One option is to diagnose the temperature change that can be attributed to each feedback explicitly, known as a warming contribution (see the methods). It is also instructive to compare the temperature changes due to a particular feedback to changes of a reference system in which the feedbacks of interest are inactive. In the radiative feedback framework, the reference system is traditionally chosen as the Planck response. The feedback factor γᵢ is then defined as the ratio of each feedback parameter to (minus) the Planck feedback λ₀: γᵢ = λᵢ/−λ₀. An advantage of this approach is that the feedback factor γᵢ is dimensionless because it is expressed relative to the reference system. It can then be used to compare the impact of very different processes, bearing in mind that its specific value depends on the reference system chosen¹²(for more details see the methods).

In addition to the approaches focused on the top of the atmosphere or the surface, it is possible to analyze the origin of three-dimensional temperature changes such as in the climate feedback response analysis method⁸⁹ (CFRAM). It has also been proposed to decompose the feedbacks in ways that differ⁵⁵, ⁹⁰ from the traditional one described in Box 1. Each methodology is adapted to a special purpose but also has its own limitations. For instance, a three-dimensional analysis can highlight the processes that are at the origin of the changes at various level in the atmosphere, but it may require model outputs that are not routinely saved by climate modeling centers. Finally, applying different methods leads to different definitions of feedbacks and ultimately differing quantitative assessment of feedback strengths.

**Limitations of the linear approach.** The standard radiative feedback framework assumes that the response of the system can be expressed as a linear function of the surface temperature. It is a very useful approximation but some processes cannot be expressed in terms of functions of single variables and the radiative feedback framework has to be adjusted to capture changes in the system not directly related to surface temperature⁹¹, ⁹². Moreover, assuming linearity in feedbacks fails in many cases, as can be expected for a system as complex as the Earth's climate. For example, the magnitude of the climate feedback parameter λ generally decreases with time in climate models after a rise in atmospheric CO₂ concentration, corresponding to increasing climate sensitivity as equilibrium is approached³², ⁷⁶, ⁹². λ may also depend on the magnitude of the perturbation⁸⁷, ⁹³⁻⁹⁵.

The non-linearity of the feedbacks can be described in different ways. A simple definition will be used here: the feedback is non-linear if the feedback factor γ is not constant. Non-linearities can be caused by several processes. The strength of the feedback can be a function of the state of the system. This state dependence can often be expressed as a time dependence when the state changes with time. Furthermore, the different processes controlling the response to a perturbation may have different time scales. Their relative contribution to local and global scale feedbacks may thus evolve leading to spatially or temporally non-constant feedback factors⁷⁶, ⁹², ⁹⁶.
<figure><img src="image_6.png" /><figcaption>**Fig. 3** Contributions of each feedback and atmospheric forcing to polar amplification. **a** Arctic (60-90N) relative to tropics (30S-30N), and **b** Antarctic (60-90S) relative to Arctic (60-90N) at year 100 of abrupt CO₂ quadrupling in climate models involved in the fifth phase of the Coupled Model Intercomparison Project (CMIP5). The feedbacks shown are the lapse rate (LR), surface albedo (A), water vapor (WV), cloud (C), and latitudinal variation in the Planck response (P', local difference from its global-mean value λ₀); the additional energetic contributions shown are the CO₂ forcing (CO₂). Atmospheric heat transport convergence (ΔAHT) and ocean heat uptake (ΔOHU) (see method section). The feedbacks are expressed as warming contributions to the total temperature change</figcaption></figure>

In polar regions, the presence of different phases of water implies that many feedback parameters display a particularly strong dependence on the state of the system near the freezing point and are thus highly non-linear. For instance, phase changes play an important role in polar clouds leading to non-linearities in the cloud feedback³², ³⁴, ³⁹, ⁴⁰. Furthermore, feedbacks related to the cryosphere generally depend on the surface area covered by snow or ice. As temperatures rise, this area decreases and the feedback strength is reduced. This is illustrated⁸⁷ in Fig. 2 for the surface albedo feedback in response to three consecutive doublings of CO₂ in the Community Climate System Model version 3 (CCSM3). At many latitudes, the value of the feedback factor is smaller for the third doubling (8 × CO₂–4 × CO₂) than it is for the first (2 × CO₂–CNTL). Between 50°S and 60°S the feedback approaches zero for the third doubling, since the Southern Ocean is already ice-free at these latitudes in the 4xCO₂ climate, and no further melting can occur. On the other hand, the value of the feedback factor increases at northern high latitudes (75°N–90°N), as the sea ice edge retreats within the central Arctic at high warming.

**Non-radiative feedbacks.** The traditional radiative feedback framework has been extended to additional processes that influence the energy balance of the Earth, offering an effective way to evaluate and compare the strength of the different feedbacks⁵⁵, ⁷⁶, ⁸⁹. This approach has also been successful for some biogeochemical and biogeophysical feedbacks⁶⁶.

In contrast, the evaluation of key non-radiative polar feedbacks is generally inconsistent among the different feedbacks and even among different studies of the same feedback. For the ice growth-thickness feedback, in analogy with the radiative feedback framework, the thickness sensitivity parameter is defined as the ratio of the sea ice thickness change to the perturbative forcing⁴⁹, but this definition has not been widely used so far. The effectiveness of the ice production-entrainment feedback can be measured⁵⁰, ⁵¹ as the ratio between the melting immediately caused by the entrainment of warmer water in the surface layer to the initial ice growth. The ice production-ocean heat storage feedback has been estimated by the ratio between the heat losses associated with sea-ice volume changes to the heat storage below the surface level⁵², ⁵³. Both quantities can be evaluated directly from observations or model results. (Supplementary Note 1).

Those definitions appear justified taken alone but the diversity of definitions and methods to quantify those feedbacks complicates their systematic evaluation and the comparison of the role of the different processes in observed changes. A common framework would thus be very helpful.

**Implications of correctly quantifying feedbacks**

The analysis and quantification of feedbacks have many potential applications. This is illustrated in this section by explaining how this can be used to understand the higher temperature changes expected in the Arctic compared to other regions, to perform a process-oriented evaluation of model behavior, and to reduce the uncertainty in projections.

**Polar amplification.** Overall, climate feedbacks are less stabilizing (i.e., feedback parameters are less negative or more positive) in polar regions than in the tropics. This explains the larger temperature changes experienced in polar regions in response to a perturbation (Fig. 3), a phenomenon referred to as polar amplification¹⁹, ²⁰, ²³, ⁹⁷. For the climate changes projected for the 21st century, polar amplification is much stronger in the Arctic than in the Antarctic. In the Arctic, the large amplification mostly results from (1) a relatively large and positive lapse rate feedback, due to a different vertical distribution of the temperature change compared to the tropics; (2) a relatively weak negative Planck response, due to smaller blackbody emissions per unit warming at lower temperatures (Stefan–Boltzmann law); and (3) a large positive surface albedo feedback, due to the loss of high albedo snow and ice-covered surfaces, as well as a contribution from atmospheric heat transport (Fig. 3a). In the Antarctic, both the weak Planck response and the positive surface albedo feedback induce polar amplification. Warming is damped relative to the Arctic due to a less positive lapse rate feedback, more negative cloud feedback, and strong ocean heat uptake in the Southern Ocean under transient warming (Fig. 3b).

Note that if feedbacks are defined using global mean rather than local surface temperature change, the Planck feedback appears strongly negative in the Arctic because the local temperature change exceeds the local mean. Additionally, the polar amplification has a large seasonal cycle, displaying over the Arctic a minimum in summer and a maximum in fall/winter. In summer, the influence of the large positive surface albedo
<figure><img src="image_7.png" /><figcaption>**Fig. 4** Amplitude of the sea ice volume seasonal cycle versus the ice production-entrainment feedback factor. The ice production-entrainment feedback factor $\gamma_{\theta}$ is defined as the ratio of the melting due to warm water entrainment to the initial ice formation$^{50}$. The values are estimated over 1990-2005 for a standardized perturbation corresponding to an increase of 10 cm of sea ice. For both observational and model datasets, the evaluation of $\gamma_{\theta}$ is performed on the basis of temperature and salinity profiles in the Southern Ocean, averaged over January–February for the period 1990-2005. Values are represented by five crosses corresponding to five sectors of the Southern Ocean$^{111}$. Results for NEMO-LIM$^{112}$, CCSM4$^{113}$, and IPSL$^{114}$ models are in colors. Estimates are given in black based on oceanic observations$^{115}$ and the sea ice volume derived from a reanalysis$^{116}$. For all datasets, the plain circles correspond to the average of $\gamma_{\theta}$ and of the amplitude of the ice volume seasonal cycle over all sectors, and $\rho$ is the correlation coefficient between these two quantities</figcaption></figure>

feedback is compensated by a strong oceanic heat uptake and negative cloud feedbacks while in fall/winter the heat released by the ocean, the contribution of lapse rate feedback and cloud feedbacks induce a large warming$^{55, 98}$.

**Origin of model biases.** In polar regions, many studies have identified model strengths and weaknesses in reproducing observations$^{1, 8, 10, 99, 100}$. Yet, perhaps the most important challenge is identifying the climate system processes that must be represented in order to consider a model realistic enough for mechanistic studies and projections. More generally, simple comparisons between model results and observations do not allow estimating the origin of model biases or their impact. Process-based model evaluation offers the possibility of exploring the causes of discrepancies more deeply, identifying the links between various variables and ultimately suggesting model improvements$^{8, 46, 101-103}$.

To illustrate this point, we compare the ice production-entrainment feedback in three existing simulations for all the sectors of the Southern Ocean with estimates derived from observations and a reanalysis (Fig. 4). A clear link is found between the value of the feedback factor and the amplitude of the seasonal cycle of ice volume: since the ice production-entrainment feedback is negative, it tends to damp the seasonal cycle; i.e., a stronger feedback corresponds to a weaker seasonal cycle. The spread across simulations and Antarctic sea ice regions in Fig. 4 stresses the large sensitivity of the feedback to the ocean properties. As many climate models suffer from large biases in their representation of the vertical structure of the Southern Ocean, they are unlikely to predict this feedback accurately. For instance, the overestimation of the amplitude of the seasonal cycle of sea ice volume in the model CCSM4 is likely related to a too weak negative feedback and improvements in the representation of

ocean properties, in particular of temperature and salinity below the surface layer, would reduce this bias.

**Uncertainties in model projections.** One justification of the development of the radiative feedback framework is to determine the processes that can explain the range of model projections for a specific scenario of future changes in radiative forcing. As expected, the models displaying the largest surface temperature changes are the ones for which the radiative feedbacks have the largest (most positive) values. The same approach can be applied to the non-radiative feedbacks investigated here as they are related to the magnitude of the response to any type of perturbation.

Although it is better if a model is able to reproduce the observations with a bias that is as small as possible, it is not clear for many variables, such as the global mean temperature, that the response to a perturbation is a function of this bias$^{104}$. In other words, there is no a priori reason to believe that a model which reproduces the present global mean temperature well will provide more reasonable projections of future climate than a model that has larger biases. Indeed, a model may have a global mean temperature close to observations due to compensations between many factors that may not necessarily balance in a projected climate$^{105}$.

The situation is distinctive for polar regions, where the feedbacks are strongly non-linear and thus state dependent. This provides an instructive way to interpret the range of model responses as a function of the value of some variables for present-day conditions$^{106, 107}$. Furthermore, it has been argued that a model displaying a more realistic mean state in polar regions will also have a better representation of key processes and thus will provide a more likely estimate of future climate changes than a model with larger biases. This has then been used to justify the selection of models based on their mean state$^{108, 109}$. This idea appears useful in principle but is hard to generalize and is subject to criticism. In particular, it is not always clear to determine how to evaluate models, which variable should be used to select models, and if the currently available model sample is adequate to apply a meaningful selection. This has led to strong debate in the community about the justification of this approach, which may artificially reduce the uncertainty range by discarding model results that are as likely as the others$^{110}$. We propose to use feedback quantification more extensively to evaluate model behavior, and foresee that this can contribute to more robust estimates of the likelihood of projections.

**A simple and consistent approach for non-radiative feedbacks**

The review above illustrates that many definitions and evaluation methods have been proposed for the various radiative and non-radiative feedbacks. Nevertheless, all the feedbacks can be described and quantified using a simple and consistent framework, based on the definition of a feedback factor $\gamma$.

For radiative feedbacks, the feedback factor $\gamma_i$ is the ratio of a particular feedback parameter to minus the Planck feedback parameter. An analogous expression can be written for any other feedback. When only one feedback is operating (see the methods for the case of multiple feedbacks), the feedback factor $\gamma$ can be quantified as the ratio between the additional changes specifically due to the feedback and the response of the full system including all the feedbacks (Total response). This additional change (Total response – Reference response) is itself computed as the difference between the response of the full system and the one of a reference system in which the feedback under consideration does
not operate (Reference response):

$$ \gamma = \frac{\text{Total response} - \text{Reference response}}{\text{Total response}} \quad (1) $$

The methodology requires explicitly identifying (1) a perturbation or a class of perturbations, (2) a response variable involved in the feedback loop, (3) the full system with all processes operating and its response to the perturbation, and (4) the reference system with the process of interest not operating and the reference system response to the perturbation. While the framework is general, a clear definition of this system is required, as the value of the feedback factor depends on the way the reference system is chosen¹². Let's examine radiative feedbacks as one example (see the methods): (1) the perturbation is the radiative forcing F, (2) the response variable is T<sub>s</sub>, (3) the full system includes one or more radiative feedbacks plus the Planck response referenced to T<sub>s</sub>, and (4) the reference system is the Planck response only.

As it is based on the same principles, analysis of non-radiative processes using this feedback factor retains the main advantages of the radiative feedback framework. First, each feedback can be associated with a well-defined conceptual model describing the mechanisms and interactions involved. This is essential in order to allow each feedback to be firmly rooted in a process-based analysis that is straightforward to apply and understand. Secondly, it is possible to evaluate the magnitude of the feedbacks using a dimensionless factor, ideally both in models and observations. This is required to assess the contributions of the different feedbacks to the total response and to compare the role of each feedback in various Earth System models to determine which is responsible for their distinctive sensitivities.

As in the example below, the feedback factor can in some cases be evaluated using observations or model outputs only, but it may also require specific additional calculations. This is illustrated using a simple model in Supplementary Note 2 for the ice growth-ice thickness feedback. In this case, potential compensations can occur between feedbacks and the interpretation of the estimates that are obtained must then take into those synergies⁶⁵, ⁸⁹.

**An example of the approach for the ice production-ocean entrainment feedback.** We illustrate the methodology with the ice production–entrainment feedback⁵⁰. For this negative feedback, (1) the perturbation is a given amount of ice production, (2) the reference variable is ice thickness, (3) the full system is the sea ice plus ocean column with the entrainment process, and (4) the reference system is the sea ice plus ocean column but without entrainment. The intensity of this feedback can then be evaluated using the ratio

$$ \gamma_{\theta} = \frac{\text{Total ice thickness changes-Ice thickness changes without entrainment}}{\text{Total ice thickness changes}} \quad (2) $$

Despite a different form, this expression is strictly equivalent to the original formulation proposed in ref. ⁵⁰ and used for Fig. 4 (see Supplementary Note 1 for the demonstration).

As the mixed layer deepens, it entrains water with increasing temperature (since temperature increases with depth) and the heat input grows. Consequently, the absolute value of the feedback factor γ<sub>θ</sub> increases with ice formation (meaning that its value decreases, since it is a negative feedback) until the end of winter (Fig. 5). This non-linear behavior can be illustrated using a simple analytical model as shown in Supplementary Note 3 and Supplementary Fig. 3.

<figure><img src="image_8.png" /><figcaption>**Fig. 5** Evolution of the ice production-entrainment feedback factor as a function of ice production. For each value of ice production, the entrainment is computed from the January–February 1990–2005 mean temperature and salinity profiles¹¹⁵ assuming a mixed layer deepening restoring the static stability of the water column after the brine release. It is shown here for a Weddell Sea location typically covered by ice in winter (near 30°W, 65°S). The strength of the pycnocline is thus evaluated in summer but it must be measured below the layer close to the surface that is warmed above freezing point temperature if sea ice completely melts, as the heat in this layer is removed quickly in fall when the temperature drops and is not involved in the ice production-ocean entrainment feedback</figcaption></figure>

In practice, it is usually not possible to completely quantify all the dependencies of a non-linear, spatially variable feedback factor. This is why it has been suggested to select a prescribed perturbation representative of the condition of interest, similarly to the classical analysis of radiative forcing in terms of doublings of atmospheric CO₂ concentrations (Fig. 2). For the ice production-ocean entrainment feedback, we propose to evaluate γ<sub>θ</sub> by considering the response to a standardized perturbation corresponding to an increase in sea ice thickness of 10 cm (Fig. 4). The number of observed profiles being much higher in summer, the feedback parameter is evaluated from data collected during this season. It is clear that the corresponding values of the feedback factors are not universally valid but they provide a standard benchmark for comparisons and analyses.

## Concluding remarks

This Perspective underlines the critical role of feedbacks in the dynamics of polar climate and the need to quantify them precisely. Feedback quantification provides a powerful tool to understand the interactions between the components of the system, to analyze model biases and to determine the origin of the differences within a set of model predictions. We have focused on some physical processes affecting the atmosphere, sea ice, ice sheets, land surfaces and ocean in polar regions. Yet, the discussion can be extended to feedbacks including biogeochemical processes.

Quantification of feedback strength is not simple as many polar feedbacks are strongly non-linear. Indeed, feedback magnitude depends on the location, the season and is a function of the climate state. We have provided here values in specific conditions for some of the feedbacks investigated. However, one single number or a range, as could be included in a table for instance, is not sufficient to fully characterize a feedback and its variations.

While the quantitative estimation of feedbacks follows well-established methodologies for radiative feedbacks, this is not the case for many other feedbacks. Nevertheless, the traditional radiative feedback analysis can be extended to define a feedback factor that can be used as a standard measure of most polar
feedbacks. This feedback factor is estimated as the relative contribution of the feedback to the total change of the system in response to a perturbation. It has several advantages: (1) it is compatible with the radiative feedback framework which has proven useful over the past three decades, (2) it is based on a clear physical interpretation of key physical processes, (3) the framework is simple to articulate, and (4) an unambiguous quantification can be obtained allowing an objective evaluation of the processes that control the response to a perturbation.

Using a feedback factor provides the clear theoretical advantages of being consistent with the feedback theory, general and easy to interpret. Nevertheless, based on data availability or on the specific goal of a study, other parameters will continue to be used to diagnose the magnitude of some feedbacks. This may of course be perfectly justified. However, the limitations and merit of such an approach compared to a more general view need to be discussed and taken into account when interpreting the results.

Even though the feedbacks discussed here are well-known, some are often only used to provide a qualitative interpretation or a narrative framework to explain the changes occurring in polar regions. In particular, the quantitative evaluation of non-radiative feedbacks in models is rare, with only few and, in some cases, no publications on this topic for some feedbacks. This Perspective motivates the use of more systematic approaches to analyze past and upcoming model results. In particular, the framework presented here allows tracing the origin of model diversity back to physical considerations. Identifying the feedbacks that are critical for correctly simulating the mean state and variability of polar climate will ultimately promote the development of targeted observational campaigns, by means of which models will be evaluated. Such campaigns already exist or are underway: the Year of Polar Prediction (YOPP) or the Multidisciplinary Drifting Observatory for the Study of Arctic Climate (MOSAiC) are two examples. In summary, advancing evaluation of feedbacks will require dedicated model experiments and careful analyses, complemented with the collection of dedicated observations that can constrain model feedbacks. Some of those elements are clearly challenging, but this will be strongly beneficial for our understanding of polar climate dynamics and of the future changes expected at high latitudes.

## Methods

**Radiative feedbacks expressed as feedback factors.** The global mean radiative balance at the TOA in response to a radiative forcing $F$ (in W m⁻²) at equilibrium can be written as

$$F + \left( \lambda_0 + \sum_i \lambda_i \right) \Delta T_s = 0 \quad (3)$$

where $\Delta T_s$ (in K) is the surface temperature change, $\lambda_0$ (~-3.2 W m⁻² K⁻¹) the Planck response and the $\lambda_i$ (in W m⁻² K⁻¹) correspond to the radiative feedback parameters related to the response of surface albedo, clouds, water vapor and vertical temperature gradient (lapse rate feedback).

In the absence of these feedbacks, the equilibrium surface temperature change in response to a doubling of CO₂ would be governed by only the Planck response, given by $\Delta T_0 = -F/\lambda$ and equal to approximately¹² 1.2 K. This temperature change is amplified or damped by individual feedbacks, depending on whether they contribute to a positive or negative additional radiative perturbation to the TOA radiative balance in response to warming. The overall equilibrium warming resulting from a CO₂ doubling is thus greater than $\Delta T_0$, likely¹²⁻¹⁴ between about 2 and 4.5 K.

The primary advantage of the feedback framework is that it allows a process-based analysis of the adjustment of the system to a radiative perturbation and a quantification of the importance of each process. Traditionally, the magnitude of each feedback is compared to that of the Planck response, giving dimensionless feedback factors $\gamma_i = \lambda_i/\lambda_0$. In turn, the surface temperature change in response to

forcing can be cast in terms of the feedback factors as

$$\Delta T_s = \Delta T_0 / \left( 1 - \sum_i \gamma_i \right) \quad (4)$$

The sum of all the feedback factors $\gamma_g = \sum_i \gamma_i = \sum_i \frac{\lambda_i}{\lambda_0}$ can also be calculated as

$$\gamma_g = \frac{\Delta T_s - \Delta T_0}{\Delta T_s} \quad (5)$$

This sum can thus be interpreted as the additional warming due to the feedbacks divided by the total temperature change.

**Radiative feedbacks expressed as warming contributions.** The formalism can be extended to assess the relative contributions of individual feedbacks to local surface warming by use of the local energy budget equation:

$$F + \left( \lambda_0 + \sum_i \lambda_i \right) \Delta T_s + \Delta \text{OHU} + \Delta \text{AHT} = 0 \quad (6)$$

where each variable is a function of latitude, $\Delta \text{OHU}$ is the change in ocean heat uptake (positive into atmospheric column) and $\Delta \text{AHT}$ represents the change in atmospheric heat flux convergence (positive into atmospheric column). Following previous studies¹⁸, ¹⁹, ⁷², we define the warming contribution as the energetic contribution (in W m⁻²) associated with a particular feedback ($\lambda_i \Delta T_s$) or atmospheric forcing ($F$, $\Delta \text{OHU}$ or $\Delta \text{AHT}$) divided by the magnitude of the global-mean Planck response $-\bar{\lambda}_0$:

$$\Delta T_s = -F/\bar{\lambda}_0 - \lambda_0 \Delta T_s/\bar{\lambda}_0 - \sum_i \lambda_i \Delta T_s/\bar{\lambda}_0 - \Delta \text{OHU}/\bar{\lambda}_0 - \Delta \text{AHT}/\bar{\lambda}_0 \quad (7)$$

where the terms on the right-hand side each represent an individual warming contribution and together sum to the total surface warming $\Delta T_s$ (with small residual ignored here); $\lambda_0 = \lambda_0 - \bar{\lambda}_0$ represents the deviation of the local Planck response from its global-mean value. Here we use local feedbacks $\lambda$ diagnosed using radiative kernels⁸⁵ averaged over years 85–115 of abrupt CO₂ quadrupling simulations of 13 models⁷²; $\Delta \text{OHU}$ is diagnosed as the anomalous net surface heat flux and $\Delta \text{AHT}$ as the difference between $\Delta \text{OHU}$ and net TOA radiation flux anomalies. Figure 3 shows the calculated warming contributions of feedbacks and forcings in polar and tropical regions.

**Feedback factor and the feedback gain when several feedbacks are operating.**

When several feedbacks are operating, their contribution to the changes in the radiative balance is additive. Using the same notations as above, the radiative balance at equilibrium is

$$F + \left( \lambda_0 + \sum_i \lambda_i \right) \Delta T_s = 0 \quad (8)$$

leading when dividing by $\lambda_0$ to

$$\frac{F}{\lambda_0} + \left( 1 + \sum_i \frac{\lambda_i}{\lambda_0} \right) \Delta T_s = 0 \quad (9)$$

and

$$-\Delta T_0 + \left( 1 - \sum_i \gamma_i \right) \Delta T_s = 0 \quad (10)$$

The contribution of the feedback factors is thus also additive.

This is not the case for feedback gain $G$

$$G = \frac{\Delta T_s}{\Delta T_0} = \frac{1}{1 - \sum_i \gamma_i} = \frac{1}{1 - \gamma_g} \quad (11)$$

As the various feedbacks are interacting to get the full response, the gain when two feedbacks are active is not the sum of the gains when the feedbacks are acting separately.

When only one feedback is acting or if only the sum of all feedbacks is considered, $\gamma_g$ can be simply evaluated by

$$\gamma_g = \frac{\Delta T_s - \Delta T_0}{\Delta T_s} \quad (12)$$

Several techniques are available⁸², ⁸⁵ to estimate the individual $\gamma_i$. They generally require to perform specific analyses to extract the contribution of a particular feedback. One example is to perform an experiment when only the investigated feedback is operating and comparing the changes $\Delta T_{si}$ in this
experiment to the one of the reference system $\Delta T_0$.

$$y_l = \frac{\Delta T_{si} - \Delta T_0}{\Delta T_{si}}. \qquad (13)$$

This solution is not the one traditionally used for radiative feedbacks but the approach can be generalized to any feedback where this alternative may be a practical solution for feedback evaluation.

Received: 19 September 2017 Accepted: 4 April 2018

Published online: 15 May 2018

## References

1. Stroeve, J. C. et al. The Arctic's rapidly shrinking sea ice cover: a research synthesis. *Clim. Change* **110**, 1005 (2012).

2. Stammerjohn, S., Massom, R., Rind, D. & Martinson, D. Regions of rapid sea ice change: an inter-hemispheric seasonal comparison. *Geophys. Res. Lett.* **39**, L06501 (2012).

3. Döscher, R., Vihma, T. & Maksimovich, E. Recent advances in understanding the Arctic climate system state and change from a sea ice perspective: a review. *Atmos. Chem. Phys.* **14**, 13751–13600 (2014).

4. Armour, K. C., Marshall, J., Scott, J., Donohoe, A. & Newsom, E. R. Southern Ocean warming delayed by circumpolar upwelling and equatorward transport. *Nat. Geos.* **9**, 549–554 (2016).

5. Jones, J. M. et al. Assessing recent trends in high-latitude Southern Hemisphere surface climate. *Nat. Clim. Change* **6**, 917–926 (2016).

6. Hobbs, W. R. et al. A review of recent changes in Southern Ocean sea ice, their drivers and forcings. *Glob. Planet. Change* **143**, 228–250 (2016).

7. Swart, N. C., Fyfe, J. C., Hawkins, E., Kay, J. E. & Jahn, A. Influence of internal variability on Arctic sea-ice trends. *Nat. Clim. Change* **5**, 86–89 (2015).

8. Notz, D. How well must climate models agree with observations? *Phil. Trans. R. Soc. A* **373**, 20140164. https://doi.org/10.1098/rsta.2014.0164 (2015).

9. Zunz, V., Goosse, H. & Massonnet, F. How does internal variability influence the ability of CMIP5 models to reproduce the recent trend in Southern Ocean sea ice extent? *Cryosphere* **7**, 451–468 (2013).

10. Flato, G. et al. Evaluation of Climate Models. *Climate Change: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change* (eds Stocker, T.F., Qin, D., Plattner, G.-K., Tignor, M., Allen, S. K., Boschung, J., Nauels, A., Xia, Y., Bex, V., & Midgley, P.M.) (Cambridge United Kingdom and New York, USA, 2013).

11. Aström, K. J. & Murray, R. M. Feedback systems: an introduction for scientists and engineers. (Princeton, New Jersey, USA, 2010) (ISBN 9781400828739).

12. Roe, G. H. Feedbacks, timescales and seeing red. *Annu. Rev. Earth. Planet. Sci.* **37**, 93–115 (2009).

13. Hansen, J. E. et al. Climate sensitivity: analysis of feedback mechanisms. *Climate processes and climate sensitivity* (eds Hansen, J. E. & Takahashi, T.) 130–163 (American Geophysical Union, Washington, DC, USA, 1984).

14. Bony, S. et al. How well do we understand and evaluate climate change feedback processes? *J. Clim.* **19**, 344–3482 (2006).

15. Wallace, J. M. & Hobbs, P. V. Atmospheric science: an introductory survey 2nd edn. *International Geophysics Series* **92**, 484 (Academic press, Burlington, MA, USA and San Diego, CA, USA, 2006).

16. Armour, K. C., Bitz, C. M. & Roe, G. H. Time-varying climate sensitivity from regional feedbacks. *J. Clim.* **26**, 4518–4534 (2013).

17. Ramanathan, V. Interactions between ice-albedo, lapse-rate and cloud-top feedbacks: an analysis of the nonlinear response of a GCM climate model. *J. Atmos. Sci.* **34**, 1885–1897 (1977).

18. Crook, J. A., Forster, P. M. & Stuber, N. Spatial patterns of modeled climate feedback and contributions to temperature response and polar amplification. *J. Clim.* **24**, 3575–3592 (2011).

19. Pithan, F. & Mauritsen, T. Arctic amplification dominated by temperature feedbacks in contemporary climate models. *Nat. Geosci.* **7**, 181–184 (2014).

20. Manabe, S. & Wetherald, R. The effects of doubling the CO₂ concentration on the climate of a general circulation model. *J. Atmos. Sci.* **32**, 3–15 (1975).

21. Dessler, A. E., Zhang, Z. & Yang, P. Water-vapor climate feedback inferred from climate fluctuations, 2003–2008. *Geophys. Res. Lett.* **35**, L20704 (2008).

22. Gordon, N. D., Jonko, A. K., Forster, P. M. & Shell, K. M. An observationally based constraint on the water-vapor feedback. *J. Geophys. Res. Atmos.* **118**, 12,435–12,443 (2013).

23. Taylor, P. C. et al. A decomposition of feedback contributions to polar warming amplification. *J. Clim.* **26**, 7023–7043 (2013).

24. Graversen, R. G. & Wang, M. Polar amplification in a coupled climate model with locked albedo. *Clim. Dyn.* **33**, 629–643 (2009).

25. Budyko, M. The effect of solar radiation variations on the climate of the Earth. *Tellus* **21**, 611–619 (1969).

26. Sellers, P. A global climate model based on the energy balance of the earth-atmosphere system. *J. Appl. Meteor.* **8**, 392–400 (1969).

27. Hall, A. The role of surface albedo feedback in climate. *J. Clim.* **17**, 1550–1568 (2004).

28. Winton, M. Surface albedo feedback estimates for the AR4 climate models. *J. Clim.* **19**, 359–365 (2006).

29. Perovich, D. K. et al. Thin and thinner: sea ice mass balance measurements during SHEBA. *J. Geophys. Res.* **108**, 8050 (2003).

30. Sturm, M. & Massom, R. A. Snow in the sea ice system: friend or foe? In *Sea Ice* (ed. Thomas, D. N.) 65–109 (Wiley-Blackwell, Oxford, United Kingdom, 2017).

31. Wetherald, R. T. & Manabe, S. Cloud feedback processes in a general circulation model. *J. Atmos. Sci.* **45**, 1397–1415 (1988).

32. Mitchell, J. F. B., Senior, C. A. & Ingram, W. J. On CO₂ and climate: a missing cloud feedback? *Nature* **341**, 132–134 (1989).

33. Vial, J., Dufresne, J.-L. & Bony, S. On the interpretation of inter-model spread in CMIP5 climate sensitivity estimates. *Clim. Dyn.* **41**, 3339–3362 (2013).

34. Zelinka, M. D., Klein, S. A. & Hartmann, D. L. Computing and partitioning cloud feedbacks using cloud property histograms. Part II: attribution to changes in cloud amount, altitude, and optical depth. *J. Clim.* **25**, 3736–3754 (2012).

35. Andrews, T., Gregory, J. M., Webb, M. J. & Taylor, K. E. Forcing, feedbacks and climate sensitivity in CMIP5 coupled atmosphere-ocean climate models. *Geophys. Res. Lett.* **39**, L09712 (2012).

36. Schweiger, A. L., Lindsay, R. W., Vavrus, S. & Francis, J. A. Relationships between Arctic sea ice and clouds during Autumn. *J. Clim.* **21**, 4799–4810 (2008).

37. Morrison, A. L., Kay, J. E., Chepfer, H., Guzman, R. & Yettella, V. Isolating the liquid cloud response to recent Arctic sea ice loss using spaceborne lidar observations. *J. Geophys. Res. Atmos.* **123**, 473–490 (2018).

38. Kay, J. E. et al. Recent advances in Arctic cloud and climate research. *Curr. Clim. Change Rep.* **2**, 159 (2016).

39. Boisvert, L. N., Wu, D. L. & Shie, C. -L. Increasing evaporation amounts seen in the Arctic between 2003 and 2013 from AIRS data. *J. Geophys. Res. Atmos.* **120**, 6865–6881 (2015).

40. Bodas-Salcedo, A., Andrews, T., Karmalkar, A. V. & Ringer, M. A. Cloud liquid water path and radiative feedbacks over the Southern Ocean. *Geophys. Res. Lett.* **43**, 10,938–10,946 (2016).

41. Wall, C. J., Kohyama, T. & Hartmann, D. L. Low-cloud, boundary layer, and sea ice interactions over the Southern Ocean during winter. *J. Clim.* **30**, 4857–4871 (2017).

42. Terai, C. R., Zelinka, M. & Klein, S. A. Constraining the low-cloud optical depth feedback at middle and high latitudes using satellite observations. *J. Geophys. Res. Atmos.* **121**, 9696–9716 (2016).

43. Ceppi, P., McCoy, D. T. & Hartmann, D. L. Observational evidence for a negative shortwave cloud feedback in mid to high latitudes. *Geophys. Res. Lett.* **43**, 1331–1339 (2016).

44. Bodas-Salcedo, A. et al. Origins of the solar radiation biases over the Southern Ocean in CMIP2 models. *J. Clim.* **27**, 41–56 (2014).

45. Tan, I., Storelvmo, T. & Zelinka, M. D. Observational constraints on mixed-phase clouds imply higher climate sensitivity. *Science* **352**, 224–227 (2016).

46. Kay, J. E. et al. Evaluating and improving cloud phase in the Community Atmosphere Model version 5 using spaceborne lidar observations. *J. Geophys. Res. Atmos.* **121**, 4162–4176 (2016).

47. Holland, M. M., Bitz, C. M. & Tremblay, B. Future abrupt reductions in the summer Arctic sea ice. *Geophys. Res. Lett.* **33**, L23503 (2006).

48. Maykut, G. A. The surface heat and mass balance. *The Geophysics of Sea Ice* (ed. Untersteiner, N.) 395–464 (Plenum Press, 1986).

49. Bitz, C. M. & Roe, G. H. A mechanism for the high rate of sea ice thinning in the Arctic Ocean. *J. Clim.* **17**, 3623–3632 (2004).

50. Martinson, D. G. Evolution of the Southern Ocean winter mixed layer and sea ice-open ocean deep-water formation and ventilation. *J. Geophys. Res. Oceans* **95**, 11641–11654 (1990).

51. Martinson, D. G. & Iannuzzi, R. A. Antarctic ocean-ice interaction: implications from ocean bulk property distributions in the Weddell Gyre. *Antarctic sea ice: physical processes, interactions and variability*, (ed. Jeffries, M.) 243–271 (American Geophysical Union, Washington, DC, USA, 1998).

52. Goosse, H. & Zunz, V. Decadal trends in the Antarctic sea ice extent ultimately controlled by ice-ocean feedback. *Cryosphere* **8**, 453–470 (2014).

53. Lecomte, O. et al. Vertical ocean heat redistribution sustaining sea-ice concentration trends in the Ross Sea. *Nat. Comm.* **8**, 258 (2017).

54. Seneviratne, S. I. et al. Investigating soil moisture-climate interactions in a changing climate: a review. *Earth-Sci. Rev.* **99**, 125–161 (2010).

55. Lainé, A., Yoshimori, M., & Abe-Ouchi, A. Surface Arctic amplification factors in CMIP5 models: land and oceanic surfaces and seasonality. *J. Clim.* **29**, 3297–3316 (2016).
56. Edwards, T. L. et al. Effect of uncertainty in surface mass balance-elevation feedback on projections of the future sea level contribution of the Greenland ice sheet. *Cryosphere* **8**, 195-208 (2014).

57. Edwards, T. L. et al. Probabilistic parameterization of the surface mass balance-elevation feedback in regional climate model simulations of the Greenland ice sheet. *Cryosphere* **8**, 181-194 (2014).

58. Schoof, C. Ice sheet grounding line dynamics: steady states, stability, and hysteresis. *J. Geophys. Res.* **112**, F03S28 (2007).

59. Docquier, D., Perichon, L. & Pattyn, F. Representing grounding line dynamics in numerical ice sheet models: recent advances and outlook. *Surv. Geophys.* **32**, 417-435 (2011).

60. Rignot, E., Mouginot, J., Morlighem, M., Seroussi, H. & Scheuchl, B. Widespread, rapid grounding line retreat of Pine Island, Thwaites, Smith, and Kohler glaciers, West Antarctica, from 1992 to 2011. *Geophys. Res. Lett.* **41**, 3502-3509 (2014).

61. Wise, M. G., Dowdeswell, J. A., Jakobsson, M. & Larter, R. D. Evidence of marine ice-cliff instability in Pine Island Bay from iceberg-keel plough marks. *Nature* **550**, 506-510 (2017).

62. Zwally, H. J. et al. Surface melt-induced acceleration of Greenland ice-sheet flow. *Science* **297**, 218-222 (2002).

63. Swingedouw, D. et al. Antarctic ice-sheet melting provides negative feedbacks on future climate warming. *Geophys. Res. Lett.* **35**, L17705 (2008).

64. Bintanja, R., van Oldenborgh, G. J., Drijfhout, S. S., Wouters, B. & Katsman, C. A. Important role for ocean warming and increased ice-shelf melt in Antarctic sea-ice expansion. *Nat. Geos.* **6**, 376-379 (2013).

65. Otto, J., Raddatz, T., Claussen, M., Brovkin, V. & Gayler, V. Separation of atmosphere-ocean-vegetation feedbacks and synergies for mid-Holocene climate. *Geophys. Res. Lett.* **36**, L09701 (2009).

66. Gregory, J. M., Jones, C. D., Cadule, P. & Friedlingstein, P. Quantifying carbon cycle feedbacks. *J. Clim.* **22**, 5232-5250 (2009).

67. Arneth, A. Terrestrial biogeochemical feedbacks in the climate system. *Nat. Geosci.* **3**, 525-532 (2010). Coauthors.

68. Schuur, E. A. G. Climate change and the permafrost carbon feedback. *Nature* **520**, 171-179 (2015). Coauthors.

69. Pabi, S., van Dijken, G. L. & Arrigo, K. R. Primary production in the Arctic Ocean, 1998-2006. *J. Geophys. Res.* **113**, C8 (2008).

70. Lengaigne, M. et al. Bio-physical feedbacks in the Arctic Ocean using an Earth system model. *Geophys. Res. Lett.* **36**, L21602 (2009).

71. Alexeev, V. A. & Jackson, C. H. Polar amplification: is atmospheric heat transport important? *Clim. Dyn.* **41**, 533-547 (2013).

72. Feldl, N., Bordoni, S. & Merlis, T. M. Coupled high-latitude climate feedbacks and their impact on atmospheric heat transport. *J. Clim.* **30**, 189-201 (2017).

73. Kay, J. E. et al. The influence of local feedbacks and northward heat transport on the equilibrium Arctic climate response to increased greenhouse gas forcing in coupled climate models. *J. Clim.* **25**, 5433-5450 (2012).

74. Roe, G. H., Feldl, N., Armour, K. C., Hwang, Y.-T. & Frierson, D. M. W. The remote impacts of climate feedbacks on regional climate predictability. *Nat. Geosci.* **8**, 135-139 (2015).

75. Cai, M. & Lu, J. Dynamical greenhouse-plus feedback and polar warming amplification. Part II: meridional and vertical asymmetries of the global warming. *Clim. Dyn.* **29**, 375-391 (2007).

76. Feldl, N. & Roe, G. H. The nonlinear and nonlocal nature of climate feedbacks. *J. Clim.* **26**, 8289-8304 (2013).

77. Zelinka, M. D. & Hartmann, D. L. Climate feedbacks and their implications for poleward energy flux changes in a warming climate. *J. Clim.* **25**, 608-624 (2011).

78. Overland, J. E. et al. Nonlinear response of mid-latitude weather to the changing Arctic. *Nat. Clim. Change* **6**, 992-999 (2016).

79. Marshall, J. et al. The ocean's role in the transient response of climate to abrupt greenhouse gas forcing. *Clim. Dyn.* **4**, 2287-2299 (2015).

80. Jungclaus, J. H., Lohmann, K. & Zanchettin, D. Enhanced 20th-century heat transfer to the Arctic simulated in the context of climate variations over the last millennium. *Clim. Past.* **10**, 2201-2213 (2014).

81. Bitz, C. M., Gent, P. R., Woodgate, R. A., Holland, M. M. & Lindsay, R. The influence of sea ice on ocean heat uptake in response to increasing CO₂. *J. Clim.* **19**, 2437-2450 (2006).

82. Gregory, J. M. et al. A new method for diagnosing radiative forcing and climate sensitivity. *Geophys. Res. Lett.* **31**, L03205 (2004).

83. Taylor, K. E. et al. Estimating shortwave radiative forcing and response in climate models. *J. Clim.* **20**, 2530-2543 (2007).

84. Soden, B. F. et al. Quantifying climate feedbacks using radiative kernels. *J. Clim.* **21**, 3504-3520 (2008).

85. Shell, K. M., Kiehl, J. T. & Shields, C. A. Using the radiative kernel technique to calculate climate feedbacks in NCAR's community atmospheric model. *J. Clim.* **21**, 2269-2282 (2008).

86. Jonko, A., Shell, K., Sanderson, B. & Danabasoglu, G. Climate feedbacks in CCSM3 under changing CO₂ forcing. Part I: adapting the linear radiative

kernel technique to feedback calculations for a broad range of forcings. *J. Clim.* **25**, 5260-5272 (2012).

87. Jonko, A. K., Shell, K. M., Sanderson, B. M. & Danabasoglu, G. Climate feedbacks in CCSM3 under changing CO₂ forcing: Part II. Variation of climate feedbacks and sensitivity with forcing. *J. Clim.* **26**, 2784-2795 (2013).

88. Sedlar, J., Shupe, M. D. & Tjernström, M. On the relationship between thermodynamic structure and cloud top, and its climate significance in the Arctic. *J. Clim.* **25**, 2374-2393 (2012).

89. Cai, M. & Lu, J. A new framework for isolating individual feedback processes in coupled general circulation climate models. Part II: method demonstrations and comparisons. *Clim. Dyn.* **32**, 887-900 (2009).

90. Feldl, N. & Roe, G. H. Four perspectives on climate feedbacks. *Geophys. Res. Lett.* **40**, 4007-4011 (2013).

91. Sherwood, S. C. et al. Adjustments in the forcing-feedback framework for understanding climate change. *Bull. Am. Meteor. Soc.* **96**, 217-228 (2015).

92. Andrews, T., Gregory, J. M. & Webb, M. J. The dependence of radiative forcing and feedback on evolving patterns of surface temperature change in climate models. *J. Clim.* **28**, 1630-1648 (2015).

93. Colman, R. A. & McAvaney, B. J. On tropospheric adjustment to forcing and climate feedbacks. *Clim. Dyn.* **36**, 1649-1658 (2011).

94. Meraner, K., Mauritsen, T. & Voigt, T. A. Robust increase in equilibrium climate sensitivity under global warming. *Geophys. Res. Lett.* **40**, 5944-5948 (2013).

95. Gregory, J. M., Andrews, T. & Good, P. The inconstancy of the transient climate response parameter under increasing CO₂. *Philos. Trans. Roy. Soc. Lond.* **373**, 20140417 (2015).

96. Rugenstein, M. A. A., Caldeira, K. & Knutti, R. Dependence of global radiative feedbacks on evolving patterns of surface heat fluxes. *Geophys. Res. Lett.* **43**, 9877-9885 (2016).

97. Holland, M. M. & Bitz, C. M. Polar amplification of climate change in coupled models. *Clim. Dyn.* **21**, 221-232 (2003).

98. Sejas, S. A. et al. Individual feedback contributions to the seasonality of surface warming. *J. Clim.* **27**, 5653-5669 (2014). 410.

99. Svensson, G. & Karlsson, J. On the Arctic wintertime climate in global climate models. *J. Clim.* **24**, 5757-5771 (2011).

100. de Boer, G., Chapman, W., Kay, J. E., Medeiros, B. & Shupe, M. D. A characterization of the present-day Arctic atmosphere in CCSM4. *J. Clim.* **25**, 2676-2695 (2012).

101. Qu, X. & Hall, A. What controls the strength of snow-albedo feedback? *J. Clim.* **20**, 3971-3981 (2007).

102. Notz, D. et al. The CMIP6 Sea-Ice Model Intercomparison Project (SIMIP): understanding sea ice through climate-model simulations. *Geosci. Model Dev.* **9**, 3427-3446 (2016).

103. Pithan, F. et al. Select strengths and biases of models in representing the Arctic winter boundary layer over sea ice: the Larcform 1 single column model intercomparison. *J. Adv. Mod. Earth Syst.* **8**, 1345-1357 (2016).

104. Hawkins, E. & Sutton, R. Connecting climate model projections of global temperature change with the real World. *Bull. Am. Meteor. Soc.* **6**, 963-980 (2016).

105. Knutti, R. Why are climate models reproducing the observed global surface warming so well? *Geophys. Res. Lett.* **35**, L18704 (2008).

106. Bitz, C. M. Some aspects of uncertainty in predicting sea ice thinning. Arctic sea ice decline: observations, projections, mechanisms, and implications. *Amer. Geophys. Union* (eds deWeaver, E., Bitz, C. M. & Tremblay, B.) 63-76 (2008).

107. Meehl, G. A. & Washington, W. M. Climate change projections in CESM1 (CAM5) compared to CCSM4. *J. Clim.* **26**, 6287-6308 (2013).

108. Wang, M. & Overland, J. E. A sea ice free summer Arctic within 30 years: an update from CMIP5 models. *Geophys. Res. Lett.* **39**, L18501 (2012).

109. Massonnet, F. et al. Constraining projections of summer Arctic sea ice. *Cryosphere* **6**, 1383-1394 (2012).

110. Eyring, V. et al. Towards improved and more routine Earth system model evaluation in CMIP. *Earth Syst. Dynam.* **7**, 813-830 (2016).

111. Gloersen, P. et al. Arctic and Antarctic sea ice, 1978-1987: satellite passive microwave observations and analysis. National Aeronautics and Space Administration, NASA SP-511 (1992).

112. Barthélemy, A., Fichefet, T., Goosse, H. & Madec, G. Modeling the interplay between sea ice formation and the oceanic mixed layer: limitations of simple brine rejection parameterizations. *Ocean Model.* **86**, 141-152 (2015).

113. Gent, P. R. et al. The community climate system model version 4. *J. Clim.* **24**, 4 973-4 991 (2011).

114. Dufresne, J. L. et al. Climate change projections using the IPSL-CM5 earth system model: from CMIP3 to CMIP. *Clim. Dyn.* **40**, 2123-2165 (2013).

115. Ridgway, K. R., Dunn, J. R. & Wilkin, J. L. Ocean interpolation by four-dimensional weighted least squares-application to the waters around Australasia. *J. Atmos. Ocean. Tech.* **19**, 1357-1375 (2002).

116. Massonnet, F. et al. A model reconstruction of the Antarctic sea ice thickness and volume changes over 1980-2008 using data assimilation. *Ocean Model.* **64**, 67-75 (2013).
117. Donohoe, A., Armour, K. C., Pendergrass, A. G. & Battisti, D. S. Shortwave and longwave radiative contributions to global warming under increasing CO₂. *Proc. Nat. Acad. Sci. 111*, 16700–16705 (2014).

## Acknowledgements

This work was supported by FNRS (Grant agreement no. T.0007.14), the PRIMAVERA EU-H2020 project (Grant agreement no. 641727), the Ministerio de Economía, Industria y Competitividad (MINECO) and Belgian Science Policy Office. K.C. Armour was supported by National Science Foundation award OCE-1523641. H. Goosse is the Research Director within the Fonds National de la Recherche Scientifique (F.R.S.-FNRS-Belgium). François Massonnet is a F.R.S-FNRS Post-Doctoral Researcher. Jennifer Kay's efforts were supported by NASA 15-CCST15-0025. A. Bodas-Salcedo was supported by the Joint UK BEIS/Defra Met Office Hadley Centre Climate Programme (GA01101). P. Kushner acknowledges the support of the Natural Sciences and Engineering Research Council of Canada's Canadian Sea Ice and Snow Evolution Network (see www.CanSISE.ca). This paper is the outcome of a workshop held in May 2016 at the Université catholique de Louvain, Belgium in the framework of the Polar Climate Prediction Initiative supported by the World Climate Research Programme.

## Author contributions

All authors shared responsibility for writing the manuscript and of the general concept that is developed. H.G. and J.E.K. conceived and supervised the study. H.S.P. and A.J. lead the preparation of Fig. 1. A.J. prepared Fig. 2 and wrote the text devoted to it. K.C.A. supplied Fig. 3 and contributed to the discussion of the links with radiative feedbacks in collaboration with all authors. O.L. and H.G. were in charge of the discussion of ice production-entrainment feedback (Figs. 4 and 5). F.M and M.V. wrote the section devoted to the sea ice growth-thickness feedback (including Supplementary Fig. 2). Supplementary Fig. 1 was performed and discussed by A.B.-S. D.D. was responsible of Table 1, with inputs of all authors.

## Additional information

**Supplementary Information** accompanies this paper at https://doi.org/10.1038/s41467-018-04173-0.

**Competing interests:** The authors declare no competing interests.

**Reprints and permission** information is available online at http://npg.nature.com/reprintsandpermissions/

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_9.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2018