An Inverse Method for Quantifying Petrological Parameters and Uncertainty in Phase Equilibrium Modelling

ABSTRACT Phase equilibrium modelling offers a powerful quantitative framework for understanding petrological processes. Yet, many studies still rely on qualitative comparisons between natural datasets and these forward modelled predictions to constrain model parameters, commonly pressure–temperature ( P – T ) conditions. Compounding this, uncertainties from the observed data or within the modelled predictions are rarely quantified, limiting confidence in the estimated P – T conditions and resulting petrological interpretations. We introduce LinaForma, an inverse modelling workflow that determines best‐fit P – T conditions (or other petrological parameters) and their associated uncertainties for a given rock system by minimizing the misfit between observed data (e.g., mineral compositions or modal proportions) and their forward modelled predictions. Uncertainty is quantified by resampling of the observed data with replacement. Diagnostic metrics identify poorly performing variables and assess the sensitivity of the inversion result to variable uncertainty. Applied to an amphibolite‐facies pelite and metabasite from the Greater Himalayan Sequence (Zanskar Himalaya, NW India), the approach proves effective across contrasting model systems that require different sets of solution models and variables, and it produces P – T estimates consistent with classical thermobarometry. The workflow offers several advantages: compatibility with outputs from any forward modelling software; flexible variable selection; systematic grid‐search inversion in multidimensional space; a robust L1‐norm misfit function resistant to outliers; and sensitivity and uncertainty analysis via bootstrap resampling. Limitations include the increasing computational demands for high‐dimensional grids ( N > 2) and the absence of explicit quantification of uncertainties inherited from the thermodynamic dataset and solution models. Alongside other emerging quantitative methods, LinaForma enables petrologists to make more informed interpretations of complex metamorphic systems and target improvements to thermodynamic datasets.

Authors

Institutions

Publication Details

Journal
Journal of Metamorphic Geology
Published
2025-12-05
DOI
https://doi.org/10.1111/jmg.70016
Primary Topic
Geological and Geochemical Analysis
Type
article
Field-Weighted Citation Impact
0.00

Funders

Controls
|||
ALL TIME
JAN
FEB
MAR
APR
MAY
JUN
JUL
AUG
SEP
OCT
article

An Inverse Method for Quantifying Petrological Parameters and Uncertainty in Phase Equilibrium Modelling

Tobermory Mackay-Champion, Ian P. Cawood
Journal of Metamorphic Geology
Geological and Geochemical Analysis
article

An Inverse Method for Quantifying Petrological Parameters and Uncertainty in Phase Equilibrium Modelling

Tobermory Mackay-Champion, Ian P. Cawood
article en

Abstract

ABSTRACT Phase equilibrium modelling offers a powerful quantitative framework for understanding petrological processes. Yet, many studies still rely on qualitative comparisons between natural datasets and these forward modelled predictions to constrain model parameters, commonly pressure–temperature ( P – T ) conditions. Compounding this, uncertainties from the observed data or within the modelled predictions are rarely quantified, limiting confidence in the estimated P – T conditions and resulting petrological interpretations. We introduce LinaForma, an inverse modelling workflow that determines best‐fit P – T conditions (or other petrological parameters) and their associated uncertainties for a given rock system by minimizing the misfit between observed data (e.g., mineral compositions or modal proportions) and their forward modelled predictions. Uncertainty is quantified by resampling of the observed data with replacement. Diagnostic metrics identify poorly performing variables and assess the sensitivity of the inversion result to variable uncertainty. Applied to an amphibolite‐facies pelite and metabasite from the Greater Himalayan Sequence (Zanskar Himalaya, NW India), the approach proves effective across contrasting model systems that require different sets of solution models and variables, and it produces P – T estimates consistent with classical thermobarometry. The workflow offers several advantages: compatibility with outputs from any forward modelling software; flexible variable selection; systematic grid‐search inversion in multidimensional space; a robust L1‐norm misfit function resistant to outliers; and sensitivity and uncertainty analysis via bootstrap resampling. Limitations include the increasing computational demands for high‐dimensional grids ( N > 2) and the absence of explicit quantification of uncertainties inherited from the thermodynamic dataset and solution models. Alongside other emerging quantitative methods, LinaForma enables petrologists to make more informed interpretations of complex metamorphic systems and target improvements to thermodynamic datasets.

Journal of Metamorphic GeologyVol. 44(1)
Planetary Science Institute (US), University of Oxford (GB), University of Hong Kong (HK)
Chinese Academy of Sciences, Natural Environment Research Council
Openalex Percentile: Top 61%
Geological and Geochemical Analysis
AI Navigator

Ask Laika to Summarize, Analyze, and Connect papers live on the map.

Summarize Papers & Methodologies

Extract key findings, datasets, and comparative methods across publications.

Benchmark Rankings & Visual Analytics

Rank top research institutions, authors, funders, topics, and journals by Field-Weighted Citation Impact (FWCI) and paper volume with instant charts.

Connect Distant Disciplines

Bridge topological clusters on the map to find hidden collaborative intersections.