Rare earth element (REE) deposits are commonly associated with carbonatites and (per)alkaline rocks where hydrothermal magmatic fluids can play a significant role in REE mobilization and deposition [1]. Thermodynamic modeling permits predicting the evolution of ore-forming fluids and can be used to test different controls on hydrothermal REE mobility including temperature, pressure, the solubility of REE minerals, aqueous REE speciation and pH evolution associated with fluid-rock interaction. Previous modeling studies either focused on REE fluoride/chloride complexation in acidic aqueous fluids [2] or near neutral/alkaline fluids associated with calcite vein formation [3]. Such models were also applied to interpret field observations in REE deposits Bayan Obo in China and Bear Lodge in Wyoming [3,4]. Recent hydrothermal calcite-fluid REE partitioning experiments provide new data to simulate the solubility of REE in calcite, REE carbonates/fluorocarbonates at high temperatures [5, 6].
We studied the competing effects controlling the mobility of REE in hydrothermal fluids between 100 and 400 °C at 500 bar. Speciation calculations were carried out in the Ca-F-CO2-Na-Cl-H2O system using the GEMS code package [7]. The properties of minerals and aqueous species were taken from the MINES thermodynamic database [3,5]. The Gallinas Mountains hydrothermal REE deposit in New Mexico was used as a field analogue to compare our models with the formation of calcite-fluorite veins hosting bastnäsite. Previous fluid inclusion studies hypothesized that the REE were transported as fluoride complexes [8] but more recent modeling studies have shown that fluoride essentially acts as a depositing ligand [2]. Here we show more detailed simulations predicting the stability of fluorite, calcite and REE minerals relevant to ore-forming processes in carbonatites and alkaline systems.
[1] Gysi et al. (2016), Econ. Geol. 111, 1241-1276; [2] Migdisov and Williams-Jones (2014), Mineral. Deposita 49, 987-997. [3] Perry and Gysi (2018), Geofluids; [4] Liu et al. (2020), Minerals 10, 495; [5] Perry and Gysi (2020), Geochim. Cosmochim. Acta 286, 177-197; [6] Gysi and Williams-Jones (2015) Chem. Geol. 392, 87-101;[7] Kulik et al. (2013), Computat. Geosci. 17, 1-24; [8] Williams-Jones et al. (2000), Econ. Geol. 95, 327-341
more »
« less
Constraints on deep, CO2-rich degassing at arc volcanoes from solubility experiments on hydrous basaltic andesite of Pavlof Volcano, Alaska Peninsula, at 300 to 1200 MPa
Abstract The solubility of CO2 in hydrous basaltic andesite was examined in fO2-controlled experiments at a temperature of 1125 °C and pressures between 310–1200 MPa. Concentrations of dissolved H2O and CO2 in experimental glasses were determined by ion microprobe calibrated on a subset of run glasses analyzed by high-temperature vacuum manometry. Assuming that the solubility of H2O in mafic melt is relatively well known, estimates of XH2Ofluid and PH2Ofluid in the saturating fluid were modeled, and by difference, values for XCO2fluid and PCO2fluid were obtained (XCO2 ~0.5–0.9); fCO2 could be then calculated from the fluid composition, temperature, and pressure. Dissolved H2O over a range of 2.3–5.5 wt% had no unequivocal influence on the dissolution of CO2 at the pressures and fluid compositions examined. For these H2O concentrations, dissolved CO2 increases with fCO2 following an empirical power-law relation: dissolved CO2 (ppmw) = 14.9−3.5+4.5[fCO2 (MPa)]0.7±0.03. The highest-pressure results plot farthest from this equation but are within its 1 standard-error uncertainty envelope. We compare our experimental data with three recent CO2-H2O solubility models: Papale et al. (2006); Iacono-Marziano et al. (2012); and Ghiorso and Gualda (2015). The Papale et al. (2006) and Iacono-Marizano et al. (2012) models give similar results, both over-predicting the solubility of CO2 in a melt of the Pavlof basaltic andesite composition across the fCO2 range, whereas the Ghiorso and Gualda (2015) model under-predicts CO2 solubility. All three solubility models would indicate a strong enhancement of CO2 solubility with increasing dissolved H2O not apparent in our results. We also examine our results in the context of previous high-pressure CO2 solubility experiments on basaltic melts. Dissolved CO2 correlates positively with mole fraction (Na+K+Ca)/Al across a compositional spectrum of trachybasalt-alkali basalt-tholeiite-icelandite-basaltic andesite. Shortcomings of current solubility models for a widespread arc magma type indicate that our understanding of degassing in the deep crust and uppermost mantle remains semi-quantitative. Experimental studies systematically varying concentrations of melt components (Mg, Ca, Na, K, Al, Si) may be necessary to identify solubility reactions, quantify their equilibrium constants, and thereby build an accurate and generally applicable solubility model.
more »
« less
- Award ID(s):
- 1664308
- NSF-PAR ID:
- 10312562
- Date Published:
- Journal Name:
- American Mineralogist
- Volume:
- 106
- Issue:
- 5
- ISSN:
- 0003-004X
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
more » « less
Abstract Plagioclase microlites in a magma nucleate and grow in response to melt supersaturation (Δ
ϕ plag ). The resultant frozen plagioclase crystal size distribution (CSD) preserves the history of decompression pathways (dP/dt ). SNGPlag is a numerical model that calculates the equilibrium composition of a decompressing magma and nucleates and grows plagioclase in response to an imposed Δϕ plag . Here, we test a new version of SNGPlag calibrated for use with basaltic andesite magmas and modeldP/dt for the ca. 12.6 ka Curacautín eruption of Llaima volcano, Chile. Instantaneous nucleation (N plag ) and growth (G plag ) rates of plagioclase were computed using the experimental results of Shea and Hammer (J Volcanol Geotherm Res 260:127–145, 10.1016/j.jvolgeores.2013.04.018, 2013) and used for SNGPlag modeling of basaltic andesite composition. MaximumN plag of 6.1 × 105 cm h−1is achieved at a Δϕ plag of 44% and the maximumG plag of 27.4 μm h−1is achieved at a Δϕ plag of 29%. Our modeled logdP/dt avg range from 2.69 ± 0.09 to 6.89 ± 0.96 MPa h−1(1σ) with an average duration of decompression from 0.87 ± 0.25 to 16.13 ± 0.29 h assuming a starting pressureP i of 110–150 MPa. These rates are similar to those derived from mafic decompression experiments for other explosive eruptions. Using assumptions for lithostatic pressure gradients (dP/dz ), we calculate ascent rates of < 1–6 m s−1. We conducted a second set of Monte Carlo simulations usingP i of 15–30 MPa to investigate the influence of shallower decompression, resulting in logdP/dt avg from 2.86 ± 0.49 to 6.00 ± 0.86 MPa h−1. ThedP/dt modeled here is two orders of magnitude lower than those calculated by Valdivia et al. (Bull Volcanol, 10.1007/s00445-021-01514-8, 2022) for the same eruption using a bubble number density meter, and suggests homogeneous nucleation raisesdP/dt by orders of magnitude in the shallow conduit. Our modeling further supports the rapid-ascent hypothesis for driving highly explosive mafic eruptions. -
Experimental determination of quartz solubility in H2O-CaCl2 solutions at 600–900 °C and 0.6–1.4 GPamore » « less
Abstract Fluid-mediated calcium metasomatism is often associated with strong silica mobility and the presence of chlorides in solution. To help quantify mass transfer at lower crustal and upper mantle conditions, we measured quartz solubility in H2O-CaCl2 solutions at 0.6–1.4 GPa, 600–900 °C, and salt concentrations to 50 mol%. Solubility was determined by weight loss of single-crystals using hydrothermal piston-cylinder methods. All experiments were conducted at salinity lower than salt saturation. Quartz solubility declines exponentially with added CaCl2 at all conditions investigated, with no evidence for complexing between silica and Ca. The decline in solubility is similar to that in H2O-CO2 but substantially greater than that in H2O-NaCl at the same pressure and temperature. At each temperature, quartz solubility at low salinity (XCaCl2 < 0.1) depends strongly on pressure, whereas at higher XCaCl2 it is nearly pressure independent. This behavior is consistent with a transition from an aqueous solvent to a molten salt near XCaCl2 ~0.1. The solubility data were used to develop a thermodynamic model of H2O-CaCl2 fluids. Assuming ideal molten-salt behavior and utilizing previous models for polymerization of hydrous silica, we derived values for the activity of H2O (aH2O), and for the CaCl2 dissociation factor (α), which may vary from 0 (fully associated) to 2 (fully dissociated). The model accurately reproduces our data along with those of previous work and implies that, at conditions of this study, CaCl2 is largely associated (<0.2) at H2O density <0.85 g/cm3. Dissociation rises isothermally with increasing density, reaching ~1.4 at 600 °C, 1.4 GPa. The variation in silica molality with aH2O in H2O-CaCl2 is nearly identical to that in H2O-CO2 solutions at 800 °C and 1.0 GPa, consistent with the absence of Ca-silicate complexing. The results suggest that the ionization state of the salt solution is an important determinant of aH2O, and that H2O-CaCl2 fluids exhibit nearly ideal molecular mixing over a wider range of conditions than implied by previous modeling. The new data help interpret natural examples of large-scale Ca-metasomatism in a wide range of lower crustal and upper mantle settings.
-
more » « less
Abstract Many lines of evidence from high P–T experiments, thermodynamic models, and natural observations suggest that slab-derived aqueous fluids, which flux mantle wedges contain variable amounts of dissolved carbon. However, constraints on the effects of H2O–CO2 fluids on mantle melting, particularly at mantle wedge P–T conditions, are limited. Here, we present new piston cylinder experiments on fertile and depleted peridotite compositions with 3.5 wt.% H2O and XCO2 [= molar CO2 / (CO2 + H2O)] of 0.04–0.17. Experiments were performed at 2–3 GPa and 1350°C to assess how temperature, peridotite fertility, and XCO2 of slab-derived fluid affects partial melting in mantle wedges. All experiments produce olivine + orthopyroxene +7 to 41 wt.% partial melt. Our new data, along with previous lower temperature data, show that as mantle wedge temperature increases, primary melts become richer in SiO2, FeO*, and MgO and poorer CaO, Al2O3, and alkalis when influenced by H2O–CO2 fluids. At constant P–T and bulk H2O content, the extent of melting in the mantle wedge is largely controlled by peridotite fertility and XCO2 of slab-fluid. High XCO2 depleted compositions generate ~7 wt.% melt, whereas, at identical P–T, low XCO2 fertile compositions generate ~30 to 40 wt.% melt. Additionally, peridotite fertility and XCO2 have significant effects on peridotite partial melt compositions. At a constant P–T–XCO2, fertile peridotites generate melts richer in CaO and Al2O3 and poorer in SiO2, MgO + FeO, and alkalis. Similar to previous experimental studies, at a constant P–T fertility condition, as XCO2 increases, SiO2 and CaO of melts systematically decrease and increase, respectively. Such distinctive effects of oxidized form of dissolved carbon on peridotite partial melt compositions are not observed if the carbon-bearing fluid is reduced, such as CH4-bearing. Considering the large effect of XCO2 on melt SiO2 and CaO concentrations and the relatively oxidized nature of arc magmas, we compare the SiO2/CaO of our experimental melts and melts from previous peridotite + H2O ± CO2 studies to the SiO2/CaO systematics of primitive arc basalts and ultra-calcic, silica-undersaturated arc melt inclusions. From this comparison, we demonstrate that across most P–T–fertility conditions predicted for mantle wedges, partial melts from bulk compositions with XCO2 ≥ 0.11 have lower SiO2/CaO than all primitive arc melts found globally, even when correcting for olivine fractionation, whereas partial melts from bulk compositions with XCO2 = 0.04 overlap the lower end of the SiO2/CaO field defined by natural data. These results suggest that the upper XCO2 limit of slab-fluids influencing primary arc magma formation is 0.04 < XCO2 < 0.11, and this upper limit is likely to apply globally. Lastly, we show that the anomalous SiO2/CaO and CaO/Al2O3 signatures observed in ultra-calcic arc melt inclusions can be reproduced by partial melting of either CO2-bearing hydrous fertile and depleted peridotites with 0 < XCO2 < 0.11 at 2–3 GPa, or from nominally CO2-free hydrous fertile peridotites at P > 3 GPa.
-
Understanding the viscosity of mantle-derived magmas is needed to model their migration mechanisms and ascent rate from the source rock to the surface. High pressure–temperature experimental data are now available on the viscosity of synthetic melts, pure carbonatitic to carbonate–silicate compositions, anhydrous basalts, dacites and rhyolites. However, the viscosity of volatile-bearing melilititic melts, among the most plausible carriers of deep carbon, has not been investigated. In this study, we experimentally determined the viscosity of synthetic liquids with ~31 and ~39 wt% SiO2, 1.60 and 1.42 wt% CO2 and 5.7 and 1 wt% H2O, respectively, at pressures from 1 to 4.7 GPa and temperatures between 1265 and 1755 °C, using the falling-sphere technique combined with in situ X-ray radiography. Our results show viscosities between 0.1044 and 2.1221 Pa·s, with a clear dependence on temperature and SiO2 content. The atomic structure of both melt compositions was also determined at high pressure and temperature, using in situ multi-angle energy-dispersive X-ray diffraction supported by ex situ microFTIR and microRaman spectroscopic measurements. Our results yield evidence that the T–T and T–O (T = Si,Al) interatomic distances of ultrabasic melts are higher than those for basaltic melts known from similar recent studies. Based on our experimental data, melilititic melts are expected to migrate at a rate ~from 2 to 57 km·yr−1 in the present-day or the Archaean mantle, respectively.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.2138/am-2021-7531
