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, excluding Eu |
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 REEi: La, Ce, Pr, Nd, Sm, Gd.
Dmz = (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 of Montel (1993) in °C |
FM |
cationic ratios |
Saturation.r
Vojtěch Janoušek, vojtech.janousek@geology.cz
Montel JM (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
Kelsey DE, Clark C, Hand M (2008) Thermobarometric modelling of zircon and monazite growth in melt-bearing systems: examples using model metapelitic and metapsammitic granulites. J Metamorph Geol 26: 199-212 doi: 10.1111/j.1525-1314.2007.00757.x
Rapp RP, Watson EB (1986) Monazite solubility and dissolution kinetics; implications for the thorium and light rare earth chemistry of felsic magmas. Contrib Mineral Petrol 94:304-316 doi: 10.1007/BF00371439
data(blatna) accessVar("blatna") mzSaturation(Xmz=0.8)