trace element geochemistry lecture 24. geochemical classification

17
Trace Element Geochemistry Lecture 24

Upload: curtis-norton

Post on 21-Jan-2016

216 views

Category:

Documents


0 download

TRANSCRIPT

Page 1: Trace Element Geochemistry Lecture 24. Geochemical Classification

Trace Element Geochemistry

Lecture 24

Page 2: Trace Element Geochemistry Lecture 24. Geochemical Classification

Geochemical Classification

Page 3: Trace Element Geochemistry Lecture 24. Geochemical Classification

The Rare Earth Elements• The rare earths and Y are strongly

electropositive. As a result, they form predominantly ionic bonds, and behave as hard charged spheres.

• The lanthanide rare earths are in the +3 valence state over a wide range of oxygen fugacities.

• In the transition metals, the s orbital of the outermost shell is filled before filling of lower electron shells is complete so the configuration of the valence electrons is similar in all the rare earth, hence all exhibit similar chemical behavior.

• Ionic radius, which decreases progressively from La3to Lu3+ (93 pm) governs their relative behavior.  

• Because of their high charge and large radii, the rare earths are incompatible elements. o The degree of incompatibility of the lanthanides varies with

atomic number. Highly charged U and Th are highly incompatible elements, as are the lightest rare earths.

o However, the heavy rare earths have sufficiently small radii that they can be accommodated to some degree in many common minerals such as Lu for Al in garnet. Eu2+ can substitute for Ca in plagioclase.

Page 4: Trace Element Geochemistry Lecture 24. Geochemical Classification

Rare Earth Diagrams• The systematic variation in lathanide rare

earth behavior is best illustrated by plotting the log of the relative abundances as a function of atomic. Relative abundances are calculated by dividing the concentration of each rare earth by its concentration in a set of normalizing values, such as the concentrations of rare earths in chondritic meteorites.o Rare earths are also refractory elements, so that their relative

abundances are the same in most primitive meteorites - and presumably (to a first approximation) in the Earth.

• Why do we use relative abundances? The abundances of even-numbered elements in the solar system are greater than those of neighboring odd-numbered elements and abundances generally decrease with increasing atomic number, leading to a saw-toothed abundance pattern. Normalizing eliminates this.

• Abundances in chondritic meteorites are generally used for normalization. However, other normalizations are possible: sediments (and waters) are often normalized to average shale.

Page 5: Trace Element Geochemistry Lecture 24. Geochemical Classification

Partition Coefficients

Page 6: Trace Element Geochemistry Lecture 24. Geochemical Classification

The Partition Coefficient• Geochemists find it convenient to define a partition or

distribution coefficient of element i between phases α and β:

• Where one phase is a liquid, the convention is the solid is placed on top:

• Incompatible elements are those with Ds/l ≪ 1. Compatible elements are those with Ds/l ≥ 1. These terms refer to partitioning between silicate melts and phases common to mantle rocks (peridotite). It is this phase assemblage that dictates whether lithophile trace elements are concentrated in the Earth’s crust, hence the significance of these terms.

Page 7: Trace Element Geochemistry Lecture 24. Geochemical Classification

Thermodynamic Basis• The chemical potentials of element i in phases α and β are

• At equilibrium:

o Since trace elements obey Henry’s Law, we can replace the activity coefficient with h.

• The left hand side is ∆G˚, so that

• and

• where C is simply a constant converting concentration units (usually ppm) to mole fraction.

• Thus the distribution coefficient is a kind of equilibrium constant.

Page 8: Trace Element Geochemistry Lecture 24. Geochemical Classification

Relationship among distribution coefficients

• In a system with three phases, α, β, and γ, if α and β are in equilibrium and α and γ are in equilibrium, then β and γ must also be in equilibrium. It follows that:

Dα-β=Dα-γ/Dα-β

• This relationship has practical use. For example, if we can determine the partition coefficient for an element between pyroxene and melt and between garnet and pyroxene, we can then calculate the garnet–melt partition coefficient for this element.

Page 9: Trace Element Geochemistry Lecture 24. Geochemical Classification

Temperature and Pressure Dependence

• In ideal solutions, the temperature dependence of the partition coefficient is the same as that of the equilibrium constant:

• ∆G can be expanded into entropy and enthalpy terms:

• In ideal solution, and assuming again that ∆V is independent of temperature and pressure, the pressure dependence is also the same as that of the equilibrium constant:

• We would predict a strong pressure dependence when the ionic radius of an element differs greatly from that of the available crystal lattice site. Thus, for example, we would predict the partition coefficient for K between pyroxene and melt would be strongly pressure-dependent since the ionic radius of K is 150 pm and is much larger that the size of the M2 site in clinopyroxene, which is normally occupied by Ca, with a radius of about 100 pm. Conversely, where the size difference is small (e.g., Mn (83 pm) substituting for Fe (78 pm)), we would expect the pressure dependence to be smaller.

• In non-ideal solutions, the T and P dependencies will be more complex because the activity coefficients (the Henry’s Law coefficients) will also depend on T and P.

• Bottom line: partition coefficient is temperature and pressure dependent.

Page 10: Trace Element Geochemistry Lecture 24. Geochemical Classification

Importance of Ionic Size and Charge

• Much of the interest in trace elements in igneous processes centers on the elements located in the lower left portion of the periodic table (K, Rb, Cs, Sr and Ba; the rare earths, Y, Zr, Nb, Hf and Ta).

• One reason for this focus of attention is that these elements are all lithophile and therefore present at relatively high abundance in the Earth’s crust and mantle. There is another reason, however: their chemical behavior is comparatively simple: their behavior in igneous systems is mainly (not entirely) a function of ionic size and charge.

• The other trace elements that receive the most attention from igneous geochemists are the first transition series elements. Though their electronic structures and bonding behavior are considerably more complex, charge and size are also important. Many of these elements, particularly Ni, Co, and Cr, have partition coefficients greater than 1 in many Mg–Fe silicate minerals. Hence the term “compatible elements” often refers to these elements.

Ionic radius (picometers) vs. ionic charge contoured for clinopyroxene/liquid partition coefficients. Cations normally present in clinopyroxene M1 and M2 sites are Ca2+, Mg2+, and Fe2+, shown by ✱ symbols. Elements whose charge and ionic radius most closely match that of the major elements have the highest partition coefficients

Page 11: Trace Element Geochemistry Lecture 24. Geochemical Classification

Quantitative Treatment• Consider the substitution reaction:

Mℓ + CaMgSi2O6 ⇄ Caℓ + MMgSi2O6

• The Gibbs free energy change of this reaction can be expressed as:

• The first term is ∆G for transferring an M2+ ion from the melt to the crystal lattice and simultaneously transferring a Ca2+ ion from the lattice site to the liquid.

• The second term is the ∆G associated with the melting of diopside, and governs the distribution of Ca between diopside and the liquid. The distribution coefficient for element M then depends on these two components of free energy:

Page 12: Trace Element Geochemistry Lecture 24. Geochemical Classification

Substitution & Lattice Strain

• If we substitute an atom that larger than one normally occupying a site in a crystal lattice, for example a Ba2+ for a Mg2+ in olivine, the lattice will strain to accommodate that atom.

• Strain is a form of stored energy.

• We can express this as?

Olivine Lattice

Page 13: Trace Element Geochemistry Lecture 24. Geochemical Classification

Strain Theory• According to the lattice strain energy theory, ∆Gexchange is dominated by the

energy associated with the lattice strain resulting from M2+ being a different size than Ca2+. Because the melt (at least at low pressure) has a far less rigid structure and is more compressible than the solid, any strain in the melt is essentially negligible compared with the strain in the solid. In other words:

• The strain energy, ∆Gstrain, may be calculated as:

• where r0 is the optimal radius of the lattice site, rM is the ionic radius of M, NA is Avogadro’s Number, and E is Young’s modulus, which is the ratio of stress applied to the resulting strain (change in dimension) and has units of pressure. It is a property of the material and is related to compressibility.

Page 14: Trace Element Geochemistry Lecture 24. Geochemical Classification

Blundy & Wood Model• We can now substitute this expression into

• and obtain

• where

• Where charge is different from the ion normally occupying the site, we much consider the coupled substitution. For example, we can balance charge by substituting both La3+ and Na+ into the Ca2+ site in diopside. In this case, our relationship is:

• The point is we can build a theoretical framework to predict how partitioning will depend on ionic size and charge (and the T and P dependency).

Page 15: Trace Element Geochemistry Lecture 24. Geochemical Classification

Blundy & Wood Model

We still depend heavily on experimental determination of partition coefficients. But by experimentally determining a few, we can use the Blundy & Wood model to predict others (as well as pressure dependence).

Page 16: Trace Element Geochemistry Lecture 24. Geochemical Classification

Dependence on Composition

Ol-liquid partition coefficient for Zn parameterized as a function of non-bridging oxygens in the melt

and T.pyx-liquid partition coefficient for

Sm as a function of liquid and solid composition

Page 17: Trace Element Geochemistry Lecture 24. Geochemical Classification

Rare Earth Partition Coefficients