Proxy data and model simulations of sediment core GeoB7920-2

The evolution of the northwest African hydrological balance throughout the Pleistocene epoch influenced the migration of prehistoric humans**1. The hydrological balance is also thought to be important to global teleconnection mechanisms during Dansgaard-Oeschger and Heinrich events**2. However, most high-resolution African climate records do not span the millennial-scale climate changes of the last glacial-interglacial cycle**1, 3, 4, 5, or lack an accurate chronology**6. Here, we use grain-size analyses of siliciclastic marine sediments from off the coast of Mauritania to reconstruct changes in northwest African humidity over the past 120,000 years. We compare this reconstruction to simulations of palaeo-humidity from a coupled atmosphere-ocean-vegetation model. These records are in good agreement, and indicate the reoccurrence of precession-forced humid periods during the last interglacial period similar to the Holocene African Humid Period. We suggest that millennial-scale arid events are associated with a reduction of the North Atlantic meridional overturning circulation and that millennial-scale humid events are linked to a regional increase of winter rainfall over the coastal regions of northwest Africa.

Theoretical model simulations for variations in thermal reaction norms

Thermal reaction norms for growth rates of six Emiliania huxleyi isolates originating from the central Atlantic (Azores, Portugal) and five isolates from the coastal North Atlantic (Bergen, Norway) were assessed. We used the template mode of variation model to decompose variations in growth rates into modes of biological interest: vertical shift, horizontal shift, and generalist-specialist variation. In line with the actual habitat conditions, isolates from Bergen (Bergen population) grew well at lower temperatures, and isolates from the Azores (Azores population) performed better at higher temperatures. The optimum growth temperature of the Azores population was significantly higher than that of the Bergen population. Neutral genetic differentiation was found between populations by microsatellite analysis. These findings indicate that E. huxleyi populations are adapted to local temperature regimes. Next to between-population variation, we also found variation within populations. Genotype-by-environment interactions resulted in the most pronounced phenotypic differences when isolates were exposed to temperatures outside the range they naturally encounter. Variation in thermal reaction norms between and within populations emphasizes the importance of using more than one isolate when studying the consequences of global change on marine phytoplankton. Phenotypic plasticity and standing genetic variation will be important in determining the potential of natural E. huxleyi populations to cope with global climate change.

Model results of ocean physical field of the reference state and the targets in NetCDF format

In a feasibility study, the potential of proxy data for the temperature and salinity during the Last Glacial Maximum (LGM, about 19 000 to 23 000 years before present) in constraining the strength of the Atlantic meridional overturning circulation (AMOC) with a general ocean circulation model was explored. The proxy data were simulated by drawing data from four different model simulations at the ocean sediment core locations of the Multiproxy Approach for the Reconstruction of the Glacial Ocean surface (MARGO) project, and perturbing these data with realistic noise estimates. The results suggest that our method has the potential to provide estimates of the past strength of the AMOC even from sparse data, but in general, paleo-sea-surface temperature data without additional prior knowledge about the ocean state during the LGM is not adequate to constrain the model. On the one hand, additional data in the deep-ocean and salinity data are shown to be highly important in estimating the LGM circulation. On the other hand, increasing the amount of surface data alone does not appear to be enough for better estimates. Finally, better initial guesses to start the state estimation procedure would greatly improve the performance of the method. Indeed, with a sufficiently good first guess, just the sea-surface temperature data from the MARGO project promise to be sufficient for reliable estimates of the strength of the AMOC.

LPJ-LMfire simulations for Europe for the Last Glacial Maximum and preindustrial control, link to netCDF files

This data sets contains LPJ-LMfire dynamic global vegetation model output covering Europe and the Mediterranean for the Last Glacial Maximum (LGM; 21 ka) and for a preindustrial control simulation (20th century detrended climate). The netCDF data files are time averages of the final 30 years of the model simulation. Each netCDF file contains four or five variables: fractional cover of 9 plant functional types (PFTs; cover), total fractional coverage of trees (treecover), population density of hunter-gatherers (foragerPD; only for the "people" simulations), fraction of the gridcell burned on 30-year average (burnedf), and vegetation net primary productivity (NPP). The model spatial resolution is 0.5-degrees For the LGM simulations, LPJ-LMfire was driven by the PMIP3 suite of eight GCMs for which LGM climate simulations were available. Also provided in this archive is the result of an LPJ-LMfire run that was forced by the average climate of all GCMs (the "GCM-mean" files), and the average of each of the individual LPJ-LMfire runs over the eight LGM scenarios individually (the "LPJ-mean" files). The model simulations are provided that include the influence of human presence on the landscape (the "people" files), and in a "world without humans" scenario (the "natural" files). Finally this archive contains the preindustrial reference simulation with and without human influence ("PI_reference_people" and "PI_reference_nat", respectively). There are therefore 22 netCDF files in this archive: 8 each of LGM simulations with and without people (total 16) and the "GCM mean" simulation (2 files) and the "LPJ mean" aggregate (2 files), and finally the two preindustrial "control" simulations ("PI"), with and without humans (2 files). In addition to the LPJ-LMfire model output (netCDF files), this archive also contains a table of arboreal pollen percent calculated from pollen samples dated to the LGM at sites throughout (lgmAP.txt), and a table containing the location of archaeological sites dated to the LGM (LGM_archaeological_site_locations.txt).

Holocene precipitation change in different monsoon sub-regions (time-slices and transient data) simulated by different global climate models, with links to model results

The recently proposed global monsoon hypothesis interprets monsoon systems as part of one global-scale atmospheric overturning circulation, implying a connection between the regional monsoon systems and an in-phase behaviour of all northern hemispheric monsoons on annual timescales (Trenberth et al., 2000). Whether this concept can be applied to past climates and variability on longer timescales is still under debate, because the monsoon systems exhibit different regional characteristics such as different seasonality (i.e. onset, peak, and withdrawal). To investigate the interconnection of different monsoon systems during the pre-industrial Holocene, five transient global climate model simulations have been analysed with respect to the rainfall trend and variability in different sub-domains of the Afro-Asian monsoon region. Our analysis suggests that on millennial timescales with varying orbital forcing, the monsoons do not behave as a tightly connected global system. According to the models, the Indian and North African monsoons are coupled, showing similar rainfall trend and moderate correlation in rainfall variability in all models. The East Asian monsoon changes independently during the Holocene. The dissimilarities in the seasonality of the monsoon sub-systems lead to a stronger response of the North African and Indian monsoon systems to the Holocene insolation forcing than of the East Asian monsoon and affect the seasonal distribution of Holocene rainfall variations. Within the Indian and North African monsoon domain, precipitation solely changes during the summer months, showing a decreasing Holocene precipitation trend. In the East Asian monsoon region, the precipitation signal is determined by an increasing precipitation trend during spring and a decreasing precipitation change during summer, partly balancing each other. A synthesis of reconstructions and the model results do not reveal an impact of the different seasonality on the timing of the Holocene rainfall optimum in the different sub-monsoon systems. They rather indicate locally inhomogeneous rainfall changes and show, that single palaeo-records should not be used to characterise the rainfall change and monsoon evolution for entire monsoon sub-systems.

Model simulation output data (CCSM3) related to uplift of Africa and Himalaya/Tibet Plateau and the influence on tropical African vegetation

Hominid evolution in the late Miocene has long been hypothesized to be linked to the retreat of the tropical rainforest in Africa. One cause for the climatic and vegetation change often considered was uplift of Africa, but also uplift of the Himalaya and the Tibetan Plateau was suggested to have impacted rainfall distribution over Africa. Recent proxy data suggest that in East Africa open grassland habitats were available to the common ancestors of hominins and apes long before their divergence and do not find evidence for a closed rainforest in the late Miocene. We used the coupled global general circulation model CCSM3 including an interactively coupled dynamic vegetation module to investigate the impact of topography on African hydro-climate and vegetation. We performed sensitivity experiments altering elevations of the Himalaya and the Tibetan Plateau as well as of East and Southern Africa. The simulations confirm the dominant impact of African topography for climate and vegetation development of the African tropics. Only a weak influence of prescribed Asian uplift on African climate could be detected. The model simulations show that rainforest coverage of Central Africa is strongly determined by the presence of elevated African topography. In East Africa, despite wetter conditions with lowered African topography, the conditions were not favorable enough to maintain a closed rainforest. A discussion of the results with respect to other model studies indicates a minor importance of vegetation-atmosphere or ocean-atmosphere feedbacks and a large dependence of the simulated vegetation response on the land surface/vegetation model.

Geochemical and palynological results of IODP Hole 302-M0004A from Lomonosov Ridge, Arctic Ocean

Several episodes of abrupt and transient warming, each lasting between 50,000 and 200,000 years, punctuated the long-term warming during the Late Palaeocene and Early Eocene (58 to 51 Myr ago) epochs**1,2. These hyperthermal events, such as the Eocene Thermal Maximum 2 (ETM2) that took place about 53.5 Myr ago**2, are associated with rapid increases in atmospheric CO2 content. However, the impacts of most events are documented only locally**3,4. Here we show, on the basis of estimates from the TEX86' proxy, that sea surface temperatures rose by 3-5 °C in the Arctic Ocean during the ETM2. Dinoflagellate fossils demonstrate a concomitant freshening and eutrophication of surface waters, which resulted in euxinia in the photic zone. The presence of palm pollen implies**5 that coldest month mean temperatures over the Arctic land masses were no less than 8 °C, in contradiction of model simulations that suggest hyperthermal winter temperatures were below freezing**6. In light of our reconstructed temperature and hydrologic trends, we conclude that the temperature and hydrographic responses to abruptly increased atmospheric CO2 concentrations were similar for the ETM2 and the better-described Palaeocene-Eocene Thermal Maximum**7,8, 55.5 Myr ago.

Stable isotope records of corals from the northern Red Sea

We present a study based on X-ray chronologies and the stable isotopic composition of fossil Porites spp. corals from the northern Gulf of Aqaba (Red Sea) covering the mid-Holocene period from 5750 to 4450 14C years BP (before present). The stable oxygen and carbon isotopic compositions of five specimens reveal regular annual periodicities. Compared with modern Porites spp. from the same environment, the average seasonal delta18O amplitude of the fossil corals is higher (by ca. 0.35-0.60?), whereas annual growth rates are lower (by ca. 3.5 to 2 mm/year). This suggests stronger seasonality of sea surface temperatures and increased variability of the oxygen isotopic composition of the sea water due to changes in the precipitation and evaporation regime during the mid-Holocene. Most likely, summer monsoon rains reached the northern end of the Red Sea at that time. Average annual coral growth rates are diminished probably due to an increased input and resuspension of terrestrial debris to the shallow marine environment during more humid conditions. Our results corroborate published reports of paleodata and model simulations suggesting a northward migration of the African monsoon giving rise to increased seasonalities during the mid-Holocene over northeastern Africa and Arabia.

Monthly Tahiti coral Sr/Ca and oxygen isotope data from IODP Hole 310-M0024A

The early last glacial termination was characterized by intense North Atlantic cooling and weak overturning circulation. This interval between ~18,000 and 14,600 years ago, known as Heinrich Stadial 1, was accompanied by a disruption of global climate and has been suggested as a key factor for the termination. However, the response of interannual climate variability in the tropical Pacific (El Niño-Southern Oscillation) to Heinrich Stadial 1 is poorly understood. Here we use Sr/Ca in a fossil Tahiti coral to reconstruct tropical South Pacific sea surface temperature around 15,000 years ago at monthly resolution. Unlike today, interannual South Pacific sea surface temperature variability at typical El Niño-Southern Oscillation periods was pronounced at Tahiti. Our results indicate that the El Niño-Southern Oscillation was active during Heinrich Stadial 1, consistent with climate model simulations of enhanced El Niño-Southern Oscillation variability at that time. Furthermore, a greater El Niño-Southern Oscillation influence in the South Pacific during Heinrich Stadial 1 is suggested, resulting from a southward expansion or shift of El Niño-Southern Oscillation sea surface temperature anomalies.

Multi-sensor global compilation of mid-Holocene sea surface temperatures, based on Mg/Ca and alkenone palaeothermometry and reconstructions obtained using planktonic foraminifera and organic-walled dinoflagellate cyst

We present and examine a multi-sensor global compilation of mid-Holocene (MH) sea surface temperatures (SST), based on Mg/Ca and alkenone palaeothermometry and reconstructions obtained using planktonic foraminifera and organic-walled dinoflagellate cyst census counts. We assess the uncertainties originating from using different methodologies and evaluate the potential of MH SST reconstructions as a benchmark for climate-model simulations. The comparison between different analytical approaches (time frame, baseline climate) shows the choice of time window for the MH has a negligible effect on the reconstructed SST pattern, but the choice of baseline climate affects both the magnitude and spatial pattern of the reconstructed SSTs. Comparison of the SST reconstructions made using different sensors shows significant discrepancies at a regional scale, with uncertainties often exceeding the reconstructed SST anomaly. Apparent patterns in SST may largely be a reflection of the use of different sensors in different regions. Overall, the uncertainties associated with the SST reconstructions are generally larger than the MH anomalies. Thus, the SST data currently available cannot serve as a target for benchmarking model simulations.

Physical oceanography during Maria S. Merian cruise MSM21/2

Flemish Pass, located at the western subpolar margin, is a passage (sill depth 1200 m) that is constrained by the Grand Banks and the underwater plateau Flemish Cap. In addition to the Deep Western Boundary Current (DWBC) pathway offshore of Flemish Cap, Flemish Pass represents another southward transport pathway for two modes of Labrador Sea Water (LSW), the lightest component of North Atlantic Deep Water carried with the DWBC. This pathway avoids potential stirring regions east of Flemish Cap and deflection into the interior North Atlantic. Ship-based velocity measurements between 2009 and 2013 at 47°N in Flemish Pass and in the DWBC east of Flemish Cap revealed a considerable southward transport of Upper LSW through Flemish Pass (15-27%, -1.0 to -1.5 Sv). About 98% of the denser Deep LSW were carried around Flemish Cap as Flemish Pass is too shallow for considerable transport of Deep LSW. Hydrographic time series from ship-based measurements show a significant warming of 0.3°C/decade and a salinification of 0.03/decade of the Upper LSW in Flemish Pass between 1993 and 2013. Almost identical trends were found for the evolution in the Labrador Sea and in the DWBC east of Flemish Cap. This indicates that the long-term hydrographic variability of Upper LSW in Flemish Pass as well as in the DWBC at 47°N is dominated by changes in the Labrador Sea, which are advected southward. Fifty years of numerical ocean model simulations in Flemish Pass suggest that these trends are part of a multidecadal cycle.

Global cooling during the Eocene Oligocene Transition

About 34 million years ago, Earth's climate shifted from a relatively ice-free world to one with glacial conditions on Antarctica characterized by substantial ice sheets. How Earth's temperature changed during this climate transition remains poorly understood, and evidence for Northern Hemisphere polar ice is controversial. Here, we report proxy records of sea surface temperatures from multiple ocean localities and show that the high-latitude temperature decrease was substantial and heterogeneous. High-latitude (45 degrees to 70 degrees in both hemispheres) temperatures before the climate transition were ~20°C and cooled an average of ~5°C. Our results, combined with ocean and ice-sheet model simulations and benthic oxygen isotope records, indicate that Northern Hemisphere glaciation was not required to accommodate the magnitude of continental ice growth during this time.

(Table 1) Bioclimatic subzones and soil characteristics along the Yamal Arctic Transect, northwestern Siberia

Sustainability of tundra vegetation under changing climate on the Yamal Peninsula, northwestern Siberia, home to the world's largest area of reindeer husbandry, is of crucial importance to the local native community. An integrated investigation is needed for better understanding of the effects of soils, climate change and grazing on tundra vegetation in the Yamal region. In this study we applied a nutrient-based plant community model - ArcVeg - to evaluate how two factors (soil organic nitrogen (SON) levels and grazing) interact to affect tundra responses to climate warming across a latitudinal climatic gradient on the Yamal Peninsula. Model simulations were driven by field-collected soil data and expected grazing patterns along the Yamal Arctic Transect (YAT), within bioclimate subzones C (high arctic), D (northern low arctic) and E (southern low arctic). Plant biomass and NPP (net primary productivity) were significantly increased with warmer bioclimate subzones, greater soil nutrient levels and temporal climate warming, while they declined with higher grazing frequency. Temporal climate warming of 2 °C caused an increase of 665 g/m**2 in total biomass at the high SON site in subzone E, but only 298 g/m**2 at the low SON site. When grazing frequency was also increased, total biomass increased by only 369 g/m**2 at the high SON site in contrast to 184 g/m**2 at the low SON site in subzone E. Our results suggest that high SON can support greater plant biomass and plant responses to climate warming, while low SON and grazing may limit plant response to climate change. In addition to the first order factors (SON, bioclimate subzones, grazing and temporal climate warming), interactions among these significantly affect plant biomass and productivity in the arctic tundra and should not be ignored in regional scale studies.

Constraints on the magnitude and patterns of ocean cooling at the Last Glacial Maximum - Supplementary material

Observation-based reconstructions of sea surface temperature from relatively stable periods in the past, such as the Last Glacial Maximum, represent an important means of constraining climate sensitivity and evaluating model simulations. The first quantitative global reconstruction of sea surface temperatures during the Last Glacial Maximum was developed by the Climate Long-Range Investigation, Mapping and Prediction (CLIMAP) project in the 1970s and 1980s. Since that time, several shortcomings of that earlier effort have become apparent. Here we present an updated synthesis of sea surface temperatures during the Last Glacial Maximum, rigorously defined as the period between 23 and 19 thousand years before present, from the Multiproxy Approach for the Reconstruction of the Glacial Ocean Surface (MARGO) project. We integrate microfossil and geochemical reconstructions of surface temperatures and include assessments of the reliability of individual records. Our reconstruction reveals the presence of large longitudinal gradients in sea surface temperature in all of the ocean basins, in contrast to the simulations of the Last Glacial Maximum climate available at present.