On Some Computational Aspects of the Variational Data Assimilation Techniques

It is important to incorporate all available observations when large-scale mathematical models arising in different fields of science and engineering are used to study various physical and chemical processes. Variational data assimilation techniques can be used in the attempts to utilize efficiently observations in a large-scale model. Variational data assimilation techniques are based on a combination of three very important components - numerical methods for solving differential equations, - splitting procedures and - optimization algorithms. It is crucial to select an optimal (or, at least, a good) combination of these three components, because models which are very expensive computationally become much more expensive (the computing time being often increased by a factor greater than 100) when a variational data assimilation technique is applied. Therefore, it is important to study the interplay between the three components of the variational data assimilation techniques as well as to apply powerful parallel computers in the computations. Some results obtained in the search for a good combination will be reported. Parallel techniques described in [1] are used in the runs related to this paper. Modules from a particular large-scale mathematical model, the Unified Danish Eulerian Model (UNI-DEM), are used in the experiments. The mathematical background of UNI-DEM is discussed in [1], [24] The ideas are rather general and can easily be applied in connection with other mathematical models.

[1]  T. E. Hull,et al.  Comparing Numerical Methods for Ordinary Differential Equations , 1972 .

[2]  Bruno Sportisse,et al.  Solving reduced chemical models in air pollution modelling , 2003 .

[3]  Peter Borrell,et al.  Transport and chemical transformation of pollutants in the troposphere : an overview of the work of EUROTRAC , 2000 .

[4]  Zahari Zlatev,et al.  Computer Treatment of Large Air Pollution Models , 1995 .

[5]  Ionel M. Navon,et al.  Second-Order Information in Data Assimilation* , 2002 .

[6]  Philip E. Gill,et al.  Practical optimization , 1981 .

[7]  E. Hairer,et al.  Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems , 2010 .

[8]  Gurij Ivanovich Marchuk,et al.  Some application of splitting-up methods to the solution of mathematical physics problems , 1968 .

[9]  J. Brandts [Review of: W. Hundsdorfer, J.G. Verwer (2003) Numerical Solution of Time-Dependent Advection-Diffusion-Reaction Equations] , 2006 .

[10]  J. M. Lewis,et al.  The use of adjoint equations to solve a variational adjustment problem with advective constraints , 1985 .

[11]  Adrian Sandu,et al.  Adjoint sensitivity analysis of regional air quality models , 2005 .

[12]  Hauke Schmidt,et al.  A four-dimensional variational chemistry data assimilation scheme for Eulerian chemistry transport modeling , 1999 .

[13]  M. C. Dodge,et al.  A photochemical kinetics mechanism for urban and regional scale computer modeling , 1989 .

[14]  Zahari Zlatev,et al.  Parallel runs of a large air pollution model on a grid of Sun computers , 2004, Math. Comput. Simul..

[15]  Adolf Ebel,et al.  Chemical Transfer and Transport Modelling , 2000 .

[16]  Zahari Zlatev,et al.  Chapter 1 - PDE systems arising in air pollution modelling and justification of the need for high speed computers , 2006 .

[17]  P. Gill,et al.  Conjugate-Gradient Methods for Large-Scale Nonlinear Optimization. , 1979 .

[18]  Zahari Zlatev,et al.  Modified Diagonally Implicit Runge–Kutta Methods , 1981 .

[19]  J. Verwer,et al.  Numerical solution of time-dependent advection-diffusion-reaction equations , 2003 .

[20]  Shengtai Li,et al.  Adjoint Sensitivity Analysis for Differential-Algebraic Equations: The Adjoint DAE System and Its Numerical Solution , 2002, SIAM J. Sci. Comput..

[21]  O. Talagrand,et al.  4D-variational data assimilation with an adjoint air quality model for emission analysis , 2000, Environ. Model. Softw..

[22]  G. Strang On the Construction and Comparison of Difference Schemes , 1968 .

[23]  J. Butcher Numerical methods for ordinary differential equations , 2003 .

[24]  Jostein K. Sundet,et al.  Model calculations of present and future levels of ozone and ozone precursors with a global and a regional model , 2001 .

[25]  Graeme Fairweather,et al.  The current state and future direction of Eulerian models in simulating the tropospheric chemistry and transport of trace species: a review , 1995 .

[26]  Adolf Ebel,et al.  Regional Modelling of Tropospheric Ozone Distribution and Budgets , 1997 .