SIMS determination of trace element partition coefficients between garnet, clinopyroxene and hydrous basaltic liquids at 2–7.5 GPa and 1080–1200°C

Trace element partition coefficients (D's) for up to 13 REE, Nb, Ta, Zr, Hf, Sr and Y have been determined by SIMS analysis of seven garnets, four clinopyroxenes, one orthopyroxene and one phlogopite crystallized from an undoped basanite and a lightly doped (200 ppm Nb, Ta and Hf) quartz tholeiite. Experiments were conducted at 2–7.5 GPa, achieving near-liquidus crystallization at relatively low temperatures of 1080–1200°C under strongly hydrous conditions (5–27 wt.% added water). Garnet and pyroxene D_REE show a parabolic pattern when plotted against ionic radius, and conform closely to the lattice strain model of Blundy and Wood (Blundy, J.D., Wood, B.J., 1994. Prediction of crystal–melt partition coefficients from elastic moduli. Nature 372, 452–454). Comparison, at constant pressure, between hydrous and anhydrous values of the strain-free partition coefficient (D₀) for the large cation sites of garnet and clinopyroxene reveals the relative importance of temperature and melt water content on partitioning. In the case of garnet, the effect of lower temperature, which serves to increase D₀, and higher water content, which serves to decrease D₀, counteract each other to the extent that water has little effect on garnet–melt D₀ values. In contrast, the effect of water on clinopyroxene–melt D₀ overwhelms the effect of temperature, such that D₀ is significantly lower under hydrous conditions. For both minerals, however, the lower temperature of the hydrous experiments tends to tighten the partitioning parabolas, increasing fractionation of light from heavy REE compared to anhydrous experiments.

Three sets of near-liquidus clinopyroxene–garnet two-mineral D values increase the range of published experimental determinations, but show significant differences from natural two-mineral D's determined for subsolidus mineral pairs. Similar behaviour is observed for the first experimental data for orthopyroxene–clinopyroxene two-mineral D's when compared with natural data. These differences are in large part of a consequence of the subsolidus equilibration temperatures and compositions of natural mineral pairs. Great care should therefore be taken when using natural mineral–mineral partition coefficients to interpret magmatic processes.

The new data for strongly hydrous compositions suggest that fractionation of Zr–Hf–Sm by garnet decreases with increasing depth. Thus, melts leaving a garnet-dominated residuum at depths of about 200 km or greater may preserve source Zr/Hf and Hf/Sm. This contrasts with melting at shallower depths where both garnet and clinopyroxene will cause Zr–Hf–Sm fractionation. Also, at shallower depths, clinopyroxene-dominated fractionation may produce a positive Sr spike in melts from spinel lherzolite, but for garnet lherzolite melting, no Sr spike will result. Conversely, clinopyroxene megacrysts with negative Sr spikes may crystallize from magmas without anomalous Sr contents when plotted on mantle compatibility diagrams. Because the characteristics of strongly hydrous silicate melt and solute-rich aqueous fluid converge at high pressure, the hydrous data presented here are particularly pertinent to modelling processes in subduction zones, where aqueous fluids may have an important metasomatic role.     

Effect of Fe on garnet–melt trace element partitioning: experiments in FCMAS and quantification of crystal-chemical controls in natural systems

Garnet–melt trace element partitioning experiments were performed in the system FeO–CaO–MgO–Al₂O₃–SiO₂ (FCMAS) at 3 GPa and 1540°C, aimed specifically at studying the effect of garnet Fe²⁺ content on partition coefficients (D^Grt/Melt). D^Grt/Melt, measured by SIMS, for trivalent elements entering the garnet X-site show a small but significant dependence on garnet almandine content. This dependence is rationalised using the lattice strain model of Blundy and Wood [Blundy, J.D., Wood, B.J., 1994. Prediction of crystal–melt partition coefficients from elastic moduli. Nature 372, 452–454], which describes partitioning of an element i with radius r_i and valency Z in terms of three parameters: the effective radius of the site r₀(Z), the strain-free partition coefficient D₀(Z) for a cation with radius r₀(Z), and the apparent compressibility of the garnet X-site given by its Young's modulus E_X(Z). Combination of these results with data in Fe-free systems [Van Westrenen, W., Blundy, J.D., Wood, B.J., 1999. Crystal-chemical controls on trace element partitioning between garnet and anhydrous silicate melt. Am. Mineral. 84, 838–847] and crystal structure data for spessartine, andradite, and uvarovite, leads to the following equations for r₀(3+) and E_X(3+) as a function of garnet composition (X) and pressure (P):

r₀(3+) [Å]=0.930X_Py+0.993X_Gr+0.916X_Alm+0.946X_Spes+1.05(X_And+X_Uv)−0.005(P [GPa]−3.0)(±0.005 Å)

E_X(3+) [GPa]=3.5×10¹²(1.38+r₀(3+) [Å])^−26.7(±30 GPa)

Accuracy of these equations is shown by application to the existing garnet–melt partitioning database, covering a wide range of P and T conditions (1.8 GPa<P<5.0 GPa; 975°C<T<1640°C). D^Grt/Melt for all 3+ elements entering the X-site (REE, Sc and Y) are predicted to within 10–40% at given P, T, and X, when D^Grt/Melt for just one of these elements is known. In the absence of such knowledge, relative element fractionation (e.g. D_Sm^Grt/Melt/D_Nd^Grt/Melt) can be predicted. As an example, we predict that during partial melting of garnet peridotite, group A eclogite, and garnet pyroxenite, r₀(3+) for garnets ranges from 0.939±0.005 to 0.953±0.009 Å. These values are consistently smaller than the ionic radius of the heaviest REE, Lu. The above equations quantify the crystal-chemical controls on garnet–melt partitioning for the REE, Y and Sc. As such, they represent a major advance en route to predicting D^Grt/Melt for these elements as a function of P, T and X.     

Partitioning coefficients between olivine and silicate melts

Variation of Nernst partition coefficients (D) between olivine and silicate melts cannot be neglected when modeling partial melting and fractional crystallization. Published natural and experimental ^{olivine/liquid}D data were examined for covariation with pressure, temperature, olivine forsterite content, and melt SiO₂, H₂O, MgO and MgO/MgO + FeO^total. Values of ^{olivine/liquid}D generally increase with decreasing temperature and melt MgO content, and with increasing melt SiO₂ content, but generally show poor correlations with other variables. Multi-element ^{olivine/liquid}D profiles calculated from regressions of D REE–Sc–Y vs. melt MgO content are compared to results of the Lattice Strain Model to link melt MgO and: D₀ (the strain compensated partition coefficient), E_M³⁺ (Young's Modulus), and r₀ (the size of the M site). Ln D₀ varies linearly with Ln MgO in the melt; E_M³⁺ varies linearly with melt MgO, with a dog-leg at ca. 1.5% MgO; and r₀ remains constant at 0.807 Å. These equations are then used to calculate ^{olivine/liquid}D for these elements using the Lattice Strain Model. These empirical parameterizations of ^{olivine/liquid}D variations yield results comparable to experimental or natural partitioning data, and can easily be integrated into existing trace element modeling algorithms. The ^{olivine/liquid}D data suggest that basaltic melts in equilibrium with pure olivine may acquire small negative Ta–Hf–Zr–Ti anomalies, but that negative Nb anomalies are unlikely to develop. Misfits between results of the Lattice Strain Model and most light rare earth and large ion lithophile partitioning data suggest that kinetic effects may limit the lower value of D for extremely incompatible elements in natural situations characterized by high cooling/crystallization rates.     

Marcasite precipitation from hydrothermal solutions

Pyrite and marcasite were precipitated by both slow addition of aqueous Fe²⁺ and SiO₃²⁻ to an H₂S solution and by mixing aqueous Fe²⁺ and Na₂S₄ solutions at 75°C. H₂S₂ or HS₂⁻ and H₂S₄ or HS₄⁻ were formed in the S₂O₃²⁻ and Na₂S₄ experiments, respectively. Marcasite formed at pH < pK₁ of the polysulfide species present (for H₂S₂, pK₁ = 5.0; for H₂S₄, pK₁ = 3.8 at 25°C). Marcasite forms when the neutral sulfane is the dominant polysulfide, whereas pyrite forms when mono-or divalent polysulfides are dominant. In natural solutions where H₂S₂ and HS₂ are likely to be the dominant polysulfides, marcasite will form only below pH 5 at all temperatures.

The pH-dependent precipitation of pyrite and marcasite may be caused by electrostatic interactions between polysulfide species and pyrite or marcasite growth surfaces: the protonated ends of H₂S₂ and HS₂ are repelled from pyrite growth sites but not from marcasite growth sites. The negative ions HS₂ and S₂²⁻ are strongly attracted to the positive pyrite growth sites. Masking of 1π_g^* electrons in the S₂ group by the protons makes HS₂ and H₂S₂ isoelectronic with AsS²⁻ and As₂²⁻, respectively ( et al., 1981). Thus, the loellingitederivative structure (marcasite) results when both ends of the polysulfide are protonated.

Marcasite occurs abundantly only for conditions below pH 5 and where H₂S₂ was formed near the site of deposition by either partial oxidation of aqueous H₂S by O₂ or by the reaction of higher oxidation state sulfur species that are reactive with H₂S at the conditions of formation e.g., S₂O₃²⁻ but not SO₄²⁻. The temperature of formation of natural marcasite may be as high as 240°C ( and , 1985), but preservation on a multimillion-year scale seems to require post-depositional temperatures of below about 160°C ( , 1973; and , 1985).     

川东北普光地区深层富钾卤水钾资源提取研究

四川盆地地下卤水资源丰富,尤其川东北地区地下卤水富含高品质钾资源。本项目以四川普光地区富钾卤水为研究对象,根据卤水组成,采用模拟计算并结合实验验证的方法,研究了高温蒸发时,氯化钠、氯化钾、光卤石、硼酸等矿物的析出阶段及特点。研究结果表明,当蒸失水率约80%时,体系中约85%的NaCl析出,同时KCl达到饱和,继续蒸发可获取钾石盐。控制总蒸失水率94%~95%时分离,体系中大于80%的钾可在这一阶段析出,且湿基中KCl品位可高达约45%;析出钾石盐后的卤水降至室温可获得NaCl、KCl、KCl·MgCl₂·6H₂O及H₃BO₃的混合物,硼的析出率可达到约80%。同时,研究表明,将蒸失水率约80%时分离石盐后的卤水直接降温可以获得湿基品位约高达70%的钾石盐矿,钾析出率约50%。综合对比分析,提出两条以钾资源开发为主的工艺路线:其一为“高温蒸发析氯化钠-高温蒸发析氯化钾-冷却析钾硼混盐”,通过加工获得KCl及H₃BO₃产品;其二为“高温蒸发析氯化钠-冷却析钾”,通过加工获得高端KCl产品。     

An experimental study of calcite dissolution rates at acidic conditions and 25 °C in the presence of NaPO3 and MgCl2

Dissolution rates of single calcite crystals were determined from sample weight loss using free-drift rotating disk techniques. Experiments were performed at 25 °C in aqueous HCl solutions over the bulk solution pH range −1 to 3 and in the presence of trace concentrations of aqueous NaPO₃ and MgCl₂. These salts were chosen for this study because aqueous magnesium and phosphate are known to strongly inhibit calcite dissolution at neutral to basic pH. Reactive solutions were undersaturated with respect to possible secondary phases. Neither an inhibition or enhancement of calcite dissolution rates was observed in the presence of aqueous MgCl₂ at pH 1 and 3. The presence of trace quantities of NaPO₃, which dissociates in solution to Na⁺ and H₂PO₄⁻, decreased the overall calcite dissolution rate at pH≤2. This contrasting behavior could be attributed to the different adsorption behavior of these dissolved species. As calcite surfaces are positively charged in acidic solutions, aqueous Mg²⁺ may not adsorb, whereas aqueous phosphate, present as either the anion H₂PO₄⁻ or the neutral species H₃PO₄⁰, readily adsorbs on calcite surfaces leading to significant dissolution inhibition.     

Liquid immiscibility and its role at the magmatic–hydrothermal transition: a summary of experimental studies

Several types of fluid immiscibility may affect the evolution of volatile-rich magmatic systems at the magmatic–hydrothermal transition. The topology of silicate–salt–H₂O systems implies that three-fluid immiscibility (silicate melt+hydrosaline melt+vapour) should be stable in a broad range of compositions and P–T conditions. The most important factor controlling the immiscibility appears to be the Coulombic properties (electric charges Z and ionic radii r) of the main network-modifying cations and the capacity for immiscibility appears to decrease in the following sequence: Mg>Ca>Sr>Ba>Li>Na>K. Liquid immiscibility is enhanced in peralkaline compositions and in the presence of nonsilicate anions such as F⁻, Cl⁻, CO₃²⁻ and BO₃³⁻. In volatile-rich magmatic systems, the H₂O is likely to react with the chloride, fluoride, borate and carbonate species and the chemical effects of high-temperature hydrolysis may be greatly enhanced by phase separation in systems with multiple immiscible fluid phases. Natural granitic magmas can thus exsolve a range of chemically and physically diverse hydrosaline liquids and the role of these fluid phases is likely to be especially significant in pegmatites and Li–F rare-metal granites.     

Water content in natural eclogite and implication for water transport into the deep upper mantle

Infrared spectroscopy and ion micro-probe measurements showed that the major constituent minerals of eclogites from the Kokchetav massif, which have been subducted to 180 km depths, contain significant amounts of water up to 870 ppm H₂O (by weight) in omphacite, 130 ppm H₂O in garnet and 740 ppm H₂O in rutile. Omphacite shows three hydroxyl absorption bands at 3440–3460, 3500–3530 and 3600–3625 cm^− 1, garnet has a single band at 3580–3630 cm^− 1 and rutile has a single sharp band at 3280 cm^− 1. The hydroxyl absorbance at these wavenumbers changes with the crystal orientation in polarized infrared radiation, indicating that the water is structurally incorporated in these minerals. The water contents in omphacite and garnet increase systematically with the metamorphic pressure of the host eclogites. The partitioning coefficient of the water content between coexisting garnet and omphacite is similar in different eclogites, D^Grt/Omp0.1–0.2, but decreases slightly at high pressure. Based on the mineral proportions of the eclogites, we estimate bulk-rock water content in the eclogites ranging from 3070 to 300 ppm H₂O (by weight). Although hydrous minerals are absent in the diamond-grade eclogite (60 kbar and 1000 °C), trace amounts of water are incorporated in the nominally anhydrous minerals such as omphacite and garnet. The presence of significant water in these minerals implies that the subducting oceanic crust can transport considerable amounts of water into the deep upper mantle beyond the stability of hydrous minerals. Such water may be stored in the deep upper mantle and have an important influence on dynamics in the Earth's interior.     

Carbonatite melt–CO2 fluid inclusions in mantle xenoliths from Tenerife, Canary Islands: a story of trapping, immiscibility and fluid–rock interaction in the upper mantle

Three types of fluid inclusions have been identified in olivine porphyroclasts in the spinel harzburgite and lherzolite xenoliths from Tenerife: pure CO₂ (Type A); carbonate-rich CO₂–SO₂ mixtures (Type B); and polyphase inclusions dominated by silicate glass±fluid±sp±silicate±sulfide±carbonate (Type C). Type A inclusions commonly exhibit a “coating” (a few microns thick) consisting of an aggregate of a platy, hydrous Mg–Fe–Si phase, most likely talc, together with very small amounts of halite, dolomite and other phases. Larger crystals (e.g. (Na,K)Cl, dolomite, spinel, sulfide and phlogopite) may be found on either side of the “coating”, towards the wall of the host mineral or towards the inclusion center. These different fluids were formed through the immiscible separations and fluid–wall-rock reactions from a common, volatile-rich, siliceous, alkaline carbonatite melt infiltrating the upper mantle beneath the Tenerife. First, the original siliceous carbonatite melt is separated from a mixed CO₂–H₂O–NaCl fluid and a silicate/silicocarbonatite melt (preserved in Type A inclusions). The reaction of the carbonaceous silicate melt with the wall-rock minerals gave rise to large poikilitic orthopyroxene and clinopyroxene grains, and smaller neoblasts. During the metasomatic processes, the consumption of the silicate part of the melt produced carbonate-enriched Type B CO₂–SO₂ fluids which were trapped in exsolved orthopyroxene porphyroclasts. At the later stages, the interstitial silicate/silicocarbonatite fluids were trapped as Type C inclusions. At a temperature above 650 °C, the mixed CO₂–H₂O–NaCl fluid inside the Type A inclusions were separated into CO₂-rich fluid and H₂O–NaCl brine. At T<650 °C, the residual silicate melt reacted with the host olivine, forming a reaction rim or “coating” along the inclusion walls consisting of talc (or possibly serpentine) together with minute crystals of NaCl, KCl, carbonates and sulfides, leaving a residual CO₂ fluid. The homogenization temperatures of +2 to +25 °C obtained from the Type A CO₂ inclusions reflect the densities of the residual CO₂ after its reactions with the olivine host, and are unrelated to the initial fluid density or the external pressure at the time of trapping. The latter are restricted by the estimated crystallization temperatures of 1000–1200 °C, and the spinel lherzolite phase assemblage of the xenolith, which is 0.7–1.7 GPa.     

Fluid and melt inclusion evidence for the origin of idiomorphic quartz crystals in topaz-bearing granite from the Salmi batholith, Karelia, Russia

Idiomorphic quartz crystals in topaz-bearing granite from the Salmi batholith contain primary inclusions of silicate melt and abundant mostly secondary aqueous fluid inclusions. Microthermometric measurements on melt inclusions give estimates for the granite solidus and liquidus of 640–680°C and 770–830°C, respectively. Using published solubility models for H₂O in granitic melts and the obtained solidus/liquidus temperatures from melt inclusions, the initial water concentration of the magma is deduced to have been approximately 3 wt.% and the minimum pressure about 2 kbar. At this initial stage, volatile-undersaturation conditions of magma were assumed. These results indicate that the idiomorphic quartz crystals are magmatic in origin and thus real phenocrysts. During subsolidus cooling and fracturing of the granite, several generations of aqueous fluid inclusions were trapped into the quartz phenocrysts. The H₂O inclusions have salinities and densities of 1–41 wt.% NaCl eq. and 0.53–1.18 g/cm³, respectively.     

Fluid evolution in the nepheline syenites of the Ditrău Alkaline Massif, Transylvania, Romania

The Ditrău Alkaline Massif is an intrusion into the Bucovina nappe system that is part of the Mesozoic crystalline zone located in Transylvania, Romania, in the Eastern Carpathians. Nepheline syenites are the most abundant rocks in the central and eastern part of the Massif, and represent the last major intrusion of the complex. Fluid inclusions in nepheline, aegirine and albite were trapped at magmatic conditions on or below the H₂O-saturated nepheline syenite solidus at about 400–600 °C and 2.5–5 kbars. Early nepheline, and to a lesser extent albite, were altered by highly saline fluids to produce cancrinite, sodalite and analcime, during this process cancrinite also trapped fluid inclusions. The fluids, in most cases, can be modeled by the H₂O–NaCl system with varying salinity; however inclusions with more complex fluid composition (containing K, Ca, CO₃, etc., in addition to NaCl) are common. Raman spectroscopic analyses of daughter minerals confirm the presence of alkali-carbonate fluids in some of the earliest inclusions in nepheline, aegirine and albite.

During crystallization, the melts exsolved a high salinity, carbonate-rich magmatic fluid that evolved to lower salinity as crystallization progressed. Phases that occur early in the paragenesis contain high-salinity inclusions while late phases contain low-salinity inclusions. The salinity trend is consistent with experimental data for the partitioning of chlorine between silicic melt and exsolved aqueous fluid at about 2.0 kbars. The activity of water (a_H2O) in the melt increases during crystallization, resulting in the formation of hydrous phases during late-stage crystallization of the nepheline syenites.     

Magmatic-to-hydrothermal crystallization in the W–Sn mineralized Mole Granite (NSW, Australia): Part I: Crystallization of zircon and REE-phosphates over three million years—a geochemical and U–Pb geochronological study

Zircon, monazite and xenotime crystallized over a temperature interval of several hundred degrees at the magmatic to hydrothermal transition of the Sn and W mineralized Mole Granite. Magmatic zircon and monazite, thought to have crystallized from hydrous silicate melt, were dated by conventional U–Pb techniques at an age of 247.6 ± 0.4 and 247.7 ± 0.5 Ma, respectively. Xenotime occurring in hydrothermal quartz is found to be significantly younger at 246.2 ± 0.5 Ma and is interpreted to represent hydrothermal growth. From associated fluid inclusions it is concluded that it precipitated from a hydrothermal brine ≤ 600 °C, which is below the accepted closure temperature for U–Pb in this mineral. These data are compatible with a two-stage crystallization process: precipitation of zircon and monazite as magmatic liquidus phases in deep crustal magma followed by complete crystallization and intimately associated Sn–W mineralization after intrusion of the shallow, sill-like body of the Mole Granite. Later hydrothermal formation of monazite in a biotite–fluorite–topaz reaction rim around a mineralized vein was dated at 244.4 ± 1.4 Ma, which distinctly postdates the Mole Granite and is possibly related to a younger hidden intrusion and its hydrothermal fluid system.

Obtaining precise age data for magmatic and hydrothermal minerals of the Mole Granite is hampered by uncertainties introduced by different corrections required for multiple highly radiogenic minerals crystallising from evolved hydrous granites, including ²³⁰Th disequilibrium due to Th/U fractionation during monazite and possibly xenotime crystallization, variable Th/U ratios of the fluids from which xenotime was precipitating, elevated contents of common lead, and post-crystallization lead loss in zircon, enhanced by the fluid-saturated environment. The data imply that monazite can also survive as a liquidus phase in protracted magmatic systems over periods of 10⁶ years. The outlined model is in agreement with prominent chemical core-rim variation of the zircon.     

Platinum solubility and partitioning in a felsic melt-vapor-brine assemblage

The solubility and partitioning of Pt in a S-free vapor - brine - rhyolite melt - Pt metal assemblage has been quantified at 800 °C, f_O2=NNO and pressures of 100 and 140 MPa. Vapor and brine were sampled at run conditions by trapping these phases as glass-hosted fluid inclusions as the melt cooled through the glass transition temperature. The vapor and brine were in equilibrium with the melt at the time of trapping and, thus, represent fluids which were sampled at the termination of each experimental run. The microthermometrically determined salinities of vapor and brine are ∼2 and ∼63 wt.% NaCl eq. and ∼9 and ∼43 wt.% NaCl eq. at 100 and 140 MPa, respectively. Platinum solubilities in vapor, brine and glass (i.e., quenched melt) were quantified by using laser ablation - inductively coupled plasma - mass spectrometry (LA-ICP-MS). Equilibrium is discussed with reference to the major and trace element concentrations of glass-hosted fluid inclusions as well as the silicate melt over run times that varied from 110 to 377 h at 140 MPa and 159 to 564 h at 100 MPa. Platinum solubility values (±1σ) in H₂O-saturated felsic melt are 0.28 ± 0.13 μg/g and 0.38 ± 0.06 μg/g at 140 and 100 MPa, respectively. Platinum solubility values () at 140 and 100 MPa, respectively, in aqueous vapor are 0.91 ± 0.29 μg/g and 0.37 ± 0.17 μg/g and in are brine 16 ± 10 μg/g and 3.3 ± 1.0 μg/g. The measured solubility data were used to calculate Nernst-type partition coefficients for Pt between vapor/melt, brine/melt and vapor/brine. The partition coefficient values () for vapor/melt, brine/melt and vapor/brine at 140 MPa are 2.9 ± 1.0, 67 ± 27, and 0.13 ± 0.05 and at 100 MPa are 1.0 ± 0.2, 6.8 ± 2.4, and 0.15 ± 0.05. The partitioning data were used to model the Pt-scavenging capacity of vapor and brine during the crystallization-driven degassing (i.e., second boiling) of a felsic silicate melt over a depth range (i.e., 3-6 km) consistent with the evolution of magmatic-hydrothermal ore deposits. Model calculations suggest that aqueous vapor and brine can scavenge sufficient quantities of Pt, and by analogy other platinum group elements (PGE), to produce economically important PGE-rich magmatic-hydrothermal ore deposits in Earth’s upper continental crust.     

Solubility of silica polymorphs in electrolyte solutions, I. Activity coefficient of aqueous silica from 25° to 250°C, Pitzer''s parameterisation

Within the framework of Pitzer's specific interaction model, interaction parameters for aqueous silica in concentrated electrolyte solutions have been derived from Marshall and co-authors amorphous silica solubility measurements. The values, at 25°C, of the Pitzer interaction parameter (λ_SiO2(aq)−i) determined in this study are the following: 0.092 (i = Na⁺), 0.032 (K⁺), 0.165 (Li⁺), 0.292 (Ca²⁺, Mg²⁺), −0.139 (SO₄²⁻), and −0.009 (NO₃⁻). A set of polynomial equations has been derived which can be used to calculate λ_SiO2(aq)−i for these ions at any temperature up to 250°C. A linear relationship between the aqueous silica-ion interaction parameters (λ_SiO2(aq)−i) and the surface electrostatic field (Z_i/r_e,i) of ions was obtained. This empirical equation can be used to estimate, in first approximation, λ_SiO2(aq)−i if no measurements are available. From this parameterisation, the calculated activity coefficient of aqueous silica is 2.52 at 25°C and 1.45 at 250°C in 5 m NaCl solution. At lower concentrations, e.g. 2 m NaCl, the activity coefficient of silica is 1.45 at 25°C and 1.2 at 250°C. Hence, in practice, it is necessary to take into account the activity coefficient of aqueous silica (λ_SiO2(aq)≠1) in hydrothermal solutions and basinal brines where the ionic strength exceeds 1. A comparison of measured [Marshall, W.L., Chen, C.-T.A., 1982. Amorphous silica solubilities, V. Prediction of solubility behaviour in aqueous mixed electrolyte solutions to 300°C. Geochim. Cosmochim. Acta 46, 289–291.] and computed amorphous silica solubility, using this parameterisation, shows a good agreement. Because the effect of individual ions on silicate and silica polymorph solubilities are additive, the present study has permitted to derive Pitzer interaction parameters that allow a precise computation of γ_SiO2(aq) in the Na---K---Ca---Mg---Cl---SO₄---HCO₃---SiO₂---H₂O system, over a large range of salt concentrations and up to temperatures of 250°C.     

Water-enhanced dynamic recrystallization and solution transfer in experimentally deformed carnallite

Cylindrical samples of polycrystalline carnallite (KMgCl₃, 6H₂O) were deformed in a triaxial apparatus at 60°C, at confining pressures between 0.1 and 31 MPa and at strain rates between 10⁻⁴ and 10⁻⁸ s⁻¹. In a number of cases, small amounts of saturated carnallite brine were added. Samples without added brine deform by intracrystalline slip, mechanical twinning, cracking, and by frictional sliding on crack surfaces. Stress-strain curves of these samples are strongly dependent on confining pressure. Addition of brine has a dramatic effect on both microstructural development and mechanical properties. Grain-boundary migration is strongly enhanced. At lower strain rates, additional intracrystalline effects start to appear, together with the onset of solution transfer. Rapid compaction in samples deformed with added brine causes high fluid pressures to develop. At higher strain rates addition of brine results in a decrease of the flow stress by a factor of two. This weakening will increase even further at strain rates below about 10⁻⁹ s⁻¹, when solution transfer becomes rate controlling. It is argued that deformation of carnallite in nature is adequately described by the flow law found for samples deformed with added brine.     

Equilibria in the system Au–Ag–S–fluid: computed and experimental data

We present results of computations on the interaction of solid-phase electrum–argentite–pyrite (weight ratios 210⁻⁵/ 210⁻³/1 and 210⁻⁵/410⁻²/1) association with Cl-containing aqueous moderately acid solutions (0.5m NaCl, pH = 3.08) at 300 °C and 500 bars. These data are a physicochemical basis for predicting the geochemical behavior of Au and Ag during the hydrothermal-metasomatic transformation of Au-Ag-pyrite. We also propose a technique of study of this process based on the phase equilibria of the subsystem Au–Ag–S with the aqueous solution at different liquid/solid (l/s) ratios, with the use of new graphic diagrams. The relationship of the composition of the solid-phase association with l/s ratio in real boundary conditions (Au = 17 ppm, mAu/mAg = 10^–3.57–10^–2.28) is shown. The maximum l/s values for complete leaching of gold and silver (l/smax = 200–800) are estimated. It has been established that argentite is the first to dissolve when mAu/mAg(s) > mAu/mAg(sol), and electrum, when mAu/mAg(s) < mAu/mAg(sol).

The experimental results showed that at 300 °C, the conversion of electrum (N_Au = 300‰) nonequilibrated with pyrite into an Au-richer form (N_Au = 730‰) and argentite follows an intricate kinetic scheme. Using the Pilling-Bedwords kinetic equation for processing data yielded the process rate constant K = 2.8(±0.5)10⁻⁵ g²cm⁻⁴day⁻¹. With this equation, the time of the complete conversion of 200 μm thick flat gold grains is 604 days. These data evidence a significant role of kinetic factors in hydrothermal-metasomatic processes involving native gold, which requires combination of thermodynamic and kinetic approaches on the construction of geologo-genetic models for hydrothermal sulfide formation.     

The effect of Ca-Tschermaks component on trace element partitioning between clinopyroxene and silicate melt

We have studied the influence of Ca-Tschermaks (Calcium Tschermaks or CaTs) content of clinopyroxene on the partitioning of trace elements between this phase and silicate melt at fixed temperature and pressure. Ion probe analyses of experiments carried out in the system Na₂O–CaO–MgO–Al₂O₃–SiO₂, at 0.1 MPa and 1218°C, produced crystal-melt partition coefficients (D) of 36 trace elements (Li, Cl, Sc, Ti, V, Cr, Fe, Co, Ge, Sr, Y, Zr, Nb, Mo, Ru, Rh, In, Sn, Sb, Ba, La, Ce, Pr, Nd, Sm, Eu, Gd, Tb, Dy, Ho, Er, Tm, Yb, Lu, Hf, Ta and W), for clinopyroxene compositions between 10 and 32 mol% CaTs. Partition coefficients for 2+ to 5+ cations show, for each charge, a near parabolic dependence of log D on ionic radius of the substituting cation, for partitioning into both the M1 and M2 sites of clinopyroxene. Fitting the results to the elastic strain model of Blundy and Wood [Blundy, J.D., Wood, B.J., 1994. Prediction of crystal-melt partition coefficients from elastic moduli. Nature 372, 452–454] we obtain results for the strain-free partition coefficients of theoretical cations (D₀), with site radius r₀, and for the site's Young's Modulus (E).

In agreement with earlier data our results show that increasing ^ivAl concentration in cpx is matched by increasing D, E_M1, E_M2 and D₀ for tri-, tetra- and pentavalent cations. The degree of fractionation between chemically similar elements (i.e. Ta/Nb, Zr/Hf) also increases. In contrast, D values for mono-, di- and hexavalent cations decrease with increasing ^ivAl in the cpx. The large suite of trace elements used has allowed us to study the effects of cation charge on D₀, r₀ and E. We have found that D₀ and r₀ decrease with increasing cation charge, e.g. r₀=0.66 Å for 4+ cations and 0.59 Å for 5+ cations substituting into M1. Values of E_M1 and E_M2 increase with cation charge as well as with increasing ^ivAl content. The increase in E_M2 is linear and close to the trend set by Hazen and Finger [Hazen, R.M., Finger, L.W., 1979. Bulk modulus-volume relationship for cation–anion polyhedra. J. Geophys. Res. 84 (10) 6723–6728] for oxides. E_M1 values are much higher and do not fit the trend predicted by the Hazen and Finger relationship.     

Conditions of pocket formation in the Oktyabrskaya tourmaline-rich gem pegmatite (the Malkhan field, Central Transbaikalia, Russia)

Coexisting melt (MI), fluid-melt (FMI) and fluid (FI) inclusions in quartz from the Oktaybrskaya pegmatite, central Transbaikalia, have been studied and the thermodynamic modeling of PVTX-properties of aqueous orthoboric-acid fluids has been carried out to define the conditions of pocket formation. At room temperature, FMI in early pocket quartz and in quartz from the coarse-grained quartz–oligoclase host pegmatite contain crystalline aggregates and an orthoboric-acid fluid. The portion of FMI in inclusion assemblages decreases and the volume of fluid in inclusions increases from the early to the late growth zones in the pocket quartz. No FMI have been found in the late growth zones. Significant variations of solid/fluid ratios in the neighboring FMI result from heterogeneous entrapment of coexisting melts and fluids by a host mineral. Raman spectroscopy, SEM EDS and EMPA indicate that the crystalline aggregates in FMI are dominated by mica minerals of the boron-rich muscovite–nanpingite CsAl₂[AlSi₃O₁₀](OH,F)₂ series as well as lepidolite. Topaz, quartz, potassium feldspar and several unidentified minerals occur in much lower amounts. Fluid isolations in FMI and FI have similar total salinity (4–8 wt.% NaCl eq.) and H₃BO₃ contents (12–16 wt.%). The melt inclusions in host-pegmatite quartz homogenize at 570–600 °C. The silicate crystalline aggregates in large inclusions in pocket quartz completely melt at 615 °C. However, even after those inclusions were significantly overheated at 650±10 °C and 2.5 kbar during 24 h they remained non-homogeneous and displayed two types: (i) glass+unmelted crystals and (ii) fluid+glass. The FMI glasses contain 1.94–2.73 wt.% F, 2.51 wt.% B₂O₃, 3.64–5.20 wt.% Cs₂O, 0.54 wt.% Li₂O, 0.57 wt.% Ta₂O₅, 0.10 wt.% Nb₂O₅, 0.12 wt.% BeO. The H₂O content of the glass could exceed 12 wt.%. Such compositions suggest that the residual melts of the latest magmatic stage were strongly enriched in H₂O, B, F, Cs and contained elevated concentrations of Li, Be, Ta, and Nb. FMI microthermometry showed that those melts could have crystallized at 615–550 °C.

Crystallization of quartz–feldspar pegmatite matrix leads to the formation of H₂O-, B- and F-enriched residual melts and associated fluids (prototypes of pockets). Fluids of different compositions and residual melts of different liquidus–solidus P–T-conditions would form pockets with various internal fluid pressures. During crystallization, those melts release more aqueous fluids resulting in a further increase of the fluid pressure in pockets. A significant overpressure and a possible pressure gradient between the neighboring pockets would induce fracturing of pockets and “fluid explosions”. The fracturing commonly results in the crushing of pocket walls, formation of new fractures connecting adjacent pockets, heterogenization and mixing of pocket fluids. Such newly formed fluids would interact with a primary pegmatite matrix along the fractures and cause autometasomatic alteration, recrystallization, leaching and formation of “primary–secondary” pockets.     

钨的气态迁移与岩浆-热液成矿作用:实验研究及其成矿学意义

低密度的水热蒸气和超临界似气流体广泛存在于中地壳至地球表面的各种地质环境中, 是成矿金属搬运和富集的重要介质。火山喷气凝结水、火山结壳和升华物、矿床的流体包裹体气相中均存在具有地质意义的W含量或含钨矿物, 表明W同样可在含水气相中溶解和迁移。本文在350 ℃~400 ℃和压力为60~200 bar的实验条件下, 测定了WO₃-H₂O体系中W在水蒸气和似气流体中的溶解度, 考察了水蒸气压力对W溶解度的影响。结果显示, W在水蒸气中的逸度(或含量)远高于依据无水体系中固体WO₃挥发性数据计算的蒸气压力, 证明气态溶质W与溶剂水蒸气之间存在促进W溶解的水合作用。经热力学方法分析, 认为可能形成了WO₃·nH₂O(g)形式的水合气体物种, 其水合数n在350 ℃、370 ℃和400 ℃时分别为1.4、1.6和2.9。因此WO₃·3H₂O(g)或H₂WO₄·2H₂O(g)及H₆WO₆(g)在温压较高的岩浆-热液或气成-热液成矿环境中(如斑岩系统)对W的气态迁移和浓集可能具有重要作用, 而在温压较低的水热蒸气条件下, W的迁移形式可能以水合数较小的WO₃·H₂O(g)(或H₂WO₄)和WO₃·2H₂O(g)(或H₂WO₄·H₂O)物种为主, 其含量或比例随水蒸气的压力而改变。某些斑岩型和脉型钨(钼)矿床常存在富气体包裹体, 伴随酸性岩浆结晶出溶的以低盐度含水蒸气占优势的岩浆流体对斑岩系统中W、Mo在高温阶段的气态迁移和矿质在花岗岩体顶部和上覆岩层的聚集具有重要意义, 之后蒸气冷凝可产生高盐度的含矿卤水或与渗流地下水混合形成低-中等盐度的成矿流体, 流体的减压沸腾(相分离)和对围岩的交代蚀变导致W、Mo等金属在不同阶段和构造-岩性部位沉淀富集。     

Gold partitioning in melt-vapor-brine systems

We used laser-ablation inductively coupled plasma mass spectrometry to measure the solubility of gold in synthetic sulfur-free vapor and brine fluid inclusions in a vapor + brine + haplogranite + magnetite + gold metal assemblage. Experiments were conducted at 800°C, oxygen fugacity buffered at Ni-NiO (NNO), and pressures ranging from 110 to 145 MPa. The wt% NaCl eq. of vapor increases from 2.3 to 19 and that of brine decreases from 57 to 35 with increasing pressure. The composition of the vapors and brines are dominated by NaCl + KCl + FeCl₂ + H₂O. Gold concentrations in vapor and brine decrease from 36 to 5 and 50 to 28 μg/g, respectively, and the calculated vapor:brine partition coefficients for gold decrease from 0.72 to 0.17 as pressure decreases from 145 to 110 MPa. These data are consistent with the thermodynamic boundary condition that the concentration of gold in the vapor and brine must approach a common value as the critical pressure is approached along the 800°C isotherm in the NaCl-KCl-FeCl₂-HCl-H₂O system.We use the equilibrium constant for gold dissolution as AuOH⁰, extrapolated from lower temperature and overlapping pressure range, to calculate expected concentrations of AuOH⁰ in our experimental vapors. These calculations suggest that a significant quantity of gold in our experimental vapors is present as a non-hydroxide species. Possible chloridogold(I) species are hypothesized based on the positively correlated gold and chloride concentrations in our experimental vapors. The absolute concentration of gold in our synthetic vapor, brine, and melt and calculated mass partition coefficients for gold between these physicochemically distinct magmatic phases suggests that gold solubility in aqueous fluids is a function of aqueous phase salinity, specifically total chloride concentration, at magmatic conditions. However, though we highlight here the effect of salinity, the combination of our data with data sets from lower temperatures evinces a significant decrease in gold solubility as temperature drops from 800°C to 600°C. This decrease in solubility has implications for gold deposition from ascending magmatic fluids.