Estimation of sea surface temperatures and salinity of ODP Site 177-1089 (Appendix A)

ODP Site 1089 is optimally located in order to monitor the occurrence of maxima in Agulhas heat and salt spillage from the Indian to the Atlantic Ocean. Radiolarian-based paleotemperature transfer functions allowed to reconstruct the climatic history for the last 450 kyr at this location. A warm sea surface temperature anomaly during Marine Isotope Stage (MIS) 10 was recognized and traced to other oceanic records along the surface branch of the global thermohaline (THC) circulation system, and is particularly marked at locations where a strong interaction between oceanic and atmospheric overturning cells and fronts occurs. This anomaly is absent in the Vostok ice core deuterium, and in oceanic records from the Antarctic Zone. However, it is present in the deuterium excess record from the Vostok ice core, interpreted as reflecting the temperature at the moisture source site for the snow precipitated at Vostok Station. As atmospheric models predict a subtropical Indian source for such moisture, this provides the necessary teleconnection between East Antarctica and ODP Site 1089, as the subtropical Indian is also the source area of the Agulhas Current, the main climate agent at our study location. The presence of the MIS 10 anomaly in the delta13C foraminiferal records from the same core supports its connection to oceanic mechanisms, linking stronger Agulhas spillover intensity to increased productivity in the study area. We suggest, in analogy to modern oceanographic observations, this to be a consequence of a shallow nutricline, induced by eddy mixing and baroclinic tide generation, which are in turn connected to the flow geometry, and intensity, of the Agulhas Current as it flows past the Agulhas Bank. We interpret the intensified inflow of Agulhas Current to the South Atlantic as responding to the switch between lower and higher amplitude in the insolation forcing in the Agulhas Current source area. This would result in higher SSTs in the Cape Basin during the glacial MIS 10, due to the release into the South Atlantic of the heat previously accumulating in the subtropical and equatorial Indian and Pacific Ocean. If our explanation for the MIS 10 anomaly in terms of an insolation variability switch is correct, we might expect that a future Agulhas SSST anomaly event will further delay the onset of next glacial age. In fact, the insolation forcing conditions for the Holocene (the current interglacial) are very similar to those present during MIS 11 (the interglacial preceding MIS 10), as both periods are characterized by a low insolation variability for the Agulhas Current source area. Natural climatic variability will force the Earth system in the same direction as the anthropogenic global warming trend, and will thus lead to even warmer than expected global temperatures in the near future.

(Table 1) Age determination of sediment core AAS-62/2

Variations in the Indian summer monsoon (ISM) intensity during the last 16.7?ka have been studied using organic carbon (Corg), d15N of sedimentary organic matter, CaCO3, sediment texture, d18OC, and Mg/Ca-derived sea surface temperature, d18O of sea water and sea surface salinity, in a 14C-dated sediment core from the eastern Arabian Sea. The d18O in water and planktonic foraminifera shells off the central west coast of India may be controlled by the ISM intensity as this area receives high precipitation and land runoff. Also, the Corg and CaCO3 contents of sediments and d15N of sedimentary organic matter may be linked to ISM-induced productivity and denitrification. The results of the present study reveal that between 16 and 15.2 ka BP, the ISM was weak with minor fluctuations and started intensifying around 15.2 ka BP, at the onset of the Bølling-Ållerød (B-A) event. The B-A event is characterized by higher water column denitrification rates comparable to the present day. The ISM signatures observed in the d18OC record of B-A event compare well with those from Timta cave of the western Himalayas and also the Asian summer monsoon signatures from the Hulu caves in China and warming signatures in Greenland Ice Sheet Project 2 (GISP2) suggesting atmospheric teleconnections through Intertropical Convergence Zone. The boundary between the Younger Dryas and the Holocene is discernible with small episodes of abrupt events of increased ISM intensity. This decrease in d18OC values at ~11.8 ka BP is contemporary with June solar insolation maximum at 30° north and the increase in methane in the GISP2 ice core supporting episodes of warmer climate and increase in ISM intensity. The ISM seems to have been most stable between 7 and 5.6 ka BP. The core exhibits periodicity of 500 years that is comparable to the Atlantic water formation and the Chinese monsoon.

Downcore rock-magnetic properties of ICDP Core 5045-1, Lake Ohrid

The bulk magnetic mineral record from Lake Ohrid, spanning the past 637 kyr, reflects large-scale shifts in hydrological conditions, and, superimposed, a strong signal of environmental conditions on glacial-interglacial and millennial timescales. A shift in the formation of early diagenetic ferrimagnetic iron sulfides to siderites is observed around 320 ka. This change is probably associated with variable availability of sulfide in the pore water. We propose that sulfate concentrations were significantly higher before ~320 ka, due to either a higher sulfate flux or lower dilution of lake sulfate due to a smaller water volume. Diagenetic iron minerals appear more abundant during glacials, which are generally characterized by higher Fe/Ca ratios in the sediments. While in the lower part of the core the ferrimagnetic sulfide signal overprints the primary detrital magnetic signal, the upper part of the core is dominated by variable proportions of high- to low-coercivity iron oxides. Glacial sediments are characterized by high concentration of high-coercivity magnetic minerals (hematite, goethite), which relate to enhanced erosion of soils that had formed during preceding interglacials. Superimposed on the glacial-interglacial behavior are millennial-scale oscillations in the magnetic mineral composition that parallel variations in summer insolation. Like the processes on glacial-interglacial timescales, low summer insolation and a retreat in vegetation resulted in enhanced erosion of soil material. Our study highlights that rock-magnetic studies, in concert with geochemical and sedimentological investigations, provide a multi-level contribution to environmental reconstructions, since the magnetic properties can mirror both environmental conditions on land and intra-lake processes.

Expanded GHOST database

In this study we review a global set of alkenone- and foraminiferal Mg/Ca-derived sea surface temperatures (SST) records from the Holocene and compare them with a suite of published Eemian SST records based on the same approach. For the Holocene, the alkenone SST records belong to the actualized GHOST database (Kim, J.-H., Schneider R.R., 2004). The actualized GHOST database not only confirms the SST changes previously described but also documents the Holocene temperature evolution in new oceanic regions such as the Northwestern Atlantic, the eastern equatorial Pacific, and the Southern Ocean. A comparison of Holocene SST records stemming from the two commonly applied paleothermometry methods reveals contrasting - sometimes divergent - SST evolution, particularly at low latitudes where SST records are abundant enough to infer systematic discrepancies at a regional scale. Opposite SST trends at particular locations could be explained by out-of-phase trends in seasonal insolation during the Holocene. This hypothesis assumes that a strong contrast in the ecological responses of coccolithophores and planktonic foraminifera to winter and summer oceanographic conditions is the ultimate reason for seasonal differences in the origin of the temperature signal provided by these organisms. As a simple test for this hypothesis, Eemian SST records are considered because the Holocene and Eemian time periods experienced comparable changes in orbital configurations, but had a higher magnitude in insolation variance during the Eemian. For several regions, SST changes during both interglacials were of a similar sign, but with higher magnitudes during the Eemian as compared to the Holocene. This observation suggests that the ecological mechanism shaping SST trends during the Holocene was comparable during the penultimate interglacial period. Although this "ecology hypothesis" fails to explain all of the available results, we argue that any other mechanism would fail to satisfactorily explain the observed SST discrepancies among proxies.

Estimates of annual mean surface temperatures in the late Quaternary of the tropical Atlantic

At least two modes of glacial-interglacial climate change have existed within the tropical Atlantic Ocean during the last 20,000 years. The first mode (defined by cold glacial and warm interglacial conditions) occurred symmetrically north and south of the equator and dominated the eastern boundary currents and tropical upwelling areas. This pattern suggests that mode 1 is driven by a glacial modification of surface winds in both hemispheres. The second mode of oceanic climate change, defined by temperature extremes centered on the deglaciation, was hemispherically asymmetrical, with the northern tropical Atlantic relatively cold and the southern tropical Atlantic relatively warm during deglaciation. A likely cause for this pattern of variation is a reduction of the presently northward cross-equatorial heat flux during deglaciation. No single mechanism accounts for all the data. Potential contributors to oceanic climate changes are linkage to high-latitude climates, modification of monsoonal winds by ice sheet and/or insolation changes, atmospheric CO2 and greenhouse effects, indirect effects of glacial meltwater, and variations in thermohaline overturn of the oceans.

(Table S1) Chronostratigraphic constrains for sediment core AND1-1B

Thirty years after oxygen isotope records from microfossils deposited in ocean sediments confirmed the hypothesis that variations in the Earth's orbital geometry control the ice ages (Hays et al., 1976, doi:10.1126/science.194.4270.1121), fundamental questions remain over the response of the Antarctic ice sheets to orbital cycles (Raymo and Huybers, 2008, doi:10.1038/nature06589). Furthermore, an understanding of the behaviour of the marine-based West Antarctic ice sheet (WAIS) during the 'warmer-than-present' early-Pliocene epoch (~5-3 Myr ago) is needed to better constrain the possible range of ice-sheet behaviour in the context of future global warming (Solomon et al., 2007). Here we present a marine glacial record from the upper 600 m of the AND-1B sediment core recovered from beneath the northwest part of the Ross ice shelf by the ANDRILL programme and demonstrate well-dated, ~40-kyr cyclic variations in ice-sheet extent linked to cycles in insolation influenced by changes in the Earth's axial tilt (obliquity) during the Pliocene. Our data provide direct evidence for orbitally induced oscillations in the WAIS, which periodically collapsed, resulting in a switch from grounded ice, or ice shelves, to open waters in the Ross embayment when planetary temperatures were up to ~3° C warmer than today ( Kim and Crowley, 2000, doi:10.1029/1999PA000459) and atmospheric CO2 concentration was as high as ~400 p.p.m.v. (van der Burgh et al., 1993, doi:10.1126/science.260.5115.1788, Raymo et al., 1996, doi:10.1016/0377-8398(95)00048-8). The evidence is consistent with a new ice-sheet/ice-shelf model (Pollard and DeConto, 2009, doi:10.1038/nature07809) that simulates fluctuations in Antarctic ice volume of up to +7 m in equivalent sea level associated with the loss of the WAIS and up to +3 m in equivalent sea level from the East Antarctic ice sheet, in response to ocean-induced melting paced by obliquity. During interglacial times, diatomaceous sediments indicate high surface-water productivity, minimal summer sea ice and air temperatures above freezing, suggesting an additional influence of surface melt (Huybers, 2006, doi:10.1126/science.1125249) under conditions of elevated CO2.

Magnetic susceptibility dataset from the Danian Contessa Highway section, Gubbio, Italy

We studied a high-resolution multiproxy data set, including magnetic susceptibility (MS), CaCO3 content, and stable isotopes (d18O and d13C), from the stratigraphic interval covering the uppermost Maastrichtian and the lower Danian, represented by the pelagic limestones of the Scaglia Rossa Formation continuously exposed in the classic sections of the Bottaccione Gorge and the Contessa Highway near Gubbio, Italy. Variations in all the proxy series are periodic and reflect astronomically forced climate changes (i.e., Milankovitch cycles). In particular, the MS proxy reflects variations in the terrigenous dust input in this pelagic, deep-marine environment. We speculate that the dust is mainly eolian in origin and that the availability and transport of dust are influenced by variations in the vegetation cover on the Maastrichtian-Paleocene African or Asian zone, which were respectively located at tropical to subtropical latitudes to the south or far to the east of the western Tethyan Umbria-Marche Basin, and were characterized by monsoonal circulation. The dynamics of monsoonal circulation are known to be strongly dependent on precession-driven and obliquity-driven changes in insolation. We propose that a threshold mechanism in the vegetation coverage may explain eccentricity-related periodicities in the terrigenous eolian dust input. Other mechanisms, both oceanic and terrestrial, that depend on the precession amplitude modulated by eccentricity, can be evoked together with the variation of dust influx in the western Tethys to explain the detected eccentricity periodicity in the d13C record. Our interpretations of the d18O and MS records suggest a warming event ~400 k.y. prior to the Cretaceous-Paleogene (K-Pg) boundary, and a period of climatic and environmental instability in the earliest Danian. Based on these multiproxy phase relationships, we propose an astronomical tuning for these sections; this leads us to an estimate of the timing and duration of several late Maastrichtian and Danian biostratigraphic and magnetostratigraphic events.

Age determination of North Atlantic sediment cores

Available overwash records from coastal barrier systems document significant variability in North Atlantic hurricane activity during the late Holocene. The same climate forcings that may have controlled cyclone activity over this interval (e.g., the West African Monsoon, El Niño-Southern Oscillation (ENSO)) show abrupt changes around 6000 yrs B.P., but most coastal sedimentary records do not span this time period. Establishing longer records is essential for understanding mid-Holocene patterns of storminess and their climatic drivers, which will lead to better forecasting of how climate change over the next century may affect tropical cyclone frequency and intensity. Storms are thought to be an important mechanism for transporting coarse sediment from shallow carbonate platforms to the deep-sea, and bank-edge sediments may offer an unexplored archive of long-term hurricane activity. Here, we develop this new approach, reconstructing more than 7000 years of North Atlantic hurricane variability using coarse-grained deposits in sediment cores from the leeward margin of the Great Bahama Bank. High energy event layers within the resulting archive are (1) broadly correlated throughout an offbank transect of multi-cores, (2) closely matched with historic hurricane events, and (3) synchronous with previous intervals of heightened North Atlantic hurricane activity in overwash reconstructions from Puerto Rico and elsewhere in the Bahamas. Lower storm frequency prior to 4400 yrs B.P. in our records suggests that precession and increased NH summer insolation may have greatly limited hurricane potential intensity, outweighing weakened ENSO and a stronger West African Monsoon-factors thought to be favorable for hurricane development.

(Table S1) Opal content of sediment core OXMZ01MV-GC31

The North American monsoon (NAM), an onshore wind shift occurring between July and September, has evolved in character during the Holocene largely due to changes in Northern Hemisphere insolation. Published paleoproxy and modeling studies suggest that prior to ~8000 cal years BP, the NAM affected a broader region than today, extending westward into the Mojave Desert of California. Holocene proxy SST records from the Gulf of California (GoC) and the adjacent Pacific provide constraints for this changing NAM climatology. Prior to ~8000 cal years BP, lower GoC SSTs would not have fueled northward surges of tropical moisture up the GoC, which presently contribute most of the monsoon precipitation to the western NAM region. During the early Holocene, the North Pacific High was further north and SSTs in the California Current off Baja California were warmer, allowing monsoonal moisture flow from the subtropical Pacific to take a more direct, northwesterly trajectory into an expanded area of the southwestern U.S. west of 114°W. A new upwelling record off southwest Baja California reveals that enhanced upwelling in the California Current beginning at ~7500 cal year BP may have triggered a change in NAM climatology, focusing the geographic expression of NAM in the southwest USA into its modern core region east of ~114°W, in Arizona and New Mexico. Holocene proxy precipitation records from the southwestern U.S. and northwestern Mexico, including lakes, vegetation/pollen, and caves are reviewed and found to be largely supportive of this hypothesis of changing Holocene NAM climatology.

(Table 1) Age determination of Barents Sea sediment cores

he separate roles of oceanic heat advection and orbital forcing on influencing early Holocene temperature variability in the eastern Nordic Seas is investigated. The effect of changing orbital forcing on the ocean temperatures is tested using the 1DICE model, and the 1DICE results are compared with new and previously published temperature reconstructions from a transect of five cores located underneath the pathway of Atlantic water, from the Faroe-Shetland Channel in the south to the Barents Sea in the north. The stronger early Holocene summer insolation at high northern latitudes increased the summer mixed layer temperatures, however, ocean temperatures underneath the summer mixed layer did not increase significantly. The absolute maximum in summer mixed layer temperatures occurred between 9 and 6 ka BP, representing the Holocene Thermal Maximum in the eastern Nordic Seas. In contrast, maximum in northward oceanic heat transport through the Norwegian Atlantic Current occurred approximately 10 ka BP. The maximum in oceanic heat transport at 10 ka BP occurred due to a major reorganization of the Atlantic Ocean circulation, entailing strong and deep rejuvenation of the Atlantic Meridional Overturning Circulation, combined with changes in the North Atlantic gyre dynamic causing enhanced transport of heat and salt into the Nordic Seas.

Multi-proxy records of Two-Yurts Lake

Within the scope of Russian-German palaeoenvironmental research, Two-Yurts Lake (TYL, Dvuh-Yurtochnoe in Russian) was chosen as the main scientific target area to decipher Holocene climate variability on Kamchatka. The 5x2 km large and 26 m deep lake is of proglacial origin and situated on the eastern flank of Sredinny Ridge at the northwestern end of the Central Kamchatka Valley, outside the direct influence of active volcanism. Here, we present results of a multi-proxy study on sediment cores, spanning about the last 7000 years. The general tenor of the TYL record is an increase in continentality and winter snow cover in conjunction with a decrease in temperature, humidity, and biological productivity after 5000-4500 cal yrs BP, inferred from pollen and diatom data and the isotopic composition of organic carbon. The TYL proxy data also show that the late Holocene was punctuated by two colder spells, roughly between 4500 and 3500 cal yrs BP and between 1000 and 200 cal yrs BP, as local expressions of the Neoglacial and Little Ice Age, respectively. These environmental changes can be regarded as direct and indirect responses to climate change, as also demonstrated by other records in the regional terrestrial and marine realm. Long-term climate deterioration was driven by decreasing insolation, while the short-term climate excursions are best explained by local climatic processes. The latter affect the configuration of atmospheric pressure systems that control the sources as well as the temperature and moisture of air masses reaching Kamchatka.

Alkenone unsaturation index and alkenones in ODP Hole 161-977A

A number of essential elements closely related to each other are involved in the Earth's climatic system. The temporal and spatial distribution of insolation determines wind patterns and the ocean's thermohaline pump. In turn, these last two are directly linked to the extension and retreat of marine and continental ice and to the chemistry of the atmosphere and the ocean. The variability of these elements may trigger, amplify, sustain or globalize rapid climatic changes. Paleoclimatic oscillations have been identified in this thesis by using fossil organic compounds synthesized by marine and terrestrial flora. High sedimentation rate deposits at the Barents and the Iberian peninsula continental margins were chosen in order to estimate the climatic changes on centennial time resolution. At the Barents margin, the sediment recovered was up to 15,000 years old (unit ''a'', from latin ''annos'') (M23258; west of the Bjørnøya island). At the Iberian margin, the sediment cores studied covered a wide range of time spans: up to 115,000 a (MD99-2343; north of the Minorca island), up to 250,000 a (ODP-977A; Alboran basin) and up to 420,000 a (MD01-2442, MD01-2443, MD01-2444, MD01-2445; close to the Tagus abyssal plain). At the northern site, inputs containing marine, continental and ancient reworked organic matter provided a detailed reconstruction of climate history at the time of the final retreat of the Barents ice sheet. At the western Barents continental slope, warm climatic conditions were observed during the early Holocene (~from 8,650 a to 5,240 a ago); in contrast, an apparent long-term cooling trend occurred in the late Holocene (~from 5,240 a to 760 a ago), in consistence with other paleoarchives from northern and southern European latitudes. The Iberian margin sites, which were never covered with large ice sheets, preserved exceptionally complete sequences of rapid events during ice ages hitherto not studied in such great detail: during the last glacial (~from 70,900 a to 11,800 a ago), the second glacial (~from 189,300 a to 127,500 a ago), the third ice age (~from 278,600 a to 244,800 a ago) and the fourth (~from 376,300 a to 337,500 a ago). In this thesis, crucial research questions were brought up concerning the severity of different glacial periods, the intensity and rates of the recorded oscillations and the long distance connections related to rapid climate change. The data obtained provide a sound basis to further research on the mechanisms involved in this rapid climate variability. An essential point of the research was the evidence that, over the past 420,000 a, at the whole Iberian margin, warm and stable long periods similar to the Holocene always ended abruptly in few centuries after a gradual deterioration of climate conditions. The detailed estimate of past climate variability provides clues to the natural end of the present warm period. Returning to an ice age in European lands would be exacerbated by a number of factors: a lack of differential solar heating between northern and southern north Atlantic latitudes, enhanced evaporation at low latitudes, and an increase in snowfall or iceberg discharges at northern regions. It must be emphasized that all climatic oscillations observed in this thesis were caused by forces of nature, i.e. the last two centuries were not taken into consideration.

Holocene experiment with coupled atmosphere-ocean-model ECHAM5/MPI-OM

Changes in the Earth's orbit lead to changes in the seasonal and meridional distribution of insolation. We quantify the influence of orbitally induced changes on the seasonal temperature cycle in a transient simulation of the last 6000 years - from the mid-Holocene to today - using a coupled atmosphere-ocean general circulation model (ECHAM5/MPI-OM) including a land surface model (JSBACH). The seasonal temperature cycle responds directly to the insolation changes almost everywhere. In the Northern Hemisphere, its amplitude decreases according to an increase in winter insolation and a decrease in summer insolation. In the Southern Hemisphere, the opposite is true. Over the Arctic Ocean, decreasing summer insolation leads to an increase in sea-ice cover. The insulating effect of sea ice between the ocean and the atmosphere leads to decreasing heat flux and favors more "continental" conditions over the Arctic Ocean in winter, resulting in strongly decreasing temperatures. Consequently, there are two competing effects: the direct response to insolation changes and a sea-ice insulation effect. The sea-ice insulation effect is stronger, and thus an increase in the amplitude of the seasonal temperature cycle over the Arctic Ocean occurs. This increase is strongest over the Barents Shelf and influences the temperature response over northern Europe. We compare our modeled seasonal temperatures over Europe to paleo reconstructions. We find better agreements in winter temperatures than in summer temperatures and better agreements in northern Europe than in southern Europe, since the model does not reproduce the southern European Holocene summer cooling inferred from the paleo reconstructions. The temperature reconstructions for northern Europe support the notion of the influence of the sea-ice insulation effect on the evolution of the seasonal temperature cycle.

(Table 1) Age determination of North Atlantic sediment cores

Stable isotope, trace metal, alkenone paleothermometry, and radiocarbon methods have been applied to sediment cores in the western subpolar North Atlantic between Hudson Strait and Cape Hatteras to reveal the history of climate in that region over the past ~11 kyr. We focus on cores from the Laurentian Fan, which is known to have rapid and continuous accumulation of hemipelagic sediment. Although results among our various proxy data are not always in agreement, the weight of the evidence (alkenone sea surface temperature (SST), d18O and abundance of Globigerinoides ruber) indicates a continual cooling of surface waters over Laurentian Fan, from about 18°C in the early Holocene to about 8°C today. Alternatively, Mg/Ca data on planktonic foraminifera indicate no systematic change in Holocene SST. The inferred long-term decrease in SST was probably driven by decreasing seasonality of Northern Hemisphere insolation. Two series of proxy data show the gradual cooling was interrupted by a two-step cold pulse that began 8500 years ago, and lasted about 700 years. Although this event is associated with the final deglaciation of Hudson Bay, there is no d18O minimum anywhere in the Labrador Sea, yet there is some evidence for it as far south as Cape Hatteras. Finally, although the 8200 year B.P. event has been implicated in decreasing North Atlantic ventilation, and hence widespread temperature depression on land and at sea, we find inconsistent evidence for a change at that time in deep ocean nutrient content at ~4 km water depth.

Stable isotope record of foraminifera from the Barents Sea

We measured the oxygen isotopic composition of planktonic and benthic foraminifera in three cores collected at key positions to reconstruct the paleoceanography of the Barents Sea: core ASV 880 on the path of the northern branch of Atlantic water inflowing from the Arctic Ocean, core ASV 1200 in the central basin near the polar front, and core ASV 1157 in the main area of brine formation. Modern seawater d18O measurements show that far from the coast, d18O variations are linearly linked to the salinity changes associated with sea ice melting. The foraminifer d18O records are dated by 14C measurements performed on mollusk shells, and they provide a detailed reconstruction of the paleoceanographic evolution of the Barents Sea during the Holocene. Four main steps were recognized: the terminal phase of the deglaciation with melting of the main glaciers, which were located on the surrounding continent and islands, the short thermal optimum from 7.8 ka B.P. to 6.8 ka B.P., a cold mid-Holocene phase with a large reduction of the inflow of Atlantic water, and the inception of the modern hydrological pattern by 4.7 ka B.P. Brine water formation was active during the whole Holocene. The paleoclimatic evolution of the Barents Sea was driven by both high-latitude summer insolation and the intensity of the Atlantic water inflow.