Japan Geoscience Union Meeting 2019

Presentation information

[J] Oral

M (Multidisciplinary and Interdisciplinary) » M-GI General Geosciences, Information Geosciences & Simulations

[M-GI33] Data-driven geosciences

Mon. May 27, 2019 3:30 PM - 5:00 PM A08 (TOKYO BAY MAKUHARI HALL)

convener:Tatsu Kuwatani(Japan Agency for Marine-Earth Science and Technology), Hiromichi Nagao(Earthquake Research Institute, The University of Tokyo), Kenta Ueki(Japan Agency for Marine-Earth Science and Technology), Masayuki Kano(Graduate school of science, Tohoku University), Chairperson:Tatsu Kuwatani(JAMSTEC), Suguru Yabe(JAMSTEC)

4:00 PM - 4:15 PM

[MGI33-07] Specifying rate constants and reaction path on the water-rock chemical reaction based on exchange Monte Carlo method and sparse modeling

*Ryosuke Oyanagi1, Tatsu Kuwatani1, Omori Toshiaki2 (1.Japan Agency for Marine-Earth Science and Technology, 2.Kobe University)

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 present a novel framework to estimate multiple reaction-rate constants from the mineral distribution pattern in a rock by using the exchange Monte Carlo method (Hukushima & Nemoto, 1997; Iba, 2001). We applied the method to the synthetic data, which includes nine unknown parameters to be estimated (Oyanagi et al., 2018). As a result, probabilistic distribution of six unknown parameters among them showed normal distributions and the proposed method could find each true parameter. However, the probabilistic distribution of others showed uniform distributions and the proposed method failed to estimate each true value because these parameters did not contribute to the observed dataset. We suggest that effective reaction paths could be extracted from mineral distribution pattern in a rock, by excluding non-effective rate constant using a sparse Modeling.


Hukushima K, Nemoto K (1996) Exchange Monte Carlo method and application to spin glass simulations. J Phys Soc Japan 65:1604–1608. doi: 10.1143/JPSJ.65.1604

Iba Y (2001) Extended Ensemble Monte Carlo. Int J Mod Phys C 12:623–656. doi: 10.1142/S0129183101001912

Oyanagi R, Okamoto A, Tsuchiya N (2018) Multiple Kinetic Parameterization in a Reactive Transport Model Using the Exchange Monte Carlo Method. Minerals 8:579. doi: 10.3390/min8120579