Water temperature reconstructions based on paired Uk'37 and Tex86 records, and paired surface and subsurface Mg/Ca records

This dataset contains the collection of available published paired Uk'37 and Tex86 records spanning multi-millennial to multi-million year time scales, as well as a collection of Mg/Ca-derived temperatures measured in parallel on surface and subsurface dwelling foraminifera, both used in the analyses of Ho and Laepple, Nature Geoscience 2016. As the signal-to-noise ratios of proxy-derived Holocene temperatures are relatively low, we selected records that contain at least the last deglaciation (oldest sample >18kyr BP).

Mg/Ca ratio of Globorotalia tumida and subsurface temperature reconstruction for ODP Site 130-806

The tropical Pacific thermocline strength, depth, and tilt are critical to tropical mean state and variability. During the early Pliocene (~3.5 to 4.5 Ma), the Eastern Equatorial Pacific (EEP) thermocline was deeper and the cold tongue was warmer than today, which resulted in an mean state with a reduced zonal sea surface temperature gradient, or El Padre. However, it is unclear whether the deep thermocline was a local feature of the EEP or a basin-wide condition with global implications. Our measurements of Mg/Ca of Globorotalia tumida in a western equatorial Pacific site indicate Pliocene subsurface temperatures warmer than today; thus, El Padre included a basin-wide thermocline that was relatively warm, deep, and weakly tilted. At ~4 Ma, thermocline steepening was coupled to cooling of the cold tongue. Since ~4 Ma, the basin-wide thermocline cooled/shoaled gradually, with implications for thermocline feedbacks in tropical dynamics and the interpretation of TEX86-derived temperatures.

Geochemical and physical sediment properties of sediment cores from the Portuguese shelf

Marine sediments from the Portuguese shelf are influenced by environmental changes in the surrounding continental and marine environment. These are largely controlled by the North Atlantic Oscillation, but additional impacts may arise from episodic tsunamis. In order to investigate these influences, a high resolution multi-proxy study has been carried out on a 5.4 m long gravity core and five box cores from the Tagus prodelta on the western Portuguese margin, incorporating geochemical (Corg/Ntotal ratios, d13Corg, d15N, d18O, Corg and CaCO3 content) and physical sediment properties (magnetic susceptibility, grain-size). Subsurface data of the five box cores indicate no major effect of early postdepositional alteration. Surface data show a higher fraction of terrigenous organic material close to the river mouth and in the southern prodelta. Gravity core GeoB 8903 covers the last 3.2 kyrs with a temporal resolution of at least 0.1 cm/yr. Very high sedimentation rates between 69 and 140 cm core depth indicate a possible disturbance of the record by the AD1755 tsunami, although no evidence for a disturbance is observed in the data. Sea surface temperature and salinity on the prodelta, the local budget of marine NO3- as well as the provenance of organic matter remained virtually constant during the past 3.2 kyrs. A positive correlation between magnetic susceptibility (MS) and North Atlantic Oscillation (NAO) is evident for the past 250 years, coinciding with a negative correlation between mean grain-size and NAO. This is assigned to a constant riverine supply of fine material with high MS, which is diluted by the riverine input of a coarser, low-MS component during NAO negative, high-precipitation phases. End-member modelling of the lithic grain-size spectrum supports this, revealing a third, coarse lithic component. The high abundance of this coarse end-member prior to 2 kyr BP is interpreted as the result of stronger bottom currents, concentrating the coarse sediment fraction by winnowing. As continental climate was more arid prior to 2 kyr BP (Subboreal), the coarse end-member may also consist of dust from local sources. A decrease in grain-size and CaCO3 content after 2 kyr BP is interpreted as a result of decreasing wind strength. The onset of a fining trend and a further decrease in CaCO3 around AD900 occurs simultaneous to climatic variations, reconstructed from eastern North Atlantic records. A strong increase in MS between AD1400 and AD1500 indicates higher lithic terrigenous input, caused by deforestation in the hinterland.

Isotopic signature of specific biomarkers derived from anaerobic methanotrophic archaea (ANME groups) and sulphate-reducing bacteria (SRB) at different cold seep provinces

Anaerobic methane-oxidizing microbial communities in sediments at cold methane seeps are important factors in controlling methane emission to the ocean and atmosphere. Here, we investigated the distribution and carbon isotopic signature of specific biomarkers derived from anaerobic methanotrophic archaea (ANME groups) and sulphate-reducing bacteria (SRB) responsible for the anaerobic oxidation of methane (AOM) at different cold seep provinces of Hydrate Ridge, Cascadia margin. The special focus was on their relation to in situ cell abundances and methane turnover. In general, maxima in biomarker abundances and minima in carbon isotope signatures correlated with maxima in AOM and sulphate reduction as well as with consortium biomass. We found ANME-2a/DSS aggregates associated with high abundances of sn-2,3-di-O-isoprenoidal glycerol ethers (archaeol, sn-2-hydroxyarchaeol) and specific bacterial fatty acids (C16:1omega5c, cyC17:0omega5,6) as well as with high methane fluxes (Beggiatoa site). The low to medium flux site (Calyptogena field) was dominated by ANME-2c/DSS aggregates and contained less of both compound classes but more of AOM-related glycerol dialkyl glycerol tetraethers (GDGTs). ANME-1 archaea dominated deeper sediment horizons at the Calyptogena field where sn-1,2-di-O-alkyl glycerol ethers (DAGEs), archaeol, methyl-branched fatty acids (ai-C15:0, i-C16:0, ai-C17:0), and diagnostic GDGTs were prevailing. AOM-specific bacterial and archaeal biomarkers in these sediment strata generally revealed very similar d13C-values of around -100 per mill. In ANME-2-dominated sediment sections, archaeal biomarkers were even more 13C-depleted (down to -120 per mill), whereas bacterial biomarkers were found to be likewise 13C-depleted as in ANME-1-dominated sediment layers (d13C: -100 per mill). The zero flux site (Acharax field), containing only a few numbers of ANME-2/DSS aggregates, however, provided no specific biomarker pattern. Deeper sediment sections (below 20 cm sediment depth) from Beggiatoa covered areas which included solid layers of methane gas hydrates contained ANME-2/DSS typical biomarkers showing subsurface peaks combined with negative shifts in carbon isotopic compositions. The maxima were detected just above the hydrate layers, indicating that methane stored in the hydrates may be available for the microbial community. The observed variations in biomarker abundances and 13C-depletions are indicative of multiple environmental and physiological factors selecting for different AOM consortia (ANME-2a/DSS, ANME-2c/DSS, ANME-1) along horizontal and vertical gradients of cold seep settings.

Foraminifer counts in core JM01-1199, and review of radiocarbon ages in the Norwegian-Barents Sea-Svalbard margin

We compare six high-resolution Holocene, sediment cores along a S-N transect on the Norwegian-Svalbard continental margin from ca 60°N to 77.4°N, northern North Atlantic. Planktonic foraminifera in the cores were investigated to show the changes in upper surface and subsurface water mass distribution and properties, including summer sea-surface temperatures (SST). The cores are located below the axis of the Norwegian Current and the West Spitsbergen Current, which today transport warm Atlantic Water to the Arctic. Sediment accumulation rates are generally high at all the core sites, allowing for a temporal resolution of 10-102 years. SST is reconstructed using different types of transfer functions, resulting in very similar SST trends, with deviations of no more than +- 1.0/1.5 °C. A transfer function based on the maximum likelihood statistical approach is found to be most relevant. The reconstruction documents an abrupt change in planktonic foraminiferal faunal composition and an associated warming at the Younger Dryas-Preboreal transition. The earliest part of the Holocene was characterized by large temperature variability, including the Preboreal Oscillations and the 8.2 k event. In general, the early Holocene was characterized by SSTs similar to those of today in the south and warmer than today in the north, and a smaller S-N temperature gradient (0.23 °C/°N) compared to the present temperature gradient (0.46 °C/°N). The southern proxy records (60-69°N) were more strongly influenced by slightly cooler subsurface water probably due to the seasonality of the orbital forcing and increased stratification due to freshening. The northern records (72-77.4°N) display a millennial-scale change associated with reduced insolation and a gradual weakening of the North Atlantic thermohaline circulation (THC). The observed northwards amplification of the early Holocene warming is comparable to the pattern of recent global warming and future climate modelling, which predicts greater warming at higher latitudes. The overall trend during mid and late Holocene was a cooling in the north, stable or weak warming in the south, and a maximum S-N SST gradient of ca 0.7 °C/°N at 5000 cal. years BP. Superimposed on this trend were several abrupt temperature shifts. Four of these shifts, dated to 9000-8000, 5500-3000 and 1000 and ~400 cal. years BP, appear to be global, as they correlate with periods of global climate change. In general, there is a good correlation between the northern North Atlantic temperature records and climate records from Norway and Svalbard.

Occurrence of late Aptian/earliest Albian radiolaria from ODP Hole 113-693A (Table 2)

Unusually well preserved Cretaceous radiolarians are observed in the subsurface sections from two drilled sites in the Weddell Sea collected during Leg 113 of the Ocean Drilling Program. Radiolarians from the lithified calcareous chalk of Hole 689B represent the first Campanian-Maestrichtian assemblage which is characterized by abundant Cromyodruppa Iconcentrica, Dictyomitra multicostata, and Protostichocapsa stocki. Abundant Pseudodictyomitra pentacolaensis and Diacanthocapsa sp. 1, on the other hand, are the main constituents of the assemblage from the latest Aptian/earliest Albian diatomite of Hole 693B. These represent the oldest and the highest-latitude reported radiolarian occurrences from the Atlantic sector of the Antarctic Ocean. The assemblages are marked by their low diversity and an absence of low- to mid-latitude zonal indices.

(Table 1) Summary of sampling points described and main petrographic, morphological and geochemical characteristics

Targeted sampling on the Dolgovskoy Mound (northern Shatsky Ridge) revealed the presence of spectacular laterally extensive and differently shaped authigenic carbonates. The sampling stations were selected based on sidescan sonar and profiler images that show patchy backscatter and irregular and discontinuous reflections in the near subsurface. The interpretation of acoustic data from the top part of the mound supports the seafloor observations and the sampling that revealed the presence of a complex subsurface plumbing system characterized by carbonates and gas. The crusts sampled consist of carbonate cemented layered hemipelagic sedimentary Unit 1 associated with several centimetres thick microbial mats. Three different carbonate morphologies were observed: (a) tabular slabs, (b) subsurface cavernous carbonates consisting of void chambers up to 20 cm**3 in size and (c) chimney and tubular conduits vertically oriented or forming a subhorizontal network in the subsurface. The methanogenic origin of the carbonates is established based on visual observations of fluids seepage structures, 13C depletion of the carbonates (d13C varying between -36.7 per mil and -27.4 per mil), and by thin carbonate layers present within the thick microbial mats. Laboratory experiments with a Hele-Shaw cell were conducted in order to simulate the gas seepage through contrasting grain size media present on the seafloor. Combined petrography, visual observations and sandbox simulations allowed a characterization of the dynamics and the structures of the plumbing system in the near subsurface. Based on sample observations and the experiments, three observed morphologies of authigenic carbonates are interpreted, respectively, as (a) Darcian porous flow through the finely laminated clayey/coccolith-rich layers, (b) gas accumulation chambers at sites where significant fluid escape was impeded by thicker clayey layers forming the laminated Unit1 and (c) focussed vertical fluid venting and subhorizontal migration of overpressured fluids released from (b). The Hele-Shaw cell experiments represent a promising tool for investigating shallow fluid flow pathways in marine systems.

Chemistry of sediment profile GeoB4906

Mineralization of organic matter and the subsequent dissolution of calcite were simulated for surface sediments of the upper continental slope off Gabon by using microsensors to measure O2, pH, pCO2 and Ca2+ (in situ), pore-water concentration profiles of NO3-, NH4+, Fe2+, and Mn2+ and SO42- (ex situ), as well as sulfate reduction rates derived from incubation experiments. The transport and reaction model CoTReM was used to simulate the degradation of organic matter by O2, [NO3]-, Fe(OH)3 and [SO4]2-, reoxidation reactions involving Fe2+ and Mn2+, and precipitation of FeS. Model application revealed an overall rate of organic matter mineralization amounting to 50 µmol C cm**-2 yr**-1, of which 77% were due to O2, 17% to [NO3]- and 3% to Fe(OH)3 and 3% to [SO4]2-. The best fit for the pH profile was achieved by adapting three different dissolution rate constants of calcite ranging between 0.01 and 0.5% d-1 and accounting for different calcite phases in the sediment. A reaction order of 4.5 was assumed in the kinetic rate law. A CaCO3 flux to the sediment was estimated to occur at a rate of 42 g m**-2 yr**-1 in the area of equatorial upwelling. The model predicts a redissolution flux of calcite amounting to 36 g m**-2 yr**-1, thus indicating that ~90% of the calcite flux to the sediment is redissolved.