Threshold contour production of rainfall intensity that induces landslides in susceptible regions of northern Turkey

This study is aimed at developing rainfall intensity contour lines to illustrate upper and lower rainfall intensity limits that trigger landslides. The first phase began by developing precipitation empirical thresholds as follows: (i) the relationship between intensity (I) and accumulative rainfall (E) and date (I-date and E-date, respectively); (ii) the antecedent rainfall up to 3, 5, 10, 15, and 30 days prior to landslide occurrence; (iii) the relationship between the I and the duration (D); and (iv) the relationship between the cumulative rainfall event and the D (ID and ED, respectively). The data recorded by two rain gauges (Rize and Rize–Pazar) in the province of Rize in Northwest Turkey were used to analyze 24 previous landslide events during the period 1985 to 2006. The second phase began by developing surface interpolation maps after deriving threshold and assumption of normality, thereby producing contours of the rainfall intensity data of the rain gauges. All thresholds were verified by no-rainfall events, and the ID threshold successfully distinguished 97 and 95% (in Rize and Rize–Pazar, respectively) of false alarm days with the highest accuracy among the thresholds (antecedent rainfall days, I-date, and E-date thresholds). Thus, a descriptive analysis was conducted using the rainfall data of selected rain gauges to test the normality of rainfall between the selected rain gauges in the study area. The linear correlation coefficient was found at 0.75, and other tests verified the normality with optimistic results. A universal or collaborative kriging interpolation (UK) in a geographic information systems’ environment was used to estimate the spatial distribution of rainfall intensity because rainfall is controlled robustly by topography and other factors. This distribution corresponds to the landslide day event and hole-effect mathematical models that are used to forecast un-sampled values in the study area. Topographic elevation, slope, and normalized difference vegetation index (NDVI) maps were integrated into the landslide day events by using rainfall intensity. The UK map was validated by a cross-validation procedure using root-mean-square error. The highest interpolated surface map accuracy was reached without using an external drift (using rainfall intensity values only) and when using an NDVI as the external drift. The optimal model was selected by comparing the measured and predicted values at approximately 78% accuracy. Ultimately, the rainfall intensity contour lines were used to illustrate the upper and lower rainfall intensity limits that trigger landslides. This study used well-documented spatio-temporal data of landslide events and introduced a primary threshold that can be validated for any future event.

[1]  Isabel F. Trigo,et al.  Shallow and deep landslides induced by rainfall in the Lisbon region (Portugal): assessment of relationships with the North Atlantic Oscillation , 2005 .

[2]  M. Rossi,et al.  Rainfall thresholds for the initiation of landslides in central and southern Europe , 2007 .

[3]  Fausto Guzzetti,et al.  Rainfall induced landslides in December 2004 in south-western Umbria, central Italy: types, extent, damage and risk assessment , 2006 .

[4]  M. Rossi,et al.  The rainfall intensity–duration control of shallow landslides and debris flows: an update , 2008 .

[5]  S. Reis,et al.  A GIS-based comparative study of frequency ratio, analytical hierarchy process, bivariate statistics , 2011 .

[6]  Işık Yilmaz,et al.  The effect of the sampling strategies on the landslide susceptibility mapping by conditional probability and artificial neural networks , 2010 .

[7]  Biswajeet Pradhan,et al.  A novel ensemble decision tree-based CHi-squared Automatic Interaction Detection (CHAID) and multivariate logistic regression models in landslide susceptibility mapping , 2014, Landslides.

[8]  J. Carpenter,et al.  Practice of Epidemiology Comparison of Random Forest and Parametric Imputation Models for Imputing Missing Data Using MICE: A CALIBER Study , 2014 .

[9]  Penelope A. Hancock,et al.  Spatial interpolation of large climate data sets using bivariate thin plate smoothing splines , 2006, Environ. Model. Softw..

[10]  J. Zêzere,et al.  Rainfall thresholds for landsliding in Lisbon Area (Portugal) , 2002 .

[11]  Karl Ropkins,et al.  openair - An R package for air quality data analysis , 2012, Environ. Model. Softw..

[12]  P. Frattini,et al.  Approaches for defining thresholds and return periods for rainfall‐triggered shallow landslides , 2009 .

[13]  Fausto Guzzetti,et al.  Lithological and seasonal control on rainfall thresholds for the possible initiation of landslides in central Italy , 2012 .

[14]  Philip K. Hopke,et al.  Comparison between Conditional Probability Function and Nonparametric Regression for Fine Particle Source Directions , 2004 .

[15]  T. Giambelluca,et al.  Comparison of geostatistical approaches to spatially interpolate month‐year rainfall for the Hawaiian Islands , 2016 .

[16]  Chris M. Mannaerts,et al.  Influence of topography on rainfall variability in Santiago Island, Cape Verde , 2014 .

[17]  P. Reichenbach,et al.  Regional hydrological thresholds for landslides and floods in the Tiber River Basin (central Italy) , 1998 .

[18]  C. Bernhofer,et al.  Comparison of spatial interpolation methods for the estimation of precipitation distribution in Distrito Federal, Brazil , 2015, Theoretical and Applied Climatology.

[19]  Menno-Jan Kraak,et al.  Principles of geographic information systems : an introductory textbook , 2000 .

[20]  A. Yalçın GIS-based landslide susceptibility mapping using analytical hierarchy process and bivariate statistics in Ardesen (Turkey): Comparisons of results and confirmations , 2008 .

[21]  Jason W. Kean,et al.  Objective definition of rainfall intensity–duration thresholds for the initiation of post-fire debris flows in southern California , 2013, Landslides.

[22]  Richard M. Iverson,et al.  Landslide triggering by rain infiltration , 2000 .

[23]  F. Bulut,et al.  Landslide susceptibility mapping using GIS and digital photogrammetric techniques: a case study from Ardesen (NE-Turkey) , 2007 .

[24]  R. Blessing,et al.  An empirical correction for absorption anisotropy. , 1995, Acta crystallographica. Section A, Foundations of crystallography.

[25]  Selçuk Reis,et al.  Designing and Developing a Province-Based Spatial Database for the Analysis of Potential Environmental Issues in Trabzon, Turkey , 2008 .

[26]  V. De Martinis,et al.  The Sustainable City: Urban Regeneration and Sustainability , 2000 .

[27]  Samuele Segoni,et al.  Integration of rainfall thresholds and susceptibility maps in the Emilia Romagna (Italy) regional-scale landslide warning system , 2015, Landslides.

[28]  Samuele Segoni,et al.  Rainfall thresholds for the forecasting of landslide occurrence at regional scale , 2012, Landslides.

[29]  C. K. Singh,et al.  DEM-based delineation for improving geostatistical interpolation of rainfall in mountainous region of Central Himalayas, India , 2017, Theoretical and Applied Climatology.

[30]  S. L. Gariano,et al.  An algorithm for the objective reconstruction of rainfall events responsible for landslides , 2015, Landslides.

[31]  G. Wieczorek,et al.  LANDSLIDES: INVESTIGATION AND MITIGATION. CHAPTER 4 - LANDSLIDE TRIGGERING MECHANISMS , 1996 .

[32]  Clayton V. Deutsch,et al.  Geostatistical Reservoir Modeling , 2002 .

[33]  S. Oliveira,et al.  Rainfall thresholds for landslide activity in Portugal: a state of the art , 2015, Environmental Earth Sciences.

[34]  S. Campus Evaluation and Prevention of Natural Risks , 2009 .

[35]  Selçuk Reis,et al.  Analyzing Land Use/Land Cover Changes Using Remote Sensing and GIS in Rize, North-East Turkey , 2008, Sensors.

[36]  M. Martina,et al.  Probabilistic rainfall thresholds for landslide occurrence using a Bayesian approach , 2012 .

[37]  A. Pasuto,et al.  Rainfall as a trigger of shallow mass movements. A case study in the Dolomites, Italy , 1998 .

[38]  J. E. Gartner,et al.  Storm rainfall conditions for floods and debris flows from recently burned areas in southwestern Colorado and southern California , 2008 .

[39]  H. Saito,et al.  Relationship between the initiation of a shallow landslide and rainfall intensity—duration thresholds in Japan , 2010 .

[40]  G. Wahba Spline models for observational data , 1990 .

[41]  Laurent Delobbe,et al.  Evaluation of radar-gauge merging methods for quantitative precipitation estimates , 2009 .

[42]  Stef van Buuren,et al.  MICE: Multivariate Imputation by Chained Equations in R , 2011 .

[43]  S. Wang,et al.  Comparison of interpolation methods for estimating spatial distribution of precipitation in Ontario, Canada , 2014 .

[44]  R. Yoshinaka,et al.  SLOPE FAILURES CAUSED BY HEAVY RAINFALL IN JAPAN , 1974 .

[45]  A. H. Thiessen PRECIPITATION AVERAGES FOR LARGE AREAS , 1911 .

[46]  C. Schlüchter,et al.  Paleoglaciations in Anatolia: A Schematic Review and First Results , 2005 .

[47]  N. Caine,et al.  The Rainfall Intensity - Duration Control of Shallow Landslides and Debris Flows , 1980 .

[48]  A. R. Mahmud,et al.  Prediction of slope failures using bivariate statistical based index of entropy model , 2012, 2012 IEEE Colloquium on Humanities, Science and Engineering (CHUSER).

[49]  C. Shope,et al.  Modeling Spatiotemporal Precipitation: Effects of Density, Interpolation, and Land Use Distribution , 2015 .

[50]  Biswajeet Pradhan,et al.  Estimation of rainfall threshold and its use in landslide hazard mapping of Kuala Lumpur metropolitan and surrounding areas , 2015, Landslides.

[51]  Hasan Arman,et al.  Rainfall Thresholds for The Initiation of Landslides in Trabzon Province of Turkey , 2010 .

[52]  T. Chai,et al.  Root mean square error (RMSE) or mean absolute error (MAE)? – Arguments against avoiding RMSE in the literature , 2014 .

[53]  Grace Wahba,et al.  Spline Models for Observational Data , 1990 .

[54]  Yuk Feng Huang,et al.  Rainfall-induced landslides in Hulu Kelang area, Malaysia , 2013, Natural Hazards.

[55]  O. F. Althuwaynee,et al.  Ionospheric TEC from the Turkish Permanent GNSS Network (TPGN) and comparison with ARMA and IRI models , 2017 .

[56]  Timothy C. Coburn,et al.  Geostatistics for Natural Resources Evaluation , 2000, Technometrics.

[57]  Davide Tiranti,et al.  Estimation of rainfall thresholds triggering shallow landslides for an operational warning system implementation , 2010 .

[58]  K. Vijay Kumar,et al.  A Novel Ensemble Based Decision Tree Model For High Dimensional Biomedicine Data , 2016 .

[59]  Fausto Guzzetti,et al.  Rainfall thresholds for the possible occurrence of landslides in Italy , 2010 .

[60]  B. Pradhan,et al.  Regional prediction of landslide hazard using probability analysis of intense rainfall in the Hoa Binh province, Vietnam , 2013, Natural Hazards.

[61]  S. L. Gariano,et al.  Rainfall thresholds for shallow landslide occurrence in Calabria, southern Italy , 2013 .

[62]  Pietro Aleotti,et al.  A warning system for rainfall-induced shallow failures , 2004 .

[63]  Applicability of R statistics in analyzing landslides spatial patterns in Northern Turkey , 2017, 2017 2nd International Conference on Knowledge Engineering and Applications (ICKEA).