*Earth Planets Space,* Vol. 62 (No. 5), pp. 483-488, 2010

doi:10.5047/eps.2010.02.005

## A simple formula for calculating porosity of magma in volcanic conduits during dome-forming eruptions

Tomofumi Kozono^{1} and Takehiro Koyaguchi^{2}

^{1}National Research Institute for Earth Science and Disaster Prevention, 3-1, Tennodai, Tsukuba, Ibaraki 305-0006, Japan

^{2}Earthquake Research Institute, University of Tokyo, 1-1-1, Yayoi, Bunkyo-ku, Tokyo 113-0032, Japan

(Received December 14, 2009; Revised February 12, 2010; Accepted February 15, 2010; Online published July 9, 2010)

**Abstract: **
We present a simple formula for analyzing factors that govern porosity of magma in dome-forming eruptions. The formula is based on a 1-dimensional steady conduit flow model with vertical gas escape, and provides the value of the porosity as a function of magma flow rate, magma properties (viscosity and permeability), and pressure. The porosity for a given pressure depends on two non-dimensional numbers ε and θ. The parameter ε represents the ratio of wall friction force to liquid-gas interaction force, and is proportional to the magma viscosity. The parameter θ represents the ratio of gravitational load to liquid-gas interaction force and is inversely proportional to the magma flow rate. Gas escape is promoted and porosity decreases with increasing ε or θ. From the possible ranges of ε and θ for typical magmatic conditions, it is inferred that the porosity is primarily determined by ε at the atmospheric pressure (near the surface), and by θ at higher pressures (in the subsurface region inside the conduit). The porosity near the surface approaches 0 owing to high magma viscosity regardless of the magnitude of the magma flow rate, whereas the subsurface porosity increases to more than 0.5 with increasing magma flow rate.

**Key words: **
Conduit flow, dome-forming eruptions, magma porosity, gas escape from magma.

Corresponding author E-mail: kozono@bosai.go.jp

[Full text] (PDF 464 KB)