Meta-modeling of a simulation chain for urban air quality

Abstract

Urban air quality simulation is an important tool to understand the impacts of air pollution. However, the simulations are often computationally expensive, and require extensive data on pollutant sources. Data on road traffic pollution, often the predominant source, can be obtained through sparse measurements, or through simulation of traffic and emissions. Modeling chains combine the simulations of multiple models to provide the most accurate representation possible, however the need to solve multiple models for each simulation increases computational costs even more. In this paper we construct a meta-modeling chain for urban atmospheric pollution, from dynamic traffic modeling to air pollution modeling. Reduced basis methods (RBM) aim to compute a cheap and accurate approximation of a physical state using approximation spaces made of a suitable sample of solutions to the model. One of the keys of these techniques is the decomposition of the computational work into an expensive one-time offline stage and a low-cost parameter-dependent online stage. Traditional RBMs require modifying the assembly routines of the computational code, an intrusive procedure which may be impossible in cases of operational model codes. We propose a non-intrusive reduced order scheme, and study its application to a full chain of operational models. Reduced basis are constructed using principal component analysis (PCA), and the concentration fields are approximated as projections onto this reduced space. We use statistical emulation to approximate projection coefficients in a non-intrusive manner. We apply a multi-level meta-modeling technique to a chain using the dynamic traffic assignment model LADTA, the emissions database COPERT IV, and the urban dispersion-reaction air quality model SIRANE to a case study on the city of Clermont-Ferrand with over 45, 000 daily traffic observations, a 47, 000-link road network, a simulation domain covering 180,text {km}^2. We assess the results using hourly NO_2 concentration observations measured at stations in the agglomeration. Computational times are reduced from nearly 3 h per simulation to under 0.1 s, while maintaining accuracy comparable to the original models. The low cost of the meta-model chain and its non-intrusive character demonstrate the versatility of the method, and the utility for long-term or many-query air quality studies such as epidemiological inquiry or uncertainty quantification.