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: In climate reanalyses for multi-decadal or longer scales with coupled atmosphere-ocean General Circulation models (CGCMs) it can be assumed that the growth of prediction errors arises chiefly from imprecisely known model parameters, which have a nonlinear relationship with the climate observations (paleoclimate proxies). Also, high-resolution CGCMs for climate analysis are extremely expensive to run, which constrains the applicability of assimilation schemes. In a model framework where we assume that model dynamic parameters account for (nearly) all forecast errors at observation times, we compare two computationally efficient iterative schemes for approximate nonlinear model parameter estimation and joint flux estimation (taking the specific shape of freshwater from melting in the Greenland ice sheet), and its physically consistent state. First, a trivial adaptation of the strong constraint incremental 4D-Var formulation leads to what we refer to as the parameter space iterative extended Kalman smoother (pIKS); a Gauss-Newton scheme. Second, a so-called parameter space fractional Kalman smoother (pFKS) is an alternative controlled-step line search, which can potentially be a more stable approach. While these iterative schemes have been used in data assimilation, we revisit them together within the context of parameter estimation in climate reanalysis, as compared to the more general 4D-Var formulation. Then, the two schemes are evaluated in numerical experiments with a simple 1D energy balance model (Ebm1D) and with a fully-coupled Community Earth System Model (CESM v1.2). Firstly, with Ebm1D the pFKS obtains a cost function similar to the adjoint method with highly reduced computational cost, while an ensemble transform Kalman filter with an m = 60 ensemble size (ETKF60) behaves slightly worse. The pIKS behaves worse than the ETKF60, but an ETKF10 (m = 10) is even worst. Accordingly, with CESM we evaluate the pKFS and the ETKF60 along with an ETKF with Gaussian Anamorphosis (ETKF-GA60). From all the options, the pFKS has the lowest cost function and seems the favored overall option under heavy computational restrictions, but the ETKF obtains better estimates of the flux term.

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: Nitrogen fixers, or diazotrophs, play a key role in the carbon and nitrogen cycle of the world oceans, but the controlling mechanisms are not comprehensively understood yet. The present study compares two paradigms on the ecological niche of diazotrophs in an Earth System Model (ESM). In our standard model configuration, which is representative for most of the state-of-the-art pelagic ecosystem models, diazotrophs take advantage of zooplankton featuring a lower food preference for diazotrophs than for ordinary phytoplankton. We compare this paradigm with the idea that diazotrophs are more competitive under oligotrophic conditions, characterized by low (dissolved, particulate, organic and inorganic) phosphorous availability. Both paradigms are supported by observational evidence and lead to a similar good agreement to the most recent and advanced observation-based nitrogen fixation estimate in our ESM framework. Further, we illustrate that the similarity between the two paradigms breaks in a RCP 8.5 anthropogenic emission scenario. We conclude that a more advanced understanding of the ecological niche of diazotrophs is mandatory for assessing the cycling of essential nutrients, especially under changing environmental conditions. Our results call for more in-situ measurements of cyanobacteria biomass if major controls of nitrogen fixation in the oceans are to be dissected.

Description: Nitrogen fixers, or diazotrophs, play a key role in the carbon and nitrogen cycle of the world oceans, but the controlling mechanisms are not comprehensively understood yet. The present study compares two paradigms on the ecological niche of diazotrophs in an Earth System Model (ESM). In our standard model configuration, which is representative for most of the state-of-the-art pelagic ecosystem models, diazotrophs take advantage of zooplankton featuring a lower food preference for diazotrophs than for ordinary phytoplankton. We compare this paradigm with the idea that diazotrophs are more competitive under oligotrophic conditions, characterized by low (dissolved, particulate, organic and inorganic) phosphorous availability. Both paradigms are supported by observational evidence and lead to a similar good agreement to the most recent and advanced observation-based nitrogen fixation estimate in our ESM framework. Further, we illustrate that the similarity between the two paradigms breaks in a RCP 8.5 anthropogenic emission scenario. We conclude that a more advanced understanding of the ecological niche of diazotrophs is mandatory for assessing the cycling of essential nutrients, especially under changing environmental conditions. Our results call for more in-situ measurements of cyanobacteria biomass if major controls of nitrogen fixation in the oceans are to be dissected.

Description: Although marine controlled source electromagnetic (CSEM) methods are effective for investigating offshore freshened groundwater (OFG) systems, interpreting the spatial extent and salinity of OFG remains challenging. Integrating CSEM resistivity models with information on sub-surface properties, such as host-rock porosity, allows for estimates of pore-water salinity. However, deterministic inversion approaches pose challenges in quantitatively analyzing these estimates as they provide only one best-fit model with no associated estimate of model parameter uncertainty. To address this limitation, we employ a trans-dimensional Markov-Chain Monte-Carlo inversion on marine CSEM data, under the assumption of horizontal stratification, collected from the Canterbury Bight, New Zealand. We integrate the resulting posterior distributions of electrical resistivity with borehole and seismic reflection data to quantify pore-water salinity with uncertainty estimates. The results reveal a low-salinity groundwater body in the center of the survey area at varying depths, hosted by consecutive silty- and fine-sand layers approximately 20 to 60 km from the coast. These observations support the previous study’s results obtained through deterministic 2-D inversion and suggest freshening of the OFG body closer to the shore within a permeable, coarse-sand layer 40 to 150 m beneath the seafloor. This implies a potential active connection between the OFG body and the terrestrial groundwater system. We demonstrate how the Bayesian approach constrains the uncertainties in resistivity models and subsequently in pore-water salinity estimates. Our findings highlight the potential of Bayesian inversions in enhancing our understanding of OFG systems, providing crucial boundary conditions for hydrogeological modeling and sustainable water resource development.

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: Nonmammalian tooth-bearing vertebrates usually replace their teeth throughout life. Much about how a replacement pattern is generated has been learned from zebrafish. However, to understand general mechanisms of tooth replacement, advantage can be taken from studying other, “nonmodel” species. We have mapped the patterns of tooth replacement in widely divergent aquatic osteichthyans using 2D charts, in which one axis is time, the other linear spacing along the tooth row. New teeth that are generated simultaneously are considered part of the same odontogenic wave. Using this approach, it appears that a similar, general pattern underlies very distinctive dentitions in distantly related species. A simple shift in spacing of odontogenic waves, or in distance between subsequent tooth positions along a row (or both), can produce dramatically different dentitions between life stages within a species, or between closely related species. Examples will be presented from salmonids, cyprinids, and cichlids. Our observations suggest that lines linking subsequent positions may have more biological significance than replacement waves (usually linking alternate positions), often used to explain the generation of patterns. The presence of a general pattern raises questions about common control mechanisms. There is now increasing evidence, at least for the zebrafish, to support a role for stem cells in continuous tooth renewal and control of replacement patterns.