Geosciences, Volume 12, Issue 7 (July 2022) – 27 articles

The majority of formulae for predicting in-channel mixing by longitudinal dispersion are based on empirical evidence from medium to large rivers, whereas small creeks and streams are under-represented despite their hydrological and ecological importance. In this study, twenty-six formulae for predicting the longitudinal dispersion coefficient (K_x) were evaluated for their applicability to small channels using field measurements and hydraulic modeling. Predicted values for K_x, following guidelines recommended in the original publications for the formulae, were compared to measured dispersion coefficients obtained from sodium-chloride plumes injected into two small channels (a concrete-lined, trapezoidal channel and a small, natural stream) based on fitting the Hayami solution to the one-dimensional advection-dispersion equation (ADE) to each plume. The predicted K_x coefficients from the formulae were also utilized to create model-simulated plumes, which were compared to those measured well downstream of the point of injection. The findings demonstrate that the predictive accuracy of the twenty-six formulae was extremely variable; none were able to predict the dispersion process in the small channels with better than ± 50% accuracy. These results show that “universal” formulae are plagued with a large degree of uncertainty and should be used with caution when applied to small channels, although more robust predictions are possible with some formulae if site-specific data are available for calibration. Full article

The objective of this work is to analyze the spatial-temporal features of the trends in the frequency and amount of erosion-hazardous precipitation in the European part of Russia (EPR) for the period 1966–2020, as a reflection of the influence of climatic changes on surface runoff from the cultivated slopes during the warm season. One hundred and fifty-nine EPR weather stations were selected for analysis based on the length of the time series and the amount of missing data. Several characteristics of erosion-hazardous precipitation were considered: the number of days with a daily precipitation of more than 12.7 mm, the number of days with a daily precipitation of 12.7 to 40 mm, the number of days with a daily precipitation of more than 40 mm, the maximum one-day precipitation. In general, it can be stated that even within the southern taiga, mixed forests, and forest-steppe ecoregion (broad-leaved forests), within which a positive increase in the frequency of erosion-hazardous precipitation was detected, there was no significant increase in the rate of washout and linear washout, which is primarily due to a more significant reduction of slope runoff and soil washout during spring snowmelt. Precipitation, the daily amount of which is more than 40 mm, as well as the maximum daily amount of precipitation, show an upward trend in the western contact zone of mixed forests and forest-steppe, on the Black Sea coast, as well as in the northern foothills of the Caucasus, where their contribution to erosion processes is likely to increase against the decrease in the number of days with precipitation of a 12.7–40 mm daily amount. Full article

Gravitational mass movements such as rockfalls, landslides, rock avalanches, or debris flows are increasingly endangering settlement areas and infrastructure facilities in the Alpine region as a result of climate change. An essential component of counteracting the dangers of such events is the construction of suitable protective structures. However, the dimensioning of these protective structures requires in-depth knowledge of the impact process on the structure. Measurements of real large mass movements such as rock avalanches fail due to the large impact forces involved. For this reason, model tests have been carried out by different institutions in different countries in recent decades. An essential aspect of the study of gravitational mass movements using model experiments is scaling experimental results to real events. Therefore, in this study, a model experiment carried out at the University of Innsbruck was recalculated in the first step using the discrete element method (DEM). Subsequently, the experimental results and the numerical DEM model were scaled to a real event using scale factors and then compared again. The aim was to show how well the results of the model tests can be scaled to describe real events of rock avalanches. Full article

The coastal area of the Kavala district, Northern Greece, is characterized by minerals enriched in rare earth elements (REE). The present study focuses on the mineralogy of the black sands from six different locations and the comprehensive mineral chemistry of the REE-bearing minerals, allanite-(Ce), epidote, monazite, thorite, zircon and titanite. Allanite-(Ce) is the most important carrier of light REE (LREE) in the studied black sands, reaching up to 23.24 wt % ΣREE. The crystal chemistry of allanite-(Ce) transitions into ferriallanite-(Ce), due to the significant involvement of Fe³⁺. High resolution backscattered electron (BSE) images were used to identify zoning that corresponds to variations in REE, Th and U. These modifications follow the exchange scheme: (Ca + (Fe³⁺, Al))₋₁(LREE, Y, Th, U + (Fe²⁺, Mg, Mn))₊₁. Epidotes may also contain up to 0.5 REE³⁺ apfu. Monazite and thorite are found as inclusions in allanite-(Ce) and are enriched in Ce, La and Nd, together with Th and U. Some zircons are enriched in Hf, while some titanites host Nb and V. Full article

Migration of dissolved salts from natural (cryopeg brines, seawater, etc.), or artificial sources can destabilize intrapermafrost gas hydrates. Salt transport patterns vary as a function of gas pressure, temperature, salinity, etc. The sensitivity of the salt migration and hydrate dissociation processes to ambient temperature and to the concentration and chemistry of saline solutions is investigated experimentally on frozen sand samples at a constant negative temperature (−6 °C). The experiments show that the ambient temperature and the solution chemistry control the critical salt concentration required for complete gas hydrate dissociation. Salt ions migrate faster from more saline solutions at higher temperatures, and the pore moisture can reach the critical salinity in a shorter time. The flux density and contents of different salt ions transported to the samples increase in the series Na₂SO₄–KCl–CaCl₂–NaCl–MgCl₂. A model is suggested to account for phase transitions of pore moisture in frozen hydrate-saturated sediments exposed to contact with concentrated saline solutions at pressures above and below the thermodynamic equilibrium, in stable and metastable conditions of gas hydrates, respectively. Full article

Pompeii Archaeological Park is the best laboratory for the study of the seismic site effects on cultural heritage: the ancient site was destroyed and buried by the 79 AD Vesuvian eruption and, furthermore, it was also affected by the 62–63 AD strong earthquake. Large sectors of the city were reconstructed after this earthquake while other parts were still under reconstruction when the fall-out and pyroclastic density currents of the eruption buried the Roman city. In order to evaluate the distribution of the damage and reconstructions due to the earthquake, detailed mappings of the structures were carried out using multidisciplinary approaches. In addition, analyses of the topographical features, subsoil stratigraphies, and geophysical surveys, responsible for local seismic amplification (site effects), allow us to define the sectors of the ancient city where the Amplification Factors (AFs) were the main ones responsible for damage. Selected areas and examples of compromised and reconstructed buildings show that the ancient topography and subsoil features (both lithological and seismic) are the main AFs. In particular, the damages caused by the 62–63 AD earthquakes seem to be mainly due to topographical factors such as steep scarps and slopes, ridges, peaks, and terraces, as well as to the major thickness of the soft sediments (loose volcanoclastic layers, paleosols, weathered lavas, and anthropogenic infillings) located over the well-lithified lavas. It is not uncommon to also have the combination of these two factors. For the first time, this multidisciplinary approach allows us to draw a seismic microzonation map for one of the most important archaeological sites of the world. Full article

Like many other lakes in the world, the interconnected Abaya and Chamo lakes in the Southern Main Ethiopian Rift are affected by rapid sediment accumulation. Although land degradation is a well-known issue in this part of the African continent, the main sediment sources, their spatial distribution and interaction in the Abaya–Chamo lakes’ basin have not yet been documented. Here, we present a systematic inventory, characterization, and spatial analysis of landslides and gullies as concentrated sediment sources, for four representative river catchments impacted by landscape rejuvenation. Using Google Earth imagery and field surveys, we mapped with high accuracy a total of 7336 gullies and 430 landslides. Recent landslides observed during the last decade were few, small and shallow, and appear to have played a minor role in the current sediment dynamics. Large landslides are old and inactive. Although they do not contribute to the current sediment budget, they contribute indirectly to landscape dynamics by favoring the occurrence of gullies. Overall, large percentages of severe to extremely degraded areas of gully erosion are located in rejuvenated landscapes, especially at the level of the old landslides. Many active gullies are connected to the river network, as such acting as the source of sediment. Our analysis is a step towards understanding the nature and control of anthropic activities on sediment production in the region. We also highlight the importance of considering the interactions between sediment sources and the connectivity of the geomorphological system. Full article

This study examines the mineralogy, texture, and weathering grades of intrusive igneous rocks from northern Pakistan, as well as their impacts on physical and strength properties. In comparison to felsic rocks, mafic and intermediate rocks have lower cumulative proportions of quartz, feldspar, and plagioclase, as well as higher specific gravity, strength (i.e., UCS and R-value), and UPV values. Similarly, samples with anhedral grain shapes, irregular boundaries, and fine to medium grain sizes (UD, ANS, and CGN) exhibited greater strength values, with compressive strengths of 121, 118, and 91 MPa and tensile strengths of 11, 9, and 12 MPa, respectively. The physical and strength properties of the investigated samples corresponded well with the weathering grades assigned to them, such as fresh (WG-I), slightly weathered (WG-II), and highly weathered (WG-III). That is, as the grade increased from WG-I to WG-III, the porosity and water absorption increased (0.28% and 0.72%, respectively), whereas the specific gravity, compressive strength, and tensile strength decreased (2.04, 20, and 2.5 MPa, respectively, for CGA). Although the presence of quartz impacts rock strength, no significant association was found between the strength and the maximum and mean grain sizes of other minerals. Full article

Sedimentary systems developed in small (<1 km²) mountain lacustrine basins represent high-resolution geological archives, able to record subtle climatic and tectonic signatures over historical times. The studied example from the Mt. Ragola ophiolitic massif in the Northern Apennines (Italy) allowed us to better understand the role of the different (neo)tectonic and climatic events on the development and distribution of large landslides and lakes/peat bogs during the last 10 kyrs. Implementing a multidisciplinary approach that includes detailed acquisition of bedrock, geomorphological, topographic, and geophysical data, we detected and mapped ridge splitting, trenches, closed depressions, double ridges, and counterscarps. These morphostructures are interpreted as relevant factors influencing the distribution of sediments in historical times by shifting the position of the local equilibrium point (i.e., erosion vs. deposition) and have been correlated to a combination of climatic (i.e., increased flood events) and tectonic (i.e., spatial–temporal clustering of seismic shocks) forcing, starting from the demise of the Little Ice Age in the mid-19th century to the present-day situation. This approach allowed us to better describe the current changes in the hydrologic cycle, reaching beyond the limits of historical instrumental data. Furthermore, allowing the recognition and dating of recent tectonic vs. gravitational deformations, it also constitutes an integrative method for assessing the local geological hazard. Full article

We present initial rock magnetic results for both lava flows and lapilli produced by the 2021 eruption of the Cumbre Vieja, La Palma (Canary Islands). Samples were taken during the eruption to minimize early alteration and weathering of the rocks and tephra. Standard procedures included progressive alternating field and thermal demagnetization, hysteresis curves, thermomagnetic experiments, progressive acquisition of isothermal remanent magnetization (IRM), and First-Order Reversal Curves (FORCs). Overall, our observations, including low to medium unblocking temperatures, isothermal remanent magnetization to 1 Tesla, and the abundance of wasp-waist hysteresis loops, strongly suggest the presence of Ti-rich titanomagnetites as the main remanence carriers in both lava flows and lapilli, in addition to some hematite as well. Whereas the former has been directly seen (SEM), hematite is elusive with nonmagnetic-based methods. Rock magnetic data, on a Day plot, also reveal that the magnetic grain size tends to be larger in the lava flows than in the lapilli. Full article

Seismic interferometry is often proposed as a cost-efficient technique for reservoir monitoring including CO₂ sequestration due to its low cost and environmental advantages over active source imaging. Although many studies have demonstrated the ability of seismic interferometry to retrieve surface waves, body wave imaging remains challenging due to their generally lower amplitudes of body waves in seismic interferometry data. An optimum data acquisition strategy can help retrieve low amplitude body waves better, however, rare attempts have been made to evaluate various data acquisition strategies. In this study, we use numerical modeling to examine three different acquisition schemes to evaluate the retrievability of P waves from seismic interferometry data. From our numerical results, we observe that (1) positing receivers beneath the attenuated weathered layer improves the data quality and signal to noise ratio, but additional processing steps including predictive deconvolution and Radom transform filter are necessary to remove the downgoing surface multiples, artifacts that are generated from this data acquisition; (2) vertical seismic profiling (VSP) alongside with the conventional surface seismic acquisition improve the target zone detection; and (3) crosswell acquisition of seismic interferometry is an ineffective means to obtain reflection events due to the non-similarity of ray paths from the noise sources meaning that the required stationary phase theory is not fulfilled. Full article

When conducting numerical upscaling, either for a fractured or a porous medium, it is important to account for anisotropy because in general, the resulting upscaled conductivity is anisotropic. Measurements made at different scales also demonstrate the existence of anisotropy of hydraulic conductivity. At the “microscopic” scale, the anisotropy results from the preferential flatness of grains, presence of shale, or variation of grain size in successive laminations. At a larger scale, the anisotropy results from preferential orientation of highly conductive geological features (channels, fracture families) or alternations of high and low conductive features (stratification, bedding, crossbedding). Previous surveys of homogenization techniques demonstrate that a wide variety of approaches exists to define and calculate the equivalent conductivity tensor. Consequently, the resulting equivalent conductivities obtained by these different methods are not necessarily equal, and they do not have the same mathematical properties (some are symmetric, others are not, for example). We present an overview of different techniques allowing a quantitative evaluation of the anisotropic equivalent conductivity for heterogeneous porous media, via numerical simulations and, in some cases, analytical approaches. New approaches to equivalent permeability are proposed for heterogeneous media, as well as discontinuous (composite) media, and also some extensions to 2D fractured networks. One of the main focuses of the paper is to explore the relations between these various definitions and the resulting properties of the anisotropic equivalent conductivity, such as tensorial or non-tensorial behavior of the anisotropic conductivity; symmetry and positiveness of the conductivity tensor (or not); dual conductivity/resistivity tensors; continuity and robustness of equivalent conductivity with respect to domain geometry and boundary conditions. In this paper, we emphasize some of the implications of the different approaches for the resulting equivalent permeabilities. Full article

Faut gouge forms at the core of the fault as the result of a slip in the upper brittle crust. Therefore, the deformation mechanisms and conditions under which the fault gouge was formed can document the stages of fault movement in the crust. We carried out a microstructural analysis on a fault gouge from a hanging-wall branch fault of the Simplon Fault Zone, a major low-angle normal fault in the European Alps. We use thin-section analysis, together with backscattered electron imaging and X-ray diffractometry (XRD), to show that a multistage history from ductile to brittle deformation within the fault gouge. We argue that this multistage deformation history is the result of continuous exhumation history from high to low temperature, along the Simplon Fault Zone. Because of the predominance of pressure solution and veining, we associated a large part of the deformation in the fault gouge with viscous-frictional behaviour that occurred at the brittle-ductile transition. Phyllosilicates and graphite likely caused fault lubrication that we suggested played a role in localizing slip along this major low-angle normal fault. Full article

At the hydrothermal stage, volcanoes are affected by internal and external processes that control their fumarolic and eruptive activity. Monitoring hydrothermal activity is challenging given the diverse nature of the processes accounting for deeper magmatic and shallow hydrothermal sources. A better understanding of these processes has commonly been achieved by combining geochemical and geophysical techniques. However, existing geochemical techniques only include the surveillance of major gas components in fumarolic emissions or major ions in cold/thermal springs. This work presents a long-term (2017–2021) surveillance of major and trace elements in fumarolic condensates from the Cratère Sud vent, a low-temperature steam-rich emission from the La Soufriere de Guadeloupe volcano. This fumarole presented a fluctuating performance, offering a unique opportunity to reveal the behavior of major and trace elements, as well as the physicochemical processes affecting magmatic and hydrothermal sources. Time-series analyses allowed for the identification of pH-related chemical fluctuations associated with (1) variable inputs of deep magmatic components at the root of the hydrothermal system, (2) pressurization episodes of the hydrothermal system with increasing fluid–rock interaction, acid gas scrubbing, and vapor scavenging of metals, and (3) the decreased hydrothermal activity, decreasing scrubbing efficiency. Variations in the volatile content (e.g., S, Sb, B, Cl, Bi, Zn, Mo, Br, Cd, Ag, Cu, and Pb), the amount of leached rock-related elements (e.g., Na, Mg, Al, Si, P, K, Ca, Ti, Cr, Mn, Fe, Rb, Sr, Y, Cs, Ba, REEs, and U), and variations in the concentration of Cl and S alone, are postulated as key parameters to monitor volcanic–hydrothermal systems in unrest, such as La Soufriere. Our results demonstrate that monitoring using condensates is a useful geochemical technique, complementing conventional methods, such as “Giggenbach” soda flasks or the so-called Multigas. Full article

The El Salvador Fault Zone (ESFZ) is part of the Central American Volcanic Arc and accommodates the oblique separation movement between the forearc sliver and the Chortis block (Caribbean Plate). In this work, a triclinic transtension model was applied to geological (fault-slip inversion, shape of volcanic calderas), seismic (focal mechanisms) and geodetic (GPS displacements) data to evaluate the characteristics of the last stages of the kinematic evolution of the arc. The El Salvador Fault Zone constitutes a large band of transtensional deformation whose direction varies between N90° E and N110° E. Its dip is about 70° S because it comes from the reactivation of a previous extensional stage. A protocol consisting of three successive steps was followed to compare the predictions of the model with the natural data. The results show a simple shear direction plunging between 20° and 50° W (triclinic flow) and a kinematic vorticity number that is mostly higher than 0.81 (simple-shearing-dominated flow). The direction of shortening of the coaxial component would be located according to the dip of the deformation band. It was concluded that this type of analytical model could be very useful in the kinematic study of active volcanic arcs, even though only information on small deformation increments is available. Full article

Multi-parameter observations, powerful for the study of lithosphere–atmosphere–ionosphere coupling (LAIC), have been performed for a recent Tokyo earthquake (EQ) with a moderate magnitude (M = 5.9) and rather larger depth (~70 km) on 7 October 2021, in the hope of predicting the next Kanto (Tokyo) huge EQ, such as the 1923 Great Kanto EQ (with a magnitude greater than 7). Various possible precursors have been searched during the two-month period of 1 September to 31 October 2021, based on different kinds of data sets: (i) ULF (ultra-low frequency) magnetic data from Kakioka, Japan, (ii) ULF/ELF (extremely low frequency) magnetic field data from the Chubu University network, (iii) meteorological data (temperature and humidity) from the Japan Meteorological Agency (JMA), (iv) AGW (atmospheric gravity wave) ERA5 data provided by the European Centre for Medium-Range Weather Forecast (ECMWF), (v) subionospheric VLF/LF (very low frequency/low frequency) data from Russia and Japan, (vi) ionosonde Japanese data, and (vii) GIM (global ionosphere map) TEC (total electron content) data. After extensive analyses of all of the above data, we have found that there are a few obvious precursors: (i) ULF/ELF electromagnetic radiation in the atmosphere, and (ii) lower ionospheric perturbations (with two independent tools from the ULF depression and subionospheric VLF anomaly) which took place just two days before the EQ. Further, ULF/ELF atmospheric electromagnetic radiation has been observed from approximately one week before the EQ until a few days after the EQ, which seems to be approximately synchronous in time to the anomalous variation in meteorological parameters (a combination of temperature and humidity, atmospheric chemical potential). On the other hand, there have been no clear anomalies detected in the stratospheric AGW activity, and in the NmF2 and TEC data for the upper F region ionosphere. So, it seems that the lithospheric origin is not strong enough to perturb the upper F region. Finally, we discuss the possible hypothesis for the LAIC process, and we can conclude that the AGW hypothesis might be ruled out, but other possible channels such as the chemical channel (radon emanation) and the associated effects might be in operation, at least, for this Tokyo EQ. Full article

Geoheritage is a geotourism resource that could lead to the sustainable development of society, and could contribute to information on geological and geomorphological characteristics, education, public awareness on geoconservation, and landscape protection. Understanding this cultural wealth is still limited, since the information available is often specialized. Communication addressed to the emotional sphere can profoundly engage people, and technical information can be transmitted more effectively. For several years, our research group dealt with geoheritage and geotourism in central Italy. This work aimed to investigate different communication methods, such as publications, announcements through social media, conferences, live performances, and field excursions aimed at the emotional involvement of the public. Furthermore, the emotional response of the public was explored through a digital survey. These communication methodologies enabled a broad public reach, without any distinction based on social origin, cultural background, or age. Narrating the landscape through different forms of art creates a personal relationship with the audience. The emotional participation of the people demonstrates the effectiveness of the method in transmitting scientific and cultural themes, such as the origin of the geological landforms, the culture and the traditions that the site generated, and the problems pertaining to its protection and conservation. Full article

In recent years, it has become a common understanding among researchers that one of the significant agents for forming a variety of submarine geomorphologies is turbidity currents that travel long distances, so-called long-runout turbidity currents. In this study, we present a mathematical model of the formation of cyclic steps due to long-runout turbidity currents based on the assumption that the highly concentrated layer near the bottom of the turbidity currents achieves an equilibrium state. The model uses the four governing equations: the momentum equation, the continuity equation, the diffusion/dispersion equation of suspended sediment, and the continuity equation of sediment on the bed (the Exner equation). We simplify the governing equations by ignoring the entrainment from ambient water since there is a distinct density interface above the highly concentrated lower layer of long-runout turbidity currents. We determine the dimensions of cyclic steps based on the solution. Agreement between the predictions and field observations turns out to be reasonable. Full article

The generalized Hoek-Brown (GHB) failure criterion can estimate the rock mass parameters required for rock mechanics–related analyses such as numerical modeling in geomechanics. The determination of GHB parameters has been developed in the field of rock mechanics. Due to the wide use of the Mohr-Coulomb criterion and the lack of an existing relationship for determining its parameters for a rock mass, equivalent Mohr-Coulomb parameters (EMC) can be derived from the GHB. To determine the differences in the use of these two criteria, we analyzed the behavior of a deep circular tunnel in nine stress states for three metamorphic rocks recovered from the Canadian Shield from rock masses that present a very blocky structure. We carried out 241 simulations using the finite element code RS2 to assess the effect of the geological strength index (GSI), in situ stress, and rock type on the deviation of wall displacement, the number of yielded elements, and the differential stress obtained by the GHB and EMC parameters. A combination of low in situ stress and high GSI yielded similar results when using both failure criteria. Full article

Salt migration may be another reason why pore-gas hydrates dissociate in permafrost, besides pressure and temperature changes. Temperature variations in frozen hydrate-saturated sediments interacting with a NaCl solution have been studied experimentally at a constant temperature, ~−6 °C typical for permafrost. The experiments with frozen sandy samples containing metastable methane hydrate show that the migration of Na+ ions in the NaCl solution and their accumulation in the sediments can induce heat-consuming hydrate dissociation and ice melting. The hydrate-saturated frozen soils cool down at higher rates than their hydrate-free counterparts and require more time to recover their initial temperature. The temperature effects in hydrate-saturated frozen sediments exposed to contact with NaCl solutions depend strongly on salt concentration. The experimental results are used to model phase changes in the pore space associated with salt-ions transport and provide insights into the reasons for temperature changes. Full article

Flood risk in Italy is a key aspect for the administrative authorities, from the national to the local level. This is especially true in Northern Italy, where the Po River, the most important river of the peninsula, and its river basin are located. In North-Western Italy, the Po Basin is described by numerous sub-basins, among which is the Tanaro River basin: here, in the last decades, floods have produced significant damage, causing an increased concern to local and regional administrations. The main goal of this study was to identify suitable sites for the construction of dams, having the function of retention basins, aiming to mitigate the flood risk in the Upper Tanaro Valley. First, using a qualitative approach, suitable sites were identified using available public data provided by regional administrations and field data obtained from geomorphological surveys, later elaborated in a Geographic Information System (GIS) environment. Several models were then produced using conventional methods to evaluate the hydrological characteristics of the study area and to assess the efficiency of each site in terms of flood water flow rate reduction: the performance was evaluated at control sections chosen in urban areas along the Upper Tanaro Valley. The results show that it is possible to find suitable locations for risk mitigation structures. These models also allowed for a rapid cost-effectiveness evaluation, which led to the definition of the best-performing site. The Upper Tanaro Valley case study here analyzed contributes to proving the importance of an integrated approach based on geomorphological, geo-hydrological, and hydraulic evaluations when dealing with the choice of a flood risk mitigation strategy. Full article

The distinctive features and fossil content of some caves from eastern Sicily (San Teodoro, Donnavilla, Fulco, Taormina, Tremilia, Spinagallo), altogether spanning from the middle Pleistocene until the beginning of the Holocene, are discussed. Although dating on vertebrate and/or invertebrate remains is available in few instances, coastal notches and marine terraces correlate with the caves, provide further chronological constraint. The San Teodoro and Spinagallo caves are the best known, whereas the Tremilia cave deserves to be better analysed. Most caves, but not the San Teodoro one (including only terrestrial faunas), testify to the transition from submarine coastal environments (documented by biogenic crusts, borings, shelly sediments), to continental conditions (vertebrate remains of the Paleoloxodon falconeri, Maccagnone and San Teodoro Faunistic Complexes). The fossil register preserved in these cavities represents a source of information useful to (1) reconstruct the palaeogeography of Sicily and its coastline, largely resulting from the interplay between tectonic and sea-level changes linked to climate fluctuations; and (2) the consequent evolution of the terrestrial biota, including the dominance of insular endemic taxa later replaced by species shared with continental Italy, after the establishment of temporary connection through the Messina Strait. Full article

Barhal Valley belongs to the Çoruh Valley System in the Kaçkar Mountains of northeastern Anatolia. This 13 km long valley is located to the south of the main weather divide and to the east of Mt. Kaçkar, with the highest peak of the mountain range being 3932 m. Today, source of an average yearly precipitation of 2000 mm of moisture is the Black Sea, situated approximately 40 km to the north of the study site. Glaciers of the Last Glacial Maximum (LGM) descended directly from Mt. Kaçkar and reached an altitude of ca. 1850 m a.s.l. (above sea level). In this study, we are exploring whether the position of Barhal Valley to the south of the main weather divide and its east–west orientation have an influence on the existence and expansion of paleoglaciers. Here, we present 32 new cosmogenic ³⁶Cl dates on erratic boulders from the Çoruh Valley System. We reconstructed three geomorphologically well-contained glacier advances in the Barhal Valley, namely at 34.0 ± 2.3 ka, 22.2 ± 2.6 ka, and 18.3 ± 1.7 ka within the time window of the global LGM. Field evidence shows that the glacier of the 18.3 ± 1.7 ka advance disappeared rapidly and that by the latest time, at 15.6 ± 1.8 ka, the upper cirques were ice-free. No evidence for Lateglacial glacier fluctuations was found, and the Neoglacial activity is restricted to the cirques with rock glaciers. A range of 2700 to 3000 m for the Equilibrium Line Altitude (ELA) at the LGM was reported based on modeling of the glacial morphology. We determined that the most likely position of the LGM ELA in the Çoruh Valley System was at 2900 m a.s.l. We suggest an alternative moisture source to the direct transport from the Black Sea for the ice accumulation in the Eastern Black See Mountains. The shift of the Polar Front and of the Siberian High Pressure System to the south during the LGM resulted in the domination of easterly airflow to the Caucasus and Kaçkar Mountains with moisture from expanded lakes in central–western Siberia and from the enlarged Aral- and Caspian Seas. Full article

This paper focuses on the study of the Vögelsberg landslide located in the municipality of Wattens (Tyrol, Austria), which reactivated in 2016, causing damages to nearby buildings and infrastructures. Since the date of reactivation, a modern monitoring system has been implemented with the installation of in-situ geodetic automated tracking total stations (ATTS), an inclinometer and two piezometers. Here, we describe two distinctive methods, the Breaks for Additive Seasonal and Trend (BFAST) and the Vector Inclination Method (VIM) used to characterize the landslide from the kinematic and geometrical point of view. The main input data, used for both methods, derive from processing a stack of several Sentinel-1 differential interferograms with the Multiple Small Baseline Subset (MSBAS) 2D and 3D algorithms. BFAST allowed highlighting the seasonality of the phenomenon from the analysis of the time series as well as the trend and the breakpoints that identify the landslide reactivation phases. These latter were then correlated with the main triggering factors such as rain and snow melting. The application of the VIM through the exploitation of the MSBAS displacement vectors allowed the reconstruction of the depth of the landslide slip surface along both the longitudinal and transversal direction and, in turn, the evaluation of the volumes of material mobilized by the landslide. The results obtained further prove that procedures for the in-depth analysis of Multi-Temporal Interferometric Synthetic Aperture Radar (MT-InSAR) data can contribute to slow-moving landslide characterization, which represents a fundamental step for landslide hazard assessment within quantitative risk analyses. Full article

To estimate the hazard posed by rock slopes, it is essential to determine the overall stability and potential detachment volume. This is mostly solved using numerical methods together with the strength reduction method (SRM). Many calculation programs do not provide a direct implementation of the Hoek–Brown (HB) criterion. Equivalent Mohr–Coulomb (MC) parameters are often used. Especially for steep rock slopes, the use of equivalent MC parameters with numerical codes and the SRM lead to poor estimates of safety factors. The problem lies in the required and often difficult estimation of a suitable range of minor principal stresses over a ‘slope height’. In the example of the stability analysis of the rock slope Vals in Tyrol/Austria, we show the differences between the application of equivalent MC parameters and a direct application of the HB criterion with apparent MC parameters. The detachment volume and stability are overestimated when applying equivalent MC parameters, as confirmed by calculations with the continuum mechanics code FLAC3D (Itasca Consulting Group). However, the SRM with HB material (i.e., apparent MC parameters) results in a safety factor that cannot be applied to HB parameters. To date, it has not been possible to determine the HB parameters for limit equilibrium via the SRM. This challenge was overcome by fitting an HB envelope to the original HB shear envelope reduced by the safety factor. The envelope is adjusted by two HB variables: GSI and D. This allows to determine the HB parameters at limit equilibrium. It helps to make more realistic predictions about the detachment mechanism and volume. Full article