Introduction

There is a growing demand for a more accurate prediction and understanding of Arctic weather and climate, which puts great expectations on the global and regional forecasting and reanalysis systems1. Currently, weather and climate models are subjected to having various systematic errors, which induce large spread in results of different climate models2 and have a significant impact on the local and mid-latitude prediction skill3. Also, contemporary global reanalysis products, often used in climate research in order to monitor climate change and variability, or to evaluate climate model systems4, have large biases in temperature, humidity and wind speed in the Arctic5,6,7.

Specifically, in the Arctic surface energy balance significant deviations between climate models are found2,4. Within the coupled atmosphere, ocean, snow, and sea-ice system, this energy balance controls the growth and melt of sea-ice, as well, as the thermal stratification of the lower atmosphere. Thus, the accurate simulation of the Arctic surface energy budget is key for improving our understanding of the rapid climate change in the Arctic, as well as, advancing long term prediction of sea-ice properties in the future.

During Arctic winter clear-sky events (CSE) over sea-ice, with the absence of solar radiation and the strong longwave radiative cooling, the air temperatures can drop to −40 °C. Studies of the winter Arctic surface energy budget show that the radiative cooling of the troposphere is balanced by the advection of heat from lower latitudes8,9. In addition, in sea-ice covered regions the ocean moderates the low temperatures during CSEs in contrast to the land surfaces, where minimum temperatures of −50 to −60 °C occur8.

In the present study, we investigate the ability of contemporary regional and global weather prediction systems and global reanalysis systems to simulate the low temperatures observed during CSEs. We use the in-situ data set from the N-ICE 2015 drift campaign10 and a pan-Arctic sea-ice surface temperature satellite product. We find a warm bias in almost all analysed model systems and show that it is induced by a missing representation of snow on sea-ice. Furthermore, the simplistic representation of sea-ice thickness and concentration in the analysed model systems contribute to inconsistencies in the simulation of sea-ice surface temperatures.

Results

Overview of the studied models and utilised observations

We use the observational in situ data set of atmospheric, snow, sea-ice and ocean observations which is available from the N-ICE 2015 campaign and taken during the four subsequent drifts of a research vessel in between January and June 2015 (Fig. 1). The meteorological conditions, the thermodynamic structure of the troposphere, and the surface energy budget during the campaign are analysed in a set of studies11,12,13. Several CSE have been observed during N-ICE campaign13. In addition, we utilise a pan-Arctic sea-ice surface temperature satellite product14 based on infrared data from AVHRR instruments in order to extend our results to the pan-Arctic scale. More details on the use of the observations are given in the Methods section.

Fig. 1
figure 1

General overview of the sea-ice cover, experimental model domains and N-ICE 2015 campaign drifts. Background: mean sea ice concentration in January–February 2015 according to OSI-SAF data; A: model domain for the AA (control) experiment, B: model domain for the AA-snow (sensitivity) experiment; coloured curves: trajectories of the different N-ICE 2015 campaign legs; inset: trajectories of N-ICE 2015 drifts and restored true colour image retrieved from MODIS radiances showing the ice cover on 3 April 2015. Leg 1: 15 January–21 February; leg 2: 24 February–19 March; leg 3: 18 April–5 June; leg 4: 7 June–22 June

We investigate the ability of regional and global weather prediction systems (HARMONIE-AROME15 configuration of the ALADIN-HIRLAM numerical weather prediction system, and IFS-HRES16), and global reanalysis systems (ERA-Interim17, ERA5, MERRA-218, JRA-5519, and NCEP-DOE Reanalysis 220) to simulate those low temperature events. For HARMONIE-AROME two configurations of the sea-ice parameterisation scheme are studied. First, the control experiment (AA) which resembles sea-ice handling in IFS-HRES and has a snow-free sea-ice layer of uniform and fixed thickness. Second, the sensitivity experiment (AA-snow) which uses extended sea-ice parameterisation scheme that includes a prognostic sea-ice thickness and snow layer model (see the Methods section for further details). General overview of the model systems discussed in the present study is provided in the Table 1.

Table 1 Overview of the numerical weather prediction systems and global reanalysis systems

Analysis of a clear-sky event from N-ICE 2015

During the winter drifts of the N-ICE 2015 campaign a large number of storms has been observed11. Since those storms are connected to the advection of warm air onto the sea-ice, they can be identified in the temperature timeseries by relatively high temperatures close to 0 °C (Fig. 2). In between those storms, periods of clear-sky conditions occurred, marked by reduced longwave downwelling radiation flux (LW\(\downarrow\)) of around \(160\ {\mathrm{W}}{{\mathrm{m}}}^{-2}\) and low temperatures of around −30 to −40 °C. We refer in the following to three clear-sky events, CSE-1, CSE-2, and CSE-3, which are marked in Fig. 2.

Fig. 2
figure 2

Observations from the N-ICE 2015 campaign compared to various reanalysis products and model experiments. a Difference between the skin temperature predicted by different models and skin temperature observed during the N-ICE 2015 campaign (background); time series of the measured skin temperature (foreground curve). b Same as a but for downwelling longwave radiation

During the CSE-1 all models showed low LW\(\downarrow\) values of 140 to \(150\ {\mathrm{W}}{{\mathrm{m}}}^{-2}\), which are consistent with the observations (see the Fig. 2) and indicate that clear-sky conditions are well simulated. However, the simulated surface temperatures (\({T}_{{\rm{skin}}}\)) are much higher than observed. While the observations show values of down to \(-4{0}^{\circ }\ {\rm{C}}\), most of the models simulate temperatures of 5 to \(1{5}^{\circ }\ {\rm{C}}\) too warm, with the strongest deviations found in the snow-free AA experiment. The sensitivity experiment AA-snow, with a prognostic sea-ice thickness and snow layer model, is the only model simulation with small deviations. For all model products, the \({T}_{{\rm{skin}}}\) is obtained from the LW\(\downarrow\) and LW\(\uparrow\), while LW\(\uparrow\) corresponds to the sea-ice covered part of the grid cell only (see also Eq. (6) in the Methods section).

During the CSE-2 and CSE-3, the models show a similar warm temperature bias as during the CSE-1, but with differences in their capabilities in simulating the low LW\(\downarrow\), which are characteristic for clear-sky conditions. During CSE-2, MERRA-2 and ERA5 show too high LW\(\downarrow\) indicating that the clear sky conditions are not well resolved. During CSE-3, the ECMWF products (ERA5, ERA-I, and IFS-HRES) and MERRA-2 do not show the resolved clear-sky conditions and provide too high LW\(\downarrow\). For this event, the temperature biases are largest in the ECMWF products with the values of up to 15 °C.

In general, the models show a consistent warm bias during the polar winter clear-sky conditions over sea-ice, when observed temperatures are usually below −25 °C. In order to better understand this discrepancy between the modelled and observed temperatures, we analyse the surface energy budget (see also the Methods section) during the CSE-1, where all models show a consistent representation of the clear sky conditions in terms of the simulated low LW\(\downarrow\) values. Averaged values of the surface energy budget are given in Table 2. The conductive heat flux \(C\) is calculated (assuming steady linear sea-ice temperature profile) by using the respective snow and sea-ice thickness properties of the models (Table 1) or, in case of AA-snow, the directly simulated daily values of the sea-ice and snow thickness. Note that, for the sake of convenience, we assumed a 1.5-m sea-ice thickness for the MERRA-2 reanalysis. We tested the sensitivity of the flux by also prescribing sea-ice thicknesses of 1 and 2 m, which resulted in conductive heat fluxes ranging from 30 to \(61\ {\mathrm{W}}{{\mathrm{m}}}^{-2}\), i.e., in the range of the values obtained from the other model systems. In order to compute the observed conductive heat flux the mean observed values for sea ice and snow thickness are used21, i.e., \({h}_{{\rm{s}}}=0.5\ {\mathrm{m}}\) and \({h}_{{\rm{i}}}=1.4\ {\mathrm{m}}\).

Table 2 Surface energy budget–Clear sky event 1

The conductive heat flux in the observations is \(16\ {\mathrm{W}}{{\mathrm{m}}}^{-2}\). The AA-snow experiment shows a conductive heat flux consistent with the observations. The AA experiment shows a much higher flux of \(61\ {\mathrm{W}}{{\mathrm{m}}}^{-2}\), which is due to the combination of a missing snow layer and a small prescribed sea-ice thickness of 0.75 m. The global reanalyses, which do not represent snow over sea-ice (NCEP-2 has a rudimentary parameterisation for snow over sea-ice, but it could be neglected. See Table 1 and Methods for further details) but have a realistic assumption of the sea-ice thickness of 1.5 to 2 m, still overestimate the conductive heat flux by a factor of more than two, with values between 28 and \(40\ {\mathrm{W}}{{\mathrm{m}}}^{-2}\). The JRA-55 and NCEP-2 have the smallest values, due to their relatively large sea-ice thickness of 2 m.

The overestimation in \(C\) from the ocean to the surface is to a large extent compensated by an increased LW\(\uparrow\). The observed value for LW\(\uparrow\) is \(182\ {\mathrm{W}}{{\mathrm{m}}}^{-2}\), while the global reanalyses show values of around \(200\ {\mathrm{W}}{{\mathrm{m}}}^{-2}\) and in the AA experiment the flux reaches \(220\ {\mathrm{W}}{{\mathrm{m}}}^{-2}\). Again, AA-snow is consistent to the observations with a value of \(179\ {\mathrm{W}}{{\mathrm{m}}}^{-2}\).

We conclude that the missing snow layer on top of the sea-ice results in the overestimated conductive heat flux from the ocean to the atmosphere, which is compensated by an increase in outgoing longwave radiation and too high surface temperatures. The resulting bias in the net surface radiation budget (\({\mathrm{LW}}{\uparrow} +\,{\mathrm{LW}}{\downarrow}\)) is about 20 to \(30\ {\mathrm{W}}{{\mathrm{m}}}^{-2}\).

A pan-Arctic view of the temperature bias

From infrared satellite observations we can estimate the pan-Arctic \({T}_{{\rm{skin}}}\) during winter clear-sky conditions in the selected period from 2015 to 2017. In the satellite observations the mean clear-sky surface temperatures are as low as −35 °C in large parts of the central Arctic and towards the Canadian coastline (Fig. 3a).

Fig. 3
figure 3

Satellite ice surface temperature product compared to various reanalysis products. a Climatological winter-time clear-sky ice surface temperature in CMEMS for years 2015–2017; b Climatological snow thickness according to the TOPAZ reanalysis for years 2015–2017; cf difference between the climatological winter-time clear-sky ice surface temperature in different reanalysis products and the climatological winter-time clear-sky ice surface temperature in CMEMS for years 2015–2017, (c) ERA5, (d) JRA-55, (e) MERRA-2, (f) NCEP-2

In the coldest areas, all reanalyses have a warm bias of about 5 to 10 °C (Fig. 3e–f). ERA5 (Fig. 3c) and MERRA-2 (Fig. 3e) have similar spatial characteristics for the surface temperature with both showing a warm bias in all sea-ice covered areas and MERRA-2 having the largest temperature biases. The JRA-55 (Fig. 3d) and NCEP-2 (Fig. 3f) reanalyses have a smaller warm bias in the central Arctic and tend to have a cold bias in the areas further away from the North pole. This cold bias is about 5 to 10 °C and is most pronounced in the Northern Barents Sea, the areas from Hudson Bay towards Baffin Bay and some parts along the Russian, Canadian and American coastlines. Note, here we derive the model simulated surface skin temperature from the LW\(\downarrow\) and LW\(\uparrow\), following the Eq. (1) in the Methods section. This is equivalent to the grid cell mean temperature of the model, and not the surface temperature of the sea-ice part of the grid cell only, and thus consistent to the infrared surface observation.

The pattern of the \({T}_{{\rm{skin}}}\) warm bias in the global reanalyses shows similar characteristics to the retrieved mean snow depth from the TOPAZ4 ocean and sea-ice reanalysis22 during CSEs in 2015 to 2017 (Fig. 3b). For example the lower snow-depths in Baffin Bay and Labrador Sea result in higher observed \({T}_{{\rm{skin}}}\) and a smaller bias, while the largest snow-depths in the central Arctic and towards the Canadian Archipelago are co-existing with the lowest temperatures and the largest temperature biases in the models.

Our findings of a temperature bias during clear sky conditions in the satellite analysis are consistent to results of the N-ICE 2015 drift analysis in the previous section. In order to further strengthen the conclusion that the snow representation is the key factor, the surface temperature bias can be formulated as a function of the misrepresentation of snow and the sea-ice thickness (see Methods section and Fig. 4). If the prescribed sea-ice thickness is 1.5 m, an error in the snow thickness of only 0.25 m can induce a warm bias of 8 °C. This could be seen in ERA5 and MERRA-2 over the central Arctic (see Fig. 3c and e) where largest warm biases of the surface temperature are in a good agreement with the pattern of the mean snow depth (Fig. 3b). A misrepresentation in the sea-ice thickness can also have a strong effect on the surface temperature in cases of only a thin snow layer. The JRA-55 and NCEP-2 reanalyses have a 2 m sea-ice thickness representation, compared to 1.5 m in ERA5 and MERRA-2. If we assume that there is a 0.5 m overestimation in sea-ice thickness in regions with no snow cover, it would lead to a cold bias of about −5 °C, and counter-balance the warm bias induced by the missing snow in areas with a snow layer of around 0.1 m (Fig. 4). This is consistent with the relatively small warm bias in central parts of the Arctic in JRA-55 and NCEP-2, and agrees with the cold bias further away from the North pole in these reanalyses (see the Fig. 3d and f). In addition, JRA-55 and NCEP-2 have a binary representation of the sea-ice cover. When sea-ice concentration is above 55% (or 50% for NCEP-2) these reanalyses assume that sea is completely ice-covered. This induces the additional spurious insulation from the ocean and, in turn, a cold surface temperature bias in regions closer to the sea-ice edge.

Fig. 4
figure 4

The dependency of the surface temperature bias on ice and snow thickness errors. The bias \(\Delta {T}_{{\rm{s}}}\) is induced by the errors of ice \(\Delta {h}_{{\rm{i}}}\) and snow thickness \(\Delta {h}_{{\rm{s}}}\). Here, we assume a simulated 1.5 m thick sea-ice without a snow layer, i.e., \({h}_{{\rm{i}}}=1.5\ {\mathrm{m}}\) and \({h}_{{\rm{s}}}=0\ {\mathrm{m}}\)

Discussion

An evaluation of the surface energy budget of current global and regional model configurations, which are used for weather forecasting and reanalysis products, highlights the importance of the simulation of snow on sea-ice. Three CSEs occurred during polar night are analysed using the N-ICE 2015 drift campaign data. The analysis shows that the models in their standard configurations are 5 to 15 °C too warm during these events. In the CSE-1 the low incoming longwave radiation in the model products is low, which indicates that the clear-sky conditions are reasonably well resolved. However, the outgoing longwave radiation is too high, as the consequence of a warm bias in the surface temperature. In turn, this bias can be attributed to the overestimation of the conductive heat flux from the ocean to the surface by 20 to \(40\ {\mathrm{W}}{{\mathrm{m}}}^{-2}\). This heat flux is regulated by the difference between the ocean and atmospheric temperatures and by the thickness of the sea-ice and snow layers. A snow layer, which has a thermal conductivity about seven times lower than that of sea-ice, is missing in most of the analysed models, and thus, is identified as the main reason for the deviations in the surface energy budget.

The warm bias in the ice surface temperatures is a characteristic feature of the analysed global reanalysis products on the pan-Arctic scale. By comparing to a multiyear product obtained from infrared satellite measurements we show that all tested systems have issues in simulating the extreme low temperatures over large parts of the Arctic. The most drastic deviations in the temperature are found over the areas with the thick snow cover according to the winter snow climatology (see the Fig. 3b–f). This strengthens the argument that the snow layer is strongly connected to the temperature bias found in the global atmospheric reanalyses. In two of the reanalyses, which have thicker sea-ice (JRA-55 and NCEP-2), the warm bias is reduced, but instead a cold bias is produced, specifically in areas where sea-ice is thinner than the prescribed 2 m and only a thin snow layer exists. In general, the characteristics of the surface temperature in the reanalysis products can be explained by the surface energy budget. We conclude that the snow component on the sea-ice improves the surface atmospheric energy budget in cold atmospheric conditions and thus is an important but often missing component in state-of-the-art reanalysis and forecasting systems. In addition, sea-ice thickness plays an important role, as well, but mainly in the areas where the snow layer is thin.

The surface energy budget is an integral part of many climate processes in the Arctic. For example, it determines the available energy for sea-ice melting and freezing, and also the thermal stratification in the lower troposphere. Thus, the accuracy of its representation can have a strong impact on the skill of climate prediction and our understanding of large-scale climate dynamics. The bias in the surface energy budget due to the misrepresentation of the snow and sea-ice layer is about \(20\ {\mathrm{W}}{{\mathrm{m}}}^{-2}\). Compared to the spread of 20 and \(60\ {\mathrm{W}}{{\mathrm{m}}}^{-2}\) between the various climate models in the winter-time net longwave radiation budget2,4, the temperature bias, which is induced by the missing snow layer and misrepresentation of sea-ice thickness, is of significant importance. This bias in the net radiation budget can induce wrong conclusions, e.g., assuming that the climate models have a “cold temperature” bias in winter4, or generally refine process studies which rely on the atmospheric reanalyses 9,23.

Due to the widespread use of the atmospheric reanalysis products for model validation, initialisation of prediction systems, forcing of ocean and sea-ice reanalyses24, etc., it is very important to take into account this temperature bias in the contemporary reanalyses. Future reanalyses could be enhanced primarily by considering a prognostic parameterisation scheme to represent snow cover over the sea-ice, and by improving the spatial characteristics of the sea-ice thickness.

Methods

Observations

We use the observational data set from the N-ICE 2015 drift north of the Svalbard Archipelago10. The N-ICE 2015 campaign consisted in total of four legs. In the present study we focus on the first two legs within the Arctic winter period, between 15 January and 19 March (Fig. 1), and analyse the meteorological observations of longwave radiation11. The skin temperature of the snow surface is obtained from the longwave upwelling LW\(\uparrow\) and downwelling LW\(\downarrow\) radiation by

$${T}_{{\rm{skin}}}={\left(\frac{{\mathrm{LW}}{\uparrow} -(1-{\varepsilon }_{{\rm{s}}}){\mathrm{LW}}{\downarrow} }{{\varepsilon }_{{\rm{s}}}\sigma }\right)}^{\frac{1}{4}},$$
(1)

where \(\sigma\) is the Stefan–Boltzmann constant and \({\varepsilon }_{s}\) is the emissivity of the snow, which we assume to be 0.9913.

In addition to the in-situ observation of N-ICE 2015, we include a satellite data set of daily sea-ice surface temperature (\({T}_{{\rm{skin}}}\)). The data product is based on infrared data from the AVHRR instrument on board the Metop-A satellite. The \({T}_{{\rm{skin}}}\) satellite data set is available from 2015 to present, as part of the Copernicus Marine Environmental Monitoring Service (CMEMS). Since satellite measurements of \({T}_{{\rm{skin}}}\) are only meaningful in case of cloud free conditions, for the satellite retrieval a cloud-mask has been used together with a bias-correction algorithm to obtain surface temperature estimates within cloud covered areas14. In the present study, however, we excluded those corrected surface temperatures, due to our focus on CSEs in winter. We define the CSEs by considering only the data points where: incoming longwave radiation in ERA5 is lower than \(160\ {\mathrm{W}}{{\mathrm{m}}}^{-2}\), the OSI SAF Ice Edge product identifies the grid cell as closed ice, there are more than eight AVHRR infrared observations in the grid cell, and the date is within the winter season (October to February). The same mask is then applied to the daily surface temperature values, in order to derive the mean \({T}_{{\rm{skin}}}\) of reanalysis and observations during clear sky conditions on a pan-Arctic scale.

Atmospheric models and global reanalysis products

We analyse a suite of global reanalyses, i.e., two configurations of the ECMWF model system (ERA-Interim25 and ERA5), the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-218), the Japanese 55-year Reanalysis (JRA-5519), and the National Centers for Environmental Prediction Reanalysis 2 (NCEP-220).

We also include the operational ECMWF global deterministic weather forecasting system (IFS-HRES16), as well as, two experimental configurations of the regional convective scale model HARMONIE-AROME15. The first of these configurations (AA) covers the same domain as the operational numerical weather prediction system AROME Arctic26 and has the same representation of the sea-ice cover. The second experimental configuration of HARMONIE-AROME (AA-snow) is similar to AA, but incorporates the parameterisations for the sea-ice mass-balance and snow layer on top of the sea-ice. The model domain for AA-snow is smaller than that one in the AA experiment (see the Fig. 1) to reduce the computational costs and to focus on the areas with extensive sea-ice cover. The AA experiment covers the time period from 10 January to 20 March 2015. The AA-snow experiment has been performed from 1 September 2014 to 20 March 2015. To allow realistic evolution of the snow cover in AA-snow, this experiment was started without snow layer on top of sea ice and snow was accumulated from the modelled precipitation. In both AA and AA-snow the initial ice thickness is 0.75 m but in AA it remains constant during the experiment and in AA-snow it evolves according to the ice mass-balance parameterisation.

Relevant features of the models, which were used to produce forecast and reanalysis products, are given in Table 1 and with a focus on the representation of the sea-ice surface in the Arctic. In short, none of the global model systems has a prognostic parameterisation of the sea-ice mass-balance nor a parameterisation of a snow layer on top of the sea-ice, which impacts the surface energy budget. NCEP-2 simulates snow on sea-ice, however, without any influence on the heat conductivity and surface albedo. In all models sea-ice concentration is updated daily by utilising satellite products. Operational configurations of HARMONIE-AROME do not simulate the sea-ice concentration and sea-ice thickness, but simulate sea-ice surface temperature26,27.

To extend the existing sea-ice parameterisation of HARMONIE-AROME with ice mass-balance calculations, we added the representation of the following processes: ice growth and melting from the bottom, and surface melting. The more complex processes of internal melting and snow-ice formation were not considered.

The interface between the sea-ice bottom and the underlying water body is governed by the balance between the sea-ice and ocean heat fluxes. When these fluxes are not in equilibrium state, the residual heat flux leads to growth or melting of the ice layer28:

$${\rho }_{{\rm{i}}}{L}_{{\rm{f}}}\frac{{\mathrm{d}}{h}_{{\rm{i}}}}{{\mathrm{d}}t}={k}_{{\rm{i}}}\left.{\frac{\partial {T}_{{\rm{i}}}}{\partial z}}\right\vert_{z={h}_{{\rm{i}}}}-{F}_{{\rm{w}}}$$
(2)

where \({\rho }_{{\rm{i}}}\) is the density of sea-ice, \({L}_{{\rm{f}}}\) is the latent heat of fusion, \({h}_{{\rm{i}}}\) is the total ice thickness, \({k}_{{\rm{i}}}\) is the thermal conductivity of sea-ice, and \({F}_{{\rm{w}}}\) is the ocean heat flux. For the AA-snow experiment, the ocean heat flux was assumed to be constant with the value of \(2\ {\mathrm{W}}{{\mathrm{m}}}^{-2}\).

The processes of surface melting is parameterised in the following way. At the first step, the ice surface temperature is calculated from the thermal balance of the surface layer of sea-ice. Then, if the obtained temperature is higher than the melting temperature of sea-ice, it is set to be equal this melting temperature and the residual heat flux \({F}_{{\rm{melt}}}\) induces the melting of the ice surface28:

$${F}_{{\rm{melt}}}=-{\rho }_{{\rm{i}}}{L}_{{\rm{f}}}\frac{{\mathrm{d}}{h}_{{\rm{i}}}}{{\mathrm{d}}t}$$
(3)

The AA-snow experiment, which uses the updated sea-ice parameterisation scheme, shows simulated snow depths between 0.45 and 0.7 m and sea-ice thickness between 1.4 and 1.7 m. This is consistent with the observed snow and sea-ice thickness of \(0.5\pm 0.2\ {\mathrm{m}}\) and \(1.4\pm 0.3\ {\mathrm{m}}\), respectively21.

We compare the model’s \({T}_{{\rm{skin}}}\) with the satellite and in-situ observations described in the previous sections. For the satellite observations \({T}_{{\rm{skin}}}\) corresponds to the model’s grid-cell mean temperature, and in case of the in-situ observations it corresponds to the \({T}_{{\rm{skin}}}\) of the sea-ice covered part only. Thus, for the satellite observation comparison we can use the Eq. (1). For the comparison with the in-situ observations we partition the longwave upwelling radiation LW\(\uparrow\) into an ocean and sea-ice covered part according to the sea-ice concentration \({\mathrm{SIC}}\) at the respective grid-cell.

$${\mathrm{LW}}{\uparrow} =({\mathrm{1}}-{\mathrm{SIC}})\cdot {\mathrm{LW}}{\uparrow }^{({\rm{ocean}})}+{\mathrm{SIC}}\cdot {\mathrm{LW}}{\uparrow }^{({\rm{sea-ice}})}$$
(4)

If we assume the ocean surface with temperature of \(-{2}^{\circ }\ {\rm{C}}\) and emissivity of 0.98, then by following the Stefan-Boltzman law we obtain

$${\mathrm{LW}}{\uparrow }^{({\rm{ocean}})}=0.98\cdot \sigma \cdot 271.1{5}^{4}\cdot$$
(5)

Thus, combining the Eqs. (1), (4), and (5), the skin temperature of the ice covered part of the grid-cell is derived as follows:

$${T}_{{\rm{skin}}}^{({\rm{sea-ice}})}={\left[\frac{{{\mathrm{SIC}}}^{-1}\left({\mathrm{LW}}{\uparrow} -({\mathrm{1}}-{\mathrm{SIC}}){\mathrm{LW}}{\uparrow }^{({\rm{ocean}})}\right)-({\mathrm{1}}-{\varepsilon }_{{\rm{s}}}){\mathrm{LW}}{\downarrow}}{{\varepsilon }_{{\rm{s}}}\sigma }\right]}^{\frac{1}{4}}$$
(6)

Surface energy budget and surface temperature

The energy budget of the surface layer, or the net-energy transfer \({F}_{{\rm{sfc}}}\) between the atmosphere and the ocean can be written as9,13

$${F}_{{\rm{sfc}}}={C}_{{\rm{t}}}\frac{\partial {T}_{{\rm{s}}}}{\partial t}-C=\left(1-{i}_{0}\right)\left(1-\alpha \right){\mathrm{SW}}{\downarrow} +{\mathrm{LW}}{\downarrow} +{\mathrm{LW}}{\uparrow} +Q$$
(7)

where \({C}_{{\rm{t}}}\) is the thermal resistance of the surface layer, \(C\) is the the conduction of heat from the ocean through the sea-ice/snow to the atmosphere, \({i}_{0}\) denotes the part of the downwelling shortwave radiation that penetrated through the surface layer, \(\alpha\) is the surface albedo, \({\mathrm{SW}}{\downarrow}\) is the surface downwelling short-wave radiation flux, LW\(\downarrow\) is the surface downwelling longwave radiation flux, LW\(\uparrow\) is the surface upwelling long-wave radiation flux, \(Q\) is the turbulent flux of sensible and latent heat. In the polar night conditions the downwelling shortwave radiation flux is very small and could be neglected from the Eq. (7). The conductive flux from the ocean to the snow surface is not directly observed. However, we can estimate it as

$$C=\frac{{k}_{{\rm{i}}}\cdot {k}_{{\rm{s}}}}{{k}_{{\rm{s}}}\cdot {h}_{{\rm{i}}}+{k}_{{\rm{i}}}\cdot {h}_{{\rm{s}}}}({T}_{{\rm{o}}}-{T}_{{\rm{skin}}})$$
(8)

where \({T}_{{\rm{o}}}\) is the surface ocean temperature, \({k}_{{\rm{s}}}\) is the thermal conductivity of snow, and \({h}_{{\rm{s}}}\) is the thickness of the snow layer on top of sea-ice. For the thermal conductivity of snow and sea-ice we assume, \({k}_{{\rm{s}}}=0.31\ {\mathrm{W}}{{\mathrm{m}}}^{-1} \, {{\mathrm{K}}}^{-1}\) and \({k}_{{\rm{i}}}=2.1\ {\mathrm{W}}{{\mathrm{m}}}^{-1} \, {{\mathrm{K}}}^{-1}\), respectively. Note that Eq. (8) implies steady linear temperature profile within the sea-ice layer. This assumption is not generally correct in case of rapid changes of weather conditions or multilayer sea-ice schemes, especially when snow and ice layers are of considerable thickness. However, taking into account minor variability of weather conditions during a single polar night CSE, Eq. (8) could provide a simple first-order estimate of the real conductive heat flux.

If we assume that a variation \(\Delta R\) in the net longwave radiation budget \(R={\mathrm{LW}}{\downarrow} +{\mathrm{LW}}{\uparrow}\) balances the changes in the conductive heat flux \(\Delta C\) induced by errors in snow thickness \(\Delta {h}_{{\rm{s}}}\) and sea ice thickness \(\Delta {h}_{{\rm{i}}}\), we can write \(\Delta R=\Delta C\). For the radiation we yield

$$\Delta R=R({T}_{{\rm{s}}}+\Delta {T}_{{\rm{s}}})-R({T}_{{\rm{s}}})=\varepsilon \sigma \cdot {({T}_{{\rm{s}}}+\Delta {T}_{{\rm{s}}})}^{4}-\varepsilon \sigma {T}_{{\rm{s}}}^{4}\approx 4\varepsilon \sigma {T}_{{\rm{s}}}^{3}\cdot \Delta {T}_{{\rm{s}}}$$
(9)

with \(\Delta {T}_{{\rm{s}}}\) being the induced temperature bias. And the variation of the conductive heat flux is defined as

$$\Delta C=C({T}_{{\rm{s}}}+\Delta {T}_{{\rm{s}}},{h}_{{\rm{i}}}+\Delta {h}_{{\rm{i}}},{h}_{{\rm{s}}}+\Delta {h}_{{\rm{s}}})-C({T}_{{\rm{s}}},{h}_{{\rm{i}}},{h}_{{\rm{s}}})$$
(10)
$$=\alpha {\prime} \cdot ({T}_{{\rm{o}}}+\Delta {T}_{{\rm{s}}}-{T}_{{\rm{skin}}})-\alpha \cdot ({T}_{{\rm{o}}}-{T}_{{\rm{skin}}})$$
(11)

where

$$\alpha \equiv \frac{{k}_{{\rm{i}}}\cdot {k}_{{\rm{s}}}}{{k}_{{\rm{s}}}\cdot {h}_{{\rm{i}}}+{k}_{{\rm{i}}}\cdot {h}_{{\rm{s}}}}$$
(12)
$$\alpha {\prime} \equiv \frac{{k}_{{\rm{i}}}\cdot {k}_{{\rm{s}}}}{{k}_{{\rm{s}}}\cdot ({h}_{{\rm{i}}}+\Delta {h}_{{\rm{i}}})+{k}_{{\rm{i}}}\cdot ({h}_{{\rm{s}}}+\Delta {h}_{{\rm{s}}})}$$
(13)

Thus, we derive the temperature bias induced by changes in snow and sea-ice thickness during clear sky conditions as

$$\Delta {T}_{{\rm{s}}}=\frac{(\alpha {\prime} -\alpha)\cdot ({T}_{{\rm{s}}}-{T}_{{\rm{o}}})}{4\varepsilon \sigma {T}_{{\rm{s}}}^{3}-\alpha {\prime} }.$$
(14)