Attribution of climatic events using data assimilation
|
|
- Sydney Chambers
- 5 years ago
- Views:
Transcription
1 Attribution of climatic events using data assimilation Alberto Carrassi Nansen Environmental and Remote Sensing Center - Norway University of Maryland College Park 24 th April 2018 A. Carrassi University of Maryland. USA, 24 th April / 33
2 The original goal: the attribution problem Using data assimilation as a tool to diagnose and evidence causal links in weather/climate events. Move toward an operational real-time event attribution system. Work main team: M. Bocquet (CEREA-ENPC, France), M. Ghil (ENS/UCLA, France/USA) and A. Hannart (OURANOS, Canada), J. Ruiz (Un. of Buenos Aires, Argentina), S. Metref (Un. of Buenos Aires, Argentina) Funding partly provided by the ANR project DADA ( ): Data Assimilation for Detection and Attribution of climate change A. Carrassi University of Maryland. USA, 24 th April / 33
3 Outline 1 Data assimilation in a nutshell 2 The attribution problem in climate science 3 Model evidence: a key quantity Definition and interpretation Contextual model evidence Computing model evidence using data assimilation 4 Model evidence: Numerical results Setup Accuracy of the DA-based estimates Low dimensional models - Application for attribution Model evidence with localization Application for model selection Application for the attribution problem (high-dim) 5 Conclusion and future directions A. Carrassi University of Maryland. USA, 24 th April / 33
4 Data assimilation in a nutshell Data assimilation in a nutshell environmental system Filtering observations forecast analysis The problem is in principle solved using a probabilistic framework The quantity of interest are the probability density functions, PDF The PDFs are evolved in time and updated at analysis times using Bayes s rules p(x y) = p(y x)p(x) p(y) model truth From Carrassi et al, t k t k+1 t k+2 time The huge dimensions of typical models and datasets hamper the use of a fully Bayesian approach in the geosciences. The Gaussian approximation is at the core of many successful DA procedures. A. Carrassi University of Maryland. USA, 24 th April / 33
5 Data assimilation in a nutshell Data assimilation in a nutshell A. Carrassi University of Maryland. USA, 24 th April / 33
6 The attribution problem in climate science The attribution problem in climate science A. Carrassi University of Maryland. USA, 24 th April / 33
7 The attribution problem in climate science The attribution problem: trends or events A. Carrassi University of Maryland. USA, 24 th April / 33
8 The attribution problem in climate science The attribution of a climatic event Do we identify a discernible human influence on a single weather or climate-related event? Implications for mitigation and adaptation policy, litigation in courts, insurance, physical understanding of the event. Define the event based on threshold exceedance of an ad-hoc climate index Derive the likelihood of the event in two model ensembles: factual world (e.g. including CO2 emission) = p 1 counterfactual world (e.g. without CO2 emission) = p 0 Derive the fraction of attributable risk (FAR) F AR = 1 p 0 p 1 A. Carrassi University of Maryland. USA, 24 th April / 33
9 The attribution problem in climate science The conventional method August 2003 European heatwave (Stott et al., 2004) They conclude: "It is very likely (> 90%) that CO2 emissions have increased the frequency of occurrence of 2003-like heatwaves by a factor of at least two." Issues with the conventional approach: It requires to run a large ensemble of model simulations, plus extrapolation using EVT for both the factual and counter-factual (very costly). Long delays in producing assessment after an event occurrence. A. Carrassi University of Maryland. USA, 24 th April / 33
10 The attribution problem in climate science The new near real-time event attribution The challenge is to design a system/procedure that allows for near real-time event attribution "The overarching challenge for the community is to move beyond research- mode case studies and to develop systems that can deliver regular, reliable and timely assessments in the aftermath of notable weather and climate- related events, typically in the weeks or months following (and not many years later as is the case with some research-mode)." Stott et al. (2013) Based on data assimilation, our approach addresses this request by relying upon meteo/climate service infrastructures (a routine for prediction and DA already in place). A. Carrassi University of Maryland. USA, 24 th April / 33
11 The attribution problem in climate science The DADA approach: justification We want to compute the FAR in the aftermath of an event of interest: F AR = 1 p0 p 1 The likelihoods p 0 /p 1 can be obtained as a by product of the data assimilation. Recognize that the likelihoods of interest are the normalization factors in the Bayesian updates. DA state update p(x y) = p(y x)p(x) p(y) = p(y) = p(y x)p(x) p(x y) They correspond to the marginal likelihood of the data and are referred to as model evidence in statistical literature. Questions: 1 Can we use DA to estimate them and in the aftermath of a climate extreme? 2 Would this be efficient? 3 How would this compare with the conventional approach? A. Carrassi University of Maryland. USA, 24 th April / 33
12 Model evidence: a key quantity Definition and interpretation Model evidence Definition & Interpretation Model evidence p(y K: M) = p(y K: ) = p(y K, y K 1,..., y ) Likelihood of (the infinite from the far past) observation sequence p(y K: ) = dx p(y K: x)p(x) It is the marginal likelihood of the data. p(x) plays the role of the prior. It depends on the underlying dynamics (the model) but it can be formulated as being dependent on any hypothesis under scrutiny. If the dynamics is ergodic and autonomous = p(x) is the invariant distribution on its attractor. p(y K: ) reflects this invariant distribution at a level of accuracy dependent on that of the observation model. We define p(y K: ) as the climatological model evidence. A. Carrassi University of Maryland. USA, 24 th April / 33
13 Model evidence: a key quantity Definition and interpretation Model evidence Applications General metric for model selection and comparison (see e.g. Metref et al., 2018). It can be used to compare the skill of several candidate models, model settings, or boundary conditions, in representing a given observed phenomenon (see e.g. Carson et al., 2018). Calibrating the state-space model s parameters (see e.g. Pulido et al., 2018). Bayesian inference (see e.g. Elsheikh et al., 2014). Quantifying the evidence supporting several theoretical hypothesis associated to a phenomena. Discriminate sources (i.e. nuclear release - Winiarek et al., 2011) A. Carrassi University of Maryland. USA, 24 th April / 33
14 Model evidence: a key quantity Definition and interpretation Statistical/climatological model evidence - Properties and issues The climatological model evidence embeds global information about the system of interest. Nevertheless: 1 It is very difficult to compute in large dimensional systems (it involves solving an integral over the state vector). 2 Its climatological character is not informative about the system s specific current conditions. 3 If the system is driven a time-dependent forcing (i.e. anthropogenic changes in greenhouse gases) its dynamics is non-autonomous the existence of an invariant p(x) is not guaranteed. A. Carrassi University of Maryland. USA, 24 th April / 33
15 Model evidence: a key quantity Contextual model evidence From climatological to contextual model evidence Proposal: move from a global to a context-dependent formulation of the likelihood p(y K: ) p(y K:1 y 0: ) and p(y K:1 y 0: ) is defined as the contextual model evidence. y 0: is used to specify the context. The informational content from past observations is propagated forward by conditioning on y 0: y K:1 is used for evidence diagnostic. The interval t K t 0 defines the evidencing window. A. Carrassi University of Maryland. USA, 24 th April / 33
16 Model evidence: a key quantity Contextual model evidence Contextual model evidence It narrows the climate perspective to the context at t 0 p(y K:1 y 0: ) = dx 0 p(y K:1 x 0 )p(x 0 y 0: ) We will show that contextual model evidence p(y K:1 y 0: ) can be estimated using DA To this purpose note that: The conditional pdf p(x 0 y 0: ) plays the role of the prior density and substitutes p(x) of the climatological model evidence. p(x 0 y 0: ) is the natural (albeit approximate) outcome of the DA machinery. The contextual model evidence is the product of the individual contextual model evidences p(y K:1 y 0: ) = K p(y k y k 1: ) k=1 A. Carrassi University of Maryland. USA, 24 th April / 33
17 Model evidence: a key quantity Computing model evidence using data assimilation Computing model evidence using data assimilation How it works Schema'c diagram of a forecast- - assimila'on cycle and contextual model evidence (CME) computa'on K- long evidencing window True/Correct model M 1 Incorrect model M 0 observa8ons EnKF- based Trajectory 1 Suppose DA (e.g. EnKF, 4DVar,...) is routinely used in a meteo/climate center using the state-of-the-art environmental model (factual model). 2 If a relevant climate event occurs then use the DA to compute the factual likelihood p 1. This comes without additional cost; it is automatically computable using the innovations and the posterior PDF obtained by the DA procedure. 3 Repeat the DA process over a specified evidencing window that encompasses the event of interest but using the counterfactual model (e.g. without the CO2 forcing). 4 This allows to compute the counterfactual likelihood p 0. 5 Use p 1 and p 0 to compute the fraction of attributable risk. A. Carrassi University of Maryland. USA, 24 th April / 33
18 Model evidence: a key quantity Computing model evidence using data assimilation How to compute contextual model evidence Method 1. Importance sampling Compute the model evidence over the window [t 0, t K ] via Importance Sampling: 1 sample N members, x i o, from p(x 0 y 0: ) and forecast from t 0 to t K 2 weight each member trajectory x i K:0 according to their likelihood p(y K:1 x i K:0 ), so as to obtain N p(y K:1 y 0: ) p(y K:1 x i 0) 1 N Importance sampling (IS): use an ensemble-da method up to t 0 = the assumption N M (i.e. ensemble size much smaller than model dimension) can be used in steps (1) and (2) Evaluation of model evidence using a particle filter proposed by Reich and Cotter, i=1 A. Carrassi University of Maryland. USA, 24 th April / 33
19 Model evidence: a key quantity Computing model evidence using data assimilation Method 2. Filtering: Linear/Gaussian case Kalman filter In the Gaussian linear case (Kalman filter - KF), the likelihood and forecast pdfs are Gaussian The Kalman filter estimate of the contextual model evidence reads p(y K:1 y 0: ) = ( K exp 1 2 y k H k x f k (2π) d R k + H k P f k HT k k=1 2 R k +H k P f k HT k R k, P f k : observation and forecast error covariances H k: observation operator. Notation remark: x 2 A = xt A 1 x. ) A. Carrassi University of Maryland. USA, 24 th April / 33
20 Model evidence: a key quantity Computing model evidence using data assimilation Method 2bis. Filtering: nonlinear/gaussian case ensemble Kalman filter Using an ensemble Kalman filter (EnKF), the Kalman filter formula for model evidence now becomes ( ) K exp 1 2 y k H k x f 2 k R k +Y k Yk p(y K:1 y 0: ) T k=1 (2π) d R k + Y k Yk T Y k = H k X k are the normalized observation anomalies, and X k the forecast anomalies. The inversion of R k + Y k Y T k can be done in ensemble subspace. The estimation is exact when the models are linear, the initial condition and observation errors are Gaussian and when the ensemble perturbation span the full range of uncertainty. A. Carrassi University of Maryland. USA, 24 th April / 33
21 Model evidence: a key quantity Computing model evidence using data assimilation Method 3. Smoothing: nonlinear/gaussian case Use a smoothing DA method to estimate p(y K:1 y 0: ) Compute the evidence with a saddle-point approximation (Laplace method): p(y K:1 y 0: ) = dx 0 p(y K:1 x 0 )p(x 0 y 0: ) ( ) dx 0 exp 1 K y k H k M k:0 (x 0 ) 2 R 2 k + ln p(x 0 y 0: ) k=1 = dx 0 exp ( L(y K:1, x 0 )) exp ( L(y K:1, x 0)) 2 x0 L(x 0 ) How to obtain x 0 characterizes the problem: 4DVar, En-smoother, etc. We shall use En4DVar and IEnKS. A. Carrassi University of Maryland. USA, 24 th April / 33
22 Model evidence: a key quantity Computing model evidence using data assimilation 3. Smoothing En4DVar and IEnKS Ensemble-4DVar (En4DVar) (see nomenclature by Lorenc, 2013) The state is parameterized in the ensemble space: x 0 = x 0 + X 0 w ( exp 1 K 2 k=1 y k H k M k:0 (x 0) 2 R k 1 2 w 2) p(y K:1 y 0: ) (2π) Kd K k=1 R k I N + K k=1 (Y k )T R 1 k Y k (Quasi-static) Iterative Ensemble Kalman Smoother (IEnKS) The initial condition and the ensemble at t 0 are sequentially updated, x 0 = x k 1 + X k 1 w k ) K exp ( 12 y k H k M k:0 (x k ) 2Rk+Y k(y k) T p(y K:1 y 0: ) k=1 (2π) d Rk + Yk (Y k )T A. Carrassi University of Maryland. USA, 24 th April / 33
23 Model evidence: Numerical results Setup Model evidence Numerical experiments Schema'c diagram of a forecast- - assimila'on cycle and contextual model evidence (CME) computa'on K- long evidencing window EnKF- based Trajectory True/Correct model M 1 Incorrect model M 0 observa8ons Forced L63 Model dx dt = σ(y x) + λ i cos θ Factual λ 1 = 0; Counterfactual λ 0 { 8 : 8} ; θ = 140. dy dt = ρx y xz + λ i sin θ L95 Model dx m = x m 1 (x m+1 x m 2 ) x m + F i m = 1,..., M dt Factual F 1 = 8; Counterfactual F 0 {5 : 11}. dz = xy βz dt A. Carrassi University of Maryland. USA, 24 th April / 33
24 Model evidence: Numerical results Accuracy of the DA-based estimates Computation of the model evidence: Accuracy assessment Comparison with two independent highly precise methods: 1 Gauss-Hermite quadrature GHQ with 32 degrees - for L63 2 Massive Monte Carlo MC with N particles (power-law extrapolation) - for L63 and L95 The most sophisticated DA method (IEnKS) provides the most accurate outcome. The En4DVar fails in a highly nonlinear setup represented here by the bi-modal L63 model. The EnKF provides sufficiently accurate estimate at a much reduced cost Suitable operationally. A. Carrassi University of Maryland. USA, 24 th April / 33
25 Model evidence: Numerical results Accuracy of the DA-based estimates Model evidence versus forcing and length of window log(p 0 ) log(p 0 ) (a) L63 (K=10) λ 0 (c) L63 (λ 0 =8) K log(p 0 ) log(p 0 ) (b) L95 (K=10) IS EnKF En-4D-Var IEnKS F (d) L95 (F 0 =11) K A. Carrassi University of Maryland. USA, 24 th April / 33
26 Model evidence: Numerical results Low dimensional models - Application for attribution Low dimensional models - Application for attribution Forced L63 with stochastic noise : dx dt = σ(y x) + λ i cos θ i + v x dy dt = ρx y xz + λ i sin θ i + v y dz = xy βz + vz dt Factual λ 1 = 0; Counter-Factual λ 0 = 0 : 40 GINI index measures the discriminating power, from 0 for random to 1 for perfect discrimination EnKF-based - CONVENTIONAL Gin Index Gin Index λ forcing σ q Model Error Gin Index Gini Index σ q λ forcing σ r Obs Error A. Carrassi University of Maryland. USA, 24 th April / 33
27 Model evidence: Numerical results Model evidence with localization Model evidence with domain localization for high-dimension model In high-dimensional systems EnKF methods require localization We derived the Local CME at each gridpoint s using the local observations y s, K { p(y K:1 s ) (2π) d 2 Σk 1 2 exp 1 } 2 (y k s H k sx f k) T Σ 1 k (y k s H k sx f k) k=1 with Σ k = H k s P f k HT k s + R k s and d the size of y k s, P f k is the prior error covariance matrix, R k s is the local obs. error cov. matrix, H k s is the local linear obs. operator. ρ loc S Based on the local CME we define the Domain localized CME (DL-CME) as an euristic global estimator { } p(y K:1) = exp w(s) ln{p(y s )}, s Γ with w(s), scalar weights inversely proportional to the localization radius. A. Carrassi University of Maryland. USA, 24 th April / 33
28 Model evidence: Numerical results Application for model selection High(er) dimensional models - Application for model selection A global atmospheric model resolving the large scale dynamic (Molteni, 2003 ) The models Res.: O(10 4 ) Vor, Div, T, Q, log(p s) True trajectory: 5 month SPEEDY (01/02-30/06/1983) Two convective relaxation times Correct param.: τ cnv = 6hs Incorrect param.: τ cnv = 9hs The observations Artificial obs.: [u, v, T, Q, p s] Frequ.: 6h Spat. distrib.: 1/4 gridpoints DA setup LETKF 50 members (Miyoshi, 2005, 2007) Hydrostat., σ-coord, spectral-transf. Convect., condens., clouds, radiat. Proba. of selection: number of successfull selection DA using all obs.; CME computed for seperate var. A. Carrassi University of Maryland. USA, 24 th April / 33
29 Model evidence: Numerical results Application for the attribution problem (high-dim) High(er) dimensional models - Application for attribution (ongoing) Factual and counterfactual models Two 5-year intervals of interest: and SPEEDY configurations: Factual Counterfactual Projet CMIP5 CMIP5 Model MPI-ESM-LR MPI-ESM-LR Institute MPI-M MPI-M Experiment historical picontrol Time freq. monthly monthly T (K) T (K) Averaged T Fact. CFact Averaged T differences Anomalies Observations Modeling realm atmos/sea atmos/sea Ensemble r1i1p1 r1i1p1 Factual perturbed: σ u = σ v = 1 m.s 1, σ T = 1 C, σ q = 1 gkg 1, σ ps = 1 hpa. Frequ.: 6h; Spat. distrib.: random 1/4 x grid; DA setup LETKF 50 members; Domain localization (h: 700 km, v: 4 km) A. Carrassi University of Maryland. USA, 24 th April / 33
30 Model evidence: Numerical results Application for the attribution problem (high-dim) Global study: DL-CME differences: p (M 0, M 1 ) = ln{ p(y K: M 1 )} ln{ p(y K: M 0 )} all obs. Periods p u v T q all Percentage of attributed cycles (K = 1) The FAR coefficient: FAR= 1 p(y K: M0) p(y K: M 1) all obs. FAR Periods FAR u 0.18 v 0.13 T 0.40 q 0.07 all 0.59 Monthly FAR averaged over the 5-year intervals A. Carrassi University of Maryland. USA, 24 th April / 33
31 Model evidence: Numerical results Application for the attribution problem (high-dim) Spatial study JJA DJF JJA DJF The CME allows for a geographical assessment of the attribution. A. Carrassi University of Maryland. USA, 24 th April / 33
32 Conclusion and future directions Conclusion We showed analytic formulae to compute model evidence using DA (Gaussian filter and smoother) and proved their validity with numerical experiments. The accuracy of this DA-based estimate scales with the sophistication of the DA method (IS EnKF IEnKS). Applications: 1 Attribution of climate-related event attribution evaluation in the aftermath of the event occurrence using DA. (Hannart et al., 2016; Carrassi et al., 2017) 2 Model selection (Metref et al., 2018) 3 Parameter estimation (Carrassi et al., 2017; Pulido et al., 2018) Open questions and future directions Model error is a central concern. 1 It can mask the factor whose casual attribution is under scrutiny confuse discrimination between factual and counter-factual worlds. 2 A DA-based method has the potential to deal with this issue taking advantage of the knowledge in model error treatment in DA. Applications with real observations is on-going. A. Carrassi University of Maryland. USA, 24 th April / 33
33 Conclusion and future directions Key references Carrassi A., M. Bocquet, A. Hannart and M. Ghil, 2017: Estimating model evidence using Data assimilation. Q. J. R. Meteorol. Soc., 143, , 2017 Carrassi, A., M. Bocquet, L. Bertino and G. Evensen. Data assimilation in the geosciences - An overview on methods, issues and perspectives. Submitted to WIREs Climate Change. Available at: arxiv: v2., Hannart A., A. Carrassi, M. Bocquet, M. Ghil, P. Naveau, M. Pulido, J. Ruiz and P. Tandeo, DADA: Data Assimilation for the Detection and Attribution of Weather- and Climate-related Events. Climatic Change, 136, , Metref S., A. Hannart, J. Ruiz, M. Bocquet, A. Carrassi, and M. Ghil, 2018: Estimating model evidence using ensemble-based data assimilation with localization - The model selection problem. Q. J. R. Meteorol. Soc. - Submitted Pulido, M., P. Tandeo, M. Bocquet, A. Carrassi and M. Lucini, Stochastic parametrization identification using ensemble Kalman filtering combined with expectation-minimization and Newton-Raphson maximum likelihood methods. In press on Tellus. A. Carrassi University of Maryland. USA, 24 th April / 33
Estimating model evidence using data assimilation
Estimating model evidence using data assimilation Application to the attribution of climate-related events Marc Bocquet 1, Alberto Carrassi 2, Alexis Hannart 3 & Michael Ghil 4,5 and the DADA Team (1 CEREA,
More informationData assimilation in the geosciences An overview
Data assimilation in the geosciences An overview Alberto Carrassi 1, Olivier Talagrand 2, Marc Bocquet 3 (1) NERSC, Bergen, Norway (2) LMD, École Normale Supérieure, IPSL, France (3) CEREA, joint lab École
More informationDetection and attribution of climate change based on data assimilation
Summer School on Dynamics, Stochastics and Predictability of the Climate System Valsavarenche, Valle d'aosta, Italy, 9-18 June, 2014 Detection and attribution of climate change based on data assimilation
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 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 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 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 informationFundamentals of Data Assimila1on
014 GSI Community Tutorial NCAR Foothills Campus, Boulder, CO July 14-16, 014 Fundamentals of Data Assimila1on Milija Zupanski Cooperative Institute for Research in the Atmosphere Colorado State University
More informationData assimilation; comparison of 4D-Var and LETKF smoothers
Data assimilation; comparison of 4D-Var and LETKF smoothers Eugenia Kalnay and many friends University of Maryland CSCAMM DAS13 June 2013 Contents First part: Forecasting the weather - we are really getting
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 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 informationDADA: Data Assimilation for the Detection and Attribution of weather- and climate-related events
DADA: Data Assimilation for the Detection and Attribution of weather- and climate-related events Alexis Hannart, Alberto Carrassi, Marc Bocquet, Michael Ghil, Philippe Naveau, Manuel Pulido, Juan Ruiz,
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 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 informationLecture 2: From Linear Regression to Kalman Filter and Beyond
Lecture 2: From Linear Regression to Kalman Filter and Beyond January 18, 2017 Contents 1 Batch and Recursive Estimation 2 Towards Bayesian Filtering 3 Kalman Filter and Bayesian Filtering and Smoothing
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 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 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 ECODYC10, Dresden 28 January 2010 Relationship
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 information4. DATA ASSIMILATION FUNDAMENTALS
4. DATA ASSIMILATION FUNDAMENTALS... [the atmosphere] "is a chaotic system in which errors introduced into the system can grow with time... As a consequence, data assimilation is a struggle between chaotic
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 informationLecture 6: Bayesian Inference in SDE Models
Lecture 6: Bayesian Inference in SDE Models Bayesian Filtering and Smoothing Point of View Simo Särkkä Aalto University Simo Särkkä (Aalto) Lecture 6: Bayesian Inference in SDEs 1 / 45 Contents 1 SDEs
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 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 informationFundamentals of Data Assimila1on
2015 GSI Community Tutorial NCAR Foothills Campus, Boulder, CO August 11-14, 2015 Fundamentals of Data Assimila1on Milija Zupanski Cooperative Institute for Research in the Atmosphere Colorado State University
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 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 informationConvergence of the Ensemble Kalman Filter in Hilbert Space
Convergence of the Ensemble Kalman Filter in Hilbert Space Jan Mandel Center for Computational Mathematics Department of Mathematical and Statistical Sciences University of Colorado Denver Parts based
More informationOn the convergence of (ensemble) Kalman filters and smoothers onto the unstable subspace
On the convergence of (ensemble) Kalman filters and smoothers onto the unstable subspace Marc Bocquet CEREA, joint lab École des Ponts ParisTech and EdF R&D, Université Paris-Est, France Institut Pierre-Simon
More informationA Global Atmospheric Model. Joe Tribbia NCAR Turbulence Summer School July 2008
A Global Atmospheric Model Joe Tribbia NCAR Turbulence Summer School July 2008 Outline Broad overview of what is in a global climate/weather model of the atmosphere Spectral dynamical core Some results-climate
More informationAn 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 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 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 informationVariational data assimilation
Background and methods NCEO, Dept. of Meteorology, Univ. of Reading 710 March 2018, Univ. of Reading Bayes' Theorem Bayes' Theorem p(x y) = posterior distribution = p(x) p(y x) p(y) prior distribution
More informationExploring and extending the limits of weather predictability? Antje Weisheimer
Exploring and extending the limits of weather predictability? Antje Weisheimer Arnt Eliassen s legacy for NWP ECMWF is an independent intergovernmental organisation supported by 34 states. ECMWF produces
More informationLecture 9. Time series prediction
Lecture 9 Time series prediction Prediction is about function fitting To predict we need to model There are a bewildering number of models for data we look at some of the major approaches in this lecture
More information1. Current atmospheric DA systems 2. Coupling surface/atmospheric DA 3. Trends & ideas
1 Current issues in atmospheric data assimilation and its relationship with surfaces François Bouttier GAME/CNRM Météo-France 2nd workshop on remote sensing and modeling of surface properties, Toulouse,
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 informationIntroduction to Data Assimilation. Saroja Polavarapu Meteorological Service of Canada University of Toronto
Introduction to Data Assimilation Saroja Polavarapu Meteorological Service of Canada University of Toronto GCC Summer School, Banff. May 22-28, 2004 Outline of lectures General idea Numerical weather prediction
More informationNumerical Weather prediction at the European Centre for Medium-Range Weather Forecasts (2)
Numerical Weather prediction at the European Centre for Medium-Range Weather Forecasts (2) Time series curves 500hPa geopotential Correlation coefficent of forecast anomaly N Hemisphere Lat 20.0 to 90.0
More informationPar$cle Filters Part I: Theory. Peter Jan van Leeuwen Data- Assimila$on Research Centre DARC University of Reading
Par$cle Filters Part I: Theory Peter Jan van Leeuwen Data- Assimila$on Research Centre DARC University of Reading Reading July 2013 Why Data Assimila$on Predic$on Model improvement: - Parameter es$ma$on
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 informationAn introduction to data assimilation. Eric Blayo University of Grenoble and INRIA
An introduction to data assimilation Eric Blayo University of Grenoble and INRIA Data assimilation, the science of compromises Context characterizing a (complex) system and/or forecasting its evolution,
More informationBayesian Statistics and Data Assimilation. Jonathan Stroud. Department of Statistics The George Washington University
Bayesian Statistics and Data Assimilation Jonathan Stroud Department of Statistics The George Washington University 1 Outline Motivation Bayesian Statistics Parameter Estimation in Data Assimilation Combined
More informationComparisons between 4DEnVar and 4DVar on the Met Office global model
Comparisons between 4DEnVar and 4DVar on the Met Office global model David Fairbairn University of Surrey/Met Office 26 th June 2013 Joint project by David Fairbairn, Stephen Pring, Andrew Lorenc, Neill
More informationConsider the joint probability, P(x,y), shown as the contours in the figure above. P(x) is given by the integral of P(x,y) over all values of y.
ATMO/OPTI 656b Spring 009 Bayesian Retrievals Note: This follows the discussion in Chapter of Rogers (000) As we have seen, the problem with the nadir viewing emission measurements is they do not contain
More informationThe analog data assimilation: method, applications and implementation
The analog data assimilation: method, applications and implementation Pierre Tandeo pierre.tandeo@imt-atlantique.fr IMT-Atlantique, Brest, France 16 th of February 2018 Before starting Works in collaboration
More informationDer SPP 1167-PQP und die stochastische Wettervorhersage
Der SPP 1167-PQP und die stochastische Wettervorhersage Andreas Hense 9. November 2007 Overview The priority program SPP1167: mission and structure The stochastic weather forecasting Introduction, Probabilities
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 informationLecture 2: From Linear Regression to Kalman Filter and Beyond
Lecture 2: From Linear Regression to Kalman Filter and Beyond Department of Biomedical Engineering and Computational Science Aalto University January 26, 2012 Contents 1 Batch and Recursive Estimation
More informationHandling nonlinearity in Ensemble Kalman Filter: Experiments with the three-variable Lorenz model
Handling nonlinearity in Ensemble Kalman Filter: Experiments with the three-variable Lorenz model Shu-Chih Yang 1*, Eugenia Kalnay, and Brian Hunt 1. Department of Atmospheric Sciences, National Central
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 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 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 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 informationComparing Local Ensemble Transform Kalman Filter with 4D-Var in a Quasi-geostrophic model
Comparing Local Ensemble Transform Kalman Filter with 4D-Var in a Quasi-geostrophic model Shu-Chih Yang 1,2, Eugenia Kalnay 1, Matteo Corazza 3, Alberto Carrassi 4 and Takemasa Miyoshi 5 1 University of
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 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 informationAnalysis Scheme in the Ensemble Kalman Filter
JUNE 1998 BURGERS ET AL. 1719 Analysis Scheme in the Ensemble Kalman Filter GERRIT BURGERS Royal Netherlands Meteorological Institute, De Bilt, the Netherlands PETER JAN VAN LEEUWEN Institute or Marine
More informationSequential Monte Carlo Samplers for Applications in High Dimensions
Sequential Monte Carlo Samplers for Applications in High Dimensions Alexandros Beskos National University of Singapore KAUST, 26th February 2014 Joint work with: Dan Crisan, Ajay Jasra, Nik Kantas, Alex
More informationIntroduction to Particle Filters for Data Assimilation
Introduction to Particle Filters for Data Assimilation Mike Dowd Dept of Mathematics & Statistics (and Dept of Oceanography Dalhousie University, Halifax, Canada STATMOS Summer School in Data Assimila5on,
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 informationStochastic Collocation Methods for Polynomial Chaos: Analysis and Applications
Stochastic Collocation Methods for Polynomial Chaos: Analysis and Applications Dongbin Xiu Department of Mathematics, Purdue University Support: AFOSR FA955-8-1-353 (Computational Math) SF CAREER DMS-64535
More informationJi-Sun Kang. Pr. Eugenia Kalnay (Chair/Advisor) Pr. Ning Zeng (Co-Chair) Pr. Brian Hunt (Dean s representative) Pr. Kayo Ide Pr.
Carbon Cycle Data Assimilation Using a Coupled Atmosphere-Vegetation Model and the LETKF Ji-Sun Kang Committee in charge: Pr. Eugenia Kalnay (Chair/Advisor) Pr. Ning Zeng (Co-Chair) Pr. Brian Hunt (Dean
More informationObservability, a Problem in Data Assimilation
Observability, Data Assimilation with the Extended Kalman Filter 1 Observability, a Problem in Data Assimilation Chris Danforth Department of Applied Mathematics and Scientific Computation, UMD March 10,
More informationLecture 7: Optimal Smoothing
Department of Biomedical Engineering and Computational Science Aalto University March 17, 2011 Contents 1 What is Optimal Smoothing? 2 Bayesian Optimal Smoothing Equations 3 Rauch-Tung-Striebel Smoother
More informationECE295, Data Assimila0on and Inverse Problems, Spring 2015
ECE295, Data Assimila0on and Inverse Problems, Spring 2015 1 April, Intro; Linear discrete Inverse problems (Aster Ch 1 and 2) Slides 8 April, SVD (Aster ch 2 and 3) Slides 15 April, RegularizaFon (ch
More informationThe priority program SPP1167 Quantitative Precipitation Forecast PQP and the stochastic view of weather forecasting
The priority program SPP1167 Quantitative Precipitation Forecast PQP and the stochastic view of weather forecasting Andreas Hense 9. November 2007 Overview The priority program SPP1167: mission and structure
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 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 informationBackground and observation error covariances Data Assimilation & Inverse Problems from Weather Forecasting to Neuroscience
Background and observation error covariances Data Assimilation & Inverse Problems from Weather Forecasting to Neuroscience Sarah Dance School of Mathematical and Physical Sciences, University of Reading
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 informationPart 1: Expectation Propagation
Chalmers Machine Learning Summer School Approximate message passing and biomedicine Part 1: Expectation Propagation Tom Heskes Machine Learning Group, Institute for Computing and Information Sciences Radboud
More informationData Assimilation: Finding the Initial Conditions in Large Dynamical Systems. Eric Kostelich Data Mining Seminar, Feb. 6, 2006
Data Assimilation: Finding the Initial Conditions in Large Dynamical Systems Eric Kostelich Data Mining Seminar, Feb. 6, 2006 kostelich@asu.edu Co-Workers Istvan Szunyogh, Gyorgyi Gyarmati, Ed Ott, Brian
More informationStability of Ensemble Kalman Filters
Stability of Ensemble Kalman Filters Idrissa S. Amour, Zubeda Mussa, Alexander Bibov, Antti Solonen, John Bardsley, Heikki Haario and Tuomo Kauranne Lappeenranta University of Technology University of
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 informationTowards a probabilistic hydrological forecasting and data assimilation system. Henrik Madsen DHI, Denmark
Towards a probabilistic hydrological forecasting and data assimilation system Henrik Madsen DHI, Denmark Outline Hydrological forecasting Data assimilation framework Data assimilation experiments Concluding
More informationData assimilation : Basics and meteorology
Data assimilation : Basics and meteorology Olivier Talagrand Laboratoire de Météorologie Dynamique, École Normale Supérieure, Paris, France Workshop on Coupled Climate-Economics Modelling and Data Analysis
More informationTropical cyclone simulations and predictions with GFDL s prototype global cloud resolving model
Tropical cyclone simulations and predictions with GFDL s prototype global cloud resolving model S.-J. Lin and GFDL model development teams NOAA/Geophysical Fluid Dynamics Laboratory Workshop on High-Resolution
More information6 Sequential Data Assimilation for Nonlinear Dynamics: The Ensemble Kalman Filter
6 Sequential Data Assimilation for Nonlinear Dynamics: The Ensemble Kalman Filter GEIR EVENSEN Nansen Environmental and Remote Sensing Center, Bergen, Norway 6.1 Introduction Sequential data assimilation
More informationCurrent Issues and Challenges in Ensemble Forecasting
Current Issues and Challenges in Ensemble Forecasting Junichi Ishida (JMA) and Carolyn Reynolds (NRL) With contributions from WGNE members 31 th WGNE Pretoria, South Africa, 26 29 April 2016 Recent trends
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 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 informationModel Uncertainty Quantification for Data Assimilation in partially observed Lorenz 96
Model Uncertainty Quantification for Data Assimilation in partially observed Lorenz 96 Sahani Pathiraja, Peter Jan Van Leeuwen Institut für Mathematik Universität Potsdam With thanks: Sebastian Reich,
More informationThe ECMWF Hybrid 4D-Var and Ensemble of Data Assimilations
The Hybrid 4D-Var and Ensemble of Data Assimilations Lars Isaksen, Massimo Bonavita and Elias Holm Data Assimilation Section lars.isaksen@ecmwf.int Acknowledgements to: Mike Fisher and Marta Janiskova
More informationData Assimilation in the Geosciences
Data Assimilation in the Geosciences An overview on methods, issues and perspectives arxiv:1709.02798v3 [physics.ao-ph] 8 Jun 2018 Alberto Carrassi, Marc Bocquet, Laurent Bertino and Geir Evensen Article
More informationProbabilistic Graphical Models
Probabilistic Graphical Models Brown University CSCI 2950-P, Spring 2013 Prof. Erik Sudderth Lecture 12: Gaussian Belief Propagation, State Space Models and Kalman Filters Guest Kalman Filter Lecture by
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 informationEnsemble Verification Metrics
Ensemble Verification Metrics Debbie Hudson (Bureau of Meteorology, Australia) ECMWF Annual Seminar 207 Acknowledgements: Beth Ebert Overview. Introduction 2. Attributes of forecast quality 3. Metrics:
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 informationThe hybrid ETKF- Variational data assimilation scheme in HIRLAM
The hybrid ETKF- Variational data assimilation scheme in HIRLAM (current status, problems and further developments) The Hungarian Meteorological Service, Budapest, 24.01.2011 Nils Gustafsson, Jelena Bojarova
More informationData assimilation as an optimal control problem and applications to UQ
Data assimilation as an optimal control problem and applications to UQ Walter Acevedo, Angwenyi David, Jana de Wiljes & Sebastian Reich Universität Potsdam/ University of Reading IPAM, November 13th 2017
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 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 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 informationImproved analyses and forecasts with AIRS retrievals using the Local Ensemble Transform Kalman Filter
Improved analyses and forecasts with AIRS retrievals using the Local Ensemble Transform Kalman Filter Hong Li, Junjie Liu, and Elana Fertig E. Kalnay I. Szunyogh, E. J. Kostelich Weather and Chaos Group
More informationImplementation and evaluation of a regional data assimilation system based on WRF-LETKF
Implementation and evaluation of a regional data assimilation system based on WRF-LETKF Juan José Ruiz Centro de Investigaciones del Mar y la Atmosfera (CONICET University of Buenos Aires) With many thanks
More informationESTIMATING CORRELATIONS FROM A COASTAL OCEAN MODEL FOR LOCALIZING AN ENSEMBLE TRANSFORM KALMAN FILTER
ESTIMATING CORRELATIONS FROM A COASTAL OCEAN MODEL FOR LOCALIZING AN ENSEMBLE TRANSFORM KALMAN FILTER Jonathan Poterjoy National Weather Center Research Experiences for Undergraduates, Norman, Oklahoma
More informationMCMC Sampling for Bayesian Inference using L1-type Priors
MÜNSTER MCMC Sampling for Bayesian Inference using L1-type Priors (what I do whenever the ill-posedness of EEG/MEG is just not frustrating enough!) AG Imaging Seminar Felix Lucka 26.06.2012 , MÜNSTER Sampling
More informationEfficient Data Assimilation for Spatiotemporal Chaos: a Local Ensemble Transform Kalman Filter
arxiv:physics/0511236v2 [physics.data-an] 29 Dec 2006 Efficient Data Assimilation for Spatiotemporal Chaos: a Local Ensemble Transform Kalman Filter Brian R. Hunt Institute for Physical Science and Technology
More information