Stream Length-gradient Hotspot and Cluster Analysis (SL-HCA) to fine-tune the detection and interpretation of knickzones on longitudinal profiles

Abstract This work presents a new approach for fine-tuning the analysis of stream longitudinal profiles. We show that applying Hotspot and Cluster Analysis (HCA), based on the Getis-Ord Gi* statistic, to the stream length-gradient (SL) index improves the visualization of anomalous values, assisting in the identification of tectonic structures and large landslides. High positive Gi* values indicate the clustering of SL anomalies (hotspots), and mirror the occurrence of knickzones on the stream long-profiles. We applied this methodology to a mountainous sector of the eastern Emilia-Romagna region, in northern Italy. Remote sensing and field surveys conducted on hotspot sites indicate that large landslides are the main process associated to over-steepened long-profile segments along streams connected to the valley slopes. Along-stream changes in bedrock resistance accounts for the main anomalies within sectors where hillslopes and valley floors are disconnected. We demonstrate that specific relationships between geometry and intensity of SL hotspots are indicative of the process responsible for the knickzone formation and, in particular that tectonic structures generally provide the longest and highest anomalies. The results of this work suggest that SL-HCA maps are more advantageous for detecting and interpreting knickzones compared with traditional SL maps, since: i) they need less input data to be computed, thus making them useful to investigate regions poorly covered by detailed geological data and/or where field surveys are difficult to be carried out and ii) the hotspot geometry can help discriminate the knickzones attributable to gravitational mass movements from litho-structural ones.

[1]  Nicola Casagli,et al.  Radiocarbon Data on Lateglacial and Holocene Landslides in the Northern Apennines , 2004 .

[2]  A. Getis The Analysis of Spatial Association by Use of Distance Statistics , 2010 .

[3]  Peter A. Rogerson,et al.  Maximum Getis–Ord Statistic Adjusted for Spatially Autocorrelated Data , 2015 .

[4]  Francesco Troiani,et al.  Types and development of stream terraces in the Marche Apennines (central Italy): a review and remarks on recent appraisals , 2012 .

[5]  Earthquake-induced ground failures in Italy from a reviewed database , 2013 .

[6]  Gustavo Lopes Queiroz,et al.  Knickpoint finder: A software tool that improves neotectonic analysis , 2015, Comput. Geosci..

[7]  Alan D. Howard,et al.  Long Profile Development of Bedrock Channels: Interaction of Weathering, Mass Wasting, Bed Erosion, and Sediment Transport , 2013 .

[8]  C. Tellini,et al.  Tectonic control on the development and distribution of large landslides in the Northern Apennines (Italy) , 2016 .

[9]  J. Azañón,et al.  Geomorphic analysis of the Sierra Cabrera, an active pop-up in the constrictional domain of conjugate strike-slip faults: The Palomares and Polopos fault zones (eastern Betics, SE Spain) , 2012 .

[10]  A. Cristea Assessment of recent tectonic evolution and geomorphic response in SE Carpathians (Romania) using hypsometric analysis , 2014 .

[11]  Y. Chen,et al.  Along-strike variations of morphotectonic features in the Western Foothills of Taiwan: tectonic implications based on stream-gradient and hypsometric analysis , 2002 .

[12]  T. Hoey,et al.  Scale Dependence of Lithological Control on Topography: Bedrock Channel Geometry and Catchment Morphometry in Western Scotland , 2010, The Journal of Geology.

[13]  L. Nibbi,et al.  Morphotectonic characterization of the quaternary intermontane basins of the Umbria-Marche Apennines (Italy) , 2014, Rendiconti Lincei.

[14]  V. Gornitz,et al.  River profiles along the Himalayan arc as indicators of active tectonics , 1983 .

[15]  José Vicente Pérez-Peña,et al.  Differentiating geology and tectonics using a spatial autocorrelation technique for the hypsometric integral , 2009 .

[16]  G. Pini Tectonosomes and Olistostromes in the Argille Scagliose of the Northern Apennines, Italy , 1999 .

[17]  D. Granger,et al.  Quantifying rock uplift rates using channel steepness and cosmogenic nuclide–determined erosion rates: Examples from northern and southern Italy , 2010 .

[18]  M. Rossi,et al.  Analysis of historical landslide time series in the Emilia‐Romagna region, northern Italy , 2010 .

[19]  M. Font,et al.  DEM and GIS analysis of the stream gradient index to evaluate effects of tectonics: The Normandy intraplate area (NW France) , 2010 .

[20]  Y. Hayakawa,et al.  GIS analysis of fluvial knickzone distribution in Japanese mountain watersheds , 2009 .

[21]  F. Pazzaglia,et al.  Knickzone propagation in the Black Hills and northern High Plains: A different perspective on the late Cenozoic exhumation of the Laramide Rocky Mountains , 2001 .

[22]  Y. Gunnell,et al.  An automated method for producing synoptic regional maps of river gradient variation: Procedure, accuracy tests, and comparison with other knickpoint mapping methods , 2011 .

[23]  Harvey M. Kelsey,et al.  Knickpoint and knickzone formation and propagation, South Fork Eel River, northern California , 2012 .

[24]  J. Keith Ord,et al.  Spatial Processes Models and Applications , 1981 .

[25]  J. T. Hack,et al.  Stream-profile analysis and stream-gradient index , 1973 .

[26]  Marco Pizziolo,et al.  Landslides in Emilia-Romagna region (Italy): strategies for hazard assessment and risk management , 2005 .

[27]  M. Pizziolo,et al.  The Susceptibility Map for Landslides with Shallow Initiation in the Emilia Romagna Region (Italy) , 2013 .

[28]  Alexander C. Whittaker,et al.  Fluvial archives, a valuable record of vertical crustal deformation , 2017 .

[29]  C. Esposito,et al.  Mutual interactions between slope-scale gravitational processes and morpho-structural evolution of central Apennines (Italy): review of some selected case histories , 2014, Rendiconti Lincei.

[30]  Veronica Tofani,et al.  Persistent Scatterers Interferometry Hotspot and Cluster Analysis (PSI-HCA) for detection of extremely slow-moving landslides , 2012 .

[31]  Massimiliano Stucchi,et al.  CPTI11, the 2011 version of the Parametric Catalogue of Italian Earthquakes , 2011 .

[32]  Salvatore Martino,et al.  Morpho-structural evolution of the valley-slope systems and related implications on slope-scale gravitational processes: New results from the Mt. Genzana case history (Central Apennines, Italy) , 2017 .

[33]  Giovanni Sarti,et al.  The Romagna Apennines, Italy: an eroded duplex , 2001 .

[34]  Aaron J. Martin,et al.  Correlations of fluvial knickzones with landslide dams, lithologic contacts, and faults in the southwestern Annapurna Range, central Nepalese Himalaya , 2012 .

[35]  J. Ord,et al.  Local Spatial Autocorrelation Statistics: Distributional Issues and an Application , 2010 .

[36]  O. Korup Rock-slope failure and the river long profile , 2006 .

[37]  José Vicente Pérez-Peña,et al.  Spatial analysis of stream power using GIS: SLk anomaly maps , 2009 .

[38]  N. Cressie The origins of kriging , 1990 .

[39]  P. Bishop,et al.  Regional analysis of bedrock stream long profiles: evaluation of Hack's SL form, and formulation and assessment of an alternative (the DS form) , 2007 .

[40]  L. Royden,et al.  An integral approach to bedrock river profile analysis , 2013 .

[41]  K. Whipple,et al.  Knickpoint initiation and distribution within fluvial networks: 236 waterfalls in the Waipaoa River, North Island, New Zealand , 2006 .

[42]  F. Troiani,et al.  The use of the Stream Length–Gradient index in morphotectonic analysis of small catchments: A case study from Central Italy , 2008 .

[43]  B. Silverman Density estimation for statistics and data analysis , 1986 .

[44]  Timothy O. Randhir,et al.  Ecohydrologic modelling of water resources and land use for watershed conservation , 2016, Ecol. Informatics.

[45]  K. Wegmann,et al.  Late Quaternary fluvial terraces of the Romagna and Marche Apennines, Italy: Climatic, lithologic, and tectonic controls on terrace genesis in an active orogen , 2009 .

[46]  Jorge Pedro Galve,et al.  Spatial analysis of stream length-gradient (SL) index for detecting hillslope processes: A case of the Gállego River headwaters (Central Pyrenees, Spain) , 2014 .