Geosciences, Volume 10, Issue 1 (January 2020) – 38 articles

Major explosive volcanic eruptions may significantly alter the global atmosphere for about 2–3 years. During that period, volcanic products (mainly H₂SO₄) with high residence time, stored in the stratosphere or, for shorter times, in the troposphere are gradually deposited onto polar ice caps. Antarctic snow may thus record acidic signals providing a history of past volcanic events. The high resolution sulphate concentration profile along a 197 m long ice core drilled at GV7 (Northern Victoria land) was obtained by Ion Chromatography on around 3500 discrete samples. The relatively high accumulation rate (241 ± 13 mm we yr ⁻¹) and the 5-cm sampling resolution allowed a preliminary counted age scale. The obtained stratigraphy covers roughly the last millennium and 24 major volcanic eruptions were identified, dated, and tentatively ascribed to a source volcano. The deposition flux of volcanic sulphate was calculated for each signature and the results were compared with data from other Antarctic ice cores at regional and continental scale. Our results show that the regional variability is of the same order of magnitude as the continental one. Full article

This article presents practical aspects of toponymics in the context of complex analysis of different kinds of data (ethnological, historical, geographical, geological, geopolitical, and sociological). The application of toponymy is here related to the city of Szczecin—a historical city and the recent capital of Western Pomerania, Poland—in order to reveal other same-named localities established on four continents. The historical, geological, and geographical backgrounds of different Szczecin locations is described, with an emphasis on natural values and geoheritage. Analysis of different kinds of data shows several toponymical units that have geographical names related to Szczecin (Stettin, Stettyn, or Settin), which are identified and associated mostly with geological or geographical units. The historical development of the locations described is related mostly to migration and exploration by people originating from Western Pomerania (historically part of Prussia and the German Empire), especially during the colonial and post-colonial ages. Full article

The behaviour of bioavailable trace metals and their stable isotopes in the modern oceans is controlled by uptake into phototrophic organisms and adsorption on and incorporation into marine authigenic minerals. Among other bioessential metals, Cd and its stable isotopes have recently been used in carbonate lithologies as novel tracer for changes in the paleo primary productivity and (bio)geochemical cycling. However, many marine sediments that were deposited during geologically highly relevant episodes and which, thus, urgently require study for a better understanding of the paleo environment are rather composed of a mixture of organic matter (OM), and detrital and authigenic minerals. In this study, we present Cd concentrations and their isotopic compositions as well as trace metal concentrations from sequential leachates of OM-rich shales of the Cryogenian basal Datangpo Formation, Yangtze Platform (South China). Our study shows variable distribution of conservative and bioavailable trace metals as well as Cd isotope compositions between sequential leachates of carbonate, OM, sulphide, and silicate phases. We show that the Cd isotope compositions obtained from OM leachates can be used to calculate the ambient Cryogenian surface seawater of the restricted Nanhua Basin by applying mass balance calculations. By contrast, early diagenetic Mn carbonates and sulphides incorporated the residual Cd from dissolved organic matter that was in isotopic equilibrium with deep/pore waters of the Nanhua Basin. Our model suggests that the Cd isotopic composition of surface seawater at that time reached values of modern oxygenated surface oceans. However, the deep water Cd isotope composition was substantially heavier than that of modern fully oxygenated oceans and rather resembles deep waters with abundant sulphide precipitation typical for modern oxygen minimum zones. This argues for incomplete recycling of Cd and other bioavailable metals shortly after the Sturtian glaciation in the redox stratified Cryogenian Nanhua Basin. Our study highlights the importance of sequential leaching procedures when dealing with impure authigenic sediments such as OM-rich carbonates, mudstones, or shales to achieve reliable trace metal concentrations and Cd isotope compositions as proxies for (bio)geochemical metal cycling in past aquatic systems. Full article

Mathematical models of different degrees of complexity, describing the motion of a snow avalanche along a path with given center line and spatially varying width, are formulated and compared. The most complete model integrates the balance equations for mass and momentum over the cross-section and achieves closure through an entrainment function based on shock theory and a modified Voellmy bed friction law where the Coulombic contribution to the bed shear stress is limited by the shear strength of the snow cover. A simplified model results from integrating these balance equations over the (time-dependent) length of the flow and postulating weak similarity of the evolving avalanche shape. On path segments of constant inclination, it can be solved for the flow depth and speed of the front in closed form in terms of the imaginary error function. Finally, the very simplest model assumes constant flow height and length. On an inclined plane, the evolution of flow depth and velocity predicted by the simplified model are close to those from the full model without entrainment and with corresponding parameters, but the simplest model with constant flow depth predicts much higher velocity values. If the friction coefficient is varied in the full model with entrainment, there can be non-monotonous behavior due to the non-linear interplay between entrainment and the limitation on the Coulomb friction. Full article

The complexity of the Krafla volcano and its geothermal system(s) has puzzled geoscientists for decades. New and old geoscientific studies are reviewed in order to shed some light on this complexity. The geological structure and history of the volcano is more complex than hitherto believed. The visible 110 ka caldera hosts, now buried, an 80 ka inner caldera. Both calderas are bisected by an ESE-WNW transverse low-density structure. Resistivity surveys show that geothermal activity has mainly been within the inner caldera but cut through by the ESE-WNW structure. The complexity of the geothermal system in the main drill field can be understood by considering the tectonic history. Isotope composition of the thermal fluids strongly suggests at least three different geothermal systems. Silicic magma encountered in wells K-39 and IDDP-1 indicates a hitherto overlooked heat transport mechanism in evolved volcanos. Basaltic intrusions into subsided hydrothermally altered basalt melt the hydrated parts, producing a buoyant silicic melt which migrates upwards forming sills at shallow crustal levels which are heat sources for the geothermal system above. This can explain the bimodal behavior of evolved volcanos like Krafla and Askja, with occasional silicic, often phreatic, eruptions but purely basaltic in-between. When substantial amounts of silicic intrusions/magma have accumulated, major basalt intrusion(s) may “ignite” them causing a silicic eruption. Full article

In the context of energy transition, new and renovated buildings often include heating and/or air conditioning energy-saving technologies based on sustainable energy sources, such as groundwater heat pumps with aquifer thermal energy storage. A new aquifer thermal energy storage system was designed and is under construction in the city of Liège, Belgium, along the Meuse River. This system will be the very first to operate in Wallonia (southern Belgium) and should serve as a reference for future shallow geothermal developments in the region. The targeted alluvial aquifer reservoir was thoroughly characterized using geophysics, pumping tests, and dye and heat tracer tests. A 3D groundwater flow heterogeneous numerical model coupled to heat transport was then developed, automatically calibrated with the state-of-the-art pilot points method, and used for simulating and assessing the future system efficiency. A transient simulation was run over a 25 year-period. The potential thermal impact on the aquifer, based on thermal needs from the future building, was simulated at its full capacity in continuous mode and quantified. While the results show some thermal feedback within the wells of the aquifer thermal energy storage system and heat loss to the aquifer, the thermal affected zone in the aquifer extends up to 980 m downstream of the building and the system efficiency seems suitable for long-term thermal energy production. Full article

Physical-based soil erosion models are playing an important role in the assessment of soil erosion, transportation, and deposition in the watershed. Most of these models were developed for data-rich areas of the world and they need a measured data to calibrate and validate their results. To apply such physical-based models, the main factor hindering is the lack of measured sediment data. The amount of sediment in the fluvial systems is the result of hydro-geomorphological processes of a watershed and the nature of stream flows. Therefore, this study aims to develop an alternative empirical model that generates the observed sediment data based on the hydro-geomorphology and nature of stream flows for Ziway Lake basin in the rift Valley of Ethiopia. By applying Soil and water Assessment Tool (SWAT), the lake basin was divided in to two sub-basins (Maki and Katar) with 26 of the watersheds within Maki. The SWAT model was calibrated and validated for both stream and sediment flow by using SUFI-2 program and its performance was assessed by using model evaluation statistics. With calibrated sediment flow rates of 26 Maki sub basins, an empirical model was developed by using its respective drainage area, average sub-basins slope, surface runoff, soil erodibility factor, stream flow rate, and average rive slopes. The applicability of the newly developed alternative model was tested by using model evaluation statistics and validated inside of Katar sub-basin. It is recommended to test the developed model in other basins to incorporate with SWAT CUP program to calibrate and validate the sediment yield at data scared area. Full article

The study of the hydrology of arid regions greatly expanded at the end of the 20th century as humans sought to reduce groundwater pollution from landfills, waste dumps and other forms of land disposal. Historically viewed as wastelands where little or no water percolated to the underlying water table, the discovery of large-scale contamination beneath arid disposal sites such as the Hanford nuclear complex in eastern Washington jumpstarted an industry in studying the hydrology of arid vadose zones and their transport behavior. These studies showed that, in spite of hyper aridity in many areas, precipitation often did infiltrate to deep water. The efforts at Yucca Mountain, Nevada to design a high-level nuclear repository stand out as one of the largest of such studies, and one that fundamentally changed our understanding of not only water flow in fractured rocks, but also of the range of our uncertainty of hydrologic processes in arid regions. In this review and commentary, we present some of the initial concepts of flow at Yucca Mountain, and the evolution in research to quantify the concepts. In light of continued stockpiling of high-level waste, and the renewed interest in opening Yucca Mountain for high-level waste, we then focus on the significant surprises and unanswered questions that remained after the end of the characterization and licensing period; questions that continue to demonstrate the challenges of a geologic repository and our uncertainty about critical processes for long-term, safe storage or disposal of some of our most toxic waste products. Full article

The strength of sand is usually characterized by the maximum value of the secant friction angle. The friction angle is a function of deformation mode, density, and stress level and is strongly correlated with dilatancy at failure. Most often, the friction angle is evaluated from results of conventional compression tests, and correlation between the friction angle of sand at triaxial compression and triaxial extension and plane strain conditions is a vital problem of soil mechanics. These correlations can be obtained from laboratory test results. The failure criteria for sand presented in literature also give the possibility of finding correlations between friction angles for different deformation modes. The general stress-dilatancy relationship obtained from the frictional state concept, with some additional assumptions, gives the possibility of finding theoretical relationships between the friction angle of sand at triaxial compression and triaxial extension and plane strain conditions. The theoretically obtained relationships presented in the paper are fully consistent with theoretical and experimental findings of soil mechanics. Full article

The coastal system, with its physical chemical and biological components, is the place where the lithosphere, atmosphere and hydrosphere intersect and interact, and in which human activity has a considerable impact on the balance of the whole ecosystem. The aim of this study is to assess the health of the coastal environment in a natural protected area in Policoro (Italy). Sand samples have been collected and analysed to detect the presence of environmental contaminants and pollutants (heavy metals). A photogrammetric survey was carried out using an unmanned aerial vehicle (UAV), and an updated orthophoto of the area was obtained using the Structure from Motion (SfM) processing technique. Geostatistical techniques have been used to assess the distribution of the contaminants in the study area. Full article

The UNESCO Global Geoparks is a good didactic instrument that allows teaching territorial elements in transversal way. The geopark provides a means to address the teaching-learning process of both cultural values and natural values, including geology. This paper aims to highlight the educational activities that have emerged from the Villuercas UNESCO Global Geopark. For this purpose, a qualitative methodology based on document analysis and fieldwork has been used. Through a system of categories, we have analyzed the existing materials on the geopark’s website, its activity book and the schools that have been part of its educational project. The geopark and educational centers collaboration has 76 educational proposals that are developed in the territory of the geopark. Since 2011, all schools in its area of action participate in its educational project with a total of 2500 students per year. The Villuercas geopark involves all students in its territory and encourages learning of a significant, transversal and collaborative type. In the same way, it has a positive impact on environmental awareness and respect for the environment and cultural resources. Full article

This paper reports the results of a field-based structural investigation of a well-exposed paleo-accretionary prism, which experienced complex deformation in a low-grade metamorphic setting. Field analyses focused on the description of structural fabrics, with the main emphasis upon parameters like the orientation, style and kinematics of foliations, folds and shear zones. We address the research to the south-westernmost part of the Alpine chain, the Ligurian Alps, where, despite their origin as turbidite sequences deposited into the closing Alpine Tethys Ocean, the Helminthoid Flysch Nappes are presently distributed in the outer part of the chain, above the foreland. The new dataset highlights different deformation patterns related to the different spatial distribution of the flysch units. This regional-scale partitioning of strain is hence associated with progressive deformation within a two-stage geodynamic evolution. Correlations among the different orogenic domains allow the proposal of a kinematic model that describes the motion of the Helminthoid Flysch from the inner to the outer part of the orogen, encompassing the shift from subduction- to collision-related Alpine geodynamic phases. Full article

Watershed-scale nonpoint source (NPS) pollution models have become important tools to understand, evaluate, and predict the negative impacts of NPS pollution on water quality. Today, there are many NPS models available for users. However, different types of models possess different form and structure as well as complexity of computation. It is difficult for users to select an appropriate model for a specific application without a clear understanding of the limitations or strengths for each model or tool. This review evaluates 14 more commonly used watershed-scale NPS pollution models to explain how and when the application of these different models are appropriate for a given effort. The models that are assessed have a wide range of capacities that include simple models used as rapid screening tools (e.g., Long-Term Hydrologic Impact Assessment (L-THIA) and Nonpoint Source Pollution and Erosion Comparison Tool (N-SPECT/OpenNSPECT)), medium-complexity models that require detail data input and limited calibration (e.g., Generalized Watershed Loading Function (GWLF), Loading Simulation Program C (LSPC), Source Loading and Management Model (SLAMM), and Watershed Analysis Risk Management Frame (WARMF)), complex models that provide sophisticated simulation for NPS pollution processes with intensive data and rigorous calibration (e.g., Agricultural Nonpoint Source pollution model (AGNPS/AnnAGNPS), Soil and Water Assessment Tool (SWAT), Stormwater Management Model (SWMM), and Hydrologic Simulation Program Fortran (HSPF)), and modeling systems that integrate various sub-models and tools, and contain the highest complexity to solve all phases of hydrologic, hydraulic, and chemical dynamic processes (e.g., Automated Geospatial Watershed Assessment Tool (AGWA), Better Assessment Science Integrating Point and Nonpoint Sources (BASINS) and Watershed Modeling System (WMS)). This assessment includes model intended use, components or capabilities, suitable land-use type, input parameter type, spatial and temporal scale, simulated pollutants, strengths and limitations, and software availability. Understanding the strengths and weaknesses of each watershed-scale NPS model will lead to better model selection for suitability and help to avoid misinterpretation or misapplication in practice. The article further explains the crucial criteria for model selection, including spatial and temporal considerations, calibration and validation, uncertainty analysis, and future research direction of NPS pollution models. The goal of this work is to provide accurate and concise insight for watershed managers and planners to select the best-suited model to reduce the harm of NPS pollution to watershed ecosystems. Full article

Urban growth, extreme climate, and mismanagement are crucial controlling factors that affect flood vulnerability at wadi catchments. Therefore, this study attempts to understand the impacts of these three factors on the flash flood vulnerability in different climatic regions in Egypt. An integrated approach is presented to evaluate the urban growth from 1984 to 2019 by using Google Images and SENTINEL-2 data, and to develop hazard maps by using a rainfall-runoff-inundation model (RRI). Annual rainfall trend analysis was performed to evaluate the temporal variability trend. The hazard maps that were created were classified into three categories (low, medium, and high) and integrated with the urban growth maps to evaluate the impacts on the flood-vulnerable areas. The results show a significant increase in urban growth resulting in an increase of prone areas for flood hazards over time. However, the degree of this hazard is mainly related to growth directions. Mismanagement affects urban growth directions in both planned and unplanned growth, whether by loss of control over unplanned growth or by deficiencies in approved plans. The rainfall analysis showed that there is no explicit relationship to increases or decreases in the flood vulnerable areas. An urban planning approach is recommended for risk reduction management based on a comprehensive study considering such factors. Full article

Among the indirect estimation approaches of soil water content in the upper layer of the soil, the “triangle method” is one of the most common that relies on the simple relationship between the optical and thermal features sensed via Earth Observation. These features are controlled by water content at the surface and within the root zone but also by meteorological forcing including air temperature and humidity, as well as solar radiation. Night- and day-time MODIS composites of land-surface temperature (LST) allowed applying a version of the triangle method that takes into account the temporal admittance of the soil. In this study, it has been applied to a long time-series of pair images to analyze the seasonal influence of the meteorological forcing on a triangle method index (or temperature–vegetation index, TVX), as well as to discuss extra challenges of the diachronic approach including seasonality effects and the variability of environmental forcing. The Imera Meridionale basin (Sicily, Italy) has been chosen to analyze the method over a time-series of 12 years. The analysis reveals that, under these specific environmental and climatic conditions (strong seasonality and rainfall out of phase with vegetation growth), Normalized Difference Vegetation Index (NDVI) and LST pairs move circularly in time within the optical vs. thermal feature space. Concordantly, the boundaries of the triangle move during the seasons. Results showed a strong correlation between TVX and rainfall normalized amplitudes of the power spectra (r² ~0.8) over the range of frequencies of the main harmonics. Full article

Since early 2018 the “Natural Hazards” Section of Geosciences journal has aimed to publish pure, experimental, or applied research that is focused on advancing methodologies, technologies, expertise, and capabilities to detect, characterize, monitor, and model natural hazards and assess their associated risks. This stream of geoscientific research has reached a high degree of specialization and represents a multi-disciplinary research realm. To inaugurate this section, the Special Issue “Key Topics and Future Perspectives in Natural Hazards Research” was launched. After a year and half since the call for papers was initially opened, the special issue is now completed with the editorial introducing the collection of 10 selected papers covering the following hot topics of natural hazards research: (i) trends in publications and research directions at international level; (ii) the role of Big Data in natural disaster management; (iii) assessment of seismic risk through the understanding and quantification of its three components (i.e., hazard, vulnerability and exposure/impact); (iv) climatic/hydro-meteorological hazards (i.e., drought, hurricanes); and (v) scientific analysis of past incidents and disaster forensics (i.e., the Oroville Dam 2017 spillway incident). The present editorial provides a summary of each paper of the collection within the current context of scientific research on natural hazards, pointing out the salient results and key messages. Full article

The present research aims to estimate the influence of site amplification on the seismic vulnerability of the historical centre of the municipality of Baranello in the Molise Region of Italy. Firstly, a structural and typological characterization of the investigated area has been done according to the EMS-98 scale. Subsequently, the vulnerability assessment of the historical buildings located there has been carried out through an appropriate survey form in order to identify the buildings which are most susceptible to seismic damage. To this purpose, the seismic event occurring in October of 2002 has been selected as a reference earthquake. Moreover, according to the AeDES form implemented by the Italian Civil Protection Department to evaluate the usability of constructions after seismic events, the calibration of the typological vulnerability curves of the built-up area has been done, and a quantitative assessment of the local seismic response has been achieved, based on the seismic motions recorded after the 2002 Molise earthquake. Finally, the local amplification factor, which negatively influences the severity of the seismic damage on the structures, has been taken into account in order to more correctly foresee the expected damage of the inspected urban sector, so as to use more appropriately the achieved results for reliable seismic risk mitigation plans. Full article

The present paper deals with the practical problem of reducing statistical uncertainty in elastic settlement analysis of shallow foundations by relying on targeted field investigation with the aim of an optimal design. In a targeted field investigation, the optimal number and location of sampling points are known a priori. As samples are taken from the material field (i.e., the ground), which simultaneously is a stress field (stresses caused by the footing), the coexistence of these two fields allows for some points in the ground to better characterize the serviceability state of structure. These points are identified herein through an extensive parametric analysis of the factors controlling the magnitude of settlement; the number of different cases considered was 3318. This is done in an advanced probabilistic framework using the Random Finite Element Method (RFEM) properly considering sampling of soil property values. In this respect, the open source RSETL2D program, which combines elastic finite element analysis with the theory of random fields, has been modified as to include the function of sampling of soil property values from the generated random fields and return the failure probability of footing against excessive settlement. Two sampling strategies are examined: (a) sampling from a single point and (b) sampling a domain (the latter refers to e.g., continuous cone penetration test data). As is shown in this work, by adopting the proper sampling strategy (defined by the number and location of sampling points), the statistical error can be significantly reduced. The error is quantified by the difference in the probability of failure comparing different sampling scenarios. Finally, from the present analysis, it is inferred that the benefit from a targeted field investigation is much greater as compared to the benefit from the use of characteristic values in a limit state design framework. Full article

Engineering methods such as soil nails, geosynthetic reinforcement, retaining structures, gabions, and shotcrete are implemented to stabilize road cut slopes along mountainous areas. However, these structures are not environmentally friendly and, particularly in Ethiopia, it is impossible to address all road problems due to financial limitations. Nowadays, soil reinforcement with plant roots is recognized as an environmentally sustainable alternative to improve shallow slope failure along mountainous transportation corridors. The aims of this study was, therefore, to conduct slope stability analysis along a road corridor by incorporating the effect of plant roots. Five plant species were selected for the analysis based on their mechanical characteristics. Namely, Eucalyptus globules (tree), Psidium guajava (shrub), Salix subserrata (shrub), Chrysopogon zizanioides, and Pennisetum macrourum (grasses). The roots’ tensile strength and soil parameters were determined through tensile strength testing and triaxial compression tests, respectively. The factor of safety of the slope was calculated by the PLAXIS-2D software. The study showed that when the slope was reinforced with plant roots, the factor of safety (FOS) improved from 22–34%. The decreasing effect of vegetation on slope stability was observed when soil moisture increased. The sensitivity analysis also indicated that: (1) as the spacing between plants decreased, the effect of vegetation on the slope increased. (2) Slope angle modification with a combination of plant roots had a significant impact on slope stabilization. Of the five-selected plant species, Salix subserrata was the promising plant species for slope stabilization as it exhibited better root mechanical properties among selected plant species. Full article

Structural pounding between adjacent, insufficiently separated buildings, or bridge segments, has been repeatedly observed during seismic excitations. Such earthquake-induced collisions may cause severe structural damage or even lead to the collapse of colliding structures. The aim of the present paper was to show the results of the study focused on determination of peak impact forces during collisions between buildings exposed to different seismic excitations. A set of different ground motion records, with various peak ground acceleration (PGA) values and frequency contents, were considered. First, pounding-involved numerical analysis was conducted for the basic parameters of colliding buildings. Then, the parametric study was carried out for different structural natural periods, structural damping ratios, gap sizes between buildings and coefficients of restitution. The results of the analysis conducted for the basic structural parameters indicate that the largest response of the analysed buildings was observed for the Duzce earthquake. The parametric study showed that the pounding-involved structural response depended substantially on all parameters considered in the analysis, and the largest response was observed for different ground motions. The results of the study presented in this paper indicate that the value of the peak impact force expected during the time of the earthquake does not depend on the PGA value of ground motion, but rather on the frequency contents of excitation and pounding scenario. It is therefore recommended that the peak impact force for buildings exposed to structural pounding during earthquakes should be determined individually for the specific structural configuration taking into account the design ground motion. Full article

The Matese and Beneventano areas coincide with the transition from the central to the southern Apennines and are characterized by both SW- and NE-dipping normal faulting seismogenic structures, responsible for the large historical earthquakes. We studied the Matese and Beneventano seismicity by means of high-precision locations of earthquakes spanning from 29 December 2013 to 4 September 2018. Events were located by using all of the available data from temporary and permanent stations in the area and a 1D computed velocity model, inverting the dataset with the Velest code. For events M > 2.8 we used P- and S-waves arrival times of the strong motion stations located in the study area. A constant value of 1.83 for V_p/V_s was computed with a modified Wadati method. The dataset consists of 2378 earthquakes, 18,715 P- and 12,295 S-wave arrival times. We computed 55 new fault plane solutions. The mechanisms show predominantly normal fault movements, with T-axis trends oriented NE–SW. Only relatively small E–W trending clusters in the eastern peripheral zones of the Apenninic belt show right-lateral strike-slip kinematics similar to that observed in the Potenza (1990–1991) and Molise (2002 and 2018) sequences. These belong to transfer zones associated with differential slab retreat of the Adriatic plate subduction beneath the Apennines. The Matese sequence (December 2013–February 2014; main shock M_w 5.0) is the most relevant part of our dataset. Hypocentral depths along the axis of the Apenninic belt are in agreement with previous seismological studies that place most of the earthquakes in the brittle upper crust. We confirm a general deepening of seismicity moving from west to the east along the Apennines. Seismicity depth is controlled by heat-flow, which is lower in the eastern side, thus causing a deeper brittle–ductile transition. Full article

All 142 known species of Curculionoidea in Eocene amber are documented, including one species of Nemonychidae, 16 species of Anthribidae, six species of Belidae, 10 species of Rhynchitidae, 13 species of Brentidae, 70 species of Curcuionidae, two species of Platypodidae, and 24 species of Scolytidae. Oise amber has eight species, Baltic amber has 118 species, and Rovno amber has 16 species. Nine new genera and 18 new species are described from Baltic amber. Four new synonyms are noted: Palaeometrioxena Legalov, 2012, syn. nov. is synonymous with Archimetrioxena Voss, 1953; Paleopissodes weigangae Ulke, 1947, syn. nov. is synonymous with Electrotribus theryi Hustache, 1942; Electrotribus erectosquamata Rheinheimer, 2007, syn. nov. is synonymous with Succinostyphlus mroczkowskii Kuska, 1996; Protonaupactus Zherikhin, 1971, syn. nov. is synonymous with Paonaupactus Voss, 1953. Keys for Eocene amber Curculionoidea are given. There are the first records of Aedemonini and Camarotini, and genera Limalophus and Cenocephalus in Baltic amber. Full article

Here we describe the molecular composition and resistance to decomposition of humic acids isolated from selected soils of the Russian Arctic and Antarctic. The degree of soil organic matter stabilization was assessed using modern instrumental methods: nuclear magnetic resonance spectroscopy (cross peak magic-angle spinning (CP/MAS) ¹³C-NMR and ¹H-¹³C heteronuclear-correlation (HETCOR)). Analysis of the humic acids showed that aromatic compounds prevail in the organic matter formed in cryoconites, located on the surfaces of the glaciers. The predominance of aliphatic fragments is revealed in the soils of the Yamal peninsula and Antarctica. This could be caused by sedimentation of fresh organic matter exhibiting low decomposition stage due to the severe climate and processes of hydrogenation in the humic acids, destruction of the C-C bonds, and formation of chains with high hydrogen content. These processes result in formation of aliphatic fragments in the humic acids. Full article

A series of six gravity cores has been used to reconstruct the depositional history of Hatton Bank (Rockall Plateau, NE Atlantic Ocean). The cores have been studied for magnetic susceptibility (MS), geochemical composition, grain size distribution, and a semi-quantitative foraminiferal association. Two main interbedded facies have been described: (i) calcareous ooze; and (ii) lithogenous silt. The study reveals prominent peaks from the MS signal, silt, Mg/Ca, Fe/Ca, Al/Ca, and Rare Earth Elements normalised by Continental Crust (REE/CC), which are sensitive indicators for Heinrich events (H1, H2, H3, H4, and H5) and ash layers. These peaks may relate to alternations in dominance of the calcareous and lithogenic facies. The sediment displays a high percentage of carbonate in interglacial layers but is lithogenic-dominated in glacial stages. The layers with prominent lithic-rich and foraminifera-poor sediments (established as Heinrich layers) may be related to a possible palaeoclimatic effect, where freshwater discharged during iceberg melting may have reduced the formation of North Atlantic Deep Water (NADW). In the study area, the mean sedimentation rates for the last glacial as ~4.2 cm ka⁻¹ and ~1.4 cm ka⁻¹ for the last ~18 ka (interglacial period) have been estimated. Besides this evidence, Fe/Ca and MS peaks may reflect the presence of basalt, either introduced through ice-rafting or transported and redistributed by bottom currents in the study area. Certain indices, including MS and Fe/Ca, are proposed as good proxies for detecting Heinrich events and ash layers in the Hatton Bank sediments and, in consequence, are parameters that can be used to infer strengthened/weakened NADW formation, according to stadials/interstadials. Moreover, we suggest that the northernmost boundary of the area with evidence of Heinrich events may be situated around 57°38′ N in the Hatton–Rockall area, at least for H4, based on the variation of the Mg/Ca and Fe/Ca curves. Full article

Natural and human-induced impacts on water resources across the globe continue to negatively impact water resources. Characterizing the hydrologic sensitivity to climatic and anthropogenic changes is problematic given the lack of monitoring networks and global-scale model uncertainties. This study presents an integrated methodology combining satellite remote sensing (e.g., GRACE, TRMM), hydrologic modeling (e.g., SWAT), and climate projections (IPCC AR5), to evaluate the impact of climatic and man-made changes on groundwater and surface water resources. The approach was carried out on two scales: regional (Morocco) and watershed (Souss Basin, Morocco) to capture the recent climatic changes in precipitation and total water storage, examine current and projected impacts on total water resources (surface and groundwater), and investigate the link between climate change and groundwater resources. Simulated (1979–2014) potential renewable groundwater resources obtained from SWAT are ~4.3 × 10⁸ m³/yr. GRACE data (2002–2016) indicates a decline in total water storage anomaly of ~0.019m/yr., while precipitation remains relatively constant through the same time period (2002–2016), suggesting human interactions as the major underlying cause of depleting groundwater reserves. Results highlight the need for further conservation of diminishing groundwater resources and a more complete understanding of the links and impacts of climate change on groundwater resources. Full article

Ice cliffs within a glacier represent a challenge for the continuity equations used in many glacier models by interrupting the validity of input parameters. In the case of Gourdon Glacier on James Ross Island, Antarctica, a ∼300–500 $\mathrm{m}$ high, almost vertical cliff, separates the outlet glacier from its main accumulation area on the plateau of the island. In 2017 and 2018 we conducted ice thickness measurements during two airborne ground penetrating radar campaigns in order to evaluate differences to older measurements from the 1990s. The observed differences are mostly smaller than the estimated error bars. In comparison to the in situ data, the published “consensus ice thickness estimate” strongly overestimates the ice thickness at the outlet. We analyse three different interpolation and ice thickness reconstruction methods. One approach additionally includes the mass input from the plateau. Differences between the interpolation methods have a minor impact on the ice discharge estimation if the used flux gates are in areas with a good coverage of in situ measurements. A much stronger influence was observed by uncertainties in the glacier velocities derived from remote sensing, especially in the direction of the velocity vector in proximity to the ice cliff. We conclude that the amount of in situ measurements should be increased for specific glacier types in order to detect biases in modeled ice thickness and ice discharge estimations. Full article

The Kola region hosts numerous Paleozoic massifs of ultrabasic alkaline rocks and carbonatites with deposits of commercially valuable metals, such as iron, tantalum, niobium, and rare earth elements. These magmatic complexes are characterized by high contents of alkaline elements at generally low contents of SiO₂ and/or Al₂O₃. In this study, we examined the precursors to the formation of the unique Paleozoic alkaline province through studying the early Precambrian stages in the evolution of the Kola collision area, from where these unique features probably originated. We mathematically modeled the changes in the chemical composition of these rocks. The obtained data can be used for metallogenic forecasting, which indicated a number of Precambrian objects in the region, namely, the Lapland Granulite Belt of the Kola region and granulite belts in Eurasia. The mathematical modeling performed during this research depicted a linear trend that defined the style of the changes in the chemical composition at the transition from the metaultrabasic-basic rocks of the Lapland granulite belt to the group of belts in Eurasia. These differences are statistically significant with respect to the obtained trend (chemical composition projected on the trend), mainly manifested as increased SiO₂ and Al₂O₃ contents with a decreasing total alkalis content, which is opposite to the indicated trends of the changing chemical composition in the Paleozoic alkaline rock units of the Kola region. We concluded that one of the reasons for the unique composition of the Paleozoic magmatism products could be a specific feature of the earlier Neoarchean stages of the tectonic-magmatic activity in the northeastern Baltic Shield, which implies a close relationship between later geological events and the early Precambrian history, at least in the study area. Full article

Depth-integrated simulations of snow avalanches have become a central part of risk analysis and mitigation. However, the common practice of applying different model parameters to mimic different avalanches is unsatisfying. In here, we analyse this issue in terms of two differently sized avalanches from the full-scale avalanche test-site Vallée de la Sionne, Switzerland. We perform depth-integrated simulations with the toolkit OpenFOAM, simulating both events with the same set of model parameters. Simulation results are validated with high-resolution position data from the GEODAR radar. Rather than conducting extensive post-processing to match radar data to the output of the simulations, we generate synthetic flow signatures inside the flow model. The synthetic radar data can be directly compared with the GEODAR measurements. The comparison reveals weaknesses of the model, generally at the tail and specifically by overestimating the runout of the smaller event. Both issues are addressed by explicitly considering deposition processes in the depth-integrated model. The new deposition model significantly improves the simulation of the small avalanche, making it starve in the steep middle part of the slope. Furthermore, the deposition model enables more accurate simulations of deposition patterns and volumes and the simulation of avalanche series that are influenced by previous deposits. Full article

Anomalies of the electric field potential gradient have been observed in the near-ground air before earthquakes in different regions of the world. Such anomalies are likely caused by radon air ionization. In this study, the impact of this precursor was estimated according to continuous observations of the electric field in Kamchatka in 1997–2002. Full article