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.

Georeferenced standardized apparent wave power along the coastline of Cabrera National Park for analysis of essential habitats of littoral species

Wind- induced exposure is one of the major forces shaping the geomorphology and biota in coastal areas. The effect of wave exposure on littoral biota is well known in marine environments (Ekebon et al., 2003; Burrows et al., 2008). In the Cabrera Archipelago National Park wave exposure has demostrated to have an effect on the spatial distribution of different stages of E.marginatus (Alvarez et al., 2010). Standarized average wave exposures during 2008 along the Cabrera Archipelago National park coast line were calculated to be applied in studies of littoral species distribution within the archipelago. Average wave exposure (or apparent wave power) was calculated for points located 50 m equidistant on the coastline following the EXA methodology (EXposure estimates for fragmented Archipelagos) (Ekebon et al., 2003). The average wave exposures were standardized from 1 to 100 (minimum and maximum in the area), showing coastal areas with different levels of mea wave exposure during the year. Input wind data (direction and intensity) from 2008 was registered at the Cabrera mooring located north of Cabrera Archipelago. Data were provided by IMEDEA (CSIC-UIB, TMMOS http://www.imedea.uib-csic.es/tmoos/boyas/). This cartography has been developed under the framework of the project EPIMHAR, funded by the National Park's Network (Spanish Ministry of Environment, Maritime and Rural Affairs, reference: 012/2007 ). Part of this work has been developed under the research programs funded by "Fons de Garantia Agrària i Pesquera de les Illes Balears (FOGAIBA)".

Core scans and lithologic analysis of Subglacial Lake Whillans sediment cores

The hydrologic system beneath the Antarctic Ice Sheet is thought to influence both the dynamics and distribution of fast flowing ice streams, which discharge most of the ice lost by the ice sheet. Despite considerable interest in understanding this subglacial network and its affect on ice flow, in situ observations from the ice sheet bed are exceedingly rare. Here we describe the first sediment cores recovered from an active subglacial lake. The lake, known as Subglacial Lake Whillans, is part of a broader, dynamic hydrologic network beneath the Whillans Ice Stream in West Antarctica. Even though "floods" pass through the lake, the lake floor shows no evidence of erosion or deposition by flowing water. By inference, these floods must have insufficient energy to erode or transport significant volumes of sediment coarser than silt. Consequently, water flow beneath the region is probably incapable of incising continuous channels into the bed and instead follows preexisting subglacial topography and surface slope. Sediment on the lake floor consists of till deposited during intermittent grounding of the ice stream following flood events. The fabrics within the till are weaker than those thought to develop in thick deforming beds suggesting subglacial sediment fluxes across the ice plain are currently low and unlikely to have a large stabilizing effect on the ice stream's grounding zone.