Hello! I’m glad these tools continue to be a useful resource for so many. I do regularly update this page when I find or develop new tools, but I don’t often check the comments posted here. If you have a question about any of these tools, feel free to reach out via email at kayla [dot] iacovino [at] nasa [dot] gov
TAS Diagram Plotter
Excel spreadsheet that allows for fast plotting onto a TAS (total alkali silica) diagram for use in naming your volcanic rock by composition. Created using coordinates as given in Le Maitre et al (2002).
Please cite as: Iacovino, Kayla, & Gouard, Cédric. (2021). TAS Diagram Plotter (3.0). Zenodo. https://doi.org/10.5281/zenodo.5907859
TAS Diagram Plotter version en espanol
Versión en español del trazador de diagramas TAS traducido del inglés por Crescencio Fernandez H.
Please cite as: Iacovino, Kayla, Gouard, Cédric, & Fernandez, Crescencio. (2021). TAS Diagram Plotter version en espanol. Zenodo. https://doi.org/10.5281/zenodo.5907888
Version 3.0 New version by Cédric Gouard. Significant improvements and new features! Logic based on rock composition will now classify your rock type even further than the fields on the TAS diagram. For example, a trachybasalt may be further identified as a hawaiite or potassic trachybasalt; a basaltic trachyandesite may be further identified as a mugearite or shoshonite. The TAS diagram is now drawn in a more sensible way from formulae, and overall user friendliness is improved.
Version 2.0 The user can now resize the TAS diagram chart and keep the region boundaries in tact. Beware that the region labels may still move around relative to their fields, however.
Plutonic TAS Diagram Plotter
Excel spreadsheet based on the Volcanic TAS Diagram Plotter (v2.0) with updated parameters and labels for plutonic rocks. Created by Paul Asimow using coordinates and labels as given in E. A. K. Middlemost (1994).
Please cite as: Iacovino, Kayla, & Asimow, Paul. (2018). Plutonic TAS Diagram Plotter (1.0). Zenodo. https://doi.org/10.5281/zenodo.5907892
Version 1.0 Based on Volcanic TAS Diagram Plotter v2.0
DensityX: Glass Density Calc
Excel spreadsheet. Input the composition of a silicate glass/melt in terms of wt% and the H2O content in terms of wt% plus the temperature and pressure of the sample. Spreadsheet outputs the density of your glass/melt in terms of g/cm3 and g/L. User inputs go in blue boxes, output is automatically generated in red boxes.
References: Iacovino & Till (2019); Lange & Carmichael (1990); Ochs & Lange (1999); Also see Hack & Thompson (2011)
Please cite as: Iacovino, K. and Till, C. B. (2019) “DensityX: A program for calculating the densities of magmatic liquids up to 1,627 °C and 30 kbar”, Volcanica, 2(1), pp. 1–10. doi: 10.30909/vol.02.01.0110.
Excel Version 4.1
4.1: Updated molecular weight of H2O from 18 to 18.02 to match python code and web app. Updated 04 Apr 2023.
4.0: Updated to match DensityX python code version 1.1.0, which includes all changes made to the program in response to reviewers of the DensityX manuscript for publication in Volcanica.
3.2: Updated reference list to include all values and reference P and T values used. Updated 10 Apr 2017.
3.1: Minor tweaks to layout to clarify user inputs and program outputs. Fixed typo in reference temperature for H2O. Updated 29 Mar 2017.
3.0: Updated temperature dependence data to refer to correct reference temperature. Calculations now work at high temperature. Updated 26 Jan 2014.
🐍 DensityX Python Library
You can download the DensityX python library locally if you like. It can also be installed via pip “pip install densityx”.
🕸️ DensityX Web Tool
Online interface to the DensityX python library. Simply upload your file and run it!
Calculate fO2 Buffer
Excel spreadsheet. Tool to calculate oxygen fugacity in terms of the common buffers when logfO2 is known or for translating between fO2 values expressed in terms of various buffers. This is useful for converting calculated oxygen fugacity values into and between delta buffer values (e.g. ΔNNO, ΔQFM…).
References:
Please cite as: Iacovino, Kayla. (2022). Calculate fO2 Buffer (1.5). Zenodo. https://doi.org/10.5281/zenodo.7387196
Version 1.5 1 Dec 2022 Fixed Re-ReO buffer (previously used log when ln should have been used). Credit Daniel Coulthard for the fix!
Version 1.4 19 July 2021 Added equation for Mn-MnO buffer. Also flagged which reactions are calibrated only for 1 bar ( these equations have no pressure dependence).
Version 1.3 3 June 2021 Added equations for Graphite-CO-CO2 (GCO), Re-ReO, Cr-Cr2O3, and MoMoO2 buffers. Added “Translate from IW” tab.
Version 1.2 28 May 2020 Added equation for Si-SiO2 buffer.
Version 1.1 Added ability to convert from a buffer value to the actual logfO2 value and into other buffers. Useful, for example, in comparing two data sets where one is expressed in terms of ΔNNO and the other in terms of ΔQFM. Updated 6 October 2016.
🐍 Python Script: Oxygen Fugacity (Redox) Buffer Plotter
Python script that plots fO2 buffer values at any given pressure and temperature range. Right click the link above and “Save Target As…” and unzip the .zip file containing a python script and library. User should edit the file ‘create_plot.py’ and run this script to create a plot.
Please cite as: Iacovino, Kayla. (2021). Oxygen Fugacity (Redox) Buffer Plotter. Zenodo. https://doi.org/10.5281/zenodo.5907875
Update 6 July 2021 Fixed QFM, HM, CoCoO, and QIF equations to use pressure in bars instead of GPa.
Ferric/Ferrous, Fe3+/FeT, fO2 Converter (Kress and Carmichael, 1991)
Excel spreadsheet. Tool for calculating the ferric/ferrous iron ratio in a melt for any given fO2 and vice versa. Also calculates new wt% and mole fraction values for FeO and Fe2O3 in your melt. Useful for determining theoretical iron speciation when total is given as FeOtotal or Fe2O3total and fO2 is known. Can also calculate or recalculate FeO and Fe2O3 wt% from given Fe3+/FeT and FeOT (or existing FeO and Fe2O3 wt%).
References: Kress & Carmichael (1991); B. R. Frost in Mineralogical Society of America “Reviews in Mineralogy” Volume 25.
Please cite as: Iacovino, Kayla. (2021). Ferric/Ferrous, Fe3+/FeT, fO2 Converter (Kress and Carmichael, 1991) (3.2). Zenodo. https://doi.org/10.5281/zenodo.5907844
Version 3.2
3.2: Corrected error in cells J5 and J6, which were previously using an improper conversion coefficient (0.8998 instead of 2). January 2021.
3.1: Added tab for calculating FeO and Fe2O3 from Fe3+/FeT. December 2020.
3.0: Added all parameters for calculations at high pressure. Redesigned user interface. Added conversion from melt comp to fO2. 12 October 2018.
2.2: Fixed bug where dropdown box to choose oxygen fugacity buffer did not work. 11 October 2016.
2.1: Fixed a critical error in the calculation of XFeO based on the ratio given by the equation in Kress & Carmichael. 15 October 2015.
WetMelts Conversion Tool
Excel spreadsheet with VisualBasic macros. Tool for converting anhydrous composition of a silicate melt in wt% plus known mol fraction of volatiles in the melt (H2O and CO2 mixture) to a hydrous melt composition in terms of both wt% and mol fraction. I created this tool as part of an independent study during my undergrad. Also see the tutorial.
WtpertoMolper Conversion Tool
Excel spreadsheet with VisualBasic macros. Tool for converting (hydrous or anhydrous) melt composition from wt% to mol%.
Oxide to Element Conversion Tool
Excel spreadsheet that will convert compositions measured in oxides (SiO2, FeO, K2O, etc) to elements (Si, Fe, K, etc), and vice versa.
Please cite as: Iacovino, Kayla. (2019). Oxide to Element Conversion Tool. Zenodo. https://doi.org/10.5281/zenodo.5907879
Water From Iron Reduction
Excel spreadsheet. Great for experimental petrologists. If a silicate melt is “fused” or “pre-melted” at the same oxygen fugacity that the experiments are ultimately run at, the ratio of ferric to ferrous iron in the melt will remain unchanged. If, however, the melt is fused in air and then run at more reducing conditions, some of the Fe2O3 will be reduced to FeO and will add water to your melt in the process. Use this tool to calculate just how much water was added to your melt due to the reduction of iron during your run. Only works if pre-run fO2 is in air (i.e. all iron is Fe2O3 before run), but calculation can be made for any experimental fO2.
References: Kress & Carmichael (1991); B. R. Frost in Mineralogical Society of America “Reviews in Mineralogy” Volume 25.
Please cite as: Iacovino, Kayla. (2020). Water From Iron Reduction (1.1). Zenodo. https://doi.org/10.5281/zenodo.5907885
Version 1.1
🐍 Link: Python Libraries for Earth Scientists
A curated list of python packages for use in the Earth Sciences by Maurizio Petrelli.
Link: Gabbrosoft Mineralogical Spreadsheets
Spreadsheets for calculating mineral formulas based on an electron microprobe analysis. The spreadsheet does all the hard work, so you don’t have to! Also find the Bulk KD calculator for calculating trace element bulk KD values given the mineral make-up of a sample, and the Magma Density calculator for calculating the density of a magma from its major element composition.
Link: Mineral Formula Calculations
Various excel spreadsheets from Carlton College. Useful tools for calculating formulas of various mineral types based on electron microprobe or EDX measurements (or anything where you have an analysis of a mineral in terms of wt% oxides).
Link: rockPTX
An excellent resource for optical petrology and mineralogy and some electron microscopy. This website is a resource put together by Frank Mazdab and contains numerous high resolution scans of thin sections (PPL and XPL) along with explanations of the mineralogy. Also lots of good tips for doing optical mineralogy.
Link: Chemical Equation Balancer
A fantastic tool that balances chemical equations (calculates reaction stoichiometry) and then allows for the calculation of moles or mass of product(s) given moles or mass of the reactant(s).
Some of my favorite models, none of which are by me but all of which are used by me!
Papale model for H2O-CO2 partitioning/solubility in magmas
Calculate isobars, fluid/melt partitioning, saturation pressures, solubilities! Calculations straight from the model of Papale P, Moretti R, Barbato D (2006) The compositional dependence of the saturation surface of H2O + CO2 fluids in silicate melts. Chemical Geology 229, 78-95.
MagmaSat for Mac
MagmaSat is a calculator that determines the saturation conditions for a mixed H2O-CO2 fluid in natural composition silicate liquids. The calculator implements the thermodynamic model of Ghiorso and Gualda (2015).
VolatileCalc
This software calculates H2O-CO2-melt equilibrium for the rhyolite and basalt systems. The software is written as Visual Basic for Applications Macros within Microsoft Excel. Calculations from the model of Newman, S., Lowenstern, J.B., 2002, VolatileCalc: a silicate melt-H2O-CO2 solution model written in Visual Basic for Excel: Computers and Geosciences vol. 28 no. 5, pp. 597-604.
SolEx
This model calculates basaltic melt and coexisting vapor compositions in the system C–O–H–S–Cl. In the forward model, total volatile inventories and melt composition are specified by the user. The parameterisation of Dixon (1997, p. 368) is used to predict the partitioning of CO2 and H2O between vapour and melt phases. An iterative procedure is employed to predict the partitioning of S and Cl components between fluid and melt phases. Melt and gas compositions and gas volume fraction are thereby modelled over pressures in the range 5–4000 bar. This approach satisfactorily reproduces independent literature data on S and Cl behaviour in basalt. SolEx is a user-friendly software package available for OS X and Windows, facilitating modelling of closed- and open-system C–O–H–S–Cl degassing in basalts. Computers and Geosciences vol. 45, pp. 87–97.
MELTS, rhyoloiteMELTS, and pMELTS
MELTS is a software package designed to facilitate thermodynamic modeling of phase equilibria in magmatic systems. It provides the ability to compute equilibrium phase relations for igneous systems. MELTS works over the temperature range 500-2000 °C and the pressure range 0-2 GPa. rhyoliteMELTS is identical to the MELTS calibration but is adjusted specifically for silicic compositions containing quartz and feldspar (P/T range equal to that of MELTS). pMELTS is optimized for bulk mantle compositions and is suitable over the temperature range 1000-2500°C and the pressure interval 1-3 GPa.
MELTS for Excel
MELTS for Excel is a macro-enabled Microsoft Excel workbook that provides a graphical user interface to set up, run, retrieve, and display calculations made with MELTS. Original publication: Gualda G.A.R., Ghiorso M.S. (2015) MELTS_Excel: A Microsoft Excel-based MELTS interface for research and teaching of magma properties and evolution. Geochemistry, Geophysics, Geosystems 16(1), 315-324 doi:10.1002/2014GC005545
Model for Viscosity of Volatile-bearing Melts
From the paper Viscosity of magmatic liquids: A model by Giordano D, Russell JK & Dingwell DB (2008). Calculate the viscosity of a melt given composition, temperature, and H2O. The link above is a version of the model in javascript. There is also an Excel spreadsheet version.
QUILF
QUILF is a program originally written in Pascal to assess equilibria among Fe-Mg-Mn-Ti oxides, pyroxenes, olivine, and quartz. The program can calculate mineral equilibria, P, T, and fO2 given the right parameters. Only works on Windows. This version of QUILF was made to run on Windows 3.1 and 95, but I have successfully run it on Windows XP. Be sure to read the QUILF Basics file, included in the zipped archive or as a standalone download. Original publication: Anderson, DJ, Lindsley, DH, and Davidson, PM (1993) A Pascal program to assess equilibria among Fe-Mg-Mn-Ti oxides, pyroxenes, olivine, and quartz. Computers and Geosciences vol. 19, no. 9, pp. 1333-1350. doi:10.1016/0098-3004(93)90033-2
I always find it difficult to find the thermodynamic data that I need when I need it. So, I’m starting a running list here. Below are some references and online tools and tables that I consult often when I need to get my fugacity on.
Fugacity Coefficient Tables
Some very handy tables that are hosted as course materials for a geoscience class at Penn State.
Patrick Barrie’s program for solving cubic equations of state
A man to whom I owe a great debt, Patrick Barrie, has made this online applet for calculating Z, V, φ, H & S using a variety of equations of state, including the Redlich-Kwong. You must know critical temperature (Tcrit), pressure (Pcrit), and the acentric factor ω.
Critical Constants and Acentric Factors for Selected Fluids
A page out of the book Molecular Thermodynamics of Fluid-Phase Equilibria, Third Edition by John M. Prausnitz – University of California, Berkeley; Rüdiger N. Lichtenthaler – University of Heidelberg; Edmundo Gomes de Azevedo – Institute Superior Técnico, Lisbon
Critical constants and accentric factors for some inorganic compounds
A page out of Perry’s Chemical Engineers’ Handbook, 1999.
Selected Properties of Molecules
This is an image of Table 2 from Holloway (1977) Fugacity and Activity of Molecular Species in Supercritical Fluids. In Fraser, D.G., Thermodynamics in Geology, 161-181.
Thermodynamic Quantities for Substances and Ions at 25°C
Link to a table of values for ΔH, ΔG, and ΔS for various substances.
Heats, free energies, and equilibrium constants of some reactions involving O2, H2, H2O, C, CO, CO2, and CH4
Wagman et al. (1945). Paper with calculations and thermodynamic data for reactions involving the above species. This is useful when an equilibrium constant of reaction (rather than equilibrium constant of formation from the elements) is required. For example, when dealing with the reaction CO + 1/2O2 = CO2. Jump to Table 8 for equilibrium constant values.
Joaquin Cortes
October 2, 2015
Hi,
A couple of years ago I wrote Pingu, which is an online tool to generate petrologic diagrams in the context of the Vhub cyberinfrastructure and I would love if you can include it as a resource in your page. In order to access the tool you need to sign for Vhub (https://vhub.org) and join the Pingu group at https://vhub.org/groups/pingu
Then the tool can be launched at https://vhub.org/tools/pingu
Thanks
Joaquin
Maryam
December 19, 2015
Hi
I’m Maryam Khosravi, a PhD student of economic geology from Department of Earth Sciences at Shiraz University. I am working on porphyry copper deposit in my PhD thesis. Could I ask you kindly to reply my question or to introduce papers and references? About this topic can you help me, if possible?
My question is that what analysis can I use for the exact determination of oxygen fugacity (fO2, quantity not quality) in alteration zones such as potassic and sercite alterations, and granitoid intrusions on porphyry copper deposits?
I’d pleasure to help me.
Sincerely,
Maryam
Grace
November 14, 2016
Hello, where do I find a reference to the TAS spreadsheet. I’m using it in an assignment and would like to reference the creator. Thank you
Kayla Iacovino
November 14, 2016
Hi Grace. The coordinates that I used to plot the TAS diagram in excel are given in Le Maitre et al (2002), which is linked above. All I did was plot those up in excel to use as a template for plotting data. You’re welcome to cite this webpage with my name (Kayla Iacovino) as the author, but it’d be good to cite Le Maitre et al (2002), too.
Franco
July 19, 2017
Hello, thanks for your TAS spreadsheet!
But if I insert more than 20 data in “YOUR DATA” table, it does not give the plot.
Ho to modify this table?
Thanks, Franco
Kayla Iacovino
August 29, 2017
Hi Franco! Great question. First, put in all the data you want to plot. Then, right click on the graph and click on “Select data…”. Now you should see a box that lists all of the Legend entries. Click on “Your Data”, and adjust the box to the right of “X values” and “Y values” to include all rows that have your data in it. For now, to help with this a bit, I’ve re-uploaded the spreadsheet with a new version that will plot data down to row 100.
Tashane Boothe
October 7, 2017
Hi,
How do i identify the samples on the graph? is there a way to change the icons or label the points on the graph?
Thanks
Jimmy
January 16, 2018
If you type a Phonolite that may reach up to 15% alkalis, the scale on the y axis change and the diagram is mess up, is there a way to fix this?
Jimmy
January 16, 2018
Is me again, the problem seems to be solved by changing the maximun value for the x axis from 78 to 80 and for the y axis from 16 to 18, some high-silica rhyolites and phonolites either change the location of the names or the point doesn’t appear, giving the impresion that the diagram is not plotting data or that the field areas change, that happened to me when i computed a data set with both high-silica rhyolites and phonolites, but the coordinates of the intersections are fixed : ). That kind of confusion is solved increasing a little bit the maximun values of the axis, the diagram end a little bit oversized, but can be easily edited.
Pablo
November 15, 2018
Hi Dr Kayla, thanks for your site!
I have tried to simulate published results of delta QFM and NNO (e.g. Moussllam et al 2016) providing the melt composition and Fe+3/FeT ratios using the tool ferric/ferrous… but I got comparatively overestimated results. Any possible error with the data I entered?
Thank you,
Dawn
November 30, 2018
Hey Kayla,
Love the TAS plotter. Could you also add in the tholeiite, alkaline magma line?
Cheers!
Sara
June 13, 2019
Hi
I am working on an olivine gabbro intrusion and need Dymek spreadsheet. Can you please help me, How i can download it? I really appreciate any help. Thank you
Abhishek Kumar Madhesiya
March 24, 2020
Hello,
I am a young research scholar in geoscience working on silicic rocks. I am looking for tectonic discrimination diagram for this system. Can you help me to find this?
Thanking you,
Best Wish.
Andy Townley
July 18, 2020
Can anyone help me with this. My research is looking at the occurrence of mafic (or so it seems at first sight) lavas on a predominantly silicic OIV. During lock down, I have been crunching EMPA and XRF data from my samples like its going out of fashion, but one thing I can find absolutely no help for is Titanite and foids (of various hues) – relatively common phases in my samples. Can anyone help me in terms of recommending recalculation spreadsheets, recommended plots, etc? Anything to enable me to make sense of these weird (but wonderful minerals and where they fit into my big scheme of things (which, to be fair is a little vague right now). Help! Thanks in advance. Best wishes
Andy Townley
MOAGINYANA LEKAUNYANE
March 13, 2021
Good Day
I am looking for the spreadsheet or a tool which can help me calculate Oxygen Fugacity Conditions using mono-mineralic compositions such as Perovskite and Monticellite Oxybarometry.