Numerical simulation of lava flows based on depth-averaged equations

Risks and damages associated with lava flows
propagation (for instance the most recent Etna eruptions) require a quantitative description of this phenomenon and a reliable forecasting of lava flow paths. Due to the high complexity of these processes, numerical solution of the complete conservation equations for real lava flows is often practically impossible. To overcome the computational difficulties, simplified models are usually adopted, including 1-D models and cellular automata. In this work we propose a simplified 2D model based on the conservation equations for lava thickness and depthaveraged velocities and temperature which result in first order partial differential equations. The proposed approach represents a good compromise between the full 3-D description and the need to decrease the computational time. The method was satisfactorily applied to reproduce some analytical solutions and to simulate a real lava flow event occurred during the 1991–93 Etna eruption.

Alcrudo, F., and F. Benkhaldoun (2001), Exact solutions to the Riemann of the shallow water equations with a step,Comput. Fluids, 30, 643–671.
Ambrosi, D. (1999), Approximation of shallow water equations by Riemann solvers, Int. J. Numer. Methods Fluids, 20, 157– 168.
Barberi, F., M. Carapezza, M. Valenza, and L. Villari (1993), The control of lava flow during the 1991–1992 eruption of Mt. Etna, J. Volcanol. Geotherm. Res., 56, 1 – 34.
Burguete, J., P. Garcia-Navarro, and R. Aliod (2002), Numerical simulation of runoff from extreme rainfall events in a mountain water catchment, Nat. Hazards Earth Syst. Sci., 2, 109– 117.
Calvari, S., M. Coltelli, M. Neri, M. Pompilio, and V. Scribano (1994), The 1991– 93 Etna eruption: Chronology and flow-field evolution, Acta Vulcanol., 4, 1 –15.
Chinnayya, A., A. LeRoux, and N. Seguin (2004), A well-balanced numerical scheme for the approximation of the shallow-water equations with topography: The resonance phenomenon, Int. J. Finite Volumes.
Costa, A., and G. Macedonio (2002), Nonlinear phenomena in fluids with temperature-dependent viscosity: An hysteresis model for magma flow in conduits, Geophys. Res. Lett., 29(10), 1402, doi:10.1029/2001GL014493.
Costa, A., and G. Macedonio (2003), Viscous heating in fluids with temperature-dependent viscosity: Implications for magma flows, Nonlinear Proc. Geophys., 10, 545– 555.
Crisp, J., and S. Baloga (1990), A model for lava flows with two thermal components, J. Geophys. Res., 95, 1255– 1270.
Ferrari, S., and F. Saleri (2004), A new two-dimensional shallow water model including pressure effects and slow varying bottom topography, ESAIM Math. Modell. Numer. Anal., 38, 211– 234.
George, D. (2004), Numerical approximation of the nonlinear shallow water equations with topography and dry beds: A Godunov-type scheme, M.S. thesis, Univ. of Wash., Seattle.
Gerbeau, J., and B. Perthame (2001), Derivation of viscous Saint-Venant system for laminar shallow water: Numerical validation, Discrete Continuos Dyn. Syst., Ser. B, 1, 89– 102.
Heinrich, P., A. Piatanesi, and H. Hebert (2001),Numerical modelling of tsunami generation and propagation from submarine slumps: The 1998 Papua New Guinea event, Geophys. J. Int., 145, 97–111.
Keszthely, L., and S. Self (1998), Some physical requirements for the emplacement of long basaltic lava flows, J. Geophys. Res., 103, 27, 447– 27,464.
Figure 4. Simulated lava thickness of the 3rd and 4th January 1992 Etna lava flow.
LeVeque, R. (1998), Balancing source terms and flux gradients in highresolution Godunov methods: The quasi-steady wave-propagation algorithm, J. Comput. Phys., 146, 346– 365.
LeVeque, R. (2002), Finite Volume Methods for Hyperbolic Problems, Cambridge Univ. Press, New York.
Monthe, L., F. Benkhaldoun and I. Elmahi (1999), Positivity preserving finite volume Roe schemes for transport-diffusion equations, Comput. Methods Appl. Mech. Eng., 178, 215– 232.
Neri, A. (1998), A local heat transfer analysis of lava cooling in the atmosphere: Application to thermal diffusion-dominated lava flows, J. Volcanol. Geotherm. Res., 81, 215–243.
Pieri, D., and S. Baloga (1986), Eruption rate, area, and length relationships for some Hawaiian lava flows, J. Volcanol. Geotherm. Res., 30, 29– 45.
Shah, Y., and J. Pearson (1974), Stability of non-isothermal flow in channels: III.Temperature-dependent pawer-law fluids with heat generation, Chem. Eng. Sci., 29, 1485–1493.