G. S. Barmina

COMAGMAT: a Fortran program to model MAGMA differentiation processes

Abstract The purpose of this paper is to present a general model (COMAGMAT) for the calculation of equilibrium temperatures and phase relations at a given extent of crystallization or melting in natural magmatic systems. This model is based on a set of empirical expressions that describe mineral-melt equilibria for major and trace elements in terms of pressure (up to 12 kbar), temperature, and liquid compositions for systems ranging from primitive basalts to dacites. These expressions are in the form of empirically calibrated mineral-melt geothermometers for Olivine, Augite, Pigeonite (Opx), Plagioclase, Ilmenite, and Magnetite used to develop the algorithm simulating multiply saturated magmatic melts. The results of the program are in the form of calculated liquid lines of descent, plus the equilibrium mineral proportions and compositions. The phase equilibria calculations form the core of a model that allows the user to simulate processes ranging from simple isobaric crystallization to in situ differentiation processes resulted from crystal settling, and polybaric fractionation.

An empirical model for the calculation of spinel-melt equilibria in mafic igneous systems at atmospheric pressure: 2. Fe-Ti oxides

Abstract In order to develop models simulating the crystallization of Fe-Ti oxides in natural lavas, we have processed published experimental data on magnetite-melt and ilmenite-melt equilibria. These data include 62 Mt-melt and 75 Ilm-melt pairs at temperatures 1040–1150 °C, oxygen fugacities from IW to NNO+2, and bulk compositions ranging from ferrobasalts to andesites and dacites. Five major cations (Fe3+, Fe2+, Ti4+, Mg2+ and Al3+) were considered for the purpose of describing Fe-Ti oxide saturation as a function of melt composition, temperature and oxygen fugacity at 1 atmosphere pressure. The empirically calibrated mineral-melt expression based on multiple linear regressions is: ln Di = a/T + blog fO2 + c + d1XNa + d2XK + d3XP, where Di represents molar distribution coefficients of the given cations between Mt/Ilm and melt; XNa, XK, and XP are the molar fractions of Na, K, and P in the melt. The empirically calibrated Mt-melt and Ilm-melt equilibria equations allowed us to develop two models for calculating crystallization temperatures of the Fe-Ti oxides in the melts with an accuracy of 10–15 °C, and compositions with an accuracy of 0.5–2 mol%. These models have been integrated into the COMAGMAT-3.5 program, improving our ability to study numerically the effects of temperature and oxygen fugacity on the stability and phase equilibria of Fe-Ti oxides. Application of this approach to the tholeiitic series of Chazhma Sill from Eastern Kamchatka (Russia) indicates oxygen fugacity conditions near NNO + 0.5. Numerical simulation of fractional crystallization of an iron-enriched basaltic andesite parent at these oxidizing conditions accurately reproduces the FeO-SiO2 relations observed in the Chazhma suite.

Simulation of primary phase relations and mineral compositions in the Partridge River intrusion, Duluth Complex, Minnesota: implications for the parent magma composition

In order to describe the composition and crystallinity of the initial (parental) magma of the Partridge River intrusion of the Keweenawan Duluth Complex, and thereby understand the mode of emplacement and solidification of the intrusion, we have applied a numerical simulation technique called geochemical thermometry (Frenkel et al. 1988). The parental magma was a low-alumina, high-Ti-P olivine tholeiite similar to typical Keweenawan low-alumina, high-Ti-P basalts associated with the Duluth Complex and from the nearby Portage Lake area of the Lake Superior region. The parental magma was emplaced as a crystal-liquid suspension, followed by chilling of an evolved, leading edge ferrodioritic liquid in the basal zone of the intrusion. The conditions of emplacement at the present crustal location were ∼1,150°C, 2 kbar, and fO2 slightly above the wustite-magnetite (WM) buffer. The main differentiation process after emplacement was the sorting and redistribution of plagioclase and olivine crystals on a local scale accompanied by less efficient convection and minor settling of olivine. Calculated crystallization sequence for the parental magma is olivine+plagioclase (1,240°C)→olivine+plagioclase+magnetite (1,146°C, WM+0.5)→olivine+plagioclase+magnetite+augite (1,140°C, WM+0.5). The calculated compositions of the cumulus olivine and plagioclase in equilibrium with the parent magma at 1,150°C are Fo66.7±1.1 and An64.5±2.5, respectively, and are similar to the estimated average composition of primary olivine (Fo69.1±2.8) and the average composition of plagioclase core (An66.3±2.8) measured in drill core samples through the intrusion (Chalokwu and Grant 1987).

The Dovyren Intrusive Complex: Problems of Petrology and Ni Sulfide Mineralization

This paper presents a review of petrological-geochemical studies at the Yoko-Dovyren Massif with an emphasis on relations between parameters of the parental magma, a model for the genesis of the lower contact zone, and the nature of Ni sulfide ore mineralization, including the evaluation of the possible ore potential. Arguments are presented in support of the conclusion that the Dovyren magma brought much intratelluric olivine of the composition Fo85–87 into the chamber, and the composition of the initial melt corresponded to gabbronorite or moderately magnesian basite with no more than 10 wt % MgO. The probable temperature of the parental magma was approximately 1200–1250°C, and the sulfur solubility did not exceed 0.10–0.12 wt % (P = 1 kbar, WM buffer). The comparison of this estimate with the average S contents in the bottom plagioperidotites (0.12±0.06 wt %) indicates that the initial magma was saturated with a sulfide phase. For the first time the problem of the composition of contaminated dunites was formulated (these rocks occur in the Layered Series and contain more magnesian olivine Fo87–92). The reason for the increase in the mg# of olivine is thought to be the partial melting and compaction of the original cumulates due to the infiltration of intercumulus melt enriched in volatile components. The volatiles were presumably provided by the thermal decomposition of carbonate xenoliths, a process that resulted in an increase in the CO2 pressure and the transfer of calcite-magnesite components of carbonates into the melt. This follows from (1) the occurrence of magnesian skarn developing after carbonates, (2) high CaO contents in olivine form the contaminated dunite, (3) the appearance of olivine-bearing pyroxenites and wehrlites in the upper part of the dunite zone, (4) correlation between the olivine and chromite composition in the contaminated and uncontaminated dunites, (5) broad variations in the oxygen isotopic composition of olivine and plagioclase from rocks of the Layered Series, (6) experimental data on the dissolution of carbonates in alkali basalt melts, and (7) analogies with isotopic-geochemical characteristics of rocks from the Jinchuan ultramafic complex. Petrological implications of the interpretation of the Dovyren chamber are discussed with reference to closed and flow-through (during an initial stage) magmatic systems. A petrological-geological model is proposed for the genesis of the Synnyr-Dovyren volcanic-plutonic complex and related Ni sulfide ore mineralization. The potential resources of Cu-Ni sulfide ores in the plagioperidotites are evaluated with regard to the still-unexposed part of the massif.

INFOREX-3.0: a database on experimental studies of phase equilibria in igneous rocks and synthetic systems II: data description and petrological applications

The purpose of the paper is to provide details of the application of the INFOREX-3.0 database, a package designed to store, retrieve, and process phase equilibria information. This most recent release of the system accesses data of 162 experimental studies, conducted from 1962 to 1994, including a total of 6174 experiments with 5188 addressed to natural igneous rocks and 986 runs carried out in synthetic systems, mostly CMAS. The total database was divided into 3893 “dry” conditions experiments, and 2281 runs performed in the presence of H2O and/or CO2: 1618 of the “wet” runs represented are water saturated. The number of 1 atm experiments (3750) is greater than the number of high-pressure runs (2474). The INFOREX database contains 8311 coexisting phase compositions: 3197 for glasses, 1247—olivine, 1429—pyroxenes, 501—spinels, 842—plagioclase. One block of the INFOREX information includes 298 liquid compositions where the Fe3+Fe2+ ratio was determined. Data for sulfur and water solubility experiments also have been systematized. The INFOREX data management system allows users to find and print out data on a specific set of mineral-melt or two mineral equilibrium experiments requested for a given range of temperatures, pressures, oxygen fugacities, and compositions in a matter of a few seconds. In addition, one can use subsets of the data to develop mineral-melt geothermometers for equilibria including olivine, plagioclase, pyroxenes, and spinels for any specific system type. Two examples illustrate the use of INFOREX for testing empirical equations proposed for the calculation of water solubility and Fe3+Fe2+ ratio in basic to acid melts.

Test of the Ballhaus–Berry–Green Ol–Opx–Sp oxybarometer and calibration of a new equation for estimating the redox state of melts saturated with olivine and spinel

Testing the Ballhaus–Berry–Green Ol–Opx–Sp oxybarometer (BBG) on independent experimental data indicates that it overestimates the oxygen fugacity by 0.6–1.3 log units under mildly reduced conditions (near the C–CO buffer) and by as much as 2–3 log units under reduced conditions (at the IW buffer and below it). A newly developed oxibarometer is suggested to minimize this effect and enhance the capabilities of redoxometry of low-pressure mineral associations, including magmatic melts undersaturated with respect to orthopyroxene (Opx). The new empirical equation of the oxybarometer is applicable to a wide range of mafic–ultramafic magmas of normal alkalinity, including terrestrial, lunar, and meteoritic systems under pressures of 0.001–25 kbar and oxygen fugacity ranging from IW–3 to NNO + 1. The derived regression fits the ΔQFM values of the calibration dataset (154 experiments) accurate to ~0.5 log units. The new oxybarometer eliminates systematic errors when redox parameters are evaluated for the reduced region (from IW–3 to C–CO) and for crystallization of magmas without Opx on the liquidus. The efficiency of the suggested model is demonstrated by its application to natural rocks: (1) low-Ti lunar basalts, (2) tholeiites from the Shatsky Rise, (3) Siberian flood basalts, (4) rocks of the layered series of the Yoko-Dovyren intrusion, and (5) mantle xenoliths collected in southern Siberia, Mongolia, China, and the southern Russian Far East. The values yielded by such oxybarometers for intrusive rocks, which underwent long-lasting cooling and postcumulus reequilibration, should be regarded with reserve.

SPINMELT-2.0: Simulation of spinel–melt equilibrium in basaltic systems under pressures up to 15 kbar: I. model formulation, calibration, and tests

The paper presents results of testing currently used models proposed to describe Cr-spinel–melt equilibrium: models of the MELTS family by M.S. Ghiorso with colleagues, the SPINMELT program by A.A. Ariskin and G.S. Nikolaev, and the “MELT–CHROMITE spinel calculator” by А.А. Pustovetov and R.L. Roeder. The new calibration of the SPINMELT model presented in this publication enables calculating a sixcomponent (Mg, Fe2+, Cr, Al, Fe3+, and Ti) composition of Cr-spinel and the

SPINMELT-2.0: Simulation of Spinel–Melt Equilibrium in Basaltic Systems under Pressures up to 15 Kbar: II. Description of the Program Package, the Topology of the Cr-spinel–Melt Model System, and Petrological Implications

T - {f_{{o_2}}}