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.

Geochemical composition of gabbroic rocks from ODP Hole 179-1105A, southwest Indian Ridge

Oxide-free olivine gabbro and gabbro, and oxide olivine gabbro and gabbro make up the bulk of the gabbroic suite recovered from Ocean Drilling Program (ODP) Leg 179 Hole 1105A, which lies 1.2 km away from Hole 735B on the eastern transverse ridge of the Atlantis II Fracture Zone, Southwest Indian Ridge. The rocks recovered during Leg 179 show striking similarities to rocks recovered from the uppermost 500 m of Hole 735B during ODP Leg 118. The rocks of the Atlantis platform were likely unroofed as part of the footwall block of a large detachment fault on the inside corner of the intersection of the Southwest Indian Ridge and the Atlantis II Transform at ~11.5 Ma. We analyzed the lithologic, geochemical, and structural stratigraphy of the section. Downhole lithologic variation allowed division of the core into 141 lithologic intervals and 4 main units subdivided on the basis of predominance of oxide gabbroic vs. oxide-free gabbroic rocks. Detailed analyses of whole-rock chemistry, mineral chemistry, microstructure, and modes of 147 samples are presented and clearly show that the gabbroic rocks are of cumulate origin. These studies also indicate that geochemistry results correlate well with downhole magnetic susceptibility and Formation MicroScanner (FMS) resistivity measurements and images. FMS images show rocks with a well-layered structure and significant numbers of mappable layer contacts or compositional contrasts. Downhole cryptic mineral and whole-rock chemical variations depict both "normal" and inverse fine-scale variations on a scale of 10 m to <2 m with significant compositional variation over a short distance within the 143-m section sampled. A Mg# shift in whole-rock or Fo contents of olivine of as much as 20-30 units over a few meters of section is not atypical of the extreme variation in downhole plots. The products of the earliest stages of basaltic differentiation are not represented by any cumulates, as the maximum Fo content was Fo78. Similarly, the extent of fractionation represented by the gabbroic rocks and scarce granophyres in the section is much greater than that represented in the Atlantis II basalts. The abundance of oxide gabbros is similar to that in Hole 735B, Unit IV, which is tentatively correlated as a similar unit or facies with the oxide gabbroic units of Hole 1105A. Oxide phases are generally present in the most fractionated gabbroic rocks and lacking in more primitive gabbroic rocks, and there is a definite progression of oxide abundance as, for example, the Mg# of clinopyroxene falls below 73-75. Coprecipitation of oxide at such early Mg#s cannot be modeled by perfect fractional crystallization. In situ boundary layer fractionation may offer a more plausible explanation for the complex juxtaposition of oxide- and nonoxide-bearing more primitive gabbroic rocks. The geochemical signal may, in part, be disrupted by the presence of mylonitic shear zones, which strike east-west and dip both to the south and north, but predominantly to the south away from the northern rift valley where they formed. Downhole deformation textures indicate increasing average strain and crystal-plastic deformation in units that contain oxides. Oxide-rich zones may represent zones of rheologic weakness in the cumulate section along which mylonitic and foliated gabbroic shear zones nucleate in the solid state at high temperature, or the oxide may be a symptom of former melt-rich zones and hypersolidus flow, as predicted during study of Hole 735B.

Spatially explicit methane inventory for Switzerland with 500 m x 500 m resolution for 2011 (NetCDF 25 MB)

We present the first high-resolution (500 m × 500 m) gridded methane (CH4) emission inventory for Switzerland, which integrates the national emission totals reported to the United Nations Framework Convention on Climate Change (UNFCCC) and recent CH4 flux studies conducted by research groups across Switzerland. In addition to anthropogenic emissions, we also include natural and semi-natural CH4 fluxes, i.e., emissions from lakes and reservoirs, wetlands, wild animals as well as uptake by forest soils. National CH4 emissions were disaggregated using detailed geostatistical information on source locations and their spatial extent and process- or area-specific emission factors. In Switzerland, the highest CH4 emissions in 2011 originated from the agricultural sector (150 Gg CH4/yr), mainly produced by ruminants and manure management, followed by emissions from waste management (15 Gg CH4/yr) mainly from landfills and the energy sector (12 Gg CH4/yr), which was dominated by emissions from natural gas distribution. Compared to the anthropogenic sources, emissions from natural and semi-natural sources were relatively small (6 Gg CH4/yr), making up only 3 % of the total emissions in Switzerland. CH4 fluxes from agricultural soils were estimated to be not significantly different from zero (between -1.5 and 0 Gg CH4/yr), while forest soils are a CH4 sink (approx. -2.8 Gg CH4/yr), partially offsetting other natural emissions. Estimates of uncertainties are provided for the different sources, including an estimate of spatial disaggregation errors deduced from a comparison with a global (EDGAR v4.2) and a European CH4 inventory (TNO/MACC). This new spatially-explicit emission inventory for Switzerland will provide valuable input for regional scale atmospheric modeling and inverse source estimation.

Major ion chemistry of snow cores along a transect in central Dronning Maud Land, Antarctica

Among the large variety of particulates in the atmosphere, calcic mineral dust particles have highly reactive surfaces that undergo heterogeneous reactions with nitrogen oxides contiguously. The association between Ca2+, an important proxy indicator of mineral dust and NO3-, a dominant anion in the Antarctic snow pack was analysed. A total of 41 snow cores (~ 1 m each) that represent snow deposited during 2008-2009 were studied along coastal-inland transects from two different regions - the Princess Elizabeth Land (PEL) and central Dronning Maud Land (cDML) in East Antarctica. Correlation statistics showed a strong association (at 99 % significance level) between NO3- and Ca2+ at the near-coastal sections of both PEL (r = 0.72) and cDML (r = 0.76) transects. Similarly, a strong association between these ions was also observed in snow deposits at the inland sections of PEL (r = 0.8) and cDML (r = 0.85). Such systematic associations between Ca2+ and NO3- is attributed to the interaction between calcic mineral dust and nitrogen oxides in the atmosphere, leading to the possible formation of calcium nitrate (Ca(NO3)2). Forward and back trajectory analyses using HYSPLIT model v. 4 revealed that Southern South America (SSA) was an important dust emitting source to the study region, aided by the westerlies. Particle size distribution showed that over 90 % of the dust was in the range < 4 µm, indicating that these dust particles reached the Antarctic region via long range transport from the SSA region. We propose that the association between Ca2+ and NO3- occurs during the long range transport due to the formation of Ca(NO3)2. The Ca(NO3)2 thus formed in the atmosphere undergo deposition over Antarctica under the influence of anticyclonic polar easterlies. However, influence of local dust sources from the nunataks in cDML evidently mask such association in the mountainous region. The study indicates that the input of dust-bound NO3- may contribute a significant fraction of the total NO3- deposited in Antarctic snow.

Mean dynamic topography and oceanographic parameters estimated from an inverse model and satellite geodesy, with link to model result in one single NetCDF file (392 MB), including the inverse data error covariance

Geostrophic surface velocities can be derived from the gradients of the mean dynamic topography-the difference between the mean sea surface and the geoid. Therefore, independently observed mean dynamic topography data are valuable input parameters and constraints for ocean circulation models. For a successful fit to observational dynamic topography data, not only the mean dynamic topography on the particular ocean model grid is required, but also information about its inverse covariance matrix. The calculation of the mean dynamic topography from satellite-based gravity field models and altimetric sea surface height measurements, however, is not straightforward. For this purpose, we previously developed an integrated approach to combining these two different observation groups in a consistent way without using the common filter approaches (Becker et al. in J Geodyn 59(60):99-110, 2012, doi:10.1016/j.jog.2011.07.0069; Becker in Konsistente Kombination von Schwerefeld, Altimetrie und hydrographischen Daten zur Modellierung der dynamischen Ozeantopographie, 2012, http://nbn-resolving.de/nbn:de:hbz:5n-29199). Within this combination method, the full spectral range of the observations is considered. Further, it allows the direct determination of the normal equations (i.e., the inverse of the error covariance matrix) of the mean dynamic topography on arbitrary grids, which is one of the requirements for ocean data assimilation. In this paper, we report progress through selection and improved processing of altimetric data sets. We focus on the preprocessing steps of along-track altimetry data from Jason-1 and Envisat to obtain a mean sea surface profile. During this procedure, a rigorous variance propagation is accomplished, so that, for the first time, the full covariance matrix of the mean sea surface is available. The combination of the mean profile and a combined GRACE/GOCE gravity field model yields a mean dynamic topography model for the North Atlantic Ocean that is characterized by a defined set of assumptions. We show that including the geodetically derived mean dynamic topography with the full error structure in a 3D stationary inverse ocean model improves modeled oceanographic features over previous estimates.