ALBERT

Description: Male secondary sexual traits are targets of inter- and/or intrasexual selection, but can vary due to a correlation with life-history traits or as by-product of adaptation to distinct environments. Trade-offs contributing to this variation may comprise conspicuousness towards conspecifics versus inconspicuousness towards predators, or between allocating resources into coloration versus the immune system. Here, we examine variation in expression of a carotenoid-based visual signal, anal-fin egg-spots, along a replicate environmental gradient in the haplochromine cichlid fish Astatotilapia burtoni. We quantified egg-spot number, area, and coloration; applied visual models to estimate the trait's conspicuousness when perceived against the surrounding tissue under natural conditions; and used the lymphocyte ratio as a measure for immune activity. We find that (i) males possess larger and more conspicuous egg-spots than females, which is likely explained by their function in sexual selection; (ii) riverine fish generally feature fewer but larger and/or more intensively colored egg-spots, which is probably to maintain signal efficiency in intraspecific interactions in long-wavelength shifted riverine light conditions; and (iii) egg-spot number and relative area correlate with immune defense, suggesting a trade-off in the allocation of carotenoids. Taken together, haplochromine egg-spots feature the potential to adapt to the respective underwater light environment, and are traded-off with investment into the immune system

Description: The ocean's potential to export carbon to depth partly depends on the fraction of primary production (PP) sinking out of the euphotic zone (i.e., the e-ratio). Measurements of PP and export flux are often performed simultaneously in the field, although there is a temporal delay between those parameters. Thus, resulting e-ratio estimates often incorrectly assume an instantaneous downward export of PP to export flux. Evaluating results from four mesocosm studies, we find that peaks in organic matter sedimentation lag chlorophyll a peaks by 2 to 15 days. We discuss the implications of these time lags (TLs) for current e-ratio estimates and evaluate potential controls of TL. Our analysis reveals a strong correlation between TL and the duration of chlorophyll a buildup, indicating a dependency of TL on plankton food web dynamics. This study is one step further toward time-corrected e-ratio estimates

Description: Stable isotope compositions can potentially be used to trace atmospheric Cd inputs to the surface ocean and anthropogenic Cd emissions to the atmosphere. Both of these applications may provide valuable insights into the effects of anthropogenic activities on the cycling of Cd in the environment. However, a lack of constraints for the Cd isotope compositions of atmospheric aerosols is currently hindering such studies. Here, we present stable Cd isotope data for aerosols collected over the Tropical Atlantic Ocean. The samples feature variable proportions of mineral dust-derived and anthropogenic Cd, yet exhibit similar isotope compositions, thus negating the distinction of these Cd sources using isotopic signatures in this region. Isotopic variability between these two atmospheric Cd sources may be identified in other areas, and thus warrants further investigation. Regardless, these data provide important initial constraints on the isotope composition of atmospheric Cd inputs to the ocean.

Description: Aim: The lives of juvenile leatherback turtles are amongst the most enigmatic of all marine mega-vertebrates. For these cryptic organisms, ocean models provide important insights into their dispersion from natal sites. Here, corroborated by fisheries bycatch data, we simulate spatio-temporal variation in hatchling dispersion patterns over five decades from the World's largest leatherback turtle nesting region. Location: Equatorial Central West Africa (3.5°N to −6°S) spanning the Gulf of Guinea in the North, Gabon and the Republic/Democratic Republic of the Congo in the South. Results: Due to dynamic oceanic conditions at these equatorial latitudes, dispersion scenarios differed significantly: (1) along the north to south gradient of the study region, (2) seasonally and (3) between years. From rookeries to the north of the equator, simulated hatchling retention rates within the Gulf of Guinea were very high (〉99%) after 6 months of drift, whilst south of the equator, retention rates were as low as c. 6% with the majority of simulated hatchlings dispersing west into the South Atlantic Ocean with the South Equatorial Current. Seasonal dispersion variability was driven by wind changes arising from the yearly north/southward migration of the intertropical convergence zone resulting in the increasing westerly dispersion of hatchlings throughout the hatching season. Annual variability in wind stress drove a long-term trend for decreased retention within the Gulf of Guinea and increased westerly dispersion into habitats in the South Atlantic Ocean. Main conclusions: Shifts in dispersion habitats arising from spatio-temporal oceanic variability expose hatchlings to different environments and threats that will influence important life history attributes such as juvenile growth/survival rates; anticipated to impact the population dynamics and size/age structure of populations into adulthood. The impacts of local and dynamic oceanic conditions thus require careful considerations, such as subregional management, when managing marine populations of conservation concern.

Description: The waning stage(s) of the Tethyan ocean(s) in the Balkans are not well understood. Controversy centres on the origin and life-span of the Cretaceous Sava Zone, which is allegedly a remnant of the last oceanic domain in the Balkan Peninsula, defining the youngest suture between Eurasia- and Adria-derived plates. In order to investigate to what extent late-Cretaceous volcanism within the Sava zone is consistent with this model, we present new age data together with trace-element and Sr–Nd–Pb isotope data for the Klepa basaltic lavas from the central Balkan Peninsula. Our new geochemical data show marked differences between the Cretaceous Klepa basalts (Sava Zone) and the rocks of other volcanic sequences from the Jurassic ophiolites of the Balkans. The Klepa basalts mostly have Sr–Nd–Pb isotopic and trace-element signatures that resemble enriched within-plate basalts, substantially different from Jurassic ophiolite basalts with MORB, BAB and IAV affinities. Trace-element modelling of the Klepa rocks indicates 2–20% polybaric melting of a relatively homogeneously metasomatised mantle source that ranges in composition from garnet lherzolite to ilmenite+apatite bearing spinel–amphibole lherzolite. Thus, the residual mineralogy is characteristic of a continental rather than oceanic lithospheric mantle source, suggesting an intracontinental within-plate origin for the Klepa basalts. Two alternative geodynamic models are internally consistent with our new findings: i) if the Sava Zone represents remnants of the youngest Neotethyan Ocean, magmatism along this zone would be situated within the forearc region and triggered by ridge subduction; ii) if the Sava Zone delimits a diffuse tectonic boundary between Adria and Europe, which had already collided in the Late Jurassic, the Klepa basalts together with a number of other magmatic centres represent volcanism related to transtensional tectonics.

Description: The upstream sources and pathways of the Denmark Strait Overflow Water and their variability have been investigated using a high-resolution model hindcast. This global simulation covers the period from 1948 to 2009 and uses a fine model mesh (1/20°) to resolve mesoscale features and the complex current structure north of Iceland explicitly. The three sources of the Denmark Strait Overflow, the shelfbreak East Greenland Current (EGC), the separated EGC, and the North Icelandic Jet, have been analyzed using Eulerian and Lagrangian diagnostics. The shelfbreak EGC contributes the largest fraction in terms of volume and freshwater transport to the Denmark Strait Overflow and is the main driver of the overflow variability. The North Icelandic Jet contributes the densest water to the Denmark Strait Overflow and shows only small temporal transport variations. During summer, the net volume and freshwater transports to the south are reduced. On interannual time scales, these transports are highly correlated with the large-scale wind stress curl around Iceland and, to some extent, influenced by the North Atlantic Oscillation, with enhanced southward transports during positive phases. The Lagrangian trajectories support the existence of a hypothesized overturning loop along the shelfbreak north of Iceland, where water carried by the North Icelandic Irminger Current is transformed and feeds the North Icelandic Jet. Monitoring these two currents and the region north of the Iceland shelfbreak could provide the potential to track long-term changes in the Denmark Strait Overflow and thus also the AMOC.

Description: The North Atlantic Current (NAC) is subject to variability on multiannual to decadal time scales, influencing the transport of volume, heat, and freshwater from the subtropical to the eastern subpolar North Atlantic (NA). Current observational time series are either too short or too episodic to study the processes involved. Here we compare the observed continuous NAC transport time series at the western flank of the Mid-Atlantic Ridge (MAR) and repeat hydrographic measurements at the OVIDE line in the eastern Atlantic with the NAC transport and circulation in the high-resolution (1/20°) ocean model configuration VIKING20 (1960–2008). The modeled baroclinic NAC transport relative to 3400 m (24.5 ± 7.1 Sv) at the MAR is only slightly lower than the observed baroclinic mean of 27.4 ± 4.7 Sv from 1993 to 2008, and extends further north by about 0.5°. In the eastern Atlantic, the western NAC (WNAC) carries the bulk of the transport in the model, while transport estimates based on hydrographic measurements from five repeated sections point to a preference for the eastern NAC (ENAC). The model is able to simulate the main features of the subpolar NA, providing confidence to use the model output to analyze the influence of the North Atlantic Oscillation (NAO). Model based velocity composites reveal an enhanced NAC transport across the MAR of up to 6.7 Sv during positive NAO phases. Most of that signal (5.4 Sv) is added to the ENAC transport, while the transport of the WNAC was independent of the NAO.

Description: While the Earth's surface has considerably warmed over the past two decades, the tropical Pacific has featured a cooling of sea surface temperatures in its eastern and central parts, which went along with an unprecedented strengthening of the equatorial trade winds, the surface component of the Pacific Walker Circulation (PWC). Previous studies show that this decadal trend in the trade winds is generally beyond the range of decadal trends simulated by climate models when forced by historical radiative forcing. There is still a debate on the origin of and the potential role that internal variability may have played in the recent decadal surface wind trend. Using a number of long control (unforced) integrations of global climate models and several observational data sets, we address the question as to whether the recent decadal to multidecadal trends are robustly classified as an unusual event or the persistent response to external forcing. The observed trends in the tropical Pacific surface climate are still within the range of the long-term internal variability spanned by the models but represent an extreme realization of this variability. Thus, the recent observed decadal trends in the tropical Pacific, though highly unusual, could be of natural origin. We note that the long-term trends in the selected PWC indices exhibit a large observational uncertainty, even hindering definitive statements about the sign of the trends.

Description: We examine the mean pathways, transit timescales, and transformation of waters flowing from the Pacific and the marginal seas through the Indian Ocean (IO) on their way toward the South Atlantic within a high-resolution ocean/sea-ice model. The model fields are analyzed from a Lagrangian perspective where water volumes are tracked as they enter the IO. The IO contributes 12.6 Sv to Agulhas leakage, which within the model is 14.1 ± 2.2 Sv, the rest originates from the South Atlantic. The Indonesian Through-flow constitutes about half of the IO contribution, is surface bound, cools and salinificates as it leaves the basin within 10–30 years. Waters entering the IO south of Australia are at intermediate depths and maintain their temperature-salinity properties as they exit the basin within 15–35 years. Of these waters, the contribution from Tasman leakage is 1.4 Sv. The rest stem from recirculation from the frontal regions of the Southern Ocean. The marginal seas export 1.0 Sv into the Atlantic within 15–40 years, and the waters cool and freshen on-route. However, the model's simulation of waters from the Gulfs of Aden and Oman are too light and hence overly influenced by upper ocean circulations. In the Cape Basin, Agulhas leakage is well mixed. On-route, temperature-salinity transformations occur predominantly in the Arabian Sea and within the greater Agulhas Current region. Overall, the IO exports at least 7.9 Sv from the Pacific to the Atlantic, thereby quantifying the strength of the upper cell of the global conveyor belt.

Description: The potential of mining seafloor massive sulfide deposits for metals such as Cu, Zn, and Au is currently debated. One key challenge is to predict where the largest deposits worth mining might form, which in turn requires understanding the pattern of subseafloor hydrothermal mass and energy transport. Numerical models of heat and fluid flow are applied to illustrate the important role of fault zone properties (permeability and width) in controlling mass accumulation at hydrothermal vents at slow spreading ridges. We combine modeled mass-flow rates, vent temperatures, and vent field dimensions with the known fluid chemistry at the fault-controlled Logatchev 1 hydrothermal field of the Mid-Atlantic Ridge. We predict that the 135 kilotons of SMS at this site (estimated by other studies) can have accumulated with a minimum depositional efficiency of 5% in the known duration of hydrothermal venting (58,200 year age of the deposit). In general, the most productive faults must provide an efficient fluid pathway while at the same time limit cooling due to mixing with entrained cold seawater. This balance is best met by faults that are just wide and permeable enough to control a hydrothermal plume rising through the oceanic crust. Model runs with increased basal heat input, mimicking a heat flow contribution from along-axis, lead to higher mass fluxes and vent temperatures, capable of significantly higher SMS accumulation rates. Nonsteady state conditions, such as the influence of a cooling magmatic intrusion beneath the fault zone, also can temporarily increase the mass flux while sustaining high vent temperatures.

Description: Diazotrophic dinitrogen (N2) fixation contributes ~76% to "new" nitrogen inputs to the sunlit open ocean, but environmental factors determining N2 fixation rates are not well constrained. Excess phosphate (phosphate-nitrate/16 〉 0) and iron availability control N2 fixation rates in the eastern tropical North Atlantic (ETNA), but it remains an open question how excess phosphate is generated within or supplied to the phosphate-depleted sunlit layer. Our observations in the ETNA region (8°N-15°N, 19°W-23°W) suggest that Prochlorococcus and Synechococcus, the two ubiquitous non-diazotrophic cyanobacteria with cellular N:P ratios higher than the Redfield ratio, create an environment of excess phosphate, which cannot be explained by diapycnal mixing, atmospheric, and riverine inputs. Thus, our results unveil a new biogeochemical niche construction mechanism by non-diazotrophic cyanobacteria for their diazotrophic phylum group members (N2 fixers). Our observations may help to understand the prevalence of diazotrophy in low-phosphate, oligotrophic regions.

Description: Oceanographic observations from the Eurasian Basin north of Svalbard collected between January and June 2015 from the N-ICE2015 drifting expedition are presented. The unique winter observations are a key contribution to existing climatologies of the Arctic Ocean, and show a ∼100 m deep winter mixed layer likely due to high sea ice growth rates in local leads. Current observations for the upper ∼200 m show mostly a barotropic flow, enhanced over the shallow Yermak Plateau. The two branches of inflowing Atlantic Water are partly captured, confirming that the outer Yermak Branch follows the perimeter of the plateau, and the inner Svalbard Branch the coast. Atlantic Water observed to be warmer and shallower than in the climatology, is found directly below the mixed layer down to 800 m depth, and is warmest along the slope, while its properties inside the basin are quite homogeneous. From late May onwards, the drift was continually close to the ice edge and a thinner surface mixed layer and shallower Atlantic Water coincided with significant sea ice melt being observed.

Description: Our study followed the seasonal cycling of soluble (SFe), colloidal (CFe), dissolved (DFe), total dissolvable (TDFe), labile particulate (LPFe) and total particulate (TPFe) iron in the Celtic Sea (NE Atlantic Ocean). Preferential uptake of SFe occurred during the spring bloom, preceding the removal of CFe. Uptake and export of Fe during the spring bloom, coupled with a reduction in vertical exchange, led to Fe deplete surface waters (〈0.2 nM DFe; 0.11 nM LPFe, 0.45 nM TDFe, 1.84 nM TPFe) during summer stratification. Below the seasonal thermocline, DFe concentrations increased from spring to autumn, mirroring NO3- and consistent with supply from remineralised sinking organic material, and cycled independently of particulate Fe over seasonal timescales. These results demonstrate that summer Fe availability is comparable to the seasonally Fe limited Ross Sea shelf, and therefore is likely low enough to affect phytoplankton growth and species composition.

Description: Current climate models disagree on how much carbon dioxide land ecosystems take up for photosynthesis. Tracking the stronger carbonyl sulfide signal could help.

Description: Oceanic dimethyl sulfide (DMS) is of interest due to its critical influence on atmospheric sulfur compounds in the marine atmosphere and its hypothesized significant role in global climate. High-resolution shipboard underway measurements of surface seawater DMS and the partial pressure of carbon dioxide (pCO2) were conducted in the Atlantic Ocean and Indian Ocean sectors of the Southern Ocean (SO), the southeast Indian Ocean, and the northwest Pacific Ocean from February to April 2014 during the 30th Chinese Antarctic Research Expedition. The SO, particularly in the region south of 58°S, had the highest mean surface seawater DMS concentration of 4.1 ± 8.3 nM (ranged from 0.1 to 73.2 nM) and lowest mean seawater pCO2 level of 337 ± 50 μatm (ranged from 221 to 411 μatm) over the entire cruise. Significant variations of surface seawater DMS and pCO2 in the seasonal ice zone (SIZ) of SO were observed, which are mainly controlled by biological process and sea ice activity. We found a significant negative relationship between DMS and pCO2 in the SO SIZ using 0.1° resolution, [DMS] seawater = -0.160 [pCO2] seawater + 61.3 (r2 = 0.594, n = 924, p 〈 0.001). We anticipate that the relationship may possibly be utilized to reconstruct the surface seawater DMS climatology in the SO SIZ. Further studies are necessary to improve the universality of this approach. Key Points: • The characteristics of surface water DMS and pCO2 distributions from the Southern Ocean to northwest Pacific Ocean are investigated • The correlations between DMS, pCO2, and environmental parameters are analyzed • Anticorrelation between DMS and pCO2 is found in the seasonal ice zone of the Southern Ocean

Description: Natural gas hydrates are considered a potential resource for gas production on industrial scales. Gas hydrates contribute to the strength and stiffness of the hydrate-bearing sediments. During gas production, the geomechanical stability of the sediment is compromised. Due to the potential geotechnical risks and process management issues, the mechanical behavior of the gas hydrate-bearing sediments needs to be carefully considered. In this study, we describe a coupling concept that simplifies the mathematical description of the complex interactions occurring during gas production by isolating the effects of sediment deformation and hydrate phase changes. Central to this coupling concept is the assumption that the soil grains form the load-bearing solid skeleton, while the gas hydrate enhances the mechanical properties of this skeleton. We focus on testing this coupling concept in capturing the overall impact of geomechanics on gas production behavior though numerical simulation of a high-pressure isotropic compression experiment combined with methane hydrate formation and dissociation. We consider a linear-elastic stress-strain relationship because it is uniquely defined and easy to calibrate. Since, in reality, the geomechanical response of the hydrate-bearing sediment is typically inelastic and is characterized by a significant shear-volumetric coupling, we control the experiment very carefully in order to keep the sample deformations small and well within the assumptions of poroelasticity. The closely coordinated experimental and numerical procedures enable us to validate the proposed simplified geomechanics-to-flow coupling, and set an important precursor toward enhancing our coupled hydro-geomechanical hydrate reservoir simulator with more suitable elastoplastic constitutive models.

Description: The Summer East Atlantic (SEA) mode is the second dominant mode of summer low-frequency variability in the Euro-Atlantic region. Using reanalysis data, we show that SEA-related circulation anomalies significantly influence temperatures and precipitation over Europe. We present evidence that part of the interannual SEA variability is forced by diabatic heating anomalies of opposing signs in the tropical Pacific and Caribbean that induce an extratropical Rossby wave train. This precipitation dipole is related to SST anomalies characteristic of the developing ENSO phases. Seasonal hindcast experiments forced with observed sea surface temperatures (SST) exhibit skill at capturing the interannual SEA variability corroborating the proposed mechanism and highlighting the possibility for improved prediction of boreal summer variability. Our results indicate that tropical forcing of the SEA likely played a role in the dynamics of the 2015 European heat wave.

Description: Polyextremophiles are present in a wide variety of extreme environments in which they must overcome various hostile conditions simultaneously such as high UVB radiation, extreme pHs and temperatures, elevated salt and heavy-metal concentration, low-oxygen pressure and scarce nutrients. High-altitude Andean lakes (HAALs; between 2000 and 4000 m) are one example of these kinds of ecosystems suffering from the highest total solar and UVB radiation on Earth where an abundant and diverse polyextremophilic microbiota was reported. In this work, we performed the first extensive isolation of UV-resistant actinobacteria from soils, water, sediments and modern stromatolites at HAALs. Based on the 16S rRNA sequence, the strains were identified as members of the genera Streptomyces, Micrococcus, Nesterenkonia, Rhodococcus, Microbacterium, Kocuria, Arthrobacter, Micromonospora, Blastococcus, Citrococcus and Brevibacterium. Most isolates displayed resistance to multiple environmental stress factors confirming their polyextremophilic nature and were able to produce effective antimicrobial compounds. HAALs constitute a largely unexplored repository of UV-resistant actinobacteria, with high potential for the biodiscovery of novel natural products.

Description: Two variants of sea-surface temperature (SST) dipole indices for the South Atlantic Ocean (SAO) has been previously described representing: (1) the South Atlantic subtropical dipole (SASD) supposedly peaking in austral summer and (2) the SAO dipole (SAOD) in winter. In this study, we present the analysis of observational data sets (1985–2014) showing the SASD and SAOD as largely constituting the same mode of ocean–atmosphere interaction reminiscent of the SAOD structure peaking in winter. Indeed, winter is the only season in which the inverse correlation between the northern and southern poles of both indices is statistically significant. The observed SASD and SAOD indices exhibit robust correlations (P ≤ 0.001) in all seasons and these are reproduced by 54 of the 63 different models of the Coupled Models Intercomparison Project analysed. Their robust correlations notwithstanding the SASD and SAOD indices appear to better capture different aspects of SAO climate variability and teleconnections

Description: Contents 670 I. 671 II. 671 III. 676 IV. 678 678 References 678 SUMMARY: Biotic interactions underlie life's diversity and are the lynchpin to understanding its complexity and resilience within an ecological niche. Algal biologists have embraced this paradigm, and studies building on the explosive growth in omics and cell biology methods have facilitated the in-depth analysis of nonmodel organisms and communities from a variety of ecosystems. In turn, these advances have enabled a major revision of our understanding of the origin and evolution of photosynthesis in eukaryotes, bacterial-algal interactions, control of massive algal blooms in the ocean, and the maintenance and degradation of coral reefs. Here, we review some of the most exciting developments in the field of algal biotic interactions and identify challenges for scientists in the coming years. We foresee the development of an algal knowledgebase that integrates ecosystem-wide omics data and the development of molecular tools/resources to perform functional analyses of individuals in isolation and in populations. These assets will allow us to move beyond mechanistic studies of a single species towards understanding the interactions amongst algae and other organisms in both the laboratory and the field.

Description: The region encompassing the Kuroshio Extension (KE) in the Northwestern Pacific Ocean (25°N–45°N and 130°E–180°E) is one of the most eddy-energetic regions of the global ocean. The three-dimensional structures and transports of mesoscale eddies in this region are comprehensively investigated by combined use of satellite data and Argo profiles. With the allocation of Argo profiles inside detected eddies, the spatial variations of structures of eddy temperature and salinity anomalies are analyzed. The results show that eddies predominantly have subsurface (near-surface) intensified temperature and salinity anomalies south (north) of the KE jet, which is related to different background stratifications between these regions. A new method based on eddy trajectories and the inferred three-dimensional eddy structures is proposed to estimate heat and salt transports by eddy movements in a Lagrangian framework. Spatial distributions of eddy transports are presented over the vicinity of the KE for the first time. The magnitude of eddy-induced meridional heat (freshwater volume) transport is on the order of 0.01 PW (103 m3/s). The eddy heat transport divergence results in an oceanic heat loss south and heat gain north of the KE, thereby reinforcing and counteracting the oceanic heat loss from air-sea fluxes south and north of the KE jet, respectively. It also suggests a poleward heat transport across the KE jet due to eddy propagation.

Description: Marine zooplankton comprise a phylogenetically and functionally diverse assemblage of protistan and metazoan consumers that occupy multiple trophic levels in pelagic food webs. Within this complex network, carbon flows via alternative zooplankton pathways drive temporal and spatial variability in production-grazing coupling, nutrient cycling, export, and transfer efficiency to higher trophic levels. We explore current knowledge of the processing of zooplankton food ingestion by absorption, egestion, respiration, excretion, and growth (production) processes. On a global scale, carbon fluxes are reasonably constrained by the grazing impact of microzooplankton and the respiratory requirements of mesozooplankton but are sensitive to uncertainties in trophic structure. The relative importance, combined magnitude, and efficiency of export mechanisms (mucous feeding webs, fecal pellets, molts, carcasses, and vertical migrations) likewise reflect regional variability in community structure. Climate change is expected to broadly alter carbon cycling by zooplankton and to have direct impacts on key species.

Description: The historical developments are reviewed that have led from a bottom-up responsibility initiative of concerned scientists to the emergence of a nationwide interdisciplinary Priority Program on the assessment of Climate Engineering (CE), funded by the German Research Foundation (DFG). Given the perceived lack of comprehensive and comparative appraisals of different CE methods, the Priority Program was designed to encompass both solar radiation management (SRM) and carbon dioxide removal (CDR) ideas, and to cover the atmospheric, terrestrial and oceanic realm. First key findings obtained by the ongoing Priority Program are summarized and reveal that compared to earlier assessments, such as the 2009 Royal Society report, more detailed investigations tend to indicate less efficiency, lower effectiveness and often lower safety. Emerging research trends are discussed in the context of the recent Paris agreement to limit global warming to less than two degrees and the associated increasing reliance on negative emission technologies. Our results show then when deployed at scales large enough to have a significant impact on atmospheric CO2, even CDR methods such as afforestation – often perceived as ‘benign’ – can have substantial side effects and may raise severe ethical, legal and governance issues. We suppose that before being deployed at climatically relevant scales, any negative-emission or climate engineering method will require careful analysis of efficiency, effectiveness and undesired side effects.

Description: Selecting appropriate indicators is essential to aggregate the information provided by climate model outputs into a manageable set of relevant metrics on which assessments of climate engineering (CE) can be based. From all the variables potentially available from climate models, indicators need to be selected that are able to inform scientists and society on the development of the Earth system under CE, as well as on possible impacts and side effects of various ways of deploying CE or not. However, the indicators used so far have been largely identical to those used in climate change assessments and do not visibly reflect the fact that indicators for assessing CE (and thus the metrics composed of these indicators) may be different from those used to assess global warming. Until now, there has been little dedicated effort to identifying specific indicators and metrics for assessing CE. We here propose that such an effort should be facilitated by a more decision-oriented approach and an iterative procedure in close interaction between academia, decision makers, and stakeholders. Specifically, synergies and trade-offs between social objectives reflected by individual indicators, as well as decision-relevant uncertainties should be considered in the development of metrics, so that society can take informed decisions about climate policy measures under the impression of the options available, their likely effects and side effects, and the quality of the underlying knowledge base.

Description: The consequences of emerging marine diseases on the evolutionary trajectories of affected host populations in the marine realm are largely unexplored. Evolution in response to natural selection depends on the genetic variation of the traits under selection and the interaction of these traits with the environment (GxE). However, in the case of diseases, genotypes of pathogens add another dimension to this interaction. Therefore, the study of disease resistance needs to be extended to the interaction of host genotype, pathogen genotype and environment (GxGxE). In the present study we used a full-sib breeding design crossing two genetically differentiated populations of the Pacific oyster Crassostrea gigas (Thunberg, 1793), to determine the influence of host genotype, pathogen genotype and temperature on disease resistance. Based on a controlled infection experiment on two early life stages, i.e. D-larvae and Pediveliger larvae at elevated and ambient water temperatures we estimated disease resistance to allopatric and sympatric Vibrio sp. by measuring survival and growth within and between genetically differentiated oyster populations. In both populations survival was higher upon infection with sympatric Vibrio sp. indicating that disease resistance has a genetic basis and is dependent on host genotype. In addition we observed a significant GxGxE effect in D-larvae, where contrary to expectations, disease resistance was higher at warm than at cold temperatures. Using thermal reaction norms, we could further show, that disease resistance is an environment dependent trait with high plasticity, which indicates the potential for a fast acclimatization to changing environmental conditions. These population specific reaction norms disappeared in hybrid crosses between both populations which demonstrates that admixture between genetically differentiated populations can influence GxGxE interactions on larger scales.

Description: Over the past 17 years, the western boundary current system of the Labrador Sea has been closely observed by maintaining the 53°N observatory (moorings and shipboard station data) measuring the top-to-bottom flow field offshore from the Labrador shelf break. Volume transports for the North Atlantic Deep Water (NADW) components were calculated using different methods, including gap filling procedures for deployment periods with suboptimal instrument coverage. On average the Deep Western Boundary Current (DWBC) carries 30.2 ± 6.6 Sv of NADW southward, which are almost equally partitioned between Labrador Sea Water (LSW, 14.9 ± 3.9 Sv) and Lower North Atlantic Deep Water (LNADW, 15.3 ± 3.8 Sv). The transport variability ranges from days to decades, with the most prominent multiyear fluctuations at interannual to near decadal time scales (±5 Sv) in the LNADW overflow water mass. These long-term fluctuations appear to be in phase with the NAO-modulated wind fluctuations. The boundary current system off Labrador occurs as a conglomerate of nearly independent components, namely, the shallow Labrador Current, the weakly sheared LSW range, and the deep baroclinic, bottom-intensified current core of the LNADW, all of which are part of the cyclonic Labrador Sea circulation. This structure is relatively stable over time, and the 120 km wide boundary current is constrained seaward by a weak counterflow which reduces the deep water export by 10–15%.

Description: Theory of local adaptation predicts that nonadapted migrants will suffer increased costs compared to local residents. Ultimately this process can result in the reduction of gene flow and culminate in speciation. Here, we experimentally investigated the relative fitness of migrants in foreign habitats, focusing on diverging lake and river ecotypes of three-spined sticklebacks. A reciprocal transplant experiment performed in the field revealed asymmetric costs of migration: whereas mortality of river fish was increased under lake conditions, lake migrants suffered from reduced growth relative to river residents. Selection against migrants thus involved different traits in each habitat but generally contributed to bidirectional reduction in gene flow. Focusing particularly on the parasitic environments, migrant fish differed from resident fish in the parasite community they harboured. This pattern correlated with both cellular phenotypes of innate immunity as well as with allelic variation at the genes of the major histocompatibility complex. In addition to showing the costs of migration in three-spined sticklebacks, this study highlights the role of asymmetric selection particularly from parasitism in genotype sorting and in the emergence of local adaptation.

Description: Trans-generational plasticity is the adjustment of phenotypes to changing habitat conditions that persist longer than the individual lifetime. Fitness benefits (adaptive TGP) are expected upon matching parent-offspring environments. In a global change scenario, several performance-related environmental factors are changing simultaneously. This lowers the predictability of offspring environmental conditions, potentially hampering the benefits of trans-generational plasticity. For the first time, we here explore how the combination of an abiotic and a biotic environmental factor in the parental generation plays out as trans-generational effect in the offspring. We fully reciprocally exposed the parental generation of the pipefish Syngnathus typhle to an immune challenge and elevated temperatures simulating a naturally occurring heatwave. Upon mating and male pregnancy, offspring were kept in ambient or elevated temperature regimes combined with a heat-killed bacterial epitope treatment. Differential gene expression (immune genes and DNA- and histone-modification genes) suggests that the combined change of an abiotic and a biotic factor in the parental generation had interactive effects on offspring performance, the temperature effect dominated over the immune challenge impact. The benefits of certain parental environmental conditions on offspring performance did not sum up when abiotic and biotic factors were changed simultaneously supporting that available resources that can be allocated to phenotypic trans-generational effects are limited. Temperature is the master regulator of trans-generational phenotypic plasticity, which potentially implies a conflict in the allocation of resources towards several environmental factors. This asks for a reassessment of trans-generational plasticity as a short-term option to buffer environmental variation in the light of climate change.

Description: Primary productivity is limited by the availability of nitrogen (N) in most of the coastal Arctic, as a large portion of N is released by the spring freshet and completely consumed during the following summer. Thus, understanding the fate of riverine nitrogen is critical to identify the link between dissolved nitrogen dynamic and coastal primary productivity to foresee upcoming changes in the Arctic seas, such as increase riverine discharge and permafrost thaw. Here, we provide a field-based study of nitrogen dynamic over the Laptev Sea shelf based on isotope geochemistry. We demonstrate that while most of the nitrate found under the surface fresh water layer is of remineralized origin, some of the nitrate originates from atmospheric input and was probably transported at depth by the mixing of brine-enriched denser water during sea-ice formation. Moreover, our results suggest that riverine dissolved organic nitrogen (DON) represents up to 6 times the total riverine release of nitrate and that about 62 to 76% of the DON is removed within the shelf waters. This is a crucial information regarding the near-future impact of climate change on primary productivity in the Eurasian coastal Arctic.

Description: The identification of native sources and vectors of introduced species informs their ecological and evolutionary history and may guide policies that seek to prevent future introductions. Population genetics provides a powerful set of tools to identify origins and vectors. However, these tools can mislead when the native range is poorly sampled or few molecular markers are used. Here, we traced the introduction of the Asian seaweed Gracilaria vermiculophylla (Rhodophyta) into estuaries in coastal western North America, the eastern United States, Europe, and northwestern Africa by genotyping more than 2,500 thalli from 37 native and 53 non-native sites at mitochondrial cox1 and 10 nuclear microsatellite loci. Overall, greater than 90% of introduced thalli had a genetic signature similar to thalli sampled from the coastline of northeastern Japan, strongly indicating this region served as the principal source of the invasion. Notably, northeastern Japan exported the vast majority of the oyster Crassostrea gigas during the 20th century. The preponderance of evidence suggests G. vermiculophylla may have been inadvertently introduced with C. gigas shipments and that northeastern Japan is a common source region for estuarine invaders. Each invaded coastline reflected a complex mix of direct introductions from Japan and secondary introductions from other invaded coastlines. The spread of G. vermiculophylla along each coastline was likely facilitated by aquaculture, fishing, and boating activities. Our ability to document a source region was enabled by a robust sampling of locations and loci that previous studies lacked and strong phylogeographic structure along native coastlines.

Description: Marine controlled source electromagnetic (CSEM) data have been collected to investigate methane seep sites and associated gas hydrate deposits at Opouawe Bank on the southern tip of the Hikurangi Margin, New Zealand. The bank is located in about 1000 m water depth within the gas hydrate stability field. The seep sites are characterized by active venting and typical methane seep fauna accompanied with patchy carbonate outcrops at the seafloor. Below the seeps, gas migration pathways reach from below the bottom-simulating reflector (at around 380 m sediment depth) toward the seafloor, indicating free gas transport into the shallow hydrate stability field. The CSEM data have been acquired with a seafloor-towed, electric multi-dipole system measuring the inline component of the electric field. CSEM data from three profiles have been analyzed by using 1-D and 2-D inversion techniques. High-resolution 2-D and 3-D multichannel seismic data have been collected in the same area. The electrical resistivity models show several zones of highly anomalous resistivities (〉50 Ωm) which correlate with high amplitude reflections located on top of narrow vertical gas conduits, indicating the coexistence of free gas and gas hydrates within the hydrate stability zone. Away from the seeps the CSEM models show normal background resistivities between ~1 and 2 Ωm. Archie's law has been applied to estimate gas/gas hydrate saturations below the seeps. At intermediate depths between 50 and 200 m below seafloor, saturations are between 40 and 80% and gas hydrate may be the dominating pore filling constituent. At shallow depths from 10 m to the seafloor, free gas dominates as seismic data and gas plumes suggest.

Description: Many species of Indo-Pacific holobenthic foraminifera have been introduced and successfully established sustainable populations in the Mediterranean Sea over the past few decades. However, known natural and anthropogenic vectors do not explain how these species were introduced long distances from their origin. We present evidence for a novel marine bioinvasion vector explaining this long-distance transport and introduction using both contemporary field and historical analyses. In 2015–2016, we found living specimens of 29 foraminiferal species in the fecal pellets of two Red Sea herbivorous rabbitfish—Siganus rivulatus and Siganus luridus in the Mediterranean. In our historical analysis, we found 34 foraminiferal species in preserved Red Sea rabbitfish specimens, dating between 1967 and 1975. In addition, we found congruent propagation patterns of the non-indigenous rabbitfish and foraminifera, lagging 4–11 yrs between discoveries, respectively. Predation of marine benthos by non-indigenous fish, followed by incomplete digestion and defecation of viable individuals, comprise the main introduction vector of these organisms into novel environments.

Description: Biological invasions are worldwide phenomena that have reached alarming levels among aquatic species. There are key challenges to understand the factors behind invasion propensity of non-native populations in invasion biology. Interestingly, interpretations cannot be expanded to higher taxonomic levels due to the fact that in the same genus, there are species that are notorious invaders and those that never spread outside their native range. Such variation in invasion propensity offers the possibility to explore, at fine-scale taxonomic level, the existence of specific characteristics that might predict the variability in invasion success. In this work, we explored this possibility from a molecular perspective. The objective was to provide a better understanding of the genetic diversity distribution in the native range of species that exhibit contrasting invasive propensities. For this purpose, we used a total of 784 sequences of the cytochrome c oxidase subunit I of mitochondrial DNA (mtDNA-COI) collected from seven Gammaroidea, a superfamily of Amphipoda that includes species that are both successful invaders (Gammarus tigrinus, Pontogammarus maeoticus, and Obesogammarus crassus) and strictly restricted to their native regions (Gammarus locusta, Gammarus salinus, Gammarus zaddachi, and Gammarus oceanicus). Despite that genetic diversity did not differ between invasive and non-invasive species, we observed that populations of non-invasive species showed a higher degree of genetic differentiation. Furthermore, we found that both geographic and evolutionary distances might explain genetic differentiation in both non-native and native ranges. This suggests that the lack of population genetic structure may facilitate the distribution of mutations that despite arising in the native range may be beneficial in invasive ranges. The fact that evolutionary distances explained genetic differentiation more often than geographic distances points toward that deep lineage divergence holds an important role in the distribution of neutral genetic diversity.

Description: Receiver functions (RF) have been used for several decades to study structures beneath seismic stations. Although most available stations are deployed on-shore, the number of ocean bottom station (OBS) experiments has increased in recent years. Almost all OBSs have to deal with higher noise levels and a limited deployment time (∼1 year), resulting in a small number of usable records of teleseismic earthquakes. Here, we use OBSs deployed as mid-aperture array in the deep ocean (4.5-5.5 km water depth) of the eastern mid-Atlantic. We use evaluation criteria for OBS data and beam forming to enhance the quality of the RFs. Although some stations show reverberations caused by sedimentary cover, we are able to identify the Moho signal, indicating a normal thickness (5-8 km) of oceanic crust. Observations at single stations with thin sediments (300-400 m) indicate that a probable sharp lithosphere-asthenosphere boundary (LAB) might exist at a depth of ∼70-80 km which is in line with LAB depth estimates for similar lithospheric ages in the Pacific. The mantle discontinuities at ∼410 km and ∼660 km are clearly identifiable. Their delay times are in agreement with PREM. Overall the usage of beam formed earthquake recordings for OBS RF analysis is an excellent way to increase the signal quality and the number of usable events.

Description: There is currently no consensus on how humans are affecting the marine nitrogen (N) cycle, which limits marine biological production and CO2 uptake. Anthropogenic changes in ocean warming, deoxygenation, and atmospheric N deposition can all individually affect the marine N cycle and the oceanic production of the greenhouse gas nitrous oxide (N2O). However, the combined effect of these perturbations on marine N cycling, ocean productivity, and marine N2O production is poorly understood. Here we use an Earth system model of intermediate complexity to investigate the combined effects of estimated 21st century CO2 atmospheric forcing and atmospheric N deposition. Our simulations suggest that anthropogenic perturbations cause only a small imbalance to the N cycle relative to preindustrial conditions (∼+5 Tg N y−1 in 2100). More N loss from water column denitrification in expanded oxygen minimum zones (OMZs) is counteracted by less benthic denitrification, due to the stratification-induced reduction in organic matter export. The larger atmospheric N load is offset by reduced N inputs by marine N2 fixation. Our model predicts a decline in oceanic N2O emissions by 2100. This is induced by the decrease in organic matter export and associated N2O production and by the anthropogenically driven changes in ocean circulation and atmospheric N2O concentrations. After comprehensively accounting for a series of complex physical-biogeochemical interactions, this study suggests that N flux imbalances are limited by biogeochemical feedbacks that help stabilize the marine N inventory against anthropogenic changes. These findings support the hypothesis that strong negative feedbacks regulate the marine N inventory on centennial time scales.

Description: When volcanic mountains slide into the sea, they trigger tsunamis. How big are these waves, and how far away can they do damage? Ritter Island provides some answers.

Description: Proxy data suggest the onset of Northern Hemisphere glaciation during the Plio-Pleistocene transition from 3.2 to 2.5 Ma resulted in enhanced climate variability at the obliquity (41 kyr) frequency. Here, we investigate the influence of the expanding Greenland ice sheet (GrIS) on the mean climate and obliquity-related variability in a series of climate model simulations. These suggest that an expanding GrIS weakens the Atlantic Meridional Overturning Circulation (AMOC) by ~1 Sv, mainly due to reduced heat loss in the Greenland-Iceland-Norwegian Sea. Moreover, the growing GrIS amplifies the Hadley circulation response to obliquity forcing driving variations in freshwater export from the tropical Atlantic and in turn variations of the AMOC. The stronger AMOC response to obliquity forcing, by about a factor of two, results in a stronger global-mean near-surface temperature response. We conclude that the AMOC response to obliquity forcing is important to understand the enhanced climate variability at the obliquity frequency during the Plio-Pleistocene transition.

Description: In this pilot study we link the yield of industrial fisheries to changes in the zooplankton mortality in an idealized way accounting for different target species (planktivorous fish—decreased zooplankton mortality; large predators—increased zooplankton mortality). This indirect approach is used in a global coupled biogeochemistry circulation model to estimate the range of the potential impact of industrial fisheries on marine biogeochemistry. The simulated globally integrated response on phytoplankton and primary production is in line with expectations—a high (low) zooplankton mortality results in a decrease (increase) of zooplankton and an increase (decrease) of phytoplankton. In contrast, the local response of zooplankton and phytoplankton depends on the region under consideration: In nutrient-limited regions, an increase (decrease) in zooplankton mortality leads to a decrease (increase) in both zooplankton and phytoplankton biomass. In contrast, in nutrient-replete regions, such as upwelling regions, we find an opposing response: an increase (decrease) of the zooplankton mortality leads to an increase (decrease) in both zooplankton and phytoplankton biomass. The results are further evaluated by relating the potential fisheries-induced changes in zooplankton mortality to those driven by CO2 emissions in a business-as-usual 21st century emission scenario. In our idealized case, the potential fisheries-induced impact can be of similar size as warming-induced changes in marine biogeochemistry.

Description: Back-arc spreading centers (BASCs) form a distinct class of ocean spreading ridges distinguished by steep along-axis gradients in spreading rate and by additional magma supplied through subduction. These characteristics can affect the population and distribution of hydrothermal activity on BASCs compared to mid-ocean ridges (MORs). To investigate this hypothesis, we comprehensively explored 600 km of the southern half of the Mariana BASC. We used water column mapping and seafloor imaging to identify 19 active vent sites, an increase of 13 over the current listing in the InterRidge Database (IRDB), on the bathymetric highs of 7 of the 11 segments. We identified both high and low (i.e., characterized by a weak or negligible particle plume) temperature discharge occurring on segment types spanning dominantly magmatic to dominantly tectonic. Active sites are concentrated on the two southernmost segments, where distance to the adjacent arc is shortest (〈40 km), spreading rate is highest (〉48 mm/yr), and tectonic extension is pervasive. Re-examination of hydrothermal data from other BASCs supports the generalization that hydrothermal site density increases on segments 〈90 km from an adjacent arc. Although exploration quality varies greatly among BASCs, present data suggest that, for a given spreading rate, the mean spatial density of hydrothermal activity varies little between MORs and BASCs. The present global database, however, may be misleading. On both BASCs and MORs, the spatial density of hydrothermal sites mapped by high-quality water-column surveys is 2–7 times greater than predicted by the existing IRDB trend of site density versus spreading rate.

Description: The evolutionary trajectory of early complex life on Earth is interpreted largely from the fossils of the Precambrian soft-bodied Ediacara Biota, which appeared and evolved during a time of dynamic biogeochemical and environmental fluctuation in the global ocean. The Ediacara Biota is historically divided into three successive Assemblages—the Avalon, the White Sea, and the Nama—which are marked by the appearance of novel biological traits and ecological strategies. In particular, the younger White Sea and Nama Assemblages record a “second wave” of ecological innovations, which included not only the development of uniquely Ediacaran body plans and ecologies, such as matground adaptations, but also the dual emergence of bilaterian-grade animals and Phanerozoic-style ecological innovations, including spatial heterogeneity, complex reproductive strategies, ecospace utilization, motility, and substrate competition. The late Ediacaran was an evolutionarily dynamic time characterized by strong environmental control over the distribution of taxa in time and space.

Description: The oceanic crustal and uppermost lithospheric mantle structure across the Gloria Fault (GF) transcurrent plate boundary between Africa and Eurasia in the Northeast Atlantic is investigated based on seismic reflection, seismic refraction and wide-angle reflection data. This experiment used 18 ocean bottom stations along an N–S 150 km long traverse together with acquisition of a multichannel seismic reflection profile. Modeling of P and S seismic waves and gravimetric anomalies allowed estimation of P- and S-wave velocities, density, Poisson's ratio and discussion of a compositional model. A five-layer model is proposed in which layers 1–3 correspond to normal sediments through typical oceanic crust layers 2 and 3. Layer 5 yielded mantle velocities above 7.9 km s−1. Layer 4 with 4 km of thickness has Vp velocities between 7.1 and 7.4 km s−1 and is clearly separated from typical oceanic crust and mantle layers. Comparison with natural analogues and published lab measurements suggest that layer 4 can be a mix of lithologies that comply with the estimated P and S velocities and computed Poisson's ratio and densities, such as, olivine cumulates, peridotite, gabbro and hydrated mantle. We favour the tectonic process that produces secondary porosity from which results serpentinization due to sea water circulation in fractures. Structural and seismic stratigraphic interpretation of the reflection profile shows that Neogene to recent tectonic deformation on this segment of the plate boundary concentrated on the southern side of the GF, that is, the Africa plate.

Description: Overriding plate topography provides constraints on subduction zone geodynamics. We investigate its evolution using fully dynamic laboratory models of subduction with techniques of stereoscopic photogrammetry and particle image velocimetry. Model results show that the topography is characterized by an area of forearc dynamic subsidence, with a magnitude scaling to 1.44–3.97 km in nature, and a local topographic high between the forearc subsided region and the trench. These topographic features rapidly develop during the slab free‐sinking phase and gradually decrease during the steady state slab rollback phase. We propose that they result from the variation of the vertical component of the trench suction force along the subduction zone interface, which gradually increases with depth and results from the gradual slab steepening during the initial transient slab sinking phase. The downward mantle flow in the nose of the mantle wedge plays a minor role in driving forearc subsidence.

Description: We provide high-resolution foraminiferal stable carbon isotope (δ13C) records from the subarctic Pacific and Eastern Equatorial Pacific (EEP) to investigate circulation dynamics between the extra-tropical and tropical North Pacific during the past 60 kyr. We measured the δ13C composition of the epibenthic foraminiferal species Cibicides lobatulus from a shallow sediment core recovered from the western Bering Sea (SO201-2-101KL; 58°52.52’N, 170°41.45’E; 630 m water depth) to reconstruct past ventilation changes close to the source region of Glacial North Pacific Intermediate Water (GNPIW). Information regarding glacial changes in the δ13C of sub-thermocline water masses in the EEP is derived from the deep-dwelling planktonic foraminifera Globorotaloides hexagonus at ODP Site 1240 (00°01.31’N, 82°27.76’W; 2921 m water depth). Apparent similarities in the long-term evolution of δ13C between GNPIW, intermediate waters in the eastern tropical North Pacific and sub-thermocline water masses in the EEP suggest the expansion of relatively 13C-depleted, nutrient-enriched, and northern-sourced intermediate waters to the equatorial Pacific under glacial conditions. Further, it appears that additional influence of GNPIW to the tropical Pacific is consistent with changes in nutrient distribution and biological productivity in surface-waters of the glacial EEP. Our findings highlight potential links between North Pacific mid-depth circulation changes, nutrient cycling, and biological productivity in the equatorial Pacific under glacial boundary conditions.

Description: A multi-mode, linear reduced-gravity model, driven by ERA-Interim monthly mean wind stress anomalies, is used to investigate interannual variability in tropical Pacific sea level as seen in satellite altimeter data. The model output is fitted to the altimeter data along the equator, in order to derive the vertical profile for the model forcing, showing that a signature from modes higher than mode six cannot be extracted from the altimeter data. It is shown that the model has considerable skill at capturing interannual sea level variability both on and off the equator. The correlation between modelled and satellite-derived sea level data exceeds 0.8 over a wide range of longitudes along the equator and readily captures the observed ENSO events. Overall, the combination of the first, second, third and fifth modes can provide a robust estimate of the interannual sea level variability, the second mode being dominant. A remarkable feature of both the model and the altimeter data is the presence of a pivot point in the western Pacific on the equator. We show that the westward displacement of the pivot point from the centre of the basin is strongly influenced by the fact that most of the wind stress variance is found in the western part of the basin. We also show that the Sverdrup transport is not fundamental to the dynamics of the recharge/discharge mechanism in our model, although the spatial structure of the wind forcing does play a role in setting the amplitude of the “warm water volume”.

Description: El Niño-Southern Oscillation (ENSO) in the Pacific is asymmetric for warm and cold events with respect to amplitude, spatial patterns and temporal evolution. Here the symmetry of the Atlantic Niño mode, which many previous studies have argued is governed by atmosphere–ocean dynamics similar to those of ENSO, is investigated using two different ocean reanalysis products. Calculation of Bjerknes feedback terms for the Pacific reveals a pronounced asymmetry between warm and cold events, though unlike most previous studies, the largest asymmetry is found in the relationship between eastern Pacific thermocline depth and SST anomalies. For the Atlantic, cold events are effectively mirror images of warm events with Bjerknes feedbacks of similar strength. The analysis supports not only the conclusion that Atlantic Niños are more symmetric than ENSO, but the hypothesis itself that the Bjerknes feedback is operative in the Atlantic given the strength of the relationship between the key variables involved.

Description: At the eastern end of the Azores-Gloria transform fault system to the southwest of Portugal, the plate boundary between Africa and Iberia is a region where deformation is accommodated over a wide tectonically-active area. The region has unleashed large earthquakes and tsunamis, including the Mw ~ 8.5 Great Lisbon earthquake of 1755. Although the source region of the 1755 earthquake is still disputed, most proposals include a source location in the vicinity of the Horseshoe Abyssal Plain (HAP), which is bounded by the 5000 m high Gorringe Bank (GB). In this study we characterise seismic activity in the region using data recorded by two local networks of ocean-bottom seismometers (OBS). The networks were deployed in the eastern HAP and at the GB. The dataset allowed the detection of 160 local earthquakes. These earthquakes cluster around the GB, to the SW of Cabo Sao Vicente, and in the HAP. Focal depths indicate deep-seated earthquakes, with depths increasing from 20-35 km (mean of 26.1 ± 7.2 km) at the GB to 15-45 km (mean 31.5 km ± 10.5 km) under the HAP. Seismic activity thus extends down to levels that are deeper than those mapped by active seismic profiling, with the majority of events occurring within the mantle. Thermal modelling suggests that temperatures of approximately 600 °C characterise the base of the seismogenic brittle lithosphere at ~45 km depth. The large source depth and thermal structure supports previous suggestions that catastrophic seismic rupture through the lithospheric mantle may indeed occur in the area.

Description: Despite a growing literature on the climate response to solar geoengineering – proposals to cool the planet by increasing the planetary albedo – there has been little published on the impacts of solar geoengineering on natural and human systems such as agriculture, health, water resources, and ecosystems. An understanding of the impacts of different scenarios of solar geoengineering deployment will be crucial for informing decisions on whether and how to deploy it. Here we review the current state of knowledge about impacts of a solar geoengineered climate and identify major research gaps. We suggest that a thorough assessment of the climate impacts of a range of scenarios of solar geoengineering deployment is needed and can build upon existing frameworks. However, solar geoengineering poses a novel challenge for climate impacts research as the manner of deployment could be tailored to pursue different objectives making possible a wide range of climate outcomes. We present a number of ideas for approaches to extend the survey of climate impacts beyond standard scenarios of solar geoengineering deployment to address this challenge. Reducing the impacts of climate change is the fundamental motivator for emissions reductions and for considering whether and how to deploy solar geoengineering. This means that the active engagement of the climate impacts research community will be important for improving the overall understanding of the opportunities, challenges and risks presented by solar geoengineering.

Description: The supply and bioavailability of iron (Fe) controls primary productivity and N2-fixation in large parts of the global ocean. An important, yet poorly quantified, source to the ocean is particulate Fe (pFe). Here we present the first combined dataset of particulate, labile-particulate (L-pFe) and dissolved Fe (dFe) from the (sub)-tropical North Atlantic. We show a strong relationship between L-pFe and dFe, indicating a dynamic equilibrium between these two phases whereby particles ‘buffer’ dFe and maintain the elevated concentrations observed. Moreover, L-pFe can increase the overall ‘available’ (L-pFe + dFe) Fe pool by up to 55%. The lateral shelf flux of this available Fe was similar in magnitude to observed soluble aerosol-Fe deposition, a comparison that has not been previously considered. These findings demonstrate that L-pFe is integral to Fe cycling and hence plays a role in regulating carbon cycling, warranting its’ inclusion in Fe budgets and biogeochemical models.

Description: Human-induced ocean warming and acidification have received increasing attention over the past decade and are considered to have substantial consequences for a broad range of marine species and their interactions. Understanding how these interactions shift in response to climate change is particularly important with regard to foundation species, such as the brown alga Fucus vesiculosus. This macroalga represents the dominant habitat former on coastal rocky substrata of the Baltic Sea, fulfilling functions essential for the entire benthic community. Its ability to withstand extensive fouling and herbivory regulates the associated community and ecosystem dynamics. This study tested the interactive effects of future warming, acidification, and seasonality on the interactions of a marine macroalga with potential foulers and consumers. F. vesiculosus rockweeds were exposed to different combinations of conditions predicted regionally for the year 2100 (+∆5°C, +∆700 μatm CO2) using multifactorial long-term experiments in novel outdoor benthic mesocosms (“Benthocosms”) over 9–12-week periods in four seasons. Possible shifts in the macroalgal susceptibility to fouling and consumption were tested using consecutive bioassays. Algal susceptibility to fouling and grazing varied substantially among seasons and between treatments. In all seasons, warming predominantly affected anti-fouling and anti-herbivory interactions while acidification had a subtle nonsignificant influence. Interestingly, anti-microfouling activity was highest during winter under warming, while anti-macrofouling and anti-herbivory activities were highest in the summer under warming. These contrasting findings indicate that seasonal changes in anti-fouling and anti-herbivory traits may interact with ocean warming in altering F. vesiculosus community composition in the future.

Description: Combined seawater radiogenic hafnium (Hf) and neodymium (Nd) isotope compositions were extracted from bulk sediment leachates and foraminifera of Site 1088, ODP Leg 177, 2082 m water depth on the Agulhas Ridge. The new data provide a continuous reconstruction of long and short-term changes in ocean circulation and continental weathering inputs since the Mid-Miocene. Due to its intermediate water depth the sediments of this core sensitively recorded changes in admixture of North Atlantic Deep Water (NADW) to the Antarctic Circumpolar Current (ACC) as a function of the strength of the Atlantic Meridional Overturning Circulation (AMOC). Nd isotope compositions (εNd) range from -7 to -11 with glacial values generally 1 to 3 units more radiogenic than during the interglacials of the Quaternary. The data reveal episodes of significantly increased AMOC strength during late Miocene and Pliocene warm periods whereas peak radiogenic εNd values mark a strongly diminished AMOC during the major intensification of Northern Hemisphere Glaciation near 2.8 Ma and in the Pleistocene after 1.5 Ma. In contrast, the Hf isotope compositions (εHf) show an essentially continuous evolution from highly radiogenic values of up to +11 during the Miocene to less radiogenic present day values (+2 to +4) during the late Quaternary. The data document a long-term transition in dominant weathering inputs, where inputs from the South America are replaced by those from Southern Africa. Moreover, radiogenic peaks provide evidence for the supply of radiogenic Hf originating from Patagonian rocks to the Atlantic sector of the Southern Ocean via dust inputs.

Description: While secondary contact between Mytilus edulis and Mytilus trossulus in North America results in mosaic hybrid zone formation, both species form a hybrid swarm in the Baltic. Despite pervasive gene flow, Baltic Mytilus species maintain substantial genetic and phenotypic differentiation. Exploring mechanisms underlying the contrasting genetic composition in Baltic Mytilus species will allow insights into processes such as speciation or adaptation to extremely low salinity. Previous studies in the Baltic indicated that only weak interspecific reproductive barriers exist and discussed the putative role of adaptation to environmental conditions. Using a combination of hydrodynamic modelling and multilocus genotyping, we investigate how oceanographic conditions influence passive larval dispersal and hybrid swarm formation in the Baltic. By combining our analyses with previous knowledge, we show a genetic transition of Baltic Mytilus species along longitude 12°-13°E, that is a virtual line between Malmö (Sweden) and Stralsund (Germany). Although larval transport only occurs over short distances (10–30 km), limited larval dispersal could not explain the position of this genetic transition zone. Instead, the genetic transition zone is located at the area of maximum salinity change (15–10 psu). Thus, we argue that selection results in weak reproductive barriers and local adaptation. This scenario could maintain genetic and phenotypic differences between Baltic Mytilus species despite pervasive introgressive hybridization.

Description: Habitat stratification by abiotic and biotic factors initiates divergence of populations and leads to ecological speciation. In contrast to fully marine waters, the Baltic Sea is stratified by a salinity gradient that strongly affects fish physiology, distribution, diversity and virulence of important marine pathogens. Animals thus face the challenge to simultaneously adapt to the concurrent salinity and cope with the selection imposed by the changing pathogenic virulence. Western Baltic spring-spawning herring (Clupea harengus) migrate to spawning grounds characterized by different salinities to which herring are supposedly adapted. We hypothesized that herring populations do not only have to cope with different salinity levels but that they are simultaneously exposed to higher-order effects that accompany the shifts in salinity, that is induced pathogenicity of Vibrio bacteria in lower saline waters. To experimentally evaluate this, adults of two populations were caught in their spawning grounds and fully reciprocally crossed within and between populations. Larvae were reared at three salinity levels, representing the spawning ground salinity of each of the two populations, or Atlantic salinity conditions resembling the phylogenetic origin of Clupea harengus. In addition, larvae were exposed to a Vibrio spp. infection. Life-history traits and gene expression analysis served as response variables. Herring seem adapted to Baltic Sea conditions and cope better with low saline waters. However, upon a bacterial infection, herring larvae suffer more when kept at lower salinities implying reduced resistance against Vibrio or higher Vibrio virulence. In the context of recent climate change with less saline marine waters in the Baltic Sea, such interactions may constitute key future stressors.

Description: Reconstructed sea surface temperatures (SSTs) derived from Mg/Ca measurements in nine encrusting coralline algal skeletons from the Aleutian archipelago in the northernmost Pacific Ocean reveal an overall increase in SST from 1665 to 2007. In the Aleutian SST reconstruction, decadal-scale variability is a transient feature present during the 1700s and early 1800s and then fully emerging post-1950. SSTs vary coherently with available instrument records of cyclone variance and vacillate in and out of coherence with multicentennial Pacific Northwest drought reconstructions as a response to SST-driven alterations of storm tracks reaching North America. These results indicate that an influence of decadal-scale variability on the North Pacific storm tracks only became apparent during the midtwentieth century. Furthermore, what has been assumed as natural variability in the North Pacific, based on twentieth century instrumental data, is not consistent with the long-term natural variability evident in reconstructed SSTs predating the anthropogenic influence.

Description: Microbathymetry data, in situ observations, and sampling along the 138200N and 138200N oceanic core complexes (OCCs) reveal mechanisms of detachment fault denudation at the seafloor, links between tectonic extension and mass wasting, and expose the nature of corrugations, ubiquitous at OCCs. In the initial stages of detachment faulting and high-angle fault, scarps show extensive mass wasting that reduces their slope. Flexural rotation further lowers scarp slope, hinders mass wasting, resulting in morphologically complex chaotic terrain between the breakaway and the denuded corrugated surface. Extension and drag along the fault plane uplifts a wedge of hangingwall material (apron). The detachment surface emerges along a continuous moat that sheds rocks and covers it with unconsolidated rubble, while local slumping emplaces rubble ridges overlying corrugations. The detachment fault zone is a set of anostomosed slip planes, elongated in the alongextension direction. Slip planes bind fault rock bodies defining the corrugations observed in microbathymetry and sonar. Fault planes with extension-parallel stria are exposed along corrugation flanks, where the rubble cover is shed. Detachment fault rocks are primarily basalt fault breccia at 138200N OCC, and gabbro and peridotite at 138300N, demonstrating that brittle strain localization in shallow lithosphere form corrugations, regardless of lithologies in the detachment zone. Finally, faulting and volcanism dismember the 138300N OCC, with widespread present and past hydrothermal activity (Semenov fields), while the Irinovskoe hydrothermal field at the 138200N core complex suggests a magmatic source within the footwall. These results confirm the ubiquitous relationship between hydrothermal activity and oceanic detachment formation and evolution.

Description: We present the first subprecessional record of seawater 87Sr/86Sr isotope ratios for a marginal Mediterranean subbasin. The sediments contained in this interval (three precessional cycles between 6.60 and 6.55 Ma) are important because they record conditions during the transition to the Messinian Salinity Crisis (MSC; 5.97 to 5.33 Ma), an event for which many details are still poorly understood. The record, derived from planktic foraminifera of the late Miocene Sorbas Basin (SE Spain), shows brief excursions with precessional cyclicity to 87Sr/86Sr ratios higher than coeval ocean 87Sr/86Sr. The hydrologic conditions required to generate the observed record are investigated using box modeling, constrained using a new paleodepth estimate (150 to 250 m) based on benthic foraminiferal assemblages. The box model results highlight the role of climate-driven interbasin density contrast as a significant driver of, or impediment to, exchange. The results are particularly significant in the context of the MSC, where 87Sr/86Sr excursions have been interpreted purely as a consequence of physical restriction. To replicate the observed temporal patterns of lithological variations and 87Sr/86Sr isotope excursions, the Sorbas Basin “box” must have a mainly positive hydrologic budget, in contrast with the Mediterranean's negative budget during the late Miocene. This result has implications for the assumption of synchronous deposition of specific sedimentary layers (sapropels) between marginal and open Mediterranean settings at subprecessional resolution. A net positive hydrologic budget in marginal Mediterranean subbasins may reconcile observations of freshwater inclusions in gypsum deposits.

Description: We report the results of a two-dimensional tomographic inversion of marine seismic refraction data from an array of ocean-bottom seismographs (OBSs), which produced an image of the crustal structure along the axial valley of the ultraslow spreading Mid-Cayman Spreading Center (MCSC). The seismic velocity model shows variations in the thickness and properties of the young oceanic crust that are consistent with the existence of two magmatic-tectonic segments along the 110 km long spreading center. Seismic wave speeds are consistent with exhumed mantle at the boundary between these two segments, but changes in the vertical gradient of seismic velocity suggest that volcanic crust occupies most of the axial valley seafloor along the seismic transect. The two spreading segments both have a low-velocity zone (LVZ) several kilometers beneath the seafloor, which may indicate the presence of shallow melt. However, the northern segment also has low seismic velocities (3 km/s) in a thick upper crustal layer (1.5–2.0 km), which we interpret as an extrusive volcanic section with high porosity and permeability. This segment hosts the Beebe vent field, the deepest known high-temperature black smoker hydrothermal vent system. In contrast, the southern spreading segment has seismic velocities as high as 4.0 km/s near the seafloor. We suggest that the porosity and permeability of the volcanic crust in the southern segment are much lower, thus limiting deep seawater penetration and hydrothermal recharge. This may explain why no hydrothermal vent system has been found in the southern half of the MCSC.

Description: Numerical Earth System Models are generic tools used to extrapolate present climate conditions into a warming future and to explore geoengineering options. Most of the current-generation models feature a simple pelagic biogeochemical model component that is embedded into a three-dimensional ocean general circulation model. The dynamics of these biogeochemical model components is essentially controlled by so-called model parameters most of which are poorly known. Here we explore the feasibility to estimate these parameters in a full-fledged three-dimensional Earth System Model by minimizing the misfit to noisy observations. The focus is on parameter identifiability. Based on earlier studies, we illustrate problems in determining a unique estimate of those parameters that prescribe the limiting effect of nutrient- and light-depleted conditions on carbon assimilation by autotrophic phytoplankton. Our results showcase that for typical models and evaluation metrics no meaningful “best” unique parameter set exists. We find very different parameter sets which are, on the one hand, equally consistent with our (synthetic) historical observations while, on the other hand, they propose strikingly differing projections into a warming climate.

Description: We present a comprehensive study showing new results from a shallow gas seep area in approximate to 40 m water depth located in the North Sea, Netherlands sector B13 that we call Dutch Dogger Bank seep area. It has been postulated that methane presumably originating from a gas reservoir in approximate to 600 m depth below the seafloor is naturally leaking to the seafloor. Our ship-based subbottom echosounder data indicate that the migrating gas is trapped in numerous gas pockets in the shallow sediments. The gas pockets are located at the boundary between the top of the Late Pliocene section and overlying fine-grained sediments, which were deposited during the early Holocene marine transgression after the last glaciation. We mapped gas emissions during three R/V Heincke cruises in 2014, 2015, and 2016 and repeatedly observed up to 850 flares in the study area. Most of them (approximate to 80%) were concentrated at five flare clusters. Our repeated analysis revealed spatial similarities of seep clusters, but also heterogeneities in emission intensities. A first calculation of the methane released from these clusters into the water column revealed a flow rate of 277 L/min (SD=140), with two clusters emitting 132 and 142 L/min representing the most significant seepage sites. Above these two flare clusters, elevated methane concentrations were recorded in atmospheric measurements. Our results illustrate the effective transport of methane via gas bubbles through a approximate to 40 m water column, and furthermore provide an estimate of the emission rate needed to allow for a contribution to the atmospheric methane concentration.

Description: Seafloor massive sulfide (SMS) deposits form on and just below the seafloor along submarine tectonic plate boundaries. The deposits form from seawater that circulates through the underlying crust, is heated, leaches metals and sulfur from the surrounding rock, and then ascends and vents at the seafloor, forming sulfide mineral accumulations rich in Cu, Zn, Pb, Au, and Ag. Hydrothermal circulation through the crust is driven by shallow magmatic heat sources along the plate boundaries. Although high temperature “black smoker” chimneys and the unique ecosystems that they support are the most recognizable features of these vent sites, the mineral deposits can take on a variety of forms, from individual chimneys of less than a meter tall to large mounds with diameters of several hundred meters. The description of the deposits as “massive” refers to the high proportion (typically over 60%) of sulfide minerals that make up the deposits. Other minerals that commonly occur in SMS deposits are sulfates (barite and anhydrite), amorphous silica, and clay minerals. At the time of writing, more than 500 sites of high temperature seafloor hydrothermal systems and related mineral deposits have been found of the seafloor.

Description: We use seismic oceanography to document and analyze oceanic thermohaline finestructure across the Tyrrhenian Sea. Multichannel seismic (MCS) reflection data were acquired during the MEDiterranean OCcidental survey in April-May 2010. We deployed along-track expendable bathythermograph probes simultaneous with MCS acquisition. At nearby locations we gathered conductivity-temperature-depth data. An autonomous glider survey added in-situ measurements of oceanic properties. The seismic reflectivity clearly delineates thermohaline finestructure in the upper 2,000 m of the water column, indicating the interfaces between Atlantic Water/Winter Intermediate Water, Levantine Intermediate Water, and Tyrrhenian Deep Water. We observe the Northern Tyrrhenian Anticyclone, a near-surface meso-scale eddy, plus laterally and vertically extensive thermohaline staircases. Using MCS we are able to fully image the anticyclone to a depth of 800 m and to confirm the horizontal continuity of the thermohaline staircases of more than 200 km. The staircases show the clearest step-like gradients in the center of the basin while they become more diffuse towards the periphery and bottom, where impedance gradients become too small to be detected by MCS. We quantify the internal wave field and find it to be weak in the region of the eddy and in the center of the staircases, while it is stronger near the coastlines. Our results indicate this is because of the influence of the boundary currents, which disrupt the formation of staircases by preventing diffusive convection. In the interior of the basin the staircases are clearer and the internal wave field weaker, suggesting that other mixing processes such as double-diffusion prevail. Synopsis We studied the internal temperature and salinity structure of the Tyrrhenian Sea (Mediterranean) using the multichannel seismic reflection method (the same used in the hydrocarbon industry). Low frequency sound (seismic) waves are produced at the surface with an explosive air source and recorded by a towed cable containing hydrophones (underwater microphones). The data are processed to reveal 'stratigraphy' that result from contrasts in density that are themselves caused by changes in temperature and salinity. In this way we can map ocean circulation in two-dimensions. We also deployed in situ oceanographic probes to measure temperature and salinity in order to corroborate and optimize the processing of the seismic data. We then quantified the internal gravity wave field by tracking the peaks of seismic trace wavelets. Our results show that the interior of the Tyrrhenian Sea is largely isolated from internal waves that are generated by a large cyclonic boundary current that contains waters from the Atlantic ocean and other parts of the Mediterranean. This isolation allows the thermohaline finestructure to form, where small scale vertical mixing processes are at play. Understanding these mixing processes will aid researchers study global ocean circulation and to add constraints that can help improve climate models.

Description: Based on an unprecedented dissolved barium (D_Ba) data set collected in the Mediterranean Sea during a zonal transect between the Lebanon coast and Gibraltar (M84/3 cruise, April 2011), we decompose the D_Ba distribution to isolate the contribution of biogeochemical processes from the impact of the oceanic circulation. We have built a simple parametric water mass analysis (Parametric Optimum Multiparameter analysis) to reconstruct the contribution of the different Mediterranean water masses to the thermohaline structure. These water mass fractions have then been used to successfully reconstruct the background vertical gradient of D_Ba reflecting the balance between the large-scale oceanic circulation and the biological activity over long time scales. Superimposed on the background field, several D_Ba anomalies have been identified. Positive anomalies are associated with topographic obstacles and may be explained by the dissolution of particulate biogenic barium (P_Ba barite) of material resuspended by the local currents. The derived dissolution rates range from 0.06 to 0.21 μmol m−2 d−1. Negative anomalies are present in the mesopelagic region of the western and eastern basins (except in the easternmost Levantine basin) as well as in the abyssal western basin. This represents the first quantification of the nonconservative component of the D_Ba signal. These mesopelagic anomalies could reflect the subtraction of D_Ba during P_Ba barite formation occurring during organic carbon remineralization. The deep anomalies may potentially reflect the transport of material toward the deep sea during winter deep convection and the subsequent remineralization. The D_Ba subtraction fluxes range from −0.07 to −1.28 μmol m−2 d−1. D_Ba-derived fluxes of P_Ba barite (up to 0.21 μmol m−2 d−1) and organic carbon (13 to 29 mmol C m−2 d−1) are in good agreement with other independent measurements suggesting that D_Ba can help constrain remineralization horizons. This study highlights the importance of quantifying the impact of the large-scale oceanic circulation in order to better understand the biogeochemical cycling of elements and to build reliable geochemical proxies.

Description: We investigated the onset and development of Cretaceous Oceanic Anoxic Event 2 (OAE2) in a newly drilled core (SN degrees 4) from the Tarfaya Basin (southern Morocco), where this interval is unusually expanded. High-resolution (centimeter-scale equivalent to centennial) analysis of bulk organic and carbonate stable isotopes and of carbonate and organic carbon content in combination with XRF scanner derived elemental distribution reveal that the ocean-climate system behaved in a highly dynamic manner prior to and during the onset of OAE2. Correlation with the latest orbital solution indicates that the main carbon isotope shift occurred during an extended minimum in orbital eccentricity (similar to 400 kyr cycle). Shorter-term fluctuations in carbonate and organic carbon accumulation and in sea level related terrigenous discharge were predominantly driven by variations in orbital obliquity. Negative excursions in organic and carbonate delta C-13 preceded the global positive delta C-13 shift marking the onset of OAE2, suggesting injection of isotopically depleted carbon into the atmosphere. The main delta C-13 increase during the early phase of OAE2 in the late Cenomanian was punctuated by a transient plateau. Maximum organic carbon accumulation occurred during the later part of the main delta C-13 increase and was associated with climate cooling events, expressed as three consecutive maxima in bulk carbonate delta O-18. The extinctions of the thermocline dwelling keeled planktonic foraminifers Rotalipora greenhornensis and Rotalipora cushmani occurred during the first and last of these cooling events and were likely associated with obliquity paced, ocean-wide expansions, and intensifications of the oxygen minimum zone, affecting their habitat space on a global scale.

Description: The potential of Coastal Ocean Alkalinization (COA), a carbon dioxide removal (CDR) climate engineering strategy that chemically increases ocean carbon uptake and storage, is investigated with an Earth system model of intermediate complexity. The CDR potential and possible environmental side effects are estimated for various COA deployment scenarios, assuming olivine as the alkalinity source in ice-free coastal waters (about 8.6% of the global ocean's surface area), with dissolution rates being a function of grain size, ambient seawater temperature and pH. Our results indicate that for a large-enough olivine deployment of small-enough grain sizes (10 μm), atmospheric CO2 could be reduced by more than 800 GtC by the year 2100. However, COA with coarse olivine grains (1000 μm) has little CO2 sequestration potential on this time scale. Ambitious CDR with fine olivine grains would increase coastal aragonite saturation Ω to levels well beyond those that are currently observed. When imposing upper limits for aragonite saturation levels (Ωlim) in the grid boxes subject to COA (Ωlim = 3.4 and 9 chosen as examples), COA still has the potential to reduce atmospheric CO2 by 265 GtC (Ωlim=3.4) to 790 GtC (Ωlim=9) and increase ocean carbon storage by 290 Gt (Ωlim=3.4) to 913 Gt (Ωlim=9) by year 2100.

Description: Transect surveys of hamlet communities (Hypoplectrus spp., Serranidae) covering 14 000 m2 across 16 reefs off La Parguera, Puerto Rico, are presented and compared with a previous survey conducted in the year 2000. The hamlet community has noticeably changed over 17 years, with a 〉 30% increase in relative abundance of the yellowtail hamlet Hypoplectrus chlorurus on the inner reefs at the expense of the other hamlet species. The data also suggest that the density of H. chlorurus has declined and that its distribution has shifted towards shallower depths. Considering that H. chlorurus has been previously identified as one of the few fish showing a positive association with seawater turbidity on the inner reefs of La Parguera and that sedimentation of terrestrial origin has increased over recent decades on these reefs, it is proposed that turbidity may constitute an important but so far overlooked ecological driver of hamlet communities.

Description: Methane hydrate-bearing sediments exist throughout the world in continental margins and in Arctic permafrost. Hydrates are ice-like compounds when dissociate due to temperature rise or reduction in fluid pressure, release gas. Because of the mechanical property changes caused by dissociation in which the loads supported by the hydrates are transferred to soil grains, these sediments may become unstable. To quantify the risk of ground instability triggered by dissociation, which may happen during operation to extract methane gas or from climate changes, a reliable predictive model is indispensable. Even though many models have been proposed, a detailed validation of the ability to model dissociation impact is still needed. This study investigated the adequacy of an spatially mobilized plane constitutive model and a modeling framework using laboratory-induced dissociation tests under shear from literature. Using laboratoryimposed temperature and pressure changes and the resulting hydrate saturation changes as input, this study was able to capture the geomechanical responses and determine the stability state of methane hydrate-bearing sediments as observed.

Description: This paper investigates new observations from the poorly understood region between the Kara and Laptev Seas in the Eastern Arctic Ocean. We discuss relevant circulation features including riverine freshwater, Atlantic-derived water, and polynya-formed dense water, emphasize Vilkitsky Strait (VS) as an important Kara Sea gateway, and analyze the role of the adjacent ∼250 km-long submarine Vilkitsky Trough (VT) for the Arctic boundary current. Expeditions in 2013 and 2014 operated closely spaced hydrographic transects and 1 year-long oceanographic mooring near VT's southern slope, and found persistent annually averaged flow of 0.2 m s−1 toward the Nansen Basin. The flow is nearly barotropic from winter through early summer and becomes surface intensified with maximum velocities of 0.35 m s−1 from August to October. Thermal wind shear is maximal above the southern flank at ∼30 m depth, in agreement with basinward flow above VT's southern slope. The subsurface features a steep front separating warm (–0.5°C) Atlantic-derived waters in central VT from cold (〈–1.5°C) shelf waters, which episodically migrates across the trough indicated by current reversals and temperature fluctuations. Shelf-transformed waters dominate above VT's slope, measuring near-freezing temperatures throughout the water column at salinities of 34–35. These dense waters are vigorously advected toward the Eurasian Basin and characterize VT as a conduit for near-freezing waters that could potentially supply the Arctic Ocean's lower halocline, cool Atlantic water, and ventilate the deeper Arctic Ocean. Our observations from the northwest Laptev Sea highlight a topographically complex region with swift currents, several water masses, narrow fronts, polynyas, and topographically channeled storms.

Description: Mixotrophs are important components of the bacterioplankton, phytoplankton, microzooplankton, and (sometimes) zooplankton in coastal and oceanic waters. Bacterivory among the phytoplankton may be important for alleviating inorganic nutrient stress and may increase primary production in oligotrophic waters. Mixotrophic phytoflagellates and dinoflagellates are often dominant components of the plankton during seasonal stratification. Many of the microzooplankton grazers, including ciliates and Rhizaria, are mixotrophic owing to their retention of functional algal organelles or maintenance of algal endosymbionts. Phototrophy among the microzooplankton may increase gross growth efficiency and carbon transfer through the microzooplankton to higher trophic levels. Characteristic assemblages of mixotrophs are associated with warm, temperate, and cold seas and with stratification, fronts, and upwelling zones. Modeling has indicated that mixotrophy has a profound impact on marine planktonic ecosystems and may enhance primary production, biomass transfer to higher trophic levels, and the functioning of the biological carbon pump.

Description: Aim Following the biogeographical approach implemented by Longhurst for the epipelagic layer, we propose here to identify a biogeochemical 3-D partition for the mesopelagic layer. The resulting partition characterizes the main deep environmental biotopes and their vertical boundaries on a global scale, which can be used as a geographical and ecological framework for conservation biology, ecosystem-based management and for the design of oceanographic investigations. Location The global ocean. Methods Based on the most comprehensive environmental climatology available to date, which is both spatially and vertically resolved (seven environmental parameters), we applied a combination of clustering algorithms (c-means, k-means, partition around medoids and agglomerative with Ward's linkage) associated with a nonparametric environmental model to identify the vertical and spatial delineation of the mesopelagic layer. Results First, we show via numerical interpretation that the vertical division of the pelagic zone varies and, hence, is not constant throughout the global ocean. Indeed, a latitudinal gradient is found between the epipelagic–mesopelagic and mesopelagic–bathypelagic vertical limits. Second, the mesopelagic layer is shown here to be composed of 13 distinguishable Biogeochemical Provinces. Each province shows a distinct range of environmental conditions and characteristic 3-D distributions. Main conclusions The historical definition of the mesopelagic zone is here revisited to define a 3-D geographical framework and characterize all the deep environmental biotopes of the deep global ocean. According to the numerical interpretation of mesopelagic boundaries, we reveal that the vertical division of the zone is not constant over the global ocean (200–1,000 m) but varies between ocean basin and with latitude. We also provide evidence of biogeochemical division of the mesopelagic zone that is spatially structured in a similar way than the epipelagic in the shallow waters but varies in the deep owing to a change of the environmental driving factors.

Description: The nature and origin of the J-magnetic anomaly along the Iberia–Newfoundland margins are controversial and its validity for plate kinematic reconstructions questioned. At present, it is interpreted as either an oceanic isochron or an edge effect of oceanic crust corresponding to lithosphere breakup. Both interpretations result in restorations that are in conflict with the current knowledge from Pyrenean and North Atlantic geology. We combine seismic interpretations and dating of magmatic additions with magnetic data to examine the nature and formation process of this anomaly and discuss its value for plate restorations. We show that the J-anomaly is the result of polygenic and multiple magmatic events occurring during and after the formation of the first oceanic crust. Therefore, we conclude that the J-anomaly cannot be used for plate kinematic studies and, more generally, we question the validity of using ill-defined magnetic anomalies outside unequivocal oceanic domains for plate reconstructions.

Description: We report a new synthesis of best estimates of the inputs of fixed nitrogen to the world ocean via atmospheric deposition and compare this to fluvial inputs and dinitrogen fixation. We evaluate the scale of human perturbation of these fluxes. Fluvial inputs dominate inputs to the continental shelf, and we estimate that about 75% of this fluvial nitrogen escapes from the shelf to the open ocean. Biological dinitrogen fixation is the main external source of nitrogen to the open ocean, i.e., beyond the continental shelf. Atmospheric deposition is the primary mechanism by which land-based nitrogen inputs, and hence human perturbations of the nitrogen cycle, reach the open ocean. We estimate that anthropogenic inputs are currently leading to an increase in overall ocean carbon sequestration of ~0.4% (equivalent to an uptake of 0.15 Pg C yr−1 and less than the Duce et al. (2008) estimate). The resulting reduction in climate change forcing from this ocean CO2 uptake is offset to a small extent by an increase in ocean N2O emissions. We identify four important feedbacks in the ocean atmosphere nitrogen system that need to be better quantified to improve our understanding of the perturbation of ocean biogeochemistry by atmospheric nitrogen inputs. These feedbacks are recycling of (1) ammonia and (2) organic nitrogen from the ocean to the atmosphere and back, (3) the suppression of nitrogen fixation by increased nitrogen concentrations in surface waters from atmospheric deposition, and (4) increased loss of nitrogen from the ocean by denitrification due to increased productivity stimulated by atmospheric inputs.

Description: The goal of this study is to computationally determine the potential distribution patterns of diffusion-driven methane hydrate accumulations in coarse-grained marine sediments. Diffusion of dissolved methane in marine gas hydrate systems has been proposed as a potential transport mechanism through which large concentrations of hydrate can preferentially accumulate in coarse-grained sediments over geologic time. Using one-dimensional compositional reservoir simulations, we examine hydrate distribution patterns at the scale of individual sand layers (1-20 m thick) that are deposited between microbially active fine-grained material buried through the gas hydrate stability zone (GHSZ). We then extrapolate to two-dimensional and basin-scale three-dimensional simulations, where we model dipping sands and multilayered systems. We find that properties of a sand layer including pore size distribution, layer thickness, dip, and proximity to other layers in multilayered systems all exert control on diffusive methane fluxes toward and within a sand, which in turn impact the distribution of hydrate throughout a sand unit. In all of these simulations, we incorporate data on physical properties and sand layer geometries from the Terrebonne Basin gas hydrate system in the Gulf of Mexico. We demonstrate that diffusion can generate high hydrate saturations (upward of 90%) at the edges of thin sands at shallow depths within the GHSZ, but that it is ineffective at producing high hydrate saturations throughout thick (greater than 10 m) sands buried deep within the GHSZ. Furthermore, we find that hydrate in fine-grained material can preserve high hydrate saturations in nearby thin sands with burial.

Description: The Norwegian Channel Ice Stream of the Scandinavian Ice Sheet extended across the northern North Sea margin during the mid- to late Quaternary, eroding older sediment from the continental shelf. Consequently, little is known about early Quaternary sedimentation on this margin. We use two- and threedimensional seismic-reflection data to investigate changing sediment volumes and sources in the northern North Sea through the Quaternary. The northern North Sea Basin was infilled during the early Quaternary by intercalated glacigenic debris-flows and contourites, which provide a record of the delivery of glacigenic sediment to the slope and the intensity of North Atlantic thermohaline circulation during early Quaternary glacial ± interglacial cycles. The infilling of the basin reduced accommodation and led to the deflection of mid- to late Quaternary sediments into the Norwegian Sea, forming the North Sea Fan. Close to the onset of the mid-Quaternary, the south-western Scandinavian Ice Sheet margin was drained by an ice stream located beneath Målùy Plateau, 60 km east of the Last Glacial Maximum Norwegian Channel Ice Stream. The southward-flowing Norwegian Sea Bottom Water current was directed into the partially filled northern North Sea Basin during the early Quaternary, and deflected progressively northwards as the basin became infilled.

Description: The impact of a subgrid-scale ice thickness distribution (ITD) and two standard ice strength formulations on simulated Arctic sea ice climate is investigated. To this end, different model configurations with and without an ITD were tuned by minimizing the weighted mean error between the simulated and observed sea ice concentration, thickness, and drift speed with an semiautomatic parameter optimization routine. The standard ITD and ice strength parameterization lead to larger errors when compared to the simple single-category model with an ice strength parameterization based on the mean ice thickness. Interestingly, the simpler ice strength formulation, which depends linearly on the mean ice thickness, also reduces the model-observation error when using an ITD. For the ice strength parameterization that makes use of the ITD, the effective ice strength depends strongly on the number of thickness categories, so that introducing more categories can lead to overall thicker ice that is more easily deformed.

Description: Dissolved rare earth element (REE) concentration data from intermediate and deep seawater form an array characterized by higher middle REE enrichments (MREE/MREE*) in the North Atlantic and a progressive increase in heavy-to-light REE ratios (HREE/LREE) as water masses age. The REEs in foraminifera are fractionated towards higher MREE/MREE* and lower HREE/LREE relative to seawater. Calculations based on a scavenging model show that the REE patterns in uncleaned core-top foraminifera resemble those adsorbed onto calcite, particulate organic material, and hydrous ferric oxides but the full extent of the REE fractionation measured in foraminifera was not reproduced by the model. However, differences in the HREE/LREE, MREE/MREE* ratios and the cerium anomaly between ocean basins are preserved and are in agreement with the seawater REE distribution. Under oxic conditions, the HREE/LREE and MREE/MREE* compositions of uncleaned foraminifera at the sediment/seawater boundary are preserved during burial but the cerium anomaly is sensitive to burial depth. In suboxic sedimentary environments, all uncleaned foraminiferal REE concentrations are elevated relative to core-top values indicating addition of REEs from pore waters. The HREE/LREE ratio is highest when sedimentation rates were greatest [Lippold et al., 2009] and when high Fe/Ca ratios in the uncleaned foraminifera indicate that Fe was mobile. In sediments that have not experienced suboxic conditions during burial, uncleaned foraminifera preserve the seawater signal taken up at the sediment/seawater interface and are therefore suggested to be a suitable archive of changes in the REE signal of past bottom waters. This article is protected by copyright. All rights reserved.

Description: First Workshop of the Carbon Dioxide Removal Model Intercomparison Project; Potsdam, Germany, 20–22 September 2016

Description: Our study presents a basin-scale 3D modeling solution, quantifying and exploring gas hydrate accumulations in the marine environment around the Green Canyon (GC955) area, Gulf of Mexico. It is the first modeling study that considers the full complexity of gas hydrate formation in a natural geological system. Overall, it comprises a comprehensive basin re-construction, accounting for depositional and transient thermal history of the basin, source rock maturation, petroleum components generation, expulsion and migration, salt tectonics and associated multi-stage fault development. The resulting 3D gas hydrate distribution in the Green Canyon area is consistent with independent borehole observations. An important mechanism identified in this study and leading to high gas hydrate saturation (〉 80 vol. %) at the base of the gas hydrate stability zone (GHSZ), is the recycling of gas hydrate and free gas enhanced by high Neogene sedimentation rates in the region. Our model predicts the rapid development of secondary intra-salt mini-basins situated on top of the allochthonous salt deposits which leads to significant sediment subsidence and an ensuing dislocation of the lower GHSZ boundary. Consequently, large amounts of gas hydrates located in the deepest parts of the basin dissociate and the released free methane gas migrates upwards to recharge the GHSZ. In total, we have predicted the gas hydrate budget for the Green Canyon area that amounts to ∼3,256 Mt of gas hydrate which is equivalent to ∼340 Mt of carbon (∼7 x 1011 m3 of CH4 at STP conditions), and consists mostly of biogenic hydrates.

Description: The ocean is an important source of nitrous oxide (N2O) to the atmosphere, yet the factors controlling N2O production and consumption in oceanic environments are still not understood nor constrained. We measured N2O concentrations and isotopomer ratios, as well as O2, nutrient and biogenic N2 concentrations, and the isotopic compositions of nitrate and nitrite at several coastal stations during two cruises off the Peru coast (~5–16°S, 75–81°W) in December 2012 and January 2013. N2O concentrations varied from below equilibrium values in the oxygen deficient zone (ODZ) to up to 190 nmol L−1 in surface waters. We used a 3-D-reaction-advection-diffusion model to evaluate the rates and modes of N2O production in oxic waters and rates of N2O consumption versus production by denitrification in the ODZ. Intramolecular site preference in N2O isotopomer was relatively low in surface waters (generally −3 to 14‰) and together with modeling results, confirmed the dominance of nitrifier-denitrification or incomplete denitrifier-denitrification, corresponding to an efflux of up to 0.6 Tg N yr−1 off the Peru coast. Other evidence, e.g., the absence of a relationship between ΔN2O and apparent O2 utilization and significant relationships between nitrate, a substrate during denitrification, and N2O isotopes, suggest that N2O production by incomplete denitrification or nitrifier-denitrification decoupled from aerobic organic matter remineralization are likely pathways for extreme N2O accumulation in newly upwelled surface waters. We observed imbalances between N2O production and consumption in the ODZ, with the modeled proportion of N2O consumption relative to production generally increasing with biogenic N2. However, N2O production appeared to occur even where there was high N loss at the shallowest stations.

Description: The stable production of high-quality fry in marine aquaculture is still hampered by unpredictable mortality caused by infectious diseases during larval rearing. Consequently, the development of new biocontrol agents is crucial for a viable aquaculture industry. The bacterial energy storage compound poly-β-hydroxybutyrate (PHB) has been shown to exhibit beneficial properties on aquatic organisms such as enhanced survival, growth, disease resistance and a controlling effect on the gastrointestinal microbiota. However, the effect of PHB on the developing immune system of fish larvae has so far not been investigated. In this study, the effect of feeding PHB-enriched Artemia nauplii on survival, growth and immune response in European sea bass (Dicentrarchus labrax) postlarvae was examined. Amorphous PHB was administered to 28-day-old sea bass postlarvae over a period of 10 days. The survival and growth performance were monitored, and the expression of 29 genes involved in immunity, growth, metabolism and stress-response was measured. While the expression of the insulin-like growth factor 1 (igf1), an indicator of relative growth, was upregulated in response to feeding PHB, the larval survival and growth performance remained unaffected. After 10 days of PHB treatment, the expression of the antimicrobial peptides dicentracin (dic) and hepcidin (hep) as well as mhc class IIa and mhc class IIb was elevated in the PHB fed postlarvae. This indicates that PHB is capable of stimulating the immune system of fish early life stages, which may be the cause of the increased resistance to diseases and robustness observed in previous studies.

Description: Methane concentrations in the water column and emissions to the atmosphere were determined for three tropical coastal lagoons surrounded by mangrove forests on the Yucatán Peninsula, Mexico. Surface water dissolved methane was sampled at different seasons over a period of 2 years in areas representing a wide range of salinities and anthropogenic impacts. The highest surface water methane concentrations (up to 8378 nM) were measured in a polluted canal associated with Terminos Lagoon. In Chelem Lagoon, methane concentrations were typically lower, except in the polluted harbor area (1796 nM). In the relatively pristine Celestún Lagoon, surface water methane concentrations ranged from 41 to 2551 nM. Methane concentrations were negatively correlated with salinity in Celestún, while in Chelem and Terminos high methane concentrations were associated with areas of known pollution inputs, irrespective of salinity. The diffusive methane flux from surface lagoon water to the atmosphere ranged from 0.0023 to 15 mmol CH4 m-2 d-1. Flux chamber measurements revealed that direct methane release as ebullition was up to 3 orders of magnitude greater than measured diffusive flux. Coastal mangrove lagoons may therefore be an important natural source of methane to the atmosphere despite their relatively high salinity. Pollution inputs are likely to substantially enhance this flux. Additional statistically rigorous data collected globally are needed to better consider methane fluxes from mangrove-surrounded coastal areas in response to sea level changes and anthropogenic pollution in order to refine projections of future atmospheric methane budgets.

Description: Global climate models exhibit large biases in the Southern Ocean. For example, in models Antarctic bottom water is formed mostly through open-ocean deep-convection rather than through shelf convection. Still, the timescale, region, and intensity of deep-convection variability vary widely among models. We investigate the physical controls of this variability in the Atlantic sector of the Southern Ocean, where most of the models simulate recurring deep-convection events. We analyzed output from eleven exemplary CMIP5 models and four versions of the Kiel Climate Model (KCM). Of several potential physical control parameters that we tested, the ones shared by all these models are: Stratification in the convection region influences the timescale of the deep-convection variability, i.e. models with a strong (weak) stratification vary on long (short) timescales. And, sea ice volume affects the modeled mean state in the Southern Ocean: large (small) sea ice volume is associated with a non-convective (convective) predominant regime.

Description: The relationships between tectonic processes, magmatism, and hydrothermal venting along ∼600 km of the slow-spreading Mariana back-arc between 12.7°N and 18.3°N reveal a number of similarities and differences compared to slow-spreading mid-ocean ridges. Analysis of the volcanic geomorphology and structure highlights the complexity of the back-arc spreading center. Here, ridge segmentation is controlled by large-scale basement structures that appear to predate back-arc rifting. These structures also control the orientation of the chains of cross-arc volcanoes that characterize this region. Segment-scale faulting is oriented perpendicular to the spreading direction, allowing precise spreading directions to be determined. Four morphologically distinct segment types are identified: dominantly magmatic segments (Type I); magmatic segments currently undergoing tectonic extension (Type II); dominantly tectonic segments (Type III); and tectonic segments currently undergoing magmatic extension (Type IV). Variations in axial morphology (including eruption styles, neovolcanic eruption volumes, and faulting) reflect magma supply, which is locally enhanced by cross-arc volcanism associated with N-S compression along the 16.5°N and 17.0°N segments. In contrast, cross-arc seismicity is associated with N-S extension and increased faulting along the 14.5°N segment, with structures that are interpreted to be oceanic core complexes—the first with high-resolution bathymetry described in an active back-arc basin. Hydrothermal venting associated with recent magmatism has been discovered along all segment types.

Description: As land-based mineral resources become increasingly difficult and expensive to acquire, the potential for mining resources from the deep seafloor has become widely discussed and debated. Exploration leases are being granted, and technologies are under development. However, the quantity and quality of the resources are uncertain, and many worry about risks to vulnerable deep-sea ecosystems. Deep-sea mining has become part of the discussion of the United Nations Sustainable Development Goals. In this article we provide a summary of benefits, costs, and uncertainties that surround this potentially attractive but contentious topic.

Description: Previous studies have estimated that mantle serpentinization reactions generate H2 at a rate of 1010–1012 mol/yr along the global mid-ocean ridge (MOR) system. Here we present results of 3-D geodynamic simulations that predict rates of additional mantle serpentinization and H2 production at oceanic transform faults (OTF). We find that the extent and rate of mantle serpentinization increases with OTF length and is maximum at intermediate slip rates of 5 to 10 cm/yr. The additional global OTF-related production of H2 is found to be between 6.1 and 10.7 × 1011 mol/yr, which is comparable to the predicted background MOR rate of 4.1–15.0 × 1011 mol H2/yr. This points to oceanic transform faults as potential sites of intense fluid-rock interaction, where chemosynthetic life could be sustained by serpentinization reactions.

Description: A global inventory of data from gas hydrate drilling expeditions is used to develop relationships between the base of structure I gas hydrate stability, top of gas hydrate occurrence, sulfate-methane transition depth, pressure (water depth), and geothermal gradients. The motivation of this study is to provide first-order estimates of the top of gas hydrate occurrence and associated thickness of the gas hydrate occurrence zone for climate-change scenarios, global carbon budget analyses, or gas hydrate resource assessments. Results from publically available drilling campaigns (21 expeditions and 52 drill sites) off Cascadia, Blake Ridge, India, Korea, South China Sea, Japan, Chile, Peru, Costa Rica, Gulf of Mexico, and Borneo reveal a first-order linear relationship between the depth to the top and base of gas hydrate occurrence. The reason for these nearly linear relationships is believed to be the strong pressure and temperature dependence of methane solubility in the absence of large difference in thermal gradients between the various sites assessed. In addition, a statistically robust relationship was defined between the thickness of the gas hydrate occurrence zone and the base of gas hydrate stability (in meters below seafloor). The relationship developed is able to predict the depth of the top of gas hydrate occurrence zone using observed depths of the base of gas hydrate stability within less than 50 m at most locations examined in this study. No clear correlation of the depth to the top and base of gas hydrate occurrences with geothermal gradient and sulfate-methane transition depth was identified.

Description: The Vestnesa Ridge comprises a 〉100 km long sediment drift located between the western continental slope of Svalbard and the Arctic mid-ocean ridges. It hosts a deep water (〉1000 m) gas hydrate and associated seafloor seepage system. Near-seafloor headspace gas compositions and its methane carbon isotopic signature along the ridge indicate a predominance of thermogenic gas sources feeding the system. Prediction of the base of the gas hydrate stability zone for theoretical pressure and temperature conditions and measured gas compositions results in an unusual underestimation of the observed bottom-simulating reflector (BSR) depth. The BSR is up to 60 m deeper than predicted for pure methane and measured gas compositions with 〉99% methane. Models for measured gas compositions with 〉4% higher-order hydrocarbons result in a better BSR approximation. However, the BSR remains 〉20 m deeper than predicted in a region without active seepage. A BSR deeper than predicted is primarily explained by unaccounted spatial variations in the geothermal gradient and by larger amounts of thermogenic gas at the base of the gas hydrate stability zone. Hydrates containing higher-order hydrocarbons form at greater depths and higher temperatures and contribute with larger amounts of carbons than pure methane hydrates. In thermogenic provinces, this may imply a significant upward revision (up to 50% in the case of Vestnesa Ridge) of the amount of carbon in gas hydrates.

Description: This paper introduces the Distribution-Independent Storm Severity Index (DI-SSI). The DI-SSI represents an approach to quantify the severity of exceptional surface wind speeds of large scale windstorms that is complementary to the SSI introduced by Leckebusch et al. While the SSI approaches the extremeness of a storm from a meteorological and potential loss (impact) perspective, the DI-SSI defines the severity in a more climatological perspective. The idea is to assign equal index values to wind speeds of the same singularity (e.g. the 99th percentile) under consideration of the shape of the tail of the local wind speed climatology. Especially in regions at the edge of the classical storm track, the DI-SSI shows more equitable severity estimates, e.g. for the extra-tropical cyclone Klaus. In order to compare the indices, their relation with the North Atlantic Oscillation is studied, which is one of the main large scale drivers for the intensity of European windstorms.

Description: Many blueschists and eclogites are inferred to have formed from oceanic basalts in subducted slabs. Knowledge of their elastic behaviour is essential for reconstructing the internal structure of subduction zones. The Cycladic Blueschist Unit, exposed on Syros Island (Greece), contains rocks belonging to an exhumed Tertiary subduction complex. They were possibly part of a subduction channel, a shear zone above the subducting slab in which exhumation is possible during subduction. Intense plastic deformation, forming crystallographic preferred orientations (CPO), accompanied blueschist and eclogite metamorphism. CPO of the constituent minerals in the collected samples was determined by time-of-flight neutron diffraction. Two samples are foliated fine-grained blueschists with strong CPO, rich in glaucophane, zoisite and phengite. Two coarser-grained eclogite samples rich in omphacite and clinozoisite, or glaucophane, have weaker CPO. Vp and Vs anisotropies were computed from the orientation distribution function and single-crystal elastic constants. All samples show velocity maxima parallel to the mineral lineation, and minima normal to the foliation, providing important constraints on orientations of seismic anisotropy in subduction channels. Vp anisotropies are up to three times higher (6.5-12%) in the blueschists than in the eclogites (3-4%), pointing to a potentially important lithological control of elastic anisotropy in subducted oceanic crust.

Description: Global-scale gene flow is an important concern in conservation biology as it has the potential to either increase or decrease genetic diversity in species and populations. Although many studies focus on the gene flow between different populations of a single species, the potential for gene flow and introgression between species is understudied, particularly in seabirds. The only well-studied example of a mixed-species, hybridizing population of petrels exists on Round Island, in the Indian Ocean. Previous research assumed that Round Island represents a point of secondary contact between Atlantic (Pterodroma arminjoniana) and Pacific species (Pterodroma neglecta and Pterodroma heraldica). This study uses microsatellite genotyping and tracking data to address the possibility of between-species hybridization occurring outside the Indian Ocean. Dispersal and gene flow spanning three oceans were demonstrated between the species in this complex. Analysis of migration rates estimated using bayesass revealed unidirectional movement of petrels from the Atlantic and Pacific into the Indian Ocean. Conversely, structure analysis revealed gene flow between species of the Atlantic and Pacific oceans, with potential three-way hybrids occurring outside the Indian Ocean. Additionally, geolocation tracking of Round Island petrels revealed two individuals travelling to the Atlantic and Pacific. These results suggest that interspecific hybrids in Pterodroma petrels are more common than was previously assumed. This study is the first of its kind to investigate gene flow between populations of closely related Procellariiform species on a global scale, demonstrating the need for consideration of widespread migration and hybridization in the conservation of threatened seabirds.

Description: Measurements of near-surface oxygen (O2) concentrations and mixed layer depth from the K1 mooring in the central Labrador Sea are used to calculate the change in column-integrated (0–1700 m) O2 content over the deep convection winter 2014/2015. During the mixed layer deepening period, November 2014 to April 2015, the oxygen content increased by 24.3 ± 3.4 mol m−2, 40% higher than previous results from winters with weaker convection. By estimating the contribution of respiration and lateral transport on the oxygen budget, the cumulative air-sea gas exchange is derived. The O2 uptake of 29.1 ± 3.8 mol m−2, driven by persistent undersaturation (≥5%) and strong atmospheric forcing, is substantially higher than predicted by standard (nonbubble) gas exchange parameterizations, whereas most bubble-resolving parameterizations predict higher uptake, comparable to our results. Generally large but varying mixed layer depths and strong heat and momentum fluxes make the Labrador Sea an ideal test bed for process studies aimed at improving gas exchange parameterizations.

Description: The taxonomy of the ocean sunfishes (Molidae) has a complicated history. Currently, three genera and four species are recognized, including two in the genus Mola (M. mola and M. ramsayi). In 2009, a genetic study revealed a potential third species, Mola species C, in Southeast Australian waters. Concentrating on this region, we obtained samples and morphological data from 27 Mola sp. C specimens, genetically confirmed the existence of this species (mtDNA D-loop and cytochrome c oxidase 1), and established its morphology across a size spectrum of 50–242 cm total length. Mola sp. C is diagnosed by clavus meristics [15–17 fin rays (13–15 principal, 2 minor), 5–7 ossicles, paraxial ossicles separate], clavus morphology (prominent smooth band back-fold, rounded clavus edge with an indent), and body scale morphology (raised conical midpoints, non-branching). This species does not develop a protruding snout, or swollen dorso- or ventrolateral ridges. Body proportions remain similar with growth. A review of the historic literature revealed that Mola sp. C is a new, hitherto undescribed species, M. tecta, which we describe and diagnose, and that it is the first proposed addition to the genus Mola in 125 years. Its core distribution is likely in the temperate waters of the Southern Hemisphere.

Description: Disentangling the processes and mechanisms underlying adaptive diversification is facilitated by the comparative study of replicate population pairs that have diverged along a similar environmental gradient. Such a setting is realized in a cichlid fish from southern Lake Tanganyika, Astatotilapia burtoni, which occurs within the lake proper as well as in various affluent rivers. Previously, we demonstrated that independent lake and stream populations show similar adaptations to the two habitat regimes. However, little is known about the evolutionary and demographic history of the A. burtoni populations in question and the patterns of genome divergence among them. Here, we apply restriction site-associated DNA sequencing (RADseq) to examine the evolutionary history, the population structure and genomic differentiation of lake and stream populations in A. burtoni. A phylogenetic reconstruction based on genome-wide molecular data largely resolved the evolutionary relationships among populations, allowing us to re-evaluate the independence of replicate lake-stream population clusters. Further, we detected a strong pattern of isolation by distance, with baseline genomic divergence increasing with geographic distance and decreasing with the level of gene flow between lake and stream populations. Genome divergence patterns were heterogeneous and inconsistent among lake-stream population clusters, which is explained by differences in divergence times, levels of gene flow and local selection regimes. In line with the latter, we only detected consistent outlier loci when the most divergent lake-stream population pair was excluded. Several of the thus identified candidate genes have inferred functions in immune and neuronal systems and show differences in gene expression between lake and stream populations.

Description: Earthquake locations along the southern Mid-Atlantic Ridge have large uncertainties due to the sparse distribution of permanent seismological stations in and around the South Atlantic Ocean. Most of the earthquakes are associated with plate tectonic processes related to the formation of new oceanic lithosphere, as they are located close to the ridge axis or in the immediate vicinity of transform faults. A local seismological network of ocean-bottom seismometers and land stations on and around the archipelago of Tristan da Cunha, allowed for the first time a local earthquake survey for one year. We relate intra-plate seismicity within the African oceanic plate segment north of the island partly to extensional stresses induced by a bordering large transform fault and to the existence of the Tristan mantle plume. The temporal propagation of earthquakes within the segment reflects the prevailing stress field. The strong extensional stresses in addition with the plume weaken the lithosphere and might hint at an incipient ridge jump. An apparently aseismic zone coincides with the proposed location of the Tristan conduit in the upper mantle southwest of the islands. The margins of this zone describe the transition between the ductile and the surrounding brittle regime. Moreover, we observe seismicity close to the islands of Tristan da Cunha and nearby seamounts, which we relate to ongoing tectono-magmatic activity.

Description: Benthic–pelagic coupling is manifested as the exchange of energy, mass, or nutrients between benthic and pelagic habitats. It plays a prominent role in aquatic ecosystems, and it is crucial to functions from nutrient cycling to energy transfer in food webs. Coastal and estuarine ecosystem structure and function are strongly affected by anthropogenic pressures; however, there are large gaps in our understanding of the responses of inorganic nutrient and organic matter fluxes between benthic habitats and the water column. We illustrate the varied nature of physical and biological benthic–pelagic coupling processes and their potential sensitivity to three anthropogenic pressures – climate change, nutrient loading, and fishing – using the Baltic Sea as a case study and summarize current knowledge on the exchange of inorganic nutrients and organic material between habitats. Traditionally measured benthic–pelagic coupling processes (e.g., nutrient exchange and sedimentation of organic material) are to some extent quantifiable, but the magnitude and variability of biological processes are rarely assessed, preventing quantitative comparisons. Changing oxygen conditions will continue to have widespread effects on the processes that govern inorganic and organic matter exchange among habitats while climate change and nutrient load reductions may have large effects on organic matter sedimentation. Many biological processes (predation, bioturbation) are expected to be sensitive to anthropogenic drivers, but the outcomes for ecosystem function are largely unknown. We emphasize how improved empirical and experimental understanding of benthic–pelagic coupling processes and their variability are necessary to inform models that can quantify the feedbacks among processes and ecosystem responses to a changing world.

Description: Greenland's bed topography is a primary control on ice flow, grounding line migration, calving dynamics, and subglacial drainage. Moreover, fjord bathymetry regulates the penetration of warm Atlantic water (AW) that rapidly melts and undercuts Greenland's marine-terminating glaciers. Here we present a new compilation of Greenland bed topography that assimilates seafloor bathymetry and ice thickness data through a mass conservation approach. A new 150 m horizontal resolution bed topography/bathymetric map of Greenland is constructed with seamless transitions at the ice/ocean interface, yielding major improvements over previous data sets, particularly in the marine-terminating sectors of northwest and southeast Greenland. Our map reveals that the total sea level potential of the Greenland ice sheet is 7.42 ± 0.05 m, which is 7 cm greater than previous estimates. Furthermore, it explains recent calving front response of numerous outlet glaciers and reveals new pathways by which AW can access glaciers with marine-based basins, thereby highlighting sectors of Greenland that are most vulnerable to future oceanic forcing.

Description: This study presents a Monte Carlo method (CMSY) for estimating fisheries reference points from catch, resilience and qualitative stock status information on data-limited stocks. It also presents a Bayesian state-space implementation of the Schaefer production model (BSM), fitted to catch and biomass or catch-per-unit-of-effort (CPUE) data. Special emphasis was given to derive informative priors for productivity, unexploited stock size, catchability and biomass from population dynamics theory. Both models gave good predictions of the maximum intrinsic rate of population increase r, unexploited stock size k and maximum sustainable yield MSY when validated against simulated data with known parameter values. CMSY provided, in addition, reasonable predictions of relative biomass and exploitation rate. Both models were evaluated against 128 real stocks, where estimates of biomass were available from full stock assessments. BSM estimates of r, k and MSY were used as benchmarks for the respective CMSY estimates and were not significantly different in 76% of the stocks. A similar test against 28 data-limited stocks, where CPUE instead of biomass was available, showed that BSM and CMSY estimates of r, k and MSY were not significantly different in 89% of the stocks. Both CMSY and BSM combine the production model with a simple stock–recruitment model, accounting for reduced recruitment at severely depleted stock sizes.