ALBERT

Description: The knowledge of tree species dependent turnover of soil organic matter (SOM) is limited, yet required to understand the carbon sequestration function of forest soil. We combined investigations of 13C and 15N and its relationship to elemental stoichiometry along soil depth gradients in 35-year old monocultural stands of Douglas fir (Pseudotsuga menziesii), black pine (Pinus nigra), European beech (Fagus sylvatica) and red oak (Quercus rubra) growing on a uniform post-mining soil. We investigated the natural abundance of 13C and 15N and the carbon:nitrogen (C:N) and oxygen:carbon (O:C) stoichiometry of litterfall and fine roots as well as SOM in the forest floor and mineral soil. Tree species had a significant effect on SOM δ13C and δ15N reflecting significantly different signatures of litterfall and root inputs. Throughout the soil profile, δ13C and δ15N were significantly related to the C:N and O:C ratio which indicates that isotope enrichment with soil depth is linked to the turnover of organic matter (OM). Significantly higher turnover of OM in soils under deciduous tree species depended to 46% on the quality of litterfall and root inputs (N content, C:N, O:C ratio), and the initial isotopic signatures of litterfall. Hence, SOM composition and turnover also depends on additional—presumably microbial driven—factors. The enrichment of 15N with soil depth was generally linked to 13C. In soils under pine, however, with limited N and C availability, the enrichment of 15N was decoupled from 13C. This suggests that transformation pathways depend on litter quality of tree species.

Description: Universität Trier (3163)

Description: Despite the importance of phosphorus (P) as a macronutrient, the factors controlling the pool sizes of organic and inorganic P (OP and IP) in soils are not yet well understood. Therefore, the aim of this study was to gain insights into the pools sizes of OP, IP and organic carbon (OC) in soils and soil particle size fractions. For this purpose, I analyzed the distribution of OP, IP, and OC among particle size fractions depending on geographical location, climate, soil depth, and land use, based on published data. The clay size fraction contained on average 8.8 times more OP than the sand size fraction and 3.9 and 3.2 times more IP and OC, respectively. The OP concentrations of the silt and clay size fraction were both negatively correlated with mean annual temperature (R2 = 0.30 and 0.31, respectively, p 〈 0.001). The OC:OP ratios of the silt and clay size fraction were negatively correlated with latitude (R2 = 0.49 and 0.34, respectively, p 〈 0.001). Yet, the OC:OP ratio of the clay size fraction changed less markedly with latitude than the OC:OP ratio of the silt and the sand size fraction. The OC concentrations of all three particle size fractions were significantly (p 〈 0.05) lower in soils converted to cropland than in adjacent soils under natural vegetation. In contrast, the OP concentration was only significantly (p 〈 0.05) decreased in the sand size fraction but not in the other two particle size fractions due to land-use change. Thus, the findings suggest that OP is more persistent in soil than OC, which is most likely due to strong sorptive stabilization of OP compounds to mineral surfaces.

Description: German Research Foundation

Description: Dissolved organic carbon (DOC) from Oa horizons has been proposed to be an important contributor for subsoil organic carbon stocks. We investigated the fate of DOC by directly injecting a DOC solution from 13C labelled litter into three soil depths at beech forest sites. Fate of injected DOC was quantified with deep drilling soil cores down to 2 m depth, 3 and 17 months after the injection. 27 ± 26% of the injected DOC was retained after 3 months and 17 ± 22% after 17 months. Retained DOC was to 70% found in the first 10 cm below the injection depth and on average higher in the topsoil than in the subsoil. After 17 months DOC in the topsoil was largely lost (− 19%) while DOC in the subsoil did not change much (− 4.4%). Data indicated a high stabilisation of injected DOC in the subsoils with no differences between the sites. Potential mineralisation as revealed by incubation experiments however, was not different between DOC injected in topsoil or subsoils underlining the importance of environmental factors in the subsoil for DOC stabilisation compared to topsoil. We conclude that stability of DOC in subsoil is primary driven by its spatial inaccessibility for microorganisms after matrix flow while site specific properties did not significantly affect stabilisation. Instead, a more fine-textured site promotes the vertical transport of DOC due to a higher abundance of preferential flow paths.

Description: Deutsche Forschungsgemeinschaft http://dx.doi.org/10.13039/501100001659

Description: Calcium (Ca) plays a crucial role for plant nutrition, soil aggregation, and soil organic matter (SOM) stabilization. Turnover and ecological functions of Ca in soils depend on soil Ca speciation. For the first time, we used synchrotron-based X-ray absorption near-edge structure (XANES) spectroscopy at the Ca K-edge (4038 eV) to investigate Ca speciation in soils. We present Ca K-edge XANES spectra of standard compounds with relevance in soils (e.g. calcite, dolomite, hydroxyapatite, anorthite, clay mineral-adsorbed Ca; Ca oxalate, formate, acetate, citrate, pectate, phytate). Calcium XANES spectra with good signal-to-noise ratios were acquired in fluorescence mode for Ca concentrations between 1 and 10 mg g−1. Most standard spectra differed markedly among each other, allowing the identification of different Ca species in soils and other environmental samples as well as Ca speciation by linear combination fitting. Calcium XANES spectra obtained for samples from different horizons of twelve temperate forest soils revealed a change from dominating lithogenic Ca to clay mineral-bound and/or organically bound Ca with advancing pedogenesis. O layer Ca was almost exclusively organically bound. With increasing SOM decomposition, shares of oxalate-bound Ca decreased. Oxalate-bound Ca was absent in calcareous, but not in silicate subsoil horizons, which can be explained by microbial decomposition in the former vs. stabilization by association to pedogenic minerals in the latter soils. Synchrotron-based Ca XANES spectroscopy is a promising novel tool to investigate the fate of Ca during pedogenesis and—when performed with high spatial resolution (µ-XANES), to study aggregation and SOM stabilization mechanisms produced by Ca.

Description: BLE Waldklimafonds

Description: Projekt DEAL

Description: Ge/Si ratios of plant phytoliths have been widely used to trace biogeochemical cycling of Si. However, until recently, information on how much of the Ge and Si transferred from soil to plants is actually stored in phytoliths was lacking. The aim of the present study is to (i) compare the uptake of Si and Ge in three grass species, (ii) localize Ge and Si stored in above-ground plant parts and (iii) evaluate the amounts of Ge and Si sequestrated in phytoliths and plant tissues. Mays (Zea mays), oat (Avena sativa) and reed canary grass (Phalaris arundinacea) were cultivated in the greenhouse on soil and sand to control element supply. Leaf phytoliths were extracted by dry ashing. Total elemental composition of leaves, phytoliths, stems and roots were measured by ICP-MS. For the localization of phytoliths and the determination of Ge and Si within leaf tissues and phytoliths scanning electron microscopy (SEM), energy dispersive x-ray spectroscopy (EDX) and laser ablation inductively coupled mass spectrometry (LA-ICP-MS) was used. The amounts of Si and Ge taken up by the species corresponded with biomass formation and decreased in the order Z. mays 〉 P. arundinacea, A. sativa. Results from LA-ICP-MS revealed that Si was mostly localized in phytoliths, while Ge was disorderly distributed within the leaf tissue. In fact, from the total amounts of Ge accumulated in leaves only 10% was present in phytoliths highlighting the role of organic matter on biogeochemical cycling of Ge and the necessity for using bulk Ge/Si instead of Ge/Si in phytoliths to trace biogeochemical cycling of Si.

Description: Tree roots penetrate the soil to several meters depth, but the role of subsoils for the supply of nutrient elements such as phosphorus (P) to the trees is poorly understood. Here, we tested the hypothesis that increased P deficiency in the topsoil results in an increased microbial recycling of P from the forest subsoil. We sampled soils from four German temperate forest sites representing a gradient in total P stocks. We analyzed the oxygen isotopic composition of HCl-extractable phosphate (δ18OP) and identified differences in P speciation with increasing soil depth using X-ray absorption near-edge structure (XANES) spectroscopy. We further determined microbial oxygen demand with and without nutrient supply at different soil depths to analyse nutrient limitation of microbial growth and used nanoscale secondary ion mass spectrometry (NanoSIMS) to visualize spatial P gradients in the rhizosphere. We found that δ18OP values in the topsoil of all sites were close to the isotopic signal imparted by biological cycling when oxygen isotopes in phosphate are exchanged by enzymatic activity. However, with increasing soil depth and increasing HCl-P concentrations, δ18Ο values continuously decreased towards values expected for primary minerals in parent material at depths below 60 cm at sites with high subsoil P stocks and below more than 2 m at sites with low subsoil P stocks, respectively. For these depths, XANES spectra also indicated the presence of apatite. NanoSIMS images showed an enrichment of P in the rhizosphere in the topsoil of a site with high P stocks, while this P enrichment was absent at a site with low P stocks and in both subsoils. Addition of C, N and P alone or in combination revealed that microbial activity in subsoils of sites with low P stocks was mostly P limited, whereas sites with high P stocks indicated N limitation or N and P co-limitation. We conclude that subsoil P resources are recycled by trees and soil microorganisms. With continued weathering of the bedrock and mobilisation of P from the weathered rocks, P cycling will proceed to greater depths, especially at sites characterised by P limitation.

Description: Rheinische Friedrich-Wilhelms-Universität Bonn (1040)

Description: Deutsche Forschungsgemeinschaft http://dx.doi.org/10.13039/501100001659

Description: Rheinische Friedrich-Wilhelms-Universität Bonn (1040)

Description: All nuclear explosions are banned by the Comprehensive Nuclear-Test-Ban Treaty. In the context of the treaty a verification regime was put into place to detect, locate, and characterize nuclear explosions at any time, by anyone and everywhere on the Earth. The International Monitoring System, which plays a key role in the verification regime, was set up by the Preparatory Commission of the Comprehensive Nuclear-Test-Ban Treaty Organization. Out of the several different monitoring techniques applied in the International Monitoring System the seismic waveform approach is the most effective and reliable technology for monitoring nuclear explosions underground. This study introduces a deterministic method of threshold monitoring that allows to asses a lower body wave magnitude limit of a potential seismic event in a certain geographical region, that can be detected by those seismic stations being part of the International Monitoring System network. The method is based on measurements of ambient seismic noise levels at the individual seismic stations along with global distance corrections terms for the body wave magnitude. The results suggest that an average global detection capability of approximately body wave magnitude 4.0 can be achieved using only stations from the primary seismic network of the International Monitoring System. The incorporation of seismic stations from the auxiliary seismic network leads to a slight improvement of the detection capability, while the use and analysis of wave arrivals from distances greater than 120∘ results in a significant improvement of the detection capability. Temporal variations in terms of hourly and monthly changes of the global detection capability can not be observed. Overall, comparisons between detection capability and manually retrieved body wave magnitudes from the Reviewed Event Bulletin suggest, that our method yields a more conservative estimation of the detection capability and that in reality detection thresholds might be even lower than estimated.

Description: Bundesanstalt für Geowissenschaften und Rohstoffe (BGR) (4230)

Description: Electrical resistivity tomography (ERT) and ground magnetic surveys were applied to characterize an old uncontrolled landfill in a former exploited sand and gravel quarry in an area to the north-west of the city of Cologne, Germany. The total magnetic field and its vertical gradient were recorded using a proton precession magnetometer to cover an area of about 43,250 m2. The magnetic data were transferred to the frequency domain and then reduced to the north magnetic pole. The amplitude of the analytical signal was calculated to define the magnetic materials within and outside the landfill. Eight ERT profiles were constructed based on the results of the magnetic survey using different electrode arrays (Wenner, dipole–dipole, and Schlumberger). In order to increase both data coverage and sensitivity and to decrease uncertainty, a non-conventional mixed array was used. The subsurface resistivity distributions were imaged using the robust (L1-norm) inversion method. The resultant inverted subsurface true resistivity data were presented in the form of 2D cross sections and 3D fence diagram. These non-invasive geophysical tools helped us to portray the covering soil, the spatial limits of the landfill, and the depth of the waste body. We also successfully detected low resistivity zones at deeper depths than expected, which probably be associated with migration pathways of the leachate plumes. The findings of the present study provide valuable information for decision makers with regards to environmental monitoring and assessment.

Description: Volkswagen Foundation http://dx.doi.org/10.13039/501100001663

Description: Universität zu Köln (1017)

Description: An experimental multi-parameter structural monitoring system has been installed on the Kurpsai dam, western Kyrgyz Republic. This system consists of equipment for seismic and strain measurements for making longer- (days, weeks, months) and shorter- (minutes, hours) term observations, dealing with, for example seasonal (longer) effects or the response of the dam to ground motion from noise or seismic events. Fibre-optic strain sensors allow the seasonal and daily opening and closing of the spaces between the dam’s segments to be tracked. For the seismic data, both amplitude (in terms of using differences in amplitudes in the Fourier spectra for mapping the modes of vibration of the dam) and their time–frequency distribution for a set of small to moderate seismic events are investigated and the corresponding phase variabilities (in terms of lagged coherency) are evaluated. Even for moderate levels of seismic-induced ground motion, some influence on the structural response can be detected, which then sees the dam quickly return to its original state. A seasonal component was identified in the strain measurements, while levels of noise arising from the operation of the dam's generators and associated water flow have been provisionally identified.

Description: Bundesministerium für Forschung und Technologie http://dx.doi.org/10.13039/501100004937

Description: Helmholtz-Zentrum Potsdam Deutsches GeoForschungsZentrum - GFZ (4217)

Description: Plant litter decomposition is a key ecological process that is mostly studied at the forest floor. However, decomposition generally starts in the canopy. In this study, we evaluated the effect of litter composition and climate on the initial phase of decomposition in the soil and two contrasting types of canopy microsites along an elevational gradient (0–2200 m a.s.l.). To this end, we incubated standard material composed by green (fast decomposing) and rooibos (slow decomposing) tea bags for three months. Tea bags were placed in soil (buried at 5 cm) and in the canopy at ca. 5 m above the ground in “micro-wetlands” (tank bromeliads) and dry crown microsites (branches). Along the elevational gradient, green tea decomposed faster than rooibos tea in all microsites and forests. Mass loss for both tea types was lowest on branches at all sites, except for green tea in a wet forest where decomposition did not significantly differ among microsites. In wet forests, decomposition did not differ between bromeliads and soil, while in a dry forest, decomposition was faster in bromeliads. We found that the effects of climatic variables [monthly average temperature (TEMP) and total precipitation (PREC) for the incubation months] on decomposition differed between microsites. Along the elevational gradient, the mass loss in soil was positively correlated with TEMP but not with PREC, whereas on branches, mass loss was negatively correlated with TEMP and positively correlated with PREC. Unlike on branches, mass loss in bromeliads slightly decreased with PREC and increased with TEMP. Our study shows that microsite conditions interact with climate (TEMP and PREC) leading to differences in the general decomposition patterns in the forest canopy.

Description: CONACYT

Description: Deutscher Akademischer Austauschdienst http://dx.doi.org/10.13039/501100001655

Description: Heinz Neumüller Stiftung

Description: Carl von Ossietzky Universität Oldenburg (3092)

Description: Selective logging is among the main causes of tropical forest degradation, but little is known about its effects on greenhouse gas (GHG) fluxes from highly weathered Ferralsol soils in Africa. We measured soil CO2, N2O, and CH4 fluxes, and their soil controlling factors at two forests that had undergone conventional selective logging and reduced-impact logging in Cameroon. Each logging system had four replicate plots, each included the disturbed strata (road, logging deck, skidding trail, and felling gap) and an undisturbed reference area. Measurements were conducted monthly from September 2016 to October 2017. Annual GHG fluxes ranged from 4.9 to 18.6 Mg CO2–C, from 1.5 to 79 kg N2O–N, and from − 4.3 to 71.1 kg CH4–C ha−1 year−1. Compared to undisturbed areas, soil CO2 emissions were reduced and soil CH4 emissions increased in skidding trails, logging decks and roads (P 〈 0.01) whereas soil N2O emissions increased in skidding trails (P = 0.03–0.05). The combined disturbed strata had 28% decrease in soil CO2 emissions, 83% increase in soil N2O emissions, and seven times higher soil CH4 emissions compared to undisturbed area (P ≤ 0.01). However, the disturbed strata represented only 4–5% of the area impacted in both logging systems, which reduced considerably the changes in soil GHG fluxes at the landscape level. Across all strata, soil GHG fluxes were regulated by soil bulk density and water-filled pore space, indicating the influence of soil aeration and gas diffusion, and by soil organic carbon and nitrogen, suggesting the control of substrate availability on microbial processes of these GHG.

Description: Deutsche Forschungsgemeinschaft http://dx.doi.org/10.13039/501100001659

Description: Deutscher Akademischer Austauschdienst http://dx.doi.org/10.13039/501100001655

Description: Georg-August-Universität Göttingen (1018)

Description: https://doi.org/10.25625/TWSFRU

Description: Coastal waters have strong gradients in dissolved organic matter (DOM) quantity and characteristics, originating from terrestrial inputs and autochthonous production. Enclosed seas with high freshwater input therefore experience high DOM concentrations and gradients from freshwater sources to more saline waters. The brackish Baltic Sea experiences such salinity gradients from east to west and from river mouths to the open sea. Furthermore, the catchment areas of the Baltic Sea are very diverse and vary from sparsely populated northern areas to densely populated southern zones. Coastal systems vary from enclosed or open bays, estuaries, fjords, archipelagos and lagoons where the residence time of DOM at these sites varies and may control the extent to which organic matter is biologically, chemically or physically modified or simply diluted with transport off-shore. Data of DOM with simultaneous measurements of dissolved organic (DO) nitrogen (N), carbon (C) and phosphorus (P) across a range of contrasting coastal systems are scarce. Here we present data from the Roskilde Fjord, Vistula and Öre estuaries and Curonian Lagoon; four coastal systems with large differences in salinity, nutrient concentrations, freshwater inflow and catchment characteristics. The C:N:P ratios of DOM of our data, despite high variability, show site specific significant differences resulting largely from differences residence time. Microbial processes seemed to have minor effects, and only in spring did uptake of DON in the Vistula and Öre estuaries take place and not at the other sites or seasons. Resuspension from sediments impacts bottom waters and the entire shallow water column in the Curonian Lagoon. Finally, our data combined with published data show that land use in the catchments seems to impact the DOC:DON and DOC:DOP ratios of the tributaries most.

Description: Academy of Sciences of Finland

Description: Deutsche Forschungsgemeinschaft http://dx.doi.org/10.13039/501100001659

Description: BONUS COCOA and Bundesministerium für Bildung und Forschung

Description: Academy of Finland

Description: Danish Research Council for Independent Research

Description: BONUS COCOA Project

Description: Leibniz-Institut für Ostseeforschung Warnemünde (IOW) (3484)

Description: Even though micropolar theories are widely applied for engineering applications such as the design of metamaterials, applications in the study of the Earth’s interior still remain limited and in particular in seismology. This is due to the lack of understanding of the required elastic material parameters present in the theory as well as the eigenfrequency $\omega _r$ which is not observed in seismic data. By showing that the general dynamic equations of the Timoshenko’s beam is a particular case of the micropolar theory we are able to connect micropolar elastic parameters to physically measurable quantities. We then present an alternative micropolar model that, based on the same physical basis as the original model, circumvents the problem of the original eigenfrequency $\omega _r$ laking in seismological data. We finally validate our model with a seismic experiment and show it is relevant to explain observed seismic dispersion curves.

Description: Westfälische Wilhelms-Universität Münster (1056)

Description: Seismic events produced by block rotations about vertical axis occur in many geodynamic contexts. In this study, we show that these rotations can be accounted for using the proper theory, namely micropolar theory, and a new asymmetric moment tensor can be derived. We then apply this new theory to the Kaikōura earthquake (2016/11/14), Mw 7.8, one of the most complex earthquakes ever recorded with modern instrumental techniques. Using advanced numerical techniques, we compute synthetic seismograms including a full asymmetric moment tensor and we show that it induces measurable differences in the waveforms proving that seismic data can record the effects of the block rotations observed in the field. Therefore, the theory developed in this work provides a full framework for future dynamic source inversions of asymmetric moment tensors.

Description: Westfälische Wilhelms-Universität Münster (1056)

Description: Two efficient and scalable numerical solution methods will be compared using exact Jacobians to solve the fully coupled Newton systems arising during fully implicit discretization of the equations for two-phase flow in porous media. These methods use algebraic multigrid (AMG) to solve the linear systems in every Newton step. The algebraic multigrid methods rely on (i) a Schur Complement Reduction (SCR-AMG) and (ii) a Constrained Pressure Residual method (CPR-AMG) to decouple elliptic and hyperbolic contributions. Both methods employ automatic differentiation (AD) to calculate exact Jacobians and are compared with finite difference (FD) approximations of the Jacobian. The superiority of AD is shown by several numerical test cases from the field of CO2 geo-sequestration comprising two- and three-dimensional examples. A weak scaling test on JUQUEEN, a BlueGene/Q supercomputer, demonstrates the efficiency and scalability of both methods. To achieve maximum comparability and reproducibility, the Portable Extensible Toolkit for Scientific Computation (PETSc) is used as framework for the comparison of all solvers.

Description: Horizon 2020 https://doi.org/10.13039/501100007601

Description: Horizon 2020 https://doi.org/10.13039/501100007601

Description: A multiphase-field approach for elasto-plastic and anisotropic brittle crack propagation in geological systems consisting of different regions of brittle and ductile materials is presented and employed to computationally study crack propagation. Plastic deformation in elasto-plastic materials such as frictional, granular or porous materials is modelled with the pressure-sensitive Drucker-Prager plasticity model. This plasticity model is combined with a multiphase-field model fulfilling the mechanical jump conditions in diffuse solid-solid interfaces. The validity of the plasticity model with phase-inherent stress and strain fields is shown, in comparison with sharp interface finite element solutions. The proposed model is capable of simulating crack formation in heterogeneous multiphase systems comprising both purely elastic and inelastic phases. We investigate the influence of different material parameters on the crack propagation with tensile tests in single- and two-phase materials. To show the applicability of the model, crack propagation in a multiphase domain with brittle and elasto-plastic components is performed.

Description: The composition of meta-gabbro and meta-basalt occurring abundant and widespread in all nappes of the nappe stack exposed in the Zermatt region of the Western Alp shows distinct patterns related to the geodynamic origin of metamorphic basic rocks. Eclogitic meta-basalts of the ophiolitic Zermatt-Saas Unit (ZSU) show enriched MORB signatures. The meta-basalts (eclogites) of the continental fragment of the Theodul Glacier Unit (TGU) derive from pre-Alpine metamorphic continental intraplate basalts. Meta-basalts (eclogites) from the continental basement of the Siviez-Mischabel nappe (SMN) derive from MORB thus a genetic relation to the TGU eclogites can be excluded. All basic igneous rocks experienced post-magmatic alteration by fluid-rock interaction ranging from processes at the seafloor, in the shallow crust, during subduction zone hydration, in the exhumation channel and late Alpine regional metamorphisms. The consequences of these alteration processes can be identified at various levels in the rock composition data. It was found that the REE data are little affected by fluid-rock alteration. Some trace elements, notably Cs, Rb, and Ba are typically massively altered relative to igneous compositions in all three groups of meta-basalts. Generally, meta-basalts from the TGU and the SMN preserved the features of the original composition whilst the ZSU meta-volcanic rocks experienced massive alteration. For the ZSU meta-volcanic rocks it is evident that Zr was gained and Y lost during high-pressure fluid-rock interaction indicating a mobile behavior of the two elements during HP-metamorphism in contrast to their behavior in hydrothermal near-surface fluid-rock interaction.

Description: Albert-Ludwigs-Universität Freiburg im Breisgau (1016)

Description: Significant sedimentation of manganese (Mn) in form of manganese oxides (MnOx) and the subsequent formation of authigenic calcium-rich rhodochrosite (Mn(Ca)CO3) were observed in the seasonally stratified hard water Lake Stechlin in north-eastern Germany. This manganese enrichment was assumed to be associated with recent eutrophication of the formerly oligotrophic lake. The mechanisms and processes involved were examined by analysing: (i) short sediment cores obtained from seven locations along a depth transect ranging from 69.5 m (the deepest point) to 38 m; (ii) sediment traps located at 20 m and 60 m water depths; (iii) water column profiles; and (iv) porewater profiles at 69.5 m and 58 m depths. Sedimentary Mn enrichment was observed at water depths below 56 m and increased to more than 25 wt% at the deepest site. Between 2010 and 2017, Mn accumulation at the deepest site was 815 g Mn m−2. Transfer of Mn from the shallower towards the deepest parts of the lake was initiated by reductive dissolution of MnOx and diffusion of dissolved Mn from the sediment to the overlying water column. Manganese was then dissipated via turbulent mixing and subsequently oxidised to MnOx before being transported towards the deepest zone. Transformation of the redeposited MnOx to Mn(Ca)CO3 favoured the final burial of Mn. We show that eutrophication and the areal spreading of anoxic conditions may intensify diagenetic processes and cause the spatial redistribution of Mn as well as its effective burial. Contrary to many previous findings, we show that increases of Mn and Mn/Fe can also be used as indicators for increasing anoxic conditions in previously oligotrophic lakes.

Description: Leibniz-Gemeinschaft http://dx.doi.org/10.13039/501100001664

Description: Leibniz-Institut für Gewässerökologie und Binnenfischerei (IGB) im Forschungsverbund Berlin e.V. (3473)

Description: For the purpose of monitoring for compliance with the Comprehensive Nuclear-Test-Ban Treaty (CTBT), the International Monitoring System (IMS) is being established that includes 40 sensor systems for atmospheric xenon radioactivity. Its purpose is to monitor the atmosphere for signatures that may indicate a nuclear explosion. Normal operational releases of radioxenon from nuclear facilities can regularly be observed by these very high-sensitive noble gas systems. Existing best estimates of releases for a generic year are unlikely to apply for any specific year at the level of individual facilities because their releases are highly variable and can change by several orders of magnitude from year to year. In this paper, best knowledge of the radioxenon emission inventory from nuclear power plants (NPPs) is collected for the calendar year 2014. The distribution function for each CTBT relevant radioxenon isotope is derived from all releases from NPPs as reported for 2014. The data of this paper can be used for developing and validating methods based on atmospheric transport modelling that are designed to enhance understanding of the impact of known sources on the IMS background observations.

Description: Universität Hamburg (1037)

Description: Even though micropolar theories are widely applied for engineering applications such as the design of metamaterials, applications in the study of the Earth’s interior still remain limited and in particular in seismology. This is due to the lack of understanding of the required elastic material parameters present in the theory as well as the eigenfrequency ωr which is not observed in seismic data. By showing that the general dynamic equations of the Timoshenko’s beam is a particular case of the micropolar theory we are able to connect micropolar elastic parameters to physically measurable quantities. We then present an alternative micropolar model that, based on the same physical basis as the original model, circumvents the problem of the original eigenfrequency ωr laking in seismological data. We finally validate our model with a seismic experiment and show it is relevant to explain observed seismic dispersion curves.

Description: Westfälische Wilhelms-Universität Münster (1056)

Description: The explosion at the Ingolstadt oil refinery was widely recorded at seismic and infrasound stations deployed throughout Central Europe, to distances of several hundred to a thousand kilometres. This study focuses on the wealth of data recorded at infrasound stations in Central and Eastern Europe, while from the many detecting seismic stations within 400 km range, only seismic and seismo-acoustic arrivals at the close-in Gräfenberg array are considered here. Most of the infrasound stations are acoustic arrays enabling us to apply array processing techniques to determine relevant wave field parameters, such as backazimuth and slowness (resp. trace velocity). These parameters not only confirm the source direction, but also put constraints on the observed arrivals’ propagation modes. Wave field parameters suggest that we observe tropospheric arrivals to about 150 km and stratospheric and/or thermospheric returns for longer distances. 1D, 2D and 3D ray tracing predict tropospheric arrivals to westerly directions up to distances of 100 km, beyond which only thermospheric returns are obtained azimuth-independent beyond 250–300 km. Stratospheric returns do not follow from any of the increasingly complex ray tracing models. Parabolic equation propagation modeling however suggests that in a number of cases stratospheric ducting may be possible. However, neither the tropospheric seismo-acoustic arrivals at the Gräfenberg array nor the various arrivals at IMS station IS26 could be modeled. Therefore, the Ingolstadt explosion along with the observed infrasonic phases provide an excellent test bed to investigate our ability in realistically forecasting atmospheric wave propagation with existing algorithms and available atmospheric models.

Description: Managed grasslands have the potential to store carbon (C) and partially mitigate climate change. However, it remains difficult to predict potential C storage under a given soil or management practice. To study C storage dynamics due to long-term (1952–2009) phosphorus (P) fertilizer and irrigation treatments in New Zealand grasslands, we measured radiocarbon (〈sup〉14〈/sup〉C) in archived soil along with observed changes in C stocks to constrain a compartmental soil model. Productivity increases from P application and irrigation in these trials resulted in very similar C accumulation rates between 1959 and 2009. The ∆〈sup〉14〈/sup〉C changes over the same time period were similar in plots that were both irrigated and fertilized, and only differed in a non-irrigated fertilized plot. Model results indicated that decomposition rates of fast cycling C (0.1 to 0.2 year〈sup〉−1〈/sup〉) increased to nearly offset increases in inputs. With increasing P fertilization, decomposition rates also increased in the slow pool (0.005 to 0.008 year〈sup〉−1〈/sup〉). Our findings show sustained, significant (i.e. greater than 4 per mille) increases in C stocks regardless of treatment or inputs. As the majority of fresh inputs remain in the soil for less than 10 years, these long term increases reflect dynamics of the slow pool. Additionally, frequent irrigation was associated with reduced stocks and increased decomposition of fresh plant material. Rates of C gain and decay highlight trade-offs between productivity, nutrient availability, and soil C sequestration as a climate change mitigation strategy.

Description: European Research Council http://dx.doi.org/10.13039/501100000781

Description: GNS Science

Description: Max Planck Institute for Biogeochemistry (2)

Description: https://github.com/ShaneStoner/Winchmore14C

Description: Freeze–thaw (FT) events exert a great physiological stress on the soil microbial community and thus significantly impact soil biogeochemical processes. Studies often show ambiguous and contradicting results, because a multitude of environmental factors affect biogeochemical responses to FT. Thus, a better understanding of the factors driving and regulating microbial responses to FT events is required. Soil chronosequences allow more focused comparisons among soils with initially similar start conditions. We therefore exposed four soils with contrasting organic carbon contents and opposing soil age (i.e., years after restoration) from a postmining agricultural chronosequence to three consecutive FT events and evaluated soil biochgeoemical responses after thawing. The major microbial biomass carbon losses occurred after the first FT event, while microbial biomass N decreased more steadily with subsequent FT cycles. This led to an immediate and lasting decoupling of microbial biomass carbon:nitrogen stoichiometry. After the first FT event, basal respiration and the metabolic quotient (i.e., respiration per microbial biomass unit) were above pre-freezing values and thereafter decreased with subsequent FT cycles, demonstrating initially high dissimilatory carbon losses and less and less microbial metabolic activity with each iterative FT cycle. As a consequence, dissolved organic carbon and total dissolved nitrogen increased in soil solution after the first FT event, while a substantial part of the liberated nitrogen was likely lost through gaseous emissions. Overall, high-carbon soils were more vulnerable to microbial biomass losses than low-carbon soils. Surprisingly, soil age explained more variation in soil chemical and microbial responses than soil organic carbon content. Further studies are needed to dissect the factors associated with soil age and its influence on soil biochemical responses to FT events.

Description: Bundesministerium für Bildung und Forschung http://dx.doi.org/10.13039/501100002347

Description: Universität zu Köln (1017)

Description: https://datenzentrum.bonares.de/data-portal.php

Description: Natural earthquakes often have very few observable foreshocks which significantly complicates tracking potential preparatory processes. To better characterize expected preparatory processes before failures, we study stick-slip events in a series of triaxial compression tests on faulted Westerly granite samples. We focus on the influence of fault roughness on the duration and magnitude of recordable precursors before large stick–slip failure. Rupture preparation in the experiments is detectable over long time scales and involves acoustic emission (AE) and aseismic deformation events. Preparatory fault slip is found to be accelerating during the entire pre-failure loading period, and is accompanied by increasing AE rates punctuated by distinct activity spikes associated with large slip events. Damage evolution across the fault zones and surrounding wall rocks is manifested by precursory decrease of seismic b-values and spatial correlation dimensions. Peaks in spatial event correlation suggest that large slip initiation occurs by failure of multiple asperities. Shear strain estimated from AE data represents only a small fraction (〈 1%) of total shear strain accumulated during the preparation phase, implying that most precursory deformation is aseismic. The relative contribution of aseismic deformation is amplified by larger fault roughness. Similarly, seismic coupling is larger for smooth saw-cut faults compared to rough faults. The laboratory observations point towards a long-lasting and continuous preparation process leading to failure and large seismic events. The strain partitioning between aseismic and observable seismic signatures depends on fault structure and instrument resolution.

Description: Helmholtz-Zentrum Potsdam - Deutsches GeoForschungsZentrum GFZ http://dx.doi.org/10.13039/501100010956

Description: Helmholtz-Zentrum Potsdam Deutsches GeoForschungsZentrum - GFZ (4217)

Description: In the past, several destructive earthquakes have occurred in the North African Atlas Mountain ranges located along the Africa–Eurasia plate boundary. Although the region is rich with impressive archaeological sites, including those in modern Tunisia, few comprehensive archaeoseismological studies have been conducted. Historic sources account at least three damaging earthquakes in the Kairouan area in central Tunisia between AD 859 and 1041. Little is known about which faults triggered these earthquakes or the size of these events. The water supply of the city of Kairouan depended on a 32-km-long aqueduct with a large bridge (now partially collapsed) at the confluence of the de Mouta and Cherichira rivers. The original bridge of Roman construction was retrofitted twice during the Aghlabid period (AD 800–903) and probably in AD 995 during the Fatimid period. The ruined section of the bridge shows damage which might be related to the AD 859 earthquake shaking. Here, we present a detailed study of the history, the status and the damage of the Cherichira aqueduct bridge using previous historic accounts and written works, a 3D laser scan model, local geological and seismological characteristics, and include results of radiocarbon dating and a timeline of events. In addition to earthquake ground motions, we consider severe flash floods on the bridge as a potential cause of the damage. We estimate the severity of such flash floods and develop a model with 18 earthquake scenarios on local reverse and strike-slip faults with magnitudes between MW 6.1 and 7.2. While a few damage patterns might be indicative of flooding, most damage can be attributed to earthquakes. It is highly probable that the earthquake in AD 859 caused enough damage to the Aghlabid bridge to render it dysfunctional; however, to resolve the question of whether another earthquake in AD 911 or 1041 caused the complete destruction of the previously retrofitted aqueduct by the Fatimids requires dating of additional sections of the bridge.

Description: Universität zu Köln (1017)

Description: Elevated soil loss and runoff rates can reduce soil fertility; therefore, soil erosion control strategies must be implemented at the hillslope and watershed scale when plantations are established and the soil is bare. In this research, we propose the use of the Vicia sativa Roth. to reduce the soil losses during the first year to control the peak of soil erosion after the plantation in tilled vineyards. To test its efficiency, rainfall simulation experiments were carried out with field surveys in The Casa Pago Gran, in the Eastern Iberian Peninsula. Our results demonstrate that soil properties (organic matter and soil bulk density) and soil erosion (soil loss, runoff and sediment concentration) were significantly different between the control plot (tilled) and with cover crops along with August 2016, 2017 and 2018 measurement periods, but not during January 2016, coinciding with the initial survey before the vetch sown. Runoff initiation was delayed in 3.7 times after 3 years (from 190 till 709 s). The runoff discharge was reduced by the Vicia sativa from 32.87 till 13.68%, the sediment concentration went down from 18.54 till 3.81 gr l−1 and the soil erosion from 3.36 to 0.29 Mg ha−1 year−1. An increase in soil bulk density was registered but did not affect the runoff generation either the soil losses, that was reduced by the plant cover. We conclude that it is necessary to include soil erosion control measures such as the use of Vicia sativa to reduce soil erosion processes during the first stages of the vineyard plantations due to the soil quality improvements and the reduction in soil and water losses.

Description: European Union Seventh Framework Programme (FP7/2007-2013)

Description: OECD (Biological Resource Management for Sustainable Agricultural Systems)

Description: Universität Trier (3163)

Description: Geospatial analysis and visualisation can be used to gain knowledge about land use structure and its changes on different spatial scales, which in turn is key to unlock the potential for sustainable land use development. This editorial provides a frame to a set of papers of the topical collection “Innovative approaches, tools and visualization techniques for analyzing land use structures and dynamics of cities and regions”, which was initiated in conjunction with the 2017 International Land Use Symposium taken place in Dresden, Germany. It first introduces current, urging land use, development and management challenges. Further on, the editorial presents the individual contributions and reflects their affiliation to the themes “Mapping and Monitoring Approaches” and “Planning, Decision Support and Participation”. Although the objectives, methods and underlying data used in the papers of this topical collection greatly vary, as pieces of a puzzle they contribute to a better analysis and understanding of current and future land use structures and dynamics of cities and regions.

Description: Leibniz-Institut für ökologische Raumentwicklung e.V. (3483)

Description: Thermal conductivity is a key parameter for many soil applications, especially for dimensioning shallow and very shallow geothermal systems based on the possible heat extraction rate and for modelling heat transfer processes around high voltage underground cables. Due to the limited purview of direct thermal conductivity measurements, for an investigation of extensive areas, usually other geophysical methods like electrical resistivity tomography measurements are applied. To derive thermal conductivity of soil from geoelectrical measurements a relation between electrical and thermal conductivity is needed. Until now only few approaches worked on a direct correlation between both conductivities. Due to the difficulties of a direct relation, within this study a modular approach of a mediate correlation between electrical and thermal conductivity was investigated. Therefore, a direct relationship between a corrected electrical conductivity and water content as well as the standard and simple thermal conductivity model of Kersten (Bull of the Univ Minnesota 28:1–227, 1949) was used. To develop this concept soil types of sand, silt loam and clay were investigated where different saturation steps and pressure loads were applied. For each configuration electrical and thermal conductivity as well as water content and bulk density was determined. To refine the results of the calculated water content a corrective factor was applied. Furthermore, bulk density as an inlet parameter of the Kersten equation was also derived based on electrical conductivity. The suggested proceeding enables the determination of thermal conductivity solely based on electrical conductivity without prior soil property information.

Description: Bundesministerium für Wirtschaft und Technologie http://dx.doi.org/10.13039/501100002765

Description: Extreme precipitation events (EPEs) are meteorological phenomena of major concern for the densely populated regions of northern and central Italy. Although statistically rare, they tend to be recurrent in autumn and share common characteristics in the large-scale dynamical evolution responsible for their generation. Past studies on EPEs have reported, as the main triggering factor, a meridionally elongated upper-level trough embedded in an incoming Rossby wave packet. In this respect, we show how the meteorological conditions leading to the devastating 1994 Piedmont flood represent a typical flow evolution for this type of extreme events. Exploiting the systematic classification of EPEs recently published by the authors and taking advantage of a new observational dataset, this article revisits the role of the large-scale flow on this and similar cases of past EPEs.

Description: Ludwig-Maximilians-Universität, Meteorologisches Institut, München, Germany

Description: In situ profiles and fixed-altitude time series of all four components of net radiation were obtained at Ny-Ålesund, Svalbard (78.9° N, 11.9° E), in the period May 04–21, 2015. Measurements were performed using adapted high-quality instrumentation classified as “secondary standard” carried by a tethered balloon system. Balloon-lifted measurements of albedo under clear-sky conditions demonstrate the local dependence on altitude and on the surface inhomogeneity of this parameter over coastal terrain of Ny-Ålesund. Depending on the surface composition within the sensor’s footprint near the coastline, the albedo over predominantly snow-covered surfaces was found to decrease to 0.548 and 0.452 at 494 m and 881 m altitude compared with 0.731 and 0.788 measured with near-surface references, respectively. Albedo profiles show an all-sky maximum at 150 m above surface level due to local surface inhomogeneity, and an averaged vertical change rate of − 0.040/100 up to 750 m aboveground level (clear sky) and − 0.034/100 m (overcast). Profiling of arctic low-level clouds reveals distinct vertical gradients in all radiative fluxes but longwave upward at the cloud top. Observed radiative cooling at the top of a partly dissolving stratus cloud with heating rates of − 40.4 to − 62.1 Kd−1 in subsequent observations is exemplified.

Description: Bundesrepublik Deutschland, Deutscher Wetterdienst, vertreten durch den Vorstand, Deutsche Meteorologische Bibliothek (4242)

Description: This paper presents a numerical simulation tool for the analysis of coupled processes related to subsurface operations. The tool combines the open-source scientific code OpenGeoSys with the reservoir simulator Eclipse enabling the coupling of thermal, hydraulic, mechanical and geochemical processes. While the coupling of multiphase flow with heat and reactive geochemical component transport has been already implemented, OpenGeoSys-Eclipse is now extended for the coupling of multiphase flow and deformation. By this, OpenGeoSys-Eclipse is capable of addressing the impact of pore pressure changes on rock stability and deformation as well as the feedback effects of geomechanical processes on multiphase flow via pore volume coupling and porosity and permeability update. The coupling is verified by several test cases of gas storage scenarios and compared with reference simulations of OpenGeoSys. The results are in good agreement regarding the general effects of geomechanical feedback on pore pressure as well as porosity and permeability changes. Differences in the results are only observed for the pore volume coupling arising from the different implementation of rock compressibility models in the two simulators. The simulations are furthermore used to investigate the relevance of addressing geomechanical feedback in numerical scenario simulations for the assessment of subsurface operations. The results show clearly, that, depending on the given storage site conditions and rock types, the feedback of deformation on pore pressure can be significant and should therefore be accounted for in the assessment.

Description: The correct choice of interface conditions and effective parameters for coupled macroscale free-flow and porous-medium models is crucial for a complete mathematical description of the problem under consideration and for accurate numerical simulation of applications. We consider single-fluid-phase systems described by the Stokes–Darcy model. Different sets of coupling conditions for this model are available. However, the choice of these conditions and effective model parameters is often arbitrary. We use large-scale lattice Boltzmann simulations to validate coupling conditions by comparison of the macroscale simulations against pore-scale resolved models. We analyse three settings (lid-driven cavity over a porous bed, infiltration problem and general filtration problem) with different geometrical configurations (channelised and staggered distributions of solid grains) and different sets of interface conditions. Effective parameters for the macroscale models (permeability tensor, boundary layer constants) are computed numerically for each geometrical configuration. Numerical simulation results demonstrate the sensitivity of the coupled Stokes–Darcy problem to the location of the sharp fluid–porous interface, the effective model parameters and the interface conditions.

Description: Bundesministerium für Bildung und Forschung https://doi.org/10.13039/501100002347

Description: Deutsche Forschungsgemeinschaft (DE)

Description: Deutsche Forschungsgemeinschaft (DE)

Description: Bundesministerium für Bildung und Forschung (DE)

Description: Geothermal district heating development has been gaining momentum in Europe with numerous deep geothermal installations and projects currently under development. With the increasing density of geothermal wells, questions related to the optimal and sustainable reservoir exploitation become more and more important. A quantitative understanding of the complex thermo-hydraulic interaction between tightly deployed geothermal wells in heterogeneous temperature and permeability fields is key for a maximum sustainable use of geothermal resources. Motivated by the geological settings of the Upper Jurassic aquifer in the Greater Munich region, we develop a computational model based on finite element analysis and gradient-free optimization to simulate groundwater flow and heat transport in hot sedimentary aquifers, and numerically investigate the optimal positioning and spacing of multi-well systems. Based on our numerical simulations, net energy production from deep geothermal reservoirs in sedimentary basins by smart geothermal multi-well arrangements provides significant amounts of energy to meet heat demand in highly urbanized regions. Our results show that taking into account heterogeneous permeability structures and a variable reservoir temperature may drastically affect the results in the optimal configuration. We demonstrate that the proposed numerical framework is able to efficiently handle generic geometrical and geological configurations, and can be thus flexibly used in the context of multi-variable optimization problems. Hence, this numerical framework can be used to assess the extractable geothermal energy from heterogeneous deep geothermal reservoirs by the optimized deployment of smart multi-well systems.

Description: Leibniz-Gemeinschaft https://doi.org/10.13039/501100001664

Description: Subsurface flow and geomechanics are often modeled with sequential approaches. This can be computationally beneficial compared with fully coupled schemes, while it requires usually compromises in numerical accuracy, at least when the sequential scheme is non-iterative. We discuss the influence of the choice of scheme on the numerical accuracy and the expected computational effort based on a comparison of a fully coupled scheme, a scheme employing a one-way coupling, and an iterative scheme using a fixed-stress split for two subsurface injection scenarios. All these schemes were implemented in the numerical simulator DuMux. This study identifies conditions of problem settings where differences due to the choice of the model approach are as important as differences in geologic features. It is shown that in particular transient and multiphase flow, effects can be causing significant deviations between non-iterative and iterative sequential schemes, which might be in the same order of magnitude as geologic uncertainty. An iterated fixed-stress split has the same numerical accuracy as a fully coupled scheme but only for a certain number of iterations which might use up the computational advantage of solving two smaller systems of equations rather than a big monolithical one.

Description: Deutsche Forschungsgemeinschaft https://doi.org/10.13039/501100001659

Description: H2020 European Research Council https://doi.org/10.13039/100010663

Description: We developed a generalized multiphase-field modeling framework for addressing the problem of brittle fracture propagation in quartz sandstones at microscopic length scale. Within this numerical approach, the grain boundaries and crack surfaces are modeled as diffuse interfaces. The two novel aspects of the model are the formulations of (I) anisotropic crack resistance in order to account for preferential cleavage planes within each randomly oriented quartz grain and (II) reduced interfacial crack resistance for incorporating lower fracture toughness along the grain boundaries that might result in intergranular crack propagation. The presented model is capable of simulating the competition between inter- and transgranular modes of fracturing based on the nature of grain boundaries, while exhibiting preferred fracturing directions within each grain. In the full parameter space, the model can serve as a powerful tool to investigate the complicated fracturing processes in heterogeneous polycrystalline rocks comprising of grains of distinct elastic properties, cleavage planes, and grain boundary attributes. We demonstrate the performance of the model through the representative numerical examples.

Description: Helmholtz-Gemeinschaft https://doi.org/10.13039/501100001656

Description: Deutsche Forschungsgemeinschaft https://doi.org/10.13039/501100001659

Description: Ministry of Baden-Wuerttemberg

Description: The given paper describes a method for automatic 3D reconstruction of bridges from cadastral footprints and airborne laser scanning point clouds. The reconstructed bridges are used to enrich 3D city models. Unlike roofs, decks of bridges are typically smooth without ridge lines or step edges. Therefore, established methods for roof reconstruction are not suitable for bridges. The standard description language for semantic city models is CityGML. This specification of the Open Geospatial Consortium assumes that surfaces are composed of planar polygons. The approximation of smooth decks by planar polygons is achieved by using a medial axis tree. Instead of the medial axis of the footprint, a modified medial axis is computed that does not consider counter bearing edges. The resulting tree represents centerline connections between all counter bearing edges and, in conjunction with filtered height values of a point cloud, serves as the basis for approximation with polygons. In addition to modeling decks, superstructures such as pylons and cables are also derived from the point cloud. For this purpose, planes carrying many superstructure points are detected using the Random Sampling Consensus Algorithm (RANSAC). Images are generated by projecting points onto these planes. Then, image processing methods are used to find connected contours that are extruded to form 3D objects. The presented method was successfully applied to all bridges of two German cities as well as to large bridges built over the Rhine River.

Description: Modelling soil erosion and sediment transport are vital to assess the impact of the flash floods. However, limited research works have studied sediment transport, especially in Egypt. This paper employs the HEC-HMS lumped hydrological model to predict the sediment load due to the flood event of 9th March 2014 in Wadi Billi, Egypt. The Modified USLE model has been used to calculate the total upland erosion, while Laursen-Copeland has been used to simulate load streams’ sediment transport potential. The Normalized Difference Vegetation Index (NDVI) has been applied over Landsat 8 image captured on 20th February 2014 using ArcMap 10.5 to determine the vegetation cover based on its spectral footprint. The resulted sedigraph showed accumulation of more than five thousand tons of sediments at the Wadi’s outlet. The results are crucial to design a suitable stormwater management system to protect the downstream urban area and to use flood water for groundwater recharge.

Description: Technische Universität Berlin (3136)

Description: https://www.hec.usace.army.mil/software/hec-hms/

Description: The differences between one classical and three state-of-the-art formulations of the mass density of humid air were quantified. Here, we present both the calculi for direct determination of the humid-air mass density employing the virial form of the thermodynamic equation of state, and a sufficiently accurate look-up-table for the quick-look determination of the humid-air mass density, which is based on the advanced Thermodynamic Equation of Seawater 2010.

Description: Leibniz-Institut für Troposphärenforschung e.V. (3489)

Description: Even though micropolar theories are widely applied for engineering applications such as the design of metamaterials, applications in the study of the Earth’s interior still remain limited and in particular in seismology. This is due to the lack of understanding of the required elastic material parameters present in the theory as well as the eigenfrequency ωr which is not observed in seismic data. By showing that the general dynamic equations of the Timoshenko’s beam is a particular case of the micropolar theory we are able to connect micropolar elastic parameters to physically measurable quantities. We then present an alternative micropolar model that, based on the same physical basis as the original model, circumvents the problem of the original eigenfrequency ωr laking in seismological data. We finally validate our model with a seismic experiment and show it is relevant to explain observed seismic dispersion curves.

Description: Westfälische Wilhelms-Universität Münster (1056)

Description: Drought is the absence or below-required supply of precipitation, runoff and or moisture for an extended time period. Modelling drought is relevant in assessing drought incidence and pattern. This study aimed to model the spatial variation and incidence of the 2018 drought in Brandenburg using GIS and remote sensing. To achieve this, we employed a Multi-Criteria Approach (MCA) by using three parameters including Precipitation, Land Surface Temperature and Normalized Difference Vegetation Index (NDVI). We acquired the precipitation data from Deutsche Wetterdienst, Land Surface Temperature and NDVI from Landsat 8 imageries on the USGS Earth Explorer. The datasets were analyzed using ArcGIS 10.7. The information from these three datasets was used as parameters in assessing drought prevalence using the MCA. The MCA was used in developing the drought model, ‘PLAN’, which was used to classify the study area into three levels/zones of drought prevalence: moderate, high and extreme drought. We went further to quantify the agricultural areas affected by drought in the study area by integrating the land use map. Results revealed that 92% of the study area was severely and highly affected by drought especially in districts of Oberhavel, Uckermark, Potsdam-Staedte, and Teltow-Flaeming. Finding also revealed that 77.54% of the total agricultural land falls within the high drought zones. We advocated for the application of drought models (such as ‘PLAN’), that incorporates flexibility (tailoring to study needs) and multi-criteria (robustness) in drought assessment. We also suggested that adaptive drought management should be championed using drought prevalence mapping.

Description: Rheinische Friedrich-Wilhelms-Universität Bonn (1040)

Description: We investigate the "macronutrient-access hypothesis", which states that the balance between stoichiometric macronutrient demand and accessible macronutrients controls nutrient assimilation by aquatic heterotrophs. Within this hypothesis, we consider bioavailable dissolved organic carbon (bDOC), reactive nitrogen (N) and reactive phosphorus (P) to be the macronutrients accessible to heterotrophic assimilation. Here, reactive N and P are the sums of dissolved inorganic N (nitrate-N, nitrite-N, ammonium-N), soluble-reactive P (SRP), and bioavailable dissolved organic N (bDON) and P (bDOP). Previous data from various freshwaters suggests this hypothesis, yet clear experimental support is missing. We assessed this hypothesis in a proof-of-concept experiment for waters from four small agricultural streams. We used seven different bDOC:reactive N and bDOC:reactive P ratios, induced by seven levels of alder leaf leachate addition. With these treatments and a stream-water specific bacterial inoculum, we conducted a 3-day experiment with three independent replicates per combination of stream water, treatment, and sampling occasion. Here, we extracted dissolved organic matter (DOM) fluorophores by measuring excitation-emission matrices with subsequent parallel factor decomposition (EEM-PARAFAC). We assessed the true bioavailability of DOC, DON, and the DOM fluorophores as the concentration difference between the beginning and end of each experiment. Subsequently, we calculated the bDOC and bDON concentrations based on the bioavailable EEM-PARAFAC fluorophores, and compared the calculated bDOC and bDON concentrations to their true bioavailability. Due to very low DOP concentrations, the DOP determination uncertainty was high, and we assumed DOP to be a negligible part of the reactive P. For bDOC and bDON, the true bioavailability measurements agreed with the same fractions calculated indirectly from bioavailable EEM-PARAFAC fluorophores (bDOC r〈sup〉2〈/sup〉 = 0.96, p 〈 0.001; bDON r〈sup〉2〈/sup〉 = 0.77, p 〈 0.001). Hence we could predict bDOC and bDON concentrations based on the EEM-PARAFAC fluorophores. The ratios of bDOC:reactive N (sum of bDON and DIN) and bDOC:reactive P (equal to SRP) exerted a strong, predictable stoichiometric control on reactive N and P uptake (R〈sup〉2〈/sup〉 = 0.80 and 0.83). To define zones of C:N:P (co-)limitation of heterotrophic assimilation, we used a novel ternary-plot approach combining our data with literature data on C:N:P ranges of bacterial biomass. Here, we found a zone of maximum reactive N uptake (C:N:P approx. 〉 114: 〈 9:1), reactive P uptake (C:N:P approx. 〉 170:21: 〈 1) and reactive N and P co-limitation of nutrient uptake (C:N:P approx. 〉 204:14:1). The “macronutrient-access hypothesis” links ecological stoichiometry and biogeochemistry, and may be of importance for nutrient uptake in many freshwater ecosystems. However, this experiment is only a starting point and this hypothesis needs to be corroborated by further experiments for more sites, by in-situ studies, and with different DOC sources.

Description: Rural and Environment Science and Analytical Services Division http://dx.doi.org/10.13039/100011310

Description: Helmholtz-Gemeinschaft http://dx.doi.org/10.13039/501100001656

Description: Helmholtz-Zentrum für Umweltforschung GmbH - UFZ (4215)

Description: The Baumgarten explosion occurred on 12 December 2017 at a gas storage site about 30 km east of Vienna, Austria. Acoustic arrivals from this accidental surface explosion were detected at dozens of stations of the AlpArray seismic network to distances up to 150 km, mainly in easterly directions. Thus it was expected that the Hungarian infrasound array PSZI located about 230 km to the east-southeast of Baumgarten would detect this acoustic wave as well. Standard progressive multichannel correlation processing and frequency-wavenumber analysis identified a signal emerging at 7:57:55 UTC from an azimuth of 296°–300° and with trace-velocity 〉 400 m/s. The extraordinarily high trace-velocity and excessive backazimuth residual, relative to the explosion site direction of 282°, however cast strong doubts on the arrival’s connection to the Baumgarten event. Accounting for the effect of non-planar geometry of the infrasound array results in a reduction of the azimuth residual by half. Additionally, 2D and 3D raytracing methods are used including the European Centre for Medium-Range Weather Forecasts (ECMWF) atmospheric model to further explain the remaining azimuth residual as well as to elucidate the large trace velocity estimates. The prevailing stratospheric winds in excess of 150 m/s are identified as the underlying cause. Including both factors the initial azimuth residual of up to 18° decreases to ~ 4°, allowing to associate the infrasound signal at PSZI with the Baumgarten event. Finally, the data from a seismic station at 30 km range is re-investigated for magnitude estimation. The local magnitude of ML 〈 1.0 explains well the scarcity of seismic observations within 50 km range, where three or four stations show signals, mainly consisting of Rg-type surface waves, but no body waves.

Description: https://www.geofon.gfz-potsdam.de/doi/network/HN

Description: https://doi.org/10.14470/UA114590

Description: https://www.geofon.gfz-potsdam.de/doi/network/SK

Description: https://www.ecmwf.int/

Description: https://www.github.com/LANL-Seismoacoustics/GeoAc

Description: Correlations between organic carbon (OC) and fine mineral particles corroborate the important role of the abundance of soil minerals with reactive surfaces to bind and increase the persistence of organic matter (OM). The storage of OM broadly consists of particulate and mineral-associated forms. Correlative studies on the impact of fine mineral soil particles on OM storage mostly combined data from differing sites potentially confounded by other environmental factors. Here, we analyzed OM storage in a soil clay content gradient of 5–37% with similar farm management and mineral composition. Throughout the clay gradient, soils contained 14 mg OC g−1 on average in the bulk soil without showing any systematic increase. Density fractionation revealed that a greater proportion of OC was stored as occluded particulate OM in the high clay soils (18–37% clay). In low clay soils (5–18% clay), the fine mineral-associated fractions had up to two times higher OC contents than high clay soils. Specific surface area measurements revealed that more mineral-associated OM was related to higher OC loading. This suggests that there is a potentially thicker accrual of more OM at the same mineral surface area within fine fractions of the low clay soils. With increasing clay content, OM storage forms contained more particulate OC and mineral-associated OC with a lower surface loading. This implies that fine mineral-associated OC storage in the studied agricultural soils was driven by thicker accrual of OM and decoupled from clay content limitations.

Description: Deutsche Forschungsgemeinschaft http://dx.doi.org/10.13039/501100001659

Description: Technische Universität München (1025)

Description: https://doi.org/10.14459/2018mp1462414

Description: The soilscapes along the southern and western coast of Crete (Greece) are dominated by coarse-grained reddish-brown slope sediments whose natural (pre-anthropogenic) configuration and properties are difficult to reconstruct due to the long history of intense land use. As a consequence, datable terrestrial sediment archives of pre-anthropogenic genesis are scarce. We present preliminary results of a study performed on an accumulation within an alluvial fan south of Stomio Bay in southwestern Crete. The studied profile is located in a small depression and is composed of a sequence of sandy to silty yellowish-brown calcareous sediments overlying coarse-grained slope sediments, including a fossil topsoil horizon. Based on macroscopic, micromorphological, geochemical, geophysical and mineralogical analysis, we interpret the fine-grained sediments to have a local aeolian origin. OSL dating indicates a final deposition phase during the early Holocene. Considering the scarcity of early Holocene terrestrial archives in Crete, the analysed profile provides valuable data for the reconstruction of landscape dynamics and paleoecological conditions as well as soil-sediment configurations during this time period. Additional research is needed to address the specific source area(s) as well as the ages of the deposition of slope sediments and formation of the fossil topsoil.

Description: Freie Universität Berlin (1008)

Description: The shale gas potential of Ediacaran and Lower Silurian shales from the Upper Yangtze platform is assessed in this study with a focus on the contributions of clay minerals and organic matter to sorption capacity. For this purpose, a multidisciplinary assessment was carried out using petrophysical, mineralogical, petrographic and geochemical methods. In terms of TOC contents (4.2%), brittle mineral contents (68.6%) and maximum gas storage capacities (0.054–0.251 mmol/g) Ediacaran shales from this study show comparable properties to other producing shale gas systems although the thermal maturity is extremely high (VRr = 3.6%). When compared to lower Silurian shales from the same region, it is evident that (1) deeper maximum burial and (2) a lack of silica-associated preservation of the pores resulted in a relatively lower mesopore volume, higher micropore volume fraction and lower overall porosity (Ediacaran shales: 1.4–4.6%; Silurian shales: 6.2–7.4%). Gas production is therefore retarded by poor interconnectivity of the pore system, which was qualitatively demonstrated by comparing experimental gas uptake kinetics. TOC content exhibits a prominent control on sorption capacity and micropore volume for both shales. However, different contributions of clay minerals to sorption capacity were identified. This can partly be attributed to different clay types but is likely also related to burial-induced recrystallisation and different origins of illite. Additionally, it was shown that variations in sorption capacity due to incorrect estimates of clay mineral contribution are in the same range as variations due to differences in thermal maturity. Article highlights Pore structure and gas storage characteristics are evaluated for the first time for Ediacaran Shales from the Upper Yangtze platform Due to a lower free gas storage capacity and diffusivity, the Ediacaran shale can be regarded as a less favorable shale gas prospect when compared to the Silurian shale Clay mineral contribution to sorption capacity is evaluated taking clay mineralogy into consideration Maturity-related changes of organic matter sorption capacity have been discussed on the basis of a compiled data set

Description: RWTH Aachen University (3131)

Description: Glacial isostatic adjustment is largely governed by the rheological properties of the Earth's mantle. Large mass redistributions in the ocean–cryosphere system and the subsequent response of the viscoelastic Earth have led to dramatic sea level changes in the past. This process is ongoing, and in order to understand and predict current and future sea level changes, the knowledge of mantle properties such as viscosity is essential. In this study, we present a method to obtain estimates of mantle viscosities by the assimilation of relative sea level rates of change into a viscoelastic model of the lithosphere and mantle. We set up a particle filter with probabilistic resampling. In an identical twin experiment, we show that mantle viscosities can be recovered in a glacial isostatic adjustment model of a simple three-layer Earth structure consisting of an elastic lithosphere and two mantle layers of different viscosity. We investigate the ensemble behaviour on different parameters in the following three set-ups: (1) global observations data set since last glacial maximum with different ensemble initialisations and observation uncertainties, (2) regional observations from Fennoscandia or Laurentide/Greenland only, and (3) limiting the observation period to 10 ka until the present. We show that the recovery is successful in all cases if the target parameter values are properly sampled by the initial ensemble probability distribution. This even includes cases in which the target viscosity values are located far in the tail of the initial ensemble probability distribution. Experiments show that the method is successful if enough near-field observations are available. This makes it work best for a period after substantial deglaciation until the present when the number of sea level indicators is relatively high.

Description: Based on the numerical weather prediction model COSMO of Germany's national meteorological service (Deutscher Wetterdienst, DWD), regional reanalysis datasets have been developed with grid spacing of up to 2 km. This development started as a fundamental research activity within the Hans-Ertel-Centre for Weather Research (HErZ) at the University of Bonn and the University of Cologne. Today, COSMO reanalyses are an established product of the DWD and have been widely used in applications on European and national German level. Successful applications of COSMO reanalyses include renewable energy assessments as well as meteorological risk estimates. The COSMO reanalysis datasets are now publicly available and provide spatio-temporal consistent data of atmospheric parameters covering both near-surface conditions and vertical profiles. This article reviews the status of the COSMO reanalyses, including evaluation results and applications. In many studies, evaluation of the COSMO reanalyses point to an overall good quality and often an added value compared to different contemporary global reanalysis datasets. We further outline current plans for the further development and application of regional reanalyses in the HErZ research group Cologne/Bonn in collaboration with the DWD.

Description: The effective radiative forcing, which includes the instantaneous forcing plus adjustments from the atmosphere and surface, has emerged as the key metric of evaluating human and natural influence on the climate. We evaluate effective radiative forcing and adjustments in 17 contemporary climate models that are participating in the Coupled Model Intercomparison Project (CMIP6) and have contributed to the Radiative Forcing Model Intercomparison Project (RFMIP). Present-day (2014) global-mean anthropogenic forcing relative to pre-industrial (1850) levels from climate models stands at 2.00 (±0.23) W m−2, comprised of 1.81 (±0.09) W m−2 from CO2, 1.08 (± 0.21) W m−2 from other well-mixed greenhouse gases, −1.01 (± 0.23) W m−2 from aerosols and −0.09 (±0.13) W m−2 from land use change. Quoted uncertainties are 1 standard deviation across model best estimates, and 90 % confidence in the reported forcings, due to internal variability, is typically within 0.1 W m−2. The majority of the remaining 0.21 W m−2 is likely to be from ozone. In most cases, the largest contributors to the spread in effective radiative forcing (ERF) is from the instantaneous radiative forcing (IRF) and from cloud responses, particularly aerosol–cloud interactions to aerosol forcing. As determined in previous studies, cancellation of tropospheric and surface adjustments means that the stratospherically adjusted radiative forcing is approximately equal to ERF for greenhouse gas forcing but not for aerosols, and consequentially, not for the anthropogenic total. The spread of aerosol forcing ranges from −0.63 to −1.37 W m−2, exhibiting a less negative mean and narrower range compared to 10 CMIP5 models. The spread in 4×CO2 forcing has also narrowed in CMIP6 compared to 13 CMIP5 models. Aerosol forcing is uncorrelated with climate sensitivity. Therefore, there is no evidence to suggest that the increasing spread in climate sensitivity in CMIP6 models, particularly related to high-sensitivity models, is a consequence of a stronger negative present-day aerosol forcing and little evidence that modelling groups are systematically tuning climate sensitivity or aerosol forcing to recreate observed historical warming.

Description: We present the result of the third Marine Ice Sheet Intercomparison project, MISMIP+. MISMIP+ is intended to be a test of ice flow models which include fast sliding marine ice streams and floating ice shelves and in particular a treatment of viscous stress that is sufficient for buttressing, where upstream ice flow is restrained by a downstream ice shelf. A set of idealized experiments test the models in circumstances where buttressing contributes to a stable steady state, and where a reduction in that buttressing causes ice stream acceleration, thinning, and grounding line retreat. We find that the most important distinction between models in this particular type of simulation is in the treatment of sliding at the bed, with other distinctions – notably the difference between the simpler and more complete treatments of englacial stress, but also the differences between numerical methods – taking a secondary role.

Description: Varved lake sediments provide long climatic records with high temporal resolution and low associated age uncertainty. Robust and detailed comparison of well-dated and annually laminated sediment records is crucial for reconstructing abrupt and regionally time-transgressive changes as well as validation of spatial and temporal trajectories of past climatic changes. The VARved sediments DAtabase (VARDA) presented here is the first data compilation for varve chronologies and associated palaeoclimatic proxy records. The current version 1.0 allows detailed comparison of published varve records from 95 lakes. VARDA is freely accessible and was created to assess outputs from climate models with high-resolution terrestrial palaeoclimatic proxies. VARDA additionally provides a technical environment that enables to explore the database of varved lake sediments using a connected data-model and can generate a state-of-the-art graphic representation of multi-site comparison. This allows to reassess existing chronologies and tephra events to synchronize and compare even distant varved lake records. Furthermore, the present version of VARDA permits to explore varve thickness data. In this paper, we report in detail on the data mining and compilation strategies for the identification of varved lakes and assimilation of high-resolution chronologies as well as the technical infrastructure of the database. Additional paleoclimate proxy data will be provided in forthcoming updates. The VARDA graph-database and user interface can be accessed online at https://varve.gfz-potsdam.de, all datasets of version 1.0 are available at http://doi.org/10.5880/GFZ.4.3.2019.003 (Ramisch et al., 2019).

Description: Various observational estimates indicate growing mass loss at Antarctica's margins but also heavier precipitation across the continent. In the future, heavier precipitation fallen on Antarctica will counteract any stronger iceberg discharge and increased basal melting of floating ice shelves driven by a warming ocean. Here, we use from nine CMIP5 models future projections, ranging from strong mitigation efforts to business-as-usual, to run an ensemble of ice-sheet simulations. We test, how the precipitation boundary condition determines Antarctica's sea-level contribution. The spatial and temporal varying climate forcings drive ice-sheet simulations. Hence, our ensemble inherits all spatial and temporal climate patterns, which is in contrast to a spatial mean forcing. Regardless of the applied boundary condition and forcing, some areas will lose ice in the future, such as the glaciers from the West Antarctic Ice Sheet draining into the Amundsen Sea. In general the simulated ice-sheet thickness grows in a broad marginal strip, where incoming storms deliver topographically controlled precipitation. This strip shows the largest ice thickness differences between the applied precipitation boundary conditions too. On average Antarctica's ice mass shrinks for all future scenarios if the precipitation is scaled by the spatial temperature anomalies coming from the CMIP5 models. In this approach, we use the relative precipitation increment per degree warming as invariant scaling constant. In contrast, Antarctica gains mass in our simulations if we apply the simulated precipitation anomalies of the CMIP5 models directly. Here, the scaling factors show a distinct spatial pattern across Antarctica. Furthermore, the diagnosed mean scaling across all considered climate forcings is larger than the values deduced from ice cores. In general, the scaling is higher across the East Antarctic Ice Sheet, lower across the West Antarctic Ice Sheet, and lowest around the Siple Coast. The latter is located on the east side of the Ross Ice Shelf.

Description: A finely laminated lake sediment record with a basal age of 11,619 ± 603 years BP was retrieved from Lake Chatyr Kol (Kyrgyz Republic). Microfacies analysis reveals the presence of seasonal laminae (varves) from the sediment basis to ~ 360 ± 40 years BP. The Chatvd19 floating varve chronology covers the time span from 360 ± 40 years BP to the base and relies on replicate varve counts on overlapping petrographic thin sections with an uncertainty of ± 5 %. The uppermost non-varved interval was chronologically constrained by 210Pb and 137Cs γ-spectrometry and interpolation based on varve thickness measurements of adjacent varved intervals with an assumed uncertainty of 10 %. Six varve types were distinguished, are described in detail and show a changing predominance of clastic-organic, clastic-calcitic or -aragonitic, calcitic-clastic, organic-clastic and clastic-diatom varves throughout the Holocene. Variations in varve thickness and the number and composition of seasonal sublayers are attributed to 1) changes in the amount of summer or winter/spring precipitation affecting local runoff and erosion and/or to 2) evaporative conditions during summer. Radiocarbon dating of bulk organic matter, daphnia remains, aquatic plant remains and Ruppia maritima seeds reveal reservoir ages with a clear decreasing trend up core from ~ 6,150 years in the early Holocene, to ~ 3,000 years in the mid-Holocene, to ~ 1,000 years and less in the late Holocene and modern times. In contrast, two radiocarbon dates from terrestrial plant remains are in good agreement with the varve-based chronology.

Description: Plankton is a massive and phylogenetically diverse group of thousands of prokaryotes, protists (unicellular eukaryotic organisms), and metazoans (multicellular eukaryotic organisms; Fig. 1). Plankton functional diversity is at the core of various ecological processes, including productivity, carbon cycling and sequestration, nutrient cycling (Falkowski 2012), interspecies interactions, and food web dynamics and structure (D'Alelio et al. 2016). Through these functions, plankton play a critical role in the health of the coastal and open ocean and provide essential ecosystem services. Yet, at present, our understanding of plankton dynamics is insufficient to project how climate change and other human-driven impacts affect the functional diversity of plankton. That limits our ability to predict how critical ecosystem services will change in the future and develop strategies to adapt to these changes.

Description: Ice flow models of the Antarctic ice sheet are commonly used to simulate its future evolution in response to different climate scenarios and inform on the mass loss that would contribute to future sea level rise. However, there is currently no consensus on estimated the future mass balance of the ice sheet, primarily because of differences in the representation of physical processes and the forcings employed. This study presents results from 18 simulations from 15 international groups focusing on the evolution of the Antarctic ice sheet during the period 2015–2100, forced with different scenarios from the Coupled Model Intercomparison Project Phase 5 (CMIP5) representative of the spread in climate model results. The contribution of the Antarctic ice sheet in response to increased warming during this period varies between −7.8 and 30.0 cm of Sea Level Equivalent (SLE). The evolution of the West Antarctic Ice Sheet varies widely among models, with an overall mass loss up to 21.0 cm SLE in response to changes in oceanic conditions. East Antarctica mass change varies between −6.5 and 16.5 cm SLE, with a significant increase in surface mass balance outweighing the increased ice discharge under most RCP 8.5 scenario forcings. The inclusion of ice shelf collapse, here assumed to be caused by large amounts of liquid water ponding at the surface of ice shelves, yields an additional mass loss of 8 mm compared to simulations without ice shelf collapse. The largest sources of uncertainty come from the ocean-induced melt rates, the calibration of these melt rates based on oceanic conditions taken outside of ice shelf cavities and the ice sheet dynamic response to these oceanic changes. Results under RCP 2.6 scenario based on two CMIP5 AOGCMs show an overall mass loss of 10 mm SLE compared to simulations done under present-day conditions, with limited mass gain in East Antarctica.

Description: The Net-Zero-2050 cluster aims for a national roadmap for net zero CO2 emissions by 2050, including integrated scenario analyses and negative emission technology assessment (see fact sheet Net-Zero-2050 Structure Project 1). This national target to substantially reduce national CO2 emissions by 2050 stems from the objective to comply with the global long-term temperature goal of well below 2°C of the Paris Agreement (UNFCCC, 2015). Within the cluster it is therefore important to decide on an approach for deriving a national remaining carbon budget from global emissions trajectories in agreement with the Paris Climate Agreement’s longterm temperature goal (UNFCCC, 2015). Allocating national carbon budgets is a balance of environmental effectiveness, equity, national capacity and ability, political feasibility, economic efficiency and technical requirements (Gignac and Matthews, 2015; Höhne et al., 2003; 2014). Given Germany’s capacity and abilities, we decided to follow a sustainable growth trajectory with a convergence phase to equal-per-capita CO2 emissions by 2035, and a net zero CO2 emissions trajectory after 2050 until the end of the century. This approach leads to a remaining Germany CO2 budget of 9 GtCO2 (from 1st January 2018 to 2050 and 2100), which we propose to be used across the Net-Zero-2050 cluster. The remaining carbon budget will serve as a target to be used in all work packages in a concerted way, either qualitatively or quantitatively, and in accordance with other work packages (see also fact sheet Net-Zero-2050 Energy Scenario Approach). The calculated budget is at the lower end of the national budget if allocated by the grandfathering approach (emissions are allocated with respect to today’s emissions shares: 5.5-13.1 GtCO2), but slightly higher than the highest estimate of an equal-per-capita remaining carbon budget (emissions are allocated with respect to Germany’s share of the global population: 3.5-8.4 GtCO2) The 9 GtCO2 national remaining CO2 budget, 6.9 GtCO2 from 1st January 2021, will need to be broken down by category (e.g. energy, land use, industrial processes, and man-made sinks and sources; see Gap Analysis Report) in order to provide a consistent approach across work packages.

Description: Doppler wind lidars (DWLs) have increasingly been used over the last decade to derive the mean wind in the atmospheric boundary layer. DWLs allow the determination of wind vector profiles with high vertical resolution and provide an alternative to classic meteorological tower observations. They also receive signals from altitudes higher than a tower and can be set up flexibly in any power-supplied location. In this work, we address the question of whether and how wind gusts can be derived from DWL observations. The characterization of wind gusts is one central goal of the Field Experiment on Sub-Mesoscale Spatio-Temporal Variability in Lindenberg (FESSTVaL). Obtaining wind gusts from a DWL is not trivial because a monostatic DWL provides only a radial velocity per line of sight, i.e., only one component of a three-dimensional vector, and measurements in at least three linearly independent directions are required to derive the wind vector. Performing them sequentially limits the achievable time resolution, while wind gusts are short-lived phenomena. This study compares different DWL configurations in terms of their potential to derive wind gusts. For this purpose, we develop a new wind retrieval method that is applicable to different scanning configurations and various time resolutions. We test eight configurations with StreamLine DWL systems from HALO Photonics and evaluate gust peaks and mean wind over 10 min at 90 m a.g.l. against a sonic anemometer at the meteorological tower in Falkenberg, Germany. The best-performing configuration for retrieving wind gusts proves to be a fast continuous scanning mode (CSM) that completes a full observation cycle within 3.4 s. During this time interval, about 11 radial Doppler velocities are measured, which are then used to retrieve single gusts. The fast CSM configuration was successfully operated over a 3-month period in summer 2020. The CSM paired with our new retrieval technique provides gust peaks that compare well to classic sonic anemometer measurements from the meteorological tower.

Description: A challenge of an energy system that nowadays more strongly depends on wind power generation is the spatial and temporal variability in winds. Nocturnal low-level jets (NLLJs) are typical wind phenomena defined as a maximum in the vertical profile of the horizontal wind speed. A NLLJ has typical core heights of 50–500 m a.g.l. (above ground level), which is in the height range of most modern wind turbines. This study presents NLLJ analyses based on new observations from Doppler wind lidars. The aim is to characterize the temporal and spatial variability in NLLJs on the mesoscale and to quantify their impacts on wind power generation. The data were collected during the Field Experiment on Submesoscale Spatio-Temporal Variability (FESSTVaL) campaign from June to August 2020 in Lindenberg and Falkenberg (Germany), located at about 6 km from each other. Both sites have seen NLLJs in about 70 % of the nights with half of them lasting for more than 3 h. Events longer than 6 h occurred more often simultaneously at both sites than shorter events, indicating the mesoscale character of very long NLLJs. Very short NLLJs of less than 1 h occurred more often in Lindenberg than Falkenberg, indicating more local influences on the wind profile. We discussed different meteorological mechanisms for NLLJ formation and linked NLLJ occurrences to synoptic weather patterns. There were positive and negative impacts of NLLJs on wind power that we quantified based on the observational data. NLLJs increased the mean power production by up to 80 % and were responsible for about 25 % of the power potential during the campaign. However, the stronger shear in the rotor layer during NLLJs can also have negative impacts. The impacts of NLLJs on wind power production depended on the relative height between the wind turbine and the core of the NLLJ. For instance, the mean increase in the estimated power production during NLLJ events was about 30 % higher for a turbine at 135 m a.g.l. compared to one at 94 m a.g.l. Our results imply that long NLLJs have an overall stronger impact on the total power production, while short events are primarily relevant as drivers for power ramps.

Description: The reaction between ozone and iodide at the sea surface is now known to be an important part of atmospheric ozone cycling, causing ozone deposition and the release of ozone-depleting reactive iodine to the atmosphere. The importance of this reaction is reflected by its inclusion in chemical transport models (CTMs). Such models depend on accurate sea surface iodide fields, but measurements are spatially and temporally limited. Hence, the ability to predict current and future sea surface iodide fields, i.e. sea surface iodide concentration on a narrow global grid, requires the development of process-based models. These models require a thorough understanding of the key processes that control sea surface iodide. The aim of this study was to explore if there are common features of iodate-to-iodide reduction amongst diverse marine phytoplankton in order to develop models that focus on sea surface iodine and iodine release to the troposphere. In order to achieve this, rates and patterns of changes in inorganic iodine speciation were determined in 10 phytoplankton cultures grown at ambient iodate concentrations. Where possible these data were analysed alongside results from previous studies. Iodate loss and some iodide production were observed in all cultures studied, confirming that this is a widespread feature amongst marine phytoplankton. We found no significant difference in log-phase, cell-normalised iodide production rates between key phytoplankton groups (diatoms, prymnesiophytes including coccolithophores and phaeocystales), suggesting that a phytoplankton functional type (PFT) approach would not be appropriate for building an ocean iodine cycling model. Iodate loss was greater than iodide formation in the majority of the cultures studied, indicating the presence of an as-yet-unidentified “missing iodine” fraction. Iodide yield at the end of the experiment was significantly greater in cultures that had reached a later senescence stage. This suggests that models should incorporate a lag between peak phytoplankton biomass and maximum iodide production and that cell mortality terms in biogeochemical models could be used to parameterise iodide production.

Description: Nitric oxide (NO) is a short-lived intermediate of the oceanic nitrogen cycle. However, our knowledge about its production and consumption pathways in oceanic environments is rudimentary. In order to decipher the major factors affecting NO photochemical production, we irradiated several artificial seawater samples as well as 31 natural surface seawater samples in laboratory experiments. The seawater samples were collected during a cruise to the western tropical North Pacific Ocean (WTNP, a N-S section from 36 to 2 degrees N along 146 to 143 degrees E with 6 and 12 stations, respectively, and a W-E section from 137 to 161 degrees E along the Equator with 13 stations) from November 2015 to January 2016. NO photoproduction rates from dissolved nitrite in artificial seawater showed increasing trends with decreasing pH, increasing temperature, and increasing salinity. In contrast, NO photoproduction rates (average: 0.5 +/- 0.2 x 10(-12) mol L-1 s(-1)) in the natural seawater samples from the WTNP did not show any correlations with pH, water temperature, salinity, or dissolved inorganic nitrite concentrations. The flux induced by NO photoproduction in the WTNP (average: 13 x 10(-12) mol M-2 S-1) was significantly larger than the NO air-sea flux density (average: 1.8 x 10(-12) Mol M-2 S-1), indicating a further NO loss process in the surface layer.

Description: Variations of the solar spectral irradiance (SSI) with the 11-year sunspot cycle have been shown to have a significant impact on temperatures and the mixing ratios of atmospheric constituents in the stratosphere and mesosphere. Uncertainties in modelling the effects of SSI variations arise from uncertainties in the empirical models reconstructing the prescribed SSI data set as well as from uncertainties in the chemistry-climate model (CCM) formulation. In this study CCM simulations with the ECHAM MESSy Atmospheric Chemistry (EMAC) model and the Community Earth System Model 1 (CESM1) – Whole Atmosphere Chemistry Climate Model (WACCM) have been performed to quantify the uncertainties of the solar responses in chemistry and dynamics that are due to the usage of five different SSI data sets or the two CCMs. We apply a two-way analysis of variance (ANOVA) to separate the influence of the SSI data sets and the CCMs on the variability of the solar response in shortwave heating rates, temperature and ozone. The ANOVA identifies the SSI data set with the strongest influence on the variability of the solar signal in shortwave heating rates in the upper mesosphere and in the upper stratosphere/lower mesosphere. The strongest influence on the variability of the solar signal in ozone and temperature is identified in the upper stratosphere/lower mesosphere. The largest influence of the CCMs on variability of the solar responses can be identified in the upper mesosphere. The solar response in the lower stratosphere also depends on the CCM used, especially in the tropics and northern hemispheric subtropics and mid latitudes, where the model dynamics modulate the solar responses.

Description: Coastal areas contribute significantly to the emissions of methane (CH4) from the ocean. In order to decipher its temporal variability in the whole water column, dissolved CH4 was measured on a monthly basis at the Boknis Eck Time-series Station (BE) located in the Eckernförde Bay (SW Baltic Sea) from 2006 to 2017. BE has a water depth of about 28 m and dissolved CH4 was measured at six water depths ranging from 0 to 25 m. In general CH4 concentrations increased with depth, indicating a sedimentary release of CH4. Pronounced enhancement of the CH4 concentrations in the bottom layer (15–25 m) was found during February, May–June and October. CH4 was not correlated with Chlorophyll a or O2 over the measurement period. Unusually high CH4 concentrations (of up to 696 nM) were sporadically observed in the upper layer (0–10 m) (e.g. in November 2013 and December 2014) and were coinciding with Major Baltic Inflow (MBI) events. Surface CH4 concentrations were always supersaturated throughout the monitoring period, indicating that the Eckernförde Bay is an intense but highly variable source of atmospheric CH4. We did not detect significant temporal trends in CH4 concentrations or emissions, despite of ongoing environmental changes such as warming and deoxygenation in the Eckernförde Bay. Overall, the CH4 variability at BE is driven by a complex interplay of various biological and physical processes.

Description: The extracellular concentration of H2O2 in surface aquatic environments is controlled by a balance between photochemical production and the microbial synthesis of catalase and peroxidase enzymes to remove H2O2 from solution. In any kind of incubation experiment, the formation rates and equilibrium concentrations of reactive oxygen species (ROSs) such as H2O2 may be sensitive to both the experiment design, particularly to the regulation of incident light, and the abundance of different microbial groups, as both cellular H2O2 production and catalase–peroxidase enzyme production rates differ between species. Whilst there are extensive measurements of photochemical H2O2 formation rates and the distribution of H2O2 in the marine environment, it is poorly constrained how different microbial groups affect extracellular H2O2 concentrations, how comparable extracellular H2O2 concentrations within large-scale incubation experiments are to those observed in the surface-mixed layer, and to what extent a mismatch with environmentally relevant concentrations of ROS in incubations could influence biological processes differently to what would be observed in nature. Here we show that both experiment design and bacterial abundance consistently exert control on extracellular H2O2 concentrations across a range of incubation experiments in diverse marine environments. During four large-scale (〉1000 L) mesocosm experiments (in Gran Canaria, the Mediterranean, Patagonia and Svalbard) most experimental factors appeared to exert only minor, or no, direct effect on H2O2 concentrations. For example, in three of four experiments where pH was manipulated to 0.4–0.5 below ambient pH, no significant change was evident in extracellular H2O2 concentrations relative to controls. An influence was sometimes inferred from zooplankton density, but not consistently between different incubation experiments, and no change in H2O2 was evident in controlled experiments using different densities of the copepod Calanus finmarchicus grazing on the diatom Skeletonema costatum (〈1 % change in [H2O2] comparing copepod densities from 1 to 10 L−1). Instead, the changes in H2O2 concentration contrasting high- and low-zooplankton incubations appeared to arise from the resulting changes in bacterial activity. The correlation between bacterial abundance and extracellular H2O2 was stronger in some incubations than others (R2 range 0.09 to 0.55), yet high bacterial densities were consistently associated with low H2O2. Nonetheless, the main control on H2O2 concentrations during incubation experiments relative to those in ambient, unenclosed waters was the regulation of incident light. In an open (lidless) mesocosm experiment in Gran Canaria, H2O2 was persistently elevated (2–6-fold) above ambient concentrations; whereas using closed high-density polyethylene mesocosms in Crete, Svalbard and Patagonia H2O2 within incubations was always reduced (median 10 %–90 %) relative to ambient waters.

Description: Cold-water corals (CWCs) constitute important deep-water ecosystems that are under increasing environmental pressure due to ocean acidification and global warming. The sensitivity of these deep-water ecosystems to environmental change is demonstrated by abundant paleorecords drilled through CWC mounds that reveal characteristic alterations between rapid formation and dormant or erosive phases. Previous studies have identified several central parameters for driving or inhibiting CWC growth such as food supply, oxygenation, and the carbon saturation state of bottom water, yet there are still large uncertainties about the relative importance of the different environmental parameters. To advance this debate we have performed a multiproxy study on a sediment core retrieved from the 25 m high Bowie Mound, located at 866 m water depth on the continental slope off southeastern Brazil, a structure built up mainly by the CWC Solenosmilia variabilis. Our results indicate a multifactorial control on CWC growth at Bowie Mound during the past ∼ 160 kyr, which reveals distinct formation pulses during northern high-latitude glacial cold events (Heinrich stadials, HSs) largely associated with anomalously strong monsoonal rainfall over the continent. The ensuing enhanced runoff elevated the terrigenous nutrient and organic-matter supply to the continental margin and likely boosted marine productivity. The dispersal of food particles towards the CWC colonies during HSs was facilitated by the highly dynamic hydraulic conditions along the continental slope that prevailed throughout glacial periods. These conditions caused the emplacement of a pronounced nepheloid layer above Bowie Mound, thereby aiding the concentration and along-slope dispersal of organic matter. Our study thus emphasizes the impact of continental climate variability on a highly vulnerable deep-marine ecosystem.

Description: Variations in the solar spectral irradiance (SSI) with the 11-year sunspot cycle have been shown to have a significant impact on temperatures and the mixing ratios of atmospheric constituents in the stratosphere and mesosphere. Uncertainties in modelling the effects of SSI variations arise from uncertainties in the empirical models reconstructing the prescribed SSI data set as well as from uncertainties in the chemistry–climate model (CCM) formulation. In this study CCM simulations with the ECHAM/MESSy Atmospheric Chemistry (EMAC) model and the Community Earth System Model 1 (CESM1)–Whole Atmosphere Chemistry Climate Model (WACCM) have been performed to quantify the uncertainties of the solar responses in chemistry and dynamics that are due to the usage of five different SSI data sets or the two CCMs. We apply a two-way analysis of variance (ANOVA) to separate the influence of the SSI data sets and the CCMs on the variability of the solar response in shortwave heating rates, temperature, and ozone. The solar response is derived from climatological differences of time slice simulations prescribing SSI for the solar maximum in 1989 and near the solar minimum in 1994. The SSI values for the solar maximum of each SSI data set are created by adding the SSI differences between November 1994 and November 1989 to a common SSI reference spectrum for near-solar-minimum conditions based on ATLAS-3 (Atmospheric Laboratory of Applications and Science-3). The ANOVA identifies the SSI data set with the strongest influence on the variability of the solar response in shortwave heating rates in the upper mesosphere and in the upper stratosphere–lower mesosphere. The strongest influence on the variability of the solar response in ozone and temperature is identified in the upper stratosphere–lower mesosphere. However, in the region of the largest ozone mixing ratio, in the stratosphere from 50 to 10 hPa, the SSI data sets do not contribute much to the variability of the solar response when the Spectral And Total Irradiance REconstructions-T (SATIRE-T) SSI data set is omitted. The largest influence of the CCMs on variability of the solar responses can be identified in the upper mesosphere. The solar response in the lower stratosphere also depends on the CCM used, especially in the tropics and northern hemispheric subtropics and mid-latitudes, where the model dynamics modulate the solar responses. Apart from the upper mesosphere, there are also regions where the largest fraction of the variability of the solar response is explained by randomness, especially for the solar response in temperature.

Description: In Kiel, in the north of Germany, marine research is rooted in a lively research community hosted mainly at Kiel University and the GEOMAR Helmholtz Centre. While the ratio of women and men is more or less balanced on all qualification levels with mainly nonpermanent junior positions, women are generally underrepresented in leading research positions. The problem of gender imbalance and inequality has been well-known for a long time. Especially in the last decade, however, manifold efforts were initiated to improve gender equality on a political and institutional level as well as within the research community itself. In our article we focus on the gender equality activities of the two large externally funded marine sciences research alliances: the Cluster of Excellence “The Future Ocean” and the Collaborative Research Centre 754 “Climate–Biogeochemistry Interactions in the Tropical Ocean”. For about a decade they offered both financial provisions and a structural framework to tackle the problem of women's underrepresentation in science and came up with innovative measures. In the following case study, we not only introduce the situation of women in marine sciences in Kiel and the structural arrangement to improve gender equality in general, but we also discuss three specific measures developed within the two collaborative research projects in detail: (i) the mentoring program via:mento_ocean for female postdocs, (ii) hiring policies integrating a gender quota for recruiting postdoctoral researchers and (iii) a code of conduct. Based on these best-practice examples we can show that progress towards gender equality has been made despite some obstacles faced when implementing the measures. This was especially the case for attracting female researchers to work in Kiel marine sciences and bringing the relevance of the topic to the surface of debates within the community. Looking at gender equality activities from a managerial point of view, we conclude from the situation in Kiel, where external funding for both research alliances ended in 2019, that even time-bound activities can initiate change. Initiatives developed by the marine sciences community were taken up by other research groups and inspired new activities at the level of the institutions

Description: The Zero Emissions Commitment (ZEC) is the change in global mean temperature expected to occur following the cessation of net CO2 emissions and as such is a critical parameter for calculating the remaining carbon budget. The Zero Emissions Commitment Model Intercomparison Project (ZECMIP) was established to gain a better understanding of the potential magnitude and sign of ZEC, in addition to the processes that underlie this metric. A total of 18 Earth system models of both full and intermediate complexity participated in ZECMIP. All models conducted an experiment where atmospheric CO2 concentration increases exponentially until 1000 PgC has been emitted. Thereafter emissions are set to zero and models are configured to allow free evolution of atmospheric CO2 concentration. Many models conducted additional second-priority simulations with different cumulative emission totals and an alternative idealized emissions pathway with a gradual transition to zero emissions. The inter-model range of ZEC 50 years after emissions cease for the 1000 PgC experiment is −0.36 to 0.29 ∘C, with a model ensemble mean of −0.07 ∘C, median of −0.05 ∘C, and standard deviation of 0.19 ∘C. Models exhibit a wide variety of behaviours after emissions cease, with some models continuing to warm for decades to millennia and others cooling substantially. Analysis shows that both the carbon uptake by the ocean and the terrestrial biosphere are important for counteracting the warming effect from the reduction in ocean heat uptake in the decades after emissions cease. This warming effect is difficult to constrain due to high uncertainty in the efficacy of ocean heat uptake. Overall, the most likely value of ZEC on multi-decadal timescales is close to zero, consistent with previous model experiments and simple theory.

Description: Many coastal oceans experience not only increased loads of nutrients but also changes in the stoichiometry of nutrient supply. Excess supply of nitrogen and stable or decreased supply of silicon lower silicon to nitrogen (Si:N) ratios, which may decrease diatom proportion in phytoplankton. To examine how Si:N ratios affect plankton community composition and food web structure, we performed a mesocosm experiment where we manipulated Si:N ratios and copepod abundance in a Baltic Sea plankton community. In high Si:N treatments, diatoms dominated. Some of them were likely spared from grazing unexpectedly resulting in higher diatom biomass under high copepod grazing. With declining Si:N ratios, dinoflagellates became more abundant under low and picoplankton under high copepod grazing. This altered plankton food web structure: under high Si:N ratios, edible diatoms were directly accessible food for copepods, while under low Si:N ratios, microzooplankton and phago-mixotrophs (mixoplankton) were a more important food source for mesograzers. The response of copepods to changes in the phytoplankton community was complex and copepod density-dependent. We suggest that declining Si:N ratios favor microzoo- and mixoplankton leading to increased complexity of planktonic food webs. Consequences on higher trophic levels will, however, likely be moderated by edibility, nutritional value or toxicity of dominant phytoplankton species.

Description: Quantification and attribution of the food web changes associated with the invasion of non-indigenous species in the marine realm often remain a challenge. One of the pelagic non-indigenous species of concern in the recent history of aquatic bioinvasions is the predatory cladoceran Cercopagis pengoi, which invaded the Baltic Sea in the early 1990s. While several studies have reported immediate declines in abundances of its potential prey, the long-term effects of C. pengoi on the food webs remain to be examined. Based on the long-term time series (1968–2018) in the Gulf of Riga (Baltic Sea), we found significant declines in abundance of the cladoceran Pleopis spp. and copepod Eurytemora affinis by 90 and 80%, respectively, are associated with the invasion of C. pengoi as well as significant alterations in seasonal abundance patterns of Pleopis spp., E. affinis and cladoceran Bosmina spp. The invasion of the non-indigenous predator has led to the changed prey abundance–temperature relationships. Special caution was taken in data preprocessing, to minimize the likelihood that observed changes in the zooplankton prey could be associated with factors other than the invasion of C. pengoi, such as temperature and storminess.

Description: The last few decades have seen dramatic changes in the hydrography and biogeochemistry of the Mediterranean Sea. The complex bathymetry and highly variable spatial and temporal scales of atmospheric forcing, convective and ventilation processes contribute to generate complex and unsteady circulation patterns and significant variability in biogeochemical systems. Part of the variability of this system can be influenced by anthropogenic contributions. Consequently, it is necessary to document details and to understand trends in place to better relate the observed processes and to possibly predict the consequences of these changes. In this context we report data from an oceanographic cruise in the Mediterranean Sea on the German research vessel Maria S. Merian (MSM72) in March 2018. The main objective of the cruise was to contribute to the understanding of long-term changes and trends in physical and biogeochemical parameters, such as the anthropogenic carbon uptake and to further assess the hydrographical situation after the major climatological shifts in the eastern and western part of the basin, known as the Eastern and Western Mediterranean Transients. During the cruise, multidisciplinary measurements were conducted on a predominantly zonal section throughout the Mediterranean Sea, contributing to the Med-SHIP and GO-SHIP long-term repeat cruise section that is conducted at regular intervals in the Mediterranean Sea to observe changes and impacts on physical and biogeochemical variables.

Description: Decreasing concentrations of dissolved oxygen in the ocean are considered one of the main threats to marine ecosystems as they jeopardize the growth of higher organisms. They also alter the marine nitrogen cycle, which is strongly bound to the carbon cycle and climate. While higher organisms in general start to suffer from oxygen concentrations 〈 ∼ 63 µM (hypoxia), the marine nitrogen cycle responds to oxygen concentration below a threshold of about 20 µM (microbial hypoxia), whereas anoxic processes dominate the nitrogen cycle at oxygen concentrations of 〈 ∼ 0.05 µM (functional anoxia). The Arabian Sea and the Bay of Bengal are home to approximately 21 % of the total volume of ocean waters revealing microbial hypoxia. While in the Arabian Sea this oxygen minimum zone (OMZ) is also functionally anoxic, the Bay of Bengal OMZ seems to be on the verge of becoming so. Even though there are a few isolated reports on the occurrence of anoxia prior to 1960, anoxic events have so far not been reported from the open northern Indian Ocean (i.e., other than on shelves) during the last 60 years. Maintenance of functional anoxia in the Arabian Sea OMZ with oxygen concentrations ranging between 〉 0 and ∼ 0.05 µM is highly extraordinary considering that the monsoon reverses the surface ocean circulation twice a year and turns vast areas of the Arabian Sea from an oligotrophic oceanic desert into one of the most productive regions of the oceans within a few weeks. Thus, the comparably low variability of oxygen concentration in the OMZ implies stable balances between the physical oxygen supply and the biological oxygen consumption, which includes negative feedback mechanisms such as reducing oxygen consumption at decreasing oxygen concentrations (e.g., reduced respiration). Lower biological oxygen consumption is also assumed to be responsible for a less intense OMZ in the Bay of Bengal. According to numerical model results, a decreasing physical oxygen supply via the inflow of water masses from the south intensified the Arabian Sea OMZ during the last 6000 years, whereas a reduced oxygen supply via the inflow of Persian Gulf Water from the north intensifies the OMZ today in response to global warming. The first is supported by data derived from the sedimentary records, and the latter concurs with observations of decreasing oxygen concentrations and a spreading of functional anoxia during the last decades in the Arabian Sea. In the Arabian Sea decreasing oxygen concentrations seem to have initiated a regime shift within the pelagic ecosystem structure, and this trend is also seen in benthic ecosystems. Consequences for biogeochemical cycles are as yet unknown, which, in addition to the poor representation of mesoscale features in global Earth system models, reduces the reliability of estimates of the future OMZ development in the northern Indian Ocean.

Description: In the current era of rapid climate change, accurate characterization of climate-relevant gas dynamics-namely production, consumption, and net emissions-is required for all biomes, especially those ecosystems most susceptible to the impact of change. Marine environments include regions that act as net sources or sinks for numerous climateactive trace gases including methane (CH4) and nitrous oxide (N2O). The temporal and spatial distributions of CH4 and N2O are controlled by the interaction of complex biogeochemical and physical processes. To evaluate and quantify how these mechanisms affect marine CH4 and N2O cycling requires a combination of traditional scientific disciplines including oceanography, microbiology, and numerical modeling. Fundamental to these efforts is ensuring that the datasets produced by independent scientists are comparable and interoperable. Equally critical is transparent communication within the research community about the technical improvements required to increase our collective understanding of marine CH4 and N2O. A workshop sponsored by Ocean Carbon and Biogeochemistry (OCB) was organized to enhance dialogue and collaborations pertaining to marine CH4 and N2O. Here, we summarize the outcomes from the workshop to describe the challenges and opportunities for near-future CH4 and N2O research in the marine environment.

Description: It is widely accepted that orbital variations are responsible for the generation of glacial cycles during the late Pleistocene. However, the relative contributions of the orbital forcing compared to CO2 variations and other feedback mechanisms causing the waxing and waning of ice sheets have not been fully understood. Testing theories of ice ages beyond statistical inferences, requires numerical modeling experiments that capture key features of glacial transitions. Here, we focus on the glacial buildup from Marine Isotope Stage (MIS) 7 to 6 covering the period from 240 to 170 ka (ka: thousand years before present). This transition from interglacial to glacial conditions includes one of the fastest Pleistocene glaciation–deglaciation events, which occurred during MIS 7e–7d–7c (236–218 ka). Using a newly developed three-dimensional coupled atmosphere–ocean–vegetation–ice sheet model (LOVECLIP), we simulate the transient evolution of Northern Hemisphere and Southern Hemisphere ice sheets during the MIS 7–6 period in response to orbital and greenhouse gas forcing. For a range of model parameters, the simulations capture the evolution of global ice volume well within the range of reconstructions. Over the MIS 7–6 period, it is demonstrated that glacial inceptions are more sensitive to orbital variations, whereas terminations from deep glacial conditions need both orbital and greenhouse gas forcings to work in unison. For some parameter values, the coupled model also exhibits a critical North American ice sheet configuration, beyond which a stationary-wave–ice-sheet topography feedback can trigger an unabated and unrealistic ice sheet growth. The strong parameter sensitivity found in this study originates from the fact that delicate mass imbalances, as well as errors, are integrated during a transient simulation for thousands of years. This poses a general challenge for transient coupled climate–ice sheet modeling, with such coupled paleo-simulations providing opportunities to constrain such parameters.

Description: In 2015, we have collected more than 60,000 scavenging amphipod specimens during two expeditions to the Clarion-Clipperton fracture Zone (CCZ), in the Northeast (NE) Pacific and to the DISturbance and re-COLonisation (DisCOL) Experimental Area (DEA), a simulated mining impact disturbance proxy in the Peru basin, Southeast (SE) Pacific. Here, we compare biodiversity patterns of the larger specimens (〉15mm) within and between these two oceanic basins. Nine scavenging amphipod species are shared between these two areas, thus indicating connectivity. We further provide evidence that disturbance proxies seem to negatively affect scavenging amphipod biodiversity, as illustrated by a reduced alpha biodiversity in the DEA (Simpson Index (D)=0.62), when compared to the CCZ (D=0.73) and particularly of the disturbance site in the DEA and the site geographically closest to it. Community compositions of the two basins differs, as evidenced by a Non-Metric Dimensional Scaling (NMDS) analysis of beta biodiversity. The NMDS also shows a further separation of the disturbance site (D1) from its neighbouring, undisturbed reference areas (D2, D3, D4 and D5) in the DEA. A single species, Abyssorchomene gerulicorbis, dominates the DEA with 60% of all individuals.

Description: The speciation of dissolved iron (DFe) in the ocean is widely assumed to consist almost exclusively of Fe(III)-ligand complexes. Yet in most aqueous environments a poorly defined fraction of DFe also exists as Fe(II), the speciation of which is uncertain. Here we deploy flow injection analysis to measure in situ Fe(II) concentrations during a series of mesocosm/microcosm/multistressor experiments in coastal environments in addition to the decay rate of this Fe(II) when moved into the dark. During five mesocosm/microcosm/multistressor experiments in Svalbard and Patagonia, where dissolved (0.2 µm) Fe and Fe(II) were quantified simultaneously, Fe(II) constituted 24 %–65 % of DFe, suggesting that Fe(II) was a large fraction of the DFe pool. When this Fe(II) was allowed to decay in the dark, the vast majority of measured oxidation rate constants were less than calculated constants derived from ambient temperature, salinity, pH, and dissolved O2. The oxidation rates of Fe(II) spikes added to Atlantic seawater more closely matched calculated rate constants. The difference between observed and theoretical decay rates in Svalbard and Patagonia was most pronounced at Fe(II) concentrations 〈2 nM, suggesting that the effect may have arisen from organic Fe(II) ligands. This apparent enhancement of Fe(II) stability under post-bloom conditions and the existence of such a high fraction of DFe as Fe(II) challenge the assumption that DFe speciation in coastal seawater is dominated by ligand bound-Fe(III) species.

Description: A strong oxygen-deficient layer is located in the upper layers of the tropical Pacific Ocean and deeper in the North Pacific. Processes related to climate change (upper-ocean warming, reduced ventilation) are expected to change ocean oxygen and nutrient inventories. In most ocean basins, a decrease in oxygen (“deoxygenation”) and an increase in nutrients have been observed in subsurface layers. Deoxygenation trends are not linear and there could be multiple influences on oxygen and nutrient trends and variability. Here oxygen and nutrient time series since 1950 in the Pacific Ocean were investigated at 50 to 300 m depth, as this layer provides critical pelagic habitat for biological communities. In addition to trends related to ocean warming the oxygen and nutrient trends show a strong influence of the Pacific Decadal Oscillation (PDO) in the tropical and the eastern Pacific, and the North Pacific Gyre Oscillation (NPGO) in particular in the North Pacific. In the Oyashio Region the PDO, the NPGO, the North Pacific Index (NPI) and an 18.6-year nodal tidal cycle overlay the long-term trend. In most eastern Pacific regions oxygen increases and nutrients decrease in the 50 to 300 m layer during the negative PDO phase, with opposite trends during the positive PDO phase. The PDO index encapsulates the major mode of sea surface temperature variability in the Pacific, and oxygen and nutrients trends throughout the basin can be described in the context of the PDO phases. El Niño and La Niña years often influence the oxygen and nutrient distribution during the event in the eastern tropical Pacific but do not have a multi-year influence on the trends.

Description: Deciphering the dynamics of dissolved oxygen in the mid-depth ocean during the last deglaciation is essential to understand the influence of climate change on modern oxygen minimum zones (OMZs). Many paleo-proxy records from the Eastern Pacific Ocean indicate an extension of oxygen depleted conditions during the deglaciation but the degree of deoxygenation has not been quantified to date. The Peruvian OMZ, one of the largest OMZs in the world, is a key area to monitor such changes in near-bottom water oxygenation in relation to changing climatic conditions. Here, we analysed the potential to use the composition of foraminiferal assemblages from the Peruvian OMZ as a quantitative redox-proxy. A multiple regression analysis was applied to a joint dataset of living (rose Bengal stained, fossilizable calcareous species) benthic foraminiferal distributions from the Peruvian continental margin. Bottom-water oxygen concentrations ([O2]BW) during sampling were used as dependant variable. The correlation was significant (R2 = 0.82; p 〈 0.05) indicating that the foraminiferal assemblages are rather governed by oxygen availability than by the deposition of particulate organic matter (R2 = 0.53; p = 0.31). We applied the regression formula to four sediment cores from the northern part of the Peruvian OMZ between 3° S and 8° S and 600 m to 1250 m water depths; thereby recording oxygenation changes at the lower boundary of the Peruvian OMZ. Each core displayed a similar trend of decreasing oxygen levels since the Last Glacial Maximum (LGM). The overall [O2]BW change from the Last Glacial Maximum and the Holocene was constrained to 30 μmol/kg at the lower boundary of the OMZ, whereas at shallower depths [O2]BW was relatively stable along the deglaciation. The deoxygenation trend was time-transgressive. It commenced at the southern core, and gradually spread to deeper waters and to the northernmost core location. This pattern indicates a gradual expansion of the OMZ during the last deglaciation, as a result of increasing surface productivity in the Eastern Equatorial Pacific and decreasing advective oxygen supply to intermediate waters off Peru.

Description: Oxygen minimum zones (OMZs) show distinct biogeochemical processes that relate to microorganisms being able to thrive under low or even absent oxygen. Microbial degradation of organic matter is expected to be reduced in OMZs, although quantitative evidence is low. Here, we present heterotrophic bacterial production (3H leucine incorporation), extracellular enzyme rates (leucine aminopeptidase/β-glucosidase) and bacterial cell abundance for various in situ oxygen concentrations in the water column, including the upper and lower oxycline, of the eastern tropical South Pacific off Peru. Bacterial heterotrophic activity in the suboxic core of the OMZ (at in situ ≤ 5 µmol O2 kg−1) ranged from 0.3 to 281 µmol C m−3 d−1 and was not significantly lower than in waters of 5–60 µmol O2 kg−1. Moreover, bacterial abundance in the OMZ and leucine aminopeptidase activity were significantly higher in suboxic waters compared to waters of 5–60 µmol O2 kg−1, suggesting no impairment of bacterial organic-matter degradation in the core of the OMZ. Nevertheless, high cell-specific bacterial production was observed in samples from oxyclines, and cell-specific extracellular enzyme rates were especially high at the lower oxycline, corroborating earlier findings of highly active and distinct micro-aerobic bacterial communities. To assess the impact of bacterial degradation of dissolved organic matter (DOM) for oxygen loss in the Peruvian OMZ, we compared diapycnal fluxes of oxygen and dissolved organic carbon (DOC) and their microbial uptake within the upper 60 m of the water column. Our data indicate low bacterial growth efficiencies of 1 %–21 % at the upper oxycline, resulting in a high bacterial oxygen demand that can explain up to 33 % of the observed average oxygen loss over depth. Our study therewith shows that microbial degradation of DOM has a considerable share in sustaining the OMZ off Peru.

Description: Small steam-driven volcanic explosions are common at volcanoes worldwide but are rarely documented or monitored; therefore, these events still put residents and tourists at risk every year. Steam-driven explosions also occur frequently (once every 2–5 years on average) at Lascar volcano, Chile, where they are often spontaneous and lack any identifiable precursor activity. Here, for the first time at Lascar, we describe the processes culminating in such a sudden volcanic explosion that occurred on October 30, 2015, which was thoroughly monitored by cameras, a seismic network, and gas (SO2 and CO2) and temperature sensors. Prior to the eruption, we retrospectively identified unrest manifesting as a gradual increase in the number of long-period (LP) seismic events in 2014, indicating an augmented level of activity at the volcano. Additionally, SO2 flux and thermal anomalies were detected before the eruption. Then, our weather station reported a precipitation event, followed by changes in the brightness of the permanent volcanic plume and (10 days later) by the sudden volcanic explosion. The multidisciplinary data exhibited short-term variations associated with the explosion, including (1) an abrupt eruption onset that was seismically identified in the 1–10 Hz frequency band, (2) the detection of a 1.7 km high white-grey eruption column in camera images, and (3) a pronounced spike in sulfur dioxide (SO2) emission rates reaching 55 kg sec−1 during the main pulse of the eruption as measured by a mini-DOAS scanner. Continuous CO2 gas and temperature measurements conducted at a fumarole on the southern rim of the Lascar crater revealed a pronounced change in the trend of the relationship between the carbon dioxide (CO2) mixing ratio and the gas outlet temperature; we believe that this change was associated with the prior precipitation event. An increased thermal anomaly inside the active crater observed through Sentinel-2 images and drone overflights performed after the steam-driven explosion revealed the presence of a fracture ~ 50 metres in diameter truncating the dome and located deep inside the active crater, which coincides well with the location of the thermal anomaly. Altogether, these observations lead us to infer that a lava dome was present and subjected to cooling and inhibited degassing. We conjecture that a precipitation event led to the short-term build-up of pressure inside the shallow dome that eventually triggered a vent-clearing phreatic explosion. This study shows the chronology of events culminating in a steam-driven explosion but also demonstrates that phreatic explosions are difficult to forecast, even if the volcano is thoroughly monitored; these findings also emphasize why ascending to the summits of Lascar and similar volcanoes is hazardous, particularly after considerable rainfall.

Description: The intraseasonal evolution of physical and biogeochemical properties during a coastal trapped wave event off central Peru is analysed using data from an extensive shipboard observational programme conducted between April and June 2017, and remote sensing data. The poleward velocities in the Peru–Chile Undercurrent were highly variable and strongly intensified to above 0.5 m s−1 between the middle and end of May. This intensification was likely caused by a first-baroclinic-mode downwelling coastal trapped wave, excited by a westerly wind anomaly at the Equator and originating at about 95∘ W. Local winds along the South American coast did not impact the wave. Although there is general agreement between the observed cross-shore-depth velocity structure of the coastal trapped wave and the velocity structure of first vertical mode solution of a linear wave model, there are differences in the details of the two flow distributions. The enhanced poleward flow increased water mass advection from the equatorial current system to the study site. The resulting shorter alongshore transit times between the Equator and the coast off central Peru led to a strong increase in nitrate concentrations, less anoxic water, likely less fixed nitrogen loss to N2 and a decrease of the nitrogen deficit compared to the situation before the poleward flow intensification. This study highlights the role of changes in the alongshore advection due to coastal trapped waves for the nutrient budget and the cumulative strength of N cycling in the Peruvian oxygen minimum zone. Enhanced availability of nitrate may impact a range of pelagic and benthic elemental cycles, as it represents a major electron acceptor for organic carbon degradation during denitrification and is involved in sulfide oxidation in sediments.

Description: Carbonyl sulfide (OCS) and carbon disulfide (CS2) are volatile sulfur gases that are naturally formed in seawater and exchanged with the atmosphere. OCS is the most abundant sulfur gas in the atmosphere, and CS2 is its most important precursor. They have gained interest due to their direct (OCS) or indirect (CS2 via oxidation to OCS) contribution to the stratospheric sulfate aerosol layer. Furthermore, OCS serves as a proxy to constrain terrestrial CO2 uptake by vegetation. Oceanic emissions of both gases contribute a major part to their atmospheric concentration. Here we present a database of previously published and unpublished, mainly ship-borne measurements in seawater and the marine boundary layer for both gases, available at https://doi.pangaea.de/10.1594/PANGAEA.905430 (Lennartz et al., 2019). The database contains original measurements as well as data digitalized from figures in publications from 42 measurement campaigns, i.e. cruises or time series stations, ranging from 1982 to 2019. OCS data cover all ocean basins except for the Arctic Ocean, as well as all months of the year, while the CS2 dataset shows large gaps in spatial and temporal coverage. Concentrations are consistent across different sampling and analysis techniques for OCS. The database is intended to support the identification of global spatial and temporal patterns and to facilitate the evaluation of model simulations.

Description: Prochlorococcus and Synechococcus are the most abundant free-living photosynthetic microorganisms in the ocean. Uncultivated lineages of these picocyanobacteria also thrive in the dimly illuminated upper part of oxygen-deficient zones (ODZs), where an important portion of ocean nitrogen (N) loss takes place via denitrification and anaerobic ammonium oxidation. Recent metagenomic studies revealed that ODZ Prochlorococcus have the genetic potential for using different N forms, including nitrate and nitrite, uncommon N sources for Prochlorococcus, but common for Synechococcus. To determine which N sources ODZ picocyanobacteria are actually using in nature, the cellular N-15 natural abundance (delta N-15) and assimilation rates of different N compounds were determined using cell sorting by flow cytometry and mass spectrometry. The natural delta N-15 of the ODZ Prochlorococcus varied from -4.0 parts per thousand to 13.0 parts per thousand (n = 9), with 50% of the values in the range of -2.1-2.6 parts per thousand. While the highest values suggest nitrate use, most observations indicate the use of nitrite, ammonium, or a mixture of N sources. Meanwhile, incubation experiments revealed potential assimilation rates of ammonium and urea in the same order of magnitude as that expected for total N in several environments including ODZs, whereas rates of nitrite and nitrate assimilation were very low. Our results thus indicate that reduced forms of N and nitrite are the dominant sources for ODZ picocyanobacteria, although nitrate might be important on some occasions. ODZ picocyanobacteria might thus represent potential competitors with anammox bacteria for ammonium and nitrite, with ammonia-oxidizing archaea for ammonium, and with nitrite-oxidizing bacteria for nitrite.

Description: The plea for using more “realistic,” community‐level, investigations to assess the ecological impacts of global change has recently intensified. Such experiments are typically more complex, longer, more expensive, and harder to interpret than simple organism‐level benchtop experiments. Are they worth the extra effort? Using outdoor mesocosms, we investigated the effects of ocean warming (OW) and acidification (OA), their combination (OAW), and their natural fluctuations on coastal communities of the western Baltic Sea during all four seasons. These communities are dominated by the perennial and canopy‐forming macrophyte Fucus vesiculosus—an important ecosystem engineer Baltic‐wide. We, additionally, assessed the direct response of organisms to temperature and pH in benchtop experiments, and examined how well organism‐level responses can predict community‐level responses to the dominant driver, OW. OW affected the mesocosm communities substantially stronger than acidification. OW provoked structural and functional shifts in the community that differed in strength and direction among seasons. The organism‐level response to OW matched well the community‐level response of a given species only under warm and cold thermal stress, that is, in summer and winter. In other seasons, shifts in biotic interactions masked the direct OW effects. The combination of direct OW effects and OW‐driven shifts of biotic interactions is likely to jeopardize the future of the habitat‐forming macroalga F. vesiculosus in the Baltic Sea. Furthermore, we conclude that seasonal mesocosm experiments are essential for our understanding of global change impact because they take into account the important fluctuations of abiotic and biotic pressures.

Description: Tephra layers produced by volcanic eruptions are widely used for correlation and dating of various deposits and landforms, for synchronization of disparate paleoenvironmental archives, and for reconstruction of magma origin. Here we present our original database TephraKam, which includes chemical compositions of volcanic glass in tephra and welded tuffs from the Kamchatka volcanic arc. The database contains 7049 major element analyses obtained by electron microprobe and 738 trace element analyses obtained by laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) on 487 samples collected in proximity of their volcanic sources in all volcanic zones in Kamchatka. The samples characterize about 300 explosive eruptions, which occurred in Kamchatka from the Pliocene until historic times. Precise or estimated ages for all samples are based on published 39Ar/40Ar dates of rocks and 14C dates of host sediments, statistical age modelling and geologic relationships with dated units. All data in TephraKam is supported by information about source volcanoes and analytical details. Using the data, we present an overview of geochemical variations of Kamchatka volcanic glasses and discuss application of this data for precise identification of tephra layers, their source volcanoes, temporal and spatial geochemical variations of pyroclastic rocks in Kamchatka.

Description: High-resolution optical and hydroacoustic seafloor data acquired in 2015 enabled the reconstruction of disturbance tracks of a past Benthic Impact Experiment that was conducted in 1989 in the Peru Basin in the course of former German environmental impact studies associated with manganese nodule mining. Based on this information, the disturbance level of the experiment regarding the plough impact and distribution and re-deposition of sediment from the evolving sediment plume was assessed qualitatively. Through this, the evolution over the 26 years of a number of the total 78 disturbance tracks could be analyzed which highlights the considerable difference between natural sedimentation in the deep-sea and sedimentation of a resettled sediment plume. Such plumes are seen as one of the most concerning impact associated with potential Mn-nodule mining. Problems in data processing became eminent while dealing with old data from the late 80s, at a time when GPS was just invented and underwater navigation was in an infant stage. However, even today the uncertainties of underwater navigation and the use of a variety of acoustical and optical sensors at different resolutions require detailed post-processing in terms of absolute geographic positioning to improve the overall accuracy of the data. In this study, a ship-based bathymetric map of the survey area was used as absolute geographic reference and a workflow was applied successfully resulting in the most accurate geo-referenced dataset of the DISCOL Experimental Area to date. The new field data were acquired with sensors attached to GEOMARs AUV Abyss and the 0.5 × 1° EM122 multibeam system of RV SONNE during cruise SO242 -1 while the old data first needed to be found and compiled before they could be digitized and properly georeferenced for the presented joined analyses.

Description: The duration and magnitude of the North Atlantic spring bloom impacts both higher trophic levels and oceanic carbon sequestration. Nutrient exhaustion offers a general explanation for bloom termination, but detail on which nutrients and their relative influence on phytoplankton productivity, community structure, and physiology is lacking. Here, we address this using nutrient addition bioassay experiments conducted across the midlatitude North Atlantic in June 2017 (late spring). In four out of six experiments, phytoplankton accumulated over 48–72 h following individual additions of either iron (Fe) or nitrogen (N). In the remaining two experiments, Fe and N were serially limiting, that is, their combined addition sequentially enhanced phytoplankton accumulation. Silicic acid (Si) added in combination with N + Fe led to further chlorophyll a (Chl a) enhancement at two sites. Conversely, addition of zinc, manganese, cobalt, vitamin B12, or phosphate in combination with N + Fe did not. At two sites, the simultaneous supply of all six nutrients, in combination with N + Fe, also led to no further Chl a enhancement, but did result in an additional 30–60% particulate carbon accumulation. This particulate carbon accumulation was not matched by a Redfield equivalent of particulate N, characteristic of high C:N organic exudates that enhance cell aggregation and sinking. Our results suggest that growth rates of larger phytoplankton were primarily limited by Fe and/or N, making the availability of these nutrients the main bottom‐up factors contributing to spring bloom termination. In addition, the simultaneous availability of other nutrients could modify bloom characteristics and carbon export efficiency.

Description: Oxygen-deficient zones (ODZs) are major sites of net natural nitrous oxide (N2O) production and emissions. In order to understand changes in the magnitude of N2O production in response to global change, knowledge on the individual contributions of the major microbial pathways (nitrification and denitrification) to N2O production and their regulation is needed. In the ODZ in the coastal area off Peru, the sensitivity of N2O production to oxygen and organic matter was investigated using 15N tracer experiments in combination with quantitative PCR (qPCR) and microarray analysis of total and active functional genes targeting archaeal amoA and nirS as marker genes for nitrification and denitrification, respectively. Denitrification was responsible for the highest N2O production with a mean of 8.7 nmol L−1 d−1 but up to 118±27.8 nmol L−1 d−1 just below the oxic–anoxic interface. The highest N2O production from ammonium oxidation (AO) of 0.16±0.003 nmol L−1 d−1 occurred in the upper oxycline at O2 concentrations of 10–30 µmol L−1 which coincided with the highest archaeal amoA transcripts/genes. Hybrid N2O formation (i.e., N2O with one N atom from NH+4 and the other from other substrates such as NO−2) was the dominant species, comprising 70 %–85 % of total produced N2O from NH+4, regardless of the ammonium oxidation rate or O2 concentrations. Oxygen responses of N2O production varied with substrate, but production and yields were generally highest below 10 µmol L−1 O2. Particulate organic matter additions increased N2O production by denitrification up to 5-fold, suggesting increased N2O production during times of high particulate organic matter export. High N2O yields of 2.1 % from AO were measured, but the overall contribution by AO to N2O production was still an order of magnitude lower than that of denitrification. Hence, these findings show that denitrification is the most important N2O production process in low-oxygen conditions fueled by organic carbon supply, which implies a positive feedback of the total oceanic N2O sources in response to increasing oceanic deoxygenation.

Description: Seamounts represent ideal systems to study the influence and interdependency of environmental gradients at a single geographic location. These topographic features represent a prominent habitat for various forms of life, including microbiota and macrobiota, spanning benthic as well as pelagic organisms. While it is known that seamounts are globally abundant structures, it still remains unclear how and to which extent the complexity of the sea floor is intertwined with the local oceanographic mosaic, biogeochemistry, and microbiology of a seamount ecosystem. Along these lines, the present study aimed to explore whether and to what extent seamounts can have an imprint on the microbial community composition of seawater and of sessile benthic invertebrates, sponges. For our high-resolution sampling approach of microbial diversity (16S rRNA gene amplicon sequencing) along with measurements of inorganic nutrients and other biogeochemical parameters, we focused on the Schulz Bank seamount ecosystem, a sponge ground ecosystem which is located on the Arctic Mid-Ocean Ridge. Seawater samples were collected at two sampling depths (mid-water, MW, and near-bed water, BW) from a total of 19 sampling sites. With a clustering approach we defined microbial microhabitats within the pelagic realm at Schulz Bank, which were mapped onto the seamount's topography and related to various environmental parameters (such as suspended particulate matter, SPM; dissolved inorganic carbon, DIC; silicate, SiO−4; phosphate, PO3−4; ammonia, NH+4; nitrate, NO2−3; nitrite, NO−2; depth; and dissolved oxygen, O2). The results of our study reveal a “seamount effect” (sensu stricto) on the microbial mid-water pelagic community at least 200 m above the sea floor. Further, we observed a strong spatial heterogeneity in the pelagic microbial landscape across the seamount, with planktonic microbial communities reflecting oscillatory and circulatory water movements, as well as processes of bentho-pelagic coupling. Depth, NO2−3, SiO−4, and O2 concentrations differed significantly between the determined pelagic microbial clusters close to the sea floor (BW), suggesting that these parameters were presumably linked to changes in microbial community structures. Secondly, we assessed the associated microbial community compositions of three sponge species along a depth gradient of the seamount. While sponge-associated microbial communities were found to be mainly species-specific, we also detected significant intra-specific differences between individuals, depending on the pelagic near-bed cluster they originated from. The variable microbial phyla (i.e. phyla which showed significant differences across varying depth, NO2−3, SiO−4, O2 concentrations, and different from local seawater communities) were distinct for every sponge species when considering average abundances per species. Variable microbial phyla included representatives of both those taxa traditionally counted for the variable community fraction and taxa counted traditionally for the core community fraction. Microbial co-occurrence patterns for the three examined sponge species Geodia hentscheli, Lissodendoryx complicata, and Schaudinnia rosea were distinct from each other. Over all, this study shows that topographic structures such as the Schulz Bank seamount can have an imprint (seamount effect sensu lato) on both the microbial community composition of seawater and sessile benthic invertebrates such as sponges by an interplay between the geology, physical oceanography, biogeochemistry, and microbiology of seamounts.

Description: The copepod Acartia tonsa is a key component of a wide range of marine ecosystems, linking energy transfer from phytoplankton to higher trophic levels, and has a central role in productivity and biogeochemistry. The interaction of end-of-century global warming and ocean acidification scenarios with testing moderate temperature effects on a seminatural copepod community is needed to understand future community functioning. Here, we deployed a mesocosm experimental set-up with a full factorial design using two temperatures (13°C and 19°C) crossed with a pCO2 gradient ranging from ambient (550 μatm) to 3000 μatm. We used the natural bacteria, phyto- and microzooplankton species composition and biomass of the Kiel Bight and tested the response of A. tonsa development, carbon growth, mortality, size and condition. The tested traits were differently affected by the interaction of temperature and acidification. Ocean acidification increased development, carbon growth, size and mortality under the warming scenario of 19°C. At 13°C mortality rates decreased, while carbon growth, size and condition increased with acidification. We conclude from our experimental approach that a single species shows a variety of responses depending on the focal functional trait. Trait-specific mesozooplankton responses need to be further investigated and compared between geographical regions, seasons and taxonomic groups.

Description: The largest and commercially appealing mineral deposits can be found in the abyssal seafloor of the Clarion-Clipperton Zone (CCZ), a polymetallic nodule province, in the NE Pacific Ocean, where experimental mining is due to take place. In anticipation of deep-sea mining impacts, it has become essential to rapidly and accurately assess biodiversity. For this reason, ophiuroid material collected during seven scientific cruises from five exploration license areas within CCZ, one area protected from mining (APEI3, Area of Particular Environmental Interest) in the periphery of CCZ and the DIS-turbance and re-COLonisation (DISCOL) Experimental Area (DEA), in the SE Pacific Ocean, was examined. Specimens were genetically analysed using a fragment of the mitochondrial cytochrome c oxidase subunit I (COI). Maximum Likelihood and Neighbour Joining trees were constructed, while four tree-based and distance-based methods of species delineation (ABGD, BINs, GMYC, mPTP) were employed to propose Secondary Species Hypotheses (SSHs) within the ophiuroids collected. The species delimitations analyses concordant results revealed the presence of 43 deep-sea brittle stars SSHs, revealing an unexpectedly high diversity and showing that the most conspicuous invertebrates in abyssal plains have been so far considerably under-estimated. The number of SSHs found in each area varied from 5 (IFREMER area) to 24 (BGR area), while 13 SSHs were represented by singletons. None of the SSHs was found to be present in all 7 areas, while the majority of species (44.2 %) had a single-area presence (19 SSHs). The most common species were Ophioleucidae sp. (Species 29), Amphioplus daleus (Species 2) and Ophiosphalma glabrum (Species 3), present in all areas except APEI3. The biodiversity patterns could be mainly attributed to POC fluxes that could explain the highest species numbers found in BGR (German contractor area) and UKSRL (UK contractor area) areas. The five exploration contract areas belong to a mesotrophic province, while in contrary the APEI3 is located in an oligotrophic province which could explain the lowest diversity as well as very low similarity with the other six study areas. Based on these results the representativeness and the appropriateness of APEI3 to meet its purpose of preserving the biodiversity of the CCZ fauna are questioned. Finally, this study provides the foundation for biogeographic and functional analyses that will provide insight into the drivers of species diversity and its role in ecosystem function.

Description: Dissolved Fe (DFe) samples from the GEOVIDE voyage (GEOTRACES GA01, May–June 2014) in the North Atlantic Ocean were analyzed using a seaFAST-pico™ coupled to an Element XR sector field inductively coupled plasma mass spectrometer (SF-ICP-MS) and provided interesting insights into the Fe sources in this area. Overall, DFe concentrations ranged from 0.09±0.01 to 7.8±0.5 nmol L−1. Elevated DFe concentrations were observed above the Iberian, Greenland, and Newfoundland margins likely due to riverine inputs from the Tagus River, meteoric water inputs, and sedimentary inputs. Deep winter convection occurring the previous winter provided iron-to-nitrate ratios sufficient to sustain phytoplankton growth and lead to relatively elevated DFe concentrations within subsurface waters of the Irminger Sea. Increasing DFe concentrations along the flow path of the Labrador Sea Water were attributed to sedimentary inputs from the Newfoundland Margin. Bottom waters from the Irminger Sea displayed high DFe concentrations likely due to the dissolution of Fe-rich particles in the Denmark Strait Overflow Water and the Polar Intermediate Water. Finally, the nepheloid layers located in the different basins and at the Iberian Margin were found to act as either a source or a sink of DFe depending on the nature of particles, with organic particles likely releasing DFe and Mn particle scavenging DFe.

Description: In the abyssal equatorial Pacific Ocean, most of the seafloor of the Clarion-Clipperton Fracture Zone (CCFZ), a 6 million km2 polymetallic nodule province, has been preempted for future mining. In light of the large environmental footprint that mining would leave and given the diversity and the vulnerability of the abyssal fauna, the International Seabed Authority has implemented a regional management plan that includes the creation of nine Areas of Particular Environmental Interest (APEIs) located at the periphery of the CCFZ. The scientific principles for the design of the APEIs were based on the best – albeit very limited – knowledge of the area. The fauna and habitats in the APEIs are unknown, as are species' ranges and the extent of biodiversity across the CCFZ. As part of the Joint Programming Initiative Healthy and Productive Seas and Oceans (JPI Oceans) pilot action “Ecological aspects of deep-sea mining”, the SO239 cruise provided data to improve species inventories, determine species ranges, identify the drivers of beta diversity patterns and assess the representativeness of an APEI. Four exploration contract areas and an APEI (APEI no. 3) were sampled along a gradient of sea surface primary productivity that spanned a distance of 1440 km in the eastern CCFZ. Between three and eight quantitative box cores (0.25 m2; 0–10 cm) were sampled in each study area, resulting in a large collection of polychaetes that were morphologically and molecularly (cytochrome c oxidase subunit I and 16S genes) analyzed. A total of 275 polychaete morphospecies were identified. Only one morphospecies was shared among all five study areas and 49 % were singletons. The patterns in community structure and composition were mainly attributed to variations in organic carbon fluxes to the seafloor at the regional scale and nodule density at the local scale, thus supporting the main assumptions underlying the design of the APEIs. However, the APEI no. 3, which is located in an oligotrophic province and separated from the CCFZ by the Clarion Fracture Zone, showed the lowest densities, lowest diversity, and a very low and distant independent similarity in community composition compared to the contract areas, thus questioning the representativeness and the appropriateness of APEI no. 3 to meet its purpose of diversity preservation. Among the four exploration contracts, which belong to a mesotrophic province, the distance decay of similarity provided a species turnover of 0.04 species km−1, an average species range of 25 km and an extrapolated richness of up to 240 000 polychaete species in the CCFZ. By contrast, nonparametric estimators of diversity predict a regional richness of up to 498 species. Both estimates are biased by the high frequency of singletons in the dataset, which likely result from under-sampling and merely reflect our level of uncertainty. The assessment of potential risks and scales of biodiversity loss due to nodule mining thus requires an appropriate inventory of species richness in the CCFZ.

Description: Climate change in Siberia is currently receiving a lot of attention because large permafrost-covered areas could provide a strong positive feedback to global warming through the release of carbon that has been sequestered there on glacial–interglacial timescales. Geological evidence and climate model experiments show that the Siberian region also played an exceptional role during glacial periods. The region that is currently known for its harsh cold climate did not experience major glaciations during the last ice age, including its severest stages around the Last Glacial Maximum (LGM). On the contrary, it is thought that glacial summer temperatures were comparable to the present day. However, evidence of glaciation has been found for several older glacial periods. We combine LGM experiments from the second and third phases of the Paleoclimate Modelling Intercomparison Project (PMIP2 and PMIP3) with sensitivity experiments using the Community Earth System Model (CESM). Together, these climate model experiments reveal that the intermodel spread in LGM summer temperatures in Siberia is much larger than in any other region of the globe and suggest that temperatures in Siberia are highly susceptible to changes in the imposed glacial boundary conditions, the included feedbacks and processes, and to the model physics of the different components of the climate model. We find that changes in the circumpolar atmospheric stationary wave pattern and associated northward heat transport drive strong local snow and vegetation feedbacks and that this combination explains the susceptibility of LGM summer temperatures in Siberia. This suggests that a small difference between two glacial periods in terms of climate, ice buildup or their respective evolution towards maximum glacial conditions can lead to strongly divergent summer temperatures in Siberia, allowing for the buildup of an ice sheet during some glacial periods, while during others, above-freezing summer temperatures preclude a multi-year snowpack from forming.

Description: In this paper we introduce a Bayesian framework, which is explicit about prior assumptions, for using model ensembles and observations together to constrain future climate change. The emergent constraint approach has seen broad application in recent years, including studies constraining the equilibrium climate sensitivity (ECS) using the Last Glacial Maximum (LGM) and the mid-Pliocene Warm Period (mPWP). Most of these studies were based on ordinary least squares (OLS) fits between a variable of the climate state, such as tropical temperature, and climate sensitivity. Using our Bayesian method, and considering the LGM and mPWP separately, we obtain values of ECS of 2.7 K (0.6–5.2, 5th–95th percentiles) using the PMIP2, PMIP3, and PMIP4 datasets for the LGM and 2.3 K (0.5–4.4) with the PlioMIP1 and PlioMIP2 datasets for the mPWP. Restricting the ensembles to include only the most recent version of each model, we obtain 2.7 K (0.7–5.2) using the LGM and 2.3 K (0.4–4.5) using the mPWP. An advantage of the Bayesian framework is that it is possible to combine the two periods assuming they are independent, whereby we obtain a tighter constraint of 2.5 K (0.8–4.0) using the restricted ensemble. We have explored the sensitivity to our assumptions in the method, including considering structural uncertainty, and in the choice of models, and this leads to 95 % probability of climate sensitivity mostly below 5 K and only exceeding 6 K in a single and most uncertain case assuming a large structural uncertainty. The approach is compared with other approaches based on OLS, a Kalman filter method, and an alternative Bayesian method. An interesting implication of this work is that OLS-based emergent constraints on ECS generate tighter uncertainty estimates, in particular at the lower end, an artefact due to a flatter regression line in the case of lack of correlation. Although some fundamental challenges related to the use of emergent constraints remain, this paper provides a step towards a better foundation for their potential use in future probabilistic estimations of climate sensitivity.