Sediment core-top GDGT data from the Arctic, the North Pacific and the Southern Ocean

TEX86 (TetraEther indeX of tetraethers consisting of 86 carbon atoms) is a sea surface temperature (SST) proxy based on the distribution of archaeal isoprenoid glycerol dialkyl glycerol tetraethers (GDGTs). In this study, we appraise the applicability of TEX86 and TEX86L in subpolar and polar regions using surface sediments. We present TEX86 and TEX86L data from 160 surface sediment samples collected in the Arctic, the Southern Ocean and the North Pacific. Most of the SST estimates derived from both TEX86 and TEX86L are anomalously high in the Arctic, especially in the vicinity of Siberian river mouths and the sea ice margin, plausibly due to additional archaeal contributions linked to terrigenous input. We found unusual GDGT distributions at five sites in the North Pacific. High GDGT-0/crenarchaeol and GDGT-2/crenarchaeol ratios at these sites suggest a substantial contribution of methanogenic and/or methanotrophic archaea to the sedimentary GDGT pool here. Apart from these anomalous findings, TEX86 and TEX86L values in the surface sediments from the Southern Ocean and the North Pacific do usually vary with overlaying SSTs. In these regions, the sedimentary TEX86-SST relationship is similar to the global calibration, and the derived temperature estimates agree well with overlaying annual mean SSTs at the sites. However, there is a systematic offset between the regional TEX86L-SST relationships and the global calibration. At these sites, temperature estimates based on the global TEX86L calibration are closer to summer SSTs than annual mean SSTs. This finding suggests that in these subpolar settings a regional TEX86L calibration may be a more suitable equation for temperature reconstruction than the global calibration.

Epibenthic foraminiferal d13C in the Recent deep Arctic Ocean

Low planktic and benthic d18O and d13C values in sediments from the Nordic seas of cold stadials of the last glaciation have been attributed to brines, formed similar to modern ones in the Arctic Ocean. To expand on the carbon isotopes of this hypothesis I investigated benthic d13C from the modern Arctic Ocean. I show that mean d13C values of live epibenthic foraminifera from the deep Arctic basins are higher than mean d13C values of upper slope epibenthic foraminifera. This agrees with mean high d13C values of dissolved inorganic carbon (DIC) in Arctic Bottom Water (ABW), which are higher than mean d13CDIC values from shallower water masses of mainly Atlantic origin. However, adjustments for oceanic 13C-Suess depletion raise subsurface and intermediate water d13CDIC values over ABW d13CDIC ones. Accordingly, during preindustrial Holocene times, the d13CDIC of ABW was as high or higher than today, but lower than the d13CDIC of younger subsurface and intermediate water. If brine-enriched water significantly ventilated ABW, brines should have had high d13CDIC values. Analogously, high-d13CDIC brines may have been formed in the Nordic seas during warm interstadials. During cold stadials, when most of the Arctic Ocean was perennially sea-ice covered, a cessation of high-d13CDIC brine rejection may have lowered d13CDIC values of ABW, and ultimately the d13CDIC in Nordic seas intermediate and deep water. So, in contrast to the idea of enhanced brine formation during cold stadials, the results of this investigation imply that a cessation of brine rejection would be more likely.

Diatom abundance in surface sediments of the Southern Ocean

Based on the quantitative analysis of diatom assemblages preserved in 274 surface sediment samples recovered in the Pacific, Atlantic and western Indian sectors of the Southern Ocean we have defined a new reference database for quantitative estimation of late-middle Pleistocene Antarctic sea ice fields using the transfer function technique. The Detrended Canonical Analysis (DCA) of the diatom data set points to a unimodal distribution of the diatom assemblages. Canonical Correspondence Analysis (CCA) indicates that winter sea ice (WSI) but also summer sea surface temperature (SSST) represent the most prominent environmental variables that control the spatial species distribution. To test the applicability of transfer functions for sea ice reconstruction in terms of concentration and occurrence probability we applied four different methods, the Imbrie and Kipp Method (IKM), the Modern Analog Technique (MAT), Weighted Averaging (WA), and Weighted Averaging Partial Least Squares (WAPLS), using logarithm-transformed diatom data and satellite-derived (1981-2010) sea ice data as a reference. The best performance for IKM results was obtained using a subset of 172 samples with 28 diatom taxa/taxa groups, quadratic regression and a three-factor model (IKM-D172/28/3q) resulting in root mean square errors of prediction (RMSEP) of 7.27% and 11.4% for WSI and summer sea ice (SSI) concentration, respectively. MAT estimates were calculated with different numbers of analogs (4, 6) using a 274-sample/28-taxa reference data set (MAT-D274/28/4an, -6an) resulting in RMSEP's ranging from 5.52% (4an) to 5.91% (6an) for WSI as well as 8.93% (4an) to 9.05% (6an) for SSI. WA and WAPLS performed less well with the D274 data set, compared to MAT, achieving WSI concentration RMSEP's of 9.91% with WA and 11.29% with WAPLS, recommending the use of IKM and MAT. The application of IKM and MAT to surface sediment data revealed strong relations to the satellite-derived winter and summer sea ice field. Sea ice reconstructions performed on an Atlantic- and a Pacific Southern Ocean sediment core, both documenting sea ice variability over the past 150,000 years (MIS 1 - MIS 6), resulted in similar glacial/interglacial trends of IKM and MAT-based sea-ice estimates. On the average, however, IKM estimates display smaller WSI and slightly higher SSI concentration and probability at lower variability in comparison with MAT. This pattern is a result of different estimation techniques with integration of WSI and SSI signals in one single factor assemblage by applying IKM and selecting specific single samples, thus keeping close to the original diatom database and included variability, by MAT. In contrast to the estimation of WSI, reconstructions of past SSI variability remains weaker. Combined with diatom-based estimates, the abundance and flux pattern of biogenic opal represents an additional indication for the WSI and SSI extent.

Registry of selected samples from the Tara Oceans Expedition (2009-2013)

The present data set is a registry of samples from the Tara Oceans Expedition (2009-2013) that were selected for publication in a special issue of the SCIENCE journal (see related references below). The registry provides details about the sampling location and methodology of each sample. Uniform resource locators (URLs) offer direct links to additional contextual environmental data and to the corresponding sequence runs used for analysis in the related literature publications in the SCIENCE journal.

Biodiversity context of selected samples from the Tara Oceans Expedition (2009-2013)

The present data set provides contextual data for samples from the Tara Oceans Expedition (2009-2013) that were selected for publication in a special issue of the SCIENCE journal (see related references below). Contextual data include various diversity indexes calculated for the sampling location using satellite and model climatologies (Darwin project, Physat) and results from the sequencing of Tara Oceans samples.

Particle flux in the NE Atlantic as recorded by sediment traps

Downward particle flux was measured using sediment traps at various depths over the Porcupine Abyssal Plain (water depth ab. 4850 m) for prolonged periods from 1989 to 1999. A strong seasonal pattern of flux was evident reaching a maximum in mid-summer. The composition of the material changed with depth, reflecting the processes of remineralisation and dissolution as the material sank through the water column. However, there was surprisingly little seasonal variation in its composition to reflect changes in the biology of the euphotic zone. Currents at the site have a strong tidal component with speeds almost always less than 15 cm/sec. In the deeper part of the water column they tend to be northerly in direction, when averaged over periods of several months. A model of upper ocean biogeochemistry forced by meteorology was run for the decade in order to provide an estimate of flux at 3000 m depth. Agreement with measured organic carbon flux is good, both in terms of the timings of the annual peaks and in the integrated annual flux. Interannual variations in the integrated flux are of similar magnitude for both the model output and sediment trap measurements, but there is no significant relationship between these two sets of estimates. No long-term trend in flux is evident, either from the model, or from the measurements. During two spring/summer periods, the marine snow concentration in the water column was assessed by time-lapse photography and showed a strong peak at the start of the downward pulse of material at 3000 m. This emphasises the importance of large particles during periods of maximum flux and at the start of flux peaks. Time lapse photographs of the seabed show a seasonal cycle of coverage of phytodetrital material, in agreement with the model output both in terms of timing and magnitude of coverage prior to 1996. However, after a change in the structure of the benthic community in 1996 no phytodetritus was evident on the seabed. The model output shows only a single peak in flux each year, whereas the measured data usually indicated a double peak. It is concluded that the observed double peak may be a reflection of lowered sediment trap efficiency when flux is very high and is dominated by large marine snow particles. Resuspension into the trap 100 m above the seabed, when compared to the primary flux at 3000 m depth (1800 mab) was lower during periods of high primary flux probably because of a reduction in the height of resuspension when the material is fresh. At 2 mab, the picture is more complex with resuspension being enhanced during the periods of higher flux in 1997, which is consistent with this hypothesis. However there was rather little relationship to flux at 3000 m in 1998. At 3000 m depth, the Flux Stability Index (FSI), which provides a measure of the constancy of the seasonal cycle of flux, exhibited an inverse relationship with flux, such that the highest flux of organic carbon was recorded during the year with the greatest seasonal variation.

Organic carbon flux and remineralization in surface sediments

Organic carbon fluxes through the sediment/water interface in the high-latitude North Atlantic were calculated from oxygen microprofiles. A wire-operated in situ oxygen bottom profiler was deployed, and oxygen profiles were also measured onboard (ex situ). Diffusive oxygen fluxes, obtained by fitting exponential functions to the oxygen profiles, were translated into organic carbon fluxes and organic carbon degradation rates. The mean Corg input to the abyssal plain sediments of the Norwegian and Greenland Seas was found to be 1.9 mg C/m**2/d. Typical values at the seasonally ice-covered East Greenland continental margin are between 1.3 and 10.9 mg C/m**2/d (mean 3.7 mg C/m**2/d), whereas fluxes on the East Greenland shelf are considerably higher, 9.1-22.5 mg C/m**2/d. On the Norwegian continental slope Corg fluxes of 3.3-13.9 mg C/m**2/d (mean 6.5 mg C/m**2/d) were found. Fluxes are considerably higher here compared to stations on the East Greenland slope at similar water depths. By repeated occupation of three sites off southern Norway in 1997 the temporal variability of diffusive O2 fluxes was found to be quite low. The seasonal signal of primary and export production from the upper water column appears to be strongly damped at the seafloor. Degradation rates of 0.004-1.1 mg C/cm**3/a at the sediment surface were calculated from the oxygen profiles. First-order degradation constants, obtained from Corg degradation rates and sediment organic carbon content, are in the range 0.03-0.6/a. Thus, the corresponding mean lifetime of organic carbon lies between 1.7 and 33.2 years, which also suggests that seasonal variations in Corg flux are small. The data presented here characterize the Norwegian and Greenland Seas as oligotrophic and relatively low organic carbon deep-sea environments.

Biomass of large protists from the Infrakingdom Rhizaria in the global ocean

The present data set is a worldwide compilation from 11 oceanographic expeditions during which an underwater vision profiler (UVP) was deployed in situ to determine the vertical distribution (biomass) of 4 taxonomic groups of plankton larger than 600 µm, belonging to the Infrakingdom Rhizaria, including Collodaria, Acantharia, Phaeodaria and other Rhizaria. Vertical distributions are binned in four layers: 0-100, 0-200, 100-500 and 0-500 m.

Magnetic susceptibility and ice-rafted debris in surface sediments of the Nordic Sea

Magnetic susceptibility and ice-rafted debris of surface sediments in the Nordic Seas were investigated to reconstruct source areas and recent transport pathways of magnetic minerals. From the distribution of magnetic susceptibility and ice-rafted debris and published data on petrographic tracers for iceberg drift, we reconstructed a counter-clockwise iceberg drift pattern during cooler phases in the Holocene, which is similar to conceptual and numerical models for Weichselian iceberg drift. The release of basaltic debris at Scoresby Sund played a significant role for the magnetic signature of stadial/interstadial events during isotope stage 3 recorded in sediment cores of the Nordic Seas.

Abundance, biomass and secondary production of benthic megafauna on the Barents Sea shelf in 2008 and 2009

Megabenthos plays a major role in the overall energy flow on Arctic shelves, but information on megabenthic secondary production on large spatial scales is scarce. Here, we estimated for the first time megabenthic secondary production for the entire Barents Sea shelf by applying a species-based empirical model to an extensive dataset from the joint Norwegian? Russian ecosystem survey. Spatial patterns and relationships were analyzed within a GIS. The environmental drivers behind the observed production pattern were identified by applying an ordinary least squares regression model. Geographically weighted regression (GWR) was used to examine the varying relationship of secondary production and the environment on a shelfwide scale. Significantly higher megabenthic secondary production was found in the northeastern, seasonally ice-covered regions of the Barents Sea than in the permanently ice-free southwest. The environmental parameters that significantly relate to the observed pattern are bottom temperature and salinity, sea ice cover, new primary production, trawling pressure, and bottom current speed. The GWR proved to be a versatile tool for analyzing the regionally varying relationships of benthic secondary production and its environmental drivers (R² = 0.73). The observed pattern indicates tight pelagic? benthic coupling in the realm of the productive marginal ice zone. Ongoing decrease of winter sea ice extent and the associated poleward movement of the seasonal ice edge point towards a distinct decline of benthic secondary production in the northeastern Barents Sea in the future.