Computation of phase equilibria by linear programming: A tool for geodynamic modeling and its application to subduction zone decarbonation

An algorithm for the construction of phase diagram sections is formulated that is well suited for geodynamic problems in which it is necessary to assess the influence of phase transitions on rock properties or the evolution and migration of fluids. The basis of the algorithm is the representation of the continuous compositional variations of solution phases by series of discrete compositions. As a consequence of this approximation the classical non-linear free energy minimization problem is trivially solved by linear programming. Phase relations are then mapped as a function of the variables of interest using bisection to locate phase boundaries. Treatment of isentropic and isothermal phase relations involving felsic and mafic silicate melts by this method is illustrated. To demonstrate the tractability of more complex problems involving mass transfer, a model for infiltration driven-decarbonation in subduction zones is evaluated. As concluded from earlier closed system models, the open-system model indicates that carbonates are likely to persist in the subducted oceanic crust beyond sub-arc depths even if the upper section of the oceanic mantle is extensively hydrated. However, in contrast to more simplistic models of slab devolatilization, the open-system model suggests slab fluid production is heterogeneous and ephemeral. Computed seismic velocity profiles, together with thermodynamic constraints, imply that for typical geothermal conditions serpentinization of the subducted mantle is unlikely to extend to N25 km depth and that the average water-content of the serpentinized mantle is b2 wt.%. D 2005 Elsevier B.V. All rights reserved.

[1]  David A. Yuen,et al.  The application of multidimensional wavelets to unveiling multi-phase diagrams and in situ physical properties of rocks , 2004 .

[2]  G. Eriksson,et al.  Theoretical Computation of Mineral Assemblages in Pyrolite and Lherzolite , 1983 .

[3]  S. Poli,et al.  Experimentally based water budgets for dehydrating slabs and consequences for arc magma generation , 1998 .

[4]  R. Carlson.  Bound water content of the lower oceanic crust estimated from modal analyses and seismic velocities of oceanic diabase and gabbro , 2003 .

[5]  Thomas H. Brown,et al.  The computation of chemical equilibrium in complex systems containing non-ideal solutions , 1987 .

[6]  J. Connolly,et al.  Petrogenetic grids for metacarbonate rocks: pressure-temperature phase-diagram projection for mixed-volatile systems , 1991 .

[7]  S. Peacock,et al.  High‐resolution models of subduction zones: Implications for mineral dehydration reactions and the transport of water into the deep mantle , 2002 .

[8]  Brian M. Smith,et al.  Cretaceous ocean crust at DSDP Sites 417 and 418: Carbon uptake from weathering versus loss by magmatic outgassing , 1989 .

[9]  S. Stein,et al.  Implications of slab mineralogy for subduction dynamics , 2001 .

[10]  D. Gubbins,et al.  Seismological constraints on the gabbro-eclogite transition in subducted oceanic crust , 1994 .

[11]  G. Bebout.  The impact of subduction-zone metamorphism on mantle-ocean chemical cycling , 1995 .

[12]  J. Connolly,et al.  Metamorphic devolatilization of subducted oceanic metabasalts: implications for seismicity, arc magmatism and volatile recycling , 2001 .

[13]  R. Powell,et al.  Calculating phase diagrams involving solid solutions via non‐linear equations, with examples using THERMOCALC , 1998 .

[14]  D. Gutzler,et al.  Post-125 Ma carbon storage associated with continent-continent collision , 1993 .

[15]  O. J. Kleppa,et al.  Thermochemistry of the high structural state plagioclases , 1980 .

[16]  K. Caldeira,et al.  Paleoatmospheric consequences of CO2 released during early Cenozoic regional metamorphism in the Tethyan orogen , 1993 .

[17]  William R. Smith,et al.  The Computation of Chemical Equilibria in Complex Systems , 1980 .

[18]  L. Stixrude,et al.  Mineralogy and elasticity of the oceanic upper mantle: Origin of the low‐velocity zone , 2005 .

[19]  R. Powell,et al.  Mixing properties and activity-composition relationships of chlorites in the system MgO-FeO-Al2O3-SiO2-H2O , 2006 .

[20]  J. Connolly.  Multivariable phase diagrams; an algorithm based on generalized thermodynamics , 1990 .

[21]  Matthias Hort,et al.  Serpentine and the subduction zone water cycle , 2004 .

[22]  Mark S. Ghiorso,et al.  Chemical mass transfer in magmatic processes IV. A revised and internally consistent thermodynamic model for the interpolation and extrapolation of liquid-solid equilibria in magmatic systems at elevated temperatures and pressures , 1995 .

[23]  Robert A. Berner,et al.  Modeling the geochemical carbon cycle , 1989 .

[24]  R. Powell,et al.  Relating formulations of the thermodynamics of mineral solid solutions: Activity modeling of pyroxenes, amphiboles, and micas , 1999 .

[25]  Mark S. Ghiorso,et al.  Chemical mass transfer in magmatic processes , 1987 .

[26]  戸高 法文,et al.  Geochemistry , 2019, Nature.

[27]  R. Powell,et al.  Thermodynamics of order-disorder in minerals: II. Symmetric formalism applied to solid solutions , 1996 .

[28]  J. Connolly,et al.  A coupled petrological–tectonic model for sedimentary basin evolution: the influence of metamorphic reactions on basin subsidence , 2001 .

[29]  G. Abers.  Seismic low-velocity layer at the top of subducting slabs: observations, predictions, and systematics , 2005 .

[30]  S. Hart,et al.  In search of a bulk-Earth composition , 1986 .

[31]  James B. Thompson,et al.  Entropy of mixing in sanidine , 1979 .

[32]  James P. Ignizio,et al.  Linear Programming in Single- and Multiple-Objective Systems , 1984 .

[33]  C. Bina,et al.  Free Energy Minimization by Simulated Annealing with Applications to Lithospheric Slabs and Mantle Plumes , 1998 .

[34]  Roger Powell,et al.  An internally consistent thermodynamic data set for phases of petrological interest , 1998 .

[35]  R. Powell,et al.  Phase relations in high-pressure metapelites in the system KFMASH (K2O–FeO–MgO–Al2O3–SiO2–H2O) with application to natural rocks , 2003 .

[36]  P. Asimow,et al.  Calculation of Peridotite Partial Melting from Thermodynamic Models of Minerals and Melts, IV. Adiabatic Decompression and the Composition and Mean Properties of Mid-ocean Ridge Basalts , 2001 .

[37]  J. Connolly,et al.  An algorithm and computer program for calculating composition phase diagrams , 1987 .

[38]  C. Sotin,et al.  Computation of seismic profiles from mineral physics: the importance of the non-olivine components for explaining the 660 km depth discontinuity , 1998 .

[39]  D. Kerrick.  Serpentinite Seduction , 2002, Science.

[40]  James A. D. Connolly,et al.  An automated strategy for calculation of phase diagram sections and retrieval of rock properties as a function of physical conditions , 2002 .

[41]  S. H. Storey,et al.  The computation of chemical equilibria , 1970 .

[42]  B. Wood,et al.  A thermodynamic model for subsolidus equilibria in the system CaOMgOAl2O3SiO2 , 1984 .

[43]  F. Dahlen.  Metamorphism of nonhydrostatically stressed rocks , 1992 .

[44]  M. Gottschalk.  Internally consistent thermodynamic data for rock-forming minerals in the system SiO 2 -TiO 2 -Al 2 O 3 -CaO-MgO-FeO-K 2 O-Na 2 O-H 2 O-CO 2 , 1996 .

[45]  M. Fuhrman,et al.  Ternary-feldspar modeling and thermometry , 1988 .

[46]  Charles H. Langmuir,et al.  The chemical composition of subducting sediment and its consequences for the crust and mantle , 1998 .

[47]  Simon M. Peacock,et al.  Subduction factory 2. Are intermediate‐depth earthquakes in subducting slabs linked to metamorphic dehydration reactions? , 2003 .

[48]  L. Rüpke,et al.  Are the regional variations in Central American arc lavas due to differing basaltic versus peridotitic slab sources of fluids , 2002 .

[49]  S. Peacock.  Large-scale hydration of the lithosphere above subducting slabs , 1993 .

[50]  E. Oelkers,et al.  SUPCRT92: a software package for calculating the standard molal thermodynamic properties of minerals, gases, aqueous species, and reactions from 1 to 5000 bar and 0 to 1000 ° C , 1992 .

[51]  Selmer M. Johnson,et al.  Chemical Equilibrium in Complex Mixtures , 1958 .

[52]  Stephan V. Sobolev,et al.  Modeling of mineralogical composition, density and elastic wave velocities in anhydrous magmatic rocks , 1994 .

[53]  T. Gerya,et al.  Inherent gravitational instability of thickened continental crust with regionally developed low- to medium-pressure granulite facies metamorphism , 2001 .

[54]  D. Kerrick.  Present and past nonanthropogenic CO2 degassing from the solid earth , 2001 .

[55]  S. M. Junaid Zaidi,et al.  Linear Programming in Single and Multiple Objective Systems , 1982 .

[56]  M. Ghiorso.  Chemical mass transfer in magmatic processes , 1985 .

[57]  M. Hirschmann,et al.  Deep global cycling of carbon constrained by the solidus of anhydrous, carbonated eclogite under upper mantle conditions , 2004 .

[58]  D. Shaw.  GEOCHEMISTRY OF PELITIC ROCKS. PART III: MAJOR ELEMENTS AND GENERAL GEOCHEMISTRY , 1956 .

[59]  J. Connolly,et al.  Subduction of ophicarbonates and recycling of CO2 and H2O , 1998 .

[60]  Konstantin V. Chudnenko,et al.  Modeling chemical mass transfer in geochemical processes; thermodynamic relations, conditions of equilibria and numerical algorithms , 1997 .

[61]  Mark S. Ghiorso,et al.  The pMELTS: A revision of MELTS for improved calculation of phase relations and major element partitioning related to partial melting of the mantle to 3 GPa , 2002 .

[62]  D. Miller,et al.  Mantle wedge water contents estimated from seismic velocities in partially serpentinized peridotites , 2003 .

[63]  J. Connolly,et al.  Metamorphic devolatilization of subducted marine sediments and the transport of volatiles into the Earth's mantle , 2001, Nature.

[64]  J. Connolly,et al.  Metamorphic controls on seismic velocity of subducted oceanic crust at 100–250 km depth , 2002 .

[65]  R. Berman,et al.  Optimized standard state and solution properties of minerals , 1996 .

[66]  P. Kelemen,et al.  Relationship between seismic P‐wave velocity and the composition of anhydrous igneous and meta‐igneous rocks , 2003 .

[67]  R. Powell,et al.  Mixing properties and activity-composition and relationships of chlorites in the system MgO-FeO-Al 2 O 3 -SiO 2 -H 2 O , 1998 .