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

L'objectif de cette thèse a été d'analyser la dynamique de l'interaction entre les ondes baroclines (ou rail des dépressions) et la variabilité basse fréquence de l'atmosphère aux moyennes latitudes. Deux approches distinctes ont été suivies pour étudier le rôle des déferlements d'ondes baroclines sur les courants-jets, l'une reposant sur les données réanalysées et l'autre sur des simulations numériques. La première partie de la thèse a plus précisément consisté à étudier le lien entre les déferlements d'ondes de Rossby et les quatre régimes de temps sur l'Atlantique Nord en hiver en utilisant les réanalyses ERA40. Le calcul des fréquences d'occurrence des déferlements d'ondes a montré que ceux-ci tendent généralement à renforcer les régimes sauf le blocage scandinave qui est détruit par du déferlement cyclonique au sud du Groenland. Ensuite, les précurseurs des transitions entre régimes de temps ont été identifiés. Le premier précurseur est relié à la propagation linéaire d'anomalies basse fréquence (période supérieure à 10 jours). Ce précurseur n'est pas systématique mais il survient durant la transition du régime zonal vers le blocage environ une semaine avant ce dernier où il prend la forme d'un train d'ondes quasi-stationnaire excité par des anomalies convectives dans l'Atlantique subtropical. Le second précurseur plus systématique intervient au niveau des interactions non-linéaires entre les tourbillons transitoires haute et basse fréquences et a pu être relié aux déferlements d'ondes. La formation et la destruction du blocage scandinave ont ensuite été plus particulièrement étudiées en analysant respectivement les transitions préférentielles du régime zonal au blocage et du blocage vers l'anticyclone groenlandais en lien avec les dépressions de surface et les déferlements d'ondes. Les dépressions de surface atteignent les mêmes intensités pendant la formation et la destruction du blocage mais ne suivent pas les mêmes trajectoires. Pendant la formation du blocage, les dépressions de surface ont des trajectoires rectilignes se dirigeant vers le nord de la Scandinavie et sont liées à un déferlement anticyclonique. Pendant la destruction du blocage, les trajectoires des dépressions de surface sont courbées sur l'Atlantique Nord en direction du Groenland et sont reliées à du déferlement cyclonique qui favorise ainsi l'apparition de l'anticyclone groenlandais. Notre analyse suggère que cette différence de comportement provient de la forme de l'écoulement basse fréquence qui n'est pas le même pendant la formation et la destruction du blocage et qui favorise un certain type de déferlement plutôt qu'un autre. Dans la seconde partie de la thèse, nous avons analysé le lien entre les températures de surface de la mer (SSTs) et le comportement du rail des dépressions avec le modèle de circulation générale de l'atmosphère Arpège-Climat en mode aquaplanète pour mieux comprendre les téléconnexions telles que l'Oscillation Arctique et/ou l'Oscillation Nord-Atlantique d'un point de vue idéalisé. / This thesis aims at analyzing the dynamics of the interaction between baroclinic waves (stormtrack) and the atmospheric low-frequency variability at midlatitudes. Two different approaches have been followed to study the impact of baroclinic wave breakings on jet-streams, one using reanalysis data and the other numerical simulations of a climate model. The first part of the PhD dealt with the link between Rossby wave breakings and the four weather regimes over the North Atlantic in winter using ERA40 reanalysis. The calculation of wave breaking frequencies showed that wave breakings tend to reinforce weather regimes except the Scandinavian blocking which is destroyed by cyclonic wave breaking south of Greenland. Then, precursors of weather regime transitions have been identified. The first precursor is linked to the linear propagation of low-frequency anomalies (period greater than ten days). This is not a systematic precursor but it occurs during the zonal to blocking transition about one week before this latter and is related to a quasi-stationary wave train excited by convective anomalies in the North Atlantic subtropics. The systematic second precursor is related to non-linear transient eddy interactions and has been linked to Rossby wave breakings. The link between the surface cyclones, Rossby wave breakings and the formation and decay of the Scandinavian blocking has been more precisely studied through the preferential transitions from the zonal weather regime to the blocking and from the blocking to the Greenland anticyclone. During the formation and decay of the blocking, surface cyclones reach the same intensities but do not follow the same trajectories. During the blocking formation, surface cyclones follow straight trajectories toward the north of Scandinavia and are linked to an anticyclonic wave breaking. Whereas during the blocking decay, surface cyclones trajectories are curved over the North Atlantic toward Greenland and are linked to a cyclonic wave breaking favouring the Greenland anticyclone formation. Our study suggests that this difference of behavior comes from the shape of the low-frequency flow which is not the same during the formation and the decay of the blocking and which can favour a particular type of wave breaking rather than another. The second part dealt with the link between sea surface temperatures (SSTs) and the storm-track in the atmospheric general circulation model Arpège-Climat in aquaplanet mode to better understand teleconnections such as the Arctic Oscillation and/or the North Atlantic Oscillation from an idealized point of view. We performed a sensitivity analysis of the eddy-driven jet variability to various stationary SST profiles.

Les précipitations sont un phénomène dont la variabilité s'étend sur une très large gamme d'échelles : de l'échelle millimétrique de la goutte d'eau (échelle microphysique) à l'échelle des circulations atmosphériques globales (échelle synoptique). Il n'existe pas de système unique capable de fournir des observations des précipitations couvrant toutes ces échelles. Les observations satellite sont celles qui actuellement résolvent le plus efficacement les grandes échelles spatiales et temporelles : de la méso-échelle à l'échelle synoptique. Dans cette thèse, nous explorons en zone tropicale les capacités des satellites à résoudre les échelles spatiales de l'ordre de 100km, jusqu'aux échelles kilométriques ; et les échelles temporelles comprises entre 24 heures et 15 minutes (afin de résoudre le cycle diurne). L'approche retenue est physico-statistique. Si les grandes échelles peuvent être résolues par des approches déterministes combinant les mesures de multiples instruments spatiaux, plusieurs facteurs limitent la pertinence des approches déterministes à fine échelle : - Les limites instrumentales en terme de résolution spatiale. - Le nombre d'instruments en orbite qui limite la fréquence d'échantillonnage des mesures. - La nature dynamique de la variabilité fine échelle. En particulier, aux fines échelles, c'est la difficulté à parfaitement localiser les structures précipitantes qui entraine les erreurs d'estimation les plus importantes. L'approche physico-statistique est ici synonyme de déterministe (pour les grandes échelles) et probabiliste (pour les fines échelles). Le premier objectif de cette thèse est de déterminer précisément la limite des échelles qui peuvent être résolues de façon déterministe. L'approche physico-statistique de l'estimation des intensités de précipitation est implémentée dans cette thèse à partir d'une méthode multicapteur déterministe pré-existante : l'algorithme TAPEER, développé dans le cadre de la mission Megha-Tropiques, qui fournit une estimation du cumul pluviométrique journalier à une résolution de 1°. C'est la génération d'ensembles désagrégés par une méthode stochastique multi-échelle qui a été retenue ici. Les ensembles sont contraints par une information fine échelle : un masque de détection des aires précipitantes dérivé des images infrarouge metosat-SG à une résolution de 3km (et avec une image toutes les 15 minutes). La génération d'ensemble permet de caractériser l'incertitude sur l'estimation à travers la dispersion des réalisations de l'ensemble. Chaque réalisation de l'ensemble est générée de façon à reproduire le plus fidèlement possible les propriétés statistiques (distribution de fréquence des intensités, autocorrélation spatiale et temporelle) des véritables champs de précipitation. Ces champs et cette technique ont un apport pour les applications hydrologiques, par exemple pour améliorer le ruissellement lié aux précipitations intenses dans les modèles. Considérer plusieurs réalisations permet de plus d'étudier la propagation des incertitudes à travers un modèle. / Rainfall variability involves a wide range of scales: from the millimeter-scale associated with microphysics to the synoptic scale of the global atmospheric circulation. No existing observation system is able to cover all these scales by itself. Satellite-based observation systems are currently the most efficient systems to resolve the large spatial and temporal scales: from mesoscale meteorology to the synoptic scale. This thesis is dedicated to the exploration of satellites ability to resolve spatial scales from 100km to 2km and temporal scales from 24h to 15 min (in order to resolve the diurnal cycle). The chosen approach is both physical and statistical (or deterministic and probabilistic). The idea is that the deterministic approach can resolve the large scales, but several factors limit its relevance when dealing with fine scales: -The limited resolution of the instruments. -The number of orbiting instruments that limits temporal sampling. -The dynamic nature of fine scale variability. At fines scales, most of the errors in rainfall estimation from satellite comes from not perfectly localizing the precipitating cells. The first objective of this thesis is to identify precisely the lowest limit in scale where the deterministic approach is appropriate. The implementation of the physical-statistical approach relies on an existing multisensor estimate of daily precipitation at a 1° resolution: the TAPEER algorithm developed as part of the Megha-Tropiques mission. The chosen method is a hybrid physical disaggregation and stochastic downscaling via a multiscale representation. The result is an ensemble of high-resolution probable realizations of the rain intensity field. The ensemble is constrained by a high resolution rain detection mask derived from meteosat-SG infrared images at 3km resolution (one image every 15 minutes). The uncertainty associated with the final estimation is handled through the ensemble dispersion. Every realization is generated so that its statistical properties (frequency distribution of the intensities, autocorrelation function) mimic those of the true rain field. The generated fields and the proposed technique contribute to hydrological applications for instance by improving the runoff associated to high precipitation rates in models. Using several realizations is a way to study uncertainty propagation through a model.

Le modèle SAFRAN-ISBA-MODCOU est évalué et sa physique améliorée. L'analyse météorologique SAFRAN est tout d'abord validée en détail. Le modèle de surface ISBA est ensuite modifié pour mieux décrire la conductivité hydraulique dans le sol. Une stratégie de calibration est définie et appliquée à l'échelle de la France. Le modèle amélioré est finalement utilisé pour évaluer les impact du changement climatique en région méditerranéenne. Un scénario d'un modèle de climat régional est désagrégé en utilisant deux méthodes différentes. L'étude montre que les incertitudes liées à la désagrégation sont importantes. L'incertitude liée au modèle d'impact est moindre, mais doit être prise en compte pour les extrêmes. Dans cette région, les débits extrêmes, et parfois, moyens, augmenteront pendant la première moitiée du XXIème siécle. A la fin du siècle, le scénario indique une baisse des débits moyens et une stabilité des débits extrêmes, ce qui conduit à une augmentation de la variabilité.

La zone Cévennes-Vivarais et son extension jusque la Méditerranée a été choisie par l'OHM-CV (Observatoire Hydrométéorologique-Méditerranéen des Cévennes-Vivarais) pour la fréquence et la représentativité des événements de pluies intenses qu'elle subit. Depuis 2002, des campagnes automnales de mesures GPS y sont menées dans le but d'améliorer notre connaissance du champ de vapeur d'eau troposphérique associé aux événements de pluie intense. Les observations GPS météorologiques sont les ZTD (délais troposphériques au zénith) et les gradients de délais. Exprimé par une composante Est-Ouest (Gew) et une composante Nord-Sud (Gns), le gradient de délai traduit l'anisotropie du champ de vapeur d'eau à proximité du site GPS.Pour le traitement des observations GPS des réseaux permanents et temporaires de la région de l'OHM-CV, une configuration optimale a été recherchée à partir de différents tests de sensibilité. Une précision de 5 mm sur les ZTD peut être obtenue avec cette stratégie d'analyse pour le réseau régional GPS de la région de l'OHM-CV ; Gew est précis à 6 mm près, alors que Gns apparaît moins précis (à 12 mm près) du fait de l'absence de satellite survolant le pôle Nord.Afin de préparer l'assimilation des observations GPS par les systèmes d'assimilation de données à méso-échelle de Météo-France (modèle de Météo-France AROME), des simulateurs d'observations GPS (ZTD, STD et gradients) ont été implémentés dans le modèle non-hydrostatique à haute résolution Méso-NH. Le laboratoire numérique que constitue le modèle Méso-NH a été un moyen efficace de quantifier la sensibilité des simulations de délais à différentes formulations de la réfractivité atmosphérique.

Pour mieux comprendre les structures de la couronne solaire de température élevée et leurs évolutions, des méthodes de vision à 3 dimensions à partir des images EUV de SOHO/EIT ont été développées. Des anaglyphes ont servi à mieux comprendre le rôle du champ magnétique et des intéractions entre champs<br />fermés et ouverts. Il est apparu que la reconnection entre des trous coronaux (CH) et des régions actives (AR) peut expliquer certaines disparitions de boucles coronales fermées. De plus, des analyses multi-instruments ont montré que la structuration des<br />CHs résulte parfois de changements dans l'activité magnétique voisine.<br />Les mesures d'incertitudes sur les reconstructions par stéréovision n'ont pas permis de valider une méthode générale pour visualiser la couronne globalement en 3D. Par contre, le principe a pu être utilisé avec succès pour la reconstruction en 3 dimensions de boucles coronales. Ainsi après analyse des paramètres physiques de 30 boucles EUV de températures intermédiaires d'une région active, il est apparu que contrairement aux boucles chaudes en X, elles étaient en équilibre hydrostatique.<br />Les boucles d'une région active émergentes sont apparues comme se détorsadant au fur et à mesure de son grandissement ce qui correspond à un transfert de l'hélicité.<br />Des corrélations entre des brusques détorsadages et des éruptions ont aussi été établis dans les cas où la torsion initiale est trop importante et permet le développement<br />d'instabilités.<br />Ces études 3D vont pouvoir permettre de mieux contraindre les bilans énergétiques pour le chauffage coronal gràce à la conservation de l'hélicité et d'améliorer la prévision de la météorologie spatiale. L'utilisation de techniques d'imageries adaptées (comme<br />le dégrillage ou le Modèle de Vision Multi-échelle) a permis de faire ressortir certains détails dans la formation de CMEs par exemple.

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.

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

Les modèles emboîtés sont largement utilisés en météorologie et en océanographie. Ils permettent un accroissement local de la résolution, dans les zones où cela semble nécessaire, via l'intégration d'un même modèle sur une hiérarchie de grilles. Dans le cas d'interaction one-way, les conditions aux frontières pour la grille fine proviennent d'une interpolation de la solution obtenue sur la grille à plus faible résolution. Dans le cas d'interaction two-way, une rétroaction de la grille fine vers la grille grossière est ajoutée. Toutefois, le problème de l'assimilation variationnelle de données dans de tels systèmes n'a pas, ou peu, été étudié à ce jour. Ces classes de méthodes, notamment l'algorithme 4D-Var, permettent d'améliorer la solution d'un modèle, jusqu'ici mono-grille, en minimisant une fonctionnelle mesurant l'écart de ce modèle aux observations présentes sur une fenêtre temporelle. Le travail présenté ici vise donc à formuler un algorithme d'assimilation 4D-Var localement multi-grille. Pour le cas général d'une grille haute résolution emboîtée localement dans une autre à plus faible résolution, nous posons les équations du système adjoint dans les deux cas d'interactions one-way et two-way. Nous montrons ainsi que la formulation adjointe fait naturellement apparaître de nouvelles interactions entre les grilles, dans le sens opposé de celles existant dans la formulation directe. De plus, nous proposons différentes variantes à ces algorithmes, réalisant un couplage faible entre les solutions des différents modèles via l'ajout d'un terme de contrôle au niveau des transferts inter-grilles. Nous présentons également l'application d'une méthode multi-grille, le Full Approximation Scheme, à l'assimilation variationnelle de données. Cette approche permet d'obtenir un algorithme d'assimilation multi-grille potentiellement très efficace. Enfin, ces méthodes sont testées sur le cas d'un modèle Saint Venant 2D. Nous constatons une réduction importante des erreurs des solutions multi-grilles, ainsi qu'une accélération de la convergence de ces algorithmes.

Le changement climatique en cours et à venir, qui se traduit au Canada comme dans la plupart des régions du monde par un climat plus chaud, peut être parfois amplifié en fonction des conditions de surface et notamment dans les zones urbanisées durant la saison estivale. Compte tenu que la moitié de l'humanité vit désormais en ville, l'urbanisation est un enjeu dans le réchauffement global. La grande métropole de Montréal au sud du Québec est une région soumise à une urbanisation propice à la formation d'îlots de chaleur urbains (ICU) qui ont une influence sur le climat régional. La présente étude analyse durant des périodes estivales l'évolution du comportement thermique de Montréal au cours des 20 dernières années (1984-2005), afin d'en dresser le bilan à l'aide d'une série d'images satellites Landsat 5 TM. Suite au traitement de dix images Landsat 5, la localisation et l'amplitude des ICU sur le territoire ont été établies. Les secteurs fortement minéralisés et/ou végétalisés ont été identifiés en explorant la relation entre les températures de surface et l'indice de végétation Normalized Difference Vegetation Index (NDVI). Ainsi on a pu constater que certaines municipalités et/ou arrondissements (ou portions du territoire) présentent des bilans thermiques de surface forts différents. Certains ont vu leur couvert végétal disparaître peu à peu à l'instar de nouvelles constructions, faisant croître la superficie au sol et l'intensité des ICU. D'autres ont en revanche concilié urbanisme et environnement en maintenant les espaces verts, ce qui a limité voire diminué la chaleur à certains endroits. Les résultats montrent aussi que les îlots de chaleur se sont amplifiés au cours des décennies quelque soit le mois ou l'année, en particulier dans la moitié sud-ouest de l'île de Montréal. Ce comportement thermique résulte donc d'une distribution des surfaces végétalisées versus non (ou peu) végétalisées. Dans cette recherche nous démontrons que la corrélation entre ces températures de surface et le NDVI est très forte avec un gradiant allant de températures fraiches pour les zones vertes vers des valeurs beaucoup plus chaudes pour les secteurs minéralisés. ______________________________________________________________________________ MOTS-CLÉS DE L’AUTEUR : Îlot de chaleur urbain, Télédétection, Infrarouge thermique, Comportement thermique, Indice de végétation NDVI, Analyse en composante principale.