Le bassin versant de la Durance assure 10% de la production hydroélectrique française, et alimente en eau potable de nombreuses villes du pourtour méditerranéen. Sa partie alpine contribue à près de la moitié des débits à sa confluence avec le Rhône. Connaître les impacts du changement climatique sur l'hydrologie de la Durance présente donc un enjeu de taille pour l'avenir de la population provençale. A cette fin, les modèles hydrologiques sont devenus des outils indispensables. Cependant, le manteau neigeux, de par sa forte variabilité spatiale, est particulièrement difficile à modéliser. Une attention particulière a été portée à la compréhension des processus nivaux, dominants dans la partie amont du bassin, et à leur représentation dans le modèle de surface continentale CLSM. L'analyse de mesures d'extension et de hauteur du manteau neigeux nous a amenés à modifier le schéma de neige de CLSM, en introduisant une hystérésis dans la relation entre ces deux variables. Nous avons ensuite évalué le potentiel devenir sous changement climatique de l'hydrosystème. Les résultats des simulations de CLSM, réalisées à l'aide de 330 scénarios climatiques, ont été confrontés à ceux de cinq autres modèles hydrologiques. Tous s'accordent sur une baisse significative des débits annuels, avec des modifications du régime hydrologique différentes selon les zones. L'incertitude reste importante sur la magnitude des changements de débits, principalement due aux scénarios climatiques. La part d'incertitude associée aux modèles hydrologiques varie selon les saisons, mettant en évidence une plus grande difficulté à représenter les processus nivaux et ceux liés à l'évapotranspiration.

No description available.

Cette thèse s'inscrit dans le développement de la plateforme EauDyssée de modélisation intégrée des hydrosystèmes régionaux, au sein du bassin pilote de la Seine. L'objectif principal est de contribuer à une meilleure simulation des tirants d'eau à l'échelle régionale afin d'améliorer la simulation des interactions nappe-rivière et de mieux quantifier les niveaux piézométriques dans les aquifères. La première partie de la thèse vise à évaluer la sensibilité d'un modèle hydraulique à la précision de la description géomorphologique des lits pour identifier le meilleur compromis entre parcimonie et réalisme et identifier les facteurs morphologiques les plus importants pour obtenir une simulation satisfaisante des tirants d'eau à l'échelle régionale. Cette étude est menée sur le Serein (affluent de l'Yonne), entre les stations limnimétriques de Dissangis et Beaumont, dans un bief bien renseigné (20 sections transversales sur 89 kms). Débits et tirants d'eau sont simulés par le modèle hydraulique HEC-RAS (équations de Saint-Venant 1D), en fonction des apports latéraux simulés par le modèle régional EauDyssée. Les résultats de cette étude montrent qu'un modèle 1D type Saint-Venant n'est pas adapté à la simulation des écoulements à l'échelle régionale. Nous avons donc développé une méthode de changement d'échelle originale, dans laquelle la modélisation fine des processus hydrauliques à haute résolution permet d'améliorer la représentation des profils d'eau en rivière et les interactions nappe-rivière simulées à l'échelle régionale par le modèle intégré EauDyssée. Cette méthodologie de changement d'échelle a été validée dans un sous bassin versant de l'Oise d'une superficie de 4500 km2, pour la période 1990-1995. Nous avons utilisé HEC-RAS pour la modélisation hydraulique d'un tronçon de l'Oise de 188 km, où 420 sections transversales sont disponibles. Le modèle permet d'interpoler des courbes de tarage simulées tous les 200m en moyenne. Ces courbes de tarage sont ensuite projetées sur les mailles rivière du modèle régional EauDyssée (résolution de 1 km), où elles permettent de simuler la fluctuation du niveau d'eau en fonction du débit à l'échelle régionale par EauDyssée. La cote de la surface libre de la rivière définissant sa charge hydraulique, ces fluctuations influencent alors les échanges entre les mailles rivière et les nappes, qui dépendent des gradients de charge verticaux entre rivière et nappe (loi de Darcy). Ce travail montre l'intérêt de l'approche pour mieux évaluer les interactions nappes-rivières à l'échelle régionale avec un faible coût de calcul. Il offre des perspectives intéressantes pour simuler des processus jusque là négligés par le modèle EauDyssée : élimination de nitrate dans les zones humides qui sont souvent situées à la zone de contact entre les nappes souterraines et la rivière, ou l'impact du changement climatique sur le fonctionnement des hydrosystèmes et plus particulièrement sur l'élimination ou le relarguage de polluants par des processus biogéochimiques, ainsi que de mieux estimer les risques d'inondation à l'échelle régionale.

L'étude hydrogéologique du bassin versant de la Bourbre a été complétée par les nécessités de création de la ville nouvelle de l'Isle d'Abeau. La gestion et la protection des ressources en eau ont été particulieremnt étudiées.

L'étude géologique et hydrogéologique de la partie centrale du massif des Bauges est traitée dans ce travail notamment la quaternaire

Notre étude, qui a débuté en juillet 1977, a été tout d'abord principalement dirigée vers le bassin versant du ruisseau du Ternèze et du massif de la Thuile à l'Est de Chambéry. Il s'agissait de préciser les relations hydrogéologiques pouvant exister entre cette petite région et les vallées de l'Isère - Combe de Savoie la limitant au Sud. L'équipement hydrologique et météorologique a été installé dans ce secteur en fonction de ces recherches à effectuer. De cette façon, une liaison était établie entre les différents travaux hydrogéologiques effectués dans le massif par G. Nicoud (Haute Vallée du Chéran) , Y. Lemordant (Plateau du Mont Revard) , M. Lepiller à l'Est (Semnoz) , Benedetti Crouzet (Bassin versant à l'Ouest du Lac d'Annecy au Nord) , et dans la plaine de l'Isère au Sud par J.C. Fourneaux. Mais nous avons rencontré quelques difficultés dans l'installation de l'équipement hydrologique de cette vallée des Aillons,qui n'a pu être réalisée Au cours de cette étude, nous allons examiner successivement la géologie de la région, son hydrogéologie ainsi qu'une étude hydrochimique se rapportant aux principales émergences.

Cette thèse étudie une large gamme de processus liés au cycle de l'eau. La sensibilité du cycle hydrologique global à la représentation de l'hydrologie du sol est analysée à l'aide d'un Modèle de Circulation Générale (MCG). Trois directions sont explorées : l'influence d'un terme de drainage, celle de la capacité en eau du sol, et celle de la variabilité sous-maille de l'hydrologie continentale, à travers une paramétrisation distribuée du ruissellement. La sensibilité de la branche atmosphérique du cycle hydrologique dépend de l'amplitude et du signe des variations d'évaporation continentale moyennes induites par les paramétrisations étudiées. Elle dépend aussi de la distribution spatiale et temporelle de ces variations. Cette distribution est particulièrement sensible à l'hydrologie du sol dans les tropiques, où elle entraine des modifications importantes de la circulation de Hadley-Walker et de la mousson indienne associée. Les paramétrisations de l'hydrologie du sol n'ont cependant qu'une influence très faible sur la convergence d'humidité moyenne sur les continents, en regard de la surestimation systématique de cette grandeur dans le MCG. La seconde partie de cette thèse présente le développement d'un Modèle de Transport Latéral de l'eau dans les grands bassins versants, le MTL, et son application pour simuler le débit de 14 grands fleuves, à partir de l'écoulement total journalier simulé par un MCG. Le MTL permet de reproduire de manière satisfaisante le retard entre le maximum de l'écoulement total sur un bassin et le maximum du débit à l'aval du bassin, si l'on tient compte des différences entre les précipitations simulées et observées, et des simplifications de l'hydrologie du sol dans le MCG.

En 2004 et 2007, neuf structures d'échappements de fluides (SEF), constituées de volcans de boue (VdB), d'évents de gaz (GS, gas seeps en anglais) et/ou de sources, ont été échantillonnées dans la partie émergée du Bassin de la Côte Est (BCE) de l'île nord de la Nouvelle-Zélande.L'étude granulométrique indique que la boue émise par les VdB, les roches encaissantes et les niveaux de décollement voisins sont composés d'argiles et de silts en majorité. L'étude de la fraction argileuse et de la roche totale par diffraction des rayons X (DRX) a révélé de nombreuses similitudes. Les volcans de boue, les roches encaissantes et les niveaux de décollement présentent ainsi le même assemblage minéralogique, à savoir smectite, illite, chlorite, kaolinite, quartz et feldspaths. Les proportions sont variables d'un échantillon à un autre mais le couple smectite-illite est toujours majoritaire.L'étude géochimique de la fraction solide indique que les échantillons sont riches en Si02, pauvres en Fe2O3, MgO, MnO et en alcalins, à quelques exceptions près. La composition en éléments majeurs s'organise entre un pôle argileux alcalin et un pôle carbonaté. Les spectres de terres rares sont similaires et caractéristiques des argiles ; Ils présentent également un faible degré de fractionnement, lié à la formation des carbonates. Ce sont donc les mêmes minéraux qui contrôlent la chimie des échantillons.L'étude géochimique de la phase liquide montre que l'eau impliquée dans les volcans de boue est d'origine marine essentiellement, et des réactions eau/roche similaires, notamment l'altération de smectite en illite. Cette étude a permis aussi d'obtenir une estimation de la température d'équilibre, comprise entre 60 et 110°C, impliquant une profondeur d'origine de 2 à 3 km, voire plus.L'étude géophysique indique qu'à l'aplomb des VdB et des deux sources chaudes étudiés, la croûte continentale a sensiblement la même épaisseur et que la profondeur de la croûte océanique en subduction avoisine les 20 km. A cette profondeur, la fusion de la péridotite n'est pas possible et la fusion résultante de la croûte continentale, responsable du volcanisme d'arc, non plus. Le gradient géothermique mesuré à TePuia est donc influencé par un autre phénomène, mais la modélisation de la lithosphère ne nous a pas permis de trouver lequel.Ces différentes études mettent en évidence des caractéristiques géochimiques, pétrographiques et minéralogiques communes aux volcans de boue de Nouvelle-Zélande. Les fluides impliqués dans ces structures proviendraient donc d'un même niveau source, recouvert du même assemblage sédimentaire. L'étude géophysique ne nous apporte aucune information à ce sujet mais permet cependant d'établir avec certitude que le régime thermique est le même du Nord de Hawke's Bay au Sud du BCE ; la région de TePuia est un cas particulier, peut-être influencée par le complexe volcanique de Matakaoa. / In 2004 and 2007, nine areas of fluids escapes structures (FES) with mud volcanoes (MVs), gas seeps (GSs) and springs, were sampled in the East Coast Basin (ECB) emerged part of New Zealand's north island.The grain size study indicates that the mud emitted by MVs, like surrounding rocks and decollement layers in the vicinity, is mainly composed of clays and silts. Clay fraction and whole rock XRD study revealed many similarities. Mud volcanoes, surrounding rocks and decollement layers have exactly the same mineral assemblage, i.e. smectite/illite, chlorite, kaolinite, quartz and feldspars. The proportions vary from one sample to another but the couple smectite/illite is always majority.The geochemical study of the solid fraction indicates that the samples are rich in Si02, low in Fe2O3, MgO, MnO and alkaline, with a few exceptions. The major element composition is organized between a pole and a pole clay alkaline carbonate. The REE patterns are similar and characteristics of clays They also have a low degree of fractionation, due to the formation of carbonates. So these are the same minerals that control the chemistry of samples.The liquid phase geochemical study shows that water involved in the MVs is mainly of marine origin, that water/rock reactions are similar, including weathering of illite in smectite, allowed us to obtain an estimate of the water equilibrium temperature between 60 and 110°C, implying an origin depth of 2-3 km.The geophysical study indicates that directly above the MVs and the two hot springs studied, the continental crust has substantially the same thickness and the depth of the subducting oceanic crust is around 20 km. At this depth, the melting of the peridotite is not possible and the resulting fusion of continental crust, responsible of arc volcanism, neither. The geothermal gradient measured at TePuia is influenced by another phenomenon, but the modeling of the lithosphere does not allow us to find which one.These studies show geochemical, petrographical and mineralogical characteristics common to mud volcanoes in New Zealand. Fluids involved in these structures therefore come from the same parent bed, covered with the same sedimentary package. The geophysical study gives us no information about it but nevertheless allows us to establish with certainty that the thermal regime is the same from Northern Hawke's Bay to the south of the ECB; TePuia region is a special case, perhaps influenced by the Matakaoa volcanic complex.

Les multifractals universels (UM) sont un outil puissant et abondement utilisé d'analyse et de simulation de champs géophysiques, comme la pluie, extrêmement variables sur une large gamme d'échelle. Ils sont basés sur le concept de cascade multiplicative qui repose sur la notion physique d'invariance d'échelle pour explorer le phénomène fondamental qu'est l'intermittence. Dans ce cadre, toute la variabilité du champ est caractérisée à l'aide de simplement trois paramètres qui ont en plus une interprétation physique. Dans cette thèse on utilise ce cadre théorique pour quantifier l'impact de la variabilité à petite échelle de la pluie en hydrologie urbaine. La première étape consiste à analyser la variabilité spatio-temporelle de données radar de précipitation à l'aide d'un modèle multifractal anisotrope simple. Divers évènements pluvieux sont analysés. Un comportement scalant a été observé sur deux gammes d'échelles séparées par une rupture à 16 km qui est discutée. Ces données sont globalement en accord avec un modèle spatio-temporel simple reposant un exposant d'anisotropie entre l'espace et de temps. Les résultats suggèrent une possible universalité des paramètres UM pour les précipitations. Cette thèse aborde également un autre aspect de l'intermittence, particulièrement important pour les longues séries temporelles pluviométriques, que sont les nombreuses mesures nulles de la pluie (c'est-à-dire un pixel où aucune pluie n'est relevée), i.e. les longues périodes sèches. L'ancienne question de la source de cette intermittence, et notamment la nécessité d'un modèle dédié, est revisitée. D'abord les effets d'un seuil sur un champ multifractal sont analysés et ensuite un « toy model » qui introduit des zéros au sein du processus de cascade et conditionnellement aux valeurs du champ est développé. Cela permet d'expliquer la plupart des comportements observés, e.g. les différences entre les statistiques évènementielles et globales. L'impact de la variabilité de la pluie est analysé à travers l'étude de la sensibilité de modèles d'hydrologie/hydraulique urbaine à la donnée de pluie. Deux bassins versants essentiellement urbains (un de 3 400 ha en Seine-Saint-Denis à proximité de Paris, et un de 900 ha à Londres) modélisés avec des modèles opérationnels semi-distribués sont pris comme cas d'études. Par ailleurs le modèle distribué Multi-Hydro (en développement au LEESU) est testé sur une portion de 145 ha du cas d'étude parisien. L'impact de la variabilité à petites échelles non mesurée des précipitations (i.e. se produisant à des échelles plus petites que 1 km en espace et 5 min en temps qui sont disponibles avec les données radar à bande C) est d'abord évalué. Ceci est réalisé par la génération d'un ensemble de pluie réaliste désagrégée en continuant stochastiquement le processus sous-jacent de cascade au-delà de l'échelle d'observation, puis la simulation de l'ensemble correspondant d'hydrographes. Il apparaît que la variabilité à petites échelles de la pluie engendre une variabilité hydrologique qui ne doit pas être négligée. De plus le modèle Multi-Hydro génère une variabilité plus importante et pas seulement au niveau du pic de débit, i.e. même pour les pluies modérées. Ces résultats mettent en lumière la nécessité d'installer des radars en bande X (dont la résolution est hectométrique) en milieu urbain. Dans un deuxième temps les outils multifractals sont employés sur les pluies et les débits simulés qui présentent aussi un comportement scalant. Il apparaît que le réseau d'assainissement transmet simplement la variabilité des précipitations sans l'atténuer, au moins en termes de statistiques multifractals / The Universal Multifractals (UM) are a powerful tool which has been extensively used to analyze and simulate geophysical fields, such as rainfall, that are extremely variable over wide range of scales. It is based on the concept of cascade phenomenology that relies on the physical notion of scale invariance to explore the fundamental phenomenon of intermittency. In this framework the whole variability of a field is characterized with the help of only three parameters that are furthermore physically meaningful. In this PhD thesis we use this theoretical framework to quantify the impacts of small scale rainfall variability in urban hydrology. The first step consists in analysing radar rainfall space-time variability with the help of a simple anisotropic multifractal model. A variety of rainfall events are analyzed. It appears that a scaling behaviour was observed on two distinct ranges of scales separated by a break at roughly 16 km that is discussed. These data sets are in overall agreement with a simple space-time scaling model relying on single anisotropy exponent between space and time. The results hint at a possible universality of the UM parameters for rainfall. This thesis also explores another facet of intermittency, which is particularly important for long time series of precipitation, that of numerous zero rainfall measurements (a pixel or a time step with no recorded rainfall), i. e. long “dry” periods. We revisit the long lasting discussion on the source of this intermittency, e.g. whether it requires a specific modelling. First the effects of a threshold on a universal multifractal field are investigated and second a toy model that introduces some zeros within the cascade process conditioned by the field value is developed. This enables to explain most of the observed behaviour, e.g. the difference between event statistics and overall statistics. The impact of rainfall variability is investigated through the analysis of the sensitivity to the rainfall input of urban hydrologic-hydraulic models. Two predominantly urban catchments (a 3 400 ha one in Seine-Saint-Denis near Paris, and a 900 ha one in London) modelled with the help of operational semi-distributed models are used as case studies. The fully distributed model Multi-Hydro (under development at LEESU) is also tested on a 147 ha portion of the Paris case study. First the impact of unmeasured small scale rainfall variability (i.e. occurring at scales smaller than 1 km in space and 5 min in time which are available with C-band radar data) is evaluated. This is achieved by generating an ensemble of realistic downscaled rainfall fields by continuing the stochastic cascade process below the observation scale and then simulating the corresponding ensemble of hydrographs. It appears that the small scale rainfall variability generates significant hydrological variability that should not be neglected. Furthermore the Multi-Hydro model generates a larger variability not only during the peak flow, but during the whole event, i.e. for moderate rain rates. These findings highlight the need to implement X-band radars (whose resolution is hectometric) in urban areas. In a second part multifractal tools are used on both rainfall and simulated discharges that also exhibit a scaling behaviour. It appears that the rainfall drainage system basically transmits the rainfall variability without damping it, at least in terms of multifractal statistics

L’hydrologie, la géochimie et la météorologie du petit bassin versant glaciaire à pergélisol continu du glacier Austre Lovén (10 km2, Spitsberg occidental, 79°N) ont été étudiées au cours des trois années hydrologiques 2010, 2011 et 2012 en vue de mieux définir le rôle des eaux souterraines de la nappe supra-pergélisol dans les flux d’eau et de solutés en sortie de bassin.Le travail de thèse a montré que la période d’écoulement survient en moyenne sur une période de 141 jours (σ = 4,5 jours) comprise entre les 9-10 mai et 26-27 octobre. La lame d’eau moyenne est de 0,940 m.a-1 (σ = 0,156 m.a-1) et se répartie à 79 % en eau de fonte (54 % de fonte de neige et de glace sur le glacier, 18 % de neige sur les versants et 7 % de neige dans la zone proglaciaire), 18 % en hauteurs de précipitation efficace et 3 % en décharge sous-glaciaire pérenne. Les incertitudes sur la quantification des composantes de l’écoulement et de la lame d’eau écoulée elle-même ne permettent pas une quantification plus précise des flux en eau de la nappe supra-pergélisol. Cette étude montrerait donc qu’en termes de flux d’eau, la composante d’eau souterraine supra-pergélisol pourrait être négligeable et ne constituer qu’un réservoir de transit de l’eau. En effet, le suivi piézométrique montre que la nappe dégèle et regèle à des profondeurs similaires malgré la présence d’épisodes de recharge. Ces recharges qui surviennent majoritairement au profit d’épisodes pluvieux importants (hauteurs cumulées sur 2 à 8 jours supérieures à 20 mm) se déchargent lentement par une alimentation de la rivière par la nappe. Les suivis thermiques du sous-sol ont enfin permis d’établir que le toit du pergélisol supposé se situe à une profondeur de 2,50 m pour une puissance de nappe pouvant atteindre plus d’1,70 m, soit 70 % de la couche active. Par ailleurs, les données thermiques montrent des températures du sous-sol à moins de 0,5°C de la température de dégel, laissant supposer une dégradation prochaine du pergélisol avec par conséquent un accroissement de l’épaisseur de l’aquifère de la nappe supra-pergélisol (couche active).Les données de géochimie montrent une augmentation de la minéralisation des eaux de surface du bassin par des contributions d’eau souterraine le long des cours d’eau, preuve que les eaux de la nappe supra-pergélisol impactent les flux de solutés aux exutoires. Les analyses en chimie élémentaire (n = 178) conduites sur les différents composants du système hydrologique montrent que les eaux du bassin possèdent un faciès sulfaté-bicarbonaté calcique à bicarbonaté-sulfaté calcique lié à l’altération des carbonates et des sulfures des moraines récentes et dans l’altération d’aluminosilicates dans les moraines plus anciennes. Les données d’isotopie du strontium (⁸⁷Sr/⁸⁶Sr  ; n = 8) corroborent un apport du calcium dissous par les aluminosilicates de par leur signature radiogénique. Les teneurs en ¹³C du CID des eaux (n = 30) sont cohérentes aussi bien avec un système ouvert sur le CO₂ du sol (à -20 ‰) et le CO₂ atmosphérique (entre -6,5 et -8 ‰) qu’avec un système fermé sur le CO₂ où la signature de δ ¹³C des carbonates marins (à 0 ‰) peut s’imposer par dissolution. Les données d’isotopie de la molécule d’eau (n = 592) ont permis d’identifier les signatures des différentes familles d’eau composant l’hydrosystème et de proposer un modèle de mélange entre 4 pôles expliquant les signatures des eaux aux exutoires du bassin versant : eaux de fonte, précipitations estivales, eaux sous-glaciaires et eaux de la nappe supra-pergélisol. Le couplage entre les teneurs en ¹⁸O et SO²₄⁻ des eaux aux exutoires confirme ce modèle de mélange, le renforçant même en termes de variations temporelles. La mesure de ³H-³He tendrait vers un temps de séjour court tandis que les premières mesures d’activité ¹⁴C (n = 15) sont faibles, renseignant une minéralisation des eaux par un carbone ancien. / The hydrology, geochemistry and meteorology of a small, Arctic glacial watershed under continuous permafrost conditions (Austre Lovén glacier catchment area of 10 km2, western Spitsbergen, 79°N) were studied during the three hydrological years 2010, 2011 and 2012 to better define the role of suprapermafrost groundwater on both water and solute fluxes at basin outlet. The runoff period occurs on average over a 141 days period (σ = 4.5 days) between May 9-10 and October 26-27. The mean total discharge is 0.940 m a-1 (σ = 0.156 m a-1) divided into 79 % of meltwater (54 % of snowmelt and icemelt from the glacier, 18 % of snowmelt from the slopes and 7 % of snowmelt from the proglacial area), 18 % of effective summer precipitation and 3% of perennial subglacial discharge. Uncertainties in quantifying the runoff components and the total discharge itself don’t allow more precise quantification of water fluxes from the suprapermafrost groundwater by water balance. This study would so show that the suprapermafrost groundwater component would be negligible in the water balance and only constitute a transit tank for water. Indeed, the potentiometric level monitoring shows that the water-table thaws and freezes at similar depths despite recharge events occurring during summer. Those recharge events mainly occur in favor of important rain events (cumulative amounts on 2 to 8 days higher that 20 mm). The water-table discharges towards rivers. The monitoring of ground temperature indicates that the frozen ground (supposed permafrost) top is located at a maximum depth of 2.50 m for a groundwater thickness reaching up to 1.70 m, or 70 % of the active layer. Thermal data show ground temperatures between 0 and - 0.5°C, close to thaw temperature, which let suppose a permafrost degradation between 2.50 and 3.40 m deep is about to occur due to climate warming. Consequently, the suprapermafrost groundwater aquifer (active layer) is expected to thicken. Geochemical data show an increase of the basin surface water mineralization by groundwater contribution along the rivers, proof of a suprapermafrost groundwater impact on the solutes fluxes at the outlets. Elementary chemical analyses (n = 178) performed on the different components of the hydrological system show that basin waters have a sulfate-bicarbonate calcium to bicarbonate-sulfate calcium type linked to carbonates dissolution and sulfides (pyrite) oxydation in the recent moraines and to aluminosilicates weathering in older moraines. Strontium isotopic data (87Sr/86Sr; n = 8) are consistent with a contribution in dissolved calcium from aluminosilicates according to their radiogenic signature. The ¹³C contents of water DIC (n = 30) are consistent just as well with a system open on the soil CO₂ (at -20 ‰) and the atmospheric CO₂ (between -6.5 et -8 ‰) than with a system closed on the CO₂ where the δ ¹³C of marine carbonates (at 0 ‰) control the δ ¹³C signature by dissolution. Isotopic data of the water molecule (n = 592) are helpful to separate the signatures of the different water masses of the hydrosystem and to propose a 4 end-members mixing model explaining the river water signature: meltwaters, summer precipitations, subglacial water and suprapermafrost groundwater. The coupling between river water contents in ¹⁸O and SO²₄⁻ confirms this mixing model, even reinforcing it in terms of temporal variations. The unique measurement of ³H-³He tends to a short residence time while the first ¹⁴C activity data (n = 15) are low, giving a water mineralization by an old carbon.