Phase Equilibria in Silicate Systems Intro Petrol EPSC212

Phase Equilibria in Silicate Systems Intro. Petrol. EPSC-212, Francis-14

Winter, J. D. ; 2001: An Introduction to Igneous and Metamorphic Petrology. Prentice Hall, Chapter 4, 2001. Fundamentals: Units Types of units commonly used: Weight units: = gms species/ 100 gms (X-Ray analytical techniques determine weight fraction element) • Molar Units: = gms species / (molecular wt. of species), normalized to 100 (most chemical phenomena are proportional to molecular proportions) • Atomic units: = (no. moles species) x (no. atoms per species), normalized to 100 atoms • Cation units: = (no. moles species) x (no. cations per species), typically normalized to some total number of cations (or anions). • Oxygen units: = (no. moles species) x (no. oxygens per species), commonly normalized to some total number of anions (closest to volume proportions) Example: Coordinates of Enstatite (Mg. Si. O 3) in: Mg 2 Si. O 4 - Si. O 2 space Weight units: 70 30 Molecular units: 50 50 Cations units: 75 25 Oxygen units 66. 7 33. 3

Volatile have an importance beyond that predicted simply by their abundance because: - Volatiles have low molecular weights: H 2 O = 18 CO 2 = 44 Si. O 4 = 92 Na. Al. Si 3 O 8 = 262 In a melt consisting of Na. Al. Si 3 O 8 clusters and H 2 O molecules: 0. 5 wt. % H 20 ~ 45 mole % H 2 O Small amounts of water produce large effects because of its low molecular wt. compared to that of a silicate magma. This effect is enhanced by the fact that at XH 2 O < 0. 3, molecular water dissociates into 2 OH’s f. H 20 ~ PH 2 O α XH 2 O. 2 H 2 O + Obridging 2 × OH

The Rock Forming Minerals: A mineral is defined as a naturally occurring crystalline phase or compound that is made up of a 3 -D ordered atomic arrangement or structure of different atoms. The dominant rock forming minerals are silicates, compounds of Si and O, because Si is the most abundant metal and oxygen the most abundant anion in the Earth's crust. Minerals can be thought of as close packings of oxygens to a first approximation, because of the large size of O, with the smaller Si and other metals occupying the interstices, or sites, between the oxygens. Mineral Name: Olivine Pyroxene Feldspar Structural Formula: Y 2 TO 4 XYT 2 O 6 WT 4 O 8 (Mg, Fe)2 Si. O 4 Ca(Mg, Fe)Si 2 O 6 (K, Na)Al. Si 3 O 8 Chemical Composition:

Rock Forming Minerals: The ratio of Si to O determines the type and abundance of the silicate mineral(s) present, and the other metals distribute themselves in sites between the oxygens according to the size and charge of their cations. The amounts of different elements in any given mineral may vary somewhat, but the type and proportion of occupied sites is fixed by the structure of the silicate mineral and can be expressed by a formula of the type: Ww Xx Y y T t Aa Where the capital letters stand for occupied sites of different co-ordination number and the small letters are small whole numbers. A = anion site W = 12 co-ordinated site X = 8 co-ordinated site Y = 6 co-ordinated site T = 4 co-ordinated site

K, Rb, Ba, Na Ca, Mn, Na Mg, Fe, Mn, Al, Ti Si, Al

Rock Forming Minerals Silicates (lithophile elements) Olivine: Y 2 TO 4 Mg 2 Si. O 4 - Fe 2 Si. O 4 Garnet: X 3 Y 2(TO 4)3 Mg 3 Al 2(Si. O 4)3 - Fe 3 Al 2(Si. O 4)3 Clinopyroxene: XYT 2 O 6 Ca. Mg. Si 2 O 6 - Ca. Fe. Si 2 O 6 Orthopyroxene: Y 2 T 2 O 6 Mg 2 Si 2 O 6 - Fe 2 Si 2 O 6 Amphibole: W 01 X 2 Y 5 T 8 O 22(OH)2 Ca 2 Mg 5 Si 8 O 22(OH)2 Na. Ca 2 Mg 4 Al(Al 2 Si 6 O 22(OH)2 Feldspar: WT 4 O 8 Ca. Al 2 Si 2 O 8 - Na. Al. Si 3 O 8 - KAl. Si 3 O 8 Quartz: TO 2 Si. O 2 Micas: W Y 2 -3 T 4 O 10(OH)2 KAl 2(Al. Si 3)O 10(OH)2 - KMg 3(Al. Si 3)O 10(OH)2 Clays & Serpentine: Y 3 T 2 O 5(OH)4 Mg 3 Si 2 O 5(OH)4 Carbonates: (lithophile elements) Calcite - Dolomite: YCO 3 Ca. CO 3 - Ca(Mg, Fe)(CO 3)2 Other Oxides: (lithophile elements) Spinel – Magnetite: Y 3 O 4 Mg. Al 2 O 4 - Fe. Cr 2 O 4 - Fe 3 O 4 Hematite – Ilmenite: Y 2 O 3 Fe 2 O 3 - Fe. Ti. O 3 Pyrite - Pyrrhotite Fe. S 2 - Fe 1 -x. S Sulfides: (chalcophile elements)

The Mineralogical Phase Rule In any chemical system at equilibrium, the following relationship holds: FDegrees of Freedom = Components - Phases + 2 F equals the minimum number of variables that must be specified in order to completely define the state of a system. F is thus the variance of the system, the number of unknowns. C equals the number of independent chemical components needed to define the composition of the system. P equals the number of physical phases present in the system, which include the number of solid minerals, plus liquid and gas phases, if present. 2 represents the variables pressure and temperature.

Single Component Systems: Si. O 2 When a solid consist of 2 coexisting minerals (phases): F = C–P+2 = 1 -2+2 = 1 Such a system is invariant at any given pressure, and thus a single component solid phase will melt at 1 unique temperature at any specified pressure. The boundary between the 2 phases in P - T space will be a univariant line with a slope approximated by: d G = - Sd. T + Vd. P = 0 d. P/d. T = S/ V This is also true for solid - liquid phase boundaries because, to a first approximation, Ho and So are constant for small changes in temperature (true for all reactions not involving a relatively compressible vapour phase).

Two Component Systems: Mg 2 Si. O 4 – Si. O 2 Pure forsterite melts at 2163 o. C at 1 atm. If extra Si. O 2 is added to the system, Si. O 2 will be present only in the melt phase, while forsterite will remain a pure phase. The temperature at which forsterite crystallizes is now an inverse function of Si. O 2 content: at equilibrium: GFo. Ol = GFo. Liq and G Fo = 0 Go. Fo + R × T × Ln(a. Fo. Ol) = Go. Fo. Liq + R × T × Ln(a. Fo. Liq) Go Fo = - R × T × Ln(a. Fo. Liq/a. Fo. Ol) or Go Fo Fo = - R × T × Ln(1 -XSi. O 2), a. Fo. Ol = 1. 0, assume activity (a) = mole fraction X) and = Ho. Fo – T × So. Fo = Ho. Fo – T × Ho. Fo / TFo If Ho and So are insensitive to small changes in T & P, then: Ho Fo – T × Ho. Fo / TFo = - R × T × Ln(1 -XSi. O 2) or Ln(1 -XSi. O 2) = Ho. Fo / R × (T - TFo) / (T×TFo) van't Hoff equation for melting point depression

Two Component Systems: Binary Phase Diagrams If a system is comprised of 2 components, then where a solid and liquid phase coexist: F =2 -2+2=2 Such a system is univariant at any given pressure, and thus the melting point of a solid will depend on the proportions of the two components. In the absence of extensive solid solution, the presence of an additional component will reduce the melting temperature of single component solid phases because the additional component typically dissolves preferentially in the liquid phase. e No solid-solution Between end-members The Eutectic point “e” of a two component system is invariant (F = 0, if pressure fixed) and is defined by the intersection of two univariant (F = 1) liquidus curves, originating from the melting temperatures of the two pure end-member phases.

Peritectic versus Eutectic Invariant points e 2 e 1 p - oliv + liq e - liq opx e 1 - liq albite + qtz opx + qtz e 2 - liq neph + albite

The Lever Rule X For bulk Composition X T 3 : Forst / liquid : b / a; Forst / whole = b /(a+b) T 2 : Forst / liquid : c / a; Forst / whole = c /(a+c) T 1 : Forst / Enst Forst / whole = d /(a+d) : d / a;

Cumulate Rocks versus Rocks that represent liquids Liquids vs Cumulates Equilibrium vs 1557 Fractional Crystallization vs Partial Melting Upon cooling to 1557 o. C, early crystallized olivine exhibits a reaction relationship with the residual liquid of composition “p” to form orthopyroxene. Either olivine or melt must disappear before cooling can continue. During partial melting, orthopyroxene begins to melt incongruently at 1557 o. C to form olivine plus a liquid of composition “p”. Orthopyroxene must be consumed before the temperature can increase.

The presence of other components in solid solution at levels that are insufficient to stabilize a separate phase destroys the invariant nature of melting. The temperature and composition of the first melt are determined by the amount of the additional component. During partial melting, these additional components are typically the first to be refined out into the melt.

Bianary Systems with extensive Solid Solution: Olivine exhibits complete solid solution between the forsterite (Mg 2 Si. O 4) and fayalite (Fe 2 Si. O 4) end-members. In Fe and Mg bearing systems, neither the olivine solid nor the olivine liquid are pure end-member components: We now have two van't Hoff equations: Ln(XFo. Liq / XFo. Ol) = Ho. Fo / R × (T - TFo) / (T × TFo) Ln(XFa. Liq / XFa. Ol) = Ho. Fa /R × (T - TFa) /(T × TFa) Because: XFa. Lig = 1 -XFo. Liq and XFa. Ol = 1 -XFo. Ol Then: Ln(XFo. Liq / XFo. Ol) = Ho. Fo / R × (T - TFo) / (T × TFo) Ln(1 -XFo. Liq / (1 -XFo. Ol)) = Ho. Fa /R × (T - TFa) / (T × TFa) The choice of any T between TFo and TFa will enable the calculation of the compositions of the coexisting olivine and liquid for that T, and thus the solidus and liquidus at any T. Exactly analogous solid solution relationships can be developed for the plagioclase series feldspars: anorthite Ca. Al 2 Si 2 O 8 - albite Na. Al. Si 3 O 8

Ternary Systems: Forsterite – Diopside – Anorthite Liquidus Projection In order to portray the magmatic phase relations of systems with more than two chemical components, we need to develop specialized projection schemes. Three component systems can be represented on a two dimensional sheet of paper, if we project only those phase relationships for which a magmatic liquid is present. P = 1 atm

Ternary Systems: Forsterite – Diopside - Anorthite A Liquid of bulk composition X cools to the olivine liquidus surface at 1600°C, at which point Forsterite begins to crystalize P = 1 atm The liquid composition moves directly away from Fo, producing a dunite cumulate, until it reaches the cotectic, at which point Diopside begins to crystallize with Forsterite. X

Forsterite – Diopside - Anorthite P = 1 atm Co-precipitation of Forsterite + Diopside causes liquid composition to move down cotectic curve, producing a wehrlite cumulate. X

Forsterite – Diopside - Anorthite P = 1 atm The liquid composition reaches the ternary eutectic at 1270°C, at which point Anorthite begins to crystallize with Diopside and Forsterite, producing a gabbroic cumulate. The composition of the liquid remains at the eutectic point until all the liquid is consumed. X

Ternary Systems The composition of the first melt of an assemblage ABD is that of invariant eutectic point e. ABD, while the composition of the first melt of assemblage DBC is that of invariant eutectic point e. DBC. The intersection of a univariant curve with the Alkemade line joining the compositions of the coexisting solid phases defines a thermal maximum along the univariant curve.

Ternary Systems with Solid Solution:

Oliv - Cpx - Qtz Liquidus Projection: The invariant point “p”, at which olivine, clinopyroxene and orthopyroxene coexist with a liquid, is a peritectic point because it lies outside of the compositional volume of the solid phases. It represents the first melt of any assemblage consisting of olivine, opx, and cpx (mantle peridotite) at 1 atm. , and is analogous in composition to a quartz-normative basalt. Similarly, the univariant curve along which olivine and orthopyroxene coexist with a liquid is a reaction curve because the tangent to the curve at any point cuts the olivine - orthopyroxene Alkemade line with a negative olivine intercept. The invariant point “e” is a eutectic and represents the composition of the first melt of an assemblage of quartz-diopside-orthopyroxene. The composition of “e” approximates that of the Earth’s continental crust.

Mantle Ocean Continent crust Si. O 2 Ti. O 2 Al 2 O 3 Mg. O Fe. O Ca. O Na 2 O K 2 O Total 45. 2 0. 7 3. 5 37. 5 8. 5 3. 1 0. 6 0. 1 99. 2 49. 4 15. 4 7. 6 10. 1 12. 5 2. 6 0. 3 99. 3 60. 3 1. 0 15. 6 3. 9 7. 2 5. 8 3. 2 2. 5 99. 5 Cations normalized to 100 cations Si Ti Al Mg Fe Ca Na K O 38. 5 0. 5 3. 6 47. 6 6. 0 2. 8 0. 9 0. 1 140. 2 46. 1 1. 0 16. 9 10. 6 7. 9 12. 5 4. 7 0. 5 153. 0 56. 4 0. 7 17. 2 5. 4 5. 6 5. 8 3. 0 161. 3 Mineralogy (oxygen units, XFe 3+ = 0. 10) Quartz Feldspar Clinopyroxene Orthopyroxene Olivine Oxides 0. 0 13. 2 6. 7 18. 3 59. 9 1. 8 0. 0 57. 3 25. 7 4. 1 9. 9 3. 0 13. 0 64. 3 5. 9 14. 7 0. 0 2. 0 Oceanic crust - MORB basalt Continental crust - granite p e

Liquidus Projections for haplo-basalts The Basalt Tetrahedron at 1 atm: The olivine - clinopyroxene plagioclase plane is a thermal divide in the haplo-basalt system at low pressures and separates natural magmas into two fundamentally different magmatic series. Sub-alkaline basaltic magmas with compositions to the Qtz-rich side of the plane fractionate towards Qtz-saturated residual liquids, such as rhyolite. Alkaline basaltic magmas with compositions to the Qtz-poor side of the plane fractionate towards residual liquids saturated in a feldspathoid, such as nepheline phonolite.

Since the dominant mineral in the mantle source of basaltic magmas is olivine, we can achieve a further simplification by projecting the liquidus of basaltic systems from the perspective of olivine:

Alkaline basalts fall to the Foid-side of the olivineclinopyroxene-plagioclase plane (1 atm thermal divide) and fractionate to foid-saturated residual liquids. Sub-alkaline basalts fall to the Quartz-side and fractionate towards quartz -saturated residual liquids. Alkaline basaltic lavas are volumetrically insignificant (~1%), but strongly enriched in highly incompatible trace elements profiles compared to sub-alkaline lavas, and low in HREE, Y, & Sc. These characteristics are generally ascribed to small degrees of partial melting at elevated pressures, leaving garnet as a phase in the refractory residue.

The Effect of Pressure 1 atm Increasing pressure shifts the olivcpx-opx peritectic point towards less Si-rich compositions. At approximately 10 kbs this invariant point moves into the oliv - cpx- opx compositional volume, and the first melt of the mantle has an olivine basalt composition. The invariant point is still a peritectic point, however, because of the extensive solid solution of cpx towards opx. At pressures exceeding 15 -20 kbs, this invariant point moves outside the simple olivine - cpx - qtz system, into the Neph-normative volume of the basalt tetrahedron. The first melt of mantle peridotite is an alkaline olivine basalt at these high pressures.

Since the dominant mineral in the mantle source of basaltic magmas is olivine, we can achieve a further simplification by projecting the liquidus of basaltic systems from the perspective of olivine: Movement of the invariant point determining the composition of the first melt with increasing pressure.