Arctic Sea Ice Decline in the 2010 s : The Increasing Role of the Ocean — Air Heat Exchange in the Late Summer

This study is focused on the specific features of ocean–air interaction in the Laptev Sea, in the late summer, on the basis of recurrent measurements during four expeditions in the 2000s and 2010s, atmospheric reanalysis products, and satellite ice concentration data. It was established that in the “icy” years, the accumulation of heat in the upper ocean layer is insignificant for the subsequent ice formation. In the “ice-free” years, the accumulated heat storage in the upper mixed layer depends on the duration of open water and the distance of the point of interest to the nearest ice edge. In a broader context, we considered possible links between the average ice area/extent in the August–September–October (ASO) period, and in the December–January–February (DJF) period, for two representative Arctic regions; that is, the Eurasian segment, defined within the bounds 60–120◦ E, 65–80◦ N, and the American segment, defined within the bounds 150◦ E–150◦ W, 65–80◦ N. Significant “seasonal memory”, characterized by the consistent change of the ice cover parameters in sequential seasons, was revealed in the Eurasian segment between 2007 and 2017. No linkage on a seasonal time scale was found in the American segment. A possible explanation for the distinguished contrast between the two geographical regions is proposed.


Introduction
The steady poleward advance of the open water area at the end of summer is in line with the existing climatic projections of the progressive transition to seasonal ice cover in the Arctic Ocean by the end of the 21st century [1].The year 2007 could be considered as a sort of tipping point between the prevailing multiyear and seasonal sea ice in the Siberian Arctic seas [2].Hydrometeorological processes in the seasonally ice-covered seas substantially differ from the seas under permanent ice cover, because of the enhanced heat exchange between the atmosphere and the upper water layer, a distinct freeze-thaw cycle, and a deeper penetration of various types of mixing.
There is growing evidence that a new state of the hydrographic regime is currently shaping under the rapidly changing conditions at the ocean surface in the East Atlantic sector of the Arctic Ocean [3][4][5][6].It is expected that this new state will be characterized by warmer and saltier well-mixed upper ocean layers, instead of a quasi-two-layer structure with a cold and fresh upper layer under the permanent ice cover.A general north-eastward shift of such structural changes has already been documented in the Barents Sea and the Western Nansen Basin, where the main driver of change is the warm Atlantic Water (AW) inflow, where its influence on the upper ocean has intensified due to the reduction of ice cover [3,5].Further to the east, the warm AW becomes sheltered from the ocean surface by high gradient transient layers, thereby losing an efficient means of delivering heat and salt upwards.In these regions, the possible driver of change is the excessive absorption of short-wave solar radiation under the conditions of increased duration of open water in summer.Whether this accumulated extra heat is rapidly lost to the atmosphere during the next fall, or the "memory" regarding the warmer surface layer is transferred to the first year ice properties, as suggested by [7], depends on multiple factors.The extended size of the ice-free area facilitates massive ice growth as soon as the surface water temperature reaches the freezing point.On the other hand, a longer duration of ice-free areas delays the onset of freezing, since more atmospheric heat can be accumulated in the upper mixed layer (UML).The rate of heat absorption depends on the prevailing atmospheric circulation patterns and the local weather conditions.Downward propagation of the absorbed heat in the water column, which ultimately determines the total heat content of the upper ocean at the end of the warm season, is controlled by density stratification and atmospheric dynamical forcing [8].This heat has to be released to the air before the freezing can start.The rate of the UML cooling in the fall primarily depends on turbulent heat flux at the ocean-air interface, because during the polar night, radiation balance is permanently negative.
In this study, we validated the outlined theoretical considerations using field measurements and atmospheric reanalysis products.We focused our investigation on the central Laptev Sea, where recurrent hydrological transects and accompanying meteorological measurements were carried out in September 2003September , 2005September , 2013, and 2015, as a part of the ongoing NABOS (Nansen and Amundsen Basins Observations System, http://nabos.iarc.uaf.edu) research program.
The paper contains four sections, including the current one.In Section 2, we describe the data and methodology used.Section 3 is devoted to the joint analysis of the available data in the context of the basic ideas formulated above.In the discussion, we extend an outcome of the performed analysis for the tentative explanation of the interannual variability in the sea ice properties in the representative regions of the Arctic Ocean, observed during the 2000-2010s.General conclusions are formulated in the final section.Supporting information on the observational and data processing techniques, as well as the comparison between the in-situ observations and reanalysis data, is presented in the Appendix A.

Data and Methods
Field data used in this study were collected during four research expeditions to the Arctic Ocean, carried out in August-September onboard the icebreaker "Kapitan Dranitsyn" (2003 and 2005), and the research vessels "Akademik Fedorov" (2013) and "Akademik Tryoshnikov" (2015).The scope of the investigations included multidisciplinary studies of the physical, chemical, and biological conditions in the Eurasian Arctic, with the major focus on the deep part of the Laptev Sea (see NABOS cruise reports at http://nabos.iarc.uaf.edu/).Hydrographic measurements were carried out at the cross-slope sections using the SeaBird-911plus CTD (Conductivity Temperature Depth) device (https://www.seabird.com/profiling/sbe-911plus-ctd/family?productCategoryId=54627473769, manufactured by SeaBird Scientific, USA).CTD accuracy is ±0.001 • C for temperature, ±0.0003 S/m for conductivity, and ±0.001 PSU for salinity.We used the vertical gradient of the potential density 0.01 kg•m 3 •m −1 as the quantitative criterion for determination of the UML depth.
Atmospheric and sea ice data used in this study included the ERA-Interim reanalysis from ECMWF (European Centre for Medium-Range Weather Forecasts).[9] and the DMSP (Defense Meteorological Satellite Platform) digitized passive microwave satellite images [10].These products are freely available on the internet.The ERA-Interim reanalysis is one of the most popular global datasets for climate studies, especially for data-sparse regions.For the Arctic region it is still not perfect [11], but it provides the best verification results for wind speed, near-surface temperature, and radiative fluxes in comparison to six other global reanalysis products, as is shown in a detailed evaluation study [12].
One of the key analyzed variables in our study is the total heat balance for the ocean surface, B, and its components: where Q S and Q L are the net shortwave and longwave radiation, and F is a total turbulent heat flux.
Since, in our study, the turbulent heat exchange at the ocean-atmosphere interface is considered relative to the ocean, the positive F values indicate heat transfer from the atmosphere to the ocean: where H and LE are, respectively, the sensible and latent turbulent heat fluxes, which are directed from the ocean to the atmosphere according to their classic definition in atmospheric sciences.The methods of calculation of the turbulent fluxes according to in-situ ship-based observations are presented in the Appendix A.
During the research cruises, meteorological measurements were permanently conducted on the route.The list of instruments used and the measured parameters is presented in the Appendix A, Table A1.In 2013 and 2015, all components of heat balance at the sea surface were calculated from the shipborne measurements.In 2005, the radiation balance was reconstructed from the ERA-Interim reanalysis data.In 2003, all components of heat balance were reconstructed from the ERA-Interim reanalysis data.For consistency, we based the analysis of meteorological conditions on ERA-Interim products in all four years.The quality of the reanalysis data was confirmed by comparison to shipborne measurements in 2005, 2013, and 2015 (see Figure A5 in the Appendix A).In cases where the discrepancy between the shipborne measurements and the ERA-Interim data appeared to be too large, we based our conclusions on the shipborne measurements, considering them to be more reliable.

Ice and Hydrometeorological Conditions in the Central Laptev Sea in Summer
Sea ice extent (the edge of the ice concentration exceeding 15%) at the peak of the summer minimum is the conventional indicator of interannual changes in the state of the Arctic sea ice cover.This indicator, however, does not consider ice thickness (age).In this sense, the four considered years belonged to different ice regimes in the Laptev Sea, that is, the dominance of perennial ice up to 2005 and the substantially increased fraction of seasonal ice in 2006 and after (see details in Discussion section).The mean September sea ice extent in the Laptev Sea according to the DMSP satellite data for four studied years is presented in Figure 1.The position of the recurrent CTD stations is shown by colored squares.Supportive additional information is given in the inset and in Table 1.

Ice Conditions
A general northward shift of the ice edge is clearly seen from the 2000s to 2010s (about 300 km between 2003 and 2015).However, the difference between 2005 and 2013 is rather small, pointing out the significant year-to-year variations.In 2003, the position of the recurrent CTD station was permanently in the ice covered waters, while in the other years, the position was in the open water for different periods of time (see the inset in Figure 1).In 2005, the open water first reached the point of interest on August 25 ± 5 days.However, for the next 20 days, the distance to the ice edge did not exceed 50 km.Considering the 10 day averaging of the used DMSP product and the relatively large error of edge determination in the melting season [13], we may expect that in 2005, between August 15 and September 15, the point was very close to, or even within, the marginal ice zone (MIZ).153   The background thermohaline structure in the upper ocean in September, during four 159 considered years, was primarily controlled by the ice conditions.This expected linkage is illustrated 160

Hydrographic Conditions
The background thermohaline structure in the upper ocean in September, during four considered years, was primarily controlled by the ice conditions.This expected linkage is illustrated by the temperature and salinity profiles in the upper 40 m depth layer at the recurrent CTD stations in Figure 2. In 2003 and in 2005, water temperature in the entire 40-m layer was negative and did not exceed −1 • C. In 2003, the maximal departure from the freezing point, up to 0.5 • C, was observed at the base of the UML.In 2005, the maximal departure was less than 0.25 • C. On the contrary, in 2013 and 2015, the departure from the freezing point in the UML exceeded 1 degree during the second CTD cast, reaching 3 degrees in 2015.However, on the 5-14 day time scale, which corresponds to the time interval between recurrent CTD casts, the effect of ice conditions on the vertical structure was less unequivocal.For example, the evolution of the UML demonstrated definite similarities in specific years with contrasting ice conditions.In "icy" 2003 and in "ice-free" 2015, the thickness and salinity of the UML increased unidirectionally between the two recurrent CTD casts, reaching about the same value of 23 m during the second cast.In 2005, the UML was 4 m shallower during the second cast.In 2013, there was no distinct seasonal UML during the first cast at all, while at the second cast, it appeared within the 11 m depth.Temperature change in the UML between the two casts in the neighboring "ice-free" years was opposite-an increase in 2013 and a decrease in 2015, while corresponding salinity changes were unidirectional (increase) in both years.These similarities and differences in the observed short-term changes over the years, which were characterized by different and similar ice conditions, point to variable atmospheric forcing as a probable factor affecting the state of the upper ocean layer in the warm season.

Weather Conditions
Synoptic scale atmospheric forcing at the ocean surface is determined by weather conditions.
The distinctive feature of the Siberian shelf seas in general (and the Laptev Sea in particular) in summer is a close neighborhood (on the scale of the baric systems) of cold ice covered surfaces, a warm continent, and relatively warm open water in between.Hence, the intensity of the ocean-air energy exchange, and even its sign, primarily depends on the characteristics of the air masses that

Weather Conditions
Synoptic scale atmospheric forcing at the ocean surface is determined by weather conditions.The distinctive feature of the Siberian shelf seas in general (and the Laptev Sea in particular) in summer is a close neighborhood (on the scale of the baric systems) of cold ice covered surfaces, a warm continent, and relatively warm open water in between.Hence, the intensity of the ocean-air energy exchange, and even its sign, primarily depends on the characteristics of the air masses that enter the affected area.This feature is a principal difference between the Polar regions and the lower latitudes.
In the mid-latitudes, there is a distinct division between the summer season, when the colder ocean gains heat from the warmer air, and the fall-winter season, when the warmer ocean releases heat to the atmosphere.In the Polar regions, even in the mid-summer season, the upper ocean layer may act as a heat absorber, as well as a heat emitter, because of the near-zero sea surface temperature (SST) close to the ice edge under the ice-free conditions [14].An air flow, which floods over an ice-free surface, may be either warmer, if originated over the continent, or colder, when coming from the ice covered regions in the high-latitude Arctic.
Considering these possible options, we assessed the synoptic conditions over the four considered years according to the ERA-Interim reanalysis data.In 2003, the position of the recurrent CTD station was permanently at the periphery of the vast cyclone occupying almost the entire Arctic Ocean (Figure 3a).Western winds forced the ice edge to the east, keeping the point in the MIZ.In 2005, between September 14 and 20, the central Laptev Sea was located in a saddle-shaped sea level pressure (SLP) field between four large scale pressure systems, centered over the Arctic Ocean and the adjacent continent (Figure 3b).Weak SLP gradients caused intermittent winds and the slow transport of air masses.Within the entire time interval between the two recurrent CTD stations in 2013, the study area remained in the pathway of warm air flow from the continent between the anticyclone over the central Arctic Ocean, and the cyclone, centered in the Kara Sea (Figure 3c).The time interval between the recurrent CTD stations in 2015 was the longest compared to the other years, that is, two and a half weeks.The study area was mainly in the cyclonic circulation field.This caused the prevalence of western winds from the Kara Sea and an inflow of the air with near-zero • C surface air temperature (SAT) (Figure 3d).

Ocean-Air Energy Exchange
Time series of total heat balance at the area of the recurrent CTD station through August-September from the ERA-Interim data (daily average values) are shown in Figure 4 (in the figure, and hereafter in the text, we consider the total turbulent heat flux F, where positive values indicate a heat transfer from the atmosphere to the ocean.See Equations ( 1) and (2).Time series of the specific components of total heat balance along the ship route are presented in the Appendix A, Figure A5.A common feature for all four years was the gradual decrease of radiation balance from the highest values (up to 100 W/m 2 ) in early August to near-zero values at the end of September (Figure 4).Such a consistent pattern is obviously explained by the seasonal cycle of the short-wave radiation.Small-scale variations of radiation balance, superimposed on this general seasonal pattern, might be linked to variable cloudiness and the difference in long-wave radiation flux at the ice-free and ice covered underlying surface.
On the other hand, the difference in the pattern of turbulent heat fluxes between 2003, 2005, 2013, and 2015 for the entire two-month time interval was very visible.In 2003 and in 2005, turbulent heat fluxes were alternating around zero until the end of September.In 2013 and in 2015, intensive ocean heat losses due to turbulent exchange were observed during the entire second part of September (mean values of the total turbulent heat flux F over the period from September 15 to September 30 at the recurrent CTD site were −49 W/m 2 on average in 2013 and −82 W/m 2 in 2015, according to the ERA-Interim reanalysis), thus making the total heat balance at the surface substantially negative.
However, there was a large difference in the turbulent heat balance at the smaller time scale in 2013 and 2015.In 2013, the mean total heat balance between the two recurrent CTD casts was positive.Both data sets (direct measurements and ERA-Interim reanalysis) consensually pointed out that the main contributors to the total heat balance in 2013 were radiation components.In 2015, the turbulent heat flux patterns between the two recurrent CTD casts were more ambiguous.According to the ERA-Interim reanalysis, the negative and positive heat fluxes alternated, yielding a small average total heat flux (Figure 4d).On the other hand, direct measurements indicated a constantly negative total heat flux along the ship route with the dominance of the turbulent heat flux (Figure A5c).Hence, in 2015, we preferred to trust the shipborne measurements and cautiously concluded the decisive contribution of turbulent fluxes to the total heat balance.
Turbulent heat fluxes are controlled by the state of the underlying surface (ice of various type or an open water), the weather conditions (SAT and wind speed), and the SAT-SST difference.A closer look at the time series of the observed meteorological parameters (see the Appendix A, Figures A2, A3  and A5) confirmed our conclusions on the dominance of turbulent heat fluxes in the total heat balance of the ocean surface in 2015.In 2013, the SAT was more than 2 • C lower than that in 2015, according to ship-based measurements, while the observed SST values in the two years were practically similar (see Table 2 for details).Hence, the observed SAT-SST difference in the Laptev sea was close to zero in 2013 and −2.5 • C in 2015.Warmer air flow from the continent in 2013 resulted in a near-zero turbulent heat flux.In 2015, the relatively cold air inflow on the warmer (over 2 • C) ice-free surface provided negative turbulent heat fluxes.This resulted in generally negative heat balances at the ocean-air interface between the two recurrent CTD casts in 2015, instead of the positive heat balance in 2013.

Heat Accumulation in the UML
In this subsection, we consider the evolution of the heat storage in the UML according to the recurrent CTD casts with linkages to the atmospheric forcing and ice conditions.The heat storage  is defined relative to the freezing point temperature:

Heat Accumulation in the UML
In this subsection, we consider the evolution of the heat storage in the UML according to the recurrent CTD casts with linkages to the atmospheric forcing and ice conditions.The heat storage C is defined relative to the freezing point temperature: where ρ is the water density, c p is the water heat capacity, h is the UML depth, T is the mean temperature of the UML, and T f is the freezing temperature that depends on the salinity S.
Long duration of open water is the prerequisite condition for building up heat storage in the UML.In 2003, the studied area was permanently covered by pack ice of variable concentrations.Despite similar radiation balance values, as in the other years (see Figure 4), the temperature in the UML (15 m thick during the first CTD cast, see Figure 2 30.79 ÷ 30.82 to 30.65 PSU.These unidirectional changes of thermohaline parameters point out that incoming atmospheric heat was spent on ice melt.A simple calculation based on Equation (2) in Reference [15]: where γ = ρ i /ρ is the relation of mean ice and sea water densities that the temperature decrease due to ice melt by −0.4 • C, within the 23-m thick layer, led to a salinity decrease by −0.16 PSU, provided that the layer was well mixed.The latter was likely to be the case due to stormy conditions experienced from September 1-3, 2003 (see Figure 4a), which provided intensive turbulent mixing.Equation (1) for the UML depth h in Reference [15], gives that temperature and salinity decrease corresponds to 10 cm of ice melt (∆h i ).Taking into account that 6 cm of ice melt may be due to the atmospheric heat flux only (26 W/m 2 , the average heat balance value between September 1 and 5, see Figure 4a), this amount appears reasonable because of the additional bottom and side melt.Hence, we may conclude that under the icy conditions in 2003, by the end of the summer season, the heat storage in the upper layer remained small (7.4 MJ/m 2 ), since the incoming atmospheric heat was almost totally spent on the ice melt.However, the presence of open water does not guarantee a noticeable warming in the UML.In 2005, the recurrent CTD station was intermittently in the open water or in the MIZ 30 days prior to the first CTD cast.It is important to note that at the time of the first cast (September 14), the ice edge moved south reaching that point.Proximity to the ice edge probably determined the lack of warming in the UML in August and in the first half of September, despite the high positive heat flux (39 W/m 2 on average between August 15 and September 14).Another explanation is that the export of ice from the north in mid-September caused the rapid cooling of the UML due to the local ice melt.This option is less probable, because there was virtually neither a temperature change (−1.56 • C) in the UML, nor a substantial change in its thickness (−3 m) between the two casts.The salinity increase in the UML during this time interval (+0.09PSU) additionally points out that there was no ice melt.This slight salinity increase may have been attributed to an entrainment of cold halocline waters in the UML, which did not change the temperature.Heat flux at the ocean surface between the two recurrent CTD casts was small (11 W/m 2 on the average, see Figure 4b), thereby resulting in the preservation of nearly the same heat content in the UML (7.4 MJ/m 2 ) as during "icy" 2003, despite the relatively long duration of open water.
Build-up of the UML heat content in the absence of sea ice could be revealed from the two recurrent CTD casts in 2013 and 2015.Vertical temperature profiles on August 27, 2013, essentially rendered the profiles in icy 2003 and 2005.Down to the 40 m depth, the temperature was quasi-uniform and close to the freezing point.Salinity contrast between the surface and the 40 m depth was less than 2 PSU.A sort of halocline between 8 and 22 m was split into several segments with small and variable vertical gradients.Hence, we could anticipate that the vertical thermohaline structure during the first CTD cast in 2013 represented an initial state of the upper ocean water column after the ice cover had just disappeared.The latter was in line with the satellite ice maps and cruise observations.During the next 11 days, a well-defined UML embedded by sharp pycnocline had formed as a result of the heat uptake from the atmosphere and the turbulent mixing [8,16].
Average heat flux at the surface between August 27 and September 7 (∆t = 11 days) was equal to 50 W/m 2 .The temperature rise (∆T) was calculated from the one-dimensional heat balance equation described in Reference [3]: applied to h = 11 m, the UML is equal to 1.28 • C, while the actual temperature increase was 1.04 • C (c p denotes the specific heat of sea water at a constant pressure).Additional warming in the UML could be attributed to the horizontal mixing with surrounding waters.As a result of the heat uptake, by September 9, the UML temperature reached a temperature of −0.26In 2015, the first CTD cast was done on September 2, when the warming season was coming to an end.Therefore, the vertical structure of the upper ocean layer represented the state after the peak of warming had already passed, under conditions when open water existed for about one month.Very sharp pycnocline, between 19 and 26 m, separated the warm and fresh UML from the cold and saline waters below.The UML temperature equaled 1.61 • C, corresponding to a 3.29 • C departure from the freezing point.After 17 days (on September 19), the pycnocline upper boundary deepened to 22 m, while the lower boundary remained at 26 m, thereby increasing the vertical gradients of the properties.Despite a weak positive average heat flux at the ocean surface (10 W/m 2 according to the ERA-Interim data), the temperature of the UML decreased by 0.15 • C, reaching 1.46 • C. As was mentioned before, in 2015, the ERA-Interim data substantially contrasted with the direct measurements, which showed a total negative heat balance due to the dominance of turbulent fluxes.The latter was in line with the observed temperature decrease in the UML.The heat content of the UML on September 19 equaled to 282.3 MJ/m 2 , which was four times larger than in 2013.

Discussion
In this study, we focused on the specific features of the ocean-air interaction in the Laptev Sea during summer, on the basis of recurrent measurements during four expeditions in 2003, 2005, 2013, and 2015.The transition to an increasingly seasonal sea ice regime after 2007 draws a distinct timeline between contrasting ice conditions before and after this year, thereby allowing the assessment of possible changes in the ocean-air interaction processes caused by this transition.
Our main goal was to assess the relative significance of various factors in the buildup of heat content of the upper ocean layer at the doorstep of the next cold season.In a general sense, the increased duration of open water in the summer provides the prerequisite conditions for heat absorption in the upper ocean layer.However, the positive radiation balance in the high Arctic culminates in the middle of the astronomical summer, while the peak of the open water area typically occurs in the first half of September.By that time, the shortwave radiation flux is declining and a turbulent heat exchange may dominate the total heat balance at the ocean surface.Thus, the rate of removal of the heat absorbed in the upper ocean layer during the fall primarily depends on the atmospheric circulation patterns and the relation of these patterns to weather conditions.Assuming an average rate of heat loss through the fall in the range 100-200 W/m 2 , we required less than a day for the UML temperature to reach the freezing point in 2003 and 2005, then from 3 to 6 days in 2013, and from 16 to 32 days in 2015.The latter meant the visible delay of the freezing onset in 2015.The 15 day delay of ice emergence (compared to 2005 and 2013) was captured by satellite imagery in 2015 (see Figure 1-inset).However, it is impossible to detect whether the first ice at this site was formed locally or was imported from the north.A crude estimation of the possible effects of a delayed freeze onset on the seasonal ice thickness may be applied using the empirical formula described in [17]: where h is the ice thickness, T − is the air temperature below zero, and τ is the time interval of calculation.Although this equation was derived for the fast ice in the Laptev Sea, it is applicable for the estimation of ice growth in the open ocean, provided that oceanic heat flux is absent.
Taking the average winter SAT in the Laptev Sea as equal to −30 • C [18], we obtain the thickness of the one-year ice after 180 days of ice growth, equal to 182 cm.Reduction of the time interval by 16-32 days gave an 8-17 cm reduction in the one-year ice thickness.Although this value is rather small, it is detectable by direct measurements.Moreover, at the shorter time intervals, the relative thinning of ice is substantially greater, which may have an effect on the ice cover properties in mid-winter, as discussed below.
The obtained results may be cautiously extended on a larger scale.Taking into account the non-uniform spatial distribution of sea ice cover at the peak of summer, at a minimum we consider two specific segments around the Arctic.We define the 'Eurasian' segment within the bounds 60-120 • E, and the 'American' segment within the bounds 150 • E-150 • W. The southern and northern boundaries of both segments were defined at 65 • N and 80 • N, respectively.The western margin of the Eurasian segment was chosen as 60 • E to exclude the Western Nansen Basin and the Barents Sea, being the regions directly affected by the warm Atlantic Water, where the substantial reduction of winter ice cover in the 2010s was well documented [3,5].Longitudinal margins of the American segment were selected in order to cover the area where the record retreat of ice cover had occurred in 2007.
Time series of the mean anomalies of ice area and ice extent, that is, the border of 15% ice concentration (shown as the fraction of the total area of the corresponding segment) during the August-September-October (ASO) period and the December-January-February (DJF) period, for selected segments, are shown in Figure 5. Anomalies were calculated relative to the mean values during the time interval 1980-2000 when the ice conditions were rather stable.The mean monthly values of both parameters for the time intervals 1980-2000 and 2007-2017 are also shown in the same graphs.A similar feature in both time series, for both segments, was the substantial increase of the ASO negative anomalies during the 2000-2010s.In the American segment, the shift to the milder ice regime started in 2003, while in the Eurasian segment, this began two years later.The maximal area of open water in the American segment was observed in 2007.In the Eurasian segment, the maximal ice decrease was observed in 2012.Comparing Figure 5 with Figure 1 in Reference [19], we argue that from 2012 onwards, the interannual variability of the ice extent in the Eurasian segment (probably together with the Barents Sea and the Western Nansen Basin, not considered here) notably shapes the pan-Arctic pattern of the interannual variability.
The DJF time series demonstrated a marked shift to the milder ice regime after 2005 in the Eurasian segment.In the American segment, there is no change from 2003 onwards compared to the span from 1980-2000.This was consistent with the pattern at the mean seasonal variability plots.The shift to the reduced values of both ice parameters in the winter months was very well pronounced in the Eurasian segment.In the American segment, there was no such shift at all-from January until April, the ice area and extent were the same in the two considered time intervals, despite the fact that the summer ice reduction in the Eurasian sector was only half that of the American segment.Hence, we could conclude that the "seasonal memory" between the ice cover properties in summer and in the following winter after 2005 was typical only for the Eurasian segment.
To understand the possible reasons behind the "Eurasian-American" difference, we considered the time series of the 3-month average (ASO) atmospheric circulation patterns over the Arctic, derived from the ERA-Interim reanalysis daily data.These patterns were presented by the spatial and temporal components of the anomaly of sea level pressure (SLP), which were determined using the standard algorithm for empirical orthogonal function decomposition (EOF).EOFs are used to distinguish the main patterns (for example, the SLP) and to indicate the prevailing localization of structures [20,21].The EOFs are mathematically orthogonal and independent of each other.The first three EOFs describe 86% of the total variance (see Figure 6, left panels) and they are well separated.The DJF time series demonstrated a marked shift to the milder ice regime after 2005 in the Eurasian segment.In the American segment, there is no change from 2003 onwards compared to the span from 1980-2000.This was consistent with the pattern at the mean seasonal variability plots.The shift to the reduced values of both ice parameters in the winter months was very well pronounced in the Eurasian segment.In the American segment, there was no such shift at all-from January until April, the ice area and extent were the same in the two considered time intervals, despite the fact that the summer ice reduction in the Eurasian sector was only half that of the American segment.Hence, we could conclude that the "seasonal memory" between the ice cover properties in summer and in the following winter after 2005 was typical only for the Eurasian segment.
To understand the possible reasons behind the "Eurasian-American" difference, we considered the time series of the 3-month average (ASO) atmospheric circulation patterns over the Arctic, derived from the ERA-Interim reanalysis daily data.These patterns were presented by the spatial and temporal components of the anomaly of sea level pressure (SLP), which were determined using the standard algorithm for empirical orthogonal function decomposition (EOF).EOFs are used to distinguish the main patterns (for example, the SLP) and to indicate the prevailing localization of structures [20,21].The EOFs are mathematically orthogonal and independent of each other.The first three EOFs describe 86% of the total variance (see Figure 6, left panels) and they are well separated.
The first EOF pattern (58% of the total variance) represents Arctic oscillation (AO) [22,23].The positive phase of the AO (corresponded to EOF1) is associated with a negative SLP anomaly close to the North Pole.Almost the entire Arctic Ocean is covered by this circulation cell.However, Ogi et al.
[ 24,25] emphasize that the summer AO pattern has a smaller meridional extension than the conventional AO.The second EOF pattern (17% of the total variance) forms a dipole containing two cells that are opposite in signs.The first zone is centered over the southern Kara Sea and it is a more The first EOF pattern (58% of the total variance) represents Arctic oscillation (AO) [22,23].The positive phase of the AO (corresponded to EOF1) is associated with a negative SLP anomaly close to the North Pole.Almost the entire Arctic Ocean is covered by this circulation cell.However, Ogi et al. [24,25] emphasize that the summer AO pattern has a smaller meridional extension than the conventional AO.The second EOF pattern (17% of the total variance) forms a dipole containing two cells that are opposite in signs.The first zone is centered over the southern Kara Sea and it is a more compact cell, while the second zone is centered over the Canadian Basin.The former is stretched in the zonal direction, with the western periphery reaching the Norwegian Sea and the eastern periphery on the Laptev Sea shelf.The second EOF was previously detected as the Pacific/North American teleconnection pattern (PNA) [26].The third EOF (11% of the total variance) is more meridional, has a highly asymmetric dipole field, and corresponds to the Arctic dipole (AD) [27][28][29].A positive phase of the AD pattern is associated with a negative SLP anomaly on the North American side of the Arctic and a positive SLP anomaly on the Siberian side of the Arctic.However, the EOF's size of area and boundaries can influence the EOF patterns [30].EOF calculations for various sizes of areas were carried out in this study.It was shown that with a shortening of the zone, the key features of the first two EOFs remained the same (with a slight shift of maximum to the northeast); the third EOF changed significantly, where it formed a saddle pattern, the so called Buell patterns [30].Time series of the corresponding principal components (PCs) are shown on the right panels of Figure 6.PC2 attracts particular attention due to the presence of a steady increasing trend after the absolute minimum in 2007.Two other significant PCs did not show any sort of persistent behavior at the decadal time intervals.The range of PC1 interannual oscillations was the largest during the 1980s-1990s (over four units), while after 2000 it reduced by two times.PC3 also did not show any notable trends.Similar to PC1, the range of PC3's interannual variability substantially decreased at the end of the record (i.e., from 2010 onwards).
Our interpretation of these features in relation to the sea ice variability described above is as follows.Large scale cyclonic circulation over the Arctic Ocean (PC1) enforces a poleward drift of the ice in the American segment, but it slows down the ice export from the Siberian seas.At the same time, this cyclonic circulation maintains an inflow of air from the continent.The temperature of the continental air is warmer than the ice or water surface in the summer.In the fall, the temperature of the air, flowing in from the cooling continent, drops down, thereby favoring heat loss at the open water surface due to a turbulent heat exchange.In the Eurasian segment, a positive PC1 retards summer warming due to an inflow of the cold air from the ice-covered ocean north of Greenland.By the end of September-early October, this cooling effect might become weaker, because by that time, the sea surface in the northern Barents Sea has already been heated by the solar radiation and it warms up the incoming air from the central basin.Therefore, the anomalous reduction of the summer sea ice area and extent in the American segment, as well as the fast cooling at the ocean surface in the fall, may be explained (at least partly) by the large scale cyclonic circulation over the Arctic Ocean.
Atmosphere 2018, 9, x FOR PEER REVIEW 13 of 23 areas were carried out in this study.It was shown that with a shortening of the zone, the key features of the first two EOFs remained the same (with a slight shift of maximum to the northeast); the third EOF changed significantly, where it formed a saddle pattern, the so called Buell patterns [30].Time series of the corresponding principal components (PCs) are shown on the right panels of Figure 6.
PC2 attracts particular attention due to the presence of a steady increasing trend after the absolute minimum in 2007.Two other significant PCs did not show any sort of persistent behavior at the decadal time intervals.The range of PC1 interannual oscillations was the largest during the 1980s-1990s (over four units), while after 2000 it reduced by two times.PC3 also did not show any notable trends.Similar to PC1, the range of PC3's interannual variability substantially decreased at the end of the record (i.e., from 2010 onwards).
Our interpretation of these features in relation to the sea ice variability described above is as However, PC1 does not provide any clue to the explanation of the revealed "seasonal memory" between the sea ice properties during the ASO and DJF periods in the Eurasian segment.We argue that an increased contribution of PC2 after 2005, which turned from negative to highly positive values in recent years, was the major driver behind this "seasonal memory".An elongated shape of the EOF2 anticyclonic cell ensured the pumping of warm and humid air from the ice-free Norwegian Sea and southern Barents Sea when the PC2 was positive.In the fall, this relatively warm air served as an "insulator" for the heat accumulated in the upper ocean during the summer, slowing down the UML cooling and delaying the onset of freezing.This insulating effect was obviously the strongest in the western part of the Eurasian segment.For example, in the fall of 2016, the anomalously high sea surface temperatures over the Barents-Kara Seas created an extremely low pan-Arctic sea ice extent throughout the fall and the early winter.In mid-November 2016, the sea ice extent even decreased for several days [31].In the American segment, the effect of an enhanced PC2 was the opposite.
Counterclockwise circulation over the Canadian Basin delivered cold air from the Greenland pack ice towards Alaska and the Bering Strait, leading to enhanced heat loss from the ocean surface to the atmosphere, facilitating faster ice formation.

Conclusion
The progressive reduction of Arctic sea ice cover alters the processes of ocean-air interaction.In the Arctic, ocean-air turbulent heat exchange is particularly significant in the late summer when the short wave radiation balance at the ocean surface is reduced, and the turbulent heat fluxes take over in the total heat balance.Local ice formation in the open water commences only after the sea surface temperature drops to the freezing point.The required time for the freezing onset depends on the amount of heat accumulated in the UML throughout the summer, as well as the rate of this heat removal.Joint analysis of the observations, atmospheric reanalysis products, and satellite data allowed for the assessment of the heat storage changes in the UML related to the variable atmospheric forcing under contrasting ice conditions in the Laptev Sea at the doorstep of the fall.It was established that in the "icy" years (2003 and 2005), the accumulation of heat in the UML was insignificant for the subsequent ice formation because the absorbed heat was mostly spent on the ice melt until the end of the summer.In the "ice-free" years (2013 and 2015), the ultimate heat storage in the UML depended on the duration of open water, the distance of the point of interest to the nearest ice edge, and on the synoptic weather patterns.We did not find any direct links between the UML heat content and the strength of the underlying pycnocline, although such links could be expected.In our opinion, the water column in the Laptev Sea has so far not undergone a significant structural transformation, However, PC1 does not provide any clue to the explanation of the revealed "seasonal memory" between the sea ice properties during the ASO and DJF periods in the Eurasian segment.We argue that an increased contribution of PC2 after 2005, which turned from negative to highly positive values in recent years, was the major driver behind this "seasonal memory".An elongated shape of the EOF2 anticyclonic cell ensured the pumping of warm and humid air from the ice-free Norwegian Sea and southern Barents Sea when the PC2 was positive.In the fall, this relatively warm air served as an "insulator" for the heat accumulated in the upper ocean during the summer, slowing down the UML cooling and delaying the onset of freezing.This insulating effect was obviously the strongest in the western part of the Eurasian segment.For example, in the fall of 2016, the anomalously high sea surface temperatures over the Barents-Kara Seas created an extremely low pan-Arctic sea ice extent throughout the fall and the early winter.In mid-November 2016, the sea ice extent even decreased for several days [31].In the American segment, the effect of an enhanced PC2 was the opposite.Counterclockwise circulation over the Canadian Basin delivered cold air from the Greenland pack ice towards Alaska and the Bering Strait, leading to enhanced heat loss from the ocean surface to the atmosphere, facilitating faster ice formation.

Conclusions
The progressive reduction of Arctic sea ice cover alters the processes of ocean-air interaction.In the Arctic, ocean-air turbulent heat exchange is particularly significant in the late summer when the short wave radiation balance at the ocean surface is reduced, and the turbulent heat fluxes take over in the total heat balance.Local ice formation in the open water commences only after the sea surface temperature drops to the freezing point.The required time for the freezing onset depends on the amount of heat accumulated in the UML throughout the summer, as well as the rate of this heat removal.Joint analysis of the observations, atmospheric reanalysis products, and satellite data allowed for the assessment of the heat storage changes in the UML related to the variable atmospheric forcing under contrasting ice conditions in the Laptev Sea at the doorstep of the fall.It was established that in the "icy" years (2003 and 2005), the accumulation of heat in the UML was insignificant for the subsequent ice formation because the absorbed heat was mostly spent on the ice melt until the end of the summer.In the "ice-free" years (2013 and 2015), the ultimate heat storage in the UML depended on the duration of open water, the distance of the point of interest to the nearest ice edge, and on the synoptic weather patterns.We did not find any direct links between the UML heat content and the strength of the underlying pycnocline, although such links could be expected.In our opinion, the water column in the Laptev Sea has so far not undergone a significant structural transformation, similar, for instance, to the Western Nansen Basin [32].Hence, the heat accumulation in summer was confined to the relatively shallow (up to 20-25m) surface layer.However, the heat capacity of this thin surface layer may be large enough to delay the onset of freezing by up to tens of days, depending on the rate of surface cooling during the fall.
In the format of discussion, we considered the possible links between the average ice area/extent in August-September-October (ASO) period, and in the December-January-February (DJF) period, for two representative Arctic regions, that is, the Eurasian segment, defined within the bounds 60-120 • E, 65-80 • N, and the American segment, defined within the bounds 150 • E-150 • W, 65-80 • N. Significant "seasonal memory", characterized by the consistent change of the ice cover parameters in sequential seasons, was retrieved in the Eurasian segment between 2007 and 2017.No linkage to the seasonal time scale was found for the American segment, despite an even stronger reduction of the ASO ice extent/area in this region from 2007 onwards.We explain the distinguished "Eurasian-American" contrast as being the result of the increasing (since 2007) intensity of a dipole circulation mode over the Arctic Ocean in the ASO period.This circulation mode is described by the second-in-order member of the Arctic SLP EOF-decomposition, containing an elongated anticyclonic cell with a longer axis, stretched eastwards over the northern margin of Eurasia, and the conjugated round-shaped cyclonic cell over the Canadian Basin.Zonal flow at the northern periphery of the anticyclonic circulation carries warm and humid air from the Nordic Seas across the Siberian shelf seas up to the western margin of the East-Siberian Sea, thereby slowing down the surface water cooling in the fall.In the American segment, on the contrary, an enhanced cyclonic circulation delivers cold air from the Greenland ice pack, thereby accelerating ocean-air heat loss and reducing the time interval before the onset of freezing.Since the presented observational data analysis was carried out for a specific region of the Arctic Ocean, our conclusions at the pan-Arctic scale should be considered as tentative, requiring additional independent testing.where ϕ q (ξ) is a universal function that describes the vertical profile of the humidity.The form of this function depends on the vertical stratification.In case of an unstable stratification, the function form Reference [41] was used, whilst for a stable stratification, the function form was taken from Reference [42].

Figure 1 .
Figure 1.Mean sea ice extent (>15% ice concentration) in September 2003, 2005, 2013, and 2015, represented by different colors.Position of the recurrent conductivity temperature depth (CTD) station is shown by the corresponding color squares.Bottom bathymetry is shown by thin black lines.Inset: the distance (km) from the position of the recurrent CTD station to the nearest ice edge.Bathymetry contours 300, 500, 1000, 2000, and 3000 m are shown by gray lines.

Atmosphere 2018, 9 ,
x FOR PEER REVIEW 5 of 23 by the temperature and salinity profiles in the upper 40 m depth layer at the recurrent CTD stations in Figure2.In 2003 and in 2005, water temperature in the entire 40-m layer was negative and did not exceed −1 °C.In 2003, the maximal departure from the freezing point, up to 0.5 °C , was observed at the base of the UML.In 2005, the maximal departure was less than 0.25 °C.On the contrary, in 2013 and 2015, the departure from the freezing point in the UML exceeded 1 degree during the second CTD cast, reaching 3 degrees in 2015.However, on the 5-14 day time scale, which corresponds to the time interval between recurrent CTD casts, the effect of ice conditions on the vertical structure was less unequivocal.For example, the evolution of the UML demonstrated definite similarities in specific years with contrasting ice conditions.In "icy" 2003 and in "ice-free" 2015, the thickness and salinity of the UML increased unidirectionally between the two recurrent CTD casts, reaching about the same value of 23 m during the second cast.In 2005, the UML was 4 m shallower during the second cast.In 2013, there was no distinct seasonal UML during the first cast at all, while at the second cast, it appeared within the 11 m depth.Temperature change in the UML between the two casts in the neighboring "ice-free" years was opposite-an increase in 2013 and a decrease in 2015, while corresponding salinity changes were unidirectional (increase) in both years.These similarities and differences in the observed short-term changes over the years, which were characterized by different and similar ice conditions, point to variable atmospheric forcing as a probable factor affecting the state of the upper ocean layer in the warm season.

Figure 2 .
Figure 2. Temperature and salinity vertical profiles at the recurrent CTD stations in different years.The first CTD cast is shown by the blue line.The second cast is shown by the red line.Freezing point temperature profile is shown by the black line.

Figure 2 .
Figure 2. Temperature and salinity vertical profiles at the recurrent CTD stations in different years.The first CTD cast is shown by the blue line.The second cast is shown by the red line.Freezing point temperature profile is shown by the black line.

Atmosphere 2018, 9 ,
x FOR PEER REVIEW 6 of 23 enter the affected area.This feature is a principal difference between the Polar regions and the lower 189 latitudes.In the mid-latitudes, there is a distinct division between the summer season, when the 190 colder ocean gains heat from the warmer air, and the fall-winter season, when the warmer ocean 191 releases heat to the atmosphere.In the Polar regions, even in the mid-summer season, the upper 192 ocean layer may act as a heat absorber, as well as a heat emitter, because of the near-zero sea surface 193 temperature (SST) close to the ice edge under the ice-free conditions[14].An air flow, which floods 194 over an ice-free surface, may be either warmer, if originated over the continent, or colder, when 195 coming from the ice covered regions in the high-latitude Arctic.196 Considering these possible options, we assessed the synoptic conditions over the four 197 considered years according to the ERA-Interim reanalysis data.In 2003, the position of the recurrent 198 CTD station was permanently at the periphery of the vast cyclone occupying almost the entire Arctic 199 Ocean (Figure 3a).Western winds forced the ice edge to the east, keeping the point in the MIZ.In 200 2005, between September 14 and 20, the central Laptev Sea was located in a saddle-shaped sea level 201 pressure (SLP) field between four large scale pressure systems, centered over the Arctic Ocean and 202 the adjacent continent (Figure 3b).Weak SLP gradients caused intermittent winds and the slow 203 transport of air masses.Within the entire time interval between the two recurrent CTD stations in 204 2013, the study area remained in the pathway of warm air flow from the continent between the 205 anticyclone over the central Arctic Ocean, and the cyclone, centered in the Kara Sea (Figure 3c).The 206 time interval between the recurrent CTD stations in 2015 was the longest compared to the other years, 207 that is, two and a half weeks.The study area was mainly in the cyclonic circulation field.This caused 208 the prevalence of western winds from the Kara Sea and an inflow of the air with near-zero °C surface 209 air temperature (SAT) (Figure 3d).

Figure 3 .
Figure 3. Horizontal distribution of mean sea level pressure (SLP) (solid lines) and mean surface air temperature (SAT) (color) during the time interval between the recurrent CTD casts in 2003 (a), 2005 (b), 2013 (c), and 2015 (d) according to the ERA-Interim reanalysis data.Frequency distribution of the wind direction and speed are presented in the polar rose plots.Position of the recurrent CTD station is shown by a yellow circle.

Figure 4 .
Figure 4. Time series of the heat balance components (color bars), total heat balance (solid lines), and wind speed (graph) at the recurrent CTD stations in 2003 (a), 2005 (b), 2013 (c), and 2015 (d), on the basis of the ERA-Interim reanalysis data.Black dotted lines indicate the timing of the recurrent CTD casts.Blue dotted lines indicate the time of the open water emergence according to satellite observations (small dots) and according to the ERA-Interim data (big dots).Positive heat flux values indicate radiative and turbulent heat transfer to the ocean.

Figure 4 .
Figure 4. Time series of the heat balance components (color bars), total heat balance (solid lines), and wind speed (graph) at the recurrent CTD stations in 2003 (a), 2005 (b), 2013 (c), and 2015 (d), on the basis of the ERA-Interim reanalysis data.Black dotted lines indicate the timing of the recurrent CTD casts.Blue dotted lines indicate the time of the open water emergence according to satellite observations (small dots) and according to the ERA-Interim data (big dots).Positive heat flux values indicate radiative and turbulent heat transfer to the ocean.
for all the discussed casts) was −1.16 • C. Five days later, it decreased by the ∆T = 0.4 • C, almost reaching the freezing point, and the UML deepened to 23 m.The mean salinity decrease (∆S) in the upper 19-20 m was equal to 0.15 PSU, from Atmosphere 2019, 10, 184 9 of 23 graphs.A similar feature in both time series, for both segments, was the substantial increase of the 386 ASO negative anomalies during the 2000-2010s.In the American segment, the shift to the milder ice 387 regime started in 2003, while in the Eurasian segment, this began two years later.The maximal area 388 of open water in the American segment was observed in 2007.In the Eurasian segment, the maximal 389 ice decrease was observed in 2012.Comparing Figure 5 with Figure 1 in Reference [19], we argue that 390 from 2012 onwards, the interannual variability of the ice extent in the Eurasian segment (probably 391 together with the Barents Sea and the Western Nansen Basin, not considered here) notably shapes 392 the pan-Arctic pattern of the interannual variability.

Figure 5 .
Figure 5.Time series of the sea ice extent anomalies (light blue bars: December, January, February (DJF); pink bars: August, September, October (ASO)) and sea ice area anomalies (blue line: DJF, red line: ASO) in the left panel, and the seasonal variation of sea ice extent (light blue bars: 1980-2000, pink bars: 2007-2017) and sea ice area (blue line: 1980-2000, red line: 2007-2017) in the right panel, for the Eurasian segment of the Arctic Ocean, 60-120° E, 65-80° N (a); the same for the American segment of the Arctic Ocean, 150° E-150° W, 65-80° N (b) according to the DMSP (Defense Meteorological Satellite Platform) digitized passive microwave satellite images [10].Anomalies are calculated as decimal fractions of the mean values over the 1980-2000 period.

Figure 5 .
Figure 5.Time series of the sea ice extent anomalies (light blue bars: December, January, February (DJF); pink bars: August, September, October (ASO)) and sea ice area anomalies (blue line: DJF, red line: ASO) in the left panel, and the seasonal variation of sea ice extent (light blue bars: 1980-2000, pink bars: 2007-2017) and sea ice area (blue line: 1980-2000, red line: 2007-2017) in the right panel, for the Eurasian segment of the Arctic Ocean, 60-120 • E, 65-80 • N (a); the same for the American segment of the Arctic Ocean, 150 • E-150 • W, 65-80 • N (b) according to the DMSP (Defense Meteorological Satellite Platform) digitized passive microwave satellite images [10].Anomalies are calculated as decimal fractions of the mean values over the 1980-2000 period.

Figure 6 .
Figure 6.Horizontal distribution of the first empirical orthogonal function (EOF1) for the mean August-September-October sea level pressure (left panel), and the time series of the corresponding principal component PC1 (right panel) (a); (b,c) the same for the second and the third EOFs.

1 ) 23 Table A1 .Figure A1 .
Figure A1.Ice conditions along the ship track from August 26-29, 2013, based on the visual observations from the bridge of the RV "Akademik Fedorov".

Figure A1 .
Figure A1.Ice conditions along the ship track from August 26-29, 2013, based on the visual observations from the bridge of the RV "Akademik Fedorov".

A. 3 .
Figures A2-A4 shows the time series of the measured onboard meteorological parameters in 2005, 2013, and 2015.Figure A5 presents a comparison between the daily-averaged heat balance components according to direct measurements at the ocean-air interface and the interpolated ERA-Interim reanalysis data.

Figures
Figures A2, A3 and A4 shows the time series of the measured onboard meteorological parameters in 2005, 2013, and 2015.Figure A5 presents a comparison between the daily-averaged heat balance components according to direct measurements at the ocean-air interface and the interpolated ERA-Interim reanalysis data.

Figure A2 .
Figure A2.Time series of the measured onboard meteorological parameters in 2005.Positive values of the heat balance components define fluxes from the atmosphere to the ocean for the radiation fluxes (Qs and Ql), and vice versa for the turbulent fluxes (H and LE).

Figure A2 .
Figure A2.Time series of the measured onboard meteorological parameters in 2005.Positive values of the heat balance components define fluxes from the atmosphere to the ocean for the radiation fluxes (Q s and Q l ), and vice versa for the turbulent fluxes (H and LE).

Figure A3 .
Figure A3.Time series of the measured onboard meteorological parameters in 2013.Designations are similar to Figure A2.

Figure A3 .
Figure A3.Time series of the measured onboard meteorological parameters in 2013.Designations are similar to Figure A2.

Figure A4 .
Figure A4.Time series of the measured onboard meteorological parameters in 2015.Designations are similar to Figure A2.

Figure A4 .
Figure A4.Time series of the measured onboard meteorological parameters in 2015.Designations are similar to Figure A2.

Figure A5 .
Figure A5.Time series of the heat balance components (color bars), total heat balance (solid lines), and the wind speed (graph) along the ship route on the basis of direct measurements (left panels) and on the basis of the ERA-Interim reanalysis data (right panels) in 2005 (a,b), 2013 (c,d), and 2015 (e,f).Black dotted lines indicate the recurrent CTD casts.Positive heat flux values indicate radiative and turbulent heat transfer to the ocean.

Figure A5 .
Figure A5.Time series of the heat balance components (color bars), total heat balance (solid lines), and the wind speed (graph) along the ship route on the basis of direct measurements (left panels) and on the basis of the ERA-Interim reanalysis data (right panels) in 2005 (a,b), 2013 (c,d), and 2015 (e,f).Black dotted lines indicate the recurrent CTD casts.Positive heat flux values indicate radiative and turbulent heat transfer to the ocean.

.
Description of the recurrent conductivity temperature depth (CTD) stations.The dates of 154 opening or closing signify the first date when the CTD station was outside or inside the border of 15% 155 ice concentration.The mean distance to the ice edge is the nearest linear path to the border of 15% ice

Table 1 .
Description of the recurrent conductivity temperature depth (CTD) stations.The dates of opening or closing signify the first date when the CTD station was outside or inside the border of 15% ice concentration.The mean distance to the ice edge is the nearest linear path to the border of 15% ice concentration in any direction.

Table 2 .
Mean values of the observed surface air temperature (SAT), sea surface temperature (SST), and their difference ∆T, averaged over the whole cruise track in the Laptev sea region (100-150 • E) for the 2013 and 2015 cruises.
• C (1.46 • C above the freezing point).The heat content by this date reached 65.8 MJ/m 2 , which was an order of magnitude higher than in 2003 and in 2005.

Table A1 .
The list of instruments used and the measured parameters during the 2013 and 2015 cruises.