ALBERT

Description: Tropical South America is one of the three main centres of the global, zonal overturning circulation of the equatorial atmosphere (generally termed the 'Walker' circulation1). Although this area plays a key role in global climate cycles, little is known about South American climate history. Here we describe sediment cores and down-hole logging results of deep drilling in the Salar de Uyuni, on the Bolivian Altiplano, located in the tropical Andes. We demonstrate that during the past 50,000 years the Altiplano underwent important changes in effective moisture at both orbital (20,000-year) and millennial timescales. Long-duration wet periods, such as the Last Glacial Maximum—marked in the drill core by continuous deposition of lacustrine sediments—appear to have occurred in phase with summer insolation maxima produced by the Earth's precessional cycle. Short-duration, millennial events correlate well with North Atlantic cold events, including Heinrich events 1 and 2, as well as the Younger Dryas episode. At both millennial and orbital timescales, cold sea surface temperatures in the high-latitude North Atlantic were coeval with wet conditions in tropical South America, suggesting a common forcing.

Description: Marine phytoplankton blooms are annual spring events that sustain active and diverse bloom-associated bacterial populations. Blooms vary considerably in terms of eukaryotic species composition and environmental conditions, but a limited number of heterotrophic bacterial lineages — primarily members of the Flavobacteriia, Alphaproteobacteria and Gammaproteobacteria — dominate these communities. In this Review, we discuss the central role that these bacteria have in transforming phytoplankton-derived organic matter and thus in biogeochemical nutrient cycling. On the basis of selected field and laboratory-based studies of flavobacteria and roseobacters, distinct metabolic strategies are emerging for these archetypal phytoplankton-associated taxa, which provide insights into the underlying mechanisms that dictate their behaviours during blooms.

Description: Owing to the turbulent nature of the ocean, mesoscale eddies are omnipresent. The impact of these transitory and approximately circular sea surface temperature fronts on the overlying atmosphere is not well known. Stationary fronts such as the Gulf Stream have been reported to lead to pronounced atmospheric changes1, 2. However, the impact of transient ocean eddies on the atmosphere has not been determined systematically, except on winds and to some extent clouds3, 4, 5, 6. Here, we examine the atmospheric conditions associated with over 600,000 individual eddies in the Southern Ocean, using satellite data. We show that ocean eddies locally affect near-surface wind, cloud properties and rainfall. The observed pattern of atmospheric change is consistent with a mechanism in which sea surface temperature anomalies associated with the oceanic eddies modify turbulence in the atmospheric boundary layer. In the case of cyclonic eddies, this modification triggers a slackening of near-surface winds, a decline in cloud fraction and water content, and a reduction in rainfall. We conclude that transient mesoscale ocean structures can significantly affect much larger atmospheric low-pressure systems that swiftly pass by at the latitudes investigated.

Description: An influence of solar irradiance variations on Earth’s surface climate has been repeatedly suggested, based on correlations between solar variability and meteorological variables1. Specifically, weaker westerly winds have been observed in winters with a less active sun, for example at the minimum phase of the 11-year sunspot cycle2, 3, 4. With some possible exceptions5, 6, it has proved difficult for climate models to consistently reproduce this signal7, 8. Spectral Irradiance Monitor satellite measurements indicate that variations in solar ultraviolet irradiance may be larger than previously thought9. Here we drive an ocean–atmosphere climate model with ultraviolet irradiance variations based on these observations. We find that the model responds to the solar minimum with patterns in surface pressure and temperature that resemble the negative phase of the North Atlantic or Arctic Oscillation, of similar magnitude to observations. In our model, the anomalies descend through the depth of the extratropical winter atmosphere. If the updated measurements of solar ultraviolet irradiance are correct, low solar activity, as observed during recent years, drives cold winters in northern Europe and the United States, and mild winters over southern Europe and Canada, with little direct change in globally averaged temperature. Given the quasiregularity of the 11-year solar cycle, our findings may help improve decadal climate predictions for highly populated extratropical regions

Description: Cold-water coral ecosystems are considered hot-spots of biodiversity and biomass production and may be a regionally important contributor to carbonate production. The impact of these ecosystems on biogeochemical processes and carbonate preservation in associated sediments were studied at Røst Reef and Traenadjupet Reef, two modern (post-glacial) cold-water coral reefs on the Mid-Norwegian shelf. Sulfate and iron reduction as well as carbonate dissolution and precipitation were investigated by combining pore-water geochemical profiles, steady state modeling, as well as solid phase analyses and sulfate reduction rate measurements on gravity cores of up to 3.25 m length. Low extents of sulfate depletion and dissolved inorganic carbon (DIC) production, combined with sulfate reduction rates not exceeding 3 nmol S cm−3 d−1, suggested that overall anaerobic carbon mineralization in the sediments was low. These data showed that the coral fragment-bearing siliciclastic sediments were effectively decoupled from the productive pelagic ecosystem by the complex reef surface framework. Organic matter being mineralized by sulfate reduction was calculated to consist of 57% carbon bound in CH2O groups and 43% carbon in -CH2- groups. Methane concentrations were below 1 μM, and failed to support the hypothesis of a linkage between the distribution of cold-water coral reefs and the presence of hydrocarbon seepage. Reductive iron oxide dissolution linked to microbial sulfate reduction buffered the pore-water carbonate system and inhibited acid-driven coral skeleton dissolution. A large pool of reactive iron was available leading to the formation of iron sulfide minerals. Constant pore-water Ca2+, Mg2+ and Sr2+ concentrations in most cores and decreasing Ca2+ and Sr2+ concentrations with depth in core 23–18 GC indicated diagenetic carbonate precipitation. This was consistent with the excellent preservation of buried coral fragments.

Description: The anaerobic oxidation of methane (AOM) with sulphate, an area currently generating great interest in microbiology, is accomplished by consortia of methanotrophic archaea (ANME) and sulphate-reducing bacteria1, 2. The enzyme activating methane in methanotrophic archaea has tentatively been identified as a homologue of methyl-coenzyme M reductase (MCR) that catalyses the methane-forming step in methanogenic archaea3, 4. Here we report an X-ray structure of the 280 kDa heterohexameric ANME-1 MCR complex. It was crystallized uniquely from a protein ensemble purified from consortia of microorganisms collected with a submersible from a Black Sea mat catalysing AOM with sulphate4. Crystals grown from the heterogeneous sample diffract to 2.1 Å resolution and consist of a single ANME-1 MCR population, demonstrating the strong selective power of crystallization. The structure revealed ANME-1 MCR in complex with coenzyme M and coenzyme B, indicating the same substrates for MCR from methanotrophic and methanogenic archaea. Differences between the highly similar structures of ANME-1 MCR and methanogenic MCR include a F430 modification, a cysteine-rich patch and an altered post-translational amino acid modification pattern, which may tune the enzymes for their functions in different biological contexts.

Description: According to small subunit ribosomal RNA (ss rRNA) sequence comparisons all known Archaea belong to the phyla Crenarchaeota, Euryarchaeota, and—indicated only by environmental DNA sequences—to the 'Korarchaeota'1, 2. Here we report the cultivation of a new nanosized hyperthermophilic archaeon from a submarine hot vent. This archaeon cannot be attached to one of these groups and therefore must represent an unknown phylum which we name 'Nanoarchaeota' and species, which we name 'Nanoarchaeum equitans'. Cells of 'N. equitans' are spherical, and only about 400 nm in diameter. They grow attached to the surface of a specific archaeal host, a new member of the genus Ignicoccus3. The distribution of the 'Nanoarchaeota' is so far unknown. Owing to their unusual ss rRNA sequence, members remained undetectable by commonly used ecological studies based on the polymerase chain reaction4. 'N. equitans' harbours the smallest archaeal genome; it is only 0.5 megabases in size. This organism will provide insight into the evolution of thermophily, of tiny genomes and of interspecies communication.

Description: In order to test the influences of ocean acidification on the ocean pelagic ecosystem, so far the largest CO2 manipulation mesocosm study (European Project on Ocean Acidification, EPOCA) was performed in Kings Bay (Kongsfjorden), Spitsbergen. During a 30 day incubation, bacterial diversity was investigated using DNA fingerprinting and clone library analysis of bacterioplankton samples. Terminal restriction fragment length polymorphism (T-RFLP) analysis of the PCR amplicons of the 16S rRNA genes revealed that general bacterial diversity, taxonomic richness and community structure were influenced by the variation of productivity during the time of incubation, but not the degree of ocean acidification. A BIOENV analysis suggested a complex control of bacterial community structure by various biological and chemical environmental parameters. The maximum apparent diversity of bacterioplankton (i.e., the number of T-RFs) in high and low pCO2 treatments differed significantly. A negative relationship between the relative abundance of Bacteroidetes and pCO2 levels was observed for samples at the end of the experiment by the combination of T-RFLP and clone library analysis. Our study suggests that ocean acidification affects the development of bacterial assemblages and potentially impacts the ecological function of the bacterioplankton in the marine ecosystem.

Description: The Atlantic meridional overturning circulation (AMOC) carries large amounts of heat into the North Atlantic influencing climate regionally as well as globally. Palaeo-records and simulations with comprehensive climate models suggest that the positive salt-advection feedback may yield a threshold behaviour of the system. That is to say that beyond a certain amount of freshwater flux into the North Atlantic, no meridional overturning circulation can be sustained. Concepts of monitoring the AMOC and identifying its vicinity to the threshold rely on the fact that the volume flux defining the AMOC will be reduced when approaching the threshold. Here we advance conceptual models that have been used in a paradigmatic way to understand the AMOC, by introducing a density-dependent parameterization for the Southern Ocean eddies. This additional degree of freedom uncovers a mechanism by which the AMOC can increase with additional freshwater flux into the North Atlantic, before it reaches the threshold and collapses: an AMOC that is mainly wind-driven will have a constant upwelling as long as the Southern Ocean winds do not change significantly. The downward transport of tracers occurs either in the northern sinking regions or through Southern Ocean eddies. If freshwater is transported, either atmospherically or via horizontal gyres, from the low to high latitudes, this would reduce the eddy transport and by continuity increase the northern sinking which defines the AMOC until a threshold is reached at which the AMOC cannot be sustained. If dominant in the real ocean this mechanism would have significant consequences for monitoring the AMOC.

Description: We present bathymetry and hydrological observations collected in the summer of 2014 from two fjord systems in southeastern Greenland with a multibeam sonar system. Our results provide a detailed bathymetric map of the fjord complex around the island of Skjoldungen in Skjoldungen Fjord and the outer part of Timmiarmiut Fjord and show far greater depths compared to the International Bathymetric Chart of the Arctic Ocean. The hydrography collected shows different properties in the fjords with the bottom water masses below 240 m in Timmiarmiut Fjord being 1–2 °C warmer than in the two fjords around Skjoldungen, but data also illustrate the influence of sills on the exchange of deeper water masses within fjords. Moreover, evidence of subglacial discharge in Timmiarmiut Fjord, which is consistent with satellite observations of ice mélange set into motion, adds to our increasing understanding of the distribution of subglacial meltwater.

Description: Iron limits phytoplankton growth and hence the biological carbon pump in the Southern Ocean1. Models assessing the impacts of iron on the global carbon cycle generally rely on dust input and sediment resuspension as the predominant sources2, 3. Although it was previously thought that most iron from deep-ocean hydrothermal activity was inaccessible to phytoplankton because of the formation of particulates4, it has been suggested that iron from hydrothermal activity5, 6, 7 may be an important source of oceanic dissolved iron8, 9, 10, 11, 12, 13. Here we use a global ocean model to assess the impacts of an annual dissolved iron flux of approximately 9×108 mol, as estimated from regional observations of hydrothermal activity11, 12, on the dissolved iron inventory of the world’s oceans. We find the response to the input of hydrothermal dissolved iron is greatest in the Southern Hemisphere oceans. In particular, observations of the distribution of dissolved iron in the Southern Ocean3 (Chever et al., manuscript in preparation; Bowie et al., manuscript in preparation) can be replicated in our simulations only when our estimated iron flux from hydrothermal sources is included. As the hydrothermal flux of iron is relatively constant over millennial timescales14, we propose that hydrothermal activity can buffer the oceanic dissolved iron inventory against shorter-term fluctuations in dust deposition.

Description: Ocean acidification resulting from the uptake of anthropogenic carbon dioxide (CO2) by the ocean is considered a major threat to marine ecosystems. Here we examined the effects of ocean acidification on microbial community dynamics in the eastern Baltic Sea during the summer of 2012 when inorganic nitrogen and phosphorus were strongly depleted. Large-volume in situ mesocosms were employed to mimic present, future and far future CO2 scenarios. All six groups of phytoplankton enumerated by flow cytometry ( 〈  20 µm cell diameter) showed distinct trends in net growth and abundance with CO2 enrichment. The picoeukaryotic phytoplankton groups Pico-I and Pico-II displayed enhanced abundances, whilst Pico-III, Synechococcus and the nanoeukaryotic phytoplankton groups were negatively affected by elevated fugacity of CO2 (fCO2). Specifically, the numerically dominant eukaryote, Pico-I, demonstrated increases in gross growth rate with increasing fCO2 sufficient to double its abundance. The dynamics of the prokaryote community closely followed trends in total algal biomass despite differential effects of fCO2 on algal groups. Similarly, viral abundances corresponded to prokaryotic host population dynamics. Viral lysis and grazing were both important in controlling microbial abundances. Overall our results point to a shift, with increasing fCO2, towards a more regenerative system with production dominated by small picoeukaryotic phytoplankton.

Description: Pliensbachian and Toarcian Ostracoda first described by Exton (Geological Paper, Carleton University, Ottawa, 79: 1–104 1979) from the Lusitanian Basin, west-central Portugal have been re-examined. As a result, a greater diversity in the Ostracoda (80 species) is now recognized. Two species are newly described (Eucytherura zambujalensis sp. nov., Ektyphocythere mediodepressa sp. nov.) from the marls and calcareous shales of the Maria Pares Hill section near the village of Zambujal. Poor preservation precludes a complete taxonomic review of the present material. Five ostracod zones are proposed; Gammacythere ubiquita–Ogmoconchella gruendeli Zone, Poly cope cerasia–Polycope cincinnata Zone, Liasina lanceolata–Ogmoconcha convexa Zone, Bairdiacypris rectangularis–Kinkelinella sermoisensis Zone, and Cytherella toarcensis-Kinkelinella costata Zone. Although the ostracod assemblages possess strong similarities to those described from Northwest Europe, some of the Zambujal assemblages are dominated by the genus Polycope. A marked faunal turnover, in association with the extinction of the Metacopina occurs in the lower Subzone of the tenuicostatum Zone of Lower Toarcian age. These faunal events are discussed in relation to changing environmental conditions.

Description: Biogeochemical models, capturing the major feedbacks of the pelagic ecosystem of the world ocean, are today often embedded into Earth system models which are increasingly used for decision making regarding climate policies. These models contain poorly constrained parameters (e.g., maximum phytoplankton growth rate), which are typically adjusted until the model shows reasonable behavior. Systematic approaches determine these parameters by minimizing the misfit between the model and observational data. In most common model approaches, however, the underlying functions mimicking the biogeochemical processes are nonlinear and non-convex. Thus, systematic optimization algorithms are likely to get trapped in local minima and might lead to non-optimal results. To judge the quality of an obtained parameter estimate, we propose determining a preferably large lower bound for the global optimum that is relatively easy to obtain and that will help to assess the quality of an optimum, generated by an optimization algorithm. Due to the unavoidable noise component in all observations, such a lower bound is typically larger than zero. We suggest deriving such lower bounds based on typical properties of biogeochemical models (e.g., a limited number of extremes and a bounded time derivative). We illustrate the applicability of the method with two real-world examples. The first example uses real-world observations of the Baltic Sea in a box model setup. The second example considers a three-dimensional coupled ocean circulation model in combination with satellite chlorophyll a.

Description: The late stage of the North East Atlantic (NEA) spring bloom was investigated during June 2005 along a transect section from 45 to 66° N between 15 and 20° W in order to characterize the contribution of siliceous and calcareous phytoplankton groups and describe their distribution in relation to environmental factors. We measured several biogeochemical parameters such as nutrients, surface trace metals, algal pigments, biogenic silica (BSi), particulate inorganic carbon (PIC) or calcium carbonate, particulate organic carbon, nitrogen and phosphorus (POC, PON and POP, respectively), as well as transparent exopolymer particles (TEP). Results were compared with other studies undertaken in this area since the JGOFS NABE program. Characteristics of the spring bloom generally agreed well with the accepted scenario for the development of the autotrophic community. The NEA seasonal diatom bloom was in the late stages when we sampled the area and diatoms were constrained to the northern part of our transect, over the Icelandic Basin (IB) and Icelandic Shelf (IS). Coccolithophores dominated the phytoplankton community, with a large distribution over the Rockall-Hatton Plateau (RHP) and IB. The Porcupine Abyssal Plain (PAP) region at the southern end of our transect was the region with the lowest biomass, as demonstrated by very low Chla concentrations and a community dominated by picophytoplankton. Early depletion of dissolved silicic acid (DSi) and increased stratification of the surface layer most likely triggered the end of the diatom bloom, leading to coccolithophore dominance. The chronic Si deficiency observed in the NEA could be linked to moderate Fe limitation, which increases the efficiency of the Si pump. TEP closely mirrored the distribution of both biogenic silica at depth and prymnesiophytes in the surface layer suggesting the sedimentation of the diatom bloom in the form of aggregates, but the relative contribution of diatoms and coccolithophores to carbon export in this area still needs to be resolved.

Description: Very short-lived substances (VSLS) contribute as source gases significantly to the tropospheric and stratospheric bromine loading. At present, an estimated 25% of stratospheric bromine is of oceanic origin. In this study, we investigate how climate change may impact the ocean- atmosphere flux of brominated VSLS, their atmospheric transport, and chemical transformations and evaluate how these changes will affect stratospheric ozone over the 21st century. Under the assumption of fixed ocean water concentrations and RCP6.0 scenario, we find an increase of the ocean- atmosphere flux of brominated VSLS of about 8-10% by the end of the 21st century compared to present day. A decrease in the tropospheric mixing ratios of VSLS and an increase in the lower stratosphere are attributed to changes in atmospheric chemistry and transport. Our model simulations reveal that this increase is counteracted by a corresponding reduction of inorganic bromine. Therefore the total amount of bromine from VSLS in the stratosphere will not be changed by an increase in upwelling. Part of the increase of VSLS in the tropical lower stratosphere results from an increase in the corresponding tropopause height. As the depletion of stratospheric ozone due to bromine depends also on the availability of chlorine, we find the impact of bromine on stratospheric ozone at the end of the 21st century reduced compared to present day. Thus, these studies highlight the different factors influencing the role of brominated VSLS in a future climate

Description: Reanalysis data sets are widely used to understand atmospheric processes and past variability, and are often used to stand in as "observations" for comparisons with climate model output. Because of the central role of water vapor (WV) and ozone (O3) in climate change, it is important to understand how accurately and consistently these species are represented in existing global reanalyses. In this paper, we present the results of WV and O3 intercomparisons that have been performed as part of the SPARC (Stratosphere–troposphere Processes and their Role in Climate) Reanalysis Intercomparison Project (S-RIP). The comparisons cover a range of timescales and evaluate both inter-reanalysis and observation-reanalysis differences. We also provide a systematic documentation of the treatment of WV and O3 in current reanalyses to aid future research and guide the interpretation of differences amongst reanalysis fields. The assimilation of total column ozone (TCO) observations in newer reanalyses results in realistic representations of TCO in reanalyses except when data coverage is lacking, such as during polar night. The vertical distribution of ozone is also relatively well represented in the stratosphere in reanalyses, particularly given the relatively weak constraints on ozone vertical structure provided by most assimilated observations and the simplistic representations of ozone photochemical processes in most of the reanalysis forecast models. However, significant biases in the vertical distribution of ozone are found in the upper troposphere and lower stratosphere in all reanalyses. In contrast to O3, reanalysis estimates of stratospheric WV are not directly constrained by assimilated data. Observations of atmospheric humidity are typically used only in the troposphere, below a specified vertical level at or near the tropopause. The fidelity of reanalysis stratospheric WV products is therefore mainly dependent on the reanalyses' representation of the physical drivers that influence stratospheric WV, such as temperatures in the tropical tropopause layer, methane oxidation, and the stratospheric overturning circulation. The lack of assimilated observations and known deficiencies in the representation of stratospheric transport in reanalyses result in much poorer agreement amongst observational and reanalysis estimates of stratospheric WV. Hence, stratospheric WV products from the current generation of reanalyses should generally not be used in scientific studies.

Description: The Sumatran subduction zone exhibits strong seismic and tsunamogenic potential with the prominent examples of the 2004, 2005 and 2007 earthquakes. Here, we invert travel time data of local earthquakes for vp and vp/vs velocity models of the central Sumatran forearc. Data were acquired by an amphibious seismometer network consisting of 52 land stations and 10 ocean bottom seismometers located on a segment of the Sumatran subduction zone that had not ruptured in a great earthquake since 1797 but witnessed recent ruptures to the north in 2005 (Nias earthquake, Mw = 8.7) and to the south in 2007 (Bengkulu earthquake, Mw = 8.5). 2D and 3D vp velocity anomalies reveal the downgoing slab and the sedimentary basins. Although the seismicity pattern in the study area appears to be strongly influenced by the obliquely subducting Investigator Fracture Zone to at least 200 km depth, the 3D velocity model shows prevailing trench parallel structures at depths of the plate interface. The tomographic model suggests a thinned crust below the basin east of the forearc islands (Nias, Pulau Batu, Siberut) at ~ 180 km distance to the trench. Vp velocities beneath the magmatic arc and the Sumatran fault zone SFZ are around 5 km/s at 10 km depth and the vp/vs ratios in the uppermost 10 km are low, indicating the presence of felsic lithologies typical for continental crust. We find moderately elevated vp/vs values of 1.85 at ~ 150 km distance to the trench in the region of the Mentawai fault. Vp/vs ratios suggest absence of large scale alteration of the mantle wedge and might explain why the seismogenic plate interface (observed as a locked zone from geodetic data) extends below the continental forearc Moho in Sumatra. Reduced vp velocities beneath the forearc basin covering the region between Mentawai Islands and the Sumatra mainland possibly reflect a reduced thickness of the overriding crust.

Description: Biological N2 fixation is an important part of the marine nitrogen cycle as it provides a source of new nitrogen that can support biological carbon export and sequestration. Research in the past decade has focused on determining the patterns of distribution and abundance of diazotrophs, defining the environmental features leading to these patterns and characterizing the factors that constrain marine N2 fixation overall. In this Review, we describe how variations in the deposition of iron from dust to different ocean basins affects the limiting nutrient for N2 fixation and the distribution of different diazotrophic species. However, many questions remain about marine N2 fixation, including the role of temperature, fixed nitrogen species, CO2 and physical forcing in controlling N2 fixation, as well as the potential for heterotrophic N2 fixation.

Description: In this study ship- and AUV-based multibeam data from the German Mn-nodule license area in the Clarion-Clipperton Zone (CCZ; eastern Pacific) are linked to ground truth data from optical imaging. Photographs obtained by an AUV enable semi-quantitative assessments of nodule coverage at a spatial resolution in the range of meters. Together with high resolution AUV bathymetry this revealed a correlation of small-scale terrain variations (〈 5 m horizontally, 〈 1 m vertically) with nodule abundance. In the presented data set, increased nodule coverage could be correlated with slopes 〉 1.8° and concave terrain. On a more regional scale, factors such as the geological setting (existence of horst and graben structures, sediment thickness, outcropping basement) and influence of bottom currents seem to play an essential role for the spatial variation of nodule abundance and the related hard substrate habitat. AUV imagery was also successfully employed to map the distribution of re-settled sediment following a disturbance and sediment cloud generation during a sampling deployment of an Epibenthic Sledge. Data from before and after the "disturbance" allows a direct assessment of the impact. Automated image processing analyzed the nodule coverage at the seafloor, revealing nodule blanketing by resettling of suspended sediment within 16 hours after the disturbance. The visually detectable impact was spatially limited to a maximum of 100m distance from the disturbance track, downstream of the bottom water current. A correlation with high resolution AUV bathymetry reveals that the blanketing pattern varies in extent by tens of meters, strictly following the bathymetry, even in areas of only slightly undulating seafloor (〈 1 m vertical change). These results highlight the importance of detailed terrain knowledge when engaging in resource assessment studies for nodule abundance estimates and defining minable areas. At the same time, it shows the importance of high resolution mapping for detailed benthic habitat studies that show a heterogeneity at scales of 10 m to 100 m. Terrain knowledge is also needed to determine the scale of the impact by seafloor sediment blanketing during mining-operations.

Description: Climate-driven changes in environmental conditions have significant and complex effects on marine ecosystems. Variability in phytoplankton elements and biochemicals can be important for global ocean biogeochemistry and ecological functions, while there is currently limited understanding on how elements and biochemicals respond to the changing environments in key coccolithophore species such as Emiliania huxleyi. We investigated responses of elemental stoichiometry and fatty acids (FAs) in a strain of E. huxleyi under three temperatures (12, 18 and 24 °C), three N : P supply ratios (molar ratios 10:1, 24:1 and 63:1) and two pCO2 levels (560 and 2400 µatm). Overall, C : N : P stoichiometry showed the most pronounced response to N : P supply ratios, with high ratios of particulate organic carbon vs. particulate organic nitrogen (POC : PON) and low ratios of PON vs. particulate organic phosphorus (PON : POP) in low-N media, and high POC : POP and PON : POP in low-P media. The ratio of particulate inorganic carbon vs. POC (PIC : POC) and polyunsaturated fatty acid proportions strongly responded to temperature and pCO2, both being lower under high pCO2 and higher with warming. We observed synergistic interactions between warming and nutrient deficiency (and high pCO2) on elemental cellular contents and docosahexaenoic acid (DHA) proportion in most cases, indicating the enhanced effect of warming under nutrient deficiency (and high pCO2). Our results suggest differential sensitivity of elements and FAs to the changes in temperature, nutrient availability and pCO2 in E. huxleyi, which is to some extent unique compared to non-calcifying algal classes. Thus, simultaneous changes of elements and FAs should be considered when predicting future roles of E. huxleyi in the biotic-mediated connection between biogeochemical cycles, ecological functions and climate change.

Description: The pre-industrial millennium is among the periods selected by the Paleoclimate Model Intercomparison Project (PMIP) for experiments contributing to the sixth phase of the Coupled Model Intercomparison Project (CMIP6) and the fourth phase of the PMIP (PMIP4). The past1000 transient simulations serve to investigate the response to (mainly) natural forcing under background conditions not too different from today, and to discriminate between forced and internally generated variability on interannual to centennial timescales. This paper describes the motivation and the experimental set-ups for the PMIP4-CMIP6 past1000 simulations, and discusses the forcing agents orbital, solar, volcanic, and land use/land cover changes, and variations in greenhouse gas concentrations. The past1000 simulations covering the pre-industrial millennium from 850 Common Era (CE) to 1849 CE have to be complemented by historical simulations (1850 to 2014 CE) following the CMIP6 protocol. The external forcings for the past1000 experiments have been adapted to provide a seamless transition across these time periods. Protocols for the past1000 simulations have been divided into three tiers. A default forcing data set has been defined for the Tier 1 (the CMIP6 past1000) experiment. However, the PMIP community has maintained the flexibility to conduct coordinated sensitivity experiments to explore uncertainty in forcing reconstructions as well as parameter uncertainty in dedicated Tier 2 simulations. Additional experiments (Tier 3) are defined to foster collaborative model experiments focusing on the early instrumental period and to extend the temporal range and the scope of the simulations. This paper outlines current and future research foci and common analyses for collaborative work between the PMIP and the observational communities (reconstructions, instrumental data).

Description: Although mesoscale ocean eddies are ubiquitous in the Southern Ocean, their average regional and seasonal association with phytoplankton has not been quantified systematically yet. To this end, we identify over 100000 mesoscale eddies with diameters of 50km and more in the Southern Ocean and determine the associated phytoplankton biomass anomalies using satellite-based chlorophyll-a (Chl) as a proxy. The mean Chl anomalies, δChl, associated with these eddies, comprising the upper echelon of the oceanic mesoscale, exceed ±10% over wide regions. The structure of these anomalies is largely zonal, with cyclonic, thermocline lifted, eddies having positive anomalies in the subtropical waters north of the Antarctic Circumpolar Current (ACC) and negative anomalies along its main flow path. The pattern is similar, but reversed for anticyclonic, thermocline deepened eddies. The seasonality of δChl is weak in subtropical waters, but pronounced along the ACC, featuring a seasonal sign switch. The spatial structure and seasonality of the mesoscale δChl can be explained largely by lateral advection, especially local eddy-stirring. A prominent exception is the ACC region in winter, where δChl is consistent with a modulation of phytoplankton light exposure caused by an eddy-induced modification of the mixed layer depth. The clear impact of mesoscale eddies on phytoplankton may implicate a downstream effect on Southern Ocean biogeochemical properties, such as mode water nutrient contents.

Description: Benthic foraminifera from Bottsand coastal lagoon, western Baltic Sea, have been studied since the mid-1960s. They were monitored annually in late autumn since 2003 at the terminal ditch of the lagoon. There were 12 different species recognised, of which three have not been recorded during earlier investigations. Dominant species showed strong interannual fluctuations and a steady increase in population densities over the last decade. Elphidium incertum, a stenohaline species of the Baltic deep water fauna, colonised the Bottsand lagoon in 2016, most likely during a period of salinities 〉19 units and water temperatures of 18 °C on average in early autumn. The high salinities probably triggered their germination from a propagule bank in the ditch bottom sediment. The new E. incertum population showed densities higher by an order of magnitude than those of the indigenous species. The latter did not decline, revealing that E. incertum used another food source or occupied a different microhabitat. Elphidium incertum survived transient periods of lower salinities in late autumn 2017, though with reduced abundances, and became a regular faunal constituent at the Bottsand lagoon.

Description: The quantitative reconstruction of past seawater salinity has yet to be achieved and the search for a direct and independent salinity proxy is ongoing. Recent culture and field studies show a significant positive correlation of Na/Ca with salinity in benthic and planktonic foraminiferal calcite. For accurate paleoceanographic reconstructions, consistent and reliable calibrations are necessary, which are still missing. In order to assess the reliability of foraminiferal Na/Ca as a direct proxy for seawater salinity, this study presents electron microprobe Na/Ca data, measured on cultured specimens of Trilobatus sacculifer. The culture experiments were conducted over a wide salinity range of 26 to 45, while temperature was kept constant. To further understand potential controlling factors of Na incorporation, measurements were also performed on foraminifera cultured at various temperatures in the range of 19.5 °C to 29.5 °C under constant salinity conditions. Foraminiferal Na/Ca ratios positively correlate with seawater salinity (Na/Caforam = 0.97 + 0.115 ⋅ Salinity, R = 0.97, p 〈 0.005). Temperature on the other hand exhibits no statistically significant relationship with Na/Ca ratios indicating salinity to be the dominant factor controlling Na incorporation. The culturing results are corroborated by measurements on T. sacculifer from Caribbean and Gulf of Guinea surface sediments. In conclusion, planktonic foraminiferal Na/Ca can be applied as a reliable proxy for reconstructing sea surface salinities, albeit species-specific calibrations might be necessary.

Description: We report here the results of total mercury (HgT) determinations along the 2014 Geotraces Geovide cruise (GA01 transect) in the North Atlantic Ocean (NA) from Lisbon (Portugal) to the coast of Labrador (Canada). HgT concentrations in unfiltered samples (HgTUNF) were log-normally distributed and ranged between 0.16 and 1.54 pmol L−1, with a geometric mean of 0.51 pmol L−1 for the 535 samples analysed. The dissolved fraction (〈 0.45 µm) of HgT (HgTF), determined on 141 samples, averaged 78 % of the HgTUNF for the entire data set, 84 % for open seawaters (below 100 m) and 91 % if the Labrador Sea data are excluded, where the primary production was high (with a winter convection down to 1400 m). HgTUNF concentrations increased eastwards and with depth from Greenland to Europe and from subsurface to bottom waters. The HgTUNF concentrations were similarly low in the subpolar gyre waters ( ∼  0.45 pmol L−1), whereas they exceeded 0.60 pmol L−1 in the subtropical gyre waters. The HgTUNF distribution mirrored that of dissolved oxygen concentration, with highest concentration levels associated with oxygen-depleted zones. The relationship between HgTF and the apparent oxygen utilization confirms the nutrient-like behaviour of Hg in the NA. An extended optimum multiparameter analysis allowed us to characterize HgTUNF concentrations in the different source water types (SWTs) present along the transect. The distribution pattern of HgTUNF, modelled by the mixing of SWTs, show Hg enrichment in Mediterranean waters and North East Atlantic Deep Water and low concentrations in young waters formed in the subpolar gyre and Nordic seas. The change in anthropogenic Hg concentrations in the Labrador Sea Water during its eastward journey suggests a continuous decrease in Hg content in this water mass over the last decades. Calculation of the water transport driven by the Atlantic Meridional Overturning Circulation across the Portugal–Greenland transect indicates northward Hg transport within the upper limb and southward Hg transport within the lower limb, with resulting net northward transport of about 97.2 kmol yr−1.

Description: The chemistry–climate model ECHAM-HAMMOZ contains a detailed representation of tropospheric and stratospheric reactive chemistry and state-of-the-art parameterizations of aerosols using either a modal scheme (M7) or a bin scheme (SALSA). This article describes and evaluates the model version ECHAM6.3-HAM2.3-MOZ1.0 with a focus on the tropospheric gas-phase chemistry. A 10-year model simulation was performed to test the stability of the model and provide data for its evaluation. The comparison to observations concentrates on the year 2008 and includes total column observations of ozone and CO from IASI and OMI, Aura MLS observations of temperature, HNO3, ClO, and O3 for the evaluation of polar stratospheric processes, an ozonesonde climatology, surface ozone observations from the TOAR database, and surface CO data from the Global Atmosphere Watch network. Global budgets of ozone, OH, NOx, aerosols, clouds, and radiation are analyzed and compared to the literature. ECHAM-HAMMOZ performs well in many aspects. However, in the base simulation, lightning NOx emissions are very low, and the impact of the heterogeneous reaction of HNO3 on dust and sea salt aerosol is too strong. Sensitivity simulations with increased lightning NOx or modified heterogeneous chemistry deteriorate the comparison with observations and yield excessively large ozone budget terms and too much OH. We hypothesize that this is an impact of potential issues with tropical convection in the ECHAM model.

Description: Lakes and reservoirs play an important role in the carbon cycle, and therefore monitoring their metabolic rates is essential. The triple oxygen-isotope anomaly of dissolved O2 [17Δ =  ln(1+δ17O) − 0.518  ×  ln(1 + δ18O)] offers a new, in situ, perspective on primary production, yet little is known about 17Δ from freshwater systems. We investigated the 17Δ together with the oxygen : argon ratio [Δ(O2 ∕ Ar)] in the subtropical Feitsui Reservoir in Taiwan from June 2014 to July 2015. Here, we present the seasonal variations in 17Δ, GP (gross production), NP (net production) and the NP ∕ GP (net to gross ratio) in association with environmental parameters. The 17Δ varied with depth and season, with values ranging between 26 and 205 per meg. The GP rates were observed to be higher (702 ± 107 mg C m−2 d−1) in winter than those (303 ± 66 mg C m−2 d−1) recorded during the summer. The overall averaged GP was 220 g C m−2 yr−1 and NP was −3 g C m−2 yr−1, implying the reservoir was net heterotrophic on an annual basis. This is due to negative NP rates from October to February (−198 ± 78 mg C m−2 d−1). Comparisons between GP rates obtained from the isotope mass balance approach and 14C bottle incubation method (14C–GP) showed consistent values on the same order of magnitude with a GP ∕ 14C–GP ratio of 1.2 ± 1.1. Finally we noted that, although typhoon occurrences were scarce, higher than average 17Δ values and GP rates were recorded after typhoon events.

Description: Oxygen optode measurements on floats and gliders suffer from a slow time response and various sources of drift in the calibration coefficients. Based on two dual-O2 Argo floats, we show how to post-correct for the effect of the optode's time response and give an update on optode in situ drift stability and in-air calibration. Both floats are equipped with an unpumped Aanderaa 4330 optode and a pumped Sea-Bird SBE63 optode. Response times for the pumped SBE63 were derived following Bittig et al. (2014) and the same methods were used to correct the time response bias. Using both optodes on each float, the time response regime of the unpumped Aanderaa optode was characterized more accurately than previously possible. Response times for the pumped SBE63 on profiling floats are in the range of 25–40 s, while they are between 60 and 95 s for the unpumped 4330 optode. Our parameterization can be employed to post-correct the slow optode time response on floats and gliders. After correction, both sensors agree to within 2–3 µmol kg−1 (median difference) in the strongest gradients (120 µmol kg−1 change over 8 min or 20 dbar) and better elsewhere. However, time response correction is only possible if measurement times are known, i.e., provided by the platform as well as transmitted and stored with the data. The O2 in-air measurements show a significant in situ optode drift of −0.40 and −0.27 % yr−1 over the available 2 and 3 years of deployment, respectively. Optode in-air measurements are systematically biased high during midday surfacings compared to dusk, dawn, and nighttime. While preference can be given to nighttime surfacings to avoid this in-air calibration bias, we suggest a parameterization of the daytime effect as a function of the Sun's elevation to be able to use all data and to better constrain the result. Taking all effects into account, calibration factors have an uncertainty of 0.1 %. In addition, in-air calibration factors vary by 0.1–0.2 % when using different reanalysis models as a reference. The overall accuracy that can be achieved following the proposed correction routines is better than 1 µmol kg−1.

Description: Conventional integration of Earth system and ocean models can accrue considerable computational expenses, particularly for marine biogeochemical applications. "Offline" numerical schemes in which only the biogeochemical tracers are time stepped and transported using a pre-computed circulation field can substantially reduce the burden and are thus an attractive alternative. One such scheme is the "transport matrix method" (TMM), which represents tracer transport as a sequence of sparse matrix–vector products that can be performed efficiently on distributed-memory computers. While the TMM has been used for a variety of geochemical and biogeochemical studies, to date the resulting solutions have not been comprehensively assessed against their "online" counterparts. Here, we present a detailed comparison of the two. It is based on simulations of the state-of-the-art biogeochemical sub-model embedded within the widely used coarse-resolution University of Victoria Earth System Climate Model (UVic ESCM). The default, non-linear advection scheme was first replaced with a linear, third-order upwind-biased advection scheme to satisfy the linearity requirement of the TMM. Transport matrices were extracted from an equilibrium run of the physical model and subsequently used to integrate the biogeochemical model offline to equilibrium. The identical biogeochemical model was also run online. Our simulations show that offline integration introduces some bias to biogeochemical quantities through the omission of the polar filtering used in UVic ESCM and in the offline application of time-dependent forcing fields, with high latitudes showing the largest differences with respect to the online model. Differences in other regions and in the seasonality of nutrients and phytoplankton distributions are found to be relatively minor, giving confidence that the TMM is a reliable tool for offline integration of complex biogeochemical models. Moreover, while UVic ESCM is a serial code, the TMM can be run on a parallel machine with no change to the underlying biogeochemical code, thus providing orders of magnitude speed-up over the online model.

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

Description: Although rising global sea levels will affect the shape of coastlines over the coming decades1, 2, the most severe and catastrophic shoreline changes occur as a consequence of local and regional-scale processes. Changes in sediment supply3 and deltaic subsidence4, 5, both natural or anthropogenic, and the occurrences of tropical cyclones4, 5 and tsunamis6 have been shown to be the leading controls on coastal erosion. Here, we use satellite images of South American mangrove-colonized mud banks collected over the past twenty years to reconstruct changes in the extent of the shoreline between the Amazon and Orinoco rivers. The observed timing of the redistribution of sediment and migration of the mud banks along the 1,500 km muddy coast suggests the dominant control of ocean forcing by the 18.6 year nodal tidal cycle7. Other factors affecting sea level such as global warming or El Niño and La Niña events show only secondary influences on the recorded changes. In the coming decade, the 18.6 year cycle will result in an increase of mean high water levels of 6 cm along the coast of French Guiana, which will lead to a 90 m shoreline retreat.

Description: The assessment of the ocean biota's role in climate climate change is often carried out with global biogeochemical ocean models that contain many components, and involve a high level of parametric uncertainty. Examination the models' fit to climatologies of inorganic tracers, after the models have been spun up to steady state, is a common, but computationally expensive procedure to assess model performance and reliability. Using new tools that have become available for global model assessment and calibration in steady state, this paper examines two different model types – a complex seven-component model (MOPS), and a very simple two-component model (RetroMOPS) – for their fit to dissolved quantities. Before comparing the models, a subset of their biogeochemical parameters has been optimised against annual mean nutrients and oxygen. Both model types fit the observations almost equally well. The simple model, which contains only nutrients and dissolved organic phosphorus (DOP), is sensitive to the parameterisation of DOP production and decay. The spatio-temporal decoupling of nitrogen and oxygen, and processes involved in their uptake and release, renders oxygen and nitrate valuable tracers for model calibration. In addition, the non-conservative nature of these tracers (with respect to their upper boundary condition) introduces the global bias as a useful additional constraint on model parameters. Dissolved organic phosphorous at the surface behaves antagonistically to phosphate, and suggests that observations of this tracer – although difficult to measure – may be an important asset for model calibration

Description: Benthic microbial methanogenesis is a known source of methane in marine systems. In most sediments, the majority of methanogenesis is located below the sulfate-reducing zone, as sulfate reducers outcompete methanogens for the major substrates hydrogen and acetate. The coexistence of methanogenesis and sulfate reduction has been shown before and is possible through the usage of noncompetitive substrates by methanogens such as methanol or methylated amines. However, knowledge about the magnitude, seasonality, and environmental controls of this noncompetitive methane production is sparse. In the present study, the presence of methanogenesis within the sulfate reduction zone (SRZ methanogenesis) was investigated in sediments (0–30 cm below seafloor, cm b.s.f.) of the seasonally hypoxic Eckernförde Bay in the southwestern Baltic Sea. Water column parameters such as oxygen, temperature, and salinity together with porewater geochemistry and benthic methanogenesis rates were determined in the sampling area "Boknis Eck" quarterly from March 2013 to September 2014 to investigate the effect of seasonal environmental changes on the rate and distribution of SRZ methanogenesis, to estimate its potential contribution to benthic methane emissions, and to identify the potential methanogenic groups responsible for SRZ methanogenesis. The metabolic pathway of methanogenesis in the presence or absence of sulfate reducers, which after the addition of a noncompetitive substrate was studied in four experimental setups: (1) unaltered sediment batch incubations (net methanogenesis), (2) 14C-bicarbonate labeling experiments (hydrogenotrophic methanogenesis), (3) manipulated experiments with the addition of either molybdate (sulfate reducer inhibitor), 2-bromoethanesulfonate (methanogen inhibitor), or methanol (noncompetitive substrate, potential methanogenesis), and (4) the addition of 13C-labeled methanol (potential methylotrophic methanogenesis). After incubation with methanol, molecular analyses were conducted to identify key functional methanogenic groups during methylotrophic methanogenesis. To also compare the magnitudes of SRZ methanogenesis with methanogenesis below the sulfate reduction zone (〉 30 cm b.s.f.), hydrogenotrophic methanogenesis was determined by 14C-bicarbonate radiotracer incubation in samples collected in September 2013. SRZ methanogenesis changed seasonally in the upper 30 cm b.s.f. with rates increasing from March (0.2 nmol cm−3 d−1) to November (1.3 nmol cm−3 d−1) 2013 and March (0.2 nmol cm−3 d−1) to September (0.4 nmol cm−3 d−1) 2014. Its magnitude and distribution appeared to be controlled by organic matter availability, C / N, temperature, and oxygen in the water column, revealing higher rates in the warm, stratified, hypoxic seasons (September–November) compared to the colder, oxygenated seasons (March–June) of each year. The majority of SRZ methanogenesis was likely driven by the usage of noncompetitive substrates (e.g., methanol and methylated compounds) to avoid competition with sulfate reducers, as was indicated by the 1000–3000-fold increase in potential methanogenesis activity observed after methanol addition. Accordingly, competitive hydrogenotrophic methanogenesis increased in the sediment only below the depth of sulfate penetration (〉 30 cm b.s.f.). Members of the family Methanosarcinaceae, which are known for methylotrophic methanogenesis, were detected by PCR using Methanosarcinaceae-specific primers and are likely to be responsible for the observed SRZ methanogenesis. The present study indicates that SRZ methanogenesis is an important component of the benthic methane budget and carbon cycling in Eckernförde Bay. Although its contributions to methane emissions from the sediment into the water column are probably minor, SRZ methanogenesis could directly feed into methane oxidation above the sulfate–methane transition zone.

Description: The balance in microbial net consumption of nitrogen and phosphorus was investigated in samples collected in two mesotrophic coastal environments: the Baltic Sea (Tvärminne field station) and the North Sea (Espegrend field station). For this, we have refined a bioassay based on the response in alkaline phosphatase activity (APA) over a matrix of combinations in nitrogen and phosphorus additions. This assay not only provides information on which element (N or P) is the primary limiting nutrient, but also gives a quantitative estimate for the excess of the secondary limiting element (P+ or N+, respectively), as well as the ratio of balanced net consumption of added N and P over short timescales (days). As expected for a Baltic Sea late spring–early summer situation, the Tvärminne assays (n =  5) indicated N limitation with an average P+ =  0.30 ± 0.10 µM-P, when incubated for 4 days. For short incubations (1–2 days), the Espegrend assays indicated P limitation, but the shape of the response surface changed with incubation time, resulting in a drift in parameter estimates toward N limitation. Extrapolating back to zero incubation time gave P limitation with N+ ≈  0.9 µM-N. The N : P ratio (molar) of nutrient net consumption varied considerably between investigated locations: from 2.3 ± 0.4 in the Tvärminne samples to 13 ± 5 and 32 ± 3 in two samples from Espegrend. Our assays included samples from mesocosm acidification experiments, but statistically significant effects of ocean acidification were not found by this method.

Description: THE rare deep-sea octopod Cirrothauma murrayi Chun 1910 was first described from a single specimen caught during the Michael Sars Expedition of 1910 (ref. 1). Until now it has been caught only four more times2. We describe here three specimens of this species that were recently caught during biological cruises of RRS Discovery (Fig. 1). All of these animals, including the Discovery ones, have been caught at depths of more than 1,500 m, except one that was dip-netted through the ice of the Arctic Ocean3.

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

Description: The Pacific sector of Antarctica, including both the Antarctic Peninsula and continental West Antarctica, has experienced substantial warming in the past 30 years. An increase in the circumpolar westerlies, owing in part to the decline in stratospheric ozone concentrations since the late 1970s, may account for warming trends in the peninsula region in austral summer and autumn. The more widespread warming in continental West Antarctica (Ellsworth Land and Marie Byrd Land) occurs primarily in austral winter and spring, and remains unexplained. Here we use observations of Antarctic surface temperature and global sea surface temperature, and atmospheric circulation data to show that recent warming in continental West Antarctica is linked to sea surface temperature changes in the tropical Pacific. Over the past 30 years, anomalous sea surface temperatures in the central tropical Pacific have generated an atmospheric Rossby wave response that influences atmospheric circulation over the Amundsen Sea, causing increased advection of warm air to the Antarctic continent. General circulation model experiments show that the central tropical Pacific is a critical region for producing the observed high latitude response. We conclude that, by affecting the atmospheric circulation at high southern latitudes, increasing tropical sea surface temperatures may account for West Antarctic warming through most of the twentieth century.

Description: Here we show the use of the 210Pb-226Ra excess method to determine the growth rate of two corals from the world's largest known cold-water coral reef, Røst Reef, north of the Arctic circle off Norway. Colonies of each of the two species that build the reef, Lophelia pertusa and Madrepora oculata, were collected alive at 350 m depth using a submersible. Pb and Ra isotopes were measured along the major growth axis of both specimens using low level alpha and gamma spectrometry and trace element compositions were studied. 210Pb and 226Ra differ in the way they are incorporated into coral skeletons. Hence, to assess growth rates, we considered the exponential decrease of initially incorporated 210Pb, as well as the increase in 210Pb from the decay of 226Ra and contamination with 210Pb associated with Mn-Fe coatings that we were unable to remove completely from the oldest parts of the skeletons. 226Ra activity was similar in both coral species, so, assuming constant uptake of 210Pb through time, we used the 210Pb-226Ra chronology to calculate growth rates. The 45.5 cm long branch of M. oculata was 31 yr with an average linear growth rate of 14.4 ± 1.1 mm yr−1 (2.6 polyps per year). Despite cleaning, a correction for Mn-Fe oxide contamination was required for the oldest part of the colony; this correction corroborated our radiocarbon date of 40 yr and a mean growth rate of 2 polyps yr−1. This rate is similar to the one obtained in aquarium experiments under optimal growth conditions.

Description: The Ocean Model Intercomparison Project (OMIP) focuses on the physics and biogeochemistry of the ocean component of Earth system models participating in the sixth phase of the Coupled Model Intercomparison Project (CMIP6). OMIP aims to provide standard protocols and diagnostics for ocean models, while offering a forum to promote their common assessment and improvement. It also offers to compare solutions of the same ocean models when forced with reanalysis data (OMIP simulations) vs. when integrated within fully coupled Earth system models (CMIP6). Here we detail simulation protocols and diagnostics for OMIP's biogeochemical and inert chemical tracers. These passive-tracer simulations will be coupled to ocean circulation models, initialized with observational data or output from a model spin-up, and forced by repeating the 1948–2009 surface fluxes of heat, fresh water, and momentum. These so-called OMIP-BGC simulations include three inert chemical tracers (CFC-11, CFC-12, SF6) and biogeochemical tracers (e.g., dissolved inorganic carbon, carbon isotopes, alkalinity, nutrients, and oxygen). Modelers will use their preferred prognostic BGC model but should follow common guidelines for gas exchange and carbonate chemistry. Simulations include both natural and total carbon tracers. The required forced simulation (omip1) will be initialized with gridded observational climatologies. An optional forced simulation (omip1-spunup) will be initialized instead with BGC fields from a long model spin-up, preferably for 2000 years or more, and forced by repeating the same 62-year meteorological forcing. That optional run will also include abiotic tracers of total dissolved inorganic carbon and radiocarbon, CTabio and 14CTabio, to assess deep-ocean ventilation and distinguish the role of physics vs. biology. These simulations will be forced by observed atmospheric histories of the three inert gases and CO2 as well as carbon isotope ratios of CO2. OMIP-BGC simulation protocols are founded on those from previous phases of the Ocean Carbon-Cycle Model Intercomparison Project. They have been merged and updated to reflect improvements concerning gas exchange, carbonate chemistry, and new data for initial conditions and atmospheric gas histories. Code is provided to facilitate their implementation.

Description: The injection of sulphur into the stratosphere by explosive volcanic eruptions is the cause of significant climate variability. Based on sulphate records from a suite of ice cores from Greenland and Antarctica, the eVolv2k database includes estimates of the magnitudes and approximate source latitudes of major volcanic stratospheric sulphur injection (VSSI) events from 500 BCE to 1900 CE, constituting an update of prior reconstructions and an extension of the record by 1000 years. The VSSI estimates incorporate improvements to the ice core records in terms of synchronization and dating, refinements to the methods used to estimate VSSI from ice core records, and includes first estimates of the random uncertainties in VSSI values. VSSI estimates for many of the largest eruptions, including Samalas (1257), Tambora (1815) and Laki (1783) are within 10% of prior estimates. A number of strong events are included in eVolv2k which are largely underestimated or not included in earlier VSSI reconstructions, including events in 540, 574, 682 and 1108 CE. The long term annual mean VSSI from major volcanic eruptions is estimated to be ∼ 0.5 Tg [S] yr−1, ∼ 50 % greater than a prior reconstruction, due to the identification of more events and an increase in the magnitude of many intermediate events. A long-term, latitudinally and monthly resolved stratospheric aerosol optical depth (SAOD) time series is reconstructed from the eVolv2k VSSI estimates, and the resulting global mean SAOD is found to be similar (within 33%) to a prior reconstruction for most of the largest eruptions. The long-term (500 BCE–900 CE) average global mean SAOD estimated from the eVolv2k VSSI estimates and including a constant "background" injection of stratospheric sulphur is ∼ 0.014, 30 % greater than a prior reconstruction. These new long-term reconstructions of past VSSI and SAOD variability give context to recent volcanic forcing, suggesting that the 20th century was a period of somewhat weaker than average volcanic forcing, with current best estimates of 20th century mean VSSI and SAOD values being 25 and 14 % less, respectively, than the mean of the 500 BCE to 1900 CE period. The reconstructed VSSI and SAOD data are available at https://doi.org/10.1594/WDCC/eVolv2k_v2〉.

Description: The Cretaceous ocean witnessed intervals of profound perturbations such as volcanic input of large amounts of CO2, anoxia, eutrophication and introduction of biologically relevant metals. Some of these extreme events were characterized by size reduction and/or morphological changes of a few calcareous nannofossil species. The correspondence between intervals of high trace metal concentrations and coccolith dwarfism suggests a negative effect of these elements on nannoplankton biocalcification processes in past oceans. In order to test this hypothesis, we explored the potential effect of a mixture of trace metals on growth and morphology of four living coccolithophore species, namely Emiliania huxleyi, Gephyrocapsa oceanica, Pleurochrysis carterae and Coccolithus pelagicus. The phylogenetic history of coccolithophores shows that the selected living species are linked to Mesozoic species showing dwarfism under excess metal concentrations. The trace metals tested were chosen to simulate the environmental stress identified in the geological record and upon known trace metal interactions with living coccolithophore algae. Our laboratory experiments demonstrated that elevated trace metal concentrations, similarly to the fossil record, affect coccolithophore algae size and/or weight. Smaller coccoliths were detected in E. huxleyi and C. pelagicus, while coccoliths of G. oceanica showed a decrease in size only at the highest trace metal concentrations. P. carterae coccolith size was unresponsive to changing trace metal concentrations. These differences among species allow discriminating the most- (P. carterae), intermediate- (E. huxleyi and G. oceanica) and least-tolerant (C. pelagicus) taxa. The fossil record and the experimental results converge on a selective response of coccolithophores to metal availability. These species-specific differences must be considered before morphological features of coccoliths are used to reconstruct paleo-chemical conditions.

Description: The Baltic Sea is a seasonally ice-covered marginal sea in northern Europe with intense wintertime ship traffic and a sensitive ecosystem. Understanding and modeling the evolution of the sea-ice pack is important for climate effect studies and forecasting purposes. Here we present and evaluate the sea-ice component of a new NEMO–LIM3.6-based ocean–sea-ice setup for the North Sea and Baltic Sea region (NEMO-Nordic). The setup includes a new depth-based fast-ice parametrization for the Baltic Sea. The evaluation focuses on long-term statistics, from a 45-year long hindcast, although short-term daily performance is also briefly evaluated. We show that NEMO-Nordic is well suited for simulating the mean sea-ice extent, concentration, and thickness as compared to the best available observational data set. The variability of the annual maximum Baltic Sea ice extent is well in line with the observations, but the 1961–2006 trend is underestimated. Capturing the correct ice thickness distribution is more challenging. Based on the simulated ice thickness distribution we estimate the undeformed and deformed ice thickness and concentration in the Baltic Sea, which compares reasonably well with observations.

Description: We developed a new method for the determination of dissolved nitric oxide (NO) in discrete seawater samples based on the combination of a purge-and-trap setup and a fluorometric detection of NO. 2,3-diaminonaphthalene (DAN) reacts with NO in seawater to form the highly fluorescent 2,3-naphthotriazole (NAT). The fluorescence intensity was linear for NO concentrations in the range from 0.14 to 19 nmol L−1. We determined a detection limit of 0.068 nmol L−1, an average recovery coefficient of 83.8 % (80.2–90.0 %), and a relative standard deviation of ±7.2 %. With our method we determined for the first time the temporal and spatial distributions of NO surface concentrations in coastal waters of the Yellow Sea off Qingdao and in Jiaozhou Bay during a cruise in November 2009. The concentrations of NO varied from below the detection limit to 0.50 nmol L−1 with an average of 0.26 ± 0.14 nmol L−1. NO surface concentrations were generally enhanced significantly during daytime, implying that NO formation processes such as NO2− photolysis are much higher during daytime than chemical NO consumption, which, in turn, lead to a significant decrease in NO concentrations during nighttime. In general, NO surface concentrations and measured NO production rates were higher compared to previously reported measurements. This might be caused by the high NO2− surface concentrations encountered during the cruise. Moreover, additional measurements of NO production rates implied that the occurrence of particles and a temperature increase can enhance NO production rates. With the method introduced here, we have a reliable and comparably easy to use method at hand to measure oceanic NO surface concentrations, which can be used to decipher both its temporal and spatial distributions as well as its biogeochemical pathways in the oceans.

Description: Zircon is a common mineral in continental crustal rocks. As it is not easily altered in processes such as erosion or transport, this mineral is often used in the reconstruction of geological processes such as the formation and evolution of the continents. Zircon can also survive under conditions of the Earth’s mantle, and rare cases of zircons crystallizing in the mantle significantly before their entrainment into magma and eruption to the surface have been reported1,2,3. Here we analyse the isotopic and trace element compositions of large zircons of gem quality from the Eger rift, Bohemian massif, and find that they are derived from the mantle. (U–Th)/He analyses suggest that the zircons as well as their host basalts erupted between 29 and 24 million years ago, but fragments from the same xenocrysts reveal U–Pb ages between 51 and 83 million years. We note a lack of older volcanism and of fragments from the lower crust, which suggests that crustal residence time before eruption is negligible and that most rock fragments found in similar basalts from adjacent volcanic fields equilibrated under mantle conditions. We conclude that a specific chemical environment in this part of the Earth’s upper mantle allowed the zircons to remain intact for about 20–60 million years.

Description: The synthetic production of monodisperse single magnetic domain nanoparticles at ambient temperature is challenging. In nature, magnetosomes--membrane-bound magnetic nanocrystals with unprecedented magnetic properties--can be biomineralized by magnetotactic bacteria. However, these microbes are difficult to handle. Expression of the underlying biosynthetic pathway from these fastidious microorganisms within other organisms could therefore greatly expand their nanotechnological and biomedical applications. So far, this has been hindered by the structural and genetic complexity of the magnetosome organelle and insufficient knowledge of the biosynthetic functions involved. Here, we show that the ability to biomineralize highly ordered magnetic nanostructures can be transferred to a foreign recipient. Expression of a minimal set of genes from the magnetotactic bacterium Magnetospirillum gryphiswaldense resulted in magnetosome biosynthesis within the photosynthetic model organism Rhodospirillum rubrum. Our findings will enable the sustainable production of tailored magnetic nanostructures in biotechnologically relevant hosts and represent a step towards the endogenous magnetization of various organisms by synthetic biology.

Description: Igneous intrusions in sedimentary basins may have a profound effect on the thermal structure and physical properties of the hosting sedimentary rocks. These include mechanical effects such as deformation and uplift of sedimentary layers, generation of overpressure, mineral reactions and porosity evolution, and fracturing and vent formation following devolatilization reactions and the generation of CO2 and CH4. The gas generation and subsequent migration and venting may have contributed to several of the past climatic changes such as the end-Permian event and the Paleocene-Eocene Thermal Maximum. Additionally, the generation and expulsion of hydrocarbons and cracking of pre-existing oil reservoirs around a hot magmatic intrusion is of significant interest to the energy industry. In this paper, we present a user-friendly 1D FEM based tool, SILLi, which calculates the thermal effects of sill intrusions on the enclosing sedimentary stratigraphy. The model is accompanied by three case studies of sills emplaced in two different sedimentary basins, the Karoo Basin in South Africa and the Vøring Basin offshore Norway. Input data for the model is the present-day well log or sedimentary column with an Excel input file and includes rock parameters such as thermal conductivity, total organic carbon (TOC) content, porosity, and latent heats. The model accounts for sedimentation and burial based on a rate calculated by the sedimentary layer thickness and age. Erosion of the sedimentary column is also included to account for realistic basin evolution. Multiple sills can be emplaced within the system with varying ages. The emplacement of a sill occurs instantaneously. The model can be applied to volcanic sedimentary basins occurring globally. The model output includes the thermal evolution of the sedimentary column through time, and the changes that take place following sill emplacement such as TOC changes, thermal maturity, and the amount of organic and carbonate-derived CO2. The TOC and vitrinite results can be readily benchmarked within the tool to present-day values measured within the sedimentary column. This allows the user to determine the conditions required to obtain results that match observables and leads to a better understanding of metamorphic processes in sedimentary basins.

Description: Past warm periods provide an opportunity to evaluate climate models under extreme forcing scenarios, in particular high ( 〉  800 ppmv) atmospheric CO2 concentrations. Although a post hoc intercomparison of Eocene ( ∼  50  Ma) climate model simulations and geological data has been carried out previously, models of past high-CO2 periods have never been evaluated in a consistent framework. Here, we present an experimental design for climate model simulations of three warm periods within the early Eocene and the latest Paleocene (the EECO, PETM, and pre-PETM). Together with the CMIP6 pre-industrial control and abrupt 4 ×  CO2 simulations, and additional sensitivity studies, these form the first phase of DeepMIP – the Deep-time Model Intercomparison Project, itself a group within the wider Paleoclimate Modelling Intercomparison Project (PMIP). The experimental design specifies and provides guidance on boundary conditions associated with palaeogeography, greenhouse gases, astronomical configuration, solar constant, land surface processes, and aerosols. Initial conditions, simulation length, and output variables are also specified. Finally, we explain how the geological data sets, which will be used to evaluate the simulations, will be developed.

Description: THE Sierra Leone Rise, located in the east equatorial Atlantic, forms a discontinuous chain of seamounts as shallow as 2 km extending with a general NE–SW trend from near the Sierra Leone coast of Africa, to the St Paul fracture zone near the Mid-Atlantic Ridge (Fig. 1). The origin of this feature has remained a topic of discussion. Sheridan et al.1 have hypothesised that the Sierra Leone Rise is a volcanic structure formed at the beginning of the opening of the Atlantic in the early Cretaceous period. The twin features of the Sierra Leone and the Ceara Rises are probably of oceanic origin and were created 80 Myr ago or later in their present-day position with respect to Africa and South America2. The Atlantic ocean exhibits several similar aseismic structures which appear symmetrically oriented with respect to the mid-oceanic ridge, such as the Walvis–Rio Grande Rise and the Iceland Faeroes–Iceland Greenland Ridges. These structures are volcanic edifices having a composition similar to that found in their associated islands3–7. Deep sea drilling of the Ceara Rise8,9 penetrated a basaltic basement of the upper Cretaceous period (Maestrichtian) (Leg 39, Site 354). Similarly, a DSDP hole (Leg 41, Site 366) on the Sierra Leone Rise, penetrated sediments of the same period, without reaching basement10. We report here the discovery of alkali-rich volcanics in an area of the Sierra Leone Rise. The sediment overlying the rock fragments is aged ∼45 Myr.

Description: The Stratospheric Sulfur and its Role in Climate (SSiRC) interactive stratospheric aerosol model intercomparison project (ISA-MIP) explores uncertainties in the processes that connect volcanic emission of sulphur gas species and the radiative forcing associated with the resulting enhancement of the stratospheric aerosol layer. The central aim of ISA-MIP is to constrain and improve interactive stratospheric aerosol models and reduce uncertainties in the stratospheric aerosol forcing by comparing results of standardized model experiments with a range of observations. In this paper we present 4 co-ordinated inter-model experiments designed to investigate key processes which influence the formation and temporal development of stratospheric aerosol in different time periods of the observational record. The "Background" (BG) experiment will focus on microphysics and transport processes under volcanically quiescent conditions, when the stratospheric aerosol is controlled by the transport of aerosols and their precursors from the troposphere to the stratosphere. The "Transient Aerosol Record" (TAR) experiment will explore the role of small- to moderate-magnitude volcanic eruptions, anthropogenic sulphur emissions and transport processes over the period 1998–2012 and their role in the warming hiatus. Two further experiments will investigate the stratospheric sulphate aerosol evolution after major volcanic eruptions. The "Historical Eruptions SO2 Emission Assessment" (HErSEA) experiment will focus on the uncertainty in the initial emission of recent large-magnitude volcanic eruptions, while the "Pinatubo Emulation in Multiple models" (PoEMS) experiment will provide a comprehensive uncertainty analysis of the radiative forcing from the 1991 Mt. Pinatubo eruption.

Description: Marine phytoplankton blooms are annual spring events that sustain active and diverse bloom-associated bacterial populations. Blooms vary considerably in terms of eukaryotic species composition and environmental conditions, but a limited number of heterotrophic bacterial lineages — primarily members of the Flavobacteriia, Alphaproteobacteria and Gammaproteobacteria — dominate these communities. In this Review, we discuss the central role that these bacteria have in transforming phytoplankton-derived organic matter and thus in biogeochemical nutrient cycling. On the basis of selected field and laboratory-based studies of flavobacteria and roseobacters, distinct metabolic strategies are emerging for these archetypal phytoplankton-associated taxa, which provide insights into the underlying mechanisms that dictate their behaviours during blooms.

Description: Shorelines exhibit long-range dependence (LRD) and have been shown in some environments to be described in the wavenumber domain by a power law characteristic of scale-independence. Recent evidence suggests that the geomorphology of barrier islands can, however, exhibit scale dependence as a result of systematic variations of the underlying framework geology. The LRD of framework geology, which influences island geomorphology and its response to storms and sea level rise, has not been previously examined. Electromagnetic induction (EMI) surveys conducted along Padre Island National Seashore (PAIS), Texas, USA, reveal that the EMI apparent conductivity σa signal and, by inference, the framework geology exhibits LRD at scales up to 101 to 102 km. Our study demonstrates the utility of describing EMI σa and LiDAR spatial series by a fractional auto-regressive integrated moving average process that specifically models LRD. This method offers a robust and compact way for quantifying the geological variations along a barrier island shoreline using three parameters (p,d,q). We discuss how ARIMA (0,d,0) models that use a single parameter d provide a quantitative measure for determining free and forced barrier island evolutionary behavior across different scales. Statistical analyses at regional, intermediate, and local scales suggest that the geologic framework within an area of paleo-channels exhibits a first order control on dune height. The exchange of sediment amongst nearshore, beach and dune in areas outside this region are scale-independent, implying that barrier islands like PAIS exhibit a combination of free and forced behaviors that affect the response of the island to sea level rise.

Description: The influence of mesoscale eddies on the flow field and the water masses, especially the oxygen distribution of the eastern tropical South Pacific, is investigated from a mooring, float, and satellite data set. Two anticyclonic (ACE1/2), one mode-water (MWE), and one cyclonic eddy (CE) are identified and followed in detail with satellite data on their westward transition with velocities of 3.2 to 6.0cms−1 from their generation region, the shelf of the Peruvian and Chilean upwelling regime, across the Stratus Ocean Reference Station (ORS;  ∼ 20°S, 85°W) to their decaying region far west in the oligotrophic open ocean. The ORS is located in the transition zone between the oxygen minimum zone and the well oxygenated South Pacific subtropical gyre. Velocity, hydrographic, and oxygen measurements at the mooring show the impact of eddies on the weak flow region of the eastern tropical South Pacific. Strong anomalies are related to the passage of eddies and are not associated with a seasonal signal in the open ocean. The mass transport of the four observed eddies across 85°W is between 1.1 and 1.8Sv. The eddy type-dependent available heat, salt, and oxygen anomalies are 8.1×1018J (ACE2), 1.0×1018J (MWE), and −8.9×1018J (CE) for heat; 25.2×1010kg (ACE2), −3.1×1010kg (MWE), and −41.5×1010kg (CE) for salt; and −3.6×1016µmol (ACE2), −3.5×1016µmol (MWE), and −6.5×1016µmol (CE) for oxygen showing a strong imbalance between anticyclones and cyclones for salt transports probably due to seasonal variability in water mass properties in the formation region of the eddies. Heat, salt, and oxygen fluxes out of the coastal region across the ORS region in the oligotrophic open South Pacific are estimated based on these eddy anomalies and on eddy statistics (gained out of 23 years of satellite data). Furthermore, four profiling floats were trapped in the ACE2 during its westward propagation between the formation region and the open ocean, which allows for conclusions on lateral mixing of water mass properties with time between the core of the eddy and the surrounding water. The strongest lateral mixing was found between the seasonal thermocline and the eddy core during the first half of the eddy lifetime.

Description: RECENT advances in 40Ar/39Ar dating1,2 have made it possible to date individual K-feldspar grains from Pleistocene tephra, a capability that greatly improves the reliability and temporal resolving power of the method. Here we apply these new techniques to the dating of a phonolite tephra from the East Eifel volcanic field in West Germany, which is sandwiched between loess and palaeosol (alfisol) deposits, and which was therefore erupted during the transition from a glacial to an interglacial period. Our age estimate for this transition is 215±4 kyr (1 σ), which has important implications for the marine δ18O timescale and for models of global climate change during the Pleistocene. The results show that single-grain dating can detect and compensate for the large quantities of xenocrystic contaminants which are found in many tephra deposits. This technique could be used to date the tephra layers found in marine sediment cores and the results could greatly enhance the reliability of the marine δ18O timescale for more rigorous Fourier analysis testing of the Milankovitch hypothesis.

Description: The eruption of Mt. Tambora in 1815 was the largest volcanic eruption of the past 500 years. The eruption had significant climatic impacts, leading to the 1816 "year without a summer", and remains a valuable event from which to understand the climatic effects of large stratospheric volcanic sulfur dioxide injections. The eruption also resulted in one of the strongest and most easily identifiable volcanic sulfate signals in polar ice cores, which are widely used to reconstruct the timing and atmospheric sulfate loading of past eruptions. As part of the Model Intercomparison Project on the climatic response to Volcanic forcing (VolMIP), five state-of-the-art global aerosol models simulated this eruption. We analyse both simulated background (no Tambora) and volcanic (with Tambora) sulfate deposition to polar regions and compare to ice core records. The models simulate overall similar patterns of background sulfate deposition, al-though there are differences in regional details and magnitude. However, the volcanic sulfate deposition varies considerably between the models with differences in timing, spatial pattern and magnitude. Mean simulated deposited sulfate on Antarctica ranges from 19 to 264 kgkm-2 and on Greenland from 31 to 194 kgkm-2, as compared to the mean ice-corederived estimates of roughly 50 kgkm-2 for both Greenland and Antarctica. The ratio of the hemispheric atmospheric sulfate aerosol burden after the eruption to the average ice sheet deposited sulfate varies between models by up to a factor of 15. Sources of this inter-model variability include differences in both the formation and the transport of sulfate aerosol. Our results suggest that deriving relationships between sulfate deposited on ice sheets and atmospheric sulfate burdens from model simulations may be associated with greater uncertainties than previously thought.

Description: Oceanic very short-lived substances (VSLSs), such as bromoform (CHBr3), contribute to stratospheric halogen loading and, thus, to ozone depletion. However, the amount, timing, and region of bromine delivery to the stratosphere through one of the main entrance gates, the Indian summer monsoon circulation, are still uncertain. In this study, we created two bromoform emission inventories with monthly resolution for the tropical Indian Ocean and west Pacific based on new in situ bromoform measurements and novel ocean biogeochemistry modeling. The mass transport and atmospheric mixing ratios of bromoform were modeled for the year 2014 with the particle dispersion model FLEXPART driven by ERA-Interim reanalysis. We compare results between two emission scenarios: (1) monthly averaged and (2) annually averaged emissions. Both simulations reproduce the atmospheric distribution of bromoform from ship- and aircraft-based observations in the boundary layer and upper troposphere above the Indian Ocean reasonably well. Using monthly resolved emissions, the main oceanic source regions for the stratosphere include the Arabian Sea and Bay of Bengal in boreal summer and the tropical west Pacific Ocean in boreal winter. The main stratospheric injection in boreal summer occurs over the southern tip of India associated with the high local oceanic sources and strong convection of the summer monsoon. In boreal winter more bromoform is entrained over the west Pacific than over the Indian Ocean. The annually averaged stratospheric injection of bromoform is in the same range whether using monthly averaged or annually averaged emissions in our Lagrangian calculations. However, monthly averaged emissions result in the highest mixing ratios within the Asian monsoon anticyclone in boreal summer and above the central Indian Ocean in boreal winter, while annually averaged emissions display a maximum above the west Indian Ocean in boreal spring. In the Asian summer monsoon anticyclone bromoform atmospheric mixing ratios vary by up to 50% between using monthly averaged and annually averaged oceanic emissions. Our results underline that the seasonal and regional stratospheric bromine injection from the tropical Indian Ocean and west Pacific critically depend on the seasonality and spatial distribution of the VSLS emissions.

Description: The calcareous tubeworm Spirorbis spirorbis is a widespread serpulid species in the Baltic Sea, where it commonly grows as an epibiont on brown macroalgae (genus Fucus). It lives within a Mg-calcite shell and could be affected by ocean acidification and temperature rise induced by the predicted future atmospheric CO2 increase. However, Spirorbis tubes grow in a chemically modified boundary layer around the algae, which may mitigate acidification. In order to investigate how increasing temperature and rising pCO2 may influence S. spirorbis shell growth we carried out four seasonal experiments in the Kiel Outdoor Benthocosms at elevated pCO2 and temperature conditions. Compared to laboratory batch culture experiments the benthocosm approach provides a better representation of natural conditions for physical and biological ecosystem parameters, including seasonal variations. We find that growth rates of S. spirorbis are significantly controlled by ontogenetic and seasonal effects. The length of the newly grown tube is inversely related to the initial diameter of the shell. Our study showed no significant difference of the growth rates between ambient atmospheric and elevated (1100 ppm) pCO2 conditions. No influence of daily average CaCO3 saturation state on the growth rates of S. spirorbis was observed. We found, however, net growth of the shells even in temporarily undersaturated bulk solutions, under conditions that concurrently favoured selective shell surface dissolution. The results suggest an overall resistance of S. spirorbis growth to acidification levels predicted for the year 2100 in the Baltic Sea. In contrast, S. spirorbis did not survive at mean seasonal temperatures exceeding 24 °C during the summer experiments. In the autumn experiments at ambient pCO2, the growth rates of juvenile S. spirorbis were higher under elevated temperature conditions. The results reveal that S. spirorbis may prefer moderately warmer conditions during their early life stages but will suffer from an excessive temperature increase and from increasing shell corrosion as a consequence of progressing ocean acidification.

Description: n the Paris Agreement in 2015 countries agreed on holding global mean surface air warming to "well below 2 degrees C above pre-industrial" levels, but the emission reduction pledges under that agreement are not ambitious enough to meet this target. Therefore, the question arises of whether restoring global warming to this target after exceeding it by artificially removing CO2 from the atmosphere is possible. One important aspect is the reversibility of ocean heat uptake and associated sea level rise, which have very long (centennial to millennial) response timescales. In this study the response of sea level rise due to thermal expansion to a 1% yearly increase of atmospheric CO2 up to a quadrupling of the pre-industrial concentration followed by a 1% yearly decline back to the pre-industrial CO2 concentration is examined using the University of Victoria Earth System Climate Model (UVic ESCM). We find that global mean thermosteric sea level (GMTSL) continues to rise for several decades after atmospheric CO2 starts to decline and does not return to pre-industrial levels for over 1000 years after atmospheric CO2 is restored to the pre-industrial concentration. This finding is independent of the strength of vertical sub-grid-scale ocean mixing implemented in the model. Furthermore, GMTSL rises faster than it declines in response to a symmetric rise and decline in atmospheric CO2 concentration partly because the deep ocean continues to warm for centuries after atmospheric CO2 returns to the pre-industrial concentration. Both GMTSL rise and decline rates increase with increasing vertical ocean mixing. Exceptions from this behaviour arise if the overturning circulations in the North Atlantic and Southern Ocean intensify beyond pre-industrial levels in model versions with lower vertical mixing, which leads to rapid cooling of the deep ocean.

Description: The climate research community uses atmospheric reanalysis data sets to understand a wide range of processes and variability in the atmosphere, yet different reanalyses may give very different results for the same diagnostics. The Stratosphere–troposphere Processes And their Role in Climate (SPARC) Reanalysis Intercomparison Project (S-RIP) is a coordinated activity to compare reanalysis data sets using a variety of key diagnostics. The objectives of this project are to identify differences among reanalyses and understand their underlying causes, to provide guidance on appropriate usage of various reanalysis products in scientific studies, particularly those of relevance to SPARC, and to contribute to future improvements in the reanalysis products by establishing collaborative links between reanalysis centres and data users. The project focuses predominantly on differences among reanalyses, although studies that include operational analyses and studies comparing reanalyses with observations are also included when appropriate. The emphasis is on diagnostics of the upper troposphere, stratosphere, and lower mesosphere. This paper summarizes the motivation and goals of the S-RIP activity and extensively reviews key technical aspects of the reanalysis data sets that are the focus of this activity. The special issue "The SPARC Reanalysis Intercomparison Project (S-RIP)" in this journal serves to collect research with relevance to the S-RIP in preparation for the publication of the planned two (interim and full) S-RIP reports.

Description: A high-resolution dynamical model (Nucleus for European Modelling of the Ocean, Mediterranean configuration – NEMO-MED12) was used to give the first simulation of the distribution of radiocarbon (14C) across the whole Mediterranean Sea. The simulation provides a descriptive overview of both the natural pre-bomb 14C and the entire anthropogenic radiocarbon transient generated by the atmospheric bomb tests performed in the 1950s and early 1960s. The simulation was run until 2011 to give the post-bomb distribution. The results are compared to available in situ measurements and proxy-based reconstructions. The radiocarbon simulation allows an additional and independent test of the dynamical model, NEMO-MED12, and its performance to produce the thermohaline circulation and deep-water ventilation. The model produces a generally realistic distribution of radiocarbon when compared with available in situ data. The results demonstrate the major influence of the flux of Atlantic water through the Strait of Gibraltar on the inter-basin natural radiocarbon distribution and characterize the ventilation of intermediate and deep water especially through the propagation of the anthropogenic radiocarbon signal. We explored the impact of the interannual variability on the radiocarbon distribution during the Eastern Mediterranean Transient (EMT) event. It reveals a significant increase in 14C concentration (by more than 60 ‰) in the Aegean deep water and at an intermediate level (value up to 10 ‰) in the western basin. The model shows that the EMT makes a major contribution to the accumulation of radiocarbon in the eastern Mediterranean deep waters.

Description: Vertebrates have evolved a sophisticated adaptive immune system that relies on an almost infinite diversity of antigen receptors that are clonally expressed by specialized immune cells that roam the circulatory system. These immune cells provide vertebrates with extraordinary antigen-specific immune capacity and memory, while minimizing self-reactivity. Plants, however, lack specialized mobile immune cells. Instead, every plant cell is thought to be capable of launching an effective immune response. So how do plants achieve specific, self-tolerant immunity and establish immune memory? Recent developments point towards a multilayered plant innate immune system comprised of self-surveillance, systemic signalling and chromosomal changes that together establish effective immunity.

Description: To survey hepatitis B virus (HBV) integration in liver cancer genomes, we conducted massively parallel sequencing of 81 HBV-positive and 7 HBV-negative hepatocellular carcinomas (HCCs) and adjacent normal tissues. We found that HBV integration is observed more frequently in the tumors (86.4%) than in adjacent liver tissues (30.7%). Copy-number variations (CNVs) were significantly increased at HBV breakpoint locations where chromosomal instability was likely induced. Approximately 40% of HBV breakpoints within the HBV genome were located within a 1,800-bp region where the viral enhancer, X gene and core gene are located. We also identified recurrent HBV integration events (in ≥4 HCCs) that were validated by RNA sequencing (RNA-seq) and Sanger sequencing at the known and putative cancer-related TERT, MLL4 and CCNE1 genes, which showed upregulated gene expression in tumor versus normal tissue. We also report evidence that suggests that the number of HBV integrations is associated with patient survival.

Description: WHILE ammonites and all other ectocochleate cephalopods became extinct, nautiloids survived relatively unchanged from the Ordovician, suggesting that they are unusually well adapted to their niche. Here we obtain high-resolution tracks of Nautilus positions and depths, combined with telemetered jet pressures, which clarify both its lifestyle and economics. Nautilus is more active in nature than in captivity1, but its energy costs are lower than projected2,3. Viewing Nautilus as 'vertic', rather than benthic, resolves this contradiction. Records show that the cost of transport is the same in any direction within a vertical plane. Living on a reef face swept by a lateral current means that vertical movements4,5 sample large areas for chemical trails. A detected trail can be followed upcurrent in the slow-moving boundary layer, but no effort is wasted on horizontal movement without good prospects for food; long-range movements are downcurrent and made by drifting. Once fed, a Nautilus can reduce its energy costs by moving to deeper, cooler waters, where a single meal can last for months.

Description: Cold-water coral (CWC) reefs are heterogeneous ecosystems comprising numerous microhabitats. A typical European CWC reef provides various biogenic microhabitats (within, on and surrounding colonies of coral species such as Lophelia pertusa, Paragorgia arborea and Primnoa resedaeformis, or formed by their remains after death). These microhabitats may be surrounded and intermixed with non-biogenic microhabitats (soft sediment, hard ground, gravel/pebbles, steep walls). To date, studies of distribution of sessile fauna across CWC reefs have been more numerous than those investigating mobile fauna distribution. In this study we quantified shrimp densities associated with key CWC microhabitat categories at the Røst Reef, Norway, by analysing image data collected by towed video sled in June 2007. We also investigated shrimp distribution patterns on the local scale (〈40 cm) and how these may vary with microhabitat. Shrimp abundances at the Røst Reef were on average an order of magnitude greater in biogenic reef microhabitats than in non-biogenic microhabitats. Greatest shrimp densities were observed in association with live Paragorgia arborea microhabitat (43 shrimp m−2, SD = 35.5), live Primnoa resedaeformis microhabitat (41.6 shrimp m−2, SD = 26.1) and live Lophelia pertusa microhabitat (24.4 shrimp m−2, SD = 18.6). In non-biogenic microhabitat, shrimp densities were 〈2 shrimp m−2. CWC reef microhabitats appear to support greater shrimp densities than the surrounding non-biogenic microhabitats at the Røst Reef, at least at the time of survey.

Description: Extensive biogeochemical transformation of organic matter takes place in the shallow continental shelf seas of Siberia. This, in combination with brine production from sea-ice formation, results in cold bottom waters with relatively high salinity and nutrient concentrations, as well as low oxygen and pH levels. Data from the SWERUS-C3 expedition with icebreaker Oden, from July to September 2014, show the distribution of such nutrient-rich, cold bottom waters along the continental margin from about 140 to 180° E. The water with maximum nutrient concentration, classically named the upper halocline, is absent over the Lomonosov Ridge at 140° E, while it appears in the Makarov Basin at 150° E and intensifies further eastwards. At the intercept between the Mendeleev Ridge and the East Siberian continental shelf slope, the nutrient maximum is still intense, but distributed across a larger depth interval. The nutrient-rich water is found here at salinities of up to ∼ 34.5, i.e. in the water classically named lower halocline. East of 170° E transient tracers show significantly less ventilated waters below about 150 m water depth. This likely results from a local isolation of waters over the Chukchi Abyssal Plain as the boundary current from the west is steered away from this area by the bathymetry of the Mendeleev Ridge. The water with salinities of ∼ 34.5 has high nutrients and low oxygen concentrations as well as low pH, typically indicating decay of organic matter. A deficit in nitrate relative to phosphate suggests that this process partly occurs under hypoxia. We conclude that the high nutrient water with salinity ∼ 34.5 are formed on the shelf slope in the Mendeleev Ridge region from interior basin water that is trapped for enough time to attain its signature through interaction with the sediment.

Description: Submarine canyons are known as one of the seafloor morphological features where living cold-water coral (CWC) communities develop in the Mediterranean Sea. We investigated the CWC community of the two westernmost submarine canyons of the Gulf of Lions canyon system: the Cap de Creus Canyon (CCC) and Lacaze-Duthiers Canyon (LDC). Coral associations have been studied through video material recorded by means of a manned submersible and a remotely operated vehicle. Video transects have been conducted and analyzed in order to obtain information on (1) coral bathymetric distribution and density patterns, (2) size structure of coral populations, and (3) coral colony position with respect to the substrate. Madrepora oculata was the most abundant CWC in both canyons, while Lophelia pertusa and Dendrophyllia cornigera mostly occurred as isolated colonies or in small patches. An important exception was detected in a vertical cliff in LDC where a large L. pertusa framework was documented. This is the first record of such an extended L. pertusa framework in the Mediterranean Sea. In both canyons coral populations were dominated by medium and large colonies, but the frequent presence of small-sized colonies also indicate active recruitment. The predominant coral orientation (90° and 135°) is probably driven by the current regime as well as by the sediment load transported by the current flows. In general, no clear differences were observed in the abundance and in the size structure of the CWC populations between CCC and LDC, despite large differences in particulate matter between canyons.

Description: As a major source for atmospheric CO2, the Peruvian upwelling region exhibits strong variability in surface fCO2 on short spatial and temporal scales. Understanding the physical processes driving the strong variability is of fundamental importance for constraining the effect of marine emissions from upwelling regions on the global CO2 budget. In this study, a frontal decay on length scales of 𝒪(10 km) was observed off the Peruvian coast following a pronounced decrease in down-frontal (equatorward) wind speed with a time lag of 9 h. Simultaneously, the sea-to-air flux of CO2 on the inshore (cold) side of the front dropped from up to 80 to 10 mmol m−2 day−1, while the offshore (warm) side of the front was constantly outgassing at a rate of 10–20 mmol m−2 day−1. Based on repeated ship transects the decay of the front was observed to occur in two phases. The first phase was characterized by a development of coherent surface temperature anomalies which gained in amplitude over 6–9 h. The second phase was characterized by a disappearance of the surface temperature front within 6 h. Submesoscale mixed-layer instabilities were present but seem too slow to completely remove the temperature gradient in this short time period. Dynamics such as a pressure-driven gravity current appear to be a likely mechanism behind the evolution of the front.

Description: Climate warming at the end of the last glaciation caused ice caps on Icelandic volcanoes to retreat. Removal of surface ice load is thought to have decreased pressures in the underlying mantle, triggering decompression melting, enhanced magma generation and increased volcanic activity1–3. Present-day climate change could have the same effect, although there may be a time lag of hundreds of years between magma generation and eruption4,5. However, in addition to increased magma generation, pressure changes associated with ice retreat should also alter the capacity for storing magma within the crust. Here we use a numerical model to evaluate the effect of the current decrease in ice load on magma storage in the crust at the Kverkfjöll volcanic system, located partially beneath Iceland’s largest ice cap. We compare the model results with radar and global positioning system measurements of surface displacement and changes in crustal stress between 2007 and 2008, during the intrusion of a deep dyke at Upptyppingar. We find that although the main component of stress recorded during dyke intrusion relates to plate extension, another component of stress is consistent with the stress field caused by the retreating ice cap. We conclude that the retreating ice cap led to enhanced capture of magma within the crust. We suggest that ice-cap retreat can promote magma storage, rather than eruption, at least in the short term.

Description: A major role in regulation of global methane fluxes has been attributed to the process of anaerobic oxidation of methane (AOM), which is performed by consortia of methanotrophic archaea and sulfate reducing bacteria. An important question remains how these energy limited, slow growing Microorganisms with generation times of 3–7 months respond to rapid natural variations in methane fluxes at cold seeps. We used an experimental flow-through column system filled with cold seep sediments naturally enriched in methanotrophic communities, to test their responses to short-term variations in methane and sulfate fluxes. At stable methane and sulfate concentrations of ∼2 mM and 28 mM, respectively, we measured constant rates of AOM and sulfate reduction (SR) for up to 160 days of incubation. When percolated with methane-free medium, the anaerobic methanotrophs ceased to produce sulfide. After a starvation phase of 40 days, the addition of methane restored former AOM and SR rates immediately. At methane concentrations between 0–2.3 mM we measured a linear correlation between methane availability, AOM and SR. At constant fluid flow velocities of 30 m yr−1, ca. 50% of the methane was consumed by the anaerobic methanotrophic (ANME) population at all concentrations tested. Reducing the sulfate concentration from 28 to 1 mM, a decrease in AOM and SR by 50% was observed, and 45% of the methane was consumed. Hence, the marine anaerobic methanotrophs(ANME) are capable of oxidizing substantial amounts of methane over a wide and variable range of fluxes of the reaction educts.

Description: The island of South Georgia is situated in the iron (Fe) depleted Antarctic Circumpolar Current of the Southern Ocean. Iron emanating from its shelf system fuels large phytoplankton blooms downstream of the island, but the actual supply mechanisms are unclear. To address this we present the first inventory of Fe, manganese (Mn) and aluminium (Al) in shelf sediments, pore waters and the water column in the vicinity of South Georgia, alongside data on zooplankton-mediated Fe cycling processes. The seafloor sediments were the main particulate Fe source to shelf bottom waters as indicated by Fe / Mn and Fe / Al ratios for shelf sediments and suspended particles in the water column. Less than 1 % of the total particulate Fe pool was leachable surface adsorbed (labile) Fe, and therefore potentially available to organisms. Pore waters formed the primary dissolved Fe (DFe) source to shelf bottom waters supplying 0.1–4 μmol DFe m−2 d−1. However, only 0.41 ± 0.26 μmol DFe m−2 d−1 was transferred to the surface mixed layer by vertical diffusive and advective mixing. Other trace metal sources to surface waters included glacial flour released by melting glaciers and zooplankton excretion processes. On average 6.5 ± 8.2 μmol m−2 d−1 of labile particulate Fe was supplied to the surface mixed layer via krill faecal pellets, with further DFe released by krill at around 1.1 ± 2.2 μmol m−2 d−1. The faecal pellets released by krill constituted of seafloor derived lithogenic material and settled algae debris, in addition to freshly ingested suspended phytoplankton specimen. The phytoplankton Fe requirement in the blooms ca. 1250 km downstream the island of South Georgia was 0.33 ± 0.11 μmol m−2 d−1, with the DFe supply by horizontal/vertical mixing, deep winter mixing and via aeolian dust estimated as ~ 0.12 μmol m−2 d−1. We suggest that additionally required DFe was provided through recycling of biogenically stored Fe following luxury Fe uptake by phytoplankton on the Fe rich shelf. This process would allow Fe to be retained in the surface mixed layer of waters downstream of South Georgia through continuous recycling and biological uptake, and facilitate the large scale blooms.

Description: In estuarine coastal systems such as the Baltic Sea, mussels suffer from low salinity which limits their distribution. Anthropogenic climate change is expected to cause further desalination which will lead to local extinctions of mussels in the low saline areas. It is commonly accepted that mussel distribution is limited by osmotic stress. However, along the salinity gradient environmental conditions for biomineralization are successively becoming more adverse as a result of reduced [Ca2+] and dissolved inorganic carbon (CT) availability. In larvae, calcification is an essential process starting during early development with formation of the prodissoconch I (PD I) shell which is completed under optimal conditions within 2 days. Experimental manipulations of seawater [Ca2+] start to impair PD I formation in Mytilus larvae at concentrations below 3 mM, which corresponds to conditions present in the Baltic at salinities below 8 g kg-1. In addition, lowering dissolved inorganic carbon to critical concentrations (〈 1 mM) similarly affected PD I size which was well correlated with calculated ΩAragonite and [Ca2+][HCO3-]/[H+] in all treatments. Comparing results for larvae from the western Baltic with a population from the central Baltic revealed significantly higher tolerance of PD I formation to lowered [Ca2+] and [Ca2+][HCO3-]/[H+] in the low saline adapted population. This may result from genetic adaptation to the more adverse environmental conditions prevailing in the low saline areas of the Baltic. The combined effects of lowered [Ca2+] and adverse carbonate chemistry represent major limiting factors for bivalve calcification and can thereby contribute to distribution limits of mussels in the Baltic Sea.

Description: The elemental stoichiometry of sea water and particulate organic matter is remarkably similar. This observation led Redfield to hypothesize that the oceanic ratio of nitrate to phosphate is controlled by the remineralization of phytoplankton biomass1. The Redfield ratio is used universally to quantitatively link the marine nitrogen and phosphorus cycles in numerous biogeochemical applications2,3,4. Yet, empirical and theoretical studies show that the ratio of nitrogen to phosphorus in phytoplankton varies greatly with taxa5,6 and growth conditions7,8,9. Here we present a dynamic five-box ecosystem model showing that non-Redfield utilization of dissolved nitrogen and phosphorus by non-nitrogen-fixing phytoplankton controls the magnitude and distribution of nitrogen fixation. In our simulations, systems dominated by rapidly growing phytoplankton with low nitrogen to phosphorus uptake ratios reduce the phosphorus available for nitrogen fixation. In contrast, in systems dominated by slow-growing phytoplankton with high nitrogen to phosphorus uptake ratios nitrogen deficits are enhanced, and nitrogen fixation is promoted. We show that estimates of nitrogen fixation are up to fourfold too high when non-Redfield uptake stoichiometries are ignored. We suggest that the relative abundance of fast- and slow-growing phytoplankton controls the amount of new nitrogen added to the ocean.

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

Description: The response of the marine carbon cycle to changes in atmospheric CO2 concentrations will be determined, in part, by the relative response of calcifying and non-calcifying organisms to global change. Planktonic foraminifera are responsible for a quarter or more of global carbonate production, therefore understanding the sensitivity of calcification in these organisms to environmental change is critical. Despite this, there remains little consensus as to whether, or to what extent, chemical and physical factors affect foraminiferal calcification. To address this, we directly test the effect of multiple controls on calcification in culture experiments and core-top measurements of Globigerinoides ruber. We find that two factors, body size and the carbonate system, strongly influence calcification intensity in life, but that exposure to corrosive bottom waters can overprint this signal post mortem. Using a simple model for the addition of calcite through ontogeny, we show that variable body size between and within datasets could complicate studies that examine environmental controls on foraminiferal shell weight. In addition, we suggest that size could ultimately play a role in determining whether calcification will increase or decrease with acidification. Our models highlight that knowledge of the specific morphological and physiological mechanisms driving ontogenetic change in calcification in different species will be critical in predicting the response of foraminiferal calcification to future change in atmospheric pCO2.

Description: The distribution of dissolved aluminium (dAl) in the water column of the North Atlantic and Labrador Sea was studied along GEOTRACES section GA01 to unravel the sources and sinks of this element. Surface water dAl concentrations were low (median of 2.5 nM) due to low aerosol deposition and removal by phytoplankton. However, surface water dAl concentrations were enhanced on the Iberian and Greenland shelves (up to 30.9 nM) due to continental inputs (rivers, glacial flour and ice melt). A negative correlation was observed between dAl in surface waters and primary production, phytoplankton community structure and biogenic opal production. The abundance of diatoms exerted a significant (p 〈 0.01) control on the surface particulate Al (pAl) to dAl ratios by decreasing dAl levels and increasing pAl levels. Dissolved Al concentrations generally increased with depth and correlated strongly with silicate (R2 〉 0.76) west of the Iberian Basin, suggesting net release of dAl at depth during remineralization of sinking biogenic opal containing particles. Enrichment of dAl at near-bottom depths was observed due to resuspension of sediments near the sediment-water interface. The highest dAl (up to 38.7 nM) concentrations were observed in Mediterranean Overflow Waters which act as a major source of dAl to mid depth waters of the eastern North Atlantic. This study clearly shows that the vertical and lateral distribution of dAl in the North Atlantic differs when compared to other regions of the North Atlantic and global ocean due to the large spatial differences both in the main source of Al, atmospheric deposition, and the main sink for Al, particle scavenging, between different oceanic regions.

Description: Autonomous underwater vehicles (AUVs) offer unique possibilities for exploring the deep seafloor in high resolution over large areas. We highlight the results from AUV-based multibeam echosounder (MBES) bathymetry / backscatter and digital optical imagery from the DISCOL area acquired during research cruise SO242 in 2015. AUV bathymetry reveals a morphologically complex seafloor with rough terrain in seamount areas and low-relief variations in sedimentary abyssal plains which are covered in Mn-nodules. Backscatter provides valuable information about the seafloor type and particularly about the influence of Mn-nodules on the response of the transmitted acoustic signal. Primarily, Mn-nodule abundances were determined by means of automated nodule detection on AUV seafloor imagery and nodule metrics such as nodules m−2 were calculated automatically for each image allowing further spatial analysis within GIS in conjunction with the acoustic data. AUV-based backscatter was clustered using both raw data and corrected backscatter mosaics. In total, two unsupervised methods and one machine learning approach were utilized for backscatter classification and Mn-nodule predictive mapping. Bayesian statistical analysis was applied to the raw backscatter values resulting in six acoustic classes. In addition, Iterative Self-Organizing Data Analysis (ISODATA) clustering was applied to the backscatter mosaic and its statistics (mean, mode, 10th, and 90th quantiles) suggesting an optimum of six clusters as well. Part of the nodule metrics data was combined with bathymetry, bathymetric derivatives and backscatter statistics for predictive mapping of the Mn-nodule density using a Random Forest classifier. Results indicate that acoustic classes, predictions from Random Forest model and image-based nodule metrics show very similar spatial distribution patterns with acoustic classes hence capturing most of the fine-scale Mn-nodule variability. Backscatter classes reflect areas with homogeneous nodule density. A strong influence of mean backscatter, fine scale BPI and concavity of the bathymetry on nodule prediction is seen. These observations imply that nodule densities are generally affected by local micro-bathymetry in a way that is not yet fully understood. However, it can be concluded that the spatial occurrence of Mn-covered areas can be sufficiently analysed by means of acoustic classification and multivariate predictive mapping allowing to determine the spatial nodule density in a much more robust way than previously possible.

Description: Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere – the "global carbon budget" – is important to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe data sets and methodology to quantify the five major components of the global carbon budget and their uncertainties. CO2 emissions from fossil fuels and industry (EFF) are based on energy statistics and cement production data, respectively, while emissions from land-use change (ELUC), mainly deforestation, are based on land-cover change data and bookkeeping models. The global atmospheric CO2 concentration is measured directly and its rate of growth (GATM) is computed from the annual changes in concentration. The ocean CO2 sink (SOCEAN) and terrestrial CO2 sink (SLAND) are estimated with global process models constrained by observations. The resulting carbon budget imbalance (BIM), the difference between the estimated total emissions and the estimated changes in the atmosphere, ocean, and terrestrial biosphere, is a measure of imperfect data and understanding of the contemporary carbon cycle. All uncertainties are reported as ±1σ. For the last decade available (2007–2016), EFF was 9.4 ± 0.5 GtC yr−1, ELUC 1.3 ± 0.7 GtC yr−1, GATM 4.7 ± 0.1 GtC yr−1, SOCEAN 2.4 ± 0.5 GtC yr−1, and SLAND 3.0 ± 0.8 GtC yr−1, with a budget imbalance BIM of 0.6 GtC yr−1 indicating overestimated emissions and/or underestimated sinks. For year 2016 alone, the growth in EFF was approximately zero and emissions remained at 9.9 ± 0.5 GtC yr−1. Also for 2016, ELUC was 1.3 ± 0.7 GtC yr−1, GATM was 6.1 ± 0.2 GtC yr−1, SOCEAN was 2.6 ± 0.5 GtC yr−1, and SLAND was 2.7 ± 1.0 GtC yr−1, with a small BIM of −0.3 GtC. GATM continued to be higher in 2016 compared to the past decade (2007–2016), reflecting in part the high fossil emissions and the small SLAND consistent with El Niño conditions. The global atmospheric CO2 concentration reached 402.8 ± 0.1 ppm averaged over 2016. For 2017, preliminary data for the first 6–9 months indicate a renewed growth in EFF of +2.0 % (range of 0.8 to 3.0 %) based on national emissions projections for China, USA, and India, and projections of gross domestic product (GDP) corrected for recent changes in the carbon intensity of the economy for the rest of the world. This living data update documents changes in the methods and data sets used in this new global carbon budget compared with previous publications of this data set (Le Quéré et al., 2016, 2015b, a, 2014, 2013). All results presented here can be downloaded from https://doi.org/10.18160/GCP-2017 (GCP, 2017).

Description: A long-term mean flow field for the subpolar North Atlantic region with a horizontal resolution of approximately 25km is created by gridding Argo-derived velocity vectors using two different topography-following interpolation schemes. The 10-day float displacements in the typical drift depths of 1000 to 1500m represent the flow in the Labrador Sea Water density range. Both mapping algorithms separate the flow field into potential vorticity (PV) conserving, i.e., topography-following contribution and a deviating part, which we define as the eddy contribution. To verify the significance of the separation, we compare the mean flow and the eddy kinetic energy (EKE), derived from both mapping algorithms, with those obtained from multiyear mooring observations. The PV-conserving mean flow is characterized by stable boundary currents along all major topographic features including shelf breaks and basin-interior topographic ridges such as the Reykjanes Ridge or the Rockall Plateau. Mid-basin northward advection pathways from the northeastern Labrador Sea into the Irminger Sea and from the Mid-Atlantic Ridge region into the Iceland Basin are well-resolved. An eastward flow is present across the southern boundary of the subpolar gyre near 52°N, the latitude of the Charlie Gibbs Fracture Zone (CGFZ). The mid-depth EKE field resembles most of the satellite-derived surface EKE field. However, noticeable differences exist along the northward advection pathways in the Irminger Sea and the Iceland Basin, where the deep EKE exceeds the surface EKE field. Further, the ratio between mean flow and the square root of the EKE, the Peclet number, reveals distinct advection-dominated regions as well as basin-interior regimes in which mixing is prevailing.

Description: Past ocean temperatures and salinities can be approximated from combined stable oxygen isotopes (δ18O) and Mg ∕ Ca measurements in fossil foraminiferal tests with varying success. To further refine this approach, we collected living planktic foraminifers by net sampling and pumping of sea surface water from the Caribbean Sea, the eastern Gulf of Mexico and the Florida Straits. Analyses of δ18O and Mg ∕ Ca in eight living planktic species (Globigerinoides sacculifer, Orbulina universa, Neogloboquadrina dutertrei, Pulleniatina obliquiloculata, Globorotalia menardii, Globorotalia ungulata, Globorotalia truncatulinoides and Globorotalia tumida) were compared to measured in situ properties of the ambient seawater (temperature, salinity and δ18Oseawater) and fossil tests of underlying surface sediments. “Vital effects” such as symbiont activity and test growth cause δ18O disequilibria with respect to the ambient seawater and a large scatter in foraminiferal Mg ∕ Ca. Overall, ocean temperature is the most prominent environmental influence on δ18Ocalcite and Mg ∕ Ca. Enrichment of the heavier 18O isotope in living specimens below the mixed layer and in fossil tests is clearly related to lowered in situ temperatures and gametogenic calcification. Mg ∕ Ca-based temperature estimates of G. sacculifer indicate seasonal maximum accumulation rates on the seafloor in early spring (March) at Caribbean stations and later in the year (May) in the Florida Straits, related to the respective mixed layer temperatures of ∼26 ∘C. Notably, G. sacculifer reveals a weak positive linear relationship between foraminiferal derived δ18Oseawater estimates and both measured in situ δ18Oseawater and salinity. Our results affirm the applicability of existing δ18O and Mg ∕ Ca calibrations for the reconstruction of past ocean temperatures and δ18Oseawater reflecting salinity due to the convincing accordance of proxy data in both living and fossil foraminifers, and in situ environmental parameters. Large vital effects and seasonally varying proxy signals, however, need to be taken into account.

Description: We propose a satellite mission that uses a near-nadir Ka-band Doppler radar to measure surface currents, ice drift and ocean waves at spatial scales of 40 km and more, with snapshots at least every day for latitudes 75 to 82°, and every few days for other latitudes. The use of incidence angles of 6 and 12° allows for measurement of the directional wave spectrum, which yields accurate corrections of the wave-induced bias in the current measurements. The instrument's design, an algorithm for current vector retrieval and the expected mission performance are presented here. The instrument proposed can reveal features of tropical ocean and marginal ice zone (MIZ) dynamics that are inaccessible to other measurement systems, and providing global monitoring of the ocean mesoscale that surpasses the capability of today's nadir altimeters. Measuring ocean wave properties has many applications, including examining wave–current interactions, air–sea fluxes, the transport and convergence of marine plastic debris and assessment of marine and coastal hazards.

Description: Gel particles, such as transparent exopolymer particles (TEP) and Coomassie stainable particles (CSP), are important organic components in the sea surface microlayer (SML). Here, we present results on the effect of different wind speeds on the accumulation and size distribution of TEP and CSP during a wind wave channel experiment in the Aeolotron. Total areas of TEP (TEPSML) and CSP (CSPSML) in the surface microlayer were exponentially related to wind speed. At wind speeds  〈  6 m s−1, accumulation of TEPSML and CSPSML occurred, decreasing at wind speeds of  〉  8 m s−1. Wind speeds  〉  8 m s−1 also significantly altered the size distribution of TEPSML in the 2–16 µm size range towards smaller sizes. The response of the CSPSML size distribution to wind speed varied through time depending on the biogenic source of gels. Wind speeds  〉  8 m s−1 decreased the slope of CSPSML size distribution significantly in the absence of autotrophic growth. For the slopes of TEP and CSP size distribution in the bulk water, no significant difference was observed between high and low wind speeds. Changes in spectral slopes between high and low wind speed were higher for TEPSML than for CSPSML, indicating that the impact of wind speed on size distribution of gel particles in the SML may be more pronounced for TEP than for CSP, and that CSPSML are less prone to aggregation during the low wind speeds. Addition of an E. huxleyi culture resulted in a higher contribution of submicron gels (0.4–1 µm) in the SML at higher wind speed ( 〉  6 m s−1), indicating that phytoplankton growth may potentially support the emission of submicron gels with sea spray aerosol.

Description: In the ocean, sinking of particulate organic matter (POM) drives carbon export from the euphotic zone and supplies nutrition to mesopelagic communities, the feeding and degradation activities of which in turn lead to export flux attenuation. Oxygen minimum zones (OMZs) with suboxic water layers (〈 5 μmol O2 kg−1) show a lower carbon flux attenuation compared to well oxygenated waters (〉 100 μmol O2 kg−1), supposedly due to reduced heterotrophic activity. This study focuses on sinking particle fluxes through hypoxic mesopelagic waters (〈 60 % μmol O2 kg−1); these represent ~ 100-times more ocean volume globally compared to suboxic waters, but have less been studied. Particle export fluxes and attenuation coefficients were determined in the Eastern Tropical North Atlantic (ETNA) using two surface tethered drifting sediment trap arrays with 7 trapping depths located between 100 and 600 m. Data on particulate matter fluxes were fitted to the normalized power function Fz = F100 (z/100)−b, with F100 being the flux at a depth (z) of 100 m and b being the attenuation coefficient. Higher b-values suggest stronger flux attenuation and are influenced by factors such as faster degradation at higher temperatures. In this study, b-values of organic carbon fluxes varied between 0.74 and 0.80 and were in the intermediate range of previous reports, but lower than expected from seawater temperatures within the upper 500 m. During this study, highest b-values were determined for fluxes of particulate hydrolysable amino acids (PHAA), followed by particulate organic phosphorus (POP), nitrogen (PN), carbon (POC), chlorophyll a, and transparent exopolymer particles (TEP), pointing to a sequential degradation of organic matter components during sinking. Our study suggests that in addition to oxygen concentration, organic matter composition co-determines transfer efficiency through the mesopelagial. The magnitude of future carbon export fluxes may therefore also depend on how organic matter quality in the surface ocean changes under influence of warming, acidification, and enhanced stratification.

Description: Well water level changes associated with magmatic unrest can be interpreted as a result of pore pressure changes in the aquifer due to crustal deformation, and so could provide constraints on the subsurface processes causing this strain. We use finite element analysis to demonstrate the response of aquifers to volumetric strain induced by pressurized magma reservoirs. Two different aquifers are invoked – an unconsolidated pyroclastic deposit and a vesicular lava flow – and embedded in an impermeable crust, overlying a magma chamber. The time-dependent, fully coupled models simulate crustal deformation accompanying chamber pressurization and the resulting hydraulic head changes as well as flow through the porous aquifer, i.e. porous flow. The simulated strain leads to centimetres (pyroclastic aquifer) to metres (lava flow aquifer) of hydraulic head changes; both strain and hydraulic head change with time due to substantial porous flow in the hydrological system. Well level changes are particularly sensitive to chamber volume, shape and pressurization strength, followed by aquifer permeability and the phase of the pore fluid. The depths of chamber and aquifer, as well as the aquifer's Young's modulus also have significant influence on the hydraulic head signal. While source characteristics, the distance between chamber and aquifer and the elastic stratigraphy determine the strain field and its partitioning, flow and coupling parameters define how the aquifer responds to this strain and how signals change with time. We find that generic analytical models can fail to capture the complex pre-eruptive subsurface mechanics leading to strain-induced well level changes, due to aquifer pressure changes being sensitive to chamber shape and lithological heterogeneities. In addition, the presence of a pore fluid and its flow have a significant influence on the strain signal in the aquifer and are commonly neglected in analytical models. These findings highlight the need for numerical models for the interpretation of observed well level signals. However, simulated water table changes do indeed mirror volumetric strain, and wells are therefore a valuable addition to monitoring systems that could provide important insights into pre-eruptive dynamics.

Description: The onset of the Palaeocene/Eocene thermal maximum (about 55 Myr ago) was marked by global surface temperatures warming by 5–7 °C over approximately 30,000 yr (ref. 1), probably because of enhanced mantle outgassing2, 3 and the pulsed release of approx1,500 gigatonnes of methane carbon from decomposing gas-hydrate reservoirs4, 5, 6, 7. The aftermath of this rapid, intense and global warming event may be the best example in the geological record of the response of the Earth to high atmospheric carbon dioxide concentrations and high temperatures. This response has been suggested to include an intensified flux of organic carbon from the ocean surface to the deep ocean and its subsequent burial through biogeochemical feedback mechanisms8. Here we present firm evidence for this view from two ocean drilling cores, which record the largest accumulation rates of biogenic barium—indicative of export palaeoproductivity—at times of maximum global temperatures and peak excursion values of delta13C. The unusually rapid return of delta13C to values similar to those before the methane release7 and the apparent coupling of the accumulation rates of biogenic barium to temperature, suggests that the enhanced deposition of organic matter to the deep sea may have efficiently cooled this greenhouse climate by the rapid removal of excess carbon dioxide from the atmosphere.

Description: The meridional Ekman volume, heat, and salt transport across two trans-Atlantic sections near 14.5° N and 11° S were estimated using in situ observations, wind products, and model data. A meridional ageostrophic velocity was obtained as the difference between the directly measured total velocity and the geostrophic velocity derived from observations. Interpreting the section mean ageostrophy to be the result of an Ekman balance, the meridional Ekman transport of 6. 2 ± 2. 3 Sv northward at 14.5° N and 11. 7 ± 2. 1 Sv southward at 11° S is estimated. The integration uses the top of the pycnocline as an approximation for the Ekman depth, which is on average about 20 m deeper than the mixed layer depth. The Ekman transport estimated based on the velocity observations agrees well with the predictions from in situ wind stress data of 6. 7 ± 3. 5 Sv at 14.5° N and 13. 6 ± 3. 3 Sv at 11° S. The meridional Ekman heat and salt fluxes calculated from sea surface temperature and salinity data or from high-resolution temperature and salinity profile data differ only marginally. The errors in the Ekman heat and salt flux calculation were dominated by the uncertainty of the Ekman volume transport estimates.

Description: Halogenated very short-lived substances (VSLSs) are naturally produced in the ocean and emitted to the atmosphere. When transported to the stratosphere, these compounds can have a significant influence on the ozone layer and climate. During a research cruise on RV Sonne in the subtropical and tropical west Indian Ocean in July and August 2014, we measured the VSLSs, methyl iodide (CH3I) and for the first time bromoform (CHBr3) and dibromomethane (CH2Br2), in surface seawater and the marine atmosphere to derive their emission strengths. Using the Lagrangian particle dispersion model FLEXPART with ERA-Interim meteorological fields, we calculated the direct contribution of observed VSLS emissions to the stratospheric halogen burden during the Asian summer monsoon. Furthermore, we compare the in situ calculations with the interannual variability of transport from a larger area of the west Indian Ocean surface to the stratosphere for July 2000–2015. We found that the west Indian Ocean is a strong source for CHBr3 (910 pmol m−2 h−1), very strong source for CH2Br2 (930 pmol m−2 h−1), and an average source for CH3I (460 pmol m−2 h−1). The atmospheric transport from the tropical west Indian Ocean surface to the stratosphere experiences two main pathways. On very short timescales, especially relevant for the shortest-lived compound CH3I (3.5 days lifetime), convection above the Indian Ocean lifts oceanic air masses and VSLSs towards the tropopause. On a longer timescale, the Asian summer monsoon circulation transports oceanic VSLSs towards India and the Bay of Bengal, where they are lifted with the monsoon convection and reach stratospheric levels in the southeastern part of the Asian monsoon anticyclone. This transport pathway is more important for the longer-lived brominated compounds (17 and 150 days lifetime for CHBr3 and CH2Br2). The entrainment of CHBr3 and CH3I from the west Indian Ocean to the stratosphere during the Asian summer monsoon is lower than from previous cruises in the tropical west Pacific Ocean during boreal autumn and early winter but higher than from the tropical Atlantic during boreal summer. In contrast, the projected CH2Br2 entrainment was very high because of the high emissions during the west Indian Ocean cruise. The 16-year July time series shows highest interannual variability for the shortest-lived CH3I and lowest for the longest-lived CH2Br2. During this time period, a small increase in VSLS entrainment from the west Indian Ocean through the Asian monsoon to the stratosphere is found. Overall, this study confirms that the subtropical and tropical west Indian Ocean is an important source region of halogenated VSLSs, especially CH2Br2, to the troposphere and stratosphere during the Asian summer monsoon.

Description: Treatment of the underwater light field in ocean biogeochemical models has been attracting increasing interest, with some models moving towards more complex parameterisations. We conduct a simple sensitivity study of a typical, highly simplified parameterisation. In our study, we vary the phytoplankton light attenuation parameter over a range constrained by data during both pre-industrial equilibrated and future climate scenario RCP8.5. In equilibrium, lower light attenuation parameters (weaker self-shading) shift net primary production (NPP) towards the high latitudes, while higher values of light attenuation (stronger shelf-shading) shift NPP towards the low latitudes. Climate forcing magnifies this relationship through changes in the distribution of nutrients both within and between ocean regions. Where and how NPP responds to climate forcing can determine the magnitude and sign of global NPP trends in this high CO2 future scenario. Ocean oxygen is particularly sensitive to parameter choice. Under higher CO2 concentrations, two simulations establish a strong biogeochemical feedback between the Southern Ocean and low-latitude Pacific that highlights the potential for regional teleconnection. Our simulations serve as a reminder that shifts in fundamental properties (e.g. light attenuation by phytoplankton) over deep time have the potential to alter global biogeochemistry.

Description: Occupying about 14 % of the world's surface, the Southern Ocean plays a fundamental role in ocean and atmosphere circulation, carbon cycling and Antarctic ice-sheet dynamics. Unfortunately, high interannual variability and a dearth of instrumental observations before the 1950s limits our understanding of how marine–atmosphere–ice domains interact on multi-decadal timescales and the impact of anthropogenic forcing. Here we integrate climate-sensitive tree growth with ocean and atmospheric observations on southwest Pacific subantarctic islands that lie at the boundary of polar and subtropical climates (52–54° S). Our annually resolved temperature reconstruction captures regional change since the 1870s and demonstrates a significant increase in variability from the 1940s, a phenomenon predating the observational record. Climate reanalysis and modelling show a parallel change in tropical Pacific sea surface temperatures that generate an atmospheric Rossby wave train which propagates across a large part of the Southern Hemisphere during the austral spring and summer. Our results suggest that modern observed high interannual variability was established across the mid-twentieth century, and that the influence of contemporary equatorial Pacific temperatures may now be a permanent feature across the mid- to high latitudes.

Description: The outer western Crimean shelf of the Black Sea is a natural laboratory to investigate effects of stable oxic versus varying hypoxic conditions on seafloor biogeochemical processes and benthic community structure. Bottom-water oxygen concentrations ranged from normoxic (175 μmol O2 L−1) and hypoxic (〈 63 μmol O2 L−1) or even anoxic/sulfidic conditions within a few kilometers' distance. Variations in oxygen concentrations between 160 and 10 μmol L−1 even occurred within hours close to the chemocline at 134 m water depth. Total oxygen uptake, including diffusive as well as fauna-mediated oxygen consumption, decreased from 15 mmol m−2 d−1 on average in the oxic zone, to 7 mmol m−2 d−1 on average in the hypoxic zone, correlating with changes in macrobenthos composition. Benthic diffusive oxygen uptake rates, comprising respiration of microorganisms and small meiofauna, were similar in oxic and hypoxic zones (on average 4.5 mmol m−2 d−1), but declined to 1.3 mmol m−2 d−1 in bottom waters with oxygen concentrations below 20 μmol L−1. Measurements and modeling of porewater profiles indicated that reoxidation of reduced compounds played only a minor role in diffusive oxygen uptake under the different oxygen conditions, leaving the major fraction to aerobic degradation of organic carbon. Remineralization efficiency decreased from nearly 100 % in the oxic zone, to 50 % in the oxic–hypoxic zone, to 10 % in the hypoxic–anoxic zone. Overall, the faunal remineralization rate was more important, but also more influenced by fluctuating oxygen concentrations, than microbial and geochemical oxidation processes.