Data assimilation : Basics and meteorology
|
|
- Clyde Brown
- 6 years ago
- Views:
Transcription
1 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 Paris 23 November
2 What is exactly Data Assimilation? The numerical algorithms. (Ensemble) Kalman Filter, Variational Assimilation, Particle Filters Assimilation vs Stability and Instability. 2
3
4 ECMWF, Technical Report 499, 2006
5 5
6 6
7 December 2007: Satellite data volumes used: around 18 millions per day quantity of satellite data used per day at ECMWF number of data used per day (millions) CONV+SAT WINDS TOTAL Year Value as of March 2010 : 25 millions per day 7
8 S. Louvel, Doctoral Dissertation, 1999
9 Physical laws governing the flow Conservation of mass Dρ/Dt + ρ divu = 0 Conservation of energy De/Dt - (p/ρ 2 ) Dρ/Dt = Q Conservation of momentum DU/Dt + (1/ρ) gradp - g + 2 Ω U = F Equation of state f(p, ρ, e) = 0 (p/ρ = rt, e = C v T) Conservation of mass of secondary components (water in the atmosphere, salt in the ocean, chemical species, ) Dq/Dt + q divu = S Physical laws available in practice in the form of a discretized (and necessarily imperfect) numerical model 9
10 European Centre for Medium-range Weather Forecasts (ECMWF, Reading, UK) Horizontal spherical harmonics triangular truncation T1279 (horizontal resolution 16 kilometres) 91 levels on the vertical (0-80 km) Dimension of state vector n Timestep = 10 minutes 10
11 Purpose of assimilation : reconstruct as accurately as possible the state of the atmospheric or oceanic flow, using all available appropriate information. The latter essentially consists of The observations proper, which vary in nature, resolution and accuracy, and are distributed more or less regularly in space and time. The physical laws governing the evolution of the flow, available in practice in the form of a discretized, and necessarily approximate, numerical model. Asymptotic properties of the flow, such as, e. g., geostrophic balance of middle latitudes. Although they basically are necessary consequences of the physical laws which govern the flow, these properties can usefully be explicitly introduced in the assimilation process. 11
12 Both observations and model are affected with some uncertainty uncertainty on the estimate. For some reason, uncertainty is conveniently described by probability distributions (don t know too well why, but it works). Assimilation is a problem in bayesian estimation. Determine the conditional probability distribution for the state of the system, knowing everything we know 12
13 Difficulties specific to assimilation of meteorological and oceanographical observations : - Very large numerical dimensions (n parameters to be estimated, p observations per 24-hour period). Difficulty aggravated in Numerical Weather Prediction by the need for the forecast to be ready in time. - Non-trivial underlying dynamics. 13
14 Bayesian estimation is however impossible in its general theoretical form in meteorological or oceanographical practice because It is impossible to explicitly describe a probability distribution in a space with dimension even as low as n 10 3, not to speak of the dimension n of present Numerical Weather Prediction models. Probability distribution of errors on data very poorly known (model errors in particular).
15 One has to restrict oneself to a much more modest goal. Two approaches exist at present Obtain some central estimate of the conditional probability distribution (expectation, mode, ), plus some estimate of the corresponding spread (standard deviations and a number of correlations). Produce an ensemble of estimates which are meant to sample the conditional probability distribution (dimension N O(10-100)).
16 Data of the form z = Γx + ζ, ζ N [µ, S] Known data vector z belongs to data space D, dimd = m, Unknown state vector x belongs to state space S, dims = n Γ known (mxn)-matrix, ζ unknown error Then conditional probability distribution is where P(x z) = N [x a, P a ] x a = (Γ T S -1 Γ) -1 Γ T S -1 [z µ] P a = (Γ T S -1 Γ) -1 Determinacy condition : rankγ = n. Requires m n. 16
17 Expressions x a = (Γ T S -1 Γ) -1 Γ T S -1 [z µ] P a = (Γ T S -1 Γ) -1 (2a) (2b) can also be obtained, from the knowledge of only µ and S, as defining the variance minimizing linear estimator of x from z, or Best Linear Unbiased Estimate (BLUE) of x from z (but significance of P a is now totally different) Expressions (2) are commonly used, and are useful, in many (especially, but non only, geophysical) applications where operator Γ is mildly nonlinear (and error ζ can be non-gaussian). In fact, going beyond (2) turns out to be very challenging in many practical problems. 17
18 Variational form. BLUE x a minimizes following scalar objective function, defined on state space S ξ S J(ξ) (1/2) [Γξ - (z-µ)] T S -1 [Γξ - (z-µ)] BLUE is invariant in any invertible linear change of coordinates, either in state or data space. 18
19 If data of the form (not restrictive) x b = x + ζ b, ζ b N [0, P b ] (3a) y = Hx + ε, ε N [0, R] (3b) ζ b and ε independent Eq. 2 become Variational form x a = x b + P b H T [HP b H T + R] -1 (y - Hx b ) P a = P b - P b H T [HP b H T + R] -1 HP b d y - Hx b is innovation vector ξ S J(ξ) = (1/2) (x b - ξ) T [P b ] -1 (x b - ξ) + (1/2) (y - Hξ) T R -1 (y - Hξ) = J b + J o 19
20 Time dimension Observation vector at time k y k = H k x k + ε k k = 0,, K E(ε k ) = 0 ; E(ε k ε jt ) R k δ kj Evolution equation x k+1 = M k x k + η k k = 0,, K-1 E(η k ) = 0 ; E(η k η jt ) Q k δ kj E(η k ε jt ) = 0 Background estimate at time 0 x b 0 = x 0 + ζb 0 E(ζ b 0 ) = 0 ; E(ζb 0 ζ b 0 T ) P b 0 E(ζ b 0ε kt ) = 0 ; E(ζ b 0η kt ) = 0 Errors uncorrelated in time 20
21 BLUE can be obtained through time-sequential algorithm called Kalman filter Background x b k and associated error covariance matrix Pb k known Analysis step x a k = xb k + Pb k H k T [H k P b k H k T + R k ] -1 (y k - H k x b k ) P a k = Pb k - Pb k H k T [H k P b k H k T + R k ] -1 H k P b k Forecast step x b k+1 = M k xa k P b k+1 = M k Pa k M k T + Q k 21
22 Costliest part of computation P b k+1 = M k P a k M k T + Q k Need for determining the temporal evolution of the uncertainty on the state of the system is the major difficulty in assimilation of meteorological and oceanographical observations 22
23 Analysis of 500-hPa geopotential for 1 December 1989, 00:00 UTC (ECMWF, spectral truncation T21, unit m. After F. Bouttier) 23
24 Temporal evolution of the 500-hPa geopotential autocorrelation with respect to point located at 45N, 35W. From top to bottom: initial time, 6- and 24-hour range. Contour interval 0.1. After F. Bouttier. 24
25 Most common approach at present is Ensemble Kalman Filter (EnKF) Uncertainty is represented, not by a covariance matrix, but by an ensemble of point estimates in state space which are meant to sample the conditional probability distribution for the state of the system (dimension N O(10-100)). Ensemble is evolved in time through the full model, which eliminates any need for linear hypothesis as to the temporal evolution. 25
26 How to update predicted ensemble with new observations? Predicted ensemble at time t : {x b i }, i = 1,, N Observation vector at same time : y = Hx + ε Gaussian approach Produce sample of probability distribution for real observed quantity Hx y i = y - ε i where ε i is distributed according to probability distribution for observation error ε. Then use Kalman formula to produce sample of analysed states x a i = xb i + Pb H T [HP b H T + R] -1 (y i - Hx b i ), i = 1,, N (2) where P b is covariance matrix of predicted ensemble {x b i }. Remark. If all errors were Gaussian, and if P b was exact covariance matrix of background error, (2) would achieve Bayesian estimation, in the sense that {x a i } would be a sample of conditional probability distribution for x, given all data up to time t. 26
27 Month-long Performance of EnKF vs. 3Dvar with WRF EnKF 3DVar (prior, solid; posterior, dotted) posterior Prior Better performance of EnKF than 3DVar also seen in both 12-h forecast and posterior analysis in terms of root-mean square difference averaged over the entire month 27 (Meng and Zhang 2007c, MWR, in review )
28 Ensemble Kalman Filter (EnKF), of which many variants have been defined and are actively studied, has now become one of the two major powerful algorithm for assimilation of meteorological and oceanographical observations. Ensemble Transform Kalman Filter (ETKF, Bishop and others) 28
29 Ensemble Kalman Filter is only exception so far, in real life assimilation, to linear (and gaussian) approach. Analysis step remains linear (and gaussian). Bayesian character of EnKF? - If everything is linear, distribution of ensemble produced by EnKF tends to underlying bayesian pdf when ensemble dimension tends to infinity (if errors are not gaussian, gaussianity is ensured by Central Limit Theorem). - If non-linearity, distribution of ensemble tends to a limit which is not the bayesian distribution (F. Le Gland) 29
30 Variational Assimilation Observation vector at time k y k = H k x k + ε k k = 0,, K E(ε k ) = 0 ; E(ε k ε jt ) R k δ kj Evolution equation x k+1 = M k x k + η k k = 0,, K-1 E(η k ) = 0 ; E(η k η jt ) Q k δ kj E(η k ε jt ) = 0 Background estimate at time 0 x b 0 = x 0 + ζb 0 E(ζ b 0 ) = 0 ; E(ζb 0 ζ b 0 T ) P b 0 E(ζ b 0ε kt ) = 0 ; E(ζ b 0η kt ) = 0 Errors uncorrelated in time 30
31 Minimise objective function (ξ 0, ξ 1,..., ξ K ) J(ξ 0, ξ 1,..., ξ K ) = (1/2) (x 0 b - ξ 0 ) T [P 0b ] -1 (x 0 b - ξ 0 ) + (1/2) Σ k=0,,k [y k - H k ξ k ] T R k -1 [y k - H k ξ k ] + (1/2) Σ k=0,,k-1 [ξ k+1 - M k ξ k ] T Q k -1 [ξ k+1 - M k ξ k ] Can include nonlinear M k and/or H k. 4D-Var. In the linear case, equivalent to Kalman Filter. 31
32 4D-Var used operationally in a number of major meteorological centres (European Centre for Medium-Range Weather Forecasts, Météo-France, UK Meteorological Office, Canadian Meteorological Centre, Japan Meteorological Agency, ). In most cases, model error is ignored. Strong constraint variational assimilation ξ 0 J(ξ 0 ) = (1/2) (x b 0 - ξ 0 )T [P b 0 ]-1 (x b 0 - ξ 0 ) + (1/2) Σ k=0,, K (y k - H k ξ k ) T R k -1 (y k - H k ξ k ) subject to ξ k+1 = M k ξ k, k = 0,, K-1 32
33 J(ξ 0 ) = (1/2) (x 0 b - ξ 0 ) T [P 0b ] -1 (x 0 b - ξ 0 ) + (1/2) Σ k [y k - H k ξ k ] T R k -1 [y k - H k ξ k ] Background is not necessary, if observations are in sufficient number to overdetermine the problem. Nor is strict linearity. How to minimize objective function with respect to initial state u = ξ 0 (u is called the control variable of the problem)? Use iterative minimization algorithm, each step of which requires the explicit knowledge of the local gradient u J ( J/ u i ) of J with respect to u. Gradient computed by adjoint method. 33
34 Courtier and Talagrand, QJRMS,
35 Courtier and Talagrand, QJRMS,
36 500-hPa geopotential field as determined by : (left) operational assimilation system of French Weather Service (3D, primitive equation) and (right) experimental variational system (2D, vorticity equation) Courtier and Talagrand, QJRMS,
37 Thépaut et al., 1993, Mon. Wea. Rev., 121,
38 Analysis increments in a 3D-Var corresponding to a height observation at the 250- hpa pressure level (no temporal evolution of background error covariance matrix) 38 Thépaut et al., 1993, Mon. Wea. Rev., 121,
39 Same as before, but at the end of a 24-hr 4D-Var Thépaut et al., 1993, Mon. Wea. Rev., 121,
40 Analysis increments in a 3D-Var corresponding to a u-component wind observation at the 1000-hPa pressure level (no temporal evolution of background error covariance matrix) Thépaut et al., 1993, Mon. Wea. Rev., 121,
41 Same as before, but at the end of a 24-hr 4D-Var Thépaut et al., 1993, Mon. Wea. Rev., 121,
42 ECMWF, Results on one FASTEX case (1997)
43 Buehner (2008) For the same numerical cost, and in meteorologically realistic situations, Ensemble Kalman Filter and Variational Assimilation produce results of similar quality. 43
44 Exact bayesian estimation? Particle filters Predicted ensemble at time t : {x b n, n = 1,, N }, each element with its own weight (probability) P(x b n) Observation vector at same time : y = Hx + ε Bayes formula P(x b n y) P(y xb n ) P(xb n ) Defines updating of weights 44
45 Bayes formula P(x b n y) P(y x b n) P(x b n) Defines updating of weights; particles are not modified. Asymptotically converges to bayesian pdf. Very easy to implement. Observed fact. For large state dimension, ensemble tends to collapse. 45
46 C. Snyder, Snyder.pdf 46
47 Problem originates in the curse of dimensionality Large dimension pdf s are very diffuse, so that very few particles (if any) are present in areas where conditional probability ( likelihood ) P(y x) is large. Bengtsson et al. (2008) and Snyder et al. (2008) evaluate that stability of filter requires the size of ensembles to increase exponentially with space dimension. 47
48 Alternative possibilities (review in van Leeuwen, 2009, Mon. Wea. Rev., ) Resampling. Define new ensemble. Simplest way. Draw new ensemble according to probability distribution defined by the updated weights. Give same weight to all particles. Particles are not modified, but particles with low weights are likely to be eliminated, while particles with large weights are likely to be drawn repeatedly. For multiple particles, add noise, either from the start, or in the form of model noise in ensuing temporal integration. Random character of the sampling introduces noise. Alternatives exist, such as residual sampling (Lui and Chen, 1998, van Leeuwen, 2003). Updated weights w n are multiplied by ensemble dimension N. Then p copies of each particle n are taken, where p is the integer part of Nw n. Remaining particles, if needed, are taken randomly from the resulting distribution. 48
49 Importance Sampling. Use a proposal density that is closer to the new observations than the density defined by the predicted particles (for instance the density defined by EnKF, after the latter has used the new observations). Independence between observations is then lost in the computation of likelihood P(y x). In particular, Guided Sequential Importance Sampling (van Leeuwen, 2002). Idea : use observations performed at time k to resample ensemble at some timestep anterior to k. Implicit Particle Filter. The evolution of each particle is guided by the observations. 49
50 van Leeuwen, 2003, Mon. Wea. Rev., 131,
51 Particle filters are actively studied (Snyder et al., van Leeuwen et al.). They are independent of any hypothesis of linearity or gaussianity. It is not clear however if their cost can be sufficiently reduced to make them usable in the large scale applications of geophysical flows. 51
52 Proportion of resources devoted to assimilation in Numerical Weather Prediction has steadily increased over time. At present at ECMWF, the cost of 24 hours of assimilation is half the global cost of the 10-day forecast (i. e., including the ensemble forecast). 52
53 What is exactly Data Assimilation? The numerical algorithms. (Ensemble) Kalman Filter, Variational Assimilation, Particle Filters Assimilation vs Stability and Instability 53
54 If there is uncertainty on the state of the system, and dynamics of the system is perfectly known, uncertainty on the state along stable modes decreases over time, while uncertainty along unstable modes increases. Stable (unstable) modes : perturbations to the basic state that decrease (increase) over time. 54
55 Observations between times -τ and 0, no background 55
56 Consequence : 4D-Var assimilation, which carries information both forward and backward in time, performed over time interval [t 0, t 1 ] over uniformly distributed noisy data. If assimilating model is perfect, estimation error is concentrated in stable modes at time t 0, and in unstable modes at time t 1. Error is smallest somewhere within interval [t 0, t 1 ]. Similar result holds true for Kalman filter (or more generally any form of sequential assimilation), in which estimation error is concentrated in unstable modes at any time. 56
57 Trevisan et al., 2010, Q. J. R. Meteorol. Soc. 57
58 Lorenz (1963) dx/dt = σ(y-x) dy/dt = ρx - y - xz dz/dt = -βz + xy with parameter values σ = 10, ρ = 28, β = 8/3 chaos 58
59
60 60
61 Pires et al., Tellus, 1996 ; Lorenz system (1963) 61
62 Pires et al., Tellus,
63 Since, after an assimilation has been performed over a period of time, uncertainty is likely to be concentrated in modes that have been unstable, it might be useful for the next assimilation, and at least in terms of cost efficiency, to concentrate corrections on the background in those modes. Actually, presence of residual noise in stable modes can be damageable for analysis and subsequent forecast. Assimilation in the Unstable Subspace (AUS) (Carrassi et al., 2007, 2008, for the case of 3D-Var) 63
64 Four-dimensional variational assimilation in the unstable subspace (4DVar-AUS) Trevisan et al., 2010, Four-dimensional variational assimilation in the unstable subspace and the optimal subspace dimension, Q. J. R. Meteorol. Soc., 136,
65 4D-Var-AUS Algorithmic implementation Define N perturbations to the current state, and evolve them according to the tangent linear model, with periodic reorthonormalization in order to avoid collapse onto the dominant Lyapunov vector (same algorithm as for computation of Lyapunov exponents). Cycle successive 4D-Var s, restricting at each cycle the modification to be made on the current state to the space spanned by the N perturbations emanating from the previous cycle (if N is the dimension of state space, that is identical with standard 4D-Var). 65
66 Experiments performed on the Lorenz (1996) model with value F = 8, which gives rise to chaos. Three values of I have been used, namely I = 40, 60, 80, which correspond to respectively N + = 13, 19 and 26 positive Lyapunov exponents. In all three cases, the largest Lyapunov exponent corresponds to a doubling time of about 2 days (with 1 day = 1/5 model time unit). Identical twin experiments (perfect model) 66
67 Observing system defined as in Fertig et al. (Tellus, 2007): At each observation time, one observation every four grid points (observation points shifted by one grid point at each observation time). Observation frequency : 1.5 hour Random gaussian observation errors with expectation 0 and standard deviation σ 0 = 0.2 ( climatological standard deviation 5.1). Sequences of variational assimilations have been cycled over windows with length τ = 1,, 5 days. Results are averaged over 5000 successive windows. 67
68 No explicit background term (i. e., with error covariance matrix) in objective function : information from past lies in the background to be updated, and in the N perturbations which define the subspace in which updating is to be made. Best performance for N slightly above number N + of positive Lyapunov exponents. 68
69 Different curves are almost identical on all three panels. Relative improvement obtained by decreasing subspace dimension N to its optimal value is largest for smaller window length τ.
70 70
71 Experiments have been performed in which an explicit background term was present, the associated error covariance matrix having been obtained as the average of a sequence of full 4D-Var s. The estimates are systematically improved, and more for full 4D-Var than for 4D-Var- AUS. But they remain qualitatively similar, with best performance for 4D-Var-AUS with N slightly above N +. 71
72 Minimum of objective function cannot be made smaller by reducing control space. Numerical tests show that minimum of objective function is smaller (by a few percent) for full 4D-Var than for 4D-Var-AUS. Full 4D-Var is closer to the noisy observations, but farther away from the truth. And tests also show that full 4D-Var performs best when observations are perfect (no noise). Results show that, if all degrees of freedom that are available to the model are used, the minimization process introduces components along the stable modes of the system, in which no error is present, in order to ensure a closer fit to the observations. This degrades the closeness of the fit to reality. The optimal choice is to restrict the assimilation to the unstable modes. 72
73 Conclusions on AUS Error concentrates in unstable modes at the end of assimilation window. It must therefore be sufficient, at the beginning of new assimilation cycle, to introduce increments only in the subspace spanned by those unstable modes. In the perfect model case (and with noisy observations), assimilation is most efficient when increments are introduced in a space with dimension slightly above the number of non-negative Lyapunov exponents. In the case of imperfect model (and of strong constraint assimilation), preliminary results lead to similar conclusions, with larger optimal subspace dimension, and less well marked optimality (not shown here). Further work necessary. In agreement with theoretical and experimental results obtained for Kalman Filter assimilation (Trevisan and Palatella, Nonlin. Processes Geophys., 2011, Ng et al., Tellus, 2011) 73
74 Lorenz (1963) dx/dt = σ(y-x) dy/dt = ρx - y - xz dz/dt = -βz + xy with parameter values σ = 10, ρ = 28, β = 8/3 chaos 74
75 Pires et al., Tellus, 1996 ; Lorenz system (1963) 75
76 Minima in the variations of objective function correspond to solutions that have bifurcated from the observed solution, and to different folds in state space. 76
77 Quasi-Static Variational Assimilation (QSVA). Increase progressively length of the assimilation window, starting each new assimilation from the result of the previous one. This should ensure, at least if observations are in a sense sufficiently dense in time, that current estimation of the system always lies in the attractive basin of the absolute minimum of objective function. 77
78 Pires et al., Tellus, 1996 ; Lorenz system (1963)
79 79
80 Swanson, Vautard and Pires, 1998, Tellus, 50A,
81 Conclusions on QSVA Quasi-Static Variational Assimilation, which consists in progressively increasing the length of the assimilation window, is capable of efficiently assimilating observations over time windows that are much longer than the length of time over which the tangent linear approximation is valid and the 4DVar objective function is quadratic. For a quasi-geostrophic model, assimilation can be usefully assimilated over periods as large as 10 days. That however requires the observations to be in a sense sufficiently dense in time. Do EnKF, EnKS, or particle filters, achieve the same accuracy? Effect of model errors? 81
82 Assimilation, which originated from the need of defining initial conditions for numerical weather forecasts, has progressively extended to many diverse applications Oceanography Atmospheric chemistry (both troposphere and stratosphere) Oceanic biogeochemistry Ground hydrology Terrestrial biosphere and vegetation cover Glaciology Magnetism (both planetary and stellar) Plate tectonics Planetary atmospheres (Mars, ) Reassimilation of past observations (mostly for climatological purposes, ECMWF, NCEP/NCAR) Identification of source of tracers Parameter identification A priori evaluation of anticipated new instruments Definition of observing systems (Observing Systems Simulation Experiments) Validation of models Sensitivity studies (adjoints) It has now become a major tool of numerical environmental science 82
83 Aubert and Fournier, Nonlin. Processes Geophys.,
84 A few of the (many) remaining problems Assimilation of images (much work being done ) Observability (what does one obtain by observing what? Data are noisy, dynamics is chaotic!) More accurate identification and quantification of errors on data, especially on assimilating model (will always require independent hypotheses) 84
85
Assimilation 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 informationEnsemble variational assimilation as a probabilistic estimator Part 1: The linear and weak non-linear case
Nonlin. Processes Geophys., 25, 565 587, 28 https://doi.org/.594/npg-25-565-28 Author(s) 28. This work is distributed under the Creative Commons Attribution 4. License. Ensemble variational assimilation
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 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 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 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 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 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 informationIntroduction to ensemble forecasting. Eric J. Kostelich
Introduction to ensemble forecasting Eric J. Kostelich SCHOOL OF MATHEMATICS AND STATISTICS MSRI Climate Change Summer School July 21, 2008 Co-workers: Istvan Szunyogh, Brian Hunt, Edward Ott, Eugenia
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 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 informationBackground Error Covariance Modelling
Background Error Covariance Modelling Mike Fisher Slide 1 Outline Diagnosing the Statistics of Background Error using Ensembles of Analyses Modelling the Statistics in Spectral Space - Relaxing constraints
More informationThe Local Ensemble Transform Kalman Filter (LETKF) Eric Kostelich. Main topics
The Local Ensemble Transform Kalman Filter (LETKF) Eric Kostelich Arizona State University Co-workers: Istvan Szunyogh, Brian Hunt, Ed Ott, Eugenia Kalnay, Jim Yorke, and many others http://www.weatherchaos.umd.edu
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 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 informationAccelerating the spin-up of Ensemble Kalman Filtering
Accelerating the spin-up of Ensemble Kalman Filtering Eugenia Kalnay * and Shu-Chih Yang University of Maryland Abstract A scheme is proposed to improve the performance of the ensemble-based Kalman Filters
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 informationIntroduction to Data Assimilation. Reima Eresmaa Finnish Meteorological Institute
Introduction to Data Assimilation Reima Eresmaa Finnish Meteorological Institute 15 June 2006 Outline 1) The purpose of data assimilation 2) The inputs for data assimilation 3) Analysis methods Theoretical
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 informationM.Sc. in Meteorology. Numerical Weather Prediction
M.Sc. in Meteorology UCD Numerical Weather Prediction Prof Peter Lynch Meteorology & Climate Cehtre School of Mathematical Sciences University College Dublin Second Semester, 2005 2006. Text for the Course
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 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 informationIn the derivation of Optimal Interpolation, we found the optimal weight matrix W that minimizes the total analysis error variance.
hree-dimensional variational assimilation (3D-Var) In the derivation of Optimal Interpolation, we found the optimal weight matrix W that minimizes the total analysis error variance. Lorenc (1986) showed
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 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 informationModel-error estimation in 4D-Var
QUARTERLY JOURNAL OF THE ROYAL METEOROLOGICAL SOCIETY Q. J. R. Meteorol. Soc. 33: 67 80 (007) Published online 3 July 007 in Wiley InterScience (www.interscience.wiley.com) DOI: 0.00/qj.94 Model-error
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 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 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 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 informationInter-comparison of 4D-Var and EnKF systems for operational deterministic NWP
Inter-comparison of 4D-Var and EnKF systems for operational deterministic NWP Project eam: Mark Buehner Cecilien Charette Bin He Peter Houtekamer Herschel Mitchell WWRP/HORPEX Workshop on 4D-VAR and Ensemble
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 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 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 informationSome Applications of WRF/DART
Some Applications of WRF/DART Chris Snyder, National Center for Atmospheric Research Mesoscale and Microscale Meteorology Division (MMM), and Institue for Mathematics Applied to Geoscience (IMAGe) WRF/DART
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 informationAssimilation Techniques (4): 4dVar April 2001
Assimilation echniques (4): 4dVar April By Mike Fisher European Centre for Medium-Range Weather Forecasts. able of contents. Introduction. Comparison between the ECMWF 3dVar and 4dVar systems 3. he current
More informationThe Structure of Background-error Covariance in a Four-dimensional Variational Data Assimilation System: Single-point Experiment
ADVANCES IN ATMOSPHERIC SCIENCES, VOL. 27, NO. 6, 2010, 1303 1310 The Structure of Background-error Covariance in a Four-dimensional Variational Data Assimilation System: Single-point Experiment LIU Juanjuan
More informationSensitivity analysis in variational data assimilation and applications
Sensitivity analysis in variational data assimilation and applications Dacian N. Daescu Portland State University, P.O. Box 751, Portland, Oregon 977-751, U.S.A. daescu@pdx.edu ABSTRACT Mathematical aspects
More informationR. E. Petrie and R. N. Bannister. Department of Meteorology, Earley Gate, University of Reading, Reading, RG6 6BB, United Kingdom
A method for merging flow-dependent forecast error statistics from an ensemble with static statistics for use in high resolution variational data assimilation R. E. Petrie and R. N. Bannister Department
More informationDATA ASSIMILATION FOR FLOOD FORECASTING
DATA ASSIMILATION FOR FLOOD FORECASTING Arnold Heemin Delft University of Technology 09/16/14 1 Data assimilation is the incorporation of measurement into a numerical model to improve the model results
More informationIntroduction to initialization of NWP models
Introduction to initialization of NWP models weather forecasting an initial value problem traditionally, initialization comprised objective analysis of obs at a fixed synoptic time, i.e. 00Z or 12Z: data
More informationRosemary Munro*, Graeme Kelly, Michael Rohn* and Roger Saunders
ASSIMILATION OF METEOSAT RADIANCE DATA WITHIN THE 4DVAR SYSTEM AT ECMWF Rosemary Munro*, Graeme Kelly, Michael Rohn* and Roger Saunders European Centre for Medium Range Weather Forecasts Shinfield Park,
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 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 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 informationInterpreting Adjoint and Ensemble Sensitivity toward the Development of Optimal Observation Targeting Strategies
Interpreting Adjoint and Ensemble Sensitivity toward the Development of Optimal Observation Targeting Strategies Brian C. Ancell 1, and Gregory J Hakim University of Washington, Seattle, WA Submitted to
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 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 informationLecture notes on assimilation algorithms
Lecture notes on assimilation algorithms Elías alur Hólm European Centre for Medium-Range Weather Forecasts Reading, UK April 18, 28 1 Basic concepts 1.1 The analysis In meteorology and other branches
More informationMultiple-scale error growth and data assimilation in convection-resolving models
Numerical Modelling, Predictability and Data Assimilation in Weather, Ocean and Climate A Symposium honouring the legacy of Anna Trevisan Bologna, Italy 17-20 October, 2017 Multiple-scale error growth
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 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 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 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 informationSimple Doppler Wind Lidar adaptive observation experiments with 3D-Var. and an ensemble Kalman filter in a global primitive equations model
1 2 3 4 Simple Doppler Wind Lidar adaptive observation experiments with 3D-Var and an ensemble Kalman filter in a global primitive equations model 5 6 7 8 9 10 11 12 Junjie Liu and Eugenia Kalnay Dept.
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 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 information6.5 Operational ensemble forecasting methods
6.5 Operational ensemble forecasting methods Ensemble forecasting methods differ mostly by the way the initial perturbations are generated, and can be classified into essentially two classes. In the first
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 informationOn the Kalman Filter error covariance collapse into the unstable subspace
Nonlin. Processes Geophys., 18, 243 250, 2011 doi:10.5194/npg-18-243-2011 Author(s) 2011. CC Attribution 3.0 License. Nonlinear Processes in Geophysics On the Kalman Filter error covariance collapse into
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 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 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 informationBred Vectors, Singular Vectors, and Lyapunov Vectors in Simple and Complex Models
Bred Vectors, Singular Vectors, and Lyapunov Vectors in Simple and Complex Models Adrienne Norwood Advisor: Eugenia Kalnay With special thanks to Drs. Kayo Ide, Brian Hunt, Shu-Chih Yang, and Christopher
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 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 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 informationWeak Constraints 4D-Var
Weak Constraints 4D-Var Yannick Trémolet ECMWF Training Course - Data Assimilation May 1, 2012 Yannick Trémolet Weak Constraints 4D-Var May 1, 2012 1 / 30 Outline 1 Introduction 2 The Maximum Likelihood
More informationParticle filters, the optimal proposal and high-dimensional systems
Particle filters, the optimal proposal and high-dimensional systems Chris Snyder National Center for Atmospheric Research Boulder, Colorado 837, United States chriss@ucar.edu 1 Introduction Particle filters
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 informationThe Impact of Background Error on Incomplete Observations for 4D-Var Data Assimilation with the FSU GSM
The Impact of Background Error on Incomplete Observations for 4D-Var Data Assimilation with the FSU GSM I. Michael Navon 1, Dacian N. Daescu 2, and Zhuo Liu 1 1 School of Computational Science and Information
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 informationVariational Data Assimilation Current Status
Variational Data Assimilation Current Status Eĺıas Valur Hólm with contributions from Mike Fisher and Yannick Trémolet ECMWF NORDITA Solar and stellar dynamos and cycles October 2009 Eĺıas Valur Hólm (ECMWF)
More informationDiagnosis of observation, background and analysis-error statistics in observation space
Q. J. R. Meteorol. Soc. (2005), 131, pp. 3385 3396 doi: 10.1256/qj.05.108 Diagnosis of observation, background and analysis-error statistics in observation space By G. DESROZIERS, L. BERRE, B. CHAPNIK
More informationEnsembles and Particle Filters for Ocean Data Assimilation
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Ensembles and Particle Filters for Ocean Data Assimilation Robert N. Miller College of Oceanic and Atmospheric Sciences
More informationIntroduction to Data Assimilation
Introduction to Data Assimilation Alan O Neill Data Assimilation Research Centre University of Reading What is data assimilation? Data assimilation is the technique whereby observational data are combined
More informationNumerical Weather prediction at the European Centre for Medium-Range Weather Forecasts
Numerical Weather prediction at the European Centre for Medium-Range Weather Forecasts Time series curves 500hPa geopotential Correlation coefficent of forecast anomaly N Hemisphere Lat 20.0 to 90.0 Lon
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 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 informationRevision of TR-09-25: A Hybrid Variational/Ensemble Filter Approach to Data Assimilation
Revision of TR-9-25: A Hybrid Variational/Ensemble ilter Approach to Data Assimilation Adrian Sandu 1 and Haiyan Cheng 1 Computational Science Laboratory Department of Computer Science Virginia Polytechnic
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 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 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 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 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 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 informationAccepted in Tellus A 2 October, *Correspondence
1 An Adaptive Covariance Inflation Error Correction Algorithm for Ensemble Filters Jeffrey L. Anderson * NCAR Data Assimilation Research Section P.O. Box 3000 Boulder, CO 80307-3000 USA Accepted in Tellus
More informationA HYBRID ENSEMBLE KALMAN FILTER / 3D-VARIATIONAL ANALYSIS SCHEME
A HYBRID ENSEMBLE KALMAN FILTER / 3D-VARIATIONAL ANALYSIS SCHEME Thomas M. Hamill and Chris Snyder National Center for Atmospheric Research, Boulder, Colorado 1. INTRODUCTION Given the chaotic nature of
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 informationTing Lei, Xuguang Wang University of Oklahoma, Norman, OK, USA. Wang and Lei, MWR, Daryl Kleist (NCEP): dual resolution 4DEnsVar
GSI-based four dimensional ensemble-variational (4DEnsVar) data assimilation: formulation and single resolution experiments with real data for NCEP GFS Ting Lei, Xuguang Wang University of Oklahoma, Norman,
More informationComparing Variational, Ensemble-based and Hybrid Data Assimilations at Regional Scales
Comparing Variational, Ensemble-based and Hybrid Data Assimilations at Regional Scales Meng Zhang and Fuqing Zhang Penn State University Xiang-Yu Huang and Xin Zhang NCAR 4 th EnDA Workshop, Albany, NY
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 informationOPERATIONAL IMPLEMENTATION OF VARIATIONAL DATA ASSIMILATION
OPERATIONAL IMPLEMENTATION OF VARIATIONAL DATA ASSIMILATION PIERRE GAUTHIER Meteorological Service of Canada Dorval, Québec, Canada 1. Introduction Over the last few years, the variational form of statistical
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 informationObservation Impact Assessment for Dynamic. Data-Driven Coupled Chaotic System
Applied Mathematical Sciences, Vol. 10, 2016, no. 45, 2239-2248 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2016.65170 Observation Impact Assessment for Dynamic Data-Driven Coupled Chaotic
More informationEnsemble Kalman Filter based snow data assimilation
Ensemble Kalman Filter based snow data assimilation (just some ideas) FMI, Sodankylä, 4 August 2011 Jelena Bojarova Sequential update problem Non-linear state space problem Tangent-linear state space problem
More informationA Comparative Study of 4D-VAR and a 4D Ensemble Kalman Filter: Perfect Model Simulations with Lorenz-96
Tellus 000, 000 000 (0000) Printed 20 October 2006 (Tellus LATEX style file v2.2) A Comparative Study of 4D-VAR and a 4D Ensemble Kalman Filter: Perfect Model Simulations with Lorenz-96 Elana J. Fertig
More informationEnsemble of Data Assimilations and Uncertainty Estimation
Ensemble of Data Assimilations and Uncertainty Estimation Massimo Bonavita ECMWF, Reading, UK Massimo.Bonavita@ecmwf.int ABSTRACT The background error covariance matrix (B) plays a fundamental role in
More information