Interactions between wind-blown snow redistribution and melt ponds in a coupled ocean–sea ice model

Introducing a parameterization of the interactions between wind-driven snow depth changes and melt pond evolution allows us to improve large scale models. In this paper we have implemented an explicit melt pond scheme and, for the first time, a wind dependant snow redistribution model and new snow thermophysics into a coupled ocean–sea ice model. The comparison of long-term mean statistics of melt pond fractions against observations demonstrates realistic melt pond cover on average over Arctic sea ice, but a clear underestimation of the pond coverage on the multi-year ice (MYI) of the western Arctic Ocean. The latter shortcoming originates from the concealing effect of persistent snow on forming ponds, impeding their growth. Analyzing a second simulation with intensified snow drift enables the identification of two distinct modes of sensitivity in the melt pond formation process. First, the larger proportion of wind-transported snow that is lost in leads directly curtails the late spring snow volume on sea ice and facilitates the early development of melt ponds on MYI. In contrast, a combination of higher air temperatures and thinner snow prior to the onset of melting sometimes make the snow cover switch to a regime where it melts entirely and rapidly. In the latter situation, seemingly more frequent on first-year ice (FYI), a smaller snow volume directly relates to a reduced melt pond cover. Notwithstanding, changes in snow and water accumulation on seasonal sea ice is naturally limited, which lessens the impacts of wind-blown snow redistribution on FYI, as compared to those on MYI. At the basin scale, the overall increased melt pond cover results in decreased ice volume via the ice-albedo feedback in summer, which is experienced almost exclusively by MYI.

The impact of variable sea ice roughness on changes in Arctic Ocean surface stress: a model study

The Arctic sea ice cover is thinning and retreating, causing changes in surface roughness that in turn modify the momentum flux from the atmosphere through the ice into the ocean. New model simulations comprising variable sea ice drag coefficients for both the air and water interface demonstrate that the heterogeneity in sea ice surface roughness significantly impacts the spatial distribution and trends of ocean surface stress during the last decades. Simulations with constant sea ice drag coefficients as used in most climate models show an increase in annual mean ocean surface stress (0.003 N/m2 per decade, 4.6%) due to the reduction of ice thickness leading to a weakening of the ice and accelerated ice drift. In contrast, with variable drag coefficients our simulations show annual mean ocean surface stress is declining at a rate of -0.002 N/m2 per decade (3.1%) over the period 1980-2013 because of a significant reduction in surface roughness associated with an increasingly thinner and younger sea ice cover. The effectiveness of sea ice in transferring momentum does not only depend on its resistive strength against the wind forcing but is also set by its top and bottom surface roughness varying with ice types and ice conditions. This reveals the need to account for sea ice surface roughness variations in climate simulations in order to correctly represent the implications of sea ice loss under global warming.

The influence of sea ice dynamics on the climate sensitivity and memory to increased Antarctic sea ice

The study analyzes the sensitivity and memory of the Southern Hemisphere coupled climate system to increased Antarctic sea ice (ASI), taking into account the persistence of the sea ice maxima in the current climate. The mechanisms involved in restoring the climate balance under two sets of experiments, which differ in regard to their sea ice models, are discussed. The experiments are perturbed with extremes of ASI and integrated for 10 yr in a large 30-member ensemble. The results show that an ASI maximum is able to persist for ; 4 yr in the current climate, followed by a negative sea ice phase. The sea ice insulating effect during the positive phase reduces heat fluxes south of 60 8 S, while at the same time these are intensified at the sea ice edge. The increased air stability over the sea ice field strengthens the polar cell while the baroclinicity increases at midlatitudes. The mean sea level pressure is reduced (increased) over high latitudes (midlatitudes), typical of the southern annular mode (SAM) positive phase. The Southern Ocean (SO) becomes colder and fresher as the sea ice melts mainly through sea ice lateral melting, the consequence of which is an increase in the ocean stability by buoyancy and mixing changes. The climate sensitivity is triggered by the sea ice insulating process and the resulting freshwater pulse (fast response), while the climate equilibrium is restored by the heat stored in the SO subsurface layers (long response). It is concluded that the time needed for the ASI anomaly to be dissipated and/or melted is shortened by the sea ice dynamical processes.

Variability and trends in Laptev Sea ice outflow between 1992-2011

Variability and trends in seasonal and interannual ice area export out of the Laptev Sea between 1992 and 2011 are investigated using satellite-based sea ice drift and concentration data. We found an average total winter (Octo- ber to May) ice area transport across the northern and east- ern Laptev Sea boundaries (NB and EB) of 3.48 × 10 5 km 2 . The average transport across the NB (2.87 × 10 5 km 2 ) is thereby higher than across the EB (0.61 × 10 5 km 2 ), with a less pronounced seasonal cycle. The total Laptev Sea ice area flux significantly increased over the last decades (0.85 × 10 5 km 2 decade − 1 , p> 0 . 95), dominated by increas- ing export through the EB (0.55 × 10 5 km 2 decade − 1 , p> 0 . 90), while the increase in export across the NB is smaller (0.3 × 10 5 km 2 decade − 1 ) and statistically not significant. The strong coupling between across-boundary SLP gradient and ice drift velocity indicates that monthly variations in ice area flux are primarily controlled by changes in geostrophic wind velocities, although the Laptev Sea ice circulation shows no clear relationship with large-scale atmospheric in- dices. Also there is no evidence of increasing wind velocities that could explain the overall positive trends in ice export. The increased transport rates are rather the consequence of a changing ice cover such as thinning and/or a decrease in con- centration. The use of a back-propagation method revealed that most of the ice that is incorporated into the Transpolar Drift is formed during freeze-up and originates from the cen- tral and western part of the Laptev Sea, while the exchange with the East Siberian Sea is dominated by ice coming from the central and southeastern Laptev Sea. Furthermore, our re- sults imply that years of high ice export in late winter (Febru- ary to May) have a thinning effect on the ice cover, which in turn preconditions the occurence of negative sea ice extent anomalies in summer.

Daily thin-ice thickness maps from modis thermal infrared imagery

Accurate knowledge of ice-production rates within the marginal ice zones of the Arctic Ocean requires monitoring of the thin-ice distribution within polynyas. The thickness of the ice layer controls the heat loss and hence the new-ice formation. An established thinice algorithm using high-resolution MODIS data allows deriving the ice-thickness distribution within polynyas. The average uncertainty is ±4.7 cm for ice thicknesses below 0.2 m. In this study, the ice-thickness distributions within the Laptev Sea polynya for the two winter seasons 2007/08 and 2008/09 are calculated. Then, a new method is applied to determine a daily MODIS thin-ice product.

Spatio-temporal variability of polynya dynamics and ice production in the Laptev Sea between the winters of 1979/80 and 2007/08

Polynyas in the Laptev Sea are examined with respect to recurrence and interannual wintertime ice production.We use a polynya classification method based on passive microwave satellite data to derive daily polynya area from long-term sea-ice concentrations. This provides insight into the spatial and temporal variability of open-water and thin-ice regions on the Laptev Sea Shelf. Using thermal infrared satellite data to derive an empirical thin-ice distribution within the thickness range from 0 to 20 cm, we calculate daily average surface heat loss and the resulting wintertime ice formation within the Laptev Sea polynyas between 1979 and 2008 using reanalysis data supplied by the National Centers for Environmental Prediction, USA, as atmospheric forcing. Results indicate that previous studies significantly overestimate the contribution of polynyas to the ice production in the Laptev Sea. Average wintertime ice production in polynyas amounts to approximately 55 km39 27% and is mostly determined by the polynya area, wind speed and associated large-scale circulation patterns. No trend in ice production could be detected in the period from 1979/80 to 2007/08.

Observations of supercooling and frazil ice formation in the Laptev Sea coastal polynya

This paper examines a hydrographic response to the wind‐driven coastal polynya activity over the southeastern Laptev Sea shelf for April–May 2008, using a combination of Environmental Satellite (Envisat) advanced synthetic aperture radar (ASAR) and TerraSAR‐X satellite imagery, aerial photography, meteorological data, and SBE‐37 salinity‐temperature‐depth and acoustic Doppler current profiler land‐fast ice edgemoored instruments. When ASAR observed the strongest end‐of‐April polynya event with frazil ice formation, the moored instruments showed maximal acoustical scattering within the surface mixed layer, and the seawater temperatures were either at or 0.02°C below freezing. We also find evidence of the persistent horizontal temperature and salinity gradients across the fast ice edge to have the signature of geostrophic flow adjustment as predicted by polynya models.

Impact of a Fram Strait cyclone on ice edge, drift, divergence, and concentration: Possibilities and limits of an observational analysis

This study aims at the determination of a Fram Strait cyclone track and of the cyclone’s impact on ice edge, drift, divergence, and concentration. A 24 h period on 13–14 March 2002 framed by two RADARSAT images is analyzed. Data are included from autonomous ice buoys, a research vessel, Special Sensor Microwave Imager (SSM/I) and QuikSCAT satellite, and the operational European Centre for Medium-Range Weather Forecasts (ECMWF) model. During this 24 h period the cyclone moved northward along the western ice edge in the Fram Strait, crossed the northern ice edge, made a left-turn loop with 150 km diameter over the sea ice, and returned to the northern ice edge. The ECMWF analysis places the cyclone track 100 km too far west over the sea ice, a deviation which is too large for representative sea ice simulations. On the east side of the northward moving cyclone, the ice edge was pushed northward by 55 km because of strong winds. On the rear side, the ice edge advanced toward the open water but by a smaller distance because of weaker winds there. The ice drift pattern as calculated from the ice buoys and the two RADARSAT images is cyclonically curved around the center of the cyclone loop. Ice drift divergence shows a spatial pattern with divergence in the loop center and a zone of convergence around. Ice concentration changes as retrieved from SSM/I data follow the divergence pattern such that sea ice concentration increased in areas of divergence and decreased in areas of convergence.

Impact of instantaneous sea ice removal in a coupled general circulation model

The impact of extreme sea ice initial conditions on modelled climate is analysed for a fully coupled atmosphere ocean sea ice general circulation model, the Hadley Centre climate model HadCM3. A control run is chosen as reference experiment with greenhouse gas concentration fixed at preindustrial conditions. Sensitivity experiments show an almost complete recovery from total removal or strong increase of sea ice after four years. Thus, uncertainties in initial sea ice conditions seem to be unimportant for climate modelling on decadal or longer time scales. When the initial conditions of the ocean mixed layer were adjusted to ice-free conditions, a few substantial differences remained for more than 15 model years. But these differences are clearly smaller than the uncertainty of the HadCM3 run and all the other 19 IPCC fourth assessment report climate model preindustrial runs. It is an important task to improve climate models in simulating the past sea ice variability to enable them to make reliable projections for the 21st century.

Impact of cyclones on sea ice in the Fram Strait: numerical simulations and observations

The sea ice export from the Arctic is of global importance due to its fresh water which influences the oceanic stratification and, thus, the global thermohaline circulation. This study deals with the effect of cyclones on sea ice and sea ice transport in particular on the basis of observations from two field experiments FRAMZY 1999 and FRAMZY 2002 in April 1999 and March 2002 as well as on the basis of simulations with a numerical sea ice model. The simulations realised by a dynamic-thermodynamic sea ice model are forced with 6-hourly atmospheric ECMWF- analyses (European Centre for Medium-Range Weather Forecasts) and 6-hourly oceanic data of a MPI-OM-simulation (Max-Planck-Institute Ocean Model). Comparing the observed and simulated variability of the sea ice drift and of the position of the ice edge shows that the chosen configuration of the model is appropriate for the performed studies. The seven observed cyclones change the position of the ice edge up to 100 km and cause an extensive decrease of sea ice coverage by 2 % up to more than 10 %. The decrease is only simulated by the model if the ocean current is strongly divergent in the centre of the cyclone. The impact is remarkable of the ocean current on divergence and shear deformation of the ice drift. As shown by sensitivity studies the ocean current at a depth of 6 m – the sea ice model is forced with – is mainly responsible for the ascertained differences between simulation and observation. The simulated sea ice transport shows a strong variability on a time scale from hours to days. Local minima occur in the time series of the ice transport during periods with Fram Strait cyclones. These minima are not caused by the local effect of the cyclone’s wind field, but mainly by the large-scale pattern of surface pressure. A displacement of the areas of strongest cyclone activity in the Nordic Seas would considerably influence the ice transport.

On the parameterization of turbulent surface fluxes over heterogeneous sea ice surfaces

Turbulent surface fluxes of momentum and sensible and latent heat as well as surface temperature, air temperature, air humidity, and wind speed were measured by the German Falcon research aircraft over the marginal ice zone (MIZ) of the northern Baltic Sea and the Fram Strait. Applying the bulk formulas and the stability functions to the measurements, the roughness lengths for momentum z0, sensible heat zT, and latent heat zq were calculated. As mean values over a wide range of sea ice conditions, we obtain z0 = 5 � 10�4 m, zT = 1 � 10�8 m, and zq = 1 � 10�7 m. These correspond to the following mean values (± standard deviations) of neutral transfer coefficients reduced to 10 m height, CDN10 = (1.9 ± 0.8) � 10�3, CHN10 = (0.9 ± 0.3) � 10�3, and CEN10 = (1.0 ± 0.2) � 10�3. An average ratio of z0/zT � 104 was observed over the range of 10�6 m < z0 < 10�2 m and differs from previously published results over compact sea ice (10�1 < z0/zT < 103). Other observational results over heterogeneous sea ice do not exist. However, our z0/zT ratio approximately agrees with observations over heterogeneous land surfaces. Flux parameterizations based on commonly used roughness lengths ratios (z0 = zT = zq) overestimate the surface heat fluxes compared to our measurements by more than 100%.

Temporal and spatial variability of surface fluxes over the ice edge zone in the northern Baltic Sea

Three land-fast ice stations (one of them was the Finnish research ice breaker Aranda) and the German research aircraft Falcon were applied to measure the turbulent and radiation fluxes over the ice edge zone in the northern Baltic Sea during the Baltic Air-Sea-Ice Study (BASIS) field experiment from 16 February to 6 March 1998. The temporal and spatial variability of the surface fluxes is discussed. Synoptic weather systems passed the experimental area in a rapid sequence and dominated the conditions (wind speed, airsurface temperature difference, cloud field) for the variability of the turbulent and radiation fluxes. At the ice stations, the largest upward sensible heat fluxes of about 100 Wm�2 were measured during the passage of a cold front when the air cooled faster (�5 K per hour) than the surface. The largest downward flux of about �200 Wm�2 occurred during warm air advection when the air temperature reached +10�C but the surface temperature remained at 0�C. Spatial variability of fluxes was observed from the small scale (scale of ice floes and open water spots) to the mesoscale (width of the ice edge zone). The degree of spatial variability depends on the synoptic situation: during melting conditions downward heat fluxes were the same over ice and open water, whereas during strong cold-air advection upward heat fluxes differed by more than 100 Wm�2. A remarkable amount of grey ice with intermediate surface temperature was observed. The ice in the Baltic Sea cannot be described by one ice type only.

The BALTIMOS (BALTEX Integrated Model System) fi eld experiments: A comprehensive atmospheric boundary layer data set for model validation over the open and ice- covered Baltic Sea

A comprehensive atmospheric boundary layer (ABL) data set was collected in eight fi eld experiments (two during each season) over open water and sea ice in the Baltic Sea during 1998–2001 with the primary objective to validate the coupled atmospheric- ice-ocean-land surface model BALTIMOS (BALTEX Integrated Model System). Measurements were taken by aircraft, ships and surface stations and cover the mean and turbulent structure of the ABL including turbulent fl uxes, radiation fl uxes, and cloud conditions. Measurement examples of the spatial variability of the ABL over the ice edge zone and of the stable ABL over open water demonstrate the wide range of ABL conditions collected and the strength of the data set which can also be used to validate other regional models.

Processes controlling surface, bottom and lateral melt of Arctic sea ice in a state of the art sea ice model

We present a modelling study of processes controlling the summer melt of the Arctic sea ice cover. We perform a sensitivity study and focus our interest on the thermodynamics at the ice–atmosphere and ice–ocean interfaces. We use the Los Alamos community sea ice model CICE, and additionally implement and test three new parametrization schemes: (i) a prognostic mixed layer; (ii) a three equation boundary condition for the salt and heat flux at the ice–ocean interface; and (iii) a new lateral melt parametrization. Recent additions to the CICE model are also tested, including explicit melt ponds, a form drag parametrization and a halodynamic brine drainage scheme. The various sea ice parametrizations tested in this sensitivity study introduce a wide spread in the simulated sea ice characteristics. For each simulation, the total melt is decomposed into its surface, bottom and lateral melt components to assess the processes driving melt and how this varies regionally and temporally. Because this study quantifies the relative importance of several processes in driving the summer melt of sea ice, this work can serve as a guide for future research priorities.

Improvement and sensitivity analysis of thermal thin-ice thickness retrievals

Considering the sea ice decline in the Arctic during the last decades, polynyas are of high research interest since these features are core areas of new ice formation. The determination of ice formation requires accurate retrieval of polynya area and thin-ice thickness (TIT) distribution within the polynya.We use an established energy balance model to derive TITs with MODIS ice surface temperatures (Ts) and NCEP/DOE Reanalysis II in the Laptev Sea for two winter seasons. Improvements of the algorithm mainly concern the implementation of an iterative approach to calculate the atmospheric flux components taking the atmospheric stratification into account. Furthermore, a sensitivity study is performed to analyze the errors of the ice thickness. The results are the following: 1) 2-m air temperatures (Ta) and Ts have the highest impact on the retrieved ice thickness; 2) an overestimation of Ta yields smaller ice thickness errors as an underestimation of Ta; 3) NCEP Ta shows often a warm bias; and 4) the mean absolute error for ice thicknesses up to 20 cm is ±4.7 cm. Based on these results, we conclude that, despite the shortcomings of the NCEP data (coarse spatial resolution and no polynyas), this data set is appropriate in combination with MODIS Ts for the retrieval of TITs up to 20 cm in the Laptev Sea region. The TIT algorithm can be applied to other polynya regions and to past and future time periods. Our TIT product is a valuable data set for verification of other model and remote sensing ice thickness data.