assimilation of multiple linearly dependent data vectors
play

Assimilation of Multiple Linearly Dependent Data Vectors Trond - PowerPoint PPT Presentation

Assimilation of Multiple Linearly Dependent Data Vectors Trond Mannseth NORCE Energy Linearly dependent data vectors Assume that we want to assimilate the data vectors { d l } L l =1 , where { d l = B l d L } L 1 l =1 and { B l } L 1 l


  1. Assimilation of Multiple Linearly Dependent Data Vectors Trond Mannseth NORCE Energy

  2. Linearly dependent data vectors Assume that we want to assimilate the data vectors { d l } L l =1 , where { d l = B l d L } L − 1 l =1 and { B l } L − 1 l =1 denotes a sequence of matrices

  3. Linearly dependent data vectors Main issue Assume that we want to assimilate the data vectors { d l } L l =1 , where { d l = B l d L } L − 1 l =1 and { B l } L − 1 l =1 denotes a sequence of matrices What is the appropriate way to assimilate such a data sequence, taking into account that some, but not necessarily all, information is used multiple times?

  4. Outline Motivation for considering linearly dependent data vectors Relation to multiple data assimilation (MDA) Brief recap of MDA condition (ensuring correct sampling in linear-Gaussian case) Generalization of MDA condition to linearly dependent data vectors (PMDA condition) PMDA condition in practice - some issues

  5. Linearly dependent data vectors—example Data grid

  6. Linearly dependent data vectors—example l = L Data grid

  7. Linearly dependent data vectors—example Multilevel data l = L l = L − 1 l = L − 2 Data . . . grid { d l = B l d L } L − 1 l =1 With multilevel data, B l denotes an averaging operator from level L to level l

  8. Linearly dependent data vectors—example Multilevel data l = L l = L − 1 l = L − 2 Data . . . grid { d l = B l d L } L − 1 l =1 With multilevel data, B l denotes an averaging operator from level L to level l Time-domain multilevel data is also a possibility

  9. Multilevel data Why bother? l = L l = L − 1 l = L − 2 Data . . . grid { d l = B l d L } L − 1 l =1

  10. Multilevel data Why bother? l = L l = L − 1 l = L − 2 Data . . . grid { d l = B l d L } L − 1 l =1 Gradually introducing more and more information, like with sequential assimilation of d 1 , d 2 , . . . , d L , can be advantageous for nonlinear problems

  11. Multilevel data Why bother? l = L l = L − 1 l = L − 2 Data . . . grid { d l = B l d L } L − 1 l =1 Gradually introducing more and more information, like with sequential assimilation of d 1 , d 2 , . . . , d L , can be advantageous for nonlinear problems Multilevel data are required in order to correspond to results from multilevel simulations

  12. Multilevel simulations . . . Sim. output grid E E E

  13. Multilevel simulations . . . and corresponding multilevel data Data grid . . . Sim. output grid E E E

  14. Outline Motivation for considering linearly dependent data vectors Relation to multiple data assimilation (MDA) Brief recap of MDA condition (ensuring correct sampling in linear-Gaussian case) Generalization of MDA condition to linearly dependent data vectors (PMDA condition) PMDA condition in practice - some issues

  15. Multiple data assimilation 1 (MDA) Brief description With MDA, the same data are assimilated multiple times. Since the data are reused, the data-error covariances must be inflated. The motivation for MDA is to improve performance on nonlinear problems by gradually introducing the available information in the data, leading to a sequence of smaller updates instead of a single large update 1 Emerick and Reynolds, Computers & Geosci 55 , 2013

  16. MDA Multiple data assimilation { d l } L l =1 { d l = d L } L − 1 l =1 Multiple use of the same information Abbreviation: MDA

  17. MDA . . . as a special case of assimilation of multiple linearly related data vectors Multiple data Assimilation of multiple assimilation linearly related data vectors { d l } L { d l } L l =1 l =1 { d l = d L } L − 1 { d l = B l d L } L − 1 l =1 l =1 Multiple use of Partially multiple use of the same information the same information Abbreviation: MDA Abbreviation: PMDA (Partially MDA)

  18. Outline Motivation for considering linearly dependent data vectors Relation to multiple data assimilation (MDA) Brief recap of MDA condition (ensuring correct sampling in linear-Gaussian case) Generalization of MDA condition to linearly dependent data vectors (PMDA condition) PMDA condition in practice - some issues

  19. MDA condition Brief recap While the motivation for MDA is to improve performance on nonlinear problems, it is desirable that it samples correctly from the posterior PDF for the parameter vector in the linear-Gaussian case

  20. MDA condition Brief recap While the motivation for MDA is to improve performance on nonlinear problems, it is desirable that it samples correctly from the posterior PDF for the parameter vector in the linear-Gaussian case. This case can be analyzed using assembled quantities, where each row corresponds to an assimilation cycle  d L  . . . 0  C L  . δ = .   . . ... .   Ξ = . .   . .   d L 0 . . . C L G L   . . Γ =   .   G L

  21. MDA condition Brief recap While the motivation for MDA is to improve performance on nonlinear problems, it is desirable that it samples correctly from the posterior PDF for the parameter vector, m , in the linear-Gaussian case. This case can be analyzed using assembled quantities, where each row corresponds to an assimilation cycle. The analysis 2 leads to an inflated assembled covariance and the MDA condition for the inflation coefficients  d L  α 1 C L . . . 0   . δ = .   . . ... .   Ξ = . .   . .   d L 0 . . . α L C L G L   . Γ = .   . � L l =1 α − 1   = 1 l G L 2 Emerick and Reynolds, Computers & Geosci 55 , 2013

  22. Outline Motivation for considering linearly dependent data vectors Relation to multiple data assimilation (MDA) Brief recap of MDA condition (ensuring correct sampling in linear-Gaussian case) Generalization of MDA condition to linearly dependent data vectors (PMDA condition) PMDA condition in practice - some issues

  23. MDA condition Slight change of notation To prepare for the description of the PMDA condition, which follows next, I use the subscript MDA for ‘MDA quantities’  d L  α 1 C L . . . 0   . δ MDA = .   . . ... .   Ξ MDA = . .   . .   d L 0 . . . α L C L G L   . . Γ MDA =   . � L l =1 α − 1   = 1 l G L

  24. MDA condition Slight change of notation To prepare for the description of the PMDA condition, which follows next, I use the subscript MDA for ‘MDA quantities’, I introduce the coefficients { λ l = α 1 / 2 } L l l =1  d L  λ 2 1 C L . . . 0   . δ MDA = .   . . ... . . .   Ξ MDA =   . .   d L λ 2 0 . . . L C L G L   . . Γ MDA =   . � − 1 = 1   � L λ 2 � l =1 l G L

  25. MDA condition Slight change of notation To prepare for the description of the PMDA condition, which follows next, I use the subscript MDA for ‘MDA quantities’, I introduce the coefficients { λ l = α 1 / 2 l =1 , I multiply the MDA condition by C − 1 } L l L  d L  λ 2 1 C L . . . 0   . δ MDA = .   . . ... . Ξ MDA = . .     . .   d L λ 2 0 . . . L C L G L   . . Γ MDA =   . � − 1 = C − 1 C − 1   � L λ 2 � l =1 L l L G L

  26. MDA condition Slight change of notation To prepare for the description of the PMDA condition, which follows next, I use the subscript MDA for ‘MDA quantities’, I introduce the coefficients { λ l = α 1 / 2 l =1 , I multiply the MDA condition by C − 1 } L L , and I l reformulate the assembled data covariance and the MDA condition slightly  d L  λ 1 C L λ 1 . . . 0   . δ MDA = .   . . ... .   Ξ MDA = . .   . .   d L 0 . . . λ L C L λ L G L   . . Γ MDA =   . l =1 ( λ l C L λ l ) − 1 = C − 1 � L   L G L

  27. MDA condition  d L   λ 1 C L λ 1 . . . 0  . . . δ MDA = . ...   Ξ MDA = . . .     . .   d L 0 . . . λ L C L λ L  G L  . Γ MDA = .   .   l =1 ( λ l C L λ l ) − 1 = C − 1 � L G L L

  28. MDA condition  d L   λ 1 C L λ 1 . . . 0  . . . δ MDA = . ...   Ξ MDA = . . .     . .   d L 0 . . . λ L C L λ L  G L  . Γ MDA = .   .   l =1 ( λ l C L λ l ) − 1 = C − 1 � L G L L d 1   . δ PMDA = .   .   d L  G 1  . Γ PMDA = .   .   G L

  29. MDA condition and PMDA condition  d L   λ 1 C L λ 1 . . . 0  . . . δ MDA = . ...   Ξ MDA = . . .     . .   d L 0 . . . λ L C L λ L  G L  . Γ MDA = .   .   l =1 ( λ l C L λ l ) − 1 = C − 1 � L G L L d 1   A 1 C 1 A T  . . . 0  . 1 δ PMDA = .   . . ... . Ξ PMDA = . .     . .   d L A L C L A T 0 . . . L  G 1  . Γ PMDA = .   . � − 1 B l = C − 1   � L l =1 B T A l C l A T � l l L G L

Recommend


More recommend