ALBERT

Description: A Benguela Niño developed in November 2010 and lasted for 5 months along the Angolan and Namibian coastlines. Maximum amplitude was reached in January 2011 with an interannual monthly Sea Surface Temperature anomaly larger than 4 °C at the Angola Benguela Front. It was the warmest event since 1995. Consistent with previous Benguela Niños, this event was generated by a relaxation of the trade winds in the western equatorial Atlantic, which triggered a strong equatorial Kelvin wave propagating eastward along the equator and then southward along the southwest African coast. In the equatorial band, the associated ocean sub-surface temperature anomaly clearly shows up in data from the PIRATA mooring array. The dynamical signature is also detected by altimetry derived Sea Surface Height and is well reproduced by an Ocean Linear Model. In contrast to previous Benguela Niños, the initial propagation of sub-surface temperature anomalies along the equator started in October and the associated warming in the Angolan Benguela Front Zone followed on as early as November 2010. The warming was then advected further south in the Northern Benguela upwelling system as far as 25°S by an anomalously strong poleward sub-surface current. Demise of the event was triggered by stronger than normal easterly winds along the Equator in April and May 2011 leading to above normal shoaling of the thermocline along the Equator and the south-west African coastline off Angola and an associated abnormal equatorward current at the Angola Benguela Front in April and May 2011.

Description: Highlights • Analytical and numerical methods are employed to investigate fluid flow in active mud volcanoes or SHHS. • The effects of conduit radius and fluid properties on the flow rate are presented. • Conduit radius of such piercement systems cannot exceed a few metres at depth. • Clasts, if not densely packed, will not affect the flow rate when they are smaller than a fifth of the conduit size. • A maximal exsolution depth between 1800 and 3200 m is inferred for CH4 and between 750 and 1000 m for CO2. Clastic eruptions involve the rapid ascension of sedimentary clasts together with fluids, gas and/or liquid phases that may further deform and brecciate the host rocks. These fluids transport the resulting mixture, called mud breccia, to the surface. Such eruptions are often associated with geological structures such as mud volcanoes, hydrothermal vent complexes and, more generally, piercement structures. They involve various processes, acting over a wide range of scales, which makes them a complex and challenging multi-phase system to model. Although piercement structures have been widely studied and discussed, only a few attempts have been made to model the dynamics of such clastic eruptions. The ongoing Lusi mud eruption, in the East Java back-arc basin, which began in May 2006, is a spectacular large scale clastic eruption. The Lusi eruptive behaviour has been extensively studied over the past decade and thus represents a unique opportunity to better understand ongoing clastic eruptions and thus fossil clastic systems. We use both analytical formulations and numerical models to investigate simple relationships between the mud breccia properties (density, viscosity, gas and clast content) and the volumetric flow rate. Our results show that the conduit radius of such piercement systems cannot exceed a few metres at depth, and that clasts, if not densely packed, will not affect the flow rate when they are smaller than a fifth of the conduit size. Using published data for the annual gas fluxes at Lusi, we infer a maximal depth at which exsolution starts. This occurs between 1800 m and 3200 m depth for methane and between 750 m and 1000 m for carbon dioxide. Based on annual gas fluxes, we estimate that the conduit radius should be no larger than 1.5 m to match the maximal mud discharge, recorded at Lusi.

Description: The stoichiometric dissociation constants of carbonic acid ( and ) were determined by measurement of all four measurable parameters of the carbonate system (total alkalinity, total dissolved inorganic carbon, pH on the total proton scale, and CO2 fugacity) in natural seawater and seawater-derived brines, with a major ion composition equivalent to that Reference Seawater, to practical salinity (SP) 100 and from 25 °C to the freezing point of these solutions and –6 °C temperature minimum. These values, reported in the total proton scale, provide the first such determinations at below-zero temperatures and for SP 〉 50. The temperature (T, in Kelvin) and SP dependence of the current and (as negative common logarithms) within the salinity and temperature ranges of this study (33 ≤ SP ≤ 100, –6 °C ≤ t ≤ 25 °C) is described by the following best-fit equations: = –176.48 + 6.14528 – 0.127714 SP + 7.396×10–5 + (9914.37 – 622.886 + 29.714 SP) T–1 + (26.05129 – 0.666812 ) lnT (σ = 0.011, n = 62), and = –323.52692 + 27.557655 + 0.154922 SP – 2.48396×10–4 + (14763.287 – 1014.819 – 14.35223 SP) T–1 + (50.385807 – 4.4630415 ) lnT (σ = 0.020, n = 62). These functions are suitable for application to investigations of the carbonate system of internal sea ice brines with a conservative major ion composition relative to that of Reference Seawater and within the temperature and salinity ranges of this study.

Description: Documenting the early tectonic and magmatic evolution of the Izu–Bonin–Mariana (IBM) arc system in the Western Pacific is critical for understanding the process and cause of subduction initiation along the current convergent margin between the Pacific and Philippine Sea plates. Forearc igneous sections provide firm evidence for seafloor spreading at the time of subduction initiation (52 Ma) and production of “forearc basalt”. Ocean floor drilling (International Ocean Discovery Program Expedition 351) recovered basement-forming, low-Ti tholeiitic basalt crust formed shortly after subduction initiation but distal from the convergent margin (nominally reararc) of the future IBM arc (Amami Sankaku Basin: ASB). Radiometric dating of this basement gives an age range (49.3–46.8 Ma with a weighted average of 48.7 Ma) that overlaps that of basalt in the present-day IBM forearc, but up to 3.3 m.y. younger than the onset of forearc basalt activity. Similarity in age range and geochemical character between the reararc and forearc basalts implies that the ocean crust newly formed by seafloor spreading during subduction initiation extends from fore- to reararc of the present-day IBM arc. Given the age difference between the oldest forearc basalt and the ASB crust, asymmetric spreading caused by ridge migration might have taken place. This scenario for the formation of the ASB implies that the Mesozoic remnant arc terrane of the Daito Ridges comprised the overriding plate at subduction initiation. The juxtaposition of a relatively buoyant remnant arc terrane adjacent to an oceanic plate was more favourable for subduction initiation than would have been the case if both downgoing and overriding plates had been oceanic.

Description: The potential of crustose coralline algae as high-resolution archives of past ocean variability in mid- to high-latitudes has only recently been recognized. Few comparisons of coralline algal proxies, such as temperature-dependent algal magnesium to calcium (Mg/Ca) ratios, with in situ-measured surface ocean data exist, even rarer are well replicated records from individual sites. We present Mg/Ca records from nine coralline algal specimens (Clathromorphum compactum) from a single site in the Gulf of Maine, North Atlantic. Sections from algal mounds were analyzed using Laser Ablation-Inductively Coupled Plasma Mass Spectrometry (LA-ICP-MS) yielding individual Mg/Ca records of up to 30 years in length. We first test intra- and intersample signal replication and show that algal Mg/Ca ratios are reproducible along several transects within individual sample specimens and between different samples from the same study site. In addition, LA-ICP-MS-derived Mg/Ca ratios are compared to electron microprobe (EMP) analyzed data on the longest-lived specimens and were found to be statistically commensurable. Second, we evaluate whether relationships between algal-based SST reconstructions and in situ temperature data can be improved by averaging Mg/Ca records from multiple algal specimens (intersample averages). We found that intersample averages yield stronger relationships to sea surface temperature (SST) data than Mg/Ca records derived from individual samples alone. Thus, Mg/Ca-based paleotemperature reconstructions from coralline algae can benefit from using multiple samples per site, and can expand temperature proxy precision from seasonal to monthly.

Description: Changes in Atlantic Meridional Overturning Circulation (AMOC) strength exert a major influence on global atmospheric circulation patterns. However, the pacing and mechanisms of low-latitude responses to high-latitude forcing are insufficiently constrained so far. To elucidate the interaction of atmospheric and oceanic forcing in tropical South America during periods of major AMOC reductions (Heinrich Stadial 1 and the Younger Dryas) we generated a high-resolution foraminiferal multi-proxy record from off the Orinoco River based on Ba/Ca and Mg/Ca ratios, as well as stable isotope measurements. The data clearly indicate a three-phased structure of HS1 based on the reconfiguration of ocean currents in the tropical Atlantic Ocean. The initial phase (HS1a) is characterized by a diminished North Brazil Current, a southward displacement of the ITCZ, and moist conditions dominating northeastern Brazil. During subsequent HS1b, the NBC was even more diminished or yet reversed and the ITCZ shifted to its southernmost position. Hence, dryer conditions prevailed in northern South America, while eastern Brazil experienced maximally wet conditions. During the final stage, HS1c, conditions are similar to HS1a. The YD represents a smaller amplitude version of HS1 with a southward-shifted ITCZ. Our findings imply that the low-latitude continental climate response to high-latitude forcing is mediated by reconfigurations of surface ocean currents in low latitudes. Our new records demonstrate the extreme sensitivity of the terrestrial realm in tropical South America to abrupt perturbations in oceanic circulation during periods of unstable climate conditions.

Description: Olivine-hosted inclusions of silicate and sulfide melts, Cr-spinel and pyroxene were studied to estimate magma composition, temperature, pressure, and fO2 at the onset and during the silicate-sulfide immiscibility in modern arc basalt from Tolbachik volcano, Kamchatka arc. We demonstrate that the olivine phenocrysts hosting sulfide and silicate melt inclusions belong to the same population. The compositions of the silicate melt inclusions in most primitive olivine (88–91 mol% Fo) represent moderately oxidized (~ QFM + 1.1) high-MgO (up to 12–12.6 wt%) and high CaO/Al2O3 (0.8–1.2) melt that has abundances and ratios of the lithophile trace elements typical of island arc magmas. The initial volatile contents in parental Tolbachik magma are estimated from the melt inclusions and mass-balance considerations to be at least 4.9 wt% H2O, 2600 ppm S, 1100 ppm Cl, 550 ppm F, and 1200 ppm CO2. These data are used to calculate the temperature (~ 1220 °C) and minimum pressure (3 kbar) at which the beginning of crystallization and exsolution of sulfide melt took place. The presence of anhydrite, especially ubiquitous in the crystallized silicate melt associated with sulfide globules, suggest that much higher sulfur abundances prior to degassing and sulfate immiscibility and/or crystallization should be expected. We tentatively considered hydrothermal accumulations of sulfur (elemental, sulfate and sulfide) in the volcanic conduit responsible for local contamination and oversaturation of the Tolbachik magma in sulfur and related sulfide immiscibility. Coexisting sulfide and sulfate can be also interpreted in favor of the magmatic sulfide oxidation and related generation of S-rich fluids. Such fluids are expected to accumulate metals released from decomposed sulfide melts and supply significant epithermal mineralization, including native gold.

Description: Coccoliths comprise a major fraction of the global carbonate sink. Therefore, changes in coccolithophores' Ca isotopic fractionation could affect seawater Ca isotopic composition, affecting interpretations of the global Ca cycle and related changes in seawater chemistry and climate. Despite this, a quantitative interpretation of coccolith Ca isotopic fractionation and a clear understanding of the mechanisms driving it are not yet available. Here, we address this gap in knowledge by developing a simple model (CaSri–Co) to track coccolith Ca isotopic fractionation during cellular Ca uptake and allocation to calcification. We then apply it to published and new δ44/40Ca and Sr/Ca data of cultured coccolithophores of the species Emiliania huxleyi and Gephyrocapsa oceanica. We identify changes in calcification rates, Ca retention efficiency and solvation–desolvation rates as major drivers of the Ca isotopic fractionation and Sr/Ca variations observed in cultures. Higher calcification rates, higher Ca retention efficiencies and lower solvation–desolvation rates increase both coccolith Ca isotopic fractionation and Sr/Ca. Coccolith Ca isotopic fractionation is most sensitive to changes in solvation–desolvation rates. Changes in Ca retention efficiency may be a major driver of coccolith Sr/Ca variations in cultures. We suggest that substantial changes in the water structure strength caused by past changes in temperature could have induced significant changes in coccolithophores' Ca isotopic fractionation, potentially having some influence on seawater Ca isotopic composition. We also suggest a potential effect on Ca isotopic fractionation via modification of the solvation environment through cellular exudates, a hypothesis that remains to be tested.

Description: Alphaproteobacterium strain MOLA1416, related to Mycoplana ramosa DSM 7292 and Chelativorans intermedius CC-MHSW-5 (93.6% 16S rRNA sequence identity) was isolated from the marine lichen, Lichina pygmaea and its chemical composition was characterized by a metabolomic network analysis using LC-MS/MS data. Twenty-five putative different compounds were revealed using a dereplication workflow based on MS/MS signatures available through GNPS (https://gnps.ucsd.edu/). In total, ten chemical families were highlighted including isocoumarins, macrolactones, erythrinan alkaloids, prodiginines, isoflavones, cyclohexane-diones, sterols, diketopiperazines, amino-acids and most likely glucocorticoids. Among those compounds, two known metabolites (13 and 26) were isolated and structurally identified and metabolite 26 showed a high cytotoxic activity against B16 melanoma cell lines with an IC50 0.6 ± 0.07 μg/mL.

Description: Highlights: • Lagrangian ocean analysis is a powerful way to analyse the output of ocean circulation models • We present a review of the Kinematic framework, available tools, and applications of Lagrangian ocean analysis • While there are unresolved questions, the framework is robust enough to be used widely in ocean modelling Abstract: Lagrangian analysis is a powerful way to analyse the output of ocean circulation models and other ocean velocity data such as from altimetry. In the Lagrangian approach, large sets of virtual particles are integrated within the three-dimensional, time-evolving velocity fields. Over several decades, a variety of tools and methods for this purpose have emerged. Here, we review the state of the art in the field of Lagrangian analysis of ocean velocity data, starting from a fundamental kinematic framework and with a focus on large-scale open ocean applications. Beyond the use of explicit velocity fields, we consider the influence of unresolved physics and dynamics on particle trajectories. We comprehensively list and discuss the tools currently available for tracking virtual particles. We then showcase some of the innovative applications of trajectory data, and conclude with some open questions and an outlook. The overall goal of this review paper is to reconcile some of the different techniques and methods in Lagrangian ocean analysis, while recognising the rich diversity of codes that have and continue to emerge, and the challenges of the coming age of petascale computing.

Description: Highlights • Comparison of encrusting and rhodolith coralline algae for paleoclimate reconstruction • Both coralline algal forms can yield SST information, but encrusting forms generally yield higher correlations to SST. • Encrusting morphologies yielded longer records due to frequent growth irregularities in rhodoliths. Abstract Coralline algae have been used for sclerochronological studies throughout the last decade. These studies have focused on two different growth morphologies of the photosynthetic coralline algae: massive crusts forming small buildups on hard substrate, and free-living branching algal nodules, known as rhodoliths. The latter are generally found on soft-substrate, where they are frequently overturned by water movement and bottom feeding organisms, leaving one side of the rhodolith partially buried in the sediment at any given time. Here we test whether either of these growth morphologies is more suitable for proxy reconstructions by comparing Mg/Ca ratios – a temperature proxy – in multiple replicates of rhodoliths of Lithothamnion glaciale and in rhodoliths as well as encrusting specimens of Clathromorphum compactum. With both species being widespread throughout the Temperate and Arctic regions, we have chosen two North Atlantic localities at Nuuk Fjord, Greenland (Subarctic), and off the southeastern coast of Newfoundland, Canada (Temperate), for this study. Two to three Mg/Ca ratio transects spanning 18 years of growth were analysed on multiple specimens with encrusting morphologies and along different sides of rhodoliths using laser ablation inductively coupled mass spectrometry and compared to remotely sensed sea surface temperature (SST) data. The length of the common time span used for comparison was limited by growth interruptions in rhodoliths. Furthermore, our comparison is based on the assumption that rhodolith growth increments are annual – an assumption that has recently been challenged by mesocosm studies. Monthly Mg/Ca values from multiple transects within each individual were compared and in samples from Nuuk fjord significant correlations were found in 4 of 4 encrusting C. compactum, 4 of 4 C. compactum rhodoliths, and 2 of 3 L. glaciale rhodoliths. In Newfoundland significant correlations were found in 6 of 6 encrusting C. compactum comparisons (average: r = 0.61, p 〈 0.001), and in 6 of 6 L. glaciale rhodolith comparisons (average: r = 0.43, p 〈 0.001) for monthly resolved time series. The monthly Mg/Ca ratios (n = 216) from each morphology were compared with instrumental Reynolds SST yielding the following correlations: encrusting C. compactum (r = 0.64, p 〈 0.001), C. compactum rhodolith (r = 0.62, p 〈 0.001) and L. glaciale (r = 0.58, p 〈 0.001). In Newfoundland both morphologies indicate a similar strength in recording SST: encrusting C. compactum (r = 0.85, p 〈 0.001) and rhodolith-forming L. glaciale (r = 0.84, p 〈 0.001). In summary, Mg/Ca ratios derived from both coralline algal growth forms can yield SST information, however, massive encrusting forms generally yield higher correlations to SST than transects measured on individual rhodoliths, which only allowed for the generation of short uninterrupted time series due to frequent growth irregularities.

Description: The subterranean estuary (STE), the subsurface mixing zone of outflowing fresh groundwater and infiltrating seawater, is an area of extensive geochemical reactions that determine the composition of groundwater that flows into coastal environments. This study examined the porewater composition of a shallow STE (〈5 m depth) in Gloucester Point, VA (USA) over two years to determine seasonal variations in dissolved organic carbon (DOC) and the reduced metabolites Fe, Mn, and sulfide. An additional aim of this study was to investigate the relative importance of salinity gradients (which have great geochemical influence in surface estuaries) versus redox gradients on STE geochemistry. Two freshwater endmembers were identified, between which redox potential and composition varied with depth-a shallow freshwater endmember was oxidizing and high in DOC, whereas a deep freshwater endmember was reducing, lower in DOC, and high in sulfide. Results showed that dissolved Fe, Mn, and sulfide varied along a redox gradient distinct from the salinity gradient, and that three-endmember mixing was required to quantify non-conservative chemical addition/removal in the STE. In addition to salinity, humic carbon was used as a quasi-conservative tracer to quantify mixing according to a three-endmember model. The vertical distributions of DOC and reduced metabolites remained approximately constant over time, but concentrations varied with season. Dissolved organic carbon concentrations were greatest in the summer, and shallow meteoric groundwater supplied the majority of DOC to the STE. In summer, there was additional evidence for shallow non-conservative addition of DOC. Dissolved Fe and Mn were highest in a subsurface plume through the middle of the STE (100-140 cm below sediment surface at the high tide line) which was characterized by higher concentrations and greater non-conservative addition in the winter. In contrast, sulfide was higher in summer at depths within the Fe and Mn plume (100-140 cm). We attribute the contrasting seasonal patterns of dissolved Fe, Mn, and sulfide to differences in microbial response to temperature changes and organic matter availability, and to competition at the ferrous-sulfidic transition zone between dissimilatory metal reduction and sulfate reduction, leading to sulfate/sulfur reducing bacteria (SRB) being more active in summer and metal reducers being more active in winter. Throughout the STE, seasonal temperature and DOC variations determined the spatial distribution and geochemical cycling of Fe, Mn, and sulfur.

Description: Highlights: • Comparison of global NEMO and FESOM configurations with emphasis on the Agulhas system. • Both models simulate a reasonable and comparable large-scale circulation. • Both models have individual strengths and weaknesses to match the observations of the WBC system. • The numerical cost of FESOM is twice the one of NEMO. Abstract: Many questions in ocean and climate modelling require the combined use of high resolution, global coverage and multi-decadal integration length. For this combination, even modern resources limit the use of traditional structured-mesh grids. Here we compare two approaches: A high-resolution grid nested into a global model at coarser resolution (NEMO with AGRIF) and an unstructured-mesh grid (FESOM) which allows to variably enhance resolution where desired. The Agulhas system around South Africa is used as a testcase, providing an energetic interplay of a strong western boundary current and mesoscale dynamics. Its open setting into the horizontal and global overturning circulations also requires global coverage. Both model configurations simulate a reasonable large-scale circulation. Distribution and temporal variability of the wind-driven circulation are quite comparable due to the same atmospheric forcing. However, the overturning circulation differs, owing each model's ability to represent formation and spreading of deep water masses. In terms of regional, high-resolution dynamics, all elements of the Agulhas system are well represented. Owing to the strong nonlinearity in the system, Agulhas Current transports of both configurations and in comparison with observations differ in strength and temporal variability. Similar decadal trends in Agulhas Current transport and Agulhas leakage are linked to the trends in wind forcing. Although the number of 3D wet grid points used in FESOM is similar to that in the nested NEMO, FESOM uses about two times the number of CPUs to obtain the same model throughput (in terms of simulated model years per day). This is feasible due to the high scalability of the FESOM code.

Description: Highlights • Widespread Holocene tephra KHG, from Kamchatka, is found as a cryptotephra in the NGRIP ice-core. • This is the first identification of tephra from the Kamchatka Peninsula in Greenland ice. • NGRIP KHG has an age of 7872 ± 50 a BP and improves age models for Kamchatka. • Existing 14C age estimates for the KHG eruption are too young. Abstract Contiguous sampling of Holocene ice from the NGRIP core, Greenland, has revealed a new rhyolitic cryptotephra that is geochemically identical to the KHG tephra, a widespread marker deposit originating from the Khangar volcano, Kamchatka. This is the first identification of tephra from the Kamchatka Peninsula in Greenland ice and the first finding of the KHG tephra outside Kamchatka. The NGRIP KHG has an age of 7872 ± 50 a BP 1950, and this date will help improve age models for Kamchatka, where existing age estimates of KHG are too young, thus highlighting the importance of locating long-range, low-concentration cryptotephra deposits in well-dated ice cores. In Greenland KHG is located close to the termination of the 8.2 ka BP cooling event that is also a climate feature in palaeo-records of Kamchatka. This tie-point therefore provides a unique opportunity to synchronise records of environmental change in distal locations.

Description: Highlights • Active hydrothermal vent ecosystems are extremely rare. • Vent ecosystems are recognized as vulnerable by international organizations. • Mineral resources at active vents would not contribute significantly to the global metal supply. • Effective networks that protect representative active vents cannot be ensured. • A prohibition on mining active vents is consistent with obligations for conservation. Abstract There is increasing interest in mining minerals on the seabed, including seafloor massive sulfide deposits that form at hydrothermal vents. The International Seabed Authority is currently drafting a Mining Code, including environmental regulations, for polymetallic sulfides and other mineral exploitation on the seabed in the area beyond national jurisdictions. This paper summarizes 1) the ecological vulnerability of active vent ecosystems and aspects of this vulnerability that remain subject to conjecture, 2) evidence for limited mineral resource opportunity at active vents, 3) non-extractive values of active vent ecosystems, 4) precedents and international obligations for protection of hydrothermal vents, and 5) obligations of the International Seabed Authority under the UN Convention on the Law of the Sea for protection of the marine environment from the impacts of mining. Heterogeneity of active vent ecosystems makes it extremely challenging to identify “representative” systems for any regional, area-based management approach to conservation. Protection of active vent ecosystems from mining impacts (direct and indirect) would set aside only a small fraction of the international seabed and its mineral resources, would contribute to international obligations for marine conservation, would have non-extractive benefits, and would be a precautionary approach.

Description: From protists to humans, all animals and plants are inhabited by microbial organisms. There is an increasing appreciation that these resident microbes influence the fitness of their plant and animal hosts, ultimately forming a metaorganism consisting of a uni- or multicellular host and a community of associated microorganisms. Research on host–microbe interactions has become an emerging cross-disciplinary field. In both vertebrates and invertebrates a complex microbiome confers immunological, metabolic and behavioural benefits; conversely, its disturbance can contribute to the development of disease states. However, the molecular and cellular mechanisms controlling the interactions within a metaorganism are poorly understood and many key interactions between the associated organisms remain unknown. In this perspective article, we outline some of the issues in interspecies interactions and in particular address the question of how metaorganisms react and adapt to inputs from extreme environments such as deserts, the intertidal zone, oligothrophic seas, and hydrothermal vents

Description: Highlights • In cold seeps of Guaymas Basin, aragonite, barite and pyrite precipitated from modified seawater. • Aragonite is highly depleted in 13C suggesting formation via anaerobic oxidation of methane. • Barite formed through mixing of reducing, Ba-rich seep fluids with a 34S-rich sulfate pool. • Pyrite framboids formed under anoxic-sulfidic water via microbial sulfate reduction. Abstract Authigenic carbonate crusts, surface muds and bivalve shell fragments have been recovered from inactive and active recently discovered cold seep sites in central Guaymas Basin. In this study, for first time, redox conditions and fluid sources involved in mineral precipitation were investigated by analyzing the mineralogy and textures of surface samples, along with skeletal contents, and C, O and S isotopes variations. The δ13C values of aragonitic bivalve shells and non-skeletal carbonate from some surface muds (1‰ to −3.7‰ V-PDB) suggest that carbonate precipitated from ambient dissolved inorganic carbon, whereas fibrous aragonite cement and non-skeletal carbonate from other sites are highly depleted in 13C (down to −47.6‰ V-PDB), suggesting formation via anaerobic oxidation of methane, characteristic of methane seepage environments. δ18O in most of the carbonates varies from +1.4‰ to +3.2‰ V-PDB, indicating that they formed from slightly modified seawater. Some non-skeletal carbonate grains from surface muds have lower δ18O values (−12.5‰ to −8.2‰ V-PDB) reflecting the influence of 18O-depleted pore water. Size distribution of pyrite framboids (mean value: 3.1 μm) scattered within diatomaceous sinter suggests formation from anoxic-sulfidic bottom waters. δ34S in pyrite is of −0.3‰ V-CDT compared to +46.6‰ V-CDT in barite, thus implying a fluid sulfate−sulfide fractionation of 21.3‰ that argues in favor of microbial sulfate reduction as the processes that mediated pyrite framboid formation, in a semi-closed system. Barite formation occurred through the mixing of reducing and Ba-rich seep fluids with a 34S-enriched sulfate pool that resulted from microbial sulfate reduction in a semi-closed system. The chemical composition of aragonite cement, barite and pyrite suggest mineral precipitation from modified seawater. Taken together, our data suggest that mineralization at the studied seep sites is controlled by the mixing of seawater with minor amounts of hydrothermal fluids, and oxygen-depleted conditions favoring anaerobic microbial processes.

Description: The submerged portions of the North-Anatolian Fault (NAF) in the Sea of Marmara and the NE-Aegean Sea are sites of large magnitude earthquakes, that leave diagnostic geological “signatures” in the sedimentary record in the form of mass-wasting deposits, turbidites, and fluid and gas escape features. This is due to the interplay of seismic-shaking, mass- and turbidity flows, sediment resuspension and fluids circulation in relatively small sub-basins with a complex paleo-oceanography, steep slopes, high rates of deformation, and diffuse fault-controlled gas and fluid seeps. To unravel the complex interrelations of these phenomena during earthquake cycles, we carried out paleoseismological studies at several key locations. Here, we report results of these studies, carried out onboard the R/V Urania over a decade, starting soon after the Mw 7.4, 1999 İzmit earthquake. Our work included high resolution mapping of active faults through multibeam bathymetry and high resolution seismic reflection profiles, multi-parameter analysis of sediment cores, as well as seafloor observations using sensors mounted on remotely-operated vehicles (ROV). The main objectives were to map active faults, determine slip-rates and earthquake recurrence times along major fault strands, and assess connections between fault deformation and fluid activity. We mapped fault geometry in the gulfs of İzmit, Gemlik and Saros, showing the trans-tensive nature of these depressions. The average slip-rates for the last ~ 10 ka was found to be 10 mm/y in the gulfs of İzmit and Saros, at the eastern and the western ends of the NAF northern strand, and 3–4 mm/yr in the Gulf of Gemlik, along the middle strand of the NAF. These rates, integrated over 10 ka of NAF activity, are smaller than those determined by the GPS geodetic measurements. Submarine paleoseismological studies in the Gulf of İzmit detected the sedimentary records of earthquakes for the last 2.4 ka, suggesting an average recurrence time of 300 years for major events. Multisensor observations and monitoring of the seafloor have shown widespread emissions of gas and fluids along the submerged part of the NAF, associated with reduced black sediments; we investigated their possible connection with the earthquake cycle.

Description: Highlights • Four of the seven seamounts northeast of the Galápagos Platform are drowned islands • The ages of the seamounts range from 5.2 Ma to 0.5 Ma • Seamount morphology changes from conical to elongate at ~1.5 Ma • The locus of volcanism appears to migrate eastward at the rate of Nazca plate motion Abstract We present new geochemical and 40Ar/39Ar analyses from seven seamounts located off the northeastern margin of the shallow Galápagos Platform. Initial volcanism at 5.2 Ma created a small island (Pico) over the current location of the hotspot with geochemically enriched lavas. There is no further record of magmatism in the study area until 3.8 to 2.5 Ma, during which four roughly conical volcanoes (Sunray, Grande, Fitzroy, and Beagle) formed through eruption of lavas derived from a depleted mantle source. Sunray, Fitzroy, and Grande were islands that existed for ~3 m.y. ending with the submergence of Fitzroy at ~0.5 Ma. The youngest seamounts, Largo and Iguana, do not appear to have been subaerial and were active at 1.3 Ma and 0.5 Ma, respectively, with the style of edifice changing from the previous large cones to E-W elongate, composite structures. The progression of magmatism suggests that Pico erupted near 91.5°W near the location of the Galápagos plume while the others formed well east of the plume center. If the locations of initial volcanism are calculated using the eastward velocity of the Nazca plate, there appears to be a progression of younger volcanism toward the east, opposite what would be expected from a fixed mantle plume source. The rate that initial volcanism moves eastward is close to the plate velocity. A combination of higher temperature and geochemical enrichment of the thickened lithosphere of the Galápagos platform could have provided a viscosity gradient at the boundary between the thick lithosphere and the thinner oceanic lithosphere to the northeast. As this boundary moved eastward with the Nazca plate, it progressively triggered shear-driven mantle upwelling and volcanism.

Description: Highlights • Melt inclusions from southern Payenia have highly variable element enrichment • Magmas formed by mixing of asthenospheric high Nb/U and lithospheric low Nb/U melts • Low Nb/U type inclusions are similar in composition to alkaline lamprophyres • Low Nb/U melts were formed by fractionation of high Nb/U melts in the SCLM • The percolative fractional crystallization involved cpx, rutile and apatite Abstract We present major and trace element compositions of melt inclusions from three alkali basalts from the Río Colorado volcanic field in the Payenia backarc province, Argentina. Modeling of diffusion profiles around the inclusions showed that most inclusions equilibrated 〈14 days after formation, indicating a short crustal residence time for the magmas and nearly direct ascent through the crust. Despite overlapping host rock isotopic compositions, the inclusions show a large variation in their degree of enrichment, and display trends that we interpret as mixing between asthenospheric OIB-type low K2O-high Nb/U melts and enriched high K2O-low Nb/U lithospheric mantle melts similar in composition to alkaline lamprophyres. The low Nb/U magmas are excessively enriched in the elements Cs, Rb, Ba, Th, U, K, Pb and Cl relative to Nb, Ta and REEs. The enriched low Nb/U components are interpreted to have formed by percolative fractional crystallization of asthenospheric high Nb/U melts in the lithospheric mantle involving crystallization of clinopyroxene, apatite and rutile. The residual fluid-rich melts either mixed directly with new batches of high Nb/U melts or metasomatized and veined the lithospheric mantle which later re-melted during continued volcanism. The major element compositions of the high K2O-low Nb/U components are distinct for the whole rocks and melt inclusions, and most enriched inclusions have lower SiO2 and higher TiO2 contents indicating derivation by melting of amphibole-bearing veins. In contrast, most wr low Nb/U basalts have higher SiO2 and lower TiO2 and were most likely formed by melting of pyroxenitic veins or peridotitic metasomatized lithospheric mantle.

Description: The spectacular Lusi eruption started in northeast Java, Indonesia, on May 29th, 2006, continuously erupting mud, water, gas, oil, and clasts ever since. Lusi provides an unprecedented opportunity to study the birth and the evolution of a large-scale and hot mud eruption. Lusi is interpreted as a hybrid between a traditional hydrocarbon-driven piercement structure (mud volcano) and a hydrothermal system fuelled by magmatic heat. Lusi is therefore an exciting natural laboratory for understanding analogue modern and palaeo-piercement systems such as mud volcanoes, sediment-hosted hydrothermal systems, and hydrothermal vent complexes. This special issue collects recent multidisciplinary work completed in the framework of the ERC-funded LUSI LAB project. These studies were conducted at and near Lusi. Contributions span across disciplines such as engineering, geochemistry, geophysics, geology and numerical modelling, including fieldwork, laboratory and theoretical approaches. The acquired results contribute to characterise the dynamics of complex interactions between volcanism and an ongoing erupting clastic system. Lusi still saves many mysteries that will be unravelled by future scientific investigations.

Description: Microplastics (〈5 mm) have been found in many fish species, from most marine environments. However, the mechanisms underlying microplastic ingestion by fish are still unclear, although they are important to determine the pathway of microplastics along marine food webs. Here we conducted experiments in the laboratory to examine microplastic ingestion (capture and swallowing) and egestion by juveniles of the planktivorous palm ruff, Seriolella violacea (Centrolophidae). As expected, fish captured preferentially black microplastics, similar to food pellets, whereas microplastics of other colours (blue, translucent, and yellow) were mostly co-captured when floating close to food pellets. Microplastics captured without food were almost always spit out, and were only swallowed when they were mixed with food in the fish's mouth. Food probably produced a ‘gustatory trap’ that impeded the fish to discriminate and reject the microplastics. Most fish (93% of total) egested all the microplastics after 7 days, on average, and 49 days at most, substantially longer than food pellets (〈2 days). No acute detrimental effects of microplastics on fish were observable, but potential sublethal effects of microplastics on the fish physiological and behavioural responses still need to be tested. This study highlights that visually-oriented planktivorous fish, many species of which are of commercial value and ecological importance within marine food webs, are susceptible to ingest microplastics resembling or floating close to their planktonic prey. Microplastic similarity to food, and food presence together with microplastics, influence the ingestion of microplastics by visually-oriented planktivorous fish.

Description: Highlights • A carbonic anhydrase gene was identified from an hydrothermal vent metagenome. • The gene product LOGACA is a dimeric α-type carbonic anhydrase. • LOGACA is highly thermostable at alkaline pH. • Thermostability correlates with secondary structure, surface charges and ion pairs. Abstract Carbonic anhydrases (CAs) are extremely fast enzymes, which have attracted much interest in the past due to their medical relevance and their biotechnological potential. An α-type CA gene was isolated from DNA derived from an active hydrothermal vent chimney, in an effort to identify novel CAs with suitable properties for CO2 capture. The gene product was recombinantly produced and characterized, revealing remarkable thermostability, also in the presence of high ionic strength alkaline conditions, which are used in some CO2 capture applications. The Tm was above 90 °C under all tested conditions. The enzyme was crystallized and the structure determined by molecular replacement, revealing a typical bacterial α-type CA non-covalent dimer, but not the disulphide mediated tetramer observed for the hyperthermophilic homologue used for molecular replacement, from Thermovibrio ammonificans. Structural comparison suggests that an increased secondary structure content, increased content of charges on the surface and ionic interactions compared to mesophilic enzymes, may be main structural sources of thermostability, as previously suggested for the homologue from Sulfurihydrogenibium yellowstonense.

Description: Highlights • Improved position accuracy of up to 0.08 cm using SAPOS®-corrected DGNSS • High-resolution snippet backscatter to detect and monitor UXO and its environment • Combined evaluation of bathymetry, backscatter intensities, and sediments Abstract The present study reports the evaluation of snippet backscatter information gathered with a high-frequency multibeam echosounder system (200–400 kHz) due to their usability to detect ammunition of different sizes in shallow coastal waters. Besides the feasibility of the snippet backscatter data, it was focused on the attainable horizontal accuracy in comparison to side-scan sonar and autonomous underwater vehicle (AUV) surveys. The data was collected in shallow coastal waters of up to 18 m water depth (Baltic Sea) close to an ammunition dumping site characterized by an almost flat seafloor covered with sand and silt sediments. The analysis of the multibeam compared to sidescan data indicates the snippet backscatter to be a promising prospective method for ammunition detection and being able to improve horizontal position accuracy of up to 0.08 m.

Description: Seafloor spreading at slow and ultraslow rates is often taken up by extension on large-offset faults called detachments, which exhume lower crustal and mantle rocks, and in some cases make up domed oceanic core complexes. The exposed footwall may reveal a characteristic pattern of spreading-parallel corrugations, the largest of which are clearly visible in multibeam bathymetric data, and whose nature and origin have been the subject of controversy. In order to tackle this debate, we use available nearbottom bathymetric surveys recently acquired with autonomous deep-sea vehicles over five corrugated detachments along the Mid-Atlantic Ridge. With a spatial resolution of 2 m, these data allow us to compare the geometry of corrugations on oceanic detachments that are characterized by differing fault zone lithologies, and accommodate varying amounts of slip. The fault surfaces host corrugations with wavelengths of 10-250 m, while individual corrugations are finite in length, typically 100-500 m. Power spectra of profiles calculated across the corrugated fault surfaces reveal a common level of roughness, and indicate that the fault surfaces are not fractal. Since systematic variation in roughness with fault offset is not evident, we propose that portions of the exposed footwalls analyzed here record constant brittle strain. We assess three competing hypotheses for corrugation formation and find that the continuous casting and varying depth to brittle-ductile transition models cannot explain the observed corrugation geometry nor available geological observations. We suggest a model involving brittle strain localization on a network of linked fractures within a zone of finite thickness is a better explanation for the observations. This model explains corrugations on oceanic detachment faults exposed at the seafloor and on normal faults in the continents, and is consistent with recently imaged corrugations on a subduction zone megathrust. Hence fracture linkage and coalescence may give rise to corrugated fault zones, regardless of earlier deformation history and tectonic setting.

Description: The relationship between fisheries and marine spatial planning (MSP) is still widely unsettled. While several scientific studies highlight the strong relation between fisheries and MSP, as well as ways in which fisheries could be included in MSP, the actual integration of fisheries into MSP often fails. In this article, we review the state of the art and latest progress in research on various challenges in the integration of fisheries into MSP. The reviewed studies address a wide range of integration challenges, starting with techniques to analyse where fishermen actually fish, assessing the drivers for fishermen's behaviour, seasonal dynamics and long-term spatial changes of commercial fish species under various anthropogenic pressures along their successive life stages, the effects of spatial competition on fisheries and projections on those spaces that might become important fishing areas in the future, and finally, examining how fisheries could benefit from MSP. This paper gives an overview of the latest developments on concepts, tools, and methods. It becomes apparent that the spatial and temporal dynamics of fish and fisheries, as well as the definition of spatial preferences, remain major challenges, but that an integration of fisheries is already possible today

Description: Highlights • The Danube deep-sea fan offers best conditions for hydrate production. • Gas production out of a hypothetical methane hydrate reservoir was simulated. • Hazard assessment to investigate the hazard of production-induced slope failures. • Factor of Safety against slope failure is not affected by the production process. • Mobilized mass could hit the production site if landslide were to happen. Methane production from gas hydrate reservoirs is only economically viable for hydrate reservoirs in permeable sediments. The most suitable known prospect in European waters is the paleo Danube deep-sea fan in the Bulgarian exclusive economic zone in the Black Sea where a gas hydrate reservoir is found 60 m below the seafloor in water depths of about 1500 m. To investigate the hazards associated with gas production-induced slope failures we carried out a slope stability analysis for this area. Screening of the area based on multibeam bathymetry data shows that the area is overall stable with some critical slopes at the inner levees of the paleo channels. Hydrate production using the depressurization method will increase the effective stresses in the reservoir beyond pre-consolidation stress, which results in sediment compaction and seafloor subsidence. The modeling results show that subsidence would locally be in the order of up to 0.4 m, but it remains confined to the immediate vicinity above the production site. Our simulations show that the Factor of Safety against slope failure (1.27) is not affected by the production process, and it is more likely that a landslide is triggered by an earthquake than by production itself. If a landslide were to happen, the mobilized sediments on the most likely failure plane could generate a landslide that would hit the production site with velocities of up to 10 m s-1. This case study shows that even in the case of production from very shallow gas hydrate reservoirs the threat of naturally occurring slope failures may be greater than that of hydrate production itself and has to be considered carefully in hazard assessments.

Description: Deep-sea ecosystems, limited by their inability to use primary production as a source of carbon, rely on other sources to maintain life. Sedimentation of organic carbon into the deep sea has been previously studied, however, the high biomass of sedimented Sargassum algae discovered during the VEMA Transit expedition in 2014/2015 to the southern North Atlantic, and its potential as a regular carbon input, has been an underestimated phenomenon. To determine the potential for this carbon flux, a literature survey of previous studies that estimated the abundance of surface water Sargassum was conducted. We compared these estimates with quantitative analyses of sedimented Sargassum appearing on photos taken with an autonomous underwater vehicle (AUV) directly above the abyssal sediment during the expedition. Organismal communities associated to Sargassum fluitans from surface waters were investigated and Sargassum samples collected from surface waters and the deep sea were biochemically analyzed (fatty acids, stable isotopes, C:N ratios) to determine degradation potential and the trophic significance within deep-sea communities. The estimated Sargassum biomass (fresh weight) in the deep sea (0.07 − 3.75 g/m2) was several times higher than that estimated from surface waters in the North Atlantic (0.024 – 0.84 g/m2). Biochemical analysis showed degradation of Sargassum occurring during sedimentation or in the deep sea, however, fatty acid and stable isotope analysis did not indicate direct trophic interactions between the algae and benthic organisms. Thus, it is assumed that components of the deep-sea microbial food web form an important link between the macroalgae and larger benthic organisms. Evaluation of the epifauna showed a diverse nano- micro-, meio, and macrofauna on surface Sargassum and maybe transported across the Atlantic, but we had no evidence for a vertical exchange of fauna components. The large-scale sedimentation of Sargassum forms an important trophic link between surface and benthic production and has to be further considered in the future as a regular carbon input to the deep-sea floor in the North Atlantic.

Description: Highlights: • First long-term study on microplastic in the marine environment • Case study based on a unique sample set from the highly human impacted Baltic Sea • Water column microplastic concentration constant over past three decades • Microplastic concentration in forage fish constant over past three decades • We hypothesise that household waste is the dominant source of Baltic marine plastics. Abstract Microplastic is considered a potential threat to marine life as it is ingested by a wide variety of species. Most studies on microplastic ingestion are short-term investigations and little is currently known about how this potential threat has developed over the last decades where global plastic production has increased exponentially. Here we present the first long-term study on microplastic in the marine environment, covering three decades from 1987 to 2015, based on a unique sample set originally collected and conserved for food web studies. We investigated the microplastic concentration in plankton samples and in digestive tracts of two economically and ecologically important planktivorous forage fish species, Atlantic herring (Clupea harengus) and European sprat (Sprattus sprattus), in the Baltic Sea, an ecosystem which is under high anthropogenic pressure and has undergone considerable changes over the past decades. Surprisingly, neither the concentration of microplastic in the plankton samples nor in the digestive tracts changed significantly over the investigated time period. Average microplastic concentration in the plankton samples was 0.21±0.15particlesm-3. Of 814 fish examined, 20% contained plastic particles, of which 95% were characterized as microplastic (〈5mm) and of these 93% were fibres. There were no significant differences in the plastic content between species, locations, or time of day the fish were caught. However, fish size and microplastic in the digestive tracts were positively correlated, and the fish contained more plastic during summer than during spring, which may be explained by increased food uptake with size and seasonal differences in feeding activity. This study highlights that even though microplastic has been present in the Baltic environment and the digestive tracts of fishes for decades, the levels have not changed in this period. This underscores the need for greater understanding of how plastic is cycled through marine ecosystems. The stability of plastic concentration and contamination over time observed here indicates that the type and level of microplastic pollution may be more closely correlated to specific human activities in a region than to global plastic production and utilization as such.

Description: Highlights: • We present the first hyperspectral image data from the deep seafloor. • The data were acquired with a new UHI in 4200 m water depth. • Supervised classification is able to detect manganese nodules and fauna. • The UHI is a promising tool for high-resolution seafloor exploration and monitoring. Abstract: Hyperspectral seafloor surveys using airborne or spaceborne sensors are generally limited to shallow coastal areas, due to the requirement for target illumination by sunlight. Deeper marine environments devoid of sunlight cannot be imaged by conventional hyperspectral imagers. Instead, a close-range, sunlight-independent hyperspectral survey approach is required. In this study, we present the first hyperspectral image data from the deep seafloor. The data were acquired in approximately 4200 m water depth using a new Underwater Hyperspectral Imager (UHI) mounted on a remotely operated vehicle (ROV). UHI data were recorded for 112 spectral bands between 378 nm and 805 nm, with a high spectral (4 nm) and spatial resolution (1 mm per image pixel). The study area was located in a manganese nodule field in the Peru Basin (SE Pacific), close to the DISCOL (DISturbance and reCOLonization) experimental area. To test whether underwater hyperspectral imaging can be used for detection and mapping of mineral deposits in potential deep-sea mining areas, we compared two supervised classification methods, the Support Vector Machine (SVM) and the Spectral Angle Mapper (SAM). The results show that SVM is superior to SAM and is able to accurately detect nodule surfaces. The UHI therefore represents a promising tool for high-resolution seafloor exploration and characterisation prior to resource exploitation.

Description: Highlights: • A High Resolution-LOPC and a FlowCAM were evaluated for ballast water monitoring. • Both instruments underestimated density compared to microscopy. • Size measurements can be affected by organism orientation and complex morphology. • Both tools might be particularly useful when working with a known community. Abstract: Many commercial ships will soon begin to use treatment systems to manage their ballast water and reduce the global transfer of harmful aquatic organisms and pathogens in accordance with upcoming International Maritime Organization regulations. As a result, rapid and accurate automated methods will be needed to monitoring compliance of ships' ballast water. We examined two automated particle counters for monitoring organisms ≥ 50 μm in minimum dimension: a High Resolution Laser Optical Plankton Counter (HR-LOPC), and a Flow Cytometer with digital imaging Microscope (FlowCAM), in comparison to traditional (manual) microscopy considering plankton concentration, size frequency distributions and particle size measurements. The automated tools tended to underestimate particle concentration compared to standard microscopy, but gave similar results in terms of relative abundance of individual taxa. For most taxa, particle size measurements generated by FlowCAM ABD (Area Based Diameter) were more similar to microscope measurements than were those by FlowCAM ESD (Equivalent Spherical Diameter), though there was a mismatch in size estimates for some organisms between the FlowCAM ABD and microscope due to orientation and complex morphology. When a single problematic taxon is very abundant, the resulting size frequency distribution curves can become skewed, as was observed with Asterionella in this study. In particular, special consideration is needed when utilizing automated tools to analyse samples containing colonial species. Re-analysis of the size frequency distributions with the removal of Asterionella from FlowCAM and microscope data resulted in more similar curves across methods with FlowCAM ABD having the best fit compared to the microscope, although microscope concentration estimates were still significantly higher than estimates from the other methods. The results of our study indicate that both automated tools can generate frequency distributions of particles that might be particularly useful if correction factors can be developed for known differences in well-studied aquatic ecosystems.

Description: Highlights • Biogeochemical processes create CO2 sources/sinks by altering seawater AT and CT • Source/sink strength depends on local seawater ‘isocapnic quotient’ (Q) • Q depends on seawater temperature and the state of the marine carbonate system • Spatiotemporal variability in Q drives heterogeneous CO2 source/sink magnitude • Future warming and CO2 emissions will modify Q and the size of CO2 sources/sinks Abstract The ocean holds a large reservoir of carbon dioxide (CO2), and mitigates climate change through uptake of anthropogenic CO2. Fluxes of CO2 between the atmosphere and surface ocean are regulated by a number of physical and biogeochemical processes, resulting in a spatiotemporally heterogeneous CO2 distribution. Determining the influence of each individual process is useful for interpreting marine carbonate system observations, and is also necessary to investigate how changes in these drivers could affect air-sea CO2 exchange. Biogeochemical processes exert an influence primarily through modifying seawater dissolved inorganic carbon (CT) and total alkalinity (AT), thus changing the seawater partial pressure of CO2 (psw). Here, we propose a novel conceptual framework through which the size of the CO2 source or sink generated by any biogeochemical process, denoted Φ, can be evaluated. This is based on the ‘isocapnic quotient’ (Q), which defines the trajectory through (AT,CT) phase space for which there is no change in psw. We discuss the limitations and uncertainties inherent in this technique, which are negligible for most practical purposes, and its links with existing, related approaches. We investigate the effect on Φ of spatiotemporal heterogeneity in Q in the present day surface ocean for several key biogeochemical processes. This leads the magnitude of the CO2 source or sink generated by processes that modify AT to vary spatiotemporally. Finally, we consider how the strength of each process as a CO2 source or sink may change in a warmer, higher-CO2 future ocean.

Description: Highlights • Complete upper Albian to early Turonian climate archive in drilled core from Tarfaya Basin. • Eccentricity pacing of mid Cretaceous OAE isotope excursions. • MCE and OAE2 associated with climate cooling and sea level fall. Abstract A 325 m long continuous succession of uppermost Albian to lower Turonian pelagic (outer shelf) deposits was recovered from a new drill site in the central part of the Tarfaya Basin (southern Morocco). Natural gamma ray wireline logging, carbonate and organic carbon content, bulk carbonate and organic carbon stable isotopes and X-ray fluorescence (XRF)-scanner derived elemental distribution data in combination with planktonic foraminiferal biostratigraphy indicate complete recovery of the Cenomanian Stage. This exceptional sediment archive allows to identify orbitally driven cyclic sedimentation patterns and to evaluate the pacing of climatic events and regional environmental change across the Albian-Cenomanian boundary (ACB), the mid-Cenomanian Event (MCE) and Oceanic Anoxic Event 2 (OAE2) in the latest Cenomanian. The deposition of organic-rich sediments in the Tarfaya Basin, likely driven by upwelling of nutrient-rich water masses, started during the latest Albian and intensified in two major steps following the MCE and the onset of OAE2. The duration and structure of the MCE and OAE2 carbon isotope excursions exhibit striking similarities, suggesting common driving mechanisms and climate-carbon cycle feedbacks. Both events were also associated with eustatic sea level falls, expressed as prominent sequence boundaries in the Tarfaya Basin. Based on the 405 kyr signal imprinted on the Natural Gamma Ray (NGR) and XRF-scanner derived Log(Zr/Rb) records, we estimate the duration of the Cenomanian Stage to be 4.8 ± 0.2 Myr.

Description: Highlights • Pb concentrations and isotope ratios presented for GEOTRACES section GA06. • Northern and southern hemisphere water masses have distinct Pb isotope ratios. • Pb isotope ratios consistent with ventilation timescales of northern water masses. • Mixing complicates interpretation of Pb distributions in southern water masses. Abstract Anthropogenic emissions have dominated marine Pb sources during the past century. Here we present Pb concentrations and isotope compositions for ocean depth profiles collected in the eastern Tropical Atlantic Ocean (GEOTRACES section GA06), to trace the transfer of anthropogenic Pb into the ocean interior. Variations in Pb concentration and isotope composition were associated with changes in hydrography. Water masses ventilated in the southern hemisphere generally featured lower 206Pb/207Pb and 208Pb/207Pb ratios than those ventilated in the northern hemisphere, in accordance with Pb isotope data of historic anthropogenic Pb emissions. The distributions of Pb concentrations and isotope compositions in northern sourced waters were consistent with differences in their ventilation timescales. For example, a Pb concentration maximum at intermediate depth (600–900 m, 35 pmol kg−1) in waters sourced from the Irminger/Labrador Seas, is associated with Pb isotope compositions (206Pb/207Pb = 1.1818–1.1824, 208Pb/207Pb = 2.4472–2.4483) indicative of northern hemispheric emissions during the 1950s and 1960s close to peak leaded petrol usage, and a transit time of ∼50–60 years. In contrast, North Atlantic Deep Water (2000–4000 m water depth) featured lower Pb concentrations and isotope compositions (206Pb/207Pb = 1.1762–1.184, 208Pb/207Pb = 2.4482–2.4545) indicative of northern hemispheric emissions during the 1910s and 1930s and a transit time of ∼80–100 years. This supports the notion that transient anthropogenic Pb inputs are predominantly transferred into the ocean interior by water mass transport. However, the interpretation of Pb concentration and isotope composition distributions in terms of ventilation timescales and pathways is complicated by (1) the chemical reactivity of Pb in the ocean, and (2) mixing of waters ventilated during different time periods. The complex effects of water mass mixing on Pb distributions is particularly apparent in seawater in the Tropical Atlantic Ocean which is ventilated from the southern hemisphere. In particular, South Atlantic Central Water and Antarctic Intermediate Water were dominated by anthropogenic Pb emitted during the last 50–100 years, despite estimates of much older average ventilation ages in this region.

Description: The larger area of and around the Early Iron Age fortress Heuneburg is focus of decades of archaeological excavations and observations. Additionally, to deciduous oak (Quercus sp.) known as the main timber during this period (Middle to Late Hallstatt c. 750–450 BCE), silver fir (Abies alba) was found recently on the Heuneburg plateau itself and its surroundings. Silver fir is now recognized being a significant source of timber during the Iron Age, in a region where its nowadays natural status is debated and its occurrence is rare. The aim of this study was to determine the possible source of the used timber, which might have been taken from the some 80 km distant Abies-rich Black Forest on the primary bed rock or from stands developed on younger geological formations nearer to the site. For this approach, radiogenic strontium isotope (87Sr/86Sr) analyses were performed on waterlogged wood and on charcoal remains buried in four archaeological contexts from Late Hallstatt period around the princely seat Heuneburg. The 87Sr/86Sr ratios were compared to those of living trees, soils and sediments with respect to the diverse geological background. The geological ground of the archaeological wood was narrowed down to molasse for three structures and to loess, moraine or mass chalk sites for one. The isotopic ratio thus points to growing sites more in the surrounding of the concerned archaeological sites and a more distant Black Forest source can be ruled out. Local evidence of Abies was additionally supported by palynological data from four archaeological sites and by anthracological data from a funeral context near the Heuneburg.

Description: Highlights • The genetic model for Algoma-type BIF is modified taking into account S-MIF results. • Metal and sulfur sources are decoupled and reflect diverse microbial metabolisms. • Sulfur deposited with oxide-facies BIF is mostly atmospheric in origin. • Little juvenile sulfur is found, despite the proximity to volcanic sources. Abstract Neoarchean Algoma-type banded iron formations (BIFs) are widely viewed as direct chemical precipitates from proximal volcanic–hydrothermal vents. However, a systematic multiple sulfur isotope study of oxide-facies BIF from a type locality in the ca. 2.74 Ga Temagami greenstone belt reveals mainly bacterial turnover of atmospheric elemental sulfur in the host basin rather than deposition of hydrothermally cycled seawater sulfate or sulfur from direct volcanic input. Trace amounts of chromium reducible sulfur that were extracted for quadruple sulfur isotope (32S–33S–34S–36S) analysis record the previously known mass-independent fractionation of volcanic SO2 in the Archean atmosphere (S-MIF) and biological sulfur cycling but only minor contributions from juvenile sulfur, despite the proximity of volcanic sources. We show that the dominant bacterial metabolisms were iron reduction and sulfur disproportionation, and not sulfate reduction, consistent with limited availability of organic matter and the abundant ferric iron deposited as Fe(OH)3. That sulfur contained in the BIF was not a direct volcanic–hydrothermal input, as expected, changes the view of an important archive of the Neoarchean sulfur cycle in which the available sulfur pools were strongly decoupled and only species produced photochemically under anoxic atmospheric conditions were deposited in the BIF-forming environment.

Description: We present a dataset of dissolved methane (CH4) in the East China Sea (ECS) during five cruises in March, May, August, October and December 2011. CH4 distribution in this region showed pronounced spatial and seasonal variability due to the complex mixing of different water masses and other variables. Surface CH4 concentrations gradually decreased from the coast to the open sea, with maximum values occurring near Changjiang estuary or outside the Hangzhou Bay. The annual mean CH4 concentration of the surface layer was 9.1 ± 1.6 nmol L-1 in the coastal area, which was nearly twice as large as that in the open sea (4.3 ± 1.3 nmol L-1). CH4 was distributed evenly from the surface to the bottom in the shelf region during March and December, while it increased gradually with depth during May and October. CH4 depth profiles exhibited various distribution features along the coast, in the middle and on the edge of continental shelf. CH4 levels at the bottom were generally higher than at the surface during all seasons, indicating obvious CH4 sources from sediments. Incubation experiments of sediment cores onboard suggested that sediment release was an important source of CH4 in the water column of the ECS. We estimated that the annual average CH4 release rate from sediments was about 1.11 μmol·m-2·d-1 on the continental shelf of the ECS. The maximum CH4 concentration and sediment-water CH4 flux both occurred in summer, which might be related with the occurrence of hypoxia in the bottom water. Surface seawater of the ECS was oversaturated with CH4 relative to the atmosphere over most of the five cruises, indicating that the ECS was a net source of atmospheric CH4. The annual mean area-weighted sea-air flux density of CH4 in the ECS was estimated to be about 10.7 μmol·m-2·d-1 in 2011. Accordingly, an area-weighted, seasonally adjusted annual rate of CH4 efflux was determined to be 2.98 × 109 mol yr-1 (∼0.05 Tg CH4 yr-1) from the ECS to the atmosphere.

Description: Highlights • Elongated fault structures are conduits for focused fluid flow. • Gas migration occurs only along a sub-set of faults across Opouawe bank. • Stress state deduced from 3D fault structures appears partially stratigraphically controlled. Abstract High-resolution 2D and 3D seismic data from Opouawe Bank, an accretionary ridge on the Hikurangi subduction margin off New Zealand, show evidence for exceptional gas migration pathways linked to the stress regime of the ridge. Although the ridge has formed by thrusting and folding in response to a sub-horizontal principal compressive stress (σ1), it is clear that local stress conditions related to uplift and extension around the apex of folding (i.e. sub-vertical σ1) are controlling shallow fluid flow. The most conspicuous structural features are parallel and horizontally-elongated extensional fractures that are perpendicular to the ridge axis. At shallower depth near the seafloor, extensional fractures evolve into more concentric structures which ultimately reach the seafloor where they terminate at gas seeps. In addition to the ridge-perpendicular extensional fractures, we also observe both ridge-perpendicular and ridge-parallel normal faults. This indicates that both longitudinal- and ridge-perpendicular extension have occurred in the past. The deepest stratigraphic unit that we image has undergone significant folding and is affected by both sets of normal faults. Shallower stratigraphic units are less deformed and only host the ridge-parallel normal faults, indicating that longitudinal extension was limited to an older phase of ridge evolution. Present-day gas migration has exploited the fabric from longitudinal extension at depth. As the gas ascends to shallower units it ‘self-generates’ its flow pathways through the more concentric structures near the seafloor. This shows that gas migration can evolve from being dependent on inherited tectonic structures at depth, to becoming self-propagating closer to the seafloor.

Description: Highlights • We reconstructed variation in nutrient utilization over the Laptev Sea throughout the Holocene. • The Holocene Siberian transgression modulated the water column structure and created unstable conditions until 4 ka. • Oceanographic conditions favourable to the onset of the Laptev Sea ‘sea-ice factory’ were reached around 2 ka. Abstract Understanding the dynamic of freshwater and sea-ice export from the Arctic is crucial to better comprehend the potential near-future climate change consequences. Here, we report nitrogen isotope data of a core from the Laptev Sea to shed light on the impact of the Holocene Siberian transgression on the summer stratification of the Laptev Sea. Our data suggest that the oceanographic setting was less favourable to sea-ice formation in the Laptev Sea during the early to mid-Holocene. It is only after the sea level reached a standstill at around 4 ka that the water column structure in the Laptev Sea became more stable. Modern-day conditions, often described as “sea-ice factory”, were reached about 2 ka ago, after the development of a strong summer stratification. These results are consistent with sea-ice reconstruction along the Transpolar Drift, highlighting the potential contribution of the Laptev Sea to the export of freshwater from the Arctic Ocean.

Description: Highlights • Epi melts have experienced no disequilibrium modification by mixing or assimilation • Melts fractionate continuously while ascending, rather than stagnating • Magma ascent is through a complex system of dykes and sills • Epi situated between compressional and extensional regime on thick island arc crust • Structural features have impact on focusing and composition of island arc magmas Abstract We present here new bathymetric, petrological and geochemical whole rock, glass and mineral data from the submarine Epi volcano in the New Hebrides (Vanuatu) island arc. The structure has previously been interpreted to be part of a larger caldera structure but new bathymetric data reveal that the volcanic cones are aligned along shear zones controlled by the local tectonic stress field parallel to the recent direction of subduction. We aim to test if there is an interaction between local tectonics and magmatism and to what extent the compositions of island arc volcanoes may be influenced by their tectonic setting. Primitive submarine Epi lavas and those from the neighbouring Lopevi and Ambrym islands originate from a depleted mantle wedge modified by addition of subduction zone components. Incompatible element ratios sensitive to fluid input (e.g., Th/Nb, Ce/Yb) in the lavas are positively correlated with those more sensitive to mantle wedge depletion (e.g., Nb/Yb, Zr/Nb) amongst the arc volcanoes suggesting that fluids or melts from the subducting sediments have a stronger impact on the more depleted compositions of the mantle wedge. The whole rock, glass and mineral major and trace element compositions and the occurrence of exclusively normally zoned clinopyroxene and plagioclase crystals combined with the absence of inversely zoned crystals and water-bearing phases in both mafic and evolved lavas suggest that the erupted melt was relatively dry compared to other subduction zone melts and has experienced little disequilibrium modification by melt mixing or assimilation. Our data also imply that differentiation of amphibole is not required to explain the incompatible element patterns but may rather result from extensive clinopyroxene fractionation in agreement with petrographic observations. Thermobarometric calculations indicate that the melts fractionated continuously during ascent, contrasting with fractionation during stagnation in an established crustal magma reservoir. We interpret the occurrence of this fractional crystallisation end-member in a relatively thick island arc crust (~30 km thickness) to result from isolated and relatively rapid ascent of melts, most likely through a complex system of dykes and sills that developed due to the tectonic positioning of Epi in a complex tectonic zone between a compressional environment in the north and an extensional setting in the south. We can show that the alignment of the cones largely depends on the local tectonic stress field at Epi that is especially influenced by a large dextral strike-slip zone, indicating that structural features have a significant impact on the location and composition of volcanic edifices.

Description: Since January 2014, the reformed Common Fisheries Policy (CFP) of the European Union is legally binding for all Member States. It prescribes the end of overfishing and the rebuilding of all stocks above levels that can produce maximum sustainable yields (MSY). This study examines the current status, exploitation pattern, required time for rebuilding, future catch, and future profitability for 397 European stocks. Fishing pressure and biomass were estimated from 2000 to the last year with available data in 10 European ecoregions and 2 wide ranging regions. In the last year with available data, 69% of the 397 stocks were subject to ongoing overfishing and 51% of the stocks were outside of safe biological limits. Only 12% of the stocks fulfilled the prescriptions of the CFP. Fishing pressure has decreased since 2000 in some ecoregions but not in others. Barents Sea and Norwegian Sea have the highest percentage (〉60%) of sustainably exploited stocks that are capable of producing MSY. In contrast, in the Mediterranean Sea, fewer than 20% of the stocks are exploited sustainably. Overfishing is still widespread in European waters and current management, which aims at maximum sustainable exploitation, is unable to rebuild the depleted stocks and results in poor profitability. This study examines four future exploitation scenarios that are compatible with the CFP. It finds that exploitation levels of 50–80% of the maximum will rebuild stocks and lead to higher catches than currently obtained, with substantially higher profits for the fishers.

Description: Highlights • Lusi crater waters represent a regional geochemical anomaly. • Erupted waters are the result of a complex mix of sedimentary and hydrothermal fluids. • Lusi is not a typical mud volcano but a sediment-hosted hydrothermal system. • The neighbouring volcanic complex feeds hydrothermal fluids for the Lusi eruption. Abstract The spectacular Lusi mud-eruption started in northeast Java the 29th of May 2006. Despite extensive research, the origin of the erupted water remains elusive and poorly constrained. Here we present a comprehensive study of the geochemistry of Lusi waters compared with those collected from surrounding areas, all collected between 2006 and 2013, including data from mud volcanoes and volcano-hosted hydrothermal springs. Within this broad context, the geochemical characteristics of the fluids expelled in the Lusi region suggest that we can classify the waters in three groups: 1) meteoric waters expelled in cold springs and artesian wells, 2) hydrothermal waters typically mixed with meteoric waters, and 3) formation water from marine sediments altered by diagenetic processes such as clay-mineral dehydration. Samples collected from the Lusi crater are Cl and Na dominated (up to 527 mM and 471.7 mM, respectively) similar to seawater indicating that altered sedimentary formation waters are predominant in this system. In addition they are enriched in Sr (up to 808.4 μM) and other elements commonly associated with hydrothermal systems, such as Li (up to 877.6 μM compared to 26 μM in seawater). Some of these elements are up to ten times enriched compared to seawater values. High-temperature fluid mineral interactions in the subsurface appear to have facilitated the transfer of Li and other mobile elements into the fluids. High temperature fluid-mineral interaction reactions are also supported by Si concentrations significantly higher compared to other sampled mud volcanoes in the island. Crater samples also show the highest δ18O values (+5‰ after correction for evaporation compared to +1‰ at the MV localities). 87Sr/86Sr ratios vary between of 0.7077 and 0.7083 and seem to reflect a general mixture of fluids from clay-mineral dehydration, carbonate recrystallization, alteration of volcanic rocks and hydrothermal imprint. Eight years of geochemical monitoring indicate that the composition of the deep-sourced Lusi fluids remain fairly constant through time. Thus our findings show that the Lusi crater waters represent a regional geochemical anomaly, and we suggest that a combination of high temperatures in the source region, and fluid-rock interactions with silicates and, possibly, carbonate-rich lithologies can explain the data. This is consistent with a model where the emitted gases migrate from a deep-seated (〉4 km) source region, likely associated with the presence of hot igneous intrusions and/or high T reactions related to the presence of neighbouring active volcanoes.

Description: Highlights • Cryptotephra study of a Holocene sedimentary record from the Chukchi Sea. • Major tephra concentration peak fingerprinted to the ∼3.6 ka Aniakchak eruption. • New electron microprobe and LA-ICP-MS glass data applicable for the Western Arctic. • Re-evaluation of the Aniakchak tephra volume. • Redeposited tephra shards map pathways of sediment transport. Abstract Developing chronologies for sediments in the Arctic Ocean and its continental margins is an important but challenging task. Tephrochronology is a promising tool for independent age control for Arctic marine sediments and here we present the results of a cryptotephra study of a Holocene sedimentary record from the Chukchi Sea. Volcanic glass shards were identified and quantified in sediment core HLY0501-01 and geochemically characterized with single-shard electron microprobe and laser ablation-inductively coupled plasma-mass spectrometry (LA-ICP-MS). This enabled us to reveal a continuous presence of glass shards with identifiable chemical compositions throughout the core. The major input of glasses into the sediments is geochemically fingerprinted to the ∼3.6 ka Aniakchak caldera II eruption (Alaska), which provides an important chronostratigraphic constraint for Holocene marine deposits in the Chukchi-Alaskan region and, potentially, farther away in the western Arctic Ocean. New findings of the Aniakchak II tephra permit a reevaluation of the eruption size and highlight the importance of this tephra as a hemispheric late Holocene marker. Other identified glasses likely originate from the late Pleistocene Dawson and Old Crow tephras while some cannot be correlated to certain eruptions. These are present in most of the analyzed samples, and form a continuous low-concentration background throughout the investigated record. A large proportion of these glasses are likely to have been reworked and brought to the depositional site by currents or other transportation agents, such as sea ice. Overall, our results demonstrate the potential for tephrochronology for improving and developing chronologies for Arctic Ocean marine records, however, at some sites reworking and redistribution of tephra may have a strong impact on the record of primary tephra deposition.

Description: Thyroid hormones (THs) are key regulators of growth, development, and metabolism in vertebrates and influence early life development of fish. TH is produced in the thyroid gland (or thyroid follicles) mainly as T4 (thyroxine), which is metabolized to T3 (3,5,3'-triiodothyronine) and T2 (3,5-diiodothyronine) by deiodinase (DIO) enzymes in peripheral tissues. The action of these hormones is mostly exerted by binding to a specific nuclear thyroid hormone receptor (THR). In this study, we i) cloned and characterized thr sequences, ii) investigated the expression pattern of the different subtypes of thrs and dios, and iii) studied how temperature affects the expression of those genes in artificially produced early life history stages of European eel (Anguilla anguilla), reared in different thermal regimes (16, 18, 20 and 22 °C) from hatch until first-feeding. We identified 2 subtypes of thr (thrα and thrβ) with 2 isoforms each (thrαA, thrαB, thrβA, thrβB) and 3 subtypes of deiodinases (dio1, dio2, dio3). All thr genes identified showed high similarity to the closely related Japanese eel (Anguilla japonica). We found that all genes investigated in this study were affected by larval age (in real time or at specific developmental stages), temperature, and/or their interaction. More specifically, the warmer the temperature the earlier the expression response of a specific target gene. In real time, the expression profiles appeared very similar and only shifted with temperature. In developmental time, gene expression of all genes differed across selected developmental stages, such as at hatch, during teeth formation or at first-feeding. Thus, we demonstrate that thrs and dios show sensitivity to temperature and are involved in and during early life development of European eel.

Description: In the deep ocean, fluxes of particulate organic carbon (POC) and calcium carbonate are positively correlated, suggesting that CaCO3 could increase sinking particle densities and/or protect the organic matter from degradation by prokaryotes, the so called “ballast effect”. Here, we used the PArticle Sinking Simulator (PASS) system to investigate the effect of increasing pressure on the biodegradation of calcifying Emiliania huxleyi aggregates. Incubations were carried out over a period of 10 days, simulating the changes in temperature and pressure in the water column of the NW Mediterranean Sea. Aggregates sinking from a depth of 200 m to 1700 m (assuming an average sinking velocity of 150 m d−1) were exposed to a natural mesopelagic prokaryotic community collected from 200 m. In contrast to previous studies, where silicifying diatom aggregates were used, the calcifying E. huxleyi aggregates were found to be more sensitive to degradation with increasing hydrostatic pressure (relative to constant atmospheric pressure). This was confirmed by changes in lipid composition which suggested increased cell lysis. Changes in particulate inorganic carbon and total alkalinity indicated that CaCO3 dissolution might have been faster under pressure. Increased hydrostatic pressure also had a positive effect on particle aggregation, which may compensate for the effect of increased cell lysis. Our results imply that in coccolithophorid-dominated sinking aggregates, the ballasting and protection effects of coccoliths may collapse throughout the water column. The increased aggregation potential with pressure observed in these controlled conditions, may balance the loss of mineral ballast to a certain extent, although this needs to be confirmed in situ.

Description: The impact of storm surge on a barrier island tends to be considered from a single cross-shore dimension, dependent on the relative elevations of the storm surge and dune crest. However, the foredune is rarely uniform and can exhibit considerable variation in height and width at a range of length scales. In this study, LiDAR data from barrier islands in Texas and Florida are used to explore how shoreline position and dune morphology vary alongshore, and to determine how this variability is altered or reinforced by storms and post-storm recovery. Wavelet analysis reveals that a power law can approximate historical shoreline change across all scales, but that stormscale shoreline change (~10 years) and dune height exhibit similar scale-dependent variations at swash and surf zone scales (〈1000 m). The in-phase nature of the relationship between dune height and storm-scale shoreline change indicates that areas of greater storm-scale shoreline retreat are associated with areas of smaller dunes. It is argued that the decoupling of storm-scale and historical shoreline change at swash and surf zone scales is also associated with the alongshore redistribution of sediment and the tendency of shorelines to evolve to a more diffusive (or straight) pattern with time. The wavelet analysis of the data for post-storm dune recovery is also characterized by red noise at the smallest scales characteristic of diffusive systems, suggesting that it is possible that small-scale variations in dune height can be repaired through alongshore recovery and expansion if there is sufficient time between storms. However, the time required for dune recovery exceeds the time between storms capable of eroding and overwashing the dune. Correlation between historical shoreline retreat and the variance of the dune at swash and surf zone scales suggests that the persistence of the dune is an important control on transgression through island migration or shoreline retreat with relative sea-level rise.

Description: Highlights: • Optimality-based modelling of microzooplankton trophic interactions in mesocsosm ecosystems. • Intraguild predation is important for modelling microzooplankton feeding interactions. • Trophic interactions structured solely by size may fail to capture feeding diversity. • Adequate representation of feeding interaction is needed for modelling ecosystem dynamics. Abstract: The zooplankton components in biogeochemical models drive top-down control of primary production and remineralisation, and thereby exert a strong impact on model performance. Who eats whom in oceanic plankton ecosystem models is often largely determined by body size. However, zooplankton of similar size can have different prey-size spectra. Thus, models with solely size-structured trophic interactions may not capture the full diversity of feeding interactions and miss important parts of zooplankton behavior. We apply an optimality-based plankton ecosystem model to analyse trophic interactions in a suite of mesocosm experiments in the Peruvian upwelling region. Sensitivity analyses reveal a dominant role of trophic structure for model performance, which cannot be compensated by parameter optimisation. The single most important aspect governing model performance is the trophic linking between dinoflagellates and ciliates. Only with a bidirectional link, i.e., both groups can prey on each other, is the model able to reproduce the differential development of the microzooplankton communities in the mesocosms. Thus, we conclude that a solely size-based trophic structure may not be appropriate to represent the most important trophic interactions in plankton ecosystems. The diversity of feeding interactions needs to be adequately represented to capture community dynamics.

Description: Highlights • Bulk rock elastic moduli of gneiss, amphibolite and marble using different techniques. • Neutron diffraction texture analysis and modeling of rock physical properties. • Measurement of seismic velocity anisotropy under pressures of up to 600 MPa. • Extrapolation of experimental data to higher pressures of 1000 MPa (crack free rock). • Comparison of modeled, experimental and extrapolated elastic anisotropy data. Abstract In this study elastic moduli of three different rock types of simple (calcite marble) and more complex (amphibolite, micaschist) mineralogical compositions were determined by modeling of elastic moduli using texture (crystallographic preferred orientation; CPO) data, experimental investigation and extrapolation. 3D models were calculated using single crystal elastic moduli, and CPO measured using time-of-flight neutron diffraction at the SKAT diffractometer in Dubna (Russia) and subsequently analyzed using Rietveld Texture Analysis. To define extrinsic factors influencing elastic behaviour, P-wave and S-wave velocity anisotropies were experimentally determined at 200, 400 and 600 MPa confining pressure. Functions describing variations of the elastic moduli with confining pressure were then used to predict elastic properties at 1000 MPa, revealing anisotropies in a supposedly crack-free medium. In the calcite marble elastic anisotropy is dominated by the CPO. Velocities continuously increase, while anisotropies decrease from measured, over extrapolated to CPO derived data. Differences in velocity patterns with sample orientation suggest that the foliation forms an important mechanical anisotropy. The amphibolite sample shows similar magnitudes of extrapolated and CPO derived velocities, however the pattern of CPO derived velocity is closer to that measured at 200 MPa. Anisotropy decreases from the extrapolated to the CPO derived data. In the micaschist, velocities are higher and anisotropies are lower in the extrapolated data, in comparison to the data from measurements at lower pressures. Generally our results show that predictions for the elastic behavior of rocks at great depths are possible based on experimental data and those computed from CPO. The elastic properties of the lower crust can, thus, be characterized with an improved degree of confidence using extrapolations. Anisotropically distributed spherical micro-pores are likely to be preserved, affecting seismic velocity distributions. Compositional variations in the polyphase rock samples do not significantly change the velocity patterns, allowing the use of RTA-derived volume percentages for the modeling of elastic moduli.

Description: Highlights • Barrier island beach and dunes can inherit morphology from framework geology. • The influence of buried paleo-channels is affected by their scale and orientation. • Non-linear interaction occurs between geology and smaller-scale coastal processes. • Framework geology is key to predicting barrier island response to sea level rise. • Fine-scale variability in dune height is greater landward of buried paleo-channels. Abstract Barrier island response and recovery to storms, and island transgression with relative sea level rise, can be influenced by the framework geology. The influence of framework geology on barrier island geomorphology has previously been examined in areas where the framework is rhythmic alongshore or consists of an isolated paleo-channel or headland. The purpose of this paper is to examine the influence of framework geology on beach and dune geomorphology at Padre Island National Seashore (PAIS), Texas, USA, where the framework geology is variable alongshore. Alongshore beach and dune morphometrics and offshore bathymetric profiles were extracted from a combined topography and bathymetry digital elevation model (DEM) using an automated approach along the ~100 km study area, and an electromagnetic induction (EMI) survey was used to map the subsurface framework geology. Wavelet decomposition, Global Wavelet (GW), and bicoherence analyses were used to test for spatial relationships between and within the extracted alongshore metrics. GW trendlines demonstrate that beach and dune morphometrics are structurally controlled. Hotspots in wavelet coherence plots between framework geology and alongshore island morphometrics indicate that the paleo-channels dissecting the island influence beach and dune morphology, with large dunes found in the area directly landward of the paleochannels. Bicoherence analysis of alongshore beach and dune morphometrics indicates that low-frequency oscillations due to framework geology interact with higher-frequency oscillations, with greater small-scale variability in the dune line directly landward of the paleo-channels. These results suggest that the paleo-channels of PAIS non-linearly influence beach and dune morphology, which in turn alters the response of the island to storms and sea level rise. It is argued that an understanding of the framework geology is key to predicting island response to sea level rise and framework geology needs to be included in barrier island models. This paper demonstrates that an irregular framework geology influences small-scale coastal processes, and creates interactions across scales that influence beach and dune morphology and affects barrier island response to storms and sea level rise.

Description: Particle aggregation and the consequent formation of marine snow alter important properties of biogenic particles (size, sinking rate, degradability), thus playing a key role in controlling the vertical flux of organic matter to the deep ocean. However, there are still large uncertainties about rates and mechanisms of particle aggregation, as well as the role of plankton community structure in modifying biomass transfer from small particles to large fast-sinking aggregates.Here we present data from a high-resolution underwater camera system that we used to observe particle size distributions and formation of marine snow (aggregates 〉0.5 mm) over the course of a 9-week in situ mesocosm experiment in the Eastern Subtropical North Atlantic. After an oligotrophic phase of almost 4 weeks, addition of nutrient-rich deep water (650 m) initiated the development of a pronounced diatom bloom and the subsequent formation of large marine snow aggregates in all 8 mesocosms. We observed a substantial time lag between the peaks of chlorophyll a and marine snow biovolume of 9-12 days, which is much longer than previously reported and indicates a marked temporal decoupling of phytoplankton growth and marine snow formation during our study. Despite this time lag, our observations revealed substantial transfer of biomass from small particle sizes (single phytoplankton cells and chains) to marine snow aggregates of up to 2.5 mm diameter (ESD), with most of the biovolume being contained in the 0.5-1 mm size range. Notably, the abundance and community composition of mesozooplankton had a substantial influence on the temporal development of particle size spectra and formation of marine snow aggregates: While higher copepod abundances were related to reduced aggregate formation and biomass transfer towards larger particle sizes, the presence of appendicularia and doliolids enhanced formation of large marine snow.Furthermore, we combined in situ particle size distributions with measurements of particle sinking velocity to compute instantaneous (potential) vertical mass flux. However, somewhat surprisingly, we did not find a coherent relationship between our computed flux and measured vertical mass flux (collected by sediment traps in 15 m depth). Although the onset of measured vertical flux roughly coincided with the emergence of marine snow, we found substantial variability in mass flux among mesocosms that was not related to marine snow numbers, and was instead presumably driven by zooplankton-mediated alteration of sinking biomass and export of small particles (fecal pellets).Altogether, our findings highlight the role of zooplankton community composition and feeding interactions on particle size spectra and formation of marine snow aggregates, with important implications for our understanding of particle aggregation and vertical flux of organic matter in the ocean.

Description: Primitive olivine-hosted melt inclusions provide information concerning the pre-eruptive volatile contents of silicate melts, but compositional changes associated with post-entrapment processes (PEP) sometimes complicate their interpretation. In particular, crystallization of the host phase along the wall of the melt inclusion and diffusion of H+ through the host promote CO2 and potentially S or other volatiles to exsolve from the melt into a separate fluid phase. Experimental rehomogenization and analysis of MI, or a combination of Raman spectroscopy, numerical modeling, and mass balance calculations are potentially effective methods to account for PEP and restore the original volatile contents of melt inclusions. In order to compare these different approaches, we studied melt inclusions from a suite of samples from Klyuchevskoy volcano (Kamchatka Arc) for which volatile compositions have been determined using experimental rehydration, Raman spectroscopy, and numerical modeling. The maximum CO2 contents of melt inclusions are in agreement (~3600-4000 ppm), regardless of the method used to correct for CO2 in the bubble, but significantly more uncertainty is observed using mass balance calculations. This uncertainty is largely due to the lack of precision associated with the petrographic method of determining bubble volumes and may also be related to the presence of daughter minerals at the glass-bubble interface.

Description: Highlights: • Ocean acidification increases phytoplankton standing stock. • This increase is more pronounced in smaller-sized taxa. • Primary consumers reac differently depending on nutrient availability. • Bacteria and micro-heterotrophs benefited under limiting conditions. • In general, heterotrophs are negatively affected at nutrient replete periods. Abstract: In situ mesocosm experiments on the effect of ocean acidification (OA) are an important tool for investigating potential OA-induced changes in natural plankton communities. In this study we combined results from various in-situ mesocosm studies in two different ocean regions (Arctic and temperate waters) to reveal general patterns of plankton community shifts in response to OA and how these changes are modulated by inorganic nutrient availability. Overall, simulated OA caused an increase in phytoplankton standing stock, which was more pronounced in smaller-sized taxa. This effect on primary producers was channelled differently into heterotroph primary consumers depending on the inorganic nutrient availability. Under limiting conditions, bacteria and micro-heterotrophs benefited with inconsistent responses of larger heterotrophs. During nutrient replete periods, heterotrophs were in general negatively affected, although there was an increase of some mesozooplankton developmental stages (i.e. copepodites). We hypothesize that changes in phytoplankton size distribution and community composition could be responsible for these food web responses.

Description: Highlights: • Identify 3 groups of gas migration structures in seismic data from the Danube Fan. • Migration structures related to shallow gas migration and flares at the seafloor. • Gas migration is controlled by lithological heterogeneity and sediment deformation. • Mass transport deposits play a role in controlling vertical migration occurrence. Abstract: A large geophysical dataset, including bathymetry, and 2D and 3D P-cable seismic data, revealed evidence of numerous gas flares near the S2 Canyon in the Danube Fan, northwestern Black Sea. This dataset allows us to investigate potential relationships between gas migration pathways, gas vents observed at the seafloor and submarine slope failures. Vertical gas migration structures as revealed in the seismics appear to be concentrated near submarine slope failure structures. Where these seismically defined features extend upwards to the seafloor, they correlate with the location of gas flares. However, not all these structures reach the seafloor, in some cases because they are capped by overlying sediments. A strong correlation is inferred between gas migration pathways, heterogeneous mass transport deposits and contacts between adjacent units of contrasting lithology. Although missing age constrains prevent a final judgement, we discuss the potential relationship between submarine slope failures and gas migration in order to determine if gas migration is a precursor to failure, or if the presence of slope failures and associated mass transport deposits facilitates the migration of gas. Our observations indicate that lithological heterogeneity, mass transport deposits and minor sediment deformation control gas migration pathways and the formation of gas chimney-like features. Gas migration is focused and gradual, resulting in gas flares where the chimney-like features extend to the seafloor, with no evidence of erosive features such as pockmarks.

Description: Oceanic basalts reflect the heterogeneities in the earth's mantle, which can be explained by five mantle end members. The HIMU end member, characterized by high time-integrated μ (238U/204Pb), is defined by the composition of lavas from the ocean islands of St. Helena, South Atlantic Ocean and Mangaia and Tubuai (Cook-Austral Islands), South Pacific Ocean. It is widely considered to be derived from a mantle reservoir that is rarely sampled and not generally involved in mixing with the other mantle components. On the other hand, the FOZO end member, located at the FOcal ZOne of oceanic volcanic rock arrays on isotope diagrams, is considered to be a widespread common component with slightly less radiogenic 206Pb/204Pb and intermediate Sr-Nd-Hf isotopic compositions. Here we present new major and trace element, Sr-Nd-Pb-Hf isotope and geochronological data from the Walvis Ridge and Richardson Seamount in the South Atlantic Ocean and the Manihiki Plateau and Eastern Chatham Rise in the southwest Pacific Ocean. Our new data, combined with literature data, document a more widespread (nearly global) distribution of the HIMU end member than previously postulated. Our survey shows that HIMU is generally associated with low-volume alkaline, carbonatitic and/or kimberlitic intraplate volcanism, consistent with derivation from low degrees of melting of CO2-rich sources. The majority of end member HIMU locations can be directly related to hotspot settings. The restricted trace element and isotopic composition (St. Helena type HIMU), but near-global distribution, point to a deep-seated, widespread reservoir, which most likely formed in the Archean. In this context we re-evaluate the origin of a widespread HIMU reservoir in an Archean geodynamic setting. We point out that the classic ocean crust recycling model cannot be applied in a plume-lid dominated tectonic setting, and instead propose that delamination of carbonatite- metasomatized subcontinental lithospheric mantle could be a suitable HIMU source.

Description: Highlights • First comprehensive data set of the seamounts from the Walvis Ridge. • The seamounts are 20–40 Myr younger than the age progressive Walvis Ridge basement. • The composition of the seamounts extends from the St. Helena HIMU to EMORB. • The seamounts are derived from a distinct source compared to the Walvis Ridge. • The temporal change from EM I to HIMU could reflect the compositional heterogeneities of the LLSVP. Abstract Volcanic activity at many oceanic volcanoes, ridges and plateaus often reawakens after hiatuses of up to several million years. Compared to the earlier magmatic phases, this late-stage (rejuvenated/post-erosional) volcanism is commonly characterized by a distinct geochemical composition. Late-stage volcanism raises two hitherto unanswered questions: Why does volcanism restart after an extended hiatus and what is the origin of this volcanism? Here we present the first 40Ar/39Ar age and comprehensive trace element and Sr–Nd–Pb–Hf isotopic data from seamounts located on and adjacent to the Walvis Ridge in the South Atlantic ocean basin. The Walvis Ridge is the oldest submarine part of the Tristan-Gough hotspot track and is famous as the original type locality for the enriched mantle one (EM I) end member. Consistent with the bathymetric data, the age data indicates that most of these seamounts are 20–40 Myr younger than the underlying or nearby Walvis Ridge basement. The trace element and isotope data reveal a distinct compositional range from the EM I-type basement. The composition of the seamounts extend from the St. Helena HIMU (high time-integrated 238U/204Pb mantle with radiogenic Pb isotope ratios) end member to an enriched (E) Mid-Ocean-Ridge Basalt (MORB) type composition, reflecting a two-component mixing trend on all isotope diagrams. The EMORB end member could have been generated through mixing of Walvis Ridge EM I with normal (N) MORB source mantle, reflecting interaction of Tristan-Gough (EM I-type) plume melts with the upper mantle. The long volcanic quiescence and the HIMU-like geochemical signature of the seamounts are unusual for classical hotspot related late-stage volcanism, indicating that these seamounts are not related to the Tristan-Gough hotspot volcanism. Two volcanic arrays in southwestern Africa (Gibeon-Dicker Willem and Western Cape province) display similar ages to the late-stage Walvis seamounts and also have HIMU-like compositions, suggesting a larger-scale event at ∼77–49 Ma. We propose that the EM I-like mantle plumes rise from the edges of the African Large Low Shear Velocity Province (LLSVP; Tristan-Gough, Discovery and Shona hotspot), whereas the HIMU-dominated intraplate lavas (St. Helena, Gibeon-Dicker Willem and Western Cape province) and the late-stage Walvis seamounts tap material from internal portions of the African LLSVP, suggesting possible lateral and/or vertical chemical zonation of the African LLSVP.

Description: The 541 ± 4 Ma-old magnesian, weakly peraluminous, calc-alkalic Donkerhoek Onanis granite is part of the ca. 6000 km2 large Donkerhoek batholith in the Southern Zone of the Damara orogen of Namibia. Linear major and trace element variations and decreasing MgO, FeO, Al2O3, CaO, K2O, Na2O, Ba and Sr concentrations with increasing SiO2 indicate that this part of the batholith represent a coherent mass and underwent fractional crystallization processes. The Donkerhoek Onanis granites are isotopically evolved (initial εNd: −4.7 to −12.3, initial 87Sr/86Sr: 0.7099–0.7157) with moderately radiogenic Pb isotope ratios (206Pb/204Pb: 17.26–18.22; 207Pb/204Pb: 15.59–15.67; 208Pb/204Pb: 37.60–38.06). Beside heterogeneities imparted by the sources, an evaluation of LREE fractionation and Nd isotope data suggests that AFC processes also modified some samples. Based on the chemical and isotope data, the Donkerhoek Onanis granites cannot be derived by partial melting of Al- and Fe-rich metasedimentary rocks of the Kuiseb formation in which they intruded. Instead, melting of meta-igneous crustal sources with Proterozoic crustal residence ages is more likely. Three igneous to meta-igneous rock suites from the area (Matchless amphibolites, Proterozoic mafic to felsic gneisses from the southern Kalahari craton basement, syn-tectonic Salem granodiorites to granites) are potential sources. An evaluation of chemical and isotope data suggests that remelting of early syn-orogenic Salem-type granites is the most likely process which would also explain the existence of ca. 563 ± 4 Ma-old zircon in the Donkerhoek Onanis granites. Comparison of the Donkerhoek Onanis granites with experimentally derived melt compositions from an intermediate igneous parent indicates temperatures between 800 and 850 °C. It is suggested that the Pan-African igneous activity in this part of the Damara Belt was a moderate-temperature intra-crustal event. Although there are some compositional similarities with juvenile granites generated in subduction zones, unradiogenic Pb isotope ratios and moderately radiogenic Sr and unradiogenic Nd isotopes suggest that reprocessed crustal rocks are more likely sources. Previously obtained high δ18O values of the Donkerhoek Onanis granites ranging from 11.8 to 13.6‰, covering the range of δ18O values obtained on Salem-type granites from the area (12.5–13.3‰) confirm this view. In contrast to igneous processes along active continental margins that produce juvenile batholiths with calc-alkaline affinities, this igneous event was not a major crust-forming episode and the Donkerhoek Onanis granites represent reprocessed crustal material.

Description: Fossil carbonate skeletons of marine organisms are archives for understanding the development and evolution of palaeo-environments. However, the correct assessment of past environment dynamics is only possible when pristine skeletons and their biogenic characteristics are unequivocally distinguishable from diagenetically-altered skeletal elements and non-biogenic features. In this study, we extend our work on diagenesis of biogenic aragonite (Casella et al. 2017) to the investigation of biogenic low-Mg calcite using brachiopod shells. We examined and compared microstructural characteristics induced by laboratory-based alteration to structural features derived from diagenetic alteration in natural environments. We used four screening methods: cathodoluminescence (CL), cryogenic and conventional field emission-scanning electron microscopy (FE-SEM), atomic force microscopy (AFM) and electron backscatter diffraction (EBSD). We base our assessments of diagenetic alteration and overprint on measurements of, a) images of optical overprint signals, b) changes in calcite crystal orientation patterns, and c) crystal co-orientation statistics. According to the screening process, altered and overprinted samples define two groups. In Group 1 the entire shell is diagenetically overprinted, whereas in Group 2 the shell contains pristine as well as overprinted parts. In the case of Group 2 shells, alteration occurred either along the periphery of the shell including the primary layer or at the interior-facing surface of the fibrous/columnar layer. In addition, we observed an important mode of the overprinting process, namely the migration of diagenetic fluids through the endopunctae corroborated by mineral formation and overprinting in their immediate vicinity, while leaving shell parts between endopunctae in pristine condition. Luminescence (CL) and microstructural imaging (FE-SEM) screening give first-order observations of the degree of overprint as they cover macro-to micron scale alteration features. For a comprehensive assessment of diagenetic overprint these screening methods should be complemented by screening techniques such as EBSD and AFM. They visualise diagenetic changes at submicron and nanoscale levels depicting the replacement of pristine nanocomposite mesocrystal biocarbonate (NMB) by inorganic rhombohedral calcite (IRC). The integration of screening methods allows for the unequivocal identification of highly-detailed alteration features as well as an assessment of the degree of diagenetic alteration.

Description: Highlights • Infra red images mosaic was acquired over an erupting clastic geysering system. • Attenuation causes of the IR spectrum are discussed. • Convective dynamics of boiling mud breccia are interpreted inside the caldera. • Pathways for fluid migration are imaged in faulted and fractured zones. Abstract The use of low-cost hand-held infrared (IR) thermal cameras based on uncooled micro-bolometer detector arrays became more widespread during the recent years. Thermal cameras have the ability to estimate temperature values without contact and therefore can be used in conditions where targets are difficult or dangerous to reach such as volcanic eruptions. Since May 2006 the Indonesian Lusi mud eruption continues to spew boiling mud, water, aqueous vapour, CO2, CH4 and covers a surface of nearly 7 km2. Here we performed surveys above and around the erupting crater using a specifically equipped remote-controlled aerial vehicle (drone). Despite the harsh logistics and the continuously varying gas concentrations, we managed to collect IR images composing mosaics to estimate the crater zone spatial and temporal thermal variations as well as that in the surrounding regions. In this manuscript we provide a) a description of the main processes that affect and control the acquisition of IR images; b) an overview of still non disclosed physical model used by the thermal camera employed during our survey (i.e. Flir I7 IR); c) a method for capturing high resolution infrared images over an erupting clastic system; d) analysis and interpretation of the acquired data and scientific results. The results show that it is possible to obtain good quality mosaics also in inaccessible areas such as erupting craters where fixed reference points are not constant, and where the presence of IR attenuation factors introduce errors in terms of temperature estimates. However the IR camera radiative transfer model (based on Lowtran model) allows the control of only some of the parameters that affect the attenuation of the IR spectrum. This results in obtained crater temperature estimates up to ∼20% lower than those measured with hand held thermometers, providing, however, very important information about the thermal gradient and the potential radioactive absorption factors. The imaged Lusi vent thermal pattern suggests the presence of shallow convective chambers inside the caldera zone; these are activated by the rise of boiling mud breccia that suddenly cools and sinks. The thermal survey conducted on the dry mud region to the NE of the crater reveals temperatures matching with those measured directly on the field with a hand held thermometer. Here the presence of hot spot anomalies and colder circular features is consistent with the migration of deeper warm fluids along a faulted and fractured area and with widespread pools distribution.

Description: Highlights • Gases collected from dry seepages and bubbling or dissolved in springs of the ophiolitic body of Kizildag (Turkey). • Large compositional variation (H2- CH4 or N2-dominated) • Hydrogen derives from low temperature serpentinization processes • Methane mainly derives from abiotic processes • Large isotopic fractionation of methane at one site due to biological oxidation Abstract We investigated the geochemical features of the gases released from the Kizildag ophiolitic complex (Hatay, Turkey). Twenty-three samples both dissolved in hyperalkaline waters and free gases (bubbling gases and dry seeps) were collected. Samples were analysed for their chemical (He, H2, O2, N2, CH4 and CO2) and isotopic (He, δ13C-CH4, δ2H-CH4, δ2H-H2) composition including the content and C-isotopic composition of C2 to C5 alkanes in free gases. Analytical results evidence H2 production through low-temperature (〈80 °C) serpentinization processes and subsequent abiogenic CH4 production through Fischer-Tropsch-type reactions. In some sample small additions of methane either of microbial or of thermogenic origin can be hypothesized. At one of the sites (Kisecik) a clear fractionation pattern due to microbial methane oxidation leading to strongly enriched isotopic values (δ13C +15‰ and δ2H −68‰) and depletion in methane concentrations has been evidenced. At the dry gas seep of Kurtbagi methane flux measurements have been made and a preliminary output estimation of about 1000 kg per year has been obtained.

Description: We present sediment pore fluid and sediment solid phase results obtained during IODP Expedition 340 from seven sites located within the Grenada Basin of the southern Lesser Antilles Volcanic Arc region. These sites are generally characterized as being low in organic carbon content and rich in calcium carbonate and volcanogenic material. In addition to the typical reactions related to organic matter diagenesis, pore fluid chemistry indicates that the diagenetic reactions fall within two broad categories; (1) reactions related to chemical exchange with volcanogenic material and (2) reactions related to carbonate dissolution, precipitation, or recrystallization. For locations dominated by reaction with volcanogenic material, these sites exhibit increases in dissolved Ca with coeval decreases in Mg. We interpret this behavior as being driven by sediment-water exchange reactions from the alteration of volcanic material that is dispersed throughout the sediment package, which likely result in formation of Mg-rich secondary authigenic clays. In contrast to this behavior, sediment sequences that exhibit decreases in Ca, Mg, Mn, and Sr with depth suggest that carbonate precipitation is an active diagenetic process affecting solute distributions. The distributions of pore fluid 87Sr/86Sr reflect these competitive diagenetic reactions between volcanic material and carbonate, which are inferred by the major cation distributions. From one site where we have solid phase 87Sr/86Sr (site U1396), the carbonate fraction is found to be generally consistent with the contemporaneous seawater isotope values. However, the 87Sr/86Sr of the non-carbonate fraction ranges from 0.7074 to 0.7052, and these values likely represent a mixture of local arc volcanic sources and trans-Atlantic eolian sources. Even at this site where there is clear evidence for diagenesis of volcanogenic material, carbonate diagenesis appears to buffer pore fluid 87Sr/86Sr from the larger changes that might be expected given the high abundance of tephra in these sediments. Part of this carbonate buffering, at this site as well as throughout the region, derives from the fact that the Sr concentration in the non-carbonate fraction is generally low (〈200 ppm), whereas the carbonate fraction has Sr concentrations approaching ~1000 ppm.

Description: Artificial neural network (ANN) is one of the most widely used methods to develop accurate predictive models based on artificial intelligence and machine learning. In the present study, the important practical aspects of developing a reliable ANN model e.g. appropriate assignment of the number of neurons, number of hidden layers, transfer function, training algorithm, dataset division and initialization of the network are discussed. As a case study, predictability of the flash point for a dataset of 740 organic compounds using ANNs was investigated via a total number of 484220ANNs to allow covering a wide range of parameters affecting the performance of an ANN. Among all studied parameters, the number of neurons or layers was found to be the most important parameters to develop a reliable ANN with low overfitting risk. To evaluate appropriate number of neurons and layers, a value of equal or greater than 10 for the ratio of the training samples to the ANN constants was suggested as a rule of thumb. More ever, a strategy for evaluation of the authentic performance of ANNs and deciding about the reliability of an ANN model was proposed which is applicable to other models developed by supervised learning. Based on the introduced considerations, an ANN model was proposed for predicting the flash point of pure organic compounds. According to the results, the new model was found to produce the lowest error compared to other available models.

Description: Molybdenum (Mo) concentrations and isotope compositions in sediments and shales are commonly used as proxies for anoxic and sulfidic (i.e., euxinic) conditions in the water column of paleo-marine systems. A basic assumption underlying this practice is that the proxy signal extracted from the geological record is controlled by long-term (order of decades to millennia) Mo scavenging in the euxinic water column rather than Mo deposition during brief episodes or events (order of weeks to months). To test whether this assumption is viable we studied the biogeochemical cycling of Mo and its isotopes in sediments of the intermittently euxinic Gotland Deep in the central Baltic Sea. Here, multiannual to decadal periods of euxinia are occasionally interrupted by inflow events during which well‑oxygenated water from the North Sea penetrates into the basin. During these events manganese (Mn) (oxyhydr)oxide minerals are precipitated in the water column, which are known to scavenge Mo. We present sediment and pore water Mo and Mo isotope data for sediment cores which were taken before and after a series of inflow events between 2014 and 2016. After seawater inflow, pore water Mo concentrations in anoxic surface sediments exceed the salinity-normalized concentration by more than two orders of magnitude and coincide with transient peaks of dissolved Mn. A fraction of the Mo liberated into the pore water is transported by diffusion in a downward direction and sequestered by organic matter within the sulfidic zone of the sediment. Diffusive flux calculations as well as a mass balance that is based on the sedimentary Mo isotope composition suggest that about equal proportions of the Mo accumulating in the basin are delivered by Mn (oxyhydr)oxide minerals during inflow events and Mo scavenging with hydrogen sulfide during euxinic periods. Since the anoxic surface sediment where Mo is released from Mn (oxyhydr)oxides are separated by several centimeters from the deeper sulfidic layers where Mo is removed, the solid phase record of Mo concentration and isotope composition would be misinterpreted if steady state Mo accumulation was assumed. Based on our observations in the Gotland Deep, we argue that short-term redox fluctuations need to be considered when interpreting Mo-based paleo-records.

Description: To analyze the global hydroclimate response during the Younger Dryas cold event, we evaluate climate model results that have been constrained with proxy-based temperatures from the North Atlantic region. We find that both the temperature and the hydroclimate response have a clear global signature. A marked cooling is simulated over the North Atlantic Ocean (more than 5 °C) and the downwind continents (2–4 °C). This response is related to the weakening of the Atlantic meridional overturning circulation under influence of meltwater discharges. The hydroclimate response is most expressed over Eurasia in a belt between 40 and 60°N, and over Northern Africa in the Sahel region. In both areas, a strong decrease in soil moisture is simulated (up to 20% reduction). In contrast, a striking increase in moisture is found over southeastern North America (15% increase), where southerly atmospheric flow brings moist air to the continent. Outside these areas that are clearly affected by the cold North Atlantic Ocean, the responses of temperature and moisture are decoupled, with different causes for these temperature and hydroclimate responses. In the tropics, the hydroclimate response is governed by the southward shift of the intertropical convergence zone (ITCZ) due to the cooling of the North Atlantic Ocean. This causes drier conditions north of the equator and wetter conditions in the Southern Hemisphere tropics. The associated changes in soil moisture are relatively gradual here, taking up to two centuries to complete, suggesting that the impact of the ITCZ shift on the tropical hydroclimate is building up. Our experiment indicates that Southern Hemisphere continents experienced a small cooling (less than 0.5 °C) during the Younger Dryas, caused by the negative radiative forcing associated with reduced atmospheric methane concentrations and enhanced dust levels. In our simulation, the bi-polar seesaw mechanism is relatively weak, so that the associated warming of the South Atlantic Ocean is not overwhelming the reduction in radiative forcing. Our results thus indicate that in the tropics and/or Southern Hemisphere, the cooling is a response to the negative radiative forcing, while the hydroclimatic changes are predominantly resulting from ITCZ variations. Consequently, when interpreting hydroclimatic proxy records from these regions, data should not be compared directly to key records from high latitudes, such as Greenland ice core stable isotope records.

Description: We present high-resolution resistivity imaging of gas hydrate pipe-like structures, as derived from marine controlled-source electromagnetic (CSEM) inversions that combine towed and ocean-bottom electric field receiver data, acquired from the Nyegga region, offshore Norway. Two-dimensional CSEM inversions applied to the towed receiver data detected four new prominent vertical resistive features that are likely gas hydrate structures, located in proximity to a major gas hydrate pipe-like structure, known as the CNE03 pockmark. The resistivity model resulting from the CSEM data inversion resolved the CNE03 hydrate structure in high resolution, as inferred by comparison to seismically constrained inversions. Our results indicate that shallow gas hydrate vertical features can be delineated effectively by inverting both ocean-bottom and towed receiver CSEM data simultaneously. The approach applied here can be utilised to map and monitor seafloor mineralisation, freshwater reservoirs, CO2 sequestration sites and near-surface geothermal systems.

Description: Reef-associated predators are thought to influence the distribution of invertebrates of surrounding sediment habitats. In this study, we analyzed the predation pressure and the distributional patterns of macro- and meiofaunal assemblages in soft sediments surrounding two coral reef sites at Rapa Nui (Easter Island), in the central South Pacific. We tested the hypothesis that reef-associated predators negatively affect sediment-dwelling invertebrates, causing macro- and meiofauna to be less abundant and diverse in soft sediments near the reefs. As expected, predation intensity was greater nearer the reef than farther away, but macro- and meiofaunal assemblages did not differ significantly with distance from the reef. Taxon richness of macro- and meiofaunal assemblages were similar irrespective of distance from the reef. Only meiofauna showed significant variation in total abundance with distance from the reef, but this trend was not consistent between the two study sites. No gradient in sediment texture was observed with distance from the reef. Underwater video recordings at one study site also revealed that local hydrodynamics cause frequent disturbance and resuspension of the upper sediment layers. Our results suggest that soft-sediment assemblages are constantly reshuffled by oceanic waves, thereby blurring the potential effects of predation on invertebrate assemblages closer to the reef.

Description: Highlights • We have tested a hypothesis where a sill intrusion is present at depth near Lusi. • We have calculated the CO2 generation following the emplacing a 150 m sill in an organic rich sequence at 4.5 km. • This scenario may provide the CO2 currently emitted from Lusi, and are consistent with geological information. Abstract The Lusi mud eruption started in 2006 and is located near the Arjuno-Welirang volcanic complex in Northeastern Java. Lusi is characterized by the eruption of aqueous vapor, CO2, and CH4 in addition to mud breccia and boiling water. However, the ultimate driving force for the eruption remains unclear. Here we investigate if Lusi could have been driven by the heat released from a deep-seated igneous sill originating from the neighboring volcanic arc. We have used a 1D thermal model to calculate the production of CO2 from thermally matured organic matter in the contact aureole of a hypothetical 150 m thick sill. The sill is tentatively emplaced at 1100 °C at 4.5 km depth within the organic-rich Eocene Ngimbang Formation. The carbon gas produced from the thermal perturbation reaches a peak of 1357 kg/m2/y CO2 equivalents shortly after sill emplacement, stressing the efficiency of organic matter transformation in contact aureoles. Our simulations show that during the first 1000 years after emplacement, 53.5 ton CO2/m2 is produced in the contact aureole. When scaled to a sill size of 150 m × 25 km2, i.e., a sill volume of 3.75 km3, the aureole has the potential to generate a total of 1350 Mt CO2 during the first 1000 years, with a peak generation of about 34 Mt CO2/y. We conclude that contact metamorphism in our hypothetical geological scenario generates CO2 in the gigaton range and represents a plausible source for the Lusi gas.

Description: Highlights: • Picophytoplankton was able to grow under a wide range of environmental conditions varying seasonally. • Picophytoplankton variability was more marked in the bay than in the lagoon. • Picophytoplankton decreased in the lagoon under enhanced anthropogenic stress. • Synechococcus, the most abundant group of the community, was frequently present in both ecosystems during the monitoring. We investigated at the single cell level during 16months (June 2012 to September 2013) the temporal distribution of picophytoplankton (picoeukaryotes, Synechococcus and Prochlorococcus) communities in two contrasted ecosystems: the Bay of Bizerte characterised by an oligotrophic regime typical of the Mediterranean Sea and the Bizerte Lagoon that exhibits a mesotrophic/eutrophic state. We aimed at depicting seasonal variations and quantifying the relationships between the environmental factors and the structure and abundance of picophytoplankton communities. Results showed that picophytoplankton groups were able to grow under a wide range of environmental conditions varying seasonally, although their abundances and contributions to the total chlorophyll biomass significantly varied and showed importance in the Bay of Bizerte. Synechococcus was the most abundant group reaching 225*103 cells·cm-3 in the Bay and 278*103 cells·cm-3 in the lagoon. This group was present all over the year in both ecosystems. Structural equation model results pointed out a different configuration regarding the picophytoplankton environmental drivers. The complexity of the configuration, i.e. number of significant links within the system, decreased under enhanced eutrophication conditions. The less exposure to anthropogenic stress, i.e. in the Bay of Bizerte, highlight a larger role of nutrient and hydrological conditions on the seasonal variations of picophytoplankton, whereas a negative effect of eutrophication on picophytoplankton communities was unveiled in the Bizerte Lagoon. We stress that such influence may be exacerbated under expected scenarios of Mediterranean warming conditions and nutrient release in coastal ecosystems.

Description: In the present study, various MD methods including Equilibrium Molecular Dynamics (EMD) and two different Non- Equilibrium Molecular Dynamics (NEMD) methods are studied and compared for evaluation of the thermal conductivity and total and partial phonon density of states spectra in silicon nanowires, as a case study. The thermal conductivity of nanowires was determined as a function of length, cross section width and temperature. According to the results obtained via various MD methods, the thermal conductivity increases by increasing the length and cross section or decreasing the MD temperature. However, it was observed that despite the same initial conditions, different MD methods could predict considerably different values for the thermal conductivity which was found to be due to the different equilibrium temperature achieved in different methods. The total phonon density of states spectra was then employed to analyze the phonon transport properties of a 45 nm SiNW simulated using various MD techniques. Two major peaks were observed at around 16.5 and 5 THz which are attributed to Si-Si bond modes. Finally partial phonon density of states was calculated to differentiate the contributions to the phonon DOS from surface atoms compared to central atoms.

Description: Major and trace elements, and volatile components have been measured in melt inclusions in olivine from fresh 2.7 Ga old komatiites from the Reliance Formation of the Belingwe Greenstone Belt, Zimbabwe. Reconstructed compositions of melt inclusions contain 20–23.5 wt% MgO and up to 0.3 wt% H2O; these compositions probably represent those of the erupted lava. In inclusions in relatively evolved (low Fo) olivines, an excess of Na2O, CaO, Li, La, Cu, Rb, Y, Sc as well as volatile components (H2O, F, Cl and S) relative to other highly incompatible elements is attributed to assimilation of seawater altered mafic material. No assimilation signature is observed for the most primitive melt inclusions hosted in the magnesium rich olivines. The primary melt composition, estimated using melt inclusions in the most magnesian olivine (Fo 93.5), contains up to 27.5 wt% MgO and ca. 0.2 wt% H2O. The presence of H2O slightly depressed the liquidus temperature to ca. 1513 °C. Our results suggest formation of the Belingwe komatiite magma at ca. 7 GPa pressure and ca. 1790 °C temperature in a mantle plume. The plume picked up water and probably chlorine through interaction with a hydrous transition mantle zone in the way similar to that previously proposed by Sobolev et al. (2016) for komatiites in Canada.

Description: The influences of North Atlantic biases on multiyear predictability of unforced surface air temperature (SAT) variability are examined in the Kiel Climate Model (KCM). By employing a freshwater flux correction over the North Atlantic to the model, which strongly alleviates both North Atlantic sea surface salinity (SSS) and sea surface temperature (SST) biases, the freshwater flux-corrected integration depicts significantly enhanced multiyear SAT predictability in the North Atlantic sector in comparison to the uncorrected one. The enhanced SAT predictability in the corrected integration is due to a stronger and more variable Atlantic Meridional Overturning Circulation (AMOC) and its enhanced influence on North Atlantic SST. Results obtained from preindustrial control integrations of models participating in the Coupled Model Intercomparison Project Phase 5 (CMIP5) support the findings obtained from the KCM: models with large North Atlantic biases tend to have a weak AMOC influence on SAT and exhibit a smaller SAT predictability over the North Atlantic sector.

Description: Most of the anthropogenic radionuclide 129I released to the marine environment from the nuclear fuel reprocessing plants (NFRP) at Sellafield (England) and La Hague (France) is transported to the Arctic Ocean via the North Atlantic Current and the Norwegian Coastal Current. 129I concentrations in seawater provides a powerful and well-established radiotracer technique to provide information about the mechanisms which govern water mass transport in the Nordic Seas and the Arctic Ocean and is gaining importance when coupled with other tracers (e.g. CFC, 236U). In this work, 129I concentrations in surface and depth profiles from the Nordic Seas and the North Atlantic (NA) Ocean collected from four different cruises between 2011 and 2012 are presented. This work allowed us to i) update information on 129I concentrations in these areas, required for the accurate use of 129I as a tracer of water masses; and ii) investigate the formation of deep water currents in the eastern part of the Nordic Seas, by the analysis of 129I concentrations and temperature-salinity (T-S) diagrams from locations within the Greenland Sea Gyre. In the Nordic Seas, 129I concentrations in seawater are of the order of 109 at·kg− 1, one or two orders of magnitude higher than those measured at the NA Ocean, not so importantly affected by the releases from the NFRP. 129I concentrations of the order of 108 atoms·kg− 1 at the Ellet Line and the PAP suggest a direct contribution from the NFRP in the NA Ocean. An increase in the concentrations in the Nordic Seas between 2002 and 2012 has been detected, which agrees with the temporal evolution of the 129I liquid discharges from the NFRPs in years prior to this. Finally, 129I profile concentrations, 129I inventories and T-S diagrams suggest that deep water formation occurred in the easternmost area of the Nordic Seas during 2012.

Description: On continental margins, upward migration of fluids from various sources and various subsurface accumulations, through the sedimentary column to the seafloor, leads to the development of cold seeps where chemical compounds are discharged into the water column. MarsiteCruise was undertaken in November 2014 to investigate the dynamics of cold seeps characterized by vigorous gas emissions in the Sea of Marmara (SoM).A previous paper published by Bourry et al. (2009) presented the gas geochemistry of three seeps sampled along three different segments in the SoM. Their findings showed that the seeps were sourced by three different reservoirs. In this paper, seventeen seeps were investigated to determine the gas sources, unravel reservoir contributions, and estimate their level of mixing. The molecular and stable isotope compositions of the gas compounds were determined to establish the empirical diagrams that usually allow to delineate source domains. The results provide insights into the complexities of source mixing within the sedimentary column of the SoM before emission of the gases into the water column. The seep gases originate from deep thermogenic or microbial hydrocarbon sources, or from a CO2-rich source. Microbial sources producing methane from primary methanogenesis have been identified in the Tekirdağand the Çinarcik basins. In addition, six different thermogenic reservoirs or six different pathways of migration are responsible for the supply of gas to the seeps on the highs and in the western basin. Five of them are undergoing biodegradation followed by secondary methanogenesis, thereby providing additional sources of microbial methane to the seeps. Overall, the gases emitted by the seventeen seeps consist of variable mixtures of different components from two or three sources.

Description: This study reports for the first time boreal to subarctic intertidal foraminiferal assemblages from saltmarshes at Borgarnes and Faskrudsfjördur on Iceland. The composition of living and dead foraminiferal assemblages was investigated along transects from the tidal flat to the highest reach of halophytic plants. The foraminiferal assemblages from Borgarnes showed 18 species in the total foraminiferal assemblage of which only 7 species were recorded in the living fauna. The assemblages were dominated by agglutinated taxa, whereas 3 calcareous species were recorded, of which only Haynesina orbicularis was found in the living fauna. The distribution limit of calcifying species corresponds to the lower boundary of the lower saltmarsh vegetation zone. Furthermore, calcareous tests showed many features of dissolution, which is an indication of a carbonate corrosive environment. The species forming the dead assemblages were mainly derived from the ambient intertidal areas and were displaced by tidal currents into the saltmarsh. The foraminiferal assemblages from Faskrudsfjördur showed two species, of which only one species was recorded in the living fauna. The assemblage was dominated by the agglutinated foraminifer Trochaminita irregularis. The foraminiferal species recorded on Iceland were the same as commonly found elsewhere in Europa. Since no species was found which is endemic to North America, Iceland is considered part of the European bio province. The foraminiferal could have been immigrated to Iceland from Europe through warm water currents, migratory birds or marine traffic since the last Ice Age.

Description: Pseudovibrio is a marine bacterial genus members of which are predominantly isolated from sessile marine animals, and particularly sponges. It has been hypothesised that Pseudovibrio spp. form mutualistic relationships with their hosts. Here, we studied Pseudovibrio phylogeny and genetic adaptations that may play a role in host colonization by comparative genomics of 31 Pseudovibrio strains, including 25 sponge isolates. All genomes were highly similar in terms of encoded core metabolic pathways, albeit with substantial differences in overall gene content. Based on gene composition, Pseudovibrio spp. clustered by geographic region, indicating geographic speciation. Furthermore, the fact that isolates from the Mediterranean Sea clustered by sponge species suggested host-specific adaptation or colonization. Genome analyses suggest that Pseudovibrio hongkongensis UST20140214-015BT is only distantly related to other Pseudovibrio spp., thereby challenging its status as typical Pseudovibrio member. All Pseudovibrio genomes were found to encode numerous proteins with SEL1 and tetratricopeptide repeats, which have been suggested to play a role in host colonization. For evasion of the host immune system, Pseudovibrio spp. may depend on type III, IV and VI secretion systems that can inject effector molecules into eukaryotic cells. Furthermore, Pseudovibrio genomes carry on average seven secondary metabolite biosynthesis clusters, reinforcing the role of Pseudovibrio spp. as potential producers of novel bioactive compounds. Tropodithietic acid, bacteriocin and terpene biosynthesis clusters were highly conserved within the genus, suggesting an essential role in survival e.g. through growth inhibition of bacterial competitors. Taken together, these results support the hypothesis that Pseudovibrio spp. have mutualistic relations with sponges.

Description: A new cyclic hexapeptide, cyclo-(Gly-Leu-Val-IIe-Ala-Phe), named bacicyclin (1), was isolated from a marine Bacillus sp. strain associated with Mytilus edulis. The sequences of the amino acid building blocks of the cyclic peptide and its structure were determined by 1D- and 2D-NMR techniques. Marfey's analysis showed that the amino acid building blocks had L-configuration in all cases except for alanine and phenylalanine, which had D-configuration. Bacicyclin (1) exhibited antibacterial activity against the clinically relevant strains Enterococcus faecalis and Staphylococcus aureus with minimal inhibitory concentration values of 8 and 12 μM, respectively. These results demonstrate the potential of marine bacteria as a promising source for the discovery of new antibiotics.

Description: Although many of the regions on and close to the mid-ocean ridges have been extensively mapped and sampled, the abyssal intraplate regions remain essentially unsampled and unmapped, leaving huge gaps in our understanding of their geologic history and present activity. Prominent bathymetric features in these intraplate regions are fracture zones. Here we present bathymetric and sampling information from a transatlantic transect along the Vema Fracture Zone (ca. 11 °N), covering crustal ages from 109 − 0 Ma on the African plate and 0–62 Ma on the South American plate. The Vema Fracture Zone is the intraplate trace of the active Vema Transform plate boundary, which offsets the present-day Mid-Atlantic Ridge by ca. 300 km left-laterally, juxtaposing zero-age crust with crust of 20 million years age. Our results show clear evidence of tectonic activity along most of the Fracture Zone, in most places likely associated with active fluid flow. Within the active Vema Transform at crustal ages of ca. 10 Ma we found clear indications of fluid flow both in the sediments and the overlying water column. This region is 〉120 km from the nearest spreading axis and increases by almost an order of magnitude the maximum off-axis distance that active hydrothermal discharge has been found in the oceanic crust. Sampling of the igneous seafloor was possible at all crustal ages and the accretionary fabric imprinted on the plate during its production was prominent everywhere. Seafloor sediments show signs of extensive bioturbation. In one area, high concentrations of spherical Mn-nodules were also found and sampled. At the end of the transect we also mapped and sampled the Puerto Rico Trough, a 〉8000 m deep basin north of the Caribbean arc. Here the seafloor morphology is more complicated and strongly influenced by transpressive tectonics.

Description: Seamounts are amongst the most common physiographic structures of the deep-ocean landscape, but remoteness and geographic complexity have limited the systematic collection of integrated and multidisciplinary data in the past. Consequently, important aspects of seamount ecology and dynamics remain poorly studied. We present a data collection of ocean currents and raw acoustic backscatter from shipboard Acoustic Doppler Current Profiler (ADCP) measurements during six cruises between 2004 and 2015 in the tropical and subtropical Northeast Atlantic to narrow this gap. Measurements were conducted at seamount locations between the island of Madeira and the Portuguese mainland (Ampère, Seine Seamount), as well as east of the Cape Verde archipelago (Senghor Seamount). The dataset includes two-minute ensemble averaged continuous velocity and backscatter profiles, supplemented by spatially gridded maps for each velocity component, error velocity and local bathymetry.

Description: Complex multifault earthquake ruptures involving secondary faults emphasize the necessity to characterize their seismogenic potential better and study their relationship with major faults to improve the seismic hazard assessment of a region. High-resolution geophysical data were interpreted to make a detailed characterization of the Averroes Fault and the North Averroes Faults, which are poorly known secondary right-lateral strike-slip faults located in the central part of the Alboran Sea (western Mediterranean). These faults appear to have evolved since the Pliocene as part of a distributed dextral strike-slip shear zone in response to local strain engendered by the diverging movement of the Carboneras Fault to the north, and the Yusuf and Alboran Ridge faults to the south. In addition, the architecture of these faults suggests that the Averroes Fault may eventually link with the Yusuf fault, thus leading to a higher seismogenic potential. Therefore, these secondary faults represent a hitherto unrecognized seismogenic hazard since they could produce earthquakes up to moment magnitude (Mw) 7.6. Our results highlight the importance of the role played by secondary faults in a specific kinematic framework. Their reciprocal linkage and their mechanical relationship with the main faults could lead to future complex fault ruptures. This information could improve fault source and earthquake models used in seismic and tsunami hazard assessment in this and similar regions.

Description: Fecundity of marine fish species is highly variable, but trade-offs between fecundity and egg quality have rarely been observed at the individual level. We investigated spatial differences in reproductive investment of individual European sprat Sprattus sprattus (Linnaeus 1758) females by determining batch fecundity, condition indices (somatic condition index and gonadosomatic index) as well as oocyte dry weight, protein content, lipid content, spawning batch energy content, and fatty acid composition. Sampling was conducted in five different spawning areas within the Baltic Sea between March and May 2012. Spawning sprat from the Kiel Bight were in a better nutritional condition compared to sprat from the Arkona Basin, Bornholm Basin, and the Gdansk Deep. These females were also producing up to twice as many oocytes, and invested significantly more energy in reproduction, than their counterparts sampled in the eastern part of the Baltic Sea. Still, oocytes produced by Kiel Bight sprat contained significantly lower fractions of the essential fatty acids 20:4 (n-6) and 22:6 (n-3). A seasonal trade-off between oocyte weight/lipid content and fecundity was found for Baltic sprat, albeit the gross energy invested into spawning remained constant. Observed spatial and seasonal differences in sprat reproductive investment may be linked to hydrographic conditions and food availability and will impact the survival probability of yolk-sac and first feeding larvae. These findings indicate that Baltic sprat is able to adapt its reproductive tactics to the highly variable pelagic habitat of the Baltic Sea.

Description: Highlights • Large seafloor depressions with diameters of up 10 km across have been mapped on the southern Chatham Rise, New Zealand. • Seismic reflection data show scarce indications for vertical fluid flow but no clear link between fluid flow and depressions. • Methane gas or methane hydrates appear to be absent on the southern Chatham Rise. • Seismic evidence suggests that vertical fluid flow was likely fuelled by polygonal faulting and silica diagenesis • The depressions are the results of erosion and sediment drift deposition of bottom currents associated with the Subtropical Front. Abstract Several giant seafloor depressions were investigated on the Chatham Rise offshore New Zealand using mainly bathymetric and seismic data, supplemented by sediment cores and reported porewater geochemistry data. The depressions have diameters of up to 11 km and occur on the southern flank of the Chatham Rise in water depths between 600 and 900 m, i.e. roughly underneath the location of the strongest thermal gradients of the Subtropical Front (STF) and characterized by eastward flowing currents. With up to 150 m of relief the depressions cut into post-Miocene deposits. Some of the depressions are partially filled with drift deposits that have similar seismic characteristics as the surrounding sediments and consist of alternations of silty muds and silts. Seismic profiles also show completely filled depressions that no longer have a bathymetric expression. Despite several pipe structures indicating vertical fluid flow, neither active fluid seepage nor indications for past fluid seepage are present at the seafloor of the Chatham Rise. Also, both pore water geochemistry and geophysical data do not show indications for an existing or past gas hydrate system in the area. Instead, seismic data suggest widespread polygonal faulting and the presence of silica diagenetic fronts. The release of mineral-bound water during silica diagenesis or fluid expulsion during sediment compaction can explain the presence of vertical fluid flow features but not the giant depressions themselves. Instead, the depressions are interpreted as the result of scouring by strong bottom currents for which fluid venting may have created the nucleation points.

Description: Highlights • We document marine forearc deformation in the Northern Chile seismic gap. • Upper-plate normal faulting off Northern Chile locally extends close to the trench. • Normal faults indicate that past earthquakes may reached the shallow plate-boundary. Abstract Seismic rupture of the shallow plate-boundary can result in large tsunamis with tragic socio-economic consequences, as exemplified by the 2011 Tohoku-Oki earthquake. To better understand the processes involved in shallow earthquake rupture in seismic gaps (where megathrust earthquakes are expected), and investigate the tsunami hazard, it is important to assess whether the region experienced shallow earthquake rupture in the past. However, there are currently no established methods to elucidate whether a margin segment has repeatedly experienced shallow earthquake rupture, with the exception of mechanical studies on subducted fault-rocks. Here we combine new swath bathymetric data, unpublished seismic reflection images, and inter-seismic seismicity to evaluate if the pattern of permanent deformation in the marine forearc of the Northern Chile seismic gap allows inferences on past earthquake behavior. While the tectonic configuration of the middle and upper slope remains similar over hundreds of kilometers along the North Chilean margin, we document permanent extensional deformation of the lower slope localized to the region 20.8°S–22°S. Critical taper analyses, the comparison of permanent deformation to inter-seismic seismicity and plate-coupling models, as well as recent observations from other subduction-zones, including the area that ruptured during the 2011 Tohoku-Oki earthquake, suggest that the normal faults at the lower slope may have resulted from shallow, possibly near-trench breaking earthquake ruptures in the past. In the adjacent margin segments, the 1995 Antofagasta, 2007 Tocopilla, and 2014 Iquique earthquakes were limited to the middle and upper-slope and the terrestrial forearc, and so are upper-plate normal faults. Our findings suggest a seismo-tectonic segmentation of the North Chilean margin that seems to be stable over multiple earthquake cycles. If our interpretations are correct, they indicate a high tsunami hazard posed by the yet un-ruptured southern segment of the seismic gap.

Description: The aggregation behaviour of polymer-coated silver nanoparticles (AgNPs) was characterized in NaCl solutions, and in two seawaters of different salinities and dissolved organic matter (DOM) contents. Representative organic coatings i.e. tannic acid (TA), alginic acid (ALG), two gum Arabic samples (GAL and GAH), branched polyethylenimine (BPEI), and non-ionic surfactants (reference material NM-300K) were selected to cover a wide range of zeta-potentials. The stability in NaCl solutions, as determined from the rate of variation in hydrodynamic size within a timeframe of one hour, followed the order BPEI ≫ NM-300K ≈ GAL ≫ ALG ≈ TA ≫ GAH. In the seawater samples the order was NM-300K ≈ GAL ≫ ALG 〉 GAH 〉 TA ≈ BPEI, and only TA, GAL and NM-300K batches behaved as expected from the NaCl experiments. Remarkably, the BPEI sample showed the largest aggregation rate in the seawater sample with the highest DOM concentration (277 μM C). The GAH sample displayed a non-monotonic variation in aggregation rate with NaCl concentration, apparently due to concomitant precipitation of AgCl. The results indicate that non-electrostatic stabilization mechanisms and DOM-coating interactions are important for the prediction of stability and persistence of polymer-coated AgNPs in seawater.

Description: Highlights • Hindon Maar Complex is a new mid-Miocene Fossil-Lagerstätte in New Zealand. • Anoxia in maar lakes allowed exquisite preservation of plant and animal fossils. • The biota is from a lake and Nothofagus/podocarp/mixed broadleaf forest ecosystem. • Fossils record high diversity at humid, warm Southern Hemisphere mid-latitudes. Abstract This paper highlights the geology, biodiversity and palaeoecology of the Hindon Maar Complex, the second Miocene Konservat-Lagerstätte to be described from New Zealand. The Lagerstätte comprises four partly eroded maar-diatreme volcanoes, with three craters filled by biogenic and highly fossiliferous lacustrine sediments. The exceptionally well-preserved and diverse biota from the site is derived from a mid-latitude Southern Hemisphere lake-forest palaeoecosystem, including many fossil taxa not previously reported from the Southern Hemisphere. The most common macrofossils are leaves of Nothofagus, but the flora also includes conifers, cycads, monocots (such as Ripogonum and palms), together with Lauraceae, Myrtaceae and Araliaceae leaves and flowers. The small maar lakes were surrounded by Nothofagus/podocarp/mixed broadleaf forest growing under humid, warm temperate to subtropical conditions. The fossil fauna comprises insects in the orders Odonata, Hemiptera, Thysanoptera, Coleoptera, Diptera, Hymenoptera and Trichoptera, and the fish assemblage includes a non-migratory species of the Southern Hemisphere Galaxias (Galaxiidae) and a significant new record of the freshwater eel Anguilla (Anguillidae). The fossil assemblage also includes the first pre-Quaternary bird feathers from New Zealand and abundant coprolites derived from fish and volant birds, presumably waterfowl. Palynomorph analysis and a 40Ar/39Ar age of 14.6 Ma obtained from basanite associated with the maar complex indicate that the Hindon Maar Complex is of mid-Miocene age (Langhian; New Zealand local stage: Lillburnian). It thus provides a new and unique perspective on Neogene terrestrial biodiversity and biogeography in the Australasian region, around the end of the mid-Miocene thermal optimum and prior to late Miocene–Pleistocene climate cooling episodes when many warm-temperate and subtropical forest components became extinct in New Zealand.

Description: MarsiteCruise was undertaken in October/November 2014 in the Sea of Marmara to gain detailed insight into the fate of fluids migrating within the sedimentary column and partially released into the water column. The overall objective of the project was to achieve a more global understanding of cold-seep dynamics in the context of a major active strike-slip fault. Five remotely operated vehicle (ROV) dives were performed at selected areas along the North Anatolian Fault and inherited faults. To efficiently detect, select and sample the gas seeps, we applied an original procedure. It combines sequentially (1) the acquisition of ship-borne multibeam acoustic data from the water column prior to each dive to detect gas emission sites and to design the tracks of the ROV dives, (2) in situ and real-time Raman spectroscopy analysis of the gas stream, and (3) onboard determination of molecular and isotopic compositions of the collected gas bubbles. The in situ Raman spectroscopy was used as a decision-making tool to evaluate the need for continuing with the sampling of gases from the discovered seep, or to move to another one. Push cores were gathered to study buried carbonates and pore waters at the surficial sediment, while CTD-Rosette allowed collecting samples to measure dissolved-methane concentration within the water column followed by a comparison with measurements from samples collected with the submersible Nautile during the Marnaut cruise in 2007. Overall, the visited sites were characterized by a wide diversity of seeps. CO2- and oil-rich seeps were found at the westernmost part of the sea in the Tekirdag Basin, while amphipods, anemones and coral populated the sites visited at the easternmost part in the Cinarcik Basin. Methane-derived authigenic carbonates and bacterial mats were widespread on the seafloor at all sites with variable size and distributions. The measured methane concentrations in the water column were up to 377 μmol, and the dissolved pore-water profiles indicated the occurrence of sulfate depleting processes accompanied with carbonate precipitation. The pore-water profiles display evidence of biogeochemical transformations leading to the fast depletion of seawater sulfate within the first 25-cm depth of the sediment. These results show that the North Anatolian Fault and inherited faults are important migration paths for fluids for which a significant part is discharged into the water column, contributing to the increase of methane concentration at the bottom seawater and favoring the development of specific ecosystems.

Description: Deep stratigraphic constrains below the Indonesian Lusi mud eruption are currently lacking due to the absence of deep wells and good quality seismic data. A collection of carbonate clasts has been sampled from the Lusi site, active since its birth in 2006. These specimens are part of a large variety of lithotypes erupted from the main crater. The carbonates analysed comprise scleractinian coral and bivalve shell fragments, probably shallow-water in origin, and clasts consisting of planktonic foraminifera-bearing mudstone, from pelagic deposits. Selected rocks were analysed using planktonic foraminifera and 87Sr/86Sr dating with the aim to constrain their age and to improve the understanding of the, so far unknown, sequence of limestone deposits inferred at this site. Based on biostratigraphy using planktonic foraminifera, one group of samples reveal to belong to the Planktonic Foraminifera Zone M5b, with an age comprised between 16.29 and 15.10 Ma (Miocene, Latest Burdighalian to Langhian). The Sr isotope-based ages of clasts analysed for 87Sr/86Sr cover a larger time window spanning from Pliocene (Zanclean and Piacenzian), Miocene (Messinian) down to Eocene (Priabonian). The Pliocene and Messinian ages are unreasonably young from what is known of the local geology and one sample provided an 87Sr/86Sr age that is ∼8 My younger compared to that obtained from the planktonic foraminifera assemblage occurring in that sample. The discrepancy suggests that this and the other unduly young samples have possibly been contaminated by geological materials with radiogenic Sr isotope composition. The minimum age of 37.18 Ma obtained by 87Sr/86Sr from a well-preserved oyster shell indicates that some of the clasts can be attributed to the deep seated Ngimbang Formation. The dating has been combined with the interpretation of seismic profiles to investigate the stratigraphy of inferred carbonates below Lusi and the PRG-1 well located some kilometres to the north east. PRG-1 borehole data are also integrated and discussed in this study. The obtained results reveal multiple implications: a) the Tuban and Kujung Formations are overlapping at Lusi site; b) the Lusi feeder conduit brecciated and mobilized to the surface carbonate lithologies buried as deep as possibly ∼3.8 km as well as even older and deeper seated lithotypes from the Ngimbang Formation; c) since the deeper carbonate samples erupted in 2006 belong to the typically not overpressured Kujung Formation, an additional overpressure generated from deeper units (Ngimbang Formation) would be required to force these fragments to surface.

Description: For the first time we present the geochemical characterization of fluids emitted from the Arjuno-Welirang volcanic complex and compare the results with those obtained sampling the neighboring spectacular Lusi eruption site (Java Island, Indonesia). The isotopic composition of the hydrothermal and cold waters from the Welirang volcanic complex indicate a meteoric origin for these springs, with values ranging from −65 to −50 and −6 to −1‰ vs V-SMOW respectively for δD and δ18O. The water erupted from the Lusi site showed clustered higher δD and δ18O isotopic values, ranging around −6 and +10 vs V-SMOW respectively. We ascribe these results to mixing between hydrothermal fluids, meteoric water, saline formation fluids, and water released during clay mineral illitization ultimately altered by additional evaporation processes. The chemical and isotopic composition of fluids emitted from fumaroles and hydrothermal springs of the Welirang volcano showed a clear magmatic signature where a CO2-dominated gas reveals δ13CCO2 ranging between −5.9 and −2.4 and helium isotope with R/Ra = 7.3. These values are very close to those measured at Lusi site (R/Ra = 7) that also have high CO2/CH4 ratio (1.7–2.2) supporting the high contribution of magmatic gases. Moreover, a great contribution of andesitic water has been recognized in the water vapour emitted from the summit fumaroles. Converging geochemical data indicate that the plumbing system of the Lusi eruption site is connected at depth with the Arjuno-Welirang volcanic complex. These data support a scenario where hydrothermal fluids from the volcanic system migrated in the sedimentary basin triggering metamorphic reactions in the organic-rich sediments that ultimately resulted in a venting system at the surface. After eleven years of incessant activity this venting system remains constantly fed by the fluids from the volcanic complex and became world known as “Lusi”, the largest ongoing clastic geysering system on Earth.

Description: An extensive survey was carried out in the Sidoarjo district (East Java, Indonesia) to investigate the gas leaking properties along fractured zones coinciding with a strike-slip system, the Watukosek Fault System (WFS) in NE Java. This structure has been the focus of attention since the beginning of the spectacular Lusi mud eruption on the 29th May 2006. The sinistral strike-slip WFS originates from the Arjuno-Welirang volcanic complex, intersects the active Lusi eruption site displaying a system of antithetic faults, and extends towards the NE of Java where mud volcanic structures reside. In the Lusi region we completed a geochemical survey along three profiles combining measurements of a) 220Rn and 222Rn activity, b) CO2 and CH4 soil gas content, c) CO2 and CH4 fluxes, and d) gas analyses. The profiles are up to ∼1.2 km long and intersect perpendicularly areas with intense fracturing and surface deformation along the WFS. The purpose was to investigate the presence and origin of soil degassing activity in potentially active fault zones. Results show that the peripheral sectors of the profiles have high 220Rn activity and reduced CO2 and CH4 fluxes and concentrations. This suggests low fluids migration that could be affected by shallow circulation. In contrast, the segments of profiles intersecting the fractured zones have the highest 222Rn activity, CO2 and CH4 flux and gas concentration values. The relationship existing among the measured parameters suggests that the WFS acts as a preferential pathway for active rise of deep fluids. The presence of such advective processes is suggested by the relatively high rate of migration needed to obtain anomalies of short-lived 222Rn in the soil pores. Gas molecular and isotopic composition reveals that all sampled localities have a mixed hydrocarbon origin implying the presence of shallow microbial and deeper thermogenic hydrocarbons. CO2 isotopic values (δ13C-CO2 ranges between −9.48‰ and 4.12‰ V-PDB) indicate the presence of mantle derived CO2 and thermo-metamorphic CO2 suggesting that elevated temperatures have a key role in this active system. The samples collected from fractured and faulted zones reveal to have gas composition similar to that obtained from Lusi crater, indicating deep origin fluids.

Description: Detailed analysis of two-dimensional seismic lines acquired in the NE Java basin has been performed to unravel the subsurface geology of the region around the Lusi mud eruption. This work revealed the existence of a system consisting of a complex set of faults, here called the Watukosek fault system, forming triangular deformation zones converging at the top of the early Miocene Carbonates. This system continues downwards with vertical individual fault segments, often bordering the steep margins of the carbonate platforms. The analysis of data includes the interpretation of seismic lines, regional structural data inferred from basement gravity maps and present-day main direction of stress. Results suggest that a possible rotation of stress direction from N40E-S40W to N-S occurred during the post-Miocene history of the Java back arc tectonic evolution. The Watukosek fault system was first generated as a tensional lineament during the E-W sinistral transpressive strain, which involved the basement. In this phase, synthetic and antithetic Riedel faults formed, the former controlling the NW-SE orientation of the structural highs represented by Oligo-Miocene carbonate platforms. As a consequence of the rotation of the main principal stress direction to a N-S direction, the Watukosek fault and similar parallel lineaments became sinistral Riedel shears, developing intense triangular deformation zones. Based on the stratigraphic position of gentle anticlinal deformations with axis corresponding to the N40E-S40W oriented triangular deformation zones, the transpressive strain linked to the N-S main stress compression occurred likely in the Late Pliocene-Early Pleistocene. The detailed examination of a) stratigraphy at the wells BJP-1 and Porong-1 as well as b) the seismo-stratigraphic architecture of the entire succession in the study area, allowed a new subsurface interpretation and revision of the stratigraphic units below Lusi. The thick Early Miocene Tuban Formation is found sandwiched between the coheval Upper Kalibeng Formation and the Early Miocene Carbonate of the Kujung Formation, which in turn overlain the older Ngimbang Formation.

Description: From its inception May 29, 2006, the Lusi mud eruption has continuously erupted fluids, boiling mud and clasts through large active vents approximately 100 m in diameter. In 2016, we conducted a Dynamic Gravity survey (DG) using a network built over four locations and two Continuous Gravity-monitoring (CG) experiments to monitor the eruption activity. The CG was done for 8 days from June 2nd to June 10th, and for 9 days from August 20th to August 29th, 2016. Atmospheric pressure and atmospheric temperature variations were recorded during each experiment to constrain potential environmental effects, and Earth and oceanic tide effects were removed from gravity signals (CG and DG). Atmospheric pressure effects were removed from CG gravity signals. At the station, closest to the hydrothermal pond, the DG survey results show a gravity increase of ∼0.009 mGal month−1, which we interpret as the growth of the mud edifice in the central area. CG-monitoring shows that gravity variations occur at a period of 12–13 h, with amplitudes reaching up to 0.020 mGal. We interpret this as relating to density variation of the rising mud mixture (fluids + coarse mud + clasts). The observed 12–13 h period variations appear to indicate that tides may have some control on the density change of rising mud mixture by triggering the release of gas trapped at depth. Our 3D gravity results around the Lusi vents show that density variations range from 100 kg m−3 to 775 kg m−3. Similarly, vent diameters better constrain density contrasts occurring within the caldera zone, which is more likely to range between 400 and 450 kg m−3, and is equivalent respectively to 27% and 31% of gas ratio change over time.

Description: Since May 2006, the active Lusi mud eruption has continuously erupted boiling mud. During the early stages of the eruption in 2006, previous gravity studies showed that the piling up of the mud constructed a large edifice that subsides within the unconsolidated sediment especially around a large area over the active crater. After ten years of continuous eruption, the size of the edifice has grown significantly over a surface framed by 10 m tall containment embankments. In 2015 and 2016, several land gravity surveys were carried out to investigate the structure of the mud edifice and the effect of local geological active features. The new residual Bouguer anomaly map, calculated for a reference density of 2670 kg m−3, shows significant changes in the local gravity field in comparison to the previous 2006-gravity map survey. The new data set shows that the gravity decrease is generally restricted along the faulted and fractured zones, around erupting vents, and in the southern part of the mud edifice. Maximum gravity variations reach 1.4 mGal in some areas of the mud edifice. In the region outside the embankment, the gravity reductions are 0.6 mGal E-W and 1.0 mGal N-S. A second vertical derivative analysis of gravity data indicates that the mud edifice continues to pile up and subside mostly in the western and southern part of the edifice. Results of a 3D forward model of a vertical cylinder shape allowed characterising the extent of compacted material along the Watukosek fault system that originates from the neighbouring volcanic arc and crosses the mud edifice. Our results support the hypothesis of local pinched volume of mud ongoing between the subsided and uplifted masses of mud. The density of compacted mud breccia material increases between 16 and 27%. Gravity data also shows that the Lusi mud edifice is built over an extended NW-SE gravity increase, interpreted as a sediment density variation within the basin, and which is parallel to the trend of the Basin.

Description: The Lusi mud eruption in East Java, Indonesia, has continued unabated for more than ten years with no end in sight. This review summarizes what has been learned about this intriguing system, from its inception to the current (2017) well-established metastable geysering system that continuously erupts mud breccia, gas, steam, and water. We discuss the initiation of Lusi, highlighting discrepancies and evidence through the published data, to build a comprehensive database that emphasizes how the results converge towards a natural scenario of this system. We argue that attempts to understand, constrain, or predict the behaviour of this system that rely on a drilling trigger can not explain subsequent observations. On the other hand, we show that a well-constrained conceptual model recognizing Lusi as a volcanically-linked hydrothermal system, has provided important insights for the documented observations over the last eleven years. The response of Lusi to the Yogykarta earthquake falls directly within the range of earthquake triggering phenomena (globally) of similar hydrothermal/geothermal systems, suggesting a natural trigger as the more likely culprit for the Lusi phenomenon. We also offer some future directions of additional scrutiny for understanding this newborn, tectonic scale, volcanic-hydrothermal complex.

Description: The spectacular Lusi mud eruption started in northeast Java, Indonesia, the 29th of May 2006 following a M6.3 strike slip seismic event. After the earthquake several mud pools aligned along a NE-SW direction appeared in the Sidoarjo district. The most prominent eruption site was named Lusi. Lusi is located ∼10 km to the NE of the northernmost cone of the Arjuno-Welirang volcanic complex with which it is connected by the Watukosek Fault System. In this study, we applied the HVSR method, which is a common tool used for site effect investigations as well as to infer buried structures and reconstruct sub-surface geology. The method is based on the ratio of the horizontal to vertical components of ground motion and it generally exhibits a peak corresponding to the fundamental frequency of the site. Spectral ratio results highlight a fundamental frequency band between 0.4 and 1.0 Hz in the Lusi neighborhood. We interpret these peaks as related to the velocity lithological contrast at depth between alluvial deposits and bluish grey clay. Our analysis also highlights the presence of a “depocenter”, characterized by fundamental frequency up to 0.3 Hz, which is interpreted as the subsidence caused by withdrawal of mud and fluids from depth (as also shown by the comparison of the HVSR results with gravimetry data). Moreover, in the area of the Lusi vent a broad-band frequency range is related to the Lusi conduit. In this paper, we show that detailed microtremor surveys could be used as a preliminary and fast approach to locate mud conduits with sufficient precision.

Description: The Lusi mud eruption in East Java, Indonesia, is an active clastic-dominated geyser and a sedimentary hosted hydrothermal system that has generated wide interest across many disciplines. This moderate-to-high enthalpy system is driven by multiphase and multicomponent processes, fluid and rock mechanics, and heat transport processes, all which present challenges in developing realistic numerical models of the underlying physics. We develop a hydrogeological conceptual model for this deep and complex hydrothermal system, and construct an appropriate 3D geological model using the available data. This geological model then serves as the basis for numerically simulations that include some of the dominant processes driving Lusi. We adopt a flexible continuum approach with an efficient numerical simulator based on the 3D geological model representing the deep structures of this hydrothermal system and geothermal reservoir by incorporating borehole information and seismic data obtained in the framework of the LUSI Lab project. The geological model is transformed into a computational grid using binary space partitioning (BSP) of the input geometry and octree refinement on the grid to perform multi-physics simulations. Thermodynamic calculations using the equation of state for a CO2-bearing aqueous NaCl solution suggest that Lusi is a two-phase flow system (Water/CO2). Finally, we present initial results from a simple hydro-mechanical multiphase numerical model that simulates processes that likely contributed to the initiation of Lusi.

Description: Clastic eruptions are the surface expression of piercement structures such as mud volcanoes or hydrothermal vent complexes and involve subsurface sediment remobilisation and fluid flow processes. During these eruptions, many different processes are involved over a wide range of temporal and physical scales, which makes it a highly challenging multi-phase and multi-processes system to model. Field studies on piercement structures rarely include monitoring and detailed descriptions of clastic eruptions, and only a few attempts have been made to model fluid flow during these events. Moreover, these models have usually only considered one or two dimensions and/or have a limited spatial resolution. In this paper, we summarise the elements that are relevant for modelling fluid flow during clastic eruption: the geometry of the system, the ascending material and the host rocks. We present the main challenges associated with the identification of processes and quantification of parameters. By analogy to magmatic systems, we suggest that the type of clastic eruptions could be controlled by the liquid-gas flow pattern in the conduit. Effusive eruptions could be explained in terms of annular flows, while slug or churn flows could be expected during explosive events. We also propose that the viscosity of liquid mud controls the presence of slug flows in the conduit. We then review the two main approaches that have been proposed to model the flow dynamics in the active conduits, Darcy and Navier-Stokes, as well as their key parameters and their validity. Finally, we discuss the limits of the previously employed models and suggest further work directions to improve our understanding of clastic eruptions.

Description: We study the local seismicity in East Java around the Arjuno-Welirang volcanic complex that is connected via the Watukosek Fault System, to the spectacular Lusi eruption site. Lusi is a sediment-hosted hydrothermal system which has been erupting since 2006. It is fed by both mantellic and hydrothermal fluids, rising and mixing with the thermogenic gases and other fluids from shallower sedimentary formations. During a period of 24 months, we observe 156 micro-seismic earthquakes with local magnitudes ranging from ML0.5 to ML1.9, within our network. The events predominantly nucleate at depths of 8–13 km below the Arjuno-Welirang volcanic complex. Despite the geological evidence of active tectonic deformation and faulting observed at the surface, practically no seismicity is observed in the sedimentary basin hosting Lusi. Although we cannot entirely rule out artifacts due to an increased detection threshold in the sedimentary basin, the deficit in significant seismicity suggests aseismic deformation beneath Lusi due to the large amount of fluids that may lubricate the fault system. An analysis of focal mechanisms of nine selected events around the Arjuno-Welirang volcanic complex indicates predominantly strike-slip faulting activity in the region SW of Lusi. This type of activity is consistent with observable features such as fault escarpment, river deviation and railroad deformation; suggesting that the Watukosek fault system extends from the volcanic complex towards the NE of Java. Our results point out that the tectonic deformation of the region is characterized by a segmented fault system being part of a broader damage zone, rather than localized along a distinct fault plane.