Sodium partitioning between clinopyroxene and silicate melts

We present a quantitative model for Na partitioning between clinopyroxene and silicate melt, applicable in the range 0.1 MPa to 6.0 GPa, 1000–2000°C over a wide span of composition. Our model is derived from the Jadeite melting equilibrium, via the relationship ΔG_ƒ(P_,T₎ = −RTlnK_ƒ = RTlnD*_Na where ΔG_ƒ(P_,T₎ is the Gibbs free energy of fusion of Jadeite at the pressure (P) and temperature (T) of interest, T is in kelvins, R is the gas constant, K_ƒ is the equilibrium constant for the the melting reaction and D*_Na is a molar partition coefficient, defined here as the molar ratio of Na in crystal to Na in melt on a six-oxygen basis. An expression for ΔG_ƒ(P_,T₎ is obtained from published experimental data on the Jadeite melting curve from 2.8 to 16.5 GPa, combined with available (or estimated) thermochemical data for jadeite crystal and Jadeite melt. The model is tested against: (1) new experimental data in the system diopside-albite from 0.1 MPa to 6.0 GPa; (2) new analyses of clinopyroxene and glass from the mid-ocean ridge basalt-pyrolite sandwich experiments of Falloon and Green (1988); and (3) published experimental data from natural and synthetic systems. The model is in good agreement with the experimental data over the entire P–Trange investigated. The only systematic deviations occur in oxidizing systems (ƒ_O2 > FMQ + 1), due to stabilization of the NaFe³ + Si₂O₆ (acmite) component in clinopyroxene, and on the diopside-albite join at 0.1 MPa. The latter can be quantitatively ascribed to nonideality of melts on the diopside-silica join in the system CaMgSi₂O₆-NaAlSi₂O₆-Si₃O₆. The activity-composition relationships adopted for crystal and melt do not require prior knowledge of the aluminum content of either phase. As D*_Na is typically within 3% relative of the weight partition coefficient (D_Na) our thermodynamic model provides a quantitative description of Na partitioning between clinopyroxene and melt in the upper mantle.