mzSaturation {GCDkit} | R Documentation |
Calculates monazite saturation temperatures for given major-element compositions and LREE contents of the magma.
mzSaturation(cats = milli, REE = filterOut(WR, c("La", "Ce", "Pr","Nd","Sm", "Gd"), 1), H2O = 3, Xmz = 0)
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 |
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
Returns a matrix 'results
' with the following components:
Dmz |
distribution coefficient |
Tmz.sat.C |
monazite saturation temperature |
Saturation.r
Vojtech Janousek, vojtech.janousek@geology.cz
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