ALBERT

Description: Ice calved from the Antarctic and Greenland Ice Sheets or tidewater glaciers ultimately melts in the ocean contributing to sea-level rise. Icebergs have also been described as biological hotspots due to their potential roles as platforms for marine mammals and birds, and as micronutrient fertilizing agents. Icebergs may be especially important in the Southern Ocean where availability of the micronutrients iron and manganese extensively limits marine primary production. Whilst icebergs have long been described as a source of iron to the ocean, their nutrient signature is poorly constrained and it is unclear if there are regional differences. Here we show that 589 ice fragments collected from floating ice in contrasting regions spanning the Antarctic Peninsula, Greenland, and smaller tidewater systems in Svalbard, Patagonia and Iceland have similar characteristic (micro)nutrient signatures with limited or no significant differences between regions. Icebergs are a minor or negligible source of macronutrients to the ocean with low concentrations of NOx (NO3 + NO2, median 0.51 µM), PO4 (median 0.04 µM), and dissolved Si (dSi, median 0.02 µM). In contrast, icebergs deliver elevated concentrations of dissolved Fe (dFe; mean 82 nM, median 12 nM) and Mn (dMn; mean 26 nM, median 2.6 nM). A tight correlation between total dissolvable Fe and Mn (R2 = 0.95) and a Mn:Fe ratio of 0.024 suggested a lithogenic origin for the majority of sediment present in ice. Total dissolvable Fe and Mn retained a strong relationship with sediment load (both R2 = 0.43, p〈0.001), whereas weaker relationships were observed for dFe, dMn and dSi. Sediment load for Antarctic ice (median 9 mg L-1, n=144) was low compared to prior reported values for the Arctic. A particularly curious incidental finding was that melting samples of ice were observed to rapidly lose their sediment load, even when sediment layers were embedded within the ice and stored in the dark. Our results demonstrated that the nutrient signature of icebergs is consistent with an atmospheric source of NOx and PO4. Conversely, high Fe and Mn, and modest dSi concentrations, are associated with englacial sediment, which experiences limited biogeochemical processing prior to release into the ocean.

Description: Artificial upwelling has been discussed as a nature-based solution to fertilize currently unproductive areas of the ocean to enhance food web productivity and atmospheric CO2 sequestration. The efficacy of this approach may be closely tied to the nutrient stoichiometry of the upwelled water, as Si-rich water upwelling should benefit the growth of diatoms, who are key players for primary production, carbon export and food web efficiency. With a mesocosm experiment in subtropical waters, we assessed the physiological and functional responses of an oligotrophic phytoplankton community to artificial upwelling under varying Si:N ratios (0.07-1.33). Deep water fertilization led to strongly enhanced primary productivity rates and net autotrophy across Si scenarios. At the community level, Si-rich upwelling temporarily increased primary production and consistently enhanced diatom growth, producing up to 10-fold higher abundances compared to Si-deficient upwelling. At the organism level, contrasting effects were observed. On the one hand, silicification and size of diatom cells remained unaffected by Si:N, which is surprising given the direct dependency of these traits on Si. On the other hand, diatom Chlorophyll a density and carbon density were strongly reduced and particulate matter C:N was elevated under Si-rich upwelling. This suggests a reduced nutritional value for higher trophic levels under high Si:N ratios. Despite these strong qualitative changes under high Si, diatom cells appeared healthy and showed high photosynthetic efficiency. Our findings reveal great physiological plasticity and adaptability in phytoplankton under artificial upwelling, with Si-dependent trade-offs between primary producer quantity and quality.

Description: Coastal upwelling areas are extraordinarily productive environments where prokaryotic communities, the principal remineralizers of dissolved organic matter (DOM), rapidly respond to phytoplankton bloom and decay dynamics. Nevertheless, the extent of variability of key microbial functions in such dynamic waters remains largely unconstrained. Our metatranscriptomics analyses of 162 marker genes encoding ecologically relevant prokaryotic functions showed distinct spatial-temporal patterns in the NW Iberian Peninsula upwelling area. Short-term (daily) changes in specific bacterial functions associated with changes in biotic and abiotic factors were superimposed on seasonal variability. Taxonomic and functional specialization of prokaryotic communities, based mostly on different resource acquisition strategies, was observed. Our results uncovered the potential influence of prokaryotic functioning on phytoplankton bloom composition and development (e.g., Cellvibrionales and Flavobacteriales increased relative gene expression related to vitamin B12 and siderophore metabolisms during Chaetoceros and Dinophyceae summer blooms). Notably, bacterial adjustments to C- or N-limitation and DMSP availability during summer phytoplankton blooms and different spatial-temporal patterns of variability in the expression of genes with different phosphate affinity indicated a complex role of resource availability in structuring bacterial communities in this upwelling system. Also, a crucial role of Cellvibrionales in the degradation of DOM (carbohydrate metabolism, TCA cycle, proteorhodopsin, ammonium, and phosphate uptake genes) during the summer phytoplankton bloom was found. Overall, this dataset revealed an intertwined mosaic of microbial interactions and nutrient utilization patterns along a spatial-temporal gradient that needs to be considered if we aim to understand the biogeochemical processes in some of the most productive ecosystems in the world´s oceans.

Description: The element silicon is everywhere! In fact, silicon is the second most abundant element in Earth’s crust. Silicon in rocks and minerals breaks down and is transported from rivers and streams into the world’s oceans. Many marine organisms need silicon as it is a crucial nutrient to build their skeletons. Silicon eventually reaches the seafloor, but its journey into the abyss is not straightforward due to biological, physical, and chemical processes. All these processes transport and transform silicon, creating a cycle that we call the marine silicon cycle. The silicon cycle is directly connected to the carbon cycle, making silicon a key player in the regulation of Earth’s climate. In this article, we discuss why we need to understand the marine silicon cycle, explain the steps that happen in the ocean, and demonstrate how the marine silicon cycle affects humans.

Description: Silicon is a crucial nutrient that can join with the element oxygen to form a substance commonly called silica. Silica, commonly known as glass, is found in rocks in the Earth’s crust and dissolves into the oceans, where organisms like algae and sponges use it to build their glassy skeletons. This process, called biosilicification, is extremely important in the silica cycle. Over time, organisms have changed the silica cycle. Today, because of these organisms, the oceans no longer contain much silica. However, when the Earth was younger and these organisms had not evolved yet, no biological processes affected silica in the oceans. The evolution of these oceanic organisms across time has removed silica from the oceans. In this article, we discuss how the evolution of silicon-using sponges, as well as tiny organisms called zooplankton and algae, have changed the amount of silica in the world’s oceans through geologic time.

Description: Riverine nutrient export is an important process in marine coastal biogeochemistry and also impacts global marine biology. The nitrogen cycle is a key player here. Internal feedbacks regulate not only nitrogen distribution, but also primary production and thereby oxygen concentrations. Phosphorus is another essential nutrient and interacts with the nitrogen cycle via different feedback mechanisms. After a previous study of the marine nitrogen cycle response to riverine nitrogen supply, we here additionally include phosphorus from river export with different phosphorus burial scenarios and study the impact of phosphorus alone and in combination with nitrogen in a global 3-D ocean biogeochemistry model. Again, we analyse the effects on near coastal and open ocean biogeochemistry. We find that the addition of bio-available riverine phosphorus alone or together with nitrogen affects marine biology on millennial timescales more than riverine nitrogen alone. Biogeochemical feedbacks in the marine nitrogen cycle are strongly influenced by the additional phosphorus. Where bio-available phosphorus is increased by river input, nitrogen concentrations increase as well, except for regions with high denitrification rates. High phosphorus burial rates decrease biological production significantly. Globally, riverine phosphorus leads to elevated primary production rates in the coastal and open oceans.

Description: Ocean alkalinity enhancement (OAE) stands as a promising carbon dioxide removal technology. Yet, this solution to climate change entails shifts in water chemistry with unknown consequences for marine fish that are critical to ecosystem health and food security. With a laboratory and mesocosm experiment, we show that early life stages of fish can be resistant to OAE. We examined metabolic rate, swimming behavior, growth and survival in Atlantic herring (Clupea harengus) and other temperate coastal fish species. Neither direct physiological nor indirect food web-mediated impacts of OAE were apparent. This was despite non-CO2-equilibrated OAE (ΔTA = +600 µmol kg-1) that induces strong perturbations (ΔpH = +0.7, pCO2 = 75 µatm) compared to alternative deployment scenarios. Whilst our results give cause for optimism regarding the large-scale application of OAE, other life history stages (embryos) and habitats (open ocean) may prove more vulnerable. Still, our study across ecological scales (organism to community) and exposure times (short- to long-term) suggests that some fish populations, including key fisheries species, may be resilient to the carbonate chemistry changes under OAE.

Description: Octocorals (Cnidaria: Anthozoa) have a global distribution and form benthic assemblages along the depth gradient, from shallow to deep waters. They often occur below SCUBA diving limits, where they can become dominant habitat builders and aggregate different taxa. During a cruise in February 2023, one octocoral specimen was collected at 1453 m depth at Kebrit Deep, in the northern Saudi Arabian Red Sea axis, an area with extremely high temperature and salinity profiles at depth. Morphological analysis coupled with DNA barcoding using two mitochondrial markers ( COI and mtMuts ), revealed that the coral belongs to Acanthogorgia , a genus of azooxanthellate octocorals known to occur from 3 to 2300 m depths in cold, temperate and tropical waters. In the Red Sea, the genus was previously only known from shallower waters. Hence, we report the deepest record of the genus Acanthogorgia from the warm and saline Red Sea basin. This finding provides novel insights on deep-water octocoral diversity in the Red Sea, a still scantily explored area of the world, while emphasizing the need for further explorations at depth.

Description: Recent studies have begun to explore the potential of enhanced benthic weathering (EBW) in the Baltic Sea as a measure for climate change mitigation. To augment the understanding of EBW under seasonally changing conditions, this study aims to investigate weathering processes under anoxia to hypoxia in corrosive bottom waters, which reflect late summer conditions in the Baltic Sea. Dunite and calcite were added to sediment cores retrieved from Eckernförde Bay (Western Baltic Sea) with a constant flow-through of deoxygenated, CO 2 -enriched Baltic Sea bottom water. The addition of both materials increased benthic alkalinity release by 2.94 μmol cm −2 d −1 (calcite) and 1.12 μmol cm −2 d −1 (dunite), compared to the unamended control experiment. These excess fluxes are significantly higher than those obtained under winter conditions. The comparison with bottom water oxygen concentrations emphasizes that highest fluxes of alkalinity were associated with anoxic phases of the experiment. An increase in Ca and Si fluxes showed that the enhanced alkalinity fluxes could be attributed to calcite and dunite weathering. First order rate constants calculated based on these data were close to rates published in previous studies conducted under different conditions. This highlights the suitability of these proxies for mineral dissolution and justifies the use of these rate constants in modeling studies investigating EBW in the Baltic Sea and areas with similar chemical conditions. Generally stable pH profiles over the course of the experiment, together with the fact that the added minerals remained on the sediment surface, suggest that corrosive bottom waters were the main driving factor for the dissolution of the added minerals. These factors have important implications for the choice of mineral and timing for EBW as a possible marine carbon dioxide removal method in seasonally hypoxic to anoxic regions of the Baltic Sea.

Description: Rhodaliids, a specific family of siphonophores inhabiting the benthic zone, have remained enigmatic due to their rarity and elusive nature. These unique organisms, primarily found in open ocean habitats, exhibit distinctive features and are characterized by complex structures. During the Red Sea Decade Expedition, two rhodaliid specimens were collected at the sea bed at water depths of 438 and 495 meters. Regardless of challenges in specimen preservation, detailed morphological analysis revealed unique characteristics. Genetic analysis, employing the 16S rRNA marker, revealed one specimen closely related to Thermopalia taraxaca , while the taxonomic traits of this specimen suggested the identification of Archangelopsis jagoa . While most of the morphological features of the second specimen were also very similar to A. jagoa , some differences in coloration suggest the discovery of either an unusual colour variant or a potential new species within this genus. Notably, this study reports the first sequencing of A. jagoa . This species demonstrated a capacity to thrive in low-oxygen environments, challenging conventional assumptions about their habitat requirements. Despite difficulties in specimen handling and genetic analysis limitations due to a lack of comprehensive data, this research sheds light on the elusive world of benthic rhodaliids.

Description: Ocean alkalinity enhancement (OAE) is considered for the long-term removal of gigatons of carbon dioxide (CO2) from the atmosphere to achieve our climate goals. Little is known, however, about the ecosystem-level changes in biogeochemical functioning that may result from the chemical sequestration of CO2 in seawater, and how stable the sequestration is. We studied these two aspects in natural plankton communities under carbonate-based, CO2-equilibrated OAE in the nutrient-poor North Atlantic. During a month-long mesocosm experiment, the majority of biogeochemical pools, including inorganic nutrients, particulate organic carbon and phosphorus as well as biogenic silica, remained unaltered across all OAE levels of up to a doubling of ambient alkalinity (+2400 µeq kg-1). Noticeable exceptions were a minor decrease in particulate organic nitrogen and an increase in the carbon to nitrogen ratio (C:N) of particulate organic matter in response to OAE. Thus, in our nitrogen limited system, nitrogen turnover processes appear more susceptible than those of other elements leading to decreased food quality and increased organic carbon storage. However, alkalinity and chemical CO2 sequestration were not stable at all levels of OAE. Two weeks after alkalinity addition, we measured a loss of added alkalinity and of the initially stored CO2 in the mesocosm where alkalinity was highest (+2400 µeq kg-1, Ωaragonite ~10). The loss rate accelerated over time. Additional tests showed that such secondary precipitation can be initiated by particles acting as precipitation nuclei and that this process can occur even at lower levels of OAE. In conclusion, on the one hand, our study under carbonate-based OAE where the carbon is already sequestered, the risk of major and sustained impacts on biogeochemical functioning may be low in the nutrient-poor ocean. On the other hand, the durability of carbon sequestration using OAE could be constrained by alkalinity loss in supersaturated waters with precipitation nuclei present. Our study provides evaluation of ecosystem impacts of an idealised OAE deployment for monitoring, reporting and verification (MRV) in an oligotrophic system. Whether biogeochemical functioning is resilient to more technically simple and economically more viable approaches that induce stronger water chemistry perturbations remains to be seen.

Description: Circulation anomalies accompanying Sudden Stratospheric Warmings (SSWs) can have a significant impact on the troposphere. This surface response is observed for some but not all SSWs, and their downward coupling is not fully understood. We use an existing classification method to separate downward- and non-propagating SSWs (d/nSSWs) in ERA5 reanalysis data for the years 1979–2019. The differences in SSW downward propagation in composites of spatial patterns clearly show that dSSWs dominate the surface regional impacts following SSWs. During dSSWs, the upper-tropospheric jet stream is significantly displaced equatorward. Wave activity analysis shows remarkable differences between d/nSSWs for planetary and synoptic-scale waves. Enhanced stratospheric planetary eddy kinetic energy (EKE) and heat fluxes around the central date of dSSWs are followed by increased synoptic-scale wave activity and even surface coupling for synoptic-scale EKE. An observed significant reduction in upper-tropospheric synoptic-scale momentum fluxes following dSSWs confirms the important role of tropospheric eddy feedbacks for coupling to the surface. Our findings emphasize the role of the lower stratosphere and synoptic-scale waves in coupling the SSW signal to the surface and agree with mechanisms suggested in earlier modeling studies.

Description: Ocean alkalinity enhancement (OAE) is considered one of the most promising approaches to actively remove carbon dioxide (CO2) from the atmosphere by accelerating the natural process of rock weathering. This approach involves introducing alkaline substances sourced from natural mineral deposits such as olivine, basalt, and carbonates or obtained from industrial waste products such as steel slags, into seawater and dispersing them over coastal areas. Some of these natural and industrial substances contain trace metals, which would be released into the oceans along with the alkalinity enhancement. The trace metals could serve as micronutrients for marine organisms at low concentrations, but could potentially become toxic at high concentrations, adversely affecting marine biota. To comprehensively assess the feasibility of OAE, it is crucial to understand how the phytoplankton, which forms the base of marine food webs, responds to ocean alkalinization and associated trace metal perturbations. In this study, we investigated the toxicity of nickel on three representative phytoplankton species across a range of Ni concentrations (from 0 to 100 µmol L-1 with 12 µmol L-1 synthetic organic ligand). The results showed that the growth of the tested species was impacted differently. The low growth inhibition and high IC50 (concentration to inhibit growth rate by 50 %) revealed that both the coccolithophore Emiliania huxleyi and the dinoflagellate Amphidinium carterae were mildly impacted by the increase in Ni concentrations while the rapid response to exposure of Ni, high growth rate inhibition, and low IC50 of Thalassiosira weissflogii indicate low tolerance to Ni in this species. In conclusion, the variability in phytoplankton sensitivity to Ni suggests that for OAE applications with Ni-rich materials caution is required and critical toxic thresholds for Ni must be avoided.

Description: The central Arctic Ocean (CAO) plays an important role in the global carbon cycle, but the current and future exchange of the climate-forcing trace gases methane (CH4) and carbon dioxide (CO2) between the CAO and the atmosphere is highly uncertain. In particular, there are very few observations of near-surface gas concentrations or direct air–sea CO2 flux estimates and no previously reported direct air–sea CH4 flux estimates from the CAO. Furthermore, the effect of sea ice on the exchange is not well understood. We present direct measurements of the air–sea flux of CH4 and CO2, as well as air–snow fluxes of CO2 in the summertime CAO north of 82.5∘ N from the Synoptic Arctic Survey (SAS) expedition carried out on the Swedish icebreaker Oden in 2021. Measurements of air–sea CH4 and CO2 flux were made using floating chambers deployed in leads accessed from sea ice and from the side of Oden, and air–snow fluxes were determined from chambers deployed on sea ice. Gas transfer velocities determined from fluxes and surface-water-dissolved gas concentrations exhibited a weaker wind speed dependence than existing parameterisations, with a median sea-ice lead gas transfer rate of 2.5 cm h−1 applicable over the observed 10 m wind speed range (1–11 m s−1). The average observed air–sea CO2 flux was −7.6 ..., and the average air–snow CO2 flux was −1.1 . Extrapolating these fluxes and the corresponding sea-ice concentrations gives an August and September flux for the CAO of −1.75 ... , within the range of previous indirect estimates. The average observed air–sea CH4 flux of 3.5 ..., accounting for sea-ice concentration, equates to an August and September CAO flux of 0.35 , lower than previous estimates and implying that the CAO is a very small (≪ 1 %) contributor to the Arctic flux of CH4 to the atmosphere.

Description: Identification of seismically active fault zones and the definition of sufficiently large respect distances from these faults which enable avoiding the damaged rock zone surrounding the ruptured ground commonly are amongst the first steps to take in the geoscientific evaluation of sites suitable for nuclear waste disposal. In this work we present a GIS-based approach, using the earthquake-epicentre locations from the instrumental earthquake record of South-Korea to identify potentially active fault zones in the country, and compare different strategies for fault zone buffer creation as originally developed for site search in the high seismicity country Japan, and the low-to-moderate seismicity countries Germany and Sweden. In order to characterize the hazard potential of the Korean fault zones, we moreover conducted slip tendency analysis, here for the first time covering the fault zones of the entire Korean Peninsula. For our analyses we used the geo-spatial information from a new version of the Geological map of South-Korea, containing the outlines of 11 rock units, which we simplified to distinguish between 4 different rock types (granites, metamorphic rocks, sedimentary rocks and igneous rocks) and the surface traces of 1,528 fault zones and 6,654 lineaments identified through years of field work and data processing, a rich geo-dataset which we will publish along with this manuscript. Our approach for identification of active fault zones was developed without prior knowledge of already known seismically active fault zones, and as a proof of concept the results later were compared to a map containing already identified active fault zones. The comparison revealed that our approach identified 16 of the 21 known seismically active faults and added 472 previously unknown potentially active faults. The 5 seismically active fault zones which were not identified by our approach are located in the NE- and SW-sectors of the Korean Peninsula, which haven’t seen much recent seismic activity, and thus are not sufficiently well covered by the seismic record. The strike directions of fault zones identified as active are in good agreement with the orientation of the current stress field of the peninsula and slip tendency analysis provided first insights into subsurface geometry such as the dip angles of both active and inactive fault zones. The results of our work are of major importance for the early-stage seismic hazard assessment that has to be conducted in support of the nuclear waste disposal siting in South-Korea. Moreover, the GIS-based methods for identification of active fault zones and buffering of respect areas around fault zone traces presented here, are applicable also elsewhere.

Description: Physical and chemical trace metal speciation are important for our understanding of metal cycling and potential toxicity to marine life. Trace metals can behave differently in diffusion processes or particle-solution interactions and have different bioavailabilities depending on their physical and chemical forms, which often depend on redox conditions. Here we investigated dissolved (〈 0.2 µm) and soluble (〈 0.02 µm) concentrations of Mn, Co, Ni, Fe, Cu, V, Mo, U, Cd, and As in oxic and suboxic deep-sea sediments of the central equatorial Pacific Ocean. Vanadium, Mo, U, As, and Cd showed no significant concentration differences between their dissolved and soluble forms, suggesting that they are present as inorganic ionic species or organic complexes in the truly dissolved or small colloidal fraction. In contrast, the colloidal fraction (〉 0.02 µm 〈 0.2 µm) of Mn, Co, Ni, and Cu increased with depth in oxic pore waters and Fe had the largest but variable colloidal pool. Soluble Mn, Co, and Ni were released in the uppermost 2-4 cm in the sediment because of reductive dissolution. The increasing colloidal fraction with depth suggests a decrease in the concentration of small organic ligands with depth, that are abundant in the surface sediment pore waters, and instead an increasing importance of larger (〉 0.02 µm) inorganic nanoparticles and colloids such as Mn and Fe (oxyhydr)oxides that control Mn, Fe, and Co cycling at depths 〉 10 cm. The distribution of Ni and Cu cannot be exclusively explained by inorganic nanoparticles and a shift from low to larger high molecular weight organic ligands might occur. These findings provide new insights into trace metal distributions in the dissolved phase, highlighting the diversity of metal complexes and the need to incorporate these in future calculations of benthic metal fluxes and ecotoxicity assessments, especially in oxic pore waters.

Description: Within the context of climate change, coastal vegetated ecosystems have the capacity for long-term carbon storage. Blue carbon refers to such carbon trapped in the oceans and coastal shelf seas. These ecosystems are under anthropogenic pressure and, to help these ecosystems to thrive and realize their carbon storage potentials, interventions require acceptance from society, in general, and adjacent coastal communities, in particular. Through a random street survey along the German coasts in 2022, quantitative and qualitative data were collected from more than 200 participants. A questionnaire comprising 50 open and closed questions was designed to assess the status quo of German coastal residents’ norms and values concerning blue carbon ecosystems. Focus was put on nature conservation and climate change perceptions. The survey results reveal that most residents along the German coast valued nature conservation while idealizing nature that is seen as “untouched” by humans. Responses regarding active interventions to improve coastal ecosystem services were diverse. Blue carbon strategies are likely to operate within this area of tension. Most respondents were aware of climate change as a threat to their home region and were in favor of an increase in action against climate change there. The respondents were familiar with CO 2 reduction and avoidance strategies. However, they were less aware of measures to remove atmospheric CO 2 and the potential of storing CO 2 in ecosystems beyond afforestation measures. Due to a lack of knowledge, no consolidated public opinions on blue carbon in coastal vegetated ecosystems could be identified, blurring societal acceptance of blue carbon strategies. While these ecosystems are particularly vulnerable to human disturbance, long-term carbon storage is essential for blue carbon. Therefore, the individual acceptance of interventions from people living in close proximity to intervention sites is key for sustained success. The present article concludes that there are possibilities to co-create knowledge and acceptance as prerequisites for blue carbon interventions to possibly become efficacious.

Description: Metaorganism research contributes substantially to our understanding of the interaction between microbes and their hosts, as well as their co-evolution. Most research is currently focused on the bacterial community, while archaea often remain at the sidelines of metaorganism-related research. Here, we describe the archaeome of a total of eleven classical and emerging multicellular model organisms across the phylogenetic tree of life. To determine the microbial community composition of each host, we utilized a combination of archaea and bacteria-specific 16S rRNA gene amplicons. Members of the two prokaryotic domains were described regarding their community composition, diversity, and richness in each multicellular host. Moreover, association with specific hosts and possible interaction partners between the bacterial and archaeal communities were determined for the marine models. Our data show that the archaeome in marine hosts predominantly consists of Nitrosopumilaceae and Nanoarchaeota, which represent keystone taxa among the porifera. The presence of an archaeome in the terrestrial hosts varies substantially. With respect to abundant archaeal taxa, they harbor a higher proportion of methanoarchaea over the aquatic environment. We find that the archaeal community is much less diverse than its bacterial counterpart. Archaeal amplicon sequence variants are usually host-specific, suggesting adaptation through co-evolution with the host. While bacterial richness was higher in the aquatic than the terrestrial hosts, a significant difference in diversity and richness between these groups could not be observed in the archaeal dataset. Our data show a large proportion of unclassifiable archaeal taxa, highlighting the need for improved cultivation efforts and expanded databases.

Description: Total alkalinity (AT) and dissolved inorganic carbon (CT) in the oceans are important properties with respect to understanding the ocean carbon cycle and its link to global change (ocean carbon sinks and sources, ocean acidification) and ultimately finding carbon-based solutions or mitigation procedures (marine carbon removal). We present a database of more than 44 400 AT and CT observations along with basic ancillary data (spatiotemporal location, depth, temperature and salinity) from various ocean regions obtained, mainly in the framework of French projects, since 1993. This includes both surface and water column data acquired in the open ocean, coastal zones and in the Mediterranean Sea and either from time series or dedicated one-off cruises. Most AT and CT data in this synthesis were measured from discrete samples using the same closed-cell potentiometric titration calibrated with Certified Reference Material, with an overall accuracy of ±4 µmol kg−1 for both AT and CT. The data are provided in two separate datasets – for the Global Ocean and the Mediterranean Sea (https://doi.org/10.17882/95414, Metzl et al., 2023), respectively – that offer a direct use for regional or global purposes, e.g., AT–salinity relationships, long-term CT estimates, and constraint and validation of diagnostic CT and AT reconstructed fields or ocean carbon and coupled climate–carbon models simulations as well as data derived from Biogeochemical-Argo (BGC-Argo) floats. When associated with other properties, these data can also be used to calculate pH, the fugacity of CO2 (fCO2) and other carbon system properties to derive ocean acidification rates or air–sea CO2 fluxes.

Description: The upper wind-driven circulation in the tropical Atlantic Ocean plays a key role in the basin-wide distribution of water mass properties and affects the transport of heat, freshwater, and biogeochemical tracers such as oxygen or nutrients. It is crucial to improve our understanding of its long-term behaviour, which largely relies on model simulations and applied forcing due to sparse observational data coverage, especially before the mid-2000s. Here, we apply two different forcing products, the Coordinated Ocean-ice Reference Experiments (CORE) v2 and the Japanese 55-year Reanalysis (JRA55-do) surface dataset, to a high-resolution ocean model. Where possible, we compare the simulated results to long-term observations. We find large discrepancies between the two simulations regarding the wind and current field. In the CORE simulation, strong, large-scale wind stress curl amplitudes above the upwelling regions of the eastern tropical North Atlantic seem to cause an overestimation of the mean and seasonal variability in the eastward subsurface current just north of the Equator. The wind stress curl of JRA55-do forcing shows much finer structures, and the JRA55-do simulation is in better agreement with the mean and intraseasonal fluctuations in the subsurface current found in observations. The northern branch of the South Equatorial Current flows westward at the surface just north of the Equator. On interannual to decadal timescales, it shows a high correlation of R=0.9 with the zonal wind stress in the CORE simulation but only a weak correlation of R=0.35 in the JRA55-do simulation. We also identify similarities between the two simulations. The strength of the eastward-flowing North Equatorial Counter Current located between 3 and 10° N covaries with the strength of the meridional wind stress just north of the Equator on interannual to decadal timescales in the two simulations. Both simulations present a comparable mean, seasonal cycle and trend of the eastward off-equatorial subsurface current south of the Equator but underestimate the current strength by half compared to observations. In both simulations, the eastward-flowing Equatorial Undercurrent weakened between 1990 and 2009. In the JRA simulation, which covers the modern period of observations, the Equatorial Undercurrent strengthened again between 2008 to 2018, which agrees with observations, although the simulation underestimates the strengthening by over a third. We propose that long-term observations, once they have reached a critical length, need to be used to test the quality of wind-driven simulations. This study presents one step in this direction.

Description: We conducted extensive sediment trap experiments in the Benguela Upwelling System (BUS) in the south-eastern Atlantic Ocean to study the influence of zooplankton on the flux of particulate organic carbon (POC) through the water column and its sedimentation. Two long term moored and sixteen short term free-floating sediment trap systems were deployed. The mooring experiments were conducted for several years and the sixteen drifters were deployed on three different research cruises between 2019 and 2021. Zooplankton was separated from the trapped material and divided into 8 different zooplankton groups. In contrast to zooplankton which actively carries POC into the traps in the form of biomass (active POC flux), the remaining fraction of the trapped material was assumed to fall passively into the traps along with sinking particles (passive POC flux). The results show, in line with other studies, that copepods dominate the active POC flux, with the active POC flux in the southern BUS (sBUS) being about three times higher than in the northern BUS (nBUS). In contrast, the differences between the passive POC fluxes in the nBUS and sBUS were small. Despite large variations, which reflected the variability within the two subsystems, the mean passive POC fluxes from the drifters and the moored traps could be described using a common POC flux attenuation equation. However, the almost equal passive POC flux, on the one hand, and large variations in the POC concentration in the surface sediments between the nBUS and sBUS, on the other hand, imply that factors others than the POC supply exert the main control on POC sedimentation in the BUS. The varying intensity of the near-bottom oxygen minimum zone (OMZ), which is more pronounced in the nBUS than in the sBUS, could in turn explain the differences in the sediments, as the lack of oxygen reduces the POC degradation. Hence, globally expanding OMZs might favour POC sedimentation in regions formerly exposed to oxygenated bottom water but bear the risk of increasing the frequency of anoxic events in the oxygen-poor upwelling systems. Apart from associated release of CH4, which is a much more potent greenhouse gas than CO2, such events pose a major threat to the pelagic ecosystem and fisheries.

Description: Nitrogen (N) is a crucial limiting nutrient for phytoplankton growth in the ocean. The main source of bioavailable N in the ocean is delivered by N2-fixing diazotrophs in the surface layer. Since field observation of N2 fixation are spatially and temporally sparse, the fundamental processes and mechanisms controlling N2 fixation are not well understood and constrained. Here, we implement benthic denitrification in an Earth System Model of intermediate complexity (UVic-ESCM 2.9) coupled to an optimality-based plankton ecosystem model (OPEM v1.1). Benthic denitrification occurs mostly in coastal upwelling regions and on shallow continental shelves, and is the largest N-loss process in the global ocean. We calibrate our model against three different combinations of observed Chl, NO3-, PO43-, O2 and N* = NO3- −16PO43- +2.9. The inclusion of N* provides a powerful constraint on biogeochemical model behavior. Our new model version including benthic denitrification simulates higher global rates of N2 fixation with a more realistic distribution extending to higher latitudes that are supported by independent estimates based on geochemical data. Oxygen deficient zone volume and water column denitrification rates are reduced in the new version, indicating that including benthic denitrification may improve global biogeochemical models that commonly overestimate anoxic zones. With the improved representation of the ocean N cycle, our new model configuration also yields better global net primary production (NPP) when compared to the independent datasets not included in the calibration. Benthic denitrification plays an important role shaping N2 fixation and NPP throughout the global ocean in our model, and should be considered when evaluating and predicting their response to environmental change.

Description: Since a pH sensor has become available that is principally suitable for use on demanding autonomous measurement platforms, the marine CO2 system can be observed independently and continuously by Biogeochemical Argo floats. This opens the potential to detect variability and long-term changes in interior ocean inorganic carbon storage and quantify the ocean sink for atmospheric CO2. In combination with a second parameter of the marine CO2 system, pH can be a useful tool to derive the surface ocean CO2 partial pressure (pCO2). The large spatiotemporal variability in the marine CO2 system requires sustained observations to decipher trends and study the impacts of short-term events (e.g., eddies, storms, phytoplankton blooms) but also puts a high emphasis on the quality control of float-based pH measurements. In consequence, a consistent and rigorous quality control procedure is being established to correct sensor offsets or drifts as the interpretation of changes depends on accurate data. By applying current standardized routines of the Argo data management to pH measurements from a pH / O2 float pilot array in the subpolar North Atlantic Ocean, we assess the uncertainties and lack of objective criteria associated with the standardized routines, notably the choice of the reference method for the pH correction (CANYON-B, LIR-pH, ESPER-NN, and ESPER-LIR) and the reference depth for this adjustment. For the studied float array, significant differences ranging between ca. 0.003 pH units and ca. 0.04 pH units are observed between the four reference methods which have been proposed to correct float pH data. Through comparison against discrete and underway pH data from other platforms, an assessment of the adjusted float pH data quality is presented. The results point out noticeable discrepancies near the surface of 〉 0.004 pH units. In the context of converting surface ocean pH measurements into pCO2 data for the purpose of deriving air–sea CO2 fluxes, we conclude that an accuracy requirement of 0.01 pH units (equivalent to a pCO2 accuracy of 10 µatm as a minimum requirement for potential future inclusion in the Surface Ocean CO2 Atlas, SOCAT, database) is not systematically achieved in the upper ocean. While the limited dataset and regional focus of our study do not allow for firm conclusions, the evidence presented still calls for the inclusion of an additional independent pH reference in the surface ocean in the quality control routines. We therefore propose a way forward to enhance the float pH quality control procedure. In our analysis, the current philosophy of pH data correction against climatological reference data at one single depth in the deep ocean appears insufficient to assure adequate data quality in the surface ocean. Ideally, an additional reference point should be taken at or near the surface where the resulting pCO2 data are of the highest importance to monitor the air–sea exchange of CO2 and would have the potential to very significantly augment the impact of the current observation network.

Description: Here we present a confocal Fe K-edge μ-XANES method (where XANES stands for X-ray absorption near-edge spectroscopy) for the analysis of Fe oxidation state in heterogeneous and one-side-polished samples. The new technique allows for an analysis of small volumes with high spatial 3D resolution of 〈100 µm3. The probed volume is restricted to that just beneath the surface of the exposed object. This protocol avoids contamination of the signal by the host material and minimizes self-absorption effects. This technique has been tested on a set of experimental glasses with a wide range of Fe3+  ΣFe ratios. The method was applied to the analysis of natural melt inclusions trapped in forsteritic to fayalitic olivine crystals of the Hekla volcano, Iceland. Our measurements reveal changes in Fe3+  ΣFe from 0.17 in basaltic up to 0.45 in dacitic melts, whereas the magnetite–ilmenite equilibrium shows redox conditions with Fe3+  ΣFe ≤0.20 (close to FMQ, fayalite–magnetite–quartz redox equilibrium) along the entire range of Hekla melt compositions. This discrepancy indicates that the oxidized nature of glasses in the melt inclusions could be related to the post-entrapment process of diffusive hydrogen loss from inclusions and associated oxidation of Fe in the melt. The Fe3+  ΣFe ratio in silicic melts is particularly susceptible to this process due to their low FeO content, and it should be critically evaluated before petrological interpretation.

Description: Sediment fluxes to the seafloor govern the fate of elements and compounds in the ocean and serve as a prerequisite for research on elemental cycling, benthic processes and sediment management strategies. To quantify these fluxes over seafloor areas, it is necessary to scale up sediment mass accumulation rates (MAR) obtained from multiple sample stations. Conventional methods for spatial upscaling involve averaging of data or spatial interpolation. However, these approaches may not be sufficiently precise to account for spatial variations of MAR, leading to poorly constrained regional sediment budgets. Here, we utilize a machine learning approach to scale up porosity and 210 Pb data from 145 and 65 stations, respectively, in the Skagerrak. The models predict the spatial distributions by considering several predictor variables that are assumed to control porosity and 210 Pb rain rates. The spatial distribution of MAR is based on the predicted porosity and existing sedimentation rate data. Our findings reveal highest MAR and 210 Pb rain rates to occur in two parallel belt structures that align with the general circulation pattern in the Skagerrak. While high 210 Pb rain rates occur in intermediate water depths, the belt of high MAR is situated closer to the coastlines due to lower porosities at shallow water depths. Based on the spatial distributions, we calculate a total MAR of 34.7 Mt yr -1 and a 210 Pb rain rate of 4.7 · 10 14 dpm yr -1 . By comparing atmospheric to total 210 Pb rain rates, we further estimate that 24% of the 210 Pb originates from the local atmospheric input, with the remaining 76% being transported laterally into the Skagerrak. The updated MAR in the Skagerrak is combined with literature data on other major sediment sources and sinks to present a tentative sediment budget for the North Sea, which reveals an imbalance with sediment outputs exceeding the inputs. Substantial uncertainties in the revised Skagerrak MAR and the literature data might close this imbalance. However, we further hypothesize that previous estimates of suspended sediment inputs into the North Sea might have been underestimated, considering recently revised and elevated estimates on coastal erosion rates in the surrounding region of the North Sea.

Description: We examine the impact of horizontal resolution and model time step on the climate of the OpenIFS version 43r3 atmospheric general circulation model. A series of simulations for the period 1979–2019 are conducted with various horizontal resolutions (i.e. ∼100, ∼50, and ∼25 km) while maintaining the same time step (i.e. 15 min) and using different time steps (i.e. 60, 30, and 15 min) at 100 km horizontal resolution. We find that the surface zonal wind bias is significantly reduced over certain regions such as the Southern Ocean and the Northern Hemisphere mid-latitudes and in tropical and subtropical regions at a high horizontal resolution (i.e. ∼25 km). Similar improvement is evident too when using a coarse-resolution model (∼100 km) with a smaller time step (i.e. 30 and 15 min). We also find improvements in Rossby wave amplitude and phase speed, as well as in weather regime patterns, when a smaller time step or higher horizontal resolution is used. The improvement in the wind bias when using the shorter time step is mostly due to an increase in shallow and mid-level convection that enhances vertical mixing in the lower troposphere. The enhanced mixing allows frictional effects to influence a deeper layer and reduces wind and wind speed throughout the troposphere. However, precipitation biases generally increase with higher horizontal resolutions or smaller time steps, whereas the surface air temperature bias exhibits a small improvement over North America and the eastern Eurasian continent. We argue that the bias improvement in the highest-horizontal-resolution (i.e. ∼25 km) configuration benefits from a combination of both the enhanced horizontal resolution and the shorter time step. In summary, we demonstrate that, by reducing the time step in the coarse-resolution (∼100 km) OpenIFS model, one can alleviate some climate biases at a lower cost than by increasing the horizontal resolution.

Description: Despite management efforts, anthropogenic nutrient enrichments continue to enhance phytoplankton blooms worldwide. Release of nitrogen and phosphorus compounds not only provides surplus of nutrients but also disbalances their stoichiometry. Declines in the relative availability of dissolved silicon might induce limitation in diatoms, major primary producers with silicified shells. We studied experimentally how nutrient enrichment and resulting decline in dissolved silicon to nitrogen ratios (Si:N) affect the structure and functioning of natural plankton communities. Nitrate was added to create a range of Si:N ratios and phosphate was supplied in Redfield ratio to nitrogen. We also manipulated copepod abundance to understand the top-down effects on communities experiencing nutrient enrichment. Nitrogen and phosphorus additions resulted in a steep phytoplankton biomass increase, followed by a post-bloom decline. Phytoplankton bloom biomass was higher in high nitrogen treatments but during the post-bloom period this trend switched. Biomass was sustained longer in high Si:N treatments, indicating that silicon limitation terminates the bloom. Many diatom species did not benefit from nitrogen and phosphorus enrichment and diatom dominance ceased below Si:N of 0.4:1. Under high grazing pressure, silicate was taken up faster suggesting that silicification is important in diatom defense. Copepods shaped plankton communities via feeding on dinoflagellates, chlorophytes and the diatom Skeletonema costatum but there was no significant effect of nitrogen and phosphorus enrichment on copepod abundance. Our results, combined with previous studies, show that while nutrient concentrations define the total phytoplankton bloom biomass, resource ratios are important in sustaining biomass and determining community structure and composition.

Description: Phytoplankton forms the base of the marine food web by transforming CO2 into organic carbon via photosynthesis. Some of the organic carbon is then transferred through the food web and exported into the deep ocean, a process known as the biological carbon pump. Despite the importance of phytoplankton for marine ecosystems and the global carbon cycle, projections of phytoplankton biomass in response to climate change differ strongly across Earth system models, illustrating uncertainty in our understanding of the underlying processes. Differences are especially large in the Southern Ocean, a region that is notoriously difficult to represent in models. Here, we argue that water column-integrated phytoplankton biomass in the Southern Ocean is projected to largely remain unchanged under climate change by the CMIP6 multi-model ensemble because of a shifting balance of bottom-up and top-down processes driven by a shoaling mixed layer depth. A shallower mixed layer is projected to improve growth conditions and consequently weaken bottom-up control. In addition to enhanced phytoplankton growth, the shoaling of the mixed layer also compresses phytoplankton closer to the surface and promotes zooplankton grazing efficiency, thus intensifying top-down control. Overall, our results suggest that while changes in bottom-up conditions stimulate enhanced growth, intensified top-down control opposes an increase in phytoplankton and becomes increasingly important for phytoplankton response under climate change in the Southern Ocean.

Description: Transparent exopolymer particles (TEP) and Coomassie stainable particles (CSP), two prominent classes of gel−like particles in the ocean primarily produced by phytoplankton, play crucial roles in ecological and biogeochemical processes, influencing microbial nutrition, growth, and particle aggregation. The distribution of these particles is intricately linked to the spatiotemporal dynamics of phytoplankton. Mesoscale cyclonic eddies (CEs) are known to stimulate phytoplankton growth and influence particle transport, but their effects on TEP and CSP remain to be determined. In the Eastern Tropical North Atlantic (ETNA), we examined three CEs: one off the Mauritanian coast during summer (Mau), one offshore during winter (Sal), and another near Brava island during winter. Mau and Brava CEs were in their intensification/maturity phase, while the Sal CE was in its decay phase. Both TEP and CSP concentrations correlated with primary productivity, but TEP increased with chlorophyll−a concentration, whereas elevated CSP coincided also with the highest abundance of pico−nanophytoplankton (〈20 µm), mainly Synechococcus. Both gels exhibited a positive correlation with bacterial biomass production, indicating their consumption by heterotrophic bacteria. TEP total area in the epipelagic waters of all CEs (Mau, Brava, and Sal) was elevated compared to surrounding waters, with on average 4, 2.5, and 1.6−fold higher values, respectively. However, no significant difference in TEP size distribution was observed within any CEs and their surroundings. Similarly, CSP total area increased in the epipelagic waters of Mau and Brava CEs, with on average 5 and 2.4−fold higher values, respectively, compared to surrounding waters. CSP particles were notably larger in these two eddies, while the Sal CE showed no significant difference from surrounding waters in CSP abundance and size. Overall, TEP and CSP exhibited distinct responses to CEs, with increased concentrations during their intensification/maturation stage and remineralization dominating during their decaying stage.

Description: Introduction Ecosystem engineers play a pivotal role in shaping habitats through their activities and presence. In shallow Baltic waters, seagrasses, patch-forming mussels, and infaunal clams modify soft bottom habitats, impacting benthic community structure. While the individual effects of these ecosystem engineers are well studied, interactions among co-occurring engineers are poorly understood.Methods We conducted a mesocosm experiment to assess the independent and combined impacts of seagrass (Zostera marina), epifaunal mussels (Mytilus spp.), and infaunal clams (Macoma balthica) on invertebrate colonization in soft sediments.Results Our findings reveal significant engineer-driven alterations in macrofaunal community structure. Combined engineer effects diverged from individual impacts, indicating potential synergies or antagonisms in sediment (re)colonization. Notably, a higher number of engineer species positively affected the diversity of settled macrofauna, with the lowest macrofaunal abundance and biomass but the highest Shannon diversity found in the presence of all three engineers.Discussion Results suggest that seagrass, mussels, and clams influence benthos through larval settlement and sediment biogeochemistry, providing insights into the distinct roles of habitat-forming organisms in shaping the benthic communities in coastal ecosystems of the Baltic Sea.

Description: Deoxygenation is tied to organic carbon (Corg) supply and utilization in marine systems. Under oxygen-depletion, bacteria maintain respiration using alternative electron acceptors such as nitrate. Since anaerobic respiration's energy yield is lower, Corg remineralization may be reduced and its residence time increased. We investigated the influence of oxygen and alternative electron acceptors' availability on Corg cycling by heterotrophic bacteria during a continuous culture experiment with Shewanella baltica, a facultative anaerobic γ-Proteobacteria in the Baltic Sea. We tested six different oxygen levels, from suboxic (〈5 µmol L-1 ) to fully oxic conditions, using media (salinity=14 g L-1 ) supplied with high (HighN) or low (LowN) inorganic nitrogen concentrations relative to glucose as labile Corg source. Our results show that suboxia limited DOC (glucose) uptake and cell growth only under LowN, while higher availability of alternative electron acceptors seemingly compensated oxygen limitation under HighN. N-loss was observed under suboxia in both nitrogen treatments. Under HighN, N-loss was highest and a C:N loss ratio of ~2.0 indicated that Corg was remineralized via denitrification. Under LowN, the C:N loss ratio under suboxia was higher (~5.5), suggesting dominance of other anaerobic respiration pathways, such as dissimilatory nitrate reduction to ammonium (DNRA). Bacterial growth efficiency was independent of oxygen concentration but higher under LowN (34±3.0%) than HighN (26±1.6%). Oxygen concentration also affected dissolved organic matter (DOM) cycling. Under oxic conditions, the release of dissolved combined carbohydrates was enhanced, and the amino acid-based degradation index (DI) pointed to more diagenetically altered DOM. Our results suggest bacterial Corg uptake in low-oxygen systems dominated by S. baltica can be limited by oxygen but compensated by high nitrate availability. Hence, suboxia diminishes Corg remineralisation only when alternative electron acceptors are lacking. Under high nitrate:Corg supply, denitrification leads to a higher N:C loss ratio, potentially counteracting eutrophication in the long run. Low nitrate:Corg supply may favour other anaerobic respiration pathways like DNRA, which sustains labile nitrogen in the system, potentially intensifying the cycle of eutrophication. Going forward, it will be crucial to establish the validity of our findings for S. baltica in natural systems with diverse organic substrates and microbial consortia.

Description: As a part of the Scientific Committee on Oceanographic Research (SCOR) Working Group #160 “Analyzing ocean turbulence observations to quantify mixing” (ATOMIX), we have developed recommendations on best practices for estimating the rate of dissipation of kinetic energy, ε, from measurements of turbulence shear using shear probes. The recommendations provided here are platform-independent and cover the conceivable range of dissipation rates in the ocean, seas, and other natural waters. They are applicable to commonly deployed platforms that include vertical profilers, fixed and moored instruments, towed profilers, submarines, self-propelled ocean gliders, and other autonomous underwater vehicles. The procedure for preparing the shear data for spectral estimation is discussed in detail, as are the quality control metrics that should accompany each estimate of ε. The methods are illustrated using a high-quality ‘benchmark’ dataset, while potential pitfalls are demonstrated with a second dataset containing common faults.

Description: The climate science community aims to improve our understanding of climate change due to anthropogenic influences on atmospheric composition and the Earth's surface. Yet not all climate interactions are fully understood and diversity in climate model experiments persists as assessed in the latest Intergovernmental Panel on Climate Change (IPCC) assessment report. This article synthesizes current challenges and emphasizes opportunities for advancing our understanding of climate change and model diversity. The perspective of this article is based on expert views from three multi-model intercomparison projects (MIPs) – the Precipitation Driver Response MIP (PDRMIP), the Aerosol and Chemistry MIP (AerChemMIP), and the Radiative Forcing MIP (RFMIP). While there are many shared interests and specialisms across the MIPs, they have their own scientific foci and specific approaches. The partial overlap between the MIPs proved useful for advancing the understanding of the perturbation-response paradigm through multi-model ensembles of Earth System Models of varying complexity. It specifically facilitated contributions to the research field through sharing knowledge on best practices for the design of model diagnostics and experimental strategies across MIP boundaries, e.g., for estimating effective radiative forcing. We discuss the challenges of gaining insights from highly complex models that have specific biases and provide guidance from our lessons learned. Promising ideas to overcome some long-standing challenges in the near future are kilometer-scale experiments to better simulate circulation-dependent processes where it is possible, and machine learning approaches for faster and better sub-grid scale parameterizations where they are needed. Both would improve our ability to adopt a smart experimental design with an optimal tradeoff between resolution, complexity and simulation length. Future experiments can be evaluated and improved with sophisticated methods that leverage multiple observational datasets, and thereby, help to advance the understanding of climate change and its impacts.

Description: We carried out measurements of the CO2 system parameters to evaluate the impact of carbonate and nutrients' chemistry on phytoplankton populations in the Gulf of Guinea (GoG). The seasonal variations of the CO2 system parameters (fCO2, DIC, pH and TA) along with nitrates and phosphates were quantified weekly at surface (between 0 and 5 m depth) (5.57 degrees N - 4.57 degrees W) in the GoG from May to December 2020. Seawater pH varied widely during the study period, ranging between 8.10-8.35 pH units; DIC and TA varied between 1810 and 2094 mu mol kg-1, and between 2051 and 2216 mu mol-1 respectively. DIC peaks coincided with the high upwelling period (August and September). For phytoplankton, a total of 60 species were found belonging to four taxonomic phyla: Bacillariophyta, Dinophyta, Chlorophyta and Dictyochophyta. The highest number of phytoplanktonic species were recorded for Bacillariophyta phylum with 36 species (60%). The phylum Dinophyta comprised 22 taxa (36%) and Chlorophyta and Dictyochophyta recorded only one species (2%). The highest specific diversities were observed in August and September with 29 and 26 taxa respectively and the lowest was found in October-November (5 taxa) and December (one taxa). Bacillariophyta and Dinophyta appeared throughout the entire study period. The only species for Chlorophyta phylum appeared in June and July and the Dictyochophyta's one in May, July and August. In general, the physical (SST, SSS) and chemical (TA, DIC, pH) parameters influenced less than 50% of the phytoplankton population in the coastal area of the GoG. Our study shows that Bacillariophyta population grows up when the physicochemical parameters' variability increase.

Description: The presented pilot for the Synthesis Product for Ocean Time Series (SPOTS) includes data from 12 fixed ship-based time-series programs. The related stations represent unique open-ocean and coastal marine environments within the Atlantic Ocean, Pacific Ocean, Mediterranean Sea, Nordic Seas, and Caribbean Sea. The focus of the pilot has been placed on biogeochemical essential ocean variables: dissolved oxygen, dissolved inorganic nutrients, inorganic carbon (pH, total alkalinity, dissolved inorganic carbon, and partial pressure of CO2), particulate matter, and dissolved organic carbon. The time series used include a variety of temporal res- olutions (monthly, seasonal, or irregular), time ranges (10–36 years), and bottom depths (80–6000 m), with the oldest samples dating back to 1983 and the most recent one corresponding to 2021. Besides having been harmo- nized into the same format (semantics, ancillary data, units), the data were subjected to a qualitative assessment in which the applied methods were evaluated and categorized. The most recently applied methods of the time- series programs usually follow the recommendations outlined by the Bermuda Time Series Workshop report (Lorenzoni and Benway, 2013), which is used as the main reference for “method recommendations by prevalent initiatives in the field”. However, measurements of dissolved oxygen and pH, in particular, still show room for improvement. Additional data quality descriptors include precision and accuracy estimates, indicators for data variability, and offsets compared to a reference and widely recognized data product for the global ocean: the GLobal Ocean Data Analysis Project (GLODAP). Generally, these descriptors indicate a high level of continuity in measurement quality within time-series programs and a good consistency with the GLODAP data product, even though robust comparisons to the latter are limited. The data are available as (i) a merged comma-separated file that is compliant with the World Ocean Circulation Experiment (WOCE) exchange format and (ii) a format dependent on user queries via the Environmental Research Division’s Data Access Program (ERDDAP) server of the Global Ocean Observing System (GOOS). The pilot increases the data utility, findability, accessibility, interoperability, and reusability following the FAIR philosophy, enhancing the readiness of biogeochemical time series. It facilitates a variety of applications that benefit from the collective value of biogeochemical time-series observations and forms the basis for a sustained time-series living data product, SPOTS, complementing relevant products for the global interior ocean carbon data (GLobal Ocean Data Analysis Project), global surface ocean carbon data (Surface Ocean CO2 Atlas; SOCAT), and global interior and surface methane and nitrous oxide data (MarinE MethanE and NiTrous Oxide product). Aside from the actual data compilation, the pilot project produced suggestions for reporting metadata, im- plementing quality control measures, and making estimations about uncertainty. These recommendations aim to encourage the community to adopt more consistent and uniform practices for analysis and reporting and to update these practices regularly. The detailed recommendations, links to the original time-series programs, the original data, their documentation, and related efforts are available on the SPOTS website. This site also pro- vides access to the data product (DOI: https://doi.org/10.26008/1912/bco-dmo.896862.2, Lange et al., 2024) and ancillary data.

Description: The field of oceanography is transitioning from data-poor to data-rich, thanks in part to increased deployment of in-situ platforms and sensors, such as those that instrument the US-funded Ocean Observatories Initiative (OOI). However, generating science-ready data products from these sensors, particularly those making biogeochemical measurements, often requires extensive end-user calibration and validation procedures, which can present a significant barrier. Openly available community-developed and -vetted Best Practices contribute to overcoming such barriers, but collaboratively developing user-friendly Best Practices can be challenging. Here we describe the process undertaken by the NSF-funded OOI Biogeochemical Sensor Data Working Group to develop Best Practices for creating science-ready biogeochemical data products from OOI data, culminating in the publication of the GOOS-endorsed OOI Biogeochemical Sensor Data Best Practices and User Guide. For Best Practices related to ocean observatories, engaging observatory staff is crucial, but having a “user-defined” process ensures the final product addresses user needs. Our process prioritized bringing together a diverse team and creating an inclusive environment where all participants could effectively contribute. Incorporating the perspectives of a wide range of experts and prospective end users through an iterative review process that included “Beta Testers’’ enabled us to produce a final product that combines technical information with a user-friendly structure that illustrates data analysis pipelines via flowcharts and worked examples accompanied by pseudo-code. Our process and its impact on improving the accessibility and utility of the end product provides a roadmap for other groups undertaking similar community-driven activities to develop and disseminate new Ocean Best Practices.

Description: We compare Holocene tree-cover changes in Europe derived from a transient MPI-ESM1.2 simulation with high spatial resolution LPJ-GUESS time-slice simulations and pollen-based quantitative reconstructions of tree cover based on the REVEALS model. The dynamic vegetation models and REVEALS agree with respect to the general temporal trends in tree cover for most parts of Europe, with a large tree cover during the mid-Holocene and a substantially smaller tree cover closer to the present time. However, the decrease in tree cover in REVEALS starts much earlier than in the models indicating much earlier anthropogenic deforestation than the prescribed land-use in the models. While LPJ-GUESS generally overestimates tree cover compared to the reconstructions, MPI-ESM indicates lower percentages of tree cover than REVEALS, particularly in Central Europe and the British Isles. A comparison of the simulated climate with chironomid-based climate reconstructions reveals that model-data mismatches in tree cover are in most cases not driven by biases in the climate. Instead, sensitivity experiments indicate that the model results strongly depend on the tuning of the models regarding natural disturbance regimes (e.g. fire and wind throw). The frequency and strength of disturbances are – like most of the parameters in the vegetation models – static and calibrated to modern conditions. However, these parameter values may not be valid during climate and vegetation states totally different from todays. In particular, the mid-Holocene natural forests were probably more stable and less sensitive to disturbances than present day forests that are heavily altered by human interventions. Our analysis highlights the fact that such model settings are inappropriate for palaeo-simulations and complicate model-data comparisons with additional challenges. Moreover, our study suggests that land-use is the main driver of forest decline in Europe during the mid- and late-Holocene.

Description: Nitrous oxide (N2O) is a long-lived potent greenhouse gas and stratospheric ozone-depleting substance, which has been accumulating in the atmosphere since the pre-industrial period. The mole fraction of atmospheric N2O has increased by nearly 25 % from 270 parts per billion (ppb) in 1750 to 336 ppb in 2022, with the fastest annual growth rate since 1980 of more than 1.3 ppb yr-1 in both 2020 and 2021. As a core component of our global greenhouse gas assessments coordinated by the Global Carbon Project (GCP), we present a global N2O budget that incorporates both natural and anthropogenic sources and sinks, and accounts for the interactions between nitrogen additions and the biochemical processes that control N2O emissions. We use Bottom-Up (BU: inventory, statistical extrapolation of flux measurements, process-based land and ocean modelling) and Top-Down (TD: atmospheric measurement-based inversion) approaches. We provide a comprehensive quantification of global N2O sources and sinks in 21 natural and anthropogenic categories in 18 regions between 1980 and 2020. We estimate that total annual anthropogenic N2O emissions increased 40 % (or 1.9 Tg N yr-1) in the past four decades (1980–2020). Direct agricultural emissions in 2020, 3.9 Tg N yr−1 (best estimate) represent the large majority of anthropogenic emissions, followed by other direct anthropogenic sources (including ‘Fossil fuel and industry’, ‘Waste and wastewater’, and ‘Biomass burning’ (2.1 Tg N yr−1), and indirect anthropogenic sources (1.3 Tg N yr−1). For the year 2020, our best estimate of total BU emissions for natural and anthropogenic sources was 18.3 (lower-upper bounds: 10.5–27.0) Tg N yr-1, close to our TD estimate of 17.0 (16.6–17.4) Tg N yr-1. For the period 2010–2019, the annual BU decadal-average emissions for natural plus anthropogenic sources were 18.1 (10.4–25.9) Tg N yr-1 and TD emissions were 17.4 (15.8–19.20 Tg N yr-1. The once top emitter Europe has reduced its emissions since the 1980s by 31 % while those of emerging economies have grown, making China the top emitter since the 2010s. The observed atmospheric N2O concentrations in recent years have exceeded projected levels under all scenarios in the Coupled Model Intercomparison Project Phase 6 (CMIP6), underscoring the urgency to reduce anthropogenic N2O emissions. To evaluate mitigation efforts and contribute to the Global Stocktake of the United Nations Framework Convention on Climate Change, we propose establishing a global network for monitoring and modeling N2O from the surface through the stratosphere. The data presented in this work can be downloaded from https://doi.org/10.18160/RQ8P-2Z4R (Tian et al. 2023).

Description: Black shale sediments from the Barremian to Aptian South Atlantic document intense and widespread burial of marine organic carbon during the initial stages of seafloor spreading between Africa and South America. The enhanced sequestration of atmospheric CO2 makes these young ocean basins potential drivers of the Early Cretaceous carbon cycle and climate perturbations. The opening of marine gateways between initially restricted basins and related circulation and ventilation changes are a commonly invoked explanation for the transient formation and disappearance of these regional carbon sinks. However, large uncertainties in paleogeographic reconstructions limit the interpretation of available paleoceanographic data and prevent any robust model-based quantifications of the proposed circulation and carbon burial changes. Here, we present a new approach to assess the principal controls on the Early Cretaceous South Atlantic and Southern Ocean circulation changes under full consideration of the uncertainties in available boundary conditions. Specifically, we use a large ensemble of 36 climate model experiments to simulate the Barremian to Albian progressive opening of the Falkland Plateau and Georgia Basin gateways with different configurations of the proto-Drake Passage, the Walvis Ridge, and atmospheric CO2 concentrations. The experiments are designed to complement available geochemical data across the regions and to test circulation scenarios derived from them. All simulations show increased evaporation and intermediate water formation at subtropical latitudes that drive a meridional overturning circulation whose vertical extent is determined by the sill depth of the Falkland Plateau. Densest water masses formed in the southern Angola Basin and potentially reached the deep Cape Basin as Walvis Ridge Overflow Water. Paleogeographic uncertainties are as important as the lack of precise knowledge of atmospheric CO2 levels for the simulated temperature and salinity spread in large parts of the South Atlantic. Overall temperature uncertainties are up to 15 °C and increase significantly with water depth. The ensemble approach reveals temporal changes in the relative importance of geographic and radiative forcings for the simulated oceanographic conditions and, importantly, nonlinear interactions between them. Progressive northward opening of the highly restricted Angola Basin increased the sensitivity of local overturning and upper ocean stratification to atmospheric CO2 concentrations due to large-scale changes in the hydrological cycle, while the chosen proto-Drake Passage depth is critical for the ocean dynamics and CO2 response in the southern South Atlantic. Finally, the simulated processes are integrated into a recent carbon burial framework to document the principal control of the regional gateway evolution on the progressive shift from the prevailing saline and oxygen-depleted subtropical water masses to the dominance of ventilated high-latitude deep waters.

Description: Ocean alkalinity enhancement (OAE) has been proposed as a carbon dioxide removal technology (CDR) allowing for long term storage of carbon dioxide in the ocean. By changing the carbonate speciation in seawater, OAE may potentially alter marine ecosystems with implications for the biological carbon pump. Using mesocosmsthe subtropical North Atlantic, we provide first empirical insights into impacts of carbonate-based OAE on the vertical flux and attenuation of sinking particles in an oligotrophic plankton community. We enhanced total alkalinity (TA) in increments of 300 μmol kg-1, reaching up to ΔTA = 2400 µmol kg-1 compared to ambient TA. We applied a pCO2-equilibrated OAE approach, i.e. dissolved inorganic carbon (DIC) was raised simultaneously with TA to maintain seawater pCO2 in equilibrium with the atmosphere, thereby keeping perturbations of seawater carbonate chemistry moderate. The vertical flux of major elements including carbon, nitrogen, phosphorus and silicon, as well as their stoichiometric ratios (e.g. carbon-to-nitrogen) remained unaffected over 29 days of OAE. The particle properties controlling the flux attenuationinking velocities and remineralization rates also remained unaffected by OAE. However, we observed abiotic mineral precipitation at high OAE levels (ΔTA = 1800 μmol kg-1 and higher) that resulted in a substantial increase in PIC formation. The associated consumption of alkalinity reduces the efficiency of CO2 removal and emphasizes the importance of maintaining OAE within a carefully defined operating range. Our findings suggest that carbon export by oligotrophic plankton communities is insensitive to OAE perturbations using a CO2 pre-equilibrated approach. The integrity of ecosystem services is a prerequisite for large-scale application and should be further tested across a variety of nutrient-regimes and for less idealized OAE approaches.

Description: An essential prerequisite for the implementation of ocean alkalinity enhancement (OAE) applications is their environmental safety. Only if it can be ensured that ecosystem health and ecosystem services are not at risk will the implementation of OAE move forward. Public opinion on OAEs will depend first and foremost on reliable evidence that no harm will be done to marine ecosystems and licensing authorities will demand measurable criteria against which environmental sustainability can be determined. In this context mesocosm experiments represent a highly valuable tool in determining the safe operating space of OAE applications. By combining realism and biological complexity with controllability and replication they provide an ideal OAE test bed and a critical stepping stone towards field applications. Mesocosm approaches can also be helpful in testing the efficacy, efficiency and permanence of OAE applications. This chapter outlines strengths and weaknesses of mesocosm approaches, illustrates mesocosm facilities and suitable experimental designs presently employed in OAE research, describes critical steps in mesocosm operation, and discusses possible approaches for alkalinity manipulation and monitoring. Building on a general treatise on each of these aspects, the chapter describes pelagic and benthic mesocosm approaches separately, given their inherent differences. The chapter concludes with recommendations for best practices in OAE-related mesocosm research.

Description: Ocean alkalinity enhancement (OAE) is an emerging strategy that aims to mitigate climate change by increasing the alkalinity of seawater. This approach involves increasing the alkalinity of the ocean to enhance its capacity to absorb and store carbon dioxide (CO2) from the atmosphere. This chapter presents an overview of the technical aspects associated with the full range of OAE methods being pursued and discusses implications for undertaking research on these approaches. Various methods have been developed to implement OAE, including the direct injection of alkaline liquid into the surface ocean; dispersal of alkaline particles from ships, platforms, or pipes; the addition of minerals to coastal environments; and the electrochemical removal of acid from seawater. Each method has its advantages and challenges, such as scalability, cost effectiveness, and potential environmental impacts. The choice of technique may depend on factors such as regional oceanographic conditions, alkalinity source availability, and engineering feasibility. This chapter considers electrochemical methods, the accelerated weathering of limestone, ocean liming, the creation of hydrated carbonates, and the addition of minerals to coastal environments. In each case, the technical aspects of the technologies are considered, and implications for best-practice research are drawn. The environmental and social impacts of OAE will likely depend on the specific technology and the local context in which it is deployed. Therefore, it is essential that the technical feasibility of OAE is undertaken in parallel with, and informed by, wider impact assessments. While OAE shows promise as a potential climate change mitigation strategy, it is essential to acknowledge its limitations and uncertainties. Further research and development are needed to understand the long-term effects, optimize techniques, and address potential unintended consequences. OAE should be viewed as complementary to extensive emission reductions, and its feasibility may be improved if it is operated using energy and supply chains with minimal CO2 emissions.

Description: The deliberate increase in ocean alkalinity (referred to as ocean alkalinity enhancement, or OAE) has been proposed as a method for removing CO2 from the atmosphere. Before OAE can be implemented safely, efficiently, and at scale several research questions have to be addressed, including (1) which alkaline feedstocks are best suited and the doses in which they can be added safely, (2) how net carbon uptake can be measured and verified, and (3) what the potential ecosystem impacts are. These research questions cannot be addressed by direct observation alone but will require skilful and fit-for-purpose models. This article provides an overview of the most relevant modelling tools, including turbulence-, regional-, and global-scale biogeochemical models and techniques including approaches for model validation, data assimilation, and uncertainty estimation. Typical bio- geochemical model assumptions and their limitations are discussed in the context of OAE research, which leads to an identification of further development needs to make models more applicable to OAE research questions. A description of typical steps in model validation is followed by proposed minimum criteria for what constitutes a model that is fit for its intended purpose. After providing an overview of approaches for sound integration of models and observations via data assimilation, the application of observing system simulation experiments (OSSEs) for observing system design is described within the context of OAE research. Criteria for model val- idation and intercomparison studies are presented. The article concludes with a summary of recommendations and potential pitfalls to be avoided.

Description: The Paris Agreement to limit global warming to well below 2 °C requires ambitious emission reduction and the balancing of remaining emissions through carbon sinks, i.e. the deployment of carbon dioxide removal (CDR). While ambitious climate protection scenarios until now consider primarily land-based CDR methods, there is growing concern about their potential to deliver sufficient CDR, and marine CDR options receive more and more interest. Based on idealized theoretical studies, Ocean Alkalinity Enhancement (OAE) appears as a promising marine CDR method. However, the knowledge base is insufficient for a robust assessment of its practical feasibility, of its side effects, social and governance aspects as well as monitoring, reporting and verification issues. A number of research efforts aim to improve this in a timely manner. We provide an overview on the current situation of developing OAE as marine CDR method, and describe the history that has led to the creation of the OAE research Best Practices Guide.

Description: Offshore freshened groundwater (OFG) is groundwater with a salinity below that of seawater that is stored in sub-seafloor sediments and rocks. OFG has been proposed as an alternative solution to relieve water scarcity in coastal regions and to enhance oil recovery. Although OFG has been documented in most continental margins, we still have a poor understanding of the extent and flow characteristics of OFG systems, and their evolution through time. In view of the general absence of appropriate field data, paleohydrogeological models have been used. The majority of these models are based on 2D approaches, and they rarely consider the future evolution of OFG systems, especially in response to predicted climate change. Here we utilize recently acquired geological, geophysical and hydrogeological data from onshore and offshore the Maltese Islands, and employ 2D and 3D numerical models, to: (i) reconstruct the evolution of the onshore-offshore groundwater system during the last 188 ka, (ii) predict the evolution of the OFG system in response to climate-related changes. We show that the mechanisms emplacing OFG include a combination of active meteoric recharge at present as well as at sea-level lowstands. The Maltese onshore-offshore groundwater system is relatively dynamic, with 23% of groundwater being preserved in the last 18 ka. The control of geology is expressed by the more prevalent distribution of OFG north of the Great Fault, which is associated to the occurrence of low permeability units, and the asymmetry of the groundwater lens during the 18 ka lowstand. A 30% decrease in recharge predicted in the coming 100 years will diminish OFG extent by 38%, whereas sea-level rise will play a negligible role. At present the estimated volume of OFG is 1 km 3 , which could potentially provide an alternative supply of potable water to the Maltese Islands for 75 years. Exploitation of OFG with minimal salinization of onshore groundwater bodies would require locating pumping wells close to the coast.

Description: Ocean alkalinity enhancement (OAE) is a proposed marine carbon dioxide removal (mCDR) approach that has the potential for large-scale uptake of significant amounts of atmospheric carbon dioxide (CO2). Removing anthropogenic legacy CO2 will be required to stabilise global surface temperatures below the 1.5–2 ∘C Paris Agreement target of 2015. In this chapter we describe the impacts of various OAE feedstocks on seawater carbonate chemistry, as well as pitfalls that need to be avoided during sampling, storage, and measurement of the four main carbonate chemistry parameters, i.e. dissolved inorganic carbon (DIC), total alkalinity (TA), pH, and CO2 fugacity (fCO2). Finally, we also discuss considerations in regard to calculating carbonate chemistry speciation from two measured parameters. Key findings are that (1) theoretical CO2 uptake potential (global mean of 0.84 mol of CO2 per mole of TA added) based on carbonate chemistry calculations is probably secondary in determining the oceanic region in which OAE would be best; (2) carbonate chemistry sampling is recommended to involve gentle pressure filtration to remove calcium carbonate (CaCO3) that might have been precipitated upon TA increase as it would otherwise interfere with a number of analyses; (3) samples for DIC and TA can be stabilised to avoid the risk of secondary CaCO3 precipitation during sample storage; and (4) some OAE feedstocks require additional adjustments to carbonate chemistry speciation calculations using available programs and routines, for instance if seawater magnesium or calcium concentrations are modified.

Description: Nitric oxide (NO) is an intermediate of various microbial nitrogen cycle processes and the open ocean and coastal areas are generally a source of NO in the atmosphere. However, our knowledge about its distribution and the main production processes in coastal areas and estuaries is rudimentary at best. To this end, dissolved NO concentrations were measured for the first time in surface waters along the lower Elbe Estuary and Hamburg Port area in July 2021. The discrete surface water samples were analyzed using a chemiluminescence detection method. The NO concentrations ranged from below the limit of detection (9.1 × 10−12 mol L−1) to 17.7 × 10−12 mol L−1, averaging at 12.5 × 10−12 mol L−1 and were supersaturated in the surface layer of both the lower Elbe Estuary and the Hamburg Port area, indicating that the study site was a source of NO to the atmosphere during the study period. On the basis of a comprehensive comparison of NO concentrations with parallel nutrient, oxygen, and nitrous oxide concentration measurements, we conclude that the observed distribution of dissolved NO was most likely resulting from microbial nitrogen transformation processes, particularly nitrification in the coastal-brackish and limnic zones of the lower Elbe Estuary and nitrifier-denitrification in the Hamburg Port area.

Description: Carbon monoxide (CO) is an atmospheric trace gas that plays a crucial role in the oxidizing capacity of the Earth’s atmosphere. Moreover, it functions as an indirect greenhouse gas, influencing the lifetimes of potent greenhouse gases such as methane. Albeit being an overall source of atmospheric CO, the role of coastal regions in the marine cycling of CO and how its budget can be affected by anthropogenic activities, remain uncertain. Here, we present the first measurements of dissolved CO in the Ria Formosa Lagoon, an anthropogenically influenced system in southern Portugal. The dissolved CO concentrations in the surface layer ranged from 0.16 to 3.1 nmol L−1 with an average concentration of 0.75 ± 0.57 nmol L−1. The CO saturation ratio ranged from 1.7 to 32.2, indicating that the lagoon acted as a source of CO to the atmosphere in May 2021. The estimated average sea-to-air flux density was 1.53 μmol m−2 d−1, mainly fueled by CO photochemical production. Microbial consumption accounted for 83 % of the CO production, suggesting that the resulting CO emissions to the atmosphere were modulated by microbial consumption in the surface waters of the Ria Formosa Lagoon. The results from an irradiation experiment with aquaculture effluent water indicated that aquaculture facilities in the Ria Formosa Lagoon seem to be a negligible source of atmospheric CO.

Description: Ocean alkalinity enhancement (OAE) is a marine carbon dioxide removal (CDR) approach. Publicly funded research projects have begun, and philanthropic funding and start-ups are collectively pushing the field forward. This rapid progress in research activities has created an urgent need to learn if and how OAE can work at scale. The Best Practices Guide to OAE research contains 7 topics broken down into 13 chapters that compare and synthesise previously published methods and offer guidance for future research.

Description: The Humboldt Current Upwelling System (HCS) is the most productive eastern boundary upwelling system (EBUS) in terms of fishery yield on the planet. EBUSs are considered hotspots of climate change with predicted expansion of mesopelagic oxygen minimum zones (OMZs) and related changes in the frequency and intensity of upwelling of nutrient-rich, low-oxygen deep water. To increase our mechanistic understanding of how upwelling impacts plankton communities and trophic links, we investigated mesozooplankton community succession and gut fluorescence, fatty acid and elemental compositions (C, N, O, P), and stable isotope (δ13C, δ15N) ratios of dominant mesozooplankton and microzooplankton representatives in a mesocosm setup off Callao (Peru) after simulated upwelling with OMZ water from two different locations and different N:P signatures (moderate and extreme treatments). An oxycline between 5 and 15 m with hypoxic conditions (〈50 µmol L−1) below ∼10 m persisted in the mesocosms throughout the experiment. No treatment effects were determined for the measured parameters, but differences in nutrient concentrations established through OMZ water additions were only minor. Copepods and polychaete larvae dominated in terms of abundance and biomass. Development and reproduction of the dominant copepod genera Paracalanus sp., Hemicyclops sp., Acartia sp., and Oncaea sp. were hindered as evident from accumulation of adult copepodids but largely missing nauplii. Failed hatching of nauplii in the hypoxic bottom layer of the mesocosms and poor nutritional condition of copepods suggested from very low gut fluorescence and fatty acid compositions most likely explain the retarded copepod development. Correlation analysis revealed no particular trophic relations between dominant copepods and phytoplankton groups. Possibly, particulate organic matter with a relatively high C:N ratio was a major diet of copepods. C:N ratios of copepods and polychaetes ranged 4.8–5.8 and 4.2–4.3, respectively. δ15N was comparatively high (∼13 ‰–17 ‰), potentially because the injected OMZ source water was enriched in δ15N as a result of anoxic conditions. Elemental ratios of dinoflagellates deviated strongly from the Redfield ratio. We conclude that opportunistic feeding of copepods may have played an important role in the pelagic food web. Overall, projected changes in the frequency and intensity of upwelling hypoxic waters may make a huge difference for copepod reproduction and may be further enhanced by varying N:P ratios of upwelled OMZ water masses.

Description: When interpreting geophysical models, we need to establish a link between the models’ physical parameters and geological units. To define these connections, it is crucial to consider and compare geophysical models with multiple, independent parameters. Particularly in complex geological scenarios, such as the rifted passive margin offshore Namibia, multi-parameter analysis and joint inversion are key techniques for comprehensive geological inferences. The models resulting from joint inversion enable the definition of specific parameter combinations, which can then be ascribed to geological units. Here we perform a user-unbiased clustering analysis of the parameters electrical resistivity and density from two models derived in a joint inversion along the Namibian passive margin. We link the resulting parameter combinations to break-up related lithology, and infer the history of margin formation. This analysis enables us to clearly differentiate two types of sediment cover. Namely, one of near-shore, thick, clastic sediments, and a second one of further offshore located, more biogenic, marine sediments. Furthermore, we clearly identify areas of interlayered massive, and weathered volcanic flows, which are usually only identified in reflection seismic studies as seaward dipping reflectors. Lastly, we find a distinct difference in the signature of the transitional crust south of- and along the supposed hot-spot track Walvis Ridge. We ascribe this contrast to an increase in magmatic activity above the volcanic centre along Walvis Ridge, and potentially a change in melt sources or depth of melting. This characterizes a rift-related southern complex, and a plume-driven Walvis Ridge regime. All of these observations demonstrate the importance of multi-parameter geophysical analysis for large-scale geological interpretations. Furthermore, our results may improve future joint inversions using direct parameter coupling, by providing a guideline for the complex passive margins parameter correlations.

Description: In this paper, we review observational and modelling results on the upwelling in the tropical Atlantic between 10∘ N and 20∘ S. We focus on the physical processes that drive the seasonal variability of surface cooling and the upward nutrient flux required to explain the seasonality of biological productivity. We separately consider the equatorial upwelling system, the coastal upwelling system of the Gulf of Guinea and the tropical Angolan upwelling system. All three tropical Atlantic upwelling systems have in common a strong seasonal cycle, with peak biological productivity during boreal summer. However, the physical processes driving the upwelling vary between the three systems. For the equatorial regime, we discuss the wind forcing of upwelling velocity and turbulent mixing, as well as the underlying dynamics responsible for thermocline movements and current structure. The coastal upwelling system in the Gulf of Guinea is located along its northern boundary and is driven by both local and remote forcing. Particular emphasis is placed on the Guinea Current, its separation from the coast and the shape of the coastline. For the tropical Angolan upwelling, we show that this system is not driven by local winds but instead results from the combined effect of coastally trapped waves, surface heat and freshwater fluxes, and turbulent mixing. Finally, we review recent changes in the upwelling systems associated with climate variability and global warming and address possible responses of upwelling systems in future scenarios.

Description: We use output from a freely-running NEMO model simulation for the equatorial Pacific to investigate the utility of linearly removing the local influence of vertical displacements of the thermocline from variations in sea surface height. We show that the resulting time series of residual sea surface height, denoted ηnlti, measures variations in near-surface heat content that are independent of the local vertical displacement of the thermocline and can arise from horizontal advection, surface heat flux and diapycnal mixing processes. We find that the variance of ηnlti and its correlation with sea surface temperature, are focused on the Niño4 region. Furthermore, ηnlti averaged over the Niño4 region is highly correlated with indices of Central Pacific El Niño Southern Oscillation (CP ENSO), and its variance in 21 year running windows shows a strong upward trend over the past 50 years, corresponding to the emergence of CP ENSO following the 1976/77 climate shift. We show that ηnlti can be estimated from observations, using satellite altimeter data and a linear multi-mode model. The time series of ηnlti, especially when estimated using the linear model, show pronounced westward propagation in the western equatorial Pacific, arguing an important role for zonal advective feedback in the dynamics of CP ENSO, in particular for cold events. We also present evidence that the role of the thermocline displacement in influencing sea surface height increased strongly after 2000 in the eastern part of the Niño4 region, at a time when CP ENSO was particularly active. Finally, the diagnostic is easy to compute and can be easily applied to mooring data or couple climate models.

Description: Barium (Ba) isotopes are a promising new tracer for riverine freshwater input to the ocean and marine biogeochemical cycling. However, many processes that affect Ba cycling at continental margins have not yet been investigated with respect to Ba isotope fractionation. Here, we present a comprehensive data set of Ba concentration and isotope data for water column, pore water and sediment samples from Kiel Bight, a seasonally stratified and hypoxic fjord in the southwestern Baltic Sea. The surface water Ba concentration and Ba isotope inventory of the water column can generally be explained by mixing of riverine freshwater and Atlantic seawater. However, the deep-water below the seasonal pycnocline (10 - 15 m water depth) is characterized by a pronounced positive Ba concentration anomaly (up to 915 nM) that is accompanied by a δ138Ba of ~+0.25 ‰, which is lighter than expected from the seawater-freshwater mixing line (Ba: 77 nM, δ138Ba: +0.32 ‰ at a salinity of 18). Pore water profiles indicate a Ba flux across the sediment-water interface, which contributes to the enrichment in isotopically light Ba in the deep-water. Pore waters of surface sediments and deep-waters are oversaturated with respect to barite. Therefore, barite dissolution is unlikely to account for the benthic Ba flux. Water column Ba concentrations closely correlate with those of the nutrients phosphate and silica, which are removed from surface waters by biological processes and recycled from the sediment by diffusion across the sediment-water interface. As nutrient-to-Ba ratios differ among sites and from those observed in open-marine systems, we propose that Ba is removed from surface waters by adsorption onto biogenic particles (rather than assimilation) and regenerated within surface sediments upon organic matter degradation. Pore water data for subsurface sediments in Kiel Bight indicate preferential transfer of isotopically heavy Ba into an authigenic phase during early diagenesis. Quantifying the burial flux associated with this authigenic Ba phase along continental margins could potentially help to settle the isotopic imbalance between known Ba source and sink fluxes in the ocean.

Description: Nitrogen (N) and phosphorus (P) biogeochemical dynamics are crucial for the regulation of the terrestrial carbon cycle. In Earth system models (ESMs) the implementation of nutrient limitations has been shown to improve the carbon cycle feedback representation and, hence, the fidelity of the response of land to simulated atmospheric CO2 rise. Here we aimed to implement a terrestrial N and P cycle in an Earth system model of intermediate complexity to improve projections of future CO2 fertilization feedbacks. The N cycle is an improved version of the Wania et al. (2012) N module, with enforcement of N mass conservation and the merger with a deep land-surface and wetland module that allows for the estimation of N2O and NO fluxes. The N cycle module estimates fluxes from three organic (litter, soil organic matter and vegetation) and two inorganic ( and ) pools and accounts for inputs from biological N fixation and N deposition. The P cycle module contains the same organic pools with one inorganic P pool; it estimates influx of P from rock weathering and losses from leaching and occlusion. Two historical simulations are carried out for the different nutrient limitation setups of the model: carbon and nitrogen (CN), as well as carbon, nitrogen and phosphorus (CNP), with a baseline carbon-only simulation. The improved N cycle module now conserves mass, and the added fluxes (NO and N2O), along with the N and P pools, are within the range of other studies and literature. For the years 2001–2015 the nutrient limitation resulted in a reduction of gross primary productivity (GPP) from the carbon-only value of 143 to 130 Pg C yr−1 in the CN version and 127 Pg C yr−1 in the CNP version. This implies that the model efficiently represents a nutrient limitation over the CO2 fertilization effect. CNP simulation resulted in a reduction of 11 % of the mean GPP and a reduction of 23 % of the vegetation biomass compared to the baseline C simulation. These results are in better agreement with observations, particularly in tropical regions where P limitation is known to be important. In summary, the implementation of the N and P cycle has successfully enforced a nutrient limitation in the terrestrial system, which has now reduced the primary productivity and the capacity of land to take up atmospheric carbon, better matching observations.

Description: State-of-the-art Earth system models typically employ grid spacings of O(100 km), which is too coarse to explicitly resolve main drivers of the flow of energy and matter across the Earth system. In this paper, we present the new ICON-Sapphire model configuration, which targets a representation of the components of the Earth system and their interactions with a grid spacing of 10 km and finer. Through the use of selected simulation examples, we demonstrate that ICON-Sapphire can (i) be run coupled globally on seasonal timescales with a grid spacing of 5 km, on monthly timescales with a grid spacing of 2.5 km, and on daily timescales with a grid spacing of 1.25 km; (ii) resolve large eddies in the atmosphere using hectometer grid spacings on limited-area domains in atmosphere-only simulations; (iii) resolve submesoscale ocean eddies by using a global uniform grid of 1.25 km or a telescoping grid with the finest grid spacing at 530 m, the latter coupled to a uniform atmosphere; and (iv) simulate biogeochemistry in an ocean-only simulation integrated for 4 years at 10 km. Comparison of basic features of the climate system to observations reveals no obvious pitfalls, even though some observed aspects remain difficult to capture. The throughput of the coupled 5 km global simulation is 126 simulated days per day employing 21 % of the latest machine of the German Climate Computing Center. Extrapolating from these results, multi-decadal global simulations including interactive carbon are now possible, and short global simulations resolving large eddies in the atmosphere and submesoscale eddies in the ocean are within reach.

Description: Global biogeochemical ocean models help to investigate the present and potential future state of the ocean, its productivity and cascading effects on higher trophic levels such as fish. They are often subjectively tuned against data sets of inorganic tracers and surface chlorophyll and only very rarely against organic components such as particulate organic carbon or zooplankton. The resulting uncertainty in biogeochemical model parameters (and parameterisations) associated with these components can explain some of the large spread of global model solutions with regard to the cycling of organic matter and its impacts on biogeochemical tracer distributions, such as oxygen minimum zones (OMZs). A second source of uncertainty arises from differences in the model spin-up length as, so far, there seems to be no agreement on the required simulation time that should elapse before a global model is assessed against observations. We investigated these two sources of uncertainty by optimising a global biogeochemical ocean model against the root-mean-squared error (RMSE) of six different combinations of data sets and different spin-up times. Besides nutrients and oxygen, the observational data sets also included phyto- and zooplankton, as well as dissolved and particulate organic phosphorus (DOP and POP, respectively). We further analysed the optimised model performance with regard to global biogeochemical fluxes, oxygen inventory and OMZ volume. Following the optimisation procedure, we evaluated the RMSE for all tracers located in the upper 100 m (except for POP, for which we considered the entire vertical domain), regardless of their consideration during optimisation. For the different optimal model solutions, we find a narrow range of the RMSE, between 14 % of the average RMSE after 10 years and 24 % after 3000 years of simulation. Global biogeochemical fluxes, global oxygen bias and OMZ volume showed a much stronger divergence among the models and over time than RMSE, indicating that even models that are similar with regard to local surface tracer concentrations can perform very differently when assessed against the global diagnostics for oxygen. Considering organic tracers in the optimisation had a strong impact on the particle flux exponent (Martin b) and may reduce much of the uncertainty in this parameter and the resulting deep particle flux. Independent of the optimisation setup, the OMZ volume showed a particularly sensitive response with strong trends over time, even after 3000 years of simulation time (despite the constant physical forcing); a high sensitivity to simulation time; and the highest sensitivity to model parameters arising from the tuning strategy setup (variation of almost 80 % of the ensemble mean). In conclusion, calibration against observations of organic tracers can help to improve global biogeochemical models even after short spin-up times; here especially, observations of deep particle flux could provide a powerful constraint. However, a large uncertainty remains with regard to global OMZ volume and its evolution over time, which can show very dynamic behaviour during the model spin-up, which renders temporal extrapolation to a final equilibrium state difficult if not impossible. Given that the real ocean shows variations on many timescales, the assumption of observations representing a steady-state ocean may require some reconsideration.

Description: Dissolved silicate (H4SiO4) is essential for the formation of the opaline skeletal structures of diatoms and other siliceous plankton. A fraction of particulate biogenic silica (bSi) formed in surface waters sinks to the seabed, where it either dissolves and returns to the water column or is permanently buried. Global silica budgets are still poorly constrained since data on benthic bSi cycling are lacking, especially on continental margins. This study describes benthic bSi cycling in the Skagerrak, a sedimentary depocenter for particles from the North Sea. Biogenic silica burial fluxes, benthic H4SiO4 fluxes to the water column and bSi burial efficiencies are reported for nine stations by evaluating data from in-situ benthic landers and sediment cores with a diagenetic reaction-transport model. The model simulates bSi contents and H4SiO4 concentrations at all sites using a novel power law to describe bSi dissolution kinetics with a small number of adjustable parameters. Our results show that, on average, 1100 mmol m-2 yr-1 of bSi rains down to the Skagerrak basin seafloor, of which 50% is released back to overlying waters, with the remainder being buried. Biogenic silica cycling in the Skagerrak is generally consistent with previously reported global trends, showing higher Si fluxes and burial efficiencies than deep-sea sites and similar values compared to other continental margins. A significant finding of this work is a molar bSi-to-organic carbon burial ratio of 0.22 in Skagerrak sediments, which is distinctively lower compared to other continental margins. We suggest that the continuous dissolution of bSi in suspended sediments transported over long distances from the North Sea leads to the apparent decoupling between bSi and organic carbon in Skagerrak sediments.

Description: Volcanic flanks subject to hydrothermal alteration become mechanically weak and gravitationally unstable, which may collapse and develop far-reaching landslides. The dynamics and trajectories of volcanic landslides are hardly preserved and challenging to determine, which is due to the steep slopes and the inherent instability. Here we analyze the proximal deposits of the 21 July 2014, landslide at Askja (Iceland), by combining high-resolution imagery from satellites and Unoccupied Aircraft Systems. We performed a Principal Component Analysis in combination with supervised classification to identify different material classes and altered rocks. We trained a maximum-likelihood classifier and were able to distinguish 7 different material classes and compare these to ground-based hyperspectral measurements that we conducted on different rock types found in the field. Results underline that the Northern part of the landslide source region is a hydrothermally altered material class, which bifurcates halfway downslope and then extends to the lake. We find that a large portion of this material is originating from a lava body at the landslide headwall, which is the persistent site of intense hydrothermal activity. By comparing the classification result to in-situ hyperspectral measurements, we were able to further identify the involved types of rocks and the degree of hydrothermal alteration. We further discuss associated effects of mechanical weakening and the relevance of the heterogeneous materials for the dynamics and processes of the landslide. As the study demonstrates the success of our approach for identification of altered and less altered materials, important implications for hazard assessment in the Askja caldera and elsewhere can be drawn.

Description: The Angolan shelf system represents a highly productive ecosystem. Throughout the year sea surface temperatures (SSTs) are cooler near the coast than further offshore. Lowest SSTs, the strongest cross-shore temperature gradient and maximum productivity occur in austral winter when seasonally prevailing upwelling favourable winds are weakest. Here, we investigate the seasonal mixed layer heat budget to analyse atmospheric and oceanic causes for heat content variability. By using different satellite and in-situ data, we derive monthly estimates of surface heat fluxes, mean horizontal advection and local heat content change. We calculate the heat budgets for the near coastal and offshore regions separately to explore processes that lead to the observed differences. The results show that the net surface heat flux warms the coastal ocean stronger than further offshore thus acting to damp spatial SST differences. Mean horizontal heat advection is dominated by meridional advection of warm water along the Angolan coast. However, its contribution to the heat budget is small. Ocean turbulence data suggests that the heat flux due to turbulent mixing across the base of the mixed layer is an important cooling term. This turbulent cooling that is strongest in shallow shelf regions is capable of explaining the observed negative cross-shore temperature gradient. The residuum of the mixed layer heat budget and uncertainties of budget terms are discussed.

Description: Carbon monoxide (CO) influences the radiative budget and oxidative capacity of the atmosphere over the Arctic Ocean, which is a source of atmospheric CO. Yet, oceanic CO cycling is understudied in this area, particu- larly in light of the ongoing rapid environmental changes. We present results from incubation experiments conducted in the Fram Strait in August–September 2019 under different environmental conditions: while lower pH did not affect CO production (GPCO) or consumption (kCO) rates, enhanced GPCO and kCO were positively correlated with coloured dis- solved organic matter (CDOM) and dissolved nitrate concen- trations, respectively, suggesting microbial CO uptake under oligotrophic conditions to be a driving factor for variability in CO surface concentrations. Both production and consump- tion of CO will likely increase in the future, but it is un- known which process will dominate. Our results will help to improve models predicting future CO concentrations and emissions and their effects on the radiative budget and the oxidative capacity of the Arctic atmosphere.

Description: The carbon cycle component of the newly developed Earth System Model of intermediate complexity CLIMBER-X is presented. The model represents the cycling of carbon through atmosphere, vegetation, soils, seawater and marine sediments. Exchanges of carbon with geological reservoirs occur through sediment burial, rock weathering and volcanic degassing. The state-of-the-art HAMOCC6 model is employed to simulate ocean biogeochemistry and marine sediments processes. The land model PALADYN simulates the processes related to vegetation and soil carbon dynamics, including permafrost and peatlands. The dust cycle in the model allows for an interactive determination of the input of the micro-nutrient iron into the ocean. A rock weathering scheme is implemented into the model, with the weathering rate depending on lithology, runoff and soil temperature. CLIMBER-X includes a simple representation of the methane cycle, with explicitly modelled natural emissions from land and the assumption of a constant residence time of CH4 in the atmosphere. Carbon isotopes 13C and 14C are tracked through all model compartments and provide a useful diagnostic for model-data comparison. A comprehensive evaluation of the model performance for present–day and the historical period shows that CLIMBER-X is capable of realistically reproducing the historical evolution of atmospheric CO2 and CH4, but also the spatial distribution of carbon on land and the 3D structure of biogeochemical ocean tracers. The analysis of model performance is complemented by an assessment of carbon cycle feedbacks and model sensitivities compared to state-of-the-art CMIP6 models. Enabling interactive carbon cycle in CLIMBER-X results in a relatively minor slow-down of model computational performance by ~20 %, compared to a throughput of ~10,000 simulation years per day on a single node with 16 CPUs on a high performance computer in a climate–only model setup. CLIMBER-X is therefore well suited to investigate the feedbacks between climate and the carbon cycle on temporal scales ranging from decades to 〉100,000 years.

Description: Introduction: The associated diverse microbiome contributes to the overall fitness of Aurelia aurita, particularly to asexual reproduction. However, how A. aurita maintains this specific microbiome or reacts to manipulations is unknown. Methods: In this report, the response of A. aurita to manipulations of its native microbiome was studied by a transcriptomics approach. Microbiome-manipulated polyps were generated by antibiotic treatment and challenging polyps with a non-native, native, and potentially pathogenic bacterium. Total RNA extraction followed by RNAseq resulted in over 155 million reads used for a de novo assembly. Results: The transcriptome analysis showed that the antibiotic-induced change and resulting reduction of the microbiome significantly affected the host transcriptome, e.g., genes involved in processes related to immune response and defense mechanisms were highly upregulated. Similarly, manipulating the microbiome by challenging the polyp with a high load of bacteria (2 × 107 cells/polyp) resulted in induced transcription of apoptosis-, defense-, and immune response genes. A second focus was on host-derived quorum sensing interference as a potential defense strategy. Quorum Quenching (QQ) activities and the respective encoding QQ-ORFs of A. aurita were identified by functional screening a cDNA-based expression library generated in Escherichia coli. Corresponding sequences were identified in the transcriptome assembly. Moreover, gene expression analysis revealed differential expression of QQ genes depending on the treatment, strongly suggesting QQ as an additional defense strategy. Discussion: Overall, this study allows first insights into A. aurita’s response to manipulating its microbiome, thus paving the way for an in-depth analysis of the basal immune system and additional fundamental defense strategies.

Description: Communicating environmental change and mitigation scenarios to stakeholders and decision-makers can be challenging. Immersive environments offer an innovative approach for knowledge transfer, allowing science-based scenarios to be discussed interactively. The use of such environments is particularly helpful for the analysis of large, multi-component geospatial datasets, as commonly employed in the classification of ecosystems. Virtual environments can play an important role in conveying and discussing the findings gathered from these geomorphometric datasets. However, textured meshes and point clouds are not always well suited for direct import to a virtual reality or the creation of a truly immersive environment, and often result in geometrical artifacts, which can be misinterpreted during the import to a game engine. Such technical hurdles may lead to viewers rejecting the experience altogether, failing to achieve a higher educational purpose. In this study, we apply an asset-based approach to create an immersive virtual representation of a coastal environment. The focus hereby is on the coastal vegetation and changes in species distribution, which could potentially be triggered by the impact of climate change. We present an easy-to-use blueprint for the game engine EPIC Unreal Engine 5. In contrast to traditional virtual reality environments, which use static textured mesh data derived from photogrammetry, this asset-based approach enables the use of dynamic and physical properties (e.g. vegetation moving due to wind or waves), which makes the virtual environment more immersive. This will help to stimulate understanding and discussion amongst different stakeholders, and will also help to foster inclusion in earth- and environmental science education.

Description: Enhanced mineral dissolution in the benthic environment is currently discussed as a potential technique for ocean alkalinity enhancement (OAE) to reduce atmospheric CO 2 levels. This study explores how biogeochemical processes affect the dissolution of alkaline minerals in surface sediments during laboratory incubation experiments. These involved introducing dunite and calcite to organic-rich sediments from the Baltic Sea under controlled conditions in an oxic environment. The sediment cores were incubated with Baltic Sea bottom water. Findings reveal that the addition of calcite increased the benthic alkalinity release from 0.4 μmol cm −2 d −1 (control) to 1.4 μmol cm −2 d −1 (calcite) as well as other weathering products such as calcium. However, these enhanced fluxes returned to lower fluxes after approximately 4 weeks yet still higher than the un-amended controls. Microbial activity appeared to be the primary driver for lowering pore water pH and thus enhanced weathering. In several sediment cores, pH profiles taken at the start of the experiments indicated activity of sulfur oxidizing Beggiatoa spp, which was verified by RNA-profiling of 16S rRNA genes. The pH profiles transitioned to those commonly associated with the activity of cable bacteria as the experiments progressed. The metabolic activity of cable bacteria would explain the significantly lower pH values (~5.6) at sediment depths of 1–3 cm, which would favor substantial calcite dissolution. However, a high abundance of cable bacteria was not reflected in 16S rRNA sequence data. Total alkalinity (TA) fluxes in these cores increased by a factor of ~3, with excess TA/calcium ratios indicating that the enhanced flux originated from calcite dissolution. The dissolution of dunite or the potential formation of secondary minerals could not be identified due to the strong natural flux of silicic acid, likely due to biogenic silica dissolution. Furthermore, no accumulation of potentially harmful metals such as nickel was observed, as highlighted as a potential risk in other studies concerning OAE. Given the complexity of sediment chemistry and changes of the benthic conditions induced by the incubation, it remains challenging to distinguish between natural and enhanced mineral weathering. Further investigation, including the identification of suitable tracers for mineral dissolution, are necessary to assess the feasibility of benthic weathering as a practical approach for OAE and climate change mitigation.

Description: The Earth climate system is out of energy balance, and heat has accumulated continuously over the past decades, warming the ocean, the land, the cryosphere, and the atmosphere. According to the Sixth Assessment Report by Working Group I of the Intergovernmental Panel on Climate Change, this planetary warming over multiple decades is human-driven and results in unprecedented and committed changes to the Earth system, with adverse impacts for ecosystems and human systems. The Earth heat inventory provides a measure of the Earth energy imbalance (EEI) and allows for quantifying how much heat has accumulated in the Earth system, as well as where the heat is stored. Here we show that the Earth system has continued to accumulate heat, with 381±61 ZJ accumulated from 1971 to 2020. This is equivalent to a heating rate (i.e., the EEI) of 0.48±0.1 W m−2. The majority, about 89 %, of this heat is stored in the ocean, followed by about 6 % on land, 1 % in the atmosphere, and about 4 % available for melting the cryosphere. Over the most recent period (2006–2020), the EEI amounts to 0.76±0.2 W m−2. The Earth energy imbalance is the most fundamental global climate indicator that the scientific community and the public can use as the measure of how well the world is doing in the task of bringing anthropogenic climate change under control. Moreover, this indicator is highly complementary to other established ones like global mean surface temperature as it represents a robust measure of the rate of climate change and its future commitment. We call for an implementation of the Earth energy imbalance into the Paris Agreement's Global Stocktake based on best available science. The Earth heat inventory in this study, updated from von Schuckmann et al. (2020), is underpinned by worldwide multidisciplinary collaboration and demonstrates the critical importance of concerted international efforts for climate change monitoring and community-based recommendations and we also call for urgently needed actions for enabling continuity, archiving, rescuing, and calibrating efforts to assure improved and long-term monitoring capacity of the global climate observing system. The data for the Earth heat inventory are publicly available, and more details are provided in Table 4.

Description: The North Pacific Subtropical Countercurrent area (STCC) is high in mesoscale eddy activities. According to the rotation direction of the eddy flow field and the sign of temperature anomaly within the eddy, they can be divided into four categories: cyclonic cold-core eddy (CCE), anticyclonic warm-core eddy (AWE), cyclonic warm-core eddy (CWE) and anticyclonic cold-core eddy (ACE). CCE and AWE are called normal eddies, and CWE and ACE are named abnormal eddies. Based on the OFES data and vector geometry automatic detection method, we find that at the sea surface, the maximum monthly number of the CCE, AWE, CWE, and ACE occurs in December (765.70 ± 52.05), January (688.20 ± 82.53), August (373.40 ± 43.09) and August (533.00 ± 56.92), respectively. The number of normal eddies is more in winter and spring, and less in summer and autumn, while abnormal eddies have the opposite distribution. The maximum rotation velocity of the four types of eddies appears in June (11.71 ± 0.75 cm/s), June (12.24 ± 0.86 cm/s), May (10.63 ± 0.99 cm/s) and June (9.97 ± 0.91 cm/s), which is fast in winter and spring. The moving speed of the four types of eddies is almost similar (about 10 ~ 11 cm/s). The amplitude of normal and abnormal eddies is both high in summer and autumn, and low in winter and spring, with larger amplitudes in normal than abnormal eddies. The eccentricity (defined as the eccentricity of the ellipse obtained by fitting the eddy boundary) of the four types of eddies is also close to each other, and their variation ranges from 0.7 to 0.8, with no apparent seasonal variation. The vertical penetration depth, which has no significant seasonal difference, is 675.13 ± 67.50 m in cyclonic eddies (CCE and CWE), which is deeper than that 622.32 ± 81.85 m in anticyclonic eddies (ACE and AWE). In addition, increasing the defined temperature threshold for abnormal eddies can significantly reduce their numbers but does not change their seasonal variation trend.

Description: Microplastic is a ubiquitous marine pollutant whose small dimensions make it biologically available to phytoplankton and zooplankton. These organisms are crucial as the basis of the marine food web and for the export of organic material in the form of faecal pellets from the surface to deeper in the water column, forming a long-term carbon sink. Previous laboratory studies have demonstrated empirically that ingestion of low density microplastics reduces the sinking rates of zooplankton faecal pellets. This study uses a complex earth system model to analyse this effect and assess its wider impacts in a changing climate. Results show that the slowing of faecal pellet sinking stimulates changes to ecosystems regionally and reduces ocean carbon uptake by about 4.4 Pg C between the years 1950-2100, 0.24% of anthropogenic emissions over this time. However, perturbation of organic particle fluxes is significant, especially in gyres, and of the order of climate change impacts over the same time period. We calculate that plastics carbon has a 3 orders of magnitude greater impact on marine ecosystems than atmospheric carbon over our centennial timescale. Large uncertainties in model parameters and simplistic model structure suggest our results should be interpreted as motivation to further investigate parameter estimation, calcification responses to pollution, and the combined effects of multiple impact mechanisms on ecosystems.

Description: For millennia, humans have gravitated towards coastlines for their resource potential and as geopolitical centres for global trade. A basic requirement ensuring water security for coastal communities relies on a delicate balance between the supply and demand of potable water. The interaction between freshwater and saltwater in coastal settings is, therefore, complicated by both natural and human-driven environmental changes at the land-sea interface. In particular, ongoing sea level rise, warming and deoxygenation might exacerbate such perturbations. In this context, an improved understanding of the nature and variability of groundwater fluxes across the land-sea continuum is timely, yet remains out of reach. The flow of terrestrial groundwater across the coastal transition zone as well as the extent of freshened groundwater below the present-day seafloor are receiving increased attention in marine and coastal sciences because they likely represent a significant, yet highly uncertain component of (bio)geochemical budgets, and because of the emerging interest in the potential use of offshore freshened groundwater as a resource. At the same time, “reverse” groundwater flux from offshore to onshore is of prevalent socio-economic interest as terrestrial groundwater resources are continuously pressured by overpumping and seawater intrusion in many coastal regions worldwide. An accurate assessment of the land-ocean connectivity through groundwater and its potential responses to future anthropogenic activities and climate change will require a multidisciplinary approach combining the expertise of geophysicists, hydrogeologists, (bio)geochemists and modellers. Such joint activities will lay the scientific basis for better understanding the role of groundwater in societal-relevant issues such as climate change, pollution and the environmental status of the coastal oceans within the framework of the United Nations Sustainable Development Goals. Here, we present our perspectives on future research directions to better understand land-ocean connectivity through groundwater, including the spatial distributions of the essential hydrogeological parameters, highlighting technical and scientific developments, and briefly discussing its societal relevance in rapidly changing coastal oceans.

Description: Seagrass meadows have a disproportionally high organic carbon (Corg) storage potential within their sediments and thus can play a role in climate change mitigation via their conservation and restoration. However, high spatial heterogeneity is observed in Corg, with wide differences seen globally, regionally, and even locally (within a seagrass meadow). Consequently, it is difficult to determine their contributions to the national remaining carbon dioxide (CO2) budget without introducing a large degree of uncertainty. To address this spatial heterogeneity, we sampled 20 locations across the German Baltic Sea to quantify Corg stocks and sources in Zostera marina seagrass-vegetated and adjacent unvegetated sediments. To predict and integrate the Corg inventory in space, we measured the physical (seawater depth, sediment grain size, current velocity at the seafloor, anthropogenic inputs) and biological (seagrass complexity) environments to determine regional and local drivers of Corg variation. Here, we show that seagrass meadows in Germany constitute a significant Corg stock, storing on average 7,785 g C/m2, 13 times greater than meadows from other parts of the Baltic Sea, and fourfold richer than adjacent unvegetated sediments. Stocks were highly heterogenous; they differed widely between (by 10-fold) and even within (by 3- to 55-fold) sites. Regionally, Corg was controlled by seagrass complexity, fine sediment fraction, and seawater depth. Autochthonous material contributed to 78% of the total Corg in seagrass-vegetated sediments, and the remaining 22% originated from allochthonous sources (phytoplankton and macroalgae). However, relic terrestrial peatland material, deposited approximately 6,000 years BP during the last deglaciation, was an unexpected and significant source of Corg. Collectively, German seagrasses in the Baltic Sea are preventing 8.14 Mt of future CO2 emissions. Because Corg is mostly produced on-site and not imported from outside the meadow boundaries, the richness of this pool may be contingent on seagrass habitat health. Disturbance of this Corg stock could act as a source of CO2 emissions. However, the high spatial heterogeneity warrants site-specific investigations to obtain accurate estimates of blue carbon and a need to consider millennial timescale deposits of Corg beneath seagrass meadows in Germany and potentially other parts of the southwestern Baltic Sea.

Description: As one of Earth's most productive marine ecosystems, the Peruvian upwelling system transports large amounts of biogenic matter from the surface to the deep ocean. Whilst particle sinking velocity is a key factor controlling the biological pump, thereby affecting carbon sequestration and O2-depletion, it has not yet been measured in this system. During a 50 d mesocosm experiment in the surface waters off the coast of Peru, we assessed particle sinking velocities and their biogeochemical and physical drivers. We further characterized the general properties of exported particles under different phytoplankton communities and nutritional states. Average sinking velocities varied between size classes and ranged from 12.8 ± 0.7 m d−1 (particles 40–100 µm) to 19.4 ± 0.7 m d−1 (particles 100–250 µm) and 34.2 ± 1.5 m d−1 (particles 250–1000 µm) (± 95 % CI). Despite a distinct plankton succession from diatoms to dinoflagellates with concomitant 5-fold drop in opal ballasting, substantial changes in sinking velocity were not observed. This illustrates the complexity of counteracting factors driving the settling behaviour of marine particles. In contrast, we found higher sinking velocities with increasing particle size and roundness and decreasing porosity. Size had by far the strongest influence among these physical particle properties, despite a high amount of unexplained variability. Our study provides a detailed analysis of the drivers of particle sinking velocity in the Peruvian upwelling system, which allows modellers to optimize local particle flux parameterization. This will help to better project oxygen concentrations and carbon sequestration in a region that is subject to substantial climate-driven changes.

Description: Atlantic herring (Clupea harengus) has a complex population structure and displays a variety of reproductive strategies. Differences in reproductive strategies among herring populations are linked to their time of spawning, as well as to their reproductive investment which can be an indicator for migratory vs. stationary behavior. These differences are reflected in the number of oocytes (fecundity) and the size of the oocytes prior spawning. We studied potential mixing of herring with different reproductive strategies during the spring spawning season on a coastal spawning ground. It has been hypothesized that both spring and autumn spawning herring co-occur on this specific spawning ground. Therefore, we investigated the reproductive traits oocyte size, fecundity, fertilization success as well as length of the hatching larvae during the spring spawning season from February to April. We used a set of 11 single nucleotide polymorphism markers (SNPs), which are associated with spawning season, to genetically identify autumn and spring spawning herring. Reproductive traits were investigated separately within these genetically distinct spawning types. Furthermore, we used multivariate analyses to identify groups with potentially different reproductive strategies within the genetic spring spawners. Our results indicate that mixing between ripe spring and autumn spawners occurs on the spawning ground during spring, with ripe autumn spawners being generally smaller but having larger oocytes than spring spawners. Within spring spawners, we found large variability in reproductive traits. A following multivariate cluster analysis indicated two groups with different reproductive investment. Comparisons with other herring populations along the Norwegian coastline suggest that the high variability can be explained by the co-occurrence of groups with different reproductive investments potentially resulting from stationary or migratory behavior. Fertilization success and the length of the hatching larvae decreased with progression of the spawning season, with strong inter-individual variation, supporting our findings. Incorporating such complex population dynamics into management strategies of this species will be essential to build its future population resilience.

Description: Limiting global warming to well below 2°C and pursuing efforts to limit it to 1.5°C, as agreed in the 2015 Paris Agreement, requires global carbon neutrality by mid-century at the latest. The corresponding carbon budget is decreasing steadily and significantly. To phase out carbon emissions in line with the specified temperature target, countries are formulating their mitigation efforts in their long-term low greenhouse gas emission development strategies (LT-LEDS). However, there are no standardized specifications for preparing these strategies, which is why the reports published to date differ widely in terms of structure and scope. To consider the multiple facets of reaching net-zero from a systemic perspective as comprehensively as possible, the authors propose the Net-Zero-2050 System: A novel, transferrable systems approach that supports the development of national endeavors toward carbon neutrality. The Net-Zero-2050 System is defined by three interconnected components: The Carbon-Emission-Based System, the surrounding Framing System and a set of system boundaries. For both systems levels, IPCC approaches were used as a basis and were then adjusted and supplemented by Net-Zero-2050. We suggest applying the Net-Zero-2050 System—beyond the project environment—in carbon emission based contexts at different levels. Especially at the national level, this would improve the comparability of the different national strategies to achieve carbon neutrality.

Description: A growing population on a planet with limited resources demands finding new sources of protein. Hence, fisheries are turning their perspectives towards mesopelagic fish, which have, so far, remained relatively unexploited and poorly studied. Large uncertainties are associated with regards to their biomass, turn-over rates, susceptibility to environmental forcing and ecological and biogeochemical role. Models are useful to disentangle sources of uncertainties and to understand the impact of different processes on the biomass. In this study, we employed two food-web models – OSMOSE and the model by Anderson et al. (2019, or A2019) – coupled to a regional physical–biogeochemical model to simulate mesopelagic fish in the Eastern Tropical South Pacific ocean. The model by A2019 produced the largest biomass estimate, 26 to 130% higher than OSMOSE depending on the mortality parameters used. However, OSMOSE was calibrated to match observations in the coastal region off Peru and its temporal variability is affected by an explicit life cycle and food web. In contrast, the model by A2019 is more convenient to perform uncertainty analysis and it can be easily coupled to a biogeochemical model to estimate mesopelagic fish biomass. However, it is based on a flow analysis that had been previously applied to estimate global biomass of mesopelagic fish but has never been calibrated for the Eastern Tropical South Pacific. Furthermore, it assumes a steady-state in the energy transfer between primary production and mesopelagic fish, which may be an oversimplification for this highly dynamic system. OSMOSE is convenient to understand the interactions of the ecosystem and how including different life stages affects the model response. The combined strengths of both models allow us to study mesopelagic fish from a holistic perspective, taking into account energy fluxes and biomass uncertainties based on primary production, as well as complex ecological interactions.

Description: The subarctic front (SAF) in the pelagic North Pacific is the northernmost front of the Kuroshio-Oyashio transition zone separating the subpolar and subtropical gyres and is marked by a strong sea surface temperature gradient. A complex interplay of e.g. variations of currents, the wind system and other forcing mechanisms causes shifts of the SAF’s position on timescales from orbital to interannual. In this study, we present proxy data from the Emperor Seamount chain, which reveal a link between long-term ENSO (El Niño/Southern Oscillation) dynamics in the tropics and shifts of the SAF. Based on sediment core SO264-45-2 from Jimmu Seamount (46°33.792’N, 169°36.072’E) located close to the modern position of the SAF, we reconstruct changes in (sub)surface temperature ((sub)SST Mg/Ca ) and δ 18 O sw-ivc (approximating salinities) via combined Mg/Ca and δ 18 O analyses of the shallow-dwelling foraminifera Globigerina bulloides and the near-thermocline-dwelling Neogloboquadrina pachyderma , biological productivity (XRF-based Ba/Ti ratios), and terrigenous input via dust (XRF-based Fe). From ~600 to ~280 ka BP we observe significantly higher SST Mg/Ca than after an abrupt change at 280 ka BP. We assume that during this time warmer water from the Kuroshio-Oyashio transition zone reached the core site, reflecting a shift of the SAF from a position at or even north of our study site prior to 280 ka BP to a position south of our study site after 280 ka BP. We propose that such a northward displacement of the SAF between 600-280 ka BP was induced by sustained La Niña-like conditions, which led to increased transport of tropical ocean heat into the Kuroshio-Oyashio transition zone via the Kuroshio Current. After ~280 ka BP, the change to more El Niño-like conditions led to less heat transfer via the Kuroshio Current with the SAF remaining south of the core location. In contrast, our productivity record shows a clear glacial-interglacial pattern that is common in the North Pacific. We assume that this pattern is connected to changes in nutrient supply or utilization, which are not primarily driven by changes of the Kuroshio and Oyashio Currents or the SAF.

Description: El Niño–Southern Oscillation (ENSO) is a major source for teleconnections, including towards the tropical North Atlantic (TNA) region, whereby TNA sea surface temperatures (SSTs) are positively correlated with ENSO in boreal spring following an ENSO event. However, the Pacific–Atlantic connection can be impacted by different ENSO characteristics, such as the amplitude, location, and timing of Pacific SST anomalies (SSTAs). Indeed, the TNA SSTAs may respond nonlinearly to strong and extreme El Niño events. However, observational data for the number of extreme ENSO events remain limited, restricting our ability to investigate the influence of observed extreme ENSO events. To overcome this issue and to further evaluate the nonlinearity of the TNA SSTA response, two coupled climate models are used, namely the Community Earth System Model version 1 – Whole Atmosphere Community Climate Model (CESM-WACCM) and the Flexible Ocean and Climate Infrastructure version 1 (FOCI). In both models the TNA SSTAs respond linearly to ENSO during extreme El Niño events but nonlinearly to extreme La Niña events for CESM-WACCM. We investigate differences by using indices for all major mechanisms that connect ENSO to the TNA and compare them with reanalysis. CESM-WACCM and FOCI overall represent the teleconnection well, including that the tropical and extratropical pathways are similar to observations. Our results also show that a large portion of the nonlinearity during La Niña is explained by the interaction between Pacific SSTAs and the overlying upper-level divergence.

Description: The aim of this work is to present the food web models developed using the Ecopath with Ecosim (EwE) software tool to describe structure and functioning of various European marine ecosystems (eastern, central and western Mediterranean Sea; Black Sea; Bay of Biscay, Celtic Sea and Iberian coast; Baltic Sea; North Sea; English Channel, Irish Sea and west Scottish Sea; and Norwegian and Barents Seas). A total of 195 Ecopath models based on 168 scientific publications, which report original, updated and modified versions, were reviewed. Seventy models included Ecosim temporal simulations while 28 implemented Ecospace spatiotemporal dynamics. Most of the models and publications referred to the western Mediterranean Sea followed by the English Channel, Irish Sea and west Scottish Sea sub-regions. In the Mediterranean Sea, the western region had the largest number of models and publications, followed by the central and eastern regions; similar trends were observed in previous literature reviews. Most models addressed ecosystem functioning and fisheries-related hypotheses while several investigated the impact of climate change, the presence of alien species, aquaculture, chemical pollution, infrastructure, and energy production. Model complexity (i.e., number of functional groups) increased over time. Main forcing factors considered to run spatial and temporal simulations were trophic interactions, fishery, and primary production. Average scores of ecosystem indicators derived from the Ecopath summary statistics were compared. Uncertainty was also investigated based on the use of the Ecosampler plug-in and the Monte Carlo routine; only one third of the reviewed publications incorporated uncertainty analysis. Only a limited number of the models included the use of the ECOIND plug-in which provides the user with quantitative output of ecological indicators. We assert that the EwE modelling approach is a successful tool which provides a quantitative framework to analyse the structure and dynamics of ecosystems, and to evaluate the potential impacts of different management scenarios.

Description: Land degradation is a cause of many social, economic, and environmental problems. Therefore identification and monitoring of high-risk areas for land degradation are necessary. Despite the importance of land degradation due to wind and water erosion in some areas of the world, the combined study of both types of erosion in the same area receives relatively little attention. The present study aims to create a land degradation map in terms of soil erosion caused by wind and water erosion of semi-dry land. We focus on the Lut watershed in Iran, encompassing the Lut Desert that is influenced by both monsoon rainfalls and dust storms. Dust sources are identified using MODIS satellite images with the help of four different indices to quantify uncertainty. The dust source maps are assessed with three machine learning algorithms encompassing the artificial neural network (ANN), random forest (RF), and flexible discriminant analysis (FDA) to map dust sources paired with soil erosion susceptibility due to water. We assess the accuracy of the maps from the machine learning results with the area under the curve (AUC) of the receiver operating characteristic (ROC) metric. The water and aeolian soil erosion maps are used to identify different classes of land degradation risks. The results show that 43 % of the watershed is prone to land degradation in terms of both aeolian and water erosion. Most regions (45 %) have a risk of water erosion and some regions (7 %) a risk of aeolian erosion. Only a small fraction (4 %) of the total area of the region had a low to very low susceptibility for land degradation. The results of this study underline the risk of land degradation for in an inhabited region in Iran. Future work should focus on land degradation associated with soil erosion from water and storms in larger regions to evaluate the risks also elsewhere.

Description: Coralline algae play important ecological roles throughout the photic zone of the world’s oceans. Recent studies have shown that attached-living coralline algae can contain records of past climate variability. So far, algal-based paleo-reconstructions are mainly available from mid- to high-latitudes, while in low latitude and temperate regions only few examples exist. Here, we investigate samples from the attached-living encrusting coralline algal species Neogoniolithon hauckii (Rothpletz) R.A Townsend & Huisman [= Neogoniolithon mamillosum (Hauck) Setchell & L.R. Mason, nom illeg.] from a temperate site in the Mediterranean Sea to assess its potential as an environmental recorder. The specimens were collected at different water depths (20 and 40 m) in the Columbretes Islands (Spain). Sclerochronological analysis of sectioned samples revealed seasonal growth patterns. Mg/Ca, Li/Ca, and Ba/Ca ratios were measured in the algal skeletons using laser ablation inductively coupled mass spectrometry (LA-ICP-MS) in ultra-high resolution. We report a mean vertical extension rate of 1.1 to 1.2 mm/year (based on analysis of Mg/Ca cycles in 40 m and 20 m samples), representing the first growth rate measurement for this species. In addition, subannual banding patterns were mapped, measured, and could be linked to high frequency variability in laser-analyzed Mg/Ca ratios. Elemental ratios analyzed in Neogoniolithon hauckii were compared to in situ water temperatures measured at the water depth of sample collection. Our results show significant positive relationships between algal Mg/Ca (R=0.55) as well as Li/Ca ratios (R=0.46) and in situ measured temperature data (40 m specimen). Ba/Ca ratios show no significant correlation to temperature and may be influenced by other factors. These data suggest potential of this species for climate reconstructions in warm-temperate regions as Neogoniolithon hauckii is not only widely distributed in the Mediterranean, but also one of the few species that may be used as a temperature archive for mesophotic coralligenous assemblages that are strongly affected by the recent anthropogenic temperature rise. Further calibration studies are needed to test the element-temperature relationships on samples with longer growth records and in different settings and water depths.

Description: Little is known about dispersal in deep-sea ecosystems, especially for sponges, which are abundant ecosystem engineers. Understanding patterns of gene flow in deep-sea sponges is essential, especially in areas where rising pressure from anthropogenic activities makes difficult to combine management and conservation. Here, we combined population genomics and oceanographic modelling to understand how Northeast Atlantic populations (Cantabrian Sea to Norway) of the deep-sea sponge Phakellia ventilabrum are connected. The analysis comprised ddRADseq derived SNP datasets of 166 individuals collected from 57 sampling stations from 17 different areas, including two Marine Protected Areas, one Special Area of Conservation and other areas with different levels of protection. The 4,017 neutral SNPs used indicated high connectivity and panmixis amongst the majority of areas (Ireland to Norway), spanning ca. 2,500-km at depths of 99–900 m. This was likely due to the presence of strong ocean currents allowing long-distance larval transport, as supported by our migration analysis and by 3D particle tracking modelling. On the contrary, the Cantabrian Sea and Roscoff (France) samples, the southernmost areas in our study, appeared disconnected from the remaining areas, probably due to prevailing current circulation patterns and topographic features, which might be acting as barriers for gene flow. Despite this major genetic break, our results suggest that all protected areas studied are well-connected with each other. Interestingly, analysis of SNPs under selection replicated results obtained for neutral SNPs. The relatively low genetic diversity observed along the study area, though, highlights the potential fragility of this species to changing climates, which might compromise resilience to future threats.

Description: Sponges harbor diverse, specific, and stable microbial communities, but at the same time, they efficiently feed on microbes from the surrounding water column. This filter-feeding lifestyle poses the need to distinguish between three categories of bacteria: food to digest, symbionts to incorporate, and pathogens to eliminate. How sponges discriminate between these categories is still largely unknown. Phagocytosis is conceivable as the cellular mechanism taking part in such discrimination, but experimental evidence is missing. We developed a quantitative in-vivo phagocytosis assay using an emerging experimental model, the sponge Halichondria panicea. We incubated whole sponge individuals with different particles, recovered the sponge (host) cells, and tracked the incorporation of these particles into the sponge cells. Fluorescence-activated cell sorting (FACS) and fluorescent microscopy were used to quantify and verify phagocytic activity, defined here as the population of sponge cells with incorporated particles. Sponges were incubated with a green microalgae to test if particle concentration in the seawater affects the percentage of phagocytic activity, and to determine the timing where the maximum of phagocytic cells are captured in a pulse-chase experiment. Lastly, we investigated the application of our phagocytic assay with other particle types (i.e., fluorescently-labeled bacteria and fluorescent beads). The percentage of sponge cells that had incorporated algae, bacteria, and beads ranged between 5 to 24 %. These phagocytic sponge cells exhibited different morphologies and sizes depending on the type of particle presented to the sponge. Particle incorporation into sponge cells was positively related to algal concentration in the seawater, suggesting that sponge cells adjust their phagocytic activity depending on the number of particles they encounter. Our results further revealed that sponge phagocytosis initiates within minutes after exposure to the particles. Fluorescent and TEM microscopy rectified algal internalization and potential digestion in sponge cells. To our knowledge, this is the first quantitative in-vivo phagocytosis assay established in sponges that could be used to further explore phagocytosis as a cellular mechanism for sponges to differentiate between different microorganisms.

Description: Nitrous oxide (N2O) is a greenhouse gas, with a global warming potential 298 times that of carbon dioxide. Estuaries can be sources of N2O, but their emission estimates have significant uncertainties due to limited data availability and high spatiotemporal variability. We investigated the spatial and seasonal variability of dissolved N2O and its emissions along the Elbe Estuary (Germany), a well-mixed temperate estuary with high nutrient loading from agriculture. During nine research cruises performed between 2017 and 2022, we measured dissolved N2O concentrations, as well as dissolved nutrient and oxygen concentrations along the estuary, and calculated N2O saturations, flux densities, and emissions. We found that the estuary was a year-round source of N2O, with the highest emissions in winter when dissolved inorganic nitrogen (DIN) loads and wind speeds are high. However, in spring and summer, N2O saturations and emissions did not decrease alongside lower riverine nitrogen loads, suggesting that estuarine in situ N2O production is an important source of N2O. We identified two hotspot areas of N2O production: the Port of Hamburg, a major port region, and the mesohaline estuary near the maximum turbidity zone (MTZ). N2O production was fueled by the decomposition of riverine organic matter in the Hamburg Port and by marine organic matter in the MTZ. A comparison with previous measurements in the Elbe Estuary revealed that N2O saturation did not decrease alongside the decrease in DIN concentrations after a significant improvement of water quality in the 1990s that allowed for phytoplankton growth to re-establish in the river and estuary. The overarching control of phytoplankton growth on organic matter and, subsequently, on N2O production highlights the fact that eutrophication and elevated agricultural nutrient input can increase N2O emissions in estuaries.

Description: The warming climate is causing a strengthening of ocean stratification. Ocean stratification, in turn, has significant impacts on physical, biogeochemical and ecological processes, such as ocean circulation, ventilation, air-sea interactions, nutrient fluxes, primary productivity and fisheries. How these processes are affected in detail by changing stratification still remains uncertain and are likely to vary locally. Here, we investigate the state and trend of different parameters characterizing the stratification of the global upper-ocean which can be derived from Argo profiles for the period 2006-2021. Among those parameters are mixed layer depth, magnitude and depth of the vertical stratification maximum. The summertime stratification maximum has increased in both hemispheres, respectively. During wintertime, the stratification maximum has intensified in the Northern Hemisphere, while changes in the Southern Hemisphere have been relatively small. Comparisons to mixed layer characteristics show that a strengthening stratification is mainly accompanied by a warming and freshening of the mixed layer. In agreement with previous observational studies, we find a large-scale mixed layer deepening that regionally contributes to the increasing stratification. Globally, the vertical stratification maximum strengthens by 7-8% and the mixed layer deepens by 4 m during 2006-2021. This hints to an ongoing de-coupling of the surface ocean from the ocean interior. The investigated changes can help determine the origin of existing model-observation discrepancies and improve predictions on climate change impact on upper-ocean ecology and biogeochemistry.

Description: Eastern boundary upwelling systems (EBUS) contribute a disproportionate fraction of the global fish catch relative to their size and are especially susceptible to global environmental change. Here we present the evolution of communities over 50 days in an in situ mesocosm 6 km offshore of Callao, Peru and in the nearby unenclosed coastal Pacific Ocean. The communities were monitored using multi-marker environmental DNA (eDNA) metabarcoding and flow cytometry. DNA extracted from weekly water samples were subjected to amplicon sequencing for four genetic loci: 1) the V1-V2 region of the 16S rRNA gene, for photosynthetic eukaryotes (via their chloroplasts) and bacteria; 2) the V9 region of the 18S rRNA gene for exploration of eukaryotes but targeting phytoplankton; 3) cytochrome oxidase I (COI), for exploration of eukaryotic taxa but targeting invertebrates, and 4) the 12S rRNA gene, targeting vertebrates. The multi-marker approach showed a divergence of communities (from microbes to fish) between the mesocosm and the unenclosed ocean. Together with the environmental information, the genetic data furthered our mechanistic understanding of the processes that are shaping EBUS communities in a changing ocean. The unenclosed ocean experienced significant variability over the course of the 50-day experiment with temporal shifts in community composition but remained dominated by organisms that are characteristic of high nutrient, upwelling conditions (e.g. diatoms, copepods, anchovies). A large directional change was found in the mesocosm community. The mesocosm community that developed was characteristic of upwelling regions when upwelling relaxes and waters stratify (e.g. dinoflagellates, nanoflagellates). The selection of dinoflagellates under the warm (coastal El Niño) and stratified conditions in the mesocosm may be an indication of how EBUS will respond under the global environmental changes (i.e. continued global warming) forecast by the IPCC.

Description: Equatorial deep jets (EDJs) are vertically alternating, stacked zonal currents that flow along the Equator in all three ocean basins at intermediate depth. Their structure can be described quite well by the sum of high-baroclinic-mode equatorial Kelvin and Rossby waves. However, the EDJ meridional width is larger by a factor of 1.5 than inviscid theory predicts for such waves. Here, we use a set of idealised model configurations representing the Atlantic Ocean to investigate the contributions of different processes to the enhanced EDJ width. Corroborated by the analysis of shipboard velocity sections, we show that widening of the EDJs by momentum loss due to irreversible mixing or other processes contributes more to their enhanced time mean width than averaging over meandering of the jets. Most of the widening due to meandering can be attributed to the strength of intraseasonal variability in the jets' depth range, suggesting that the jets are meridionally advected by intraseasonal waves. A slightly weaker connection to intraseasonal variability is found for the EDJ widening by momentum loss. These results enhance our understanding of the dynamics of the EDJs and, more generally, of equatorial waves in the deep ocean.

Description: Tropical peat swamps are essential ecosystems, which provide numerous services, and also serve as a rich source of dissolved organic carbon (DOC), hydrogen ions and trace elements to peat draining rivers. However, not much is known about trace element export from tropical peat swamps. We investigated trace element dynamics in rivers and estuaries draining tropical peat swamps on Borneo, and examined the influence of estuarine processes as well as dissolved organic carbon (DOC) on the distribution and concentration of trace elements. Our results indicate acidic conditions (pH = 3.3) and high DOC concentration (3500 µmol L −1 ) at salinities〈1. We observed an initial release of trace elements at low salinity (0.05〈S〈 0.5), followed by scavenging to particles at intermediate salinities (0.5〈S〈10) due to an increasing ionic strength and pH. Peak concentrations (µmol kg −1 ) of Al (24.9), Si (96.2), Mn (4.9), Cu (0.035) and Ni (0.047) were observed during the dry season (July), and Fe concentrations (43.2) were highest during the wet season (December). We used the NICA-Donnan model to investigate the combined impact of DOC and pH on the formation of solid iron hydroxide (Fe(OH) 3 (s)). The Maludam river was predicted to be supersaturated for Fe hydroxides and the results affirmed our model prediction. The output showed Fe and Cu had a strong affinity for DOC and to a lesser extent Al and Ni in the conditions prevailing at the study sites. Statistical analyses also indicated strong correlation between Cu and Ni (r 2 = 0.97, 0.94 and 0.82) in Maludam, Sebuyau and Belait rivers and estuaries, respectively. The results obtained in this study are comparable to values published for southeast Asia and other continents for pristine peat draining rivers.

Description: The main goal of the study is to establish the spatial and temporal distribution of pyroclastic material from large explosive eruptions of the volcanoes of Kamchatka, the Kuril, and Aleutian Islands to create a generalized tephrochronological model and reveal patterns of explosive activity in this region. This paper presents new data on the composition of volcanic ash (tephra) found in the Pleistocene deposits of the northwestern Pacific from the eastern slope of the Detroit Rise (northwestern part of the Imperial Ridge), 450–550 km east of the Kamchatka Peninsula. Eleven layers and lenses of tephra aged from 28 to 245 ka, which were previously unknown, were studied in the core Lv63-4-2. Their stratigraphic position and age were determined based on age models developed in this study. Based on the geochemical composition of volcanic glass (determined using an electron microprobe), seven layers were correlated with tephra from several cores in the northwestern Pacific and the Bering Sea. The obtained results supplement the information on large explosive eruptions of volcanoes in the region and their periods of activity. They also allow the development of a generalized tephrochronological model of Quaternary deposits, which is necessary for stratigraphic correlation, and of paleooceanological and paleogeographic reconstructions.

Description: The abyssal plains are vast areas without large scale relief that occupy much of the ocean floor. Although long considered relatively featureless, they are now known to display substantial biological heterogeneity across different spatial scales. Ecological research in these regions benefits increasingly from non-destructive visual sampling of epifaunal organisms with imaging technology. We analysed images from ultra-high-definition towed camera transects at depths of around 3500 m across three stations (100–130 km apart) in the Bering Sea, to ask whether the density and distribution of visible epifauna indicated any substantial heterogeneity. We identified 71 different megafaunal taxa, of which 24 occurred at only one station. Measurements of the two most abundant faunal elements, the holothurian Elpidia minutissima and two xenophyophores morphotypes (the more common identifiable as Syringammina limosa), indicated significant differences in local densities and patchy aggregations that were strikingly dissimilar among stations. One station was dominated by xenophyophores, one was relatively depauperate in both target taxa as well as other identified megafauna, and the third station was dominated by Elpidia. This is an unexpected level of variation within comparable transects in a well-mixed oceanic basin, reinforcing the emerging view that abyssal habitats encompass biological heterogeneity at similar spatial scales to terrestrial continental realms.

Description: Ocean data synthesis products for specific biogeochemical essential ocean variables have the potential to facilitate today’s biogeochemical ocean data usage and comply with the Findable Accessible Interoperable and Reusable (FAIR) data principles. The products constitute key outputs from the Global Ocean Observation System, laying the observational foundation for information and services regarding climate and environmental status of the ocean. Using the Framework of Ocean Observing (FOO) readiness level concept, we present an evaluation framework for biogeochemical data synthesis products, which enables a systematic assessment of each product’s maturity. A new criteria catalog provides the foundation for assigning scores to the nine FOO readiness levels. As an example, we apply the assessment to four existing biogeochemical essential ocean variables data products. In descending readiness level order these are: The Surface Ocean CO2 Atlas (SOCAT); the Global Ocean Data Analysis Project (GLODAP); the MarinE MethanE and NiTrous Oxide (MEMENTO) data product and the Global Ocean Oxygen Database and ATlas (GO2DAT). Recognizing that the importance of adequate and comprehensive data from the essential ocean variables will grow, we recommend using this assessment framework to guide the biogeochemical data synthesis activities in their development. Moreover, we envision an overarching cross-platform FAIR biogeochemical data management system that sustainably supports the products individually and creates an integrated biogeochemical essential ocean variables data synthesis product; in short a system that provides truly comparable and FAIR data of the entire biogeochemical essential ocean variables spectrum.

Description: Deoxygenation of the ocean has been occurring over the last half century, particularly in poorly ventilated coastal waters. In coastal and estuarine environments, both the water column and sediments play key roles in controlling oxygen variability. In this study, we focus on controls of oxygen concentration in Bedford Basin (BB), a 70 m deep, seasonally hypoxic semi-enclosed fjord on the West Atlantic coast in Nova Scotia. The basin is connected to the Scotian Shelf via a narrow 20 m deep sill that restricts the resupply of bottom water. Hypoxia was recorded seasonally in 2018, 2019 and 2021 with minimum oxygen concentration of 5, 6.7 and 2.7 μM, respectively. Using a 1-D benthic-pelagic coupled model we investigate oxygen consumption and resupply processes during these years. The model was constrained with weekly water column measurements of temperature, salinity, chlorophyll-a fluorescence and dissolved oxygen from a monitoring station in the central basin together with seasonal measurements of benthic diffusive oxygen uptake. Our model suggested that 29-81%, and up to 36% of bottom water re-oxygenation occurred during the winter mixing period and through summer/fall intrusions of Scotian Shelf water, respectively. Occasional shelf water intrusions occurred rapidly, on a timescale of a few hours, and delivered equivalent amounts of oxygen as winter mixing and were sufficient to end bottom water hypoxia. Collectively, these mechanisms supplied the majority of the oxygen delivered to the bottom water. Oxygen supply to bottom waters during periods of water column stratification accounted for 19-36% of the annual flux. The mean benthic uptake was 12 ± 8 mmol m -2 d -1 and contributed ~20% of the total oxygen consumption below the sill depth. In 2021, sea surface temperature (SST) was unusually high and likely resulted in 50% less bottom water oxygenation compared to 2018 and 2019 due to increased stratification; SST in BB was found to be increasing at a rate of 0.11 ± 0.02 °C/year. Climate control on water column stratification are discussed and numerical experiments are used to compare the effects of different water column mixing scenarios on bottom water oxygenation.

Description: One of the major problems in the volcanic surveillance is how data from several techniques can be correlated and used to discriminate between possible precursors of volcanic eruptions and changes related to non-eruptive processes. Gas chemical surveys and measurements of SO 2 emission rates performed in the past (2006–2019) at Lastarria volcano in Northern Chile have revealed a persistent increment of magmatic sourced gas emissions since late November 2012, following a 13 years period of intense ground uplift. In this work, we provide new insights into the gas-chemical evolution of Lastarria’s fumarolic discharges obtained from direct sampling (2006–2019) and SO 2 emission rates using UV camera and DOAS instruments (2018–2019) and link these to pre-existing information on ground deformation (1998–2016) in order to determine the origin of observed degassing and ground deformation processes. We revise the four mechanisms originally proposed as alternatives by Lopez et al. (Geosphere, 2018, 14 (3), 983–1007) to explain the changes observed in the fluid geochemistry and ground deformation between 2009 and 2012, in order to explain major changes in gas-geochemistry over an extended period between 1998 and 2019. We hypothesize that a continuous sequence of processes explains the evolution in the fluid geochemistry of fumarolic discharges. Two mechanisms are responsible of the changes in the gas composition during the studied period, corresponding to a 1) deep magma chamber (7–15 km depth) pressurized by volatile exsolution (1998–2020), which is responsible of the large-scale deformation; followed by 2) a crystallization-induced degassing (2001–2020) and pressurization of the hydrothermal system (2003-early November 2012), where the former process induced the changes in the gas composition from hydrothermal-dominated to magmatic-dominated, whereas the last produced the small-scale deformation at Lastarria volcano. The changes in the gas composition since late November 2012, which were strongly dominated by magmatic volatiles, produced two consecutive processes: 1) acidification (late November 2012–2020) and 2) depletion (2019–2020) of the hydrothermal system. In this work we have shown that a long-term surveillance of the chemistry of fluid discharges provides valuable insights into underlying magmatic/volcanic processes, and consequently, for forecasting future eruptions.

Description: The Peruvian upwelling system is a highly productive ecosystem with a large oxygen minimum zone (OMZ) close to the surface. In this work, we carried out a mesocosm experiment off Callao, Peru, with the addition of water masses from the regional OMZ collected at two different sites simulating two different upwelling scenarios. Here, we focus on the pelagic remineralization of organic matter by the extracellular enzyme activity of leucine aminopeptidase (LAP) and alkaline phosphatase activity (APA). After the addition of the OMZ water, dissolved inorganic nitrogen (N) was depleted, but the standing stock of phytoplankton was relatively high, even after N depletion (mostly 〉 4 µg chlorophyll a L−1). During the initial phase of the experiment, APA was 0.6 nmol L−1 h−1 even though the PO concentration was 〉 0.5 µmol L−1. Initially, the dissolved organic phosphorus (DOP) decreased, coinciding with an increase in the PO concentration that was probably linked to the APA. The LAP activity was very high, with most of the measurements in the range of 200–800 nmol L−1 h−1. This enzyme hydrolyzes terminal amino acids from larger molecules (e.g., peptides or proteins), and these high values are probably linked to the highly productive but N-limited coastal ecosystem. Moreover, the experiment took place during a rare coastal El Niño event with higher than normal surface temperatures, which could have affected enzyme activity. Using a nonparametric multidimensional scaling analysis (NMDS) with a generalized additive model (GAM), we found that biogeochemical variables (e.g., nutrient and chlorophyll-a concentrations) and phytoplankton and bacterial communities explained up to 64 % of the variability in APA. The bacterial community best explained the variability (34 %) in LAP. The high hydrolysis rates for this enzyme suggest that pelagic N remineralization, likely driven by the bacterial community, supported the high standing stock of primary producers in the mesocosms after N depletion.