Matlab Codes for Volcanology - Magma Viscosity

By Tom D Pering

University of Sheffield, England

Download (M)

Licensed according to this deed.

Published on

Abstract

A simple Matlab function for calculating the viscosity of magma. It is based on the method of Shaw (1972). Viscosities of magmatic silicate liquids: an empirical method of prediction. American Journal of Science 272, pp. 870-893

It requires the inputs:
SiO2
TiO2
Al2O3
Fe2O3
MnO
MgO
CaO
Na2O
K2O
P2O5
H2O - water content
mTc - magma temperature in celsius

The outputs are:
viscosity_gas = viscosity of the magma with the water content

viscosity_no_gas = viscosity of the magma with no gas

totalweight = total weight of all the inputs

During creation I also used the excel spreadsheet by John D. Winter as a guide (Available at: http://www.whitman.edu/geology/winter/) - which is called "Magma Viscosity spreadsheet". As per my magma density code, the excel spreadsheet is probably more user friendly, however, it is useful to have this file in Matlab. It is also useful to use in combination with modelling work.

If you have any questions please don't hesitate to contact me.

Cite this work

Researchers should cite this work as follows:

  • Tom D Pering (2013), "Matlab Codes for Volcanology - Magma Viscosity," https://vhub.org/resources/2792.

    BibTex | EndNote

Tags

  1. MATLAB
  2. magma viscosity