Japan Geoscience Union Meeting 2019

Presentation information

[E] Poster

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

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

Wed. May 29, 2019 5:15 PM - 6:30 PM Poster Hall (International Exhibition Hall8, Makuhari Messe)

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)

[MGI30-P03] Accurate estimation of posterior error covariance in a 4D-Var inverse analysis

*Yosuke Niwa1,2, Yosuke Fujii2,3 (1.National Institute for Environmental Studies, 2.Meteorological Research Institute, Japan Meteorological Agency, 3.The Institute of Statistical Mathematics)

Keywords:inversion, data assimilation, posterior error covariance

The four-dimensional variational method (4D-Var) is one prominent data assimilation/inversion method. However, it is a non-trivial task for a conventional 4D-Var to estimate a posterior error variance-covariance matrix. This study developed a method to estimate a posterior error variance-covariance matrix especially focusing on its off-diagonal elements, i.e., covariance. Off-diagonal elements are usually difficult to be estimated but provide error correlations that are beneficial for interpreting optimized parameter variations. The method was developed within a 4D-Var framework using a quasi-Newton method with the Broyden–Fletcher–Goldfarb–Shanno (BFGS) formula. This study newly introduces numerical techniques that ensure conjugacy among the set of increment vector pairs in the BFGS formula. Through application to an inverse problem of atmospheric CO2, the developed method was evaluated with three diagnostic measures and compared against existing methods. The evaluation revealed that the developed method could provide accurate estimates of the posterior variance-covariance matrix, in terms not only of the diagonal but also of the off-diagonal elements. Although far more expensive than optimal state estimation, the computational efficiency was found reasonable for practical use, especially in conjunction with an ensemble approach.