A search for evolutionary period variations of Cepheids using the Harvard plate stacks: FF Aql

Photographic plates of the Harvard Observatory stacks have been used to obtain 1261 brightness estimates for the low-amplitude Cepheid FF Aql (P = 4.47 days). Combined with published visual, photoelectric, and CCD observations, these data enable the construction of an O-C diagram spanning 122 years. The resulting O-C diagram is parabola-shaped, enabling the determination of quadratic brightness elements for the first time. The evolutionary rate of increase of the period, dP/dt = 0.072 ± 0.011 s/year, is in agreement with theoretical computations for the third crossing of the instability strip. The available data, reduced using the technique of by Eddington and Plakidis, reveal the presence of small, random period fluctuations, ? = 0.0061^d ± 0.0044^d, that do not distort the evolutionary trend of the O-C residuals.     

A search for evolutionary period variations of Cepheids using the harvard plate stacks: II Carinae

We obtained 380 photographic brightness estimates for the long-period classical Cepheid II Car (P = 64.4 days) from the plate stacks of the Harvard College Observatory. Combined with the published photoelectric and CCD observations, these data enable the construction of an O-C diagram covering 120 years. The O-C diagram is parabolic, making it possible for the first time to determine the quadratic light-curve elements and calculate the evolutionary rate of increase of the period, dP/dt = 719(±15) s/year (log(Ṗ/P) = −4.64), in agreement with theoretical computations for the third crossing of the instability strip. The available data, reduced using the method described by Eddington and Plakidis, demonstrate the presence of small random period fluctuations that do not distort the evolutionary changes of the O-C residuals.     

The stability of sphene; experimental redetermination and geologic implications

Hydrothermal experiments on the reaction: rutile + calcite + quartz = sphene + CO₂, suggest equilibrium at: P_fluid = 2 kbar, 500 ± 5°C, X_CO2 = 0.5 ± 0.006; P_fluid = 3.45 kbar, 535 ± 10°C, X_CO2 = 0.5 ± 0.006; P_fluid = 5 kbar, <580°C, X_CO2 = 0.5 ± 0.006; P_fluid = 5 kbar, <635°C, X_CO2 = 0.97 ± 0.035.The resultant stability field of sphene is much larger than that suggested by Schuiling and Vink (1967). Calculation of the reaction in P-T space directly from free energy data indicates a discrepancy in Δ_r of about 1.5 kcal/gfw when compared to our experimental data. In addition, P-T extrapolation of the reaction using thermochemical data disagrees with the constraints of the two experimental brackets, indicating possible errors in the entropy of one or more solid phases and/or in CO₂ free energy data.The combination of these experimental results with published data on reactions involving additional minerals allows the recognition of sphene and/or rutile-bearing assemblages which can be used as indicators of P-T-X_CO2 conditions at various grades of metamorphism. In particular, mineral assemblages useful in delimiting P-T-X_CO2 conditions are: muscovite-sphene (H₂O-rich fluids), diopside-rutile-quartz (high temperatures or a CO₂-rich fluid at low temperatures), epidote-rutile-quartz (low temperatures, H₂O-rich conditions), prehnite-rutile (very low temperatures, very H₂O-rich conditions).     

High-pressure X-ray study of LiCrSi2O6 clinopyroxene and the general compressibility trends for Li-clinopyroxenes

High-pressure single-crystal X-ray diffraction measurements of synthetic LiCrSi₂O₆ clinopyroxene (with space group P2₁/c) were performed in a diamond-anvil cell up to 7.970 GPa. No phase transition has been observed within the pressure range investigated, but the elastic behavior at lower pressures (up to ~2.5 GPa) is affected by an anomalous softening due to the proximity of the phase transition to the HT-C2/c phase at 330 K and at ambient pressure. A third-order Birch–Murnaghan equation of state fitted to the compression data above 2.5 GPa yields a bulk modulus K _T0 = 93(2) GPa and its first derivative K′ = 8.8(6). The structural data measured up to 7.970 GPa confirm that the space group P2₁/c is maintained throughout the whole pressure range investigated. The atomic parameters, obtained from the integrated diffraction intensities, suggest that the Li coordination polyhedron changes its coordination number from 5 to 6 at 6–7 GPa by means of the approach of the bridging O atom, related to the increased kinking of the B tetrahedral chain. Furthermore, at higher pressures, the structural evolution of LiCrSi₂O₆ provides evidence in the variation of kinking angles and bond lengths of a potential phase transition above 8 GPa to the HP-C2/c space group. A comparison of the Li-clinopyroxenes (M1 = Cr, Al, Sc, Ga, Mg + Fe) previously investigated and our sample shows that their elastic behavior and structural mechanisms of compression are analogous.     

Physical properties of calcium aluminates from vibrational spectroscopy

Heat capacity (C_V) and entropy (S) as a function of temperature were calculated for phases in the CaO-Al₂O₃ system from vibrational spectra using a quasi-harmonic model. Calculated values of C_V at 298 K for the calcium aluminates may have uncertainties as small as ±1%, based on comparison with published calorimetric data for CaO, Al₂O₃, and CaAl₂O₄. For hibonite (CaAl₁₂O₁₉), we predict C_P as 519.3 J/mol-K and S as 391.7 J/mol-K, at 298 K. For grossite (CaAl₄O₇), calculated values of C_P as 195.9 J/mol-K and of S as 172.0 J/mol-K are slightly smaller than the available calorimetric data at 298 K, consistent with calorimetric data having been obtained from samples containing ∼10 wt% hibonite impurities. Thermal conductivity at 298 K (k₀) is predicted from peak widths of the vibrational modes using the damped harmonic oscillator model of a phonon gas. Calculations of k₀ for CaO and Al₂O₃ differ from the measurements by 17% and 5%, respectively. The discrepancy for lime is larger due to uncertainties in its peak widths. This comparison suggests that our results for the calcium aluminates should be more accurate than conventional measurements of k₀ which are commonly uncertain by ∼25%.     

Placing constraints on phase equilibria and thermophysical properties in the system MgO-SiO2 by a thermodynamically consistent vibrational method

We use a lattice vibrational technique to derive thermophysical and thermochemical properties and phase equilibria in the system MgO-SiO₂ at pressures and temperatures relevant to Earth’s mantle. The technique is based on an extension of Kieffer’s model to incorporate details of the phonon spectrum, and it includes treatment of intrinsic anharmonicity. We use a least squares inversion technique applied to available experimental data, and show that it results in an accurate representation of thermodynamic properties and sound wave velocities of high-pressure phases in the system MgSiO₃. The vibrational method has been validated against laboratory experimental data in the temperature range between 0 and 2500 K and at pressures between 1 bar and 30 GPa. The technique results in a phase diagram consistent with the majority of thermophysical and thermochemical data. It is shown that intrinsic anharmonicity affects significantly slopes and positions of the phase boundaries. Our analysis indicates inconsistencies in a number of data sets of thermophysical properties for stishovite, majorite and ortho-enstatite necessitating new measurements. For akimotoite elasticity data at high-pressure and high-temperature conditions and 1 bar heat capacity measurements are needed. For stishovite elasticity measurements are necessary to reconcile elasticity data with V-P-T measurements. Additionally V-P-T measurements at pressures higher than 10 GPa are needed, which should be reconciled with V-P-T data at lower pressures. Raman and infrared spectroscopic data are necessary for both clino-enstatite and majorite. Additionally structural data are needed to resolve the discrepancy between values for the degree of disorder in majorite. Volume expansion data for ortho-enstatite are needed and effects causing differences in measured volume expansion should be elucidated.     

Determination of the solubility products of sodium carbonate minerals and an application to trona deposition in Lake Magadi (Kenya)

The ion-interaction model of PITZER (1973), is very effective in deriving stability relationships at high concentrations for the system Na-Cl-HCO₃-CO₃-OH-H₂O. The solubility products of the main sodium carbonates have been calculated from solubility data between 5 and 50°C. The stability diagram in log p_co2 — temperature coordinates and the invariant points deduced from the newly determined data are in good agreement with the most recent measurements.These results are used to calculate the activities of the major dissolved species in Lake Magadi brines (Kenya). The thermodynamic treatment confirms the main conclusions reached earlier by Eugster (1970, 1980) mainly from field observations. Trona precipitation occurs at equilibrium while natron is likely to form when the temperature decreases below 25°C. After the salt deposition the CO₂ supply from the atmosphere is too slow to allow equilibrium between the atmosphere and the brines. In the next stages of evaporative concentration thermonatrite and halite precipitate. The deposition of the latter salts along with the observed HCO^?₃ depletion suggest that fractional crystallization is likely to control trona deposition.     

Methane: An equation of state with application to the ternary system H2O-CO2-CH4

The following hardsphere modified Redlich-Kwong (HSMRK) equation of state was obtained by least squares fitting to available P-V-T data for methane (P in bars; T in Kelvins; v in cm³ mol^?1; b = 60.00 cm³ mol^?1; R = 83.14 cm³barmol^?1K^?1): $\text{P}\text{RT(1 + y + y}{}^2\text{?y}{}^3\text{v(1?y)}{}^3\text{)}\text{-}\text{c(T) + d(T)}\text{v}\text{+}\text{e(T)}\text{v}{}^2\text{/v(v + b)T}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}$$\text{y =}\text{b}\text{4v}$c(T) = 13.403 × 10⁶ + (9.28 × 10⁴)T + 2.7 T²d(T) = 5.216 × 10⁹ ? (6.8 × 10⁶)T + (3.28 × 10³)T²e(T) = (?2.3322 × 10¹¹) + (6.738 × 10⁸)T + (3.179 × 10⁵)T² For the P-T range of experimental data used in the fit (50 to 8600 bars and from 320 to 670 K), calculated volumes and fugacity coefficients for CH₄ relative to experimentally determined volumes and fugacity coefficients have average percent deviations of 0.279 and 1.373, respectively. The HSMRK equation, which predicts linear isochores over a wide P-T range, should yield reasonable estimates of fugacity coefficients for CH₄ to pressures and temperatures well outside the P-T range of available P-V-T data. Calculations for the system H₂O-CO₂-CH₄, using the HSMRK equations for H₂O and CO₂ of Kerrick and Jacobs (1981) and the HSMRK equation for CH₄ of this study, indicate that compared to the binary H₂O-CO₂ system, small amounts of CH₄ in the ternary system H₂O-CO₂-CH₄ slightly increases the activity of H₂O, and significantly decreases the activity of CO₂.     

Systematics of metal-silicate partitioning for many siderophile elements applied to Earth’s core formation

Superliquidus metal-silicate partitioning was investigated for a number of moderately siderophile (Mo, As, Ge, W, P, Ni, Co), slightly siderophile (Zn, Ga, Mn, V, Cr) and refractory lithophile (Nb, Ta) elements. To provide independent constrains on the effects of temperature, oxygen fugacity and silicate melt composition, isobaric (3 GPa) experiments were conducted in piston cylinder apparatus at temperature between 1600 and 2600 °C, relative oxygen fugacities of IW−1.5 to IW−3.5, and for silicate melt compositions ranging from basalt to peridotite. The effect of pressure was investigated through a combination of piston cylinder and multi-anvil isothermal experiments between 0.5 and 18 GPa at 1900 °C. Oxidation states of siderophile elements in the silicate melt as well as effect of carbon saturation on partitioning are also derived from these results. For some elements (e.g. Ga, Ge, W, V, Zn) the observed temperature dependence does not define trends parallel to those modeled using metal-metal oxide free energy data. We correct partitioning data for solute interactions in the metallic liquid and provide a parameterization utilized in extrapolating these results to the P-T-X conditions proposed by various core formation models. A single-stage core formation model reproduces the mantle abundances of several siderophile elements (Ni, Co, Cr, Mn, Mo, W, Zn) for core-mantle equilibration at pressures from 32 to 42 GPa along the solidus of a deep peridotitic magma ocean (∼3000 K for this pressure range) and oxygen fugacities relevant to the FeO content of the present-day mantle. However, these P-T-fO₂ conditions cannot produce the observed concentrations of Ga, Ge, V, Nb, As and P. For more reducing conditions, the P-T solution domain for single stage core formation occurs at subsolidus conditions and still cannot account for the abundances of Ge, Nb and P. Continuous core formation at the base of a magma ocean at P-T conditions constrained by the peridotite liquidus and fixed fO₂ yields concentrations matching observed values for Ni, Co, Cr, Zn, Mn and W but underestimates the core/mantle partitioning observed for other elements, notably V, which can be reconciled if accretion began under reducing conditions with progressive oxidation to fO₂ conditions consistent with the current concentration of FeO in the mantle as proposed by Wade and Wood (2005). However, neither oxygen fugacity path is capable of accounting for the depletions of Ga and Ge in the Earth’s mantle. To better understand core formation, we need further tests integrating the currently poorly-known effects of light elements and more complex conditions of accretion and differentiation such as giant impacts and incomplete equilibration.     

Proposed method for groundwater vulnerability mapping in carbonate (karstic) aquifers: the COP method

The ‘COP method’ has been developed for the assessment of intrinsic vulnerability of carbonate aquifers in the frame of the European COST Action 620. This method uses the properties of overlying layers above the water table (O factor), the concentration of flow (C factor) and precipitation (P factor) over the aquifer, as the parameters to assess the intrinsic vulnerability of groundwater. This method considers karst characteristics, such as the presence of swallow holes (C factor) and their catchment areas as well as karstic landforms, as factors which decrease the natural protection provided by overlying layers (O factor). The P factor allows for consideration of the spatial and temporal variability of precipitation, which is considered the transport agent of contamination. Two carbonate aquifers in the South of Spain, Sierra de Líbar (a conduit flow system) and Torremolinos (a diffuse flow system), have been selected for the application and validation of the method and the results have been compared with three methods widely applied in different aquifers around the world (AVI, GOD and DRASTIC). Comparisons with these methods and validation tools (hydrogeological data and tracer test) show the advantages of the COP method in the assessment of vulnerability of karstic groundwaters.     

Vibrational spectroscopy of beryllium aluminosilicates: Heat capacity calculations from band assignments

Far-infrared, mid-IR, and Raman powder spectra were measured on six phases (bromellite, chrysoberyl, phenakite, bertrandite, beryl, and euclase) in the system BeO-Al₂O₃-SiO₂-H₂O. A single-crystal absorption spectrum of IR fundamentals in beryl is also presented, which more closely resembles the powder absorption spectrum than it does absorption spectra calculated from single-crystal reflection data. Assignments of the SiO₄ and BeO₄ internal vibrations are made in accordance with each mineral's symmetry and composition and by comparison to structural analogs. Heat capacities C _v calculated for these partial band assignments agree with C _v derived from experimental C _p for all six phases, provided that Kieffer's (1979c) model is slightly modified to correctly enumerate both Si-O and Be-O stretching modes in the high frequency region (>750 cm^?1). Si-O stretching bands were found to out-number Be-O stretching modes in the high-energy region of the vibrational spectra with two exceptions: (1) For those phases containing oxygen ions not coordinated to silicon, vibrations occurring at v>1,080 cm^?1 that are attributable to Be-O (H) stretching must be treated separately in the model in order to calculate C _v accurately. (2) Minerals consisting entirely of interlocking Si and Be tetrahedra (i.e., phases without Al or OH) can be modeled by one optic continuum representing all optical modes. These results, along with the occurrence of very low energy lattice vibrations for Be-silicates within Al, suggests that although Be-O bonds are generally weaker than neighboring Si-O bonds, Be mimics the network-forming characteristic of Si to a limited extent.     

The calculated geometry of silicate liquid immiscibility

Calculations on a new empirical solution model, constructed as a Kohler oxide solution with special oxide coefficients and nearly symmetrical binary interaction parameters, reproduce the simple geometry found for silicate liquid immiscibility in synthetic systems. Both ordinary regular and subregular solutions are unsuitable for these calculations because, as components are added, these solutions result in complicated liquid immiscibility with many extraneous solvii. The geometry of immiscibility is determined by spinodal calculation followed by a graphical construction. Spinodal calculation produces a vector R which must be roughly parallel to tie lines and is used to compare the real with modelled partitioning. The results agree very well with the available experimental work for SiO₂-M(1)O-M(2)O systems M = Fe, Mg, Ca, Mn, Zn, Ba, Pb) and reasonably well for the quartz-fayalite-leucite section with minor Fe₂O₃, TiO₂ and P₂O₅. The new model predicts that the flattened liquidus of diopside, in the diopside-leucite-quartz and diopside-nepheline-quartz systems, is due to the metastable extension of the diopside-quartz melt solvus, only 100–150°C below the liquidus of diopside. Preliminary attempts to extend the coefficient model to natural examples of magma immiscibility are not very successful.     

Ab initio calculations of 17O and nT NMR parameters (nT=31P, 29Si) in H3 TOTH3 dimers and T 3O9 trimeric rings

NMR shieldings (σ) and electric field gradients (eq) are calculated using ab initio methods at the O and T nuclei (where T=P, Si) in two different types of molecules-TH₃ dimers, i.e. H₃SiOSiH₃ and H₃POPH ₃ ²⁺ , and TO₄ trimeric rings, i.e., Si₃O ₉ ^6- and P₃O ₉ ^3- , which serve as models for assessing the effects of polymerization, bond length and bond angle variation on the NMR properties of polymerized silicates and phosphates. In agreement with earlier ab initio studies on H₃SiOSiH₃ we confirm that σ(²⁹Si), σ(³¹P), σ(¹⁷O) and eq(¹⁷O) all decrease as θ(SiOSi) decreases in the range from 180° to 100°. However, correction for artifacts due to distant core electrons leads to a considerably reduced value for the anisotropy in σ ^O, bringing it into better agreement with estimated experimental values. The qualitative change in σ(²⁹Si) with θ(SiOSi) can be understood on the basis of changes in the energies of the highest energy occupied MO's and consequent variations in their contributions to the paramagnetic part of the shielding. For H₃POPH ₃ ²⁺ we calculate a larger value of eq^O than for the analog Si compound but the same type of variation of σ(¹⁷O) with θ(TOT). The change in σ(³¹P) with θ(POP) is, however, calculated to be much smaller than in the Si case and a maximum is predicted for intermediate angles. For the trimeric rings we obtain energy optimized geometries in good agreement with x-ray structural data, with T-O terminal distances systematically shorter than the T-O bridging distances. Calculated σ(T) anisotropies are also in good agreement with experiment and can be simply related to the calculated structure. After correction for distant core effects we obtain a change in σ(³¹P) between PO ₄ ^3- and P₃O ₉ ^3- in reasonable agreement with experiment.     

Pyrope–knorringite garnets: overview of experimental data and natural parageneses

This paper gives an analytical overview of the experimental data obtained by different authors at high P and T in the model system MgO–Al₂O₃–SiO₂–Cr₂O₃ (MASCr). A set of four simple polynomial equations is proposed for the temperature and pressure dependence of chromium content in garnet and spinel in the assemblage Gar + Opx + Es and Gar + Fo + Opx + Sp.From the first equation, one can estimate the minimum pressure at a given temperature which is required for the formation of peridotite garnets of uncertain paragenesis with a known knorringite content. A combination of the second and third equations helps estimate P and T from the chromium content of garnet and spinel from assemblages containing both minerals. If the spinel composition is unknown, but there is reason to assign garnet to a spinel-bearing paragenesis, the fourth equation is applicable for estimating pressure at given temperature.Originally, the proposed garnet–spinel geothermobarometer was developed for a harzburgite paragenesis. However, it is applicable to garnets with CaO/Cr₂O₃ < 0.903 (including lherzolitic ones), that is, those within the Pyr–Kn–Uv triangle of the reciprocal quaternary diagram Pyr–Cros–Uv–Kn.Using the above equations and an empirical P_CG geobarometer (Grütter et al., 2006), comparative geothermobarometric estimates were obtained for a set of garnet and garnet–spinel inclusions in diamonds and intergrowths with diamond, as well as garnet inclusions in spinel. If garnet has CaO/Cr₂O₃ = 0.35–0.40, the results are in good accord. For Cr-richest and Ca-poorest garnets, the P_CG barometer shows pressures 10–15% higher compared with our estimates.     

The crystal and magnetic structure of Li-aegirine LiFe3+Si2O6: a temperature-dependent study

Single crystals of Li-aegirine LiFe³⁺Si₂O₆ were synthesized at 1573?K and 3?GPa, and a polycrystalline sample suitable for neutron diffraction was produced by ceramic sintering at 1223?K. LiFe³⁺Si₂O₆ is monoclinic, space group C2/c, a=9.6641(2)?Å, b= 8.6612(3)?Å, c=5.2924(2)?Å, β=110.12(1)^° at 300?K as refined from powder neutron data. At 229?K Li-aegirine undergoes a phase transition from C2/c to P2₁ /c. This is indicated by strong discontinuities in the temperature variation of the lattice parameters, especially for the monoclinic angle β and by the appearance of Bragg reflections (hkl) with h+k≠2n. In the low-temperature form two non-equivalent Si-sites with 〈SiA–O〉=1.622?Å and 〈SiB–O〉=1.624?Å at 100?K are present. The bridging angles of the SiO₄ tetrahedra O3–O3–O3 are 192.55(8)° and 160.02(9)° at 100?K in the two independent tetrahedral chains in space group P2₁ /c, whereas it is 180.83(9)° at 300?K in the high-temperature C2/c phase, i.e. the chains are nearly fully expanded. Upon the phase transition the Li-coordination changes from six to five. At 100?K four Li–O bond lengths lie within 2.072(4)–2.172(3)?Å, the fifth Li–O bond length is 2.356(4)?Å, whereas the Li–O3?A bond lengths amount to 2.796(4)?Å. From ⁵⁷Fe Mössbauer spectroscopic measurements between 80 and 500?K the structural phase transition is characterized by a small discontinuity of the quadrupole splitting. Temperature-dependent neutron powder diffraction experiments show first occurrence of magnetic reflections at 16.5?K in good agreement with the point of inflection in the temperature-dependent magnetization of LiFe³⁺Si₂O₆. Distinct preordering phenomena can be observed up to 35?K. At the magnetic phase transition the unit cell parameters exhibit a pronounced magneto-striction of the lattice. Below T _N Li-aegirine shows a collinear antiferromagnetic structure. From our neutron powder diffraction experiments we extract a collinear antiferromagnetic spin arrangement within the a–c plane.     

The reaction forsterite+cordierite=aluminous orthopyroxene+spinel in the system MgO-Al2O3-SiO2

Reversal experiments at 1,150–1,300°C on the reaction forsterite+cordierite=aluminous orthopyroxene+spinel in the system MgO-Al₂O₃-SiO₂ show the equilibrium to have a negativedT/dP. The slope andT-P location of this equilibrium have been modelled using available heat capacity data and various structural models which explore the configurational entropy contributions to the totalΔS. The experimental data are consistent with the aluminous orthopyroxene model of Ganguly and Ghose (1979) where limited Al disorder occurs between theM1 andM2 sites, Al-Si mixing occurs on the tetrahedralB site with the ‘aluminum avoidance’ principle maintained, and Mg-Al disorder occurs in spinel with an interchange enthalpy of 9–12 kcal mol^?1. Additionally, Al-Si disordering which occurs in the indialite structure of cordierite is inconsistent with the experimental data and all pyroxene and spinel mixing models; consequently, Si and Al in anhydrous cordierites to 1,300°C in the system MgO-Al₂O₃-SiO₂ must be largely ordered.     

P-T-X data on P21/c-clinopyroxenes and their displacive phase transitions

The P2₁/c clinopyroxene kanoite (ideally MnMgSi₂O₆) was studied as a function of pressure and temperature using powder X-ray diffraction, differential scanning calorimetry (DSC) and optical methods. The temperature of the P2₁/c to high-temperature (HT) C2/c transition ranges from 425?°C in endmember MnMgSi₂O₆ to 125?°C in natural samples with an aegirine component. Compiling pigeonite and clinoenstatite–clinoferrosilite literature data, the temperature of the transformation was found to decrease linearly with M2 cation size. A synchrotron powder diffraction study in a heated diamond-anvil cell (DAC) yielded compression and thermal expansion data for low kanoite of composition Mn_1.2Mg_0.4Fe_0.4Si₂O₆. The high-pressure (HP) phase transition from P2₁/c to HPC2/c was reversed at 5.8 GPa at 417?°C. The high-temperature phase transition from P2₁/c to HTC2/c was rather indistinct and occurred at approximately 530?°C and 0.76 GPa. In a separate experiment, the HT transition was observed optically in a hydrothermal DAC between 0.0 and 0.4 GPa. The in-situP-T data of both experiments yielded an increase in transition temperature with increasing pressure (approx. 149?°C/GPa) and suggest a change in character of the transition from first order to continuous with increasing pressure. The data indicate that the HTC2/c and HPC2/c polymorphs are distinct phases with different stability fields. Since the high-temperature and the high-pressure polymorphs of kanoite were shown to be isotypic with other low-Ca clinopyroxenes such as the (Mg,Fe)SiO₃ series, the conclusions we draw from this study are valid for all clinopyroxenes with small (<0.88 Å) M1 and M2 cation sizes. The petrologic implications of these conclusions for the occurrence of “clinoenstatite” in the Alpe Arami peridotite are discussed.     

Metamorphic equilibria in the siliceous dolomite system: 6 kbar experimental data and geologic implications

Hydrothermal experiments with H₂O-CO₂ fluids at P_fluid = 6 kbar yielded the following quilibrium conditions for reactions important in metamorphosed siliceous dolomites (T = °C; X = X_co2): (3) dolomite + 2 quartz = diopside + 2 CO₂T = 620 ± 8X = 0.73 ± 0.03 (5) 5 dolomite + 8 quartz + H₂O = tremolite + 3 calcite + 7 CO₂T = 600 ± 5 550 ±5 540±5 500±5X = 0.66 ± 0.03 0.21 ± 0.03 0.21 ± 0.04 0.06 ± 0.02 (7) 3 dolomite + 4 quartz + H₂O = talc + 3 calcite + 3 CO₂T = 550±5 500±5 450 ±5X = 0.25 ± 0.05 0.07 ± 0.02 0.03 ± 0.02 (8) 2 dolomite + talc + 4 quartz = tremolite + 4 CO₂T = 550 ± 5 540 ±5 500 ± 5X = 0.22 ± 0.03 0.21 ± 0.02 0.06 ± 0.02 A thermodynamically self-consistent 6 kbar T-X_CO2, topology results by extrapolating equilibria from experimental brackets using a modified Redlich-Kwong equation for activities in H₂O-CO₂ mixtures. This topology restricts the assemblage talc + calcite to a narrow stability band in T-X_CO2 space at X_CO2 < 0.55 and T < 590°C. Accordingly, the occurrence of talc + calcite in pure siliceous dolomites metamorphosed at P_fluid = 6 kbar implies correspondingly water-rich fluids.     

Zn and Pb solubility and speciation in aqueous chloride fluids at T-P parameters corresponding to granitoid magma degassing and crystallization

The solubility of all possible Zn and Pb species in aqueous chloride fluids was evaluated by means of thermodynamic simulations in systems ZnO(PbO)-aqueous solution of NaCl (KCl, NaCl + HCl) within broad ranges of temperature (600–900°C), pressure (0.7–5 kbar), and chloride concentrations, under parameters corresponding to the crystallization and degassing of granitoid magmas in the Earth’s crust. Our simulation results demonstrate that the addition of Cl to the fluid phase in the form of Na(K)Cl and HCl significantly increases the concentrations of Cl-bearing Zn and Pb complexes and the total concentration of the metals in the solutions in equilibrium with the solid oxides. In Zn-bearing fluids, the Zn(OH) ₂ ⁰ , ZnOH⁺, and Zn(OH) ₃ ^? —hydroxyl complexes and the ZnCl ₂ ⁰ , and ZnCl⁺ chlorocomplexes, which are predominant at low Cl concentrations (C_Cl < 0.05–0.1 m) give way to ZnCl ₄ ^2? with increasing C_Cl, which becomes the predominant Zn species of the fluid at C_Cl > 0.1–0.5 m throughout the whole temperature range in question and pressures higher than 1 kbar. For Pb-bearing fluids, the T-P-X region dominated by the Pb(OH) ₂ ⁰ , and Pb(OH) ₃ ^? hydroxyl complexes is remarkably wider than the analogous region for Zn, particularly at elevated temperatures (≥700°C) in alkaline solutions. An increase in C_Cl is associated with an increase in the concentration and changes in the speciation of Pb chlorocomplexes: PbCl ₂ ⁰ → PbCl ₃ ^? → PbCl ₄ ^2? . The concentrations of Zn and Pb chlorocomplexes increase with increasing pressure, decreasing temperature, and decrease pH with the addition of HCl to the system. It is demonstrated that the solubility of ZnO at any given T-P-X in alkaline solutions with low chloride concentrations are lower than the solubility of PbO. The Zn concentration increases more significantly than with the Pb concentration with increasing C_Cl and decreasing pH, so that the Zn concentration in acidic solutions is higher than the Pb concentration over broad ranges of temperature, pressure, and Cl concentration. Chloride complexes of Zn (ZnCl ₂ ⁰ , and ZnCl ₄ ^2? ) and Pb (PbCl ₂ ⁰ , and PbCl ₃ ^? are proved to be predominant within broad T-P-X-pH ranges corresponding to the parameters under which magmatic fluid are generated. Our simulation results confirm the hypothesis that chlorocomplexes play a leading role in Zn and Pb distribution between aqueous chloride fluids and granitic melts. These simulation results are consistent with experimental data on the Zn and Pb distribution coefficients (D(Zn)^f/m and D(Pb)^f/m, respectively) between aqueous chloride fluids and granitic melts that demonstrated that (1) D(Zn)^f/m and D(Pb)^f/m increase with increasing Na and K chloride concentrations in the aqueous fluid, (2) both D(Zn)^f/m and D(Pb)^f/m drastically increase when HCl is added to the fluid, and (3) (D(Zn)^f/m is higher than D(Pb)^f/m at any given T-P-X parameters. The experimentally established decrease in D(Zn)^f/m and D(Pb)^f/m with increasing pressure (at unchanging temperature and Cl concentration) is likely explained by an increase in the alkalinity of the aqueous chloride fluid in equilibrium with granite melt and, correspondingly, a decrease in the Zn and Pb solubility in this fluid.