Data on mineral composition of volcanic glass and radiocarbon age of Upper Quaternary sediments from the Sea of Okhotsk

Lithological horizons have been distinguished in sediments cores from different parts of the Sea of Okhotsk based on primary descriptions of sediments and smear slides, and analyses of contents of both calcium carbonate and organic carbon, and opal. Sediment lithology has been correlated with oxygen isotope records and the standard isotope scale and radiocarbon data by AMS method for three cores studied in detail. This allowed to determine in detail periods of carbonaceous and diatomaceous ooze accumulation in the Sea of Okhotsk. Changes in magnetic susceptibility and grain size composition of sediments have been also compared with oxygen-isotope curves and radiocarbon datings. Obtained results confirm that variations in magnetic susceptibility are related with oxygen-isotope stages and influenced by climatic changes. Tephra interlayers K0, TR, K2, K3 have been identified by mineralogical analyses in all studied cores. Stratigraphic location of these tephra interlayers in detailed studied cores and their radiocarbon ages (8.1, 8.05, 26.8, and about 60 ka, respectively) provided base correlation between the interlayers and volcanic eruptions on the Kamchatka Peninsula and the Kuril Islands. This allows to use the former ones as time markers for deep-sea sediments of the Sea of Okhotsk. New lithostratigraphic and tephrochronologic data obtained allowed to correlate Upper Quaternary sediments from the Sea of Okhotsk.

O, C and Sr isotopic data for plankton and benthic foraminifera from DSDP Legs 26, 32, 36, and 71 sediments

A detailed d18O and d13C stratigraphy has been generated from analysis of well-preserved Albian - Early Maastrichtian foraminifera from Deep Sea Drilling Project (DSDP) Sites 511 and 327 (Falkland Plateau; ~58°S - 62°S paleolatitude) in the southern South Atlantic, and Cenomanian and Coniacian - Santonian foraminifera from DSDP Site 258 (Naturaliste Plateau; ~58°S paleolatitude) in the southern Indian Ocean. These results, when combined with previously published Maastrichtian stable isotope data from Ocean Drilling Program (ODP) Site 690 (Weddell Sea, ~65°S paleolatitude), provide new insight into the climatic and oceanographic history of the southern high latitudes during Middle-Late Cretaceous time. The planktonic foraminifer d18O curves reveal a gradual warming of surface waters from the Albian through the Cenomanian followed by extremely warm surface waters from the Turonian through the early Campanian. Long-term cooling of surface waters began in the late early Campanian and continued through the end of the Maastrichtian. The benthic foraminifer d18O record generally parallels changes in the oxygen isotopic curves defined by shallow-dwelling planktonic foraminifera. The vertical oxygen and carbon isotopic gradients were relatively low during the Albian - Cenomanian, high from the Turonian - Early Campanian, and then low during the late Campanian and Maastrichtian.

Results of Heinrich event 1 simulation from the University of Victoria Earth System-Climate Model (UVic ESCM) and the Community Climate System Model (CCSM3)

We investigated changes in tropical climate and vegetation cover associated with abrupt climate change during Heinrich Event 1 (HE1, ca. 17.5 ka BP) using two different global climate models: the University of Victoria Earth System-Climate Model (UVic ESCM) and the Community Climate System Model version 3 (CCSM3). Tropical South American and African pollen records suggest that the cooling of the North Atlantic Ocean during HE1 influenced the tropics through a southward shift of the rain belt. In this study, we simulated the HE1 by applying a freshwater perturbation to the North Atlantic Ocean. The resulting slowdown of the Atlantic Meridional Overturning Circulation was followed by a temperature seesaw between the Northern and Southern Hemispheres, as well as a southward shift of the tropical rain belt. The shift and the response pattern of the tropical vegetation around the Atlantic Ocean were more pronounced in the CCSM3 than in the UVic ESCM simulation. For tropical South America, opposite changes in tree and grass cover were modeled around 10° S in the CCSM3 but not in the UVic ESCM. In tropical Africa, the grass cover increased and the tree cover decreased around 15° N in the UVic ESCM and around 10° N in the CCSM3. In the CCSM3 model, the tree and grass cover in tropical Southeast Asia responded to the abrupt climate change during the HE1, which could not be found in the UVic ESCM. The biome distributions derived from both models corroborate findings from pollen records in southwestern and equatorial western Africa as well as northeastern Brazil.

ELPA (European Leaf Physiognomic Approach): Grid data set of the leaf physiognomic composition of the extant European hardwood vegetation

Physiognomic traits of plant leaves such as size, shape or margin are decisively affected by the prevailing environmental conditions of the plant habitat. On the other hand, if a relationship between environment and leaf physiognomy can be shown to exist, vegetation represents a proxy for environmental conditions. This study investigates the relationship between physiognomic traits of leaves from European hardwood vegetation and environmental parameters in order to create a calibration dataset based on high resolution grid cell data. The leaf data are obtained from synthetic chorologic floras, the environmental data comprise climatic and ecologic data. The high resolution of the data allows for a detailed analysis of the spatial dependencies between the investigated parameters. The comparison of environmental parameters and leaf physiognomic characters reveals a clear correlation between temperature related parameters (e.g. mean annual temperature or ground frost frequency) and the expression of leaf characters (e.g. the type of leaf margin or the base of the lamina). Precipitation related parameters (e.g. mean annual precipitation), however, show no correlation with the leaf physiognomic composition of the vegetation. On the basis of these results, transfer functions for several environmental parameters are calculated from the leaf physiognomic composition of the extant vegetation. In a next step, a cluster analysis is applied to the dataset in order to identify "leaf physiognomic communities". Several of these are distinguished, characterised and subsequently used for vegetation classification. Concerning the leaf physiognomic diversity there are precise differences between each of these "leaf physiognomic classes". There is a clear increase of leaf physiognomic diversity with increasing variability of the environmental parameters: Northern vegetation types are characterised by a more or less homogeneous leaf physiognomic composition whereas southern vegetation types like the Mediterranean vegetation show a considerable higher leaf physiognomic diversity. Finally, the transfer functions are used to estimate palaeo-environmental parameters of three fossil European leaf assemblages from Late Oligocene and Middle Miocene. The results are compared with results obtained from other palaeo-environmental reconstructing methods. The estimates based on a direct linear ordination seem to be the most realistic ones, as they are highly consistent with the Coexistence Approach.

Treeline dynamics in Siberia under changing climates as inferred from an individual-based model for Larix, supplementary material

Siberian boreal forests are expected to expand northwards in the course of global warming. However, processes of the treeline ecotone transition, as well astiming and related climate feedbacks are still not understood. Here, we present 'Larix Vegetation Simulator' LAVESI, an individual-based spatially-explicit model that can simulate Larix gmelinii (RUPR.) RUPR. stand dynamics in an attempt to improve our understanding about past and future treeline movements under changing climates. The relevant processes (growth, seed production and dispersal, establishment and mortality) are incorporated and adjusted to observation data mainly gained from the literature. Results of a local sensitivity analysis support the robustness of the model's parameterization by giving relatively small sensitivity values. We tested the model by simulating tree stands under modern climate across the whole Taymyr Peninsula, north-central Siberia (c. 64-80° N; 92-119° E). We find tree densities similar to observed forests in the northern to mid-treeline areas, but densities are overestimated in the southern parts of the simulated region. Finally, from a temperature-forcing experiment, we detect that the responses of tree stands lag the hypothetical warming by several decades, until the end of 21st century. With our simulation experiments we demonstrate that the newly-developed model captures the dynamics of the Siberian latitudinal treeline.

Pollen data and pollen-based climate reconstruction for the KTK1 pollen record from Lake Kotokel

In this study a radiocarbon-dated pollen record from Lake Kotokel (52°47' N, 108°07' E, 458 m a.s.l.) located in southern Siberia east of Lake Baikal was used to derive quantitative characteristics of regional vegetation and climate from about 15 kyr BP (1 kyr = 1000 cal. yr) until today. Quantitative reconstruction of the late glacial vegetation and climate dynamics suggests that open steppe and tundra communities predominated in the study area prior to ca. 13.5 kyr BP and again during the Younger Dryas interval, between 12.8 and 11.6 kyr BP. The pollen-based climate reconstruction suggests lower-than-present mean January (~ -38 °C) and July (~ 12 °C) temperatures and annual precipitation (~ 270-300 mm) values during these time intervals. Boreal woodland replaced the primarily open landscape around Kotokel three times at about 14.8-14.7 kyr BP, during the Allerød Interstadial between 13.3-12.8 kyr BP and with the onset of the Holocene interglacial between 11.5 and 10.5 kyr BP, presumably in response to a noticeable increase in precipitation, and in July and January temperatures. The maximal spread of the boreal forest (taiga) communities in the region is associated with a warmer and wetter-than-present climate (Tw ~ 17-18 °C, Tc ~ -19 °C, Pann ~ 500-550 mm) that occurred ca. 10.8-7.3 kyr BP. During this time interval woody vegetation covered more than 50 % of the area within a 21x21 km window around the lake. The pollen-based best modern analogue reconstruction suggests a decrease in woody cover percentages and in all climatic variables about 7-6.5 kyr BP. Our results demonstrate a gradual decrease in precipitation and mean January temperature towards their present-day values in the region around Lake Kotokel since that time.

ROCADA: Romanian daily gridded climatic dataset (1961-2013) V1.0

Daily records of nine meteorological variables covering the interval 1961-2013 were used in order to create a state-of-the-art homogenized climatic dataset over Romania at a spatial resolution of 0.1°. All meteorological stations with full data records, as well as stations with up to 30 % missing data, were used for the following variables: air pressure (150 stations); minimum, maximum, and average air temperature (150 stations); soil temperature (127 stations); precipitation (188 stations); sunshine hours (135 stations); cloud cover (104 stations); relative humidity (150 stations). For each parameter, the data series were first homogenized with the software MASH (Multiple Analysis of Series for Homogenization); then, the data series were gridded by means of the software MISH (Meteorological Interpolation based on Surface Homogenized Data).

Grain-size, lithic grains, foraminifera-derived and dinocyst-derived data of sediment core MD99-2281

The last glacial period was punctuated by abrupt climatic events with extrema known as Heinrich and Dansgaard-Oeschger events. These millennial events have been the subject of many paleoreconstructions and model experiments in the past decades, but yet the hydrological processes involved remain elusive. In the present work, high-resolution analyses were conducted on the 12-42 ka BP section of core MD99-2281 retrieved southwest of the Faeroe Islands, and combined with analyses conducted in two previous studies (Zumaque et al., 2012; Caulle et al., 2013). Such a multiproxy approach, coupling micropaleontological, geochemical and sedimentological analyses, allows us to track surface, subsurface, and deep hydrological processes occurring during these rapid climatic changes. Records indicate that the coldest episodes of the studied period (Greenland stadials and Heinrich stadials) were characterized by a strong stratification of surface waters. This surface stratification seems to have played a key role in the dynamics of subsurface and deep-water masses. Indeed, periods of high surface stratification are marked by a coupling of subsurface and deep circulations which sharply weaken at the beginning of stadials, while surface conditions progressively deteriorate throughout these cold episodes; conversely, periods of decreasing surface stratification (Greenland interstadials) are characterized by a coupling of surface and deep hydrological processes, with progressively milder surface conditions and gradual intensification of the deep circulation, while the vigor of the subsurface northward Atlantic flow remains constantly high. Our results also reveal different and atypical hydrological signatures during Heinrich stadials (HSs): while HS1 and HS4 exhibit a "usual" scheme with reduced overturning circulation, a relatively active North Atlantic circulation seems to have prevailed during HS2, and HS3 seems to have experienced a re-intensification of this circulation during the middle of the event. Our findings thus bring valuable information to better understand hydrological processes occurring in a key area during the abrupt climatic shifts of the last glacial period.

Inventory of SINOPS project data sets

During the SINOPS project, an optimal state of the art simulation of the marine silicon cycle is attempted employing a biogeochemical ocean general circulation model (BOGCM) through three particular time steps relevant for global (paleo-) climate. In order to tune the model optimally, results of the simulations are compared to a comprehensive data set of 'real' observations. SINOPS' scientific data management ensures that data structure becomes homogeneous throughout the project. Practical work routine comprises systematic progress from data acquisition, through preparation, processing, quality check and archiving, up to the presentation of data to the scientific community. Meta-information and analytical data are mapped by an n-dimensional catalogue in order to itemize the analytical value and to serve as an unambiguous identifier. In practice, data management is carried out by means of the online-accessible information system PANGAEA, which offers a tool set comprising a data warehouse, Graphical Information System (GIS), 2-D plot, cross-section plot, etc. and whose multidimensional data model promotes scientific data mining. Besides scientific and technical aspects, this alliance between scientific project team and data management crew serves to integrate the participants and allows them to gain mutual respect and appreciation.