Les régimes dynamiques dans les basses couches de l'atmosphère à l'île de la Réunion sont conditionnés par l'action du relief et du rayonnement sur l'écoulement synoptique. L'île est située en permanence dans le flux des alizés de sud-est et son relief élevé, culminant à 3000 m dans le centre de l'île et 2600 m au sud, constitue un obstacle important pour l'écoulement moyen. Le relief, le chauffage différentiel des pentes et le contraste thermique avec l'océan influent sur les échanges locaux entre la couche limite marine et la troposphère libre.L'analyse des phénomènes complexes de l'écoulement atmosphérique au niveau de l'île a pour but initial de caractériser les transferts de polluants émis localement. Cette étude a deux applications : * La première rentre dans le cadre préparatoire à la mise en place, à 2200 m d'altitude sur le Piton Maïdo, de l'observatoire atmosphérique de l'OPAR, à l'horizon 2011. L'objectif est de comprendre les circulations locales induites par le relief et les transports associés afin de discriminer d'éventuelles pollutions par les sources locales sur le signal qui sera mesuré in situ au sommet du Maïdo.* La seconde vise à étudier la diffusion des panaches volcaniques du Piton de la Fournaise. L'éruption majeure d'avril 2007 du volcan réunionnais a montré que des panaches pouvaient générer d'importantes pollutions dans diverses parties de l'île allant jusqu'à poser de réels problèmes environnementaux et de santé publique.Ces applications s'appuient sur une étude par modélisation numérique à haute résolution des écoulements atmosphériques dans les basses couches au niveau de l'île, au moyen du modèle météorologique de recherche MésoNH, permettant de conduire un ensemble de simulations sur cas idéalisés puis sur cas réels avec diffusion de traceurs passifs.

La modélisation par simulation des grandes échelles (Large-Eddy Simulation - LES) des processus physiques régissant la couche limite atmosphérique (CLA) demeure complexe de part la difficulté des modèles à capter l'évolution de la turbulence entre différentes conditions de stratification. De ce fait, l'étude LES du cycle diurne complet de la CLA comprenant des situations convectives la journée et des conditions stables la nuit est très peu documenté. La simulation de la couche limite stable où la turbulence est faible, intermittente et qui est caractérisée par des structures turbulentes de petite taille est tout particulièrement compliquée. En conséquence, la capacité de la LES à bien reproduire les conditions météorologiques de la CLA, notamment en situation stable, est étudiée à l'aide du code de mécanique des fluides développé par EDF R&D, Code_Saturne. Dans une première étude, le modèle LES est validé sur un cas de couche limite convective quasi stationnaire sur terrain homogène. L'influence des modèles sous-maille de Smagorinsky, Germano-Lilly, Wong-Lilly et WALE (Wall-Adapting Local Eddy-viscosity) ainsi que la sensibilité aux méthodes de paramétrisation sur les champs moyens, les flux et les variances est discutées. Dans une seconde étude le cycle diurne complet de la CLA pendant la campagne de mesure Wangara est modélisé. L'écart aux mesures étant faible le jour, ce travail se concentre sur les difficultés rencontrées la nuit à bien modéliser la couche limite stable. L'impact de différents modèles sous-maille ainsi que la sensibilité au coefficient de Smagorinsky ont été analysés. Par l'intermédiaire d'un couplage radiatif réalisé en LES, les répercussions du rayonnement infrarouge et solaire sur le jet de basse couche nocturne et le gradient thermique près de la surface sont exposées. De plus l'adaptation de la résolution du domaine à l'intensité de la turbulence et la forte stabilité atmosphérique durant l'expérience Wangara sont commentées. Enfin un examen des oscillations numériques inhérentes à Code_Saturne est réalisé afin d'en limiter les effets

Le transport de la neige par le vent est une composante importante de l'interaction entre l'atmosphère et la cryosphère. En zone de montagne, il influence la distribution temporelle et spatiale de la couverture neigeuse au cours de l'hiver et a en premier lieu des conséquences sur le danger d'avalanche. La modélisation numérique de ce phénomène permet d'étudier les interactions complexes entre le manteau neigeux et le vent et d'en estimer les conséquences de manière distribuée. Dans ce contexte, cette thèse décrit le développement et l'évaluation d'un modèle couplé atmosphère/manteau neigeux dédié à l'étude du transport de la neige par le vent en zone de montagne reposant sur le modèle atmosphérique Meso-NH et le modèle détaillé de manteau neigeux Crocus. Le transport de la neige par le vent a été étudié sur le site expérimental du Col du Lac Blanc (massif des Grandes Rousses, France). Une base de données d'épisodes de transport couvrant dix hivers a tout d'abord été utilisée pour déterminer les caractéristiques principales de ces épisodes. Des simulations avec le modèle Crocus (non couplé à Meso-NH) ont ensuite montré qu'il était nécessaire de tenir compte des transformations mécaniques des grains de neige induites par le vent afin de simuler une évolution réaliste de la vitesse seuil de transport. Le site expérimental a également été le siège de deux campagnes de mesures en 2011 et 2012 visant à collecter de données de validation pour le modèle. Elles renseignent sur les conditions météorologiques près de la surface, sur les quantités de neige transportées et sur la localisation des zones d'érosion et de dépôt de la neige grâce à l'utilisation d'un laser terrestre. Le modèle de transport de neige par le vent Meso-NH/Crocus a été développé. Il intègre le transport de la neige en saltation et en suspension turbulente ainsi que la sublimation des particules de neige transportée. Un schéma à deux moments permet de simuler l'évolution spatiale et temporelle de la distribution en taille des particules. L'utilisation d'un schéma de couche limite de surface à l'interface entre Meso-NH et Crocus s'est révélé nécessaire pour représenter les forts gradients de concentration en particules de neige observés près de la surface. Meso-NH/Crocus est le premier modèle couplé atmosphère/manteau neigeux capable de simuler de manière interactive le transport de la neige par le vent en zone alpine. Meso-NH/Crocus a été évalué en relief réel grâce aux données collectées lors de la première campagne de mesure en 2011. La simulation d'un épisode de transport sans chute de neige simultanée montre que le modèle reproduit de manière satisfaisante les principales structures d'un écoulement en relief complexe ainsi que les profils verticaux de vitesse de vent et de flux de particules de neige en suspension près de la surface. En revanche, la résolution horizontale de 50 m est insuffisante pour reproduire avec précision la localisation des zones d'érosion et de dépôt autour du Col du Lac Blanc. La prise en compte de la sublimation réduit la quantité de neige déposée de l'ordre de 5%.Les techniques de descente d'échelle dynamique (grid nesting) ont ensuite été utilisées pour simuler un second épisode de transport avec chute de neige. L'augmentation de la résolution horizontale intensifie les contrastes de vitesse de vent entre versants au vent et sous le vent. En revanche, elle modifie peu les quantités et les structures spatiales des précipitations solides autour du Col du Lac Blanc. Lorsqu'il est activé, le transport devient la principale source d'hétérogénéités des accumulations neigeuses / Blowing and drifting snow are crucial components of the interaction between the cryosphere and the atmosphere. In mountainous areas, it affects the temporal and spatial distribution of snow depth throughout the winter season and influences avalanche formation. Numerical modeling offers a solution for studying the complex interaction between the snowpack and the wind field and to assess the related processes in a spatially distributed way. In this context, this PhD describes the development and the validation of a coupled snow/atmosphere model which is dedicated to the study of blowing and drifting snow in alpine terrain. The coupled model consists in the atmospheric model Meso-NH and the detailed snowpack model Crocus. Blowing and drifting snow have been monitored at the Col du Lac Blanc (Grandes Rousses range, French Alps) experimental site. A database consisting of blowing snow events observed over 10 years allowed us to identify the main features of these events. Numerical simulations using Crocus illustrated the necessity of taking the wind-dependence of snow grain characteristics into account in order to simulate satisfactorily the occurrence of blowing snow events. We also carried out two measurement campaigns at our experimental site in 2011 and 2012 in order to collect validation data for the model. This includes measurements of vertical profiles of wind speed and snow particle fluxes near the surface and the mapping of areas of erosion and deposition using terrestrial laser scanning. The coupled Meso-NH/Crocus model has been developed in order to account for blowing and drifting snow. It simulates snow transport in saltation and in turbulent suspension and includes the sublimation of suspended snow particles. In the atmosphere, a double-moment scheme allows the model to simulate the spatial and temporal evolution of the snow particle size distribution. The implementation of a surface boundary layer scheme at the interface between Meso-NH and Crocus turned out to be necessary to reproduce the strong vertical gradient of snow particle concentration near the surface. Meso-NH/Crocus is the first coupled snow-atmosphere model that can simulate snow transport in alpine terrain in an interactive way.Meso-NH/Crocus has been evaluated against data collected near Col du Lac Blanc during the first measurement campaign in 2011. The simulation of a blowing snow event without concurrent snowfall showed that the model captures the main structures of atmospheric flow in complex terrain, the vertical profile of wind speed and the snow particle fluxes. However, the horizontal resolution of 50 m is found to be insufficient to simulate the location of areas of snow erosion and deposition observed around Col du Lac Blanc. Blowing snow sublimation leads to a reduction in snow deposition of approximately 5%.We used downscaling techniques (grid nesting) to simulate a second blowing event with concurrent snowfall. The increase in horizontal resolution enhanced the contrast of wind speed between windward and leeward slopes. However, it only slightly affects the amount and the spatial pattern of snow precipitation around Col du Lac Blanc. When activated, blowing and drifting snow are the main sources of spatial variability of snow accumulation

Afin de mieux comprendre les mécanismes liés à l’érosion du sol sous l’effet du vent, le transport de particules solides dans un écoulement de couche limite turbulente à l’échelle d’une soufflerie est étudié à l’aide de simulations numériques. La présence d’une ou plusieurs collines Gaussiennes au sol permet d’étudier les effets de la topographie sur le transport, le dépôt et la réémission de particules solides. L’écoulement du fluide porteur est résolu par la Simulation des Grandes Échelles (SGE). Des modèles de paroi pour la vitesse du fluide sont implémentés afin de mieux représenter l’écoulement proche d’une colline. Le mouvement des particules est pris en compte par un suivi Lagrangien. Des modèles d’envol et de rebond sont développés et utilisés pour prendre en compte l’émission et l’impact au sol des particules. Dans la première partie, l’écoulement au-dessus de collines transversales est simulé et validé par des comparaisons avec différentes expériences. Selon Oke [1988], l’écoulement dans la canopée urbaine peut être schématiquement caractérisé par différents régimes en fonction du positionnement relatif des obstacles. Ce concept est appliqué au cas des dunes, assimilées à des collines dans notre étude. L’accent est mis sur la zone de recirculation (ZR) formée derrière ces collines. Les variations de la ZR sont examinées en fonction de différents paramètres dont la configuration des collines et le nombre de Reynolds. De plus, une étude portant sur la sous couche rugueuse est effectuée de façon à déterminer l’effet de la rugosité due à la couche de particules solides au sol. La seconde partie du travail porte sur la simulation des particules au-dessus des collines. L’objectif est l’amélioration des modélisations concernant l’envol, le rebond et le couplage entre le fluide et les particules. Un premier travail de validation est réalisé en utilisant le modèle complet de transport des particules solides. En particulier, l’évolution du flux d’émission des particules, estimé par le modèle d’envol, en fonction du nombre de Shields, donne des résultats comparables aux modèles classiques de saltation et aux expériences de la littérature. Au-dessus des collines, le transport des particules solides est étudié par des profils de concentration et de vitesse moyenne. Pour analyser les résultats, deux cartographies sont réalisées. La première donne l’intensité des événements locaux et instantanés qui seraient à l’origine de l’évacuation des particules piégées au sein de la ZR. La seconde montre la distribution des particules déposées au sol. Ces résultats permettent d’identifier des zones sujettes à l’érosion et à l’accumulation autour des collines. Enfin, les flux des particules piégées et déposées à l’intérieur de la ZR sont quantifiés et comparés aux flux des particules émises en amont. Ces flux, bien que faibles par rapport au flux entrant, contribueraient aux migrations des dunes et à l’avancée des déserts. / The transport of solid particles inside a laboratory-scale turbulent boundary-layer is studied by numerical simulations, to obtain a better understanding of the mechanisms associated with wind erosion of soil. The presence of one or several Gaussian hills allows a study of the topographic effects on the transport, deposition and re-emission of solid particles. The carrier fluid motion is resolved in a Large Eddy Simulation (LES). Wall models are implemented to better account for the effects of turbulent flow near the terrain. Particle trajectories are calculated using a Lagrangian tracking. Take-off and rebound models are developed in order to take into account particle emissions and impacts at the wall. In the first part, the flow over transversal Gaussian hills is simulated and validated by comparison with different experiments. According to Oke [1988], the flow inside an urban canopy can be schematically characterised into different flow regimes depending on the relative localisation of the obstacles at the ground. This concept is applied to the case of sand dunes, assimilated to 2D hills in this study. The focus is on the recirculation zone (RZ) on the lee side, which has the characteristic of increasing the residence time and the interaction fluid/particle in general, particle trapping and deposition in particular. The variations of RZ with different hill geometries and Reynolds numbers are examined. A study on the roughness sublayer is conducted in order to determine the roughness effects due to the layer of solid particles on the wall. The second part of the work is devoted to the simulation of solid particle transport over the Gaussian hills. The objective is to improve the modelling of particle take-off, rebound and the two-way coupling between the fluid and the particle. A first work of validation is conducted by using the complete model of solid particle transport developed in this thesis. In particular, the evolution of particle emission flux predicted by the take-off model is in accordance with classical saltation models and experiments from the literature. Over the Gaussian hills, analysis of particle transport is conducted using concentration and mean velocity fields. Two mappings are realised. The first indicates the intensity of the local and instantaneous flow structures that arguably regulate the re-entrainment of particles trapped inside the RZ. The second shows the accumulation of particles on the wall. These results highlight zones prone to wind erosion and particle deposition around the hills. Last but not least, the fluxes of particle trapping and deposition inside the RZ are quantified and compared to the incoming flux from upstream. These fluxes, albeit relatively weak in comparison to the incoming one, contribute potentially to dune migrations and desertification.

Cette thèse s'inscrit dans un projet de modélisation numérique de la dispersion atmosphérique de polluants à travers le code de mécanique des fluides numérique Code_Saturne. L'objectif est de pouvoir simuler la dispersion atmosphérique de polluants en environnement complexe, c'est-à-dire autour de centrales, sites industriels ou en milieu urbain. Dans ce contexte, nous nous concentrons sur la modélisation de la dispersion des polluants à micro-échelle, c'est-à-dire pour des distances de l'ordre de quelques mètres à quelques kilomètres et correspondant à des échelles de temps de l'ordre de quelques dizaines de secondes à quelques dizaines de minutes : on parle de modélisation en champ proche. L’approche suivie dans ces travaux de recherche suit une formulation hybride eulérienne/lagrangienne, où les champs dynamiques moyens relatifs au fluide porteur (pression, vitesse, température, turbulence) sont calculés via une approche eulérienne et sont ensuite fournis au solveur lagrangien. Ce type de formulation est couramment utilisé dans la littérature atmosphérique pour son efficacité numérique. Le modèle lagrangien stochastique considéré dans nos travaux est le Simplified Langevin Model (SLM), développé par Pope (1985,2000). Ce modèle appartient aux méthodes communément appelées méthodes PDF (Probability Density Function), et, à notre connaissance, n'a pas été exploité auparavant dans le contexte de la dispersion atmosphérique. Premièrement, nous montrons que le SLM respecte le critère dit de mélange homogène (Thomson, 1987). Ce critère, essentiel pour juger de la bonne qualité d'un modèle lagrangien stochastique, correspond au fait que si des particules sont initialement uniformément réparties dans un fluide incompressible, alors elles doivent le rester. Nous vérifions le bon respect du critère de mélange homogène pour trois cas de turbulence inhomogène représentatifs d'une large gamme d'applications pratiques : une couche de mélange, un canal plan infini, ainsi qu'un cas de type atmosphérique mettant en jeu un obstacle au sein d'une couche limite neutre. Nous montrons que le bon respect du critère de mélange homogène réside simplement en la bonne introduction du terme de gradient de pression en tant que terme de dérive moyen dans le modèle de Langevin (Pope, 1987; Minier et al., 2014; Bahlali et al., 2018c). Nous discutons parallèlement de l'importance de la consistance entre champs eulériens et lagrangiens dans le cadre de telles formulations hybrides eulériennes/lagrangiennes. Ensuite, nous validons le modèle dans le cas d'un rejet de polluant ponctuel et continu, en conditions de vent uniforme et turbulence homogène. Dans ces conditions, nous disposons en effet d'une solution analytique nous permettant une vérification précise. Nous observons que dans ce cas, le modèle lagrangien discrimine bien les deux différents régimes de diffusion de champ proche et champ lointain, ce qui n'est pas le cas d'un modèle eulérien à viscosité turbulente (Bahlali et al., 2018b).Enfin, nous travaillons sur la validation du modèle sur plusieurs campagnes expérimentales en atmosphère réelle, en tenant compte de la stratification thermique de l'atmosphère et de la présence de bâtiments. Le premier programme expérimental considéré dans nos travaux concerne le site du SIRTA (Site Instrumental de Recherche par Télédétection Atmosphérique), dans la banlieue sud de Paris, et met en jeu une stratification stable de la couche limite atmosphérique. La seconde campagne étudiée est l'expérience MUST (Mock Urban Setting Test). Réalisée aux Etats-Unis, dans le désert de l'Utah, cette expérience a pour but de représenter une ville idéalisée, au travers d'un ensemble de lignées de conteneurs. Deux rejets ont été simulés et analysés, respectivement en conditions d'atmosphère neutre et stable (Bahlali et al., 2018a) / This Ph.D. thesis is part of a project that aims at modeling pollutant atmospheric dispersion with the Computational Fluid Dynamics code Code_Saturne. The objective is to simulate atmospheric dispersion of pollutants in a complex environment, that is to say around power plants, industrial sites or in urban areas. In this context, the focus is on modeling the dispersion at micro-scale, that is for distances of the order of a few meters to a few kilometers and corresponding to time scales of the order of a few tens of seconds to a few tens of minutes: this is also called the near field area. The approach followed in this thesis follows a hybrid Eulerian/Lagrangian formulation, where the mean dynamical fields relative to the carrier fluid (pressure, velocity, temperature, turbulence) are calculated through an Eulerian approach and are then provided to the Lagrangian solver. This type of formulation is commonly used in the atmospheric literature for its numerical efficiency. The Lagrangian stochastic model considered in our work is the Simplified Langevin Model (SLM), developed by Pope (1985,2000). This model belongs to the methods commonly referred to as PDF (Probability Density Function) methods, and, to our knowledge, has not been used before in the context of atmospheric dispersion. First, we show that the SLM meets the so-called well-mixed criterion (Thomson, 1987). This criterion, essential for any Lagrangian stochastic model to be regarded as acceptable, corresponds to the fact that if particles are initially uniformly distributed in an incompressible fluid, then they must remain so. We check the good respect of the well-mixed criterion for three cases of inhomogeneous turbulence representative of a wide range of practical applications: a mixing layer, an infinite plane channel, and an atmospheric-like case involving an obstacle within a neutral boundary layer. We show that the good respect of the well-mixed criterion lies simply in the good introduction of the pressure gradient term as the mean drift term in the Langevin model (Pope, 1987; Minier et al., 2014; Bahlali et al., 2018c). Also, we discuss the importance of consistency between Eulerian and Lagrangian fields in the framework of such Eulerian/Lagrangian hybrid formulations. Then, we validate the model in the case of continuous point source pollutant dispersion, under uniform wind and homogeneous turbulence. In these conditions, there is an analytical solution allowing a precise verification. We observe that in this case, the Lagrangian model discriminates well the two different near- and far-field diffusion regimes, which is not the case for an Eulerian model based on the eddy-viscosity hypothesis (Bahlali et al., 2018b).Finally, we work on the validation of the model on several experimental campaigns in real atmosphere, taking into account atmospheric thermal stratification and the presence of buildings. The first experimental program considered in our work has been conducted on the `SIRTA' site (Site Instrumental de Recherche par Télédétection Atmosphérique), in the southern suburb of Paris, and involves a stably stratified surface layer. The second campaign studied is the MUST (Mock Urban Setting Test) experiment. Conducted in the United States, in Utah's desert, this experiment aims at representing an idealized city, through several ranges of containers. Two cases are simulated and analyzed, respectively corresponding to neutral and stable atmospheric stratifications (Bahlali et al., 2018a)

Le transport de neige par le vent occasionne des nuisances qui entravent certaines activités du secteur des transports et du tourisme. Il déplace les dépôts de neige et modifie la structure du manteau neigeux. A ce titre, il participe à l'augmentation du risque d'avalanche. Ces différents aspects justifient l'activité de recherche menée depuis plusieurs années au Cemagref de Grenoble et dans laquelle cette étude trouve naturellement sa place. Ainsi les travaux précédents ont contribué à la modélisation physique et numérique de prévision du transport de neige par le vent. Cela dit, le transport éolien de particules met en jeu différents phénomènes complexes que les connaissances actuelles ne suffisent pas à décrire avec précision. Dans ce contexte, des travaux expérimentaux ont été conduits à la fois en soufflerie et sur le terrain. Ces derniers ont permis des premières investigations sur le comportement des fluctuations de concentration en régime permanent, en particulier par le biais de traitements d'images et d'analyses spectrales. Cette thèse a également porté sur la variabilité des profils de concentration en conditions naturelles, c'est-à-dire dans un écoulement turbulent instationnaire. Tous les résultats in-situ ont nécessité l'usage d'instruments adaptés dont certains ont du être réétalonnés de manière spécifique.

Le relief des vallées modifie la structure et l'évolution de la couche limite atmosphérique (CLA) par comparaison avec le sol plat et ainsi affecte le transport et le mélange de polluants. L'objet de cette étude est alors de comprendre, au moyen de l'observation et de la simulation numérique, l'évolution de la structure de l'atmosphère (en particulier la CLA) dans les vallées. A cette fin, des données issues du programme POVA, relatif à la qualité de l'air dans les vallées alpines, sont analysées. Un modèle de type simulation des grandes échelles est également adapté (et évalué au moyen des données expérimentales). Les processus d'échanges verticaux couplant la CLA convective à l'atmosphère libre (AL) sont ensuite détaillés sur sol plat. Les résultats de cette étude permettent de caractériser et de quantifier les spécificités des écoulements en vallée et les mécanismes d'entraînement à l'interface entre la CLA et l'AL sur sol plat.

L'objectif de cette thèse est d'étudier la capacité des modèles météorologiques à prévoir des épisodes de cisaillements de vent dans les basses couches de l'atmosphère sur une zone limitée à un aéroport et d'examiner l'apport pour la modélisation d'observations locales à haute fréquence. Nous avons choisi l'aéroport international de Nice, régulièrement soumis à des variations rapides de la direction et de l'intensité du vent selon l'horizontale dans la CLA, appelées aussi renverses. Un pro leur de vent et trois anémomètres sont installés sur les pistes de l'aéroport. Au début de l'année 2009, une campagne de mesures incluant un lidar vent à balayage et un anémomètre sonique s'est déroulée sur l'aéroport fournissant des observations complémentaires. L'ensemble des mesures à haute fréquence temporelle et des simulations numériques obtenues avec le modèle de recherche Méso-NH à 2.5 km de résolution, a fourni une vision de l'enchaînement complexe des écoulements conduisant à des cisaillements de vent d'origine différente. Cette complémentarité a aussi permis d'estimer la capacité du modèle numérique à reproduire les cisaillements de vent. Pour les trois situations étudiées, il reproduit la structure horizontale et verticale de l'écoulement malgré des erreurs de placement spatio-temporel. Bien que les écoulements locaux participent à la mise en place des conditions nécessaires au cisaillement de vent, c'est l'écoulement de méso-échelle (ondes piégées ou talweg d'altitude) qui va déterminer la position du phénomène. Nous avons réalisé des comparaisons avec le modèle opérationnel de Météo-France AROME ainsi que des tests de sensibilité pour étudier l'in uence des conditions de couplage et de la résolution. Nous avons, en particulier, augmenté la résolution horizontale de 2.5 km à 500 m sur un domaine centré sur l'aéroport de Nice sur les situations étudiées. Une résolution de 500 m permet d'améliorer la représentation d'écoulements locaux et de variations locales du vent mais n'améliore pas la position des cisaillements de vent par rapport à une échelle plus grossière. L'extension horizontale limitée du domaine à haute résolution augmente la sensibilité aux conditions aux limites de grande échelle. Pour améliorer les prévisions et contraindre le modèle numérique vers les observations disponibles sur le site d'étude, un système l'assimilation de données basé sur le 'nudging' et permettant de prendre en compte des données à haute fréquence temporelle, le "nudging direct et rétrograde" (BFN pour 'Back and Forth Nudging'), a été mis en place. Nous avons appliqué cet algorithme aux équations de Lorenz pour con rmer le comportement de cette méthode par rapport à des résultats publiés antérieurement avec d'autres méthodes d'assimilation de données. Les résultats encourageants, ont conduit à l'introduction du BFN dans Méso-NH. Nous avons mis en place des simulations avec assimilation de données simulées dans des conditions idéalisées qui ont montré une réponse cohérente du modèle numérique à l'introduction de pro ls verticaux de vent.

L'étude se concentre sur l'analyse de la dynamique de la Couche Limite Atmosphérique stable et hivernale en vallée, notamment dans des situations conduisant à des épisodes de pollution particulaire. L'analyse est effectuée à partir de simulations numériques de type LES et de mesures locales. Dans une première partie, on étudie les vents catabatiques prenant naissance sur les pentes des vallées encaissées. On montre qu'ils sont instationnaires, inhomogènes et turbulents, que leur diffusivité turbulente évolue quadratiquement avec le nombre de Froude et décroît avec la stratification ambiante. On s'intéresse ensuite au champ d'ondes internes généré par ces vents. On trouve que sa fréquence ne dépend que de la stratification et est indépendante des caractéristiques des vents qui lui donnent naissance. Dans une seconde partie, on s'intéresse au système de vents de vallées et aux inversions thermiques dans les vallées Grenobloises. Les conditions météorologiques conduisant aux épisodes de pollution sont étudiées et leur lien avec les mécanismes précédents est explicité. On montre que ces épisodes se déroulent toujours dans un contexte anticyclonique, sont induits par la présence d'une inversion et que leur évolution est liée à celle des régimes de temps. Pendant ces épisodes, le système de vents local est toujours similaire, indépendant du régime synoptique et constitué de vents thermiques, dont l'organisation spatiale est gouvernée par la géométrie du site. Ces courants sont contenus dans la couche d'inversion qui persiste pendant toute la durée de l'épisode et n'est pas détruite si l'énergie solaire est insuffisante. Le seuil énergétique correspondant est mis en évidence.

C'est a l'échelle de la turbulence que se produit l'entraînement d'air sec du flux d'Harmattan a l'intérieur du flux humide de Mousson. Ce processus de petite échelle est analysé dans le cadre du programme AMMA (Analyse Multidisciplinaires de la Mousson Africaine) qui vise a mieux renseigner et prévoir la mousson de l'Afrique de l'ouest. Les mesures atmosphériques réalisées par l'avion de recherche français ATR-42 et la modélisation des grands tourbillons (LES) sont complémentaires et permettent de documenter la structure verticale, moyenne et turbulente de la couche limite Sahélienne, de d'écrire les intrusions d'air sec et leur contribution dans les transferts turbulents, de quantifier le processus d'entraînement, de tester les paramétrisations existantes et de les relier aux conditions de grande échelle, et aux caractéristiques de la couche limite et de ses interfaces, notamment du cisaillement de vent Mousson/Harmattan.