TEXTURAL AND MICROANALYSIS OF IGNEOUS ROCKS: TOOLS ...

calculated melt concentrations for Sr and Ba (within their analytical uncertainty).
Likewise, Bindeman et al. [7] showed that variations of ∼150 ◦ C produce < 10%
differences for particular partition coefficients, which were often within error of
the respective D values.
Blundy and Wood [11] presented an alternative model for quantitative predic-
tion of D values based upon thermodynamic principles and lattice strain theory
(i.e., [15]). Specifically, they noted the dominant roles of cation size, lattice site
size, and elasticity of the lattice site related to cation substitution in minerals.
The elastic response of a lattice site to cation substitution is measured by Young’s
Modulus (E), and the value of E plag is dependent upon plagioclase An content
and the charge of the substituent cation [11, 12]. Equation 4.2 (Equation 2 of
Blundy and Wood [11]) is the basis of their predictive D model that I apply to
plagioclase,
Di = D z+
o
∗ exp
−4πEplagNA
ro
2 (r2 o − r2 i ) + 1
3 (r3 i − r3 o)
RT
(4.2)
where Di of a trace element (i) is a function of the strain compensated partition
coefficient for the divalent (D 2+
0
) or trivalent (D3+ 0 ) substituent cation, Young’s
Modulus for the plagioclase M lattice site (E plag), NA (Avogadro’s Number), the
radius of the plagioclase M site (ro), the radius of the substituent cation (ri),
the gas constant (R), and crystallization temperature (T). For a more thorough
discussion of this model see Blundy and Wood [11, 12].
4.3.6.1 Partition coefficients for Elan Bank Plagioclase Crystals
I used equation 4.2 and the parameterizations of Blundy and Wood [11, 12]
to calculate D values for Rb, Sr, Y, Ba, La, Ce, Pr, Nd, and Sm, as partitioning
183