(Table 2) Survivorship of transplanted seedlings of Picea glauca inside and outside three tree islands in 1994, 1995, 2009 and after 15 years

The northern boundary of boreal forest and the ranges of tree species are expected to shift northward in response to climate warming, which will result in a decrease in the albedo of areas currently covered by tundra vegetation, an increase in terrestrial carbon sequestration, and an alteration of biodiversity in the current Low Arctic. Central to the prediction of forest expansion is an increase in the reproductive capacity and establishment of individual trees. We assessed cone production, seed viability, and transplanted seedling success of Picea glauca (Moench.) Voss. (white spruce) in the early 1990s and again in the late 2000s at four forest stand sites and eight tree island sites (clonal populations beyond present treeline) in the Mackenzie Delta region of the Northwest Territories, Canada. Over the past 20 years, average temperatures in this region have increased by 0.9 °C. This area has the northernmost forest-tundra ecotone in North America and is one of the few circumpolar regions where the northern limit of conifer trees reaches the Arctic Ocean. We found that cone production and seed viability did not change between the two periods of examination and that both variables decreased northward across the forest-tundra ecotone. Nevertheless, white spruce individuals at the northern limit of the forest-tundra ecotone produced viable seeds. Furthermore, transplanted seedlings were able to survive in the northernmost sites for 15 years, but there were no signs of natural regeneration. These results indicate that if climatic conditions continue to ameliorate, reproductive output will likely increase, but seedling establishment and forest expansion within the forest-tundra of this region is unlikely to occur without the availability of suitable recruitment sites. Processes that affect the availability of recruitment sites are likely to be important elsewhere in the circumpolar ecotone, and should be incorporated into models and predictions of climate change and its effects on the northern forest-tundra ecotone.

Geochemistry on Holocene sediments from the Mediterranean

The enhanced accumulation of organic matter in Eastern Mediterranean sapropels and their unusually low d15N values have been attributed to either enhanced nutrient availability which led to elevated primary production and carbon sequestration or to enhanced organic matter preservation under anoxic conditions. In order to evaluate these two hypothesis we have determined Ba/Al ratios, amino acid composition, N and organic C concentrations and d15N in sinking particles, surface sediments, eight spatially distributed core records of the youngest sapropel S1 (10-6 ka) and older sapropels (S5, S6) from two locations. These data suggest that (i) temporal and spatial variations in d15N of sedimentary N are driven by different degrees of diagenesis at different sites rather than by changes in N-sources or primary productivity and (ii) present day TOC export production would suffice to create a sapropel like S1 under conditions of deep-water anoxia. This implies that both enhanced TOC accumulation and d15N depletion in sapropels were due to the absence of oxygen in deep waters. Thus preservation plays a major role for the accumulation of organic-rich sediments casting doubt on the need of enhanced primary production for sapropel formation.

Age-depth relation and red colour reflectance of ODP Leg 177 sites

The Antarctic Polar Front is an important biogeochemical divider in the Southern Ocean. Laminated diatom mat deposits record episodes of massive flux of the diatom Thalassiothrix antarctica beneath the Antarctic Polar Front and provide a marker for tracking the migration of the Front through time. Ocean Drilling Program Sites 1091, 1093 and 1094 are the only deep piston cored record hitherto sampled from the sediments of the circumpolar biogenic opal belt. Mapping of diatom mat deposits between these sites indicates a glacial-interglacial front migration of up to 6 degrees of latitude in the early/mid Pleistocene. The mid-Pleistocene transition marks a stepwise minimum 7° northward migration of the locus of the Polar Front sustained for about 450 kyr until an abrupt southward return to a locus similar to its modern position and further south than any mid-Pleistocene locus. This interval from a "900 ka event" that saw major cooling of the oceans and a d13C minimum through to the 424 ka Mid-Brunhes Event at Termination V is also seemingly characterised by 1) sustained decreased carbonate in the sub-tropical south Atlantic, 2) reduced strength of Antarctic deep meridional circulation, 3) lower interglacial temperatures and lower interglacial atmospheric CO2 levels (by some 30 per mil) than those of the last 400 kyr, evidencing less complete deglaciation. This evidence is consistent with a prolonged period lasting 450 kyr of only partial ventilation of the deep ocean during interglacials and suggests that the mechanisms highlighted by recent hypotheses linking mid-latitude atmospheric conditions to the extent of deep ocean ventilation and carbon sequestration over glacial-interglacial cycles are likely in operation during the longer time scale characteristic of the mid-Pleistocene transition. The cooling that initiated the "900 ka event" may have been driven by minima in insolation amplitude related to eccentricity modulation of precession that also affected low latitude climates as marked by threshold changes in the African monsoon system. The major thresholds in earth system behaviour through the mid-Pleistocene transition were likely governed by an interplay of the 100 kyr and 400 kyr eccentricity modulation of precession.

Total nitrogen from solid phase in the Jena Experiment (Main Experiment up to 30cm depth, year 2008)

This data set contains measurements of total nitrogen from the main experiment plots of a large grassland biodiversity experiment (the Jena Experiment; see further details below). In the main experiment, 82 grassland plots of 20 x 20 m were established from a pool of 60 species belonging to four functional groups (grasses, legumes, tall and small herbs). In May 2002, varying numbers of plant species from this species pool were sown into the plots to create a gradient of plant species richness (1, 2, 4, 8, 16 and 60 species) and functional richness (1, 2, 3, 4 functional groups). Plots were maintained by bi-annual weeding and mowing. Soil sampling and analysis: Stratified soil sampling was performed in April 2008 to a depth of 30 cm. Three independent samples per plot were taken using a split tube sampler with an inner diameter of 4.8 cm (Eijkelkamp Agrisearch Equipment, Giesbeek, the Netherlands). Soil samples were segmented to a depth resolution of 5 cm in the field, giving six depth subsamples per core, and made into composite samples per depth. Sampling locations were less than 30 cm apart from sampling locations in other years. Samples were dried at 40°C. All soil samples were passed through a sieve with a mesh size of 2 mm. Because of much higher proportions of roots in the soil, the samples were further sieved to 1 mm according to common root removal methods. No additional mineral particles were removed by this procedure. Total nitrogen concentration was analyzed on ball-milled subsamples (time 4 min, frequency 30 s-1) by an elemental analyzer at 1150°C (Elementaranalysator vario Max CN; Elementar Analysensysteme GmbH, Hanau, Germany).

Total nitrogen from solid phase in the Jena Experiment (Block 2 Main Experiment up to 100cm depth, year 2002)

This data set contains measurements of total nitrogen from the main experiment plots of a large grassland biodiversity experiment (the Jena Experiment; see further details below). In the main experiment, 82 grassland plots of 20 x 20 m were established from a pool of 60 species belonging to four functional groups (grasses, legumes, tall and small herbs). In May 2002, varying numbers of plant species from this species pool were sown into the plots to create a gradient of plant species richness (1, 2, 4, 8, 16 and 60 species) and functional richness (1, 2, 3, 4 functional groups). Plots were maintained by bi-annual weeding and mowing. Soil sampling and analysis: Stratified soil sampling to a depth of 1m was performed before sowing in April 2002. Three independent samples per plot were taken of all plots in block 2 using a motor-driven soil column cylinder (Cobra, Eijkelkamp, 8.3 cm in diameter). Soil samples were dried at 40°C and segmented to a depth resolution of 5 cm giving 20 depth subsamples per core. All samples were analyzed independently. All soil samples were passed through a sieve with a mesh size of 2 mm. Rarely present visible plant remains were removed using tweezers. Total nitrogen concentration was analyzed on ball-milled subsamples (time 4 min, frequency 30 s-1) by an elemental analyzer at 1150°C (Elementaranalysator vario Max CN; Elementar Analysensysteme GmbH, Hanau, Germany).

Total nitrogen from solid phase in the Jena Experiment (Main Experiment up to 30cm depth, year 2004)

This data set contains measurements of total nitrogen from the main experiment plots of a large grassland biodiversity experiment (the Jena Experiment; see further details below). In the main experiment, 82 grassland plots of 20 x 20 m were established from a pool of 60 species belonging to four functional groups (grasses, legumes, tall and small herbs). In May 2002, varying numbers of plant species from this species pool were sown into the plots to create a gradient of plant species richness (1, 2, 4, 8, 16 and 60 species) and functional richness (1, 2, 3, 4 functional groups). Plots were maintained by bi-annual weeding and mowing. Soil sampling and analysis: Stratified soil sampling was performed in April 2004 to a depth of 30 cm. Three independent samples per plot were taken using a split tube sampler with an inner diameter of 4.8 cm (Eijkelkamp Agrisearch Equipment, Giesbeek, the Netherlands). Soil samples were segmented to a depth resolution of 5 cm in the field, giving six depth subsamples per core, and made into composite samples per depth. Sampling locations were less than 30 cm apart from sampling locations in other years. Samples were dried at 40°C. All soil samples were passed through a sieve with a mesh size of 2 mm. Because of much higher proportions of roots in the soil, the samples were further sieved to 1 mm according to common root removal methods. No additional mineral particles were removed by this procedure. Total nitrogen concentration was analyzed on ball-milled subsamples (time 4 min, frequency 30 s-1) by an elemental analyzer at 1150°C (Elementaranalysator vario Max CN; Elementar Analysensysteme GmbH, Hanau, Germany).

(Tabls S1) Bulk geochemistry for sediment samples collected from the Cenomanian-Turonian boundary section at Pont d'Issole

Oceanic Anoxic Event 2 (OAE2), spanning the Cenomanian-Turonian boundary (CTB), represents one of the largest perturbations in the global carbon cycle in the last 100 Myr. The d13Ccarb, d13Corg, and d18O chemostratigraphy of a black shale-bearing CTB succession in the Vocontian Basin of France is described and correlated at high resolution to the European CTB reference section at Eastbourne, England, and to successions in Germany, the equatorial and midlatitude proto-North Atlantic, and the U.S. Western Interior Seaway (WIS). Delta13C (offset between d13Ccarb and d13Corg) is shown to be a good pCO2 proxy that is consistent with pCO2 records obtained using biomarker d13C data from Atlantic black shales and leaf stomata data from WIS sections. Boreal chalk d18O records show sea surface temperature (SST) changes that closely follow the Delta13C pCO2 proxy and confirm TEX86 results from deep ocean sites. Rising pCO2 and SST during the Late Cenomanian is attributed to volcanic degassing; pCO2 and SST maxima occurred at the onset of black shale deposition, followed by falling pCO2 and cooling due to carbon sequestration by marine organic productivity and preservation, and increased silicate weathering. A marked pCO2 minimum (~25% fall) occurred with a SST minimum (Plenus Cold Event) showing >4°C of cooling in ~40 kyr. Renewed increases in pCO2, SST, and d13C during latest Cenomanian black shale deposition suggest that a continuing volcanogenic CO2 flux overrode further drawdown effects. Maximum pCO2 and SST followed the end of OAE2, associated with a falling nutrient supply during the Early Turonian eustatic highstand.

Stable isotope record, opal accumulation rate and geochemistry of middle Miocene sediments of IODP Site 321-U1338

During the middle Miocene, Earth's climate transitioned from a relatively warm phase (Miocene climatic optimum) into a colder mode with re-establishment of permanent ice sheets on Antarctica, thus marking a fundamental step in Cenozoic cooling. Carbon sequestration and atmospheric CO2 drawdown through increased terrestrial and/or marine productivity have been proposed as the main drivers of this fundamental transition. We integrate high-resolution (1-3 k.y.) benthic stable isotope data with XRF-scanner derived biogenic silica and carbonate accumulation estimates in an exceptionally well-preserved sedimentary archive, recovered at Integrated Ocean Drilling Program Site U1338, to reconstruct eastern equatorial Pacific productivity variations and to investigate temporal linkages between high- and low-latitude climate change over the interval 16-13 Ma. Our records show that the climatic optimum (16.8-14.7 Ma) was characterized by high amplitude climate variations, marked by intense perturbations of the carbon cycle. Episodes of peak warmth at (southern hemisphere) insolation maxima coincided with transient shoaling of the carbonate compensation depth and enhanced carbonate dissolution in the deep ocean. A switch to obliquity-paced climate variability after 14.7 Ma concurred with a general improvement in carbonate preservation and the onset of stepwise global cooling, culminating with extensive ice growth over Antarctica at ~13.8 Ma. We find that two massive increases in opal accumulation at ~14.0 and ~13.8 Ma occurred just before and during the final and most prominent cooling step, supporting the hypothesis that enhanced siliceous productivity in the eastern equatorial Pacific contributed to CO2 drawdown.

Ba/Sr ratio in marine barite from five DSDP and ODP holes

The Paleocene-Eocene Thermal Maximum (PETM), ca. 55 Ma, was a period of extreme global warming caused by rapid emission of greenhouse gases. It is unknown what ended this episode of greenhouse warming, but high oceanic export productivity over thousands of years (as indicated by high accumulation rates of barium, Ba) may have been a factor in ending this warm period by carbon sequestration. However, Ba has a short oceanic residence time (~10 k.y.), so a prolonged global increase in Ba accumulation rates requires an increase in input of Ba to the ocean, increasing barite saturation. We use a novel proxy for barite saturation (Sr/Ba in marine barite) to demonstrate that the seawater saturation state with respect to barite did not change across the PETM. The observations of increased barite burial, no change in saturation, and the short residence time can be reconciled if Ba burial decreased at continental margin and shelf sites due to widespread occurrence of suboxic conditions, leading to Ba release into the water column, combined with increased biological export production at some pelagic sites, resulting in Ba sink reorganization.

(Table S1) Diatom abundance in coastal upwelling systems

Coastal upwelling systems account for approximately half of global ocean primary production and contribute disproportionately to biologically driven carbon sequestration. Diatoms, silica-precipitating microalgae, constitute the dominant phytoplankton in these productive regions, and their abundance and assemblage composition in the sedimentary record is considered one of the best proxies for primary production. The study of the sedimentary diatom abundance (SDA) and total organic carbon content (TOC) in the five most important coastal upwelling systems of the modern ocean (Iberia-Canary, Benguela, Peru-Humboldt, California and Somalia-Oman) reveals a global-scale positive relationship between diatom production and organic carbon burial. The analysis of SDA in conjunction with environmental variables of coastal upwelling systems such as upwelling strength, satellite-derived net primary production and surface water nutrient concentrations shows different relations between SDA and primary production on the regional scale. At the global-scale, SDA appears modulated by the capacity of diatoms to take up silicic acid, which ultimately sets an upper limit to global export production in these ocean regions.

Total nitrogen from solid phase in the Jena Experiment (Block 2 Main Experiment up to 100cm depth, year 2007)

This data set contains measurements of total nitrogen from the main experiment plots of a large grassland biodiversity experiment (the Jena Experiment; see further details below). In the main experiment, 82 grassland plots of 20 x 20 m were established from a pool of 60 species belonging to four functional groups (grasses, legumes, tall and small herbs). In May 2002, varying numbers of plant species from this species pool were sown into the plots to create a gradient of plant species richness (1, 2, 4, 8, 16 and 60 species) and functional richness (1, 2, 3, 4 functional groups). Plots were maintained by bi-annual weeding and mowing. Stratified soil sampling to a depth of 1m was repeated in April 2007 (as had been done before sowing in April 2002). Three independent samples per plot were taken of all plots in block 2 using a motor-driven soil column cylinder (Cobra, Eijkelkamp, 8.3 cm in diameter). Soil samples were dried at 40°C and segmented to a depth resolution of 5 cm giving 20 depth subsamples per core. All samples were analyzed independently. All soil samples were passed through a sieve with a mesh size of 2 mm. Because of much higher proportions of roots in the soil, the samples in 2007 were further sieved to 1 mm according to common root removal methods. No additional mineral particles were removed by this procedure. Total nitrogen concentration was analyzed on ball-milled subsamples (time 4 min, frequency 30 s-1) by an elemental analyzer at 1150°C (Elementaranalysator vario Max CN; Elementar Analysensysteme GmbH, Hanau, Germany).

Total nitrogen from solid phase in the Jena Experiment (Main Experiment up to 30cm depth, year 2006)

This data set contains measurements of total nitrogen from the main experiment plots of a large grassland biodiversity experiment (the Jena Experiment; see further details below). In the main experiment, 82 grassland plots of 20 x 20 m were established from a pool of 60 species belonging to four functional groups (grasses, legumes, tall and small herbs). In May 2002, varying numbers of plant species from this species pool were sown into the plots to create a gradient of plant species richness (1, 2, 4, 8, 16 and 60 species) and functional richness (1, 2, 3, 4 functional groups). Plots were maintained by bi-annual weeding and mowing. Soil sampling and analysis: Stratified soil sampling was performed in April 2006 to a depth of 30 cm. Three independent samples per plot were taken using a split tube sampler with an inner diameter of 4.8 cm (Eijkelkamp Agrisearch Equipment, Giesbeek, the Netherlands). Soil samples were segmented to a depth resolution of 5 cm in the field, giving six depth subsamples per core, and made into composite samples per depth. Sampling locations were less than 30 cm apart from sampling locations in other years. Samples were dried at 40°C. All soil samples were passed through a sieve with a mesh size of 2 mm. Because of much higher proportions of roots in the soil, the samples were further sieved to 1 mm according to common root removal methods. No additional mineral particles were removed by this procedure. Total nitrogen concentration was analyzed on ball-milled subsamples (time 4 min, frequency 30 s-1) by an elemental analyzer at 1150°C (Elementaranalysator vario Max CN; Elementar Analysensysteme GmbH, Hanau, Germany).

Total nitrogen from solid phase in the Jena Experiment (Main Experiment up to 30cm depth, year 2002)

This data set contains measurements of total nitrogen from the main experiment plots of a large grassland biodiversity experiment (the Jena Experiment; see further details below). In the main experiment, 82 grassland plots of 20 x 20 m were established from a pool of 60 species belonging to four functional groups (grasses, legumes, tall and small herbs). In May 2002, varying numbers of plant species from this species pool were sown into the plots to create a gradient of plant species richness (1, 2, 4, 8, 16 and 60 species) and functional richness (1, 2, 3, 4 functional groups). Plots were maintained by bi-annual weeding and mowing. Soil sampling and analysis: Stratified soil sampling was performed before sowing in April 2002. Five independent samples per plot were taken using a split tube sampler with an inner diameter of 4.8 cm (Eijkelkamp Agrisearch Equipment, Giesbeek, the Netherlands). Soil samples were dried at 40°C and then segmented to a depth resolution of 5 cm giving six depth subsamples per core. All samples were analyzed independently and averaged values per depth layer are reported. Sampling locations were less than 30 cm apart from sampling locations in other years. Subsequently, samples were dried at 40°C. All soil samples were passed through a sieve with a mesh size of 2 mm. Rarely present visible plant remains were removed using tweezers. Total nitrogen concentration was analyzed on ball-milled subsamples (time 4 min, frequency 30 s-1) by an elemental analyzer at 1150°C (Elementaranalysator vario Max CN; Elementar Analysensysteme GmbH, Hanau, Germany).

Biogeochemistry of the Eastern Weddell Gyre

The Southern Ocean (SO) plays a key role in modulating atmospheric CO2 via physical and biological processes. However, over much of the SO, biological activity is iron-limited. New in situ data from the Antarctic zone south of Africa in a region centered at -20°E - 25°E reveal a previously overlooked region of high primary production, comparable in size to the northwest African upwelling region. Here, sea ice together with enclosed icebergs is channeled by prevailing winds to the eastern boundary of the Weddell Gyre, where a sharp transition to warmer waters causes melting. This cumulative melting provides a steady source of iron, fuelling an intense phytoplankton bloom that is not fully captured by monthly satellite production estimates. These findings imply that future changes in sea-ice cover and dynamics could have a significant effect on carbon sequestration in the SO.

Eddy covariance and environmental data of boreal bogs plant species

In boreal bogs plant species are low in number, but they differ greatly in their growth forms and photosynthetic properties. We assessed how ecosystem carbon (C) sink dynamics were affected by seasonal variations in photosynthetic rate and leaf area of different species. Photosynthetic properties (light-response parameters), leaf area development and areal cover (abundance) of the species were used to quantify species-specific net and gross photosynthesis rates (PN and PG, respectively), which were summed to express ecosystem-level PN and PG. The ecosystem-level PG was compared with a gross primary production (GPP) estimate derived from eddy covariance measurements (EC). Species areal cover rather than differences in photosynthetic properties determined the species with the highest PG of both vascular plants and Sphagna. Species-specific contributions to the ecosystem PG varied over the growing season, which in turn determined the seasonal variation in ecosystem PG. The upscaled growing-season PG estimate, 230 g C/m**2, agreed well with the GPP estimated by the EC, 243 g C/m**2. Sphagna were superior to vascular plants in ecosystem-level PG throughout the growing season but had a lower PN. PN results indicated that areal cover of the species together with their differences in photosynthetic parameters shape the ecosystem-level C balance. Species with low areal cover but high photosynthetic efficiency appear to be potentially important for the ecosystem C sink. Results imply that functional diversity may increase the stability of C sink of boreal bogs.