Articles | Volume 22, issue 7
Atmos. Chem. Phys., 22, 4763–4778, 2022
Atmos. Chem. Phys., 22, 4763–4778, 2022
Research article
12 Apr 2022
Research article | 12 Apr 2022

Ship-based estimates of momentum transfer coefficient over sea ice and recommendations for its parameterization

Ship-based estimates of momentum transfer coefficient over sea ice and recommendations for its parameterization
Piyush Srivastava1,a, Ian M. Brooks1, John Prytherch2, Dominic J. Salisbury1, Andrew D. Elvidge3, Ian A. Renfrew3, and Margaret J. Yelland4 Piyush Srivastava et al.
  • 1School of Earth & Environment, University of Leeds, Leeds, LS2 9JT, UK
  • 2Department of Meteorology, Stockholm University, Stockholm, Sweden
  • 3School of Environmental Sciences, University of East Anglia, Norwich, NR4 7TJ, UK
  • 4National Oceanography Centre, Southampton, SO14 3ZH, UK
  • anow at: Centre of Excellence in Disaster and Mitigation and Management, Indian Institute of Technology, Roorkee, India

Correspondence: Piyush Srivastava (


A major source of uncertainty in both climate projections and seasonal forecasting of sea ice is inadequate representation of surface–atmosphere exchange processes. The observations needed to improve understanding and reduce uncertainty in surface exchange parameterizations are challenging to make and rare. Here we present a large dataset of ship-based measurements of surface momentum exchange (surface drag) in the vicinity of sea ice from the Arctic Clouds in Summer Experiment (ACSE) in July–October 2014, and the Arctic Ocean 2016 experiment (AO2016) in August–September 2016. The combined dataset provides an extensive record of momentum flux over a wide range of surface conditions spanning the late summer melt and early autumn freeze-up periods, and a wide range of atmospheric stabilities. Surface exchange coefficients are estimated from in situ eddy covariance measurements. The local sea-ice fraction is determined via automated processing of imagery from ship-mounted cameras. The surface drag coefficient, CD10n, peaks at local ice fractions of 0.6–0.8, consistent with both recent aircraft-based observations and theory. Two state-of-the-art parameterizations have been tuned to our observations, with both providing excellent fits to the measurements.

1 Introduction

The Arctic region is changing rapidly. Surface temperatures are rising at a rate more than twice the planetary average, a process known as Arctic Amplification (Serreze and Barry, 2011; Cohen et al., 2014; Stuecker et al., 2018; Dai et al., 2019). Such rapid warming is drastically altering the physical landscape of the Arctic, most visibly the dramatic reduction in sea-ice extent (Onarheim et al., 2018), thickness and age (Ricker et al., 2017; Kwok, 2018), and has the potential to impact a host of biological and chemical processes (Howes et al., 2015; Lehnherr et al., 2018). Changes in the Arctic may also impact lower latitudes via modification of weather patterns and ocean circulation (Cohen et al., 2014; Overland et al., 2016).

Although climate models robustly reproduce Arctic amplification, they have been less successful in making accurate seasonal forecasts of sea-ice extent (Stroeve et al., 2014) or even capturing the observed sea-ice decline over the past decades (Stroeve et al., 2012). There is also large inter-model variability in projections of future climate over varying timescales (Hodson et al., 2013; Stroeve et al., 2014; Zampieri et al., 2018). A major source of uncertainty in models is the representation of turbulence-driven surface exchanges (Bourassa et al., 2013; Vihma et al., 2014; Tsamados et al., 2014; LeMone et al., 2018). Turbulent exchange is a subgrid-scale process parameterized in terms of resolved model variables and surface transfer coefficients. A lack of observational data in high-latitude environments has resulted in great uncertainty in the parameterization of the transfer coefficients of momentum (CD), heat (CH) and moisture (CE). Here, we focus on the parameterization of the momentum transfer (drag) coefficient, CD.

The exchange of momentum between the atmosphere and sea ice directly affects the dynamical evolution of both the atmospheric boundary layer and the sea ice. The exchange is partly dependent on physical properties of the surface. With ongoing sea ice loss and the increasing spatial extent of the Arctic Ocean's marginal ice zone (MIZ) (Strong and Rigor, 2013; Rolph et al., 2020), the nature of this exchange is subject to change, implying that improved understanding of the physical processes is critical. Recent studies have shown that the model reproduction of future sea-ice thickness and extent (Rae et al., 2014; Tsamados et al., 2014), the near-surface atmosphere (Rae et al., 2014; Renfrew et al., 2019) and the polar ocean (Stössel et al., 2008; Roy et al., 2015) are all sensitive to the parameterization of surface momentum exchange over sea ice.

Most models have rather simplified approaches for parameterizing the transfer coefficients over sea ice: prescribing either a constant value for equivalent neutral transfer coefficients for all sea ice, or two different values, corresponding to the MIZ and pack ice conditions along with empirical ice morphological parameters (Notz, 2012; Lüpkes et al., 2013; Elvidge et al., 2016). They then typically utilize a classical “mosaic” or “flux-averaging” approach, where fluxes are estimated separately over sea ice and open water for each grid box and an “effective” turbulent flux is calculated as the weighted average using the fractions of open water and sea ice (Claussen, 1990; Vihma, 1995).

For models that assume a fixed CD10n over ice, the flux-averaging method leads to a monotonically increasing CD10n across the MIZ; this is not supported by observations (Hartman et al., 1994; Mai et al., 1996; Schröder et al., 2003; Andreas et al., 2010; Elvidge et al., 2016), which indicate a peak at ice fractions of 50 %–80 %. The value of sea-ice concentration at which CD10n peaks depends upon the ice morphology (Elvidge et al., 2016). It arises because of the contribution of form drag at the edges of floes, leads, melt ponds and ridges (Arya, 1973, 1975; Andreas et al., 2010; Lüpkes et al., 2012, Lüpkes and Gryanik, 2015; Elvidge et al., 2016).

Andreas et al. (2010) suggested a simple empirically based parameterization of CD10n in terms of a quadratic function of ice concentration. Based on theoretical considerations (Arya, 1973, 1975; Hanssen-Bauer and Gjessing, 1988; Garbrecht et al., 2002; Birnbaum and Lüpkes, 2002; Lüpkes and Birnbaum, 2005), Lüpkes et al. (2012; L2012 hereafter) developed a physically based hierarchical parameterization for CD10n, which, at its lowest level of complexity, requires only ice fraction as the independent variable. The L2012 parameterization scheme qualitatively reproduces the observed peak in CD10n over the MIZ.

Recently, Elvidge et al. (2016; hereafter E2016) used aircraft measurements over the Arctic MIZ to develop a dataset of 195 independent estimates of CD10n over the MIZ, more than doubling the number of observations previously available. Their observations were consistent with the theory of L2012; however, they found a large variation in CD10ni (CD10n for 100 % ice cover) demonstrating that this depends strongly on ice morphology – as also found by Castellani et al. (2014), who applied bulk parameterizations to ice morphology data based on laser altimetry. E2016 recommended modified values of key parameters in the L2012 scheme and, subsequently, this scheme with these settings has been implemented in the Met Office Unified Model (MetUM). Renfrew et al. (2019) demonstrated that this new scheme significantly reduced biases and root-mean-square errors in the simulated wind speed, air temperature and momentum flux over, and downstream of, the MIZ; in addition to having widespread impacts throughout the Arctic and Antarctic via, for example, mean sea level pressure. The new scheme became part of the operational forecasting system at the Met Office in September 2018 and part of the latest climate model configuration (in GL8). However, at present a constant value of CD10ni is used, a known limitation in the veracity of surface momentum exchange over sea ice.

At present, the complexity of physically based parameterizations of momentum exchange over sea ice exceeds the parameterization constraints provided by observations. In other words, despite recent progress, we are still lacking the observational datasets required for further parameterization development. Here, we utilize a large dataset of ship-based measurements of surface momentum exchange – made as part of the Arctic Clouds in Summer Experiment (ACSE) in July–October 2014, and the Arctic Ocean 2016 expedition (AO2016) in August–September 2016 – to study momentum exchange over heterogeneous sea ice. We investigate the relationship between surface drag and sea-ice concentration within the existing framework suggested by L2012, E2016, and its recent extension by Lüpkes and Gryanik (2015; hereafter L2015) using over 500 new estimates of surface drag and local sea-ice concentration measurements derived from on-board imagery, over varying sea-ice conditions and a range of near-surface atmospheric stabilities.

2 Parameterization background

The surface flux of momentum is τ=-ρu*2=ρCDU2, where ρ is the air density, u* is the friction velocity, and U is the wind speed at a reference height. The drag coefficient, CD, is derived from Monin–Obukhov similarity theory (MOST, Monin and Obukhov, 1954) as

(1) C D = κ 2 ln z / z 0 - ψ m z / L + ψ m z 0 / L - 2

Here, κ is the von Kármán constant, z is the reference height at which the transfer coefficient is evaluated, z0 is the aerodynamic roughness length, L is the Obukhov length, and ψm is an integrated stability correction function (Stull, 1988). Here the small term, ψmz0/L is neglected.

Over land surfaces, the aerodynamic roughness length is in general taken as constant depending upon the surface characteristics, whereas over the open water the roughness length varies with wind speed and is typically parameterized using a type of Charnock relation (Charnock, 1955). When the surface consists of a mix of ice and open water, an effective turbulent flux over the area is usually calculated by taking a weighted average over the fraction of open water and sea ice (Vihma, 1995):

(2) C D 10 n = 1 - A C D 10 nw + A C D 10 ni .

Here, CD10nw and CD10ni are the neutral transfer coefficients for momentum over water and ice surfaces respectively, and A is the fraction of the surface covered by ice. Over sea ice, an additional drag contribution, the form drag, CD10nf, is generated owing to air-flow pressure against the edges of floes, leads and melt ponds (Andreas et al., 2010; L2012; L2015; E2016). The overall equivalent neutral drag coefficient is then given by

(3) C D 10 n = 1 - A C D 10 nw + A C D 10 ni + C D 10 nf .

Lüpkes et al. (2012) proposed a hierarchical parameterization for CD10n in which form drag, at its lowest level of complexity, is parameterized as a function of ice fraction only:

(4) C D 10 nf = A h f D i S c 2 c e 2 ln 2 h f / z 0 w ln 2 10 / z 0 w .

Here, Di is the characteristic length scale of the floe, hf is the freeboard height, Sc is the sheltering function, and ce is the effective resistance coefficient. Lüpkes et al. (2012) provided simplified forms for these parameters, either in terms of ice fraction or as constants:

(5) S c = 1 - exp - s D w h f ,




(9) A * = 1 1 - D min / D max 1 / β .

For operational purposes, L2012 suggested optimal values of the parameters used in the above expressions (Table 1). E2016 evaluated the L2012 scheme with in situ aircraft measurements and found that slightly modifying values of these key parameters (Table 1) represented the behaviour of CD10n well.

The L2012 scheme assumes that the wind profile is always adjusted to the local surface. However, this assumption is not necessarily valid where the surface conditions change over small spatial scales, and the fetch over the local surface is insufficient for the wind profile to come into equilibrium with its characteristics. To overcome this issue in the existing schemes, L2015 suggested a fetch-dependent parameterization of the form drag component of the total drag at arbitrary height:

(10) C Dnf = C Dnf , w 1 - A + C Dnf , i A ,

where CDnf,w and CDnf,i are, respectively, the neutral form-drag coefficients related to the fetch over open water and over ice, and are expressed as

(11) C Dnf , k = A h f D i S c , k 2 c e 2 ln 2 h f / e z 0 , k ln 2 z p / z 0 , k , where k = i , w .

Thus, the L2015 scheme incorporates two form drag contributions and both are weighted by their respective surface fractions. Equation (11) differs from the formulation in L2012 (Eq. 4) only by the inclusion of the Eulerian number, e, in the logarithmic term of the numerator, consistent with previous work by other groups (e.g. Hanssen-Bauer and Gjessing, 1988), and is valid for any reference height, zp. Here we evaluate L2012 and L2015 against in situ estimates of CD10n to assess the impact of including form drag for both water and ice surfaces. We do not evaluate the higher levels of complexity in L2015. The values of various parameters used in the L2015 parameterization, both as originally published and tuned to our observations, are given in Table 1, along with those for L2012 and E2016.

Table 1Parameter settings for the form drag component of the L2012 scheme (Lüpkes et al., 2012, rows 1–4): as recommended in L2012, E2016A and E2016B, P2021–L2012; and the L2015 scheme (Lüpkes and Gryanik, 2015, rows 5–8): as recommended in L2015, and fit to ACSE + AO2016, ACSE only, AO2016 only. The L2012 variants use: Dmin= 8 m and Dmax= 300 m, whereas the L2015 variants use Dmin= 300 m. The primary tuning parameter is the effective resistance coefficient, ce, whereas β has a second-order effect on the shape of the curve. Sc and hf were tunable parameters in L2012, but found by L2015 to have a marginal impact and set as constants for simplicity. P2021–L2015 is the proposed parameterization which is presented with L2015 anchored at the observed values of CD10nw and CD10ni.

Download Print Version | Download XLSX

3 Measurement and methods

3.1 Field measurements

We utilize data from two field campaigns, the Arctic Cloud in Summer Experiment (ACSE, Tjernström et al., 2015, 2019; Achtert et al., 2020), part of the Swedish–Russian–US Arctic Ocean Investigation on Climate-Cryosphere-Carbon (SWERUS-C3), and the Arctic-Ocean 2016 (AO2016) expedition. Both ACSE and AO2016 were carried out on board the Swedish icebreaker Oden. The ACSE cruise took place between 5 July and 5 October 2014, starting and ending in Tromsø, Norway, and working around the Siberian shelf, through the Kara, Laptev, East Siberian and Chukchi seas (Fig. 1). There was a change of crew and science team in Utqiaġvik (formerly Barrow), Alaska, on 20 August. The AO2016 expedition took place between 8 August and 19 September 2016 in the central Arctic Ocean, starting from, and returning to, Longyearbyen, Svalbard (Fig. 1).

3.2 Surface turbulence and meteorological measurements

Turbulent fluxes were measured with an eddy covariance system installed at the top of Oden's foremast, 20.3 m above the waterline. On ACSE this consisted of a Metek USA-100 sonic anemometer with heated sensing heads, a Li-COR Li-7500 open path gas analyser, and an Xsens MTi-700-G motion sensing package installed at the base of the anemometer. The ship's absolute heading and velocity were obtained from its navigation system. The Metek sonic anemometer failed at the start of AO2016 and was replaced with a Gill R3 sonic anemometer. The raw turbulent wind components, at 20 Hz, were corrected for platform motion following Edson et al. (1998) and Prytherch et al. (2015). Corrections for flow distortion of the mean wind were derived from a computational fluid dynamics model (Moat et al., 2015). Turbulent fluxes of heat, momentum and moisture were estimated by the eddy covariance technique over 30 min averaging intervals.

Figure 1The cruise tracks of ACSE (leg 1, 5 July to 18 August 2014 (magenta) and leg 2, 21 August to 5 October (black)) and AO2016 (8 August to 19 September 2016, (blue)) with (a) sea-ice fraction from in situ imagery and (b) CD10n, for each 30 min flux period shown. The sea-ice extent from AMSR2 on 7 August 2014 (white) and 2016 (cyan) – about midway through ACSE and at the start of AO2016 – are shown for reference and to give an indication of the variability between years.

Mean temperature (T) and relative humidity (RH) at the mast top were measured with an aspirated sensor – a Rotronic T/ RH sensor during ACSE and a Vaisala HMP-110 during AO2016. Additional T, RH and pressure (P) measurements were made by a Vaisala PTU300 sensor on the seventh deck of the ship. Pressure at the mast top was obtained by height-adjusting the measurement from the seventh deck. The surface skin temperature was obtained from two Heitronics KT15 infra-red surface temperature sensors, mounted above the bridge and viewing the surface on either side of the ship. Digital imagery of the surface around the ship was obtained from 2 Mobotix M24 IP cameras mounted on the port and starboard beam rails above the bridge, approximately 25 m above the surface. Images were recorded at 1 min intervals during ACSE and at 15 s intervals during AO2016. Profiles of atmospheric thermodynamic structure and winds were obtained from Vaisala RS92 radiosondes, launched every 6 h throughout both cruises. During ACSE a Radiometer Physics HATPRO scanning microwave radiometer provided additional retrievals of lower-atmosphere temperature profiles every 5 min.

3.3 Estimation of turbulence parameters and data screening

The transfer coefficient of momentum is computed as

(12) C D 10 n = u * U 10 n 2 ,

where u* is the measured friction velocity, and U10n is the 10 m equivalent neutral wind speed corresponding to the 10 m wind speed U10, and determined using Monin–Obukhov similarity theory and the Businger–Dyer stability correction function fm (Businger et al., 1971) as U10n=U10fm.

A total of 3421 and 1555 individual half-hourly flux estimates were obtained during ACSE and AO2016 respectively. Data were removed from the analysis if they failed a set of flux quality control criteria (Foken and Wichura, 1996) resulting in a subset of 1804 (247) flux estimates. Additional quality control criteria were applied to filter out data unreliable for analysis of transfer coefficients:

  • The relative wind direction was restricted to ±120 from bow-on, where the flow is clear of the ship's superstructure.

  • Data points where the stability parameter, z/L, was greater than 1 or less than −2 were removed to avoid the effects of strong stability and instability.

  • Sign constancy between the turbulent heat flux and mean gradients was enforced. The Richardson number, RiB, and z/L should always have the same sign, whereas the sensible heat flux, H, should have the opposite sign to z/L. Inconsistencies may arise owing to combined measurement uncertainties where the temperature gradient or H is small.

  • Data were also removed where the 10 m wind speed was less than 3 m s−1.

After initial quality control we have a total of 1403 and 162 half-hourly flux estimates. For ACSE data, an additional quality control criterion was applied based on the boundary-layer profiles. Working around the MIZ, ACSE experienced multiple warm air advection events. These result in strong near-surface air mass modification and the formation of very low or surface-based temperature inversions (Tjernström et al., 2019), which can exhibit significant spatial and temporal variability. Temperature profiles from the HATPRO microwave radiometer, bias-corrected through extensive comparisons with 6-hourly radiosondes, were used to detect surface inversions with 5 min temporal resolution following Tjernström et al. (2019). The profiles were classified as surface-based inversions, low-level inversions (inversion base height < 200 m), and well-mixed boundary layers (inversion base height > 200 m). The flux periods with mixed surface-based and low-level inversions were discarded from the analysis on the basis that the change in near-surface thermodynamic structure was likely to compromise the quality of the flux-profile relationship upon which the calculation of CD10n depends. Following this step, we were left with 1051 data points from the ACSE campaign. No high-frequency profile measurements were available from AO2016; however, operating much further from the ice edge, AO2016 was not subject to the frequent warm air advection and air mass modification events seen during ACSE.

3.4 Determination of sea-ice concentration

Estimates of sea-ice fraction are drawn from two sources: (i) a local estimate of ice fraction determined from digital imagery from the ship; and (ii) daily ice fractions, derived from satellite-based Advanced Microwave Scanning Radiometer (AMSR2) passive microwave measurements (Spreen et al., 2008).

Our local raw imagery consists of high-definition (2048×1536) images of the surface to port and starboard, obtained from Mobotix MX-M24M IP cameras mounted above the ship's bridge, 25 m above the surface. Additional images from a camera pointed over the bow are used for visual inspection while selecting the periods when Oden was in the ice, but are not processed because the ship dominates the near field of the image.

On-board imagery can provide local surface properties including sea ice and melt pond fractions with a spatial resolution of order metres on a time base matched to the flux-averaging time (Weissling et al., 2009). The large volume of imagery sampled requires automated image-processing techniques to estimate ice properties (Perovich et al., 2002; Renner et al., 2014; Miao et al., 2015; Webster et al., 2015; Wright and Polashenski, 2018). Here we use the Open Source Sea-ice Processing (OSSP) algorithm of Wright and Polashenski (2018). The surface properties obtained during each 30 min interval are averaged to give the sea-ice and melt-pond fractions. The image processing methodology is described in Appendix A. Limitations on quality and availability of imagery resulted in a further 206 (36) flux estimates for the ACSE (AO2016) datasets being discarded. After all quality control criterion are applied and flux estimates matched with robust estimates of the local ice fraction, we retain a total of 542 flux estimates: 416 from ACSE and 126 from AO2016. Initially, melt ponds are treated as open water; the impact of this is examined in Appendix B.

The area of each image is approximately 34 225 m2; the total area imaged within each 30 min averaging interval varies with the number of images passing quality control, and the ship's movement, but is up to a maximum of approximately 2 km2 during ACSE and 6.7 km2 during AO2016.

Satellite-based sea-ice products are widely used to prescribe ice concentration in operational forecast models and have been used to assess the dependence of in situ flux measurements as a function of ice fraction (e.g. Prytherch et al., 2017). However, they have significant uncertainties when related to in situ flux measurements due to their relatively coarse temporal and spatial resolution (Weissling et al., 2009), resulting in a mismatch between the satellite footprint and that of the surface flux measurement, and the times of the measurements. The AMSR2 satellite measurements used here provide daily sea-ice concentration on a 6.25 km grid, whereas eddy covariance flux estimates are for 30 min periods and have a footprint of the order of a few hundred metres to a kilometre. The AMSR2 estimates are interpolated spatially to the locations of each flux measurement.

4 Results

4.1 Atmospheric conditions during ACSE and AO2016

Figures 2 and 3 show the meteorological and surface conditions during ACSE and AO2016. The first half of ACSE was dominated by relatively low winds, and surface temperatures close to 0 C when in the ice; much warmer temperatures are associated with open coastal waters. The second half of the cruise experienced higher, and more variable winds, associated with multiple low-pressure systems. Temperatures first fell to the freezing point of salt water on day of year (DoY) 218, although Sotiropoulou et al. (2016) identified the start of freeze up as DoY 241. AO2016 saw a shift from relatively low surface air pressure, and mostly low winds to higher pressure and more variable winds with frequent occurrence of high winds around DoY 237.

Out of the total of 542 flux estimates, we have 184, 282 and 76 flux estimates in stable (z/L> 0.01), near-neutral (−0.01<z/L< 0.01) and unstable conditions (z/L<0.01) respectively. This distribution in static stability augments the limited datasets already available over the marginal ice zone, which have been predominantly in unstable conditions (e.g. E2016).

4.2 Ice surface characteristics

Figure 4 shows the variation over time of the ice, melt-pond, and open-water fractions determined from the on-board imagery and from AMSR2. Note that the variations represent geographic variability along the cruise tracks, as well as temporal changes in ice conditions. During the early phase of the ACSE campaign, to mid-July (DoY 185 to 196), the ice encountered was mostly old ice (Tjernström et al., 2019), with an average ice concentration of about 70 %; melt ponds and open water had 11 % and 18 % coverage respectively. In this phase the average concentration from AMSR2 was 95 % – larger than the sum of average local ice and melt-pond concentrations (81 %). From late July to early August (DoY 209 to 229), the average local ice concentration was 56 %, with melt-pond and open-water fractions of about 17 % and 27 %. Here the AMSR2 concentration was lower than that from imagery, at 44 %.

Figure 2Time series of (a) 10 m neutral wind speed, U10n and surface pressure, Psurf, (secondary axis), (b) friction velocity, u* (c) surface temperature, Tsurf and relative humidity, RH, (secondary axis), (d) air temperature, Tair, (e) 10 m equivalent neutral drag coefficient, CD10n, (f) ice fraction from AMSR2 satellite, (g) Monin–Obukhov stability parameter, z/L, for z/L> 0 (stable) and (h) z/L< 0 (unstable), for ACSE data. The grey dots are 30 min flux periods from the whole cruise, whereas the black dots correspond to the flux data points that pass quality control. In panels (c) and (d), the dashed red lines show Tair=Tsurf= 0 C.


Figure 3As Fig. 2 but for AO2016.


During this period, warm continental air from Siberia flowed northward across the Oden's track causing a rapid melting of ice (Tjernström et al., 2015, 2019). From late August to mid-September (DoY 239 to 255), the average local concentration declined to 27.6 % and the surface was mostly characterized by large areas of open water (63 %). The AMSR2 ice fraction was 54 %, approximately twice that from the imagery. During late September (DoY 268 to 272), a large variation in the surface conditions was observed and often the ice concentration was higher than 90 % owing to the presence of newly formed thin ice, nilas and pancake ice. During the AO2016 campaign, the surface was mostly characterized by old and thick ice, with intermittent patches of thin ice and melt ponds, reflecting the more northerly cruise location. The average ice concentrations from imagery and AMSR2 were found to be about 80 % and 90 % respectively.

Figure 4Time series of ice, melt-pond and open-water fractions (white, blue and black symbols respectively) from the local imagery, and ice fraction (green) from AMSR2, interpolated to the ship location. (a) is for ACSE and (b) for AO2016.


Figure 5 shows a direct comparison of the ice fraction from the in situ imagery and AMSR2. There is a broad correspondence, but a very high degree of scatter, and AMSR2 tends to overestimate the local sea-ice fraction; the correlation coefficient, mean absolute bias and root-mean-square error are 0.64, 0.21 and 0.28 respectively. It is clear from the ice concentration time series, however, that the bias between AMSR2 and the local ice fraction varies over time and appears to be related to the surface conditions of melt or freeze up, in particular when changes are rapid. The largest difference between ice fractions from both projects was found during the early freeze-up season where there is extensive very thin ice.

Figure 5A comparison of the ice fraction derived from the local imagery and from AMSR2 for both field campaigns. The linear regression (AAMSR= 0.584AImg+ 0.321) and 1:1 lines are shown in blue and black respectively.


4.3 Variation of momentum transfer coefficient with sea-ice concentration

We first assess the variability of surface drag with the sea-ice fraction using local ice concentration from the onboard imagery (Fig. 6). The range and median values of CD10nover sea ice (AImg> 0) are similar to those of previous studies (Banke and Smith, 1971; Overland et al., 1985; Guest and Davidson, 1987; Castellani et al., 2014; E2016). The peak CD10n is found at the 0.6–0.8 ice-fraction bin, consistent with L2012 and E2016). The median values of CD10n in both datasets agree well for high ice fractions (Fig. 6b and c); however, there are insufficient AO2016 data for AImg< 0.5 to make a robust comparison with ACSE. Given the good general agreement between ACSE and AO2016, we will only consider the joint dataset from here on.

Figure 6CD10n as a function of ice fraction, as derived from local imagery (AImg) for (a) the joint ACSE and AO2016 datasets (n= 542), (b) the ACSE dataset (n= 416) and (c) the AO2016 dataset (n= 126). The boxes show the interquartile range and the bin median (black squares) for bins of width = 0.2 and plotted at the mean ice fraction for the bin; the number of data points in each bin is noted at the median level. Whiskers indicate the range of the estimates, excluding any outliers, which are plotted individually if present. Parameterization schemes are overlain as indicated, with each curve anchored at the observed median values of CD10nw (AImg= 0) and CD10ni (defined here as AImg> 0.8) for each dataset.


The measurements are compared with the L2012, L2015 and E2016 parameterization schemes. Note that these all require specified values of CD10n over open water (AImg= 0) and solid ice (AImg= 1) (see Eq. 2); these vary with conditions, dramatically so for AImg= 1, as demonstrated by E2016. Here, we follow E2016 and fix the values of CD10nw and CD10ni used in the parameterizations to the measurements, using the observed median values at AImg= 0 and AImg> 0.8 respectively for each dataset. Note that A> 0.8 is used, as opposed to A= 1, as navigational consideration meant that the ship rarely operated in regions with an ice fraction of 1. The mean ice fraction in this bin is 0.89. Here, we do not adjust any of the other tuneable parameters in these parameterizations.

L2012 overestimates the observations of all but the lowest ice concentrations. E2016a and E2016b – which follow L2012 with settings tuned to measurements over the MIZ from Fram Strait and the Barents Sea – correspond well with the observations, with only a slight overestimation of the peak values. L2015, which accounts for form drag over water as well as over ice, is a close match to the observations for AImg> 0.6 but overestimates the 0.2–0.4 and 0.4–0.6 bins and peaks at too low an ice concentration. Note that we will tune the L2015 scheme using our measurements in Sect. 4.4.

The median value of CD10n at A= 0 was 1.65×10-3, which is higher than those typically found over the open ocean (Smith, 1980; Large and Yeager, 2009). This may be a result of the open water measurements being made under fetch-limited conditions close to the ice edge, or within regions of open water within the pack ice, where an under-developed wave state may result in higher drag (Drennan et al., 2003). We cannot, however, exclude the possibility that they result from an incomplete correction for flow distortion over the ship (Yelland et al., 1998, 2002), or that the flux footprint includes flow over nearby ice that is not visible in the imagery.

Figure 7 shows CD10n as a function of the AMSR2 ice fraction. There is broad agreement with the values in Fig. 6a at low and high ice concentrations, but there is no peak in CD10n at intermediate concentrations. Instead, the measurements with higher drags have moved to either lower or higher ice fraction bins. This is consistent with neither our in situ imagery nor previous aircraft-based studies, suggesting that it might be a limitation of the AMSR2 imagery.

Figure 7As Fig. 6a but for the ice fraction derived from the AMSR2 satellite.


4.4 Updating parameterizations using local sea-ice concentration measurements

L2015 extended the parameterization of L2012 to explicitly represent the impact of fetch dependence over heterogeneous surfaces in a physically consistent manner. To date, this scheme has been unconstrained by observational data. Here, we validate the scheme and provide recommendations for its tuneable parameters based on the joint ACSE and AO2016 datasets. E2016 pointed out that variation in the morphological parameters β and ce in L2012 could explain the variability of CD10n within concentration bins. Reducing the values of β and ce from those suggested by L2012 resulted in a better fit to their data.

In the fetch-dependent L2015 parameterization, increasing β (sea-ice morphology exponent, Eqs. 8 and 9) results in decreasing CD10n, mostly at high ice concentrations, whereas increasing ce (the effective resistance coefficient) increases CD10n at all concentrations. Here we have adjusted the L2015 values of β and ce to optimize the fit to our measurements. The revised values of the coefficients are given in Table 1. For a consistent comparison, similar tuning is applied to L2012.

Figure 8 shows CD10n plotted against AImg along with the tuned L2015 and L2012 schemes, both anchored to the observed values of CD10n at AImg= 0 and AImg> 0.8 for the joint dataset.

Figure 8CD10n as a function of the local ice fraction (AImg) for the joint ACSE and AO2016 (n= 542) dataset. Here, the proposed parameterization (red line, P2021–L2015) is presented with L2015 anchored at the observed values of CD10nw (AImg= 0) and CD10ni(AImg> 0.8) with the coefficients β and ce shown in Table 1. For comparison, the L2012 scheme (green line, P2021–L2012) is also tuned with the curve anchored at the same values of CD10ni and CD10nw with β and ce shown in Table 1.


Both L2012 and L2015 provide an excellent fit to the data, passing close to the median observed values at all ice fractions. The fitted curve for the joint dataset (P2021–L2015) works equally well for the individual datasets (Fig. S1, Supplement).

In the analysis above we have considered CD10n as a function of ice fraction – no distinction is made between melt ponds and open water. However, there are uncertainties in the surface classification, in particular for the determination of melt-pond fraction. Thin ice and shallow melt ponds can appear very similar in colour, and potentially be misclassified by the image-processing algorithm. An assessment of the sensitivity of the fitting of the L2015 to the presence and treatment of melt ponds (see Appendix B) shows that they have little impact.

Melt ponds are explicitly included in the L2012 and L2015 parameterizations in their more complex levels of implementation, where the edges of melt ponds provide a source of form drag. Tsamados et al. (2014) modelled the different contributions to the total drag using L2012 implemented within the CICE sea ice model (Hunke and Lipscomb, 2010). They found that melt ponds made a negligible contribution to the drag except over the oldest, thickest ice just north of the Canadian archipelago, consistent with our observations.

Both L2012 and L2015 can be tuned to provide excellent fits to the observations (Fig. 8). Even without tuning to this dataset, the differences between L2015, E2016A and E2016B are modest and all lie within the interquartile range of the observed CD10n at all ice fractions (Fig. 6a). The largest source of uncertainty in the application of these schemes is the value of the drag coefficient at 100 % ice fraction, CD10ni, which must be prescribed, and is strongly dependent on ice morphology. Table 2 lists values of the neutral drag coefficient for very high ice fractions (0.8–1 from this study, and 0.9–1 or 1 in previous studies) reported in the literature. The best estimates (mean or median values) vary by a factor of more than 4. As discussed in previous studies, CD10ni depends on the sea-ice morphology and so prescribing this as one value is a drastic simplification.

Table 2Overview of neutral drag coefficients based on in situ eddy covariance measurements over “complete” sea-ice cover (CD10ni) from this and previous studies. The values are taken from the literature, and so vary as to whether the mean, median or a range of values is shown. The definition of “complete” sea ice covers a range of ice fractions (0.8–1, 0.9–1.0, and exactly 1) depending on the study. N is the number of data points in this bin where specified. The 2nd column provides the interquartile range or the range from 1 to +1 s, where s is the standard deviation. Note Guest and Davidson (1987) uses the turbulence dissipation method. Overland (1985) compiles values from a variety of previous studies in various locations, as well as new data; thus, the values reproduced here are a compilation by morphology.

Download Print Version | Download XLSX

5 Conclusions

An extensive set of measurements of drag coefficients over sea ice, obtained during two research cruises within the Arctic Ocean, has been utilized to evaluate the dependence of drag on the ice fraction. The final dataset consists of 542 estimates of drag coefficients along with estimates of the local ice fraction obtained from high-resolution imagery of the surface around the ship. The measurements cover a wide geographic area, summer melt and early autumn freeze up, and a range of surface conditions from thick multiyear ice, through melting ice with melt ponds, to newly formed thin and pancake ice, and near-surface stability conditions of −2<z/L< 1, a much wider range than E2016. This wide range of conditions means that the results should be broadly representative of much of the Arctic sea ice region.

The dependence of CD10n on the ice fraction is evaluated in the context of state-of-the-art parameterization schemes (Lüpkes et al., 2012; Lüpkes and Gryanik, 2015). The most recent of these (Lüpkes and Gryanik, 2015) attempts to account for the impact of short fetch over ice/water over spatially highly heterogeneous surfaces. When tuned to the observations, the parameterizations provide an excellent representation of CD10n as a function of ice fraction:

The main conclusions are

  • The data support the existence of a negatively skewed distribution of CD10n with ice concentration, with a peak value for fractions of 0.6–0.8, consistent with the predicted behaviour from Lüpkes et al. (2012) and observations of Elvidge et al. (2016).

  • When tuned to our measurements, both L2012 and L2015 provide an excellent fit to the observed variation of CD10n with the ice fraction. The impact of small-scale surface heterogeneity and the influence of fetch is likely to increase with increasing contrast in the skin temperatures of the ice and water surfaces, and thus play a greater role in the winter.

  • Melt ponds had no significant impact on the drag coefficient over the study area. The optimal fit of the L2015 parameterization to the measurements had little sensitivity to the uncertainty in partitioning of melt ponds to the ice or water fractions when estimating the local ice fraction, and there was little sensitivity to the presence of melt ponds at all for the conditions observed.

  • When evaluated against the AMSR2 retrieval of the ice fraction, the behaviour of CD10n is not consistent with in situ observations, for example, no peak is seen at intermediate ice fractions. This is likely a result of several factors: a mismatch in the spatial scale between the in situ flux footprint (of order 100 s of metres to 1 km) and the satellite footprint (6.5 km); potential spatial offsets in location matching resulting from the low temporal resolution of the satellite data (daily retrievals) combined with drifting of the ice; and the high scatter and varying mean bias between the in situ and satellite estimates of the ice fraction. The mean bias in particular displays temporal and spatial coherence that suggests a dependence upon surface conditions. This finding cautions against the use of comparatively low-resolution remote-sensing products when evaluating parameterizations.

Atmospheric stability may also play a role here, as it will affect how rapidly the atmospheric surface layer adjusts to changes in surface properties. L2015 incorporates stability effects at the higher levels of parameterization complexity, but not within the simplest complexity level used here. A much larger dataset, including the details of surface heterogeneity, would be required to evaluate the details of both stability and fetch dependencies.

Sea ice and climate models are starting to incorporate components of form drag within their surface exchange schemes for sea ice (e.g. Tsamados et al., 2014; Renfrew et al., 2019). But at present, most do not use all the components of the more complex versions of schemes such as L2012 or L2015. Instead, they tend to rely on the simplest versions where drag is only a function of the ice fraction. In operational forecast models, where only a prescribed ice concentration from a satellite retrieval may be available, this seems appropriate, but within more complex coupled weather and climate prediction models there is the potential for using output from the sea-ice model to adjust the drag coefficient (E2016; Renfrew et al., 2019). The skill of parameterization is strongly dependent on the accurate representation of the drag at 100 % ice fraction, CD10ni, which varies significantly with ice morphology (Lüpkes et al., 2012; Lüpkes and Gryanik, 2015; Elvidge et al., 2016, 2021). Tackling the representation of CD10ni should be the next challenge in improving air-ice surface drag in weather and climate models.

Appendix A: Image processing and evaluation of the local ice fraction

A total of ∼500 000 images of the surface around the ship were obtained over the two cruises; thus, this required an automated approach to estimating the local ice fraction. Here we use the Open Source Sea-ice Processing (OSSP) algorithm of Wright and Polashenski (2018).

A1 (a) Pre-processing

Of the images available for each flux period, a subset of visibly good images was selected for further processing. The rejection of images was due to the presence of dense fog, moisture or ice on the camera lens, strong surface reflection of direct sunlight, or insufficient illumination. The selected subsets consist of 10 to 60 images in each flux period (e.g. Fig. A1a). These images are first corrected for lens distortion. The lens-specific distortion coefficients and intrinsic parameters were determined using the Computer Vision System Toolbox of MATLAB. The corrected images (2048×1536 pixels) were then cropped to select a region within ∼200 m of the ship (2009×1111 pixels) – e.g. see Fig. A1b.

A2 (b) Training and implementation of the algorithm

The success of any machine learning-based algorithm depends upon the quality of the training dataset. As the ice conditions varied substantially throughout the campaigns, extensive training data were needed to cover the wide range of conditions. The initial training images selected were from the first and last images from each flux period. Additional images were added iteratively depending upon the performance of the algorithm on randomly selected images. After multiple trials, we settled on three different training datasets for (i) images with a visibly large ice fraction, (ii) images with a large open water fraction and (iii) images showing newly formed thin ice. Our approach was to generate a training dataset that could be utilized equally on imagery from other campaigns, while keeping the number of discrete training datasets as small as possible. The training data, identifying ice, water and melt-ponds, were generated based on user classification of the training images via a Graphical User Interface (GUI), and thus depends upon the ability of the user to identify the surface features correctly.

A3 (c) Post-processing

The OSSP algorithm produces an indexed image having pixel-wise information about surface features (open water, melt ponds, ice) for each input image (e.g. Fig. A1c). As the images were necessarily taken at an oblique angle, the indexed images need to be orthorectified to derive the correct fractions of ice, melt pond and water. Orthorectification of imagery is a process by which pixel elements of an oblique image are restored to their true vertical perspective position. The angular separation of each

pixel (after correction for lens distortion) was determined from a laboratory calibration of the cameras. The angle from the horizon (the horizontal) in the images and the height of the camera above the surface then allow the location of each pixel on the surface to be calculated. The masked images were interpolated onto a regular xy grid after orthorectification and area fractions of ice and melt ponds were estimated as a fraction of the total number of pixels for each category (e.g. Fig. A1d). The average fractions of ice and water for a flux period are then calculated by taking an average over all the images in that period. Only flux periods having more than 30 available images are included in the analysis.

Each orthorectified image has an area of approximately 34 225 m2; the total area included in each 30 min average varies with ship manoeuvres and the number of images passing quality control. With a maximum in-ice ship speed of 5 m s−1, the 1 min images from ACSE do not overlap, providing a maximum area of 2.05 km−2. For the 15 s imagery during AO2016, images overlap by 75 m at a ship speed of 5 m s−1, giving a maximum area of approximately 6.7 km−2. It is implicitly assumed that the ice fraction determined along the ship track is representative of that within the flux footprint, i.e. that the ice structure is more or less homogeneous within the footprint.

Figure A1An example of the image processing workflow. Panel (a) is an example raw image; (b) shows the image corrected for the lens distortion, where the region of focus is shown by the yellow rectangle; (c) shows the image after processing by the OSSP algorithm where the masking colours – white, blue and black – represent ice, melt pond or submerged ice and open water areas respectively; (d) shows the orthorectified image showing the true distance of each surface feature away from the camera.


Appendix B: Sensitivity to melt ponds

Here we investigate the sensitivity of our tuning of L2015 to the melt pond fraction. We reclassify 50 % and 100 % of melt ponds as ice instead of water (Fig. B1) and the L2015 function is re-fitted to the revised ice fractions and compared with our original fit.

The reclassification of melt ponds as ice has the effect of moving some drag estimates into higher ice fraction bins, slightly increasing the median value of CD10n at high ice fractions. The refitted L2015 functions reflect this slightly higher drag at high ice fraction but are essentially unchanged for A< 0.5. Note that even when A> 0.5 the change in the L2015 functions is very small compared with the variation in CD10n within each ice concentration bin. We further investigate the sensitivity of the parameterization to the presence of melt ponds by a simple sub-setting of the data by the melt-pond fraction. In all cases the melt-pond fraction is < 0.6. Figure B2 shows CD10n with the ice fraction for cases where the melt-pond fraction is <  0.3 (Fig. B2a) and < 0.1 (Fig. B2b). The L2015 function is fitted to these subsets of data and compared with that of the full dataset. The revised fits differ negligibly from that of the full dataset, suggesting that CD10n might not be strongly dependent on the extent of melt ponds. In short, the sensitivity of the parameterization to the treatment of melt ponds is negligible.

Figure B1Parameterization sensitivity to the melt-pond (MP) fraction. Panel (a) re-classifies half of melt ponds to sea ice (A' = A + 0.5  MP), whereas (b) re-classifies all of the melt points to sea ice (A' = A + MP), where A is the sea-ice fraction. The curves show the L2015 parameterization, tuned to the original ice fraction observations (red; P2021–L2015) and tuned to the adjusted ice fraction observations (black).


Figure B2(a) MP < 30 % (b) MP < 10 %. In each panel, the red curve is the fitting curve obtained for the joint ACSE and AO2016 data, the same as shown in Fig. 8 (P2021–L2015), and the black lines are the fitting curve obtained for the data shown in the “respective” panels.


Data availability

The data used here are publicly available from the Centre for Environmental Data Analysis (CEDA) archives as separate files for each cruise. ACSE cruise data are available from CEDA archives (, Brooks et al., 2022a) and the AO2016 data are available from CEDA archive (, Brooks et al., 2022b).


The supplement related to this article is available online at:

Author contributions

All authors contributed to the design of the study. PS analysed the datasets and wrote the manuscript with contributions from all co-authors. IMB, JP, and DJS collected the data during ACSE; IMB and JP collected the data during AO2016. JP processed the flux data. DJS and PS processed the surface imagery.

Competing interests

The contact author has declared that neither they nor their co-authors have any competing interests.


Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


We would like to thank the captains and crews of the icebreaker Oden, along with the technical and logistical support staff of the Swedish Polar Research Secretariat, for their assistance throughout the ACSE and AO2016 cruises. We thank Michael Tjernström, Ola Persson, Matthew Shupe, Barbara Brooks, Joseph Sedlar, and Georgia Sotiropoulou for their contributions to the ACSE measurement campaign.

Financial support

This research has been supported by the UK Natural Environment Research Council (NERC; grant nos. NE/S000453/1 and NE/S000690/1). Margaret J. Yelland was also supported by NERC (grant nos. NE/N018095/1 and NE/V013254/1). The contribution of Ian M. Brooks, John Prytherch, and Dominic J. Salisbury to the ACSE cruise was funded by NERC (grant no. NE/K011820/1). Participation in the AO2016 cruise was supported by the Swedish Polar Research Secretariat. John Prytherch was also supported by the Knut and Alice Wallenberg Foundation (grant no. 2016-0024). Piyush Srivastava was supported by the Department of Science and Technology (DST), Government of India (grant no. DST/INSPIRE/04/2019/003125).

Review statement

This paper was edited by Peter Haynes and reviewed by Christof Lüpkes and one anonymous referee.


Achtert, P., O'Connor, E. J., Brooks, I. M., Sotiropoulou, G., Shupe, M. D., Pospichal, B., Brooks, B. J., and Tjernström, M.: Properties of Arctic liquid and mixed-phase clouds from shipborne Cloudnet observations during ACSE 2014, Atmos. Chem. Phys., 20, 14983–15002,, 2020. 

Arya, S. P. S.: Contribution of form drag on pressure ridges to the air stress on Arctic ice, J. Geophys. Res., 78, 7092–7099,, 1973. 

Arya, S. P. S.: A drag partition theory for determining the large-scale roughness parameter and wind stress on the Arctic pack ice, J. Geophys. Res., 80, 3447–3454,, 1975. 

Andreas, E. L., Horst, T. W., Grachev, A. A., Persson, P. O. G., Fairall, C. W., Guest, P. S., and Jordan, R. E.: Parametrizing turbulent exchange over summer sea ice and the marginal ice zone, Q. J. Roy. Meteor. Soc., 136, 927–943,, 2010. 

Banke, E. G. and Smith, S. D.: Wind stress over ice and over water in the Beaufort Sea, J. Geophys. Res., 76, 7368–7374,, 1971. 

Birnbaum, G. and Lüpkes, C.: A new parameterization of surface drag in the marginal sea ice zone, Tellus A, 54,, 107–123,, 2002. 

Bourassa, M. A., Gille, S. T., Bitz, C., Carlson, D., Cerovecki, I., Clayson, C. A., Cronin, M. F., Drennan, W. M., Fairall, C. W., Hoffman, R. N., Magnusdottir, G., Pinker, R. T., Renfrew, I. A., Serreze, M., Speer, K., Talley, L. D., and Wick, G. A.: High-Latitude Ocean and Sea Ice Surface Fluxes: Challenges for Climate Research, B. Am. Meteorol. Soc., 94, 403–423,, 2013. 

Brooks, I. M., Prytherch, J., and Srivastava, P.: CANDIFLOS: Surface fluxes from ACSE measurement campaign on icebreaker Oden, 2014, NERC EDS Centre for Environmental Data Analysis [data set],, 2022a. 

Brooks, I. M., Prytherch, J., and Srivastava, P.: CANDIFLOS: Surface fluxes from AO2016 measurement campaign on icebreaker Oden, 2016, NERC EDS Centre for Environmental Data Analysis [data set],, 2022b. 

Businger, J. A., Wyngaard, J. C., Izumi, Y., and Bradley, E. F.: Flux Profile Relationships in the Atmospheric Surface Layer, J. Atmos. Sci., 28, 181–189,<0181:FPRITA>2.0.CO;2, 1971. 

Castellani, G., Lüpkes, C., Hendricks, S., and Gerdes, R.: Variability of Arctic sea-ice topography and its impact on the atmospheric surface drag, J. Geophys. Res.-Oceans, 119, 6743–6762,, 2014. 

Charnock, H.: Wind stress over a water surface, Q. J. Roy. Meteor. Soc., 81, 639–640,, 1955. . 

Claussen, M.: Area-averaging of surface fluxes in a neutrally stratified, horizontally inhomogeneous atmospheric boundary layer, Atmos. Environ., 24, 1349–1360,, 1990. 

Cohen, J., Screen, J. A., Furtado, J. C., Barlow, M., Whittleston, D., Coumou, D., Francis, J., Dethloff, K., Entekhabi, D., Overland, J., and Jones, J.: Recent Arctic amplification and extreme mid-latitude weather, Nat. Geosci., 7, 627–637,, 2014. 

Dai, A., Luo, D., Song, M., and Liu, J.: Arctic amplification is caused by sea-ice loss under increasing CO2, Nat. Commun., 10, 121,, 2019. 

Drennan, W. M., Graber, H. C., Hauser, D., and Quentin, C.: On the wave age dependence of wind stress over pure wind seas, J. Geophys. Res., 108, 8062,, 2003. 

Edson, J. B., Hinton, A. A., Prada, K. E., Hare, J. E., and Fairall, C. W.: Direct covariance flux estimates from mobile platforms at sea, J. Atmos. Ocean. Tech., 15, 547–562,<0547:DCFEFM>2.0.CO;2, 1998. 

Elvidge, A. D., Renfrew, I. A., Weiss, A. I., Brooks, I. M., Lachlan-Cope, T. A., and King, J. C.: Observations of surface momentum exchange over the marginal ice zone and recommendations for its parametrisation, Atmos. Chem. Phys., 16, 1545–1563,, 2016. 

Elvidge, A. D., Renfrew, I. A., Brooks, I. M., Srivastava, P., Yelland, M. J., Prytherch, J.: Surface heat and moisture exchange in the marginal ice zone: Observations and a new parameterization scheme for weather and climate models, J. Geophys. Res., 126, e2021JD034827,, 2021. 

Foken, T. and Wichura, B.: Tools for quality assessment of surface based flux measurements 1, Agr. Forest Meteorol., 78, 83–105,, 1996. 

Garbrecht, T., Lüpkes, C., Hartmann, J., and Wolff, M.: Atmospheric drag coefficients over sea ice – validation of a parameterisation concept, Tellus A, 54, 205–219,, 2002. 

Guest, P. S. and Davidson, K. L.: The Effect of Observed Ice Conditions on the Drag Coefficient in the Summer East Greenland Sea Marginal Ice Zone, J. Geophys. Res., 92, 6943–6954,, 1987. 

Hanssen-Bauer, I. and Gjessing, Y. T.: Observations and model calculations of aerodynamic drag on sea ice in the Fram Strait, Tellus A, 40, 151–161,, 1988. 

Hartmann, J., Kottmeier, C., Wamser, C., and Augstein, E.: Aircraft measured atmospheric momentum, heat and radiation fluxes over Arctic sea ice, The Polar Oceans and their Role in Shaping the Global Environment, 85, 443–454,, 1994. 

Hodson, D. L. R., Keeley, S. P. E., West, A., Ridley, J., Hawkins, E., and Hewitt, H. T.: Identifying uncertainties in Arctic climate change projections, Clim. Dynam., 40, 2849–2865,, 2013. 

Howes, E. L., Joos, F., Eakin, M., and Gattuso, J.-P.: An updated synthesis of the observed and projected impacts of climate change on the chemical, physical and biological processes in the oceans, Front. Mar. Sci., 2, 36,, 2015. 

Hunke, E. C., and Lipscomb, W. H.: CICE: The Los Alamos sea ice model documentation and software user's manual version 4.1, Los Alamos National Laboratory Tech. Rep. LA-CC-06-012, 76 pp.,'s_manual.pdf (last access: 20 January 2020), 2010. 

Kwok, R.: Arctic Sea ice thickness, volume, and multiyear ice coverage: losses and coupled variability (1958–2018). Environ. Res. Lett., 13, 105005,, 2018. 

Large, W. G. and Yeager, S. G.: The global climatology of an interannually varying air–sea flux data set, Clim. Dynam., 33, 341–364,, 2009. 

Lehnherr, I., St. Louis, V. L., Sharp, M., Gardner, A., Smol, J. P., Schiff, S. L., Muir, D. C. G., Mortimer, C. A., Michelutti, N., Tarnocai, C., St. Pierre, K. A., Emmerton, C. A., Wiklund, J. A., Köck, G., Lamoureux, S. F., and Talbot, C. H.: The world's largest High Arctic lake responds rapidly to climate warming, Nat. Commun., 9, 1290,, 2018. 

LeMone, M. A., Angevine, W. M., Bretherton, C. S., Chen, F., Dudhia, J., Fedorovich, F., Katsaros, K. B., Lenschow, D. H., Mahrt, L., Patton, E. G., Sun, J., Tjernström, M., and Weil, J.: 100 Years of Progress in Boundary Layer Meteorology, Meteor. Mon., 59, 9.1–9.85,, 2018. 

Lüpkes, C. and Birnbaum, G.: Surface drag in the Arctic marginal sea-ice zone: A comparison of different parameterisation concepts, Bound.-Lay. Meteorol., 117, 179–211,, 2005. 

Lüpkes, C. and Gryanik, V. M.: A stability-dependent parametrization of transfer coefficients for momentum and heat over polar sea ice to be used in climate models, J. Geophys. Res.-Atmos., 120, 552–581,, 2015. 

Lüpkes, C., Gryanik, V. M., Hartmann, J., and Andreas, E. L.: A parametrization, based on sea ice morphology, of the neutral atmospheric drag coefficients for weather prediction and climate models, J. Geophys. Res., 117, D13112,, 2012. 

Lüpkes, C., Gryanik, V. M., Rösel, A., Birnbaum, G., and Kaleschke, L.: Effect of sea ice morphology during Arctic summer on atmospheric drag coefficients used in climate models, Geophys. Res. Lett., 40, 446–451,, 2013. 

Mai, S., Wamser, C., and Kottmeier, C.: Geometric and aerodynamic roughness of sea ice, Bound.-Lay. Meteorol., 77, 233–248,, 1996. 

Miao, X., Xie, H., Ackley, S., Perovich, D., and Ke, C.: Object based detection of Arctic sea ice and melt ponds using high spatial resolution aerial photographs, Cold Reg. Sci. Technol., 119, 211–222,, 2015. 

Moat, B. I., Yelland, M. J., and Brooks, I. M.: Airflow distortion at instrument sites on the ODEN during the ACSE project (National Oceanography Centre Internal Document, 17) Southampton, GB. National Oceanography Centre 114 pp., available at: (last access: 16 January 2020), 2015. 

Monin, A. S. and Obukhov, A. M.: Osnovnye zakonomernosti turbulentnogo peremeshivanija v prizemnom sloe atmosfery [Basic Laws of Turbulent Mixing in the Atmosphere Near the Ground], Trudy geofiz. inst. AN SSSR 24(151), 163–187, 1954. 

Notz, D.: Challenges in simulating sea ice in Earth System Models, Wiley Interdiscip, Rev. Clim. Change, 3, 509–526,, 2012. 

Onarheim, I. H., Eldevik, T., Smedsrud, L. H., and Stroeve, J. C.: Seasonal and regional manifestation of Arctic Sea ice loss, J. Climate, 31, 4917–4932,, 2018. 

Overland, J. E.: Atmospheric boundary-layer structure and drag coefficients over sea ice, J. Geophys. Res., 90, 9029–9049,, 1985. 

Overland, J. E., Dethloff, K., Francis, J. A., Hall, R. J., Hanna, E., Kim, S., Screen, J. A., Shepherd, T. G., and Vihma, T.: Nonlinear response of mid-latitude weather to the changing Arctic, Nat. Clim. Change, 6, 992–999,, 2016. 

Perovich, D. K., Tucker, W. B., and Ligett, K. A.: Aerial observations of the evolution of ice surface conditions during summer, J. Geophys. Res., 107, 8048,, 2002. 

Prytherch, J., Yelland, M. J., Brooks, I. M., Tupman, D. J., Pascal, R. W., Moat, B. I., and Norris, S. J.: Motion-correlated flow distortion and wave-induced biases in air–sea flux measurements from ships, Atmos. Chem. Phys., 15, 10619–10629,, 2015. 

Prytherch, J., Brooks, I. M., Crill, P. M., Thornton, B. F., Salisbury, D. J., Tjernström, M., Anderson, L. G., Geibel, M. C., and Humborg, C.: Direct determination of the air-sea CO2 gas transfer velocity in Arctic sea ice regions, Geophys. Res. Lett., 44, 3770–3778,, 2017. 

Rae, J. G. L., Hewitt, H. T., Keen, A. B., Ridley, J. K., Edwards, J. M., and Harris, C. M.: A sensitivity study of the sea ice simulation in the global coupled climate model, HadGEM3, Ocean Model., 74, 60–76,, 2014. 

Renner, A. H. H., Gerland, S., Haas, C., Spreen, G., Beckers, J. F., Hansen, E., Nicolaus, M., and Goodwin, H.: Evidence of Arctic sea ice thinning from direct observations, Geophys. Res. Lett., 41, 5029–5036,, 2014. 

Renfrew, I. A., Elvidge, A. D., and Edwards, J. M: Atmospheric sensitivity to marginal-ice zone drag: local and global responses, Q. J. Roy. Meteor. Soc., 145, 1165–1179,, 2019. 

Ricker, R., Hendricks, S., Kaleschke, L., Tian-Kunze, X., King, J., and Haas, C.: A weekly Arctic sea-ice thickness data record from merged CryoSat-2 and SMOS satellite data, The Cryosphere, 11, 1607–1623,, 2017. 

Rolph, R. J., Feltham, D. L., and Schröder, D.: Changes of the Arctic marginal ice zone during the satellite era, The Cryosphere, 14, 1971–1984,, 2020. 

Roy, F., Chevallier, M., Smith, G., Dupont, F., Garric, G., Lemieux, J.-F., Lu, Y., and Davidson, F.: Arctic sea ice and freshwater sensitivity to the treatment of the atmosphere-ice ocean surface layer, J. Geophys. Res.-Oceans., 120, 4392–4417,, 2015. 

Schröder, D., Vihma, T., Kerber, A., and Brümmer, B.: On the parameterisation of Turbulent Surface Fluxes Over Heterogeneous Sea Ice Surfaces, J. Geophys. Res., 108, 3195,, 2003. 

Serreze, M. C. and Barry, R. G.: Processes and impacts of Arctic amplification: A research synthesis, Global Planet. Change, 77, 85–96,, 2011. 

Smith, S.: Wind stress and heat flux over the ocean in gale force winds, J. Phys. Oceanogr., 10, 709–726,<0709:WSAHFO>2.0.CO;2, 1980. 

Sotiropoulou, G., Tjernström, M., Sedlar, J., Achtert, P., Brooks, B. J., Brooks, I. M., Persson, P. O. G., Prytherch, J., Salisbury, D. J., Shupe, M. D., and Johnston, P. E.: Atmospheric Conditions during the Arctic Clouds in Summer Experiment (ACSE): Contrasting Open Water and Sea Ice Surfaces during Melt and Freeze-Up Seasons, J. Climate, 29, 8721–8744,, 2016. 

Spreen, G., Kaleschke, L., and Heygster, G.: Sea ice remote sensing using AMSR-E 89-GHz channels, J. Geophys. Res., 113, C02S03,, 2008. 

Stössel, A., Cheon, W.-G., and Vihma, T.: Interactive momentum flux forcing over sea ice in a global ocean GCM, J. Geophys. Res., 113, C05010,, 2008. 

Stroeve, J. C., Kattsov, V., Barrett, A., Serreze, M., Pavlova, T., Holland, M., and Meier, W. N.: Trends in Arctic sea ice extent from CMIP5, CMIP3 and observations, Geophys, Res. Letts., 39, L16502,, 2012. 

Stroeve, J. C., Hamilton, L., Bitz, C., and Blanchard-Wigglesworth, E.: Predicting September sea ice: Ensemble skill of the SEARCH sea ice outlook 2008–2013, Geophys. Res. Lett., 41, 2411–2418,, 2014. 

Strong, C. and Rigor, I. G.: Arctic marginal ice zone trending wider in summer and narrower in winter, Geophys. Res. Lett., 40, 4864–4868,, 2013. 

Stuecker, M. F., Bitz, C. M., Armour, K. C., Proistosescu, C., Kang, S. M., Xie, S. P., Kim, D., McGregor, S., Zhang, W., Zhao, S., Cai, W., Dong, Y., and Jin, F. F.: Polar amplification dominated by local forcing and feedbacks, Nat. Clim. Change, 8, 1076–1081,, 2018. 

Stull, R. B.: An introduction to boundary layer meteorology, Kluwer Academic Publishers, Dordrecht,, 1988. 

Tjernström, M., Shupe, M. D., Brooks, I. M., Persson, P. O. G., Prytherch, J., Salisbury, D. J., Sedlar, J., Achtert, P., Brooks, B. J., Johnston, P. E., and Sotiropoulou, G.: Warm-air advection, air mass transformation and fog causes rapid ice melt, Geophys. Res. Lett., 42, 5594–5602,, 2015. 

Tjernström, M., Shupe, M. D., Prytherch, J., Achtert, P., Brooks, I. M., and Sedlar, J.: Arctic summer air-mass transformation, surface inversions and the surface energy budget, J. Climate, 32, 769–789,, 2019. 

Tsamados, M., Feltham, D. L., Schroeder, D., Flocco, D., Farrell, S. L., Kurtz, N., Laxon, S. L., and Bacon, S.: Impact of Variable Atmospheric and Oceanic Form Drag on Simulations of Arctic Sea Ice, J. Phys. Oceanogr., 44, 1329–1353,, 2014. 

Vihma, T.: Subgrid Parameterization of Surface Heat and Momentum Fluxes over Polar Oceans, J. Geophys. Res., 100, 22625–22646,, 1995. 

Vihma, T., Pirazzini, R., Fer, I., Renfrew, I. A., Sedlar, J., Tjernström, M., Lüpkes, C., Nygård, T., Notz, D., Weiss, J., Marsan, D., Cheng, B., Birnbaum, G., Gerland, S., Chechin, D., and Gascard, J. C.: Advances in understanding and parameterization of small-scale physical processes in the marine Arctic climate system: a review, Atmos. Chem. Phys., 14, 9403–9450,, 2014. 

Webster, M. A., Rigor, I. G., Perovich, D. K., Richter-menge, J. A., Polashenski, C. M., and Light, B.: Seasonal evolution of melt ponds on Arctic sea ice, J. Geophys. Res., 120, 5968–5982,, 2015. 

Weissling, B., Ackley, S., Wagner, P., and Xie, H.: EISCAM – Digital image acquisition and processing for sea ice parameters from ships, Cold Reg. Sci. Technol., 57, 49–60,, 2009. 

Wright, N. C. and Polashenski, C. M.: Open-source algorithm for detecting sea ice surface features in high-resolution optical imagery, The Cryosphere, 12, 1307–1329,, 2018. 

Yelland, M. J., Moat, B. I., Taylor, P. K., Pascal, R. W., Hutchings, J., and Cornell, V. C.: Wind stress measurements from the open ocean corrected for airflow distortion by the ship, J. Phys. Ocean., 28, 1511–1526,<1511:WSMFTO>2.0.CO;2, 1998. 

Yelland, M. J., Moat, B. I., Pascal, R. W., and Berry, D. I.: CFD model estimates of the airflow distortion over research ships and the impact on momentum flux measurements. J. Atmos. Ocean Tech., 19, 1477–1499,<1477:CMEOTA>2.0.CO;2, 2002. 

Zampieri, L., Goessling, H. F., and Jung, T.: Bright prospects for Arctic sea ice prediction on subseasonal time scales, Geophys. Res. Letts., 45, 9731–9738,, 2018. 

Short summary
The parameterization of surface turbulent fluxes over sea ice remains a weak point in weather forecast and climate models. Recent theoretical developments have introduced more extensive physics but these descriptions are poorly constrained due to a lack of observation data. Here we utilize a large dataset of measurements of turbulent fluxes over sea ice to tune the state-of-the-art parameterization of wind stress, and compare it with a previous scheme.
Final-revised paper