Ground Instability Research Articles

Correlation Between Geotechnical Indexes and Landslide Occurrence in Southwestern Cyprus Using GIS and Machine Learning

Landslides are significantly influenced by the properties of the geological materials. As such, effective landslide susceptibility and hazard assessment necessitates the use of carefully selected and well-organized spatial data on geology and ground characteristics. The present study explores the spatial correlation between landslides and geotechnical indexes pertinent to the problem of slope stability. For this purpose, a geodatabase containing geotechnical properties was created for the southwestern part of the island of Cyprus, an area noted for its frequent ground instability issues and the availability of a comprehensive landslide database. Then, statistical correlations are established between the recorded landslides (active and inactive) in the region and key geotechnical, geological and geomorphological factors. The analysis was initially performed using the Frequency Ratio method, followed by two advanced machine learning techniques, namely Random Forest and Shapley Additive Explanations. The analysis results reveal that for weak argillaceous geomaterials, the clay content and plasticity index constitute high importance variables, followed by the geomorphological factors such as the slope angle. In rocky formations with clear rock mass structure, as main factors emerge the Geological Strength Index and the uniaxial compressive strength. The strong correlation identified between geotechnical factors and the spatial distribution of landslides underscores the potential benefit of integrating these variables in landslide susceptibility and hazard analysis methodologies. This adjusted emphasis provides clearer insights into the relationship between geological properties and landslide occurrences, which is crucial for developing more accurate predictive models and effective mitigation strategies.

Permafrost instability negates the positive impact of warming temperatures on boreal radial growth

Climate warming can alleviate temperature and nutrient constraints on tree growth in boreal regions, potentially enhancing boreal productivity. However, in permafrost environments, warming also disrupts the physical foundation on which trees grow, leading to leaning trees or "drunken" forests. Tree leaning might reduce radial growth, undermining potential benefits of warming. Here, we found widespread radial growth reductions in southern latitude boreal forests since the 1980s. At mid latitudes, radial growth increased from ~1980 to ~2000 but showed recent signs of decline afterward. Increased growth was evident since the 1980 s at higher latitudes, where radial growth appears to be temperature limited. However, recent changes in permafrost stability, and the associated increased frequency of tree leaning events, emerged as a significant stressor, leading to reduced radial growth in boreal trees at the highest latitudes, where permafrost is extensive. We showed that trees growing in unstable permafrost sites allocated more nonstructural carbohydrate reserves to offset leaning which compromised radial growth and potential carbon uptake benefits of warming. This higher allocation of resources in drunken trees is needed to build the high-density reaction wood, rich in lignin, that is required to maintain a vertical position. With continued climate warming, we anticipate widespread reductions in radial growth in boreal forests, leading to lower carbon sequestration. These findings enhance our understanding of how climate warming and indirect effects, such as ground instability caused by warming permafrost, will affect boreal forest productivity in the future.

Natural Hazard Assessment in the Southeastern Margin of the Ría de Arosa (Pontevedra, Spain) Using GIS Techniques

The characterization of natural hazards in coastal environments is of great necessity, especially in the current context of global climate change and increasing population concentrations. This research focuses on a multi-hazard analysis of the main geotechnical, geomorphological, hydrological, and lithological risks in the southeastern margin of the Ría de Arosa using Geographic Information System techniques. The integration of geotechnical characterization maps and natural hazard maps has allowed for the identification of areas with a high susceptibility to natural disasters, which is crucial for territorial planning and management in the context of growing urban pressure and global climate change. The results indicate that poorly consolidated surface formations, especially in transitional areas such as dunes and marshes, are particularly vulnerable. Additionally, areas with higher lithological competence have been identified, where slope changes contribute to ground instability. This analysis provides valuable tools for decision-making and the implementation of risk management policies, promoting sustainable development, the protection of coastal ecosystems, and the prevention of risks from urban planning and civil engineering activities in the Ría de Arosa.

Numerical evaluation of ground source heat pumps in a thawing permafrost region

Permafrost degradation poses significant environmental and geological challenges in Arctic and subarctic regions, particularly in areas like Umiujaq, Canada. The warming climate leads to thawing permafrost, causing ground instability, disrupting hydrology, and impacting local built environment. This study evaluates the use of Ground Source Heat Pump (GSHP) operation for mitigating ground subsidence in a permafrost region using a two-dimensional Thermo-Hydro-Mechanical (THM) coupled finite element analysis considering the ground poro-elastic and poro-plastic responses. The research uses a single-well scenario to demonstrate the interactions among thermal, hydraulic, and mechanical processes. The impact of GSHP operation under different temperature management strategies, including a scenario with a constant GSHP temperature of −5 °C throughout the year is numerically investigated. Results indicate that GSHP operation exacerbates ground deformation near the borehole, particularly during winter months. However, maintaining GSHP operation throughout the entire year can mitigate extreme subsidence fluctuations, leading to a more stable subsurface environment. While GSHP systems provide effective thermal regulation, their operation can introduce mechanical stresses that potentially disturb the ground close to the borehole. Therefore, careful design, operation, and further research are essential to balance thermal benefits with ground stability in permafrost regions.

Assessing landslide susceptibility in northern Morocco: A geostatistical mapping approach in Al Hoceima-Ajdir

Landslides are a significant global risk due to their impact on nature, society, and the economy. The Al Hoceima-Ajdir region is highly susceptible to geohazards such as earthquakes and ground instabilities caused by local tectonic movements and specific climatic and geomorphological conditions. This study aims to assess landslide susceptibility in Al Hoceima-Ajdir using a digital mapping approach that integrates the geostatistical Frequency Ratio tool, geomatic tools, and field observations. The frequency ratio method, a widely used technique in geospatial analysis, evaluates the relative frequency of landslides based on the spatial distribution of conditioning factors such as slope, lithology and land use – land cover. This method calculates the ratio of the frequency of landslides within specific classes of these factors to their respective frequencies in the entire study area. This research introduces a novel susceptibility map highlighting major threats to several natural and touristic sites crucial to the region's socio-economic activities. The susceptibility map was developed by analyzing 40 historical landslides along various predisposing factors, including topographic, geo-hydrological, and anthropogenic parameters. The results identify areas prone to landslides and classify their susceptibility into five levels ranging from very high to very low, with each level indicating a different degree of risk and requiring a specific set of preventive measures. The reliability of this map was evaluated using the Receiver Operating Characteristics (ROC) method, demonstrating a high accuracy of 89%, further validated by field visits. This high-precision susceptibility map is significant for pre- and post-risk management, enhancing prevention and mitigation strategies. It also protects society, the economy, and the environment by providing a reliable regional planning and hazard management tool. The findings of this study offer a significant advancement in the methodology for landslide susceptibility assessment and present a model that can be adapted to other regions facing similar geohazards. The research underscores the importance of integrating geostatistical tools with field data to produce accurate and actionable susceptibility maps, benefitting the global scientific community by improving hazard prediction and management strategies.

Soil Instability and Building Damage in the Urban District, Morocco: Chefchaouen City as Case Study

Ground instability is one of the most significant issues that affect the social, economic, and environmental levels in the northern Rif area and many other parts of Morocco. To understand and analyze this natural hazard of terrain instability, Chefchaouen City is chosen as a case study. based on the inventory of ground movements carried out by the Geotechmed Laboratory to identify and characterize cases of instability in this city and their relationship with the geological structure and soil characteristics is the aim of this research. As well as, the inventory of constructions damaged by instability in all municipality districts. The results indicated that the soils in the administrative district were reliable. Several factors contributed to the 2016 landslide in the Ain Haouzi District, including high rainfall, water resurgence, and human activity. The damage inventory of 83 houses showed how the geographical distribution of the severely damaged homes was placed on top of the distribution of landslides in the surrounding areas. The management plan needs to adapt to the geomorphological characteristics and behavior of the subsoil. That will help to avoid any future material damage and provide better guidance for the planning and the city's expansion development, with its remarkable urban sprawl.

Are the Italian microzonation level 2 abacuses applicable in the Friuli Venezia Giulia (Italy) plain? Comparison between the national abacuses and the numerically simulated amplification factors and between the derived elastic response spectra

Seismic shaking of an area is strongly affected by the local geology. The so-called local site effects must be considered for the estimation of seismic effects on structures and urban planning. Thus, the seismic microzonation is the process aimed at identifying and mapping the subsoil local response in a given area, usually at urban/municipality scale and in terms of ground shaking parameters and susceptibility to ground instabilities. In Italy, for areas that can be schematised as a 1D subsoil model (e.g. alluvial plain), a simplified approach is proposed to quantify the seismic amplification (amplification factor, AF). This approach consists of tables of correspondences, called seismic abacuses, available for the whole national area as well as for some regional territories, and derived for simplified subsoil models. In this work, the results of the comparison between the AF values retrieved from national abacuses applied in the Friuli Venezia Giulia (Italy) plain municipalities and those from 1D numerical simulations are presented. In general, the abacuses underestimate the local seismic site effects a part for sites with a shallow bedrock. No correlations/trends were identified between the AF derived from abacuses and those from numerical simulations. Moreover, considering the elastic acceleration response spectra, it emerges that in the 49.5% of the FVG analysed sites the abacuses approach, even though it underestimates the real seismic response, is a more suitable approximation compared to the soil class simplified approach proposed by the Italian regulation. Finally, what emerges is that the limit of 30 m, as indicated in the Italian regulation, to consider a deep or shallow bedrock seams underestimated, and the AFs are not correlated with the seismic bedrock depth when it is higher than 100 m.

Spatial and temporal evolution of mining-induced sinkholes in the Gavorrano area (Tuscany, Italy): Insights from remote sensing and field data

The failure of underground cavities and the formation of sinkholes is a major risk at both active and dismissed mining areas. However, while ground deformation related to most active mines is constantly monitored today, the ground instabilities around historical and dismissed mining areas often remain disregarded. Here, we combined a multidisciplinary dataset of historical mining reports and maps, aerial images, field measurements and InSAR observations to reconstruct the spatial and temporal evolution of ground deformation around a sinkhole prone area in the mining district of Gavorrano (Tuscany, Italy). We identified the three sinkholes of Monte Calvo, Valsecchi, and Ravi, with the latter two never reported in literature. The three sinkholes have a spatial correlation with the underground cavities and galleries of major mines. InSAR observations between 20162022 also revealed that an area of ∼ 700 m × 400 m around the Monte Calvo sinkhole is still subsiding with rates of ∼5 mm/yr. Conversely, no clear evidence of deformation is observed at Valsecchi, Ravi, and the nearby city of Gavorrano. Our observations suggest that the sinkhole activity in Gavorrano has been induced by the mining activity in the area and that subsidence is still ongoing at Monte Calvo.

Biomorphological traits and leaf dry matter content are important to plant persistence in a highly unstable volcanic ground

The colonisation of newly formed territories by plants during primary succession is a crucial stage in the formation of ecosystems. Adaptations of species to withstand harsh environment allow them to survive and form pioneer communities. We analysed ten categorical and quantitative plant traits to examine their role in plant resistance to substrate instability in primary volcanic habitats. The research questions were: 1) which plant traits enable plants to persist under the most unstable ground conditions?; 2) what are the ecological strategies of species that grow under the condition of stress, such as low mineral nutrition, coupled with intense disturbance, such as unstable ground? The research was conducted on the Tolbachinsky Dol plateau in Kamchatka, Russia, on loose volcanic sediments of the 1975 eruption. At altitudes ranging from 700 m a.s.l. to 1000 m a.s.l., 40 sample plots were established along the gradients of altitude and degrees of ground instability. The species composition of vascular plants and the percentage cover of species, as well as other habitat characteristics, were assessed. Linear models were used to investigate the relationship between traits and ground instability. The principal component analysis was used to identify groups of characteristic traits. Leaf dry matter content (LDMC) had the largest explained variance. The LDMC community weighted means decreased with increasing disturbance (ground instability). A whole set of categorical traits, i.e. plant functional types, characterised the adaptation of plants to unstable ground. These characteristics include life-form (according to Serebryakov classification), density of shoots arrangement, morphology of underground organs, and, to a lesser extent, dispersal mode. The study analysed 45 species, and it was found that only four of them were adapted to perform under conditions of high substrate instability. The study showed that the mutual interaction of a suite of biomorphological traits and leaf density characteristics formed a syndrome of pioneer species adapted to unstable ground. The ecological strategies of the majority of species in the studied communities were characterised by the predominance of the stress-tolerant strategy, explained by the extremely poor mineral nutrition in the primary habitats. Under conditions of both low nutrient availability and high disturbance intensity, species with a mixed stress-tolerant-ruderal strategy are favoured. In the gradient of increasing disturbance, the role of the ruderal axis increased while stress tolerance decreased. However, the role of ecological strategy may be weakened by the presence of other important adaptive traits.

Mental illness in former residents affected by the extraction of rock salt from a mining company in Maceió, Alagoas State, Brazil

This study aimed to investigate the presence of mental illness in victims of soil instability in neighborhoods affected by rock salt extraction from a mining company located in the city of Maceió, Alagoas, Brazil. It is a quantitative, descriptive-analytical, and cross-sectional study. The sample was intentional and non-probabilistic and consisted of 158 participants, with a 0.50 "power of effect" (moderate size effect) and a 99.9% confidence level. A link was made available on social media with an invitation to participate in the online survey. Most participants were adult, female, black/mixed-race, with a high level of education, and reported a perceived worsening in monthly income, physical and mental health, and in the self-perception of their lives after the relocation of their homes. Most participants had depressive symptoms (87.34%), moderate and severe anxiety (55.7%), and positive screening for common mental disorders (77.22%). By analyzing the current perception of the changes in the lives of former residents before and after being relocated due to the ground instability, it was identified that, after relocation, those with greater psychiatric demand and suicidal ideation had a higher prevalence of depressive symptoms, symptoms of severe anxiety, and positive screening for common mental disorders. This research highlights the importance of elaborating policies, strategies and actions to prevent impacts on mental health, recovery, and health rehabilitation in scenarios of socio-environmental disasters. Implementing preventive measures and offering adequate psychological and social support can reduce the detrimental effects of such disasters on the affected population.

Ground investigation: a tool for landslide risk mitigation of infrastructures

Recent monitoring experiences for bridges in Central and South Italy have shown that it is not common for bridges to be directly damaged by landslides phenomena in their vicinities, either because their siting was very carefully selected at the origin or because appropriate protection measures have been taken to reduce geological risk to acceptable levels. In the first case, periodic visual inspections of the manufacts may suffice to monitor possible consequences of ground instability phenomena, but in the second case monitoring of the infrastructure and of the ground is the only possible tool to evaluate the existence of a landslide-structure interaction which can influence its performance and safety. Ground investigation is a fundamental activity for landslide monitoring and analysis, and this is specifically true when a landslide-infrastructure interaction problem needs to be investigated and assessed. By discussing the results from case studies recently inspected, this paper shows how the knowledge of the physical processes and the monitoring of selected parameters are the first irreplaceable tools to mitigate risks associated with instability phenomena interacting with sensitive infrastructures.

A novel model for multi-risk ranking of buildings at city level based on open data: the test site of Rome, Italy

In the context of population concentration in large cities, assessing the risks posed by geological hazards to enhance urban resilience is becoming increasingly important. This study introduces a robust and replicable procedure for assessing ground instability hazards and associated physical risks. Specifically, our comprehensive model integrates spatial hazard assessments, multi-satellite InSAR data, and physical features of the built environment to rank and prioritize assets facing multiple risks, with a focus on ground instabilities. The model generates risk scores based on hazard probability, potential damage, and displacement rates, aiding decision-makers in identifying high-risk buildings and implementing appropriate mitigation measures to reduce economic losses. The procedure was tested in Rome, Italy, where the analysis revealed that 60% of the examined buildings (90 × 103) are at risk of ground instability. Specifically, 33%, 22%, and 5% exhibit the highest multi-risk score for sinkholes, landslides, and subsidence, respectively. Landslide risk prevails among residential structures, while retail and office buildings face a higher risk of subsidence and sinkholes. Notably, our study identified a positive correlation between mitigation expenses and the multi-risk scores of nearby buildings, highlighting the practical implications of our findings for urban planning and risk management strategies.

InSAR time series and LSTM model to support early warning detection tools of ground instabilities: mining site case studies

Early alarm systems can activate vital precautions for saving lives and the economy threatened by natural hazards and human activities. Interferometric synthetic aperture radar (InSAR) products generate valuable ground motion data with high spatial and temporal resolutions. Integrating the InSAR products and forecasting models make possible to set up early alarm systems to monitor vulnerable areas. This study proposes a technical support to early warning detection tools of ground instabilities using machine learning and InSAR time series that is capable of forecasting regions affected by potential collapses. A long short-term memory (LSTM) model is tailored to predict ground movements in three forecast ranges (i.e., SAR observations): 3, 4, and 5 multistep. A contribution of the proposed strategy is utilizing adjacent time series to decrease the possibility of falsely detecting safe regions as significant movements. The proposed tool offers ground motion-based outcomes that can be interpreted and utilized by experts to activate early alarms to reduce the consequences of possible failures in vulnerable infrastructures, such as mining areas. Three case studies in Spain, Brazil, and Australia, where fatal incidents happened, are analyzed by the proposed early alert detector to illustrate the impact of chosen temporal and spatial ranges. Since most early alarm systems are site dependent, we propose a general tool to be interpreted by experts for activating reliable alarms. The results show that the proposed tool can identify potential regions before collapse in all case studies. In addition, the tool can suggest an optimum selection of InSAR temporal (i.e., number of images) and spatial (i.e., adjacent measurement points) combinations based on the available SAR images and the characteristics of the study area.

Effect of Particle Size Distributions (PSDs) on ground responses induced by tunnelling in dense coarse-grained soils: A DEM investigation

Tunnelling-induced ground responses in coarse-grained soils often exhibit complex behaviours such as soil arching and progressive ground instability. However, the effect of particle size distributions (PSDs) on ground responses in tunnelling remains largely unexplored. This study provides new insights into PSD’s effect on tunnelling ground response and the underlying mechanism based on granular skeleton quantification. The coarse-grained soils with 9 typical PSDs were modelled by the discrete element method (DEM) with varied granular skeleton structures. A representative volume element (RVE) assembly method, with a novel iteration-based specimen generation method, is developed to build ground models incorporating millions of soil grains, where an eccentric displacement boundary was embedded for reproducing the tunnel volume losses. The results demonstrated that the change in skeleton structures from fine-grains to coarse-grains dominated leads to the acceleration of soil arch emergence and the deceleration of soil arch collapse-rebuilding. With an increase in ded, the strengthened soil arch leads to a more evident deformation localization in ground, and results in more significant hysteresis of ground instability from the tunnel to the surface. The effect of PSDs on ground responses becomes prominent as the overburden thickness decreases and the tunnel volume loss increases, and a decrease in ded would increase the risk of overall failure in tunnelling.

Exploiting satellite SAR for archaeological prospection and heritage site protection

ABSTRACT Optical and Synthetic Aperture Radar (SAR) remote sensing has a long history of use and reached a good level of maturity in archaeological and cultural heritage applications, yet further advances are viable through the exploitation of novel sensor data and imaging modes, big data and high-performance computing, advanced and automated analysis methods. This paper showcases the main research avenues in this field, with a focus on archaeological prospection and heritage site protection. Six demonstration use-cases with a wealth of heritage asset types (e.g. excavated and still buried archaeological features, standing monuments, natural reserves, burial mounds, paleo-channels) and respective scientific research objectives are presented: the Ostia-Portus area and the wider Province of Rome (Italy), the city of Wuhan and the Jiuzhaigou National Park (China), and the Siberian “Valley of the Kings” (Russia). Input data encompass both archive and newly tasked medium to very high-resolution imagery acquired over the last decade from satellite (e.g. Copernicus Sentinels and ESA Third Party Missions) and aerial (e.g. Unmanned Aerial Vehicles, UAV) platforms, as well as field-based evidence and ground truth, auxiliary topographic data, Digital Elevation Models (DEM), and monitoring data from geodetic campaigns and networks. The novel results achieved for the use-cases contribute to the discussion on the advantages and limitations of optical and SAR-based archaeological and heritage applications aimed to detect buried and sub-surface archaeological assets across rural and semi-vegetated landscapes, identify threats to cultural heritage assets due to ground instability and urban development in large metropolises, and monitor post-disaster impacts in natural reserves.

InterpolatiON of InSAR Time series for the dEtection of ground deforMatiOn eVEnts (ONtheMOVE): application to slow-moving landslides

The aim of this work is to develop an innovative methodology to analyse the time series (TS) of interferometric satellite data. TS are important tools for the ground displacement monitoring, mostly in areas in which in situ instruments are scarce. The proposed methodology allows to classify the trend of TS in three classes (uncorrelated, linear, non-linear) and to obtain the parameters of non-linear time series to characterise the magnitude and timing of changes of ground instabilities. These parameters are the beginning and end of the non-linear deformation break(s), the length of the event(s) in days, and the quantification of the cumulative displacement in mm. The methodology was tested on two Sentinel-1 datasets (2014–2020) covering the Alpine and Apennine sectors of the Piemonte region, an area prone to slow-moving slope instabilities. The results were validated at the basin scale (Pellice-Chisone and Piota basin) and at a local scale (Brenvetto, Champlas du Col and Casaleggio Boiro landslides) comparing with in situ monitoring system measurements, possible triggering factors (rainfall, snow) and already-collected events of the territory. The good correlation of the results has proven that the methodology can be a useful tool to local and regional authorities for risk planning and management of the area, also in terms of near real-time monitoring of the territory both at local and regional scale.

Bayesian Monte Carlo Inversion of InSAR Time Series Deformation Induced by Wastewater Injection: A Case Study in West Texas

AbstractWastewater disposal can induce detectable surface uplift, which may cause ground instability and threaten infrastructure. The distributions of local hydro‐geomechanical parameters, especially Young's modulus and hydraulic conductivity, play an essential role in these geohazards. To constrain these parameters, we have inverted spatio‐temporal deformation measured by Interferometric Synthetic Aperture Radar (InSAR) and injection information using a Bayesian Monte Carlo approach with a poroelastic finite element model. Sentinel‐1A/B imagery from 2014 to 2020 is processed to track the spatio‐temporal deformation in Winkler county, West Texas, USA. The posterior distribution of subsurface effective volumes reveals under‐reported volumes in the well near the deformation center. In addition, the inversion results provide better constraints for the parameters than those solely obtained based on the cumulative spatial deformation or temporal development of the deformation center.

The mechanisms and drivers of lithification in slag‐dominated artificial ground

AbstractUnconsolidated artificial ground is an ever‐increasing feature on the Earth's surface but it poses various challenges such as pollutant release and ground instability. The process of lithification could be an important factor in changing the properties of artificial ground and ameliorating these challenges. In this study, a lithified deposit of a furnace slag associated with a former iron and steel works in Scotland was analysed to determine the mechanisms and drivers of lithification. Scanning electron microscope analysis showed that Ca leached from around the edges of clasts of slag through reaction of the chemically unstable slag with water from an adjacent water body. Dissolution of Ca (and OH−) from the slag caused the water in contact with the slag to become hyperalkaline, facilitating ingassing and hydroxylation of CO2 from the atmosphere (fingerprinted through carbon isotope analysis). Reaction of the dissolved Ca and CO2 led to precipitation of calcite. Scanning electron microscope analysis showed the calcite is distributed between slag clasts, forming rims around clasts and cementing clasts together into a solid rock‐like mass. Understanding the mechanisms and drivers of lithification in artificial ground will be important, given its widespread nature particularly in urban areas where artificial ground is the substrate of most development.

Appraisal of Ancient Quarries and WWII Air Raids as Factors of Subsidence in Rome: A Geomatic Approach

Ancient mining and quarrying activities left anthropogenic geomorphologies that have shaped the natural landscape and affected environmental equilibria. The artificial structures and their related effects on the surrounding environment are analyzed here to characterize the quarrying landscape in the southeast area of Rome in terms of its dimensions, typology, state of preservation and interface with the urban environment. The increased occurrence of sinkhole events in urban areas has already been scientifically correlated to ancient cavities under increasing urban pressure. In this scenario, additional interacting anthropogenic factors, such as the aerial bombardments perpetrated during the Second World War, are considered here. These three factors have been investigated by employing a combined geomatic methodology. Information on air raids has been organized in vector archives. A dataset of historical aerial photographs has been processed into Digital Surface Models and orthomosaics to reconstruct the quarry landscape and its evolution, identify typologies of exploitation and forms of collapse and corroborate the discussion concerning the induced historical and recent subsidence phenomena, comparing these outputs with photogrammetric products obtained from recent satellite data. Geological and urbanistic characterization of the study area allowed a better connection between these historical and environmental factors. In light of the information gathered so far, SAR interferometric products allowed a preliminary interpretation of ground instabilities surrounding historical quarries, air raids and recent subsidence events. Various sub-areas of the AOI where the presence of the considered factors also corresponds to areas in slight subsidence in the SAR velocity maps have been highlighted. Bivariate hotspot analysis allowed substantiating the hypothesis of a spatial correlation between these multiple aspects.

A ground instabilities inventory between Tetouan and Jebha (Morocco): mapping, description and analysis

A ground instabilities inventory between Tetouan and Jebha (Morocco): mapping, description and analysis