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.

Single non-normalized data of electron probe analyses of all glass shard samples from the Seward Peninsula and the Lipari obsidian reference standard

Permafrost degradation influences the morphology, biogeochemical cycling and hydrology of Arctic landscapes over a range of time scales. To reconstruct temporal patterns of early to late Holocene permafrost and thermokarst dynamics, site-specific palaeo-records are needed. Here we present a multi-proxy study of a 350-cm-long permafrost core from a drained lake basin on the northern Seward Peninsula, Alaska, revealing Lateglacial to Holocene thermokarst lake dynamics in a central location of Beringia. Use of radiocarbon dating, micropalaeontology (ostracods and testaceans), sedimentology (grain-size analyses, magnetic susceptibility, tephra analyses), geochemistry (total nitrogen and carbon, total organic carbon, d13Corg) and stable water isotopes (d18O, dD, d excess) of ground ice allowed the reconstruction of several distinct thermokarst lake phases. These include a pre-lacustrine environment at the base of the core characterized by the Devil Mountain Maar tephra (22 800±280 cal. a BP, Unit A), which has vertically subsided in places due to subsequent development of a deep thermokarst lake that initiated around 11 800 cal. a BP (Unit B). At about 9000 cal. a BP this lake transitioned from a stable depositional environment to a very dynamic lake system (Unit C) characterized by fluctuating lake levels, potentially intermediate wetland development, and expansion and erosion of shore deposits. Complete drainage of this lake occurred at 1060 cal. a BP, including post-drainage sediment freezing from the top down to 154 cm and gradual accumulation of terrestrial peat (Unit D), as well as uniform upward talik refreezing. This core-based reconstruction of multiple thermokarst lake generations since 11 800 cal. a BP improves our understanding of the temporal scales of thermokarst lake development from initiation to drainage, demonstrates complex landscape evolution in the ice-rich permafrost regions of Central Beringia during the Lateglacial and Holocene, and enhances our understanding of biogeochemical cycles in thermokarst-affected regions of the Arctic.

Baseline characteristics of climate, permafrost, and land cover from Samoylov Island, Lena River Delta, Siberia

Samoylov Island is centrally located within the Lena River Delta at 72° N, 126° E and lies within the Siberian zone of continuous permafrost. The landscape on Samoylov Island consists mainly of late Holocene river terraces with polygonal tundra, ponds and lakes, and an active floodplain. The island has been the focus of numerous multidisciplinary studies since 1993, which have focused on climate, land cover, ecology, hydrology, permafrost and limnology. This paper aims to provide a framework for future studies by describing the characteristics of the island's meteorological parameters (temperature, radiation and snow cover), soil temperature, and soil moisture. The land surface characteristics have been described using high resolution aerial images in combination with data from ground-based observations. Of note is that deeper permafrost temperatures have increased between 0.3 to 1.3 °C over the last five years. However, no clear warming of air and active layer temperatures is detected since 1998, though winter air temperatures during recent years have not been as cold as in earlier years.

(Table S1) Radiocarbon dates of organic matter from Lena Delta ice wedges

Relative to the past 2,000 years, the Arctic region has warmed significantly over the past few decades. However, the evolution of Arctic temperatures during the rest of the Holocene is less clear. Proxy reconstructions, suggest a long-term cooling trend throughout the mid- to late Holocene, whereas climate model simulations show only minor changes or even warming. Here we present a record of the oxygen isotope composition of permafrost ice wedges from the Lena River Delta in the Siberian Arctic. The isotope values, which reflect winter season temperatures, became progressively more enriched over the past 7,000 years, reaching unprecedented levels in the past five decades. This warming trend during the mid- to late Holocene is in opposition to the cooling seen in other proxy records. However, most of these existing proxy records are biased towards summer temperatures. We argue that the opposing trends are related to the seasonally different orbital forcing over this interval. Furthermore, our reconstructed trend as well as the recent maximum are consistent with the greenhouse gas forcing and climate model simulations, thus reconciling differing estimates of Arctic and northern high-latitude temperature evolution during the Holocene.

Benthic foraminiferal record from Santa Barbara basin

The oxygen minimum zone (OMZ) of the late Quaternary California margin experienced abrupt and dramatic changes in strength and depth in response to changes in intermediate water ventilation, ocean productivity, and climate at orbital through millennial time scales. Expansion and contraction of the OMZ is exhibited at high temporal resolution (107-126 year) by quantitative benthic foraminiferal assemblage changes in two piston cores forming a vertical profile in Santa Barbara Basin (569 m, basin floor; 481 m, near sill depth) to 34 and 24 ka, respectively. Variation in the OMZ is quantified by new benthic foraminiferal groupings and new dissolved oxygen index based on documented relations between species and water-mass oxygen concentrations. Foraminiferal-based paleoenvironmental assessments are integrated with principal component analysis, bioturbation, grain size, CaCO3, total organic carbon, and d13C to reconstruct basin oxygenation history. Fauna responded similarly between the two sites, although with somewhat different magnitude and taxonomic expression. During cool episodes (Younger Dryas and stadials), the water column was well oxygenated, most strongly near the end of the glacial episode (17-16 ka; Heinrich 1). In contrast, the OMZ was strong during warm episodes (Bølling/Allerød, interstadials, and Pre-Boreal). During the Bølling/Allerød, the OMZ shoaled to <360 m of contemporaneous sea level, its greatest vertical expansion of the last glacial cycle. Assemblages were then dominated by Bolivina tumida, reflecting high concentrations of dissolved methane in bottom waters. Short decadal intervals were so severely oxygen-depleted that no benthic foraminifera were present. The middle to late Holocene (6-0 ka) was less dysoxic than the early Holocene.

AMS radiocarbon dates on ten monospecific samples of planktonic foraminifer Globigerinoides ruber of sediment core GeoB8509-2 (Table 1)

Aeolian and fluvial sediment transport to the Atlantic Ocean offshore Mauritania were reconstructed based on grain-size distributions of the carbonate-free silt fraction of three marine sediment records of Cap Timiris Canyon to monitor the climatic evolution of present-day arid north-western Africa. During the late Pleistocene, predominantly coarse-grained particles, which are interpreted as windborne dust, characterise glacial dry climate conditions with a low sea level and extended sand seas that reach onto the exposed continental shelf off Mauritania. Subsequent particle fining and the abrupt decrease in terrigenous supply are attributed to humid climate conditions and dune stabilisation on the adjacent African continent with the onset of the Holocene humid period. Indications for an ancient drainage system, which was discharging fluvial mud offshore via Cap Timiris Canyon, are provided by the finest end member for early to mid Holocene times. However, in comparison to the Senegal and Niger River further south, the river system connecting Cap Timiris Canyon with the Mauritanian hinterland was starved during the late Holocene and is non-discharging under present-day arid climate conditions.

Percentages of sand-sized material and biogenic components in sediment core GIK12392-1

The paleo-oceanography of the southeastern North Atlantic Ocean during the last 150,000 yr has been studied using biogenous and terrigenous components of hemipelagic sediments sampled close to the northwest African continental margin. Variations of oxygen isotope ratios in shells of benthic calcareous foraminifers in two cores allow the assignment of absolute ages to these cores (in the best case at 1000 yr increments). The uncorrected bulk sedimentation rates of the longest core range from 3.4 to 7.6 cm/ 1000 yr during Interglacial conditions, and from 6.5 to 9.9 cm/1000 yr during Glacial conditions; all other cores have given results of the same order of magnitude, but with generally increasing values towards the continental edge. The distribution of sediment components allow us to make inferences about paleo-oceanographic changes in this region. Frequencies of biogenic components from benthic organisms, oxygen isotope ratios measured in benthic calcareous foraminiferal shells, the total carbonate contents of the sediment and distributions of biogenic components from planktonic organisms often fluctuate in concert. However, all fluctuations which can be attributed to changes of the bottom water masses (North Atlantic Deep Water) seem to precede by several thousand years those which can be linked to changes of the surface water mass distributions or to changes of the climate over the neighboring land masses. Late Quaternary planktonic foraminiferal assemblages in the cores from the northwest African continental margin can be defined satisfactorily in the way that distributions of assemblages found in sediment surface samples from the northeast Atlantic Ocean have been explained. The distributions of assemblages in the northwest African cores can also be used to estimate past sea surface temperatures and salinities. The downcore record of these estimates reveals two warm periods during the last 150,000 yr, the lower one corresponding to the oxygen isotope stage 5 e (equivalent to the Eemian proper in Europe), the upper one to the younger half of the Holocene. Winter surface water temperatures during oxygen isotope stages 6, 4, 3, and 2 are remarkably constant in most cores, while summer sea surface temperatures during stage 3 reach values comparable to those of the warm periods during the Late Holocene and Eemian. Estimated winter sea surface temperatures range from > 16 °C to < 11°C, the summer sea surface temperatures from > 22 °C to < 15 °C during the last 150,000 yr. Estimates of the winter sea surface salinities fluctuate between 36.6? and 35.5?, the higher values being restricted to the warm periods since the penultimate Glacial. Estimates for sea surface temperatures and salinities for two cores from the center of today's coastal upwelling region show less pronounced fluctuations than the record of the open ocean cores in the case of the station 12379 off Cape Barbas, more pronounced in the case of station 12328 off Cape Blanc. Seasonal differences between winter and summer sea surface temperatures derived from the estimated temperatures are today more pronounced in the boundary region of the ocean to the continent than further away from the continent. The differences are generally higher during warm climatic periods of the last 150,000 yr than during cooler ones. The abundance of terrigenous grains in the coarse fractions generally decreases with increasing distance from the continental edge, and also from south to north. The dominant portion of the terrigenous detritus is carried out into the ocean during the relatively cool climatic periods (stage 6, 4, later part of stage 3, stage 2 and oldest part of stage 1). The enhanced precision of dating combined with the stratigraphic resolution of these high deposition rate cores make it clear that the peaks of the terrigenous input off this part of the northwest African continental margin occur simultaneously with times of rapid sea level fluctuations resulting from large volume changes of the large Glacial ice sheets.

(Table S1) Age determination of sediment cores from the Labrador Sea

High resolution flow speed reconstructions of two core sites located on Gardar Drift in the northeast Atlantic Basin and Orphan Knoll in the northwest Atlantic Basin reveal a long-term decrease in flow speed of Northeast Atlantic Deep Water (NEADW) after 6,500 years. Benthic foraminiferal oxygen isotopes of sites currently bathed in NEADW show a 0.2per mil depletion after 6,500 years, shortly after the start of the development of a carbon isotope gradient between NEADW and Norwegian Sea Deep Water. We consider these changes in near-bottom flow vigor and benthic foraminiferal isotope records to mark a significant reorganization of the Holocene deep ocean circulation, and attribute the changes to a weakening of NEADW flow during the mid to late Holocene that allowed the shoaling of Lower Deep Water and deeper eastward advection of Labrador Sea Water into the northeast Atlantic Basin.

(Table S1) TEX86 and BIT record from sediment core PBBC-1

A molecular organic geochemical proxy (TEX86) for sea surface temperature (SST) is compared with a foraminifera-based SST proxy (Mg/Ca) in a decadal-resolution marine sedimentary record spanning the last 1000 years from the Gulf of Mexico. We assess the relative strengths of the organic and inorganic paleoceanographic techniques for reconstructing high-resolution SST variability during recent climate events, including the Little Ice Age (LIA) and the Medieval Warm Period (MWP). SST estimates based on the molecular organic proxy TEX86 show a similar magnitude and pattern of SST variability to foraminiferal Mg/Ca-SST estimates but with some important differences. For instance, both proxies show a cooling (1°C-2°C) of Gulf of Mexico SSTs during the LIA. During the MWP, however, Mg/Ca-SSTs are similar to near-modern SSTs, while TEX86 indicates SSTs that were cooler than modern. Using the respective SST calibrations for each proxy results in TEX86-SST estimates that are 2°C-4°C warmer than Mg/Ca-SST throughout the 1000 year record. We interpret the TEX86-SST as a summer-weighted SST signal from the upper mixed layer, whereas the Mg/Ca-SST better reflects the mean annual SST. Downcore differences in the SST estimates between the two proxies (DeltaT = TEX86 - Mg/Ca) are interpreted in the context of varying seasonality and/or changing water column temperature gradients.

(Table 1) Mean annual ground temperatures and description of permafrost boreholes in Russia

The results of the International Permafrost Association's International Polar Year Thermal State of Permafrost (TSP) project are presented based on field measurements from Russia during the IPY years (2007-09) and collected historical data. Most ground temperatures measured in existing and new boreholes show a substantial warming during the last 20 to 30 years. The magnitude of the warming varied with location, but was typically from 0.5°C to 2°C at the depth of zero annual amplitude. Thawing of Little Ice Age permafrost is ongoing at many locations. There are some indications that the late Holocene permafrost has begun to thaw at some undisturbed locations in northeastern Europe and northwest Siberia. Thawing of permafrost is most noticeable within the discontinuous permafrost domain. However, permafrost in Russia is also starting to thaw at some limited locations in the continuous permafrost zone. As a result, a northward displacement of the boundary between continuous and discontinuous permafrost zones was observed. This data set will serve as a baseline against which to measure changes of near-surface permafrost temperatures and permafrost boundaries, to validate climate model scenarios, and for temperature reanalysis.