4. DATA ASSIMILATION FUNDAMENTALS
|
|
- Alban Jacobs
- 6 years ago
- Views:
Transcription
1 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 destruction of knowledge and its restoration by new observations." Leith (1993) ATMOSPHERIC STATE OBSERVATIONS ASSIMILATED MODEL REALITY TIME Ross Bannister, EO and DA, QUEST ES Page 3 of 51
2 The vector notation for fields and data and the need for an a- priori The 'state vector', x The 'observation vector', y u v θ p q λ1 λ n φ 1 φ m l 1 l L u zonal wind field v meridional wind field θ potential temperature p pressure q specific humidity λ longitude φ latitude l vertical level NWP models: > 10 7 elements (5 n m L) y 1 y y 3 y N NWP models: ~10 6 elements The 'forward model', h y = h [x ] + ε No. of obs. << No. of (unknown) elements in x. This is an under-constrained (and inexact) inverse problem. Need to fill-in the missing information with prior knowledge. An 'a-priori' state (a.k.a. 'first guess', 'background', 'forecast') is needed to make the assimilation problem well posed. Ross Bannister, EO and DA, QUEST ES Page 4 of 51
3 Vectors and matrices Vector/matrix notation is a powerful and compact way of dealing with large volumes of data. A matrix operator acts on an input vector to give an output vector, e.g. ( ) x () ( 1 A 11 A 1 A x (1) 1 1N x () (1) x () A = Ax = 1 A A N x (1) x () A N1 A N A NN N x (1) N )( ) Matrix products do not commute in general, e.g. x (3) = ABCx (1) (1) CBAx Some matrices can be inverted (must be 'square' and non-singluar), e.g. ( ) x (1) ( 1 (A 1 ) 11 (A 1 ) 1 (A 1 ) 1 1N x (1) (A = 1 ) 1 (A 1 ) (A 1 ) N x (1) (A 1 ) N N1 (A 1 ) N (A 1 ) NN N ) ( ) x () x () x () (matrices are complicated to invert, ie (A 1 ) ij A 1 ij.) The matrix transpose make rows into columns and columns into rows (also for vectors), e.g. ( ( A 11 A 1 A 1N A 11 A 1 A N1 A A = 1 A A N ), A T A = 1 A A N A N1 A N A NN A 1N A N A ( ) NN x 1 x x =, x T = (x 1, x, x N ) x N ) The inner product ('scalar' or 'dot' product), e.g. () x T x (1) = x () 1 x (1) 1 + x () x (1) + +x () The outer product (a matrix), e.g. x () x (1) T = ( N x (1) N ) x () 1 x (1) 1 x () 1 x (1) x () x (1) 1 x () x (1) x () 1 x (1) N x (1) N x () N x (1) x () x () N x (1) 1 x () N x (1) N = scalar Ross Bannister, EO and DA, QUEST ES Page 5 of 51
4 Early data assimilation ("objective analysis") The method of "successive corrections" Bergthorsson & Doos (1955), Cressman (1959) Analysis obs. (obs-rich regions). Analysis a-priori (obs-poor regions). Simple scheme to develop. Computationally cheap. Use a-priori in absence of observations. obs. grid-point Analysis is a linear combination of nearby observations, and an a-priori. Poor account of error statistics of obs and a-priori. Direct observations only. Anomalous spreading of obs. information. No multivariate relations (e.g. geostrophy). Ross Bannister, EO and DA, QUEST ES Page 6 of 51
5 'Optimal' Interpolation (OI) Introduced in the 1970s - a more powerful formulation of data assimilation xa = xb + K (y h [x B]) K is a rectangular matrix operator (the 'gain matrix'). K depends upon the error covariance matrices B and R, and linearization H. K = BH T (HBH T + R) 1. The Best Linear Unbiased Estimator (BLUE). R : observation error covariance matrix. Describes the error statistics of the observations (see later). B : background error covariance matrix. Describes the error statistics of the a-priori state (see later). H : linearized observation operator. h [x B + δx ] h [x B] + Hδx. Account taken of a-priori and obs. error statistics. Allows assimilation of some indirect obs. Use a-priori in absence of obs. Works as an inverse model. Too expensive for single global solution. Difficult to know B. No consistency with the equations of motion. Ross Bannister, EO and DA, QUEST ES Page 7 of 51
6 Types of errors 1. Random errors Data Arising from Obs. Noise Forecast Stochastic processes in model, init. conds. Assim. Input data E.g. repeated measurement of temperature: truth mean spread (ref. 'error statistics'). Systematic errors Data Arising from Obs. E.g. reading errors Forecast Model formulation, init. conds. Assim. Input data, formulation E.g. As before but with biased thermometer: truth mean spread Biases should be corrected where possible. 3. Representativeness error Data Arising from Assim. Unresolved variability E.g. Interpolation of model grid values to location of observation (forward operator h [x ]): V 1 obs V V 11 V 1 Ross Bannister, EO and DA, QUEST ES Page 8 of 51
7 Probability distribution functions Error statistics are described by a probability density function (PDF). PDFs of random and representativeness errors are often expressed together. They are often approximated by the normal (Gaussian) distribution. A scalar (ie a single piece of information), y P (y) A vector (multiple pieces of information), x x P ( x ) σ x x 1 xn ( ) x = y (y y ) P (y) exp σ Mean : y Variance : σ ( Std. dev. : σ) y ( P (x ) exp 1 (x x ) B 1 (x x )) Mean Covariance : x : B = ( x 1 ) δx 1 δx 1 δx δx 1 δx N δx δx 1 δx δx δx N δx N δx 1 δx N δx δx N = δx δx T where δx = x x Ross Bannister, EO and DA, QUEST ES Page 9 of 51
8 Example of covariances: forecast as the a-priori xb is a forecast and so the equations of motion will influence strongly the covariance patterns. Example: Geostrophic error covariances Geostrophic balance: v = 1 ρf p x u = 1 ρf p y Courtesy, Univ. of Washington Pressure-pressure covariances assumption: δp i δp j = σ exp r ij L L 750 km Ross Bannister, EO and DA, QUEST ES Page 30 of 51
9 Variational data assimilation The 'method of least squares' - simple version J (x ) = (x xb) + (y h [x ]) J : cost function (a scalar) xb : a-priori (background) state y : observations h [x ] : observation operator (forward model) x : variable x A = x min J = Carl Fredrich Gauss "analysis" The 'method of least squares' - considering error statistics J (x ) = 1 (x 1 (y x B) T B 1 (x h [x ]) T R 1 (y xb) + h [x ]) B : background error covariance matrix R : observation error covariance matrix This is the form used in operational weather forecasting, deriving satellite retrievals, etc. Non-Euclidean L norm. Assumes perfect forward model, unbiased data. This is consistent with a Gaussian model of error statistics (next slide). 'Var.' is efficient enough to solve the global problem. Ross Bannister, EO and DA, QUEST ES Page 31 of 51
10 The Bayesian view of data assimilation Bayes' Theorem P(y, x) = P(x y)p(y ) P(x P(x, y) = P(y x)p(x ) y) = P(y x)p(x ) P(y ) P (x ) P (y x) Rev. Thomas Bayes P (x ( y) exp 1 (x x B) T B 1 (x B)) x exp ( 1 (h [x ] y) T R 1 (h [x ] y )) ( exp 1 (x x B) T B 1 (x xb) + 1 (h [x ] y) T R 1 (h )) [x ] y Maximum likelihood Minimum penalty, J J [x ] = 1 (x x B) T B 1 (x xb) + 1 (h [x ] y) T R 1 (h [x ] y) xa = x min J = "analysis" Ross Bannister, EO and DA, QUEST ES Page 3 of 51
11 Minimising the cost function The problem reduces to a (badly conditioned) optimisation problem in dimensional phase space. J [x ] = 1 (x x B) T B 1 (x xb) + 1 (h [x ] y) T R 1 (h [x ] y) x xb xa x 1 Descent algorithms minimize J iteratively. They need the local gradient, x J of the cost function at each iteration. The adjoint method is used to compute the adjoint. 1 The curvature (a.k.a. inverse Hessian, ) at x indicates the error statistics of the analysis. A very badly conditioned problem. ( xj) 1 Ross Bannister, EO and DA, QUEST ES Page 33 of 51 A
12 Algebraic minimization of the cost function Under simplified conditions the cost function can be minimized algebraically. Assume that the linearization of the forward model is reasonable h [x ] h [x B] + H (x xb) J [x ] = 1 (x x B) T B 1 (x xb) + 1 (H (x x B) (y h [x B]) ) T R 1 (H (x xb) (y h [x B]) ) 1. Calculate the gradient vector ( x J = ) J / x 1 J / x = J / x N B 1 (x xb) + H T R 1 (h [x ] y). The special x that has zero gradient minimizes J (this cost function is quadratic and convex) xa x J xa = 0 = xb + (B 1 + H T R 1 H) 1 H T R 1 (y h [x B ]) = xb + BH T (R + HBH T ) 1 (y h [x B ]) This is the OI formula with the BLUE! Ross Bannister, EO and DA, QUEST ES Page 34 of 51
13 Types of data assimilation Sequential data assimilation methods Data assimilation performed at each batch of observations. Model forecast made between batches (for background). E.g. OI, KF*, EnKF*, etc. Explicit formula used for analysis. Very expensive. 1d-Var Data assimilation performed for vertical profile only, where satellite makes observations. Used as a 'pre-main-assimilation' step to produce vertical profiles of model quantities (retrievals) from satellite radiances. * Kalman Filter (KF) and Ensemble Kalman Filter (EnKF). Ross Bannister, EO and DA, QUEST ES Page 35 of 51
14 3d-Var Data assimilation performed every 6 hours. 6 hour model forecast between analysis times (for background). Adequate for re-analysis. Relatively cheap. Observations within ±3 hours are not at analysis time. No dynamical constraint used. 3d-Var. time window 4d-Var Data assimilation performed every 6 (1) hours. 6 (1) hour model forecast between analysis times (for background). Used (e.g.) by ECMWF and Met Office for operational weather forecasting. Model used as a dynamical constraint. Observations are compared to the model trajectory at the correct time. Perfect model assumption. Expensive, but not unfeasible. 4d-Var. time window t = 3 t = 0 t = +3 t = +6 t = 3 t = 0 t = +3 t = +6 3dFGAT '3d First Guess at Time' is half-way between 3d and 4d-Var. 'Strong constraint' - it is possible to use 4d-Var with a model under the 'weak-constraint' formulation. Ross Bannister, EO and DA, QUEST ES Page 36 of 51
15 The 4d-Var cost function 4d-Var. time window yt = 1 yt = yt = 3 J [x 0] = 1 (x 0 xb) T B 1 (x 0 xb) + 1 t (h t [x t] yt) T R 1 t (h t [x t] yt) t = 3 t = 0 t = +3 t = +6 The observation vector comprises subvectors, y for time t. The observation operator ht acts on model state xt. Vary x0 in the minimization - the state at the start of the 4d-Var. cycle. Future states in the cycle are computed with the forecast model, xt = M x. Forward model is the composite operator h. Important issues: t t t M 0 t 0 0 Tangent linear model (and its adjoint) needed and can be difficult to find. Forecast model can be highly non-linear (e.g. sensitive dependence in model's convection scheme - on/off 'switches'). Ross Bannister, EO and DA, QUEST ES Page 37 of 51
16 Assimilation of sequences of satellite images in 4d-Var (Courtesy Samantha Pullen, Met Office) Sequence of observed brightness temperatures Sequence of simulated brightness temperatures Ross Bannister, EO and DA, QUEST ES Page 38 of 51
17 '4d-Var.' demonstration with a double pendulum θ1 m 1 ( x = ) θ θ 1 θ1 l 1 m l θ θ l 3 L = T V d d t ( L θ i) = L θi m 3 V = gm 1 l 1 cos θ1 gm l cos θ1 gm 3 (l cos θ1 + l 3 cos θ) T = ½ m 1 (x1 + y1) + ½ m (x + y) + ½ m 3 (x3 + y3) Demonstrate '4d-Var' with an OSSE - 'Observation System Simulation Experiment'. Also known as a 'twin experiment'. Choose a set of initial conditions and run the model (truth). Add random noise to generate pseudo-observations. Forget the truth and try to recover it by assimilating the observations. Use observations of and only (no observations of and ). θ1 θ θ 1 θ Ross Bannister, EO and DA, QUEST ES Page 39 of 51
18 OSSE demonstration (double pendulum) - truth run 4 3 θ1 Angle θ 1 Rate of change angle time 'Truth' Observations Ross Bannister, EO and DA, QUEST ES Page 40 of 51
19 OSSE demonstration (double pendulum) - '4d-Var.' run 4 3 θ1 Angle θ 1 Rate of change angle time 'Truth' Observations '4d-Var.' analysis Ross Bannister, EO and DA, QUEST ES Page 41 of 51
20 OSSE demonstration (double pendulum) - 'obs insertion' run 4 3 θ1 Angle θ 1 Rate of change angle time 'Truth' Observations 'Insertion' run Ross Bannister, EO and DA, QUEST ES Page 4 of 51
21 Issues with data assimilation Data assimilation is a computer intensive process. For one cycle, 4d-Var. can use up to 100 times more computer power than the forecast. The B-matrix (forecast error covariance matrix in Var.) is difficult to deal with. Assimilation process is very sensitive to B. Least well-known part of data assimilation. In operational data assimilation, B is a matrix. Need to model the B-matrix - use technique of 'control variable transforms'. In reality B is flow dependent. Practically, B is quasi-static. Data assimilation replies on optimality. Issues of suboptimality arise if: Actual error distributions are non-gaussian, B or R are inappropriate. Forward models are inaccurate or are non-linear. Data have biases. Cost function has not converged adequately (in Var.). Assimilation can introduce undesirable imbalances. Quantities not constrained by observations can be poor (e.g. diagnosed quantities): Precipitation. Vertical velocity, etc. Ross Bannister, EO and DA, QUEST ES Page 43 of 51
Fundamentals 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 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 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 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 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 informationOptimal Interpolation ( 5.4) We now generalize the least squares method to obtain the OI equations for vectors of observations and background fields.
Optimal Interpolation ( 5.4) We now generalize the least squares method to obtain the OI equations for vectors of observations and background fields. Optimal Interpolation ( 5.4) We now generalize the
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; 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 informationThe Inversion Problem: solving parameters inversion and assimilation problems
The Inversion Problem: solving parameters inversion and assimilation problems UE Numerical Methods Workshop Romain Brossier romain.brossier@univ-grenoble-alpes.fr ISTerre, Univ. Grenoble Alpes Master 08/09/2016
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 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 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 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 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 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 informationBrian J. Etherton University of North Carolina
Brian J. Etherton University of North Carolina The next 90 minutes of your life Data Assimilation Introit Different methodologies Barnes Analysis in IDV NWP Error Sources 1. Intrinsic Predictability Limitations
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 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 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 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 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 informationIntroduction to data assimilation and least squares methods
Introduction to data assimilation and least squares methods Eugenia Kalnay and many friends University of Maryland October 008 (part 1 Contents (1 Forecasting the weather - we are really getting better!
More informationTangent-linear and adjoint models in data assimilation
Tangent-linear and adjoint models in data assimilation Marta Janisková and Philippe Lopez ECMWF Thanks to: F. Váňa, M.Fielding 2018 Annual Seminar: Earth system assimilation 10-13 September 2018 Tangent-linear
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 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 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 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 informationParallel Algorithms for Four-Dimensional Variational Data Assimilation
Parallel Algorithms for Four-Dimensional Variational Data Assimilation Mie Fisher ECMWF October 24, 2011 Mie Fisher (ECMWF) Parallel 4D-Var October 24, 2011 1 / 37 Brief Introduction to 4D-Var Four-Dimensional
More informationEnsemble prediction and strategies for initialization: Tangent Linear and Adjoint Models, Singular Vectors, Lyapunov vectors
Ensemble prediction and strategies for initialization: Tangent Linear and Adjoint Models, Singular Vectors, Lyapunov vectors Eugenia Kalnay Lecture 2 Alghero, May 2008 Elements of Ensemble Forecasting
More informationObservation impact on data assimilation with dynamic background error formulation
Observation impact on data assimilation with dynamic background error formulation ALEXANDER BECK alexander.beck@univie.ac.at Department of, Univ. Vienna, Austria Thanks to: Martin Ehrendorfer, Patrick
More informationData assimilation concepts and methods March 1999
Data assimilation concepts and methods March 1999 By F. Bouttier and P. Courtier Abstract These training course lecture notes are an advanced and comprehensive presentation of most data assimilation methods
More informationBasic concepts of data assimilation
Basic concepts of data assimilation Numerical Weather Prediction NWP is an initial-boundary value problem: - given an estimate of the present state of the atmosphere (initial conditions), and appropriate
More informationData assimilation for large state vectors
Data assimilation for large state vectors Frédéric Chevallier Laboratoire des Sciences du Climat et de l Environnement Gif-sur-Yvette France Outline Introduction Basic concepts Prior information Analytical
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 informationRepresentation of inhomogeneous, non-separable covariances by sparse wavelet-transformed matrices
Representation of inhomogeneous, non-separable covariances by sparse wavelet-transformed matrices Andreas Rhodin, Harald Anlauf German Weather Service (DWD) Workshop on Flow-dependent aspects of data assimilation,
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 informationDRAGON ADVANCED TRAINING COURSE IN ATMOSPHERE REMOTE SENSING. Inversion basics. Erkki Kyrölä Finnish Meteorological Institute
Inversion basics y = Kx + ε x ˆ = (K T K) 1 K T y Erkki Kyrölä Finnish Meteorological Institute Day 3 Lecture 1 Retrieval techniques - Erkki Kyrölä 1 Contents 1. Introduction: Measurements, models, inversion
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 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 informationComparison of of Assimilation Schemes for HYCOM
Comparison of of Assimilation Schemes for HYCOM Ashwanth Srinivasan, C. Thacker, Z. Garraffo, E. P. Chassignet, O. M. Smedstad, J. Cummings, F. Counillon, L. Bertino, T. M. Chin, P. Brasseur and C. Lozano
More informationThe University of Reading
The University of Reading Radial Velocity Assimilation and Experiments with a Simple Shallow Water Model S.J. Rennie 2 and S.L. Dance 1,2 NUMERICAL ANALYSIS REPORT 1/2008 1 Department of Mathematics 2
More informationUsing Observations at Different Spatial. Scales in Data Assimilation for. Environmental Prediction. Joanne A. Waller
UNIVERSITY OF READING DEPARTMENT OF MATHEMATICS AND STATISTICS Using Observations at Different Spatial Scales in Data Assimilation for Environmental Prediction Joanne A. Waller Thesis submitted for the
More informationDirect assimilation of all-sky microwave radiances at ECMWF
Direct assimilation of all-sky microwave radiances at ECMWF Peter Bauer, Alan Geer, Philippe Lopez, Deborah Salmond European Centre for Medium-Range Weather Forecasts Reading, Berkshire, UK Slide 1 17
More informationMathematical Concepts of Data Assimilation
Mathematical Concepts of Data Assimilation N.K. Nichols 1 Introduction Environmental systems can be realistically described by mathematical and numerical models of the system dynamics. These models can
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 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 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 informationSatellite Observations of Greenhouse Gases
Satellite Observations of Greenhouse Gases Richard Engelen European Centre for Medium-Range Weather Forecasts Outline Introduction Data assimilation vs. retrievals 4D-Var data assimilation Observations
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 informationQuantifying observation error correlations in remotely sensed data
Quantifying observation error correlations in remotely sensed data Conference or Workshop Item Published Version Presentation slides Stewart, L., Cameron, J., Dance, S. L., English, S., Eyre, J. and Nichols,
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 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 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 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 informationDemonstration and Comparison of of Sequential Approaches for Altimeter Data Assimilation in in HYCOM
Demonstration and Comparison of of Sequential Approaches for Altimeter Data Assimilation in in HYCOM A. Srinivasan, E. P. Chassignet, O. M. Smedstad, C. Thacker, L. Bertino, P. Brasseur, T. M. Chin,, F.
More informationReducing the Impact of Sampling Errors in Ensemble Filters
Reducing the Impact of Sampling Errors in Ensemble Filters Jeffrey Anderson NCAR Data Assimilation Research Section The National Center for Atmospheric Research is sponsored by the National Science Foundation.
More 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 informationLocalization in the ensemble Kalman Filter
Department of Meteorology Localization in the ensemble Kalman Filter Ruth Elizabeth Petrie A dissertation submitted in partial fulfilment of the requirement for the degree of MSc. Atmosphere, Ocean and
More 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 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 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 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 informationLinearized Physics (LP): Progress and Issues
Linearized Physics (LP): Progress and Issues Philippe Lopez ECMWF, Reading, UK with thanks to Marta Janisková, Alan Geer and Peter Bauer ECMWF-JCSDA Workshop on Assimilating Satellite Observations of Clouds
More informationNumerical Weather Prediction: Data assimilation. Steven Cavallo
Numerical Weather Prediction: Data assimilation Steven Cavallo Data assimilation (DA) is the process estimating the true state of a system given observations of the system and a background estimate. Observations
More informationFour-Dimensional Ensemble Kalman Filtering
Four-Dimensional Ensemble Kalman Filtering B.R. Hunt, E. Kalnay, E.J. Kostelich, E. Ott, D.J. Patil, T. Sauer, I. Szunyogh, J.A. Yorke, A.V. Zimin University of Maryland, College Park, MD 20742, USA Ensemble
More 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 informationQuantifying observation error correlations in remotely sensed data
Quantifying observation error correlations in remotely sensed data Conference or Workshop Item Published Version Presentation slides Stewart, L., Cameron, J., Dance, S. L., English, S., Eyre, J. and Nichols,
More informationNonlinear variational inverse problems
6 Nonlinear variational inverse problems This chapter considers highly nonlinear variational inverse problems and their properties. More general inverse formulations for nonlinear dynamical models will
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 informationSimple Examples. Let s look at a few simple examples of OI analysis.
Simple Examples Let s look at a few simple examples of OI analysis. Example 1: Consider a scalar prolem. We have one oservation y which is located at the analysis point. We also have a ackground estimate
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 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 informationMasahiro Kazumori, Takashi Kadowaki Numerical Prediction Division Japan Meteorological Agency
Development of an all-sky assimilation of microwave imager and sounder radiances for the Japan Meteorological Agency global numerical weather prediction system Masahiro Kazumori, Takashi Kadowaki Numerical
More informationCoupled atmosphere-ocean data assimilation in the presence of model error. Alison Fowler and Amos Lawless (University of Reading)
Coupled atmosphere-ocean data assimilation in the presence of model error Alison Fowler and Amos Lawless (University of Reading) Introduction Coupled DA methods are being developed to initialise forecasts
More informationA Penalized 4-D Var data assimilation method for reducing forecast error
A Penalized 4-D Var data assimilation method for reducing forecast error M. J. Hossen Department of Mathematics and Natural Sciences BRAC University 66 Mohakhali, Dhaka-1212, Bangladesh I. M. Navon Department
More informationEstimating Observation Impact in a Hybrid Data Assimilation System: Experiments with a Simple Model
Estimating Observation Impact in a Hybrid Data Assimilation System: Experiments with a Simple Model NOAA / NCEP / EMC College Park, MD 20740, USA 10 February 2014 Overview Goal Sensitivity Theory Adjoint
More informationAspects of the practical application of ensemble-based Kalman filters
Aspects of the practical application of ensemble-based Kalman filters Lars Nerger Alfred Wegener Institute for Polar and Marine Research Bremerhaven, Germany and Bremen Supercomputing Competence Center
More informationCurrent Limited Area Applications
Current Limited Area Applications Nils Gustafsson SMHI Norrköping, Sweden nils.gustafsson@smhi.se Outline of talk (contributions from many HIRLAM staff members) Specific problems of Limited Area Model
More informationConditioning of the Weak-Constraint Variational Data Assimilation Problem for Numerical Weather Prediction. Adam El-Said
THE UNIVERSITY OF READING DEPARTMENT OF MATHEMATICS AND STATISTICS Conditioning of the Weak-Constraint Variational Data Assimilation Problem for Numerical Weather Prediction Adam El-Said Thesis submitted
More informationLocalization and Correlation in Ensemble Kalman Filters
Localization and Correlation in Ensemble Kalman Filters Jeff Anderson NCAR Data Assimilation Research Section The National Center for Atmospheric Research is sponsored by the National Science Foundation.
More informationON DIAGNOSING OBSERVATION ERROR STATISTICS WITH LOCAL ENSEMBLE DATA ASSIMILATION
ON DIAGNOSING OBSERVATION ERROR STATISTICS WITH LOCAL ENSEMBLE DATA ASSIMILATION J. A. Waller, S. L. Dance, N. K. Nichols University of Reading 1 INTRODUCTION INTRODUCTION Motivation Only % of observations
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 informationHow is balance of a forecast ensemble affected by adaptive and non-adaptive localization schemes?
1 How is balance of a forecast ensemble affected by adaptive and non-adaptive localization schemes? Ross Bannister Stefano Migliorini, Laura Baker, Ali Rudd NCEO, DIAMET, ESA University of Reading 2 Outline
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 informationFLORIDA STATE UNIVERSITY COLLEGE OF ARTS AND SCIENCES ADAPTIVE OBSERVATIONS IN A 4D-VAR FRAMEWORK APPLIED TO THE NONLINEAR BURGERS EQUATION MODEL
FLORIDA STATE UNIVERSITY COLLEGE OF ARTS AND SCIENCES ADAPTIVE OBSERVATIONS IN A 4D-VAR FRAMEWORK APPLIED TO THE NONLINEAR BURGERS EQUATION MODEL By MD. JAKIR HOSSEN A Thesis submitted to the Department
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 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 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 informationVariational data assimilation of lightning with WRFDA system using nonlinear observation operators
Variational data assimilation of lightning with WRFDA system using nonlinear observation operators Virginia Tech, Blacksburg, Virginia Florida State University, Tallahassee, Florida rstefane@vt.edu, inavon@fsu.edu
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 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 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 informationEnsemble Kalman Filter
Ensemble Kalman Filter Geir Evensen and Laurent Bertino Hydro Research Centre, Bergen, Norway, Nansen Environmental and Remote Sensing Center, Bergen, Norway The Ensemble Kalman Filter (EnKF) Represents
More informationA 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 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 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 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 information[POLS 8500] Review of Linear Algebra, Probability and Information Theory
[POLS 8500] Review of Linear Algebra, Probability and Information Theory Professor Jason Anastasopoulos ljanastas@uga.edu January 12, 2017 For today... Basic linear algebra. Basic probability. Programming
More informationData Assimilation in Geosciences: A highly multidisciplinary enterprise
Data Assimilation in Geosciences: A highly multidisciplinary enterprise Adrian Sandu Computational Science Laboratory Department of Computer Science Virginia Tech Data assimilation fuses information from
More informationNumerical Methods for PDE-Constrained Optimization
Numerical Methods for PDE-Constrained Optimization Richard H. Byrd 1 Frank E. Curtis 2 Jorge Nocedal 2 1 University of Colorado at Boulder 2 Northwestern University Courant Institute of Mathematical Sciences,
More information