Dans cette thèse, nous analysons les phénomènes météorologiques régionaux révélés par les récentes missions martiennes. L'approche est numérique et observationnelle : nous avons construit un modèle météorologique méso-échelle complet pour la planète Mars, dont les applications potentielles sont nombreuses, et réalisé une cartographie à haute résolution de la pression de surface en nous appuyant sur les données du spectro-imageur OMEGA à bord du satellite Mars Express. Nous nous concentrons notamment sur les intenses circulations atmosphériques au voisinage des volcans, canyons et cratères martiens. Nous analysons également l'activité météorologique régionale des plaines martiennes, bien plus marquée que ne l'indique l'analyse synoptique. Les mouvements convectifs dans la couche limite martienne sont étudiés. Mars apparaît comme un environnement propice à d'intenses phénomènes météorologiques régionaux.

L'objet de ce travail de thèse est l'étude en soufflerie du sillage lointain des éoliennes à axe horizontal. La complexité phénoménologique du sillage des éoliennes fait que les mécanismes le régissant sont généralement traités d'un point de vue stationnaire : les principaux paramètres (déficit de vitesse, production de turbulence...) sont alors moyennés dans le temps. Néanmoins, considérer les instationnarités du sillage d'une éolienne placée en écoulement atmosphérique permet d'observer un phénomène appelé meandering, traduisant un battement aléatoire du sillage. Ce travail est construit en deux grandes parties : la vision stationnaire et la vision instationnaire du sillage d'un modèle simplifié d'éolienne, basé sur le principe du disque de Froude, placé dans une couche limite atmosphérique (CLA) modélisée en soufflerie à l'échelle 1/400. La première partie est composée de plusieurs études paramétriques sur l'influence du point de fonctionnement d'une éolienne, sur sa hauteur de mât ou encore sur le type de CLA. Une éolienne en situation de dérapage fait également l'objet d'études paramétriques. Celles-ci visent à étudier le comportement du sillage dans différentes situations et il est alors montré les limites des lois empiriques présentes dans la littérature, concernant notamment l'influence de la turbulence ambiante. La seconde partie se focalise sur la vision instationnaire, jusqu'ici rarement considérée. En utilisant l'anémométrie par fil chaud, les résultats montrent le rôle des grandes échelles de la turbulence atmosphérique sur l'apparition du meandering. L'amplitude du battement et les dimensions du sillage instantané sont ensuite appréhendées de façon quantitative par mesures PIV, montrant l'influence du point de fonctionnement de l'éolienne et de l'intensité de turbulence ambiante.

La couche-limite atmosphérique turbulente stable, particulièrement en zone de relief, n'est pas totalement comprise. Elle est, donc, mal représentée par les modèles atmosphériques. En présence de pente et d'un refroidissement du sol, l'augmentation locale de masse volumique génère un écoulement catabatique. En région de montagne, le maximum de vent est généralement enregistré à une hauteur (z_j) de 1-10 m. Le jet de paroi engendre un changement de signe du flux de qdm ainsi qu'une variabilité du flux de chaleur sensible proche du sol. Ces variabilités de flux turbulents contredisent l'applicabilité de la théorie des similitudes de Monin-Obukhov (TSMO), pourtant utilisée de manière universelle dans les modèles atmosphériques. Si la TSMO est discutable pour les cas très stables, c'est en présence de pente qu'elle n'est naturellement plus valide puisqu'elle néglige le couplage entre les équations de vitesse et de température. Il est donc impossible de représenter correctement un écoulement catabatique (z_j O(1m)) par un modèle utilisant la TSMO et avec une résolution verticale de l'ordre de la hauteur du maximum de vent. L'objectif du travail de thèse est d'apporter une contribution dans la compréhension et la modélisation de ces écoulements.Afin de compléter les observations peu nombreuses, une campagne de mesure a été réalisée sur une pente raide (20-40 deg) : la pente ouest du Grand Colon (chaîne de Belledonne, Alpes). Les analyses spectrales témoignent de la sensibilité de l'écoulement local aux perturbations externes, même faibles. Les caractéristiques turbulentes classiques sont observées à haute-fréquence alors que des comportements moins standards sont observés aux fréquences intermédiaires ou basses et expliqués par la présence de perturbations turbulentes d'énergie du même ordre de grandeur que l'injection locale. Les cospectres montrent un comportement propre aux écoulements catabatiques: recouvrement progressif selon z des corrélations croisées <0 et >0. La TSMO est mise en défaut pour l'écoulement observé et une solution alternative est utilisée pour estimer les flux en surface, permettant une bonne description de la vitesse de frottement.Le modèle 1D de surface de ISBA (Météo-France) est modifié pour répondre à la modélisation des écoulements catabatiques. Dans un premier temps, le modèle est validé sur un cas standard: en comparaison avec un modèle de Prandtl adapté. Dans un second temps, les données in-situ sont modélisées, d'abord en fournissant des profils de diffusivités effectives puis en utilisant un modèle modifié de turbulence d'ordre 1.5. Les modélisations 1D représentent correctement les champs moyens de vitesse et température mais montrent cependant des comportements trop diffusifs. Le modèle de longueur de mélange est principalement remis en cause, y compris en utilisant des paramétrisations adaptées.Des simulations LES 3D réalistes (Meso-NH, Météo-France) sont effectuées à haute résolution pour représenter le cas d'étude. Ces modélisations représentent finement les variabilités spatiales de l'écoulement catabatique. Cependant, des biais sont engendrés principalement par l'utilisation de la TSMO en condition aux limites de surface. Malgré la forte résolution spatiale, l'utilisation de la TSMO repousse à seulement z=2 m la perception des termes sources de l'écoulement catabatique par le modèle, alors que la source de l'écoulement atteint son maximum précisément en surface. Les modèles analytiques d'écoulement catabatique (de type Prandtl, qui pourraient aisément être intégrés en conditions aux limites) nécessitent de connaître "a priori" les profils de diffusivité. Ceci implique l'utilisation d'un modèle de turbulence. Le couplage du modèle 1D de surface (précédemment modifié et validé "off-line") est donc proposé pour répondre au manque de description de la physique par les CaL classiques de surface. Le travail préliminaire du couplage est présenté et des solutions sont proposées en perspective. / The stable atmospheric boundary layer, particularly in complex terrain, is not yet fully understood and it is thus still inadequately modelled. A surface cooling of a sloping terrain generates katabatic wind due to local density increase. This flow behaves as a wall-bounded turbulent jet, often simply modelled by a local balance between the buoyancy force and the turbulent friction. In mountainous regions, the wind maximum is typically observed at a height (z_j) of 1-10 m above the ground. The wall-bounded jet is responsible for a momentum-flux sign change and a heat-flux variability close to the ground. Those turbulent-flux variabilities are fully conflicting with the aplicability of the Monin-Obukhov similarity theory (MOST), which is nevertheless universally used in the atmospheric models to provide the surface boundary condition. If the MOST is already questionable for the very stable cases, it is obviously not valide over sloping surfaces because it neglectes the coupling of the wind and temperature equations, which constitutes the katabatic source. Hence, it is not possible to adequatly model a katabatic flow (z_j O(1m)) using the MOST, especially with a vertical resolution of the order of magnitude of z_j. The aim of the this PhD work is thus to improve the current understanding and modelling capacity of the katabatic winds.Since data sets of turbulent-katabatic-flow measurements are still scarce, a new field campain was carried out on a steep slope (20-40 deg): the west face of the Grand-Colon mountain (Belledonne ridge, French Alps). The experimental setup was mainly composed of a 6m mast with four sonic-anemometer levels (1, 2, 4 and 6m) to measure the turbulence on both sides of the katabatic jet. The spectral analysis shows the hight sensitivity of the local flow to external perturbations, even when these are weak. The hight-frequency subrange shows a classical behaviour (energy-injection frequency, inertial subrange), but the spectra of the intermediate and low-frequency subranges are less typical: turbulent perturbations with an energy of the order of magnitude of the local injection are present. A specific cospectra behaviour of the katabatic flows is shown: negative and positive cross-correlations overlap gradually, increasing z. The MOST fails in representing the observed flow and a surface-flux alternative estimation is succesfully used to describe the friction vellocity.The 1D surface model of ISBA (Météo-France) is modified to model katabatic flows. The model is firstly validated with a standard calibrated Prandtl model (with variable eddy difusivity). Secondly, the field data are modelled both with a prescribed effective diffusivity (from data) and using the 1.5-order turbulence scheme. The mean velocity and temperature fields are well reproduced, but it appears that the model is over-diffusive (which generates excessive fluxes), even when an adapted mixing-length is used.Realistic 3D LES simulations (Meso-NH, Météo-France) are computed with high resolution to model the field data. Spatial flow variabilities over sloping terrain are finely represented, but are biased, mainly due to the using of MOST for the surface boundary counditions. The using of MOST shifts the start of the katabatic source detection by the atmospheric model to a height of 2 m, while the katabatic source reaches its maximum at the surface. Analytical katabatic models (of the Prandtl type, which could be easily used to feed surface boundary counditions) need an "apriori" definition of the eddy and heat diffusivities. Currently, the general definition of these diffusivities is only possible by the use of turbulent models that include closures. The coupling of the previously-presented 1D surface model (validated off-line) is suggested to overcome the lack of physics description in the classic surface boundary counditions. Preliminary work on this coupling is developed and perspective solutions are proposed.

Cette thèse porte sur l'étude de la transition d'après-midi (TA) de la couche limite atmosphérique (CLA). La TA est une période complexe du cycle diurne, de par son manque de stationnarité et parce qu'elle est régie par un ensemble de forçages qui faiblissent et dont les rôles respectifs changent par rapport à ce qu'ils ont pu être durant la période convective. Il résulte une moins bonne compréhension de cette période que des régimes quasi-stationnaires convectifs, neutres ou stables. La TA pourrait conditionner l'établissement du brouillard nocturne, influencer le développement de la couche limite du lendemain et être une phase clé du cycle diurne pour la ventilation des espèces en trace vers l'atmosphère libre. Par conséquent, de meilleures description et compréhension de la TA pourraient améliorer les modèles météorologiques et de qualité de l'air. Les principaux objectifs de la thèse ont été de mieux comprendre comment la turbulence décroît pendant la TA en mettant l'accent sur l'évolution de la structure verticale de la turbulence lorsque les transferts d'énergie en surface diminuent progressivement, et sur le rôle que joue la surface à ce moment particulier du cycle diurne. Le projet international BLLAST (Boundary Layer Late Afternoon and Sunset Turbulence) a servi de cadre à cette thèse. La campagne de mesures qui s'est déroulée en été 2011 sur le Plateau de Lannemezan, au nord de la chaîne Pyrénéenne, a fourni les nombreuses observations utilisées dans cette thèse. Une approche complémentaire aux observations est l'utilisation d'un modèle atmosphérique permettant de résoudre explicitement les échelles de la turbulence (LES). / This thesis focuses on the afternoon transition (AT) of the atmospheric boundary layer (ABL). This transitional period is complex, because it is non-stationary and most of the forcings, though smaller than during the previous convective period, may come into play. Thus, this transitional period is less understood than the quasi-stationary convective, neutral or stable regimes. Yet, the AT could impact on the nocturnal fog set up, influence the boundary-layer development on the following day and play a crucial role in the transport and dispersion of pollutants and trace species towards the free troposphere. Therefore, a better understanding of the AT could improve the meteorological models. One of the main objectives of the thesis is to improve the knowledge of the decay of turbulence during the AT, when the surface energy transfers are gradually decreasing. We especially put emphasis on the evolution of the ABL turbulence vertical structure. This work is in the core of the BLLAST (Boundary Layer Late Afternoon and Sunset Turbulence) international project. A field campaign took place in summer 2011 in France, on the northern side of the Pyrenean foothills (at "Plateau de Lannemezan"), providing numerous observations used in the thesis. In addition, we also used a Large-Eddy Simulation (LES) with which the turbulent scales can be explicitly resolved. Thanks to this approach, the evolution during the AT of the scale definitions, scaling laws, turbulence characteristics and of the role of the boundary layer processes are now better understood. Thanks to the set of observations, an overview of the BLLAST ATs has been done, according to the various meteorological conditions as well as the structure and evolution of the mean and turbulence structure of the ABLs. A large variability of the AT duration was observed, depending on the surface characteristics and atmospheric conditions. The evolution of the turbulent characteristics during the AT has been studied at the surface and higher in the ABL by means of aircraft measurements and LES data. The study points out the LES ability to reproduce the turbulence evolution throughout the afternoon.

Ce travail de recherche concerne l’interaction entre une couche limite atmosphérique et une canopée (représentant un couvert forestier). J’ai étudié le problème complexe de production et d’évolution de grosses structures turbulentes au dessus de couverts homogènes et hétérogènes, moyennement denses. J’ai abordé ce problème en mettant en œuvre les outils de la simulation numérique des grosses structures (LES) et du calcul haute performance (HPC). Les résultats numériques obtenus, reproduisent correctement les principales caractéristiques de cet écoulement, telles qu’elles sont rapportées dans la littérature : la formation d’une première génération de structures cohérentes en rouleaux, orientées transversalement par rapport à la direction de l’écoulement principal, puis la réorganisation et la déformation de ces structures qui évoluent vers une forme en fer à cheval. Les résultats obtenus au dessus d’un couvert discontinu (représentant une clairière ou une coupure de combustible dans une forêt), ont été comparés avec des données expérimentales collectées dans une soufflerie. Ceux-ci confirment l’existence d’un niveau élevé de turbulence au sein même du couvert végétal à une distance égale à 8 fois la hauteur de canopée. Cette zone, (appelée « Enhance Gust Zone » dans la littérature), est par ailleurs caractérisée par l’existence d’un pic local du facteur de dissymétrie (« skewness factor »).Le transport d’un contaminant passif émis par le feuillage a été également étudié, dans deux configurations, en supposant que la concentration à la surface du feuillage pouvait être considérée comme constante (source infinie) ou variable (source finie) en fonction du temps. Les résultats montrent un impact significatif de cette hypothèse sur la dynamique et le niveau des concentrations relevées dans l’atmosphère. / The work is dedicated to the investigation of the interaction between an Atmospheric Boundary Layer and a canopy (representing a forest cover). We have focused our attention to the complex problem of the generation and transformation of turbulent vortices over homogeneous, heterogeneous and sparse canopy. This problem has been studied using Large Eddy Simulation (LES) approach and High Performance Computing (HPC) technique.The numerical results reproduced correctly all the main characteristics of this flow, as reported in the literature: the formation of a first generation of coherent structures aligned transversally with the wind flow direction, the reorganization and the deformation of these vortex tubes into horse-shoe structures. The results obtained with the introduction of a discontinuity in the canopy (reproducing a clearing or a fuel break in a forest) are compared with the experimental data collected in a wind tunnel. In this case, the results confirmed the existence of a strong turbulence activity inside the canopy at a distance equal to 8 times the height of the canopy, referenced in the literature as the Enhance Gust Zone (EGZ) characterized by a local peak of the skewness factor. Then, the process of passive scalar transport from a forest canopy into a clear atmosphere is studied for two cases, i.e., when the concentration held by the forest canopy is either constant or variable. While this difference has little influence on the concentration patterns, results show that it has an important influence on the concentration magnitude as well as on the dynamics of the total concentration in the atmosphere.

L'ozone troposphérique, polluant secondaire affectant la santé des êtres vivants, concerne particulièrement le département des Alpes-Maritimes au cours de la saison photochimique. Espace littoral montagneux, ce territoire est largement dominé durant l'été notamment par des conditions anticycloniques permettant aux brises thermiques de s'exprimer. Ce régime de vent est au cœur de la problématique de la pollution par l'ozone puisqu'il entraîne fréquemment une accumulation des polluants primaires et secondaires au fil des jours au sein de la même masse d'air. Bien qu'étant un département faiblement industrialisé, les Alpes-Maritimes sont victimes d'un fort ensoleillement qui permet aux polluants primaires, émis majoritairement par le trafic routier près du littoral et par la végétation dans l'arrière-pays, de produire de l'ozone.<br />A travers les données de pollution issues du réseau de surveillance de la qualité de l'air d'AtmoPACA ainsi qu'à partir de très nombreuses mesures de terrain, l'objectif est de mieux appréhender les relations entre la variabilité spatiale et temporelle de l'ozone et celle des conditions météorologiques à différentes échelles. Après avoir détaillé l'historique des mesures d'ozone et de dioxyde d'azote disponibles dans le département, une première approche à macro-échelle est menée entre les réanalyses du NCEP et les niveaux de pollution par l'ozone dans neuf stations de mesures des Alpes-Maritimes. Ce premier niveau d'analyse permet de définir les configurations météorologiques générales caractérisant un épisode de pollution par l'ozone. La présence d'une dorsale anticyclonique sur l'Europe de l'ouest, associée à de faibles vitesses de vent, de faibles taux d'humidité relative en surface et d'une faible vorticité relative, provoque une dégradation de la qualité de l'air dans le département.<br /> Un second niveau d'analyse est alors abordé : il s'agit de préciser à méso-échelle et à micro-échelle les conditions météorologiques propices à de fortes concentrations d'ozone. Pour cela des campagnes de mesures itinérantes d'ozone sont effectuées dans l'ensemble du département ; une importante base de données est constituée sur la commune de Nice et dans ses alentours. Le recours au modèle météorologique RAMS permet alors de mieux comprendre la variabilité spatiale et temporelle de l'ozone induite par les conditions climatiques locales. Les variables météorologiques les plus corrélées aux concentrations d'ozone sont la vitesse du vent en surface, l'énergie cinétique turbulente, la hauteur de la couche limite atmosphérique et l'humidité relative. Tout indique que moins le volume d'air dans lequel les polluants primaires sont émis est important, et moins l'intensité du brassage de l'air est forte, alors plus les concentrations d'ozone sont élevées. Il semblerait également que ces conditions locales du temps aient plus d'influence sur les niveaux de pollution par l'ozone que la configuration synoptique à macro-échelle. Bien qu'étant nécessaire, la présence d'un anticyclone sur l'Europe de l'ouest n'est pas suffisante pour expliquer le comportement de l'ozone localement dans les Alpes-Maritimes.<br />Parallèlement à ces différentes approches, les origines spatiales de la pollution photochimique affectant cet espace littoral montagneux sont recherchées. L'advection d'importantes quantités d'ozone par la brise de mer en journée indique que ce polluant s'accumule au dessus de la mer ; excepté ce phénomène d'accumulation induit par des aller-retour de la masse d'air guidée par l'alternance entre la brise de mer et la brise de terre, l'origine des fortes concentrations d'ozone au dessus de la mer n'est pas clairement comprise. Il semblerait qu'un schéma de recirculation des masses d'air en trois dimensions permette la création d'un empilement de couches stratifiées sur la mer.

La pollution atmosphérique à l'échelle régionale est le bilan des interactions entre des processus très différents: les émissions , la chimie, le transport, le mélange et le dépôt des espèces gazeuses. La prévision de la qualité de l'air nécessite donc l'utilisation de modèles, qui prennent en compte les émissions par le biais de cadastres. Les concentrations de polluants simulées dépendent fortement des émissions utilisées. Or, les cadastres utilisés pour les représenter comportent de grandes incertitudes. Comme il est difficile actuellement d'affiner leurs méthodes de construction, il reste la possibilité d'ajouter de l'information aux cadastres existant. L'optimisation des émissions utilise l'information contenue dans les mesures pour obtenir le cadastre qui minimise la différence entre les concentrations simulées et mesurées.<br /><br />Une méthode d'inversion des émissions anthropiques à l'échelle régionale, utilisant les mesures de routine et basée sur le modèle CHIMERE et son adjoint, a été élaborée et validée. Une technique de krigeage permet d'utiliser de façon optimale les informations disponibles dans l'espace des concentrations. L'enchaînement de cycles krigeage-optimisation améliore la qualité des résultats. Une technique d'agrégation spatiale dynamique est utilisée pour réduire la dimension du problème.<br /><br />Les émissions de NOx du cadastre élaboré par AIRPARIF pour l'Ile-de-France ont été inversées pendant les étés 1998 et 1999, les épisodes de la campagne ESQUIF étant étudiés en détail. L'optimisation corrige des écarts importants entre concentrations mesurées et simulées. Cependant, d'une façon générale, le niveau de fiabilité des résultats diminue avec la densité du réseau de mesure. Les résultats présentant le plus haut niveau de confiance concernent donc les flux d'émission les plus intenses d'Ile-de-France. Les corrections apportées à la masse moyenne émise dans l'ensemble du domaine et aux profils temporels correspondant sont en accord avec l'estimation de l'incertitude sur le cadastre pour les NOx (15%) obtenue lors de la campagne ESQUIF.

La prévision des phénomènes de pollution atmosphérique est l'objectif central de la campagne ESCOMPTE-2001, qui s'est déroulée durant l'été 2001 dans la région très hétérogène de Marseille/Fos/Berre. Cet objectif nécessite une bonne compréhension et prise en compte, par les modèles numériques de physico-chimie, des processus physiques intervenant dans la Couche Limite Atmosphérique (CLA), au sein de laquelle s'effectuent le transport et la diffusion des polluants émis en surface.<br />Dans le cadre de la campagne ESCOMPTE-2001, ce travail de thèse est consacré à l'étude de la basse troposphère en période de brise de mer, situation météorologique défavorable à la qualité de l'air des zones côtières. Il a permis notamment de mettre en évidence une oscillation de l'intensité de la brise de mer et des compétitions de brises de mer locales et régionales, qui modifient le temps d'advection de l'air marin sur la surface terrestre et ont donc des répercussions importantes sur le développement de la CLA et sa concentration en polluants. Ces travaux s'appuient principalement sur le réseau de quatre radars profileurs de vent UHF mis en place sur la zone côtière de Marseille/Fos/Berre, permettant une description en trois dimensions et en continu dans le temps des écoulements de brises de mer et de la CLA.<br />Pour les besoins de cette étude phénoménologique, des développements méthodologiques sont apportés sur la mesure des propriétés turbulentes de la CLA par les radars UHF (termes du bilan de l'énergie cinétique turbulente) et sur leur utilisation en réseau pour l'étude de la trajectographie de panaches de polluants.

Les Modèles de Chimie-Transport (CTM) simulent aujourd'hui de façon assez réaliste les concentrations des espèces responsables de la pollution photochimique estivale. Cependant, la complexité des phénomènes mis en jeu et la variabilité spatiale et temporelle des émissions de polluants sont telles que, quel que soit le modèle utilisé, il est impossible de reproduire parfaitement ces concentrations. Lorsqu'un paramètre interne ou d'entrée du modèle est mal décrit une erreur importante peut etre observée certains jours. Le modèle CHIMERE est un CTM couvrant la majeur partie de l'Europe avec une résolution d'une cinquantaine de kilomètres. Il offre la possibilité de zoomer et de simuler de façon plus détaillée les concentrations des polluants sur des régions clefs, telles que l'Ile-de-France et la région de Berre et de Marseille.<br /><br />Les simulations de ce modèle ont été compaés avec des observations de surface et des données aéroportées de la campagne d'Etude et Simulation de la QUalité de l'Air en Ile-de-France (ESQUIF) à l'échelle européenne et régionale. Cette comparaison a permis ainsi de quantifier l'erreur globale commise sur les concentrations d'ozone et de son précurseur, le dioxyde d'azote. Différentes méthodes (Interpolation Statistique, Krigeage) ont été testées et adaptées au cas de la pollution dans le but de corriger cette erreur. Elles ont été comparées et validées de manière objective. Il a été montré que la combinaison des observations de surface et des simulations du modèle permet d'obtenir efficacement des cartes tridimensionnelles (i.e. des analyses) de concentration des deux polluants, plus réalistes que les simulations brutes.<br /><br />Les expériences menées depuis trois ans dans le cadre du projet PIONEER (Prévisibilité et Incertitude de l'Ozone à l'échelle Européenne et Régionale) montrent que l'erreur de prévision peut éventuellement se propager d'une région vers une autre. Les analyses des concentrations d'ozone, produites à l'échelle européenne ont été également utilisées pour réinitialiser le modèle de prévision. L'objectif était alors de savoir s'il est possible d'améliorer aussi les prévisions à court termes en utilisant de meilleurs états initiaux que les prévisions de la veille.

Avec 70% de la surface de la planète recouverte par les mers et les océans, la présence de particules salées, produites majoritairement par le déferlement de vagues, représente une composante majeure dans le cycle géochimique de l'atmosphère et dans le bilan radiatif terrestre. Nos connaissances doivent progresser vers une description et une modélisation des flux actuels plus précises afin d'être capable d'anticiper les changements éventuels de la composition de l'atmosphère. Ce travail de thèse s'inscrit dans le cadre général d'une meilleure caractérisation des particules injectées dans l'atmosphère par le déferlement des vagues. Un des objectifs est d'etudier la validité des modèles paramétriques et numériques pour prédire la variation spatio-temporelle des particules d'aérosols à l'échelle locale et régionale. Une première partie du travail est consacrée à l'étude de la variation spatiale des concentrations d'aérosols dans une zone côtière méditerranéenne. Pour ce faire, un couplage entre un code météorologique méso-échelle (RAMS) et un modèle déterministe de distributions d'aérosols développé par le LSEET (Medex)a été réalisé. Les prédictions du couplage sont confrontées aux données enregistrées lors d'une campagne expérimentale menée en mai 2007 à la station de mesures du laboratoire sur l'île de Porquerolles. Une deuxième partie du travail s'est concentrée sur l'amélioration du modèle paramétrique Medex. A partir, de données enregistrées en mai 2007, les prédictions de Medex ont été affinées, en particulier par correction de l'influence saisonnière. Dans une troisième partie, on s'est intéressé aux paramètres qui influencent les variations de concentrations en aérosols lors d'un épisode de Mistral. Il est montré que des paramètres autres que la vitesse du vent interviennent, tel que le fetch et la hauteur de la couche limite marine. Une dernière partie est consacrée au développement d'un modèle numérique de transport des aérosols. Les premiers résultats permettent de montrer l'influence de la stabilité atmosphérique sur les profils horizontaux et verticaux d'aérosols.