Editor’s Choice Articles

Landslides are one of the most destructive and most recurring natural calamities in the Himalayan region. Their occurrence leads to immense damage to infrastructure and loss of land, human lives, and livestock. One of the most affected regions is the Bhutan Himalayas, where the majority of the landslides are rainfall-induced. The present study aims to determine the hazard and risk associated with rainfall-induced landslides for the Phuentsholing region located in the southwestern part of the Bhutan Himalayas. The work involves developing a landslide risk map using hazard and vulnerability maps utilizing landslide records from 2004 to 2014. The landslide hazard map was generated by determining spatial and temporal probabilities for the study region. The spatial probability was computed by analyzing the landslide contributing factors like geology, slope, elevation, rainfall, and vegetation based on comprehensive field study and expertise about the area. The contributing factors were divided into various classes and the percentage of landslide occurrence under each class was calculated to understand its contributing significance. Thereafter, a weighted linear combination approach was used in a GIS environment to develop the spatial probability map which was multiplied with temporal probabilities based on regional rainfall thresholds already determined for the region. Consequently, vulnerability assessment was conducted using key elements at risk (population, land use/land cover, proximity to road, proximity to stream) and the weights were provided based on expert judgment and comprehensive field study. Finally, risk was determined and the various regions in the study area were categorized as high, medium, and low risk. Such a study is necessary for low-economic countries like Bhutan which suffers from unavailability of extensive data and research. The study is conducted for a specific region but can be extended to other areas around the investigated area. The tool can serve as an indicator for the civil authorities to analyze the risk posed by landslides due to the rapid infrastructure development in the region. Full article

Rockfalls evolve rapidly and unpredictably in mountain environments and can cause considerable losses to human societies, structures, economical activities, and also natural and historical heritage. Rockfall risk analyses are complex and multi-scale processes involving several disciplines and techniques. This complexity is due to the main features of rockfall phenomena, which are extremely variable over space and time. Today, a considerable number of methods exists for protecting land, as well as assessing and managing the risk level. These methodologies are often very different from each other, depending on the data required, the purposes of the analysis, and the reference scale adopted, i.e., the analysis level of detail. Nevertheless, several questions still remain open with reference to each phase of the hazard and risk process. This paper is devoted to a general overview of existing risk estimation methodologies and a critical analysis of some open questions with the aim of highlighting possible further research topics. A typical risk assessment framework is exemplified by analyzing a real case study. Each step of the process is treated at both the detailed and the large scale in order to highlight the main characteristics of each level of detail. Full article

Deep slope movements and, eventually, slope failure, have been often interpreted to be due to slope-vegetation-atmosphere interaction on slopes formed of clayey materials in the Italian Southern-Eastern Apennines, as reported in the literature. Such slopes are generally formed of flysch, within which clay is the main lithotype. Such clays are characterized by a disturbed meso-fabric, as an effect of the intense tectonics. The paper presents the results of coupled hydromechanical numerical analyses of the slope-vegetation-atmosphere interaction for a clay slope representative for the geomechanical scenario where such climate-induced deep slope movements have been repeatedly recorded. In the analyses, different model initialization procedures and parameter values were adopted. The comparison of the numerical results with the site data is aimed at assessing the effects of the soil-vegetation-atmosphere interaction taking place in the top strata of the slope, on the stress-strain conditions across the whole slope, and on the slope stability. The comparison between the numerical results of the analyses carried out entailing different initialization stages are intended to evaluate the influence of such a stage on the model predictions. It is found that only when the slope model initialization accounts for the slope loading history, developed over geological time, the numerical predictions get close to the site observations. In such case, the numerical results confirm that deep movements consequent to progressive failure may take place in clay slopes due to the slope-vegetation-atmosphere interaction. Full article

The Vigan-Aggao Fault is a 140-km-long complex active fault system consisting of multiple traces in the westernmost part of the Philippine Fault Zone (PFZ) in northern Luzon, the Philippines. In this paper, its traces, segmentation, and oblique left-lateral strike-slip motion are determined from horizontal and vertical displacements measured from over a thousand piercing points pricked from displaced spurs and streams observed from Google Earth Pro satellite images. This work marks the first instance of the extensive use of Google Earth as a tool in mapping and determining the kinematics of active faults. Complete 3D image coverage of a major thoroughgoing active fault system is freely and easily accessible on the Google Earth Pro platform. It provides a great advantage to researchers collecting morphotectonic displacement data, especially where access to aerial photos covering the entire fault system is next to impossible. This tool has not been applied in the past due to apprehensions on the positional measurement accuracy (mainly of the vertical component). The new method outlined in this paper demonstrates the applicability of this tool in the detailed mapping of active fault traces through a neotectonic analysis of fault-zone features. From the sense of motion of the active faults in northern Luzon and of the major bounding faults in central Luzon, the nature of deformation in these regions can be inferred. An understanding of the kinematics is critical in appreciating the distribution and the preferred mode of accommodation of deformation by faulting in central and northern Luzon resulting from oblique convergence of the Sunda Plate and the Philippine Sea Plate. The location, extent, segmentation patterns, and sense of motion of active faults are critical in coming up with reasonable estimates of the hazards involved and identifying areas prone to these hazards. The magnitude of earthquakes is also partly dependent on the type and nature of fault movement. With a proper evaluation of these parameters, earthquake hazards and their effects in different tectonic settings worldwide can be estimated more accurately. Full article

The distinction between turbidites, contourites and hemipelagites in modern and ancient deep-water systems has long been a matter of controversy. This is partly because the processes themselves show a degree of overlap as part of a continuum, so that the deposit characteristics also overlap. In addition, the three facies types commonly occur within interbedded sequences of continental margin deposits. The nature of these end-member processes and their physical parameters are becoming much better known and are summarised here briefly. Good progress has also been made over the past decade in recognising differences between end-member facies in terms of their sedimentary structures, facies sequences, ichnofacies, sediment textures, composition and microfabric. These characteristics are summarised here in terms of standard facies models and the variations from these models that are typically encountered in natural systems. Nevertheless, it must be acknowledged that clear distinction is not always possible on the basis of sedimentary characteristics alone, and that uncertainties should be highlighted in any interpretation. A three-scale approach to distinction for all deep-water facies types should be attempted wherever possible, including large-scale (oceanographic and tectonic setting), regional-scale (architecture and association) and small-scale (sediment facies) observations. Full article

Hydrovolcanism is a type of volcanism where magma and water interact either explosively or non-explosively. The less frequently used term, hydromagmatism, includes all the processes responsible for magma and water interaction in a magmatic system. Hydrovolcanism is commonly used as a synonym for phreatomagmatism. However, in recent years phreatomagmatism appears more in association with volcanic eruptions that occur in shallow subaqueous or terrestrial settings and commonly involves molten fuel-coolant interaction (MFCI) driven processes. Here a revised and reviewed classification scheme is suggested on the basis of the geo-environment in which the magma-water interaction takes place and the explosivity plus mode of energy transfer required to generate kinetic energy to produce pyroclasts. Over the past decade researchers have focused on the role hydrovolcanism/phreatomagmatism plays in the formation of maar craters, the evolution of diatremes and the signatures of magma—water interaction in the geological record. In the past five years, lithofacies-characterization is the most common approach to studying hydrovolcanism. By far mafic monogenetic volcanic fields generated the greatest number of research results. Significant knowledge gaps are identified, especially in developing tools to identify the textural signatures hydrovolcanism leave behind on eruptive products and exploring the role of hydrovolcanism in the growth of intermediate and silicic small volume volcanoes. 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

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

2I/Borisov is the second ever interstellar object (ISO). It is very different from the first ISO ’Oumuamua by showing cometary activities, and hence provides a unique opportunity to study comets that are formed around other stars. Here we present early imaging and spectroscopic follow-ups to study its properties, which reveal an (up to) 5.9 km comet with an extended coma and a short tail. Our spectroscopic data do not reveal any emission lines between 4000–9000 Angstrom; nevertheless, we are able to put an upper limit on the flux of the C2 emission line, suggesting modest cometary activities at early epochs. These properties are similar to comets in the solar system, and suggest that 2I/Borisov—while from another star—is not too different from its solar siblings. Full article

In the framework of this work, studies on the Yamal crater formed as a result of a cryogenic eruption of a water-gas fluid were carried out. The structure and variations of the composition of the geochemical field along the section of the upper horizons of permafrost are considered on the basis of field work, including the drilling of boreholes near the crater. The obtained regularities of the distribution of chemical elements, and gases between the mineral component of the soil and meltwater, suggest that permafrost at the site of the funnel are the remains of a sub-lake paleo-talik, from which mineralized water and gases were expulsed into the yet unfrozen reservoir that previously existed at this place. The component composition of gases suggests that they are products of biochemical processes similar to those that occur in modern peatlands. The δ¹³C value for methane extracted from the sediments of the near-contact zone of the Yamal crater was found to be −76‰. The predominance of high molecular weight normal alkanes in frozen bitumen indicates that the original organic substrate which was buried contained remains of higher vegetation. The Yamal funnel was formed by the sediment’s “explosion” while the water-gas fluid was released. The volume of the ejected sediment amounted to about 220 thousand m³. Full article

The fossil record of Tenebrionidae (excluding the Quartenary) is presented. In total, 122 fossil species, clearly belonging to the family, are known; some beetles were determined only to genus; 78 genera are listed in the fossil record, including 29 extinct genera. The great diversity of tenebrionids occurs in the Lower Cretaceous Lagerstätte of China (Yixian Formation), Middle Paleocene of France (Menat), Lower Eocene deposits of Germany (Geiseltal), Upper Eocene Baltic amber (Eastern Europe), Upper Eocene deposits of Florissant Formation (USA) and Miocene (Dominican amber). Tenebrionids of the following major lineages, including seven subfamilies, are currently known in the fossil record. These include the lagrioid branch (Lagriinae, Nilioninae), pimelioid branch (Pimeliinae), and tenebrioid branch (Alleculinae, Tenebrioninae, Diaperinae, Stenochiinae). The importance of the fossil record for evolutionary reconstructions and phylogenetic patterns is discussed. The oldest Jurassic and Early Cretaceous darkling beetles of the tenebrionoid branch consist of humid-adapted groups from the extant tribes Alleculini, Ctenopodiini (Alleculinae), and Alphitobiini (Tenebrioninae). Thus, paleontological evidence suggests that differentiation of the family started at least by the Middle Jurassic but does not indicate that xerophilic darkling beetles differentiated much earlier than mesophilic groups. Full article

The complex cyclic shear stress path experienced by the soil during an earthquake, which could also induce liquefaction phenomena, can be approximated in the laboratory only by using sophisticated testing apparatuses. Cyclic triaxial tests have been widely used, especially for coarse grained soils, as in this study. In the framework of the design for the seismic retrofitting of the ‘‘Ritiro viaduct’’ foundations along the A20 motorway connecting Messina with Palermo (Italy), a soil liquefaction study was also carried out. With this aim, a detailed geological and geotechnical characterization of the area was performed by in situ and laboratory tests, including seismic dilatometer Marchetti tests (SDMTs), the combined resonant column (RCT) and cyclic loading torsional shear tests (CLTSTs), and undrained cyclic loading triaxial tests (CLTxTs). In particular, the paper presents the results of cyclic triaxial tests carried out on isotropically consolidated specimens of a sandy soil. The seismic retrofitting works include the reinforcement of the foundation and replacement of the decks with newly designed type and structural schemes, mixed steel, and concrete with continuous girder. During the investigation, data were acquired for the characterization of materials, for the definition of degradation phenomena with the relative identification of possible causes, and for the estimation of the residual performance characteristics of the building. The structural campaign of investigations necessary to determine all of the key parameters useful for a correct definition of the residual performance capabilities of the work was divided into two phases: One in situ and one in the laboratory. Full article

Forecasting estuarine circulation is a hot topic, especially in densely populated regions, like Santos (Brazil). This paper aims to improve a water-level forecasting system for the Santos estuary, particularly the physical forcing determining the residual tide, which in extreme cases increase the predicting errors. The MOHID hydrodynamic model was implemented with a nested downscaling approach. All automatic procedures to provide a high-resolution real-time forecast system are managed by the AQUASAFE software. Water-level observation and prediction datasets (2016–2017) of five tide gauges in the Santos channel were analyzed, resulting in distinct model configurations, aiming to minimize forecasting inaccuracies. Current MOHID open boundary reference solutions were modified: the astronomical solution was updated from FES2012 to FES2014 whereas the meteorological component (Copernicus Marine Environment Monitoring Service (CMEMS) global solution) time resolution was altered from daily to hourly data. Furthermore, the correlation between significant wave height with positive residual tide events was identified. The model validation presented a minimum Root Mean Square Error (RMSE) of 12.5 cm. Despite FES2014 solution improvements at the bay entrance, errors increase in inner stations were maintained, portraying the need for better bathymetric data. The use of a CMEMS hourly resolution decreased the meteorological tide errors. A linear regression method was developed to correct the residual tide through post-processing, under specific wave height conditions. Overall, the newest implementation increased the water-level forecast accuracy, particularly under extreme events. Full article

Dramatic changes in sea ice have been observed in both poles in recent decades. However, the observational period for sea ice is short, and the climate models tasked with predicting future change in sea ice struggle to capture the current Antarctic trends. Paleoclimate archives, from marine sedimentary records and coastal Antarctic ice cores, provide a means of understanding sea ice variability and its drivers over decadal to centennial timescales. In this study, we collate published records of Antarctic sea ice over the past 2000 years (2 ka). We evaluate the current proxies and explore the potential of combining marine and ice core records to produce multi-archive reconstructions. Despite identifying 92 sea ice reconstructions, the spatial and temporal resolution is only sufficient to reconstruct circum-Antarctic sea ice during the 20th century, not the full 2 ka. Our synthesis reveals a 90 year trend of increasing sea ice in the Ross Sea and declining sea ice in the Bellingshausen, comparable with observed trends since 1979. Reconstructions in the Weddell Sea, the Western Pacific and the Indian Ocean reveal small negative trends in sea ice during the 20th century (1900–1990), in contrast to the observed sea ice expansion in these regions since 1979. Full article

Some (ultra)high-pressure metamorphic rocks that formed during continental collision preserve relict minerals, indicating a two-stage evolution: first, subduction to mantle depths and exhumation to the lower-crustal level (with simultaneous cooling), followed by intensive heating that can be characterized by a β-shaped pressure–temperature–time (P–T–t) path. Based on a two-dimensional (2D) coupled petrological–thermomechanical tectono-magmatic numerical model, we propose a possible sequence of tectonic stages that could lead to these overprinting metamorphic events along an orogenic β-shaped P–T–t path: the subduction and exhumation of continental crust, followed by slab retreat that leads to extension and subsequent asthenospheric upwelling. During the last stage, the exhumed crustal material at the crust–mantle boundary undergoes heating from the underlying hot asthenospheric mantle. This slab rollback scenario is further compared numerically with the classical continental collision scenario associated with slab breakoff, which is often used to explain the late heating impulse in the collisional orogens. The mantle upwelling occurring in the experiments with slab breakoff, which is responsible for the heating of the exhumed crustal material, is not related to the slab breakoff but can be caused either by slab bending before slab breakoff or by post-breakoff exhumation of the subducted crust. Our numerical modeling predictions align well with a variety of orogenic P–T–t paths that have been reported from many Phanerozoic collisional orogens, such as the Variscan Bohemian Massif, the Triassic Dabie Shan, the Cenozoic Northwest Himalaya, and some metamorphic complexes in the Alps. Full article

This work is focused on the landslide susceptibility assessment, applied to Mauritius Island. The study area is a volcanic island located in the western part of the Indian Ocean and it is characterized by a plateau-like morphology interrupted by three rugged mountain areas. The island is severely affected by geo-hydrological hazards, generally triggered by tropical storms and cyclones. The landslide susceptibility analysis was performed through an integrated approach based on morphometric analysis and preliminary Geographical Information System (GIS)-based techniques, supported by photogeological analysis and geomorphological field mapping. The analysis was completed following a mixed heuristic and statistical approach, integrated using GIS technology. This approach led to the identification of eight landslide controlling factors. Hence, each factor was evaluated by assigning appropriate expert-based weights and analyzed for the construction of thematic maps. Finally, all the collected data were mapped through a cartographic overlay process in order to realize a new zonation of landslide susceptibility. The resulting map was grouped into four landslide susceptibility classes: low, medium, high, and very high. This work provides a scientific basis that could be effectively applied in other tropical areas showing similar climatic and geomorphological features, in order to develop sustainable territorial planning, emergency management, and loss-reduction measures. Full article

Earthquake occurrence is usually unpredictable apart from sites in the vicinity of volcanoes. It is not easy to measure displacements caused by seismic phenomena using classical geodetic methods, which are based on point survey. Therefore, the surveying of ground movements caused by seismic events should be carried out continuously. Nowadays, remote sensing data and InSAR are often applied to monitor ground displacements in areas affected by seismicity. The effects of severe nearby mining-induced earthquakes have been discussed in the paper. The earthquakes occurred in 2017 and had a magnitude of 4.7 and 4.8. The distance between the epicenters of the mining-induced earthquakes was around 1.6 km. The aim of the investigation has been to analyze the spatio-temporal distribution of ground movements caused by the two tremors using the InSAR technique. Superposition of surface displacement has been studied in time and space. The main scientific aim has been to prove that in the areas where high-energy tremors occur, ground movements overlap. Due to proximity between the epicenters, the mining-induced earthquakes caused the formation of a large subsidence trough with the dimension of approximately 1.2 km × 4.2 km and total subsidence of ca. 116 mm. Two-time phases of subsidence were determined with temporal overlapping. The subsidence analysis has enhanced the cognition of the impact of mining-induced seismicity on the kinematics of surface changes. Moreover, the present work supports the thesis that InSAR is a valuable and adequately accurate technique to monitor ground displacements caused by mining induced earthquakes. Full article

Over the past two decades, La Soufrière volcano in Guadeloupe has displayed a growing degassing unrest whose actual source mechanism still remains unclear. Based on new measurements of the chemistry and mass flux of fumarolic gas emissions from the volcano, here we reveal spatio-temporal variations in the degassing features that closely relate to the 3D underground circulation of fumarolic fluids, as imaged by electrical resistivity tomography, and to geodetic-seismic signals recorded over the past two decades. Discrete monthly surveys of gas plumes from the various vents on La Soufrière lava dome, performed with portable MultiGAS analyzers, reveal important differences in the chemical proportions and fluxes of H₂O, CO₂, H₂S, SO₂ and H_2, which depend on the vent location with respect to the underground circulation of fluids. In particular, the main central vents, though directly connected to the volcano conduit and preferentially surveyed in past decades, display much higher CO₂/SO₂ and H₂S/SO₂ ratios than peripheral gas emissions, reflecting greater SO₂ scrubbing in the boiling hydrothermal water at 80–100 m depth. Gas fluxes demonstrate an increased bulk degassing of the volcano over the past 10 years, but also a recent spatial shift in fumarolic degassing intensity from the center of the lava dome towards its SE–NE sector and the Breislack fracture. Such a spatial shift is in agreement with both extensometric and seismic evidence of fault widening in this sector due to slow gravitational sliding of the southern dome sector. Our study thus provides an improved framework to monitor and interpret the evolution of gas emissions from La Soufrière in the future and to better forecast hazards from this dangerous andesitic volcano. Full article

The relative impact of water and wind on total erosion was investigated by means of an experimental-empirical study. Wind erosion and water erosion were measured at five different sites: (1) Mediterranean fallow, (2) Mediterranean orchard, (3) wheat field, (4) vineyard and (5) sand substrate. Mean erosion rates ranged from 1.55 to 618 g·m⁻²·h⁻¹ for wind and from 0.09 to 133.90 g·m⁻²·h⁻¹ for rain eroded material over all tested sites. Percentages (%) of eroded sediment for wind and rain, respectively, were found to be 2:98 on Mediterranean fallow, 11:89 on Mediterranean orchard, 3:97 on wheat field, 98:2 on vineyard and 99:1 on sand substrate. For the special case of soil surface crust destroyed by goat trampling, the measured values emphasize a strong potential impact of herding on total soil erosion. All sites produced erosion by wind and rain, and relations show that both erosive forces may have an impact on total soil erosion depending on site characteristics. The results indicate a strong need to focus on both wind and water erosion particularly concerning soils and substrates in vulnerable environments. Measured rates show a general potential erosion depending on recent developments of land use and climate change and may raise awareness of scientist, farmers and decision makers about potential impact of both erosive forces. Knowledge about exact relationship is key for an adapted land use management, which has great potential to mitigate degradation processes related to climate change. Full article

The Sistan Basin has been recognized as one of the most active dust sources and windiest desert environments in the world. Although the dust activity in Sistan maximizes during the summer, rare but intense dust storms may also occur in the winter. This study aims to elucidate the atmospheric dynamics related to dust emission and transport, dust-plume characteristics, and impacts on aerosol properties and air quality during an intense dust storm over Sistan in February 2019. The dust storm was initiated by strong northerly winds (~20 ms⁻¹) associated with the intrusion of a cold front from high latitudes. The upper-level potential vorticity (PV)-trough evolved into a cut-off low in the mid and upper troposphere and initiated unstable weather over Afghanistan and northern Pakistan. At the surface, density currents emanating from deep convective clouds and further strengthened by downslope winds from the mountains, caused massive soil erosion. The passage of the cold front reduced the temperature by ~10 °C and increased the atmospheric pressure by ~10 hPa, while the visibility was limited to less than 200 m. The rough topography played a major role in modulating the atmospheric dynamics, wind field, dust emissions, and transport pathways. Meso-NH model simulates large amounts of columnar mass dust loading (> 20 g m⁻²) over Sistan, while the intense dust plume was mainly traveling below 2 km and increased the particulate matter (PM₁₀) concentrations up to 1800 µg m⁻³ at Zabol. The dust storm was initially moving in an arc-shaped pathway over the Sistan Basin and then it spread away. Plumes of dust covered a large area in southwest Asia, reaching the northern Arabian Sea, and the Thar desert one to two days later, while they strongly affected the aerosol properties at Karachi, Pakistan, by increasing the aerosol optical depth (AOD > 1.2) and the coarse-mode fraction at ~0.7. Full article

Local scour is a dynamic process evolving during the lifetime of bridges as a result of the changes in hydrologic and hydraulic conditions. Current approaches for scour risk assessment are generally based on the evaluation of the equilibrium scour depth for a flood event with a prefixed return period. The temporal evolution of the bridge-pier scour process is usually disregarded, by assuming that equilibrium conditions are always attained, regardless of the flood properties. However, recent studies have highlighted the importance of accounting for the contribution of multiple flood events and their exact hydrograph shape. This study aims at quantifying the epistemic uncertainty related to the modelling of the temporal evolution of scour under multiple consecutive flood events in clear-water conditions. A simple numerical case study is considered, using a Markovian framework to describe probabilistically the progression of scour. Well-known time-dependent scour models are used to estimate the temporal evolution of the scour-depth under each flood hydrograph, and the scour estimates are compared with those obtained using widely employed equilibrium scour formulas. Results show that the expected scour depth is influenced by the parameters used to describe the flood hydrograph and that the probability distribution of the scour depth is highly sensitive to the choice of the time-dependent scour model. The uncertainty in the scour estimation stemming from the formula adopted in this study for describing the temporal evolution of the scour depth can be higher than those related to the formula adopted for equilibrium scour. Full article

Despite landslide inventories being compiled throughout the world every year at different scales, limited efforts have been made to critically compare them using various techniques or by different investigators. Event-based landslide inventories indicate the location, distribution, and detected boundaries of landslides caused by a single event, such as an earthquake or a rainstorm. Event-based landslide inventories are essential for landslide susceptibility mapping, hazard modeling, and further management of risk mitigation. In Nepal, there were several attempts to map landslides in detail after the Gorkha earthquake. Particularly after the main event on 25 April 2015, researchers around the world mapped the landslides induced by this earthquake. In this research, we compared four of these published inventories qualitatively and quantitatively using different techniques. Two principal methodologies, namely the cartographical degree of matching and frequency area distribution (FAD), were optimized and applied to evaluate inventory maps. We also showed the impact of using satellite imagery with different spatial resolutions on the landslide inventory generation by analyzing matches and mismatches between the inventories. The results of our work give an overview of the impact of methodology selection and outline the limitations and advantages of different remote sensing and mapping techniques for landslide inventorying. Full article

The San Joaquin Valley and Tulare basins in California’s Central Valley have intensive agricultural activity and groundwater demand that has caused significant subsidence and depletion of water resources in the past. We measured groundwater pumping-induced land subsidence in the southern Central Valley from March 2015 to May 2017 using Sentinel-1 interferometric synthetic aperture radar (InSAR) data. The InSAR measurements provided fine spatial details of subsidence patterns and displayed a superposition of secular and seasonal variations that were coherent across our study region and correlated with precipitation variability and changes in freshwater demand. Combining InSAR and Global Positioning System (GPS) data, precipitation, and in situ well records showed a broad scale slowdown/cessation of long term subsidence in the wetter winter of 2017, likely reflecting the collective response of the Central Valley aquifer system to heavier-than-usual precipitation. We observed a very good temporal correlation between the Gravity Recovery and Climate Experiment (GRACE) satellite groundwater anomaly (GWA) variation and long-term subsidence records, regardless of local hydrogeology and mechanical properties. This indicates the subsidence from satellite geodesy is a very useful indicator for tracking groundwater storage change. With the continuing acquisition of Sentinel-1 and other satellites, we anticipate decadal-scale subsidence records with a spatial resolution of tens to hundreds of meters will be available in the near future to be combined with basin-averaged GRACE measurements to improve our estimate of time-varying groundwater change. Full article

Waste management (WM) is a demanding undertaking in all countries, with important implications for human health, environmental preservation, sustainability and circular economy. The method of sanitary landfilling for final disposal of waste remains a generally accepted and used method but the available scientific evidence on the waste-related environmental and health effects is not conclusive. Comparative studies of various WM methods (landfilling, incineration, composting etc.) show that among the municipal solid waste (MSW) treatment and disposal technological options, sanitary landfilling or open dumping is popular in most countries because of the relative low cost and low-technical requirement. The European Union (EU) Directive on waste landfills has introduced specific goals for reducing the volume of disposed waste and very strict requirements for landfilling and landfill sites. Evaluation of the impact of landfills on the environment is a crucial topic in the literature and has received increased attention recently, given growing environmental concerns. The main goal of this survey was to conduct a comprehensive assessment of possible impacts of MSW landfills on the environment. The main conclusion of the overall assessment of the literature is that the disposal of MSW in landfills entails a number of environmental risks but with respect to the current situation and rich style of living adopted in industrially developed countries, the idea of WM systems functioning without landfilling—at least in the foreseeable future within one generation—seems to be somewhat unreal. The results also provided important information of landfills as a source of environmental risk. Results of this research may have an important impact on landfill management and the disposal of waste. From the literature review, it is evident that even if high levels of waste avoidance, reuse and recycling are achieved, some waste materials will always need to be forwarded for disposal. Full article

Satellite gravimetry data acquired by the Gravity Recovery and Climate Experiment (GRACE) allows to derive the temporal evolution in ice mass for both the Antarctic Ice Sheet (AIS) and the Greenland Ice Sheet (GIS). Various algorithms have been used in a wide range of studies to generate Gravimetric Mass Balance (GMB) products. Results from different studies may be affected by substantial differences in the processing, including the applied algorithm, the utilised background models and the time period under consideration. This study gives a detailed description of an assessment of the performance of GMB algorithms using actual GRACE monthly solutions for a prescribed period as well as synthetic data sets. The inter-comparison exercise was conducted in the scope of the European Space Agency’s Climate Change Initiative (CCI) project for the AIS and GIS, and was, for the first time, open to everyone. GMB products generated by different groups could be evaluated and directly compared against each other. For the period from 2003-02 to 2013-12, estimated linear trends in ice mass vary between −99 Gt/yr and −108 Gt/yr for the AIS and between −252 Gt/yr and −274 Gt/yr for the GIS, respectively. The spread between the solutions is larger if smaller drainage basins or gridded GMB products are considered. Finally, findings from the exercise formed the basis to select the algorithms used for the GMB product generation within the AIS and GIS CCI project. Full article

Landslides over steep slopes, floods along rivers plains and debris flows across valleys are hydrogeological phenomena typical for mountain regions. Such events are generally triggered by rainfall, which can have large variability in terms of both its intensity and volume. Furthermore, terrain predisposition and the presence of some disturbances, such as wildfires, can have an adverse effect on the potential risk. Modelling the complex interaction between these components is not a simple task and cannot always be carried out using instability thresholds that only take into account the characteristics of the rainfall events. In some particular cases, external factors can modify the existing delicate equilibrium on the basis of which stability thresholds are defined. In particular, events such as wildfires can cause the removal of vegetation coverage and the modification of the soil terrain properties. Therefore, wildfires can effectively reduce the infiltration capacity of the terrain and modify evapotranspiration. As a result, key factors for slope stability, such as the trend of the degree of saturation of the terrain, can be strongly modified. Thus, studying the role of wildfire effects on the terrain’s hydrological balance is fundamental to establish the critical conditions that can trigger potential slope failures (i.e., shallow landslides and possible subsequent debris flows). In this work, we investigate the consequences of wildfire on the stability of slopes through a hydrological model that takes into account the wildfire effects and compare the results to the current stability thresholds. Two case studies in the Ardenno (IT) and Ronco sopra Ascona (CH) municipalities were chosen for model testing. The aim of this paper is to propose a quantitative analysis of the two cases studies, taking into account the role of fire in the slope stability assessment. The results indicate how the post-fire circumstances strongly modify the ability of the terrain to absorb rainfall water. This effect results in a persistently drier terrain until a corner point is reached, after which the stability of the slope could be undermined by a rainfall event of negligible intensity. Full article

We digitize surface rupture maps and compile observational data from 67 publications on ten of eleven historical, surface-rupturing earthquakes in Australia in order to analyze the prevailing characteristics of surface ruptures and other environmental effects in this crystalline basement-dominated intraplate environment. The studied earthquakes occurred between 1968 and 2018, and range in moment magnitude (Mw) from 4.7 to 6.6. All earthquakes involved co-seismic reverse faulting (with varying amounts of strike-slip) on single or multiple (1–6) discrete faults of ≥ 1 km length that are distinguished by orientation and kinematic criteria. Nine of ten earthquakes have surface-rupturing fault orientations that align with prevailing linear anomalies in geophysical (gravity and magnetic) data and bedrock structure (foliations and/or quartz veins and/or intrusive boundaries and/or pre-existing faults), indicating strong control of inherited crustal structure on contemporary faulting. Rupture kinematics are consistent with horizontal shortening driven by regional trajectories of horizontal compressive stress. The lack of precision in seismological data prohibits the assessment of whether surface ruptures project to hypocentral locations via contiguous, planar principal slip zones or whether rupture segmentation occurs between seismogenic depths and the surface. Rupture centroids of 1–4 km in depth indicate predominantly shallow seismic moment release. No studied earthquakes have unambiguous geological evidence for preceding surface-rupturing earthquakes on the same faults and five earthquakes contain evidence of absence of preceding ruptures since the late Pleistocene, collectively highlighting the challenge of using mapped active faults to predict future seismic hazards. Estimated maximum fault slip rates are 0.2–9.1 m Myr⁻¹ with at least one order of uncertainty. New estimates for rupture length, fault dip, and coseismic net slip can be used to improve future iterations of earthquake magnitude—source size—displacement scaling equations. Observed environmental effects include primary surface rupture, secondary fracture/cracks, fissures, rock falls, ground-water anomalies, vegetation damage, sand-blows/liquefaction, displaced rock fragments, and holes from collapsible soil failure, at maximum estimated epicentral distances ranging from 0 to ~250 km. ESI-07 intensity-scale estimates range by ± 3 classes in each earthquake, depending on the effect considered. Comparing Mw-ESI relationships across geologically diverse environments is a fruitful avenue for future research. Full article

Alakol Lake is one of the largest hydrologically closed lake located in Balkash-Alakol River Basin in southeast Kazakhstan. Having a coastline approximately at 490 km, Alakol Lake has faced multiple threats due to both natural and anthropogenic factors as a result of tectonic movements, geology, wind-wave conditions, growing tourism activities, fishing, and transport, etc. The present study aims to investigate the historical trends in coastline changes along Alakol Lake in Kazakhstan and estimate its change rate by using remote sensing data in particular scale-space images Landsat-5 TM, 7 ETM+, 8 OLI, and Sentinel-2A. Based on Landsat and Sentinel data, the modified normalized difference water index was calculated to demonstrate the coastline changes along Alakol Lake between 1990 and 2018. Moreover, the monitoring and analysis of coastline dynamics is based on the main morphometric characteristics of Alakol Lake including water surface area, coastline length, geomorphology of the coast, etc. Our results reveal that there is a continuous coastline retreat, depending on the coast types. For example, in the case of the denudation coasts, a land inundation was from 120 to 270 m between 1990 and 2018. In the case of the accumulative coast (mainly northeast, north, and northwest coasts) a land inundation was from 200 to 900 m. A vast area of agricultural land around Alakol Lake become flooded and lost. This study demonstrates the importance of monitoring coastline dynamics because it provides essential information for understanding the coastal response to contemporary nature and anthropogenic impacts. Full article

Persistent non-explosive passive degassing is a common characteristic of active volcanoes. Distinct periodic components in measurable parameters of gas release have been widely identified over timescales ranging from seconds to months. The development and implementation of high temporal resolution gas measurement techniques now enables the robust quantification of high frequency processes operating on timescales comparable to those detectable in geophysical datasets. This review presents an overview of the current state of understanding regarding periodic volcanic degassing, and evaluates the methods available for detecting periodicity, e.g., autocorrelation, variations of the Fast Fourier Transform (FFT), and the continuous wavelet transform (CWT). Periodicities in volcanic degassing from published studies were summarised and statistically analysed together with analyses of literature-derived datasets where periodicity had not previously been investigated. Finally, an overview of current knowledge on drivers of periodicity was presented and discussed in the framework of four main generating categories, including: (1) non-volcanic (e.g., atmospheric or tidally generated); (2) gas-driven, shallow conduit processes; (3) magma movement, intermediate to shallow storage zone; and (4) deep magmatic processes. Full article

The first 10 days of September 2015 were marked by intense dust activity over the Middle East and the Arabian Peninsula. This study examines the atmospheric conditions at the origin of the large dust storms during this period. We particularly investigate the atmospheric dynamics leading to the development of a large dry cyclone over Iraq on 31 August 2015 which in turn generated an intense dust storm that affected most of the countries around the Arabian Gulf and lasted for 5 days. We found that the cyclone developed over Northwest Iraq as a transfer to low levels of a cut-off low which had formed two days earlier at upper levels over Turkey. Large dust loads exceeding 250 tons were emitted and moved southeast in a cyclonic shape toward the Arabian sea. The second large dust storm on 6-8 September 2015 occurred over Syria and affected all the coastal countries on the eastern side of the Mediterranean Sea. It was associated with the occurrence of a series of density currents over northeast Syria emanating from deep convection over the mountainous border between Syria and Turkey. The unusual development of deep convection over this area was associated with a blocking high and interaction with orography. Both the cut-off high and the cut-off low occurred during a period characterized by a meandering polar jet and an enhanced subtropical jet causing unstable weather over mid-latitudes which in turn led to highly polluted atmosphere by natural dust in the affected countries. Full article

Expansive soil shows dual swell–shrink which is not suitable for construction. Several mitigating techniques exist to counteract the problem promulgated by expansive clayey soils. This paper explored the potential mecho-chemical reinforcement of expansive clayey soil to mitigate the effect of upward swelling pressure and heave. The polypropylene fiber is randomly distributed in the soil for mechanical stabilization, and the industrial residual silica fume is used as a chemical stabilizer. The experimental analysis was made in three phases which involved tests on mechanically-reinforced expansive soil, using randomly distributed polypropylene fibers with different percentages (0.25%, 0.50%, and 1.00%), and which were 12 mm length. The second phase of experiments was carried out on chemical stabilized expansive soil with different percentages (2%, 4%, and 8%) of silica, and the next phase of the experiment focused on the combination of mecho-chemical stabilization of the expansive soil with different combinations of silica (i.e., 2%, 4%, and 8%) and polypropylene fibers (i.e., 0.25%, 0.50%, and 1.00%). Maximum dry density (MDD), optimum moisture content (OMC), liquid limit (LL), plastic limit (PL), plastic index (PI), grain size, and constant volume swelling pressure tests were performed on unreinforced and reinforced expansive soil, to investigate the effects of polypropylene fiber and silica fume on the engineering properties of expansive clayey soil. The experimental results illustrate that the inclusion of polypropylene fiber has a significant effect on the upward swelling pressure and expansion property of expansive soil. The reduction in the upward swelling pressure and expansion is a function of fiber content. These results also indicated that the use of silica fume caused a reduction in upward swelling potential, and its effect was considerably more than the influence of fiber. Full article

We provide a dataset of the landslides induced by the 2016 Pedernales megathrust earthquake, Ecuador (Mw 7.8, focal depth of 20 km) and compare their spatial distribution with mapped bedrock lithology, horizontal peak ground acceleration (PGA-h) and the macroseismic intensity based on earthquake-induced environmental effects (ESI-07). We studied 192 coseismic landslides (classified as coherent, disrupted and lateral spreads) located in the epicentral area, defined by the VII to IX_ESI-07 isoseismals. Based on our findings, lahar deposits, tuffs and volcanoclastic units are the most susceptible to landslides occurrence. Alluvial plains with fluvial loose fine sand are the most susceptible setting for lateral spreading, with a maximum intensity of IX_ESI-07. The coherent landslides are frequently found in altered shale and siltstone geological units with moderate slopes (8°–16°), with typical intensity ranging between VII and VIII_ESI-07. Our analysis draws a typical framework for slope movements triggered by subduction earthquakes in Ecuador. The most dangerous setting is the coastal region, a relatively highly urbanized area located near the epicenter and where liquefaction can trigger massive lateral spreading events. Coherent and disrupted landslides, dominating the more internal hilly region, can be triggered also in moderate slope settings (i.e., less than 10°). Indeed, the regression analysis between seismic intensity, PGA-h and landslide occurrence shows that most of the events occurred at PGA-h values between 0.4 g and 1.2 g, at a distance of 30 to 50 km from the rupture plane. Our database suggests that lithology and hillslope geometry are the main geological/geomorphological factors controlling coseismic landslides occurrence; while the distance from the rupture plane plays a significant role on determining the landslide size. Finally, we underline that coseismically-triggered landslides are among the most common environmental effects occurring during large subduction events that can be effectively used to properly evaluate the earthquake macroseismic field. The landslide inventory we compiled is suitable for assessing the vulnerability of physical environment from subduction earthquakes in Ecuador, and offers a primary data source for future worldwide analysis. Full article

Morphologically diverse organo-sedimentary structures (including microbial mats and stromatolites) provide a palaeobiological record through more than three billion years of Earth history. Since understanding much of the Archaean fossil record is contingent upon proving the biogenicity of such structures, mechanistic interpretations of well-preserved fossil microbialites can reinforce our understanding of their biogeochemistry and distinguish unambiguous biological characteristics in these structures, which represent some of the earliest records of life. Mechanistic morphogenetic understanding relies upon the analysis of geomicrobiological experiments. Herein, we report morphological-biogeochemical comparisons between micromorphologies observed in growth experiments using photosynthetic mats built by the cyanobacterium Coleofasciculus chthonoplastes (formerly Microcoleus) and green anoxygenic phototrophic Chloroflexus spp. (i.e., Coleofasciculus–Chloroflexus mats), and Precambrian organo-sedimentary structures, demonstrating parallels between them. In elevated ambient concentrations of Cu (toxic to Coleofasciculus), Coleofasciculus–Chloroflexus mats respond by forming centimetre-scale pinnacle-like structures (supra-lamina complexities) associated with large quantities of EPS at their surfaces. µPIXE mapping shows that Cu and other metals become concentrated within surficial sheath-EPS-Chloroflexus-rich layers, producing density-differential micromorphologies with distinct fabric orientations that are detectable using X-ray computed micro-tomography (X-ray µCT). Similar micromorphologies are also detectable in stromatolites from the 3.481 Ga Dresser Formation (Pilbara, Western Australia). The cause and response link between the presence of toxic elements (geochemical stress) and the development of multi-layered topographical complexities in organo-sedimentary structures may thus be considered an indicator of biogenicity, being an indisputably biological and predictable morphogenetic response reflecting, in this case, the differential responses of Coleofasciculus and Chloroflexus to Cu. Growth models for microbialite morphogenesis rely upon linking morphology to intrinsic (biological) and extrinsic (environmental) influences. Since the pinnacles of Coleofasciculus–Chloroflexus mats have an unambiguously biological origin linked to extrinsic geochemistry, we suggest that similar micromorphologies observed in ancient organo-sedimentary structures are indicative of biogenesis. An identical Coleofasciculus–Chloroflexus community subjected to salinity stress also produced supra-lamina complexities (tufts) but did not produce identifiable micromorphologies in three dimensions since salinity seems not to negatively impact either organism, and therefore cannot be used as a morphogenetic tool for the interpretation of density-homogeneous micro-tufted mats—for example, those of the 3.472 Ga Middle Marker horizon. Thus, although correlative microscopy is the keystone to confirming the biogenicity of certain Precambrian stromatolites, it remains crucial to separately interrogate each putative trace of ancient life, ideally using three-dimensional analyses, to determine, where possible, palaeoenvironmental influences on morphologies. Widespread volcanism and hydrothermal effusion into the early oceans likely concentrated toxic elements in early biomes. Morphological diversity in fossil microbialites could, therefore, reflect either (or both of) differential exposure to ambient fluids enriched in toxic elements and/or changing ecosystem structure and tolerance to elements through evolutionary time—for example, after incorporation into enzymes. Proof of biogenicity by deducing morphogenesis (i.e., a process preserved in the fossil record) overcomes many of the shortcomings inherent to the proof of biogenicity by descriptions of morphology alone. Full article

Unmanned Aerial Vehicles (UAV) are a rapidly evolving tool in geosciences and are increasingly deployed for studying the dynamic processes of the earth’s surface. To assess the potential of autonomous low-cost UAVs for the mapping and monitoring of alpine glaciers, we conducted multiple aerial surveys on the Kanderfirn in the Swiss Alps in 2017 and 2018 using open hardware and software of the Paparazzi UAV project. The open-source photogrammetry software OpenDroneMap was tested for the generation of high-resolution orthophotos and digital surface models (DSMs) from aerial imagery and cross-checked with the well-established proprietary software Pix4D. Accurately measured ground control points served for the determination of the geometric accuracy of the orthophotos and DSMs. A horizontal (xy) accuracy of 0.7–1.2 m and a vertical (z) accuracy of 0.7–2.1 m was achieved for OpenDroneMap, compared to a xy-accuracy of 0.3–0.5 m and a z-accuracy of 0.4–0.5 m obtained for Pix4D. Based on the analysis and comparison of different orthophotos and DSMs, surface elevation, roughness and brightness changes from 3 June to 29 September 2018 were quantified. While the brightness of the glacier surface decreased linearly over the ablation season, the surface roughness increased. The mean DSM-based elevation change across the glacier tongue was 8 m, overestimating the measured melting and surface lowering at the installed ablation stakes by about 1.5 m. The presented results highlight that self-built fixed-wing UAVs in tandem with open-source photogrammetry software are an affordable alternative to commercial remote-sensing platforms and proprietary software. The applied low-cost approach also provides great potential for other regions and geoscientific disciplines. Full article

Along with rapid population growth in Vietnam, there is an increasing dependence on groundwater for various activities. An Giang province is known to be one of the agricultural intensification areas of The Vietnamese Mekong Delta (VMD). This study aimed to evaluate the spatiotemporal variation of groundwater quality for a period of ten years from 2009 to 2018 in An Giang. The weighted groundwater quality index (GWQI) was developed based on the fuzzy analytic hierarchy process (Fuzzy-AHP) for assigning weighted parameters. The results show that that shallow wells in the Northeast and Southeast regions of An Giang were mostly categorized under “bad water” quality with high arsenic (As) concentration over the years partly due to huge amounts of sediment deposition in monsoon season. Overall, the reason for the poor groundwater quality in An Giang was the combined effect of both natural and human activities. On the other hand, we detected high values of GWQI links with high As concentration in areas where people extract more groundwater for irrigation. Temporal variation of GWQI suggested that groundwater quality at eight wells has improved from 2009 to 2018 in the wet season as compared to the dry season. The reason behind the improvement of groundwater quality during wet season was the decrease in river discharge, which causes less deposition of suspended solids near the flood plains. Moreover, the filling of unused wells can reduce the movement of pollutants from unused wells to groundwater aquifers. Although there was not sufficient evidence to show the relationship between As and sediment concentration, the temporal reduction trend in river discharge and suspended solids was detected in An Giang. The understanding of groundwater quality can help policymakers protect and manage limited water resources in the long-term. Full article

The Temple of Athena is one of the main sacred areas of the Greek–Roman settlement of Poseidonia-Paestum (southern Italy). Several archaeological excavations were carried out here between the late nineteenth and early twentieth century. Unfortunately, the locations of these excavations are only approximately known, as are the geomorphology and stratigraphy of the temple area. A multidisciplinary study, including stratigraphic, geomorphological, archaeological, and sedimentological investigations, remote sensing, and electromagnetic and geoelectrical tests, was therefore carried out, shedding new light on the geomorphology and stratigraphy of the SW and W temple sectors. The geophysical data obtained revealed anomalies in the subsoil that probably correspond to ancient structures and the cutting of the travertine deposits around the temple. The position and extension of the trenches of the early archaeological excavations were also established. Full article

Rock brittleness is pivotal in the development of the unconventional reservoirs. However, the existence of various methods of calculating the brittleness index (BI) such as the mineral-based brittleness index (MBI), the log-based brittleness index (LBI), and the elastic-based brittleness index (EBI) lead to inconclusive estimations of the brittleness index. Hence, in this work, the existing correlations are applied on prolific unconventional plays in the U.S. such as the Marcellus, Bakken, Niobrara, and Chattanooga Formation to examine the various BI methods. A detailed comparison between the MBI, LBI, and EBI has also been conducted. The results show that a universal correlation cannot be derived in order to define brittleness since it is a function of lithology. Correlation parameters vary significantly from one shale play to another. Nevertheless, an overall trend shows that abundant quartz and carbonates content yield high brittleness values, while the high clay content and porosity lower the rock brittleness. Full article

Recent and historical studies of earthquake-induced liquefaction, as well as paleoliquefaction studies, demonstrate the potential usefulness of liquefaction data in the assessment of the earthquake potential of seismic sources. Paleoliquefaction studies, along with other paleoseismology studies, supplement historical and instrumental seismicity and provide information about the long-term behavior of earthquake sources. Paleoliquefaction studies focus on soft-sediment deformation features, including sand blows and sand dikes, which result from strong ground shaking. Most paleoliquefaction studies have been conducted in intraplate geologic settings, but a few such studies have been carried out in interplate settings. Paleoliquefaction studies provide information about timing, location, magnitude, and recurrence of large paleoearthquakes, particularly those with moment magnitude, M, greater than 6 during the past 50,000 years. This review paper presents background information on earthquake-induced liquefaction and resulting soft-sediment deformation features that may be preserved in the geologic record, best practices used in paleoliquefaction studies, and application of paleoliquefaction data in earthquake source characterization. The paper concludes with two examples of regional paleoliquefaction studies—in the Charleston seismic zone and the New Madrid seismic zone in the southeastern and central United States, respectively—which contributed to seismic source models used in earthquake hazard assessment. Full article

Landslides are one of the most devastating and commonly recurring natural hazards in the Indian Himalayas. They contribute to infrastructure damage, land loss and human casualties. Most of the landslides are primarily rainfall-induced and the relationship has been well very well-established, having been commonly defined using empirical-based models which use statistical approaches to determine the parameters of a power-law equation. One of the main drawbacks using the traditional empirical methods is that it fails to reduce the uncertainties associated with threshold calculation. The present study overcomes these limitations by identifying the precipitation condition responsible for landslide occurrence using an algorithm-based model. The methodology involves the use of an automated tool which determines cumulated event rainfall–rainfall duration thresholds at various exceedance probabilities and the associated uncertainties. The analysis has been carried out for the Kalimpong Region of the Darjeeling Himalayas using rainfall and landslide data for the period 2010–2016. The results signify that a rainfall event of 48 hours with a cumulated event rainfall of 36.7 mm can cause landslides in the study area. Such a study is the first to be conducted for the Indian Himalayas and can be considered as a first step in determining more reliable thresholds which can be used as part of an operational early-warning system. Full article

Climate reanalyses provide key information to calibrate proxy records in regions with scarce direct observations. The climate reanalysis used to perform a proxy calibration should accurately reproduce the local climate variability. Here we present a regional scale evaluation of meteorological parameters using ERA-Interim and ERA5 reanalyses compared to in-situ observations from 13 automatic weather stations (AWS), located in the southern Antarctic Peninsula and Ellsworth Land, Antarctica. Both reanalyses seem to perform better in the escarpment area (>1000 m a.s.l) than on the coast. A significant improvement is observed in the performance of ERA5 over ERA-Interim. ERA5 is highly accurate, representing the magnitude and variability of near-surface air temperature and wind regimes. The higher spatial and temporal resolution provided by ERA5 reduces significantly the cold coastal biases identified in ERA-Interim and increases the accuracy representing the wind direction and wind speed in the escarpment. The slight underestimation in the wind speed obtained from the reanalyses could be attributed to an interplay of topographic factors and the effect of local wind regimes. Three sites in this region are highlighted for their potential for ice core studies. These sites are likely to provide accurate proxy calibrations for future palaeoclimatic reconstructions. Full article

This study examines the history of the paradigm concerning a lunar (or solar-system-wide) terminal cataclysm (also called “Late Heavy Bombardment” or LHB), a putative, brief spike in impacts at ~3.9 Ga ago, preceded by low impact rates. We examine origin of the ideas, why they were accepted, and why the ideas are currently being seriously revised, if not abandoned. The paper is divided into the following sections: Overview of paradigm. Pre-Apollo views (1949–1969). Initial suggestions of cataclysm (ca. 1974). Ironies. Alternative suggestions, megaregolith evolution (1970s). Impact melt rocks “establish” cataclysm (1990). Imbrium redux (ca. 1998). Impact melt clasts (early 2000s). Dating of front-side lunar basins? Dynamical models “explain” the cataclysm (c. 2000s). Asteroids as a test case. Impact melts predating 4.0 Ga ago (ca. 2008–present.). Biological issues. Growing doubts (ca. 1994–2014). Evolving Dynamical Models (ca. 2001–present). Connections to lunar origin. Dismantling the paradigm (2015–2018). “Megaregolith Evolution Model” for explaining the data. Conclusions and new directions for future work. Full article

This paper examined three different types of recycled materials, such as beer green glass, waste tile, and asphalt, which will be used in different mixtures in order to prepare concrete specimens and, more specifically, their effect on concrete strength and how the petrographic characteristics of various recycled materials influenced the durability of C25/30 strength class concrete. Particular emphasis was placed on the effect of artificial microroughness of glassy and smooth surfaces of recycled materials on their final concrete strength. The concrete strength values do not show great variance, but their limited differences have been qualitatively interpreted by a new promising petrographic methodology, including the study of the surface texture of the used aggregate materials. Concretes are produced with constant volume proportions, workability, mixing, and curing conditions while using different sizes of each aggregate type. The aggregates were mixed both in dry and water saturated states in concretes. Concretes that are made by a mixture of beer green glass with quartz primer, as well as of tile with quartz primer, presented the optimum possible results of the compressive strength. Full article

Dust storms present numerous hazards to human society and are particularly significant to people living in the Dust Belt which stretches from the Sahara across the Middle East to northeast Asia. This paper presents a review of dust storm variability and trends in frequency on decadal timescales from three Dust Belt settlements with long-term (>50 years) meteorological records: Nouakchott, Mauritania; Zabol, Iran, and Minqin, China. The inhabitants of each of these settlements have experienced a decline in dust storms in recent decades, since the late 1980s at Nouakchott, since 2004 at Zabol, and since the late 1970s at Minqin. The roles of climatic variables and human activities are assessed in each case, as drivers of periods of high dust storm frequency and subsequent declines in dust emissions. Both climatic and human variables have been important but overall the balance of research conclusions indicates natural processes (precipitation totals, wind strength) have had greater impact than human action, in the latter case both in the form of mismanagement (abandoned farmland, water management schemes) and attempts to reduce wind erosion (afforestation projects). Understanding the drivers of change in dust storm dynamics at the local scale is increasingly important for efforts to mitigate dust storm hazards as climate change projections suggest that the global dryland area is likely to expand in the twenty-first century, along with an associated increase in the risk of drought and dust emissions. Full article

This paper summarizes current understanding of the processes that determine the dynamics of the subsea permafrost–hydrate system existing in the largest, shallowest shelf in the Arctic Ocean; the East Siberian Arctic Shelf (ESAS). We review key environmental factors and mechanisms that determine formation, current dynamics, and thermal state of subsea permafrost, mechanisms of its destabilization, and rates of its thawing; a full section of this paper is devoted to this topic. Another important question regards the possible existence of permafrost-related hydrates at shallow ground depth and in the shallow shelf environment. We review the history of and earlier insights about the topic followed by an extensive review of experimental work to establish the physics of shallow Arctic hydrates. We also provide a principal (simplified) scheme explaining the normal and altered dynamics of the permafrost–hydrate system as glacial–interglacial climate epochs alternate. We also review specific features of methane releases determined by the current state of the subsea-permafrost system and possible future dynamics. This review presents methane results obtained in the ESAS during two periods: 1994–2000 and 2003–2017. A final section is devoted to discussing future work that is required to achieve an improved understanding of the subject. Full article

The paper addresses the fields of combined close-range photogrammetry and terrestrial laser scanning in the light of ship modelling. The authors pointed out precision and measurement accuracy due to their possible complex application for ship hulls inventories. Due to prescribed vitality of every ship structure, it is crucial to prepare documentation to support the vessel processes. The presented methods are directed, combined photogrammetric techniques in ship hull inventory due to submarines. The class of photogrammetry techniques based on high quality photos are supposed to be relevant techniques of the inventories’ purpose. An innovative approach combines these methods with Terrestrial Laser Scanning. The process stages of data acquisition, post-processing, and result analysis are presented and discussed due to market requirements. Advantages and disadvantages of the applied methods are presented. Full article

In the present paper, we analyse ground tilt and seismicity at Campi Flegrei caldera and Ischia Island, two volcanic areas located in the south of Italy. These areas have been well studied for many years from a petrological, volcanological and geophysical view point. Moreover, due to the high seismic and volcanic risk for the populations living there, they are continuously monitored by networks of geophysical and geochemical sensors. We summarize the most important results that we obtained so far, concerning the observations of relationships between seismic activity and ground tilt anomalies, focusing on the time interval 2015–2018. First, we present a detailed description of the tiltmeter and seismic networks in both the investigated areas, as well as their development and improvement over time that has enabled high quality data collection. From the joint analysis of the seismic and borehole tiltmeter signals, we often notice concurrence between tilt pattern variations and the occurrence of seismicity. Moreover, the major tilt anomalies appear to be linked with the rate and energy of volcano-tectonic earthquakes, as well as with exogenous phenomena like solid Earth tides and hydrological cycles. The analysis that we present has potential applicability to other volcanic systems. Our findings show how the joint use tilt and seismic data can contribute to better understanding of the dynamics of volcanoes. Full article

This brief note aims to describe the history, from its early original idea, of the new macroseismic scale: The Environmental Seismic Intensity Scale 2007 (ESI 2007). It can be used together with other existing scales or alone when needed for measuring the intensity of an earthquake on the basis of the primary and secondary effects of a seismic event on the natural environment. These effects could be the major sources of earthquake hazards, as recently proved. This note also aims to contribute to the understanding of processes that induced the researcher to develop an idea, to pursue it, and bring it to its end, first through the help of valuable Italian researchers and then through the constructive exchange of ideas with researchers of different cultural backgrounds operating almost everywhere in the world. This note is sponsored and approved by the International Union for Quaternary Research (INQUA), and the Environmental Seismic Intensity scale (ESI-07) was published in 2007 after a revision process of about eight years. Full article

In 2011 Brazil experienced the worst disaster in the country’s history. There were 918 deaths and thousands made homeless in the mountainous region of Rio de Janeiro State due to several landslides triggered by heavy rainfalls. This area constantly suffers high volumes of rain and episodes of landslides. Due to these experiences, we used the MaCumBa (Massive CUMulative Brisk Analyser) software to identify rainfall intensity–duration thresholds capable of triggering landslides in the most affected municipalities of this region. More than 3000 landslides and rain data from a 10-year long dataset were used to define the thresholds and one year was used to validate the results. In this work, a set of three thresholds capable of defining increasing alert levels (moderate, high and very high) has been defined for each municipality. Results show that such thresholds may be used for early alerts. In the future, the same methodology can be replicated to other Brazilian municipalities with different datasets, leading to more accurate warning systems. Full article

Depleted gas reservoirs are appealing targets for carbon dioxide (CO ${}_2$ ) sequestration because of their storage capacity, proven seal, reservoir characterization knowledge, existing infrastructure, and potential for enhanced gas recovery. Low abandonment pressure in the reservoir provides additional voidage-replacement potential for CO ${}_2$ and allows for a low surface pump pressure during the early period of injection. However, the injection process poses several challenges. This work aims to raise awareness of key operational challenges related to CO ${}_2$ injection in low-pressure reservoirs and to provide a new approach to assessing the phase behavior of CO ${}_2$ within the wellbore. When the reservoir pressure is below the CO ${}_2$ bubble-point pressure, and CO ${}_2$ is injected in its liquid or supercritical state, CO ${}_2$ will vaporize and expand within the well-tubing or in the near-wellbore region of the reservoir. This phenomenon is associated with several flow assurance problems. For instance, when CO ${}_2$ transitions from the dense-state to the gas-state, CO ${}_2$ density drops sharply, affecting the wellhead pressure control and the pressure response at the well bottom-hole. As CO ${}_2$ expands with a lower phase viscosity, the flow velocity increases abruptly, possibly causing erosion and cavitation in the flowlines. Furthermore, CO ${}_2$ expansion is associated with the Joule–Thomson (IJ) effect, which may result in dry ice or hydrate formation and therefore may reduce CO ${}_2$ injectivity. Understanding the transient multiphase phase flow behavior of CO ${}_2$ within the wellbore is crucial for appropriate well design and operational risk assessment. The commonly used approach analyzes the flow in the wellbore without taking into consideration the transient pressure response of the reservoir, which predicts an unrealistic pressure gap at the wellhead. This pressure gap is related to the phase transition of CO ${}_2$ from its dense state to the gas state. In this work, a new coupled approach is introduced to address the phase behavior of CO ${}_2$ within the wellbore under different operational conditions. The proposed approach integrates the flow within both the wellbore and the reservoir at the transient state and therefore resolves the pressure gap issue. Finally, the energy costs associated with a mitigation process that involves CO ${}_2$ heating at the wellhead are assessed. Full article

On 15th November 2017, the Pohang earthquake (M_w 5.4) had strong ground shaking that caused severe liquefaction and lateral spreading across the Heunghae Basin, around Pohang city, South Korea. Such liquefaction is a rare phenomenon during small or moderate earthquakes (M_W < 5.5). There are only a few examples around the globe, but more so in the Korean Peninsula. In this paper, we present the results of a systematic survey of the secondary ground effects—i.e., soil liquefaction and ground cracks—developed during the earthquake. Most of the liquefaction sites are clustered near the epicenter and close to the Heunghae fault. Based on the geology, tectonic setting, distribution, and clustering of the sand boils along the southern part of the Heunghae Basin, we propose a geological model, suggesting that the Heunghae fault may have acted as a barrier to the propagation of seismic waves. Other factors like the mountain basin effect and/or amplification of seismic waves by a blind thrust fault could play an important role. Liquefaction phenomenon associated with the 2017 Pohang earthquake emphasizes that there is an urgent need of liquefaction potential mapping for the Pohang city and other areas with a similar geological setting. In areas underlain by extensive unconsolidated basin fill sediments—where the records of past earthquakes are exiguous or indistinct and there is poor implementation of building codes—future earthquakes of similar or larger magnitude as the Pohang earthquake are likely to occur again. Therefore, this represents a hazard that may cause significant societal and economic threats in the future. Full article