mzSaturation {GCDkit}R Documentation

Monazite saturation (Montel 1993)

Description

Calculates monazite saturation temperatures for given major-element compositions and LREE contents of the magma.

Usage

mzSaturation(cats = milli, 
    REE = filterOut(WR, c("La", "Ce", "Pr","Nd","Sm", "Gd"), 1), 
    H2O = 3, Xmz = 0)

Arguments

cats

numeric matrix; whole-rock data recast to millications

REE

numeric matrix with LREE concentrations - only complete set of La-Gd

H2O

assumed water contents of the magma

Xmz

mole fractions of the REE-phosphates in monazite

Details

This function uses saturation model of Montel (1993). The formulae are as follows:

LREE = sum(REEi/at.weight(REEi))/Xmz

where REE\_i: La, Ce, Pr, Nd, Sm, Gd.

Dmz = 100 (Na + K + 2Ca)/Al * 1/(Al + Si)

T.calc = 13318/(9.5 + 2.34 Dmz + 0.3879 sqrt(H2O) - ln(LREE)) - 273.15

Value

Returns a matrix 'results' with the following components:

Dmz

distribution coefficient

Tmz.sat.C

monazite saturation temperature

Plugin

Saturation.r

Author(s)

Vojtech Janousek, vojtech.janousek@geology.cz

References

Montel J M (1993) A model for monazite/melt equilibrium and application to the generation of granitic magmas. Chem Geol 110: 127-146 doi: 10.1016/0009-2541(93)90250-M


[Package GCDkit version 4.0 Index]