Afin d’étudier les variations climatiques enregistrées par les carottes de glace, il est nécessaire d’avoir des datations précises à la fois pour les phases gaz et glace. Le but de ce travail de thèse a été d’améliorer les chronologies des carottes de glace, couvrant les derniers 800 000 ans, au travers de nouvelles mesures de la composition isotopique de l’air (δ15N, δ18Oatm et δO2/N2) piégé dans la glace d’EPICA Dôme C (EDC) et de l’utilisation de l’outil de datation "Datice". Le premier résultat important de cette thèse a été la production de la chronologie cohérente pour les carottes de glace ("Antarctic Ice Core Chronology", AICC2012) pour EDC, Vostok, EPICA Droning Maud Land (EDML), TALos Dôme ICE core (TALDICE) et NorthGRIP. Sur cette nouvelle chronologie la théorie du see-saw bipolaire reste valable. AICC2012 donne un âge pour la Terminaison II en accord avec les autres archives climatiques. De plus, les durées des périodes interglaciaires restent inchangées par rapport à la chronologie EDC3. Lors de la construction d’AICC2012 nous avons mis en évidence plusieurs points nécessitant des améliorations. Nous avons donc procédé à l’amélioration de Datice dans le but d’intégrer correctement les contraintes issues du comptage des couches annuelles et leurs erreurs. Ces améliorations conduisent à des chronologies cohérentes tout en respectant les hypothèses sous-jacentes à Datice. Nous proposons aussi une nouvelle formulation de l’erreur associée à la fonction d’amincissement à partir de l’analyse des propriétés mécaniques de la glace dans le cas d’EDC. Pour finir, les nouvelles mesures du δO2/N2 et du δ18Oatm effectuées sur de la glace bien conservée d’EDC nous ont permis de définir de nouvelles contraintes d’âge. La comparaison de ces traceurs mesurés à Vostok, EDC et Dôme F sur le MIS 5 a permis de mettre en évidence une possible influence de paramètres climatiques locaux sur le δO2/N2. L’analyse du retard entre le δ18Oatm et la précession sur les derniers 800 ka montre des variations de ce dernier. Nous suggérons que ce retard est augmenté par l’occurence d’évènements de Heinrich à certaines périodes. Les résultats de cette thèse sont à prendre en compte pour le prochain exercice de datation cohérente pour les carottes de glace. / In order to study the climate variations recorded by ice cores, it is necessary to have precise chronologies for the ice and gas phases. The aim of this work has been to improve ice cores chronologies, covering the last 800 000 years, through new measurements of the isotopic composition of the air δ15N, δ18Oatm et δO2/N2) trapped in EPICA Dome C (EDC) ice core and the use of the Datice dating tool.The first important result of this PhD has been the production of the Antarctic Ice Core Chronology (AICC2012), common for EDC, Vostok, EPICA Droning Maud Land (EDML), TALos Dome ICE core (TALDICE) and NorthGRIP ice cores. The bipolar see-saw theory is still valid on the new chronology. The AICC2012 chronology gives an age for Termination II in good agreement with other climate archives. Moreover, the duration of interglacial periods is unchanged compared to EDC3. While building the AICC2012 chronology, we have pointed out several limitations. Since then, we have improved Datice in order to correctly integrate constraints deduced from layer counting and their associated uncertainties. These improvements permit to build coherent chronologies respecting the underlying hypotheses of Datice. Moreover, we propose a new parameterization of the uncertainty associated with the background thinning function based on ice mechanical properties of EDC ice core. Finally, we were able to deduce new age constraints thanks to the new measurements of δO2/N2 and δ18Oatm performed on well-conserved ice from EDC. A multi-proxy comparison of Vostok, EDC and Dome F ice cores over MIS 5 has highlighted a possible influence of local climatic parameters on δO2/N2. The analysis of the delay between δ18Oatm and precession shows some variability over the last 800 ka. We propose that the delay between δ18Oatm and precession is increased during periods associated with Heinrich events. The results obtained during this PhD should be used for the next ice core coherent chronology.

Les glaces polaires, en plus d'améliorer notre compréhension du climat, donnent accès aux retombées de béryllium-10 (10Be), isotope d'origine cosmogénique formé par l'interaction du rayonnement cosmique et des hautes couches de l'atmosphère. Son taux de production dépend de l'intensité du rayonnement cosmique primaire, qui est modulé par l'activité solaire et le champ magnétique terrestre. Le 10Be contient donc des informations sur ces deux paramètres. Il a permis, entre autres, l'amélioration des chronologies des carottes de glace grâce à des marqueurs stratigraphiques absolus liés à des évènements remarquables du champ géomagnétique comme l'excursion de Laschamp ou l'inversion de Brunhes-Matuyama.EPICA Dome C (75° 06' S, 123° 21' E) est une carotte de glace de 3270 m de long forée en Antarctique Est. Elle offre un enregistrement climatique complet durant les 800 000 dernières années (kyr BP). Dans le cadre de cette thèse, 2200 échantillons de 10Be ont été mesurés entre 2384 m (269 kyr BP) et 2627 m (355 kyr BP) de profondeur. Cette séquence continue offre, d'une part, la possibilité d'étudier l'activité solaire durant la période interglaciaire du Stade Isotopique Marin (MIS) 9.3, où la résolution atteint ~20 ans. Nos résultats contrastent avec ceux durant l'Holocène, avec l'absence remarquée du cycle de de Vries (210 ans) dans notre profil de 10Be. D'autre part, cette séquence nous a permis de vérifier que l'estimation classique de l'accumulation de neige obtenue par les modèles de datation de carottes de glace est correcte à 20% près lors de la succession des cycles glaciaires – interglaciaires.Ces mesures ont été combinées avec celles précédemment effectuées sous la responsabilité de Grant Raisbeck (publication en préparation). Ceci permet de disposer d'un profil continu et détaillé (résolution de 20 à 250 ans) entre 200 et 800 kyr BP sur la carotte de glace EPICA Dome C. Les variations de flux de 10Be observées se comparent bien aux changements d'intensité du champ géomagnétique enregistrés dans les sédiments marins. À partir de cette observation, il est possible de proposer une synchronisation de ces profils afin de les placer sur une échelle d'âge commune. Les déphasages observés entre l'augmentation de température en Antarctique (augmentation du D à EDC) et la hausse du niveau marin global (baisse du 18O marin des archives sédimentaires) n'excèdent pas 3200 ans excepté à la Terminaison VII. La faible résolution des enregistrements de 18O dans les sédiments marins rend l'analyse des déphasages délicate. Les résultats obtenus encouragent à analyser des périodes plus récentes autour de la Terminaison II, incluant les excursions géomagnétiques de Blake et d'Iceland Basin. Ceci permettrait en effet de contraindre plus fortement la synchronisation entre le signal paléomagnétique dans les sédiments marins et le flux de 10Be à EDC. / The polar ice cores, in addition to improving our understanding of the climate, give access to beryllium-10 (10Be) fallout, an isotope of cosmogenic origin created by the interaction of Galactic Cosmic Rays with the upper atmosphere. Its production rate depends on the intensity of the primary cosmic rays, which are modulated by solar activity and the Earth's magnetic field. 10Be therefore provides information on these two parameters. It has allowed, among others, the improvement of ice cores chronologies thanks to absolute stratigraphic markers linked to remarkable events of the geomagnetic field as the Laschamp excursion or the Brunhes-Matuyama reversal.EPICA Dome C (75° 06' S, 123° 21' E) is a 3270~m ice core drilled in East Antarctica. It offers a complete climate record over the last 800 000 years (kyr BP). In the framework of the PhD, 2200 10Be samples were measured between 2384 m (269 kyr BP) and 2627 m (355 kyr BP) deep. This continuous sequence provides, on the one hand, the opportunity to study the solar activity during the interglacial period of the Marine Isotope Stage (MIS) 9.3, where the resolution reaches ~20 years. Our results contrast with those during the Holocene, with the noted absence of the de Vries cycle (210 years) in our 10Be profile. On the other hand, this sequence allowed us to verify that the classical estimate of snow accumulation obtained by the dating models of ice cores is correct to within 20% during the succession of glacial – interglacial cycles.These measurements were combined with those previously performed under the responsibility of Grant Raisbeck (publication in preparation). This allows to dispose of a continuous and detailed profile (resolution: 20 to 250 years) between 200 and 800 kyr BP on the EPICA Dome C ice core. The observed 10Be flux variations compare well with changes in the intensity of the geomagnetic field recorded in marine sediments. From this observation, it is possible to propose a synchronization of these profiles in order to place them on a common age scale. The observed phase shifts between the rise of temperature in Antarctica (increase of D at EDC) and the rise of global sea level (decrease of marine 18O from marine sediments) do not exceed 3200, except at the Termination VII. The low resolution of the 18O records from marine sediments makes it difficult to analyze the phase shifts. The obtained results encourage to analyze more recent periods around Termination II, including the Blake and Iceland Basin geomagnetic excursions. Indeed, this would allow to constrain more strongly the synchronization between the paleomagnetic signal in marine sediments and the 10Be flux at EDC.

L'Antarctique est encore un continent à explorer pour comprendre le climat terrestre, passé, présent et futur. Grâce à l'analyse des carottages, on peut remonter l'histoire du climat terrestre sur plusieurs centaines de milliers d'années. Grâce à l'étude de sa topographie, on peut remonter à la variabilité de son climat, connaître sa dynamique : écoulement de la glace, vêlages d'icebergs. . . Grâce au suivi de la topographie, c'est-à-dire la variation de la hauteur de glace, on peut estimer les pertes et les gains de masses ainsi que la contribution à la hausse du niveau de la mer, pour déduire l'impact sur les océans et les littoraux. Pour explorer ce continent, nous utilisons dans cette thèse l'instrument embarqué sur satellite qu'est l'altimètre. Une onde radar ou laser est envoyée sur le sol et la réflexion de ce signal est enregistrée permettant d'extraire des paramètres pertinents et la hauteur de glace. Les satellites ont pris le relais des expéditions célèbres du XIXème siècle, notamment depuis le lancement du premier altimètre couvrant jusqu'à 82æ S ce continent, embarqué sur la mission ERS-1. Depuis 2013, SARAL marque un changement dans les caractéristiques instrumentales par rapport aux missions précédentes, car elle émet dans une fréquence appartenant à la bande Ka (36.75 Ghz) au lieu des bandes Ku (13.6 Ghz) et S (3.2 Ghz) ce qui implique une interaction avec la surface différente qu'il faut pouvoir comprendre. Avec près de trois ans d'observations, ce manuscrit regroupe les travaux sur le traitement des données altimétriques dans un but de valorisation de la mission SARAL. Nous présentons le contexte géographique et le contexte technique pour montrer les difficultés du traitement des données . En utilisant la méthode des points de croisement que nous expliciterons, nous présentons l'importance de la comparaison de deux missions concomitantes, ENVISAT et ICESat afin d'estimer la profondeur de pénétration en bande Ku de l'onde radar et corriger de cet effet, étude que l'on pourra répéter avec SARAL une fois la mission ICESat-2 lancée. Nous avons adapté une chaîne de calibration et de validation des données sur glace continentale qui nous permet d'établir des diagnostics, des statistiques et des suivis temporels. La pente de la surface joue un rôle prépondérant dans la dégradation de la précision des données. Nous mettons en place un processus de sélection des données altimétriques afin d'augmenter cette précision. Grâce à cette sélection et l'application de la correction de marée océanique, nous quantifions avec des métriques aux points de croisement la performance de la mesure de SARAL sur glace continentale, et nous réduisons l'imprécision de la mesure de plus de 90%. Enfin, grâce aux sorties de notre chaîne de validation, nous décrivons les observations de SARAL sur près de trois ans afin d'apporter une analyse préliminaire cohérente avec la mission précédente qu'est ENVISAT. Nous soulevons les améliorations à apporter pour aller plus loin dans la physique de la mesure. Que ce soit la comparaison avec d'autres missions, l'utilisation de modèles ou des modifications dans l'algorithme d'extraction des paramètres du signal appelé retracking, SARAL apporte de nouvelles perspectives afin que l'évolution de la calotte Antarctique soit connue avec une précision toujours croissante. / Antarctica still remains a fascinating place to be explored. With ice corings, the past Earth climate can be retraced. Studying its surface, its climate variability and its dynamic are better known : ice velocity, iceberg calvings. . . Through the height variations monitoring, ice gains and losses are estimated, leading to the contribution to the sea-level rise, from which is deduced the impact over coastal areas or the oceans. To explore this continent, we use the sensor called altimeter : a radar or laser wave is sent from the satellite to the surface and the reflected signal is recorded. From this signal we extract relevant parameters and the height. The satellites have observed this area, notably since the altimetric mission ERS-1 launching, until 82æ S in 1991, following the still famous explorations from the beginning of the twentieth century. SARAL, launched in February 2013, innovates because of the major change in the frequency used, the Ka-band (36.75 Ghz) instead of the Ku-band (13.6Ghz) and the S-band (3.2 Ghz) implying a different interaction between the radar wave and the snowpack that needs to be investigated. Using almost three years of observations, we focus on the altimetric signal processing and its validation. The geographic and the technical aspects are introduced and we show the limitations in the altimetric data processing. Using the crossover method (explained in greater detail later), we compare two simultaneous missions, ENVISAT and ICESat in order to better constrain the penetration effect of the radar-wave into the snowpack and correct it. The method that is used will be of great interest to compare SARAL with the future mission ICESat- 2. A new calibration and validation tool has been implemented, allowing a long-term survey of the Antarctic area, providing statistics, diagnostics and temporal series. The slope effect is the major limitation in the precision assessment of the mission. Thanks to a new way of selecting the data and the oceanic tide correction, we quantify with metrics computed at crossover points the accuracy of altimetric data over the Antarctic ice sheet. We describe in the last part the SARAL observations to give a preliminary analysis in agreement with the former mission ENVISAT. Other ways of improvement are presented, like future altimetric missions, modifications in the so-called retracking algorithm that extracts the relevant parameters or even the use of models. SARAL raises perspectives to estimate with a growing precision the evolution of the Antarctic continent, and this thesis is a state of the art about the different possible processings needed to do so.

Dans le contexte actuel du réchauffement climatique, l'une des principales sources d'incertitude pour l'élévation du niveau de la mer est la contribution de la calotte Antarctique. L'étendue et les conditions météorologiques extrêmes de ce continent font de la télédétection spatiale un moyen utile pour son suivi sur le long terme. Les observations satellites altimétriques et radiométriques dans la gamme des micro-ondes rendent compte de l'évolution des propriétés du manteau neigeux de la calotte. L'altimétrie radar, par des mesures répétées de l'élévation de la topographie de surface, permet de quantifier les variations de volume sur l'ensemble du continent. Cependant, la pénétration de l'onde radar dans la neige affecte négativement cette quantification. Les méthodes proposées pour minimiser les erreurs de pénétration sont toutes basées sur des relations avec le coefficient de rétrodiffusion radar. La compréhension des variations annuelles et inter-annuelles du coefficient de rétrodiffusion est nécessaire pour améliorer la précision de l'estimation de l'élévation de la surface donc du bilan de volume de la calotte. Cette thèse a pour objectif d'étudier le coefficient de rétrodiffusion mesuré par les altimètres sur l'ensemble du continent, sujet qui jusqu'à aujourd'hui a reçu peu d'attention. Les altimètres radars embarqués à bord d'ENVISAT (bandes S et Ku) et de SARAL/AltiKa (bande Ka) ont des sensibilités différentes aux propriétés de la neige. Nous nous sommes intéressés aux caractéristiques annuelles et inter-annuelles des coefficients de rétrodiffusion dans ces trois bandes. Une étude de sensibilité a été réalisée avec un modèle électromagnétique afin de déterminer les propriétés du manteau neigeux qui dominent le signal saisonnier. On montre que le signal saisonnier est sensible à la densité et la rugosité de surface dans la bande S, à la température de la neige dans la bande Ka et à l'une ou à l'autre de ces variables selon la région dans la bande Ku. Les caractéristiques saisonnières du coefficient de rétrodiffusion sont ensuite comparées à celles des températures de brillance acquises par les radiomètres à bord de SARAL et de SSM/I. Les résultats indiquent une influence significative de la rugosité de surface sur les températures de brillance de la bande Ka, influence souvent considérée négligeable dans la modélisation de la température de brillance. Cette étude apporte une meilleure connaissance de la dynamique saisonnière des propriétés de proche surface de la calotte Antarctique. Elle fournit de nouveaux indices pour développer dans le futur des algorithmes robustes de correction de l'erreur de pénétration. Elle met également en lumière l'importance des missions altimétriques multi-fréquences et les possibilités qu'offrent le signal de la bande S pour l'étude des variabilités saisonnières de la rugosité de surface. En définitive, la rugosité de surface est un paramètre important à prendre en compte pour obtenir de meilleures estimations et modélisations des coefficients de rétrodiffusion et des températures de brillance. / In the context of global climate changes, the Antarctic ice sheet contribution to sea-level rise is one of the main uncertainty sources. The extent and extreme meteorological conditions of this continent render remote sensing a useful tool for long term monitoring. Altimetry and radiometry observations in the microwave range reveal variations of the volume of the ice sheet and surface properties of the snowpack. Radar altimeters, provide repeated observations of the surface topography elevation, which allow the quantification of volume variations of the ice sheet. However, the penetration of radar waves in dry and cold snowpack adversely affects the estimated surface elevation. Approaches to minimize the penetration error are all based on a relationship with the backscattering coefficient. Understanding the annual and interannual variations of the backscattering coefficient is thus a key issue in order to improve the estimation accuracy of the surface elevation and to refine the ice-sheet volume trend. This thesis aims at studying the backscattering coefficients acquired by radar altimeters, which until now have received little attention. Radar altimeters on board ENVISAT (S and Ku bands) and SARAL/AltiKa (Ka band) have different sensitivities to the snowpack properties. The annual and interannual variations of the backscattering coefficient at the three bands is investigated. Sensitivity tests are carried out with an electromagnetic model to determine the prevailing snowpack properties that drive the signal. The seasonal signal is sensitive to surface density and roughness at S band, to snow temperature at Ka band and to either snow surface density and roughness or temperature depending on the location on the continent at Ku band. The seasonal signal of the backscattering coefficient is then compared with that of the brightness temperature measured by radiometers on SARAL and SSM/I. The results show a significant influence of surface roughness on brightness temperatures at Ka band, which has often been neglected in brightness temperature modeling studies. This thesis provides a better understanding of the seasonal dynamics of the near surface properties of the Antarctic ice sheet. It also provides new clues to build a more robust corrections of the penetration errors in the future. It highlights the importance of multi-frequency altimetry missions and the potential of the S band to study the seasonal variability in surface roughness. In summary, surface roughness is an important property which should be taken into account for a better modeling of backscattering coefficient and brightness temperature.

La concentration en isotopes de l'eau dans les glaces polaires et les sédiments marins est un précieux indicateur des températures passées si l'on comprend bien les échanges isotopiques entre l'océan, l'atmosphère et la cryosphère. La dynamique et la capacité de stockage des calottes glaciaires demeurent cependant méconnues à cause de la rareté des prélèvements de glace. Cette thèse, basée sur la modélisation de la glace en 3D, développe une méthode de traçage pour suivre les particules de glace afin de retrouver leurs conditions de déposition et de prédire la fine stratigraphie des calottes du Groenland et de l'Antarctique. La méthode est validée en comparant la prédiction du modèle avec les couches de glace aux sites de forage polaire profonds. En utilisant le modèle d'évolution du Groenland de l'UBC (Marshall et Clarke 1997) et les traceurs aux sites de GRIP, Dye 3 et Camp Century, j'étudie le climat lors de la précédente période inter-glaciaire, il y a 127000 ans, et trouve que la fonte de glace au Groenland provoqua une montée de 3.5 à 4.5 m du niveau des mers. La méthode est aussi appliquée aux conditions de déposition de la glace prélevée en Antarctique à Vostok, Dome Fuji et Dome C en combinant les traceurs au modèle de calotte du LGGE (Ritz et al. 2001). Je teste ainsi les hypothèses d'interprétation des enregistrements climatiques et examine les changements d'altitude et de morphologie de la calotte antarctique. Ayant prédit la fine stratigraphie des calottes polaires, j'en déduis la première estimation glaciologique de leur composition isotopique moyenne pour le présent et le passé et boucle ainsi le budget en isotope de l'eau de notre planète.

Ce travail traite de la simulation numérique du climat des grandes calottes de glace, en particulier des calottes de l'Antarctique et du Groenland, toujours existantes, dans des conditions climatiques différentes, à l'aide de modèles de circulation générale de l'atmosphère (MCGA). Le MCGA à grille variable LMDz a été adapté aux spécificités du climat polaire et validé pour le climat actuel. L'approche d'une grille variable, qui permet d'utiliser le MCGA à haute résolution spatiale (autour de 100 km) sur la région d'intérêt à un coût numérique raisonnable, a été validée en analysant la dynamique atmosphérique au bord de la région ciblée à l'aide d'un schéma de suivi des cyclones individuels. Des simulations du climat du Dernier Maximum Glaciaire (DMG) ont été faites pour le Groenland et l'Antarctique et analysées en tenant compte des archives glaciaires disponibles. Une explication possible des différences entre les deux méthodes principales de reconstruction des paléotempératures - l'analyse des isotopes de l'eau et la mesure directe de la température de la glace dans le trou de forage - au centre du Groenland a pu être proposée. Cette explication est basée sur des changements de paramètres climatiques locaux. C'est la première fois que l'approche de grille variable a été utilisée dans un MCGA pour des simulations du climat polaire à l'échelle de quelques années. Les simulations paléoclimatiques faites avec LMDz sont à une résolution spatiale inégalée à ce jour. Finalement, le climat du DMG, simulé par plusieurs MCGA dans le cadre du projet international PMIP (Paleoclimate Modelling Intercomparison Programme), a été analysé, et des implications des résultats pour l'interprétation des enregistrements glaciaires ont été discutées.

Le système climatique terrestre est forme de plusieurs composantes qui interagissent étroitement (océans, atmosphère, calottes de glace, biosphère et lithosphère). La modélisation numérique permet d' étudier ces interactions. Dans ce travail, nous avons modélisé l' évolution de la calotte de glace Antarctique, le plus grand réservoir d'eau douce sur Terre. Nous avons utilisé le modèle 3D développé au LGGE. Il tient compte du couplage thermomécanique entre les vitesses et les températures dans la glace et de l' écoulement spécifique de la glace flottante. Le modèle a été amélioré sur plusieurs points : la loi de déformation de la glace, la loi de frottement sous les fleuves de glace, le calcul de la chaleur de déformation et la paramétrisation de la fusion sous les plates-formes de glace flottante. De plus, des modifications numériques ont été apportées. Différentes expériences ont montré une grande asymétrie de comportement de la calotte en englacement et déglacement. Une simulation des quatre derniers cycles glaciaires-interglaciaires a permis de valider les modifications effectuées sur le modèle. Les résultats indiquent que, dans le futur, la contribution de l'Antarctique à l'augmentation du niveau des mers pour les prochains siècles devrait être négligeable, puis s'accentuer ensuite. Nous avons aussi montré la forte sensibilité des résultats à la carte de socle utilisée en entrée du modèle. Enfin, la résolution du modèle a été augment ée en passant d'une maille de 40 km à 20 km, ce qui permet de mieux prendre en compte les structures fines de la topographie.

L'objectif principal du travail présenté ici est d'une part la détermination de la relation entre la température locale de l'air et la composition isotopique de la neige à la station de Vostok (Antarctique de l'Est) et d'autre part la reconstruction de variabilité climatique dans cette région au cours des 200 dernières années. On a produit des données expérimentales concernant la composition isotopique et l'accumulation de la neige à partir de 8 puits creusés à Vostok dont 2 couvrent les 200 dernières années et qui ont été réalisés pendant les saisons d'été 1998/99, 1999/2000 et 2001/02. Le régime météorologique est documenté par les observations instrumentales et les radiosondages de l'atmosphère (depuis Décembre 1957) archivés à l'Institut Recherche Arctique et Antarctique (St. Pétersbourg, Russie). Nous avons aussi utilisé les données d'accumulation de la neige venant du réseau des balises installé à proximité de la station de Vostok en 1970, documenté la distribution spatiale de la teneur isotopique de la neige de surface, mesuré le nivellement de la surface, et étudié la variation de la composition isotopique des précipitations et celle de la neige transportée par vent, qui ont été collectés de Décembre 1999 à Décembre 2000. En comparant les données météorologiques et les données de l'accumulation de la neige, nous avons établi que la plupart des précipitations au centre de l'Antarctique était formée pendant des conditions de "ciel clair" (diamond dust). Par ailleurs les données de sondages aérologiques suggèrent que la température moyenne de condensation à Vostok est très proche de la température de la limite haute de la couche d'inversion. Pour un site donné, on a mis en évidence que la variabilité temporelle de la composition isotopique de la neige est influencé par le relief de la surface de neige (microrelief, "méso-dunes" et probablement méga-dunes). Cet effet diminue considérablement le rapport "signal sur bruit". On met en évidence que les variations saisonnières des teneurs isotopiques des précipitations suivent bien celles de la température locale. L'excès en deutérium varie au cours de l'année et il est, quant à lui, influencé par la source d'évaporation. Pour les 40 dernières années et pour l'échelle de temps multi-décennale, une relation linéaire significative est observée entre la moyenne de la composition isotopique de la neige et la température de l'air près de la surface. Mais pour des échelles de temps plus courtes (10 ans) les variations isotopiques apparaissent dominées par les conditions de la source. Pendant les deux derniers siècles, l'accumulation et la composition isotopique de la neige présentent des fluctuations avec une période de 50 ans environ, qui ont des similarités avec l'indice de l'Oscillation Décennale du Pacifique. Si cela est vérifié, cette relation impliquerait une téléconnection climatique de l'Antarctique centrale avec Pacifique tropical. Par ailleurs, l'excès en deutérium suggère un changement brusque vers 1963 dans la circulation atmosphérique qui alimente la région de Vostok en l'humidité.

L'atmosphère des hautes latitudes sud demeure très peu documentée à l'heure actuelle. Or, la compréhension des systèmes biogéochimiques à ces latitudes est nécessaire tant pour l'interprétation des archives glaciaires que pour évaluer l'importance de ces régions dans le climat actuel et futur, en liaison notamment avec les émissions océaniques de diméthylsulfure (DMS). Des études de distributions en taille de l'aérosol sur la base côtière antarctique de Dumont d'Urville ont permis de mieux cerner les problèmes de fractionnement de l'aérosol marin par rapport à la composition de l'eau de mer. En hiver, un déficit de sulfate par rapport au sodium ou au chlore est observé. La chronologie et l'ampleur de ce phénomène ont pu être précisés, dans le but de quantifier le nss-sulfate hivernal. Un second phénomène de fractionnement du sel de mer, lié à la remobilisation de HCI et traduit par un déficit de chlore par rapport au sodium, se révèle peu marqué sur le site suite à la combinaison de plusieurs paramètres, parmi lesquels l'influence locale d'émissions ornithogéniques alcalines. Depuis fin 1998, des mesures de DMS et de DMSO complètent les prélèvements d'aérosols continus depuis 1991. Après trois années de mesures, il apparaît une très grande variabilité interannuelle des teneurs en DMS liée à une variabilité importante de la source océanique locale en été. Si les niveaux de DMS paraissent fortement dépendants de cette source à une échelle locale, les espèces oxydées (MSA, sulfate) sont plutôt représentatives en été des émissions océaniques à une échelle régionale. En hiver, leurs teneurs sont très faibles (15 à 20 ng/m3 de nss-sulfate, hors perturbation volcanique, dont une moitié serait d'origine biogénique). Enfin, des études de photochimie en été ont confirmé l'importance d'un puits hétérogène du DMSO qui semble avoir une dépendance photochimique et s'inscrit comme étant une voie de production privilégiée de MSA particulaire.

Le mercure est un métal toxique émis dans l'environnement par des sources naturelles et anthropiques. Le cycle naturel de ce polluant, a fortiori la réponse des écosystèmes à ces perturbations par l'Homme, est encore mal quantifié. L'air des névés polaires constitue une archive de la composition passée de l'atmosphère qui permet d'améliorer la connaissance du cycle du mercure. En juillet 2009, à NEEM, Groenland du Nord, nous avons mesuré les concentrations en mercure élémentaire gazeux dans l'air du névé entre 30 et 76 mètres de profondeur, et, grâce à des travaux de modélisation des processus de diffusion, reconstruit les concentrations atmosphériques en Hg° depuis les années 1950, à une échelle spatiale représentative de l'Arctique. Le signal reconstruit présente des concentrations maximales dans les années 1970, synchrones avec le maximum d'émissions anthropiques supposé, puis une décroissance très peu marquée, suggérant une persistance de cette contamination dans l'Arctique. Le plateau antarctique revêt également un intérêt paléoenvironnemental dans la perspective de caractériser l'impact des perturbations anthropiques aux hautes latitudes sud. Néanmoins, la réactivité du mercure dans cette région est méconnue, et sa caractérisation est nécessaire. Nous avons mené une campagne de mesure en janvier 2009 sur le site de Dôme C, Antarctique de l'Est et observé une réactivité marquée du mercure dans la basse troposphère et dans l'air interstitiel du manteau neigeux, mettant en évidence l'existence de processus d'oxydation puissants dans ces deux matrices. De fortes concentrations en mercure ont aussi été mesurées dans la neige de surface.