ALBERT

Description: The present study describes Rossby wave packet (RWP) properties in the upper-troposphere and lower-stratosphere (UTLS) with the use of Global Navigation Satellite System radio occultation (GNSS-RO) measurements. This global study covering both hemisphere's extratropics is the first to tackle medium and synoptic-scale waves with GNSS-RO. We use one decade of GNSS-RO temperature and pressure data from the CHAMP, COSMIC, GRACE, Metop-A, Metop-B, SAC-C and TerraSAR-X missions; combining them into one gridded dataset for the years 2007–2016. Our approach to extract RWP anomalies and their envelope uses Fourier and Hilbert transforms over longitude without pre- or post-processing the data. Our study is purely based on observations, only using ERA-Interim winds to provide information about the background wind regimes. The RWP structures that we obtain in the UTLS agree well with theory and earlier studies, in terms of coherent phase/group propagation, zonal scale and distribution over latitudes. Furthermore, we show that RWP pressure anomalies maximize around the tropopause, while RWP temperature anomalies maximize right above tropopause height with a contrasting minimum right below. RWP activity follows the zonal-mean tropopause during all seasons. RWP anomalies in the lower stratosphere are dynamically coupled to the upper troposphere. They are part of the same system with a quasi-barotropic structure across the UTLS. RWP activity often reaches up to 20 km height and occasionally higher, defying the Charney–Drazin criterion. We note enhanced amplitude and upward propagation of RWP activity during sudden stratospheric warmings. We provide observational support for improvements in RWP diagnostics and wave trend analysis in models and reanalyses. Wave quantities follow the tropopause, and diagnosing them on fixed pressure levels (which the tropopause does not follow) can lead to aliasing. Our novel approach analysing GNSS-RO pressure anomalies provides wave signals with better continuity and coherence across the UTLS and the stratosphere, compared to temperature anomalies. Thus, RWP vertical propagation is much easier to analyse with pressure data.

Description: The eastern boundary region of the southeastern Pacific Ocean hosts one of the world's most dynamic and productive upwelling systems with an associated oxygen minimum zone (OMZ). The variability in downward export fluxes in this region, with strongly varying surface productivity, upwelling intensities and water column oxygen content, is however poorly understood. Thorium-234 (234Th) is a powerful tracer to study the dynamics of export fluxes of carbon and other elements, yet intense advection and diffusion in nearshore environments impact the assessment of depth-integrated 234Th fluxes when not properly evaluated. Here we use vessel-mounted acoustic Doppler current profiler (VmADCP) current velocities, satellite wind speed and in situ microstructure measurements to determine the magnitude of advective and diffusive fluxes over the entire 234Th flux budget at 25 stations from 11 to 16∘ S in the Peruvian OMZ. Contrary to findings along the GEOTRACES P16 eastern section, our results showed that weak surface wind speed during our cruises induced low upwelling rates and minimal upwelled 234Th fluxes, whereas vertical diffusive 234Th fluxes were important only at a few shallow shelf stations. Horizontal advective and diffusive 234Th fluxes were negligible because of small alongshore 234Th gradients. Our data indicated a poor correlation between seawater 238U activity and salinity. Assuming a linear relationship between the two would lead to significant underestimations of the total 234Th flux by up to 40 % in our study. Proper evaluation of both physical transport and variability in 238U activity is thus crucial in coastal 234Th flux studies. Finally, we showed large temporal variations on 234Th residence times across the Peruvian upwelling zone and cautioned future carbon export studies to take these temporal variabilities into consideration while evaluating carbon export efficiency.

Description: Rapid monsoon changes since the last deglaciation remain poorly constrained due to the scarcity of geological archives. Here we present a high-resolution scanning X-ray fluorescence (XRF) analysis of a 13.5 m terrace succession on the western Chinese Loess Plateau (CLP) to infer rapid monsoon changes since the last deglaciation. Our results indicate that Rb∕Sr and Zr∕Rb are sensitive indicators of chemical weathering and wind sorting, respectively, which are further linked to the strength of the East Asian summer monsoon (EASM) and the East Asian winter monsoon (EAWM). During the last deglaciation, two cold intervals of the Heinrich event 1 and Younger Dryas were characterized by intensified winter monsoon and weakened summer monsoon. The EAWM gradually weakened at the beginning of the Holocene, while the EASM remained steady till 9.9 ka and then grew stronger. Both the EASM and EAWM intensities were relatively weak during the Middle Holocene, indicating a mid-Holocene climatic optimum. Rb∕Sr and Zr∕Rb exhibit an antiphase relationship between the summer and winter monsoon changes on a centennial timescale during 16–1 ka. Comparison of these monsoon changes with solar activity and North Atlantic cooling events reveals that both factors can lead to abrupt changes on a centennial timescale in the Early Holocene. During the Late Holocene, North Atlantic cooling became the major forcing of centennial monsoon events.

Description: We developed a coupling scheme for the Community Earth System Model version 1.2 (CESM1.2) and the Model of Early Diagenesis in the Upper Sediment of Adjustable complexity (MEDUSA), and explored the effects of the coupling on solid components in the upper sediment and on bottom seawater chemistry by comparing the coupled model's behaviour with that of the uncoupled CESM having a simplified treatment of sediment processes. CESM is a fully coupled atmosphere–ocean–sea-ice–land model and its ocean component (the Parallel Ocean Program version 2; POP2) includes a biogeochemical component (the Biogeochemical Elemental Cycling model; BEC). MEDUSA was coupled to POP2 in an offline manner so that each of the models ran separately and sequentially with regular exchanges of necessary boundary condition fields. This development was done with the ambitious aim of a future application for long-term (spanning a full glacial cycle; i.e. ∼105 years) climate simulations with a state-of-the-art comprehensive climate model including the carbon cycle, and was motivated by the fact that until now such simulations have been done only with less-complex climate models. We found that the sediment–model coupling already had non-negligible immediate advantages for ocean biogeochemistry in millennial-timescale simulations. First, the MEDUSA-coupled CESM outperformed the uncoupled CESM in reproducing an observation-based global distribution of sediment properties, especially for organic carbon and opal. Thus, the coupled model is expected to act as a better “bridge” between climate dynamics and sedimentary data, which will provide another measure of model performance. Second, in our experiments, the MEDUSA-coupled model and the uncoupled model had a difference of 0.2 ‰ or larger in terms of δ13C of bottom water over large areas, which implied a potentially significant model uncertainty for bottom seawater chemical composition due to a different way of sediment treatment. For example, an ocean model that does not treat sedimentary processes depending on the chemical composition of the ambient water can overestimate the amount of remineralization of organic matter in the upper sediment in an anoxic environment, which would lead to lighter δ13C values in the bottom water. Such a model uncertainty would be a fundamental issue for paleo model–data comparison often relying on data derived from benthic foraminifera.

Description: The Sea Surface Microlayer (SML) is known to be enriched by trace metals relative to the underlying water and harbor diverse microbial communities (i.e., neuston). However, the processes linking metals and biota in the SML are not yet fully understood. The metal (Cd, Co, Cu, Fe, Ni, Mo, V, Zn and Pb) concentrations in aerosol samples in the SML (dissolved and total fractions) and in subsurface waters (SSWs; dissolved fraction at ∼1 m depth) from the western Mediterranean Sea were analyzed in this study during a cruise in May–June 2017. The composition and abundance of the bacterial community in the SML and SSW, the primary production, and Chl a in the SSW were measured simultaneously at all stations during the cruise. Residence times in the SML of metals derived from aerosol depositions were highly variable and ranged from minutes for Fe (3.6±6.0 min) to a few hours for Cu (5.8±6.2 h). Concentrations of most of the dissolved metals in both the SML and SSW were positively correlated with the salinity gradient and showed the characteristic eastward increase in the surface waters of the Mediterranean Sea (MS). In contrast, the total fraction of some reactive metals in the SML (i.e., Cu, Fe, Pb and Zn) showed a negative correlation with salinity and a positive correlation with microbial abundance, which might be associated with microbial uptake. Our results show a strong negative correlation between the dissolved and total Ni concentration and heterotrophic bacterial abundance in the SML and SSW, but we cannot ascertain whether this correlation reflects a toxicity effect or is the result of some other process.

Description: Proxy records represent an invaluable source of information for reconstructing past climatic variations, but they are associated with considerable uncertainties. For a systematic quantification of these reconstruction errors, however, knowledge is required not only of their individual sources but also of their auto-correlation structure as this determines the timescale dependence of their magnitude, an issue that has been often ignored until now. Here a spectral approach to uncertainty analysis is provided for paleoclimate reconstructions obtained from single sediment proxy records. The formulation in the spectral domain rather than the time domain allows for an explicit demonstration and quantification of the timescale dependence that is inherent in any proxy-based reconstruction uncertainty. This study is published in two parts. In this first part, the theoretical concept is presented, and analytic expressions are derived for the power spectral density of the reconstruction error of sediment proxy records. The underlying model takes into account the spectral structure of the climate signal, seasonal and orbital variations, bioturbation, sampling of a finite number of signal carriers, and uncorrelated measurement noise, and it includes the effects of spectral aliasing and leakage. The uncertainty estimation method, based upon this model, is illustrated by simple examples. In the second part of this study, published separately, the method is implemented in an application-oriented context, and more detailed examples are presented.

Description: Enhanced summer insolation during the early and mid-Holocene drove increased precipitation and widespread expansion of vegetation across the Sahara during the African humid period (AHP). While changes in atmospheric dynamics during this time have been a major focus of palaeoclimate modelling efforts, the transient nature of the shift back to the modern desert state at the end of this period is less well understood. Reconstructions reveal a spatially and temporally complex end of the AHP, with an earlier end in the north than in the south and in the east than in the west. Some records suggest a rather abrupt end, whereas others indicate a gradual decline in moisture availability. Here we investigate the end of the AHP based on a transient simulation of the last 7850 years with the comprehensive Earth system model MPI-ESM1.2. The model largely reproduces the time-transgressive end of the AHP evident in proxy data, and it indicates that it is due to the regionally varying dynamical controls on precipitation. The impact of the main rain-bringing systems, i.e. the summer monsoon and extratropical troughs, varies spatially, leading to heterogeneous seasonal rainfall cycles that impose regionally different responses to the Holocene insolation decrease. An increase in extratropical troughs that interact with the tropical mean flow and transport moisture to the western Sahara during the mid-Holocene delays the end of the AHP in that region. Along the coast, this interaction maintains humid conditions for a longer time than further inland. Drying in this area occurs when this interaction becomes too weak to sustain precipitation. In the lower latitudes of west Africa, where the rainfall is only influenced by the summer monsoon dynamics, the end of the AHP coincides with the retreat of the monsoonal rain belt. The model results clearly demonstrate that non-monsoonal dynamics can also play an important role in forming the precipitation signal and should therefore not be neglected in analyses of north African rainfall trends.

Description: The supereruption of Los Chocoyos (14.6°N, 91.2°W) in Guatemala ∼84kyr ago was one of the largest volcanic events of the past 100000 years. Recent petrologic data show that the eruption released very large amounts of climate-relevant sulfur and ozone-destroying chlorine and bromine gases (523±94Mt sulfur, 1200±156Mt chlorine, and 2±0.46Mt bromine). Using the Earth system model (ESM) of the Community Earth System Model version 2 (CESM2) coupled with the Whole Atmosphere Community Climate Model version 6 (WACCM6), we simulated the impacts of the sulfur- and halogen-rich Los Chocoyos eruption on the preindustrial Earth system. Our simulations show that elevated sulfate burden and aerosol optical depth (AOD) persists for 5 years in the model, while the volcanic halogens stay elevated for nearly 15 years. As a consequence, the eruption leads to a collapse of the ozone layer with global mean column ozone values dropping to 50DU (80% decrease) and leading to a 550% increase in surface UV over the first 5 years, with potential impacts on the biosphere. The volcanic eruption shows an asymmetric-hemispheric response with enhanced aerosol, ozone, UV, and climate signals over the Northern Hemisphere. Surface climate is impacted globally due to peak AOD of 〉6, which leads to a maximum surface cooling of 〉6K, precipitation and terrestrial net primary production decrease of 〉25%, and sea ice area increases of 40% in the first 3 years. Locally, a wetting (〉100%) and strong increase in net primary production (NPP) (〉700%) over northern Africa is simulated in the first 5 years and related to a southward shift of the Intertropical Convergence Zone (ITCZ) to the southern tropics. The ocean responds with pronounced El Niño conditions in the first 3 years that shift to the southern tropics and are coherent with the ITCZ change. Recovery to pre-eruption ozone levels and climate takes 15 years and 30 years, respectively. The long-lasting surface cooling is sustained by an immediate increase in the Arctic sea ice area, followed by a decrease in poleward ocean heat transport at 60°N which lasts up to 20 years. In contrast, when simulating Los Chocoyos conventionally by including sulfur and neglecting halogens, we simulate a larger sulfate burden and AOD, more pronounced surface climate changes, and an increase in column ozone. By comparing our aerosol chemistry ESM results to other supereruption simulations with aerosol climate models, we find a higher surface climate impact per injected sulfur amount than previous studies for our different sets of model experiments, since the CESM2(WACCM6) creates smaller aerosols with a longer lifetime, partly due to the interactive aerosol chemistry. As the model uncertainties for the climate response to supereruptions are very large, observational evidence from paleo archives and a coordinated model intercomparison would help to improve our understanding of the climate and environment response.

Description: Electromagnetic (EM) geophysical methods are well equipped to distinguish electrical resistivity contrasts between freshwater-saturated and seawater-saturated formations. Beneath the semi-arid, rapidly urbanizing island of Malta, offshore groundwater is an important potential resource but it is not known whether the regional mean sea-level aquifer (MSLA) extends offshore. To address this uncertainty, land-based alongshore and across-shore time-domain electromagnetic (TDEM) responses were acquired with the G-TEM instrument (Geonics Ltd., Mississauga, ON, Canada) and used to map the onshore structure of the aquifer. 1-D inversion results suggest that the onshore freshwater aquifer resides at 4–24 m depth, underlain by seawater-saturated formations. The freshwater aquifer thickens with distance from the coastline. We present 2D and 3D electromagnetic forward modeling based on finite-element (FE) analysis to further constrain the subsurface geometry of the onshore freshwater body. We interpret the high resistivity zones that as brackish water-saturated bodies are associated with the mean sea-level aquifer. Generally, time-domain electromagnetic (TDEM) results provide valuable onshore hydrogeological information, which can be augmented with marine and coastal transition-zone measurements to assess potential hydraulic continuity of terrestrial aquifers extending offshore.

Description: The geographic position, highly fluctuating sea temperatures and hypersalinity make Persian Gulf an extreme environment. Although this unique environment has high biodiversity dominated by invertebrates, its potential in marine biodiscovery has largely remained untapped. Herein, we aimed at a detailed analysis of the metabolome and bioactivity profiles of the marine sponge Axinella sinoxea collected from the northeast coast of the Persian Gulf in Iran. The crude extract and its Kupchan subextracts were tested in multiple in-house bioassays, and the crude extract and its CHCl3-soluble portion showed in vitro antibacterial activity against Methicillin-resistant Staphylococcus aureus (MRSA) and Enterococcus faecium (Efm). A molecular networking (MN)-based dereplication strategy by UPLC-MS/MS revealed the presence of phospholipids and steroids, while 1H NMR spectroscopy indicated the presence of additional metabolites, such as diketopiperazines (DKPs). Integrated MN and 1H NMR analyses on both the crude and CHCl3 extracts combined with an antibacterial activity-guided isolation approach afforded eight metabolites: a new diketopiperazine, (-)-cyclo(L-trans-Hyp-L-Ile) (8); a known diketopiperazine, cyclo(L-trans-Hyp-L-Phe) (7); two known phospholipids, 1-O-hexadecyl-sn-glycero-3-phosphocholine (1) and 1-O-octadecanoyl-sn-glycero-3-phosphocholine (2); two known steroids, 3β-hydroxycholest-5-ene-7,24-dione (3) and (22E)-3β-hydroxycholesta-5,22-diene-7,24-dione (4); two known monoterpenes, loliolide (5) and 5-epi-loliolide (6). The chemical structures of the isolates were elucidated by a combination of NMR spectroscopy, HRMS and [α]D analyses. All compounds were tested against MRSA and Efm, and compound 3 showed moderate antibacterial activity against MRSA (IC50 value 70 μg/mL). This is the first study that has dealt with chemical and bioactivity profiling of A. sinoxea leading to isolation and characterization of pure sponge metabolites.

Description: Leaf litter contributes to the functioning of aquatic ecosystems through allochthonous inputs of carbon, nitrogen, and other elements. Here, we examine leaf litter nutrient inputs and decomposition associated with four plant species using a mesocosm approach. Native sycamore fig Ficus sycomorus L., and silver cluster–leaf Terminalia sericea Burch. ex DC. decomposition dynamics were compared to invasive tickberry Lantana camara L. and guava Psidium guajava L., whereby phosphate, nitrate, nitrite, silicate, and ammonium releases were quantified over time. Leaf inputs significantly reduced pH, with reductions most marked by invasive L. camara. Conductivity was heightened by all leaf input treatments, except native T. sericea. Leaf inputs significantly affected all nutrient levels monitored in the water over time, except for silicate. In particular, leaf litter from invasive L. camara drove significantly increased nutrient concentrations compared to other native plant species, whilst effects of invasive P. guajava were less statistically clear. The end weights of the leaf litter demonstrated decomposition differences among the species types, following a decreasing order of P. guajava 〉 T. sericea 〉 F. sycomorus 〉 L. camara, further suggesting high organic inputs from invasive L. camara. The study results highlight that differential leaf litter decomposition rates of four plant species can play a significant role in nutrient release, in turn altering aquatic ecosystem productivity. However, these effects likely depend on species-specific differences, rather than between invasive–native species generally. Shifting terrestrial plant communities may alter aquatic community composition, but specific effects are likely associated with leaf traits

Description: We investigate the changes in terrestrial natural methane emissions between the Last Glacial Maximum (LGM) and preindustrial (PI) periods by performing time-slice experiments with a methane-enabled version of MPI-ESM, the Max Planck Institute Earth System Model. We consider all natural sources of methane except for emissions from wild animals and geological sources, i.e. emissions from wetlands, fires, and termites. Changes are dominated by changes in tropical wetland emissions, with mid-to-high-latitude wetlands playing a secondary role, and all other natural sources being of minor importance. The emissions are determined by the interplay of vegetation productivity, a function of CO2 and temperature; source area size, affected by sea level and ice sheet extent; and the state of the West African monsoon, with increased emissions from northern Africa during strong monsoon phases. We show that it is possible to explain the difference in atmospheric methane between LGM and PI purely by changes in emissions. As emissions more than double between LGM and PI, changes in the atmospheric lifetime of CH4, as proposed in other studies, are not required.

Description: A compilation of the published literature on dust content in terrestrial and marine sediment cores was synchronized with pollen data and speleothem growth phases on the Greenland Ice Core Chronology 2005 (GICC05) time axis. Aridity patterns for eight key areas of the global climate system have been reconstructed for the last 60 000 years. These records have different time resolutions and different dating methods, i.e. different types of stratigraphy. Nevertheless, all regions analysed in this study show humid conditions during early Marine Isotope Stage 3 (MIS3) and the early Holocene or deglaciation, but not always at the same time. Such discrepancies have been interpreted as regional effects, although stratigraphic uncertainties may affect some of the proposed interpretations. In comparison, most of the MIS2 interval becomes arid in all of the Northern Hemisphere records, but the peak arid conditions of the Last Glacial Maximum (LGM) and Heinrich event 1 differ in duration and intensity among regions. In addition, we also compare the aridity synthesis with modelling results using a global climate model (GCM). Indeed, geological archives and GCMs show agreement on the aridity pattern for the Holocene or deglaciation, for the LGM and for late MIS3.

Description: Terrestrial plants are known to “garden” the microbiota of their rhizosphere via released metabolites (that can attract beneficial microbes and deter pathogenic microbes). Such a “gardening” capacity is also known to be dynamic in plants. Although microbial “gardening” has been recently demonstrated for seaweeds, we do not know whether this capacity is a dynamic property in any aquatic flora like in terrestrial plants. Here, we tested the dynamic microbial “gardening” capacity of seaweeds using the model invasive red seaweed Agarophyton vermiculophyllum. Following an initial extraction of surface-associated metabolites (immediately after field collection), we conducted a long-term mesocosm experiment for 5 months to test the effect of two different salinities (low = 8.5 and medium = 16.5) on the microbial “gardening” capacity of the alga over time. We tested “gardening” capacity of A. vermiculophyllum originating from two different salinity levels (after 5 months treatments) in settlement assays against three disease causing pathogenic bacteria and seven protective bacteria. We also compared the capacity of the alga with field-collected samples. Abiotic factors like low salinity significantly increased the capacity of the alga to deter colonization by pathogenic bacteria while medium salinity significantly decreased the capacity of the alga over time when compared to field-collected samples. However, capacity to attract beneficial bacteria significantly decreased at both tested salinity levels when compared to field-collected samples. Dynamic microbial “gardening” capacity of a seaweed to attract beneficial bacteria and deter pathogenic bacteria is demonstrated for the first time. Such a dynamic capacity as found in the current study could also be applicable to other aquatic host–microbe interactions. Our results may provide an attractive direction of research towards manipulation of salinity and other abiotic factors leading to better defended A. vermiculophyllum towards pathogenic bacteria thereby enhancing sustained production of healthy A. vermiculophyllum in farms

Description: Ascidians and their associated microbiota are prolific producers of bioactive marine natural products. Recent culture-independent studies have revealed that the tunic of the solitary ascidian Ciona intestinalis (sea vase) is colonized by a diverse bacterial community, however, the biotechnological potential of this community has remained largely unexplored. In this study, we aimed at isolating the culturable microbiota associated with the tunic of C. intestinalis collected from the North and Baltic Seas, to investigate their antimicrobial and anticancer activities, and to gain first insights into their metabolite repertoire. The tunic of the sea vase was found to harbor a rich microbial community, from which 89 bacterial and 22 fungal strains were isolated. The diversity of the tunic-associated microbiota differed from that of the ambient seawater samples, but also between sampling sites. Fungi were isolated for the first time from the tunic of Ciona. The proportion of bioactive extracts was high, since 45% of the microbial extracts inhibited the growth of human pathogenic bacteria, fungi or cancer cell lines. In a subsequent bioactivity-and metabolite profiling-based approach, seven microbial extracts were prioritized for in-depth chemical investigations. Untargeted metabolomics analyses of the selected extracts by a UPLC-MS/MS-based molecular networking approach revealed a vast chemical diversity with compounds assigned to 22 natural product families, plus many metabolites that remained unidentified. This initial study indicates that bacteria and fungi associated with the tunic of C. intestinalis represent an untapped source of putatively new marine natural products with pharmacological relevance.

Description: Here we present the first radiometric age and geochemical (major and trace element and isotope) data for samples from the Hydrographer Ridge, a back arc volcano of the Kurile Island Arc, and a newly discovered chain of volcanoes (“Sonne Volcanoes”) on the northwestern continental slope of the Kurile Basin on the opposite side of the arc. The 40Ar/39Ar age and geochemical data show that Hydrographer Ridge (3.2–3.3 Ma) and the “Sonne Volcanoes” (25.3–25.9 Ma) have very similar trace element and isotope characteristics to those of the Kurile Island Arc, indicating derivation from a common magma source. We conclude that the age of the “Sonne Volcanoes” marks the time of opening of the Kurile Basin, implying slow back arc spreading rates of 1.3–1.8 cm/y. Combined with published data from the Kurile fore arc, our data suggest that the processes of subduction, Kurile Basin opening and frontal arc extension occurred synchronously and that extension in the rear part and in the frontal part of the Kurile Island Arc must have been triggered by the same mechanism.

Description: Massive fouling by the invasive ascidian Ciona intestinalis in Prince Edward Island (PEI, Canada) has been causing devastating losses to the local blue mussel farms. In order to gain first insights into so far unexplored factors that may contribute to the invasiveness of C. intestinalis in PEI, we undertook comparative microbiome and metabolome studies on specific tissues from C. intestinalis populations collected in invaded (PEI) and native regions (Helgoland and Kiel, Germany). Microbial community analyses and untargeted metabolomics revealed clear location- and tissue-specific patterns showing that biogeography and the sampled tissue shape the microbiome and metabolome of C. intestinalis. Moreover, we observed higher microbial and chemical diversity in C. intestinalis from PEI than in the native populations. Bacterial OTUs specific to C. intestinalis from PEI included Cyanobacteria (e.g., Leptolyngbya sp.) and Rhodobacteraceae (e.g., Roseobacter sp.), while populations from native sampling sites showed higher abundances of e.g., Firmicutes (Helgoland) and Epsilonproteobacteria (Kiel). Altogether 121 abundant metabolites were putatively annotated in the global ascidian metabolome, of which 18 were only detected in the invasive PEI population (e.g., polyketides and terpenoids), while six (e.g., sphingolipids) or none were exclusive to the native specimens from Helgoland and Kiel, respectively. Some identified bacteria and metabolites reportedly possess bioactive properties (e.g., antifouling and antibiotic) that may contribute to the overall fitness of C. intestinalis. Hence, this first study provides a basis for future studies on factors underlying the global invasiveness of Ciona species

Description: The Global Ocean Data Analysis Project (GLODAP) is a synthesis effort providing regular compilations of surface-to-bottom ocean biogeochemical data, with an emphasis on seawater inorganic carbon chemistry and related variables determined through chemical analysis of seawater samples. GLODAPv2.2020 is an update of the previous version, GLODAPv2.2019. The major changes are data from 106 new cruises added, extension of time coverage to 2019, and the inclusion of available (also for historical cruises) discrete fugacity of CO2 (fCO2) values in the merged product files. GLODAPv2.2020 now includes measurements from more than 1.2 million water samples from the global oceans collected on 946 cruises. The data for the 12 GLODAP core variables (salinity, oxygen, nitrate, silicate, phosphate, dissolved inorganic carbon, total alkalinity, pH, CFC-11, CFC-12, CFC-113, and CCl4) have undergone extensive quality control with a focus on systematic evaluation of bias. The data are available in two formats: (i) as submitted by the data originator but updated to WOCE exchange format and (ii) as a merged data product with adjustments applied to minimize bias. These adjustments were derived by comparing the data from the 106 new cruises with the data from the 840 quality-controlled cruises of the GLODAPv2.2019 data product using crossover analysis. Comparisons to empirical algorithm estimates provided additional context for adjustment decisions; this is new to this version. The adjustments are intended to remove potential biases from errors related to measurement, calibration, and data-handling practices without removing known or likely time trends or variations in the variables evaluated. The compiled and adjusted data product is believed to be consistent to better than 0.005 in salinity, 1 % in oxygen, 2 % in nitrate, 2 % in silicate, 2 % in phosphate, 4 µmol kg−1 in dissolved inorganic carbon, 4 µmol kg−1 in total alkalinity, 0.01–0.02 in pH (depending on region), and 5 % in the halogenated transient tracers. The other variables included in the compilation, such as isotopic tracers and discrete fCO2, were not subjected to bias comparison or adjustments.

Description: Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere in a changing climate – the “global carbon budget” – is important to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe and synthesize data sets and methodology to quantify the five major components of the global carbon budget and their uncertainties. Fossil CO2 emissions (EFOS) are based on energy statistics and cement production data, while emissions from land-use change (ELUC), mainly deforestation, are based on land use and land-use change data and bookkeeping models. Atmospheric CO2 concentration is measured directly and its growth rate (GATM) is computed from the annual changes in concentration. The ocean CO2 sink (SOCEAN) and terrestrial CO2 sink (SLAND) are estimated with global process models constrained by observations. The resulting carbon budget imbalance (BIM), the difference between the estimated total emissions and the estimated changes in the atmosphere, ocean, and terrestrial biosphere, is a measure of imperfect data and understanding of the contemporary carbon cycle. All uncertainties are reported as ±1σ. For the last decade available (2010–2019), EFOS was 9.6 ± 0.5 GtC yr−1 excluding the cement carbonation sink (9.4 ± 0.5 GtC yr−1 when the cement carbonation sink is included), and ELUC was 1.6 ± 0.7 GtC yr−1. For the same decade, GATM was 5.1 ± 0.02 GtC yr−1 (2.4 ± 0.01 ppm yr−1), SOCEAN 2.5 ±  0.6 GtC yr−1, and SLAND 3.4 ± 0.9 GtC yr−1, with a budget imbalance BIM of −0.1 GtC yr−1 indicating a near balance between estimated sources and sinks over the last decade. For the year 2019 alone, the growth in EFOS was only about 0.1 % with fossil emissions increasing to 9.9 ± 0.5 GtC yr−1 excluding the cement carbonation sink (9.7 ± 0.5 GtC yr−1 when cement carbonation sink is included), and ELUC was 1.8 ± 0.7 GtC yr−1, for total anthropogenic CO2 emissions of 11.5 ± 0.9 GtC yr−1 (42.2 ± 3.3 GtCO2). Also for 2019, GATM was 5.4 ± 0.2 GtC yr−1 (2.5 ± 0.1 ppm yr−1), SOCEAN was 2.6 ± 0.6 GtC yr−1, and SLAND was 3.1 ± 1.2 GtC yr−1, with a BIM of 0.3 GtC. The global atmospheric CO2 concentration reached 409.85 ± 0.1 ppm averaged over 2019. Preliminary data for 2020, accounting for the COVID-19-induced changes in emissions, suggest a decrease in EFOS relative to 2019 of about −7 % (median estimate) based on individual estimates from four studies of −6 %, −7 %, −7 % (−3 % to −11 %), and −13 %. Overall, the mean and trend in the components of the global carbon budget are consistently estimated over the period 1959–2019, but discrepancies of up to 1 GtC yr−1 persist for the representation of semi-decadal variability in CO2 fluxes. Comparison of estimates from diverse approaches and observations shows (1) no consensus in the mean and trend in land-use change emissions over the last decade, (2) a persistent low agreement between the different methods on the magnitude of the land CO2 flux in the northern extra-tropics, and (3) an apparent discrepancy between the different methods for the ocean sink outside the tropics, particularly in the Southern Ocean. This living data update documents changes in the methods and data sets used in this new global carbon budget and the progress in understanding of the global carbon cycle compared with previous publications of this data set (Friedlingstein et al., 2019; Le Quéré et al., 2018b, a, 2016, 2015b, a, 2014, 2013). The data presented in this work are available at https://doi.org/10.18160/gcp-2020 (Friedlingstein et al., 2020).

Description: Seamounts are abundant and prominent features on the deep-sea floor and intersperse with the nodule fields of the Clarion-Clipperton Fracture Zone (CCZ). There is a particular interest in characterising the fauna inhabiting seamounts in the CCZ because they are the only other ecosystem in the region to provide hard substrata besides the abundant nodules on the soft-sediment abyssal plains. It has been hypothesised that seamounts could provide refuge for organisms during deep-sea mining actions or that they could play a role in the (re-)colonisation of the disturbed nodule fields. This hypothesis is tested by analysing video transects in both ecosystems, assessing megafauna composition and abundance. Nine video transects (ROV dives) from two different license areas and one Area of Particular Environmental Interest in the eastern CCZ were analysed. Four of these transects were carried out as exploratory dives on four different seamounts in order to gain first insights into megafauna composition. The five other dives were carried out in the neighbouring nodule fields in the same areas. Variation in community composition observed among and along the video transects was high, with little morphospecies overlap along intra-ecosystem transects. Despite the observation of considerable faunal variations within each ecosystem, differences between seamounts and nodule fields prevailed, showing significantly different species associations characterising them, thus calling into question their use as a possible refuge area.

Description: The eastern tropical South Pacific (ETSP) represents one of the most productive areas in the ocean that is characterized by a pronounced oxygen minimum zone (OMZ). Particulate organic matter (POM) that sinks out of the euphotic zone is supplied to the anoxic sediments and utilized by microbial communities. The degradation of POM is associated with dissolved organic matter (DOM) production and reworking. The release of recalcitrant DOM to the overlying waters may represent an important organic matter escape mechanism from remineralization within sediments but received little attention in OMZ regions so far. Here, we combine measurements of dissolved organic carbon (DOC) and dissolved organic nitrogen (DON) with DOM optical properties in the form of chromophoric (CDOM) and fluorescent (FDOM) DOM from pore waters and near-bottom waters of the ETSP off Peru. We evaluate diffusion–driven fluxes and net in situ fluxes of DOC and DON in order to investigate processes affecting DOM cycling at the sediment–water interface along a transect 12° S. To our knowledge, these are the first data for sediment release of DON and pore water CDOM and FDOM for the ETSP off Peru. Pore-water DOC and DON accumulated with increasing sediment depth, suggesting an imbalance between DOM production and remineralization within sediments. High DON accumulation resulted in very low pore water DOC / DON ratios (〉 1) which could be caused by either an "imbalance" in DOC and DON remineralization, or to the presence of an additional nitrogen source. Diffusion driven fluxes of DOC and DON exhibited high spatial variability. They varied from 0.2–0.1 mmol m−2 d−1 to 2.52–1.3 mmol m−2 d−1 and from −0.042–0.02 mmol m−2 d−1 to 3.32–1.7 mmol m−2 d−1, respectively. Generally low net in situ DOC and DON fluxes as well as steepening of spectral slope (S) of CDOM and accumulation of humic-like FDOM at the near-bottom waters over time indicated active microbial DOM utilization at the sediment–water interface, potentially stimulated by nitrate (NO3−) and nitrite (NO2−). The microbial DOC utilization rates, estimated in our study, may be sufficient to support denitrification rates of 0.2–1.4 mmol m−2 d−1, suggesting that sediment release of DOM contributes substantially to nitrogen loss processes in the ETSP off Peru.

Description: Marine aggregates are the vector for biogenically bound carbon and nutrients from the euphotic zone to the interior of the oceans. To improve the representation of this biological carbon pump in the global biogeochemical HAMburg Ocean Carbon Cycle (HAMOCC) model, we implemented a novel Microstructure, Multiscale, Mechanistic, Marine Aggregates in the Global Ocean (M4AGO) sinking scheme. M4AGO explicitly represents the size, microstructure, heterogeneous composition, density and porosity of aggregates and ties ballasting mineral and particulate organic carbon (POC) fluxes together. Additionally, we incorporated temperature-dependent remineralization of POC. We compare M4AGO with the standard HAMOCC version, where POC fluxes follow a Martin curve approach with (i) linearly increasing sinking velocity with depth and (ii) temperature-independent remineralization. Minerals descend separately with a constant speed. In contrast to the standard HAMOCC, M4AGO reproduces the latitudinal pattern of POC transfer efficiency, as recently constrained by Weber et al. (2016). High latitudes show transfer efficiencies of ≈0.25±0.04, and the subtropical gyres show lower values of about 0.10±0.03. In addition to temperature as a driving factor for remineralization, diatom frustule size co-determines POC fluxes in silicifier-dominated ocean regions, while calcium carbonate enhances the aggregate excess density and thus sinking velocity in subtropical gyres. Prescribing rising carbon dioxide (CO2) concentrations in stand-alone runs (without climate feedback), M4AGO alters the regional ocean atmosphere CO2 fluxes compared to the standard model. M4AGO exhibits higher CO2 uptake in the Southern Ocean compared to the standard run, while in subtropical gyres, less CO2 is taken up. Overall, the global oceanic CO2 uptake remains the same. With the explicit representation of measurable aggregate properties, M4AGO can serve as a test bed for evaluating the impact of aggregate-associated processes on global biogeochemical cycles and, in particular, on the biological carbon pump.

Description: With increasing pressure to extract minerals from the deep-sea bed, understanding the ecological and evolutionary processes that limit the spatial distribution of species is critical to assessing ecosystem resilience to mining impacts. The aim of our study is to gain a better knowledge about the abyssal isopod crustacean fauna of the central Pacific manganese nodule province (Clarion–Clipperton Fracture Zone, CCZ). In total, we examined 22 epibenthic sledge (EBS) samples taken at five abyssal areas located in the central northern Pacific including four contracting areas and one Area of Particular Environmental Interest (APEI3). Additional samples come from the DISturbance and reCOLonization experiment (DISCOL) area situated in the Peru Basin, southeastern Pacific. Using an integrative approach that combined morphological and genetic methods with species delimitation analyses (SDs) we assessed patterns of species range size, diversity, and community composition for four different isopod families (Munnopsidae Lilljeborg, 1864; Desmosomatidae Sars, 1897; Haploniscidae Hansen, 1916; and Macrostylidae Hansen, 1916) displaying different dispersal capacities as adults. Isopods are brooders, so their distribution and connectivity cannot be explained by larval dispersal but rather by adult locomotion. In particular, our objectives were to (1) identify potential differences in the distributional ranges of isopod families relative to their locomotory potential and to (2) evaluate the representativeness of the APEI for the preservation of regional biodiversity in the CCZ following mining disturbances. From 619 specimens, our SD analysis could distinguish 170 species, most of which were new to science (94.1 %). We found that increased locomotory ability correlated with higher species diversity with 9 species of Macrostylidae, 23 of Haploniscidae, 52 of Desmosomatidae, and 86 of Munnopsidae. This is supported by family-level rarefaction analyses. As expected, we found the largest species ranges in the families with swimming abilities, with a maximum recorded species range of 5245 and 4480 km in Munnopsidae and Desmosomatidae, respectively. The less motile Haploniscidae and Macrostylidae had maximal species ranges of 1391 and 1440 km, respectively. Overall, rarefaction analyses indicated that species richness did not vary much between areas, but the real number of species was still not sufficiently sampled. This is also indicated by the large proportion of singletons (40.5 %) found in this study. The investigated contractor areas in the CCZ were more similar in species composition and had a higher proportion of shared species between each other than the closely located APEI3 and the distantly located DISCOL area. In fact, the DISCOL area, located in the Peru Basin, had more species in common with the core CCZ areas than APEI3. In this regard, APEI3 does not appear to be representative as serving as a reservoir for the fauna of the investigated contractor areas, at least for isopods, as it has a different species composition. Certainly, more data from other APEIs, as well as preservation reference zones within contractor areas, are urgently needed in order to assess their potential as resources of recolonization of impacted seabed.

Description: Northern peatlands have been a persistent natural carbon sink since the Last Glacial Maximum. The continued growth and expansion of these carbon-rich ecosystems could offset a large portion of anthropogenic carbon emissions before the end of the present interglacial period. Here we used an impeded drainage model and gridded data on the depth to bedrock and the fraction of histosol-type soils to evaluate the limits to the growth of northern peatland carbon stocks. Our results show that the potential carbon stock in northern peatlands could reach a total of 875±125 Pg C before the end of the present interglacial, which could, as a result, remove 330±200 Pg C of carbon from the atmosphere. We argue that northern peatlands, together with the oceans, will potentially play an important role in reducing the atmospheric carbon dioxide concentration over the next 5000 years.

Description: The Ligurian Basin is located in the Mediterranean Sea to the north-west of Corsica at the transition from the Western Alpine orogen to the Apennine system and was generated by the south-eastward trench retreat of the Apennines–Calabrian subduction zone. Late-Oligocene-to-Miocene rifting caused continental extension and subsidence, leading to the opening of the basin. Yet it remains unclear if rifting caused continental break-up and seafloor spreading. To reveal its lithospheric architecture, we acquired a 130 km long seismic refraction and wide-angle reflection profile in the Ligurian Basin. The seismic line was recorded in the framework of SPP2017 4D-MB, a Priority Programme of the German Research Foundation (DFG) and the German component of the European AlpArray initiative, and trends in a NE–SW direction at the centre of the Ligurian Basin, roughly parallel to the French coastline. The seismic data were recorded on the newly developed GEOLOG recorder, designed at GEOMAR, and are dominated by sedimentary refractions and show mantle Pn arrivals at offsets of up to 70 km and a very prominent wide-angle Mohorovičić discontinuity (Moho) reflection. The main features share several characteristics (e.g. offset range, continuity) generally associated with continental settings rather than documenting oceanic crust emplaced by seafloor spreading. Seismic tomography results are complemented by gravity data and yield a ∼ 6–8 km thick sedimentary cover and the seismic Moho at 11–13 km depth below the sea surface. Our study reveals that the oceanic domain does not extend as far north as previously assumed. Whether Oligocene–Miocene extension led to extremely thinned continental crust or exhumed subcontinental mantle remains unclear. A low grade of mantle serpentinisation indicates a high rate of syn-rift sedimentation. However, rifting failed before oceanic spreading was initiated, and continental crust thickens towards the NE within the northern Ligurian Basin.

Description: We present a new near-global coupled biogeochemical ocean-circulation model configuration. The configuration features a horizontal discretization with a grid spacing of less than 11 km in the Southern Ocean and gradually coarsens in meridional direction to more than 200 km at 64∘ N, where the model is bounded by a solid wall. The underlying code framework is the Geophysical Fluid Dynamics Laboratory (GFDL)'s Modular Ocean Model coupled to the Biogeochemistry with Light, Iron, Nutrients and Gases (BLING) ecosystem model of Galbraith et al. (2010). The configuration is unique in that it features both a relatively equilibrated oceanic carbon inventory and an eddying ocean circulation based on a realistic model geometry/bathymetry – a combination that has been precluded by prohibitive computational cost in the past. Results from a simulation with climatological forcing and a sensitivity experiment with increasing winds suggest that the configuration is sufficiently equilibrated to explore Southern Ocean carbon uptake dynamics on decadal timescales. The configuration is dubbed MOMSO, a Modular Ocean Model Southern Ocean configuration.

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: The project MarParCloud (Marine biological production, organic aerosol Particles and marine Clouds: a process chain) aims to improve our understanding of the genesis, modification and impact of marine organic matter (OM) from its biological production, to its export to marine aerosol particles and, finally, to its ability to act as ice-nucleating particles (INPs) and cloud condensation nuclei (CCN). A field campaign at the Cape Verde Atmospheric Observatory (CVAO) in the tropics in September–October 2017 formed the core of this project that was jointly performed with the project MARSU (MARine atmospheric Science Unravelled). A suite of chemical, physical, biological and meteorological techniques was applied, and comprehensive measurements of bulk water, the sea surface microlayer (SML), cloud water and ambient aerosol particles collected at a ground-based and a mountain station took place. Key variables comprised the chemical characterization of the atmospherically relevant OM components in the ocean and the atmosphere as well as measurements of INPs and CCN. Moreover, bacterial cell counts, mercury species and trace gases were analyzed. To interpret the results, the measurements were accompanied by various auxiliary parameters such as air mass back-trajectory analysis, vertical atmospheric profile analysis, cloud observations and pigment measurements in seawater. Additional modeling studies supported the experimental analysis. During the campaign, the CVAO exhibited marine air masses with low and partly moderate dust influences. The marine boundary layer was well mixed as indicated by an almost uniform particle number size distribution within the boundary layer. Lipid biomarkers were present in the aerosol particles in typical concentrations of marine background conditions. Accumulation- and coarse-mode particles served as CCN and were efficiently transferred to the cloud water. The ascent of ocean-derived compounds, such as sea salt and sugar-like compounds, to the cloud level, as derived from chemical analysis and atmospheric transfer modeling results, denotes an influence of marine emissions on cloud formation. Organic nitrogen compounds (free amino acids) were enriched by several orders of magnitude in submicron aerosol particles and in cloud water compared to seawater. However, INP measurements also indicated a significant contribution of other non-marine sources to the local INP concentration, as (biologically active) INPs were mainly present in supermicron aerosol particles that are not suggested to undergo strong enrichment during ocean–atmosphere transfer. In addition, the number of CCN at the supersaturation of 0.30 % was about 2.5 times higher during dust periods compared to marine periods. Lipids, sugar-like compounds, UV-absorbing (UV: ultraviolet) humic-like substances and low-molecular-weight neutral components were important organic compounds in the seawater, and highly surface-active lipids were enriched within the SML. The selective enrichment of specific organic compounds in the SML needs to be studied in further detail and implemented in an OM source function for emission modeling to better understand transfer patterns, the mechanisms of marine OM transformation in the atmosphere and the role of additional sources. In summary, when looking at particulate mass, we see oceanic compounds transferred to the atmospheric aerosol and to the cloud level, while from a perspective of particle number concentrations, sea spray aerosol (i.e., primary marine aerosol) contributions to both CCN and INPs are rather limited.

Description: During the Late Pleistocene and Holocene retreat of paleo-ice sheets in North America and Europe, vast proglacial lakes existed along the land terminating margins. These proglacial lakes impacted ice sheet dynamics by imposing boundary conditions analogous to a marine terminating margin. Such lacustrine boundary conditions cause changes in the ice sheet geometry, stress balance and frontal ablation and therefore affect the mass balance of the entire ice sheet. Despite this, dynamically evolving proglacial lakes have rarely been considered in detail in ice sheet modeling endeavors. In this study, we describe the implementation of an adaptive lake boundary in the Parallel Ice Sheet Model (PISM), which we call PISM-LakeCC. We test our model with a simplified glacial retreat setup of the Laurentide Ice Sheet (LIS). By comparing the experiments with lakes to control runs with no lakes, we show that the presence of proglacial lakes locally enhances the ice flow, which leads to a lowering of the ice sheet surface. In some cases, this also results in an advance of the ice margin and the emergence of ice lobes. In the warming climate, increased melting on the lowered ice surface drives the glacial retreat. For the LIS, the presence of lakes triggers a process similar to marine ice sheet instability, which caused the collapse of the ice saddle over Hudson Bay. In the control experiments without lakes, Hudson Bay is still glaciated when the climate reaches present-day (PD) conditions. The results of our study demonstrate that glacio-lacustrine interactions play a significant role in the retreat of land terminating ice sheet margins.

Description: Characterizing the variability across timescales is important for understanding the underlying dynamics of the Earth system. It remains challenging to do so from palaeoclimate archives since they are more often than not irregular, and traditional methods for producing timescale-dependent estimates of variability, such as the classical periodogram and the multitaper spectrum, generally require regular time sampling. We have compared those traditional methods using interpolation with interpolation-free methods, namely the Lomb–Scargle periodogram and the first-order Haar structure function. The ability of those methods to produce timescale-dependent estimates of variability when applied to irregular data was evaluated in a comparative framework, using surrogate palaeo-proxy data generated with realistic sampling. The metric we chose to compare them is the scaling exponent, i.e. the linear slope in log-transformed coordinates, since it summarizes the behaviour of the variability across timescales. We found that, for scaling estimates in irregular time series, the interpolation-free methods are to be preferred over the methods requiring interpolation as they allow for the utilization of the information from shorter timescales which are particularly affected by the irregularity. In addition, our results suggest that the Haar structure function is the safer choice of interpolation-free method since the Lomb–Scargle periodogram is unreliable when the underlying process generating the time series is not stationary. Given that we cannot know a priori what kind of scaling behaviour is contained in a palaeoclimate time series, and that it is also possible that this changes as a function of timescale, it is a desirable characteristic for the method to handle both stationary and non-stationary cases alike.

Description: The Parallel Ice Sheet Model (PISM) is applied to the Antarctic Ice Sheet over the last two glacial cycles (≈210 000 years) with a resolution of 16 km. An ensemble of 256 model runs is analyzed in which four relevant model parameters have been systematically varied using full-factorial parameter sampling. Parameters and plausible parameter ranges have been identified in a companion paper (Albrecht et al., 2020) and are associated with ice dynamics, climatic forcing, basal sliding and bed deformation and represent distinct classes of model uncertainties. The model is scored against both modern and geologic data, including reconstructed grounding-line locations, elevation–age data, ice thickness, surface velocities and uplift rates. An aggregated score is computed for each ensemble member that measures the overall model–data misfit, including measurement uncertainty in terms of a Gaussian error model (Briggs and Tarasov, 2013). The statistical method used to analyze the ensemble simulation results follows closely the simple averaging method described in Pollard et al. (2016). This analysis reveals clusters of best-fit parameter combinations, and hence a likely range of relevant model and boundary parameters, rather than individual best-fit parameters. The ensemble of reconstructed histories of Antarctic Ice Sheet volumes provides a score-weighted likely range of sea-level contributions since the Last Glacial Maximum (LGM) of 9.4±4.1 m (or 6.5±2.0×106km3 ), which is at the upper range of most previous studies. The last deglaciation occurs in all ensemble simulations after around 12 000 years before present and hence after the meltwater pulse 1A (MWP1a). Our ensemble analysis also provides an estimate of parametric uncertainty bounds for the present-day state that can be used for PISM projections of future sea-level contributions from the Antarctic Ice Sheet.

Description: Uncertainties in projections from Earth system models (ESMs) are associated to a large degree with the imperfect representation of the marine plankton ecosystem, in particular the physiology of primary and secondary producers. Here we describe the implementation of an optimality-based plankton-ecosystem model (OPEM) with variable C:N:P stoichiometry in the University of Victoria ESM (UVic) and the behaviour of two calibrated reference configurations, which differ in the assumed temperature dependence of diazotrophs. Predicted tracer distributions of oxygen and dissolved inorganic nutrients are similar to those of an earlier fixed-stoichiometry model (Keller et al., 2012). Compared to the classic fixed-stoichiometry model, OPEM is closer to recent satellite-based estimates of net community production (NCP), despite overestimating net primary production (NPP), can better reproduce deep-ocean gradients in the NO3:PO43− ratio, and partially explains observed patterns of particulate C:N:P in the surface ocean. Allowing diazotrophs to grow (but not necessarily fix N2) at similar temperatures as other phytoplankton results in a better representation of surface Chl and NPP in the Arctic and Antarctic Oceans. Deficiencies of our calibrated OPEM configurations may serve as a magnifying glass for shortcomings in global biogeochemical models and hence guide future model development. The overestimation of NPP at low latitudes indicates the need for improved representations of temperature effects on biotic processes, as well as phytoplankton community composition, which may be represented by locally-varying parameters based on suitable trade-offs. Discrepancies between observed and predicted vertical gradients in particulate C:N:P ratios suggest the need to include preferential P remineralisation, which could also benefit the representation of N2 fixation. While OPEM yields a much improved distribution of surface N* (NO3− 16·PO43− + 2.9 mmol m−3), it still fails to reproduce observed N* in the Arctic, possibly related to a mis-representation of the phytoplankton community there and the lack of benthic denitrification in the model. Coexisting ordinary and diazotrophic phytoplankton can exert strong control on N* in our simulations, which questions the interpretation of N* as reflecting the balance of N2 fixation and denitrification.

Description: Simulations of the glacial–interglacial history of the Antarctic Ice Sheet provide insights into dynamic threshold behavior and estimates of the ice sheet's contributions to global sea-level changes for the past, present and future. However, boundary conditions are weakly constrained, in particular at the interface of the ice sheet and the bedrock. Also climatic forcing covering the last glacial cycles is uncertain, as it is based on sparse proxy data. We use the Parallel Ice Sheet Model (PISM) to investigate the dynamic effects of different choices of input data, e.g., for modern basal heat flux or reconstructions of past changes of sea level and surface temperature. As computational resources are limited, glacial-cycle simulations are performed using a comparably coarse model grid of 16 km and various parameterizations, e.g., for basal sliding, iceberg calving, or for past variations in precipitation and ocean temperatures. In this study we evaluate the model's transient sensitivity to corresponding parameter choices and to different boundary conditions over the last two glacial cycles and provide estimates of involved uncertainties. We also discuss isolated and combined effects of climate and sea-level forcing. Hence, this study serves as a “cookbook” for the growing community of PISM users and paleo-ice sheet modelers in general. For each of the different model uncertainties with regard to climatic forcing, ice and Earth dynamics, and basal processes, we select one representative model parameter that captures relevant uncertainties and motivates corresponding parameter ranges that bound the observed ice volume at present. The four selected parameters are systematically varied in a parameter ensemble analysis, which is described in a companion paper.

Description: The surface mass balance scheme dEBM (diurnal Energy Balance Model) provides a novel interface between the atmosphere and land ice for Earth system modeling, which is based on the energy balance of glaciated surfaces. In contrast to empirical schemes, dEBM accounts for changes in the Earth’s orbit and atmospheric composition. The scheme only requires monthly atmospheric forcing (precipitation, temperature, shortwave and longwave radiation, and cloud cover). It is also computationally inexpensive, which makes it particularly suitable to investigate the ice sheets' response to long-term climate change. After calibration and validation, we analyze the surface mass balance of the Greenland Ice Sheet (GrIS) based on climate simulations representing two warm climate states: a simulation of the mid-Holocene (approximately 6000 years before present) and a climate projection based on an extreme emission scenario which extends to the year 2100. The former period features an intensified summer insolation while the 21st century is characterized by reduced outgoing longwave radiation. Specifically, we investigate whether the temperature–melt relationship, as used in empirical temperature-index methods, remains stable under changing insolation and atmospheric composition. Our results indicate that the temperature–melt relation is sensitive to changes in insolation on orbital timescales but remains mostly invariant under the projected warming climate of the 21st century.

Description: Mastermind is famous two-player game. The first player (codemaker) chooses a secret code which the second player (codebreaker) is supposed to crack within a minimum number of code guesses (queries). Therefore, the codemaker’s duty is to help the codebreaker by providing a well-defined error measure between the secret code and the guessed code after each query. We consider a variant, called Yes-No AB-Mastermind, where both secret code and queries must be repetition-free and the provided information by the codemaker only indicates if a query contains any correct position at all. For this Mastermind version with n positions and 𝑘≥𝑛 colors and ℓ:=𝑘+1−𝑛, we prove a lower bound of ∑𝑘𝑗=ℓlog2𝑗 and an upper bound of 𝑛log2𝑛+𝑘 on the number of queries necessary to break the secret code. For the important case 𝑘=𝑛, where both secret code and queries represent permutations, our results imply an exact asymptotic complexity of Θ(𝑛log𝑛) queries.

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: Ice rises and ice rumples are locally grounded features found in coastal Antarctica and are surrounded by otherwise freely floating ice shelves. An ice rise has an independent flow regime, whereas the flow regime of an ice rumple conforms to that of the ice shelf and merely slows the flow of ice. In both cases, local highs in the bathymetry are in contact with the ice shelf from below, thereby regulating the large-scale ice flow, with implications for the upstream continental grounding line position. This buttressing effect, paired with the suitability of ice rises as a climate archive, necessitates a better understanding of the transition between ice rise and ice rumple, their evolution in response to a change in sea level, and their dynamic interaction with the surrounding ice shelf. We investigate this behaviour using a three-dimensional full Stokes ice flow model with idealised ice rises and ice rumples. The simulations span end-member basal friction scenarios of almost stagnant and fully sliding ice at the ice–bed interface. We analyse the coupling with the surrounding ice shelf by comparing the deviations between the non-local full Stokes surface velocities and the local shallow ice approximation (SIA). Deviations are generally high at the ice divides and small on the lee sides. On the stoss side, where ice rise and ice shelf have opposing flow directions, deviations can be significant. Differences are negligible in the absence of basal sliding where the corresponding steady-state ice rise is larger and develops a fully independent flow regime that is well described by SIA. When sea level is increased, and a transition from ice rise to ice rumple is approached, the divide migration is more abrupt the higher the basal friction. In each scenario, the transition occurs after the stoss-side grounding line has moved over the bed high and is positioned on a retrograde slope. We identify a hysteretic response of ice rises and ice rumples to changes in sea level, with grounded area being larger in a sea-level-increase scenario than in a sea-level-decrease scenario. This hysteresis shows not only irreversibility following an equal increase and subsequent decrease in sea level but also that the perturbation history is important when the ice rise or ice rumple geometry is not known. The initial grounded area needs to be carefully considered, as this will determine the formation of either an ice rise or an ice rumple, thereby causing different buttressing effects.

Description: Ice shelves surrounding the Antarctic perimeter moderate ice discharge towards the ocean through buttressing. Ice-shelf evolution and integrity depend on the local surface accumulation, basal melting and on the spatially variable ice-shelf viscosity. These components of ice-shelf mass balance are often poorly constrained by observations and introduce uncertainties in ice-sheet projections. Isochronal radar stratigraphy is an observational archive for the atmospheric, oceanographic and ice-flow history of ice shelves. Here, we predict the stratigraphy of locally accumulated ice on ice shelves with a kinematic forward model for a given atmospheric and oceanographic scenario. This delineates the boundary between local meteoric ice (LMI) and continental meteoric ice (CMI). A large LMI to CMI ratio hereby marks ice shelves whose buttressing strength is more sensitive to changes in atmospheric precipitation patterns. A mismatch between the steady-state predictions of the kinematic forward model and observations from radar can highlight inconsistencies in the atmospheric and oceanographic input data or be an indicator for a transient ice-shelf history not accounted for in the model. We discuss pitfalls in numerical diffusion when calculating the age field and validate the kinematic model with the full Stokes ice-flow model Elmer/Ice. The Roi Baudouin Ice Shelf (East Antarctica) serves as a test case for this approach. There, we find a significant east–west gradient in the LMI / CMI ratio. The steady-state predictions concur with observations on larger spatial scales (〉10 km), but deviations on smaller scales are significant, e.g., because local surface accumulation patterns near the grounding zone are underestimated in Antarctic-wide estimates. Future studies can use these mismatches to optimize the input data or to pinpoint transient signatures in the ice-shelf history using the ever growing archive of radar observations of internal ice stratigraphy.

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 Atlantic Multidecadal Oscillation (AMO) is coherently linked to climate variations over many parts of the globe. Despite recent achievements, the mechanism by which the AMO influences regional precipitation over South America is not well understood. In this study, we isolate the atmospheric response to the AMO using a water isotope-enabled version of the Community Earth System Model version 1.2 (iCESM1.2) and determine its influence on (sub)tropical South American regional precipitation and atmospheric circulation. The results suggest an interhemispheric seesaw in Hadley circulation strength and that the section of the Atlantic Hadley cell is marked by a stronger upward air component south of the equator during the cold AMO phase. We also find that the precipitation anomalies over (sub)tropical South America during AMO phases are mainly related to changes in the Atlantic Intertropical Convergence Zone (ITCZ) core strength, where in the cold (warm) AMO phase the core region strengthens (weakens) from February to July, while from July to November the core region weakens (strengthens). Our results stress the importance of acknowledging the dynamics of season- and regional-dependent ITCZ responses, as they are sufficient to produce observed AMO-related signals even in the absence of marked changes in the ITCZ position.

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: Military munitions from World War I and II dumped at the seafloor are a threat to the marine environment and its users. Decades of saltwater exposure make the explosives fragile and difficult to dispose of. If required, the munition is blast-in-place. In August 2019, 42 ground mines were detonated in a controlled manner underwater during a NATO maneuver in the German Natura2000 Special Area of Conservation Fehmarn Belt, the Baltic Sea. In June 2020, four detonation craters were investigated with a multibeam echosounder for the first time. This dataset is represented here as maps of bathymetry, slope angle, and height difference to the surrounding. The circular craters were still clearly visible a year after the detonation. The diameter and depth of the structures were between 7.5–12.6 m and 0.7–2.2 m, respectively. In total, about 321 m2 of the seafloor was destroyed along the track line.

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: 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: Ectoenzymatic activity, prokaryotic heterotrophic abundances and production were determined in the Mediterranean Sea. Sampling was carried out in the sub-surface, the deep chlorophyll maximum layer (DCM), the core of the Levantine intermediate waters and in the deeper part of the mesopelagic layers. Michaelis–Menten kinetics were assessed using a large range of concentrations of fluorogenic substrates (0.025 to 50 µM). As a consequence, Km (Michaelis–Menten half-saturation constant) and Vm (maximum hydrolysis velocity) parameters were determined for both low- and high-affinity enzymes for alkaline phosphatase, aminopeptidase (LAP) and β-glucosidase (βGLU). Based on the constant derived from the high-LAP-affinity enzyme (0.025–1 µM substrate concentration range), in situ hydrolysis of N proteins contributed 48 % ± 30 % to the heterotrophic bacterial nitrogen demand within the epipelagic layers and 180 % ± 154 % in the Levantine intermediate waters and the upper part of the mesopelagic layers. The LAP hydrolysis rate was higher than bacterial N demand only within the deeper layer and only when considering the high-affinity enzyme. Based on a 10 % bacterial growth efficiency, the cumulative hydrolysis rates of C proteins and C polysaccharides contributed on average 2.5 % ± 1.3  % to the heterotrophic bacterial carbon demand in the epipelagic layers sampled (sub-surface and DCM). This study clearly reveals potential biases in current and past interpretations of the kinetic parameters for the three enzymes tested based on the fluorogenic-substrate concentration used. In particular, the LAP / βGLU enzymatic ratios and some of the depth-related trends differed between the use of high and low concentrations of fluorogenic substrates.

Description: It is widely accepted that the commensal gut microbiota contributes to the health and well-being of its host. The solitary tunicate Ciona intestinalis emerges as a model organism for studying host–microbe interactions taking place in the gut, however, the potential of its gut-associated microbiota for marine biodiscovery remains unexploited. In this study, we set out to investigate the diversity, chemical space, and pharmacological potential of the gut-associated microbiota of C. intestinalis collected from the Baltic and North Seas. In a culture-based approach, we isolated 61 bacterial and 40 fungal strains affiliated to 33 different microbial genera, indicating a rich and diverse gut microbiota dominated by Gammaproteobacteria. In vitro screening of the crude microbial extracts indicated their antibacterial (64% of extracts), anticancer (22%), and/or antifungal (11%) potential. Nine microbial crude extracts were prioritized for in-depth metabolome mining by a bioactivity- and chemical diversity-based selection procedure. UPLC-MS/MS-based metabolomics combining automated (feature-based molecular networking and in silico dereplication) and manual approaches significantly improved the annotation rates. A high chemical diversity was detected where peptides and polyketides were the predominant classes. Many compounds remained unknown, including two putatively novel lipopeptides produced by a Trichoderma sp. strain. This is the first study assessing the chemical and pharmacological profile of the cultivable gut microbiota of C. intestinalis

Description: An oxygen decrease of the intermediate-depth low-oxygen zones (300 to 700 m) is seen in time series for selected tropical areas for the period 1960 to 2008 in the eastern tropical Atlantic, the equatorial Pacific and the eastern tropical Indian Ocean. These nearly 5-decade time series were extended to 68 years by including rare historic data starting in 1950 and more recent data. For the extended time series between 1950 and 2018, the deoxygenation trend for the layer 300 to 700 m is similar to the deoxygenation trend seen in the shorter time series. Additionally, temperature, salinity, and nutrient time series in the upper-ocean layer (50 to 300 m) of these areas were investigated since this layer provides critical pelagic habitat for biological communities. Due to the low amount of data available, the results are often not statistically significant within the 95 % confidence interval but nevertheless indicate trends worth discussing. Generally, oxygen is decreasing in the 50 to 300 m layer, except for an area in the eastern tropical South Atlantic. Nutrients also showed long-term trends in the 50 to 300 m layer in all ocean basins and indicate overlying variability related to climate modes. Nitrate increased in all areas. Phosphate also increased in the Atlantic Ocean and Indian Ocean areas, while it decreased in the two areas of the equatorial Pacific Ocean. Silicate decreased in the Atlantic and Pacific areas but increased in the eastern Indian Ocean. Hence, oxygen and nutrients show trends in the tropical oceans, though nutrients trends are more variable between ocean areas than the oxygen trends; therefore, we conclude that those trends are more dependent on local drivers in addition to a global trend. Different positive and negative trends in temperature, salinity, oxygen and nutrients indicate that oxygen and nutrient trends cannot be completely explained by local warming.

Description: We resolve a previously unrecognized shallow subducting seamount from a re-processed multichannel seismic depth image crossing the 1994 M7.8 Java tsunami earthquake slip area. Seamount subduction is related to the uplift of the overriding plate by lateral shortening and vertical thickening, causing pronounced back-thrusting at the landward slope of the forearc high and the formation of splay faults branching off the landward flank of the subducting seamount. The location of the seamount in relation to the 1994 earthquake hypocentre and its co-seismic slip model suggests that the seamount acted as a seismic barrier to the up-dip co-seismic rupture propagation of this moderate size earthquake. The wrapping of the co-seismic slip contours around the seamount indicates that it diverted rupture propagation, documenting the control of forearc structures on seismic rupture.

Description: The main objective of the present study was to assess the effects of sulfur (S) nutrition on plant growth, overall quality, secondary metabolites, and antibacterial and radical scavenging activities of hydroponically grown lettuce cultivars. Three lettuce cultivars, namely, Pazmanea RZ (green butterhead, V1), Hawking RZ (green multi-leaf lettuce, V2), and Barlach RZ (red multi-leaf, V3) were subjected to two S-treatments in the form of magnesium sulfate (+S) or magnesium chloride (−S). Significant differences were observed under −S treatments, especially among V1 and V2 lettuce cultivars. These responses were reflected in the yield, levels of macro- and micro-nutrients, water-soluble sugars, and free inorganic anions. In comparison with the green cultivars (V1 and V2), the red-V3 cultivar revealed a greater acclimation to S starvation, as evidenced by relative higher plant growth. In contrast, the green cultivars showed higher capabilities in production and superior quality attributes under +S condition. As for secondary metabolites, sixteen compounds (e.g., sesquiterpene lactones, caffeoyl derivatives, caffeic acid hexose, 5-caffeoylquinic acid (5-OCQA), quercetin and luteolin glucoside derivatives) were annotated in all three cultivars with the aid of HPLC-DAD-MS-based untargeted metabolomics. Sesquiterpene lactone lactucin and anthocyanin cyanidin 3-O-galactoside were only detected in V1 and V3 cultivars, respectively. Based on the analyses, the V3 cultivar was the most potent radical scavenger, while V1 and V2 cultivars exhibited antibacterial activity against Staphylococcus aureus in response to S provision. Our study emphasizes the critical role of S nutrition in plant growth, acclimation, and nutritional quality. The judicious-S application can be adopted as a promising antimicrobial prototype for medical applications.

Description: Six-degree-of-freedom (6DoF) sensors measure translation along three axes and rotation around three axes. These collocated measurements make it possible to fully describe building motion without the need for an external reference point. This is an advantage for building health monitoring, which uses interstory drift and building eigenfrequencies to monitor stability. In this paper, IMU50 6DoF sensors are characterized to determine their suitability for building health monitoring. The sensors are calibrated using step table methods and by comparison with earth’s rotation and gravity. These methods are found to be comparable. The sensor’s self-noise is examined through the power spectral density and the Allan deviation of data recorded in a quiet environment. The effect of temperature variation is tested between 14 and 50 ∘C. It appears that the self-noise of the rotation components increases while the self-noise of the acceleration components decreases with temperature. The comparison of the sensor self-noise with ambient building signal and higher amplitude shaking shows that these sensors are in general not sensitive enough for ambient signal building health monitoring in the frequency domain, but could be useful for monitoring interstory drift and building motion during, for example, strong earthquake shaking in buildings similar to those examined here.

Description: The Persian Gulf is a unique and biologically diverse marine environment dominated by invertebrates. In continuation of our research interest in the chemistry and biological activity of marine sponges from the Persian Gulf, we selected the excavating sponge Cliona celata for detailed metabolome analyses, in vitro bioactivity screening, and chemical isolation studies. A UPLC-MS/MS (MS2) molecular-networking-based dereplication strategy allowed annotation and structural prediction of various diketopiperazines (DKPs) and etzionin-type diketopiperazine hydroxamates (DKPHs) in the crude sponge extract. The molecular-networking-guided isolation approach applied to the crude extract afforded the DKPH etzionin (1) and its two new derivatives, clioetzionin A (2) and clioetzionin B (3). Another new modified DKP (4) was identified by MS/MS analyses but could not be isolated in sufficient quantities to confirm its structure. The chemical characterization of the purified DKPHs 1–3 was performed by a combination of 1D and 2D NMR spectroscopy, HRMS, HRMS/MS, and [α]D analyses. Compounds 1 and 2 exhibited broad antibacterial, antifungal, and anticancer activities, with IC50 values ranging from 19.6 to 159.1 µM. This is the first study investigating the chemical constituents of a C. celata specimen from the Persian Gulf. It is also the first report of full spectroscopic data of etzionin based on extensive spectroscopic analyses

Description: During the early part of the last glacial termination (17.2-15 ka) and coincident with a ∼ 35 ppm rise in atmospheric CO2, a sharp 0.3‰-0.4‰ decline in atmospheric δ13CO2 occurred, potentially constraining the key processes that account for the early deglacial CO2 rise. A comparable δ13C decline has also been documented in numerous marine proxy records from surface and thermocline-dwelling planktic foraminifera. The δ13C decline recorded in planktic foraminifera has previously been attributed to the release of respired carbon from the deep ocean that was subsequently transported within the upper ocean to sites where the signal was recorded (and then ultimately transferred to the atmosphere). Benthic δ13C records from the global upper ocean, including a new record presented here from the tropical Pacific, also document this distinct early deglacial δ13C decline. Here we present modeling evidence to show that rather than respired carbon from the deep ocean propagating directly to the upper ocean prior to reaching the atmosphere, the carbon would have first upwelled to the surface in the Southern Ocean where it would have entered the atmosphere. In this way the transmission of isotopically light carbon to the global upper ocean was analogous to the ongoing ocean invasion of fossil fuel CO2. The model results suggest that thermocline waters throughout the ocean and 500-2000m water depths were affected by this atmospheric bridge during the early deglaciation.

Description: Ocean warming and acidification pose serious threats to cold-water corals (CWCs) and the surrounding habitat. Yet, little is known about the role of natural short-term and seasonal environmental variability, which could be pivotal to determine the resilience of CWCs in a changing environment. Here, we provide continuous observational data of the hydrodynamic regime (recorded using two benthic landers) and point measurements of the carbonate and nutrient systems from five Lophelia pertusa reefs in the Langenuen Fjord, southwestern Norway, from 2016 to 2017. In this fjord setting, we found that over a tidal (〈24 h) cycle during winter storms, the variability of measured parameters at CWC depths was comparable to the intra-annual variability, demonstrating that single point measurements are not sufficient for documenting (and monitoring) the biogeochemical conditions at CWC sites. Due to seasonal and diurnal forcing, parts of the reefs experienced temperatures up to 4 °C warmer (i.e., 〉12 °C) than the mean conditions and high CT concentrations of 20 µmol kg−1 over the suggested threshold for healthy CWC reefs (i.e., 〉2170 µmol kg−1). Combined with hindcast measurements, our findings indicate that these shallow fjord reefs may act as an early hotspot for ocean warming and acidification. We predict that corals in Langenuen will face seasonally high temperatures (〉18 °C) and hypoxic and corrosive conditions within this century. Therefore, these fjord coral communities could forewarn us of the coming consequences of climate change on CWC diversity and function.

Description: The northern margin of the Ligurian Basin shows notable seismicity at the Alpine front, including frequent magnitude 4 events. Seismicity decreases offshore towards the Basin centre and Corsica, revealing a diffuse distribution of low magnitude earthquakes. We analyse data of the amphibious AlpArray seismic network with focus on the offshore component, the AlpArray OBS network, consisting of 24 broadband ocean bottom seismometers deployed for eight months, to reveal the seismicity and depth distribution of micro-earthquakes beneath the Ligurian Sea. Two clusters occurred between ~10 km to ~16 km depth below sea surface, within the lower crust and uppermost mantle. Thrust faulting focal mechanisms indicate compression and an inversion of the Ligurian Basin, which is an abandoned Oligocene rift basin. The Basin inversion is suggested to be related to the Africa-Europe plate convergence. The locations and focal mechanisms of seismicity suggest reactivation of pre-existing rift structures. Slightly different striking directions of faults in the basin centre compared to faults further east and hence away from the abandoned rift may mimic the counter-clockwise rotation of the Corsica-Sardinia block during ~20–16 Ma. The observed cluster events support the hypothesis of strengthening of crust and uppermost mantle during rifting related extension and thinning of continental crust.

Description: Previous studies with Baltic Sea phytoplankton combining elevated seawater temperature with CO2 revealed the importance of size trait‐based analyses, in particular dividing the plankton into edible (〉5 and 〈100 μm) and inedible ([removed]100 μm) size classes for mesozoopankton grazers. While the edible phytoplankton responded predominantly negative to warming and the inedible group stayed unaffected or increased, independent from edibility most phytoplankton groups gained from CO2. Because the ratio between edible and inedible taxa changes profoundly over seasons, we investigated if community responses can be predicted according to the prevailing composition of edible and inedible groups. We experimentally explored the combined effects of elevated temperatures and CO2 concentrations on a late‐summer Baltic Sea community. Total phytoplankton significantly increased in response to elevated CO2 in particular in combination with temperature, driven by a significant gain of the inedible 〈5 μm fraction and large filamentous cyanobacteria. Large flagellates disappeared. The edible group was low as usual in summer and decreased with both factors due to enhanced copepod grazing and overall decline of small flagellates. Our results emphasize that the responses of summer communities are complex, but can be predicted by the composition and dominance of size classes and groups.

Description: Oceans play a major role on the exchange of carbon with the atmosphere and thereby on past climates with glacial/interglacial variations of the CO2 concentration. The melting of ice sheets during deglaciations lets the sea level rise which leads to the flooding of coastal land areas resulting in the transfer of terrestrial organic matter to the ocean. However, the consequences of such fluxes on the ocean biogeochemical cycle and uptake/release of CO2 are poorly constrained. Moreover, this potentially important exchange of carbon at the land-sea interface is not represented in most Earth System Models. We present here the implementation of terrestrial organic matter fluxes into the ocean at the transiently changing land-sea interface in the Max Planck Institute for Meteorology Earth System Model (MPI-ESM) and investigate their effect on the biogeochemistry during the last deglaciation. Our results show that during the deglaciation, most of the terrestrial organic matter inputs to the ocean occurs during Meltwater Pulse 1a (between 15–14 ka) which leads to additional 21.2 GtC of terrestrial origin (mostly originating from wood and humus). Although this additional organic matter input is relatively small in comparison to the global ocean inventory (0.06 %) and thus doesn’t have an impact on the global CO2 flux, the terrestrial organic matter fluxes initiate oceanic outgassing at regional hotspots like in Indonesia for a few hundred years. Finally, sensitivity experiments highlight that terrestrial organic matter fluxes are the drivers of oceanic outgassing in flooded coastal regions during Meltwater Pulse 1a. Furthermore, the magnitude of outgassing is rather insensitive to higher carbon to nutrients ratios of the terrestrial organic matter. Our results provide a first estimate of the importance of terrestrial organic matter fluxes in a transient deglaciation simulation. Moreover, our model development is an important step towards a fully coupled carbon cycle in an Earth System Model applicable for simulations of glacial/interglacial cycles.

Description: In this study, we analyzed a large seismological dataset from temporary and permanent networks in the southern and eastern Alps to establish high-precision hypocenters and 1-D VP and VP/VS models. The waveform data of a subset of local earthquakes with magnitudes in the range of 1–4.2 ML were recorded by the dense, temporary SWATH-D network and selected stations of the AlpArray network between September 2017 and the end of 2018. The first arrival times of P and S waves of earthquakes are determined by a semi-automatic procedure. We applied a Markov chain Monte Carlo inversion method to simultaneously calculate robust hypocenters, a 1-D velocity model, and station corrections without prior assumptions, such as initial velocity models or earthquake locations. A further advantage of this method is the derivation of the model parameter uncertainties and noise levels of the data. The precision estimates of the localization procedure is checked by inverting a synthetic travel time dataset from a complex 3-D velocity model and by using the real stations and earthquakes geometry. The location accuracy is further investigated by a quarry blast test. The average uncertainties of the locations of the earthquakes are below 500 m in their epicenter and ∼ 1.7 km in depth. The earthquake distribution reveals seismicity in the upper crust (0–20 km), which is characterized by pronounced clusters along the Alpine frontal thrust, e.g., the Friuli-Venetia (FV) region, the Giudicarie–Lessini (GL) and Schio-Vicenza domains, the Austroalpine nappes, and the Inntal area. Some seismicity also occurs along the Periadriatic Fault. The general pattern of seismicity reflects head-on convergence of the Adriatic indenter with the Alpine orogenic crust. The seismicity in the FV and GL regions is deeper than the modeled frontal thrusts, which we interpret as indication for southward propagation of the southern Alpine deformation front (blind thrusts).

Description: The recent development of ambient noise tomography, in combination with the increasing number of permanent seismic stations and dense networks of temporary stations operated during passive seismic experiments, provides a unique opportunity to build the first high-resolution 3-D shear wave velocity (vS) model of the entire crust of the Bohemian Massif (BM). This paper provides a regional-scale model of velocity distribution in the BM crust. The velocity model with a cell size of 22 km is built using a conventional two-step inversion approach from Rayleigh wave group velocity dispersion curves measured at more than 400 stations. The shear velocities within the upper crust of the BM are ∼0.2 km s−1 higher than those in its surroundings. The highest crustal velocities appear in its southern part, the Moldanubian unit. The Cadomian part of the region has a thinner crust, whereas the crust assembled, or tectonically transformed in the Variscan period, is thicker. The sharp Moho discontinuity preserves traces of its dynamic development expressed in remnants of Variscan subductions imprinted in bands of crustal thickening. A significant feature of the presented model is the velocity-drop interface (VDI) modelled in the lower part of the crust. We explain this feature by the anisotropic fabric of the lower crust, which is characterised as vertical transverse isotropy with the low velocity being the symmetry axis. The VDI is often interrupted around the boundaries of the crustal units, usually above locally increased velocities in the lowermost crust. Due to the north-west–south-east shortening of the crust and the late-Variscan strike-slip movements along the north-east–south-west oriented sutures preserved in the BM lithosphere, the anisotropic fabric of the lower crust was partly or fully erased along the boundaries of original microplates. These weakened zones accompanied by a velocity increase above the Moho (which indicate an emplacement of mantle rocks into the lower crust) can represent channels through which portions of subducted and later molten rocks have percolated upwards providing magma to subsequently form granitoid plutons.

Description: Accurate subsurface velocity models are crucial for geological interpretations based on seismic depth images. Seismic reflection tomography is an effective iterative method to update and refine a preliminary velocity model for depth imaging. Based on residual move-out analysis of reflectors in common image point gathers an update of the velocity is estimated by a ray-based tomography. To stabilize the tomography, several preconditioning strategies exist. Most critical is the estimation of the depth error to account for the residual move-out of the reflector in the common image point gathers. Because the depth errors for many closely spaced image gathers must be picked, manual picking is extremely time-consuming, human biased, and not reproducible. Data-driven picking algorithms based on coherence or semblance analysis are widely used for hyperbolic or linear events. However, for complex-shaped depth events, pure data-driven picking is difficult. To overcome this, the warping method named Non-Rigid Matching is used to estimate a depth error displacement field. Warping is used, e.g., to merge photographic images or to match two seismic images from time-lapse data. By calculating the displacements between an offset to its neighbouring offset in the common image point domain, a locally smooth-shaped displacement field is defined for each data sample. Depending on the complexity of the subsurface, sample tracking through the displacement field along predefined horizons or on a simple regular grid yields discrete depth error values for the tomography. The application to a multi-channel seismic line across the Sunda subduction zone offshore Lombok island, Indonesia, illustrates the approach and documents the advantages of the method to estimate a detailed velocity structure in a complex tectonic regime. By incorporating the warping scheme into the reflection tomography, we demonstrate an increase in the velocity resolution and precision by improving the data-driven accuracy of depth error picks with arbitrary shapes. This approach offers the possibility to use the full capacities of tomography and further leads to more accurate interpretations of complex geological structures.

Description: Direct comparison between paleo oceanic δ13C records and model results facilitates assessing simulated distributions and properties of water masses in the past. To accomplish this, we include a new representation of the stable carbon isotope 13C into the HAMburg Ocean Carbon Cycle model (HAMOCC), the ocean biogeochemical component of the Max Planck Institute Earth System Model (MPI-ESM). 13C is explicitly resolved for all existing oceanic carbon pools. We account for fractionation during air-sea gas exchange and for biological fractionation εp associated with photosynthetic carbon fixation during phytoplankton growth. We examine two εp parameterisations of different complexity: εpPopp varies with surface dissolved CO2 concentration (Popp et al., 1989), while εpLaws additionally depends on local phytoplankton growth rates (Laws et al., 1995). When compared to observations of δ13C in dissolved inorganic carbon (DIC), both parameterisations yield similar performance. However, with regard to δ13C in particulate organic carbon εpPopp shows a considerably improved performance than εpLaws, because the latter results in a too strong preference for 12C. The model also well reproduces the oceanic 13C Suess effect, i.e. the intrusion of the isotopically light anthropogenic CO2 into the ocean, based on comparison to other existing 13C models and to observation-based oceanic carbon uptake estimates over the industrial period. We further apply the approach of Eide et al. (2017a), who derived the first global oceanic 13C Suess effect estimate based on observations, to our model data that has ample spatial and temporal coverage. With this we are able to analyse in detail the underestimation of 13C Suess effect by this approach as it has been noted by Eide et al. (2017a). Based on our model we find underestimations of 13C Suess effect at 200 m by 0.24 ‰ in the Indian Ocean, 0.21 ‰ in the North Pacific, 0.26 ‰ in the South Pacific, 0.1 ‰ in the North Atlantic and 0.14 ‰ in the South Atlantic. We attribute the major sources of the underestimation to two assumptions in Eide et al. (2017a)'s approach: a spatially-constant preformed component of δ13CDIC in year 1940 and neglecting 13C Suess effect in CFC-12 free water.

Description: A hierarchy of global 1/4° (ORCA025) and Atlantic Ocean 1/20° nested (VIKING20X) ocean/sea-ice models is described. It is shown that the eddy-rich configurations performed in hindcasts of the past 50–60 years under CORE and JRA55-do atmospheric forcings realistically simulate the large-scale horizontal circulation, the distribution of the mesoscale, overflow and convective processes, and the representation of regional current systems in the North and South Atlantic. The representation, and in particular the long-term temporal evolution, of the Atlantic Meridional Overturning Circulation (AMOC) strongly depends on numerical choices for the application of freshwater fluxes. The interannual variability of the AMOC instead is highly correlated among the model experiments and also with observations, including the 2010 minimum observed by RAPID at 26.5° N pointing at a dominant role of the forcing. Regional observations in western boundary current systems at 53° N, 26.5° N and 11° S are explored in respect to their ability to represent the AMOC and to monitor the temporal evolution of the AMOC. Apart from the basin-scale measurements at 26.5° N, it is shown that in particular the outflow of North Atlantic Deepwater at 53° N is a good indicator of the subpolar AMOC trend during the recent decades, if the latter is provided in density coordinates. The good reproduction of observed AMOC and WBC trends in the most reasonable simulations indicate that the eddy-rich VIKING20X is capable in representing realistic forcing-related and ocean-intrinsic trends.

Description: Current trends in the global climate facilitate the displacement of numerous marine species from their native distribution ranges to higher latitudes when facing warming conditions. In this work, we analyzed occurrences of a circumtropical reef fish, the spotfin burrfish, Chilomycterus reticulatus (Linnaeus, 1958), in the Madeira Archipelago (NE Atlantic) between 1898 and 2021. In addition to available data sources, we performed an online survey to assess the distribution and presence of this species in the Madeira Archipelago, along with other relevant information, such as size class and year of the first sighting. In total, 28 valid participants responded to the online survey, georeferencing 119 C. reticulatus sightings and confirming its presence in all archipelago islands. The invasiveness of the species was screened using the Aquatic Species Invasiveness Screening Kit. Five assessments rated the fish as being of medium risk of establishing a local population and becoming invasive. Current temperature trends might have facilitated multiple sightings of this thermophilic species in the Madeira Archipelago. The present study indicates an increase in C. reticulatus sightings in the region. This underlines the need for updated comprehensive information on species diversity and distribution to support informed management and decisions. The spread of yet another thermophilic species in Madeiran waters provides further evidence of an ongoing tropicalization, emphasizing the need for monitoring programs and the potential of citizen science in complementing such programs.

Description: This paper presents a modelling study on the fate of CHBr3 and its product gases in the troposphere within the context of tropical deep convection. A cloud-scale case study was conducted along the west coast of Borneo, where several deep convective systems were triggered on the afternoon and early evening of 19 November 2011. These systems were sampled by the Falcon aircraft during the field campaign of the SHIVA project and analysed using a simulation with the cloud-resolving meteorological model C-CATT-BRAMS at 2x2 km resolution that represents the emissions, transport by large-scale flow, convection, photochemistry, and washout of CHBr3 and its product gases (PGs). We find that simulated CHBr3 mixing ratios and the observed values in the boundary layer and the outflow of the convective systems agree. However, the model underestimates the background CHBr3 mixing ratios in the upper troposphere, which suggests a missing source at the regional scale. An analysis of the simulated chemical speciation of bromine within and around each simulated convective system during the mature convective stage reveals that 〉 85% of the bromine derived from CHBr3 and its PGs is transported vertically to the point of convective detrainment in the form of CHBr3 and that the remaining small fraction is in the form of organic PGs, principally insoluble brominated carbonyls produced from the photo-oxidation of CHBr3. The model simulates that within the boundary layer and free troposphere, the inorganic PGs are only present in soluble forms, i.e. HBr, HOBr, and BrONO2, and, consequently, within the convective clouds, the inorganic PGs are almost entirely removed by wet scavenging. We find that HBr is the most abundant PG in background lower-tropospheric air and that this prevalence of HBr is a result of the relatively low background tropospheric ozone levels at the regional scale. Contrary to a previous study in a different environment, for the conditions in the simulation, the insoluble Br-2 species is hardly formed within the convective systems and therefore plays no significant role in the vertical transport of bromine. This likely results from the relatively small quantities of simulated inorganic bromine involved, the presence of HBr in large excess compared to HOBr and BrO, and the relatively efficient removal of soluble compounds within the convective column.

Description: Improved quantification techniques of natural sources is needed to explain variations in atmospheric methane. In polar regions, high uncertainties in current estimates of methane release from the seabed remain. We present two unique 10 and 3 months long time-series of bottom water measurements of physical and chemical parameters from two autonomous ocean observatories deployed at separate intense seabed methane seep sites (91 and 246 m depth) offshore Western Svalbard from 2015 to 2016. Results show high short term (100–1000 nmol L-1 within hours) and seasonal variation, as well as higher (2–7 times) methane concentrations compared to previous measurements. Rapid variability is explained by uneven distribution of seepage and changing ocean current directions. No overt influence of tidal hydrostatic pressure or water temperature variations on methane concentration was observed, but an observed negative correlation with temperature at the 246 site fits with hypothesized seasonal blocking of lateral methane pathways in the sediments. Negative correlation between bottom water methane concentration/variability and wind forcing, concomitant with signs of weaker water column stratification, indicates increased potential for methane release to the atmosphere in fall/winter. We highlight uncertainties in methane inventory estimates based on discrete water sampling and present new information about short- and long-term methane variability which can help constrain future estimates of seabed methane seepage.

Description: Earth system and climate modelling involves the simulation of processes on a wide range of scales and within and across various compartments of the Earth system. In practice, component models are often developed independently by different research groups, adapted by others to their special interests and then combined using a dedicated coupling software. This procedure not only leads to a strongly growing number of available versions of model components and coupled setups but also to model- and high-performance computing (HPC)-system-dependent ways of obtaining, configuring, building and operating them. Therefore, implementing these Earth system models (ESMs) can be challenging and extremely time consuming, especially for less experienced modellers or scientists aiming to use different ESMs as in the case of intercomparison projects. To assist researchers and modellers by reducing avoidable complexity, we developed the ESM-Tools software, which provides a standard way for downloading, configuring, compiling, running and monitoring different models on a variety of HPC systems. It should be noted that ESM-Tools is not a coupling software itself but a workflow and infrastructure management tool to provide access to increase usability of already existing components and coupled setups. As coupled ESMs are technically the more challenging tasks, we will focus on coupled setups, always implying that stand-alone models can benefit in the same way. With ESM-Tools, the user is only required to provide a short script consisting of only the experiment-specific definitions, while the software executes all the phases of a simulation in the correct order. The software, which is well documented and easy to install and use, currently supports four ocean models, three atmosphere models, two biogeochemistry models, an ice sheet model, an isostatic adjustment model, a hydrology model and a land-surface model. Compared to previous versions, ESM-Tools has lately been entirely recoded in a high-level programming language (Python) and provides researchers with an even more user-friendly interface for Earth system modelling. ESM-Tools was developed within the framework of the Advanced Earth System Model Capacity project, supported by the Helmholtz Association.

Description: Ocean science and hydroacoustic seafloor mapping rely on accurate navigation underwater. By exploiting terrain information provided by a multibeam echosounder system, it is possible to significantly improve map quality. This article presents an algorithm capable of improving map quality and accuracy by aligning consecutive pings to tiles that are matched pairwise. A globally consistent solution is calculated from these matches. The proposed method has the potential to be used online in addition to other navigation solutions, but is mainly targeted for post processing. The algorithm was tested using different parameter settings on an AUV and a ship-based dataset. The ship-based dataset is publicly available as a benchmark. The original accurate navigation serving as a ground truth, alongside trajectories that include an artificial drift, are available. This allows quantitative comparisons between algorithms and parameter settings.

Description: The past and future evolution of the Antarctic Ice Sheet is largely controlled by interactions between the ocean and floating ice shelves. To investigate these interactions, coupled ocean and ice sheet model configurations are required. Previous modelling studies have mostly relied on high-resolution configurations, limiting these studies to individual glaciers or regions over short timescales of decades to a few centuries. We present a framework to couple the dynamic ice sheet model PISM (Parallel Ice Sheet Model) with the global ocean general circulation model MOM5 (Modular Ocean Model) via the ice shelf cavity model PICO (Potsdam Ice-shelf Cavity mOdel). As ice shelf cavities are not resolved by MOM5 but are parameterized with the PICO box model, the framework allows the ice sheet and ocean components to be run at resolutions of 16 km and 3∘ respectively. This approach makes the coupled configuration a useful tool for the analysis of interactions between the Antarctic Ice Sheet and the global ocean over time spans of the order of centuries to millennia. In this study, we describe the technical implementation of this coupling framework: sub-shelf melting in the ice sheet component is calculated by PICO from modelled ocean temperatures and salinities at the depth of the continental shelf, and, vice versa, the resulting mass and energy fluxes from melting at the ice–ocean interface are transferred to the ocean component. Mass and energy fluxes are shown to be conserved to machine precision across the considered component domains. The implementation is computationally efficient as it introduces only minimal overhead. Furthermore, the coupled model is evaluated in a 4000 year simulation under constant present-day climate forcing and is found to be stable with respect to the ocean and ice sheet spin-up states. The framework deals with heterogeneous spatial grid geometries, varying grid resolutions, and timescales between the ice and ocean component in a generic way; thus, it can be adopted to a wide range of model set-ups.

Description: The skill of global ocean biogeochemical models, and the earth system models in which they are embedded, can be improved by systematic calibration of the parameter values against observations. However, such tuning is seldom undertaken as these models are computationally very expensive. Here we investigate the performance of DFO-LS, a local, derivative-free optimisation algorithm which has been designed for computationally expensive models with irregular model–data misfit landscapes typical of biogeochemical models. We use DFO-LS to calibrate six parameters of a relatively complex global ocean biogeochemical model (MOPS) against synthetic dissolved oxygen, phosphate and nitrate “observations” from a reference run of the same model with a known parameter configuration. The performance of DFO-LS is compared with that of CMA-ES, another derivative-free algorithm that was applied in a previous study to the same model in one of the first successful attempts at calibrating a global model of this complexity. We find that DFO-LS successfully recovers five of the six parameters in approximately 40 evaluations of the misfit function (each one requiring a 3000-year run of MOPS to equilibrium), while CMA-ES needs over 1200 evaluations. Moreover, DFO-LS reached a “baseline” misfit, defined by observational noise, in just 11–14 evaluations, whereas CMA-ES required approximately 340 evaluations. We also find that the performance of DFO-LS is not significantly affected by observational sparsity, however fewer parameters were successfully optimised in the presence of observational uncertainty. The results presented here suggest that DFO-LS is sufficiently inexpensive and robust to apply to the calibration of complex, global ocean biogeochemical models.

Description: An increasing number of dead zoning (hypoxia) has been reported as a consequence of declining levels of dissolved oxygen in coastal oceans all over the globe. Despite substantial efforts a quantitative description of hypoxia up to a level enabling reliable predictions has not been achieved yet for most regions of societal interest. This does also apply to Eckernförde Bight (EB) situated in the Baltic Sea, Germany. The aim of this study is to dissect underlying mechanisms of hypoxia in EB, to identify key sources of uncertainties, and to explore the potential of existing monitoring programs to predict hypoxia by developing and documenting a workflow that may be applicable to other regions facing similar challenges. Our main tool is an ultra-high spatially resolved general ocean circulation model based on a code framework of proven versatility in that it has been applied to various regional and even global simulations in the past. Our model configuration features a spacial horizontal resolution of 100 m (unprecedented in the underlying framework which is used in both global and regional applications) and includes an elementary representation of the biogeochemical dynamics of dissolved oxygen. In addition, we integrate artificial “clocks” that measure the residence time of the water in EB along with timescales of (surface) ventilation. Our approach relies on an ensemble of hindcast model simulations, covering the period from 2000 to 2018, designed to cover a range of poorly known model parameters for vertical background mixing (diffusivity) and local oxygen consumption within EB. Feed-forward artificial neural networks are used to identify predictors of hypoxia deep in EB based on data at a monitoring site at the entrance of EB. Our results consistently show that the dynamics of low (hypoxic) oxygen concentrations in bottom waters deep inside EB is, to first order, determined by the following antagonistic processes: (1) the inflow of low-oxygenated water from the Kiel Bight (KB) – especially from July to October – and (2) the local ventilation of bottom waters by local (within EB) subduction and vertical mixing. Biogeochemical processes that consume oxygen locally are apparently of minor importance for the development of hypoxic events. Reverse reasoning suggests that subduction and mixing processes in EB contribute, under certain environmental conditions, to the ventilation of the KB by exporting recently ventilated waters enriched in oxygen. A detailed analysis of the 2017 fish-kill incident highlights the interplay between westerly winds importing hypoxia from KB and ventilating easterly winds which subduct oxygenated water.

Description: Understanding trophic interactions is essential for the prediction and measurement of structure and function in aquatic environments. Communities in these ecosystems may be shaped by variables such as predator diversity, prey density and emergent multiple predator effects (MPEs), which are likely to influence trophic dynamics. In this study, we examined the effect of key predatory fish in floodplain wetlands, namely Oreochromis mossambicus and Enteromius paludinosus, towards Chironomidae prey, using a comparative functional response (FR) approach. We used single predator species as well as intra- and interspecific paired species to contrast FRs under multiple predator scenarios. Attack rate and handling time estimates from single predator FRs were used to predict multiple predators’ feeding rates, which were compared to observe multiple predators’ feeding rates to quantify potential MPEs. From single fish trials, each species displayed a significant Type II FR, characterized by high feeding rates at low prey densities. Oreochromis mossambicus had a steeper (initial slope, i.e., higher attack rate) and higher (asymptote of curve, i.e., shorter handling time and higher maximum feeding rate) FR, whereas E. paludinosus exhibited lower-magnitude FRs (i.e., lower attack rate, longer handling time and lower feeding rate). In multiple predator scenarios, feeding rates were well-predicted by those of single predators, both in conspecific and interspecific pairs, and thus we did not find evidence for antagonistic or synergistic MPEs. Predator–prey interactions in wetland systems can have significant consequences on the structure and dynamics of ecological communities. In turn, this could have destabilizing effects on resources in tropical wetlands. These results, although experimental, help us understand how trophic interaction among conspecific or interspecific fish species in Austral tropical wetlands might influence their aquatic prey species. This will help us to understand food web dynamics better.

Description: Seaweed endophytic (algicolous) fungi are talented producers of bioactive natural products. We have previously isolated two strains of the endophytic fungus, Pyrenochaetopsis sp. FVE-001 and FVE-087, from the thalli of the brown alga Fucus vesiculosus. Initial chemical studies yielded four new decalinoylspirotetramic acid derivatives with antimelanoma activity, namely pyrenosetins A–C (1–3) from Pyrenochaetopsis sp. strain FVE-001, and pyrenosetin D (4) from strain FVE-087. In this study, we applied a comparative metabolomics study employing HRMS/MS based feature-based molecular networking (FB MN) on both Pyrenochaetopsis strains. A higher chemical capacity in production of decalin derivatives was observed in Pyrenochaetopsis sp. FVE-087. Notably, several decalins showed different retention times despite the same MS data and MS/MS fragmentation pattern with the previously isolated pyrenosetins, indicating they may be their stereoisomers. FB MN-based targeted isolation studies coupled with antimelanoma activity testing on the strain FVE-087 afforded two new stereoisomers, pyrenosetins E (5) and F (6). Extensive NMR spectroscopy including DFT computational studies, HR-ESIMS, and Mosher’s ester method were used in the structure elucidation of compounds 5 and 6. The 3′R,5′R stereochemistry determined for compound 6 was identical to that previously reported for pyrenosetin C (3), whose stereochemistry was revised as 3′S,5′R in this study. Pyrenosetin E (5) inhibited the growth of human malignant melanoma cells (A-375) with an IC50 value of 40.9 μM, while 6 was inactive. This study points out significant variations in the chemical repertoire of two closely related fungal strains and the versatility of FB MN in identification and targeted isolation of stereoisomers. It also confirms that the little-known fungal genus Pyrenochaetopsis is a prolific source of complex decalinoylspirotetramic acid derivatives.

Description: Climate-relevant trace gas air-sea exchange exerts an important control on air quality and climate, especially in remote regions of the planet such as the Southern Ocean. It is clear that polar regions exhibit seasonal trends in productivity and biogeochemical cycling, but almost all of the measurements there are skewed to summer months. If we want to understand how the Southern Ocean effects the balance of climate through trace gas air-sea exchange, it is essential to expand our measurement database over greater temporal and spatial scales, including all seasons. Therefore, in this study, we report measured concentrations of dimethylsulphide (DMS, and related sulphur compounds) and isoprene in the Atlantic sector of the Southern Ocean during the winter to understand the spatial and temporal distribution in comparison to current knowledge and climatological calculations for the Southern Ocean. The observations of isoprene are the first in the winter season in the Southern Ocean. We find that concentrations and fluxes of DMS and isoprene in the investigated area are generally lower than those presented or calculated in currently used climatologies and models. More data is urgently needed to better interpolate climatological values and validate process-oriented models, as well as to explore how finer measurement resolution, both spatially and temporally, can influence air-sea flux calculations.

Description: The Leeuwin Current flowing southward along West Australia is an important conduit for the poleward heat transport and interocean water exchange between the tropical and the subantarctic ocean areas. Its past development, and its relationship to Southern Ocean change and to Australian ecosystem response, however is largely unknown. We here reconstruct sea surface and thermocline temperatures and salinities from foraminiferal-based Mg/Ca and stable oxygen isotopes from offshore southwest and southeast Australia reflecting the Leeuwin Current dynamics over the last 60 kyrs. Its variability resembles the biomass burning development in Australasia from ~60–20 ka BP implying that climate-modulated changes related to the Leeuwin Current most likely affected Australian vegetational and fire regimes. In particular during ~60–43 ka BP, warmest thermocline temperatures point to a strongly developed Leeuwin Current during Antarctic cool periods when the Antarctic Circumpolar Current weakened. The pronounced centennial-scale variations in Leeuwin Current strength appear in line with the migrations of the Southern Hemisphere frontal system and are captured by prominent changes in the Australian megafauna biomass. We argue that the concerted action of a rapidly changing Leeuwin Current, the ecosystem response in Australia, and human interference since ~50 BP enhanced the ecological stress on the Australian megafauna until a tipping point was reached at ~43 ka BP, after which faunal recuperation no longer took place. While being weakest during the last glacial maximum, the deglacial Leeuwin Current intensified at times of poleward migrations of the Subtropical Front. During the Holocene, the thermocline off South Australia was considerably shallower compared to the short-term glacial and deglacial periods of Leeuwin Current intensification.

Description: Machine learning spatial modeling is used for mapping the distribution of deep-sea polymetallic nodules (PMN). However, the presence and influence of spatial autocorrelation (SAC) have not been extensively studied. SAC can provide information regarding the variable selection before modeling, and it results in erroneous validation performance when ignored. ML models are also problematic when applied in areas far away from the initial training locations, especially if the (new) area to be predicted covers another feature space. Here, we study the spatial distribution of PMN in a geomorphologically heterogeneous area of the Peru Basin, where SAC of PMN exists. The local Moran’s I analysis showed that there are areas with a significantly higher or lower number of PMN, associated with different backscatter values, aspect orientation, and seafloor geomorphological characteristics. A quantile regression forests (QRF) model is used using three cross-validation (CV) techniques (random-, spatial-, and cluster-blocking). We used the recently proposed “Area of Applicability” method to quantify the geographical areas where feature space extrapolation occurs. The results show that QRF predicts well in morphologically similar areas, with spatial block cross-validation being the least unbiased method. Conversely, random-CV overestimates the prediction performance. Under new conditions, the model transferability is reduced even on local scales, highlighting the need for spatial model-based dissimilarity analysis and transferability assessment in new areas.

Description: We describe and test a new model of biological marine silicate cycling, implemented in the Kiel Marine Biogeochemical Model version 3 (KMBM3), embedded in the University of Victoria Earth System Climate Model (UVic ESCM) version 2.9. This new model adds diatoms, which are a key component of the biological carbon pump, to an existing ecosystem model. This new model combines previously published parameterizations of a diatom functional type, opal production and export with a novel, temperature-dependent dissolution scheme. Modelled steady-state biogeochemical rates, carbon and nutrient distributions are similar to those found in previous model versions. The new model performs well against independent ocean biogeochemical indicators and captures the large-scale features of the marine silica cycle to a degree comparable to similar Earth system models. Furthermore, it is computationally efficient, allowing both fully coupled, long-timescale transient simulations and “offline” transport matrix spinups. We assess the fully coupled model against modern ocean observations, the historical record starting from 1960 and a business-as-usual atmospheric CO2 forcing to the year 2300. The model simulates a global decline in net primary production (NPP) of 1.4 % having occurred since the 1960s, with the strongest declines in the tropics, northern midlatitudes and Southern Ocean. The simulated global decline in NPP reverses after the year 2100 (forced by the extended RCP8.5 CO2 concentration scenario), and NPP returns to 98 % of the pre-industrial rate by 2300. This recovery is dominated by increasing primary production in the Southern Ocean, mostly by calcifying phytoplankton. Large increases in calcifying phytoplankton in the Southern Ocean offset a decline in the low latitudes, producing a global net calcite export in 2300 that varies only slightly from pre-industrial rates. Diatom distribution moves southward in our simulations, following the receding Antarctic ice front, but diatoms are outcompeted by calcifiers across most of their pre-industrial Southern Ocean habitat. Global opal export production thus drops to 75 % of its pre-industrial value by 2300. Model nutrients such as phosphate, silicate and nitrate build up along the Southern Ocean particle export pathway, but dissolved iron (for which ocean sources are held constant) increases in the upper ocean. This different behaviour of iron is attributed to a reduction of low-latitude NPP (and consequently, a reduction in both uptake and export and particle, including calcite scavenging), an increase in seawater temperatures (raising the solubility of particulate iron) and stratification that “traps” the iron near the surface. These results are meant to serve as a baseline for sensitivity assessments to be undertaken with this model in the future.

Description: The region around the town Albstadt, SW Germany, was struck by four damaging earthquakes with magnitudes greater than 5 during the last century. These earthquakes occurred along the Albstadt Shear Zone (ASZ), which is characterized by more or less continuous microseismicity. As there are no visible surface ruptures that may be connected to the fault zone, we study its characteristics by its seismicity distribution and faulting pattern. We use the earthquake data of the state earthquake service of Baden-Württemberg from 2011 to 2018 and complement it with additional phase picks beginning in 2016 at the AlpArray and StressTransfer seismic networks in the vicinity of the ASZ. This extended data set is used to determine new minimum 1-D seismic vp and vs velocity models and corresponding station delay times for earthquake relocation. Fault plane solutions are determined for selected events, and the principal stress directions are derived. The minimum 1-D seismic velocity models have a simple and stable layering with increasing velocity with depth in the upper crust. The corresponding station delay times can be explained well by the lateral depth variation of the crystalline basement. The relocated events align about north–south with most of the seismic activity between the towns of Tübingen and Albstadt, east of the 9∘ E meridian. The events can be separated into several subclusters that indicate a segmentation of the ASZ. The majority of the 25 determined fault plane solutions feature an NNE–SSW strike but NNW–SSE-striking fault planes are also observed. The main fault plane associated with the ASZ dips steeply, and the rake indicates mainly sinistral strike-slip, but we also find minor components of normal and reverse faulting. The determined direction of the maximum horizontal stress of 140–149∘ is in good agreement with prior studies. Down to ca. 7–8 km depth SHmax is bigger than SV; below this depth, SV is the main stress component. The direction of SHmax indicates that the stress field in the area of the ASZ is mainly generated by the regional plate driving forces and the Alpine topography.

Description: Serum calcium isotopes (δ44/42Ca) have been suggested as a non-invasive and sensitive Ca balance marker. Quantitative δ44/42Ca changes associated with Ca flux across body compartment barriers relative to the dietary Ca and the correlation of δ44/42CaSerum with bone histology are unknown. We analyzed Ca and δ44/42Ca by mass-spectrometry in rats after two weeks of standard-Ca-diet (0.5%) and after four subsequent weeks of standard- and of low-Ca-diet (0.25%). In animals on a low-Ca-diet net Ca gain was 61 ± 3% and femur Ca content 68 ± 41% of standard-Ca-diet, bone mineralized area per section area was 68 ± 15% compared to standard-Ca-diet. δ44/42Ca was similar in the diets, and decreased in feces and urine and increased in serum in animals on low-Ca-diet. δ44/42CaBone was higher in animals on low-Ca-diet, lower in the diaphysis than the metaphysis and epiphysis, and unaffected by gender. Independent of diet, δ44/42CaBone was similar in the femora and ribs. At the time of sacrifice, δ44/42CaSerum inversely correlated with intestinal Ca uptake and histological bone mineralization markers, but not with Ca content and bone mineral density by µCT. In conclusion, δ44/42CaBone was bone site specific, but mechanical stress and gender independent. Low-Ca-diet induced marked changes in feces, serum and urine δ44/42Ca in growing rats. δ44/42CaSerum inversely correlated with markers of bone mineralization.

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: We developed a new version of the Alfred Wegener Institute Climate Model (AWI-CM3), which has higher skills in representing the observed climatology and better computational efficiency than its predecessors. Its ocean component FESOM2 (Finite-volumE Sea ice-Ocean Model) has the multi-resolution functionality typical of unstructured-mesh models while still featuring a scalability and efficiency similar to regular-grid models. The atmospheric component OpenIFS (CY43R3) enables the use of the latest developments in the numerical-weather-prediction community in climate sciences. In this paper we describe the coupling of the model components and evaluate the model performance on a variable-resolution (25-125 km) ocean mesh and a 61 km atmosphere grid, which serves as a reference and starting point for other ongoing research activities with AWI-CM3. This includes the exploration of high and variable resolution and the development of a full Earth system model as well as the creation of a new sea ice prediction system. At this early development stage and with the given coarse to medium resolutions, the model already features above-CMIP6-average skills (where CMIP6 denotes Coupled Model Intercomparison Project phase 6) in representing the climatology and competitive model throughput. Finally we identify remaining biases and suggest further improvements to be made to the model.

Description: Corals are globally important calcifiers that exhibit complex responses to anthropogenic warming and acidification. Although coral calcification is supported by high seawater pH, photosynthesis by the algal symbionts of zooxanthellate corals can be promoted by elevated pCO2. To investigate the mechanisms underlying corals’ complex responses to global change, three species of tropical zooxanthellate corals (Stylophora pistillata, Pocillopora damicornis, and Seriatopora hystrix) and one species of asymbiotic cold-water coral (Desmophyllum pertusum, syn. Lophelia pertusa) were cultured under a range of ocean acidification and warming scenarios. Under control temperatures, all tropical species exhibited increased calcification rates in response to increasing pCO2. However, the tropical species’ response to increasing pCO2 flattened when they lost symbionts (i.e., bleached) under the high-temperature treatments—suggesting that the loss of symbionts neutralized the benefit of increased pCO2 on calcification rate. Notably, the cold-water species that lacks symbionts exhibited a negative calcification response to increasing pCO2, although this negative response was partially ameliorated under elevated temperature. All four species elevated their calcifying fluid pH relative to seawater pH under all pCO2 treatments, and the magnitude of this offset (Δ[H+]) increased with increasing pCO2. Furthermore, calcifying fluid pH decreased along with symbiont abundance under thermal stress for the one species in which calcifying fluid pH was measured under both temperature treatments. This observation suggests a mechanistic link between photosymbiont loss (‘bleaching’) and impairment of zooxanthellate corals’ ability to elevate calcifying fluid pH in support of calcification under heat stress. This study supports the assertion that thermally induced loss of photosymbionts impairs tropical zooxanthellate corals’ ability to cope with CO2-induced ocean acidification.

Description: Microbial co-cultivation is a promising approach for the activation of biosynthetic gene clusters (BGCs) that remain transcriptionally silent under artificial culture conditions. As part of our project aiming at the discovery of marine-derived fungal agrochemicals, we previously used four phytopathogens as model competitors in the co-cultivation of 21 marine fungal strains. Based on comparative untargeted metabolomics analyses and anti-phytopathogenic activities of the co-cultures, we selected the co-culture of marine Cosmospora sp. with the phytopathogen Magnaporthe oryzae for in-depth chemical studies. UPLC-MS/MS-based molecular networking (MN) of the co-culture extract revealed an enhanced diversity of compounds in several molecular families, including isochromanones, specifically induced in the co-culture. Large scale co-cultivation of Cosmospora sp. and M. oryzae resulted in the isolation of five isochromanones from the whole co-culture extract, namely the known soudanones A, E, D (1-3) and their two new derivatives, soudanones H-I (4-5), the known isochromans, pseudoanguillosporins A and B (6, 7), naphtho-γ-pyrones, cephalochromin and ustilaginoidin G (8, 9), and ergosterol (10). Their structures were established by NMR, HR-ESIMS, FT-IR, electronic circular dichroism (ECD) spectroscopy, polarimetry ([α]D), and Mosher’s ester reaction. Bioactivity assays revealed antimicrobial activity of compounds 2 and 3 against the phytopathogens M. oryzae and Phytophthora infestans, while pseudoanguillosporin A (6) showed the broadest and strongest anti-phytopathogenic activity against Pseudomonas syringae, Xanthomonas campestris, M. oryzae and P. infestans. This is the first study assessing the anti-phytopathogenic activities of soudanones.

Description: Large amounts of methane (CH4) could be released as a result of the gradual or abrupt thawing of Arctic permafrost due to global warming. Once available, this potent greenhouse gas is emitted into the atmosphere or transported laterally into aquatic ecosystems via hydrologic connectivity at the surface or via groundwaters. While high northern latitudes contribute up to 5 % of total global CH4 emissions, the specific contribution of Arctic rivers and streams is largely unknown. We analyzed high-resolution continuous CH4 concentrations measured between 15 and 17 June 2019 (late freshet) in a ∼120 km transect of the Kolyma River in northeast Siberia. The average partial pressure of CH4 (pCH4) in tributaries (66.8–206.8 µatm) was 2–7 times higher than in the main river channel (28.3 µatm). In the main channel, CH4 was up to 1600 % supersaturated with respect to atmospheric equilibrium. Key sites along the riverbank and at tributary confluences accounted for 10 % of the navigated transect and had the highest pCH4 (41 ± 7 µatm) and CH4 emissions (0.03 ± 0.004 ) compared to other sites in the main channel, contributing between 14 % to 17 % of the total CH4 flux in the transect. These key sites were characterized by warm waters (T〉14.5 ∘C) and low specific conductivities (κ〈88 µS cm−1). The distribution of CH4 in the river could be linked statistically to T and κ of the water and to their proximity to the shore z, and these parameters served as predictors of CH4 concentrations in unsampled river areas. The abundance of CH4-consuming bacteria and CH4-producing archaea in the river was similar to those previously detected in nearby soils and was also strongly correlated to T and κ. These findings imply that the source of riverine CH4 is closely related with sites near land. The average total CH4 flux density in the river section was 0.02 ± 0.006 , equivalent to an annual CH4 flux of 1.24×107 g CH4 yr−1 emitted during a 146 d open water season. Our study highlights the importance of high-resolution continuous CH4 measurements in Arctic rivers for identifying spatial and temporal variations, as well as providing a glimpse of the magnitude of riverine CH4 emissions in the Arctic and their potential relevance to regional CH4 budgets.

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: Two linear proline-rich peptides (1–2), bearing an N-terminal pyroglutamate, were isolated from the marine bacterium Microbacterium sp. V1, associated with the marine sponge Petrosia ficiformis, collected in the volcanic CO2 vents in Ischia Island (South Italy). Peptide production was triggered at low temperature following the one strain many compounds (OSMAC) method. Both peptides were detected together with other peptides (3–8) via an integrated, untargeted MS/MS-based molecular networking and cheminformatic approach. The planar structure of the peptides was determined by extensive 1D and 2D NMR and HR-MS analysis, and the stereochemistry of the aminoacyl residues was inferred by Marfey’s analysis. Peptides 1–8 are likely to arise from Microbacterium V1 tailor-made proteolysis of tryptone. Peptides 1 and 2 were shown to display antioxidant properties in the ferric-reducing antioxidant power (FRAP) assay.

Description: Several strains of a Gram-negative, anaerobic photoautotrophic, motile, rod-shaped bacterium, designated as B14B, A-7R, and A-7Y were isolated from biofilms of low-mineralized soda lakes in central Mongolia and Russia (southeast Siberia). They had lamellar stacks as photosynthetic structures and bacteriochlorophyll a as the major photosynthetic pigment. The strains were found to grow at 25–35 °C, pH 7.5–10.2 (optimum, pH 9.0), and with 0–8% (w/v) NaCl (optimum, 0%). In the presence of sulfide and bicarbonate, acetate, butyrate, yeast extract, lactate, malate, pyruvate, succinate, and fumarate promoted growth. The DNA G + C content was 62.9–63.0 mol%. While the 16S rRNA gene sequences confirmed that the new strains belonged to the genus Ectothiorhodospira of the Ectothiorhodospiraceae, comparison of the genome nucleotide sequences of strains B14B, A-7R, and A-7Y revealed that the new isolates were remote from all described Ectothiorhodospira species both in dDDH (19.7–38.8%) and in ANI (75.0–89.4%). The new strains are also genetically differentiated by the presence of a nitric oxide reduction pathway that is lacking from all other Ectiothiorhodospiraceae. We propose to assign the isolates to the new species, Ectothiorhodospira lacustris sp. nov., with the type strain B14BT (=DSM 116064T = KCTC 25542T = UQM 41491T).

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: Marine heatwaves (MHWs) are extreme events characterized by abnormally high sea surface temperatures, and they have significant impacts on marine ecosystems and human society. The rapid and accurate forecasting of MHWs is crucial for preventing and responding to the impacts they can lead to. However, the research on relevant forecasting methods is limited, and a dedicated forecasting system specifically tailored for the South China Sea (SCS) region has yet to be reported. This study proposes a novel forecasting system utilizing U-Net and ConvLSTM models to predict MHWs in the SCS. Specifically, the U-Net model is used to forecast the intensity of MHWs, while the ConvLSTM model is employed to predict the probability of their occurrence. The indication of an MHW relies on both the intensity forecasted by the U-Net model exceeding threshold T and the occurrence probability predicted by the ConvLSTM model surpassing threshold P. Incorporating sensitivity analysis, optimal thresholds for T are determined as 0.9 °C, 0.8 °C, 1.0 °C, and 1.0 °C for 1-, 3-, 5-, and 7-day forecast lead times, respectively. Similarly, optimal thresholds for P are identified as 0.29, 0.30, 0.20, and 0.28. Employing these thresholds yields the highest forecast accuracy rates of 0.92, 0.89, 0.88, and 0.87 for the corresponding forecast lead times. This innovative approach gives better predictions of MHWs in the SCS, providing invaluable reference information for marine management authorities to make well-informed decisions and issue timely MHW warnings.

Description: Mesoscale eddies are ubiquitous in the ocean, yet our understanding of their evolutions, particularly eddy merging processes, remains enigmatic. In this study, the merging processes of two cyclonic–cyclonic and two anticyclonic–anticyclonic eddies are analyzed in the Subtropical Northwest Pacific using satellite remote sensing altimetry data. The results reveal that, as eddies approach each other, their contours become connected, leading to the formation of multi-core eddies. Simultaneously, the merging process prompts substantial exchanges of energy and vorticity, resulting in the dissipation of one eddy and the emergence of a more extensive merged eddy. Throughout the merging process, the eddy contours elongate upwards along the centerline (the line connecting eddy centers) and there are distinct changes in both the horizontal and vertical morphology of the eddies. Notably, after the merging, the eddies distinctly exhibit intensified signals of sea surface temperature and vertical temperature anomaly, an outcome of their transformative fusion. The findings of this study significantly enhance our understanding of mesoscale eddy dynamics, particularly in the intricate eddy merging process. However, it is important to note that, due to limitations in vertical observational data, this study does not provide a quantitative portrayal of the vertical mechanisms of eddy merging, which also underscores a pivotal avenue for future research in the field.

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: Biogeochemical cycling of carbon (C) and nitrogen (N) in the ocean depends on both the composition and activity of underlying biological communities and on abiotic factors. The Southern Ocean is encircled by a series of strong currents and fronts, providing a barrier to microbial dispersion into adjacent oligotrophic gyres. Our study region straddles the boundary between the nutrient-rich Southern Ocean and the adjacent oligotrophic gyre of the South Indian Ocean, providing an ideal region to study changes in microbial productivity. Here, we measured the impact of C- and N- uptake on microbial community diversity, contextualized by hydrographic factors and local physico-chemical conditions across the Southern Ocean and South Indian Ocean. We observed that contrasting physico-chemical characteristics led to unique microbial diversity patterns, with significant correlations between microbial alpha diversity and primary productivity (PP). However, we detected no link between specific PP (PP normalized by chlorophyll a concentration) and microbial alpha and beta diversity. Prokaryotic alpha and beta diversity were correlated with biological N2 fixation, itself a prokaryotic process, and we detected measurable N2 fixation to 60° S. While regional water masses have distinct microbial genetic fingerprints in both the eukaryotic and prokaryotic fractions, PP and N2 fixation vary more gradually and regionally. This suggests that microbial phylogenetic diversity is more strongly bounded by physical oceanographic features, while microbial activity responds more to chemical factors. We conclude that concomitant assessments of microbial diversity and activity is central in understanding the dynamics and complex responses of microorganisms to a changing ocean environment.

Description: We present an ensemble of last glacial inception (LGI) simulations for the Northern Hemisphere that captures a significant fraction of inferred ice volume changes within proxy uncertainties. This ensemble was performed with LCice 1.0, a coupled ice sheet and climate model, varying parameters of both climate and ice sheet components, as well as the coupling between them. Certain characteristics of the spatiotemporal pattern of ice growth and subsequent retreat in both North America (NA) and Eurasia (EA) are sensitive to parameter changes while others are not. We find that the initial inception of ice over NA and EA is best characterized by the nucleation of ice at high-latitude and high-elevation sites. Subsequent spreading and merger along with large-scale conversion of snowfields dominate in different sectors. The latter plays an important role in the merging of eastern and western ice regions in NA. The inception peak ice volume in the ensemble occurs approximately at 111 ka and therefore lags the summer 60∘ N insolation minimum by more than 3 kyr. Ice volumes consistently peak earlier over EA than NA. The inception peak in North America is characterized by a merged Laurentide and Cordilleran ice sheet, with the Davis Strait covered in ice in ∼80 % of simulations. Ice also bridges Greenland and Iceland in all runs by 114 ka and therefore blocks the Denmark Strait. This latter feature would thereby divert the East Greenland Current and Denmark Strait overflow with a potentially significant impact on ocean circulation. The Eurasian ice sheet at its inception peak varies across ensemble runs between a continuous ice sheet and multiple smaller ice caps. In both continents, the colder high latitudes (i.e. Ellesmere and Svalbard) tend to grow ice through the entire simulation (to 102 ka), while lower latitudes lose ice after ∼110 ka. We find temperature decreases over the initial phases of the inception lead to the expansion of NA ice sheet area and that subsequent precipitation increases contribute to its thickening. EA ice sheet area also expands with decreasing temperatures, but sea ice limits any increases in precipitation, leading to an earlier retreat away from the EA maximum ice sheet volume. We also examine the extent to which the capture of both LGI ice growth and retreat constrains the coupled ice–climate model sensitivity to changing atmospheric pCO2. The 55-member sub-ensemble that meets our criteria for “acceptable” ice growth and retreat has an equilibrium climate sensitivity lower bound that is 0.3 ∘C higher than that of the full ensemble. This suggests some potential value of fully coupled ice–climate modelling of the last glacial inception to constrain future climate change.

Description: Open-ocean oxygen minimum zones (OMZs) occur in regions with high biological productivity and weak ventilation. They restrict marine habitats and alter biogeochemical cycles. Global models generally show a large model–data misfit with regard to oxygen. Reliable statements about the future development of OMZs and the quantification of their interaction with climate change are currently not possible. One of the most intense OMZs worldwide is located in the Arabian Sea (AS). We give an overview of the main model deficiencies with a detailed comparison of the historical state of 10 climate models from the 5th Coupled Model Intercomparison Project (CMIP5) that present our present-day understanding of physical and biogeochemical processes. Most of the models show a general underestimation of the OMZ volume in the AS compared to observations that is caused by an overly shallow layer of oxygen-poor water in the models. The deviation of oxygen values in the deep AS is the result of oxygen levels that are too high simulated in the Southern Ocean formation regions of Indian Ocean Deep Water in the models compared to observations and uncertainties in the deepwater mass transport from the Southern Ocean northward into the AS. Differences in simulated water mass properties and ventilation rates of Red Sea Water and Persian Gulf Water cause different mixing in the AS and thus influence the intensity of the OMZ. These differences in ventilation rates also point towards variations in the parameterizations of the overflow from the marginal seas among the models. The results of this study are intended to foster future model improvements regarding the OMZ in the AS.