Data from the EPICA Dome C ice core EDC

The Antarctic Vostok ice core provided compelling evidence of the nature of climate, and of climate feedbacks, over the past 420,000 years. Marine records suggest that the amplitude of climate variability was smaller before that time, but such records are often poorly resolved. Moreover, it is not possible to infer the abundance of greenhouse gases in the atmosphere from marine records. Here we report the recovery of a deep ice core from Dome C, Antarctica, that provides a climate record for the past 740,000 years. For the four most recent glacial cycles, the data agree well with the record from Vostok. The earlier period, between 740,000 and 430,000 years ago, was characterized by less pronounced warmth in interglacial periods in Antarctica, but a higher proportion of each cycle was spent in the warm mode. The transition from glacial to interglacial conditions about 430,000 years ago (Termination V) resembles the transition into the present interglacial period in terms of the magnitude of change in temperatures and greenhouse gases, but there are significant differences in the patterns of change. The interglacial stage following Termination V was exceptionally long - 28,000 years compared to, for example, the 12,000 years recorded so far in the present interglacial period. Given the similarities between this earlier warm period and today, our results may imply that without human intervention, a climate similar to the present one would extend well into the future.

(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.

Species occurrence and richness of lichen, bryophytes and vascular plants near Abisko, Sweden and Toolik lake, Alaska

Little is known about the impact of changing temperature regimes on composition and diversity of cryptogam communities in the Arctic and Subarctic, despite the well-known importance of lichens and bryophytes to the functioning and climate feedbacks of northern ecosystems. We investigated changes in diversity and abundance of lichens and bryophytes within long-term (9-16 years) warming experiments and along natural climatic gradients, ranging from Swedish subarctic birch forest and subarctic/subalpine tundra to Alaskan arctic tussock tundra. In both Sweden and Alaska, lichen diversity responded negatively to experimental warming (with the exception of a birch forest) and to higher temperatures along climatic gradients. Bryophytes were less sensitive to experimental warming than lichens, but depending on the length of the gradient, bryophyte diversity decreased both with increasing temperatures and at extremely low temperatures. Among bryophytes, Sphagnum mosses were particularly resistant to experimental warming in terms of both abundance and diversity. Temperature, on both continents, was the main driver of species composition within experiments and along gradients, with the exception of the Swedish subarctic birch forest where amount of litter constituted the best explanatory variable. In a warming experiment in moist acidic tussock tundra in Alaska, temperature together with soil ammonium availability were the most important factors influencing species composition. Overall, dwarf shrub abundance (deciduous and evergreen) was positively related to warming but so were the bryophytes Sphagnum girgensohnii, Hylocomium splendens and Pleurozium schreberi; the majority of other cryptogams showed a negative relationship to warming. This unique combination of intercontinental comparison, natural gradient studies and experimental studies shows that cryptogam diversity and abundance, especially within lichens, is likely to decrease under arctic climate warming. Given the many ecosystem processes affected by cryptogams in high latitudes (e.g. carbon sequestration, N2-fixation, trophic interactions), these changes will have important feedback consequences for ecosystem functions and climate.

Stable isotope ratios of foraminifera from various DSDP/ODP sites

The equator to high southern latitude sea surface and vertical temperature gradients are reconstructed from oxygen isotope values of planktonic and benthic foraminifers for the following five time intervals: late Paleocene, early Eocene, early middle Eocene, late Eocene, and early Oligocene. Paleotemperatures are calculated using standard oxygen isotope/temperature equations with adjustments to account for (1) variations in sea water delta18O related to changes in global ice volume over time and (2) latitudinal gradients in surface water delta18O. These reconstructions indicate that sea-surface temperatures (SST) of the Southern Oceans in the early Eocene were as high as 15°C, whereas temperatures during the late Paleocene and early middle Eocene reached maximum levels of 10°-12°C. By the late Eocene and early Oligocene high latitude SST had declined to 6 and 4°C, respectively. For most of the early Paleogene, low latitude sub-tropical temperatures remained constant and well within the range of Holocene temperatures (24°-25°C) but by the late Eocene and early Oligocene declined to values in the range of 18° to 22°C. The late Paleogene apparent decline in tropical temperatures, however, might be artificial because of dissolution of near-surface foraminifera tests which biased sediment assemblages toward deeper-dwelling foraminifera. Moreover, according to recent plate reconstructions, it appears that the majority of sites upon which the late Eocene and early Oligocene tropical temperatures were previously established were located either in or near regions likely to have been influenced by upwelling. Global deepwater temperature on average paralleled southern ocean SST for most of the Paleogene. We speculate based on the overall timing and character of marine sea surface temperature variation during the Paleogene that some combination of both higher levels of greenhouse gases and increased heat transport was responsible for the exceptional high-latitude warmth of the early Eocene.

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.

(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.

Sea surface temperature reconstruction based on alkenones of sediment core M35003-4

Evidence for abrupt climate changes on millennial and shorter timescales is widespread in marine and terrestrial climate records (Dansgard et al., 1993, doi:10.1038/364218a0; Bond et al., 1993, doi:10.1038/365143a0; Charles et al., 1996, doi:10.1016/0012-821X(96)00083-0, Bard et al., 1997, doi:10.1038/385707a0). Rapid reorganization of ocean circulation is considered to exert some control over these changes (Broecker et al., 1985, doi:10.1038/315021a0), as are shifts in the concentrations of atmospheric greenhouse gases (Broecker, 1994, doi:10.1038/372421a0). The response of the climate system to these two influences is fundamentally different: slowing of thermohaline overturn in the North Atlantic Ocean is expected to decrease northward heat transport by the ocean and to induce warming of the tropical Atlantic (Crowley, 1992, doi:10.1029/92PA01058; Manabe and Stouffer, 1997, doi:10.1029/96PA03932), whereas atmospheric greenhouse forcing should cause roughly synchronous global temperature changes (Manabe et al., 1991, doi:10.1175/1520-0442(1991)004<0785:TROACO>2.0.CO;2). So these two mechanisms of climate change should be distinguishable by the timing of surface-water temperature variations relative to changes in deep-water circulation. Here we present a high-temporal-resolution record of sea surface temperatures from the western tropical North Atlantic Ocean which spans the past 29,000 years, derived from measurements of temperature-sensitive alkenone unsaturation in sedimentary organic matter. We find significant warming is documented for Heinrich event H1 (16,900-15,400 calendar years bp) and the Younger Dryas event (12,900-11,600 cal. yr bp), which were periods of intense cooling in the northern North Atlantic. Temperature changes in the tropical and high-latitude North Atlantic are out of phase, suggesting that the thermohaline circulation was the important trigger for these rapid climate changes.

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.

Seawater carbonate chemistry and particulate organic particles during a semicontinuous batch culture experiment with Trichodesmium IMS101, 2007

Diazotrophic (N2-fixing) cyanobacteria provide the biological source of new nitrogen for large parts of the ocean. However, little is known about their sensitivity to global change. Here we show that the single most important nitrogen fixer in today's ocean, Trichodesmium, is strongly affected by changes in CO2 concentrations. Cell division rate doubled with rising CO2 (glacial to projected year 2100 levels) prompting lower carbon, nitrogen and phosphorus cellular contents, and reduced cell dimensions. N2 fixation rates per unit of phosphorus utilization as well as C:P and N:P ratios more than doubled at high CO2, with no change in C:N ratios. This could enhance the productivity of N-limited oligotrophic oceans, drive some of these areas into P limitation, and increase biological carbon sequestration in the ocean. The observed CO2 sensitivity of Trichodesmium could thereby provide a strong negative feedback to atmospheric CO2 increase.

The Lena River Delta - Land cover classification of tundra environments based on Landsat 7 ETM+ data and its application for upscaling of methane emissions

The Lena River Delta, situated in Northern Siberia (72.0 - 73.8° N, 122.0 - 129.5° E), is the largest Arctic delta and covers 29,000 km**2. Since natural deltas are characterised by complex geomorphological patterns and various types of ecosystems, high spatial resolution information on the distribution and extent of the delta environments is necessary for a spatial assessment and accurate quantification of biogeochemical processes as drivers for the emission of greenhouse gases from tundra soils. In this study, the first land cover classification for the entire Lena Delta based on Landsat 7 Enhanced Thematic Mapper (ETM+) images was conducted and used for the quantification of methane emissions from the delta ecosystems on the regional scale. The applied supervised minimum distance classification was very effective with the few ancillary data that were available for training site selection. Nine land cover classes of aquatic and terrestrial ecosystems in the wetland dominated (72%) Lena Delta could be defined by this classification approach. The mean daily methane emission of the entire Lena Delta was calculated with 10.35 mg CH4/m**2/d. Taking our multi-scale approach into account we find that the methane source strength of certain tundra wetland types is lower than calculated previously on coarser scales.

Time series of air chemistry measurements at Georg-von-Neumayer station, Dronning Maud Land, Antarctica in the years 1982 to 1991

The first Air Chemistry Observatory at the German Antarctic station Georg von Neumayer (GvN) was operated for 10 years from 1982 to 1991. The focus of the established observational programme was on characterizing the physical properties and chemical composition of the aerosol, as well as on monitoring the changing trace gas composition of the background atmosphere, especially concerning greenhouse gases. The observatory was designed by the Institut für Umweltphysik, University of Heidelberg (UHEIIUP). The experiments were installed inside the bivouac lodge, mounted on a sledge and put upon a snow hill to prevent snow accumulation during blizzards. All experiments were under daily control and daily performance protocols were documented. A ventilated stainless steel inlet stack (total height about 3-4 m above the snow surface) with a 50% aerodynamic cut-off diameter around 7-10 µm at wind velocities between 4-10 m/s supplied all experiments with ambient air. Contamination free sampling was realized by several means: (i) The Air Chemistry Observatory was situated in a clean air area about 1500 m south of GvN. Due to the fact that northern wind directions are very rare, contamination from the base can be excluded for most of the time. (ii) The power supply (20 kW) is provided by a cable from the main station, thus no fuel-driven generator is operated in the very vicinity. (iii) Contamination-free sampling is controlled by the permanently recorded wind velocity, wind direction and by condensation particle concentration. Contamination was indicated if one of the following criteria were given: Wind direction within a 330°-30° sector, wind velocity <2.2 m/s or >17.5 m/s, or condensation particle concentrations >2500/cm**3 during summer, >800/cm**3 during spring/autumn and >400/cm**3 during winter. If one or a definable combination of these criteria were given, high volume aerosol sampling and part of the trace gas sampling were interrupted. Starting at 1982 through 1991-01-14 surface ozone was measured with an electrochemical concentration cell (ECC). Surface ozone mixing ratio are given in ppbv = parts per 10**9 by volume. The averaging time corresponds to the given time intervals in the data sheet. The accuracy of the values are better than ±1 ppbv and the detection limit is around 1.0 ppbv. Aerosols were sampled on two Whatman 541 cellulose filters in series and analyzed by ion chromatography at the UHEI-IUP. Generally, the sampling period was seven days but could be up to two weeks on occasion. The air flow was around 100 m**3/h and typically 10000-20000 m**3 of ambient air was forced through the filters for one sample. Concentration values are given in nanogram (ng) per 1 m**3 air at standard pressure and temperature (1013 mbar, 273.16 K). Uncertainties of the values were approximately ±10% to ±15% for the main components MSA, chloride, nitrate, sulfate and sodium, and between ±20% and ±30% for the minor species bromide, ammonium, potassium, magnesium and calcium.

Barite Accumulation rates across Paleocene-Eocene Thermal Maximum (PETM)

Barite accumulation rates (BAR) have been measured from 12 DSDP/ODP site globally (DSDP site 525, 549 and ODP site 690, 738, 1051, 1209, 1215, 1220, 1221, 1263,1265 and 1266A) to reconstruct the export production across Paleocene Eocene Thermal Maximum (PETM) around 55.9 million year ago. Our results suggesting a general increase in export productivity. We propose that changes in marine ecosystems, resulting from high atmospheric partial pressure of CO2 and ocean acidification, led to enhanced carbon export from the photic zone to depth, thereby increasing the efficiency of the biological pump. We estimate that an annual carbon export flux out of the euphotic zone and into the deep ocean waters could have amounted to about 15 Gt during the PETM. About 0.4% of this carbon is expected to have entered the refractory dissolved organic pool, where it could be sequestered from the atmosphere for tens of thousands of years. Our estimates are consistent with the amount of carbon redistribution expected for the recovery from the PETM.