A data-driven method for improving the correlation estimation in serial ensemble Kalman filter
|
|
- Candice Barber
- 5 years ago
- Views:
Transcription
1 A data-driven method for improving the correlation estimation in serial ensemble Kalman filter Michèle De La Chevrotière, 1 John Harlim 2 1 Department of Mathematics, Penn State University, 2 Department of Mathematics and Department of Meteorology, Penn State University The 7th EnKF Data Assimilation Workshop May 16
2 Outline 1. Correlation statistics in EnKF 2. Learning algorithm for improving the correlation estimation 3. Serial EnKF with transformed correlations 4. Numerical experiments on the Lorenz-96 model. Conclusion Michèle De La Chevrotière Data-driven method for correlation estimation in serial EnKF 2 / 18
3 Correlation statistics in EnKF Michèle De La Chevrotière Data-driven method for correlation estimation in serial EnKF 3 / 18
4 EnKF in linear and Gaussian settings In the linear and Gaussian, and under certain assumptions (controllability, observability, uncorrelated noises), the Kalman posterior error covariance matrix P a will converge to the fixed point Ĉ of Ĉ = (I KH)(F ĈF + Q), where K = ĈH (HĈH + R) 1 is the Kalman gain, F is the deterministic linear forecast model with noise covariance matrix Q, H is the linear observation operator, and R is the observation error covariance matrix. Michèle De La Chevrotière Data-driven method for correlation estimation in serial EnKF 4 / 18
5 EnKF in linear and Gaussian settings In the linear and Gaussian, and under certain assumptions (controllability, observability, uncorrelated noises), the Kalman posterior error covariance matrix P a will converge to the fixed point Ĉ of Ĉ = (I KH)(F ĈF + Q), where K = ĈH (HĈH + R) 1 is the Kalman gain, F is the deterministic linear forecast model with noise covariance matrix Q, H is the linear observation operator, and R is the observation error covariance matrix. EnKF construct an ensemble of analysis solutions {x a,k } K k=1 so that the ensemble mean and covariance match the Kalman solution. Michèle De La Chevrotière Data-driven method for correlation estimation in serial EnKF 4 / 18
6 EnKF in linear and Gaussian settings In the linear and Gaussian, and under certain assumptions (controllability, observability, uncorrelated noises), the Kalman posterior error covariance matrix P a will converge to the fixed point Ĉ of Ĉ = (I KH)(F ĈF + Q), where K = ĈH (HĈH + R) 1 is the Kalman gain, F is the deterministic linear forecast model with noise covariance matrix Q, H is the linear observation operator, and R is the observation error covariance matrix. EnKF construct an ensemble of analysis solutions {x a,k } K k=1 so that the ensemble mean and covariance match the Kalman solution. Here x a R N, and K is the ensemble size. Michèle De La Chevrotière Data-driven method for correlation estimation in serial EnKF 4 / 18
7 EnKF in linear and Gaussian settings In the linear and Gaussian, and under certain assumptions (controllability, observability, uncorrelated noises), the Kalman posterior error covariance matrix P a will converge to the fixed point Ĉ of Ĉ = (I KH)(F ĈF + Q), where K = ĈH (HĈH + R) 1 is the Kalman gain, F is the deterministic linear forecast model with noise covariance matrix Q, H is the linear observation operator, and R is the observation error covariance matrix. EnKF construct an ensemble of analysis solutions {x a,k } K k=1 so that the ensemble mean and covariance match the Kalman solution. Here x a R N, and K is the ensemble size. The scheme approximates both prior and posterior covariances by the ensemble covariance matrix, e.g., P a 1 K 1 XXT, where X is the matrix R N K of analysis ensemble perturbations from the ensemble mean (the k th column of X is X a,k = x a,k x a,k ). Michèle De La Chevrotière Data-driven method for correlation estimation in serial EnKF 4 / 18
8 EnKF in linear and Gaussian settings In the linear and Gaussian, and under certain assumptions (controllability, observability, uncorrelated noises), the Kalman posterior error covariance matrix P a will converge to the fixed point Ĉ of Ĉ = (I KH)(F ĈF + Q), where K = ĈH (HĈH + R) 1 is the Kalman gain, F is the deterministic linear forecast model with noise covariance matrix Q, H is the linear observation operator, and R is the observation error covariance matrix. EnKF construct an ensemble of analysis solutions {x a,k } K k=1 so that the ensemble mean and covariance match the Kalman solution. Here x a R N, and K is the ensemble size. The scheme approximates both prior and posterior covariances by the ensemble covariance matrix, e.g., P a 1 K 1 XXT, where X is the matrix R N K of analysis ensemble perturbations from the ensemble mean (the k th column of X is X a,k = x a,k x a,k ). In the linear and Gaussian setting, the EnKF solution is optimal with finite ensemble if x a,k t x a t N (0, Ĉ) in the limit as t. Michèle De La Chevrotière Data-driven method for correlation estimation in serial EnKF 4 / 18
9 Motivation Nonlinear and non-gaussian case: estimating the correlation ρ = D 1/2 CD 1/2, D = diag(c), is much more difficult because ρ is time-dependent and does not equilibrate to an invariant ρeq. Furthermore, the EnKF sample correlation, rk, is such that r (t) lim rk (t) ρ(t), t. K Miche le De La Chevrotie re Data-driven method for correlation estimation in serial EnKF / 18
10 Motivation Nonlinear and non-gaussian case: estimating the correlation ρ = D 1/2 CD 1/2, D = diag(c), is much more difficult because ρ is time-dependent and does not equilibrate to an invariant ρeq. Furthermore, the EnKF sample correlation, rk, is such that r (t) lim rk (t) ρ(t), t. K Since estimating the true correlation ρ is very difficult, then instead, we try to estimate r. Miche le De La Chevrotie re Data-driven method for correlation estimation in serial EnKF / 18
11 Motivation Nonlinear and non-gaussian case: estimating the correlation ρ = D 1/2 CD 1/2, D = diag(c), is much more difficult because ρ is time-dependent and does not equilibrate to an invariant ρeq. Furthermore, the EnKF sample correlation, rk, is such that r (t) lim rk (t) ρ(t), t. K Since estimating the true correlation ρ is very difficult, then instead, we try to estimate r. How can we transform rk into r? Miche le De La Chevrotie re Data-driven method for correlation estimation in serial EnKF / 18
12 Motivation Nonlinear and non-gaussian case: estimating the correlation ρ = D 1/2 CD 1/2, D = diag(c), is much more difficult because ρ is time-dependent and does not equilibrate to an invariant ρeq. Furthermore, the EnKF sample correlation, rk, is such that r (t) lim rk (t) ρ(t), t. K Since estimating the true correlation ρ is very difficult, then instead, we try to estimate r. How can we transform rk into r? The idea: I. Find a linear map L that transforms the undersampled correlation matrix rk into a sample of improved correlation matrix: LT rk r, Miche le De La Chevrotie re L RN N M. Data-driven method for correlation estimation in serial EnKF / 18
13 Motivation Nonlinear and non-gaussian case: estimating the correlation ρ = D 1/2 CD 1/2, D = diag(c), is much more difficult because ρ is time-dependent and does not equilibrate to an invariant ρeq. Furthermore, the EnKF sample correlation, rk, is such that r (t) lim rk (t) ρ(t), t. K Since estimating the true correlation ρ is very difficult, then instead, we try to estimate r. How can we transform rk into r? The idea: I. Find a linear map L that transforms the undersampled correlation matrix rk into a sample of improved correlation matrix: L RN N M State Variables State Variables State Variables Miche le De La Chevrotie re State Variables State Variables State Variables r r ^ rl ^ r rl rl^ LTrKLTrK LTrK rk rk 1 Observation Variables 1 Observation Variables Observation Variables rk r L r EnKF L= LT rk r, State Variables State Variables State Variables Data-driven method for correlation estimation in serial EnKF / 18
14 Motivation Nonlinear and non-gaussian case: estimating the correlation ρ = D 1/2 CD 1/2, D = diag(c), is much more difficult because ρ is time-dependent and does not equilibrate to an invariant ρeq. Furthermore, the EnKF sample correlation, rk, is such that r (t) lim rk (t) ρ(t), t. K Since estimating the true correlation ρ is very difficult, then instead, we try to estimate r. How can we transform rk into r? The idea: I. Find a linear map L that transforms the undersampled correlation matrix rk into a sample of improved correlation matrix: L RN N M State Variables State Variables State Variables State Variables State Variables State Variables r r ^ rl ^ r rl rl^ LTrKLTrK LTrK rk rk 1 Observation Variables 1 Observation Variables Observation Variables rk r L r EnKF L= LT rk r, State Variables State Variables State Variables II. Implement the map L into the serial LLS EnKF of Anderson (03). Miche le De La Chevrotie re Data-driven method for correlation estimation in serial EnKF / 18
15 Learning algorithm for improving the correlation estimation Michèle De La Chevrotière Data-driven method for correlation estimation in serial EnKF 6 / 18
16 A minimization problem For nonlinear filtering problems, where h(x) R N R M, we consider the sample cross correlation: r K xy = D 1/2 X ( 1 K 1 XY )D 1/2 Y R N M. [The K th column of Y is Y k = h(x a,k ) y a, where y a = K k=1 h(x a,k )/K] Michèle De La Chevrotière Data-driven method for correlation estimation in serial EnKF 7 / 18
17 A minimization problem For nonlinear filtering problems, where h(x) R N R M, we consider the sample cross correlation: r K xy = D 1/2 X ( 1 K 1 XY )D 1/2 Y R N M. [The K th column of Y is Y k = h(x a,k ) y a, where y a = K k=1 h(x a,k )/K] Onward: suppress the subscript xy (write r K xy as r K ; ˆr = lim K r K ). Michèle De La Chevrotière Data-driven method for correlation estimation in serial EnKF 7 / 18
18 A minimization problem For nonlinear filtering problems, where h(x) R N R M, we consider the sample cross correlation: r K xy = D 1/2 X ( 1 K 1 XY )D 1/2 Y R N M. [The K th column of Y is Y k = h(x a,k ) y a, where y a = K k=1 h(x a,k )/K] Onward: suppress the subscript xy (write r K xy as r K ; ˆr = lim K r K ). We seek to find a linear map L such that for every pair (i, j), L(, i, j) takes the sample correlations r K (, y j) to the limiting correlation ˆr(x i, y j). Michèle De La Chevrotière Data-driven method for correlation estimation in serial EnKF 7 / 18
19 A minimization problem For nonlinear filtering problems, where h(x) R N R M, we consider the sample cross correlation: r K xy = D 1/2 X ( 1 K 1 XY )D 1/2 Y R N M. [The K th column of Y is Y k = h(x a,k ) y a, where y a = K k=1 h(x a,k )/K] Onward: suppress the subscript xy (write r K xy as r K ; ˆr = lim K r K ). We seek to find a linear map L such that for every pair (i, j), L(, i, j) takes the sample correlations r K (, y j) to the limiting correlation ˆr(x i, y j). Inspired by Anderson (12), we formulate: min (r K (, y j) L(, i, j) ˆr(x i, y j)) 2 p(r K ˆr)p(ˆr) dr K dˆr. L(,i,j) [ 1,1] [ 1,1] Michèle De La Chevrotière Data-driven method for correlation estimation in serial EnKF 7 / 18
20 A minimization problem For nonlinear filtering problems, where h(x) R N R M, we consider the sample cross correlation: r K xy = D 1/2 X ( 1 K 1 XY )D 1/2 Y R N M. [The K th column of Y is Y k = h(x a,k ) y a, where y a = K k=1 h(x a,k )/K] Onward: suppress the subscript xy (write r K xy as r K ; ˆr = lim K r K ). We seek to find a linear map L such that for every pair (i, j), L(, i, j) takes the sample correlations r K (, y j) to the limiting correlation ˆr(x i, y j). Inspired by Anderson (12), we formulate: min (r K (, y j) L(, i, j) ˆr(x i, y j)) 2 p(r K ˆr)p(ˆr) dr K dˆr. L(,i,j) [ 1,1] [ 1,1] In practice, can use any reliable correlation statistics r L to approximate the limiting correlation ˆr, so that r L p(ˆr). (e.g. EnKF with large ensemble L). Michèle De La Chevrotière Data-driven method for correlation estimation in serial EnKF 7 / 18
21 Learning the map L 1. Train offline: Generate an OSSE assimilation analysis ensemble {x a,k t with L 1 and T assimilation cycles. } T,L t,k=1 Michèle De La Chevrotière Data-driven method for correlation estimation in serial EnKF 8 / 18
22 Learning the map L 1. Train offline: Generate an OSSE assimilation analysis ensemble {x a,k t with L 1 and T assimilation cycles. 2. Compute {r L t } T t=1, so that r L t p(ˆr). } T,L t,k=1 Michèle De La Chevrotière Data-driven method for correlation estimation in serial EnKF 8 / 18
23 Learning the map L 1. Train offline: Generate an OSSE assimilation analysis ensemble {x a,k t with L 1 and T assimilation cycles. 2. Compute {r L t } T t=1, so that r L t p(ˆr). 3. Compute {r K t } T t=1, subsampling K out of L members, so that r K t p(r K ˆr = r L t ). } T,L t,k=1 Michèle De La Chevrotière Data-driven method for correlation estimation in serial EnKF 8 / 18
24 Learning the map L 1. Train offline: Generate an OSSE assimilation analysis ensemble {x a,k t with L 1 and T assimilation cycles. 2. Compute {r L t } T t=1, so that r L t p(ˆr). 3. Compute {r K t } T t=1, subsampling K out of L members, so that r K t p(r K ˆr = r L t ). 4. Monte Carlo approximation to the minimization problem: } T,L t,k=1 or 1 min L(,i,j) T T t=1 (r K t (, y j) L(, i, j) r L t (x i, y j) ) 2. row of A u entry of b LSP min u Au b 2 when u = (A A) 1 A b Michèle De La Chevrotière Data-driven method for correlation estimation in serial EnKF 8 / 18
25 Learning the map L 1. Train offline: Generate an OSSE assimilation analysis ensemble {x a,k t with L 1 and T assimilation cycles. 2. Compute {r L t } T t=1, so that r L t p(ˆr). 3. Compute {r K t } T t=1, subsampling K out of L members, so that r K t p(r K ˆr = r L t ). 4. Monte Carlo approximation to the minimization problem: } T,L t,k=1 or 1 min L(,i,j) T T t=1 (r K t (, y j) L(, i, j) r L t (x i, y j) ) 2. row of A u entry of b LSP min u Au b 2 when u = (A A) 1 A b Offline training, little computational overhead, almost no tuning (K, T, L). Michèle De La Chevrotière Data-driven method for correlation estimation in serial EnKF 8 / 18
26 Learning the map L 1. Train offline: Generate an OSSE assimilation analysis ensemble {x a,k t with L 1 and T assimilation cycles. 2. Compute {r L t } T t=1, so that r L t p(ˆr). 3. Compute {r K t } T t=1, subsampling K out of L members, so that r K t p(r K ˆr = r L t ). 4. Monte Carlo approximation to the minimization problem: } T,L t,k=1 or 1 min L(,i,j) T T t=1 (r K t (, y j) L(, i, j) r L t (x i, y j) ) 2. row of A u entry of b LSP min u Au b 2 when u = (A A) 1 A b Offline training, little computational overhead, almost no tuning (K, T, L). L is a scalar: the map found can be viewed as a type of data-driven localization function (L d ). Michèle De La Chevrotière Data-driven method for correlation estimation in serial EnKF 8 / 18
27 Serial EnKF with transformed correlations Michèle De La Chevrotière Data-driven method for correlation estimation in serial EnKF 9 / 18
28 Adapted serial LLS EnKF (Anderson, 13) LLS EnKF assimilates observations j = 1,..., M one at a time: 1 for j = 1 to M do /* Loop over observations */ 2 y a j = y b j + (Py b j y j + R jj ) 1 Py b j y j (yj o yb j) ; 3 y a,k j = y a Rjj j + R jj (y b,k +Py b j y j y b j j); 4 y k j = ya,k j y b,k j ; /* compute the obs update */ x a,k = x b,k + P xy b j y P k b y j y j j 6 end ; /* regresses the update onto x */ The straightforward modification is to replace in line : P b xy j = D 1/2 X r(, y j)d 1/2 Y D 1/2 X L(,, j) r(, y j )D 1/2 Y, where D X, D Y and r are estimated form the EnKF with ensemble size K. Michèle De La Chevrotière Data-driven method for correlation estimation in serial EnKF / 18
29 Numerical experiments on the Lorenz-96 model Michèle De La Chevrotière Data-driven method for correlation estimation in serial EnKF 11 / 18
30 Linear direct observations Maps L and L d found by regressing onto ETKF L00 N = 40 state variables, (top) and 40 (bottom) observations observation time steps 1 (left) and (right) 40 obs obs L., L., L., L., obs obs L., L., L., L., n = 1 n = (a) ensemble members. n = 1 n = (b) 40 ensemble members. Localized structure of the mappings; L d (dashed black) is an envelope for the mappings L. Michèle De La Chevrotière Data-driven method for correlation estimation in serial EnKF 12 / 18
31 Linear direct observations En En nn = = 1 (c) L En n=en =1 nen n= En En En En= n En n =1 nen = En En En En obs obs En En En 1. EnEn En obs obs En En obs0. 1 En En En obs obs En En obs obs Filtering results with adapted LLS EnKF + covariance inflation nen =En n= (d) Ld. Time mean RMSE normalized by the RMSE of the ETKF using 00 ensemble members shown for the serial EnKF using L and Ld, plotted against inflation factor and K Miche le De La Chevrotie re Data-driven method for correlation estimation in serial EnKF 13 / 18
32 Linear direct observations Filtering results Best tuned inflation values observations, observation time step n = 1, comparison of 3 localization methods Time Mean of RMSE Time Mean of Spread ETKF L Ld GC Obs error Time mean of RMSE (left) and spread (right) as a function of ensemble size GC is more performant for K > but is extensively tuned! RMSE.0 2. Spread Michèle De La Chevrotière Data-driven method for correlation estimation in serial EnKF 14 /
33 Linear direct observations Regression using different products observations, observation time step n = 1, all use best tuned inflation values Time Mean of RMSE Time Mean of Spread ETKF L (against ETKF) Ld (against ETKF) L (against GC) Ld (against GC) GC 1 Ensemble Size 1 Ensemble Size Two choices of regression products: (1) ETKF with 00 ensemble members and (2) serial EnKF with GC localization with optimal half-width. Michèle De La Chevrotière Data-driven method for correlation estimation in serial EnKF 1 / 18
34 Linear indirect observations Maps L and L d found by regressing onto ETKF L00 observations, ensemble members obs L., L., L., n = 1 n = n = Michèle De La Chevrotière Data-driven method for correlation estimation in serial EnKF 16 / 18
35 Time M Linear indirect observations Filtering results with localization and best tuned inflation values observations, observation time step n = Time M Time Mean of RMSE Time Mean of Spread ETKF L Ld GC Obs error Time Mean of Spread Ensemble Size Ensemble Size e Mean of Spread Time mean of RMSE (left) and spread (right) Michèle De La Chevrotière Data-driven method for correlation estimation in serial EnKF 17 / 18
36 Summary A data-driven method for improving the correlation estimation in serial ensemble Kalman filter is introduced. The method finds a linear map that transforms, at each assimilation cycle, the poorly estimated sample correlation into a sample of improved correlation. This map is obtained from an offline training procedure with almost no tuning as the solution of a linear regression problem that uses appropriate sample correlation statistics obtained from historical data assimilation product. In numerical tests with the Lorenz-96 model for ranges of cases of linear and nonlinear observation models, the proposed scheme improves the filter estimates, especially when ensemble size is small relative to the dimension of the state space. Reference: M. De La Chevrotière, J. Harlim, A data-driven method for improving the correlation estimation in serial ensemble Kalman filter (16), Mon. Wea. Rev., submitted. Michèle De La Chevrotière Data-driven method for correlation estimation in serial EnKF 18 / 18
An Efficient Ensemble Data Assimilation Approach To Deal With Range Limited Observation
An Efficient Ensemble Data Assimilation Approach To Deal With Range Limited Observation A. Shah 1,2, M. E. Gharamti 1, L. Bertino 1 1 Nansen Environmental and Remote Sensing Center 2 University of Bergen
More informationDART_LAB Tutorial Section 2: How should observations impact an unobserved state variable? Multivariate assimilation.
DART_LAB Tutorial Section 2: How should observations impact an unobserved state variable? Multivariate assimilation. UCAR 2014 The National Center for Atmospheric Research is sponsored by the National
More informationLocalization and Correlation in Ensemble Kalman Filters
Localization and Correlation in Ensemble Kalman Filters Jeff Anderson NCAR Data Assimilation Research Section The National Center for Atmospheric Research is sponsored by the National Science Foundation.
More informationGaussian Filtering Strategies for Nonlinear Systems
Gaussian Filtering Strategies for Nonlinear Systems Canonical Nonlinear Filtering Problem ~u m+1 = ~ f (~u m )+~ m+1 ~v m+1 = ~g(~u m+1 )+~ o m+1 I ~ f and ~g are nonlinear & deterministic I Noise/Errors
More informationA Note on the Particle Filter with Posterior Gaussian Resampling
Tellus (6), 8A, 46 46 Copyright C Blackwell Munksgaard, 6 Printed in Singapore. All rights reserved TELLUS A Note on the Particle Filter with Posterior Gaussian Resampling By X. XIONG 1,I.M.NAVON 1,2 and
More informationData assimilation in high dimensions
Data assimilation in high dimensions David Kelly Courant Institute New York University New York NY www.dtbkelly.com February 12, 2015 Graduate seminar, CIMS David Kelly (CIMS) Data assimilation February
More informationAdvances and Challenges in Ensemblebased Data Assimilation in Meteorology. Takemasa Miyoshi
January 18, 2013, DA Workshop, Tachikawa, Japan Advances and Challenges in Ensemblebased Data Assimilation in Meteorology Takemasa Miyoshi RIKEN Advanced Institute for Computational Science Takemasa.Miyoshi@riken.jp
More informationData Assimilation Research Testbed Tutorial
Data Assimilation Research Testbed Tutorial Section 3: Hierarchical Group Filters and Localization Version 2.: September, 26 Anderson: Ensemble Tutorial 9//6 Ways to deal with regression sampling error:
More informationData Assimilation Research Testbed Tutorial
Data Assimilation Research Testbed Tutorial Section 2: How should observations of a state variable impact an unobserved state variable? Multivariate assimilation. Single observed variable, single unobserved
More informationEnsemble Kalman Filter
Ensemble Kalman Filter Geir Evensen and Laurent Bertino Hydro Research Centre, Bergen, Norway, Nansen Environmental and Remote Sensing Center, Bergen, Norway The Ensemble Kalman Filter (EnKF) Represents
More informationA new Hierarchical Bayes approach to ensemble-variational data assimilation
A new Hierarchical Bayes approach to ensemble-variational data assimilation Michael Tsyrulnikov and Alexander Rakitko HydroMetCenter of Russia College Park, 20 Oct 2014 Michael Tsyrulnikov and Alexander
More informationAsynchronous data assimilation
Ensemble Kalman Filter, lecture 2 Asynchronous data assimilation Pavel Sakov Nansen Environmental and Remote Sensing Center, Norway This talk has been prepared in the course of evita-enkf project funded
More informationAddressing the nonlinear problem of low order clustering in deterministic filters by using mean-preserving non-symmetric solutions of the ETKF
Addressing the nonlinear problem of low order clustering in deterministic filters by using mean-preserving non-symmetric solutions of the ETKF Javier Amezcua, Dr. Kayo Ide, Dr. Eugenia Kalnay 1 Outline
More informationCorrecting biased observation model error in data assimilation
Correcting biased observation model error in data assimilation Tyrus Berry Dept. of Mathematical Sciences, GMU PSU-UMD DA Workshop June 27, 217 Joint work with John Harlim, PSU BIAS IN OBSERVATION MODELS
More informationState and Parameter Estimation in Stochastic Dynamical Models
State and Parameter Estimation in Stochastic Dynamical Models Timothy DelSole George Mason University, Fairfax, Va and Center for Ocean-Land-Atmosphere Studies, Calverton, MD June 21, 2011 1 1 collaboration
More informationReview of Covariance Localization in Ensemble Filters
NOAA Earth System Research Laboratory Review of Covariance Localization in Ensemble Filters Tom Hamill NOAA Earth System Research Lab, Boulder, CO tom.hamill@noaa.gov Canonical ensemble Kalman filter update
More informationKalman Filter and Ensemble Kalman Filter
Kalman Filter and Ensemble Kalman Filter 1 Motivation Ensemble forecasting : Provides flow-dependent estimate of uncertainty of the forecast. Data assimilation : requires information about uncertainty
More informationPractical Aspects of Ensemble-based Kalman Filters
Practical Aspects of Ensemble-based Kalman Filters Lars Nerger Alfred Wegener Institute for Polar and Marine Research Bremerhaven, Germany and Bremen Supercomputing Competence Center BremHLR Bremen, Germany
More informationLocal Ensemble Transform Kalman Filter
Local Ensemble Transform Kalman Filter Brian Hunt 11 June 2013 Review of Notation Forecast model: a known function M on a vector space of model states. Truth: an unknown sequence {x n } of model states
More informationEnsemble Data Assimilation and Uncertainty Quantification
Ensemble Data Assimilation and Uncertainty Quantification Jeff Anderson National Center for Atmospheric Research pg 1 What is Data Assimilation? Observations combined with a Model forecast + to produce
More informationErgodicity in data assimilation methods
Ergodicity in data assimilation methods David Kelly Andy Majda Xin Tong Courant Institute New York University New York NY www.dtbkelly.com April 15, 2016 ETH Zurich David Kelly (CIMS) Data assimilation
More informationData assimilation with and without a model
Data assimilation with and without a model Tim Sauer George Mason University Parameter estimation and UQ U. Pittsburgh Mar. 5, 2017 Partially supported by NSF Most of this work is due to: Tyrus Berry,
More informationEstimating Observation Impact in a Hybrid Data Assimilation System: Experiments with a Simple Model
Estimating Observation Impact in a Hybrid Data Assimilation System: Experiments with a Simple Model Rahul Mahajan NOAA / NCEP / EMC - IMSG College Park, MD 20740, USA Acknowledgements: Ron Gelaro, Ricardo
More informationEstimating Observation Impact in a Hybrid Data Assimilation System: Experiments with a Simple Model
Estimating Observation Impact in a Hybrid Data Assimilation System: Experiments with a Simple Model NOAA / NCEP / EMC College Park, MD 20740, USA 10 February 2014 Overview Goal Sensitivity Theory Adjoint
More informationConvergence of Square Root Ensemble Kalman Filters in the Large Ensemble Limit
Convergence of Square Root Ensemble Kalman Filters in the Large Ensemble Limit Evan Kwiatkowski, Jan Mandel University of Colorado Denver December 11, 2014 OUTLINE 2 Data Assimilation Bayesian Estimation
More informationObjective localization of ensemble covariances: theory and applications
Institutionnel Grand Public Objective localization of ensemble covariances: theory and applications Yann Michel1, B. Me ne trier2 and T. Montmerle1 Professionnel (1) Me te o-france & CNRS, Toulouse, France
More informationAbstract 2. ENSEMBLE KALMAN FILTERS 1. INTRODUCTION
J5.4 4D ENSEMBLE KALMAN FILTERING FOR ASSIMILATION OF ASYNCHRONOUS OBSERVATIONS T. Sauer George Mason University, Fairfax, VA 22030 B.R. Hunt, J.A. Yorke, A.V. Zimin, E. Ott, E.J. Kostelich, I. Szunyogh,
More informationFinite-size ensemble Kalman filters (EnKF-N)
Finite-size ensemble Kalman filters (EnKF-N) Marc Bocquet (1,2), Pavel Sakov (3,4) (1) Université Paris-Est, CEREA, joint lab École des Ponts ParisTech and EdF R&D, France (2) INRIA, Paris-Rocquencourt
More informationEnKF Review. P.L. Houtekamer 7th EnKF workshop Introduction to the EnKF. Challenges. The ultimate global EnKF algorithm
Overview 1 2 3 Review of the Ensemble Kalman Filter for Atmospheric Data Assimilation 6th EnKF Purpose EnKF equations localization After the 6th EnKF (2014), I decided with Prof. Zhang to summarize progress
More informationAssimilation des Observations et Traitement des Incertitudes en Météorologie
Assimilation des Observations et Traitement des Incertitudes en Météorologie Olivier Talagrand Laboratoire de Météorologie Dynamique, Paris 4èmes Rencontres Météo/MathAppli Météo-France, Toulouse, 25 Mars
More informationReducing the Impact of Sampling Errors in Ensemble Filters
Reducing the Impact of Sampling Errors in Ensemble Filters Jeffrey Anderson NCAR Data Assimilation Research Section The National Center for Atmospheric Research is sponsored by the National Science Foundation.
More informationConvective-scale data assimilation in the Weather Research and Forecasting model using a nonlinear ensemble filter
Convective-scale data assimilation in the Weather Research and Forecasting model using a nonlinear ensemble filter Jon Poterjoy, Ryan Sobash, and Jeffrey Anderson National Center for Atmospheric Research
More informationForecasting and data assimilation
Supported by the National Science Foundation DMS Forecasting and data assimilation Outline Numerical models Kalman Filter Ensembles Douglas Nychka, Thomas Bengtsson, Chris Snyder Geophysical Statistics
More informationEnKF and Catastrophic filter divergence
EnKF and Catastrophic filter divergence David Kelly Kody Law Andrew Stuart Mathematics Department University of North Carolina Chapel Hill NC bkelly.com June 4, 2014 SIAM UQ 2014 Savannah. David Kelly
More informationJuli I. Rubin. NRC Postdoctoral Research Associate Naval Research Laboratory, Monterey, CA
Development of the Ensemble Navy Aerosol Analysis Prediction System and its application of the Data Assimilation Research Testbed in Support of Aerosol Forecasting Juli I. Rubin NRC Postdoctoral Research
More informationBayesian Inverse problem, Data assimilation and Localization
Bayesian Inverse problem, Data assimilation and Localization Xin T Tong National University of Singapore ICIP, Singapore 2018 X.Tong Localization 1 / 37 Content What is Bayesian inverse problem? What is
More informationGenerating climatological forecast error covariance for Variational DAs with ensemble perturbations: comparison with the NMC method
Generating climatological forecast error covariance for Variational DAs with ensemble perturbations: comparison with the NMC method Matthew Wespetal Advisor: Dr. Eugenia Kalnay UMD, AOSC Department March
More informationSmoothers: Types and Benchmarks
Smoothers: Types and Benchmarks Patrick N. Raanes Oxford University, NERSC 8th International EnKF Workshop May 27, 2013 Chris Farmer, Irene Moroz Laurent Bertino NERSC Geir Evensen Abstract Talk builds
More informationWhat do we know about EnKF?
What do we know about EnKF? David Kelly Kody Law Andrew Stuart Andrew Majda Xin Tong Courant Institute New York University New York, NY April 10, 2015 CAOS seminar, Courant. David Kelly (NYU) EnKF April
More informationData assimilation in high dimensions
Data assimilation in high dimensions David Kelly Kody Law Andy Majda Andrew Stuart Xin Tong Courant Institute New York University New York NY www.dtbkelly.com February 3, 2016 DPMMS, University of Cambridge
More informationAdaptive ensemble Kalman filtering of nonlinear systems
Adaptive ensemble Kalman filtering of nonlinear systems Tyrus Berry George Mason University June 12, 213 : Problem Setup We consider a system of the form: x k+1 = f (x k ) + ω k+1 ω N (, Q) y k+1 = h(x
More informationOrganization. I MCMC discussion. I project talks. I Lecture.
Organization I MCMC discussion I project talks. I Lecture. Content I Uncertainty Propagation Overview I Forward-Backward with an Ensemble I Model Reduction (Intro) Uncertainty Propagation in Causal Systems
More informationThe Ensemble Kalman Filter:
p.1 The Ensemble Kalman Filter: Theoretical formulation and practical implementation Geir Evensen Norsk Hydro Research Centre, Bergen, Norway Based on Evensen, Ocean Dynamics, Vol 5, No p. The Ensemble
More informationEnsemble Kalman Filter potential
Ensemble Kalman Filter potential Former students (Shu-Chih( Yang, Takemasa Miyoshi, Hong Li, Junjie Liu, Chris Danforth, Ji-Sun Kang, Matt Hoffman), and Eugenia Kalnay University of Maryland Acknowledgements:
More informationData assimilation with and without a model
Data assimilation with and without a model Tyrus Berry George Mason University NJIT Feb. 28, 2017 Postdoc supported by NSF This work is in collaboration with: Tim Sauer, GMU Franz Hamilton, Postdoc, NCSU
More informationA variance limiting Kalman filter for data assimilation: I. Sparse observational grids II. Model error
A variance limiting Kalman filter for data assimilation: I. Sparse observational grids II. Model error Georg Gottwald, Lewis Mitchell, Sebastian Reich* University of Sydney, *Universität Potsdam Durham,
More informationEnKF and filter divergence
EnKF and filter divergence David Kelly Andrew Stuart Kody Law Courant Institute New York University New York, NY dtbkelly.com December 12, 2014 Applied and computational mathematics seminar, NIST. David
More informationRobust Ensemble Filtering With Improved Storm Surge Forecasting
Robust Ensemble Filtering With Improved Storm Surge Forecasting U. Altaf, T. Buttler, X. Luo, C. Dawson, T. Mao, I.Hoteit Meteo France, Toulouse, Nov 13, 2012 Project Ensemble data assimilation for storm
More informationLocalization in the ensemble Kalman Filter
Department of Meteorology Localization in the ensemble Kalman Filter Ruth Elizabeth Petrie A dissertation submitted in partial fulfilment of the requirement for the degree of MSc. Atmosphere, Ocean and
More informationRelative Merits of 4D-Var and Ensemble Kalman Filter
Relative Merits of 4D-Var and Ensemble Kalman Filter Andrew Lorenc Met Office, Exeter International summer school on Atmospheric and Oceanic Sciences (ISSAOS) "Atmospheric Data Assimilation". August 29
More informationThe Ensemble Kalman Filter:
p.1 The Ensemble Kalman Filter: Theoretical formulation and practical implementation Geir Evensen Norsk Hydro Research Centre, Bergen, Norway Based on Evensen 23, Ocean Dynamics, Vol 53, No 4 p.2 The Ensemble
More informationEnvironment Canada s Regional Ensemble Kalman Filter
Environment Canada s Regional Ensemble Kalman Filter May 19, 2014 Seung-Jong Baek, Luc Fillion, Kao-Shen Chung, and Peter Houtekamer Meteorological Research Division, Environment Canada, Dorval, Quebec
More informationEnsemble square-root filters
Ensemble square-root filters MICHAEL K. TIPPETT International Research Institute for climate prediction, Palisades, New Yor JEFFREY L. ANDERSON GFDL, Princeton, New Jersy CRAIG H. BISHOP Naval Research
More informationEnsemble Kalman Filters for WRF-ARW. Chris Snyder MMM and IMAGe National Center for Atmospheric Research
Ensemble Kalman Filters for WRF-ARW Chris Snyder MMM and IMAGe National Center for Atmospheric Research Preliminaries Notation: x = modelʼs state w.r.t. some discrete basis, e.g. grid-pt values y = Hx
More informationShort tutorial on data assimilation
Mitglied der Helmholtz-Gemeinschaft Short tutorial on data assimilation 23 June 2015 Wolfgang Kurtz & Harrie-Jan Hendricks Franssen Institute of Bio- and Geosciences IBG-3 (Agrosphere), Forschungszentrum
More informationPerformance of ensemble Kalman filters with small ensembles
Performance of ensemble Kalman filters with small ensembles Xin T Tong Joint work with Andrew J. Majda National University of Singapore Sunday 28 th May, 2017 X.Tong EnKF performance 1 / 31 Content Ensemble
More informationLocal Ensemble Transform Kalman Filter: An Efficient Scheme for Assimilating Atmospheric Data
Local Ensemble Transform Kalman Filter: An Efficient Scheme for Assimilating Atmospheric Data John Harlim and Brian R. Hunt Department of Mathematics and Institute for Physical Science and Technology University
More informationA Spectral Approach to Linear Bayesian Updating
A Spectral Approach to Linear Bayesian Updating Oliver Pajonk 1,2, Bojana V. Rosic 1, Alexander Litvinenko 1, and Hermann G. Matthies 1 1 Institute of Scientific Computing, TU Braunschweig, Germany 2 SPT
More informationImplications of the Form of the Ensemble Transformation in the Ensemble Square Root Filters
1042 M O N T H L Y W E A T H E R R E V I E W VOLUME 136 Implications of the Form of the Ensemble Transformation in the Ensemble Square Root Filters PAVEL SAKOV AND PETER R. OKE CSIRO Marine and Atmospheric
More informationLagrangian Data Assimilation and Manifold Detection for a Point-Vortex Model. David Darmon, AMSC Kayo Ide, AOSC, IPST, CSCAMM, ESSIC
Lagrangian Data Assimilation and Manifold Detection for a Point-Vortex Model David Darmon, AMSC Kayo Ide, AOSC, IPST, CSCAMM, ESSIC Background Data Assimilation Iterative process Forecast Analysis Background
More informationFour-Dimensional Ensemble Kalman Filtering
Four-Dimensional Ensemble Kalman Filtering B.R. Hunt, E. Kalnay, E.J. Kostelich, E. Ott, D.J. Patil, T. Sauer, I. Szunyogh, J.A. Yorke, A.V. Zimin University of Maryland, College Park, MD 20742, USA Ensemble
More informationFundamentals of Data Assimilation
National Center for Atmospheric Research, Boulder, CO USA GSI Data Assimilation Tutorial - June 28-30, 2010 Acknowledgments and References WRFDA Overview (WRF Tutorial Lectures, H. Huang and D. Barker)
More informationGaussian Process Approximations of Stochastic Differential Equations
Gaussian Process Approximations of Stochastic Differential Equations Cédric Archambeau Centre for Computational Statistics and Machine Learning University College London c.archambeau@cs.ucl.ac.uk CSML
More informationSystematic strategies for real time filtering of turbulent signals in complex systems
Systematic strategies for real time filtering of turbulent signals in complex systems Statistical inversion theory for Gaussian random variables The Kalman Filter for Vector Systems: Reduced Filters and
More informationHow 4DVAR can benefit from or contribute to EnKF (a 4DVAR perspective)
How 4DVAR can benefit from or contribute to EnKF (a 4DVAR perspective) Dale Barker WWRP/THORPEX Workshop on 4D-Var and Ensemble Kalman Filter Intercomparisons Sociedad Cientifica Argentina, Buenos Aires,
More informationMethods of Data Assimilation and Comparisons for Lagrangian Data
Methods of Data Assimilation and Comparisons for Lagrangian Data Chris Jones, Warwick and UNC-CH Kayo Ide, UCLA Andrew Stuart, Jochen Voss, Warwick Guillaume Vernieres, UNC-CH Amarjit Budiraja, UNC-CH
More informationA Moment Matching Particle Filter for Nonlinear Non-Gaussian. Data Assimilation. and Peter Bickel
Generated using version 3.0 of the official AMS L A TEX template A Moment Matching Particle Filter for Nonlinear Non-Gaussian Data Assimilation Jing Lei and Peter Bickel Department of Statistics, University
More informationAn ETKF approach for initial state and parameter estimation in ice sheet modelling
An ETKF approach for initial state and parameter estimation in ice sheet modelling Bertrand Bonan University of Reading DARC Seminar, Reading October 30, 2013 DARC Seminar (Reading) ETKF for ice sheet
More informationGaussian Process Approximations of Stochastic Differential Equations
Gaussian Process Approximations of Stochastic Differential Equations Cédric Archambeau Dan Cawford Manfred Opper John Shawe-Taylor May, 2006 1 Introduction Some of the most complex models routinely run
More information(Extended) Kalman Filter
(Extended) Kalman Filter Brian Hunt 7 June 2013 Goals of Data Assimilation (DA) Estimate the state of a system based on both current and all past observations of the system, using a model for the system
More informationEnKF Localization Techniques and Balance
EnKF Localization Techniques and Balance Steven Greybush Eugenia Kalnay, Kayo Ide, Takemasa Miyoshi, and Brian Hunt Weather Chaos Meeting September 21, 2009 Data Assimilation Equation Scalar form: x a
More informationP 1.86 A COMPARISON OF THE HYBRID ENSEMBLE TRANSFORM KALMAN FILTER (ETKF)- 3DVAR AND THE PURE ENSEMBLE SQUARE ROOT FILTER (EnSRF) ANALYSIS SCHEMES
P 1.86 A COMPARISON OF THE HYBRID ENSEMBLE TRANSFORM KALMAN FILTER (ETKF)- 3DVAR AND THE PURE ENSEMBLE SQUARE ROOT FILTER (EnSRF) ANALYSIS SCHEMES Xuguang Wang*, Thomas M. Hamill, Jeffrey S. Whitaker NOAA/CIRES
More informationThe Kalman Filter. Data Assimilation & Inverse Problems from Weather Forecasting to Neuroscience. Sarah Dance
The Kalman Filter Data Assimilation & Inverse Problems from Weather Forecasting to Neuroscience Sarah Dance School of Mathematical and Physical Sciences, University of Reading s.l.dance@reading.ac.uk July
More informationEfficient Data Assimilation for Spatiotemporal Chaos: a Local Ensemble Transform Kalman Filter
Efficient Data Assimilation for Spatiotemporal Chaos: a Local Ensemble Transform Kalman Filter arxiv:physics/0511236 v1 28 Nov 2005 Brian R. Hunt Institute for Physical Science and Technology and Department
More information(Toward) Scale-dependent weighting and localization for the NCEP GFS hybrid 4DEnVar Scheme
(Toward) Scale-dependent weighting and localization for the NCEP GFS hybrid 4DEnVar Scheme Daryl Kleist 1, Kayo Ide 1, Rahul Mahajan 2, Deng-Shun Chen 3 1 University of Maryland - Dept. of Atmospheric
More informationEnsemble Data Assimila.on for Climate System Component Models
Ensemble Data Assimila.on for Climate System Component Models Jeffrey Anderson Na.onal Center for Atmospheric Research In collabora.on with: Alicia Karspeck, Kevin Raeder, Tim Hoar, Nancy Collins IMA 11
More informationEnsemble of Data Assimilations methods for the initialization of EPS
Ensemble of Data Assimilations methods for the initialization of EPS Laure RAYNAUD Météo-France ECMWF Annual Seminar Reading, 12 September 2017 Introduction Estimating the uncertainty in the initial conditions
More informationHierarchical Bayes Ensemble Kalman Filter
Hierarchical Bayes Ensemble Kalman Filter M Tsyrulnikov and A Rakitko HydroMetCenter of Russia Wrocław, 7 Sep 2015 M Tsyrulnikov and A Rakitko (HMC) Hierarchical Bayes Ensemble Kalman Filter Wrocław, 7
More informationEnhancing information transfer from observations to unobserved state variables for mesoscale radar data assimilation
Enhancing information transfer from observations to unobserved state variables for mesoscale radar data assimilation Weiguang Chang and Isztar Zawadzki Department of Atmospheric and Oceanic Sciences Faculty
More informationRelationship between Singular Vectors, Bred Vectors, 4D-Var and EnKF
Relationship between Singular Vectors, Bred Vectors, 4D-Var and EnKF Eugenia Kalnay and Shu-Chih Yang with Alberto Carrasi, Matteo Corazza and Takemasa Miyoshi 4th EnKF Workshop, April 2010 Relationship
More informationEnsemble 4DVAR for the NCEP hybrid GSI EnKF data assimilation system and observation impact study with the hybrid system
Ensemble 4DVAR for the NCEP hybrid GSI EnKF data assimilation system and observation impact study with the hybrid system Xuguang Wang School of Meteorology University of Oklahoma, Norman, OK OU: Ting Lei,
More informationPATTERN RECOGNITION AND MACHINE LEARNING CHAPTER 13: SEQUENTIAL DATA
PATTERN RECOGNITION AND MACHINE LEARNING CHAPTER 13: SEQUENTIAL DATA Contents in latter part Linear Dynamical Systems What is different from HMM? Kalman filter Its strength and limitation Particle Filter
More informationError (E) Covariance Diagnostics. in Ensemble Data Assimilation. Kayo Ide, Daryl Kleist, Adam Kubaryk University of Maryland
Error (E) Covariance Diagnostics in Ensemble Data Assimilation Kayo Ide, Daryl Kleist, Adam Kubaryk University of Maryland 7 th EnKF Data Assimilation Workshop May 2016 Outline Part I. Background: E-diagnostics
More informationIntroduction to Ensemble Kalman Filters and the Data Assimilation Research Testbed
Introduction to Ensemble Kalman Filters and the Data Assimilation Research Testbed Jeffrey Anderson, Tim Hoar, Nancy Collins NCAR Institute for Math Applied to Geophysics pg 1 What is Data Assimilation?
More informationDART Tutorial Sec'on 5: Comprehensive Filtering Theory: Non-Iden'ty Observa'ons and the Joint Phase Space
DART Tutorial Sec'on 5: Comprehensive Filtering Theory: Non-Iden'ty Observa'ons and the Joint Phase Space UCAR The Na'onal Center for Atmospheric Research is sponsored by the Na'onal Science Founda'on.
More informationToward improved initial conditions for NCAR s real-time convection-allowing ensemble. Ryan Sobash, Glen Romine, Craig Schwartz, and Kate Fossell
Toward improved initial conditions for NCAR s real-time convection-allowing ensemble Ryan Sobash, Glen Romine, Craig Schwartz, and Kate Fossell Storm-scale ensemble design Can an EnKF be used to initialize
More informationThe Canadian approach to ensemble prediction
The Canadian approach to ensemble prediction ECMWF 2017 Annual seminar: Ensemble prediction : past, present and future. Pieter Houtekamer Montreal, Canada Overview. The Canadian approach. What are the
More informationSome ideas for Ensemble Kalman Filter
Some ideas for Ensemble Kalman Filter Former students and Eugenia Kalnay UMCP Acknowledgements: UMD Chaos-Weather Group: Brian Hunt, Istvan Szunyogh, Ed Ott and Jim Yorke, Kayo Ide, and students Former
More informationLagrangian Data Assimilation and Its Application to Geophysical Fluid Flows
Lagrangian Data Assimilation and Its Application to Geophysical Fluid Flows Laura Slivinski June, 3 Laura Slivinski (Brown University) Lagrangian Data Assimilation June, 3 / 3 Data Assimilation Setup:
More informationModel error and parameter estimation
Model error and parameter estimation Chiara Piccolo and Mike Cullen ECMWF Annual Seminar, 11 September 2018 Summary The application of interest is atmospheric data assimilation focus on EDA; A good ensemble
More informationEnsemble forecasting and flow-dependent estimates of initial uncertainty. Martin Leutbecher
Ensemble forecasting and flow-dependent estimates of initial uncertainty Martin Leutbecher acknowledgements: Roberto Buizza, Lars Isaksen Flow-dependent aspects of data assimilation, ECMWF 11 13 June 2007
More informationFour-dimensional ensemble Kalman filtering
Tellus (24), 56A, 273 277 Copyright C Blackwell Munksgaard, 24 Printed in UK. All rights reserved TELLUS Four-dimensional ensemble Kalman filtering By B. R. HUNT 1, E. KALNAY 1, E. J. KOSTELICH 2,E.OTT
More informationEvolution of Forecast Error Covariances in 4D-Var and ETKF methods
Evolution of Forecast Error Covariances in 4D-Var and ETKF methods Chiara Piccolo Met Office Exeter, United Kingdom chiara.piccolo@metoffice.gov.uk Introduction Estimates of forecast error covariances
More informationApplication of the Ensemble Kalman Filter to History Matching
Application of the Ensemble Kalman Filter to History Matching Presented at Texas A&M, November 16,2010 Outline Philosophy EnKF for Data Assimilation Field History Match Using EnKF with Covariance Localization
More informationCan hybrid-4denvar match hybrid-4dvar?
Comparing ensemble-variational assimilation methods for NWP: Can hybrid-4denvar match hybrid-4dvar? WWOSC, Montreal, August 2014. Andrew Lorenc, Neill Bowler, Adam Clayton, David Fairbairn and Stephen
More informationOptimal Localization for Ensemble Kalman Filter Systems
Journal December of the 2014 Meteorological Society of Japan, Vol. Á. 92, PERIÁÑEZ No. 6, pp. et 585 597, al. 2014 585 DOI:10.2151/jmsj.2014-605 Optimal Localization for Ensemble Kalman Filter Systems
More information4DEnVar. Four-Dimensional Ensemble-Variational Data Assimilation. Colloque National sur l'assimilation de données
Four-Dimensional Ensemble-Variational Data Assimilation 4DEnVar Colloque National sur l'assimilation de données Andrew Lorenc, Toulouse France. 1-3 décembre 2014 Crown copyright Met Office 4DEnVar: Topics
More informationMaximum Likelihood Ensemble Filter Applied to Multisensor Systems
Maximum Likelihood Ensemble Filter Applied to Multisensor Systems Arif R. Albayrak a, Milija Zupanski b and Dusanka Zupanski c abc Colorado State University (CIRA), 137 Campus Delivery Fort Collins, CO
More informationA Stochastic Collocation based. for Data Assimilation
A Stochastic Collocation based Kalman Filter (SCKF) for Data Assimilation Lingzao Zeng and Dongxiao Zhang University of Southern California August 11, 2009 Los Angeles Outline Introduction SCKF Algorithm
More informationABSTRACT. Sabrina Rainwater Doctor of Philosophy, Professor Brian R. Hunt Department of Mathematics and
ABSTRACT Title of dissertation: NONLINEAR AND MULTIRESOLUTION ERROR COVARIANCE ESTIMATION IN ENSEMBLE DATA ASSIMILATION Sabrina Rainwater Doctor of Philosophy, 2012 Dissertation directed by: Professor
More information