Academic literature on the topic 'Phénomènes de la calotte polaire'

Lorsque le champ magnétique interplanétaire est dirigé vers le Nord, les régions polaires de la magnétosphère terrestre, situées à haute altitude, sont peuplées d'électrons accélérés. L'ionosphère polaire voit quant à elle l'apparition d'aurores. Grâce aux observations magnétosphériques de la mission multi-satellites Cluster et à notre modélisation du transport de particules le long des lignes de force du champ magnétique terrestre B, nous avons montré, au cours de cette thèse, que les électrons à l'origine de ces aurores proviennent de la couche frontière de la magnétosphère. Localisée à très haute altitude, cette couche sépare les milieux terrestre et interplanétaire. Accélérés par des champs électriques parallèles à B, situés à des altitudes aussi bien inférieures que supérieures à celles des satellites, ces électrons précipitent dans l'ionosphère avec des énergies de l'ordre du keV. Ils sont accompagnés par des échappements d'ions ionosphériques et l'ensemble des particules transporte des courants alignés le long de B et dirigés vers la magnétosphère, équivalents à quelques µA/m² aux altitudes ionosphériques. Le système de courant se referme par des courants alignés dirigés vers la Terre, d'intensité comparable, portés par des faisceaux d'électrons étroits et variables. Très alignés le long de B, ces électrons s'échappent de l'ionosphère avec des énergies de ~ 50 eV et déclenchent des instabilités faisceau-plasma à l'origine d'ondes électrostatiques large bande observées par Cluster et reproduites par les modèles.

Lorsque le champ magnétique interplanétaire est dirigé vers le Nord, les régions polaires de la magnétosphère terrestre, situées à haute altitude, sont peuplées d'électrons accélérés. L'ionosphère polaire voit quant à elle l'apparition d'aurores. Grâce aux observations magnétosphériques de la mission multi-satellites Cluster et à notre modélisation du transport de particules le long des lignes de force du champ magnétique terrestre B, nous avons montré, au cours de cette thèse, que les électrons à l'origine de ces aurores proviennent de la couche frontière de la magnétosphère. Localisée à très haute altitude, cette couche sépare les milieux terrestre et interplanétaire. Accélérés par des champs électriques parallèles à B, situés à des altitudes aussi bien inférieures que supérieures à celles des satellites, ces électrons précipitent dans l'ionosphère avec des énergies de l'ordre du keV. Ils sont accompagnés par des échappements d'ions ionosphériques et l'ensemble des particules transporte des courants alignés le long de B et dirigés vers la magnétosphère, équivalents à quelques µA/m² aux altitudes ionosphériques. Le système de courant se referme par des courants alignés dirigés vers la Terre, d'intensité comparable, portés par des faisceaux d'électrons étroits et variables. Très alignés le long de B, ces électrons s'échappent de l'ionosphère avec des énergies de ~ 50 eV et déclenchent des instabilités faisceau-plasma à l'origine d'ondes électrostatiques large bande observées par Cluster et reproduites par les modèles.

L'objectif de cette thèse est de mieux comprendre les variations temporelles des signaux altimétriques sur la calotte polaire Antarctique. Nous exploitons les observations effectuées par l'altimètre à bord d'ENVISAT entre janvier 2003 et décembre 2007. Ces observations s'étendent jusqu'à 82°Sud, ce qui permet de couvrir environ 80% du continent Antarctique. Pendant la période d'étude deux fréquences sont exploitables (Bande S, soit 3.2GHz, et Bande Ku soit 13.6GHz), ce qui permet de mieux cerner la sensibilité de la mesure aux variations d'état du manteau neigeux couvrant la calotte polaire Antarctique. L'état de la surface neigeuse en Antarctique varie dans l'espace comme dans le temps, en réponse aux événements météorologiques de proche surface, et tout particulièrement les vents. Selon l'état de la surface, l'onde électromagnétique émise par l'altimètre va pénétrer plus ou moins profondément dans le manteau neigeux, modifiant la forme de l'écho enregistré lorsque cette onde retourne au satellite. La précision des mesures de topographie de la calotte polaire dépend alors de la connaissance des propriétés de surface, ainsi que de la compréhension des mécanismes contrôlant l'interaction des ondes électromagnétiques avec la surface neigeuse. Nous étudions dans cette thèse comment la forme de l'écho altimétrique évolue lorsque les conditions de vent à la surface de la calotte polaire changent, et affectent l'état de la surface. Nous commençons par étudier la variabilité spatiale et temporelle des observations altimétriques d'une part, et des intensités de vent d'autre part (réanalyses des modèles atmosphériques globaux NCEP/NCAR, NCEP/DOE et ERA-Interim). Nous en identifions les échelles de variations spatiales et temporelles caractéristiques. Cela nous permet de concevoir les traitements à appliquer afin d'homogénéiser les différentes sources d'information. Nous constatons que les variations de nos séries temporelles décrivent essentiellement un signal annuel, toujours maximal en été austral pour les variations de vents. En revanche les variations annuelles de la forme d'onde ne sont pas homogènes à l'échelle de la calotte polaire, mais mettent en évidence deux régions principales. Nous utilisons ensuite les propriétés de la forme d'onde altimétrique, et leurs variations temporelles pour reconstruire les variations d'intensité de vents. Nous obtenons de très bonnes cohérences, ce qui atteste la sensibilité de la mesure altimétrique aux changements d'état de la surface induits par les variations de vent. Nous utilisons en dernier lieu un modèle régional (MAR/LGGE), qui prend en compte les interactions des vents avec l'état de la surface avec une résolution spatiale plus fine. Cela nous permet de confirmer les résultats déjà obtenus, et de mettre en valeur les spécificités du modèle MAR. Les méthodes mises en jeu sont empiriques et les résutats obtenus qualitatifs. Nous avons toutefois identifié deux régions majeures, réagissant différemment aux variations de vents, que nous expliquons par l'existence de structures de surface particulières, et un impact différent des vents sur l'état de la surface.

L'objectif de cette thèse est de mieux comprendre les variations temporelles des signaux altimétriques sur la calotte polaire Antarctique. Nous exploitons les observations effectuées par l'altimètre à bord d'ENVISAT entre janvier 2003 et décembre 2007. Ces observations s'étendent jusqu'à 82°Sud, ce qui permet de couvrir environ 80% du continent Antarctique. Pendant la période d'étude deux fréquences sont exploitables (Bande S, soit 3. 2GHz, et Bande Ku soit 13. 6GHz), ce qui permet de mieux cerner la sensibilité de la mesure aux variations d'état du manteau neigeux couvrant la calotte polaire Antarctique. L'état de la surface neigeuse en Antarctique varie dans l'espace comme dans le temps, en réponse aux événements météorologiques de proche surface, et tout particulièrement les vents. Selon l'état de la surface, l'onde électromagnétique émise par l'altimètre va pénétrer plus ou moins profondément dans le manteau neigeux, modifiant la forme de l'écho enregistré lorsque cette onde retourne au satellite. La précision des mesures de topographie de la calotte polaire dépend alors de la connaissance des propriétés de surface, ainsi que de la compréhension des mécanismes contrôlant l'interaction des ondes électromagnétiques avec la surface neigeuse. Nous étudions dans cette thèse comment la forme de l'écho altimétrique évolue lorsque les conditions de vent à la surface de la calotte polaire changent, et affectent l'état de la surface. Nous commençons par étudier la variabilité spatiale et temporelle des observations altimétriques d'une part, et des intensités de vent d'autre part (réanalyses des modèles atmosphériques globaux NCEP/NCAR, NCEP/DOE et ERA-Interim). Nous en identifions les échelles de variations spatiales et temporelles caractéristiques. Cela nous permet de concevoir les traitements à appliquer afin d'homogénéiser les différentes sources d'information. Nous constatons que les variations de nos séries temporelles décrivent essentiellement un signal annuel, toujours maximal en été austral pour les variations de vents. En revanche les variations annuelles de la forme d'onde ne sont pas homogènes à l'échelle de la calotte polaire, mais mettent en évidence deux régions principales. Nous utilisons ensuite les propriétés de la forme d'onde altimétrique, et leurs variations temporelles pour reconstruire les variations d'intensité de vents
This work aims at improving our understanding of the altimetric time series acquired over the Antarctic Ice Sheet. Dual frequency data (S Band - 3. 2GHz and Ku Band - 13. 6GHz) from thealtimeter onboard the ENVISAT satellite are used, during a five year time period from january2003 until december 2007. These data cover around 80% of the surface of the Antarctic continent,up to 82°S. Having data in two different frequencies is valuable when it comes to better estimatethe altimeter sensitivity regarding snow surface property changes. Over the Antarctic ice sheet, snow surface changes with respect to space and time, beingaffected by meteorological conditions close to the surface, and especially winds. The altimetricwave penetrates more or less deeply beneath the surface, depending on snow surface and subsurfaceproperties. As a result, when the wave comes back to the satellite, the recorded signal, namedwaveform, is more or less distorted. The accuracy of the ice sheet topographic changes computedthanks to satellite altimetric techniques depends on our knowledge of the processes inducing thisdistortion. The purpose of the present work is to better understand the effect of changing windconditions on altimetric data. Winds in Antarctica are indeed famous for their strength and theirimpact on the snow surface state. First, spatial and temporal variability of the altimetric data on the one hand, and of wind speedreanalysis fields (from ERA-Interim, NCEP/NCAR and NCEP/DOE projects) on the other handare studied. We estimate spatial and temporal typical length scales for all datasets. As a result, weare able to smooth the data, so that all datasets have the same spatial and temporal caractericticlength scales. Furthermore, we note that our time series are well described by an annual signal. This annual cycle shows that whereas wind speed would always be maximum in austral winter,altimetric seasonal cycles have very different behaviors depending on the location. .

Les prochaines modèles climatiques comprendront un modèle de calotte polaire afin de tenir en compte la dynamique de la glace et les interactions glace-océans dans ses projections. D'une part, l'océan Austral (SO) pilote l'accélération des glaciers de l'Antarctique via une augmentation de la fonte basale des ice shelves. D'autre part, l'accélération de la décharge de glace de l'Antarctic Ice Sheet (AIS) contribue à la montée du niveau de la mer et est susceptible de devenir le plus grand contributeur de la cryosphère d'ici la fin du siècle. En outre, l'adoucissement relié, peut avoir des répercussions importantes sur la glace de mer et sur la formation des eaux profondes. Cependant, on ne sait pas encore comment les modèles d'océan et de calotte polaire des futurs systèmes couplés vont représenter les interactions glace-océan, causes et conséquences du déséquilibre de masse de AIS. Ici, dans ce travail, les différents aspects des modèles de océan et calotte polaire ont été étudiés. Une première étape de cette thèse a été concentrée à la représentation des flux d'eau douce glaciaires dans les modèles océaniques actuels. Basé dans estimations glaciologiques, la fonte basal des ice shelves a été répartie dans une grille de ORCA025, et les taux de production d'icebergs ont été appliqués dans une version améliorée du modèle d'iceberg NEMO-ICB. Cette étude préliminaire a été utilisé pour produire une climatologie d'eau de fonte provenant des icebergs, valable pour forcer les modèles de océan actuels. Ce travail montre l'importance de représenter les flux d'eau de fonte des icebergs lors de la modélisation de la glace de mer, qui peut être obtenu en utilisant notre climatologie. Ces améliorations ont été pris en compte dans l'étude de la réponse du modèle de océan a la perte de masse de AIS. Cette étude considère une perturbation réaliste de l'eau douce glaciaire aussi près que possible de sa représentation dans les futurs modèles couplés ice-sheet/océan. Selon nos résultats, jusqu'à 50% des changements récents de volume de la glace de mer pourrait être causée par le bilan masse de l'AIS. Le forçage en eau douce glaciaire semble être cruciale pour représenter correctement les interactions glace-océan et projeter la glace de mer dans les futurs systèmes couplés. Cependant, l'estimation de l'apport d'eau douce glaciaire dans les modèles climatiques futurs sera fortement affecté par la capacité des modèles de calotte polaire de reproduire les migrations de grounding line des glacières de "marine ice sheets". Les modèles de calotte polaire actuels présentent grandes incertitudes liées aux différents réglages. Dans le contexte des futurs modèles climatiques, nous avons étudié la sensibilité des retraites de la grounding line produites par l'océan à l'application de deux lois de frottement différentes et deux différentes approximations du stress glacier. Les modèle réagit de façon presque similaire aux approximations SSA ou SSA *. Par contre, les différences dans la contribution du glacier à l'élévation du niveau de la mer peuvent être jusqu'à 50% en fonction de la loi de frottement considéré. La loi de friction Schoof, la plus physique, est nettement plus réactif aux perturbations océaniques que la loi Weertman, et devrait être pris en compte dans les systèmes couplés futurs. Ce travail souligne que les incertitudes liées aux estimations des modèles de la calotte glaciaire de migrations de grounding line peuvent contribuer non seulement à des incertitudes du futur niveau de la mer, mais aussi de la glace de mer à travers des interactions glace-océan dans les futures models climatiques. Tel conclusion suggère la nécessité d'améliorer la représentation de la fonte basal des ice shelves et le frottement du glacier, afin d'améliorer les projections climatiques des modèles climatiques, dans lequel la distribution spatiale et saisonnière des eau douce glaciaires peut jouer un rôle important en établir la glace de mer
The next generation of climate models will include an ice-sheet model in order to improve the ice sheet mass balance projections by accounting for the ice dynamics and ice-oceans interactions. On the one hand, the Southern Ocean (SO) is indeed driving the acceleration of the Antarctic outlet glaciers via an increase in the basal melting of the ice shelves. On the other hand, the increasing ice discharge from Antarctic Ice Sheet (AIS) contributes to the current sea level rise and is likely to become the largest cryospheric contributor to sea level rise by the end of the current century. In addition, the related freshening may have significant implications on future sea-ice cover and on bottom water formation. However, it is not clear yet how the ocean and ice-sheet components of future coupled systems will account for the ice-ocean interactions, which are both causes and consequences of the AIS mass imbalance. Here in this work, different aspects of the standalone ocean and ice-sheet components have been investigated. A first step of this thesis has been focused in the representation of the glacial freshwater fluxes in current ocean models. Based on recent glaciological estimates, the ice shelf basal melting fluxes have been spatially distributed in an ORCA025 grid, and the calving rates have been applied into an improved version of the NEMO-ICB iceberg model. This preliminary study has been used to produce a monthly iceberg meltwater climatology, to be used to force current ocean models. This work shows the importance of representing the iceberg meltwater fluxes when modeling sea ice, which can be inexpensively achieve by using our climatology. The improvements in the representation of the glacial freshwater fluxes have been considered in the study of the ocean model response to the Antarctic mass imbalance. This study considers a realistic perturbation in the glacial freshwater forcing as close as possible as it will be represented in future ice-sheet/ocean models. According to our results, up to 50% of the recent Antarctic sea ice volume changes might be caused by the observed decadal AIS mass imbalance rate. Glacial freshwater forcing appears to be crucial to correctly represent the ice-ocean interactions and projecting sea ice cover of future coupled systems. However, the estimation of the glacial freshwater input in future climate models will be strongly dependent upon the capacity of ice-sheet models to reproduce the grounding line migrations of marine ice sheet glaciers. Current ice-sheet models present large uncertainties related to different parametrizations. In the context of the future climate models, we have studied the sensitivity of ocean-driven grounding line retreats to the application of two different friction laws and two different englacial stress approximations. The model responses almost indistinctively to either the SSA or the SSA* englacial stress approximations. However, differences in the contribution of the glacier to the sea level rise can be up to 50% depending on the friction law considered. The more physically constrained Schoof friction law is significantly more reactive to the ocean perturbations than Weertman law and should be considered in future coupled systems. This work underlines that uncertainties related to the ice sheet model estimates of grounding line migrations may not only contribute to uncertainties in sea level projections, but also the sea ice cover through the ice-ocean interaction in future ocean models.This conclusion suggests the need for improving the representation of both the ice shelf basal melting and the glacier interaction with the bedrock, in order to improve the climate projections of future climate models, in which the spatial and seasonal distribution of the glacial freshwater fluxes may play an important role in setting the sea ice cover

La calotte polaire antarctique est constituée d'une épaisse couche de glace surmontée d'un névé d'une centaine de mètres d'épaisseur. Si les observations satellites, avec une couverture quasi-globale du continent, ont permis depuis bientôt 40 ans de fournir de nombreuses informations sur sa topographie et quelques propriétés de sa surface, la connaissance de son manteau neigeux est encore très limitée. Les précipitations qui se déposent sont balayées par les vents, s'accumulent pour former des couches. En s'enfouissant, les grains de neige qui les constituent subissent un métamorphisme qui modifie leur taille et leurs caractéristiques. Pourtant, Les taux d'accumulations ou les tailles des grains de neige, dont la signature est contenue dans les observations de la proche surface, sont des indicateurs climatiques clés, jusqu'alors mal connus.

Depuis 2002 et le lancement de ENVISAT, on dispose d'un altimètre radar qui couvre 80 \% de la calotte polaire Antarctique, dont la particularité est d'acquérir des signaux à deux fréquences différentes (bande S à 3.2 GHz et bande Ku à 13.6 GHz). Ces deux ondes pénètrent dans le manteau neigeux sur plusieurs mètres et ont des sensibilités aux propriétés de la neige différentes. Ainsi, l'idée de cette thèse est d'utiliser cette double information pour retrouver les propriétés du manteau neigeux.

On se propose de résoudre cette problématique par une analyse et une modélisation des signaux altimétriques bi-fréquences sur la calotte polaire, puis par leur inversion. On se penche tout d'abord sur quelques études de cas pour estimer la sensibilité des signaux aux différentes propriétés de la neige: i/ On montre tout d'abord que le signal altimétrique est sensible à la rugosité de la surface à différentes échelles, puis ii/ que le signal altimétrique est sujet à des variations saisonnières causées par la densification de la neige en surface, et enfin iii/ que les ondes radars sont réfléchies par des strates en profondeur.

Un modèle de l'interaction de l'onde avec le manteau neigeux est réalisé simultanément aux deux fréquences, afin de permettre une comparaison de ces signaux entre eux. Les résultats du modèle sont utilisés pour expliquer les variations saisonnières précédemment observées. Finalement, les paramètres du manteau neigeux sont estimés à l'échelle de la calotte polaire antarctique. Les tailles de grains retrouvées présentent un grossissement vers l'intérieur du continent. La densité montre des variations saisonnières de plusieurs g.cm3 notamment sur les côtes antarctiques. Certaines régions présentent un état de surface de la neige particulièrement lisse (Dronning Maud Land, par exemple).

La donnée in situ de l'état de surface de la neige étant quasi inexistante sur les calottes polaires, on développe finalement un protocole de mesure de la rugosité de la neige, qui est testé sur un glacier du Spitzberg.

L'évolution des calottes polaires est régie à la fois par une dynamique d'écoulement complexe et par des mécanismes tel le glissement à la base, la température de la glace ou le bilan de masse en surface. De plus, de nombreuses boucles de rétroactions sont constatées entre les différents phénomènes impliquées. Tout ceci rend la modélisation de cette évolution complexe. Malgré tout, un certain nombre de modèles ont été développés dans cette optique. Ceux-ci font tous intervenir des paramètres influents qui dans certains cas sont peu ou pas connus. Ils nécessitent donc d'être correctement spécifiés. L'assimilation de données peut permettre une meilleure estimation de ces paramètres grâce à l'utilisation d'observations qui sont peu nombreuses en glaciologie. Dans cette thèse, nous nous intéressons à la mise en place de systèmes d'assimilation performants pour deux problèmes inverses concernant l'évolution des calottes polaires. Pour mieux nous concentrer sur ce point, nous avons travaillé avec un modèle d'évolution de calotte simplifié (appelé Winnie) qui, cependant, représente bien la plupart des processus complexes de la dynamique de la glace, et permet de travailler à différentes échelles de temps. Dans un premier temps, nous mettons en place une approche 4D-Var pour la reconstruction de l'évolution d'un paramètre climatique influant sur l'évolution d'une calotte sur une échelle de temps typique de 20 000 ans. Elle nécessite notamment l'écriture du code adjoint du modèle. Dans un second temps, nous nous intéressons au problème du spin-up. Ce problème de calibration du modèle pour des simulations à échelle de temps courtes (pas plus de 100 ans) consiste plus particulièrement en la reconstruction conjointe de l'état initial, de la topographie du socle rocheux et des paramètres de glissement basal. Nous développons ici une approche filtre de Kalman d'ensemble pour résoudre ce problème.

Le cristal de glace est l'un des matériaux naturels les plus anisotropes. L'analyse des glaces polaires indique que les cristaux s'orientent selon des directions privilégiées et que la fabrique (distribution des orientations cristallines des grains) du polycristal de glace est le résultat de l'histoire des déformations que ce polycristal a subi. Le comportement macroscopique du polycristal dépend de la fabrique et peut être lui aussi fortement anisotrope. Il a été montré que cette anisotropie est propice à créer des perturbations de la stratigraphie et qu'elle influence fortement l'écoulement de la calotte. L'objectif de cette étude est de construire un modèle d'écoulement pour la glace polaire anisotrope et l'évolution de sa fabrique. Dans ce but, nous modélisons le monocristal de glace comme un matériau continu orthotrope de révolution autour de son axe c et nous le comparons à un modèle de grain à plans de glissement. La fabrique est décrite de manière continue par le tenseur d'orientation d'ordre 2 et une fonction de fermeture pour le tenseur d'orientation d'ordre 4. Par homogénéisation, en supposant des contraintes ou des vitesses de déformation uniformes, nous obtenons des solutions analytiques pour le comportement du polycristal de glace et l'évolution de sa fabrique. A partir de ces solutions, nous adoptons une loi de comportement orthotrope linéaire pour le polycristal et une équation d'évolution pour le tenseur d'orientation du second ordre. Ces équations permettent de bien reproduire les résultats du modèle auto-cohérent utilisant une description discrète de la fabrique, avec un gain important au niveau du temps de calcul et du nombre de variables nécessaires pour décrire la fabrique. Ces équations sont implantées ensuite dans un code aux éléments finis utilisé pour simuler l'écoulement d'une calotte polaire présentant une anisotropie induite évolutive. Par des tests synthétiques, nous montrons l'influence de l'anisotropie sur l'écoulement de la calotte. Enfin, nous présentons une analyse des données de fabrique du forage de Dôme C à la lumière de notre modèle
The ice crystal is one of the most anisotropic natural materials. Observations of polar ice show that the crystals are orientated along very marked privileged directions and that the fabric of the ice polycristal is the result of the strain history that the polycristal has undergone. The macroscopic polycristal behaviour depends on the fabric and can be strongly anisotropic. It has been shown that anisotropy can introduce the occurrence of stratigraphy disturbances and influences strongly the flow of the ice sheet. The aim of this work is to build an ice flow model for anisotropic ice and the evolution of its fabric. To this end, we model the ice monocrystal as a transversely isotropic continuous medium around its c-axis and we compare it with a slip-planes grain model. The ice fabric is described in a continuous way by the second order orientation tensor and a closure approximation for the fourth ordre orientation tensor. By homogenization, assuming uniform stresses or strain-rates, we obtain analytical solutions for the polycrystal behaviour and the evolution of its fabric. From these solutions, we derive a linear orthotropic flow law for the ice polycrystal and an equation for the evolution of the second order orientation tensor. These equations are able to reproduce the results of the self-consistent model using a discrete fabric description, with an important improvement as regards the computing time and the number of variables necessary to describe the fabric. These equations are then implemented in a finite element code, used for simulating the flow of an ice sheet with strain-induced evolving anisotropy. By performing synthetic tests, we show the influence of anisotropy on the ice flow. Finally, through our model results, we propose an interpretation of the Dome-C ice-core fabric data

D'importants dépôts de sulfates ont été découverts dans certaines régions équatoriales de Mars : les canyons de Valles Marineris, les terrains chaotiques, et les plaines de Meridiani. Afin d'expliquer la chimie de ces dépôts, ainsi que la formation d'un tel volume de sédiments de même composition, sur des surfaces aussi larges et variées, l'une des hypothèses les plus récentes suggère une origine glaciaire. Le but de cette thèse est d'évaluer la plausibilité de ce mode de formation, par une étude morphologique et minéralogique comparative des sulfates trouvés sur la calotte polaire Nord et sur les terrains chaotiques (Aram Chaos). L'analyse des sédiments superficiels de la calotte polaire montre que ceux-ci sont riches en gypse. Ces particules riches en sulfates proviennent directement de la calotte polaire, s'accumulent en surface par ablation de la glace, et sont transportées par le vent dans les dunes circum-polaires. L'analyse des affleurements observés à l'équateur sur Aram Chaos montre des dépôts stratifiés clairs, semblables à ceux de Valles Marineris et de Meridiani Planum ; ils contiennent un mélange de sulfates monohydratés, d'oxydes de fer et de minéraux hydratés. Ces dépôts pourraient également avoir une origine glaciaire car, durant les périodes de haute obliquité, d'anciennes accumulations de glace ont pu se former aux basses latitudes. De plus, l'observation d'analogues terrestres en Antarctique révèle que des sulfates se forment effectivement dans la glace par des réactions post-dépôt, et peuvent être compatibles avec la chimie des minéraux observés à l'équateur martien, si l'accumulation de la glace coïncide avec le dépôt de matériaux volcaniques.

L'observation de l'orientation des grains des carottes forées en Antarctique et au Groenland montre que la glace polaire acquiert une fabrique représentative de l'histoire des déformations qu'elle a subies. Le comportement viscoplastique du cristal de glace étant fortement anisotrope, il en résulte une anisotropie macroscopique du polycristal qui dépend des orientations préférentielles des grains qui le composent. A l'échelle du polycristal, le comportement viscoplastique est obtenu par homogénéisation sur tous les cristaux, en supposant leur comportement linéaire orthotrope de révolution et une répartition homogène des contraintes dans tout le polycristal (modéle statique). La fabrique du polycristal est décrite par une fonction continue de distribution des orientations (ODF). Avec ces hypothèses, nous avons obtenu des résultats analytiques pour le comportement orthotrope d'un polycristal et l'évolution de sa fabrique, qui permettent de représenter la fabrique par une ODF ne dépendant que de trois paramètres. A l'échelle de la calotte polaire, les champs de vitesses et fabriques de l'écoulement stationnaire sont calculés par une méthode couplée, qui alterne : - le calcul des vitesses, pour un champ de fabriques donné, par la méthode des éléments finis, - le calcul de la fabrique, pour un champ de vitesses donné, par une méthode qui utilise les lignes de courant de l'écoulement, jusqu'à ce que la convergence soit atteinte, à la fois sur les vitesses et la fabrique. Nous étudions l'influence du champ de températures et du comportement du grain sur la formation des fabriques au sein des calottes polaires. Le modèle est appliqué à la ligne d'écoulement GRIP-GISP2 (Groenland Central) et les résultats sont comparès aux mesures de terrain.