Development of a data assimilation algorithm. (English) Zbl 1142.65394
Summary: 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
\(\bullet \) numerical methods for solving differential equations,
\(\bullet \) splitting procedures and
\(\bullet \) 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 by V.N. Alexandrov, W. Owczarz, P.G. Thomsen and Z. Zlatev [Parallel runs of a large air pollution model on a grid of Sun computers, Mathematics and Computers in Simulation 65, 557–577 (2004)] 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 [loc. cit. ; 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.
\(\bullet \) numerical methods for solving differential equations,
\(\bullet \) splitting procedures and
\(\bullet \) 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 by V.N. Alexandrov, W. Owczarz, P.G. Thomsen and Z. Zlatev [Parallel runs of a large air pollution model on a grid of Sun computers, Mathematics and Computers in Simulation 65, 557–577 (2004)] 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 [loc. cit. ; 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.
MSC:
| 65L99 | Numerical methods for ordinary differential equations |
Keywords:
scientific models; splitting techniques; systems of PDEs; systems of ODEs; stiffness; accuracy; data assimilationReferences:
| [1] | Alexandrov, V. N.; Owczarz, W.; Thomsen, P. G.; Zlatev, Z., Parallel runs of a large air pollution model on a grid of Sun computers, Mathematics and Computers in Simulation, 65, 557-577 (2004) |
| [2] | Zlatev, Z., Computer Treatment of Large Air Pollution Models (1995), Kluwer Academic Publishers: Kluwer Academic Publishers Dordrecht, Boston, London · Zbl 0852.65058 |
| [3] | Elbern, H.; Schmidt, H., A four-dimensional variational chemistry data assimilation scheme for Eulerian chemistry transport modelling, Journal of Geophysical Research, 104, 18583-18598 (1999) |
| [4] | Elbern, H.; Schmidt, H.; Talagrand, O.; Ebel, A., 4D-variational data assimilation with an adjoint air quality model for emission analysis, Environmental Modelling & Software, 15, 539-548 (2000) |
| [5] | Sandu, A.; Daescu, D. N.; Carmichael, G. R.; Chai, T., Adjoint sensitivity analysis of regional air quality models, Journal of Computational Physics, 204, 222-252 (2005) · Zbl 1061.92061 |
| [6] | Le Dimet, F.-X.; Navon, I. M.; Daescu, D. N., Second order information in data assimilation, Monthly Weather Review, 130, 629-648 (2002) |
| [7] | Lewis, J. M.; Derber, J. C., The use of adjoint equations to solve a variational adjustment problem with advective constraints, Tellus, 37A, 309-322 (1985) |
| [8] | Vreugdenhil, C. B., Chapter 2: Linear central difference methods, (Vreugdenhil, C. B.; Koren, B., Numerical Methods for Advection-Diffusion Problems. Numerical Methods for Advection-Diffusion Problems, Notes in Numerical Fluid Mechanics, vol. 45 (1993), Vieweg: Vieweg Braunschwieg, Wiesbaden), 27-54 · Zbl 0805.76042 |
| [9] | NAG Library Fortran Manual: E04 - minimizing and maximizing a function, Numerical Algorithms Group (NAG), Banbury Road 7, Oxford, England. http://www.nag.co.uk, 2004; NAG Library Fortran Manual: E04 - minimizing and maximizing a function, Numerical Algorithms Group (NAG), Banbury Road 7, Oxford, England. http://www.nag.co.uk, 2004 |
| [10] | P.E. Gill, W. Murray, Safeguarded steplength algorithms for optimization using descent methods, NPL Report NAC 37, National Physical Laboratory, Teddington, Middlesex, UK, 1973; P.E. Gill, W. Murray, Safeguarded steplength algorithms for optimization using descent methods, NPL Report NAC 37, National Physical Laboratory, Teddington, Middlesex, UK, 1973 |
| [11] | Gill, P. E.; Murray, W., Newton-type methods for unconstrained and linearly constrained optimization, Mathematical Programming, 7, 311-350 (1974), Teddington, Middlesex, UK, 1973 · Zbl 0297.90082 |
| [12] | P.E. Gill, W. Murray, Minimization subject to bounds on the variables, NPL Report NAC 72, National Physical Laboratory, Teddington, Middlesex, UK, 1976; P.E. Gill, W. Murray, Minimization subject to bounds on the variables, NPL Report NAC 72, National Physical Laboratory, Teddington, Middlesex, UK, 1976 |
| [13] | P.E. Gill, W. Murray, Conjugate-gradient methods for large-scale nonlinear optimization, Technical Report SOL 79-15, Department of Operations Research, Stanford University, Stanford, CA, USA, 1979; P.E. Gill, W. Murray, Conjugate-gradient methods for large-scale nonlinear optimization, Technical Report SOL 79-15, Department of Operations Research, Stanford University, Stanford, CA, USA, 1979 |
| [14] | Gill, P. E.; Murray, W.; Wright, M. H., Practical Optimization (1981), Academic Press: Academic Press New York · Zbl 0503.90062 |
| [15] | Zlatev, Z.; Dimov, I., (Computational and Environmental Challenges in Environmental Modelling. Computational and Environmental Challenges in Environmental Modelling, Studies in Computational Mathematics, vol. 13 (2006), Elsevier Science: Elsevier Science Amsterdam) · Zbl 1120.65103 |
| [16] | Djouad, R.; Sportisse, B., Solving reduced chemical models in air pollution modelling, Applied Numerical Mathematics, 44, 49-61 (2003) · Zbl 1038.92037 |
| [17] | Ebel, A., Chemical transfer and transport modelling, (Borrell, P.; Borrell, P. M., Transport and Chemical Transformation of Pollutants in the Troposphere (2000), Springer: Springer Berlin), 85-128 |
| [18] | Ebel, A.; Memmesheimer, M.; Jacobs, H., Regional modelling of tropospheric ozone distributions and budgets, (Varotsos, C., Atmospheric Ozone Dynamics: I. Global Environmental Change. Atmospheric Ozone Dynamics: I. Global Environmental Change, NATO ASI Series, vol. 53 (1997), Springer: Springer Berlin), 37-57 |
| [19] | Gery, M. W.; Whitten, G. Z.; Killus, J. P.; Dodge, M. C., A photochemical kinetics mechanism for urban and regional modeling, Journal of Geophysical Research, 94, 12925-12956 (1989) |
| [20] | Jonson, J. E.; Sundet, J.; Tarrason, L., Model calculations of present and future levels of ozone and ozone precursors with a global and a regional model, Atmospheric Environment, 35, 525-537 (2001) |
| [21] | Peters, L. K.; Berkowitz, C. M.; Carmichael, G. R.; Easter, R. C.; Fairweather, G.; Ghan, S. J.; Hales, J. M.; Leung, L. R.; Pennell, W. R.; Potra, F. A.; Saylor, R. D.; Tsang, T. T., The current state and future direction of Eulerian models in simulating tropospheric chemistry and transport of trace species: A review, Atmospheric Environment, 29, 189-222 (1995) |
| [22] | D. Simpson, H. Fagerli, J.E. Jonson, S.G. Tsyro, P. Wind, J.-P. Tuovinen, Transboundary acidification, eutrophication and ground level ozone in Europe, Part I. Unified EMEP model description, EMEP/MSC-W Status Report 1/2003, Norwegian Meteorological Institute, Oslo, Norway, 2003; D. Simpson, H. Fagerli, J.E. Jonson, S.G. Tsyro, P. Wind, J.-P. Tuovinen, Transboundary acidification, eutrophication and ground level ozone in Europe, Part I. Unified EMEP model description, EMEP/MSC-W Status Report 1/2003, Norwegian Meteorological Institute, Oslo, Norway, 2003 |
| [23] | Hairer, E.; Wanner, G., Solving Ordinary Differential Equations, II: Stiff and Differential-Algebraic Problems (1991), Springer: Springer Berlin, Heidelberg, New York, London · Zbl 0729.65051 |
| [24] | Hundsdorfer, W.; Verwer, J. G., Numerical Solution of Time-Dependent Advection-Diffusion-Reaction Equations (2003), Springer: Springer Berlin · Zbl 1030.65100 |
| [25] | Lambert, J. D., Numerical Methods for Ordinary Differential Systems (1991), Wiley: Wiley Chichester, New York, Brisbane, Toronto, Singapore · Zbl 0745.65049 |
| [26] | Zlatev, Z., Modified diagonally implicit Runge-Kutta methods, SIAM Journal on Scientific and Statistical Computing, 2, 321-334 (1981) · Zbl 0475.65040 |
| [27] | Marchuk, G. I., Some application of splitting-up methods to the solution of mathematical physics problems, Applik. Mat., 13, 2, 103-132 (1968) · Zbl 0159.44702 |
| [28] | Strang, G., On the construction and comparison of difference schemes, SIAM Journal on Numerical Analysis, 5, 505-517 (1968) · Zbl 0184.38503 |
| [29] | Cao, Y.; Shengtai, L.; Petzold, L.; Serban, R., Adjoint sensitivity analysis for differential-algebraic equations: The adjoint DAE system and its numerical solution, SIAM Journal on Scientific Computing, 24, 1076-1089 (2003) · Zbl 1034.65066 |
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.
