GPS IN ANTARCTICA

extracts from Marie-Noelle Bouin PhD. Thesis, 1999 ©







Antarctic French base GPS permanent station (DUM1)









1- Antarctic Ice Sheet Melting and Post Glacial Rebound

The present day behaviour of the antarctic ice cap is a matter of great importance since it is related to, among other things, the short term evolution of the relative sea level. Is the antarctic ice cap melting, as suggested by the recent observation of gigantic icebergs leaving the peninsula ? Is it on the contrary going o thicken in the next years due to increase of precipitation generated by global temperature rise ?
Paradoxically, one don't get direct measurements of ice cap thickness variations, precise enough and widely distributed enough to answer this question. But, since changes in the mass of the ice cap do influence the crustal surface position underneath, informations on ice cap evolution can be inffered from measurements of crustal position at given sites. The recent developpment of GPS (Global Positioning System), a spatial geodetic technique with centimetric position accuracy, enables us to actually do those measurements.
The motions of chosen sites at the surface of the Antarctic crust are not only a response to present day ice unloading, but also to the full history of ice thickening or thining over the past thousands of year since the last galcial maximum 18 000 years agoo. The late effects come from the visous mantle (Post-Glacial-Rebound), when the instantaneous effects come from it's elastic component. Therefore, measuring present day crustal motions in Antarctica, give information not only on the last deglaciation but also on the crustal and mantle parameters involved in the process.

2- Geodetical Measurements of Motions and Deformations

All those motions (viscous or elastic) have very small amplitude, usually less then a centimeter per year. Nevertheless, using the permanent GPS stations operated at the permanently occupied antarctic bases on the continent, it is possible to detect such small motions after a few years of measurements.

Tectonic Plates in Southern Hemisphere The Network of Permanent GPS Sites in Antarctica


We choose to use data from permanent IGS (International GPS System) sites, because they cover the all continent since 1995, although they are only a few of them. The situation considerably improved from 1995 (only 4 stations) to 1999 (already 9 stations). We contributed to this network by installing a GPS permanent station at the french base of Dumont d'Urville (Terre Adélie) in December 1997. This was done in the framework of the french antarctic programm, funded and coordinated by the french institute for polar research (IFRTP), in collaboration with the french institute for geography (IGN).



3- Different Expected Motions

The expected displacements are thus made of a superposition of different effects with different time constants and spatial distribution.


4- Description of Geodetic Solution

The projection of our measurements in the global reference system (ITRF) was conducted using different techniques : CATREF software (used for combining different data sets to realize the ITRF), and GLOBK (Kalman filtering of epoch measurements), after the initial GPS measurements processing done with the GAMIT (GPS At MIT) software. We used those different techniques to assess the level of coherence and internal precision of the data set. The ultimate solution is produced by combining daily coordinates produced by free network solutions, by the mean of daily 7 parameters (Helmert) transformations to the reference frame. We used both ITRF96 and ITRF97 and produce comparisons. In order to constrain this refernce frame we were conducted to use a certain number of additional stations, outside the Antarctic plate itself.


The network of permanent GPS stations used in our computations


5- Obtained Motions and Velocities

5a- Temporal Series at Dumont d'Urville (Terre Adélie)
The carefull analysis of temporal series allow to obtain horizontal and vertical velocities for stations with at least 3 years of data. Some of them exhibit non-linear behaviour for different reasons, either the state of the antenna and/or the mounting device including the monument or true ground motion.
The temporal series for Dumont d'Urville Show the effect of the Magnitude Mw=8.1 earthquake, which took place under the Balleny Islands (600 kms away), march 25th, 1998. At this date, the position of the station was ofset by approximately 1 cm. This displacement is slightely less than the modelled co-seismic displacement using Okada's formule in elastic homogeneous half-infinite space, but quite coherent in direction. Of course, it is the velocity corrected for this "step displacement" which is used to quantify the solid rotation pole of the Antarctic plate.

5b- Horizontal Velocities and Antarctic Plate Rotation Pole
All the stations located on the Antarctic plate exhibit very coherent velocities very much in agreement with the rotation of a rigid block : the undeformed Antarctic plate. We estimate the pole of rotation located by 62° North latitude and 147° West longitude, with an angular rotation speed of 0.26 °/Myr.
On the figure, the red arrows stand for our solution, when the green arrows depict the velocities corresponding to a rigid rotation around a given pole (green dot circled in red), and the blue arrows depict the residual velocities : the difference between the measured velocities and the predicted velocities based on the rigid rotation. One realize the very good agreement with the model by looking at the blue arrows (residual) which are below a couple of mm/yr everywhere. Therefore, it can be stated that the Antarctic plate is rigid at this level (2mm/yr) of precision. It should be noted that O'higgins at the tip end of the peninsula is an exeption.
One can also realize that the estimated pole, differs from the Nuvel-1a (brown dot on the figure). An independant computation (Heflin et al., 1999) using Australian data also suggest a new rotation pole for Australia. The absolute poles for those 2 plates differ significantly from the Nuvel-1a estimates, but the relative motion between Antarctica and Australia remains unchanged whether one looks at the Nuvel-1a pair of poles or our new estimates. Since this relative motion is the only strong constraint (it is given by the magnetic anomalies and the transform fault azimuth at the ridge between the two plates), we iestimate this fact is a good validation of our results.

5c- Comparison with DORIS System
The french positioning system DORIS also give displacement velocities for antarctic stations. (Crétaux et al., 1998). In this case also a pole of rotation can be determined, and residuals with respect to this rigid rotation. The DOIS pole is almost identical to the Nuvel-1a pole, but residuals are more important and reach 5 mm/yr on 3 stations. This may indicate that DORIS velocities are either less precise than GPS velocities or affected by local "errors" like unstable benchmarks. The fact that One finds a pole similar to Nuvel-1a indicates that DORIS can produce a reasonnably good reference system, independant from GPS.


Temporal series at Dumont
The cross (resp. the circles) stand for positions before (resp. after) the earthquake.




horizontal velocities and position of rotation pole and residuals with respect to this rigid rotation



DORIS horizontal velocities, and the estimated pole, and residual velocities


6- Deglaciation and Instantaneous Elastic Rebound on the Peninsula

La seule station où l'on obtient une vitesse résiduelle horizontale non nulle après retrait du mouvement correspondant à la rotation rigide de la plaque tectonique antarctique et l'influence du séisme des Iles Balleny, est celle de O'Higgins située à l'extrémité de la Péninsule. La très bonne cohérence de ce résidu d'une solution à l'autre, et son amplitude de 8 mm/an significative par rapport aux incertitudes nous ont orientés vers une interprétation glaciologique. La Péninsule est la région la plus chaude de tout le continent antarctique, soit la meilleure candidate à une débâcle rapide. Cette vitesse horizontale peut être associée à un retrait élastique instantané vers un centre de décharge glaciaire situé près de la base de la Péninsule (plate-formes de Larsen ou de Filchner-Ronne, sur lesquelles les taux de vêlage sont élevés actuellement). En effet, la croute continentale réagit comme une éponge élastique qui se regonfle (mouvement vertical vers le haut) et qui se rétracte (mouvement horizontal en arrière) quand on arrete d'appuyer dessus (la masse de glace qui diminue). L'amplitude des mouvements horizontaux d'origine flexurale est d'environ 1/3 de celle des mouvements verticaux de rebond élastique correspondants. Ce résidu horizontal de 8 mm/an serait donc cohérente avec le mouvement vertical positif de 25 mm/an que l'on observe sur cette station de O'Higgins.


7- Vertical Velocities : Towards a general uplift, more important than expected ?

Parmi les vitesses verticales sur le continent antarctique, la plus élevée, dans le sens d'une surrection importante, est celle observée sur cette même station de O'Higgins, avec 25 +/- 12 mm/an. Ce mouvement est trop important, comparé aux autres vitesses verticales des stations antarctiques, pour ne traduire que du rebond visqueux. Il indique certainement une remontée élastique d'au moins 10 à 15 mm/an, correspondant à une évacuation glaciaire actuelle importante au voisinage de la Péninsule. L'application d'une modélisation simple à partir de disques représentant la charge de glace, permet d'évaluer la diminution de l'épaisseur de glace nécessaire pour provoquer un rebond vertical de 1 cm/an à O'Higgins. Si cette diminution est uniforme sur l'ensemble de la Péninsule, il faudrait qu'elle atteigne 65 à 70 cm/an, chiffre important par rapport aux évaluations courantes des modèles, mais qui reste possible, et même comparable à des valeurs locales (par exemple celles du scénario J92 de (Jacobs et al., 1992). Un tel mouvement constitue une indication importante et nouvelle d'une accélération sensible de l'évacuation glaciaire sur la Péninsule.
Paradoxalement, il est difficile d'obtenir des mesures directes fiables de l'ablation des plate-formes, qui soient représentatives du comportement glaciaire sur l'ensemble de la Péninsule. On y a cependant constaté depuis le début des années 90 une accélération du vêlage, avec la disparition rapide de plusieurs grandes plates-formes glaciaires, mais les glaciologues hésitent à en déduire une tendance à long terme.
Cette vitesse de remontée importante à O'Higgins vient donc confirmer des observations glaciologiques. Les études de Rignot, 1998b montrent également des vitesses de retrait des lignes d'échouage sur les plates-formes de la Péninsule beaucoup plus rapides que les résultats des modèles, qui servent eux-mêmes de base aux bilans de l'équilibre actuel de la calotte antarctique. On ne dispose pas encore de mesures globales de vitesse de remontée de la croûte d'origine géodésique venant confirmer les résultats de cette étude. Localement, a la station japonaise de Syowa (Kamimuna, 1998), les résultats de plus de 20 ans de mesures de nivellement et de marégraphie, relayées par du GPS, montrent des vitesses de remontée continue de 10 mm/an environ.
Notre résultat constitue donc une première validation géodésique de ces observations glaciologiques. Le fait de ne disposer que d'une vitesse GPS sur la Péninsule empêche des interprétations plus détaillées. L'installation mi-1998 d'une seconde station GPS permanente à la station de Palmer devrait permettre d'obtenir d'ici quelques années une autre vitesse verticale sur la côte Ouest de la Péninsule. Par la suite, il faudrait disposer de données plus nombreuses, couvrant si possible l'ensemble de la Péninsule, qui devraient fournir une cartographie plus complète du mouvement de rebond élastique, et donc des zones principales d'évacuation glaiciaire. La fiabilité des vitesses verticales GPS dépend en grande partie de la façon dont les données ont été obtenues.
Les variations saisonnières très importantes, détectées sur la composante verticale de O'Higgins, confirment la présence d'un mouvement de rebond élastique. Leur période est de 365 jours exactement, avec une amplitude de 5 cm. Le mécanisme (encore incertain) justifiant des variations d'une telle ampleur devrait faire appel à des variations très importantes de l'évacuation glaciaire sur la Péninsule entre l'été et l'hiver.
Les vitesses verticales GPS sur le reste des stations antarctiques, bien que reflétant un mouvement moins important, et avec des incertitudes du même ordre de grandeur que les vitesses elles-mêmes, indiquent une remontée globale du continent. Un rebond de 10 mm/an environ se situe plutôt vers la limite supérieure des estimations des modèles de rebond post-glaciaires, montrant ainsi que la relaxation visqueuse provoquée par la dernière déglaciation sur l'Antarctique est encore effective. Le fait de ne pas détecter de divergence horizontale claire accompagnant ce mouvement vertical peut être attribué à l'incertitude sur les vitesses résiduelles horizontales par rapport à la rotation de plaque tectonique rigide. Une vitesse inférieure à 2 mm/an, conforme aux prédictions des modèles ICE-3G, LC79 ou D91, passerait inaperçue parmi les résidus horizontaux du même ordre de grandeur.

Les vitesses verticales mesurées en Antarctique :
  • en rouge les vitesses des stations sur la plaque Antarctique
  • en vert les vitesses des stations sur les autres plaques.


vitesses en plan : les flèches pointant toutes vers l'extérieur indiquent une surrection générale




vitesses en coupe : altitude en fonction de la latitude. les barres grises verticales indiquent l'incertitude sur les vitesses



8- Conclusion

Nous avons montré dans cette étude que le traitement de 4 années de données GPS continues sur l'Antarctique et ses environs, avec une méthode adaptée, permettait d'obtenir des vitesses exploitables, sur la composante verticale comme horizontale.

Les enjeux de ce dernier point sont particulièrement importants, puisque la fonte totale de la partie Ouest de la calotte implique une augmentation du niveau des mers de plus de 6m. L'année 1998 a été celle du GPS en Antarctique, puisque 4 nouvelles stations permanentes fournissent désormais des données accessibles. L'intérêt de la géodésie dans la détection des variations glaciaires n'a jamais été aussi évident, c'est à partir de maintenant que cette nouvelle direction de recherche prend tout son sens.




© all results displayed here are extracted from M.N. Bouin PhD. thesis, conducted at "laboratoire de géologie de l'ENS", under co-direction of Ch. Vigny (ENS) et C. Boucher (IGN), defended October 14th, 1999, in front of the jury composed by : F. Barlier (président), K. Lambeck (reviewer), M. Kasser (reviewer), F. Remy (examinor), B. Ambrosius (examinor), C. Boucher, et Ch. Vigny.





Retour au sommaire