Evaluation of uncertainties in solid–aqueous–gas chemical equilibrium calculations

•This algorithm predicts uncertainties in equilibrium compositions.•A dataset of uncertainties on thermodynamic parameters is tabulated.•This dataset is consistent with the SUPCRT92 database.•Methods of nonlinear optimization are applied.•This code can be used as a core for calculating equilibrium compositions.

Thermodynamic calculations are traditionally carried out under the assumption of specified input parameters. Errors associated to the results are not often estimated. Here, we propose a novel algorithm that propagates the uncertainty intervals on thermodynamic constants to the uncertainty in chemical equilibrium compositions. The computing uses a dataset of uncertainties on thermodynamic parameters for minerals, solution species and gases consistent with the SUPCRT92 database. Also the algorithm of nonlinear optimization is thoroughly described and realized on a base of the CRONO software. This code can be incorporated into reactive mass transport models as a core for calculating equilibrium compositions. The performance of the algorithm is tested in an experimental system involving Mont Terri's Opalinus Clay interacting with pore water. Its effectiveness is also evaluated against Monte Carlo simulations and Latin Hypercube sampling.