Relative abundance of dinoflagellate cysts in surface sediments and dinoflagellate counts from three sediment cores from the Atlantic sector of the Southern Ocean

In this study we investigate the potential of organic-walled dinoflagellate cysts (dinocysts) as tools for quantifying past sea-surface temperatures (SST) in the Southern Ocean. For this purpose, a dinocyst reference dataset has been formed, based on 138 surface sediment samples from different circum-Antarctic environments. The dinocyst assemblages of these samples are composed of phototrophic (gonyaulacoid) and heterotrophic (protoperidinioid) species that provide a broad spectrum of palaeoenvironmental information. The relationship between the environmental parameters in the upper water column and the dinocyst distribution patterns of individual species has been established using the statistical method of Canonical Correspondence Analysis (CCA). Among the variables tested, summer SST appeared to correspond to the maximum variance represented in the dataset. To establish quantitative summer SST reconstructions, a Modern Analogue Technique (MAT) has been performed on data from three Late Quaternary dinocyst records recovered from locations adjacent to prominent oceanic fronts in the Atlantic sector of the Southern Ocean. These dinocyst time series exhibit periodic changes in the dinocyst assemblage during the last two glacial/interglacial-cycles. During glacial conditions the relative abundance of protoperidinioid cysts was highest, whereas interglacial conditions are characterised by generally lower cyst concentrations and increased relative abundance of gonyaulacoid cysts. The MAT palaeotemperature estimates show trends in summer SST changes following the global oxygen isotope signal and a strong correlation with past temperatures of the last 140,000 years based on other proxies. However, by comparing the dinocyst results to quantitative estimates of summer SSTs based on diatoms, radiolarians and foraminifer-derived stable isotope records it can be shown that in several core intervals the dinocyst-based summer SSTs appeared to be extremely high. In these intervals the dinocyst record seems to be highly influenced by selective degradation, leading to unusual temperature ranges and to unrealistic palaeotemperatures. We used the selective degradation index (kt-index) to determine those intervals that have been biased by selective degradation in order to correct the palaeotemperature estimates. We show that after correction the dinocyst based SSTs correspond reasonably well with other palaeotemperature estimates for this region, supporting the great potential of dinoflagellate cysts as a basis for quantitative palaeoenvironmental studies.

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.

(Table 1) Surface water chemistry of northern Victoria Land lakes

Concentrations of major ions, silicate and nutrients (total N and P) were measured in samples of surface water from 28 lakes in ice-free areas of northern Victoria Land (East Antarctica). Sixteen lakes were sampled during austral summers 2001/02, 2003/04, 2004/05 and 2005/06 to assess temporal variation in water chemistry. Although samples showed a wide range in ion concentrations, their composition mainly reflected that of seawater. In general, as the distance from the sea increased, the input of elements from the marine environment (through aerosols and seabirds) decreased and there was an increase in nitrate and sulfate concentrations. Antarctic lakes lack outflows and during the austral summer the melting and/or ablation of ice cover, water evaporation and leaching processes in dry soils determine a progressive increase in water ion concentrations. During the five-year monitoring survey, no statistically significant variation in the water chemistry were detected, except for a slight (hardly significant) increase in TN concentrations. However, Canonical Correspondence Analysis (CCA) indicated that other factors besides distance from the sea, the presence of nesting seabirds, the sampling time and percentage of ice cover affect the composition of water in Antarctic cold desert environments.

Environmental variables and multivariate seabed classification maps via k means clustering of Potter Cove, Antarctica, link to input and result files in GeoTIFF format

This study subdivides the Potter Cove, King George Island, 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 includes in total 42 different environmental variables, interpolated based on samples taken during Australian summer seasons 2010/2011 and 2011/2012. The statistical errors of several interpolation methods (e.g. IDW, Indicator, Ordinary and Co-Kriging) with changing settings have been compared and the most reasonable method has been applied. 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 and 4, 7, 10 as well as 12 were identified as reasonable numbers for clustering the Potter Cove. Especially the results of 10 and 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.

Palynology and dinoflagellate counts of surface sediments, and sea surface conditions, of sites in the North Pacific

Palynological analyses were performed on 53 surface sediment samples from the North Pacific Ocean, including the Bering and Okhotsk Seas (37-64°N, 144°E-148°W), in order to document the relationships between the dinocyst distribution and sea-surface conditions (temperatures, salinities, primary productivity and sea-ice cover). Samples are characterized by concentrations ranging from 18 to 143816 cysts/cm**3 and the occurrence of 32 species. A canonical correspondence analysis (CCA) was carried out to determine the relationship between environmental variables and the distribution of dinocyst taxa. The first and second axes represent, respectively, 47% and 17.8% of the canonical variance. Axis 1 is positively correlated with all parameters except to the sea-ice and primary productivity in August, which are on the negative side. Results indicate that the composition of dinocyst assemblages is mostly controlled by temperature and that all environmental variables are correlated together. The CCA distinguishes 3 groups of dinocysts: the heterotrophic taxa, the genera Impagidinium and Spiniferites as well as the cyst of Pentapharsodinium dalei and Operculodinium centrocarpum. Five assemblage zones can be distinguished: 1) the Okhotsk Sea zone, which is associated to temperate and eutrophic conditions, seasonal upwellings and Amur River discharges. It is characterized by the dominance of O. centrocarpum, Brigantedinium spp. and Islandinium minutum; 2) the Western Subarctic Gyre zone with subpolar and mesotrophic conditions due to the Kamchatka Current and Alaska Stream inflows. Assemblages are dominated by Nematosphaeropsis labyrinthus, Pyxidinopsis reticulata and Brigantedinium spp.; 3) the Bering Sea zone, depicting a subpolar environment, influenced by seasonal upwellings and inputs from the Anadyr and Yukon Rivers. It is characterized by the dominance of I. minutum and Brigantedinium spp.; 4) the Alaska Gyre zone with temperate conditions and nutrient-enriched surface waters, which is dominated by N. labyrinthus and Brigantedinium spp. and 5) the Kuroshio Extension-North Pacific-Subarctic Current zone characterized by a subtropical and oligotrophic environment, which is dominated by O. centrocarpum, N. labyrinthus and warm taxa of the genus Impagidinium. Transfer functions were tested using the modern analog technique (MAT) on the North Pacific Ocean (= 359 sites) and the entire Northern Hemisphere databases ( = 1419 sites). Results confirm that the updated Northern Hemisphere database is suitable for further paleoenvironmental reconstructions, and the best results are obtained for temperatures with an accuracy of +/-1.7 °C.

Diatom abundance in surface sediments of the northern North Pacific

In order to map the modern distribution of diatoms and to establish a reliable reference data set for paleoenvironmental reconstruction in the northern North Pacific, a new data set including the relative abundance of diatom species preserved in a total of 422 surface sediments was generated, which covers a broad range of environmental variables characteristic of the subarctic North Pacific, the Sea of Okhotsk and the Bering Sea between 30° and 70°N. The biogeographic distribution patterns as well as the preferences in sea surface temperature of 38 diatom species and species groups are documented. A Q-mode factor analysis yields a three-factor model representing assemblages associated with the Arctic, Subarctic and Subtropical water mass, indicating a close relationship between the diatom composition and the sea surface temperatures. The relative abundance pattern of 38 diatom species and species groups was statistically compared with nine environmental variables, i.e. the summer sea surface temperature and salinity, annual surface nutrient concentration (nitrate, phosphate, silicate), summer and winter mixed layer depth and summer and winter sea ice concentrations. Canonical Correspondence Analysis (CCA) indicates 32 species and species groups have strong correspondence with the pattern of summer sea surface temperature. In addition, the total diatom flux data compiled from ten sediment traps reveal that the seasonal signals preserved in the surface sediments are mostly from spring through autumn. This close relationship between diatom composition and the summer sea surface temperature will be useful in deriving a transfer function in the subarctic North Pacific for the quantitative paleoceanographic and paleoenvironmental studies. The relative abundance of the sea-ice indicator diatoms Fragilariopsis cylindrus and F. oceanica of >20% in the diatom composition is used to represent the winter sea ice edge in the Bering Sea. The northern boundary of the distribution of F. doliolus in the open ocean is suggested to be an indicator of the Subarctic Front, while the abundance of Chaetoceros resting spores may indicate iron input from nearby continents and shelves and induced productivity events in the study area.

Ostracode relative abundances and geochemistry in surface sediments from Lake Ohrid

Understanding the ecology of bioindicators such as ostracods is essential in order to reconstruct past environmental and climate change from analysis of fossil assemblages preserved in lake sediment cores. Knowledge of the ecology of ancient Lake Ohrid's ostracod fauna is very limited and open to debate. In advance of the Ohrid ICDP-Drilling project, which has potential to generate high-resolution long-term paleoenvironmental data of global importance in paleoclimate research, we sampled Lake Ohrid and a wide range of habitat types in its surroundings to assess 1) the composition of ostracod assemblages in lakes, springs, streams, and short-lived seasonal water bodies, 2) the geographical distribution of ostracods, and 3) the ecological characteristics of individual ostracod species. In total, 40 species were collected alive, and seven species were preserved as valves and empty carapaces. Of the 40 ostracod species, twelve were endemic to Lake Ohrid. The most common genus in the lake was Candona, represented by 13 living species, followed by Paralimnocythere, represented by five living species. The most frequent species was Cypria obliqua. Species with distinct distributions included Heterocypris incongruens, Candonopsis kingsleii, and Cypria lacustris. The most common species in shallow, flooded areas was H. incongruens, and the most prominent species in ditches was C. kingsleii. C. lacustris was widely distributed in channels, springs, lakes, and rivers. Statistical analyses were performed on a "Lake Ohrid" dataset, comprising the subset of samples from Lake Ohrid alone, and an "entire" dataset comprising all samples collected. The unweighted pair group mean average (UPGMA) clustering was mainly controlled by species-specific depth preferences. Canonical Correspondence Analysis (CCA) with forward selection identified water depth, water temperature, and pH as variables that best explained the ostracod distribution in Lake Ohrid. The lack of significance of conductivity and dissolved oxygen in CCA of Ohrid data highlight the uniformity across the lake of the well-mixed waters. In the entire area, CCA revealed that ostracod distribution was best explained by water depth, salinity, conductivity, pH, and dissolved oxygen. Salinity was probably selected by CCA due to the presence of Eucypris virens and Bradleystrandesia reticulata in short-lived seasonal water bodies. Water depth is an important, although indirect, influence on ostracod species distribution which is probably associated with other factors such as sediment texture and food supply. Some species appeared to be indicators for multiple environmental variables, such as lake level and water temperature.

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.

Benthic foraminifera asemblages and isotopic composition of surface sediment samples off Morocco

The influence of different primary productivity regimes on live (Rose Bengal stained) and dead benthic foraminiferal distribution, as well as on the stable carbon isotopic composition of foraminiferal tests, was investigated in sediment surface samples (0-1 cm) from the upwelling region off Morocco between Cape Ghir (31°N) and Cape Yubi (27°N). A combination of factor analysis, detrended correspondence analysis (DCA) and canonical correspondence analysis (CCA) was applied to the benthic foraminiferal data sets. Five major assemblages for both the live and dead fauna were revealed by factor analysis. In the cape regions organic matter fluxes are enhanced by high chlorophyll-a concentrations in the overlying surface waters. Here, benthic foraminiferal faunas are characterized by identical live and dead assemblages, high standing stocks, and low species delta13C values, indicating constant year-round high productivity. Bulimina marginata dominates the unique fauna at the shallowest station off Cape Ghir indicating highest chlorophyll-a concentrations. Off both capes, the succession of the Bulimina aculeata/Uvigerina mediterranea assemblage, the Sphaeroidina bulloides/Gavelinopsis translucens assemblage, and the Hoeglundina elegans assemblage from the shelf to the deep sea reflects the decrease in chlorophyll-a concentrations, hence the export flux. In contrast, the area between the capes is characterized by differently composed live and dead assemblages, low standing stocks, and less depleted delta13C values, thus reflecting low primary productivity. High foraminiferal numbers of Epistominella exigua, Eponides pusillus, and Globocassidulina subglobosa in the dead fauna indicate a seasonally varying primary productivity signal. Significantly lower mean delta13C values were recorded in Bulimina mexicana, Cibicidoides kullenbergi, H. elegans, U. mediterranea and Uvigerina peregrina. Cibicidoides wuellerstorfi is a faithful recorder of bottom water delta13C in the Canary Islands regions. The mean delta13C signal of this species is not significantly influenced by constant high organic matter fluxes. The species-specific offset between live and dead specimens is the same.

(Appendix) Analysis of benthic foraminifera from ODP Leg 181 and DSDP Leg 90

Canonical correspondence analysis indicates that the distribution of Neogene benthic foraminiferal faunas (>63 µm) in seven DSDP and ODP sites (500-4500 m water depth) east of New Zealand (38-51°S, 170°E-170°W) is most strongly influenced by depth (water mass stratification), and secondly by age (palaeoceanographic changes influencing faunal composition and biotic evolution). Stratigraphic faunal changes are interpretted in terms of the pulsed sequential development of southern, and later northern, polar glaciation and consequent cooling of bottom waters, increased vertical and lateral stratification of ocean water masses, and increased overall and seasonal surface water productivity. Oligocene initiation of the Antarctic Circumpolar Current and Deep Western Boundary Current (DWBC), flowing northwards past New Zealand, resulted in extensive hiatuses throughout the Southwest Pacific, some extending through into the Miocene. Planktic foraminiferal fragmentation index values indicate that carbonate dissolution was significant at abyssal depths throughout most of the Neogene, peaking at upper abyssal depths in the late Miocene (11-7 Ma), with the lysocline progressively deepened thereafter. Miocene abyssal faunas are dominated by Globocassidulina subglobosa and Oridorsalis umbonatus, with increasing Epistominella exigua after 16 Ma at upper abyssal depths. Peak abundances of Epistominella umbonifera indicate increased input of cold Southern Component Water to the DWBC at 7-6 Ma. Faunal association changes imply establishment of the modern Oxygen Minimum Zone (upper Circumpolar Deep Water) in the latest Miocene. Significant latitudinal differences between the benthic foraminiferal faunas at lower bathyal depths indicate the existence of an oceanic front along the Chatham Rise (location of present Subtropical Front), since the early late Miocene at least, with more pulsed productivity (higher E. exigua) along the south side. Modern Antarctic Intermediate Water faunal associations were established north of the Chatham Rise at 10-9 Ma, and south of it at 3-1.5 Ma. Middle-upper bathyal faunas on the Campbell Plateau are dominated by reticulate bolivinids during the early and middle Miocene, indicative of sustained productivity above relatively sluggish, suboxic bottom waters. Faunal changes and hiatuses indicate increased current vigour over the Campbell Plateau from the latest Miocene on. Surface water productivity (food supply) appears to have increased in three steps (at times of enhanced global cooling) marked by substantially increased relative abundance of: (1) Abditodentrix pseudothalmanni, Alabaminella weddellensis, Cassidulina norvangi (16-15 Ma, increased pulsed productivity); (2) Bulimina marginata f. aculeata, Nonionella auris, Trifarina angulosa, Uvigerina peregrina (3-1.5 Ma, increased overall productivity); and (3) Cassidulina carinata (1-0.5 Ma, increased overall productivity). Three intervals of deep-sea benthic foraminiferal taxonomic turnover are recognised (16-15, 11.5-10, 2-0.5 Ma) corresponding to intervals of enhanced global cooling and possible productivity changes. The late Pliocene-middle Pleistocene extinction, associated with increasing Northern Hemisphere glaciation, culminating in the middle Pleistocene climatic transition, was more significant in the study area than the earlier Neogene turnovers.

Late Pliocene dinoflagellate cyst and diatom analysis from a high resolution sequence in DSDP Site 594, Chatham Rise, south west Pacific

A high resolution mixed carbonate and siliciclastic sequence from DSDP Site 594 contains a detailed record of climate change in the late Pliocene. The sequence can be accurately dated by the LAD of Nitzschia weaveri, the LAD of Thalassiosira insigna, the LAD of T. vulnifica and the LAD of T. kolbei diatom datums. Carbonate content and delta18O signatures provide added resolution and place the sequence between isotope stage 100 and 92. The sequence contains well-preserved and diverse dinoflagellate cyst floras. Use of principal component (PCA) and canonical correspondence analyses (CCA) identifies changes in the assemblages that principally reflect warming and cooling trends. Species association with warmer climates included Impagidinium patulum, I. paradoxum and I. sp. cf. paradoxum while those from cooler climates include Invertecysta tabulata and I. velorum. CCA is shown to be a valuable method of determining the past environmental preferences of extinct species such as I. tabulata.

Artificial oxygen fluxes measured by the eddy correlation method using stirring-sensitive oxygen microsensor and oxygen optodes in a flume experiment

In the last decade, the aquatic eddy correlation (EC) technique has proven to be a powerful approach for non-invasive measurements of oxygen fluxes across the sediment water interface. Fundamental to the EC approach is the correlation of turbulent velocity and oxygen concentration fluctuations measured with high frequencies in the same sampling volume. Oxygen concentrations are commonly measured with fast responding electrochemical microsensors. However, due to their own oxygen consumption, electrochemical microsensors are sensitive to changes of the diffusive boundary layer surrounding the probe and thus to changes in the ambient flow velocity. The so-called stirring sensitivity of microsensors constitutes an inherent correlation of flow velocity and oxygen sensing and thus an artificial flux which can confound the benthic flux determination. To assess the artificial flux we measured the correlation between the turbulent flow velocity and the signal of oxygen microsensors in a sealed annular flume without any oxygen sinks and sources. Experiments revealed significant correlations, even for sensors designed to have low stirring sensitivities of ~0.7%. The artificial fluxes depended on ambient flow conditions and, counter intuitively, increased at higher velocities because of the nonlinear contribution of turbulent velocity fluctuations. The measured artificial fluxes ranged from 2 - 70 mmol m**-2 d**-1 for weak and very strong turbulent flow, respectively. Further, the stirring sensitivity depended on the sensor orientation towards the flow. Optical microsensors (optodes) that should not exhibit a stirring sensitivity were tested in parallel and did not show any significant correlation between O2 signals and turbulent flow. In conclusion, EC data obtained with electrochemical sensors can be affected by artificial flux and we recommend using optical microsensors in future EC-studies. Flume experiments were conducted in February 2013 at the Institute for Environmental Sciences, University of Koblenz-Landau Landau. Experiments were performed in a closed oval-shaped acrylic glass flume with cross-sectional width of 4 cm and height of 10 cm and total length of 54 cm. The fluid flow was induced by a propeller driven by a motor and mean flow velocities of up to 20 cm s-1 were generated by applying voltages between 0 V and 4 V DC. The flume was completely sealed with an acrylic glass cover. Oxygen sensors were inserted through rubber seal fittings and allowed positioning the sensors with inclinations to the main flow direction of ~60°, ~95° and ~135°. A Clark type electrochemical O2 microsensor with a low stirring sensitivity (0.7%) was tested and a fast-responding needle-type O2 optode (PyroScience GmbH, Germany) was used as reference as optodes should not be stirring sensitive. Instantaneous three-dimensional flow velocities were measured at 7.4 Hz using stereoscopic particle image velocimetry (PIV). The velocity at the sensor tip was extracted. The correlation of the fluctuating O2 sensor signals and the fluctuating velocities was quantified with a cross-correlation analysis. A significant cross-correlation is equivalent to a significant artificial flux. For a total of 18 experiments the flow velocity was adjusted between 1.7 and 19.2 cm s**-1, and 3 different orientations of the electrochemical sensor were tested with inclination angles of ~60°, ~95° and ~135° with respect to the main flow direction. In experiments 16-18, wavelike flow was induced, whereas in all other experiments the motor was driven by constant voltages. In 7 experiments, O2 was additionally measured by optodes. Although performed simultaneously with the electrochemical sensor, optode measurements are listed as separate experiments (denoted by the attached 'op' in the filename), because the velocity time series was extracted at the optode tip, located at a different position in the flume.

Metadata und statistic analysis of archael and bacterial sequences originating from sediments of the Håkon Mosby mud volcano (Z1-Z3, all habitats)

DNA extraction was carried out as described on the MICROBIS project pages (http://icomm.mbl.edu/microbis ) using a commercially available extraction kit. We amplified the hypervariable regions V4-V6 of archaeal and bacterial 16S rRNA genes using PCR and several sets of forward and reverse primers (http://vamps.mbl.edu/resources/primers.php). Massively parallel tag sequencing of the PCR products was carried out on a 454 Life Sciences GS FLX sequencer at Marine Biological Laboratory, Woods Hole, MA, following the same experimental conditions for all samples. Sequence reads were submitted to a rigorous quality control procedure based on mothur v30 (doi:10.1128/AEM.01541-09) including denoising of the flow grams using an algorithm based on PyroNoise (doi:10.1038/nmeth.1361), removal of PCR errors and a chimera check using uchime (doi:10.1093/bioinformatics/btr381). The reads were taxonomically assigned according to the SILVA taxonomy (SSURef v119, 07-2014; doi:10.1093/nar/gks1219) implemented in mothur and clustered at 98% ribosomal RNA gene V4-V6 sequence identity. V4-V6 amplicon sequence abundance tables were standardized to account for unequal sampling effort using 1000 (Archaea) and 2300 (Bacteria) randomly chosen sequences without replacement using mothur and then used to calculate inverse Simpson diversity indices and Chao1 richness (doi:10.2307/4615964). Bray-Curtis dissimilarities (doi:10.2307/1942268) between all samples were calculated and used for 2-dimensional non metric multidimensional scaling (NMDS) ordinations with 20 random starts (doi:10.1007/BF02289694). Stress values below 0.2 indicated that the multidimensional dataset was well represented by the 2D ordination. NMDS ordinations were compared and tested using Procrustes correlation analysis (doi:10.1007/BF02291478). All analyses were carried out with the R statistical environment and the packages vegan (available at: http://cran.r-project.org/package=vegan), labdsv (available at: http://cran.r-project.org/package=labdsv), as well as with custom R scripts. Operational taxonomic units at 98% sequence identity (OTU0.03) that occurred only once in the whole dataset were termed absolute single sequence OTUs (SSOabs; doi:10.1038/ismej.2011.132). OTU0.03 sequences that occurred only once in at least one sample, but may occur more often in other samples were termed relative single sequence OTUs (SSOrel). SSOrel are particularly interesting for community ecology, since they comprise rare organisms that might become abundant when conditions change.16S rRNA amplicons and metagenomic reads have been stored in the sequence read archive under SRA project accession number SRP042162.