No description available.

Les régions polaires représentent actuellement les régions mondiales où les changements liés au réchauffement climatique se manifestent de manière particulière et souvent spectaculaire. Notre capacité actuelle à comprendre les modes naturels de la variabilité climatique dans ces zones, ainsi que les impacts directs ou indirects de l'activité anthropique sur ces modes, reste encore relativement limitée. Ainsi, la nappe glaciaire Antarctique et ses composantes (atmosphère, océan, glace de mer) restent encore mal représentées dans les modèles climatiques actuels, contribuant notamment à de fortes incertitudes sur les projections climatiques futures et l'évolution associée du niveau des océans. La poursuite des études sur l'équilibre ou non (et l'évolution future) du bilan de masse (surface et total) des différents secteurs du continent Antarctique, par rapport au réchauffement climatique actuel, s'avère donc primordiale. Dans le cadre de ce manuscrit, je propose de contribuer à réduire les incertitudes du Bilan de Masse de Surface (BMS) Antarctique, permettant d'aborder une optimisation de la cartographie actuelle de ce dernier que ce soit avec les modèles climatiques ou les cartes d'interpolation de BMS.

Dans ce travail, nous avons adapté et mis en oeuvre le modèle de neige CROCUS pour le calcul et la reconstruction du bilan de masse des glaciers alpins. Le bilan de masse est calculé au pas horaire à partir des données de température, vitesse du vent, humidité, précipitations, rayonnement infra-rouge, rayonnement solaire direct et diffus et nébulosité. Il est spatialisé sur l'ensemble de la surface du glacier avec un pas horizontal de 200 m. Les données météorologiques pour la validation du modèle proviennent du modèle SAFRAN. La glace de glacier a été introduite dans CROCUS (albédo, rugosité), et la chaîne SAFRAN/CROCUS a été validée à partir de mesures de terrain et de reconstructions photogrammétriques du bilan de masse de surface et d'images satellite de la ligne de neige. Le bilan de masse pour la période 1981-2004 est reconstruit pour les glaciers de St Sorlin et d'Argentière, en bon accord avec les données disponibles.<br /><br />Le modèle permet de tester la sensibilité du bilan de masse des glaciers aux paramètres météorologiques. Pour des amplitudes de variations des paramètres météorologiques compatibles avec un changement climatique réaliste (prévisions climatiques de type IPCC), la température se révèle le paramètre dominant, et la sensibilité de l'altitude de la ligne d'équilibre est de 125 m.K-1, ou même 160 m.K-1 si la variation concomitante du rayonnement infra-rouge est prise en compte. L'impact de la variabilité des paramètres et des combinaisons de paramètres entre eux n'est pas négligeable par rapport à l'impact d'un changement de la moyenne météorologique (tendance).<br /><br />Des simulations longues de bilan de masse ont été effectuées sur la période 1800-2165, en utilisant des résultats de simulations climatiques réalisées par plusieurs modèles de circulation générale dans le cadre de l'IPCC AR4. D'après le modèle, le glacier de St Sorlin devrait disparaître avant 2100.

La zone non saturée (ZNS) des sols contient de l'eau plus ou moins fortement retenue par succion capillaire. Concrètement, cette succion modifie la pression interne de l'eau capillaire et donc modifie ses propriétés thermodynamiques. De ce fait, la spéciation des solutés aqueux, la solubilité des gaz et des solides sont modifiées. Pour permettre de faire des bilans de masse à l'échelle des profils, on a utilisé les courbes de rétention en eau du sol, qui relient un volume d'eau porale à une gamme de succion capillaire. A cette dernière gamme sont associées les propriétés de la succion moyenne. Certains modèles de courbes de rétention sont spécialement adaptés au domaine des fortes succions et ils sont ici étudiés en détail. Leur utilisation permet de constater que les volumes retenus à des succions géochimiquement significatives peuvent atteindre plusieurs dizaines de litres par mètre cube de sol et ainsi participer au transfert de masse en profondeur. Plusieurs faits d'observations pourraient être re-interprétés par cette approche capillaire, comme les successions anormales de paragénèses dans certains régolites.<br />Des premières études de bilans de masse ont été réalisées en s'appuyant sur un modèle en double porosité, avec un domaine macroporal de percolation, et un domaine microporal de rétention. Ces deux zones sont couplées par des échanges diffusifs. Le modèle utilisé est celui des colonnes de transfert-réactif 1Dde PhreeqC, un des modèles géochimiques les plus courants. Les caractéristiques géochimiques des solutions capillaires sont calculées avec Thermo-ZNS, un logiciel de calculs thermodynamiques développé pour la ZNS au BRGM. Les résultats de simulation des colonnes mettent en évidence l'effet de la capillarité sur la stabilité des minéraux, la vitesse des réactions et la complexation de surface, consistant avec les observations naturelles de terrain.

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 / Albedo is defined as the ratio of reflected to incident radiation over the solar spectrum and is a key parameter in the surface energy balance of snow and ice. This parameter is highly variable both temporally and spatially; thus remote sensing is an ideally suited approach for the retrieval of albedo data.The albedo value depends on both physical properties of the target and the characteristics of the incident radiation. Furthermore, the physical parameters linked with the albedo concept vary in consideration with spectral and angular ranges. Measurements of bi-directional reflectance over natural snow have shown that the Lambertian hypothesis may lead to significant error when estimating albedo from remote sensing data.Detailed knowledge of the angular distribution of radiation reflected by snow allows for the development of a new method to retrieve albedo values for mountainous, snow/ice covered areas. This method takes into account multiple reflections on mountainous areas, anisotropy of radiation reflected by snow and ice, and albedo spectral variations with surface physical properties. The method is applied to visible and near-infrared terrestrial photographs (spatial resolution 10 m) and MODIS data (spatial resolution 250 m). The accuracy of the method is evaluated at ±10% on the retrieved albedo value using concurrent field measurements at theSaint Sorlin Glacier (Grandes Rousses, France) during the summers of 2008 and 2009. The method is used to retrieve albedo data for this glacier from 2000 to 2009. Results indicate that the albedo of the ablation area of the Saint Sorlin Glacier has not shown any decreasing trend over this decade, in opposition to results presented for the Morteratsch Glacier (Switzerland). In addition, the minimal value over the summer period of the whole glacier averaged albedo is highly correlated to the specific annual mass-balance.Albedo data from MODIS and terrestrial photographs are then assimilated into the snow model CROCUS. This assimilation allows for an estimation of the spatialized mass-balance of the Saint Sorlin Glacier over the five studied hydrological years. Root mean square error is evaluated to 0.5 m w.e. For this study, we have used mid-scale meteorological data from SAFRAN. A brief analysis of the contribution from the atmospheric fluxes to the surface energy balance shows that, for the time period considered in this study, the shortwave radiation budget is the main process determining the surface energy balance. Furthermore, variability in shortwave radiation budget explains the major part of the daily variability in surface energy balance.The methods developed in this work are readily applicable to other temperate glaciers. They allows spatialized mass-balance reconstruction on a decadal scale and lead to improved quantification of the physical processes controlling mass-balance in temperate glaciers

Sur la Cordillère Royale en Bolivie, nous avons reconstitué les bilans de masse de 21 glaciers représentatifs, par photogrammétrie numérique. Les MNT's ont été construit à partir des clichés aériens de l'IGM-Bolivie de 1956, 1963, 1975, 1983, 1997 et 2006, avec une précision décimétrique ou métrique.<br />L'analyse des résultats sur le glacier du Zongo montre qu'un tiers de la surface du glacier (zone d'ablation) contrôle 80% du bilan de masse annuel. La série du Zongo constitue la série de bilan de masse glaciaire la plus longue et la plus précise des Andes tropicales.<br />Nous avons analysé les bilans de masse des 21 glaciers depuis 1956, lorsqu'on s'affranchit des tendances de ces glaciers, les variations temporelles du bilan de masse montrent un signal commun d'origine climatique.<br />Nous avons montré ainsi que 88% de la variance des bilans cumulés sont expliqués par l'exposition et l'altitude. Cette corrélation et les restitutions photogrammétriques ont permis de calculer les bilans de masse et la perte de surfaces de 376 glaciers au cours de la période 1963-2006 (43% de perte de leur volume et 48% de perte de leur surface depuis 1975).<br />Enfin, notre étude évalue la contribution des glaciers au système d'alimentation en eau potable de la ville de La Paz. En moyenne sur la période 1975-2006, la fonte glaciaire a contribué à hauteur de 18% et 15% sur la période récente 1997-2006 (12% en saison humide et 27% en saison sèche). Nous avons enfin calculé les débits attendus en cas de disparition totale des glaciers (12% en moins à l'année, 9% en moins en saison humide et 25% en moins en saison sèche).

La qualité globale des eaux continentales (lacs et réservoirs) continue à se détériorer dans de nombreuses régions du Mexique. Le réservoir Cointzio, situé dans la partie sud du Plateau Central Mexicain, sur la ceinture "Trans-Mexicain Volcanique Belt " (TMVB), ne fait pas exception. Ce réservoir tropical turbide est un système monomictique chaud (surface = 6 km2, capacité = 66 Mm3, temps de séjour ~ 1 an). Il est stratégique pour l'approvisionnement en eau potable de la ville de Morelia, capitale de l'état du Michoacán, et pour l'irrigation en aval pendant la saison sèche. Le réservoir est menacé par l'accumulation des sédiments et des nutriments provenant d'eaux non traitées dans le bassin versant amont. La haute teneur en particules d'argile très fines et le manque d'installations de traitement de l'eau conduisent à de graves épisodes d'eutrophisation (jusqu'à 70 µg chl. a L-1) ainsi qu'à des périodes saisonnières de forte turbidité (profondeur de Secchi < 30 cm) et d'anoxie (de Mai à Octobre). Sur la base de mesures intensives de terrain réalisées en 2009 (échantillonnage dans le bassin versant, sédiments déposés, profils verticaux de l'eau, bilan entrées et sorties) nous avons présenté une étude intégrée du fonctionnement hydrodynamique et biogéochimique du réservoir Cointzio. Les profils verticaux de température (T), de concentration en sédiment en suspension (TSS), d'oxygène dissous (DO), de chlorophylle a, de carbone et de nutriments ont été utilisés pour évaluer les cycles internes au réservoir. Pour compléter l'analyse des données de terrain, nous avons examiné la capacité de différents modèles numériques 1DV (Aquasim modèle biogéochimique couplé avec le modèle de mélange k-ε) à (i) reproduire les principaux cycles biogéochimiques dans le réservoir Cointzio et (ii) préciser les scénarios de réduction des éléments nutritifs (P and N) pour limiter l'eutrophisation dans les prochaines décennies. Le modèle k-ε s'est avéré le plus robuste pour reproduire les conditions hydrodynamiques, en particulier la stratification modérée liée à la très grande turbidité et au régime soutenu de vent thermique. Le modèle Aquasim s'est avéré pertinent pour reproduire les principaux cycles de l'oxygène dissous DO, des nutriments et de la chlorophylle a au cours de l'année 2009. Les différentes simulations ont souligné l'impact négatif à long terme du réchauffement climatique. A la fin du siècle (2090), une augmentation de la température de l'air atteignant 4.4° C a été prédite à partir de modèles de circulation mondiale. Couplé avec une année hydrologique sèche, une telle situation pourrait conduire à des conditions anoxiques sévères et à des blooms importants de chlorophylle a (jusqu'à de 94 µg L-1). Diverses simulations ont montré que la réduction drastique de l'apport de nutriments (à 90%) serait nécessaire pour réduire de façon significative les teneurs en chlorophylle a. Si de telles mesures d'atténuation sont adoptées, le pic maximal de chlorophylle a se stabiliserait à 40 µg L-1, au lieu de 94 mg L-1, après une période de cinq ans d'efforts. À notre connaissance, cette étude fournit la première application numérique de modèles k-ε et AQUASIM pour simuler les niveaux d'eutrophisation élevés dans un réservoir tropical très turbide.

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

Les glaciers intègrent naturellement des informations hydrologiques et énergétiques dans des régions climatiques variées du globe et sont sensibles à des changements minimes de bilan d’énergie de surface. Parmi les paramètres mesurés in-situ traduisant les variations climatiques, le bilan de masse de surface est la variable directement reliée aux conditions atmosphériques qui contrôlent les processus d’accumulation et d’ablation. Un grand nombre de modèles, de complexités diverses, existent pour simuler l’évolution du bilan de masse. Cependant, des questions perdurent quant à l’utilisation de l’approche la plus appropriée, notamment sur de longues périodes de temps. Basés sur de longues séries de mesures réalisées sur quatre glaciers des Alpes françaises (issues du service d’observation GLACIOCLIM), ces travaux de thèse ont pour but d’étudier la sensibilité du bilan de masse aux variables atmosphériques et topographiques, afin d’identifier celles essentielles pour contraindre les modèles.En s’appuyant tout d’abord sur une approche empirique, nos résultats indiquent que pour chaque glacier, la variabilité temporelle de l’ablation s’explique principalement par la variabilité de la température tandis que la variabilité spatiale de la fonte de la neige et de la glace est largement gouvernée par l’insolation. Néanmoins, à l’échelle du glacier, considérer l’insolation en plus de la température ne conduit pas à une amélioration de la performance des modèles empiriques et un modèle type degré-jour apparaît suffisant pour simuler son bilan de masse. L’utilisation du modèle à base physique SURFEX/ISBA-Crocus a permis d’étudier la sensibilité du bilan de masse estival à chacune des variables du bilan d’énergie de surface. Le modèle s’avère très performant pour la simulation des bilans spatialisés, mais nécessite un forçage météorologique précis et distribué spatialement. En particulier, le vent, souvent complexe à représenter dans les modèles atmosphériques,apparaît comme un facteur prépondérant pour modéliser correctement les processus de fonte. Avec des forçages météorologiques bien contraints, ces deux approches (empirique et physique) se révèlent très performantes quant à la modélisation du bilan de masse estival. Toutefois, nos résultats montrent que les bilans annuels restent très conditionnés par les bilans hivernaux du fait de la rétroaction de l’albédo de surface. Une attention particulière a donc été portée sur la simulation de l’accumulation.A partir des données de l’observatoire GLACIOCLIM, notre étude n’a pas permis d’établir de relation significative entre la variabilité spatio-temporelle des bilans hivernaux et les variables morpho-topographiques.Pour cette raison, une campagne spécifique basée sur des mesures LiDAR a été mise en œuvre au Col du Midi (Massif du Mont-Blanc) afin d’affiner le réseau de mesures et d’améliorer la compréhension de la variabilité spatio-temporelle de la couverture neigeuse. Nos résultats montrent la complexité d’obtenir des mesures d’accumulation par différence de MNT, du fait des vitesses d’émergence/submergence. / Glaciers naturally record hydrological and energetic information into varied climatic regions of the world, and are sensitive to small changes in their surface energy balance. The glacier surface mass balance is directly linked to the atmospheric variables which control accumulation and ablation processes. Many models, with different complexities, enable simulations of mass balance evolution, but some questions arise about the best approach to use, especially over long time periods. Based on long measurement series on four glaciers in the French Alps (GLACIOCLIM observatory), this thesis aims at studying the mass balance sensitivity to climatic variations, to identify relevant meteorological variables for modelling. Based on an empirical approach performed on each glacier, temperature was found to be the main driver of temporal ice/snow ablation variability while solar radiations strongly influence the spatial distribution of summer mass balance. However, at glacier scale, to include solar radiation in melt models does not improvethe performances and a classical degree-day model is sufficient to simulate glacier-wide mass balance.The snowpack model SURFEX-ISBA-Crocus was then used to study summer mass balance sensitivity to all surface energy balance fluxes. Results indicate good performances of the model to simulate mass balance at each stake, but it requires accurate meteorological forcing. In particular wind, generally difficult to assessby atmospheric models, appears as a key factor for an accurate ablation modelling. Our work emphasizes that both physical and empirical approaches are very efficient for ablation modeling when forced with accurate meteorological data. Yet, annual mass balances remain very sensitive to wintermass balance due to the surface albedo feedback. For this reason, simulating accumulation processes remainsas important as ablation modelling. Using GLACIOCLIM data measurements, our study failed to find relevant relationship between winter mass balance and topographical variables such as slope, aspect, elevation, due to insufficient spatio-temporal resolution. To fill this gap of measurements, we set up a dedicated field campaign to improve our understanding of the spatio-temporal variability of the snow depth distribution on Alpineglaciers. This year-round campaign was performed using LiDAR acquisitions at Col du Midi (Mont Blancmassif). Results underline the complexity of accurate accumulation measurements from DEM differences due to submergence/emergence velocities.

La qualité globale des eaux continentales (lacs et réservoirs) continue à se détériorer dans de nombreuses régions du Mexique. Le réservoir Cointzio, situé dans la partie sud du Plateau Central Mexicain, sur la ceinture “Trans-Mexicain Volcanique Belt » (TMVB), ne fait pas exception. Ce réservoir tropical turbide est un système monomictique chaud (surface = 6 km2, capacité = 66 Mm3, temps de séjour ~ 1 an). Il est stratégique pour l'approvisionnement en eau potable de la ville de Morelia, capitale de l'état du Michoacán, et pour l'irrigation en aval pendant la saison sèche. Le réservoir est menacé par l'accumulation des sédiments et des nutriments provenant d'eaux non traitées dans le bassin versant amont. La haute teneur en particules d'argile très fines et le manque d'installations de traitement de l'eau conduisent à de graves épisodes d'eutrophisation (jusqu'à 70 µg chl. a L-1) ainsi qu'à des périodes saisonnières de forte turbidité (profondeur de Secchi < 30 cm) et d'anoxie (de Mai à Octobre). Sur la base de mesures intensives de terrain réalisées en 2009 (échantillonnage dans le bassin versant, sédiments déposés, profils verticaux de l'eau, bilan entrées et sorties) nous avons présenté une étude intégrée du fonctionnement hydrodynamique et biogéochimique du réservoir Cointzio. Les profils verticaux de température (T), de concentration en sédiment en suspension (TSS), d'oxygène dissous (DO), de chlorophylle a, de carbone et de nutriments ont été utilisés pour évaluer les cycles internes au réservoir. Pour compléter l'analyse des données de terrain, nous avons examiné la capacité de différents modèles numériques 1DV (Aquasim modèle biogéochimique couplé avec le modèle de mélange k-ε) à (i) reproduire les principaux cycles biogéochimiques dans le réservoir Cointzio et (ii) préciser les scénarios de réduction des éléments nutritifs (P and N) pour limiter l'eutrophisation dans les prochaines décennies. Le modèle k-ε s'est avéré le plus robuste pour reproduire les conditions hydrodynamiques, en particulier la stratification modérée liée à la très grande turbidité et au régime soutenu de vent thermique. Le modèle Aquasim s'est avéré pertinent pour reproduire les principaux cycles de l'oxygène dissous DO, des nutriments et de la chlorophylle a au cours de l'année 2009. Les différentes simulations ont souligné l'impact négatif à long terme du réchauffement climatique. A la fin du siècle (2090), une augmentation de la température de l'air atteignant 4.4° C a été prédite à partir de modèles de circulation mondiale. Couplé avec une année hydrologique sèche, une telle situation pourrait conduire à des conditions anoxiques sévères et à des blooms importants de chlorophylle a (jusqu'à de 94 µg L-1). Diverses simulations ont montré que la réduction drastique de l'apport de nutriments (à 90%) serait nécessaire pour réduire de façon significative les teneurs en chlorophylle a. Si de telles mesures d'atténuation sont adoptées, le pic maximal de chlorophylle a se stabiliserait à 40 µg L-1, au lieu de 94 mg L-1, après une période de cinq ans d'efforts. À notre connaissance, cette étude fournit la première application numérique de modèles k-ε et AQUASIM pour simuler les niveaux d'eutrophisation élevés dans un réservoir tropical très turbide. / La qualité globale des eaux continentales (lacs et réservoirs) continue à se détériorer dans de nombreuses régions du Mexique. Le réservoir Cointzio, situé dans la partie sud du Plateau Central Mexicain, sur la ceinture “Trans-Mexicain Volcanique Belt » (TMVB), ne fait pas exception. Ce réservoir tropical turbide est un système monomictique chaud (surface = 6 km2, capacité = 66 Mm3, temps de séjour ~ 1 an). Il est stratégique pour l'approvisionnement en eau potable de la ville de Morelia, capitale de l'état du Michoacán, et pour l'irrigation en aval pendant la saison sèche. Le réservoir est menacé par l'accumulation des sédiments et des nutriments provenant d'eaux non traitées dans le bassin versant amont. La haute teneur en particules d'argile très fines et le manque d'installations de traitement de l'eau conduisent à de graves épisodes d'eutrophisation (jusqu'à 70 µg chl. a L-1) ainsi qu'à des périodes saisonnières de forte turbidité (profondeur de Secchi < 30 cm) et d'anoxie (de Mai à Octobre). Sur la base de mesures intensives de terrain réalisées en 2009 (échantillonnage dans le bassin versant, sédiments déposés, profils verticaux de l'eau, bilan entrées et sorties) nous avons présenté une étude intégrée du fonctionnement hydrodynamique et biogéochimique du réservoir Cointzio. Les profils verticaux de température (T), de concentration en sédiment en suspension (TSS), d'oxygène dissous (DO), de chlorophylle a, de carbone et de nutriments ont été utilisés pour évaluer les cycles internes au réservoir. Pour compléter l'analyse des données de terrain, nous avons examiné la capacité de différents modèles numériques 1DV (Aquasim modèle biogéochimique couplé avec le modèle de mélange k-ε) à (i) reproduire les principaux cycles biogéochimiques dans le réservoir Cointzio et (ii) préciser les scénarios de réduction des éléments nutritifs (P and N) pour limiter l'eutrophisation dans les prochaines décennies. Le modèle k-ε s'est avéré le plus robuste pour reproduire les conditions hydrodynamiques, en particulier la stratification modérée liée à la très grande turbidité et au régime soutenu de vent thermique. Le modèle Aquasim s'est avéré pertinent pour reproduire les principaux cycles de l'oxygène dissous DO, des nutriments et de la chlorophylle a au cours de l'année 2009. Les différentes simulations ont souligné l'impact négatif à long terme du réchauffement climatique. A la fin du siècle (2090), une augmentation de la température de l'air atteignant 4.4° C a été prédite à partir de modèles de circulation mondiale. Couplé avec une année hydrologique sèche, une telle situation pourrait conduire à des conditions anoxiques sévères et à des blooms importants de chlorophylle a (jusqu'à de 94 µg L-1). Diverses simulations ont montré que la réduction drastique de l'apport de nutriments (à 90%) serait nécessaire pour réduire de façon significative les teneurs en chlorophylle a. Si de telles mesures d'atténuation sont adoptées, le pic maximal de chlorophylle a se stabiliserait à 40 µg L-1, au lieu de 94 mg L-1, après une période de cinq ans d'efforts. À notre connaissance, cette étude fournit la première application numérique de modèles k-ε et AQUASIM pour simuler les niveaux d'eutrophisation élevés dans un réservoir tropical très turbide.