Climate at high-obliquity

The question of climate at high obliquity is raised in the context of both exoplanet studies (e.g. habitability) and paleoclimates studies (evidence for low-latitude glaciation during the Neoproterozoic and the ”Snowball Earth” hypothesis). States of high obliquity, φ, are distinctive in that, for φ ≥54◦, the poles receive more solar radiation in the annual mean than the Equator, opposite to the present day situation. In addition, the seasonal cycle of insolation is extreme, with the poles alternatively “facing” the sun and sheltering in the dark for months. The novelty of our approach is to consider the role of a dynamical ocean in controlling the surface climate at high obliquity, which in turn requires understanding of the surface winds patterns when temperature gradients are reversed. To address these questions, a coupled ocean-atmosphere-sea ice GCM configured on an aquaplanet is employed. Except for the absence of topography and modified obliquity, the set-up is Earth-like. Two large obliquities φ, 54◦ and 90◦, are compared to today’s Earth value, φ=23.5◦. Three key results emerge at high obliquity: 1) despite reversed temper- ature gradients, mid-latitudes surface winds are westerly and trade winds exist at the equator (as for φ=23.5◦) although the westerlies are confined to the summer hemisphere, 2) a habitable planet is possible with mid-latitude temperatures in the range 300-280 K and 3) a stable climate state with an ice cap limited to the equatorial region is unlikely. We clarify the dynamics behind these features (notably by an analysis of the potential vorticity structure and conditions for baroclinic instability of the atmosphere). Interestingly, we find that the absence of a stable partially glaciated state is critically linked to the absence of ocean heat transport during winter, a feature ultimately traced back to the high seasonality of baroclinic instability conditions in the atmosphere.

Mean climate and variability of the atmosphere and ocean on an aquaplanet

Numerical experiments are described that pertain to the climate of a coupled atmosphere–ocean–ice system in the absence of land, driven by modern-day orbital and CO2 forcing. Millennial time-scale simulations yield a mean state in which ice caps reach down to 55° of latitude and both the atmosphere and ocean comprise eastward- and westward-flowing zonal jets, whose structure is set by their respective baroclinic instabilities. Despite the zonality of the ocean, it is remarkably efficient at transporting heat meridionally through the agency of Ekman transport and eddy-driven subduction. Indeed the partition of heat transport between the atmosphere and ocean is much the same as the present climate, with the ocean dominating in the Tropics and the atmosphere in the mid–high latitudes. Variability of the system is dominated by the coupling of annular modes in the atmosphere and ocean. Stochastic variability inherent to the atmospheric jets drives variability in the ocean. Zonal flows in the ocean exhibit decadal variability, which, remarkably, feeds back to the atmosphere, coloring the spectrum of annular variability. A simple stochastic model can capture the essence of the process. Finally, it is briefly reviewed how the aquaplanet can provide information about the processes that set the partition of heat transport and the climate of Earth.

An event-based approach to understanding decadal fluctuations in the Atlantic meridional overturning circulation

Many previous studies have shown that unforced climate model simulations exhibit decadal-scale fluctuations in the Atlantic meridional overturning circulation (AMOC), and that this variability can have impacts on surface climate fields. However, the robustness of these surface fingerprints across different models is less clear. Furthermore, with the potential for coupled feedbacks that may amplify or damp the response, it is not known whether the associated climate signals are linearly related to the strength of the AMOC changes, or if the fluctuation events exhibit nonlinear behaviour with respect to their strength or polarity. To explore these questions, we introduce an objective and flexible method for identifying the largest natural AMOC fluctuation events in multicentennial/multimillennial simulations of a variety of coupled climate models. The characteristics of the events are explored, including their magnitude, meridional coherence and spatial structure, as well as links with ocean heat transport and the horizontal circulation. The surface fingerprints in ocean temperature and salinity are examined, and compared with the results of linear regression analysis. It is found that the regressions generally provide a good indication of the surface changes associated with the largest AMOC events. However, there are some exceptions, including a nonlinear change in the atmospheric pressure signal, particularly at high latitudes, in HadCM3. Some asymmetries are also found between the changes associated with positive and negative AMOC events in the same model. Composite analysis suggests that there are signals that are robust across the largest AMOC events in each model, which provides reassurance that the surface changes associated with one particular event will be similar to those expected from regression analysis. However, large differences are found between the AMOC fingerprints in different models, which may hinder the prediction and attribution of such events in reality.

Atmospheric and oceanic contributions to irreducible forecast uncertainty of Arctic surface climate

Uncertainty of Arctic seasonal to interannual predictions arising from model errors and initial state uncertainty has been widely discussed in the literature, whereas the irreducible forecast uncertainty (IFU) arising from the chaoticity of the climate system has received less attention. However, IFU provides important insights into the mechanisms through which predictability is lost, and hence can inform prioritization of model development and observations deployment. Here, we characterize how internal oceanic and surface atmospheric heat fluxes contribute to IFU of Arctic sea ice and upper ocean heat content in an Earth system model by analyzing a set of idealized ensemble prediction experiments. We find that atmospheric and oceanic heat flux are often equally important for driving unpredictable Arctic-wide changes in sea ice and surface water temperatures, and hence contribute equally to IFU. Atmospheric surface heat flux tends to dominate Arctic-wide changes for lead times of up to a year, whereas oceanic heat flux tends to dominate regionally and on interannual time scales. There is in general a strong negative covariance between surface heat flux and ocean vertical heat flux at depth, and anomalies of lateral ocean heat transport are wind-driven, which suggests that the unpredictable oceanic heat flux variability is mainly forced by the atmosphere. These results are qualitatively robust across different initial states, but substantial variations in the amplitude of IFU exist. We conclude that both atmospheric variability and the initial state of the upper ocean are key ingredients for predictions of Arctic surface climate on seasonal to interannual time scales.

A reversal of climatic trends in the North Atlantic since 2005

In the mid-1990s the North Atlantic subpolar gyre warmed rapidly, which had important climate impacts, such as increased hurricane numbers, and changes to rainfall over Africa, Europe and North America. Evidence suggests that the warming was largely due to a strengthening of the ocean circulation, particularly the Atlantic Meridional Overturning Circulation (AMOC). Since the mid-1990s direct and indirect measurements have suggested a decline in the strength of the ocean circulation, which is expected to lead to a reduction in northward heat transport. Here we show that since 2005 a large volume of the upper North Atlantic Ocean has cooled significantly by approximately -0.45C or 1.5x10^22 J, reversing the previous warming trend. By analysing observations and a state-of-the-art climate model, we show that this cooling is consistent with a reduction in the strength of the ocean circulation and heat transport, linked to record low densities in the deep Labrador Sea. The low density in the deep Labrador Sea is primarily due to deep ocean warming since 1995, but a long-term freshening also played a role. The observed upper ocean cooling since 2005 is not consistent with the hypothesis that anthropogenic aerosols directly drive Atlantic temperatures.

Temperature Dynamics Investigation At Small And Shallow Lakes Using Hydrodynamic Model

A three-dimensional time-dependent hydrodynamic and heat transport model of Lake Binaba, a shallow and small dam reservoir in Ghana, emphasizing the simulation of dynamics and thermal structure has been developed. Most numerical studies of temperature dynamics in reservoirs are based on one- or two-dimensional models. These models are not applicable for reservoirs characterized with complex flow pattern and unsteady heat exchange between the atmosphere and water surface. Continuity, momentum and temperature transport equations have been solved. Proper assignment of boundary conditions, especially surface heat fluxes, has been found crucial in simulating the lake’s hydrothermal dynamics. This model is based on the Reynolds Average Navier-Stokes equations, using a Boussinesq approach, with a standard k − ε turbulence closure to solve the flow field. The thermal model includes a heat source term, which takes into account the short wave radiation and also heat convection at the free surface, which is function of air temperatures, wind velocity and stability conditions of atmospheric boundary layer over the water surface. The governing equations of the model have been solved by OpenFOAM; an open source, freely available CFD toolbox. As its core, OpenFOAM has a set of efficient C++ modules that are used to build solvers. It uses collocated, polyhedral numerics that can be applied on unstructured meshes and can be easily extended to run in parallel. A new solver has been developed to solve the hydrothermal model of lake. The simulated temperature was compared against a 15 days field data set. Simulated and measured temperature profiles in the probe locations show reasonable agreement. The model might be able to compute total heat storage of water bodies to estimate evaporation from water surface.

State of the Arctic October 2006

This State of the Arctic Report presents a review of recent data by an international group of scientists who developed a consensus on the information content and reliability. The report highlights data primarily from 2000 to 2005 with a first look at winter 2006, providing an update to some of the records of physical processes discussed in the Arctic Climate Impact Assessment (ACIA, 2004, 2005). Of particular note: • Atmospheric climate patterns are shifting (Fig. 1). The late winter/spring pattern for 2000–2005 had new hot spots in northeast Canada and the East Siberian Sea relative to 1980–1999. Late winter 2006, however, shows a return to earlier climate patterns, with warm temperatures in the extended region near Svalbard. • Ocean salinity and temperature profiles at the North Pole and in the Beaufort Sea, which changed abruptly in the 1990s, show that conditions since 2000 have relaxed toward the pre-1990 climatology, although 2001–2004 has seen an increase in northward ocean heat transport through Bering Strait (Fig. 2), which is thought to impact sea ice loss. • Sea ice extent continues to decrease. The sea ice extent in September 2005 was the minimum observed in summer during the satellite era (beginning in 1979), marking an unprecedented series of extreme ice extent minima beginning in 2002 (Fig. 3). The sea ice extent in March 2006 was also the minimum observed in winter during the satellite era. • Tundra vegetation greenness increased, primarily due to an increase in the abundance of shrubs. Boreal forest vegetation greenness decreased, possibly due to drought conditions (Fig. 4). • There is increasing interest in the stability of the Greenland ice sheet. The velocity of outlet glaciers increased in 2005 relative to 2000 and 1995, but uncertainty remains with regard to the total mass balance. • Permafrost temperatures continue to increase. However, data on changes in the active layer thickness (the relatively thin layer of ground between the surface and permafrost that undergoes seasonal freezing and thawing) are less conclusive. While some of the sites show a barely noticeable increasing trend in the thickness of the active layer, most of them do not. • Globally, 2005 was the warmest year in the instrumental record (beginning in 1880), with the Arctic providing a large contribution toward this increase. Many of the trends documented in the ACIA are continuing, but some are not. Taken collectively, the observations presented in this report indicate that during 2000–2005 the Arctic system showed signs of continued warming. However, there are a few indications that certain elements may be recovering and returning to recent climatological norms (for example, the central Arctic Ocean and some wind patterns). These mixed tendencies further illustrate the sensitivity and complexity of the Arctic physical system. They underline the importance of maintaining and expanding efforts to observe and better understand this important component of the climate system to provide accurate predictions of its future state.

Wind-driven cross-equatorial flow in the Indian ocean

Máster en Oceanografía

Structuring of polymer surface by evaporation of sessile microdrops

In this thesis, we investigated the evaporation of sessile microdroplets on different solid substrates. Three major aspects were studied: the influence of surface hydrophilicity and heterogeneity on the evaporation dynamics for an insoluble solid substrate, the influence of external process parameters and intrinsic material properties on microstructuring of soluble polymer substrates and the influence of an increased area to volume ratio in a microfluidic capillary, when evaporation is hindered. In the first part, the evaporation dynamics of pure sessile water drops on smooth self-assembled monolayers (SAMs) of thiols or disulfides on gold on mica was studied. With increasing surface hydrophilicity the drop stayed pinned longer. Thus, the total evaporation time of a given initial drop volume was shorter, since the drop surface, through which the evaporation occurs, stays longer large. Usually, for a single drop the volume decreased linearly with t1.5, t being the evaporation time, for a diffusion-controlled evaporation process. However, when we measured the total evaporation time, ttot, for multiple droplets with different initial volumes, V0, we found a scaling of the form V0 = attotb. The more hydrophilic the substrate was, the more showed the scaling exponent a tendency to an increased value up to 1.6. This can be attributed to an increasing evaporation rate through a thin water layer in the vicinity of the drop. Under the assumption of a constant temperature at the substrate surface a cooling of the droplet and thus a decreased evaporation rate could be excluded as a reason for the different scaling exponent by simulations performed by F. Schönfeld at the IMM, Mainz. In contrast, for a hairy surface, made of dialkyldisulfide SAMs with different chain lengths and a 1:1 mixture of hydrophilic and hydrophobic end groups (hydroxy versus methyl group), the scaling exponent was found to be ~ 1.4. It increased to ~ 1.5 with increasing hydrophilicity. A reason for this observation can only be speculated: in the case of longer hydrophobic alkyl chains the formation of an air layer between substrate and surface might be favorable. Thus, the heat transport to the substrate might be reduced, leading to a stronger cooling and thus decreased evaporation rate. In the second part, the microstructuring of polystyrene surfaces by drops of toluene, a good solvent, was investigated. For this a novel deposition technique was developed, with which the drop can be deposited with a syringe. The polymer substrate is lying on a motorized table, which picks up the pendant drop by an upward motion until a liquid bridge is formed. A consecutive downward motion of the table after a variable delay, i.e. the contact time between drop and polymer, leads to the deposition of the droplet, which can evaporate. The resulting microstructure is investigated in dependence of the processes parameters, i.e. the approach and the retraction speed of the substrate and the delay between them, and in dependence of the intrinsic material properties, i.e. the molar mass and the type of the polymer/solvent system. The principal equivalence with the microstructuring by the ink-jet technique was demonstrated. For a high approach and retraction speed of 9 mm/s and no delay between them, a concave microtopology was observed. In agreement with the literature, this can be explained by a flow of solvent and the dissolved polymer to the rim of the pinned droplet, where polymer is accumulated. This effect is analogue to the well-known formation of ring-like stains after the evaporation of coffee drops (coffee-stain effect). With decreasing retraction speed down to 10 µm/s the resulting surface topology changes from concave to convex. This can be explained with the increasing dissolution of polymer into the solvent drop prior to the evaporation. If the polymer concentration is high enough, gelation occurs instead of a flow to the rim and the shape of the convex droplet is received. With increasing delay time from below 0 ms to 1s the depth of the concave microwells decreases from 4.6 µm to 3.2 µm. However, a convex surface topology could not be obtained, since for longer delay times the polymer sticks to the tip of the syringe. Thus, by changing the delay time a fine-tuning of the concave structure is accomplished, while by changing the retraction speed a principal change of the microtopolgy can be achieved. We attribute this to an additional flow inside the liquid bridge, which enhanced polymer dissolution. Even if the pendant drop is evaporating about 30 µm above the polymer surface without any contact (non-contact mode), concave structures were observed. Rim heights as high as 33 µm could be generated for exposure times of 20 min. The concave structure exclusively lay above the flat polymer surface outside the structure even after drying. This shows that toluene is taken up permanently. The increasing rim height, rh, with increasing exposure time to the solvent vapor obeys a diffusion law of rh = rh0  tn, with n in the range of 0.46 ~ 0.65. This hints at a non-Fickian swelling process. A detailed analysis showed that the rim height of the concave structure is modulated, unlike for the drop deposition. This is due to the local stress relaxation, which was initiated by the increasing toluene concentration in the extruded polymer surface. By altering the intrinsic material parameters i.e. the polymer molar mass and the polymer/solvent combination, several types of microstructures could be formed. With increasing molar mass from 20.9 kDa to 1.44 MDa the resulting microstructure changed from convex, to a structure with a dimple in the center, to concave, to finally an irregular structure. This observation can be explained if one assumes that the microstructuring is dominated by two opposing effects, a decreasing solubility with increasing polymer molar mass, but an increasing surface tension gradient leading to instabilities of Marangoni-type. Thus, a polymer with a low molar mass close or below the entanglement limit is subject to a high dissolution rate, which leads to fast gelation compared to the evaporation rate. This way a coffee-rim like effect is eliminated early and a convex structure results. For high molar masses the low dissolution rate and the low polymer diffusion might lead to increased surface tension gradients and a typical local pile-up of polymer is found. For intermediate polymer masses around 200 kDa, the dissolution and evaporation rate are comparable and the typical concave microtopology is found. This interpretation was supported by a quantitative estimation of the diffusion coefficient and the evaporation rate. For a different polymer/solvent system, polyethylmethacrylate (PEMA)/ethylacetate (EA), exclusively concave structures were found. Following the statements above this can be interpreted with a lower dissolution rate. At low molar masses the concentration of PEMA in EA most likely never reaches the gelation point. Thus, a concave instead of a convex structure occurs. At the end of this section, the optically properties of such microstructures for a potential application as microlenses are studied with laser scanning confocal microscopy. In the third part, the droplet was confined into a glass microcapillary to avoid evaporation. Since here, due to an increased area to volume ratio, the surface properties of the liquid and the solid walls became important, the influence of the surface hydrophilicity of the wall on the interfacial tension between two immiscible liquid slugs was investigated. For this a novel method for measuring the interfacial tension between the two liquids within the capillary was developed. This technique was demonstrated by measuring the interfacial tensions between slugs of pure water and standard solvents. For toluene, n-hexane and chloroform 36.2, 50.9 and 34.2 mN/m were measured at 20°C, which is in a good agreement with data from the literature. For a slug of hexane in contact with a slug of pure water containing ethanol in a concentration range between 0 and 70 (v/v %), a difference of up to 6 mN/m was found, when compared to commercial ring tensiometry. This discrepancy is still under debate.

Impact of suborbital climate changes in the North Atlantic on ice sheet dynamics at the Mid-Pleistocene Transition

[1] Early and Mid-Pleistocene climate, ocean hydrography and ice sheet dynamics have been reconstructed using a high-resolution data set (planktonic and benthicδ18O time series, faunal-based sea surface temperature (SST) reconstructions and ice-rafted debris (IRD)) record from a high-deposition-rate sedimentary succession recovered at the Gardar Drift formation in the subpolar North Atlantic (Integrated Ocean Drilling Program Leg 306, Site U1314). Our sedimentary record spans from late in Marine Isotope Stage (MIS) 31 to MIS 19 (1069–779 ka). Different trends of the benthic and planktonic oxygen isotopes, SST and IRD records before and after MIS 25 (∼940 ka) evidence the large increase in Northern Hemisphere ice-volume, linked to the cyclicity change from the 41-kyr to the 100-kyr that occurred during the Mid-Pleistocene Transition (MPT). Beside longer glacial-interglacial (G-IG) variability, millennial-scale fluctuations were a pervasive feature across our study. Negative excursions in the benthicδ18O time series observed at the times of IRD events may be related to glacio-eustatic changes due to ice sheets retreats and/or to changes in deep hydrography. Time series analysis on surface water proxies (IRD, SST and planktonicδ18O) of the interval between MIS 31 to MIS 26 shows that the timing of these millennial-scale climate changes are related to half-precessional (10 kyr) components of the insolation forcing, which are interpreted as cross-equatorial heat transport toward high latitudes during both equinox insolation maxima at the equator.

Amplified inception of European Little Ice Age by sea ice-ocean-atmosphere feedbacks

The inception of the Little Ice Age (~1400–1700 AD) is believed to have been driven by an interplay of external forcing and climate system internal variability. While the hemispheric signal seems to have been dominated by solar irradiance and volcanic eruptions, the understanding of mechanisms shaping the climate on a continental scale is less robust. In an ensemble of transient model simulations and a new type of sensitivity experiments with artificial sea ice growth, the authors identify a sea ice–ocean–atmosphere feedback mechanism that amplifies the Little Ice Age cooling in the North Atlantic–European region and produces the temperature pattern suggested by paleoclimatic reconstructions. Initiated by increasing negative forcing, the Arctic sea ice substantially expands at the beginning of the Little Ice Age. The excess of sea ice is exported to the subpolar North Atlantic, where it melts, thereby weakening convection of the ocean. Consequently, northward ocean heat transport is reduced, reinforcing the expansion of the sea ice and the cooling of the Northern Hemisphere. In the Nordic Seas, sea surface height anomalies cause the oceanic recirculation to strengthen at the expense of the warm Barents Sea inflow, thereby further reinforcing sea ice growth. The absent ocean–atmosphere heat flux in the Barents Sea results in an amplified cooling over Northern Europe. The positive nature of this feedback mechanism enables sea ice to remain in an expanded state for decades up to a century, favoring sustained cold periods over Europe such as the Little Ice Age. Support for the feedback mechanism comes from recent proxy reconstructions around the Nordic Seas.

A Numerical Investigation of the Gulf Stream and Its Meanders in Response to Cold Air Outbreaks

The three-dimensional Princeton Ocean Model is used to examine the modification of the Gulf Stream and its meanders by cold air outbreaks. Two types of Gulf Stream meanders are found in the model. Meanders on the shoreward side of the Gulf Stream are baroclinically unstable. They are affected little by the atmospheric forcing because their energy source is stored at the permanent thermocline, well below the influence of the surface forcing. Meanders on the seaward side of the stream are both barotropically and baroclinically unstable. The energy feeding these meanders is stored at the surface front separating the Gulf Stream and the Sargasso Seal which is greatly reduced in case of cold air outbreaks. Thus, meanders there reduce strength and also seem to slow their downstream propagation due to the southward Ekman flow. Heat budget calculations suggest two almost separable processes. The oceanic heal released to the atmosphere during these severe cooling episodes comes almost exclusively from the upper water column. Transport of heat by meanders from the Gulf Stream to the shelf, though it is large, does not disrupt the principal balance. It is balanced nicely with the net heat transport in the downstream direction.

A multi-model comparison of Atlantic multidecadal variability

A multi-model analysis of Atlantic multidecadal variability is performed with the following aims: to investigate the similarities to observations; to assess the strength and relative importance of the different elements of the mechanism proposed by Delworth et al. (J Clim 6:1993–2011, 1993) (hereafter D93) among coupled general circulation models (CGCMs); and to relate model differences to mean systematic error. The analysis is performed with long control simulations from ten CGCMs, with lengths ranging between 500 and 3600 years. In most models the variations of sea surface temperature (SST) averaged over North Atlantic show considerable power on multidecadal time scales, but with different periodicity. The SST variations are largest in the mid-latitude region, consistent with the short instrumental record. Despite large differences in model configurations, we find quite some consistency among the models in terms of processes. In eight of the ten models the mid-latitude SST variations are significantly correlated with fluctuations in the Atlantic meridional overturning circulation (AMOC), suggesting a link to northward heat transport changes. Consistent with this link, the three models with the weakest AMOC have the largest cold SST bias in the North Atlantic. There is no linear relationship on decadal timescales between AMOC and North Atlantic Oscillation in the models. Analysis of the key elements of the D93 mechanisms revealed the following: Most models present strong evidence that high-latitude winter mixing precede AMOC changes. However, the regions of wintertime convection differ among models. In most models salinity-induced density anomalies in the convective region tend to lead AMOC, while temperature-induced density anomalies lead AMOC only in one model. However, analysis shows that salinity may play an overly important role in most models, because of cold temperature biases in their relevant convective regions. In most models subpolar gyre variations tend to lead AMOC changes, and this relation is strong in more than half of the models.

Mid-Pliocene palynomorph record of the North Atlantic

The mid-Pliocene was an episode of prolonged global warmth and strong North Atlantic thermohaline circulation, interrupted briefly at circa 3.30 Ma by a global cooling event corresponding to marine isotope stage (MIS) M2. Paleoceanographic changes in the eastern North Atlantic have been reconstructed between circa 3.35 and 3.24 Ma at Deep Sea Drilling Project Site 610 and Integrated Ocean Drilling Program Site 1308. Mg/Ca ratios and d18O from Globigerina bulloides are used to reconstruct the temperature and relative salinity of surface waters, and dinoflagellate cyst assemblages are used to assess variability in the North Atlantic Current (NAC). Our sea surface temperature data indicate warm waters at both sites before and after MIS M2 but a cooling of ~2-3°C during MIS M2. A dinoflagellate cyst assemblage overturn marked by a decline in Operculodinium centrocarpum reflects a southward shift or slowdown of the NAC between circa 3.330 and 3.283 Ma, reducing northward heat transport 23-35 ka before the global ice volume maximum of MIS M2. This will have established conditions that ultimately allowed the Greenland ice sheet to expand, leading to the global cooling event at MIS M2. Comparison with an ice-rafted debris record excludes fresh water input via icebergs in the northeast Atlantic as a cause of NAC decline. The mechanism causing the temporary disruption of the NAC may be related to a brief reopening of the Panamanian Gateway at about this time.

Dinoflagellate cyst assemblages, oxygen isotope records of foraminifera, and alkenone and Mg/Ca-based SST estimates for the time interval 3.40-3.18 Ma (Late Pliocene) from five North Atlantic and Caribbean sediment cores

The early Late Pliocene (3.6 to ~3.0 million years ago) is the last extended interval in Earth's history when atmospheric CO2 concentrations were comparable to today's and global climate was warmer. Yet a severe global glaciation during marine isotope stage (MIS) M2 interrupted this phase of global warmth ~3.30 million years ago, and is seen as a premature attempt of the climate system to establish an ice-age world. Our geochemical and palynological records from five marine sediment cores along a Caribbean to eastern North Atlantic transect show that increased Pacific-to-Atlantic flow via the Central American Seaway weakened the North Atlantic Current (NAC) and attendant northward heat transport prior to MIS M2. The consequent cooling of the northern high latitude oceans permitted expansion of the Greenland ice sheet during MIS M2, despite near-modern atmospheric CO2 concentrations. Before and after MIS M2, heat transport via the NAC was crucial in maintaining warm climates comparable to those predicted for the end of this century.