Japan Geoscience Union Meeting 2019

Presentation information

[J] Oral

S (Solid Earth Sciences ) » S-TT Technology & Techniques

[S-TT46] Seismic Big Data Analysis Based on the State-of-the-Art of Bayesian Statistics

Mon. May 27, 2019 1:45 PM - 3:15 PM A08 (TOKYO BAY MAKUHARI HALL)

convener:Hiromichi Nagao(Earthquake Research Institute, The University of Tokyo), Aitaro Kato(Earthquake Research Institute, the University of Tokyo), Takuto Maeda(Graduate School of Science and Technology, Hirosaki University), Keisuke Yano(The university of Tokyo), Chairperson:Nana Yoshimitsu(Earthquake Research Institute, The University of Tokyo), Takeru Matsuda(Graduate School of Information Science and Technology, The University of Tokyo)

2:15 PM - 2:30 PM

[STT46-03] Uncertainty quantification based on 4DVar data assimilation for massive simulation models

*Shin-ichi Ito1, Masayuki Kano2, Hiromichi Nagao1 (1.The University of Tokyo, 2.Tohoku University)

Keywords:Bayesian statistics, Data assimilation, Uncertainty quantification

Data assimilation (DA) is a fundamental computational technique that integrates numerical simulation models and observation data on the basis of Bayesian statistics. Originally developed for meteorology, especially weather forecasting, DA is now an accepted technique in various scientific fields. One key issue that remains controversial is the implementation of DA in massive simulation models under the constraints of limited computation time and resources. This work proposes a DA technique based on a four dimensional variational (4DVar) method for massive autonomous models that produces optimum estimates and their uncertainties within reasonable computation time and resource constraints. The uncertainties are given as several elements of an inverse Hessian matrix, which is the covariance matrix of a normal distribution that approximates the target posterior probability density function in the neighborhood of the optimum. The proposed method based on a second-order adjoint method allows us to directly evaluate the elements of the inverse Hessian matrix without computing all of its elements. Some numerical tests assuming a massive simulation model confirm that the proposed method works well and enables us to evaluate the uncertainties of the parameters involved in the assumed model faster than the conventional DA approaches. Such information regarding uncertainty is valuable, as it can be used to optimize the design of experiments.