In 1974, newly available satellite observations unveiled the presence of a giant ice-free area within the Antarctic ice pack, which persisted throughout the winter, and formed again in the next two winters. Subsequent research showed that deep convective overturning kept the waters ice-free, through the massive release of heat rising from the deep sea. While the polynya has aroused continued interest among climate scientists, it has not reappeared since 1976. Here we use model experiments to show that deep convection in the Southern Ocean, common in current generation climate models, is highly sensitive to anthropogenic forcing, and ceases in many models when forced by a high emissions climate change scenario. The slowdown in deep ventilation follows from the gradual freshening of polar surface waters, a trend which is borne out by observations over recent decades. Our results suggest that deep convection in the Southern Ocean will be less common in future, and may have already been significantly reduced compared to the pre-industrial period, with important consequences for ocean circulation and climate. / En 1974, des observations satellite nouvellement disponibles révélèrent la présence d'une géante surface d'eau libre au sein de la glace de mer entourant l'Antarctique, qui persista tout au long de l'hiver et réapparut les deux hivers suivants. Les recherches qui suivirent montrèrent que les eaux étaient maintenues libres de glace par la convection profonde, permettant à une grande quantité de chaleur de remonter des profondeurs pour être ensuite libérée dans l'atmosphère. Si la polynya continue de susciter l'intérêt des climatologues, elle n'est cependant pas réapparue depuis 1976. Nous utilisons ici des expériences de modélisation pour montrer que la convection profonde dans l'Océan Austral, commune dans les modèles de climat actuels, est fortement sensible au forçage anthropique, et cesse dans beaucoup de modèles quand ceux-ci sont forcés par un scénario de fortes émissions. Le ralentissement de la ventilation profonde résulte de la baisse progressive de la salinité des eaux de surface, une tendance corroborée par les observations des dernières décennies. Nos résultats suggèrent que la convection profonde dans l'Océan Austral sera moins fréquente dans le futur, et a peut-être déjà été significativement affaiblie relativement à la période préindustrielle, avec d'importantes conséquences pour la circulation océanique et le climat.

Current sea-ice models evolve towards higher spatial resolutions and levels of convergence of their numerical solution, which refines their simulated ice deformation fields, while still employing the values of the dynamic parameters used in earlier models. In this work, we investigate the effect of spatial resolution and convergence on the choice of optimal model parameters by comparing drag coefficients and ice compressive and shear strength sensitivity studies conducted with a viscous-plastic sea-ice model of the Arctic. We observe that convergence impacts the optimal parameters greatly and that a low converged solution induces ice drift biases that cannot be rectified by changing the parameter values. The performed experiments do not allow drawing conclusions on the effect of spatial resolution but indicate that varying the resolution of the model wind forcing alters the optimal parameters if this forcing is derived using a non-linear spatial interpolation scheme. / Les modèles de glace de mer évoluent aujourd'hui vers une résolution spatiale et un niveau de convergence numérique accrus, ce qui rafine les champs de déformation simulés, mais emploient toujours d'anciennes valeurs de paramètres dynamiques. Dans cette étude, nous examinons l'effet de la résolution et de la convergence sur les paramètres dynamiques optimaux en étudiant la sensibilité d'un modèle visco-plastique de glace de mer à la valeur des coefficients de friction et de résistance en compression et tension de la glace. Nous observons que le niveau de convergence affecte largement les paramètres optimaux et qu'une convergence insuffisante induit des biais sur la dérive des glaces qui ne peuvent être rectifiés en changeant leur valeur. Les expériences effectuées ne permettent pas de tirer de conclusions quant à l'effet de la résolution du modèle mais indiquent que la résolution des vents de forçage affecte les paramètres optimaux s'ils sont dérivés d'une méthode d'interpolation spatiale non-linéaire.

Ice algae are an important component of the carbon cycle in the Arctic, and can therefore have an impact on climate. I investigate the dynamics of an ice algae bloom by coupling an algae-nutrient model [Lavoie et al., 2005] with a thermodynamic sea ice model [Huwald et al., 2005]. The sea ice component is a more sophisticated model than what has been used in past Arctic ice algae model studies. To validate the model, I simulate an algal bloom at the base of the ice over a season and compare with data from the Resolute area in the Canadian Arctic Archipelago. Results suggest that bloom dynamics are strongly related to the ice growth/melt rate, with ice melt being the trigger for bloom decline. Being able to accurately model physical conditions is essential before ice algae can be accurately modeled, and some recommendations for improvement are discussed. / Les algues de glace sont une composante importante du cycle du carbone dans l'Arctique et peuvent avoir une influence sur le climat. J'explore la dynamique d'une efflorescence algale par le couplage d'un modèle d'algues /nutriments [Lavoie et al., 2005] avec un modèle de la thermodynamique de la glace de mer [Huwald et al., 2005]. Le modèle de glace est plus sophistiqué que les modèles utilisés dans les études précédentes sur les algues de glace Arctiques. Pour valider le modèle, je simule une efflorescence algale à la base de la glace pendant une saison et on compare avec des données de la région de Resolute dans l'Archipel Canadien Arctique. Les résultats indiquent que la dynamique de l'efflorescence est fortement liée au taux de croissance et de fonte, la fonte étant le déclencheur du déclin de l'efflorescence. La capacité de modéliser avec précision les conditions physiques est essentiel pour simuler correctement les algues de glaces, et quelques améliorations sont suggérées.

In this thesis an analysis of the variability of the liquid freshwater (FW) export from the Arctic Ocean on annual and seasonal timescales is presented. Due to missing long-term observations, the variability of the liquid FW export is not well known or understood. Model simulations are therefore currently the only way to study the variability of the FW export from the Arctic. / To investigate the role of the atmospheric forcing for the variability of the liquid FW export, a model simulation for 1950-2007 from the University of Victoria Earth System Climate Model (UVic ESCM) is analyzed. It is shown that large-scale atmospheric circulation changes generally control the variability of the FW export through changes in the FW storage in the Beaufort Gyre. These changes have a large influence on the variability of the FW export through the Canadian Arctic Archipelago (CAA), whereas the Fram Strait FW export is also influenced by changes in the FW storage in the Eurasian basin. / In order to better understand the differences between the mechanisms driving the export variability through Fram Strait and the CAA, passive dye tracers are added to the ocean module of a state-of-the-art global general circulation model, the Community Climate System Model Version 3 (CCSM3). These tracers allow the identification of FW from different sources, and therefore the individual investigation of the export variability of FW from individual sources. It is shown that the Fram Strait FW export is made up mainly of Eurasian runoff and Pacific FW, whereas the FW exported through the CAA comes primarily from Pacific FW and North American runoff. The variability of the FW exports from individual sources is largely in phase in the CAA, as the CAA FW export is mainly driven by velocity anomalies, not FW concentration anomalies. In Fram Strait on the other hand, FW concentration anomalies contribute as much to the FW export variability as velocity anomalies. The variability of the Fram Strait FW concentrations from the two main FW sources is not in phase, as Pacific FW and Eurasian runoff have different pathways to Fram Strait and their variability is governed by different mechanisms. Whereas the Eurasian runoff export depends strongly on the release of FW from the Eurasian shelf during years with an anticyclonic circulation anomaly (negative Vorticity index), the variability of the Pacific export is mainly controlled by changes in the Pacific FW stored in the Beaufort Gyre, with increased export during years with a cyclonic circulation anomaly (positive Vorticity index). A high vertical resolution of the ocean model is found to be important to resolve the role of FW concentration changes for the Fram Strait FW export variability. / The model simulation also shows that in contrast to the interannual variability, the seasonal variability of the Fram Strait FW export is driven almost entirely by the seasonal cycle of sea-ice melt, with a smaller influence of velocity changes or advected FW concentration changes. The disappearance of the summer sea-ice cover in the Arctic during the 21st century might therefore affect the seasonal cycle of the Fram Strait FW export. / Cette thèse de doctorat présente une étude de la variabilité du flux d'eau douce de l'océan Arctique vers l'Atlantique Nord. Parce qu'il existe peu d'observations sur ce flux d'eau douce, sa variabilité n'est pas bien connue. Par conséqeuent, des simulations numériques sont nécessaires pour l'étude du flux d'eau douce. / Premièrement, nous avons utilisé le "Earth Sytem ClimateModel" de l'université de Victoria (UVic ESCM) pour analyser la variabilité du flux d'eau douce pour la période 1950-2007. Nos résultats indiquent que la circulation atmosphérique détermine la variabilité du flux d'eau douce de l'Arctique par son influence sur le tourbillon de Beaufort. Les changements de la circulation dans le tourbillon de Beaufort ont une grande influence sur le flux d'eau douce par l'archipel canadien. Le flux d'eau douce par le détroit de Fram est aussi influencé par les changements de la circulation dans le tourbillon de Beaufort, mais également par les changements de la circulation océanique dans le bassin européen / Par la suite, nous avons implanté des traceurs, représentant les différents types d'eau douce dans l'Arctique, dans un autre modèle climatique: le "Community Climate System Model Version 3" (CCSM3). Ces traceurs nous permettent d'analyser la variabilité du flux d'eau douce en détails. Nous remarquons que la plupart de l'eau douce exportée par le détroit de Fram provient des fleuves eurasiens et de l'océan Pacifique (passant du Pacifique à l'Arctique par le détroit de Béring). Par contraste, le flux d'eau douce par l'archipel canadien est principalement composé d'eau douce provenant de l'océan Pacifique et des fleuves d'Amérique du Nord. Les variablités associées au flux d'eau douce provenant de différentes sources par l'archipel canadien sont en phase. Ceci est dû au fait que la variabilité du flux d'eau douce est controllée par la vitesse de l'eau dans l'archipel canadien. Par contre, la variabilité du flux d'eau douce par le détroit de Fram est controllée par la vitesse de l'eau et aussi par la concentration en eau douce. En outre, le flux d'eau douce provenant de sources différentes ne sont pas en phase dans le détroit de Fram parce que le trajet de l'eau douce provenant des fleuves eurasiens et celui de l'océan Pacifique sont différents et leur variabilités sont controllées par des mécanismes différents. Nous remarquons que le flux d'eau douce provenant des fleuves eurasiens par le détroit de Fram dépend du transport d'eau douce du plateau eurasien. L'eau douce quitte le plateau eurasien quand la circulation atmosphérique est anti-cyclonique (l'index de vorticité est positif). Par contraste, le flux d'eau douce provenant de l'océan Pacifique est plus fort quand le tourbillon de Beaufort est réduit, et cette situation se produit lorsque la circulation atmosphérique est cyclonique (l'index de vorticité est dans ce cas negatif). De plus, nous remarquons qu'une haute résolution spatiale est nécessaire pour représenter la / Par ailleurs, les simulations numériques avec le CCSM3 révèlent que la variabilité saisonnière du flux d'eau douce est déterminée par la fonte de la glace de mer dans le détroit de Fram. La variabilité saisonnière de la vitesse a peu d'influence sur la variabilité saisonnière du flux d'eau douce, et ne fait que retarder d'un mois le maximum et minimum du flux d'eau douce. Parce que la fonte de la glace de mer dans le détroit de Fram détermine la variabilité saisonnière du flux d'eau douce par le détroit de Fram, la disparition de la glace de mer dans l'Arctique en été au cours du 21e siècle (selon les prédictions des modèles climatiques) pourrait changer la variabilité saisonnière du flux d'eau douce.

Western boundary intensification is a major feature of the ocean circulation. It is understood that the change in the Coriolis parameter with latitude is the principal component responsible for this phenomenon. However, traditional arguments for western boundary intensification also call for an interplay between forcing and dissipation. In this thesis, a new geometry, namely the periodic meridional channel, is compared to the closed basin setting in decaying, two-dimensional, beta-plane turbulence. By isolating the effects of the western and eastern boundaries, persistent western boundary intensification is shown to be possible as a balance between nonlinear effects and Rossby wave reflections. The conditions under which this occurs are explained, and particular attention is given to the importance of resolution in resolving the Rossby wave dynamics at the boundaries. A discussion of the regimes that do not exhibit western boundary intensification is also included. / L'intensification des frontières ouest est une caractéristique majeure de la circulation océanique. Il est reconnu que la variation du paramètre de Coriolis avec la latitude est l'élément principalement responsable de ce phénomène. Néanmoins, les arguments traditionnels pour l'intensification des frontières ouest demandent aussi la présence d'une force externe et de dissipation. Dans cette thèse, une nouvelle géométrie, c'est-à-dire le canal méridional périodique, est comparée à l'environnement du bassin fermé dans la turbulence libre sur le plan beta en deux dimensions. L'isolation des effets aux frontières ouest permet l'émergence d'une intensification des frontières ouest persistante en tant qu'équilibre entre les effets non linéaires et la réflexion des ondes de Rossby. Les conditions sous lesquelles ce processus est réalisé sont expliquées, et une attention particulière est portée à l'importance de la résolution pour résoudre la dynamique des ondes de Rossby aux frontières. Une analyse des régimes qui ne manifestent pas d'intensification des frontières ouest est aussi incluse.

The idea that basin-like dynamics may influence or control the Antarctic Circumpolar Current (ACC) is investigated with idealized analytic and numerical models. A simple 2-layer analytic model is developed to predict the transport evolution with the wind stress amplitude. At very low forcing, a non-zero minimum is predicted. This is followed by two distinct dynamical regimes for stronger forcing: a linearly increasing Stommel regime and a saturation regime in which the transport ceases to increase. The vertical distribution of the flow obtained using the geometry of the geostrophic contours (or characteristics) is key to predicting the occurrence of this transport saturation. Many eddy-permitting numerical simulations in large domains are carried over a wide range of parameters. The simulations using a reference zonal wind stress profile agree qualitatively with the analytic model. However, quantitative discrepancies are observed in the saturation regime: i) when a topographic continental ridge is added along the western boundary and ii) when the bottom drag is varied. When a continental ridge is added, eddy fluxes associated with zonal jets enhance the bottom layer recirculation and lower the saturation transport values. When the bottom drag is increased, the lower layer recirculation is suppressed, and this increases the saturation transport values. Experiments investigating the relative roles of the wind stress and wind stress curl in Drake Passage latitudes are also carried out. It is found that the transport is increased when adding a significant constant wind stress. In this regime dominated by the wind stress itself, there is an offset between the numerical results and what is predicted by the analytic model. The vertical momentum flux by mesoscale eddies can be used to distinguish between different regimes: an upward momentum transfer is observed when the dynamics is dominated by the wind stress curl and a downward flux is observed when it is not. In the regime where the wind stress curl dominates, Sverdrup circulation applies over most of the domain --- even in absence of meridional barriers. Also in this regime, transport is saturated, as suggested by the analytic model. The analytic model is also generalized to a continuous stratification and numerical experiments varying the vertical resolution are carried out to test its robustness. These simulations show that the 2-layer and 5-layer models give equivalent results when inertial effects are weak. However, in the 5-layer simulations, topographically-driven inertial recirculations blocking Drake Passage reduce the transport when inertial effects are strong. This behavior disapears, however, when realistic topography is used. In this context, the numerical results agree well with the predictions of the analytic model. It is also found that when the wind stress curl dominates, meridional walls play an important role in the dynamics at weak forcing but become less and less important as the forcing increases. / Dans cette thèse, on étudié l'idée qu'une circulation de type bassin peut influencer et contrôler la dynamique du Courant Circumpolaire Antarctique (CCA) à l'aide de modèles analytique et numérique. Dans un premier temps, on développe un modèle analytique simple à deux couches pour estimer l'évolution du transport en fonction de l'amplitude du vent appliqué à la surface. À très faible amplitude, ce modèle prévoit un transport minimum non-nul. Deux régimes dynamique distincts succèdent à ce minimum: un régime de type "Stommel", dans lequel le transport augmente linéairement et un régime de "saturation" dans lequel le transport plafonne. On utilise la géométrie des "contours géostrophiques", pierre angulaire de la théorie, pour obtenir la distribution verticale de la circulation et estimer l'occurrence de ce régime de saturation.On effectue ensuite un grand nombre de simulations numériques à haute résolution spatiale, en variant la plupart des paramètres du modèle afin de tester la théorie analytique. On définit un profil de vent "référence" soufflant vers l'est suivant une fonction sin^2(y), où y est la latitude. Les simulations utilisant ce profil de référence correspondent qualitativement aux prévisions de la théorie analytique. Par contre, on observe des différences quantitatives dans le régime saturation: i) lorsqu'un plateau continental est ajouté à la frontière ouest et ii) lorsqu'on varie le coefficient de friction au fond. Lorsqu'on ajoute un plateau continental, les flux de tourbillons associés aux jets longitudinaux favorisent la circulation abyssale et baissent ainsi les valeurs de saturation du transport. Lorsqu'on augmente le coefficient de friction, la circulation abyssale est supprimée, ce qui augmente les valeurs de saturation du transport.Dans les expériences où l'on rajoute un vent constant au profil de référence, un décalage est observé entre les résultats des simulations numériques et les prévisions du modèle analytique. Cela définit un nouveau régime où le vent lui-même est fort en comparaison à son rotationnel. Le flux vertical de quantité de mouvement des tourbillons méso-échelle peut être utilisé afin de distinguer les différents régimes. En effet, ce flux est orienté vers le haut lorsque la dynamique est dominée par le rotationnel tandis qu'il est vers le bas lorsque le vent lui-même domine la dynamique. Dans le régime dominé par le rotationnel, une circulation de type Sverdrup est observée dans l'ensemble du bassin, même en absence de péninsules. De plus, le transport y est saturé, tel que suggéré par la théorie analytique.Le modèle analytique est ensuite généralisé au cas où la stratification est continue. On effectue des simulations numériques où la résolution verticale est variée afin de tester cette théorie généralisée. Ces simulations montrent que les modèles à 2 et 5 couches donnent des résultats similaires lorsque les effets d'inertie sont faibles. Par contre, d'intenses circulations générées au-dessus de la topographie bloquent le détroit de Drake et réduisent le transport lorsque les effets d'inertie sont importants. Cependant, ce comportement disparaît lorsqu'une topographie plus réaliste est utilisée. Dans ce contexte, les simulations numériques correspondent aux prévisions du modèle analytique. De plus, on observe que dans un régime où le rotationnel est dominant, les péninsules jouent un rôle important dans la dynamique lorsque le vent est faible tandis que leur rôle devient de plus en plus négligeable lorsque le vent augmente.

A Lagrangian dynamic and thermodynamic iceberg drift model was developed, coded and validated against observations. First, the model was used to create a climatology (1979-2000) of iceberg drift in the Southern Ocean. The simulation reproduced the main patterns of motion and the northernmost extent of Antarctic icebergs as inferred from satellite and ship observations. The model was then used to hindcast 29 giant iceberg drift tracks in what was the first study of an iceberg model ability to reproduce the motion of individual icebergs around the Antarctic continent on timescales of years. The shape and timing of twelve of the twenty-nine tracks was successfully modeled with a model error in the 0.9-50% range. In six cases, the shape of the observed drift track was reproduced but the timing was off, and in the remaining eleven simulation the icebergs moved in the wrong direction. The model error was found to be independent of simulation length suggesting that the error was due to inaccuracies in the forcing data rather than in the physics of the model. In particular, model performance deteriorated in coastal areas and in the southern portions of the Weddell and Ross sea, highlighting the need for higher resolution forcing data in these regions. The model accuracy would benefit from a better definition of the Antarctic coastline, a better representation of Katabatic winds off the continent and a forcing ocean model which would include a dynamic and thermodynamic sea-ice component. / Un modèle Lagrangien dynamique-thermodynamique pour la dérive d'icebergs a été développé, codé et validé à l'aide d'observations. Premièrement, nous avons produit, à l'aide du modèle, une climatologie (1979-2000) de la dérive d'icebergs dans l'Océan du Sud. Les principales tendances du mouvement des icebergs simulés sont en accord avec les observations satellitaires et les mesures in-situ. Le modèle simule bien la limite septentrionale des icebergs d'Antarctique. Nous avons ensuite simulé vingt-neuf trajectoires individuelles d'icebergs géants. C'est la première fois qu'une telle étude est menée pour des icebergs observés autour de l'Antarctique et sur une échelle de plusieurs années. Dans douze cas, le tracé et le minutage de la trajectoire observée a été reproduit avec succès (erreur de 0.9-50%). Six simulations avaient des erreurs de temps mais non de trajet et dans les onze simulations restantes, l'iceberg a dérivé dans la mauvaise direction. Il a été établi que l'erreur du modèle était indépendante de la durée de la simulation, suggérant que l'erreur était due au champ de forçage plutôt qu'aux équations physiques du modèle. En particulier, une détérioration de la qualité des résultats a été observée dans les régions côtières et dans les parties sud des mers de Ross et de Weddell; soulignant ainsi le besoin d'améliorer le champ de forçage dans ces régions. D'autres moyens d'augmenter la précision du modèle seraient, entre autre, une meilleure définition de la géographie côtière de l'Antarctique, une meilleure représentation des vents catabatiques et un modèle océanique incluant une composante de glace dynamique et thermodynamique.

In this thesis we investigate the sensitivity of the Arctic and Antarctic sea-ice cover and global ocean properties to wind-stress and radiative forcings from 1500 to 2000. In a first step, the conversion of the granular sea-ice model (GRAN) (Tremblay and Mysak, 1997) from Cartesian to spherical coordinates is presented. The GRAN is coupled to the Earth System Climate Model of the University of Victoria as an application. The sea-ice thickness and drift in the Arctic Ocean agree well with observations. To further validate the model, the results of thermodynamic component are compared with the Surface Heat Budget of the Arctic Ocean (SHEBA) datasets which were collected between autumn 1997 and autumn 1998. The simulated growth rate is larger and the melt rate is smaller than observed. In a second step, the model is used to investigate which forcings had a dominant effect on the sea-ice cover in both polar regions during the Little Ice Age (LIA), defined as the period between 1500 and 1850, and the industrial period (1850 - 2000). Three different reconstructed wind-stress fields which take into account the North Atlantic Oscillation, one general circulation model wind-stress field, and three different radiative forcings are used (i.e., volcanic activity, insolation changes, greenhouse gas changes). The annual surface air temperature anomalies for the Northern Hemisphere, which are used as model validation, show good agreement with reconstructed temperature anomalies. The simulated sea-ice area and volume in the Northern Hemisphere were larger during the LIA as compared to the present. The comparison between wind-driven and radiatively-driven changes shows that both forcings result in equal magnitude changes in the case of ice volume; for ice area, the wind-driven part is twice as large as the radiatively-driven part. The simulations suggest that the main radiative forcing before 1850 was volcanic forcing, whereas after 1850 the greenhouse gas changes dominated the / Nous étudions dans cette thèse de doctorat l'effet du vent et du forçage radiatif sur la glace de mer de l'Arctique et de l'Antarctique et les proprétés de l'eau de mer. Cette étude se consentre essentiellement sur la période comprise entre 1500 et 2000. En premier lieu, la conversion du modèle de glace de mer granulaire (GRAN) des coordonnées cartésiennes aux coordonnées sphériques est présentée. Comme application, le GRAN est couplé au 'Earth System Climate Model' de l'Université de Victoria. La dérive et l'épaisseur de la glace dans l'océan Arctique concordent avec les observations. Le module thermodynamique est validé en comparant ses résultats aux données du 'Surface Heat Budget of the Arctic Ocean' (SHEBA). Ces données ont été prises entre l'automne 1997 et l'automne 1998. La croissance de glace simulée par le modèle est supérieure et la fonte inférieure à celles observées. En deuxième lieu, le modèle est utilisé pour étudier quels forçages ont un effet dominant sur la couverture de glace de mer des deux régions polaires durant la petite ère glacière (PEG), la période comprise entre 1500 et 1850, et la ère industrielle (1850 - 2000). Trois reconstructions différentes du vent modulée par l'oscillation nord-Atlantique, les vents obtenus d'un modèle de circulation général et trois forçages radiatifs distincts (i.e., activité volcanique, changements de l'insolation et changements des concentrations de gaz à effet de serre) sont utilisés pour forcer le modèle. Les anomalies de la temperature de surface annuelle pour l'hémisphère Nord, qui sont utilisées pour valider le modèle, montrent un bon accord avec les reconstructions des anomalies de temperature. Le volume et l'étendue de la couverture de glace simulés dans l'hémisphère Nord sont plus élevés pendant la PEG que les valeurs présentes. Les résultats obtenus montrent que les changements de volume causés soit par le vent soit par les forçages ra

The Eocene-Oligocene boundary marked the transition from a very warm to a much cooler global climate, and the first appearance of large ice sheets in Antarctica. The opening of the Drake Passage and hence the establishment of the Antarctic Circumpolar Current (ACC) has been proposed as a mechanism for the dramatic Southern Hemisphere cooling at this time. Previous modelling studies explicitly addressing the impact of the Drake Passage used models with simple representations of an atmosphere, and did not consider the Panamanian Seaway, which was open at the time, allowing Atlantic and Pacific waters to exchange. Here, a coupled, fully-dynamic climate model is used to simulate the ocean circulation and climate impacts of the Drake Passage as modulated by the presence or absence of the Ismuth of Panama. With the Panama Ismuth, closing the Drake Passage does not lead to a shutdown of northern overturning, as commonly observed in simpler models, and the temperature impact on Antarctica is similar to prior estimates. In this case, Antarctic cooling reflects the thermal isolation of Antarctica, as warm subtropical waters cannot pass the geostrophic barrier of the Drake Passage. Without the Panama Isthmus, a salinity threshold is reached in the North Atlantic, leading to the cessation of northern overturning and dominance of southern overturning. In the case with no Panama Isthmus, the cooling of surface air temperatures over Antarctica caused by the Drake Passage is 3 × larger than the case with a Panama Isthmus. The larger temperature change reflects the additional onset of inter- hemispheric overturning, transferring heat from the Southern Hemisphere to the Northern Hemisphere, as well as an atmospheric radiative feedback driven by changes in cloud distributions. / La frontère Éocène-Oligocène marque la transition d'un climat chaud à un climat globallement plus frais, et l'apparition de larges calottes glacière en Antarctique. L'ouverture du Passage de Drake et donc l' établissement du Courant Circumpolaire Antarctique a ét é proposé comme un mécanisme permettant le refroidissement drastique de l'hémisèphre sud à cette époque. Des expériences passées de modélisation traitant explicitement des impacts du Passage de Drake utilisent des modèles ne représentant l'atmosphère que de façon simple sans considérer la route maritime du Panama qui ouverte à cette époque, permettait des échanges entre l'Atlantique et le Pacifique. Dans cette étude, un modèle couplé et complètement dynamique est utilisé pour simuler les impacts du Passage de Drake sur la circulation océanique ainsi que sur le climat, tout cela modulé par la présence ou l'absence de l'Isthme de Panama. Avec un Isthme de Panama, la fermeture du Passage de Drake ne résulte pas en l'arrêt de la Circulation Méridienne du Nord comme comunément observé avec de modèles simples malgrés le fait que la réaction de température soit similaire aux estimations passées. Dans ce cas-ci le refroidissement est causé par une isolation thermale de l'Antarctique résultant en l'impossibilité des eaux subtropicales chaudes de franchir la barrière géostropique du Passage de Drake. Sans l'Isthme de Panama, un seuil de salinité est atteint dans l'Atlantique Nord conduisant a la céssation de la Circulation M ́eridienne du Nord ainsi que la domination de celle du sud. Dans le cas ou l'Isthme de Panama est enlevé, le refroisdissement de la température atmopshérique de surface sur l'Antarctique résultant de l'ouverture du Passage de Drake est 3 × plus importante que dans le cas ou l'Isthme de Panama est présent. Cette réaction de température plus intense est le résultat de l'enclenchement de la Circulation Méridiennedu Nord qui transfère de la chaleur de l'Hémisphère Sud à l'Hémisphère Nord et d'une réaction radiative de l'atmosphère résultant de changements dans la distribution des nuages.

The Antarctic Circumpolar Current (ACC), long recognized for its importance in world ocean and climate dynamics, is close to purely zonal in many regions. This means the majority of meridional transport of heat or other scalars is largely achieved by transient mesoscale eddies. Many long-running climate models are still unable to resolve the mesoscale; their effects are parameterized, most commonly using the GM90 parameterization which relates mass transport to the slopes of isopycnal layers via an eddy diffusivity. In this project, the ACC was modelled using a two-layer quasi-geostrophic model with idealistic topography. A series of simulations were carried out with a combination of topography heights and forcing strengths. Time-mean statistics were used to fit several forms of the eddy diffusivity. A paradox was found in the cross-jet diffusivity. In a global sense, it increased as the strength of forcing and jet velocities increased, consistent with higher eddy kinetic energy. On the other hand, a local definition was found to be a local minimum in the jet core where the velocity peaked, consistent with research suggesting the ACC cores suppress mixing. In addition, a skew diffusivity was fitted and found to be globally non-zero for the higher topography, corresponding to a net eastwards flux. For the ridge simulations, large-scale variation in the cross and skew eddy diffusivity was noted and attributed to the effect of local bottom topography on baroclinic instability. This motivated a matrix definition of the eddy diffusivity, which included the bottom topography gradient as an additional parameter to estimate along and cross jet mass fluxes. In all forms the cross-jet diffusivity showed a large dependence on the density gradient. / Le rôle du Courant Circumpolaire Antarctique (CCA) dans la dynamique mondiale des océans est reconnu depuis longtemps. L'écoulement du CCA est presque de purement zonal dans de nombreuses régions, ce qui signifie que la majorité du transport traversant les jets est effectuée par les remous méso-échelles. La plupart des modèles climatiques de longue durée sont encore incapables de résoudre la méso-échelle et les effets de ces remous doivent être paramétrés. Actuellement, le modèle GM90, qui estime le transport de masse à partir des pistes isopycnales en utilisant une diffusivité turbulente est le plus fréquemment employé. Dans ce projet, le CCA a été représenté par un modèle quasi-géostrophique de deux couches avec de la topographie idéaliste. Une série de simulations a été réalisée en changeant la hauteur de la topographie et le niveau du forçage. Dans le but d'évaluer plusieurs formes de la diffusivité turbulente, une analyse des statistiques temporelles a été réalisée. Dans un sens global, la diffusivité turbulente normale dans la direction méridienne, s'est augmentée avec le gradient de densité. Une définition locale, cependant, était un minimum dans le noyau de jet, en accord avec des recherches qui suggèrent la suppression du transport traversant le CCA. Une diffusivité turbulente dans la direction du débit moyen, a été également estimée, montrant un flux net vers l'est pour les topographies et les forçages supérieurs. Des répartitions spatiales ont été notées dans les champs des diffusivités turbulentes dans la direction normale et parallèle à l'écoulement moyen. Celles ont motivé une matrice de diffusion turbulente afin d'essayer de prendre en compte l'effet de la topographie sur les flux de masse. Cette forme a été appliquée globalement et localement, et a permis de réduire la variation spatiale observée.