IR calibrations for water determination in olivine,r-GeO<Subscript>2</Subscript>, and SiO<Subscript>2</Subscript> polymorphs

Mineral-specific IR absorption coefficients were calculated for natural and synthetic olivine, SiO₂ polymorphs, and GeO₂ with specific isolated OH point defects using quantitative data from independent techniques such as proton–proton scattering, confocal Raman spectroscopy, and secondary ion mass spectrometry. Moreover, we present a routine to detect OH traces in anisotropic minerals using Raman spectroscopy combined with the “Comparator Technique”. In case of olivine and the SiO₂ system, it turns out that the magnitude of ε for one structure is independent of the type of OH point defect and therewith the peak position (quartz ε = 89,000 ± 15,000  \textl \textmol_{\textH2\textO}^-1 \textcm^-2\text{l}\,\text{mol}_{{\text{H}_2}\text{O}}^{-1}\,\text{cm}^{-2}), but it varies as a function of structure (coesite ε = 214,000 ± 14,000  \textl \textmol_{\textH2\textO}^-1 \textcm^-2\text{l}\,\text{mol}_{{\text{H}_2}\text{O}}^{-1}\,\text{cm}^{-2}; stishovite ε = 485,000 ± 109,000  \textl \textmol_{\textH2\textO}^-1 \textcm^-2\text{l}\,\text{mol}_{{\text{H}_2}\text{O}}^{-1}\,\text{cm}^{-2}). Evaluation of data from this study confirms that not using mineral-specific IR calibrations for the OH quantification in nominally anhydrous minerals leads to inaccurate estimations of OH concentrations, which constitute the basis for modeling the Earth’s deep water cycle.     

The H<Subscript>2</Subscript>O solubility of alkali basaltic melts: an experimental study

Experiments were conducted to determine the water solubility of alkali basalts from Etna, Stromboli and Vesuvius volcanoes, Italy. The basaltic melts were equilibrated at 1,200°C with pure water, under oxidized conditions, and at pressures ranging from 163 to 3,842 bars. Our results show that at pressures above 1 kbar, alkali basalts dissolve more water than typical mid-ocean ridge basalts (MORB). Combination of our data with those from previous studies allows the following simple empirical model for the water solubility of basalts of varying alkalinity and fO₂ to be derived: \textH ₂ \textO( \textwt% ) = \text H ₂ \textO_\textMORB ( \textwt% ) + ( 5.84 ×10 ^{- 5} *\textP - 2.29 ×10 ^{- 2} ) ×( \textNa₂ \textO + \textK₂ \textO )( \textwt% ) + 4.67 ×10 ^{- 2} ×\Updelta \textNNO - 2.29 ×10 ^{- 1} {\text{H}}_{ 2} {\text{O}}\left( {{\text{wt}}\% } \right) = {\text{ H}}_{ 2} {\text{O}}_{\text{MORB}} \left( {{\text{wt}}\% } \right) + \left( {5.84 \times 10^{ - 5} *{\text{P}} - 2.29 \times 10^{ - 2} } \right) \times \left( {{\text{Na}}_{2} {\text{O}} + {\text{K}}_{2} {\text{O}}} \right)\left( {{\text{wt}}\% } \right) + 4.67 \times 10^{ - 2} \times \Updelta {\text{NNO}} - 2.29 \times 10^{ - 1} where H₂O_MORB is the water solubility at the calculated P, using the model of Dixon et al. (1995). This equation reproduces the existing database on water solubilities in basaltic melts to within 5%. Interpretation of the speciation data in the context of the glass transition theory shows that water speciation in basalt melts is severely modified during quench. At magmatic temperatures, more than 90% of dissolved water forms hydroxyl groups at all water contents, whilst in natural or synthetic glasses, the amount of molecular water is much larger. A regular solution model with an explicit temperature dependence reproduces well-observed water species. Derivation of the partial molar volume of molecular water using standard thermodynamic considerations yields values close to previous findings if room temperature water species are used. When high temperature species proportions are used, a negative partial molar volume is obtained for molecular water. Calculation of the partial molar volume of total water using H₂O solubility data on basaltic melts at pressures above 1 kbar yields a value of 19 cm³/mol in reasonable agreement with estimates obtained from density measurements.     

An Equation of State for Hypersaline Water in Great Salt Lake,Utah, USA

Great Salt Lake (GSL) is one of the largest and most saline lakes in the world. In order to accurately model limnological processes in GSL, hydrodynamic calculations require the precise estimation of water density (ρ) under a variety of environmental conditions. An equation of state was developed with water samples collected from GSL to estimate density as a function of salinity and water temperature. The ρ of water samples from the south arm of GSL was measured as a function of temperature ranging from 278 to 323 degrees Kelvin (^oK) and conductivity salinities ranging from 23 to 182 g L⁻¹ using an Anton Paar density meter. These results have been used to develop the following equation of state for GSL (σ = ± 0.32 kg m⁻³):

H<Subscript>2</Subscript>O storage capacity of olivine and low-Ca pyroxene from 10 to 13 GPa: consequences for dehydration melting above the transition zone

The onset of hydrous partial melting in the mantle above the transition zone is dictated by the H₂O storage capacity of peridotite, which is defined as the maximum concentration that the solid assemblage can store at P and T without stabilizing a hydrous fluid or melt. H₂O storage capacities of minerals in simple systems do not adequately constrain the peridotite water storage capacity because simpler systems do not account for enhanced hydrous melt stability and reduced H₂O activity facilitated by the additional components of multiply saturated peridotite. In this study, we determine peridotite-saturated olivine and pyroxene water storage capacities at 10–13 GPa and 1,350–1,450°C by employing layered experiments, in which the bottom ~2/3 of the capsule consists of hydrated KLB-1 oxide analog peridotite and the top ~1/3 of the capsule is a nearly monomineralic layer of hydrated Mg# 89.6 olivine. This method facilitates the growth of ~200-μm olivine crystals, as well as accessory low-Ca pyroxenes up to ~50 μm in diameter. The presence of small amounts of hydrous melt ensures that crystalline phases have maximal H₂O contents possible, while in equilibrium with the full peridotite assemblage (melt + ol + pyx + gt). At 12 GPa, olivine and pyroxene water storage capacities decrease from ~1,000 to 650 ppm, and ~1,400 to 1,100 ppm, respectively, as temperature increases from 1,350 to 1,450°C. Combining our results with those from a companion study at 5–8 GPa (Ardia et al., in prep.) at 1,450°C, the olivine water storage capacity increases linearly with increasing pressure and is defined by the relation C_{\textH2 \textO}^\textolivine ( \textppm ) = 57.6( ±16 ) ×P( \textGPa ) - 169( ±18 ). C_{{{\text{H}}_{2} {\text{O}}}}^{\text{olivine}} \left( {\text{ppm}} \right) = 57.6\left( { \pm 16} \right) \times P\left( {\text{GPa}} \right) - 169\left( { \pm 18} \right). Adjustment of this trend for small increases in temperature along the mantle geotherm, combined with experimental determinations of D_{\textH2 \textO}^{\textpyx/olivine} D_{{{\text{H}}_{2} {\text{O}}}}^{\text{pyx/olivine}} from this study and estimates of D_{\textH2 \textO}^{\textgt/\textolivine} D_{{{\text{H}}_{2} {\text{O}}}}^{{{\text{gt}}/{\text{olivine}}}} , allows for estimation of peridotite H₂O storage capacity, which is 440 ± 200 ppm at 400 km. This suggests that MORB source upper mantle, which contains 50–200 ppm bulk H₂O, is not wet enough to incite a global melt layer above the 410-km discontinuity. However, OIB source mantle and residues of subducted slabs, which contain 300–1,000 ppm bulk H₂O, can exceed the peridotite H₂O storage capacity and incite localized hydrous partial melting in the deep upper mantle. Experimentally determined values of D_{\textH2 \textO}^{\textpyx/\textolivine} D_{{{\text{H}}_{2} {\text{O}}}}^{{{\text{pyx}}/{\text{olivine}}}} at 10–13 GPa have a narrow range of 1.35 ± 0.13, meaning that olivine is probably the most important host of H₂O in the deep upper mantle. The increase in hydration of olivine with depth in the upper mantle may have significant influence on viscosity and other transport properties.     

Complexation of Pb(II) by Chloride Ions in Aqueous Solutions

Lead chloride formation constants at 25°C were derived from analysis of previous spectrophotometrically generated observations of lead speciation in a variety of aqueous solutions (HClO₄–HCl and NaCl–NaClO₄ mixtures, and solutions of MgCl₂ and CaCl₂). Specific interaction theory analysis of these formation constants produced coherent estimates of (a) PbCl⁺, \textPbCl₂⁰ {\text{PbCl}}_{2}^{0} , and PbCl₃⁻ formation constants at zero ionic strength, and (b) well-defined depictions of the dependence of these formation constants on ionic strength. Accompanying examination of a recent IUPAC critical assessment of lead formation constants, in conjunction with the spectrophotometrically generated formation constants presented in this study, revealed significant differences among various subsets of the IUPAC critically selected data. It was found that these differences could be substantially reduced through reanalysis of the formation constant data of one of the subsets. The resulting revised lead chloride formation constants are in good agreement with the formation constants derived from the earlier spectrophotometrically generated data. Combining these data sets provides an improved characterization of lead chloride complexation over a wide range of ionic strengths:

Paragenesis of unusual Fe-cordierite (sekaninaite)-bearing paralava and clinker from the Kuznetsk coal basin,Siberia, Russia

Sekaninaite (X_Fe > 0.5)-bearing paralava and clinker are the products of ancient combustion metamorphism in the western part of the Kuznetsk coal basin, Siberia. The combustion metamorphic rocks typically occur as clinker beds and breccias consisting of vitrified sandstone–siltstone clinker fragments cemented by paralava, resulting from hanging-wall collapse above burning coal seams and quenching. Sekaninaite–Fe-cordierite (X_Fe = 95–45) is associated with tridymite, fayalite, magnetite, ± clinoferrosilite and ±mullite in paralava and with tridymite and mullite in clinker. Unmelted grains of detrital quartz occur in both rocks (<3 vol% in paralavas and up to 30 vol% in some clinkers). Compositionally variable siliceous, K-rich peraluminous glass is <30% in paralavas and up to 85% in clinkers. The paralavas resulted from extensive fusion of sandstone–siltstone (clinker), and sideritic/Fe-hydroxide material contained within them, with the proportion of clastic sediments ≫ ferruginous component. Calculated dry liquidus temperatures of the paralavas are 1,120–1,050°C and 920–1,050°C for clinkers, with calculated viscosities at liquidus temperatures of 10^1.6–7.0 and 10^7.0–9.8 Pa s, respectively. Dry liquidus temperatures of glass compositions range between 920 and 1,120°C (paralava) and 920–960°C (clinker), and viscosities at these temperatures are 10^9.7–5.5 and 10^8.8–9.7 Pa s, respectively. Compared with worldwide occurrences of cordierite–sekaninaite in pyrometamorphic rocks, sekaninaite occurs in rocks with X_Fe (mol% FeO/(FeO + MgO)) > 0.8; sekaninaite and Fe-cordierite occur in rocks with X_Fe 0.6–0.8, and cordierite (X_Fe < 0.5) is restricted to rocks with X_Fe < 0.6. The crystal-chemical formula of an anhydrous sekaninaite based on the refined structure is | \textK_0.02 |(\textFe_1.54^{2 +} \textMg_0.40 \textMn_0.06 )_{\Upsigma 2.00}^M [(\textAl_1.98 \textFe_0.02^{2 +} \textSi_1.00 )_{\Upsigma 3.00}^T1 (\textSi_3.94 \textAl_2.04 \textFe_0.02^{2 +} )_{\Upsigma 6.00}^T2 \textO₁₈ ]. \left| {{\text{K}}_{0.02} } \right|({\text{Fe}}_{1.54}^{2 + } {\text{Mg}}_{0.40} {\text{Mn}}_{0.06} )_{\Upsigma 2.00}^{M} [({\text{Al}}_{1.98} {\text{Fe}}_{0.02}^{2 + } {\text{Si}}_{1.00} )_{\Upsigma 3.00}^{T1} ({\text{Si}}_{3.94} {\text{Al}}_{2.04} {\text{Fe}}_{0.02}^{2 + } )_{\Upsigma 6.00}^{T2} {\text{O}}_{18} ].     

The notion of “distinguishability” between bulk elastic parameters on the basis of the Gibbs deformation energy

The present work aims in discussing a principle that distinguishes between elastic parameters sets, $ \{ \Upphi \} \equiv \{ K_{0} , \, K^{\prime}, \, V_{0} ,\ldots\} The present work aims in discussing a principle that distinguishes between elastic parameters sets, { \Upphi } o { K₀ ,  K^￠,  V₀ ,?} \{ \Upphi \} \equiv \{ K_{0} , \, K^{\prime}, \, V_{0} ,\ldots\} , on the basis of an energetic criterion: once a reference set, { \Upphi_R } \{ \Upphi_{R} \} , is given, another one can be fixed, { \Upphi _min } \left\{ {\Upphi_{ \min } } \right\} , so that they are as close as possible to each other, but yield non-equivalent deformation energy curves \Updelta G({ \Upphi } )_\textdeform \Updelta G(\{ \Upphi \} )_{\text{deform}} , i.e. they give \Updelta G({ \Upphi_R } )_\textdeform \Updelta G(\{ \Upphi_{R} \} )_{\text{deform}} and \Updelta G({ \Upphi _min } )_\textdeform \Updelta G(\{ \Upphi_{ \min } \} )_{\text{deform}} such that | \Updelta G({ \Upphi _min } )_\textdeform - \Updelta G({ \Upphi_R } )_\textdeform | 3 1×s[\Updelta G_\textdeform ]. \left| {\Updelta G(\{ \Upphi_{ \min } \} )_{\text{deform}} - \Updelta G(\{ \Upphi_{R} \} )_{\text{deform}} } \right| \ge 1\times \sigma [\Updelta G_{\text{deform}} ]. ΔG _deform, calculated using the equation of state (EoS), and its uncertainty σ[ΔG _deform], obtained by a propagation of the errors affecting { \Upphi } \{ \Upphi \} are crucial to fix which mineral assemblage forms at P–T conditions and allow one to assess the reliability of such a prediction. We explore some properties related to the principle introduced, using the average values of the elastic parameters found in literature and related uncertainties for di-octahedral mica, olivine, garnet and clinopyroxene. Two elementary applications are briefly discussed: the effect of refining V ₀ in fitting EoSs to P–V experimental data, in the case of garnet and omphacite, and the phengite 3T–2M ₁ relative stability, controlled by pressure.     

A 57Fe M?ssbauer spectroscopic study of sogdianite: an example of a symmetric electric field gradient around Fe3+

Sogdianite, a double-ring silicate of composition ( \textZr_{0. 7 6} \textTi_{0. 3 8}^{4 +} \textFe_{0. 7 3}^{3 +} \textAl_0.13 )_{\Upsigma = 2} ( \square _{1. 1 5} \textNa_{0. 8 5} )_{\Upsigma = 2} \textK[\textLi ₃ \textSi _{1 2} \textO ₃₀ ] ( {\text{Zr}}_{0. 7 6} {\text{Ti}}_{0. 3 8}^{4 + } {\text{Fe}}_{0. 7 3}^{3 + } {\text{Al}}_{0.13} )_{\Upsigma = 2} \left( {\square_{ 1. 1 5} {\text{Na}}_{0. 8 5} } \right)_{\Upsigma = 2} {\text{K}}[{\text{Li}}_{ 3} {\text{Si}}_{ 1 2} {\text{O}}_{ 30} ] from Dara-i-Pioz, Tadjikistan, was studied by the combined application of ⁵⁷Fe M?ssbauer spectroscopy and electronic structure calculations. The M?ssbauer spectrum confirms published microprobe and X-ray single-crystal diffraction results that indicate that Fe³⁺ is located at the octahedral A-site and that no Fe²⁺ is present. Both the measured and calculated quadrupole splitting, ΔE _Q, for Fe³⁺ are virtually 0 mm s⁻¹. Such a value is unusually small for a silicate and it is the same as the ΔE _Q value for Fe³⁺ in structurally related sugilite. This result is traced back to the nearly regular octahedral coordination geometry corresponding to a very symmetric electric field gradient around Fe³⁺. A crystal chemical interpretation for the regular octahedral geometry and the resulting low ΔE _Q value for Fe³⁺ in the M?ssbauer spectrum of sogdianite is that structural strain is largely “taken up” by weak Li–O bonds permitting highly distorted LiO₄ tetrahedra. Weak Li–O bonding allows the edge-shared more strongly bonded Fe³⁺O₆ octahedra to remain regular in geometry. This may be a typical property for all double-ring silicates with tetrahedrally coordinated Li.     

Thermal expansion of plagioclase feldspars

The volume thermal expansion coefficient and the anisotropy of thermal expansion were determined for nine natural feldspars with compositions, in terms of albite (NaAlSi₃O₈, Ab) and anorthite (CaAl₂Si₂O₈, An), of Ab₁₀₀, An₂₇Ab₇₃, An₃₅Ab₆₅, An₄₆Ab₅₄, An₆₀Ab₄₀, An₇₈Ab₂₂, An₈₉Ab₁₁, An₉₆Ab₄ and An₁₀₀ by high resolution powder diffraction with a synchrotron radiation source. Unit-cell parameters were determined from 124 powder patterns of each sample, collected over the temperature range 298–935 K. The volume thermal expansion coefficient of the samples determined by a linear fit of V/V ₀ = α(T − T ₀) varies with composition (X _An in mol %) as:

Structural relaxation and crystal field stabilization in Cr3+-containing oxides and silicates

The effect of crystal structure relaxation in oxygen-based Cr³⁺-containing minerals on the crystal field stabilization energy (CFSE) is considered. It is shown that the dependence of \textCFSE_{\textCr ³⁺} {\text{CFSE}}_{{{\text{Cr}}^{ 3+ } }} , which is found from optical absorption spectra, on the average interatomic distances is described by the power function with a negative exponent c \mathord

TitaniQ under pressure: the effect of pressure and temperature on the solubility of Ti in quartz

Quartz and rutile were synthesized from silica-saturated aqueous fluids between 5 and 20 kbar and from 700 to 940°C in a piston-cylinder apparatus to explore the potential pressure effect on Ti solubility in quartz. A systematic decrease in Ti-in-quartz solubility occurs between 5 and 20 kbar. Titanium K-edge X-ray absorption near-edge structure (XANES) measurements demonstrate that Ti⁴⁺ substitutes for Si⁴⁺ on fourfold tetrahedral sites in quartz at all conditions studied. Molecular dynamic simulations support XANES measurements and demonstrate that Ti incorporation onto fourfold sites is favored over interstitial solubility mechanisms. To account for the P–T dependence of Ti-in-quartz solubility, a least-squares method was used to fit Ti concentrations in quartz from all experiments to the simple expression

HTP21/c–C2/c phase transition and kinetics of Fe2+–Mg order–disorder of an Fe-poor pigeonite: implications for the cooling history of ureilites

A natural Ca-poor pigeonite (Wo₆En₇₆Fs₁₈) from the ureilite meteorite sample PCA82506-3, free of exsolved augite, was studied by in situ high-temperature single-crystal X-ray diffraction. The sample, monoclinic P2₁/c, was annealed up to 1,093°C to induce a phase transition from P2₁/c to C2/c symmetry. The variation with increasing temperature of the lattice parameters and of the intensity of the b-type reflections (h + k = 2n + 1, present only in the P2₁/c phase) showed a displacive phase transition P2₁/c to C2/c at a transition temperature T _Tr = 944°C, first order in character. The Fe–Mg exchange kinetics was studied by ex situ single-crystal X-ray diffraction in a range of temperatures between the closure temperature of the Fe–Mg exchange reaction and the transition temperature. Isothermal disordering annealing experiments, using the IW buffer, were performed on three crystals at 790, 840 and 865°C. Linear regression of ln k _D versus 1/T yielded the following equation: ln k_\textD = - 3717( ±416)/T(K) + 1.290( ±0.378);    (R² = 0.988) \ln \,k_{\text{D}} = - 3717( \pm 416)/T(K) + 1.290( \pm 0.378);\quad (R^{2} = 0.988) . The closure temperature (T _c) calculated using this equation was ∼740(±30)°C. Analysis of the kinetic data carried out taking into account the e.s.d.'s of the atomic fractions used to define the Fe–Mg degree of order, performed according to Mueller’s model, allowed us to retrieve the disordering rate constants C ₀ K _dis⁺ for all three temperatures yielding the following Arrhenius relation: ln( C₀ K_\textdis ⁺ ) = ln K₀ - Q/(RT) = 20.99( ±3.74) - 26406( ±4165)/T(K);    (R² = 0.988) \ln \left( {C_{0} K_{\text{dis}}^{ + } } \right) = \ln \,K_{0} - Q/(RT) = 20.99( \pm 3.74) - 26406( \pm 4165)/T(K);\quad (R^{2} = 0.988) . An activation energy of 52.5(±4) kcal/mol for the Fe–Mg exchange process was obtained. The above relation was used to calculate the following Arrhenius relation modified as a function of X _Fe (in the range of X _Fe = 0.20–0.50): ln( C₀ K_\textdis ⁺ ) = (21.185 - 1.47X_\textFe ) - \frac(27267 - 4170X_\textFe )T(K) \ln \left( {C_{0} K_{\text{dis}}^{ + } } \right) = (21.185 - 1.47X_{\text{Fe}} ) - {\frac{{(27267 - 4170X_{\text{Fe}} )}}{T(K)}} . The cooling time constant, η = 6 × 10⁻¹ K⁻¹ year⁻¹ calculated on the PCA82506-3 sample, provided a cooling rate of the order of 1°C/min consistent with the extremely fast late cooling history of the ureilite parent body after impact excavation.     

Iron in monticellite as an oxygen barometer for kimberlite magmas

Monticellite is a common magmatic mineral in the groundmass of kimberlites. A new oxygen barometer for kimberlite magmas is calibrated based on the Fe content of monticellite, CaMgSiO₄, in equilibrium with kimberlite liquids in experiments at 100 kPa from 1,230 to 1,350°C and at logfO₂ from NNO-4.1 to NNO+5.3 (where NNO is the nickel–nickel oxide buffer). The XFe_Mtc/XFe_liq was found to decrease with increasing fO₂, consistent with only Fe²⁺ entering the monticellite structure. Although the XFe-in-monticellite varies with temperature and composition, these dependencies are small compared to that with fO₂. The experimental data were fitted by weighted least square regression to the following relationship: \Updelta \textNNO = \frac{ log[ 0.858( ±0.021)\fracX_\textFe^\textLiq X_\textFe^\textMtc ] - 0.139( ±0.022) }0.193( ±0.004) \Updelta {\text{NNO}} = \frac{{\left\{ {\log \left[ {0.858( \pm 0.021)\frac{{X_{\text{Fe}}^{\text{Liq}} }}{{X_{\text{Fe}}^{\text{Mtc}} }}} \right] - 0.139( \pm 0.022)} \right\}}}{0.193( \pm 0.004)} where ΔNNO is the fO₂ relative to that of the Nickel-bunsenite (NNO) buffer and XFe_liq/XFe_Mtc is the ratio of mole fraction of Fe in liquid and Fe-in-monticellite (uncertainties at 2σ). The application of this oxygen barometer to natural kimberlites from both the literature and our own investigations, assuming the bulk rock FeO is that of their liquid FeO, revealed a range in fO₂ from NNO-3.5 to NNO+1.7. A range of Mg/(Mg + Fe²⁺) (Mg#) for kimberlite melts of 0.46–0.88 was derived from the application of the experimentally determined monticellite-liquid Kd Fe²⁺–Mg to natural monticellites. The range in Mg# is broader and less ultramafic than previous estimates of kimberlites, suggesting an evolution under a wide range of petrologic conditions.     

High-pressure electronic absorption spectroscopy of natural and synthetic Cr<Superscript>3+</Superscript>-bearing clinopyroxenes

Comparison of polarized optical absorption spectra of natural Ca-rich diopsides and synthetic NaCrSi2O6 and LiCrSi2O6 clinopyroxenes, evidences as vivid similarities, as noticeable differences. The similarities reflect the fact that in all cases Cr3+ enters the small octahedral M1-site of the clinopyroxene structure. The differences are due to some iron content in the natural samples causing broad intense near infrared bands of electronic spin-allowed dd transitions of Fe2+(M1, M2) and intervalence Fe2+/Fe3+ charge-transfer transition, and by different symmetry and different local crystal fields strength of Cr3+ in the crystal structures. The positions of the spin-allowed bands of Cr3+, especially of the low energy one caused by the electronic 4 A 2g → 2 T 1g transition, are found to be in accordance with mean M1–O distances. The local relaxation parameter ε calculated for limCr 3+ → 0 from the spectra and interatomic á Cr - O ñ \left\langle {Cr - O} \right\rangle and á Mg - O ñ \left\langle {Mg - O} \right\rangle distances yields a very high value, 0.96, indicating that in the clinopyroxene structure the local lattice relaxation around the “guest” ion, Cr3+, deviates greatly from the “diffraction” value, ε = 0, than in any other known Cr3+-bearing systems studied so far. Under pressure the spin-allowed bands of Cr3+ shift to higher energies and decrease in intensity quite in accordance with the crystal field theoretical expectations, while the spin-forbidden absorption lines remain practically unshifted, but also undergo a strong weakening. There is no evident dependence of the Racah parameter B of Cr3+ reflecting the covalence of the oxygen-chromium bond under pressure: within the uncertainty of determination it may be regarded as practically constant. The values of CrO6 octahedral modulus, k\textpoly\textloc k_{\text{poly}}^{\text{loc}} , derived from high-pressure spectra of natural chromium diopside and synthetic NaCrSi2O6 kosmochlor are very close, ~203 and ~196 GPa, respectively, being, however, nearly twice higher than that of MgO6 octahedron in diopside, 105(4) GPa, obtained by Thompson and Downs (2008). Such a strong stiffening of the structural octahedron, i.e. twice higher value of k\textCr3 + \textloc k_{{{\text{Cr}}^{3 + } }}^{\text{loc}} comparing with that of k\textMg2 + \textloc k_{{{\text{Mg}}^{2 + } }}^{\text{loc}} , may be caused by simultaneous substitution of Ca2+ by larger Na+ in the neighboring M2 sites at so-called jadeite-coupled substitution Mg2+ + Ca2+ → Cr3+ + Na+. It is also remarkable that the values of CrO6 octahedral modulus of NaCrSi2O6 kosmochlor obtained here are nearly twice larger than that of 90(16) GPa, evaluated by high-pressure X-ray structural refinement by Origlieri et al. (2003). Taking into account that the overall compressibility of the clinopyroxene structure should mainly be due to the compressibility of M1- and M2-sites, our k\textCr3 + \textloc k_{{{\text{Cr}}^{3 + } }}^{\text{loc}} -value, ~196 GPa, looks much more consistent with the bulk modulus value, 134(1) GPa.     

Practical Estimates of Tensile Strength and Hoek–Brown Strength Parameter m i of Brittle Rocks

By applying the Griffith stress criterion of brittle failure, one can find that the uniaxial compressive strength (σ_c) of rocks is eight times the value of the uniaxial tensile strength (σ_t). The Griffith strength ratio is smaller than what is normally measured for rocks, even with the consideration of crack closure. The reason is that Griffith’s theories address only the initiation of failure. Under tensile conditions, the crack propagation is unstable so that the tensile crack propagation stress (σ_cd)_t and the peak tensile strength σ_t are almost identical to the tensile crack initiation stress (σ_ci)_t. On the other hand, the crack growth after crack initiation is stable under a predominantly compressive condition. Additional loading is required in compression to bring the stress from the crack initiation stress σ_ci to the peak strength σ_c. It is proposed to estimate the tensile strength of strong brittle rocks from the strength ratio of R = \fracs_\textc | s_\textt | = 8\fracs_\textc s_\textci . R = {\frac{{\sigma_{\text{c}} }}{{\left| {\sigma_{\text{t}} } \right|}}} = 8{\frac{{\sigma_{\text{c}} }}{{\sigma_{\text{ci}} }}}. The term \fracs_\textc s_\textci {\frac{{\sigma_{\text{c}} }}{{\sigma_{\text{ci}} }}} accounts for the difference of crack growth or propagation in tension and compression in uniaxial compression tests. \fracs_c s_ci {\frac{{\sigma_{c} }}{{\sigma_{ci} }}} depends on rock heterogeneity and is larger for coarse grained rocks than for fine grained rocks. σ_ci can be obtained from volumetric strain measurement or acoustic emission (AE) monitoring. With the strength ratio R determined, the tensile strength can be indirectly obtained from | s_\textt | = \fracs_\textc R = \fracs_\textci 8. \left| {\sigma_{\text{t}} } \right| = {\frac{{\sigma_{\text{c}} }}{R}} = {\frac{{\sigma_{\text{ci}} }}{8}}. It is found that the predicted tensile strengths using this method are in good agreement with test data. Finally, a practical estimate of the Hoek–Brown strength parameter m _i is presented and a bi-segmental or multi-segmental representation of the Hoek–Brown strength envelope is suggested for some brittle rocks. In this fashion, the rock strength parameters like σ_t and m _i, which require specialty tests such as direct tensile (or Brazilian) and triaxial compression tests for their determination, can be reasonably estimated from uniaxial compression tests.     

Five new E′ centers and their 29Si hyperfine structures in electron-irradiated α-quartz

Single-crystal electron paramagnetic resonance (EPR) spectra of fast-electron-irradiated quartz, after annealing at 120 and 200°C, reveal five new E′ type centers, herein labeled E ₅^￠ , E ₆^￠ , E ₇^￠ , E ₈^￠ , \textand E ₉^￠ E_{ 5}^{\prime } ,\,E_{ 6}^{\prime } ,\,E_{ 7}^{\prime } ,\,E_{ 8}^{\prime } ,\,{\text{and}}\,E_{ 9}^{\prime } . Centers E ₅^￠ , E ₇^￠ , \textand E ₉^￠ E_{ 5}^{\prime } ,\,E_{ 7}^{\prime } ,\,{\text{and}}\,E_{ 9}^{\prime } are characterized by the orientations of the unique principal g and A(²⁹Si) axes close to a short Si–O bond direction, hence representing new variants of the well-established E ₁^￠ E_{ 1}^{\prime } center. Centers E ₆^￠ E_{ 6}^{\prime } and E ₈^￠ E_{ 8}^{\prime } have the orientations of the unique principal g and A(²⁹Si) axes approximately along a long Si–O bond direction, similar to the E ₂^￠ E_{ 2}^{\prime } centers. Therefore, these new E′ type centers apparently arise from the removal of different oxygen atoms and represent variable local distortions around the oxygen vacancies.     

Experimental and computational study of trace element distribution between orthopyroxene and anhydrous silicate melt: substitution mechanisms and the effect of iron

Although orthopyroxene (Opx) is present during a wide range of magmatic differentiation processes in the terrestrial and lunar mantle, its effect on melt trace element contents is not well quantified. We present results of a combined experimental and computational study of trace element partitioning between Opx and anhydrous silicate melts. Experiments were performed in air at atmospheric pressure and temperatures ranging from 1,326 to 1,420°C in the system CaO–MgO–Al₂O₃–SiO₂ and subsystem CaO–MgO–SiO₂. We provide experimental partition coefficients for a wide range of trace elements (large ion lithophile: Li, Be, B, K, Rb, Sr, Cs, Ba, Th, U; rare earth elements, REE: La, Ce, Nd, Sm, Y, Yb, Lu; high field strength: Zr, Nb, Hf, Ta, Ti; transition metals: Sc, V, Cr, Co) for use in petrogenetic modelling. REE partition coefficients increase from $ D_{\text{La}}^{{\text{Opx}} {\hbox{-}} {\text{melt}}} \sim 0.0005 Although orthopyroxene (Opx) is present during a wide range of magmatic differentiation processes in the terrestrial and lunar mantle, its effect on melt trace element contents is not well quantified. We present results of a combined experimental and computational study of trace element partitioning between Opx and anhydrous silicate melts. Experiments were performed in air at atmospheric pressure and temperatures ranging from 1,326 to 1,420°C in the system CaO–MgO–Al₂O₃–SiO₂ and subsystem CaO–MgO–SiO₂. We provide experimental partition coefficients for a wide range of trace elements (large ion lithophile: Li, Be, B, K, Rb, Sr, Cs, Ba, Th, U; rare earth elements, REE: La, Ce, Nd, Sm, Y, Yb, Lu; high field strength: Zr, Nb, Hf, Ta, Ti; transition metals: Sc, V, Cr, Co) for use in petrogenetic modelling. REE partition coefficients increase from $ D_{\text{La}}^{{\text{Opx}} {\hbox{-}} {\text{melt}}} \sim 0.0005 $ D_{\text{La}}^{{\text{Opx}} {\hbox{-}} {\text{melt}}} \sim 0.0005 to $ D_{\text{Lu}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.109 $ D_{\text{Lu}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.109 , D values for highly charged elements vary from $ D_{\text{Th}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.0026 $ D_{\text{Th}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.0026 through $ D_{\text{Nb}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.0033 $ D_{\text{Nb}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.0033 and $ D_{\text{U}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.0066 $ D_{\text{U}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.0066 to $ D_{\text{Ti}}^{{\text{Opx}} {\hbox{-}} {\text{melt}}} \sim 0.058 $ D_{\text{Ti}}^{{\text{Opx}} {\hbox{-}} {\text{melt}}} \sim 0.058 , and are all virtually independent of temperature. Cr and Co are the only compatible trace elements at the studied conditions. To elucidate charge-balancing mechanisms for incorporation of REE into Opx and to assess the possible influence of Fe on Opx-melt partitioning, we compare our experimental results with computer simulations. In these simulations, we examine major and minor trace element incorporation into the end-members enstatite (Mg₂Si₂O₆) and ferrosilite (Fe₂Si₂O₆). Calculated solution energies show that R²⁺ cations are more soluble in Opx than R³⁺ cations of similar size, consistent with experimental partitioning data. In addition, simulations show charge balancing of R³⁺ cations by coupled substitution with Li⁺ on the M1 site that is energetically favoured over coupled substitution involving Al–Si exchange on the tetrahedrally coordinated site. We derived best-fit values for ideal ionic radii r ₀, maximum partition coefficients D ₀, and apparent Young’s moduli E for substitutions onto the Opx M1 and M2 sites. Experimental r ₀ values for R³⁺ substitutions are 0.66–0.67 ? for M1 and 0.82–0.87 ? for M2. Simulations for enstatite result in r ₀ = 0.71–0.73 ? for M1 and ~0.79–0.87 ? for M2. Ferrosilite r ₀ values are systematically larger by ~0.05 ? for both M1 and M2. The latter is opposite to experimental literature data, which appear to show a slight decrease in $ r_{0}^{{{\text{M}}2}} $ r_{0}^{{{\text{M}}2}} in the presence of Fe. Additional systematic studies in Fe-bearing systems are required to resolve this inconsistency and to develop predictive Opx-melt partitioning models for use in terrestrial and lunar magmatic differentiation models.     

The carbon dioxide solubility in alkali basalts: an experimental study

Experiments were conducted to determine CO₂ solubilities in alkali basalts from Vesuvius, Etna and Stromboli volcanoes. The basaltic melts were equilibrated with nearly pure CO₂ at 1,200°C under oxidizing conditions and at pressures ranging from 269 to 2,060 bars. CO₂ solubility was determined by FTIR measurements. The results show that alkalis have a strong effect on the CO₂ solubility and confirm and refine the relationship between the compositional parameter Π devised by Dixon (Am Mineral 82:368–378, 1997) and the CO₂ solubility. A general thermodynamic model for CO₂ solubility in basaltic melts is defined for pressures up to 2 kbars. Based on the assumption that O²⁻ and CO₃²⁻ mix ideally, we have:

Structure and properties of loaded silica contacts during pressure solution: impedance spectroscopy measurements under hydrothermal conditions

In order to investigate directly the structure and properties of grain boundaries in silicate materials undergoing pressure solution, in situ measurements of these properties are required. We report electrical impedance spectroscopy measurements, performed, under hydrothermal conditions, on individual glass–glass and glass-quartz contacts undergoing pressure solution. Resulting estimates of the average grain boundary diffusivity product ( Z = Dd_\textav C^* Z = D\delta_{\text{av}} C^{*} ) for silica transport and of the average grain boundary fluid film thickness ( d_\textav \delta_{\text{av}} ) fall in the ranges 6.3 ± 1.4 × 10⁻¹⁸ m³ s⁻¹ and 350 ± 210 nm, respectively. However, the average values for Z and d_\textav \delta_{\text{av}} obtained were likely influenced by cracking and irregular dissolution of the dissolving contact surfaces, rather than representing uniformly wetted grain boundary properties. Post-mortem SEM observations indicate that the contact surfaces were internally rough. Taken together, our data support the notion that during pressure solution of quartz, grain boundary diffusion is rapid, and interface processes (dissolution and precipitation) are more likely to be rate-limiting than diffusion.     

Infrared spectroscopic characterization of dehydration and accompanying phase transition behaviors in NAT-topology zeolites

Relative humidity ( P_{\textH 2 \textO} P_{{{\text{H}}_{ 2} {\text{O}}}} , partial pressure of water)-dependent dehydration and accompanying phase transitions in NAT-topology zeolites (natrolite, scolecite, and mesolite) were studied under controlled temperature and known P_{\textH 2 \textO} P_{{{\text{H}}_{ 2} {\text{O}}}} conditions by in situ diffuse-reflectance infrared Fourier transform spectroscopy and parallel X-ray powder diffraction. Dehydration was characterized by the disappearance of internal H₂O vibrational modes. The loss of H₂O molecules caused a sequence of structural transitions in which the host framework transformation path was coupled primarily via the thermal motion of guest Na⁺/Ca²⁺ cations and H₂O molecules. The observation of different interactions of H₂O molecules and Na⁺/Ca²⁺ cations with host aluminosilicate frameworks under high- and low- P_{\textH 2 \textO} P_{{{\text{H}}_{ 2} {\text{O}}}} conditions indicated the development of different local strain fields, arising from cation–H₂O interactions in NAT-type channels. These strain fields influence the Si–O/Al–O bond strength and tilting angles within and between tetrahedra as the dehydration temperature is approached. The newly observed infrared bands (at 2,139 cm⁻¹ in natrolite, 2,276 cm⁻¹ in scolecite, and 2,176 and 2,259 cm⁻¹ in mesolite) result from strong cation–H₂O–Al–Si framework interactions in NAT-type channels, and these bands can be used to evaluate the energetic evolution of Na⁺/Ca²⁺ cations before and after phase transitions, especially for scolecite and mesolite. The 2,176 and 2,259 cm⁻¹ absorption bands in mesolite also appear to be related to Na⁺/Ca²⁺ order–disorder that occur when mesolite loses its Ow4 H₂O molecules.