Introduction
La loi sur l'air du 30 décembre 1996 oblige les zones urbaines françaises de plus de 100 000 habitants à surveiller, analyser et informer la population sur l'état de la qualité de l'air qu'elle respire. Durant les périodes estivales, la ville de Grenoble et ses alentours sont exposés à des périodes de fortes concentrations d'ozone (O3). L'Association pour le contrôle et la préservation de l'air dans la région grenobloise (ASCOPARG) a donc lancé un projet de modélisation du nuage photochimique estival en collaboration avec l'École
Polytechnique Fédérale de Lausanne (EPFL) sur trois ans. En 1999, la principale activité a été d'organiser la campagne GRENOble PHOTochemistry (GRENOPHOT) qui s'est déroulée de début juillet à mi-août sur l'ensemble de l'agglomération. La ville de Grenoble, qui se trouve au carrefour de trois vallées en forme de Y, plus communément appelé Y grenoblois, se démarque par ses spécificités dans la problématique de la qualité de l'air :
• épisodes fréquents de concentrations d'ozone dépassant la valeur horaire limite de 90 ppb ;
• topographie complexe du site avec des vallées étroites retenant les masses d'air ;
• circulation spécifique des masses d'air dans le Y grenoblois.
La production d'O3 est une fonction non linéaire dépendant des oxydes d'azote (NOx) et des composés organiques volatils (COV). La non-linéarité du processus de formation d'O3 en fonction de ses précurseurs implique qu'il n'existe pas de relation simple entre les émissions et l'O3 [1]. Pour comprendre sa formation, l'emploi d'un modèle eulérien à méso-échelle [2] est nécessaire.
Dans ce travail, le modèle METeorological PHOtochemistry MODel (METPHOMOD) [3], couplé au mécanisme chimique Reactiv Atmospheric Chemistry Mechanism (RACM) [4], simule la formation d'ozone sur le Y grenoblois. Ces résultats sont utilisés pour étudier la dispersion spatio-temporelle du panache des photo-oxydants et leur distribution verticale en fonction du régime de vent. Cette étude se situe dans le cadre de l'effort d'amélioration des performances des modèles photochimiques à l'aide de mesures représentatives 3D. Elle représente une première sur la région grenobloise d'une part par l'apport de mesures dans les trois vallées, d'autre part par leur utilisation pour valider le modèle METPHOMOD sur le cas de base, enfin pour mieux comprendre l'évolution des processus photochimiques en terrain complexe.
Campagne de mesure GRENOPHOT (25-26-27 juillet 1999)
Le but de la campagne GRENOPHOT est d'évaluer la production d'ozone et d'autres polluants dans le Y grenoblois situé au pied des Alpes, tant d'un point de vue expérimental que numérique. Cette campagne de mesures a vu se déployer , de début juillet à mi-août, un réseau de mesures (au sol et en altitude) renforcé dans le sud de l'agglomération, un système DOAS (Difterential Optical Absorption Spectroscopy) [5], deux LIDAR (Light lntensity Detection And Ranging) [6], un profileur de vent, et trois jours de mesures aéroportées réparties se sont ajoutées aux mesures habituelles du réseau ASCOP A RG (Figure 1, p. 280).
Glossaire
• ASCOPARG : Association pour le contrôle et la préservation de l'air en région grenobloise .
• ASL (Above Sea Level) : unité prenant le niveau de la mer comme référence.
• CITEPA : Centre interprofessionnel technique d'études de la pollution atmosphérique .
• CML : couche de mélange convective.
• COV : composés organiques volatils.
• DEM : Digital Elevation Model.
• DOAS : Differential Optical Absorption Spectroscopy.
• EPFL : École Polytechnique Fédérale de Lausanne.
• GRENOPHOT : GRENOble PHOTochemistry. http://lpas .epfl.ch/MOD/study/asco/
• GTOP030 : base de données de la topographie mondiale . http://edcdaac.usgs.gov/gtop030/gtop030.html
• ISM : Institut Suisse de Météorologie .
• LIDAR: Light Intensity Detection And Ranging.
• METPHOMOD : METeorological PHOtochemical MODel. http://www .giub.unibe.ch/klimetlmetphomod/
• NOx : somme des oxydes d'azote (NO + N02) .
• POI : période d'observation intensive .
• ppb (part per billion) : unité de concentration s'exprimant en milliardième de volume.
Cette unité est principalement utilisée lorsqu'oncompare les concentrations d'un même polluantà différentes altitudes . Dans un paque t d'air àmouvement vertical, la concentration d'un polluant, exprimée en ppb, ne varie pas, contrairement aux concentrations exprimées en microgrammes par mètre cube.
À 9° C et 950 hPa, 1 ppb O3 =1,95 Ilg/m3
À - 8° C et 653 hPa, 1 ppb O3 =1,42 Ilg/m3
• RACM : Regional Atmospheric Chemistry Mechanism.
• USGS/LCS: United States Geological Survey's/Land Cover System http://edcdaac.usgs.gov/glcc/glcc.html
La première période d'observation intensive (POI) a eu lieu du 25 au 27 juillet sous l'anticyclone des Açores bien installé sur les Alpes dans un flux de nord faiblissant. Le 26 juillet, à 17 heures, heure locale (HL), un fort panache d'03 atteignant des concentrations de 95 ppb (maximum mesuré au cours de GRENOPHOT) s'est développé à 20 km au sud de Grenoble. La seconde POI s'est déroulée du 1er au 3 août au cours d'une situation météorologique moins stable ayant généré des niveaux d'ozone inférieurs à 90 ppb. La première POI a été prise comme le cas de référence de la partie de simulation pour GRENOPHOT . Cette période ensoleillée a conduit à des niveaux élevés de photo-oxydants au sud de l'agglomération grenobloise [7].
Modélisation
Description du modèle
Le modèle employé pour cette étude est le modèle METPHOMOD [3]. C'est un modèle eulérien tridimensionnel utilisé pour le calcul d'épisodes de forte pollution en photo-oxydants dans les régions montagneuses. Il se compose d'un module météorologique pronostique :dynamique des fluides atmosphériques (champs de vents tridimensionnels, champs de pression hydrostatique ou non hydrostatique, densité de l'air, température et humidité) , turbulence atmosphérique (modèle k-t [8]), radiation solaire [9], équilibre du rayonnement infrarouge [10], transfert atmosphérique avec le sol et la végétation (implémentation d'un module en cinq couches) [11], nuages, pluie et neige. Il contient aussi un module chimique traitant la chimie en phase gazeuse, les émissions et les dépositions. METPHOMOD utilise un maillage cartésien avec une grille régulière.
Le module chimique en phase gazeuse utilisé est le module RACM. Il traite 237 réactions chimiques incluant 17 espèces inorganiques stables, 4 espèces inorganiques intermédiaires, 32 espèces organiques stables et 24 espèces organiques intermédiaires [4].
Mise en place : modélisation et nesting
La technique de nesting mise en place sur la ville de Grenoble implique un premier calcul sur un grand domaine, contenant le plus petit domaine, centré sur la zone d'étude. Les résultats obtenus sur le grand domaine sont alors utilisés comme condition de bord sur le petit domaine. Cette approche présente l'avantage de s'affranchir des conditions de bords pour la météorologie et la chimie (ce qui est important pour les photo-oxydants) sur le petit domaine de calcul.
Domaines de calcul
À l'horizontale, le grand domaine de simulation est un carré de 198 km de côté, orienté nord-sud. Il contient 33 x 33 cellules de 6 x 6 km. Ce domaine couvre l'ensemble de la région Rhône-Alpes et ses principales cités. Quant au petit domaine de simulation, il couvre un rectangle de 68 x 78 km de côté, centré sur Grenoble et orienté nord-sud, avec 34 x 39 cellules de 2 x 2 km. La figure 2, p. 281, montre les domaines de simulation et les isocontours de terrain.
À la verticale, un système de coordonnées cartésiennes avec 26 niveaux est utilisé pour le grand et le petit domaine. Les calculs météorologiques et chimiques sont effectués sur les mêmes grilles.
La base de données GTOP030 [12] a été utilisée pour les données de topographie. Il s'agit d'un modèle global digital d'élévation (DEM) avec une résolution de l'ordre de 1 x 1 km, rendu disponible par l'United States Geological Survey's (USGS) Earth Resources Observation System (EROS). Les données d'occupation des sols proviennent aussi de USGS EROS [13] avec une résolution de 1 x 1 km dans la classification USGS/Land Cover System (LCS). Dans le module de sol, 12 classes sont prises en compte et paramétrées selon le type de sol, l'albédo, la capacité calorifique du sol, la diffusivité du sol, l'humidité relative de l'air, la résistance d'évaporation du sol et des plantes.
Figure 1. Répartition spatiale des différents systèmes de mesures sur l'agglomération grenobloise au cours de la campagne de mesures intensives. (Reproduit avec l'autorisation de l'ASCOPARG).
Network of measurement systems over the Grenoble area during the intensive measurement campaign. (Reproduced with permission from ASCOPARG)
Figure 2. Domaines de simulation numérique pour GRENOPHOT
GRENOPHOT simulation domains.
Les émissions
Pour la campagne GRENOPHOT, deux inventaires d'émissions horaires ont été réalisés par le Centre interprofessionnel technique d'études de la pollution atmosphérique (CITEPA) [14] : le premier pour la grande grille avec une résolution de 6 x 6 km, et le second pour la ville de Grenoble avec une résolution de 1 x 1 km. Ces cadastres font la distinction entre trois types d'émissions : les émissions surfaciques (neuf secteurs), linéaires (autoroutes et routes nationales) et ponctuelles, selon deux types de jours (semaine et dimanche). Chacun de ces inventaires inclut les émissions de NOx, de méthane (CH4), de monoxyde de carbone (CO) et 23 spéciations de composés organiques volatils non méthaniques (COVNM) [15]. Ces émissions ont été assemblées suivant les 32 classes de COV [4] nécessaires à l'emploi du module chimique RACM. Sur le grand domaine de calcul, les émissions les plus fortes se localisent sur les villes de Lyon et de Grenoble ainsi que le long des autoroutes partant de ces deux grandes métropoles.
Conditions initiales et de bords
Calcul météorologique
Le calcul météorologique commence à 2 h 00 HL le 24 juillet. Il est initialisé en tout point du domaine à l'horizontale et pour chaque niveau vertical à partir des calculs effectués par le modèle météorologique de prévision de l'Institut Suisse de Météorologie (ISM) [16], pour les paramètres : vitesse du vent, température potentielle et humidité de l'air. Au cours de la simulation, les paramètres extraits du calcul ISM sont interpolés linéairement par METPHOMOD pour forcer les bords. Au sommet du modèle, la pression est fixée par les valeurs ISM et les autres paramètres sont soumis à l'hypothèse du gradient nul. Les prévisions du modèle de l'ISM montrent un régime de nord avec des vents de 4 à 5 m/s, faiblissant dans la troposphère libre le 24 juillet. À partir du 26 juillet, ce régime stable bascule progressivement vers un régime de vent de sud-ouest devenant fort dans la nuit du 27 au 28 juillet.
Calcul chimique
Il est important de disposer de bonnes conditions de bords pour l'03, les COV, le CO et les NOx. La condition de bord pour l'ozone est la plus délicate à fixer car on peut avoir plusieurs dizaines de ppb de différence entre deux épisodes de pollution. En utilisant la technique du nesting, l'influence de la région Rhône-Alpes et de ses grandes cités sur Grenoble est mieux prise en compte. Pour les quatre bords du grand domaine de calcul, les concentrations standard trouvées dans la troposphère libre ont été utilisées comme valeurs d'initialisation (Tableau 1). Le modèle commence le calcul 24 heures avant l'épisode pour atténuer l'influence des conditions initiales. Dans notre cas, il commence le 24 juillet à 2 h 00 HL.
Comparaison des résultats de simulation avec les mesures
Le modèle simule dans le temps la variation spatiale des photo-oxydants au cours de la POI du 27 juillet. L'ensemble de ces comparaisons portent sur la région de Grenoble et ses contreforts alpins, avec une analyse plus poussée dans la partie sud (Figure 3, ci-contre).
Météorologie
Au cours de la période d'observation, le temps synoptique sur la France a été dominé par l'anticyclone des Açores avec des vents d'altitude de nord faiblissant et tournant au sud-ouest la journée du 27 juillet. Les profils de vent réalisés avec le profileur de vent, à Vif, ont montré une stratification verticale de l'atmosphère [17] suivant trois couches : la couche synoptique située au-dessus de 2 500 à 3 000 rn, la couche thermique directement influencée par la topographie du site entre le sol et 1 200 à 2 400 m, et une couche intermédiaire, plus difficile à caractériser, variant entre 1 200 et 3 000 m.
Au cours de cette POI, la couche synoptique est caractérisée par un vent de nord avec une vitesse de 16 m/s le 25 juillet , faiblissant à 9 m/s le 27 juillet à midi en ayant tourné progressivement au sud-ouest, soit une bascule de 135 degrés. Les mesures réalisées avec le profileur de vent ont mis en évidence ces trois couches, reproduites par le modèle. La figure 4, p. 284, illustre la difficulté pour reproduire correctement la complexité des champs de vent dans une zone comme Grenoble. Néanmoins, ces comparaisons montrent la stratification verticale au niveau de la direction du vent ; elle est plus nuancée pour la vitesse et montre que le modèle, à ce moment-là, sous-estime la vitesse du vent au niveau de la première couche.
|
Espèces chimiques |
Concentrations (ppb) |
|
O3 |
60 |
|
NO |
0,1 |
|
NO2 |
0,1 |
|
HONO |
3,85 10-6 |
|
HNO3 |
9,12 10-6 |
|
HNO4 |
1,89 10-6 |
|
N2O5 |
3,75 10-6 |
|
NO3 |
1,91 10-7 |
|
HO2 |
2,19 10-7 |
|
CO |
4,40 102 |
|
HCHO |
3,0 |
|
ALD |
3,0 |
|
KET |
3,0 |
|
MGLY |
1,79 10-6 |
|
PAN |
0,2 |
|
MO2 |
1,16 10-6 |
|
ACO3 |
2,25 10-9 |
|
ONIT |
5,52 10-7 |
|
HC3 |
3,0 |
|
HC5 |
3,0 |
|
HC8 |
2,6 |
|
ETE |
1,5 |
|
OLT |
1,0 |
|
OLI |
0,6 |
|
TOL |
1,3 |
|
XYL |
1,5 10-3 |
|
DCB |
2,20 10-6 |
|
CSL |
1,91 10-7 |
|
H2O2 |
2,47 10-7 |
|
ISO |
2,0 10-11 |
|
CH4 |
1,7 103 |
Tableau 1. Concentrations de fond des espèces chimiques en phase gazeuse prises comme condition initiale et condition de bord sur le grand domaine de calcul.
Background of chemical species concentrations used for the big simulated domain.
Figure 3. Classification des stations de mesures pour les comparaisons par type. Les deux stations rurales sont aussi des stations de moyenne altitude (Saint-Dizier 1 000 m, Saint-Barthélemy 620 m).
Ground stations typology
Les résultats du modèle montrent une interaction complexe entre le vent à l'échelle synoptique et les systèmes de vents locaux , qui se développent dans les trois vallées du Y grenoblois. Les vents locaux se développent dans la couche thermique située entre le sol et 1 200 à 2 400 m. L'après-midi du 27 juillet, le modèle montre un vent de sud-ouest dans la troposphère libre et au sommet des massifs montagneux. Cependant, ce régime synoptique n'a pas d'influence sur les vents locaux de la couche thermique, dominé par les effets de canalisation des vallées qui transportent les masses d'air de la ville vers le sud de l'agglomération, et les vents de pente qui stockent les masses d'air sur les pentes de Belledonne avant de les pousser en fin de journée vers le sud (Figure 5, p. 285). La vitesse et la direction du vent simulées sont comparées dans le temps avec les mesures de quatre stations ponctuelles situées autour de la ville (Figure 6, p. 286). À Saint-Barthélemy, le cycle diurne (vents catabatiques-anabatiques) est bien mis en évidence par les mesures et la simulation. Pour la station de l'IGA, les résultats de simulation sont moins bons. Ce résultat est dû à la situation géographique de convergence de fond de vallée de cette station qui est une zone de turbulence. Des mesures traditionnelles de température ont été réalisées pour ces stations et les comparaisons avec les simulations montrent un bon accord (Figure 7, p. 287). Ces résultats montrent la capacité du modèle à retrouver le régime thermique et à simuler correctement les systèmes de vents de vallée, paramètres primordiaux pour les processus photochimiques.
Chimie
Distribution horizontale de l'ozone
La figure 8, p. 288, montre la distribution spatiale des concentrations d'O3 sur la couche de surface du modèle le matin et en fin d'après-midi. Le matin, la destruction de l'O3 en surface par titration avec le NO dans les zones émettrices de ce polluant précurseur est mise en évidence. L'après-midi, le panache d'O3
Figure 4. Comparaison des profils de vents, directions (degrés) et vitesses (m/s), réalisés par le profileur de vent et calculés par le modèle, à Vif, le 25 juillet à 12 h 00 HL.
Comparison between windprofiler and model calculalions for wind direction (°) and wind speed (m/s), at Vif, the 25th July at 12 h 00 Local Time (LT).
produit par les polluants précurseurs (NOx et COV) émis par l'agglomération est transporté au sud de la ville. C'est là qu'il atteint son maximum à environ 20 km du centre urbain. Les masses d'air coincées entre le massif du Vercors et le massif de l'Oisans, du fait des vents faibles, permettent aux polluants primaires de réagir ensemble pour former le maximum d'O3. Ensuite, ce panache d'O3 est poussé par les vents de vallée au sud de la ville, avant d'être dilué par le transport et la turbulence, et finalement alimenter la couche réservoir en début de nuit.
Pour faciliter les comparaisons des résultats de simulation avec les mesures, la figure 3, p. 283, suggère une classification des stations de mesures suivant trois types :
• type 1 rural : stations de Saint-Barthélemy et de Saint-Dizier (niveaux d'O3 élevés la nuit) ;
• type 2 semi-urbain : stations de Vif, Champ-sur-Drac, Voreppe et Versoud (niveaux d'O3 moyennement élevés la nuit) ;
• type 3 urbain : stations de Villeneuve et Fontaine (niveaux d'O3 nocturnes inférieurs à 10 ppb).
Les stations de type 2 et 3 se caractérisent par de faibles niveaux d'O3 nocturnes car elles sont proches des sources de pollution. Ces concentrations, basses durant la nuit, sont essentiellement dues à la destruction de l'O3 par les NO,( En revanche, les stations de type 1, loin des sources de pollution, affichent des niveaux d'O3 élevés la nuit car la déposition sèche est le processus de destruction de l'O3 le plus important. Dans la journée, les stations de type 3, soumises à de fortes émissions de NOx. affichent des maximums d'O3 assez faibles. Les fortes concentrations d'O3 indiquées par les stations de type 1 et 2 situées au sud de Grenoble identifient la zone où se développe le panache photochimique.
Le maximum d'O3, compris entre 90 et 100 ppb, apparaît entre 14 h 00 et 15 h 00 HL le 27 juillet dans la vallée du Drac où se trouvent les stations de Champ-sur-Drac, Vif et Saint-Barthélemy (Tableau 2, p. 287). Les stations de type 2 enregistrent 90 ppb d'O3 à 14 h 00 HL. À 16 h 00 HL, le maximum d'O3 de 95 ppb est capturé par Saint-Barthélemy, station de type 1, soit 1 heure après avoir enregistré le maximum d'O3 en fond de vallée (Figure 9, p. 289). Les concentrations nocturnes d'03 simulées sont dans l'ensemble en accord avec les mesures. Le modèle simule correctement les concentrations élevées d'O3 pour les stations de type 1. L'O3 stocké dans la couche réservoir est mis en évidence au cours des trois jours et atteint une concentration de 70 ppb durant la nuit du 26 au 27 juillet. Les faibles concentrations résultant des réactions avec les NOx sont bien reproduites en milieu urbain (Figure 8, p. 288).
Figure 5. Vents simulés dans la couche de surface du modèle, le 26 juillet à 9 h 00 et 17 h 00 HL, et le 27 juillet à 10 h 00 et 17 h 00 HL.
Wind in the lowest layer on 26th July at 9 h 00 and 17 h 00 LT, and 27th July al 10 h 00 and 17 h 00 LT.
Distribution horizontale des NOx
La figure 10, p. 290, montre que les NOx sont observés près des sources d'émissions (centres urbains et grands axes routiers). De ce fait, leur concentration est plus difficile à modéliser car les NOx vont être grandement influencés par des effets locaux mal considérés à méso-échelle. La figure 11, p. 291, montre un assez bon accord entre les mesures et les valeurs simulées pour les stations de type 1 et 2. Le stations de type 3 présentent des valeurs nocturnes simulées plus élevées que les mesures (Figure 9, p. 289). Le modèle sous-estime la hauteur de mélange la nuit sur les zones urbaines, car la formulation de calcul de la turbulence est mal
Figure 6. Comparaison des directions (degrés) et des vitesses de vent (mis) dans la couche de surface du modèle (trait plein) avec les mesures (cercles) des stations ponctuelles
Comparison between ground measurements and mode ! calculallons !rom 2s1n to 2 ?1h July for w1nd direction (0 and wind speed (mis). Solid lines are mode! results and circles are measurements.
Figure 7. Comparaison de la température (°C) dans la couche de surface du modèle (trait plein) avec les mesures (cercles) des stations ponctuelles
Comparison between ground measurements and model calculations from 25th to 27th July of temperature (°C). Solid lines are model results and circles are measurements.
Tableau 2. Concentrations maximales horaires d'03, de NO, de N02 et delormaldëhyde (HCHO). en ppb, mesurées et simulées à Vif le 26 et le 27 juillet dans le panache de la ville
Highest O3, NO, NO2 and HCHO concentrations measured and simulated in the Grenoble plume on 26th and 27th July
|
26-07-99 |
27-07-99 |
||||||
|
Espèces chimiques |
mesures |
simulations |
pic |
Espèces chimiques |
mesures |
simulations |
pic |
|
O3 |
79,4 ppb |
79,8 ppb |
oui |
O3 |
90,5 ppb |
88,3 ppb |
oui |
|
NO |
2,7 ppb |
2,21 ppb |
oui |
NO |
- |
0,72 ppb |
- |
|
NO2 |
8,6 ppb |
8,98 ppb |
oui |
NO2 |
7,4 ppb |
3,79 ppb |
oui |
|
HCHO |
5,9 ppb |
5,70 ppb |
oui |
HCHO |
7,2 ppb |
5,69 ppb |
non |
adaptée à la ville. Les polluants émis sont moins dilués que dans la réalité et restent alors prisonniers dans le premier niveau vertical du modèle [18].
Distribution verticale de l'ozone
Sur le site de Grenoble, les mouvements verticaux et le mélange des masses d'air vont conditionner le maximum d'O3. On a montré qu'au cours des épisodes photochimiques une stratification verticale du vent se met en place. C'est pourquoi on s'attend à trouver une stratification verticale pour la chimie.
La figure 12, p. 292, compare l'évolution verticale dans le temps de l'O3 mesurée par LIDAR à Vif et la simulation. La couche de mélange convective (CMC) monte d'un jour à l'autre car les températures augmentent entre le 25 et le 27 juillet avec peu d'humidité. Le modèle calcule une hauteur de mélange de 1 300 m le 25 juillet et de 2 100 m le 27 juillet. Les mesures LIDAR montrent une CMC un peu plus haute
Figure 8. Concentrations d'O3 simulées dans la couche de surface du modèle:(a) 25 juillet 17 h 00 HL, (b) 25 juillet 20 h 00 HL, (c) 27 juillet 9 h 00 HL, (d) 27 juillet 18 h 00 HL
Simulated O3 values [ppb] in the lowest layer on 25th July at (a) 17 h 00, (b) 20 h 00 and 27th July at (c) 9 h OO, and (d) 18 h 00 LT.
s'élevant jusqu'à 2 400 m le 27 juillet. Cette différence peut être due à un effet de la topographie complexe du site de Grenoble. Le modèle et le LIDAR mettent ainsi en évidence une CMC dont le sommet est plus haut que celui mesuré pour les régions de plaine [17]. La nuit tombée , la CMC est ramenée à sa valeur nocturne, et il se forme en altitude la couche résiduelle appelée aussi couche réservoir (Figure 12, p. 292) qui va stocker une partie de 1'O3 produit au cours de la journée dans la CMC. Entre le 26 et le 27 juillet, le modèle et le LIDAR mettent en évidence cette couche réservoir (19) entre 1 200 et 2 000 m avec un maximum d'O3 de 70 ppb. Dans la couche stable nocturne, entre le sol et 600 m, l'O3 est rapidement détruit en réagissant avec la molécule NO. Cette destruction commence dans la couche de surface à 2 h 00 HL, le 26 juillet, et s'accentue à partir de 5 h 00 jusqu'à 8 h 00 HL avec l'augmentation des émissions de NOx dues majoritairement au trafic. Entre 8 h 00 et 11 h 30 HL, le brassage de la couche réservoir avec la CMC explique l'augmentation des concentrations d'O3 proches du sol avant l'arrivée du panache produit par l'agglomération.
Figure 9. Comparaison des concentrations d'O3 (ppb) simulées (trait plein) dans la couche de surface du modèle avec les mesures (cercles) des stations ponctuelles
Comparison between ground measurements and mode ! simulations from 25th to 27th July for O3 (ppb). Solid lines are model results and circles are measurements.
Au niveau de la couche de surface, les profils verticaux d'03 affichent systématiquement une discontinuité entre la mesure au sol du DOAS et le premier point de mesures LIDAR (Figure 13, p. 293). Les simulations réalisées sans tenir compte du flux de déposition sèche pour l'03 montrent que cette discontinuité, vue par les mesures, est bien réelle et correspond à la déposition sèche. Ce dépôt sec peut être estimé à 10 ppb pour les 100 premiers mètres au-dessus du sol. La simulation montre que c'est vers 15 h 00 HL que le dépôt sec de 1'O3 est le plus élevé, c'est-à-dire quand le maximum d'O3 est atteint.
Figure 10. Concentrations de NOx simulées dans la couche de surface du modèle :(a) 25 juillet 11 h 00 HL,(b) 25 juillet 16 h 00 HL, (c) 27 juillet 9 h 00 HL, (d) 27 juillet 17 h 00 HL.
Modeled NOx values [ppb] on the surface on 25th July al (a) 11 h 00, (b) 16 h 00 and 27th July at (c) 9 h 00, and (d) 17 h 000 LT.
Le panache d'O3 produit par la ville est mis en évidence à Vif par une répartition verticale sur la hauteur de la CMC les après-midi des 25, 26 et 27 juillet. L'O3 généré par les NOx et les COV principalement émis au centre-ville, est transporté horizontalement par les vents thermiques poussant le panache vers le sud. Pour estimer la production de l'agglomération, un calcul a été effectué pour le petit domaine sans émissions anthropiques. Le tableau 3, p. 294, montre les maximums d'O3 et le ΔO3 entre le cas de base et le scénario zéro émissions pour le 26 et le 27 juillet à Saint-Barthélemy. Grâce à ces résultats, la production de l'agglomération peut être estimée à 32 ppb dans ce cas, soit un tiers de la valeur totale d'O3 mesurée et simulée.
Figure 11. Comparaison des concentrations de NOx simulées (trait plein) dans la couche de surface du modèle avec les mesures (cercles) à différentes stations ponctuelles
Comparison between ground measurements and model simulations from 25th to 27th July for NOx, Solid fines are model results and circles are measurements
Figure 12. Comparaison de l'évolution verticale temporelle des concentrations en O3 simulées et mesurées par le LIDAR, à Vif, entre le sol et 3 500 m asl. Le LIDAR a une résolution verticale de 10 m et une résolution temporelle de 6 min, moyennée sur 1 heure.
Measured and simulated vertical profile of O3 [ppb] between the 25th and the 27th July. Measured vertical O3 values [ppb] were performed by the EPFL mobile tropospheric ozone UV-DIAL LIDAR with a vertical resolution of 10 m and a temporal resolution of 6 min, averaged over 1 hour.
Conclusion
Les modèles eulériens pronostiques comme METPHOMOD utilisent un grand nombre de données en entrée. Les données météorologiques synoptiques et les concentrations chimiques de fond ne sont souvent pas disponibles et doivent être estimées, ce qui conduit à une incertitude pour les concentrations simulées. De plus, la résolution spatiale et temporelle des émissions devra être améliorée dans le futur pour générer moins d'erreurs et d'incertitudes.
METPHOMOD devrait aussi bénéficier d'un module pour la chimie hétérogène et les aérosols, et mieux prendre en compte la couche de surface et l'influence des nuages sur la radiation solaire.
En dépit de ces futures améliorations, les résultats de la campagne de mesures, associés à la modélisation sur la région de Grenoble, ont permis d'identifier une stratification verticale dynamique différente de la stratification verticale chimique lors d'épisodes photochimiques estivaux sur Grenoble. Les comparaisons spatiales et temporelles des simulations avec les mesures montrent que le modèle est capable de reproduire correctement la situation réelle et le panache au sud de la ville. La modélisation du cas de base a permis de valider le modèle sur un épisode de trois jours (25, 26 et 27 juillet 1999). L'accord qualitatif entre les mesures au sol disponibles et les résultats
Comparaison des profils de concentrations en O3 (ppb) réalisés avec le LIDAR (la première valeur de chaque profil est la valeur DOAS au sol) et calculés par le modèle, à Vif, le 26 juillet à différentes heures de la journée :(a) 9 h 00, (b) 14 h 00 et (c) 16 h 00 HL. Dans la figure bla mesure réalisée par avion est ajoutée el valide la mesure LIDAR
Comparison belween O3 LIDAR vertical profile using DOAS measurement for ground value and simulation on 261h July at Vif site :(a) 9 h 00.(b) 14 h 00 and (c) 16 h 00 LT. On the figure b, the airerait measurement validates the LIDAR.
|
Maximum O3 calculé (cas de base) |
ΔO3 maximum (cas de base-simulation zéro émissions) |
|
|
26-07-99 |
84,2 ppb |
32,5 ppb |
|
27-07-99 |
93,7 ppb |
31,8 ppb |
Tableau 3. Maximums d'O3 pour le cas de base et Δ03 entre le cas de base et le scénario zéro émissions pour le 26 et le 27 juillet
Highest O3 values between base case and Δ03 (base case-zero emissions) simulations on 26th and 27th July.
du modèle sont bons pour l'O3 et peuvent encore être améliorés pour les polluants primaires tels que les NOx, surtout en zone urbaine. Les profils verticaux LIDAR et la modélisation mettent en évidence une CMC supérieure à 2 100 m au cours de la journée , et la formation de la couche réservoir d'O3 entre 1 200 et 2 000 m pendant la nuit. Le maximum d'O3 est atteint le 27, car le régime de vents dans la cuvette grenobloise n'est pas influencé par la situation synoptique avec un vent soufflant du sud-ouest. Les résultats montrent que la ville génère environ un tiers du maximum d'03. C'est sur cette quantité d'03 produit en deux à trois jours au cours de périodes de smog estival qu'il faut agir, car elle amène aux dépassements des seuils réglementaires.
Sachant que le problème du smog estival peut seulement être résolu par une baisse massive des polluants précurseurs, et que les modèles permettent une interprétation quantitative de l'influence des émissions des polluants précurseurs. les bons résultats de METPHOMOD suggèrent son emploi pour le calcul de scénarios de réduction des émissions afin d'appuyer le développement de stratégies d'abattement pour l'ozone telles que l'évaluation des nouveaux Plans de déplacements urbains (POU).
Les auteurs tiennent vivement à remercier l'ASCOPARG pour avoir organisé ce projet avec le Laboratoire de Pollution de l'Air de l'EPFL et pour son soutien financier. Nous remercions aussi l'ensemble du personnel ASCOPARG qui a participé à la campagne de mesures pour son enthousiasme et sa disponibilité ainsi que les sociétés DEGREANE et METAIR.