Multiple kinetic parameterization in a reactive transport model using the exchange Monte Carlo method

Ryosuke Oyanagi, Atsushi Okamoto, Noriyoshi Tsuchiya

Research output: Contribution to journalArticlepeer-review

5 Citations (Scopus)


Water–rock interaction in surface and subsurface environments occurs in complex multicomponent systems and involves several reactions, including element transfer. Such kinetic information is obtained by fitting a forward model into the temporal evolution of solution chemistry or the spatial pattern recorded in the rock samples, although geochemical and petrological data are essentially sparse and noisy. Therefore, the optimization of kinetic parameters sometimes fails to converge toward the global minimum due to being trapped in a local minimum. In this study, we simultaneously present a novel framework to estimate multiple reaction-rate constants and the diffusivity of aqueous species from the mineral distribution pattern in a rock by using the reactive transport model coupled with the exchange Monte Carlo method. Our approach can estimate both the maximum likelihood and error of each parameter. We applied the method to the synthetic data, which were produced using a model for silica metasomatism and hydration in the olivine–quartz–H2O system. We tested the robustness and accuracy of our method over a wide range of noise intensities. This methodology can be widely applied to kinetic analyses of various kinds of water–rock interactions.

Original languageEnglish
Article number579
Issue number12
Publication statusPublished - 2018 Dec


  • Data-driven approach
  • Exchange Monte Carlo method
  • Geochemical inverse problem
  • Markov chain Monte Carlo
  • Optimization
  • Reaction kinetics
  • Reaction rate
  • Reactive transport modeling


Dive into the research topics of 'Multiple kinetic parameterization in a reactive transport model using the exchange Monte Carlo method'. Together they form a unique fingerprint.

Cite this