Assamblages and cluster analysis from DSDP Hole 22-214 of the eastern Indian Ocean

Multivariate analyses of latest Pliocene through Holocene benthic foraminifera from 61 samples from Deep-Sea Drilling Project (DSDP) Site 214, eastem Indian Ocean were carried out. The 46 highest ranked species were used in R-mode factor analysis which has enabled to the identification of three environmentally significant assemblages at Site 214. Assemblage 1 is characterized by Uvigerina hispido-costata, Osangularia culter , Gavelinopsis lobatulus, Cibicides wuellerstorfi and Karreriella baccata as principal species. This assemblage is inferred to reflect high-energy, well-oxygenated and probably low-organic carbon deep-sea environment at Site 214. Assemblage 2 is defined principally by Globocassidulina pacifica and U. proboscidea and is considered to indicate an organic carbon-rich environment which resulted from high surface productivity irrespective of dissolved oxygen content. Assemblage 3 is marked by Oridorsalis umbonatus, Textularia lythostrota, Hoeglundina elegans, Pyrgo murrhina, and Pullenia quinqueloba as principal species. This assemblage is inferred to indicate a low-organic carbon environment with high pore water oxygen concentration leading to better preservation of deep-sea sediments.

Sea surface and the sea floor regionalization of the Southern Ocean by multivariate cluster analysis, links to ArcGIS project files

This study subdivides the Weddell Sea, Antarctica, into seafloor regions using multivariate statistical methods. These regions are categories used for comparing, contrasting and quantifying biogeochemical processes and biodiversity between ocean regions geographically but also regions under development within the scope of global change. The division obtained is characterized by the dominating components and interpreted in terms of ruling environmental conditions. The analysis uses 28 environmental variables for the sea surface, 25 variables for the seabed and 9 variables for the analysis between surface and bottom variables. The data were taken during the years 1983-2013. Some data were interpolated. The statistical errors of several interpolation methods (e.g. IDW, Indicator, Ordinary and Co-Kriging) with changing settings have been compared for the identification of the most reasonable method. The multivariate mathematical procedures used are regionalized classification via k means cluster analysis, canonical-correlation analysis and multidimensional scaling. Canonical-correlation analysis identifies the influencing factors in the different parts of the cove. Several methods for the identification of the optimum number of clusters have been tested. For the seabed 8 and 12 clusters were identified as reasonable numbers for clustering the Weddell Sea. For the sea surface the numbers 8 and 13 and for the top/bottom analysis 8 and 3 were identified, respectively. Additionally, the results of 20 clusters are presented for the three alternatives offering the first small scale environmental regionalization of the Weddell Sea. Especially the results of 12 clusters identify marine-influenced regions which can be clearly separated from those determined by the geological catchment area and the ones dominated by river discharge.

Ladderane lipid analysis and pore water and sediment chemistry of sediment cores from the continental shelf and slope off northwest Africa

The presence and abundance of anaerobic ammonium-oxidizing (anammox) bacteria was investigated in continental shelf and slope sediments (300-3000 m water depth) off northwest Africa in a combined approach applying quantitative polymerase chain reaction (q-PCR) analysis of anammox-specific 16S rRNA genes and anammox-specific ladderane biomarker lipids. We used the presence of an intact ladderane monoether lipid with a phosphocholine (PC) headgroup as a direct indicator for living anammox bacteria and compared it with the abundance of ladderane core lipids derived from both living and dead bacterial biomass. All investigated sediments contained ladderane lipids, both intact and core lipids, in agreement with the presence of anammoxspecific 16S rRNA gene copies, indicating that anammox occurs at all sites. Concentrations of ladderane core lipids in core top sediments varied between 0.3 and 97 ng g**-1 sediment, with the highest concentrations detected at the sites located on the shelf at shallower water depths between 300 and 500 m. In contrast, the C20 [3]-ladderane monoether-PC lipid was most abundant in a core top sediment from 1500 m water depth. Both anammox-specific 16S rRNA gene copy numbers and the concentration of the C20 [3]-ladderane monoether-PC lipid increased downcore in sediments located at greater water depths, showing highest concentrations of 1.2 x 10**8 copies g**-1 sediment and 30 pg g**-1 sediment, respectively, at the deepest station of 3000 m water depth. This suggests that the relative abundance of anammox bacteria is higher in sediments at intermediate to deep water depths where carbon mineralization rates are lower but where anammox is probably more important than denitrification.

(Table S2) Relative abundance of archaeal 16S rRNA genes in sediment cores

Deep drilling into the marine sea floor has uncovered a vast sedimentary ecosystem of microbial cells (Parkes et al., 1994, doi:10.1038/371410a0; D'Hondt et al., 2004, doi:10.1126/science.1101155). Extrapolation of direct counts of stained microbial cells to the total volume of habitable marine subsurface sediments suggests that between 56 Pg (Parkes et al., 1994, doi:10.1038/371410a0) and 303 Pg (Whitman et al., 1998) of cellular carbon could be stored in this largely unexplored habitat. From recent studies using various culture-independent techniques, no clear picture has yet emerged as to whether Archaea or Bacteria are more abundant in this extensive ecosystem (Schippers et al., doi:10.1038/nature03302; Inagaki et al., doi:10.1073/pnas.0511033103 ; Mauclaire et al., doi:10.1111/j.1472-4677.2004.00035.x; Biddle et al., doi:10.1073/pnas.0600035103). Here we show that in subsurface sediments buried deeper than 1 m in a wide range of oceanographic settings at least 87% of intact polar membrane lipids, biomarkers for the presence of live cells (Biddle et al., doi:10.1073/pnas.0600035103; Sturt et al., 2004, doi:10.1002/rcm.1378), are attributable to archaeal membranes, suggesting that Archaea constitute a major fraction of the biomass. Results obtained from modified quantitative polymerase chain reaction and slot-blot hybridization protocols support the lipid-based evidence and indicate that these techniques have previously underestimated archaeal biomass. The lipid concentrations are proportional to those of total organic carbon. On the basis of this relationship, we derived an independent estimate of amounts of cellular carbon in the global marine subsurface biosphere. Our estimate of 90 Pg of cellular carbon is consistent, within an order of magnitude, with previous estimates, and underscores the importance of marine subsurface habitats for global biomass budgets.

Pore water measurements of sediment cores from the Helgoland mud area, North Sea

Iron reduction in subseafloor sulfate-depleted and methane-rich marine sediments is currently a subject of interest in subsurface geomicrobiology. While iron reduction and microorganisms involved have been well studied in marine surface sediments, little is known about microorganisms responsible for iron reduction in deep methanic sediments. Here, we used quantitative PCR (Q-PCR)-based 16S rRNA gene copy numbers and pyrosequencing-based relative abundances of bacteria and archaea to investigate covariance between distinct microbial populations and specific geochemical profiles in the top 5 m of sediment cores from the Helgoland mud area, North Sea. We found that gene copy numbers of bacteria and archaea were specifically higher around the peak of dissolved iron in the methanic zone (250-350 cm. The higher copy numbers at these depths were also reflected by the relative sequence abundances of members of the candidate division JS1, methanogenic and Methanohalobium/ANME-3 related archaea. The distribution of these populations was strongly correlated to the profile of pore-water Fe2+ while that of Desulfobacteraceae corresponded to the pore-water sulfate profile. Furthermore, specific JS1 populations also strongly co-varied with the distribution of Methanosaetaceae in the methanic zone. Our data suggest that the interplay among JS1 bacteria, methanogenic archaea and Methanohalobium/ANME-3-related archaea may be important for iron reduction and methane cycling in deep methanic sediments of the Helgoland mud area and perhaps in other methane-rich depositional environments. .

Results of bulk sediment X-ray diffraction analysis and quantification of mineral phases based on the RockJock and on the QUAX quantitative analysis

We have performed quantitative X-ray diffraction (qXRD) analysis of 157 grab or core-top samples from the western Nordic Seas between (WNS) ~57°-75°N and 5° to 45° W. The RockJock Vs6 analysis includes non-clay (20) and clay (10) mineral species in the <2 mm size fraction that sum to 100 weight %. The data matrix was reduced to 9 and 6 variables respectively by excluding minerals with low weight% and by grouping into larger groups, such as the alkali and plagioclase feldspars. Because of its potential dual origins calcite was placed outside of the sum. We initially hypothesized that a combination of regional bedrock outcrops and transport associated with drift-ice, meltwater plumes, and bottom currents would result in 6 clusters defined by "similar" mineral compositions. The hypothesis was tested by use of a fuzzy k-mean clustering algorithm and key minerals were identified by step-wise Discriminant Function Analysis. Key minerals in defining the clusters include quartz, pyroxene, muscovite, and amphibole. With 5 clusters, 87.5% of the observations are correctly classified. The geographic distributions of the five k-mean clusters compares reasonably well with the original hypothesis. The close spatial relationship between bedrock geology and discrete cluster membership stresses the importance of this variable at both the WNS-scale and at a more local scale in NE Greenland.

Operator and replicability bias in comparative taphonomic studies

The operator effect is a well-known methodological bias already quantified in some taphonomic studies. However, the replicability effect, i.e., the use of taphonomic attributes as a replicable scientific method, has not been taken into account to the present. Here, we quantified for the first time this replicability bias using different multivariate statistical techniques, testing if the operator effect is related to the replicability effect. We analyzed the results reported by 15 operators working on the same dataset. Each operator analyzed 30 biological remains (bivalve shells) from five different sites, considering the attributes fragmentation, edge rounding, corrasion, bioerosion and secondary color. The operator effect followed the same pattern reported in previous studies, characterized by a worse correspondence for those attributes having more than two levels of damage categories. However, the effect did not appear to have relation with the replicability effect, because nearly all operators found differences among sites. Despite the binary attribute bioerosion exhibited 83% of correspondence among operators it was the taphonomic attributes that showed the highest dispersion among operators (28%). Therefore, we conclude that binary attributes (despite showing a reduction of the operator effect) diminish replicability, resulting in different interpretations of concordant data. We found that a variance value of nearly 8% among operators, was enough to generate a different taphonomic interpretation, in a Q-mode cluster analysis. The results reported here showed that the statistical method employed influences the level of replicability and comparability of a study and that the availability of results may be a valid alternative to reduce bias.

Geochemical provenance analysis of fine-grained sediment of sediment core AND1-1

Bulk chemical fine-grained sediment compositions from southern Victoria Land glacimarine sediments provide significant constraints on the reconstruction of sediment provenance models in the McMurdo Sound during Late Cenozoic time. High-resolution (~ 1 ka) geochemical data were obtained with a non-destructive AVAATECH XRF Core Scanner (XRF-CS) on the 1285 m long ANDRILL McMurdo Ice Shelf Project (MIS) sediment core AND-1B. This data set is complemented by high-precision chemical analyses (XRF and ICP-OES) on discrete samples. Statistical analyses reveal three geochemical facies which are interpreted to represent the following sources for the sediments recovered in the AND-1B core: 1) local McMurdo Volcanic Group (MVG) rocks, 2) Transantarctic Mountain rocks west of Ross Island (W TAM), and 3) Transantarctic Mountain rocks from more southerly areas (S TAM). Data indicate in combination with other sediment facies analyses (McKay et al., 2009, doi:10.1130/B26540.1) and provenance scenarios (Talarico and Sandroni, 2009, doi:10.1016/j.gloplacha.2009.04.007) that diamictites at the drill site are largely dominated by local sources (MVG) and are interpreted to indicate cold polar conditions with dry-based ice. MVG is interpreted to indicate cold polar condition with dry-based ice. A mixture of MVG and W TAM is interpreted to represent polar conditions and the S TAM facies is interpreted to represent open-marine conditions. Down-core variations in geochemical facies in the AND-1B core are interpreted to represent five major paleoclimate phases over the past 14 Ma. Cold polar conditions with major MVG influence occur below 1045 mbsf and above 120 mbsf. A section of warmer climate conditions with extensive peaks of S TAM influence characterizes the rest of the core, which is interrupted by a section from 525 to 855 mbsf of alternating influences of MVG and W TAM.

Distribution of deep-sea foraminifera in surface sediments east of New Zealand

Paleobathymetric assessments of fossil foraminiferal faunas play a significant role in the analysis of the paleogeographic, sedimentary, and tectonic histories of New Zealand's Neogene marine sedimentary basins. At depths >100 m, these assessments often have large uncertainties. This study, aimed at improving the precision of paleodepth assessments, documents the present-day distribution of deep-sea foraminifera (>63 µm) in 66 samples of seafloor sediment at 90-700 m water depth (outer shelf to mid-abyssal), east of New Zealand. One hundred and thirty-nine of the 465 recorded species of benthic foraminifera are new records for the New Zealand region. Characters of the foraminiferal faunas which appear to provide the most useful information for estimating paleobathymetry are, in decreasing order of reliability: relative abundance of common benthic species; benthic species associations; upper depth limits of key benthic species; and relative abundance of planktic foraminifera. R mode cluster analysis on the quantitative census data of the 58 most abundant species of benthic foraminifera produced six species associations within three higher level clusters: (1) calcareous species most abundant at mid-bathyal to outer shelf depths (<1000 m); (2) calcareous species most abundant at mid-bathyal and greater depths (>600 m); (3) agglutinated species mostly occurring at deep abyssal depths (>3000 m). A detrended correspondence analysis ordination plot exhibits a strong relationship between these species associations and bathymetry. This is manifest in the bathymetric ranges of the relative abundance peaks of many of the common benthic species (e.g., Abditodentrix pseudothalmanni 500-2800 m, Bolivina robusta 200-650 m, Bulimina marginata f. marginata 20-600 m, B. marginata f. aculeata 400-3000 m, Cassidulina norvangi 1000-4500 m, Epistominella exigua 1000-4700 m, and Trifarina angulosa 10-650 m), which should prove useful in paleobathymetric estimates. The upper depth limits of 28 benthic foraminiferal species (e.g., Fursenkoina complanata 200 m, Bulimina truncana 450 m, Melonis affinis 550 m, Eggerella bradyi 750 m, and Cassidulina norvangi 1000 m) have potential to improve the precision of paleobathymetric estimates based initially on the total faunal composition. The planktic percentage of foraminiferal tests increases from outer shelf to upper abyssal depths followed by a rapid decline within the foraminiferal lysocline (below c. 3600 m). A planktic percentage <50% is suggestive of shelf depths, and >50% is suggestive of bathyal or abyssal depths above the CCD. In the abyssal zone there is dramatic taphonomic loss of most agglutinated tests (except some textulariids) at burial depths of 0.1-0.2 m, which negates the potential usefulness of these taxa in paleobathymetric assessments.

Late Quaternary benthic foraminiferal record of the Bounty Trough, east of New Zealand

The Bounty Trough, east of New Zealand, lies along the southeastern edge of the present-day Subtropical Front (STF), and is a major conduit via the Bounty Channel, for terrigenous sediment supply from the uplifted Southern Alps to the abyssal Bounty Fan. Census data on 65 benthic foraminiferal faunas (>63 µm) from upper bathyal (ODP 1119), lower bathyal (DSDP 594) and abyssal (ODP 1122) sequences, test and refine existing models for the paleoceanographic and sedimentary history of the trough through the last 150 ka (marine isotope stages, MIS 6-1). Cluster analysis allows recognition of six species groups, whose distribution patterns coincide with bathymetry, the climate cycles and displaced turbidite beds. Detrended canonical correspondence analysis and comparisons with modern faunal patterns suggest that the groups are most strongly influenced by food supply (organic carbon flux), and to a lesser extent by bottom water oxygen and factors relating to sediment type. Major faunal changes at upper bathyal depths (1119) probably resulted from cycles of counter-intuitive seaward-landward migrations of the Southland Front (SF) (north-south sector of the STF). Benthic foraminiferal changes suggest that lower nutrient, cool Subantarctic Surface Water (SAW) was overhead in warm intervals, and higher nutrient-bearing, warm neritic Subtropical Surface Water (STW) was overhead in cold intervals. At lower bathyal depths (594), foraminiferal changes indicate increased glacial productivity and lowered bottom oxygen, attributed to increased upwelling and inflow of cold, nutrient-rich, Antarctic Intermediate Water (AAIW) and shallowing of the oxygen-minimum zone (upper Circum Polar Deep Water, CPDW). The observed cyclical benthic foraminiferal changes are not a result of associations migrating up and down the slope, as glacial faunas (dominated by Globocassidulina canalisuturata and Eilohedra levicula at upper and lower bathyal depths, respectively) are markedly different from those currently living in the Bounty Trough. On the abyssal Bounty Fan (1122), faunal changes correlate most strongly with grain size, and are attributed to varying amounts of mixing of displaced and in-situ faunas. Most of the displaced foraminifera in turbiditic sand beds are sourced from mid-outer shelf depths at the head of the Bounty Channel. Turbidity currents were more prevalent during, but not restricted to, glacial intervals.