1 / 58

Errors, Uncertainties in Data Assimilation

Errors, Uncertainties in Data Assimilation. François-Xavier LE DIMET Université Joseph Fourier+INRIA Projet IDOPT, Grenoble, France. Acknowlegment. Pierre Ngnepieba ( FSU) Youssuf Hussaini ( FSU) Arthur Vidard ( ECMWF) Victor Shutyaev ( Russ. Acad. Sci.) Junqing Yang ( LMC , IDOPT).

mave
Télécharger la présentation

Errors, Uncertainties in Data Assimilation

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Errors, Uncertainties in Data Assimilation François-Xavier LE DIMET Université Joseph Fourier+INRIA Projet IDOPT, Grenoble, France

  2. Acknowlegment • Pierre Ngnepieba ( FSU) • Youssuf Hussaini ( FSU) • Arthur Vidard ( ECMWF) • Victor Shutyaev ( Russ. Acad. Sci.) • Junqing Yang ( LMC , IDOPT)

  3. Prediction: What information is necessary ? • Model • law of conservation mass, energy • Laws of behaviour • Parametrization of physical processes • Observations in situ and/or remote • Statistics • Images

  4. Forecast.. • Produced by the integration of the model from an initial condition • Problem : how to link together heterogeneous sources of information • Heterogeneity in : • Nature • Quality • Density

  5. Basic Problem • U and V control variables, V being and error on the model • J cost function • U* and V* minimizes J

  6. Optimality System • P is the adjoint variable. • Gradients are couputed by solving the adjoint model then an optimization method is performed.

  7. Remark on statistical information • Statistical information is included in the assimilation • In the norm of the discrepancy between the solution of the model ( approximation of the inverse of the covariance matrix) • In the background term ( error covariance matrix)

  8. Remarks: • This method is used since May 2000 for operational prediction at ECMWF and MétéoFrance, Japanese Meteorological Agency ( 2005) with huge models ( 10 millions of variable. • The Optimality System countains all the available information • The O.S. should be considered as a « Generalized Model » • Only the O.S. makes sense.

  9. Errors • On the model • Physical approximation (e.g. parametrization of subgrid processes) • Numerical discretization • Numerical algorithms ( stopping criterions for iterative methods • On the observations • Physical measurement • Sampling • Some « pseudo-observations », from remote sensing, are obtained by solving an inverse problem.

  10. Sensitivity of the initial condition with respect to errors on the models and on the observations. • The prediction is highly dependant on the initial condition. • Models have errors • Observations have errors. • What is the sensitivity of the initial condition to these errors ?

  11. Optimality System : including errors on the model and on the observation

  12. Second order adjoint

  13. Models and Data • Is it necessary to improve a model if data are not changed ? • For a given model what is the « best » set of data? • What is the adequation between models and data?

  14. A simple numerical experiment. • Burger’s equation with homegeneous B.C.’s • Exact solution is known • Observations are without error • Numerical solution with different discretization • The assimilation is performed between T=0 and T=1 • Then the flow is predicted at t=2.

  15. Partial Conclusion • The error in the model is introduced through the discretization • The observations remain the same whatever be the discretization • It shows that the forecast can be downgraded if the model is upgraded. • Only the quality of the O.S. makes sense.

  16. Remark 1 • How to improve the link between data and models? • C is the operator mapping the space of the state variable into the space of observations • We considered the liear case.

  17. Remark 2 : ensemble prediction • To estimate the impact of uncertainies on the prediction several prediction are performed with perturbed initial conditions • But the initial condition is an artefact : there is no natural error on it . The error comes from the data throughthe data assimilation process • If the error on the data are gaussian : what about the initial condition?

  18. Because D.A. is a non linear process then the initial condition is no longer gaussian

  19. Control of the error

  20. Choice of the base

  21. Remark . • The model has several sources of errors • Discretization errors may depends on the second derivative : we can identify this error in a base of the first eigenvalues of the Laplacian • The systematic error may depends be estimated using the eigenvalues of the correlation matrix

  22. Numerical experiment • With Burger’s equation • Laplacian and covariance matrix have considered separately then jointly • The number of vectors considered in the correctin term varies

  23. An application in oceanographyin A. Vidard’s Ph.D. • Shallow water on a square domain with a flat bottom. • An bias term is atted into the equation and controlled

  24. RMS ot the sea surface height with or without control of the bias

  25. An application in hydrology(Yang Junqing ) • Retrieve the evolution of a river • With transport+sedimentation

  26. Physical phenomena • fluid and solid transport • different time scales

  27. N 2D sedimentation modeling 1. Shallow-water equations 2. Equation of constituent concentration 3. Equation of the riverbed evolution

  28. Semi-empirical formulas • Bed load function : • Suspended sediment transport rate : are empirical constants

  29. An example of simulation Initial river bed • Domain : • Space step : 2 km in two directions • Time step : 120 seconds Simulated evolution of river bed (50 years)

  30. Model error estimation controlled system • model • cost function • optimality conditions • adjoint system(to calculate the gradient)

  31. Reduction of the size of the controlled problem • Change the space bases Suppose is a base of the phase space and is time-dependent base function on [0, T], so that then the controlled variables are changed to with controlled space size

More Related