Geowissenschaften
Refine
Year of publication
Document Type
- Article (593)
- Doctoral Thesis (139)
- Book (25)
- Contribution to a Periodical (21)
- Conference Proceeding (19)
- Working Paper (19)
- Part of Periodical (10)
- Diploma Thesis (8)
- diplomthesis (7)
- Preprint (6)
Is part of the Bibliography
- no (872)
Keywords
- climate change (13)
- Climate change (5)
- Geochemistry (5)
- Klimaänderung (5)
- Atmospheric chemistry (4)
- Boden (4)
- Deutschland (4)
- Klima (4)
- Modellierung (4)
- Stratosphäre (4)
Institute
- Geowissenschaften (872)
- Senckenbergische Naturforschende Gesellschaft (59)
- Präsidium (49)
- Biodiversität und Klima Forschungszentrum (BiK-F) (45)
- Geographie (26)
- Biowissenschaften (15)
- Medizin (11)
- Physik (7)
- Institut für Ökologie, Evolution und Diversität (6)
- Biochemie und Chemie (5)
Der Paläontologe und das Meer : Eberhard Gischler erforscht das Klimagedächtnis von Korallenriffen
(2018)
Constraining the architecture of complex 3D volcanic plumbing systems within active rifts, and their impact on rift processes, is critical for examining the interplay between faulting, magmatism and magmatic fluids in developing rift segments. The Natron basin of the East African Rift System provides an ideal location to study these processes, owing to its recent magmatic-tectonic activity and ongoing active carbonatite volcanism at Oldoinyo Lengai. Here, we report seismicity and fault plane solutions from a 10-month temporary seismic network spanning Oldoinyo Lengai, Naibor Soito volcanic field and Gelai volcano. We locate 6827 earthquakes with ML -0.85 to 3.6, which are related to previous and ongoing magmatic and volcanic activity in the region, as well as regional tectonic extension. We observe seismicity down to ~17 km depth north and south of Oldoinyo Lengai and shallow seismicity (3 - 10 km) beneath Gelai, including two swarms. The deepest seismicity (~down to 20 km) occurs above a previously imaged magma body below Naibor Soito. These seismicity patterns reveal a detailed image of a complex volcanic plumbing system, supporting potential lateral and vertical connections between shallow- and deep-seated magmas, where fluid and melt transport to the surface is facilitated by intrusion of dikes and sills. Focal mechanisms vary spatially. T-axis trends reveal dominantly WNW-ESE extension near Gelai, while strike-slip mechanisms and a radial trend in P-axes are observed in the vicinity of Oldoinyo Lengai. These data support local variations in the state of stress, resulting from a combination of volcanic edifice loading and magma-driven stress changes imposed on a regional extensional stress field. Our results indicate that the southern Natron basin is a segmented rift system, in which fluids preferentially percolate vertically and laterally in a region where strain transfers from a border fault to a developing magmatic rift segment.
The complete elastic stiffness tensor of thiourea has been determined from thermal diffuse scattering (TDS) using high-energy photons (100 keV). Comparison with earlier data confirms a very good agreement of the tensor coefficients. In contrast with established methods to obtain elastic stiffness coefficients (e.g. Brillouin spectroscopy, inelastic X-ray or neutron scattering, ultrasound spectroscopy), their determination from TDS is faster, does not require large samples or intricate sample preparation, and is applicable to opaque crystals. Using high-energy photons extends the applicability of the TDS-based approach to organic compounds which would suffer from radiation damage at lower photon energies.
Strain localization in the lithosphere and the formation, evolution, and maintenance of resulting plate boundaries play a crucial role in plate tectonics and thermo‐chemical mantle convection. Previously activated lithospheric deformation zones often appear to maintain a “memory” of weakening, leading to tectonic inheritance within plate reorganizations including the Wilson cycle. Different mechanisms have been proposed to explain such strain localization, but it remains unclear which operates on what spatio‐temporal scales, and how to best incorporate them in large‐scale mantle convection models. Here, we analyze two candidates, (1), grain‐size sensitive rheology and, (2), damage‐style parameterizations of yield, stress which are sometimes used to approximate the former. Grain‐size reduction due to dynamic recrystallization can drive localization in the ductile domain, and grain growth provides a time‐dependent rheological hardening component potentially enabling the preservation of rheological heterogeneities. We compare the dynamic weakening and hardening effects as well as the timescales of strength evolution for a composite rheology including grain‐size dynamics with a pseudo‐plastic rheology including damage‐ (or “strain”‐) dependent weakening. We explore the implications of different proposed grain‐size evolution laws, and test to which extent strain‐dependent rheologies can mimic the weakening and hardening effects of the more complex micro‐physical behavior. Such an analysis helps to better understand the parallels and differences between various strain‐localization modeling approaches used in different tectonics and geodynamics communities. More importantly, our results contribute to efforts to identify the key ingredients of strain‐localization and damage hysteresis within plate tectonics and how to represent those in planetary‐scale modeling.
The Ceboruco is a 2280 m high stratovolcano located in Nayarit State, Mexico. Despite its last eruption which occurred in 1870, it is the most active volcano in the area, showing volcanicearthquake activity together with ongoing vapor emissions. The magnetotelluric survey was carried out in November 2016. It was part of a geothermal project (CeMIEGeo-P24) and focused on the determination of the electrical conductivity distribution in the subsurface of the volcano.
The Magnetotelluric Apparent Resistivity Tensor, as introduced by Brown (2016), can be decomposed into an amplitude and a phase tensor. The fundamental physics behind those new tensors were presented in Hering et al. (2019), using canonical models in 1-D (isotropic and anisotropic) and 2-D resistivity environments. Here, the tensors are introduced for a high-quality data set, where their interpretational benefits become very obvious. Additionally, results from an isotropic 3-D inversion are presented and compared to an alternative 3-D anisotropic forward model.
During the first two days of August 2016 a seismic crisis occurred on Brava, Cape Verde, which – according to observations based on a local seismic network – was characterized by more than thousand volcano–seismic signals. Brava is considered an active volcanic island, although it has not experienced any historic eruptions. Seismicity significantly exceeded the usual level during the crisis. We report on results based on data from a temporary seismic–array deployment on the neighbouring island of Fogo at a distance of about 35 km. The array was in operation from October 2015 to December 2016 and recorded a total of 1343 earthquakes, 355 thereof were localized. On 1 and 2 August we observed 54 earthquakes, 25 of which could be located beneath Brava. We further evaluate the observations with regards to possible precursors to the crisis and its continuation. Our analysis shows a migration of seismicity around Brava, but no distinct precursory pattern. However, the observations suggest that similar earthquake swarms commonly occur close to Brava. The results further confirm the advantages of seismic arrays as tools for the remote monitoring of regions with limited station coverage or access.
During the first two days of August 2016 a seismic crisis occurred on Brava, Cabo Verde, which – according to observations based on a local seismic network – was characterized by more than a thousand volcano-seismic signals. Brava is considered an active volcanic island, although it has not experienced any historic eruptions. Seismicity significantly exceeded the usual level during the crisis. We report on results based on data from a temporary seismic-array deployment on the neighbouring island of Fogo at a distance of about 35 km. The array was in operation from October 2015 to December 2016 and recorded a total of 1343 earthquakes in the region of Fogo and Brava; 355 thereof were localized. On 1 and 2 August we observed 54 earthquakes, 25 of which could be located beneath Brava. We further evaluate the observations with regards to possible precursors to the crisis and its continuation. Our analysis shows a significant variation in seismicity around Brava, but no distinct precursory pattern. However, the observations suggest that similar earthquake swarms commonly occur close to Brava. The results further confirm the advantages of seismic arrays as tools for the remote monitoring of regions with limited station coverage or access.
Rodrigues Ridge connects the Réunion hotspot track with the Central Indian Ridge (CIR) and has been suggested to represent the surface expression of a sub-lithospheric flow channel. From global earthquake catalogues, the seismicity in the region has been associated mainly with events related to the fracture zones at the CIR. However, some segments of the CIR appear void of seismic events. Here, we report on the seismicity recorded at a temporary array of 10 seismic stations operating on Rodrigues Island from September 2014 to June 2016. The array analysis was performed in the time domain by time shifting and stacking the complete waveforms. Event distances were estimated based on a 1-D velocity model and the travel time differences between S and P wave arrivals. We detected and located 63 new events that were not reported by the global networks. Most of the events (51) are located off the CIR and can be classified as intraplate earthquakes. Local magnitudes varied between 1.6 and 3.7. Four seismic clusters were observed that occurred to the west of the spreading segment of the CIR. The Rodrigues Ridge appeared to be aseismic during the period of operation. The lack of seismic activity along both Rodrigues Ridge and the sections of the CIR to the east of Rodrigues may be explained by partially molten upper-mantle material, possibly in relation to the proposed material flow between the Réunion plume and the CIR.
Tetra-auricupride, ideally AuCu, represents the only species showing the coexistence of Au with an elevated level of Pt, as in the case of a detrital grain studied structurally for the first time, from an ophiolite-associated placer at Bolshoy Khailyk, western Sayans, Russia. We infer that tetra-auricupride can incorporate as much as ~30 mol. % of a “PtCu” component, apparently without significant modification of the unit cell. The unit-cell parameters of platiniferous tetra-auricupride are: a 2.790(1) Å, c 3.641(4) Å, with c/a = 1.305, which are close to those reported for ordered AuCu(I) in the system Au–Cu, and close also to the cell parameters of tetraferroplatinum (PtFe), which both appear to crystallize in the same space group, P4/mmm. These intermetallic compounds and natural alloys are thus isostructural. The closeness of their structures presumably allows Pt to replace Au atoms so readily. The high extent of Cu + Au enrichment is considered to be a reflection of geochemical evolution and buildup in levels of the incompatible Cu and Au with subordinate Pt in a remaining volume of melt at low levels of fO2 and fS2 in the system.
In the nineteenth century, two Neolithic axe-heads were reported from the Michelsberg enclosure system at Kapellenberg. The recent identification of an unusually large tumulus, from which the axe-heads were almost certainly once recovered, reveals that socio-political hierarchisation, linked to the emergence of high-ranking elites in Brittany and the Paris Basin during the fifth millennium cal BC, may have extended into Central Europe.
Wetlands such as bogs, swamps, or freshwater marshes are hotspots of biodiversity. For 5.1 million km2 of inland wetlands, the dynamics of area and water storage, which strongly impact biodiversity and ecosystem services, were simulated using the global hydrological model WaterGAP. For the first time, the impacts of both human water use and man‐made reservoirs (WUR) and future climate change (CC) on wetlands around the globe were quantified. WUR impacts are concentrated in arid/semiarid regions, where WUR decreased mean wetland water storage by more than 5% on 8.2% of the mean wetland area during 1986–2005 (Am), with highest decreases in groundwater depletion area. Using output of three climate models, CC impacts on wetlands were quantified, distinguishing unavoidable impacts [i.e., at 2 °C global warming (GW)] from avoidable impacts (difference between 3 °C and 2 °C impacts). Even unavoidable CC impacts are projected to be much larger than WUR impacts, also in arid/semiarid regions. On most wetland area with reliable estimates, avoidable CC impacts are more than twice as large as unavoidable impacts. In case of 2 °C GW, half of Am is estimated to be unaffected by mean storage changes of more than 5%, but only one third in case of 3 °C GW. Temporal variability of water storage will increase for most wetlands. Wetlands in dry regions will be affected the most, particularly by water storage decreases in the dry season. Different from wealthier countries, low‐income countries will dominantly suffer from a decrease in wetland water storage due to CC.
In November 2016, magnetotelluric (MT) data were collected at the Ceboruco Volcano in cooperation with the Centro de Sismología y Volcanología de Occidente (SisVoc, Universidad de Guadalajara, Mexico). The Ceboruco is a 2280 m high stratovolcano, located in Nayarit State, Mexico. It is placed in the central part of the Tepic-Zacoalco Rift (TZR), which constitutes the north-western end of the Trans-Mexican Volcanic Belt. Together with Chapala and Colima (in the Jalisco Block), they form the triple rift system developed as a consequence of the ongoing subduction of the Rivera and Cocos oceanic plates beneath the North American continental crust. Although its last eruption occurred in 1870, it is the most active volcano in the area, showing volcanic-earthquake activity together with ongoing vapor emissions. The survey was part of a geothermal project (CeMIEGeo-P24) and focused on the determination of electrical conductivity properties to characterize the deep structure and the geothermal potential of the Volcano. Frequency dependent magnetotelluric response functions were calculated from 25 broadband MT stations, which covered an area of 10 x 10 km2 including its crater, calderas and foreland. The results were interpreted using anisotropic 3-D forward modelling and isotropic 3-D inversion approaches, considering strong topographical effects. The final resistivity model implies a highly conductive layer, reaching from near-surface to approximately 2 km depth, which might be related to a hydrothermal system. Here, mineralized fluids and clay minerals can cause high conductivities around 1 S/m. For longer periods, the principal axes of the MT response tensors (phase tensor, apparent resistivity tensor) are in good agreement with the strike direction of the underlying rift system. However, they are not rendered by the isotropic inversion. Thus the data suggest an anisotropic electrical conductivity at greater depth with its principal axis determined by the response tensors.
Anthropogenic climate change is expected to impact ecosystem structure, biodiversity and ecosystem services in Africa profoundly. We used the adaptive Dynamic Global Vegetation Model (aDGVM), which was originally developed and tested for Africa, to quantify sources of uncertainties in simulated African potential natural vegetation towards the end of the 21st century. We forced the aDGVM with regionally downscaled high‐resolution climate scenarios based on an ensemble of six general circulation models (GCMs) under two representative concentration pathways (RCPs 4.5 and 8.5). Our study assessed the direct effects of climate change and elevated CO2 on vegetation change and its plant‐physiological drivers. Total increase in carbon in aboveground biomass in Africa until the end of the century was between 18% to 43% (RCP4.5) and 37% to 61% (RCP8.5) and was associated with woody encroachment into grasslands and increased woody cover in savannas. When direct effects of CO2 on plants were omitted, woody encroachment was muted and carbon in aboveground vegetation changed between –8 to 11% (RCP 4.5) and –22 to –6% (RCP8.5). Simulated biome changes lacked consistent large‐scale geographical patterns of change across scenarios. In Ethiopia and the Sahara/Sahel transition zone, the biome changes forecast by the aDGVM were consistent across GCMs and RCPs. Direct effects from elevated CO2 were associated with substantial increases in water use efficiency, primarily driven by photosynthesis enhancement, which may relieve soil moisture limitations to plant productivity. At the ecosystem level, interactions between fire and woody plant demography further promoted woody encroachment. We conclude that substantial future biome changes due to climate and CO2 changes are likely across Africa. Because of the large uncertainties in future projections, adaptation strategies must be highly flexible. Focused research on CO2 effects, and improved model representations of these effects will be necessary to reduce these uncertainties.
Precipitation extremes with devastating socioeconomic consequences within the South American Monsoon System (SAMS) are expected to become more frequent in the near future. The complexity in SAMS behavior, however, poses severe challenges for reliable future projections. Thus, robust paleomonsoon records are needed to constrain the high spatiotemporal variability in the response of SAMS rainfall to different climatic drivers. This study uses Ti/Ca ratios from X‐ray fluorescence scanning of a sediment core retrieved off eastern Brazilian to trace precipitation changes over the past 322 Kyr. The results indicate that despite the spatiotemporal complexity of the SAMS, insolation forcing is the primary pacemaker of variations in the monsoonal system. Additional modulation by atmospheric pCO2 suggests that SAMS intensity over eastern Brazil will be suppressed by rising CO2 emissions in the future. Lastly, our record reveals an unprecedented strong and persistent wet period during Marine Isotope Stage 6 driven by anomalously strong trade winds.
The complex magnetotelluric (MT) apparent resistivity tensor can be decomposed into two real tensors, the apparent resistivity and the resistivity phase tensors, which represent relationships between the observed electric field at a point on the Earth's surface and an associated apparent current density. We explain the differences between these tensors and conventional estimates of apparent resistivity and phase for simple resistivity environments and demonstrate, using canonical models in 1‐D and 2‐D environments, that both tensors are more sensitive to vertical and horizontal resistivity gradients than their conventional counterparts. The properties of the new tensors are explained using electromagnetic induction theory and the effects of associated charges at resistivity boundaries. We introduce a new way to plot tensor ellipses, which brings significant improvements to the interpretation of MT data, using appropriate visualization software. The apparent resistivity tensor gives information about the magnitude and direction of apparent resistivity subsurface structures and has a strong response to vertical resistivity contrasts. The resistivity phase tensor is highly sensitive to vertical boundaries and the associated fields in the TM mode. It is also free from static distortions under the same conditions implied for the conventional phase tensor. These findings have prompted a study in the potential of the new tensors for 3‐D inversions. The results from a 3‐D inversion of a canonical oblique conductor straddling two quarter spaces show distinct improvements in resolving the boundaries of the conductor and open a promising field for future studies.
The endemic argan woodlands cover large parts of South Morocco and create a characteristic landscape with areas of sparsely vegetated and bare soil surfaces between single trees. This unique ecosystem has been under extensive agrosilvopastoral management for centuries and is now at risk of degradation caused by overgrazing and increasing scarcity and variability of rainfall.
To investigate susceptibility to wind erosion, we conducted an experimental–empirical study including wind tunnel tests and a drone‐generated digital elevation model and quantified wind‐erodible material on five different associated surface types by means of sediment catchers. The highest emission flux was measured on freshly ploughed surfaces (1875 g m–2 h–1), while older ploughed areas with a re‐established crust produced a much lower emission flux (795 g m–2 h–1). Extensive tillage may have been a sustainable practice for generations, but increasing drought and uncertainty of rainfall now lead to an acute risk of severe soil erosion and dust production. The typical crusted surfaces characterized by residual rock fragment accumulation and wash processes produced the second highest emission flux (1,354 g m–2 h–1). Material collected from tree‐shaded areas (933 g m–2 h–1) was revealed to be a considerable source of organic material, possibly affecting substrate conditions positively on a larger regional scale. The lowest flux was measured on rock fragment‐covered surfaces (301 g m–2 h–1).
The data show that open argan woodland may be a considerable source for wind erosion and dust production, depending on surface characteristics strongly related to management. An adapted management must include the conservation of argan trees to offer a promising approach to prevent severe wind erosion and dust production and mitigate possible impacts of land‐use change and climate change related shifts in wind and rainfall patterns.
A realistic simulation of the atmospheric boundary layer (ABL) depends on an accurate representation of the land–atmosphere coupling. Land surface temperature (LST) plays an important role in this context and the assimilation of LST can lead to improved estimates of the boundary layer and its processes. We assimilated synthetic satellite LST retrievals derived from a nature run as truth into a fully coupled, state‐of‐the‐art land–atmosphere numeric weather prediction model. As assimilation system a local ensemble transform Kalman filter was used and the control vector was augmented by the soil temperature and humidity. To evaluate the concept of the augmented control vector, two‐day case‐studies with different control vector settings were conducted for clear‐sky periods in March and August 2017. These experiments with hourly LST assimilation were validated against the nature run and overall, the RMSE of atmospheric and soil temperature of the first‐guess (and analysis) were reduced. The temperature estimate of the ABL was particularly improved during daytime as was the estimate of the soil temperature during the whole diurnal cycle. The best impact of LST assimilation on the soil and the ABL was achieved with the augmented control vector. Through the coupling between the soil and the atmosphere, the assimilation of LST can have a positive impact on the temperature forecast of the ABL even after 15 hr because of the memory of the soil. These encouraging results motivate further work towards the assimilation of real satellite LST retrievals.
The continental expression of global cooling during the Miocene Climate Transition in Central Asia is poorly documented, as the tectonically active setting complicates the correlation of Neogene regional and global climatic developments. This study presents new geochemical data (CaSO4 content, carbonate δ13C and δ18O) from the endorheic alluvial‐lacustrine Aktau succession (Ili Basin, south‐east Kazakhstan) combined with findings from the previously published facies evolution. Time series analysis revealed long‐eccentricity forcing of the paleohydrology throughout the entire succession, split into several facies‐dependent segments. Orbital tuning, constrained by new laser ablation U‐Pb dates and a preexisting magnetostratigraphy, places the succession in a 5.0 Ma long interval in the middle to late Miocene (15.6 to 10.6 Ma). The long‐term water accumulation in the Ili Basin followed the timing of the Miocene Climate Transition, suggesting increased precipitation in the catchment area in response to climate cooling and stronger westerly winds. This was paced by minima of the 2.4 Ma eccentricity cycle, which favored the establishment of a discharge playa (~14.3 Ma) and a perennial lake (12.6 to 11.8 Ma). Furthermore, low obliquity amplitudes (nodes) caused a transient weakening of the westerlies at ~13.7 to 13.5 Ma and at ~12.7 Ma, resulting in negative hydrological budgets and salinization. Flooding of the windward Ili Basin coeval with aridification in the leeward basins suggests that the Tian Shan was a climate boundary already in the middle Miocene. Our results emphasize the impact of climate fluctuations on the westerlies' strength and thus on Central Asian hydrology.
An easy-to-use model to evaluate conductivities at high and middle latitudes in the height range 70–100 km is presented. It is based on electron density profiles obtained with the EISCAT VHF radar during 11 years and on the neutral atmospheric model MSIS95. The model uses solar zenith angle, geomagnetic activity and season as input parameters. It was mainly constructed to study the properties of Schumann resonances that depend on such conductivity profiles.
Remote sensing data are essential for monitoring the Earth’s surface waters, especially since the amount of publicly available in-situ data is declining. Satellite altimetry provides valuable information on the water levels and variations of lakes, reservoirs and rivers. In combination with satellite imagery, the derived time series allow the monitoring of lake storage changes and river discharge. However, satellite altimetry is limited in terms of its spatial resolution due to its measurement geometry, only providing information in the nadir direction beneath the satellite’s orbit. In a case study in the Mississippi River Basin (MRB), this study investigates the potential and limitations of past and current satellite missions for the monitoring of basin-wide storage changes. For that purpose, an automated target detection is developed and the extracted lake surfaces are merged with the satellites’ tracks. This reveals that the current altimeter configuration misses about 80% of all lakes larger than 0.1 km2 in the MRB and 20% of lakes larger than 10 km2, corresponding to 30% and 7% of the total water area, respectively. Past altimetry configurations perform even more poorly. From the larger water bodies represented by a global hydrology model, at least 91% of targets and 98% of storage changes are captured by the current altimeter configuration. This will improve significantly with the launch of the planned Surface Water and Ocean Topography (SWOT) mission.
Increases in water demand often result in unsustainable water use, leaving insufficient amounts of water for the environment. Therefore, water-saving strategies have been introduced to the environmental policy agenda in many (semi)-arid regions. As many such interventions failed to reach their objectives, a comprehensive tool is needed to assess them. We introduced a constructive framework to assess the proposed strategies by estimating five key components of the water balance in an area: (1) Demand; (2) Availability; (3) Withdrawal; (4) Depletion and (5) Outflow. The framework was applied to assess the Urmia Lake Restoration Program (ULRP) which aimed to increase the basin outflow to the lake to reach 3.1 × 109 m3 yr−1. Results suggested that ULRP could help to increase the Outflow by up to 57%. However, successful implementation of the ULRP was foreseen to be impeded because of three main reasons: (i) decreasing return flows; (ii) increased Depletion; (iii) the impact of climate change. Decreasing return flows and increasing Depletion were expected due to the introduction of technologies that increase irrigation efficiency, while climate change could decrease future water availability by an estimated 3–15%. We suggest that to reach the intervention target, strategies need to focus on reducing water depletion rather than water withdrawals. The framework can be used to comprehensively assess water-saving strategies, particularly in water-stressed basins.
Ice nucleating particle concentrations of the past: insights from a
600 year old Greenland ice core
(2020)
Ice nucleating particles (INPs) affect the microphysics in cloud and precipitation processes. Hence, they modulate the radiative properties of clouds. However, atmospheric INP concentrations of the past are basically unknown. Here, we present INP measurements from an ice core in Greenland, which dates back to the year 1370. In total 135 samples were analyzed with the FRIDGE droplet freezing assay in the temperature range from −14 ◦C to −35 ◦C. The sampling frequency was set to 1 in 10 years from 1370 to 1960. From 1960 to 1990 the frequency was increased to 1 sample per year. Additionally, a number of special events were probed, including volcanic episodes. The typical time coverage of a sample was on the order of a few months. Historical atmospheric INP concentrations were estimated with a conversion factor, which depends on the snow accumulation rate of the ice core, particle dry deposition velocity and the wet scavenging ratio. Typical atmospheric INP concentrations were on the order of 0.1 L -1 at −25 ◦C. The INP variability was found to be about 1 – 2 orders of magnitude. Yet, the short-term variability from samples over a seasonal cycle was considerably lower. INP concentrations were significantly correlated to chemical tracers derived from continuous flow analysis (CFA) and ion chromatography (IC) over a broad range of nucleation temperatures. The highest correlation coefficients were found for the particle concentration (dp > 1.2 µm). The correlation is higher for the seasonal samples, where INP concentrations follow a clear annual pattern, highlighting the importance of the annual dust input in Greenland from East Asian deserts during spring. Scanning electron microscopy (SEM) of single particles retrieved from selected samples found particles of soil origin to be the dominant fraction, verifying the significance of mineral dust particles as INPs. Overall, the concentrations compare reasonably well to present day INP concentrations, albeit they are on the lower side. However, we found that the INP concentration at medium supercooled temperatures differed before and after 1960. Average INP concentrations at −23 ◦C, −24 ◦C, −25 ◦C, −26 ◦C and −28 ◦C were significantly higher (and more variable) in the modern day period, which could indicate a potential anthropogenic impact or some post-coring contamination of the topmost, very porous firn.
Ice-nucleating particle concentrations of the past: insights from a 600-year-old Greenland ice core
(2020)
Ice-nucleating particles (INPs) affect the microphysics in cloud and precipitation processes. Hence, they modulate the radiative properties of clouds. However, atmospheric INP concentrations of the past are basically unknown. Here, we present INP measurements from an ice core in Greenland, which dates back to the year 1370. In total 135 samples were analyzed with the FRIDGE droplet freezing assay in the temperature range from −14 to −35 ∘C. The sampling frequency was set to 1 in 10 years from 1370 to 1960. From 1960 to 1990 the frequency was increased to one sample per year. Additionally, a few special events were probed, including volcanic episodes. The typical time coverage of a sample was on the order of a few months. Historical atmospheric INP concentrations were estimated with a conversion factor, which depends on the snow accumulation rate of the ice core, particle dry deposition velocity, and wet scavenging ratio. Typical atmospheric INP concentrations were on the order of 0.1 L−1 at −25 ∘C. The INP variability was found to be about 1–2 orders of magnitude. Yet, the short-term variability from samples over a seasonal cycle was considerably lower. INP concentrations were significantly correlated to some chemical tracers derived from continuous-flow analysis (CFA) and ion chromatography (IC) over a broad range of nucleation temperatures. The highest correlation coefficients were found for the particle concentration (spherical diameter dp > 1.2 µm). The correlation is higher for a time period of seasonal samples, where INP concentrations follow a clear annual pattern, highlighting the importance of the annual dust input in Greenland from East Asian deserts during spring. Scanning electron microscopy (SEM) analysis of selected samples found mineral dust to be the dominant particle fraction, verifying their significance as INPs. Overall, the concentrations compare reasonably well to present-day INP concentrations, albeit they are on the lower side. However, we found that the INP concentration at medium supercooled temperatures differed before and after 1960. Average INP concentrations at −23, −24, −25, −26, and −28 ∘C were significantly higher (and more variable) in the modern-day period, which could indicate a potential anthropogenic impact, e.g., from land-use change.
In this study, we use simulations from seven global vegetation models to provide the first multi‐model estimate of fire impacts on global tree cover and the carbon cycle under current climate and anthropogenic land use conditions, averaged for the years 2001–2012. Fire globally reduces the tree covered area and vegetation carbon storage by 10%. Regionally, the effects are much stronger, up to 20% for certain latitudinal bands, and 17% in savanna regions. Global fire effects on total carbon storage and carbon turnover times are lower with the effect on gross primary productivity (GPP) close to 0. We find the strongest impacts of fire in savanna regions. Climatic conditions in regions with the highest burned area differ from regions with highest absolute fire impact, which are characterized by higher precipitation. Our estimates of fire‐induced vegetation change are lower than previous studies. We attribute these differences to different definitions of vegetation change and effects of anthropogenic land use, which were not considered in previous studies and decreases the impact of fire on tree cover. Accounting for fires significantly improves the spatial patterns of simulated tree cover, which demonstrates the need to represent fire in dynamic vegetation models. Based upon comparisons between models and observations, process understanding and representation in models, we assess a higher confidence in the fire impact on tree cover and vegetation carbon compared to GPP, total carbon storage and turnover times. We have higher confidence in the spatial patterns compared to the global totals of the simulated fire impact. As we used an ensemble of state‐of‐the‐art fire models, including effects of land use and the ensemble median or mean compares better to observational datasets than any individual model, we consider the here presented results to be the current best estimate of global fire effects on ecosystems.
Aim: Plant life‐forms characterize key morphological strategies that enable large‐scale comparisons of plant communities. This study applies Raunkiær's plant life‐form concept that was developed for temperate climate to a subtropical island flora, in parts, dominated by summer aridity. We quantify how plant life‐form patterns as well as patterns of important plant functional traits (PFTs) relate to important climate and topographic characteristics.
Location: La Palma, Canary Islands.
Taxon: Flora of La Palma.
Methods: We assigned each native plant species a plant life‐form, that is, phanerophyte, chamaephyte, hemicryptophyte, geophyte and therophyte, as well as PFTs (succulence and N‐fixer). We used stacked species distribution models to assess occurrence probability for each species using the Atlantis database (500 m × 500 m grid). We related richness and percentage values for each plant life‐form and PFT to climate and topography.
Results: Plant life‐forms and PFTs showed a clear pattern within geographic but also climate space, while topography had a minor effect. Phanerophytes mainly contributed to the flora in humid areas. Chamaephytes and hemicryptophytes most strongly contributed to the summit scrub flora and, to some degree, also to the arid coastal regions. Geophytes and therophytes were mainly found in dry coastal regions. N‐fixers contributed mainly to warm‐arid and cool‐arid regions, while succulent species were mainly found in arid coastal regions.
Main conclusions: Raunkiær's plant life‐form concept can be comprehensively transferred to a subtropical island flora by adapting to local unfavourable growing conditions, that is, aridity. Using the strong environmental gradients offered by our study island, we identify substantial climate‐driven variation in patterns of plant life‐forms and PFTs that might be used for large‐scale comparisons in macroecological studies. The growth strategies reflected in Raunkiær's plant life‐forms suggest differences in species establishment and coexistence dynamics within different parts of the island's climate space.
We apply seismic full waveform inversion to SH‐ and Love‐wave data for investigating the near‐surface lithology at an archaeological site. We evaluate the resolution of the applied full waveform inversion algorithm through ground truthing in the form of an excavation and sediment core studies. Thereby, we investigate the benefits of full waveform inversion in comparison with other established methods of near‐surface prospecting in terms of resolution capabilities and interpretation security. The study is performed in a presumed harbour area of the ancient Thracian city of Ainos. The exemplary target is the source of a linear magnetic anomaly oriented perpendicular to the coast, which was found in a previous magnetic gradiometry survey, suggesting a mole. The SH‐wave full waveform inversion recovered a subsurface SH‐wave velocity model with submeter resolution showing lateral and vertical velocity variation between 40 and 150 m/s. To tame the non‐linearity of the full waveform inversion, a sequential inversion of frequency bands has to be combined with time‐windowing in order to separate the Love wave from the reflected SH wavefield. We compare the full waveform inversion results with multichannel analysis of surface waves, standard seismic reflection imaging, electrical resistivity tomography and electromagnetic induction. It turns out that the respective depth sections are correlated to a certain degree with the full waveform inversion results. However, the structural resolution of the other geophysical methods is significantly lower than for the full waveform inversion. An exception is the reflection seismic imaging, which shows the same resolution as full waveform inversion but can only be interpreted together with the full waveform inversion–based velocity model. An archaeological excavation as well as coring data allows ground truthing and a direct understanding of the geophysical structures. The results show that the target was a sort of near‐surface trench of about 3–4 m width and 0.8 m to 1.0 m depth, filled with silty sediment, which differs from the layered surrounding in colour and composition. The ground truthing revealed that only SH‐wave full waveform inversion and seismic reflection imaging could image the trench and sediment structure with satisfying lateral and depth resolution. We emphasize that the velocity distribution from SH‐wave full waveform inversion agrees closely with the excavated subsurface structures, and that the discovered changes in seismic velocity correlate with changes in the sand content in the respective sediment facies sequences. The study demonstrated that SH‐wave full waveform inversion is capable to image structural and lithological changes in the near subsurface at scales as low as 0.5 m, thus providing the high resolution needed for archaeological and geoarchaeological prospection.
Opportunities and challenges for paleoaltimetry in "small" orogens: insights from the European Alps
(2020)
Many stable isotope paleoaltimetry studies have focused on paleoelevation reconstructions of orogenic plateaus such as the Tibetan or Andean Plateaus. We address the opportunities and challenges of applying stable isotope paleoaltimetry to “smaller” orogens. We do this using a high‐resolution isotope tracking general circulation model (ECHAM5‐wiso) and explore the precipitation δ18O (δ18Op) signal of Cenozoic paleoclimate and topographic change in the European Alps. Results predict a maximum δ18Op change of 4–5‰ (relative to present day) during topographic development of the Alps. This signal of topographic change has the same magnitude as changes in δ18Op values resulting from Pliocene and Last Glacial Maximum global climatic change. Despite the similar magnitude of the isotopic signals resulting from topographic and paleoclimate changes, their spatial patterns across central Europe differ. Our results suggest that an integration of paleoclimate modeling, multiproxy approaches, and low‐elevation reference proxy records distal from an orogen improve topographic reconstructions.
Global‐scale gradient‐based groundwater models are a new endeavor for hydrologists who wish to improve global hydrological models (GHMs). In particular, the integration of such groundwater models into GHMs improves the simulation of water flows between surface water and groundwater and of capillary rise and thus evapotranspiration. Currently, these models are not able to simulate water table depth adequately over the entire globe. Unsatisfactory model performance compared to well observations suggests that a higher spatial resolution is required to better represent the high spatial variability of land surface and groundwater elevations. In this study, we use New Zealand as a testbed and analyze the impacts of spatial resolution on the results of global groundwater models. Steady‐state hydraulic heads simulated by two versions of the global groundwater model G3M, at spatial resolutions of 5 arc‐minutes (9 km) and 30 arc‐seconds (900 m), are compared with observations from the Canterbury region. The output of three other groundwater models with different spatial resolutions is analyzed as well. Considering the spatial distribution of residuals, general patterns of unsatisfactory model performance remain at the higher resolutions, suggesting that an increase in model resolution alone does not fix problems such as the systematic overestimation of hydraulic head. We conclude that (1) a new understanding of how low‐resolution global groundwater models can be evaluated is required, and (2) merely increasing the spatial resolution of global‐scale groundwater models will not improve the simulation of the global freshwater system.
Barriereinseln und Außensände, geformt durch eine Kombination aus Wind, Wellen, Strömung und Küstenlängstransport gelten als morphologisch hoch aktive Küstenbereiche und variieren häufig in Ursprung, Genese und Entwicklung. Sie besitzen durch ihre dissipative Wirkungsweise eine bedeutende Schutzfunktion für rückwärtige Inseln, Halligen und Festlandbereiche und bilden vor der Küste Nordfrieslands die westliche Außengrenze des Wattenmeeres.
Ziel der nachfolgenden Studie ist es, anhand hochauflösender Georadarmessungen und sedimentologischer Daten aus Bohrungen die Landschaftsgeschichte an der Westküste Amrums sowie die Entstehung und interne sedimentäre Architektur des der Insel westlich angelagerten Kniepsandes zu untersuchen und die Prozesse, die zur Genese geführt haben, durch Datierungen zeitlich einzuordnen. Auf Grundlage der gewonnenen Daten wurden zwei stratigraphische Modelle entwickelt, welche die geologisch-geomorphologischen Prozesse und Sedimentationsbedingungen im Westküstenvorfeld erklären.
Ein Modell zeigt die Landschaftsentwicklung an der Westküste Amrums und beschreibt die Sedimentationsbedingungen, die im Vorfeld der Westküste herrschten, bevor der Kniepsand an die Insel heranwanderte. Auf der Landoberfläche des ertrinkenden saaleeiszeitlichen Geestkerns wurden zu Beginn der Flandrischen Transgression feinkörnige Sedimente eines Mischund Schlickwatts abgelagert. Es ist davon auszugehen, dass der damals noch weit vor der Küste Amrums liegende Kniepsand eine Barriere bildete und so an der heute hochenergetischen Westküste für strömungsberuhigte Sedimentationsbedingungen sorgte. Durch Erosion am Geestkern bildeten sich in unterschiedlichen Höhenpositionen fossile Kliffs, die dem damaligen Meeresspiegelstand entsprechen.
Ein weiteres Modell beschreibt den Andockmechanismus des Kniepsandes an die Insel Amrum. Durch die Anlagerung des ehemaligen Außensandes und den damit einhergehenden Sedimentinput wurden die Bedingungen für eine großflächige Dünenbildung geschaffen.
Als erste Ergebnisse der derzeit laufenden paläomagnetischen Untersuchungen im Rhein-Main-Gebiet werden die Befunde von einigen stratigraphisch bedeutenden Profilen mitgeteilt. Während die Remanenzwerte aus altpleistozänen Tonen der Kelsterbacher Terrasse mit einer Ausnahme revers sind, konnte im Autobahneinschnitt Abenheim (NW Worms) und im Lößprofil der Ziegeleigrube Bad Soden am Taunus der Ubergang von reverser zu normaler Magnetisierung gefunden werden. Die Grenze zwischen MATUYAMA- und BRUNHES-Epoche liegt im Profil Bad Soden unter dem 6., das Jaramillo-Event unter dem 7. fossilen Bt-Horizont.
Die über dem letztinterglazialen Boden (Parabraunerde bzw. pseudovergleyte Parabraunerde) folgenden würmzeitlichen äolischen Sedimente können durch kennzeichnende Bodenhorizonte in 3 Abschnitte (Alt-, Mittel- und Jungwürm) gegliedert werden. Charakteristisch für das Altwürm sind Humuszonen, für das Mittelwürm neben einigen „Naßböden" braune Verlehmungszonen bis zu einer Mächtigkeit von 1,1 m und für das Jungwürm mehrere schwach ausgebildete, geringmächtige Verbraunungszonen und „Naßböden". Den wichtigsten Leithorizont des Jungwürms bildet das Kärlicher Tuffbändchen, das in jüngster Zeit auch in Nordhessen aufgefunden wurde. Abschließend wird das in Hessen aufgestellte Gliederungsschema mit den Würmlöß-Gliederungen in anderen Teilen Europas verglichen und eine Parallelisierung versucht.
Der Auffassung, das Substrat der Lockerbraunerde auf dem Oberwald des Vogelsberges sei ausschließlich ein äolisches Sediment der Jüngeren Tundrenzeit, wird widersprochen. Verschiedene Befunde, vor allem 14C-Datierungen, lassen vielmehr den Schluß zu, daß große Teile des Substrates vieler Lockerbraunerden holozäne anthropogene Kolluvien sind. Bekräftigt wird dagegen der Befund gleichen (jungtundrenzeitlichen) Alters für Lockerbraunerde-Substrat und Deckschutt (Decksediment).
Aus dem Gebiet der weichselzeitlichen Vereisung in Polen werden allgemein verbreitete geringmächtige periglaziale Deckschichten beschrieben. Sie zeigen in der Regel eine äolische Beeinflussung und unterscheiden sich dadurch vom Liegenden. Es handelt sich um spätglaziale Bildungen, wie sie aus dem Jungmoränengebiet der DDR seit langem bekannt sind. Ähnliche Substrate wurden auch im nördlichen Alpenvorland gefunden.
Im Lößprofil der Ziegelei Glos bei Lisieux ist eine Abfolge fossiler Boden aufgeschlossen, die sich gut mit der hessischen Lößstratigraphie verbinden läßt. Über dem letztinterglazialen Boden liegen nach einer Diskordanz der Lohner Boden, der E1-, der E2- und der E4-Naß-boden. Letzterer wird häufig als Äquivalent des „Sol de Kesselt" angesehen. Somit zeigt sich auch für dieses Gebiet, daß dieser Boden stratigraphisch nicht dem Lohner Boden entsprechen kann.
Un tuf volcanique a ete mis en evidence dans trois profils de loess de Hesbaye (Rocourt et Lixhe) et du Limbourg neerlandais (Nagelbeek). Ce tuf correspond probablement ä l'Eltviller Tuff connu en Hesse (SEMMEL 1967). En consequence, cette decouverte permet de comparer directement une partie de la sequence des loess du Pleistocene superieur propre ä chaque region
In der Monte Cavallo-Gruppe fand Verf. in Stauseesedimenten fossile Holzstücke, für die die 14 C-Datierung ein Alter von 29 350 ± 460 Jahren vor 1950 n. Chr. ergab. Die schluffigen Ablagerungen, die in etwa 900 m, fast 80 m über dem heutigen Flußbett des T. Caltea aufgeschlossen sind, ruhen auf Schottern und werden von Moräne überlagert.
Aufgrund seiner stratigraphischen Lage und seines 14C-datierten Alters kann das Holz (Picea abies bzw. Larix) einem zeitlichen Äquivalent des Paudorf-Interstadials zugeordnet werden. Dieser Fund von Großresten ermöglicht somit eine erste absolute Datierung des Paudorf-Interstadials für die Südalpen und Oberitalien.
Das Cranium eines fossilen Hominiden des Formenkreises Homo sapiens sapiens wurde relativ-geologisch sowie absolut durch Radiokohlenstoff und Aminosäuren auf ungefähr 31 000 Jahre B.P. datiert. Andere absolute sowie relative Daten an Mollusken und Mammutzähnen in überlagernden jüngeren Straten datieren auf 18 000 — 21000 und 16 000 Jahren B.P. Geomorphologische und geophysikalische Datierungen stimmen somit gut überein. Er ist der älteste datierte und früheste Bewohner Zentraleuropas, der dem Homo sapiens sapiens angehört.
Aus Franken wird die Entwicklung quartärer Hohlformen beschrieben, deren Rekonstruktion mit Hilfe lößstratigraphischer Methoden (fossile Böden, Tuffbänder, Umlagerungszonen etc.) möglich ist. Bei vielen Formen zeigt sich, daß sie bereits größere Vorläuferformen präwürmzeitlichen Alters hatten. Die Entwicklung während des Würms läßt sich an manchen Beispielen in besonders instruktiver Weise verfolgen. Zu Beginn des Würms, im unteren Mittelwürm und im unteren Jungwürm dominierte zeitweise die Abtragung und Verlagerung. Im oberen Mittelwürm sowie im oberen Jungwürm herrschte äolische Lößsedimentation vor. Diese Ergebnisse stimmen gut mit den bereits aus anderen mitteleuropäischen Lößgebieten bekannten Befunden überein. Mit dem Trockental-System von Helmstadt wird die Entwicklung von Hohlformen beschrieben, deren Anlage bis in das ältere Pleistozän zurückreicht.
In vorliegender Untersuchung wurde der Rißlöß zwischen der 1. und 2. fossilen Parabraunerde anhand schwächerer Bodenbildungen und eingeschalteter Abtragungsphasen zu gliedern versucht. Im jüngeren Riß herrschte starke Lößsedimentation vor, wobei es in mindestens 6 kaltfeuchten Abschnitten zur Ausbildung schwacher periglazialer Naßböden kam. Die Naßbodenserie wurde als Bruchköbeler Böden (B) bezeichnet. Im jüngsten Rißlöß ist wenige dm unter dem Eemboden als tephrochronologischer Leithorizont der Krifteler Tuff (vgl. SEMMEL 1968) eingeschaltet. Den mittleren Profilbereich im Rißlöß zeichnen feuchtere Klimaabschnitte mit starken Verschwemmungsphasen aus, die in den meisten Profilen zu erheblichen Diskordanzen geführt haben. An der Basis der wenigen kompletten Rißlöß-Profile treten vorwiegend in Hessen über der zumeist gekappten 2. fossilen Parabraunerde maximal zwei Schwarzerden auf, die von SEMMEL (1968) als Weilbacher Humuszonen bezeichnet werden. Unmittelbar über diesen Schwarzerden folgt die Ostheimer Zone, eine Fließerde aus aufgearbeitetem Solumaterial der liegenden Böden. Insgesamt zeigt die aus den Rißböden rekonstruierte Klimaabfolge — neben geringfügigen Abweichungen — überraschende Parallelen zur paläopedologisch-klimatischen Gliederung der Würmkaltzeit.
Crustal structure at the western end of the North Anatolian Fault Zone from deep seismic sounding
(2001)
The first deep seismic sounding experiment in Northwestern Anatolia was carried out in October 1991 as part of the "German - Turkish Project on Earthquake Prediction Research" in the Mudurnu area of the North Anatolian Fault Zone. The experiment was a joint enterprise by the Institute of Meteorology and Geophysics of Frankfurt University, the Earthquake Research Institute (ERI) in Ankara, and the Turkish Oil Company (TPAO). Two orthogonal profiles, each 120 km in length with a crossing point near Akyazi, were covered in succession by 30 short period tape recording seismograph stations with 2 km station spacing. 12 shots, with charge sizes between 100 and 250 kg, were fired and 342 seismograms out of 360 were used for evaluation. By coincidence an M b = 4.5 earthquake located below Imroz Island was also recorded and provided additional information on Moho and the sub-Moho velocity. A ray tracing method orginally developed by Weber (1986) was used for travel time inversion. From a compilation of all data two generalized crustal models were derived, one with velocity gradients within the layers and one with constant layer velocities. The latter consists of a sediment cover of about 2 km with V p » 3.6 km/s, an upper crystalline crust down to 13 km with V p » 5.9 km/s, a middle crust down to 25 km depth with V p » 6.5 km/s, a lower crust down to 39 km Moho depth with V p » 7.0 km/s and V p » 8.05 km/s below the Moho. The structure of the individual profiles differs slightly. The thickest sediment cover is reached in the Izmit-Sapanca-trough and in the Akyazi basin. Of particular interest is a step of about 4 km in the lower crust near Lake Sapanca and probably an even larger one in the Moho (derived from the Imroz earthquake data). After the catastrophic earthquake of Izmit on 17 August 1999 this significant heterogeneity in crustal structure appears in a new light with regard to the possible cause of the Izmit earthquake. Heterogeneities in structure are frequently also heterogeneities in strength and stress that impede or even lock rupture. The Izmit earthquake is discussed in relation to a large stepover or jog at the North Anatolian Fault.
Aim: Predicting future changes in species richness in response to climate change is one of the key challenges in biogeography and conservation ecology. Stacked species distribution models (S‐SDMs) are a commonly used tool to predict current and future species richness. Macroecological models (MEMs), regression models with species richness as response variable, are a less computationally intensive alternative to S‐SDMs. Here, we aim to compare the results of two model types (S‐SDMS and MEMs), for the first time for more than 14,000 species across multiple taxa globally, and to trace the uncertainty in future predictions back to the input data and modelling approach used.
Location: Global land, excluding Antarctica.
Taxon: Amphibians, birds and mammals.
Methods: We fitted S‐SDMs and MEMs using a consistent set of bioclimatic variables and model algorithms and conducted species richness predictions under current and future conditions. For the latter, we used four general circulation models (GCMs) under two representative concentration pathways (RCP2.6 and RCP6.0). Predicted species richness was compared between S‐SDMs and MEMs and for current conditions also to extent‐of‐occurrence (EOO) species richness patterns. For future predictions, we quantified the variance in predicted species richness patterns explained by the choice of model type, model algorithm and GCM using hierarchical cluster analysis and variance partitioning.
Results: Under current conditions, species richness predictions from MEMs and S‐SDMs were strongly correlated with EOO‐based species richness. However, both model types over‐predicted areas with low and under‐predicted areas with high species richness. Outputs from MEMs and S‐SDMs were also highly correlated among each other under current and future conditions. The variance between future predictions was mostly explained by model type.
Main conclusions: Both model types were able to reproduce EOO‐based patterns in global terrestrial vertebrate richness, but produce less collinear predictions of future species richness. Model type by far contributes to most of the variation in the different future species richness predictions, indicating that the two model types should not be used interchangeably. Nevertheless, both model types have their justification, as MEMs can also include species with a restricted range, whereas S‐SDMs are useful for looking at potential species‐specific responses.
This study reports and discusses new radiometric ages, petrographical and volcanological observations and whole rock geochemical data of the rocks of the Rudnik Mts. volcano-intrusive complex. The complex hosts a Pb-Zn-Ag deposit and belongs to the Serbo-Macedonian metallogenetic belt. Two distinct igneous events are distignuished. The first occurred >30 Ma and was mainly characterized by extrusive and shallow intrusive dacites and andesites and was unrelated to mineralization. The second igneous event occurred <23 Ma and was highly heterogeneous in terms of volcanic products and petrographic varieties, but with predominance of quartzlatites. The dacite-andesites (first event) and the quartzlatites (second event) are geochemically similar and display a calc-alkaline affinity and highly incompatible element enriched patterns on spider diagrams, but the younger quartzlatites are richer in K2O, Rb and Ba and poorer in Sr. This is taken as evidence that mixing between an ultrapotassic lamprophyre/lamproite magma and an acid calc-alkaline (dacite-like) magma was essential petrogenetic processes during the second event. The proposed simplified volcanological model suggests that this mixing was responsible for triggering strongly explosive volcanic activity as well as for providing conditions for active hydrothermal and mineralization processes. The observed link between a specific magmatic phase and ore deposit formation can be a general phennomenon in the Balkans, and must be addressed by further and more advanced studies.