Your search

Abstract The Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6) is the primary effort of CMIP6 (Coupled Model Intercomparison Project?Phase 6) focusing on ice sheets, designed to provide an ensemble of process-based projections of the ice-sheet contribution to sea-level rise over the twenty-first century. However, the behavior of the Antarctic Ice Sheet beyond 2100 remains largely unknown: several instability mechanisms can develop on longer time scales, potentially destabilizing large parts of Antarctica. Projections of Antarctic Ice Sheet evolution until 2300 are presented here, using an ensemble of 16 ice-flow models and forcing from global climate models. Under high-emission scenarios, the Antarctic sea-level contribution is limited to less than 30 cm sea-level equivalent (SLE) by 2100, but increases rapidly thereafter to reach up to 4.4 m SLE by 2300. Simulations including ice-shelf collapse lead to an additional 1.1 m SLE on average by 2300, and can reach 6.9 m SLE. Widespread retreat is observed on that timescale in most West Antarctic basins, leading to a collapse of large sectors of West Antarctica by 2300 in 30%?40% of the ensemble. While the onset date of retreat varies among ice models, the rate of upstream propagation is highly consistent once retreat begins. Calculations of sea-level contribution including water density corrections lead to an additional ?10% sea level and up to 50% for contributions accounting for bedrock uplift in response to ice loading. Overall, these results highlight large sea-level contributions from Antarctica and suggest that the choice of ice sheet model remains the leading source of uncertainty in multi-century projections.

Oceanic mesoscale eddy mixing plays a crucial role in Earth’s climate system by redistributing heat, salt, and carbon. For many ocean and climate models, mesoscale eddies still need to be parameterized. This is often done via an eddy diffusivity K , which sets the strength of turbulent downgradient tracer fluxes. A well-known effect is the modulation of K in the presence of background potential vorticity (PV) gradients, which suppresses cross-PV gradient mixing. Topographic slopes can induce such suppression through topographic PV gradients. However, this effect has received little attention, and topographic effects are often not included in parameterizations for K . In this study, we show that it is possible to describe the effect of topography on K analytically in a barotropic framework, using a simple stochastic representation of eddy–eddy interactions. We obtain an analytical expression for the depth-averaged K as a function of the bottom slope, which we validate against diagnosed eddy diffusivities from a numerical model. The obtained analytical expression can be generalized to any constant barotropic PV gradient. Moreover, the expression is consistent with empirical parameterizations for eddy diffusivity over topography from previous studies and provides a physical rationalization for these parameterizations. The new expression helps to understand how eddy diffusivities vary across the ocean, and thus how mesoscale eddies impact ocean mixing processes.

During the mid-Pliocene warm period (mPWP; 3.264–3.025 Ma), atmospheric CO2 concentrations were approximately 400 ppm, and the Antarctic Ice Sheet was substantially reduced compared to today. Antarctica is surrounded by the Southern Ocean, which plays a crucial role in the global oceanic circulation and climate regulation. Using results from the Pliocene Model Intercomparison Project (PlioMIP2), we investigate Southern Ocean conditions during the mPWP with respect to the pre-industrial period. We find that the mean sea surface temperature (SST) warming in the Southern Ocean is 2.8 °C, while global mean SST warming is 2.4 °C. The enhanced warming is strongly tied to a dramatic decrease in sea ice cover over the mPWP Southern Ocean. We also see a freshening of the ocean (sub)surface, driven by an increase in precipitation over the Southern Ocean and Antarctica. The warmer and fresher surface leads to a highly stratified Southern Ocean that can be related to weakening of the deep abyssal overturning circulation. Sensitivity simulations show that the decrease in sea ice cover and enhanced warming is largely a consequence of the reduction in the Antarctic Ice Sheet. In addition, the mPWP geographic boundary conditions are responsible for approximately half of the increase in mPWP SST warming, sea ice loss, precipitation, and stratification increase over the Southern Ocean. From these results, we conclude that a strongly reduced Antarctic Ice Sheet during the mPWP has a substantial influence on the state of the Southern Ocean and exacerbates the changes that are induced by a higher CO2 concentration alone. This is relevant for the long-term future of the Southern Ocean, as we expect melting of the western Antarctic Ice Sheet in the future, an effect that is not currently taken into account in future projections by Coupled Model Intercomparison Project (CMIP) ensembles.

Polar warming, ice melt and strong precipitation events are strongly affected by episodic poleward advection of warm and moist air (Woods and Caballero 2016 J. Clim. 29 4473–85; Wille et al 2019 Nat. Geosci. 12 911–6), which, in turn, is linked to variability in poleward moisture transport (PMT) (Nash et al 2018 J. Geophys. Res. Atmos. 123 6804–21). However, processes governing regional impacts of PMT as well as long-term trends remain largely unknown. Here we use an ensemble of state-of-the-art global climate models in standardized scenario simulations (1850–2100) to show that both the Arctic and the Antarctic exhibit distinct geographical patterns of PMT-related warming. Specifically, years with high PMT experience considerable warming over subarctic Eurasia and West-Antarctica (Raphael et al 2016 Bull. Am. Meteorol. Soc. 97 111–21), whereas precipitation is distributed more evenly over the polar regions. The warming patterns indicate preferred routes of atmospheric rivers (Woods and Caballero 2016 J. Clim. 29 4473–85), which may regionally enhance atmospheric moisture content, cloud cover, and downward longwave radiative heating in years with comparatively high PMT (Scott et al 2019 J. Clim. 32 665–84). Trend-analyses reveal that the link between PMT-variability and regional precipitation patterns will weaken in both polar regions. Even though uncertainties associated with intermodel differences are considerable, the advection of warm and moist air associated with PMT-variability is likely to increasingly cause mild conditions in both polar regions, which in the Arctic will reinforce sea-ice melt. Similarly, the results suggest that warm years in West-Antarctica disproportionally contribute to ice sheet melt (Trusel et al 2015 Nat. Geosci. 8 927–32), enhancing the risk of ice-sheet instabilities causing accelerated and sudden sea-level rise.

Abstract The additional water from the Antarctic ice sheet and ice shelves due to climate-induced melt can impact ocean circulation and global climate. However, the major processes driving melt are not adequately represented in Coupled Model Intercomparison Project phase 6 (CMIP6) models. Here, we analyze a novel multi-model ensemble of CMIP6 models with consistent meltwater addition to examine the robustness of the modeled response to meltwater, which has not been possible in previous single-model studies. Antarctic meltwater addition induces a substantial weakening of open-ocean deep convection. Additionally, Antarctic Bottom Water warms, its volume contracts, and the sea surface cools. However, the magnitude of the reduction varies greatly across models, with differing anomalies correlated with their respective mean-state climatology, indicating the state-dependency of the climate response to meltwater. A better representation of the Southern Ocean mean state is necessary for narrowing the inter-model spread of response to Antarctic meltwater.

The Antarctic Ice Sheet represents the largest source of uncertainty in future sea level rise projections, with a contribution to sea level by 2100 ranging from −5 to 43 cm of sea level equivalent under high carbon emission scenarios estimated by the recent Ice Sheet Model Intercomparison for CMIP6 (ISMIP6). ISMIP6 highlighted the different behaviors of the East and West Antarctic ice sheets, as well as the possible role of increased surface mass balance in offsetting the dynamic ice loss in response to changing oceanic conditions in ice shelf cavities. However, the detailed contribution of individual glaciers, as well as the partitioning of uncertainty associated with this ensemble, have not yet been investigated. Here, we analyze the ISMIP6 results for high carbon emission scenarios, focusing on key glaciers around the Antarctic Ice Sheet, and we quantify their projected dynamic mass loss, defined here as mass loss through increased ice discharge into the ocean in response to changing oceanic conditions. We highlight glaciers contributing the most to sea level rise, as well as their vulnerability to changes in oceanic conditions. We then investigate the different sources of uncertainty and their relative role in projections, for the entire continent and for key individual glaciers. We show that, in addition to Thwaites and Pine Island glaciers in West Antarctica, Totten and Moscow University glaciers in East Antarctica present comparable future dynamic mass loss and high sensitivity to ice shelf basal melt. The overall uncertainty in additional dynamic mass loss in response to changing oceanic conditions, compared to a scenario with constant oceanic conditions, is dominated by the choice of ice sheet model, accounting for 52 % of the total uncertainty of the Antarctic dynamic mass loss in 2100. Its relative role for the most dynamic glaciers varies between 14 % for MacAyeal and Whillans ice streams and 56 % for Pine Island Glacier at the end of the century. The uncertainty associated with the choice of climate model increases over time and reaches 13 % of the uncertainty by 2100 for the Antarctic Ice Sheet but varies between 4 % for Thwaites Glacier and 53 % for Whillans Ice Stream. The uncertainty associated with the ice–climate interaction, which captures different treatments of oceanic forcings such as the choice of melt parameterization, its calibration, and simulated ice shelf geometries, accounts for 22 % of the uncertainty at the ice sheet scale but reaches 36 % and 39 % for Institute Ice Stream and Thwaites Glacier, respectively, by 2100. Overall, this study helps inform future research by highlighting the sectors of the ice sheet most vulnerable to oceanic warming over the 21st century and by quantifying the main sources of uncertainty.

As the climate warms, the grounded ice sheet and floating ice shelves surrounding Antarctica are melting and releasing additional freshwater into the Southern Ocean. Nonetheless, almost all existing coupled climate models have fixed ice sheets and lack the physics required to represent the dominant sources of Antarctic melt. These missing ice dynamics represent a key uncertainty that is typically unaccounted for in current global climate change projections. Previous modelling studies that have imposed additional Antarctic meltwater have demonstrated regional impacts on Southern Ocean stratification, circulation, and sea ice, as well as remote changes in atmospheric circulation, tropical precipitation, and global temperature. However, these previous studies have used widely varying rates of freshwater forcing, have been conducted using different climate models and configurations, and have reached differing conclusions on the magnitude of meltwater–climate feedbacks. The Southern Ocean Freshwater Input from Antarctica (SOFIA) initiative brings together a team of scientists to quantify the climate system response to Antarctic meltwater input along with key aspects of the uncertainty. In this paper, we summarize the state of knowledge on meltwater discharge from the Antarctic ice sheet and ice shelves to the Southern Ocean and explain the scientific objectives of our initiative. We propose a series of coupled and ocean–sea ice model experiments, including idealized meltwater experiments, historical experiments with observationally consistent meltwater input, and future scenarios driven by meltwater inputs derived from stand-alone ice sheet models. Through coordinating a multi-model ensemble of simulations using a common experimental design, open data archiving, and facilitating scientific collaboration, SOFIA aims to move the community toward better constraining our understanding of the climate system response to Antarctic melt.

Clouds have a large effect on the radiation budget and represent a major source of uncertainty in climate models. Supercooled liquid clouds can exist at temperatures as low as 235 K, and the radiative effect of these clouds depends on the complex refractive index (CRI) of liquid water. Laboratory measurements have demonstrated that the liquid-water CRI is temperature-dependent, but corroboration with field measurements is difficult. Here we present measurements of the downwelling infrared radiance and in-situ measurements of supercooled liquid water in a cloud at temperatures as low as 240 K, made at South Pole Station in 2001. These results demonstrate that including the temperature dependence of the liquid-water CRI is essential for accurate calculations of radiative transfer through supercooled liquid clouds. Furthermore, we show that when cloud properties are retrieved from infrared radiances (using the spectral range 500–1,200 cm−1) spurious ice may be retrieved if the 300 K CRI is used for cold liquid clouds (∼240 K). These results have implications for radiative transfer in climate models as well as for retrievals of cloud properties from infrared radiance spectra.

Antarctic sea ice has paradoxically become more extensive over the past four decades despite a warming climate. The regional expression of this trend has been linked to changes in vertical redistribution of ocean heat and large-scale wind-field shifts. However, the short length of modern observations has hindered attempts to attribute this trend to anthropogenic forcing or natural variability. Here, we present two new decadal-resolution records of sea ice and sea surface temperatures that document pervasive regional climate heterogeneity in Indian Antarctic sea-ice cover over the last 2,000 years. Data assimilation of our marine records in a climate model suggests that the reconstructed dichotomous regional conditions were driven by the multi-decadal variability of the El Niño Southern Oscillation and Southern Annular Mode (SAM). For example, during an El Niño/SAM– combination, the northward sea-ice transport was reduced while heat advection from the subtropics to the Southern Ocean increased, which resulted in reduced sea-ice extent in the Indian sector as sea ice was compacted along the Antarctic coast. Our results therefore indicate that natural variability is large in the Southern Ocean and suggest that it has played a crucial role in the recent sea-ice trends and their decadal variability in this region.

Open-ocean polynyas effectively couple the ocean and atmosphere through large ice-free areas within the sea-ice cover, release vast quantities of oceanic heat, and impact deep ocean ventilation. Changes in polynya activity, particularly in the Weddell Sea, may be key to longer time-scale climate fluctuations, feedbacks and abrupt change. While changes in the occurrence of Weddell Sea polynyas are generally attributed to changes in the atmospheric surface forcing, the role of internal ocean dynamics for polynya variability is not well-resolved. In this study we employ a global coupled ocean-sea ice model with a repeating annual atmospheric cycle to explore changes in Weddell Sea water mass properties, stratification and ocean circulation driven by open-ocean polynyas. During the 1300-year long simulation, two large polynyas occur in the central Weddell Sea. Our results suggest that Weddell polynyas may be triggered without inter-annual changes in the atmospheric forcing. This highlights the role of ocean processes in preconditioning and triggering open-ocean polynyas on multi-centennial time-scales. The simulated polynyas form due to internal ocean-sea ice dynamics associated with a slow build-up and subsequent release of subsurface heat. A strong stratification and weak vertical mixing is necessary for building the subsurface heat reservoir. Once the water column turns unstable, enhanced vertical mixing of warm and saline waters into the surface layer causes efficient sea ice melt and the polynya appears. Subsequent, vigorous deep convection is maintained through upwelling of warm deep water leading to enhanced bottom water formation. We find a cessation of simulated deep convection and polynya activity due to long-term cooling and freshening of the subsurface heat reservoir. As subsurface waters in the Southern Ocean are now becoming warmer and saltier, we speculate that larger and more persistent Weddell polynyas could become more frequent in the future.

We examine the response of the Community Earth System Model Versions 1 and 2 (CESM1 and CESM2) to abrupt quadrupling of atmospheric CO2 concentrations (4xCO2) and to 1% annually increasing CO2 concentrations (1%CO2). Different estimates of equilibrium climate sensitivity (ECS) for CESM1 and CESM2 are presented. All estimates show that the sensitivity of CESM2 has increased by 1.5 K or more over that of CESM1. At the same time the transient climate response (TCR) of CESM1 and CESM2 derived from 1%CO2 experiments has not changed significantly—2.1 K in CESM1 and 2.0 K in CESM2. Increased initial forcing as well as stronger shortwave radiation feedbacks are responsible for the increase in ECS seen in CESM2. A decomposition of regional radiation feedbacks and their contribution to global feedbacks shows that the Southern Ocean plays a key role in the overall behavior of 4xCO2 experiments, accounting for about 50% of the total shortwave feedback in both CESM1 and CESM2. The Southern Ocean is also responsible for around half of the increase in shortwave feedback between CESM1 and CESM2, with a comparable contribution arising over tropical ocean. Experiments using a thermodynamic slab-ocean model (SOM) yield estimates of ECS that are in remarkable agreement with those from fully coupled Earth system model (ESM) experiments for the same level of CO2 increase. Finally, we show that the similarity of TCR in CESM1 and CESM2 masks significant regional differences in warming that occur in the 1%CO2 experiments for each model.

To investigate the role of tides in Weddell Sea ocean-ice shelf melt interactions, and resulting consequences for ocean properties and sea ice interactions, we develop a regional ocean-sea ice model configuration, with time-varying ocean boundary and atmospheric forcing, including the deep open ocean (at 2.5–4 km horizontal resolution), the southwestern continental shelf (≈2.5 km), and the adjacent cavities of eastern Weddell, Larsen, and Filchner-Ronne ice shelves (FRIS, 1.5–2.5 km). Simulated circulation, water mass, and ice shelf melt properties compare overall well with available open ocean and cavity observational knowledge. Tides are shown to enhance the kinetic energy of the time-varying flow in contact with the ice shelves, thereby increasing melt. This dynamically driven impact of tides on net melting is to almost 90% compensated by cooling through the meltwater that is produced but not quickly exported from regions of melting in the Weddell Sea cold-cavity regime. The resulting systematic tide-driven enhancement of both produced meltwater and its refreezing on ascending branches of, especially the FRIS, cavity circulation acts to increase net ice shelf melting (by 50% in respect to the state without tides, ≈50 Gt yr−1). In addition, tides also increase the melt-induced FRIS cavity circulation, and the meltwater export by the FRIS outflow. Simulations suggest attendant changes on the open-ocean southwestern continental shelf, characterized by overall freshening and small year-round sea ice thickening, as well as in the deep southwestern Weddell Sea in the form of a marked freshening of newly formed bottom waters.

Current global warming is causing significant changes in snowfall in polar regions, directly impacting the mass balance of the ice caps. The only water supply in Antarctica, precipitation, is poorly estimated from surface measurements. The onboard cloud-profiling radar of the CloudSat satellite provided the first real opportunity to estimate solid precipitation at continental scale. Based on CloudSat observations, we propose to explore the vertical structure of precipitation in Antarctica over the 2007–2010 period. A first division of this data set following a topographical approach (continent vs. peripheral regions, with a 2,250 m topographical criterion) shows a high snowfall rate (275 mm yr at 1,200 m above ground level) with low relative seasonal variation ( ) over the peripheral areas. Over the plateau, the snowfall rate is low (34 mm yr at 1,200 m above ground level) with a much larger relative seasonal variation ( ). A second study that follows a geographical division highlights the average vertical structure of precipitation and temperature depending on the regions and their interactions with topography. In particular, over ice shelves, we see a strong dependence of the distribution of snowfall on the sea ice coverage. Finally, the relationship between precipitation and temperature is analyzed and compared with a simple analytical relationship. This study highlights that precipitation is largely dependent on the advection of air masses along the topographic slopes with an average vertical wind of 0.02 m s . This provides new diagnostics to evaluate climate models with a three-dimensional approach of the atmospheric structure of precipitation.

While observed mesospheric polar nitric acid enhancements have been attributed to energetic particle precipitation through ion cluster chemistry in the past, this phenomenon is not reproduced in current whole-atmosphere chemistry-climate models. We investigate such nitric acid enhancements resulting from energetic electron precipitation events using a recently developed variant of the Whole Atmosphere Community Climate Model (WACCM) that includes a sophisticated ion chemistry tailored for the D-layer of the ionosphere (50–90 km), namely, WACCM-D. Using the specified dynamics mode, that is, nudging dynamics in the troposphere and stratosphere to meteorological reanalyses, we perform a 1-year-long simulation (July 2009–June 2010) and contrast WACCM-D with the standard WACCM. Both WACCM and WACCM-D simulations are performed with and without forcing from medium-to-high energy electron precipitation, allowing a better representation of the energetic electrons penetrating into the mesosphere. We demonstrate the effects of the strong particle precipitation events which occurred during April and May 2010 on nitric acid and on key ion cluster species, as well as other relevant species of the nitrogen family. The 1-year-long simulation allows the event-related changes in neutral and ionic species to be placed in the context of their annual cycle. We especially highlight the role played by medium-to-high energy electrons in triggering ion cluster chemistry and ion-ion recombinations in the mesosphere and lower thermosphere during the precipitation event, leading to enhanced production of nitric acid and raising its abundance by 2 orders of magnitude from 10−4 to a few 10−2 ppb.

The establishment of a modern-like monsoon climate in East Asia by the early Miocene was a complex process forced by several factors, and previous studies paid less attention to global cooling. Here we investigate this process using climate modeling by considering changes in topography and global cooling under the early Miocene boundary conditions. Using early Miocene paleogeography and an atmospheric CO2 concentration of 560 ppmv, our model results indicate that a nonzonal climate pattern has appeared in East China but that this climate exhibits weak precipitation and wind seasonality. Such seasonality strengthens as the concentration of atmospheric CO2 decreases from 560 to 420 ppmv and resembles the modern-like condition after the growth of the Antarctic ice sheet. Although the development of East Asian topography can further strengthen this seasonality, our results indicate that global cooling is also pivotal for the establishment of a modern-like monsoon climate in East Asia.

An idealized eddy-resolving numerical model, with topographic features common to the southern Weddell Sea, is constructed to study mechanisms through which warm deep water enters a wide continental shelf with a trough. The open ocean, represented by a 1700 m deep channel, is connected to a 400 m deep shelf with a continental slope. The shelf is narrow (50 km) in the east but widens to 300 km at the center of the model domain. Over the narrow shelf, the slope front is balanced by wind-driven Ekman downwelling and counteracting eddy overturning, favoring on-shelf transport of warm water in summer scenarios when fresher surface water is present. Over the wide shelf, the Ekman downwelling ceases, and the mesoscale eddies relax the front. Inflow of warm water is sensitive to along-shelf salinity gradients and is most efficient when denser water over the wide shelf favors up-slope eddy transport along isopycnals of the V-shaped slope front. Inflow along the eastern side of the trough cannot penetrate the sill region due to potential vorticity constraints, while along the western trough flank, eddy-induced inflow crosses the sill and reaches the ice front. The warm inflow into the trough is sensitive to the density of the outflowing dense shelf water. For weaker winds, absence of the dense water outflow leads to a reversal of the trough circulation and a strong inflow of warm water, while for stronger winds, baroclinic effects become less important and the inflow is similar to experiments including dense water outflow.

We review recent progress in understanding the role of sea ice, land surface, stratosphere, and aerosols in decadal-scale predictability and discuss the perspectives for improving the predictive capabilities of current Earth system models (ESMs). These constituents have received relatively little attention because their contribution to the slow climatic manifold is controversial in comparison to that of the large heat capacity of the oceans. Furthermore, their initialization as well as their representation in state-of-the-art climate models remains a challenge. Numerous extraoceanic processes that could be active over the decadal range are proposed. Potential predictability associated with the aforementioned, poorly represented, and scarcely observed constituents of the climate system has been primarily inspected through numerical simulations performed under idealized experimental settings. The impact, however, on practical decadal predictions, conducted with realistically initialized full-fledged climate models, is still largely unexploited. Enhancing initial-value predictability through an improved model initialization appears to be a viable option for land surface, sea ice, and, marginally, the stratosphere. Similarly, capturing future aerosol emission storylines might lead to an improved representation of both global and regional short-term climatic changes. In addition to these factors, a key role on the overall predictive ability of ESMs is expected to be played by an accurate representation of processes associated with specific components of the climate system. These act as “signal carriers,” transferring across the climatic phase space the information associated with the initial state and boundary forcings, and dynamically bridging different (otherwise unconnected) subsystems. Through this mechanism, Earth system components trigger low-frequency variability modes, thus extending the predictability beyond the seasonal scale.

This article investigates the annual cycle observed in the Antarctic baseline aerosol scattering coefficient, total particle number concentration, and particle number size distribution (PNSD), as measured at Troll Atmospheric Observatory. Mie theory shows that the annual cycles in microphysical and optical aerosol properties have a common cause. By comparison with observations at other Antarctic stations, it is shown that the annual cycle is not a local phenomenon, but common to central Antarctic baseline air masses. Observations of ground-level ozone at Troll as well as backward plume calculations for the air masses arriving at Troll demonstrate that the baseline air masses originate from the free troposphere and lower stratosphere region, and descend over the central Antarctic continent. The Antarctic summer PNSD is dominated by particles with diameters <100 nm recently formed from the gas-phase despite the absence of external sources of condensible gases. The total particle volume in Antarctic baseline aerosol is linearly correlated with the integral insolation the aerosol received on its transport pathway, and the photooxidative production of particle volume is mostly limited by photooxidative capacity, not availability of aerosol precursor gases. The photooxidative particle volume formation rate in central Antarctic baseline air is quantified to 207 ±4 μm3/(MJ m). Further research is proposed to investigate the applicability of this number to other atmospheric reservoirs, and to use the observed annual cycle in Antarctic baseline aerosol properties as a benchmark for the representation of natural atmospheric aerosol processes in climate models.

An ongoing challenge in attributing anthropogenic climate change is to distinguish anthropogenic and natural changes of atmospheric composition, e.g. concerning atmospheric aerosol and its climate effects. Aerosol properties measured at pristine locations, to the extend they still exist, can serve as a climate model benchmark for verifying the representation of natural aerosol processes in the model.

We have examined changes in climate which result from the sudden termination of geoengineering after 50 years of offsetting a 1% per annum increase in CO2 concentrations by a reduction of solar radiation, as simulated by 11 different climate models in experiment G2 of the Geoengineering Model Intercomparison Project. The models agree on a rapid increase in global-mean temperature following termination accompanied by increases in global-mean precipitation rate and decreases in sea-ice cover. There is no agreement on the impact of geoengineering termination on the rate of change of global-mean plant net primary productivity. There is a considerable degree of consensus for the geographical distribution of temperature change following termination, with faster warming at high latitudes and over land. There is also considerable agreement regarding the distribution of reductions in Arctic sea-ice, but less so for the Antarctic. There is much less agreement regarding the patterns of change in precipitation and net primary productivity, with a greater degree of consensus at higher latitudes.