550 Geowissenschaften
Refine
Year of publication
Document Type
- Article (888)
- Doctoral Thesis (194)
- Contribution to a Periodical (32)
- Book (26)
- Working Paper (22)
- Part of Periodical (20)
- 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 (7)
- Klima (7)
- Klimaänderung (7)
- Modellierung (7)
- COSMO-CLM (6)
- Klimawandel (6)
- precipitation (6)
- Atmospheric chemistry (5)
- Deutschland (5)
Institute
- Geowissenschaften (775)
- Geowissenschaften / Geographie (135)
- Biodiversität und Klima Forschungszentrum (BiK-F) (61)
- Geographie (61)
- Senckenbergische Naturforschende Gesellschaft (54)
- Präsidium (44)
- Extern (28)
- Biowissenschaften (22)
- Institut für Ökologie, Evolution und Diversität (8)
- Physik (8)
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.
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.
Forschung Frankfurt : das Wissenschaftsmagazin der Goethe-Universität. 2020, Nr. 2 ; Klimakrise
(2020)
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.
Die Bewertung der Nitrataustragsgefährdung (NAG) landwirtschaftlich genutzter Flächen in Wasserschutzgebieten (WSG) erfolgte bislang auf Basis bodenkundlicher Kartierungen und wurde seit 1996 nach einem im Staatsanzeiger für das Land Hessen veröffentlichten Merkblatt des ehemaligen Hessischen Landesamtes für Bodenforschung im Rahmen der Muster-Wasserschutzgebietsverordnung geregelt (HLfB 1996, HMUJFG 1996). Infolge der Verfügbarkeit hochauflösender Bodendaten in Form der „Bodenflächendaten 1: 5.000, landwirtschaftliche Nutzfläche“ (BFD5L) wird die Ermittlung der Nitrataustragsgefährdung landwirtschaftlich genutzter Flächen neu geregelt. Die BFD5L liefert Auswertungen der Bodenschätzungsdaten zur Feldkapazität des Wurzelraums sowie weiterer relevanter Parameter, die zur Bewertung der Nitrataustragsgefährdung herangezogen werden können.
Um die Eignung der BFD5L-Daten zur Ermittlung der Nitrataustragsgefährdung zu überprüfen, wurden in den Jahren 2009 bis 2012 bodenkundliche Vergleichskartierungen im Rahmen eines Pilotvorhabens im Wasserschutzgebiet Eschollbrücken/Pfungstadt in Südhessen, im Wassereinzugsgebiet der Quelle Meineringhausen bei Korbach, im Wasserschutzgebiet des Tiefbrunnens Spieß der Gemeinde Bad Emstal sowie im WSG Quelle Ohmes der Stadt Kirtorf durchgeführt. Ziel war es, die Umsetzungsmöglichkeiten bei der Nutzung der BFD5LDaten in organisatorischer und technischer Hinsicht zu erproben und das bisherige Verfahren zu überarbeiten (PETER & MILLER 2009, PETER & MILLER 2010a und 2010b, PETER & MILLER 2012).
Die Ergebnisse der Vergleichskartierungen zeigen, dass sich die Daten der BFD5L grundsätzlich für die Ermittlung der Nitrataustragsgefährdung in Wasserschutzgebieten eignen. Lediglich für Flächen, für die nach den bislang im System BFD5L enthaltenen Methoden keine Kennwerte abgeleitet werden können sowie für Sonderstandorte, muss die Nitrataustragsgefährdung durch bodenkundliche Geländearbeiten ermittelt werden.
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.
Anthropogenic climate change is expected to impact ecosystem structure, biodiversity and ecosystem services in Africa profoundly. We used the adaptive Dynamic Global Vegetation Model (aDGVM), which was originally developed and tested for Africa, to quantify sources of uncertainties in simulated African potential natural vegetation towards the end of the 21st century. We forced the aDGVM with regionally downscaled high‐resolution climate scenarios based on an ensemble of six general circulation models (GCMs) under two representative concentration pathways (RCPs 4.5 and 8.5). Our study assessed the direct effects of climate change and elevated CO2 on vegetation change and its plant‐physiological drivers. Total increase in carbon in aboveground biomass in Africa until the end of the century was between 18% to 43% (RCP4.5) and 37% to 61% (RCP8.5) and was associated with woody encroachment into grasslands and increased woody cover in savannas. When direct effects of CO2 on plants were omitted, woody encroachment was muted and carbon in aboveground vegetation changed between –8 to 11% (RCP 4.5) and –22 to –6% (RCP8.5). Simulated biome changes lacked consistent large‐scale geographical patterns of change across scenarios. In Ethiopia and the Sahara/Sahel transition zone, the biome changes forecast by the aDGVM were consistent across GCMs and RCPs. Direct effects from elevated CO2 were associated with substantial increases in water use efficiency, primarily driven by photosynthesis enhancement, which may relieve soil moisture limitations to plant productivity. At the ecosystem level, interactions between fire and woody plant demography further promoted woody encroachment. We conclude that substantial future biome changes due to climate and CO2 changes are likely across Africa. Because of the large uncertainties in future projections, adaptation strategies must be highly flexible. Focused research on CO2 effects, and improved model representations of these effects will be necessary to reduce these uncertainties.
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 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.
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.
A realistic simulation of the atmospheric boundary layer (ABL) depends on an accurate representation of the land–atmosphere coupling. Land surface temperature (LST) plays an important role in this context and the assimilation of LST can lead to improved estimates of the boundary layer and its processes. We assimilated synthetic satellite LST retrievals derived from a nature run as truth into a fully coupled, state‐of‐the‐art land–atmosphere numeric weather prediction model. As assimilation system a local ensemble transform Kalman filter was used and the control vector was augmented by the soil temperature and humidity. To evaluate the concept of the augmented control vector, two‐day case‐studies with different control vector settings were conducted for clear‐sky periods in March and August 2017. These experiments with hourly LST assimilation were validated against the nature run and overall, the RMSE of atmospheric and soil temperature of the first‐guess (and analysis) were reduced. The temperature estimate of the ABL was particularly improved during daytime as was the estimate of the soil temperature during the whole diurnal cycle. The best impact of LST assimilation on the soil and the ABL was achieved with the augmented control vector. Through the coupling between the soil and the atmosphere, the assimilation of LST can have a positive impact on the temperature forecast of the ABL even after 15 hr because of the memory of the soil. These encouraging results motivate further work towards the assimilation of real satellite LST retrievals.
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.
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.
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 (PBAPs) account for a large fraction of coarse particulate matter, and fungal spores are among the most abundant PBAPs in this area as well as in other vegetated continental regions. Furthermore, PBAPs 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 as well as the microclimatic factors temperature and air humidity in parallel to 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 using 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 from 62 % to 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, whereas 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.
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.
An easy-to-use model to evaluate conductivities at high and middle latitudes in the height range 70–100 km is presented. It is based on electron density profiles obtained with the EISCAT VHF radar during 11 years and on the neutral atmospheric model MSIS95. The model uses solar zenith angle, geomagnetic activity and season as input parameters. It was mainly constructed to study the properties of Schumann resonances that depend on such conductivity profiles.
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).
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).
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.
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.
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.
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.
In this study, we use simulations from seven global vegetation models to provide the first multi‐model estimate of fire impacts on global tree cover and the carbon cycle under current climate and anthropogenic land use conditions, averaged for the years 2001–2012. Fire globally reduces the tree covered area and vegetation carbon storage by 10%. Regionally, the effects are much stronger, up to 20% for certain latitudinal bands, and 17% in savanna regions. Global fire effects on total carbon storage and carbon turnover times are lower with the effect on gross primary productivity (GPP) close to 0. We find the strongest impacts of fire in savanna regions. Climatic conditions in regions with the highest burned area differ from regions with highest absolute fire impact, which are characterized by higher precipitation. Our estimates of fire‐induced vegetation change are lower than previous studies. We attribute these differences to different definitions of vegetation change and effects of anthropogenic land use, which were not considered in previous studies and decreases the impact of fire on tree cover. Accounting for fires significantly improves the spatial patterns of simulated tree cover, which demonstrates the need to represent fire in dynamic vegetation models. Based upon comparisons between models and observations, process understanding and representation in models, we assess a higher confidence in the fire impact on tree cover and vegetation carbon compared to GPP, total carbon storage and turnover times. We have higher confidence in the spatial patterns compared to the global totals of the simulated fire impact. As we used an ensemble of state‐of‐the‐art fire models, including effects of land use and the ensemble median or mean compares better to observational datasets than any individual model, we consider the here presented results to be the current best estimate of global fire effects on ecosystems.
Aim: Plant life‐forms characterize key morphological strategies that enable large‐scale comparisons of plant communities. This study applies Raunkiær's plant life‐form concept that was developed for temperate climate to a subtropical island flora, in parts, dominated by summer aridity. We quantify how plant life‐form patterns as well as patterns of important plant functional traits (PFTs) relate to important climate and topographic characteristics.
Location: La Palma, Canary Islands.
Taxon: Flora of La Palma.
Methods: We assigned each native plant species a plant life‐form, that is, phanerophyte, chamaephyte, hemicryptophyte, geophyte and therophyte, as well as PFTs (succulence and N‐fixer). We used stacked species distribution models to assess occurrence probability for each species using the Atlantis database (500 m × 500 m grid). We related richness and percentage values for each plant life‐form and PFT to climate and topography.
Results: Plant life‐forms and PFTs showed a clear pattern within geographic but also climate space, while topography had a minor effect. Phanerophytes mainly contributed to the flora in humid areas. Chamaephytes and hemicryptophytes most strongly contributed to the summit scrub flora and, to some degree, also to the arid coastal regions. Geophytes and therophytes were mainly found in dry coastal regions. N‐fixers contributed mainly to warm‐arid and cool‐arid regions, while succulent species were mainly found in arid coastal regions.
Main conclusions: Raunkiær's plant life‐form concept can be comprehensively transferred to a subtropical island flora by adapting to local unfavourable growing conditions, that is, aridity. Using the strong environmental gradients offered by our study island, we identify substantial climate‐driven variation in patterns of plant life‐forms and PFTs that might be used for large‐scale comparisons in macroecological studies. The growth strategies reflected in Raunkiær's plant life‐forms suggest differences in species establishment and coexistence dynamics within different parts of the island's climate space.
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.
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.
Opportunities and challenges for paleoaltimetry in "small" orogens: insights from the European Alps
(2020)
Many stable isotope paleoaltimetry studies have focused on paleoelevation reconstructions of orogenic plateaus such as the Tibetan or Andean Plateaus. We address the opportunities and challenges of applying stable isotope paleoaltimetry to “smaller” orogens. We do this using a high‐resolution isotope tracking general circulation model (ECHAM5‐wiso) and explore the precipitation δ18O (δ18Op) signal of Cenozoic paleoclimate and topographic change in the European Alps. Results predict a maximum δ18Op change of 4–5‰ (relative to present day) during topographic development of the Alps. This signal of topographic change has the same magnitude as changes in δ18Op values resulting from Pliocene and Last Glacial Maximum global climatic change. Despite the similar magnitude of the isotopic signals resulting from topographic and paleoclimate changes, their spatial patterns across central Europe differ. Our results suggest that an integration of paleoclimate modeling, multiproxy approaches, and low‐elevation reference proxy records distal from an orogen improve topographic reconstructions.
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.
Barriereinseln und Außensände, geformt durch eine Kombination aus Wind, Wellen, Strömung und Küstenlängstransport gelten als morphologisch hoch aktive Küstenbereiche und variieren häufig in Ursprung, Genese und Entwicklung. Sie besitzen durch ihre dissipative Wirkungsweise eine bedeutende Schutzfunktion für rückwärtige Inseln, Halligen und Festlandbereiche und bilden vor der Küste Nordfrieslands die westliche Außengrenze des Wattenmeeres.
Ziel der nachfolgenden Studie ist es, anhand hochauflösender Georadarmessungen und sedimentologischer Daten aus Bohrungen die Landschaftsgeschichte an der Westküste Amrums sowie die Entstehung und interne sedimentäre Architektur des der Insel westlich angelagerten Kniepsandes zu untersuchen und die Prozesse, die zur Genese geführt haben, durch Datierungen zeitlich einzuordnen. Auf Grundlage der gewonnenen Daten wurden zwei stratigraphische Modelle entwickelt, welche die geologisch-geomorphologischen Prozesse und Sedimentationsbedingungen im Westküstenvorfeld erklären.
Ein Modell zeigt die Landschaftsentwicklung an der Westküste Amrums und beschreibt die Sedimentationsbedingungen, die im Vorfeld der Westküste herrschten, bevor der Kniepsand an die Insel heranwanderte. Auf der Landoberfläche des ertrinkenden saaleeiszeitlichen Geestkerns wurden zu Beginn der Flandrischen Transgression feinkörnige Sedimente eines Mischund Schlickwatts abgelagert. Es ist davon auszugehen, dass der damals noch weit vor der Küste Amrums liegende Kniepsand eine Barriere bildete und so an der heute hochenergetischen Westküste für strömungsberuhigte Sedimentationsbedingungen sorgte. Durch Erosion am Geestkern bildeten sich in unterschiedlichen Höhenpositionen fossile Kliffs, die dem damaligen Meeresspiegelstand entsprechen.
Ein weiteres Modell beschreibt den Andockmechanismus des Kniepsandes an die Insel Amrum. Durch die Anlagerung des ehemaligen Außensandes und den damit einhergehenden Sedimentinput wurden die Bedingungen für eine großflächige Dünenbildung geschaffen.
Als erste Ergebnisse der derzeit laufenden paläomagnetischen Untersuchungen im Rhein-Main-Gebiet werden die Befunde von einigen stratigraphisch bedeutenden Profilen mitgeteilt. Während die Remanenzwerte aus altpleistozänen Tonen der Kelsterbacher Terrasse mit einer Ausnahme revers sind, konnte im Autobahneinschnitt Abenheim (NW Worms) und im Lößprofil der Ziegeleigrube Bad Soden am Taunus der Ubergang von reverser zu normaler Magnetisierung gefunden werden. Die Grenze zwischen MATUYAMA- und BRUNHES-Epoche liegt im Profil Bad Soden unter dem 6., das Jaramillo-Event unter dem 7. fossilen Bt-Horizont.
Die über dem letztinterglazialen Boden (Parabraunerde bzw. pseudovergleyte Parabraunerde) folgenden würmzeitlichen äolischen Sedimente können durch kennzeichnende Bodenhorizonte in 3 Abschnitte (Alt-, Mittel- und Jungwürm) gegliedert werden. Charakteristisch für das Altwürm sind Humuszonen, für das Mittelwürm neben einigen „Naßböden" braune Verlehmungszonen bis zu einer Mächtigkeit von 1,1 m und für das Jungwürm mehrere schwach ausgebildete, geringmächtige Verbraunungszonen und „Naßböden". Den wichtigsten Leithorizont des Jungwürms bildet das Kärlicher Tuffbändchen, das in jüngster Zeit auch in Nordhessen aufgefunden wurde. Abschließend wird das in Hessen aufgestellte Gliederungsschema mit den Würmlöß-Gliederungen in anderen Teilen Europas verglichen und eine Parallelisierung versucht.
Der Auffassung, das Substrat der Lockerbraunerde auf dem Oberwald des Vogelsberges sei ausschließlich ein äolisches Sediment der Jüngeren Tundrenzeit, wird widersprochen. Verschiedene Befunde, vor allem 14C-Datierungen, lassen vielmehr den Schluß zu, daß große Teile des Substrates vieler Lockerbraunerden holozäne anthropogene Kolluvien sind. Bekräftigt wird dagegen der Befund gleichen (jungtundrenzeitlichen) Alters für Lockerbraunerde-Substrat und Deckschutt (Decksediment).
Aus dem Gebiet der weichselzeitlichen Vereisung in Polen werden allgemein verbreitete geringmächtige periglaziale Deckschichten beschrieben. Sie zeigen in der Regel eine äolische Beeinflussung und unterscheiden sich dadurch vom Liegenden. Es handelt sich um spätglaziale Bildungen, wie sie aus dem Jungmoränengebiet der DDR seit langem bekannt sind. Ähnliche Substrate wurden auch im nördlichen Alpenvorland gefunden.
Im Lößprofil der Ziegelei Glos bei Lisieux ist eine Abfolge fossiler Boden aufgeschlossen, die sich gut mit der hessischen Lößstratigraphie verbinden läßt. Über dem letztinterglazialen Boden liegen nach einer Diskordanz der Lohner Boden, der E1-, der E2- und der E4-Naß-boden. Letzterer wird häufig als Äquivalent des „Sol de Kesselt" angesehen. Somit zeigt sich auch für dieses Gebiet, daß dieser Boden stratigraphisch nicht dem Lohner Boden entsprechen kann.
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
In der Monte Cavallo-Gruppe fand Verf. in Stauseesedimenten fossile Holzstücke, für die die 14 C-Datierung ein Alter von 29 350 ± 460 Jahren vor 1950 n. Chr. ergab. Die schluffigen Ablagerungen, die in etwa 900 m, fast 80 m über dem heutigen Flußbett des T. Caltea aufgeschlossen sind, ruhen auf Schottern und werden von Moräne überlagert.
Aufgrund seiner stratigraphischen Lage und seines 14C-datierten Alters kann das Holz (Picea abies bzw. Larix) einem zeitlichen Äquivalent des Paudorf-Interstadials zugeordnet werden. Dieser Fund von Großresten ermöglicht somit eine erste absolute Datierung des Paudorf-Interstadials für die Südalpen und Oberitalien.
Das Cranium eines fossilen Hominiden des Formenkreises Homo sapiens sapiens wurde relativ-geologisch sowie absolut durch Radiokohlenstoff und Aminosäuren auf ungefähr 31 000 Jahre B.P. datiert. Andere absolute sowie relative Daten an Mollusken und Mammutzähnen in überlagernden jüngeren Straten datieren auf 18 000 — 21000 und 16 000 Jahren B.P. Geomorphologische und geophysikalische Datierungen stimmen somit gut überein. Er ist der älteste datierte und früheste Bewohner Zentraleuropas, der dem Homo sapiens sapiens angehört.
Aus Franken wird die Entwicklung quartärer Hohlformen beschrieben, deren Rekonstruktion mit Hilfe lößstratigraphischer Methoden (fossile Böden, Tuffbänder, Umlagerungszonen etc.) möglich ist. Bei vielen Formen zeigt sich, daß sie bereits größere Vorläuferformen präwürmzeitlichen Alters hatten. Die Entwicklung während des Würms läßt sich an manchen Beispielen in besonders instruktiver Weise verfolgen. Zu Beginn des Würms, im unteren Mittelwürm und im unteren Jungwürm dominierte zeitweise die Abtragung und Verlagerung. Im oberen Mittelwürm sowie im oberen Jungwürm herrschte äolische Lößsedimentation vor. Diese Ergebnisse stimmen gut mit den bereits aus anderen mitteleuropäischen Lößgebieten bekannten Befunden überein. Mit dem Trockental-System von Helmstadt wird die Entwicklung von Hohlformen beschrieben, deren Anlage bis in das ältere Pleistozän zurückreicht.
In vorliegender Untersuchung wurde der Rißlöß zwischen der 1. und 2. fossilen Parabraunerde anhand schwächerer Bodenbildungen und eingeschalteter Abtragungsphasen zu gliedern versucht. Im jüngeren Riß herrschte starke Lößsedimentation vor, wobei es in mindestens 6 kaltfeuchten Abschnitten zur Ausbildung schwacher periglazialer Naßböden kam. Die Naßbodenserie wurde als Bruchköbeler Böden (B) bezeichnet. Im jüngsten Rißlöß ist wenige dm unter dem Eemboden als tephrochronologischer Leithorizont der Krifteler Tuff (vgl. SEMMEL 1968) eingeschaltet. Den mittleren Profilbereich im Rißlöß zeichnen feuchtere Klimaabschnitte mit starken Verschwemmungsphasen aus, die in den meisten Profilen zu erheblichen Diskordanzen geführt haben. An der Basis der wenigen kompletten Rißlöß-Profile treten vorwiegend in Hessen über der zumeist gekappten 2. fossilen Parabraunerde maximal zwei Schwarzerden auf, die von SEMMEL (1968) als Weilbacher Humuszonen bezeichnet werden. Unmittelbar über diesen Schwarzerden folgt die Ostheimer Zone, eine Fließerde aus aufgearbeitetem Solumaterial der liegenden Böden. Insgesamt zeigt die aus den Rißböden rekonstruierte Klimaabfolge — neben geringfügigen Abweichungen — überraschende Parallelen zur paläopedologisch-klimatischen Gliederung der Würmkaltzeit.
Crustal structure at the western end of the North Anatolian Fault Zone from deep seismic sounding
(2001)
The first deep seismic sounding experiment in Northwestern Anatolia was carried out in October 1991 as part of the "German - Turkish Project on Earthquake Prediction Research" in the Mudurnu area of the North Anatolian Fault Zone. The experiment was a joint enterprise by the Institute of Meteorology and Geophysics of Frankfurt University, the Earthquake Research Institute (ERI) in Ankara, and the Turkish Oil Company (TPAO). Two orthogonal profiles, each 120 km in length with a crossing point near Akyazi, were covered in succession by 30 short period tape recording seismograph stations with 2 km station spacing. 12 shots, with charge sizes between 100 and 250 kg, were fired and 342 seismograms out of 360 were used for evaluation. By coincidence an M b = 4.5 earthquake located below Imroz Island was also recorded and provided additional information on Moho and the sub-Moho velocity. A ray tracing method orginally developed by Weber (1986) was used for travel time inversion. From a compilation of all data two generalized crustal models were derived, one with velocity gradients within the layers and one with constant layer velocities. The latter consists of a sediment cover of about 2 km with V p » 3.6 km/s, an upper crystalline crust down to 13 km with V p » 5.9 km/s, a middle crust down to 25 km depth with V p » 6.5 km/s, a lower crust down to 39 km Moho depth with V p » 7.0 km/s and V p » 8.05 km/s below the Moho. The structure of the individual profiles differs slightly. The thickest sediment cover is reached in the Izmit-Sapanca-trough and in the Akyazi basin. Of particular interest is a step of about 4 km in the lower crust near Lake Sapanca and probably an even larger one in the Moho (derived from the Imroz earthquake data). After the catastrophic earthquake of Izmit on 17 August 1999 this significant heterogeneity in crustal structure appears in a new light with regard to the possible cause of the Izmit earthquake. Heterogeneities in structure are frequently also heterogeneities in strength and stress that impede or even lock rupture. The Izmit earthquake is discussed in relation to a large stepover or jog at the North Anatolian Fault.
This study reports and discusses new radiometric ages, petrographical and volcanological observations and whole rock geochemical data of the rocks of the Rudnik Mts. volcano-intrusive complex. The complex hosts a Pb-Zn-Ag deposit and belongs to the Serbo-Macedonian metallogenetic belt. Two distinct igneous events are distignuished. The first occurred >30 Ma and was mainly characterized by extrusive and shallow intrusive dacites and andesites and was unrelated to mineralization. The second igneous event occurred <23 Ma and was highly heterogeneous in terms of volcanic products and petrographic varieties, but with predominance of quartzlatites. The dacite-andesites (first event) and the quartzlatites (second event) are geochemically similar and display a calc-alkaline affinity and highly incompatible element enriched patterns on spider diagrams, but the younger quartzlatites are richer in K2O, Rb and Ba and poorer in Sr. This is taken as evidence that mixing between an ultrapotassic lamprophyre/lamproite magma and an acid calc-alkaline (dacite-like) magma was essential petrogenetic processes during the second event. The proposed simplified volcanological model suggests that this mixing was responsible for triggering strongly explosive volcanic activity as well as for providing conditions for active hydrothermal and mineralization processes. The observed link between a specific magmatic phase and ore deposit formation can be a general phennomenon in the Balkans, and must be addressed by further and more advanced studies.
This paper provides new geochemical data focusing on valuable elements in the coal, parting, and floor samples in the No. 5 coal seam of the Taiyuan Formation from the Wujiawan mine, Datong coalfield, northern China. The minerals mainly consist of kaolinite, calcite, and pyrite, as well as trace amounts of quartz and illite. The No. 5 coal is enriched in Li, Ga, high field strength elements (HFSEs), and rare earth elements and yttrium (REY) when compared with world hard coals. Of particular interest is the high average concentration of Li (67.66 μg/g), which is around seven times higher than the value for world hard coals. Lithium, Ga, and HFSEs have strong inorganic affinities, whereas REY have organic affinities. The main carrier of Li, Ga, and HFSEs is aluminosilicate minerals, while REY appear to occur with organophosphorus. These HFSEs are enriched, both in the parting and in the adjacent coal samples. This suggests that these elements are likely to leach out during the diagenetic process. The distribution patterns of REY, along with the ratio of Al2O3/TiO2 and the figure of Zr/TiO2 vs. Nb/Y are suggestive of their derivation from felsic parent material. In the northern and eastern part of the Datong coalfield, there are several regions where the Li content is higher than the mineable grade, in particular in the northern Datong coalfield where there is a mine with an Li content of 294.6 μg/g. This is significantly higher than the mineable grade. Therefore, there is a potential for financially viable recovery of Li in these coals of the Datong coalfield.
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).
The Early Permian coal is of great value in the Tengxian Coalfield, Shandon Province, Eastern China. This work deals with the new data focusing on mineralogical characteristics in the Early Permian Shanxi Formation No. 3 coal from the Jinyuan Mine. The Jinyuan coal is a low ash and highly volatile A bituminous coal. Minerals in the No. 3 coal mainly comprise of kaolinite, ankerite, illite, calcite, siderite, and quartz, with varying compositions of trace amounts of pyrite, jarosite, bassanite, anatase, and rutile. According to mineral assemblage in the coal plies, three Types (A to C) can be identified in the No. 3 coal. The dominant minerals in Type A are poorly-ordered kaolinite, illite, quartz, pyrite, and jarosite. Type B is mainly composed of well-ordered kaolinite, illite, siderite, ankerite, and calcite. Type C, with just one sample (JY-3-7c), which contains high proportions of calcite (54%) and ankerite (34%). Terrigenous minerals are elevated in coal plies that typically have relatively high contents of ash yield. The formation of syngenetic pyrite was generally due to seawater, while the sulphate minerals (jarosite and coquimbite) were derived from the oxidation of pyrite. Epigenetic vein-like or fracture-fillings carbonate minerals (ankerite, calcite, and siderite), kaolinite, and pyrite, as well as authigenic quartz were derived from the influx of hydrothermal fluids during different periods, from the authigenic to epigenetic. The paragonite in the coal may have been formed by the precipitated from Na-rich hydrothermal fluids. No effects of magmatic intrusion on mineralogy were investigated in this research.
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.
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
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.
The accurate knowledge of the groundwater storage variation (ΔGWS) is essential for reliable water resource assessment, particularly in arid and semi-arid environments (e.g., Australia, the North China Plain (NCP)) where water storage is significantly affected by human activities and spatiotemporal climate variations. The large-scale ΔGWS can be simulated from a land surface model (LSM), but the high model uncertainty is a major drawback that reduces the reliability of the estimates. The evaluation of the model estimate is then very important to assess its accuracy. To improve the model performance, the terrestrial water storage variation derived from the Gravity Recovery And Climate Experiment (GRACE) satellite mission is commonly assimilated into LSMs to enhance the accuracy of the ΔGWS estimate. This study assimilates GRACE data into the PCRaster Global Water Balance (PCR-GLOBWB) model. The GRACE data assimilation (DA) is developed based on the three-dimensional ensemble Kalman smoother (EnKS 3D), which considers the statistical correlation of all extents (spatial, temporal, vertical) in the DA process. The ΔGWS estimates from GRACE DA and four LSM simulations (PCR-GLOBWB, the Community Atmosphere Biosphere Land Exchange (CABLE), the Water Global Assessment and Prognosis Global Hydrology Model (WGHM), and World-Wide Water (W3)) are validated against the in situ groundwater data. The evaluation is conducted in terms of temporal correlation, seasonality, long-term trend, and detection of groundwater depletion. The GRACE DA estimate shows a significant improvement in all measures, notably the correlation coefficients (respect to the in situ data) are always higher than the values obtained from model simulations alone (e.g., ~0.15 greater in Australia, and ~0.1 greater in the NCP). GRACE DA also improves the estimation of groundwater depletion that the models cannot accurately capture due to the incorrect information of the groundwater demand (in, e.g., PCR-GLOBWB, WGHM) or the unavailability of a groundwater consumption routine (in, e.g., CABLE, W3). In addition, this study conducts the inter-comparison between four model simulations and reveals that PCR-GLOBWB and CABLE provide a more accurate ΔGWS estimate in Australia (subject to the calibrated parameter) while PCR-GLOBWB and WGHM are more accurate in the NCP (subject to the inclusion of anthropogenic factors). The analysis can be used to declare the status of the ΔGWS estimate, as well as itemize the possible improvements of the future model development.
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.
Measurements of halogenated trace gases in 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, maximizing measurement precision. However, individual compounds might be affected during the sampling and 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 high-resolution sampler (HIRES) air sampling unit, which is part of the Civil Aircraft for the Regular Investigation of the atmosphere Based on an Instrument Container (CARIBIC) instrument package. The HIRES unit holds 88 lightweight stainless-steel cylinders that are pressurized in flight to 4.5 bar using metal bellows pumps. The HIRES unit was first deployed in 2010 but has up to now not been used for regular halocarbon analysis with the exception of chloromethane analysis. The sample collection unit was tested for the sampling and storage effects of 28 halogenated compounds. The focus was on compound stability in the stainless-steel canisters during storage of up to 5 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 parts per billion by volume (ppbv). Most of the investigated (hydro)chlorofluorocarbons and long-lived hydrofluorocarbons were found to be stable over a storage time of up to 5 weeks and were unaltered by ozone being present during pressurization. Some compounds such as dichloromethane, trichloromethane, and tetrachloroethene started to decrease in the canisters after a storage time of more than 2 weeks or exhibited lowered mixing ratios in samples pressurized with ozone present. A few compounds such as tetrachloromethane and tribromomethane were found to be unstable in the HIRES stainless-steel canisters independent of ozone levels. Furthermore, growth was observed during storage for some species, namely for HFC-152a, HFC-23, and Halon 1301.
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.
Micromorphology is a suitable method to study the contents and stratigraphic relationships of pit fills. Within the ramparts of Corneşti-Iarcuri, fill layers of a pit were sampled. Th e pit fill was macroscopically divided into primary and secondary fill due to striking differences. These differences could be verified and concretized micromorphologically.
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 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 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.
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.
Inappropriate land management leads to soil loss with destruction of the land’s resource and sediment input into the receiving river. Part of the sediment budget of a catchment is the estimation of soil loss. In the Ruzizi catchment in the Eastern Democratic Republic of the Congo (DRC), only limited research has been conducted on soil loss mainly dealing with local observations on geomorphological forms or river load measurements; a regional quantification of soil loss is missing so far. Such quantifications can be calculated using the Universal Soil Loss Equation (USLE). It is composed of four factors: precipitation (R), soil (K), topography (LS), and vegetation cover (C). The factors can be calculated in different ways according to the characteristics of the study area. In this paper, different approaches for calculating the single factors are reviewed and validated with field work in two sub-catchments of Ruzizi River supplying the water for the reservoirs of Ruzizi I and II hydroelectric dams. It became obvious that the (R)USLE model provides the best results with revised R and LS factors. C factor calculations required to conduct a supervised classification using the Maximum Likelihood Procedure. Different C factor values were assigned to the land cover classes. The calculations resulted in a soil loss rate for the predominantly occurring Ferralsols and Leptosols of around 576 kt/yr in both catchments, when 2016 landcover and precipitation are used. This represents an area-normalized value of 40.4 t/ha/yr for Ruzizi I and 50.5 t/ha/yr for Ruzizi II due to different landcover in the two sub-catchments. The mean value for the whole study area is 47.8 t/ha/yr or even 27.1 t/ha/yr when considering land management techniques like terracing on the slopes (P factor). This work has shown that the (R)USLE model can serve as an easy to handle tool for soil loss quantification when comprehensive field work results are sparse. The model can be implemented in Geographic Information Systems (GIS) with free data; hence, a validation is crucial. It becomes apparent that the use of high resolution Sentinel 2a MSI data as the basis for C factor calculations is an appropriate method for considering heterogeneous Land Use Land Cover (LULC) patterns. To transfer the approach to other regions, the calculation of factor R needs to be modified
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.
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.
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. ...
Under the Kigali Amendment to the Montreal Protocol, new controls are being implemented to reduce emissions of HFC-23 (CHF3), a by-product during the manufacture of HCFC-22 (CHClF2). Starting in 2015, China and India, who dominate global HCFC-22 production (75% in 2017), set out ambitious programs to reduce HFC-23 emissions. Here, we estimate that these measures should have seen global emissions drop by 87% between 2014 and 2017. Instead, atmospheric observations show that emissions have increased and in 2018 were higher than at any point in history (15.9 ± 0.9 Gg yr−1). Given the magnitude of the discrepancy between expected and observation-inferred emissions, it is likely that the reported reductions have not fully materialized or there may be substantial unreported production of HCFC-22, resulting in unaccounted-for HFC-23 by-product emissions. The difference between reported and observation-inferred estimates suggests that an additional ~309 Tg CO2-equivalent emissions were added to the atmosphere between 2015 and 2017.
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.
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.
Invasive plant species are increasingly altering species composition and the functioning of ecosystems from a local to a global scale. The grass species Pennisetum setaceum has recently raised concerns as an invader on different archipelagos worldwide. Among these affected archipelagos are the Canary Islands, which are a hotspot of endemism. Consequently, conservation managers and stakeholders are interested in the potential spreading of this species in the archipelago. We identify the current extent of the suitable habitat for P. setaceum on the island of La Palma to assess how it affects island ecosystems, protected areas (PAs), and endemic plant species richness. We recorded in situ occurrences of P. setaceum from 2010 to 2018 and compiled additional ones from databases at a 500 m × 500 m resolution. To assess the current suitable habitat and possible distribution patterns of P. setaceum on the island, we built an ensemble model. We projected habitat suitability for island ecosystems and PAs and identified risks for total as well as endemic plant species richness. The suitable habitat for P. setaceum is calculated to cover 34.7% of the surface of La Palma. In open ecosystems at low to mid elevations, where native ecosystems are already under pressure by land use and human activities, the spread of the invader will likely lead to additional threats to endemic plant species. Forest ecosystems (e.g., broadleaved evergreen and coniferous forests) are not likely to be affected by the spread of P. setaceum because of its heliophilous nature. Our projection of suitable habitat of P. setaceum within ecosystems and PAs on La Palma supports conservationists and policymakers in prioritizing management and control measures and acts as an example for the potential threat of this graminoid invader on other islands.
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.
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.
Enhanced aridification of Central Asia driven by the combined effects of orogenic surface uplift, Paratethys retreat, changes in atmospheric moisture transport and global cooling is one of the most prominent Cenozoic climate change events of the Northern Hemisphere. Deciphering regional long-term patterns of Central Asian hydrology is, therefore, a key element in understanding the role of Northern Hemisphere mid-latitude drying in the global hydrological system. Pedogenic carbonates record information of the paleoecosystem, the paleogeography, hydrology, tectonic and climatic conditions as well as the sedimentary regime during the time in which they formed. The calcrete-bearing paleosols in continental sedimentary basins yield the possibility for understanding these processes, mechanisms and controlling factors. This study characterizes long-term paleoenvironmental conditions between the late Oligocene and early Miocene in SE Kazakhstan based on stable isotopes, elemental geochemistry and laser ablation U-Pb geochronology from alluvial, fluvial and pedogenic deposits. The main topics addressed in this thesis are paleosol and calcrete formation processes in the light of geochemistry, tectonic and climatic influences on both, as well as the paleoenvironmental and hydrological conditions which led to these processes.
The sedimentological succession of the 14 km-long Kendyrlisai Valley in the Ili Basin, SE Kazakhstan in Central Asia, yields exceptional well-developed paleosols which provide an insight into the relationship between tectonics and sedimentation as well as soil formation processes. The valley accommodates more than 600 m of siliciclastic sediments deposited in a fluvial to alluvial environment in the lower part which grade into lacustrine to fluvial deposition with minor gypsum accumulation in the upper part. The yet undated sediments of Kendyrlisai Valley are compared with a biostratigraphically dated and well described succession in an adjacent location. The correlation of the two successions revealed a proximal-distal relation with the Kendyrlisai Valley succession represents the basin margin position with coarse alluvial-fluvial sediments. The calcrete nodules in Kendyrlisai Valley paleosols give the opportunity of U-Pb dating with LA-ICP-MS which uses the variation in both U and Pb to directly date the calcite. The U-Pb dating revealed a late Oligocene-early Miocene age for the investigated section. The analysis of facies and geochemistry of the paleosol profiles in Kendyrlisai Valley gives insight into calcretization stages, pedogenic processes and landscape stability. Kendyrlisai Valley paleosols show a low to moderate weathering intensity with the highest weathering intensity between 0.5 and 0.8 m depth in the paleosol profile. The comparison between acid leachable and non-acid leachable fraction indicate an incorporation of Ca and Ba in incipient calcrete calcite, whereas clay minerals adsorb Mg and Sr. During the evolution from early to more mature calcrete stages, i.e. calcrete nodules, Ba is lost by dissolution and subsequently adsorbed on clay minerals. The nodules consist almost exclusively of pure calcium carbonate with rare substitution of Ca by Mg indicating a constant supply of Ca by weathering and solution during calcrete formation. The occurrence of Mg-bearing clay minerals indicates high evaporative conditions with alkaline waters (pH >9).
Sampled calcrete nodule cross sections reveal more than one weathering and calcrete formation phase with a higher supply of Ca, Mg and Fe during early phases. The organic and inorganic carbon isotope composition of the calcrete nodules reflects C3 metabolism under occasionally moisture stress, resulting in higher δ13C values and lower respiration rates. This study also presents calculated atmospheric pCO2 values from the calcrete nodules with the equation from Cerling et al. (1999). The average calculated CO2 concentration for Kendyrlisai Valley paleosols is 313 ± 110 ppmv. The comparison with data from other studies during the late Oligocene–early Miocene reveal lower values for Kendyrlisai Valley paleosols, which may stem from an underestimation of the plant-respired CO2.
The knowledge of the variability within the paleosols and calcrete nodules enables a reliable paleoclimatic interpretation for the Kendyrlisai Valley and beyond the Ili Basin. Sedimentary facies and geochemical weathering indices suggest an increased surface and groundwater discharge fed by orographically-enhanced precipitation in the Tien Shan hinterland. In contrast, pedogenic stable isotope data and elevated rates of Mg fixation in clay minerals mirror enhanced rates of evaporation in the vadose zone due to protracted aridification. This study posits that pronounced surface uplift of the Tien Shan Mountains during the Oligocene–Miocene transition promoted regionally increased orographic precipitation and the development of fluvial discharge systems.
The comparison with other studies in adjacent basins creates a precipitation pattern for Central Asia during the late Oligocene to early Miocene. The westerlies supplied Central Asia with atmospheric moisture from the West, possibly from the Eastern Paratethys. The uplifting Tien Shan Mountain ranges captured this moisture on their luv-side, resulting in a pronounced orographic rainfall in the adjacent Ili and Issyk Kul Basins. The Tarim Basin and the Valley of Lakes experienced drier conditions due to the rain shadow effect on the lee-side of the Tien Shan and Altai Mountains. In course of this findings, the thesis highlights the crucial role the Tien Shan Mountains play and had been playing in former times as an orographic barrier for the distribution of atmospheric moisture across Central Asia.
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.
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.
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.
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.
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 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.
The exchange of heat, momentum, and mass in the atmosphere over mountainous terrain is controlled by synoptic-scale dynamics, thermally driven mesoscale circulations, and turbulence. This article reviews the key challenges relevant to the understanding of exchange processes in the mountain boundary layer and outlines possible research priorities for the future. The review describes the limitations of the experimental study of turbulent exchange over complex terrain, the impact of slope and valley breezes on the structure of the convective boundary layer, and the role of intermittent mixing and wave–turbulence interaction in the stable boundary layer. The interplay between exchange processes at different spatial scales is discussed in depth, emphasizing the role of elevated and ground-based stable layers in controlling multi-scale interactions in the atmosphere over and near mountains. Implications of the current understanding of exchange processes over mountains towards the improvement of numerical weather prediction and climate models are discussed, considering in particular the representation of surface boundary conditions, the parameterization of sub-grid-scale exchange, and the development of stochastic perturbation schemes.
Throughout mankind’s history, the need to secure and protect the home settlement was an essential one. This holds especially true for the city of Ainos (modern Enez) in Turkish Thrace. Due to its continuous settlement history since the 7th/6th century BC, several different types of city walls were built—sometimes even on top of each other—several of which have been preserved over time. To decipher the construction style, the course and the age of a buried city wall segment in the southern part of the former city, a geoscientific multi-proxy approach including magnetic gradiometer (MG) and electrical resistivity tomography (ERT) measurements in combination with granulometrical, sedimentological and microfaunistical investigations on sediment cores was applied. We were able to (1) present reasonable arguments for its Hellenistic age; (2) reveal the course of this wall segment and extrapolate it further north into a less studied area; and (3) demonstrate that in this near-coastal area, the former swampy terrain had been consolidated for constructing the wall. Our multi-proxy approach serves as a valuable example for investigating buried structures in archaeological contexts, avoiding a less-economical, time-consuming, or even forbidden excavation.
Atmospheric aerosols and their effect on clouds are thought to be important for anthropogenic radiative forcing of the climate, yet remain poorly understood1. Globally, around half of cloud condensation nuclei originate from nucleation of atmospheric vapours2. It is thought that sulfuric acid is essential to initiate most particle formation in the atmosphere3,4, and that ions have a relatively minor role5. Some laboratory studies, however, have reported organic particle formation without the intentional addition of sulfuric acid, although contamination could not be excluded6,7. Here we present evidence for the formation of aerosol particles from highly oxidized biogenic vapours in the absence of sulfuric acid in a large chamber under atmospheric conditions. The highly oxygenated molecules (HOMs) are produced by ozonolysis of α-pinene. We find that ions from Galactic cosmic rays increase the nucleation rate by one to two orders of magnitude compared with neutral nucleation. Our experimental findings are supported by quantum chemical calculations of the cluster binding energies of representative HOMs. Ion-induced nucleation of pure organic particles constitutes a potentially widespread source of aerosol particles in terrestrial environments with low sulfuric acid pollution.
Atmospheric nanoaerosols have extensive effects on the Earth’s climate and human health. This cumulative work focuses on the development and characterization of instrumentation for measuring various parameters of atmospheric nanoaerosols, and its use to understand new particle formation from organic precursors. The principal research question is, how the chemical composition of nanoaerosol particles can be measured and how atmospheric chemistry influences aerosol processes, especially new particle formation and growth. Therefore, nanoaerosols are investigated under various aspects. More specifically, an instrument is developed to analyze nanoparticles, and field as well as chamber studies are conducted.
The main project is the instrument development of the Thermal Desorption Differential Mobility Analyzer (TD-DMA, project 1, Wagner et al. (2018)). This instrument analyzes the chemical composition of small aerosol particles. By characterization and testing in chamber experiments, it is proven to be suitable for the analysis of freshly nucleated particles.
The second project (Wagner et al. (2017)) applies a broad spectrum of aerosol measurement instruments for the characterization of aerosol particles produced by a skyscraper blasting. A comprehensive picture of the particle population emitted by the demolition is obtained.
Project 3 (K¨urten et al. (2016)) is also an ambient aerosol measurement, focusing of new particle formation in a rural area in central Germany, and the ability of a negative nitrate CI-APi-TOF to detect various substances in atmosphere. Project 4 (Heinritzi et al. (2016)) is a characterization of the negative nitrate CI-APi-TOF used in projects 1, 3, 5, 6, 7 and 8. The following projects focus on understanding new particle formation from atmospherically abundant organic precursors. Key instruments comprise the negative nitrate CI-APiTOF for gas-phase measurements of the nucleating species, and various sizing and counting instruments for quantifying the particle formation and growth. Project 5 (Kirkby et al. (2016)) shows that biogenic organic compounds formed from alpha-pinene can nucleate on their own without the influence of e.g. sulfuric acid. Project 6 (Tr¨ostl et al. (2016)) describe the subsequent growth of these particles. Project 7 (Stolzenburg et al. (2018)) covers the temperature dependence of this growth and in project 8 (Heinritzi et al. (2018)), the suppressing influence of isoprene on the new particle formation is assessed.
We present a study characterizing aerosol particles resulting from a skyscraper blasting. High mass concentrations with a maximum of 844.9 μg m-3 were present for a short time period of approximately 15 minutes. They result in a day mean of 32.6 μg m-3 compared to a 27.6 μg m-3 background not exceeding the 50 μg m-3 EU maximum permissive value. The increase in particle number concentration was less pronounced with a maximum concentration of 6.9 ⋅ 104 cm-3 compared to the local background value of 1.8 ⋅ 104 cm-3. The size-resolved number concentration shows a single mode of ultrafine particles at approximately 93 nm. The spatial distribution of deposited dust was investigated with Bergerhoff glass collection vessels, showing a decrease with distance. In the deposited dust samples the concentrations of twelve metals was determined, non of them exceeded the regional background concentrations significantly. The chemical composition of individual particles emitted by the demolition was studied by Scanning Electron Microscopy. They were mainly concrete and steel particles, with 60% calcium carbonates, 19% calcium sulfates, 19% silicates and 2% steel. In energy-dispersive X-Ray Spectroscopy, no fibers like asbestos were observed. Using a broad spectrum of instruments and methods, we obtain comprehensive characterization of the particles emitted by the demolition.
The overarching goal of the thesis was to create a holistic predictive framework, a vegetation model, by improving the representations of and interactions between the biosphere, hydrosphere, atmosphere and pedosphere. Vegetation models rep- resent a crucial component of Earth system model since the properties of the land surface, via interactions with the atmosphere, can have extremely large climatic effect. Yet, there remains great uncertainty associated with the dynamics of the vegetated land surface. Various vegetation models have been critiqued for numerous reasons including overly simplistic representations of vegetation, prescribed vegetation, poor representations of diversity, inaccurate representations of competition, non-transparent model calibration, and poor responses to drought. The purpose of the creation of this "next generation" model was to address deficiencies common to current vegetation modelling paradigms.
The representation of the biosphere within this framework was improved via two separate development axes. First, ecological realism was improved by integrating concepts from community assembly theory, co-existence theory, and evolutionary theory. Explicitly, rather than defining teleonomic rules to define plant behaviour the process of natural selection is modelled. By modelling the pro- cess of natural selection and its affect on relative fitness, myriad "rules" which continually adapt to biotic and abiotic conditions "come out" as a consequence of the modelled dynamics rather than being "put in". In aDGVM2 (adaptive Dynamic Global Vegetation model 2) communities of plants and their trait values evolve through time, this evolution is constrained by trade-offs between traits. Poorly performing individuals are more likely to die and produce fewer copies of themselves, this results in a filtering of trait values. Further, the community and species’ trait values can evolve through successive generations via reproduction, mutation and crossover which we approximate by using a genetic optimisation algorithm. Thus, a plant community consisting of individuals and species with potentially novel and diverse trait values is assembled iteratively through time.
We tested the assertion that improved integration of concepts from community assembly, evolutionary, and co-existence theory could address limitations of DGVMs in Chapter 2. We demonstrated that such an approach does indeed allow diverse communities of plants to emerge from the modelling framework. We showed that the position of the emergent communities in trait space differed along abiotic gradients and that, in simulations where reproductive isolation was simulated, communities emerged which were composed of multiple co-existing clusters in trait-space. Simulated trait values of co-existing strategies emerging from aDGVM2 were often multimodal, indicative of the emergence of multiple life- history strategies.
Second, to successfully model how natural selection forms a community requires accurate representation of how resource availability affects fitness. In the majority of dynamic global vegetation models (DGVMs) there is no real representation of plant hydraulics with plant water availability being calculated as a simple function of relative soil moisture content and root fractions across a number of soil layers. Worryingly, a number vegetation models appear to under represent the magnitude of these observed responses to drought. This was deficiency was ad- dressed in Chapter 3 by designing a simplified version of the cohesion tension theory of sapwood ascent where elements determining plant conductances are considered in series and implementing a set of trait trade-offs which influence a plant’s hydraulic strategy whereby hydraulic safety trades-off against xylem and leaf conductivity.
Interactions between the biosphere, pedosphere, and hydrosphere can also potentially mediate water resource availability and thus fitness. In the majority of DGVMs the volume of soil explored and explorable by plant roots in fixed glob- ally and usually constrained to a depth not greater than 3m. However, we know that soils can have a strong effect on vegetation distributions, that soil depth is not constant globally, and that plants root to variable depths.
In Chapter 4 I explored interactions between soil depth, plant rooting and the emergent properties of communities and highlighted the importance of considering interactions between the biosphere, hydrosphere, pedosphere, and fire. Here I demonstrated that, in addition to fire and precipitation, edaphic constraints on the volume of soil explorable by plant roots (e.g. by shallow soils, lateritic layers, anoxic conditions due to water logging, toxicity resulting from heavy metal concentrations) can affect the process of plant community assembly, alter the mean values of multiple traits in communities, and the trait diversity of communities.
...
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.
Historically, the expansion of soy plantations has been a major driver of land-use/cover change (LUCC) in Brazil. While a series of recent public actions and supply-chain commitments reportedly curbed the replacement of forests by soy, the expansion of the agricultural commodity still poses a considerable threat to the Amazonian and Cerrado biomes. Identification of areas under high risk of soy expansion is thus paramount to assist conservation efforts in the region. We mapped the areas suitable for undergoing transition to soy plantations in the Legal Amazon with a machine-learning approach adopted from the ecological modeling literature. Simulated soy expansion for the year 2014 exhibited favorable validation scores compared to other LUCC models. We then used our model to simulate how potential future infrastructure improvements would affect the 2014 probabilities of soy occurrence in the region. In addition to the 2.3 Mha of planted soy in the Legal Amazon in 2014, our model identified another 14.7 Mha with high probability of soy conversion in the region given the infrastructure conditions at that time. Out of those, pastures and forests represented 9.8 and 0.4 Mha, respectively. Under the new infrastructure scenarios simulated, the Legal Amazonian area under high risk of soy conversion increased by up to 2.1 Mha (14.6%). These changes led to up to 11.4 and 51.4% increases in the high-risk of conversion areas of pastures and forests, respectively. If conversion occurs in the identified high-risk areas, at least 4.8 Pg of CO2 could be released into the atmosphere, a value that represents 10 times the total CO2 emissions of Brazil in 2014. Our results highlight the importance of targeting conservation policies and enforcement actions, including the Soy Moratorium, to mitigate future forest cover loss associated with infrastructure improvements in the region.
We evaluate the near-surface representation of thermally driven winds in the Swiss Alps in a numerical weather prediction model at km-scale resolution. In addition, the influence of grid resolution (2.2 km and 1.1 km), topography filtering, and land surface datasets on the accuracy of the simulated valley winds is investigated. The simulations are evaluated against a comprehensive set of surface observations for an 18-day fair-weather summer period in July 2006. The episode is characterized by strong diurnal wind systems and the formation of shallow convection over the mountains, which transitions to precipitating convection in some areas. The near-surface winds (10 m above ground level) follow a typical diurnal pattern with strong daytime up-valley flow and weaker nighttime down-valley flow. At a 2.2 km resolution the valley winds are poorly simulated for most stations, while at a 1.1 km resolution the diurnal cycle of the valley winds is well represented in most large (e.g., Rhein valley at Chur and Rhone valley at Visp) and medium-sized valleys (e.g., Linth valley at Glarus). In the smaller valleys (e.g., Maggia valley at Cevio), the amplitude of the valley wind is still significantly underestimated, even at a 1.1 km resolution. Detailed sensitivity experiments show that the use of high-resolution land surface datasets, for both the soil characteristics as well as for the land cover, and reduced filtering of the topography are essential to achieve good performance at a 1.1 km resolution
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.
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.
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.
To quantify water flows between groundwater (GW) and surface water (SW) as well as the impact of Abstract. To quantify water flows between groundwater (GW) and surface water (SW) as well as the impact of capillary rise on evapotranspiration by global hydrological models (GHMs), it is necessary to replace the bucket-like linear GW reservoir model typical for hydrological models with a fully integrated gradient-based GW flow model. Linear reservoir models can only simulate GW discharge to SW bodies, provide no information on the location of the GW table and assume that there is no GW flow among grid cells. A gradient-based GW model simulates not only GW storage but also hydraulic head, which together with information on SW table elevation enables the quantification of water flows from GW to SW and vice versa. In addition, hydraulic heads are the basis for calculating lateral GW flow among grid cells and capillary rise.
G³M is a new global gradient-based GW model with a spatial resolution of 5' that will replace the current linear GW reservoir in the 0.5° WaterGAP Global Hydrology Model (WGHM). The newly developed model framework enables inmemory coupling to WGHM while keeping overall runtime relatively low, allowing sensitivity analyses and data assimilation. This paper presents the G³M concept and specific model design decisions together with results under steady-state naturalized conditions, i.e. neglecting GW abstractions. Cell-specific conductances of river beds, which govern GW-SW interaction, were determined based on the 30'' steady-state water table computed by Fan et al. (2013). Together with an appropriate choice for the effective elevation of the SW table within each grid cell, this enables a reasonable simulation of drainage from GW to SW such that, in contrast to the GW model of de Graaf et al. (2015, 2017), no additional drainage based on externally provided values for GW storage above the floodplain is required in G³M. Comparison of simulated hydraulic heads to observations around the world shows better agreement than de Graaf et al. (2015). In addition, G³M output is compared to the output of two established macro-scale models for the Central Valley, California, and the continental United States, respectively. As expected, depth to GW table is highest in mountainous and lowest in flat regions. A first analysis of losing and gaining rivers and lakes/wetlands indicates that GW discharge to rivers is by far the dominant flow, draining diffuse GW recharge, such that lateral flows only become a large fraction of total diffuse and focused recharge in case of losing rivers and some areas with very low GW recharge. G³M does not represent losing rivers in some dry regions. This study presents the first steps towards replacing the linear GW reservoir model in a GHM while improving on recent efforts, demonstrating the feasibility of the approach and the robustness of the newly developed framework.
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.
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.
Analysis of stratospheric transport from an observational point of view is frequently realized by evaluation of 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 and so far primarily achieved by assumptions about dynamics and spectra themselves. This paper introduces a modified version of an inverse method to infer age spectra from mixing ratios of short-lived trace gases. For a full description of transport seasonality the formulation includes an imposed seasonal cycle to gain multimodal spectra. The EMAC model simulation used for a proof of concept 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. 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 time scales. Results indicate that the modified inverse age spectra match the annual and seasonal pulse age spectra well on global scale beyond 1.5 years 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 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. As the inverse method aims for future implementation on in situ observational data, possible critical factors for this purpose are delineated finally.
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.
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.
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.
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.
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.
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.
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.