550 Geowissenschaften
Refine
Year of publication
Document Type
- Article (896) (remove)
Has Fulltext
- yes (896)
Keywords
- Climate change (7)
- climate change (7)
- COSMO-CLM (6)
- Palaeoclimate (6)
- Atmospheric chemistry (5)
- Biogeochemistry (5)
- Palaeoceanography (5)
- precipitation (5)
- loess (4)
- Clausius–Clapeyron scaling (3)
Institute
- Geowissenschaften (546)
- Geowissenschaften / Geographie (100)
- Biodiversität und Klima Forschungszentrum (BiK-F) (59)
- Senckenbergische Naturforschende Gesellschaft (51)
- Extern (22)
- Biowissenschaften (21)
- Geographie (15)
- Präsidium (12)
- Institut für Ökologie, Evolution und Diversität (9)
- Physik (6)
Analyzing the impact of streamflow drought on hydroelectricity production: a global-scale study
(2021)
Electricity production by hydropower is negatively affected by drought. To understand and quantify risks of less than normal streamflow for hydroelectricity production (HP) at the global scale, we developed an HP model that simulates time series of monthly HP worldwide and thus enables analyzing the impact of drought on HP. The HP model is based on a new global hydropower database (GHD), containing 8,716 geo-localized plant records, and on monthly streamflow values computed by the global hydrological model WaterGAP with a spatial resolution of 0.5°. The GHD includes 44 attributes and covers 91.8% of the globally installed capacity. The HP model can reproduce HP trends, seasonality, and interannual variability that was caused by both (de)commissioning of hydropower plants and hydrological variability. It can also simulate streamflow drought and its impact on HP reasonably well. Global risk maps of HP reduction were generated for both 0.5° grid cells and countries, revealing that 67 out of the 134 countries with hydropower suffer, in 1 out of 10 years, from a reduction of more than 20% of mean annual HP and 18 countries from a reduction of more than 40%. The developed HP model enables advanced assessments of drought impacts on hydroelectricity at national to international levels.
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.
AirCore samplers have been increasingly used to capture vertical profiles of trace gases reaching from the ground up to about 30 km, in order to validate remote sens- ing instruments and to investigate transport processes in the stratosphere. When deployed to a weather balloon, accu- rately attributing the trace gas measurements to the sampling altitudes is nontrivial, especially in the stratosphere. In this paper we present the CO-spiking experiment, which can be deployed to any AirCore on any platform in order to evalu- ate different computational altitude attribution processes and to experimentally derive the vertical resolution of the profile by injecting small volumes of signal gas at predefined GPS altitudes during sampling. We performed two CO-spiking flights with an AirCore from the Goethe University Frankfurt (GUF) deployed to a weather balloon in Traînou, France, in June 2019. The altitude retrieval based on an instantaneous pressure equilibrium assumption slightly overestimates the sampling altitudes, especially at the top of the profiles. For these two flights our altitude attribution is accurate within 250 m below 20 km. Above 20 km the positive bias becomes larger and reaches up to 1.2 km at 27 km altitude. Differences in descent velocities are shown to have a major impact on the altitude attribution bias. We parameterize the time lag between the theoretically attributed altitude and the actual CO-spike release altitude for both flights together and use it to empirically correct our AirCore altitude retrieval. Regard- ing the corrected profiles, the altitude attribution is accurate within ±120 m throughout the profile. Further investigations are needed in order to test for the scope of validity of this correction parameter regarding different ambient conditions and maximum flight altitudes. We derive the vertical resolu- tion from the CO spikes of both flights and compare it to the modeled vertical resolution. The modeled vertical resolution is too optimistic compared to the experimentally derived res- olution throughout the profile, albeit agreeing within 220 m. All our findings derived from the two CO-spiking flights are strictly bound to the GUF AirCore dimensions. The newly introduced CO-spiking experiment can be used to test differ- ent combinations of AirCore configurations and platforms in future studies.
Acesta excavata (Fabricius, 1779) is a slow growing bivalve from the Limidae family and is often found associated with cold-water coral reefs along the European continental margin. Here we present the compositional variability of frequently used proxy elemental ratios (Mg/ Ca, Sr/Ca, Na/Ca) measured by laser-ablation mass spectrometry (LA-ICP-MS) and com- pare it to in-situ recorded instrumental seawater parameters such as temperature and salin- ity. Shell Mg/Ca measured in the fibrous calcitic shell section was overall not correlated with seawater temperature or salinity; however, some samples show significant correlations with temperature with a sensitivity that was found to be unusually high in comparison to other marine organisms. Mg/Ca and Sr/Ca measured in the fibrous calcitic shell section display significant negative correlations with the linear extension rate of the shell, which indicates strong vital effects in these bivalves. Multiple linear regression analysis indicates that up to 79% of elemental variability is explicable with temperature and salinity as independent pre- dictor values. Yet, the overall results clearly show that the application of Element/Ca (E/Ca) ratios in these bivalves to reconstruct past changes in temperature and salinity is likely to be complicated due to strong vital effects and the effects of organic material embedded in the shell. Therefore, we suggest to apply additional techniques, such as clumped isotopes, in order to exactly determine and quantify the underlying vital effects and possibly account for these. We found differences in the chemical composition between the two calcitic shell lay- ers that are possibly explainable through differences of the crystal morphology. Sr/Ca ratios also appear to be partly controlled by the amount of magnesium, because the small magne- sium ions bend the crystal lattice which increases the space for strontium incorporation. Oxi- dative cleaning with H2O2 did not significantly change the Mg/Ca and Sr/Ca composition of the shell. Na/Ca ratios decreased after the oxidative cleaning, which is most likely a leaching effect and not caused by the removal of organic matter.
The future physiology of marine phytoplankton will be impacted by a range of changes in global ocean conditions, including salinity regimes that vary spatially and on a range of short- to geological timescales. Coccolithophores have global ecological and biogeochemical significance as the most important calcifying marine phytoplankton group. Previous research has shown that the morphology of their exoskeletal calcified plates (coccoliths) responds to changing salinity in the most abundant coccolithophore species, Emiliania huxleyi. However, the extent to which these responses may be strain-specific is not well established. Here we investigated the growth response of six strains of E. huxleyi under low (ca. 25) and high (ca. 45) salinity batch culture conditions and found substantial variability in the magnitude and direction of response to salinity change across strains. Growth rates declined under low and high salinity conditions in four of the six strains but increased under both low and high salinity in strain RCC1232 and were higher under low salinity and lower under high salinity in strain PLYB11. When detailed changes in coccolith and coccosphere size were quantified in two of these strains that were isolated from contrasting salinity regimes (coastal Norwegian low salinity of ca. 30 and Mediterranean high salinity of ca. 37), the Norwegian strain showed an average 26% larger mean coccolith size at high salinities compared to low salinities. In contrast, coccolith size in the Mediterranean strain showed a smaller size trend (11% increase) but severely impeded coccolith formation in the low salinity treatment. Coccosphere size similarly increased with salinity in the Norwegian strain but this trend was not observed in the Mediterranean strain. Coccolith size changes with salinity compiled for other strains also show variability, strongly suggesting that the effect of salinity change on coccolithophore morphology is likely to be strain specific. We propose that physiological adaptation to local conditions, in particular strategies for plasticity under stress, has an important role in determining ecotype responses to salinity.
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.
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.
This article presents the findings from systematically reviewing 26 empirical research studies published from 2005 to 2014 on the use of GIS for learning and teaching. By employing methods of narrative synthesis and qualitative content analysis, the study gives evidence about the state of knowledge of competence‐based GIS education. The results explain what factors and variables effect GIS learning in terms of technology use, major subject contents, learning contexts, and didactic and pedagogical aspects. They also show what facets of knowledge, process skills, and affect the research literature has investigated. The analysis of the type and quality of the methods used indicates that current GIS education research is a heterogeneous field that needs a systematic research framework for future efforts, according to empirical education research.
Good quality data on precipitation are a prerequisite for applications like short-term weather forecasts, medium-term humanitarian assistance, and long-term climate modelling. In Sub-Saharan Africa, however, the meteorological station networks are frequently insufficient, as in the Cuvelai-Basin in Namibia and Angola. This paper analyses six rainfall products (ARC2.0, CHIRPS2.0, CRU-TS3.23, GPCCv7, PERSIANN-CDR, and TAMSAT) with respect to their performance in a crop model (APSIM) to obtain nutritional scores of a household’s requirements for dietary energy and further macronutrients. All products were calibrated to an observed time series using Quantile Mapping. The crop model output was compared against official yield data. The results show that the products (i) reproduce well the Basin’s spatial patterns, and (ii) temporally agree to station records (r = 0.84). However, differences exist in absolute annual rainfall (range: 154 mm), rainfall intensities, dry spell duration, rainy day counts, and the rainy season onset. Though calibration aligns key characteristics, the remaining differences lead to varying crop model results. While the model well reproduces official yield data using the observed rainfall time series (r = 0.52), the products’ results are heterogeneous (e.g., CHIRPS: r = 0.18). Overall, 97% of a household’s dietary energy demand is met. The study emphasizes the importance of considering the differences among multiple rainfall products when ground measurements are scarce.
The Global South is facing severe challenges in ensuring livelihood security due to climate change impacts, environmental degradation and population growth as well as changing lifestyles. These complex problems cannot be solely solved by single scientific disciplines – they require transdisciplinary research (TDR). Stakeholders from civil society, the corporate sector, government and science need to pool their knowledge to find solutions for sustainable transformations. In Namibia, we have been involved in TDR projects on water supply, and sanitation services as well as livestock management in rangeland systems. In this paper, we review two TDR projects that differ in multiple ways and hence allow us to carve out structural differences and critically discuss research outcomes, lessons learned and the challenge of North–South collaborations. Our review builds upon published and unpublished project documents as well as expert interviews with Namibian and German researchers who were involved in the projects. Our results show that TDR can be put into practice in different ways, depending on the research focus and the period available. The TDR phases of problem framing, inter- and transdisciplinary integration were implemented with different tools and foci points. We discuss the role of project length and funding conditions for project success and outcome generation. In addition, we critically consider the role of Namibian and German researchers in these international collaborations. The conclusions we draw touch upon the points of preparatory research funding, the equal acknowledgement of Global South contributions to joint research projects and the explicit handling of TDR components in project work. Significance: • The current social-ecological challenges are complex and require TDR as a mode of knowledge coproduction, particularly in a development context. • Inter- and transdisciplinary integration are critical processes for a project to be successful and require the allocation of adequate time and monetary resources. • Longer-term projects with a funded preparatory research phase constitute a structural model for TDR as project outcomes can evolve over time. • Global South researchers carry a hidden burden in international collaborations that has to be adequately acknowledged upfront in project planning and final products.
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.
2019 hat Frankfurt einen Rekord geknackt: Am 25. Juli wurde im Stadtteil Westend die Tageshöchsttemperatur von 40,2 Grad Celsius erreicht. Damit war Frankfurt der heißeste Ort Hessens seit Beginn der Wetteraufzeichnungen im Jahr 1881. Ein Superlativ, der zeigt: Der Klimawandel stellt auch die Städte vor große Herausforderungen. Gesucht werden Wege, um mit seinen Folgen zurechtzukommen.
Aus der Redaktion
(2020)
Mit dem Klimawandel könnte das Grundwasser künftig weltweit in vielen Regionen knapp werden. Das zeigt eine große internationale Studie, die vom Doktoranden Robert Reinecke und der Geographin Prof. Petra Döll initiiert wurde. Vielerorts könnte das zu Wassermangel führen oder einen bestehenden Wassermangel verschärfen. Gleichzeitig werden andere Regionen unter steigenden Grundwasserständen zu leiden haben. Ackerflächen könnten durch Vernässung verloren gehen.
Schon zu früheren Zeiten in der Erdgeschichte gab es Warmzeiten durch starke Treibhauseffekte, mit tropischen Temperaturen in weiten Teilen der Erde, hohem Meeresspiegel und massivem Artensterben. Das belegen Daten aus der Paläoklimatologie. Wenn man heutige Klimamodelle auf solche geologischen Warmzeiten anwendet, kann man sie testen und verbessern. So verhilft die Paläoklimatologie zu einem präziseren Blick in unsere Klimazukunft.
Forschung Frankfurt : das Wissenschaftsmagazin der Goethe-Universität. 2020, Nr. 2 ; Klimakrise
(2020)
Wetlands such as bogs, swamps, or freshwater marshes are hotspots of biodiversity. For 5.1 million km2 of inland wetlands, the dynamics of area and water storage, which strongly impact biodiversity and ecosystem services, were simulated using the global hydrological model WaterGAP. For the first time, the impacts of both human water use and man‐made reservoirs (WUR) and future climate change (CC) on wetlands around the globe were quantified. WUR impacts are concentrated in arid/semiarid regions, where WUR decreased mean wetland water storage by more than 5% on 8.2% of the mean wetland area during 1986–2005 (Am), with highest decreases in groundwater depletion area. Using output of three climate models, CC impacts on wetlands were quantified, distinguishing unavoidable impacts [i.e., at 2 °C global warming (GW)] from avoidable impacts (difference between 3 °C and 2 °C impacts). Even unavoidable CC impacts are projected to be much larger than WUR impacts, also in arid/semiarid regions. On most wetland area with reliable estimates, avoidable CC impacts are more than twice as large as unavoidable impacts. In case of 2 °C GW, half of Am is estimated to be unaffected by mean storage changes of more than 5%, but only one third in case of 3 °C GW. Temporal variability of water storage will increase for most wetlands. Wetlands in dry regions will be affected the most, particularly by water storage decreases in the dry season. Different from wealthier countries, low‐income countries will dominantly suffer from a decrease in wetland water storage due to CC.
Anthropogenic climate change is expected to impact ecosystem structure, biodiversity and ecosystem services in Africa profoundly. We used the adaptive Dynamic Global Vegetation Model (aDGVM), which was originally developed and tested for Africa, to quantify sources of uncertainties in simulated African potential natural vegetation towards the end of the 21st century. We forced the aDGVM with regionally downscaled high‐resolution climate scenarios based on an ensemble of six general circulation models (GCMs) under two representative concentration pathways (RCPs 4.5 and 8.5). Our study assessed the direct effects of climate change and elevated CO2 on vegetation change and its plant‐physiological drivers. Total increase in carbon in aboveground biomass in Africa until the end of the century was between 18% to 43% (RCP4.5) and 37% to 61% (RCP8.5) and was associated with woody encroachment into grasslands and increased woody cover in savannas. When direct effects of CO2 on plants were omitted, woody encroachment was muted and carbon in aboveground vegetation changed between –8 to 11% (RCP 4.5) and –22 to –6% (RCP8.5). Simulated biome changes lacked consistent large‐scale geographical patterns of change across scenarios. In Ethiopia and the Sahara/Sahel transition zone, the biome changes forecast by the aDGVM were consistent across GCMs and RCPs. Direct effects from elevated CO2 were associated with substantial increases in water use efficiency, primarily driven by photosynthesis enhancement, which may relieve soil moisture limitations to plant productivity. At the ecosystem level, interactions between fire and woody plant demography further promoted woody encroachment. We conclude that substantial future biome changes due to climate and CO2 changes are likely across Africa. Because of the large uncertainties in future projections, adaptation strategies must be highly flexible. Focused research on CO2 effects, and improved model representations of these effects will be necessary to reduce these uncertainties.
Precipitation extremes with devastating socioeconomic consequences within the South American Monsoon System (SAMS) are expected to become more frequent in the near future. The complexity in SAMS behavior, however, poses severe challenges for reliable future projections. Thus, robust paleomonsoon records are needed to constrain the high spatiotemporal variability in the response of SAMS rainfall to different climatic drivers. This study uses Ti/Ca ratios from X‐ray fluorescence scanning of a sediment core retrieved off eastern Brazilian to trace precipitation changes over the past 322 Kyr. The results indicate that despite the spatiotemporal complexity of the SAMS, insolation forcing is the primary pacemaker of variations in the monsoonal system. Additional modulation by atmospheric pCO2 suggests that SAMS intensity over eastern Brazil will be suppressed by rising CO2 emissions in the future. Lastly, our record reveals an unprecedented strong and persistent wet period during Marine Isotope Stage 6 driven by anomalously strong trade winds.
The complex magnetotelluric (MT) apparent resistivity tensor can be decomposed into two real tensors, the apparent resistivity and the resistivity phase tensors, which represent relationships between the observed electric field at a point on the Earth's surface and an associated apparent current density. We explain the differences between these tensors and conventional estimates of apparent resistivity and phase for simple resistivity environments and demonstrate, using canonical models in 1‐D and 2‐D environments, that both tensors are more sensitive to vertical and horizontal resistivity gradients than their conventional counterparts. The properties of the new tensors are explained using electromagnetic induction theory and the effects of associated charges at resistivity boundaries. We introduce a new way to plot tensor ellipses, which brings significant improvements to the interpretation of MT data, using appropriate visualization software. The apparent resistivity tensor gives information about the magnitude and direction of apparent resistivity subsurface structures and has a strong response to vertical resistivity contrasts. The resistivity phase tensor is highly sensitive to vertical boundaries and the associated fields in the TM mode. It is also free from static distortions under the same conditions implied for the conventional phase tensor. These findings have prompted a study in the potential of the new tensors for 3‐D inversions. The results from a 3‐D inversion of a canonical oblique conductor straddling two quarter spaces show distinct improvements in resolving the boundaries of the conductor and open a promising field for future studies.
The endemic argan woodlands cover large parts of South Morocco and create a characteristic landscape with areas of sparsely vegetated and bare soil surfaces between single trees. This unique ecosystem has been under extensive agrosilvopastoral management for centuries and is now at risk of degradation caused by overgrazing and increasing scarcity and variability of rainfall.
To investigate susceptibility to wind erosion, we conducted an experimental–empirical study including wind tunnel tests and a drone‐generated digital elevation model and quantified wind‐erodible material on five different associated surface types by means of sediment catchers. The highest emission flux was measured on freshly ploughed surfaces (1875 g m–2 h–1), while older ploughed areas with a re‐established crust produced a much lower emission flux (795 g m–2 h–1). Extensive tillage may have been a sustainable practice for generations, but increasing drought and uncertainty of rainfall now lead to an acute risk of severe soil erosion and dust production. The typical crusted surfaces characterized by residual rock fragment accumulation and wash processes produced the second highest emission flux (1,354 g m–2 h–1). Material collected from tree‐shaded areas (933 g m–2 h–1) was revealed to be a considerable source of organic material, possibly affecting substrate conditions positively on a larger regional scale. The lowest flux was measured on rock fragment‐covered surfaces (301 g m–2 h–1).
The data show that open argan woodland may be a considerable source for wind erosion and dust production, depending on surface characteristics strongly related to management. An adapted management must include the conservation of argan trees to offer a promising approach to prevent severe wind erosion and dust production and mitigate possible impacts of land‐use change and climate change related shifts in wind and rainfall patterns.
A realistic simulation of the atmospheric boundary layer (ABL) depends on an accurate representation of the land–atmosphere coupling. Land surface temperature (LST) plays an important role in this context and the assimilation of LST can lead to improved estimates of the boundary layer and its processes. We assimilated synthetic satellite LST retrievals derived from a nature run as truth into a fully coupled, state‐of‐the‐art land–atmosphere numeric weather prediction model. As assimilation system a local ensemble transform Kalman filter was used and the control vector was augmented by the soil temperature and humidity. To evaluate the concept of the augmented control vector, two‐day case‐studies with different control vector settings were conducted for clear‐sky periods in March and August 2017. These experiments with hourly LST assimilation were validated against the nature run and overall, the RMSE of atmospheric and soil temperature of the first‐guess (and analysis) were reduced. The temperature estimate of the ABL was particularly improved during daytime as was the estimate of the soil temperature during the whole diurnal cycle. The best impact of LST assimilation on the soil and the ABL was achieved with the augmented control vector. Through the coupling between the soil and the atmosphere, the assimilation of LST can have a positive impact on the temperature forecast of the ABL even after 15 hr because of the memory of the soil. These encouraging results motivate further work towards the assimilation of real satellite LST retrievals.
The continental expression of global cooling during the Miocene Climate Transition in Central Asia is poorly documented, as the tectonically active setting complicates the correlation of Neogene regional and global climatic developments. This study presents new geochemical data (CaSO4 content, carbonate δ13C and δ18O) from the endorheic alluvial‐lacustrine Aktau succession (Ili Basin, south‐east Kazakhstan) combined with findings from the previously published facies evolution. Time series analysis revealed long‐eccentricity forcing of the paleohydrology throughout the entire succession, split into several facies‐dependent segments. Orbital tuning, constrained by new laser ablation U‐Pb dates and a preexisting magnetostratigraphy, places the succession in a 5.0 Ma long interval in the middle to late Miocene (15.6 to 10.6 Ma). The long‐term water accumulation in the Ili Basin followed the timing of the Miocene Climate Transition, suggesting increased precipitation in the catchment area in response to climate cooling and stronger westerly winds. This was paced by minima of the 2.4 Ma eccentricity cycle, which favored the establishment of a discharge playa (~14.3 Ma) and a perennial lake (12.6 to 11.8 Ma). Furthermore, low obliquity amplitudes (nodes) caused a transient weakening of the westerlies at ~13.7 to 13.5 Ma and at ~12.7 Ma, resulting in negative hydrological budgets and salinization. Flooding of the windward Ili Basin coeval with aridification in the leeward basins suggests that the Tian Shan was a climate boundary already in the middle Miocene. Our results emphasize the impact of climate fluctuations on the westerlies' strength and thus on Central Asian hydrology.
Bioaerosols are considered to play a relevant role in atmospheric processes, but their sources, properties and spatiotemporal distribution in the atmosphere are not yet well characterized. In the Amazon Basin, primary biological aerosol particles (PBAP) account for a large fraction of coarse particulate matter, and fungal spores are among the most abundant PBAP there as well as in other vegetated continental regions. furthermore, PBAP could also be important ice nuclei in Amazonia. Measurement data on the release of fungal spores under natural conditions, however, are sparse. Here we present an experimental approach to analyze and quantify the spore release from fungi and other spore producing organisms under natural and laboratory conditions. For measurements under natural conditions, the samples were kept in their natural environment and a setup was developed to estimate the spore release numbers and sizes together with the microclimatic factors temperature and air humidity, as well as the mesoclimatic parameters net radiation, rain, and fog occurrence. For experiments in the laboratory, we developed a cuvette to assess the particle size and number of newly released fungal spores under controlled conditions, simultaneously measuring temperature and relative humidity inside the cuvette. Both approaches were combined with bioaerosol sampling techniques to characterize the released particles by microscopic methods. For fruiting bodies of the basidiomycetous species, Rigidoporus microporus, the model species for which these techniques were tested, the highest frequency of spore release occurred in the range of 62 and 96 % relative humidity. The results obtained for this model species reveal characteristic spore release patterns linked to environmental or experimental conditions, indicating that the moisture status of the sample may be a regulating factor, while temperature and light seem to play a minor role for this species. The presented approach enables systematic studies aimed at the quantification and validation of spore emission rates and inventories, which can be applied to a regional mapping of cryptogamic organisms under given environmental conditions.
Bioaerosols are considered to play a relevant role in atmospheric processes, but their sources, properties, and spatiotemporal distribution in the atmosphere are not yet well characterized. In the Amazon Basin, primary biological aerosol particles (PBAPs) account for a large fraction of coarse particulate matter, and fungal spores are among the most abundant PBAPs in this area as well as in other vegetated continental regions. Furthermore, PBAPs could also be important ice nuclei in Amazonia. Measurement data on the release of fungal spores under natural conditions, however, are sparse. Here we present an experimental approach to analyze and quantify the spore release from fungi and other spore-producing organisms under natural and laboratory conditions. For measurements under natural conditions, the samples were kept in their natural environment and a setup was developed to estimate the spore release numbers and sizes as well as the microclimatic factors temperature and air humidity in parallel to the mesoclimatic parameters net radiation, rain, and fog occurrence. For experiments in the laboratory, we developed a cuvette to assess the particle size and number of newly released fungal spores under controlled conditions, simultaneously measuring temperature and relative humidity inside the cuvette. Both approaches were combined with bioaerosol sampling techniques to characterize the released particles using microscopic methods. For fruiting bodies of the basidiomycetous species, Rigidoporus microporus, the model species for which these techniques were tested, the highest frequency of spore release occurred in the range from 62 % to 96 % relative humidity. The results obtained for this model species reveal characteristic spore release patterns linked to environmental or experimental conditions, indicating that the moisture status of the sample may be a regulating factor, whereas temperature and light seem to play a minor role for this species. The presented approach enables systematic studies aimed at the quantification and validation of spore emission rates and inventories, which can be applied to a regional mapping of cryptogamic organisms under given environmental conditions.
An easy-to-use model to evaluate conductivities at high and middle latitudes in the height range 70–100 km is presented. It is based on electron density profiles obtained with the EISCAT VHF radar during 11 years and on the neutral atmospheric model MSIS95. The model uses solar zenith angle, geomagnetic activity and season as input parameters. It was mainly constructed to study the properties of Schumann resonances that depend on such conductivity profiles.
The turnover time of terrestrial carbon (τ) controls the global carbon cycle – climate feedback and, yet, is poorly simulated by the current Earth System Models (ESMs). In this study, by assessing apparent carbon turnover time as the ratio between carbon stocks and fluxes, we provide a new, updated ensemble of diagnostic terrestrial carbon turnover times and associated uncertainties on a global scale using multiple, state-of-the-art, observation-based datasets of soil organic carbon stock (Csoil), vegetation biomass (Cveg) and gross primary productivity (GPP). Using this new ensemble, we estimated the global average τ to be 42$% &' years when the full soil depth is considered, longer than the previous estimates of 23$) &* years. Only considering the top 1 m (assuming maximum active layer depth is up to 1 meter) of soil carbon in circumpolar regions yields a global τ of 35$) &' years. Csoil in circumpolar regions account for two thirds of the total uncertainty in global τ estimates, whereas Csoil in non-circumpolar contributes merely 9.38%. GPP (2.25%) and Cveg (0.05%) contribute even less to the total uncertainty. Therefore, the high uncertainty in Csoil is the main factor behind the uncertainty in global τ, as reflected in the larger range of full-depth Csoil (3152-4372 PgC). The uncertainty is especially high in circumpolar regions with a behaviour of ESMs which could contribute to uncertainty reductions in future projections of the carbon cycle - climate feedback. The dataset of the terrestrial turnover time ensemble (DOI: 10.17871/bgitau.201911) is openly available from the data portal: https://doi.org/10.17871/bgitau.201911 (Fan et al., 2019) uncertainty of 50% and the spatial correlations among different datasets are also low compared to other regions. Overall, we argue that current global datasets do not support robust estimates of τ globally, for which we need clarification on variations of Csoil with soil depth and stronger estimates of Csoil in circumpolar regions. Despite the large variation in both magnitude and spatial patterns of τ, we identified robust features in the spatial patterns of τ that emerge regardless of soil depth and differences in data sources of Csoil, Cveg and GPP. Our findings show that the latitudinal gradients of τ are consistent across different datasets and soil depth. Furthermore, there is a strong consensus on the negative correlation between τ and temperature along latitude that is stronger in temperate zones (30ºN-60ºN) than in subtropical and tropical zones (30ºS30ºN). The identified robust patterns can be used to infer the response of τ to climate and for constraining contemporaneous behaviour of ESMs which could contribute to uncertainty reductions in future projections of the carbon cycle - climate feedback. The dataset of the terrestrial turnover time ensemble (DOI:10.17871/bgitau.201911) is openly available from the data portal: https://doi.org/10.17871/bgitau.201911 (Fan et al., 2019).
The turnover time of terrestrial ecosystem carbon is an emergent ecosystem property that quantifies the strength of land surface on the global carbon cycle–climate feedback. However, observation- and modeling-based estimates of carbon turnover and its response to climate are still characterized by large uncertainties. In this study, by assessing the apparent whole ecosystem carbon turnover times (τ) as the ratio between carbon stocks and fluxes, we provide an update of this ecosystem level diagnostic and its associated uncertainties in high spatial resolution (0.083∘) using multiple, state-of-the-art, observation-based datasets of soil organic carbon stock (Csoil), vegetation biomass (Cveg) and gross primary productivity (GPP). Using this new ensemble of data, we estimated the global median τ to be 43+7−7 yr (median+difference to percentile 75−difference to percentile 25) when the full soil is considered, in contrast to limiting it to 1 m depth. Only considering the top 1 m of soil carbon in circumpolar regions (assuming maximum active layer depth is up to 1 m) yields a global median τ of 37+3−6 yr, which is longer than the previous estimates of 23+7−4 yr (Carvalhais et al., 2014). We show that the difference is mostly attributed to changes in global Csoil estimates. Csoil accounts for approximately 84 % of the total uncertainty in global τ estimates; GPP also contributes significantly (15 %), whereas Cveg contributes only marginally (less than 1 %) to the total uncertainty. The high uncertainty in Csoil is reflected in the large range across state-of-the-art data products, in which full-depth Csoil spans between 3362 and 4792 PgC. The uncertainty is especially high in circumpolar regions with an uncertainty of 50 % and a low spatial correlation between the different datasets (0.2<r<0.5) when compared to other regions (0.6<r<0.8). These uncertainties cast a shadow on current global estimates of τ in circumpolar regions, for which further geographical representativeness and clarification on variations in Csoil with soil depth are needed. Different GPP estimates contribute significantly to the uncertainties of τ mainly in semiarid and arid regions, whereas Cveg causes the uncertainties of τ in the subtropics and tropics. In spite of the large uncertainties, our findings reveal that the latitudinal gradients of τ are consistent across different datasets and soil depths. The current results show a strong ensemble agreement on the negative correlation between τ and temperature along latitude that is stronger in temperate zones (30–60∘ N) than in the subtropical and tropical zones (30∘ S–30∘ N). Additionally, while the strength of the τ–precipitation correlation was dependent on the Csoil data source, the latitudinal gradients also agree among different ensemble members. Overall, and despite the large variation in τ, we identified robust features in the spatial patterns of τ that emerge beyond the differences stemming from the data-driven estimates of Csoil, Cveg and GPP. These robust patterns, and associated uncertainties, can be used to infer τ–climate relationships and for constraining contemporaneous behavior of Earth system models (ESMs), which could contribute to uncertainty reductions in future projections of the carbon cycle–climate feedback. The dataset of τ is openly available at https://doi.org/10.17871/bgitau.201911 (Fan et al., 2019).
Remote sensing data are essential for monitoring the Earth’s surface waters, especially since the amount of publicly available in-situ data is declining. Satellite altimetry provides valuable information on the water levels and variations of lakes, reservoirs and rivers. In combination with satellite imagery, the derived time series allow the monitoring of lake storage changes and river discharge. However, satellite altimetry is limited in terms of its spatial resolution due to its measurement geometry, only providing information in the nadir direction beneath the satellite’s orbit. In a case study in the Mississippi River Basin (MRB), this study investigates the potential and limitations of past and current satellite missions for the monitoring of basin-wide storage changes. For that purpose, an automated target detection is developed and the extracted lake surfaces are merged with the satellites’ tracks. This reveals that the current altimeter configuration misses about 80% of all lakes larger than 0.1 km2 in the MRB and 20% of lakes larger than 10 km2, corresponding to 30% and 7% of the total water area, respectively. Past altimetry configurations perform even more poorly. From the larger water bodies represented by a global hydrology model, at least 91% of targets and 98% of storage changes are captured by the current altimeter configuration. This will improve significantly with the launch of the planned Surface Water and Ocean Topography (SWOT) mission.
Increases in water demand often result in unsustainable water use, leaving insufficient amounts of water for the environment. Therefore, water-saving strategies have been introduced to the environmental policy agenda in many (semi)-arid regions. As many such interventions failed to reach their objectives, a comprehensive tool is needed to assess them. We introduced a constructive framework to assess the proposed strategies by estimating five key components of the water balance in an area: (1) Demand; (2) Availability; (3) Withdrawal; (4) Depletion and (5) Outflow. The framework was applied to assess the Urmia Lake Restoration Program (ULRP) which aimed to increase the basin outflow to the lake to reach 3.1 × 109 m3 yr−1. Results suggested that ULRP could help to increase the Outflow by up to 57%. However, successful implementation of the ULRP was foreseen to be impeded because of three main reasons: (i) decreasing return flows; (ii) increased Depletion; (iii) the impact of climate change. Decreasing return flows and increasing Depletion were expected due to the introduction of technologies that increase irrigation efficiency, while climate change could decrease future water availability by an estimated 3–15%. We suggest that to reach the intervention target, strategies need to focus on reducing water depletion rather than water withdrawals. The framework can be used to comprehensively assess water-saving strategies, particularly in water-stressed basins.
Ice nucleating particle concentrations of the past: insights from a
600 year old Greenland ice core
(2020)
Ice nucleating particles (INPs) affect the microphysics in cloud and precipitation processes. Hence, they modulate the radiative properties of clouds. However, atmospheric INP concentrations of the past are basically unknown. Here, we present INP measurements from an ice core in Greenland, which dates back to the year 1370. In total 135 samples were analyzed with the FRIDGE droplet freezing assay in the temperature range from −14 ◦C to −35 ◦C. The sampling frequency was set to 1 in 10 years from 1370 to 1960. From 1960 to 1990 the frequency was increased to 1 sample per year. Additionally, a number of special events were probed, including volcanic episodes. The typical time coverage of a sample was on the order of a few months. Historical atmospheric INP concentrations were estimated with a conversion factor, which depends on the snow accumulation rate of the ice core, particle dry deposition velocity and the wet scavenging ratio. Typical atmospheric INP concentrations were on the order of 0.1 L -1 at −25 ◦C. The INP variability was found to be about 1 – 2 orders of magnitude. Yet, the short-term variability from samples over a seasonal cycle was considerably lower. INP concentrations were significantly correlated to chemical tracers derived from continuous flow analysis (CFA) and ion chromatography (IC) over a broad range of nucleation temperatures. The highest correlation coefficients were found for the particle concentration (dp > 1.2 µm). The correlation is higher for the seasonal samples, where INP concentrations follow a clear annual pattern, highlighting the importance of the annual dust input in Greenland from East Asian deserts during spring. Scanning electron microscopy (SEM) of single particles retrieved from selected samples found particles of soil origin to be the dominant fraction, verifying the significance of mineral dust particles as INPs. Overall, the concentrations compare reasonably well to present day INP concentrations, albeit they are on the lower side. However, we found that the INP concentration at medium supercooled temperatures differed before and after 1960. Average INP concentrations at −23 ◦C, −24 ◦C, −25 ◦C, −26 ◦C and −28 ◦C were significantly higher (and more variable) in the modern day period, which could indicate a potential anthropogenic impact or some post-coring contamination of the topmost, very porous firn.
Ice-nucleating particle concentrations of the past: insights from a 600-year-old Greenland ice core
(2020)
Ice-nucleating particles (INPs) affect the microphysics in cloud and precipitation processes. Hence, they modulate the radiative properties of clouds. However, atmospheric INP concentrations of the past are basically unknown. Here, we present INP measurements from an ice core in Greenland, which dates back to the year 1370. In total 135 samples were analyzed with the FRIDGE droplet freezing assay in the temperature range from −14 to −35 ∘C. The sampling frequency was set to 1 in 10 years from 1370 to 1960. From 1960 to 1990 the frequency was increased to one sample per year. Additionally, a few special events were probed, including volcanic episodes. The typical time coverage of a sample was on the order of a few months. Historical atmospheric INP concentrations were estimated with a conversion factor, which depends on the snow accumulation rate of the ice core, particle dry deposition velocity, and wet scavenging ratio. Typical atmospheric INP concentrations were on the order of 0.1 L−1 at −25 ∘C. The INP variability was found to be about 1–2 orders of magnitude. Yet, the short-term variability from samples over a seasonal cycle was considerably lower. INP concentrations were significantly correlated to some chemical tracers derived from continuous-flow analysis (CFA) and ion chromatography (IC) over a broad range of nucleation temperatures. The highest correlation coefficients were found for the particle concentration (spherical diameter dp > 1.2 µm). The correlation is higher for a time period of seasonal samples, where INP concentrations follow a clear annual pattern, highlighting the importance of the annual dust input in Greenland from East Asian deserts during spring. Scanning electron microscopy (SEM) analysis of selected samples found mineral dust to be the dominant particle fraction, verifying their significance as INPs. Overall, the concentrations compare reasonably well to present-day INP concentrations, albeit they are on the lower side. However, we found that the INP concentration at medium supercooled temperatures differed before and after 1960. Average INP concentrations at −23, −24, −25, −26, and −28 ∘C were significantly higher (and more variable) in the modern-day period, which could indicate a potential anthropogenic impact, e.g., from land-use change.
In this study, we use simulations from seven global vegetation models to provide the first multi‐model estimate of fire impacts on global tree cover and the carbon cycle under current climate and anthropogenic land use conditions, averaged for the years 2001–2012. Fire globally reduces the tree covered area and vegetation carbon storage by 10%. Regionally, the effects are much stronger, up to 20% for certain latitudinal bands, and 17% in savanna regions. Global fire effects on total carbon storage and carbon turnover times are lower with the effect on gross primary productivity (GPP) close to 0. We find the strongest impacts of fire in savanna regions. Climatic conditions in regions with the highest burned area differ from regions with highest absolute fire impact, which are characterized by higher precipitation. Our estimates of fire‐induced vegetation change are lower than previous studies. We attribute these differences to different definitions of vegetation change and effects of anthropogenic land use, which were not considered in previous studies and decreases the impact of fire on tree cover. Accounting for fires significantly improves the spatial patterns of simulated tree cover, which demonstrates the need to represent fire in dynamic vegetation models. Based upon comparisons between models and observations, process understanding and representation in models, we assess a higher confidence in the fire impact on tree cover and vegetation carbon compared to GPP, total carbon storage and turnover times. We have higher confidence in the spatial patterns compared to the global totals of the simulated fire impact. As we used an ensemble of state‐of‐the‐art fire models, including effects of land use and the ensemble median or mean compares better to observational datasets than any individual model, we consider the here presented results to be the current best estimate of global fire effects on ecosystems.
Aim: Plant life‐forms characterize key morphological strategies that enable large‐scale comparisons of plant communities. This study applies Raunkiær's plant life‐form concept that was developed for temperate climate to a subtropical island flora, in parts, dominated by summer aridity. We quantify how plant life‐form patterns as well as patterns of important plant functional traits (PFTs) relate to important climate and topographic characteristics.
Location: La Palma, Canary Islands.
Taxon: Flora of La Palma.
Methods: We assigned each native plant species a plant life‐form, that is, phanerophyte, chamaephyte, hemicryptophyte, geophyte and therophyte, as well as PFTs (succulence and N‐fixer). We used stacked species distribution models to assess occurrence probability for each species using the Atlantis database (500 m × 500 m grid). We related richness and percentage values for each plant life‐form and PFT to climate and topography.
Results: Plant life‐forms and PFTs showed a clear pattern within geographic but also climate space, while topography had a minor effect. Phanerophytes mainly contributed to the flora in humid areas. Chamaephytes and hemicryptophytes most strongly contributed to the summit scrub flora and, to some degree, also to the arid coastal regions. Geophytes and therophytes were mainly found in dry coastal regions. N‐fixers contributed mainly to warm‐arid and cool‐arid regions, while succulent species were mainly found in arid coastal regions.
Main conclusions: Raunkiær's plant life‐form concept can be comprehensively transferred to a subtropical island flora by adapting to local unfavourable growing conditions, that is, aridity. Using the strong environmental gradients offered by our study island, we identify substantial climate‐driven variation in patterns of plant life‐forms and PFTs that might be used for large‐scale comparisons in macroecological studies. The growth strategies reflected in Raunkiær's plant life‐forms suggest differences in species establishment and coexistence dynamics within different parts of the island's climate space.
The complexity of atmospherical processes has always yielded a multitude of ways of knowing about the weather. What has been lacking in the historiography of meteorology so far is a way to formulate differences between forms of knowledge in a way that does not privilege modern scientific structures, but focuses instead on the epistemological category of causality. Using causality as ground of comparison for different knowledge claims, I shall argue, may enable researchers to investigate meteorological knowledge across time periods, perhaps even geographical regions, in a more symmetrical manner. This review demonstrates this approach as a means to organize a large set of historical meteorological writings from German countries between 1750 and 1850. Three distinct forms of knowledge (Semiotics, Physics, and Organics of the weather) during that time and in that region are suggested and will be described. While a bibliography with a national perspective from the 1880s was the basis for the selection of historical sources, such a setup proved awkward even to contemporaries. In addition, the bibliography came with a number of biases and shortcomings that will be critically reviewed.
We apply seismic full waveform inversion to SH‐ and Love‐wave data for investigating the near‐surface lithology at an archaeological site. We evaluate the resolution of the applied full waveform inversion algorithm through ground truthing in the form of an excavation and sediment core studies. Thereby, we investigate the benefits of full waveform inversion in comparison with other established methods of near‐surface prospecting in terms of resolution capabilities and interpretation security. The study is performed in a presumed harbour area of the ancient Thracian city of Ainos. The exemplary target is the source of a linear magnetic anomaly oriented perpendicular to the coast, which was found in a previous magnetic gradiometry survey, suggesting a mole. The SH‐wave full waveform inversion recovered a subsurface SH‐wave velocity model with submeter resolution showing lateral and vertical velocity variation between 40 and 150 m/s. To tame the non‐linearity of the full waveform inversion, a sequential inversion of frequency bands has to be combined with time‐windowing in order to separate the Love wave from the reflected SH wavefield. We compare the full waveform inversion results with multichannel analysis of surface waves, standard seismic reflection imaging, electrical resistivity tomography and electromagnetic induction. It turns out that the respective depth sections are correlated to a certain degree with the full waveform inversion results. However, the structural resolution of the other geophysical methods is significantly lower than for the full waveform inversion. An exception is the reflection seismic imaging, which shows the same resolution as full waveform inversion but can only be interpreted together with the full waveform inversion–based velocity model. An archaeological excavation as well as coring data allows ground truthing and a direct understanding of the geophysical structures. The results show that the target was a sort of near‐surface trench of about 3–4 m width and 0.8 m to 1.0 m depth, filled with silty sediment, which differs from the layered surrounding in colour and composition. The ground truthing revealed that only SH‐wave full waveform inversion and seismic reflection imaging could image the trench and sediment structure with satisfying lateral and depth resolution. We emphasize that the velocity distribution from SH‐wave full waveform inversion agrees closely with the excavated subsurface structures, and that the discovered changes in seismic velocity correlate with changes in the sand content in the respective sediment facies sequences. The study demonstrated that SH‐wave full waveform inversion is capable to image structural and lithological changes in the near subsurface at scales as low as 0.5 m, thus providing the high resolution needed for archaeological and geoarchaeological prospection.
Opportunities and challenges for paleoaltimetry in "small" orogens: insights from the European Alps
(2020)
Many stable isotope paleoaltimetry studies have focused on paleoelevation reconstructions of orogenic plateaus such as the Tibetan or Andean Plateaus. We address the opportunities and challenges of applying stable isotope paleoaltimetry to “smaller” orogens. We do this using a high‐resolution isotope tracking general circulation model (ECHAM5‐wiso) and explore the precipitation δ18O (δ18Op) signal of Cenozoic paleoclimate and topographic change in the European Alps. Results predict a maximum δ18Op change of 4–5‰ (relative to present day) during topographic development of the Alps. This signal of topographic change has the same magnitude as changes in δ18Op values resulting from Pliocene and Last Glacial Maximum global climatic change. Despite the similar magnitude of the isotopic signals resulting from topographic and paleoclimate changes, their spatial patterns across central Europe differ. Our results suggest that an integration of paleoclimate modeling, multiproxy approaches, and low‐elevation reference proxy records distal from an orogen improve topographic reconstructions.
Global‐scale gradient‐based groundwater models are a new endeavor for hydrologists who wish to improve global hydrological models (GHMs). In particular, the integration of such groundwater models into GHMs improves the simulation of water flows between surface water and groundwater and of capillary rise and thus evapotranspiration. Currently, these models are not able to simulate water table depth adequately over the entire globe. Unsatisfactory model performance compared to well observations suggests that a higher spatial resolution is required to better represent the high spatial variability of land surface and groundwater elevations. In this study, we use New Zealand as a testbed and analyze the impacts of spatial resolution on the results of global groundwater models. Steady‐state hydraulic heads simulated by two versions of the global groundwater model G3M, at spatial resolutions of 5 arc‐minutes (9 km) and 30 arc‐seconds (900 m), are compared with observations from the Canterbury region. The output of three other groundwater models with different spatial resolutions is analyzed as well. Considering the spatial distribution of residuals, general patterns of unsatisfactory model performance remain at the higher resolutions, suggesting that an increase in model resolution alone does not fix problems such as the systematic overestimation of hydraulic head. We conclude that (1) a new understanding of how low‐resolution global groundwater models can be evaluated is required, and (2) merely increasing the spatial resolution of global‐scale groundwater models will not improve the simulation of the global freshwater system.
Barriereinseln und Außensände, geformt durch eine Kombination aus Wind, Wellen, Strömung und Küstenlängstransport gelten als morphologisch hoch aktive Küstenbereiche und variieren häufig in Ursprung, Genese und Entwicklung. Sie besitzen durch ihre dissipative Wirkungsweise eine bedeutende Schutzfunktion für rückwärtige Inseln, Halligen und Festlandbereiche und bilden vor der Küste Nordfrieslands die westliche Außengrenze des Wattenmeeres.
Ziel der nachfolgenden Studie ist es, anhand hochauflösender Georadarmessungen und sedimentologischer Daten aus Bohrungen die Landschaftsgeschichte an der Westküste Amrums sowie die Entstehung und interne sedimentäre Architektur des der Insel westlich angelagerten Kniepsandes zu untersuchen und die Prozesse, die zur Genese geführt haben, durch Datierungen zeitlich einzuordnen. Auf Grundlage der gewonnenen Daten wurden zwei stratigraphische Modelle entwickelt, welche die geologisch-geomorphologischen Prozesse und Sedimentationsbedingungen im Westküstenvorfeld erklären.
Ein Modell zeigt die Landschaftsentwicklung an der Westküste Amrums und beschreibt die Sedimentationsbedingungen, die im Vorfeld der Westküste herrschten, bevor der Kniepsand an die Insel heranwanderte. Auf der Landoberfläche des ertrinkenden saaleeiszeitlichen Geestkerns wurden zu Beginn der Flandrischen Transgression feinkörnige Sedimente eines Mischund Schlickwatts abgelagert. Es ist davon auszugehen, dass der damals noch weit vor der Küste Amrums liegende Kniepsand eine Barriere bildete und so an der heute hochenergetischen Westküste für strömungsberuhigte Sedimentationsbedingungen sorgte. Durch Erosion am Geestkern bildeten sich in unterschiedlichen Höhenpositionen fossile Kliffs, die dem damaligen Meeresspiegelstand entsprechen.
Ein weiteres Modell beschreibt den Andockmechanismus des Kniepsandes an die Insel Amrum. Durch die Anlagerung des ehemaligen Außensandes und den damit einhergehenden Sedimentinput wurden die Bedingungen für eine großflächige Dünenbildung geschaffen.
Als erste Ergebnisse der derzeit laufenden paläomagnetischen Untersuchungen im Rhein-Main-Gebiet werden die Befunde von einigen stratigraphisch bedeutenden Profilen mitgeteilt. Während die Remanenzwerte aus altpleistozänen Tonen der Kelsterbacher Terrasse mit einer Ausnahme revers sind, konnte im Autobahneinschnitt Abenheim (NW Worms) und im Lößprofil der Ziegeleigrube Bad Soden am Taunus der Ubergang von reverser zu normaler Magnetisierung gefunden werden. Die Grenze zwischen MATUYAMA- und BRUNHES-Epoche liegt im Profil Bad Soden unter dem 6., das Jaramillo-Event unter dem 7. fossilen Bt-Horizont.
Die über dem letztinterglazialen Boden (Parabraunerde bzw. pseudovergleyte Parabraunerde) folgenden würmzeitlichen äolischen Sedimente können durch kennzeichnende Bodenhorizonte in 3 Abschnitte (Alt-, Mittel- und Jungwürm) gegliedert werden. Charakteristisch für das Altwürm sind Humuszonen, für das Mittelwürm neben einigen „Naßböden" braune Verlehmungszonen bis zu einer Mächtigkeit von 1,1 m und für das Jungwürm mehrere schwach ausgebildete, geringmächtige Verbraunungszonen und „Naßböden". Den wichtigsten Leithorizont des Jungwürms bildet das Kärlicher Tuffbändchen, das in jüngster Zeit auch in Nordhessen aufgefunden wurde. Abschließend wird das in Hessen aufgestellte Gliederungsschema mit den Würmlöß-Gliederungen in anderen Teilen Europas verglichen und eine Parallelisierung versucht.
Der Auffassung, das Substrat der Lockerbraunerde auf dem Oberwald des Vogelsberges sei ausschließlich ein äolisches Sediment der Jüngeren Tundrenzeit, wird widersprochen. Verschiedene Befunde, vor allem 14C-Datierungen, lassen vielmehr den Schluß zu, daß große Teile des Substrates vieler Lockerbraunerden holozäne anthropogene Kolluvien sind. Bekräftigt wird dagegen der Befund gleichen (jungtundrenzeitlichen) Alters für Lockerbraunerde-Substrat und Deckschutt (Decksediment).
Aus dem Gebiet der weichselzeitlichen Vereisung in Polen werden allgemein verbreitete geringmächtige periglaziale Deckschichten beschrieben. Sie zeigen in der Regel eine äolische Beeinflussung und unterscheiden sich dadurch vom Liegenden. Es handelt sich um spätglaziale Bildungen, wie sie aus dem Jungmoränengebiet der DDR seit langem bekannt sind. Ähnliche Substrate wurden auch im nördlichen Alpenvorland gefunden.
Im Lößprofil der Ziegelei Glos bei Lisieux ist eine Abfolge fossiler Boden aufgeschlossen, die sich gut mit der hessischen Lößstratigraphie verbinden läßt. Über dem letztinterglazialen Boden liegen nach einer Diskordanz der Lohner Boden, der E1-, der E2- und der E4-Naß-boden. Letzterer wird häufig als Äquivalent des „Sol de Kesselt" angesehen. Somit zeigt sich auch für dieses Gebiet, daß dieser Boden stratigraphisch nicht dem Lohner Boden entsprechen kann.
Un tuf volcanique a ete mis en evidence dans trois profils de loess de Hesbaye (Rocourt et Lixhe) et du Limbourg neerlandais (Nagelbeek). Ce tuf correspond probablement ä l'Eltviller Tuff connu en Hesse (SEMMEL 1967). En consequence, cette decouverte permet de comparer directement une partie de la sequence des loess du Pleistocene superieur propre ä chaque region
In der Monte Cavallo-Gruppe fand Verf. in Stauseesedimenten fossile Holzstücke, für die die 14 C-Datierung ein Alter von 29 350 ± 460 Jahren vor 1950 n. Chr. ergab. Die schluffigen Ablagerungen, die in etwa 900 m, fast 80 m über dem heutigen Flußbett des T. Caltea aufgeschlossen sind, ruhen auf Schottern und werden von Moräne überlagert.
Aufgrund seiner stratigraphischen Lage und seines 14C-datierten Alters kann das Holz (Picea abies bzw. Larix) einem zeitlichen Äquivalent des Paudorf-Interstadials zugeordnet werden. Dieser Fund von Großresten ermöglicht somit eine erste absolute Datierung des Paudorf-Interstadials für die Südalpen und Oberitalien.
Das Cranium eines fossilen Hominiden des Formenkreises Homo sapiens sapiens wurde relativ-geologisch sowie absolut durch Radiokohlenstoff und Aminosäuren auf ungefähr 31 000 Jahre B.P. datiert. Andere absolute sowie relative Daten an Mollusken und Mammutzähnen in überlagernden jüngeren Straten datieren auf 18 000 — 21000 und 16 000 Jahren B.P. Geomorphologische und geophysikalische Datierungen stimmen somit gut überein. Er ist der älteste datierte und früheste Bewohner Zentraleuropas, der dem Homo sapiens sapiens angehört.
Aus Franken wird die Entwicklung quartärer Hohlformen beschrieben, deren Rekonstruktion mit Hilfe lößstratigraphischer Methoden (fossile Böden, Tuffbänder, Umlagerungszonen etc.) möglich ist. Bei vielen Formen zeigt sich, daß sie bereits größere Vorläuferformen präwürmzeitlichen Alters hatten. Die Entwicklung während des Würms läßt sich an manchen Beispielen in besonders instruktiver Weise verfolgen. Zu Beginn des Würms, im unteren Mittelwürm und im unteren Jungwürm dominierte zeitweise die Abtragung und Verlagerung. Im oberen Mittelwürm sowie im oberen Jungwürm herrschte äolische Lößsedimentation vor. Diese Ergebnisse stimmen gut mit den bereits aus anderen mitteleuropäischen Lößgebieten bekannten Befunden überein. Mit dem Trockental-System von Helmstadt wird die Entwicklung von Hohlformen beschrieben, deren Anlage bis in das ältere Pleistozän zurückreicht.
In vorliegender Untersuchung wurde der Rißlöß zwischen der 1. und 2. fossilen Parabraunerde anhand schwächerer Bodenbildungen und eingeschalteter Abtragungsphasen zu gliedern versucht. Im jüngeren Riß herrschte starke Lößsedimentation vor, wobei es in mindestens 6 kaltfeuchten Abschnitten zur Ausbildung schwacher periglazialer Naßböden kam. Die Naßbodenserie wurde als Bruchköbeler Böden (B) bezeichnet. Im jüngsten Rißlöß ist wenige dm unter dem Eemboden als tephrochronologischer Leithorizont der Krifteler Tuff (vgl. SEMMEL 1968) eingeschaltet. Den mittleren Profilbereich im Rißlöß zeichnen feuchtere Klimaabschnitte mit starken Verschwemmungsphasen aus, die in den meisten Profilen zu erheblichen Diskordanzen geführt haben. An der Basis der wenigen kompletten Rißlöß-Profile treten vorwiegend in Hessen über der zumeist gekappten 2. fossilen Parabraunerde maximal zwei Schwarzerden auf, die von SEMMEL (1968) als Weilbacher Humuszonen bezeichnet werden. Unmittelbar über diesen Schwarzerden folgt die Ostheimer Zone, eine Fließerde aus aufgearbeitetem Solumaterial der liegenden Böden. Insgesamt zeigt die aus den Rißböden rekonstruierte Klimaabfolge — neben geringfügigen Abweichungen — überraschende Parallelen zur paläopedologisch-klimatischen Gliederung der Würmkaltzeit.
Crustal structure at the western end of the North Anatolian Fault Zone from deep seismic sounding
(2001)
The first deep seismic sounding experiment in Northwestern Anatolia was carried out in October 1991 as part of the "German - Turkish Project on Earthquake Prediction Research" in the Mudurnu area of the North Anatolian Fault Zone. The experiment was a joint enterprise by the Institute of Meteorology and Geophysics of Frankfurt University, the Earthquake Research Institute (ERI) in Ankara, and the Turkish Oil Company (TPAO). Two orthogonal profiles, each 120 km in length with a crossing point near Akyazi, were covered in succession by 30 short period tape recording seismograph stations with 2 km station spacing. 12 shots, with charge sizes between 100 and 250 kg, were fired and 342 seismograms out of 360 were used for evaluation. By coincidence an M b = 4.5 earthquake located below Imroz Island was also recorded and provided additional information on Moho and the sub-Moho velocity. A ray tracing method orginally developed by Weber (1986) was used for travel time inversion. From a compilation of all data two generalized crustal models were derived, one with velocity gradients within the layers and one with constant layer velocities. The latter consists of a sediment cover of about 2 km with V p » 3.6 km/s, an upper crystalline crust down to 13 km with V p » 5.9 km/s, a middle crust down to 25 km depth with V p » 6.5 km/s, a lower crust down to 39 km Moho depth with V p » 7.0 km/s and V p » 8.05 km/s below the Moho. The structure of the individual profiles differs slightly. The thickest sediment cover is reached in the Izmit-Sapanca-trough and in the Akyazi basin. Of particular interest is a step of about 4 km in the lower crust near Lake Sapanca and probably an even larger one in the Moho (derived from the Imroz earthquake data). After the catastrophic earthquake of Izmit on 17 August 1999 this significant heterogeneity in crustal structure appears in a new light with regard to the possible cause of the Izmit earthquake. Heterogeneities in structure are frequently also heterogeneities in strength and stress that impede or even lock rupture. The Izmit earthquake is discussed in relation to a large stepover or jog at the North Anatolian Fault.
This study reports and discusses new radiometric ages, petrographical and volcanological observations and whole rock geochemical data of the rocks of the Rudnik Mts. volcano-intrusive complex. The complex hosts a Pb-Zn-Ag deposit and belongs to the Serbo-Macedonian metallogenetic belt. Two distinct igneous events are distignuished. The first occurred >30 Ma and was mainly characterized by extrusive and shallow intrusive dacites and andesites and was unrelated to mineralization. The second igneous event occurred <23 Ma and was highly heterogeneous in terms of volcanic products and petrographic varieties, but with predominance of quartzlatites. The dacite-andesites (first event) and the quartzlatites (second event) are geochemically similar and display a calc-alkaline affinity and highly incompatible element enriched patterns on spider diagrams, but the younger quartzlatites are richer in K2O, Rb and Ba and poorer in Sr. This is taken as evidence that mixing between an ultrapotassic lamprophyre/lamproite magma and an acid calc-alkaline (dacite-like) magma was essential petrogenetic processes during the second event. The proposed simplified volcanological model suggests that this mixing was responsible for triggering strongly explosive volcanic activity as well as for providing conditions for active hydrothermal and mineralization processes. The observed link between a specific magmatic phase and ore deposit formation can be a general phennomenon in the Balkans, and must be addressed by further and more advanced studies.
This paper provides new geochemical data focusing on valuable elements in the coal, parting, and floor samples in the No. 5 coal seam of the Taiyuan Formation from the Wujiawan mine, Datong coalfield, northern China. The minerals mainly consist of kaolinite, calcite, and pyrite, as well as trace amounts of quartz and illite. The No. 5 coal is enriched in Li, Ga, high field strength elements (HFSEs), and rare earth elements and yttrium (REY) when compared with world hard coals. Of particular interest is the high average concentration of Li (67.66 μg/g), which is around seven times higher than the value for world hard coals. Lithium, Ga, and HFSEs have strong inorganic affinities, whereas REY have organic affinities. The main carrier of Li, Ga, and HFSEs is aluminosilicate minerals, while REY appear to occur with organophosphorus. These HFSEs are enriched, both in the parting and in the adjacent coal samples. This suggests that these elements are likely to leach out during the diagenetic process. The distribution patterns of REY, along with the ratio of Al2O3/TiO2 and the figure of Zr/TiO2 vs. Nb/Y are suggestive of their derivation from felsic parent material. In the northern and eastern part of the Datong coalfield, there are several regions where the Li content is higher than the mineable grade, in particular in the northern Datong coalfield where there is a mine with an Li content of 294.6 μg/g. This is significantly higher than the mineable grade. Therefore, there is a potential for financially viable recovery of Li in these coals of the Datong coalfield.
Tropical cyclones (TC) represent a substantial threat to life and property for Caribbean and adjacent populations. The prospective increase of TC magnitudes, expressed in the 15th chapter of the IPCC AR5 report, entails a rising probability of ecological and social disasters, which were tragically exemplified by several severe Caribbean TC strikes during the past 20 years. Modern IPCC-grade climate models, however, still lack the required spatial and temporal resolution to accurately consider the underlying boundary conditions that modulate long-time TC patterns beyond the Instrumental Era. It is thus necessary to provide a synoptic mechanistic understanding regarding the origin of such long-time patterns, in order to predict reliable changes of TC magnitude and frequency under future climate scenarios. Caribbean TC records are still rare and often lack the necessary continuity and resolution to overcome these limitations. Here, we report on an annually-resolved sedimentary archive from the bottom of the Great Blue Hole (Lighthouse Reef, Belize). The TC record encompasses 1885 years and extends all existing site-specific TC archives both in terms of resolution and duration. We identified a likely connection between long-term TC patterns and climate phenomena responses to Common Era climate variations and offer a conceptual and comparative view considering several involved tropospheric and oceanographic control mechanisms such as the El-Niño-Southern-Oscillation, the North Atlantic Oscillation and the Atlantic Multidecadal Oscillation. These basin-scaled climate modes exercise internal control on TC activity by modulating the thermodynamic environment (sea-surface temperature and vertical wind shear stress dynamics) for enhanced/suppressed TC formation both on millennial (primary) and multi-decadal (secondary) time scales. We interpret the beginning of the Medieval Warm Period (MWP) as an important time interval of the Common Era record and suspect that the southward migration of the intertropical convergence zone (ITCZ) caused, in combination with extensive hydro-climate changes, a shift in the tropical Atlantic TC regime. The TC activity in the south-western Caribbean changed in general from a stable and less active stage (100–900 CE) to a more active and variable state (1,100 CE-modern).
The Early Permian coal is of great value in the Tengxian Coalfield, Shandon Province, Eastern China. This work deals with the new data focusing on mineralogical characteristics in the Early Permian Shanxi Formation No. 3 coal from the Jinyuan Mine. The Jinyuan coal is a low ash and highly volatile A bituminous coal. Minerals in the No. 3 coal mainly comprise of kaolinite, ankerite, illite, calcite, siderite, and quartz, with varying compositions of trace amounts of pyrite, jarosite, bassanite, anatase, and rutile. According to mineral assemblage in the coal plies, three Types (A to C) can be identified in the No. 3 coal. The dominant minerals in Type A are poorly-ordered kaolinite, illite, quartz, pyrite, and jarosite. Type B is mainly composed of well-ordered kaolinite, illite, siderite, ankerite, and calcite. Type C, with just one sample (JY-3-7c), which contains high proportions of calcite (54%) and ankerite (34%). Terrigenous minerals are elevated in coal plies that typically have relatively high contents of ash yield. The formation of syngenetic pyrite was generally due to seawater, while the sulphate minerals (jarosite and coquimbite) were derived from the oxidation of pyrite. Epigenetic vein-like or fracture-fillings carbonate minerals (ankerite, calcite, and siderite), kaolinite, and pyrite, as well as authigenic quartz were derived from the influx of hydrothermal fluids during different periods, from the authigenic to epigenetic. The paragonite in the coal may have been formed by the precipitated from Na-rich hydrothermal fluids. No effects of magmatic intrusion on mineralogy were investigated in this research.
We present novel measurements of five short-lived brominated source gases (CH2Br2, CHBr3, CH2ClBr, CHCl2Br and CHClBr2) obtained using a gas chromatograph-mass spectrometer system on board the High Altitude and Long Range Research Aircraft (HALO). The instrument is extremely sensitive due to the use of chemical ionisation, allowing detection limits in the lower parts per quadrillion (10-15) range. Data from three campaigns using the HALO aircraft are presented, where the Upper Troposphere/Lower Stratosphere (UTLS) of the Northern Hemisphere mid to high latitudes were sampled during winter and during late summer to early fall. We show that an observed decrease with altitude in the stratosphere is consistent with the relative lifetimes of the different compounds. Distributions of the five source gases and total organic bromine just below the tropopause shows an increase in mixing ratio with latitude, in particular during polar winter. This increase in mixing ratio is explained by increasing lifetimes at higher latitudes during winter. As the mixing ratio at the extratropical tropopause are generally higher than those derived for the tropical tropopause, extratropical troposphere-to-stratosphere transport will result in elevated levels of organic bromine in comparison to air transported over the tropical tropopause. The observations are compared to model estimates using different emission scenarios. A scenario which has emissions most strongly concentrated to low latitudes cannot reproduce the observed latitudinal distributions and will tend to overestimate bromine input through the tropical tropopause from CH2Br2 and CHBr3. Consequently, the scenario also overestimates the amount of brominated organic gases in the stratosphere. The two scenarios with the highest overall emissions of CH2Br2 tend to overestimate mixing ratios at the tropical tropopause but are in much better agreement with extratropical tropopause values, showing that not only total emissions but also latitudinal distributions in the emissions are of importance. While an increase in tropopause values with latitude is reproduced with all emission scenarios during winter, the simulated extratropical tropopause values are on average lower than the observations during late summer to fall. We show that a good knowledge of the latitudinal distribution of tropopause mixing ratios and of the fractional contributions of tropical and extratropical air is needed to derive stratospheric inorganic bromine in the lowermost stratosphere from observations. Depending on the underlying emission scenario, differences of a factor 2 in reactive bromine derived from observations and model outputs are found for the lowermost stratosphere, based on source gas injection. We conclude that a good representation of the contributions of different source regions is required in models for a robust assessment of the role of short-lived halogen source gases on ozone depletion in the UTLS.
Background: Point of care devices for performing targeted coagulation substitution in bleeding patients have become increasingly important in recent years. New on the market is the Quantra® from HemoSonics (LC, Charlottesville, VA, US). It uses sonorheometry, a sonic estimation of elasticity via resonance (SEER), a novel ultrasound-based technology that measures viscoelastic properties of whole blood. Several studies have already shown the comparability with devices already established on the market such as the ROTEM® (TEM International GmbH, Munich, Germany).
Objective: In contrast to existing studies, the planned study will be the first prospective interventional study using the new Quantra® system in a cardiac surgical patient cohort. The aim is to investigate the non-inferiority between an already existing coagulation algorithm, based on ROTEM®/Multiplate®, and a new algorithm based on the Quantra®, for the treatment of coagulopathic cardiac surgical patients.
Methods: The study is divided into two phases. In an initial observation phase, whole blood samples of 20 patients will be analyzed using both ROTEM®/Multiplate® and Quantra® obtained at three defined points of time (prior to surgery, after completion of cardiopulmonary bypass, on arrival in the intensive care unit). The obtained threshold values will be used to create an algorithm for hemotherapy. In a second intervention phase, the new algorithm will be tested against an algorithm used routineously for years at our department for non-inferiority.
Results: The main objective of the examination is the cumulative loss of blood within 24 hours after surgery. Statistical calculations based on literature and in-house data suggest that the new algorithm is not inferior if the difference in cumulative blood loss is < 150ml/24 h.
Conclusions: Because of the comparability of the Quantra® sonorheometry system with ROTEM® rotational thromboelastometric measurement methods, the existing hemotherapy treatment algorithm can be adapted to the Quantra device with a proof of non-inferiority. Clinical Trial: International Registered Report Identifier (IRRID): clinicaltrials.gov: NCT03902275
Inclusions of breyite (previously known as walstromite-structured CaSiO3) in diamond are usually interpreted as retrogressed CaSiO3 perovskite trapped in the transition zone or the lower mantle. However, the thermodynamic stability field of breyite does not preclude its crystallization together with diamond under upper-mantle conditions (6–10 GPa). The possibility of breyite forming in subducted sedimentary material through the reaction CaCO3 + SiO2 = CaSiO3 + C + O2 was experimentally evaluated in the CaO–SiO2–C–O2 ± H2O system at 6–10 GPa, 900–1500 ∘C and oxygen fugacity 0.5–1.0 log units below the Fe–FeO (IW) buffer. One experimental series was conducted in the anhydrous subsystem and aimed at determining the melting temperature of the aragonite–coesite (or stishovite) assemblage. It was found that melting occurs at a lower temperature (∼1500 ∘C) than the decarbonation reaction, which indicates that breyite cannot be formed from aragonite and silica under anhydrous conditions and an oxygen fugacity above IW – 1. In the second experimental series, we investigated partial melting of an aragonite–coesite mixture under hydrous conditions at the same pressures and redox conditions. The melting temperature in the presence of water decreased strongly (to 900–1200 ∘C), and the melt had a hydrous silicate composition. The reduction of melt resulted in graphite crystallization in equilibrium with titanite-structured CaSi2O5 and breyite at ∼1000 ∘C. The maximum pressure of possible breyite formation is limited by the reaction CaSiO3 + SiO2 = CaSi2O5 at ∼8 GPa. Based on the experimental results, it is concluded that breyite inclusions found in natural diamond may be formed from an aragonite–coesite assemblage or carbonate melt at 6–8 GPa via reduction at high water activity.
The accurate knowledge of the groundwater storage variation (ΔGWS) is essential for reliable water resource assessment, particularly in arid and semi-arid environments (e.g., Australia, the North China Plain (NCP)) where water storage is significantly affected by human activities and spatiotemporal climate variations. The large-scale ΔGWS can be simulated from a land surface model (LSM), but the high model uncertainty is a major drawback that reduces the reliability of the estimates. The evaluation of the model estimate is then very important to assess its accuracy. To improve the model performance, the terrestrial water storage variation derived from the Gravity Recovery And Climate Experiment (GRACE) satellite mission is commonly assimilated into LSMs to enhance the accuracy of the ΔGWS estimate. This study assimilates GRACE data into the PCRaster Global Water Balance (PCR-GLOBWB) model. The GRACE data assimilation (DA) is developed based on the three-dimensional ensemble Kalman smoother (EnKS 3D), which considers the statistical correlation of all extents (spatial, temporal, vertical) in the DA process. The ΔGWS estimates from GRACE DA and four LSM simulations (PCR-GLOBWB, the Community Atmosphere Biosphere Land Exchange (CABLE), the Water Global Assessment and Prognosis Global Hydrology Model (WGHM), and World-Wide Water (W3)) are validated against the in situ groundwater data. The evaluation is conducted in terms of temporal correlation, seasonality, long-term trend, and detection of groundwater depletion. The GRACE DA estimate shows a significant improvement in all measures, notably the correlation coefficients (respect to the in situ data) are always higher than the values obtained from model simulations alone (e.g., ~0.15 greater in Australia, and ~0.1 greater in the NCP). GRACE DA also improves the estimation of groundwater depletion that the models cannot accurately capture due to the incorrect information of the groundwater demand (in, e.g., PCR-GLOBWB, WGHM) or the unavailability of a groundwater consumption routine (in, e.g., CABLE, W3). In addition, this study conducts the inter-comparison between four model simulations and reveals that PCR-GLOBWB and CABLE provide a more accurate ΔGWS estimate in Australia (subject to the calibrated parameter) while PCR-GLOBWB and WGHM are more accurate in the NCP (subject to the inclusion of anthropogenic factors). The analysis can be used to declare the status of the ΔGWS estimate, as well as itemize the possible improvements of the future model development.
Stability of halocarbons in whole air samples from the upper
troposphere and lowermost stratosphere
(2019)
Measurements of halogenated hydrocarbons of ambient air frequently rely on canister sampling followed by offline laboratory analysis. This allows for a large number of compounds to be analysed under stable conditions, maximising measurement precision. However, individual compounds might be affected during storage of canister samples. In order to assess halocarbon stability in whole air samples from the upper troposphere and lowermost stratosphere, we performed stability tests using the air sampling unit High REsolution Sampler (HIRES) which is part of the CARIBIC (Civil Aircraft for the Regular Investigation of the Atmosphere Based on an Instrument Container) instrument package. HIRES holds 88 light-weight stainless steel cylinders that are pressurized in flight to 4.5 bar using metal bellows pumps. The HIRES sampling unit was first deployed in 2010, but has up to now not been used for regular halocarbon analysis with exception of chloromethane. The sample collection unit was tested for sampling and storage effects of 28 halogenated compounds. The focus was on compound stability in the stainless steel canisters during storage of up to five weeks and on the influence of ozone, since flights take place in the upper troposphere and lowermost stratosphere with ozone mixing ratios of up to several hundred ppbV. Most of the investigated (hydro)chlorofluorocarbons and long-lived hydrofluorocarbons were found to be stable over a storage time of up to five weeks and were unaltered by ozone being present during pressurization. Some compounds such as for example dichloromethane, trichloromethane and tetrachloroethene started to decrease in the canisters after a storage time of more than two weeks or exhibited lowered mixing ratios in samples pressurized with ozone present. Few compounds such as for example tetrachloromethane and tribromomethane were found to be not stable in the HIRES stainless-steel canisters independent of ozone levels. Also growth was observed during storage, namely for HFC-152a and HFC-23.
Measurements of halogenated trace gases in ambient air frequently rely on canister sampling followed by offline laboratory analysis. This allows for a large number of compounds to be analysed under stable conditions, maximizing measurement precision. However, individual compounds might be affected during the sampling and storage of canister samples. In order to assess halocarbon stability in whole-air samples from the upper troposphere and lowermost stratosphere, we performed stability tests using the high-resolution sampler (HIRES) air sampling unit, which is part of the Civil Aircraft for the Regular Investigation of the atmosphere Based on an Instrument Container (CARIBIC) instrument package. The HIRES unit holds 88 lightweight stainless-steel cylinders that are pressurized in flight to 4.5 bar using metal bellows pumps. The HIRES unit was first deployed in 2010 but has up to now not been used for regular halocarbon analysis with the exception of chloromethane analysis. The sample collection unit was tested for the sampling and storage effects of 28 halogenated compounds. The focus was on compound stability in the stainless-steel canisters during storage of up to 5 weeks and on the influence of ozone, since flights take place in the upper troposphere and lowermost stratosphere with ozone mixing ratios of up to several hundred parts per billion by volume (ppbv). Most of the investigated (hydro)chlorofluorocarbons and long-lived hydrofluorocarbons were found to be stable over a storage time of up to 5 weeks and were unaltered by ozone being present during pressurization. Some compounds such as dichloromethane, trichloromethane, and tetrachloroethene started to decrease in the canisters after a storage time of more than 2 weeks or exhibited lowered mixing ratios in samples pressurized with ozone present. A few compounds such as tetrachloromethane and tribromomethane were found to be unstable in the HIRES stainless-steel canisters independent of ozone levels. Furthermore, growth was observed during storage for some species, namely for HFC-152a, HFC-23, and Halon 1301.
A list of authors and their affiliations appears at the end of the paper New-particle formation is a major contributor to urban smog, but how it occurs in cities is often puzzling. If the growth rates of urban particles are similar to those found in cleaner environments (1–10 nanometres per hour), then existing understanding suggests that new urban particles should be rapidly scavenged by the high concentration of pre-existing particles. Here we show, through experiments performed under atmospheric conditions in the CLOUD chamber at CERN, that below about +5 degrees Celsius, nitric acid and ammonia vapours can condense onto freshly nucleated particles as small as a few nanometres in diameter. Moreover, when it is cold enough (below −15 degrees Celsius), nitric acid and ammonia can nucleate directly through an acid–base stabilization mechanism to form ammonium nitrate particles. Given that these vapours are often one thousand times more abundant than sulfuric acid, the resulting particle growth rates can be extremely high, reaching well above 100 nanometres per hour. However, these high growth rates require the gas-particle ammonium nitrate system to be out of equilibrium in order to sustain gas-phase supersaturations. In view of the strong temperature dependence that we measure for the gas-phase supersaturations, we expect such transient conditions to occur in inhomogeneous urban settings, especially in wintertime, driven by vertical mixing and by strong local sources such as traffic. Even though rapid growth from nitric acid and ammonia condensation may last for only a few minutes, it is nonetheless fast enough to shepherd freshly nucleated particles through the smallest size range where they are most vulnerable to scavenging loss, thus greatly increasing their survival probability. We also expect nitric acid and ammonia nucleation and rapid growth to be important in the relatively clean and cold upper free troposphere, where ammonia can be convected from the continental boundary layer and nitric acid is abundant from electrical storms.
Diamond formation in the Earth has been extensively discussed in recent years on the basis of geochemical analysis of natural materials, high-pressure experimental studies, or theoretical aspects. Here, we demonstrate experimentally for the first time, the spontaneous crystallization of diamond from CH4-rich fluids at pressure, temperature and redox conditions approximating those of the deeper parts of the cratonic lithospheric mantle (5-7 GPa) without using diamond seed crystals or carbides. In these experiments the fluid phase is nearly pure methane, even though the oxygen fugacity was significantly above metal saturation. We propose several previously unidentified mechanisms that may promote diamond formation under such conditions and which may also have implications for the origin of sublithospheric diamonds. These include the hydroxylation of silicate minerals like olivine and pyroxene, H2 incorporation into these phases and the "etching" of graphite by H2 and CH4 and reprecipitation as diamond. This study also serves as a demonstration of our new high-pressure experimental technique for obtaining reduced fluids, which is not only relevant for diamond synthesis, but also for investigating the metasomatic origins of diamond in the upper mantle, which has further implications for the deep carbon cycle.
The frequency of extreme events has changed, having a direct impact on human lives. Regional climate models help us to predict these regional climate changes. This work presents an atmosphere–ocean coupled regional climate system model (RCSM; with the atmospheric component COSMO-CLM and the ocean component NEMO) over the European domain, including three marginal seas: the Mediterranean, North, and Baltic Sea. To test the model, we evaluate a simulation of more than 100 years (1900–2009) with a spatial grid resolution of about 25 km. The simulation was nested into a coupled global simulation with the model MPI-ESM in a low-resolution configuration, whose ocean temperature and salinity were nudged to the ocean–ice component of the MPI-ESM forced with the NOAA 20th Century Reanalysis (20CR). The evaluation shows the robustness of the RCSM and discusses the added value by the coupled marginal seas over an atmosphere-only simulation. The coupled system is stable for the complete 20th century and provides a better representation of extreme temperatures compared to the atmosphere-only model. The produced long-term dataset will help us to better understand the processes leading to meteorological and climate extremes.
The most frequently used boundary-layer turbulence parameterization in numerical weather prediction (NWP) models are turbulence kinetic energy (TKE) based-based schemes. However, these parameterizations suffer from a potential weakness, namely the strong dependence on an ad-hoc quantity, the so-called turbulence length scale. The physical interpretation of the turbulence length scale is difficult and hence it cannot be directly related to measurements or large eddy simulation (LES) data. Consequently, formulations for the turbulence length scale in basically all TKE schemes are based on simplified assumptions and are model-dependent. A good reference for the independent evaluation of the turbulence length scale expression for NWP modeling is missing. Here we propose a new turbulence length scale diagnostic which can be used in the gray zone of turbulence without modifying the underlying TKE turbulence scheme. The new diagnostic is based on the TKE budget: The core idea is to encapsulate the sum of the molecular dissipation and the cross-scale TKE transfer into an effective dissipation, and associate it with the new turbulence length scale. This effective dissipation can then be calculated as a residuum in the TKE budget equation (for horizontal sub-domains of different sizes) using LES data. Estimation of the scale dependence of the diagnosed turbulence length scale using this novel method is presented for several idealized cases.
Metasomatic evolution of coesite-bearing diamondiferous eclogite from the Udachnaya Kimberlite
(2020)
A coesite-bearing diamondiferous eclogite from the Udachnaya kimberlite (Daldyn field, Siberian craton) has been studied to trace its complex evolution recorded in rock-forming and minor mineral constituents. The eclogite sample is composed of rock-forming omphacite (60 vol%), garnet (35 vol%) and quartz/coesite (5 vol%) and contains intergranular euhedral zoned olivine crystals, up to 200 µm long, coexisting with phlogopite, orthopyroxene, clinopyroxene (secondary), K-feldspar, plagioclase, spinel, sodalite and djerfisherite. Garnet grains are zoned, with a relatively homogeneous core and a more magnesian overgrowth rim. The rim zones further differ from the core in having higher Zr/Y (6 times that in the cores), ascribed to interaction with, or precipitation from, a kimberlite-related melt. Judging by pressure-temperature estimates (~1200 °C; 6.2 GPa), the xenolith originated at depths of ~180–200 km at the base of the continental lithosphere. The spatial coexistence of olivine, orthopyroxene and coesite/quartz with K-Na-Cl minerals in the xenolith indicates that eclogite reacted with a deep-seated kimberlite melt. However, Fe-rich olivine, orthopyroxene and low-pressure minerals (sodalite and djerfisherite) likely result from metasomatic reaction at shallower depths during transport of the eclogite by the erupting kimberlite melt. Our results demonstrate that a mixed eclogitic-peridotitic paragenesis, reported previously from inclusions in diamond, can form by interaction of eclogite and a kimberlite-related melt.
Inappropriate land management leads to soil loss with destruction of the land’s resource and sediment input into the receiving river. Part of the sediment budget of a catchment is the estimation of soil loss. In the Ruzizi catchment in the Eastern Democratic Republic of the Congo (DRC), only limited research has been conducted on soil loss mainly dealing with local observations on geomorphological forms or river load measurements; a regional quantification of soil loss is missing so far. Such quantifications can be calculated using the Universal Soil Loss Equation (USLE). It is composed of four factors: precipitation (R), soil (K), topography (LS), and vegetation cover (C). The factors can be calculated in different ways according to the characteristics of the study area. In this paper, different approaches for calculating the single factors are reviewed and validated with field work in two sub-catchments of Ruzizi River supplying the water for the reservoirs of Ruzizi I and II hydroelectric dams. It became obvious that the (R)USLE model provides the best results with revised R and LS factors. C factor calculations required to conduct a supervised classification using the Maximum Likelihood Procedure. Different C factor values were assigned to the land cover classes. The calculations resulted in a soil loss rate for the predominantly occurring Ferralsols and Leptosols of around 576 kt/yr in both catchments, when 2016 landcover and precipitation are used. This represents an area-normalized value of 40.4 t/ha/yr for Ruzizi I and 50.5 t/ha/yr for Ruzizi II due to different landcover in the two sub-catchments. The mean value for the whole study area is 47.8 t/ha/yr or even 27.1 t/ha/yr when considering land management techniques like terracing on the slopes (P factor). This work has shown that the (R)USLE model can serve as an easy to handle tool for soil loss quantification when comprehensive field work results are sparse. The model can be implemented in Geographic Information Systems (GIS) with free data; hence, a validation is crucial. It becomes apparent that the use of high resolution Sentinel 2a MSI data as the basis for C factor calculations is an appropriate method for considering heterogeneous Land Use Land Cover (LULC) patterns. To transfer the approach to other regions, the calculation of factor R needs to be modified
Accretion rates of Holocene tropical coral reefs in three areas in the Atlantic, Pacific and Indian Oceans have been quantified in 79 dated core sections in 34 reef cores from Belize, the Maldives and French Polynesia. Holocene vertical reef accretion rate averages 5.05 m/kyr and has decreased during the past 10 kyr. Accretion rates in branched and massive coral facies are statistically similar. Reef accretion rate is positively correlated with the rate of sea‐level rise, that is the degree of creation of accommodation space, and with climate as expressed in a Holocene sea surface temperature anomaly. Accommodation space is also created by subsidence, but at a rate one to two orders of magnitude lower than that created by glacio‐eustasy (0.04 to 0.16 m/kyr). Lagoonal background sedimentation in adjacent reef lagoons averages 0.89 m/kyr as measured in 72 dated core sections in 28 cores. Lagoonal carbonate sedimentation on top of underlying mangrove peat usually starts after a considerable hiatus of ca 3 kyr on average. The lagoonal background sedimentation rate increased during the Holocene, probably due to deepening. The differences between vertical reef accretion and lagoonal background sedimentation rates are a major factor in the production of the widely known saucer shapes typical of tropical reefs and carbonate platforms, that is the creation of unfilled accommodation space. Reef core recovery, used as a proxy for reef consolidation, and core depth exhibit a statistically negative correlation based on data from 326 core barrels. Recovery and marine cement abundance (average volume 8.6%) also decrease from windward to leeward core positions. These observations are presumably a result of both a decrease in the rate of sea‐level rise that is the increase in time available for submarine cementation during the Holocene and the amount of flushing of reef interstices by marine waters.
Diurnal valley winds frequently form over complex topography, particularly under fair weather conditions, and have a significant impact on the local weather and climate. Since diurnal valley winds result from complex and multi-scale interactions, their representation in numerical weather prediction models is challenging. Better understanding of these local winds based on observations is crucial to improve the accuracy of the forecasts. This study investigates the diurnal evolution of the three-dimensional mean wind structure in a deep Alpine valley, the Rhone valley at Sion, using data from a radar wind profiler and a surface weather station operated continuously from 1 September 2016 to 17 July 2017. In particular, the wind profiler data was analyzed for a subset of days on which fair weather conditions allowed for the full development of thermally driven winds. A pronounced diurnal cycle of the wind speed, as well as a reversal of the wind direction twice per day is documented for altitudes up to about 2 km above ground level (AGL) in the warm season and less than 1 km AGL in winter. The diurnal pattern undergoes significant changes during the course of the year. Particularly during the warm-weather months of May through to September, a low-level wind maximum occurs, where mean maximum up-valley velocities of 8–10 m s−1 are found between 15–16 UTC at altitudes around 200 m AGL. In addition, during nighttime, a down-valley jet with maximum wind speeds of 4–8 m s−1 around 1 km AGL is found. A case study of a three-day period in September 2016 illustrates the occurrence of an elevated layer of cross-valley flow around 1–1.5 km AGL.
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.
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.
Convective shower characteristics simulated with the convection-permitting climate model COSMO-CLM
(2019)
This paper evaluates convective precipitation as simulated by the convection-permitting climate model (CPM) Consortium for Small-Scale Modeling in climate mode (COSMO-CLM) (with 2.8 km grid-spacing) over Germany in the period 2001–2015. Characteristics of simulated convective precipitation objects like lifetime, area, mean intensity, and total precipitation are compared to characteristics observed by weather radar. For this purpose, a tracking algorithm was applied to simulated and observed precipitation with 5-min temporal resolution. The total amount of convective precipitation is well simulated, with a small overestimation of 2%. However, the simulation underestimates convective activity, represented by the number of convective objects, by 33%. This underestimation is especially pronounced in the lowlands of Northern Germany, whereas the simulation matches observations well in the mountainous areas of Southern Germany. The underestimation of activity is compensated by an overestimation of the simulated lifetime of convective objects. The observed mean intensity, maximum intensity, and area of precipitation objects increase with their lifetime showing the spectrum of convective storms ranging from short-living single-cell storms to long-living organized convection like supercells or squall lines. The CPM is capable of reproducing the lifetime dependence of these characteristics but shows a weaker increase in mean intensity with lifetime resulting in an especially pronounced underestimation (up to 25%) of mean precipitation intensity of long-living, extreme events. This limitation of the CPM is not identifiable by classical evaluation techniques using rain gauges. The simulation can reproduce the general increase of the highest percentiles of cell area, total precipitation, and mean intensity with temperature but fails to reproduce the increase of lifetime. The scaling rates of mean intensity and total precipitation resemble observed rates only in parts of the temperature range. The results suggest that the evaluation of coarse-grained (e.g., hourly) precipitation fields is insufficient for revealing challenges in convection-permitting simulations.
The sources and critical enrichment processes for granite related tin ores are still not well understood. The Erzgebirge represents one of the classical regions for tin mineralization. We investigated the four largest plutons from the Western Erzgebirge (Germany) for the geochemistry of bulk rocks and autocrystic zircons and relate this information to their intrusion ages. The source rocks of the Variscan granites were identified as high-grade metamorphic rocks based on the comparison of Hf-O isotope data on zircons, the abundance of xenocrystic zircon ages as well as Nd and Hf model ages. Among these rocks, restite is the most likely candidate for later Variscan melts. Based on the evolution with time, we could reconstruct enrichment factors for tin and tungsten starting from the protoliths (575 Ma) that were later converted to high-grade metamorphic rocks (340 Ma) and served as sources for the older biotite granites (323–318 Ma) and the tin granites (315–314 Ma). This evolution involved a continuous enrichment of both tin and tungsten with an enrichment factor of ~15 for tin and ~7 for tungsten compared to the upper continental crust (UCC). Ore level concentrations (>10–100 times enrichment) were achieved only in the greisen bodies and dykes by subsequent hydrothermal processes.
Video and image data are regularly used in the field of benthic ecology to document biodiversity. However, their use is subject to a number of challenges, principally the identification of taxa within the images without associated physical specimens. The challenge of applying traditional taxonomic keys to the identification of fauna from images has led to the development of personal, group, or institution level reference image catalogues of operational taxonomic units (OTUs) or morphospecies. Lack of standardisation among these reference catalogues has led to problems with observer bias and the inability to combine datasets across studies. In addition, lack of a common reference standard is stifling efforts in the application of artificial intelligence to taxon identification. Using the North Atlantic deep sea as a case study, we propose a database structure to facilitate standardisation of morphospecies image catalogues between research groups and support future use in multiple front-end applications. We also propose a framework for coordination of international efforts to develop reference guides for the identification of marine species from images. The proposed structure maps to the Darwin Core standard to allow integration with existing databases. We suggest a management framework where high-level taxonomic groups are curated by a regional team, consisting of both end users and taxonomic experts. We identify a mechanism by which overall quality of data within a common reference guide could be raised over the next decade. Finally, we discuss the role of a common reference standard in advancing marine ecology and supporting sustainable use of this ecosystem.
Here we present a comprehensive attempt to correlate aragonitic Na∕Ca ratios from Desmophyllum pertusum (formerly known as Lophelia pertusa), Madrepora oculata and a caryophylliid cold-water coral (CWC) species with different seawater parameters such as temperature, salinity and pH. Living CWC specimens were collected from 16 different locations and analyzed for their Na∕Ca ratios using solution-based inductively coupled plasma-optical emission spectrometry (ICP-OES) measurements.
The results reveal no apparent correlation with salinity (30.1–40.57 g kg−1) but a significant inverse correlation with temperature (−0.31±0.04 mmolmol−1∘C−1). Other marine aragonitic organisms such as Mytilus edulis (inner aragonitic shell portion) and Porites sp. exhibit similar results highlighting the consistency of the calculated CWC regressions. Corresponding Na∕Mg ratios show a similar temperature sensitivity to Na∕Ca ratios, but the combination of two ratios appears to reduce the impact of vital effects and domain-dependent geochemical variation. The high degree of scatter and elemental heterogeneities between the different skeletal features in both Na∕Ca and Na∕Mg, however, limit the use of these ratios as a proxy and/or make a high number of samples necessary. Additionally, we explore two models to explain the observed temperature sensitivity of Na∕Ca ratios for an open and semi-enclosed calcifying space based on temperature-sensitive Na- and Ca-pumping enzymes and transport proteins that change the composition of the calcifying fluid and consequently the skeletal Na∕Ca ratio.
The Arctic Svalbard Archipelago hosts the world’s northernmost cold-water ‘carbonate factories’ thriving here despite of presumably unfavourable environmental conditions and extreme seasonality. Two contrasting sites of intense biogenic carbonate production, the rhodolith beds in Mosselbukta in the north of the archipelago and the barnacle-mollusc dominated carbonate sediments accumulating in the strong hydrodynamic regime of the Bjørnøy-Banken south of Spitsbergen, were the targets of the RV Maria S. Merian cruise 55 in June 2016. By integrating data from physical oceanography, marine biology, and marine geology, the present contribution characterises the environmental setting and biosedimentary dynamics of these two polar carbonate factories. Repetitive CTD profiling in concert with autonomous temperature/salinity loggers on a long-term settlement platform identified spatiotemporal patterns in the involved Atlantic and Polar water masses, whereas short-term deployments of a lander revealed fluctuations of environmental variables in the rhodolith beds in Mosselbukta and at same depth (46 m) at Bjørnøy-Banken. At both sites, dissolved inorganic nutrients in the water column were found depleted (except for elevated ammonium concentrations) and show an overall increase in concentration and N:P ratios toward deeper waters. This indicates that a recycling system was fuelling primary production after the phytoplankton spring bloom at the time of sampling in June 2016. Accordingly, oxygen levels were found elevated and carbon dioxide concentrations (pCO2) markedly reduced, on average only half the expected equilibrium values. Backed up by seawater stable carbon and oxygen isotope signatures, this is interpreted as an effect of limited air-sea gas exchange during seasonal ice cover in combination with a boost in community photosynthesis during the spring phytoplankton bloom. The observed trends are enhanced by the onset of rhodophyte photosynthesis in the rhodolith beds during the polar day upon retreat of sea-ice. Potential adverse effects of ocean acidification on the local calcifier community are thus predicted to be seasonally buffered by the marked drop in pCO2 during the phase of sea-ice cover and spring phyto-plankton bloom, but this effect will diminish should the seasonal sea-ice formation continue to decline. Among the 25 macrobenthos taxa identified from images captured by the lander’s camera system, all but three species were calcifiers contributing to the carbonate production. Biodiversity was found to be much higher in Mosselbukta (21 taxa) compared to Bjørnøy-Banken (8 taxa), which is considered as a result of enhanced habitat diversity provided in the rhodolith beds by the bioengineering crustose alga Lithothamnion glaciale. Filter-feeding activity of selected key species did reveal group-specific but no common activity patterns. Biotic disturbance of the filtering activity was common, in contrast to abiotic factors, with hermit crabs representing the primary trigger. Motion tracking of rhodoliths revealed a high frequency of dislocation, triggered not by abiotic factors but by the activity of benthic invertebrates, in particular echinoids ploughing below or moving over the rhodoliths. The echinoid Strongylocentrotus sp. is the most abundant component of the associated fauna, thereby considerably contributing both to carbonate production and to grazing bioerosion. Together, these results portray a high degree of seasonal as well as short-term dynamics in environmental conditions that despite many similarities support distinctly different communities and biodiversity patterns in the calcifying macrobenthos at the two studied polar carbonate factories.
The exchange of heat, momentum, and mass in the atmosphere over mountainous terrain is controlled by synoptic-scale dynamics, thermally driven mesoscale circulations, and turbulence. This article reviews the key challenges relevant to the understanding of exchange processes in the mountain boundary layer and outlines possible research priorities for the future. The review describes the limitations of the experimental study of turbulent exchange over complex terrain, the impact of slope and valley breezes on the structure of the convective boundary layer, and the role of intermittent mixing and wave–turbulence interaction in the stable boundary layer. The interplay between exchange processes at different spatial scales is discussed in depth, emphasizing the role of elevated and ground-based stable layers in controlling multi-scale interactions in the atmosphere over and near mountains. Implications of the current understanding of exchange processes over mountains towards the improvement of numerical weather prediction and climate models are discussed, considering in particular the representation of surface boundary conditions, the parameterization of sub-grid-scale exchange, and the development of stochastic perturbation schemes.
Throughout mankind’s history, the need to secure and protect the home settlement was an essential one. This holds especially true for the city of Ainos (modern Enez) in Turkish Thrace. Due to its continuous settlement history since the 7th/6th century BC, several different types of city walls were built—sometimes even on top of each other—several of which have been preserved over time. To decipher the construction style, the course and the age of a buried city wall segment in the southern part of the former city, a geoscientific multi-proxy approach including magnetic gradiometer (MG) and electrical resistivity tomography (ERT) measurements in combination with granulometrical, sedimentological and microfaunistical investigations on sediment cores was applied. We were able to (1) present reasonable arguments for its Hellenistic age; (2) reveal the course of this wall segment and extrapolate it further north into a less studied area; and (3) demonstrate that in this near-coastal area, the former swampy terrain had been consolidated for constructing the wall. Our multi-proxy approach serves as a valuable example for investigating buried structures in archaeological contexts, avoiding a less-economical, time-consuming, or even forbidden excavation.
Atmospheric aerosols and their effect on clouds are thought to be important for anthropogenic radiative forcing of the climate, yet remain poorly understood1. Globally, around half of cloud condensation nuclei originate from nucleation of atmospheric vapours2. It is thought that sulfuric acid is essential to initiate most particle formation in the atmosphere3,4, and that ions have a relatively minor role5. Some laboratory studies, however, have reported organic particle formation without the intentional addition of sulfuric acid, although contamination could not be excluded6,7. Here we present evidence for the formation of aerosol particles from highly oxidized biogenic vapours in the absence of sulfuric acid in a large chamber under atmospheric conditions. The highly oxygenated molecules (HOMs) are produced by ozonolysis of α-pinene. We find that ions from Galactic cosmic rays increase the nucleation rate by one to two orders of magnitude compared with neutral nucleation. Our experimental findings are supported by quantum chemical calculations of the cluster binding energies of representative HOMs. Ion-induced nucleation of pure organic particles constitutes a potentially widespread source of aerosol particles in terrestrial environments with low sulfuric acid pollution.
We present a study characterizing aerosol particles resulting from a skyscraper blasting. High mass concentrations with a maximum of 844.9 μg m-3 were present for a short time period of approximately 15 minutes. They result in a day mean of 32.6 μg m-3 compared to a 27.6 μg m-3 background not exceeding the 50 μg m-3 EU maximum permissive value. The increase in particle number concentration was less pronounced with a maximum concentration of 6.9 ⋅ 104 cm-3 compared to the local background value of 1.8 ⋅ 104 cm-3. The size-resolved number concentration shows a single mode of ultrafine particles at approximately 93 nm. The spatial distribution of deposited dust was investigated with Bergerhoff glass collection vessels, showing a decrease with distance. In the deposited dust samples the concentrations of twelve metals was determined, non of them exceeded the regional background concentrations significantly. The chemical composition of individual particles emitted by the demolition was studied by Scanning Electron Microscopy. They were mainly concrete and steel particles, with 60% calcium carbonates, 19% calcium sulfates, 19% silicates and 2% steel. In energy-dispersive X-Ray Spectroscopy, no fibers like asbestos were observed. Using a broad spectrum of instruments and methods, we obtain comprehensive characterization of the particles emitted by the demolition.
A twentieth century-long coupled atmosphere-ocean regional climate simulation with COSMO-CLM (Consortium for Small-Scale Modeling, Climate Limited-area Model) and NEMO (Nucleus for European Modelling of the Ocean) is studied here to evaluate the added value of coupled marginal seas over continental regions. The interactive coupling of the marginal seas, namely the Mediterranean, the North and the Baltic Seas, to the atmosphere in the European region gives a comprehensive modelling system. It is expected to be able to describe the climatological features of this geographically complex area even more precisely than an atmosphere-only climate model. The investigated variables are precipitation and 2 m temperature. Sensitivity studies are used to assess the impact of SST (sea surface temperature) changes over land areas. The different SST values affect the continental precipitation more than the 2 m temperature. The simulated variables are compared to the CRU (Climatic Research Unit) observational data, and also to the HOAPS/GPCC (Hamburg Ocean Atmosphere Parameters and Fluxes from Satellite Data, Global Precipitation Climatology Centre) data. In the coupled simulation, added skill is found primarily during winter over the eastern part of Europe. Our analysis shows that, over this region, the coupled system is dryer than the uncoupled system, both in terms of precipitation and soil moisture, which means a decrease in the bias of the system. Thus, the coupling improves the simulation of precipitation over the eastern part of Europe, due to cooler SST values and in consequence, drier soil.
Historically, the expansion of soy plantations has been a major driver of land-use/cover change (LUCC) in Brazil. While a series of recent public actions and supply-chain commitments reportedly curbed the replacement of forests by soy, the expansion of the agricultural commodity still poses a considerable threat to the Amazonian and Cerrado biomes. Identification of areas under high risk of soy expansion is thus paramount to assist conservation efforts in the region. We mapped the areas suitable for undergoing transition to soy plantations in the Legal Amazon with a machine-learning approach adopted from the ecological modeling literature. Simulated soy expansion for the year 2014 exhibited favorable validation scores compared to other LUCC models. We then used our model to simulate how potential future infrastructure improvements would affect the 2014 probabilities of soy occurrence in the region. In addition to the 2.3 Mha of planted soy in the Legal Amazon in 2014, our model identified another 14.7 Mha with high probability of soy conversion in the region given the infrastructure conditions at that time. Out of those, pastures and forests represented 9.8 and 0.4 Mha, respectively. Under the new infrastructure scenarios simulated, the Legal Amazonian area under high risk of soy conversion increased by up to 2.1 Mha (14.6%). These changes led to up to 11.4 and 51.4% increases in the high-risk of conversion areas of pastures and forests, respectively. If conversion occurs in the identified high-risk areas, at least 4.8 Pg of CO2 could be released into the atmosphere, a value that represents 10 times the total CO2 emissions of Brazil in 2014. Our results highlight the importance of targeting conservation policies and enforcement actions, including the Soy Moratorium, to mitigate future forest cover loss associated with infrastructure improvements in the region.
We evaluate the near-surface representation of thermally driven winds in the Swiss Alps in a numerical weather prediction model at km-scale resolution. In addition, the influence of grid resolution (2.2 km and 1.1 km), topography filtering, and land surface datasets on the accuracy of the simulated valley winds is investigated. The simulations are evaluated against a comprehensive set of surface observations for an 18-day fair-weather summer period in July 2006. The episode is characterized by strong diurnal wind systems and the formation of shallow convection over the mountains, which transitions to precipitating convection in some areas. The near-surface winds (10 m above ground level) follow a typical diurnal pattern with strong daytime up-valley flow and weaker nighttime down-valley flow. At a 2.2 km resolution the valley winds are poorly simulated for most stations, while at a 1.1 km resolution the diurnal cycle of the valley winds is well represented in most large (e.g., Rhein valley at Chur and Rhone valley at Visp) and medium-sized valleys (e.g., Linth valley at Glarus). In the smaller valleys (e.g., Maggia valley at Cevio), the amplitude of the valley wind is still significantly underestimated, even at a 1.1 km resolution. Detailed sensitivity experiments show that the use of high-resolution land surface datasets, for both the soil characteristics as well as for the land cover, and reduced filtering of the topography are essential to achieve good performance at a 1.1 km resolution
The synthesis of polynitrogen compounds is of fundamental importance due to their potential as environmentally-friendly high energy density materials. Attesting to the intrinsic difficulties related to their formation, only three polynitrogen ions, bulk stabilized as salts, are known. Here, magnesium and molecular nitrogen are compressed to about 50 GPa and laser-heated, producing two chemically simple salts of polynitrogen anions, MgN4 and Mg2N4. Single-crystal X-ray diffraction reveals infinite anionic polythiazyl-like 1D N-N chains in the crystal structure of MgN4 and cis-tetranitrogen N44− units in the two isosymmetric polymorphs of Mg2N4. The cis-tetranitrogen units are found to be recoverable at atmospheric pressure. Our results respond to the quest for polynitrogen entities stable at ambient conditions, reveal the potential of employing high pressures in their synthesis and enrich the nitrogen chemistry through the discovery of other nitrogen species, which provides further possibilities to design improved polynitrogen arrangements.
Abiotic formation of n-alkane hydrocarbons has been postulated to occur within Earth's crust. Apparent evidence was primarily based on uncommon carbon and hydrogen isotope distribution patterns that set methane and its higher chain homologues apart from biotic isotopic compositions associated with microbial production and closed system thermal degradation of organic matter. Here, we present the first global investigation of the carbon and hydrogen isotopic compositions of n-alkanes in volcanic-hydrothermal fluids hosted by basaltic, andesitic, trachytic and rhyolitic rocks. We show that the bulk isotopic compositions of these gases follow trends that are characteristic of high temperature, open system degradation of organic matter. In sediment-free systems, organic matter is supplied by surface waters (seawater, meteoric water) circulating through the reservoir rocks. Our data set strongly implies that thermal degradation of organic matter is able to satisfy isotopic criteria previously classified as being indicative of abiogenesis. Further considering the ubiquitous presence of surface waters in Earth’s crust, abiotic hydrocarbon occurrences might have been significantly overestimated.
To quantify water flows between groundwater (GW) and surface water (SW) as well as the impact of Abstract. To quantify water flows between groundwater (GW) and surface water (SW) as well as the impact of capillary rise on evapotranspiration by global hydrological models (GHMs), it is necessary to replace the bucket-like linear GW reservoir model typical for hydrological models with a fully integrated gradient-based GW flow model. Linear reservoir models can only simulate GW discharge to SW bodies, provide no information on the location of the GW table and assume that there is no GW flow among grid cells. A gradient-based GW model simulates not only GW storage but also hydraulic head, which together with information on SW table elevation enables the quantification of water flows from GW to SW and vice versa. In addition, hydraulic heads are the basis for calculating lateral GW flow among grid cells and capillary rise.
G³M is a new global gradient-based GW model with a spatial resolution of 5' that will replace the current linear GW reservoir in the 0.5° WaterGAP Global Hydrology Model (WGHM). The newly developed model framework enables inmemory coupling to WGHM while keeping overall runtime relatively low, allowing sensitivity analyses and data assimilation. This paper presents the G³M concept and specific model design decisions together with results under steady-state naturalized conditions, i.e. neglecting GW abstractions. Cell-specific conductances of river beds, which govern GW-SW interaction, were determined based on the 30'' steady-state water table computed by Fan et al. (2013). Together with an appropriate choice for the effective elevation of the SW table within each grid cell, this enables a reasonable simulation of drainage from GW to SW such that, in contrast to the GW model of de Graaf et al. (2015, 2017), no additional drainage based on externally provided values for GW storage above the floodplain is required in G³M. Comparison of simulated hydraulic heads to observations around the world shows better agreement than de Graaf et al. (2015). In addition, G³M output is compared to the output of two established macro-scale models for the Central Valley, California, and the continental United States, respectively. As expected, depth to GW table is highest in mountainous and lowest in flat regions. A first analysis of losing and gaining rivers and lakes/wetlands indicates that GW discharge to rivers is by far the dominant flow, draining diffuse GW recharge, such that lateral flows only become a large fraction of total diffuse and focused recharge in case of losing rivers and some areas with very low GW recharge. G³M does not represent losing rivers in some dry regions. This study presents the first steps towards replacing the linear GW reservoir model in a GHM while improving on recent efforts, demonstrating the feasibility of the approach and the robustness of the newly developed framework.
In global hydrological models, groundwater (GW) is typically represented by a bucket-like linear groundwater reservoir. Reservoir models, however, (1) can only simulate GW discharge to surface water (SW) bodies but not recharge from SW to GW, (2) provide no information on the location of the GW table, and (3) assume that there is no GW flow among grid cells. This may lead, for example, to an underestimation of groundwater resources in semiarid areas where GW is often replenished by SW or to an underestimation of evapotranspiration where the GW table is close to the land surface. To overcome these limitations, it is necessary to replace the reservoir model in global hydrological models with a hydraulic head gradient-based GW flow model.
We present G3M, a new global gradient-based GW model with a spatial resolution of 5′ (arcminutes), which is to be integrated into the 0.5∘ WaterGAP Global Hydrology Model (WGHM). The newly developed model framework enables in-memory coupling to WGHM while keeping overall runtime relatively low, which allows sensitivity analyses, calibration, and data assimilation. This paper presents the G3M concept and model design decisions that are specific to the large grid size required for a global-scale model. Model results under steady-state naturalized conditions, i.e., neglecting GW abstractions, are shown. Simulated hydraulic heads show better agreement to observations around the world compared to the model output of de Graaf et al. (2015). Locations of simulated SW recharge to GW are found, as is expected, in dry and mountainous regions but areal extent of SW recharge may be underestimated. Globally, GW discharge to rivers is by far the dominant flow component such that lateral GW flows only become a large fraction of total diffuse and focused recharge in the case of losing rivers, some mountainous areas, and some areas with very low GW recharge. A strong sensitivity of simulated hydraulic heads to the spatial resolution of the model and the related choice of the water table elevation of surface water bodies was found. We suggest to investigate how global-scale groundwater modeling at 5′ spatial resolution can benefit from more highly resolved land surface elevation data.
Latitudinal and bathymetrical species richness patterns in the NW Pacific and adjacent Arctic Ocean
(2019)
Global scale analyses have recently revealed that the latitudinal gradient in marine species richness is bimodal, peaking at low-mid latitudes but with a dip at the equator; and that marine species richness decreases with depth in many taxa. However, these overall and independently studied patterns may conceal regional differences that help support or qualify the causes in these gradients. Here, we analysed both latitudinal and depth gradients of species richness in the NW Pacific and its adjacent Arctic Ocean. We analysed 324,916 distribution records of 17,414 species from 0 to 10,900 m depth, latitude 0 to 90°N, and longitude 100 to 180°N. Species richness per c. 50 000 km2 hexagonal cells was calculated as alpha (local average), gamma (regional total) and ES50 (estimated species for 50 records) per latitudinal band and depth interval. We found that average ES50 and gamma species richness decreased per 5° latitudinal bands and 100 m depth intervals. However, average ES50 per hexagon showed that the highest species richness peaked around depth 2,000 m where the highest total number of species recorded. Most (83%) species occurred in shallow depths (0 to 500 m). The area around Bohol Island in the Philippines had the highest alpha species richness (more than 8,000 species per 50,000 km2). Both alpha and gamma diversity trends increased from the equator to latitude 10°N, then further decreased, but reached another peak at higher latitudes. The latitudes 60–70°N had the lowest gamma and alpha diversity where there is almost no ocean in our study area. Model selection on Generalized Additive Models (GAMs) showed that the combined effects of all environmental predictors produced the best model driving species richness in both shallow and deep sea. The results thus support recent hypotheses that biodiversity, while highest in the tropics and coastal depths, is decreasing at the equator and decreases with depth below ~2000 m. While we do find the declines of species richness with latitude and depth that reflect temperature gradients, local scale richness proved poorly correlated with many environmental variables. This demonstrates that while regional scale patterns in species richness may be related to temperature, that local scale richness depends on a greater variety of variables.