Reply to “Oxygen isotope heterogeneity of the mantle beneath the Canary Islands: a discussion of the paper of Gurenko et al.”

The comment by Day et al. (Contrib Mineral Petrol, 2012) (1) discusses the validity of the previously obtained oxygen isotope data for El Hierro and La Palma (Canary Island) olivines, (2) questions the approach by Gurenko et al. (Contrib Mineral Petrol 162:349–363, 2011) of using weakly correlated variations of δ¹⁸O_olivine values with X _px (proportion of pyroxenite-derived melt in the parental magma), and (3) provides reasons why oxygen isotope data by secondary ion mass spectrometry (SIMS) “offer sensitive means for detecting melt-crust interactions.” We respond these comments and report a new set of oxygen isotope measurements performed by SIMS and single-grain laser fluorination methods. These measurements confirm our previous data and conclusions and demonstrate the ability of the SIMS technique to analyze O isotopes in terrestrial samples with 2-sigma uncertainty better than ±0.25 ‰.     

Oxygen isotope heterogeneity of the mantle beneath the Canary Islands: a discussion of the paper of Gurenko et al.

Gurenko et al. (Contrib Mineral Petrol 162:349–363, 2011) report laser-assisted fluorination (LF) and secondary ionization mass spectrometry (SIMS) ¹⁸O/¹⁶O datasets for olivine grains from the Canary Islands of Gran Canaria, Tenerife, La Gomera, La Palma and El Hierro. As with prior studies of oxygen isotopes in Canary Island lavas (e.g. Thirlwall et al. Chem Geol 135:233–262, 1997; Day et al. Geology 37:555–558, 2009, Geochim Cosmochim Acta 74:6565–6589, 2010), these authors find variations in δ¹⁸O_ol (~4.6–6.0 ‰) beyond that measured for mantle peridotite olivine (Mattey et al. Earth Planet Sci Lett 128:231–241, 1994) and interpret this variation to reflect contributions from pyroxenite-peridotite mantle sources. Furthermore, Gurenko et al. (Contrib Mineral Petrol 162:349–363, 2011) speculate that δ¹⁸O_ol values for La Palma olivine grains measured by LF (Day et al. Geology 37:555–558, 2009, Geochim Cosmochim Acta 74:6565–6589, 2010) may be biased to low values due to the presence of altered silicate, possibly serpentine. The range in δ¹⁸O_ol values for Canary Island lavas are of importance for constraining their origin. Gurenko et al. (Contrib Mineral Petrol 162:349–363, 2011) took a subset (39 SIMS analyses from 13 grains from a single El Hierro lava; EH4) of a more extensive dataset (321 SIMS analyses from 110 grains from 16 Canary Island lavas) to suggest that δ¹⁸O_ol is weakly correlated (R ² = 0.291) with the parameter used by Gurenko et al. (Earth Planet Sci Lett 277:514–524, 2009) to describe the estimated weight fraction of pyroxenite-derived melt (Xpx). With this relationship, end-member δ¹⁸O values for HIMU-peridotite (δ¹⁸O = 5.3 ± 0.3 ‰) and depleted pyroxenite (δ¹⁸O = 5.9 ± 0.3 ‰) were defined. Although the model proposed by Gurenko et al. (Contrib Mineral Petrol 162:349–363, 2011) implicates similar pyroxenite-peridotite mantle sources to those proposed by Day et al. (Geology 37:555–558, 2009, Geochim Cosmochim Acta 74:6565–6589, 2010) and Day and Hilton (Earth Planet Sci Lett 305:226–234, 2011), there are significant differences in the predicted δ¹⁸O values of end member components in the two models. In particular, Day et al. (Geochim Cosmochim Acta 74:6565–6589, 2010) proposed a mantle source for La Palma lavas with low-δ¹⁸O (<5 ‰), rather than higher-δ¹⁸O (c.f. the HIMU-peridotite composition of Gurenko et al. in Contrib Mineral Petrol 162:349–363, 2011). Here we question the approach of using weakly correlated variations in δ¹⁸O_ol and the Xpx parameter to define mantle source oxygen isotope compositions, and provide examples of why this approach appears flawed. We also provide reasons why the LF datasets previously published for Canary Island lavas remain robust and discuss why LF and SIMS data may provide complementary information on oxygen isotope variations in ocean island basalts (OIB), despite unresolved small-scale uncertainties associated with both techniques.     

Cordierite IV: structural heterogeneity and energetics of Mg–Fe solid solutions

The local structural heterogeneity and energetic properties of 22 natural Mg–Fe cordierites, ideal formula (Mg,Fe)₂Al₄Si₅O₁₈·x(H₂O,CO₂), were investigated at length scales given by powder infrared spectroscopy (IR) and also by published electronic absorption spectra. The studied samples have iron mole fractions from X_Fe = 0.06 to 0.82 and cover most of the Mg–Fe cordierite binary. Variations in wavenumbers and line widths of the IR bands were determined as a function of composition. Most modes shift linearly to lower wavenumbers with increasing X_Fe, except those at high wavenumbers located between 900 and 1,200 cm^-1. They are vibrations that have a large internal (Si,Al)O₄ character and are not greatly affected by Mg–Fe exchange on the octahedral site. The lower wavenumber modes can be best characterized as lattice vibrations having mixed character. The systematics of the wavenumber shifts suggest small continuous variations in the "average" cordierite structure with Mg–Fe exchange and are consistent with an ideal volume of mixing, V^mix= 0, behavior (Boberski and Schreyer 1990). IR line broadening was measured using the autocorrelation function for three wavenumber regions in order to determine the range of structural heterogeneity between roughly 2 and 100 Å (0.2–10.0 nm) in the solid solution. In order to do this, an empirical correction was first made to account for the effect that small amounts of channel Na have on the phonon systematics. The results show that between 1,200 and 540 cm^-1 the line widths of the IR bands broaden slightly and linearly with increasing X_Fe. Between 350 and 125 cm^-1 nonlinear behavior was observed and it may be related to dynamic effects. These results suggest minimal excess elastic enthalpies of mixing for Mg–Fe cordierite solid solutions. Channel Na should affect measurably the thermodynamic properties of natural cordierites as evidenced by variations in the IR spectra of Na-containing samples. Occluded H₂O (Class I) and CO₂ should have little interaction with the framework and can be considered nearly "free" molecules. They should not give rise to measurable structural heterogeneity in the framework. The contribution of the crystal field stabilization energy (CFSE) of octahedral Fe²⁺ to the energetics of Mg–Fe cordierites was also investigated using published electronic absorption spectra (Khomenko et al. 2001). Two bands are observed between 8,000 and 10,500 cm^-1 and they represent electronic dd-excitations of octahedral Fe²⁺ derived from the ⁵T_2g ⁵E_g transition. They shift to higher wavenumbers with increasing X_Mg in cordierite. An analysis gives slightly asymmetric excess -CFSE across the Mg–Fe cordierite join with a maximum of about –550 J/mole towards iron-rich compositions.Editorial responsibility: J. Hoefs     

Silicon and oxygen self diffusion in enstatite polycrystals: the Milke et al. (2001) rim growth experiments revisited

Milke et al. (Contrib Mineral Petrol 142:15–26, 2001) studied the diffusion of Si, Mg and O in synthetic polycrystalline enstatite reaction rims. The reaction rims were grown at 1,000°C and 1 GPa at the contacts between forsterite grains with normal isotopic compositions and a quartz matrix extremely enriched in ¹⁸O and ²⁹Si. The enstatite reaction rim grew from the original quartz-forsterite interface in both directions producing an inner portion, which replaced forsterite and an outer portion, which replaced quartz. Here we present new support for this statement, as the two portions of the rim are clearly distinguished based on crystal orientation mapping using electron backscatter diffraction (EBSD). Milke et al. (Contrib Mineral Petrol 142:15–26, 2001) used the formalism of LeClaire (J Appl Phys 14:351–356, 1963) to derive the coefficient of silicon grain boundary diffusion from stable isotope profiles across the reaction rims. LeClaires formalism is designed for grain boundary tracer diffusion into an infinite half space with fixed geometry. A fixed geometry is an undesired limitation in the context of rim growth. We suggest an alternative model, which accounts for simultaneous layer growth and superimposed silicon and oxygen self diffusion. The effective silicon bulk diffusivity obtained from our model is approximately equal within both portions of the enstatite reaction rim: D _Si,En ^eff =1.0–4.3×10^–16 m² s^–1. The effective oxygen diffusion is relatively slow in the inner portion of the reaction rim, D _O,En ^eff =0.8–1.4×10^–16 m² s^–1, and comparatively fast, D _O,En ^eff =5.9–11.6×10^–16 m² s^–1, in its outer portion. Microstructural evidence suggests that transient porosity and small amounts of fluid were concentrated at the quartz-enstatite interface during rim growth. This leads us to suspect that the presence of an aqueous fluid accelerated oxygen diffusion in the outer portion of the reaction rim. In contrast, silica diffusion does not appear to have been affected by the spatial variation in the availability of an aqueous fluid.

---

     

Paragonite stability at 700°C in the presence of H<Subscript>2</Subscript>O–NaCl fluids: constraints on H<Subscript>2</Subscript>O activity and implications for high pressure metamorphism

We carried out reversed piston-cylinder experiments on the equilibrium paragonite = jadeite + kyanite + H₂O at 700°C, 1.5–2.5 GPa, in the presence of H₂O-NaCl fluids. Synthetic paragonite and jadeite and natural kyanite were used as starting materials. The experiments were performed on four different nominal starting compositions: X(H₂O)=1.0, 0.90, 0.75 and 0.62. Reaction direction and extent were determined from the weight change in H₂O in the capsule, as well as by optical and scanning electron microscopy (SEM). At X(H₂O)=1.0, the equilibrium lies between 2.25 and 2.30 GPa, in good agreement with the 2.30–2.45 GPa reversal of Holland (Contrib Miner Petrol 68:293–301, 1979). Lowering X(H₂O) decreases the pressure of paragonite breakdown to 2.10–2.20 GPa at X(H₂O)=0.90 and 1.85–1.90 GPa at X(H₂O)=0.75. The experiments at X(H₂O) = 0.62 yielded the assemblage albite + corundum at 1.60 GPa, and jadeite + kyanite at 1.70 GPa. This constrains the position of the isothermal paragonite–jadeite–kyanite–albite–corundum–H₂O invariant point in the system Na₂O–Al₂O₃–SiO₂–H₂O to be at 1.6–1.7 GPa and X(H₂O)~0.65±0.05. The data indicate that H₂O activity, a(H₂O), is 0.75–0.86, 0.55–0.58, and <0.42 at X(H₂O)=0.90, 0.75, and 0.62, respectively. These values approach X(H₂O)², and agree well with the a(H₂O) model of Aranovich and Newton (Contrib Miner Petrol 125:200–212, 1996). Our results demonstrate that the presence or absence of paragonite can be used to place limits on a(H₂O) in high-pressure metamorphic environments. For example, nearly pure jadeite and kyanite from a metapelite from the Sesia Lanzo Zone formed during the Eo-Alpine metamorphic event at 1.7–2.0 GPa, 550–650°C. The absence of paragonite requires a fluid with low a(H₂O) of 0.3–0.6, which could be due to the presence of saline brines.     

An ion microprobe study of anhydrous oxygen diffusion in anorthite: a comparison with hydrothermal data and some geological implications

The diffusion rate of ¹⁸O tracer atoms in anorthite (An₉₇Ab₀₃) under anhydrous conditions has been measured using SIMS techniques. The tracer source was ¹⁸O₂ 98.4% gas at 1 bar, in the temperature range 1300° C–850° C. The measured diffusion constants are D ₀=1 _–0.6 ⁺¹ ×10^–9 m²s^–1 Q=236±8 kJ mol^–1 Comparison of these values with published data for ¹⁸O diffusion in anorthite under hydrothermal conditions shows that dry oxygen diffusivities are orders of magnitude lower than equivalent wet values at similar temperatures. The effect of these differences on oxygen isotope equilibration during cooling is discussed.     

Experimental study of a low-alkali tholeiite at 1–5 kbar: optimal condition for the crystallization of high-An plagioclase in hydrous arc tholeiite

We conducted melting experiments on a low-alkali tholeiite (SiO₂ ~52 wt%, MgO ~6.5 wt%, CaO/Na₂O~4.4, Al₂O₃/SiO₂ ~0.33) under both H₂O-undersaturated and H₂O-saturated conditions to investigate the effect of H₂O on the Ca–Na partitioning between plagioclase and melt. Experiments were performed in the temperature and pressure ranges of 1,000–1,300°C and 1–5 kbar, respectively, with varying H₂O contents of 0–12wt%. Redox condition was 0–2 log unit above NNO (nickel–nickel oxide) buffer. Temperature-bulk H₂O diagrams for the low-alkali tholeiite are constructed at 1, 2, and 5 kbar, and compositions of near-liquidus plagioclase and coexisting melt are determined. To exclude the effect of melt composition (CaO/Na₂O and Al₂O₃/SiO₂ ratios) on plagioclase composition and to reveal the effect of H₂O on An (=100×Ca/(Ca+Na)) content and (=(Ca/Na)_pl/(Ca/Na)_melt), we focused on the composition of near-liquidus plagioclases which crystallized from melts with nearly constant CaO/Na₂O and Al₂O₃/SiO₂ ratios. Our experimental results show that, at each experimental pressure, An content of the near-liquidus plagioclase and the K_D^Ca-Na almost linearly increases as H₂O content in melt increases. Each of the An content and the variations in a low-alkali tholeiitic system (CaO/Na₂O~4.0–4.5, Al₂O₃/SiO₂ ~0.27–0.33) can be described by one equation using temperature, pressure, and melt H₂O content as parameters. An content and of liquidus plagioclase increases with increasing melt H₂O and with decreasing pressure, elucidating that nearly H₂O-saturated conditions of 2–3 kbar is optimal for the crystallization of the most An-rich plagioclase (>An₈₈). We suggest this pressure condition of 2–3 kbar, corresponding to depth of 7–11 km, plays an important role for the origin of An-rich plagioclase in H₂O-rich low-alkali tholeiite. At pressures more than ca. 4 kbar, crystallization of liquidus Ca-rich clinopyroxene decreases the CaO/Na₂O ratio of liquid, thus prohibiting the crystallization of high-An plagioclase from hydrous tholeiite.     

Potential protonation sites in the Al<Subscript>2</Subscript>SiO<Subscript>5</Subscript> polymorphs based on polarized FTIR spectroscopy and properties of the electron density distribution

Potential protonation sites for, kyanite, sillimanite, and andalusite, located in a mapping of the (3, −3) critical points displayed by their L(r) = −∇²ρ(r) distributions, are compared with polarized single-crystal FTIR spectra of kyanite and sillimanite determined earlier and with andalusite measured in this study. For andalusite, seven peaks were observed when the electric vector, E, is parallel to [100]: four intense ones at 3,440, 3,460, 3,526, and 3,597 cm⁻¹ and three weaker ones at 3,480, 3,520, and 3,653 cm⁻¹. Six peaks, three intense ones at 3,440, 3,460, and 3,526 cm⁻¹ and three weaker ones at 3,480, 3,520, and 3,653 cm⁻¹ when E parallels [010]. No peaks were observed when E is parallel to [001]. The concentration of water in andalusite varies between 110 and 168 ppm by weight % H₂O. Polarized FTIR spectra indicate that the OH vector is parallel to (001) in andalusite and sillimanite and in kyanite. Examination of the L(r) (3, −3) critical points in comparison with the polarized FTIR indicates that H prefers to bond to the oxygen atoms O1 and O2 in andalusite and O2 and O4 in sillimanite which correspond to the underbonded oxygen atoms and those with the largest L(r) maxima. In kyanite, comparison of the FTIR spectrum and the critical points indicates that H will preferentially bond to the two 4-coordinated O2 and O6 atoms.     

Investigation of discharge and groundwater contamination characteristics around the Karasu Spring, Muş Area, Turkey

The study area is 35 km east of the city center of Mu, Turkey. The rock units in the area include Paleozoic Bitlis Metamorphites, Cenozoic flysch and Solhan Volcanities and Pliocene formations. A long-term discharge change was monitored in the Karasu spring. The amount of water discharged by the spring during June 3 and October 11, 2002 was 7.18×10⁶ m³ while its storage capacity during May 26 and November 11, 2002 was 4.42×10⁶ m³.NH₄ – N, NO₂ – N, NO₃ – N, total PO₄ – P, SO₄ and total organic carbon (TOC) contamination were determined in the recharge area of the spring and other well and spring waters in the vicinity. NH₄ – N, NO₂ – N, NO₃ – N, total PO₄ – P, SO₄ and total organic carbon contamination around the Budak, ivbai, Gölbai, Sazlikbai and Altinova villages are well above the drinking water standards of the Turkish Standard Institute (TSE) and World Health Organization (WHO) (1984). Similarly, water samples of K₄, K₅ and K₆ branches of the Karasu spring and S₁, S₂, S₃ and S₄ points in the swamp area yielded concentrations much more than limit values.     

Diamondiferous peridotite xenoliths from the Argyle (AK1) lamproite pipe,Western Australia

This paper describes a suite of peridotite xenoliths. some carrying diamonds at high grades, from the richly diamondiferous early Proterozoic (1180 Ma) Argyle (AK1) lamproite pipe, in northwestern Australia. The peridotites are mostly coarse garnet lherzolites but also include garnet harzburgite, chromite — garnet peridotite, a garnet wehrlite, and an altered spinel peridotite with extremely Cr-rich chromite. In all cases the garnet has been replaced by a kelyphite-like, symplectic intergrowth of Alrich pyroxenes, Al-spinel and secondary silicates. The peridotites have refractory compositions characterized by high Mg/(Mg+Fe) and depletion in lithophile elements (Al₂O₃ and CaO < 1%, Na₂O0.03%) and high field strength cations such as Ti, Zr, Y, and Yb. Olivines have high Mg/(Mg+Fe) (Mg ^91–93 ) and, like olivine inclusions in diamonds from the Argyle pipe, contain detectable amounts of Cr₂O₃ (0.03%–0.07%) but have very low CaO contents (typically 0.04%–0.05%). Enstatites (Mg ^92–94 ) have comparatively high Cr₂O₃ (0.2%–0.45%) and Na₂O (up to 0.18%) but very low Al₂O₃ contents (0.5%–0.7%). Diopsides (Mg ^92–94 , Ca/(Ca+Mg+Fe)=0.37–0.43) are Cr-rich (0.7%–1.9% Cr₂O₃) and have low Al₂O₃ (0.7%–2.2%) and Na₂O (0.5%–1.6%) contents. Many have high K₂O contents, typically 0.1%–0.4% but up to 1.3% K₂O in one xenolith. The chromite coexisting with former garnet is Mg-and Cr-rich [Mg/(Mg+Fe²⁺)=0.68–0.72, Cr/(Cr+Al)=0.72–0.79] whereas chromite in the spinel peridotite is even more Cr-rich (65% Cr₂O₃, Cr/(Cr+Al)=0.85, resembling inclusions in diamond. One highly serpentinized former garnet peridotite contains a Cr-rich (up to 13% Cr₂O₃) titanate resembling armalcolite but containing significant K₂O (1%–2.5%), CaO (0.6%–2.2%), ZrO₂ (0.1%–0.8%), SrO (0.1%–0.3%), and BaO (up to 0.58%): this appears to have formed as an overprint of the primary mineralogy. Temperatures and pressures estimated from coexisting pyroxenes and reconstructed garnet compositions indicate that the garnet lherzolites equilibrated at 1140°–1290° C and 5.0–5.9 GPa (160–190 km depth), within the stability field of diamond. Oxygen fugacties within the diamond forming environment are estimated from spinel-bearing assemblages to be reducing, with f _O2 between MW and IW. The presence of significant K in the diopsides from the peridotite xenoliths and in diopsides from heavy mineral concentrate from the Argyle pipe implies metasomatic enrichment of the subcontinental lithosphere within the diamond stability field. The P-T conditions estimated for the Argyle peridotites demonstrate that diamondiferous lamproite magmas incorporate mantle xenoliths from similar depths to kimberlites in cratonic settings, and imply that Proterozoic cratonized orogenic belts can have lithospheric roots of comparable thickness to beneath Archaean cratons. These roots lie at the base of the lithosphere within the stability field of diamond. The xenoliths, the calcic nature of chrome pyropes from heavy mineral concentrate, and the diamond inclusion assemblage indicate that the lighosphere beneath the Western Australian lamproites is mostly depleted lherozolite rather than the harzburgite commonly found beneath Archaean cratons. Nevertheless, the dominance of eclogitic paragenesis inclusions in Argyle diamonds indicates a significant proportion of diamondiferous eclogite is also present. The form, mineral inclusion assemblage, and the C-isotopic composition of diamonds in the peridotite xenoliths suggest that disaggregated diamondiferous peridotites are the source of the planar octahedral diamonds which constitute a minor component of the Argyle production. These diamonds are believed to have formed from mantle carbon in reduced, refractory peridotite (Iherzolite-harzburgite) in contrast to the predominant strongly ¹³C-depleted eclogitic suite diamonds which contain a recycled crustal carbon component. The source region of the lamproites has undergone long-term (2 Ga) enrichment in incompatible elements.     

The temperature dependence of water solubility in enstatite

The solubility of water in pure enstatite was measured on samples synthesized under water-saturated conditions at 15 kbar and temperatures ranging from 700 to 1,100°C. Polarized FTIR measurements on millimetre-sized, clear crystals showed that water solubility increases strongly with temperature, from 101 ppm by weight at 700°C to 269 ppm by weight at 1,100°C. The position and shape of the infrared bands hardly changes with temperature, with one notable exception: a band close to 3,380 cm^–1 is present in samples synthesized between 700 and 1,000°C, while this band is absent from samples synthesized at 1,100°C. This effect appears to be very reproducible and points towards a slight change in the crystal structure of enstatite between 1,000 and 1,100°C at 15 kbar. The water solubility data of this study as well as those of Rauch and Keppler (Contrib Mineral Petrol 143:525–536, 2002) can be reproduced by the equation where K is water solubility, is water fugacity, A is 0.01354 ppm/bar, V^solid=12.1 cm³/mol is the volume change of enstatite during incorporation of water, and H^{1 bar}=-4,563 J/mol is the reaction enthalpy at 1 bar. This equation predicts the following behaviour of water solubility in enstatite as a function of pressure and temperature: (1) water solubility increases with pressure up to a maximum around 80 kbar; (2) water solubility decreases with temperature at 1 bar; and (3) water solubility increases with temperature between 10 and 100 kbar. If the observed temperature dependence for enstatite were representative for other upper mantle minerals as well, it would have the following implications: (1) Lateral temperature gradients in the upper mantle could cause major variations in water contents at the same depth; in particular, hot mantle plumes may scavenge water from the surrounding shallow upper mantle. (2) The scavenging of water by hot plumes could be a major factor in increasing the mobility of plumes. (3) The predicted temperature dependence of water solubility at the base of the upper mantle may allow plumes to bypass the transition zone water filter postulated by Bercovici and Karato (Nature 425:39–44, 2003).     

The role of hydrogen bonding in the thermal expansion and dehydration of brushite,di-calcium phosphate dihydrate

The unit-cell and atomic parameters of perdeuterated brushite have been extracted from Rietveld analysis of neutron powder diffraction data within the temperature range 4.2 to 470 K. The thermal expansion of brushite is anisotropic, with the largest expansion along the b axis due principally to the effect of the O(1)···D(4) and O(3)···D(2) hydrogen bonds. Expansion along the c axis, influenced by the Ow1···D(5) interwater hydrogen bond, is also large. The high temperature limits for the expansion coefficients for the unit-cell edges a, b and c are 9.7(5) × 10^–6, 3.82(9) × 10^–5 and 5.54(5) × 10^–5 K^–1, respectively, and for the cell volume it is 9.7(1) × 10^–5 K^–1. The angle displays oscillatory variation, and empirical data analysis results in = 1.28(3) × 10^–6sin(0.0105 T) K^–1, within this temperature range. The evolution of the thermal expansion tensor of brushite has been calculated between 50 T 400 K. At 300 K the magnitudes of the principal axes are ₁₁ = 50(6) × 10^{–6 K–1}, ₂₂ = 26.7(7) × 10^–6 K^–1 and ₃₃ = 7.0(5) × 10^–6 K^–1. The intermediate axis, ₂₂, is parallel to b, and using IRE convention for the tensor orthonormal basis, the axes ₁₁ and ₃₃ have directions equal to (–0.228, 0, –0.974) and (–0.974, 0, 0.228) respectively. Under the conditions of these experiments, the onset of dehydration occurred at temperatures above 400 K. Bond valence analysis combined with assessments of the thermal evolution of the bonding within brushite suggests that dehydration is precipitated through instabilities in the chemical environment of the second water molecule.     

Crystal structure refinement of hendricksite,A Zn- and Mn-rich trioctahedral potassium mica: A contribution to the crystal chemistry of zinc-bearing minerals

Summary The crystal structure of hendricksite, a trioctahedral mica of biotite type, characterized by high Zn²⁺ and Mn²⁺ contents has been refined by least square methods. The structural formula is: (K_0.89Na_0.10Ba_0.04)(Mg_1.57Zn_0.54Mn _0.40 ²⁺ Fe _0.25 ²⁺ Al_0.07Ti_0.07Cr_0.01)(Si_2.92Al_1.08)O₁₀ (OH)₂. The space group isC2/m and the cell parameters are:a=5.340(2) Å,b=9.524(2) Å,c=10.235(3) Å, =100.07(2)^o, the cell volume isV=497.98 ų. The final unweightedR=0.072. Average cation-anion distances in polyhedra are: T–O=1.659 Å, M(1)–O=2.093 Å, M(2)–O=2.088 Å, A–O_long=3.316 Å and A–O_short=3.004 Å; A is the alkaline cation. The rotation angle of tetrahedra is =6.7°. The analysis of electron densities, of the dimensions and distorsions of polyhedra shows that Zn²⁺ is exclusively in octahedral sites; there is no order between six-fold coordinated cations. A comparison between the structural features of hendricksite and those of the two main end-members of biotites, phlogopite and annite, is presented.The effect of the strong covalence of Zn–O bonds is particularly visible on the dimensions and orientations of the thermal ellipsoids of octahedral sites M(1) and M(2) which contain zinc. In all the published structures of trioctahedral micas, the ellipsoids of cationic sites are uniaxial positive, elongated parallel toc ^*. In hendricksite, this is observed only for the two zinc-free sites (T and A; in the octahedra M(1) and M(2), which contain zinc, the ellipsoids are approximately uniaxial negative, flattened parallel toa, which is a unique situation.Zinc which habitually favours the tetrahedral coordinations with oxygen, enters the octahedra only, i.e. the chemically anisotropic sites, in hendricksite. The strong polarizability of Zn²⁺ is proposed to explain this behaviour.An examination of the behaviour of Zn²⁺ in other compounds shows that this situation is general, zinc favours chemically anisotropic sites and specially those adjacent to OH or H₂O.

---

Affinement de la structure cristalline de la hendricksite, mica trioctaédrique potassique riche en Zn et Mn; une contribution à la connaissance cristallochimique des minéraux zincifères

Résumé On a affiné par moindres carrés la structure de la hendricksite, mica trioctaédrique de type biotite, caractérisé par une teneur élevée en Zn²⁺ et Mn²⁺. La formule structurale de ce mica est: (K_0m89Na_0,10Ba_0,04)(Mg_1,57Zn_0,54Mn _0,40 ²⁺ Fe _0,25 ²⁺ Al_0,07Ti_0,07Cr_0,01)(Si_2,92Al_1,08)O₁₀(OH)₂. Le groupe spatial estC2/m et les paramètres de la maille:a=5,340(2) Å,b=9,254(2) Å,c=10,235(3) Å, =100,07(2)°; le volume de la maille estV=497,98 ų. Le résidu final non-pondéré estR=0,072. Les distances cation-anion moyennes dans les polyèdres sont les suivantes: T–O=1,659 Å, M(1)–O=2,093 Å, M(2)–O=2,088 Å, A–O_long=3,316 Å et A–O_court=3,004 Å où A désigne le cation alcalin. L'angle de rotation tétraédrique, =6,7°, est très semblable à celui de la phlogopite. L'analyse des densités électroniques, des dimensions et distorsions des polyèdres montre que Zn²⁺ est exclusivement en coordinance octaédrique et qu'il n'y a pas d'ordre entre les cations hexacoordonnés. On présente une comparaison des caractères structuraux de la hendricksite avec ceux des deux principaux pôles des biotites, la phlogopite et l'annite.L'effet de la forte covalence de la liaison Zn–O est particulièrement visible sur les dimensions et orientations des ellipsoides d'agitation thermique des deux sites octaédriques, sites zincifères. Dans toutes les structures de micas trioctaédriques publiées, les ellipsoides des sites cationiques sont uniaxes positifs, allongés parallèlement àc ^*, ce qui s' observe effectivement dans les deux sites non-zincifères (T et A) de la hendricksite, par contre, dans les octaèdres M(1) et M(2), qui contiennent le zinc, les ellipsoides sont approximativement uniaxes négatifs, applatis parallèlement àa, ce qui est une situation unique.Le zinc, qui se fixe généralement en sites tétraédriques dans les structures de type oxyde, occupe les sites octaédriques, c'est-à-dire les sites chimiquement anisotropes dans la hendricksite. La forte polarisabilité de Zn²⁺ est proposée pour expliquer ce comportement.Un examen du comportement de Zn²⁺ dans d'autres phases montre que cette situation est tout à fait générale, le zinc privilégiant les sites chimiquement anisotropes et en particulier ceux adjacents à OH où H₂O.

     

Heat capacity,entropy and phase equilibria of stishovite

The low-temperature heat capacity (C _P) of stishovite (SiO₂) synthesized with a multi-anvil device was measured over the range of 5–303 K using the heat capacity option of a physical properties measurement system (PPMS) and around ambient temperature using a differential scanning calorimeter (DSC). The entropy of stishovite at standard temperature and pressure calculated from DSC-corrected PPMS data is 24.94 J mol⁻¹ K⁻¹, which is considerably smaller (by 2.86 J mol⁻¹ K⁻¹) than that determined from adiabatic calorimetry (Holm et al. in Geochimica et Cosmochimica Acta 31:2289–2307, 1967) and about 4% larger than the recently reported value (Akaogi et al. in Am Mineral 96:1325–1330, 2011). The coesite–stishovite phase transition boundary calculated using the newly determined entropy value of stishovite agrees reasonably well with the previous experimental results by Zhang et al. (Phys Chem Miner 23:1–10, 1996). The calculated phase boundary of kyanite decomposition reaction is most comparable with the experimental study by Irifune et al. (Earth Planet Sci Lett 77:245–256, 1995) at low temperatures around 1,400 K, and the calculated slope in this temperature range is mostly consistent with that determined by in situ X-ray diffraction experiments (Ono et al. in Am Mineral 92:1624–1629, 2007).     

Hydroxyl in mantle olivine xenocrysts from the Udachnaya kimberlite pipe

The incorporation of hydrogen in mantle olivine xenocrysts from the Udachnaya kimberlite pipe was investigated by Fourier-transform infrared spectroscopy and secondary ion mass spectrometry (SIMS). IR spectra were collected in the OH stretching region on oriented single crystals using a conventional IR source at ambient conditions and in situ at temperatures down to −180°C as well as with IR synchrotron radiation. The IR spectra of the samples are complex containing more than 20 strongly polarized OH bands in the range 3,730–3,330 cm⁻¹. Bands at high energies (3,730–3,670 cm⁻¹) were assigned to inclusions of serpentine, talc and the 10 Å phase. All other bands are believed to be intrinsic to olivine. The corresponding point defects are (a) associated with vacant Si sites (3,607 cm⁻¹ E || a, 3,597 E || a, 3,571 cm⁻¹ E || c, 3,567 E || c, and 3,556 E || b), and (b) with vacant M1 sites (most of the bands polarized parallel to a). From the pleochroic behavior and position of the OH bands associated with the vacant M1 sites, we propose two types of hydrogen—one bonded to O1 and another to O2, so that both OH vectors are strongly aligned parallel to a. The O2–H groups may be responsible for the OH bands at higher wavenumbers than those for the O1–H groups. The multiplicity of the corresponding OH bands in the spectra can be explained by different chemical environments and by slightly different distortions of the M1 sites in these high-pressure olivines. Four samples were investigated by SIMS. The calculated integral molar absorption coefficient using the IR and SIMS results of 37,500±5,000 L mol H₂O cm⁻² is within the uncertainties slightly higher than the value determined by Bell et al. (J Geophys Res 108(B2):2105–2113, 2003) (28,450±1,830 L mol H₂O cm⁻²). The reason for the difference is the different distributions of the absorption intensity of the spectra of both studies (mean wavenumber 3,548 vs. 3,570 cm⁻¹). Olivine samples with a mean wavenumber of about 3,548 cm⁻¹ should be quantified with the absorption coefficient as determined in this study; those containing more bands at higher wavenumber (mean wavenumber 3,570 cm⁻¹) should be quantified using the value determined by Bell et al. (J Geophys Res 108(B2):2105–2113, 2003).

Oxygen isotope fractionation between calcite and tremolite: an experimental study

Oxygen isotope partitioning between calcite and tremolite was experimentally calibrated in the presence of small amounts of a supercritical CO₂–H₂O fluid at temperatures from 520 to 680° C and pressures from 3 to 10 kbar. The experiments were carried out within the stability field of the calcite-tremolite assemblage based on phase equilibrium relationships in the system CaO–MgO–SiO₂–CO₂–H₂O, so that decomposition of calcite and tremolite was avoided under the experimental conditions. Appropriate proportions of carbon dioxide to water were used to meet this requirement. Large weight ratios of mineral to fluid were employed in order to make the isotopic exchange between calcite and tremolite in the presence of a fluid close to that without fluid. The data processing method for isotopic exchange in a three-phase system has been applied to extrapolate partial equilibrium data to equilibrium values. The determined fractionation factors between calcite (Cc) and tremolite (Tr) are expressed as:10³1n _Cc-Tr=3.80 × 10⁶/T ²-1.67By combining the present data with the experimental calibrations of Clayton et al. (1989) on the calcite-quartz system, we obtain the fractionation for the quartztremolite system: 10³1n _Qz-Tr=4.18 × 10⁶/T ²-1.67Our experimental calibrations are in good agreement with the theoretical calculations of Hoffbauer et al. (1994) and the empirical estimates of Bottinga and Javoy (1975) based on isotopic data from naturall assemblages. At 700 C good agreement also exists between our experimental data and theoretical values calculated by Zheng (1993b). With decreasing temperature, however, an increasing difference between these data appears.Retrograde isotopic reequilibration by oxygen diffusion may be common for amphibole relative to diopside in metamorphic rocks. However, isotopic equilibrium in amphibole can be preserved in cases of rapid cooling.     

An infrared and <Superscript>1</Superscript>H MAS NMR investigation of strong hydrogen bonding in ussingite,Na<Subscript>2</Subscript>AlSi<Subscript>3</Subscript>O<Subscript>8</Subscript>(OH)

The mineral ussingite, Na₂AlSi₃O₈(OH), an interrupted tectosilicate, has strong hydrogen bonding between OH and the other nonbridging oxygen atom in the structure. Infrared spectra contain a strongly polarized, very broad OH-stretching band with an ill-defined maximum between 1500 and 1800 cm^–1, and a possible OH librational bending mode at 1295 cm^–1. The IR spectra confirm the orientation of the OH vector within the triclinic unit cell as determined from X-ray refinement (Rossi et al. 1974). There are three distinct bands in the ¹H NMR spectrum of ussingite: a predominant band at 13.5 ppm (TMS) representing 90% of the structural hydrogen, a second band at 15.9 ppm corresponding to 8% of the protons, and a third band at 11.0 ppm accounting for the remaining 2% of structural hydrogen. From the correlation between hydrogen bond length and ¹H NMR chemical shift (Sternberg and Brunner 1994), the predominant hydrogen bond length (H...O) was calculated to be 1.49 Å, in comparison to the hydrogen bond length determined from X-ray refinement (1.54 Å). The population of protons at 15.9 ppm is consistent with 5–8% Al–Si disorder. Although the ussingite crystal structure and composition are similar to those of low albite, the bonding environment of OH in low albite and other feldspars, as characterized through IR and ¹H NMR, is fundamentally different from the strong hydrogen bonding found in ussingite.     

Thermodynamic properties,low-temperature heat-capacity anomalies,and single-crystal X-ray refinement of hydronium jarosite, (H<Subscript>3</Subscript>O)Fe<Subscript>3</Subscript>(SO<Subscript>4</Subscript>)<Subscript>2</Subscript>(OH)<Subscript>6</Subscript>

Crystals of hydronium jarosite were synthesized by hydrothermal treatment of Fe(III)–SO₄ solutions. Single-crystal XRD refinement with R₁=0.0232 for the unique observed reflections (|F_o| > 4_F) and wR₂=0.0451 for all data gave a=7.3559(8) Å, c=17.019(3) Å, V^o=160.11(4) cm³, and fractional positions for all atoms except the H in the H₃O groups. The chemical composition of this sample is described by the formula (H₃O)_0.91Fe_2.91(SO₄)₂[(OH)_5.64(H₂O)_0.18]. The enthalpy of formation (H^o_f) is –3694.5 ± 4.6 kJ mol^–1, calculated from acid (5.0 N HCl) solution calorimetry data for hydronium jarosite, -FeOOH, MgO, H₂O, and -MgSO₄. The entropy at standard temperature and pressure (S^o) is 438.9±0.7 J mol^–1 K^–1, calculated from adiabatic and semi-adiabatic calorimetry data. The heat capacity (C_p) data between 273 and 400 K were fitted to a Maier-Kelley polynomial C_p(T in K)=280.6 + 0.6149T–3199700T^–2. The Gibbs free energy of formation is –3162.2 ± 4.6 kJ mol^–1. Speciation and activity calculations for Fe(III)–SO₄ solutions show that these new thermodynamic data reproduce the results of solubility experiments with hydronium jarosite. A spin-glass freezing transition was manifested as a broad anomaly in the C_p data, and as a broad maximum in the zero-field-cooled magnetic susceptibility data at 16.5 K. Another anomaly in C_p, below 0.7 K, has been tentatively attributed to spin cluster tunneling. A set of thermodynamic values for an ideal composition end member (H₃O)Fe₃(SO₄)₂(OH)₆ was estimated: G^o_f= –3226.4 ± 4.6 kJ mol^–1, H^o_f=–3770.2 ± 4.6 kJ mol^–1, S^o=448.2 ± 0.7 J mol^–1 K^–1, C_p (T in K)=287.2 + 0.6281T–3286000T^–2 (between 273 and 400 K).     

Mid-Pleistocene lavas from the Seguam volcanic center,central Aleutian arc: closed-system fractional crystallization of a basalt to rhyodacite eruptive suite

In contrast to adjacent volcanic centers of the modern central Aleutian arc, Seguam Island developed on strongly extended arc crust. K-Ar dates indicate that mid-Pleistocene, late-Pleistocene, and Holocene eruptive phases constitute Seguam. This study focuses on the petrology of the mid-Pleistocene, 1.07–07 Ma, Turf Point Formation (TPF) which is dominated by an unusual suite of porphyritic basalt and basaltic andesite lavas with subordinate phenocryst-poor andesite to rhyodacite lavas. Increasing whole-rock FeO^*/MgO from basalt to dacite, the anhydrous Plag+Ol+Cpx±Opx±Mt phenocryst assemblage, groundmass pigeonite, and the reaction Ol+Liq=Opx preserved in the mafic lavas indicate a tholeiitic affinity. Thermometry and comparison to published phase equilibria suggests that most TPF basalts crystallized Plag+Ol+Cpx±Mt at 1160°C between about 3–5 kb (±1–2% H₂O), andesites crystallized Plag+Cpx+Opx±Mt at 1000°C between 3–4 kb with 3–5% H₂O, and dacites crystallized Plag +Cpx±Opx±Mt at 1000°C between 1–2 kb with 2–3% H₂O. All lavas crystallized at f _{o 2} close to the NNO buffer. Mineral compositions and textures indicate equilibrium crystallization of the evolved lavas; petrographic evidence of open-system mixing or assimilation is rare. MgO, CaO, Al₂O₃, Cr, Ni, and Sr abundances decrease and K₂O, Na₂O, Rb, Ba, Zr, and Pb increase with increasing SiO₂ (50–71%). LREE enrichment [(Ce/Yb)_n=1.7±0.2] characterizes most TPF lavas; total REE contents increase and Eu anomalies become more negative with increasing SiO₂. Relative to other Aleutian volcanic centers, TPF basalts and basaltic andesites have lower K₂O, Na₂O, TiO₂, Rb, Ba, Sr, Zr, Y, and LREE abundances. ⁸⁷Sr/⁸⁶Sr ratios (0.70361–0.70375) and ratios of ²⁰⁶Pb/²⁰⁴Pb (18.88–18.97), ²⁰⁷Pb/²⁰⁴Pb (15.58–15.62), ²⁰⁸Pb/²⁰⁴Pb (38.46–38.55) are the highest measured for any suite of lavas in the oceanic portion of the Aleutian arc. Conversely, Nd values (+5.8 to+6.7) are among the lowest from the Aleutians. Sr, Nd, and Pb ratios are virtually constant from basalt through rhyodacite, whereas detectable isotopic heterogenity is observed at most other Aleutian volcanic centers. Major and trace element, REE, and Sr, Nd, and Pb isotopic compositions are consistent with the basaltic andesitic, andesitic, dacitic, and rhyodacitic liquids evolving from TPF basaltic magma via closed-system fractional crystallization alone. Fractionation models suggest that removal of 80 wt% cumulate (61% Plag, 17% Cpx, 12% Opx, 7% Ol, and 3% Mt) can produce 20 wt% rhyodacitic residual liquid per unit mass of parental basaltic liquid. Petrologic and physical constraints favor segregation of small batches of basalt from a larger mid-crustal reservoir trapped below a low-density upper crustal lid. In these small magma batches, the degree of cooling, crystallization, and fractionation are functions of the initial mass of basaltic magma segregated, the thermal state of the upper crust, and the magnitude of extension. Tholeiitic magmas erupted at Seguam evolved by substantially different mechanisms than did calc-alkaline lavas erupted at the adjacent volcanic centers of Kanaga and Adak on unextended arc crust. These variable differentiation mechanisms and liquid lines of descent reflect contrasting thermal and mechanical conditions imposed by the different tectonic environments in which these centers grew. At Seguam, intra-arc extension promoted eruption of voluminous basalt and its differentiates, unmodified by interaction with lower crustal or upper mantle wallrocks.     

Basalt contamination by continental crust: Some experiments and models

Chemical interaction between molten basalt and felsic minerals of the continental crust (quartz, K-feldspar, and oligoclase) was examined in static and dynamic experiments at 1,200°–1,400° C. Under circumstances of continuous stirring at 1,400°, -quartz dissolves in tholeiite melt at a rate of 3.3×10^–6 g/s per cm² of contact area; at 1,300°, the solution rate is 1.5×10^–6 g/cm{cm2}s. The feldspars are molten at the experimental conditions, and interact with contacting basalt melt by diffusion in the liquid state. This is a complex process characterized by rapid initial diffusion of alkalies to establish a distribution between felsic melt and basalt similar to that observed in cases of actual two-liquid equilibrium (both alkalies reach concentrations in the felsic melt 1.5–3 times those in the basalt). Alkali diffusion may be uphill or downhill, depending on which direction of net flux is required to produce a two-liquid type distribution. Once this distribution is attained, subsequent diffusion of all melt species is slow and apparently limited by the diffusivity of SiO₂, which is 10^–9-l0^–10 cm²/s at 1,200° C. Interdiffusion experiments involving molten basalt and synthetic granite confirm the behavior illustrated by the feldspar/basalt results, and give similar SiO₂ diffusivities.The solution rates and interdiffusion data can be used to model basalt contamination processes likely to occur in the continental crust. For the restricted case of solid quartzitic xenoliths, the uptake of SiO₂ in a well-mixed basalt magma is quite fast: appreciable SiO₂ contamination may occur over exposure times of only days to years. If basalt magma induces local melting of crustal rocks, the assimilation process becomes one of liquid-state interdiffusion. In this case, the varying diffusivities of ions and their differing preferences for silicic relative to basaltic melts can produce marked selective contamination effects. Selective contamination of ascending basaltic magmas is particularly likely in the case of K₂O, which may be introduced in substantial amounts even when other elements remain unaffected. The Na₂O content of mantle-derived magmas is buffered against contamination by crustal materials, and K₂O is buffered against further increases once it reaches a level of 1–1.5 wt.%.