Analyse de la fonte glaciaire et nivale dans les Andes tropicales à partir d’un bilan d’énergie : Glacier de l’Antisana, Equateur (0°28’S)

Antisana Glacier 15 basin surface mass balance (SMB) (0.28 km2;0°28'S,78° 09'W, 5750 m), located in the tropical Andes of Ecuador, between 2000 and 2008 has been modeled using two models: an empirical positive degree-day (PDD) and other physical surface energy balance approach..As a first stage, melting was calculated using daily temperature and the cumulative precipitation while sublimation was computed as wind-speed function. The PDD was applied at daily-time step and calibrated at 4900 m a.s.l. The period time to calibrate was March 2002 - August 2003 and validated was January - November 2005. A significant link between melting rates and temperature was revealed when model discriminates snow and ice surface. This relationship is explained by the link between the short-wave radiation and air temperature. However, this relationship disappears under intense wind-speed events, however this events show low air temperature, therefore they are not affect the final results. The PDD model was applied for modeling SMB, ablation, snow-line position and equilibrium line altitude. The results show this model type fits to SMB measures over Antisana glacier. However, it should not apply to other tropical regions, particularly where there are high sublimation rates (pronounced dry seasons) or where glaciers are located under the isotherm 0 ° C.A second stage, a distributed surface energy balance (DBES) model over glacier 15a of the Antisana was applied to quantify ice and snow ablation. For this task was used a model adapted to tropical conditions from Adélie Land-Antarctica BES model. Initially, the model was applied at 4900 m.a.s.l. , with a slope of 28º and NW (45 °) azimuth. For computing DBES, the glacier was divided in 18 altitude strips ( 50m/level). The model generated slops and azimuths from a digital elevation model (DEM). The meteorological vertical gradients were measured in the study zone. This values was used into de model as constants during both periods calibration and validation. The DBES model reproduced the spatial variation and temporal of SMB between 4850 m.a.s.l. and 5700 m.a.s.l.. The unpublished results of SEB distributed on the glacier 15a allowed to analyze the physical processes that to generate the spatio-temporal variation of the ablation.The DBES results show that the model applied in a former research was based on wrong assumptions and inaccurate. The results of DBES and PDD models show that the accumulated precipitation measured at the bottom of the glacier are greatly underestimated. It implies the need for new protocols for improving the precipitation measurement. Both the former research and SEB confirm the crucial role of shortwave radiation and albedo values over glacier melt. It puts also in evidence that snow metamorphism in the study zone occurs faster than other latitudes, associated with accelerated degradation of surface albedo. Finally, this research explains the reason why the PDD model provides good quality results in this region, which contradicted the results of previous study. This research contributes a better understanding about sensibility of SMB with temperature variations in the inner tropics.Ce travail a permis de modéliser le bilan de masse superficiel sur le bassin versant du Glacier 15 du massif Antisana (0.28 km2; 0°28'S, 78°09'W), dans les Andes tropical de l'Equateur. Le bilan de masse superficiel (MBS) entre 2000 et 2008 sur le glacier 15a de l'Antisana a été modélisé en utilisant deux approches parallèles : un modèle empirique de type degré-jour et un modèle physique de type bilan d'énergie.Dans un premier temps, la fonte a été calculée à l'aide de la température moyenne et des cumuls de précipitations alors que le calcul de la sublimation a été réalisé grâce à la vitesse du vent. Ce modèle a été appliqué au pas de temps journalier et a été calé à 4900 m d'altitude entre mars 2002 et août 2003 et validé entre janvier et novembre 2005. Un lien significatif entre la température et la fonte apparaît lorsqu'une distinction entre neige et glace est effectuée. La relation entre la fusion et la température est expliquée essentiellement par le lien significatif existant entre le bilan radiatif de courtes longueurs d'onde et la température, car le rayonnement solaire est aussi le principal moteur de la fonte. Néanmoins, cette relation disparaît dès lors que le vent devient intense. Ce point a cependant peu de conséquences sur les calculs effectués, car les périodes ventées sont marquées par des températures faibles. Le modèle a été appliqué pour évaluer dans quelle mesure cette approche est appropriée pour modéliser le bilan de masse superficiel, l'ablation, et les altitudes de la ligne de neige et de la ligne d'équilibre. Les résultats montrent qu'un modèle empirique de type degré-jour permet une modélisation précise du MBS à l'échelle d'un glacier en Equateur. Néanmoins, ce modèle ne devrait pas être appliqué à d'autres régions tropicales, en particulier là où la sublimation est importante (grâce à une saison sèche prononcée) ou lorsque les glaciers sont situés au-dessus de l'altitude moyenne de l'isotherme 0°C.Dans un deuxième temps, un modèle de bilan d'énergie superficiel (BES) distribué développé initialement pour des applications en Terre d'Adélie (Antarctique) a été appliqué sur le site du glacier 15a de l'Antisana pour quantifier l'ablation de la glace et de la neige. Ce modèle a été adapté pour les conditions tropicales puis forcé à l'aide des mesures météorologiques obtenues à proximité et sur glacier 15a. Le modèle a tout d'abord été appliqué à une altitude de référence à 4900 m pour une pente de 28° et avec un azimut NO (45°). La spatialisation du bilan d'énergie a ensuite été effectuée en discrétisant le glacier par tranches de 50 m d'altitude, puis en calculant la pente et l'azimut moyen de chaque tranche à l'aide de un modèle numérique de terrain (MNT) précis du glacier. Les données d'entrées ont été distribuées sur la langue du glacier 15a en supposant que les forçages météorologiques varient selon des gradients verticaux constants. Ces gradients ont été établis à partir de mesures de terrain. Le modèle a été calé à partir des données de bilan de masse mensuel et d'albédo journalier mesurées sur le glacier 15a entre mars 2002 et août 2003, puis validé à partir de mesures obtenues entre janvier 2005 et décembre 2008. Le modèle reproduit précisément les variations spatiales et temporelles du bilan de masse entre 4850 m et 5700 m d'altitude. Les résultats inédits de bilan d'énergie de surface distribués sur le glacier 15a ont ainsi permis d'analyser les processus physiques à l'origine des variations spatio-temporelles de l'ablation. Enfin, cette étude explique pour la première fois pourquoi les modèles de type degré-jour donnent systématiquement des résultats de bonne qualité dans cette région, venant contredire les résultats d'études antérieures. Cette thèse constitue ainsi une étape importante dans la compréhension de la sensibilité des bilans de masse aux variations de température dans cette région