QUILF: A pascal program to assess equilibria among FeMgMnTi oxides, pyroxenes, olivine, and quartz

doi:10.1016/0098-3004(93)90033-2

Computers & Geosciences

Volume 19, Issue 9, October 1993, Pages 1333-1350

https://doi.org/10.1016/0098-3004(93)90033-2 Get rights and content

Abstract

Program QUILF assesses equilibria among Ti-magnetite, ilmenite, augite, pigeonite, orthopyroxene, olivine, and quartz (or subassemblages thereof). Oxide and silicate equilibria are related through the QUIIF equilibrium:

Depending on the assemblage, QUILF can provide information on temperature, pressure, oxygen fugacity, and the activities of SiO₂, TiO₂, and Fe^o at which the phases were last in equilibrium. For many low-variance assemblages, the system is overdetermined; thus quantitative information can be extracted even if one phase is altered or has reequilibrated. QUILF equilibria can reduce the uncertainties in temperature and oxygen fugacity as determined from coexisting ilmenite and Ti-magnetite alone.

QUILF is written in Turbo Pascal for IBM PC and compatibles. The compiled program is approximately 210 kbyte; it also requires two data files that total approximately 10 kbyte. The compositions of the phases first must be projected into seven-component space (CaOMgOMnOFeOFe₂O₃TiO₂SiO₂) before they are used in program QUILF; routines to accomplish this are included in the program, and also are available as separate BASIC programs.

References (20)

D.J. Andersen
Internally consistent solution models for FeMgMnTi oxides
D.J. Andersen et al.
Internally consistent solution models for FeMgMnTi oxides: Part II. FeMgTi oxides and olivine
Am. Mineralogist
(1991)
D.J. Andersen et al.
Internally consistent solution models for FeMgMnTi oxides: FeTi oxides
Am. Mineralogist
(1988)
S.R. Bohlen et al.
Feldspar and oxide thermometry of granulites in the Adirondack Highlands
Contr. Mineralogy and Petrology
(1977)
G.M. Brown
Pyroxenes from the early and middle stages of fractionation of the Skaergaard intrusion, east Greenland
Mineral. Magazine
(1957)
P.M. Davidson
Thermodynamic analysis of quadrilateral pyroxenes. Part I. Derivation of the ternary nonconvergent site-disorder model
Contr. Mineralogy and Petrology
(1985)
P.M. Davidson et al.
Thermodynamic analysis of pyroxene-olivine-quartz equilibria in the system CaO-MgO-FeO-SiO₂
Am. Mineralogist
(1989)
P.M. Davidson et al.
CaFeMg olivines: phase relations and a solution model
Contr. Mineralogy and Petrology
(1984)
B.R. Frost et al.
Equilibria among FeTi oxides, pyroxenes, olivine, and quartz: Part II. Application
Am. Mineralogist
(1992)
B.R. Frost et al.
Fe-Ti oxide silicate equilibria: assemblages with fayalitic olivine
Am. Mineralogist
(1988)

There are more references available in the full text version of this article.

Cited by (653)

On the origin of metallic iron in eucrite breccias: Effects of impact shock and mixing on the surface of (4) Vesta
2024, Icarus
We performed scanning electron microscopy and electron microprobe analysis on the eucritic meteorites Errachidia 002 and Northwest Africa 11911, with a focus on the metallic iron disseminated in them, in order to constrain the processes involved in its formation, and discuss its implication for their parent body collisional history. Errachidia 002 is a polymict eucrite that consists of basaltic and noritic (diogenite) clasts set in a fragmental matrix. Fe-metals in this meteorite contain <0.04 wt% Ni, and 0.06 wt% Co, and they are mostly associated with silica in pyroxene. This suggests formation by decomposition of troilite, which generated S-vapors (S₂) that subsequently reacted with FeO in pyroxene during an impact event, producing pure Fe-metal and silica, with escape of SO₂ from the rock. The presence of late-stage phases in proximity of a few Fe-metal grains, as well as low Fe/Mn coupled with high Fe/Mg in pyroxene indicate a different scenario in which pure Fe-metal in this meteorite formed by reduction from a residual melt. NWA 11911 is also a polymict eucrite composed mainly of unequilibrated and equilibrated eucritic clasts set in a clastic groundmass. Here, metallic iron is either enclosed in the fragmental matrix, or within a vitrophyric impact melt clast. Its Ni and Co contents are ∼2.8 and 0.6 wt%, implying an exogenous origin, namely contamination by an iron-bearing impactor. Our results point out the major role impacts played in controlling the surface features and mineral compositions of early formed bodies in the solar system such as (4) Vesta.
Reduced post-collisional magmas for Tengchong volcanoes, southeastern Tibetan Plateau
2024, Journal of Asian Earth Sciences
Rapidly cooled and well-preserved magnetite and ilmenite in young volcanic rocks provide good samples for constraining the magmatic redox conditions. We report chemical compositions of magnetite and ilmenite from the youngest Tengchong volcanoes at the southeastern margin of the Tibetan Plateau, to constrain the oxygen fugacity and temperatures for the Tengchong magmas. Oxygen fugacity values (△FMQ) relative to the fayalite-magnetite-quartz (FMQ) buffer are −1.8 ∼ −0.7 (averaging −1.1 ± 0.4) for Dayingshan; −1.3 ∼ −0.7（averaging −1.1 ± 0.3） for Maanshan; and −1.3 for Heikongshan. The temperatures for magnetite-ilmenite pairs of Tengchong lavas range from 789 to 1090 °C. The relatively reducing oxygen fugacity values (△FMQ = −1.8 ∼ −0.7) of the Tengchong magmas are significantly lower than those oxidizing arc magmas (△FMQ = +0.5 ∼ +3.5). The reduced magmas at Tengchong may exemplify other Quaternary post-collisional lavas of the Tibetan Plateau. The crystallization temperatures from mineral-liquid equilibrium are 1160–1192 °C for olivine, 1049–1110 °C for plagioclase, 978–1122 °C for clinopyroxene, 789–1090 °C for Fe-Ti oxides, 743–792 °C for zircon. The high Fe-Ti oxide temperature and Ti-in-zircon temperature indicate warm magmatic conditions (>750 °C) prior to eruptions of Tengchong magmas.
Petrogenesis of oxide-apatite mineralization associated with proterozoic anorthosite massifs at Lac Mirepoix, Quebec, Canada: A multi-injection model for Fe-Ti-P mineralization in the Central Grenville Province
2023, Ore Geology Reviews
The Lac Mirepoix Fe-Ti-P mineralization presents several mineralized lenses that outcrop near the border of the 1080 (±2) Ma Vanel and the 1016 (±2) Ma Mattawa Anorthosite massifs, in the Central Grenville Province, Quebec, Canada. The mineralization is (hemo)-ilmenite-dominated, accompanied by magnetite and apatite. It is subdivided in three different zones due to the appearance of different cumulate phases: zone I comprises mainly massive oxide (>70 % hemo-ilmenite ± magnetite) layers hosted in anorthosite. Towards the center of the mineralization (zone II), massive oxide layers are less common whereas apatite-bearing cumulates appear, forming massive nelsonite (50–70 % magnetite ± ilmenite and 25–30 % apatite) and oxide apatite norite (OAN, 15–25 % hemo-ilmenite ± magnetite and 8–20 % apatite). Finally, zone III is marked by the alternance of OAN layers (10–25 m), richer in magnetite, in addition to (hemo)-ilmenite and apatite, and an absence of massive oxide and nelsonite. In-situ trace element analysis of plagioclase, apatite and oxides by laser ablation ICP-MS, reveals cryptic variations related to magma differentiation and multiple injections of ferrodiorite parental magmas of similar composition. In-situ U-Pb dating of zircon from the OAN mineralization itself indicates two different crystallization ages between zone III (1048 ± 8 Ma) and zone I (964 ± 9 Ma), favouring a model of multiple injections of similar Fe-Ti-P-rich melts of similar composition over a period of 80 Mys rather than in-situ crystallization of a single intrusion, which is supported by trace-element geochemistry. The Lac Mirepoix mineralization records the following fractional crystallization sequence of a high-Ti-P magma, residual after the anorthosite formation: first, massive oxides of hemo-ilmenite crystallized (high Ti/Fe) by oxide settling, with primitive compositions, similar to the world-class Lac Tio deposit. The oxide-apatite mineralization (nelsonite and OAN) crystallized from the residual liquid (lower Ti/Fe, evolved compositions) in which magnetite and apatite were more abundant, similar to nearby mineralization in the area (e.g. Lac à l’Orignal Fe-Ti-P deposit). The observed decrease in Ti content of the evolving melt is supported by the liquid line of descent of several ferrodiorite dykes within the host-rocks and mineralization. Finally, we propose a model of multiple injections of residual Fe-Ti-P-rich liquids drained from or filter-pressed within diapirs of plagioclase-rich mushes that were emplaced over 80 Mys apart along the same crustal detachment zone.
Magmatic differentiation within the transcrustal magma reservoir of Develidağ volcano, Central Anatolia (Turkey)
2023, Journal of Asian Earth Sciences
The polygenetic activity of the Develidağ volcano from Miocene to Pliocene (∼5.5 Ma) is represented by basaltic/dacitic calk-alkaline and Na-alkaline magmas on the SE rim of Central Anatolian Volcanic Province (Kayseri, Turkey). Ubiquitous resorption patterns among An-rich (> 75 X_An) plagioclase and complex pyroxene with common microlithic to porphyritic texture indicate the complex evolutionary history of magmas. Common normal zoning trend among the crystallizing phases, decreasing major element oxides, and systematic enrichment of Th, Zr, with LILE against increasing silica among rocks attests to fractional crystallization and crustal assimilation processes. Estimated intensive reservoir parameters based on mineral-only, mineral pair, and mineral-liquid equilibria demonstrate that partial crystallization of magmas spanned a broad P (kbar) / T (°C) gradient over the NNO buffer. Least-square mass balance (39–71% FC), trace element (50% AFC; r = 0.3–0.5), and pMELTS/rhyolite-MELTS (9.23–0.5 kbar, 1227–850 °C) models indicate the transition between ascent and cooling-driven fractionation, involvement of crustal material with the mixing between less evolved and evolved melts along the vertically extensive magmatic reservoir system. Micro- to mesoscale reconstruction of magmatic reservoir processes along the formation of physico-chemically distinct magma series provides a better understanding of eruptive processes and the tectonomagmatic evolution of a transcrustal reservoir system before the transtensional regime at the Central Anatolian Volcanic Province.
Pre-eruptive storage conditions and magmatic evolution of the Bora-Baricha-Tullu Moye volcanic system, Main Ethiopian Rift
2023, Lithos
Bora-Baricha-Tullu Moye is a Late Quaternary volcanic system in the Main Ethiopian Rift, characterised by products of both explosive and effusive volcanic eruptions. The petrological and geochemical characteristics of the volcanic products are investigated using a combination of petrography, major and trace element whole rock analyses and in-situ major element analyses of phenocryst phases, matrix glass and melt inclusions. The bulk rock compositions vary from basalt to peralkaline rhyolite (comendite and pantellerite), and the chemical variability can largely be explained by fractional crystallisation processes with minor crustal assimilation and magma mixing. The dominant mineral phases such as clinopyroxenes and feldspars show a tendency for Fe and Na enrichment respectively from the basalts towards the pantellerites. The comendite and pantellerite deposits show systematic variations towards more evolved glass and mineral composition with the stratigraphy. The combination of thermometry (i.e., clinopyroxene-liquid, feldspar-liquid, olivine-liquid and clinopyroxene-only) and barometry (i.e., clinopyroxene-liquid and clinopyroxene-only) modelling suggests that the basaltic magmas are stored at high temperature (1070–1190 °C) at mid-to-deep-crustal levels (∼7–29 km). The peralkaline rhyolite melts are stored at lower temperature (i.e., 805–900 °C for comendite; 700–765 °C for pantellerite) at shallow crustal levels (∼4 km). The conditions of pre-eruptive storage as recorded in the comendite and pantellerite rocks in combination with stratigraphic constraints, suggests a progressive temporal evolution of the magma reservoirs to cooler storage temperatures.
Contrasting magmatic evolutions of the Three Sister Volcanoes reflect increased heat flow, crustal melting and silicic magmatism in the Central Oregon Cascade Arc
2023, Chemical Geology
The Three Sisters stretch, one of the most magmatically active segments of the Oregon Cascade Arc, includes the Three Sisters volcanoes proper (North, Middle and South) as well as the shields and cinder cones of the Mount Bachelor Chain and McKenzie Pass. The Three Sisters record a magmatic transition from fundamentally basaltic andesitic activity in older North Sister products to more diverse activity at Middle Sister, and unusually silicic activity in a Cascade volcano, including rhyolite lava, at South Sister.
North Sister magmas reflect magmatic evolution of basaltic andesite through crystal fractionation, mostly of olivine, and incorporation of crustal melts. Middle Sister, although mostly composed of basaltic andesite and andesite, produced dacitic magma through assimilation fractional crystallization processes involving incorporation of silicic melts. South Sister developed a subjacent magma chamber under the southern half of the volcano, where andesitic magmas evolved, dominantly through fractional crystallization, to dacite and rhyolite. A prominent compositional gap at South Sister between dacite and rhyolite and depletion of middle and heavy rare earth elements in rhyolite suggest that amphibole was prominently involved in fractionation processes involving separation of silicic magma from a gabbroic or amphibolitic crystal mush. Silicic magma eventually formed a stagnant cap in the system, where it evolved through small scale fractional crystallization. Rhyolite magma was finally heated to elevated temperatures prior to eruption, presumably through influx of less evolved magma into the base of the stagnant silicic cap and was erupted in late rhyolite lava of Rock Mesa and the Devil's Hill Chain.
South Sister lies along strike of the Northwest Rift of Newberry Volcano, which forms the northwestern terminus of the High Lava Plains trend of south-central Oregon. The increased magmatic activity of the Three Sisters stretch and unusual abundance of silicic products, including rhyolite, of South Sister may reflect the arrival and intersection of the High Plains magmatic anomaly with the Cascade Arc.

View all citing articles on Scopus

View full text

QUILF: A pascal program to assess equilibria among FeMgMnTi oxides, pyroxenes, olivine, and quartz

Abstract

Internally consistent solution models for FeMgMnTi oxides

Internally consistent solution models for FeMgMnTi oxides: Part II. FeMgTi oxides and olivine

Am. Mineralogist

Internally consistent solution models for FeMgMnTi oxides: FeTi oxides

Am. Mineralogist

Feldspar and oxide thermometry of granulites in the Adirondack Highlands

Contr. Mineralogy and Petrology

Pyroxenes from the early and middle stages of fractionation of the Skaergaard intrusion, east Greenland

Mineral. Magazine

Thermodynamic analysis of quadrilateral pyroxenes. Part I. Derivation of the ternary nonconvergent site-disorder model

Contr. Mineralogy and Petrology

Thermodynamic analysis of pyroxene-olivine-quartz equilibria in the system CaO-MgO-FeO-SiO2

Am. Mineralogist

CaFeMg olivines: phase relations and a solution model

Contr. Mineralogy and Petrology

Equilibria among FeTi oxides, pyroxenes, olivine, and quartz: Part II. Application

Am. Mineralogist

Fe-Ti oxide silicate equilibria: assemblages with fayalitic olivine

Am. Mineralogist

Thermodynamic analysis of pyroxene-olivine-quartz equilibria in the system CaO-MgO-FeO-SiO₂