Hydro-morphodynamic modelling of a volcano-induced sediment-laden outburst flood at Sólheimajökull, Iceland

Outburst floods are one of the most catastrophic natural hazards for populations and infrastructure. Such high-magnitude sudden onset floods generally comprise of an advancing intense kinematic water wave that can induce considerable sediment transport. The exploration and investigation of sediment-laden outburst floods cannot be limited solely to water flow but must also include the flood-induced sediment transport. Understanding the complex flow-bed interaction process in large (field) scale outburst floods is still limited, not least due to a lack of well-constrained field data, but also because consensus on appropriate modelling schemes has yet to be decided. In recent years, attention has focussed on the numerical models capable of describing the process of erosion, transport and deposition in such flows and they are now at a point at which they provide useful quantitative data. Although the "exact" measure of bed change is still unattainable the numerical models enhance and improve insights into large outburst flood events. In this study, a volcano-induced jökulhlaup or glacial outburst flood (GLOF) at Sólheimajökull, Iceland is reproduced by novel 2D hydro-morphodynamic model that considers both bedload and suspended load based on shallow water theory. The simulation of sediment-laden outburst flood is shown to perform well, with further insights into the flow-bed interaction behaviour obtained from the modelling output. These results are beneficial to flood risk management and hazard prevention and mitigation. In summary, the modelling outputs show that (1) the quantity of bed erosion and deposition are sensitive to the sediment gain size, yet, the influences are not so significant when considering flow discharge; (2) finer resolution of topography increases the computational time significantly yet the results are not affected correspondingly; (3) the bed changes simulated by the present model achieves reasonably good agreement with those by the commercial Delft3D; (4) the flood is accelerated by about 30% due to the incorporation of sediment transport; (5) the rapid sediment-laden outburst flood causes a rapid morphological change and considerable amount of erosion and deposition, and the total erosion and deposition volumes increase simultaneously and tend to an approximate constant value; (6) and the peak erosion rate and deposition rate occurs at the peak flow. Spatial distribution of bed erosion and deposition in the river channel after the GLOF