ALBERT

Description: Accurate and reliable biodiversity estimates of marine zooplankton are a prerequisite to understand how changes in diversity can affect whole ecosystems. Species identification in the deep sea is significantly impeded by high numbers of new species and decreasing numbers of taxonomic experts, hampering any assessment of biodiversity. We used in parallel morphological, genetic, and proteomic characteristics of specimens of calanoid copepods from the abyssal South Atlantic to test if proteomic fingerprinting can accelerate estimating biodiversity. We cross-validated the respective molecular discrimination methods with morphological identifications to establish COI and proteomic reference libraries, as they are a pre-requisite to assign taxonomic information to the identified molecular species clusters. Due to the high number of new species only 37% of the individuals could be assigned to species or genus level morphologically. COI sequencing was successful for 70% of the specimens analysed, while proteomic fingerprinting was successful for all specimens examined. Predicted species richness based on morphological and molecular methods was 42 morphospecies, 56 molecular operational taxonomic units (MOTUs) and 79 proteomic operational taxonomic units (POTUs), respectively. Species diversity was predicted based on proteomic profiles using hierarchical cluster analysis followed by application of the variance ratio criterion for identification of species clusters. It was comparable to species diversity calculated based on COI sequence distances. Less than 7% of specimens were misidentified by proteomic profiles when compared with COI derived MOTUs, indicating that unsupervised machine learning using solely proteomic data could be used for quickly assessing species diversity.

Description: Plankton is a massive and phylogenetically diverse group of thousands of prokaryotes, protists (unicellular eukaryotic organisms), and metazoans (multicellular eukaryotic organisms; Fig. 1). Plankton functional diversity is at the core of various ecological processes, including productivity, carbon cycling and sequestration, nutrient cycling (Falkowski 2012), interspecies interactions, and food web dynamics and structure (D'Alelio et al. 2016). Through these functions, plankton play a critical role in the health of the coastal and open ocean and provide essential ecosystem services. Yet, at present, our understanding of plankton dynamics is insufficient to project how climate change and other human-driven impacts affect the functional diversity of plankton. That limits our ability to predict how critical ecosystem services will change in the future and develop strategies to adapt to these changes.

Description: This study investigates the mesoscale dynamics involved in the 8–11 October 2008 unseasonably strong African dust episode, during which dust was transported to the Iberian Peninsula (IP). We employ observational datasets and a high-resolution Weather Research and Forecasting model coupled with Chemistry simulations. The analysis shows that during 0900–1200 UTC 9 October, a mesoscale convective system developed over the Atlas Mountains and resulted in a southwestward propagating convective cold pool outflow on the southern foothills of the Anti-Atlas, which lifted dust from the source region. Between 1200 and 1800 UTC 9 October, new moist convection was enhanced over the Atlas Mountains due to intensifying confluence among a heat low, moist southwesterly Atlantic sea-breeze front, and northeasterly flow associated with the convective cold pool near western Algeria. This new moist convection intensified the strength of the convective cold pool outflow and haboob, both of which continued propagating southwestward. At 1200 UTC 10 October, the low-pressure system migrated poleward on the southern slopes of the Anti-Atlas Mountains in association with a mountain-plains solenoidal circulation due to the daytime differential heating between the southern slopes of the Anti-Atlas and nearby atmosphere. The deepening low-pressure and strengthening Atlantic sea-breeze redirected an equatorward advancing dust plume into the poleward direction. The dust plume ultimately crossed the Saharan Atlas Mountains on 11 October and finally impacted the IP. Key Points: - WRF-Chem simulation of an unseasonably strong haboob on the southern slopes of the Atlas Mountains - The equatorward-advancing dust plume was recirculated in the poleward direction by an Atlantic sea-breeze front - The Atlantic sea-breeze front and an intensified upper-level cutoff vortex are instrumental for dust transport over the Iberian Peninsula

Description: An extensive, reprocessed two‐dimensional (2D) seismic data set was utilized together with available well data to study the Tiddlybanken Basin in the southeastern Norwegian Barents Sea, which is revealed to be an excellent example of base salt rift structures, evaporite accumulations and evolution of salt structures. Late Devonian–early Carboniferous NE‐SW regional extensional stress affected the study area and gave rise to three half‐grabens that are separated by a NW‐SE to NNW‐SSE trending horst and an affiliated interference transfer zone. The arcuate nature of the horst is believed to be the effect of pre‐existing Timanian basement grain, whereas the interference zone formed due to the combined effect of a Timanian (basement) lineament and the geometrical arrangement of the opposing master faults. The interference transfer zone acted as a physical barrier, controlling the facies distribution and sedimentary thickness of three‐layered evaporitic sequences (LES). During the late Triassic, the northwestern part of a salt wall was developed due to passive diapirism and its evolution was influenced by halite lithology between the three‐LES. The central and southeastern parts of the salt wall did not progress beyond the pedestal stage due to lack of halite in the deepest evaporitic sequence. During the Triassic–Jurassic transition, far‐field stresses from the Novaya Zemlya fold‐and‐thrust belt reactivated the pre‐salt Carboniferous rift structures. The reactivation led to the development of the Signalhorn Dome, rejuvenated the northwestern part of the salt wall and affected the sedimentation rates in the southeastern broad basin. The salt wall together with the Signalhorn Dome and the Carboniferous pre‐salt structures were again reactivated during post‐Early Cretaceous, in response to regional compressional stresses. During this main tectonic inversion phase, the northwestern and southeastern parts of the salt wall were rejuvenated; however, salt reactivation was minimized towards the interference transfer zone beneath the centre of the salt wall.

Description: Despite the implication of aerosols for the radiation budget, there are persistent differences in data for the aerosol optical depth (τ) for 1998–2019. This study presents a comprehensive evaluation of the large-scale spatio-temporal patterns of mid-visible τ from modern data sets. In total, we assessed 94 different global data sets from eight satellite retrievals, four aerosol-climate model ensembles, one operational ensemble product, two reanalyses, one climatology and one merged satellite product. We include the new satellite data SLSTR and aerosol-climate simulations from the Coupled Model Intercomparison Project Phase 6 (CMIP6) and the Aerosol Comparisons between Observations and Models Phase 3 (AeroCom-III). Our intercomparison highlights model differences and observational uncertainty. Spatial mean τ for 60°N – 60°S ranges from 0.124 to 0.164 for individual satellites, with a mean of 0.14. Averaged τ from aerosol-climate model ensembles fall within this satellite range, but individual models do not. Our assessment suggests no systematic improvement compared to CMIP5 and AeroCom-I. Although some regional biases have been reduced, τ from both CMIP6 and AeroCom-III are for instance substantially larger along extra-tropical storm tracks compared to the satellite products. The considerable uncertainty in observed τ implies that a model evaluation based on a single satellite product might draw biased conclusions. This underlines the need for continued efforts to improve both model and satellite estimates of τ, for example, through measurement campaigns in areas of particularly uncertain satellite estimates identified in this study, to facilitate a better understanding of aerosol effects in the Earth system. Key Points: - Present-day patterns in aerosol optical depth differ substantially between 94 modern global data sets - The range in spatial means from individual satellites is −11% to +17% of the multi-satellite mean - Spatial means from climate model intercomparison projects fall within the satellite range but strong regional differences are identified

Description: Solar radiation received at the Earth's surface (Rs) is comprised of two components, the direct radiation (Rd) and the diffuse radiation (Rf). Rd, the direct beam from the sun, is essential for concentrated solar power generation. Rf, scattered by atmospheric molecules, aerosols, or cloud droplets, has a fertilization effect on plant photosynthesis. But how Rd and Rf change diurnally is largely unknown owing to the lack of long-term measurements. Taking advantage of 22 years of homogeneous hourly surface observations over China, this study documents the climatological means and evolutions in the diurnal cycles of Rd and Rf since 1993, with an emphasis on their implications for solar power and agricultural production. Over the solar energy resource region, we observe a loss of Rd which is relatively large near sunrise and sunset at low solar elevation angles when the sunrays pass through the atmosphere on a longer pathway. However, the concentrated Rd energy covering an average 10-hr period around noon during a day is relatively unaffected. Over the agricultural crop resource region, the large amounts of clouds and aerosols scattering more of the incoming light result in Rf taking the main proportion of Rs during the whole day. Rf resources and their fertilization effect in the main crop region of China further enhances since 1993 over almost all hours of the day. Key Points: - The loss of direct radiation over China since 1993 is relatively large at sunrise and sunset with little effect on solar power generation - The diffuse component dominates solar radiation normally near sunrise and sunset, but for the whole day over the main sown area of China - The diffuse fraction is further enhanced in the main sown area of China over almost all hours of the day since 1993

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: 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: 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: 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: 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: 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: 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: 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: 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: 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: 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: 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: 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: 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: 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: 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: 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: 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: 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: Most arcs show systematic temporal and spatial variations in magmatism with clear shifts in igneous rock compositions between those of the magmatic front (MF) and those in the backarc (BA). It is unclear if similar magmatic polarity is seen for extensional continental arcs. Herein, we use geochemical and isotopic characteristics coupled with zircon U‐Pb geochronology to identify the different magmatic style of the Iran convergent margin, an extensional system that evolved over 100 Myr. Our new and compiled U‐Pb ages indicate that major magmatic episodes for the NE Iran BA occurred at 110–80, 75–50, 50–35, 35–20, and 15–10 Ma. In contrast to NE Iran BA magmatic episodes, compiled data from MF display two main magmatic episodes at 95–75 and 55–5 Ma, indicating more continuous magmatism for the MF than for the BA. We show that Paleogene Iran serves as a useful example of a continental arc under extension. Our data also suggest that there is not a clear relationship between the subduction velocity of Neotethyan Ocean beneath Iran and magmatic activity in Iran. Our results imply that the isotopic compositions of Iran BA igneous rocks do not directly correspond to the changes in tectonic processes or geodynamics, but other parameters such as the composition of lithosphere and melt source(s) should be considered. In addition, changes in subduction zone dynamics and contractional versus extensional tectonic regimes influenced the composition of MF and BA magmatic rocks. These controls diminished the geochemical and isotopic variations between the magmatic front and backarc.

Description: The use of animal-attached data loggers to quantify animal movement has increased in popularity and application in recent years. High-resolution tri-axial acceleration and magnetometry measurements have been fundamental in elucidating fine-scale animal movements, providing information on posture, traveling speed, energy expenditure, and associated behavioral patterns. Heading is a key variable obtained from the tandem use of magnetometers and accelerometers, although few field investigations have explored fine-scale changes in heading to elucidate differences in animal activity (beyond the notable exceptions of dead-reckoning). This paper provides an overview of the value and use of animal heading and a prime derivative, angular velocity about the yaw axis, as an important element for assessing activity extent with potential to allude to behaviors, using “free-ranging” Loggerhead turtles (Caretta caretta) as a model species. We also demonstrate the value of yaw rotation for assessing activity extent, which varies over the time scales considered and show that various scales of body rotation, particularly rate of change of yaw, can help resolve differences between fine-scale behavior-specific movements. For example, oscillating yaw movements about a central point of the body's arc implies bouts of foraging, while unusual circling behavior, indicative of conspecific interactions, could be identified from complete revolutions of the longitudinal axis. We believe this approach should help identification of behaviors and “space-state” approaches to enhance our interpretation of behavior-based movements, particularly in scenarios where acceleration metrics have limited value, such as for slow-moving animals.

Description: Key Points:  AUV mapping along the north rift identifies hummocky mounds with pillows, channelized, and inflated flows formed during the 2015 eruption.  Impulsive sounds formed by lava/seawater interaction track mound growth over a 28 day period with an average extrusion rate of 22-45 m3s-1.  The sounds record the history of flow advancement and inflation from multiple eruptive centers, and is used to infer volcanic activity style. Quantifying eruption dynamics in submarine environments is challenging. During the 2015 eruption of Axial Seamount, the formation of hummocky mounds along the north rift was accompanied by tens‐of‐thousands of impulsive acoustic signals generated by the interaction of lava and seawater. A catalog of these sounds was integrated with detailed seafloor mapping to better understand eruptive processes in time and space. Mounds grew over a period of 28 days with average extrusion rates of 22 to 45 m3s‐1. The most distant mounds, ~ 9.5 to 15.5 km down rift from the caldera, grew primarily over the first few days of the eruption. The focus of eruptive activity then retreated ~5 km toward the caldera where it was sustained. Mounds are constructed as a series of superimposed lobes formed through alternating periods of flow inflation, generating up to 30‐m‐thick hummocks, and periods of flow advancement, with 〈0.02 ms‐1 average speeds typically observed.

Description: Oceanic transform faults and fracture zones represent major bathymetric features that keep the records of past and present strike‐slip motion along conservative plate boundaries. Although they play an important role in ridge segmentation and evolution of the lithosphere, their structural characteristics, and their variation in space and time, are poorly understood. To address some of the unknowns, we conducted interdisciplinary geophysical studies in the equatorial Atlantic Ocean, the region where some of the most prominent transform discontinuities have been developing. Here we present the results of the data analysis in the vicinity of the Chain Fracture Zone (FZ), on the South American Plate. The crustal structure across the Chain FZ, at the contact between ~10 and 24 Ma oceanic lithosphere, is sampled along seismic reflection and refraction profiles. We observe that the crustal thickness within and across the Chain FZ ranges from ~4.6‐5.9 km, which compares with the observations reported for slow‐slipping transform discontinuities globally. We attribute this presence of close to normal oceanic crustal thickness within fracture zones to the mechanism of lateral dike propagation, previously considered to be valid only in fast‐slipping environments. Furthermore, the combination of our results with other datasets enabled us to extend the observations to morpho‐tectonic characteristics on a regional scale. Our broader view suggests that the formation of the transverse ridge is closely associated with a global plate reorientation that was also responsible for the propagation and for shaping lower‐order Mid‐Atlantic Ridge segmentation around the equator.

Description: Aim: Global biodiversity and ecosystems are highly impacted by anthropogenic activities, such as climate change and introduction of non-indigenous species. As numerous species from the Ponto-Caspian region have established in the North and Baltic Seas, as well as in the Laurentian Great Lakes, there have been large number of studies examining environmental tolerance of these species to determine their future potential to spread. However, many of those studies were conducted only on adult stages, while neglecting the possibility that early life history stages might not be equally resilient. Location: Northern European, Ponto-Caspian and North American regions. Methods: To determine whether juveniles would demonstrate the same environmental tolerance as their parents, we examined the salinity tolerance of adults and juveniles of one Northern European (Gammarus salinus), one Ponto-Caspian (Pontogammarus maeoticus) and one North American species (Gammarus tigrinus). Additionally, we compared our study to that of Paiva et al. (Global Change Biology, 24, 2018, 2708), who tested the salinity tolerance of the same species using only adults. Results: Our study determined that both adults and juveniles of all three species tolerated wide ranges of salinity, with juveniles of G. salinus tolerating only slightly narrower salinity range than their parents, while those of P. maeoticus and G. tigrinus much narrower range. Additionally, we determined better survival and higher growth rates of juveniles of G. salinus in higher salinities and better survival of P. maeoticus in lower salinities. Main conclusions: Based on juvenile salinity tolerance, our study further supported findings of Paiva et al. (2018), where Northern European species may be adapted to marine, while Ponto-Caspian to lower saline and freshwater environments. The North American species is probably adapted to intermediate salinities. As juveniles do not tolerate the same salinity stress as adults, we emphasize the importance of testing all life history stages when predicting species resilience to environmental stressors.

Description: Neogloboquadrina pachyderma is the dominant species of planktonic foraminifera found in polar waters and is therefore invaluable for paleoceanographic studies of the high latitudes. However, the geochemistry of this species is complicated due to the development of a thick calcite crust in its final growth stage and at greater depths within the water column. We analyzed the in situ Mg/Ca and δ18O in discrete calcite zones using LA‐ICP‐MS, EPMA and SIMS within modern N. pachyderma shells from the highly dynamic Fram Strait and the seasonally isothermal/isohaline Irminger Sea. Here we compare shell geochemistry to the measured temperature, salinity and δ18Osw in which the shells calcified to better understand the controls on N. pachyderma geochemical heterogeneity. We present a relationship between Mg/Ca and temperature in N. pachyderma lamellar calcite that is significantly different than published equations for shells that contained both crust and lamellar calcite. We also document highly variable SIMS δ18O results (up to a 3.3‰ range in single shells) on plankton tow samples which we hypothesize is due to the granular texture of shell walls. Finally, we document that the δ18O of the crust and lamellar calcite of N. pachyderma from an isothermal/isohaline environment are indistinguishable from each other, indicating that shifts in N. pachyderma δ18O are primarily controlled by changes in environmental temperature and/or salinity rather than differences in the sensitivities of the two calcite types to environmental conditions.

Description: This study proposes a new process‐based framework to characterize and classify runoff events of various magnitudes occurring in a wide range of catchments. The framework uses dimensionless indicators that characterize space–time dynamics of precipitation events and their spatial interaction with antecedent catchment states, described as snow cover, distribution of frozen soils, and soil moisture content. A rigorous uncertainty analysis showed that the developed indicators are robust and regionally consistent. Relying on covariance‐ and ratio‐based indicators leads to reduced classification uncertainty compared to commonly used (event‐based) indicators based on absolute values of metrics such as duration, volume, and intensity of precipitation events. The event typology derived from the proposed framework is able to stratify events that exhibit distinct hydrograph dynamics even if streamflow is not directly used for classification. The derived typology is therefore able to capture first‐order controls of event runoff response in a wide variety of catchments. Application of this typology to about 180,000 runoff events observed in 392 German catchments revealed six distinct regions with homogeneous event type frequency that match well regions with similar behavior in terms of runoff response identified in Germany. The detected seasonal pattern of event type occurrence is regionally consistent and agrees well with the seasonality of hydroclimatic conditions. The proposed framework can be a useful tool for comparative analyses of regional differences and similarities of runoff generation processes at catchment scale and their possible spatial and temporal evolution.

Description: We investigate how the atmosphere is affected by the cold sea surface temperature (SST) anomalies of the Benguela Niñas using reanalysis data and a high-resolution atmospheric model. A composite analysis of reanalysis data based on five Benguela Niña events (5 years out of 39 years for 1979–2017) reveals that the rainfall along the Angolan coast is reduced significantly and other anomalies of precipitation are detected over the African continent and equatorial Atlantic. Those anomalies can be explained by the anomalies of vertically-integrated moisture flux (VIMF) and its divergence and convergence. Additionally, the Namibian low-level cloud and sea level pressure (SLP) are enhanced around the ABFZ by the Benguela Niñas. A simulation of the model forced by the cold SST anomalies of the Benguela Niñas reproduces the reduction (enhancement) of rainfall and VIMF divergent (convergent) anomaly over the Angolan coastal region (Gulf of Guinea). However, the other anomalies over the African continent are not significant. It is suggested that the effects of the Benguela Niñas are limited along southwestern coast of Africa. Composite analysis of reanalyses of rainfall anomalies associated with Benguela Niña events shows pattern over the other regions which might be induced by larger-scale atmospheric anomaly. This large-scale atmospheric anomaly can be linked with the South Atlantic Anticyclone and/or forced by the teleconnection from the tropical Pacific, for instance. The increment of the Namibian low-level cloud and SLP is simulated by the numerical experiments consistently. The low-level cloud becomes more frequent between 950 and 900 hPa and the radiative cooling by longwave radiation is reinforced during the Benguela Niñas events. In contrast, the cloud formation around 850 hPa is reduced and radiative cooling is weakened. It is indicated that this change in the cooling rate possibly induces the strengthening in the inversion layer over the cold SST anomalies.

Description: Nitrous oxide (N2O) is a strong greenhouse gas and an ozone depleting agent. In marine environments, N2O is produced biologically via ammonium oxidation, nitrite, and nitrate reduction. The relative importance of these principle production pathways is strongly influenced by oxygen availability. We conducted 15N tracer experiments of N2O production in parallel with measurements of N2O concentration and natural abundance isotopes/isotopomers in Saanich Inlet, a seasonally anoxic fjord, to investigate how temporal and vertical oxygen gradients regulate N2O production pathways and rates. In April, June, and August 2018, the depth of the oxic‐anoxic interface (dissolved oxygen = 2.5 μmol L−1 isoline) progressively deepened from 110 to 160 m. Within the oxygenated and suboxic water column, N2O supersaturation coincided with peak ammonium oxidation activity. Conditions in the anoxic deep water were potentially favorable to N2O production from nitrate and nitrite reduction, but N2O undersaturation was observed indicating that N2O consumption exceeded rates of production. In October, tidal mixing introduced oxygenated water from outside the inlet, displacing the suboxic and anoxic deep water. This oxygenation event stimulated N2O production from ammonium oxidation and increased water column N2O supersaturation while inhibiting nitrate and nitrite reduction to N2O. Results from 15N tracer incubation experiments and natural abundance isotopomer measurements both implicated ammonium oxidation as the dominant N2O production pathway in Saanich Inlet, fueled by high ammonium fluxes (0.6–3.5 nmol m−2 s−1) from the anoxic depths. Partial denitrification contributed little to water column N2O production because of low availability of nitrate and nitrite.

Description: In a widely‐held conception, the biological carbon pump (BCP) is equal to the export of organic matter out of the euphotic zone. Using global ocean‐atmosphere model experiments we show that the change in export production is a poor measure of the biological pump's feedback to the atmosphere. The change in global true oxygen utilization (TOU), an integrative measure of the imprint of the BCP on marine oxygen, however, is in good agreement with the net change in the biogenic air‐sea flux of oxygen. Since TOU correlates very well with apparent oxygen utilization (AOU) in our experiments, we propose to measure the change of AOU from data of global float programs to monitor the feedback of the BCP to the atmosphere. For the current ocean we estimate that BCP changes effect a CO2 uptake by the ocean in the range of 0.07 to 0.14 GtC/yr.

Description: Rapid increases in upper 700‐m Indian Ocean heat content (IOHC) since the 2000s have focused attention on its role during the recent global surface warming hiatus. Here, we use ocean model simulations to assess distinct multidecadal IOHC variations since the 1960s and explore the relative contributions from wind stress and buoyancy forcing regionally and with depth. Multidecadal wind forcing counteracted IOHC increases due to buoyancy forcing from the 1960s to the 1990s. Wind and buoyancy forcing contribute positively since the mid‐2000s, accounting for the drastic IOHC change. Distinct timing and structure of upper ocean temperature changes in the eastern and western Indian Ocean are linked to the pathway how multidecadal wind forcing associated with the Interdecadal Pacific Oscillation is transmitted and affects IOHC through local and remote winds. Progressive shoaling of the equatorial thermocline—of importance for low‐frequency variations in Indian Ocean Dipole occurrence—appears to be dominated by multidecadal variations in wind forcing.

Description: The Svalbard margin represents one of the northernmost gas hydrate provinces worldwide. Vestnesa Ridge (VR) and Svyatogor Ridge (SR) west of Svalbard are two prominent sediment drifts showing abundant pockmarks and sites of seismic chimney structures. Some of these sites at VR are associated with active gas venting and were the focus of drilling and coring with the seafloor‐deployed MARUM‐MeBo70 rig. Understanding the nature of fluid migration and gas hydrate distribution requires (amongst other parameters) knowledge of the thermal regime and in situ gas and pore‐fluid composition. In situ temperature data were obtained downhole at a reference site at VR defining a geothermal gradient of ~78 mK m‐1 (heat flow ~95 mW m‐2). Additional heat‐probe data were obtained to describe the thermal regime of the pockmarks. The highest heat flow values were systematically seen within pockmark depressions and were uncorrelated to gas venting occurrences. Heat flow within pockmarks is typically ~20 mW m‐2 higher than outside pockmarks. Using the downhole temperature data and gas compositions from drilling we model the regional base of the gas hydrate stability zone (BGHSZ). Thermal modeling including topographic effects suggest a BGHSZ up to 40 m deeper than estimated from seismic data. Uncertainties in sediment properties (velocity and thermal conductivity) are only partially explaining the mismatch. Capillary effects due to small sediment grain sizes may shift the free gas occurrence above the equilibrium BGHSZ. Changes in gas composition or pore fluid salinity at greater depth may also explain the discrepancy in observed and modeled BGHSZ.

Description: The neodymium (Nd) isotopic composition of seawater is a valuable tool for the reconstruction of past water mass provenance and hence deep water geometry. A meaningful interpretation of Nd isotope down‐core records requires knowledge of potential variations of water mass end member characteristics. While often assumed temporally constant, recent investigations revealed glacial‐interglacial variability of the northern and southern Nd isotope end members in the Atlantic. These new constraints have a strong influence on the interpretation of the Atlantic deep water mass evolution, yet the processes responsible for the end member shifts remain uncertain. Here we combine a new compilation of Atlantic Nd isotope reconstructions of the early Holocene with the Nd‐enabled Bern3D model to quantify the recently proposed hypothesis of a northern Nd isotope end member shift during the early Holocene. We achieve the best model‐data fit with a strong increase of the Nd flux in the northern high latitudes by a factor of 3 to 4, which lowers the northern end member signature by about 1 ε‐unit. Our findings thus agree with the rationale that glacially weathered material entered the northern Northwest Atlantic after the ice sheets retreated late in the deglaciation and released substantial amounts of unradiogenic Nd as suggested previously. Further, we find that variations in the strength of the Atlantic Meridional Overturning Circulation (AMOC) cannot reproduce the observed Nd isotope excursions of the compiled data, ruling out an early Holocene AMOC “overshoot.”

Description: Large amounts of methane, a potent greenhouse gas, are stored in hydrates beneath the seafloor. Sea level changes can trigger massive methane release into the ocean. It is not clear, however, whether surficial seafloor processes can cause comparable discharge. Previously, fluid migration was difficult to study due to a lack of spatially dense seismic and thermal observations. Here we examine a gas hydrate site at Four‐Way‐Closure Ridge off SW Taiwan using a high‐resolution 3‐D seismic cube, together with bottom‐simulating reflections (BSRs) mapped in the cube, a thermal probe data set, and 3‐D thermal modeling results. We document, on a scale of tens of meters, the interaction between surficial sedimentary processes, fluid flow, and a dynamic gas hydrate system. Fluid migrates upward through dipping permeable strata in the limb, the slope basin, and along thrust faults and ridge‐top normal faults. The seismic data also reveal several double BSRs that underlie seabed sedimentary sliding and depositional features. Abrupt changes in subsurface pressure and temperature due to the rapid seabed sedimentary processes can cause a rapid shift of the base of the gas hydrate stability zone. This shift may be either downward or upward and would result in the accumulation or dissociation of hydrate in sediments sandwiched by the double BSRs, respectively. We propose that dynamic surficial processes on the seafloor together with shallow focused fluid flow affect hydrate distribution and saturation at depth and may even result in methane expulsion into the ocean if such localized features are common along convergent plate boundaries.

Description: Key points: • Detectable self-potential and deformation anomalies result from poroelastic responses of volcanicaquifers by subsurface pressurization • Self-potential amplitudes and their polarity are sensitive to magmatic stressing and pressure sourceorientation (dike vs. sill) • Our multi-physics approach provides new insights into pre-eruptive processe Pre‐eruptive electrical signals at active volcanoes are generally interpreted in terms of electrokinetic processes. Spatio‐temporal self‐potential (SP) signals can be caused by strain‐induced fluid flow in volcanic aquifers, however, previous studies lack the quantitative assessments of these phenomena and the underpinning poroelastic responses. Here we use Finite‐Element Analysis to study poroelastic responses induced by subsurface stressing from sill and dike sources by jointly solving for ground displacements, pore pressure and SP signals. We evaluate the influence of pressure source orientation on the poroelastic response in two different volcanic aquifers (pyroclastic and lava flow) to provide insights on emergent geodetic and SP signals and their sensitivity to governing parameters. Strain‐induced SP amplitudes deduced from a reference parameter set vary in both aquifer models and are of negative polarity (‐0.35 mV and ‐22.6 mV) for a pressurized dike and of positive polarity (+4 mV and +20 mV) for a pressurized sill. Importantly, we find uniquely different SP and ground displacement patterns from either sill or dike intrusions. Our study shows that SP signals are highly sensitive to the subsurface Young's modulus, streaming potential coupling coefficient and electrical conductivity of the poroelastic domains. For the set of parameters tested, the dike model predicts SP amplitudes of up to ‐947 mV which are broadly representative of recorded amplitudes from active volcanoes. Our study demonstrates that electrokinetic processes reflect magma‐induced stress and strain variations and highlights the potential of joint geodetic and SP studies to gain new insights on causes of volcanic unrest.

Description: The modern history of North Atlantic sea surface temperature shows variability coinciding with changes in air temperature and rainfall over the Northern Hemisphere. There is a debate about this variability and, in particular, whether it is internal to the ocean‐atmosphere system or is forced by external factors (natural and anthropogenic). Here we present a temperature record, obtained using the Sr/Ca ratio measured in a skeleton of a sclerosponge, that shows agreement with the instrumental record over the past 150 years as well as multidecadal temperature variability over the last 600 years. Comparison with climate simulations of the last millennium shows that large cooling events recorded, in the sclerosponge, are consistent with natural (primarily volcanic activity) and anthropogenic forcings. There are, however, multidecadal periods not connected to current estimates of external forcing over the last millennium allowing for alternative explanations, such as internally driven changes in ocean and atmospheric circulation.

Description: Mixotrophic organisms are increasingly recognized as important components of ecosystems, but the factors controlling their nutrition pathways (in particular their autotrophy-heterotrophy balance) are little known. Both autotrophy and heterotrophy are expected to respond to density-dependent mechanisms but not necessarily in the same direction and/or strength. We hypothesize that the autotrophy-heterotrophy balance of mixotrophic organisms might therefore be a function of population densities. To investigate this relationship, we sampled mixotrophic jellyfish holobionts (host, Mastigias papua etpisoni; symbiont, Cladocopium sp.) in a marine lake (Palau, Micronesia) on six occasions (from 2010 to 2018). Over this period, population densities varied similar to 100 fold. We characterized the nutrition of the holobionts using the delta C-13 and delta N-15 signatures as well as C:N ratios. delta C-13 values increased and delta N-15 values decreased with increasing population densities (respectively, R-2 = 0.86 and 0.70, P 〈 0.05). Although less distinct, C:N ratios increased with increasing population densities (R-2 = 0.59, 0.1 〉 P 〉 0.05). This indicates that the autotrophy-heterotrophy balance tends toward autotrophy when population densities increase. We propose that the availability of zooplanktonic prey is the main driver of this pattern. These results demonstrate that the autotrophy-heterotrophy balance of mixotrophic jellyfishes can be tightly regulated by density-dependent mechanisms.

Description: Ice nucleating particles (INPs) affect the radiative properties of cold clouds. Knowledge concerning their concentration above ground level and their potential sources is scarce. Here we present the first highly temperature resolved ice nucleation spectra of airborne samples from an aircraft campaign during late winter in 2018. Most INP spectra featured low concentration levels (〈3 · 10−4 L−1 at −15°C). However, we also found INP concentrations of up to 1.8·10−2 L−1 at −15°C and freezing onsets as high as −7.5°C for samples mainly from the marine boundary layer. Shape and onset temperature of the ice nucleation spectra of those samples as well as heat sensitivity hint at biogenic INP. Colocated measurements additionally indicate a local marine influence rather than long-range transport. Our results suggest that even in late winter above 80°N a local marine source for biogenic INP, which can efficiently nucleate ice at high temperatures, is present.

Description: Investigation of active faults is crucial for the seismic hazard assessment and, in the case of volcanic belts, it provides a deeper understanding of the interactions between volcanism and tectonic faulting. In this study, we report the results of the first paleoseismological and tephrochronological investigation undertaken on Holocene faulting in Kamchatka's volcanic belts. The studied trenches and additional excavations are located along the axial fault zone of the Eastern Volcanic Front, where the earlier dated tephra layers provide a robust age control of the faulting events. Electron microprobe analysis of glass from 22 tephra samples permitted correlations among the disparate tephra profiles for constructing a summary tephra sequence. The latter, together with published geochronological data, allowed the construction of a Bayesian age model. Detailed examination of the tephra layers deformed by faulting allowed us to reconstruct and date six faulting events with the offsets of 1 to 20 cm indicating paleoearthquakes with magnitudes of Mw 〈 5.4. Holocene crustal seismicity of the Eastern Volcanic Front manifests temporal clustering rather than a uniform flux of events. However, no correlation between dated seismic events and the largest Holocene eruptions of proximal volcanoes was observed.

Description: Many marine microbial eukaryotes combine photosynthetic with phagotrophic nutrition, but incomplete understanding of such mixotrophic protists, their functional diversity, and underlying physiological mechanisms limits the assessment and modeling of their roles in present and future ocean ecosystems. We developed an experimental system to study responses of mixotrophic protists to availability of living prey and light, and used it to characterize contrasting physiological strategies in two stramenopiles in the genus Ochromonas. We show that oceanic isolate CCMP1393 is an obligate mixotroph, requiring both light and prey as complementary resources. Interdependence of photosynthesis and heterotrophy in CCMP1393 comprises a significant role of mitochondrial respiration in photosynthetic electron transport. In contrast, coastal isolate CCMP2951 is a facultative mixotroph that can substitute photosynthesis by phagotrophy and hence grow purely heterotrophically in darkness. In contrast to CCMP1393, CCMP2951 also exhibits a marked photoprotection response that integrates non‐photochemical quenching and mitochondrial respiration as electron sink for photosynthetically produced reducing equivalents. Facultative mixotrophs similar to CCMP2951 might be well adapted to variable environments, while obligate mixotrophs similar to CCMP1393 appear capable of resource efficient growth in oligotrophic ocean environments. Thus, the responses of these phylogenetically close protists to the availability of different resources reveals niche differentiation that influences impacts in food webs and leads to opposing carbon cycle roles.

Description: Numerical simulations of ocean biogeochemical cycles need to adequately represent particle sinking velocities (SV). For decades, Stokes' Law estimating particle SV from density and size has been widely used. But while Stokes' Law holds for small, smooth, and rigid spheres settling at low Reynolds number, it fails when applied to marine aggregates complex in shape, structure, and composition. Minerals and zooplankton can alter phytoplankton aggregates in ways that change their SV, potentially improving the applicability of Stokes' models. Using rolling cylinders, we experimentally produced diatom aggregates in the presence and absence of minerals and/or microzooplankton. Minerals and to a lesser extent microzooplankton decreased aggregate size and roughness and increased their sphericity and compactness. Stokes' Law parameterized with a fractal porosity modeled adequately size‐SV relationships for mineral‐loaded aggregates. Phytoplankton‐only aggregates and those exposed to microzooplankton followed the general Navier‐Stokes drag equation suggesting an indiscernible effect of microzooplankton and a drag coefficient too complex to be calculated with a Stokes' assumption. We compared our results with a larger data set of ballasted and nonballasted marine aggregates. This confirmed that the size‐SV relationships for ballasted aggregates can be simulated by Stokes' models with an adequate fractal porosity parameterization. Given the importance of mineral ballasting in the ocean, our findings could ease biogeochemical model parameterization for a significant pool of particles in the ocean and especially in the mesopelagic zone where the particulate organic matter : mineral ratio decreases. Our results also reinforce the importance of accounting for porosity as a decisive predictor of marine aggregate SV.

Description: Ocean surface partial pressure of carbon dioxide (pCO2) is a key factor controlling air–sea CO2 fluxes. Most surface pCO2 data are collected with relatively large and complex air–water equilibrators coupled to stand‐alone infrared analyzers installed on Ships of OPportunity (SOOP‐CO2). This approach has proven itself through years of successful deployments, but expansion and sustainability of the future measurement network faces challenges in terms of certification, autonomy, and maintenance, which motivates development of new systems. Here, we compare performance of three underway pCO2 measurement systems (General Oceanics, SubCtech, and Pro‐Oceanus), including a recently developed compact flow‐through, sensor‐based system. The systems were intercompared over a period of 34 days during two crossings of the subpolar North Atlantic Ocean. With a mean difference from the General Oceanics system of −5.7 ± 4.0 μatm (Pro‐Oceanus) and −4.7 ± 2.9 μatm (SubCtech) during the 1st crossing, our results indicate potential for good agreement between the systems. The study highlighted the challenge of assuring accuracy over long periods of time, particularly seen in a worse agreement during the 2nd crossing, and revealed a number of sources of systematic errors. These can influence accuracy of the measurements, agreement between systems and include slow response of membrane‐based systems to pCO2 changes, “within‐ship” respiration due to biofouling, and bias in measurement of the temperature of equilibration. These error sources can be controlled or corrected for, however, if unidentified, their magnitude can be significant relative to accuracy criteria assigned to the highest‐quality data in global databases. The advantages of the compact flow‐through system are presented along with a discussion of future solutions for improving data quality.

Description: The Atlantic Meridional Overturning Circulation (AMOC) is a key mechanism of heat, freshwater, and carbon redistribution in the climate system. The precept that the AMOC has changed abruptly in the past, notably during and at the end of the last ice age, and that it is “very likely” to weaken in the coming century due to anthropogenic climate change is a key motivation for sustained observations of the AMOC. This paper reviews the methodology and technology used to observe the AMOC and assesses these ideas and systems for accuracy, shortcomings, potential improvements, and sustainability. We review hydrographic techniques and look at how these traditional techniques can meet modern requirements. Transport mooring arrays (TMAs) provide the “gold standard” for sustained AMOC observing, utilizing dynamic height, current meter, and other instrumentation and techniques to produce continuous observations of the AMOC. We consider the principle of these systems and how they can be sustained and improved into the future. Techniques utilizing indirect measurements, such as satellite altimetry, coupled with in situ measurements, such as the Argo float array, are also discussed. Existing technologies that perhaps have not been fully exploited for estimating AMOC are reviewed and considered for this purpose. Technology is constantly evolving, and we look to the future of technology and how it can be deployed for sustained and expanded AMOC measurements. Finally, all of these methodologies and technologies are considered with a view to a sustained and sustainable future for AMOC observation.

Description: The Calabrian subduction zone is one of the narrowest arcs on Earth and a key area to understand the geodynamic evolution of the Mediterranean and other marginal seas. Here in the Ionian Sea, the African plate subducts beneath Eurasia. Imaging the boundary between the downgoing slab and the upper plate along the Calabrian subduction zone is important for assessing the potential of the subduction zone to generate mega‐thrust earthquakes and was the main objective of this study. Here we present and analyze the results from a 380 km long, wide‐angle seismic profile spanning the complete subduction zone, from the deep Ionian Basin and the accretionary wedge to NE Sicily, with additional constraints offered by 3‐D Gravity modeling and the analysis of earthquake hypocenters. The velocity model for the wide‐angle seismic profile images thin oceanic crust throughout the basin. The Calabrian backstop extends underneath the accretionary wedge to about 100 km SE of the coast. The seismic model was extended in depth using earthquake hypocenters. The combined results indicate that the slab dip increases abruptly from 2‐3° to 60‐70° over a distance of ≤50 km underneath the Calabrian backstop. This abrupt steepening is likely related to the roll‐back geodynamic evolution of the narrow Calabrian slab which shows great similarity to the shallow and deep geometry of the Gibraltar slab.

Description: A method using a linear shallow water model is presented for decomposing the temporal variability of the barotropic streamfunction in a high‐resolution ocean model. The method is based in the vertically‐averaged momentum equations and is applied to the time series of annual mean streamfunction from the model configuration VIKING20 for the northern North Atlantic. An important result is the role played by the nonlinear advection terms in VIKING20 for driving transport. The method is illustrated by examining how the Gulf Stream transport in the recirculation region responds to the winter North Atlantic Oscillation (NAO). While no statistically significant response is found in the year overlapping with the winter NAO index, there is a tendency for the Gulf Stream transport to increase as the NAO becomes more positive. This becomes significant in lead years 1 and 2 when the mean flow advection (MFA) and eddy momentum flux (EMF) contributions, associated with nonlinear momentum advection, dominate. Only after 2 years, does the potential energy (PE) term, associated with the density field, start to play a role and it is only after 5 years that the transport dependence on the NAO ceases to be significant. It is also shown that the PE contribution to the transport streamfunction has significant memory of up to 5 years in the Labrador and Irminger Seas. However, it is only around the northern rim of these seas that VIKING20 and the transport reconstruction exhibit similar memory. This is due to masking by the MFA and EMF contributions.

Description: A new autotrophic hydrogen‐oxidizing Chromatiaceae bacterium, namely bacterium CTD079, was enriched from a water column sample at 1500 m water depth in the southern Pacific Ocean. Based on the phylogeny of 16S rRNA genes it was closely related to a scaly snail endosymbiont (99.2% DNA sequence identity) whose host so far is only known to colonize hydrothermal vents along the Indian ridge. The average nucleotide identity between the genomes of CTD079 and the snail endosymbiont was 91%. The observed differences likely reflect adaptations to their specific habitats. For example, CTD079 encodes additional enzymes like the formate dehydrogenase increasing the organism's spectrum of energy generation pathways. Other additional physiological features of CTD079 included the increase of viral defense strategies, secretion systems and specific transporters for essential elements. These important genome characteristics suggest an adaptation to life in the open ocean.

Description: Eutrophication and climate change will affect habitats of species and more generally, the structure and functioning of ecosystems. We used a three‐dimensional, coupled hydrodynamic‐biogeochemical model to investigate potential future changes in size and location of potential habitats of marine species during the 21st century in a large, eutrophicated brackish sea (the Baltic Sea, northern Europe). We conducted scenario projections under the combined impact of nutrient load and climate change. Possible future changes of the eutrophication state of this sea were also assessed through two policy‐relevant indicators. The results imply a physiologically more stressful environment for marine species by the end of the 21st century: volumes of higher salinity water become more hypoxic/anoxic and the volumes of low salinity, oxic water increase. For example, these results impact and reduce cod reproductive habitats. The decrease is mainly climate change induced in the Baltic basins less directly influenced by inflows of saline, oxic water to the Baltic Sea (E Gotland and Gdansk Basins). In basins more directly influenced by such inflows (Arkona and Bornholm Basins), the combined effect from climate change and nutrient loads is of importance. The results for the eutrophication state indicators clearly indicate a more eutrophic sea than at present without a rigorous nutrient reduction policy, that is, the necessity to implement the Baltic Sea Action Plan. The multidisciplinary, multiscenario assessment strategy presented here provides a useful concept for the evaluation of impacts from cumulative stresses of changing climate and socioeconomic pressures on future eutrophication indicators and habitats of marine species.

Description: There is a strong economic interest in commercial deep‐sea mining of polymetallic nodules and therefore a need to define suitable preservation zones in the abyssal plain of the Clarion Clipperton Fracture Zone (CCZ). However, besides ship‐based multibeam data, only sparse continuous environmental information is available over large geographic scales. We test the potential of modelling meiofauna abundance and diversity on high taxonomic level on large geographic scale using a random forest approach. Ship‐based multibeam bathymetry and backscatter signal are the only sources for 11 predictor variables, as well as the modelled abundance of polymetallic nodules on the seafloor. Continuous meiofauna predictions have been combined with all available environmental variables and classified into classes representing abyssal habitats using k‐means clustering. Results show that ship‐based, multibeam‐derived predictors can be used to calculate predictive models for meiofauna distribution on a large geographic scale. Predicted distribution varies between the different meiofauna response variables. To evaluate predictions, random forest regressions were additionally computed with 1,000 replicates, integrating varying numbers of sampling positions and parallel samples per site. Higher numbers of parallel samples are especially useful to smoothen the influence of the remarkable variability of meiofauna distribution on a small scale. However, a high number of sampling positions is even more important, integrating a greater amount of natural variability of environmental conditions into the model. Synthesis and applications. Polymetallic nodule exploration contractors are required to define potential mining and preservation zones within their licence area. The biodiversity and the environment of preservation zones should be representative of the sites that will be impacted by mining. Our predicted distributions of meiofauna and the derived habitat maps are an essential first step to enable the identification of areas with similar ecological conditions. In this way, it is possible to define preservation zones not only based on expert opinion and environmental proxies but also integrating evidence from the distribution of benthic communities.

Description: Equatorial deep jets (EDJ) are vertically stacked, downward propagating zonal currents that alternate in direction with depth. In the tropical Atlantic, they have been shown to influence both surface conditions and tracer variability. Despite their importance, the EDJ are absent in most ocean models. Here we show that EDJ can be generated in an idealized ocean model when the model is driven only by the convergence of the meridional flux of intraseasonal zonal momentum diagnosed from a companion model run driven by steady wind forcing, corroborating the recent theory that intraseasonal momentum flux convergence maintains the EDJ. Additionally, the EDJ in our model nonlinearly generate mean zonal currents at intermediate depths that show similarities in structure to the observed circulation in the deep equatorial Atlantic, indicating their importance for simulating the tropical ocean mean state.

Description: Mesoscale eddies can be strengthened by the absorption of submesoscale eddies resulting from mixed-layer baroclinic instabilities. This is shown for mesoscale eddies in the Agulhas Current system by investigating the kinetic energy cascade with a spectral and a coarse-graining approach in two model simulations of the Agulhas region. One simulation resolves mixed-layer baroclinic instabilities and one does not. When mixed-layer baroclinic instabilities are included, the largest submesoscale near-surface fluxes occur in winter-time in regions of strong mesoscale activity for upscale as well as downscale directions. The forward cascade at the smallest resolved scales occurs mainly in frontogenetic regions in the upper 30 m of the water column. In the Agulhas ring path, the forward cascade changes to an inverse cascade at a typical scale of mixed-layer eddies (15 km). At the same scale, the largest sources of the upscale flux occur. After the winter, the maximum of the upscale flux shifts to larger scales. Depending on the region, the kinetic energy reaches the mesoscales in spring or early summer aligned with the maximum of mesoscale kinetic energy. This indicates the importance of submesoscale flows for the mesoscale seasonal cycle. A case study shows that the underlying process is the mesoscale absorption of mixed-layer eddies. When mixed-layer baroclinic instabilities are not included in the simulation, the open-ocean upscale cascade in the Agulhas ring path is almost absent. This contributes to a 20 %-reduction of surface kinetic energy at mesoscales larger than 100 km when submesoscale dynamics are not resolved by the model.