Home > Research > Publications & Outputs > The influence of cross-sectional channel geomet...
View graph of relations

The influence of cross-sectional channel geometry on rheology and flux estimates for active lava flows

Research output: Contribution to Journal/MagazineJournal articlepeer-review

Published
Article number829
<mark>Journal publication date</mark>06/2014
<mark>Journal</mark>Bulletin of Volcanology
Volume76
Number of pages15
Publication StatusPublished
Early online date26/06/14
<mark>Original language</mark>English

Abstract

Lava rheology and effusion rate are critical factors in determining the evolution of lava flows. However, direct and accurate field measurements are difficult to carry out, and estimates are usually based on measurements of the flows surface velocity and assumptions of sub-surface geometry. Using numerical flow models, we show that the potential for error due to geometry uncertainty is minimized if a semi-elliptical cross sectional channel shape is assumed. Flow is simulated for isothermal Newtonian, temperature-dependent Newtonian, and isothermal power-law rheology lavas. For isothermal Newtonian lava, we find that error in channel shape alone can make apparent viscosity estimates 3.5 times too large (e.g. for inappropriate use of the Jeffreys equation on a narrow semi-elliptical channel). For a temperature-dependent rheology, using
an analytical approximation for Newtonian flow in a semi-elliptical geometry yields apparent viscosity and flux values that are more accurate than estimates which assume a rectangular geometry, for all channel shapes considered, including rectangular channels. Viscosity calculations for real channels on Mauna Loa and Mount Etna show that for a Newtonian rheology, a semi-elliptical analytical solution gives an approximation 3 times closer to the actual viscosity than a rectangle with the same depth while, if the lava
is shear-thinning (power law exponent m = 0.6), a rectangular approximation is 15% more accurate. Our results can be used to bracket possible viscosity and flux estimates when channel topography is poorly constrained.