Despite the large number of recent advances and developments in landslide susceptibility mapping (LSM) there is still a lack of studies focusing on specific aspects of LSM model sensitivity. For ...example, the influence of factors such as the survey scale of the landslide conditioning variables (LCVs), the resolution of the mapping unit (MUR) and the optimal number and ranking of LCVs have never been investigated analytically, especially on large data sets. In this paper we attempt this experimentation concentrating on the impact of model tuning choice on the final result, rather than on the comparison of methodologies. To this end, we adopt a simple implementation of the random forest (RF), a machine learning technique, to produce an ensemble of landslide susceptibility maps for a set of different model settings, input data types and scales. Random forest is a combination of Bayesian trees that relates a set of predictors to the actual landslide occurrence. Being it a nonparametric model, it is possible to incorporate a range of numerical or categorical data layers and there is no need to select unimodal training data as for example in linear discriminant analysis. Many widely acknowledged landslide predisposing factors are taken into account as mainly related to the lithology, the land use, the geomorphology, the structural and anthropogenic constraints. In addition, for each factor we also include in the predictors set a measure of the standard deviation (for numerical variables) or the variety (for categorical ones) over the map unit. As in other systems, the use of RF enables one to estimate the relative importance of the single input parameters and to select the optimal configuration of the classification model. The model is initially applied using the complete set of input variables, then an iterative process is implemented and progressively smaller subsets of the parameter space are considered. The impact of scale and accuracy of input variables, as well as the effect of the random component of the RF model on the susceptibility results, are also examined. The model is tested in the Arno River basin (central Italy). We find that the dimension of parameter space, the mapping unit (scale) and the training process strongly influence the classification accuracy and the prediction process. This, in turn, implies that a careful sensitivity analysis making use of traditional and new tools should always be performed before producing final susceptibility maps at all levels and scales.
Background:
Thyroid nodules are a common finding in the general population, and their detection is increasing with the widespread use of ultrasound (US). Thyroid cancer is found in 5–15% of cases ...depending on sex, age, and exposure to other risk factors. Some US parameters have been associated with increased risk of malignancy. However, no characteristic seems sufficiently reliable in isolation to diagnose malignancy. The objective of this meta-analysis was to evaluate the diagnostic performance of US features for thyroid malignancy in patients with unselected thyroid nodules and nodules with indeterminate fine-needle aspiration (FNA) cytology.
Methods:
Electronic databases were reviewed for studies published prior to July 2012 that evaluated US features of thyroid nodules and reported postoperative histopathologic diagnosis. A manual search of references of review and key articles, and previous meta-analyses was also performed. A separate meta-analysis was performed including only nodules with indeterminate cytology. Analyzed features were solid structure, hypoechogenicity, irregular margins, absence of halo, microcalcifications, central vascularization, solitary nodule, heterogeneity, taller than wide shape, and absence of elasticity.
Results:
Fifty-two observational studies (12,786 nodules) were included. Nine studies included nodules with indeterminate cytology as a separate category, comprising 1851 nodules. In unselected nodules, all US features were significantly associated with malignancy with an odds ratio varying from 1.78 to 35.7, and microcalcifications, irregular margins, and a taller than wide shape had high specificities (Sp; 87.8%, 83.1%, 96.6%) and positive likelihood ratios (LHR; 3.26, 2.99, 8.07). Absence of elasticity was the single feature with the best diagnostic performance (sensitivity 87.9%, Sp 86.2%, and positive LHR 6.39). The presence of central vascularization was the most specific US feature in nodules with indeterminate cytology (Sp 96% and positive LHR 2.13).
Conclusions:
US features in isolation do not provide reliable information to select nodules that should have a FNA performed. A combination of US characteristics with higher likelihood ratios and consequently with higher post-test probabilities of malignancy—microcalcifications, or a taller than wide shape, or irregular margins, or absence of elasticity—will probably identify nodules with an increased risk for malignancy. Further studies are required to standardize elastography techniques and evaluate outcomes, especially in nodules with an indeterminate cytology.
This paper concerns a regional scale warning system for landslides that relies on a decisional algorithm based on the comparison between rainfall recordings and statistically defined thresholds. The ...latter were based on the total amount of rainfall, which was cumulated considering different time intervals: 1-, 2- and 3-day cumulates took into account the critical rainfall influencing shallow movements, whilst a variable time interval cumulate (up to 240 days) was used to consider the triggering of deep-seated landslides in low permeability terrains. A prototypal version of the model was initially set up to define statistical thresholds. Then, thresholds were calibrated using a database of past georegistered and dated landslides. A validation procedure showed that the calibration highly improves the results and therefore the model was integrated in the regional warning system of Emilia Romagna (Italy) for civil protection purposes. The proposed methodology could be easily implemented in other similar regions and countries where a sufficiently organised meteorological network is present.
Classification and regression problems are a central issue in geosciences. In this paper, we present Classification and Regression Treebagger (ClaReT), a tool for classification and regression based ...on the random forest (RF) technique. ClaReT is developed in Matlab and has a simple graphic user interface (GUI) that simplifies the model implementation process, allows the standardization of the method, and makes the classification and regression process reproducible. This tool performs automatically the feature selection based on a quantitative criterion and allows testing a large number of explanatory variables. First, it ranks and displays the parameter importance; then, it selects the optimal configuration of explanatory variables; finally, it performs the classification or regression for an entire dataset. It can also provide an evaluation of the results in terms of misclassification error or root mean squared error. We tested the applicability of ClaReT in two case studies. In the first one, we used ClaReT in classification mode to identify the better subset of landslide conditioning variables (LCVs) and to obtain a landslide susceptibility map (LSM) of the Arno river basin (Italy). In the second case study, we used ClaReT in regression mode to produce a soil thickness map of the Terzona catchment, a small sub-basin of the Arno river basin. In both cases, we performed a validation of the results and a comparison with other state-of-the-art techniques. We found that ClaReT produced better results, with a more straightforward and easy application and could be used as a valuable tool to assess the importance of the variables involved in the modeling.
In this paper, the updating of the landslide inventory of Tuscany region is presented. To achieve this goal, satellite SAR data processed with persistent scatter interferometry (PSI) technique have ...been used. The updating leads to a consistent reduction of unclassified landslides and to an increasing of active landslides. After the updating, we explored the characteristics of the new inventory, analysing landslide distribution and geomorphological features. Several maps have been elaborated, as sliding index or landslide density map; we also propose a density-area map to highlight areas with different landslide densities and sizes. A frequency-area analysis has been performed, highlighting a classical negative power-law distribution. We also explored landslide frequency for lithology, soil use and several morphological attributes (elevation, slope gradient, slope curvature), considering both all landslides and classified landslide types (flows, falls and slides).
Within the framework of FP7, an EU-funded SafeLand project, a questionnaire was prepared to collect information about the use of remote sensing for landslide study and to evaluate its actual ...application in landslide detection, mapping and monitoring. The questionnaire was designed using a Google form and was disseminated among end-users and researchers involved in landslide studies in Europe. In total, 49 answers from 17 different European countries were collected. The outcomes showed that landslide detection and mapping is mainly performed with aerial photos, often associated with optical and radar imagery. Concerning landslide monitoring, satellite radars prevail over the other types of data. Remote sensing is mainly used for detection/mapping and monitoring of slides, flows and lateral spreads with a preferably large scale of analysis (1:5000-1:25 000). All the compilers integrate remote sensing data with other thematic data, mainly geological maps, landslide inventory maps and DTMs and derived maps. According to the research and working experience of the compilers, remote sensing is generally considered to have a medium effectiveness/reliability for landslide studies. The results of the questionnaire can contribute to an overall sketch of the use of remote sensing in current landslide studies and show that remote sensing can be considered a powerful and well-established instrument for landslide mapping, monitoring and hazard analysis.
HIRESSS (HIgh REsolution Slope Stability Simulator) is a physically based distributed slope stability simulator for analyzing shallow landslide triggering conditions in real time and on large areas ...using parallel computational techniques. The physical model proposed is composed of two parts: hydrological and geotechnical. The hydrological model receives the rainfall data as dynamical input and provides the pressure head as perturbation to the geotechnical stability model that computes the factor of safety (FS) in probabilistic terms. The hydrological model is based on an analytical solution of an approximated form of the Richards equation under the wet condition hypothesis and it is introduced as a modeled form of hydraulic diffusivity to improve the hydrological response. The geotechnical stability model is based on an infinite slope model that takes into account the unsaturated soil condition. During the slope stability analysis the proposed model takes into account the increase in strength and cohesion due to matric suction in unsaturated soil, where the pressure head is negative. Moreover, the soil mass variation on partially saturated soil caused by water infiltration is modeled. The model is then inserted into a Monte Carlo simulation, to manage the typical uncertainty in the values of the input geotechnical and hydrological parameters, which is a common weak point of deterministic models. The Monte Carlo simulation manages a probability distribution of input parameters providing results in terms of slope failure probability. The developed software uses the computational power offered by multicore and multiprocessor hardware, from modern workstations to supercomputing facilities (HPC), to achieve the simulation in reasonable runtimes, compatible with civil protection real time monitoring. A first test of HIRESSS in three different areas is presented to evaluate the reliability of the results and the runtime performance on large areas.
We propose an original approach to develop rainfall thresholds to be used in civil protection warning systems for the occurrence of landslides at regional scale (i.e. tens of thousands of ...kilometres), and we apply it to Tuscany, Italy (23 000 km2). Purpose-developed software is used to define statistical intensity-duration rainfall thresholds by means of an automated and standardized analysis of rainfall data. The automation and standardization of the analysis brings several advantages that in turn have a positive impact on the applicability of the thresholds to operational warning systems. Moreover, the possibility of defining a threshold in very short times compared to traditional analyses allowed us to subdivide the study area into several alert zones to be analysed independently, with the aim of setting up a specific threshold for each of them. As a consequence, a mosaic of several local rainfall thresholds is set up in place of a single regional threshold. Even if pertaining to the same region, the local thresholds vary substantially and can have very different equations. We subsequently analysed how the physical features of the test area influence the parameters and the equations of the local thresholds, and found that some threshold parameters can be put in relation with the prevailing lithology. In addition, we investigated the possible relations between effectiveness of the threshold and number of landslides used for the calibration. A validation procedure and a quantitative comparison with some literature thresholds showed that the performance of a threshold can be increased if the areal extent of its test area is reduced, as long as a statistically significant landslide sample is present. In particular, we demonstrated that the effectiveness of a warning system can be significantly enhanced if a mosaic of site-specific thresholds is used instead of a single regional threshold.
In this paper, we present preliminary results of the IPL project No. 198 “Multi-scale rainfall triggering models for Early Warning of Landslides (MUSE).” In particular, we perform an assessment of ...the geotechnical and hydrological parameters affecting the occurrence of landslides. The aim of this study is to improve the reliability of a physically based model high resolution slope stability simulator (HIRESSS) for the forecasting of shallow landslides. The model and the soil characterization have been tested in Northern Tuscany (Italy), along the Apennine chain, an area that is historically affected by shallow landslides. In this area, the main geotechnical and hydrological parameters controlling the shear strength and permeability of soils have been determined by in situ measurements integrated by laboratory analyses. Soil properties have been statistically characterized to provide more refined input data for the slope stability model. Finally, we have tested the ability of the model to predict the occurrence of shallow landslides in response to an intense meteoric precipitation.
We set up an early warning system for rainfall-induced landslides in Tuscany (23 000 km2). The system is based on a set of state-of-the-art intensity–duration rainfall thresholds (Segoni et al., ...2014b) and makes use of LAMI (Limited Area Model Italy) rainfall forecasts and real-time rainfall data provided by an automated network of more than 300 rain gauges. The system was implemented in a WebGIS to ease the operational use in civil protection procedures: it is simple and intuitive to consult, and it provides different outputs. When switching among different views, the system is able to focus both on monitoring of real-time data and on forecasting at different lead times up to 48 h. Moreover, the system can switch between a basic data view where a synoptic scenario of the hazard can be shown all over the region and a more in-depth view were the rainfall path of rain gauges can be displayed and constantly compared with rainfall thresholds. To better account for the variability of the geomorphological and meteorological settings encountered in Tuscany, the region is subdivided into 25 alert zones, each provided with a specific threshold. The warning system reflects this subdivision: using a network of more than 300 rain gauges, it allows for the monitoring of each alert zone separately so that warnings can be issued independently. An important feature of the warning system is that the visualization of the thresholds in the WebGIS interface may vary in time depending on when the starting time of the rainfall event is set. The starting time of the rainfall event is considered as a variable by the early warning system: whenever new rainfall data are available, a recursive algorithm identifies the starting time for which the rainfall path is closest to or overcomes the threshold. This is considered the most hazardous condition, and it is displayed by the WebGIS interface. The early warning system is used to forecast and monitor the landslide hazard in the whole region, providing specific alert levels for 25 distinct alert zones. In addition, the system can be used to gather, analyze, display, explore, interpret and store rainfall data, thus representing a potential support to both decision makers and scientists.