Sous l'effet de charges en surface, la surface terrestre s'affaisse pour retrouver un état d'équilibre isostatique. Inversement, quand ces charges viennent à diparaître, la surface remonte pour retrouver sa configuration initiale. Dans le cas d'une calotte glaciaire, la taille mais aussi la vitesse à laquelle peut évoluer la masse de glace est à l'origine d'une réponse hautement spécifique et particulièrement de comportements transitoires associés. Après avoir introduit la notion d'isostasie, nous décrivons quelles peuvent être les conséquences à la fois de la taille et de la vitesse d'évolution de la charge sur la réponse terrestre. Cela permet de clarifier la notion de "contexte glaciaire" en matière de charge à la surface de la Terre. Ensuite sont présentés les nombreux impacts que cette réponse isostatique peut avoir en retour sur la dynamique glaciaire, ainsi que les façons d'en rendre compte. Cette première partie aura montré tout l'intérêt de developper un modèle de Terre élaboré dont la physique à la base s'avère plus rigoureuse que dans les paramétrisations usuelles utilisées jusqu'à maintenant dans les modèles de glace. Les caractéristiques majeures ainsi que les équations à la base de ce modèle sont alors présentées en détail et révèlent le degré de sophistication de cette nouvelle approche. La richesse de données issues de mesures sur les zones de rebond postglaciaire comme la Fennoscandie sont un excellent moyen de valider pareil modèle. En comparant les résultats du rebond simulé avec ces données il est possible d'estimer le réalisme du modèle. Bien que l'ensemble de ces données ne soit pas simultanément reproduit de manière parfaite, le modèle montre cependant un comportement suffisamment réaliste pour pouvoir par la suite être couplé à un modèle de glace. Les résultats de cette incorporation sont alors comparés à ceux obtenus avec les paramétrisations usuelles au cours d'une simulation de la calotte Antarctique au cours du dernier cycle glaciaire. Ces résultats confirment le relativement bon comportement de notre modèle complet de Terre par rapport aux autres approches, ce qui tendrait à discréditer l'emploi de certaines de ces paramétrisations de la réponse terrestre dans les modèles évolutifs de calottes glaciaires.

Le bilan de masse en surface (noté BMS ; l'accumulation de neige diminuée de l'ablation) de la calotte glaciaire antarctique est sensible aux paramètres climatiques et contribue directement aux variations du niveau moyen des mers. Il est donc important, dans le cadre de la prévision du changement climatique, de développer des outils capables de simuler les processus physiques régissant le bilan de masse en surface antarctique. L'approche développée dans cette thèse consiste à utiliser une cascade de modèles atmosphériques allant de la grande échelle vers l'échelle locale. Ainsi, un modèle climatique régional (Modèle atmosphérique régional, MAR), forcé par des réanalyses du Centre européen pour les prévisions météorologiques à moyen terme (CEPMMT), fournit à un modèle diagnostique de désagrégation physique des précipitations les champs météorologiques nécessaires calculés à l'échelle régionale (typiquement, de résolution 40 km). Dans un premier temps, il est montré que le bilan de masse en surface généré par le MAR est conforme aux observations dans la plupart des régions. Toutefois, le ruissellement est surestimé ; ce problème disparaît en introduisant une dépendance de l'albédo avec la distance zénithale de l'astre solaire. Dans un second temps, il est montré que malgré la relative simplicité des paramétrisations physiques du désagrégateur, la connaissance du relief de fine échelle (de résolution 5 km) permet d'améliorer la variabilité spatiale de la précipitation, et, par conséquent, du BMS, sur les régions côtières de l'Antarctique. La validation est menée à l'aide, notamment, de mesures de hauteurs de neige délivrées par des stations météorologiques automatiques. Sur le site côtier de Law Dome, le gradient d'accumulation nette est davantage dû au forage orographique subi par la précipitation qu'au processus de chasse-neige. Le modèle de désagrégation sous-estime fortement la précipitation sur le plateau Antarctique, où les nuages stratosphériques polaires associés au refroidissement radiatif pourraient jouer un rôle dans la génération de la précipitation pendant la nuit polaire.