1En Bourgogne, durant la période estivale, comme dans la plupart des régions françaises, les indices de qualité de l’air élaborés et communiqués par les AASQA (Associations Agréées de Surveillance de la Qualité de l’Air) sont majoritairement influencés par les concentrations atmosphériques en ozone. En effet ce polluant estival, à l’exception des abords des zones industrielles, est celui qui a le plus d’effet sur la santé humaine (Mégie, 1993). Ainsi, la prévision des concentrations en ozone est devenue un enjeu important en terme de santé publique.
2Pour effectuer ces prévisions, il faut prendre en compte le fait que l’ozone est un polluant secondaire issu de processus complexes. Interviennent dans sa formation des précurseurs (oxydes d’azote (Nox) et Composés Organiques Volatils Non Méthaniques (COVNM)) transformés lors de conditions météorologiques favorables (température et ensoleillement élevés, vents faibles). Ainsi, en Scandinavie, la diminution des émissions de précurseurs amorcée depuis 1990 commence à avoir des effets positifs sur la pollution à l’ozone (Solberg et al., 2005).
3Une collaboration entre l’université de Bourgogne, Météo-France et ATMOSF’air a été initiée à cette intention. L’année 2003 a servi de période test pour douze stations situées dans quatre agglomérations (Dijon, Chalon-sur-Saône, Macon et Montceau-les-Mines) ainsi qu’à Saint-Brisson, station rurale du massif du Morvan (Figure 1). L’approche par station permet de confronter prévisions et observations et de prendre en compte les différences importantes observées quant aux concentrations d’ozone entre centres et périphéries des villes (Bavay et Roussel, 1992 ; Richard et al., 2000 ; Jimenez et al., 2005). La période se limite à la saison chaude et ensoleillée (1er avril – 30 septembre) pour laquelle les concentrations d’ozone sont maximales, en particulier quand prédominent des conditions anticycloniques (Bavay et Roussel, 1992 ; Mégie et al., 1993 ; Richard et al., 2000).
Figure 1 : Stations ATMOSF’air, Météo France et mailles ALADIN.


4Les différentes prévisions de l’ozone peuvent être effectuées selon deux familles de méthodes :
51) purement numérique : les prévisions sont issues de modèles de prévisions météorologiques comprenant un module chimie-transport. Ces approches sont développées dans plusieurs régions françaises, par exemple pour l’île-de-France (Vautard et al., 2000) ou la Provence (Durand et Cros, 2004) ;
62) statistico-dynamique : des outils statistiques sont développés en aval de modèles numériques (en l’occurrence des modèles de prévisions météorologiques) intégrant ou non la chimie. La composante statistique de ces modèles s’appuie sur des régressions multiples, des analyses discriminantes, des réseaux de neurones, des méthodes de recherche de cas passés similaires ou des arbres de décisions (Fromage et Gilibert 1997 ; Ghiaus, 2005 ; McMillan et al., 2005).
7La méthode statistico-dynamique, connue et testée en météorologie opérationnelle sous la dénomination d’adaptations statistiques de sortie de modèle, est la base retenue pour cette étude de performance des prévisions locales de concentrations d’ozone.
8Le prédictand est la valeur à prévoir: la concentration maximale d’ozone du lendemain (J + 1). Pour chaque jour, nous avons sélectionné la valeur maximale d’ozone entre 9 h et 14 h UTC, tranche horaire où se situent généralement les valeurs les plus élevées (Houzé, 1999). Ce rythme diurne des concentrations en ozone est associé à l’insolation mais aussi à l’épaisseur de la couche limite (Bavay et Roussel, 1992). En outre, ce choix horaire permet une prévision à partir de 16h légales en été (14 h UTC), horaire compatible avec le souci opérationnel de pouvoir déclencher, le cas éventuel, les procédures d’information des autorités et de la population.
- 1 Les COVNM (Composés Organiques Volatiles Non Méthaniques) ne sont pas mesurés en Bourgogne. Dans (...)
9Nous ne présenterons ici que les prédicteurs chimiques mesurés par ATMOSF’air1 dans les 12 stations (Figure 1) et sélectionnés par les modèles :
101) l’ozone du jour (J) permet de prendre en compte un éventuel phénomène d’accumulation (Vergez, 1999) ;
112) les oxydes d’azote (NOx) constituent l’un des principaux précurseurs de la pollution photochimique (Jacob, 2000).
12Pour chacun, deux valeurs sont considérées :
131) la moyenne 9 h - 14 h UTC du jour J, qui minimise l’impact des variations rapides (haute fréquence) ;
142) le maximum 9 h - 14 h UTC du jour J, étant donné qu’un changement de temps rapide en cours de journée peut parfois rendre la moyenne peu pertinente.
15Ils sont tous issus ou dérivés d’ALADIN, version à domaine limité du modèle ARPEGE de Météo-France. ALADIN-France permet de gagner un facteur 2 par rapport à la résolution spatiale la plus forte d’ARPEGE auquel il est couplé (Pailleux et al., 2000). La version utilisée est donc la base opérationnelle des prévisions atmosphériques, lancée chaque nuit à 00h UTC. Les simulations sont archivées selon des échéances séparées de 3 h. Nous en avons retenues deux : 27 et 36 heures qui permettent de documenter : 1) la fin de nuit (échéance 27,03 h UTC., soit 5 h légales en été, moment où les températures comme les concentrations d’ozone sont minimales) ; 2) le début d’après midi (échéance 36, 12 h UTC., soit 14 h légales, moment de la journée ou températures et ozone atteignent leur valeurs maximales).
16Nous utilisons, pour trois points de grille centrés sur Dijon, Chalon-sur-Saône et Saint-Brisson (Figure 1), des valeurs prévues au niveau du sol (température, pression, flux solaires et thermiques, précipitations liquides de convection et précipitations liquides à grande échelle) ainsi que sur 11 niveaux (2, 10, 20, 50, 100, 250, 500, 750, 1 000, 1 250, 1 500 m) au-dessus du sol (température, composantes horizontales et verticales du vent, humidité relative). Enfin, nous intégrons la nébulosité du sol à 1 500 m, de 1 500 m à 5 km, au-delà de 5 km ainsi que totale.
17À partir de ces variables, nous en avons calculé l’amplitude diurne, ou encore la différence entre deux niveaux d’altitude afin de prendre en compte la stratification verticale dans les basses couches. L’ensemble représente au total 170 prédicteurs. A noter des valeurs manquantes dans les données qui interdisent les prévisions (cas de la fin juin 2003 par exemple). Pour Dijon et Chalon-sur-Saône, la pertinence des prévisions météorologiques est validée par les valeurs réellement observées aux stations de Dijon-Longvic pour le point de grille « Dijon », Chalon-Champforgeuil et Rully (rayonnement global) pour le point de grille « Chalon » (Figure 1).
18La régression linéaire multiple consiste à expliquer un phénomène (ici les concentrations maximales d’ozone à J + 1), par la conjonction des variations de plusieurs phénomènes. Les variables sont sélectionnées pas à pas (stepwise). Il s’agit d’un algorithme de régression utilisant les coefficients de corrélation partielle.
19Les équations statistiques utiles à la prévision doivent être construites avec des données utilisées pour caler les paramètres qui diffèrent des données servant à évaluer les résultats. On sépare les observations en jeu d’apprentissage d’une part (1999-2002) et jeu de validation d’autre part (2003). Sur la période d’apprentissage, nous avons utilisé la procédure de test dite Jack-knife, s’apparentant à la validation croisée dans la mesure où l’on retire tour à tour chaque observation de l’échantillon d’apprentissage. Ceci permet une validation à l’aide de données différentes de celles qui ont servi à l’ajustement de la droite. Cette méthode permet d’évaluer la performance réelle d’un modèle de prévision avant qu’il ne soit appliqué (Thiria et al., 1997). Nous obtenons, sur la période d’apprentissage (1999-2002), et ce pour chaque station, une équation de régression sélectionnant un maximum de 5 prédicteurs parmi les 170 soumis, sachant que ces 5 prédicteurs peuvent varier selon la station. Ce nombre maximal imposé privilégie le dynamisme des modèles sur leur performance en termes de corrélations.
20Les séries de prévisions ainsi obtenues ont des variances trop faibles, ce pour l’ensemble des stations ATMOSF’air. Pour pallier ce problème, classique en adaptation statistique de sortie de modèle météorologique, le correctif K (Ward et Folland, 1991) renforce le « dynamisme » des prévisions car augmente l’écart type de la prévision en l’ajustant à celui de l’observation de la période d’apprentissage (1999-2002).
21À partir des températures horaires relevées dans les stations de Météo-France (Dijon-Longvic ou Chalon-Champforgeuil) nous avons calculé la « climatologie » de la période d’apprentissage en effectuant une moyenne sur les 4 ans (Figure 2a). Celle-ci se caractérise par l’addition du cycle diurne (en ordonnées) et du régime saisonnier (en abscisses) sur lesquels se surimposent des fluctuations météorologiques indépendantes du calendrier et de l’heure. Ces fluctuations ne sont pas totalement lissées en raison de la brièveté de l’échantillon (4 ans). Par exemple, le hasard a fait, qu’à Chalon-Champforgeuil, la première décade d’août fut plutôt fraîche lors des étés 1999-2002 (Figure 2a).
22Les températures relevées sur la période de test (Figure 2b) présentent les mêmes caractères avec toutefois :
231) des fluctuations horaires et surtout d’un jour à l’autre plus marquées (il ne s’agit pas d’une moyenne) ;
242) des températures globalement bien plus élevées, caractéristiques du printemps et de l’été 2003 sur une bonne part de l’Europe de l’ouest (Grazzini et al., 2003).
25La période dite « de canicule » s’individualise bien avec des températures dépassant chaque jour 35 °C entre les 3 et 13 août, 38 °C entre les 5 et 13 août (excepté le 8) et culminant à 39,2 °C le 12 août (Figure 2b).
Figure 2 : Températures horaires en avril-septembre à Châlon-Champforgeuol (en °C).
a : période d’apprentissage : moyenne 1999-2002.
b : période de test : 2003.


26Les mêmes calculs sont effectués sur les valeurs d’ozone relevées dans 12 stations ATMOSF’air. Nous présentons ceux de Chalon-Champforgeuil, station située à proximité d’une station météorologique. Sur la période d’apprentissage (Figure 3a), les concentrations moyennes d’ozone sont généralement inférieures à 40 µg/m3 la nuit. Elles sont comprises entre 80 et 120 µg/m3 lors des après-midi d’été. Ces valeurs relativement faibles signifient que la pollution est globalement faible. Ceci n’exclue pas des épisodes marqués qui, n’intervenant pas les mêmes jours d’une année à l’autre, n’apparaissent pas de matière nette sur une moyenne. Cette forte variabilité inter annuelle explique que sur un échantillon court (4 ans) les fluctuations soient fortes d’un jour à l’autre. Ainsi n’y a-t-il par exemple pas d’autre interprétation que la brièveté de la période d’apprentissage pour expliquer les faibles concentrations observées vers le 15 juillet (Figure 3a).
27Le printemps et l’été 2003 (Figure 3b) se caractérisent par la présence de plusieurs séquences où les concentrations en ozone sont élevées (fin avril, première quinzaine de juin, mi-juillet, première quinzaine d’août, ...). Ces épisodes d’ozone sont à mettre en relation avec les conditions météorologiques exceptionnelles (fort ensoleillement et température élevée). Ainsi, lors de la canicule (dimanche 3 au mercredi 13 août), la valeur seuil légale d’information et de recommandation de 180 µg/m3 a été dépassée lors de 6 journées à la station de Chalon-Champforgeuil (6, 7, 9, 10, 12 et 13 août).
Figure 3 : Concentrations horaires en avril-septembre à Chalon-Champforgeuil (en µg/m3).
a : période d’apprentissage : moyenne 1999-2002.
b : période de test: 2003.


28Du point de vue du prévisionniste, ces conditions, tant d’ozone que météorologiques, signifient que la période de test comporte des valeurs inconnues en période d’apprentissage. L’enjeu est bien de savoir si, dans de telles conditions, les modèles statistico-dynamiques, à l’image des modèles dynamiques (Vautard et al., 2005) sont susceptibles de fournir des prévisions réalistes, gage de leur robustesse.
29Nous détaillons trois équations dédiées à des stations « périurbaines » de Saône-et-Loire :
30Chalon-Champforgeuil: O3max Pré = 44.5 + (0.46 O3max) + (4.41 T1) – (2.84 t1) – (O.36 H11) + (1.91 10-6 F_Th)
31Macon-Champlevert: O3max Pré = 55.7 + (0.44 O3moy) + (4.49 T1) – (2.72 t1) – (O.37 H11) + (1.9 10-6 F_Th)
32Montceau-Ecluse: O3max Pré = 73.6 + (0.52 O3max) + (1.28 T1) – (O.45 H8) – (1.07 U10) + (1.91 10-6 F_Th)
33Avec:
34O3max Pré : concentration maximale d’ozone prévue pour le lendemain (J+1)
35O3max: concentration maximale mesurée à la station par ATMOSF’air le jour J entre 9 h et 14 h UTC
36O3moy: concentration moyenne mesurée à la station par ATMOSF’air le jour J entre 9 h et 14 h UTC
37T1: température ALADIN échéance 36 h, soit J + 1 12 h UTC, à 2 m
38t1: température ;ALADIN échéance 27 h, soit J + 1 3 h UTC, à 2 m
39H11: humidité ALADIN échéance 36h, soit J + 1 12 h UTC, à 1 500 m
40H8: humidité ALADIN échéance 36h, soit J + 1 12 h UTC, à 750 m
41F_Th: flux thermique ; ALADIN échéance 36h, soit J + 1 12 h UTC
42U10 : composante zonale du vent ALADIN échéance 36 h, soit J + 1 12 h UTC, à 1 250 m
43Les prédicteurs sélectionnés varient peu d’une station à l’autre. Deux types de prédicteurs n’entrent que dans certains modèles : le vent et la température nocturne et quatre types de prédicteurs sont systématiquement pris en compte :
441) les concentrations en ozone du jour ;
452) la température à 2 m prévue pour le lendemain à 12 h UTC ;
463) l’humidité relative de l’air prévue à 12 h UTC (à différents niveaux d’altitude selon la station) ;
474) le flux thermique prévu à 12 h UTC.
48La présence de l’ozone du jour dans les modèles traduit l’inertie des concentrations d’un jour à l’autre (auto corrélation de rang 1) tout particulièrement en phase d’accumulation. Celle de la température et du rayonnement à 12 h UTC. prend en compte les conditions de formation de l’ozone. Les coefficients affectés à ces trois types de prédicteurs sont positifs. Aux autres prédicteurs (humidité, vent, température nocturne) sont associés des coefficients négatifs. Les phénomènes de destruction, de dépôt et de dispersion de l’ozone sont ainsi également pris en compte, mais potentiellement de manière un peu moins efficace
49Plusieurs indices statistiques donnent une mesure de la performance de ces modèles. Nous avons retenu le coefficient de détermination (R2) qui indique la part commune de variance des concentrations maximales d’ozone observées et prévues (Figure 4). On retiendra que les performances sont globalement très bonnes (minimum de 0.63 pour Dijon-Pasteur et maximum de 0.78 pour les deux stations mâconnaises) et qu’elles diffèrent peu d’une station à l’autre.
Figure 4 : Coefficients de détermination des 12 modèles.


50Pour Chalon-Champforgeuil, une des stations où le R2 est le moins élevé, les grandes phases (épisodes comme jours à faibles concentrations) sont globalement bien prévues (Figure 5). De surcroît les niveaux prévus de pollution sont conformes aux observations. Ce second point est fondamental et très encourageant car il montre l’aptitude de notre démarche à prévoir les valeurs extrêmes avec de bons ordres de grandeur : les pics de pollution en ozone sont donc prévisibles quantitativement par des méthodes statistico-dynamiques, ce qui était loin d’être un résultat acquis au départ, la prévision des extrêmes étant toujours la plus délicate.
Figure 5 : Concentrations quotidiennes maximales en ozone observées et prévues pour Chalon-Champforgeuil (2 avril - 30 septembre 2003)


Noir fin: observé ; Grisé épais : prévu
51Sur les différentes stations l’erreur moyenne de prévision est de l’ordre de 20 µg/m3. Dans le détail, le recouvrement entre séries observée et prévue est imparfait. Un zoom effectué sur la canicule et les jours qui l’encadrent montre un décalage : les prévisions ont tendance à suivre les observations avec un jour de retard, c’est très net les (1, 3, 4, 8, 11 et 14 août (Figure 6). Pour tester l’inertie induite par le prédicteur « ozone du jour », nous avons calculé les corrélations entre les deux courbes (observé et prévu) avec un décalage d’un jour. Les résultats montrent que les corrélations entre l’ozone prévu à J + 1 et l’ozone observé sont de même ordre selon que l’on considère l’ozone du lendemain (J + 1) ou l’ozone du jour (J). Cela signifie que le prédicteur « ozone du jour » joue un rôle déterminant dans la prévision et restreint la réactivité des modèles. Ce problème peut théoriquement être important dans les situations de changement de temps rapides (c’est le cas du 14 août), malgré l’introduction du facteur correctif K (cf. IV-B).
Figure 6 : Modèles de prévision avec et sans ozone du jour comme prédicteur (Chalon-Champforgeuil, 30 juillet au 17 août 2003)


Noir fin: observé ; Grisé épais : prévu ; Barbules : prévu sans ozone du jour.
52Pour discriminer le rôle de l’ozone du jour dans les prévisions nous avons refait tourner les modèles en supprimant ce prédicteur et ne conservant que les quatre prédicteurs issus des prévisions météorologiques. Nous sommes donc dans une situation virtuelle d’une atmosphère sans mémoire du point de vue de la chimie de l’ozone. Le déphasage a disparu. Il est donc imputable à la trop grande inertie liée à l’ozone du jour. En revanche, ce modèle se caractérise par une perte de variance expliquée : R2 chute de 0.62 à 0.51 pour Chalon-Champforgeuil, et une diminution comparable est observée pour les autres stations. Il s’avère également incapable de prévoir des valeurs réalistes, comme en témoignent les sous-estimations systématiques associées aux épisodes d’ozone où des phénomènes d’accumulation se développent, en particulier en fin d’épisode. Supprimer le prédicteur « ozone du jour », se traduit ainsi par une dégradation de la prévision les 31 juillet, 7, 9, 10, 12 et 13 août (Figure 6). Il est donc indispensable de considérer les concentrations en ozone du jour pour prévoir celles du lendemain. Ceci est cohérent compte tenu de la chimie de l’ozone, qui, du fait de la complexité des réactions peut se former pendant plusieurs jours à la suite d’émissions de précurseurs (Derwent et al., 2005).
53Une autre erreur systématique peut être liée aux différences de concentrations, en précurseurs de l’ozone et en ozone, observée entre les jours ouvrés et les week-ends et jours fériés. Ces différences se traduisent généralement par des concentrations en ozone anormalement élevées dans les villes les week-ends (Jimenez et al., 2005). Ce phénomène est observé en Bourgogne (Houzé, 2000) mais n’est pas suffisamment fort pour que cela soit judicieux de l’intégrer dans les modèles.
54Les autres erreurs de prévisions ont un caractère non systématique. Il s’agit d’erreurs isolées (dans le temps) mais pouvant être importantes (les plus fortes atteignent 60 µg/m3 en positif comme en négatif). De juillet à septembre, des erreurs de prévisions importantes, à savoir supérieures à |30| µg/m3, interviennent lors de 18 journées sur 92 (soit 20 % environ). Il s’agit de 12 sous-estimations (12, 16, 23, 29 juillet, 2, 7, 10, 13, 19, 27 août, 13, 21 septembre) et de 6 surestimations (2, 24 juillet, 15, 20 août, 3, 5 septembre). Nous en analysons une de chaque: la surestimation pour le jeudi 14 août et la sous-estimation pour le samedi 16 août. Lors de ces deux journées, les modèles statistico-dynamiques donnent de piètres résultats dans la plupart des 12 stations. Ceux de Chalon-Champforgeuil (tableau 1) et de Dijon-Balzac (tableau 2) révèlent l’ampleur des erreurs (de l’ordre de 50 µg/m3). Pour tenter d’expliquer ces erreurs, il est logique de se pencher sur les paramètres météorologiques prévus en amont et ;nous confrontons les observations météorologiques effectuées dans les stations Météo-France de Chalon-Champforgeuil et de Dijon-Longvic d’une part, avec les prévisions ALADIN sur les points de grilles correspondants.
Tableau 1 : Chalon-sur-Saône : observations et prévisions pour les 14 et 16 août 2003


T° (°C) et humidité relative (%) à 12 h UTC à Chalon-Champforgeuil. Rayonnement global (W/m²) cumulé depuis le lever du jour à 12 h UTC à Rully. T°, humidité et rayonnement global prévu à 12 h UTC sur maille ALADIN Chalon-sur-Saône. Ozone (µg/m3) prévu et observé à Chalon-Champforgeuil.
Tableau 2 : Dijon : observations et prévisions pour les 14 et 16 août 2003


T° (°C) et humidité relative (%) à 12 h UTC à Dijon-Longvic. Rayonnement global (W/m²) cumulé depuis le lever du jour à 12 h UTC à Dijon-Longvic. T°, humidité et rayonnement global prévu à 12 h UTC sur maille ALADIN Dijon. Ozone (µg/m3) prévu et observé à Dijon-Longvic.
55Les prévisions pour le jeudi 14 août sont largement surestimées. Ceci s’explique 1/ par l’inertie du modèle (à Chalon-Champforgeuil (Figure 6) les concentrations record enregistrées le mercredi 13 août (217 _g/m3) tirent la prévision vers le haut) et 2/ par l’erreur de prévision météorologique (à Chalon-sur-Saône comme à Dijon le temps est nettement moins chaud et ensoleillé que prévu (tableaux 1 et 2)). Les prévisions pour le samedi 16 août sont, a contrario, largement sous-estimées. La journée est nettement plus chaude et ensoleillée que prévu. L’humidité relative est bien plus faible, à Chalon-sur-Saône (tableau 1) comme à Dijon (tableau 2). Notons que ces erreurs météorologiques a priori importantes ne sont pas si rares. Il s’agit dans la plupart des cas d’erreurs associés à un changement de temps qui intervient quelques heures plus tôt (cas de la dégradation du 14 août) ou plus tard (cas de l’amélioration du 16 août) que prévu. Toutefois, cette analyse n’est pas surprenante, car les méthodes statistiques de sortie de modèle dynamique ne peuvent bien sur pas gommer les erreurs ou approximations des prévisions atmosphériques.
56Enjeu important en matière d’effets sur la santé, la prévision des pics de concentrations d’ozone dans l’air ambiant a bénéficié en Bourgogne d’un projet ambitieux de recherche en collaboration entre l’université, Météo-France et ATMOSF’air. Basées sur une méthode d’adaptation statistique de sortie du modèle atmosphérique ALADIN, les prévisions de concentrations maximales journalières d’ozone donnent des résultats très encourageants. Compte tenu de la difficulté classique à prévoir des valeurs extrêmes (qui sont en matière de concentrations polluantes les seules réellement intéressantes de manière à favoriser la prévention) et du caractère tout à fait exceptionnel des conditions atmosphériques de la période de validation (printemps-été 2003, considérée actuellement comme la canicule référence), il est rare et donc remarquable, pour un premier essai, d’afficher un taux de réussite de prévisions supérieur à 80% et une grande robustesse de la méthode. Les résultats obtenus démontrent donc, sans oublier qu’en amont les prévisions météorologiques gardent certaines limites, la faisabilité, en Bourgogne mais sans doute sous d’autres cieux, d’un système opérationnel de prévision des pics de concentration d’ozone. Ceci souligne aussi la nécessité absolue de disposer des modèles dynamiques de prévisions atmosphériques les plus fins tout comme d’excellents réseaux de mesure de la qualité de l’air, données fondamentales sans lesquelles la mise en œuvre de ce type de prévisions resterait à l’état de projet.