1 / 8

Basic concepts

Maximum Likelihood Estimation and Simplified Kalman Filter tecniques for real time Data Assimilation. Basic concepts. State vector : atmospheric state of the forecast model True state : best representation of the reality

wyanet
Download Presentation

Basic concepts

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Maximum Likelihood Estimation and Simplified Kalman Filter tecniques for real time Data Assimilation

  2. Basic concepts State vector : atmospheric state of the forecast model True state : best representation of the reality Background (or priori) state : estimate of the true state calculated from the analysis Observations vector y :observations data used for a given analysis Innovation vector Observation matrix H (or transition matrix): linear operator

  3. Modelling of errors Background error’s covariance matrix: Background unbiased error: Observation unbiased error: Observation error’s covariance matrix: Note: observation and backround errors are mutually uncorrelated

  4. The Kalman Filter Equations for the analysis and propagation of the model state for each observation time step Innovation Kalman Gain Too large inverted matrix Equations for the analysis and propagation covariance matrices

  5. On line estimation of covariance parameters Maximum likelihood approach suggested by Dee for Hirlam data assimilation [1] Dee, D.P., 1991: Simplification of the Kalman filter for meteorological data assimilation, Q.J.R. Meteorol. Soc.117,365-384. [2] Dee, D.P., 1995: On-line estimation oferror covariance parameters for atmospheric data assimilation, Mon. Wea. Rev.123,1128-1145. Estimation tecniques imposing the optimatility condition of the filter for the innovation process at first tested for a 28-variable atmospheric model [3] Mehra, R.K., 1970: On the identification of variances and adaptive Kalman filtering. IEEE Trans. Automat. Contr. Vol. AC15,175-184. [4] Todini, E., 1978: Mutually interactive state-parameter (MISP) estimation. In Chao-Lin Chu (ed) Applications of Kalman Filter to Hydrology, Hydraulics, and Water Resources. Stochastic Hydraulics Program Dept. Civil Eng. University of Pittsburgh (Penn). [5] Lorenz, E, 1965: A study of the predicatibility of a 28-variable atmospheric model. Tellus XVII, 321-333

  6. Maximum likelihood method We parametrize the aspected value of the innovation as follows We introduce the pdf which it is assumed to be Gaussian We find the estimation of the parameters maximizing the pdf or minimizing the log function of pdf

  7. Application to HIRVDA Data • Selected observational data for each type of observations • Choose a appropiated sub-domain, because too many observations made the dimension of S too large to invert • Calculate optimal parameter in each sub-domain and used to estimated the background error standard deviations used for HIRLAM forecast.

  8. Estimations of statistics for observations and model errors Assuming the following conditions: • The rank of matrices Q and R are equal • The model error is stationary (Q is time indipendent) • For each time step k Appling the last conditions we can calculate the recursive form for R and Q. For his retrictive conditions at first this method will be tested for a 28-variable meteorological model.

More Related