Simulation of error cycling
|
|
- Cody Walton
- 6 years ago
- Views:
Transcription
1 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
2 Motivations and questions EDA and innovations for diagnosing contributions (background errors, observation errors, model errors) in error cycling over one week. Revisit formalism & some previous EDA experiments in the litterature : Respective global amplitudes of these 3 error sources? Evolution of error contributions during the cycling?
3 Contents EDA context at Météo-France Quasi-linear expansion of forecast errors Old background errors / Recent observation errors Observation errors / Model errors
4 Contents EDA context at Météo-France Quasi-linear expansion of forecast errors Old background errors / Recent observation errors Observation errors / Model errors
5 Ensemble of perturbed Data Assimilations (EDA) : simulation of error cycling ε b ε f = Mε a + ε m ε a = (I-KH)ε b + Kε o ε o t i -6h t i t i +6h (e.g. Houtekamer et al 1996, Fisher 2003, Berre et al 2006)
6 Global operational EDA at Météo-France 25 members, T479 (40 km) L105, Arpege 4D-Var (minim T149), 6h cycle. 4D-Var analysis perturbations : observation perturbations (drawn from R, incl. spatial corr. for AMVs), background perturbations (evolved analysis pertbs and model pertbs). Multiplicative inflation of forecast perturbations, using innovation-based diagnostics. Spatially filtered variances for observation QC and minimisation, wavelet-filtered 3D correlations. An EDA is also being developed at mesoscale (AROME, oper 2018).
7 Dynamics of background error variances with EDA Standard deviations of surface pressure (hpa) (2/2/2010) (superimposed with MSLP analysis, in hpa). (e.g. Berre et al 2007, Raynaud et al 2012)
8 Dynamics of horizontal correlations from ensemble and wavelet filtering Length scales (km) for wind near 500 hpa (28/2/2010), superimposed with geopotential (Berre, Varella et Desroziers 2015)
9 Titre de la diapo Dynamics of vertical correlations from ensemble and wavelet filtering Vertical correlations of temperature between 850 & 870 hpa (28/2/2010) (Berre, Varella and Desroziers 2015)
10 Contents EDA context at Météo-France Quasi-linear expansion of forecast errors Old background errors / Recent observation errors Observation errors / Model errors
11 Expansion of forecast error contributions (quasi-linear framework) At a given step t 0 of the cycling : ε a 0 = (I-K 0 H 0 ) ε b 0 + K 0 εo 0 ε f 0 = M 0 εa 0 + εm 0 = M 0 (I-K 0 H 0 ) ε b 0 + M 0 K 0 εo 0 + εm 0 At a later step t i (e.g. one week later) : ε f i = T i+1 ε b 0 + Σ T i-j M j K j εo j + Σ T i-j ε m j j 0 j 0 where T i-j = Π M k (I-K k H k ) = cycling operator k (over k successive analysis/forecast steps from t j to t i ). (El Ouaraini and Berre 2011)
12 Age of error contributions Error contributions over one week, from t 0 to t i : OLD CYCLES t 0 RECENT CYCLES t j t i date of cycle (t) ε b 0 old background errors ε o j εm j recent obs. errors recent model errors ε f i Full expansion of forecast errors at t i : ε f i = T i+1 ε b 0 + Σ T i-j M j K j εo j + Σ T i-j ε m j where T i-j = Π M k (I-K k H k ) = cycling operator k j 0 j 0 (over k successive analysis/forecast steps from t j to t i ). (El Ouaraini and Berre 2011)
13 Several processes in error cycling Expansion of forecast errors : ε f i = T i+1 ε b 0 + Σ T i-j M j K j εo j + Σ T i-j ε m j where T i-j = Π M k (I-K k H k ) = cycling operator. k j 0 j 0 Old background errors ε b 0 : repeted analysis damping & model propagation. Recent observation errors ε o j : filtering (K) & propagation (M) ; damping & propagation (T) ; accumulation (Σ). Recent model errors ε m j : damping & propagation ; accumulation.
14 Links with EDA and innovations Same equation for errors and perturbations : ε f i = T i+1 ε b 0 + Σ T i-j M j K j εo j + Σ T i-j ε m j j 0 j 0 Simulation of error cycling : cycle observation and model perturbations added to deterministic system (e.g. with 4D-Var and non linear forecasts included). Observation and model perturbations are often (+/- indirectly) derived from innovation-based estimates of covariances. Amplitudes of some estimated (accumulated) error contributions can be diagnosed and compared using EDA and innovations.
15 Experimental framework Revisit some sensitivity experiments with a previous EDA configuration : no model perturbations (offline tuning of horizontally averaged variances) ; same K in all compared configurations (static spectral correlations, flow-dependent variances). Comparison «cold start» / «warm start» ensembles : diagnosis / evolution of old background perturbations. EDA configurations without (recent) model perturbations : to diagnose the 2 other contributions (old background errors, recent observation errors). to estimate model error variance by comparison with innovation-based estimates of forecast errors. Offline diagnosis of global amplitudes (variances) of perturbations.
16 Contents EDA context at Météo-France Quasi-linear expansion of forecast errors Old background errors / Recent observation errors Observation errors / Model errors
17 Contribution of recent observation errors ε f i = T i+1 ε b 0 + Σ T i-j M j K j εo j + Σ T i-j ε m j j 0 j 0 Contribution of recent observation errors to forecast errors : ε fo i = Σ T i-j M j K j ε o j j 0 Simulation and evolution using EDA : * use ε b 0 =0 (cold start = start from unperturbed background), then cycle with ε o j ( j 0 ) (and with εm j = 0). * compute evolution of EDA spread : σ( ε fo t = Σ T t-j M t K t ε o t ).
18 Accumulation of recent observation errors Evolution of σ( ε fo t = Σ T t-j M j K j ε o j ) j 0 at successive analysis steps (t 0, t 1, t i ) σ(ε fo t ) (vorticity (10-5 s -1 ), 500 hpa) (El Ouaraini and Berre 2011, see also Fisher et al 2005) Observation error contribution increases (accumulation), then it stabilizes after 3-5 days of cycling, due to analysis damping : T t-j = Π M k (I-K k H k ). k
19 Accumulation of recent observation errors Evolution of ε fo t = Σ T t-j M j K j ε o j j 0 at successive analysis steps (t 0, t 1, t 2, t i ) : ε fo 0 = M 0 K 0 εo 0 ε fo 1 = M 1 K 1 εo 1 + T 1 εfo 0 ε fo 2 = M 2 K 2 εo 2 + T 1 εfo 1 + T 2 εfo 0. ε fo t = M t K t ε o t + T 1 ε fo t-1 + T 2 εfo t T 12 εfo t-12 + Σ T k εfo t-k k>12 T is mainly damping, so old obs. error contrib. (with time-lag k>12) become negligible, and σ(ε fo t ) stabilizes ~ at step t = 12 (3 days).
20 Contribution of old background errors ε f i = T i+1 ε b 0 + Σ T i-j M j K j εo j + Σ T i-j ε m j Contribution of «old» background errors : ε fb i = T i+1 εb 0 with T i+1 = Π M k (I-K k H k ) Simulation and evolution using EDA : use warm start = ensemble started 6 days before t 0 : * use ε b 0 =Σ T 0-j M j K j εo j (= contribution of old observation errors), j <0 j 0 i k=0 j 0 then cycle with ε o j ( j 0 ) (and with εm j = 0). * compute evolution of sqrt of σ²(ε fb t ) = σ²(εfb t + εfo t ) σ²(εfo t ) (~sqrt of EDA variance difference between warm start and cold start)
21 Evolution of old background error contribution σ(ε fb t = T t+1 εb 0 ) t with T t+1 = Π M k (I-K k H k ) k=0 σ(ε fo t ) (vorticity (10-5 s -1 ), 500 hpa) ε b 0 = Σ T 0-j M j K j εo j (contribution of old observation errors) j <0 σ(ε fb 0 ) > σ(εfo 0 ) (for vorticity). σ(ε fb i = T i+1 εb 0 ) decreases ~ linearly / analysis damping.
22 Contribution of old+recent observation errors σ(ε fb t + εfo t ) σ(ε fb t = T t+1 εb 0 ) σ(ε fo t ) (vorticity (10-5 s -1 ), 500 hpa) Total contribution of (old+recent) obs errors is stable : damping of old observation errors is compensated by accumulation of recent obs errors.
23 Dependence on amplitude of old background errors σ(ε fb t = T t+1 εb 0 ) (vorticity (10-5 s -1 ), 500 hpa) ε b 0 = Σ T 0-j M j K j εo j (old observation errors) j <0 ε b 0 = B1/2 η (old model errors + old observation errors) σ(ε fb t = T t+1 εb 0 ) decreases ~ linearly / analysis damping.
24 First conclusions / formalism + experiments Quasi-linear expansion of forecast errors / contributions of old background errors, recent observation errors, recent model errors. Compare cold/warm start EDA (without recent model perturbations) to diagnose error contrib. of old background and recent observations. Old background errors vanish (~linearly) after 3-5 days of cycling. Compensated by accumulation of recent observation errors. This may help for interpretation/diagnosis of total forecast errors (~ recent observation & model errors).
25 Contents EDA context at Météo-France Quasi-linear expansion of forecast errors Old background errors / Recent observation errors Observation errors / Model errors
26 Innovation covariances background error covariance E[ (y Hx b )(y Hx b ) T ] = R + HBH T separation
27 Covariances of analysis residuals y δy = y H(x b ) (innovation) δy H δx = H(x b ) H(x a ) (increment) ε o E[H δx δy T ] = HB*H T δx x a ε a E[(y H(x a )) δy T ] = R* x b ε b x* (true state) (Desroziers et al 2005)
28 Estimation of model error contributions (to forecast errors) t 0 t j t i date of cycle ε b 0 ε o j εm j ε f i Error ε f i of forecast issued from cycle t i : ε f i = T i+1 ε b 0 + Σ T i-j M j K j εo j + Σ T i-j ε m j [ T i-j = Π M k (I-K k H k ) ] j 0 j 0 k obs. errors model errors accumulated from t 0 to t i Compare σ δy ( εf i = ε fb0 i + Σ T i-j M j K j εo j + Σ T i-j ε m j ) (innovation-based) with σ EDA ( ε fb0 i + Σ T i-j M j K j εo j ) (using unperturbed model from t 0 to t i ).
29 Total forecast error versus observation error contribution σ δy (εf i = ε fb0 i +Σ T i-j M j K j εo j + Σ T i-j εm j ) σ EDA ( ε fb0 i +Σ T i-j M j K j εo j ) Pressure (hpa) Standard deviation of forecast errors (aircraft observations of temperature (K)) (Raynaud et al 2012)
30 Total forecast error versus observation error contribution σ 2 δy ( εfb0 i +Σ T i-j M j K j εo j + Σ T i-j εm j ) ~ σ2 ( Σ T i-j M j K j ε o j ) + σ2 ( Σ T i-j ε m j ) If cov(ε fb0 i, εfm i ) ~ 0, i.e. if the one-week period is large enough to neglect time correlations between old background errors & recently accumulated model errors, in addition to : σ 2 EDA ( εfb0 i +Σ T i-j M j K j εo j ) = σ2 EDA ( Σ T i-j M j K j εo j ) if σ2 (ε fb0 i ) ~ 0
31 Quantification of model error accumulated during cycling σ δy/eda ( εfm i = Σ T i-j ε m j ) σ EDA ( ε fo i = Σ T i-j M j K j ε o j ) Pressure (hpa) Standard deviation of error contributions aircraft observations of temperature (K))
32 Possible guidance for model error representation? Variance of model error accumulated over 3-5 days can be diagnosed : model error contrib. ~ one half of 6h forecast error variance ; obs. error contrib. ~ other half. Variance of model error accumulated over 6h: similar formalism, but assumptions/diagnostics on temporal correlations required. Importance of innovation-based estimates (R and B). Possible model error representations (Q 1/2 η, SPPT, SKEB, etc) may be compared/adjusted with such diagnostics.
33 Conclusions EDA for estimation of flow-dependent covariances, combined with spatial filtering methods. EDA (and innovations) for diagnosing contributions to error cycling : obs. errors contribute ~ within the last 3-5 days of DA cycling. old background errors vanish ~ after 3-5 days of DA cycling. model error contrib. are similar in amplitude to obs. error contrib. LAM : contribution of LBC errors over 1/3 of ALADIN-France domain due to advection in the cycling (El Ouaraini et al 2015). Extension to diagnostics of spatial structures, pursue comparison EDA vs innovations, etc.
34 Thank you for your attention
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
More informationLoï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
More informationThe ECMWF Hybrid 4D-Var and Ensemble of Data Assimilations
The Hybrid 4D-Var and Ensemble of Data Assimilations Lars Isaksen, Massimo Bonavita and Elias Holm Data Assimilation Section lars.isaksen@ecmwf.int Acknowledgements to: Mike Fisher and Marta Janiskova
More informationEnsemble 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
More informationFiltering 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
More informationRecent 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
More informationRecent 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
More informationReport 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
More informationAROME-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
More informationData Assimilation at CHMI
Data Assimilation at CHMI progress report from September 2008 Alena Trojáková Czech Hydrometeorological Institute Na Šabatce 17, 143 06 Prague 4 Contents 1 Introduction 3 2 Evaluation of the assimilation
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 informationDiagnosis of observation, background and analysis-error statistics in observation space
Q. J. R. Meteorol. Soc. (2005), 131, pp. 3385 3396 doi: 10.1256/qj.05.108 Diagnosis of observation, background and analysis-error statistics in observation space By G. DESROZIERS, L. BERRE, B. CHAPNIK
More informationModel error and parameter estimation
Model error and parameter estimation Chiara Piccolo and Mike Cullen ECMWF Annual Seminar, 11 September 2018 Summary The application of interest is atmospheric data assimilation focus on EDA; A good ensemble
More informationNumerical Weather prediction at the European Centre for Medium-Range Weather Forecasts (2)
Numerical Weather prediction at the European Centre for Medium-Range Weather Forecasts (2) Time series curves 500hPa geopotential Correlation coefficent of forecast anomaly N Hemisphere Lat 20.0 to 90.0
More informationImpact 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 *,
More informationEnsemble Kalman Filter based snow data assimilation
Ensemble Kalman Filter based snow data assimilation (just some ideas) FMI, Sodankylä, 4 August 2011 Jelena Bojarova Sequential update problem Non-linear state space problem Tangent-linear state space problem
More informationAssimilation of SEVIRI data II
Assimilation of SEVIRI data II report from LACE stay in Budapest, 23/1/07-17/2/07 Scientific supervisors: Roger Randriamampianina Gergely Bölöni Alena Trojáková Czech Hydrometeorological Institute Na Šabatce
More informationObservation error specifications
Observation error specifications Gérald Desroziers, with many contributions Météo-France, CNRS 4 av. G. Coriolis, 357 Toulouse Cédex, France gerald.desroziers@meteo.fr ABSTRACT The aim of this paper is
More informationLAM-EPS activities in LACE
LAM-EPS activities in LACE Martin Belluš with contributions of A. Trojáková, M. Szűcs, Ch. Wittmann, F. Weidle, Y. Wang, C. Wastl, S. Taşcu, J. Cedilnik and E. Keresturi AROME-EPS STOCHASTIC ALADIN-LAEF
More informationComparing 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
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 informationChanges 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,
More informationON DIAGNOSING OBSERVATION ERROR STATISTICS WITH LOCAL ENSEMBLE DATA ASSIMILATION
ON DIAGNOSING OBSERVATION ERROR STATISTICS WITH LOCAL ENSEMBLE DATA ASSIMILATION J. A. Waller, S. L. Dance, N. K. Nichols University of Reading 1 INTRODUCTION INTRODUCTION Motivation Only % of observations
More informationImproved 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
More informationDoppler radial wind spatially correlated observation error: operational implementation and initial results
Doppler radial wind spatially correlated observation error: operational implementation and initial results D. Simonin, J. Waller, G. Kelly, S. Ballard,, S. Dance, N. Nichols (Met Office, University of
More informationRepresenting model error in ensemble data assimilation
Nonlin. Processes Geophys., 21, 971 985, 214 www.nonlin-processes-geophys.net/21/971/214/ doi:1.5194/npg-21-971-214 Author(s) 214. CC Attribution 3. License. Representing model error in ensemble 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 informationBackground Error Covariance Modelling
Background Error Covariance Modelling Mike Fisher Slide 1 Outline Diagnosing the Statistics of Background Error using Ensembles of Analyses Modelling the Statistics in Spectral Space - Relaxing constraints
More 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 informationData assimilation in mesoscale modeling and numerical weather prediction Nils Gustafsson
Data assimilation in mesoscale modeling and numerical weather prediction Nils Gustafsson Croatian USA Workshop on Mesometeorology June 2012 Perspective: What are the important issues for development of
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 informationAn 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
More informationHybrid variational-ensemble data assimilation. Daryl T. Kleist. Kayo Ide, Dave Parrish, John Derber, Jeff Whitaker
Hybrid variational-ensemble data assimilation Daryl T. Kleist Kayo Ide, Dave Parrish, John Derber, Jeff Whitaker Weather and Chaos Group Meeting 07 March 20 Variational Data Assimilation J Var J 2 2 T
More informationThe 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
More informationThe satellite winds in the operational NWP system at Météo-France
The satellite winds in the operational NWP system at Météo-France Christophe Payan CNRM UMR 3589, Météo-France/CNRS 13th International Winds Workshop, Monterey, USA, 27 June 1st July 2016 Outline Operational
More informationThe Role of Diagnostics in Numerical Weather Prediction
The Role of Diagnostics in Numerical Weather Prediction Mark Rodwell (with the support of ECMWF colleagues and external collaborators) Using ECMWF s Forecasts (UEF2016) 8 June, ECMWF Reading 1 Ensemble
More informationDA experiences and plans in Tunisia
ALADIN Data Assimilation basic kit Working Days DA experiences and plans in Tunisia Anis Satouri Lisbon 22-03-2017 Models Installed on the HP Server ]]] ALADIN AROME HARMONIE Model version CY38t1_bf.03
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 informationDiagnostics of the prediction and maintenance of Euro-Atlantic blocking
Diagnostics of the prediction and maintenance of Euro-Atlantic blocking Mark Rodwell, Laura Ferranti, Linus Magnusson Workshop on Atmospheric Blocking 6-8 April 2016, University of Reading European Centre
More informationSummary 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
More informationERA-CLIM: Developing reanalyses of the coupled climate system
ERA-CLIM: Developing reanalyses of the coupled climate system Dick Dee Acknowledgements: Reanalysis team and many others at ECMWF, ERA-CLIM project partners at Met Office, Météo France, EUMETSAT, Un. Bern,
More informationLAM-EPS activities in RC LACE
LAM-EPS activities in RC LACE Martin Belluš with contributions of M. Imrišek, Ch. Wittmann, C. Wastl, I. O. Plenković, M. Szűcs and E. Keresturi RC LACE Predictability Area - two main subjects: AROME-EPS
More informationUse of analysis ensembles in estimating flow-dependent background error variance
Use of analysis ensembles in estimating flow-dependent background error variance Lars Isaksen, Mike Fisher and Judith Berner European Centre for Medium-range Weather Forecasts Shinfield Park, Reading,
More informationSensitivity analysis in variational data assimilation and applications
Sensitivity analysis in variational data assimilation and applications Dacian N. Daescu Portland State University, P.O. Box 751, Portland, Oregon 977-751, U.S.A. daescu@pdx.edu ABSTRACT Mathematical aspects
More information4DEnVar: link with 4D state formulation of variational assimilation and different possible implementations
QuarterlyJournalof theoyalmeteorologicalsociety Q J Meteorol Soc 4: 97 October 4 A DOI:/qj35 4DEnVar: lin with 4D state formulation of variational assimilation and different possible implementations Gérald
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 informationThe Canadian approach to ensemble prediction
The Canadian approach to ensemble prediction ECMWF 2017 Annual seminar: Ensemble prediction : past, present and future. Pieter Houtekamer Montreal, Canada Overview. The Canadian approach. What are the
More 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 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 informationUse of randomization to diagnose the impact of observations on analyses and forecasts
Q. J. R. Meteorol. Soc. (2005), 131, pp. 2821 2837 doi: 10.1256/qj.04.151 Use of randomization to diagnose the impact of observations on analyses and forecasts By GÉRALD DESROZIERS, PIERRE BROUSSEAU and
More informationDevelopment and research of GSI based hybrid EnKF Var data assimilation for HWRF to improve hurricane prediction
Development and research of GSI based hybrid EnKF Var data assimilation for HWRF to improve hurricane prediction Xuguang Wang, Xu Lu, Yongzuo Li School of Meteorology University of Oklahoma, Norman, OK,
More informationNonlinear 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
More informationQuantifying observation error correlations in remotely sensed data
Quantifying observation error correlations in remotely sensed data Conference or Workshop Item Published Version Presentation slides Stewart, L., Cameron, J., Dance, S. L., English, S., Eyre, J. and Nichols,
More informationLAM-EPS activities in RC LACE
LAM-EPS activities in RC LACE Martin Belluš with contributions of M. Szűcs, M. Dian, Ch. Wittmann, F. Weidle, Y. Wang, C. Wastl, S. Taşcu, R. Pomaga and E. Keresturi 1 Overview of activities (since last
More informationHow 4DVAR can benefit from or contribute to EnKF (a 4DVAR perspective)
How 4DVAR can benefit from or contribute to EnKF (a 4DVAR perspective) Dale Barker WWRP/THORPEX Workshop on 4D-Var and Ensemble Kalman Filter Intercomparisons Sociedad Cientifica Argentina, Buenos Aires,
More informationEvolution of Forecast Error Covariances in 4D-Var and ETKF methods
Evolution of Forecast Error Covariances in 4D-Var and ETKF methods Chiara Piccolo Met Office Exeter, United Kingdom chiara.piccolo@metoffice.gov.uk Introduction Estimates of forecast error covariances
More 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 informationComparison 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
More informationEnsemble of Data Assimilations and Uncertainty Estimation
Ensemble of Data Assimilations and Uncertainty Estimation Massimo Bonavita ECMWF, Reading, UK Massimo.Bonavita@ecmwf.int ABSTRACT The background error covariance matrix (B) plays a fundamental role in
More informationApplication and verification of ECMWF products 2012
Application and verification of ECMWF products 2012 Instituto Português do Mar e da Atmosfera, I.P. (IPMA) 1. Summary of major highlights ECMWF products are used as the main source of data for operational
More informationThe 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
More informationQuantifying observation error correlations in remotely sensed data
Quantifying observation error correlations in remotely sensed data Conference or Workshop Item Published Version Presentation slides Stewart, L., Cameron, J., Dance, S. L., English, S., Eyre, J. and Nichols,
More informationObservation impact on data assimilation with dynamic background error formulation
Observation impact on data assimilation with dynamic background error formulation ALEXANDER BECK alexander.beck@univie.ac.at Department of, Univ. Vienna, Austria Thanks to: Martin Ehrendorfer, Patrick
More informationBackground 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
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 informationObjective localization of ensemble covariances: theory and applications
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
More informationUse of satellite data in ALADIN HARMONIE/Norway
Use of satellite data in ALADI HARMOIE/orway Andrea Storto and Roger Randriamampianina orwegian Meteorological Institute, Oslo orway andrea.storto@met.no, rogerr@met.no Abstract The ALADI/HARMOIE three-dimensional
More informationCurrent Limited Area Applications
Current Limited Area Applications Nils Gustafsson SMHI Norrköping, Sweden nils.gustafsson@smhi.se Outline of talk (contributions from many HIRLAM staff members) Specific problems of Limited Area Model
More informationProgress towards better representation of observation and background errors in 4DVAR
Progress towards better representation of observation and background errors in 4DVAR Niels Bormann 1, Massimo Bonavita 1, Peter Weston 2, Cristina Lupu 1, Carla Cardinali 1, Tony McNally 1, Kirsti Salonen
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 informationAssimilation of the IASI data in the HARMONIE data assimilation system
Assimilation of the IASI data in the HARMONIE data assimilation system Roger Randriamampianina Acknowledgement: Andrea Storto (met.no), Andrew Collard (ECMWF), Fiona Hilton (MetOffice) and Vincent Guidard
More informationEnsemble-variational assimilation with NEMOVAR Part 2: experiments with the ECMWF system
Ensemble-variational assimilation with NEMOVAR Part 2: experiments with the ECMWF system La Spezia, 12/10/2017 Marcin Chrust 1, Anthony Weaver 2 and Hao Zuo 1 1 ECMWF, UK 2 CERFACS, FR Marcin.chrust@ecmwf.int
More informationCurrent Issues and Challenges in Ensemble Forecasting
Current Issues and Challenges in Ensemble Forecasting Junichi Ishida (JMA) and Carolyn Reynolds (NRL) With contributions from WGNE members 31 th WGNE Pretoria, South Africa, 26 29 April 2016 Recent trends
More informationEnsemble aerosol forecasts and assimila1on at ECMWF
Ensemble aerosol forecasts and assimila1on at ECMWF Angela Benede*, Miha Razinger, Luke Jones & Jean- Jacques Morcre
More informationA new operational convection-permitting NWP system for tropical cyclone forecasting in the SW Indian Ocean
A new operational convection-permitting NWP system for tropical cyclone forecasting in the SW Indian Ocean Olivier BOUSQUET, Soline BIELLI, David BARBARY, Marie-Dominique LEROUX, Christelle BARTHE, Dominique
More informationSatellite Observations of Greenhouse Gases
Satellite Observations of Greenhouse Gases Richard Engelen European Centre for Medium-Range Weather Forecasts Outline Introduction Data assimilation vs. retrievals 4D-Var data assimilation Observations
More 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) 1 AIMS 2 OBSERVATION ERRORS
More informationMet Office convective-scale 4DVAR system, tests and improvement
Met Office convective-scale 4DVAR system, tests and improvement Marco Milan*, Marek Wlasak, Stefano Migliorini, Bruce Macpherson Acknowledgment: Inverarity Gordon, Gareth Dow, Mike Thurlow, Mike Cullen
More informationRetrospective and near real-time tests of GSIbased EnKF-Var hybrid data assimilation system for HWRF with airborne radar data
Retrospective and near real-time tests of GSIbased EnKF-Var hybrid data assimilation system for HWRF with airborne radar data Xuguang Wang, Xu Lu, Yongzuo Li University of Oklahoma, Norman, OK In collaboration
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 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 informationRecent 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.
More informationBias Correction of Satellite Data in GRAPES-VAR
Bias Correction of Satellite Data in GRAPES-VAR Wei Han Chinese Academy of Meteorological Sciences, CMA ITSC15, 2006-10, Italy Outline Status of GRAPES-3DVAR Main components of GRAPES Usage of satellite
More informationRecent activities related to EPS (operational aspects)
Recent activities related to EPS (operational aspects) Junichi Ishida and Carolyn Reynolds With contributions from WGE members 31th WGE Pretoria, South Africa, 26 29 April 2016 GLOBAL 2 Operational global
More informationEFSO and DFS diagnostics for JMA s global Data Assimilation System: their caveats and potential pitfalls
EFSO and DFS diagnostics for JMA s global Data Assimilation System: their caveats and potential pitfalls Daisuke Hotta 1,2 and Yoichiro Ota 2 1 Meteorological Research Institute, Japan Meteorological Agency
More informationICON. 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
More informationEnsemble 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
More informationModelling 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
More informationComparing Variational, Ensemble-based and Hybrid Data Assimilations at Regional Scales
Comparing Variational, Ensemble-based and Hybrid Data Assimilations at Regional Scales Meng Zhang and Fuqing Zhang Penn State University Xiang-Yu Huang and Xin Zhang NCAR 4 th EnDA Workshop, Albany, NY
More 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 informationCross-validation methods for quality control, cloud screening, etc.
Cross-validation methods for quality control, cloud screening, etc. Olaf Stiller, Deutscher Wetterdienst Are observations consistent Sensitivity functions with the other observations? given the background
More informationGSI Tutorial Background and Observation Errors: Estimation and Tuning. Daryl Kleist NCEP/EMC June 2011 GSI Tutorial
GSI Tutorial 2011 Background and Observation Errors: Estimation and Tuning Daryl Kleist NCEP/EMC 29-30 June 2011 GSI Tutorial 1 Background Errors 1. Background error covariance 2. Multivariate relationships
More informationSevere weather warnings at the Hungarian Meteorological Service: Developments and progress
Severe weather warnings at the Hungarian Meteorological Service: Developments and progress István Ihász Hungarian Meteorological Service Edit Hágel Hungarian Meteorological Service Balázs Szintai Department
More informationEnsemble-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
More informationQUALITY CONTROL OF WINDS FROM METEOSAT 8 AT METEO FRANCE : SOME RESULTS
QUALITY CONTROL OF WINDS FROM METEOSAT 8 AT METEO FRANCE : SOME RESULTS Christophe Payan Météo France, Centre National de Recherches Météorologiques, Toulouse, France Astract The quality of a 30-days sample
More informationLatest thoughts on stochastic kinetic energy backscatter - good and bad
Latest thoughts on stochastic kinetic energy backscatter - good and bad by Glenn Shutts DARC Reading University May 15 2013 Acknowledgments ECMWF for supporting this work Martin Leutbecher Martin Steinheimer
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 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 informationAssimilating cloud and precipitation: benefits and uncertainties
Assimilating cloud and precipitation: benefits and uncertainties Alan Geer Thanks to: Katrin Lonitz, Peter Lean, Richard Forbes, Cristina Lupu, Massimo Bonavita, Mats Hamrud, Philippe Chambon, Fabrizio
More informationThe assimilation of AMSU-A radiances in the NWP model ALADIN. The Czech Hydrometeorological Institute Patrik Benáček 2011
The assimilation of AMSU-A radiances in the NWP model ALADIN The Czech Hydrometeorological Institute Patrik Benáček 2011 Outline Introduction Sensor AMSU-A Set-up of model ALADIN Set-up of experiments
More information