ALBERT

Description: The distribution of the main water masses in the Atlantic Ocean are investigated with the Optimal Multi-Parameter (OMP) method. The properties of the main water masses in the Atlantic Ocean are described in a companion article; here these definitions are used to map out the general distribution of those water masses. Six key properties, including conservative (potential temperature and salinity) and non-conservative (oxygen, silicate, phosphate and nitrate), are incorporated into the OMP analysis to determine the contribution of the water masses in the Atlantic Ocean based on the GLODAP v2 observational data. To facilitate the analysis the Atlantic Ocean is divided into four vertical layers based on potential density. Due to the high seasonal variability in the mixed layer, this layer is excluded from the analysis. Central waters are the main water masses in the upper/central layer, generally featuring high potential temperature and salinity and low nutrient concentrations and are easily distinguished from the intermediate water masses. In the intermediate layer, the Antarctic Intermediate Water (AAIW) from the south can be detected to ~30°N, whereas the Subarctic Intermediate Water (SAIW), having similarly low salinity to the AAIW flows from the north. Mediterranean Overflow Water (MOW) flows from the Strait of Gibraltar as a high salinity water. NADW dominates the deep and overflow layer both in the North and South Atlantic. In the bottom layer, AABW is the only natural water mass with high silicate signature spreading from the Antarctic to the North Atlantic. Due to the change of water mass properties, in this work we renamed to North East Antarctic Bottom Water NEABW north of the equator. Similarly, the distributions of Labrador Sea Water (LSW), Iceland Scotland Overflow Water (ISOW), and Denmark Strait Overflow Water (DSOW) forms upper and lower portion of NADW, respectively roughly south of the Grand Banks between ~50 and 66°N. In the far south the distributions of Circumpolar Deep Water (CDW) and Weddell Sea Bottom Water (WSBW) are of significance to understand the formation of the AABW.

Description: The northeast Atlantic encompasses archetypal examples of volcanic rifted margins. Twenty-five years after the last ODP (Ocean Drilling Program) leg on these volcanic margins, the reasons for excess melting are still disputed with at least three competing hypotheses being discussed. We are proposing a new drilling campaign that will constrain the timing, rates of volcanism, and vertical movements of rifted margins. This will allow us to parameterise geodynamic models that can distinguish between the hypotheses. Furthermore, the drilling-derived data will help us to understand the role of breakup magmatism as a potential driver for the Palaeocene–Eocene thermal maximum (PETM) and its influence on the oceanographic circulation in the earliest phase of the northeast Atlantic Ocean formation. Tackling these questions with a new drilling campaign in the northeast Atlantic region will advance our understanding of the long-term interactions between tectonics, volcanism, oceanography, and climate and the functioning of subpolar northern ecosystems and climate during intervals of extreme warmth.

Description: Recent surveys of marine microbial diversity have identified a previously unrecognized lineage of diplonemid protists as being among the most diverse heterotrophic eukaryotes in global oceans. Despite their monophyly (and assumed importance), they lack a formal taxonomic description, and are informally known as deep-sea pelagic diplonemids (DSPDs) or marine diplonemids. Recently, we documented morphology and molecular sequences from several DSPDs, one of which is particularly widespread and abundant in environmental sequence data. To simplify the communication of future work on this important group, here we formally propose to erect the family Eupelagonemidae to encompass this clade, as well as a formal genus and species description for the apparently most abundant phylotype, Eupelagonema oceanica, for which morphological information and single-cell amplified genome data are currently available. © 2018 International Society of Protistologists

Description: Spores of the dinoflagellate Chytriodinium are known to infest copepod eggs causing their lethality. Despite the potential to control the population of such an ecologically important host, knowledge about Chytriodinium parasites is limited: we know little about phylogeny, parasitism, abundance, or geographical distribution. We carried out genome sequence surveys on four manually isolated sporocytes from the same sporangium, which seemed to be attached to a copepod nauplius, to analyze the phylogenetic position of Chytriodinium based on SSU and concatenated SSU/LSU rRNA gene sequences, and also characterize two genes related to the plastidial heme pathway, hemL and hemY. The results suggest the presence of a cryptic plastid in Chytriodinium and a photosynthetic ancestral state of the parasitic Chytriodinium/Dissodinium clade. Finally, by mapping Tara Oceans V9 SSU amplicon data to the recovered SSU rRNA gene sequences from the sporocytes, we show that globally, Chytriodinium parasites are most abundant within the pico/nano- and mesoplankton of the surface ocean and almost absent within microplankton, a distribution indicating that they generally exist either as free-living spores or host-associated sporangia. © 2018 International Society of Protistologists

Description: The El Niño Southern Oscillation (ENSO) with its warm (El Niño) and cold (La Niña) phase has strong impacts on marine ecosystems off Peru. This influence extends from changes in nutrient availability to productivity and oxygen levels. While several studies have demonstrated the influence of ENSO events on biological productivity, less is known about their impact on oxygen concentrations. In situ observations along the Peruvian and Chilean coast have shown a strong water column oxygenation during the 1997/1998 strong El Niño event. These observations suggest a deepening of the oxygen minimum zone (OMZ) along the continental shelf. However, due to reduced spatial coverage of the existing in situ observations, no studies have yet demonstrated the OMZ response to El Niño events in the whole Eastern Tropical South Pacific (ETSP). Furthermore, most studies have focused on El Niño events. Much less attention was given to the oxygen dynamics under La Niña influence. Here, we provide a comprehensive analysis of the ENSO influence on OMZ dynamics. Interannual variability of the OMZ during the period 1990–2010 is derived from a regional coupled physical-biogeochemical model forced with realistic atmospheric and lateral boundary conditions. Our results show a reduction of the vertical extent and a deepening of suboxic waters (SW) during the El Niño phase. During the La Niña phase, there is a vertical expansion of SW. These fluctuations in OMZ extent are due to changes in oxygen supply into its core depth mainly from lateral margins. During the El Niño phase, the enhanced lateral oxygen supply from the subtropics is the main reason for the reduction of SW in both coastal and offshore regions. During the La Niña phase, the oxygenated subtropical waters are blocked by the poleward transport along the southern margin of the OMZ. Consequently, oxygen concentrations within the OMZ are reduced and suboxic conditions expand during La Niña. The detailed analysis of transport pathways presented here provides new insights into how ENSO variability affects the oxygen-sensitive marine biogeochemistry of the ETSP.

Description: Oxygen minimum zones (OMZs) in the open ocean occur below the surface in regions of weak ventilation and high biological productivity. Very low levels of dissolved oxygen affect marine life and alter biogeochemical cycles. One of the most intense but least understood OMZs in the world is located in the Arabian Sea in a depth range between 300 to 1000 m. Within the last decades observations suggest a decreasing oxygen trend. Thus, an improved understanding of the crucial processes is necessary for a reliable assessment of the future development of the Arabian Sea OMZ. This study uses a combination of observational data as well as reanalysis velocity fields from the ocean model Hycom (Hybrid Coordinate Ocean Model) to explore the ventilation dynamics of the Arabian Sea OMZ. Our results show that the OMZ features a strong seasonal cycle with regional differences that is correlated with the monsoon system: In the eastern basin, the OMZ is strongest during the winter monsoon with a core thickness of 1000 m depth and oxygen values of less than 5 µmol/kg. Ventilation during that phase is dominated by Persian Gulf water, that clockwise circles the perimeter of the basin and enters the OMZ from the north. During the summer monsoon ventilation from the southeast leads to higher oxygen values indicating a reverse flow along the Indian coast in the intermediate layer compared to the southeastward surface currents. The seasonal cycle in the western basin has the same seasonality as the one in the eastern basin with a core thickness of 900 m during the winter monsoon. The oxygen supply during the summer monsoon is weaker compared to the eastern basin and correlates with the ventilation of Persian Gulf (Red Sea) water during the summer monsoon (autumn inter-monsoon) phase. As the interior exchange between the eastern and western basin is weak, the more pronounced OMZ in the eastern basin is explained by prolonged ventilation time scales. For the eastern (western) basin Persian Gulf water needs 2–3 (1–2) years and Red Sea water 7–8 (3–4) years to ventilate the OMZ.

Description: Ice sheet numerical modeling is an important tool to estimate the dynamic contribution of the Antarctic ice sheet to sea level rise over the coming centuries. The influence of initial conditions on ice sheet model simulations, however, is still unclear. To better understand this influence, an initial state intercomparison exercise (initMIP) has been developed to compare, evaluate, and improve initialization procedures and estimate their impact on century-scale simulations. initMIP is the first set of experiments of the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6), which is the primary Coupled Model Intercomparison Project Phase 6 (CMIP6) activity focusing on the Greenland and Antarctic ice sheets. Following initMIP-Greenland, initMIP-Antarctica has been designed to explore uncertainties associated with model initialization and spin-up and to evaluate the impact of changes in external forcings. Starting from the state of the Antarctic ice sheet at the end of the initialization procedure, three forward experiments are each run for 100 years: a control run, a run with a surface mass balance anomaly, and a run with a basal melting anomaly beneath floating ice. This study presents the results of initMIP-Antarctica from 25 simulations performed by 16 international modeling groups. The submitted results use different initial conditions and initialization methods, as well as ice flow model parameters and reference external forcings. We find a good agreement among model responses to the surface mass balance anomaly but large variations in responses to the basal melting anomaly. These variations can be attributed to differences in the extent of ice shelves and their upstream tributaries, the numerical treatment of grounding line, and the initial ocean conditions applied, suggesting that ongoing efforts to better represent ice shelves in continental-scale models should continue.

Description: We present here the first results, for the preindustrial and mid-Holocene climatological periods, of the newly developed isotope-enhanced version of the fully coupled Earth system model MPI-ESM, called hereafter MPI-ESM-wiso. The water stable isotopes H162O, H182O and HDO have been implemented into all components of the coupled model setup. The mid-Holocene provides the opportunity to evaluate the model response to changes in the seasonal and latitudinal distribution of insolation induced by different orbital forcing conditions. The results of our equilibrium simulations allow us to evaluate the performance of the isotopic model in simulating the spatial and temporal variations of water isotopes in the different compartments of the hydrological system for warm climates. For the preindustrial climate, MPI-ESM-wiso reproduces very well the observed spatial distribution of the isotopic content in precipitation linked to the spatial variations in temperature and precipitation rate. We also find a good model–data agreement with the observed distribution of isotopic composition in surface seawater but a bias with the presence of surface seawater that is too 18O-depleted in the Arctic Ocean. All these results are improved compared to the previous model version ECHAM5/MPIOM. The spatial relationships of water isotopic composition with temperature, precipitation rate and salinity are consistent with observational data. For the preindustrial climate, the interannual relationships of water isotopes with temperature and salinity are globally lower than the spatial ones, consistent with previous studies. Simulated results under mid-Holocene conditions are in fair agreement with the isotopic measurements from ice cores and continental speleothems. MPI-ESM-wiso simulates a decrease in the isotopic composition of precipitation from North Africa to the Tibetan Plateau via India due to the enhanced monsoons during the mid-Holocene. Over Greenland, our simulation indicates a higher isotopic composition of precipitation linked to higher summer temperature and a reduction in sea ice, shown by positive isotope–temperature gradient. For the Antarctic continent, the model simulates lower isotopic values over the East Antarctic plateau, linked to the lower temperatures during the mid-Holocene period, while similar or higher isotopic values are modeled over the rest of the continent. While variations of isotopic contents in precipitation over West Antarctica between mid-Holocene and preindustrial periods are partly controlled by changes in temperature, the transport of relatively 18O-rich water vapor near the coast to the western ice core sites could play a role in the final isotopic composition. So, more caution has to be taken about the reconstruction of past temperature variations during warm periods over this area. The coupling of such a model with an ice sheet model or the use of a zoomed grid centered on this region could help to better describe the role of the water vapor transport and sea ice around West Antarctica. The reconstruction of past salinity through isotopic content in sea surface waters can be complicated for regions with strong ocean dynamics, variations in sea ice regimes or significant changes in freshwater budget, giving an extremely variable relationship between the isotopic content and salinity of ocean surface waters over small spatial scales. These complicating factors demonstrate the complexity of interpreting water isotopes as past climate signals of warm periods like the mid-Holocene. A systematic isotope model intercomparison study for further insights on the model dependency of these results would be beneficial.

Description: Probabilistic spatial reconstructions of past climate states are valuable to quantitatively study the climate system under different forcing conditions because they combine the information contained in a proxy synthesis into a comprehensible product. Unfortunately, they are subject to a complex uncertainty structure due to complicated proxy–climate relations and sparse data, which makes interpolation between samples difficult. Bayesian hierarchical models feature promising properties to handle these issues, like the possibility to include multiple sources of information and to quantify uncertainties in a statistically rigorous way. We present a Bayesian framework that combines a network of pollen and macrofossil samples with a spatial prior distribution estimated from a multi-model ensemble of climate simulations. The use of climate simulation output aims at a physically reasonable spatial interpolation of proxy data on a regional scale. To transfer the pollen data into (local) climate information, we invert a forward version of the probabilistic indicator taxa model. The Bayesian inference is performed using Markov chain Monte Carlo methods following a Metropolis-within-Gibbs strategy. Different ways to incorporate the climate simulations into the Bayesian framework are compared using identical twin and cross-validation experiments. Then, we reconstruct the mean temperature of the warmest and mean temperature of the coldest month during the mid-Holocene in Europe using a published pollen and macrofossil synthesis in combination with the Paleoclimate Modelling Intercomparison Project Phase III mid-Holocene ensemble. The output of our Bayesian model is a spatially distributed probability distribution that facilitates quantitative analyses that account for uncertainties.

Description: Heteroplasmy occurs when copies of an organellar genome (plastid or mitochondrial) differ from one another either within a cell or among cells within an individual. This phenomenon was first discovered in plastids over 100 years ago, though ‘heteroplasmy’ was not formally defined until decades later. Mitochondrial and plastid heteroplasmy have since been discovered in diverse taxa, including numerous plants, particularly those with the gynodioecious breeding system. Though heteroplasmy can arise through mutations in organellar genomes, biparental inheritance of organelles often generates heteroplasmy. This article is intended to summarise the cytoploid nature of organellar genomes and different aspects of mitochondrial and plastid heteroplasmy, from its discovery, to the development of detection methods. We will also discuss the effects of heteroplasmy on individuals, populations, and species, with regard to fitness, cytonuclear interactions, and genome evolution. Finally, we highlight a few research questions which we consider in need of further investigation to understand the prevalence and effects of heteroplasmy in natural settings.

Description: Large calderas are among the Earth's major volcanic features. They are associated with large magma reservoirs and elevated geothermal gradients. Caldera-forming eruptions result from the withdrawal and collapse of the magma chambers and produce large-volume pyroclastic deposits and later-stage deformation related to post-caldera resurgence and volcanism. Unrest episodes are not always followed by an eruption; however, every eruption is preceded by unrest. The Campi Flegrei caldera (CFc), located along the eastern Tyrrhenian coastline in southern Italy, is close to the densely populated area of Naples. It is one of the most dangerous volcanoes on Earth and represents a key example of an active, resurgent caldera. It has been traditionally interpreted as a nested caldera formed by collapses during the 100–200 km3 Campanian Ignimbrite (CI) eruption at ∼39 ka and the 40 km3 eruption of the Neapolitan Yellow Tuff (NYT) at ∼15 ka. Recent studies have suggested that the CI may instead have been fed by a fissure eruption from the Campanian Plain, north of Campi Flegrei. A MagellanPlus workshop was held in Naples, Italy, on 25–28 February 2017 to explore the potential of the CFc as target for an amphibious drilling project within the International Ocean Discovery Program (IODP) and the International Continental Drilling Program (ICDP). It was agreed that Campi Flegrei is an ideal site to investigate the mechanisms of caldera formation and associated post-caldera dynamics and to analyze the still poorly understood interplay between hydrothermal and magmatic processes. A coordinated onshore–offshore drilling strategy has been developed to reconstruct the structure and evolution of Campi Flegrei and to investigate volcanic precursors by examining (a) the succession of volcanic and hydrothermal products and related processes, (b) the inner structure of the caldera resurgence, (c) the physical, chemical, and biological characteristics of the hydrothermal system and offshore sediments, and (d) the geological expression of the phreatic and hydromagmatic eruptions, hydrothermal degassing, sedimentary structures, and other records of these phenomena. The deployment of a multiparametric in situ monitoring system at depth will enable near-real-time tracking of changes in the magma reservoir and hydrothermal system.

Description: The international endeavour to retrieve a continuous ice core, which spans the middle Pleistocene climate transition ca. 1.2–0.9 Myr ago, encompasses a multitude of field and model-based pre-site surveys. We expand on the current efforts to locate a suitable drilling site for the oldest Antarctic ice core by means of 3-D continental ice-sheet modelling. To this end, we present an ensemble of ice-sheet simulations spanning the last 2 Myr, employing transient boundary conditions derived from climate modelling and climate proxy records. We discuss the imprint of changing climate conditions, sea level and geothermal heat flux on the ice thickness, and basal conditions around previously identified sites with continuous records of old ice. Our modelling results show a range of configurational ice-sheet changes across the middle Pleistocene transition, suggesting a potential shift of the West Antarctic Ice Sheet to a marine-based configuration. Despite the middle Pleistocene climate reorganisation and associated ice-dynamic changes, we identify several regions conducive to conditions maintaining 1.5 Myr (million years) old ice, particularly around Dome Fuji, Dome C and Ridge B, which is in agreement with previous studies. This finding strengthens the notion that continuous records with such old ice do exist in previously identified regions, while we are also providing a dynamic continental ice-sheet context.

Description: While nonlinear stochastic partial differential equations arise naturally in spatiotemporal modeling, inference for such systems often faces two major challenges: sparse noisy data and ill-posedness of the inverse problem of parameter estimation. To overcome the challenges, we introduce a strongly regularized posterior by normalizing the likelihood and by imposing physical constraints through priors of the parameters and states. We investigate joint parameter-state estimation by the regularized posterior in a physically motivated nonlinear stochastic energy balance model (SEBM) for paleoclimate reconstruction. The high-dimensional posterior is sampled by a particle Gibbs sampler that combines a Markov chain Monte Carlo (MCMC) method with an optimal particle filter exploiting the structure of the SEBM. In tests using either Gaussian or uniform priors based on the physical range of parameters, the regularized posteriors overcome the ill-posedness and lead to samples within physical ranges, quantifying the uncertainty in estimation. Due to the ill-posedness and the regularization, the posterior of parameters presents a relatively large uncertainty, and consequently, the maximum of the posterior, which is the minimizer in a variational approach, can have a large variation. In contrast, the posterior of states generally concentrates near the truth, substantially filtering out observation noise and reducing uncertainty in the unconstrained SEBM.

Description: We present a software toolbox that allows the efficient collection, management and maintenance of larger paleoceanographic data sets. The program combines a graphical user interface (GUI) with a simple document-based database and functionality for visualization, stratigraphy (visual alignment and radiocarbon calibration), age modelling and efficient ensemble time-series generation to create larger homogenous data compilations. Data can be imported from Excel or text files, are stored locally in netCDF format and can be easily exchanged between collaborating scientists. Within a data collection, data can be imported either to proxy-specific sub-collections or to a multi-proxy (“miscellaneous”) sub-collection that allows the import and management of any downcore data. A single age model is shared among all proxies of a core within a collection. The stand-alone software can be used with Windows and macOS and does not require web access. Installers of the current version for both Windows 10 and macOS including the C++ code can be downloaded from https://www.marum.de/Stefan-Mulitza/PaleoDataView.html (last access: 5 December 2019) along with a detailed user guide.

Description: Combining ocean general circulation models with proxy data via data assimilation is a means to obtain estimates of past ocean states that are consistent with model physics as well as with proxy data. The climate during the Last Glacial Maximum (LGM, 19–23 ka) was substantially different from today. Even though boundary conditions are comparatively well known, the large-scale patterns of the ocean circulation during this time remain uncertain. Previous efforts to combine ocean models with proxy data have shown dissimilar results regarding the state of the ocean, in particular of the Atlantic Meridional Overturning Circulation. Here, we present a new LGM ocean state estimate that extents previous estimates by using global benthic as well as planktic data on the oxygen isotopic composition of calcite. It is further constrained by global seasonal and annual sea surface temperature (SST) reconstructions. The estimate shows an Atlantic Ocean that is similar to the Late Holocene Atlantic Ocean but with a reduced formation of Antarctic Bottom Water, in contrast to results of previous studies. The results indicate that SST and oxygen isotopic data alone do not require the presence of a shallower North Atlantic Deep Water and a more extensive Antarctic Bottom Water, and highlight the need for more proxy data of different types to obtain reliable ocean state estimates. Additional adjoint sensitivity experiments reveal that data from the deep North Atlantic and from the global deep Southern Ocean are most important to constrain the Atlantic Meridional Overturning Circulation.

Description: The species composition of many groups of marine plankton appears well predicted by sea surface temperature (SST). Consequently, fossil plankton assemblages have been widely used to reconstruct past SST. Most applications of this approach make use of the highest possible taxonomic resolution. However, not all species are sensitive to temperature, and their distribution may be governed by other parameters. There are thus reasons to question the merit of including information about all species, both for transfer function performance and for its effect on reconstructions. Here we investigate the effect of species selection on planktonic foraminifera transfer functions. We assess species importance for transfer function models using a random forest technique and evaluate the performance of models with an increasing number of species. Irrespective of using models that use the entire training set (weighted averaging) or models that use only a subset of the training set (modern analogue technique), we find that the majority of foraminifera species does not carry useful information for temperature reconstruction. Less than one-third of the species in the training set is required to provide a temperature estimate with a prediction error comparable to a transfer function that uses all species in the training set. However, species selection matters for paleotemperature estimates. We find that transfer function models with a different number of species but with the same error may yield different reconstructions of sea surface temperature when applied to the same fossil assemblages. This ambiguity in the reconstructions implies that fossil assemblage change reflects a combination of temperature and other environmental factors. The contribution of the additional factors is site and time specific, indicating ecological and geological complexity in the formation of the sedimentary assemblages. The possibility of obtaining multiple different reconstructions from a single sediment record presents a previously unrecognized source of uncertainty for sea surface temperature estimates based on planktonic foraminifera assemblages. This uncertainty can be evaluated by determining the sensitivity of the reconstructions to species pruning.

Description: Although narrow‐azimuth towed‐streamer data provides good image quality for structural interpretation, it is generally accepted that for wide‐azimuth marine surveys seabed receivers deliver superior seismic reflection measurements and seismically derived reservoir attributes. However, seabed surveys are not widely used due to the higher acquisition costs when compared to streamer acquisition. In recent years, there have been significant engineering efforts to automate receiver deployment and retrieval in order to minimize the cost differential and conduct cost‐efficient seabed receiver seismic surveys. These engineering efforts include industrially engineered nodes, nodes‐on‐a‐rope deployment schemes and even robotic nodes, which swim to and from the deployment location. This move to automation is inevitable, leading to robotization of seismic data acquisition for exploration and development activities in the oil and gas industry. We are developing a robotic‐based technology, which utilizes autonomous underwater vehicles as seismic sensors without the need of using a remotely operated vehicle for deployment and retrieval. In this paper, we describe the autonomous underwater vehicle evolution throughout the project years from initial heavy and bulky nodes to fully autonomous light and flexible underwater receivers. Results obtained from two field pilot tests using different generations of autonomous underwater vehicles indicate that the seismic coupling, and navigation based on underwater acoustics are very reliable and robust.

Description: Regional climate modeling bridges the gap between the coarse resolution of current global climate models and the regional-to-local scales, where the impacts of climate change are of primary interest. Here, we present a review of the added value of the regional climate modeling approach within the scope of paleoclimate research and discuss the current major challenges and perspectives. Two time periods serve as an example: the Holocene, including the Last Millennium, and the Last Glacial Maximum. Reviewing the existing literature reveals the benefits of regional paleo climate modeling, particularly over areas with complex terrain. However, this depends largely on the variable of interest, as the added value of regional modeling arises from a more realistic representation of physical processes and climate feedbacks compared to global climate models, and this affects different climate variables in various ways. In particular, hydrological processes have been shown to be better represented in regional models, and they can deliver more realistic meteorological data to drive ice sheet and glacier modeling. Thus, regional climate models provide a clear benefit to answer fundamental paleoclimate research questions and may be key to advance a meaningful joint interpretation of climate model and proxy data.

Description: A marine seismic method based on continuous source and receiver wavefields has been developed. The method requires continuous recording of the seismic data. The source that may consist of multiple source elements can emit signals continuously while moving. The ideal source wavefield to be used with this method should be as white as possible both in a temporal and a spatial sense to avoid deep notches in the spectrum enabling a stable multi‐dimensional deconvolution. White noise has such properties. However, equipment that can generate white noise does not exist. In order to generate a continuous source wavefield that is approaching the properties of white noise using existing equipment onboard marine seismic vessels, individual air‐guns can be triggered with short randomized time intervals in a near‐continuous fashion. The main potential benefits with the method are to reduce the environmental impact of marine seismic surveys and to improve acquisition efficiency. The peak sound pressure levels are significantly reduced by triggering one air‐gun at a time compared to conventional marine seismic sources. Sound exposure levels are also reduced in most directions. Since the method is based on continuous recording of seismic data and the air‐guns are triggered based on time and not based on position, there are less vessel speed limitations compared to conventional marine seismic data acquisition. Also, because the source wavefield is spread out in time, the wavefields emitted from source elements in different cross‐line positions can be designed such that the emitted wavefield is spatially white in this direction. This means that source elements in multiple cross‐line positions can be operated simultaneously, potentially improving the cross‐line sampling and/or the acquisition efficiency.

Description: Basalt weathering is one of many relevant processes balancing the global carbon cycle via land–ocean alkalinity fluxes. The CO2 consumption by weathering can be calculated using alkalinity and is often scaled with runoff and/or temperature. Here, it is tested if the surface age distribution of a volcanic system derived by geological maps is a useful proxy for changes in alkalinity production with time. A linear relationship between temperature normalized alkalinity fluxes and the Holocene area fraction of a volcanic field was identified using information from 33 basalt volcanic fields, with an r2=0.93. This relationship is interpreted as an aging function and suggests that fluxes from Holocene areas are ∼10 times higher than those from old inactive volcanic fields. However, the cause for the decrease with time is probably a combination of effects, including a decrease in alkalinity production from material in the shallow critical zone as well as a decline in hydrothermal activity and magmatic CO2 contribution. The addition of fresh reactive material on top of the critical zone has an effect in young active volcanic settings which should be accounted for, too. A comparison with global models suggests that global alkalinity fluxes considering Holocene basalt areas are ∼60 % higher than the average from these models imply. The contribution of Holocene areas to the global basalt alkalinity fluxes is today however only ∼5 %, because identified, mapped Holocene basalt areas cover only ∼1 % of the existing basalt areas. The large trap basalt proportion on the global basalt areas today reduces the relevance of the aging effect. However, the aging effect might be a relevant process during periods of globally intensive volcanic activity, which remains to be tested.

Description: Dynamic vegetation models simulate global vegetation in terms of fractional coverage of a few plant functional types (PFTs). Although these models often share the same concept, they differ with respect to the number and kind of PFTs, complicating the comparability of simulated vegetation distributions. Pollen-based vegetation reconstructions are initially only available in the form of time series of individual taxa that are not distinguished in the models. Thus, to evaluate simulated vegetation distributions, the modelling results and pollen-based vegetation reconstructions have to be converted into a comparable format. The classical approach is the method of biomisation, but hitherto PFT-based biomisation methods were only available for individual models. We introduce and evaluate a simple, universally applicable technique to harmonise PFT distributions by assigning them into nine mega-biomes, using only assumptions on the minimum PFT cover fractions and few bioclimatic constraints (based on the 2 m temperature). These constraints mainly follow the limitation rules used in the classical biome models (here BIOME4). We test the method for six state-of-the-art dynamic vegetation models that are included in Earth system models based on pre-industrial, mid-Holocene and Last Glacial Maximum simulations. The method works well, independent of the spatial resolution or the complexity of the models. Large biome belts (such as tropical forest) are generally better represented than regionally confined biomes (warm–temperate forest, savanna). The comparison with biome distributions inferred via the classical biomisation approach of forcing biome models (here BIOME1) with the simulated climate states shows that the PFT-based biomisation is even able to keep up with the classical method. However, as the new method considers the PFT distributions actually calculated by the Earth system models, it allows for a direct comparison and evaluation of simulated vegetation distributions which the classical method cannot do. Thereby, the new method provides a powerful tool for the evaluation of Earth system models in general.

Description: The sedimentary stable nitrogen isotope compositions of bulk organic matter (δ15Nbulk) and silicon isotope composition of diatoms (δ30SiBSi) both mainly reflect the degree of past nutrient utilization by primary producers. However, in ocean areas where anoxic and suboxic conditions prevail, the δ15Nbulk signal ultimately recorded within the sediments is also influenced by water column denitrification causing an increase in the subsurface δ15N signature of dissolved nitrate (δ15NO3−) upwelled to the surface. Such conditions are found in the oxygen minimum zone off Peru, where at present an increase in subsurface δ15NO3− from North to South along the shelf is observed due to ongoing denitrification within the pole-ward flowing subsurface waters, while the δ30Si signature of silicic acid (δ30Si(OH)4) at the same time remains unchanged. Here, we present three new δ30SiBSi records between 11° S and 15° S and compare these to previously published δ30SiBSi and δ15Nbulk records from Peru covering the past 600 years. We present a new approach to calculate past subsurface δ15NO3− signatures based on the correlation of δ30SiBSi and δ15Nbulk signatures at a latitudinal resolution for different time periods. Our results show source water δ15NO3− compositions during the last 200 years, the Current Warm Period (CWP) and during short-term arid events prior to that, which are close to modern values increasing southward from 7 to 10 ‰ (between 11° S and 15° S). In contrast, humid conditions during the Little Ice Age (LIA) reflect consistently low δ15NO3− values between 6 and 7.5‰. Furthermore, we are able to relate the short-term variability in both isotope compositions to changes in the ratio of nutrients (NO3− : Si(OH)4) taken up by different dominating phytoplankton groups (diatoms and non-siliceous phytoplankton) under the variable climatic conditions of the past 600 years.

Description: The open ocean is a major source of nitrous oxide (N2O), an atmospheric trace gas attributable to global warming and ozone depletion. Intense sea-to-air N2O fluxes occur in major oceanic upwelling regions such as the eastern tropical South Pacific (ETSP). The ETSP is influenced by the El Niño–Southern Oscillation that leads to inter-annual variations in physical, chemical, and biological properties in the water column. In October 2015, a strong El Niño event was developing in the ETSP; we conduct field observations to investigate (1) the N2O production pathways and associated biogeochemical properties and (2) the effects of El Niño on water column N2O distributions and fluxes using data from previous non-El Niño years. Analysis of N2O natural abundance isotopomers suggested that nitrification and partial denitrification (nitrate and nitrite reduction to N2O) were occurring in the near-surface waters; indicating that both pathways contributed to N2O effluxes. Higher-than-normal sea surface temperatures were associated with a deepening of the oxycline and the oxygen minimum layer. Within the shelf region, surface N2O supersaturation was nearly an order of magnitude lower than that of non-El Niño years. Therefore, a significant reduction of N2O efflux (75 %–95 %) in the ETSP occurred during the 2015 El Niño. At both offshore and coastal stations, the N2O concentration profiles during El Niño showed moderate N2O concentration gradients, and the peak N2O concentrations occurred at deeper depths during El Niño years; this was likely the result of suppressed upwelling retaining N2O in subsurface waters. At multiple stations, water-column inventories of N2O within the top 1000 m were up to 160 % higher than those measured in non-El Niño years, indicating that subsurface N2O during El Niño could be a reservoir for intense N2O effluxes when normal upwelling is resumed after El Niño.

Description: We investigate the role of the tropics, the stratosphere, and atmosphere‐ocean coupling for seasonal forecasts of strong, potentially damaging, Northern Hemisphere extratropical winter wind storm frequencies. This is done by means of relaxation experiments with the European Centre for Medium‐Range Weather Forecasts model, which allow us to prescribe perfect forecasts for specific parts of the coupled atmosphere‐ocean system. We find that perfect predictions of the Northern Hemisphere stratosphere significantly enhance winter storm predictive skill between eastern Greenland and Northern Europe. Correct seasonal predictions of the occurrence of stratospheric sudden warmings play a decisive role. The importance of correctly predicting the tropics and of two‐way atmosphere‐ocean coupling, both for forecasting stratospheric sudden warming risk and, correspondingly, severe winter storm frequency, is noted.

Description: The early Eocene (56 to 48 million years ago) is inferred to have been the most recent time that Earth's atmospheric CO2 concentrations exceeded 1000 ppm. Global mean temperatures were also substantially warmer than present day. As such, study of early Eocene climate provides insight into how a super-warm Earth system behaves and offers an opportunity to 10 evaluate climate models under conditions of high greenhouse gas forcing. The Deep Time Model Intercomparison Project (DeepMIP) is a systematic model-model and model-data intercomparison of three early Paleogene time slices: latest Paleocene, Paleocene-Eocene thermal maximum and early Eocene climatic optimum. A previous article outlined the model experimental design for climate model simulations. In this article, we outline the methodologies to be used for the compilation and analysis of climate proxy data, primarily proxies for temperature and CO2. This paper establishes the protocols for a concerted and 15 coordinated effort to compile the climate proxy records across a wide geographic range. The resulting climate "atlas" will be used to constrain and evaluate climate models for the three selected time intervals, and provide insights into the mechanisms that control these warm climate states. We provide version 0.1 of this database, in anticipation that this will be expanded in subsequent publications.

Description: The northern part of the South China Sea is characterized by widespread occurrence of bottom simulating reflectors (BSR) indicating the presence of marine gas hydrate. Because the area covers both a tectonically inactive passive margin and the termination of a subduction zone, the influence of tectonism on the dynamics of gas hydrate systems can be studied in this region. Geophysical data show that there are multiple thrust faults on the active margin while much fewer and smaller faults exist in the passive margin. This tectonic difference matches with a difference in the geophysical characteristics of the gas hydrate systems. High hydrate saturation derived from ocean bottom seismometer data and controlled source electromagnetic data and conspicuous high‐amplitude reflections in P‐Cable 3D seismic data above the BSR are found in the anticlinal ridges of the active margin. In contrast all geophysical evidence for the passive margin points to normal to low hydrate saturations. Geochemical analyses of gas samples collected at seep sites on the active margin show methane with heavy δ13C isotope composition, while gas collected at the passive margin shows light carbon isotope composition. Thus, we interpret the passive margin as a typical gas hydrate province fuelled by biogenic production of methane and the active margin gas hydrate system as a system that is fuelled not only by biogenic gas production but also by additional advection of thermogenic methane from the subduction system.

Description: The mechanisms controlling the variability of oxygen levels in the ocean are poorly quantified. We focus here on the impact of wind synoptic variability associated with tropical convective regions and extra‐tropical storms. Removing the wind higher frequencies of variability (2 days – 1 month) in an atmosphere reanalysis used to force an ocean model decreases wind stress by up to 20% in the tropics and 50% in the mid‐latitudes, weakening wind‐driven ocean circulation by 20%. Oxygen levels decrease by up to 10 mmol.m‐3 in tropical oceans and 30 mmol.m‐3 in subtropical gyres mainly due to changes in advective processes. While a large part of the tropical oxygen anomaly has local origins, changes in oxygen levels in the subtropical gyres modulate tropical oxygen distribution. Our study suggests that the “storminess” of the ocean is an important parameter that could determine the future evolution of poorly oxygenated regions.

Description: The Benguela Upwelling system (BUS) is the most productive of all eastern boundary upwelling ecosystems and it hosts a well‐developed oxygen minimum zone. As such, the BUS is a potential hotspot for production of N2O, a potent greenhouse gas derived from microbially‐driven decay of sinking organic matter. Yet, the extent at which near‐surface waters emit N2O to the atmosphere in the BUS is highly uncertain. Here we present the first high‐resolution surface measurements of N2O across the northern part of the BUS (nBUS). We found strong gradients with a three‐fold increase in N2O concentrations near the coast as compared with open ocean waters. Our observations show enhanced sea‐to‐air fluxes of N2O (up to 1.67 nmol m−2 s−1) in association with local upwelling cells. Based on our data we suggest that the nBUS can account for 13% of the total coastal upwelling source of N2O to the atmosphere.

Description: Hydroxylamine (NH 2 OH), a short-lived intermediate in the nitrogen cycle, is a potential precursor of nitrous oxide (N 2 O) in the ocean. However, measurements of NH 2 OH in the ocean are sparse. Here we present a data set of depth profiles of NH 2 OH from the equatorial Atlantic Ocean and the eastern tropical South Pacific and compare it to N 2 O, nitrate, and nitrite profiles under varying oxygen conditions. The presence of NH 2 OH in surface waters points toward surface nitrification in the upper 100 m. Overall, we found a ratio of 1:3 between NH 2 OH and N 2 O in open ocean areas when oxygen concentrations were 〉50 μmol/L. In the equatorial Atlantic Ocean and the open ocean eastern tropical South Pacific, where nitrification is the dominant N 2 O production pathway, stepwise multiple regressions demonstrated that N 2 O, NH 2 OH, and nitrate concentrations were highly correlated, suggesting that NH 2 OH is a potential indicator for nitrification.

Description: Radiative forcing from volcanic aerosol impacts surface temperatures; however, the background climate state also affects the response. A key question thus concerns whether constraining forcing estimates is more important than constraining initial conditions for accurate simulation and attribution of posteruption climate anomalies. Here we test whether different realistic volcanic forcing magnitudes for the 1815 Tambora eruption yield distinguishable ensemble surface temperature responses. We perform a cluster analysis on a superensemble of climate simulations including three 30-member ensembles using the same set of initial conditions but different volcanic forcings based on uncertainty estimates. Results clarify how forcing uncertainties can overwhelm initial-condition spread in boreal summer due to strong direct radiative impact, while the effect of initial conditions predominate in winter, when dynamics contribute to large ensemble spread. In our setup, current uncertainties affecting reconstruction-simulation comparisons prevent conclusions about the magnitude of the Tambora eruption and its relation to the “year without summer.”

Description: Sea surface salinity is one of the most important parameters to reconstruct in paleoclimatology, reflecting amongst others the hydrological cycle, paleo-density, ice volume, and regional and global circulation of water masses. Recent culture studies and a Red Sea field study revealed a significant positive relation between salinity and Na incorporation within benthic and planktonic foraminiferal shells. However, these studies reported varying partitioning of Na between and within the same species. The latter could be associated with ontogenetic variations, most likely spine loss. Varying Na concentrations were observed in different parts of foraminiferal shells, with especially spines and regions close to the primary organic sheet being enriched in Na. In this study, we unravel the Na composition of different components of the planktonic foraminiferal shell wall using Electron Probe Micro Analysis (EPMA) and solution-ICP-MS. A model is presented to interpret EPMA data for spines and spine bases to quantitatively assess differences in composition and contribution to whole shell Na/Ca signals. The same model can also be applied to other spatial inhomogeneities observed in foraminiferal shell chemistry, like elemental (e.g. Mg, Na, S) banding and/or hotspots. The relative contribution of shell calcite, organic linings, spines and spine bases to whole shell Na chemistry is considered quantitatively. This study shows that whereas the high Na areas may be susceptible to taphonomy, the Na chemistry of the shell itself seems relatively robust. Comparing both shell and spine Na/Ca values with salinity shows that shell chemistry records salinity, albeit with a very modest slope.

Description: Increasing global energy demands have led to the ongoing intensification of hydrocarbon extraction from marine areas. Hydrocarbon extractive activities pose threats to native marine biodiversity, such as noise, light, and chemical pollution, physical changes to the sea floor, invasive species, and greenhouse gas emissions. Here, we assessed at a global scale the spatial overlap between offshore hydrocarbon activities and marine biodiversity (〉25,000 species, nine major ecosystems, and marine protected areas), and quantify the changes over time. We discovered that two-thirds of global offshore hydrocarbon activities occur in areas within the top 10% for species richness, range rarity, and proportional range rarity values globally. Thus, while hydrocarbon activities are undertaken in less than one percent of the ocean's area, they overlap with approximately 85% of all assessed species. Of conservation concern, 4% of species with the largest proportion of their range overlapping hydrocarbon activities are range restricted, potentially increasing their vulnerability to localized threats such as oil spills. While hydrocarbon activities have extended to greater depths since the mid-1990s, we found that the largest overlap is with coastal ecosystems, particularly estuaries, saltmarshes and mangroves. Furthermore, in most countries where offshore hydrocarbon exploration licensing blocks have been delineated, they do not overlap with marine protected areas (MPAs). Although this is positive in principle, many countries have far more licensing block areas than protected areas, and in some instances, MPA coverage is minimal. These findings suggest the need for marine spatial prioritization to help limit future spatial overlap between marine conservation priorities and hydrocarbon activities. Such prioritization can be informed by the spatial and quantitative baseline information provided here. In increasingly shared seascapes, prioritizing management actions that set both conservation and development targets could help minimize further declines of biodiversity and environmental changes at a global scale.

Description: Genetic data have great potential for improving fisheries management by identifying the fundamental management units—that is, the biological populations—and their mixing. However, so far, the number of practical cases of marine fisheries management using genetics has been limited. Here, we used Atlantic cod in the Baltic Sea to demonstrate the applicability of genetics to a complex management scenario involving mixing of two genetically divergent populations. Specifically, we addressed several assumptions used in the current assessment of the two populations. Through analysis of 483 single nucleotide polymorphisms (SNPs) distributed across the Atlantic cod genome, we confirmed that a model of mechanical mixing, rather than hybridization and introgression, best explained the pattern of genetic differentiation. Thus, the fishery is best monitored as a mixed-stock fishery. Next, we developed a targeted panel of 39 SNPs with high statistical power for identifying population of origin and analyzed more than 2,000 tissue samples collected between 2011 and 2015 as well as 260 otoliths collected in 2003/2004. These data provided high spatial resolution and allowed us to investigate geographical trends in mixing, to compare patterns for different life stages and to investigate temporal trends in mixing. We found similar geographical trends for the two time points represented by tissue and otolith samples and that a recently implemented geographical management separation of the two populations provided a relatively close match to their distributions. In contrast to the current assumption, we found that patterns of mixing differed between juveniles and adults, a signal likely linked to the different reproductive dynamics of the two populations. Collectively, our data confirm that genetics is an operational tool for complex fisheries management applications. We recommend focussing on developing population assessment models and fisheries management frameworks to capitalize fully on the additional information offered by genetically assisted fisheries monitoring.

Description: Prediction and Research Moored Array in the Tropical Atlantic (PIRATA) is a multinational program initiated in 1997 in the tropical Atlantic to improve our understanding and ability to predict ocean-atmosphere variability. PIRATA consists of a network of moored buoys providing meteorological and oceanographic data transmitted in real time to address fundamental scientific questions as well as societal needs. The network is maintained through dedicated yearly cruises, which allow for extensive complementary shipboard measurements and provide platforms for deployment of other components of the Tropical Atlantic Observing System. This paper describes network enhancements, scientific accomplishments and successes obtained from the last 10 years of observations, and additional results enabled by cooperation with other national and international programs. Capacity building activities and the role of PIRATA in a future Tropical Atlantic Observing System that is presently being optimized are also described.

Description: Parasites are one of the strongest selective agents in nature. They select for hosts that evolve counter‐adaptive strategies to cope with infection. Helminth parasites are special because they can modulate their hosts’ immune responses. This phenomenon is important in epidemiological contexts, where coinfections may be affected. How different types of hosts and helminths interact with each other is insufficiently investigated. We used the three‐spined stickleback (Gasterosteus aculeatus) – Schistocephalus solidus model to study mechanisms and temporal components of helminth immune modulation. Sticklebacks from two contrasting populations with either high resistance (HR) or low resistance (LR) against S. solidus, were individually exposed to S. solidus strains with characteristically high growth (HG) or low growth (LG) in G. aculeatus. We determined the susceptibility to another parasite, the eye fluke Diplostomum pseudospathaceum, and the expression of 23 key immune genes at three time points after S. solidus infection. D. pseudospathaceum infection rates and the gene expression responses depended on host and S. solidus type and changed over time. Whereas the effect of S. solidus type was not significant after three weeks, T regulatory responses and complement components were upregulated at later time points if hosts were infected with HG S. solidus. HR hosts showed a well orchestrated immune response, which was absent in LR hosts. Our results emphasize the role of regulatory T cells and the timing of specific immune responses during helminth infections. This study elucidates the importance to consider different coevolutionary trajectories and ecologies when studying host‐parasite interactions.

Description: Seafloor massive sulphide (SMS) deposits are of increasing economic interest in order to satisfy the relentless growth in worldwide metal demand. The Trans‐Atlantic Geotraverse (TAG) hydrothermal field at 26°N on the Mid‐Atlantic Ridge hosts several such deposits. This study presents new controlled source electromagnetic (CSEM), bathymetric and magnetic results from the TAG field. Potential SMS targets were selected based on their surface expressions in high‐resolution bathymetric data. High‐resolution reduced‐to‐the‐pole magnetic data show negative anomalies beneath and surrounding the SMS deposits, revealing large areas of hydrothermal alteration. CSEM data, sensitive to the electrical conductivity of SMS mineralization, further reveal a maximum thickness of up to 80 m and conductivities of up to 5 S/m. SMS samples have conductivities of up to a few thousand S/m, suggesting that remotely inferred conductivities represent an average of metal sulphide ores combined with silicified and altered host basalt that likely dominates at greater depths.

Description: The climate of the last two millennia was characterised by decadal to multi‐centennial variations which were recorded in terrestrial records and had important societal impacts. The cause of these climatic events is still under debate but changes in the North Atlantic circulation have often been proposed to play an important role. In this review we compile available high‐resolution paleoceanographic datasets from the northern North Atlantic and Nordic Seas. The records are grouped into regions related to modern ocean conditions and their variability is discussed. We additionally discuss our current knowledge from modelling studies, with a specific focus on the dynamical changes that are not well inferred from the proxy records. An illustration is provided through the analysis of two climate model ensembles and an individual simulation of the last millennium. This review thereby provides an up‐to‐date paleo‐perspective on the North Atlantic multidecadal to multi‐centennial ocean variability across the last two millennia.

Description: In the course of the ongoing global intensification and diversification of human pressures, the study of variation patterns of biological traits along environmental gradients can provide relevant information on the performance of species under shifting conditions. The pronounced salinity gradient, co‐occurrence of multiple stressors, and accelerated rates of change make the Baltic Sea and its transition to North Sea a suitable region for this type of study. Focusing on the bladderwrack Fucus vesiculosus, one of the main foundation species on hard‐bottoms of the Baltic Sea, we analyzed the phenotypic variation among populations occurring along 2,000 km of coasts subjected to salinities from 4 to 〉30 and a variety of other stressors. Morphological and biochemical traits, including palatability for grazers, were recorded at 20 stations along the Baltic Sea and four stations in the North Sea. We evaluated in a common modeling framework the relative contribution of multiple environmental drivers to the observed trait patterns. Salinity was the main and, in some cases, the only environmental driver of the geographic trait variation in F. vesiculosus. The decrease in salinity from North Sea to Baltic Sea stations was accompanied by a decline in thallus size, photosynthetic pigments, and energy storage compounds, and affected the interaction of the alga with herbivores and epibiota. For some traits, drivers that vary locally such as wave exposure, light availability or nutrient enrichment were also important. The strong genetic population structure in this macroalgae might play a role in the generation and maintenance of phenotypic patterns across geographic scales. In light of our results, the desalination process projected for the Baltic Sea could have detrimental impacts on F. vesiculosus in areas close to its tolerance limit, affecting ecosystem functions such as habitat formation, primary production, and food supply.

Description: Key Points: • A regional ocean model is used to examine multidecadal shelf temperature changes on the Agulhas Bank • There are distinct shelf temperature regime changes in 1966 and 1996 • These regime shifts are caused by changes in coastal upwelling linked to large-scale wind variability The Agulhas Bank is an important area for the spawning of small pelagic fish and other species. Here, within a NEMO ocean model, we investigate changes in temperature over the Bank on multidecadal time scales. In agreement with previous observational studies, a shift to colder temperatures is found in 1997. The model also simulates an earlier shift from colder to warmer temperatures in 1966. These shifts are coastally confined and shown, using a climatologically forced model run as a control, to be driven by a north‐south migration in the large‐scale wind belts, rather than by changes in downward heat fluxes or changes in the Agulhas Current itself. The zonal wind changes on the Agulhas Bank show a significant relationship with the Southern Annular Mode, showing some promise for future predictability of cold and warm regimes on the Agulhas Bank. Thus, while the Agulhas Current has been shown in previous work to have a large impact on intra‐annual and interannual temperature variability, this work shows that multidecadal variability in temperature on the shelf is likely to be wind forced.

Description: Gravitational sinking of photosynthetically fixed particulate organic carbon (POC) constitutes a key component of the biological carbon pump. The fraction of POC leaving the surface ocean depends on POC sinking velocity (SV) and remineralization rate (Cremin), both of which depend on plankton community structure. However, the key drivers in plankton communities controlling SV and Cremin are poorly constrained. In fall 2014, we conducted a 6 weeks mesocosm experiment in the subtropical NE Atlantic Ocean to study the influence of plankton community structure on SV and Cremin. Oligotrophic conditions prevailed for the first 3 weeks, until nutrient‐rich deep water injected into all mesocosms stimulated diatom blooms. SV declined steadily over the course of the experiment due to decreasing CaCO3 ballast and – according to an optical proxy proposed herein – due to increasing aggregate porosity mostly during an aggregation event after the diatom bloom. Furthermore, SV was positively correlated with the contribution of picophytoplankton to the total phytoplankton biomass. Cremin was highest during a Synechococcus bloom under oligotrophic conditions and in some mesocosms during the diatom bloom after the deep‐water addition while it was particularly low during harmful algal blooms. The temporal changes were considerably larger in Cremin (max. 15‐fold) than in SV (max. 3‐fold). Accordingly, estimated POC transfer efficiency to 1000 m was mainly dependent on how the plankton community structure affected Cremin. Our approach revealed key players and interactions in the plankton food web influencing POC export efficiency thereby improving our mechanistic understanding of the biological carbon pump. Key points Sinking velocity was higher during oligotrophy than during blooms which is linked to ballast, porosity, and phytoplankton size structure Remineralization was highly variable but tended to be higher during Synechococcus or diatom blooms and lower during harmful algal blooms Plankton community structure had a considerably larger influence on particle remineralization rate than on sinking velocity

Description: The physiological processes driving the rapid rates of calcification in larval bivalves are poorly understood. Here, we use a calcification substrate-limited approach (low dissolved inorganic carbon, CT) and mRNA sequencing to identify proteins involved in bicarbonate acquisition during shell formation. As a secondary approach, we examined expression of ion transport and shell matrix proteins (SMPs) over the course of larval development and shell formation. We reared four families of Mytilus edulis under ambient (ca. 1865 µmol/kg) and low CT (ca. 941 µmol/kg) conditions and compared expression patterns at six developmental time points. Larvae reared under low CT exhibited a developmental delay, and a small subset of contigs was differentially regulated between ambient and low CT conditions. Of particular note was the identification of one contig encoding an anion transporter (SLC26) which was strongly upregulated (2.3–2.9 fold) under low CT conditions. By analyzing gene expression profiles over the course of larval development, we are able to isolate sequences encoding ion transport and SMPs to enhance our understanding of cellular pathways underlying larval calcification processes. In particular, we observe the differential expression of contigs encoding SLC4 family members (sodium bicarbonate cotransporters, anion exchangers), calcium-transporting ATPases, sodium/calcium exchangers, and SMPs such as nacrein, tyrosinase, and transcripts related to chitin production. With a range of candidate genes, this work identifies ion transport pathways in bivalve larvae and by applying comparative genomics to investigate temporal expression patterns, provides a foundation for further studies to functionally characterize the proteins involved in larval calcification.

Description: Mid‐ocean ridge axes are marked by segmentation of the axes and underlying magmatic systems. Fine‐scale segmentation has mainly been studied along fast‐spreading ridges. Here we offer insight into the third‐ and fourth‐order segmentation of intermediate‐spreading ridges and their temporal evolution. The Alarcón Rise and the Endeavour Segment have similar spreading rates (49 and 52.5 mm/year, respectively) but contrasting morphologies that vary from an axial high with a relatively narrow axial summit trough to an axial valley. One‐meter resolution bathymetry acquired by autonomous underwater vehicles, lava geochemistry, and ages from sediment cores is combined with available seismic reflection profiles to analyze variations in (1) geometry and orientation of the axial summit trough or valley, (2) seafloor depth near the axis, and (3) distribution of hydrothermal vents, (4) lava chemistry, and (5) flow ages between contiguous axes. Along both intermediate‐spreading segments, third‐ and fourth‐order discontinuities and associated segments are similar in dimension to what has been observed along fast‐spreading ridges. The Alarcón Rise and the Endeavour Segment also allow the study of the evolution of fine‐scale segmentation over periods of 300 to 4,000 years. Comparison between old and young axes reveals that the evolution of fine‐scale segmentation depends on the intensity of the magmatic activity. High magmatic periods are associated with rapid evolution of third‐order segments, while low magmatic activity periods, dominated by tectonic deformation and/or hydrothermal activity, are associated with little to no change in segmentation.

Description: The subterranean flow of water through sand barriers between coastal lagoons and the sea, driven by a positive hydraulic gradient, is a net new pathway for solute transfer to the sea. On the sea side of sand barriers, seawater circulation in the swash-zone generates a flux of recycled and new solutes. The significance and temporal variability of these vectors to the French Mediterranean Sea is unknown, despite lagoons constituting ~ 50% of the coastline. A one-dimensional Ra-224(ex)/Ra-223 reactive-transport model was used to quantify water flow between a coastal lagoon (La Palme) and the sea over a 6-month period. Horizontal flow between the lagoon and sea decreased from ~ 85 cm d(-1) during May 2017 (0.3 m(3) d(-1) m(-1) of shoreline) to ~ 20 cm d(-1) in July and was negligible in the summer months thereafter due to a decreasing hydraulic gradient. Seawater circulation in the swash-zone varied from 10 to 52 cm d(-1) (0.4-2.1 m(3) d(-1) m(-1)), driven by short-term changes in the prevailing wind and wave regimes. Both flow paths supply minor dissolved silica fluxes on the order of ~ 3-10 mmol Si d(-1) m(-1). Lagoon-sea water exchange supplies a net dissolved inorganic carbon (DIC) flux (320-1100 mmol C d(-1) m(-1)) two orders of magnitude greater than seawater circulation and may impact coastal ocean acidification. The subterranean flow of water through sand barriers represents a significant source of new DIC, and potentially other solutes, to the Mediterranean Sea during high lagoon water-level periods and should be considered in seasonal element budgets.

Description: We study the structure and tectonics of the collision zone between the Nazca Ridge (NR) and the Peruvian margin constrained by seismic, gravimetric, bathymetric, and natural seismological data. The NR was formed in an on-ridge setting, and it is characterized by a smooth and broad shallow seafloor (swell) with an estimated buoyancy flux of ~7 Mg/s. The seismic results show that the NR hosts an oceanic lower crust 10–14 km thick with velocities of 7.2–7.5 km/s suggesting intrusion of magmatic material from the hot spot plume to the oceanic plate. Our results show evidence for subduction erosion in the frontal part of the margin likely enhanced by the collision of the NR. The ridge-trench collision zone correlates with the presence of a prominent normal scarp, a narrow continental slope, and (uplifted) shelf. In contrast, adjacent of the collision zone, the slope does not present a topographic scarp and the continental slope and shelf become wider and deeper. Geophysical and geodetic evidence indicate that the collision zone is characterized by low seismic coupling at the plate interface. This is consistent with vigorous subduction erosion enhanced by the subducting NR causing abrasion and increase of fluid pore pressure at the interplate contact. Furthermore, the NR has behaved as a barrier for rupture propagation of megathrust earthquakes (e.g., 1746 Mw 8.6 and 1942 Mw 8.1 events). In contrast, for moderate earthquakes (e.g., 1996 Mw 7.7 and 2011 Mw 6.9 events), the NR has behaved as a seismic asperity nucleating at depths 〉20 km.

Description: Due to its remoteness, the deep-sea floor remains an understudied ecosystem of our planet. The patchiness of existing data sets makes it difficult to draw conclusions about processes that apply to a wider area. In our study we show how different settings and processes determine sediment heterogeneity on small spatial scales. We sampled solid phase and porewater from the upper 10 m of an approximately 7.4×13 km2 area in the Peru Basin, in the southeastern equatorial Pacific Ocean, at 4100 m water depth. Samples were analyzed for trace metals, including rare earth elements and yttrium (REY), as well as for particulate organic carbon (POC), CaCO3, and nitrate. The analyses revealed the surprisingly high spatial small-scale heterogeneity of the deep-sea sediment composition. While some cores have the typical green layer from Fe(II) in the clay minerals, this layer is missing in other cores, i.e., showing a tan color associated with more Fe(III) in the clay minerals. This is due to varying organic carbon contents: nitrate is depleted at 2–3 m depth in cores with higher total organic carbon contents but is present throughout cores with lower POC contents, thus inhibiting the Fe(III)-to-Fe(II) reduction pathway in organic matter degradation. REY show shale-normalized (SN) patterns similar to seawater, with a relative enrichment of heavy REY over light REY, positive LaSN anomaly, negative CeSN anomaly, and positive YSN anomaly and correlate with the Fe-rich clay layer and, in some cores, also correlate with P. We therefore propose that Fe-rich clay minerals, such as nontronite, as well as phosphates, are the REY-controlling phases in these sediments. Variability is also seen in dissolved Mn and Co concentrations between sites and within cores, which might be due to dissolving nodules in the suboxic sediment, as well as in concentration peaks of U, Mo, As, V, and Cu in two cores, which might be related to deposition of different material at lower-lying areas or precipitation due to shifting redox boundaries.

Description: Global climatologies of the seawater CO2 chemistry variables are necessary to assess the marine carbon cycle in depth. The climatologies should adequately capture seasonal variability to properly address ocean acidification and similar issues related to the carbon cycle. Total alkalinity (A(T)) is one variable of the seawater CO2 chemistry system involved in ocean acidification and frequently measured. We used the Global Ocean Data Analysis Project version 2.2019 (GLODAPv2) to extract relationships among the drivers of the A(T) variability and A(T) concentration using a neural network (NNGv2) to generate a monthly climatology. The GLODAPv2 quality-controlled dataset used was modeled by the NNGv2 with a root-mean-squared error (RMSE) of 5.3 mu mol kg(-1). Validation tests with independent datasets revealed the good generalization of the network. Data from five ocean time-series stations showed an acceptable RMSE range of 3-6.2 mu mol kg(-1). Successful modeling of the monthly A(T) variability in the time series suggests that the NNGv2 is a good candidate to generate a monthly climatology. The climatological fields of A(T) were obtained passing through the NNGv2 the World Ocean Atlas 2013 (WOA13) monthly climatologies of temperature, salinity, and oxygen and the computed climatologies of nutrients from the previous ones with a neural network. The spatiotemporal resolution is set by WOA13: 1 degrees x 1 degrees in the horizontal, 102 depth levels (0-5500 m) in the vertical and monthly (0-1500 m) to annual (1550-5500 m) temporal resolution. The product is distributed through the data repository of the Spanish National Research Council (CSIC; https://doi.org/10.20350/digitalCSIC/8644, Broullon et al., 2019).

Description: Does disease resistance evolution in vitro reflect resistance evolution in vivo? Hernandez and Koskella conducted serial passage experiments of the plant pathogenic bacterium Pseudomonas syringae and two lytic bacteriophages in high‐nutrient medium (in vitro) and in a tomato plant (in vivo). High levels of bacterial resistance to phages evolved in vitro but not in vivo, suggesting that high costs and low benefits of resistance explain the observed pattern.

Description: Rationale: Nitrous oxide (N2O) is an atmospheric trace gas regulating Earth's climate, and is a key intermediate of many nitrogen cycling processes in aquatic ecosystems. Laser-based technology for N2O concentration and isotopic/isotopomeric analyses has potential advantages, which include high analytical specificity, low sample size requirement and reduced cost. Methods: An autosampler with a purge-and-trap module is coupled to a cavity ring-down spectrometer to achieve automated and high-throughput measurements of N2O concentrations, N2O isotope ratios (δ15Nbulk and δ18O values) and position-specific isotopomer ratios (δ15Nα and δ15Nβ values). The system provides accuracy and precision similar to those for measurements made by traditional isotope ratio mass spectrometry (IRMS) techniques. Results: The sample sizes required were 0.01–1.1 nmol-N2O. Measurements of four N2O isotopic/isotopomeric references were cross-calibrated with those obtained by IRMS. With a sample size of 0.50 nmol-N2O, the measurement precision (1σ) for δ15Nα, δ15Nβ, δ15Nbulk and δ18O values was 0.61, 0.33, 0.41 and 0.43‰, respectively. Correction schemes were developed for sample size-dependent isotopic/isotopomeric deviations. The instrumental system demonstrated consistent measurements of dissolved N2O concentrations, isotope/isotopomer ratios and production rates in seawater. Conclusions: The coupling of an autosampler with a purge-and-trap module to a cavity ring-down spectrometer not only significantly reduces sample size requirements, but also offers comprehensive investigation of N2O production pathways by the measurement of natural abundance and tracer level isotopes and isotopomers. Furthermore, the system can perform isotopic analyses of dissolved and solid nitrogen-containing samples using N2O as the analytical proxy.

Description: Distinct differences were observed in geochemical signatures in sediments from two sites drilled in the upper plate of the Costa Rica margin during Integrated Ocean Drilling Program (IODP) Expedition 334. The upper 80 m at Site U1379, located on the outer shelf, show pore water non‐steady state conditions characteristic of a declining methane flux. These contrast with analyses of the upper sediment layers at the middle slope site (U1378) that reflect steady state conditions. Distinct carbonate‐rich horizons up to 11 meters thick were recovered between 63 and 310 meters below seafloor at Site U1379 but were not found at Site U1378. The carbonates and dissolved inorganic carbon from Site U1379 have a depleted carbon stable isotope signal (up to ‐25‰) that indicates anaerobic methane oxidation. This inference is further supported by distinct δ34S‐pyrite and magnetic susceptibility records that reveal fluctuations of the sulfate‐methane transition in response to methane flux variations. Tectonic reconstructions of this margin document a marked subsidence event after arrival of the Cocos Ridge, 2.2 ± 0.2 million years ago (Ma), followed by increased sedimentation rates and uplift. As the seafloor at Site U1379 rose from ~2000 m to the present water depth of ~126 m, the site moved out of the gas hydrate stability zone (GHSZ) at ~1.1 Ma, triggering upward methane advection, methane oxidation, and the onset of massive carbonate formation. Younger carbonate occurrences and the non‐steady state pore profiles at Site U1379 reflect continued episodic venting likely modulated by changes in the underlying methane reservoir.

Description: The flow (flux) of climate-critical gases, such as carbon dioxide (CO2), between the ocean and the atmosphere is a fundamental component of our climate and an important driver of the biogeochemical systems within the oceans. Therefore, the accurate calculation of these air–sea gas fluxes is critical if we are to monitor the oceans and assess the impact that these gases are having on Earth's climate and ecosystems. FluxEngine is an open-source software toolbox that allows users to easily perform calculations of air–sea gas fluxes from model, in situ, and Earth observation data. The original development and verification of the toolbox was described in a previous publication. The toolbox has now been considerably updated to allow for its use as a Python library, to enable simplified installation, to ensure verification of its installation, to enable the handling of multiple sparingly soluble gases, and to enable the greatly expanded functionality for supporting in situ dataset analyses. This new functionality for supporting in situ analyses includes user-defined grids, time periods and projections, the ability to reanalyse in situ CO2 data to a common temperature dataset, and the ability to easily calculate gas fluxes using in situ data from drifting buoys, fixed moorings, and research cruises. Here we describe these new capabilities and demonstrate their application through illustrative case studies. The first case study demonstrates the workflow for accurately calculating CO2 fluxes using in situ data from four research cruises from the Surface Ocean CO2 ATlas (SOCAT) database. The second case study calculates air–sea CO2 fluxes using in situ data from a fixed monitoring station in the Baltic Sea. The third case study focuses on nitrous oxide (N2O) and, through a user-defined gas transfer parameterisation, identifies that biological surfactants in the North Atlantic could suppress individual N2O sea–air gas fluxes by up to 13 %. The fourth and final case study illustrates how a dissipation-based gas transfer parameterisation can be implemented and used. The updated version of the toolbox (version 3) and all documentation is now freely available.

Description: Dissolved organic carbon (DOC) plays critical roles in marine carbon cycling, but its sources and sinks remain uncertain. In this study, we monitored DOC exudation rates of Sargassum natans under visible light (lambda 〉 390 nm) and solar radiation. DOC release rates ranged from 7 to 10 mu g C g(biomass)(-1) hr(-1) (wet weight) under visible light, but increased to 23 to 41 mu g C g(biomass)(-1) hr(-1) when exposed to natural sunlight. Results indicate that DOC released by Sargassum could amount to 0.3 to 1.2 Tg C/year, potentially contributing significantly to the marine DOC pool in the Gulf of Mexico and Western North Atlantic. We employed the Folin-Ciocalteu phenolic content method, nuclear magnetic resonance (NMR) spectroscopy, and ultrahigh resolution Fourier transform ion cyclotron resonance mass spectrometry (FT-ICR MS) to characterize the diverse pool of organic compounds exuded from Sargassum. Results from these complementary methods showed that Sargassum release large quantities of phlorotannins, a class of polyphenols that have very similar properties to terrestrial DOC. These phlorotannins and their oxygenated phenolic derivatives exhibit a high hydrogen deficiency and functionalization (i.e., 4 to 6 oxygen atoms per aromatic ring), representing 5 to 18% of the released DOC isolated by solid phase extraction. Thus, Sargassum is the largest biological source of open ocean polyphenols recorded to date. The amount of polyphenolic DOC released by Sargassum challenges previous beliefs that all polyphenols found within the oceans are remnants of terrestrial organic matter, although the stability of phlorotannins and their derivatives needs to be further evaluated.

Description: We present a three-dimensional gas hydrate systems model of the southern Hikurangi subduction margin in eastern New Zealand. The model integrates thermal and microbial gas generation, migration, and hydrate formation. Modeling these processes has improved the understanding of factors controlling hydrate distribution. Three spatial trends of concentrated hydrate occurrence are predicted. The first trend (I) is aligned with the principal deformation front in the overriding Australian plate. Concentrated hydrate deposits are predicted at or near the apexes of anticlines and to be mainly sourced from focused migration and recycling of microbial gas generated beneath the hydrate stability zone. A second predicted trend (II) is related to deformation in the subducting Pacific plate associated with former Mesozoic subduction beneath Gondwana and the modern Pacific-Australian plate boundary. This trend is enhanced by increased advection of thermogenic gas through permeable layers in the subducting plate and focused migration into the Neogene basin fill above Cretaceous-Paleogene structures. The third trend (III) follows the northern margin of the Hikurangi Channel and is related to the presence of buried strata of the Hikurangi Channel system. The predicted trends are consistent with pronounced seismic reflection anomalies related to free gas in the pore space and strength of the bottom-simulating reflection. However, only trend I is also associated with clear and widespread seismic indications of concentrated gas hydrate. Total predicted hydrate masses at the southern Hikurangi Margin are between 52,800 and 69,800 Mt. This equates to 3.4–4.5 Mt hydrate/km2, containing 6.33 × 108–8.38 × 108 m3/km2 of methane.

Description: Global climatic changes during the last Glacial and Deglacial have been related to variations of the Atlantic Meridional Overturning Circulation (AMOC). Here, we present new and refined 231Pa/230Th down‐core profiles extending back to 30 ka BP from the northwestern Atlantic along the Atlantic Deep Western Boundary Current (DWBC), which is the main component of the southward deep backflow of the AMOC. Besides the well‐known Bermuda Rise records, available high‐resolution 231Pa/230Th data in the northwestern Atlantic are still sparse. Our new records along with reconstructions of deep water provenance from Nd isotopes constrain the timing and magnitude of past changes in AMOC from an additional northwestern Atlantic region forming a depth transect between 3000 and 4760 m water depth. Our extended and improved dataset confirms the weakening of the AMOC during deglacial cold spells such as Heinrich Event 1 and the Younger Dryas interrupted by a reinvigoration during the Bølling‐Allerød interstadial as seen in the prominent 231Pa/230Th records from the Bermuda Rise. However, in contrast to the Bermuda Rise records we find a clearly reduced circulation strength during the Last Glacial Maximum in the deep Atlantic.

Description: Gelatinous zooplankton (Cnidaria, Ctenophora, and Urochordata, namely, Thaliacea) are ubiquitous members of plankton communities linking primary production to higher trophic levels and the deep ocean by serving as food and transferring “jelly‐carbon” (jelly‐C) upon bloom collapse. Global biomass within the upper 200 m reaches 0.038 Pg C, which, with a 2–12 months life span, serves as the lower limit for annual jelly‐C production. Using over 90,000 data points from 1934 to 2011 from the Jellyfish Database Initiative as an indication of global biomass (JeDI: http://jedi.nceas.ucsb.edu, http://www.bco‐dmo.org/dataset/526852), upper ocean jelly‐C biomass and production estimates, organism vertical migration, jelly‐C sinking rates, and water column temperature profiles from GLODAPv2, we quantitatively estimate jelly‐C transfer efficiency based on Longhurst Provinces. From the upper 200 m production estimate of 0.038 Pg C year−1, 59–72% reaches 500 m, 46–54% reaches 1,000 m, 43–48% reaches 2,000 m, 32–40% reaches 3,000 m, and 25–33% reaches 4,500 m. This translates into ~0.03, 0.02, 0.01, and 0.01 Pg C year−1, transferred down to 500, 1,000, 2,000, and 4,500 m, respectively. Jelly‐C fluxes and transfer efficiencies can occasionally exceed phytodetrital‐based sediment trap estimates in localized open ocean and continental shelves areas under large gelatinous blooms or jelly‐C mass deposition events, but this remains ephemeral and transient in nature. This transfer of fast and permanently exported carbon reaching the ocean interior via jelly‐C constitutes an important component of the global biological soft‐tissue pump, and should be addressed in ocean biogeochemical models, in particular, at the local and regional scale.

Description: Sea surface salinity (SSS) is an important variable in the global ocean circulation. However, decadal to interdecadal changes in SSS are not well understood due to the lack of instrumental data. Here we reconstruct SSS from a paired, bimonthly resolved coral δ18O and Sr/Ca record from La Reunion Island that extends from 1913 to 1995. Coral Sr/Ca correlates with regional sea surface temperature (SST) back to 1966, when instrumental coverage is good, while coral δ18O does not. The slope of the monthly (annual mean) coral Sr/Ca-SST regression is −0.040 mmol/mol per 1 °C (−0.068 mmol/mol per 1 °C) consistent with published estimates of the Sr/Ca-SST relationship. Coral Sr/Ca suggest a warming of 0.39 °C since 1913. δ18O seawater is calculated by subtracting the temperature component from measured coral δ18O, using coral Sr/Ca as well as historical SST products. The derived δ18O seawater reconstructions are correlated (r 〉 0.6), and all show a significant shift in the midtwentieth century (−0.17‰ to −0.19‰), indicating a freshening of SSS by 0.7 psu. However, the timing of this shift depends on the temperature component and varies from 1947 (δ18O seawater calculated with historical SST) to the late 1950s (δ18O seawater calculated with coral Sr/Ca). Coral Sr/Ca shows warm temperature anomalies in the mid-1950s, while historical SST products show warm anomalies from 1940 to 1945 followed by cooling in the 1950s, a pattern typical for the World War II bias. This suggests that historical SST may bias reconstructions of δ18O seawater and SSS from corals.

Description: Climate reconstructions are means to extract the signal from uncertain paleo-observations, so-called proxies. It is essential to evaluate these reconstructions to understand and quantify their uncertainties. Similarly, comparing climate simulations and proxies requires approaches to bridge the temporal and spatial differences between both and to address their specific uncertainties. One way to achieve these two goals is so-called pseudoproxies. These are surrogate proxy records within the virtual reality of a climate simulation. They in turn depend on an understanding of the uncertainties of the real proxies including the noise characteristics disturbing the original environmental signal. Common pseudoproxy approaches so far concentrate on data with high temporal resolution over the last approximately 2000 years. Here we provide a simple but flexible noise model for potentially low-resolution sedimentary climate proxies for temperature on millennial timescales, the code for calculating a set of pseudoproxies from a simulation, and one example of pseudoproxies. The noise model considers the influence of other environmental variables, a dependence on the climate state, a bias due to changing seasonality, modifications of the archive (for example bioturbation), potential sampling variability, and a measurement error. Model, code, and data allow us to develop new ways of comparing simulation data with proxies on long timescales. Code and data are available at https://doi.org/10.17605/OSF.IO/ZBEHX (Bothe et al., 2018).

Description: Atmospheric deposition is an important source of micronutrients to the ocean, but atmospheric deposition fluxes remain poorly constrained in most ocean regions due to the limited number of field observations of wet and dry atmospheric inputs. Here we present the distribution of dissolved aluminium (dAl), as a tracer of atmospheric inputs, in surface waters of the Atlantic Ocean along GEOTRACES sections GA01, GA06, GA08, and GA10. We used the surface mixed layer concentrations of dAl to calculate atmospheric deposition fluxes using a simple steady state model. We have optimized the aerosol Al fractional solubility, dAl residence time within the surface mixed layer and depth of the surface mixed layer for each separate cruise to calculate the atmospheric deposition fluxes. We calculated the lowest deposition fluxes of 0.15 ± 0.1 and 0.27 ± 0.13 g m−2 yr−1 for the South and North Atlantic Ocean (〉 40° S and 〉 40° N), respectively, and highest fluxes of 2.67 ± 1.96 and 3.82 ± 2.72 g m−2 yr−1 for the South East Atlantic and tropical Atlantic Ocean, respectively. Overall, our estimations are comparable to atmospheric dust deposition model estimates and reported field-based atmospheric deposition estimates. We note that our estimates diverge from atmospheric dust deposition model flux estimates in regions influenced by riverine Al inputs and in upwelling regions. As dAl is a key trace element in the GEOTRACES Programme, the approach presented in this study allows calculations of atmospheric deposition fluxes at high spatial resolution for remote ocean regions.

Description: Aim: Numerous regions worldwide are highly impacted by anthropogenic activities and globalization, with climate change and species introductions being among the greatest stressors to biodiversity and ecosystems. A main donor region of non‐indigenous species (NIS) for numerous European water bodies, as well as in the North American Great Lakes is the Ponto‐Caspian region (i.e., Black, Azov and Caspian Seas), with some of those species having significant impact on local communities and ecosystem functioning. Location: Northern European, Ponto‐Caspian and North American regions. Methods: To determine environmental tolerance of native species and related NIS under current and future global warming scenarios of the Baltic Sea, we conducted common garden experiments to test temperature tolerance of three euryhaline gammarid species: one Baltic (Gammarus oceanicus), one Ponto‐Caspian (Pontogammarus maeoticus) and one North American species (Gammarus tigrinus) in two different salinities. Results: Our results determined that mortality of P. maeoticus in all temperature treatments (i.e., increased, control, and decreased) at the end of both experiments (i.e., conducted in salinities of 10 and 16 g/kg) was lower when compared to mortality of G. oceanicus and (c) G. tigrinus. The highest mortality was observed for G. oceanicus, reaching 100% in both experiments in the increased temperature treatment. Main conclusions: Due to the high environmental tolerance of the Ponto‐Caspian species tested in this study, as well as the fact that Ponto‐Caspian species evolved in environmentally variable habitats and currently inhabit warmer waters than species from North America and Northern Europe, we suggest that species from the Ponto‐Caspian region may benefit from global warming when invading new areas. Those new invasions may, in the best case scenario, increase biodiversity of the Baltic Sea. However, if notorious invaders arrive, they may have a significant impact on local communities and ecosystem functioning.

Description: Reconstruction of early Cenozoic deep‐water circulation is one of the keys to modelling Earth's greenhouse‐to‐icehouse surface evolution, but it has long been hampered by the paucity of information from the central South Pacific. To help overcome this knowledge gap, we present new micropalaeontological data from dredged carbonates (R/V Sonne Expedition SO193) at several eastern volcanic salients of the Manihiki Plateau. Interestingly, despite appreciable longitudinal separations among the dredged sites, ages indicated by the foraminiferal assemblages are consistently around the Middle Eocene [including mixed Turonian (Late Cretaceous)/Eocene at a single site], suggesting widespread post‐Eocene cessation of the pelagic sedimentation. By integrating with independent seismic and chronostratigraphic data (Deep Sea Drilling Project Leg 33) for large‐scale erosion of top‐Eocene–Oligocene sedimentary units on the eastern Manihiki Plateau, our results can be viewed as novel physical evidence for the intensification of central South Pacific deep‐water circulation since the Eocene/Oligocene climatic transition.

Description: The nature of the Ionian Sea crust has been the subject of scientific debate for more than 30 years, mainly because seismic imaging of the deep crust and upper mantle of the Ionian Abyssal Plain (IAP) has not been conclusive to date. The IAP is sandwiched between the Calabrian and Hellenic subduction zones in the central Mediterranean. To univocally confirm the proposed oceanic nature of the IAP crust as a remnant of the Tethys ocean and to confute its interpretation as a strongly thinned part of the African continental crust, a NE-SW oriented 131 km long seismic refraction and wide-angle reflection profile consisting of eight ocean bottom seismometers and hydrophones was acquired in 2014. A P-wave velocity model developed from travel time forward modelling is refined by gravimetric data and synthetic modelling of the seismic data. A roughly 6km thick crust with velocities ranging from 5.1km/s to 7.2km/s, top to bottom, can be traced throughout the IAP. In the vicinity of the Medina Seamounts at the southern IAP boundary, the crust thickens to about 9km and seismic velocities decrease to 6.8km/s at the crust-mantle boundary. The seismic velocity distribution and depth of the crust-mantle boundary in the IAP document its oceanic nature, and support the interpretation of the IAP as a remnant of the Tethys oceanic lithosphere formed during the Permian and Triassic period.

Description: We describe the development of the “Paleoclimate PLASIM-GENIE emulator” PALEO-PGEM and its application to derive a downscaled high-resolution spatiotemporal description of the climate of the last five million years. The 5-million-year time frame is interesting for a range of paleo-environmental questions, not least because it encompasses the evolution of humans. However, the choice of time-frame was primarily pragmatic; tectonic changes can be neglected to first order, so that it is reasonable to consider climate forcing restricted to the Earth's orbital configuration, ice-sheet state and the concentration of atmosphere CO2. The approach uses the Gaussian process emulation of the singular value decomposition of boundary-condition ensembles of the intermediate complexity atmosphere-ocean GCM PLASIM-GENIE. Spatial fields of bioclimatic variables of surface air temperature (warmest and coolest seasons) and precipitation (wettest and driest seasons) are emulated at 1,000 year intervals, driven by time-series of scalar boundary-condition forcing (CO2, orbit and ice-volume), and assuming the climate is in quasi-equilibrium. Paleoclimate anomalies at climate model resolution are interpolated onto the observed modern climatology to produce a high-resolution spatiotemporal paleoclimate reconstruction of the Pliocene-Pleistocene.

Description: The interaction between the atmosphere, specifically the North Atlantic Oscillation (NAO), and the North Atlantic ocean circulation on sub‐decadal timescale is analyzed in a subset of models participating in the Coupled Model Intercomparison Project phase 5 (CMIP5). From preindustrial control runs of at least 500 years length, we derive anomaly patterns in the atmospheric and ocean circulation and of air‐sea heat exchange. All models simulate a distinct dipolar oceanic overturning anomaly at the sub‐decadal timescale, with centers at 30° N and 55° N. The dipolar overturning anomaly goes along with marked anomalies in the North Atlantic sea surface temperature and gyre circulation. Lag‐regression analyses demonstrate, with relatively small ensemble spread, how the atmosphere and the ocean circulation interact. The dipolar anomalies in the overturning are forced by NAO‐related wind stress curl anomalies. Anomalous surface heat fluxes in concert with anomalous vertical motions drive a meridional dipolar heat content anomaly in the upper ocean, and it is this dipolar heat content anomaly which carries the coupled system from one phase of the sub‐decadal cycle to the other by reversing the tendencies in the overturning circulation. The coupled sub‐decadal variability derived from the CMIP5 models is characterized by three elements: a wind‐driven part steering the dipolar overturning anomaly, surface heat flux anomalies that support a heat build‐up in the subpolar gyre region, and the heat storage memory which is instrumental in the phase reversal of the NAO.

Description: Long-term measurements of volcanic gas emissions conducted during the recent decade suggest that under certain conditions the magnitude or chemical composition of volcanic emissions exhibits periodic variations with a period of about two weeks. A possible cause of such a periodicity can be attributed to the Earth tidal potential. The phenomenology of such a link has been debated for long, but no quantitative model has yet been proposed. The aim of this paper is to elucidate whether a causal link from the tidal forcing to variation in the volcanic degassing can be traced analytically. We model the response of a simplified magmatic system to the local tidal gravity variations and derive a periodical vertical magma displacement in the conduit with an amplitude of 0.1–1 m, depending on geometry and physical state of the magmatic system. We find that while the tide-induced vertical magma displacement has presumably no significant direct effect on the volatile solubility, the differential magma flow across the radial conduit profile may result in a significant increase of the bubble coalescence rate in a depth of several kilometres by up to several ten percent. Because bubble coalescence facilitates separation of gas from magma and thus enhances volatile degassing, we argue that the derived tidal variation may propagate to a manifestation of varying volcanic degassing behaviour. The presented model provides a first basic framework which establishes an analytical understanding of the link between the Earth tides and volcanic degassing.

Description: Ground-based atmospheric observations of CO2, δ(O2∕N2), N2O, and CH4 were used to make estimates of the air–sea fluxes of these species from the Lüderitz and Walvis Bay upwelling cells in the northern Benguela region, during upwelling events. Average flux densities (±1σ) were 0.65±0.4 µmol m−2 s−1 for CO2, −5.1±2.5 µmol m−2 s−1 for O2 (as APO), 0.61±0.5 nmol m−2 s−1 for N2O, and 4.8±6.3 nmol m−2 s−1 for CH4. A comparison of our top-down (i.e., inferred from atmospheric anomalies) flux estimates with shipboard-based measurements showed that the two approaches agreed within ±55 % on average, though the degree of agreement varied by species and was best for CO2. Since the top-down method overestimated the flux density relative to the shipboard-based approach for all species, we also present flux density estimates that have been tuned to best match the shipboard fluxes. During the study, upwelling events were sources of CO2, N2O, and CH4 to the atmosphere. N2O fluxes were fairly low, in accordance with previous work suggesting that the evasion of this gas from the Benguela is smaller than for other eastern boundary upwelling systems (EBUS). Conversely, CH4 release was quite high for the marine environment, a result that supports studies that indicated a large sedimentary source of CH4 in the Walvis Bay area. These results demonstrate the suitability of atmospheric time series for characterizing the temporal variability of upwelling events and their influence on the overall marine greenhouse gas (GHG) emissions from the northern Benguela region.

Description: Rationale Boron (B) is an essential micronutrient in plants and its isotope variations are used to gain insights into plant metabolism, which is important for crop plant cultivation. B isotope variations were used to trace intraplant fractionation mechanisms in response to the B concentration in the irrigation water spanning the range from B depletion to toxic levels. Methods A fully validated analytical procedure based on MC‐ICP‐MS, sample decomposition and B matrix separation was applied to study B isotope fractionation. The validation was accomplished by establishing a complete uncertainty budget and by applying reference materials, yielding expanded measurement uncertainties of 0.8 ‰ for pure boric acid solutions and ≤ 1.5 ‰ for processed samples. With this validated procedure SI traceable B isotope amount ratios were determined in plant reference materials for the first time. Results The B isotope compositions of irrigation water and bell pepper samples suggest passive diffusion of the heavy 11B isotope into the roots during low to high B concentrations while uptake of the light 10B isotope was promoted during B depletion, probably by active processes. A systematic enrichment of the heavy 11B isotope in higher located plant parts was observed (average Δ11Bleaf‐roots = 20.3 ± 2.8 ‰ (1 SD)), possibly by a facilitated transport of the heavy 11B to growing meristems by B transporters. Conclusions B isotopes can be used to identify plant metabolism in response to the B concentration in the irrigation water and during intraplant B transfer. The large B isotope fractionation within the plants demonstrates the importance of biological B cycling for the global B cycle.

Description: Recent studies using reanalysis data and complex models suggest that the Tropics influence midlatitude blocking. Here, the influence of tropical precipitation anomalies is investigated further using a dry dynamical model driven by specified diabatic heating anomalies. The model uses a quasi‐realistic setup based on idealized orography and an idealized representation of the land‐ocean thermal contrast. Results concerning the El Niño Southern Oscillation and the Madden‐Julian Oscillation are mostly consistent with previous studies and emphasize the importance of tropical dynamics for driving the variability of blocking at midlatitudes. It is also shown that a common bias in models of the Coupled Model Intercomparison Project Phase 5 (CMIP5), namely, excessive tropical precipitation, leads to an underestimation of midlatitude blocking in our model, also a common bias in the CMIP5 models. The strongest blocking anomalies associated with the tropical precipitation bias are found over Europe, where the underestimation of blocking in CMIP5 models is also particularly strong.

Description: Key Points: - Rapid subsurface oceanographic change in the tropical W Atlantic reflect shifting Subtropical Gyre - Subsurface warming responds to deglacial AMOC perturbations (Heinrich Stadials 2, 1, and the Younger Dryas) - Southward propagation of Salinity Maximum Water during Northern Hemisphere cold spells shift the mixing zone of tropical and subtropical waters During times of deglacial Atlantic Meridional Overturning Circulation (AMOC) perturbations, the tropical Atlantic experienced considerable warming at subsurface levels. Coupled ocean‐atmosphere simulations corroborate the tight teleconnection between the tropical Atlantic and climate change at high northern latitudes, but still underestimate the relevance of the subsurface N Atlantic Subtropical Gyre (STG) for heat and salt storage and its sensitivity to rapid climatic change. We here reconstruct vertical and lateral temperature and salinity gradients in the tropical W Atlantic and the Caribbean over the last 30 kyrs, based on planktic deep and shallow dwelling foraminiferal Mg/Ca and δ18O‐records. The rapid and large amplitude subsurface changes illustrate a dynamic STG associated with abrupt shifts of North Atlantic hydrographic and atmospheric regimes. During full glacial conditions, the STG has been shifted southward while intensified Ekman‐downwelling associated to strengthened trade winds fostered the formation of warm and saline Salinity Maximum Water (SMW). The southward propagation of SMW was facilitated by the glacially eastward deflected North Brazil Current. During periods of significant AMOC perturbations (Heinrich Stadials 1, and the Younger Dryas), extreme subsurface warming by ~6°C led to diminished lateral subsurface temperature gradients. Coevally, a deep thermocline suggests that SMW fully occupied the subsurface tropical W Atlantic and that the STG reached its southernmost position. During the Holocene, modern‐like conditions gradually developed with the northward retreat of SMW and the development of a strong thermocline ridge between the Subtropical Gyre and the tropical W Atlantic.

Description: Quaternary East Asian winter monsoon (EAWM) evolution has long been attributed to high-latitude Northern Hemisphere climate change. However, it cannot explain the distinct relationships of the EAWM in the northern and southern East Asian marginal sea in paleoclimatic records. Here we present an EAWM record of the northern East China Sea over the past 300 ka and a transient climate simulation with the Kiel Climate Model through the Holocene. Both proxy record and simulation suggest anticorrelated long-term EAWM evolution between the northern East China Sea and the South China Sea. We suggest that this spatial discrepancy of EAWM can be interpreted as El Niño–Southern Oscillation (ENSO)-like controlling, which generates cyclonic/anticyclonic wind anomalies in the northern/southern East Asian marginal sea. This research explains much of the controversy in nonorbital scale variability of Quaternary EAWM records in the East Asian marginal sea and supports a potent role of tropical forcing in East Asian winter climate change.

Description: The atmospheric CO2 concentration increased by about 20ppm from 6000BCE to the pre-industrial period (1850CE). Several hypotheses have been proposed to explain mechanisms of this CO2 growth based on either ocean or land carbon sources. Here, we apply the Earth system model MPI-ESM-LR for two transient simulations of climate and carbon cycle dynamics during this period. In the first simulation, atmospheric CO2 is prescribed following ice-core CO2 data. In response to the growing atmospheric CO2 concentration, land carbon storage increases until 2000BCE, stagnates afterwards, and decreases from 1CE, while the ocean continuously takes CO2 out of the atmosphere after 4000BCE. This leads to a missing source of 166Pg of carbon in the ocean-land-atmosphere system by the end of the simulation. In the second experiment, we applied a CO2 nudging technique using surface alkalinity forcing to follow the reconstructed CO2 concentration while keeping the carbon cycle interactive. In that case the ocean is a source of CO2 from 6000 to 2000BCE due to a decrease in the surface ocean alkalinity. In the prescribed CO2 simulation, surface alkalinity declines as well. However, it is not sufficient to turn the ocean into a CO2 source. The carbonate ion concentration in the deep Atlantic decreases in both the prescribed and the interactive CO2 simulations, while the magnitude of the decrease in the prescribed CO2 experiment is underestimated in comparison with available proxies. As the land serves as a carbon sink until 2000BCE due to natural carbon cycle processes in both experiments, the missing source of carbon for land and atmosphere can only be attributed to the ocean. Within our model framework, an additional mechanism, such as surface alkalinity decrease, for example due to unaccounted for carbonate accumulation processes on shelves, is required for consistency with ice-core CO2 data. Consequently, our simulations support the hypothesis that the ocean was a source of CO2 until the late Holocene when anthropogenic CO2 sources started to affect atmospheric CO2.

Description: New marine geophysical data acquired across the partly ice‐covered northern East Greenland continental margin highlight a complex interaction between tectonic and magmatic events. Breakup‐related lava flows are imaged in reflection seismic data as seaward dipping reflectors (SDRs), which are found to decrease in size both northwards and southwards from a central point at 75° N. We provide evidence that the magnetic anomaly pattern in the shelf area is related to volcanic phases and not to the presence of oceanic crust. The remnant magnetization of the individual lava flows is used to deduce a relative timing of the emplacement of the volcanic wedges. We find that the SDRs have been emplaced over a period of 2‐4 Ma progressively from north to south and from landward to seaward. The new data indicate a major post‐middle Eocene magmatic phase around the landward termination of the West Jan Mayen Fracture Zone. This post‐40 Ma volcanism likely was associated with the progressive separation of the Jan Mayen microcontinent from East Greenland. The break‐up of the Greenland Sea started at several isolated seafloor spreading cells whose location was controlled by rift structures and led to the present‐day segmentation of the margin. The original rift basins were subsequently connected by steady‐state seafloor spreading that propagated southwards, from the Greenland Fracture Zone to the Jan Mayen Fracture Zone. Key Points Polyphase Cenozoic volcanic rifting and consecutive emplacement of breakup‐related lava flows units along the northern East Greenland margin Breakup along restricted margin segments is followed by north to south directed progressive opening of the Greenland Sea Widespread post‐middle Eocene (〈 40 Ma) offshore magmatism, associated with the breakup of the Jan Mayen microcontinent from East Greenland

Description: The stable isotopic composition of particulate organic carbon (δ13CPOC) in the surface waters of the global ocean can vary with the aqueous CO2 concentration ([CO2(aq)]) and affects the trophic transfer of carbon isotopes in the marine food web. Other factors such as cell size, growth rate and carbon concentrating mechanisms decouple this observed correlation. Here, the variability in δ13CPOC is investigated in surface waters across the south subtropical convergence (SSTC) in the Atlantic Ocean, to determine carbon isotope fractionation (ϵp) by phytoplankton and the contrasting mechanisms of carbon uptake in the subantarctic and subtropical water masses. Our results indicate that cell size is the primary determinant of δ13CPOC across the Atlantic SSTC in summer. Combining cell size estimates with CO2 concentrations, we can accurately estimate "p within the varying surface water masses in this region. We further utilize these results to investigate future changes in "p with increased anthropogenic carbon availability. Our results suggest that smaller cells, which are prevalent in the subtropical ocean, will respond less to increased [CO2(aq)] than the larger cells found south of the SSTC and in the wider Southern Ocean. In the subantarctic water masses, isotopic fractionation during carbon uptake will likely increase, both with increasing CO2 availability to the cell, but also if increased stratification leads to decreases in average community cell size. Coupled with decreasing δ13C of [CO2(aq)] due to anthropogenic CO2 emissions, this change in isotopic fractionation and lowering of δ13CPOC may propagate through the marine food web, with implications for the use of δ13CPOC as a tracer of dietary sources in the marine environment.

Description: The dense AlpArray network allows studying seismic wave propagation with high spatial resolution. Here we introduce an array approach to measure arrival angles of teleseismic Rayleigh waves. The approach combines the advantages of phase correlation as in the two-station method with array beamforming to obtain the phase-velocity vector. 20 earthquakes from the first two years of the AlpArray project are selected, and spatial patterns of arrival-angle deviations across the AlpArray are shown in maps, depending on period and earthquake location. The cause of these intriguing spatial patterns is discussed. A simple wave-propagation modelling example using an isolated anomaly and a Gaussian beam solution suggests that much of the complexity can be explained as a result of wave interference after passing a structural anomaly along the wave paths. This indicates that arrival-angle information constitutes useful additional information on the Earth structure, beyond what is currently used in inversions.

Description: The amount of additional future temperature change following a complete cessation of CO2 emissions is a measure of the unrealized warming to which we are committed due to CO2 already emitted to the atmosphere. This "zero emissions commitment" (ZEC) is also an important quantity when estimating the remaining carbon budget - a limit on the total amount of CO2 emissions consistent with limiting global mean temperature at a particular level. In the recent IPCC Special Report on Global Warming of 1.5 degrees C, the carbon budget framework used to calculate the remaining carbon budget for 1.5 degrees C included the assumption that the ZEC due to CO2 emissions is negligible and close to zero. Previous research has shown significant uncertainty even in the sign of the ZEC. To close this knowledge gap, we propose the Zero Emissions Commitment Model Intercomparison Project (ZECMIP), which will quantify the amount of unrealized temperature change that occurs after CO2 emissions cease and investigate the geophysical drivers behind this climate response. Quantitative information on ZEC is a key gap in our knowledge, and one that will not be addressed by currently planned CMIP6 simulations, yet it is crucial for verifying whether carbon budgets need to be adjusted to account for any unrealized temperature change resulting from past CO2 emissions. We request only one top-priority simulation from comprehensive general circulation Earth system models (ESMs) and Earth system models of intermediate complexity (EMICs) - a branch from the 1% CO2 run with CO2 emissions set to zero at the point of 1000 PgC of total CO2 emissions in the simulation - with the possibility for additional simulations, if resources allow. ZECMIP is part of CMIP6, under joint sponsorship by C4MIP and CDR-MIP, with associated experiment names to enable data submissions to the Earth System Grid Federation. All data will be published and made freely available.

Description: During R/V Meteor cruise 141/1, pore fluids of near surface sediments were investigated to find indications for hydrothermal activity in the Terceira Rift (TR), a hyper‐slow spreading center in the Central North Atlantic Ocean. To date, submarine hydrothermal fluid venting in the TR has only been reported for the D. João de Castro seamount, which presently seems to be inactive. Pore fluids sampled close to a volcanic cone at 2800 m water depth show an anomalous composition with Mg, SO4, and total alkalinity (TA) concentrations significantly higher than seawater and a nearby reference core. The most straightforward way of interpreting these deviations is the dissolution of the hydrothermally formed mineral caminite (MgSO4 0.25Mg(OH)2 0.2H2O). This interpretation is corroborated by a thorough investigation of fluid isotope systems (δ26Mg, δ30Si, δ34S, δ44/42Ca, and 87Sr/86Sr). Caminite is known from mineral assemblages with anhydrite, and forms in hydrothermal recharge zones only under specific conditions such as high fluid temperatures and in altered oceanic crust, which are conditions generally met at the TR. We hypothesize that caminite was formed during hydrothermal activity and is now dissolving during the waning state of the hydrothermal system, so that caminite mineralization is shifted out of its stability zone. Ongoing fluid circulation through the basement is transporting the geochemical signal via slow advection towards the seafloor.

Description: Mining of polymetallic nodules in abyssal seafloor sediments promises to address the growing worldwide demand for metallic minerals. Given that prospective mining operations are likely to have profound impacts on deep seafloor communities, industrial investment has been accompanied by scientific involvement for the assessment of baseline conditions and provision of guidelines for environmentally sustainable mining practices. Benthic meiofaunal communities were studied in four prospective mining areas of the Clarion–Clipperton Zone (CCZ) in the eastern Pacific Ocean, arranged in a southeast–northwest fashion coinciding with the productivity gradient in the area. Additionally, samples were collected from the Area of Particular Environmental Interest no. 3 (APEI-3) in the northwest of the CCZ, where mining will be prohibited and which should serve as a “source area” for the biota within the larger CCZ. Total densities in the 0–5 cm upper layer of the sediment were influenced by sedimentary characteristics, water depth and nodule density at the various sampling locations, indicating the importance of nodules for meiofaunal standing stock. Nematodes were the most abundant meiobenthic taxon, and their assemblages were typically dominated by a few genera (generally 2–6) accounting for 40 %–70 % of all individuals, which were also widely spread along the CCZ and shared among all sampled license areas. However, almost half of the communities consisted of rare genera, each contributing less than 5 % to the overall abundances and displaying a distribution which was usually restricted to a single license area. The same observations (dominant and widely spread versus rare and scattered) could be made for the species of one of the dominant genera, Halalaimus, implying that it might be mainly these rare genera and species that will be vulnerable to mining-induced changes in their habitat.

Description: Chloromethane (CH3Cl) is the most important natural input of reactive chlorine to the stratosphere, contributing about 16 % to stratospheric ozone depletion. Due to the phase-out of anthropogenic emissions of chlorofluorocarbons, CH3Cl will largely control future levels of stratospheric chlorine. The tropical rainforest is commonly assumed to be the strongest single CH3Cl source, contributing over half of the global annual emissions of about 4000 to 5000 Gg (1 Gg = 109 g). This source shows a characteristic carbon isotope fingerprint, making isotopic investigations a promising tool for improving its atmospheric budget. Applying carbon isotopes to better constrain the atmospheric budget of CH3Cl requires sound information on the kinetic isotope effects for the main sink processes: the reaction with OH and Cl in the troposphere. We conducted photochemical CH3Cl degradation experiments in a 3500 dm3 smog chamber to determine the carbon isotope effect (ε=k13C/k12C−1 ) for the reaction of CH3Cl with OH and Cl. For the reaction of CH3Cl with OH, we determined an ε value of (−11.2±0.8) ‰ (n=3) and for the reaction with Cl we found an ε value of (−10.2±0.5 ) ‰ (n=1), which is 5 to 6 times smaller than previously reported. Our smaller isotope effects are strongly supported by the lack of any significant seasonal covariation in previously reported tropospheric δ13C(CH3Cl) values with the OH-driven seasonal cycle in tropospheric mixing ratios. Applying these new values for the carbon isotope effect to the global CH3Cl budget using a simple two hemispheric box model, we derive a tropical rainforest CH3Cl source of (670±200) Gg a−1, which is considerably smaller than previous estimates. A revision of previous bottom-up estimates, using above-ground biomass instead of rainforest area, strongly supports this lower estimate. Finally, our results suggest a large unknown CH3Cl source of (1530±200) Gg a−1.

Description: Since 2009, unprecedented comprehensive long-term gas hydrate observations have become available from Ocean Networks Canada's NEPTUNE cabled ocean observatory at the northern Cascadia margin. Several experiments demonstrate the scientific importance of permanent power and Internet connectivity to the ocean floor as they have advanced the field of gas hydrate related research. One example is the cabled crawler Wally at Barkley Canyon, enabling live in situ exploration of the hydrate mounds and its associated benthic communities through the crawler's mobility and permanent accessibility throughout the year. Another example is a bubble-imaging sonar at Clayoquot Slope, revealing the strong relationship between ebullition of natural gas and tidal pressure, without apparent correlation to earthquakes, storms, or temperature fluctuations, in year-long continuous recordings. Finally, regular observatory maintenance cruises allow additional science sampling including echo-sounder surveys to extend the observatory footprint. Long-term trends in the data are not yet apparent but can also become evident from continuous measurements, as ocean observatories such as NEPTUNE are built for a 25-year lifetime, and expansion of the observatory networks makes these findings comparable and testable. Plain Language Summary Natural gas near the ocean floor creates a rapidly changing environment where it is important to collect data continuously in order to determine the magnitude, speed, and potential mechanism of change. This long-standing challenge of year-round access to the deep ocean has been tackled by Ocean Networks Canada through cabling the northern Cascadia seafloor, providing power and Internet communication-ideal for power-hungry instruments, large data volumes, and real-time access. The presence of gas influences the shape of the seafloor, animal activity, and potential escape of methane, a potent greenhouse gas. A seafloor crawler Wally was operated around deep canyon mounds of gas hydrate (a solid gas-water composite) since 2009 and helped discover environmental changes influencing sea life. Further along the continental slope, an acoustic sonar monitored rising methane bubbles where the bubbling appears to be controlled neither by earthquakes, winter storms, nor subtle temperature changes but actually strongly by tidal pressure. Regular maintenance of the observatory by ship allows more data to be collected near the cabled seafloor sites, extending the observations to a larger area. Ocean observatories are built to last decades and therefore more data for more research can be collected, potentially detecting relatively slow processes as well.

Description: Nitric oxide (NO) is a short-lived compound of the marine nitrogen cycle; however, our knowledge about its oceanic distribution and turnover is rudimentary. Here we present the measurements of dissolved NO in the surface and bottom layers at 75 stations in the Bohai Sea (BS) and the Yellow Sea (YS) in June 2011. Moreover, NO photoproduction rates were determined at 27 stations in both seas. The NO concentrations in the surface and bottom layers were highly variable and ranged from below the limit of detection (i.e., 32 pmol L−1) to 616 pmol L−1 in the surface layer and 482 pmol L−1 in the bottom layer. There was no significant difference (p〉0.05) between the mean NO concentrations in the surface (186±108 pmol L−1) and bottom (174±123 pmol L−1) layers. A decreasing trend of NO in bottom-layer concentrations with salinity indicates a NO input by submarine groundwater discharge. NO in the surface layer was supersaturated at all stations during both day and night and therefore the BS and YS were a persistent source of NO to the atmosphere at the time of our measurements. The average flux was about 4.5×10−16 mol cm−2 s−1 and the flux showed significant positive relationship with the wind speed. The accumulation of NO during daytime was a result of photochemical production, and photoproduction rates were correlated to illuminance. The persistent nighttime NO supersaturation pointed to an unidentified NO dark production. NO sea-to-air flux densities were much lower than the NO photoproduction rates. Therefore, we conclude that the bulk of the NO produced in the mixed layer was rapidly consumed before its release to the atmosphere.

Description: There is a converging body of evidence supporting a measurable slowdown of the Atlantic Meridional Overturning Circulation (AMOC) as climate warms and Northern Hemisphere ice sheets inexorably shrink. Within this context, we assess the variability of the AMOC during the Holocene based on a marine sediment core retrieved from the deep northwest Atlantic, which sensitively recorded large‐scale deglacial transitions in deep water circulation. While there is a diffuse notion of Holocene variability in Labrador and Nordic Seas overturning, we report a largely invariable deep water circulation for the last ~11,000 years, even during the meltwater pulse associated with the 8.2‐ka event. Sensitivity tests along with high‐resolution 231Pa/230Th data constrain the duration and the magnitude of possible Holocene AMOC variations. The generally constant baseline during the Holocene suggests attenuated natural variability of the large‐scale AMOC on submillennial timescales and calls for compensating effects involving the upstream components of North Atlantic Deep Water.

Description: The Arctic marine climate system is changing rapidly, which is seen in the warming of the ocean and atmosphere, decline of sea ice cover, increase in river discharge, acidification of the ocean, and changes in marine ecosystems. Socio-economic activities in the coastal and marine Arctic are simultaneously changing. This calls for the establishment of a marine Arctic component of the Pan-Eurasian Experiment (MA-PEEX). There is a need for more in situ observations on the marine atmosphere, sea ice, and ocean, but increasing the amount of such observations is a pronounced technological and logistical challenge. The SMEAR (Station for Measuring Ecosystem–Atmosphere Relations) concept can be applied in coastal and archipelago stations, but in the Arctic Ocean it will probably be more cost-effective to further develop a strongly distributed marine observation network based on autonomous buoys, moorings, autonomous underwater vehicles (AUVs), and unmanned aerial vehicles (UAVs). These have to be supported by research vessel and aircraft campaigns, as well as various coastal observations, including community-based ones. Major manned drifting stations may occasionally be comparable to terrestrial SMEAR flagship stations. To best utilize the observations, atmosphere–ocean reanalyses need to be further developed. To well integrate MA-PEEX with the existing terrestrial–atmospheric PEEX, focus is needed on the river discharge and associated fluxes, coastal processes, and atmospheric transports in and out of the marine Arctic. More observations and research are also needed on the specific socio-economic challenges and opportunities in the marine and coastal Arctic, and on their interaction with changes in the climate and environmental system. MA-PEEX will promote international collaboration; sustainable marine meteorological, sea ice, and oceanographic observations; advanced data management; and multidisciplinary research on the marine Arctic and its interaction with the Eurasian continent.

Description: Phytoplankton calcifiers contribute to global carbon cycling through their dual formation of calcium carbonate and particulate organic carbon (POC). The carbonate might provide an efficient export pathway for the associated POC to the deep ocean, reducing the particles' exposure to biological degradation in the upper ocean and increasing the particle settling rate. Previous work has suggested ballasting of POC by carbonate might increase in a warming climate, in spite of increasing carbonate dissolution rates, because calcifiers benefit from the widespread nutrient limitation arising from stratification. We compare the biogeochemical responses of three models containing (1) a single mixed phytoplankton class, (2) additional explicit small phytoplankton and calcifiers, and (3) additional explicit small phytoplankton and calcifiers with a prognostic carbonate ballast model, to two rapid changes in atmospheric CO2. The first CO2 scenario represents a rapid (151-year) transition from a stable icehouse climate (283.9 ppm) into a greenhouse climate (1263 ppm); the second represents a symmetric rapid transition from a stable greenhouse climate into an icehouse climate. We identify a slope change in the global net primary production response with a transition point at about 3.5 ∘C global mean sea surface temperature change in all models, driven by a combination of physical and biological changes. We also find that in both warming and cooling scenarios, the application of a prognostic carbonate ballast model moderates changes in carbon export production, suboxic volume, and nitrate sources and sinks, reducing the long-term model response to about one-third that of the calcifier model without ballast. Explicit small phytoplankton and calcifiers, and carbonate ballasting, increase the physical separation of nitrate sources and sinks through a combination of phytoplankton competition and lengthened remineralization profile, resulting in a significantly higher global nitrate inventory in this model compared to the single phytoplankton type model (15 % and 32 % higher for icehouse and greenhouse climates). Higher nitrate inventory alleviates nitrate limitation, increasing phytoplankton sensitivity to changes in physical limitation factors (light and temperature). This larger sensitivity to physical forcing produces stronger shifts in ocean phosphate storage between icehouse and greenhouse climates. The greenhouse climate is found to hold phosphate and nitrate deeper in the ocean, despite a shorter remineralization length scale than the icehouse climate, because of the longer residence times of the deep water masses. We conclude the global biogeochemical impact of calcifiers extends beyond their role in global carbon cycling, and that the ecological composition of the global ocean can affect how ocean biogeochemistry responds to climate forcing.

Description: Natural products and derivatives thereof are of considerable importance in the discovery of new pharmaceuticals, for example, for the treatment of cancer, diabetes, inflammation diseases, and infection diseases caused by bacteria, fungi, viruses, or parasites. The great biodiversity of marine microorganisms is reflected in their huge chemical diversity, which provides a rich source of biologically active compounds. An increasing interest in marine microorganisms as promising producers of new compounds with potential medical applications has raised increasing interest in the sustainable exploration of marine microbial resources for the discovery of new antibiotics, which is highlighted. The bottlenecks in the development of drugs using the large marine natural product pipeline are also discussed.

Description: The coastal upwelling regime off Peru in December 2012 showed considerable vertical concentration gradients of dissolved nitrous oxide (N2O) across the top few meters of the ocean. The gradients were predominantly downward, i.e., concentrations decreased toward the surface. Ignoring these gradients causes a systematic error in regionally integrated gas exchange estimates, when using observed concentrations at several meters below the surface as input for bulk flux parameterizations – as is routinely practiced. Here we propose that multi-day near-surface stratification events are responsible for the observed near-surface N2O gradients, and that the gradients induce the strongest bias in gas exchange estimates at winds of about 3 to 6 m s−1. Glider hydrographic time series reveal that events of multi-day near-surface stratification are a common feature in the study region. In the same way as shorter events of near-surface stratification (e.g., the diurnal warm layer cycle), they preferentially exist under calm to moderate wind conditions, suppress turbulent mixing, and thus lead to isolation of the top layer from the waters below (surface trapping). Our observational data in combination with a simple gas-transfer model of the surface trapping mechanism show that multi-day near-surface stratification can produce near-surface N2O gradients comparable to observations. They further indicate that N2O gradients created by diurnal or shorter stratification cycles are weaker and do not substantially impact bulk emission estimates. Quantitatively, we estimate that the integrated bias for the entire Peruvian upwelling region in December 2012 represents an overestimation of the total N2O emission by about a third, if concentrations at 5 or 10 m depth are used as surrogate for bulk water N2O concentration. Locally, gradients exist which would lead to emission rates overestimated by a factor of two or more. As the Peruvian upwelling region is an N2O source of global importance, and other strong N2O source regions could tend to develop multi-day near-surface stratification as well, the bias resulting from multi-day near-surface stratification may also impact global oceanic N2O emission estimates.

Description: The trait‐based approach is increasingly used in plankton ecology to understand diversity, community dynamics, and biogeography. While on the global scale phytoplankton traits are fairly well established, zooplankton traits are only beginning to be understood. One taxa‐transcending aspect of zooplankton diversity is the distinction between ambush and active feeding strategies. We present a global‐scale empirical estimate of feeding strategy derived from copepod abundance observations, which for the first time suggests a distinct trait biogeography with ambush feeding as the dominant feeding strategy at higher, but not at lower latitudes. To explain this trait biogeography, we develop a minimalist trade‐off based model of feeding strategies based on encounter rates between zooplankton predators and their phyto‐ and zooplankton prey. Encounter rates are governed by the two traits, size and motility, that trade off against predation risk. Coupled to a three‐dimensional dynamic green ocean model, our idealized encounter model captures the observed feeding strategy biogeography. In the model, this pattern arises from competing dominant food chains within the food web and is shaped by a trophic trait cascade of active vs. passive feeding in adjacent trophic levels. The dominant feeding strategy structures the pathways and efficiency of energy and biomass transfer through the model food web, with consequences for primary production, export and higher trophic levels. Understanding feeding strategies is therefore important for fisheries, biogeochemical cycling, and long‐term predictions of ecosystem dynamics and functioning by global dynamic green ocean models.

Description: As one of the most abundant and ubiquitous representatives of marine and brackish coastal macrophytobenthos communities, the genus Ulva is not only an important primary producer but also of ecological and morphogenetic interest to many scientists. Ulva mutabilis became an important model organism to study morphogenesis and mutualistic interactions of macroalgae and microorganisms. Here, we report that our collections of Ulva compressa Linnaeus (1753) from Germany are conspecific with the type strains of the model organism U. mutabilis Føyn (1958), which were originally collected at Olhão on the south coast of Portugal and have from that time on been maintained in culture as gametophytic and parthenogenetic lab strains. Different approaches were used to test conspecificity: (i) comparisons of vegetative and reproductive features of cultured material of U. mutabilis and German U. compressa demonstrated a shared morphological pattern; (ii) gametes of U. compressa and U. mutabilis successfully mated and developed into fertile sporophytic first‐generation offspring; (iii) molecular phylogenetics and species delimitation analyses based on the Generalized Mixed Yule‐Coalescent method showed that U. mutabilis isolates (sl‐G[mt+]) and (wt‐G[mt‐]) and U. compressa belong to a unique Molecular Operational Taxonomic Unit. According to these findings, there is sufficient evidence that U. mutabilis and U. compressa should be regarded as conspecific.

Description: The coastal waters of the Baltic Sea are subject to high variations in environmental conditions, triggered by natural and anthropogenic causes. Thus, in situ measurements of water parameters can be strategic for our understanding of the dynamics in shallow water habitats. In this study we present the results of a monitoring program at low water depths (1–2.5 m), covering 13 stations along the Baltic coast of Schleswig-Holstein, Germany. The provided dataset consists of records for dissolved inorganic nutrient concentrations taken twice a month and continuous readings at 10 min intervals for temperature, salinity and oxygen content. Data underwent quality control procedures and were flagged. On average, a data availability of 〉90 % was reached for the monitoring period within 2016–2018. The obtained monitoring data reveal great temporal and spatial variabilities of key environmental factors for shallow water habitats in the southwestern Baltic Sea. Therefore the presented information could serve as realistic key data for experimental manipulations of environmental parameters as well as for the development of oceanographic, biogeochemical or ecological models.

Description: Five initialization and ensemble generation methods are investigated with respect to their impact on the prediction skill of the German decadal prediction system "Mittelfristige Klimaprognose" (MiKlip). Among the tested methods, three tackle aspects of model‐consistent initialization using the ensemble Kalman filter (EnKF), the filtered anomaly initialization (FAI) and the initialization method by partially coupled spin‐up (MODINI). The remaining two methods alter the ensemble generation: the ensemble dispersion filter (EDF) corrects each ensemble member with the ensemble mean during model integration. And the bred vectors (BV) perturb the climate state using the fastest growing modes. The new methods are compared against the latest MiKlip system in the low‐resolution configuration (Preop‐LR), which uses lagging the climate state by a few days for ensemble generation and nudging toward ocean and atmosphere reanalyses for initialization. Results show that the tested methods provide an added value for the prediction skill as compared to Preop‐LR in that they improve prediction skill over the eastern and central Pacific and different regions in the North Atlantic Ocean. In this respect, the EnKF and FAI show the most distinct improvements over Preop‐LR for surface temperatures and upper ocean heat content, followed by the BV, the EDF and MODINI. However, no single method exists that is superior to the others with respect to all metrics considered. In particular, all methods affect the Atlantic Meridional Overturning Circulation in different ways, both with respect to the basin‐wide long‐term mean and variability, and with respect to the temporal evolution at the 26° N latitude.

Description: Ground-breaking measurements from the ocean observatories initiative Irminger Sea surface mooring (60°N, 39°30′W) are presented that provide the first in situ characterization of multiwinter surface heat exchange at a high latitude North Atlantic site. They reveal strong variability (December 2014 net heat loss nearly 50% greater than December 2015) due primarily to variations in frequency of intense short timescale (1–3 days) forcing. Combining the observations with the new high resolution European Centre for Medium Range Weather Forecasts Reanalysis 5 (ERA5) atmospheric reanalysis, the main source of multiwinter variability is shown to be changes in the frequency of Greenland tip jets (present on 15 days in December 2014 and 3 days in December 2015) that can result in hourly mean heat loss exceeding 800 W/m2. Furthermore, a new picture for atmospheric mode influence on Irminger Sea heat loss is developed whereby strongly positive North Atlantic Oscillation conditions favor increased losses only when not outweighed by the East Atlantic Pattern.

Description: In order to understand the effect of global change on marine fishes, it is imperative to quantify the effects on fundamental parameters such as survival and growth. Larval survival and recruitment of the Atlantic cod (Gadus morhua) were found to be heavily impaired by end-of-century levels of ocean acidification. Here, we analysed larval growth among 35–36 days old surviving larvae, along with organ development and ossification of the skeleton. We combined CO2treatments (ambient: 503 µatm, elevated: 1,179 µatm) with food availability in order to evaluate the effect of energy limitation in addition to the ocean acidification stressor. As expected, larval size (as a proxy for growth) and skeletogenesis were positively affected by high food availability. We found significant interactions between acidification and food availability. Larvae fed ad libitum showed little difference in growth and skeletogenesis due to the CO2 treatment. Larvae under energy limitation were significantly larger and had further developed skeletal structures in the elevated CO2 treatment compared to the ambient CO2 treatment. However, the elevated CO2 group revealed impairments in critically important organs, such as the liver, and had comparatively smaller functional gills indicating a mismatch between size and function. It is therefore likely that individual larvae that had survived acidification treatments will suffer from impairments later during ontogeny. Our study highlights important allocation trade-off between growth and organ development, which is critically important to interpret acidification effects on early life stages of fish.

Description: Eddy covariance measurements show gas transfer velocity suppression at medium to high wind speed. A wind-wave interaction described by the transformed Reynolds number is used to characterize environmental conditions favoring this suppression. We take the transformed Reynolds number parameterization to review the two most cited wind speed gas transfer velocity parameterizations: Nightingale et al. (2000) and Wanninkhof (1992, 2014). We propose an algorithm to adjust k values for the effect of gas transfer suppression and validate it with two directly measured dimethyl sulfide (DMS) gas transfer velocity data sets that experienced gas transfer suppression. We also show that the data set used in the Nightingale 2000 parameterization experienced gas transfer suppression. A compensation of the suppression effect leads to an average increase of 22% in the k vs. u relationship. Performing the same correction for Wanninkhof 2014 leads to an increase of 9.85 %. Additionally, we applied our gas transfer suppression algorithm to global air-sea flux climatologies of CO2 and DMS. The global application of gas transfer suppression leads to a decrease of 11% in DMS outgassing. We expect the magnitude of Reynolds suppression on any global air-sea gas exchange to be about 10%

Description: Industrial interest in deep-sea mineral extraction began decades ago and today it is at an all-time high, accelerated by global demand for metals. Several seafloor ecosystem disturbance experiments were performed beginning in the 1970’s, including the DISturbance and reCOLonization experiment (DISCOL) conducted in the Peru Basin in 1989. A large seafloor disturbance was created by repeatedly plowing the seafloor over an area of ~ 10.8 km2. Though a number of studies in abyssal mining regions have evaluated megafaunal biodiversity and ecosystem responses, few have included quantitative and detailed data on fishes or scavengers despite their ecological importance as top predators. We used towed camera transects and baited camera data to evaluate the fish community at the DISCOL site. The abyssal fish community was relatively diverse with 16 taxa dominated by Ipnops meadi. Fish density was lower in ploughed habitat during the several years following disturbance but thereafter increased over time in part due to changes in regional environmental conditions. 26 years post disturbance there were no differences in overall total fish densities between reference and experimental areas, but the dominant fish, I. meadi, still exhibited much lower densities in ploughed habitat suggesting only partial fish community recovery. The scavenging community was dominated by eelpouts (Pachycara spp), hermit crabs (Probeebei mirabilis) and shrimp. The large contribution of hermit crabs appears unique amongst abyssal scavenger studies worldwide. The abyssal fish community at DISCOL was similar to that in the more northerly Clarion Clipperton Zone, though some species have only been observed at DISCOL thus far. Also, further species level identifications are required to refine this assessment. Additional studies across the polymetallic nodule provinces of the Pacific are required to further evaluate the environmental drivers of fish density and diversity and species biogeographies, which will be important for the development of appropriate management plans aimed at minimizing human impact from deep-sea mining.

Description: Genetic divergence among populations arises through natural selection or drift and is counteracted by connectivity and gene flow. In sympatric populations, isolating mechanisms are thus needed to limit the homogenizing effects of gene flow to allow for adaptation and speciation. Chromosomal inversions act as an important mechanism maintaining isolating barriers, yet their role in sympatric populations and divergence with gene flow is not entirely understood. Here, we revisit the question of whether inversions play a role in the divergence of connected populations of the marine fish Atlantic cod (Gadus morhua), by exploring a unique data set combining whole-genome sequencing data and behavioural data obtained with acoustic telemetry. Within a confined fjord environment, we find three genetically differentiated Atlantic cod types belonging to the oceanic North Sea population, the western Baltic population and a local fjord-type cod. Continuous behavioural tracking over 4 year revealed temporally stable sympatry of these types within the fjord. Despite overall weak genetic differentiation consistent with high levels of gene flow, we detected significant frequency shifts of three previously identified inversions, indicating an adaptive barrier to gene flow. In addition, behavioural data indicated that North Sea cod and individuals homozygous for the LG12 inversion had lower fitness in the fjord environment. However, North Sea and fjord-type cod also occupy different depths, possibly contributing to prezygotic reproductive isolation and representing a behavioural barrier to gene flow. Our results provide the first insights into a complex interplay of genomic and behavioural isolating barriers in Atlantic cod and establish a new model system towards an understanding of the role of genomic structural variants in adaptation and diversification.

Description: The Bjerknes feedback is the dominant positive feedback in the equatorial ocean basins. To examine the seasonality, symmetry, and stationarity of the Pacific and Atlantic Bjerknes feedbacks we decompose them into three feedback elements that relate thermocline depth, sea surface temperature (SST), and western basin wind stress variability to each other. We partition feedback elements into composites associated with positive or negative anomalies. Using robust regression, we diagnose the strength of each composite. For the recent period 1993‐2012, composites of the Pacific Bjerknes feedback elements agree well with previous work. Positive composites are generally stronger than negative composites, and all feedback elements are weakest in late boreal spring. In the Atlantic, differences between positive and negative composites are less consistent across feedback elements. Specifically, wind variability seems to play a less important role in shaping atmosphere‐ocean coupling in the Atlantic when compared to the Pacific. However, a clear seasonality emerges: Feedback elements are generally strong in boreal summer and, for the negative composites, again in boreal winter. The Atlantic Bjerknes feedback is dominated by subsurface‐surface coupling. Applying our analysis to overlapping 25‐year periods for 1958‐2009 shows that the strengths of feedback elements in both ocean basins vary on decadal time scales. While the overall asymmetry of the Pacific Bjerknes feedback is robust, the strength and symmetry of Atlantic feedback elements vary considerably between decades. Our results indicate that the Atlantic Bjerknes feedback is non‐stationary on decadal time scales.

Description: Due to strong mean state‐biases most coupled models are unable to simulate equatorial Atlantic variability. Here, we use the Kiel Climate Model to assess the impact of bias reduction on the seasonal prediction of equatorial Atlantic sea surface temperature (SST). We compare a standard experiment (STD) with an experiment that employs surface heat flux correction to reduce the SST bias (FLX) and, in addition, apply a correction for initial errors in SST. Initial conditions for both experiments are generated in partially coupled mode, and seasonal hindcasts are initialized at the beginning of February, May, August and November for 1981–2012. Surface heat flux correction generally improves hindcast skill. Hindcasts initialized in February have the least skill, even though the model bias is not particularly strong at that time of year. In contrast, hindcasts initialized in May achieve the highest skill. We argue this is because of the emergence of a closed Bjerknes feedback loop in boreal summer in FLX that is a feature of observations but is missing in STD.