550 Geowissenschaften
Refine
Year of publication
Document Type
- Article (907)
- Doctoral Thesis (194)
- Contribution to a Periodical (32)
- Book (26)
- Working Paper (22)
- Part of Periodical (21)
- Conference Proceeding (18)
- Part of a Book (9)
- Diploma Thesis (8)
- diplomthesis (7)
- Other (5)
- Report (5)
- Master's Thesis (4)
- Review (4)
- Periodical (3)
- Bachelor Thesis (2)
- Habilitation (2)
- Preprint (2)
Language
Keywords
- climate change (11)
- Climate change (8)
- Klima (8)
- Klimaänderung (7)
- Modellierung (7)
- COSMO-CLM (6)
- Klimawandel (6)
- Palaeoclimate (6)
- precipitation (6)
- Atmospheric chemistry (5)
Institute
- Geowissenschaften (778)
- Geowissenschaften / Geographie (148)
- Biodiversität und Klima Forschungszentrum (BiK-F) (62)
- Geographie (61)
- Senckenbergische Naturforschende Gesellschaft (57)
- Präsidium (44)
- Extern (28)
- Biowissenschaften (22)
- Institut für Ökologie, Evolution und Diversität (10)
- Physik (8)
Accretion rates of Holocene tropical coral reefs in three areas in the Atlantic, Pacific and Indian Oceans have been quantified in 79 dated core sections in 34 reef cores from Belize, the Maldives and French Polynesia. Holocene vertical reef accretion rate averages 5.05 m/kyr and has decreased during the past 10 kyr. Accretion rates in branched and massive coral facies are statistically similar. Reef accretion rate is positively correlated with the rate of sea‐level rise, that is the degree of creation of accommodation space, and with climate as expressed in a Holocene sea surface temperature anomaly. Accommodation space is also created by subsidence, but at a rate one to two orders of magnitude lower than that created by glacio‐eustasy (0.04 to 0.16 m/kyr). Lagoonal background sedimentation in adjacent reef lagoons averages 0.89 m/kyr as measured in 72 dated core sections in 28 cores. Lagoonal carbonate sedimentation on top of underlying mangrove peat usually starts after a considerable hiatus of ca 3 kyr on average. The lagoonal background sedimentation rate increased during the Holocene, probably due to deepening. The differences between vertical reef accretion and lagoonal background sedimentation rates are a major factor in the production of the widely known saucer shapes typical of tropical reefs and carbonate platforms, that is the creation of unfilled accommodation space. Reef core recovery, used as a proxy for reef consolidation, and core depth exhibit a statistically negative correlation based on data from 326 core barrels. Recovery and marine cement abundance (average volume 8.6%) also decrease from windward to leeward core positions. These observations are presumably a result of both a decrease in the rate of sea‐level rise that is the increase in time available for submarine cementation during the Holocene and the amount of flushing of reef interstices by marine waters.
The frequency of extreme events has changed, having a direct impact on human lives. Regional climate models help us to predict these regional climate changes. This work presents an atmosphere–ocean coupled regional climate system model (RCSM; with the atmospheric component COSMO-CLM and the ocean component NEMO) over the European domain, including three marginal seas: the Mediterranean, North, and Baltic Sea. To test the model, we evaluate a simulation of more than 100 years (1900–2009) with a spatial grid resolution of about 25 km. The simulation was nested into a coupled global simulation with the model MPI-ESM in a low-resolution configuration, whose ocean temperature and salinity were nudged to the ocean–ice component of the MPI-ESM forced with the NOAA 20th Century Reanalysis (20CR). The evaluation shows the robustness of the RCSM and discusses the added value by the coupled marginal seas over an atmosphere-only simulation. The coupled system is stable for the complete 20th century and provides a better representation of extreme temperatures compared to the atmosphere-only model. The produced long-term dataset will help us to better understand the processes leading to meteorological and climate extremes.
The sources and critical enrichment processes for granite related tin ores are still not well understood. The Erzgebirge represents one of the classical regions for tin mineralization. We investigated the four largest plutons from the Western Erzgebirge (Germany) for the geochemistry of bulk rocks and autocrystic zircons and relate this information to their intrusion ages. The source rocks of the Variscan granites were identified as high-grade metamorphic rocks based on the comparison of Hf-O isotope data on zircons, the abundance of xenocrystic zircon ages as well as Nd and Hf model ages. Among these rocks, restite is the most likely candidate for later Variscan melts. Based on the evolution with time, we could reconstruct enrichment factors for tin and tungsten starting from the protoliths (575 Ma) that were later converted to high-grade metamorphic rocks (340 Ma) and served as sources for the older biotite granites (323–318 Ma) and the tin granites (315–314 Ma). This evolution involved a continuous enrichment of both tin and tungsten with an enrichment factor of ~15 for tin and ~7 for tungsten compared to the upper continental crust (UCC). Ore level concentrations (>10–100 times enrichment) were achieved only in the greisen bodies and dykes by subsequent hydrothermal processes.
Unter dem Namen #climonomics – EU-Klimakonferenz für Schüler*innen: A Friday for Future veranstaltete das PolECulE-Projekt, ein Kooperationsprojekt der Englisch- und Politikdidaktik, am 25. Oktober eine Veranstaltung für Schüler*innen. Hier wurde die aktuelle EU-Klimapolitik im Rahmen eines Rollenspiels erörtert und debattiert. Nachdem die Hessenschau, FAZ, der Deutschlandfunk und die Frankfurter Rundschau ausgiebig über das Projekt berichteten, schreibt auch Charlotte Wittich (17), angehende Abiturientin am Heinrich-von-Gagern-Gymnasium, von ihren Erfahrungen.
Parabens and sorbic acid are commonly used as food preservatives due to their antimicrobial effect. However, their use in foods for infants and young children is not permitted in the European Union. Previous studies found these compounds in some gel-filled baby teethers, whereby parabens, which are well-known as endocrine disruptors, were identified in the polymer-based chewing surface consisting of ethylene-vinyl acetate (EVA). To assess the exposure of infants and young children to these products, the application of parabens in teethers should be thoroughly investigated. Therefore, the present study aimed to apply a representative migration test procedure combined with an accurate analytical method to examine gel-filled baby teethers without elaborate sample preparation, high costs, and long processing times. Accordingly, solid-phase extraction (SPE), in combination with a stable isotope dilution assay (SIDA) and subsequent gas chromatography–mass spectrometry (GC–MS) for analysis of methyl-, ethyl-, and n-propylparaben (MeP, EtP, and n-PrP), was found to be well-suited, with recoveries ranging from 93 to 99%. The study compared the release of these parabens from intact teether surfaces into water and saliva simulant under real-life conditions, with total amounts of detected parabens found to be in the range of 101–162 µg 100 mL−1 and 57–148 µg 100 mL−1, respectively. Furthermore, as a worst-case scenario, the release into water was examined using a long-term migration study.
A twentieth century-long coupled atmosphere-ocean regional climate simulation with COSMO-CLM (Consortium for Small-Scale Modeling, Climate Limited-area Model) and NEMO (Nucleus for European Modelling of the Ocean) is studied here to evaluate the added value of coupled marginal seas over continental regions. The interactive coupling of the marginal seas, namely the Mediterranean, the North and the Baltic Seas, to the atmosphere in the European region gives a comprehensive modelling system. It is expected to be able to describe the climatological features of this geographically complex area even more precisely than an atmosphere-only climate model. The investigated variables are precipitation and 2 m temperature. Sensitivity studies are used to assess the impact of SST (sea surface temperature) changes over land areas. The different SST values affect the continental precipitation more than the 2 m temperature. The simulated variables are compared to the CRU (Climatic Research Unit) observational data, and also to the HOAPS/GPCC (Hamburg Ocean Atmosphere Parameters and Fluxes from Satellite Data, Global Precipitation Climatology Centre) data. In the coupled simulation, added skill is found primarily during winter over the eastern part of Europe. Our analysis shows that, over this region, the coupled system is dryer than the uncoupled system, both in terms of precipitation and soil moisture, which means a decrease in the bias of the system. Thus, the coupling improves the simulation of precipitation over the eastern part of Europe, due to cooler SST values and in consequence, drier soil.
Often in climate system studies, linear and symmetric statistical measures are applied to quantify interactions among subsystems or variables. However, they do not allow identification of the driving and responding subsystems. Therefore, in this study, we aimed to apply asymmetric measures from information theory: the axiomatically proposed transfer entropy and the first principle-based information flow to detect and quantify climate interactions. As their estimations are challenging, we initially tested nonparametric estimators like transfer entropy (TE)-binning, TE-kernel, and TE k-nearest neighbor and parametric estimators like TE-linear and information flow (IF)-linear with idealized two-dimensional test cases along with their sensitivity on sample size. Thereafter, we experimentally applied these methods to the Lorenz-96 model and to two real climate phenomena, i.e., (1) the Indo-Pacific Ocean coupling and (2) North Atlantic Oscillation (NAO)–European air temperature coupling. As expected, the linear estimators work for linear systems but fail for strongly nonlinear systems. The TE-kernel and TE k-nearest neighbor estimators are reliable for linear and nonlinear systems. Nevertheless, the nonparametric methods are sensitive to parameter selection and sample size. Thus, this work proposes a composite use of the TE-kernel and TE k-nearest neighbor estimators along with parameter testing for consistent results. The revealed information exchange in Lorenz-96 is dominated by the slow subsystem component. For real climate phenomena, expected bidirectional information exchange between the Indian and Pacific SSTs was detected. Furthermore, expected information exchange from NAO to European air temperature was detected, but also unexpected reversal information exchange. The latter might hint to a hidden process driving both the NAO and European temperatures. Hence, the limitations, availability of time series length and the system at hand must be taken into account before drawing any conclusions from TE and IF-linear estimations.
The goal of limited area models (LAMs) is to downscale coarse-gridded general circulation model output to represent small-scale features of weather and climate. The LAM needs information from the driving coarse-gridded model passing through its lateral boundaries. The treatment of this information transfer causes inconsistencies between driving and nested models and, subsequently, issues in regional weather and climate simulations. This work examines errors arising from choices taken by the modeler (temporal update frequency of boundary data, spatial resolution jump, and numerical lateral boundary formulation) systematically in an idealized simulation environment. So-called Big-Brother Experiments were performed with the LAM COSMO-CLM (0.11° grid spacing). A baroclinic wave in a zonal channel was simulated over flat terrain with and without a Gaussian hill. The results reveal that the quality of the driving data, here represented by simulations only differing from the LAM simulations by reduced spatial resolution, dominates the performance of the nested model. Consequently, at the simulated mesoscale, the performance of the nested small-scale model simulations is weakly sensitive to the numerical lateral boundary formulation (Davies relaxation or the newly implemented, computationally less demanding Mesinger Eta-model formulation). The performance sensitivity to boundary update frequency and resolution jump is small when at least 6-hourly updates and a resolution jump factor of maximally six is used. Gaussian hill LAM simulations illustrated the strength of downscaling; they can represent small-scale features missing in the coarse-scale driving simulations. In the idealized simulation experiments, spectral nudging is not advisable as it imprints the driving models deficits on the nested simulation.
In der vorliegenden Arbeit wurde ein neuer optischer Aufbau für das Laserlabor der Abteilung Kristallographie im FB 11 an der Goethe-Universität Frankfurt beschrieben. Mit Hilfe dieses Aufbaus konnten verschiedene spektroskopische Methoden genutzt werden, um die - von Druck und Temperatur abhängige - Phasenstabilität von Calcium- und Eisencarbonaten zu untersuchen. Mit Hilfe von Raman-Spektroskopie konnte das Phasendiagramm von Calciumcarbonat (CaCO3) teilweise neu bestimmt werden. Fluoreszenzuntersuchungen an dotierten CaCO3 Proben ergaben, dass sich Europium-dotierter Calcit zunächst in eine amorphe Form umwandelt, bevor er bei ca. 15 GPa in eine amorphe 'aragonitische' Form umgewandelt wird. Die Umwandlung ist nicht reversibel. Laserheizexperimente bei 18.5 GPa an dotiertem Siderit (FeCO3) führten zur Bildung eines neuen Hochdruck-Hochtemperatur FeCO3 -Polymorphs. Die Strukturlösung erfolgte mit Hilfe von Röntgendaten, die am Deutschen Elektronen-Synchrotron (DESY) in Hamburg gewonnen wurden. Schließlich wurde eine neue Methode zur Bestimmung von Temperaturen in Laserheizexperimenten beschrieben. Sie beruht auf der Abschwächung eines Fluoreszenzsignals durch die Temperatur, welche durch die Wechselwirkung eines Heizlasers mit der Probe erzeugt wird.
The synthesis of polynitrogen compounds is of fundamental importance due to their potential as environmentally-friendly high energy density materials. Attesting to the intrinsic difficulties related to their formation, only three polynitrogen ions, bulk stabilized as salts, are known. Here, magnesium and molecular nitrogen are compressed to about 50 GPa and laser-heated, producing two chemically simple salts of polynitrogen anions, MgN4 and Mg2N4. Single-crystal X-ray diffraction reveals infinite anionic polythiazyl-like 1D N-N chains in the crystal structure of MgN4 and cis-tetranitrogen N44− units in the two isosymmetric polymorphs of Mg2N4. The cis-tetranitrogen units are found to be recoverable at atmospheric pressure. Our results respond to the quest for polynitrogen entities stable at ambient conditions, reveal the potential of employing high pressures in their synthesis and enrich the nitrogen chemistry through the discovery of other nitrogen species, which provides further possibilities to design improved polynitrogen arrangements.
In dieser Dissertation wird die Parametrisierung von subgitterskaligen (SGS) Prozessen in Atmosphärenmodellen untersucht. Die Arbeit befasst sich mit den stochastisch angetriebenen Flachwassergleichungen, im ersten Teil in einer räumlichen Dimension und im zweiten Teil in zwei Dimensionen. Die Einteilung in aufgelöste und SGS-Variable erfolgt in beiden Fällen über lokale räumliche Mittel der Ursprungsvariable und deren Abweichungen vom lokalen Mittel.
Im eindimensionalen Fall liegt zwischen den Variablen eine deutliche Separation der charakteristischen Zeitskalen vor, wodurch die Anwendung der stochastischen Moden Reduktion (SMR) ermöglicht wird. Die SMR generiert ein reduziertes Modell der aufgelösten Variable mit einer stochastischen SGS-Parametrisierung, im Folgenden auch Schließung genannt. Die SMR-Schließung basiert auf den Grundgleichungen des Flachwassermodells und ist numerisch effizient einsetzbar, da sie nur eine geringe Anzahl von benachbarten Zellen koppelt. Sie verbessert die Ergebnisse des reduzierten Modells und übertrifft die Ergebnisse zweier zum Vergleich untersuchter empirischer stochastischer Schließungen. Den größten Zugewinn liefert sie im Energiespektrum, insbesondere für kleine Skalen. Das Ergebnis der SMR-Schließung kann verbessert werden, indem die Amplitude der stochastischen Schließungskomponente gedämpft wird. Die SMR-Schließung ist skalenabhängig im Sinne der räumlichen Modellauflösung. Untersucht wird die Schließung bei Halbierung und Viertelung der räumlichen Auflösung, wo sie ihre Überlegenheit gegenüber den empirischen Schließungen wiederholt bestätigt.
Im Unterschied zum eindimensionalen Fall ist in zwei Dimensionen auch die Corioliskraft enthalten und eine räumliche Divergenz der Schwerewellen möglich. Zwischen der aufgelösten und der SGS-Variable kommt es erneut zu einer Separation der charakteristischen Zeitskalen. Die Separation ist allerdings weniger stark ausgeprägt als im eindimensionalen Fall. Grund hierfür ist das Auftreten einer lang korrelierten geostrophisch balancierten Mode, welche auch auf die SGS-Variable projiziert. Das Vorgehen zur Bestimmung der SMR-Schließung für das zweidimensionale Modell verläuft analog zum eindimensionalen Fall. Es werden die Ergebnisse des hoch aufgelösten Referenzmodells und zweier Modelle ohne SGS-Schließung verglichen.
The impact of precipitation in shallow cumulus convection on the moisture variance and third-order moments of moisture is investigated with the help of large-eddy simulations. Three idealized simulations based on the Rain in Cumulus over the Ocean field experiment are analyzed: one nonprecipitating, on a smaller domain, and two precipitating cases, on a larger domain with different initial profiles of moisture. Results show that precipitation and the associated cloud organization lead to increased generation of higher-order moments (HOM) of moisture compared to the nonprecipitating case. To understand the physical mechanism and the role of individual processes in this increase, budgets of HOM of moisture are studied. Microphysics directly decreases the generation of HOM of moisture, but this effect is not dominant. The gradient production term is identified as the main source term in the HOM budgets. The influence of the gradient production term on moisture variance is further examined separately in cloud active and nonactive regions. The main contribution to the gradient production term comes from the smaller cloud active region because of the stronger moisture flux. Further analyses of the horizontal and vertical cross sections of moisture fluctuations show that the precipitation-induced downdrafts and updrafts are the main mechanism for the generation of moisture variance. The variance increase is linked to shallow dry downdraft regions with horizontal divergence in the subcloud layer, moist updrafts with horizontal convergence in the bulk cloud layer, and finally wider areas of horizontal divergence in the cloud inversion layer.
Motions on planetary spatial scales in the atmosphere are governed by the planetary geostrophic equations. However, little attention has been paid to the interaction between the baroclinic and barotropic flows within the planetary geostrophic scaling. This is the focus of the present study, which utilizes planetary geostrophic equations for a Boussinesq fluid supplemented by a novel evolution equation for the barotropic flow. The latter is affected by meridional momentum flux due to baroclinic flow and drag by the surface wind. The barotropic wind, on the other hand, affects the baroclinic flow through buoyancy advection. Via a relaxation towards a prescribed buoyancy profile the model produces realistic major features of the zonally symmetric wind and temperature fields. We show that there is considerable cancellation between the barotropic and the baroclinic surface zonal mean zonal winds. Linear and nonlinear model responses to steady diabatic zonally asymmetric forcing are investigated, and the arising stationary waves are interpreted in terms of analytical solutions. We also study the problem of baroclinic instability on the sphere within the present model.
Convective shower characteristics simulated with the convection-permitting climate model COSMO-CLM
(2019)
This paper evaluates convective precipitation as simulated by the convection-permitting climate model (CPM) Consortium for Small-Scale Modeling in climate mode (COSMO-CLM) (with 2.8 km grid-spacing) over Germany in the period 2001–2015. Characteristics of simulated convective precipitation objects like lifetime, area, mean intensity, and total precipitation are compared to characteristics observed by weather radar. For this purpose, a tracking algorithm was applied to simulated and observed precipitation with 5-min temporal resolution. The total amount of convective precipitation is well simulated, with a small overestimation of 2%. However, the simulation underestimates convective activity, represented by the number of convective objects, by 33%. This underestimation is especially pronounced in the lowlands of Northern Germany, whereas the simulation matches observations well in the mountainous areas of Southern Germany. The underestimation of activity is compensated by an overestimation of the simulated lifetime of convective objects. The observed mean intensity, maximum intensity, and area of precipitation objects increase with their lifetime showing the spectrum of convective storms ranging from short-living single-cell storms to long-living organized convection like supercells or squall lines. The CPM is capable of reproducing the lifetime dependence of these characteristics but shows a weaker increase in mean intensity with lifetime resulting in an especially pronounced underestimation (up to 25%) of mean precipitation intensity of long-living, extreme events. This limitation of the CPM is not identifiable by classical evaluation techniques using rain gauges. The simulation can reproduce the general increase of the highest percentiles of cell area, total precipitation, and mean intensity with temperature but fails to reproduce the increase of lifetime. The scaling rates of mean intensity and total precipitation resemble observed rates only in parts of the temperature range. The results suggest that the evaluation of coarse-grained (e.g., hourly) precipitation fields is insufficient for revealing challenges in convection-permitting simulations.
In November 2016, magnetotelluric (MT) data were collected at the Ceboruco Volcano in cooperation with the Centro de Sismología y Volcanología de Occidente (SisVoc, Universidad de Guadalajara, Mexico). The Ceboruco is a 2280 m high stratovolcano, located in Nayarit State, Mexico. It is placed in the central part of the Tepic-Zacoalco Rift (TZR), which constitutes the north-western end of the Trans-Mexican Volcanic Belt. Together with Chapala and Colima (in the Jalisco Block), they form the triple rift system developed as a consequence of the ongoing subduction of the Rivera and Cocos oceanic plates beneath the North American continental crust. Although its last eruption occurred in 1870, it is the most active volcano in the area, showing volcanic-earthquake activity together with ongoing vapor emissions. The survey was part of a geothermal project (CeMIEGeo-P24) and focused on the determination of electrical conductivity properties to characterize the deep structure and the geothermal potential of the Volcano. Frequency dependent magnetotelluric response functions were calculated from 25 broadband MT stations, which covered an area of 10 x 10 km2 including its crater, calderas and foreland. The results were interpreted using anisotropic 3-D forward modelling and isotropic 3-D inversion approaches, considering strong topographical effects. The final resistivity model implies a highly conductive layer, reaching from near-surface to approximately 2 km depth, which might be related to a hydrothermal system. Here, mineralized fluids and clay minerals can cause high conductivities around 1 S/m. For longer periods, the principal axes of the MT response tensors (phase tensor, apparent resistivity tensor) are in good agreement with the strike direction of the underlying rift system. However, they are not rendered by the isotropic inversion. Thus the data suggest an anisotropic electrical conductivity at greater depth with its principal axis determined by the response tensors.
Deriving stratospheric age of air spectra using an idealized set of chemically active trace gases
(2019)
Analysis of stratospheric transport from an observational point of view is frequently realized by evaluation of the mean age of air values from long-lived trace gases. However, this provides more insight into general transport strength and less into its mechanism. Deriving complete transit time distributions (age spectra) is desirable, but their deduction from direct measurements is difficult. It is so far primarily based on model work. This paper introduces a modified version of an inverse method to infer age spectra from mixing ratios of short-lived trace gases and investigates its basic principle in an idealized model simulation. For a full description of transport seasonality the method includes an imposed seasonal cycle to gain multimodal spectra. An ECHAM/MESSy Atmospheric Chemistry (EMAC) model simulation is utilized for a general proof of concept of the method and features an idealized dataset of 40 radioactive trace gases with different chemical lifetimes as well as 40 chemically inert pulsed trace gases to calculate pulse age spectra. It is assessed whether the modified inverse method in combination with the seasonal cycle can provide matching age spectra when chemistry is well-known. Annual and seasonal mean inverse spectra are compared to pulse spectra including first and second moments as well as the ratio between them to assess the performance on these timescales. Results indicate that the modified inverse age spectra match the annual and seasonal pulse age spectra well on global scale beyond 1.5 years of mean age of air. The imposed seasonal cycle emerges as a reliable tool to include transport seasonality in the age spectra. Below 1.5 years of mean age of air, tropospheric influence intensifies and breaks the assumption of single entry through the tropical tropopause, leading to inaccurate spectra, in particular in the Northern Hemisphere. The imposed seasonal cycle wrongly prescribes seasonal entry in this lower region and does not lead to a better agreement between inverse and pulse age spectra without further improvement. Tests with a focus on future application to observational data imply that subsets of trace gases with 5 to 10 species are sufficient for deriving well-matching age spectra. These subsets can also compensate for an average uncertainty of up to ±20 % in the knowledge of chemical lifetime if a deviation of circa ±10 % in modal age and amplitude of the resulting spectra is tolerated.
An accelerating Brewer-Dobson circulation (BDC) is a robust signal of climate change in model predictions but has been questioned by trace gas observations. We analyze stratospheric mean age of air and the full age spectrum as measures for the BDC and its trend. Age of air is calculated with the Chemical Lagrangian Model of the Stratosphere (CLaMS) driven by ERA-Interim, JRA-55 and MERRA-2 reanalysis data to assess the robustness of the representation of the BDC in current generation meteorological reanalyses. We find that climatological mean age significantly depends on the reanalysis, with JRA-55 showing the youngest and MERRA-2 the oldest mean age. Consideration of the age spectrum indicates that the older age for MERRA-2 is related to a stronger spectrum tail, likely related to weaker tropical upwelling and stronger recirculation. Seasonality of stratospheric transport is robustly represented in reanalyses, with similar mean age variations and age spectrum peaks. Long-term changes over 1989–2015 turn out to be similar for the reanalyses with mainly decreasing mean age accompanied by a shift of the age spectrum peak towards shorter transit times, resembling the forced response in climate model simulations to increasing greenhouse gas concentrations. For the shorter periods 1989–2001 and 2002–2015 age of air changes are less robust. Only ERA-Interim shows the hemispheric dipole pattern in age changes during 2002–2015 as viewed by recent satellite observations. Consequently, the representation of decadal variability of the BDC in current generation reanalyses appears less robust and a major uncertainty of modelling the BDC.
The endemic argan tree (Argania spinosa) populations in South Morocco are highly degraded due to their use as a biomass resource in dry years and illegal firewood extraction. The intensification and expansion of agricultural land lead to a retreat of the wooded area, while the remaining argan open woodlands are often overgrazed. Thus, canopy-covered areas decrease while areas without vegetation cover between the argan trees increase. In total, 36 rainfall simulation experiments as well as 60 infiltration measurements were conducted to investigate the potential difference between tree-covered areas and bare intertree areas. In addition, 60 soil samples were taken under the trees and in the intertree areas parallel to the contour lines. Significant differences using a t-test were found between tree and intertree areas for the studied parameters Ksat, Kh, pH, electric conductivity, percolation stability, total C-content, total N-content, K-content, Na-content, and Mg-content. Surface runoff and soil losses were not as conclusive but showed similar trends. The results showed that argan trees influence the soil underneath significantly, while the soil in intertree areas is less protected and more degraded. It is therefore reasonable to assume further degradation of the soil when intertree areas extend further due to lack of rejuvenation of argan trees.
In global hydrological models, groundwater (GW) is typically represented by a bucket-like linear groundwater reservoir. Reservoir models, however, (1) can only simulate GW discharge to surface water (SW) bodies but not recharge from SW to GW, (2) provide no information on the location of the GW table, and (3) assume that there is no GW flow among grid cells. This may lead, for example, to an underestimation of groundwater resources in semiarid areas where GW is often replenished by SW or to an underestimation of evapotranspiration where the GW table is close to the land surface. To overcome these limitations, it is necessary to replace the reservoir model in global hydrological models with a hydraulic head gradient-based GW flow model.
We present G3M, a new global gradient-based GW model with a spatial resolution of 5′ (arcminutes), which is to be integrated into the 0.5∘ WaterGAP Global Hydrology Model (WGHM). The newly developed model framework enables in-memory coupling to WGHM while keeping overall runtime relatively low, which allows sensitivity analyses, calibration, and data assimilation. This paper presents the G3M concept and model design decisions that are specific to the large grid size required for a global-scale model. Model results under steady-state naturalized conditions, i.e., neglecting GW abstractions, are shown. Simulated hydraulic heads show better agreement to observations around the world compared to the model output of de Graaf et al. (2015). Locations of simulated SW recharge to GW are found, as is expected, in dry and mountainous regions but areal extent of SW recharge may be underestimated. Globally, GW discharge to rivers is by far the dominant flow component such that lateral GW flows only become a large fraction of total diffuse and focused recharge in the case of losing rivers, some mountainous areas, and some areas with very low GW recharge. A strong sensitivity of simulated hydraulic heads to the spatial resolution of the model and the related choice of the water table elevation of surface water bodies was found. We suggest to investigate how global-scale groundwater modeling at 5′ spatial resolution can benefit from more highly resolved land surface elevation data.
In global hydrological models, groundwater storages and flows are generally simulated by linear reservoir models. Recently, the first global gradient-based groundwater models were developed in order to improve the representation of groundwater-surface water interactions, capillary rise, lateral flows and human water use impacts. However, the reliability of model outputs is limited by a lack of data as well as model assumptions required due to the necessarily coarse spatial resolution. The impact of data quality is presented by showing the sensitivity of a groundwater model to changes in the only available global hydraulic conductivity data-set. To better understand the sensitivity of model output to uncertain spatially distributed parameter inputs, we present the first application of a global sensitivity method for a global-scale groundwater model using nearly 2000 steady-state model runs of the global gradient-based groundwater model G3M. By applying the Morris method in a novel domain decomposition approach that identifies global hydrological response units, spatially distributed parameter sensitivities are determined for a computationally expensive model. Results indicate that globally simulated hydraulic heads are equally sensitive to hydraulic conductivity, groundwater recharge and surface water body elevation, though parameter sensitivities vary regionally. For large areas of the globe, rivers are simulated to be either losing or gaining, depending on the parameter combination, indicating a high uncertainty of simulating the direction of flow between the two compartments. Mountainous and dry regions show a high variance in simulated head due to numerical difficulties of the model, limiting the reliability of computed sensitivities in these regions. This instability is likely caused by the uncertainty in surface water body elevation. We conclude that maps of spatially distributed sensitivities can help to understand complex behaviour of models that incorporate data with varying spatial uncertainties. The findings support the selection of possible calibration parameters and help to anticipate challenges for a transient coupling of the model.
Understanding new particle formation and growth is important because of the strong impact of these processes on climate and air quality. Measurements to elucidate the main new particle formation mechanisms are essential; however, these mechanisms have to be implemented in models to estimate their impact on the regional and global scale. Parameterizations are computationally cheap ways of implementing nucleation schemes in models, but they have their limitations, as they do not necessarily include all relevant parameters. Process models using sophisticated nucleation schemes can be useful for the generation of look-up tables in large-scale models or for the analysis of individual new particle formation events. In addition, some other important properties can be derived from a process model that implicitly calculates the evolution of the full aerosol size distribution, e.g., the particle growth rates. Within this study, a model (SANTIAGO – Sulfuric acid Ammonia NucleaTIon And GrOwth model) is constructed that simulates new particle formation starting from the monomer of sulfuric acid up to a particle size of several hundred nanometers. The smallest sulfuric acid clusters containing one to four acid molecules and a varying amount of base (ammonia) are allowed to evaporate in the model, whereas growth beyond the pentamer (five sulfuric acid molecules) is assumed to be entirely collision-controlled. The main goal of the present study is to derive appropriate thermodynamic data needed to calculate the cluster evaporation rates as a function of temperature. These data are derived numerically from CLOUD (Cosmics Leaving OUtdoor Droplets) chamber new particle formation rates for neutral sulfuric acid–water–ammonia nucleation at temperatures between 208 and 292 K. The numeric methods include an optimization scheme to derive the best estimates for the thermodynamic data (dH and dS) and a Monte Carlo method to derive their probability density functions. The derived data are compared to literature values. Using different data sets for dH and dS in SANTIAGO detailed comparison between model results and measured CLOUD new particle formation rates is discussed.
Understanding new particle formation and growth is important because of the strong impact of these processes on climate and air quality. Measurements to elucidate the main new particle formation mechanisms are essential; however, these mechanisms have to be implemented in models to estimate their impact on the regional and global scale. Parameterizations are computationally cheap ways of implementing nucleation schemes in models but they have their limitations, as they do not necessarily include all relevant parameters. Process models using sophisticated nucleation schemes can be useful for the generation of look-up tables in large scale models or for the analysis of individual new particle formation events. In addition, some other important properties can be derived from a process model that implicitly calculates the evolution of the full aerosol size distribution, e.g., the particle growth rates. Within this study, a model (SANTIAGO, Sulfuric acid Ammonia NucleaTIon And GrOwth model) is constructed that simulates new particle formation starting from the monomer of sulfuric acid up to a particle size of several hundred nanometers. The smallest sulfuric acid clusters containing one to four acid molecules and varying amount of base (ammonia) are allowed to evaporate in the model, whereas growth beyond the pentamer (5 sulfuric acid molecules) is assumed to be entirely collision-controlled. The main goal of the present study is to derive appropriate thermodynamic data needed to calculate the cluster evaporation rates as a function of temperature. These data are derived numerically from CLOUD (Cosmics Leaving OUtdoor Droplets) chamber new particle formation rates for neutral sulfuric acid-water-ammonia nucleation at temperatures between 208 K and 292 K. The numeric methods include an optimization scheme to derive the best estimates for the thermodynamic data (dH and dS) and a Monte Carlo method to derive their probability density functions. The derived data are compared to literature values. Using different data sets for dH and dS in SANTIAGO detailed comparison between model results and measured CLOUD new particle formation rates is discussed.
The complex magnetotelluric (MT) apparent resistivity tensor can be decomposed into two real tensors, the apparent resistivity and the resistivity phase tensors, which represent relationships between the observed electric field at a point on the Earth's surface and an associated apparent current density. We explain the differences between these tensors and conventional estimates of apparent resistivity and phase for simple resistivity environments and demonstrate, using canonical models in 1‐D and 2‐D environments, that both tensors are more sensitive to vertical and horizontal resistivity gradients than their conventional counterparts. The properties of the new tensors are explained using electromagnetic induction theory and the effects of associated charges at resistivity boundaries. We introduce a new way to plot tensor ellipses, which brings significant improvements to the interpretation of MT data, using appropriate visualization software. The apparent resistivity tensor gives information about the magnitude and direction of apparent resistivity subsurface structures and has a strong response to vertical resistivity contrasts. The resistivity phase tensor is highly sensitive to vertical boundaries and the associated fields in the TM mode. It is also free from static distortions under the same conditions implied for the conventional phase tensor. These findings have prompted a study in the potential of the new tensors for 3‐D inversions. The results from a 3‐D inversion of a canonical oblique conductor straddling two quarter spaces show distinct improvements in resolving the boundaries of the conductor and open a promising field for future studies.
Tetra-auricupride, ideally AuCu, represents the only species showing the coexistence of Au with an elevated level of Pt, as in the case of a detrital grain studied structurally for the first time, from an ophiolite-associated placer at Bolshoy Khailyk, western Sayans, Russia. We infer that tetra-auricupride can incorporate as much as ~30 mol. % of a “PtCu” component, apparently without significant modification of the unit cell. The unit-cell parameters of platiniferous tetra-auricupride are: a 2.790(1) Å, c 3.641(4) Å, with c/a = 1.305, which are close to those reported for ordered AuCu(I) in the system Au–Cu, and close also to the cell parameters of tetraferroplatinum (PtFe), which both appear to crystallize in the same space group, P4/mmm. These intermetallic compounds and natural alloys are thus isostructural. The closeness of their structures presumably allows Pt to replace Au atoms so readily. The high extent of Cu + Au enrichment is considered to be a reflection of geochemical evolution and buildup in levels of the incompatible Cu and Au with subordinate Pt in a remaining volume of melt at low levels of fO2 and fS2 in the system.
Recently, carbonates have attracted a lot of attention, due to the recognition of their importance in the global carbon cycle. This was enabled by improvement of the experimental techniques that allow for investigating the stability, structure, and physical properties of materials and high-pressures and high-temperatures, that is, they allow for investigating minerals and geochemical processes at the conditions occurring deep inside Earth. Although a lot of research has been focused on carbonates, there are still some open questions regarding their structure and physical properties at such extreme conditions. The aim of this thesis is to establish a deeper understanding of the nature of the phase transitions in carbonates by studying how do the atoms building up the crystal structure vibrate, that is lattice dynamics. The methodology adapted in this study is a combination of experimental and computational methods which allows for a very thorough examination of the problem. The computational approach allows to determine parameters that are elusive or tedious to measure, and the experimental results provide a solid benchmark for the calculations. This tandem of methods has been widely used for investigating lattice dynamics of various materials. In this study it was used to elucidate the structure and properties of carbonates in the deep Earth conditions
The main objective of this PhD work is to assess the impact of fine-scale air-sea interaction on the performance of a regional climate prediction model in marginal sea regions. Focus is on the North and Baltic Seas, the largest marginal sea area in the mid-latitudes. Motivation for this work is to better understand the interaction between the different components of the climate system, namely atmosphere, ocean and sea-ice. In addition to that, the sea regions of interest, the North and Baltic Seas, are orographically complex and cannot be resolved by a global ocean model. The ice coverage on the Baltic Sea is underestimated in the stand-alone atmospheric model COSMO-CLM due to the low water freezing temperature value assumed, which is not applicable for such brackish water body. To fulfil the thesis goal, a new regional coupled atmosphere-ocean-ice system was developed for these two seas, named COSMO-CLM/NEMO. The two-way coupling system involves active feedback from both component models: the limited-area climate model COSMO-CLM and the regional ocean model NEMO-NORDIC.
The coupled system COSMO-CLM/NEMO for the North and Baltic Seas was used to study the impact of sea surface temperature and sea ice on the atmosphere on diffrent topics. The long term impact of the North and Baltic Seas was studied through 15- year long simulations driven by European Center for Medium-Range Weather Forecasts (ECMWF) Interim reanalysis (ERA-Interim) data. Furthermore, to see whether the marginal sea modelling can advance the simulation of extreme climate events, the coupled model was used to reproduce six extreme snowband phenomena over the Baltic Sea in simulations driven by ERA-interim data. Last but not least, the role of the North and Baltic Sea model in improving long-term regional climate prediction was examined. Two sets of experiments with coupled and uncoupled models, each set has five independent decadal hindcasts forced by global climate model, were carried out.
All results were compared with observations and the stand-alone atmospheric model COSMO-CLM results. In all experiments, COSMO-CLM/NEMO showed good agreement with observations. Improvements compared with the uncoupled COSMO-CLM were also found. Coupling was found to affect the air temperature not only around the coupled sea region but also inland. The convective snowbands over the Baltic Sea were successfully reproduced by the coupled model. The high contrast of temperature in the air column, as well as considerably high amounts of surface heat fluxes exchanged between air and sea could not be simulated by COSMO-CLM without the help of reanalysis data. The coupled model also provided better forecasts in decadal scales compared with the uncoupled model and the global model. The added predictability came from the initialized regional seas and better simulated sea surface temperatures by the ocean model.
The impact of the North and Baltic Seas on the climate of the surrounding regions is in certain phases dominated by the North Atlantic Oscillation (NAO) activity. In this thesis, the relation between the NAO and the marginal sea influences was studied. It is confirmed by this study that, in strong phases, the NAO can overpower the impact of the local seas. During dominant phases of NAO, the European climate is mainly governed by large-scale circulation. On the other hand, the local seas play an important role in determining the European climate when NAO is in weak phases.
The added value of the coupled model raises promising perspectives for research in this field. It points to a potential benefit of using the coupled atmosphere-ocean-ice system for climate prediction in the region surrounding the North and Baltic Seas. Along with that, it is still a challenge to complete the model representation of the climate system by adding more climate components (such as a hydrological model). Further improvement of the coupled system can be achieved by coupling for a larger sea region, or by trying to reduce remaining low performance of the coupled model in some areas with a better configuration of the current system.
Das Wissen über die Wolkenmikrophysik und die Wechselwirkung zwischen Niederschlag und Aerosol ist ein wichtiger Baustein zur Optimierung von Klima- und Wettermodellen. Ein Großteil des Niederschlags in den mittleren Breiten fällt aus Mischphasenwolken, die aus unterkühlten Tröpfchen und Eispartikeln bestehen. Die Eispartikel bilden sich an speziellen Aerosolpartikeln, die als Eiskeime (INP) wirken können. Die Wahrscheinlichkeit eines Aerosols als Eiskeim zu wirken, nimmt mit abnehmender Temperatur und steigender Wassersättigung zu. Mineralstaubpartikel sind die häufigsten Eiskeime, die ab Temperaturen ≤−15°C aktiv sind, biologische Partikel wirken schon bei wärmeren Bedingungen. Große Wissenslücken bestehen noch bei der globalen Konzentration von Eiskeimen, inklusive deren geographischer und jahreszeitlicher Variabilität.
Im Zentrum der Experimente, die für diese Arbeit durchgeführt wurden, steht der Eiskeimzähler FRIDGE (Frankfurt Ice Deposition Freezing Experiment). Je nach Aufbau und Anwendung des Instruments werden zwei verschiedene Ansätze zur Aktivierung von Eiskeimen verfolgt. Die ursprüngliche und namensgebende Methode in Form einer Vakuum-Diffusionskammer wurde zur Untersuchung von Eisnukleation via Depositionsgefrieren (an INP_D) entwickelt (Klein et al., 2010). Danielczok (2015) nutzte einige Bestandteile des Analysegeräts, um auch Immersionsgefrieren (an INP_I) in Form von Tröpfchengefrieren à la Vali (1971) zu studieren. In der vorliegenden Arbeit wurde diese Anwendungsmöglichkeit von FRIDGE weiterentwickelt. Ein zentraler Schritt war dabei die präzise Charakterisierung des Gefrierverhaltens von Tröpfchen aus Reinstwasser ohne zusätzliches Aerosol. Die Einbeziehung dieses sogenannten Hintergrundgefrierens, das für jedes Instrument und Messverfahren spezifisch ist, ermöglichte es, die Minimaltemperatur, für die zuverlässige Ergebnisse produziert werden können, von −22°C auf −29°C herabzusetzen. Der dadurch hinzugewonnene Temperaturbereich ist für Eisnukleation in Mischphasenwolken äußerst relevant.
Beide Anwendungsmethoden wurden im Rahmen des Fifth International Workshop on Ice Nucleation – Phase 2 (FIN-02) sowie bei einer weiteren Kampagne zur Messung von Eisnukleation an Cellulosepartikeln mit über zwanzig anderen Eiskeimzählern verglichen. Mit FRIDGE als Diffusionskammer wurde für die Mehrheit der untersuchten Aerosoltypen eine zufriedenstellende Übereinstimmung mit den anderen Instrumenten beobachtet. Die Experimente mit gefrierenden Tröpfchen in FRIDGE erzielten ausnahmslos Ergebnisse, die inmitten der Werte der anderen Instrumente lagen. Die erfolgreiche Validierung – besonders der neuen Anwendungsmethode – war das erste Ziel dieser Arbeit und die notwendige Voraussetzung für die anschließenden Feldmessungen.
Atmosphärische Eiskeimkonzentrationen wurden in mehrwöchigen Feldmesskampagnen an drei sehr unterschiedlichen Orten und atmosphärischen Bedingungen untersucht: an der Hochalpinen Forschungsstation Jungfraujoch (JFJ), am Storm Peak Laboratory (SPL) in den Rocky Mountains und am Cyprus Atmospheric Observatory (CAO) in Zypern. Am JFJ wurde eine INP-Konzentration beobachtet, die um den Faktor 20 niedriger als an den anderen beiden Stationen war. Der Grund dafür war, dass sich das Jungfraujoch die meiste Zeit der Messungen in der freien Troposphäre befand. Dementsprechend waren die Bedingungen an der Station von aerosolpartikelarmer Luft mit wenigen Eiskeimen geprägt. An zwei Standorten wurde Mineralstaub als ein Parameter, der die lokale INP-Konzentration positiv beeinflusst, identifiziert. Sowohl am JFJ als auch am CAO erhöhte ferntransportierter Saharastaub die INP-Konzentration.
Die Kombination der zwei Analysemethoden, die Aerosolproben aus derselben Luft entweder in trockenem Ausgangszustand oder in Tröpfchen suspendiert untersuchen, offenbarte eine interessante Eigenschaft der INP. Es herrschte eine offensichtliche Parallelität von INP_D- und INP_I- Konzentrationen. Bei fast allen Messungen waren INP_I etwa 10-mal häufiger als INP_D. Die Aussage trifft gleichermaßen auf sehr niedrige Konzentrationen am JFJ wie auch auf hohe Konzentrationen am SPL und CAO zu. Die einzige Ausnahme bilden Cellulosepartikel. Daraus lässt sich schließen, dass INP_D und INP_I dieselben Partikel sind, die je nach Aktivierungskonditionen Eisnukleation unterschiedlich effektiv hervorrufen.
Wetlands such as bogs, swamps, or freshwater marshes are hotspots of biodiversity. For 5.1 million km2 of inland wetlands, the dynamics of area and water storage, which strongly impact biodiversity and ecosystem services, were simulated using the global hydrological model WaterGAP. For the first time, the impacts of both human water use and man‐made reservoirs (WUR) and future climate change (CC) on wetlands around the globe were quantified. WUR impacts are concentrated in arid/semiarid regions, where WUR decreased mean wetland water storage by more than 5% on 8.2% of the mean wetland area during 1986–2005 (Am), with highest decreases in groundwater depletion area. Using output of three climate models, CC impacts on wetlands were quantified, distinguishing unavoidable impacts [i.e., at 2 °C global warming (GW)] from avoidable impacts (difference between 3 °C and 2 °C impacts). Even unavoidable CC impacts are projected to be much larger than WUR impacts, also in arid/semiarid regions. On most wetland area with reliable estimates, avoidable CC impacts are more than twice as large as unavoidable impacts. In case of 2 °C GW, half of Am is estimated to be unaffected by mean storage changes of more than 5%, but only one third in case of 3 °C GW. Temporal variability of water storage will increase for most wetlands. Wetlands in dry regions will be affected the most, particularly by water storage decreases in the dry season. Different from wealthier countries, low‐income countries will dominantly suffer from a decrease in wetland water storage due to CC.
State-of-the-art climate models contain, to a significant degree, empirical components. In particular, subgrid-scale (SGS) parameterizations are usually highly tuned against observations or high-resolution model data. While this enables the models to minimize the error during hindcasts, it is not guaranteed that it yields a benefit for climate projections because of climate change. In this thesis the Fluctuation-Dissipation theorem (FDT) is used to update the statistics of the system in the presence of an external forcing. If the empirical parameters are tuned objectively to the data (i.e., they depend on the statistics of the data), then they might be updated with the FDT. This ansatz is tested within a framework of a semi-empirical model (SEM) based on the leading variance patterns of a quasigeostrophic three-layer model (QG3LM) and supplemented by a purely data-driven parameterization. We show that the FDT is able to successfully update the tuning parameters of the data-driven SGS closure, resulting in a systematic improvement in model performance in comparison to an untreated SEM. Ideally, SGS parameterizations should contain little to no tuning parameters. Thus, complementary to the FDT approach we investigate a stochastic SGS closure constrained by first principles that is calculated using the stochastic mode reduction (SMR). The SMR allows for an analytic derivation of the SGS closure from the model equations while requiring only minimal tuning. We successfully apply the SMR to the QG3LM and construct the reduced stochastic model (RSM). Furthermore, we show that the RSM is more robust against an external forcing than the SEM. Additionally, we find that, under appropriate conditions, the FDT is able to update the empirical parts of the RSM. Yet, only for the response in mean streamfunction the RSM provides useful results, while the response in covariance of the streamfunction is incorrect for most cases. Nevertheless, we obtain a remarkably accurate response in both moments for the RSM in an idealized setting. In combination with the results of the FDT study this indicates that the considered RSM is too low dimensional and encourages us to investigate the response of larger RSMs in the future.
The continental expression of global cooling during the Miocene Climate Transition in Central Asia is poorly documented, as the tectonically active setting complicates the correlation of Neogene regional and global climatic developments. This study presents new geochemical data (CaSO4 content, carbonate δ13C and δ18O) from the endorheic alluvial‐lacustrine Aktau succession (Ili Basin, south‐east Kazakhstan) combined with findings from the previously published facies evolution. Time series analysis revealed long‐eccentricity forcing of the paleohydrology throughout the entire succession, split into several facies‐dependent segments. Orbital tuning, constrained by new laser ablation U‐Pb dates and a preexisting magnetostratigraphy, places the succession in a 5.0 Ma long interval in the middle to late Miocene (15.6 to 10.6 Ma). The long‐term water accumulation in the Ili Basin followed the timing of the Miocene Climate Transition, suggesting increased precipitation in the catchment area in response to climate cooling and stronger westerly winds. This was paced by minima of the 2.4 Ma eccentricity cycle, which favored the establishment of a discharge playa (~14.3 Ma) and a perennial lake (12.6 to 11.8 Ma). Furthermore, low obliquity amplitudes (nodes) caused a transient weakening of the westerlies at ~13.7 to 13.5 Ma and at ~12.7 Ma, resulting in negative hydrological budgets and salinization. Flooding of the windward Ili Basin coeval with aridification in the leeward basins suggests that the Tian Shan was a climate boundary already in the middle Miocene. Our results emphasize the impact of climate fluctuations on the westerlies' strength and thus on Central Asian hydrology.
Un tuf volcanique a ete mis en evidence dans trois profils de loess de Hesbaye (Rocourt et Lixhe) et du Limbourg neerlandais (Nagelbeek). Ce tuf correspond probablement ä l'Eltviller Tuff connu en Hesse (SEMMEL 1967). En consequence, cette decouverte permet de comparer directement une partie de la sequence des loess du Pleistocene superieur propre ä chaque region
Within the world’s oceans, regionally distinct ecological niches develop due to differences in water temperature, nutrients, food availability, predation and light intensity. This results in differences in the vertical dispersion of planktonic foraminifera on the global scale. Understanding the controls on these modern-day distributions is important when using these organisms for paleoceanographic reconstructions. As such, this study constrains modern depth habitats for the northern equatorial Indian Ocean, for 14 planktonic foraminiferal species (G. ruber, G. elongatus, G. pyramidalis, G. rubescens, T. sacculifer, G. siphonifera, G. glutinata, N. dutertrei, G. bulloides, G. ungulata, P. obliquiloculata, G. menardii, G. hexagonus, G. scitula) using stable isotopic signatures (δ18O and δ13C) and Mg/Ca ratios. We evaluate two aspects of inferred depth habitats: (1) the significance of the apparent calcification depth (ACD) calculation method/equations and (2) regional species-specific ACD controls. Through a comparison with five global, (sub)tropical studies we found the choice of applied equation and δ18Osw significant and an important consideration when comparing with the published literature. The ACDs of the surface mixed layer and thermocline species show a tight clustering between 73–109 m water depth coinciding with the deep chlorophyll maximum (DCM). Furthermore, the ACDs for the sub-thermocline species are positioned relative to secondary peaks in the local primary production. We surmise that food source plays a key role in the relative living depths for the majority of the investigated planktonic foraminifera within this oligotrophic environment of the Maldives and elsewhere in the tropical oceans.
Here we present a comprehensive attempt to correlate aragonitic Na∕Ca ratios from Desmophyllum pertusum (formerly known as Lophelia pertusa), Madrepora oculata and a caryophylliid cold-water coral (CWC) species with different seawater parameters such as temperature, salinity and pH. Living CWC specimens were collected from 16 different locations and analyzed for their Na∕Ca ratios using solution-based inductively coupled plasma-optical emission spectrometry (ICP-OES) measurements.
The results reveal no apparent correlation with salinity (30.1–40.57 g kg−1) but a significant inverse correlation with temperature (−0.31±0.04 mmolmol−1∘C−1). Other marine aragonitic organisms such as Mytilus edulis (inner aragonitic shell portion) and Porites sp. exhibit similar results highlighting the consistency of the calculated CWC regressions. Corresponding Na∕Mg ratios show a similar temperature sensitivity to Na∕Ca ratios, but the combination of two ratios appears to reduce the impact of vital effects and domain-dependent geochemical variation. The high degree of scatter and elemental heterogeneities between the different skeletal features in both Na∕Ca and Na∕Mg, however, limit the use of these ratios as a proxy and/or make a high number of samples necessary. Additionally, we explore two models to explain the observed temperature sensitivity of Na∕Ca ratios for an open and semi-enclosed calcifying space based on temperature-sensitive Na- and Ca-pumping enzymes and transport proteins that change the composition of the calcifying fluid and consequently the skeletal Na∕Ca ratio.
The Arctic Svalbard Archipelago hosts the world’s northernmost cold-water ‘carbonate factories’ thriving here despite of presumably unfavourable environmental conditions and extreme seasonality. Two contrasting sites of intense biogenic carbonate production, the rhodolith beds in Mosselbukta in the north of the archipelago and the barnacle-mollusc dominated carbonate sediments accumulating in the strong hydrodynamic regime of the Bjørnøy-Banken south of Spitsbergen, were the targets of the RV Maria S. Merian cruise 55 in June 2016. By integrating data from physical oceanography, marine biology, and marine geology, the present contribution characterises the environmental setting and biosedimentary dynamics of these two polar carbonate factories. Repetitive CTD profiling in concert with autonomous temperature/salinity loggers on a long-term settlement platform identified spatiotemporal patterns in the involved Atlantic and Polar water masses, whereas short-term deployments of a lander revealed fluctuations of environmental variables in the rhodolith beds in Mosselbukta and at same depth (46 m) at Bjørnøy-Banken. At both sites, dissolved inorganic nutrients in the water column were found depleted (except for elevated ammonium concentrations) and show an overall increase in concentration and N:P ratios toward deeper waters. This indicates that a recycling system was fuelling primary production after the phytoplankton spring bloom at the time of sampling in June 2016. Accordingly, oxygen levels were found elevated and carbon dioxide concentrations (pCO2) markedly reduced, on average only half the expected equilibrium values. Backed up by seawater stable carbon and oxygen isotope signatures, this is interpreted as an effect of limited air-sea gas exchange during seasonal ice cover in combination with a boost in community photosynthesis during the spring phytoplankton bloom. The observed trends are enhanced by the onset of rhodophyte photosynthesis in the rhodolith beds during the polar day upon retreat of sea-ice. Potential adverse effects of ocean acidification on the local calcifier community are thus predicted to be seasonally buffered by the marked drop in pCO2 during the phase of sea-ice cover and spring phyto-plankton bloom, but this effect will diminish should the seasonal sea-ice formation continue to decline. Among the 25 macrobenthos taxa identified from images captured by the lander’s camera system, all but three species were calcifiers contributing to the carbonate production. Biodiversity was found to be much higher in Mosselbukta (21 taxa) compared to Bjørnøy-Banken (8 taxa), which is considered as a result of enhanced habitat diversity provided in the rhodolith beds by the bioengineering crustose alga Lithothamnion glaciale. Filter-feeding activity of selected key species did reveal group-specific but no common activity patterns. Biotic disturbance of the filtering activity was common, in contrast to abiotic factors, with hermit crabs representing the primary trigger. Motion tracking of rhodoliths revealed a high frequency of dislocation, triggered not by abiotic factors but by the activity of benthic invertebrates, in particular echinoids ploughing below or moving over the rhodoliths. The echinoid Strongylocentrotus sp. is the most abundant component of the associated fauna, thereby considerably contributing both to carbonate production and to grazing bioerosion. Together, these results portray a high degree of seasonal as well as short-term dynamics in environmental conditions that despite many similarities support distinctly different communities and biodiversity patterns in the calcifying macrobenthos at the two studied polar carbonate factories.
Here we present a comprehensive attempt to correlate aragonitic Na / Ca ratios from Lophelia pertusa, Madrepora oculata and a caryophylliid cold-water coral (CWC) species with different seawater parameters such as temperature, salinity and pH. Living CWC specimens were collected from 16 different locations and analyzed for their Na / Ca content using solution-based inductively coupled plasma-optical emission spectrometry (ICP-OES) measurements. The results reveal no apparent correlation with salinity (30.1–40.57 g/kg) but a significant inverse correlation with temperature (−0.31 mmol/mol/°C). Other marine aragonitic organisms such as Mytilus edulis (inner aragonitic shell portion) and Porites sp. exhibit similar results highlighting the consistency of the calculated CWC regressions. Corresponding Na / Mg ratios show a similar temperature sensitivity to Na / Ca ratios, but the combination of two ratios appear to reduce the impact of vital effects and domain-dependent geochemical variation. The high degree of scatter and elemental heterogeneities between the different skeletal features in both Na / Ca and Na / Mg however limit the use of these ratios as a proxy and/or make a high number of samples necessary. Additionally, we explore two models to explain the observed temperature sensitivity of Na / Ca ratios for an open and semi-enclosed calcifying space based on temperature sensitive Na and Ca pumping enzymes and transport proteins that change the composition of the calcifying fluid and consequently the skeletal Na / Ca ratio.
Bioaerosols are considered to play a relevant role in atmospheric processes, but their sources, properties and spatiotemporal distribution in the atmosphere are not yet well characterized. In the Amazon Basin, primary biological aerosol particles (PBAP) account for a large fraction of coarse particulate matter, and fungal spores are among the most abundant PBAP there as well as in other vegetated continental regions. furthermore, PBAP could also be important ice nuclei in Amazonia. Measurement data on the release of fungal spores under natural conditions, however, are sparse. Here we present an experimental approach to analyze and quantify the spore release from fungi and other spore producing organisms under natural and laboratory conditions. For measurements under natural conditions, the samples were kept in their natural environment and a setup was developed to estimate the spore release numbers and sizes together with the microclimatic factors temperature and air humidity, as well as the mesoclimatic parameters net radiation, rain, and fog occurrence. For experiments in the laboratory, we developed a cuvette to assess the particle size and number of newly released fungal spores under controlled conditions, simultaneously measuring temperature and relative humidity inside the cuvette. Both approaches were combined with bioaerosol sampling techniques to characterize the released particles by microscopic methods. For fruiting bodies of the basidiomycetous species, Rigidoporus microporus, the model species for which these techniques were tested, the highest frequency of spore release occurred in the range of 62 and 96 % relative humidity. The results obtained for this model species reveal characteristic spore release patterns linked to environmental or experimental conditions, indicating that the moisture status of the sample may be a regulating factor, while temperature and light seem to play a minor role for this species. The presented approach enables systematic studies aimed at the quantification and validation of spore emission rates and inventories, which can be applied to a regional mapping of cryptogamic organisms under given environmental conditions.
Latitudinal and bathymetrical species richness patterns in the NW Pacific and adjacent Arctic Ocean
(2019)
Global scale analyses have recently revealed that the latitudinal gradient in marine species richness is bimodal, peaking at low-mid latitudes but with a dip at the equator; and that marine species richness decreases with depth in many taxa. However, these overall and independently studied patterns may conceal regional differences that help support or qualify the causes in these gradients. Here, we analysed both latitudinal and depth gradients of species richness in the NW Pacific and its adjacent Arctic Ocean. We analysed 324,916 distribution records of 17,414 species from 0 to 10,900 m depth, latitude 0 to 90°N, and longitude 100 to 180°N. Species richness per c. 50 000 km2 hexagonal cells was calculated as alpha (local average), gamma (regional total) and ES50 (estimated species for 50 records) per latitudinal band and depth interval. We found that average ES50 and gamma species richness decreased per 5° latitudinal bands and 100 m depth intervals. However, average ES50 per hexagon showed that the highest species richness peaked around depth 2,000 m where the highest total number of species recorded. Most (83%) species occurred in shallow depths (0 to 500 m). The area around Bohol Island in the Philippines had the highest alpha species richness (more than 8,000 species per 50,000 km2). Both alpha and gamma diversity trends increased from the equator to latitude 10°N, then further decreased, but reached another peak at higher latitudes. The latitudes 60–70°N had the lowest gamma and alpha diversity where there is almost no ocean in our study area. Model selection on Generalized Additive Models (GAMs) showed that the combined effects of all environmental predictors produced the best model driving species richness in both shallow and deep sea. The results thus support recent hypotheses that biodiversity, while highest in the tropics and coastal depths, is decreasing at the equator and decreases with depth below ~2000 m. While we do find the declines of species richness with latitude and depth that reflect temperature gradients, local scale richness proved poorly correlated with many environmental variables. This demonstrates that while regional scale patterns in species richness may be related to temperature, that local scale richness depends on a greater variety of variables.
Increasing atmospheric CO2 stimulates photosynthesis which can increase net primary production (NPP), but at longer timescales may not necessarily increase plant biomass. Here we analyse the four decade-long CO2-enrichment experiments in woody ecosystems that measured total NPP and biomass. CO2 enrichment increased biomass increment by 1.05 ± 0.26 kg C m−2 over a full decade, a 29.1 ± 11.7% stimulation of biomass gain in these early-secondary-succession temperate ecosystems. This response is predictable by combining the CO2 response of NPP (0.16 ± 0.03 kg C m−2 y−1) and the CO2-independent, linear slope between biomass increment and cumulative NPP (0.55 ± 0.17). An ensemble of terrestrial ecosystem models fail to predict both terms correctly. Allocation to wood was a driver of across-site, and across-model, response variability and together with CO2-independence of biomass retention highlights the value of understanding drivers of wood allocation under ambient conditions to correctly interpret and predict CO2 responses.
Cenozoic lignite deposits are widespread across Europe, Asia, America, Australia, and Indonesia. These deposits were the subject of numerous studies on changes in regional/global paleoclimates, paleobotany, paleoenvironment, and basin evolutions, which led to the formation of these lignites. In some of these Cenozoic lignite deposit basins, a succession of pale and dark lignite layers has been described in the Miocene Lower Rhine Basin in Germany, the Oligo-Miocene Gippsland Basin in southeastern Australia, and several Mio-Pliocene basins in southwestern China. Furthermore, pale and dark lithotypes in lignite seams also have been found in some Pliocene lignite deposit basins from Slovenia, Serbia, and Poland. The widespread cyclic occurrence of pale and dark layers in lignite basins might represent alternating depositional conditions related to the changes in plant communities, the regional/global climate, the tectonic setting, the Asian monsoon, and orbital periodicity during peat formation. ...
Abiotic formation of n-alkane hydrocarbons has been postulated to occur within Earth's crust. Apparent evidence was primarily based on uncommon carbon and hydrogen isotope distribution patterns that set methane and its higher chain homologues apart from biotic isotopic compositions associated with microbial production and closed system thermal degradation of organic matter. Here, we present the first global investigation of the carbon and hydrogen isotopic compositions of n-alkanes in volcanic-hydrothermal fluids hosted by basaltic, andesitic, trachytic and rhyolitic rocks. We show that the bulk isotopic compositions of these gases follow trends that are characteristic of high temperature, open system degradation of organic matter. In sediment-free systems, organic matter is supplied by surface waters (seawater, meteoric water) circulating through the reservoir rocks. Our data set strongly implies that thermal degradation of organic matter is able to satisfy isotopic criteria previously classified as being indicative of abiogenesis. Further considering the ubiquitous presence of surface waters in Earth’s crust, abiotic hydrocarbon occurrences might have been significantly overestimated.
During my PhD, I was applying the clumped isotope technique to modern brachiopods and fossil belemnites, and I conducted methodological work. Carbonate clumped isotope thermometry is a tool to reconstruct carbonate precipitation temperatures. In contrast to oxygen isotope thermometry, i.e., the δ18O-thermometer, the carbonate clumped isotope thermometer does not require an estimate for the oxygen isotope composition of the seawater, as it considers the fractionation of isotopes exclusively amongst carbonate isotopologues. The ∆47 value of a carbonate expresses the abundance of the 13C–18O bond bearing carbonate isotopologue, within the carbonate, relative to its random distribution. In thermodynamic equilibrium, the ∆47 value of a given carbonate is solely a function of the carbonate precipitation temperature. However, kinetic isotope fractionations, i.e., vital effects, driven by diffusion, pH or incomplete oxygen isotope exchange between water and dissolved inorganic carbonate species can cause the carbonate to be precipitated with isotopic compositions that are offset from those predicted for thermodynamic equilibrium.
Brachiopods serve as important geochemical archives of past climate conditions. To investigate the nature and significance of kinetic controls on brachiopod shell δ18O and ∆47 values, in collaboration with the BASE-LiNE Earth ITN, I analysed the bulk and clumped isotope compositions of eighteen modern brachiopod shells, collected from different geographic locations and water depths that cover a substantial range of growth temperatures. Growth temperatures and seawater δ18O values for each brachiopod were independently determined. Most of the analysed brachiopods exhibit combined offsets from clumped and oxygen isotope equilibrium, and there is a significant negative correlation between the offset values. The observed correlation slope between offset ∆47 and offset δ18O point to the importance of kinetic effects associated with Knudsen diffusion and incomplete hydration and hydroxylation of CO2 (aq), occurring during biomineralisation. The correlations between the growth rates of the analysed brachiopods and both the offset ∆47 and the offset δ18O values provide further arguments for the presence of kinetic effects. In conclusion, the oxygen and clumped isotope composition of modern brachiopod shells are affected by growth rate-induced kinetic effects that hinder their use for palaeoceanography.
Video and image data are regularly used in the field of benthic ecology to document biodiversity. However, their use is subject to a number of challenges, principally the identification of taxa within the images without associated physical specimens. The challenge of applying traditional taxonomic keys to the identification of fauna from images has led to the development of personal, group, or institution level reference image catalogues of operational taxonomic units (OTUs) or morphospecies. Lack of standardisation among these reference catalogues has led to problems with observer bias and the inability to combine datasets across studies. In addition, lack of a common reference standard is stifling efforts in the application of artificial intelligence to taxon identification. Using the North Atlantic deep sea as a case study, we propose a database structure to facilitate standardisation of morphospecies image catalogues between research groups and support future use in multiple front-end applications. We also propose a framework for coordination of international efforts to develop reference guides for the identification of marine species from images. The proposed structure maps to the Darwin Core standard to allow integration with existing databases. We suggest a management framework where high-level taxonomic groups are curated by a regional team, consisting of both end users and taxonomic experts. We identify a mechanism by which overall quality of data within a common reference guide could be raised over the next decade. Finally, we discuss the role of a common reference standard in advancing marine ecology and supporting sustainable use of this ecosystem.
We present novel measurements of five short-lived brominated source gases (CH2Br2, CHBr3, CH2ClBr, CHCl2Br and CHClBr2) obtained using a gas chromatograph-mass spectrometer system on board the High Altitude and Long Range Research Aircraft (HALO). The instrument is extremely sensitive due to the use of chemical ionisation, allowing detection limits in the lower parts per quadrillion (10-15) range. Data from three campaigns using the HALO aircraft are presented, where the Upper Troposphere/Lower Stratosphere (UTLS) of the Northern Hemisphere mid to high latitudes were sampled during winter and during late summer to early fall. We show that an observed decrease with altitude in the stratosphere is consistent with the relative lifetimes of the different compounds. Distributions of the five source gases and total organic bromine just below the tropopause shows an increase in mixing ratio with latitude, in particular during polar winter. This increase in mixing ratio is explained by increasing lifetimes at higher latitudes during winter. As the mixing ratio at the extratropical tropopause are generally higher than those derived for the tropical tropopause, extratropical troposphere-to-stratosphere transport will result in elevated levels of organic bromine in comparison to air transported over the tropical tropopause. The observations are compared to model estimates using different emission scenarios. A scenario which has emissions most strongly concentrated to low latitudes cannot reproduce the observed latitudinal distributions and will tend to overestimate bromine input through the tropical tropopause from CH2Br2 and CHBr3. Consequently, the scenario also overestimates the amount of brominated organic gases in the stratosphere. The two scenarios with the highest overall emissions of CH2Br2 tend to overestimate mixing ratios at the tropical tropopause but are in much better agreement with extratropical tropopause values, showing that not only total emissions but also latitudinal distributions in the emissions are of importance. While an increase in tropopause values with latitude is reproduced with all emission scenarios during winter, the simulated extratropical tropopause values are on average lower than the observations during late summer to fall. We show that a good knowledge of the latitudinal distribution of tropopause mixing ratios and of the fractional contributions of tropical and extratropical air is needed to derive stratospheric inorganic bromine in the lowermost stratosphere from observations. Depending on the underlying emission scenario, differences of a factor 2 in reactive bromine derived from observations and model outputs are found for the lowermost stratosphere, based on source gas injection. We conclude that a good representation of the contributions of different source regions is required in models for a robust assessment of the role of short-lived halogen source gases on ozone depletion in the UTLS.
Stability of halocarbons in whole air samples from the upper
troposphere and lowermost stratosphere
(2019)
Measurements of halogenated hydrocarbons of ambient air frequently rely on canister sampling followed by offline laboratory analysis. This allows for a large number of compounds to be analysed under stable conditions, maximising measurement precision. However, individual compounds might be affected during storage of canister samples. In order to assess halocarbon stability in whole air samples from the upper troposphere and lowermost stratosphere, we performed stability tests using the air sampling unit High REsolution Sampler (HIRES) which is part of the CARIBIC (Civil Aircraft for the Regular Investigation of the Atmosphere Based on an Instrument Container) instrument package. HIRES holds 88 light-weight stainless steel cylinders that are pressurized in flight to 4.5 bar using metal bellows pumps. The HIRES sampling unit was first deployed in 2010, but has up to now not been used for regular halocarbon analysis with exception of chloromethane. The sample collection unit was tested for sampling and storage effects of 28 halogenated compounds. The focus was on compound stability in the stainless steel canisters during storage of up to five weeks and on the influence of ozone, since flights take place in the upper troposphere and lowermost stratosphere with ozone mixing ratios of up to several hundred ppbV. Most of the investigated (hydro)chlorofluorocarbons and long-lived hydrofluorocarbons were found to be stable over a storage time of up to five weeks and were unaltered by ozone being present during pressurization. Some compounds such as for example dichloromethane, trichloromethane and tetrachloroethene started to decrease in the canisters after a storage time of more than two weeks or exhibited lowered mixing ratios in samples pressurized with ozone present. Few compounds such as for example tetrachloromethane and tribromomethane were found to be not stable in the HIRES stainless-steel canisters independent of ozone levels. Also growth was observed during storage, namely for HFC-152a and HFC-23.
Ecophysiological studies on Antarctic cryptophytes to assess whether climatic changes such as ocean acidification and enhanced stratification affect their growth in Antarctic coastal waters in the future are lacking so far. This is the first study that investigated the combined effects of increasing availability of pCO2 (400 and 1000 µatm) and irradiance (20, 200 and 500 μmol photons m−2 s −1) on growth, elemental composition and photophysiology of the Antarctic cryptophyte Geminigera cryophila. Under ambient pCO2, this species was characterized by a pronounced sensitivity to increasing irradiance with complete growth inhibition at the highest light intensity. Interestingly, when grown under high pCO2 this negative light effect vanished and it reached highest rates of growth and particulate organic carbon production at the highest irradiance compared to the other tested experimental conditions. Our results for G. cryophila reveal beneficial effects of ocean acidification in conjunction with enhanced irradiance on growth and photosynthesis. Hence, cryptophytes such as G. cryophila may be potential winners of climate change, potentially thriving better in more stratified and acidic coastal waters and contributing in higher abundance to future phytoplankton assemblages of coastal Antarctic waters.
Highlights
• Subcrustal earthquakes detected beneath Fogo volcano, Cape Verde.
• At the focal depth of 40 km temperatures are likely too high for brittle failure.
• The earthquakes may originate from magma injection into a deep subcrustal reservoir.
• This observation indicates a distinct magma supply system of Fogo volcano.
Abstract
Fogo volcano belongs to the Cape Verde hotspot and its most recent eruption occurred from November 2014 to February 2015. From January to December 2016 we operated a temporary seismic network and array on Fogo and were able to locate 289 earthquakes in total. Array analysis shows that most of the events occur within the crust at distances >25 km near the neighboring island of Brava. However, on 15th August 2016 the network recorded an isolated cluster of >20 earthquakes, 13 of which could be located beneath the southern part of Fogo. The differences between S- and P-wave arrival times at steep incidence clearly indicate focal depths between approximately 38 and 44 km whereas receiver-function analyses place the Moho discontinuity at depths between 11 and 14 km. Thus, the earthquakes are located well within the upper mantle directly beneath Fogo. In view of the elevated upper-mantle temperatures within a hotspot regime, we propose that fracturing induced by magmatic injection is the most likely cause for the observed deep earthquakes.
Diamonds growing in the Earth’s mantle often trap inclusions of fluids that are highly saline in composition. These fluids are thought to emerge from deep in subduction zones and may also be involved in the generation of some of the kimberlite magmas. However, the source of these fluids and the mechanism of their transport into the mantle lithosphere are unresolved. Here, we present experimental results showing that alkali chlorides are stable solid phases in the mantle lithosphere below 110 km. These alkali chlorides are formed by the reaction of subducted marine sediments with peridotite and show identical K/Na ratios to fluid inclusions in diamond. At temperatures >1100°C and low pressures, the chlorides are unstable; here, potassium is accommodated in mica and melt. The reaction of subducted sediments with peridotite explains the occurrence of Mg carbonates and the highly saline fluids found in diamonds and in chlorine-enriched kimberlite magmas.
2019 hat Frankfurt einen Rekord geknackt: Am 25. Juli wurde im Stadtteil Westend die Tageshöchsttemperatur von 40,2 Grad Celsius erreicht. Damit war Frankfurt der heißeste Ort Hessens seit Beginn der Wetteraufzeichnungen im Jahr 1881. Ein Superlativ, der zeigt: Der Klimawandel stellt auch die Städte vor große Herausforderungen. Gesucht werden Wege, um mit seinen Folgen zurechtzukommen.
Precipitation extremes with devastating socioeconomic consequences within the South American Monsoon System (SAMS) are expected to become more frequent in the near future. The complexity in SAMS behavior, however, poses severe challenges for reliable future projections. Thus, robust paleomonsoon records are needed to constrain the high spatiotemporal variability in the response of SAMS rainfall to different climatic drivers. This study uses Ti/Ca ratios from X‐ray fluorescence scanning of a sediment core retrieved off eastern Brazilian to trace precipitation changes over the past 322 Kyr. The results indicate that despite the spatiotemporal complexity of the SAMS, insolation forcing is the primary pacemaker of variations in the monsoonal system. Additional modulation by atmospheric pCO2 suggests that SAMS intensity over eastern Brazil will be suppressed by rising CO2 emissions in the future. Lastly, our record reveals an unprecedented strong and persistent wet period during Marine Isotope Stage 6 driven by anomalously strong trade winds.
The complexity of atmospherical processes has always yielded a multitude of ways of knowing about the weather. What has been lacking in the historiography of meteorology so far is a way to formulate differences between forms of knowledge in a way that does not privilege modern scientific structures, but focuses instead on the epistemological category of causality. Using causality as ground of comparison for different knowledge claims, I shall argue, may enable researchers to investigate meteorological knowledge across time periods, perhaps even geographical regions, in a more symmetrical manner. This review demonstrates this approach as a means to organize a large set of historical meteorological writings from German countries between 1750 and 1850. Three distinct forms of knowledge (Semiotics, Physics, and Organics of the weather) during that time and in that region are suggested and will be described. While a bibliography with a national perspective from the 1880s was the basis for the selection of historical sources, such a setup proved awkward even to contemporaries. In addition, the bibliography came with a number of biases and shortcomings that will be critically reviewed.
WaterGAP (Water - Global Assessment and Prognosis) is a tool for modeling global water use and water availability. It participates among other models in the ISIMIP initiative (The Inter-Sectoral Impact Model Intercomparison Project). As part of this initiative, the water temperature should be calculated by participating hydrological models because it plays a vital role in many chemical, physical and biological processes. Therefore, the subject of this master thesis is to implement the physically based surface water temperature computation after VAN BEEK ET AL. (2012) and WANDERS ET AL. (2019) into WaterGAP and compare the results to the statistical regression approach by PUNZET ET AL. (2012). The computation is validated with observed water temperature data obtained from the GEMStat water quality database. The results are good for arctic and temperate latitudes. Surface water temperatures for tropical rivers are overestimated, most likely due to the overestimation of precipitation temperatures, incoming radiation and groundwater temperatures. The comparison with the regression model by PUNZET ET AL. (2012) shows matching results. The regression model even matches with WaterGAP results for most of the simulations of the future under climate change conditions, where the regression model should stop working due to changing environmental parameters. Several assumptions had to be made in order to implement the water temperature calculation in Water-GAP. These include, e.g., discharge temperatures for power plant cooling water, precipitation and surface runoff temperatures. For model improvements, perhaps three different values for the different regions of the world should be used to cool down the precipitation and surface runoff. The model could also be improved by refining the ice formation calculation, especially for the conditions when the ice melts, breaks up and is transported downstream. Furthermore, the feedback to the river channel roughness could be implemented if ice has formed. The WaterGAP model upgraded with the water temperature calculation will help the ISIMIP initiative in the future.
Abstract
The mineralogy, chemical composition, and physical properties of cratonic mantle eclogites with oceanic crustal protoliths can be modified by secondary processes involving interaction with fluids and melts, generated in various slab lithologies upon subduction (auto‐metasomatism) or mantle metasomatism after emplacement into the cratonic lithosphere. Here we combine new and published data to isolate these signatures and evaluate their effects on the chemical and physical properties of eclogite. Mantle metasomatism involving kimberlite‐like, ultramafic carbonated melts (UM carbonated melts) is ubiquitous though not pervasive, and affected between ~20% and 40% of the eclogite population at the various localities investigated here, predominantly at ~60–150 km depth, overlapping cratonic midlithospheric seismic discontinuities. Its hallmarks include lower jadeite component in clinopyroxene and grossular component in garnet, an increase in bulk‐rock MgO ± SiO2, and decrease in FeO and Al2O3 contents, and LREE‐enrichment accompanied by higher Sr, Pb, Th, U, and in part Zr and Nb, as well as lower Li, Cu ± Zn. This is mediated by addition of a high‐temperature pyroxene from a UM carbonated melt, followed by redistribution of this component into garnet and clinopyroxene. As clinopyroxene‐garnet trace‐element distribution coefficients increase with decreasing garnet grossular component, clinopyroxene is the main carrier of the metasomatic signatures. UM carbonated melt‐metasomatism at >130–150 km has destroyed the diamond inventory at some localities. These mineralogical and chemical changes contribute to low densities, with implications for eclogite gravitational stability, but negligible changes in shear‐wave velocities, and, if accompanied by H2O‐enrichment, will enhance electrical conductivities compared to unenriched eclogites.
Plain Language Summary
Oceanic crust formed at spreading ridges is recycled in subduction zones and undergoes metamorphism to eclogite. Some of this material is captured in the overlying lithospheric mantle, where it is exhumed by passing magmas. Having formed in spreading ridges, these eclogites have proven invaluable archives for the onset of plate tectonics, for the construction of cratons during subduction/collision, as probes of the convecting mantle from which their precursors formed, and as generators of heterogeneity upon recycling into Earth's convecting mantle. During subduction and until exhumation, interaction with fluids and melts (called metasomatism) can change the mineralogy, chemical composition, and physical properties of mantle eclogites, complicating their interpretation, but a comprehensive study of these effects is lacking so far. We investigated mantle eclogites from ancient continents (cratons) around the globe in order to define hallmarks of metasomatism by subduction‐related fluids and small‐volume ultramafic carbonated mantle melts. We find that the latter is pervasive and occurs predominantly at midlithospheric depths where seismic discontinuities are detected, typically causing diamond destruction and a reduction in density. This has consequences for their gravitational stability and for the interpretation of shearwave velocities in cratons.
Metasomatic evolution of coesite-bearing diamondiferous eclogite from the Udachnaya Kimberlite
(2020)
A coesite-bearing diamondiferous eclogite from the Udachnaya kimberlite (Daldyn field, Siberian craton) has been studied to trace its complex evolution recorded in rock-forming and minor mineral constituents. The eclogite sample is composed of rock-forming omphacite (60 vol%), garnet (35 vol%) and quartz/coesite (5 vol%) and contains intergranular euhedral zoned olivine crystals, up to 200 µm long, coexisting with phlogopite, orthopyroxene, clinopyroxene (secondary), K-feldspar, plagioclase, spinel, sodalite and djerfisherite. Garnet grains are zoned, with a relatively homogeneous core and a more magnesian overgrowth rim. The rim zones further differ from the core in having higher Zr/Y (6 times that in the cores), ascribed to interaction with, or precipitation from, a kimberlite-related melt. Judging by pressure-temperature estimates (~1200 °C; 6.2 GPa), the xenolith originated at depths of ~180–200 km at the base of the continental lithosphere. The spatial coexistence of olivine, orthopyroxene and coesite/quartz with K-Na-Cl minerals in the xenolith indicates that eclogite reacted with a deep-seated kimberlite melt. However, Fe-rich olivine, orthopyroxene and low-pressure minerals (sodalite and djerfisherite) likely result from metasomatic reaction at shallower depths during transport of the eclogite by the erupting kimberlite melt. Our results demonstrate that a mixed eclogitic-peridotitic paragenesis, reported previously from inclusions in diamond, can form by interaction of eclogite and a kimberlite-related melt.
Schon zu früheren Zeiten in der Erdgeschichte gab es Warmzeiten durch starke Treibhauseffekte, mit tropischen Temperaturen in weiten Teilen der Erde, hohem Meeresspiegel und massivem Artensterben. Das belegen Daten aus der Paläoklimatologie. Wenn man heutige Klimamodelle auf solche geologischen Warmzeiten anwendet, kann man sie testen und verbessern. So verhilft die Paläoklimatologie zu einem präziseren Blick in unsere Klimazukunft.
Diamond formation in the Earth has been extensively discussed in recent years on the basis of geochemical analysis of natural materials, high-pressure experimental studies, or theoretical aspects. Here, we demonstrate experimentally for the first time, the spontaneous crystallization of diamond from CH4-rich fluids at pressure, temperature and redox conditions approximating those of the deeper parts of the cratonic lithospheric mantle (5-7 GPa) without using diamond seed crystals or carbides. In these experiments the fluid phase is nearly pure methane, even though the oxygen fugacity was significantly above metal saturation. We propose several previously unidentified mechanisms that may promote diamond formation under such conditions and which may also have implications for the origin of sublithospheric diamonds. These include the hydroxylation of silicate minerals like olivine and pyroxene, H2 incorporation into these phases and the "etching" of graphite by H2 and CH4 and reprecipitation as diamond. This study also serves as a demonstration of our new high-pressure experimental technique for obtaining reduced fluids, which is not only relevant for diamond synthesis, but also for investigating the metasomatic origins of diamond in the upper mantle, which has further implications for the deep carbon cycle.
The turnover time of terrestrial ecosystem carbon is an emergent ecosystem property that quantifies the strength of land surface on the global carbon cycle–climate feedback. However, observation- and modeling-based estimates of carbon turnover and its response to climate are still characterized by large uncertainties. In this study, by assessing the apparent whole ecosystem carbon turnover times (τ) as the ratio between carbon stocks and fluxes, we provide an update of this ecosystem level diagnostic and its associated uncertainties in high spatial resolution (0.083∘) using multiple, state-of-the-art, observation-based datasets of soil organic carbon stock (Csoil), vegetation biomass (Cveg) and gross primary productivity (GPP). Using this new ensemble of data, we estimated the global median τ to be 43+7−7 yr (median+difference to percentile 75−difference to percentile 25) when the full soil is considered, in contrast to limiting it to 1 m depth. Only considering the top 1 m of soil carbon in circumpolar regions (assuming maximum active layer depth is up to 1 m) yields a global median τ of 37+3−6 yr, which is longer than the previous estimates of 23+7−4 yr (Carvalhais et al., 2014). We show that the difference is mostly attributed to changes in global Csoil estimates. Csoil accounts for approximately 84 % of the total uncertainty in global τ estimates; GPP also contributes significantly (15 %), whereas Cveg contributes only marginally (less than 1 %) to the total uncertainty. The high uncertainty in Csoil is reflected in the large range across state-of-the-art data products, in which full-depth Csoil spans between 3362 and 4792 PgC. The uncertainty is especially high in circumpolar regions with an uncertainty of 50 % and a low spatial correlation between the different datasets (0.2<r<0.5) when compared to other regions (0.6<r<0.8). These uncertainties cast a shadow on current global estimates of τ in circumpolar regions, for which further geographical representativeness and clarification on variations in Csoil with soil depth are needed. Different GPP estimates contribute significantly to the uncertainties of τ mainly in semiarid and arid regions, whereas Cveg causes the uncertainties of τ in the subtropics and tropics. In spite of the large uncertainties, our findings reveal that the latitudinal gradients of τ are consistent across different datasets and soil depths. The current results show a strong ensemble agreement on the negative correlation between τ and temperature along latitude that is stronger in temperate zones (30–60∘ N) than in the subtropical and tropical zones (30∘ S–30∘ N). Additionally, while the strength of the τ–precipitation correlation was dependent on the Csoil data source, the latitudinal gradients also agree among different ensemble members. Overall, and despite the large variation in τ, we identified robust features in the spatial patterns of τ that emerge beyond the differences stemming from the data-driven estimates of Csoil, Cveg and GPP. These robust patterns, and associated uncertainties, can be used to infer τ–climate relationships and for constraining contemporaneous behavior of Earth system models (ESMs), which could contribute to uncertainty reductions in future projections of the carbon cycle–climate feedback. The dataset of τ is openly available at https://doi.org/10.17871/bgitau.201911 (Fan et al., 2019).
The turnover time of terrestrial carbon (τ) controls the global carbon cycle – climate feedback and, yet, is poorly simulated by the current Earth System Models (ESMs). In this study, by assessing apparent carbon turnover time as the ratio between carbon stocks and fluxes, we provide a new, updated ensemble of diagnostic terrestrial carbon turnover times and associated uncertainties on a global scale using multiple, state-of-the-art, observation-based datasets of soil organic carbon stock (Csoil), vegetation biomass (Cveg) and gross primary productivity (GPP). Using this new ensemble, we estimated the global average τ to be 42$% &' years when the full soil depth is considered, longer than the previous estimates of 23$) &* years. Only considering the top 1 m (assuming maximum active layer depth is up to 1 meter) of soil carbon in circumpolar regions yields a global τ of 35$) &' years. Csoil in circumpolar regions account for two thirds of the total uncertainty in global τ estimates, whereas Csoil in non-circumpolar contributes merely 9.38%. GPP (2.25%) and Cveg (0.05%) contribute even less to the total uncertainty. Therefore, the high uncertainty in Csoil is the main factor behind the uncertainty in global τ, as reflected in the larger range of full-depth Csoil (3152-4372 PgC). The uncertainty is especially high in circumpolar regions with a behaviour of ESMs which could contribute to uncertainty reductions in future projections of the carbon cycle - climate feedback. The dataset of the terrestrial turnover time ensemble (DOI: 10.17871/bgitau.201911) is openly available from the data portal: https://doi.org/10.17871/bgitau.201911 (Fan et al., 2019) uncertainty of 50% and the spatial correlations among different datasets are also low compared to other regions. Overall, we argue that current global datasets do not support robust estimates of τ globally, for which we need clarification on variations of Csoil with soil depth and stronger estimates of Csoil in circumpolar regions. Despite the large variation in both magnitude and spatial patterns of τ, we identified robust features in the spatial patterns of τ that emerge regardless of soil depth and differences in data sources of Csoil, Cveg and GPP. Our findings show that the latitudinal gradients of τ are consistent across different datasets and soil depth. Furthermore, there is a strong consensus on the negative correlation between τ and temperature along latitude that is stronger in temperate zones (30ºN-60ºN) than in subtropical and tropical zones (30ºS30ºN). The identified robust patterns can be used to infer the response of τ to climate and for constraining contemporaneous behaviour of ESMs which could contribute to uncertainty reductions in future projections of the carbon cycle - climate feedback. The dataset of the terrestrial turnover time ensemble (DOI:10.17871/bgitau.201911) is openly available from the data portal: https://doi.org/10.17871/bgitau.201911 (Fan et al., 2019).
We apply seismic full waveform inversion to SH‐ and Love‐wave data for investigating the near‐surface lithology at an archaeological site. We evaluate the resolution of the applied full waveform inversion algorithm through ground truthing in the form of an excavation and sediment core studies. Thereby, we investigate the benefits of full waveform inversion in comparison with other established methods of near‐surface prospecting in terms of resolution capabilities and interpretation security. The study is performed in a presumed harbour area of the ancient Thracian city of Ainos. The exemplary target is the source of a linear magnetic anomaly oriented perpendicular to the coast, which was found in a previous magnetic gradiometry survey, suggesting a mole. The SH‐wave full waveform inversion recovered a subsurface SH‐wave velocity model with submeter resolution showing lateral and vertical velocity variation between 40 and 150 m/s. To tame the non‐linearity of the full waveform inversion, a sequential inversion of frequency bands has to be combined with time‐windowing in order to separate the Love wave from the reflected SH wavefield. We compare the full waveform inversion results with multichannel analysis of surface waves, standard seismic reflection imaging, electrical resistivity tomography and electromagnetic induction. It turns out that the respective depth sections are correlated to a certain degree with the full waveform inversion results. However, the structural resolution of the other geophysical methods is significantly lower than for the full waveform inversion. An exception is the reflection seismic imaging, which shows the same resolution as full waveform inversion but can only be interpreted together with the full waveform inversion–based velocity model. An archaeological excavation as well as coring data allows ground truthing and a direct understanding of the geophysical structures. The results show that the target was a sort of near‐surface trench of about 3–4 m width and 0.8 m to 1.0 m depth, filled with silty sediment, which differs from the layered surrounding in colour and composition. The ground truthing revealed that only SH‐wave full waveform inversion and seismic reflection imaging could image the trench and sediment structure with satisfying lateral and depth resolution. We emphasize that the velocity distribution from SH‐wave full waveform inversion agrees closely with the excavated subsurface structures, and that the discovered changes in seismic velocity correlate with changes in the sand content in the respective sediment facies sequences. The study demonstrated that SH‐wave full waveform inversion is capable to image structural and lithological changes in the near subsurface at scales as low as 0.5 m, thus providing the high resolution needed for archaeological and geoarchaeological prospection.
Remote sensing data are essential for monitoring the Earth’s surface waters, especially since the amount of publicly available in-situ data is declining. Satellite altimetry provides valuable information on the water levels and variations of lakes, reservoirs and rivers. In combination with satellite imagery, the derived time series allow the monitoring of lake storage changes and river discharge. However, satellite altimetry is limited in terms of its spatial resolution due to its measurement geometry, only providing information in the nadir direction beneath the satellite’s orbit. In a case study in the Mississippi River Basin (MRB), this study investigates the potential and limitations of past and current satellite missions for the monitoring of basin-wide storage changes. For that purpose, an automated target detection is developed and the extracted lake surfaces are merged with the satellites’ tracks. This reveals that the current altimeter configuration misses about 80% of all lakes larger than 0.1 km2 in the MRB and 20% of lakes larger than 10 km2, corresponding to 30% and 7% of the total water area, respectively. Past altimetry configurations perform even more poorly. From the larger water bodies represented by a global hydrology model, at least 91% of targets and 98% of storage changes are captured by the current altimeter configuration. This will improve significantly with the launch of the planned Surface Water and Ocean Topography (SWOT) mission.
Forschung Frankfurt : das Wissenschaftsmagazin der Goethe-Universität. 2020, Nr. 2 ; Klimakrise
(2020)
Aus der Redaktion
(2020)
This article presents the findings from systematically reviewing 26 empirical research studies published from 2005 to 2014 on the use of GIS for learning and teaching. By employing methods of narrative synthesis and qualitative content analysis, the study gives evidence about the state of knowledge of competence‐based GIS education. The results explain what factors and variables effect GIS learning in terms of technology use, major subject contents, learning contexts, and didactic and pedagogical aspects. They also show what facets of knowledge, process skills, and affect the research literature has investigated. The analysis of the type and quality of the methods used indicates that current GIS education research is a heterogeneous field that needs a systematic research framework for future efforts, according to empirical education research.
Tropical cyclones (TC) represent a substantial threat to life and property for Caribbean and adjacent populations. The prospective increase of TC magnitudes, expressed in the 15th chapter of the IPCC AR5 report, entails a rising probability of ecological and social disasters, which were tragically exemplified by several severe Caribbean TC strikes during the past 20 years. Modern IPCC-grade climate models, however, still lack the required spatial and temporal resolution to accurately consider the underlying boundary conditions that modulate long-time TC patterns beyond the Instrumental Era. It is thus necessary to provide a synoptic mechanistic understanding regarding the origin of such long-time patterns, in order to predict reliable changes of TC magnitude and frequency under future climate scenarios. Caribbean TC records are still rare and often lack the necessary continuity and resolution to overcome these limitations. Here, we report on an annually-resolved sedimentary archive from the bottom of the Great Blue Hole (Lighthouse Reef, Belize). The TC record encompasses 1885 years and extends all existing site-specific TC archives both in terms of resolution and duration. We identified a likely connection between long-term TC patterns and climate phenomena responses to Common Era climate variations and offer a conceptual and comparative view considering several involved tropospheric and oceanographic control mechanisms such as the El-Niño-Southern-Oscillation, the North Atlantic Oscillation and the Atlantic Multidecadal Oscillation. These basin-scaled climate modes exercise internal control on TC activity by modulating the thermodynamic environment (sea-surface temperature and vertical wind shear stress dynamics) for enhanced/suppressed TC formation both on millennial (primary) and multi-decadal (secondary) time scales. We interpret the beginning of the Medieval Warm Period (MWP) as an important time interval of the Common Era record and suspect that the southward migration of the intertropical convergence zone (ITCZ) caused, in combination with extensive hydro-climate changes, a shift in the tropical Atlantic TC regime. The TC activity in the south-western Caribbean changed in general from a stable and less active stage (100–900 CE) to a more active and variable state (1,100 CE-modern).
Mit dem Klimawandel könnte das Grundwasser künftig weltweit in vielen Regionen knapp werden. Das zeigt eine große internationale Studie, die vom Doktoranden Robert Reinecke und der Geographin Prof. Petra Döll initiiert wurde. Vielerorts könnte das zu Wassermangel führen oder einen bestehenden Wassermangel verschärfen. Gleichzeitig werden andere Regionen unter steigenden Grundwasserständen zu leiden haben. Ackerflächen könnten durch Vernässung verloren gehen.
The endemic argan woodlands cover large parts of South Morocco and create a characteristic landscape with areas of sparsely vegetated and bare soil surfaces between single trees. This unique ecosystem has been under extensive agrosilvopastoral management for centuries and is now at risk of degradation caused by overgrazing and increasing scarcity and variability of rainfall.
To investigate susceptibility to wind erosion, we conducted an experimental–empirical study including wind tunnel tests and a drone‐generated digital elevation model and quantified wind‐erodible material on five different associated surface types by means of sediment catchers. The highest emission flux was measured on freshly ploughed surfaces (1875 g m–2 h–1), while older ploughed areas with a re‐established crust produced a much lower emission flux (795 g m–2 h–1). Extensive tillage may have been a sustainable practice for generations, but increasing drought and uncertainty of rainfall now lead to an acute risk of severe soil erosion and dust production. The typical crusted surfaces characterized by residual rock fragment accumulation and wash processes produced the second highest emission flux (1,354 g m–2 h–1). Material collected from tree‐shaded areas (933 g m–2 h–1) was revealed to be a considerable source of organic material, possibly affecting substrate conditions positively on a larger regional scale. The lowest flux was measured on rock fragment‐covered surfaces (301 g m–2 h–1).
The data show that open argan woodland may be a considerable source for wind erosion and dust production, depending on surface characteristics strongly related to management. An adapted management must include the conservation of argan trees to offer a promising approach to prevent severe wind erosion and dust production and mitigate possible impacts of land‐use change and climate change related shifts in wind and rainfall patterns.
In partially molten regions inside the earth melt buoyancy may trigger upwelling of both solid and fluid phases, i.e. diapirism. If the melt is allowed to move separately with respect to the matrix, melt perturbations may evolve into solitary porosity waves. While diapirs may form on a wide range of scales, porosity waves are restricted to sizes of a few times the compaction length. Thus, the size of a partially molten perturbation controls whether a diapir or a porosity wave will emerge. We study the transition from diapiric rise to solitary porosity waves by solving the two-phase flow equations of conservation of mass and momentum in 2D with porosity dependent matrix viscosity. We systematically vary the initial size of a porosity perturbation from 1 to 100 times the compaction length. If the perturbation is much larger than a regular solitary wave, its Stokes velocity is large and therefore faster than the segregating melt. Consequently, the fluid is not able to form a porosity wave and a diapir emerges. For small perturbations solitary waves emerge, either with a positive or negative vertical matrix velocity inside. In between the diapir and solitary wave regimes we observe a third regime of solitary wave induced focusing of melt. In these cases, diapirism is dominant but the fluid is still fast enough to locally build up small solitary waves which rise slightly faster than the diapir and form finger like structures at the front of the diapir. In our numerical simulations the width of these fingers is controlled by the compaction length or the grid size, whichever is larger. In cases where the compaction length becomes similar to or smaller than the grid size the finger-like leading solitary porosity waves are no more properly resolved, and too big and too fast waves may be the result. Therefore, one should be careful in large scale two-phase flow modelling with melt focusing especially when compaction length and grid size are of similar order.
Diurnal valley winds frequently form over complex topography, particularly under fair weather conditions, and have a significant impact on the local weather and climate. Since diurnal valley winds result from complex and multi-scale interactions, their representation in numerical weather prediction models is challenging. Better understanding of these local winds based on observations is crucial to improve the accuracy of the forecasts. This study investigates the diurnal evolution of the three-dimensional mean wind structure in a deep Alpine valley, the Rhone valley at Sion, using data from a radar wind profiler and a surface weather station operated continuously from 1 September 2016 to 17 July 2017. In particular, the wind profiler data was analyzed for a subset of days on which fair weather conditions allowed for the full development of thermally driven winds. A pronounced diurnal cycle of the wind speed, as well as a reversal of the wind direction twice per day is documented for altitudes up to about 2 km above ground level (AGL) in the warm season and less than 1 km AGL in winter. The diurnal pattern undergoes significant changes during the course of the year. Particularly during the warm-weather months of May through to September, a low-level wind maximum occurs, where mean maximum up-valley velocities of 8–10 m s−1 are found between 15–16 UTC at altitudes around 200 m AGL. In addition, during nighttime, a down-valley jet with maximum wind speeds of 4–8 m s−1 around 1 km AGL is found. A case study of a three-day period in September 2016 illustrates the occurrence of an elevated layer of cross-valley flow around 1–1.5 km AGL.
The most frequently used boundary-layer turbulence parameterization in numerical weather prediction (NWP) models are turbulence kinetic energy (TKE) based-based schemes. However, these parameterizations suffer from a potential weakness, namely the strong dependence on an ad-hoc quantity, the so-called turbulence length scale. The physical interpretation of the turbulence length scale is difficult and hence it cannot be directly related to measurements or large eddy simulation (LES) data. Consequently, formulations for the turbulence length scale in basically all TKE schemes are based on simplified assumptions and are model-dependent. A good reference for the independent evaluation of the turbulence length scale expression for NWP modeling is missing. Here we propose a new turbulence length scale diagnostic which can be used in the gray zone of turbulence without modifying the underlying TKE turbulence scheme. The new diagnostic is based on the TKE budget: The core idea is to encapsulate the sum of the molecular dissipation and the cross-scale TKE transfer into an effective dissipation, and associate it with the new turbulence length scale. This effective dissipation can then be calculated as a residuum in the TKE budget equation (for horizontal sub-domains of different sizes) using LES data. Estimation of the scale dependence of the diagnosed turbulence length scale using this novel method is presented for several idealized cases.
The weather of the atmospheric boundary layer significantly affects our life on Earth. Thus, a realistic modelling of the atmospheric boundary layer is crucial. Hereby, the processes of the atmospheric boundary layer depend on an accurate representation of the land-atmosphere coupling in the model. In this context the land surface temperature (LST) plays an important role. In this thesis, it is examined if the assimilation of LST can lead to improved estimates of the boundary layer and its processes.
To properly assimilate the LST retrievals, a suitable model equivalent in the weather prediction model is necessary. In the weather forecast model of the German Weather Service used here, the LST is modelled without a vegetation temperature. To compensate for this deficit, two different vegetation parameterizations were investigated and the better one, a conductivity scheme, was implemented. In order to make optimal use of the influence of the assimilation of the LST observation on the model system, it is useful to pass on the information of the observation to land and atmosphere already in the assimilation step. For that reason, a fully coupled land-atmosphere prediction model was used. Therefore, the existing control vector of the assimilation system, a local ensemble transform Kalman filter, was extended by the soil temperature and moisture. In two-day case studies in March and August 2017, different configurations of the augmented assimilation system were evaluated based on observing system simulation experiments (OSSE).
LST was assimilated hourly over two days in the weakly and strongly coupled assimilation system. In addition, every six hours a free 24-hour forecast was simulated. The experiments were validated with the simulated truth (a high-resolution model run) and compared against an experiment without assimilation. It was shown that the prediction of the boundary layer temperature, especially during the day, and the prediction of the soil temperature, during the whole day and night, could be improved.
The best impact of LST assimilation was achieved with the fully coupled system. The humidity variables of the model benefited only partially from the LST assimilation. For this reason, covariances in the model ensemble were investigated in more detail. To check their compatibility with the high-resolution model run the ensemble consistency score was introduced. It was found that the covariances between the LST and the temperatures of the high-resolution model run were better represented in the ensemble than those between the LST and the humidity variables.
Drought is understood as both a lack of water (i.e., a deficit as compared to some requirement) and an anomaly in the condition of one or more components of the hydrological cycle. Most drought indices, however, only consider the anomaly aspect, i.e., how unusual the condition is. In this paper, we present two drought hazard indices that reflect both the deficit and anomaly aspects. The soil moisture deficit anomaly index, SMDAI, is based on the drought severity index, DSI, but is computed in a more straightforward way that does not require the definition of a mapping function. We propose a new indicator of drought hazard for water supply from rivers, the streamflow deficit anomaly index, QDAI, which takes into account the surface water demand of humans and freshwater biota. Both indices are computed and analyzed at the global scale, with a spatial resolution of roughly 50 km, for the period 1981-2010, using monthly time series of variables computed by the global water resources and the model WaterGAP2.2d. We found that the SMDAI and QDAI values are broadly similar to values of purely anomaly-based indices. However, the deficit anomaly indices provide more differentiated, spatial and temporal patterns that help to distinguish the degree of the actual drought hazard to vegetation health or the water supply. QDAI can be made relevant for stakeholders with different perceptions about the importance of ecosystem protection, by adapting the approach for computing the amount of water that is required to remain in the river for the well being of the river ecosystem. Both deficit anomaly indices are well suited for inclusion in local or global drought risk studies.
This article presents the findings from systematically reviewing 26 empirical research studies published from 2005 to 2014 on the use of GIS for learning and teaching. By employing methods of narrative synthesis and qualitative content analysis, the study gives evidence about the state of knowledge of competence-based GIS education. The results explain what factors and variables effect GIS learning in terms of technology use, major subject contents, learning contexts, and didactic and pedagogical aspects. They also show what facets of knowledge, process skills, and affect the research literature has investigated. The analysis of the type and quality of the methods used indicates that current GIS education research is a heterogeneous field that needs a systematic research framework for future efforts, according to empirical education research.
Extensive black shale deposits formed in the Early Cretaceous South Atlantic, supporting the notion that this emerging ocean basin was a globally important site of organic carbon burial. The magnitude of organic carbon burial in marine basins is known to be controlled by various tectonic, oceanographic, hydrological, and climatic processes acting on different temporal and spatial scales, the nature and relative importance of which are poorly understood for the young South Atlantic. Here we present new bulk and molecular geochemical data from an Aptian–Albian sediment record recovered from the deep Cape Basin at Deep Sea Drilling Project (DSDP) Site 361, which we combine with general circulation model results to identify driving mechanisms of organic carbon burial. A multi-million year decrease (i.e. Early Aptian–Albian) in organic carbon burial, reflected in a lithological succession of black shale, gray shale, and red beds, was caused by increasing bottom water oxygenation due to abating tectonic restriction via South Atlantic-Southern Ocean gateways. These results emphasize basin evolution and ocean gateway development as a decisive primary control on enhanced organic carbon preservation in the Cape Basin at geological time scales (>1 Myr). The Early Aptian black shale sequence comprises alternations of shales with high (>5%) and relatively low (~3%) organic carbon content of marine sources, the former being deposited during the global Oceanic Anoxic Event (OAE) 1a, as well as during repetitive events before and after OAE 1a. In all cases, these short-term events of enhanced organic carbon burial coincided with strong influxes of sediments derived from the proximal African continent, indicating closely coupled climate–land–ocean interactions. Supported by our model results, we propose that fluctuations in weathering-derived nutrient input from the southern African continent, linked to fluctuations in pCO2 and/or orbitally driven humidity/aridity, were the underlying drivers of short-term organic carbon burial in the deep Cape Basin. These results suggest that deep marine environments of emerging ocean basins responded sensitively and directly to short term fluctuations in riverine nutrient fluxes. We explain this relationship by the lack of wide and mature continental shelf seas that could have acted as a barrier or filter for nutrient transfer from the continent into the deep ocean.
Analysing the composition of ambient ultrafine particles (UFP) is a challenging task due to the low mass and chemical complexity of small particles, yet it is a prerequisite for the identification ofparticle sources and the assessment of potential health risks. Here, we show the molecular characterization of UFP, based on cascade impactor (Nano-MOUDI) 10samples that were collected at an air quality monitoring station nearby one of Europe`s largest airports in Frankfurt, Germany. At this station, particle-size-distribution measurements show enhanced number concentration of particles smaller than 50nm during airport operating hours. We sampled the lower UFP fraction (0.010-0.018 μm; 0.018-0.032 μm; 0.032-0.056 μm) when the air masses arrived from the airport. We developed an optimized filter extraction procedure, used ultra-high performance liquid chromatography (UHPLC) for compound separation, and a heated electrospray ionization (HESI) source with an 15Orbitrap high-resolution mass spectrometer (HRMS) as a detector for organic compounds. A non-target screening detected ~200 organic compounds in the UFP fraction with sample-to-blank ratios larger than five. We identified the largest signals as homologous series of pentaerythritol esters (PEE) and trimethylolpropane esters (TMPE), which are base stocks of aircraft lubrication oils. We unambiguously attribute the majority of detected compounds to jet engine lubrication oils by matching retention times, high-resolution/accurate mass (HR/AM) measurements, and comparing MS/MS fragmentation patterns between both ambient samples and commercially available jet oils. For each UFP stage, we created molecular fingerprints to visualize the complex chemical composition ofthe organic fraction and their average carbon oxidation state. These graphs underline the presence of the homologous series of PEE and TMPE, and the appearance of jet oil additives (e.g. tricresyl phosphate (TCP)). Targeted screening on TCP confirmed the absence of the harmful tri-orthoisomer, while we identified a thermal transformation product of TMPE-based lubrication oil (trimethylolpropane phosphate (TMP-P)). Even though a quantitative determination of the identified compounds is limited, the presented method enables the qualitative detection of molecular markers for jet engine lubricants in UFP and thus strongly improves the source apportionment of UFP near airports.
High-pressure single-crystal to 20 GPa and powder diffraction measurements to 50 GPa, show that the structure of Pb2SnO4 strongly distorts on compression with an elongation of one axis. A structural phase transition occurs between 10 GPa and 12 GPa, with a change of space group from Pbam to Pnam. The resistivity decreases by more than six orders of magnitude when pressure is increased from ambient conditions to 50 GPa. This insulator-to-semiconductor transition is accompanied by a reversible appearance change from transparent to opaque. Density functional theory-based calculations show that at ambient conditions the channels in the structure host the stereochemically-active Pb 6s2 lone electron pairs. On compression the lone electron pairs form bonds between Pb2+ ions. Also provided is an assignment of irreducible representations to the experimentally observed Raman bands.
Climate controls the broad-scale distribution of vegetation and change in climate will alter the vegetation distribution, biome boundaries, biodiversity, phenology and supply of ecosystem services. A better understanding of the consequences of climate change is required, particularly in under-investigated regions such as tropical Asia, i.e., South and South-east Asia, which is a host to 7 of the 36 global biodiversity hotspots. Conservation strategies would also require an in-depth understanding of the response of vegetation to climate change. Therefore, the main objective of this thesis was to investigate the impact of climate change and rising CO2 vegetation in tropical Asia. Dynamic global vegetation model (DGVMs) are the well-known tools to investigate vegetation-climate interactions and climate change impacts on ecosystems. In this thesis, I used a complex trait-based DGVM called adaptive dynamic vegetation model version 2 (aDGVM2).
In Chapter 1, I presented a brief background of the phytogeography and discussed the exiting knowledge gap on vegetation-climate interactions in the region. One major disadvantage for available DGVMs studies for the tropical Asia is that most of them have used fixed plant functional types (PFTs) and do not explicitly represent the distinct varieties of vegetation type of the region such as Asian savannas. In Chapter 2, I discussed at great length to improve DGVMs for South Asia and discussed ways to include them in the model for better representation of region vegetation-climate interaction.
I upgraded the current version of aDGVM2 and added a new vegetation type i.e., C3 grasses, and modified the sub-module to simulate photosynthesis for each individual plants to aDGVM2. In chapter 3, I used this updated version of aDGVM2 to simulate the current and future vegetation distribution in South Asia under RCP4.5 and RCP8.5 (RCP: representative concentration pathway). The model predicted an increase in biomass, canopy cover, and tree height under the presence of CO2 fertilization, which triggered transitions towards tree-dominated biomes by the end of the 21st century under both RCPs. I found that vegetation along the Western Ghats and the Himalayas are more susceptible to change due to climate change and open biomes such as grassland and savanna are prone to woody encroachment.
In Chapter 4, the study domain was extended to include South-east Asia to verify if the model configuration used in Chapter 3 can also simulate vegetation patterns in tropical Asia. The aDGVM2 simulations showed a robust trend of increasing vegetation biomass and transitions from small deciduous vegetation to taller evergreen vegetation across most of tropical Asia. Shifts in plant phenology also affect ecosystem carbon cycles and ecosystem feedback to climate, yet the quantification of such impacts remains challenging. The study showed increased biomass due to CO2 fertilization, indicates that the region can remain a carbon sink given there is no other resource limitation. However, nutrient limitations on CO2 fertilization effects were not included in the study, and carbon sink potential has to be seen with caution.
In Chapter 5, I focused on Asian savannas, which have been mismanaged since the colonial era due to misinterpretation as a degraded forest. I proposed a biome classification scheme to distinguish between degraded forest or woodland and savanna based on the abundance of grass biomass and canopy cover. I found that considering vegetation systems as woodland or degraded forest could easily be mistaken as a potential for forest restoration within a tree-centric perspective. This would put approximately 35% to 40% of a unique savanna biome at risk. Although projected woody encroachments may imply a transition toward the forest that benefits climate mitigation. This raises potential conflicts of interest between biodiversity conservation in open ecosystems, i.e., savanna and active afforestation, to enhance carbon sequestration. Proper management strategies should be taken into account to maintain a balance for both objective
In conclusion, the model predicted that vegetation in South and South-East Asia would significantly shift towards tree-dominated biomes due to CO2-induced fertilization of C3-photosynthesis. The simulation under fixed CO2 and rising CO2 scenarios clearly showed that rising level of atmospheric CO2 is responsible for most of the predicted change in biome properties. This study is an important step towards understanding ecosystems of South and Southeast Asia, specifically savannas. The aDGVM2 can serve as tools to inform decision making for climate adaptation and mitigation for savanna. The thesis, thus contributes to our ability to improve conservation strategies to mitigate the consequences of climate change.
The Ceboruco is a 2280 m high stratovolcano located in Nayarit State, Mexico. Despite its last eruption which occurred in 1870, it is the most active volcano in the area, showing volcanicearthquake activity together with ongoing vapor emissions. The magnetotelluric survey was carried out in November 2016. It was part of a geothermal project (CeMIEGeo-P24) and focused on the determination of the electrical conductivity distribution in the subsurface of the volcano.
The Magnetotelluric Apparent Resistivity Tensor, as introduced by Brown (2016), can be decomposed into an amplitude and a phase tensor. The fundamental physics behind those new tensors were presented in Hering et al. (2019), using canonical models in 1-D (isotropic and anisotropic) and 2-D resistivity environments. Here, the tensors are introduced for a high-quality data set, where their interpretational benefits become very obvious. Additionally, results from an isotropic 3-D inversion are presented and compared to an alternative 3-D anisotropic forward model.
A list of authors and their affiliations appears at the end of the paper New-particle formation is a major contributor to urban smog, but how it occurs in cities is often puzzling. If the growth rates of urban particles are similar to those found in cleaner environments (1–10 nanometres per hour), then existing understanding suggests that new urban particles should be rapidly scavenged by the high concentration of pre-existing particles. Here we show, through experiments performed under atmospheric conditions in the CLOUD chamber at CERN, that below about +5 degrees Celsius, nitric acid and ammonia vapours can condense onto freshly nucleated particles as small as a few nanometres in diameter. Moreover, when it is cold enough (below −15 degrees Celsius), nitric acid and ammonia can nucleate directly through an acid–base stabilization mechanism to form ammonium nitrate particles. Given that these vapours are often one thousand times more abundant than sulfuric acid, the resulting particle growth rates can be extremely high, reaching well above 100 nanometres per hour. However, these high growth rates require the gas-particle ammonium nitrate system to be out of equilibrium in order to sustain gas-phase supersaturations. In view of the strong temperature dependence that we measure for the gas-phase supersaturations, we expect such transient conditions to occur in inhomogeneous urban settings, especially in wintertime, driven by vertical mixing and by strong local sources such as traffic. Even though rapid growth from nitric acid and ammonia condensation may last for only a few minutes, it is nonetheless fast enough to shepherd freshly nucleated particles through the smallest size range where they are most vulnerable to scavenging loss, thus greatly increasing their survival probability. We also expect nitric acid and ammonia nucleation and rapid growth to be important in the relatively clean and cold upper free troposphere, where ammonia can be convected from the continental boundary layer and nitric acid is abundant from electrical storms.
Background: Point of care devices for performing targeted coagulation substitution in bleeding patients have become increasingly important in recent years. New on the market is the Quantra® from HemoSonics (LC, Charlottesville, VA, US). It uses sonorheometry, a sonic estimation of elasticity via resonance (SEER), a novel ultrasound-based technology that measures viscoelastic properties of whole blood. Several studies have already shown the comparability with devices already established on the market such as the ROTEM® (TEM International GmbH, Munich, Germany).
Objective: In contrast to existing studies, the planned study will be the first prospective interventional study using the new Quantra® system in a cardiac surgical patient cohort. The aim is to investigate the non-inferiority between an already existing coagulation algorithm, based on ROTEM®/Multiplate®, and a new algorithm based on the Quantra®, for the treatment of coagulopathic cardiac surgical patients.
Methods: The study is divided into two phases. In an initial observation phase, whole blood samples of 20 patients will be analyzed using both ROTEM®/Multiplate® and Quantra® obtained at three defined points of time (prior to surgery, after completion of cardiopulmonary bypass, on arrival in the intensive care unit). The obtained threshold values will be used to create an algorithm for hemotherapy. In a second intervention phase, the new algorithm will be tested against an algorithm used routineously for years at our department for non-inferiority.
Results: The main objective of the examination is the cumulative loss of blood within 24 hours after surgery. Statistical calculations based on literature and in-house data suggest that the new algorithm is not inferior if the difference in cumulative blood loss is < 150ml/24 h.
Conclusions: Because of the comparability of the Quantra® sonorheometry system with ROTEM® rotational thromboelastometric measurement methods, the existing hemotherapy treatment algorithm can be adapted to the Quantra device with a proof of non-inferiority. Clinical Trial: International Registered Report Identifier (IRRID): clinicaltrials.gov: NCT03902275
Groundwater is the largest source of accessible freshwater with its dynamics having significantly changed due to human withdrawals, and being projected to continue to as a result of climate change. The pumping of groundwater has led to lowered water tables, decreased base flow, and depletion.
Global hydrological models (GHMs) are used to simulate the global freshwater cycle, assessing impacts of changes in climate and human freshwater use. Currently, groundwater is commonly represented by a bucket-like linear storage component in these models. Bucket models, however, cannot provide information on the location of the groundwater table. Due to this limitation, they can only simulate groundwater discharge to surface water bodies but not recharge from surface water to groundwater and calculate no lateral and vertical groundwater flow whatsoever among grid cells. For instance this may lead to an underestimation of groundwater resources in semiarid areas, where groundwater is often replenished by surface water. In order to overcome these limitations it is necessary to replace the linear groundwater model in GHMs with a hydraulic head gradient-based groundwater flow model
This thesis presents the newly developed global groundwater model G3M and its coupling to the GHM WaterGAP spanning over 70,000 lines of newly developed code. Development and validation of the modeling software are discussed along with numerical challenges. Based on the newly developed software, a global natural equilibrium groundwater model is presented showing better agreements with observations than previous models. Groundwater discharge to rivers is found to be the most dominant flow component globally, compared to flows to other surface water bodies and lateral flows. Furthermore, first global maps of the distribution of gaining and losing surface water bodies are displayed.
For the purpose of determining the uncertainty in model outcomes a sensitivity study is conducted with an innovative approach through applying a global sensitivity analysis for a computationally complex model. First global maps of spatially distributed parameter sensitivities are presented. The results at hand indicate that globally simulated hydraulic heads are equally sensitive to hydraulic conductivity, groundwater recharge and surface water body elevation, even though parameter sensitivities do vary regionally.
A high resolution model of New Zealand is developed to further understand the involved uncertainties connected to the spatial resolution of the global model. This thesis finds that a new understanding is necessary how these models can be evaluated and that a simple increase in spatial resolution is not improving the model performance when compared to observations.
Alongside the assessment of the natural equilibrium, the concept of a fully coupled transient model as integrated storage component replacing the former model in the hydrological model WaterGAP is discussed. First results reveal that the model shows reasonable response to seasonal variability although it contains persistent head trends leading to global overestimates of water table depth due to an incomplete coupling. Nonetheless, WaterGAP-G3M is already able to show plausible long term storage trends for areas that are known to be affected by groundwater depletion. In comparison with two established regional models in the Central Valley the coupled model shows a highly promising simulation of storage declines.
Global‐scale gradient‐based groundwater models are a new endeavor for hydrologists who wish to improve global hydrological models (GHMs). In particular, the integration of such groundwater models into GHMs improves the simulation of water flows between surface water and groundwater and of capillary rise and thus evapotranspiration. Currently, these models are not able to simulate water table depth adequately over the entire globe. Unsatisfactory model performance compared to well observations suggests that a higher spatial resolution is required to better represent the high spatial variability of land surface and groundwater elevations. In this study, we use New Zealand as a testbed and analyze the impacts of spatial resolution on the results of global groundwater models. Steady‐state hydraulic heads simulated by two versions of the global groundwater model G3M, at spatial resolutions of 5 arc‐minutes (9 km) and 30 arc‐seconds (900 m), are compared with observations from the Canterbury region. The output of three other groundwater models with different spatial resolutions is analyzed as well. Considering the spatial distribution of residuals, general patterns of unsatisfactory model performance remain at the higher resolutions, suggesting that an increase in model resolution alone does not fix problems such as the systematic overestimation of hydraulic head. We conclude that (1) a new understanding of how low‐resolution global groundwater models can be evaluated is required, and (2) merely increasing the spatial resolution of global‐scale groundwater models will not improve the simulation of the global freshwater system.
Aquatic ecosystems are globally contaminated with microplastics (MP). However, comparative data on MP levels in freshwater systems is still scarce. Therefore, the aim of this study is to quantify MP abundance in water and sediment of the German river Elbe using visual, spectroscopic (Fourier-transform infrared spectroscopy) and thermo analytical (pyrolysis gas chromatography mass spectrometry) methods. Samples from eleven German sites along the German part of the Elbe were collected, both in the water and sediment phase, in order to better understand MP sinks and transport mechanisms. MP concentrations differed between the water and sediment phase. Sediment concentrations (mean: 3,350,000 particles m−3, 125–5000 μm MP) were in average 600,000-fold higher than water concentrations (mean: 5.57 particles m−3, 150–5000 μm MP). The abundance varied between the sampling sites: In sediments, the abundance decreased in the course of the river while in water samples no such clear trend was observed. This may be explained by a barrage retaining sediments and limiting tidal influence in the upstream parts of the river. Particle shape differed site-specifically with one site having exceptionally high quantities of spheres, most probably due to industrial emissions of PS-DVB resin beads. Suspended MP consisted predominantly of polyethylene and polypropylene whereas sediments contained a higher diversity of polymer types. Determined MP concentrations correspond well to previous results from other European rivers. In a global context, MP levels in the Elbe relate to the lower (water) to middle section (sediment) of the global range of MP concentrations determined for rivers worldwide. This highlights that elevated MP levels are not only found in single countries or continents, but that MP pollution is an issue of global concern.
During the first two days of August 2016 a seismic crisis occurred on Brava, Cabo Verde, which – according to observations based on a local seismic network – was characterized by more than a thousand volcano-seismic signals. Brava is considered an active volcanic island, although it has not experienced any historic eruptions. Seismicity significantly exceeded the usual level during the crisis. We report on results based on data from a temporary seismic-array deployment on the neighbouring island of Fogo at a distance of about 35 km. The array was in operation from October 2015 to December 2016 and recorded a total of 1343 earthquakes in the region of Fogo and Brava; 355 thereof were localized. On 1 and 2 August we observed 54 earthquakes, 25 of which could be located beneath Brava. We further evaluate the observations with regards to possible precursors to the crisis and its continuation. Our analysis shows a significant variation in seismicity around Brava, but no distinct precursory pattern. However, the observations suggest that similar earthquake swarms commonly occur close to Brava. The results further confirm the advantages of seismic arrays as tools for the remote monitoring of regions with limited station coverage or access.
During the first two days of August 2016 a seismic crisis occurred on Brava, Cape Verde, which – according to observations based on a local seismic network – was characterized by more than thousand volcano–seismic signals. Brava is considered an active volcanic island, although it has not experienced any historic eruptions. Seismicity significantly exceeded the usual level during the crisis. We report on results based on data from a temporary seismic–array deployment on the neighbouring island of Fogo at a distance of about 35 km. The array was in operation from October 2015 to December 2016 and recorded a total of 1343 earthquakes, 355 thereof were localized. On 1 and 2 August we observed 54 earthquakes, 25 of which could be located beneath Brava. We further evaluate the observations with regards to possible precursors to the crisis and its continuation. Our analysis shows a migration of seismicity around Brava, but no distinct precursory pattern. However, the observations suggest that similar earthquake swarms commonly occur close to Brava. The results further confirm the advantages of seismic arrays as tools for the remote monitoring of regions with limited station coverage or access.
One of the theoretical tensions that has arisen from Anthropocene studies is what Dipesh Chakrabarty has called the 'two figures of the human', and the question of which of these two figures of the human inheres in the concept of the Anthropocene more. On the one hand, the Human is conceived as the universal reasoning subject upon whom political rights and equality are based, and on the other hand, humankind is the collection of all individuals of our species, with all of the inequalities, differences, and variability inherent in any species category. This chapter takes up Deborah Coen's argument that Chakrabarty's claim of the 'incommensurability' of these two figures of the human ignores the way both were constructed within debates over how to relate local geophysical specificities to theoretical generalities. This chapter examines two cases in the history of science. The first is Martin Rudwick's historical exploration of how geologists slowly gained the ability to use fossils and highly local stratigraphic surveys to reconstruct the history of the Earth in deep time, rather than resort to speculative cosmological theory. The second is Coen's own history of imperial, Austrian climate science, a case where early nineteenth-century assumptions about the capriciousness of the weather gave way to theories of climate informed by thermodynamics and large-scale data collection.
Increases in water demand often result in unsustainable water use, leaving insufficient amounts of water for the environment. Therefore, water-saving strategies have been introduced to the environmental policy agenda in many (semi)-arid regions. As many such interventions failed to reach their objectives, a comprehensive tool is needed to assess them. We introduced a constructive framework to assess the proposed strategies by estimating five key components of the water balance in an area: (1) Demand; (2) Availability; (3) Withdrawal; (4) Depletion and (5) Outflow. The framework was applied to assess the Urmia Lake Restoration Program (ULRP) which aimed to increase the basin outflow to the lake to reach 3.1 × 109 m3 yr−1. Results suggested that ULRP could help to increase the Outflow by up to 57%. However, successful implementation of the ULRP was foreseen to be impeded because of three main reasons: (i) decreasing return flows; (ii) increased Depletion; (iii) the impact of climate change. Decreasing return flows and increasing Depletion were expected due to the introduction of technologies that increase irrigation efficiency, while climate change could decrease future water availability by an estimated 3–15%. We suggest that to reach the intervention target, strategies need to focus on reducing water depletion rather than water withdrawals. The framework can be used to comprehensively assess water-saving strategies, particularly in water-stressed basins.
Strain localization in the lithosphere and the formation, evolution, and maintenance of resulting plate boundaries play a crucial role in plate tectonics and thermo‐chemical mantle convection. Previously activated lithospheric deformation zones often appear to maintain a “memory” of weakening, leading to tectonic inheritance within plate reorganizations including the Wilson cycle. Different mechanisms have been proposed to explain such strain localization, but it remains unclear which operates on what spatio‐temporal scales, and how to best incorporate them in large‐scale mantle convection models. Here, we analyze two candidates, (1), grain‐size sensitive rheology and, (2), damage‐style parameterizations of yield, stress which are sometimes used to approximate the former. Grain‐size reduction due to dynamic recrystallization can drive localization in the ductile domain, and grain growth provides a time‐dependent rheological hardening component potentially enabling the preservation of rheological heterogeneities. We compare the dynamic weakening and hardening effects as well as the timescales of strength evolution for a composite rheology including grain‐size dynamics with a pseudo‐plastic rheology including damage‐ (or “strain”‐) dependent weakening. We explore the implications of different proposed grain‐size evolution laws, and test to which extent strain‐dependent rheologies can mimic the weakening and hardening effects of the more complex micro‐physical behavior. Such an analysis helps to better understand the parallels and differences between various strain‐localization modeling approaches used in different tectonics and geodynamics communities. More importantly, our results contribute to efforts to identify the key ingredients of strain‐localization and damage hysteresis within plate tectonics and how to represent those in planetary‐scale modeling.
Ice-nucleating particle concentrations of the past: insights from a 600-year-old Greenland ice core
(2020)
Ice-nucleating particles (INPs) affect the microphysics in cloud and precipitation processes. Hence, they modulate the radiative properties of clouds. However, atmospheric INP concentrations of the past are basically unknown. Here, we present INP measurements from an ice core in Greenland, which dates back to the year 1370. In total 135 samples were analyzed with the FRIDGE droplet freezing assay in the temperature range from −14 to −35 ∘C. The sampling frequency was set to 1 in 10 years from 1370 to 1960. From 1960 to 1990 the frequency was increased to one sample per year. Additionally, a few special events were probed, including volcanic episodes. The typical time coverage of a sample was on the order of a few months. Historical atmospheric INP concentrations were estimated with a conversion factor, which depends on the snow accumulation rate of the ice core, particle dry deposition velocity, and wet scavenging ratio. Typical atmospheric INP concentrations were on the order of 0.1 L−1 at −25 ∘C. The INP variability was found to be about 1–2 orders of magnitude. Yet, the short-term variability from samples over a seasonal cycle was considerably lower. INP concentrations were significantly correlated to some chemical tracers derived from continuous-flow analysis (CFA) and ion chromatography (IC) over a broad range of nucleation temperatures. The highest correlation coefficients were found for the particle concentration (spherical diameter dp > 1.2 µm). The correlation is higher for a time period of seasonal samples, where INP concentrations follow a clear annual pattern, highlighting the importance of the annual dust input in Greenland from East Asian deserts during spring. Scanning electron microscopy (SEM) analysis of selected samples found mineral dust to be the dominant particle fraction, verifying their significance as INPs. Overall, the concentrations compare reasonably well to present-day INP concentrations, albeit they are on the lower side. However, we found that the INP concentration at medium supercooled temperatures differed before and after 1960. Average INP concentrations at −23, −24, −25, −26, and −28 ∘C were significantly higher (and more variable) in the modern-day period, which could indicate a potential anthropogenic impact, e.g., from land-use change.
Ice nucleating particle concentrations of the past: insights from a
600 year old Greenland ice core
(2020)
Ice nucleating particles (INPs) affect the microphysics in cloud and precipitation processes. Hence, they modulate the radiative properties of clouds. However, atmospheric INP concentrations of the past are basically unknown. Here, we present INP measurements from an ice core in Greenland, which dates back to the year 1370. In total 135 samples were analyzed with the FRIDGE droplet freezing assay in the temperature range from −14 ◦C to −35 ◦C. The sampling frequency was set to 1 in 10 years from 1370 to 1960. From 1960 to 1990 the frequency was increased to 1 sample per year. Additionally, a number of special events were probed, including volcanic episodes. The typical time coverage of a sample was on the order of a few months. Historical atmospheric INP concentrations were estimated with a conversion factor, which depends on the snow accumulation rate of the ice core, particle dry deposition velocity and the wet scavenging ratio. Typical atmospheric INP concentrations were on the order of 0.1 L -1 at −25 ◦C. The INP variability was found to be about 1 – 2 orders of magnitude. Yet, the short-term variability from samples over a seasonal cycle was considerably lower. INP concentrations were significantly correlated to chemical tracers derived from continuous flow analysis (CFA) and ion chromatography (IC) over a broad range of nucleation temperatures. The highest correlation coefficients were found for the particle concentration (dp > 1.2 µm). The correlation is higher for the seasonal samples, where INP concentrations follow a clear annual pattern, highlighting the importance of the annual dust input in Greenland from East Asian deserts during spring. Scanning electron microscopy (SEM) of single particles retrieved from selected samples found particles of soil origin to be the dominant fraction, verifying the significance of mineral dust particles as INPs. Overall, the concentrations compare reasonably well to present day INP concentrations, albeit they are on the lower side. However, we found that the INP concentration at medium supercooled temperatures differed before and after 1960. Average INP concentrations at −23 ◦C, −24 ◦C, −25 ◦C, −26 ◦C and −28 ◦C were significantly higher (and more variable) in the modern day period, which could indicate a potential anthropogenic impact or some post-coring contamination of the topmost, very porous firn.
Rodrigues Ridge connects the Réunion hotspot track with the Central Indian Ridge (CIR) and has been suggested to represent the surface expression of a sub-lithospheric flow channel. From global earthquake catalogues, the seismicity in the region has been associated mainly with events related to the fracture zones at the CIR. However, some segments of the CIR appear void of seismic events. Here, we report on the seismicity recorded at a temporary array of 10 seismic stations operating on Rodrigues Island from September 2014 to June 2016. The array analysis was performed in the time domain by time shifting and stacking the complete waveforms. Event distances were estimated based on a 1-D velocity model and the travel time differences between S and P wave arrivals. We detected and located 63 new events that were not reported by the global networks. Most of the events (51) are located off the CIR and can be classified as intraplate earthquakes. Local magnitudes varied between 1.6 and 3.7. Four seismic clusters were observed that occurred to the west of the spreading segment of the CIR. The Rodrigues Ridge appeared to be aseismic during the period of operation. The lack of seismic activity along both Rodrigues Ridge and the sections of the CIR to the east of Rodrigues may be explained by partially molten upper-mantle material, possibly in relation to the proposed material flow between the Réunion plume and the CIR.
Inclusions of breyite (previously known as walstromite-structured CaSiO3) in diamond are usually interpreted as retrogressed CaSiO3 perovskite trapped in the transition zone or the lower mantle. However, the thermodynamic stability field of breyite does not preclude its crystallization together with diamond under upper-mantle conditions (6–10 GPa). The possibility of breyite forming in subducted sedimentary material through the reaction CaCO3 + SiO2 = CaSiO3 + C + O2 was experimentally evaluated in the CaO–SiO2–C–O2 ± H2O system at 6–10 GPa, 900–1500 ∘C and oxygen fugacity 0.5–1.0 log units below the Fe–FeO (IW) buffer. One experimental series was conducted in the anhydrous subsystem and aimed at determining the melting temperature of the aragonite–coesite (or stishovite) assemblage. It was found that melting occurs at a lower temperature (∼1500 ∘C) than the decarbonation reaction, which indicates that breyite cannot be formed from aragonite and silica under anhydrous conditions and an oxygen fugacity above IW – 1. In the second experimental series, we investigated partial melting of an aragonite–coesite mixture under hydrous conditions at the same pressures and redox conditions. The melting temperature in the presence of water decreased strongly (to 900–1200 ∘C), and the melt had a hydrous silicate composition. The reduction of melt resulted in graphite crystallization in equilibrium with titanite-structured CaSi2O5 and breyite at ∼1000 ∘C. The maximum pressure of possible breyite formation is limited by the reaction CaSiO3 + SiO2 = CaSi2O5 at ∼8 GPa. Based on the experimental results, it is concluded that breyite inclusions found in natural diamond may be formed from an aragonite–coesite assemblage or carbonate melt at 6–8 GPa via reduction at high water activity.
Motivated by the question of whether and how wave–wave interactions should be implemented into atmospheric gravity-wave parametrizations, the modulation of triadic gravity-wave interactions by a slowly varying and vertically sheared mean flow is considered for a non-rotating Boussinesq fluid with constant stratification. An analysis using a multiple-scale WKBJ (Wentzel–Kramers–Brillouin–Jeffreys) expansion identifies two distinct scaling regimes, a linear off-resonance regime, and a nonlinear near-resonance regime. Simplifying the near-resonance interaction equations allows for the construction of a parametrization for the triadic energy exchange which has been implemented into a one-dimensional WKBJ ray-tracing code. Theory and numerical implementation are validated for test cases where two wave trains generate a third wave train while spectrally passing through resonance. In various settings, of interacting vertical wavenumbers, mean-flow shear, and initial wave amplitudes, the WKBJ simulations are generally in good agreement with wave-resolving simulations. Both stronger mean-flow shear and smaller wave amplitudes suppress the energy exchange among a resonantly interacting triad. Experiments with mean-flow shear as strong as in the vicinity of atmospheric jets suggest that internal gravity-wave dynamics are dominated in such regions by wave modulation. However, triadic gravity-wave interactions are likely to be relevant in weakly sheared regions of the atmosphere.