Japan Geoscience Union Meeting 2019

Presentation information

[J] Poster

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

[M-GI33] Data-driven geosciences

Mon. May 27, 2019 1:45 PM - 3:15 PM Poster Hall (International Exhibition Hall8, Makuhari Messe)

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)

[MGI33-P02] Bayesian statistical evaluation method for detrital zircon geochronology

*Tan Furukawa1, Tatsuki Tsujimori1 (1.Tohoku University)

Keywords:Bayesian statistics, detrital zircon geochronology, MCMC algorithm

Over the past two decades, detrital zircon U-Pb geochronology has become an essential method not only in sedimentary provenance analyses but also in characterization of sedimentary and metamorphic strata. Although the age data has been conventionally evaluated by using relative frequency plots (here after called histograms) and density curves, the lack of statistical verification can cause miss reading of data. In order to overcome the statistical issues raised, this study proposes new statistical evaluation method for detrital zircon geochronology.
The new method is based on Bayesian statistics, adapted a theory proposed by Riihimäki and Vehtari (2014) that provides a flexible modeling for density curve using the Logistic Gaussian Process (LGP)
For introducing non-negative consistency for density curve, the logistic density transformation is performed by p(x)=exp(f(x))/∫exp(f(s))ds; x is age, p(x) is density curve, and f(x) is Gaussian Process; f(x)∼GP(0,κ(x,x')). The covariance matrix κ(x,x') is calculated by Gaussian kernel with hyperparameters ρ and σ. For the actual computation, the logistic density transformation has integral difficulties, then, M discretized softmax function is used for approximation for it.
The likelihood function is L(p|ρ,σ)=yTp. p is M discretized p(x), and y is vector of yi, which is the number of observations that fall within the i'th M discretized age subregion.
We have used programming language R and Stan for the MCMC computation of the statistical model above. When the number of detrital zircon data is several hundred or less, the hyperparameters of Gaussian process ρ and σ do not converge, and Bayesian estimation can not be performed. For this reason, the author estimated most plausible ρ and σ by WAIC. Using zircon data from our own data and the literature, in this presentation, we will demonstrate powerfulness of our new statistical evaluation.