ALBERT

Description: An ensemble-based data assimilation framework for a coupled ocean–atmosphere model is applied to investigate the influence of assimilating different types of ocean observations on the ocean and atmosphere simulation. The data assimilation is performed with the parallel data assimilation framework (PDAF) for the climate model AWI-CM. Observations of the ocean, namely satellite sea-surface temperature (SST) and temperature and salinity profiles, are assimilated into the ocean component. The atmospheric state is only influenced by the model dynamics. Different assimilation scenarios were carried out with different combinations of observations to investigate to what extent the assimilation into the coupled model leads to a better estimation of the state of the ocean as well as the atmosphere. The influence of the data assimilation is assessed by comparing the ocean prediction with dependent and independent ocean observations. For the atmosphere, the assimilation result is compared with the ERA-Interim atmospheric reanalysis data. The ocean temperature and salinity are improved by all the assimilation scenarios in the coupled system. The assimilation leads to a response of the atmosphere throughout the troposphere and impacts the global atmospheric circulation. Globally the temperature and wind speed are improved in the atmosphere on average.

Description: To evaluate the present sea-ice changes in a longer-term perspective the knowledge of sea-ice variability on pre-industrial and geological time scales is essential. For the interpretation of proxy reconstructions it is necessary to understand the recent signals of different sea-ice proxies from various regions. We present 260 new sediment surface samples collected in the (sub-) Arctic Oceans that were analysed for specific sea-ice (IP25) and open-water phytoplankton biomarkers (brassicasterol, dinosterol, HBI III). This new biomarker dataset was combined with 615 previously published biomarker surface samples into a pan-Arctic database. The resulting pan-Arctic biomarker and sea-ice index (PIP25) database shows a spatial distribution correlating well with the diverse modern sea-ice concentrations. We find correlations of PBIP25, PDIP25 and PIIIIP25 with spring and autumn concentrations. Similar correlations with modern sea-ice concentrations are observed in Baffin Bay. However, the correlations of the PIP25 indices with modern sea-ice concentrations differ in Fram Strait from those of the (sub-) Arctic dataset, which is likely caused by region-specific differences in sea-ice variability, nutrient availability and other environmental conditions. The extended (sea-ice) biomarker database strengthens the validity of biomarker sea-ice reconstructions in different Arctic regions and shows how different sea-ice proxies combined may resolve specific seasonal sea-ice signals.

Description: 〈jats:title〉Abstract〈/jats:title〉〈jats:p〉 〈jats:list〉 〈jats:list-item〉〈jats:p〉Social–ecological systems (SES) exhibit complex cause‐and‐effect relationships. Capturing, interpreting, and responding to signals that indicate changes in ecosystems is key for sustainable management in SES. Breaks in this signal–response chain, when feedbacks are missing, will allow change to continue until a point when abrupt ecological surprises may occur.〈/jats:p〉〈/jats:list-item〉 〈jats:list-item〉〈jats:p〉In these situations, societies and local ecosystems can often become uncoupled. In this paper, we demonstrate how the red loop–green loop (RL–GL) concept can be used to uncover missing feedbacks and to better understand past social–ecological dynamics. Reinstating these feedbacks in order to recouple the SES may ultimately create more sustainable systems on local scales.〈/jats:p〉〈/jats:list-item〉 〈jats:list-item〉〈jats:p〉The RL–GL concept can uncover missing feedbacks through the characterization of SES dynamics along a spectrum of human resource dependence. Drawing on diverse qualitative and quantitative data sources, we classify SES dynamics throughout the history of Jamaican coral reefs along the RL–GL spectrum. We uncover missing feedbacks in red‐loop and red‐trap scenarios from around the year 600 until now. The Jamaican coral reef SES dynamics have moved between all four dynamic states described in the RL–GL concept: green loop, green trap, red loop and red trap.〈/jats:p〉〈/jats:list-item〉 〈jats:list-item〉〈jats:p〉We then propose mechanisms to guide the current unsustainable red traps back to more sustainable green loops, involving mechanisms of seafood trade and ecological monitoring. By gradually moving away from seafood exports, Jamaica may be able to return to green‐loop dynamics between the local society and their locally sourced seafood. We discuss the potential benefits and drawbacks of this proposed intervention and give indications of why an export ban may insure against future missing feedbacks and could prolong the sustainability of the Jamaican coral reef ecosystem.〈/jats:p〉〈/jats:list-item〉 〈jats:list-item〉〈jats:p〉Our approach demonstrates how the RL–GL approach can uncover missing feedbacks in a coral reef SES, a way the concept has not been used before. We advocate for how the RL–GL concept in a feedback setting can be used to synthesize various types of data and to gain an understanding of past, present and future sustainability that can be applied in diverse social–ecological settings.〈/jats:p〉〈/jats:list-item〉 〈/jats:list〉 〈/jats:p〉〈jats:p〉A free 〈jats:ext-link xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="http://onlinelibrary.wiley.com/doi/10.1002/pan3.10092/suppinfo"〉Plain Language Summary〈/jats:ext-link〉 can be found within the Supporting Information of this article.〈/jats:p〉

Description: 〈jats:title〉Abstract〈/jats:title〉〈jats:p〉 〈jats:list〉 〈jats:list-item〉〈jats:p〉Herbivory is a key process on coral reefs, which, through grazing of algae, can help sustain coral‐dominated states on frequently disturbed reefs and reverse macroalgal regime shifts on degraded ones.〈/jats:p〉〈/jats:list-item〉 〈jats:list-item〉〈jats:p〉Our understanding of herbivory on reefs is largely founded on feeding observations at small spatial scales, yet the biomass and structure of herbivore populations is more closely linked to processes which can be highly variable across large areas, such as benthic habitat turnover and fishing pressure. Though our understanding of spatiotemporal variation in grazer biomass is well developed, equivalent macroscale approaches to understanding bottom‐up and top‐down controls on herbivory are lacking.〈/jats:p〉〈/jats:list-item〉 〈jats:list-item〉〈jats:p〉Here, we integrate underwater survey data of fish abundances from four Indo‐Pacific island regions with herbivore feeding observations to estimate grazing rates for two herbivore functions, cropping (which controls turf algae) and scraping (which promotes coral settlement by clearing benthic substrate), for 72 coral reefs. By including a range of reef states, from coral to algal dominance and heavily fished to remote wilderness areas, we evaluate the influences of benthic habitat and fishing on the grazing rates of fish assemblages.〈/jats:p〉〈/jats:list-item〉 〈jats:list-item〉〈jats:p〉Cropping rates were primarily influenced by benthic condition, with cropping maximized on structurally complex reefs with high substratum availability and low macroalgal cover. Fishing was the primary driver of scraping function, with scraping rates depleted at most reefs relative to remote, unfished reefs, though scraping did increase with substratum availability and structural complexity.〈/jats:p〉〈/jats:list-item〉 〈jats:list-item〉〈jats:p〉Ultimately, benthic and fishing conditions influenced herbivore functioning through their effect on grazer biomass, which was tightly correlated to grazing rates. For a given level of biomass, we show that grazing rates are higher on reefs dominated by small‐bodied fishes, suggesting that grazing pressure is greatest when grazer size structure is truncated.〈/jats:p〉〈/jats:list-item〉 〈jats:list-item〉〈jats:p〉Stressors which cause coral declines and clear substrate for turf algae will likely stimulate increases in cropping rates, in both fished and protected areas. In contrast, scraping functions are already impaired at reefs inhabited by people, particularly where structural complexity has collapsed, indicating that restoration of these key processes will require scraper biomass to be rebuilt towards wilderness levels.〈/jats:p〉〈/jats:list-item〉 〈/jats:list〉 〈/jats:p〉〈jats:p〉A free 〈jats:ext-link xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="http://onlinelibrary.wiley.com/doi/10.1111/1365-2435.13457/suppinfo"〉Plain Language Summary〈/jats:ext-link〉 can be found within the Supporting Information of this article.〈/jats:p〉

Description: This article has been accepted for publication in Geophysical Journal International ©:The Author(s) 2020. Published by Oxford University Press on behalf of the Royal Astronomical Society. All rights reserved.Uploaded in accordance with the publisher's self-archiving policy. All rights reserved

Description: The application of a physics-based earthquake simulator to Central Italy allowed the compilation of a synthetic seismic catalogue spanning 100 000 yr, containing more than 300 000 M ≥ 4.0 simulated earthquakes, without the limitations that real catalogues suffer in terms of completeness, homogeneity and time duration. The seismogenic model upon which we applied the simulator code was derived from version 3.2.1 of the Database of Individual Seismogenic Sources (DISS; http://diss.rm.ingv.it/diss/), selecting, and modifying where appropriate, all the fault systems that are recognized in the portion of Central Italy considered in this study, with a total of 54 faults. Besides tectonic stress loading and static stress transfer as in the previous versions, the physical model on which the latest version of our simulation algorithm is based also includes the Rate and State constitutive law that helps to reproduce Omori’s law. One further improvement in our code was also the introduction of trapezoidalshaped faults that perform better than known faults. The resulting synthetic seismic catalogue exhibits typical magnitude, space and time features which are comparable to those in real observations. These features include the total seismic moment rate, the earthquake magnitude distribution, and the short- and medium-term earthquake clustering. A typical aspect of the observed seismicity in Central Italy, aswell as across thewhole Italian landmass and elsewhere, is the occurrence of earthquake sequences characterized by multiple main shocks of similar magnitude. These sequences are different from the usual earthquake clusters and aftershock sequences, since they have at least two main shocks of similar magnitude. Therefore, special attentionwas devoted to verifyingwhether the simulated catalogue includes this notable aspect. For this purpose, we developed a computer code especially for this work to count the number of multiple events contained in a seismic catalogue under a quantitative definition. We found that the last version of the simulator code produces a slightly larger number of multiple events than the previous versions, but not as large as in the real catalogue. A possible reason for this drawback is the lack of components such as pore-pressure changes due to fluid-diffusion in the adopted physical model.

Description: Published

Description: 526–542

Description: 6T. Studi di pericolosità sismica e da maremoto

Description: JCR Journal

Description: 3D earthquake locations, focal mechanisms and stress tensor distribution in a 16- month interval covering the 2018 Mt. Etna flank eruption, enabled us to investigate the relationship between magma intrusion and structural response of the volcano and shed light on the dynamic processes affecting the instability of Mt. Etna. The magma intrusion likely caused tension in the flanks of the volcano, leading to significant ground deformation and redistribution of stress on the neighbouring faults at the edge of Mt. Etna's unstable sector, encouraging the ESE sliding of the eastern flank of the volcano. Accordingly, FPSs of the post-eruptive events show strike slip faulting mechanisms, under a stress regime characterized by a maximum compressive σ1, NE-SW oriented. In this perspective, any flank eruption could temporarily enhance the sliding process of both the southern and eastern flanks of the volcano.

Description: Published

Description: 334-344

Description: 5V. Processi eruttivi e post-eruttivi

Description: JCR Journal

Description: Deep-towed geophysical surveys require precise knowledge of navigational parameters such as instrument position and orientation because navigational uncertainties reflect in the data and therefore in the inferred geophysical properties of the subseafloor. We address this issue for the case of electrical conductivity inferred from controlled source electromagnetic data. We show that the data error is laterally variable due to irregular motion during deep towing, but also due to lateral variations in conductivity, including those resulting from topography. To address this variability and quantify the data error prior to inversion, we propose a 2-D perturbation study. Our workflow enables stable and geologically reliable results for multicomponent and multifrequency inversions. An error estimation workflow is presented, which comprises the assessment of navigational uncertainties, perturbation of navigational parameters, and forward modelling of electric field amplitudes for a homogeneous and then a heterogeneous subseafloor conductivity model. Some navigational uncertainties are estimated from variations of direct measurements. Other navigational parameters required for inversion are derived from the measured quantities and their error is calculated by means of error propagation. Some navigational parameters show direct correlation with the measured electric fields. For example, the antenna dip correlates with the vertical electric field and the depth correlates with the horizontal electric field. For the perturbation study each standard deviation is added to the navigational parameters. Forward models are run for each perturbation. Amplitude deviations are summed in quadrature with the stacking error for a total, laterally varying, data error. The error estimation is repeated for a heterogeneous subseafloor model due to the large conductivity range (several orders of magnitude), which affects the forward model. The approach enables us to utilize data from several components (multiple electric fields, frequencies and receivers) in the inversion to constrain the final model and reduce ambiguity. The final model is geologically reasonable, in this case enabling the identification of conductive metal sulphide deposits on the seafloor.

Description: Marine heatwaves have been observed worldwide and are expected to increase in both frequency and intensity due to climate change. Such events may cause ecosystem reconfigurations arising from species range contraction or redistribution, with ecological, economic and social implications. Macrophytes such as the brown seaweed Fucus vesiculosus and the seagrass Zostera marina are foundation species in many coastal ecosystems of the temperate northern hemisphere. Hence, their response to extreme events can potentially determine the fate of associated ecosystems. Macrophyte functioning is intimately linked to the maintenance of photosynthesis, growth and reproduction, and resistance against pathogens, epibionts and grazers. We investigated morphological, physiological, pathological and chemical defence responses of western Baltic Sea F. vesiculosus and Z. marina populations to simulated near‐natural marine heatwaves. Along with (a) the control, which constituted no heatwave but natural stochastic temperature variability (0HW), two treatments were applied: (b) two late‐spring heatwaves (June, July) followed by a summer heatwave (August; 3HW) and (c) a summer heatwave only (1HW). The 3HW treatment was applied to test whether preconditioning events can modulate the potential sensitivity to the summer heatwave. Despite the variety of responses measured in both species, only Z. marina growth was impaired by the accumulative heat stress imposed by the 3HW treatment. Photosynthetic rate, however, remained high after the last heatwave indicating potential for recovery. Only epibacterial abundance was significantly affected in F. vesiculosus. Hence both macrophytes, and in particular F. vesiculosus, seem to be fairly tolerant to short‐term marine heatwaves at least at the intensities applied in this experiment (up to 5°C above mean temperature over a period of 9 days). This may partly be due to the fact that F. vesiculosus grows in a highly variable environment, and may have a high phenotypic plasticity.

Description: Drill sites in the southern Bay of Bengal at 3°N 91°E (International Ocean Discovery Program Expedition 362) have sampled for the first time a complete section of the Nicobar Fan and below to the oceanic crust. This generally overlooked part of the Bengal–Nicobar Fan System may provide new insights into uplift and denudation rates of the Himalayas and Tibetan Plateau. The Nicobar Fan comprises sediment gravity‐flow deposits, mostly turbidites, that alternate with hemipelagite drapes and pelagite intervals of varying thicknesses. The decimetre‐thick to metre‐thick oldest pre‐fan sediments (limestones/chalks) dated at 69 Ma are overlain by volcanic material and slowly accumulated pelagites (0.5 g cm−2 kyr−1). At Expedition 362 Site U1480, terrigenous input began in the early Miocene at ca 22.5 Ma as muds, overlain by very thin‐bedded and thin‐bedded muddy turbidites at ca 19.5 Ma. From 9.5 Ma, sand content and sediment supply sharply increase (from 1–5 to 10–50 g cm−2 kyr−1). Despite the abundant normal faulting in the Nicobar Fan compared with the Bengal Fan, it offers a better‐preserved and more homogeneous sedimentary record with fewer unconformities. The persistent connection between the two fans ceased at 0.28 Ma when the Nicobar Fan became inactive. The Nicobar Fan is a major sink for Himalaya‐derived material. This study presents integrated results of International Ocean Discovery Program Expedition 362 with older Deep Sea Drilling Project/Ocean Drilling Program/International Ocean Discovery Program sites that show that the Bengal–Nicobar Fan System experienced successive large‐scale avulsion processes that switched sediment supply between the Bengal Fan (middle Miocene and late Pleistocene) and the Nicobar Fan (late Miocene to early Pleistocene). A quantitative analysis of the submarine channels of the Nicobar Fan is also presented, including their stratigraphic frequency, showing that channel size/area and abundance peaked at ca 2 to 3 Ma, but with a distinct low at 3 to 7 Ma: the intervening stratigraphic [sub]unit was a time of reduced sediment accumulation rates.

Description: The generation of reliable age models for palaeoenvironmental and archaeological records in the Eurasian Arctic is often problematic when using conventional dating techniques. Tephrochronology can potentially improve the chronologies of such records and synchronise disparate sedimentary archives. However, to date, systematic tephra studies are lacking for this region. This paper presents the first cryptotephra data from the White Sea region (northwestern Russia) based on a peat core spanning the past ~1800 years. We identify seven geochemical glass populations that derive from six Icelandic volcanoes and correlate four of them to north European tephra isochrons; these include Askja ad 1875, the basaltic component of the ad 877 Landnám tephra, and tephras BTD‐15 (c. ad 1750–1650) and SL‐2/SB‐2 (ad 803–767) from unknown eruptions of Katla and Snæfellsjökull, respectively. The remaining three populations originate from Grímsvötn, Hekla and Katla; however, their attribution to individual eruptions remains ambiguous. These findings highlight the potential to extend the Late Holocene tephrochronological framework of northern Europe to the west Eurasian Arctic. The detection of at least three basaltic tephras in the core suggests that basaltic shards can be transported over larger distances than previously known and that peatlands are well suited to preserve such components.

Description: Climate engineering (CE) measures are increasingly discussed when dealing with the adverse impacts of climate change. While much research has focused on individual methods, few studies attempt to compare and rank the effectiveness of these measures. Furthermore, model uncertainties are seldom acknowledged and lesser still, estimated when CE scenarios are assessed. In this work, we quantify the variance in outcomes due to poorly constrained model parameters under several idealized CE scenarios. The four scenarios considered are (1) warming under the high emission scenario Representative Concentration Pathway 8.5 without CE applied and the same emission scenario with (2) afforestation,(3) solar radiation management, and (4) artificial ocean alkalinization. By considering the parametric uncertainty in model outputs, we demonstrate the problems with comparing these scenarios using a single parameter setting. Using statistical emulation, we estimate the probability distributions of several model outcomes. Based on such distributions, we suggest an approach to ranking the effectiveness of the scenarios considered according to their probability of avoiding climate thresholds.

Description: A new estimate of Agulhas leakage transport is calculated using profiling floats and drifters. Since Richardson's seminal estimate of 15 Sv in 2007, the number of floats and drifters passing through the Agulhas Current has quadrupled. Within uncertainties we find the same leakage percentages as Richardson, with 34% of drifters leaking at the surface and 21% of floats leaking at 1,000 m depth. We find that the drifters tend to follow a northward leakage pathway via the Benguela Current compared to the northwestward leakage pathway of the floats along the Agulhas Ring corridor. We simulate the isobaric and profiling behavior of the floats and drifters using two high resolution models and two offline Lagrangian tracking tools, quantifying for the first time the sampling biases associated with the observations. We find that the isobaric bias cannot be robustly simulated but likely causes an underestimate of observed leakage by one or two Sverdrups. The profiling behavior of the floats causes no significant bias in the leakage. Fitting a simulated vertical leakage profile to the observed leakage percentages from the floats and drifters and using the mean Agulhas transport observed by a moored array at 34°S we find an improved Agulhas leakage transport of 21.3 Sv, with an estimated error of 4.7 Sv. Our new leakage transport is higher primarily because we account for leakage at depths down to 2,000 m, while Richardson considered only the top 1,000 m of the water column.

Description: Upwelling ocean currents associated with oxygen minimum zones (OMZs) supply nutrients fuelling intense marine productivity. Perturbations in the extent and intensity of OMZs are projected in the future, but it is currently uncertain how this will impact fluxes of redox‐sensitive trace metal micronutrients to the surface ocean. Here we report seawater concentrations of Fe, Mn, Co, Cd, and Ni alongside the redox indicator iodide/iodate in the Peruvian OMZ during the 2015 El Niño event. The El Niño drove atypical upwelling of oxygen‐enriched water over the Peruvian Shelf, resulting in oxidized iodine and strongly depleted Fe (II), total dissolved Fe, and reactive particulate Fe concentrations relative to non‐El Niño conditions. Observations of Fe were matched by the redox‐sensitive micronutrients Co and Mn, but not by non‐redox‐sensitive Cd and Ni. These observations demonstrate that oxygenation of OMZs significantly reduces water column inventories of redox‐sensitive micronutrients, with potential impacts on ocean productivity. Plain Language Summary Some trace metals, including iron, are essential micronutrients for phytoplankton growth. However, the solubility of iron is very low under oxygenated conditions. Consequently, restricted iron availability in oxygen‐rich seawater can limit phytoplankton growth in the ocean, including in the Eastern Tropical South Pacific. Under typical conditions, depleted oxygen on the South American continental shelf is generally thought to enhance iron supply to the ocean, fuelling phytoplankton productivity in overlying waters. However, the impact of changes in oxygenation, which are predicted to occur in the future, are not known. The 2015 El Niño event led to unusually high oxygen on the Peruvian shelf, offering a system‐scale test on how oxygen influences seawater iron concentrations. We show that El Niño‐driven oxygenation resulted in marked decreases in iron and other metals sensitive to oxygen (cobalt and manganese), whilst metals not sensitive to oxygen (cadmium and nickel) were unaffected. The measured reductions in iron may have led to decreased phytoplankton productivity.

Description: The causes of the seasonal cycle of vertical turbulent cooling at the base of the mixed layer are assessed using observations from moored buoys in the tropical Atlantic Intertropical Convergence Zone (ITCZ) (4°N, 23°W) and trade wind (15°N, 38°W) regions together with mixing parameterizations and a one-dimensional model. At 4°N the parameterized turbulent cooling rates during 2017–2018 and 2019 agree with indirect estimates from the climatological mooring heat budget residual: both show mean cooling of 25–30 W m (Formula presented.) during November–July, when winds are weakest and the mixed layer is thinnest, and 0–10 W m (Formula presented.) during August–October. Mixing during November–July is driven by variability on multiple time scales, including subdiurnal, near-inertial, and intraseasonal. Shear associated with tropical instability waves (TIWs) is found to generate mixing and monthly mean cooling of 15–30 W m (Formula presented.) during May–July in 2017 and 2019. At 15°N the seasonal cycle of turbulent cooling is out of phase compared to 4°N, with largest cooling of up to 60 W m (Formula presented.) during boreal fall. However, the relationships between wind speed, mixed layer depth, and turbulent mixing are similar: weaker mean winds and a thinner mixed layer in the fall are associated with stronger mixing and turbulent cooling of SST. These results emphasize the importance of seasonal modulations of mixed layer depth at both locations and shear from TIWs at 4°N.

Description: Submarine groundwater discharge (SGD) into coastal areas is a common global phenomenon and is rapidly gaining scientific interest due to its influence on marine ecology, the coastal sedimentary environment and its potential as a future freshwater resource. We conducted an integrated study of hydroacoustic surveys combined with geochemical porewater and water column investigations at a well‐known groundwater seep site in Eckernförde Bay (Germany). We aim to better constrain the effects of shallow gas and SGD on high frequency multibeam backscatter data and to present acoustic indications for submarine groundwater discharge. Our high‐quality hydroacoustic data reveal hitherto unknown internal structures within the pockmarks in Eckernförde Bay. Using precisely positioned sediment core samples, our hydroacoustic‐geochemical approach can differentiate intra‐pockmark regimes that were formerly assigned to pockmarks of a different nature. We demonstrate that high‐frequency multibeam data, in particular the backscatter signals, can be used to detect shallow free gas in areas of enhanced groundwater advection in muddy sediments. Intriguingly, our data reveal relatively small (typically 〈15 m across) pockmarks within the much larger, previously mapped, pockmarks. The small pockmarks, which we refer to as “intra‐pockmarks”, have formed due to the localized ascent of gas and groundwater; they manifest themselves as a new type of ‘eyed’ pockmarks, revealed by their acoustic backscatter pattern. Our data suggest that, in organic‐rich muddy sediments, morphological lows combined with a strong multibeam backscatter signal can be indicative of free shallow gas and subsequent advective groundwater flow.

Description: A major surface circulation feature of the Arctic Ocean is the Transpolar Drift (TPD), a current that transports river‐influenced shelf water from the Laptev and East Siberian Seas toward the center of the basin and Fram Strait. In 2015, the international GEOTRACES program included a high‐resolution pan‐Arctic survey of carbon, nutrients, and a suite of trace elements and isotopes (TEIs). The cruises bisected the TPD at two locations in the central basin, which were defined by maxima in meteoric water and dissolved organic carbon concentrations that spanned 600 km horizontally and ~25‐50 m vertically. Dissolved TEIs such as Fe, Co, Ni, Cu, Hg, Nd, and Th, which are generally particle‐reactive but can be complexed by organic matter, were observed at concentrations much higher than expected for the open ocean setting. Other trace element concentrations such as Al, V, Ga, and Pb were lower than expected due to scavenging over the productive East Siberian and Laptev shelf seas. Using a combination of radionuclide tracers and ice drift modeling, the transport rate for the core of the TPD was estimated at 0.9 ± 0.4 Sv (106 m3 s‐1). This rate was used to derive the mass flux for TEIs that were enriched in the TPD, revealing the importance of lateral transport in supplying materials beneath the ice to the central Arctic Ocean and potentially to the North Atlantic Ocean via Fram Strait. Continued intensification of the Arctic hydrologic cycle and permafrost degradation will likely lead to an increase in the flux of TEIs into the Arctic Ocean.

Description: Cichlid fishes’ famous diversity in body coloration is accompanied by a highly diverse and complex visual system. Although cichlids possess an unusually high number of seven cone opsin genes, they express only a subset of these during their ontogeny, accounting for their astonishing interspecific variation in visual sensitivities. Much of this diversity is thought to have been shaped by natural selection as cichlids inhabit a variety of habitats with distinct light environments. Also, sexual selection might have contributed to the observed visual diversity, and sexual dimorphism in coloration potentially co‐evolved with sexual dimorphism in opsin expression. We investigated sex‐specific opsin expression of several cichlids from Africa and the Neotropics and collected and integrated datasets on sex‐specific body coloration, species‐specific visual sensitivities, lens transmission and habitat light properties for some of them. We comparatively analyzed this wide range of molecular and ecological data, illustrating how integrative approaches can address specific questions on the factors and mechanisms driving diversification, and the evolution of cichlid vision in particular. We found that both sexes expressed opsins at the same levels ‐ even in sexually dimorphic cichlid species – which argues against coevolution of sexual dichromatism and differences in sex‐specific visual sensitivity. Rather, a combination of environmental light properties and body coloration shaped the diversity in spectral sensitivities among cichlids. We conclude that although cichlids are particularly colorful and diverse and often sexually dimorphic, it would appear that natural rather than sexual selection is a more powerful force driving visual diversity in this hyper‐diverse lineage.

Description: Key Points: - The Kamchatka arc lavas show across-arc variations in chalcophile elements, suggesting that the amount of fluid decreases with depth - Slab-derived fluids have a negligible contribution to the Li budget of the Kamchatka arc lavas - CKD lavas have high U/Th, Li/Y, La/Sm and B/Nb ratios, indicating that lawsonite breakdown reaction dominates the water release Chalcophile elements and lithium (Li) isotopes were measured on lavas from a 220 km transect across the Kamchatka arc in order to investigate the fluid variations below arc volcanoes and to trace the geochemical behaviour of Li in convergent plate margins. From the Eastern Volcanic Front (EVF), through the Central Kamchatka Depression (CKD), into the Sredinny Range (SR) volcanic zones, chalcophile element ratios (e.g., As/Ce and Sb/Ce) show clear across‐arc variations, decreasing (e.g., As/Ce: 0.20 to 0.03 and Sb/Ce: 0.013 to 0.002) with increasing depth above the slab (110 to 400 km). This clearly indicates a gradually decreasing influx of slab‐derived fluids added to the mantle wedge as the slab subducts below Kamchatka. In addition, the anomalously high U/Th, La/Sm and B/Nb ratios in the CKD lavas suggest lawsonite breakdown reaction dominates the fluid release in this area. However, Li/Y (0.07 to 1.78) and δ7Li (+1.8 to +5.4‰, with an exception of +8.6‰ in CKD) show limited variations and values similar to the MORB mantle. A dehydration model suggests that slab‐derived fluids, which are characterized by high Li concentration and high δ7Li, do not control the Li budget in Kamchatka arc lavas. Therefore, the isotopic heavy Li from slab‐derived fluids likely equilibrates in the sub‐arc mantle, which acts as a buffer for Li systematics. In addition, based on the Li isotopic signatures of Klyuchevskoy volcano, our study demonstrates insignificant Li isotopic fractionation during mantle melting and subsequent differentiation.

Description: Diverse microbial ecosystems underpin life in the sea. Among these microbes are many unicellular eukaryotes that span the diversity of the eukaryotic tree of life. However, genetic tractability has been limited to a few species, which do not represent eukaryotic diversity or environmentally relevant taxa. Here, we report on the development of genetic tools in a range of protists primarily from marine environments. We present evidence for foreign DNA delivery and expression in 13 species never before transformed and for advancement of tools for eight other species, as well as potential reasons for why transformation of yet another 17 species tested was not achieved. Our resource in genetic manipulation will provide insights into the ancestral eukaryotic lifeforms, general eukaryote cell biology, protein diversification and the evolution of cellular pathways.

Description: Communities are shaped by scale dependent processes. To study the diversity and variation of microbial communities across scales, the invasive and widespread seaweed Agarophyton vermiculophyllum presents a unique opportunity. We characterized pro‐ and eukaryotic communities associated with this holobiont across its known distribution range, which stretches over the northern hemisphere. Our data reveal that community composition and diversity in the holobiont vary at local but also larger geographic scales. While processes acting at the local scale (i.e., within population) are the main structuring drivers of associated microbial communities, changes in community composition also depend on processes acting at larger geographic scales. Interestingly, the largest analysed scale (i.e., native and non‐native ranges) explained variation in the prevalence of predicted functional groups, which could suggest a functional shift in microbiota occurred over the course of the invasion process. While high variability in microbiota at the local scale supports A. vermiculophyllum to be a generalist host, we also identified a number of core taxa. These geographically independent holobiont members imply that cointroduction of specific microbiota may have additionally promoted the invasion process.

Description: Atlantic cod (Gadus morhua) is one of the most important fish species in northern Europe for several reasons including its predator status in marine ecosystems, its historical role in fisheries, its potential in aquaculture and its strong public profile. However, due to over-exploitation in the North Atlantic and changes in the ecosystem, many cod populations have been reduced in size and genetic diversity. Cod populations in the Baltic Proper, Kattegat and North Sea have been analyzed using a species specific single nucleotide polymorphism (SNP) array. Using a subset of 8,706 SNPs, moderate genetic differences were found between subdivisions in three traditionally delineated cod management stocks: Kattegat, western and eastern Baltic. However, an FST measure of population differentiation based on allele frequencies from 588 outlier loci for 2 population groups, one including 5 western and the other 4 eastern Baltic populations, indicated high genetic differentiation. In this paper, differentiation has been demonstrated not only between, but also within western and eastern Baltic cod stocks for the first time, with salinity appearing to be the most important environmental factor influencing the maintenance of cod population divergence between the western and eastern Baltic Sea.

Description: The Endeavour Segment of the Juan de Fuca Ridge is well known for its abundance of hydrothermal vents and chimneys. One-meter scale multibeam mapping data collected by an autonomous undersea vehicle revealed 572 chimneys along the central 14 km of the segment, although only 47 are named and known to be active. Hydrothermal deposits are restricted to the axial graben and the near-rims of the graben above a seismically mapped axial magma lens. The sparse eruptive activity on the segment during the last 4,300 years has not buried inactive chimneys, as occurs at more magmatically robust mid-ocean ridges.

Description: Although the core velocity of the Atlantic North Equatorial Undercurrent (NEUC) is low (0.1−0.3 m s−1), it has been suggested to act as an important oxygen supply route towards the oxygen minimum zone in the eastern tropical North Atlantic. For the first time, the intraseasonal to interannual NEUC variability and its impact on oxygen are investigated based on shipboard and moored velocity observations around 5°N, 23°W. In contrast to previous studies that were mainly based on models or hydrographic data, we find hardly any seasonal cycle of NEUC transports in the central Atlantic. The NEUC transport variability is instead dominated by sporadic intraseasonal events. Only some of these events are associated with high oxygen levels suggesting an occasional eastward oxygen supply by NEUC transport events. Nevertheless, they likely contribute to the local oxygen maximum in the mean shipboard section along 23°W at the NEUC core position.

Description: The development and application of geochemical techniques to identify redox conditions in modern and ancient aquatic environments has intensified over recent years. Iron (Fe) speciation has emerged as one of the most widely used procedures to distinguish different redox regimes in both the water column and sediments, and is the main technique used to identify oxic, ferruginous (anoxic, Fe(II) containing) and euxinic (anoxic, sulfidic) water column conditions. However, an international sediment reference material has never been developed. This has led to concern over the consistency of results published by the many laboratories that now utilise the technique. Here, we report an interlaboratory comparison of four Fe speciation reference materials for palaeoredox analysis, which span a range of compositions and reflect deposition under different redox conditions. We provide an update of extraction techniques used in Fe speciation and assess the effects of both test portion mass, and the use of different analytical procedures, on the quantification of different Fe fractions in sedimentary rocks. While atomic absorption spectroscopy and inductively coupled plasma‐optical emission spectrometry produced comparable Fe measurements for all extraction stages, the use of ferrozine consistently underestimated Fe in the extraction step targeting mixed ferrous–ferric minerals such as magnetite. We therefore suggest that the use of ferrozine is discontinued for this Fe pool. Finally, we report the combined data of four independent Fe speciation laboratories to characterise the Fe speciation composition of the reference materials. These reference materials are available to the community to provide an essential validation of in‐house Fe speciation measurements.

Description: The Sea of Galilee in northeast Israel is a freshwater lake filling a morphological depression along the Dead Sea Fault. It is located in a tectonically complex area, where a N-S main fault system intersects secondary fault patterns non-univocally interpreted by previous reconstructions. A set of multiscale geophysical, geochemical and seismological data, reprocessed or newly collected, was analysed to unravel the interplay between shallow tectonic deformations and geodynamic processes. The result is a neotectonic map highlighting major seismogenic faults in a key region at the boundary between the Africa/Sinai and Arabian plates. Most active seismogenic displacement occurs along NNW-SSE oriented transtensional faults. This results in a left-lateral bifurcation of the Dead Sea Fault forming a rhomb-shaped depression we named the Capharnaum Trough, located off-track relative to the alleged principal deformation zone. Low-magnitude (ML = 3–4) epicentres accurately located during a recent seismic sequence are aligned along this feature, whose activity, depth and regional importance is supported by geophysical and geochemical evidence. This case study, involving a multiscale/multidisciplinary approach, may serve as a reference for similar geodynamic settings in the world, where unravelling geometric and kinematic complexities is challenging but fundamental for reliable earthquake hazard assessments.

Description: Resting egg production is considered the most common form of dormancy in aquatic invertebrates. Given that many taxa at least partially terminate resting egg state using environmental cues, knowledge on environmental drivers of hatching success is important, particularly within the context of climate change and environmental degradation. Fairy shrimp (anostracans) are temporary wetland specialists that are reliant on resting egg production for population persistence. Temporary wetlands are common in many arid regions projected to experience increases in temperature, and in areas often compromised by human‐mediated activities. In this study, we assessed the combined effects of light and temperature on the hatching success of Streptocephalus cafer (Anostraca) dormant eggs from temporary wetlands in an arid environment. Both temperature and light altered hatching success, with emergent effects evident. Light caused a significant threefold increase in hatching success overall, while temperature effects were non‐linear, with hatching optimised at 27°C, and especially under light conditions. These results are discussed within the context of shifting climates and disturbances to temporary wetland ecosystems.

Description: Bathymetry (seafloor depth), is a critical parameter providing the geospatial context for a multitude of marine scientific studies. Since 1997, the International Bathymetric Chart of the Arctic Ocean (IBCAO) has been the authoritative source of bathymetry for the Arctic Ocean. IBCAO has merged its efforts with the Nippon Foundation-GEBCO-Seabed 2030 Project, with the goal of mapping all of the oceans by 2030. Here we present the latest version (IBCAO Ver. 4.0), with more than twice the resolution (200 × 200 m versus 500 × 500 m) and with individual depth soundings constraining three times more area of the Arctic Ocean (∼19.8% versus 6.7%), than the previous IBCAO Ver. 3.0 released in 2012. Modern multibeam bathymetry comprises ∼14.3% in Ver. 4.0 compared to ∼5.4% in Ver. 3.0. Thus, the new IBCAO Ver. 4.0 has substantially more seafloor morphological information that offers new insights into a range of submarine features and processes; for example, the improved portrayal of Greenland fjords better serves predictive modelling of the fate of the Greenland Ice Sheet.

Description: Most molluscs possess shells, constructed from a vast array of microstructures and architectures. The fully formed shell is composed of calcite or aragonite. These CaCO3 crystals form complex biocomposites with proteins, which although typically less than 5% of total shell mass, play significant roles in determining shell microstructure. Despite much research effort, large knowledge gaps remain in how molluscs construct and maintain their shells, and how they produce such a great diversity of forms. Here we synthesize results on how shell shape, microstructure, composition and organic content vary among, and within, species in response to numerous biotic and abiotic factors. At the local level, temperature, food supply and predation cues significantly affect shell morphology, whilst salinity has a much stronger influence across latitudes. Moreover, we emphasize how advances in genomic technologies [e.g. restriction site-associated DNA sequencing (RAD-Seq) and epigenetics] allow detailed examinations of whether morphological changes result from phenotypic plasticity or genetic adaptation, or a combination of these. RAD-Seq has already identified single nucleotide polymorphisms associated with temperature and aquaculture practices, whilst epigenetic processes have been shown significantly to modify shell construction to local conditions in, for example, Antarctica and New Zealand. We also synthesize results on the costs of shell construction and explore how these affect energetic trade-offs in animal metabolism. The cellular costs are still debated, with CaCO3 precipitation estimates ranging from 1–2 J/mg to 17–55 J/mg depending on experimental and environmental conditions. However, organic components are more expensive (~29 J/mg) and recent data indicate transmembrane calcium ion transporters can involve considerable costs. This review emphasizes the role that molecular analyses have played in demonstrating multiple evolutionary origins of biomineralization genes. Although these are characterized by lineage-specific proteins and unique combinations of co-opted genes, a small set of protein domains have been identified as a conserved biomineralization tool box. We further highlight the use of sequence data sets in providing candidate genes for in situ localization and protein function studies. The former has elucidated gene expression modularity in mantle tissue, improving understanding of the diversity of shell morphology synthesis. RNA interference (RNAi) and clustered regularly interspersed short palindromic repeats - CRISPR-associated protein 9 (CRISPR-Cas9) experiments have provided proof of concept for use in the functional investigation of mollusc gene sequences, showing for example that Pif (aragonite-binding) protein plays a significant role in structured nacre crystal growth and that the Lsdia1 gene sets shell chirality in Lymnaea stagnalis. Much research has focused on the impacts of ocean acidification on molluscs. Initial studies were predominantly pessimistic for future molluscan biodiversity. However, more sophisticated experiments incorporating selective breeding and multiple generations are identifying subtle effects and that variability within mollusc genomes has potential for adaption to future conditions. Furthermore, we highlight recent historical studies based on museum collections that demonstrate a greater resilience of molluscs to climate change compared with experimental data. The future of mollusc research lies not solely with ecological investigations into biodiversity, and this review synthesizes knowledge across disciplines to understand biomineralization. It spans research ranging from evolution and development, through predictions of biodiversity prospects and future-proofing of aquaculture to identifying new biomimetic opportunities and societal benefits from recycling shell products.

Description: Climate change has led to a ~ 40% reduction in summer Arctic sea-ice cover extent since the 1970s. Resultant increases in light availability may enhance phytoplankton production. Direct evidence for factors currently constraining summertime phytoplankton growth in the Arctic region is however lacking. GEOTRACES cruise GN05 conducted a Fram Strait transect from Svalbard to the NE Greenland Shelf in summer 2016, sampling for bioessential trace metals (Fe, Co, Zn, Mn) and macronutrients (N, Si, P) at ~ 79°N. Five bioassay experiments were conducted to establish phytoplankton responses to additions of Fe, N, Fe + N and volcanic dust. Ambient nutrient concentrations suggested N and Fe were deficient in surface seawater relative to typical phytoplankton requirements. A west-to-east trend in the relative deficiency of N and Fe was apparent, with N becoming more deficient towards Greenland and Fe more deficient towards Svalbard. This aligned with phytoplankton responses in bioassay experiments, which showed greatest chlorophyll-a increases in + N treatment near Greenland and + N + Fe near Svalbard. Collectively these results suggest primary N limitation of phytoplankton growth throughout the study region, with conditions potentially approaching secondary Fe limitation in the eastern Fram Strait. We suggest that the supply of Atlantic-derived N and Arctic-derived Fe exerts a strong control on summertime nutrient stoichiometry and resultant limitation patterns across the Fram Strait region.

Description: In continental settings, seismic failure is generally restricted to crustal depth. Crustal structure is therefore an important proxy to evaluate seismic hazard of continental fault systems. Here we present a seismic velocity model across the Gibraltar Arc System, from the Eurasian Betics Range (South Iberian margin), across offshore East Alboran and Pytheas (African margin) basins, and ending onshore in North Morocco. Our results reveal the nature and configuration of the crust supporting the coexistence of three different crustal domains: the continental crust of the Betics, the continental crust of the Pytheas Basin (south Alboran Basin) and onshore Morocco, and a distinct domain formed of magmatic arc crust under the East Alboran Basin. The magmatic arc under the East Alboran Basin is characterized by a velocity structure containing a relatively high‐velocity lower crust (~7 km/s) bounded at the top and base by reflections. The lateral extension of this crust is mapped integrating a second perpendicular wide‐angle seismic profile along the Eastern Alboran basin, together with basement samples, multibeam bathymetry, and a grid of deep‐penetrating multichannel seismic profiles. The transition between crustal domains is currently unrelated to extensional and magmatic processes that formed the basin. The abrupt transition zones between the different crustal domains support that they are bounded by crustal‐scale active fault systems that reactivate inherited structures. Seismicity in the area is constrained to upper‐middle crust depths, and most earthquakes nucleate outside of the magmatic arc domain. Key Points New velocity model reveals the lithospheric structure under the Betics (South Iberia), the Alboran Basin and the North African margin The East Alboran Basin is floored by magmatic arc crust, while the southern area of the Alboran Basin is floored by continental crust Seismic activity is constrained to the upper‐middle continental crust. Crustal domains are likely bounded by active faults

Description: We constrain the lithospheric mantle density of the North China Craton (NCC) at both in situ and standard temperature‐pressure (STP) conditions from gravity data. The lithosphere‐asthenosphere boundary (LAB) depth is constrained by our new thermal model, which is based on a new regional heat flow data set and a recent regional crustal model NCcrust. The new thermal model shows that the thermal lithosphere thickness is 〈120 km in most of the NCC, except for the northern and southern parts with the maximum depth of 170 km. The gravity calculations reveal a highly heterogeneous density structure of the lithospheric mantle with in situ and STP values of 3.22–3.29 and 3.32–3.40 g/cm3, respectively. Thick and reduced‐density cratonic‐type lithosphere is preserved mostly in the southern NCC. Most of the Eastern Block has a thin (90–140 km) and high‐density lithospheric mantle. Most of the Western Block has a high‐density lithospheric mantle and a thin (80–110 km) lithosphere typical of Phanerozoic regions, which suggests that the Archean lithosphere is no longer present there. We conclude that in almost the entire NCC the lithosphere has lost its cratonic characteristics by geodynamic processes that include, but are not limited to, the Paleozoic closure of the Paleo‐Asian Ocean in the north, the Mesozoic Yangtze Craton flat subduction in the south, the Mesozoic Pacific subduction in the east, the Cenozoic remote response to the Indian‐Eurasian collision in the west, and the Cenozoic extensional tectonics (possibly associated with the slab roll‐back) in the center.

Description: Most animals, including sponges (Porifera), have species-specific microbiomes. Which genetic or environmental factors play major roles structuring the microbial community at the intraspecific level in sponges is, however, largely unknown. In this study, we tested whether geographic location or genetic structure of conspecific sponges influences their microbial assembly. For that, we used three sponge species with different rates of gene flow, and collected samples along their entire distribution range (two from the Mediterranean and one from the Southern Ocean) yielding a total of 393 samples. These three sponge species have been previously analysed by microsatellites or single nucleotide polymorphisms, and here we investigate their microbiomes by amplicon sequencing of the microbial 16S rRNA gene. The sponge Petrosia ficiformis, with highly isolated populations (low gene flow), showed a stronger influence of the host genetic distance on the microbial composition than the spatial distance. Host-specificity was therefore detected at the genotypic level, with individuals belonging to the same host genetic cluster harbouring more similar microbiomes than distant ones. On the contrary, the microbiome of Ircinia fasciculata and Dendrilla antarctica - both with weak population structure (high gene flow) - seemed influenced by location rather than by host genetic distance. Our results suggest that in sponge species with high population structure, the host genetic cluster influence the microbial community more than the geographic location.

Description: Observations and reanalysis products are used to investigate the substantial weakening in the southeastern tropical Atlantic sea‐surface temperature (SST) variability since 2000. Relative to 1982‐1999, the March‐April‐May SST variability in the Angola‐Benguela area (ABA) has decreased by more than 30 %. Both equatorial remote forcing and local forcing are known to play an important role in driving SST variability in the ABA. Compared to 1982‐1999, since 2000 equatorial remote forcing had less influence on ABA SSTs whereas local forcing has become more important. In particular, the robust correlation that existed between the equatorial zonal wind stress and the ABA SSTs has substantially weakened, suggesting less influence of Kelvin waves on ABA SSTs. Moreover, the strong correlation linking the South Atlantic Anticyclone and the ABA SSTs has reduced. Finally, multidecadal surface warming of the ABA could also have played a role in the weakening of the interannual SST variability.

Description: More than half of Earth’s freshwater resources are held by the Antarctic Ice Sheet, which thus represents by far the largest potential source for global sea-level rise under future warming conditions1. Its long-term stability determines the fate of our coastal cities and cultural heritage. Feedbacks between ice, atmosphere, ocean, and the solid Earth give rise to potential nonlinearities in its response to temperature changes. So far, we are lacking a comprehensive stability analysis of the Antarctic Ice Sheet for different amounts of global warming. Here we show that the Antarctic Ice Sheet exhibits a multitude of temperature thresholds beyond which ice loss is irreversible. Consistent with palaeodata2 we find, using the Parallel Ice Sheet Model3,4,5, that at global warming levels around 2 degrees Celsius above pre-industrial levels, West Antarctica is committed to long-term partial collapse owing to the marine ice-sheet instability. Between 6 and 9 degrees of warming above pre-industrial levels, the loss of more than 70 per cent of the present-day ice volume is triggered, mainly caused by the surface elevation feedback. At more than 10 degrees of warming above pre-industrial levels, Antarctica is committed to become virtually ice-free. The ice sheet’s temperature sensitivity is 1.3 metres of sea-level equivalent per degree of warming up to 2 degrees above pre-industrial levels, almost doubling to 2.4 metres per degree of warming between 2 and 6 degrees and increasing to about 10 metres per degree of warming between 6 and 9 degrees. Each of these thresholds gives rise to hysteresis behaviour: that is, the currently observed ice-sheet configuration is not regained even if temperatures are reversed to present-day levels. In particular, the West Antarctic Ice Sheet does not regrow to its modern extent until temperatures are at least one degree Celsius lower than pre-industrial levels. Our results show that if the Paris Agreement is not met, Antarctica’s long-term sea-level contribution will dramatically increase and exceed that of all other sources.

Description: The quasi‐biennial oscillation (QBO) of the equatorial zonal wind leads to zonally symmetric temperature variations in the stratosphere that descend downward. Here we investigate the QBO‐induced temperature anomalies in the tropical tropopause layer (TTL) and detect pronounced longitudinal variations of the signal. In addition, the QBO temperature anomalies show a strong seasonal variability. The magnitude of these seasonal and longitudinal QBO variations is comparable to the magnitude of the well‐known zonal mean QBO signal in the TTL. At the cold point tropopause, the strongest QBO variations of around ±1.6 K are found over regions of active convection such as the West Pacific and Africa during boreal winter. The weakest QBO variations of ±0.25 K are detected over the East Pacific during boreal summer, while the zonal mean signal ranges around ±0.7 K. The longitudinal variations are associated with enhanced convective activity that occurs during QBO cold phases and locally enhances the cold anomalies.

Description: We present a mechanism for self‐sustained ocean circulation changes that cause abrupt temperature changes over Greenland in a multimillennial climate model simulation with glacial CO2 concentrations representative of Marine Isotope Stage 3. The Atlantic meridional overturning circulation (AMOC) and the subpolar gyre (SPG) oscillate on millennial time scales. When the AMOC is strong, the SPG is weak and contracted; when the AMOC is weak, the SPG is strong and extensive. The coupling between the two systems via wind‐driven and density‐driven feedbacks is key to maintaining the oscillations. The SPG controls the transport of heat and salt into the deep‐water formation sites and thus controls the AMOC strength. The strength and location of the deep‐water formation affect the density‐driven part of the SPG and thus control the mean strength and extent of the SPG. This mechanism supports the hypothesis that coupled ocean‐ice‐atmosphere interactions could have triggered abrupt glacial climate change.

Description: Past changes in ocean 14C disequilibria have been suggested to reflect the Southern Ocean control on global exogenic carbon cycling. Yet, the volumetric extent of the glacial carbon pool and the deglacial mechanisms contributing to release remineralized carbon, particularly from regions with enhanced mixing today, remain insufficiently constrained. Here, we reconstruct the deglacial ventilation history of the South Indian upwelling hotspot near Kerguelen Island, using high-resolution 14C-dating of smaller-than-conventional foraminiferal samples and multi-proxy deep-ocean oxygen estimates. We find marked regional differences in Southern Ocean overturning with distinct South Indian fingerprints on (early de-)glacial atmospheric CO2 change. The dissipation of this heterogeneity commenced 14.6 kyr ago, signaling the onset of modern-like, strong South Indian Ocean upwelling, likely promoted by rejuvenated Atlantic overturning. Our findings highlight the South Indian Ocean’s capacity to influence atmospheric CO2 levels and amplify the impacts of inter-hemispheric climate variability on global carbon cycling within centuries and millennia.

Description: Benthic fauna refers to all fauna that live in or on the seafloor, which researchers typically divide into size classes meiobenthos (32/64 µm–0.5/1 mm), macrobenthos (250 µm–1 cm), and megabenthos (〉1 cm). Benthic fauna play important roles in bioturbation activity, mineralization of organic matter, and in marine food webs. Evaluating their role in these ecosystem functions requires knowledge of their global distribution and biomass. We therefore established the BenBioDen database, the largest open-access database for marine benthic biomass and density data compiled so far. In total, it includes 11,792 georeferenced benthic biomass and 51,559 benthic density records from 384 and 600 studies, respectively. We selected all references following the procedure for systematic reviews and meta-analyses, and report biomass records as grams of wet mass, dry mass, or ash-free dry mass, or carbon per m2 and as abundance records as individuals per m2. This database provides a point of reference for future studies on the distribution and biomass of benthic fauna.

Description: It is fundamentally important for many animal ecologists to quantify the costs of animal activities, although it is not straightforward to do so. The recording of triaxial acceleration by animal-attached devices has been proposed as a way forward for this, with the specific suggestion that dynamic body acceleration (DBA) be used as a proxy for movement-based power. Dynamic body acceleration has now been validated frequently, both in the laboratory and in the field, although the literature still shows that some aspects of DBA theory and practice are misunderstood. Here, we examine the theory behind DBA and employ modelling approaches to assess factors that affect the link between DBA and energy expenditure, from the deployment of the tag, through to the calibration of DBA with energy use in laboratory and field settings. Using data from a range of species and movement modes, we illustrate that vectorial and additive DBA metrics are proportional to each other. Either can be used as a proxy for energy and summed to estimate total energy expended over a given period, or divided by time to give a proxy for movement-related metabolic power. Nonetheless, we highlight how the ability of DBA to predict metabolic rate declines as the contribution of non-movement-related factors, such as heat production, increases. Overall, DBA seems to be a substantive proxy for movement-based power but consideration of other movement-related metrics, such as the static body acceleration and the rate of change of body pitch and roll, may enable researchers to refine movement-based metabolic costs, particularly in animals where movement is not characterized by marked changes in body acceleration.

Description: We conduct the seismic signal analysis on vintage and recently collected multichannel seismic reflection profiles from the Ionian Basin to characterize the deep basin Messinian evaporites. These evaporites were deposited in deep and marginal Mediterranean sedimentary basins as a consequence of the “salinity crisis” between 5.97 and 5.33 Ma, a basin‐wide oceanographic and ecological crisis whose origin remains poorly understood. The seismic markers of the Messinian evaporites in the deep Mediterranean basins can be divided in two end‐members, one of which is the typical “trilogy” of gypsum and clastics (Lower Unit – LU), halite (Mobile Unit – MU) and upper anhydrite and marl layers (Upper Unit – UU) traced in the Western Mediterranean Basins. The other end‐member is a single MU unit subdivided in seven sub‐units by clastic interlayers located in the Levant Basin. The causes of these different seismic expressions of the Messinian salinity crisis (MSC) appear to be related to a morphological separation between the two basins by the structural regional sill of the Sicily Channel. With the aid of velocity analyses and seismic imaging via prestack migration in time and depth domains, we define for the first time the seismic signature of the Messinian evaporites in the deep Ionian Basin, which differs from the known end‐members. In addition, we identify different evaporitic depositional settings suggesting a laterally discontinuous deposition. With the information gathered we quantify the volume of evaporitic deposits in the deep Ionian Basin as 500,000 km3 ± 10%. This figure allows us to speculate that the total volume of salts in the Mediterranean basin is larger than commonly assumed. Different depositional units in the Ionian Basin suggest that during the MSC it was separated from the Western Mediterranean by physical thresholds, from the Po Plain/Northern Adriatic Basin, and the Levant Basin, likely reflecting different hydrological and climatic conditions. Finally, the evidence of erosional surfaces and V‐shaped valleys at the top of the MSC unit, together with sharp evaporites pinch out on evaporite‐free pre‐Messinian structural highs, suggest an extreme Messinian Stage 3 base level draw down in the Ionian Basin. Such evidence should be carefully evaluated in the light of Messinian and post‐Messinian vertical crustal movements in the area. The results of this study demonstrates the importance of extracting from seismic data the Messinian paleotopography, the paleomorphology and the detailed stratal architecture in the in order to advance in the understanding of the deep basins Messinian depositional environments. Highlights First description of a new type of deepwater Messinian salt giant in the Ionian Sea. First quantification of the Messinian salt volume in the Ionian Sea. New seismic evidence of erosional surfces and Lago Mare deposits in the deep Ionian Basin. Further evidence of sea level lowering during the Messinian Salinity Crisis. Evidence for a different, physically separated deepwater Messinian salt basins in the Mediterranean.

Description: Prochlorococcus and Synechococcus are the most abundant free-living photosynthetic microorganisms in the ocean. Uncultivated lineages of these picocyanobacteria also thrive in the dimly illuminated upper part of oxygen-deficient zones (ODZs), where an important portion of ocean nitrogen (N) loss takes place via denitrification and anaerobic ammonium oxidation. Recent metagenomic studies revealed that ODZ Prochlorococcus have the genetic potential for using different N forms, including nitrate and nitrite, uncommon N sources for Prochlorococcus, but common for Synechococcus. To determine which N sources ODZ picocyanobacteria are actually using in nature, the cellular N-15 natural abundance (delta N-15) and assimilation rates of different N compounds were determined using cell sorting by flow cytometry and mass spectrometry. The natural delta N-15 of the ODZ Prochlorococcus varied from -4.0 parts per thousand to 13.0 parts per thousand (n = 9), with 50% of the values in the range of -2.1-2.6 parts per thousand. While the highest values suggest nitrate use, most observations indicate the use of nitrite, ammonium, or a mixture of N sources. Meanwhile, incubation experiments revealed potential assimilation rates of ammonium and urea in the same order of magnitude as that expected for total N in several environments including ODZs, whereas rates of nitrite and nitrate assimilation were very low. Our results thus indicate that reduced forms of N and nitrite are the dominant sources for ODZ picocyanobacteria, although nitrate might be important on some occasions. ODZ picocyanobacteria might thus represent potential competitors with anammox bacteria for ammonium and nitrite, with ammonia-oxidizing archaea for ammonium, and with nitrite-oxidizing bacteria for nitrite.

Description: The immune system affects senescence (declines in probabilities of survival or reproduction with age), by shaping late age vulnerability to chronic inflammatory diseases and infections. It is also a dynamic interactive system that must balance competing demands across the life course. Thus, immune system function remains an important frontier in understanding the evolution of senescence. Here, we review our expanding mechanistic understanding of immune function over the life course, in the context of theoretical predictions from life-history evolution. We are especially interested in stage- and sex-dependent costs and benefits of investment in the immune system, given differential life-history priorities of the life stages and sexes. We introduce the costs likely to govern immune allocation across the life course. We then discuss theoretical expectations for differences between the sexes and their likely consequences in terms of how the immune system is both modulated by and may modulate senescence, building on information from life-history theory, experimental immunology and demography. We argue that sex differences in immune function provide a potentially powerful probe of selection pressures on the immune system across the life course. In particular, differences in 'competing' and 'caring' between the sexes have evolved across the tree of life, providing repeated instances of divergent selection pressures on immune function occurring within the same overall bauplan. We conclude by detailing an agenda for future research, including development of theoretical predictions of the differences between the sexes under an array of existing models for sex differences in immunity, and empirical tests of such predictions across the tree of life. A free Plain Language Summary can be found within the Supporting Information of this article.

Description: The plea for using more “realistic,” community‐level, investigations to assess the ecological impacts of global change has recently intensified. Such experiments are typically more complex, longer, more expensive, and harder to interpret than simple organism‐level benchtop experiments. Are they worth the extra effort? Using outdoor mesocosms, we investigated the effects of ocean warming (OW) and acidification (OA), their combination (OAW), and their natural fluctuations on coastal communities of the western Baltic Sea during all four seasons. These communities are dominated by the perennial and canopy‐forming macrophyte Fucus vesiculosus—an important ecosystem engineer Baltic‐wide. We, additionally, assessed the direct response of organisms to temperature and pH in benchtop experiments, and examined how well organism‐level responses can predict community‐level responses to the dominant driver, OW. OW affected the mesocosm communities substantially stronger than acidification. OW provoked structural and functional shifts in the community that differed in strength and direction among seasons. The organism‐level response to OW matched well the community‐level response of a given species only under warm and cold thermal stress, that is, in summer and winter. In other seasons, shifts in biotic interactions masked the direct OW effects. The combination of direct OW effects and OW‐driven shifts of biotic interactions is likely to jeopardize the future of the habitat‐forming macroalga F. vesiculosus in the Baltic Sea. Furthermore, we conclude that seasonal mesocosm experiments are essential for our understanding of global change impact because they take into account the important fluctuations of abiotic and biotic pressures.

Description: Over the past decades, coralline algae have increasingly been used as archives of palaeoclimate information due to their seasonal growth bands and their vast distribution from high latitudes to the tropics. Traditionally, these reconstructions have been performed mainly on high latitude species, limiting the geographical area of their potential use. Here we assess the use of temperate crustose fossil coralline algae from shallow water habitats for palaeoenvironmental reconstruction to generate records of past climate change. We determine the potential of three different species of coralline algae, Lithothamnion minervae, Lithophyllum stictaeforme and Mesophyllum philippii, with different growth patterns, as archives for pH (δ11B) and temperature (Mg/Ca) reconstruction in the Mediterranean Sea. Mg concentration is driven by temperature but modulated by growth rate, which is controlled by species‐specific and intraspecific growth patterns. L. minervae is a good temperature recorder, showing a moderate warming trend in specimens from 11.37 cal ka BP (from 14.2 ± 0.4°C to 14.9 ± 0.15°C) to today. In contrast to Mg, all genera showed consistent values of boron isotopes (δ11B) suggesting a common control on boron incorporation. The recorded δ11B in modern and fossil coralline specimens is in agreement with literature data about early Holocene pH, opening new perspectives of coralline‐based, high‐resolution pH reconstructions in deep time.

Description: We report on the geochemistry of hydrocarbons and pore waters down to 62.5 mbsf, collected by drilling with the MARUM‐MeBo70 and by gravity coring at the Lunde pockmark in the Vestnesa Ridge. Our data document the origin and transformations of volatiles feeding gas emissions previously documented in this region. Gas hydrates are present where a fracture network beneath the pockmark focusses migration of thermogenic hydrocarbons characterized by their C1/C2+ and stable isotopic compositions (δ2H‐CH4, δ13C‐CH4). Measured geothermal gradients (~80°C km‐1) and known formation temperatures (〉70°C) suggest that those hydrocarbons are formed at depths 〉800 mbsf. A combined analytical/modeling approach, including concentration and isotopic mass balances, reveals that pockmark sediments experience diffuse migration of thermogenic hydrocarbons. However, at sites without channeled flow this appears to be limited to depths 〉 ~50 mbsf. At all sites we document a contribution of microbial methanogenesis to the overall carbon cycle that includes a component of secondary carbonate reduction (CR) – i.e. reduction of dissolved inorganic carbon (DIC) generated by anaerobic oxidation of methane (AOM) in the uppermost methanogenic zone. AOM and CR rates are spatially variable within the pockmark and are highest at high‐flux sites. These reactions are revealed by δ13C‐DIC depletions at the sulfate‐methane interface at all sites. However, δ13C‐CH4 depletions are only observed at the low methane flux sites because changes in the isotopic composition of the overall methane pool are masked at high‐flux sites. 13C‐depletions of TOC suggest that at seeps sites, methane‐derived carbon is incorporated into de novo synthesized biomass.

Description: Sympatric speciation occurs without geographical barriers and is thought to often be driven by ecological specialization of individuals that eventually diverge genetically and phenotypically. Distinct morphologies between sympatric populations occupying different niches have been interpreted as such differentiating adaptive phenotypes, yet differences in performance and thus likely adaptiveness between them were rarely tested. Here, we investigated if divergent body shapes of two sympatric crater lake cichlid species from Nicaragua, one being a shore‐associated (benthic) species while the other prefers the open water zones (limnetic), affect cruising (Ucrit) and sprinting (Usprint) swimming abilities ‐ performances particularly relevant to their respective lifestyles. Furthermore, we investigated species differences in oxygen consumption (MO2) across different swimming speeds and compare gene expression in gills and white muscle at rest and during exercise. We found a superior cruising ability in the limnetic Amphilophus zaliosus compared to the benthic A. astorquii, while sprinting was not different, suggesting that their distinct morphologies affect swimming performance. Increased cruising swimming ability in A. zaliosus was linked to a higher oxygen demand during activity (but not rest), indicating different metabolic rates during exercise ‐ a hypothesis supported by coinciding gene expression patterns of gill transcriptomes. We identified differentially expressed genes linked to swimming physiology, regulation of swimming behaviour and oxygen intake. A combination of physiological and morphological differences may thus underlie adaptations to these species’ distinct niches. This complex ecological specialization likely resulted in morphological and physiological trade‐offs that contributed to the rapid establishment and maintenance of divergence with gene flow.

Description: The duration and magnitude of the North Atlantic spring bloom impacts both higher trophic levels and oceanic carbon sequestration. Nutrient exhaustion offers a general explanation for bloom termination, but detail on which nutrients and their relative influence on phytoplankton productivity, community structure, and physiology is lacking. Here, we address this using nutrient addition bioassay experiments conducted across the midlatitude North Atlantic in June 2017 (late spring). In four out of six experiments, phytoplankton accumulated over 48–72 h following individual additions of either iron (Fe) or nitrogen (N). In the remaining two experiments, Fe and N were serially limiting, that is, their combined addition sequentially enhanced phytoplankton accumulation. Silicic acid (Si) added in combination with N + Fe led to further chlorophyll a (Chl a) enhancement at two sites. Conversely, addition of zinc, manganese, cobalt, vitamin B12, or phosphate in combination with N + Fe did not. At two sites, the simultaneous supply of all six nutrients, in combination with N + Fe, also led to no further Chl a enhancement, but did result in an additional 30–60% particulate carbon accumulation. This particulate carbon accumulation was not matched by a Redfield equivalent of particulate N, characteristic of high C:N organic exudates that enhance cell aggregation and sinking. Our results suggest that growth rates of larger phytoplankton were primarily limited by Fe and/or N, making the availability of these nutrients the main bottom‐up factors contributing to spring bloom termination. In addition, the simultaneous availability of other nutrients could modify bloom characteristics and carbon export efficiency.

Description: Considerable uncertainty remains over how increasing atmospheric CO2 and anthropogenic climate changes are affecting open‐ocean marine ecosystems from phytoplankton to top predators. Biological time series data are thus urgently needed for the world's oceans. Here, we use the carbon stable isotope composition of tuna to provide a first insight into the existence of global trends in complex ecosystem dynamics and changes in the oceanic carbon cycle. From 2000 to 2015, considerable declines in δ13C values of 0.8‰–2.5‰ were observed across three tuna species sampled globally, with more substantial changes in the Pacific Ocean compared to the Atlantic and Indian Oceans. Tuna recorded not only the Suess effect, that is, fossil fuel‐derived and isotopically light carbon being incorporated into marine ecosystems, but also recorded profound changes at the base of marine food webs. We suggest a global shift in phytoplankton community structure, for example, a reduction in 13C‐rich phytoplankton such as diatoms, and/or a change in phytoplankton physiology during this period, although this does not rule out other concomitant changes at higher levels in the food webs. Our study establishes tuna δ13C values as a candidate essential ocean variable to assess complex ecosystem responses to climate change at regional to global scales and over decadal timescales. Finally, this time series will be invaluable in calibrating and validating global earth system models to project changes in marine biota.

Description: Key Points Calibration of XRF core scanning data highlights the need for careful examination of sediment properties such as porosity/water Grain size and water content in the sediment trigger systematic artifacts in the signal intensity of light elements (e.g. Si and Al) Known terrigenous flux proxies (e.g Ti/Ca, Fe/Ca) are influenced by sea level variations X‐ray fluorescence (XRF) core scanning of marine and lake sediments has been extensively used to study changes in past environmental and climatic processes over a range of timescales. The interpretation of XRF‐derived element ratios in paleoclimatic and paleoceanographic studies primarily considers differences in the relative abundances of particular elements. Here we present new XRF core scanning data from two long sediment cores in the Andaman Sea in the northern Indian Ocean and show that sea level related processes influence terrigenous inputs based proxies such as Ti/Ca, Fe/Ca, and elemental concentrations of the transition metals (e.g. Mn). Zr/Rb ratios are mainly a function of changes in median grain size of lithogenic particles and often covary with changes in Ca concentrations that reflect changes in biogenic calcium carbonate production. This suggests that a common process (i.e. sea level) influences both records. The interpretation of lighter element data (e.g. Si and Al) based on low XRF counts is complicated as variations in mean grain size and water content result in systematic artifacts and signal intensities not related to the Al or Si content of the sediments. This highlights the need for calibration of XRF core scanning data based on discrete sample analyses and careful examination of sediment properties such as porosity/water content for reliably disentangling environmental signals from other physical properties. In the case of the Andaman Sea, reliable extraction of a monsoon signal will require accounting for the sea level influence on the XRF data.

Description: Weddell Sea-derived Antarctic Bottom Water (AABW) is one of the most important deep water masses in the Southern Hemisphere occupying large portions of the deep Southern Ocean (SO) today. While substantial changes in SO-overturning circulation were previously suggested, the state of Weddell Sea AABW export during glacial climates remains poorly understood. Here we report seawater-derived Nd and Pb isotope records that provide evidence for the absence of Weddell Sea-derived AABW in the Atlantic sector of the SO during the last two glacial maxima. Increasing delivery of Antarctic Pb to regions outside the Weddell Sea traced SO frontal displacements during both glacial terminations. The export of Weddell Sea-derived AABW resumed late during glacial terminations, coinciding with the last major atmospheric CO2 rise in the transition to the Holocene and the Eemian. Our new records lend strong support for a previously inferred AABW overturning stagnation event during the peak Eemian interglacial.

Description: Idealized models or emulators of volcanic aerosol forcing have been widely used to reconstruct the spatiotemporal evolution of past volcanic forcing. However, existing models, including the most recently developed Easy Volcanic Aerosol (EVA; Toohey et al., doi: 10.5194/gmd‐2016‐83), (i) do not account for the height of injection of volcanic SO urn:x-wiley:jgrd:media:jgrd55987:jgrd55987-math-0001; (ii) prescribe a vertical structure for the forcing; and (iii) are often calibrated against a single eruption. We present a new idealized model, EVA_H, that addresses these limitations. Compared to EVA, EVA_H makes predictions of the global mean stratospheric aerosol optical depth that are (i) similar for the 1979–1998 period characterized by the large and high‐altitude tropical SO urn:x-wiley:jgrd:media:jgrd55987:jgrd55987-math-0002 injections of El Chichón (1982) and Mount Pinatubo (1991); (ii) significantly improved for the 1998–2015 period characterized by smaller eruptions with a large variety of injection latitudes and heights. Compared to EVA, the sensitivity of volcanic forcing to injection latitude and height in EVA_H is much more consistent with results from climate models that include interactive aerosol chemistry and microphysics, even though EVA_H remains less sensitive to eruption latitude than the latter models. We apply EVA_H to investigate potential biases and uncertainties in EVA‐based volcanic forcing data sets from phase 6 of the Coupled Model Intercomparison Project (CMIP6). EVA and EVA_H forcing reconstructions do not significantly differ for tropical high‐altitude volcanic injections. However, for high‐latitude or low‐altitude injections, our reconstructed forcing is significantly lower. This suggests that volcanic forcing in CMIP6 last millenium experiments may be overestimated for such eruptions.

Description: Along the margins of continental ice sheets, lakes formed in isostatically depressed basins duringglacial retreat. Their shorelines and extent are sensitive to the ice margin and the glacial history of the region.Proglacial lakes, in turn, also impact the glacial isostatic adjustment due to loading, and ice dynamics by posing amarine‐like boundary condition at the ice margin. In this study we present a tool that efficiently identifies lake basinsand the corresponding maximum water level for a given ice sheet and topography reconstruction. This algorithm,called the LakeCC model, iteratively checks the whole map for a set of increasing water levels and fills isolated basinsuntil they overflow into the ocean. We apply it to the present‐day Great Lakes and the results show good agreement(∼1−4%) with measured lake volume and depth. We then apply it to two topography reconstructions of NorthAmerica between the Last Glacial Maximum and the present. The model successfully reconstructs glacial lakes suchas Lake Agassiz, Lake McConnell and the predecessors of the Great Lakes. LakeCC can be used to judge the quality ofice sheet reconstructions.

Description: Species distribution models (SDMs) correlate species occurrences with environmental predictors, and can be used to forecast distributions under future climates. SDMs have been criticized for not explicitly including the physiological processes underlying the species response to the environment. Recently, new methods have been suggested to combine SDMs with physiological estimates of performance (physiology-SDMs). In this study, we compare SDM and physiology-SDM predictions for select marine species in the Mediterranean Sea, a region subjected to exceptionally rapid climate change. We focused on six species and created physiology-SDMs that incorporate physiological thermal performance curves from experimental data with species occurrence records. We then contrasted projections of SDMs and physiology-SDMs under future climate (year 2100) for the entire Mediterranean Sea, and particularly the ‘warm’ trailing edge in the Levant region. Across the Mediterranean, we found cross-validation model performance to be similar for regular SDMs and physiology-SDMs. However, we also show that for around half the species the physiology-SDMs substantially outperform regular SDM in the warm Levant. Moreover, for all species the uncertainty associated with the coefficients estimated from the physiology-SDMs were much lower than in the regular SDMs. Under future climate, we find that both SDMs and physiology-SDMs showed similar patterns, with species predicted to shift their distribution north-west in accordance with warming sea temperatures. However, for the physiology-SDMs predicted distributional changes are more moderate than those predicted by regular SDMs. We conclude, that while physiology-SDM predictions generally agree with the regular SDMs, incorporation of the physiological data led to less extreme range shift forecasts. The results suggest that climate-induced range shifts may be less drastic than previously predicted, and thus most species are unlikely to completely disappear with warming climate. Taken together, the findings emphasize that physiological experimental data can provide valuable supplemental information to predict range shifts of marine species.

Description: Back‐arc basins open in response to subduction processes, which cause extension in the upper plate, usually along trench‐parallel spreading axes. However, global seismic databases reveal that the majority of seismic events in the Lau Basin occur along transcurrent (strike‐slip) rather than extensional faults. To better characterize active deformation in this region we compared Centroid Moment Tensors (CMTs), calculated for large (Mw 〉5), shallow (〈30 km) seismic events to the orientations of seafloor lineaments mapped throughout the Lau Basin. Ship‐based multibeam was combined with vertical gravity gradient data to provide comprehensive coverage to create the lineament map. By comparing the possible focal planes of the CMTs to the orientations of the lineaments, the most likely fault plane solutions were selected, thus classifying the faults and establishing the nature of the highly variable stress regimes in the basin. We resolved the strike, dip and dip direction of 308 faults, and classified 258 additional structures by fault type. The analysis highlights a stress regime that is dominated by a combination of left‐lateral and right‐lateral strike‐slip faults, large‐scale transcurrent motion along rigid crustal‐scale fault zones, and non‐rigid diffuse deformation along pre‐existing seafloor structures, with extension mainly limited to the tips of propagating rifts and spreading centers. By resolving many of the uncertain motions on the mapped lineaments of the Lau Basin, the CMT analysis addresses a number of questions concerning basin‐scale stress regimes and microplate development, complementing GPS measurements and providing a more complete picture of the complexities of back‐arc basin development.

Description: Although offshore freshened groundwater (OFG) systems have been documented in numerous continental margins worldwide, their geometry, controls and emplacement dynamics remain poorly constrained. Here we integrate controlled-source electromagnetic, seismic reflection and borehole data with hydrological modelling to quantitatively characterise a previously unknown OFG system near Canterbury, New Zealand. The OFG system consists of one main, and two smaller, low salinity groundwater bodies. The main body extends up to 60 km from the coast and a seawater depth of 110 m. We attribute along-shelf variability in salinity to permeability heterogeneity due to permeable conduits and normal faults, and to recharge from rivers during sea level lowstands. A meteoric origin of the OFG and active groundwater migration from onshore are inferred. However, modelling results suggest that the majority of the OFG was emplaced via topographically-driven flow during sea level lowstands in the last 300 ka. Global volumetric estimates of OFG will be significantly revised if active margins, with steep coastal topographies like the Canterbury margin, are considered.

Description: The Middle Eocene Climatic Optimum (MECO) was a gradual warming event and carbon cycle perturbation that occurred between 40.5 and 40.1 Ma. A number of characteristics, including greater-than-expected deep-sea carbonate dissolution, a lack of globally-coherent negative δ 13 40 C excursion in marine carbonates, a duration longer than the characteristic timescale of carbon-cycle recovery, and the absence of a clear trigger mechanism, challenge our current understanding of the Earth system and its regulatory feedbacks. This makes the MECO one of the most enigmatic events in the Cenozoic, dubbed a middle Eocene ‘carbon cycle conundrum’. Here we use boron isotopes in planktic foraminifera to better constrain pCO2 changes over the event. Over the MECO itself, we find that pCO2 rose by only 0.55-0.75 doublings, thus requiring a much more modest carbon injection than previously indicated by the alkenone δ 13 C-pCO2 proxy. In addition, this rise in pCO2 was focused around the peak of the 400 kyr warming trend. Before this, considerable global carbonate δ 18 O change was asynchronous with any coherent ocean pH (and hence pCO2) excursion. This finding suggests that middle Eocene climate (and perhaps a nascent cryosphere) was highly sensitive to small changes in radiative forcing.

Description: All multicellular organisms are genetic mosaics owing to somatic mutations. The accumulation of somatic genetic variation in clonal species undergoing asexual (or clonal) reproduction may lead to phenotypic heterogeneity among autonomous modules (termed ramets). However, the abundance and dynamics of somatic genetic variation under clonal reproduction remain poorly understood. Here we show that branching events in a seagrass (Zostera marina) clone or genet lead to population bottlenecks of tissue that result in the evolution of genetically differentiated ramets in a process of somatic genetic drift. By studying inter-ramet somatic genetic variation, we uncovered thousands of single nucleotide polymorphisms that segregated among ramets. Ultra-deep resequencing of single ramets revealed that the strength of purifying selection on mosaic genetic variation was greater within than among ramets. Our study provides evidence for multiple levels of selection during the evolution of seagrass genets. Somatic genetic drift during clonal propagation leads to the emergence of genetically unique modules that constitute an elementary level of selection and individuality in long-lived clonal species.

Description: Predation is a pervasive force that structures food webs and directly influences ecosystem functioning. The relative body sizes of predators and prey may be an important determinant of interaction strengths. However, studies quantifying the combined influence of intra‐ and interspecific variation in predator–prey body size ratios are lacking. We use a comparative functional response approach to examine interaction strengths between three size classes of invasive bluegill and largemouth bass toward three scaled size classes of their tilapia prey. We then quantify the influence of intra‐ and interspecific predator–prey body mass ratios on the scaling of attack rates and handling times. Type II functional responses were displayed by both predators across all predator and prey size classes. Largemouth bass consumed more than bluegill at small and intermediate predator size classes, while large predators of both species were more similar. Small prey were most vulnerable overall; however, differential attack rates among prey were emergent across predator sizes. For both bluegill and largemouth bass, small predators exhibited higher attack rates toward small and intermediate prey sizes, while larger predators exhibited greater attack rates toward large prey. Conversely, handling times increased with prey size, with small bluegill exhibiting particularly low feeding rates toward medium–large prey types. Attack rates for both predators peaked unimodally at intermediate predator–prey body mass ratios, while handling times generally shortened across increasing body mass ratios. We thus demonstrate effects of body size ratios on predator–prey interaction strengths between key fish species, with attack rates and handling times dependent on the relative sizes of predator–prey participants. Considerations for intra‐ and interspecific body size ratio effects are critical for predicting the strengths of interactions within ecosystems and may drive differential ecological impacts among invasive species as size ratios shift.

Description: One key contribution to the wide range of 1.5 degrees C carbon budgets among recent studies is the non-CO2 climate forcing scenario uncertainty. Based on a partitioning of historical non-CO2 forcing, we show that currently there is a net negative non-CO2 forcing from fossil fuel combustion (FFC), and a net positive non-CO2 climate forcing from land-use change (LUC) and agricultural activities. We perform a set of future simulations in which we prescribed a 1.5 degrees C temperature stabilisation trajectory, and diagnosed the resulting 1.5 degrees C carbon budgets. Using the historical partitioning, we then prescribed adjusted non-CO2 forcing scenarios consistent with our model's simulated decrease in FFC CO2 emissions. We compared the diagnosed carbon budgets from these adjusted scenarios to those resulting from the default RCP scenario's non-CO2 forcing, and to a scenario in which proportionality between future CO2 and non-CO2 forcing is assumed. We find a wide range of carbon budget estimates across scenarios, with the largest budget emerging from the scenario with assumed proportionality of CO2 and non-CO2 forcing. Furthermore, our adjusted-RCP scenarios produce carbon budgets that are smaller than the corresponding default RCP scenarios. Our results suggest that ambitious mitigation scenarios will likely be characterised by an increasing contribution of non-CO2 forcing, and that an assumption of continued proportionality between CO2 and non-CO2 forcing would lead to an overestimate of the remaining carbon budget. Maintaining such proportionality under ambitious fossil fuel mitigation would require mitigation of non-CO2 emissions at a rate that is substantially faster than found in the standard RCP scenarios.

Description: Key points:  First insights into the crustal structure of the northeastern Lau Basin, along a 290 km transect at 17°20’S.  Crust in southern Fonualei Rift and Spreading Center was created by extension of arc crust and variable amount of magmatism.  Magmatic underplating is present in some parts of the southern Niuafo’ou Microplate The northeastern Lau Basin is one of the fastest opening and magmatically most active back‐arc regions on Earth. Although the current pattern of plate boundaries and motions in this complex mosaic of microplates is reasonably understood, the internal structure and evolution of the back‐arc crust are not. We present new geophysical data from a 290 km long east‐west oriented transect crossing the Niuafo’ou Microplate (back‐arc), the Fonualei Rift and Spreading Centre (FRSC) and the Tofua Volcanic Arc at 17°20’S. Our P‐wave tomography model and density modelling suggests that past crustal accretion inside the southern FRSC was accommodated by a combination of arc crustal extension and magmatic activity. The absence of magnetic reversals inside the FRSC supports this and suggests that focused seafloor spreading has until now not contributed to crustal accretion. The back‐arc crust constituting the southern Niuafo’ou Microplate reveals a heterogeneous structure comprising several crustal blocks. Some regions of the back‐arc show a crustal structure similar to typical oceanic crust, suggesting they originate from seafloor spreading. Other crustal blocks resemble a structure that is similar to volcanic arc crust or a ‘hydrous’ type of oceanic crust that has been created at a spreading center influenced by slab‐derived water at distances 〈 50 km to the arc. Throughout the back‐arc region we observe a high‐velocity (Vp 7.2‐7.5 km s‐1) lower crust, which is an indication for magmatic underplating, which is likely sustained by elevated upper mantle temperatures in this region.

Description: The glass sponge Aphrocallistes vastus contributes to the formation of large reefs unique to the Northeast Pacific Ocean. These habitats have tremendous filtration capacity that facilitates flow of carbon between trophic levels. Their sensitivity and resilience to climate change, and thus persistence in the Anthropocene, is unknown. Here we show that ocean acidification and warming, alone and in combination have significant adverse effects on pumping capacity, contribute to irreversible tissue withdrawal, and weaken skeletal strength and stiffness of A. vastus. Within one month sponges exposed to warming (including combined treatment) ceased pumping (50–60%) and exhibited tissue withdrawal (10–25%). Thermal and acidification stress significantly reduced skeletal stiffness, and warming weakened it, potentially curtailing reef formation. Environmental data suggests conditions causing irreversible damage are possible in the field at +0.5 °C above current conditions, indicating that ongoing climate change is a serious and immediate threat to A. vastus, reef dependent communities, and potentially other glass sponges.

Description: High‐temperature hydrothermal venting has been discovered on all modern mid‐ocean ridges at all spreading rates. Although significant strides have been made in understanding the underlying processes that shape such systems, several first‐order discrepancies between model predictions and observations remain. One key paradox is that numerical experiments consistently show entrainment of cold ambient seawater in shallow high permeability ocean crust causing a temperature drop that is difficult to reconcile with high vent temperatures. We investigate this conundrum using a thermo‐hydro‐chemical model that couples hydrothermal fluid flow with anhydrite‐ and pyrite‐forming reactions in the shallow subseafloor. The models show that precipitation of anhydrite in warming seawater and in cooling hydrothermal fluids during mixing results in the formation of a chimney‐like subseafloor structure around the upwelling, high‐temperature plume. The establishment of such anhydrite‐sealed zones reduces mixing between the hydrothermal fluid and seawater and results in an increase in vent temperature. Pyrite subsequently precipitates close to the seafloor within the anhydrite chimney. Although anhydrite thus formed may be dissolved when colder seawater circulates through the crust away from the spreading axis, the inside pyrite walls would be preserved as veins in present‐day metal deposits, thereby preserving the history of hydrothermal circulation through shallow oceanic crust.

Description: Sponges are the oldest known extant animal-microbe symbiosis. These ubiquitous benthic animals play an important role in marine ecosystems in the cycling of dissolved organic matter (DOM), the largest source of organic matter on Earth. The conventional view on DOM cycling through microbial processing has been challenged by the interaction between this efficient filter-feeding host and its diverse and abundant microbiome. Here we quantify, for the first time, the role of host cells and microbial symbionts in sponge heterotrophy. We combined stable isotope probing and nanoscale secondary ion mass spectrometry to compare the processing of different sources of DOM (glucose, amino acids, algal-produced) and particulate organic matter (POM) by a high-microbial abundance (HMA) and low-microbial abundance (LMA) sponge with single-cell resolution. Contrary to common notion, we found that both microbial symbionts and host choanocyte (i.e. filter) cells and were active in DOM uptake. Although all DOM sources were assimilated by both sponges, higher microbial biomass in the HMA sponge corresponded to an increased capacity to process a greater variety of dissolved compounds. Nevertheless, in situ feeding data demonstrated that DOM was the primary carbon source for both the LMA and HMA sponge, accounting for ~90% of their heterotrophic diets. Microbes accounted for the majority (65–87%) of DOM assimilated by the HMA sponge (and ~60% of its total heterotrophic diet) but 〈5% in the LMA sponge. We propose that the evolutionary success of sponges is due to their different strategies to exploit the vast reservoir of DOM in the ocean.

Description: The time‐domain controlled source electromagnetic method is a geophysical prospecting tool applied to image the subsurface resistivity distribution on land and in the marine environment. In its most general setup, a square‐wave current is fed into a grounded horizontal electric dipole, and several electric and magnetic field receivers at defined offsets to the imposed current measure the electromagnetic response of the Earth. In the marine environment, the application often uses only inline electric field receivers that, for a 50% duty‐cycle current waveform, include both step‐on and step‐off signals. Here, forward and inverse 1D modelling is used to demonstrate limited sensitivity towards shallow resistive layers in the step‐off electric field when transmitter and receivers are surrounded by conductive seawater. This observation is explained by a masking effect of the direct current signal that flows through the seawater and primarily affects step‐off data. During a step‐off measurement, this direct current is orders of magnitude larger than the inductive response at early and intermediate times, limiting the step‐off sensitivity towards shallow resistive layers in the seafloor. Step‐on data measure the resistive layer at times preceding the arrival of the direct current signal leading to higher sensitivity compared to step‐off data. Such dichotomous behaviour between step‐on and step‐off data is less obvious in onshore experiments due to the lack of a strong overlying conductive zone and corresponding masking effect from direct current flow. Supported by synthetic 1D inversion studies, we conclude that time‐domain controlled source electromagnetic measurements on land should apply both step‐on and step‐off data in a combined inversion approach to maximise signal‐to‐noise ratios and utilise the sensitivity characteristics of each signal. In an isotropic marine environment, step‐off electric fields have inferior sensitivity towards shallow resistive layers compared to step‐on data, resulting in an increase of non‐uniqueness when interpreting step‐off data in a single or combined inversion.

Description: Marine food webs are highly compartmentalized, and characterizing the trophic niches among consumers is important for predicting how impact from human activities affects the structuring and functioning of marine food webs. Biomarkers such as bulk stable isotopes have proven to be powerful tools to elucidate trophic niches, but they may lack in resolution, particularly when spatiotemporal variability in a system is high. To close this gap, we investigated whether carbon isotope (delta C-13) patterns of essential amino acids (EAAs), also termed delta(13)C(AA)fingerprints, can characterize niche differentiation in a highly dynamic marine system. Specifically, we tested the ability of delta(13)C(AA)fingerprints to differentiate trophic niches among six functional groups and ten individual species in the Baltic Sea. We also tested whether fingerprints of the common zooplanktivorous fishes, herring and sprat, differ among four Baltic Sea regions with different biochemical conditions and phytoplankton assemblages. Additionally, we investigated how these results compared to bulk C and N isotope data for the same sample set. We found significantly different delta(13)C(AA)fingerprints among all six functional groups. Species differentiation was in comparison less distinct, due to partial convergence of the species' fingerprints within functional groups. Herring and sprat displayed region-specific delta(13)C(AA)fingerprints indicating that this approach could be used as a migratory marker. Niche metrics analyses showed that bulk isotope data had a lower power to differentiate between trophic niches than delta(13)C(AA)fingerprinting. We conclude that delta(13)C(AA)fingerprinting has a strong potential to advance our understanding of ecological niches, and trophic linkages from producers to higher trophic levels in dynamic marine systems. Given how management practices of marine resources and habitats are reshaping the structure and function of marine food webs, implementing new and powerful tracer methods are urgently needed to improve the knowledge base for policy makers.

Description: Ecosystem connectivity is an essential consideration for marine spatial planning of competing interests in the deep sea. Immobile, adult communities are connected through freely floating larvae, depending on new recruits for their health and to adapt to external pressures. We hypothesize that the vertical swimming ability of deep-sea larvae, before they permanently settle at the bottom, is one way larvae can control dispersal. We test this hypothesis with more than 3x108 simulated particles with a range of active swimming behaviours embedded within the currents of a high-resolution ocean model. Despite much stronger horizontal ocean currents, vertical swimming of simulated larvae can have an order of magnitude impact on dispersal. These strong relationships between larval dispersal, pathways, and active swimming demonstrate that lack of data on larval behaviour traits is a serious impediment to modelling deep-sea ecosystem connectivity; this uncertainty greatly limits our ability to develop ecologically coherent marine protected area networks.

Description: Key Points: - Fluid flow is focused along Nootka Fault traces resulting in shallow bright spots - Two seafloor mounds are the result of basaltic intrusions in the Nootka Fault zone - Gas hydrates occur at the Nootka Slope and are imaged seismically as bottom- simulating reflectors suggesting a regional heat-flow of ~80 mW/m2 along the slope Abstract Geophysical and geochemical data indicate there is abundant fluid expulsion in the Nootka fault zone (NFZ) between the Juan de Fuca and Explorer plates and the Nootka continental slope. Here we combine observations from 〉 20 years of investigations to demonstrate the nature of fluid‐flow along the NFZ, which is the seismically most active region off Vancouver Island. Seismicity reaching down to the upper mantle is linked to near‐seafloor manifestation of fluid flow through a network of faults. Along the two main fault traces, seismic reflection data imaged bright spots 100 – 300 m below seafloor that lie above changes in basement topography. The bright spots are conformable to sediment layering, show opposite‐to‐seafloor reflection polarity, and are associated with frequency‐reduction and velocity push‐down indicating the presence of gas in the sediments. Two seafloor mounds ~15 km seaward of the Nootka slope are underlain by deep, non‐conformable high amplitude reflective zones. Measurements in the water column above one mound revealed a plume of warm water, and bottom‐video observations imaged hydrothermal vent system biota. Pore fluids from a core at this mound contain predominately microbial methane (C1) with a high proportion of ethane (C2) yielding C1/C2 ratios 〈 500 indicating a possible slight contribution from a deep source. We infer the reflective zones beneath the two mounds are basaltic intrusions that create hydrothermal circulation within the overlying sediments. Across the Nootka continental slope, gas hydrate related bottom‐simulating reflectors are widespread and occur at depths indicating heat‐flow values of 80 – 90 mW/m2.

Description: The distribution of dissolved silicon isotopes (δ30Si) was examined along the US GEOTRACES East Pacific Zonal Transect (GP16) extending from Peru to Tahiti (10°S and 15°S latitude). Surface waters in the subtropical gyre displayed high δ30Si due to strong utilization of silicic acid (DSi). In contrast, surface waters close to the Peruvian coast where upwelling prevailed were less depleted and only moderately fractionated. δ30Si of water masses along the transect were compared with the results of an Optimum Multiparameter Analysis that quantified the fractional contributions of endmember water masses in each sample. Strong admixture of intermediate waters obscured the expected heavy isotopic signatures of Subantarctic Mode Water and Antarctic Intermediate Water. Isotope values were nearly homogenous below 2000 m (Average: +1.3 ± 0.1 ‰, 1 s.d.) despite the 25 μmol kg‐1 range in the DSi content among water masses. This homogeneity confirms prior observations and model results that predict nearly constant δ30Si values of +1.0 to +1.2 ‰ for Pacific deep waters with [DSi] 〉 100 μmol kg‐1. Waters above the East Pacific Rise (EPR) influenced by hydrothermal activity showed a small increase in [DSi] together with dissolved iron, but overall stations close to the EPR were slightly depleted in [DSi] (3 to 6 μmol kg‐1) with no significant shift in δ30Si compared to adjacent waters. Hydrothermal [DSi] appears to precipitate within the conduit of the EPR or upon contact with cold seawater resulting in a negligible influence of hydrothermal fluids on δ30Si in deep water. Key Points Surface waters have a large range in dissolved silicon isotopes covering nutrient‐rich coastal upwelling to oligotrophic waters Deep water masses with DSi concentrations 〉 100 μmol kg‐1 show homogenous silicon isotope signatures despite up to 25 μmol kg‐1 differences in [DSi] Hydrothermal fluids have a negligible effect on Si isotope distributions in the deep Pacific

Description: Reconstructions of the global mean annual temperature evolution during the Holocene yield conflicting results. One temperature reconstruction shows global cooling during the late Holocene. The other reconstruction reveals global warming. Here we show that both a global warming mode and a cooling mode emerge when performing a spatio-temporal analysis of annual temperature variability during the Holocene using data from a transient climate model simulation. The warming mode is most pronounced in the tropics. The simulated cooling mode is determined by changes in the seasonal cycle of Arctic sea-ice that are forced by orbital variations and volcanic eruptions. The warming mode dominates in the mid-Holocene, whereas the cooling mode takes over in the late Holocene. The weighted sum of the two modes yields the simulated global temperature trend evolution. Our findings have strong implications for the interpretation of proxy data and the selection of proxy locations to compute global mean temperatures.

Description: Serpentine alteration is recognized as an important process for element cycling, however, related silicon fluxes are unknown. Pore fluids from serpentinite seamounts sampled in the Mariana forearc region during IODP Expedition 366 were investigated for their Si, B, and Sr isotope signatures (δ 30 Si, δ 11 B, and 87 Sr/ 86 Sr, respectively) to study serpentinization in the mantle wedge and shallow serpentine alteration to authigenic clays by seawater. While serpentinization in the mantle wedge caused no significant Si isotope fractionation, implying closed system conditions, serpentine alteration by seawater led to the formation of authigenic phyllosilicates, causing the highest natural fluid δ 30 Si values measured to date (up to +5.2 ± 0.2‰). Here we show that seafloor alteration of serpentinites is a source of Si to the ocean with extremely high fluid δ 30 Si values, which can explain anomalies in the marine Si budget like in the Cascadia Basin and which has to be considered in future investigations of the global marine Si cycle.

Description: Every year, about four percent of the plastic waste generated worldwide ends up in the ocean. What happens to the plastic there is poorly understood, though a growing body of evidence suggests it is rapidly spreading throughout the global ocean. The mechanisms of this spread are straightforward for buoyant larger plastics that can be accurately modelled using Lagrangian particle models. But the fate of the smallest size fractions (the microplastics) are less straightforward, in part because they can aggregate in sinking marine snow and faecal pellets. This biologically-mediated pathway is suspected to be a primary surface microplastic removal mechanism, but exactly how it might work in the real ocean is unknown. We search the parameter space of a new microplastic model embedded in an earth system model to show that biological uptake can significantly shape global microplastic inventory and distributions and even account for the budgetary “missing” fraction of surface microplastic, despite being an inefficient removal mechanism. While a lack of observational data hampers our ability to choose a set of “best” model parameters, our effort represents a first tool for quantitatively assessing hypotheses for microplastic interaction with ocean biology at the global scale.

Description: What controls subduction megathrust seismogenesis downdip of the mantle wedge corner (MWC)? We propose that, in the region of the 2010 Mw=8.8 Maule, Chile, earthquake, serpentine minerals derived from the base of the hydrated mantle wedge exert a dominant control. Based on modeling, we predict that the megathrust fault zone near the MWC contains abundant lizardite/chrysotile‐rich serpentinite that transforms to antigorite‐rich serpentinite at greater depths. From the MWC at 32–40 km depth to at least 55 km, the predominantly velocity‐strengthening megathrust accommodated dynamic propagation of the 2010 rupture but with small slip and negative stress drop. The downdip distribution of interplate aftershocks exhibits a gap around the MWC that can be explained by the velocity‐strengthening behavior of lizardite/chrysotile. Interspersed velocity‐weakening and dynamic weakening antigorite‐rich patches farther downdip may be responsible for increased abundance of aftershocks and possibly for some of the high‐frequency energy radiation during the 2010 rupture.

Description: Key Points: • Observational transport time series of the Atlantic Subtropical Cells reveals dominant seasonal variability for horizontal branches • On time scales longer than ~5 years, interior thermocline layer transport convergence modulates equatorial sea surface temperature anomalies • Western boundary current and interior transport anomalies are partly compensating each other at thermocline level on all time scales The shallow meridional overturning cells of the Atlantic Ocean, the subtropical cells (STCs), consist of poleward Ekman transport at the surface, subduction in the subtropics, equatorward flow at thermocline level and upwelling along the equator and at the eastern boundary. In this study, we provide the first observational estimate of transport variability associated with the horizontal branches of the Atlantic STCs in both hemispheres based on Argo float data and supplemented by reanalysis products. Thermocline layer transport convergence and surface layer transport divergence between 10°N and 10°S are dominated by seasonal variability. Meridional thermocline layer transport anomalies at the western boundary and in the interior basin are anti‐correlated and partially compensate each other at all resolved time scales. It is suggested that the seesaw‐like relation is forced by the large‐scale off‐equatorial wind stress changes through low‐baroclinic‐mode Rossby wave adjustment. We further show that anomalies of the thermocline layer interior transport convergence modulate sea surface temperature (SST) variability in the upwelling regions along the equator and at the eastern boundary at time scales longer than 5 years. Phases of weaker (stronger) interior transport are associated with phases of higher (lower) equatorial SST. At these time scales, STC transport variability is forced by off‐equatorial wind stress changes, especially by those in the southern hemisphere. At shorter time scales, equatorial SST anomalies are, instead, mainly forced by local changes of zonal wind stress.

Description: Stable oxygen isotope records from central Greenland suggest disproportionally large long‐term surface air temperature (SAT) variability during the Last Glacial Maximum (LGM) relative to preindustrial times. Large perturbations in mean atmospheric circulation and its variability forced by extensive Northern Hemisphere ice sheet coverage have been suggested as cause for the enhanced Greenland SAT variability. Here, we assess the factors driving Greenland SAT variability during the LGM by means of dedicated climate model simulations and find remote forcing from the Pacific of critical importance. Atmospheric teleconnections from the Interdecadal Pacific Oscillation (IPO), a multidecadal oscillation of sea surface temperature in the Pacific Ocean, strongly intensify under LGM conditions, driving enhanced surface wind variability over Greenland, which in turn amplifies SAT variability by anomalous atmospheric heat transport. A major role of the IPO in forcing Greenland SAT variability also is supported by a number of models from the Paleoclimate Modeling Intercomparison Project Phase III

Description: Gravitationally consistent solutions of the Sea Level Equation from leakage‐corrected monthly‐mean GFZ RL06 Gravity Recovery and Climate Experiment (GRACE) and GRACE Follow‐On (GRACE‐FO) Stokes coefficients reveal that barystatic sea level averaged over the whole global ocean was rising by 1.72 mm a−1 during the period April 2002 until August 2016. This rate refers to a truely global ocean averaging domain that includes all polar and semienclosed seas. The result corresponds to 2.02 mm a−1 mean barystatic sea level rise in the open ocean with a 1,000 km coastal buffer zone as obtained from a direct spatial integration of monthly GRACE data. The bias of +0.3 mm a−1 is caused by below‐average barystatic sea level rise in close proximity to coastal mass losses induced by the smaller gravitational attraction of the remaining continental ice and water masses. Alternative spherical harmonics solutions from CSR, JPL, and TU Graz reveal open‐ocean rates between 1.94 and 2.08 mm a−1, thereby demonstrating that systematic differences among the processing centers are much reduced in the latest release. We introduce in this paper a new method to approximate spatial leakage from the differences of two differently filtered global gravity fields. A globally constant and time‐invariant scale factor required to obtain full leakage from those filter differences is found to be 3.9 for GFZ RL06 when filtered with DDK3, and lies between 3.9 and 4.4 for other processing centers. Spatial leakage is estimated for every month in terms of global grids, thereby providing also valuable information of intrabasin leakage that is potentially relevant for hydrologic and hydrometeorological applications.

Description: We use observations from novel biogeochemical profiling floats deployed by the Southern Ocean Carbon and Climate Observations and Modeling program to estimate annual net community production (ANCP; associated with carbon export) from the seasonal drawdown of mesopelagic oxygen and surface nitrate in the Southern Ocean. Our estimates agree with previous observations in showing an increase in ANCP in the vicinity of the polar front (∼3 mol C m−2 y−1), compared to lower rates in the subtropical zone (≤ 1 mol C m−2 y−1) and the seasonal ice zone (〈2 mol C m−2 y−1). Paradoxically, the increase in ANCP south of the subtropical front is associated with elevated surface nitrate and silicate concentrations, but decreasing surface iron. We hypothesize that iron limitation promotes silicification in diatoms, which is evidenced by the low silicate to nitrate ratio of surface waters around the Antarctic polar front. High diatom silicification increases the ballasting effect of particulate organic carbon and overall ANCP in this region. A model-based assessment of our methods shows a good agreement between ANCP estimates based on oxygen and nitrate drawdown and the modeled downward organic carbon flux at 100 m. This agreement supports the presumption that net biological consumption is the dominant process affecting the drawdown of these chemical tracers and that, given sufficient data, ANCP can be inferred from observations of oxygen and/or nitrate drawdown in the Southern Ocean.

Description: This paper presents the first detailed study of a late Pleistocene marine tephra sequence from the NW Pacific, downwind from the Kamchatka volcanic arc. Sediment core SO201-2-40, located on the Meiji Rise similar to 400 km offshore the peninsula, includes 25 tephras deposited within the last 215 ka. Volcanic glass from the tephras was characterized using single-shard electron microprobe analysis and laser ablation inductively coupled mass spectrometry. The age of tephras was derived from a new age model based on paleomagnetic and paleoclimate studies. Geochemical correlation of distal tephras to Kamchatkan pyroclastic deposits allowed the identification of tephras from the Karymsky, Gorely, Opala and Shiveluch eruptive centers. Three of these tephras were also correlated to other marine and terrestrial sites and hence are identified as the best markers for the north-west Pacific region. These are an early Holocene tephra from the Karymsky caldera (similar to 8.7 ka) and two tephras falling into the Marine Isotope Stage (MIS) 6 glacial time: an MIS 6.4 tephra from Shiveluch (similar to 141 ka) and the MIS 6.5 Rauchua tephra (similar to 175 ka) from Karymsky. The data presented in this study can be used in paleovolcanological and paleoceanographic reconstructions.

Description: Ecosystems are controlled by 'bottom-up' (resources) and 'top-down' (predation) forces. Viral infection is now recognized as a ubiquitous top-down control of microbial growth across ecosystems but, at the same time, cell death byviral predation influences, and is influenced by, resource availability. In this Review, we discuss recent advances in understanding the biogeochemical impact of viruses, focusing on how metabolic reprogramming of host cells during lytic viral infection alters the flow of energy and nutrients in aquatic ecosystems. Our synthesis revealed several emerging themes. First, viral infection transforms host metabolism, in part through virus-encoded metabolic genes; the functions performed by these genes appear to alleviate energetic and biosynthetic bottlenecks to viral production. Second, viral infection depends on the physiological state of the host cell and on environmental conditions, which are challenging to replicate in the laboratory. Last, metabolic reprogramming of infected cells and viral lysis alter nutrient cycling and carbon export in the oceans, although the net impacts remain uncertain. This Review highlights the need for understanding viral infection dynamics in realistic physiological and environmental contexts to better predict their biogeochemical consequences.

Description: Ballast water has been identified as a leading vector for introduction of non‐indigenous species. Recently, the International Maritime Organization implemented management standards—D‐2—where all large, commercial ships trading internationally are required to adopt an approved treatment system using technologies such as ultraviolet radiation or chlorination. However, current management regulations are based only on the total abundance of viable taxa transported (i.e. total propagule pressure), largely ignoring species richness (i.e. colonization pressure). To determine the efficacy of chlorine treatment in reducing invasion risks and changes in transported biological communities inside ballast tanks, we used DNA metabarcoding‐based approaches to estimate colonization pressure (here, the number of species/operational taxonomic units [OTUs] introduced) and relative propagule pressure (relative abundance of each species/OTU) of zooplankton communities in control and chlorine treated tanks during four transatlantic voyages. Our study demonstrated that transport itself did not significantly reduce colonization pressure of zooplankton species, nor did chlorine treatment. Chlorine treatment altered community structure by reducing relative propagule pressure of some taxa such as Mollusca and Rotifera, while increasing relative propagule pressure of some Oligohymenophorea and Copepoda species. Synthesis and applications. Chlorine treatment may not reduce invasion risks as much as previously thought. Reduction in total propagule pressure does not mean reduction in abundance of all species equally. While some taxa might experience drastically reduced abundance, others might not change at all or increase due to hatching from dormant stages initiated by chlorine exposure. Therefore, management strategies should consider changes in total propagule pressure and colonization pressure when forecasting risk of new invasions. We therefore recommend adopting new approaches, such as DNA metabarcoding‐based methods, to assess the whole biodiversity discharged from ballast water. As species responses to chlorine treatment are variable and affected by concentration, we also recommend a combination of different technologies to reduce introduction risks of aquatic organisms.

Description: Facilitated by the intensification of global trading, the introduction and dispersal of species to areas in which they are historically non-native is nowadays common. From an evolutionary standpoint, invasions are paradoxical: not only non-native environments could be different from native ones for which introduced individuals would be ill-adapted, but also small founding population size should be associated with reduced adaptive potential. As such, biological invasions are considered valuable real-time evolutionary experiments. Here, we investigated the population structure and adaptive potential of the highly invasive topmouth gudgeon (Pseudorasbora parva) across Europe and East Asia. We RAD-sequenced 301 specimens from sixteen populations and three distinct within-catchment invaded regions as well as two locations in the native range. With 13,785 single nucleotide polymorphisms, we provide conclusive evidence for a genome-wide signature of two distinct invasion events, in Slovakia and Turkey, each originating from a specific area in the native range. A third invaded area, in France, appears to be the result of dispersal within the invasive range. Few loci showed signs of selection, the vast majority of which being identified in the Slovakian region. Functional annotation suggests that faster early stage development, resistance to pollution and immunocompetence contribute to the invasion success of the local habitats. By showing that populations in the invasive range have different evolutionary histories, our study reinforces the idea that populations, rather than species, are the units to consider in invasion biology.

Description: Indian Summer Monsoon (ISM) rainfall has a direct effect on the livelihoods of two billion people in the Indian-subcontinent. Yet, our understanding of the drivers of multi-decadal variability of the ISM is far from being complete. In this context, large-scale forcing of ISM rainfall variability with multi-decadal resolution over the last two millennia is investigated using new records of sea surface salinity (δ18Ow) and sea surface temperatures (SSTs) from the Bay of Bengal (BoB). Higher δ18Ow values during the Dark Age Cold Period (1550 to 1250 years BP) and the Little Ice Age (700 to 200 years BP) are suggestive of reduced ISM rainfall, whereas lower δ18Ow values during the Medieval Warm Period (1200 to 800 years BP) and the major portion of the Roman Warm Period (1950 to 1550 years BP) indicate a wetter ISM. This variability in ISM rainfall appears to be modulated by the Atlantic Multi-decadal Oscillation (AMO) via changes in large-scale thermal contrast between the Asian land mass and the Indian Ocean, a relationship that is also identifiable in the observational data of the last century. Therefore, we suggest that inter-hemispheric scale interactions between such extra tropical forcing mechanisms and global warming are likely to be influential in determining future trends in ISM rainfall.

Description: Rivers are a major supplier of particulate and dissolved material to the ocean, but their role as sources of bio-essential dissolved iron (dFe) is thought to be limited due to rapid, efficient Fe removal during estuarine mixing. Here, we use trace element and radium isotope data to show that the influence of the Congo River margin on surface Fe concentrations is evident over 1000 km from the Congo outflow. Due to an unusual combination of high Fe input into the Congo-shelf-zone and rapid lateral transport, the Congo plume constitutes an exceptionally large offshore dFe flux of 6.8 ± 2.3 × 108 mol year−1. This corresponds to 40 ± 15% of atmospheric dFe input into the South Atlantic Ocean and makes a higher contribution to offshore Fe availability than any other river globally. The Congo River therefore contributes significantly to relieving Fe limitation of phytoplankton growth across much of the South Atlantic.

Description: The Atlantic Ocean overturning circulation is important to the climate system because it carries heat and carbon northward, and from the surface to the deep ocean. The high salinity of the subpolar North Atlantic is a prerequisite for overturning circulation, and strong freshening could herald a slowdown. We show that the eastern subpolar North Atlantic underwent extreme freshening during 2012 to 2016, with a magnitude never seen before in 120 years of measurements. The cause was unusual winter wind patterns driving major changes in ocean circulation, including slowing of the North Atlantic Current and diversion of Arctic freshwater from the western boundary into the eastern basins. We find that wind-driven routing of Arctic-origin freshwater intimately links conditions on the North West Atlantic shelf and slope region with the eastern subpolar basins. This reveals the importance of atmospheric forcing of intra-basin circulation in determining the salinity of the subpolar North Atlantic.

Description: Numerical models are a suitable tool to quantify impacts of predicted climate change on complex ecosystems but are rarely used to study effects on benthic macroalgal communities. Fucus vesiculosus L. is a habitat-forming macroalga in the Baltic Sea and alarming shifts from the perennial Fucus community to annual filamentous algae are reported. We developed a box model able to simulate the seasonal growth of the Baltic Fucus-grazer-epiphyte system. This required the implementation of two state variables for Fucus biomass in units of carbon (C) and nitrogen (N). Model equations describe relevant physiological and ecological processes, such as storage of C and N assimilates by Fucus, shading effects of epiphytes or grazing by herbivores on both Fucus and epiphytes, but with species-specific rates and preferences. Parametrizations of the model equations and required initial conditions were based on measured parameters and process rates in the near-natural Kiel Outdoor Benthocosm (KOB) experiments during the Biological Impacts of Ocean Acidification project. To validate the model, we compared simulation results with observations in the KOB experiment that lasted from April 2013 until March 2014 under ambient and climate-change scenarios, that is, increased atmospheric temperature and partial pressure of carbon dioxide. The model reproduced the magnitude and seasonal cycles of Fucus growth and other processes in the KOBs over 1 yr under different scenarios. Now having established the Fucus model, it will be possible to better highlight the actual threat of climate change to the Fucus community in the shallow nearshore waters of the Baltic Sea.

Description: Rationale To detect the small changes in past pH, the boron isotope ratio of coral carbonates, expressed as the δ11B value, needs to be both precise and accurate (2sd〈〈1‰). Boron measurements by Multi‐collector Inductively Coupled Plasma Mass Spectrometry (MC‐ICPMS) requires the boron to be carefully purified before analysis, which is time consuming, and requires specialist training. Here, we use the prepFAST‐MC that enables the automatic extraction of B (up to 25 ng load) from a CaCO3 matrix. Methods Samples were purified using the prepFAST‐MC automated system with a ~25‐μL column of Amberlite IRA743 resin. Boron isotope measurements were performed byMC‐ICPMS. The effect of matrix load, speed of sample load onto the column, and blank contamination were tested to evaluate the effect on the purification process. The optimised protocol was tested on various standards and samples of aragonite corals. Results The blank contribution for the approach is ~60 pg and is negligible given our sample size (〈0.2% sample size). Efficiency of matrix removal is demonstrated with the addition of up to 1.6 mg of dissolved low‐B calcium carbonate to NIST SRM 951 with no impact on the accuracy of δ11B values. The Japanese Geological Survey Porites reference material JCp‐1, boric acid standard NIST SRM 951, and seawater, all processed on the prepFAST‐MC, give δ11B values within error of literature values (δ11BJCp‐1 = 24.31 ±0.20‰ (2sd, n=20); δ11BNIST 951 = ‐0.02 ±0.15‰ (2sd, n=13) and δ11BSeawater = 39.50 ±0.06 ‰ (2sd, n=2)). Results obtained from the coral Siderastrea siderea purified with the prepFAST‐MC show an average offset from the manual ion exchange protocols of Δδ11B=0.01 ±0.28‰ (2sd, n=12). Conclusions Our study demonstrates the capacity of the prepFAST‐MC to generate accurate and reproducible δ11B values for a range of material, without fractionation, with efficient matrix removal and with negligible blank contribution.

Description: A site at the gas hydrate stability limit was investigated offshore northwestern Svalbard to study methane transport in sediment. The site was characterized by chemosynthetic communities (sulfur bacteria mats, tubeworms) and gas venting. Sediments were sampled with in‐situ porewater collectors and by gravity coring followed by analyses of porewater constituents, sediment and carbonate geochemistry, and microbial activity, taxonomy, and lipid biomarkers. Sulfide and alkalinity concentrations showed concentration maxima in near‐surface sediments at the bacterial mat and deeper maxima at the gas vent site. Sediments at the periphery of the chemosynthetic field were characterized by two sulfate‐methane transition zones (SMTZ) at ~204 and 45 cm depth, where activity maxima of microbial anaerobic oxidation of methane (AOM) with sulfate were found. Amplicon sequencing and lipid biomarker indicate that AOM at the SMTZs was mediated by ANME‐1 archaea. A 1D numerical transport reaction model suggests that the deeper SMTZ‐1 formed on centennial scale by vertical advection of methane, while the shallower SMTZ‐2 could only be reproduced by non‐vertical methane injections starting on decadal scale. Model results were supported by age distribution of authigenic carbonates, showing youngest carbonates within SMTZ‐2. We propose that non‐vertical methane injection was induced by increasing blockage of vertical transport or formation of sediment fractures. Our study further suggests that the methanotrophic response to the non‐vertical methane injection was commensurate with new methane supply. This finding provides new information about for the response time and efficiency of the benthic methane filter in environments with fluctuating methane transport.

Description: Heme b is an iron-containing cofactor in hemoproteins that participates in the fundamental processes of photosynthesis and respiration in phytoplankton. Heme b concentrations typically decline in waters with low iron concentrations but due to lack of field data, the distribution of heme b in particulate material in the ocean is poorly constrained. Here we report particulate heme b distributions across the Atlantic Ocean (59.9°N to 34.6°S). Heme b concentrations in surface waters ranged from 0.10 to 33.7 pmol L−1 (median = 1.47 pmol L−1, n = 974) and were highest in regions with a high biomass. The ratio of heme b to particulate organic carbon (POC) exhibited a mean value of 0.44 μmol heme b mol−1 POC. We identified the ratio of 0.10 µmol heme b mol−1 POC as the cut-off between heme b replete and heme b deficient (anemic) phytoplankton. By this definition, we observed anemic phytoplankton populations in the Subtropical South Atlantic and Irminger Basin. Comparison of observed and modelled heme b suggested that heme b could account for between 0.17–9.1% of biogenic iron. Our large scale observations of heme b relative to organic matter provide further evidence of the impact of changes in iron supply on phytoplankton iron status across the Atlantic Ocean.

Description: Ocean bottom pressure (OBP) variability serves as a proxy of ocean mass variability, the knowledge of which is needed in geophysical applications. The question of how well it can be modeled by the present general ocean circulation models on time scales in excess of 1 day is addressed here by comparing the simulated OBP variability with the observed one. To this end, a new multiyear data set is used, obtained with an array of bottom pressure gauges deployed deeply along a transect across the Southern Ocean. We present a brief description of OBP data and show large‐scale correlations over several thousand kilometers at all time scales using daily and monthly averaged data. Annual and semiannual cycles are weak. Close to the Agulhas Retroflection, signals of up to 30 cm equivalent water height are detected. Further south, signals are mostly intermittent and noisy. It is shown that the models simulate consistent patterns of bottom pressure variability on monthly and longer scales except for areas with high mesoscale eddy activity, where high resolution is needed to capture the variability due to eddies. Furthermore, despite good agreement in the amplitude of variability, the in situ and simulated OBP show only modest correlation.

Description: Background: The giant squid (Architeuthis dux; Steenstrup, 1857) is an enigmatic giant mollusc with a circumglobal distribution in the deep ocean, except in the high Arctic and Antarctic waters. The elusiveness of the species makes it difficult to study. Thus, having a genome assembled for this deep-sea-dwelling species will allow several pending evolutionary questions to be unlocked. Findings: We present a draft genome assembly that includes 200 Gb of Illumina reads, 4 Gb of Moleculo synthetic long reads, and 108 Gb of Chicago libraries, with a final size matching the estimated genome size of 2.7 Gb, and a scaffold N50 of 4.8 Mb. We also present an alternative assembly including 27 Gb raw reads generated using the Pacific Biosciences platform. In addition, we sequenced the proteome of the same individual and RNA from 3 different tissue types from 3 other species of squid (Onychoteuthis banksii, Dosidicus gigas, and Sthenoteuthis oualaniensis) to assist genome annotation. We annotated 33,406 protein-coding genes supported by evidence, and the genome completeness estimated by BUSCO reached 92%. Repetitive regions cover 49.17% of the genome. Conclusions: This annotated draft genome of A. dux provides a critical resource to investigate the unique traits of this species, including its gigantism and key adaptations to deep-sea environments.

Description: Parasites are arguably among the strongest drivers of natural selection, constraining hosts to evolve resistance and tolerance mechanisms. Although, the genetic basis of adaptation to parasite infection has been widely studied, little is known about how epigenetic changes contribute to parasite resistance and eventually, adaptation. Here, we investigated the role of host DNA methylation modifications to respond to parasite infections. In a controlled infection experiment, we used the three-spined stickleback fish, a model species for host-parasite studies, and their nematode parasite Camallanus lacustris. We showed that the levels of DNA methylation are higher in infected fish. Results furthermore suggest correlations between DNA methylation and shifts in key fitness and immune traits between infected and control fish, including respiratory burst and functional trans-generational traits such as the concentration of motile sperm. We revealed that genes associated with metabolic, developmental and regulatory processes (cell death and apoptosis) were differentially methylated between infected and control fish. Interestingly, genes such as the neuropeptide FF receptor 2 and the integrin alpha 1 as well as molecular pathways including the Th1 and Th2 cell differentiation were hypermethylated in infected fish, suggesting parasite-mediated repression mechanisms of immune responses. Altogether, we demonstrate that parasite infection contributes to genome-wide DNA methylation modifications. Our study brings novel insights into the evolution of vertebrate immunity and suggests that epigenetic mechanisms are complementary to genetic responses against parasite-mediated selection.

Description: Slower warming in the deep ocean encourages a perception that its biodiversity is less exposed to climate change than that of surface waters. We challenge this notion by analysing climate velocity, which provides expectations for species’ range shifts. We find that contemporary (1955–2005) climate velocities are faster in the deep ocean than at the surface. Moreover, projected climate velocities in the future (2050–2100) are faster for all depth layers, except at the surface, under the most aggressive GHG mitigation pathway considered (representative concentration pathway, RCP 2.6). This suggests that while mitigation could limit climate change threats for surface biodiversity, deep-ocean biodiversity faces an unavoidable escalation in climate velocities, most prominently in the mesopelagic (200–1,000 m). To optimize opportunities for climate adaptation among deep-ocean communities, future open-ocean protected areas must be designed to retain species moving at different speeds at different depths under climate change while managing non-climate threats, such as fishing and mining.

Description: Exploitation and degradation of the mysterious layer between the sunlit ocean surface and the abyss jeopardize fish stocks and the climate.