Objective localization of ensemble covariances: theory and applications
|
|
|
- Austen Stevenson
- 5 years ago
- Views:
Transcription
1 Institutionnel Grand Public Objective localization of ensemble covariances: theory and applications Yann Michel1, B. Me ne trier2 and T. Montmerle1 Professionnel (1) Me te o-france & CNRS, Toulouse, France (2) NCAR, Boulder, CO, USA 25 Feb Noir & blanc 4th International Symposium on Data Assimilation. RIKEN AICS, Kobe, Japan. Papeterie Y. Michel et. al., ISDA 2015 Objective localization 1
2 Introduction: Ensemble data assimilation Early days variationnal data assimilation considered modelled B estimated from NMC method. EnKF uses Monte Carlo sampling (ensembles) to better estimate covariances. We can also use ensembles in variationnal data assimilation, but there is sampling noise. Perturbed Backgrounds Cycle n (Perturbed) Observations Perturbed Backgrounds Cycle n+1 (Perturbed) Observations Cycle n+2 Perturbed Backgrounds Assimilation (Perturbed) Forecasts Assimilation (Perturbed) Forecasts Perturbed Analyses Perturbed Analyses Y. Michel et. al., ISDA 2015 Objective localization 2
3 Introduction: Localization Covariance localization through a Schur product: B = L B B ij = L ij Bij Sample correlation, localization and localized correlation for specic humidity in the boundary layer (AROME-France - 30 members) Y. Michel et. al., ISDA 2015 Objective localization 3
4 Introduction: Localization Covariance localization through a Schur product: B = L B B ij = L ij Bij Sample correlation, localization and localized correlation for specic humidity in the boundary layer (AROME-France - 30 members) Y. Michel et. al., ISDA 2015 Objective localization 3
5 Introduction: Localization Covariance localization through a Schur product: B = L B B ij = L ij Bij Sample correlation, localization and localized correlation for specic humidity in the boundary layer (AROME-France - 30 members) Y. Michel et. al., ISDA 2015 Objective localization 3
6 Introduction: Localization Covariance localization through a Schur product: B = L B B ij = L ij Bij Sample correlation, localization and localized correlation for specic humidity in the boundary layer (AROME-France - 30 members) Y. Michel et. al., ISDA 2015 Objective localization 3
7 Introduction: Localization Covariance localization through a Schur product: B = L B B ij = L ij Bij Sample correlation, localization and localized correlation for specic humidity in the boundary layer (AROME-France - 30 members) Y. Michel et. al., ISDA 2015 Objective localization 3
8 Introduction: Localization Covariance localization through a Schur product: B = L B B ij = L ij Bij Sample correlation, localization and localized correlation for specic humidity in the boundary layer (AROME-France - 30 members) Y. Michel et. al., ISDA 2015 Objective localization 3
9 Introduction: Localization Covariance localization through a Schur product: B = L B B ij = L ij Bij Sample correlation, localization and localized correlation for specic humidity in the boundary layer (AROME-France - 30 members) Y. Michel et. al., ISDA 2015 Objective localization 3
10 Introduction: Localization Covariance localization through a Schur product: B = L B B ij = L ij Bij Sample correlation, localization and localized correlation for specic humidity in the boundary layer (AROME-France - 30 members) Y. Michel et. al., ISDA 2015 Objective localization 3
11 Introduction: Localization Covariance localization through a Schur product: B = L B B ij = L ij Bij Sample correlation, localization and localized correlation for specic humidity in the boundary layer (AROME-France - 30 members) Y. Michel et. al., ISDA 2015 Objective localization 3
12 Introduction: Localization Covariance localization through a Schur product: B = L B B ij = L ij Bij Sample correlation, localization and localized correlation for specic humidity in the boundary layer (AROME-France - 30 members) Y. Michel et. al., ISDA 2015 Objective localization 3
13 Introduction: Localization Covariance localization through a Schur product: B = L B B ij = L ij Bij Sample correlation, localization and localized correlation for specic humidity in the boundary layer (AROME-France - 30 members) Y. Michel et. al., ISDA 2015 Objective localization 3
14 Introduction: Localization Covariance localization through a Schur product: B = L B B ij = L ij Bij Sample correlation, localization and localized correlation for specic humidity in the boundary layer (AROME-France - 30 members) Y. Michel et. al., ISDA 2015 Objective localization 3
15 Introduction: Localization Covariance localization through a Schur product: B = L B B ij = L ij Bij Sample correlation, localization and localized correlation for specic humidity in the boundary layer (AROME-France - 30 members) Y. Michel et. al., ISDA 2015 Objective localization 3
16 Introduction: Localization Covariance localization through a Schur product: B = L B B ij = L ij Bij Sample correlation, localization and localized correlation for specic humidity in the boundary layer (AROME-France - 30 members) Y. Michel et. al., ISDA 2015 Objective localization 3
17 Introduction Aim of this study: Development of a general theory for the linear filtering of background error sample covariances. (Application to the spatial filtering of variances). Application to the spatial localization of covariances. Illustration with the global ARPEGE model and the convective scale AROME model. Main results: Dependency on the ensemble size; Dependency on the variable; Vertical dependency of horizontal and vertical localizations. Y. Michel et. al., ISDA 2015 Objective localization 4
18 Outline 1 Introduction 2 Theory 3 Horizontal localization 4 Vertical localization 5 Conclusion and perspectives Y. Michel et. al., ISDA 2015 Objective localization 5
19 Outline 1 Introduction 2 Theory 3 Horizontal localization 4 Vertical localization 5 Conclusion and perspectives Y. Michel et. al., ISDA 2015 Objective localization 6
20 Theory: how to get started (I) The sample covariance matrix B is estimated using standard unbiased formulae: B = 1 N 1 XXT where X is the matrix of perturbations. Asymptotic convergence with ensemble size B = lim B N Sampling error is B e = B B Estimate using a linear filter is B = F B + f Filtering error is B e = B B Y. Michel et. al., ISDA 2015 Objective localization 7
21 Theory: how to get started (II) From our statistical knowledge Unbiased sampling error E[ B e ] = 0 Covariance of sampling error has known expression Cov( B ij, e B kl) e = 1 N Ξ ijkl 1 N B ij B kl 1 + N (N 1) ( B ik B jl + B il B jk) = 1 N 1 ( B ik B jl + B il B jk) (under Gaussianity) From the linear filtering theory Optimal filter has unbiased error E[ B e ] = 0 Orthogonality relationship E[ B e ij B kl ] = 0 Y. Michel et. al., ISDA 2015 Objective localization 8
22 Theory: end result L ij = (N 1)2 N(N 3) (N 1) (N+1)(N 2) N E[ Ξijij ] N 1 E[ṽ + i ṽ j ] (N 2)(N 3) E[ B2 ij ] N(N 2)(N 3) E[ B2 ij ] ( ) N 1 E[ṽ i ṽ j ] E[ B2 ij ] (under Gaussianity) ( ) (N 1) (N+1)(N 2) N 1 1 E[ C 2 ij ] (at small separation distance) and where: N is ensemble size, Ξ ijij is the sample fourth order moment, ṽ j is the sample variance at j, B 2 ij is the sample covariance squared, C 2 ij is the sample correlation squared. Y. Michel et. al., ISDA 2015 Objective localization 9
23 Ergodic assumptions E is replaced by averaging over sampling points. Sampling is randomly distributed over the domain. Binning is isotropic and based on distance. Y. Michel et. al., ISDA 2015 Objective localization 10
24 Ergodic assumptions E is replaced by averaging over sampling points. Sampling is randomly distributed over the domain. Binning is isotropic and based on distance. Y. Michel et. al., ISDA 2015 Objective localization 10
25 Ergodic assumptions E is replaced by averaging over sampling points. Sampling is randomly distributed over the domain. Binning is isotropic and based on distance. Y. Michel et. al., ISDA 2015 Objective localization 10
26 Ergodic assumptions E is replaced by averaging over sampling points. Sampling is randomly distributed over the domain. Binning is isotropic and based on distance. Y. Michel et. al., ISDA 2015 Objective localization 10
27 Outline 1 Introduction 2 Theory 3 Horizontal localization 4 Vertical localization 5 Conclusion and perspectives Y. Michel et. al., ISDA 2015 Objective localization 11
28 Results: horizontal localizations for the convective scale AROME model (I) Raw horizontal localizations are well approximated by Gaspari and Cohn s compactly supported function. Y. Michel et. al., ISDA 2015 Objective localization 12
29 Results: horizontal localizations length-scales for the convective scale AROME model (II) Wind Temperature Specific Humidity Localization length-scale increases with ensemble size. Local maxima at the top of boundary layer and at the tropopause. Y. Michel et. al., ISDA 2015 Objective localization 13
30 Outline 1 Introduction 2 Theory 3 Horizontal localization 4 Vertical localization 5 Conclusion and perspectives Y. Michel et. al., ISDA 2015 Objective localization 14
31 Results: vertical localizations Comparison between variables for the global ARPEGE model (ens. size=30) Wind Temperature Specific Humidity Y. Michel et. al., ISDA 2015 Objective localization 15
32 Results: vertical localizations Comparison between variables for the global ARPEGE model (ens. size=50) Wind Temperature Specific Humidity Y. Michel et. al., ISDA 2015 Objective localization 16
33 Results: vertical localizations Comparison between variables for the global ARPEGE model (ens. size=90) Wind Temperature Specific Humidity Y. Michel et. al., ISDA 2015 Objective localization 17
34 Results: vertical localizations Comparison between variables for the global ARPEGE model (ens. size=90) Wind Temperature Specific Humidity Mostly similar between variables except when there are negative correlations (e.g. for T) - as localization follows squared correlation. Y. Michel et. al., ISDA 2015 Objective localization 17
35 Conclusion and perspectives Theory Localization is a linear filter Optimal (in the RMSE sense) filtering gives orthogonality relationships Sampling noise on covariances has known statistical structure Optimal localization of background error covariances, depending only on quantities estimated over the ensemble. Ménétrier, B., Montmerle, T., Michel, Y. and Berre, L. (2014) Linear Filtering of Sample Covariances for Ensemble-Based Data Assimilation. Part I: Optimality Criteria and Application to Variances Filtering and Covariances Localization. Monthly Weather Review. Y. Michel et. al., ISDA 2015 Objective localization 18
36 Conclusion and perspectives Application Localization length-scale increases with ensemble size. There are differences between variables: lobes on the vertical for T, larger length-scales on the horizontal at the tropopause (T, U ) and at the top of the boundary layer (U ). We will test these localizations in the 3D/4D-EnVar schemes being developped for AROME model. Ménétrier, B., Montmerle, T., Michel, Y. and Berre, L. (2014) Linear Filtering of Sample Covariances for Ensemble-Based Data Assimilation. Part II: Application to a Convective-Scale NWP Model. Monthly Weather Review. Y. Michel et. al., ISDA 2015 Objective localization 19
Ensemble of Data Assimilations methods for the initialization of EPS
Ensemble of Data Assimilations methods for the initialization of EPS Laure RAYNAUD Météo-France ECMWF Annual Seminar Reading, 12 September 2017 Introduction Estimating the uncertainty in the initial conditions
Loïk Berre Météo-France (CNRM/GAME) Thanks to Gérald Desroziers
Estimation and diagnosis of analysis/background errors using ensemble assimilation Loïk Berre Météo-France (CNRM/GAME) Thanks to Gérald Desroziers Outline 1. Simulation of the error evolution 2. The operational
Variational ensemble DA at Météo-France Cliquez pour modifier le style des sous-titres du masque
Cliquez pour modifier le style du titre Variational ensemble DA at Météo-France Cliquez pour modifier le style des sous-titres du masque L. Berre, G. Desroziers, H. Varella, L. Raynaud, C. Labadie and
4DEnVar. 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
Representation 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,
Recent achievements in the data assimilation systems of ARPEGE and AROME-France
Recent achievements in the data assimilation systems of ARPEGE and AROME-France P. Brousseau and many colleagues from (CNRM/GMAP) 38th EWGLAM and 23 SRNWP Meeting Rome, 04 October 2016 Meteo-France NWP
EnKF 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
Inhomogeneous Background Error Modeling and Estimation over Antarctica with WRF-Var/AMPS
Inhomogeneous Background Error Modeling and Estimation over Antarctica with WRF-Var/AMPS Yann MICHEL 1 Météo-France, CNRM/GMAP 2 NCAR, MMM/DAG 10 th Annual WRF Users Workshop 23 th June 2009 Yann MICHEL
Report on the Joint SRNWP workshop on DA-EPS Bologna, March. Nils Gustafsson Alex Deckmyn.
Report on the Joint SRNWP workshop on DA-EPS Bologna, 22-24 March Nils Gustafsson Alex Deckmyn http://www.smr.arpa.emr.it/srnwp/ Purpose of the workshop On the one hand, data assimilation techniques require
A 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
Comparison of ensemble and NMC type of background error statistics for the ALADIN/HU model
Comparison of ensemble and NMC type of background error statistics for the ALADIN/HU model Kristian Horvath horvath@cirus.dhz.hr Croatian Meteorological and Hydrological Service supervised by Bölöni Gergely
Filtering of variances and correlations by local spatial averaging. Loïk Berre Météo-France
Filtering of variances and correlations by local spatial averaging Loïk Berre Météo-France Outline 1. Contrast between two extreme approaches in Var/EnKF? 2.The spatial structure of sampling noise and
Localization 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.
Environment 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
Thunderstorm-Scale EnKF Analyses Verified with Dual-Polarization, Dual-Doppler Radar Data
Thunderstorm-Scale EnKF Analyses Verified with Dual-Polarization, Dual-Doppler Radar Data David Dowell and Wiebke Deierling National Center for Atmospheric Research, Boulder, CO Ensemble Data Assimilation
Relative 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
Impact of an improved Backgrund-error covariance matrix with AROME 3Dvar system over Tunisia
Institut National de la Météorologie Joint 28th ALADIN Workshop & HIRLAM All Staff Meeting Impact of an improved Backgrund-error covariance matrix with AROME 3Dvar system over Tunisia Wafa KHALFAOUI *,
Assimilation of IASI reconstructed radiances from Principal Components in AROME model
Assimilation of IASI reconstructed radiances from Principal Components in AROME model J. Andrey, V. Guidard, N. Fourrie and J.-F. Mahfouf CNRM-GAME (Météo-France and CNRS) October 3 th, 215 Outline 1 Background
DATA 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
New developments for the use of microphysical variables for the assimilation of IASI radiances in convective scale models
New developments for the use of microphysical variables for the assimilation of IASI radiances in convective scale models Pauline Martinet a, Nadia Fourrié a, Florence Rabier a, Lydie Lavanant b, Antonia
Enhancing 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
Data 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:
Model Error in the Forecast Ensemble System at Météo-France (PEARP)
1. Model Error in the Forecast Ensemble System at Météo-France (PEARP) M. Boisserie L. Descamps P. Arbogast CNRM, Météo-France, Toulouse. 2 Introduction For a long time, it was assumed that model error
Localization 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
Reducing 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.
Comparing the SEKF with the DEnKF on a land surface model
Comparing the SEKF with the DEnKF on a land surface model David Fairbairn, Alina Barbu, Emiliano Gelati, Jean-Francois Mahfouf and Jean-Christophe Caret CNRM - Meteo France Partly funded by European Union
The ECMWF Hybrid 4D-Var and Ensemble of Data Assimilations
The Hybrid 4D-Var and Ensemble of Data Assimilations Lars Isaksen, Massimo Bonavita and Elias Holm Data Assimilation Section lars.isaksen@ecmwf.int Acknowledgements to: Mike Fisher and Marta Janiskova
The 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
Kalman 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
EnKF 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
Ensemble 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
Inter-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
A 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
1. 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,
Gaussian 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
Data assimilation in the geosciences An overview
Data assimilation in the geosciences An overview Alberto Carrassi 1, Olivier Talagrand 2, Marc Bocquet 3 (1) NERSC, Bergen, Norway (2) LMD, École Normale Supérieure, IPSL, France (3) CEREA, joint lab École
Variational 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
An Ensemble Kalman Filter for NWP based on Variational Data Assimilation: VarEnKF
An Ensemble Kalman Filter for NWP based on Variational Data Assimilation: VarEnKF Blueprints for Next-Generation Data Assimilation Systems Workshop 8-10 March 2016 Mark Buehner Data Assimilation and Satellite
Ensemble forecasting and flow-dependent estimates of initial uncertainty. Martin Leutbecher
Ensemble forecasting and flow-dependent estimates of initial uncertainty Martin Leutbecher acknowledgements: Roberto Buizza, Lars Isaksen Flow-dependent aspects of data assimilation, ECMWF 11 13 June 2007
The Evolution of Dispersion Spectra and the Evaluation of Model Differences in an Ensemble Estimation of Error Statistics for a Limited-Area Analysis
3456 M O N T H L Y W E A T H E R R E V I E W VOLUME 134 The Evolution of Dispersion Spectra and the Evaluation of Model Differences in an Ensemble Estimation of Error Statistics for a Limited-Area Analysis
Numerical 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
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)
Improved structure functions for 3D VAR
Improved structure functions for 3D VAR HIRLAM All Staff Meeting Sofia, Bulgaria, 15 18 May, 2006 Magnus Lindskog, Nils Gustafsson, Martin Ridal och Per Dahlgren Swedish Meteorological and Hydrological
Background error modelling: climatological flow-dependence
Background error modelling: climatological flow-dependence Yann MICHEL NCAR/MMM/B Meeting 16 th April 2009 1 Introduction 2 A new estimate of lengthscales 3 Climatological flow-dependence Yann MICHEL B
Estimating 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
State Space Models for Wind Forecast Correction
for Wind Forecast Correction Valérie 1 Pierre Ailliot 2 Anne Cuzol 1 1 Université de Bretagne Sud 2 Université de Brest MAS - 2008/28/08 Outline 1 2 Linear Model : an adaptive bias correction Non Linear
DATA ASSIMILATION OF VIS / NIR REFLECTANCE INTO THE DETAILED SNOWPACK MODEL SURFEX/ISBA- CROCUS
DATA ASSIMILATION OF VIS / NIR REFLECTANCE INTO THE DETAILED SNOWPACK MODEL SURFEX/ISBA- CROCUS Luc CHARROIS 1,2 PhD 2013 / 2016 E. Cosme 2, M. Dumont 1, M. Lafaysse 1, S. Morin 1, Q. Libois 2, G. Picard
T2.3: Use of ensemble information in ocean analysis and development of efficient 4D-Var
T2.3: Use of ensemble information in ocean analysis and development of efficient 4D-Var A. Weaver 1,B.Ménétrier 1,J.Tshimanga 1 and A. Vidard 2 1 CERFACS, Toulouse 2 INRIA/LJK, Grenoble December 10, 2015
Modelling of background error covariances for the analysis of clouds and precipitation
Modelling of background error covariances for the analysis of clouds and precipitation Thibaut Montmerle, Yann Michel and Benjamin Ménétrier Météo-France/CNRM-GAME 42 av. G. Coriolis, 31057 Toulouse, France
New Applications and Challenges In Data Assimilation
New Applications and Challenges In Data Assimilation Met Office Nancy Nichols University of Reading 1. Observation Part Errors 1. Applications Coupled Ocean-Atmosphere Ensemble covariances for coupled
Improved 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
Ensemble-variational assimilation with NEMOVAR Part 2: experiments with the ECMWF system
Ensemble-variational assimilation with NEMOVAR Part 2: experiments with the ECMWF system Toulouse, 20/06/2017 Marcin Chrust 1, Hao Zuo 1 and Anthony Weaver 2 1 ECMWF, UK 2 CERFACS, FR Marcin.chrust@ecmwf.int
Simulation of error cycling
Simulation of error cycling Loïk BERRE, Météo-France/CNRS ISDA, Reading, 21 July 2016 with inputs from R. El Ouaraini, L. Raynaud, G. Desroziers, C. Fischer Motivations and questions EDA and innovations
Review 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
Comparing Variational, Ensemble-based and Hybrid Data Assimilations at Regional Scales
Comparing Variational, Ensemble-based and Hybrid Data Assimilations at Regional Scales Meng Zhang and Fuqing Zhang Penn State University Xiang-Yu Huang and Xin Zhang NCAR 4 th EnDA Workshop, Albany, NY
A Comparison between the 4DVAR and the Ensemble Kalman Filter Techniques for Radar Data Assimilation
NOVEMBER 2005 C A Y A E T A L. 3081 A Comparison between the 4DVAR and the Ensemble Kalman Filter Techniques for Radar Data Assimilation A. CAYA, J.SUN, AND C. SNYDER National Center for Atmospheric Research,*
Ensemble 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
Impact of IASI assimilation in convective scale model AROME
Impact of IASI assimilation in convective scale model AROME Vincent GUIDARD Pierre BROUSSEAU, Nadia FOURRIE, Florence RABIER Météo-France and CNRS / CNRM-GAME 2 Overview 1. What is AROME? 2. IASI in next
A general hybrid formulation of the background-error covariance matrix for ensemble-variational ocean data assimilation
.. A general hybrid formulation of the background-error covariance matrix for ensemble-variational ocean data assimilation Anthony Weaver Andrea Piacentini, Serge Gratton, Selime Gürol, Jean Tshimanga
Generating 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
P 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
ICON. Limited-area mode (ICON-LAM) and updated verification results. Günther Zängl, on behalf of the ICON development team
ICON Limited-area mode (ICON-LAM) and updated verification results Günther Zängl, on behalf of the ICON development team COSMO General Meeting, Offenbach, 07.09.2016 Outline Status of limited-area-mode
Bayesian Inverse problem, Data assimilation and Localization
Bayesian Inverse problem, Data assimilation and Localization Xin T Tong National University of Singapore ICIP, Singapore 2018 X.Tong Localization 1 / 37 Content What is Bayesian inverse problem? What is
Particle Filters. Outline
Particle Filters M. Sami Fadali Professor of EE University of Nevada Outline Monte Carlo integration. Particle filter. Importance sampling. Degeneracy Resampling Example. 1 2 Monte Carlo Integration Numerical
Tropical Pacific Ocean model error covariances from Monte Carlo simulations
Q. J. R. Meteorol. Soc. (2005), 131, pp. 3643 3658 doi: 10.1256/qj.05.113 Tropical Pacific Ocean model error covariances from Monte Carlo simulations By O. ALVES 1 and C. ROBERT 2 1 BMRC, Melbourne, Australia
Recent experience at Météo-France on the assimilation of observations at high temporal frequency Cliquez pour modifier le style du titre
Recent experience at Météo-France on the assimilation of observations at high temporal frequency Cliquez pour modifier le style du titre J.-F. Mahfouf, P. Brousseau, P. Chambon and G. Desroziers Météo-France/CNRS
What s new for ARPEGE/ALADIN since Utrecht 2009?
What s new for ARPEGE/ALADIN since Utrecht 2009? E. Bazile with contributions from GMAP/OBS and GMAP/PROC ASM HIRLAM ALADIN 12-16 April, 2010 Krakow Outline QPF performance of ARPEGE/ALADIN with the new
Data assimilation with sea ice models
SIPN workshop NCAR, Boulder 1-2 April 2014 Data assimilation with sea ice models François Massonnet brine content mixed layer heat content heat content concentration thickness brine content mixed layer
University of Oklahoma, Norman Oklahoma Monthly Weather Review Submitted June 2017 Accepted November 2017
Development of a Hybrid En3DVar Data Assimilation System and Comparisons with 3DVar and EnKF for Radar Data Assimilation with Observing System Simulation Experiments Rong Kong 1,2, Ming Xue 1,2, and Chengsi
Snow-atmosphere interactions at Dome C, Antarctica
Snow-atmosphere interactions at Dome C, Antarctica Éric Brun, Vincent Vionnet CNRM/GAME (Météo-France and CNRS) Christophe Genthon, Delphine Six, Ghislain Picard LGGE (CNRS and UJF)... and many colleagues
Convective-scale data assimilation in the Weather Research and Forecasting model using a nonlinear ensemble filter
Convective-scale data assimilation in the Weather Research and Forecasting model using a nonlinear ensemble filter Jon Poterjoy, Ryan Sobash, and Jeffrey Anderson National Center for Atmospheric Research
Changes in the Arpège 4D-VAR and LAM 3D-VAR. C. Fischer With contributions by P. Brousseau, G. Kerdraon, J.-F. Mahfouf, T.
Changes in the Arpège 4D-VAR and LAM 3D-VAR C. Fischer With contributions by P. Brousseau, G. Kerdraon, J.-F. Mahfouf, T. Montmerle Content Arpège 4D-VAR Arome-France Other applications: Aladin Overseas,
Recent 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
Assimilation des Observations et Traitement des Incertitudes en Météorologie
Assimilation des Observations et Traitement des Incertitudes en Météorologie Olivier Talagrand Laboratoire de Météorologie Dynamique, Paris 4èmes Rencontres Météo/MathAppli Météo-France, Toulouse, 25 Mars
The Big Leap: Replacing 4D-Var with 4D-EnVar and life ever since
The Big Leap: Replacing 4D-Var with 4D-EnVar and life ever since Symposium: 20 years of 4D-Var at ECMWF 26 January 2018 Mark Buehner 1, Jean-Francois Caron 1 and Ping Du 2 1 Data Assimilation and Satellite
Convergence of Square Root Ensemble Kalman Filters in the Large Ensemble Limit
Convergence of Square Root Ensemble Kalman Filters in the Large Ensemble Limit Evan Kwiatkowski, Jan Mandel University of Colorado Denver December 11, 2014 OUTLINE 2 Data Assimilation Bayesian Estimation
Recent progress in convective scale Arome NWP system and on-going research activities
Recent progress in convective scale Arome NWP system and on-going research activities P. Brousseau, P. Chambon, G. Faure, R. Honnert, A. Mary, N. Merlet, Y. Seity, B. Vié, E. Wattrelot (presented by F.
Radar data assimilation using a modular programming approach with the Ensemble Kalman Filter: preliminary results
Radar data assimilation using a modular programming approach with the Ensemble Kalman Filter: preliminary results I. Maiello 1, L. Delle Monache 2, G. Romine 2, E. Picciotti 3, F.S. Marzano 4, R. Ferretti
Heterogeneous Convective-Scale Background Error Covariances with the Inclusion of Hydrometeor Variables
2994 M O N T H L Y W E A T H E R R E V I E W VOLUME 139 Heterogeneous Convective-Scale Background Error Covariances with the Inclusion of Hydrometeor Variables YANN MICHEL Météo-France, CNRM-GAME and CNRS,
Summary of activities with SURFEX data assimilation at Météo-France. Jean-François MAHFOUF CNRM/GMAP/OBS
Summary of activities with SURFEX data assimilation at Météo-France Jean-François MAHFOUF CNRM/GMAP/OBS Outline Status at Météo-France in 2008 Developments undertaken during 2009-2011 : Extended Kalman
The 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
4. 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
On Ensemble and Particle Filters for Numerical Weather Prediction
On Ensemble and Particle Filters for Numerical Weather Prediction Roland Potthast Hendrik Reich, Christoph Schraff Andreas Rhodin, Ana Fernandez Alexander Cress, Dora Foring Annika Schomburg, Africa Perianez
Nonlinear Balance at Mesoscale: Balanced vertical motion with respect to model convection temperature tendencies
Motivations and Methodology Results Summary Nonlinear Balance at Mesoscale: Balanced vertical motion with respect to model convection temperature tendencies Christian Pagé Peter Zwack page.christian@uqam.ca
Convergence of the Ensemble Kalman Filter in Hilbert Space
Convergence of the Ensemble Kalman Filter in Hilbert Space Jan Mandel Center for Computational Mathematics Department of Mathematical and Statistical Sciences University of Colorado Denver Parts based
AROME-EPS development in Météo-France
AROME-EPS development in Météo-France Int.Conf.Ens.Methods Toulouse, 15 Nov 2012 francois.bouttier@meteo.fr collaborators: Olivier Nuissier, Laure Raynaud, Benoît Vié Ensemble setup Boundary conditions
Advances 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
Meteo-France operational land surface analysis for NWP: current status and perspectives
Meteo-France operational land surface analysis for NWP: current status and perspectives Camille Birman, Jean-François Mahfouf, Yann Seity, Eric Bazile, Françoise Taillefer, Zied Sassi, Nadia Fourrié, Vincent
State of the art of wind forecasting and planned improvements for NWP Helmut Frank (DWD), Malte Mülller (met.no), Clive Wilson (UKMO)
State of the art of wind forecasting and planned improvements for NWP Helmut Frank (DWD), Malte Mülller (met.no), Clive Wilson (UKMO) thanks to S. Bauernschubert, U. Blahak, S. Declair, A. Röpnack, C.
ELEG 3143 Probability & Stochastic Process Ch. 6 Stochastic Process
Department of Electrical Engineering University of Arkansas ELEG 3143 Probability & Stochastic Process Ch. 6 Stochastic Process Dr. Jingxian Wu wuj@uark.edu OUTLINE 2 Definition of stochastic process (random
A Comparison between the 3/4DVAR and Hybrid Ensemble-VAR Techniques for Radar Data Assimilation ABSTRACT
36 th AMS CONFERENCE ON RADAR METEOROLOGY, 16-2 SEPTEMBER 213, BRECKENRIDGE, COLORADO A Comparison between the 3/4DVAR and Hybrid Ensemble-VAR Techniques for Radar Data Assimilation Hongli Wang *1, Xiang-Yu
Ensemble aerosol forecasts and assimila1on at ECMWF
Ensemble aerosol forecasts and assimila1on at ECMWF Angela Benede*, Miha Razinger, Luke Jones & Jean- Jacques Morcre
Adaptive Localization: Proposals for a high-resolution multivariate system Ross Bannister, HRAA, December 2008, January 2009 Version 3.
Adaptive Localization: Proposals for a high-resolution multivariate system Ross Bannister, HRAA, December 2008, January 2009 Version 3.. The implicit Schur product 2. The Bishop method for adaptive localization
How 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
Estimating 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 Rahul Mahajan NOAA / NCEP / EMC - IMSG College Park, MD 20740, USA Acknowledgements: Ron Gelaro, Ricardo
Modal view of atmospheric predictability
Modal view of atmospheric predictability Nedjeljka Žagar University of Ljubljana, Ljubljana, Slovenia Based on Žagar, N., R. Buizza and J. Tribbia, J. Atmos. Sci., 0, and Žagar, N., J. Anderson, N. Collins,
A Note on the Particle Filter with Posterior Gaussian Resampling
Tellus (6), 8A, 46 46 Copyright C Blackwell Munksgaard, 6 Printed in Singapore. All rights reserved TELLUS A Note on the Particle Filter with Posterior Gaussian Resampling By X. XIONG 1,I.M.NAVON 1,2 and
Performance of ensemble Kalman filters with small ensembles
Performance of ensemble Kalman filters with small ensembles Xin T Tong Joint work with Andrew J. Majda National University of Singapore Sunday 28 th May, 2017 X.Tong EnKF performance 1 / 31 Content Ensemble
Can 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
A Non-Gaussian Ensemble Filter Update for Data Assimilation
4186 M O N T H L Y W E A T H E R R E V I E W VOLUME 138 A Non-Gaussian Ensemble Filter Update for Data Assimilation JEFFREY L. ANDERSON NCAR Data Assimilation Research Section,* Boulder, Colorado (Manuscript