Methane and sulphate consumption in different habitats of the Håkon Mosby Mud Volcano (HMMV)

The Hakon Mosby Mud Volcano is a highly active methane seep hosting different chemosynthetic communities such as thiotrophic bacterial mats and siboglinid tubeworm assemblages. This study focuses on in situ measurements of methane fluxes to and from these different habitats, in comparison to benthic methane and oxygen consumption rates. By quantifying in situ oxygen, methane, and sulfide fluxes in different habitats, a spatial budget covering areas of 10-1000 -m diameter was established. The range of dissolved methane efflux (770-2 mmol m-2 d-1) from the center to the outer rim was associated with a decrease in temperature gradients from 46°C to < 1°C m-1, indicating that spatial variations in fluid flow control the distribution of benthic habitats and activities. Accordingly, total oxygen uptake (TOU) varied between the different habitats by one order of magnitude from 15 mmol m-2 d-1 to 161 mmol m-2 d-1. High fluid flow rates appeared to suppress benthic activities by limiting the availability of electron acceptors. Accordingly, the highest TOU was associated with the lowest fluid flow and methane efflux. This was most likely due to the aerobic oxidation of methane, which may be more relevant as a sink for methane as previously considered in submarine ecosystems.

Seasonal variations in surface metabolite composition of Fucus vesiculosus and Fucus serratus sampled at outer Kiel Fjord (August 2012 - July 2013)

The chemical composition of surface associated metabolites of two Fucus species (Fucus vesiculosus and Fucus serratus) was analysed by means of gas chromatography-mass spectrometry (GC-MS) to describe temporal patterns in chemical surface composition. Method: The two perennial brown macroalgae F. vesiculosus and F. serratus were sampled monthly at Bülk, outer Kiel Fjord, Germany (54°27'21 N / 10°11'57 E) over an entire year (August 2012 - July 2013). Per month and species six non-fertile Fucus individuals were collected from mixed stands at a depth of 0.5 m under mid water level. For surface extraction approx. 50 g of the upper 5-10 cm apical thalli tips were cut off per species. The surface extraction of Fucus was performed according to the protocol of de Nys and co-workers (1998) with minor modifications (see Rickert et al. 2015). GC/EI-MS measurements were performed with a Waters GCT premier (Waters, Manchester, UK) coupled to an Agilent 6890N GC equipped with a DB-5 ms 30 m column (0.25 mm internal diameter, 0.25 mM film thickness, Agilent, USA). The inlet temperature was maintained at 250°C and samples were injected in split 10 mode. He carrier gas flow was adjusted to 1 ml min-1. Alkanes were used for referencing of retention times. For further details (GC-MS sample preparation and analysis) see the related publication (Rickert et al. submitted to PLOS ONE).

222Radon flux map for Europe in netCDF format

A high-resolution 222Radon (222Rn) flux map for Europe was developed, based on a parameterization of 222Rn production and transport in the soil. The 222Rn exhalation rate is parameterized based on soil properties, uranium content, and modelled soil moisture from two different land-surface reanalysis data sets. Spatial variations in exhalation rates are primarily determined by the uranium content of the soil, but also influenced by soil texture and local water table depth. Temporal variations are related to soil moisture variations as the molecular diffusion in the unsaturated soil zone depends on available air-filled pore space. Monthly 222Rn exhalation rates from European soils were calculated with a nominal spatial resolution of 0.083° x 0.083°. The two realizations of the 222Rn flux map, based on the different soil moisture data sets, both realistically reproduce the observed seasonality in the fluxes but yield considerable differences for absolute flux values. The mean 222Rn flux from soils in Europe is estimated to be 10 mBq/m**2/s (ERA-Interim/Land soil moisture) or 15 mBq/m**2/s (GLDAS-Noah soil moisture) for the period 2006-2010. The 222Rn flux maps for Europe are available for the application in atmospheric transport studies, e.g to evaluate the performance of atmospheric transport models.

Geochemistry on Holocene sediments from the Mediterranean

The enhanced accumulation of organic matter in Eastern Mediterranean sapropels and their unusually low d15N values have been attributed to either enhanced nutrient availability which led to elevated primary production and carbon sequestration or to enhanced organic matter preservation under anoxic conditions. In order to evaluate these two hypothesis we have determined Ba/Al ratios, amino acid composition, N and organic C concentrations and d15N in sinking particles, surface sediments, eight spatially distributed core records of the youngest sapropel S1 (10-6 ka) and older sapropels (S5, S6) from two locations. These data suggest that (i) temporal and spatial variations in d15N of sedimentary N are driven by different degrees of diagenesis at different sites rather than by changes in N-sources or primary productivity and (ii) present day TOC export production would suffice to create a sapropel like S1 under conditions of deep-water anoxia. This implies that both enhanced TOC accumulation and d15N depletion in sapropels were due to the absence of oxygen in deep waters. Thus preservation plays a major role for the accumulation of organic-rich sediments casting doubt on the need of enhanced primary production for sapropel formation.

Chemical composition of basalts and minerals from basalts of the Bouvet triple junction

This study focuses on mafic volcanic rocks from the Bouvet triple junction, which fall into six geochemically distinct groups: (1) N-MORB, the most widespread type, encountered throughout the study area. (2) Subalkaline volcanics, hawaiites and mugearites strongly enriched in lithophile elements and radiogenic isotopes and composing the Bouvet volcanic rise, and compositionally similar basalts and basaltic andesites from the Spiess Ridge, generated in a deeper, fertile mantle region. (3) Relatively weakly enriched basalts, T-MORB derived by the mixing of Type 1 and 2 melts and exposed near the axes of the Mid-Atlantic, Southwest Indian, and America-Antarctic Ridges. (4) Basalts with a degree of trace lithophile element enrichment similar to the Spiess Ridge and Bouvet Island rocks, but higher in K, P, Ti, and Cr. These occur within extensional structures: the rift valley of the Southwest Indian Ridge, grabens of the East Dislocation Zone, and the linear rise between the Spiess Ridge and Bouvet volcano. Their parental melts presumably separated from plume material that spread from the main channels and underwent fluid-involving differentiation in the mantle. (5) A volcanic suite ranging from basalt to rhyolite, characterized by low concentrations of lithophile elements, particularly TiO2, and occurring on the Shona Seamount and other compressional features within the Antarctic and South American plates near the Bouvet triple junction. Unlike Types 1 to 4, which display tholeiitic differentiation trends, this suite is calc-alkaline. Its parental melts were presumably related to the plume material as well but, subsequently, they underwent a profound differentiation involving fluids and assimilated surrounding rocks in closed magma chambers in the upper mantle. Alternatively, the Shona Seamount might be a fragment of an ancient oceanic island arc. (6) Enriched basalts, distinguished from the other enriched rock types in very high P and radiogenic isotope abundances and composing a tectonic uplift near the junction of the three rifts. It thus follows that the main factors responsible for the compositional diversity of volcanic rocks in this region include (i) mantle source heterogeneity, (ii) plume activity, (iii) an intricate geodynamic setup at the triple junction giving rise to stresses in adjacent plate areas, and (iv) the geological prehistory. The slow spreading rate and ensuing inefficient mixing of the heterogeneous mantle material result in strong spatial variations in basaltic compositions.

(Table 1) Basin area, average late-summer speedup, number of sinks and runoff of West Greenland Ice Sheet catchments

We use interferometric synthetic aperture radar observations recorded in a land-terminating sector of western Greenland to characterise the ice sheet surface hydrology and to quantify spatial variations in the seasonality of ice sheet flow. Our data reveal a non-uniform pattern of late-summer ice speedup that, in places, extends over 100 km inland. We show that the degree of late-summer speedup is positively correlated with modelled runoff within the 10 glacier catchments of our survey, and that the pattern of late-summer speedup follows that of water routed at the ice sheet surface. In late-summer, ice within the largest catchment flows on average 48% faster than during winter, whereas changes in smaller catchments are less pronounced. Our observations show that the routing of seasonal runoff at the ice sheet surface plays an important role in shaping the magnitude and extent of seasonal ice sheet speedup.

Growth rates of late Miocene corals from Crete (Greece)

Modern scleractinian corals are classical components of marine shallow warm water ecosystems. Their occurrence and diversity patterns in the geological record have been widely used to infer past climates and environmental conditions. Coral skeletal composition data reflecting the nature of the coral environment are often affected by diagenetic alteration. Ghost structures of annual growth rhythms are, however, often well preserved in the transformed skeleton. We show that these relicts represent a valuable source of information on growth conditions of fossil corals. Annual growth bands were measured in massive hemispherical Porites of late Miocene age from the island of Crete (Greece) that were found in patch reefs and level bottom associations of attached mixed clastic environments as well as isolated carbonate environments. The Miocene corals grew slowly, about 2-4 mm/yr, compatible with present-day Porites from high-latitude reefs. Slow annual growth of the Miocene corals is in good agreement with the position of Crete at the margin of the Miocene reef belt. Within a given time slice, extension rates were lowest in level bottom environments and highest in attached inshore reef systems. Because sea surface temperatures (SSTs) can be expected to be uniform within a time slice, spatial variations in extension rates must reflect local variations in light levels (low in the level bottom communities) and nutrients (high in the attached reef systems). During the late Miocene (Tortonian-early Messinian), maximum linear extension rates remained remarkably constant within seven chronostratigraphic units, and if the relationship of SSTs and annual growth rates observed for modern massive Indo-Pacific Porites spp. applies to the Neogene, minimum (winter) SSTs were 20°-21°C. Although our paleoclimatic record has a low resolution, it fits the trends revealed by global data sets. In the near future we expect this new and easy to use Porites thermometer to add important new information to our understanding of Neogene climate.

Amino acids from ODP Holes 117-722B and 177-724C

Sedimentary d15N records are valuable archives of ocean history but they are often modified during early diagenesis. Here we quantify the effect of early diagenetic enrichment on sedimentary N-isotope composition in order to obtain the pristine signal of reactive N assimilated in the euphotic zone. This is possible by using paired data of d15N and amino acid composition of sediment samples, which can be applied to estimate the degree of organic matter degradation. We determined d15N and amino acid composition in coeval sediments from Ocean Drilling Program (ODP) Hole 772 B in the central Arabian Sea and from Hole 724 C situated on the Oman Margin in the western Arabian Sea coastal upwelling area. The records span the last 130 kyr and include two glacial-interglacial cycles. These new data are used in conjunction with data available for surface sediments that cover a wide range of organic matter degradation states, and with other cores from the northern and eastern Arabian Sea to explore spatial variations in the isotopic signal. In order to reconstruct pristine N values we apply the relationship between organic matter degradation and 15N enrichment in surface sediments to correct the core records for early diagenetic enrichment. Reconstructed d15N values suggest a significant role of N2-fixation during glacial stages. An evaluation of two preservation indices based on amino acid composition (Reactivity Index, RI; Jennerjahn and Ittekkot, 1997; and the Degradation Index, DI; Dauwe et al., 1999) in both recent sediments and core samples suggests that the RI is more suitable than the DI in correcting Arabian Sea d15N records for early diagenetic enrichment.

Rare-earth elements and isotope ratios at DSDP Leg 82 Holes

We report 48 analyses of rare-earth elements (REE) and 15 143Nd/144Nd and 87Sr/86Sr analyses for basalts from the eight holes drilled during Leg 82. Discrete and distinct REE patterns and 143Nd/144Nd ratios characterize the eight holes, with little variation observed downhole except in Holes 561 and 558, thus suggesting dominantly long-term temporal and large-scale spatial variations in the mantle source of these basalts beneath the Mid-Atlantic Ridge over the last 35 Ma of its spreading activity. There is a good inverse correlation between 143Nd/144Nd and (La/Sm)EF with one exception in Hole 558 (approximately 35 Ma), the latter suggesting a recent (35 Ma) light REE depletion event, perhaps caused by dynamic or fractional melting. Short-term temporal and small-scale spatial mantle source variability is also evident in Hole 561 (approximately 18 Ma), which has rapid fluctuations in REE patterns and 143Nd/144Nd ratios (suggesting rapid transfer of magma from the time of melting) and is evidence contrary to the presence of a well-mixed magma chamber at this particular site and time. The mantle source variations noted can be interpreted within two extreme models. The first model invokes a convecting mantle depleted in large ion lithophile elements (LILE) and containing lumps (or veins) of LILE-enriched material of various shapes and sizes, passively and randomly distributed throughout. A second more restrictive model considers the interaction of fixed mantle plumes and the LILE-depleted asthenosphere flowing towards a migrating Mid- Atlantic Ridge (MAR) axis. With the exception of Hole 558 and the uncertainties of reconstructions of absolute plate movements in the region, the observed variations can be explained by two hot spots; the nearly ridge-centered Azores hot spot (plume) and another hot spot located beneath the African plate that may be affecting the source of basalts currently erupting at the MAR axis at 35°N and which, in the past, would have produced the New England chain of seamounts on the North American plate and (later) the Atlantis-Great Meteor chain on the African plate. Basalts erupted south of the Hayes Fracture Zone have not been affected by either of these two hot spots over the last 35 Ma and appear to have been continuously derived from the LILE-depleted source. Subaxial flow downridge from the Azores plume appears to have started 9 Ma, on the basis of the southward converging V-shaped time-transgressive ridges branching from the Pico and Corves Island, or not earlier than 16 Ma, on the basis of the geochemical results. Variations within Hole 558 remains unexplained by the latter model, unless we hypothesize a third hot spot.

A Late Pleistocene sea level stack

Late Pleistocene sea level has been reconstructed from ocean sediment core data using a wide variety of proxies and models. However, the accuracy of individual reconstructions is limited by measurement error, local variations in salinity and temperature, and assumptions particular to each technique. Here we present a sea level stack (average) which increases the signal-to-noise ratio of individual reconstructions. Specifically, we perform principal component analysis (PCA) on seven records from 0-430 ka and five records from 0-798 ka. The first principal component, which we use as the stack, describes ~80 % of the variance in the data and is similar using either five or seven records. After scaling the stack based on Holocene and Last Glacial Maximum (LGM) sea level estimates, the stack agrees to within 5 m with isostatically adjusted coral sea level estimates for Marine Isotope Stages 5e and 11 (125 and 400 ka, respectively). When we compare the sea level stack with the d18O of benthic foraminifera, we find that sea level change accounts for about ~40 % of the total orbital-band variance in benthic d18O, compared to a 65 % contribution during the LGM-to-Holocene transition. Additionally, the second and third principal components of our analyses reflect differences between proxy records associated with spatial variations in the d18O of seawater.

Coral calcifying fluid pH dictates response to ocean acidification

Ocean acidification driven by rising levels of CO2 impairs calcification, threatening coral reef growth. Predicting how corals respond to CO2 requires a better understanding of how calcification is controlled. Here we show how spatial variations in the pH of the internal calcifying fluid (pHcf) in coral (Stylophora pistillata) colonies correlates with differential sensitivity of calcification to acidification. Coral apexes had the highest pHcf and experienced the smallest changes in pHcf in response to acidification. Lateral growth was associated with lower pHcf and greater changes with acidification. Calcification showed a pattern similar to pHcf, with lateral growth being more strongly affected by acidification than apical. Regulation of pHcf is therefore spatially variable within a coral and critical to determining the sensitivity of calcification to ocean acidification.

In-situ Iron isotope measurements and mineral compositions from ODP Hole 206-1256D and IODP Hole 335-U1256D

In-situ Fe isotope measurements have been carried out to estimate the impact of the hydrothermal metamorphic overprint on the Fe isotopic composition of Fe-Ti-oxides and Fe-sulfides of the different lithologies of the drilled rocks from IODP Hole 1256D (eastern equatorial Pacific; 15 Ma crust formed at the East Pacific Rise). Most igneous rocks normally have a very restricted range in their 56Fe/54Fe ratio. In contrast, Fe isotope compositions of hot fluids (> 300 °C) from mid-ocean-ridge spreading centers define a narrow range that is shifted to lower delta 56Fe values by 0.2 per mil - 0.5 per mil as compared to igneous rocks. Therefore, it is expected that mineral phases that contain large amounts of Fe are especially affected by the interaction with a fluid that fractionates Fe isotopes during exsolution/precipitation of those minerals. We have used a femtosecond UV-Laser ablation system to determine mineral 56Fe/54Fe ratios of selected samples with a precision of < 0.1 per mil (2 sigma level) at micrometer-scale. We have found significant variations of the delta 56Fe (IRMM-014) values in the minerals between different samples as well as within samples and mineral grains. The overall observed scale of delta 56Fe (magnetite) in 1256D rocks ranges from - 0.12 to + 0.64 per mil, and of delta 56Fe (ilmenite) from - 0.77 to + 0.01 per mil. Pyrite in the lowermost sheeted dike section is clearly distinguishable from the other investigated lithological units, having positive delta 56Fe values between + 0.29 and + 0.56 per mil, whereas pyrite in the other samples has generally negative delta 56Fe values from - 1.10 to - 0.59 permil. One key observation is that the temperature dependent inter-mineral fractionations of Fe isotopes between magnetite and ilmenite are systematically shifted towards higher values when compared to theoretically expected values, while synthesized, well equilibrated magnetite-ilmenite pairs are compatible with the theoretical predictions. Theoretical considerations including beta-factors of different aqueous Fe-chlorides and Rayleigh-type fractionations in the presence of a hydrous, chlorine-bearing fluid can explain this observation. The disagreement between observed and theoretical equilibrium fractionation, the fact that magnetite, in contrast to ilmenite shows a slight downhole trend in the delta 56Fe values, and the observation of small scale heterogeneities within single mineral grains imply that a general re-equilibration of the magnetite-ilmenite pairs is overprinted by kinetic fractionation effects, caused by the interaction of magnetite/ilmenite with hydrothermal fluids penetrating the upper oceanic crust during cooling, or incomplete re-equilibration at low temperatures. Furthermore, the observation of significant small-scale variations in the 56Fe/54Fe ratios of single minerals in this study highlights the importance of high spatial-resolution-analyses of stable isotope ratios for further investigations.