Dans cette thèse, nous présentons un modèle granulaire décrivant la dynamique d'un assemblage de floes rigides de tailles et de formes diverses, soumis aux forces de traînée dues aux courants atmosphérique et océanique. Ce modèle est basé sur les équations des moments linéaire et angulaire pour décrire la dynamique régulière des floes et sur la résolution de problèmes linéaires de complémentarité pour traiter les collisions entre les floes. Entre les collisions, le mouvement d'un floe isolé satisfait la conservation des équations des moments linéaire et angulaire écrites à partir des formulations classiques des traînées dues au vent et à l'océan. Nous décrivons les collisions entre les floes comme des événements instantanés et les traitons avant qu'elles n'entraînent une interpénétration. Cela implique la notion d'impulsion de contact et la mise sous la forme de problèmes linéaires de complémentarité basés sur la condition de Signorini pour la non interpénétration et la loi de Coulomb. La nature du contact est représentée à travers un coefficient de friction et un coefficient de restitution décrivant la perte d'énergie cinétique durant la collision. Dans cette présente version du modèle, le coefficient de restitution est fixé. Le modèle a été validé en utilisant des données obtenues du mouvement de disques de bois évoluant en bassin de test aussi bien qu'en comparant le comportement des floes simulés avec un comportement attendu dans des scénarios classiques de dérive de glace et de collisions entre des solides rigides. Les résultats de simulations comprenant différents assemblages contenant des floes de tailles et de formes variées, soumis à différents scénarios de forçage, sont aussi discutés. Ils montrent tout le potentiel de notre approche sans qu'une analyse détaillée et complète n'ait encore été proposée. / In this thesis, we present a model describing the dynamics of a population of ice floes with arbitrary shapes and sizes, which are exposed to atmospheric and oceanic skin drag. The granular model presented is based on simplified momentum equations for ice floe motion between collisions and on the resolution of linear complementarity problems to deal with ice floe collisions. Between collisions, the motion of an individual ice floe satisfies the linear and angular momentum conservation equations, with classical formula applied to account for atmospheric and oceanic skin drag. To deal with collisions, before they lead to interpenetration, we included a linear complementarity problem based on the Signorini condition and Coulombs law. The nature of the contact is described through a constant coefficient of friction, as well as a coefficient of restitution describing the loss of kinetic energy during the collision. In the present version of our model, this coefficient is fixed. The model was validated using data obtained from the motion of interacting artificial wood floes in a test basin. The results of simulations comprising few hundreds of ice floes of various shapes and sizes, exposed to different forcing scenarios, and under different configurations, are also discussed. They show that the progressive clustering of ice floes as the result of kinetic energy dissipation during collisions is well captured, and suggest a collisional regimes of floe dispersion at small scales, different from a large-scale regime essentially driven by wind forcing.

A quase inacessibilidade de grandes partes do Oceano Austral torna o conhecimento da espessura do gelo marinho limitado. Esta informação é essencial para a determinação do balanço de massa deste componente da criosfera. Na transição do inverno para a primavera de 2006, uma equipe de pesquisadores, coletou uma série de perfis de espessura de gelo marinho no norte e noroeste do mar de Weddell. Eles estavam a bordo do navio de pesquisa alemão Polarstern do Alfred-Wegener-Institut, Helmholtz-Zentrum für Polar- und Meeresforschung (AWI) e utilizaram um sistema composto por um sensor eletromagnético, altímetro laser e um Sistema de Posicionamento Global Diferencial (DGPS) aerotransportado por helicóptero. Simultaneamente, a Agência Espacial Europeia (ESA), adquiriu imagens ENVISAT ASAR WSM da área de pesquisa. O objetivo principal desta tese é determinar a relação existente entre o retroespalhamento medido pelo Radar de Abertura Sintética e a espessura do gelo marinho obtida pelo HEM (Helicopter-borne ElectroMagnetic system) usando os dados citados acima. Utilizamos os programas de computador NEST®, MATLAB®, EXCEL®, ArcGIS®/ArcMAP®, para processar, analisar e selecionar as imagens, para determinar a relação entre o retroespalhamento e as medidas, quase-tempo-coincidentes, de espessura do gelo. Projetamos as trajetórias dos voos sobre as imagens obtidas nas mesmas datas e extraímos os pixels referentes aos locais onde foram medidas as espessuras de gelo marinho. Apropriamos os valores de espessura do gelo para a área de cada pixel sobre o qual se referiam. Uma análise estatística determinou que o parâmetro que melhor representa a espessura do gelo dentro da área do pixel é a média. A regressão linear é a melhor forma de ajuste das relações entre o valor de retroespalhamento do pixel e a espessura do gelo marinho contido na área deste pixel. O coeficiente de correlação linear de Pearson, resultante de análise paramétrica, indica uma forte correlação (0,75) entre retroespalhamento e espessura do gelo marinho. Porém, a análise não paramétrica de Spearman resultou em um coeficiente de correlação baixo (0,06) o que pode indicar que os dados analisados são compostos por duas populações distintas (e.g., gelo de primeiro ano e plurianual). Contudo, a análise não paramétrica de Kolmogorov-Smirnov aventa a possibilidade de que não tenhamos amostrado toda população. Esta seria a razão de não haver valores de espessuras relativas ao intervalo entre -9.21dB e -1.35dB, o que poderia induzir ao baixo valor do coeficiente de correlação na análise de Spearman. Aplicamos a equação linear: y=0,6345x+12,015 às imagens e pudemos separá-las em doze classes: uma para água e onze para gelo marinho, estas com intervalos de um metro. Isto possibilitou a observação de importantes feições oceanográficas como: canais de águas abertas; cristas de compressão, decaimento do gelo marinho, deslocamento de icebergs (como indicadores de correntes, marés e ventos), liberação de gelo por geleiras de maré, desprendimento de iceberg da plataforma de gelo Larsen C, esteiras de ondas e vórtices oceânicos. Com base no exposto, podemos afirmar que, estatisticamente, nossos resultados são robustos e significantes, com nível de confiança entre 95% e 99%. A equação que propomos é um primeiro passo para inferir-se a espessura do gelo marinho a partir de coeficientes de retroespalhamento SAR. / The almost inaccessibility of large parts of the Southern Ocean makes the knowledge on the sea-ice thickness limited. This information is essential for determining the mass balance of this cryosphere component. During the transition from winter to spring 2006, a researcher team performed several sea-ice thickness profiles in the north and northwest of the Weddell Sea. They were aboard the German research vessel Polarstern from the Alfred-Wegener-Institute, Helmholtz-Zentrum für Polar- und Meeresforschung (AWI) and utilizing for those measurements a Helicopter-borne Electromagnetic system (HEM). The HEM is composed of an electromagnetic sensor, laser altimeter, and a Differential Global Positioning System (DGPS). Simultaneously, the European Space Agency (ESA) ENVISAT ASAR WSM acquired images from the research area. The main objective of this thesis is to determine the relationship between the Synthetic Aperture Radar backscattering and the HEM’s sea-ice thickness measurements using data listed above. We used the computer programs: NEST®, MATLAB®, EXCEL®, ArcGIS® / ArcMap®, to process, analyze and select images, in order to determine the relationship between the backscattering and the quasi-time-coincident ice thickness measurements. We projected the flights trajectories on the images obtained on the same dates, extracting the pixels pertaining to the places where the sea-ice thicknesses were determined. We ascribed ice thickness values for each area covered by the pixel. A statistical analysis determined that the best ice thickness parameter within a pixel area is its mean. Linear regression is the best way to adjust the relationship between the pixel backscatter value and the sea-ice thickness contained within the pixel area. The Pearson linear correlation coefficient, resulting from parametric analysis, indicates a strong correlation (0.75) between backscatter and sea-ice thickness. However, the nonparametric Spearman analysis resulted in a low correlation coefficient (0.06), which may indicate that the analyzed data consist of two distinct populations (e.g., first-year and multi-year ice). However, the Kolmogorov-Smirnov nonparametric analysis brought up the possibility that we just have not sampled the entire population. This could explain the no existence of sea-ice thicknesses values on the interval from -9.21dB to -1.35dB, which could lead to the low correlation coefficient in the Spearman analysis. We applied the linear equation: =,+, to the images, separating them into twelve classes: one for water and eleven to sea-ice, the latter in one-meter thickness intervals. The results enabled the observation of important oceanographic features such as open water channels, pressure ridges, sea-ice decay, icebergs motion (as indicators of currents, tides and winds), glaciers discharge, iceberg calving from Larsen C ice shelf, wakes and oceanic eddies. Based on the foregoing, we can say that, statistically, our results are robust and significantly, with a confidence level from 95% to 99%. The proposed equation is a first step to inferring sea-ice thickness from SAR backscatter coefficients.

A quase inacessibilidade de grandes partes do Oceano Austral torna o conhecimento da espessura do gelo marinho limitado. Esta informação é essencial para a determinação do balanço de massa deste componente da criosfera. Na transição do inverno para a primavera de 2006, uma equipe de pesquisadores, coletou uma série de perfis de espessura de gelo marinho no norte e noroeste do mar de Weddell. Eles estavam a bordo do navio de pesquisa alemão Polarstern do Alfred-Wegener-Institut, Helmholtz-Zentrum für Polar- und Meeresforschung (AWI) e utilizaram um sistema composto por um sensor eletromagnético, altímetro laser e um Sistema de Posicionamento Global Diferencial (DGPS) aerotransportado por helicóptero. Simultaneamente, a Agência Espacial Europeia (ESA), adquiriu imagens ENVISAT ASAR WSM da área de pesquisa. O objetivo principal desta tese é determinar a relação existente entre o retroespalhamento medido pelo Radar de Abertura Sintética e a espessura do gelo marinho obtida pelo HEM (Helicopter-borne ElectroMagnetic system) usando os dados citados acima. Utilizamos os programas de computador NEST®, MATLAB®, EXCEL®, ArcGIS®/ArcMAP®, para processar, analisar e selecionar as imagens, para determinar a relação entre o retroespalhamento e as medidas, quase-tempo-coincidentes, de espessura do gelo. Projetamos as trajetórias dos voos sobre as imagens obtidas nas mesmas datas e extraímos os pixels referentes aos locais onde foram medidas as espessuras de gelo marinho. Apropriamos os valores de espessura do gelo para a área de cada pixel sobre o qual se referiam. Uma análise estatística determinou que o parâmetro que melhor representa a espessura do gelo dentro da área do pixel é a média. A regressão linear é a melhor forma de ajuste das relações entre o valor de retroespalhamento do pixel e a espessura do gelo marinho contido na área deste pixel. O coeficiente de correlação linear de Pearson, resultante de análise paramétrica, indica uma forte correlação (0,75) entre retroespalhamento e espessura do gelo marinho. Porém, a análise não paramétrica de Spearman resultou em um coeficiente de correlação baixo (0,06) o que pode indicar que os dados analisados são compostos por duas populações distintas (e.g., gelo de primeiro ano e plurianual). Contudo, a análise não paramétrica de Kolmogorov-Smirnov aventa a possibilidade de que não tenhamos amostrado toda população. Esta seria a razão de não haver valores de espessuras relativas ao intervalo entre -9.21dB e -1.35dB, o que poderia induzir ao baixo valor do coeficiente de correlação na análise de Spearman. Aplicamos a equação linear: y=0,6345x+12,015 às imagens e pudemos separá-las em doze classes: uma para água e onze para gelo marinho, estas com intervalos de um metro. Isto possibilitou a observação de importantes feições oceanográficas como: canais de águas abertas; cristas de compressão, decaimento do gelo marinho, deslocamento de icebergs (como indicadores de correntes, marés e ventos), liberação de gelo por geleiras de maré, desprendimento de iceberg da plataforma de gelo Larsen C, esteiras de ondas e vórtices oceânicos. Com base no exposto, podemos afirmar que, estatisticamente, nossos resultados são robustos e significantes, com nível de confiança entre 95% e 99%. A equação que propomos é um primeiro passo para inferir-se a espessura do gelo marinho a partir de coeficientes de retroespalhamento SAR. / The almost inaccessibility of large parts of the Southern Ocean makes the knowledge on the sea-ice thickness limited. This information is essential for determining the mass balance of this cryosphere component. During the transition from winter to spring 2006, a researcher team performed several sea-ice thickness profiles in the north and northwest of the Weddell Sea. They were aboard the German research vessel Polarstern from the Alfred-Wegener-Institute, Helmholtz-Zentrum für Polar- und Meeresforschung (AWI) and utilizing for those measurements a Helicopter-borne Electromagnetic system (HEM). The HEM is composed of an electromagnetic sensor, laser altimeter, and a Differential Global Positioning System (DGPS). Simultaneously, the European Space Agency (ESA) ENVISAT ASAR WSM acquired images from the research area. The main objective of this thesis is to determine the relationship between the Synthetic Aperture Radar backscattering and the HEM’s sea-ice thickness measurements using data listed above. We used the computer programs: NEST®, MATLAB®, EXCEL®, ArcGIS® / ArcMap®, to process, analyze and select images, in order to determine the relationship between the backscattering and the quasi-time-coincident ice thickness measurements. We projected the flights trajectories on the images obtained on the same dates, extracting the pixels pertaining to the places where the sea-ice thicknesses were determined. We ascribed ice thickness values for each area covered by the pixel. A statistical analysis determined that the best ice thickness parameter within a pixel area is its mean. Linear regression is the best way to adjust the relationship between the pixel backscatter value and the sea-ice thickness contained within the pixel area. The Pearson linear correlation coefficient, resulting from parametric analysis, indicates a strong correlation (0.75) between backscatter and sea-ice thickness. However, the nonparametric Spearman analysis resulted in a low correlation coefficient (0.06), which may indicate that the analyzed data consist of two distinct populations (e.g., first-year and multi-year ice). However, the Kolmogorov-Smirnov nonparametric analysis brought up the possibility that we just have not sampled the entire population. This could explain the no existence of sea-ice thicknesses values on the interval from -9.21dB to -1.35dB, which could lead to the low correlation coefficient in the Spearman analysis. We applied the linear equation: =,+, to the images, separating them into twelve classes: one for water and eleven to sea-ice, the latter in one-meter thickness intervals. The results enabled the observation of important oceanographic features such as open water channels, pressure ridges, sea-ice decay, icebergs motion (as indicators of currents, tides and winds), glaciers discharge, iceberg calving from Larsen C ice shelf, wakes and oceanic eddies. Based on the foregoing, we can say that, statistically, our results are robust and significantly, with a confidence level from 95% to 99%. The proposed equation is a first step to inferring sea-ice thickness from SAR backscatter coefficients.

A complexa interação que ocorre entre os processos oceânicos e atmosféricos no Oceano Austral afeta a circulação oceânica global em diferentes camadas. O Mar de Weddell e o Mar de Ross possuem reconhecida importância na formação da Água de Fundo Antártica (AABW). O objetivo principal deste trabalho é caracterizar as massas de água dos Mares de Weddell e Ross, através dos dados mais recentes de reanálise oceânica SODA (Simple Ocean Data Assimilation). Através da técnica de separação de massas de água Análise Multiparamétrica Ótima (AMO) foi possível a identificação de 3 principais massas de água no Mar de Ross: Água Profunda Circumpolar Superior (UCDW), Água Profunda Circumpolar Inferior (LCDW) e Água de Plataforma de Baixa Salinidade (LSSW). A UCDW foi a que apresentou a maior variabilidade, não atingindo a Plataforma de gelo do MR durante os anos de 1950-1974. No Mar de Weddell foi possível a identificação das seguintes massas de água: Água Profunda Cálida (WDW), Água Profunda do Mar de Weddell (WSDW) e Água de Fundo do Mar de Weddell (WSBW). A WDW atingiu valores >70% à 800m. A WSDW possui em seu núcleo valores > 90% entre 2000 e 3500m. A WSBW, apresenta ~100% em profundidades > 4000m. / The complex interaction that occurs between the oceanic and atmospheric processes in the Southern Ocean affects global ocean circulation in different layers. The Weddell and Ross Seas have recognized importance in the formation of Antarctic Bottom Water (AABW). This work aims to characterize the water masses of the Weddell and Ross Seas, using the latest ocean data reanalysis SODA (Simple Ocean Data Assimilation). Through the water masses separation technique, Optimum Multiparameter Analysis (OMP), it was possible to identify three main water masses in Ross Sea: Upper Circumpolar Deep Water (UCDW), Lower Circumpolar Deep Water (LCDW) and Low Salinity Shelf Water (LSSW). UCDW showed the greatest variability, not reaching the Ross Sea Ice Shelf during the years 1950-1974. It was possible to identify the following water masses in Weddell Sea: Warm Deep Water (WDW), Weddell Sea Deep Water (WSDW) and Weddell Sea Bottom Water (WSBW). WDW reached values up to 70% in 800m. WSDW has in its core values > 90% between 2000 and 3500m. WSBW presents a contribution up to 100% at depths > 4000m.

Sea ice plays a key role in the global climate system. Indeed, through the albedo effect it reflects significant solar radiation away from the oceans, while it also plays a key role in the momentum and heat transfer between the atmosphere and ocean by acting as an insulating layer between the two. Furthermore, as more sea ice melts due to climate change, additional fresh water is released into the upper oceans, affecting the global circulation of the ocean as a whole. While there has been significant effort in recent decades, the ability to simulate sea ice has lagged behind other components of the climate system and most Earth System Models fail to capture the observed losses of Arctic sea ice, which is largely attributed to our inability to resolve sea ice dynamics. The most widely accepted model for sea ice dynamics is the Viscous- Plastic (VP) rheology, which leads to a very non-linear set of partial differential equations that are known to be intrinsically difficult to solve numerically. This work builds on recent advances in solving these equations with a Jacobian-Free Newton- Krylov (JFNK) solver. We present an improved JFNK solver, where a fully second order discretization is achieved via the Crank Nicolson scheme and consistency is improved via a novel approach to the rheology term. More importantly, we present a significant improvement to the Jacobian approximation used in the Newton iterations, and partially form the action of the matrix by expressing the linear and nearly linear terms in closed form and approximating the remaining highly non-linear term with a second order approximation of its Gateaux derivative. This is in contrast with the previous approach which used a first order approximation for the Gateaux derivative of the whole functional. Numerical tests on synthetic equations confirm the theoretical convergence rate and demonstrate the drastic improvements seen by using a second order approximation in the Gateaux derivative. To produce a fast and efficient solver for VP sea ice dynamics, the improved JFNK solver is then coupled with a non- oscillatory, central differencing scheme for transporting sea ice as well as a novel method for tracking the ice domain using a level set method. Two idealized test cases are then presented and simulation results discussed, demonstrating the solver’s ability to efficiently produce Viscous-Plastic, physically motivated solutions. / Graduate

Le réchauffement régional de l'ouest de la Péninsule Antarctique (WAP) combiné à la diminution attendue de glace de mer et à l'apparition printanière du trou d'ozone pourrait modifier la composition et la structure de la communauté microbienne. De plus, ces variations environnementales pourraient modifier le potentiel de la WAP en tant que puits de CO2. Dans ce contexte, cette thèse visait à évaluer les effets combinés du changement climatique sur la production primaire et sur la composition et la structure de la communauté microbienne de la WAP. Cette thèse visait également à évaluer le rôle de la structure, de la composition, de la production primaire et de la respiration de la communauté microbienne sur les échanges de CO2 entre l'atmosphère et l'océan. Cette étude a premièrement permis de décrire les variations de l'étendue de glace de mer, de l'épaisseur de la couche d'ozone et de la température de surface de l'eau dans la WAP au cours des 30 dernières années (1972-2007) et notamment d'observer le retrait de plus en plus précoce de la glace de mer en relation avec le réchauffement des eaux de la WAP. L'évolution de ces paramètres environnementaux offre une nouvelle fenêtre temporelle de production primaire. Ainsi, cette étude a permis de montrer que la production primaire annuelle a augmenté de 1997 à 2007, et ceci, en relation avec l'anomalie de glace de l'hiver précédent. En effet, la production primaire journalière était négativement et positivement corrélée avec, respectivement, l'étendue de glace de mer et la température de l'eau de septembre à novembre et de février à mars, suggérant que le réchauffement régional de la WAP favorise plus de production primaire durant le printemps et l'automne. En revanche, le retrait précoce de la glace de mer en coïncidence avec l'apparition printanière du trou d'ozone a provoqué l'augmentation de la photoinhibition au printemps (avec 11,6 ± 2,8 % de la production primaire journalière en moyenne). En conséquence, le changement climatique régional de la WAP a, à la fois, un effet positif et un effet négatif sur la production primaire. Cette étude a également permis de décrire la dynamique de la communauté microbienne marine dans l'archipel de Melchior (dans la WAP) de l'automne au printemps 2006. En raison des conditions environnementales extrêmes, l'abondance et la biomasse de la communauté microbienne étaient faibles durant l'automne et l'hiver et dominées par les petites cellules (< 2 µm) et donc par un réseau trophique microbien. En effet, la biomasse phytoplanctonique était faible durant l'automne et l'hiver (avec une concentration moyenne en chlorophylle a, Chl-a, de 0,3 et 0,13 µg l-1, respectivement). La biomasse phytoplanctonique a augmenté au printemps (avec un maximum de Chl-a de 1,13 µg l-1), mais, en dépit des conditions de croissance favorables, est restée faible et le phytoplancton était toujours dominé par de petites cellules (2-20 µm) et donc par le réseau trophique microbien ou multivore. De plus, la disparition précoce de glace de mer durant le printemps 2006 a exposé les eaux de la WAP à de fortes radiations ultraviolettes B (RUVB, 280-320 nm), qui ont eu un effet négatif sur la communauté microbienne des eaux de surface. Cette étude a également mis en évidence la relation existante entre les échanges CO2 et d'O2 entre l'atmosphère et l'océan dans la WAP et la biomasse, la composition, la production primaire et la respiration de la communauté microbienne. Il existait tout d'abord une relation positive entre la concentration en Chl-a et la proportion de diatomées dans la communauté phytoplanctonique. De plus, il existait une corrélation négative significative entre la Chl-a et le ΔpCO2. La production primaire nette de la communauté (NCP) était principalement contrôlée par la production primaire et était négativement et positivement reliée avec le ΔpCO2 et le pourcentage de saturation de l'O2, respectivement, suggérant que la production primaire joue un rôle majeur dans les échanges de CO2 et d'O2 entre l'atmosphère et l'océan dans la WAP. Par ailleurs, le ΔpCO2 moyen au cours des trois années étudiées était de -20,04 ± 44,3 µatm, menant à un puits de CO2 potentiel durant l'été et l'automne dans la région. Le sud de la WAP était un puits potentiel de CO2 (-43,60 ± 39,06 µatm) durant l'automne alors que le nord de la WAP était principalement une source potentielle de CO2 durant l'été ou l'automne (-4,96 ± 37,6 et 21,71 ± 22,39 µatm, respectivement). Les plus fortes concentrations en Chl-a mesurées dans le sud de la WAP pourraient expliquer cette distribution spatiale. / Regional warming in the western Antarctic Peninsula (WAP), along with the expected decrease in sea-ice cover and the seasonal ozone layer breakdown could modify the composition and the structure of the microbial community. In addition, these environmental changes could modify the potential of the WAP as a CO2 sink. In this context, this thesis aimed at evaluating the combined effects of regional climatic changes on the primary production and the composition and structure of the microbial community in the WAP. In a second time, this thesis aimed at evaluating the role of the microbial community structure, composition, primary production and respiration on air-sea CO2 gas exchanges.First, the variations in sea-ice cover, stratospheric ozone layer thickness and sea surface temperature over the last 30 years (1972-2007) were described. Related to the warming of WAP waters, the retreat of sea-ice was happening earlier each decade in the WAP. The observed changes in these environmental parameters offer a new temporal window for primary production. Indeed, the annual primary production increased from 1997 to 2007, in relation with the sea-ice cover anomaly for the previous winter. In addition, daily primary production was negatively and positively correlated to, respectively, sea-ice cover and sea-water temperature from September to November and from February to March, suggesting that regional warming favoured more primary production during spring and fall. On the contrary, the early retreat of sea-ice in spring, in coincidence with the spring ozone layer breakdown, led to an increase in photoinhibition (with an average of 11.6 ± 2.8 % of the daily primary production being photoinhibited). Therefore, regional climatic changes in the WAP had both a positive and a negative impact on primary production.The microbial community variability was also described in the Melchior Archipelago (in the WAP) from fall to spring 2006. Because of the extreme environmental conditions, the microbial community abundance and biomass were low in fall and winter and the community was dominated by small cells (< 2 µm), hence by a microbial food-web. Indeed, phytoplanktonic biomass was low during fall and winter (with respective chlorophyll a concentration, Chl-a, of 0.3 and 0.13 µg l-1). Phytoplankton biomass increased in spring (with a maximum Chl-a of 1.13 µg l-1) but, despite favourable growth conditions, phytoplankton was still dominated by small cells (2-20 µm), hence by a microbial or multivorous food-web. In addition, the early retreat of sea-ice in the spring 2006 exposed the WAP waters to strong ultraviolet B radiations (UVBR, 280-320 nm) that had a negative impact on the microbial community in surface waters.Finally, the relationship between air-sea CO2 and O2 exchanges in the WAP with the phytoplankton community biomass and composition and with the microbial community primary production and respiration was described. A positive relationship existed between Chl-a and the proportion of diatoms in the phytoplankton community. In addition, a negative relationship existed between Chl-a and ΔpCO2. The net community production (NCP) was mainly controlled by primary production and was negatively and positively related to ΔpCO2 and the %O2 saturation, respectively, suggesting that primary production was the main driver of air-sea CO2 and O2 gas exchanges in the WAP. In addition, the average ΔpCO2 for the summers and falls 2002 to 2004 was -20.04 ± 44.3 µatm, leading to a potential CO2 sink during this period in the WAP. The southern WAP was a potential CO2 sink (-43.60 ± 39.06 µatm) during fall while the northern part of the Peninsula was mainly a potential CO2 source during summer and fall (-4.96 ± 37.6 and 21.71 ± 22.39 µatm, respectively). The higher Chl-a concentrations measured in the southern WAP may explain this spatial distribution.

Thesis: S.M., Joint Program in Oceanography/Applied Ocean Science and Engineering (Massachusetts Institute of Technology, Department of Earth, Atmospheric, and Planetary Sciences; and the Woods Hole Oceanographic Institution), 2019 / Cataloged from PDF version of thesis. / Includes bibliographical references (pages 53-58). / The Arctic surface air temperature has warmed nearly twice as much as the global mean since the mid-20th century. Arctic sea ice has also been declining rapidly in recent decades. There is still discussion about how much of this Arctic amplification is caused by local factors, such as changes in surface albedo, versus remote factors, such as changes in heat transport from the midlatitudes. This thesis focuses mainly on the role of poleward heat transport on Arctic amplification. Most of the previous studies on this topic have defined ocean heat transport as the zonally averaged ocean heat transport at 65°N or 70°N, which ignores the physical pathways of heat into the Arctic and may include recirculation of heat in the North Atlantic. In this thesis, we define the ocean heat transport as the heat transport across five sections surrounding the Arctic, to create a closed domain in the Arctic. / Previous studies on Arctic amplification have used either a single model run or have compared results from a multi-model ensemble. While the multi-model ensemble approach may potentially average out biases in individual models, the ensemble spread confounds the model differences and the internal climate variability. In this thesis, we investigate the Arctic amplification in the Community Earth System Model version 1 (CESMi) Large Ensemble. The CESMI Large Ensemble includes 40 members that use the same model and external forcing, but different initializations. This simulates different climate trajectories that can occur in a given atmosphere-ocean-land-cryosphere system. We find that CESMI Large Ensemble projects a large increase towards the end of the 21st century in ocean heat transport into the Arctic, and that the increase in ocean heat transport is significantly correlated with Arctic amplification. / The main contributor to the increase in ocean heat transport is the increase across the Barents Sea Opening. The increase in Barents Sea Opening ocean heat transport is highly correlated with the decrease in sea ice in the Barents-Kara Sea region. We propose that this is because the increase in ocean heat transport melts the ice at the sea ice margin, which results in increased surface heat flux from the ocean and further local feedback through decreased surface albedo and increased cloud coverage. We also find that while the changes in atmosphere heat transport into the Arctic circle at 66.5 N are on the same order as the changes in ocean heat transport, they are not correlated with Arctic amplification. / by Laura Elizabeth Fleming. / S.M. / S.M. Joint Program in Oceanography/Applied Ocean Science and Engineering (Massachusetts Institute of Technology, Department of Earth, Atmospheric, and Planetary Sciences; and the Woods Hole Oceanographic Institution)

Antarctic ice shelves control the ice flux from the continent to the ocean. As such, they play a major role in the stability of the ice sheet and its potential contribution to sea level rise, especially in the context of global change. Below some of these ice shelves, marine ice can be found which is a product of the Deep Thermohaline Circulation. Due to its specific genetic process, marine ice has intrinsic physical (grain size, ice fabric, bubble content, ) and chemical (impurities, water stable isotopes) properties, that differ from those of 'meteoric ice' formed on the continent through snow metamorphism or 'sea ice' resulting from sea water freezing at the ocean-atmosphere surface. Until now however, the effect of these specific properties on marine ice rheology is still very poorly understood.<p><p>The principal objective being to include realistic mechanical parameters for marine ice in ice shelf flow models, uniaxial compression experiments have been performed on various types of marine ice samples. Technical developments are an important component of this thesis has they were necessary to equip the laboratory with the appropriate tools (pneumatic rig, automatic ice fabric data handling).<p><p>Results from experimental compression on isotropic marine ice show that it represents the higher boundary for meteoric ice viscosity throughout the whole temperature range, thereby validating Cuffey and Paterson's relationship with an enhancement factor equals to 1.<p><p>Marine ice is however often quite anisotropic, showing elongated crystals and wide single maximum fabric, that should impact its mechanical properties. Experiments on pre-oriented marine ice samples have therefore been carried out combining the study of epsilon_{oct} vs. tau_{oct} with a thorough analysis of microstructural data 'before' and 'after' the experiment. <p><p>Depending on the orientation of the sample in the applied stress field and on the intensity of the latter, anisotropic marine ice can be harder or softer than its isotropic counterpart, with n=4 often observed in Glen's flow law. Associating the experimental geometrical settings to potential natural equivalent, results suggest that anisotropic marine ice would strengthen ice shelf flow in most areas (for a same given temperature), apart from suturing areas between individual ice streams as they merge to form the ice shelf, where it could become weaker than meteoric ice in certain circumstances.<p><p>Finally, preliminary sensitivity studies, using a simple ice shelf model with our experimental parameters of Glen's flow law have allowed us to discuss the potential impact of rift location, rift size and thermal regime in the ice shelf behavior. / Doctorat en Sciences / info:eu-repo/semantics/nonPublished

Thesis (Ph. D.)--Joint Program in Oceanography/Applied Ocean Science and Engineering (Massachusetts Institute of Technology, Dept. of Earth, Atmospheric, and Planetary Sciences; and the Woods Hole Oceanographic Institution), 2006. / Includes bibliographical references. / Two years of moored oceanographic and automatic weather station data which span the winter ice seasons of 2001-2003 within Marguerite Bay on the western Antarctic Peninsula (wAP) shelf were collected as part of the Southern Ocean Global Ocean Ecosystems Dynamics program. In order to characterize the ice environment in the region, a novel methodology is developed for determining ice coverage, draft and velocity from moored upward-looking acoustic Doppler current profiler data. A linear momentum balance shows the importance of internal ice stresses in the observed motion of the ice pack. Strong inertial, not tidal, motions were observed in both the sea ice and upper ocean. Estimates of upward diapycnal fluxes of heat and salt from the Upper Circumpolar Deep Water to the surface mixed layer indicate almost no contribution from double diffusive convection. A one-dimensional vertical mixed layer model adapted for investigation of mixing beneath an ice-covered ocean indicates that the initial wind event, rather than subsequent inertial shear, causes the majority of the mixing. This work points towards episodic wind-forced shear at the base of the mixed layer coupled with static instability from brine rejection due to ice production as a major factor in mixing on the wAP shelf. / by Jason Hyatt. / Ph.D.

The atmospheric boundary layer in the Arctic is essential for the understanding of climate change and improving regional weather prediction. The aim of this study is to investigate to which degree wind speed profiles retrieved in the Arctic agree with well known wind profile concepts and understand which local impact factors influence the wind speed profile. As part of the Nansen Legacy project, scientists from the University Centre in Svalbard and the University of Bergen installed two wind lidars onboard the research vessel “Kronprins Haakon” during the “Winter Process Cruise” in February 2021. Wind speed profiles were collected over a period of two weeks. They were manually classified into three categories based on their shape. The ideally shaped profiles were fitted against the wind profile power law to identify the exponent, α, for use in the Arctic marine boundary layer. α was found to be 4-5 times smaller than the conventionally applied α = 1/7 for profiles retrieved over open water, which was associated with unstable atmospheric conditions. Additionally, α was found to be considerably larger than 1/7 when sea ice was present, which was associated with stable conditions. A dependency on wind speed was also found. These results underline the importance of adjusting the exponent in order to ac- curately model the wind speed in the Arctic marine boundary layer. The results might be important for optimizing potential wind energy production, which is of great interest with the increasing human activ- ity in the Arctic. Reversed profiles (wind speed maxima closest to the surface) were mainly measured over open ocean and during low wind speeds and were speculated to be related to swell conditions. Pro- files containing a maxima in low levels were primarily measured during stable atmospheric conditions when sea ice was present. Future research in Arctic conditions would benefit from extending wind speed measurements to even lower levels and including stability measurements for an even deeper analysis.