Premium
Improvement of retrospective optimal interpolation by incorporating eigen‐decomposition and covariance inflation
Author(s) -
Song HyoJong,
Lim GyuHo
Publication year - 2011
Publication title -
quarterly journal of the royal meteorological society
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.744
H-Index - 143
eISSN - 1477-870X
pISSN - 0035-9009
DOI - 10.1002/qj.911
Subject(s) - covariance , weighting , algorithm , computer science , interpolation (computer graphics) , mathematics , rank (graph theory) , divergence (linguistics) , statistics , mathematical optimization , artificial intelligence , medicine , motion (physics) , linguistics , philosophy , combinatorics , radiology
Retrospective optimal interpolation (ROI) is a method that is used to minimize cost functions with multiple minima without using adjoint models. We address two weaknesses associated with the cost‐effective formulation of ROI and offer possible solutions. The first weakness of the cost‐effective ROI formulation is that the error tolerance would become large in practical application due to computation costs. When the error tolerance is large, accuracy‐saturated modes do not extract information from new incoming observations even though they are likely to be flawed. To address this problem, we modify the existing ‘reduced‐resolution’ formulation by using eigen‐decomposition of the background error covariance at each analysis step. We refer to the modified algorithm as a ‘reduced‐rank’ algorithm. This modification allows us to deal with larger error variances while analysing the same number of control variables, because eigen‐decomposition steeply reorders the error variance in descending order. As a result, when the reduced‐rank algorithm is applied, the number of analysed control variables becomes smaller than when the reduced‐resolution algorithm is used with the same error tolerance. The second weakness is the underestimation of the trailing mode of background error covariance that is projected onto the future observation space. This originates from errors in control variables extracted from analysis procedures. To prevent the occurrence of filter divergence due to this underestimation, we reduce the weighting of the observational increments in the analysis. By implicitly assuming that the rate of being projected onto the future observation space of the trailing eigenmodes is similar to that of the leading eigenmodes, we develop a method that inflates the observation error covariance and consequently improves analysis quality in the Lorenz 40‐variable and the 960‐variable model experiments. Copyright © 2011 Royal Meteorological Society