Local scale air quality model system for diagnostic and forecasting simulations using the finite element method

Air pollution is an important topic with a great social impact; it is related with public health, environment and ecology, and climate change. Scientists have developed several models in the last thirty years, and regional air quality operational systems are used routinely by governments and agencies. Efforts have also been done to simulate the air quality in the local scale; main models are Gaussian and Puff models, that are based on a Lagrangian approach. In contrast with these models, in this thesis we have developed a system using an Eulerian approach. This model is specifically designed for regions with complex orography where the Lagrangian models have problems computing the trajectory of the particles. This model can be used for diagnostic or prediction simulations. Air quality operational systems depend on the orography, meteorological data, and emission data. Air quality models use processors to incorporate these data into the model. The data can come from numerical weather prediction systems, experimental data, or databases. In this thesis we have developed processors, specifically designed for the local scale, to incorporate these data into our system. To incorporate the orography, we have developed a mesh generation algorithm suitable for complex terrain discretization; it also allows to insert layers that can match the regional models. A wind field model has also been used; it can interpolate a three-dimensional wind field from some station measurements using a log-linear vertical wind profile, or can interpolate it from a numerical weather prediction system. Once an interpolated wind field is computed, a mass-consistent model is applied to ensure null divergence and impermeability in the terrain. The wind field is modified to take into account the injection of the pollutants into the atmosphere. Briggs studied the trajectory of the plume rise giving some empirical equations that will be used in our model. Briggs' equations describe the trajectory in a plane; our model will modify this trajectory adapting it to the ambient wind field. This modification allows the plume rise to surround the mountains or channel into the valleys. The transport and reaction of pollutants in the atmosphere is then computed using an splitting method, so the transport and the chemical reactions are computed independently. To solve the transport of pollutants we have used a finite element method stabilized using least squares. The chemical reaction is simulated using simplified models such as RIVAD, or more complex ones such as CB05. To obtain more accurate results we have used adaptation. An error indicator has been used to adapt the mesh to the solution. To adapt the mesh to the concentration distribution of all the species is very demanding, for this reason we have used a multimesh method where every chemical specie has its own mesh where we solve the transport and the chemical reactions are simulated in a common mesh. The system developed in this thesis has diagnostic and forecasting capabilities. For this reason we present two different applications. The first one is a diagnostic application in La Palma island (Spain), where wind measurements are given, and SO2 and NO2 emissions from a stack are considered. The topography of the island is real, from a digital elevation model, but the wind field measurements, and the stack location and emissions, are simulated. The second application is a forecasting application data from the CMAQ benchmark test. It is located in the surrounding of Pineville Kentucky. In this application we have used all the data from CMAQ and the chemical reaction model CB05.La contaminació atmosfèrica té gran impacte social; està relacionada ambla salut pública, l'ecologia, el medi ambient i el canvi climàtic. En els últims trenta anys, els científics han desenvolupat diversos models que els governs i les agències mediambientals utilitzen diàriament. També s'han fet esforços per simular la qualitat de l'aire en l'escala local. Els models principals són els models Gaussians i Puff, que es basen en una descripció Lagrangiana. En contrast amb aquests models, en aquesta tesi s'ha desenvolupat un sistema que utilitza una descripció Euleriana. Aquest model està dissenyat específicament per a les regions amb topografia complexa on els models Lagrangians tenen problemes calculant la trajectòria de les partícules. El model presentat en aquesta tesi pot ser utilitzat tant pel diagnòstic com per la predicció. Els sistemes operatius utilitzats actualment depenen de l'orografia, de dades meteorològiques i de dades d'emissió. Per tal d'incorporar aquestes dades, els models de qualitat de l'aire han desenvolupat diferents preprocessadors. Les dades poden venir dels sistemes numèrics de predicció meteorològics, de dades experimentals o de bases de dades. En aquesta tesi hem desenvolupat preprocessadors dissenyats específicament per a l'àmbit local, per tal d'incorporar aquestes dades al nostre sistema. Per incorporar l'orografia, hem desenvolupat un algoritme de generació de malles adequat per terreny complex; l'algoritme també ens permet inserir capes que poden coincidir amb la dels models regionals. S'ha desenvolupat un model de càlcul de camp de vent; a partir de les dades s'interpola un camp de vent tridimensional fent servir un perfil log-lineal vertical del vent, o s'interpola a partir dels resultats d'un sistema de predicció meteorològica. Quan s'ha calculat el vent interpolat, es fa servir un model de massa consistent per obtenir un camp de vent definitiu on s'ha imposat divergència zero i impermeabilitat del terreny. El camp de vent s'ha de modificar per tenir en compte la injecció dels contaminants a l'atmosfera. Briggs va estudiar la trajectòria de l'elevació de la ploma i va trobar equacions empíriques que utilitzarem en el nostre model. Les equacions de Briggs descriuen la trajectòria de la ploma en un pla vertical; el nostre model modificarà aquesta trajectòria adaptant-la al camp de vent ambiental. Aquesta modificació permet que l'elevació de la ploma rodegi les muntanyes o s'acanali a les valls. El transport i la reacció de contaminants a l'atmosfera es calcula utilitzant un mètode de "splitting", de manera que el transport i la reacció química es calculen de forma independent. Per resoldre el transport de contaminants, hem utilitzat el mètode d'elements finits estabilitzat amb mínims quadrats. La reacció química es simula mitjançant models simplificats com el model RIVAD, o més complexes com el CB05. Per obtenir resultats més precisos hem adaptat la malla a la solució utilitzant un indicador d'error. Haver d'adaptar la malla a la distribució de la concentració de totes les espècies contaminants és molt exigent i per aquest motiu hem utilitzat un mètode "multimesh" on cada espècie química té la seva pròpia malla on resolem el transport i les reaccions químiques es simulen en una malla comú. El sistema desenvolupat en aquesta tesi té capacitats de diagnòstic i pronòstic. Per aquesta raó es presenten dues aplicacions diferents. La primera és una aplicació de diagnòstic a l'illa de la Palma (Espanya), on es disposa de mesures de vent, i de les emissions de SO2 i NO2 d'una xemeneia. La topografia de l'illa és real, a partir d'un model digital del terreny, però les dades del camp de vent, la ubicació de la xemeneia i el valor de les emissions són simulades. La segona aplicació és una aplicació de predicció fent servir els resultats del benchmark del CMAQ. Es simula una zona a Pineville Kentucky. En aquesta aplicació hem utilitzat totes les dades del CMAQ i el model químic CB05.Postprint (published version