Are polynyas self‐sustaining?

Abstract

In this paper, 441 Conductivity Temperature Depth (CTD) casts from the North Water (NOW) Polynya study were used to calculate geostrophic currents between the 10 and 200 dbar surface during April, May and June 1998. Results for April and May indicated a surface intensified southward flow of 10 to 15 cm s^–1 with a small return flow along the Greenland coast in agreement with inferred currents described by Melling et al. (2001) and surface ice drifts found by Wilson et al. (2001). Southward transports at this time were 0.4–0.55 Sv in April and May. In June, however, surface currents diminished markedly: southward transports declined to 0.1–0.35 Sv, coincident with a decrease in directly measured winds over the polynya and in the surface barometric pressure difference between Grise Fjord and the Carey Islands that was used as a surrogate for the local north wind speed. There was no evident decrease in air pressure difference between Resolute and Grise Fjord, indicative of the strength of the north wind over the eastern Arctic in general. The results are consistent with present thinking that the NOW Polynya is primarily a latent heat polynya, forced by dominant north winds. The idea, broached here, is that the polynya creates its own microclimate which sustains the polynya's ice‐free condition after its initial formation. The mechanism is identified by an anomalous low pressure region associated with surface buoyancy flux in the polynya and is pursued through the application of a simple geostrophic adjustment model that suggests two self‐sustaining mechanisms. Firstly, the frontal intrusion of the cold ambient terrestrial air mass drives a significant surface wind that transports frazil ice to the edge of the polynya before it can congeal. Secondly, rotation at these high latitudes restricts the penetration of the front into the polynya, essentially insulating the centre from freezing temperatures.     

Physical properties and isotopic characteristics of landfast sea ice around the North Water (NOW) Polynya region

Abstract

An ice core sampling program was conducted during the North Water (NOW) Polynya Project 1998 Experiment in northern Baffin Bay during April‐May 1998. The physical properties of snow and sea ice as well as the microstructure and stable isotopic composition of first‐year landfast sea ice near the polynya were investigated. The thickness of sea ice at the sampling sites ranged between 147 and 194 cm with thinner snow cover during the period between mid‐April and late May. The ice was characterized as typical first‐year landfast sea ice, being composed of a thin granular ice layer at the top and an underlying columnar ice layer towards the bottom of the ice. The samples obtained at a site closer to the ice edge of the polynya contained a thin granular ice layer originating from frazil ice near the ice bottom. Formation of frazil ice was considered to be caused by turbulent processes induced by winds, waves and currents forced from the polynya and also mixing with water masses produced at the polynya.     

On the boundary-layer flow over a Canadian Archipelago polynya

A numerical model is presented for predicting boundary-layer parameters for flow over a polynya, i.e., an isolated lead of open water that recurs annually at the same geographical location during Arctic winters. As the flow encounters a polynya, it experiences a sudden change from a cold, dry, rough, snow-covered ice surface to a warmer, wetter, and smoother open water surface. The present model includes both the heat balance equation and the vapor conservation equation in addition to the usual mass continuity equation, x-momentum equation, and the turbulent energy equation. In addition, the buoyancy term is added to the energy equation to account for the buoyancy force introduced by the underlying warmer surface. To close the system, the present model uses the Glushko-type mixing length relationship. During the numerical calculations, a fully implicit finite-difference method was used, and stable numerical solutions were obtained for a fetch over 1 km. Comparison of model results with measured results over a polynya in the Canadian Archipelago show good agreement.     

The Arctic as a trigger for glacial terminations

We propose a hypothesis to explain the very abrupt terminations that end most of the glacial episodes. During the last glaciation, the buildup and southerly expansion of large continental ice-sheets in the Northern Hemisphere and extensive cover of sea ice in the N. Pacific and the N. Atlantic imposed a much more zonal climatic circulation system than exists today. We hypothesize that this, in combination with the frigid (dry) polar air led to a significant decrease in freshwater runoff into the Arctic Ocean. In addition the freshwater contribution of the fresher Pacific water was completely eliminated by the emergence of the Bering Strait (sill depth 50 m). As the Arctic freshwater input was depleted, regions of the Arctic Ocean lost surface stability and eventually overturned, bringing warmer deep water to the surface where it melted the overlying sea ice. This upwelled water was quickly cooled and sank as newly formed deep water. For sustained overturn events, such as might have occurred during the peak of very large glacial periods (i.e. the last glacial maximum), the voluminous deep water formed would eventually overflow into the Nordic Seas and North Atlantic necessitating an equally voluminous rate of return flow of warmer surface waters from the North Atlantic thus breaking down the Arctic's zonal isolation, melting the expansive NA sea ice cover and initiating oceanic heating of the atmosphere over the ice-sheets bordering the NA. We suggest that the combined effect of these overturn-induced events in concert with a Milankovitch warming cycle, was sufficient to drive the system to a termination. We elaborate on this proposed sequence of events, using the model for the formation of the Weddell Sea polynya as proposed by Martinson et al. (1981) and various, albeit sparse, data sets from the circum-Arctic region to apply and evaluate this hypothesis to the problem of glacial terminations.     

Heat budget of the upper Arctic Ocean under a warming climate

The heat budget of the upper Arctic Ocean is examined in an ensemble of coupled climate models under idealised increasing CO₂ scenarios. All of the experiments show a strong amplification of surface air temperatures but a smaller increase in sea surface temperature than the rest of the world as heat is lost to the atmosphere as the sea-ice cover is reduced. We carry out a heat budget analysis of the Arctic Ocean in an ensemble of model runs to understand the changes that occur as the Arctic becomes ice free in summer. We find that as sea-ice retreats heat is lost from the ocean surface to the atmosphere contributing to the amplification of Arctic surface temperatures. Furthermore, heat is mixed upwards into the mixed layer as a result of increased upper ocean mixing and there is increased advection of heat into the Arctic as the ice edge retreats. Heat lost from the upper Arctic Ocean to the atmosphere is therefore replenished by mixing of warmer water from below and by increased advection of warm water from lower latitudes. The ocean is therefore able to contribute more to Arctic amplification.     

Time-dependent response of a zonally averaged ocean–atmosphere–sea ice model to Milankovitch forcing

An ocean–atmosphere–sea ice model is developed to explore the time-dependent response of climate to Milankovitch forcing for the time interval 5–3 Myr BP. The ocean component is a zonally averaged model of the circulation in five basins (Arctic, Atlantic, Indian, Pacific, and Southern Oceans). The atmospheric component is a one-dimensional (latitudinal) energy balance model, and the sea-ice component is a thermodynamic model. Two numerical experiments are conducted. The first experiment does not include sea ice and the Arctic Ocean; the second experiment does. Results from the two experiments are used to investigate (1) the response of annual mean surface air and ocean temperatures to Milankovitch forcing, and (2) the role of sea ice in this response. In both experiments, the response of air temperature is dominated by obliquity cycles at most latitudes. On the other hand, the response of ocean temperature varies with latitude and depth. Deep water formed between 45°N and 65°N in the Atlantic Ocean mainly responds to precession. In contrast, deep water formed south of 60°S responds to obliquity when sea ice is not included. Sea ice acts as a time-integrator of summer insolation changes such that annual mean sea-ice conditions mainly respond to obliquity. Thus, in the presence of sea ice, air temperature changes over the sea ice are amplified, and temperature changes in deep water of southern origin are suppressed since water below sea ice is kept near the freezing point.     

The climatic response of the Arctic Ocean to Soviet river diversions

A numerical model is constructed to evaluate the effect of river diversions on the circulation of the Arctic Ocean, including the climatically important response in the extent of sea ice. The ocean model solves the primitive equations of motion in finite-difference form for the irregular geometry of the Arctic Ocean and Greenland/Norwegian Sea, using 110 km horizontal grid spacing and up to 13 unevenly spaced levels in the vertical. Annual mean atmospheric conditions and river discharges are specified from observations. The presence of sea ice is diagnosed on the basis of model ocean temperature; and the effects of sea ice on the surface fluxes of momentum, heat, and salt are included in a simplified way. Lateral exchanges at the southernmost boundary are held near observed values but respond to circulation changes in the Greenland/Norwegian Sea. Three equilibrium solutions are obtained by eighty-year integrations from simple initial conditions: the first with inflow from all rivers, the second with one-third of the inflow diverted from four major rivers (the Ob, Yenesei, Dvina, and Pechora), and the third with total diversion from those rivers. The middle case corresponds to maximal diversions which are either planned or envisioned by the Soviet Union over the next fifty years, whereas the final extreme case is run in the event that model sensitivity is low relative to that of nature.The control integration gives a good simulation of known water masses and currents. In the Central Arctic, for example, the model correctly predicts a strong shallow halocline, a relatively warm intermediate layer of Atlantic origin, and a temperature jump across the deep Lomonosov Ridge. The overall pattern of surface salinity and the margin of the pack ice are also properly simulated.When runoff into the marginal Kara and Barents Seas is diverted, either in part or in full, almost no effect on the halocline results in the Central Arctic. In particular, deep convection does not develop in the Eurasian Basin, the possibility of which was suggested by Aagaard and Coachman (1975). The vertical stability within the two marginal seas is considerably decreased by the total diversion of four rivers, but not to the point of convective overturning. The surface currents in this area change to confine the water with increased salinity to the shelf region. At deeper levels, an increased salinity tongue spreads into the deep basins of the ice-free Greenland/Norwegian Sea, where existing deep convection is slightly enhanced. As a result, there is some additional heat loss from the Atlantic layer before it enters the Central Arctic. The ice extent remains nearly the same as before within the Kara and Barents Seas. In fact, since modified bottom currents over the continental shelf bring in less heat from the Greenland Sea, an increased thickness of sea ice may result there, in spite of reduced vertical stability. These model responses are generally in agreement with those suggested by Micklin (1981) and by Soviet investigations of the effect of river diversions. These annualmean results should be regarded as tentative, pending confirmation by studies which include the seasonal cycles of runoff and atmospheric forcing.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

On the relationship between spatial patterns of sea‐ice type and the mechanisms which create and maintain the North Water (NOW) polynya

Abstract

Polynyas represent polar oceanic areas with anomalous low sea‐ice concentrations. The North Water (NOW) Polynya refers to a region at the northern end of Baffin Bay which encompasses three separate polynyas. This paper examines the spatial patterns of sea‐ice cover within the NOW region during the winter, spring and fall of 1998 in the context of polynya formation and maintenance mechanisms. To accomplish this a sea‐ice classification scheme for RADARSAT‐1 ScanSAR imagery, obtained between 21 January and 7 December 1998, was developed and implemented within a Geographic Information System (GIS).

The results identify a clear and consistent spatial structure of sea‐ice cover throughout the winter, spring and fall of 1998. Temporally, the polynya opened southward along the Canadian coast and westward away from the Greenland coast. Comparison with parallel oceanographic, atmospheric and ice motion studies suggested that the polynya was primarily controlled by a latent heat mechanism with the exception of the west Greenland coast between Whale Sound and Cape York. The underlying mechanism used to explain the polynya's occurrence along this location is delayed ice formation during freeze‐up and a resultant thinner winter ice cover causing earlier spring ablation than surrounding areas. Arguments for oceanic and/or atmospheric sensible heat contributions are made.     

Simulation of the Atlantic meridional overturning circulation in an atmosphere–ocean global coupled model. Part I: a mechanism governing the variability of ocean convection in a preindustrial experiment

A preindustrial climate experiment was conducted with the third version of the CNRM global atmosphere–ocean–sea ice coupled model (CNRM-CM3) for the Intergovernmental Panel on Climate Change Fourth Assessment Report (IPCC AR4). This experiment is used to investigate the main physical processes involved in the variability of the North Atlantic ocean convection and the induced variability of the Atlantic meridional overturning circulation (MOC). Three ocean convection sites are simulated, in the Labrador, Irminger and Greenland–Iceland–Norwegian (GIN) Seas in agreement with observations. A mechanism linking the variability of the Arctic sea ice cover and convection in the GIN Seas is highlighted. Contrary to previous suggested mechanisms, in CNRM-CM3 the latter is not modulated by the variability of freshwater export through Fram Strait. Instead, the variability of convection is mainly driven by the variability of the sea ice edge position in the Greenland Sea. In this area, the surface freshwater balance is dominated by the freshwater input due to the melting of sea ice. The ice edge position is modulated either by northwestward geostrophic current anomalies or by an intensification of northerly winds. In the model, stronger than average northerly winds force simultaneous intense convective events in the Irminger and GIN Seas. Convection interacts with the thermohaline circulation on timescales of 5–10 years, which translates into MOC anomalies propagating southward from the convection sites.     

Impacts of atmosphere–sea ice–ocean interaction on Southern Ocean deep convection in a climate system model

Deep convection in polar oceans plays a critical role in the variability of global climate. In this study, we investigate potential impacts of atmosphere–sea ice–ocean interaction on deep convection in the Southern Ocean (SO) of a climate system model (CSM) by changing sea ice–ocean stress. Sea ice–ocean stress plays a vital role in the horizontal momentum exchange between sea ice and the ocean, and can be parameterized as a function of the turning angle between sea ice and ocean velocity. Observations have shown that the turning angle is closely linked to the sea-ice intrinsic properties, including speed and roughness, and it varies spatially. However, a fixed turning angle, i.e., zero turning angle, is prescribed in most of the state-of-the-art CSMs. Thus, sensitivities of SO deep convection to zero and non-zero turning angles are discussed in this study. We show that the use of a non-zero turning angle weakens open–ocean deep convection and intensifies continental shelf slope convection. Our analyses reveal that a non-zero turning angle first induces offshore movement of sea ice transporting to the open SO, which leads to sea ice decrease in the SO coastal region and increase in the open SO. In the SO coastal region, the enhanced sea-ice divergence intensifies the formation of denser surface water descending along continental shelf by enhanced salt flux and reduced freshwater flux, combined with enhanced Ekman pumping and weakened stratification, contributing to the occurrence and intensification of continental shelf slope convection. On the other hand, the increased sea ice in the open SO weakens the westerlies, enhances sea-level pressure, and increases freshwater flux, whilst oceanic cyclonic circulation slows down, sea surface temperature and sea surface salinity decrease in the open SO response to the atmospheric changes. Thus, weakened cyclonic circulation, along with enhanced freshwater flux, reduced deep–ocean heat content, and increased stability of sea water, dampens the open–ocean deep convection in the SO, which in turn cools the sea surface temperature, increases sea-level pressure, and finally increases sea-ice concentration, providing a positive feedback. In the CSM, the use of a non-zero turning angle has the capability to reduce the SO warm bias. These results highlight the importance of an accurate representation of sea ice–ocean coupling processes in a CSM.     

The role of quasi-stationary flaw polynya in formation of dense shelf waters in the wintertime and their subsequent slope cascading (the Laptev Sea case study)

Shelf areas in the region of the Severnaya Zemlya Archipelago in the Laptev Sea are characterized by existing quasi-stationary flaw polynya that periodically opens throughout the entire wintertime under the action of strong offshore winds, which occur during the passage of cyclones. In periods of the open water surface, a near-surface turbulent layer or forced convection layer is formed in the flaw polynya; the water in the layer formed undergoes intense salinization and its dense increases due to active volumetric frazil ice production. As a result of the gravity force action, intense three-dimensional convective circulation develops in the underlying layers. It leads to a fast convective adjustment of the entire water column, especially, in the late winter, when residual stratification in the area of polynya is weakened with the total action of salinization due to the background static ΣMs _back and periodical local frazil ice formation ΣMs _f . On the whole for the entire winter period ΣMs _f is 3.4 times greater than ΣMs _back, although, during one month, probable lifetime of polynya with open water surface is several days. However, in these periods, salt fluxes with frazil ice production exceed background salt fluxes in the congelation polynya and background salt fluxes under heavy ice (limiting the polynya) 10–80 times. Spreading outside the polynia, dense shelf waters form in the area of polynya mesoscale baroclinic circulation, first generating intense shelf cascading, then intense slope cascading, which is of a local and random character. Some estimates of elements of baroclinic circulation of a convective origin in the area of polynia were obtained from the laboratory modeling results and are confirmed by field observation data.     

Variability of Fram Strait sea ice export: causes, impacts and feedbacks in a coupled climate model

Analyses of a 500-year control integration of the global coupled atmosphere–sea ice–ocean model ECHAM5.0/MPI-OM show a high variability in the ice export through Fram Strait on interannual to decadal timescales. This variability is mainly determined by variations in the sea level pressure gradient across Fram Strait and thus geostrophic wind stress. Ice thickness anomalies, formed at the Siberian coast and in the Chukchi Sea, propagate across the Arctic to Fram Strait and contribute to the variability of the ice export on a timescale of about 9 years. Large anomalies of the ice export through Fram Strait cause fresh water signals, which reach the Labrador Sea after 1–2 years and lead to significant changes in the deep convection. The associated anomalies in ice cover and ocean heat release have a significant impact on air temperature in the Labrador Sea and on the large-scale atmospheric circulation. This affects the sea ice transport and distribution in the Arctic again. Sensitivity studies, simulating the effect of large ice exports through Fram Strait, show that the isolated effect of a prescribed ice/fresh water anomaly is very important for the climate variability in the Labrador Sea. Thus, the ice export through Fram Strait can be used for predictability of Labrador Sea climate up to 2 years in advance.     

Impact of a realistic river routing in coupled ocean–atmosphere simulations of the Last Glacial Maximum climate

The presence of large ice sheets over North America and North Europe at the Last Glacial Maximum (LGM) strongly impacted Northern hemisphere river pathways. Despite the fact that such changes may significantly alter the freshwater input to the ocean, modified surface hydrology has never been accounted for in coupled ocean–atmosphere general circulation model simulations of the LGM climate. To reconstruct the LGM river routing, we use the ICE-5G LGM topography. Because of the uncertainties in the extent of the Fennoscandian ice sheet in the Eastern part of the Kara Sea, we consider two more realistic river routing scenarios. The first scenario is characterised by the presence of an ice dammed lake south of the Fennoscandian ice sheet, and corresponds to the ICE-5G topography. This lake is fed by the Ob and Yenisei rivers. In the second scenario, both these rivers flow directly into the Arctic Ocean, which is more consistent with the latest QUEEN ice sheet margin reconstructions. We study the impact of these changes on the LGM climate as simulated by the IPSL_CM4 model and focus on the overturning thermohaline circulation. A comparison with a classical LGM simulation performed using the same model and modern river basins as designed in the PMIP2 exercise leads to the following conclusions: (1) The discharge into the North Atlantic Ocean is increased by 2,000 m³/s between 38° and 54°N in both simulations that contain LGM river routing, compared to the classical LGM experiment. (2) The ice dammed lake is shown to have a weak impact, relative to the classical simulation, both in terms of climate and ocean circulation. (3) In contrast, the North Atlantic deep convection and meridional overturning are weaker than during the classical LGM run if the Ob and Yenisei rivers flow directly into the Arctic Ocean. The total discharge into the Arctic Ocean is increased by 31,000 m³/s, relative to the classical LGM simulation. Consequentially, northward ocean heat transport is weaker, and sea ice more extensive, in better agreement with existing proxy data.     

一个海洋-海冰耦合模式的北极气候模拟

This paper evaluates the simulation of Arctic sea ice states using an ocean-ice coupled model that employs LASG/IAP(the State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics/the Institute of Atmospheric Physics) Climate Ocean Model(LICOM) and the sea-ice model from the Bergen Climate Model(BCM).It is shown that the coupled model can reasonably reproduce the major characteristics of the mean state,annual cycle,and interannual variability of the Arctic sea ice concentration.The coupled model also shows biases that were generally presented in other models,such as the underestimation of summer sea ice concentration and thickness as well as the unsatisfactory sea ice velocity.Sensitivity experiments indicate that the insufficient performance of the ocean model at high latitudes may be the main reason for the biases in the coupled model.The smoother and the fake "island",which had to be used due to the model’s grid in the North Pole region,likely caused the ocean model’s weak performance.Sea ice model thermodynamics are also responsible for the sea ice simulation biases.Therefore,both the thermodynamic module of the sea ice component and the model grid of the ocean component need to be further improved.     

Thinner Sea Ice Contribution to the Remarkable Polynya Formation North of Greenland in August 2018

In August 2018, a remarkable polynya was observed off the north coast of Greenland, a perennial ice zone where thick sea ice cover persists. In order to investigate the formation process of this polynya, satellite observations, a coupled iceocean model, ocean profiling data, and atmosphere reanalysis data were applied. We found that the thinnest sea ice cover in August since 1978(mean value of 1.1 m, compared to the average value of 2.8 m during 1978-2017) and the modest southerly wind caused by a positive North Atlantic Oscillation(mean value of 0.82, compared to the climatological value of-0.02) were responsible for the formation and maintenance of this polynya. The opening mechanism of this polynya differs from the one formed in February 2018 in the same area caused by persistent anomalously high wind. Sea ice drift patterns have become more responsive to the atmospheric forcing due to thinning of sea ice cover in this region.     

Simulation of sea ice in FGOALS-g2: Climatology and late 20th century changes

Sea ice is an important component in the Earth’s climate system. Coupled climate system models are indispensable tools for the study of sea ice, its internal processes, interaction with other components, and projection of future changes. This paper evaluates the simulation of sea ice by the Flexible Global Ocean-Atmosphere-Land System model Grid-point Version 2 (FGOALS-g2), in the fifth phase of the Coupled Model Inter-comparison Project (CMIP5), with a focus on historical experiments and late 20th century simulation. Through analysis, we find that FGOALS-g2 produces reasonable Arctic and Antarctic sea ice climatology and variability. Sea ice spatial distribution and seasonal change characteristics are well captured. The decrease of Arctic sea ice extent in the late 20th century is reproduced in simulations, although the decrease trend is lower compared with observations. Simulated Antarctic sea ice shows a reasonable distribution and seasonal cycle with high accordance to the amplitude of winter-summer changes. Large improvement is achieved as compared with FGOALS-g1.0 in CMIP3. Diagnosis of atmospheric and oceanic forcing on sea ice reveals several shortcomings and major aspects to improve upon in the future: (1) ocean model improvements to remove the artificial island at the North Pole; (2) higher resolution of the atmosphere model for better simulation of important features such as, among others, the Icelandic Low and westerly wind over the Southern Ocean; and (3) ocean model improvements to accurately receive freshwater input from land, and higher resolution for resolving major water channels in the Canadian Arctic Archipelago.     

A box model study of the Greenland Sea,Norwegian Sea,and Arctic Ocean

We develop a simple dynamical system model of the Arctic Ocean and marginal seas by applying the Martinson, Killworth and Gordon box model of a high-latitude two-layer ocean to four regions connected together: the Greenland Sea, the Norwegian Sea, the Arctic Ocean, and the Greenland Gyre. The latter is a small convective region embedded in the northwest corner of the Norwegian Sea. The model for each region consists of a thermodynamic ice layer that covers two layers of saline water which can, under specific conditions, become statically unstable and hence create a state of active overturning. The system is forced by monthly mean atmospheric temperatures in the four regions, by continental runoffs and by inflows from adjacent oceans. The model predicts the ice thickness, and the temperature and salinity of the water in the upper layer of the four regions. Also determined are the water temperature and salinity of the lower layer in the Arctic Ocean box. The convective state of any given region, i.e. whether it is in an active overturning mode or not, is also determined as a continuous function of time. The different output variables of the model, which are the response to climatological forcing conditions, compare favourably with observed data. In the control run, the Arctic Ocean region is characterized by continuous ice cover, the Greenland Sea and Greenland Gyre have ice cover only during winter, and the Norwegian Sea region never forms an ice cover. Another feature of the control run is the winter time occurrence of convective overturning in the upper 200 m in the Greenland Gyre region. The model is also used for different anomaly experiments: a positive air temperature anomaly which represents a global warming of the earth, a negative salt anomaly in the Norwegian Sea which simulates the great salinity anomaly of the 1960s and 1970s, and an increase in the ice flux through Fram Strait which parameterizes anomalous ice production in the Arctic.     

Coupled climate modelling of ocean circulation changes during ice age inception

Freshening of high latitude surface waters can change the large-scale oceanic transport of heat and salt. Consequently, atmospheric and sea ice perturbations over the deep water production sites excite a large-scale response establishing an oceanic "teleconnection" with time scales of years to centuries. To study these feedbacks, a coupled atmosphere-ocean-sea ice model consisting of a two dimensional atmospheric energy and moisture balance model (EMBM) coupled to a thermodynamic sea ice model and an ocean general circulation model is utilised. The coupled model reproduces many aspects of the present oceanic circulation. We also investigate the climate impact of changes in fresh water balance during an ice age initiation. In this experiment part of the precipitation over continents is stored within continental ice sheets. During the buildup of ice sheets the oceanic stratification in the North Atlantic is weakened by a reduced continental run-off leading to an enhanced thermohaline circulation. Under these conditions salinity is redistributed such that deep water is more saline than under present conditions. Once the ice sheets built up, we simulate an ice age climate without net fresh water storage on the continents. In this case the coupled model reproduces the shallow and weak overturning cell, an ice edge advance insulating the upper ocean, and many other aspects of the glacial circulation.     

Sensitivity of open-water ice growth and ice concentration evolution in a coupled atmosphere-ocean-sea ice model

A coupled atmosphere-ocean-sea ice model is applied to investigate to what degree the area-thickness distribution of new ice formed in open water affects the ice and ocean properties. Two sensitivity experiments are performed which modify the horizontal-to-vertical aspect ratio of open-water ice growth. The resulting changes in the Arctic sea-ice concentration strongly affect the surface albedo, the ocean heat release to the atmosphere, and the sea-ice production. The changes are further amplified through a positive feedback mechanism among the Arctic sea ice, the Atlantic Meridional Overturning Circulation (AMOC), and the surface air temperature in the Arctic, as the Fram Strait sea ice import influences the freshwater budget in the North Atlantic Ocean. Anomalies in sea-ice transport lead to changes in sea surface properties of the North Atlantic and the strength of AMOC. For the Southern Ocean, the most pronounced change is a warming along the Antarctic Circumpolar Current (ACC), owing to the interhemispheric bipolar seasaw linked to AMOC weakening. Another insight of this study lies on the improvement of our climate model. The ocean component FESOM is a newly developed ocean-sea ice model with an unstructured mesh and multi-resolution. We find that the subpolar sea-ice boundary in the Northern Hemisphere can be improved by tuning the process of open-water ice growth, which strongly influences the sea ice concentration in the marginal ice zone, the North Atlantic circulation, salinity and Arctic sea ice volume. Since the distribution of new ice on open water relies on many uncertain parameters and the knowledge of the detailed processes is currently too crude, it is a challenge to implement the processes realistically into models. Based on our sensitivity experiments, we conclude a pronounced uncertainty related to open-water sea ice growth which could significantly affect the climate system sensitivity.     

Atmospheric Response To Spatial Variations In Concentration And Size Of Polynyas In The Southern Ocean Sea-Ice Zone

Although it is well known that sea-ice regions are important components of the Earth's climate system, the exchanges of energy between ocean, ice and atmosphere are not well understood. The majority of past observational and modelling studies of atmosphere-surface interactions over sea-ice regions were primarily concerned with airflow over a single, isolated area of open water. The more realistic situations of multiple polynyas within a sea-ice field and different areal concentrations of sea ice were studied here. Spatial structure of the atmospheric boundary layer in response to this surface was simulated using a high-resolution numerical model. A sea-ice concentration of 80%, typical of the Southern Ocean sea-ice zone, was maintained within a 100-km wide domain. The effects of three polynya characteristics were assessed: their horizontal extent; local concentration of sea ice (LCI); and their arrangement with ice floes. Over polynyas of all sizes distinct plumes of upward heat flux, their width and height closely linked to polynya width, resulted in mixed layers 600 to 1000 m deep over and downwind of the polynyas, their depth increasing with polynya width. Mean surface heat flux (MSHF) increased with size in polynyas less than 30 km wide. The air-to-ice MSHF over the first 10 km of sea-ice downwind of each polynya and the domain-average surface heat flux increased linearly with polynya width. Turbulent kinetic energy plumes occurred over all polynyas, their heights and widths increasing with polynya widths. Downward flux of high momentum air in the plumes caused increased wind speeds over polynyas in the layer from about 300–1000 m above the surface, the depth varying directly with polynya width. MSHFs decreased as LCIs increased. The arrangement of polynyas had relatively little effect on the overall depth of the modified layer but did influence the magnitude and spatial structure of vertical heat transfer. In the two-polynya case the MSHF over the polynyas was larger when they were closer together. Although the MSHF over the sea ice between the polynyas decreased in magnitude as their separation increased, the percentage of the polynya-to-air heat recaptured by this ice floe increased fivefold.