ALBERT

Description: Eutrophication-driven harmful algal blooms (HABs) can have secondary effects on larval fishes that rely on estuaries as nurseries. However, few studies worldwide have quantified these effects despite the global rise in eutrophication. This study presents a novel approach using biochemical body condition analyses to evaluate the impact of HABs on the growth and body condition of the larvae of an estuarine resident fish. Recurrent phytoplankton blooms of Heterosigma akashiwo occur in the warm-temperate Sundays Estuary on the southeast coast of South Africa. The response in body condition and assemblage structure on larval estuarine roundherring (Gilchristella aestuaria) was measured in conjunction with bloom conditions, water quality and zooplanktonic prey and predators. Larvae and early juveniles were sampled during varying intensity levels, duration and frequency of hypereutrophic blooms. This study demonstrated that extensive HABs could significantly impact larval roundherring, G. aestuaria, by decreasing larval nutritional condition and limiting their growth, resulting in poor grow-out into the juvenile phase. Poor condition and growth may likely affect recruitment success to adult populations, and since G. aestuaria is an important forage fish and zooplanktivore, poor recruitment will hold consequences for estuarine food webs.

Description: The Observing Air–Sea Interactions Strategy (OASIS) is a new United Nations Decade of Ocean Science for Sustainable Development programme working to develop a practical, integrated approach for observing air–sea interactions globally for improved Earth system (including ecosystem) forecasts, CO2 uptake assessments called for by the Paris Agreement, and invaluable surface ocean information for decision makers. Our “Theory of Change” relies upon leveraged multi-disciplinary activities, partnerships, and capacity strengthening. Recommendations from 〉40 OceanObs’19 community papers and a series of workshops have been consolidated into three interlinked Grand Ideas for creating #1: a globally distributed network of mobile air–sea observing platforms built around an expanded array of long-term time-series stations; #2: a satellite network, with high spatial and temporal resolution, optimized for measuring air–sea fluxes; and #3: improved representation of air–sea coupling in a hierarchy of Earth system models. OASIS activities are organized across five Theme Teams: (1) Observing Network Design & Model Improvement; (2) Partnership & Capacity Strengthening; (3) UN Decade OASIS Actions; (4) Best Practices & Interoperability Experiments; and (5) Findable–Accessible–Interoperable–Reusable (FAIR) models, data, and OASIS products. Stakeholders, including researchers, are actively recruited to participate in Theme Teams to help promote a predicted, safe, clean, healthy, resilient, and productive ocean.

Description: Strong anisotropy of seismic velocity in the Earth’scrust poses serious challenges for seismic imaging. Where in situ seismic properties are not available the anisotropy can be determined from velocity analysis of surface and borehole seismic profiles. This is well established for dense, long-offset reflection seismic data. However, it is unknown how applicable this approach is for sparse seismic reflection data with low fold and short offsets in anisotropic metamorphic rocks. Here we show that anisotropy parameters can be determined from a sparse 3D data set at the COSC-1 borehole site in the Swedish Caledonides and that the results agree well with the seismic anisotropy parameters determined from seismic laboratory measurements on core samples. Applying these anisotropy parameters during 3D seismic imaging improves the seismic image of the high amplitude reflections especially in the vicinity of the lower part of the borehole. Strong reflections in the resulting seismic data show good correlation with the borehole-derived lithology. Our results aid the interpretation and extrapolation of the seismic stratigraphy of the Lower Seve Nappe in Jämtland and other parts in the Caledonides.

Description: The southern boundary of the Cayman Trough in the Caribbean is marked by the Swan Islands transform fault (SITF), which also represents the ocean-continent transition of the Honduras continental margin. This is one of the few places globally where a transform continental margin is currently active. The CAYSEIS experiment acquired an ∼165-km-long seismic refraction and gravity profile (P01) running across this transform margin, and along the ridge-axis of the Mid-Cayman Spreading Centre (MCSC) to the north. This profile reveals not only the crustal structure of an actively evolving transform continental margin, that juxtaposes Mesozoic-age continental crust to the south against zero-age ultraslow spread oceanic crust to the north, but also the nature of the crust and uppermost mantle beneath the ridge-transform intersection (RTI). The traveltimes of arrivals recorded by ocean-bottom seismographs (OBSs) deployed along-profile have been inverse and forward modelled, in combination with gravity modelling, to reveal an ∼25-km-thick continental crust that has been continuously thinned over a distance of ∼65 km to ∼10 km adjacent to the SITF, where it is juxtaposed against ∼3-4-km-thick oceanic crust. This thinning is primarily accommodated within the lower crust. Since Moho reflections are only sparsely observed, and, even then, only by a few OBSs located on the continental margin, the 7.5 km s-1 velocity contour is used as a proxy to locate the crust-mantle boundary along-profile. Along the MCSC, the crust-mantle boundary appears to be a transition zone, at least at the seismic wavelengths used for CAYSEIS data acquisition. Although the traveltime inversion only directly constrains the upper crust at the SITF, gravity modelling suggests that it is underlain by a higher density (〉3000 kg m-3) region spanning the width (∼15 km) of its bathymetric expression, that may reflect a broad region of metasomatism, mantle hydration or melt-depleted lithospheric mantle. At the MCSC ridge-axis to the north, the oceanic crust appears to be forming in zones, where each zone is defined by the volume of its magma supply. The ridge tip adjacent to the SITF is currently in a magma rich phase of accretion. However, there is no evidence for melt leakage into the transform zone. The width and crustal structure of the SITF suggests its motion is currently predominantly orthogonal to spreading. Comparison to CAYSEIS Profile P04, located to the west and running across-margin and through 10 Ma MCSC oceanic crust, suggests that, at about this time, motion along the SITF had a left-lateral transtensional component, that accounts for its apparently broad seabed appearance westwards.

Description: Diatoms often dominate phytoplankton in temperate, polar and upwelling regions. Decreases in silicate availability or silicon to nitrogen (Si:N) ratios may induce silicon limitation in diatoms and lower their proportion within phytoplankton communities. The effects of such changes on the nutritional quality of phytoplankton are not well understood. To examine how changing Si:N ratios affect plankton nutritional value, we applied a range of Si:N ratios on a natural plankton community and manipulated grazing pressure to assess top-down effects of copepod selective grazing. Diatom proportion in phytoplankton increased with increasing Si:N ratios and so did phytoplankton nutritional quality in terms of major fatty acid concentrations, such as polyunsaturated fatty acids, docosahexaenoic (DHA) and eicosapentaenoic (EPA) acids. However, stoichiometric quality (carbon to nitrogen and carbon to phosphorus ratios), DHA:EPA and omega 3:6 (omega 3:omega 6) ratios declined with increasing Si:N ratios, suggesting that proportions between essential compounds in copepod diet may be more favorable in lowered Si:N ratios. Copepods had a negative effect on DHA contents, DHA:EPA and omega 3:omega 6 ratios, indicating possible selective grazing on more nutritious plankton. Our findings show that declining silicate concentrations can affect stoichiometric and biochemical quality of phytoplankton, which copepods can also moderate by selective grazing.

Description: Temperature and dehydration stress are two major co-occurring environmental stressors threatening the physiology, biochemistry, and ecology of insects. As such, understanding adaptive responses to desiccation stress is critical for predicting climate change impacts, particularly its influence on insect invasions. Here, we assessed water balance and desiccation resistance of the invasive Tuta absoluta (Meyrick, 1917) (Lepidoptera: Gelechiidae), and infer how eco-physiology shapes its niche. We measured basal body water and lipid content, water loss rates (WLRs), and desiccation resistance in larvae (second to fourth instars) and adults. Body -water, -lipid, and WLRs significantly varied across life stages. Second instars recorded the lowest while fourth instars exhibited the highest body water and lipid content. Adult body water and lipid content were higher than second and third instars and lower than fourth instars while proportion of body water and lipid contents were highest in adults and second larval instars respectively. Water loss rates were significantly highest in fourth-instar larvae compared to other life stages, but differences among stages were less apparent at longer exposure durations (48 h). Desiccation resistance assays showed that second instars had greatest mortality while fourth-instar larvae and adults were the most desiccation tolerant. Our results show that T. absoluta fourth-instar larvae and adults are the most resilient developmental stages and potentially contribute most to the invasion success of the pest in arid environments. Incorporation of these species-specific eco-physiological traits in predictive models can help refine invasive species potential spread under changing climates.

Description: The fishery for Northern Atlantic cod (Gadus morhua) off Newfoundland and Labrador, Eastern Canada, presents the most spectacular case of an exploited stock crashed in a few decades by an industrial bottom trawl fishery under a seemingly sophisticated management regime after half a millennium of sustainable fishing. The fishery, which had generated annual catches of 100000 to 200000 tonnes from the beginning of the 16th century to the 1950s, peaked in 1968 at 810000 tonnes, followed by a devastating collapse and closure 24 years later. Since then, stock recovery may have been hindered by premature openings, with vessels targeting the remains of the cod population. Previous research paid little attention towards using multicentury time series to inform sustainable catches and recovery plans. Here, we show that a simple stock assessment model can be used to model the cod population trajectory for the entire period from 1508 to 2019 for which catch estimates are available. The model suggests that if fishing effort and mortality had been stabilized in the 1980s, precautionary annual yields of about 200000 tonnes could have been sustained. Our analysis demonstrates the value of incorporating prior knowledge to counteract shifting baseline effects on reference points and contemporary perceptions of historical stock status.

Description: The freshwater sponge Ephydatia muelleri and its Chlorella-like algal partner is an emerging model for studying animal: algal endosymbiosis. The sponge host is a tractable laboratory organism, and the symbiotic algae are easily cultured. We took advantage of these traits to interrogate questions about mechanisms that govern the establishment of durable intracellular partnerships between hosts and symbionts in facultative symbioses. We modified a classical experimental approach to discern the phagocytotic mechanisms that might be co-opted to permit persistent infections, and identified genes differentially expressed in sponges early in the establishment of endosymbiosis. We exposed algal-free E. muelleri to live native algal symbionts and potential food items (bacteria and native heat-killed algae), and performed RNA-Seq to compare patterns of gene expression among treatments. We found a relatively small but interesting suite of genes that are differentially expressed in the host exposed to live algal symbionts, and a larger number of genes triggered by host exposure to heat-killed algae. The upregulated genes in sponges exposed to live algal symbionts were mostly involved in endocytosis, ion transport, metabolic processes, vesicle-mediated transport, and oxidation–reduction. One of the host genes, an ATP-Binding Cassette transporter that is downregulated in response to live algal symbionts, was further evaluated for its possible role in the establishment of the symbiosis. We discuss the gene expression profiles associated with host responses to living algal cells in the context of conditions necessary for long-term residency within host cells by phototrophic symbionts as well as the genetic responses to sponge phagocytosis and immune-driven pathways.

Description: Marine sponges play a major ecological role in recycling resources on coral reef ecosystems. The cycling of resources may largely depend on the stability of the host-microbiome interactions and their susceptibility to altered environmental conditions. Given the current coral to algal phase shift on coral reefs, we investigated whether the sponge-associated bacterial communities of four sponge species, with either high or low microbial abundances (HMA and LMA), remain stable at two reefs sites with different coral to algae cover ratios. Additionally, we assessed the bacterial community composition of two of these sponge species before and after a reciprocal transplantation experiment between the sites. An overall stable bacterial community composition was maintained across the two sites in all sponge species, with a high degree of host-specificity. Furthermore, the core bacterial communities of the sponges remained stable also after a 21-day transplantation period, although a minor shift was observed in less abundant taxa (〈 1%). Our findings support the conclusion that host identity and HMA-LMA status are stronger traits in shaping bacterial community composition than habitat. Nevertheless, long-term microbial monitoring of sponges along with benthic biomass and water quality assessments are needed for identifying ecosystem tolerance ranges and tipping points in ongoing coral reef phase shifts.

Description: The maximum sustainable yield (MSY) concept is widely considered to be outdated and misleading. In response, fisheries scientists have developed models that often diverge radically from the first operational version of the concept. We show that the original MSY concept was deeply rooted in ecology and that going back to that version would be beneficial for fisheries, not least because the various substitutes have not served us well.

Description: Food webs are central entities mediating processes and external pressures in marine ecosystems. They are essential to understand and predict ecosystem dynamics and provision of ecosystem services. Paradoxically, utilization of food web knowledge in marine environmental conservation and resource management is limited. To better understand the use of knowledge and barriers to incorporation in management, we assess its application related to the management of eutrophication, chemical contamination, fish stocks, and non-indigenous species. We focus on the Baltic, a severely impacted, but also intensely studied and actively managed semi-enclosed sea. Our assessment shows food web processes playing a central role in all four areas, but application varies strongly, from formalized integration in management decisions, to support in selecting indicators and setting threshold values, to informal knowledge explaining ecosystem dynamics and management performance. Barriers for integration are complexity of involved ecological processes and that management frameworks are not designed to handle such information. We provide a categorization of the multi-faceted uses of food web knowledge and benefits of future incorporation in management, especially moving towards ecosystem-based approaches as guiding principle in present marine policies and directives. We close with perspectives on research needs to support this move considering global and regional change.

Description: Many coastal oceans experience not only increased loads of nutrients but also changes in the stoichiometry of nutrient supply. Excess supply of nitrogen and stable or decreased supply of silicon lower silicon to nitrogen (Si:N) ratios, which may decrease diatom proportion in phytoplankton. To examine how Si:N ratios affect plankton community composition and food web structure, we performed a mesocosm experiment where we manipulated Si:N ratios and copepod abundance in a Baltic Sea plankton community. In high Si:N treatments, diatoms dominated. Some of them were likely spared from grazing unexpectedly resulting in higher diatom biomass under high copepod grazing. With declining Si:N ratios, dinoflagellates became more abundant under low and picoplankton under high copepod grazing. This altered plankton food web structure: under high Si:N ratios, edible diatoms were directly accessible food for copepods, while under low Si:N ratios, microzooplankton and phago-mixotrophs (mixoplankton) were a more important food source for mesograzers. The response of copepods to changes in the phytoplankton community was complex and copepod density-dependent. We suggest that declining Si:N ratios favor microzoo- and mixoplankton leading to increased complexity of planktonic food webs. Consequences on higher trophic levels will, however, likely be moderated by edibility, nutritional value or toxicity of dominant phytoplankton species.

Description: Quantification and attribution of the food web changes associated with the invasion of non-indigenous species in the marine realm often remain a challenge. One of the pelagic non-indigenous species of concern in the recent history of aquatic bioinvasions is the predatory cladoceran Cercopagis pengoi, which invaded the Baltic Sea in the early 1990s. While several studies have reported immediate declines in abundances of its potential prey, the long-term effects of C. pengoi on the food webs remain to be examined. Based on the long-term time series (1968–2018) in the Gulf of Riga (Baltic Sea), we found significant declines in abundance of the cladoceran Pleopis spp. and copepod Eurytemora affinis by 90 and 80%, respectively, are associated with the invasion of C. pengoi as well as significant alterations in seasonal abundance patterns of Pleopis spp., E. affinis and cladoceran Bosmina spp. The invasion of the non-indigenous predator has led to the changed prey abundance–temperature relationships. Special caution was taken in data preprocessing, to minimize the likelihood that observed changes in the zooplankton prey could be associated with factors other than the invasion of C. pengoi, such as temperature and storminess.

Description: The copepod Acartia tonsa is a key component of a wide range of marine ecosystems, linking energy transfer from phytoplankton to higher trophic levels, and has a central role in productivity and biogeochemistry. The interaction of end-of-century global warming and ocean acidification scenarios with testing moderate temperature effects on a seminatural copepod community is needed to understand future community functioning. Here, we deployed a mesocosm experimental set-up with a full factorial design using two temperatures (13°C and 19°C) crossed with a pCO2 gradient ranging from ambient (550 μatm) to 3000 μatm. We used the natural bacteria, phyto- and microzooplankton species composition and biomass of the Kiel Bight and tested the response of A. tonsa development, carbon growth, mortality, size and condition. The tested traits were differently affected by the interaction of temperature and acidification. Ocean acidification increased development, carbon growth, size and mortality under the warming scenario of 19°C. At 13°C mortality rates decreased, while carbon growth, size and condition increased with acidification. We conclude from our experimental approach that a single species shows a variety of responses depending on the focal functional trait. Trait-specific mesozooplankton responses need to be further investigated and compared between geographical regions, seasons and taxonomic groups.

Description: Dinoflagellates possess many cellular characteristics with unresolved evolutionary histories. These include nuclei with greatly expanded genomes and chromatin packaged using histone-like proteins and dinoflagellate-viral nucleoproteins instead of histones, highly reduced mitochondrial genomes with extensive RNA editing, a mix of photosynthetic and cryptic secondary plastids, and tertiary plastids. Resolving the evolutionary origin of these traits requires understanding their ancestral states and early intermediates. Several early-branching dinoflagellate lineages are good candidates for such reconstruction, however these cells tend to be delicate and environmentally sparse, complicating such analyses. Here, we employ transcriptome sequencing from manually isolated and microscopically documented cells to resolve the placement of two cells of one such genus, Abedinium, collected by remotely operated vehicle in deep waters off the coast of Monterey Bay, CA. One cell corresponds to the only described species, Abedinium dasypus, whereas the second cell is distinct and formally described as Abedinium folium, sp. nov. Abedinium has classically been assigned to the early-branching dinoflagellate subgroup Noctilucales, which is weakly supported by phylogenetic analyses of small subunit ribosomal RNA, the single characterized gene from any member of the order. However, an analysis based on 221 proteins from the transcriptome places Abedinium as a distinct lineage, separate from and basal to Noctilucales and the rest of the core dinoflagellates. The transcriptome also contains evidence of a cryptic plastid functioning in the biosynthesis of isoprenoids, iron–sulfur clusters, and heme, a mitochondrial genome with all three expected protein-coding genes (cob, cox1, and cox3), and the presence of some but not all dinoflagellate-specific chromatin packaging proteins.

Description: The Law of the Sea as well as regional and national laws and agreements require exploited populations or stocks to be managed so that they can produce maximum sustainable yields. However, exploitation level and stock status are unknown for most stocks because the data required for full stock assessments are missing. This study presents a new method (AMSY) that estimates relative population size when no catch data are available using time-series of catch-per-unit-effort or other relative abundance indices as the main input. AMSY predictions for relative stock size were not significantly different from the “true” values when compared with simulated data. Also, they were not significantly different from relative stock size estimated by data-rich models in 88% of the comparisons within 140 real stocks. Application of AMSY to 38 data-poor stocks showed the suitability of the method and led to the first assessments for 23 species. Given the lack of catch data as input, AMSY estimates of exploitation come with wide margins of uncertainty which may not be suitable for management. However, AMSY seems to be well suited for estimating productivity as well as relative stock size and may, therefore, aid in the management of data-poor stocks.

Description: Diurnal vertical migration (DVM) is a widespread phenomenon in the upper ocean, but it remains unclear to what degree it also involves passively transported micro- and meso-zooplankton. These organisms are difficult to monitor by in situ sensing and observations from discrete samples are often inconclusive. Prime examples of such ambiguity are planktonic foraminifera, where contradictory evidence for DVM continues to cast doubt on the stability of species vertical habitats, which introduces uncertainties in geochemical proxy interpretation. To provide a robust answer, we carried out highly replicated randomised sampling with 41 vertically resolved plankton net hauls taken within 26 hours in a confined area of 400 km2 in the tropical North Atlantic, where DVM in larger plankton occurs. Manual enumeration of planktonic foraminifera cell density consistently reveals the highest total cell concentrations in the surface mixed layer (top 50 m) and analysis of cell density in seven individual species representing different shell sizes, life strategies and presumed depth habitats reveals consistent vertical habitats not changing over the 26 hours sampling period. These observations robustly reject the existence of DVM in planktonic foraminifera in a setting where DVM occurs in other organisms.

Description: Sponges are one of the most dominant organisms in marine ecosystems. One reason for their success is their association with microorganisms that are besides the host itself responsible for the chemical defence. Sponge abundances have been increasing on coral reefs in the Western Indian Ocean (WIO) and are predicted to increase further with rising anthropogenic impacts on coral reefs. However, there is a paucity of information on chemical ecology of sponges from the WIO and their prokaryotic community composition. We used a combination of Illumina sequencing and a predictive metagenomic analysis to (i) assess the prokaryotic community composition of sponges from Zanzibar, (ii) predict the presence of KEGG metabolic pathways responsible for bioactive compound production and (iii) relate their presence to the degree of observed chemical defence in their respective sponge host. We found that sponges from Zanzibar host diverse prokaryotic communities that are host species-specific. Sponge-species and respective specimens that showed strong chemical defences in previous studies were also predicted to be highly enriched in various pathways responsible for secondary metabolite production. Hence, the combined sequencing and predictive metagenomic approach proved to be a useful indicator for the metabolic potential of sponge holobionts

Description: Transparent exopolymer particles (TEP) are polysaccharide-rich microgels that are prevalent in the marine environment and have important roles in the aggregation of organic matter and carbon export from the euphotic zone. TEP are readily colonized by bacteria and utilized by specialized taxa, such as Alteromonadaceae. However, bacterial community composition specifically attached to natural TEP remains largely unknown. In this study, we isolated individual TEP from Plymouth Sound (UK) and performed DNA sequencing of the TEP-attached bacterial communities. We also sampled the cognate bulk seawater total bacterial communities for comparison. The bacterial communities associated with individual TEP showed distinct differences compared to the total bulk bacterioplankton communities, with Alteromonadaceae significantly more abundant on TEP. The TEP-associated Alteromonadaceae consisted of two operational taxonomic units that were closely related to Marinobacter and Glaciecola, both previously associated with biogenic aggregates and microgel-rich habitats. This study provides novel insight into marine bacterial–microgel interactions.

Description: There is a recognized need for new methods with modest data requirements to provide preliminary estimates of stock status for data-limited stocks (e.g. Rudd and Thorson, 2018). Froese et al. (2018) provide such a method, which derives estimates of relative stock size from length frequency (LF) data of exploited stocks. They show that their length-based Bayesian biomass estimation method (LBB) can reproduce the “true” parameters used in simulated data and can approximate the relative stock size as estimated independently by more data-demanding methods in 34 real stocks. However, in a comment on LBB, Hordyk et al. (2019) claim (i) that the master equation of LBB is incomplete because it does not correct for the pile-up effect caused by aggregating length measurements into length classes or “bins”, (ii) that LBB is highly sensitive to equilibrium assumptions and wrongly uses maximum observed length (Lmax) for guidance in setting a prior for the estimation of asymptotic length (Linf), and (iii) that the default prior used by LBB for the ratio between natural mortality and somatic growth rate (M/K) of 1.5 (SD = 0.15) is inadequate for many exploited species. These comments are addressed below

Description: On 25 December 2016, a Mw 7.6 earthquake broke a portion of the Southern Chilean subduction zone south of Chiloé Island, located in the central part of the Mw 9.5 1960 Valdivia earthquake. This region is characterized by repeated earthquakes in 1960 and historical times with very sparse interseismic activity due to the subduction of a young (~15 Ma), and therefore hot, oceanic plate. We estimate the co-seismic slip distribution based on a kinematic finite fault source model, and through joint inversion of teleseismic body waves and strong motion data. The coseismic slip model yields a total seismic moment of 3.94×1020 Nm that occurred over ~30 s, with the rupture propagating mainly downdip, reaching a peak-slip of ~4.2 m. Regional moment tensor inversion of stronger aftershocks reveals thrust type faulting at depths of the plate interface. The fore- and aftershock seismicity is mostly related to the subduction interface with sparse seismicity in the overriding crust. The 2016 Chiloé event broke a region with increased locking and most likely broke an asperity of the 1960 earthquake. The updip limit of the main event, aftershocks, foreshocks and interseismic activity are spatially similar, located ~15 km offshore and parallel to Chiloé Islands west coast. The coseismic slip model of the 2016 Chiloé earthquake suggests a peak slip of 4.2 m that locally exceeds the 3.38 m slip deficit that has accumulated since 1960. Therefore, the 2016 Chiloé earthquake possibly released strain that has built up prior to the 1960 Valdivia earthquake.

Description: We used a molecular phylogenetic approach to investigate species delimitation and diversification in the northeastern Atlantic and Mediterranean musseldrills of the Ocinebrina aciculata complex, based on molecular data from topotypical material of many of the nominal taxa. The complex is shown to consist of at least five species: Ocinebrina aciculata (Lamarck, 1822) from the Atlantic and western Mediterranean; O. cf. corallina (Scacchi, 1836) from the central Mediterranean Sea; O. reinai Bonomolo & Crocetta, 2012 from the Tyrrhenian Sea; O. corallinoides Pallary, 1912 from the Gulf of Gabès; and O. aegeensis n. sp. currently known from the Aegean Sea only. The new species is differentiated from the other taxa by very subtle morphological diagnostic features, although it is clearly identified by genetic distance and apomorphic DNA-sequence characters. The identity of Murex corallinus Scacchi, 1836 (type species of Ocinebrina Jousseaume, 1880) could not be defined with certainty, pending genetic comparison of specimens of the â €? large Tyrrhenian morphotype' (corresponding to the neotype, but not assayed herein) with the assayed â €? small morphotype'.

Description: Marine sponges are early-branched metazoans known to harbor dense and diverse microbial communities. Yet the role of the so far uncultivable alphaproteobacterial lineages that populate these sessile invertebrates remains unclear. We applied a sequence composition-dependent binning approach to assemble one Rhodospirillaceae genome from the Spongia officinalis microbial metagenome and contrast its functional features with those of closely related sponge-associated and free-living genomes. Both symbiotic and free-living Rhodospirillaceae shared a suite of common features, possessing versatile carbon, nitrogen, sulfur and phosphorus metabolisms. Symbiotic genomes could be distinguished from their free-living counterparts by the lack of chemotaxis and motility traits, enrichment of genes required for the uptake and utilization of organic sulfur compounds—particularly taurine—, higher diversity and abundance of ABC transporters, and a distinct repertoire of genes involved in natural product biosynthesis, plasmid stability, cell detoxification and oxidative stress remediation. These sessile symbionts may more effectively contribute to host fitness via nutrient exchange, and also host detoxification and chemical defense. Considering the worldwide occurrence and high diversity of sponge-associated Rhodospirillaceae verified here using a tailored in silico approach, we suggest that these organisms are not only relevant to holobiont homeostasis but also to nutrient cycling in benthic ecosystems.

Description: Continental hyperextension during magma-poor rifting at the Deep Galicia Margin is characterised by a complex pattern of faulting, thin continental fault blocks, and the serpentinisation, with local exhumation, of mantle peridotites along the S-reflector, interpreted as a detachment surface. In order to understand fully the evolution of these features, it is important to image seismically the structure and to model the velocity structure to the greatest resolution possible. Travel-time tomography models have revealed the long-wavelength velocity structure of this hyperextended domain, but are often insufficient to match accurately the short-wavelength structure observed in reflection seismic imaging. Here we demonstrate the application of two-dimensional (2D) time-domain acoustic full-waveform inversion to deep water seismic data collected at the Deep Galicia Margin, in order to attain a high resolution velocity model of continental hyperextension. We have used several quality assurance procedures to assess the velocity model, including comparison of the observed and modelled waveforms, checkerboard tests, testing of parameter and inversion strategy, and comparison with the migrated reflection image. Our final model exhibits an increase in the resolution of subsurface velocities, with particular improvement observed in the westernmost continental fault blocks, with a clear rotation of the velocity field to match steeply dipping reflectors. Across the S-reflector there is a sharpening in the velocity contrast, with lower velocities beneath S indicative of preferential mantle serpentinisation. This study supports the hypothesis that normal faulting acts to hydrate the upper mantle peridotite, observed as a systematic decrease in seismic velocities, consistent with increased serpentinisation. Our results confirm the feasibility of applying the full-waveform inversion method to sparse, deep water crustal datasets.

Description: Overfishing and rapid environmental shifts pose severe challenges to the resilience and viability of marine fish populations. To develop and implement measures that enhance species’ adaptive potential to cope with those pressures while, at the same time, ensuring sustainable exploitation rates is part of the central goal of fisheries management. Here, we argue that a combination of biophysical modelling and population genomic assessments offer ideal management tools to define stocks, their physical connectivity and ultimately, their short-term adaptive potential. To date, biophysical modelling has often been confined to fisheries ecology whereas evolutionary hypotheses remain rarely considered. When identified, connectivity patterns are seldom explored to understand the evolution and distribution of adaptive genetic variation, a proxy for species’ evolutionary potential. Here, we describe a framework that expands on the conventional seascape genetics approach by using biophysical modelling and population genomics. The goals are to identify connectivity patterns and selective pressures, as well as putative adaptive variants directly responding to the selective pressures and, ultimately, link both to define testable hypotheses over species response to shifting ecological conditions and overexploitation.

Description: This study presents a new method (LBB) for the analysis of length frequency data from commercial catches. LBB works for species that grow throughout their lives, such as most commercially-important fish and invertebrates, and requires no input in addition to length frequency data. It estimates asymptotic length, length at first capture, relative natural mortality, and relative fishing mortality. Standard fisheries equations can then be used to approximate current exploited biomass relative to unexploited biomass. In addition, these parameters allow the estimation of length at first capture that would maximize catch and biomass for a given fishing effort, and estimation of a proxy for the relative biomass capable of producing maximum sustainable yields. Relative biomass estimates of LBB were not significantly different from the “true” values in simulated data and were similar to independent estimates from full stock assessments. LBB also presents a new indicator for assessing whether an observed size structure is indicative of a healthy stock. LBB results will obviously be misleading if the length frequency data do not represent the size composition of the exploited size range of the stock or if length frequencies resulting from the interplay of growth and mortality are masked by strong recruitment pulses.

Description: In this study, the complex frequency-shifted perfectly matched layer (CFS-PML) in stretching Cartesian coordinates, is successfully applied to three-dimensional (3D) frequency-domain marine controlled-source electromagnetic (CSEM) field modelling. The Dirichlet boundary, which is usually used within the traditional framework of EM modeling algorithms, assumes the electric or magnetic field values are zero at the boundaries. This requires the boundaries be sufficiently far away from the sources in the area of interest. To mitigate the boundary artifacts, a large modelling area may be necessary even though cell sizes are allowed to grow toward the boundaries due to the diffusion of the electromagnetic wave propagation. Compared with the conventional Dirichlet boundary, the PML boundary is preferred as the modelling area of interest could be restricted to the target region and only a few absorbing layers surrounding can effectively depress the artificial boundary effect without losing the numerical accuracy. Furthermore, for joint inversion of seismic and marine CSEM data, if we used the PML for CSEM field simulation instead of the conventional Dirichlet, the modeling area for these two different geophysical data collected from the same survey area could the same, which is convenient for joint inversion grid matching. We apply the CFS-PML boundary to 3D marine CSEM modelling by using the staggered finite-difference (SFD) discretization. Numerical test indicates that the modeling algorithm using the CFS-PML also shows good accuracy compared to the Dirichlet. Furthermore, the modeling algorithm using the CFS-PML shows advantages in computational time and memory saving than that using the Dirichlet boundary. For the 3D example in this study, the memory saving using the PML is nearly 42 % and the time saving is around 48% compared to using the Dirichlet.

Description: The BONUS symposium “Science delivery for sustainable use of the Baltic Sea living resources” held in Tallinn, Estonia, in October 2017 was an opportunity for the presentation and discussion of 107 papers that examined the state and dynamics of living resources of the Baltic Sea, and associated management challenges. The symposium included a half-day stakeholder panel discussion that addressed the main challenges related to sustainable management and matching research and policy/management needs. Based on the five symposium papers published in this Special Issue as well as the stakeholder panel discussion, it can be concluded that (i) new observations about the feeding ecology of clupeids supports a more complete understanding of trophic interactions in the pelagic realm and improved calibration of ecosystem models, (ii) to safequard sustainable and diverse fisheries resources, one should take into account the specific local characteristics of the fish community, (iii) to safeguard sustainable use of marine resources and mitigate cross-sectoral and transboundary conflicts, a risk-based approach should be adopted, and (iv) incorporation of scientific advice into management faces several obstacles including the reality that not all readily available knowledge is currently being incorporated into the decision-making process.

Description: Fisheries and marine ecosystem-based management requires a holistic understanding of the dynamics of fish communities and their responses to changes in environmental conditions. Environmental conditions can simultaneously shape the spatial distribution and the temporal dynamics of a population, which together can trigger changes in the functional structure of communities. Here, we developed a comprehensive framework based on complementary multivariate statistical methodologies to simultaneously investigate the effects of environmental conditions on the spatial, temporal and functional dynamics of species assemblages. The framework is tested using survey data collected during more than 4000 fisheries hauls over the Baltic Sea between 2001 and 2016. The approach revealed the Baltic fish community to be structured into three sub-assemblages along a strong and temporally stable salinity gradient decreasing from West to the East. Additionally, we highlight a mismatch between species and functional richness associated with a lower functional redundancy in the Baltic Proper compared with other sub-areas, suggesting an ecosystem more susceptible to external pressures. Based on a large dataset of community data analysed in an innovative and comprehensive way, we could disentangle the effects of environmental changes on the structure of biotic communities-key information for the management and conservation of ecosystems.

Description: The European and American eels spawn in the subtropical convergence zone (STCZ) in the Sargasso Sea, a dynamic and relatively productive area that is strongly influenced by front and eddy formations and subducted high-saline water masses. To understand how the physical and biological environments may affect the early life history of eels, we conducted a detailed bio-physical investigation of the water column at a site of high eel larvae abundance. Diel measurements and sampling in the upper 300 m revealed strong variations in hydrographic conditions and mean depths of different taxonomic groups; however, characteristics patterns of distribution were apparent. Most species showed diel vertical migrations, ascending about 20-30 m at night, whereas examples of night-time downward migration were also seen. European eel larvae were among the species showing more extensive diel vertical migration: their population mean depth changed from 160 m at day to 100 m at night where abundance peaked at 45 m depth. Distribution and migration of eel larvae corresponded to patterns observed for small hydrozoans, supporting a proposed predator-prey linkage. The study demonstrates the diverse and vertically strongly structured plankton community of STCZ where larvae of eel and other fish find a wide range of potential niches.

Description: An indoor mesocosm experiment was carried out to investigate the combined effects of ocean acidification and warming on the species composition and biogeochemical element cycling during a winter/spring bloom with a natural phytoplankton assemblage from the Kiel fjord, Germany. The experimental setup consisted of a "Control" (ambient temperature of similar to 4.8 degrees C and similar to 535 +/- 25 mu atm pCO(2)), a "High-CO2" treatment (ambient temperature and initially 1020 +/- 45 mu atm pCO(2)) and a "Greenhouse" treatment (similar to 8.5 degrees C and initially 990 +/- 60 mu atm pCO(2)). Nutrient replete conditions prevailed at the beginning of the experiment and light was provided at in situ levels upon reaching pCO(2) target levels. A diatom-dominated bloom developed in all treatments with Skeletonema costatum as the dominant species but with an increased abundance and biomass contribution of larger diatom species in the Greenhouse treatment. Conditions in the Greenhouse treatment accelerated bloom development with faster utilization of inorganic nutrients and an earlier peak in phytoplankton biomass compared to the Control and High CO2 but no difference in maximum concentration of particulate organic matter (POM) between treatments. Loss of POM in the Greenhouse treatment, however, was twice as high as in the Control and High CO2 treatment at the end of the experiment, most likely due to an increased proportion of larger diatom species in that treatment. We hypothesize that the combination of warming and acidification can induce shifts in diatom species composition with potential feedbacks on biogeochemical element cycling.

Description: Light and nutrients are essential resources for phytoplankton growth and considered to shape the size structure and other morphometric traits (surface:volume ratio, deviation from spherical shape) of phytoplankton communities. If morphometric traits influence the growth and resource use, shifts by one of the two factors should influence the capability to utilize the other factor. We performed a two-step experiment, where a natural phytoplankton community was first exposed to three different light levels (supposed to be limiting, saturating and slightly inhibiting for the majority of species) and grown until stationary phase. Then, the pre-conditioned communities were split into two nutrient treatments (control and saturating nutrient pulse) and again grown until stationary phase under the medium light intensity. During the experimental light phase, community mean cell-size increased with light, but surface:volume ratio and deviation from spherical shape decreased. Moreover, in response to the following nutrient pulse, the low light pre-conditioned communities showed the highest initial growth rates in response to the nutrient pulse. The high light pre-conditioned communities showed the highest conversion of the nutrient pulse into biomass during the stationary phase. These results demonstrate how the imprint of one environmental factor on trait distribution influences the ability to cope with another.

Description: Aurelia aurita (Linneaus, 1758) is a cosmopolitan scyphozoan, probably the most investigated jellyfish in temperate and highly productive coastal ecosystems. Despite a prominent top-down control in plankton food webs, a mechanistic understanding of A. aurita population dynamics and trophic interactions has been barely addressed. Here we develop a food web dynamic model to assess A. aurita role in the seasonal plankton dynamics of the Kiel Fjord, southwestern Baltic Sea. The model couples low trophic level dynamics, based on a classical Nutrient Phytoplankton Zooplankton Detritus (NPZD) model, to a stage-resolved copepod model (referencing Pseudocalanus sp.) and a jellyfish model (A. aurita ephyra and medusa) as consumers and predators, respectively. Simulations showed the relevance of high abundances of A. aurita, which appear related with warm winter temperatures, promoting a shift from a copepod-dominated food web to a ciliate and medusa dominated one. The model captured the intraspecific competition triggered by the medusae abundance and characterized by a negative relationship between population density and individual size/weight. Our results provide a mechanistic understanding of an emergent trait such as size shaping the food web functioning, driving predation rates and population dynamics of A. aurita, driving its sexual reproductive strategy at the end of the pelagic phase.

Description: Atlantic herring (Clupea harengus) is a benthic spawner, therefore its eggs are prone to encounter different water conditions during embryonic development, with bottom waters often depleted of oxygen and enriched in CO2. Some Atlantic herring spawning grounds are predicted to be highly affected by ongoing Ocean Acidification and Warming with water temperature increasing by up to +3°C and CO2 levels reaching ca. 1000 μatm (RCP 8.5). Although many studies investigated the effects of high levels of CO2 on the embryonic development of Atlantic herring, little is known about the combination of temperature and ecologically relevant levels of CO2. In this study, we investigated the effects of Ocean Acidification and Warming on embryonic metabolic and developmental performance such as mitochondrial function, respiration, hatching success (HS) and growth in Atlantic herring from the Oslo Fjord, one of the spawning grounds predicted to be greatly affected by climate change. Fertilized eggs were incubated under combinations of two PCO2 conditions (400 μatm and 1100 μatm) and three temperatures (6, 10 and 14°C), which correspond to current and end-of-the-century conditions. We analysed HS, oxygen consumption (MO2) and mitochondrial function of embryos as well as larval length at hatch. The capacity of the electron transport system (ETS) increased with temperature, reaching a plateau at 14°C, where the contribution of Complex I to the ETS declined in favour of Complex II. This relative shift was coupled with a dramatic increase in MO2 at 14°C. HS was high under ambient spawning conditions (6–10°C), but decreased at 14°C and hatched larvae at this temperature were smaller. Elevated PCO2 increased larval malformations, indicating sub-lethal effects. These results indicate that energetic limitations due to thermally affected mitochondria and higher energy demand for maintenance occur at the expense of embryonic development and growth.

Description: The ocean load in glacial isostatic adjustment (GIA) modelling is represented by the so-called sea level equation (SLE). The SLE describes the mass redistribution of water between ice sheets and oceans on a deforming Earth. Despite various teams independently investigating GIA, there has been no systematic intercomparison among the numerical solvers of the SLE through which the methods may be validated. The goal of this paper is to present a series of synthetic examples designed for testing and comparing the numerical implementations of the SLE in GIA modelling. The 10 numerical codes tested combine various temporal and spatial parametrizations. The time-domain or Laplace-domain discretizations are used to solve the SLE through time, while spherical harmonics, finite differences or finite elements parametrize the GIA-related field variables spatially. The surface ice-water load and solid Earth’s topography are represented spatially either on an equiangular grid, a Gauss–Legendre or an equiarea grid with icosahedron-shaped spherical pixels. Comparisons are made in a series of five benchmark examples with an increasing degree of complexity. Due to the complexity of the SLE, there is no analytical solution to it. The accuracy of the numerical implementations is therefore assessed by the differences of the individual solutions with respect to a reference solution. While the benchmark study does not result in GIA predictions for a realistic loading scenario, we establish a set of agreed-upon results that can be extended in the future by including more complex case studies, such as solutions with realistic loading scenarios, the rotational feedback in the linear-momentum equation, and by considering a 3-D viscosity structure of the Earth’s mantle. The test computations performed so far show very good agreement between the individual results and their ability to capture the main features of sea-surface variation and the surface vertical displacement. The differences found can often be attributed to the different approximations inherent in the various algorithms. This shows the accuracy that can be expected from different implementations of the SLE, which helps to assess differences noted in the literature between predictions for realistic loading cases.

Description: Seafloor massive sulphides (SMSs) are regarded as a potential future resource to satisfy the growing global demand of metals including copper, zinc and gold. Aside from mining and retrieving profitable amounts of massive sulphides from the seafloor, the present challenge is to detect and delineate significant SMS accumulations, which are generally located near mid-ocean ridges and along submarine volcanic arc and backarc spreading centres. Currently, several geophysical technologies are being developed to detect and quantify SMS occurrences that often exhibit measurable contrasts in their physical parameters compared to the surrounding host rock. Here, we use a short, fixed-offset controlled source electromagnetic (CSEM) system and a coincident-loop transient electromagnetic (TEM) system, which in theory allow the detection of SMS in the shallow seafloor due to a significant electrical conductivity contrast to their surroundings. In 2016, CSEM and TEM experiments were carried out at several locations near the Trans- Atlantic Geotraverse hydrothermal field to investigate shallow occurrences of massive sulphides below the seafloor. Measurements were conducted in an area that contains distinct SMS sites located several kilometres off-axis from the Mid-Atlantic ridge, some of which are still connected to hydrothermal activity and others where hydrothermal activity has ceased. Based on the quality of the acquired data, both experiments were operationally successful. However, the data analysis indicates bias caused by three-dimensional (3D) effects of the rough bathymetry in the study area and, thus, data interpretation remains challenging. Therefore, we study the influence of 3D bathymetry for marine CSEM and TEM experiments, focusing on shallow 3D conductors located beneath mound-like structures.We analyse synthetic inversion models for attributes associated with 3D distortions of CSEM and TEM data that are not sufficiently accounted for in conventional 1D (TEM) and 2D (CSEM) interpretation schemes. Before an adequate quantification of SMS in the region is feasible, these 3D effects need to be studied to avoid over/underestimation of SMS using the acquired EM data. The sensitivity of CSEM and TEM to bathymetry is investigated by means of 3D forward modelling, followed by 1D (TEM) and 2D (CSEM) inversion of the synthetic data using realistic error conditions. Subsequently, inversion models of the synthetic 3D data are analysed and compared to models derived from the measured data to illustrate that 3D distortions are evident in the recorded data sets.

Description: Infections arising from multidrug-resistant pathogenic bacteria are spreading rapidly throughout the world and threaten to become untreatable. The origins of resistance are numerous and complex, but one underlying factor is the capacity of bacteria to rapidly export drugs through the intrinsic activity of efflux pumps. In this Review, we describe recent advances that have increased our understanding of the structures and molecular mechanisms of multidrug efflux pumps in bacteria. Clinical and laboratory data indicate that efflux pumps function not only in the drug extrusion process but also in virulence and the adaptive responses that contribute to antimicrobial resistance during infection. The emerging picture of the structure, function and regulation of efflux pumps suggests opportunities for countering their activities.

Description: The transition from benthos to plankton requires multiple adaptations, yet so far it remains unclear how these are acquired in the course of the transition. To investigate this process, we analyzed the genetic diversity and distribution patterns of a group of foraminifera of the genus Bolivina with a tychopelagic mode of life (same species occurring both in benthos and plankton). We assembled a global sequence data set for this group from single-cell DNA extractions and occurrences in metabarcodes from pelagic environmental samples. The pelagic sequences all cluster within a single monophyletic clade within Bolivina. This clade harbors three distinct genetic lineages, which are associated with incipient morphological differentiation. All lineages occur in the plankton and benthos, but only one lineage exhibits no limit to offshore dispersal and has been shown to grow in the plankton. These observations indicate that the emergence of buoyancy regulation within the clade preceded the evolution of pelagic feeding and that the evolution of both traits was not channeled into a full transition into the plankton. We infer that in foraminifera, colonization of the planktonic niche may occur by sequential cooptation of independently acquired traits, with holoplanktonic species being recruited from tychopelagic ancestors

Description: The sea surface microlayer (SML) is the uppermost layer of the water column that links the ocean and atmosphere. It accumulates a variety of biogenic surface-active and buoyant substances, including gelatinous material, such as transparent exopolymer particles (TEP) and Coomassie stainable particles (CSP), potentially affecting air–sea exchange processes. Here, we studied the influence of the annual cycle of phytoplankton production on organic matter (OM) accumulation in the SML relative to the subsurface water (SSW). Sampling was performed monthly from April 2012 to November 2013 at the Boknis Eck Time Series Station (Baltic Sea). For SML sampling, we used the Garrett screen, while SSW samples were collected by Niskin bottles at 1 m depth. Samples were analyzed for carbohydrates, amino acids, TEP, CSP, chlorophyll a (SSW only) and bacterial abundance. Our data showed that the SML reflected the SSW during most parts of the year, with changes mainly responding to bloom formation and decay. OM composition during phytoplankton blooms clearly differed from periods of higher bacterial abundance. Of all components investigated, only the enrichment of total carbohydrates in the SML was inversely related to the wind speed indicating that wind-driven mixing also affected the accumulation of OM in the SML during our study.

Description: Naturally produced polybrominated diphenyl ethers (PBDEs) pervade the marine environment and structurally resemble toxic man-made brominated flame retardants. PBDEs bioaccumulate in marine animals and are likely transferred to the human food chain. However, the biogenic basis for PBDE production in one of their most prolific sources, marine sponges of the order Dysideidae, remains unidentified. Here, we report the discovery of PBDE biosynthetic gene clusters within sponge-microbiome-associated cyanobacterial endosymbionts through the use of an unbiased metagenome-mining approach. Using expression of PBDE biosynthetic genes in heterologous cyanobacterial hosts, we correlate the structural diversity of naturally produced PBDEs to modifications within PBDE biosynthetic gene clusters in multiple sponge holobionts. Our results establish the genetic and molecular foundation for the production of PBDEs in one of the most abundant natural sources of these molecules, further setting the stage for a metagenomic-based inventory of other PBDE sources in the marine environment.

Description: Processes linked with the genesis, evolution and emplacement of silicic complexes in arcs are still poorly constrained. Of particular interest are the depth of magma production, the relative contribution of crystal fractionation versus crustal partial melting and the timescales involved. The Soufrière Volcanic Complex (SVC) on St Lucia is one of the largest silicic centres in the Lesser Antilles arc. Here we present the results of a detailed mineralogical study, including in situ Sr isotopes in plagioclase and in situ δ18O in dated zircons, of both SVC and Pre-SVC volcanic rocks to place constraints on the processes intrinsic to the development and evolution of the silicic complex. These data suggest that the production of silicic magma in the SVC occurs in two stages. The first stage involves differentiation of mafic magma by crustal assimilation and mineral fractionation in the middle–lower crust of the arc to produce magmas with intermediate compositions. These intermediate magmas are water-rich (∼7 wt %) and have high 87Sr/86Sr, Ba, Sr and La/Sm (∼5) compared with Pre-SVC lavas. Near-constant trace element and isotopic compositions throughout the SVC lifespan indicate that the same process was persistent over the last 600 kyr. In the second stage, the intermediate magmas are transferred to a shallower and more differentiated chamber (∼6 km depth). During ascent, any crystals or xenocrysts residual from stage one in the deeper chamber become fully resorbed and the magma crystallizes calcic amphibole microphenocrysts, followed by anorthite-rich plagioclase close to or at the water saturation depth. During mixing upon recharge within the shallow chamber, anorthite-rich plagioclase from the recharging magma is partially resorbed; so are the crystals in equilibrium with the resident differentiated magma. The recharge event probably causes chamber-wide convection. Mixing is thought to trigger eruption of the silicic complex magmas.

Description: Eukaryotic phytoplankton exhibit an enormous species richness, displaying a range of phylogenetic, morphological and physiological diversity. Yet, until recently, very little was known about the diversity, genetic variation and evolutionary processes within species and populations. An approach to explore this diversity and to understand evolution of phytoplankton is to use population genetics as a conceptual framework and methodology. Here, we discuss the patterns, processes and questions that population genetic studies have revealed in eukaryotic phytoplankton. First, we describe the main biological processes generating genetic variation. We specifically discuss the importance of life-cycle complexity for genetic and phenotypic diversity and consider how such diversity can be maintained during blooms when rapid asexual proliferation dominates. Next, we explore how genetic diversity is partitioned over time and space, with a focus on the processes shaping this structure, in particular selection and genetic exchange. Our aim is also to show how population genetics can be used to make inferences about realized dispersal and sexual recombination, as these processes are so difficult to study directly. Finally, we highlight important open questions and suggest promising avenues for future studies that will be made possible by new sequencing technologies

Description: Marine spatial planning (MSP) is considered a valuable tool in the ecosystem-based management of marine areas. Predictive modelling may be applied in the MSP framework to obtain spatially explicit information about biodiversity patterns. The growing number of statistical approaches used for this purpose implies the urgent need for comparisons between different predictive techniques. In this study, we evaluated the performance of selected machine learning and regression-based methods that were applied for modelling fish community indices. We hypothesized that habitat features can influence fish assemblage and investigated the effect of environmental gradients on demersal fish diversity (species richness and Shannon–Weaver Index). We used fish data from the Baltic International Trawl Surveys (2001–2014) and maps of six potential predictors: bottom salinity, depth, seabed slope, growth season bottom temperature, seabed sediments and annual mean bottom current velocity. We compared the performance of six alternative modelling approaches: generalized linear models, generalized additive models, multivariate adaptive regression splines, support vector machines, boosted regression trees and random forests. We applied repeated 10-fold cross-validation, using accuracy as the measure of model quality. Finally, we selected random forest as the best performing algorithm and implemented it for the spatial prediction of fish diversity from the Baltic Proper to the Kattegat. To obtain information on the data reliability and confidence of the developed models, which are essential for MSP, we estimated the uncertainty of predictions with standard deviation of predictions obtained from all the trees in the ensemble random forest method. We showed how state-of-the-art predictive techniques, based on easily available data and simple Geographic Information System tools, can be used to obtain reliable spatial information about fish diversity. Our comparative work highlighted the potential of machine learning method to reduce prediction error in modelling of demersal fish diversity in the framework of MSP.

Description: The Mozambique Ridge, a prominent basement high in the southwestern Indian Ocean, consists of four major geomorphological segments associated with numerous phases of volcanic activity in the Lower Cretaceous. The nature and origin of the Mozambique Ridge have been intensely debated with one hypothesis suggesting a Large Igneous Province origin. High-resolution seismic reflection data reveal a large number of extrusion centres with a random distribution throughout the southern Mozambique Ridge and the nearby Transkei Rise. Intrabasement reflections emerge from the extrusion centres and are interpreted to represent massive lava flow sequences. Such lava flow sequences are characteristic of eruptions leading to the formation of continental and oceanic flood basalt provinces, hence supporting a Large Igneous Province origin of the Mozambique Ridge. We observe evidence for widespread post-sedimentary magmatic activity that we correlate with a southward propagation of the East African Rift System. Based on our volumetric analysis of the southern Mozambique Ridge we infer a rapid sequential emplacement between ∼131 and ∼125 Ma, which is similar to the short formation periods of other Large Igneous Provinces like the Agulhas Plateau.

Description: The Nifonea submarine volcano rises 1000 m above the seafloor of the Vate Trough back-arc basin behind the New Hebrides island arc. This large volcanic edifice has a caldera of ∼8 km diameter and is connected to two ∼20 km long volcanic rift zones in the back-arc basin. We present new chemical and isotope data for volcanic glasses and whole-rocks from both the volcano and its rift zones. Lavas from Nifonea volcano show an evolution towards more incompatible element enrichment, with the most enriched lavas being the youngest eruption products on the caldera floor. These are products of significant fractional crystallization, show minor contamination by hydrothermal fluids (〈0·3%) and reflect mixing of melts derived from depleted upper mantle and melts from an enriched source similar to those occurring in the North Fiji Basin. The enrichment in Nb of these lavas is comparable with that of some lavas from the New Hebrides island arc (e.g. Mota Lava island), where these coexist with typical island arc basalts. The lavas erupted along the rift zones in the Vate Trough back-arc basin are relatively depleted in incompatible elements, indicating melting of depleted upper mantle with a minor addition of a sediment-derived fluid. Our observations suggest that the mantle beneath Vate Trough is heterogeneous on a small scale (〈20 km) and that the occurrence of these enriched and fertile mantle portions has a stronger control on melting processes than the influx from the subducting slab, as all samples were recovered at a similar distance from the trench.

Description: The influence of winter on the selection of dominant taxa for the phytoplankton spring bloom was studied in batch culture experiments. Different natural phytoplankton assemblages from different phases of the temperate zone winter were exposed to varying periods of darkness (0, 6/7, 13 and 19 weeks) followed by a re-exposure to saturating light intensity for 14 days to experimentally simulate the onset of spring. The results showed that dark incubation has a strong effect on shaping the phytoplankton community composition. Many taxa disappeared in the absolute darkness. Dark survival ability might be an important contributing factor for the success of diatoms in spring. Different phytoplankton starting assemblages were dominated by the same bloom-forming diatoms, Skeletonema marinoi and Thalassosira spp., after dark incubation for only 6 weeks, irrespective of the high dissimilarities between phytoplankton communities. The growth capacity of surviving phytoplankton is almost unimpaired by darkness. Similar growth rates as that before darkness could be resumed for the surviving taxa with a potential lag time of 1–7 days dependent on taxon and the duration of darkness.

Description: The general positive effect of warmer winters on the abundance of small-sized zooplankton in the following spring and early summer has been reported from different parts of the Baltic Sea, but the mechanism of this link is not clear. Although causal links cannot be deduced with confidence from observational data, sufficiently detailed analyses can nevertheless provide insights to the potential mechanisms. We present an example of such an analysis, scrutinizing the effects of winter and spring hydroclimate on the abundance of small-sized dominant calanoid copepods (Eurytemora affinis and Acartia spp.), using data from 2080 zooplankton samples collected over 55 years (1957–2012) from a shallow coastal habitat (Pärnu Bay, Gulf of Riga) in the Baltic Sea. Our results indicated that the milder winters brought about higher abundances, and reduced seasonality of small-sized copepods, whereas ambient sea surface temperature (SST) mostly affected the relative abundance of adult stages. The sliding window correlation tests revealed temporal shifts in the effects of controlling variables: with the continuous increase in SST, the effect of winter temperature on the abundance of Acartia spp. weakened. In contrast, E. affinis was consistently affected by SST, but the effect of winter temperature was more pronounced during the period of on average colder winters.

Description: Marine sponges (phylum Porifera) are a diverse, phylogenetically deep-branching clade known for forming intimate partnerships with complex communities of microorganisms. To date, 16S rRNA gene sequencing studies have largely utilised different extraction and amplification methodologies to target the microbial communities of a limited number of sponge species, severely limiting comparative analyses of sponge microbial diversity and structure. Here, we provide an extensive and standardised dataset that will facilitate sponge microbiome comparisons across large spatial, temporal and environmental scales. Samples from marine sponges (n = 3569 specimens), seawater (n = 370), marine sediments (n = 65) and other environments (n = 29) were collected from different locations across the globe. This dataset incorporates at least 269 different sponge species, including several yet unidentified taxa. The V4 region of the 16S rRNA gene was amplified and sequenced from extracted DNA using standardised procedures. Raw sequences (total of 1.1 billion sequences) were processed and clustered with a) a standard protocol using QIIME closed-reference picking resulting in 39,543 Operational Taxonomic Units (OTU) at 97% sequence identity, b) a de novo protocol using Mothur resulting in 518,246 OTUs, and c) a new high-resolution Deblur protocol resulting in 83,908 unique bacterial sequences. Abundance tables, representative sequences, taxonomic classifications and metadata are provided. This dataset represents a comprehensive resource of sponge-associated microbial communities based on 16S rRNA gene sequences that can be used to address overarching hypotheses regarding host-associated prokaryotes, including host-specificity, convergent evolution, environmental drivers of microbiome structure and the sponge-associated rare biosphere.

Description: Phytoplankton cell or colony sizes range from 〈1 µm to several cm, i.e. 4–5 orders of magnitude in linear dimensions, which is roughly equivalent to the log-size span within terrestrial vegetation. It is commonplace to assume that smaller phytoplankton have an advantage in growth related traits while larger ones are more resistant to losses. However, the current state of literature calls for a more differentiated view. It is still controversial, whether smaller phytoplankton have higher maximal growth rates (µmax) or if there is a peak of µmax at intermediate size (102 µm3 cell volume). Smaller phytoplankton have an advantage in nutrient acquisition at low concentrations while larger phytoplankton have an advantage in utilizing nutrient pulses and exploiting vertical gradients. At equal density, larger phytoplankton experience bigger sinking losses. Small phytoplankton (〈5–10 µm) are more affected mostly from grazing by protists and tunicates, while larger phytoplankton are more affected by copepod and krill grazing. Size spectra within the most important higher taxa show some conspicuous differences between marine and lake phytoplankton, e.g. the absence of very large diatoms (〉105 µm3) in lake phytoplankton and the absence of large (〉103 µm3) green algae in marine plankton. Overall, size is one of the most important traits for the performance of phytoplankton, but it is overly simplistic to equate small size with metabolic advantages

Description: Adaptive radiation is thought to be responsible for the evolution of a great portion of the past and present diversity of life. Instances of adaptive radiation, characterized by the rapid emergence of an array of species as a consequence to their adaptation to distinct ecological niches, are important study systems in evolutionary biology. However, because of the rapid lineage formation in these groups, and occasional gene flow between the participating species, it is often difficult to reconstruct the phylogenetic history of species that underwent an adaptive radiation. In this study, we present a novel approach for species-tree estimation in rapidly diversifying lineages, where introgression is known to occur, and apply it to a multimarker data set containing up to 16 specimens per species for a set of 45 species of East African cichlid fishes (522 individuals in total), with a main focus on the cichlid species flock of Lake Tanganyika. We first identified, using age distributions of most recent common ancestors in individual gene trees, those lineages in our data set that show strong signatures of past introgression. This led us to formulate three hypotheses of introgression between different lineages of Tanganyika cichlids: the ancestor of Boulengerochromini (or of Boulengerochromini and Bathybatini) received genomic material from the derived H-lineage; the common ancestor of Cyprichromini and Perissodini experienced, in turn, introgression from Boulengerochromini and/or Bathybatini; and the Lake Tanganyika Haplochromini and closely related riverine lineages received genetic material from Cyphotilapiini. We then applied the multispecies coalescent model to estimate the species tree of Lake Tanganyika cichlids, but excluded the lineages involved in these introgression events, as the multispecies coalescent model does not incorporate introgression. This resulted in a robust species tree, in which the Lamprologini were placed as sister lineage to the H-lineage (including the Eretmodini), and we identify a series of rapid splitting events at the base of the H-lineage. Divergence ages estimated with the multispecies coalescent model were substantially younger than age estimates based on concatenation, and agree with the geological history of the Great Lakes of East Africa. Finally, we formally tested the three hypotheses of introgression using a likelihood framework, and find strong support for introgression between some of the cichlid tribes of Lake Tanganyika. [Adaptive radiation; Cichlidae; introgression; Lake Tanganyika; species network.]

Description: To the Editor: We congratulate Halsey et al. for their important discussion of the random nature of P values in an earlier issue of this journal1. P values from identical experiments can differ greatly in a way that is surprising to many2, 3, 4, 5, 6. The failure to appreciate this…

Description: Ocean acidification is an escalating environmental issue and associated changes in the ocean carbonate system have implications for many calcifying organisms. The present study followed the growth of Sepia officinalis from early-stage embryos, through hatching, to 7-week-old juveniles. Responses of cuttlefish to elevated pCO(2) (hypercapnia) were investigated to test the impacts of near-future and extreme ocean acidification conditions on growth, developmental time, oxygen consumption, and yolk utilization as proxies for individual fitness. We further examined gross morphological characteristics of the internal calcareous cuttlebone to determine whether embryonically secreted shell lamellae are impacted by environmental hypercapnia. Embryonic growth was reduced and hatching delayed under elevated pCO(2), both at environmentally relevant levels (0.14 kPa pCO(2) similar to predicted ocean conditions in 2100) and extreme conditions (0.40 kPa pCO(2)). Comparing various metrics from control and intermediate treatments generally showed no significant difference in experimental measurements. Yet, results from the high pCO(2) treatment showed significant changes compared with controls and revealed a consistent general trend across the three treatment levels. The proportion of animal mass contributed by the cuttlebone increased in both elevated pCO(2) treatments. Gross cuttlebone morphology was affected under such conditions and cuttlebones of hypercapnic individuals were proportionally shorter. Embryonic shell morphology was maintained consistently in all treatments, despite compounding hypercapnia in the perivitelline fluid; however, post-hatching, hypercapnic animals developed denser cuttlebone laminae in shorter cuttlebones. Juvenile cuttlefish in acidified environments thus experience lower growth and yet increased calcification of their internal shell. The results of this study support recent findings that early cuttlefish life stages are more vulnerable towards hypercapnia than juveniles and adults, which may have negative repercussions on the biological fitness of cuttlefish hatchlings in future oceans.

Description: Annual growth zones in cod otoliths from the eastern Baltic stock are less discrete than in other cod stocks leading to biased age reading, which recently led to a failure of age-based assessment in the eastern Baltic cod stock. In this study, we explored the applicability of minor and trace element patterns in cod otoliths for age determination. By first identifying elements of interest in a stock without ageing problems, western Baltic cod, we then tested their applicability on another stock without ageing problems, North Sea cod, and finally applied this knowledge to estimate age of eastern Baltic cod. In western Baltic cod, matching patterns with respect to occurrence of minima and maxima in both otolith opacity and element concentrations were found for Cu, Zn, and Rb, and inverse patterns with Mg and Mn. No match was found for Pb, Ba, and Sr. In the test stock, the North Sea cod, the same patterns in Cu, Zn, Rb, Mg, and Mn signals occurred. All eastern Baltic cod with low visual contrast between growth zones exhibited clearly defined synchronous cycles in Cu, Zn, Rb and Pb. Using a combined finite differencing method and structural break models approach, the statistical significance of the local profile minima were identified, based on which their age could be estimated. Despite extensive environmental differences between the three areas examined, the element concentrations of Cu, Zn, and Rb were strongly correlated in all individuals with similar correlations in all three areas, suggesting that the incorporation mechanisms are the same for these elements and independent of environmental concentrations.

Description: Global warming and ocean acidification are among the most important stressors for aquatic ecosystems in the future. To investigate their direct and indirect effects on a near-natural plankton community, a multiple-stressor approach is needed. Hence, we set up mesocosms in a full-factorial design to study the effects of both warming and high CO2 on a Baltic Sea autumn plankton community, concentrating on the impacts on microzooplankton (MZP). MZP abundance, biomass, and species composition were analysed over the course of the experiment. We observed that warming led to a reduced time-lag between the phytoplankton bloom and an MZP biomass maximum. MZP showed a significantly higher growth rate and an earlier biomass peak in the warm treatments while the biomass maximum was not affected. Increased pCO2 did not result in any significant effects on MZP biomass, growth rate, or species composition irrespective of the temperature, nor did we observe any significant interactions between CO2 and temperature. We attribute this to the high tolerance of this estuarine plankton community to fluctuations in pCO2, often resulting in CO2 concentrations higher than the predicted end-of-century concentration for open oceans. In contrast, warming can be expected to directly affect MZP and strengthen its coupling with phytoplankton by enhancing its grazing pressure.

Description: It is currently under debate whether organisms that regulate their acid–base status under environmental hypercapnia demand additional energy. This could impair animal fitness, but might be compensated for via increased ingestion rates when food is available. No data are yet available for dominant Calanus spp. from boreal and Arctic waters. To fill this gap, we incubated Calanus glacialis at 390, 1120, and 3000 µatm for 16 d with Thalassiosira weissflogii (diatom) as food source on-board RV Polarstern in Fram Strait in 2012. Every 4 d copepods were subsampled from all CO2 treatments and clearance and ingestion rates were determined. During the SOPRAN mesocosm experiment in Bergen, Norway, 2011, we weekly collected Calanus finmarchicus from mesocosms initially adjusted to 390 and 3000 µatm CO2 and measured grazing at low and high pCO2. In addition, copepods were deep frozen for body mass analyses. Elevated pCO2 did not directly affect grazing activities and body mass, suggesting that the copepods did not have additional energy demands for coping with acidification, neither during long-term exposure nor after immediate changes in pCO2. Shifts in seawater pH thus do not seem to challenge these copepod species.

Description: In a recent publication (Froese et al., ICES Journal of Marine Science; doi:10.1093/icesjms/fsv122), we presented a critique of the balanced harvesting (BH) approach to fishing. A short section dealt with the size-spectrum models used to justify BH, wherein we pointed out the lack of realism of these models, which mostly represented ecosystems as consisting of a single cannibalistic species. Andersen et al. (ICES Journal of Marine Science; doi:10.1093/icesjms/fsv211) commented on our paper and suggested that we criticized size-spectrum models in general and that we supposedly made several erroneous statements. We stress that we only referred to the size-spectrum models that we cited, and we respond to each supposedly erroneous statement. We still believe that the size-spectrum models used to justify BH were highly unrealistic and not suitable for evaluating real-world fishing strategies. We agree with Andersen et al. that BH is unlikely to be a useful guiding principle for ecosystem-based fisheries management, for many reasons. The use of unrealistic models is one of them.

Description: The quest for the ideal therapeutic target in chronic kidney disease (CKD) has been riddled with many obstacles stemming from the molecular complexity of the disease and its co-morbidities. Recent advances in omics technologies and the resulting amount of available data encompassing genomics, proteomics, peptidomics, transcriptomics and metabolomics has created an opportunity for integrating omics datasets to build a comprehensive and dynamic model of the molecular changes in CKD for the purpose of biomarker and drug discovery. This article reviews relevant concepts in omics data integration using systems biology, a mathematical modelling method that globally describes a biological system on the basis of its modules and the functional connections that govern their behaviour. The review describes key databases and bioinformatics tools, as well as the challenges and limitations of the current state of the art, along with practical application to CKD therapeutic target discovery. Moreover, it describes how systems biology and visualization tools can be used to generate clinically relevant molecular models with the capability to identify specific disease pathways, recognize key events in disease development and track disease progression.

Description: This paper applies nonlinear Bayesian inversion to marine controlled source electromagnetic (CSEM) data collected near two sites of the Integrated Ocean Drilling Program (IODP) Expedition 311 on the northern Cascadia Margin to investigate subseafloor resistivity structure related to gas hydrate deposits and cold vents. The Cascadia margin, off the west coast of Vancouver Island, Canada, has a large accretionary prism where sediments are under pressure due to convergent plate boundary tectonics. Gas hydrate deposits and cold vent structures have previously been investigated by various geophysical methods and seabed drilling. Here, we invert time-domain CSEM data collected at Sites U1328 and U1329 of IODP Expedition 311 using Bayesian methods to derive subsurface resistivity model parameters and uncertainties. The Bayesian information criterion is applied to determine the amount of structure (number of layers in a depth-dependent model) that can be resolved by the data. The parameter space is sampled with the Metropolis–Hastings algorithm in principal-component space, utilizing parallel tempering to ensure wider and efficient sampling and convergence. Nonlinear inversion allows analysis of uncertain acquisition parameters such as time delays between receiver and transmitter clocks as well as input electrical current amplitude. Marginalizing over these instrument parameters in the inversion accounts for their contribution to the geophysical model uncertainties. One-dimensional inversion of time-domain CSEM data collected at measurement sites along a survey line allows interpretation of the subsurface resistivity structure. The data sets can be generally explained by models with 1 to 3 layers. Inversion results at U1329, at the landward edge of the gas hydrate stability zone, indicate a sediment unconformity as well as potential cold vents which were previously unknown. The resistivities generally increase upslope due to sediment erosion along the slope. Inversion results at U1328 on the middle slope suggest several vent systems close to Bullseye vent in agreement with ongoing interdisciplinary observations.

Description: The TAMMAR segment of the Mid-Atlantic Ridge forms a classic propagating system centred about two degrees south of the Kane Fracture Zone. The segment is propagating to the south at a rate of 14 mm yr−1, 15 per cent faster than the half-spreading rate. Here, we use seismic refraction data across the propagating rift, sheared zone and failed rift to investigate the crustal structure of the system. Inversion of the seismic data agrees remarkably well with crustal thicknesses determined from gravity modelling. We show that the crust is thickened beneath the highly magmatic propagating rift, reaching a maximum thickness of almost 8 km along the seismic line and an inferred (from gravity) thickness of about 9 km at its centre. In contrast, the crust in the sheared zone is mostly 4.5–6.5 km thick, averaging over 1 km thinner than normal oceanic crust, and reaching a minimum thickness of only 3.5 km in its NW corner. Along the seismic line, it reaches a minimum thickness of under 5 km. The PmP reflection beneath the sheared zone and failed rift is very weak or absent, suggesting serpentinisation beneath the Moho, and thus effective transport of water through the sheared zone crust. We ascribe this increased porosity in the sheared zone to extensive fracturing and faulting during deformation. We show that a bookshelf-faulting kinematic model predicts significantly more crustal thinning than is observed, suggesting that an additional mechanism of deformation is required. We therefore propose that deformation is partitioned between bookshelf faulting and simple shear, with no more than 60 per cent taken up by bookshelf faulting.

Description: High-throughput techniques based on restriction site-associated DNA sequencing (RADseq) are enabling the low-cost discovery and genotyping of thousands of genetic markers for any species, including non-model organisms, which is revolutionizing ecological, evolutionary and conservation genetics. Technical differences among these methods lead to important considerations for all steps of genomics studies, from the specific scientific questions that can be addressed, and the costs of library preparation and sequencing, to the types of bias and error inherent in the resulting data. In this Review, we provide a comprehensive discussion of RADseq methods to aid researchers in choosing among the many different approaches and avoiding erroneous scientific conclusions from RADseq data, a problem that has plagued other genetic marker types in the past.

Description: Aquatic invertebrate communities are influenced by interactions between the abiotic and biotic environment at multiple spatial and temporal scales. Studies of mesozooplankton community patterns in relation to spatial and temporal scales are rare. In this study, we examined scale-specific variability of mesozooplankton in the shallow coastal Baltic Sea and related this variability to key environmental proxies. Seasonality defined the majority of variability in taxonomic composition and abundance patterns, as well as in aggregated parameters of zooplankton. However, these properties also varied spatially at a large, 100-km scale. The variability in all properties except taxonomic composition was negligible at the smaller spatial scale. Taxonomic richness increased until moderate levels of total abundance, whereas peak blooms were always characterized by higher dark diversity. Shannon diversity was unrelated to total abundance. Observed spatio-temporal patterns were strongly related to abiotic forcing and uncoupled from phytoplankton standing stock and primary production. Results show the importance of seasonality over spatial variability and abiotic factors over phytoplankton variability for sub-boreal brackish coastal mesozooplankton at the spatial scales studied. Information loss from spatial generalization can be larger for taxonomic occurrences and rare species than for species abundances and aggregated community parameters such as total abundance or taxonomic richness.

Description: Little is known about the effects that subducting an oceanic large igneous province (LIP) has on the petrogenesis of submarine arc volcanoes and their geochemical composition. The southern Kermadec arc represents a rare example where an LIP—the Hikurangi Plateau—is currently subducting and where its effect on mantle composition, element recycling and arc volcanism can be studied. We present mineral chemistry and whole-rock major and trace element, and Sr–Nd–Pb isotope data from samples recovered from the southern Kermadec arc volcanoes Rumble II East and Rumble II West, together with shipboard gravity and magnetic measurements. The Rumble II volcanoes (including a volcanic cone ∼10 km further west) form an ∼23 km long arc–backarc transect located ∼250 km north of New Zealand above the subducting Hikurangi Plateau. Although only a short distance apart, rocks from the two volcanoes have different mineral and whole-rock geochemical compositions. Lavas from Rumble II East are predominantly basaltic and contain primitive olivine phenocrysts (≤Fo91), high-Mg# clinopyroxene (≤96) and anorthitic plagioclase (≤An97). Geochemically these lavas are very diverse and cover a spectrum from low Th/Yb (〈0·15) at high Ba/Th (〉1014) to higher Th/Yb (〉0·15) at lower Ba/Th (〈844). This spectrum, together with 206Pb/204Pb and 143Nd/144Nd in the range of 18·74–18·83 and 0·51309–0·51298 respectively (at similar to slightly elevated 87Sr/86Sr), suggests a mantle wedge that has undergone previous melt extraction and significant fluid addition from the subducting Pacific Plate and that contains sediment and HIMU-type Hikurangi Plateau components. The geochemistry of the sediment–HIMU-type components is exemplified in an olivine pyroxenite (e.g. 206Pb/204Pb = 20·02; 87Sr/86Sr = 0·70516; 143Nd/144Nd = 0·5126). We propose that the olivine pyroxenite formed through melt or fluid–rock metasomatism and represents the first direct evidence of a near Moho arc mantle rock that shows the imprint from a subducting HIMU-type (Hikurangi) seamount. Conversely, lavas from Rumble II West and the cone ∼10 km to the west are generally more silica rich than Rumble II East lavas and mainly contain plagioclase with less ortho- and clinopyroxene + olivine phenocrysts. The low Ba/Th (〈470) and 206Pb/204Pb (〈18·74), a range of 143Nd/144Nd (0·51297–0·51307) and elevated Th/Yb (0·13–0·39) in these lavas can best be explained by minor sediment input into a less depleted mantle wedge. In addition, the geochemical composition of the Rumble II West lavas does not require involvement of a Hikurangi component, placing a spatial limit on Hikurangi material influencing regional melt generation beneath the backarc. Supported by a gravity model requiring two distinct magma chambers, the different geochemical compositions of Rumble II East and West lavas are inconsistent with a shared magma plumbing system. The different geochemical compositions of lavas from the two Rumble II volcanoes furthermore demonstrate that across-arc geochemical heterogeneities can occur within a few kilometres and may originate from both a geochemically heterogeneous mantle wedge and Moho transition layer, recording inherited geochemical heterogeneities beneath the volcanoes.

Description: Recent studies have analysed valuable compilations of data for the size-scaling of phytoplankton traits, but these cannot be employed directly in most large-scale modelling studies, which typically do not explicitly resolve the relevant trait values. Although some recent large-scale modelling studies resolve species composition and sorting within communities, most do not account for the observed flexible response of phytoplankton communities, such as the dynamic acclimation often observed in laboratory experiments. In order to derive a simple yet flexible model of phytoplankton growth that can be useful for a wide variety of ocean modelling applications, we combine two trade-offs, one for growth and the other for nutrient uptake, under the optimality assumption, i.e. that intracellular resources are dynamically allocated to maximize growth rate. This yields an explicit equation for growth as a function of nutrient concentration and daily averaged irradiance. We furthermore show how with this model effective Monod parameter values depend on both the underlying trait values and environmental conditions. We apply this new model to two contrasting time-series observation sites, including idealized simulations of size diversity. The flexible model responds differently compared with an inflexible control, suggesting that acclimation by individual species could impact models of plankton diversity.

Description: In the Baltic Sea, two genetically distinct cod populations occur, the eastern and the western Baltic cod. Since 2006, cod abundance has increased substantially in the Arkona Basin (SD 24), the potential mixing area between the two stocks management areas, presumably due to spill-over from the eastern stock. In this study, the spatio-temporal dynamics of stock mixing were analysed using shape analysis of archived otoliths. Further, the impact of eastern cod immigration on recruitment in the western Baltic Sea was investigated using hydrographic drift modelling. The percentage of eastern Baltic cod in the Arkona Basin increased from ca. 30% before 2005 to 〉80% in recent years. Geographic patterns in stock mixing with a pronounced east–west trend suggest that immigration occurs north of Bornholm, but propagates throughout the Arkona Basin. The immigration cannot be attributed to spawning migration, as no seasonal trend in stock mixing was observed. Based on environmental threshold levels for egg survival and time-series of hydrography data, the habitat suitable for successful spawning of eastern cod was estimated to range between 20 and 50% of the maximum possible habitat size, limited by primarily low salinity. Best conditions occurred irregularly in May–end June, interspersed with years where successful spawning was virtually impossible. Using a coupled hydrodynamic modelling and particle-tracking approach, the drift and survival of drifters representing eastern cod eggs was estimated. On average, 19% of the drifters in the Arkona Basin survive to the end of the yolk-sac stage, with mortality primarily after bottom contact due to low salinity. The general drift direction of the surviving larvae was towards the east. Therefore, it is the immigration of eastern cod, rather than larval transport, that contributes to cod recruitment in the western Baltic Sea.

Description: A recent optimality-based model for phytoplankton growth and diazotrophy was applied at two stations located in the oligotrophic western and the ultra-oligotrophic eastern subtropical North Atlantic. Contrary to the common view that diazotrophy is favoured by nitrogen (N) depletion relative to the Redfield equivalent of phosphorus (P), we find that optimality-based diazotrophy could explain N fixation in both regions in spite of relatively high N:P supply ratios. This is possible because the availability of an additional source of N for diazotrophs makes them strong competitors for P under oligotrophic conditions. The best reproduction of observations, especially of net primary production, is only achieved with preferential remineralization of P relative to N and atmospheric deposition. In line with observations, a higher rate of nitrogen fixation is predicted for the eastern site, owing to a larger niche for diazotrophs resulting from stronger oligotrophy and lower N:P supply ratios due to weaker atmospheric N deposition. Because the competitive advantage of diazotrophs under nutrient starvation diminishes with increasing supply N:P ratio, the predicted increase of atmospheric N deposition due to anthropogenic activity could negatively affect N2 fixation in the Atlantic Ocean.

Description: Remarkable technological advances have revealed ever more properties and behaviours of individual microorganisms, but the novel data generated by these techniques have not yet been fully exploited. In this Opinion article, we explain how individual-based models (IBMs) can be constructed based on the findings of such techniques and how they help to explore competitive and cooperative microbial interactions. Furthermore, we describe how IBMs have provided insights into self-organized spatial patterns from biofilms to the oceans of the world, phage–CRISPR dynamics and other emergent phenomena. Finally, we discuss how combining individual-based observations with IBMs can advance our understanding at both the individual and population levels, leading to the new approach of microbial individual-based ecology (μIBE).

Description: Bacteria use quorum sensing to orchestrate gene expression programmes that underlie collective behaviours. Quorum sensing relies on the production, release, detection and group-level response to extracellular signalling molecules, which are called autoinducers. Recent work has discovered new autoinducers in Gram-negative bacteria, shown how these molecules are recognized by cognate receptors, revealed new regulatory components that are embedded in canonical signalling circuits and identified novel regulatory network designs. In this Review we examine how, together, these features of quorum sensing signal–response systems combine to control collective behaviours in Gram-negative bacteria and we discuss the implications for host–microbial associations and antibacterial therapy.

Description: Mesoscale ocean eddies are ubiquitous coherent rotating structures of water with radial scales on the order of 100 kilometers. Eddies play a key role in the transport and mixing of momentum and tracers across the World Ocean. We present a global daily mesoscale ocean eddy dataset that contains ~45 million mesoscale features and 3.3 million eddy trajectories that persist at least two days as identified in the AVISO dataset over a period of 1993–2014. This dataset, along with the open-source eddy identification software, extract eddies with any parameters (minimum size, lifetime, etc.), to study global eddy properties and dynamics, and to empirically estimate the impact eddies have on mass or heat transport. Furthermore, our open-source software may be used to identify mesoscale features in model simulations and compare them to observed features. Finally, this dataset can be used to study the interaction between mesoscale ocean eddies and other components of the Earth System.

Description: Microorganisms produce a wealth of structurally diverse specialized metabolites with a remarkable range of biological activities and a wide variety of applications in medicine and agriculture, such as the treatment of infectious diseases and cancer, and the prevention of crop damage. Genomics has revealed that many microorganisms have far greater potential to produce specialized metabolites than was thought from classic bioactivity screens; however, realizing this potential has been hampered by the fact that many specialized metabolite biosynthetic gene clusters (BGCs) are not expressed in laboratory cultures. In this Review, we discuss the strategies that have been developed in bacteria and fungi to identify and induce the expression of such silent BGCs, and we briefly summarize methods for the isolation and structural characterization of their metabolic products.

Description: Cultivable Actinobacteria are the largest source of microbially derived bioactive molecules. The high demand for novel antibiotics highlights the need for exploring novel sources of these bacteria. Microbial symbioses with sessile macro-organisms, known to contain bioactive compounds likely of bacterial origin, represent an interesting and underexplored source of Actinobacteria. We studied the diversity and potential for bioactive-metabolite production of Actinobacteria associated with two marine lichens ( Lichina confinis and L. pygmaea ; from intertidal and subtidal zones) and one littoral lichen ( Roccella fuciformis ; from supratidal zone) from the Brittany coast (France), as well as the terrestrial lichen Collema auriforme (from a riparian zone, Austria). A total of 247 bacterial strains were isolated using two selective media. Isolates were identified and clustered into 101 OTUs (98% identity) including 51 actinobacterial OTUs. The actinobacterial families observed were: Brevibacteriaceae, Cellulomonadaceae, Gordoniaceae, Micrococcaceae, Mycobacteriaceae, Nocardioidaceae, Promicromonosporaceae, Pseudonocardiaceae, Sanguibacteraceae and Streptomycetaceae. Interestingly, the diversity was most influenced by the selective media rather than lichen species or the level of lichen thallus association. The potential for bioactive-metabolite biosynthesis of the isolates was confirmed by screening genes coding for polyketide synthases types I and II. These results show that littoral lichens are a source of diverse potentially bioactive Actinobacteria.

Description: Large amounts of the greenhouse gas methane are released from the seabed to the water column1, where it may be consumed by aerobic methanotrophic bacteria2. The size and activity of methanotrophic communities, which determine the amount of methane consumed in the water column, are thought to be mainly controlled by nutrient and redox dynamics3–7. Here, we report repeated measurements of methanotrophic activity and community size at methane seeps west of Svalbard, and relate them to physical water mass properties and modelled ocean currents. We show that cold bottom water, which contained a large number of aerobic methanotrophs, was displaced by warmer water with a considerably smaller methanotrophic community within days. Ocean current simulations using a global ocean/sea-ice model suggest that this water mass exchange is consistent with short-term variations in the meandering West Spitsbergen Current. We conclude that the shift from an offshore to a nearshore position of the current can rapidly and severely reduce methanotrophic activity in the water column. Strong fluctuating currents are common at many methane seep systems globally, and we suggest that they affect methane oxidation in the water column at other sites, too.

Description: Climate change is altering oceanic conditions in a complex manner, and the concurrent amendment of multiple properties will modify environmental stress for primary producers. So far, global modelling studies have focused largely on how alteration of individual properties will affect marine life. Here, we use global modelling simulations in conjunction with rotated factor analysis to express model projections in terms of regional trends in concomitant changes to biologically influential multi-stressors. Factor analysis demonstrates that regionally distinct patterns of complex oceanic change are evident globally. Preliminary regional assessments using published evidence of phytoplankton responses to complex change reveal a wide range of future responses to interactive multi-stressors with 〈20–300% shifts in phytoplankton physiological rates, and many unexplored potential interactions. In a future ocean, provinces will encounter different permutations of change that will probably alter the dominance of key phytoplankton groups and modify regional productivity, ecosystem structure and biogeochemistry. Consideration of regionally distinct multi-stressor patterns can help guide laboratory and field studies as well as the interpretation of interactive multi-stressors in global models.

Description: Reflection on the rapidly growing field of ocean acidification research highlights priorities for future research on the changing ocean.

Description: Thin oceanic crust is formed by decompression melting of the upper mantle at mid-ocean ridges, but the origin of the thick and buoyant continental crust is enigmatic. Juvenile continental crust may form from magmas erupted above intraoceanic subduction zones, where oceanic lithosphere subducts beneath other oceanic lithosphere. However, it is unclear why the subduction of dominantly basaltic oceanic crust would result in the formation of andesitic continental crust at the surface. Here we use geochemical and geophysical data to reconstruct the evolution of the Central American land bridge, which formed above an intra-oceanic subduction system over the past 70Myr. We find that the geochemical signature of erupted lavas evolved from basaltic to andesitic about 10Myr ago - coincident with the onset of subduction of more oceanic crust that originally formed above the Galápagos mantle plume. We also find that seismic P-waves travel through the crust at velocities intermediate between those typically observed for oceanic and continental crust. We develop a continentality index to quantitatively correlate geochemical composition with the average P-wave velocity of arc crust globally. We conclude that although the formation and evolution of continents may involve many processes, melting enriched oceanic crust within a subduction zone - a process probably more common in the Archaean - can produce juvenile continental crust.

Description: The stomachs of most vertebrates operate at an acidic pH of 2 generated by the gastric H+/K+-ATPase located in parietal cells. The acidic pH in stomachs of vertebrates is believed to aid digestion and to protect against environmental pathogens. Little attention has been placed on whether acidic gastric pH regulation is a vertebrate character or a deuterostome ancestral trait. Here, we report alkaline conditions up to pH 10.5 in the larval digestive systems of ambulacraria (echinoderm + hemichordate), the closest relative of the chordate. Microelectrode measurements in combination with specific inhibitors for acid-base transporters and ion pumps demonstrated that the gastric alkalization machinery in sea urchin larvae is mainly based on direct H+ secretion from the stomach lumen and involves a conserved set of ion pumps and transporters. Hemichordate larvae additionally utilized HCO 3- transport pathways to generate even more alkaline digestive conditions. Molecular analyses in combination with acidification experiments supported these findings and identified genes coding for ion pumps energizing gastric alkalization. Given that insect larval guts were also reported to be alkaline, our discovery raises the hypothesis that the bilaterian ancestor utilized alkaline digestive system while the vertebrate lineage has evolved a strategy to strongly acidify their stomachs.

Description: Symbiotic relationships between phytoplankton and N-2-fixing microorganisms play a crucial role in marine ecosystems. The abundant and widespread unicellular cyanobacteria group A (UCYN-A) has recently been found to live symbiotically with a haptophyte. Here, we investigated the effect of nitrogen (N), phosphorus (P), iron (Fe) and Saharan dust additions on nitrogen (N-2) fixation and primary production by the UCYN-A-haptophyte association in the subtropical eastern North Atlantic Ocean using nifH expression analysis and stable isotope incubations combined with single-cell measurements. N-2 fixation by UCYN-A was stimulated by the addition of Fe and Saharan dust, although this was not reflected in the nifH expression. CO2 fixation by the haptophyte was stimulated by the addition of ammonium nitrate as well as Fe and Saharan dust. Intriguingly, the single-cell analysis using nanometer scale secondary ion mass spectrometry indicates that the increased CO2 fixation by the haptophyte in treatments without added fixed N is likely an indirect result of the positive effect of Fe and/or P on UCYN-A N-2 fixation and the transfer of N-2-derived N to the haptophyte. Our results reveal a direct linkage between the marine carbon and nitrogen cycles that is fuelled by the atmospheric deposition of dust. The comparison of single-cell rates suggests a tight coupling of nitrogen and carbon transfer that stays balanced even under changing nutrient regimes. However, it appears that the transfer of carbon from the haptophyte to UCYN-A requires a transfer of nitrogen from UCYN-A. This tight coupling indicates an obligate symbiosis of this globally important diazotrophic association.

Description: Olivine major and trace element compositions from 12 basalts from the southern Payenia volcanic province in Argentina have been analyzed by electron microprobe and laser ablation inductively coupled plasma mass spectrometry. The olivines have high Fe/Mn and low Ca/Fe and many fall at the end of the global olivine array, indicating that they were formed from a pyroxene-rich source distinct from typical mantle peridotite. The olivines with the highest Fe/Mn have higher Zn/Fe, Zn and Co and lower Co/Fe than the olivines with lower Fe/Mn, also suggesting contributions from a pyroxene-rich source. Together with whole-rock radiogenic isotopes and elemental concentrations, the samples indicate mixing between two mantle sources: (1) a pyroxene-rich source with EM-1 ocean island basalt type trace element and isotope characteristics; (2) a peridotitic source with more radiogenic Pb that was metasomatized by subduction-zone fluids and/or melts. The increasing contributions from the pyroxene-rich source in the southern Payenia basalts are correlated with an increasing Fe-enrichment, which caused the olivines to have lower forsterite contents at a given Ni content. Al-in-olivine crystallization temperatures measured on olivine–spinel pairs are between 1155 and 1243°C and indicate that the magmas formed at normal upper mantle (asthenospheric) temperatures of ∼1350°C. The pyroxene-rich material is interpreted to have been brought up from the deeper parts of the upper mantle by vigorous asthenospheric upwelling caused by break-off of the Nazca slab south of Payenia during the Pliocene and roll-back of the subducting slab beneath Payenia. The pyroxene-rich mantle mixed with peridotitic metasomatized South Atlantic mantle in the mantle wedge beneath Payenia.

Description: Quasi-decadal variability in solar irradiance has been suggested to exert a substantial effect on Earth’s regional climate. In the North Atlantic sector, the 11-year solar signal has been proposed to project onto a pattern resembling the North Atlantic Oscillation (NAO), with a lag of a few years due to ocean-atmosphere interactions. The solar/NAO relationship is, however, highly misrepresented in climate model simulations with realistic observed forcings. In addition, its detection is particularly complicated since NAO quasi-decadal fluctuations can be intrinsically generated by the coupled ocean-atmosphere system. Here we compare two multi-decadal ocean-atmosphere chemistry-climate simulations with and without solar forcing variability. While the experiment including solar variability simulates a 1–2-year lagged solar/NAO relationship, comparison of both experiments suggests that the 11-year solar cycle synchronizes quasi-decadal NAO variability intrinsic to the model. The synchronization is consistent with the downward propagation of the solar signal from the stratosphere to the surface.

Description: Dissolved oxygen in sea water affects marine habitats and biogeochemical cycles. Oceanic zones with oxygen deficits represent 7% of the volume and 8% of the area of the oceans, and are thought to be expanding. One of the most pronounced lies in the region off Peru, where mesoscale activity in the form of fronts and eddies is strong. Here, we study the dynamics of the Peruvian oxygen minimum zone in a Lagrangian framework, using a coupled physical-biogeochemical numerical model and finite-size Lyapunov exponent fields, to evaluate the role of mesoscale activity. We find that, at depths between 380 and 600 m, mesoscale structures have two distinct roles. First, their mean positions and paths delimit and maintain the oxygen minimum zone boundaries. Second, their high-frequency fluctuations inject oxygen across the oxygen minimum zone boundaries and eddy fluxes are one order of magnitude higher than mean oxygen fluxes. We conclude that these eddy fluxes contribute to the ventilation of the Peruvian oxygen minimum zone.

Description: The symbiotic association of corals and unicellular algae of the genus Symbiodinium in the southern Persian/Arabian Gulf (PAG) display an exceptional heat tolerance, enduring summer peak temperatures of up to 36 °C. As yet, it is not clear whether this resilience is related to the presence of specific symbiont types that are exclusively found in this region. Therefore, we used molecular markers to identify the symbiotic algae of three Porites species along 〉1000 km of coastline in the PAG and the Gulf of Oman and found that a recently described species, Symbiodinium thermophilum, is integral to coral survival in the southern PAG, the world’s hottest sea. Despite the geographic isolation of the PAG, we discovered that representatives of the S. thermophilum group can also be found in the adjacent Gulf of Oman providing a potential source of thermotolerant symbionts that might facilitate the adaptation of Indian Ocean populations to the higher water temperatures expected for the future. However, corals from the PAG associated with S. thermophilum show strong local adaptation not only to high temperatures but also to the exceptionally high salinity of their habitat. We show that their superior heat tolerance can be lost when these corals are exposed to reduced salinity levels common for oceanic environments elsewhere. Consequently, the salinity prevailing in most reefs outside the PAG might represent a distribution barrier for extreme temperature-tolerant coral/Symbiodinium associations from the PAG.

Description: The comb jelly Mnemiopsis leidyi is considered to be a successful invasive species, partly due to its high reproduction potential. However, due to the absence of direct carbon measurements of eggs, specific reproduction rates remain uncertain. We show that egg carbon is 0.22 ± 0.02 µg C and up to 21 times higher than previously extrapolated. With maximum rates of 11 232 eggs ind−1 day−1, largest animals in northern Europe invest ∼10% day−1 of their body carbon into reproduction. The comb jelly M. leidyi has received wide public and scientific attention during the last decades due to its commonly observed formation of bloom abundances in native and invaded areas (e.g. Costello et al., 2012; Riisgård et al., 2012). One of the traits suggested to be responsible for M. leidyi's invasion success is its high fecundity. At the northern end of its distribution range in native areas, M. leidyi has been shown to produce up to 9380 and 14 233 eggs ind−1 day−1 (Kremer, 1976a; Graham et al., 2009), with similar rates of 9910 eggs ind−1 day−1 for the native southern population in Biscayne Bay, FL, USA (Baker and Reeve, 1974). Within invaded European waters, rates of up to 3000 and 12 000 eggs ind−1 day−1 have been recorded for northern and southern populations, respectively (Zaika and Revkov, 1994; Javidpour et al., 2009). Since M. leidyi is a simultaneous hermaphrodite and fertilized eggs are produced on a daily basis during favorable conditions (Jaspers, 2012), M. leidyi can circumvent the Allee effect and efficiently seed new populations even from few founding individuals. Although …

Description: New marine geophysical data recorded across the Tonga-Kermadec subduction zone are used to image deformation and seismic velocity structures of the forearc and Pacific Plate where the Louisville Ridge seamount chain subducts. Due to the obliquity of the Louisville Ridge to the trench and the fast 128 mm yr−1 south–southwest migration of the ridge-trench collision zone, post-, current and pre-seamount subduction deformation can be investigated between 23°S and 28°S. We combine our interpretations from the collision zone with previous results from the post- and pre-collision zones to define the along-arc variation in deformation due to seamount subduction. In the pre-collision zone the lower-trench slope is steep, the mid-trench slope has ∼3-km-thick stratified sediments and gravitational collapse of the trench slope is associated with basal erosion by subducting horst and graben structures on the Pacific Plate. This collapse indicates that tectonic erosion is a normal process affecting this generally sediment starved subduction system. In the collision zone the trench-slope decreases compared to the north and south, and rotation of the forearc is manifest as a steep plate boundary fault and arcward dipping sediment in a 12-km-wide, ∼2-km-deep mid-slope basin. A ∼3 km step increase in depth of the middle and lower crustal isovelocity contours below the basin indicates the extent of crustal deformation on the trench slope. At the leading edge of the overriding plate, upper crustal P-wave velocities are ∼4.0 km s−1 and indicate the trench fill material is of seamount origin. Osbourn Seamount on the outer rise has extensional faulting on its western slope and mass wasting of the seamount provides the low Vp material to the trench. In the post-collision zone to the north, the trench slope is smooth, the trench is deep, and the crystalline crust thins at the leading edge of the overriding plate where Vp is low, ∼5.5 km s−1. These characteristics are attributed to a greater degree of extensional collapse of the forearc in the wake of seamount subduction. The northern end of a seismic gap lies at the transition from the smooth lower-trench slope of the post-collision zone, to the block faulted and elevated lower-trench slope in the collision zone, suggesting a causative link between the collapse of the forearc and seismogenesis. Along the forearc, the transient effects of a north-to-south progression of ridge subduction are preserved in the geomorphology, whereas longer-term effects may be recorded in the ∼80 km offset in trench strike at the collision zone itself.

Description: Global warming was reported to cause growth reductions in tropical shallow water corals in both, cooler and warmer, regions of the coral species range. This suggests regional adaptation with less heat-tolerant populations in cooler and more thermo-tolerant populations in warmer regions. Here, we investigated seasonal changes in the in situ metabolic performance of the widely distributed hermatypic coral Pocillopora verrucosa along 12 degrees latitudes featuring a steep temperature gradient between the northern (28.5 degrees N, 21-27 degrees C) and southern (16.5 degrees N, 28-33 degrees C) reaches of the Red Sea. Surprisingly, we found little indication for regional adaptation, but strong indications for high phenotypic plasticity: Calcification rates in two seasons (winter, summer) were found to be highest at 28-29 degrees C throughout all populations independent of their geographic location. Mucus release increased with temperature and nutrient supply, both being highest in the south. Genetic characterization of the coral host revealed low inter-regional variation and differences in the Symbiodinium clade composition only at the most northern and most southern region. This suggests variable acclimatization potential to ocean warming of coral populations across the Red Sea: high acclimatization potential in northern populations, but limited ability to cope with ocean warming in southern populations already existing at the upper thermal margin for corals

Description: The Murray Ridge/Dalrymple Trough system forms the boundary between the Indian and Arabian plates in the northern Arabian Sea. Geodetic constraints from the surrounding con- tinents suggest that this plate boundary is undergoing oblique extension at a rate of a few millimetres per year. We present wide-angle seismic data that constrains the composition of the Ridge and of adjacent lithosphere beneath the Indus Fan. We infer that Murray Ridge, like the adjacent Dalrymple Trough, is underlain by continental crust, while a thin crustal section beneath the Indus Fan represents thinned continental crust or exhumed serpentinized mantle that forms part of a magma-poor rifted margin. Changes in crustal structure across the Murray Ridge and Dalrymple Trough can explain short-wavelength gravity anomalies, but a long-wavelength anomaly must be attributed to deeper density contrasts that may result from a large age contrast across the plate boundary. The origin of this fragment of continental crust remains enigmatic, but the presence of basement fabrics to the south that are roughly parallel to Murray Ridge suggests that it separated from the India/Seychelles/Madagascar block by extension during early breakup of Gondwana.

Description: Over the past 60 years, both average daily precipitation intensity and extreme precipitation have increased in many regions1, 2, 3. Part of these changes, or even individual events4, 5, have been attributed to anthropogenic warming6, 7. Over the Black Sea and Mediterranean region, the potential for extreme summertime convective precipitation has grown8 alongside substantial sea surface temperature increase. A particularly devastating convective event experienced in that region was the July 2012 precipitation extreme near the Black Sea town of Krymsk9. Here we study the effect of sea surface temperature (SST) increase on convective extremes within the region, taking the Krymsk event as a showcase example. We carry out ensemble sensitivity simulations with a convection-permitting atmospheric model and show the crucial role of SST increase in the extremeness of the event. The enhancement of lower tropospheric instability due to the current warmer Black Sea allows deep convection to be triggered, increasing simulated precipitation by more than 300% relative to simulations with SSTs characteristic of the early 1980s. A highly nonlinear precipitation response to incremental SST increase suggests that the Black Sea has exceeded a regional threshold for the intensification of convective extremes. The physical mechanism we identify indicates that Black Sea and Mediterranean coastal regions may face abrupt amplifications of convective precipitation under continued SST increase, and illustrates the limitations of thermodynamical bounds for estimating the temperature scaling of convective extremes.

Description: To date, the parameters that determine the rupture area of great subduction zone earthquakes remain contentious. On 1 April 2014, the Mw 8.1 Iquique earthquake ruptured a portion of the well-recognized northern Chile seismic gap but left large highly coupled areas un-ruptured. Marine seismic reflection and swath bathymetric data indicate that structural variations in the subducting Nazca Plate control regional-scale plate-coupling variations, and the limited extent of the 2014 earthquake. Several under-thrusting seamounts correlate to the southward and up-dip arrest of seismic rupture during the 2014 Iquique earthquake, thus supporting a causal link. By fracturing of the overriding plate, the subducting seamounts are likely further responsible for reduced plate-coupling in the shallow subduction zone and in a lowly coupled region around 20.5°S. Our data support that structural variations in the lower plate influence coupling and seismic rupture offshore Northern Chile, whereas the structure of the upper plate plays a minor role.

Description: Oceanic dissolved organic carbon (DOC) is an important carbon pool, similar in magnitude to atmospheric CO2, but the fate of its oldest forms is not well understood. Hot hydrothermal circulation may facilitate the degradation of otherwise un-reactive dissolved organic matter, playing an important role in the long-term global carbon cycle. The oldest, most recalcitrant forms of DOC, which make up most of oceanic DOC, can be recovered by solid-phase extraction. Here we present measurements of solid-phase extractable DOC from samples collected between 2009 and 2013 at seven vent sites in the Atlantic, Pacific and Southern oceans, along with magnesium concentrations, a conservative tracer of water circulation through hydrothermal systems. We find that magnesium and solid-phase extractable DOC concentrations are correlated, suggesting that solid-phase extractable DOC is almost entirely lost from solution through mineralization or deposition during circulation through hydrothermal vents with fluid temperatures of 212-401°C. In laboratory experiments, where we heated samples to 380°C for four days, we found a similar removal efficiency. We conclude that thermal degradation alone can account for the loss of solid-phase extractable DOC in natural hydrothermal systems, and that its maximum lifetime is constrained by the timescale of hydrothermal cycling, at about 40 million years

Description: Cold-water corals (CWC) are widely distributed around the world forming extensive reefs at par with tropical coral reefs. They are hotspots of biodiversity and organic matter processing in the world’s deep oceans. Living in the dark they lack photosynthetic symbionts and are therefore considered to depend entirely on the limited flux of organic resources from the surface ocean. While symbiotic relations in tropical corals are known to be key to their survival in oligotrophic conditions, the full metabolic capacity of CWC has yet to be revealed. Here we report isotope tracer evidence for efficient nitrogen recycling, including nitrogen assimilation, regeneration, nitrification and denitrification. Moreover, we also discovered chemoautotrophy and nitrogen fixation in CWC and transfer of fixed nitrogen and inorganic carbon into bulk coral tissue and tissue compounds (fatty acids and amino acids). This unrecognized yet versatile metabolic machinery of CWC conserves precious limiting resources and provides access to new nitrogen and organic carbon resources that may be essential for CWC to survive in the resource-depleted dark ocean.

Description: The interoceanic transfer of seawater between the Indian Ocean and the Atlantic, ‘Agulhas leakage’, forms a choke point for the overturning circulation in the global ocean. Here, by combining output from a series of high-resolution ocean and climate models with in situ and satellite observations, we construct a time series of Agulhas leakage for the period 1870–2014. The time series demonstrates the impact of Southern Hemisphere westerlies on decadal timescales. Agulhas leakage shows a correlation with the Atlantic Multi-decadal Oscillation on multi-decadal timescales; the former leading by 15 years. This is relevant for climate in the North Atlantic

Description: Sea-level rise1 is one of the most pressing aspects of anthropogenic global warming with far-reaching consequences for coastal societies. However, sea-level rise did2, 3, 4, 5, 6, 7 and will strongly vary from coast to coast8, 9, 10. Here we investigate the long-term internal variability effects on centennial projections of dynamic sea level (DSL), the local departure from the globally averaged sea level. A large ensemble of global warming integrations has been conducted with a climate model, where each realization was forced by identical CO2 increase but started from different atmospheric and oceanic initial conditions. In large parts of the mid- and high latitudes, the ensemble spread of the projected centennial DSL trends is of the same order of magnitude as the globally averaged steric sea-level rise, suggesting that internal variability cannot be ignored when assessing twenty-first-century DSL trends. The ensemble spread is considerably reduced in the mid- to high latitudes when only the atmospheric initial conditions differ while keeping the oceanic initial state identical; indicating that centennial DSL projections are strongly dependent on ocean initial conditions.

Description: Aquatic ecosystems experience large natural variation in elemental composition of carbon (C), nitrogen (N) and phosphorus (P), which is further enhanced by human activities. Primary producers typically reflect the nutrient ratios of their resource, whose stoichiometric composition can vary widely in conformity to environmental conditions. In contrast, C to nutrient ratios in consumers are largely constrained within a narrow range, termed homeostasis. In comparison to crustacean zooplankton, less is known about the ability of protozoan grazers and rotifer species to maintain stoichiometric balance. In this study, we used laboratory experiments with a primary producer (Nannochloropsis sp.), three different species of protozoan grazers and one mesozooplankton species: two heterotrophic dinoflagellates (Gyrodinium dominans and Oxyrrhis marina), a ciliate (Euplotes sp.) and a rotifer (Brachionus plicatilis) to test the stoichiometric response to five nutrient treatments. We showed that the dependency of zooplankton C:N:P ratios on C: nutrient ratios of their food source varies among species. Similar to the photoautotroph, the two heterotrophic dinoflagellates weakly regulated their internal stoichiometry. In contrast, the strength of stoichiometric regulation increased to strict homeostasis in both the ciliate and the rotifer, similar to crustacean zooplankton. Our study further shows that ciliate and rotifer growth can be constrained by imbalanced resource supply. It also indicates that these key primary consumers have the potential to trophically upgrade poor stoichiometric autotrophic food quality for higher trophic levels.

Description: Restructure data-gathering and evaluation networks to address climate change, energy, food, health and water provision, say Yonglong Lu and colleagues.

Description: Global mean surface warming has stalled since the end of the twentieth century1, 2, but the net radiation imbalance at the top of the atmosphere continues to suggest an increasingly warming planet. This apparent contradiction has been reconciled by an anomalous heat flux into the ocean3, 4, 5, 6, 7, 8, induced by a shift towards a La Niña-like state with cold sea surface temperatures in the eastern tropical Pacific over the past decade or so. A significant portion of the heat missing from the atmosphere is therefore expected to be stored in the Pacific Ocean. However, in situ hydrographic records indicate that Pacific Ocean heat content has been decreasing9. Here, we analyse observations along with simulations from a global ocean–sea ice model to track the pathway of heat. We find that the enhanced heat uptake by the Pacific Ocean has been compensated by an increased heat transport from the Pacific Ocean to the Indian Ocean, carried by the Indonesian throughflow. As a result, Indian Ocean heat content has increased abruptly, which accounts for more than 70% of the global ocean heat gain in the upper 700 m during the past decade. We conclude that the Indian Ocean has become increasingly important in modulating global climate variability.

Description: An index of water-circulation strength in the North Atlantic Ocean has been derived from sea-level measurements. This provides fresh evidence of the ocean's leading role in multidecadal climate variability. See Letter p.508

Description: About 60 years ago, the critical depth hypothesis was proposed to describe the occurrence of spring phytoplankton blooms and emphasized the role of stratification for the timing of onset. Since then, several alternative hypotheses appeared focusing on the role of grazing and mixing processes such as turbulent convection or wind activity. Surprisingly, the role of community composition—and thus the distribution of phytoplankton traits—for bloom formation has not been addressed. Here, we discuss how trait variability between competing species might influence phytoplankton growth during the onset of the spring bloom. We hypothesize that the bloom will only occur if there are species with a combination of traits fitting to the environmental conditions at the respective location and time. The basic traits for formation of the typical spring bloom are high growth rates and photoadaptation to low light conditions, but other traits such as nutrient kinetics and grazing resistance might also be important. We present concise ideas on how to test our theoretical considerations experimentally. Furthermore, we suggest that future models of phytoplankton blooms should include both water column dynamics and variability of phytoplankton traits to make realistic projections instead of treating the phytoplankton bloom as an aggregate community phenomenon.

Description: Increasingly, spatial geochemical zonation, present as geographically distinct, subparallel trends, is observed along hotspot tracks, such as Hawaii and the Galapagos. The origin of this zonation is currently unclear. Recently zonation was found along the last B70 Myr of the Tristan-Gough hotspot track. Here we present new Sr–Nd–Pb–Hf isotope data from the older parts of this hotspot track (Walvis Ridge and Rio Grande Rise) and re-evaluate published data from the Etendeka and Parana flood basalts erupted at the initiation of the hotspot track. We show that only the enriched Gough, but not the less-enriched Tristan, component is present in the earlier (70–132 Ma) history of the hotspot. Here we present a model that can explain the temporal evolution and origin of plume zonation for both the Tristan-Gough and Hawaiian hotspots, two end member types of zoned plumes, through processes taking place in the plume sources at the base of the lower mantle.