Introduction to ensemble forecasting. Eric J. Kostelich
|
|
- Julia Reed
- 5 years ago
- Views:
Transcription
1 Introduction to ensemble forecasting Eric J. Kostelich SCHOOL OF MATHEMATICS AND STATISTICS MSRI Climate Change Summer School July 21, 2008
2 Co-workers: Istvan Szunyogh, Brian Hunt, Edward Ott, Eugenia Kalnay, Jim Yorke and many others! Thanks to: Dave Kuhl Papers, preprints, and codes: eric MSRI Lecture #2 E. Kostelich MATHEMATICS AND STATISTICS 2 / 32
3 Principal papers Preprints: Initial papers: E. Ott et al., Tellus A 56 (2004), I. Szunyogh et al., Tellus A 57 (2005), Refined mathematical implementation: B. R. Hunt, E. K., I. Szunyogh, Physica D 230 (2007) Results with real data: I. Szunyogh, E.K. et al., Tellus A 60 (2008) MSRI Lecture #2 E. Kostelich MATHEMATICS AND STATISTICS 3 / 32
4 Recap from last time In a chaotic process, every point is sensitive Uncertainties in initial conditions grow exponentially (at least for awhile) The weather is chaotic (as far as anyone can tell) The uncertainty in the global weather vector roughly doubles every 2 days Forecast horizon: about 2 weeks MSRI Lecture #2 E. Kostelich MATHEMATICS AND STATISTICS 4 / 32
5 Relevant U. S. organizations The National Oceanographic and Atmospheric Administration (NOAA) is a division of the Department of Commerce The National Centers for Environmental Prediction (NCEP) is the division of NOAA responsible for developing and maintaining weather forecast models Spectrum of models: Global Forecast System (GFS), Regional Spectral Model (RSM), etc. Model data is distributed to local Weather Service offices, which generate public forecast products MSRI Lecture #2 E. Kostelich MATHEMATICS AND STATISTICS 5 / 32
6 Other important modeling efforts NASA develops and maintains its own forecast models International agreements to share forecasts and observations (NCEP, UK Met Office, ECMWF, Canada, Japan, Brazil, etc.) Research community: Weather Research and Forecasting model (WRF) NOAA and the U. S. Navy develop and maintain ocean models Private sector efforts: AccuWeather, airlines, etc. MSRI Lecture #2 E. Kostelich MATHEMATICS AND STATISTICS 6 / 32
7 What do we want to predict? The best long-term forecast is climatology (the mean is the maximum likelihood estimate) Prior to the mid 1960s, the starting initial condition was climatology The U. S. Weather Service defines normal as the average Example: in Phoenix, Arizona, tomorrow s weather will be sunny with 96% probability Exceptional weather often is of greatest interest MSRI Lecture #2 E. Kostelich MATHEMATICS AND STATISTICS 7 / 32
8 What is data assimilation? The process by which empirical measurements are incorporated into a forecast model to refine an estimate of the initial condition The distinction between variables and parameters is a matter of definition Operational weather forecast centers perform data assimilation steps 4 times per day (0Z, 6Z, 12Z, 18Z) Real-time constraints: NCEP allows 20 minutes MSRI Lecture #2 E. Kostelich MATHEMATICS AND STATISTICS 8 / 32
9 Measures of forecast quality One objective measure of goodness: forecast observations A 72-hour forecast today is as accurate as a 36-hour forecast in 1985 Holy grail: 7-day forecasts that are as accurate as 3-day forecasts are now MSRI Lecture #2 E. Kostelich MATHEMATICS AND STATISTICS 9 / 32
10 Many applications besides weather Controls (e.g., airplane autopilots) Ocean and climate models (obviously) Biological models (e.g., Tim Sauer & Steve Schiff) Parameter estimation MSRI Lecture #2 E. Kostelich MATHEMATICS AND STATISTICS 10 / 32
11 Some fundamental problems Naive approach: direct insertion Difficulty: there are usually many more grid points than available measurements Does not account for errors in the measurement Does not exploit correlations between nearby grid points The variables in the model are not necessarily the ones that can be easily measured MSRI Lecture #2 E. Kostelich MATHEMATICS AND STATISTICS 11 / 32
12 Example: Global Forecast System Principal variables in the GFS: natural logarithm of surface pressure virtual temperature divergence and vorticity of the wind field Principal measurements: barometric pressure sensible temperature relative humidity wind speed and direction satellite radiances (complicated!) MSRI Lecture #2 E. Kostelich MATHEMATICS AND STATISTICS 12 / 32
13 Typical 6-hour land surface dataset: 31,310 locations MSRI Lecture #2 E. Kostelich MATHEMATICS AND STATISTICS 13 / 32
14 Typical 6-hour surface marine dataset: 2,642 locations MSRI Lecture #2 E. Kostelich MATHEMATICS AND STATISTICS 14 / 32
15 Typical 6-hour satellite dataset: 53,842 locations MSRI Lecture #2 E. Kostelich MATHEMATICS AND STATISTICS 15 / 32
16 The observation space For these reasons, data assimilation is done in the observation space Given a vector of observations y, interpolate the model state x to the same locations The interpolation operator is denoted H The innovation is y H(x) MSRI Lecture #2 E. Kostelich MATHEMATICS AND STATISTICS 16 / 32
17 Basic idea: Weighted least squares Observations: y R p, y = H(x t ) + ε Observation errors: E(ε) = 0, E(εε T ) = R Model forecast ( background ): x R n, x b = x t + η Model errors: E(η) = 0, E(ηη T ) = P b Goal: minimize the objective function J(x) = [y H(x)] T R 1 [y H(x)]+(x x b ) T P 1 b (x x b) Minimization produces an analysis x a with associated covariance P a MSRI Lecture #2 E. Kostelich MATHEMATICS AND STATISTICS 17 / 32
18 Simplest assumptions The observation errors ε are normally distributed with mean 0 and covariance R Model errors similarly: N(0,P b ) When the underlying model is linear, it can be shown the the minimizer x a of J is unique, unbiased and has minimum variance among all linear estimators Weather models are linear enough over 6-hour intervals, but there is no guarantee of optimality MSRI Lecture #2 E. Kostelich MATHEMATICS AND STATISTICS 18 / 32
19 The dimensionality problem Must evaluate J(x) = [y H(x)] T R 1 [y H(x)]+(x x b ) T P 1 b (x x b) where y R p, x R n Current NCEP operations: p 1.75 million and n 3 billion We need R 1 (p p) and P 1 b (n n) MSRI Lecture #2 E. Kostelich MATHEMATICS AND STATISTICS 19 / 32
20 The computational complexity problem Inversion of a k k matrix is an O(k 3 ) algorithm If a matrix takes 1 sec to invert, then a matrix takes sec R is nearly diagonal if observation errors are mostly uncorrelated P b is not diagonal Computing P b (t + t) from P a (t) requires integration of the tangent linear model MSRI Lecture #2 E. Kostelich MATHEMATICS AND STATISTICS 20 / 32
21 Complexity reduction strategies Localization: Try to do the minimzation over smaller regions of the globe Estimate and precompute P 1 b : Assume that the forecast uncertainty is approximately constant from one day to the next. (Used in all current operational DA systems) Thin the observations and use only the most important ones MSRI Lecture #2 E. Kostelich MATHEMATICS AND STATISTICS 21 / 32
22 Each strategy has drawbacks Assuming P b constant ignores the errors of the day Generally regarded as one of the key impediments to better forecasts The result of sequential assimilation of observations depends on the order of processing Must assure continuity at the boundaries of the smaller regions MSRI Lecture #2 E. Kostelich MATHEMATICS AND STATISTICS 22 / 32
23 The Local Ensemble Transform Kalman Filter (LETKF) Addresses many of these problems Exploits the geometry of uncertainty in chaotic processes to lower the dimension but still account for errors of the day Assimilates all the data at once Uses localization and sets of observations that vary slowly in space to help assure continuity Permits efficient implementation on massively parallel computers MSRI Lecture #2 E. Kostelich MATHEMATICS AND STATISTICS 23 / 32
24 The geometry of forecast uncertainty The size of a typical high- or low-pressure system is about 1000 km 1000 km ( Texas) The GFS, when run at medium (T62) resolution, contains about 3000 grid-point variables in Texas-sized regions Suppose we run k statistically equivalent forecasts What are the singular values of the resulting 3000 k forecast matrix X F? MSRI Lecture #2 E. Kostelich MATHEMATICS AND STATISTICS 24 / 32
25 Correlation and dimensionality Over most Texas-sized regions, one solution looks much like another The columns of X F tend to be highly correlated MSRI Lecture #2 E. Kostelich MATHEMATICS AND STATISTICS 25 / 32
26 Correlation and dimensionality Over most Texas-sized regions, one solution looks much like another The columns of X F tend to be highly correlated so the SVD of X F yields a good rank-r approximation even when r k MSRI Lecture #2 E. Kostelich MATHEMATICS AND STATISTICS 25 / 32
27 Key empirical finding This was a key finding by D. J. Patil et al. PRL 86 (2001), GFS at T62 resolution: 3000 grid variables over typical Texas-sized region Typical ensemble of 100 k 200 forecasts generates a 3000 k forecast matrix X F whose first r singular vectors, 40 r 80, yield an excellent approximation of the forecast uncertainty MSRI Lecture #2 E. Kostelich MATHEMATICS AND STATISTICS 26 / 32
28 The ensemble dimension The ensemble dimension (E-dimension) of an n k matrix is E (s 1 + s s k ) 2 s s s2 k Measures the eccentricity of the ellipse of forecast uncertainty MSRI Lecture #2 E. Kostelich MATHEMATICS AND STATISTICS 27 / 32
29 Example: s 1 = 3.78, s 2 = 3.60, E dim = MSRI Lecture #2 E. Kostelich MATHEMATICS AND STATISTICS 28 / 32
30 Example: s 1 = 19.24, s 2 = 4.35, E dim = MSRI Lecture #2 E. Kostelich MATHEMATICS AND STATISTICS 29 / 32
31 Example: s 1 = 83.65, s 2 = 4.33, E dim = MSRI Lecture #2 E. Kostelich MATHEMATICS AND STATISTICS 30 / 32
32 The key idea behind the LETKF If the E-dimension is much less than the dimension of the overall space, then the distribution is flat The ensemble forecast uncertainty over a typical synoptic region resembles a pancake (at least for short intervals) Reduce the dimensionality of the problem by changing coordinates to the r-dimensional subspace containing most of the forecast uncertainty The dynamics reduces the uncertainty in the remaining directions MSRI Lecture #2 E. Kostelich MATHEMATICS AND STATISTICS 31 / 32
33 Next lecture Outline of the Kalman filter Mathematical details of how we accomplish the dimension reduction Results with operational models and real observations MSRI Lecture #2 E. Kostelich MATHEMATICS AND STATISTICS 32 / 32
The 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 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 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 informationDevelopment of the Local Ensemble Transform Kalman Filter
Development of the Local Ensemble Transform Kalman Filter Istvan Szunyogh Institute for Physical Science and Technology & Department of Atmospheric and Oceanic Science AOSC Special Seminar September 27,
More informationAsynchronous data assimilation
Ensemble Kalman Filter, lecture 2 Asynchronous data assimilation Pavel Sakov Nansen Environmental and Remote Sensing Center, Norway This talk has been prepared in the course of evita-enkf project funded
More 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 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 informationAdaptive ensemble Kalman filtering of nonlinear systems
Adaptive ensemble Kalman filtering of nonlinear systems Tyrus Berry George Mason University June 12, 213 : Problem Setup We consider a system of the form: x k+1 = f (x k ) + ω k+1 ω N (, Q) y k+1 = h(x
More informationAssimilating Nonlocal Observations using a Local Ensemble Kalman Filter
Tellus 000, 000 000 (0000) Printed 16 February 2007 (Tellus LATEX style file v2.2) Assimilating Nonlocal Observations using a Local Ensemble Kalman Filter Elana J. Fertig 1, Brian R. Hunt 1, Edward Ott
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 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 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 informationThe WMO Observation Impact Workshop. lessons for SRNWP. Roger Randriamampianina
The WMO Observation Impact Workshop - developments outside Europe and lessons for SRNWP Roger Randriamampianina Hungarian Meteorological Service (OMSZ) Outline Short introduction of the workshop Developments
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 informationBred vectors: theory and applications in operational forecasting. Eugenia Kalnay Lecture 3 Alghero, May 2008
Bred vectors: theory and applications in operational forecasting. Eugenia Kalnay Lecture 3 Alghero, May 2008 ca. 1974 Central theorem of chaos (Lorenz, 1960s): a) Unstable systems have finite predictability
More informationSatellite Retrieval Assimilation (a modified observation operator)
Satellite Retrieval Assimilation (a modified observation operator) David D. Kuhl With Istvan Szunyogh and Brad Pierce Sept. 21 2009 Weather Chaos Group Meeting Overview Introduction Derivation of retrieval
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 informationESTIMATING CORRELATIONS FROM A COASTAL OCEAN MODEL FOR LOCALIZING AN ENSEMBLE TRANSFORM KALMAN FILTER
ESTIMATING CORRELATIONS FROM A COASTAL OCEAN MODEL FOR LOCALIZING AN ENSEMBLE TRANSFORM KALMAN FILTER Jonathan Poterjoy National Weather Center Research Experiences for Undergraduates, Norman, Oklahoma
More 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 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 informationMultivariate Correlations: Applying a Dynamic Constraint and Variable Localization in an Ensemble Context
Multivariate Correlations: Applying a Dynamic Constraint and Variable Localization in an Ensemble Context Catherine Thomas 1,2,3, Kayo Ide 1 Additional thanks to Daryl Kleist, Eugenia Kalnay, Takemasa
More informationAbstract 2. ENSEMBLE KALMAN FILTERS 1. INTRODUCTION
J5.4 4D ENSEMBLE KALMAN FILTERING FOR ASSIMILATION OF ASYNCHRONOUS OBSERVATIONS T. Sauer George Mason University, Fairfax, VA 22030 B.R. Hunt, J.A. Yorke, A.V. Zimin, E. Ott, E.J. Kostelich, I. Szunyogh,
More 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 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 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 informationWe honor Ed Lorenz ( ) who started the whole new science of predictability
Atmospheric Predictability: From Basic Theory to Forecasting Practice. Eugenia Kalnay Alghero, May 2008, Lecture 1 We honor Ed Lorenz (1917-2008) who started the whole new science of predictability Ed
More informationLocal Predictability of the Performance of an. Ensemble Forecast System
Local Predictability of the Performance of an Ensemble Forecast System Elizabeth Satterfield Istvan Szunyogh University of Maryland, College Park, Maryland To be submitted to JAS Corresponding author address:
More informationWeight interpolation for efficient data assimilation with the Local Ensemble Transform Kalman Filter
QUARTERLY JOURNAL OF THE ROYAL METEOROLOGICAL SOCIETY Published online 17 December 2008 in Wiley InterScience (www.interscience.wiley.com).353 Weight interpolation for efficient data assimilation with
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 informationEnsemble Assimilation of Global Large-Scale Precipitation
Ensemble Assimilation of Global Large-Scale Precipitation Guo-Yuan Lien 1,2 in collaboration with Eugenia Kalnay 2, Takemasa Miyoshi 1,2 1 RIKEN Advanced Institute for Computational Science 2 University
More informationObservation Bias Correction with an Ensemble Kalman Filter
Tellus 000, 000 000 (0000) Printed 10 April 2007 (Tellus LATEX style file v2.2) Observation Bias Correction with an Ensemble Kalman Filter Elana J. Fertig 1, Seung-Jong Baek 2, Brian R. Hunt 1, Edward
More informationThe Local Ensemble Transform Kalman Filter and its implementation on the NCEP global model at the University of Maryland
The Local Ensemble Transform Kalman Filter and its implementation on the NCEP global model at the University of Maryland Istvan Szunyogh (*), Elizabeth A. Satterfield (*), José A. Aravéquia (**), Elana
More informationComparison between Local Ensemble Transform Kalman Filter and PSAS in the NASA finite volume GCM: perfect model experiments
Comparison between Local Ensemble Transform Kalman Filter and PSAS in the NASA finite volume GCM: perfect model experiments Junjie Liu 1*, Elana Judith Fertig 1, and Hong Li 1 Eugenia Kalnay 1, Brian R.
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 informationThe Local Ensemble Transform Kalman Filter and its implementation on the NCEP global model at the University of Maryland
The Local Ensemble Transform Kalman Filter and its implementation on the NCEP global model at the University of Maryland Istvan Szunyogh (*), Elizabeth A. Satterfield (*), José A. Aravéquia (**), Elana
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 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 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 informationWeight interpolation for efficient data assimilation with the Local Ensemble Transform Kalman Filter
Weight interpolation for efficient data assimilation with the Local Ensemble Transform Kalman Filter Shu-Chih Yang 1,2, Eugenia Kalnay 1,3, Brian Hunt 1,3 and Neill E. Bowler 4 1 Department of Atmospheric
More informationOverview of data assimilation in oceanography or how best to initialize the ocean?
Overview of data assimilation in oceanography or how best to initialize the ocean? T. Janjic Alfred Wegener Institute for Polar and Marine Research Bremerhaven, Germany Outline Ocean observing system Ocean
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 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 informationEstimating observation impact without adjoint model in an ensemble Kalman filter
QUARTERLY JOURNAL OF THE ROYAL METEOROLOGICAL SOCIETY Q. J. R. Meteorol. Soc. 134: 1327 1335 (28) Published online in Wiley InterScience (www.interscience.wiley.com) DOI: 1.12/qj.28 Estimating observation
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 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 informationQuarterly Journal of the Royal Meteorological Society !"#$%&'&(&)"&'*'+'*%(,#$,-$#*'."(/'*0'"(#"(1"&)23$)(4#$2#"( 5'$*)6!
!"#$%&'&(&)"&'*'+'*%(,#$,-$#*'."(/'*0'"(#"("&)$)(#$#"( '$*)! "#$%&'()!!"#$%&$'()*+"$,#')+-)%.&)/+(#')0&%&+$+'+#')+&%(! *'&$+,%-./!0)! "! :-(;%/-,(;! '/;!?$@A-//;B!@
More informationObservability, a Problem in Data Assimilation
Observability, Data Assimilation with the Extended Kalman Filter 1 Observability, a Problem in Data Assimilation Chris Danforth Department of Applied Mathematics and Scientific Computation, UMD March 10,
More informationUncertainty in Operational Atmospheric Analyses. Rolf Langland Naval Research Laboratory Monterey, CA
Uncertainty in Operational Atmospheric Analyses 1 Rolf Langland Naval Research Laboratory Monterey, CA Objectives 2 1. Quantify the uncertainty (differences) in current operational analyses of the atmosphere
More informationAnalysis sensitivity calculation in an Ensemble Kalman Filter
Analysis sensitivity calculation in an Ensemble Kalman Filter Junjie Liu 1, Eugenia Kalnay 2, Takemasa Miyoshi 2, and Carla Cardinali 3 1 University of California, Berkeley, CA, USA 2 University of Maryland,
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 informationComparison between Local Ensemble Transform Kalman Filter and PSAS in the NASA finite volume GCM perfect model experiments
Nonlin. Processes Geophys., 15, 645 659, 2008 Author(s) 2008. This work is distributed under the Creative Commons Attribution 3.0 License. Nonlinear Processes in Geophysics Comparison between Local Ensemble
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 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 informationMonthly Weather Review The Hybrid Local Ensemble Transform Kalman Filter
Monthly Weather Review The Hybrid Local Ensemble Transform Kalman Filter --Manuscript Draft-- Manuscript Number: Full Title: Article Type: Corresponding Author: Corresponding Author's Institution: First
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 informationMesoscale meteorological models. Claire L. Vincent, Caroline Draxl and Joakim R. Nielsen
Mesoscale meteorological models Claire L. Vincent, Caroline Draxl and Joakim R. Nielsen Outline Mesoscale and synoptic scale meteorology Meteorological models Dynamics Parametrizations and interactions
More informationStatus of Atmospheric Winds in Relation to Infrasound. Douglas P. Drob Space Science Division Naval Research Laboratory Washington, DC 20375
Status of Atmospheric Winds in Relation to Infrasound Douglas P. Drob Space Science Division Naval Research Laboratory Washington, DC 20375 GOT WINDS? Douglas P. Drob Space Science Division Naval Research
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 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 informationWRF-LETKF The Present and Beyond
November 12, 2012, Weather-Chaos meeting WRF-LETKF The Present and Beyond Takemasa Miyoshi and Masaru Kunii University of Maryland, College Park miyoshi@atmos.umd.edu Co-investigators and Collaborators:
More informationP3.11 A COMPARISON OF AN ENSEMBLE OF POSITIVE/NEGATIVE PAIRS AND A CENTERED SPHERICAL SIMPLEX ENSEMBLE
P3.11 A COMPARISON OF AN ENSEMBLE OF POSITIVE/NEGATIVE PAIRS AND A CENTERED SPHERICAL SIMPLEX ENSEMBLE 1 INTRODUCTION Xuguang Wang* The Pennsylvania State University, University Park, PA Craig H. Bishop
More informationEnsemble Kalman Filter potential
Ensemble Kalman Filter potential Former students (Shu-Chih( Yang, Takemasa Miyoshi, Hong Li, Junjie Liu, Chris Danforth, Ji-Sun Kang, Matt Hoffman), and Eugenia Kalnay University of Maryland Acknowledgements:
More 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 informationSTRONGLY COUPLED ENKF DATA ASSIMILATION
STRONGLY COUPLED ENKF DATA ASSIMILATION WITH THE CFSV2 Travis Sluka Acknowledgements: Eugenia Kalnay, Steve Penny, Takemasa Miyoshi CDAW Toulouse Oct 19, 2016 Outline 1. Overview of strongly coupled DA
More informationNumerical Weather Prediction Chaos, Predictability, and Data Assimilation
July 23, 2013, DA summer school, Reading, UK Numerical Weather Prediction Chaos, Predictability, and Data Assimilation Takemasa Miyoshi RIKEN Advanced Institute for Computational Science Takemasa.Miyoshi@riken.jp
More informationNATS 101 Section 13: Lecture 25. Weather Forecasting Part II
NATS 101 Section 13: Lecture 25 Weather Forecasting Part II NWP s First Baby Steps: Mid-Twentieth Century It wasn t until the development of computers in the 1940s and 1950s that NWP could be even attempted.
More informationarxiv: v1 [physics.ao-ph] 23 Jan 2009
A Brief Tutorial on the Ensemble Kalman Filter Jan Mandel arxiv:0901.3725v1 [physics.ao-ph] 23 Jan 2009 February 2007, updated January 2009 Abstract The ensemble Kalman filter EnKF) is a recursive filter
More informationOptimal Localization for Ensemble Kalman Filter Systems
Journal December of the 2014 Meteorological Society of Japan, Vol. Á. 92, PERIÁÑEZ No. 6, pp. et 585 597, al. 2014 585 DOI:10.2151/jmsj.2014-605 Optimal Localization for Ensemble Kalman Filter Systems
More 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 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 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 informationEnsemble-based Data Assimilation of TRMM/GPM Precipitation Measurements
January 16, 2014, JAXA Joint PI Workshop, Tokyo Ensemble-based Data Assimilation of TRMM/GPM Precipitation Measurements PI: Takemasa Miyoshi RIKEN Advanced Institute for Computational Science Takemasa.Miyoshi@riken.jp
More informationASSESMENT OF THE SEVERE WEATHER ENVIROMENT IN NORTH AMERICA SIMULATED BY A GLOBAL CLIMATE MODEL
JP2.9 ASSESMENT OF THE SEVERE WEATHER ENVIROMENT IN NORTH AMERICA SIMULATED BY A GLOBAL CLIMATE MODEL Patrick T. Marsh* and David J. Karoly School of Meteorology, University of Oklahoma, Norman OK and
More informationCoupled Global-Regional Data Assimilation Using Joint States
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Coupled Global-Regional Data Assimilation Using Joint States Istvan Szunyogh Texas A&M University, Department of Atmospheric
More informationClimate Variables for Energy: WP2
Climate Variables for Energy: WP2 Phil Jones CRU, UEA, Norwich, UK Within ECEM, WP2 provides climate data for numerous variables to feed into WP3, where ESCIIs will be used to produce energy-relevant series
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 informationJ5.8 ESTIMATES OF BOUNDARY LAYER PROFILES BY MEANS OF ENSEMBLE-FILTER ASSIMILATION OF NEAR SURFACE OBSERVATIONS IN A PARAMETERIZED PBL
J5.8 ESTIMATES OF BOUNDARY LAYER PROFILES BY MEANS OF ENSEMBLE-FILTER ASSIMILATION OF NEAR SURFACE OBSERVATIONS IN A PARAMETERIZED PBL Dorita Rostkier-Edelstein 1 and Joshua P. Hacker The National Center
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 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 informationData assimilation : Basics and meteorology
Data assimilation : Basics and meteorology Olivier Talagrand Laboratoire de Météorologie Dynamique, École Normale Supérieure, Paris, France Workshop on Coupled Climate-Economics Modelling and Data Analysis
More informationA Local Ensemble Kalman Filter for Atmospheric Data Assimilation
Tellus 000, 000 000 (0000) Printed 1 April 2004 (Tellus LATEX style file v2.2) A Local Ensemble Kalman Filter for Atmospheric Data Assimilation By EDWARD OTT 1, BRIAN R. HUNT 2, ISTVAN SZUNYOGH 3, ALEKSEY
More informationJ9.4 ERRORS OF THE DAY, BRED VECTORS AND SINGULAR VECTORS: IMPLICATIONS FOR ENSEMBLE FORECASTING AND DATA ASSIMILATION
J9.4 ERRORS OF THE DAY, BRED VECTORS AND SINGULAR VECTORS: IMPLICATIONS FOR ENSEMBLE FORECASTING AND DATA ASSIMILATION Shu-Chih Yang 1, Matteo Corazza and Eugenia Kalnay 1 1 University of Maryland, College
More informationA local ensemble transform Kalman filter data assimilation system for the NCEP global model
Tellus 000, 000 000 (0000) Printed 14 December 2006 (Tellus LATEX style file v2.2) A local ensemble transform Kalman filter data assimilation system for the NCEP global model By ISTVAN SZUNYOGH 1, ERIC
More informationTropical Cyclone Ensemble Data Assimilation
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Tropical Cyclone Ensemble Data Assimilation Istvan Szunyogh Texas A&M University, Department of Atmospheric Sciences, 3150
More informationInterpretation of two error statistics estimation methods: 1 - the Derozier s method 2 the NMC method (lagged forecast)
Interpretation of two error statistics estimation methods: 1 - the Derozier s method 2 the NMC method (lagged forecast) Richard Ménard, Yan Yang and Yves Rochon Air Quality Research Division Environment
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 informationIntroduction to Meteorology and Weather Forecasting
Introduction to Meteorology and Weather Forecasting ENVI1400 : Meteorology and Forecasting : lecture 1 2 040909 ENVI1400 : Meteorology and Forecasting : lecture 1 3 040914 ENVI1400 : Meteorology and Forecasting
More informationParameter Estimation in EnKF: Surface Fluxes of Carbon, Heat, Moisture and Momentum
Parameter Estimation in EnKF: Surface Fluxes of Carbon, Heat, Moisture and Momentum *Ji-Sun Kang, *Eugenia Kalnay, *Takemasa Miyoshi, + Junjie Liu, # Inez Fung, *Kayo Ide *University of Maryland, College
More informationLATE REQUEST FOR A SPECIAL PROJECT
LATE REQUEST FOR A SPECIAL PROJECT 2016 2018 MEMBER STATE: Italy Principal Investigator 1 : Affiliation: Address: E-mail: Other researchers: Project Title: Valerio Capecchi LaMMA Consortium - Environmental
More informationDIAGNOSING OBSERVATION ERROR STATISTICS FOR NUMERICAL WEATHER PREDICTION
DIAGNOSING OBSERVATION ERROR STATISTICS FOR NUMERICAL WEATHER PREDICTION J. Waller, S. Dance, N. Nichols (University of Reading) D. Simonin, S. Ballard, G. Kelly (Met Office) EMS Annual Meeting: European
More informationGenerating climatological forecast error covariance for Variational DAs with ensemble perturbations: comparison with the NMC method
Generating climatological forecast error covariance for Variational DAs with ensemble perturbations: comparison with the NMC method Matthew Wespetal Advisor: Dr. Eugenia Kalnay UMD, AOSC Department March
More information4D-Var or Ensemble Kalman Filter? TELLUS A, in press
4D-Var or Ensemble Kalman Filter? Eugenia Kalnay 1 *, Hong Li 1, Takemasa Miyoshi 2, Shu-Chih Yang 1, and Joaquim Ballabrera-Poy 3 1 University of Maryland, College Park, MD, 20742-2425 2 Numerical Prediction
More informationA 3DVAR Land Data Assimilation Scheme: Part 2, Test with ECMWF ERA-40
A 3DVAR Land Data Assimilation Scheme: Part 2, Test with ECMWF ERA-40 Lanjun Zou 1 * a,b,c Wei Gao a,d Tongwen Wu b Xiaofeng Xu b Bingyu Du a,and James Slusser d a Sino-US Cooperative Center for Remote
More informationMADIS Airlines for America Briefing
MADIS Airlines for America Briefing Meteorological Assimilated Data Ingest System (MADIS) FPAW Briefing Steve Pritchett NWS Aircraft Based Observations Program Manager MADIS Defined MADIS is a meteorological
More informationGuo-Yuan Lien*, Eugenia Kalnay, and Takemasa Miyoshi University of Maryland, College Park, Maryland 2. METHODOLOGY
9.2 EFFECTIVE ASSIMILATION OF GLOBAL PRECIPITATION: SIMULATION EXPERIMENTS Guo-Yuan Lien*, Eugenia Kalnay, and Takemasa Miyoshi University of Maryland, College Park, Maryland 1. INTRODUCTION * Precipitation
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 informationJi-Sun Kang. Pr. Eugenia Kalnay (Chair/Advisor) Pr. Ning Zeng (Co-Chair) Pr. Brian Hunt (Dean s representative) Pr. Kayo Ide Pr.
Carbon Cycle Data Assimilation Using a Coupled Atmosphere-Vegetation Model and the LETKF Ji-Sun Kang Committee in charge: Pr. Eugenia Kalnay (Chair/Advisor) Pr. Ning Zeng (Co-Chair) Pr. Brian Hunt (Dean
More informationThe benefits and developments in ensemble wind forecasting
The benefits and developments in ensemble wind forecasting Erik Andersson Slide 1 ECMWF European Centre for Medium-Range Weather Forecasts Slide 1 ECMWF s global forecasting system High resolution forecast
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 informationThe Use of a Self-Evolving Additive Inflation in the CNMCA Ensemble Data Assimilation System
The Use of a Self-Evolving Additive Inflation in the CNMCA Ensemble Data Assimilation System Lucio Torrisi and Francesca Marcucci CNMCA, Italian National Met Center Outline Implementation of the LETKF
More information