A strategy to optimize the use of retrievals in data assimilation
|
|
- Jean Sims
- 5 years ago
- Views:
Transcription
1 A strategy to optimize the use of retrievals in data assimilation R. Hoffman, K. Cady-Pereira, J. Eluszkiewicz, D. Gombos, J.-L. Moncet, T. Nehrkorn, S. Greybush 2, K. Ide 2, E. Kalnay 2, M. J. Hoffman 3, and R. J. Wilson 4 AER, 2 U. Maryland, 3 John Hopkins, 4 GFDL Friday, 6 May, 20 Ross N. Hoffman, Vice President, R&D Atmospheric and Environmental Research, 20 Atmospheric and Environmental Research (AER). All rights reserved.
2 Introduction Is it better to use radiances or retrievals in data assimilation? Using retrievals reduces complexity in the H-operator, reduces data volume, allows arbitrary cloud clearing and retrieval methods and makes the assimilation system more modular. Radiances are available sooner, have uncorrelated errors, and are independent of the prior. Using EOFs from the retrieval scheme can reduce data volume and reduces vertical interpolation errors. The averaging kernel (AK) concept allows removing the influence of the prior (allowing for interactive retrievals) and rotating to a representation where the obs. errors are uncorrelated. 20 Atmospheric and Environmental Research (AER). All rights reserved. 2
3 Examples Work in progress from on-going Mars data assimilation experiments using GFDL Mars GCM TES and MCS radiances and retrievals LETKF ensemble Kalman filter method Examples from TES retrievals OSS forward model Rodgers type iterative physical retrieval scheme Inverse Methods for Atmospheric Sounding: Theory and Practice Temperature profiles only [for now] 20 Atmospheric and Environmental Research (AER). All rights reserved. 3
4 Retrievals Retrievals, ˆx (22 of 887 total retrievals) 22 retrievals sampled from 6 hour period Pressure (hpa) Temperature (K) 20 Atmospheric and Environmental Research (AER). All rights reserved. 4
5 Mars Retrieval Locations and Topography 887 retrieval loca7ons 2"#$)$% 6+7$)89% 3"++'$%!"#)0'#)$% Hellas Plani7a 4'++"$% 5"$)0% TES data at the northern fall equinox (Ls=80) of Mars Year 24 (MY24) 20 Atmospheric and Environmental Research (AER). All rights reserved. 5
6 Forward model vertical grid Retrieval data pressure levels (linear) Retrieval data pressure levels (log) Temperature (x a ) profile (linear) p ln(p) T a Temperature (K) pressure levels Pressure level index 20 Atmospheric and Environmental Research (AER). All rights reserved. 6
7 Temperature EOFs Temperature EOFs EOF EOF 2 EOF 3 EOF 4 EOF 5 EOF 6 Pressure level 0 0 N Eigenvalue Variance Sum % 89.6% % 95.% % 98.2% 4 25.% 99.3% % 99.6% % 99.8% From the prior covariance S a 20 Atmospheric and Environmental Research (AER). All rights reserved. 7
8 Averaging kernel method (after Rodgers) ˆx = Ax + ( I! A) x a + G y " y. Here ˆx is the retrieval, x is the true state vector, and x a is the prior. A is the averaging kernel.! y includes measurement error and forward model error which is unbiased with covariance matrix, S!. G y is the sensitivity of retrieval to the radiances Here S a is the prior covariance and the retrieval error covariance is determined from K =!F G y =!ˆx!y = ŜK T S # " = S a K T (KS a K T + S " ) # Ŝ! = K T S "! K + S a! And is the Jacobian of the forward problem, F, evaluated at the solution.!x Finally ( ) A = G y K = Ŝ K T S! " K ˆx 20 Atmospheric and Environmental Research (AER). All rights reserved. 8
9 Interface: required data for assimilation Inputs to the retrieval Prior mean and covariance and radiance error covariance Outputs from the retrieval Retrieval and Jacobian x a, S a, and S! OR, if the channel set is large Prior mean Retrieval, posterior covariance, and averaging kernel ˆx, ˆx and K x a Ŝ, and A 20 Atmospheric and Environmental Research (AER). All rights reserved. 9
10 Jacobian of the forward problem F 0 Jacobian #374, P 0 =0.hPa K =!F!x Pressure (hpa) Wavenumber (cm ) 20 Atmospheric and Environmental Research (AER). All rights reserved. 0
11 Posterior correlations from Ŝ 0 Correlations for Ŝ # Pressure (hpa) Pressure (hpa) 20 Atmospheric and Environmental Research (AER). All rights reserved.
12 Data assimilation interface Within the data assimilation, define the new observation as and the new observation operator by We will now be comparing observed and simulated quantities with the same degree of smoothing. Now the observation increments are unbiased if the ŷ A = ˆx! ( I! A) x a, y A = Ax. ŷ A! y A = G y " y,! y S m = G y S! G y T = AŜ are unbiased, and have covariance 20 Atmospheric and Environmental Research (AER). All rights reserved. 2
13 Averaging kernel, A Rows of A R (eq ) for # 374, P 0 =0.hPa Pressure (hpa) Atmospheric and Environmental Research (AER). All rights reserved. 3
14 Old vs. New observations ˆx ŷ A = ˆx! I! A Retrievals, ˆx ŷ a eq(5) ( ) x a Pressure (hpa) 0 0 Pressure (hpa) Temperature (K) Atmospheric and Environmental Research (AER). All rights reserved. 4
15 Forecast as prior (interactive retrievals) Since the observation increment ŷ A! y A = G y " y, does not depend on the prior, we may use the ensemble background mean as and the ensemble sample covariance as S a. x a But G y does depend on S a, the prior covariance. 20 Atmospheric and Environmental Research (AER). All rights reserved. 5
16 Rotated averaging kernel method S m!/2 If we scale by we effectively rotate to a space where the observation errors are unbiased, uncorrelated, and have unit variance. The new observation is ŷ A = S m!/2 and the new observation operator by ( ˆx! ( I! A) x a ), y A = S m!/2 Ax = A R x, S m is covariance of the errors of the new observakons where A R = S m!/2 A. R is for rotated 20 Atmospheric and Environmental Research (AER). All rights reserved. 6
17 Rotated averaging kernel A R Rows of A R (eq ) for # 374, P 0 =0.hPa Rows of A R (eq ) for # 374, P 0 =0.hPa A Pressure (hpa) 0 0 A R = S m!/2 A Atmospheric and Environmental Research (AER). All rights reserved. 7
18 New vs. rotated observations!/2 ŷ A = ˆx! I! A ŷ A = S m ( ) x a ŷ a eq(5) ( ˆx! I! A ) ŷ a2 eq(9) ( ) x a Pressure (hpa) 0 0 Pressure (hpa) Atmospheric and Environmental Research (AER). All rights reserved. 8
19 EOF analysis of the retrieval method Subtract from to get a version in terms of anomalies, We now project this version into EOF coefficient space. For justification, the filtering step is normally done within the retrieval algorithm initially and at each step of iteration. If we define then we have with x = Ax + I! A ( ) x ˆx = Ax + ( I! A) x a + G y " y x ˆ! = A x! + ( I " A) x a! + G y # y.!a = E T AE ˆ! =! A! + ( I " A! )! a + #!! " = E T G y! y ˆ! = E T ( ˆx " x )! a = E T ( x a " x ) 20 Atmospheric and Environmental Research (AER). All rights reserved. 9
20 Rotated EOF averaging kernel method! a S! = E T S m E S! "/2 has covariance given by. We scale by to rotate to a space where the observation errors are unbiased, uncorrelated, and have unit variance. The new observation as ŷ A = S ( "/2! ˆ! " ( I "!A )! ) a The new observation operator by y A =! A R where we have defined the rotated EOF-space AK, ( x! x ), "/2!A R = S!!AE T. Or work in EOF coefficient space:!a R = S! "/2!A y A =! A R! 20 Atmospheric and Environmental Research (AER). All rights reserved. 20
21 EOF projected (rotated) averaging kernel ÃE T for #374, P 0 =0.hPa Rows of à R (eq 28) for #374, P 0 =0.hPa!AE T !A R = S! "/2!AE T Pressure (hpa) Atmospheric and Environmental Research (AER). All rights reserved. 2
22 Old vs. EOF-projected observations ˆx Retrievals, ˆx ŷa = S ( "/2! ˆ! " ( I "!A )! ) a EOF ŷ a3 eq(26) 2 EOF 2 Pressure (hpa) 0 0 EOF # 3 4 EOF 3 EOF 4 5 EOF Temperature (K) 6 EOF Atmospheric and Environmental Research (AER). All rights reserved. 22
23 Summary Convert standard retrievals into "observations" with expected errors that should be zero mean, uncorrelated, and unit variance, and independent of the background or prior. Define a corresponding obs-function (or H-operator) that is a weighted sum of the temperatures on the radiative transfer model vertical grid. No changes to the assimilation method are needed, except to interpolate to the radiative transfer model vertical grid and to calculate the weighted sum. Projecting onto EOFs used by the retrieval can reduce the number of observations Based on ideas in Rodgers' book, "Inverse Methods for Atmospheric Sounding: Theory and Practice. Next step: Try it out in our Mars LETKF data assimilation system. For details visit h<p://arxiv.org/abs/ Atmospheric and Environmental Research (AER). All rights reserved. 23
24 Vertical interpolation We want to determine the EOF coefficients that best fit the temperatures from the dynamical model. Define an interpolation V that transforms x from the retrieval grid to the dynamical model grid. Then the best fit solution is E T V T VE! = E T V T (x b " V x) 20 Atmospheric and Environmental Research (AER). All rights reserved. 24
Introduction 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 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 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 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 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 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 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 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 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 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 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 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 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 informationData Assimilation Development for the FV3GFSv2
Data Assimilation Development for the FV3GFSv2 Catherine Thomas 1, 2, Rahul Mahajan 1, 2, Daryl Kleist 2, Emily Liu 3,2, Yanqiu Zhu 1, 2, John Derber 2, Andrew Collard 1, 2, Russ Treadon 2, Jeff Whitaker
More informationEnsemble Kalman Filter potential
Ensemble Kalman Filter potential Former students (Shu-Chih( Yang, Takemasa Miyoshi, Hong Li, Junjie Liu, Chris Danforth, Ji-Sun Kang, Matt Hoffman), and Eugenia Kalnay University of Maryland Acknowledgements:
More informationData Assimilation: 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 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 informationIntroduction to Data Assimilation. Saroja Polavarapu Meteorological Service of Canada University of Toronto
Introduction to Data Assimilation Saroja Polavarapu Meteorological Service of Canada University of Toronto GCC Summer School, Banff. May 22-28, 2004 Outline of lectures General idea Numerical weather prediction
More informationNumerical Weather Prediction: 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 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 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 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 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 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 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 informationUncertainty in Models of the Ocean and Atmosphere
Uncertainty in Models of the Ocean and Atmosphere Robert N. Miller College of Oceanic and Atmospheric Sciences Oregon State University Corvallis, OR 97330 Uncertainty in Models of the Ocean and Atmosphere
More informationInverse Theory. COST WaVaCS Winterschool Venice, February Stefan Buehler Luleå University of Technology Kiruna
Inverse Theory COST WaVaCS Winterschool Venice, February 2011 Stefan Buehler Luleå University of Technology Kiruna Overview Inversion 1 The Inverse Problem 2 Simple Minded Approach (Matrix Inversion) 3
More informationData Assimilation Research Testbed Tutorial
Data Assimilation Research Testbed Tutorial Section 2: How should observations of a state variable impact an unobserved state variable? Multivariate assimilation. Single observed variable, single unobserved
More 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 informationSpatial and inter-channel observation error characteristics for AMSU-A and IASI and applications in the ECMWF system
Spatial and inter-channel observation error characteristics for AMSU-A and IASI and applications in the system Niels Bormann, Andrew Collard, Peter Bauer Outline 1) Estimation of observation errors AMSU-A
More informationComparison of 3D-Var and LETKF in an Atmospheric GCM: SPEEDY
Comparison of 3D-Var and LEKF in an Atmospheric GCM: SPEEDY Catherine Sabol Kayo Ide Eugenia Kalnay, akemasa Miyoshi Weather Chaos, UMD 9 April 2012 Outline SPEEDY Formulation Single Observation Eperiments
More informationRecent Advances in EnKF
Recent Advances in EnKF Former students (Shu-Chih( Yang, Takemasa Miyoshi, Hong Li, Junjie Liu, Chris Danforth, Ji-Sun Kang, Matt Hoffman, Steve Penny, Steve Greybush), and Eugenia Kalnay University of
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 informationDART Tutorial Part II: How should observa'ons impact an unobserved state variable? Mul'variate assimila'on.
DART Tutorial Part II: How should observa'ons impact an unobserved state variable? Mul'variate assimila'on. UCAR The Na'onal Center for Atmospheric Research is sponsored by the Na'onal Science Founda'on.
More informationError (E) Covariance Diagnostics. in Ensemble Data Assimilation. Kayo Ide, Daryl Kleist, Adam Kubaryk University of Maryland
Error (E) Covariance Diagnostics in Ensemble Data Assimilation Kayo Ide, Daryl Kleist, Adam Kubaryk University of Maryland 7 th EnKF Data Assimilation Workshop May 2016 Outline Part I. Background: E-diagnostics
More 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 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 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 informationRelationship between Singular Vectors, Bred Vectors, 4D-Var and EnKF
Relationship between Singular Vectors, Bred Vectors, 4D-Var and EnKF Eugenia Kalnay and Shu-Chih Yang with Alberto Carrasi, Matteo Corazza and Takemasa Miyoshi ECODYC10, Dresden 28 January 2010 Relationship
More informationData Assimilation Research Testbed Tutorial
Data Assimilation Research Testbed Tutorial Section 3: Hierarchical Group Filters and Localization Version 2.: September, 26 Anderson: Ensemble Tutorial 9//6 Ways to deal with regression sampling error:
More 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 informationRecent Data Assimilation Activities at Environment Canada
Recent Data Assimilation Activities at Environment Canada Major upgrade to global and regional deterministic prediction systems (now in parallel run) Sea ice data assimilation Mark Buehner Data Assimilation
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 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 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 informationEnhancing information transfer from observations to unobserved state variables for mesoscale radar data assimilation
Enhancing information transfer from observations to unobserved state variables for mesoscale radar data assimilation Weiguang Chang and Isztar Zawadzki Department of Atmospheric and Oceanic Sciences Faculty
More informationPlans for the Assimilation of Cloud-Affected Infrared Soundings at the Met Office
Plans for the Assimilation of Cloud-Affected Infrared Soundings at the Met Office Ed Pavelin and Stephen English Met Office, Exeter, UK Abstract A practical approach to the assimilation of cloud-affected
More informationExtended Kalman Filter Tutorial
Extended Kalman Filter Tutorial Gabriel A. Terejanu Department of Computer Science and Engineering University at Buffalo, Buffalo, NY 14260 terejanu@buffalo.edu 1 Dynamic process Consider the following
More informationP1.12 MESOSCALE VARIATIONAL ASSIMILATION OF PROFILING RADIOMETER DATA. Thomas Nehrkorn and Christopher Grassotti *
P1.12 MESOSCALE VARIATIONAL ASSIMILATION OF PROFILING RADIOMETER DATA Thomas Nehrkorn and Christopher Grassotti * Atmospheric and Environmental Research, Inc. Lexington, Massachusetts Randolph Ware Radiometrics
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 informationStatistical Validation Assessment Model for Advanced Sounding Products
Statistical Validation Assessment Model for Advanced Sounding Products N. Pougatchev T. August, X. Calbet, T. Hultberg, O. Oduleye P. Schlüssel, B. StillerK. St. Germain and G. Bingham Hyperspectral Workshop
More informationDART_LAB Tutorial Section 5: Adaptive Inflation
DART_LAB Tutorial Section 5: Adaptive Inflation UCAR 14 The National Center for Atmospheric Research is sponsored by the National Science Foundation. Any opinions, findings and conclusions or recommendations
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 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 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 informationAssimilation of cloud/precipitation data at regional scales
Assimilation of cloud/precipitation data at regional scales Thomas Auligné National Center for Atmospheric Research auligne@ucar.edu Acknowledgments to: Steven Cavallo, David Dowell, Aimé Fournier, Hans
More informationData Assimilation for Tropospheric CO. Avelino F. Arellano, Jr. Atmospheric Chemistry Division National Center for Atmospheric Research
Data Assimilation for Tropospheric CO Avelino F. Arellano, Jr. Atmospheric Chemistry Division National Center for Atmospheric Research Caveat: Illustrative rather than quantitative, applied rather than
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 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 informationAnalysis Methods in Atmospheric and Oceanic Science
Analysis Methods in Atmospheric and Oceanic Science Copyright 2016 University of Maryland. 1 AOSC 652 Root Finding & Function Minimization Week 7, Day 2 12 Oct 2016 2 Method #3: Newton-Raphson Root is
More informationDynamic Inference of Background Error Correlation between Surface Skin and Air Temperature
Dynamic Inference of Background Error Correlation between Surface Skin and Air Temperature Louis Garand, Mark Buehner, and Nicolas Wagneur Meteorological Service of Canada, Dorval, P. Quebec, Canada Abstract
More informationQuarterly Journal of the Royal Meteorological Society !"#$%&'&(&)"&'*'+'*%(,#$,-$#*'."(/'*0'"(#"(1"&)23$)(4#$2#"( 5'$*)6!
!"#$%&'&(&)"&'*'+'*%(,#$,-$#*'."(/'*0'"(#"("&)$)(#$#"( '$*)! "#$%&'()!!"#$%&$'()*+"$,#')+-)%.&)/+(#')0&%&+$+'+#')+&%(! *'&$+,%-./!0)! "! :-(;%/-,(;! '/;!?$@A-//;B!@
More informationIASI Level 2 Processing
IASI Level Processing Peter Schlüssel EUMESA Soundings from High Spectral esolution Observations Madison 06-08 May 003 EUM.EPS.SYS..03.00 Issue Slide: Outline IASI instrument properties IASI Level processor
More informationQuarterly Journal of the Royal Meteorological Society
Quarterly Journal of the Royal Meteorological Society Effects of sequential or simultaneous assimilation of observations and localization methods on the performance of the ensemble Kalman filter Journal:
More information(Toward) Scale-dependent weighting and localization for the NCEP GFS hybrid 4DEnVar Scheme
(Toward) Scale-dependent weighting and localization for the NCEP GFS hybrid 4DEnVar Scheme Daryl Kleist 1, Kayo Ide 1, Rahul Mahajan 2, Deng-Shun Chen 3 1 University of Maryland - Dept. of Atmospheric
More information16. Data Assimilation Research Team
16. Data Assimilation Research Team 16.1. Team members Takemasa Miyoshi (Team Leader) Shigenori Otsuka (Postdoctoral Researcher) Juan J. Ruiz (Visiting Researcher) Keiichi Kondo (Student Trainee) Yukiko
More informationCIMSS Hyperspectral IR Sounding Retrieval (CHISR) Processing & Applications
CIMSS Hyperspectral IR Sounding Retrieval (CHISR) Processing & Applications Jun Li @, Elisabeth Weisz @, Jinlong Li @, Hui Liu #, Timothy J. Schmit &, Jason Otkin @ and many other CIMSS collaborators @Cooperative
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 informationSimple Linear Regression for the MPG Data
Simple Linear Regression for the MPG Data 2000 2500 3000 3500 15 20 25 30 35 40 45 Wgt MPG What do we do with the data? y i = MPG of i th car x i = Weight of i th car i =1,...,n n = Sample Size Exploratory
More informationEnvironment Canada s Regional Ensemble Kalman Filter
Environment Canada s Regional Ensemble Kalman Filter May 19, 2014 Seung-Jong Baek, Luc Fillion, Kao-Shen Chung, and Peter Houtekamer Meteorological Research Division, Environment Canada, Dorval, Quebec
More 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 informationChallenges in Inverse Modeling and Data Assimilation of Atmospheric Constituents
Challenges in Inverse Modeling and Data Assimilation of Atmospheric Constituents Ave Arellano Atmospheric Chemistry Division National Center for Atmospheric Research NCAR is sponsored by the National Science
More informationCoupled Ocean-Atmosphere Assimilation
Coupled Ocean-Atmosphere Assimilation Shu-Chih Yang 1, Eugenia Kalnay 2, Joaquim Ballabrera 3, Malaquias Peña 4 1:Department of Atmospheric Sciences, National Central University 2: Department of Atmospheric
More informationImproved Use of AIRS Data at ECMWF
Improved Use of AIRS Data at ECMWF A.D. Collard, A.P. McNally European Centre for Medium-Range Weather Forecasts, Reading, U.K. W.W. Wolf QSS Group, Inc., NOAA Science Center, 5200 Auth Road, Camp Springs
More information(Statistical Forecasting: with NWP). Notes from Kalnay (2003), appendix C Postprocessing of Numerical Model Output to Obtain Station Weather Forecasts
35 (Statistical Forecasting: with NWP). Notes from Kalnay (2003), appendix C Postprocessing of Numerical Model Output to Obtain Station Weather Forecasts If the numerical model forecasts are skillful,
More information8. Kalman Filtering. ATM 552 Notes Kalman Filtering: Chapter 8 page 231
AM 552 Notes Kalman Filtering: Chapter 8 page 231 8. Kalman Filtering In this section we will discuss the technique called Kalman filtering. Kalman filtering is a technique for estimating the state of
More informationEnhanced Use of Radiance Data in NCEP Data Assimilation Systems
Enhanced Use of Radiance Data in NCEP Data Assimilation Systems John C. Derber*, Paul VanDelst #, XiuJuan Su &, Xu Li &, Kozo Okamoto % and Russ Treadon* Introduction *NOAA/NWS/NCEP/EMC # CIMSS/UW-Madison
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 informationRetrieval of upper tropospheric humidity from AMSU data. Viju Oommen John, Stefan Buehler, and Mashrab Kuvatov
Retrieval of upper tropospheric humidity from AMSU data Viju Oommen John, Stefan Buehler, and Mashrab Kuvatov Institute of Environmental Physics, University of Bremen, Otto-Hahn-Allee 1, 2839 Bremen, Germany.
More informationConfiguration of All-sky Microwave Radiance Assimilation in the NCEP's GFS Data Assimilation System
Configuration of All-sky Microwave Radiance Assimilation in the NCEP's GFS Data Assimilation System Yanqiu Zhu 1, Emily Liu 1, Rahul Mahajan 1, Catherine Thomas 1, David Groff 1, Paul Van Delst 1, Andrew
More informationWegCenter/UniGraz Technical Report for FFG-ALR No. 2/2007
Wegener Center for Climate and Global Change University of Graz Leechgasse 25, A-8010 Graz, Austria Atmospheric Remote Sensing and Climate System Research Group MULTICLIM Project Report WegCenter/UniGraz
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 informationTwo-dimensional characterization of atmospheric prole retrievals from limb sounding observations
Journal of Quantitative Spectroscopy & Radiative Transfer 86 (2004) 45 71 www.elsevier.com/locate/jqsrt Two-dimensional characterization of atmospheric prole retrievals from limb sounding observations
More informationVariable localization in an Ensemble Kalman Filter: application to the carbon cycle data assimilation
1 Variable localization in an Ensemble Kalman Filter: 2 application to the carbon cycle data assimilation 3 4 1 Ji-Sun Kang (jskang@atmos.umd.edu), 5 1 Eugenia Kalnay(ekalnay@atmos.umd.edu), 6 2 Junjie
More informationAircraft Validation of Infrared Emissivity derived from Advanced InfraRed Sounder Satellite Observations
Aircraft Validation of Infrared Emissivity derived from Advanced InfraRed Sounder Satellite Observations Robert Knuteson, Fred Best, Steve Dutcher, Ray Garcia, Chris Moeller, Szu Chia Moeller, Henry Revercomb,
More informationApplications of an ensemble Kalman Filter to regional ocean modeling associated with the western boundary currents variations
Applications of an ensemble Kalman Filter to regional ocean modeling associated with the western boundary currents variations Miyazawa, Yasumasa (JAMSTEC) Collaboration with Princeton University AICS Data
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 informationRecent developments for CNMCA LETKF
Recent developments for CNMCA LETKF Lucio Torrisi and Francesca Marcucci CNMCA, Italian National Met Center Outline Implementation of the LETKF at CNMCA Treatment of model error in the CNMCA-LETKF The
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 informationAdvances and Challenges in Ensemblebased Data Assimilation in Meteorology. Takemasa Miyoshi
January 18, 2013, DA Workshop, Tachikawa, Japan Advances and Challenges in Ensemblebased Data Assimilation in Meteorology Takemasa Miyoshi RIKEN Advanced Institute for Computational Science Takemasa.Miyoshi@riken.jp
More informationUniv. of Maryland-College Park, Dept. of Atmos. & Oceanic Science. NOAA/NCEP/Environmental Modeling Center
The Tangent Linear Normal Mode Constraint in GSI: Applications in the NCEP GFS/GDAS Hybrid EnVar system and Future Developments Daryl Kleist 1 David Parrish 2, Catherine Thomas 1,2 1 Univ. of Maryland-College
More informationNonlinear State Estimation! Particle, Sigma-Points Filters!
Nonlinear State Estimation! Particle, Sigma-Points Filters! Robert Stengel! Optimal Control and Estimation, MAE 546! Princeton University, 2017!! Particle filter!! Sigma-Points Unscented Kalman ) filter!!
More informationCalculation and Application of MOPITT Averaging Kernels
Calculation and Application of MOPITT Averaging Kernels Merritt N. Deeter Atmospheric Chemistry Division National Center for Atmospheric Research Boulder, Colorado 80307 July, 2002 I. Introduction Retrieval
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 informationAssimilation of Satellite Infrared Brightness Temperatures and Doppler Radar Observations in a High-Resolution OSSE
Assimilation of Satellite Infrared Brightness Temperatures and Doppler Radar Observations in a High-Resolution OSSE Jason Otkin and Becky Cintineo University of Wisconsin-Madison, Cooperative Institute
More informationA data-driven method for improving the correlation estimation in serial ensemble Kalman filter
A data-driven method for improving the correlation estimation in serial ensemble Kalman filter Michèle De La Chevrotière, 1 John Harlim 2 1 Department of Mathematics, Penn State University, 2 Department
More informationA study on the spread/error relationship of the COSMO-LEPS ensemble
4 Predictability and Ensemble Methods 110 A study on the spread/error relationship of the COSMO-LEPS ensemble M. Salmi, C. Marsigli, A. Montani, T. Paccagnella ARPA-SIMC, HydroMeteoClimate Service of Emilia-Romagna,
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 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 informationUpdate on the KENDA project
Christoph Schraff Deutscher Wetterdienst, Offenbach, Germany and many colleagues from CH, D, I, ROM, RU Km-scale ENsemble-based Data Assimilation : COSMO priority project Local Ensemble Transform Kalman
More information