An internally consistent dataset with uncertainties and correlations: 3. Applications to geobarometry, worked examples and a computer program

Abstract This paper provides methods and a description of a Pascal computer program, thermocalc, for various thermodynamic calculations using the thermodynamic dataset presented in earlier papers in this series (Holland & Powell, 1985; Powell & Holland, 1985). The dataset involves uncertainties on the thermodynamic parameters and therefore allows uncertainties to be calculated on results, for example in geothermometry and geobarometry. Recommendations are made for the uncertainties on activities to be used in calculations on rocks, particular emphasis being placed on preventing underestimates of these uncertainties at small mole fractions. Apposite examples of phase diagram and rock calculations are presented with ouput from thermocalc, demonstrating the utility of the program. Of the rock calculations, the most valuable are considered to be those involving simultaneous combination ‘least squares’of calculated conditions for a set of reactions applicable to a rock. This set of reactions involves the independent reactions which can be written between the end-members in the minerals in a rock and in the thermodynamic dataset. In contrast to an approach based on specific geothermometers and geobarometers, this approach maximizes the benefit of having an internally consistent thermodynamic dataset. thermocalc is available in IBM PC and Mac versions, from Roger Powell for A$25 or Tim Holland for £10 per version.

[1]  Joseph V. Smith,et al.  Investigations concerning the Breakdown of Albite at Depth in the Earth , 1967, The Journal of Geology.

[2]  D. Kerrick,et al.  Experimental and thermodynamic study of equilibria in the system CaO-MgO-SiO 2 -H 2 O-CO 2 , 1975 .

[3]  E. Ghent.  Plagioclase-garnet-Al 2 SiO 5 -quartz; a potential geobarometer-geothermometer , 1976 .

[4]  J. R. Goldsmith.  The melting and breakdown reactions of anorthite at high pressures and temperatures , 1980 .

[5]  E. Ghent,et al.  Geobarometry and geothermometry of plagioclase-biotite-garnet-muscovite assemblages , 1981 .

[6]  R. C. Newton,et al.  Thermodynamics of the Garnet—Plagioclase—Al2SiO5—Quartz Geobarometer , 1981 .

[7]  R. C. Newton,et al.  Thermodynamic calibration of geobarometers based on the assemblages garnet-plagioclase-orthopyroxene (clinopyroxene)-quartz , 1982 .

[8]  N. Harris,et al.  Geobarometry, Geothermometry, and Late Archean Geotherms from the Granulite Facies Terrain of South India , 1982, The Journal of Geology.

[9]  N. Chatterjee,et al.  The system CaO-Al2O3-SiO2-H2O: new phase equilibria data, some calculated phase relations, and their petrological applications , 1984 .

[10]  J. Ferry,et al.  Constraints on the thermodynamic mixing properties of plagioclase feldspars , 1984 .

[11]  R. Powell,et al.  An internally consistent thermodynamic dataset with uncertainties and correlations: 1. Methods and a worked example , 1985 .

[12]  D. Jenkins,et al.  The high—low albite relations revealed by reversal of degree of order at high pressures , 1985 .

[13]  R. Powell,et al.  An internally consistent thermodynamic dataset with uncertainties and correlations: 2. Data and results , 1985 .

[14]  E. Salje,et al.  Domain wall formation in minerals: I. theory of twin boundary shapes in Na-feldspar , 1985 .

[15]  R. Powell.  Geothermometry and geobarometry: a discussion , 1985, Journal of the Geological Society.