On the stability of magmatic cordierite and new thermobarometric equations for cordierite-saturated liquids

In this work, we have reviewed a large compositional dataset (571 analyses) for natural and experimental glasses to understand the physico-chemical and compositional conditions of magmatic cordierite crystallization. Cordierite crystallizes in peraluminous liquids (A/CNK ≥1) at temperatures ≥750 °C, pressures ≤700 MPa, variable H₂O activity (0.1–1.0) and relatively low fO₂ conditions (≤NNO ? 0.5). In addition to A/CNK ratio ≥1, a required condition for cordierite crystallization is a Si + Al cation value of the rhyolite liquid of 4 p8O (i.e. calculated on the 8 oxygen anhydrous basis), which is consistent with low Fe³⁺ contents and the absence or low content of non-bridging oxygens (NBO). This geochemical condition is strongly supported by the rare, if not unique, structure of cordierite where the tetrahedral framework is composed almost exclusively of Si and Al cations the sum of which is equal to 4 p8O [i.e. (Mg,Fe)_8/9Al_16/9Si_20/9O₈], indicating that aluminium (and cordierite) saturation is limited by rhyolite liquids with Al = 4 ? Si. Indeed, synthetic or natural systems with Al > 4 ? Si always show metastable glass-in-glass separation or crystallization of refractory minerals such as corundum (Al_16/3O₈) and aluminosilicates (Al_16/5Si_8/5O₈). Multivariate regression analyses of literature data for experimental glasses coexisting with magmatic cordierite produced two empirical equations to independently calculate the T (±13 °C; ME, maximum error = 29 °C) and P (±16 %; ME% = 27 %) conditions of cordierite saturation. The greatest influence on the two equations is exerted by H₂O_melt and Al concentrations, respectively. Testing of these equations with other thermobarometric constraints (e.g. feldspar-liquid, GASP, Grt–Bt and Grt–Crd equilibria) and thermodynamic models (NCKFMASHTO and NCKFMASH systems) was successfully performed for Crd-bearing rhyolites and residual enclaves from San Vincenzo (Tuscany, Italy), Morococala Field (Bolivia) and El Hoyazo (Spain). The reliability of each calculated P–T pair was graphically evaluated using the minimum and maximum P–T–H₂O relationships for peraluminous rhyolite liquids modified after the metaluminous relationships in this work. Both P–T calculations and checking can be easily performed with the attached user-friendly spreadsheet (i.e. Crd-sat_TB).