JpGU-AGU Joint Meeting 2017

Presentation information

[EE] Oral

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

[M-GI28] [EE] Data assimilation: A fundamental approach in geosciences

Mon. May 22, 2017 10:45 AM - 12:15 PM 301B (International Conference Hall 3F)

convener:Shin'ya Nakano(The Institute of Statistical Mathematics), Yosuke Fujii(Meteorological Research Institute, Japan Meteorological Agency), SHINICHI MIYAZAKI(Graduate School of Science, Kyoto University), Takemasa Miyoshi(RIKEN Advanced Institute for Computational Science), Chairperson:Shin'ya Nakano(The Institute of Statistical Mathematics)

12:00 PM - 12:15 PM

[MGI28-12] Data assimilation for massive autonomous systems based on a second-order adjoint method

Shin-ichi Ito1, *Hiromichi Nagao1,2, Akinori Yamanaka3, Yuhki Tsukada4, Toshiyuki Koyama4, Masayuki Kano1, Junya Inoue5 (1.Earthquake Research Institute, The University of Tokyo, 2.Graduate School of Information Science and Technology, The University of Tokyo, 3.Graduate School of Mechanical Systems Engineering, Tokyo University of Agriculture and Technology, 4.Graduate School of Engineering, Nagoya University, 5.Research Center for Advanced Science and Technology)

Keywords:data assmilation, adjoint method, phase-field model, uncertainty quantification, Bayesian statistics

We propose an adjoint-based data assimilation 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 diagonal 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. Conventional algorithms for deriving the inverse Hessian matrix require O(CN2+N3) computations and O(N2) memory, where N is the number of degrees of freedom of a given autonomous system and C is the number of computations needed to simulate time series of suitable length. The proposed method using a second-order adjoint method allows us to directly evaluate the diagonal elements of the inverse Hessian matrix without computing all of its elements. This drastically reduces the number of computations to O(C) and the amount of memory to O(N) for each diagonal element. The proposed method is validated through numerical tests using a massive two-dimensional Kobayashi phase-field model. We confirm that the proposed method correctly reproduces the parameter and initial state assumed in advance, and successfully evaluates the uncertainty of the parameter.