# Bilans de masse du Glacier Blanc entre 1952, 1981 et 2002 obtenus par modèles numériques de terrain

Emmanuel Thibert, Joël Faure, Christian Vincent

► To cite this version:

Emmanuel Thibert, Joël Faure, Christian Vincent. Bilans de masse du Glacier Blanc entre 1952, 1981 et 2002 obtenus par modèles numériques de terrain. La Houille Blanche - Revue internationale de l'eau, 2005, 2, pp.72-78. 10.1051/lhb:200502010. insu-00421361

HAL Id: insu-00421361

https://insu.hal.science/insu-00421361v1

Submitted on 22 Nov 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY-SA 4.0 - Attribution - ShareAlike - International
License
<figure><img src="image_2.png" /></figure>

# Bilans de masse du Glacier Blanc entre 1952, 1981 et 2002 obtenus par modèles numériques de terrain

Emmanuel Thibert, Joël Faure & Christian Vincent

To cite this article: Emmanuel Thibert, Joël Faure & Christian Vincent (2005) Bilans de masse du Glacier Blanc entre 1952, 1981 et 2002 obtenus par modèles numériques de terrain, La Houille Blanche, 91:2, 72-78, DOI: 10.1051/lhb:200502010

To link to this article: https://doi.org/10.1051/lhb:200502010

<figure><img src="image_3.png" /></figure>

© Société Hydrotechnique de France, 2005

<figure><img src="image_4.png" /></figure>

Published online: 01 Jun 2007.

<figure><img src="image_5.png" /></figure>

Submit your article to this journal

<figure><img src="image_6.png" /></figure>

Article views: 15

<figure><img src="image_7.png" /></figure>

View related articles

<figure><img src="image_8.png" /></figure>

Citing articles: 1 View citing articles

<figure><img src="image_9.png" /></figure>
# Bilans de masse du Glacier Blanc entre 1952, 1981 et 2002 obtenus par modèles numériques de terrain

*Mass balances of the Glacier Blanc between 1952, 1981 and 2002 calculated from digital elevation models*

**Emmanuel Thibert**

Parc National des Ecrins, Domaine de Charance, 05000 Gap
CEMAGREF-ETNA, BP 76, 38402 Saint Martin d'Hères cedex
emmanuel.thibert@cemagref.fr

**Joël Faure**

Parc National des Ecrins, Domaine de Charance, 05000 Gap

**Christian Vincent**

LGGE, CNRS-UJF, BP 96, 38402 Saint Martin d'Hères cedex

Three Digital Elevation Models (DEM) of the Glacier Blanc have been calculated using photogrammetric restitutions of aerial photographs realised in 1952, 1981 and 2002 by IGN (the french administration for geography) and SINTEGRA (private company, Meylan, France). The aerial triangulation and the orientation of the stereographic images based on GPS control points lead to an absolute orientation whose residual standard deviation is between 0.25 to 1.25 m in planimetry and between 0.4 and 1.55 m in altitude (Lambert III, IGN69). Spot heights and 5 meter contour lines have been drawn (120 points/ha) which enables to build a Triangular Irregular Network (TIN) according to Delauney algorithm. This vector network has been linearly interpolated to obtain a DEM which is a raster grid of 10 meter spatial resolution. The subtraction of those DEM yields altitudinal variations on the surface of the glacier from which cumulative mass balance has been deduced. Results show a mass balance increasing of 2.1 m between 1952 and 1981 and decreasing of 10.8 m between 1981 and 2002. Results are compared with those of neighbour glaciers.

# I ■ INTRODUCTION

## ● I.1 Objectifs

Les glaciers de montagne sont maintenant largement recon-
nus comme d'excellents indicateurs des changements climati-
ques survenus lors des derniers siècles (IPCC, 2001 ; Haeberli,
1995). En particulier, les variations de masse des glaciers peu-
vent être utilisées pour témoigner du réchauffement du climat
au cours du XXe siècle (Vincent et al., 2004). Contrairement
aux relevés des variations de longueur des glaciers qui
sont les données les plus anciennes et les plus répandues
(Oerlemans et al., 1998) mais qui fournissent un signal clima-
tique déphasé, filtré (coupure des hautes fréquences) et convo-

lué à l'écoulement du glacier (Paterson, 1994 ; Johannesson,
1989), les variations de masse sont des signaux climatiques
directs. En effet, la plupart des glaciers des Alpes étant tem-
pérés, ils constituent d'un point de vue thermodynamique
des systèmes « ouverts » (pouvant échanger de la matière et
de l'énergie) dont la température est au point de fusion. Les
: changes d'énergie se manifestent donc directement par une
variation de masse (ablation, accumulation interne). Le bilan
de masse enregistre également, via l'accumulation, les préci-
pitations de neige de l'hiver. Le suivi du bilan masse permet
donc d'accéder aux variations de ces paramètres climatiques.
Malheureusement, très peu de glaciers font l'objet d'un tel
suivi sur une longue période (Haeberli et al., 1998). Il est
néanmoins possible d'accéder au bilan de masse cumulé entre
## Glacier Blanc
Plan photogrammétrique

Etat en octobre 2002

<figure><img src="image_10.png" /><figcaption>Figure 1 : Plan photogrammétrique du Glacier Blanc dressé à partir de la restitution des clichés SINTEGRA du 2 octobre 2002.</figcaption></figure>

2 dates données en utilisant des documents comme les cartes topographiques et des photographies aériennes.

De tels documents sont disponibles sur le Glacier Blanc pour réaliser un bilan géométrique du glacier et d'en tirer le bilan de masse. Ce glacier (44°56'N, 6°23'E, Massif des Ecrins, France) est le plus grand glacier du Massif des Ecrins. En octobre 2002, il couvre une surface est de 534 ha, et s'étend sur 5,9 km de long (distance curviligne). Il prend naissance à 4 014 m au Dôme des Ecrins sur le versant Nord-Est de la Barre des Ecrins (4 102 m), versant dont la pente moyenne est de 30°. Un vaste plateau (pente moyenne inférieure à 10°) orienté Est entre 3 000 et 3 400 m d'altitude constitue sur 350 ha la plus grande partie du bassin d'accumulation. En dessous de 3 000 m, la langue d'ablation orientée Sud-Est (pente moyenne 20°) se termine au front à 2 315 m d'altitude (fig. 1).

Le Glacier Blanc est par sa taille le plus important gla-cier méridional de l'Europe de l'Ouest ayant fait l'objet de suivis plus ou moins réguliers. Les relevés annuels de fluctuations du front les plus anciens remontent à 1889-1907 (Commission Française des Glaciers ; Allix 1927). Entre 1921 et 1965, ces relevés sont assurés par le Service RTM de l'Administration des Eaux et Forêts, et sont associés, à partir de 1922, à des mesures de vitesse sur deux profils transversaux situés sur la langue d'ablation (Mougin, 1934 ; Poncet, 1966). Un troisième profil est installé en amont de la chute de sérac en 1950. Interrompus en 1966, ces relevés de vitesse et de position du front sont repris en 1983 par le laboratoire de Glaciologie (CNRS) jusqu'en 1995 (date où le glacier, en retrait, disparaît du profil aval), puis par le Parc National des Ecrins (Thibert et al., 2003 ; fig. 2). A partir de 1975 et jusqu'en 1991, le Cemagref débute des campagnes périodiques (3 ans) de photogrammétrie aérien-nes qui couvrent uniquement la langue d'ablation du glacier (Burnet et Mura, 1985). En 1968, le bassin supérieur fait

<figure><img src="image_11.png" /><figcaption>Figure 2 : Cumul des variations de longueur et surface du Glacier Blanc entre 1815 et 2003 (d'après Mougin, 1934 ; Letreguilly, 1983 ; Reynaud, 1995 ; Thibert et al., 2003 ; et cette étude).</figcaption></figure>

l'objet d'une campagne de sismique (Gluck, 1969) qui a pour objet de déterminer l'épaisseur de glace sur 2 profils du bassin d'accumulation. Des forages sont également réalisés sur la langue d'ablation pour déterminer quelques épaisseurs (Benoist, 1972). Une synthèse de ces travaux est effectuée par Letréguilly (1983) et Reynaud (1995). Les fluctuations du bilan de masse ont été reconstituées par télédétection sur la période 1985-2000 par Rabatel et al. (2002). Depuis 2000, les bilans de masse sont mesurés dans la zone d'accumulation à l'aide de carottage et de perches d'ablation (Thibert et al., 2003).
## I.2 Principe de calcul du bilan à partir de Modèles Numériques de Terrain (MNT)

Ce calcul du bilan de masse entre 2 états du glacier nécessite la connaissance de l'altitude des points du glacier et également des rives rocheuses qui ont été délaissées et/ou recouvertes entre ces 2 états. La fonction altitude $Z(x, y)$ étant connue sur cette surface $S_0$ et exprimée dans un système géoréférencé (coordonnées planimétriques Lambert III et altitude IGN69 dans notre cas), la variation de masse du glacier s'exprimera à partir de la différence $\Delta z$ d'altitude entre les 2 états par la relation :

$$ \Delta m = \iint_{S_0} \Delta z(x, y) \rho(x, y) dxdy = \sum_{i} \sum_{j} \Delta h_{ij} \rho_{ij} s_{pixel} $$

où $\rho(x, y)$ désigne la densité moyenne du matériau (glace, névé ou neige) sur l'épaisseur en question. Cette variation de masse peut être également écrite de manière discrète lorsque la fonction altitude n'est connue que sur une grille de points, ce qui sera le cas avec un modèle numérique de terrain qui fournit la matrice $\Delta h_{ij}$ pour chaque pixel de surface $s_{pixel}$. Dans les deux cas, la variation de masse s'écrit comme le produit de la surface $S_0$ par la valeur moyenne du produit de la variation d'altitude par la densité. Le bilan global net entre les deux états est donc donné directement par :

$$ \Delta b = \frac{\Delta m}{S_0} = \langle \rho \Delta h \rangle_{S_0} $$

où $\langle \rho \Delta h \rangle_{S_0}$ désigne la valeur moyenne sur la surface $S_0$ du produit de la densité par la variation d'altitude. Il est important de noter ici que la valeur moyenne $\langle \rho \Delta h \rangle$ ne fournit le bilan net que si la surface recouvre l'ensemble du glacier. Si il ne s'agit que d'une section du glacier, il faut appliquer l'équation de continuité (conservation de la masse) et tenir compte des flux de glace à l'entrée et à la sortie de la section pour relier le bilan à la variation de hauteur moyenne de la section (Reynaud et al., 1986).

## II ■ MÉTHODE

Les documents disponibles sur le glacier Blanc pour réaliser ce type de travail sont : La carte au 1/10 000 de 1904 de Flusin et al. (1905), la carte du Service Géographique des Armées au 1/20 000 de 1928 (SGA, 1933), puis les couvertures aériennes de l'Institut Géographique National (IGN 1952, 1967, 1981, 1993) et enfin une couverture de SINTEGRA réalisée en 2002 spécialement pour cette étude. En complément de la couverture SINTEGRA de 2002 servant de référence, notre choix s'est porté sur :

- les clichés IGN datant de 1952, première couverture aérienne du massif des Ecrins,

- les clichés IGN de 1981 qui présentent l'intérêt d'être contemporain du maximum de bilan de masse mesuré à cette époque sur les glaciers de Sarennes et Saint Sorlin (Vincent, 2002).

L'extension des prises de vues réalisées par SINTEGRA a été choisie de manière à couvrir les terrains occupés par le glacier Blanc tels qu'ils figurent sur la carte de Flusin et al.

<figure><img src="image_12.png" /><figcaption>Figure 3 : Orthophotoplan du glacier Blanc obtenu à partir des prises de vues du 2 octobre 2002. Les croix blanches sont les points utilisés pour le calage des 3 couvertures aériennes de 1952, 1981 et 2002. Le trait noir indique la limite des terrains restitués pour obtenir le modèle numérique de terrain de référence.</figcaption></figure>

(1905) et de pouvoir ainsi calculer ultérieurement un bilan de masse à partir de ce document (fig. 3).

### II.1 Paramètres des prises de vue

Les clichés de 1952 proviennent de la mission La Grave — Chorges F3435-3438 du 1er août 1952, soit 6 clichés sur deux axes de vol, échelle au centre 1/25 000. Les clichés de 1981 proviennent de la mission IFN05 P (inventaire forestier national, département des Hautes Alpes) du 29 juillet 1981. Deux axes de vol couvrent le glacier en 8 clichés dont l'échelle au centre est d'environ 1/17 000. La couverture SINTEGRA a été effectuée en fin de période d'ablation, le 3 octobre 2002 sur un axe de vol comportant 5 clichés au 1/20 000.

### II.2 Stéréopréparation

La campagne de stéréopréparation a été effectuée en août 2001. Compte tenu de l'antériorité des prises de vues IGN, des points remarquable du relief (gendarmes, blocs de rochers, toiture des refuges) ont été choisis d'après leur visibilité sur les anciens clichés. En complément, des points ont été implantés (croix de peinture blanche et bâches) en prévision de la couverture aérienne SINTEGRA de 2002. Dix-huit points d'appui ont été utilisés (fig. 3). Ils ont été rattachés au système planimétrique NTF — Lambert III et altimétrique NGF — IGN69 par géodésie GPS différentiel bi-fréquence (statique rapide post-traité). La précision est de l'ordre de quelques cm en planimétrie et d'une dizaine de centimètres en altimétrie relativement aux points géodésiques utilisés (IGN, réseau REGAL).

### II.3 Calage des clichés

Ce travail ainsi que la restitution ont été réalisés par SINTEGRA sur un restituteur de type analytique (stéréorestituteur M/s LEICA DSR 15). Les couvertures étant
<table><caption>Tableau 1. Erreur de calage des clichés σ<sub>xy</sub> (planimétrique) et σ<sub>z</sub> (altimétrique), erreur altimétrique résultante ΔZ pour le MNT par rapport aux points d'appui, incertitude sur les variations d'altitude ΔH entre le MNT (1952, 1981) et le MNT de référence (2002) et incertitude résultante Δb sur le bilan de masse calculé par rapport à l'année 2002.</caption><thead><tr><td>Année/précision</td><td>σ<sub>xy</sub></td><td>σ<sub>z</sub></td><td>ΔZ</td><td>ΔH</td><td>Δb</td></tr></thead><tbody><tr><td>1952</td><td>1,25 m</td><td>1,55 m</td><td>2,0 m</td><td>2,1 m</td><td>1,9 m eq. eau</td></tr><tr><td>1981</td><td>0,25 m</td><td>0,40 m</td><td>0,50 m</td><td>0,96 m</td><td>0,9 m eq. eau</td></tr><tr><td>2002</td><td>0,45 m</td><td>0,55 m</td><td>0,72 m</td><td>—</td><td>—</td></tr></tbody></table>

parfois constituées de plusieurs couples stéréo sur 2 axes de vol, et le nombre de points d'appui étant réduit, il a d'abord été réalisé une aérotriangulation (calages des clichés entre eux) puis un calage de cet ensemble par rapport au réseau géodésique (calcul en bloc par un programme de compensation type moindres carrés). La précision en planimétrie et altimétrie (σ<sub>xy</sub>σ<sub>z</sub>) de l'orientation absolue des clichés ainsi obtenue est indiquée au tableau 1.

## II.4 Restitution photogrammétrique

Afin d'établir les Modèles Numériques de Terrain (MNT) du glacier, les courbes de niveau ont été filées tous les 5 mètres (sur le glacier, les rives rocheuses et les terrains déglacés depuis 1904). Des détails bien visibles par l'opérateur dans le stéréomodèle ont été également pointés pour fournir des points côtés entre les courbes de niveau (certaines crevasses et séracs, gendarmes, crêtes, talwegs et blocs rocheux). Avec la densité de points obtenue (environment 120 points/ha), il est alors possible d'obtenir un plan photogrammétrique au 1/5 000 et un modèle numérique de terrain jusqu'à 2,5 mètres de résolution spatiale. La figure 1 illustre le plan obtenu à partir des clichés d'octobre 2002. L'orthophotoplan issu des mêmes clichés est représenté à la figure 3.

<figure><img src="image_13.png" /><figcaption>Figure 4 : Réseau de triangle irrégulier (TIN) construit sur les courbes de niveau et les points côtés. Aucun côté de triangle ne coupe les courbes de niveau.</figcaption></figure>

## II.5 Modèle numérique de terrain

Le MNT vecteur construit est de type réseau de triangle irrégulier (Triangular Irregular Network = TIN) s'appuyant sur les angles des courbes de niveau et les points côtés (cf. fig. 4). L'algorithme retenu pour construire ce réseau est celui de Delauney (les angles des triangles sont maximisés). Les algorithmes classiques de triangulation spécifiques aux courbes de niveau (triangulation contrainte) n'ont pas été retenus pour pouvoir inclure les points côtés dans ce réseau. Néanmoins, la densité des courbes de niveau est telle que tous les segments de courbes de niveau constituent le côté d'un triangle. Par conséquent, aucune courbe de niveau n'est coupée par le côté d'un triangle. Ceci assure la continuité entre les altitudes fournies par les courbes de niveau et le réseau de triangle. Le résultat du calcul de triangulation a donc de ce point de vue la même propriété que celui des algorithmes de triangulation contraintes et incorpore en plus les points côtés.

Le réseau de triangle ainsi obtenu a été interpolé sur une grille de 10 m (quadrillage Lambert III) pour obtenir le modèle numérique d'altimétrie (530 ×420 pixels). L'interpolation est linéaire : l'altitude de chaque pixel est celle du point du triangle intercepté par le centre du pixel. Ce MNT est donc un fichier de type image (raster) dont l'attribut est l'altitude (IGN 69). La soustraction des MNT donne les variations d'altitude et donc le bilan de volume

<figure><img src="image_14.png" /><figcaption>Figure 5 : Carte des variations d'altitude du glacier entre 1952 et 2002 représenté en 256 niveaux de gris. En noir figurent les courbes d'isovaleur des pertes d'altitudes à 20 m (entre –5 et –85 m). La perte moyenne d'altitude du glacier est de 11,2 mètres.</figcaption></figure>
<figure><img src="image_15.png" /><figcaption>Figure 6: Tracé des courbes de niveau obtenues par le stéréorestituteur (en trait noir continu) et par le modèle numérique de terrain à 10 m de résolution (en pointillé). Dans les zones au relief complexe, comme ici dans la chute de séracs, les courbes de niveau du MNT ont tendance à simplifier le relief à petite échelle (en dessous de 20 m) en supprimant les crêtes et talwegs.</figcaption></figure>

entre deux états du glacier. La figure 5 représente les variations d'altitude entre 1952 et 2002 et les courbes d'isovaleur de variation d'altitude à 20 m.

### III ■ RÉSULTATS

#### ● III.1 Fidélité du MNT au plan photogrammétrique

Elle a été évaluée en comparant les courbes de niveau calculées d'après le MNT à celles obtenues directement par le stéréo-restituteur. Les discordances de tracé sont inférieures à 1 m dans la plupart des cas. Seul dans les zones où le relief est complexe à petite échelle (en dessous de 20 m), le maillage d'interpolation du MNT (10 m) moyenne le relief. Les courbes de niveau du MNT ont un rayon de courbure qui est plus élevé que celui des courbes données par le stéréo-restituteur (fig. 6). Les crêtes et talwegs de petite échelle (< 20 m) ont tendance à être éliminés. La densité de point permettrait d'augmenter la résolution du MNT jusqu'à un mètre dans ces zones ou d'adopter une grille à pas variable. Cependant, l'information recherchée étant le bilan global net sur l'ensemble du glacier, nous avons conservé le maillage à 10 mètres.

#### ● III.2 Incertitude sur les altitudes

L'erreur globale $\Delta Z$ associée à l'altitude fournie par le MNT $Z(x, y)$ provient a) de l'erreur de calage altimétrique des clichés, b) de l'erreur de calage planimétrique qui induit dans les zones à relief une erreur d'altitude du premier ordre et c) de la résolution de la grille du MNT. Cette dernière erreur semble négligeable compte tenue de la bonne fidélité au plan photogrammétrique du MNT à 10 m. Les deux premières contributions peuvent s'écrire :

$$ \Delta Z = (\Delta z)_{calage} + \left(\frac{\partial Z}{\partial x}\right)_y \Delta x + \left(\frac{\partial Z}{\partial y}\right)_x \Delta y = \sigma_z + \sigma_{xy} \text{grad}Z $$

où $\sigma_z$ désigne l'erreur de calage altimétrique des clichés, $\sigma_{xy}$ l'erreur de calage planimétrique et $\text{grad}Z$ le gradient moyen du MNT sur le glacier. Avec un gradient moyen de 38 % pour les trois années étudiées, cette erreur peut se calculer en fonction des erreurs de calage de chaque ensemble de clichés. Il en résulte l'erreur $\Delta H$ sur les variations d'altitude entre les MNT de 1952 et 1981 par rapport au MNT de référence de 2002, et l'erreur résultante $\Delta h$ sur le bilan de masse (tab. 1).

#### ● III.3 Calcul du bilan de masse

La différence d'altitude donnée par les MNT doit être convertie en lame d'eau moyennant une hypothèse sur la densité. On a ici simplement pris comme facteur 0,9 sur l'ensemble du glacier pour passer à l'équivalent en eau.

Nous avons tenu compte des dates précoces des prises de vue de l'IGN (1er août 1952 et 29 juillet 1981) pour faire une correction de lame d'eau correspondant à l'ablation jusqu'à début octobre, mois de la prise de vue des clichés de référence de 2002. Pour ceci, nous avons calculé le cumul de température positive sur août et septembre (degrés jours positifs) pour chaque point du glacier en utilisant les normales saisonnières de la station MétéoFrance de Pelvoux les Claux (1 270 m).

Le MNT fournissant l'altitude, le cumul de température est calculé en utilisant le gradient adiabatique de 5,5 $10^{-3} \,^\circ\text{C/m}$. Nous utilisons 2 facteurs de sensibilité correspondant au matériau glace en dessous de 3 000 m (0,7 $\text{cm}/^\circ\text{C/jour}$, Torinesi et al., 2001) et au névé au dessus de 3 000 m (0,5 $\text{cm}/^\circ\text{C/jour}$, Thibert et al., 2003). Cette correction donne à titre d'exemple environ –3,8 m d'eau à 2 250 m et –1,2 m d'eau à 3 100 m pour la période août-septembre. La lame d'eau de fonte pondérée sur l'ensemble du glacier est de –1,4 m pour les deux années 1952 et 1981. Les résultats sont rassemblés au tableau 2.

<table><caption>Tableau 2. Bilan volumique et bilan de masse du glacier Blanc pour les périodes 1952-2002 et 1981-2002.</caption><thead><tr><td>Période</td><td>Bilan volumique (m)</td><td>Bilan de masse (m eq. Eau)</td><td>Correction d'ablation (m eq. Eau)</td><td>bilan corrigé (m eq. Eau)</td></tr></thead><tbody><tr><td>1952-2002</td><td>-11,2</td><td>-10,1</td><td>1,4</td><td>-8,7 +/- 1,9</td></tr><tr><td>1981-2002</td><td>-13,5</td><td>-12,2</td><td>1,4</td><td>-10,8 +/- 0,9</td></tr></tbody></table>
IV ■ DISCUSSION

En terme de variation du stock d'eau, ces chiffres cor-
respondent à une perte de 52 millions de m³ d'eau entre
1952 et 2002. Le glacier blanc a regagné 12,4 millions de
m³ entre 1952 et 1981. La perte de volume d'eau se chiffre
à 63,4 millions de m³ entre 1981 et 2002. La surface du
glacier de 599 ha en 1952, s'abaisse à 587 ha en 1981 (alors
que son volume a augmenté), et passe à 534 ha en 2002.

Les résultats obtenus pour le glacier Blanc sont compa-
rés aux mesures de bilan cumulé et de bilan obtenu égale-
ment par photogrammétrie sur les glaciers de Sarennes, St
Sorlin et Argentière (tab. 3 et fig. 7, d'après Vincent, 2002 et
Michel et al., 2003).

Sur la période 1981-2002, le bilan moyen annuel est de
– 0,51 m eq. eau/an au glacier Blanc, identique à celui du
glacier d’Argentière, mais environ deux fois moins élevé
qu’au glacier de Sarennes (– 0,98 m eq. eau/an), le glacier
de St Sorlin ayant un comportement intermédiaire avec
– 0,77 m eq. eau/an.

Sur la période 1952-1981, le glacier Blanc voit son bilan augmenter de 2,1 m eq. eau. Les bilans sont également positifs pour Argentière (+ 5,9 m eq. eau) et St Sorlin (+ 0,4 m eq. eau), alors que Sarennes est en déficit (– 10,9 m eq. eau).

Malgré sa position méridionale, le glacier Blanc semble avoir un comportement proche de celui du glacier d'Argentière. Plusieurs éléments peuvent être avancés pour expliquer cette similitude : tout comme Argentière, le glacier Blanc est un glacier de grande taille (534 ha et 5,9 km de long en 2002) disposant d'un vaste bassin d'accumulation à haute altitude (80 % de la surface est entre 3 000 et 4 000 m). Ce bassin est situé sous les plus hauts sommets du massif et l'effet orographique doit probablement forcer la pluvio-métrie. Les mesures d'accumulation effectuées sur le bas-sin indiquent en effet, à 3 100 m, deux fois et demi plus de précipitation que ce qui est mesuré par MétéoFrance à la station de Pelvoux les Claux à 1 270 m (Thibert et al., 2003). Ainsi, les glaciers qui ont des zones d'accumulation

élevées sont beaucoup moins sensibles aux variations du
bilan d'énergie estival. Et, ce sont d'ailleurs principalement
ces conditions estivales qui ont affecté les glaciers alpins au
cours du XXe siècle et, en particulier, au cours des 20 der-
nières années (Vincent et al., 2004). Une autre hypothèse
pouvant expliquer une plus faible variabilité du bilan serait
l'existence d'une zone d'accumulation froide dans la partie
supérieure de la face Nord Est des Ecrins. Son réchauffe-
ment aurait alors absorbé une partie de l'énergie qui n'aurait
pas contribué à la fusion de la glace. Mais cette hypothèse
n'est absolument pas vérifiée pour l'instant.

<figure><img src="image_16.png" /><figcaption>Figure 7 : Comparaison sur la période 1952-2002 des bilans de masse cumulés obtenus à partir des mesures de terrain (méthode glaciologique : glaciers de Sarennes, St Sorlin et Argentière) et par photogrammétrie (glaciers Blanc, St Sorlin et Argentière). Données d'après Vincent (2002), Michel et al. (2003) et cette étude.</figcaption></figure>

<table><caption>Tableau 3 : Comparaison des bilans de masse mesurés à Sarennes, St Sorlin, Argentière et glacier Blanc.</caption>
   <thead>
    <tr>
     <td>
      Période/bilan (m eq. eau)
     </td>
     <td>
      Sarennes
     </td>
     <td>
      St Sorlin
     </td>
     <td>
      Glacier Blanc
     </td>
     <td>
      Argentière
     </td>
    </tr>
   </thead>
   <tbody>
    <tr>
     <td>
      1952-2002
     </td>
     <td>
      -31,5
     </td>
     <td>
      -15,8
     </td>
     <td>
      -8,7 +/- 1,9
     </td>
     <td>
      -4,8 (1949)
     </td>
    </tr>
    <tr>
     <td>
      1952-1981
     </td>
     <td>
      -10,9
     </td>
     <td>
      +0,4
     </td>
     <td>
      +2,1 +/- 2,1
     </td>
     <td>
      +5,9 (1949)
     </td>
    </tr>
    <tr>
     <td>
      1981-2002
     </td>
     <td>
      -20,6
     </td>
     <td>
      -16,2
     </td>
     <td>
      -10,8 +/- 0,9
     </td>
     <td>
      -10,7
     </td>
    </tr>
   </tbody>
  </table>

REMERCIEMENTS

Les auteurs tiennent à remercier R. Chevalier et C. Roger
pour leur aide lors de la campagne de stéréopréparation,
ainsi que B. Barré et I. Prévitali de la société SINTEGRA
pour la qualité des travaux de restitution. Cette étude a été
financée par le Parc National des Ecrins (Ministère de l'Eco-
logie et du Développement Durable).

BIBLIOGRAPHIE

ALLIX A., (1927). — Mémoire sur les observations glaciologiques faites en Dauphiné jusqu'en 1924. *Etudes glaciologiques*, Tome VI, p. 1-138. Imprimerie Nationale, Paris.

BENOIST J.-P. (1972). — Déformation d'une langue glaciaire (Glacier Blanc), thèse de troisième cycle, Université de Grenoble, publication n° 162 du laboratoire de glaciologie du CNRS, 69 p.
BURNET R. et MURAT R. (1985). — Surveillance des glaciers par photogrammétrie aérienne, Cemagref, division nivologie, Grenoble, 6 p.

FLUSIN G., JACOB C. (1905). — Carte des glaciers Blanc et Noir (levés de terrain de 1904).

GLUCK S. (1969). — Epaisseur du glacier Blanc (Massif de l'Oisans) dans sa partie supérieure. *Comptes Rendus de l'Académie des Sciences*, Série D, n° 268, p. 1583-1585.

HAEBERLI W. (1995). — Glacier fluctuations and climate change detection. *Geogr. Fis. Dinam. Quat.*, 18, 191-199.

HAEBERLI W., HOELZLE M., SUTER S. & FRAUENFELDER R. (1998). — Fluctuation of glaciers, 1990-1995, Vol. VII, UNESCO, Paris.

INSTITUT GÉOGRAPHIQUE NATIONAL (1952). — Mission La Grave-Chorges, (1967) Mission F XXIV-36. (1981) Mission IFN 05 P. (1993) Mission FD 05.

INTERGOUVERNEMENTAL PANEL ON CLIMATE CHANGE (IPCC) (2001). — Climate Change 2001, Contribution of Working Group I to the Third Assessment Report of the Intergouvernemental Panel on Climate Change, Cambridge Univ. Press, New York.

JOHANNESSON T., RAYMOND C. & WADDINGTON E. (1989). — Time-scale for adjustment of glaciers to changes in mass balance. *J. Glaciolo.*, 35, 355-369.

LETRÉGUILLY A. (1983). — Etude des mesures effectuées sur le gla-cier Blanc (Massif des Ecrins), Mémoire de DEA, Laboratoire de Glaciologie (CNRS), Université Joseph Fourier, Grenoble.

MICHEL J.-P., GAY M. et VALLA F. (2003). — Bilan de masse du glacier de Sarennes, 54e saison 2001-2002, SHF, Grenoble.

MOUGIN M. (1934). — Glacier du Pelvoux et Vénéon, p. 71-88 ; Variations du volume des langues terminales : p. 203-210 ; Cinématique glaciaire : p. 252-255, Observations sur la neige et les glaciers des massifs méridionaux du Pelvoux : p. 272-290. Etudes glaciologiques, Paris, Imprimerie Nationale, Tome VII.

OERLEMANS J. et al. (1998). — Modelling the response of glaciers to climate warming. *Clim. Dyn.*, 14, 267-274.

PATERSON W.S.B. (1994). — The physics of glaciers, 3rd ed., Elsevier Sci., New York.

PONCET A. (1966). — Observations glaciologiques sur les gla-ciers Blanc et Noir. Centre technique forestier de Grenoble, 20 décembre 1966, 13 p.

RABATEL A., J.-P. DEDIEU, ET L., REYNAUD (2002). — Reconstitution des fluctuations du bilan de masse du Glacier Blancs (Massif des Écrins, France) entre 1987 et 2000 par télédétection optique (imagerie Spot et Landsat). La Houille Blanche, n° 6/7, p 64.

REYNAUD L., VALLON M. ET LETREGUILLY A. (1986). — Mass balance measurements : problems and two new methods of determining variations. *Journal of Glaciology*, Vol. 32, n° 112.

REYNAUD L. (1995). — Rapport d'étude sur les suivis glaciologiques dans le Parc National des Ecrins, Laboratoire de Glaciologie et de Géophysique de l'Environnement, Parc National des Ecrins, 11 p.

SERVICE GÉOGRAPHIQUE DE L'ARMÉE (1933). — Carte 1/20 000, 34-36, Saint-Christophe en Oisans, feuilles 3 et 4.

THIBERT E., FAUREJ., VINCENT C. ET REYNAUD L. (2003). — Monitoring of glacier Blanc, french Alps. Alpine Glaciological Meeting, Grenoble, France.

TORINESI O., LETRÉGUILLY A. & VALLA F. (2001). — A century reconstruction of the mass balance of Glacier de Sarennes, French Alps. *Journal of Glaciology*, Vol. 48, 142-148.

VINCENT C. (2002). — Influence of climate change over 20th century on four french glacier mass balances. *Journal of Geophysical Research*, Vol. 107, No. D19, 4375.

VINCENT C., KAPPENBERGER G., VALLA F., BAUDER A., FUNK M. & LE MEUR E. (2004). — Ice ablation as evidence of climate change in the Alps over the 20th century. *Journal of Geophysical Research*, Vol. 109, No. D10, 10104.
SOCIETE HYDROTECHNIQUE DE FRANCE

25 rue des Favorites - F 75015 PARIS - Tél. 01-42-50-91-03 - Fax 01-42-50-59-83 - mail shf@shf.asso.fr - www.shf.asso.fr

SIRET 784 309 056 00025 - Code APE 731Z - Code TVA FR23784309056 - CCF PARIS VAUGIRARD 30056-00073-00735402180-23 - CCP Paris : 6148 20 N

<table><thead><tr><th></th><th>Euros</th><th>Nb ex</th></tr></thead><tbody><tr><td colspan="3"><strong>COLLECTION : EAU ET ENVIRONNEMENT</strong></td></tr><tr><td>Crues et Laves torrentielles (1992)</td><td>75,00</td><td></td></tr><tr><td>La pluie, source de vie, choc de pollution (1993)</td><td>30,00</td><td></td></tr><tr><td>Glaciologie et nivologie - Etat des recherches et des connaissances à la fin du XXè s.(1995)</td><td>30,00</td><td></td></tr><tr><td>Le risque crue en région parisienne (1997)</td><td>60,00</td><td></td></tr><tr><td>La gestion des risques liés aux inondations (1999)</td><td>75,00</td><td></td></tr><tr><td>Protection des eaux souterraines (mai 2000)</td><td>75,00</td><td></td></tr><tr><td>Gestion des sédiments, de la source à la mer (mars 2001)</td><td>75,00</td><td></td></tr><tr><td>Imagerie satellitaire et radar, au service de l'eau (juin 2001)</td><td>75,00</td><td></td></tr><tr><td>Forêts et eau (septembre 2001)</td><td>75,00</td><td></td></tr><tr><td>Rencontres de Grenoble (novembre 2001)</td><td>30,00</td><td></td></tr><tr><td>Variations climatiques et Hydrologie (décembre 2001)</td><td>75,00</td><td></td></tr><tr><td>Hidroenergia 2002 - Conférence sur la Petite Hydroélectricité (juillet 2002)</td><td>100,00</td><td></td></tr><tr><td>Le SAGE, une chance pour la politique de l'eau, (mai 2003)</td><td>55,00</td><td></td></tr><tr><td>Quels modèles physiques pour le 21è siècle ? (octobre 2003)</td><td>60,00</td><td></td></tr><tr><td>Etiages et crues extrêmes régionaux en Europe : perspectives historiques (janvier 2004)</td><td>60,00</td><td></td></tr><tr><td>Crues méditerrannéennes : vulnérabilité, aléa, prévision et gestion de crise (juin 2004)</td><td>60,00</td><td></td></tr><tr><td colspan="3"><strong>COLLECTION : MECANIQUE DES FLUIDES</strong></td></tr><tr><td>3ème journées Cavitation (1996)</td><td>75,00</td><td></td></tr><tr><td>Hydrodynamique des procédés industriels diphasiques (1998)</td><td>60,00</td><td></td></tr><tr><td>Interactions mécaniques entre fluides et structures (novembre 1999)</td><td>60,00</td><td></td></tr><tr><td>Microfluidique – microécoulements liquides et gazeux et leurs applications (décembre 2002)</td><td>90,00</td><td></td></tr><tr><td>Modélisation des écoulements diphasiques (novembre 2004) - CDrom</td><td>40,00</td><td></td></tr><tr><td>Microfluidique – μFlu’04 (décembre 2004) - CDRom</td><td>40,00</td><td></td></tr><tr><td colspan="3"><strong>COLLECTION : JOURNEES DE L'HYDRAULIQUE</strong></td></tr><tr><td>13èmes. Influence des activités de l'homme sur le cycle hydrométéorologique (1974)</td><td>45,00</td><td></td></tr><tr><td>14èmes. La mécanique des Fluides et l'environnement (1976)</td><td>45,00</td><td></td></tr><tr><td>15èmes. L'Hydrotechnique au service d'une politique de l'Eau (1978)</td><td>45,00</td><td></td></tr><tr><td>16èmes. Rejets de chaleur à l'atmosphère (1980)</td><td>45,00</td><td></td></tr><tr><td>17èmes. L'assainissement de demain (1982)</td><td>épuisé</td><td></td></tr><tr><td>18èmes. L'Hydraulique et la maîtrise du Littoral (1984)</td><td>45,00</td><td></td></tr><tr><td>20èmes. Machines hydrauliques (1989)</td><td>45,00</td><td></td></tr><tr><td>21èmes. Les Eaux souterraines et la gestion des eaux (1991)</td><td>45,00</td><td></td></tr><tr><td>22èmes. L'Avenir de l'Eau (1992)</td><td>45,00</td><td></td></tr><tr><td>23èmes. Crues et Inondations (1994)</td><td>épuisé</td><td></td></tr><tr><td>24èmes L'Eau, l'Homme et la Nature (1996)</td><td>épuisé</td><td></td></tr><tr><td>25èmes L'école française de l'eau au service du développement mondial (1998)</td><td>45,00</td><td></td></tr><tr><td>26èmes Les logiciels de mécanique des fluides au service de l'industrie et de l'environnement (2000)</td><td>30,00</td><td></td></tr><tr><td>27èmes Eau et Economie (2002)</td><td>45,00</td><td></td></tr><tr><td>28èmes L'Eau et le Monde vivant (2004)</td><td>45,00</td><td></td></tr><tr><td colspan="3"><strong>COLLECTION SHF - ADEME : Petite Hydroélectricité</strong></td></tr><tr><td>Guide général (1999)</td><td>75,00</td><td></td></tr><tr><td>Intégration dans l'environnement (1998)</td><td>45,00</td><td></td></tr><tr><td>Spécifications techniques propres aux Equipements (1997)</td><td>75,00</td><td></td></tr><tr><td>Actualité d'une énergie renouvelable (actes du colloque, 1998)</td><td>45,00</td><td></td></tr><tr><td colspan="3"><strong>CAHIERS DES CHARGES - TYPES :</strong></td></tr><tr><td colspan="3"><strong>Pour la fourniture de turbines hydrauliques : (édition française ou édition anglaise)</strong></td></tr><tr><td>1-turbines Francis, Pelton et Kaplan, + spécifications propres aux turbines Francis (1972)</td><td>90,00</td><td></td></tr><tr><td>2-spécifications propres aux turbines Pelton (1972)</td><td>38,00</td><td></td></tr><tr><td>3-spécifications propres aux turbines Kaplan (1973)</td><td>38,00</td><td></td></tr><tr><td>4-groupes du type Bulbe (1974)</td><td>90,00</td><td></td></tr><tr><td colspan="3"><strong>Pour la fourniture et le montage des conduites forcées en acier et de leurs accessoires (1972)</strong></td></tr><tr><td>Données techniques annexes (1968)</td><td>90,00</td><td></td></tr><tr><td></td><td>60,00</td><td></td></tr><tr><td colspan="3"><strong>COLLECTION : HISTOIRE ET HOUILLE BLANCHE</strong></td></tr><tr><td>Actes du 2ème congrès de la Houille Blanche (1914)</td><td>150,00</td><td></td></tr><tr><td>Actes du 3ème congrès de la Houille Blanche (1925)</td><td>150,00</td><td></td></tr><tr><td>T.V.A. 5,5 % incluse pour la France (Frais de port en sus)</td><td></td><td>Tarif au 1/01/2005</td></tr></tbody></table>

BON DE COMMANDE :

NOM : ORGANISME

ADRESSE :

MAIL :
Nouvelles de l'AIRH : sommaires du Journal of Hydraulic Research

JHR 2005 VOL 43(2)

* Decision support model for operation of multi-purpose water resources systems, O. R. DOLLING and E A. VARAS

* Elaboration of management tool of a reservoir dam on the Sebou river (Morocco) using an implicit hydraulic model, M. IGOUZAL and A. MASLOUHI

* Sediment transport capacity in rivers, S. YANG

* Energy losses and threshold conditions for choking in channel contractions, B. WU and A. MOLINAS

* Investigation of near wall velocity in 3-D smooth channel flows, S. Yang, SIOW-YONG and J.A. MCCORQUODALE

* Curvilinear turbulence modelling of open channel flow, J.J.R. WILLIAMS

* Correlation and distribution of shear stress for turbulent flow in a smooth rectangular open channel, A. BILGIL

* Floating breakwaters under regular and irregular wave forcing: reflection and transmission characteristics, E. KOUTANDOS, P. PRINOS and X. GIRONELLA

* Fall velocities of natural sediment particles: a simple mathematical presentation of the fall velocity law, K. SHE, L. TRIM and D. POPE

* Inverse modeling of unsaturated heterogeneous media by two-phase pressure derivative matching, K. MASUMOTO

* Multi-level fracture network model and FE solution for ground water flow in rock mass, C. JUNRUI, L. SHOUYI and WUYANQING

JHR ISSUE Vol. 43 (3) 2005

* Du Boys and sediment transport, W.H. HAGER

* Physical modelling of the flow field in an undular tidal bore, H. CHANSON MIAHR

* Influence of bed roughness on sediment suspension: experimental and theoretical studies, B.S. MAZUMDER, K. GHOSHAL and D.C. DALAL

* Velocity measurements in a developed open channel flow in the presence of an upstream perturbation, RAM BALACHANDAR and V.C. PATEL

* Success of booster chlorination for water supply networks with genetic algorithms, OSMAN N. OZDEMIR and M. ERKAN UÇANER

* Turbulence particle models for tracking free surfaces, SONGDONG SHAO and HITOSHI GOTOH

* Three-dimensional numerical simulation of wind-driven flows in closed channels and basins, FRANCESCO CIOFFI, FRANCESCO GALLERANO and ENRICO NAPOLI

* Modelling suspended sediment transport using an integrated numerical and ANNs model, B. LIN and M.M. NAMIN

* Parameter identification using optimization techniques in open-channel inverse problems, HÉLÈNE ROUX and DENIS DARTUS

* Theoretical and experimental study on flow distribution at fracture intersections, HUYUNJIN, MAO GENHAI, CHENG WEIPING and ZHANG JUNJUN

JHR ISSUE Vol. 43 (4) 2005

* Finite volume modeling of variable density shallow-water flow equations for a well-mixed estuary: application to the Río Maipo estuary in central Chile, B.LOOSE, Y. NIÑO C, C. ESCAURIAZA M.

* Two-dimensional simulation of a thermally stratified reservoir with high sediment-laden inflow, D.L. YOUNG, Q.H. LIN, K. MURUGESAN

* SPH simulation of solitary wave interaction with a curtain-type breakwater, SONGDONG SHAO

* High-resolution TVD schemes in finite volume method for hydraulic shock wave modelling, GWO-FONG LIN, JIHN-SUNG LAI, WEN-DAR GUO

* Two-dimensional and line jets in a weak cross-flow, J.F. HUANG, M.J. DAVIDSON, R.I. NOKES

* On incompressible turbulent flow: partial average based theory and applications, G. GAO and Y. YONG

* Vertically integrated numerical model to simulate breach formation by flow overtopping earth-filled dams: hydraulics, R. FRENETTE, A. MUNTEANU

* Discharge equations for venturimeter and orificemeter, PRABHATA K. SWAMEE

* Modified log-wake law for zero-pressure-gradient turbulent boundary layers, JUNKE GUO, PIERRE Y. JULIEN, ROBERT N. MERONEY