skip to main content


Title: Stability of sodalite relative to nepheline in NaCl–H2O brines at 750 °C: Implications for hydrothermal formation of sodalite
ABSTRACT Formation of the feldspathoid sodalite (Na6Al6Si6O24·2NaCl) by reaction of nepheline (NaAlSiO4) with NaCl-bearing brines was investigated at 3 and 6 kbar and at a constant temperature of 750 °C to determine the brine concentration at which sodalite forms with variation in pressure. The reaction boundary was located by reaction-reversal experiments in the system NaAlSiO4–NaCl–H2O at a brine concentration of 0.16 ± 0.08 XNaCl [= molar ratio NaCl/(NaCl + H2O)] at 3 kbar and at a brine concentration of 0.35 ± 0.03 XNaCl at 6 kbar. Characterization of the sodalite using both X-ray diffraction and infrared spectroscopy after treatment in these brines indicated no obvious evidence of water or hydroxyl incorporation into the cage structure of sodalite. The data from this study were combined with earlier results by Wellman (1970) and Sharp et al. (1989) at lower (1–1.5 kbar) and higher (7–8 kbar) pressures, respectively, on sodalite formation from nepheline and NaCl which models as a concave-down curve in XNaCl – P space. In general, sodalite buffers the concentration of neutral aqueous NaCl° in the brine to relatively low values at P < 4 kbar, but NaCl° increases rapidly at higher pressures. Thermochemical modeling of these data was done to determine the activity of the aqueous NaCl° relative to a 1 molal (m) standard state, demonstrating very low activities (<0.2 m, or 1.2 wt.%) of NaCl° at 3 kbar and lower, but rising to relatively high activities (>20 m, or 54 wt.%) of NaCl° at 6 kbar or higher. The results from this study place constraints on the concentration of NaCl° in brines coexisting with nepheline and sodalite and, because of the relative insensitivity of this reaction to temperature, can provide a convenient geobarometer for those localities where the fluid compositions that formed nepheline and sodalite can be determined independently.  more » « less
Award ID(s):
1725053
NSF-PAR ID:
10162732
Author(s) / Creator(s):
;
Date Published:
Journal Name:
The Canadian Mineralogist
Volume:
58
Issue:
1
ISSN:
0008-4476
Page Range / eLocation ID:
3 to 18
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. The granitic water-saturated solidus (G-WSS) is the lower temperature limit of magmatic mineral crystallization. The accepted water-saturated solidus for granitic compositions was largely determined >60 years ago1. More recent advances in experimental petrology, improved analytical techniques, and recent observations that granitic systems can remain active or spend a significant proportion of their lives at conditions below the traditional G-WSS2–5 necessitate a careful experimental investigation of the near-solidus regions of granitic systems. Natural and synthetic starting materials were melted at 10 kbar and 900°C with 48 wt% H2O to produce hydrous glasses for subsequent experiments at lower PT conditions used to locate the G-WSS. We performed crystallization experiments and melting experiments at temperatures ranging from 575 to 800°C and 1, 6, 8, and 10 kbar on 12 granitoid compositions. First, we ran a series of isothermal crystallization experiments along each isobar at progressively lower temperatures until runs completely crystallized to identify apparent solidus temperatures. Geochemical analyses of quenched glass compositions demonstrate that progressive crystallization drives all starting compositions towards silica-rich, water-saturated rhyolitic/granitic melts (e.g., ~7578 wt% SiO2). After identifying the apparent solidus temperatures at which the various compositions crystallized, we then ran series of reversal-type melting experiments. With the goal of producing rocks with hydrous equilibrium microstructures, we crystallized compositions at temperatures ~10°C below the apparent solidus identified in crystallization experiments, and then heated isobarically to conditions that produced ~20% melt during the crystallization experiments. Importantly, crystallization experiments and heating experiments at the same PT conditions produced similar proportions of melt, crystals, and vapor. A time-series of experiments 230 days at PT conditions previously identified to produce ~10% to 20% melt did not reveal any kinetic effects on melt crystallization. Experiments at 6 to 10 kbar crystallized/melted at temperatures close to the published G-WSS. However, at lower pressures where the published G-WSS is strongly curved in PT space, all compositions investigated contained melt to temperatures ~75 to 100°C below the accepted G-WSS. The similarity of crystallization temperatures for the higher-pressure experiments to previously published results, similar phase proportions in melting and crystallization experiments, and the lack of kinetic effects on crystallization collectively suggest that our lower pressure constraints on the G-WSS are accurate. The new experimental results demonstrating that the lower-pressure G-WSS is significantly lower than unanimously accepted estimates will help us to better understand the storage conditions, evolution, and potential for eruption in mid- to upper-crustal silicic magmatic systems. (1) Tuttle, O.; Bowen, N. Origin of Granite in the Light of Experimental Studies in the System NaAlSi3O8–KAlSi3O8–SiO2–H2O; Geological Society of America Memoirs; Geological Society of America, 1958; Vol. 74. https://doi.org/10.1130/MEM74. (2) Rubin, A. E.; Cooper, K. M.; Till, C. B.; Kent, A. J. R.; Costa, F.; Bose, M.; Gravley, D.; Deering, C.; Cole, J. Rapid Cooling and Cold Storage in a Silicic Magma Reservoir Recorded in Individual Crystals. Science 2017, 356 (6343), 1154–1156. https://doi.org/10.1126/science.aam8720. (3) Andersen, N. L.; Jicha, B. R.; Singer, B. S.; Hildreth, W. Incremental Heating of Bishop Tuff Sanidine Reveals Preeruptive Radiogenic Ar and Rapid Remobilization from Cold Storage. Proceedings of the National Academy of Sciences 2017, 114 (47), 12407–12412. https://doi.org/10.1073/pnas.1709581114. (4) Ackerson, M. R.; Mysen, B. O.; Tailby, N. D.; Watson, E. B. Low-Temperature Crystallization of Granites and the Implications for Crustal Magmatism. Nature 2018, 559 (7712), 94–97. https://doi.org/10.1038/s41586-018-0264-2. (5) Glazner, A. F.; Bartley, J. M.; Coleman, D. S.; Lindgren, K. Aplite Diking and Infiltration: A Differentiation Mechanism Restricted to Plutonic Rocks. Contributions to Mineralogy and Petrology 2020, 175 (4). https://doi.org/10.1007/s00410-020-01677-1. 
    more » « less
  2. Abstract The composition of clinopyroxene and clinopyroxene-liquid (Cpx-Liq) pairs are frequently used to calculate crystallization/equilibration pressures in igneous systems. While canonical uncertainties are often assigned to calculated pressures based on fits to calibration or test datasets, the sources of these uncertainties (and thus ways to reduce them) have not been rigorously assessed. We show that considerable uncertainties in calculated pressures arise from analytical error associated with Electron Probe Microanalyser (EPMA) measurements of Cpx. Specifically, low X-ray counts during analysis of elements with concentrations <1 wt% resulting from insufficient count times and/or low beam currents yield highly imprecise measurements (1σ errors of 10–40% for Na2O). Low analytical precision propagates into the calculation of pressure-sensitive mineral components such as jadeite. Using Monte Carlo approaches, we demonstrate that elemental variation resulting from analytical precision alone generates pressures spanning ~4 kbar (~15 km) for a single Cpx and ~6 kbar for a single Cpx-Liq pair using popular barometry expressions. In addition, analytical uncertainties in mineral compositions produce highly correlated arrays between pressure and temperature that have been previously attributed to transcrustal magma storage. Before invoking such geological interpretations, a more mundane origin from analytical imprecision must be ruled out. Most importantly, low analytical precision does not just affect the application of barometers to natural systems; it has also affected characterization of Cpx in experimental products used to calibrate and test barometers. The impact of poor precision on each individual measurement is often magnified by the small number of measurements made within experimental charges, meaning that low analytical precision and true variability in mineral compositions have not been sufficiently mediated by averaging multiple EPMA analyses. We compile the number of Cpx measurements performed in N = 307 experiments used to calibrate existing barometers, and N = 490 new experiments, finding ~45% of experiment charges were characterized by ≤5 individual Cpx analyses. Insufficient characterization of the true composition of experimental phases likely accounts for the fact that all Cpx-based barometers exhibit large errors (± 3 kbar) when tested using global experimental datasets. We suggest specific changes to analytical and experimental protocols, such as increased count times and/or higher beam currents when measuring low concentration elements in relatively beam resistant Cpx in experiments and natural samples. We also advocate for increasing the number of analyses per experimental charge, resolving interlaboratory analytical offsets and improving data reporting. Implementing these changes is essential to produce a more robust dataset to calibrate and test the next generation of more precise and accurate Cpx-based barometers. In turn, this will enable more rigorous investigation of magma storage geometries in a variety of tectonic settings (e.g. distinguishing true transcrustal storage vs. storage in discrete reservoirs). 
    more » « less
  3. Abstract

    Using molten‐salt synthetic techniques, NaNbO3(Space groupPbcm; No. 57) was prepared in high purity at a reaction time of 12 hours and a temperature of 900°C. All NaNbO3products were prepared from stoichiometric ratios of Nb2O5and Na2CO3together with the addition of a salt flux introduced at a 10:1 molar ratio of salt to NaNbO3, that is, using the Na2SO4, NaF, NaCl, and NaBr salts. A solid‐state synthesis was performed in the absence of a molten salt to serve as a control. The reaction products were all found to be phase pure through powder X‐ray diffraction, for example, with refined lattice constants ofa = 5.512(5) Å,b = 5.567(3) Å, andc = 15.516(8) Å from the Na2SO4salt reaction. The products were characterized using UV‐Vis diffuse reflectance spectroscopy to have a bandgap size of ~3.5 eV. The particles sizes were analyzed by scanning electron microscopy (SEM) and found to be dependent upon the flux type used, from ~<1 μm to >10 μm in length, with overall surface areas that could be varied from 0.66 m2/g (for NaF) to 1.55 m2/g (for NaBr). Cubic‐shaped particle morphologies were observed for the metal halide salts with the set of exposed (100)/(010)/(001) crystal facets, while a truncated octahedral morphology formed in the sodium sulfate salt reaction with predominantly the set of (110)/(101)/(011) crystal facets. The products were found to be photocatalytically active for hydrogen production under UV‐Vis irradiation, with the aid of a 1 wt% Pt surface cocatalyst. The platinized NaNbO3particles were suspended in an aqueous 20% methanol solution and irradiated by UV‐Vis light (λ > 230 nm). After 6 hours of irradiation, the average total hydrogen production varied with the particle morphologies and sizes, with 753 µmol for Na2SO4, 334 µmol for NaF, 290 µmol for NaCl, 81 µmol for NaBr, and 249 µmol for the solid‐state synthesized NaNbO3. These trends show a clear relationship to particle sizes, with smaller particles showing higher photocatalytic activity in the order of NaF > NaCl > NaBr. Furthermore, the particle morphologies obtained from the Na2SO4flux showed even higher photocatalytic activity, though having a relatively similar overall surface area, owing to the higher activity of the (110) crystal facets. The apparent quantum yield (100 mW/cm2,λ = 230 to 350 nm, pH = 7) was measured to be 3.7% for NaNbO3prepared using the NaF flux, but this was doubled to 6.8% when prepared using the Na2SO4flux. Thus, these results demonstrate the powerful utility of flux synthetic techniques to control particle sizes and to expose higher‐activity crystal facets to boost their photocatalytic activities for molecular hydrogen production.

     
    more » « less
  4. Abstract Scapolites are pervasive rock-forming aluminosilicates that are found in metamorphic, igneous, and hydrothermal environments; nonetheless, the stability field of Cl-rich scapolite is not well constrained. This experimental study investigated two reactions involving Cl-rich scapolite. First, the anhydrous reaction 1 of plagioclase + halite + calcite to form scapolite [modeled as: 3 plagioclase (Ab80An20) + 0.8 NaCl + 0.2 CaCO3 = scapolite (Ma80Me20)] was investigated to determine the effect of the Ca-rich meionite (Me = Ca4Al6Si6O24CO3) component on the Na end-member marialite (Ma = Na4Al3Si9O24Cl). Second, the effect of water on this reaction was investigated using the hydrothermally equivalent reaction 2, H2O + scapolite (Ma80Me20) = 3 plagioclase (Ab80An20) + CaCO3 + liquid, where the liquid is assumed to be a saline-rich hydrous-silicate melt. Experiments were conducted with synthetic phases over the range of 500–1030 °C and 0.4–2.0 GPa. For reaction 1, intermediate composition scapolite shows a wide thermal stability and is stable relative to plagioclase + halite + calcite at temperatures above 750 °C at 0.4 GPa and 760 °C at 2.0 GPa. For reaction 2, intermediate scapolite appears to be quite tolerant of water; it forms at a minimum bulk salinity [XNaCl = molar ratio of NaCl/(NaCl+H2O)] of the brine of approximately 0.2 XNaCl at 830 and 680 °C at pressures of 2.0 and 1.5 GPa, respectively. Based on the study done by Almeida and Jenkins (2017), pure marialite is very intolerant of water when compared to intermediate composition scapolite. Compositional changes in the scapolite and plagioclase were characterized by X-ray diffraction and electron microprobe analysis and found to shift from the nominal bulk compositions to the observed compositions of Ma85Me15 for scapolite and to Ab91An09 for plagioclase. These results were used to model the phase equilibria along the marialitemeionite join in temperature-composition space. This study demonstrates that a small change in the scapolite composition from end-member marialite to Ma85Me15 expands the stability field of marialite significantly, presumably due to the high entropy of mixing in scapolite, as well as increases its tolerance to water. This supports the much more common presence of intermediate scapolites in hydrothermal settings than either end-member meionite or marialite as is widely reported in the literature. 
    more » « less
  5. Abstract The Bishop Tuff (BT), erupted from the Long Valley caldera in California, displays two types of geochemical gradients with temperature: one is related to magma mixing, whereas the other is found in the high-SiO2 rhyolite portion of the Bishop Tuff and is characterized by twofold or lower concentration variations in minor and trace elements that are strongly correlated with temperature. It is proposed that the latter zonation, which preceded phenocryst growth, developed as a result of mineral–melt partitioning between interstitial melt and surrounding crystals in a parental mush, from which variable melt fractions were segregated. To test this hypothesis, trends of increasing vs decreasing element concentrations with temperature (as a proxy for melt fraction), obtained from published data on single-clast pumice samples from the high-SiO2 rhyolite portion of the Bishop Tuff, were used to infer their relative degrees of incompatibility vs compatibility between crystals and melt in the parental mush. Relative compatibility values (RCVi) for all elements i, defined as the concentration slope with temperature divided by average concentration, are shown to be linearly correlated with their respective bulk partition coefficients (bulk Di). Mineral–melt partition coefficients from the literature were used to constrain the average stoichiometry of the crystallization/melting reaction in the parental mush: 32 % quartz + 34 % plagioclase + 31 % K-feldspar + 1·60 % biotite + 0·42 % titanomagnetite + 0·34 % ilmenite + 0·093 % allanite + 0·024 % zircon + 0·025 % apatite = 100 % liquid. The proportions of tectosilicates in the reaction (i.e. location of eutectic) are consistent with depths of melt segregation of ~400–550 MPa and an activity of H2O of ~0·4–0·6. Temperatures of <770–780 °C are constrained by allanite in the reaction. Evidence that a fluid phase was present in the parental mush is seen in the decreasing versus increasing H2O and CO2 contents with temperature in the segregated interstitial melt that formed the high-SiO2 rhyolite portion of the Bishop Tuff. The presence of an excess fluid phase, which strongly partitions CO2 relative to the melt, is required to explain the compatible behavior of CO2, whereas the fluid abundance must have been low to explain the incompatible behavior of H2O. Calculated degassing paths for interstitial melts, which segregated from the parental mush and ascended to shallower depths to grow phenocrysts, match published volatile analyses in quartz-hosted melt inclusions and constrain fluid abundances in the mush to be ≤1 wt%. The source of volatiles in the parental mush, irrespective of whether it formed by crystallization or partial melting, must have been primarily from associated basalts, as granitoid crust is too volatile poor. Approximately twice as much basalt as rhyolite is needed to provide the requisite volatiles. The determination of bulk Di for several elements gives the bulk composition of the parental leucogranitic mush and shows that it is distinct from Mesozoic Sierran arc granitoids, as expected. Collectively, the results from this study provide new constraints for models of the complex, multi-stage processes throughout the Plio-Quaternary, involving both mantle-derived basalt and pre-existing crust, that led to the origin of the parental body to the Bishop Tuff. 
    more » « less