Your search

During the Quaternary, ice sheets experienced several retreat–advance cycles, strongly influencing climate patterns. In order to properly simulate these phenomena, it is preferable to use physics-based models instead of parameterizations to estimate the surface mass balance (SMB), which strongly influences the evolution of the ice sheet. To further investigate the potential of these SMB models, this work evaluates the BErgen Snow SImulator (BESSI), a multi-layer snow model with high computational efficiency, as an alternative to providing the SMB for the Earth system model iLOVECLIM for multi-millennial simulations as in paleostudies. We compare the behaviors of BESSI and insolation temperature melt (ITM), an existing SMB scheme of iLOVECLIM during the Last Interglacial (LIG). Firstly, we validate the two SMB models using the regional climate model Mod- èle Atmosphérique Régional (MAR) as forcing and reference for the present-day climate over the Greenland and Antarctic ice sheets. The evolution of the SMB over the LIG (130–116 ka) is computed by forcing BESSI and ITM with transient climate forcing obtained from iLOVECLIM for both ice sheets. For present-day climate conditions, both BESSI and ITM exhibit good performance compared to MAR despite a much simpler model setup. While BESSI performs well for both Antarctica and Greenland for the same set of parame- ters, the ITM parameters need to be adapted specifically for each ice sheet. This suggests that the physics embedded in BESSI allows better capture of SMB changes across varying climate conditions, while ITM displays a much stronger sen- sitivity to its tunable parameters. The findings suggest that BESSI can provide more reliable SMB estimations for the iLOVECLIM framework to improve the model simulations of the ice sheet evolution and interactions with climate for multi-millennial simulations.

Water stable isotope records in polar ice cores have been largely used to reconstruct past local temperatures and other climatic information such as evaporative source region conditions of the precipitation reaching the ice core sites. However, recent studies have identified post-depositional processes taking place at the ice sheet's surface, modifying the original precipitation signal and challenging the traditional interpretation of ice core isotopic records. In this study, we use a combination of existing and new datasets of precipitation, snow surface, and subsurface isotopic compositions (δ18O and deuterium excess (d-excess)); meteorological parameters; ERA5 reanalyses; outputs from the isotope-enabled climate model ECHAM6-wiso; and a simple modelling approach to investigate the transfer function of water stable isotopes from precipitation to the snow surface and subsurface at Dome C in East Antarctica. We first show that water vapour fluxes at the surface of the ice sheet result in a net annual sublimation of snow, from 3.1 to 3.7 mm w.e. yr−1 (water equivalent) between 2018 and 2020, corresponding to 12 % to 15 % of the annual surface mass balance. We find that the precipitation isotopic signal cannot fully explain the mean, nor the variability in the isotopic composition observed in the snow, from annual to intra-monthly timescales. We observe that the mean effect of post-depositional processes over the study period enriches the snow surface in δ18O by 3.0 ‰ to 3.3 ‰ and lowers the snow surface d-excess by 3.4 ‰ to 3.5 ‰ compared to the incoming precipitation isotopic signal. We also show that the mean isotopic composition of the snow subsurface is not statistically different from that of the snow surface, indicating the preservation of the mean isotopic composition of the snow surface in the top centimetres of the snowpack. This study confirms previous findings about the complex interpretation of the water stable isotopic signal in the snow and provides the first quantitative estimation of the impact of post-depositional processes on the snow isotopic composition at Dome C, a crucial step for the accurate interpretation of isotopic records from ice cores.

Future climate and sea level projections depend sensitively on the response of the Antarctic Ice Sheet to ocean-driven melting and the resulting freshwater fluxes into the Southern Ocean. Circumpolar Deep Water (CDW) transport across the Antarctic continental shelf and into cavities beneath ice shelves is increasingly recognised as a crucial heat source for ice shelf melt. Quantifying past changes in the temperature of CDW is therefore of great benefit for modelling ice sheet response to past warm climates, for validating paleoclimate models, and for putting recent and projected changes in CDW temperature into context. Here we compile the available bottom water temperature reconstructions representative of CDW over the past 800 kyr. Estimated interglacial warming reached anomalies of +0.6 +/- 0.4 degrees C (MIS 11) and +0.5 +/- 0.5 degrees C (MIS 5) relative to present. Glacial cooling typically reached anomalies of ca. -1.5 to -2 degrees C, therefore maintaining positive thermal forcing for ice shelf melt even during glacials in the Amundsen Sea region of West Antarctica. Despite high variance amongst a small number of records and poor (4 kyr) temporal resolution, we find persistent and close relationships between our estimated CDW temperature and Southern Ocean sea surface temperature, Antarctic surface air temperature, and global deep-water temperature reconstructions at glacial-cycle timescales. Given the important role that CDW plays in connecting the world's three main ocean basins and in driving Antarctic Ice Sheet mass loss, additional temperature reconstructions targeting CDW are urgently needed to increase temporal and spatial resolution and to decrease uncertainty in past CDW temperatures - whether for use as a boundary condition, for model validation, or for understanding past oceanographic changes.

Earth’s obliquity and eccentricity cycles are strongly imprinted on Earth’s climate and widely used to measure geological time. However, the record of these imprints on the oxygen isotope record in deep-sea benthic foraminifera (δ18Ob) shows contradictory signals that violate isotopic principles and cause controversy over climate-ice sheet interactions. Here, we present a δ18Ob record of high fidelity from International Ocean Drilling Program (IODP) Site U1406 in the northwest Atlantic Ocean. We compare our record to other records for the time interval between 28 and 20 million years ago, when Earth was warmer than today, and only Antarctic ice sheets existed. The imprint of eccentricity on δ18Ob is remarkably consistent globally whereas the obliquity signal is inconsistent between sites, indicating that eccentricity was the primary pacemaker of land ice volume. The larger eccentricity-paced early Antarctic ice ages were vulnerable to rapid termination. These findings imply that the self-stabilizing hysteresis effects of large land-based early Antarctic ice sheets were strong enough to maintain ice growth despite consecutive insolation-induced polar warming episodes. However, rapid ice age terminations indicate that resistance to melting was weaker than simulated by numerical models and regularly overpowered, sometimes abruptly.

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.

Ice streams regulate most ice mass loss in Antarctica. Determining ice stream response to warmer conditions during the Pliocene could provide insights into their future behaviour, but this is hindered by a poor representation of subglacial topography in ice-sheet models. We address this limitation using a high-resolution model for Dronning Maud Land (East Antarctica). We show that contrary to dynamic thinning of the region’s ice streams following ice-shelf collapse, the largest ice stream, Jutulstraumen, thickens by 700 m despite lying on a retrograde bed slope. We attribute this counterintuitive thickening to a shallower Pliocene subglacial topography and inherent high lateral stresses at its flux gate. These conditions constrict ice drainage and, combined with increased snowfall, allow ice accumulation upstream. Similar stress balances and increased precipitation projections occur across 27% of present-day East Antarctica, and understanding how lateral stresses regulate ice-stream discharge is necessary for accurately assessing Antarctica’s future sea-level rise contribution.

New Zealand was among the last habitable places on earth to be colonized by humans. Charcoal records indicate that wildfires were rare prior to colonization and widespread following the 13th- to 14th-century Māori settlement, but the precise timing and magnitude of associated biomass-burning emissions are unknown, as are effects on light-absorbing black carbon aerosol concentrations over the pristine Southern Ocean and Antarctica. Here we used an array of well-dated Antarctic ice-core records to show that while black carbon deposition rates were stable over continental Antarctica during the past two millennia, they were approximately threefold higher over the northern Antarctic Peninsula during the past 700 years. Aerosol modelling demonstrates that the observed deposition could result only from increased emissions poleward of 40° S—implicating fires in Tasmania, New Zealand and Patagonia—but only New Zealand palaeofire records indicate coincident increases. Rapid deposition increases started in 1297 (±30 s.d.) in the northern Antarctic Peninsula, consistent with the late 13th-century Māori settlement and New Zealand black carbon emissions of 36 (±21 2 s.d.) Gg y−1 during peak deposition in the 16th century. While charcoal and pollen records suggest earlier, climate-modulated burning in Tasmania and southern Patagonia, deposition in Antarctica shows that black carbon emissions from burning in New Zealand dwarfed other preindustrial emissions in these regions during the past 2,000 years, providing clear evidence of large-scale environmental effects associated with early human activities across the remote Southern Hemisphere.

The dominant feature of large-scale mass transfer in the modern ocean is the Atlantic meridional overturning circulation (AMOC). The geometry and vigour of this circulation influences global climate on various timescales. Palaeoceanographic evidence suggests that during glacial periods of the past 1.5 million years the AMOC had markedly different features from today1; in the Atlantic basin, deep waters of Southern Ocean origin increased in volume while above them the core of the North Atlantic Deep Water (NADW) shoaled2. An absence of evidence on the origin of this phenomenon means that the sequence of events leading to global glacial conditions remains unclear. Here we present multi-proxy evidence showing that northward shifts in Antarctic iceberg melt in the Indian–Atlantic Southern Ocean (0–50° E) systematically preceded deep-water mass reorganizations by one to two thousand years during Pleistocene-era glaciations. With the aid of iceberg-trajectory model experiments, we demonstrate that such a shift in iceberg trajectories during glacial periods can result in a considerable redistribution of freshwater in the Southern Ocean. We suggest that this, in concert with increased sea-ice cover, enabled positive buoyancy anomalies to ‘escape’ into the upper limb of the AMOC, providing a teleconnection between surface Southern Ocean conditions and the formation of NADW. The magnitude and pacing of this mechanism evolved substantially across the mid-Pleistocene transition, and the coeval increase in magnitude of the ‘southern escape’ and deep circulation perturbations implicate this mechanism as a key feedback in the transition to the ‘100-kyr world’, in which glacial–interglacial cycles occur at roughly 100,000-year periods.

Understanding climate proxy records that preserve physical characteristics of past climate is a prerequisite to reconstruct long-term climatic conditions. Water stable isotope ratios (δ18O) constitute a widely used proxy in ice cores to reconstruct temperature and climate. However, the original climate signal is altered between the formation of precipitation and the ice, especially in low-accumulation areas such as the East Antarctic Plateau. Atmospheric conditions under which the isotopic signal is acquired at Aurora Basin North (ABN), East Antarctica, are characterized with the regional atmospheric model Modèle Atmosphérique Régional (MAR). The model shows that 50% of the snow is accumulated in less than 24 days year−1. Snowfall occurs throughout the year and intensifies during winter, with 64% of total accumulation between April and September, leading to a cold bias of −0.86°C in temperatures above inversion compared to the annual mean of −29.7°C. Large snowfall events are associated with high-pressure systems forcing warm oceanic air masses toward the Antarctic interior, which causes a warm bias of +2.83°C. The temperature-δ18O relationship, assessed with the global atmospheric model ECHAM5-wiso, is primarily constrained by the winter variability, but the observed slope is valid year-round. Three snow δ18O records covering 2004–2014 indicate that the anomalies recorded in the ice core are attributable to the occurrence of warm winter storms bringing precipitation to ABN and support the interpretation of δ18O in this region as a marker of temperature changes related to large-scale atmospheric conditions, particularly blocking events and variations in the Southern Annular Mode.

Atmospheric CO2 concentrations (pCO2) varied on millennial timescales in phase with Antarctic temperature during the last glacial period. A prevailing view has been that carbon release and uptake by the Southern Ocean dominated this millennial-scale variability in pCO2. Here, using Earth System Model experiments with an improved parameterization of ocean vertical mixing, we find a major role for terrestrial and oceanic carbon releases in driving the pCO2 trend. In our simulations, a change in Northern Hemisphere insolation weakens the Atlantic Meridional Overturning Circulation (AMOC) leading to increasing pCO2 and Antarctic temperatures. The simulated rise in pCO2 is caused in equal parts by increased CO2 outgassing from the global ocean due to a reduced biological activity and changed ventilation rates, and terrestrial carbon release as a response to southward migration of the Intertropical Convergence Zone. The simulated terrestrial release of carbon could explain stadial declines in organic carbon reservoirs observed in recent ice core δ13C measurements. Our results show that parallel variations in Antarctic temperature and pCO2 do not necessitate that the Southern Ocean dominates carbon exchange; instead, changes in carbon flux from the global ocean and land carbon reservoirs can explain the observed pCO2 (and δ13C) changes.

In central Antarctica, where accumulation rates are very low, summer sublimation of surface snow is a key element of the surface mass balance, but its fingerprint in isotopic composition of water (δD, δ18O, and δ17O) remains unclear. In this study, we examined the influence of summer sublimation on δD, δ18O, and δ17O in precipitation using data sets of isotopic composition of precipitation at various sites on the inland East Antarctica. We found unexpectedly low δ18O values in the summer precipitation, decoupled from surface air temperatures. This feature can be explained by the combined effects of weak or nonexistent temperature inversion and moisture recycling associated with sublimation-condensation processes in summer. Isotopic fractionation during the moisture-recycling process also explains the observed high values of d-excess and 17O-excess in summer precipitation. Our results suggest that the local cycle of sublimation-condensation in summer is an important process for the isotopic composition of surface snow, water vapor, and consequently precipitation on inland East Antarctica.

The glacimarine environment of the Antarctic Peninsula region is one of the fastest warming places on Earth today, but details of changes in the recent past remain unknown. Large distances and widespread variability separate late Holocene palaeoclimate reconstructions in this region. This study focuses on a marine sediment core collected from ca. 2000 m below sea level in the Central Bransfield Strait that serves as a key for understanding changes in this region. The core yielded a high sedimentation rate and therefore provides an exceptional high-resolution sedimentary record composed of hemipelagic sediment, with some turbidites. An age model has been created using radiocarbon dates that span the Late Holocene: 3560 cal yr BP to present. This chronostratigraphic framework was used to establish five units, which are grouped into two super-units: a lower super-unit (3560–1600 cal yr BP) and an upper super-unit (1600 cal yr BP–present), based on facies descriptions, laser particle size analysis, x-ray analysis, multi-sensor core logger data, weight percentages and isotopic values of total organic carbon and nitrogen. We interpret the signal contained within the upper super-unit as an increase in surface water irradiance and/or shortening of the sea-ice season and the five units are broadly synchronous with climatic intervals across the Antarctic Peninsula region. While the general trends of regional climatic periods are represented in the Bransfield Basin core we have examined, each additional record that is obtained adds variability to the known history of the Antarctic Peninsula, rather than clarifying specific trends. Keywords: Antarctic Peninsula; palaeoclimate; Holocene; marine; isotopes.

The Integrated Ocean Drilling Program Expedition 318 to the Wilkes Land margin of Antarctica recovered a sedimentary succession ranging in age from lower Eocene to the Holocene. Excellent stratigraphic control is key to understanding the timing of paleoceanographic events through critical climate intervals. Drill sites recovered the lower and middle Eocene, nearly the entire Oligocene, the Miocene from about 17 Ma, the entire Pliocene and much of the Pleistocene. The paleomagnetic properties are generally suitable for magnetostratigraphic interpretation, with well-behaved demagnetization diagrams, uniform distribution of declinations, and a clear separation into two inclination modes. Although the sequences were discontinuously recovered with many gaps due to coring, and there are hiatuses from sedimentary and tectonic processes, the magnetostratigraphic patterns are in general readily interpretable. Our interpretations are integrated with the diatom, radiolarian, calcareous nannofossils and dinoflagellate cyst (dinocyst) biostratigraphy. The magnetostratigraphy significantly improves the resolution of the chronostratigraphy, particularly in intervals with poor biostratigraphic control. However, Southern Ocean records with reliable magnetostratigraphies are notably scarce, and the data reported here provide an opportunity for improved calibration of the biostratigraphic records. In particular, we provide a rare magnetostratigraphic calibration for dinocyst biostratigraphy in the Paleogene and a substantially improved diatom calibration for the Pliocene. This paper presents the stratigraphic framework for future paleoceanographic proxy records which are being developed for the Wilkes Land margin cores. It further provides tight constraints on the duration of regional hiatuses inferred from seismic surveys of the region.

The chemical properties of the mid-depth and deep Southern Ocean are diagnostic of the mechanisms of abrupt changes in the global ocean throughout the late Pleistocene, because the regional water mass conversion and mixing help determine global ocean gradients. Here we define continuous time series of Southern Ocean vertical gradients by differencing the records from two high deposition rate deep sea sedimentary sequences that span the last several ice age cycles. The inferred changes in vertical carbon and oxygen isotopic gradients were dominated by variability on the millennial scale, and they followed closely the abrupt climate events of the high latitude Northern Hemisphere. In particular, the stadial events of at least the last 200kyr were characterized by enhanced mid-deep gradients in both δ13C (dissolved inorganic carbon) and δ18O (temperature). Interstadial events, conversely, featured reduced vertical gradients in both properties. The glacial terminations represented exceptions to this pattern of variability, as the vertical carbon isotopic gradient flattened dramatically at times of peak warmth in the Southern Ocean surface waters and with little or no corresponding change δ18O gradient. The available evidence suggests that properties of the upper layer of the Southern Ocean (Antarctic Intermediate Water) were influenced by an atmospherically mediated teleconnection to high latitude Northern Hemisphere.

Holocene climate variability in the southeast Atlantic sector of the Southern Ocean and Antarctic is assessed and quantified through integration of available marine sediment core and Antarctic ice core data. We use summer sea surface temperature (SSST) and sea ice presence (SIP) reconstructions from two marine sediment cores recovered north (50 °S) and south (53.2 °S) of the present day Antarctic Polar Front (APF), as well as an atmospheric temperature and sea ice proxy from the EPICA ice core from Dronning Maud Land (EDML). We find reasonably good agreement in the timing of climate evolution in the analyzed series. Almost all records show a gradual glacial-to-Holocene climate transition, interrupted by the Antarctic cold reversal around 13 000 cal yr BP, and early Holocene climatic optimum (HCO) at about 11 000 cal yr BP. During the early HCO, the seasonal ice cover retreats to south of 53 °S; it then readvances in the course of the mid- to late Holocene. The maximum winter sea ice edge position during the recent 10 000 years varied mainly within 51–53 °S, with sporadic growth to north of 50 °S, a position similar to that during the last glacial. The onset of the Neoglacial period after ca 4000 yr BP is associated with a steepening of the SSST gradient between the marine core sites, strengthening of the westerlies and cooling in the inland ice sheet. The agreement in timing between elevated SSST during the early HCO and decreased deuterium excess in EDML and other ice cores from different locations in the East Antarctic suggests that the retreat of sea ice during the early HCO and weakening of the APF was a general feature of the East Antarctic climate during that time.

A quantitative biochronological study by Cody et al. (2008) integrates comprehensive diatom biostratigraphy, magnetostratigraphy, and tephrostratigraphy from 32 Neogene sections around the Southern Ocean and Antarctic continental margin. A recent method, known as Constrained Optimization (CONOP), which can be viewed as a multidimensional version of graphic correlation, is applied to that very interesting database. The goal of the present paper is to discuss some theoretical aspects of quantitative biochronology and to compare the constrained optimization with the deterministic method called Unitary Associations (UAM), a graph theoretical model. We illustrate the fact that the UAM is an extremely powerful and unique theory allowing an in-depth analysis of the internal conflicting inter-taxon stratigraphic relationships, inherent to any complex biostratigraphical database.

The geographic extent of cooling associated with the Antarctic Cold Reversal is unclear. Dating of glacial moraines in New Zealand suggests that the cooling extended into the southern mid-latitudes, possibly as a result of the northward migration of the southern subtropical front.

We investigated deep water changes in the Southern Ocean during the last glacial inception, in relationship to surface hydrology and global climatology, to better understand the mechanisms of the establishment of a glacial ocean circulation. Changes in benthic foraminiferal δ13C from three high-resolution cores are compared and indicate decoupled intermediate and deep water changes in the Southern Ocean. From the comparison with records from the North Atlantic, South Atlantic, and the Southern Ocean, we show that the early southern deep water δ13C drop observed at the MIS 5.5–5.4 transition occurred before any significant reduction of North Atlantic Deep Water ventilation. We propose that this drop is linked to the northward expansion of poorly ventilated Antarctic Bottom Water (AABW) mass in the Southern Ocean. Associated with an early cooling in the high southern latitudes, the westerly winds and surface oceanic fronts would migrate equatorward, thus weakening the upwelling of Circumpolar Deep Waters. Reduced heat brought to Antarctic surface waters would enhance sea ice formation during winters and the deep convection of cold and poorly ventilated AABW.

The East Antarctic Ice Sheet is the largest, highest, coldest, driest, and windiest ice sheet on Earth. Understanding of the surface mass balance (SMB) of Antarctica is necessary to determine the present state of the ice sheet, to make predictions of its potential contribution to sea level rise, and to determine its past history for paleoclimatic reconstructions. However, SMB values are poorly known because of logistic constraints in extreme polar environments, and they represent one of the biggest challenges of Antarctic science. Snow accumulation is the most important parameter for the SMB of ice sheets. SMB varies on a number of scales, from small-scale features (sastrugi) to ice-sheet-scale SMB patterns determined mainly by temperature, elevation, distance from the coast, and wind-driven processes. In situ measurements of SMB are performed at single points by stakes, ultrasonic sounders, snow pits, and firn and ice cores and laterally by continuous measurements using ground-penetrating radar. SMB for large regions can only be achieved practically by using remote sensing and/or numerical climate modeling. However, these techniques rely on ground truthing to improve the resolution and accuracy. The separation of spatial and temporal variations of SMB in transient regimes is necessary for accurate interpretation of ice core records. In this review we provide an overview of the various measurement techniques, related difficulties, and limitations of data interpretation; describe spatial characteristics of East Antarctic SMB and issues related to the spatial and temporal representativity of measurements; and provide recommendations on how to perform in situ measurements.

The termination of the last ice age (Termination 1; T1) is crucial for our understanding of global climate change and for the validation of climate models. There are still a number of open questions regarding for example the exact timing and the mechanisms involved in the initiation of deglaciation and the subsequent interhemispheric pattern of the warming. Our study is based on a well-dated and high-resolution alkenone-based sea surface temperature (SST) record from the SE-Pacific off southern Chile (Ocean Drilling Project Site 1233) showing that deglacial warming at the northern margin of the Antarctic Circumpolar Current system (ACC) began shortly after 19,000 years BP (19 kyr BP). The timing is largely consistent with Antarctic ice-core records but the initial warming in the SE-Pacific is more abrupt suggesting a direct and immediate response to the slowdown of the Atlantic thermohaline circulation through the bipolar seesaw mechanism. This response requires a rapid transfer of the Atlantic signal to the SE-Pacific without involving the thermal inertia of the Southern Ocean that may contribute to the substantially more gradual deglacial temperature rise seen in Antarctic ice-cores. A very plausible mechanism for this rapid transfer is a seesaw-induced change of the coupled ocean–atmosphere system of the ACC and the southern westerly wind belt. In addition, modelling results suggest that insolation changes and the deglacial CO2 rise induced a substantial SST increase at our site location but with a gradual warming structure. The similarity of the two-step rise in our proxy SSTs and CO2 over T1 strongly demands for a forcing mechanism influencing both, temperature and CO2. As SSTs at our coring site are particularly sensitive to latitudinal shifts of the ACC/southern westerly wind belt system, we conclude that such latitudinal shifts may substantially affect the upwelling of deepwater masses in the Southern Ocean and thus the release of CO2 to the atmosphere as suggested by the conceptual model of [Toggweiler, J.R., Rusell, J.L., Carson, S.R., 2006. Midlatitude westerlies, atmospheric CO2, and climate change during ice ages. Paleoceanography 21. doi:10.1029/2005PA001154].