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: Marine methane hydrate is an ice-like substance that is stable in sediment around marine continental margins where water depths are greater than ~450–700 m. The release of methane due to melting of hydrates is considered to be a mechanism for past global carbon-cycle perturbations and could exacerbate ongoing anthropogenic climate change. Increases in bottom-water temperature at the landward limit of marine hydrate around continental margins, where vulnerable hydrate exists at or below the seabed, cause methane to vent into the ocean. However, this setting represents only ~3.5% of the global hydrate reservoir. The potential for methane from hydrate in deeper water to reach the atmosphere was considered negligible. Here we use three-dimensional (3D) seismic imagery to show that, on the Mauritanian margin, methane migrated at least 40 km below the base of the hydrate stability zone and vented through 23 pockmarks at the shelf break, probably during warmer Quaternary interglacials. We demonstrate that, under suitable circumstances, some of the 96.5% of methane bound in deeper water distal hydrates can reach the seafloor and vent into the ocean beyond the landward limit of marine hydrate. This reservoir should therefore be considered for estimating climate change-induced methane release during a warming world.

Description: Methylmercury is a potent toxin threatening the global population mainly through the consumption of marine fish. Hydrothermal venting directly delivers natural mercury to the ocean, yet its global flux remains poorly constrained. To determine the extent to which anthropogenic inputs have increased oceanic mercury levels, it is crucial to estimate natural mercury levels. Here we combine observations of vent fluids, plume waters, seawater and rock samples to quantify the release of mercury from the Trans-Atlantic Geotraverse hydrothermal vent at the Mid-Atlantic Ridge. The majority (67–95%) of the mercury enriched in the vent fluids (4,966 ± 497 pmol l −1 ) is rapidly diluted to reach background seawater levels (0.80 pmol l −1 ). A small Hg fraction (2.6–10%) is scavenged to the Trans-Atlantic Geotraverse mound rocks. Scaling up our findings and previous work, we propose a mercury flux estimate of 1.5–64.7 t per year from mid-ocean ridges. This hydrothermal flux is small in comparison to anthropogenic inputs. This suggests that most of the mercury present in the ocean must be of anthropogenic origin and that the implementation of emissions reduction measures outlined in the Minamata Convention could effectively reduce mercury levels in the global ocean and subsequently in marine fish.

Description: The Mediterranean Sea has been sampled irregularly by research vessels in the past, mostly by national expeditions in regional waters. To monitor the hydrographic, biogeochemical and circulation changes in the Mediterranean Sea, a systematic repeat oceanographic survey programme called Med-SHIP was recommended by the Mediterranean Science Commission (CIESM) in 2011, as part of the Global Ocean Ship-based Hydrographic Investigations Program (GO-SHIP). Med-SHIP consists of zonal and meridional surveys with different frequencies, where comprehensive physical and biogeochemical properties are measured with the highest international standards. The first zonal survey was done in 2011 and repeated in 2018. In addition, a network of meridional (and other key) hydrographic sections were designed: the first cycle of these sections was completed in 2016, with three cruises funded by the EU project EUROFLEETS2. This paper presents the physical and chemical data of the meridional and key transects in the Western and Eastern Mediterranean Sea collected during those cruises.

Description: Large explosive volcanic eruptions from island arcs pour pyroclastic currents into marine basins, impacting ecosystems and generating tsunamis that threaten coastal communities and infrastructures. Risk assessments require robust records of such highly hazardous events, which is challenging as most of the products lie buried under the sea. Here we report the discovery by IODP Expedition 398 of a giant rhyolitic pumice deposit emplaced 520 ± 10 ky ago at water depths of 200 to 1000 m during a high-intensity, shallow submarine eruption of ancestral Santorini Volcano. Pyroclastic currents discharged into the sea transformed into turbidity currents and slurries, forming a 〉89 ± 8 km 3 volcaniclastic megaturbidite up to 150 m thick in the surrounding marine basins, while breaching of the sea surface by the eruption column laid down veneers of ignimbrite on three islands. The eruption is one of the largest recorded on the South Aegean Volcanic Arc, and highlights the hazards from submarine explosive eruptions.

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: 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: Photosynthesis fuels primary production at the base of marine food webs. Yet, in many surface ocean ecosystems, diel-driven primary production is tightly coupled to daily loss. This tight coupling raises the question: which top-down drivers predominate in maintaining persistently stable picocyanobacterial populations over longer time scales? Motivated by high-frequency surface water measurements taken in the North Pacific Subtropical Gyre (NPSG), we developed multitrophic models to investigate bottom-up and top-down mechanisms underlying the balanced control of Prochlorococcus populations. We find that incorporating photosynthetic growth with viral- and predator-induced mortality is sufficient to recapitulate daily oscillations of Prochlorococcus abundances with baseline community abundances. In doing so, we infer that grazers in this environment function as the predominant top-down factor despite high standing viral particle densities. The model-data fits also reveal the ecological relevance of light-dependent viral traits and non-canonical factors to cellular loss. Finally, we leverage sensitivity analyses to demonstrate how variation in life history traits across distinct oceanic contexts, including variation in viral adsorption and grazer clearance rates, can transform the quantitative and even qualitative importance of top-down controls in shaping Prochlorococcus population dynamics.

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: Climate change is opening the Arctic Ocean to increasing human impact and ecosystem changes. Arctic fjords, the region’s most productive ecosystems, are sustained by a diverse microbial community at the base of the food web. Here we show that Arctic fjords become more prokaryotic in the picoplankton (0.2–3 µm) with increasing water temperatures. Across 21 fjords, we found that Arctic fjords had proportionally more trophically diverse (autotrophic, mixotrophic, and heterotrophic) picoeukaryotes, while subarctic and temperate fjords had relatively more diverse prokaryotic trophic groups. Modeled oceanographic connectivity between fjords suggested that transport alone would create a smooth gradient in beta diversity largely following the North Atlantic Current and East Greenland Current. Deviations from this suggested that picoeukaryotes had some strong regional patterns in beta diversity that reduced the effect of oceanographic connectivity, while prokaryotes were mainly stopped in their dispersal if strong temperature differences between sites were present. Fjords located in high Arctic regions also generally had very low prokaryotic alpha diversity. Ultimately, warming of Arctic fjords could induce a fundamental shift from more trophic diverse eukaryotic- to prokaryotic-dominated communities, with profound implications for Arctic ecosystem dynamics including their productivity patterns.

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: The sailboat Seaexplorer collected underway sea surface partial pressure of CO2 (pCO2) data for 129 days (2018–2021), including an Antarctic circumnavigation. By comparing ensembles of data-driven air-sea CO2 fluxes computed with and without sailboat data and applying a detection algorithm, we show that these sailboat observations significantly increase the regional carbon uptake in the North Atlantic and decrease it in the Southern Ocean. While compensating changes in both basins limit the global effect, the Southern Ocean–particularly frontal regions (40°S–60°S) during summertime—exhibited the largest air-sea CO2 flux changes, averaging 20% of the regional mean. Assessing the sensitivity of the air-sea CO2 flux to measurement uncertainty, the results stay robust within the expected random measurement uncertainty (± 5 μatm) but remain undetectable with a measurement offset of 5 µatm. We thus conclude that sailboats fill essential measurement gaps in remote ocean regions.

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: We present chromosome-level genome assemblies from representative species of three independently evolved seagrass lineages: Posidonia oceanica, Cymodocea nodosa, Thalassia testudinum and Zostera marina. We also include a draft genome of Potamogeton acutifolius, belonging to a freshwater sister lineage to Zosteraceae. All seagrass species share an ancient whole-genome triplication, while additional whole-genome duplications were uncovered for C. nodosa, Z. marina and P. acutifolius. Comparative analysis of selected gene families suggests that the transition from submerged-freshwater to submerged-marine environments mainly involved fine-tuning of multiple processes (such as osmoregulation, salinity, light capture, carbon acquisition and temperature) that all had to happen in parallel, probably explaining why adaptation to a marine lifestyle has been exceedingly rare. Major gene losses related to stomata, volatiles, defence and lignification are probably a consequence of the return to the sea rather than the cause of it. These new genomes will accelerate functional studies and solutions, as continuing losses of the ‘savannahs of the sea’ are of major concern in times of climate change and loss of biodiversity.

Description: The relationship between energy reserves of cold-water corals (CWCs) and their physiological performance remains largely unknown. In addition, it is poorly understood how the energy allocation to different metabolic processes might change with projected decreasing food supply to the deep sea in the future. This study explores the temporal and spatial variations of total energy reserves (proteins, carbohydrates and lipids) of the CWC Desmophyllum dianthus and their correlation with its calcification rate. We took advantage of distinct horizontal and vertical physico-chemical gradients in Comau Fjord (Chile) and examined the changes in energy reserves over one year in an in situ reciprocal transplantation experiment (20 m vs. 300 m and fjord head vs. mouth). Total energy reserves correlated positively with calcification rates. The fast-growing deep corals had higher and less variable energy reserves, while the slower-growing shallow corals showed pronounced seasonal changes in energy reserves. Novel deep corals (transplanted from shallow) were able to quickly increase both their calcification rates and energy reserves to similar levels as native deep corals. Our study shows the importance of energy reserves in sustaining CWC growth in spite of aragonite undersaturated conditions (deep corals) in the present, and potentially also future ocean.

Description: Caldera-forming eruptions of silicic volcanic systems are among the most devastating events on Earth. By contrast, post-collapse volcanic activity initiating new caldera cycles is generally considered less hazardous. Formed after Santorini’s latest caldera-forming eruption of ~1600 bce , the Kameni Volcano in the southern Aegean Sea enables the eruptive evolution of a recharging multi-cyclic caldera to be reconstructed. Kameni’s eruptive record has been documented by onshore products and historical descriptions of mainly effusive eruptions dating back to 197 bce . Here we combine high-resolution seismic reflection data with cored lithologies from International Ocean Discovery Program Expedition 398 at four sites to determine the submarine architecture and volcanic history of intra-caldera deposits from Kameni. Our shore-crossing analysis reveals the deposits of a submarine explosive eruption that produced up to 3.1 km 3 of pumice and ash, which we relate to a historical eruption in 726 ce . The estimated volcanic explosivity index of magnitude 5 exceeds previously considered worst-case eruptive scenarios for Santorini. Our finding that the Santorini caldera is capable of producing large explosive eruptions at an early stage in the caldera cycle implies an elevated hazard potential for the eastern Mediterranean region, and potentially for other recharging silicic calderas.

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: A natural plankton community from oligotrophic subtropical waters of the Atlantic near Gran Canaria, Spain, was subjected to varying degrees of ocean alkalinity enhancement (OAE) to assess the potential physiological effects, in the context of the application of ocean carbon dioxide removal (CDR) techniques. We employed 8.3 m3 mesocosms with a sediment trap attached to the bottom, creating a gradient in total alkalinity (TA). The lowest point on this gradient was 2400 μmol · L-1, which corresponded to the natural alkalinity of the environment, and the highest point was 4800 μmol · L-1. Over the course of the 33-day experiment, the plankton community exhibited two distinct phases. In phase-I (days 5–20), a notable decline in the photosynthetic efficiency (Fv/Fm) was observed. This change was accompanied by substantial reductions in the abundances of picoeukaryotes, small size nanoeukaryotes (nanoeukaryotes-1), and microplankton. The cell viability of picoeukaryotes, as indicated by fluorescein-di-acetate hydrolysis by cellular esterases (FDA- green fluorescence), slightly increased by the end of phase-I whilst the viability of nanoeukaryotes 1 and Synechococcus spp . did not change. Reactive oxygen species levels (ROS-green fluorescence) showed no significant changes for any of the functional groups. In contrast, in phase-II (days 21–33), a pronounced community response was observed. Increases in Fv/Fm in the intermediate OAE treatments of ∆900 to ∆1800 μmol · L-1 and in chlorophyll-a (Chl-a), chlorophyll-c2 (Chl-c2) , fucoxanthin and divinyl-Chl-a were attributed to the emergence of blooms of large size nanoeukaryotes (nanoeukaryotes-2) from the genera Chrysochromulina, as well as picoeukaryotes. Synechococcus spp. also flourished towards the end of this phase. In parallel, we observed a total 20 % significant change in the metaproteome of the phytoplankton community. This is considered a significant alteration in protein expression, having substantial impacts on cellular functions and the physiology of the organisms. Medium levels of ∆TA showed more upregulated and less downregulated proteins than higher ∆TA treatments. Under these conditions, cell viability significantly increased in pico and nanoeukaryotes-1 in intermediate alkalinity levels, while in Synechococcus spp., nanoeukaryotes-2 and microplankton remained stable. ROS levels did not significantly change in any functional group. The pigment ratios DD+DT : FUCO, and DD+DT : Chl-a increased in medium ∆TA treatments, supporting the idea of nutrient deficiency alleviation and the absence of physiological stress. Taken all data together, this study shows that there is minimal evidence indicating a harmful impact of high alkalinity on the plankton community. The OAE treatments did not result in physiological fitness impairment, thus OAE did not cause cellular stress in the phytoplankton community studied.

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: In September 2021, volcanic aerosol (mainly freshly formed sulfate plumes) originating from the eruption of Cumbre Vieja on La Palma, Canary Islands, Spain, crossed Cabo Verde at altitudes below 2 km. On 24 September 2021, an extraordinary large aerosol optical depth (AOD) close to 1 (daily mean at 500 nm) was observed at Mindelo, Cabo Verde. This event provided favorable conditions to obtain lidar-derived profiles of extinction and backscatter coefficients, lidar ratio, and depolarization ratio at 355, 532 and 1064 nm in the sulfate aerosol plume. A novel feature of the lidar system operated at Mindelo is the availability of extinction, lidar ratio and depolarization measurements at 1064 nm in addition to the standard wavelengths of 355 and 532 nm. Having measurements of these parameters at all three wavelengths is a major advantage for the aerosol characterization and in aerosol typing efforts as the lidar ratio and the particle linear depolarization ratio are key parameters for this purpose. In this article, we present the key results of the lidar observations obtained on one specific day, namely on 24 September 2021 at 04:38-05:57 UTC, including the first ever measurements of the particle extinction coefficient, the lidar ratio and the depolarization ratio at 1064 nm for volcanic sulfate, and discuss the findings in terms of aerosol optical properties and mass concentrations by comparison with a reference observation (16 September 2021) representing the typical background conditions before the start of the eruptions. We found an unusual high particle extinction coefficient of 721 +/- 51, 549 +/- 38 and 178 +/- 13 Mm - 1 , as well as an enhanced lidar ratio of 66.9 +/- 10.1, 60.2 +/- 9.2 and 30.8 +/- 8.7 sr at 355, 532 and 1064 nm, respectively, in the sulfate-dominated planetary boundary layer (PBL). The particle linear depolarization ratio was 〈= 0.9 % at all respective wavelengths. It is the first time that lidar-derived intensive aerosol optical properties could be derived for volcanic sulfate at all three wavelengths, and thus it is a highly valuable data set for global aerosol characterization. The lidar analysis also revealed a sulfate-related AOD of about 0.35 +/- 0.03 at 532 nm of the total PBL-related AOD of 0.43. The rest of the AOD contribution was caused by a lofted Saharan dust layer extending from 1.4 to 5 km and leading to a total AOD of 0.79 at 532 nm. Volcanic ash contribution to the observed aerosol plumes could be mostly excluded based on trajectory analysis and the observed optical properties. Peak mass concentration was 178.5 +/- 44.6 mu g m - 3 in the volcanic-influenced and sulfate-dominated polluted PBL, showing the hazardous potential of such sulfate plumes to significantly worsen local air quality even at remote locations.

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: 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: Thriving in both epipelagic and mesopelagic layers, Rhizaria are biomineralizing protists, mixotrophs or flux-feeders, often reaching gigantic sizes. In situ imaging showed their contribution to oceanic carbon stock, but left their contribution to element cycling unquantified. Here, we compile a global dataset of 167,551 Underwater Vision Profiler 5 Rhizaria images, and apply machine learning models to predict their organic carbon and biogenic silica biomasses in the uppermost 1000 m. We estimate that Rhizaria represent up to 1.7% of mesozooplankton carbon biomass in the top 500 m. Rhizaria biomass, dominated by Phaeodaria, is more than twice as high in the mesopelagic than in the epipelagic layer. Globally, the carbon demand of mesopelagic, flux-feeding Phaeodaria reaches 0.46 Pg C y −1 , representing 3.8 to 9.2% of gravitational carbon export. Furthermore, we show that Rhizaria are a unique source of biogenic silica production in the mesopelagic layer, where no other silicifiers are present. Our global census further highlights the importance of Rhizaria for ocean biogeochemistry.

Description: Observations of spatio-temporal variability of the deep ocean are rare and little is known about occurrence of deep ocean mesoscale dynamics. Here, we make use of 2.5 years of time series data from three distributed sensor arrays, which acquired high-resolution temperature, pressure and sound speed data of the bottom layer offshore northern Chile. Estimating salinity and density from the direct observations enable access to the full spectrum of hydrographic variability from a multi-hourly to annual time scale and with average inter-station distances of less than 1 km. Analyses revealed interannual warming over the continental slope of 0.002 °C yr−1–0.003 °C yr−1, and could trace periodic hydrographic anomalies, likely related to coastal-trapped waves, as far as to the lower continental slope. A concurrent change in the shape of the warm anomalies and the rate of deep-sea warming that occurs with the crossing of the deep-sea trench suggests that the abyssal part of the eastern boundary current system off Chile does not extend past the deep sea trench. Furthermore, the comparison of anomaly timing and shape in between stations implies southwards flow over the mid to lower continental slope, centred closer to the trench.

Description: Enhancing ocean productivity by artificial upwelling is evaluated as a nature-based solution for food security and climate change mitigation. Fish production is intended through diatom-based plankton food webs as these are assumed to be short and efficient. However, our findings from mesocosm experiments on artificial upwelling in the oligotrophic ocean disagree with this classical food web model. Here, diatoms did not reduce trophic length and instead impaired the transfer of primary production to crustacean grazers and small pelagic fish. The diatom-driven decrease in trophic efficiency was likely mediated by changes in nutritional value for the copepod grazers. Whilst diatoms benefitted the availability of essential fatty acids, they also caused unfavorable elemental compositions via high carbon-to-nitrogen ratios (i.e. low protein content) to which the grazers were unable to adapt. This nutritional imbalance for grazers was most pronounced in systems optimized for CO 2 uptake through carbon-to-nitrogen ratios well beyond Redfield. A simultaneous enhancement of fisheries production and carbon sequestration via artificial upwelling may thus be difficult to achieve given their opposing stoichiometric constraints. Our study suggest that food quality can be more critical than quantity to maximize food web productivity during shorter-term fertilization of the oligotrophic ocean.

Description: Turbulent mixing in the ocean, lakes and reservoirs facilitates the transport of momentum, heat, nutrients, and other passive tracers. Turbulent fluxes are proportional to the rate of turbulent kinetic energy dissipation per unit mass, ε. A common method for ε measurements is using microstructure profilers with shear probes. Such measurements are now widespread, and a non-expert practitioner will benefit from best practice guidelines and benchmark datasets. As a part of the Scientific Committee on Oceanographic Research (SCOR) working group on “Analysing ocean turbulence observations to quantify mixing” (ATOMIX), we compiled a collection of five benchmark data of ε from measurements of turbulence shear using shear probes. The datasets are processed using the ATOMIX recommendations for best practices documented separately. Here, we describe and validate the datasets. The benchmark collection is from different types of instruments and covers a wide range of environmental conditions. These datasets serve to guide the users to test their ε estimation methods and quality-assurance metrics, and to standardize their data for archiving.

Description: The evolution of the northern hemispheric climate during the last glacial period was beset by quasi-episodic iceberg discharge events from the Laurentide ice sheet, known as Heinrich events (HEs). The paleo record places most HEs into the cold stadial of the Dansgaard-Oeschger cycle. However, not every Dansgaard-Oeschger cycle is associated with a HE, revealing a complex interplay between the two modes of glacial variability. Here, using a coupled ice sheet-solid earth model, we introduce a mechanism that explains the synchronicity of HEs and Dansgaard-Oeschger cycles. Unlike earlier studies, our mechanism does not require a trigger during the stadial. Instead, the atmospheric warming signal during the interstadial of the Dansgaard-Oeschger cycle causes enhanced ice stream thickening that leads to the HE during the late interstadial. We demonstrate that this mechanism reproduces the key HE characteristics and provides an explanation for synchronous HEs from different regions of the Laurentide ice sheet.

Description: There has been extensive research into the nonlinear responses of the Earth system to astronomical forcing during the last glacial cycle. However, the speed and spatial geometry of ice sheet expansion to its largest extent at the Last Glacial Maximum 21 thousand years ago remains uncertain. Here we use an Earth system model with interactive ice sheets to show that distinct initial North American (Laurentide) ice sheets at 38 thousand years ago converge towards a configuration consistent with the Last Glacial Maximum due to feedbacks between atmospheric circulation and ice sheet geometry. Notably, ice advance speed and spatial pattern in our model are controlled by the amount of summer snowfall, which is dependent on moisture transport pathways from the North Atlantic warm pool linked to ice sheet geometry. The consequence of increased summer snowfall on the surface mass balance of the ice sheet is not only the direct increase in accumulation but the indirect reduction in melt through the snow/ice–albedo feedback. These feedbacks provide an effective mechanism for ice growth for a range of initial ice sheet states and may explain the rapid North American ice volume increase during the last ice age and potentially driving growth during previous glacial periods.

Description: In this paper we describe the implementation of the carbon isotopes 13C and 14C (radiocarbon) into the marine biogeochemistry model REcoM3. The implementation is tested in long-term equilibrium simulations where REcoM3 is coupled with the ocean general circulation model FESOM2.1, applying a low-resolution configuration and idealized climate forcing. Focusing on the carbon-isotopic composition of dissolved inorganic carbon (δ13CDIC and Δ14CDIC), our model results are largely consistent with reconstructions for the pre-anthropogenic period. Our simulations also exhibit discrepancies, e.g. in upwelling regions and the interior of the North Pacific. Some of these differences are due to the limitations of our ocean circulation model setup, which results in a rather shallow meridional overturning circulation. We additionally study the accuracy of two simplified modelling approaches for dissolved inorganic 14C, which are faster (15 % and about a factor of five, respectively) than the complete consideration of the marine radiocarbon cycle. The accuracy of both simplified approaches is better than 5 %, which should be sufficient for most studies of Δ14CDIC.

Description: Laboratory experiments showed that the isotopic fractionation of δ13C and of δ18O during calcite formation of planktic foraminifera are species-specific functions of ambient CO concentration. This effect became known as the carbonate ion effect (CIE), whose role for the interpretation of marine sediment data will be investigated here in an in-depth analysis of the 13C cycle. For this investigation, we constructed new 160 kyr long mono-specific stacks of changes in both δ13C and δ18O from either the planktic foraminifera Globigerinoides ruber (rub) or Trilobatus sacculifer (sac) from 112 and 40 marine records, respectively, from the wider tropics (latitudes below 38°). Both mono-specific time series Δ(δ13Crub) and Δ(δ13Csac) are very similar to each other, and a linear regression through a scatter plot of both data sets has a slope of ∼ 0.99 – although the laboratory-based CIE for both species differs by a factor of nearly 2, implying that they should record distinctly different changes in δ13C, if we accept that the carbonate ion concentration changes on glacial–interglacial timescales. For a deeper understanding of the 13C cycle, we use the Solid Earth version of the Box model of the Isotopic Carbon cYCLE (BICYLE-SE) to calculate how surface-ocean CO should have varied over time in order to be able to calculate the potential offsets which would by caused by the CIE quantified in culture experiments. Our simulations are forced with atmospheric reconstructions of CO2 and δ13CO2 derived from ice cores to obtain a carbon cycle which should at least at the surface ocean be as close as possible to expected conditions and which in the deep ocean largely agrees with the carbon isotope ratio of dissolved inorganic carbon (DIC), δ13CDIC, as reconstructed from benthic foraminifera. We find that both Δ(δ13Crub) and Δ(δ13Csac) agree better with changes in simulated δ13CDIC when ignoring the CIE than those time series which were corrected for the CIE. The combination of data- and model-based evidence for the lack of a role for the CIE in Δ(δ13Crub) and Δ(δ13Csac) suggests that the CIE as measured in laboratory experiments is not directly transferable to the interpretation of marine sediment records. The much smaller CIE-to-glacial–interglacial-signal ratio in foraminifera δ18O, when compared to δ13C, prevents us from drawing robust conclusions on the role of the CIE in δ18O as recorded in the hard shells of both species. However, theories propose that the CIE in both δ13C and δ18O depends on the pH in the surrounding water, suggesting that the CIE should be detectable in neither or both of the isotopes. Whether this lack of role of the CIE in the interpretation of planktic paleo-data is a general feature or is restricted to the two species investigated here needs to be checked with further data from other planktic foraminiferal species.

Description: An increasing number of climate model simulations is becoming available for the transition from the Last Glacial Maximum to the Holocene. Assessing the simulations' reliability requires benchmarking against environmental proxy records. To date, no established method exists to compare these two data sources in space and time over a period with changing background conditions. Here, we develop a new algorithm to rank simulations according to their deviation from reconstructed magnitudes and temporal patterns of orbital and millennial-scale temperature variations. The use of proxy forward modeling allows for accounting for non-climatic processes that affect the temperature reconstructions. It further avoids the need to reconstruct gridded fields or regional mean temperature time series from sparse and uncertain proxy data. First, we test the reliability and robustness of our algorithm in idealized experiments with prescribed deglacial temperature histories. We quantify the influence of limited temporal resolution, chronological uncertainties, and non-climatic processes by constructing noisy pseudo-proxies. While model–data comparison results become less reliable with increasing uncertainties, we find that the algorithm discriminates well between simulations under realistic non-climatic noise levels. To obtain reliable and robust rankings, we advise spatial averaging of the results for individual proxy records. Second, we demonstrate our method by quantifying the deviations between an ensemble of transient deglacial simulations and a global compilation of sea surface temperature reconstructions. The ranking of the simulations differs substantially between the considered regions and timescales, which suggests that optimizing for agreement with the temporal patterns of a small set of proxies might be insufficient for capturing the spatial structure of the deglacial temperature variability. We attribute the diversity in the rankings to more regionally confined temperature variations in reconstructions than in simulations, which could be the result of uncertainties in boundary conditions, shortcomings in models, or regionally varying characteristics of reconstructions such as recording seasons and depths. Future work towards disentangling these potential reasons can leverage the flexible design of our algorithm and its demonstrated ability to identify varying levels of model–data agreement. Additionally, the algorithm can be applied to variables like oxygen isotopes and climate transitions such as the penultimate deglaciation and the last glacial inception.

Description: Oxygen isotopes in biogenic silica (δ18OBSi) from lake sediments allow for quantitative reconstruction of past hydroclimate and proxy-model comparison in terrestrial environments. The signals of individual records have been attributed to different factors, such as air temperature (Tair), atmospheric circulation patterns, hydrological changes, and lake evaporation. While every lake has its own local set of drivers of δ18O variability, here we explore the extent to which regional or even global signals emerge from a series of paleoenvironmental records. This study provides a comprehensive compilation and combined statistical evaluation of the existing lake sediment δ18OBSi records, largely missing in other summary publications (i.e. PAGES network). For this purpose, we have identified and compiled 71 down-core records published to date and complemented these datasets with additional lake basin parameters (e.g. lake water residence time and catchment size) to best characterize the signal properties. Records feature widely different temporal coverage and resolution, ranging from decadal-scale records covering the past 150 years to records with multi-millennial-scale resolution spanning glacial–interglacial cycles. The best coverage in number of records (N = 37) and data points (N = 2112) is available for Northern Hemispheric (NH) extratropical regions throughout the Holocene (roughly corresponding to Marine Isotope Stage 1; MIS 1). To address the different variabilities and temporal offsets, records were brought to a common temporal resolution by binning and subsequently filtered for hydrologically open lakes with lake water residence times 〈 100 years. For mid- to high-latitude (〉 45° N) lakes, we find common δ18OBSi patterns among the lake records during both the Holocene and Common Era (CE). These include maxima and minima corresponding to known climate episodes, such as the Holocene Thermal Maximum (HTM), Neoglacial Cooling, Medieval Climate Anomaly (MCA) and the Little Ice Age (LIA). These patterns are in line with long-term air temperature changes supported by previously published climate reconstructions from other archives, as well as Holocene summer insolation changes. In conclusion, oxygen isotope records from NH extratropical lake sediments feature a common climate signal at centennial (for CE) and millennial (for Holocene) timescales despite stemming from different lakes in different geographic locations and hence constitute a valuable proxy for past climate reconstructions.

Description: We present transient simulations of the last glacial inception using the Earth system model CLIMBER-X with dynamic vegetation, interactive ice sheets, and visco-elastic solid Earth responses. The simulations are initialized at the middle of the Eemian interglacial (125 kiloyears before present, ka) and run until 100 ka, driven by prescribed changes in Earth's orbital parameters and greenhouse gas concentrations from ice core data. CLIMBER-X simulates a rapid increase in Northern Hemisphere ice sheet area through MIS5d, with ice sheets expanding over northern North America and Scandinavia, in broad agreement with proxy reconstructions. While most of the increase in ice sheet area occurs over a relatively short period between 119 and 117 ka, the larger part of the increase in ice volume occurs afterwards with an almost constant ice sheet extent. We show that the vegetation feedback plays a fundamental role in controlling the ice sheet expansion during the last glacial inception. In particular, with prescribed present-day vegetation the model simulates a global sea level drop of only ∼ 20 m, compared with the ∼ 35 m decrease in sea level with dynamic vegetation response. The ice sheet and carbon cycle feedbacks play only a minor role during the ice sheet expansion phase prior to ∼ 115 ka but are important in limiting the deglaciation during the following phase characterized by increasing summer insolation. The model results are sensitive to climate model biases and to the parameterization of snow albedo, while they show only a weak dependence on changes in the ice sheet model resolution and the acceleration factor used to speed up the climate component. Overall, our simulations confirm and refine previous results showing that climate–vegetation–cryosphere feedbacks play a fundamental role in the transition from interglacial to glacial states characterizing Quaternary glacial cycles.

Description: Weather causes extremes in photovoltaic and wind power production. Here we present a comprehensive climatology of anomalies in photovoltaic and wind power production associated with weather patterns in Europe considering the 2019 and potential 2050 installations, and hourly to ten-day events. To that end, we performed kilometer-scale numerical simulations of hourly power production for 23 years and paired the output with a weather classification which allows a detailed assessment of weather-driven spatio-temporal production anomalies. Our results highlight the dependency of low-power production events on the installed capacities and the event duration. South-shifted Westerlies (Anticyclonic South-Easterlies) are associated with the lowest hourly (ten-day) extremes for the 2050 (both) installations. Regional power production anomalies can differ from the ones in the European mean. Our findings suggest that weather patterns can serve as indicators for expected photovoltaic and wind power production anomalies and may be useful for early warnings in the energy sector.

Description: Phytoplankton form the base of the marine food web by transforming CO2 into organic carbon via photosynthesis. Despite the importance of phytoplankton for marine ecosystems and global carbon cycling, 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 total (depth-integrated) phytoplankton biomass in the Southern Ocean is projected to largely remain unchanged under climate change by the Coupled Model Intercomparison Project Phase 6 (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 on average to improve growth conditions, consequently weaken bottom-up control, and confine phytoplankton closer to the surface. An increase in the phytoplankton concentration promotes zooplankton grazing efficiency, thus intensifying top-down control. However, large differences across the model ensemble exist, with some models simulating a decrease in surface phytoplankton concentrations. To reduce uncertainties in projections of surface phytoplankton concentrations, we employ an emergent constraint approach using the observed sensitivity of surface chlorophyll concentration, taken as an observable proxy for phytoplankton, to seasonal changes in the mixed-layer depth as an indicator for future changes in surface phytoplankton concentrations. The emergent constraint reduces uncertainties in surface phytoplankton concentration projections by around one-third and increases confidence that surface phytoplankton concentrations will indeed rise due to shoaling mixed layers under global warming, thus favouring intensified 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 the phytoplankton response to climate change in the Southern Ocean.

Description: To better understand possible reasons for the diverse modeling results and large discrepancies of the detected solar fingerprints, we took one step back and assessed the "initial" solar signals in the middle atmosphere based on large ensemble simulations with multiple climate models — FOCI, EMAC, and MPI-ESM-HR. Consistent with previous work, we find that the 11-year solar cycle signals in the short wave heating rate (SWHR) and ozone anomalies are robust and statistically significant in all three models. These "initial" solar cycle signals in SWHR, ozone, and temperature anomalies are sensitive to the strength of the solar forcing. Correlation coefficients of the solar cycle with the SWHR, ozone, and temperature anomalies linearly increase along with the enhancement of the solar cycle amplitude, and this reliance becomes more complex when the solar cycle amplitude exceeds a certain threshold. In addition, the cold bias in the tropical stratopause of EMAC dampens the subsequent results of the "initial" solar signal. The warm pole bias in MPI-ESM-HR leads to a weak polar night jet (PNJ), which may limit the top-down propagation of the initial solar signal. Although FOCI simulated a so-called top-down response as revealed in previous studies in a period with large solar cycle amplitudes, its warm bias in the tropical upper stratosphere results in a positive bias in PNJ and can lead to a "reversed" response in some extreme cases. We suggest a careful interpretation of the single model result and further re-examination of the solar signal based on more climate models.

Description: In the framework of a changing climate, it is useful to devise methods capable of effectively assessing and monitoring the changing landscape of air-sea CO2 fluxes. In this study, we developed an integrated machine learning tool to objectively classify and track marine carbon biomes under seasonally and interannually changing environmental conditions. The tool was applied to the monthly output of a global ocean biogeochemistry model at 0.25° resolution run under atmospheric forcing for the period 1958–2018. Carbon biomes are defined as regions having consistent relations between surface CO2 fugacity (fCO2) and its main drivers (temperature, dissolved inorganic carbon, alkalinity). We detected carbon biomes by using an agglomerative hierarchical clustering (HC) methodology applied to spatial target-driver relationships, whereby a novel adaptive approach to cut the HC dendrogram based on the compactness and similarity of the clusters was employed. Based only on the spatial variability of the target-driver relationships and with no prior knowledge on the cluster location, we were able to detect well-defined and geographically meaningful carbon biomes. A deep learning model was constructed to track the seasonal and interannual evolution of the carbon biomes, wherein a feed-forward neural network was trained to assign labels to detected biomes. We find that the area covered by the carbon biomes responds robustly to seasonal variations in environmental conditions. A seasonal alternation between different biomes is observed over the North Atlantic and Southern Ocean. Long-term trends in biome coverage over the 1958–2018 period, namely a 10 % expansion of the subtropical biome in the North Atlantic and a 10 % expansion of the subpolar biome in the Southern Ocean, are suggestive of long-term climate shifts. Our approach thus provides a framework that can facilitate the monitoring of the impacts of climate change on the ocean carbon cycle and the evaluation of carbon cycle projections across Earth System Models.

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: 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: 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: 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 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: 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: 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: 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: 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: 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: The South Atlantic Meridional Overturning Circulation initiative began as a grassroots effort to study the South Atlantic Ocean and its impact on climate. In striving towards this goal, it has also become a platform for the empowerment of women and international scientists.

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: Phytoplankton forms the base of aquatic food webs and element cycling in diverse aquatic systems. The fate of phytoplankton-derived organic matter, however, often remains unresolved as it is controlled by complex, interlinked remineralization and sedimentation processes. We here investigate a rarely considered control mechanism on sinking organic matter fluxes: fungal parasites infecting phytoplankton. We demonstrate that bacterial colonization is promoted 3.5-fold on fungal-infected phytoplankton cells in comparison to non-infected cells in a cultured model pathosystem (diatom Synedra, fungal microparasite Zygophlyctis, and co-growing bacteria), and even ≥17-fold in field-sampled populations (Planktothrix, Synedra, and Fragilaria). Additional data obtained using the Synedra–Zygophlyctis model system reveals that fungal infections reduce the formation of aggregates. Moreover, carbon respiration is 2-fold higher and settling velocities are 11–48% lower for similar-sized fungal-infected vs. non-infected aggregates. Our data imply that parasites can effectively control the fate of phytoplankton-derived organic matter on a single-cell to single-aggregate scale, potentially enhancing remineralization and reducing sedimentation in freshwater and coastal systems.

Description: In search for critical elements, polymetallic nodules at the deep abyssal seafloor are targeted for mining operations. Nodules efficiently scavenge and retain several naturally occurring uranium-series radioisotopes, which predominantly emit alpha radiation during decay. Here, we present new data on the activity concentrations of thorium-230, radium-226, and protactinium-231, as well as on the release of radon-222 in and from nodules from the NE Pacific Ocean. In line with abundantly published data from historic studies, we demonstrate that the activity concentrations for several alpha emitters are often higher than 5 Bq g −1 at the surface of the nodules. These observed values can exceed current exemption levels by up to a factor of 1000, and even entire nodules commonly exceed these limits. Exemption levels are in place for naturally occurring radioactive materials (NORM) such as ores and slags, to protect the public and to ensure occupational health and radiation safety. In this context, we discuss three ways of radiation exposure from nodules, including the inhalation or ingestion of nodule fines, the inhalation of radon gas in enclosed spaces and the potential concentration of some radioisotopes during nodule processing. Seen in this light, inappropriate handling of polymetallic nodules poses serious health risks.

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.

Description: Recent studies indicate that mantle plumes, which transfer material and heat from the earth’s interior to its surface, represent multifaceted upwellings. The Tristan-Gough hotspot track (South Atlantic), which formed above a mantle plume, documents spatial geochemical zonation in two distinct sub-tracks since ~70 Ma. The origin and the sudden appearance of two distinct geochemical flavors is enigmatic, but could provide insights into the structural evolution of mantle plumes. Sr–Nd–Pb–Hf isotope data from the Late Cretaceous Rio Grande Rise and adjacent Jean Charcot Seamount Chain (South American Plate), which represent the counterpart of the older Tristan-Gough volcanic track (African Plate), extends the bilateral-zonation to ~100 Ma. Our results support recent numerical models, demonstrating that mantle plumes can split into distinct upper mantle conduits, and provide evidence that these plumelets formed at the plume head-to-plume tail transition. We attribute the plume zonation to sampling the geochemically-graded margin of the African Large Low-Shear-Velocity Province.

Description: Changes in the Atlantic Meridional Overturning Circulation (AMOC) represent a crucial component of Northern Hemisphere climate variability. In modelling studies decadal overturning variability has been attributed to the intensity of deep winter convection in the Labrador Sea. This linkage is challenged by transport observations at sections across the subpolar gyre. Here we report simulations with an eddy-rich ocean model which captures the observed concentration of downwelling in the northeastern Atlantic and the negligible impact of interannual variations in Labrador Sea convection during the last decade. However, the exceptionally cold winters in the Labrador Sea during the first half of the 1990s induced a positive AMOC anomaly of more than 20%, mainly by augmenting the downwelling in the northeastern North Atlantic. The remote effect of excessive Labrador Sea buoyancy forcing is related to rapid spreading of mid-depth density anomalies into the Irminger Sea and their entrainment into the deep boundary current off Greenland.

Description: Global warming may alter the dynamics of infectious diseases by affecting important steps in the transmission of pathogens and parasites. In trematode parasites, the emergence of cercarial stages from their hosts is temperature-dependent, being highest around a thermal optimum. If environmental temperatures exceed this optimum as a consequence of global warming, this may affect cercarial transmission. However, our knowledge of cercarial emergence patterns of species from high temperature environments is currently very limited. Here, we investigated the effect of temperature on the emergence of two common trematode species from an abundant mud snail Pirenella cingulata in the Persian Gulf, the warmest sea on Earth. Infected snails were incubated in the laboratory at 6 temperatures from 10 to 40°C for 3 days. We found an optimal temperature for cercarial emergence of 32.0°C and 33.5°C for Acanthotrema tridactyla and Cyathocotylidae gen. sp., respectively, which are the warmest recorded thermal optima for any aquatic trematode species. Emergence of both species dropped at 40°C, suggesting upper thermal limits to emergence. Overall, Persian Gulf trematodes may be among the most heat-tolerant marine trematode species, indicating a potential for dispersing to regions that will continue to warm in the future.

Description: Knowledge of the characteristics of natural climate variability is vital when assessing the range of plausible future climate trajectories in the next decades to centuries. The reliable detection of climate fluctuations on multidecadal to centennial timescales depends on proxy reconstructions and model simulations, as the instrumental record extends back only a few decades in most parts of the world. Systematic comparisons between model-simulated and proxy-based inferences of natural variability, however, often seem contradictory. Locally, simulated temperature variability is consistently smaller on multidecadal and longer timescales than is indicated by proxy-based reconstructions, implying that climate models or proxy interpretations might have deficiencies. In contrast, at global scales, studies found agreement between simulated and proxy reconstructed temperature variations. Here we review the evidence regarding the scale of natural temperature variability during recent millennia. We identify systematic reconstruction deficiencies that may contribute to differing local and global model–proxy agreement but conclude that they are probably insufficient to resolve such discrepancies. Instead, we argue that regional climate variations persisted for longer timescales than climate models simulating past climate states are able to reproduce. This would imply an underestimation of the regional variability on multidecadal and longer timescales and would bias climate projections and attribution studies. Thus, efforts are needed to improve the simulation of natural variability in climate models accompanied by further refining proxy-based inferences of variability.

Description: As Earth's atmospheric temperatures and human populations increase, more people are becoming vulnerable to natural and human-induced disasters. This is particularly true in Central America, where the growing human population is experiencing climate extremes (droughts and floods), and the region is susceptible to geological hazards, such as earthquakes and volcanic eruptions, and environmental deterioration in many forms (soil erosion, lake eutrophication, heavy metal contamination, etc.). Instrumental and historical data from the region are insufficient to understand and document past hazards, a necessary first step for mitigating future risks. Long, continuous, well-resolved geological records can, however, provide a window into past climate and environmental changes that can be used to better predict future conditions in the region. The Lake Izabal Basin (LIB), in eastern Guatemala, contains the longest known continental records of tectonics, climate, and environmental change in the northern Neotropics. The basin is a pull-apart depression that developed along the North American and Caribbean plate boundary ∼ 12 Myr ago and contains 〉 4 km of sediment. The sedimentological archive in the LIB records the interplay among several Earth System processes. Consequently, exploration of sediments in the basin can provide key information concerning: (1) tectonic deformation and earthquake history along the plate boundary; (2) the timing and causes of volcanism from the Central American Volcanic Arc; and (3) hydroclimatic, ecologic, and geomicrobiological responses to different climate and environmental states. To evaluate the LIB as a potential site for scientific drilling, 65 scientists from 13 countries and 33 institutions met in Antigua, Guatemala, in August 2022 under the auspices of the International Continental Scientific Drilling Program (ICDP) and the US National Science Foundation (NSF). Several working groups developed scientific questions and overarching hypotheses that could be addressed by drilling the LIB and identified optimal coring sites and instrumentation needed to achieve the project goals. The group also discussed logistical challenges and outreach opportunities. The project is not only an outstanding opportunity to improve our scientific understanding of seismotectonic, volcanic, paleoclimatic, paleoecologic, and paleobiologic processes that operate in the tropics of Central America, but it is also an opportunity to improve understanding of multiple geological hazards and communicate that knowledge to help increase the resilience of at-risk Central American communities.

Description: Abyssal seafloor communities cover more than 60% of Earth’s surface. Despite their great size, abyssal plains extend across modest environmental gradients compared to other marine ecosystems. However, little is known about the patterns and processes regulating biodiversity or potentially delimiting biogeographical boundaries at regional scales in the abyss. Improved macroecological understanding of remote abyssal environments is urgent as threats of widespread anthropogenic disturbance grow in the deep ocean. Here, we use a new, basin-scale dataset to show the existence of clear regional zonation in abyssal communities across the 5,000 km span of the Clarion–Clipperton Zone (northeast Pacific), an area targeted for deep-sea mining. We found two pronounced biogeographic provinces, deep and shallow-abyssal, separated by a transition zone between 4,300 and 4,800 m depth. Surprisingly, species richness was maintained across this boundary by phylum-level taxonomic replacements. These regional transitions are probably related to calcium carbonate saturation boundaries as taxa dependent on calcium carbonate structures, such as shelled molluscs, appear restricted to the shallower province. Our results suggest geochemical and climatic forcing on distributions of abyssal populations over large spatial scales and provide a potential paradigm for deep-sea macroecology, opening a new basis for regional-scale biodiversity research and conservation strategies in Earth’s largest biome.

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: 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: Increasing Greenland Ice Sheet–melting is anticipated to impact watermass transformation in the subpolar North Atlantic and ultimately the meridional overturning circulation. Complex ocean and climate models are widely applied to predict magnitude and timing of related impacts under projected future climate. We discuss the role of the ocean mean state, subpolar gyre circulation, mesoscale eddies and atmospheric coupling in shaping the response of the subpolar North Atlantic Ocean to enhanced Greenland runoff. In a suite of eight dedicated 60 to 100-year long model experiments with and without atmospheric coupling, with eddy processes parameterized and explicitly simulated, with regular and significantly enlarged Greenland runoff, we find (1) a major impact by the interactive atmosphere in enabling a compensating temperature feedback, (2) a non-negligible influence by the ocean mean state biased towards greater stability in the coupled simulations, both of which making the Atlantic Merdional Overturning Circulation less susceptible to the freshwater perturbation applied, and (3) a more even spreading of the runoff tracer in the subpolar North Atlantic and enhanced inter-gyre exchange with the subtropics in the strongly eddying simulations. Overall, our experiments demonstrate the important role of mesoscale ocean dynamics and atmosphere feedbacks in projections of the climate system response to enhanced Greenland Ice Sheet–melting and hence underline the necessity to advance scale-aware eddy parameterizations for next-generation climate models.

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: The Palaeocene–Eocene Thermal Maximum (PETM) was a global warming event of 5–6 °C around 56 million years ago caused by input of carbon into the ocean and atmosphere. Hydrothermal venting of greenhouse gases produced in contact aureoles surrounding magmatic intrusions in the North Atlantic Igneous Province have been proposed to play a key role in the PETM carbon-cycle perturbation, but the precise timing, magnitude and climatic impact of such venting remains uncertain. Here we present seismic data and the results of a five-borehole transect sampling the crater of a hydrothermal vent complex in the Northeast Atlantic. Stable carbon isotope stratigraphy and dinoflagellate cyst biostratigraphy reveal a negative carbon isotope excursion coincident with the appearance of the index taxon Apectodinium augustum in the vent crater, firmly tying the infill to the PETM. The shape of the crater and stratified sediments suggests large-scale explosive gas release during the initial phase of vent formation followed by rapid, but largely undisturbed, diatomite-rich infill. Moreover, we show that these vents erupted in very shallow water across the North Atlantic Igneous Province, such that volatile emissions would have entered the atmosphere almost directly without oxidation to CO 2 and at the onset of the PETM.

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: Using population genomics, we reveal the global dispersal history of eelgrass (Zostera marina) from its origin in the Northwest Pacific Ocean, across to the East Pacific, Atlantic and the Mediterranean Sea. Striking differences in genetic diversity among the locations reflect past glaciations and repeated bottlenecks during Z. marina’s worldwide colonization.

Description: Armed conflicts have, in addition to severe impacts on human lives and infrastructure, also impacts on the environment, which needs to be assessed and documented. On September the 26th 2022, unknown perpetrators deliberately ruptured the two gas pipelines Nord Stream 1 and 2 with four coordinated explosions near a major chemical munition dump site near the Danish island of Bornholm in the Baltic Sea. While the massive release of natural gas into atmosphere raised serious concerns concerning the contribution to climate change—this paper assesses the overlooked direct impact of the explosions on the marine ecosystem. Seals and porpoises within a radius of four km would be at high risk of being killed by the shockwave, while temporary impact on hearing would be expected up to 50 km away. As the Baltic Proper population of harbour porpoises ( Phocoena phocoena ) is critically endangered, the loss or serious injury of even a single individual is considered a significant impact on the population. The rupture moreover resulted in the resuspension of 250000 metric tons of heavily contaminated sediment from deep-sea sedimentary basin for over a week, resulting in unacceptable toxicological risks towards fish and other biota in 11 km 3 water in the area for more than a month.

Description: The presence of anti-nutritive compounds like glucosinolates (GSLs) in the rapeseed meal severely restricts its utilization as animal feed. Therefore, reducing the GSL content to 〈 18 µmol/g dry weight in the seeds is a major breeding target. While candidate genes involved in the biosynthesis of GSLs have been described in rapeseed, comprehensive functional analyses are missing. By knocking out the aliphatic GSL biosynthesis genes BnMYB28 and BnCYP79F1 encoding an R2R3 MYB transcription factor and a cytochrome P450 enzyme, respectively, we aimed to reduce the seed GSL content in rapeseed. After expression analyses on single paralogs, we used an ethyl methanesulfonate (EMS) treated population of the inbred winter rapeseed 'Express617' to detect functional mutations in the two gene families. Our results provide the first functional analysis by knock-out for the two GSL biosynthesis genes in winter rapeseed. We demonstrate that independent knock-out mutants of the two genes possessed significantly reduced seed aliphatic GSLs, primarily progoitrin. Compared to the wildtype Express617 control plants (36.3 µmol/g DW), progoitrin levels were decreased by 55.3% and 32.4% in functional mutants of BnMYB28 (16.20 µmol/g DW) and BnCYP79F1 (24.5 µmol/g DW), respectively. Our study provides a strong basis for breeding rapeseed with improved meal quality in the future.

Description: Globalization challenges sustainability by intensifying the ecological and economic impacts of biological invasions. These impacts may be unevenly distributed worldwide, with costs disproportionately incurred by a few regions. We identify economic cost distributions of invasions among origin and recipient countries and continents, and determine socio-economic and biodiversity-related predictors of cost dynamics. Using data filtered from the InvaCost database, which inevitably includes geographic biases in cost reporting, we found that recorded costly invasive alien species have originated from almost all regions, most frequently causing impacts to Europe. In terms of cost magnitude, reported monetary costs predominantly resulted from species with origins in Asia impacting North America. High reported cost linkages (flows) between species' native countries and their invaded countries were related to proxies of shared environments and shared trade history. This pattern can be partly attributed to the legacy of colonial expansion and trade patterns. The characterization of 'sender' and 'receiver' regions of invasive alien species and their associated cost can contribute to more sustainable economies and societies while protecting biodiversity by informing biosecurity planning and the prioritization of control efforts across invasion routes.

Description: Heinrich-type ice-sheet surges are one of the prominent signals of glacial climate variability. They are characterised as abrupt, quasi-periodic episodes of ice-sheet instabilities during which large numbers of icebergs are released from the Laurentide ice sheet. The mechanisms controlling the timing and occurrence of Heinrich-type ice-sheet surges remain poorly constrained to this day. Here, we use a coupled ice sheet–solid Earth model to identify and quantify the importance of boundary forcing for the surge cycle length of Heinrich-type ice-sheet surges for two prominent ice streams of the Laurentide ice sheet – the land-terminating Mackenzie ice stream and the marine-terminating Hudson ice stream. Both ice streams show responses of similar magnitude to surface mass balance and geothermal heat flux perturbations, but Mackenzie ice stream is more sensitive to ice surface temperature perturbations, a fact likely caused by the warmer climate in this region. Ocean and sea-level forcing as well as different frequencies of the same forcing have a negligible effect on the surge cycle length. The simulations also highlight the fact that only a certain parameter space exists under which ice-sheet oscillations can be maintained. Transitioning from an oscillatory state to a persistent ice streaming state can result in an ice volume loss of up to 30 % for the respective ice stream drainage basin under otherwise constant climate conditions. We show that Mackenzie ice stream is susceptible to undergoing such a transition in response to all tested positive climate perturbations. This underlines the potential of the Mackenzie region to have contributed to prominent abrupt climate change events of the last deglaciation.

Description: Subpolar overturning in the North Atlantic Ocean shows substantial seasonality, with a maximum in late spring, a minimum in early winter, and a total range of about 9 Sv, according to observations from the OSNAP array between 2014 and 2020. Understanding the variability of the Atlantic Meridional Overturning Circulation is essential for better predictions of our changing climate. Here we present an updated time series (August 2014 to June 2020) from the Overturning in the Subpolar North Atlantic Program. The 6-year time series allows us to observe the seasonality of the subpolar overturning and meridional heat and freshwater transports. The overturning peaks in late spring and reaches a minimum in early winter, with a peak-to-trough range of 9.0 Sv. The overturning seasonal timing can be explained by winter transformation and the export of dense water, modulated by a seasonally varying Ekman transport. Furthermore, over 55% of the total meridional freshwater transport variability can be explained by its seasonality, largely owing to overturning dynamics. Our results provide the first observational analysis of seasonality in the subpolar North Atlantic overturning and highlight its important contribution to the total overturning variability observed to date.

Description: Understanding the relationship between surface marine ecosystems and the export of carbon to depth by sinking organic particles is key to representing the effect of ecosystem dynamics and diversity, and their evolution under multiple stressors, on the carbon cycle and climate in models. Recent observational technologies have greatly increased the amount of data available, both for the abundance of diverse plankton groups and for the concentration and properties of particulate organic carbon in the ocean interior. Here we use synthetic model data to test the potential of using machine learning (ML) to reproduce concentrations of particulate organic carbon within the ocean interior based on surface ecosystem and environmental data. We test two machine learning methods that differ in their approaches to data-fitting, the random forest and XGBoost methods. The synthetic data are sampled from the PlankTOM12 global biogeochemical model using the time and coordinates of existing observations. We test 27 different combinations of possible drivers to reconstruct small (POCS) and large (POCL) particulate organic carbon concentrations. We show that ML can successfully be used to reproduce modelled particulate organic carbon over most of the ocean based on ecosystem and modelled environmental drivers. XGBoost showed better results compared to random forest thanks to its gradient boosting trees' architecture. The inclusion of plankton functional types (PFTs) in driver sets improved the accuracy of the model reconstruction by 58 % on average for POCS and by 22 % for POCL. Results were less robust over the equatorial Pacific and some parts of the high latitudes. For POCS reconstruction, the most important drivers were the depth level, temperature, microzooplankton and PO4, while for POCL it was the depth level, temperature, mixed-layer depth, microzooplankton, phaeocystis, PO4 and chlorophyll a averaged over the mixed-layer depth. These results suggest that it will be possible to identify linkages between surface environmental and ecosystem structure and particulate organic carbon distribution within the ocean interior using real observations and to use this knowledge to improve both our understanding of ecosystem dynamics and of their functional representation within models.

Description: The southern African climate is strongly impacted by climate change. Precipitation is a key variable in this region, as it is linked to agriculture and water supply. Simulations with a regional atmospheric model over the past decades and the 21st century display a decrease in the past precipitation over some coastal areas of South Africa and an increase over the rest of southern Africa. However, precipitation is projected to decrease over the whole southern part of the domain in the future. This study shows that the Agulhas Current system, including the current and the leakage, which surrounds the continent in the east and south, impacts this precipitation trend. A reduction in the strength of the Agulhas Current is linked to a reduction in precipitation along the southeast coast. The Agulhas leakage, the part of the Agulhas Current that leaves the system and flows into the South Atlantic, impacts winter precipitation in the southwest of the continent. A more intense Agulhas leakage is linked to a reduction in precipitation in this region.

Description: The planktonic realm from bacteria to zooplankton provides the baseline for pelagic aquatic food webs. However, multiple trophic levels are seldomly included in time series studies, hampering a holistic understanding of the influence of seasonal dynamics and species interactions on food web structure and biogeochemical cycles. Here, we investigated plankton community composition, focusing on bacterio-, phyto- and large mesozooplankton, and how biotic and abiotic factors correlate at the Linnaeus Microbial Observatory (LMO) station in the Baltic Sea from 2011 to 2018. Plankton communities structures showed pronounced dynamic shifts with recurring patterns. Summarizing the parts of the planktonic microbial food web studied here to total carbon, a picture emerges with phytoplankton consistently contributing 〉 39% while bacterio- and large mesozooplankton contributed ~ 30% and ~ 7%, respectively, during summer. Cyanophyceae, Actinobacteria, Bacteroidetes, and Proteobacteria were important groups among the prokaryotes. Importantly, Dinophyceae, and not Bacillariophyceae, dominated the autotrophic spring bloom whereas Litostomatea (ciliates) and Appendicularia contributed significantly to the consumer entities together with the more traditionally observed mesozooplankton, Copepoda and Cladocera. Our findings of seasonality in both plankton composition and carbon stocks emphasize the importance of time series analyses of food web structure for characterizing the regulation of biogeochemical cycles and appropriately constraining ecosystem models.

Description: Nutrient availability limits phytoplankton growth throughout much of the global ocean. Here we synthesize available experimental data to identify three dominant nutrient limitation regimes: nitrogen is limiting in the stratified subtropical gyres and in the summertime Arctic Ocean, iron is most commonly limiting in upwelling regions, and both nutrients are frequently co-limiting in regions in between the nitrogen and iron limited systems. Manganese can be co-limiting with iron in parts of the Southern Ocean, whilst phosphate and cobalt can be co-/serially limiting in some settings. Overall, an analysis of experimental responses showed that phytoplankton net growth can be significantly enhanced through increasing the number of different nutrients supplied, regardless of latitude, temperature, or trophic status, implying surface seawaters are often approaching nutrient co-limitation. Assessments of nutrient deficiency based on seawater nutrient concentrations and nutrient stress diagnosed via molecular biomarkers showed good agreement with experimentally-assessed nutrient limitation, validating conceptual and theoretical links between nutrient stoichiometry and microbial ecophysiology.

Description: In the hadal zone of the ocean (6–11 km), the characteristics of sinking marine snow particles and their attached microbial communities remain elusive, despite their potential importance for benthic life thriving at extreme pressures (60–110 MPa). Here, we used simulation experiments to explore how increasing pressure levels modify the microbial degradation, organic matter composition, and microbiome of sinking diatom aggregates. Individual aggregates were incubated in rotating tanks in which pressure was incrementally increased to simulate a descent from surface to hadal depth within 20 days. Incubations at atmospheric pressure served as controls. With increasing pressure, microbial respiration and diatom degradation decreased gradually and ceased completely at 60 MPa. Dissolved organic carbon leaked substantially from the aggregates at ≥40 MPa, while diatom lipid and pigment contents decreased moderately. Bacterial abundance remained stable at 〉40 MPa, but bacterial community composition changed significantly at 60–100 MPa. Thus, pressure exposure reduces microbial degradation and transforms both organic matter composition and microbiomes of sinking particles, which may seed hadal sediments with relatively fresh particulate organic matter and putative pressure-tolerant microbes.

Description: Marine diazotrophs convert dinitrogen (N-2) gas into bioavailable nitrogen (N), supporting life in the global ocean. In 2012, the first version of the global oceanic diazotroph database (version 1) was published. Here, we present an updated version of the database (version 2), significantly increasing the number of in situ diazotrophic measurements from 13 565 to 55 286. Data points for N-2 fixation rates, diazotrophic cell abundance, and nifH gene copy abundance have increased by 184 %, 86 %, and 809 %, respectively. Version 2 includes two new data sheets for the nifH gene copy abundance of non-cyanobacterial diazotrophs and cell-specific N2 fixation rates. The measurements of N-2 fixation rates approximately follow a log-normal distribution in both version 1 and version 2. However, version 2 considerably extends both the left and right tails of the distribution. Consequently, when estimating global oceanic N-2 fixation rates using the geometric means of different ocean basins, version 1 and version 2 yield similar rates (43-57 versus 45-63 TgNyr (-1); ranges based on one geometric standard error). In contrast, when using arithmetic means, version 2 suggests a significantly higher rate of 223 +/- 30 TgNyr (-1) (mean +/- standard error; same hereafter) compared to version 1 (74 +/- 7 TgNyr (-1)). Specifically, substantial rate increases are estimated for the South Pacific Ocean (88 +/- 23 versus 20 +/- 2 TgNyr 1), primarily driven by measurements in the southwestern subtropics, and for the North Atlantic Ocean (40 +/- 9 versus 10 +/- 2 TgNyr (-1)). Moreover, version 2 estimates the N-2 fixation rate in the Indian Ocean to be 35 +/- 14 TgNyr (-1), which could not be estimated using version 1 due to limited data availability. Furthermore, a comparison of N-2 fixation rates obtained through different measurement methods at the same months, locations, and depths reveals that the conventional N-15(2) bubble method yields lower rates in 69% cases compared to the new N-15(2) dissolution method. This updated version of the database can facilitate future studies in marine ecology and biogeochemistry. The database is stored at the Figshare repository (https://doi.org/10.6084/m9.figshare.21677687; Shao et al., 2022).

Description: Eutrophication is accelerating the recent expansion of oxygen-depleted coastal marine environments. Several bolivinid foraminifera are abundant in these oxygen-depleted settings, and take up nitrate through the pores in their shells for denitrification. This makes their pore density a possible nitrate proxy. This study documents three aspects related to the porosity of bolivinids. 1. A new automated image analysis technique to determine the number of pores in bolivinids is tested. 2. The pore patterns of Bolivina spissa from five different ocean settings are analysed. The relationship between porosity, pore density and mean pore size significantly differs between the studied locations. Their porosity is mainly controlled by the size of the pores at the Gulf of Guayaquil (Peru), but by the number of pores at other studied locations. This might be related to the presence of a different cryptic Bolivina species in the Gulf of Guayaquil. 3. The pore densities of closely related bolivinids in core-top samples are calibrated as a bottom-water nitrate proxy. Bolivina spissa and Bolivina subadvena showed the same correlation between pore density and bottom-water nitrate concentrations, while the pore density of Bolivina argentea and Bolivina subadvena accumeata is much higher.

Description: Stable water isotopes in polar ice cores are widely used to reconstruct past temperature variations over several orbital climatic cycles. One way to calibrate the isotope–temperature relationship is to apply the present-day spatial relationship as a surrogate for the temporal one. However, this method leads to large uncertainties because several factors like the sea surface conditions or the origin and transport of water vapor influence the isotope–temperature temporal slope. In this study, we investigate how the sea surface temperature (SST), the sea ice extent, and the strength of the Atlantic Meridional Overturning Circulation (AMOC) affect these temporal slopes in Greenland and Antarctica for Last Glacial Maximum (LGM, ∼ 21 000 years ago) to preindustrial climate change. For that, we use the isotope-enabled atmosphere climate model ECHAM6-wiso, forced with a set of sea surface boundary condition datasets based on reconstructions (e.g., GLOMAP) or MIROC 4m simulation outputs. We found that the isotope–temperature temporal slopes in East Antarctic coastal areas are mainly controlled by the sea ice extent, while the sea surface temperature cooling affects the temporal slope values inland more. On the other hand, ECHAM6-wiso simulates the impact of sea ice extent on the EPICA Dome C (EDC) and Vostok sites through the contribution of water vapor from lower latitudes. Effects of sea surface boundary condition changes on modeled isotope–temperature temporal slopes are variable in West Antarctica. This is partly due to the transport of water vapor from the Southern Ocean to this area that can dampen the influence of local temperature on the changes in the isotopic composition of precipitation and snow. In the Greenland area, the isotope–temperature temporal slopes are influenced by the sea surface temperatures near the coasts of the continent. The greater the LGM cooling off the coast of southeastern Greenland, the greater the transport of water vapor from the North Atlantic, and the larger the temporal slopes. The presence or absence of sea ice very near the coast has a large influence in Baffin Bay and the Greenland Sea and influences the slopes at some inland ice core stations. The extent of the sea ice far south slightly influences the temporal slopes in Greenland through the transport of more depleted water vapor from lower latitudes to this area. The seasonal variations of sea ice distribution, especially its retreat in summer, influence the isotopic composition of the water vapor in this region and the modeled isotope–temperature temporal slopes in the eastern part of Greenland. A stronger LGM AMOC decreases LGM-to-preindustrial isotopic anomalies in precipitation in Greenland, degrading the isotopic model–data agreement. The AMOC strength modifies the temporal slopes over inner Greenland slightly and by a little on the coasts along the Greenland Sea where the changes in surface temperature and sea ice distribution due to the AMOC strength mainly occur.

Description: An international, multidisciplinary research group is proposing the “NICA-BRIDGE” drilling project, within the framework of the International Continental Scientific Drilling Program (ICDP). The project goal is to conduct scientific drilling in Lake Nicaragua and Lake Managua (Nicaragua, Central America) to obtain long lacustrine sediment records to (a) extend the neotropical paleoclimate record back to the Pliocene, making it one of the longest continental tropical climate archives in the world, and to (b) provide geological data on the long-term complex interplay among tectonics, volcanism, sea-level dynamics, climate change, and biosphere. The lakes are the two largest in Central America, and they are located in a trench-parallel half graben that hosts the volcanic front, which developed during or prior to the Pliocene, as a consequence of subduction-related tectonic activity. The lakes are uniquely suited for multidisciplinary scientific investigation as their long, con- tinuous sediment records (several Myr) will facilitate the study of (1) terrestrial and marine basin development at the southern Central American margin, (2) alternating lacustrine and marine environments in response to tec- tonic and climatic changes, (3) the longest record of tropical climate proxies, (4) the evolution of (and transition between) the Miocene to Pliocene/Pleistocene and Pleistocene to present volcanic arcs, which were separated by slab rollback, (5) the significance of the lakes as hot spots for endemism, and (6) the Great American Biotic Interchange at this strategic location, i.e., the N–S and reverse migration of fauna after the land bridge between the Americas was established. The planned ICDP project offers an opportunity to explore these topics through continent-based seismolog- ical, volcanological, paleoclimatological, paleoecological, and paleoenvironmental studies, combined with an International Ocean Discovery Program (IODP) drill project to explore its oceanic continuation. In preparation of this drilling project, an ICDP workshop was held in Montelimar, Nicaragua, on 2–5 March 2020 to develop drilling strategies and refine scientific questions, objectives, and hypotheses. The workshop was organized and hosted by the principal investigators and the Instituto Nicaragüense de Estudios Territoriales (INETER), with funding from the ICDP. Forty-five researchers from 12 countries participated in the workshop, including representatives from ICDP. During the workshop, previous research data on the study lakes, including new recent surveys, were reviewed, and a three-phase strategy for the proposed research was developed. The aim of Phase 0 is to complement the pre-site surveys where we identified the need for further data. In Phase I, with ICDP support, we will obtain sediment cores ∼ 100 m long, which will allow us to investigate many of the scientific questions. Based on the data from those drill cores, coring locations will be identified for a future Phase II, which we envisage as a combined ICDP/IODP project to collect deep drill cores in the lakes and the offshore Sandino Basin in order to extend Phase I results to much deeper time. The Sandino Basin is the oceanic continuation of the depression in which the studied lakes are located, and complementary marine drilling will improve the understanding of the evolution of this complex margin.

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: Ocean ecosystems are at the forefront of the climate and biodiversity crises, yet we lack a unified approach to assess their state and inform sustainable policies. This blueprint is designed around research capabilities and cross-sectoral partnerships. We highlight priorities including integrating basin-scale observation, modelling and genomic approaches to understand Atlantic oceanography and ecosystem connectivity; improving ecosystem mapping; identifying potential tipping points in deep and open ocean ecosystems; understanding compound impacts of multiple stressors including warming, acidification and deoxygenation; enhancing spatial and temporal management and protection. We argue that these goals are best achieved through partnerships with policy-makers and community stakeholders, and promoting research groups from the South Atlantic through investment and engagement. Given the high costs of such research (€800k to €1.7M per expedition and €30–40M for a basin-scale programme), international cooperation and funding are integral to supporting science-led policies to conserve ocean ecosystems that transcend jurisdictional borders.