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: Spatial predictions of total organic carbon (TOC) concentrations and stocks are crucial for understanding marine sediments’ role as a significant carbon sink in the global carbon cycle. In this study, we present a geospatial prediction of TOC concentrations and stocks at a 5 x 5 arc minute grid scale, using a deep learning model — a novel machine learning approach based on a new compilation of over 22,000 global TOC measurements and a new set of predictors, such as seafloor lithologies, grain size distribution, and an alpha-chlorophyll satellite data. In our study, we compared the predictions and discuss the limitations from various machine learning methods. Our findings reveal that the neural network approach outperforms methods such as k Nearest Neighbors and random forests, which tend to overfit to the training data, especially in highly heterogeneous and complex geological settings. We provide estimates of mean TOC concentrations and total carbon stock in both continental shelves and deep sea settings across various marine regions and oceans. Our model suggests that the upper 10 cm of oceanic sediments harbors approximately 171 Pg of TOC stock and has a mean TOC concentration of 0.68 %. Furthermore, we introduce a standardized methodology for quantifying predictive uncertainty using Monte Carlo dropout and present a map of information gain, that measures the expected increase in model knowledge achieved through in-situ sampling at specific locations which is pivotal for sampling strategy planning.

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: 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: 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: Current global warming results in rising sea-water temperatures, and the loss of sea ice in arctic and subarctic oceans impacts the community composition of primary producers with cascading effects on the food web and potentially on carbon export rates. This study analyzes metagenomic shotgun and diatom rbcL amplicon-sequencing data from sedimentary ancient DNA (sedaDNA) of the subarctic western Bering Sea that records phyto- and zooplankton community changes over the last glacial–interglacial cycle, including the last interglacial period (Eemian). Our data show that interglacial and glacial plankton communities differ, with distinct Eemian and Holocene plankton communities. The generally warm Holocene period is dominated by pico-sized cyanobacteria and bacteria-feeding heterotrophic protists, while the Eemian period is dominated by eukaryotic pico-sized chlorophytes and Triparmaceae. In contrast, the glacial period is characterized by micro-sized phototrophic protists, including sea-ice associated diatoms in the family Bacillariaceae and co-occurring diatom-feeding crustaceous zooplankton. Our deep-time record of plankton community changes reveals a long-term decrease in phytoplankton cell size coeval with increasing temperatures, and resembling community changes in the currently warming Bering Sea. The phytoplankton community in the warmer-than-present Eemian period is distinct from modern communities and limits the use of the Eemian as an analog for future climate scenarios. However, under enhanced future warming, the expected shift towards the dominance of small-sized phytoplankton and heterotrophic protists might result in an increased productivity, whereas the community’s potential of carbon export will be decreased, thereby weakening the subarctic Bering Sea’s function as an effective carbon sink.

Description: Fjord-like estuaries are hotspots of biogeochemical cycling due to steep physicochemical gradients. The spatiotemporal distribution of nitrous oxide (N2O) within many of these systems is poorly described, especially in the southern hemisphere. The goal of this study is to describe the spatiotemporal distribution of N2O within a southern hemisphere fjord-like estuary, describe the main environmental drivers of this distribution, the air/sea flux of N2O, and the main drivers of N2O production. Cruises were undertaken in Macquarie Harbour, Tasmania to capture N2O concentrations and water column physicochemical profiles in winter (July 2022), spring (October 2022), summer (February 2023), and autumn (April 2023). N2O samples were collected at one depth at system end members, and at 5 depths at 4 stations within the harbour. Results indicate that N2O is consistently supersaturated (reaching 170 % saturation) below the system’s freshwater lens where oxygen concentrations are often hypoxic, but infrequently anoxic. In the surface lens, levels of N2O saturation vary with estimated river flow and with proximity to the system’s main freshwater endmember. The linear relationship between AOU and ΔN2O saturation indicates that nitrification is the process generating N2O in the system. When river flow was high (July and October 2022), surface water N2O was undersaturated (as low as 70 %) throughout most of the harbour. When river flow was low (February and April 2023) N2O was observed to be supersaturated at most stations. Calculated air/sea fluxes of N2O indicated that the system is generally a source of N2O to the atmosphere under weak river flow conditions and a sink during strong river flow conditions. The diapycnal flux was a minor contributor to surface water N2O concentrations, and subhalocline N2O is intercepted by the riverine surface lens and transported out of the system to the ocean during strong river flow conditions. In a changing climate, Western Tasmania is expected to receive higher winter rainfall and lower summer rainfall which may augment the source and sink dynamics of this system by enhancing the summer / autumn efflux of N2O to the atmosphere. This study is the first to report observations of N2O distribution, generation processes, and estimated diapycnal / surface N2O fluxes from this system.