Geowissenschaften
Refine
Year of publication
Document Type
- Article (593)
- Doctoral Thesis (139)
- Book (25)
- Contribution to a Periodical (21)
- Conference Proceeding (19)
- Working Paper (19)
- Part of Periodical (10)
- Diploma Thesis (8)
- diplomthesis (7)
- Preprint (6)
- Part of a Book (5)
- Report (5)
- Master's Thesis (4)
- Review (3)
- Bachelor Thesis (2)
- Habilitation (2)
- Other (2)
- Periodical (2)
Is part of the Bibliography
- no (872)
Keywords
- climate change (13)
- Climate change (5)
- Geochemistry (5)
- Klimaänderung (5)
- Atmospheric chemistry (4)
- Boden (4)
- Deutschland (4)
- Klima (4)
- Modellierung (4)
- Stratosphäre (4)
Institute
- Geowissenschaften (872)
- Senckenbergische Naturforschende Gesellschaft (59)
- Präsidium (49)
- Biodiversität und Klima Forschungszentrum (BiK-F) (45)
- Geographie (26)
- Biowissenschaften (15)
- Medizin (11)
- Physik (7)
- Institut für Ökologie, Evolution und Diversität (6)
- Biochemie und Chemie (5)
[Nachruf] Arno Semmel
(2010)
In partially molten regions inside the Earth, melt buoyancy may trigger upwelling of both solid and fluid phases, i.e., diapirism. If the melt is allowed to move separately with respect to the matrix, melt perturbations may evolve into solitary porosity waves. While diapirs may form on a wide range of scales, porosity waves are restricted to sizes of a few times the compaction length. Thus, the size of a partially molten perturbation in terms of compaction length controls whether material is dominantly transported by porosity waves or by diapirism. We study the transition from diapiric rise to solitary porosity waves by solving the two-phase flow equations of conservation of mass and momentum in 2D with porosity-dependent matrix viscosity. We systematically vary the initial size of a porosity perturbation from 1.8 to 120 times the compaction length. If the perturbation is of the order of a few compaction lengths, a single solitary wave will emerge, either with a positive or negative vertical matrix flux. If melt is not allowed to move separately to the matrix a diapir will emerge. In between these end members we observe a regime where the partially molten perturbation will split up into numerous solitary waves, whose phase velocity is so low compared to the Stokes velocity that the whole swarm of waves will ascend jointly as a diapir, just slowly elongating due to a higher amplitude main solitary wave. Only if the melt is not allowed to move separately to the matrix will no solitary waves build up, but as soon as two-phase flow is enabled solitary waves will eventually emerge. The required time to build them up increases nonlinearly with the perturbation radius in terms of compaction length and might be too long to allow for them in nature in many cases.
Auf der Bad Homburg Conference 2021 wurden ausgewählte Fragen der Klimapolitik aus verschiedenen Perspektiven von internationalen Expertinnen und Experten aus Wissenschaft und Zivilgesellschaft, Wirtschaft und Politik diskutiert. Der UniReport hat einige Stimmen zur Konferenz eingeholt, die jeweils wichtige Erkenntnisse, aber auch Streitpunkte und offene Fragen benennen.
In this thesis the Quadrupole Electrical Resistivity Tomography (QERT) method is presented as a new measurement concept for profile-based geoelectric field measurements. The concept is based on a tensorial formulation of the apparent resistivity in order to make three-dimensional statements about the underground conductivity structure. For a simple application of the method a number of similarities to the classical dipole-dipole method were made, such as the presentation of the measurement data in a pseudo-section. The added value of the method compared to the classical profile-based methods is especially the differentiation of lateral structures. Anomalies, which are located laterally to the profile, can be detected with respect to their position (left-right) as well as their conductivity contrast. For the practical implementation of the concept a measuring device was developed and constructed, the CR Device. The device uses 64 channels for simultaneous signal recording of voltage and current time series with up to 1 kHz sampling rate. The current injection is freely programmable and allows any survey design. The measurement of the voltages is performed against a common reference (CR) electrode and thus allows the reconstruction of any dipole voltage by difference formation. A complementary, Matlab-based software package completes the measuring system. An evaluation module allows the raw data of the CR device to be read in, processed and displayed in a suitable form. An inversion module allows the inversion of measurement data into a three-dimensional subsurface model. With a modeling module, measurements over any subsurface situation can be simulated and subsequently analysed. A field measurement on a volcanic maar in the Eifel region, Germany, demonstrates the benefits of the method. A QERT profile was set-up tangentially to a conductive anomaly in the centre of the maar. The measurement data were successfully inverted into a geologically coherent 3D resistivity model.
Analysing the composition of ambient ultrafine particles (UFP) is a challenging task due to the low mass and chemical complexity of small particles, yet it is a prerequisite for the identification ofparticle sources and the assessment of potential health risks. Here, we show the molecular characterization of UFP, based on cascade impactor (Nano-MOUDI) 10samples that were collected at an air quality monitoring station nearby one of Europe`s largest airports in Frankfurt, Germany. At this station, particle-size-distribution measurements show enhanced number concentration of particles smaller than 50nm during airport operating hours. We sampled the lower UFP fraction (0.010-0.018 μm; 0.018-0.032 μm; 0.032-0.056 μm) when the air masses arrived from the airport. We developed an optimized filter extraction procedure, used ultra-high performance liquid chromatography (UHPLC) for compound separation, and a heated electrospray ionization (HESI) source with an 15Orbitrap high-resolution mass spectrometer (HRMS) as a detector for organic compounds. A non-target screening detected ~200 organic compounds in the UFP fraction with sample-to-blank ratios larger than five. We identified the largest signals as homologous series of pentaerythritol esters (PEE) and trimethylolpropane esters (TMPE), which are base stocks of aircraft lubrication oils. We unambiguously attribute the majority of detected compounds to jet engine lubrication oils by matching retention times, high-resolution/accurate mass (HR/AM) measurements, and comparing MS/MS fragmentation patterns between both ambient samples and commercially available jet oils. For each UFP stage, we created molecular fingerprints to visualize the complex chemical composition ofthe organic fraction and their average carbon oxidation state. These graphs underline the presence of the homologous series of PEE and TMPE, and the appearance of jet oil additives (e.g. tricresyl phosphate (TCP)). Targeted screening on TCP confirmed the absence of the harmful tri-orthoisomer, while we identified a thermal transformation product of TMPE-based lubrication oil (trimethylolpropane phosphate (TMP-P)). Even though a quantitative determination of the identified compounds is limited, the presented method enables the qualitative detection of molecular markers for jet engine lubricants in UFP and thus strongly improves the source apportionment of UFP near airports.
Extensive black shale deposits formed in the Early Cretaceous South Atlantic, supporting the notion that this emerging ocean basin was a globally important site of organic carbon burial. The magnitude of organic carbon burial in marine basins is known to be controlled by various tectonic, oceanographic, hydrological, and climatic processes acting on different temporal and spatial scales, the nature and relative importance of which are poorly understood for the young South Atlantic. Here we present new bulk and molecular geochemical data from an Aptian–Albian sediment record recovered from the deep Cape Basin at Deep Sea Drilling Project (DSDP) Site 361, which we combine with general circulation model results to identify driving mechanisms of organic carbon burial. A multi-million year decrease (i.e. Early Aptian–Albian) in organic carbon burial, reflected in a lithological succession of black shale, gray shale, and red beds, was caused by increasing bottom water oxygenation due to abating tectonic restriction via South Atlantic-Southern Ocean gateways. These results emphasize basin evolution and ocean gateway development as a decisive primary control on enhanced organic carbon preservation in the Cape Basin at geological time scales (>1 Myr). The Early Aptian black shale sequence comprises alternations of shales with high (>5%) and relatively low (~3%) organic carbon content of marine sources, the former being deposited during the global Oceanic Anoxic Event (OAE) 1a, as well as during repetitive events before and after OAE 1a. In all cases, these short-term events of enhanced organic carbon burial coincided with strong influxes of sediments derived from the proximal African continent, indicating closely coupled climate–land–ocean interactions. Supported by our model results, we propose that fluctuations in weathering-derived nutrient input from the southern African continent, linked to fluctuations in pCO2 and/or orbitally driven humidity/aridity, were the underlying drivers of short-term organic carbon burial in the deep Cape Basin. These results suggest that deep marine environments of emerging ocean basins responded sensitively and directly to short term fluctuations in riverine nutrient fluxes. We explain this relationship by the lack of wide and mature continental shelf seas that could have acted as a barrier or filter for nutrient transfer from the continent into the deep ocean.
Extensive black shale deposits formed in the Early Cretaceous South Atlantic, supporting the notion that this emerging ocean basin was a globally important site of organic carbon burial. The magnitude of organic carbon burial in marine basins is known to be controlled by various tectonic, oceanographic, hydrological, and climatic processes acting on different temporal and spatial scales, the nature and relative importance of which are poorly understood for the young South Atlantic. Here we present new bulk and molecular geochemical data from an Aptian–Albian sediment record recovered from the deep Cape Basin at Deep Sea Drilling Project (DSDP) Site 361, which we combine with general circulation model results to identify driving mechanisms of organic carbon burial. A multimillion-year decrease (i.e., Early Aptian–Albian) in organic carbon burial, reflected in a lithological succession of black shale, gray shale, and red beds, was caused by increasing bottom water oxygenation due to abating hydrographic restriction via South Atlantic–Southern Ocean gateways. These results emphasize basin evolution and ocean gateway development as a decisive primary control on enhanced organic carbon preservation in the Cape Basin at geological timescales (> 1 Myr). The Early Aptian black shale sequence comprises alternations of shales with high (> 6 %) and relatively low (∼ 3.5 %) organic carbon content of marine sources, the former being deposited during the global Oceanic Anoxic Event (OAE) 1a, as well as during repetitive intervals before and after OAE 1a. In all cases, these short-term intervals of enhanced organic carbon burial coincided with strong influxes of sediments derived from the proximal African continent, indicating closely coupled climate–land–ocean interactions. Supported by our model results, we show that fluctuations in weathering-derived nutrient input from the southern African continent, linked to changes in orbitally driven humidity and aridity, were the underlying drivers of repetitive episodes of enhanced organic carbon burial in the deep Cape Basin. These results suggest that deep marine environments of emerging ocean basins responded sensitively and directly to short-term fluctuations in riverine nutrient fluxes. We explain this relationship using the lack of wide and mature continental shelf seas that could have acted as a barrier or filter for nutrient transfer from the continent into the deep ocean.
In the Central German Uplands, Fagus sylvatica and Picea abies have been particularly affected by climate change. With the establishment of beech forests about 3000 years ago and pure spruce stands 500 years ago, they might be regarded as ‘neophytes’ in the Hessian forests. Palaeoecological investigations at wetland sites in the low mountain ranges and intramontane basins point to an asynchronous vegetation evolution in a comparatively small but heterogenous region. On the other hand, palynological data prove that sustainably managed woodlands with high proportions of Tilia have been persisting for several millennia, before the spread of beech took place as a result of a cooler and wetter climate and changes in land management. In view of increasingly warmer and drier conditions, Tilia cordata appears especially qualified to be an important silvicultural constituent of the future, not only due to its tolerance towards drought, but also its resistance to browsing, and the ability to reproduce vegetatively. Forest managers should be encouraged to actively promote the return to more stress-tolerant lime-dominated woodlands, similar to those that existed in the Subboreal chronozone.
The ICON single-column mode
(2021)
The single-column mode (SCM) of the ICON (ICOsahedral Nonhydrostatic) modeling framework is presented. The primary purpose of the ICON SCM is to use it as a tool for research, model evaluation and development. Thanks to the simplified geometry of the ICON SCM, various aspects of the ICON model, in particular the model physics, can be studied in a well-controlled environment. Additionally, the ICON SCM has a reduced computational cost and a low data storage demand. The ICON SCM can be utilized for idealized cases—several well-established cases are already included—or for semi-realistic cases based on analyses or model forecasts. As the case setup is defined by a single NetCDF file, new cases can be prepared easily by the modification of this file. We demonstrate the usage of the ICON SCM for different idealized cases such as shallow convection, stratocumulus clouds, and radiative transfer. Additionally, the ICON SCM is tested for a semi-realistic case together with an equivalent three-dimensional setup and the large eddy simulation mode of ICON. Such consistent comparisons across the hierarchy of ICON configurations are very helpful for model development. The ICON SCM will be implemented into the operational ICON model and will serve as an additional tool for advancing the development of the ICON model.
Role of the tropical atlantic for the interhemispheric heat transport during the last deglaciation
(2021)
Abstract
During the last deglaciation abrupt millennial-scale perturbations of the Atlantic Meridional Overturning Circulation massively altered the interhemispheric heat distribution affecting, for example, continental ice volume and hydroclimate. If and how the related cross-equatorial heat transport was controlled by the interplay between the southward-flowing Brazil Current (BC) and northward-flowing North Brazil Current (NBC) remains controversial. To assess the role of tropical heat transport during the last deglaciation, we obtained a high-resolution foraminiferal Mg/Ca-based sea surface temperature (SST) record from the BC domain at 21.5°S. The data reveal a yet undocumented warming of at least 4.6°C of the BC during Heinrich Stadial 1 at ∼16 ka indicating massive oceanic heat accumulation in the tropical South Atlantic. Simultaneously, a strongly diminished NBC prevented the release of this excess heat into the northern tropics. The observed magnitude of heat accumulation substantially exceeds numerical model simulations, stressing the need to further scrutinize atmospheric and oceanic heat transport during extreme climatic events.
Plain Language Summary
The Atlantic overturning circulation underwent abrupt millennial-scale perturbations. Such phases of sluggish oceanic circulation resulted in a substantial reduction of northward heat transport. As a consequence, substantial cooling occurred in the Northern Hemisphere and warming occurred in the Southern Hemisphere with severe effects on tropical precipitation. The distribution of heat within the western tropical Atlantic is accomplished by the southward-flowing BC and the northward-flowing NBC. By reconstructing SSTs for the interval between 20,000 and 10,000 yr before present, we assess the role of both currents in the interhemispheric heat transport during weak Atlantic overturning. We found that a sluggish overturning circulation resulted in anomalous southward heat transport by the BC in concert with a weak NBC, which lead to a yet undocumented warming of at least 4.6°C in the western tropical South Atlantic. This warming significantly exceeds reconstructions based on numerical simulations. This points to the need to further improve our understanding of changes in the cross-equatorial oceanic and atmospheric heat transport in response to rapid changes in ocean circulation, in particular as a significant weakening of the Atlantic overturning circulation is predicted in the wake of anthropogenic climate change.
Production and use of many synthetic halogenated trace gases are regulated internationally due to their contribution to stratospheric ozone depletion or climate change. In many applications they have been replaced by shorter-lived compounds, which have become measurable in the atmosphere as emissions increased. Non-target monitoring of trace gases rather than targeted measurements of well-known substances is needed to keep up with such changes in the atmospheric composition. We regularly deploy gas chromatography (GC) coupled to time-of-flight mass spectrometry (TOF-MS) for analysis of flask air samples and in situ measurements at the Taunus Observatory, a site in central Germany. TOF-MS acquires data over a continuous mass range that enables a retrospective analysis of the dataset, which can be considered a type of digital air archive. This archive can be used if new substances come into use and their mass spectrometric fingerprint is identified. However, quantifying new replacement halocarbons can be challenging, as mole fractions are generally low, requiring high measurement precision and low detection limits. In addition, calibration can be demanding, as calibration gases may not contain sufficiently high amounts of newly measured substances or the amounts in the calibration gas may have not been quantified. This paper presents an indirect data evaluation approach for TOF-MS data, where the calibration is linked to another compound which could be quantified in the calibration gas. We also present an approach to evaluate the quality of the indirect calibration method, select periods of stable instrument performance and determine well suited reference compounds. The method is applied to three short-lived synthetic halocarbons: HFO-1234yf, HFO-1234ze(E), and HCFO-1233zd(E). They represent replacements for longer-lived hydrofluorocarbons (HFCs) and exhibit increasing mole fractions in the atmosphere.
The indirectly calibrated results are compared to directly calibrated measurements using data from TOF-MS canister sample analysis and TOF-MS in situ measurements, which are available for some periods of our dataset. The application of the indirect calibration method on several test cases can result in uncertainties of around 6 % to 11 %. For hydro(chloro-)fluoroolefines (denoted H(C)FOs), uncertainties up to 23 % are achieved. The indirectly calculated mole fractions of the investigated H(C)FOs at Taunus Observatory range between measured mole fractions at urban Dübendorf and Jungfraujoch stations in Switzerland.
Paläoklimarekonstruktionen, die es sich zum Ziel gesetzt haben, Klima-Mensch Interaktionen auf lange Zeitreihen betrachtet zu erforschen, nehmen begünstigt durch die aktuell intensiv geführte Klimadebatte, einen immer größer werdenden Stellenwert in der öffentlichen und wissenschaftlichen Wahrnehmung ein. Denn trotz aller wissenschaftlicher Fortschritte, die in den vergangenen Jahrzehnten im Bereich der modernen Klimaforschung gemacht wurden, bleibt die zuverlässige Vorhersage und Modellierung von zukünftigen Klimaveränderungen noch immer eine der größten Herausforderungen unser heutigen Zeit. Betrachtet man die Karibik exemplarisch in diesem Rahmen, dann prognostizieren viele Modellrechnungen, infolge steigender Ozeantemperaturen, ein deutlich häufigeres Auftreten von tropischen Stürmen und Hurrikanen sowie eine Verschiebung hin zu höheren Sturmstärken. Dieser Trend stellt für die Karibik und viele daran angrenzende Staaten eine der größten Gefahren des modernen Klimawandels dar, den es wissenschaftlich über einen langen Zeitrahmen zu erforschen gilt.
Klimaprognosen stützen sich meist vollständig auf hoch-aufgelöste instrumentelle Datensätze. Diese sind aber alle durch einen wesentlichen Aspekt limitiert. Aufgrund ihrer eingeschränkten Verfügbarkeit (~150 Jahre) fehlt ihnen die erforderliche Tiefe, um die auf langen Zeitskalen operierenden Prozesse der globalen Klimadynamik adäquat abbilden zu können. Betrachtet man das Holozän in seiner Gesamtheit, so wurde die globale Klimadynamik über die vergangenen ~11,700 Jahre von periodisch auftretenden Prozessen und Abläufen gesteuert. Diese wirken grundsätzlich über Zeiträume von mehreren Jahrzehnten, teilweise Jahrhunderten und in einigen Fällen sogar Jahrtausenden. Viele dieser natürlichen Prozesse, können in der kurzen Instrumentellen Ära nicht gänzlich identifiziert und angemessen in Klimamodellen berücksichtig werden. Die alleinige Berücksichtigung der Instrumentellen Ära bietet daher nur eine eingeschränkte Perspektive, um die Ursachen und Abläufe von vergangenen sowie mögliche Folgen von zukünftigen Klimaveränderungen zu verstehen. Um diese Einschränkung zu überwinden, ist es somit erforderlich, dass die geowissenschaftliche Forschung mit Proxymethoden ein zusammenfassendes und mechanistisches Verständnis über alle Holozänen Klimaveränderungen erlangt.
Wenn man sich diese Limitierung, die ansteigenden Ozeantemperaturen und das in der Karibik in den vergangen 20 Jahren vermehrte Auftreten von starken tropischen Zyklonen ins Gedächtnis ruft, ist es nachvollziehbar, dass im Rahmen dieser Doktorarbeit ein zwei Jahrtausende langer und jährlich aufgelöster Klimadatensatz erarbeitet werden soll, der spät Holozäne Variationen von Ozeanoberflächenwasser-temperaturen (SST) und daraus resultierende lang-zeitliche Veränderungen in der Häufigkeit tropischer Zyklone widerspiegelt. In Zentralamerika wird das Ende der Maya Hochkultur (900-1100 n.Chr.) mit drastischen Umweltveränderungen (z.B. Dürren) assoziiert, die während der Mittelalterlichen Warmzeit (MWP; 900-1400 n.Chr.) durch eine globale Klimaveränderung hervorgerufen wurde. Die aus einem „Blue Hole“ abgeleiteten Informationen über Klimavariationen der Vergangenheit können als Referenz für die gegenwärtige Klimakriese verwendet werden.
Als „Blue Hole“ wird eine Karsthöhle bezeichnet, die sich subaerisch während vergangener Meeresspiegeltiefstände im karbonatischen Gerüst eines Riffsystems gebildet hat und in Folge eines Meeresspiegelanstiegs vollständig überflutet wurde. In einigen wenigen marinen „Blue Holes“ treten anoxische Bodenwasserbedingungen auf. Die in diesen anoxischen Karsthöhlen abgelagerten Abfolgen mariner Sedimente können als einzigartiges Klimaarchiv verwendet werden, da sie aufgrund des Fehlens von Bioturbation eine jährliche Schichtung (Warvierung) aufweisen.
In dieser kumulativen Dissertation über das „Great Blue Hole“ werden die Ergebnisse eines 3-jährigen Forschungsprojekts vorgestellt, dass das Ziel verfolgte einen wissenschaftlich herausragenden spät Holozänen Klimadatensatz für die süd-westliche Karibik zu erzeugen. Beim „Great Blue Hole“ handelt es sich um ein weltweit einzigartiges marines Sedimentarchiv für diverse spät Holozäne Klima-veränderungen, das im Zuge dieser Dissertation sowohl nach paläoklimatischen als auch nach sedimentologischen Fragestellungen untersucht wurde. Die vorliegende Doktorarbeit befasst sich im Einzelnen mit (1) der Ausarbeitung eines jährlich aufgelösten Archives für tropische Zyklone, (2) der Entwicklung eines jährlich aufgelösten SST Datensatzes und (3) einer kompositionellen Quantifizierung der sedimentären Abfolgen sowie einer faziell-stratigraphischen Charakterisierung von Schönwetter-Sedimenten und Sturmlagen. Zu jedem dieser drei Aspekte, wurde jeweils ein Fachartikel bei einer anerkannten wissenschaftlichen Fachzeitschrift mit „peer-review“ Verfahren veröffentlicht.
Der insgesamt 8.55 m lange Sedimentbohrkern („BH6“), der für diese Dissertation untersucht wurde, stammt vom Boden des 125 m tiefen und 320 m breiten „Great Blue Holes“, das sich in der flachen östlichen Lagune des 80 km vor der Küste von Belize (Zentralamerika) gelegenen „Lighthouse Reef“ Atolls befindet. Durch seine besondere Geomorphologie wirkt das, innerhalb des atlantischen „Hurrikan Gürtels“ positionierte, „Great Blue Hole“ wie eine gigantische Sedimentfalle. Die unter Schönwetter-Bedingungen kontinuierlich abgelagerten Abfolgen feinkörniger karbonatischer Sedimente, werden von groben Sturmlagen unterbrochen, die auf „over-wash“ Prozesse von tropischen Zyklonen zurückzuführen sind.
...
Abstract
The mineralogy, chemical composition, and physical properties of cratonic mantle eclogites with oceanic crustal protoliths can be modified by secondary processes involving interaction with fluids and melts, generated in various slab lithologies upon subduction (auto‐metasomatism) or mantle metasomatism after emplacement into the cratonic lithosphere. Here we combine new and published data to isolate these signatures and evaluate their effects on the chemical and physical properties of eclogite. Mantle metasomatism involving kimberlite‐like, ultramafic carbonated melts (UM carbonated melts) is ubiquitous though not pervasive, and affected between ~20% and 40% of the eclogite population at the various localities investigated here, predominantly at ~60–150 km depth, overlapping cratonic midlithospheric seismic discontinuities. Its hallmarks include lower jadeite component in clinopyroxene and grossular component in garnet, an increase in bulk‐rock MgO ± SiO2, and decrease in FeO and Al2O3 contents, and LREE‐enrichment accompanied by higher Sr, Pb, Th, U, and in part Zr and Nb, as well as lower Li, Cu ± Zn. This is mediated by addition of a high‐temperature pyroxene from a UM carbonated melt, followed by redistribution of this component into garnet and clinopyroxene. As clinopyroxene‐garnet trace‐element distribution coefficients increase with decreasing garnet grossular component, clinopyroxene is the main carrier of the metasomatic signatures. UM carbonated melt‐metasomatism at >130–150 km has destroyed the diamond inventory at some localities. These mineralogical and chemical changes contribute to low densities, with implications for eclogite gravitational stability, but negligible changes in shear‐wave velocities, and, if accompanied by H2O‐enrichment, will enhance electrical conductivities compared to unenriched eclogites.
Plain Language Summary
Oceanic crust formed at spreading ridges is recycled in subduction zones and undergoes metamorphism to eclogite. Some of this material is captured in the overlying lithospheric mantle, where it is exhumed by passing magmas. Having formed in spreading ridges, these eclogites have proven invaluable archives for the onset of plate tectonics, for the construction of cratons during subduction/collision, as probes of the convecting mantle from which their precursors formed, and as generators of heterogeneity upon recycling into Earth's convecting mantle. During subduction and until exhumation, interaction with fluids and melts (called metasomatism) can change the mineralogy, chemical composition, and physical properties of mantle eclogites, complicating their interpretation, but a comprehensive study of these effects is lacking so far. We investigated mantle eclogites from ancient continents (cratons) around the globe in order to define hallmarks of metasomatism by subduction‐related fluids and small‐volume ultramafic carbonated mantle melts. We find that the latter is pervasive and occurs predominantly at midlithospheric depths where seismic discontinuities are detected, typically causing diamond destruction and a reduction in density. This has consequences for their gravitational stability and for the interpretation of shearwave velocities in cratons.
Der Paläontologe und das Meer : Eberhard Gischler erforscht das Klimagedächtnis von Korallenriffen
(2018)
Constraining the architecture of complex 3D volcanic plumbing systems within active rifts, and their impact on rift processes, is critical for examining the interplay between faulting, magmatism and magmatic fluids in developing rift segments. The Natron basin of the East African Rift System provides an ideal location to study these processes, owing to its recent magmatic-tectonic activity and ongoing active carbonatite volcanism at Oldoinyo Lengai. Here, we report seismicity and fault plane solutions from a 10-month temporary seismic network spanning Oldoinyo Lengai, Naibor Soito volcanic field and Gelai volcano. We locate 6827 earthquakes with ML -0.85 to 3.6, which are related to previous and ongoing magmatic and volcanic activity in the region, as well as regional tectonic extension. We observe seismicity down to ~17 km depth north and south of Oldoinyo Lengai and shallow seismicity (3 - 10 km) beneath Gelai, including two swarms. The deepest seismicity (~down to 20 km) occurs above a previously imaged magma body below Naibor Soito. These seismicity patterns reveal a detailed image of a complex volcanic plumbing system, supporting potential lateral and vertical connections between shallow- and deep-seated magmas, where fluid and melt transport to the surface is facilitated by intrusion of dikes and sills. Focal mechanisms vary spatially. T-axis trends reveal dominantly WNW-ESE extension near Gelai, while strike-slip mechanisms and a radial trend in P-axes are observed in the vicinity of Oldoinyo Lengai. These data support local variations in the state of stress, resulting from a combination of volcanic edifice loading and magma-driven stress changes imposed on a regional extensional stress field. Our results indicate that the southern Natron basin is a segmented rift system, in which fluids preferentially percolate vertically and laterally in a region where strain transfers from a border fault to a developing magmatic rift segment.
The complete elastic stiffness tensor of thiourea has been determined from thermal diffuse scattering (TDS) using high-energy photons (100 keV). Comparison with earlier data confirms a very good agreement of the tensor coefficients. In contrast with established methods to obtain elastic stiffness coefficients (e.g. Brillouin spectroscopy, inelastic X-ray or neutron scattering, ultrasound spectroscopy), their determination from TDS is faster, does not require large samples or intricate sample preparation, and is applicable to opaque crystals. Using high-energy photons extends the applicability of the TDS-based approach to organic compounds which would suffer from radiation damage at lower photon energies.
Strain localization in the lithosphere and the formation, evolution, and maintenance of resulting plate boundaries play a crucial role in plate tectonics and thermo‐chemical mantle convection. Previously activated lithospheric deformation zones often appear to maintain a “memory” of weakening, leading to tectonic inheritance within plate reorganizations including the Wilson cycle. Different mechanisms have been proposed to explain such strain localization, but it remains unclear which operates on what spatio‐temporal scales, and how to best incorporate them in large‐scale mantle convection models. Here, we analyze two candidates, (1), grain‐size sensitive rheology and, (2), damage‐style parameterizations of yield, stress which are sometimes used to approximate the former. Grain‐size reduction due to dynamic recrystallization can drive localization in the ductile domain, and grain growth provides a time‐dependent rheological hardening component potentially enabling the preservation of rheological heterogeneities. We compare the dynamic weakening and hardening effects as well as the timescales of strength evolution for a composite rheology including grain‐size dynamics with a pseudo‐plastic rheology including damage‐ (or “strain”‐) dependent weakening. We explore the implications of different proposed grain‐size evolution laws, and test to which extent strain‐dependent rheologies can mimic the weakening and hardening effects of the more complex micro‐physical behavior. Such an analysis helps to better understand the parallels and differences between various strain‐localization modeling approaches used in different tectonics and geodynamics communities. More importantly, our results contribute to efforts to identify the key ingredients of strain‐localization and damage hysteresis within plate tectonics and how to represent those in planetary‐scale modeling.
The Ceboruco is a 2280 m high stratovolcano located in Nayarit State, Mexico. Despite its last eruption which occurred in 1870, it is the most active volcano in the area, showing volcanicearthquake activity together with ongoing vapor emissions. The magnetotelluric survey was carried out in November 2016. It was part of a geothermal project (CeMIEGeo-P24) and focused on the determination of the electrical conductivity distribution in the subsurface of the volcano.
The Magnetotelluric Apparent Resistivity Tensor, as introduced by Brown (2016), can be decomposed into an amplitude and a phase tensor. The fundamental physics behind those new tensors were presented in Hering et al. (2019), using canonical models in 1-D (isotropic and anisotropic) and 2-D resistivity environments. Here, the tensors are introduced for a high-quality data set, where their interpretational benefits become very obvious. Additionally, results from an isotropic 3-D inversion are presented and compared to an alternative 3-D anisotropic forward model.
During the first two days of August 2016 a seismic crisis occurred on Brava, Cape Verde, which – according to observations based on a local seismic network – was characterized by more than thousand volcano–seismic signals. Brava is considered an active volcanic island, although it has not experienced any historic eruptions. Seismicity significantly exceeded the usual level during the crisis. We report on results based on data from a temporary seismic–array deployment on the neighbouring island of Fogo at a distance of about 35 km. The array was in operation from October 2015 to December 2016 and recorded a total of 1343 earthquakes, 355 thereof were localized. On 1 and 2 August we observed 54 earthquakes, 25 of which could be located beneath Brava. We further evaluate the observations with regards to possible precursors to the crisis and its continuation. Our analysis shows a migration of seismicity around Brava, but no distinct precursory pattern. However, the observations suggest that similar earthquake swarms commonly occur close to Brava. The results further confirm the advantages of seismic arrays as tools for the remote monitoring of regions with limited station coverage or access.
During the first two days of August 2016 a seismic crisis occurred on Brava, Cabo Verde, which – according to observations based on a local seismic network – was characterized by more than a thousand volcano-seismic signals. Brava is considered an active volcanic island, although it has not experienced any historic eruptions. Seismicity significantly exceeded the usual level during the crisis. We report on results based on data from a temporary seismic-array deployment on the neighbouring island of Fogo at a distance of about 35 km. The array was in operation from October 2015 to December 2016 and recorded a total of 1343 earthquakes in the region of Fogo and Brava; 355 thereof were localized. On 1 and 2 August we observed 54 earthquakes, 25 of which could be located beneath Brava. We further evaluate the observations with regards to possible precursors to the crisis and its continuation. Our analysis shows a significant variation in seismicity around Brava, but no distinct precursory pattern. However, the observations suggest that similar earthquake swarms commonly occur close to Brava. The results further confirm the advantages of seismic arrays as tools for the remote monitoring of regions with limited station coverage or access.
Rodrigues Ridge connects the Réunion hotspot track with the Central Indian Ridge (CIR) and has been suggested to represent the surface expression of a sub-lithospheric flow channel. From global earthquake catalogues, the seismicity in the region has been associated mainly with events related to the fracture zones at the CIR. However, some segments of the CIR appear void of seismic events. Here, we report on the seismicity recorded at a temporary array of 10 seismic stations operating on Rodrigues Island from September 2014 to June 2016. The array analysis was performed in the time domain by time shifting and stacking the complete waveforms. Event distances were estimated based on a 1-D velocity model and the travel time differences between S and P wave arrivals. We detected and located 63 new events that were not reported by the global networks. Most of the events (51) are located off the CIR and can be classified as intraplate earthquakes. Local magnitudes varied between 1.6 and 3.7. Four seismic clusters were observed that occurred to the west of the spreading segment of the CIR. The Rodrigues Ridge appeared to be aseismic during the period of operation. The lack of seismic activity along both Rodrigues Ridge and the sections of the CIR to the east of Rodrigues may be explained by partially molten upper-mantle material, possibly in relation to the proposed material flow between the Réunion plume and the CIR.
Tetra-auricupride, ideally AuCu, represents the only species showing the coexistence of Au with an elevated level of Pt, as in the case of a detrital grain studied structurally for the first time, from an ophiolite-associated placer at Bolshoy Khailyk, western Sayans, Russia. We infer that tetra-auricupride can incorporate as much as ~30 mol. % of a “PtCu” component, apparently without significant modification of the unit cell. The unit-cell parameters of platiniferous tetra-auricupride are: a 2.790(1) Å, c 3.641(4) Å, with c/a = 1.305, which are close to those reported for ordered AuCu(I) in the system Au–Cu, and close also to the cell parameters of tetraferroplatinum (PtFe), which both appear to crystallize in the same space group, P4/mmm. These intermetallic compounds and natural alloys are thus isostructural. The closeness of their structures presumably allows Pt to replace Au atoms so readily. The high extent of Cu + Au enrichment is considered to be a reflection of geochemical evolution and buildup in levels of the incompatible Cu and Au with subordinate Pt in a remaining volume of melt at low levels of fO2 and fS2 in the system.
In the nineteenth century, two Neolithic axe-heads were reported from the Michelsberg enclosure system at Kapellenberg. The recent identification of an unusually large tumulus, from which the axe-heads were almost certainly once recovered, reveals that socio-political hierarchisation, linked to the emergence of high-ranking elites in Brittany and the Paris Basin during the fifth millennium cal BC, may have extended into Central Europe.
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.
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.
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.
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.
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.
Aim: Predicting future changes in species richness in response to climate change is one of the key challenges in biogeography and conservation ecology. Stacked species distribution models (S‐SDMs) are a commonly used tool to predict current and future species richness. Macroecological models (MEMs), regression models with species richness as response variable, are a less computationally intensive alternative to S‐SDMs. Here, we aim to compare the results of two model types (S‐SDMS and MEMs), for the first time for more than 14,000 species across multiple taxa globally, and to trace the uncertainty in future predictions back to the input data and modelling approach used.
Location: Global land, excluding Antarctica.
Taxon: Amphibians, birds and mammals.
Methods: We fitted S‐SDMs and MEMs using a consistent set of bioclimatic variables and model algorithms and conducted species richness predictions under current and future conditions. For the latter, we used four general circulation models (GCMs) under two representative concentration pathways (RCP2.6 and RCP6.0). Predicted species richness was compared between S‐SDMs and MEMs and for current conditions also to extent‐of‐occurrence (EOO) species richness patterns. For future predictions, we quantified the variance in predicted species richness patterns explained by the choice of model type, model algorithm and GCM using hierarchical cluster analysis and variance partitioning.
Results: Under current conditions, species richness predictions from MEMs and S‐SDMs were strongly correlated with EOO‐based species richness. However, both model types over‐predicted areas with low and under‐predicted areas with high species richness. Outputs from MEMs and S‐SDMs were also highly correlated among each other under current and future conditions. The variance between future predictions was mostly explained by model type.
Main conclusions: Both model types were able to reproduce EOO‐based patterns in global terrestrial vertebrate richness, but produce less collinear predictions of future species richness. Model type by far contributes to most of the variation in the different future species richness predictions, indicating that the two model types should not be used interchangeably. Nevertheless, both model types have their justification, as MEMs can also include species with a restricted range, whereas S‐SDMs are useful for looking at potential species‐specific responses.
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.
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.
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 crystal structure of MgCO3-II has long been discussed in the literature where DFT-based model calculations predict a pressure-induced transition of the carbon atom from the sp2 to the sp3 type of bonding. We have now determined the crystal structure of iron-bearing MgCO3-II based on single-crystal X-ray diffraction measurements using synchrotron radiation. We laser-heated a synthetic (Mg0.85Fe0.15)CO3 single crystal at 2500 K and 98 GPa and observed the formation of a monoclinic phase with composition (Mg2.53Fe0.47)C3O9 in the space group C2/m that contains tetrahedrally coordinated carbon, where CO44− tetrahedra are linked by corner-sharing oxygen atoms to form three-membered C3O96− ring anions. The crystal structure of (Mg0.85Fe0.15)CO3 (magnesium iron carbonate) at 98 GPa and 300 K is reported here as well. In comparison with previous structure-prediction calculations and powder X-ray diffraction data, our structural data provide reliable information from experiments regarding atomic positions, bond lengths, and bond angles.
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.
Vegetation responds to drought through a complex interplay of plant hydraulic mechanisms, posing challenges for model development and parameterization. We present a mathematical model that describes the dynamics of leaf water-potential over time while considering different strategies by which plant species regulate their water-potentials. The model has two parameters: the parameter λ describing the adjustment of the leaf water potential to changes in soil water potential, and the parameter Δψww describing the typical ‘well-watered’ leaf water potentials at non-stressed (near-zero) levels of soil water potential. Our model was tested and calibrated on 110 time-series datasets containing the leaf- and soil water potentials of 66 species under drought and non-drought conditions. Our model successfully reproduces the measured leaf water potentials over time based on three different regulation strategies under drought. We found that three parameter sets derived from the measurement data reproduced the dynamics of 53% of an drought dataset, and 52% of a control dataset [root mean square error (RMSE) < 0.5 MPa)]. We conclude that, instead of quantifying water-potential-regulation of different plant species by complex modeling approaches, a small set of parameters may be sufficient to describe the water potential regulation behavior for large-scale modeling. Thus, our approach paves the way for a parsimonious representation of the full spectrum of plant hydraulic responses to drought in dynamic vegetation models.
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.
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.
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.