Training course: Data assimilation
from Monday, 10 May 2021 (08:45) to Friday, 14 May 2021 (17:00)
-
Monday, 10 May 202108:45 Online course opensOnline course opens08:45 - 09:0509:05 Introduction - Sarah Keeley (ECMWF) Andy Brown (ECMWF)Introduction
- Sarah Keeley (ECMWF)
- Andy Brown (ECMWF)
09:05 - 10:05In this session we will sort out: general house keeping for the course, such as computing accounts; have an overview of ECMWF research and introduce ourselves to one another.10:05 BreakBreak10:05 - 10:4010:40 Overview of assimilation methods - Massimo Bonavita (ECMWF)Overview of assimilation methods- Massimo Bonavita (ECMWF)
10:40 - 11:40The goal of the ECMWF Earth System data assimilation is to provide an accurate and physically coherent description of the state of the atmosphere, ocean, sea ice and land surface as an initial point for our forecasts. This requires blending in a statistically optimal way information from a huge variety of observations and our prior knowledge about the physical laws of the Earth system, which is encapsulated in our models. In this lecture we will lay the general conceptual framework on how to achieve this from a Bayesian perspective. We will then highlight the approximations and hypotheses which are required to make the assimilation problem computationally tractable and which underlie the practical data assimilation algorithms which will be described in detail in this training course. By the end of lecture you should be able to: • understand the basics of how a geophysical data assimilation system works; • understand the main approximations and hypotheses which are required to build practical data assimilation algorithms for large geophysical systems11:40 Lunch breakLunch break11:40 - 12:5012:50 Assimilation Algorithms: (1) Basic concepts - Sebastien Massart (ECMWF)Assimilation Algorithms: (1) Basic concepts- Sebastien Massart (ECMWF)
12:50 - 13:50This lecture will explain the basic concepts of the assimilation algorithms. The terminology used in the next lectures will be introduced. Simple examples will conduce towards the formulation of the optimal minimum-variance analysis. The optimal interpolation method will finally be presented. By the end of the lecture the participants should be able to: • Recognize the notations used for the rest of the week • Solve the optimal minimum-variance analysis problem • Apply the optimal interpolation method13:50 BreakBreak13:50 - 14:1514:15 Assimilation Algorithms: (2) 3D-Var - Sebastien Massart (ECMWF)Assimilation Algorithms: (2) 3D-Var- Sebastien Massart (ECMWF)
14:15 - 15:15This lecture will present the 3D-Var assimilation algorithm. This algorithm is based in the formulation of a cost function to minimize. Minimization methods will be presented together with some information on how to improve their efficiency. By the end of the lecture the participants should be able to: • Recognize the 3D-Var cost function • Explain the various terms of the cost function • Question the efficiency of methods designed to find the minimum of the cost function15:15 Comfort breakComfort break15:15 - 15:3015:30 Conventional and actively sensed observations - Sean Healy (ECMWF)Conventional and actively sensed observations- Sean Healy (ECMWF)
15:30 - 16:30This lecture will introduce how observations are an essential part of the data assimilation system. It will focus on in situ (also called conventional) observations, from surface stations, drifters, aircraft and radiosondes. They are important both for direct use in the data assimilation system and for diagnostics. Radiosonde and surface observations also help to control the biases in the assimilation system. However, they are diverse and hey can be complex, so close attention to quality control, observation uncertainty and (in some cases) bias correction is needed to optimise their use. The use of new BUFR format high resolution radiosonde data will also be presented. The lecture will also introduce the actively sensed satellite observations used for data assimilation at ECMWF: radio occultation data, scatterometer winds, and altimeter wind/significant wave height. By the end of the lecture the student should be able to: • understand how in situ and actively sensed observations are used in data assimilation, including bias aspects and observation uncertainty aspects. • appreciate the diverse and complex range of in situ observations used in modern NWP. • understand how radio occultation data, scatterometer winds and altimeter data are used in data assimilation. -
Tuesday, 11 May 202108:45 Online course opensOnline course opens08:45 - 09:0509:05 Analysis of radiance observations - Tony McNally (ECMWF)Analysis of radiance observations
- Tony McNally (ECMWF)
09:05 - 10:05The primary purpose of this lecture is to explore the implications of the fact that satellites can only measure radiation at the top of the atmosphere and do not measure the geophysical variables we require for NWP (e.g. temperature, humidity and wind). The link between the atmospheric variables and the measured radiances is the radiative transfer equation - the key elements of which are discussed. It is shown how - with careful frequency selection - satellite measurements can be made for which the relationship to geophysical variables is greatly simplified. Despite these simplifications, it is shown that the extraction of detailed profile information from downward looking radiance measurements is a formally ill posed inverse problem. Data assimilation is introduced as the solution to this inverse problem, where background information and satellite observations are combined to produce a best or optimal estimate of the atmospheric state. The main elements of the assimilation scheme (such as the chain of observation operators for radiances) and its key statistical inputs are examined. In particular, it is shown that incorrect specification of observation errors (R) and background errors (B) can severely limit the successful exploitation of satellite data. By the end of this lecture you will: • understand exactly what a satellite actually measures (radiance) • appreciate the complex relationship between what is measured and what we wish to know for NWP • how information is extracted from satellite measurements in data assimilation10:05 BreakBreak10:05 - 10:4010:40 Assimilation Algorithms: (3) 4D-Var - Sebastien Massart (ECMWF)Assimilation Algorithms: (3) 4D-Var- Sebastien Massart (ECMWF)
10:40 - 11:4011:40 Lunch breakLunch break11:40 - 12:5012:50 Assimilation Algorithms: (4) Ensemble Kalman filters - Massimo Bonavita (ECMWF)Assimilation Algorithms: (4) Ensemble Kalman filters- Massimo Bonavita (ECMWF)
12:50 - 13:50The aim of this lecture is to introduce the concept of the EnKF in the context of atmospheric data assimilation. Strengths and weaknesses of the algorithm will be discussed and results of the ECMWF implementation will be presented. By the end of the lecture the participants should be able to: • Describe the basic EnKF algorithm and its connections with the Kalman Filter; • Discuss some of the advantages and the limitations of EnKF algorithms with respect to more established variational algorithms; • Be aware of recent developments in hybrid variational-EnKF data assimilation13:50 BreakBreak13:50 - 14:1514:15 Bias correction methods - Niels Bormann (ECMWF)Bias correction methods- Niels Bormann (ECMWF)
14:15 - 15:15In this lecture, the variational bias correction scheme (VarBC) as used at ECMWF is explained. VarBC replaced the tedious job of estimating observation bias off-line for each satellite instrument or in-situ network by an automatic self-adaptive system. This is achieved by making the bias estimation an integral part of the ECMWF variational data assimilation system, where now both the initial model state and observation bias estimates are updated simultaneously. By the end of the session you should be able to realize that: • many observations are biased, and that the characteristics of bias varies widely between types of instruments • separation between model bias and observation bias is often difficult • the success of an adaptive system implicitly relies on a redundancy in the underlying observing system.15:15 Comfort breakComfort break15:15 - 15:3015:30 Tangent linear and adjoints - Marcin Chrust (ECMWF)Tangent linear and adjoints- Marcin Chrust (ECMWF)
15:30 - 16:30The goal of this lecture is to familiarise the student with the notion of tangent linear and adjoint models, and their use in variational data assimilation. A general overview of the current use of tangent linear and adjoint models in the ECMWF system will also be provided. Theoretical definitions and practical examples of tangent liner and adjoint models will be given. The student will be invited to work some simple tangent linear and adjoint derivations together with the instructor. A brief introduction to automatic differentiation software will also be given. By the end of the session you should be able to: • define what tangent linear and adjoint models are • derive tangent linear and adjoint equations for a simple nonlinear equation • describe the use of tangent linear and adjoint codes within the ECMWF's 4D-VAR system.16:30 BreakBreak16:30 - 16:4516:45 Practical session: Tangent linear and adjoints - Sebastien Massart (ECMWF) Marcin Chrust (ECMWF)Practical session: Tangent linear and adjoints- Sebastien Massart (ECMWF)
- Marcin Chrust (ECMWF)
16:45 - 17:45 -
Wednesday, 12 May 202108:45 Online course opensOnline course opens08:45 - 09:0509:05 Background error modelling in data assimilation - Elias Holm (ECMWF)Background error modelling in data assimilation
- Elias Holm (ECMWF)
09:05 - 10:05The background error is central to the performance of the analysis system and tells how much confidence to put in the best available forecast which is to be updated with new observations. The lecture will review how background errors are estimated and represented for current variational algorithms.10:05 BreakBreak10:05 - 10:4010:40 Data assimilation diagnostics - forecast sensitivity - Bruce Ingleby (ECMWF)Data assimilation diagnostics - forecast sensitivity- Bruce Ingleby (ECMWF)
10:40 - 11:4011:40 Lunch breakLunch break11:40 - 12:5012:50 Model error in data assimilation - Patrick Laloyaux (ECMWF)Model error in data assimilation- Patrick Laloyaux (ECMWF)
12:50 - 13:50In this lecture, the impact of model error on variational data assimilation will be presented. This lecture will introduce weak-constraint 4D-Var as a way to account for model error in the data assimilation process. Several examples of results from simplified implementations in the IFS will be shown. By the end of the lecture the participants should be able to: • describe the impact of model error on the data assimilation process, • explain the difficulties in properly accounting for model error in data assimilation.13:50 BreakBreak13:50 - 14:1514:15 Assimilation Algorithms: (5) Hybrid data assimilation methds - Massimo Bonavita (ECMWF)Assimilation Algorithms: (5) Hybrid data assimilation methds- Massimo Bonavita (ECMWF)
14:15 - 15:1515:15 Comfort breakComfort break15:15 - 15:3015:30 Practical session: DA experiments with OOPS - Massimo Bonavita (ECMWF) Marcin Chrust (ECMWF) Sebastien Massart (ECMWF) Patrick Laloyaux (ECMWF)Practical session: DA experiments with OOPS- Massimo Bonavita (ECMWF)
- Marcin Chrust (ECMWF)
- Sebastien Massart (ECMWF)
- Patrick Laloyaux (ECMWF)
15:30 - 17:15 -
Thursday, 13 May 202108:45 Online course opensOnline course opens08:45 - 09:0509:05 Parametrization and data assimilation - Philippe Lopez (ECMWF)Parametrization and data assimilation
- Philippe Lopez (ECMWF)
09:05 - 10:05This one-hour lecture will identify the challenges associated with the use of physical parametrizations in the context of four-dimensional variational data assimilation (4D-Var). The importance of the linearity constraint in 4D-Var and the methods to address it will be detailed. The set of linearized physical parametrizations used at ECMWF will be briefly presented. Examples of the use of physical parametrizations in variational data assimilation and its impact on forecast quality will be given. By the end of the lecture, the students should be able: • to tell why physical parametrizations are needed in data assimilation. • to recognize the importance of the regularization of the linearized code10:05 BreakBreak10:05 - 10:4010:40 Data Assimilation and Machine Learning - Alan Geer (ECMWF)Data Assimilation and Machine Learning- Alan Geer (ECMWF)
10:40 - 11:4011:40 Lunch breakLunch break11:40 - 12:5012:50 Land data assimilation - Patricia de Rosnay (ECMWF)Land data assimilation- Patricia de Rosnay (ECMWF)
12:50 - 13:50The aim of these sessions is to understand the role of land surface data assimilation on medium range weather forecasts. We will give an overview of the different approaches used to assimilate land surface data and to initialise model variables in NWP. We will present the current observing systems and describe the land data assimilation structure within ECMWF system. By the end of the session you should be able to: • identify the different observations used for snow and soil moisture data assimilation • define land surface data assimilation approaches used for NWP • describe the role of land surface data assimilation on medium-range weather forecasts13:50 BreakBreak13:50 - 14:1514:15 Gather Town Practical Session: DA experiments with OOPS - Marcin Chrust (ECMWF) Patrick Laloyaux (ECMWF) Massimo Bonavita (ECMWF)Gather Town Practical Session: DA experiments with OOPS- Marcin Chrust (ECMWF)
- Patrick Laloyaux (ECMWF)
- Massimo Bonavita (ECMWF)
14:15 - 15:1515:15 Comfort breakComfort break15:15 - 15:3015:30 Practical Session: DA experiments with OOPS - Massimo Bonavita (ECMWF) Marcin Chrust (ECMWF) Patrick Laloyaux (ECMWF) Sebastien Massart (ECMWF)Practical Session: DA experiments with OOPS- Massimo Bonavita (ECMWF)
- Marcin Chrust (ECMWF)
- Patrick Laloyaux (ECMWF)
- Sebastien Massart (ECMWF)
15:30 - 16:30 -
Friday, 14 May 202108:45 Online course opensOnline course opens08:45 - 09:0509:05 Data assimilation of atmospheric composition - Nicolas Bousserez (ECMWF)Data assimilation of atmospheric composition
- Nicolas Bousserez (ECMWF)
09:05 - 10:05At ECMWF, atmospheric composition data are assimilated into the IFS as part of the Copernicus Atmospheric Monitoring Service. On a global scale, atmospheric composition represents the full state of the global atmosphere covering phenomena such as desert dust plumes, long-range transport of atmospheric pollutants or ash plumes from volcanic eruptions, but also variations and long-term changes in the background concentrations of greenhouse gases. The aim of this lecture is to give an overview of the work that is carried out at ECMWF regarding the assimilation of atmospheric composition data, and to address why this is of interest and which special challenges are faced when assimilating atmospheric composition data. By the end of the session you should: • have some understanding of the work carried out at ECMWF to assimilate data of atmospheric composition10:05 BreakBreak10:05 - 10:4010:40 Ocean data assimilation - Hao Zuo (ECMWF)Ocean data assimilation- Hao Zuo (ECMWF)
10:40 - 11:40This lecture provides an overview of a typical ocean data assimilation system for initialization and re-analyses application. The lecture uses as an example the ECMWF ocean data assimilation system, which is based the NEMOVAR (3Dvar FGAT). This will be used to discuss design of the assimilation cycle, formulation of error covariances, observations assimilated and evaluation procedure, among others. By the end of the lecture students should be able to: • describe the different components involved in an ocean data assimilation system • list the commonalities and differences between ocean and atmosphere data assimilation • describe the basics of the physical ocean observing system • explain the essential multivariate relationships between ocean variables • identify the limitations of the existing systems.11:40 Lunch breakLunch break11:40 - 12:5012:50 Reanalysis methods - Hans Hersbach (ECMWF)Reanalysis methods- Hans Hersbach (ECMWF)
12:50 - 13:50The aim of this session is to understand how data assimilation can improve our knowledge of past weather over long time-scales. We will present recent advances that help capture changes over time in observing system networks, and project this variation in information content into uncertainty estimates of the reanalysis products. We will also discuss the applications of reanalysis, which generally put weather events into the climate context. By the end of the session you should be able to: • explain what are the goals of data assimilation in a reanalysis data assimilation system • list the key aspects that require particular attention in reanalysis, as compared to numerical weather prediction • describe the most common problems in reanalysis products13:50 BreakBreak13:50 - 14:1514:15 Coupled data assimilation opportunities and challenges - Phil Browne (ECMWF)Coupled data assimilation opportunities and challenges- Phil Browne (ECMWF)
14:15 - 15:15At ECMWF, we are striving to move towards an Earth System approach to our data assimilation techniques. We currently have models not only of the atmosphere, but of the ocean, the land surface, sea ice, waves, and atmospheric composition. These systems interact with each other in different ways and all need to be initialised through the incorporation of observational data. The aim of this lecture is to recognise the benefits and challenges associated with data assimilation in coupled models. By the end of the lecture the participants should be able to: • Recall the challenges associated with variational data assimilation in systems with different timescales and computer codes. • Describe the benefits of having more consistently balanced coupled systems from coupled data assimilation. • Explain the differences between weakly and strongly coupled data assimilation approaches. • Discuss the various methods that are in use at ECWMF and explain the planned developments of the systems.15:15 Comfort breakComfort break15:15 - 15:3015:30 Final discussion and Q&A - Elias Holm Patrick Laloyaux (ECMWF) Sebastien Massart (ECMWF) Massimo Bonavita (ECMWF) Marcin Chrust (ECMWF)Final discussion and Q&A- Elias Holm
- Patrick Laloyaux (ECMWF)
- Sebastien Massart (ECMWF)
- Massimo Bonavita (ECMWF)
- Marcin Chrust (ECMWF)
15:30 - 16:30