Modelling of Cascading Hazards in Naturally- Dammed Lakes
Total Page:16
File Type:pdf, Size:1020Kb
MODELLING OF CASCADING HAZARDS IN NATURALLY- DAMMED LAKES ANISH RATNA SHAKYA June, 2020 SUPERVISORS: Prof. Dr, C. V., Westen Ir, G. N., Parodi MODELLING OF CASCADING HAZARDS IN NATURALLY- DAMMED LAKES ANISH RATNA SHAKYA Enschede, The Netherlands, June, 2020 Thesis submitted to the Faculty of Geo-Information Science and Earth Observation of the University of Twente in partial fulfilment of the requirements for the degree of Master of Science in Geo-information Science and Earth Observation. Specialization: Natural Hazards and Disaster Risk Reduction SUPERVISORS: Prof. Dr, C. V., Westen Ir, G. N., Parodi THESIS ASSESSMENT BOARD: Prof. Dr, N., Kerle (Chair) Dr, M., Mergili (External Examiner, BOKU University, Austria) DISCLAIMER This document describes work undertaken as part of a programme of study at the Faculty of Geo-Information Science and Earth Observation of the University of Twente. All views and opinions expressed therein remain the sole responsibility of the author, and do not necessarily represent those of the Faculty. ABSTRACT Naturally-dammed lake outburst floods are sudden, catastrophic events releasing high volume of water and sediments with high peak discharges in the mountainous region. The lake outburst results from the interaction of the underlying series of natural events affecting one upon another, also known as cascading hazards. Insights of the cascading hazards can be understood by modelling them in an integrative way. Among several processes involved in the lake outburst flood, this study focusses on integrated modelling of three major cascading processes: slope stability, impulse wave and flood. The main objective of this research is to evaluate the outburst flood hazard modelling of naturally-dammed lakes by integrated process chain modelling under different scenarios which aim to support the disaster risk mitigation planning. For this, Lake Sarez and Lake Rivakkul were chosen as two case study areas from the mountainous region of Tajikistan. Three scenarios were designed for Lake Sarez and two scenarios were designed for Lake Rivakkul. In Scenario 1 of Lake Sarez, the volume and depth of potential rockslide were predicted from a physically- based slope stability model, Scoops3D using the geotechnical parameters assigned from the past literature. The outputs were then coupled with the Heller-Hager model which estimated the impulse wave characteristics. The optimistic and pessimistic range of parameter values from past studies were chosen in Scenario 2 and Scenario 3 respectively. These parameter values were assigned to Heller-Hager model to estimate the wave run-up heights, overtopping volume and duration. Simple hydrographs were constructed from the outputs of the Heller-Hager model which could be coupled with the flood model. For the Lake Rivakkul, an ideal rotational slip was assumed for the potential rockslide in Scenario 1 from which the volume and depth were predicted. Heller-Hager was used to estimate the wave run-up heights. In Scenario 2, the dam impounding Lake Rivakkul was assumed to breach and completely empty the lake. The peak discharge and time to peak were estimated from the empirical equations to construct the input hydrograph for the two-dimensional flood modelling in the HEC-RAS. Terrain data was extracted from SRTM DEM (30 m) with the computational mesh size of 50 m. A full momentum scheme with the timestep of 5 seconds and a uniform roughness coefficient of 0.06 was used. The resulting flood was routed at 24 km and 60 km downstream of the lake outburst till Rivak village and Khorugh respectively. Scenario 1 and Scenario 3 generated the run-up heights that overtopped the Usoi dam impounding Lake Sarez in the expense of several limitations. However, the Usoi dam was considered stable by numerous studies so downstream flood modelling was not done but the hydrograph was constructed from the overtopping volume and duration demonstrating its plausibility of coupling with the flood model. Scenario 1 of Lake Rivakkul did not generate enough run-up height to overtop the dam therefore, it was not able to couple with the flood model. In Scenario 2 of Lake Rivakkul, the peak flood reached at the Rivak village and Khorugh at 2.25 hours and 4.45 hours respectively. The average inundation depths at these places reached about 3.6 m and 4 m accordingly. The results of the chain modelling were represented only as the first-pass hazard assessment. Uncertainties and limitations of each process models along with the challenges of modelling the lake outburst hazard chain were highlighted. Modelling of the lake outburst hazard cascade will inherit some uncertainties and quantification of uncertainties are recommended for a range of acceptable results. To overcome these challenges in the data-scarce mountainous environment, recommendations were made regarding the selection of advanced numerical models complemented by field investigations and the future use of physically-based multi-hazard models which can analyse the entire event within a single run. Keywords: Lake outburst, modelling, mountainous, rockslide, impulse waves, flood i ACKNOWLEDGEMENTS “Challenges are what make life interesting. And overcoming them is what makes life meaningful”. - Buddha First of all, I am grateful to my parents and brother. Their relentless love and support drove me to this achievement. I would like to extend my sincere gratitude to my first supervisor, Professor Dr. Cees van Westen for providing me with all support and knowledge, since the phase of choosing the research topic to the internship and to polishing the final product of this thesis. I am also extremely thankful to my second supervisor, Ir. Gabriel Parodi for his crucial guidance and insights. It is with their mentoring and constructive feedback that I could complete this work. I would like to acknowledge the Multi-hazard risk assessment project at district level in Tajikistan for providing me with essential thematic data for the study. I am grateful to all the teachers of NHR domain from whom I gained indispensable skills and knowledge to excel in my career. My gratitude goes to ITC Student Affairs for encouraging us to go through the tough times. To all my international friends from NHR and other domains as well, I am thankful for all the memories that we have shared. We faced a lot of challenges and most importantly, we overcame together and shared our cultures. Finally, I would like to thank my Nepalese friends in ITC, with whom I got to share my experiences and got to enjoy some quality time together. Their company made me feel homely. ii TABLE OF CONTENTS 1. INTRODUCTION .............................................................................................................................................. 7 1.1. Background ...................................................................................................................................................................7 1.2. Literature review ..........................................................................................................................................................8 1.3. Problem statement and scientific significance ..................................................................................................... 12 1.4. Research objectives and research questions ........................................................................................................ 13 1.5. Organization of thesis and workflow.................................................................................................................... 13 2. STUDY AREA AND DATA COLLECTION ........................................................................................... 15 2.1. Case study area 1: Lake Sarez ................................................................................................................................. 16 2.2. Case study area 2: Lake Rivakkul ........................................................................................................................... 20 3. DESCRIPTION OF THE MODELS USED ............................................................................................. 23 3.1. Slope stability model: Scoops3D ........................................................................................................................... 24 3.2. Impulse wave generation, propagation and run-up model: Heller-Hager model ......................................... 25 3.3. Flood routing model: HEC-RAS 5.0.7 ................................................................................................................. 27 3.4. Coupling of models .................................................................................................................................................. 28 4. HAZARD MODELLING OF LAKE SAREZ .......................................................................................... 30 4.1. Scenario design ......................................................................................................................................................... 30 4.2. Input data preparation ............................................................................................................................................. 31 4.3. Results and discussion ............................................................................................................................................. 38 5. HAZARD MODELLING OF LAKE RIVAKKUL ................................................................................