Development of a data assimilation algorithm

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 (for example, in order to obtain more reliable initial values). 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 will 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 of numerical methods, splitting techniques and optimization algorithms will be reported. Parallel techniques described in [V.N. Alexandrov, W. Owczarz, P.G. Thomsen, Z. Zlatev, Parallel runs of a large air pollution model on a grid of Sun computers, Mathematics and Computers in Simulation, 65 (2004) 557-577] are used in the runs. 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 [V.N. Alexandrov, W. Owczarz, P.G. Thomsen, Z. Zlatev, Parallel runs of a large air pollution model on a grid of Sun computers, Mathematics and Computers in Simulation, 65 (2004) 557-577, Z. Zlatev, Computer Treatment of Large Air Pollution Models, Kluwer Academic Publishers, Dordrecht, Boston, London, 1995]. The ideas are rather general and can easily be applied in connection with other mathematical models.

[1]  I. H. Öğüş,et al.  NATO ASI Series , 1997 .

[2]  J. Lambert Numerical Methods for Ordinary Differential Systems: The Initial Value Problem , 1991 .

[3]  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 .

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

[19]  Philip E. Gill,et al.  Newton-type methods for unconstrained and linearly constrained optimization , 1974, Math. Program..

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

[21]  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 .

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

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

[24]  Z. Zlatev Modified Diagonally Implicit Runge-Kutta Methods , 1980 .

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

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

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