MDTanaliza: understanding digital elevation models when facing gravity-driven flows in a hazard assessment context

The Digital Elevation Model (DEM) is one of the most important inputs when working on computer simulations to define areas potentially affected by gravity-driven flow hazards. Sometimes there is only one DEM available in the study area, although nowadays the availability of several different DEMs of varying resolutions and precision is more frequent, the characteristics of which are not necessarily familiar for the users. This availability may increase during emergency events in high-risk areas where updated and more precise DEMs may be calculated every few months using less expensive and more advanced methods. In order to design gravity-driven flow hazard maps, it is critical to understand the changes that may be taking place in computer simulations that use various DEMs to get a general understanding, for example, of the drainage network representation or how resolution changes may make the flow paths differ and how this may affect the interpretation of existing hazard maps that may be enhanced with the new DEM and cartography-base maps. To support researchers in these tasks, a simple and flexible ANSI C software tool called MDTanaliza has been developed. The aim of this tool is to analyze all possible gravity-driven flow paths, recreating former or new morphologies and calculating various different basic morphometric parameters in order to understand and compare the DEMs available in a study area. We have also developed some innovative approaches to improve the gravity-driven flow path analysis: the Filling Depression function to work with non-depressionless DEM and to reduce the time computing; the Pathway Force Interaction function that allows the previous gravity-driven flow paths to be considered in the next path calculation; and the Restricted Multiflow function to control the dispersion in the Multi-Directional flow model. As an example, the summit area of the Cotopaxi volcano has been used.

[1]  D. Montgomery,et al.  Digital elevation model grid size, landscape representation, and hydrologic simulations , 1994 .

[2]  K. Beven,et al.  THE PREDICTION OF HILLSLOPE FLOW PATHS FOR DISTRIBUTED HYDROLOGICAL MODELLING USING DIGITAL TERRAIN MODELS , 1991 .

[3]  K. Jones.  A comparison of algorithms used to compute hill slope as a property of the DEM , 1998 .

[4]  Berthold K. P. Horn,et al.  Hill shading and the reflectance map , 1981, Proceedings of the IEEE.

[5]  M. D. Fleming,et al.  Machine Processing of Landsat MSS Data and DMA Topographic Data for Forest Cover Type Mapping , 1979 .

[6]  A. Harris,et al.  FLOWGO: a kinematic thermo-rheological model for lava flowing in a channel , 2001 .

[7]  J. Vallance,et al.  OBJECTIVE DELINEATION OF LAHAR-INUNDATION HAZARD ZONES , 1998 .

[8]  Santiago Beguería,et al.  Validation and Evaluation of Predictive Models in Hazard Assessment and Risk Management , 2006 .

[9]  P. Burrough,et al.  Principles of geographical information systems , 1998 .

[10]  P. Mothes,et al.  Character, Stratigraphy and Magnitude of Historical Lahars of Cotopaxi Volcano (Ecuador) , 2004 .

[11]  P. Pilesjö,et al.  Estimating slope from raster data: a test of eight different algorithms in flat, undulating and steep terrain , 2011 .

[12]  Robert D’Ercole,et al.  Percepción del riesgo asociado al volcán Cotopaxi y vulnerabilidad en el Valle de Los Chillos (Ecuador) , 2009 .

[13]  P. Mothes.  Lahars of Cotopaxi Volcano, Ecuador: hazard and risk evaluation , 1992 .

[14]  Friedrich Barnikel.  The value of historical documents for hazard zone mapping , 2004 .

[15]  P. Mothes,et al.  The 1877 lahar deposits on the eastern flank of Cotopaxi volcano , 2007 .

[16]  Vern Manville,et al.  The sensitivity of a volcanic flow model to digital elevation model accuracy: experiments with digitised map contours and interferometric SAR at Ruapehu and Taranaki volcanoes, New Zealand , 2003 .

[17]  L. Blair.  A prediction. , 1995, Hospitals & health networks.

[18]  P. Julien,et al.  Two‐Dimensional Water Flood and Mudflow Simulation , 1993 .

[19]  Alicia García,et al.  Assessment and Modelling of Lava Flow Hazard on Lanzarote (Canary Islands) , 2001 .

[20]  S. Schilling,et al.  Laharz_py: GIS tools for automated mapping of lahar inundation hazard zones , 2014 .

[21]  Carl Gerstenecker,et al.  Validation of Digital Elevation Models around Merapi Volcano, Java, Indonesia , 2005 .

[22]  M. Smyth.  Movement and emplacement mechanisms of the Rio Pita volcanic debris avalanche, and its role in the evolution of Cotopaxi volcano , 1991 .

[23]  S. K. Jenson,et al.  Extracting topographic structure from digital elevation data for geographic information-system analysis , 1988 .

[24]  Qiming Zhou,et al.  Error Analysis on Grid-Based Slope and Aspect Algorithms , 2004 .

[25]  M. Berti,et al.  DFLOWZ: A free program to evaluate the area potentially inundated by a debris flow , 2014, Comput. Geosci..

[26]  L. Wang,et al.  An efficient method for identifying and filling surface depressions in digital elevation models for hydrologic analysis and modelling , 2006, Int. J. Geogr. Inf. Sci..

[27]  John F. O'Callaghan,et al.  The extraction of drainage networks from digital elevation data , 1984, Comput. Vis. Graph. Image Process..

[28]  Joseph. Wood,et al.  The geomorphological characterisation of Digital Elevation Models , 1996 .

[29]  R. Cioni,et al.  Evidence for lahar-triggering mechanisms in complex stratigraphic sequences: the post-twelfth century eruptive activity of Cotopaxi Volcano, Ecuador , 2013, Bulletin of Volcanology.

[30]  David John Unwin,et al.  Introductory Spatial Analysis , 1982 .

[31]  Nathan J. Wood,et al.  Reducing risk from lahar hazards: concepts, case studies, and roles for scientists , 2014, Journal of Applied Volcanology.

[32]  S. Schilling.  Copies of this report can be purchased from: , 1996 .

[33]  D. Tarboton.  A new method for the determination of flow directions and upslope areas in grid digital elevation models , 1997 .

[34]  R. Cioni,et al.  Lahar hazard assessment in the southern drainage system of Cotopaxi volcano, Ecuador: Results from multiscale lahar simulations , 2014 .

[35]  José Manuel Marrero,et al.  Virtual tools for volcanic crisis management, and evacuation decision support: applications to El Chichón volcano (Chiapas, México) , 2013, Natural Hazards.

[36]  Ramón Ortiz,et al.  Automatic GIS-based system for volcanic hazard assessment , 2007 .

[37]  Wayne D. Iverson,et al.  VIEWIT: computation of seen areas, slope, and aspect for land-use planning , 1975 .

[38]  M. Ruiz,et al.  Relationship between volcanic ash fallouts and seismic tremor: quantitative assessment of the 2015 eruptive period at Cotopaxi volcano, Ecuador , 2016, Bulletin of Volcanology.

[39]  E. Jordan,et al.  Estimation by photogrammetry of the glacier recession on the Cotopaxi Volcano (Ecuador) between 1956 and 1997 / Estimation par photogrammétrie de la récession glaciaire sur le Volcan Cotopaxi (Equateur) entre 1956 et 1997 , 2005 .

[40]  P. Mothes,et al.  The rhyolitic–andesitic eruptive history of Cotopaxi volcano, Ecuador , 2008 .

[41]  Anselm Smolka,et al.  Natural disasters and the challenge of extreme events: risk management from an insurance perspective , 2006, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.