This work is concerned with the reconstruction of unknown atmospheric distributions of species' concentrations and the simultaneous identification of unknown physical parameters. The physical and chemical processes are represented by a simplified chemical transport model. The second focus of the contribution lies on the development and the application of the adaptive techniques to the resulting inverse problem. The underlying theoretical framework is the Dual Weighted Residual (DWR) method for goal-oriented mesh optimization.Index Terms-Atmospheric modeling, inverse problems, goal-oriented adaptivity, mesh optimization.
IntroductionIn this contribution, we showcase the methods stemming from optimal control theory for modeling the transport of chemical species in the Earth's atmosphere by remote sensing techniques such as satellite imaging. The goal of the presented work is the reconstruction of three-dimensional time-dependent concentration fields of chemical constituents in the Earth's atmosphere using indirect measurements such as satellite-based observations. For the latter case, the observations are generally given in form of vertical column densities (VCD).In the last years, significant effort has been spent on refining models to appropriately describe processes affecting the atmospheric constituent transport. The "general" chemical transport model is given by a system of weakly coupled parabolic equations. The coupling between the equations usually occur due to the chemical transformations, affecting the species' concentrations. The use of such models implies that the chemical mechanisms are precisely known. However, there are many situations in which these transformations cannot be modeled properly since the chemical mechanism is only poorly understood. Furthermore, depending on the form of the given measurements and the unknown parameters, the use of such a general model could be not appropriate. To address this issue, we propose a simplified model which allows to reconstruct the constituent concentrations together with unknown sources and sinks, including chemical transformations. However, other parameters can also be the goal of the computations. E.g. in the event that there are many candidates for modeling the chemical processes of a species, the reconstructed source/sink terms could also serve as a guideline for choosing the "right" mechanism.The resulting problems are solved using the Galerkin finite element method. On uniformly refined meshes this method leads to the same algebraic systems as the discretization using finite differences. However, the Galerkin discretization offers the possibility for the systematic a priori and a posteriori error analysis. Moreover, the a posteriori error analysis gives the opportunity to optimize the meshes used in the solution process and in the error control.To this end, we shortly describe the Dual Weighted Residual Method originally proposed in [1] and [2]. The method provides reliable a posteriori error estimators based on local residuals weighted by s...