Low Antarctic continental climate sensitivity due to high ice sheet orography

The Antarctic continent has not warmed in the last seven decades, despite a monotonic increase in the atmospheric concentration of greenhouse gases. In this paper, we investigate whether the high orography of the Antarctic ice sheet (AIS) has helped delay warming over the continent. To that end, we contrast the Antarctic climate response to CO2-doubling with present-day orography to the response with a flattened AIS. To corroborate our findings, we perform this exercise with two different climate models. We find that, with a flattened AIS, CO2-doubling induces more latent heat transport toward the Antarctic continent, greater moisture convergence over the continent and, as a result, more surface-amplified condensational heating. Greater moisture convergence over the continent is made possible by flattening of moist isentropic surfaces, which decreases humidity gradients along the trajectories on which extratropical poleward moisture transport predominantly occurs, thereby enabling more moisture to reach the pole. Furthermore, the polar meridional cell disappears when the AIS is flattened, permitting greater CO2-forced warm temperature advection toward the Antarctic continent. Our results suggest that the high elevation of the present AIS plays a significant role in decreasing the susceptibility of the Antarctic continent to CO2-forced warming.


INTRODUCTION
The impacts of anthropogenic greenhouse gas loading of the atmosphere are global, but nowhere are these changes more consequential than over the polar regions. Changes in polar climate have important impacts globally: ice sheet melt and polar ocean warming result in sea level rise 1 ; warming and shifts in sea ice seasonality impact keystone species, like the Southern Ocean krill, which, in turn, impact the productivity of fisheries worldwide 2 ; surface-amplified polar warming can shift the jet stream and storm tracks, effectively altering midlatitude climate variability 3 ; and hemispherically asymmetric changes in polar ice cover, such as during the last glacial period, impact the position and intensity of the tropical rain belts 4 .
The polar regions, however, present a conundrum for understanding present-day climate change, as recent polar warming differs markedly between the Arctic and Antarctic. Over the past several decades, the Arctic has warmed rapidly 5 , a phenomenon linked to sea ice decline 6,7 and associated radiative feedbacks [8][9][10] . Indeed, the central Arctic has warmed by nearly 6 K over the last three decades 11 . Over the same period, Antarctic sea ice area has modestly expanded 12 , and warming has been nearly non-existent over much of the Antarctic ice sheet (AIS) ( Fig. 1; also see ref. 13 ).
The orography of the AIS, which towers nearly 4 km above sea level at its highest, is possibly the most obvious factor which could account for weak (or non-existent) warming over the Antarctic continent. The presence of the AIS has a substantial impact on the mean state dynamics and thermodynamics of the Southern Hemisphere (see, e.g., [14][15][16][17][18][19] ). Global climate model (GCM) experiments in which Antarctic continental orography is flattened exhibit weaker baroclinicity over the Southern Ocean, greater baroclinicity over the Antarctic continent, and, consequently, more frequent incursion of midlatitude eddies over the Antarctic plateau 14,20,21 . At the same time, (equatorward) katabatic flow away from the Antarctic continent ceases when Antarctic orography is flattened, and the Southern hemispheric polar cell vanishes 15,16 .
Taken together, these factors imply a greater (moist and dry) atmospheric energy transport toward the Antarctic continent when orography is flattened 22 . In a warmer CO 2 -forced world, such factors would be expected to further enhance energy transport toward the Antarctic continent, as a consequence of increased atmospheric sensible and latent heat content. These considerations suggest that the presence or absence of Antarctic orography may considerably affect the local Antarctic continental warming in response to increased CO 2 .
Recently, a single coarse-resolution GCM was used to investigate the impact of AIS orography onwarming over the Antarctic continent with CO 2 -doubling 23 . While reporting that the CO 2forced lapse rate feedback over the Antarctic continent increased significantly when Antarctic orography was flattened, implying an increase in surface-amplified warming, that study did not address which dynamical mechanisms may trigger such an increase, nor the robustness of that increase in the lapse rate feedback.
In the present study, we address these questions. We use two state-of-the-art GCMs-with different physical parameterizations and different equilibrium climate sensitivities-to show that when CO 2 is doubled, the orography of the present-day AIS shields the Antarctic continent from dynamic extrapolar transport processes that might otherwise warm it. When the orography of the AIS is flattened, we show that these shielding mechanisms are no longer at play, and winter season surface-amplified warming with CO 2doubling is significantly greater. As expected, we find that such a dynamically driven increase in surface-amplified warming augments the lapse rate feedback over the Antarctic continent. Our results imply that the high orography of Antarctica constitutes a significant local dynamical barrier to CO 2 -induced surface warming, both today and in the near future, provided orographic changes remain small.

Antarctic response to CO 2 -doubling with flattened orography
The key finding of this study is that in response to atmospheric CO 2 -doubling, warming over the Antarctic continent is significantly greater with flattened orography than with present-day orography in both CESM1.1 and CCSM4.0 ( Fig. 2a; compare red and blue bars). Warming is greater in the annual mean (0.9 K greater in CCSM4.0 and 0.7 K greater in CESM1.1) because of significantly greater warming over winter (June-July-August (JJA); 1.5 K and 1.2 K greater warming in CCSM4.0 and CESM1.1, respectively) and the shoulder seasons (March-April-May (MAM) and September-October-November (SON); 1.4 K and 1.1 K greater warming in MAM in CCSM4.0 and CESM1.1, respectively, and 0.8 K and 0.4 K greater warming in SON in CCSM4.0 and CESM1.1, respectively). In both models, CO 2 -induced warming in the present-day and flattened orography runs is comparable in summer (DJF).
The amplified warming over the Antarctic continent with flattened orography is also evident in Fig. 2b, where we compare Antarctic amplification, defined as the ratio between surface warming over the Antarctic continent and warming over the globe, in the present-day and flattened Antarctic orography experiments. In the annual mean, we find greater amplification of CO 2 -forced warming over the Antarctic continent with flattened orography than with present-day orography: the Antarctic amplification ratio increases from 1.5 to 1.9 in CESM1.1, and from 1.4 to 1.5 in CCSM4.0. Seasonally, greater Antarctic amplification with flattened orography is most prominent in winter (JJA), when the Antarctic amplification ratio increases from 1.6 to 2.2 in CESM1.1, and from 1.4 to 1.8 in CCSM4.0. Increased CO 2 -forced Antarctic amplification with flattened Antarctic orography is also evident in the shoulder seasons, though it is stronger in fall (MAM; increases from 1.7 to 2.3 in CESM1.1, and from 1.5 to 1.8 in CCSM4.0) and weaker in spring (SON; increases from 1.4 to 1.7 in CESM1.1, and from 1.3 to 1.4 in CCSM4.0). We note that while both models are in qualitative agreement, there are considerable quantitative discrepancies, implying that single model studies must be treated with caution.
Since it is clear that both greater Antarctic amplification and greater CO 2 -forced warming over the Antarctic continent with flattened orography are most prominent in winter (JJA), for the remainder of this study we focus on this season, as we analyze the Antarctic continental warming and bring to light the mechanisms responsible for it.
In Fig. 3, we present the spatial pattern of winter season Antarctic surface warming in CESM1.1 and CCSM4.0. In both models, increased CO 2 -forced surface warming with flattened Antarctic orography is greatest over the West Antarctic Ice Sheet and Peninsula, and smallest over regions of the East Antarctic. In CESM1.1, West Antarctic warming extends into the Western Pacific Ocean Sector of the East Antarctic, while in CCSM4.0, warming extends in the opposite direction into the Weddell Sea Sector. Over some areas of the West Antarctic, winter season warming is nearly 5 K greater-in both models-when Antarctic orography is flattened. The spatial pattern of Antarctic surface warming in the winter season resembles that in the annual mean (compare Fig. 3 to Supplementary Fig. 1), suggesting that winter season processes drive the annual mean signal.
To understand why CO 2 -forced warming over the Antarctic continent is greater with flattened orography than with presentday orography, we now examine how (moist and dry) transport processes into the Antarctic respond differently to CO 2 -doubling in the two cases.
Increased moisture advection with flattened orography It is well-known that when forced with CO 2 , nearly all climate models respond with an increase in poleward latent heat transport 24 : this is clearly seen in our runs with CESM1.1 and CCSM4.0 (Fig. 4, solid blue lines). Such increased poleward latent  heat transport is expected on theoretical grounds, since atmospheric specific humidity depends exponentially on temperature through the Clausius-Clapeyron relation, which enhances the meridional gradient of precipitable water in a warmer world 25 . Conversely, poleward dry static energy transport declines in many models with CO 2 forcing 24 , including CESM1.1 and CCSM4.0 (Fig. 4, pink lines), since polar-amplified warming weakens the meridional temperature gradient.
When Antarctic orography is flattened, we find a larger increase in CO 2 -forced latent heat transport toward the Antarctic continent than when orography is at its present-day height. This is evident in both CESM1.1 and CCSM4.0 (Fig. 4, panels a and b, respectively; compare solid and dashed blue lines), corresponding to a 0.05 PW increase in latent heat transport across 70S in both models; this represents a 100% increase in the poleward latent heat transport response to CO 2 -doubling in CCSM4.0, and a 50% increase in CESM1.1.
To understand why latent heat transport at high Southern latitudes increases more with CO 2 forcing when Antarctic orography is flattened, we consider the spatial structure of the moist isentropes (isolines of both the equivalent potential temperature, θ E , and the MSE) in both flattened and present-day Antarctic orography experiments ( Fig. 5). In the latitude-height plane, flattened Antarctic orography is associated with flattened moist isentropes over the pole, compared to present-day Antarctic orography where isentropes are more sloped (compare, for example, the 270 K moist isentrope in Fig. 5, panels a, c with those in panels b, d); this is true whether atmospheric CO 2 is at preindustrial concentrations (  In the extratropics, moist isentropes effectively act as isolines of the poleward moisture transport streamfunction 26,27 , particularly for moisture transport to the Antarctic continent 28 . The moist static energy (MSE), the sum of atmospheric latent heat, sensible heat, and geopotential, is conserved along a moist isentrope, and can be written as where L v is the latent heat of fusion, q is the specific humidity, C p is the specific heat of dry air at constant pressure, T is the temperature, g is the acceleration due to gravity, and Φ is the height relative to the geoid. When moist isentropes are flatter (as when Antarctic orography is flattened), air parcels gain less geopotential as they flow poleward, and therefore retain more moisture. This is also evident in the change in specific humidity with CO 2 -doubling (Fig. 5, colors): the (poleward) gradient of this specific humidity change along a moist isentropic surface, ∇(Δq)| θE , is greater with present-day Antarctic orography (panels a, c) than with flattened Antarctic orography (panels b, d), indicating that more moisture is lost along this trajectory with present-day orography than with flattened orography (compare, for example, the gradient of the change in specific humidity, denoted by the blue colors, along the 270 K moist isentrope in Fig. 5a, c with those in 5b, d).
As poleward latent heat transport towards the Antarctic continent increases more in response to CO 2 -doubling when Antarctic orography is flattened (see Fig. 4), condensational heating in the boundary layer over Antarctica also increases more: this is seen by comparing the left and right columns in Fig. 6, where condensational heating has been calculated in the lowest 200 hPa of the atmospheric column. Flattening of moisture transport trajectories (by flattening moist isentropes, see Fig. 6) results in less condensational depletion of poleward-moving moisture at the edge of the continent, thereby permitting more moisture to condense over the continental interior. Similarly, when Antarctic orography is flattened, CO 2 -forced precipitation increases are smaller over the edge of the Antarctic continent, and significantly greater over the continental interior (see Supplementary Fig. 2). In both CCSM4.0 and CESM1.1 flattened orography experiments, the increase in condensational heating in the boundary layer with CO 2 -doubling approximately occurs over the same regions where CO 2 -forced temperature increases are greatest (particularly over the West Antarctic in both models; compare Fig. 3a, b to Fig. 6a, b), indicating that condensational heating plays a key role in warming these regions.
Increased sensible heat advection with flattened orography While greater moisture transport amplifies CO 2 -forced warming over the Antarctic continent when Antarctic orography is flattened, dry processes also contribute to greater warming. Previous studies on how flattened Antarctic orography impacts the mean state climate suggest a range of local changes in the atmospheric circulation, including weakening of the polar vortex 14,17 , slowing of the eddy-driven jet 14,16,22 , increased baroclinicity closer to the pole 14,21 and a cessation of katabatic flow over the continent 15,16,18 . In both CCSM4.0 and CESM1.1, the latter is associated with the disappearance of the polar cell in the mass overturning streamfunction when orography is flattened (Fig. 7, compare contours in panels a, c with b, d; also see refs. 15,22 ).
The presence (or absence) of the polar cell determines how meridional temperature advection responds to CO 2 -doubling over the high Southern latitudes. When Antarctic orography is at present-day height, the temperature advection response is northward in the lower troposphere (Fig. 7a, c, colors), as warm temperature anomalies are advected away from the Antarctic continent by the lower branch of the polar cell. When Antarctic orography is flattened, on the other hand, the meridional temperature advection response is towards the continent (i.e., southward; Fig. 7b, d, colors) in the lower troposphere, as warmer boundary layer air from areas where sea ice has retreated is more readily advected poleward when the polar cell is absent.
Previous studies have shown that CO 2 -forced GCMs simulate similar advection of lower tropospheric warm anomalies from newly ice-free ocean regions in the Arctic to adjacent land regions: while summer season sea ice decline and ocean mixed layer warming in the Arctic directly cause (delayed) surface-amplified warming over the Arctic Ocean in fall and winter 7 , warming over boreal land regions is primarily through warm air advection from the polar oceans 29 . Therefore, when Antarctic orography is flattened, there is greater similarity between the poles in how advection of boundary layer air from newly ice-free ocean waters warms land areas in winter.
Greater surface-amplified warming with flattened orography As a result of both the dry and moist transport processes described above, Antarctic continent surface temperatures warm more with CO 2 -doubling when Antarctic orography is flattened (recall Figs. 2 and 3).
The overlying atmosphere also warms with CO 2 -doubling, and that warming is also more surface-amplified with flattened Antarctic orography than with present-day orography (Fig. 8, compare red with black lines). As we have shown, such surfaceamplified warming is a result of greater atmospheric energy transport convergence (both dry and moist) into the lower troposphere when Antarctic orography is flattened, combined with strong atmospheric static stability in winter which prevents  Fig. 4; also see ref. 30 ).
Greater surface-amplified warming with CO 2 -doubling when Antarctic orography is flattened also leads to a stronger (positive) lapse rate feedback over the Antarctic continent itself: 1.2 W/m 2 /K (0.8 W/m 2 /K) with flattened Antarctic orography compared to 0.5 W/m 2 /K (0.3 W/m 2 /K) with present-day Antarctic orography in CESM1.1 (CCSM4.0), which we have computed using the radiative kernel method 31 with a CAM3-derived temperature kernel 32 . The positive lapse rate feedback in the high latitudes is one of the most important radiative factors responsible for polaramplified warming in response to CO 2 forcing, particularly that of the Arctic 9,10,30,33,34 . Its increase when Antarctic orography is flattened is consistent with the increase in Antarctic (continental) amplification reported here (recall Fig. 2).

DISCUSSION
In this study, we have used two state-of-the-art GCMs to show that Antarctic ice sheet orography robustly decreases climate sensitivity over the Antarctic continent, and that a flattened Antarctic ice sheet would experience significantly greater surface warming than the present-day Antarctic ice sheet. When CO 2 is doubled, we have shown that both moist and dry dynamic processes are more efficient at warming a flattened Antarctic continent. These dynamic processes differ from those that occur when Antarctica is at its present-day elevation because the presence of Antarctic orography significantly alters the atmospheric circulation at high Southern latitudes (see, e.g., [14][15][16]; differences in the unperturbed circulation due to differences in orography result in different (moist and dry transport) responses to CO 2 -doubling.
Over the Antarctic continent, our results generally agree with those of ref. 23 , who also found greater surface-amplified warming over the Antarctic continent (and, therefore, a more positive lapse rate feedback) with CO 2 -doubling when Antarctic orography is flattened. In our study, we have specifically identified a number of atmospheric processes which act to increase CO 2 -forced surfaceamplified warming (and, hence, the lapse rate feedback) over the Antarctic continent when its orography is flattened. We suggest that these mechanisms are robust because they depend upon changes in atmospheric circulation with Antarctic orography flattening that are consistent across a broad range of GCMs of varying complexity (see, e.g., 14,17,19,21,22 ). Unlike 23 , we do not find that annual mean Antarctic surface warming with CO 2 -doubling exceeds that over the Arctic when Antarctic orography was flattened in either of the models we studied (see Supplementary Fig. 3). In fact, our two models disagree on the impact of Antarctic orography flattening on Arctic warming: CCSM4.0 shows more Arctic warming with CO 2 -doubling when Antarctic orography is flattened, while CESM1.1 shows less. Therefore, the Arctic response to a flattened AIS does not appear to be robust. In addition, we find that the remote (atmospheric and oceanic) energy transport response to CO 2 -doubling when Antarctic orography is flattened, reported by 23 , is also not robust (see Supplementary Fig. 4). Specifically, CCSM4.0 shows a greater increase in ocean heat transport into the Arctic with CO 2 -doubling when Antarctic orography is flattened, whereas CESM1.1 shows a smaller increase; these transport responses are consistent with increased Arctic warming in the former and decreased Arctic Antarctic orography is at its present-day height and b, d when Antarctic orography is flattened to 10% of its present-day height. The contours of moist entropy are shown for the corresponding pre-industrial experiment (i.e., PI or PI_FA; in black) and the CO 2 -doubling experiment (i.e., 2 × CO 2 or 2 × CO 2 _FA; in red).
In this study, we've identified Antarctic orography as one of (likely) several factors that reduces the magnitude of CO 2 -forced warming over the Antarctic continent. Further factors may also be responsible for a weak climate change signal over the Antarctic continent, including heat uptake by the Southern Ocean 37 and the meagre decline of Antarctic sea ice (in both observations and future GCM projections; see refs. 12,38 ). More research is necessary to identify the extent to which warming or cooling over the open Southern Ocean and sea ice zones impacts surface temperatures over high-elevation regions of the Antarctic ice sheet, as the dynamics of the extratropical atmosphere suggests that these regions should be relatively isolated from adjacent low-elevation regions 28,39,40 . Our experiments, for example, show that differences in sea ice loss and warming over the marginal ice zone are not primarily responsible for greater warming over the flattened Antarctic continent. In fact, CESM1.1 simulates less Antarctic sea ice retreat with CO 2 -doubling when Antarctica is flattened whereas CCSM4.0 simulates more Antarctic sea ice retreat (see Supplementary Fig. 5), even though both models simulate greater warming over the continent (recall Figs. 2 and 3).
For the sake of completeness, we also remind the reader that other factors have been suggested to have contributed to the absence of warming over the Antarctic continent in the past several decades, but these suggestions have been shown to be off target. For instance, on the basis of idealized model runs 41 , it was suggested that the formation of the ozone hole over the South pole, which has caused the majority of the positive trends in the Southern Annual Mode with accompanying surface wind stress changes 42 , could have contributed to colder SSTs and increased sea ice extent around Antarctica 43 . However, numerous studies with realistic models have convincingly shown that ozone depletion (which increases the amount of shortwave radiation reaching the surface) has in fact contributed to warmer SSTs and decreased sea ice extent [44][45][46][47][48][49] .
Other studies have highlighted the existence of a negative greenhouse effect over the Antarctic continent: the instantaneous outgoing longwave radiation at the top-of-atmosphere increases, rather than decreases, with higher levels of atmospheric CO 2 50-52 . That negative greenhouse effect, which only occurs in some months of the year, owes its existence to the tropospheric temperature inversion over the extremely cold Antarctic surface, and is enhanced by the relative absence of tropospheric water vapor in that region 53 . While that negative greenhouse effect has indeed been observed from satellites (as documented by 53 ), it dissipates rapidly following abrupt CO 2 -quadrupling in fully coupled GCMs due to fast stratospheric adjustments 54 , and does not result in a cooling of the Antarctic surface. Furthermore, Flanner et al. 55 show that the net surface longwave radiative impact of greenhouse gases will always tend to heat the surface at high latitudes because of the local temperature inversion, regardless of whether the greenhouse effect is positive or negative at the top-of-atmosphere. In the present study, we do find that the net (downward) surface longwave flux with CO 2doubling is greater when Antarctic orography is flattened. However, we are leery to attribute the surface-amplified warming with flat orography to this factor: analysis of surface radiative kernels indicates that anomalies in the downward longwave flux at the surface primarily arise as a consequence of surface temperature anomalies, rather being the cause of those anomalies 56 . A thorough assessment of instantaneous radiative forcing, and of the accompanying rapid adjustments, is outside the scope of the present study, as we are here interested in the Antarctic surface climate response to CO 2 -doubling at quasi-equilibrium, not in the details of the radiative forcing and adjustment immediately following the doubling of CO 2 .
We conclude by acknowledging some important caveats on the results presented here. First, we have only used two models in our study, and a more extensive model intercomparison may identify additional mechanisms that contribute to a reduction of Antarctic continental warming under CO 2 -forcing. Furthermore, we reiterate that although models generally agree on the local dynamic atmospheric impacts of Antarctic orography flattening, there is little intermodel agreement regarding the remote impacts of Antarctic orography flattening. For example, some GCMs simulate a global mean cooling when Antarctic orography is flattened 22 , while other GCMs simulate a global mean warming 57 , suggesting a large intermodel spread in the global radiative forcing that accompanies Antarctic orography changes. We expect that major difference across models in the global radiative forcing associated with flattening of Antarctic orography may also result in differences in the global CO 2 -doubling response when Antarctic orography is flattened, as we have noted above. Despite these caveats, we believe that the mechanisms we proposed in this study are likely to be robust, in that they depend on changes in local atmospheric dynamics circa the Antarctic continent, which are known to be consistent across a range of models.

GCM experiments
We utilize two state-of-the-art, fully-coupled GCMs: the Community Climate System Model, version 4.0 (CCSM4.0; see ref. 58 ), and the Community Earth System Model, version 1.1 (CESM1.1; see ref. 59 ). The primary difference between the two GCMs is in the atmosphere model component: CCSM4.0 utilizes version 4 of the Community Atmosphere Model (CAM4 60 ) at 2°spatial resolution, while CESM1.1 utilizes version 5 Fig. 7 Mass overturning streamfunction and temperature advection change with CO 2 -doubling. Mean meridional mass overturning streamfunction (contours; in kg per sec) and change in temperature advection with CO 2 -doubling (colors; in Km s −1 ) in the a, b CESM1.1 and c, d CCSM4.0 when a, c Antarctic orography is at its present-day height, 2 × CO 2 minus PI, and when b, d Antarctic orography is flattened to 10% of its present-day height, 2 × CO 2 _FA minus PI_FA. The mass overturning streamfunction is shown for the corresponding pre-industrial experiment (PI or PI_FA).
(CAM5 61 ) at 1°spatial resolution. CAM4 and CAM5 implement different atmospheric cloud, convection, boundary layer (including gravity wave drag), and precipitation parameterizations; as a result, the mean state climate and climate sensitivity of these two models differ significantly 62 .
Other than the atmosphere model, all other components are nearly identical between CCSM4.0 and CESM1.1. Both models utilize the fully dynamic CICE4 sea ice component 63 and the POP2 ocean component 64 with an isopycnal mixing parameterization to capture transport from subgrid scale eddies 65 . Both ocean and sea ice are run on the same (nominally) 1°spatial resolution grid with the North Pole singularity centered over Greenland. Also, the land component for both CCSM4.0 and CESM1.1 is the CLM4 66 , which is run at the same resolution as the respective atmosphere model (i.e., 2°for CCSM4.0 and 1°for CESM1.1).
For both models, we perform a set of four runs: (1) a pre-industrial control experiment with present-day Antarctic orography (PI); (2) a preindustrial control experiment with flattened Antarctic orography (PI_FA); (3) a CO 2 -doubling experiment with present-day Antarctic orography (2 × CO 2 ); and (4) a CO 2 -doubling experiment with flattened Antarctic orography (2 × CO 2 _FA). For both models, the pre-industrial control experiment with present-day Antarctic orography (PI) was obtained from a millennium-length, publicly available, equilibrated run (see Table 1, Component Set and Configuration). Each of the other experiments were branched from this equilibrated pre-industrial control run, and integrated for a total of 230 model years with the modifications noted: CO 2 -doubling (atmospheric CO 2 concentrations abruptly doubled from 284 to 568 ppm at the beginning of the experiment), Antarctic orography flattening (reduction of Antarctic orography to 10% of its present-day elevation at the beginning of the experiment), or both CO 2 -doubling and orography flattening. The details of each experiment are shown in Table 1.
Annual and seasonal climatologies are calculated using the final 30 years of each experiment. The net top-of-atmosphere anomalies for all experiments are given in Supplementary Table 1. For each model, the CO 2 -forced response is computed as the difference between the 2 × CO 2 and PI runs, either with present-day Antarctic orography (i.e., 2 × CO 2 minus PI) or with flattened Antarctic orography (i.e., 2 × CO 2 _FA minus PI_FA); the former shows the climate response to CO 2 -doubling when Antarctic orography is at its present-day elevation, and the latter shows the climate response to CO 2 -doubling with flattened Antarctic orography.

DATA AVAILABILITY
Global climate model data generated for this study are available at https://doi.org/ 10.5281/zenodo.3998429.