Japan Geoscience Union Meeting 2015

Presentation information

Oral

Symbol M (Multidisciplinary and Interdisciplinary) » M-TT Technology & Techniques

[M-TT40] New frontier of data analysis in geoscience: Data-driven approach

Thu. May 28, 2015 9:00 AM - 10:45 AM 201A (2F)

Convener:*Tatsu Kuwatani(Graduate School of Environmental Sciences, The University of Tokyo), Takeshi Komai(none), Hideaki Miyamoto(The University Museum, The University of Tokyo), Katsuaki Koike(Laboratory of Environmental Geosphere Engineering, Department of Urban Management, Graduate School of Engineering, Kyoto University), Takane Hori(R&D Center for Earthquake and Tsunami, Japan Agency for Marine-Earth Science and Technology), Hiromichi Nagao(Earthquake Research Institute, The University of Tokyo), Chair:Kenta Ueki(Japan Agency for Marine-Earth Science and Technology, Department of Solid Earth Geochemistry), Tatsu Kuwatani(Graduate School of Environmental Sciences, The University of Tokyo)

9:30 AM - 9:45 AM

[MTT40-09] Bayesian analysis of nonlinear dynamics in heterogeneous reactions -- toward understanding of rock-forming processes --

*Toshiaki OMORI1, Tatsu KUWATANI2, Atsushi OKAMOTO2, Koji HUKUSHIMA3 (1.Kobe University, 2.Tohoku University, 3.The University of Tokyo)

Keywords:data analysis in geoscience, data-driven approach, Bayesian statistics, nonlinear dynamics, heterogeneous reaction, rock-forming processes

Heterogeneous reactions, which occur at interface between two phases, play important roles in rock-forming processes. In heterogeneous reactions, dynamics are intrinsically nonlinear, because surface area governs actual reaction rates and changes temporally. In general, even static surface area, not to mention the temporal change of surface area, is very difficult to be measured. Therefore, it has still been difficult to determine kinetic parameters of heterogeneous reactions even by well-designed laboratory experiments, because we can only use nonlinear and incomplete data sets.

In this study, we develop novel methodologies which estimate kinetic parameters from nonlinear dynamics of heterogeneous reactions using a Bayesian probabilistic approach. With an increasing demand of data-driven analysis, the Bayesian estimation has been widely applied to various fields in the natural sciences including physics, medicine, brain science, and earth sciences (e.g. Watanabe et al., 2009; Naruse et al., 2009; Shouno and Okada, 2010; Omori et al., 2013; Kuwatani et al., 2014; Omori, 2014). In this study, we first formulate probabilistic models of nonlinear dynamics of heterogeneous systems based on state space model. Sequential Monte Carlo Method and EM algorithms are employed for the proposed state-space model in order to simultaneously estimate time course of heterogeneous systems and kinetic parameters underlying nonlinear dynamics. The proposed method provides a very suitable framework for inversion analysis of nonlinear dynamics from incomplete datasets.

In this study, we focus a simple heterogeneous reaction, in which the solid reactant changes to solid product via intermediate product dissolved in fluid phase. It can be regarded as one of the most fundamental and substantial heterogeneous reactions for solid-fluid interactions. Using our proposed method, the reaction rate constants of dissolution and precipitation as well as the temporal changes of solid reactants and products were successfully estimated only by the observable temporal change of the dissolved intermediate product. The proposed method can be potentially applied to actual laboratory experiments of heterogeneous kinetics in various fields of earth and planetary sciences.