Les études sur l'Antarctique sont motivées par l'interprétation des données que renferme la glace ancienne sur les climats passés et sur les possibles interactions entre la calotte et un éventuel changement climatique. Les radiomètres micro-ondes sont des instruments bien adaptés pour l'étude de la neige, car ils sont indépendants des conditions nuageuses et d'éclairement. Lorsque la neige est sèche, ils sont à même de fournir des informations en profondeur. Les caractéristiques du manteau neigeux (température, densité et taille de grains) sont liées aux processus de dépôt et leur évolution dépend des conditions climatiques. Ce sont ces caractéristiques, déterminantes vis à vis de la signature spectrale de la neige, que nous cherchons à évaluer à partir des données de télédétection. Les données utilisées proviennent du Scanning Multichannel Microwave Radiometer à bord du satellite Nimbus 7 (1978-86). Dans un premier temps, ces mesures ont été comparées à des données de terrain. A une forte stratification du manteau neigeux correspond une grande différence entre les polarisations verticale et horizontale, surtout pour les fréquences 6.6 et 10.7 GHz. On a trouvé aussi une forte corrélation entre la valeur moyenne de la taille des grains sur 0-2 mètres et la différence de comportement entre les fréquences 18 et 6.6 GHz (gradient en fréquence). Dans une seconde étape, nous avons modélisé l'émissivité microonde de la neige. Le modèle est basé sur la résolution des équations de Maxwell au travers de la théorie des fortes fluctuations (Stogryn 1986). Une solution analytique a été développée dans le cas d'un milieu uniforme en profondeur. Lorsque les caractéristiques de la neige changent en fonction de la profondeur, la solution fait appel à une résolution numérique. Le manteau neigeux est considéré isotherme, stratifié horizontalement et isotrope à l'intérieur de chaque strate. Les résultats du modèle font apparaitre plus clairement les contributions de la densité de la neige, de la taille des grains et surtout de la stratification sur les signatures spectrales.

La télédétection satellitaire est en mesure de fournir des informations globales sur les calottes polaires et d'en effectuer un suivi régulier dans le temps. Pour être facilement interprétées, afin d'en déduire les caractéristiques du manteau neigeux observé depuis l'espace (taille, forme de grains, rugosité de surface ... ), les données satellitaires doivent être validées et inversées à l'aide de paramétrisations simples. Il est alors indispensable de définir des paramètres robustes et simples de la taille et de la forme des grains de neige à partir de leur observation. Un moyen d'observation et de mesure est le traitement d'images gui a permis d' établir ces relations et qui a permis un traitement automatique d'un grand nombre de données indépendamment de l'observateur. Un autre problème glaciologique, est l'interprétation des données issues de l'analyse des gaz piégés dans les bulles d'air. Cette étude implique, en particulier, la datation de la glace dans le névé lors de la fermeture des pores, indispensable pour déterminer dans la glace profonde l'âge du gaz par rapport à la glace qui l'emprisonne. En effet la datation des gaz contenus dans les bulles des carottes de glace s'effectue à partir de l'âge de la glace, et par conséquent, elle néces ite la connaissance de la différence d'âge entre les gaz et la glace lors de la fermeture des pores. Ces datations nécessitent l'utilisation des modèles de densification, de fermeture des pores, et de diffusion des gaz dans le névé. Les carottages réalisés dans les régions centrales de l'Antarctique et du Groenland permettent d'obtenir des données expérimentales essentielles pour caractériser le névé polaire. Le traitement automatique des images issues de ces données a fourni les paramètres de structure introduits dans ces modèles. Un autre aspect étudié au Laboratoire de Glaciologie est l'augmentation de la taille moyenne des cristaux de glace avec la profondeur ou l'âge de la glace. Les analyses des lames minces de glace permettant de déterminer ces lois de grossissement des grains ont été réalisées jusgu'à présent manuellement. Le traitement automatique d' images 2D a permis d' analyser les lois de grossissement des grains sur d' autres sites et a permis d' étudier les microstructures en terme d' évolution des distributions de taille et en terme d'évolution de la topologie des grains.

Le nitrate, produit de la fin de chaîne de réaction des oxydes d'azotes de l'atmosphère (NOx = NO +NO2), est l'un des ions le plus abondant de la neige et de la glace polaire. Ses rapports isotopiques stables (δ18O, δ15N et Δ17O) ont été abondamment utilisés pour contraindre ses sources et les chemins réactionnels. De plus, le nitrate archivé dans les carottes de glace profondes peut apporter de nouvelles contraintes sur les conditions climatiques passées. Cependant, le dépôt de nitrate dans les régions polaires à faible accumulation est réversible en raison des processus post-dépôts, compliquant l'interprétation des enregistrements. Actuellement, il existe des enregistrements de nitrate issus de carottes de glace profonds couvrant de l'information climatique sur plusieurs milliers d'années dont leur interprétation dépend d'une quantification précise ces phénomènes post-dépôts. Nous avons étudié expérimentalement le transfert d'excès-17O de l'ozone durant la réaction en phase gaz de NO2 + O3  NO3 + O2, qui est une réaction importante de la chimie nocturne de formation du nitrate. De cette étude nous avons déterminé la fonction de transfert du 17O donnée par : ∆17O(O3*) = (1.23 ± 0.19) × ∆17O(O3)bulk + (9.02 ± 0.99). Nous avons aussi évalué la distribution intramoléculaire des isotopes de l'oxygène de l'ozone et observé que l'excès d'enrichissement résidait de manière prépondérante sur les atomes terminaux de l'ozone. Ces résultats auront une implication importante sur la compréhension de la formation du nitrate via les mécanismes d'oxydation des précurseurs NOx. L'impact de la photolyse sur les concentrations et les compositions isotopiques stables du nitrate est étudié dans ce travail de thèse sur la base d'étude de laboratoire et de terrain. Une étude de laboratoire a été conduite en irradiant de la neige naturelle de Dôme C avec une lampe UV à xénon et en utilisant différents filtres UV (280 nm, 305 nm et 320 nm). Sur la base des mesures des rapports isotopiques de l'oxygène et de l'azote, la dépendance aux longueurs d'onde des fractionnements isotopiques a été déterminée. En conséquence, en présence de lumière UV de haute énergie, le fractionnement isotopique est décalé vers des valeurs moins négatives et vice versa. Sur la base des fractionnements isotopiques obtenus en laboratoire, nous avons dérivé un décalage apparent de la valeur du zéro point d'énergie (ZPE) qui apporte une meilleure contrainte sur la section efficace d'absorption du 15NO3-. Ce décalage apparent est obtenu en minimisant les écarts entre les observations et les fractionnements isotopiques calculés basés sur un modèle de décalage ZPE, modèle qui inclut outre le décalage ZPE, le changement des largeurs, de l'asymétrie et de l'amplitude des sections efficaces d'absorption lors de la substitution isotopique. Nous avons validé le nouveau ZPE apparent en conduisant une étude de terrain à Dôme C, Antarctique. Dans cette étude, un dispositif expérimental a été construit sur le site et l'effet du rayonnement solaire UV sur la photolyse du nitrate de la neige investigué. Cette étude était basée sur la comparaison de deux puits remplis par de la neige soufflée homogénéisée dont l'un des deux puits n'était pas soumis aux rayonnements UV. Le fractionnement isotopique de 15N pour la neige exposée aux UV (-67.9 ± 12 ‰) est en bon accord avec le modèle de décalage ZPE estimé au cours de ce travail de thèse (-55.4 ‰). Ces valeurs sont aussi dans la gamme des fractionnements isotopiques apparents observée précédemment au Dôme C. Nous pensons que l'inclusion des ces nouvelles connaissances dans un modèle prédisant l'enrichissement du 15N dans les carottes de glace permettra une interprétation quantitative de l'information préservée dans la glace.

L'objectif de cette étude est d'évaluer à l'aide des techniques de phylogénie moléculaire la diversité bactérienne non cultivable dans le manteau neigeux du Mont Blanc (MtBl) contenant de la poussière saharienne déposée au cours de quatre événements pendant la période de 2006 à 2009. Le but final est de découvrir les bactéries qui pourraient être impliquées dans l'établissement du microbiote neigeux.

La neige est un matériau poreux dont la microstructure change en permanence. L'ensemble de ces transformations, qui prend le nom de ''métamorphisme", est susceptible d'affecter les propriétés thermiques, mécaniques et électromagnétiques de la neige au niveau macroscopique. En particulier, les échanges d'énergie et de matière à l'intérieur du manteau neigeux et entre la neige et l'atmosphère sont fortement influencés par l'évolution au cours du temps de la microstructure de la neige. Une représentation adéquate du métamorphisme dans les modèles de manteau neigeux s'avère donc cruciale. La microstructure d'un matériau poreux peut être raisonnablement décrite en se servant d'un nombre réduit de variables. En effet, la masse volumique, la surface spécifique (SSA) et la distribution de courbure permettent de caractériser la microstructure d'un matériau. Cependant, dans le cas de la neige cette approche n'en est qu'à ses débuts et n'a pas encore été appliquée de façon systématique. Des variables semi-empiriques, difficiles à mesurer et dépourvues de lien direct avec d'autres propriétés physiques, sont encore largement utilisées dans les modèles détaillés de manteau neigeux. Ce travail de thèse s'inscrit dans cette tentative de représenter la microstructure de la neige au cours du temps à l'aide de variables bien définies et mesurables sur le terrain. Parmi ces variables, nous nous sommes attachés notamment à la SSA, qui constitue une grandeur essentielle pour l'étude du manteau neigeux et de son évolution temporelle. Différentes lois d'évolution de la SSA ont été étudiées, à partir de relations empiriques basées sur des ajustements de données expérimentales jusqu'aux modèles physiques qui représentent le flux de la vapeur d'eau entre les grains de neige. Ces lois ont été dans un premier temps testées à l'aide d'un modèle simplifié de manteau neigeux et puis introduites directement dans le modèle SURFEX/ISBA-Crocus. Pour ce faire, la SSA dans Crocus a été transformée en variable prognostique, en remplaçant d'autres variables semi-empiriques préexistantes. Les différentes formulations de l'évolution temporelle de la SSA ont été comparées à des mesures de terrain, acquises lors de deux campagnes à Summit (Groenland) et au Col de Porte (France). Ces mesures ont été effectuées en utilisant de nouvelles techniques optiques et ont permis d'obtenir un riche jeu de données avec une grande résolution verticale. Les résultats montrent que les différentes formulations sont comparables et reproduisent bien les mesures, avec un écart quadratique moyen entre les valeurs de SSA simulées et observées inférieur à 10 m^2/kg. Enfin, nous avons contribué à faire le pont entre la microstructure de la neige et ses propriétés macroscopiques. En particulier, nous nous sommes intéressés au lien entre, d'une part, la SSA et, d'autre part, les propriétés mécaniques et optiques. Dans le premier cas, nous avons investigué la corrélation entre la SSA et la résistance à l'enfoncement mesurée avec un Snow Micro Pen (SMP). Les résultats encore préliminaires semblent indiquer que la SSA peut être dérivée de la masse volumique et de grandeurs micro-mécaniques estimées à partir du signal du SMP avec un modèle statistique. Dans le deuxième cas, nous avons simulé l'albédo de surface à Summit à partir des profils mesurés de masse volumique et de SSA et du contenu en impuretés. Les résultats de cette étude ont démontré que l'albédo spectral peut être correctement simulé à l'aide d'un modèle de transfert radiatif et l'énergie absorbée par le manteau neigeux peut être estimée avec une précision d'environ 1%.

La prédétermination de la hauteur de départ des avalanches représente un défi majeur pour l'évaluation du risque en montagne. Cette hauteur constitue en effet un ingrédient d'entrée important des procédures de zonage et de cartographie du risque. Dans cette thèse, nous présentons un formalisme rigoureux dans lequel les distributions de hauteur de départ d'avalanche sont exprimées à travers un couplage des facteurs mécaniques et météorologiques. Le critère de stabilité du système plaque - couche fragile est étudié en utilisant une analyse mécanique par éléments finis prenant en compte l'hétérogénéité spatiale des propriétés mécaniques. Considérant qu'une avalanche ne peut se produire que si la hauteur de chute de neige dépasse une hauteur critique correspondant au critère de stabilité, les distributions de hauteur de départ obtenues à partir du modèle mécanique sont couplées avec la distribution des chutes de neige extrêmes sur 3 jours. Nous montrons que ce modèle couplé est capable de reproduire des données de terrain de 369 avalanches naturelles de plaque à La Plagne (France). Non seulement la queue de la distribution en loi puissance, correspondant à des épaisseurs de plaque élevées, mais aussi le corps de la distribution pour les plaques moins épaisses, sont bien reproduits par le modèle. Les avalanches petites à moyennes semblent être essentiellement contrôlées par la mécanique, tandis que les grosses avalanches et l'exposant de la loi puissance associé, sont influencés par un couplage mécanique - météorologique fort. Par ailleurs, nous démontrons que la distribution obtenue est fortement dépendante de l'espace, et, en utilisant les processus max-stables permettant une interpolation spatiale rigoureuse, notre modèle couplé est utilisé pour obtenir des cartes de hauteur de départ d'avalanche pour différentes périodes de retour sur l'ensemble des Alpes françaises.

La région de l'archipel canadien, située en Arctique, connaît actuellement d'importants changements climatiques, se traduisant notamment par une augmentation des températures, une réduction de l'étendue de la banquise marine et du couvert nival terrestre ou encore une perte de masse significative des calottes glaciaires disséminées sur les îles de l'archipel. Parmi ces calottes glaciaires, la calotte Barnes, située en Terre de Baffin, ne fait pas exception comme le montrent les observations satellitaires qui témoignent d'une importante perte de masse ainsi que d'une régression de ses marges, sur les dernières décennies. Bien que les calottes glaciaires de l'archipel canadien ne représentent que quelques dizaines de centimètres d'élévation potentielle du niveau des mers, leur perte de masse est une composante non négligeable de l'augmentation actuelle du niveau des mers. Les projections climatiques laissent à penser que cette contribution pourrait rester significative dans les décennies à venir. Cependant, afin d'estimer les évolutions futures de ces calottes glaciaires et leur impact sur le climat ou le niveau des mers, ils est nécessaire de caractériser les processus physiques tels que les modifications du bilan de masse de surface. Cette connaissance est actuellement très limitée du fait notamment du sous-échantillonnage des régions arctiques en terme de stations météorologiques permanentes. Une autre particularité de certaines calottes de l'archipel canadien, et de la calotte Barnes en particulier, est de présenter un processus d'accumulation de type glace surimposée, ce phénomène étant à prendre en compte dans l'étude des processus de surface. Pour palier au manque de données, l'approche retenue a été d'utiliser des données de télédétection, qui offrent l'avantage d'une couverture spatiale globale ainsi qu'une bonne répétitivité temporelle. En particulier les données acquises dans le domaine des micro-ondes passives est d'un grand intérêt pour l'étude de surfaces enneigées. En complément de ces données, la modélisation du manteau neigeux, tant d'un point de vue des processus physiques que de l'émission électromagnétique permet d'avoir accès à une compréhension fine des processus de surface tels que l'accumulation de la neige, la fonte, les transferts d'énergie et de matière à la surface, etc. Ces différents termes sont regroupés sous la notion de bilan de masse de surface. L'ensemble du travail présenté dans ce manuscrit a donc consisté à développer des outils permettant d'améliorer la connaissance des processus de surface des calottes glaciaires du type de celles que l'on rencontre dans l'archipel canadien, l'ensemble du développement méthodologique ayant été réalisé sur la calotte Barnes à l'aide du schéma de surface SURFEX-CROCUS pour la modélisation physique et du modèle DMRT-ML pour la partie électromagnétique. Les résultats ont tout d'abord permis de mettre en évidence une augmentation significative de la durée de fonte de surface sur la calotte Barnes (augmentation de plus de 30% sur la période 1979-2010), mais aussi sur la calotte Penny, elle aussi située en Terre de Baffin et qui présente la même tendance (augmentation de l'ordre de 50% sur la même période). Ensuite, l'application d'une chaîne de modélisation physique contrainte par diverses données de télédétection a permis de modéliser de manière réaliste le bilan de masse de surface de la dernière décennie, qui est de +6,8 cm/an en moyenne sur la zone sommitale de la calotte, qui est une zone d'accumulation. Enfin, des tests de sensibilité climatique sur ce bilan de masse ont permis de mettre en évidence un seuil à partir duquel cette calotte voit disparaître sa zone d'accumulation. Les modélisations effectuées suggèrent que ce seuil a de fortes chances d'être atteint très prochainement, pour une augmentation de température moyenne inférieure à 1°C, ce qui aurait pour conséquence une accélération de la perte de masse de la calotte.

L'albédo, fraction de rayonnement réfléchi dans le spectre solaire, est une variable clef du bilan énergétique des surfaces enneigées et englacées. Cette grandeur possède une forte variabilité spatio-temporelle ce qui fait de la télédétection un outil adapté pour son étude. L'albédo dépend à la fois des propriétés physiques du milieu considéré et des caractéristiques du rayonnement incident. Les différentes grandeurs liées à l'albédo sont fonction des domaines angulaires et spectraux des radiations considérées. Les mesures de répartition angulaire du rayonnement réfléchi par la neige ont montré que l'hypothèse lambertienne pouvait conduire à des erreurs non négligeables lors de la détermination de l'albédo par télédétection. La connaissance des caractéristiques de la répartition angulaire du rayonnement réfléchi par la neige permet de développer une nouvelle méthode de détermination de l'albédo en zones montagneuses. Cette méthode prend en compte les effets liés à la forte variabilité topographique des terrains de montagne, à l'anisotropie du rayonnement réfléchi par la neige et par la glace ainsi que les variations spectrales de l'albédo en fonction des propriétés physiques de la surface. Elle a été appliquée à deux types de données : des photographies terrestres visibles et proche infrarouges (résolution spatiale 10 m) et des images MODIS (résolution spatiale 250 m). L'incertitude sur la valeur de l'albédo ainsi déterminée est évaluée à ±10% grâce aux mesures de terrain effectuées sur le glacier de Saint Sorlin (massif des Grandes Rousses, France). L'étude des cartes d'albédo issues de dix années (2000-2009) d'images MODIS montre qu'il n'y a pas de décroissance marquée de la valeur de l'albédo en zone d'ablation au contraire de ce qui a été prouvé pour le glacier du Morteratsch (Suisse). De plus, il existe une corrélation très élevée entre la valeur minimale de la moyenne de l'albédo sur le glacier, i.e. l'albédo moyen du glacier le jour où la ligne de neige est proche de la ligne d'équilibre, et la valeur du bilan de masse annuel spécifique. L'assimilation des données d'albédo obtenues grâce aux images MODIS et aux photographies terrestres dans le modèle de neige CROCUS permet une bonne estimation du bilan de masse spatialisé du glacier de Saint Sorlin (rmse=0.5 m w.e. pour les cinq années hydrologiques étudiées). Les forçages météorologiques utilisés pour cette étude sont de moyenne échelle. L'analyse succincte de la contribution des différents flux atmosphériques au bilan d'énergie de surface montre qu'en zone d'ablation comme en zone d'accumulation, le bilan radiatif net courtes longueurs d'ondes constitue la source principale d'énergie et que la variabilité de ce flux explique la majeure partie de la variabilité journalière de la somme des flux atmosphériques. Appliquées à d'autres glaciers, ces méthodes permettraient de savoir si les conclusions établies pour notre seul glacier d'étude sont valables pour d'autres glaciers. Elles rendraient également possibles la reconstruction du bilan de masse spatialisé sur 10 ans d'autres glaciers et potentiellement une meilleure quantification des processus physiques mis en jeu dans le bilan de masse de ces glaciers tempérés

Le bilan d'énergie de surface du Plateau Antarctique est essentiellement contrôlé par les propriétés physiques des premiers centimètres du manteau neigeux. Or l'évolution de cette neige de surface est complexe car elle dépend de processus fondamentalement imbriqués: vitesse de métamorphisme, profils de température, pénétration du rayonnement solaire, précipitations, transport de neige par le vent, etc. L'objectif de ces travaux de thèse est d'étudier ces diverses composantes et leur couplage afin de simuler l'évolution de la densité de la neige et de la taille de grain (surface spécifique) sur le Plateau Antarctique. Pour représenter de manière physique l'absorption de l'énergie solaire à l'intérieur du manteau, un modèle de transfert radiatif à fine résolution spectrale (TARTES) a été implémenté dans le modèle de manteau neigeux détaillé Crocus. TARTES permet de calculer le profil vertical d'absorption d'énergie dans un manteau stratifié dont les caractéristiques sont connues. Parmi elles, la forme des grains, explicitement prise en compte dans TARTES, a été peu étudiée jusqu'à présent. C'est pourquoi une méthode de détermination expérimentale de la forme optique des grains est proposée et appliquée à un grand nombre d'échantillons de neige. Cette méthode, basée sur des mesures optiques, des simulations TARTES, et l'inférence bayésienne, a permis de déterminer la forme la plus adéquate pour simuler les propriétés optiques de la neige, et a mis en évidence le fait que représenter la neige par un ensemble de particules sphériques conduisait à surestimer la profondeur de pénétration du rayonnement d'environ 30%. L'impact de l'absorption en profondeur du rayonnement sur les profils de température dans le manteau et son métamorphisme est ensuite étudié par des approches analytique et numérique, mettant en valeur la sensibilité des profils aux propriétés de la neige proche de la surface. En particulier, la densité de la neige sur les premiers centimètres est cruciale pour le bilan d'énergie du manteau car elle impacte à la fois la profondeur de pénétration du rayonnement et la conductivité thermique du manteau. Puisque le modèle Crocus tient compte de ce couplage entre propriétés optiques et physiques du manteau, il est utilisé pour estimer l'influence des conditions météorologiques sur la variabilité temporelle des propriétés physiques de la neige de surface à Dôme C. Ces simulations sont évaluées au regard d'un jeu de données collectées lors de missions de terrain et de mesures automatiques de l'albédo spectral et de la pénétration du rayonnement dans la neige. Ces observations mettent en évidence le rôle primordial des précipitations dans les variations rapides de taille de grain en surface et l'existence d'un cycle saisonnier de cette taille de grain. Ces variations sont bien simulées par Crocus lorsque le forçage atmosphérique qui lui est imposé est adéquat. En particulier, l'impact du vent sur l'évolution du manteau est fondamental car il contrôle la densité de surface par le biais du transport de neige. Ce transport est aussi à l'origine de la variabilité spatiale des propriétés de la neige observée à Dôme C. C'est pourquoi une modélisation stochastique de l'érosion et du transport de neige dans Crocus est proposée. En plus d'expliquer la variabilité spatiale de la densité et de la taille de grain, elle permet de reproduire celle de l'accumulation annuelle ainsi que les variations rapides de hauteur de neige liées à des épisodes de vent. Ces travaux ont permis une meilleure représentation des processus physiques qui contrôlent les variations des propriétés de la neige de surface à Dôme C, tout en soulignant le rôle primordial du vent, dont l'impact sur le manteau est particulièrement complexe à simuler. / The surface energy balance of the Antarctic Plateau is mainly governed by the physical properties of the snowpack in the topmost centimeters, whose evolution is driven by intricated processes such as: snow metamorphism, temperature profiles variations, solar radiation penetration, precipitation, snow drift, etc. This thesis focuses on the interactions between all these components and aims at simulating the evolution of snow density and snow grain size (specific surface area) on the Antarctic Plateau. To physically model the absorption of solar radiation within the snowpack, a radiative transfer model with high spectral resolution (TARTES) is implemented in the detailed snowpack model Crocus. TARTES calculates the vertical profile of absorbed radiation in a layered snowpack whose characteristics are given. These characteristics include snow grain shape, a parameter that has been seldom studied. For this reason, an experimental method to estimate the optical grain shape is proposed and applied to a large number of snow samples. This method, which combines optical measurements, TARTES simulations and Bayesian inference, is used to estimate the optimal shape to be used in snow optical models. In addition, it highlights that representing snow as a collection of spherical particles results in overestimation of the penetration depth of solar radiation. The influence of the penetration of solar radiation on the snow temperature profiles is then investigated with analytical and numerical tools. The results point out the high sensitivity of the temperature profiles to surface snow physical properties. In particular, the density of the topmost centimeters of the snowpack is critical for the energy budget of the snowpack because it impacts both the effective thermal conductivity and the penetration depth of light. To simulate the evolution of snow physical properties at Dome C by taking into account their interdependence with snow optical properties, the model Crocus is used, driven by meteorological data. These simulations are evaluated against a set of data collected during field missions as well as automatic measurements of snow spectral albedo and penetration depth. These observations highlight the influence of weather conditions on the temporal variability of surface snow properties. They show the existence of a slow decrease of snow grain size at the surface during summer. Rapid changes are also observed, essentially due to precipitation. These variations are well simulated by Crocus when forced by an appropriate atmospheric forcing. In particular, the impact of wind on the evolution of the snowpack is crucial because it controls the surface density through snow transport. This transport is also responsible for the spatial variability of snow properties observed at Dome C. That is why a stochastic representation of snow erosion and transport in Crocus is proposed. It explains well the observations of the spatial variability of density and grain size, and reproduces the variability of the annual accumulation as well as rapid changes in snow height resulting from drift events. This study improves our understanding of the physical processes which drive the properties of snow close to the surface on the Antarctic Plateau, and also points out the critical role of wind, the impact of which is very difficult to account for in models yet.

Les environnements montagneux sont un lieu privilégié d'échange d'eau et d'énergie. Les rivières de montagne alimentent en eau 40% de la population mondiale et sont sujettes à une pression démographique et climatique important. Dans ce contexte, la compréhension des processus météorologiques, hydrologiques et hydrogéologiques est fondamentale pour la gestion globale de la ressource en eau. L'étude, présentée dans ce manuscrit de thèse, se positionne au sein des environnements montagneux où l'hydrologie est influencée par le couvert neigeux saisonnier et par les glaciers, et propose une approche de modélisation interdisciplinaire afin d'améliorer la compréhension des processus en jeu.Aujourd'hui, si les modèles sont capables de simuler le débit sur les rivières de montagnes jaugées sous influence nivale et glaciaire, un certain nombre d'incertitudes persistent quant à l'utilisation de tels modèles hors de leur conditions de validation (en réponse à un climat différent ou sur un domaine non-jaugé). La principale source d'incertitude est liée au manque de connaissance des précipitations en montagne, dont la mesure est rare et incertaine. C'est pourtant la principale composante du bilan hydrologique. A cet égard, nous proposons d'exploiter l'information fournie par la géométrie du couvert neigeux et des glaciers, en tant que “pluviomètres géants” à l'échelle de ces réservoirs, dans un modèle hydrologique à réservoirs conceptuels reposant sur la notion de bassin versant.L'information, hydrologique, nivale et glaciaire est évaluée dans un cadre de calibration multi-objectifs. Les résultats montrent que, dans cette configuration, la validation conjointe du modèle hydrologique par le débit journalier, le bilan de masse glaciaire annuel et la hauteur de neige locale journalier permet de réduire fortement l'incertitude sur le forçage météorologique journalier et d'améliorer la robustesse du modèle. Ce résultat préliminaire nous a permis de reconstruire, en conséquence, le bilan de masse local annuel à l'échelle des glaciers.Par ailleurs, la représentation des glaciers au sein d'un modèle hydrologique pose un certain nombre de défis, surtout dans la perspective de simuler les processus hydrologiques à l'échelle pluri-annuelle. En particulier, la prise en compte de l'évolution de la géométrie des glaciers au sein d'un modèle hydrologique est balbutiante. A cet égard, nous proposons, dans cette étude, des axes d'amélioration de la représentation des glaciers au sein d'un modèle hydrologique par un angle d'investigation géomorphologique. Cette approche a permis d'élaborer un modèle probabiliste permettant de décrire les surfaces englacées au sein d'un bassin versant selon une courbe de niveau. / Mountainous environments are a privileged place of water and energy exchange. Mountainous rivers feed about 40% of the world population and are subjected to climate change and a growing demography. In this context, the comprehension of meteorological, hydrological and hydrogeological processes is essential for a better overall management of water resource. This PhD study is focused squarely on the mountainous environments where hydrology is influenced by snow cover and glaciers, and introduces a multidisciplinary modeling approach in order to improve our comprehension of the process involved.Today, hydrological models are able to simulate gauged mountainous river streamflows under the influence of snow and glaciers but some uncertainties remain when applying such models out of their calibration phase (for instance in response to a different climate or on a ungauged basin). The main uncertainty source is the lack of knowledge of mountainous precipitations, whose measure is sparse and uncertain. It remains the principal component of the hydrological budget though. In this study, we suggest using the meteorological information provided by snow cover and glaciers as “giant pluviometers” to their reservoir scales, with a conceptual reservoir model associated with the concept of watershed.The information provided by hydrology processes, snow and glaciers is assessed in a multi-objective calibration phase. Results show that, in this configuration, the joint validation of the hydrological model by daily streamflow, annual mass balance and daily local snow depth reduces significantly the uncertainty on the meteorological forcing and improves the model robustness. This preliminary result has motivated, consequently, the local annual mass balance of the glaciers.Also, the representation of glaciers in a hydrological model raises a certain amount of issues, especially in the perspective of simulation long-term hydrological processes. In particular, the consideration of the evolution of the glacier geometry is at an early stage. To this end, we propose, in this study, ways of improvements for the representation of glaciers from a geomorphological perspective. This approach allowed us to build a probabilistic model able to describe the glaciated surfaces within a watershed according to a given topographic contour line.