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

Mixed layer heat and salinity budgets during the onset of the 2011 Atlantic cold tongue

The mixed layer (ML) temperature and salinity changes in the central tropical Atlantic have been studied by a dedicated experiment (Cold Tongue Experiment (CTE)) carried out from May to July 2011. The CTE was based on two successive research cruises, a glider swarm, and moored observations. The acquired in situ data sets together with satellite, reanalysis, and assimilation model data were used to evaluate box-averaged ML heat and salinity budgets for two subregions: (1) the western equatorial Atlantic cold tongue (ACT) (23°-10°W) and (2) the region north of the ACT. The strong ML heat loss in the ACT region during the CTE was found to be the result of the balance of warming due to net surface heat flux and cooling due to zonal advection and diapycnal mixing. The northern region was characterized by weak cooling and the dominant balance of net surface heat flux and zonal advection. A strong salinity increase occurred at the equator, 10°W, just before the CTE. During the CTE, ML salinity in the ACT region slightly increased. Largest contributions to the ML salinity budget were zonal advection and the net surface freshwater flux. While essential for the ML heat budget in the ACT region, diapycnal mixing played only a minor role for the ML salinity budget. In the region north of the ACT, the ML freshened at the beginning of the CTE due to precipitation, followed by a weak salinity increase. Zonal advection changed sign contributing to ML freshening at the beginning of the CTE and salinity increase afterward.

(Table 1) Sea-surface temperatures, intermediate-water temperatures, and latitudinal gradient at three time slices at the end of the Cretaceous

Climatic and oceanographic variations during the last 2 m.y. of the Maastrichtian inferred from high-resolution (10 k.y.) stable isotope analysis of the mid-latitude South Atlantic Deep Sea Drilling Project Site 525 reveal a major warm pulse followed by rapid cooling prior to the Cretaceous-Tertiary boundary. Between 66.85 and 65.52 Ma, cool but fluctuating temperatures average 9.9 and 15.4°C in intermediate and surface waters, respectively. This interval is followed by an abrupt short-term warming between 65.45 and 65.11 Ma, which increased temperatures by 2-3°C in intermediate waters, and decreased the vertical thermal gradient to an average of 2.7°C. This warm pulse may be linked to increased atmospheric pCO2, increased poleward heat transport, and the switch of an intermediate water source from high to low-middle latitudes. During the last 100 k.y. of the Maastrichtian, intermediate and surface temperatures decreased by an average of 2.1 and 1.4°C, respectively, compared to the maximum temperature between 65.32 and 65.24 Ma.

Daily ice production and polynya area in the Laptev Sea region during 1979-2009

Polynyas in the Laptev Sea are examined with respect to recurrence and interannual wintertime ice production.We use a polynya classification method based on passive microwave satellite data to derive daily polynya area from long-term sea-ice concentrations. This provides insight into the spatial and temporal variability of open-water and thin-ice regions on the Laptev Sea Shelf. Using thermal infrared satellite data to derive an empirical thin-ice distribution within the thickness range from 0 to 20 cm, we calculate daily average surface heat loss and the resulting wintertime ice formation within the Laptev Sea polynyas between 1979 and 2008 using reanalysis data supplied by the National Centers for Environmental Prediction, USA, as atmospheric forcing. Results indicate that previous studies significantly overestimate the contribution of polynyas to the ice production in the Laptev Sea. Average wintertime ice production in polynyas amounts to approximately 55 km3 ± 27% and is mostly determined by the polynya area, wind speed and associated large-scale circulation patterns. No trend in ice production could be detected in the period from 1979/80 to 2007/08.

Siliceous microplankton in surface sediments and the water column of the Red Sea

We studied the siliceous microplankton assemblages (mainly diatoms) from plankton tows (mesh size 20 µm) and surface sediment samples collected along a N-S transect in the northern Red Sea (28-21°N). In addition, we analyzed differences/similarities between plankton and sediment assemblages within a brine-filled basin (the southern basin) of the Shaban Deep and compared these assemblages with those from outside the brine. Plankton samples revealed the overwhelming dominance of diatoms over other siliceous groups. Diatoms accounted for ca. 97% of all biosiliceous particles at 120-20 m (vs. 2.9% silicoflagellates and 0.4% radiolarians), and ca. 94% at 200-120 m (vs. 4.5% silicoflagellates and 1.6% radiolarians). In general, a marine, warm-water (tropical/subtropical) diatom assemblage characterizes the plankton samples. Representatives of the Nitzschia bicapitata group are by far the most abundant contributors at both depth intervals (average=43%), ranging from ca. 30% in the North to ca. 60% in the South. Biogenic opal content in non-brine surface sediments is very low, (below 0.2 wt.% SiO2); and concentration of siliceous microorganisms is also low and of the order of 5*10**3-10**4 microorganisms/g dry sediment. Diatoms are the main contributors to the opal signal in the 20-40 µm fraction, while they share dominance with radiolarians in the >40 µm fraction. Total diatom concentrations average 1.2*10**4 valves/g in the 20-40 µm fraction and 4*10**3 valves/g in the >40 µm fraction. Robust taxa of warm water affinity (Alveus marinus, Azpeitia neocrenulata, Azpeitia nodulifera and Roperia tesselata) characterize the surface sediments. In contrast, biogenic opal content in brine surface sediment samples is much higher than in the non-brine samples, ranging from 2.8 to 3.8 wt.% SiO2, and concentration of siliceous microorganisms is 3-4 orders of magnitude higher. In addition here, diatoms dominate the opal signal. The taxa found in these samples are a mixture of non-brine and plankton samples, and fragile forms (e.g., N. bicapitata group, Neodelphineis indica) are well preserved in these sediments. Thus, brine sediments in this region seem to offer a great potential for palaeoenvironmental studies.

Temperature time series and physical properties of snow samples of a high-arctic permafrost site on Svalbard, Norway

Independent measurements of radiation, sensible and latent heat fluxes and the ground heat flux are used to describe the annual cycle of the surface energy budget at a high-arctic permafrost site on Svalbard. During summer, the net short-wave radiation is the dominant energy source, while well developed turbulent processes and the heat flux in the ground lead to a cooling of the surface. About 15% of the net radiation is consumed by the seasonal thawing of the active layer in July and August. The Bowen ratio is found to vary between 0.25 and 2, depending on water content of the uppermost soil layer. During the polar night in winter, the net long-wave radiation is the dominant energy loss channel for the surface, which is mainly compensated by the sensible heat flux and, to a lesser extent, by the ground heat flux, which originates from the refreezing of the active layer. The average annual sensible heat flux of -6.9 W/m**2 is composed of strong positive fluxes in July and August, while negative fluxes dominate during the rest of the year. With 6.8 W/m**2, the latent heat flux more or less compensates the sensible heat flux in the annual average. Strong evaporation occurs during the snow melt period and particularly during the snow-free period in summer and fall. When the ground is covered by snow, latent heat fluxes through sublimation of snow are recorded, but are insignificant for the average surface energy budget. The near-surface atmospheric stratification is found to be predominantly unstable to neutral, when the ground is snow-free, and stable to neutral for snow-covered ground. Due to long-lasting near-surface inversions in winter, an average temperature difference of approximately 3 K exists between the air temperature at 10 m height and the surface temperature of the snow.

Combined effects of short-term ocean acidification and heat shock in a benthic copepod Tigriopus japonicus Mori

Warming of the world's oceans is predicted to have many negative effects on organisms as they have optimal thermal windows. In coastal waters, however, both temperatures and pCO2 (pH) exhibit diel variations, and biological performances are likely to be modulated by physical and chemical environmental changes. To understand how coastal zooplankton respond to the combined impacts of heat shock and increased pCO2, the benthic copepod Tigriopus japonicus were treated at temperatures of 24, 28, 32 and 36 °C to simulate natural coastal temperatures experienced in warming events, when acclimated in the short term to either ambient (LC, 390 µatm) or future CO2 (HC, 1000 µatm). HC and heat shock did not induce any mortality of T. japonicus, though respiration increased up to 32 °C before being depressed at 36 °C. Feeding rate peaked at 28 °C but did not differ between CO2 treatments. Expression of heat shock proteins (hsps mRNA) was positively related to temperature, with no significant differences between the CO2 concentrations. Nauplii production was not affected across all treatments. Our results demonstrate that T. japonicus responds more sensitively to heat shocks rather than to seawater acidification; however, ocean acidification may synergistically act with ocean warming to mediate the energy allocation of copepods.

Last Interglacial synthesis of high-latitude temperature: temperature anomalies and associated errors for 4 time slices

The Last Interglacial (LIG, 129-116 thousand of years BP, ka) represents a test bed for climate model feedbacks in warmer-than-present high latitude regions. However, mainly because aligning different palaeoclimatic archives and from different parts of the world is not trivial, a spatio-temporal picture of LIG temperature changes is difficult to obtain. Here, we have selected 47 polar ice core and sub-polar marine sediment records and developed a strategy to align them onto the recent AICC2012 ice core chronology. We provide the first compilation of high-latitude temperature changes across the LIG associated with a coherent temporal framework built between ice core and marine sediment records. Our new data synthesis highlights non-synchronous maximum temperature changes between the two hemispheres with the Southern Ocean and Antarctica records showing an early warming compared to North Atlantic records. We also observe warmer than present-day conditions that occur for a longer time period in southern high latitudes than in northern high latitudes. Finally, the amplitude of temperature changes at high northern latitudes is larger compared to high southern latitude temperature changes recorded at the onset and the demise of the LIG. We have also compiled four data-based time slices with temperature anomalies (compared to present-day conditions) at 115 ka, 120 ka, 125 ka and 130 ka and quantitatively estimated temperature uncertainties that include relative dating errors. This provides an improved benchmark for performing more robust model-data comparison. The surface temperature simulated by two General Circulation Models (CCSM3 and HadCM3) for 130 ka and 125 ka is compared to the corresponding time slice data synthesis. This comparison shows that the models predict warmer than present conditions earlier than documented in the North Atlantic, while neither model is able to produce the reconstructed early Southern Ocean and Antarctic warming. Our results highlight the importance of producing a sequence of time slices rather than one single time slice averaging the LIG climate conditions.

First-year and second-year snow properties on sea ice in the Weddell Sea during POLARSTERN cruise ANT-XXII/2 (ISPOL) during the drift

Observations of snow properties, superimposed ice, and atmospheric heat fluxes have been performed on first-year and second-year sea ice in the western Weddell Sea, Antarctica. Snow in this region is particular as it does usually survive summer ablation. Measurements were performed during Ice Station Polarstern (ISPOL), a 5-week drift station of the German icebreaker RV Polarstern. Net heat flux to the snowpack was 8 W/m**2, causing only 0.1 to 0.2 m of thinning of both snow cover types, thinner first-year and thicker second-year snow. Snow thinning was dominated by compaction and evaporation, whereas melt was of minor importance and occurred only internally at or close to the surface. Characteristic differences between snow on first-year and second-year ice were found in snow thickness, temperature, and stratigraphy. Snow on second-year ice was thicker, colder, denser, and more layered than on first-year ice. Metamorphism and ablation, and thus mass balance, were similar between both regimes, because they depend more on surface heat fluxes and less on underground properties. Ice freeboard was mostly negative, but flooding occurred mainly on first-year ice. Snow and ice interface temperature did not reach the melting point during the observation period. Nevertheless, formation of discontinuous superimposed ice was observed. Color tracer experiments suggest considerable meltwater percolation within the snow, despite below-melting temperatures of lower layers. Strong meridional gradients of snow and sea-ice properties were found in this region. They suggest similar gradients in atmospheric and oceanographic conditions and implicate their importance for melt processes and the location of the summer ice edge.

Degradation rates and productivity signals of organic-walled dinoflagellate cyst species in surface sediment samples

To understand the role of the ocean within the global carbon cycle, detailed information is required on key-processes within the marine carbon cycle; bio-production in the upper ocean, export of the produced material to the deep ocean and the storage of carbon in oceanic sediments. Quantification of these processes requires the separation of signals of net primary production and the rate of organic matter decay as reflected in fossil sediments. This study examines the large differences in degradation rates of organic-walled dinoflagellate cyst species to separate these degradation and productivity signals. For this, accumulation rates of cyst species known to be resistant (R-cysts) or sensitive (S-cysts) to aerobic degradation of 62 sites are compared to mean annual chlorophyll-a, sea-surface temperature, sea-surface salinity, nitrate and phosphate concentrations of the upper waters and deep-water oxygen concentrations. Furthermore, the degradation of sensitive cysts, as expressed by the degradation constant k and reaction time t, has been related to bottom water [O2]. The studied sediments were taken from the Arabian Sea, north-western African Margin (North Atlantic), western-equatorial Atlantic Ocean/Caraibic, south-western African margin (South Atlantic) and Southern Ocean (Atlantic sector). Significant relationships are observed between (a) accumulation rates of R-cysts and upper water chlorophyll-a concentrations, (b) accumulation rates of S-cysts and bottom water [O2] and (c) degradation rates of S-cysts (kt) and bottom water [O2]. Relationships that are extremely weak or are clearly insignificant on all confidence intervals are between (1) S-cyst accumulation rates and chlorophyll-a concentrations, sea-surface temperature (SST), sea-surface salinity (SSS), phosphate concentrations (P) and nitrate concentrations (N), (2) between R-cyst accumulation rates and bottom water [O2], SST, SSS, P and N, and between (3) kt and water depth. Co-variance is present between the parameters N and P, N, P and chlorophyll-a, oxygen and water depth. Correcting for this co-variance does not influence the significance of the relationship given above. The possible applicability of dinoflagellate cyst degradation to estimate past net primary production and deep ocean ventilation is discussed.

(Fig. 1D, 3E) Sea surface temperature and alkenone C37:4 for sediment core MD01-2444 during the penultimate glaciation

Past glacials can be thought of as natural experiments in which variations in boundary conditions influenced the character of climate change. However, beyond the last glacial, an integrated view of orbital- and millennial-scale changes and their relation to the record of glaciation has been lacking. Here, we present a detailed record of variations in the land-ocean system from the Portuguese margin during the penultimate glacial and place it within the framework of ice-volume changes, with particular reference to European ice-sheet dynamics. The interaction of orbital- and millennial-scale variability divides the glacial into an early part with warmer and wetter overall conditions and prominent climate oscillations, a transitional mid-part, and a late part with more subdued changes as the system entered a maximum glacial state. The most extreme event occurred in the mid-part and was associated with melting of the extensive European ice sheet and maximum discharge from the Fleuve Manche river. This led to disruption of the meridional overturning circulation, but not a major activation of the bipolar seesaw. In addition to stadial duration, magnitude of freshwater forcing, and background climate, the evidence also points to the influence of the location of freshwater discharges on the extent of interhemispheric heat transport.

Eocene-Oligocene sea surface temperature and pCO2 estimates

The long-term cooling trend of the Cenozoic is punctuated by shorter-term climatic events, such as the inception of permanent ice sheets on Antarctica at the Eocene?Oligocene Transition (~33.7 Ma). Taking advantage of the excellent state of preservation of coccolith calcite in equatorial Atlantic deep-sea cores, we unveil progressive tropical warming in the Atlantic Ocean initiated 4 million years prior to Antarctic glaciation. Warming preceding glaciation may appear counterintuitive, but we argue that this long-term climatic precursor to the EOT reinforced cooling of austral high latitudes via the redistribution of heat at the surface of the oceans. We discuss this new prominent paleoceanographic and climatic feature in the context of overarching pCO2 decline and the establishment of an Antarctic circumpolar current.