Model-based changes in global annual mean surface temperature change (Delta T_g) and radiative forcing due to land ice albedo changes (Delta R_[LI]) over the last 5 Myr, supplementary material

It is still an open question how equilibrium warming in response to increasing radiative forcing - the specific equilibrium climate sensitivity S - depends on background climate. We here present palaeodata-based evidence on the state dependency of S, by using CO2 proxy data together with a 3-D ice-sheet-model-based reconstruction of land ice albedo over the last 5 million years (Myr). We find that the land ice albedo forcing depends non-linearly on the background climate, while any non-linearity of CO2 radiative forcing depends on the CO2 data set used. This non-linearity has not, so far, been accounted for in similar approaches due to previously more simplistic approximations, in which land ice albedo radiative forcing was a linear function of sea level change. The latitudinal dependency of ice-sheet area changes is important for the non-linearity between land ice albedo and sea level. In our set-up, in which the radiative forcing of CO2 and of the land ice albedo (LI) is combined, we find a state dependence in the calculated specific equilibrium climate sensitivity, S[CO2,LI], for most of the Pleistocene (last 2.1 Myr). During Pleistocene intermediate glaciated climates and interglacial periods, S[CO2,LI] is on average ~ 45 % larger than during Pleistocene full glacial conditions. In the Pliocene part of our analysis (2.6-5 Myr BP) the CO2 data uncertainties prevent a well-supported calculation for S[CO2,LI], but our analysis suggests that during times without a large land ice area in the Northern Hemisphere (e.g. before 2.82 Myr BP), the specific equilibrium climate sensitivity, S[CO2,LI], was smaller than during interglacials of the Pleistocene. We thus find support for a previously proposed state change in the climate system with the widespread appearance of northern hemispheric ice sheets. This study points for the first time to a so far overlooked non-linearity in the land ice albedo radiative forcing, which is important for similar palaeodata-based approaches to calculate climate sensitivity. However, the implications of this study for a suggested warming under CO2 doubling are not yet entirely clear since the details of necessary corrections for other slow feedbacks are not fully known and the uncertainties that exist in the ice-sheet simulations and global temperature reconstructions are large.

Tab. 2: Grain-sizes of cobble on the surface of sanderwurzel

yResults of 13 field investigations between 1966 and 1990 of the southwestern to eastern margin of Kötlujökull and its proglacial area are summarized with respect to sandar and their formation. Generally, the results are based on sedimentological examinations in the field and laboratory, on analyses of aerial photographs, and investigations of the glacier slope. The methods permitted a more detailed reconstruction of sandar evolution in the proglacial area of Kötlujökull since 1945, of tendencies in development and of single data going back until the last decades of the 19th century. Accordingly, there existed special periods of "flachsander"-formations with raised coarsegrained "sanderwurzels" resultant from the outbreak of subglacial meltwater tunneloutlets and other periods with "hochsander-"formations by supraglacial drainage. At present the belts of hochsanders in front of the glacier come up to more than 4 m in thickness and 1000 m in width, therefore containing perhaps more sediment direct in front of Kötlujökull than the old belts of flachsanderwurzels. In one case the explosion-like subglacial meltwater outburst combined with the genesis of a sanderwurzel could be observed for a time and is thoroughly discussed. The event is referred to the outburst of a sub- to inglacial meltwater body being under extreme hydrostatic press ures which is combined with the genesis of a new subglacial tunneloutlet as a new flachsander. Often these outbursts led to the destruction of a morainic belt more than 1000 m in width. Presumably the whole event was finished in not more than a few days. In addition to a characteristic pear-shaped form and water-moved stones up to diameters of 1 m the wurzels possess a single "main-channel" with rectangular cross-sections as far as 4 m deep and 50 m wide just as small flat channels resembling fish bones in connection with the main channel. Presumably, they have been active only in the last stage of wurzel formation. With regard to the subglacial tunnel gates long-living L-meltwater outlets are distinguished from short-living K-meltwater outlets. These are always combined with a raised coarse-grained sanderwurzel, but its meltwater discharge is generally decreasing and ceases after some years, whereas the discharge of L-meltwater outlets continues unchanged for long times (except seasonal differences). The material of flachsanders is preponderantly composed of mugearitic and andesitic cobble extending at least for some kilometres from the glacier margin, whereas the hochsanders correspond to medium to coarse sands without clay and without alternations into the direction of flow. The hochsander fans are covered with small braidet channels. Their sedimentary structures are determined by the short time changing of supraglacial meltwater discharge and the upper flow regime combined with the development of antidunes, which rule the channel-flows during the main activity periods in summer. Unlike the subglacial drainage the supraglacial drainage led to only weak effects of erosion on the glacier foreland. So the hochsanders refilled depressions of morainic areas or grew up on older flachsanderwurzels. Whereas all large flachsanders developed in front of approximate stationary glacier margins, the evolution of coherent belts of hochsanders were combined with progressive glacier fronts. On the other hand, there was obviously no evolution at all of large sandar in front of back-melting margins of Kötlujökull. Based on examinations of the glacier surface and on analyses of aerial photographs the different types of sandar are referred to different structures of the glacier snout. Finally chances of surviving of sandar in the proglacial area of Kötlujökull are shortly discussed just as the possibility of an application of the Islandic research results on Pleistocene sandar in northern Germany.

Plio-Pleistocene Mg/Ca and TEX86 sea surface temperature records of ODP Sites 184-1143 and 165-999,

The western warm pools of the Atlantic and Pacific oceans are a critical source of heat and moisture for the tropical climate system. Over the past five million years, global mean temperatures have cooled by 3-4 °C. Yet, current reconstructions of sea surface temperatures indicate that temperature in the warm pools has remained stable during this time. This stability has been used to suggest that tropical sea-surface temperatures are controlled by some sort of thermostat-like regulation. Here we reconstruct sea surface temperatures in the South China Sea, Caribbean Sea and western equatorial Pacific Ocean for the past five million years, using a combination of the Mg/Ca, TEXH86-and Uk'37 surface temperature proxies. Our data indicate that during the period of Pliocene warmth from about 5 to 2.6 million years ago, the western Pacific and western Atlantic warm pools were about 2 °C warmer than today. We suggest that the apparent lack of warming seen in the previous reconstructions was an artefact of low seawater Mg/Ca ratios in the Pliocene oceans. Taking this bias into account, our data indicate that tropical sea surface temperatures did change in conjunction with global mean temperatures. We therefore conclude that the temperature of the warm pools of the equatorial oceans during the Pliocene was not limited by a thermostat-like mechanism.

Diatom abundance and paleo-temperature reconstruction for sediment cores from the southeast Norwegian Sea

Although the pulsating nature and the abruptness of the last deglaciation are well documented in marine and land records, very few marine records have so far been able to capture the high-frequency climatic changes recorded in the Greenland ice core Dye 3. We studied high-resolution sediment cores from SE Norwegian Sea, which display a detailed climatic record during the last deglaciation comparable to that of Dye 3. Accelerator mass spectrometry age control of the cores enables us to correlate this record in detail with continental records. The results indicate that the surface waters of the SE Norwegian Sea were seasonally ice free after 13,400 B.P. The Bølling/Allerød interstadial complex (13,200-11,200 B.P.) was a climatically unstable period with changing Arctic-Subarctic conditions. This period was punctuated by four progressively more severe sea surface temperature (SST) minima: between 12,900-12,800 B.P. (BCP I); 12,500-12,400 B.P. (BCP II); 12,300-12,000 B.P. (OD I); and 11,800-11,500 B.P. (OD II). The Younger Dryas (YD) (11,200-10,200 B.P.) represents the severest and most prolonged cold episode of this series of climatic deteriorations. It was bounded by very rapid SST changes and characterized by Arctic-Polar conditions. The first true warm Atlantic water incursion to the SE Norwegian Sea took place around 10,100 B.P., followed by a brief cooler condition between 9900-9600 B.P. (YD II). The early Holocene climatic optimum occurred between 8000-5000 B.P. A conceptual model is proposed where meltwater fluxes are suggested to cause the observed instability in the SST record.

Temperature reconstruction of sediment cores from the northeastern Arabian Sea

In order to reconstruct the monsoonal variability during the late Holocene we investigated a complete, annually laminated sediment record from the oxygen minimum zone (OMZ) off Pakistan for oxygen isotopes of planktic foraminifera and alkenone-derived sea surface temperatures (SST). Significant SST changes of up to 3°C which cannot be explained by changes in the alkenone-producing coccolithophorid species (inferred from the Gephyrocapsa oceanica / Emiliania huxleyi ratio) suggest that SST changes are driven by changes in the monsoon strength. Our high-(decadal)-resolution data indicate that the late Holocene in the northeastern Arabian Sea was not characterized by a stable uniform climate, as inferred from the Greenland ice cores, but by variations in the dominance of the SW monsoon conditions with significant effects on temperatures. Highest SST fluctuations of up to 3.0°C and 2.5°C were observed for the time interval from 4600 to 3300 years B.P. and during the past 500 years. The significant, short-term SST changes during the past 500 years might be related to climatic instabilities known from the northern latitudes ("Little Ice Age") and confirm global effects. Surface salinity values, reconstructed from delta18O records after correction for temperature-related oxygen isotope fractionation, suggest that in general, the past 5000 years were characterized by higher-than-recent evaporation and more intense SW monsoon conditions. However, between 4600 and 3700 years B.P., evaporation dropped, SW monsoon weakened, and NE monsoon conditions were comparatively enhanced. For the past 1500 years we infer strongly fluctuating monsoon conditions. Comparisons of reconstructed salinity records with ice accumulation data from published Tibetan ice core and Tibetan tree ring width data reveal that during the past 2000 years, enhanced evaporation in the northeastern Arabian Sea correlates with periods of increased ice accumulation in Tibet, and vice versa. This suggests a strong climatic relationship between both monsoon-controlled areas.

Temperature and salinity reconstruction for the Last Interglacial Period in the North Atlantic Ocean

Eight deep-sea sediment cores from the North Atlantic Ocean ranging from 31° to 72°N are studied to reconstruct the meridional gradients in surface hydrographic conditions during the interval of minimum ice volume within the last interglacial period. Using benthic foraminiferal ?18O measurements and estimates of Sea Surface Temperature (SST) and Sea Surface Salinity (SSS), we show that summer SSTs and SSSs decreased gradually during the interval of minimum ice volume at high-latitude sites (52°-72°N) whereas they were stable or increased during the same time period at low-latitude sites (31°-41°N). This increase in meridional gradients of SSTs and SSSs may have been due to changes in the latitudinal distribution of summer and annual-average insolation and associated oceanic and atmospheric feedbacks. These trends documented for the Eemian ice volume minimum period are similar to corresponding changes observed during the Holocene and may have had a similar origin.

Clay mineraology and grain size composition of surface sediments from South Pacific Ocean

The reconstruction of low-latitude ocean-atmosphere interactions is one of the major issues of (paleo-)environmental studies. The trade winds, extending over 20° to 30° of latitude in both hemispheres, between the subtropical highs and the intertropical convergence zone, are major components of the atmospheric circulation and little is known about their long-term variability on geological time-scales, in particular in the Pacific sector. We present the modern spatial pattern of eolian-derived marine sediments in the eastern equatorial and subtropical Pacific (10°N to 25°S) as a reference data set for the interpretation of SE Pacific paleo-dust records. The terrigenous silt and clay fractions of 75 surface sediment samples have been investigated for their grain-size distribution and clay-mineral compositions, respectively, to identify their provenances and transport agents. Dust delivered to the southeast Pacific from the semi- to hyper-arid areas of Peru and Chile is rather fine-grained (4-8 µm) due to low-level transport within the southeast trade winds. Nevertheless, wind is the dominant transport agent and eolian material is the dominant terrigenous component west of the Peru-Chile Trench south of ~ 5°S. Grain-size distributions alone are insufficient to identify the eolian signal in marine sediments due to authigenic particle formation on the sub-oceanic ridges and abundant volcanic glass around the Galapagos Islands. Together with the clay-mineral compositions of the clay fraction, we have identified the dust lobe extending from the coasts of Peru and Chile onto Galapagos Rise as well as across the equator into the doldrums. Illite is a very useful parameter to identify source areas of dust in this smectite-dominated study area.

Diatom biostratigraphy and sea ice concentration reconstruction from core DA06-139G, Vaigat Strait in Disko Bugt, West Greenland

A diatom-based sea-ice concentration (SIC) transfer function is developed using 72 surface samples from west of Greenland and around Iceland, and through comparison with the associated modern SIC. Canonical correspondence analysis on surface sediment diatoms and monthly average of SIC reveals that April SIC is the most important environmental factor controlling the distribution of diatoms in the area, and permits the development of a diatom-based SIC transfer function. The consistency between reconstructed SIC based on diatoms from West Greenland and the instrumental and documentary data during the last ~75 years demonstrates that the diatom-based SIC reconstruction is reliable for studying the palaeoceanography off West Greenland. Relatively warm conditions with strong influence of the Irminger Current (IC) are indicated for the early part of the record (~5000-3860 cal. yr BP), corresponding in time to the latest part of the Holocene Thermal Maximum. The April SIC oscillated around the mean value between 3860 and 1510 cal. yr BP and was above mean afterwards, particularly during the time interval 1510-1120 cal. yr BP and after 650 cal. yr BP, indicating more extensive sea-ice cover in Disko Bugt. A high degree of consistency between the reconstructed April SIC and changes in the diatom species suggests that the sea-ice condition in Disko Bugt is strongly influenced by variations in the relative strength of two components of the West Greenland Current, i.e. the cold East Greenland Current and the relatively warm IC.

Sea surafce temperature reconstruction for sediment core SU81-18

Reconstructing the impact of Heinrich events outside the main belt of ice rafting is crucial to understanding the underlying causes of these abrupt climatic events. A high-resolution study of a marine sediment core from the Iberian margin demonstrates that this midlatitude area was strongly affected both by cooling and advection of low-salinity arctic water masses during the last three Heinrich events. These paleoclimatic time series reveal the internal complexity of each of the last three Heinrich events and illustrate the value of parallel studies of the organic and inorganic fractions of the sediments.

Geochemistry, sediment components, and paleoproductivity reconstruction for the Miocene to early Pliocene sediments

We reconstruct paleoproductivity at three sites in the Atlantic Ocean (Ocean Drilling Program Sites 982, 925, and 1088) to investigate the presence and extent of the late Miocene to early Pliocene 'biogenic bloom' from 9 to 3 Ma. Our approach involves construction of multiple records including benthic foraminiferal and CaCO3 accumulation rates, Uvigerina counts, dissolution proxies, and geochemical tracers for biogenic and detrital fluxes. This time interval also contains the so-called late Miocene carbon isotope shift, a well-known decrease in benthic foraminiferal d13C values. We find that the timing of paleoproductivity maxima differs among the three sites. At Site 982 (North Atlantic), benthic foraminifera and CaCO3 accumulation were both at a maximum at ~5 Ma, with smaller peaks at ~6 Ma. The paleoproductivity maximum was centered earlier (~6.6-6.0 Ma) in the tropical Atlantic (Site 925). In the South Atlantic (Site 1088), paleoproductivity increased even earlier, between 8.2 Ma and 6.2 Ma, and remained relatively high until ~5.4 Ma. We note that there is some overlap between the interval of maximum productivity between Sites 925 and 1088, as well as the minor productivity increase at Site 982. We conclude that the paleoproductivity results support hypotheses aiming to place the biogenic bloom into a global context of enhanced productivity. In addition, we find that at all three sites the d13C shift is accompanied by carbonate dissolution. This observation is consistent with published studies that have sought a relationship between the late Miocene carbon isotope shift and carbonate preservation.

Automated composite depth scale reconstruction of ODP Leg 138 and 154 sites

A composite section, which reconstructs a continuous stratigraphic record from cores of multiple nearby holes, and its associated composite depth scale are important tools for analyzing sediment recovered from a drilling site. However, the standard technique for creating composite depth scales on drilling cruises does not correct for depth distortion within each core. Additionally, the splicing technique used to create composite sections often results in a 10-15% offset between composite depths and measured drill depths. We present a new automated compositing technique that better aligns stratigraphy across holes, corrects depth offsets, and could be performed aboard ship. By analyzing 618 cores from seven Ocean Drilling Program (ODP) sites, we estimate that ?80% of the depth offset in traditional composite depth scales results from core extension during drilling and extraction. Average rates of extension are 12.4 ± 1.5% for calcareous and siliceous cores from ODP Leg 138 and 8.1 ± 1.1% for calcareous and clay-rich cores from ODP Leg 154. Also, average extension decreases as a function of depth in the sediment column, suggesting that elastic rebound is not the dominant extension mechanism.