The Summer Surface Energy Budget of the Ice-Free Area of Northern James Ross Island and Its Impact on the Ground Thermal Regime

: Despite the key role of the surface energy budget in the global climate system, such investigations are rare in Antarctica. In this study, the surface energy budget measurements from the largest ice-free area on northern James Ross Island, in Antarctica, were obtained. The components of net radiation were measured by a net radiometer, while sensible heat ﬂux was measured by a sonic anemometer and ground heat ﬂux by heat ﬂux plates. The surface energy budget was compared with the rest of the Antarctic Peninsula Region and selected places in the Arctic and the impact of surface energy budget components on the ground thermal regime was examined. Mean net radiation on James Ross Island during January–March 2018 reached 102.5 W m − 2 . The main surface energy budget component was the latent heat ﬂux, while the sensible heat ﬂux values were only 0.4 W m − 2 lower. Mean ground heat ﬂux was only 0.4 Wm -2 , however, it was negative in 47% of January–March 2018, while it was positive in the rest of the time. The ground thermal regime was a ﬀ ected by surface energy budget components to a depth of 50 cm. The strongest relationship was found between ground heat ﬂux and ground surface temperature. Further analysis conﬁrmed that active layer refroze after a sequence of three days with negative ground heat ﬂux even in summer months. Daily mean net radiation and ground heat ﬂux were signiﬁcantly reduced when cloud amount increased, while the inﬂuence of snow cover on ground surface temperature was negligible.


Introduction
Polar regions with their vast ice-sheets are important climatic factors influencing the climate of the whole planet. The white surface of glaciers in Antarctica and in the Arctic has an albedo of up to 90% [1], thereby reflecting incoming shortwave radiation and reducing the amount of energy available to the Earth. In the Antarctic, the net amount of incoming solar radiation further affects the grade of turbulent fluxes, atmospheric circulation on a large scale and, as a consequence, the energy budget of the whole planet.
More than 99% of Antarctica is covered by ice [2], but the ice-free area is further predicted to enlarge [3][4][5], especially in the Antarctic Peninsula (AP) region. Despite the fact that warming in the AP region has decelerated or even been absent since the 2000s, large decadal [6] and spatial variability in air temperature [7] was found to be a typical feature in this area. Therefore, it does not necessarily indicate a return to a colder climate from a long-term perspective [8].
Most studies have analysed the surface energy budget over sea-ice, snow and ice-covered areas (e.g., [9][10][11][12]). Ice-free areas, such as in the AP region, on the other hand, have received less attention. Choi et al. [13] reported that the mean monthly net radiation on an ice-free area on King George Island 1.
To quantify the summer surface energy budget over the ice-free area on the Ulu Peninsula; 2.
To identify the connection between components of the surface energy budget and verify if the energy fluxes differ between the western and eastern side of the AP; and 3.
To determine the effect of snow cover and cloudiness on the surface energy budget and ground thermal regime.

Study Site and Instrumentation
The Ulu Peninsula, covering 552 km 2 (see [26]), is the largest ice-free area in the AP region. JRI lies to the south-east of the tip of the AP (Figure 1a), with a prevailing easterly flow often modified by the AP orography in relation to the southern and south-western winds [27]. Consequently, the local climate on JRI is characterized by a 4-5 • C lower mean air temperature compared to the north-western AP [28]. The study site was situated at 10 meters above sea level (m a.s.l.) on a Holocene marine terrace in the northern part of the Ulu Peninsula ( Figure 1). The Johann Gregor Mendel Station is located approximately 100 m to the north-east of the meteorological tower and a sonic anemometer, while the seashore is about 150 m to the north. The area is composed of dry sandy ground (about 10% moisture) underlain by continuous permafrost with an active layer of thickness between 50 and 65 cm [23]. temperature from sonic anemometer were also included in the latent heat flux calculation. The ground surface temperature was defined as the ground temperature at a depth of 2 cm. and wind rose representing wind frequencies with respect to their wind speed from study site for period from 23 January to 21 March, 2018 (c). Table 1. Overview of meteorological parameters ("Parameter") measured on James Ross Island during 23 January-20 March, 2018. In the "Height/Depth" column, positive numbers refer to height above the ground, while negative numbers refer to depth.  The measurements were carried out between 23 January and 21 March, 2018 (58 days). Details on the instruments used for measurements are in Table 1 and they were all temporarily installed at one site. For the latent heat flux calculation [29], we used the measurements of 2-m air temperature, relative humidity and surface ground temperature measurements. As only relative humidity was measured, 2-m specific humidity was calculated with the use of air temperature and atmospheric pressure data. Furthermore, vertical kinematic eddy heat flux, Monin-Obukhov length and virtual temperature from sonic anemometer were also included in the latent heat flux calculation. The ground surface temperature was defined as the ground temperature at a depth of 2 cm.

Data Processing and Analysis
The surface energy budget of an ice-free surface can be written as: where Q net is the net radiation, Q G is the ground heat flux, Q H is the sensible heat flux, Q E is the latent heat flux and C is the residual. The Q net has the following components: where Q Snet is the net shortwave flux consisting of incoming (Q S↓ ) and outgoing (Q S↑ ) shortwave radiation, while Q Lnet is the net longwave flux, which can be divided into incoming (Q L↓ ) and outgoing (Q L↑ ) longwave radiation. All of the measured data were processed as 30-min averages. The Q S↓ measurements were further corrected with the use of a more precise CM11 pyranometer (Kipp and Zonen). Simple linear regression was used for adding the relative deviation to CNR4 measurement [30]. Surface albedo was calculated daily as the ratio of Q S↑ and Q S↓ in the 3 h closest to the maximum of Q S↓ in order to avoid measurement unreliability under large solar zenith angles and low flux magnitudes.
Ground heat flux measurement was corrected using the calorimetric method [31,32] as follows: where G i-1 is the heat flux density at the top of a layer, here at a depth of 2 cm, δz i is the layer thickness (3 cm), C i is the volumetric heat capacity for the layer (1.17 MJ m −3 K −1 defined for the area by Hrbáček et al. [23]), ∂T i /∂t is the rate of a 30-min change in the mean layer temperature and G i is the heat flux density at the bottom of the layer at a depth of 5 cm. The sensible heat flux was calculated from sonic anemometer using the eddy covariance technique (e.g., [33]). However, it needs to be noted that when the term "sensible heat flux" is used in this study, it is the buoyancy flux, which was calculated instead of true sensible heat flux due to the absence of fast-response moisture measurements [34]. Sensible heat fluxes of 30-min average were calculated only if Atmosphere 2020, 11, 877 5 of 18 more than 54% of the measurement were available. Consequently, <2% of sensible heat flux data were missing. Latent heat was calculated as in Langer et al. [29]. Briefly, latent heat flux was estimated as: where ρ air is the density of air, L lg is the latent heat of vaporization of water, r a is atmospheric resistance, q is specific humidity, z m is measurement height and z 0 is roughness length. Roughness length was originally chosen to be 10 −3 m over bare ground and 10 −4 m over a snow-covered surface [35,36], but due to unrealistically high values of latent heat flux, they were both lowered by one order to 10 −4 and 10 −5 m, respectively. The atmospheric resistance r a was calculated as in formula D4 in Langer et al. 2011 [29]. Specific humidity at the surface q(z 0 ) was derived from the surface temperature using the Magnus formula [37]. Latent heat flux measurements were only available for 92% of the study period. For all surface energy budget components, daily and weekly means were calculated if <20% were missing, hence the loss of daily mean fluxes on 8 February and 12-14 March, 2018 ( Figure 2).  The data were mostly analysed and compared as daily means; 30-min averages were only examined in Table A1 and day-long case studies in order to determine the influence of snow and cloudiness on the surface energy budget (Section 3.4). For the correlation analysis (Section 3.3), Spearman's rank correlation coefficient was used due to non-normal distributions of some data. Simple linear regression was employed for calculating regression function in Section 3.3. The  The data were mostly analysed and compared as daily means; 30-min averages were only examined in Table A1 and day-long case studies in order to determine the influence of snow and cloudiness on the surface energy budget (Section 3.4). For the correlation analysis (Section 3.3), Spearman's rank correlation coefficient was used due to non-normal distributions of some data. Simple linear regression was employed for calculating regression function in Section 3.3. The following convention was used: the net radiation components were considered positive if energy transport was directed towards the ground surface; the turbulent and ground heat fluxes were denoted as positive if energy transport was directed away from the ground surface.
For the case studies, it was necessary to determine a cloudy and a clear-sky day. We only had three categories: clear sky, cloudy and overcast. Cloudiness was derived from a combination of hemispherical photos and Q S↓ measurements. We calculated the mean daily ratio of Q S↓ to extra-terrestrial radiation (simply approximated from solar constant and solar zenith angle) in hourly intervals. We then estimated that if the daily mean ratio was above 0.6, the sky was (mostly) clear. If the ratio was below 0.4, the day was overcast. Ratios between 0.4 and 0.6 meant a cloudy day.

Measurement Accuracy
The residual term C in our study was 36.1 ± 104.0 W m −2 (Table A1), indicating that errors prevented surface energy budget closure. Unlike the measurements from KGI, where C was not determined but was contained within Q H [17] or Q G [13], we tried to limit the uncertainty by calculating unmeasured Q E . Non-zero C has been almost unavoidable in surface energy budget studies from polar regions [29,38], even though errors should have been cancelled out unless there was a systematic error. In our case, 35% of Q net was allocated to neither Q G nor the turbulent fluxes due to a combination of multiple measurement errors.
Firstly, the measurement accuracy of the CNR4 for longwave radiation of ±10% might have increased the values of Q net (Table 2). Moreover, Westermann et al. [38] pointed out that uncertainty about albedo under conditions with patches of snow and bare-ground could be up to 0.1, which would further lower or raise Q net by 3.4%.
Secondly, the eddy covariance technique requires stationarity. However, complex orography around the measuring site was prone to the development of local circulation patterns [36]. That could have led to systematic underestimation of Q H up to 25% [39], while the use of buoyancy flux as Q H means that another 10% might need to be added to Q H [34].
The error in Q E was difficult to assess because no comparable measurements existed. Nevertheless, Q E calculating depended strongly on the correct choice of z 0 . In this study, using the values for snow-covered surfaces and bare ground according to current knowledge [35,36] could lead to an increase of Q E up to 33%. On the other hand, as both Prosek et al. [17] and Choi et al. [13] reported mean Q E~QH , it is likely that the error is smaller. Even though the maximum uncertainty in Q G estimation caused by instrument accuracy could have been up to ±20%, the mean value of Q G was below 1 W m −2 . As a result, the error in Q G would not contribute significantly to C.
Finally, individual energy budget components were measured on different spatial scales. Q G was derived from a point measurement and radiation sensors had a footprint area of several square meters. For the sonic anemometer at a height of 2 m, 60% of the flux might have come from an area with a radius of 75 m or larger [40]. During high wind speeds, the turbulent eddies (and derived energy fluxes) measured by the sonic anemometer were coming from the Bohemian Stream Catchment (see Figure 1 and Section 3.1), where the ground was waterlogged or moister than close to the meteorological tower. Consequently, more energy was allocated to Q E in the area from which the turbulence was coming, influencing the residual term C.
Considering the above-mentioned errors, the magnitude of C could be explained by the uncertainty inherited from methods of measurements, the landscape diversity around the measuring site and the choice of parameters in the calculation. Essential for future experiment concept and better accuracy Atmosphere 2020, 11, 877 7 of 18 of eddy-covariance method could be two-level humidity measurement or krypton hygrometer for measuring water vapour fluctuations with a high-frequency sampling rate.

Meteorological Conditions and Energy Budget Components
From 23 January to 21 March, 2018, the mean daily air temperature on the Ulu Peninsula was −2.5 • C with a range from −13.9 • C to 5.3 • C. Daily mean air temperatures often changed from negative to positive, however, they were <0 • C after 25 February, with the exception of 28 February and 11 March ( Figure 2d). On the other hand, the daily mean ground surface temperature was mostly positive before 28 February and mostly negative afterwards. The daily mean ground temperature at 50 cm was −0.2 • C, varying between −1.3 and 0.3 • C only. In 21% of the measurement period, the wind blew from south-southwest (Figure 1c), where the Bohemian Stream catchment lies. The highest mean daily wind speed of 13.6 ms −1 was measured on 15 March, 2018 during south-westerly wind (Figure 2e).
Mean albedo was 0.11 when the ground surface was bare, while it rose to 0.54 on the days with snow cover (Figure 2b). On 8 and 9 February, albedo increased to over 0.70. Snow cover also intensified after 13 March and during 2-5 March, when 3-cm deep snow cover led to the albedo reaching 0.64. The snow cover periods corresponded well to a decrease in mean daily Q net , which was mostly lower than 40 W m −2 during the intensified snow cover period and only reached 60 W m −2 on 3 March.
Mean Q net was 102.5 W m −2 and mean Q L↓ (264.4 W m −2 ) was by 45% higher than mean Q S↓ (182.7 W m −2 ). Mean Q L↑ was 311.1 Wm −2 , while mean Q S↑ was 33.4 Wm −2 (Table A1). The highest energy budget component, apart from the radiation fluxes, was Q E with a mean value of 33.2 W m −2 , closely followed by Q H (32.8 W m −2 ). In both cases, standard deviations were about 1.5times the mean value of the flux. The mean Q G was only 0.5 W m −2 due to the averaging-out of positive and negative Q G (Figure 2c). The residual term C reached 36.1 W m −2 as a mean value for the study period (more information in Table A1 and Figure A1).
Both mean Q H and Q E depended on wind direction ( Figure 3). Their mean values were the lowest when the wind was from the east and south-south-east (146.5 • to 168.5 • ), while they were the highest for north-easterly and west-north-westerly wind (281.5 • to 303.5 • ). When the wind was westerly, both mean Q H and Q E increased with increasing wind speed (data not shown). Mean Q H was by 11.1 Wm −2 higher than mean Q E for the south-south-westerly wind direction (191.5 • to 213.5 • ), while it was lower for easterly wind (Q H : 2.5 Wm −2 , Q E : 13.4 Wm −2 ) and for west-south-westerly to north-west-northerly wind (236.5 • to 348.5 • , the difference between mean Q H and mean Q E of up to 19.6 Wm −2 ). The 95% confidence intervals for both Q H and Q E were the lowest for south-south-westerly wind direction (191.5 • to 213.5 • ) and the highest for northerly and north-north-westerly wind direction for Q H and Q E , respectively. It is, therefore, clear that in cases of low wind velocity (wind from the northerly sector), higher uncertainty on the closure of surface energy budget occurred and caused higher variation of both turbulent fluxes in the related wind directions.

Ground Temperature Response to the Largest Surface Energy Fluxes
The ground thermal regime was significantly influenced by individual components of the surface energy budget. Mean daily QE was higher than QH in 63% of the study period ( Figure 2). The highest mean daily QH was 100.1 W m -2 on 1 February (Figure 2). Between 31 January and 1 February, the daily mean ground surface temperature increased by 1.7 °C. Mean daily QE was the highest on 24 February (116.2 W m -2 ), when QH reached only 50.3 Wm -2 . A concurrent increase in air and ground surface temperature before 24 February pointed to active layer thawing ( Figure 4). Between 24 and 25 February, the ground surface temperature decreased by -5.8 °C, indicating that intense evaporation might have led to energy depletion in ground surface layers, which resulted in the drop in surface temperature.

Ground Temperature Response to the Largest Surface Energy Fluxes
The ground thermal regime was significantly influenced by individual components of the surface energy budget. Mean daily Q E was higher than Q H in 63% of the study period ( Figure 2). The highest mean daily Q H was 100.1 W m −2 on 1 February (Figure 2). Between 31 January and 1 February, the daily mean ground surface temperature increased by 1.

Ground Temperature Response to the Largest Surface Energy Fluxes
The ground thermal regime was significantly influenced by individual components of the surface energy budget. Mean daily QE was higher than QH in 63% of the study period ( Figure 2). The highest mean daily QH was 100.1 W m -2 on 1 February (Figure 2). Between 31 January and 1 February, the daily mean ground surface temperature increased by 1.7 °C. Mean daily QE was the highest on 24 February (116.2 W m -2 ), when QH reached only 50.3 Wm -2 . A concurrent increase in air and ground surface temperature before 24 February pointed to active layer thawing (Figure 4). Between 24 and 25 February, the ground surface temperature decreased by -5.8 °C, indicating that intense evaporation might have led to energy depletion in ground surface layers, which resulted in the drop in surface temperature.  Daily mean Q G was the highest on 11 March, when it reached 23.9 W m −2 . The Q G on 11 March led to warming of the whole column from the ground surface to a depth of about 30 cm (Figure 4), with daily mean ground surface temperatures rising by 5.3 • C between 10 and 11 March, while the increase was by 0.2 • C at a depth of 30 cm. On 1 February, 24 February and 11 March, the change in daily mean ground temperature at a depth of 50 cm was never by more than 0.1 • C day −1 , even though between 1 and 4 February there was a gradual increase of ground temperature at a depth of 50 cm observed.
Active layer refreezing at the end of the thawing season occurred during 26-28 February (Figure 4). It is clear from Figure 2 that, from 25 to 27 February, daily mean values of Q G were negative, resulting in a mean heat loss of −5.9 W m −2 . Such a negative mean three-day value of Q G had not been seen since the start of the measurement. It was also connected to the lowest three-day value of Q net since 23 January.

Longer-Term Influence of Surface Energy Budget Components on Ground Temperature Variation
Relationships between ground temperatures and surface energy budget components were the strongest for ground surface temperature and decreased with depth (Table 3). For Q net , the highest correlation coefficients were 0.71 for ground surface temperature and ground temperature at a depth of 30 cm. For Q G , the difference between the strength of correlations with ground surface temperature and ground temperature at 30 cm was nearly 0.2. For both Q net and Q G , correlation coefficients were the highest without a lag for ground surface temperature. Correlation coefficients were close to zero for a depth of 75 cm, which was connected to a nearly isothermal temperature regime with no daily and seasonal variation. Table 3. Daily average correlation coefficients between ground surface temperatures (T surf ), ground temperatures at depth of 30 cm (GT 30 ), 50 cm (GT 50 ) and 75 cm (GT 75 ) and net radiation (Q net ), ground heat flux (Q G ) and air temperature (T air ) on James Ross Island during 23 January to 21 March 2018. The labels "0-3 days" correspond to the lag of ground temperature time series. Daily mean Q G also had several significant and strong relationships with other fluxes (Table 4). For instance, a significant anticorrelation between Q G and albedo meant that snow cover led to more negative Q G , hence the flux from deeper levels towards the ground surface. Furthermore, there was a strong positive relationship between Q G and air temperature. Daily mean Q G could be expressed by the following regression functions: Q G = 5.7 + 2.1T air (5) and

Variable
respectively. Consequently, an increase in the summer air temperature of 1 • C would lead to an increase in summer Q G on JRI by 2.1 W m −2 . The relationship between air temperature and ground surface temperature was even stronger than the relationship between air temperature and Q G , as the correlation coefficient reached 0.95 (Table 3). The reason is that they both primarily respond to Q net . Consistent with Fourier [42], the lag between Q G and ground temperature increased with depth. The correlations were quite high probably due to lack of vegetation cover and dry soil, enabling fast heat transport between air and the ground.

Impact of Snow Cover and Cloudiness
The impact of snow cover and cloudiness was studied with the use of day-long case studies representing contrasting situations that occurred during the summer season. On 4 March (Figure 5b), lower-than-average Q L↑ (−284.4 W m −2 ) connected to low surface temperature (−2.7 • C) compensated for losses due to reflection from snow cover. Q L↓ was, for most of the time, 0-5 W m −2 higher compared with the average daily cycle (Figure 5a) due to the cloudy sky. Q net was below average apart from a short period between 1300 and 1400 local time when Q S↓ increased to 940 W m −2 because of a sudden reduction in cloud cover. Daily mean Q G on 4 March was −12.7 W m −2 , compared with the mean value of 0.4 W m −2 ( Table 1). The decrease between 4 and 5 March was −4.1 • C day −1 for air temperature and −4.9 • C day −1 for ground surface temperature (Figure 4). Consequently, both the surface and deeper ground levels were losing heat, despite the isolating effect of snow. Atmosphere 2020, 11, x FOR PEER REVIEW 11 of 18 The effect of cloudiness was studied on 18 and 15 February, when albedo <0.1 confirmed snowfree conditions. On the overcast day (18 February, Figure 6a), QL ↓ was, for most of the time, 10-25 W m -2 higher compared with the average daily cycle (Figure 5a). Daily mean Qnet was 44.8 W m -2 lower than on the averaged day. QG responded to cloudiness by negative values apart from an increase to 10.5 W m -2 around 1400. Consequently, the mean heat loss of the ground on 18 February was only -11.3 W m -2 , causing a drop of 2.4 °C in daily mean ground temperature at 10 cm depth (Figure 4). On the other hand, on 15 February, the sky was mostly clear apart from one episode of clouds evidenced by sudden decreases in QS ↓ at 1000 (Figure 6b). The large value of QS ↓ led to greater daily amplitudes of both Qnet and QG with daily means of 171.8 W m -2 and 6.4 W m -2 , respectively. The supply of heat to the ground was reflected in ground temperature increases between 15 and 16 February from the surface down to 20 cm depth (Figure 4). A sudden increase in QE at 1430 was connected to an increase of air temperature above 0°C (Figure 6b). The effect of cloudiness was studied on 18 and 15 February, when albedo <0.1 confirmed snow-free conditions. On the overcast day (18 February, Figure 6a), Q L↓ was, for most of the time, 10-25 W m −2 higher compared with the average daily cycle (Figure 5a). Daily mean Q net was 44.8 W m −2 lower than on the averaged day. Q G responded to cloudiness by negative values apart from an increase to 10.5 W m −2 around 1400. Consequently, the mean heat loss of the ground on 18 February was only −11.3 W m −2 , causing a drop of 2.4 • C in daily mean ground temperature at 10 cm depth (Figure 4). On the other hand, on 15 February, the sky was mostly clear apart from one episode of clouds evidenced by sudden decreases in Q S↓ at 1000 (Figure 6b). The large value of Q S↓ led to greater daily amplitudes of both Q net and Q G with daily means of 171.8 W m −2 and 6.4 W m −2 , respectively. The supply of heat to the ground was reflected in ground temperature increases between 15 and 16 February from the surface down to 20 cm depth (Figure 4). A sudden increase in Q E at 1430 was connected to an increase of air temperature above 0 • C (Figure 6b).

Surface Energy Budget Components in Polar Regions
The values of Qnet, QH, QE and QG from extant surface energy budget studies were only compared for February (Southern Hemisphere) and August (Northern Hemisphere) due to different lengths of measurement periods (Figure 7). Qnet on northern JRI for this month was 46-54% higher than on KGI [13] or dry tundra on Svalbard [43] or Axel Heiberg Island in northern Canada [44]. The only exception was the study from KGI by Prosek et al. [17] where Qnet was only about 10% lower than on JRI. This can be attributed to higher cloudiness on the western side of the AP and a consequent decrease in QS ↓, as can be evidenced by a lower ratio of Qnet to QS ↓ in this study (0.56) than in the study from KGI (0.62) by Choi et al. [13].
The value of QG on JRI was also slightly higher than on KGI [17], however, in both cases, QG comprised about 5% of Qnet on average. The amount of QG received by the soil decreases during summer in the ice-free polar regions (e.g., [13,24,29,38]), for instance by -0.46 W m -2 per 10 days on KGI [17]. Therefore, it is possible that around summer solstice, the soil on JRI might have received up to 5 W m -2 . The ratio of QG to Qnet on JRI was 2-4 times lower than in dry Arctic tundra, where summer QG usually formed 10-20% of Qnet (e.g., [23,45,46]), which could have been connected to differences

Surface Energy Budget Components in Polar Regions
The values of Q net , Q H , Q E and Q G from extant surface energy budget studies were only compared for February (Southern Hemisphere) and August (Northern Hemisphere) due to different lengths of measurement periods (Figure 7). Q net on northern JRI for this month was 46-54% higher than on KGI [13] or dry tundra on Svalbard [43] or Axel Heiberg Island in northern Canada [44]. The only exception was the study from KGI by Prosek et al. [17] where Q net was only about 10% lower than on JRI. This can be attributed to higher cloudiness on the western side of the AP and a consequent decrease in Q S↓, as can be evidenced by a lower ratio of Q net to Q S↓ in this study (0.56) than in the study from KGI (0.62) by Choi et al. [13].
in ground surface properties such as vegetation cover and polar day occurrence. In Choi et al. [13], QG formed nearly 40% of Qnet on KGI, but the value was most likely overestimated, as it also contained the residual term. Both on JRI and KGI [17] in February, QH was comparable to QE. Similar to dry Arctic tundra in late summer, QH on JRI also formed about 30% of Qnet (e.g., [44,47]). Summer QE in the Arctic were higher than QH by 20% to 100% [40] with the difference increasing as the autumn approached (Table  4 in [24,48]). That also seemed to be the case on JRI, as February mean QE was about 20% higher than QH.
Mean QE was higher than QH ( Figure 3) when the wind was from the westerly sector, while it was vice versa for the south-south-westerly wind. For the westerly sector, there was a clear influence of fetch from the Prince Gustav Channel, which was often free of sea ice. For the south-south-westerly wind, it is likely that measured QE would be higher than calculated QE, as the ground is waterlogged in that direction (Figure 1).

Surface Energy Budget Impact on Ground Thermal Regime
The main heat source for ground warming on JRI was Qnet and QS ↓ , as evidenced by the cooccurrence of the most intensive ground warming with the largest Qnet ( Figure 2 and Figure 4). For the most part, Qnet supplied heat to QG, both turbulent fluxes and ground warming. However, in 53% of the time, the ground was also contributing via QG to warming the air, which also occasionally happened in the dry Arctic (western Greenland) [45]. On 11 March, Qnet and QH were the main energy sources, as the ground was cooler than the air above ( Figure 2). On that day, the available energy was penetrating well via the frozen ground due to higher thermal conductivity of frozen compared to unfrozen ground, which led to slight ground warming (Figure 4). From the surface energy budget components, the closest relationship of ground temperature was found to be with QG.
The energy transfer via the ground was similar on the eastern and western sides of the AP, as demonstrated, for instance, by the correlations between Qnet and ground temperatures on JRI (this study) and KGI [17]. At the measurement site on northern JRI, there was bare ground, while Prosek et al. [17] carried out measurements over regularly flooded ground covered with a 3-cm-thick vegetation cover consisting of 30% moss and 70% grass. It is, therefore, clear that the insulating effect of vegetation cover acted to delay the energy transfer, which might have consequently operated at the level of days rather than hours.
The specific climate conditions and the absence of vegetation cover on JRI study site strongly influenced the ground thermal regime, as the ground temperature at a depth of 10 cm was about 4 °C lower than for the corresponding period in Arctic lichen tundra (Figure 7 in [49]). The ground The value of Q G on JRI was also slightly higher than on KGI [17], however, in both cases, Q G comprised about 5% of Q net on average. The amount of Q G received by the soil decreases during summer in the ice-free polar regions (e.g., [13,24,29,38]), for instance by −0.46 W m −2 per 10 days on KGI [17]. Therefore, it is possible that around summer solstice, the soil on JRI might have received up to 5 W m −2 . The ratio of Q G to Q net on JRI was 2-4 times lower than in dry Arctic tundra, where summer Q G usually formed 10-20% of Q net (e.g., [23,45,46]), which could have been connected to differences in ground surface properties such as vegetation cover and polar day occurrence. In Choi et al. [13], Q G formed nearly 40% of Q net on KGI, but the value was most likely overestimated, as it also contained the residual term.
Both on JRI and KGI [17] in February, Q H was comparable to Q E . Similar to dry Arctic tundra in late summer, Q H on JRI also formed about 30% of Q net (e.g., [44,47]). Summer Q E in the Arctic were higher than Q H by 20% to 100% [40] with the difference increasing as the autumn approached (Table 4 in [24,48]). That also seemed to be the case on JRI, as February mean Q E was about 20% higher than Q H .
Mean Q E was higher than Q H (Figure 3) when the wind was from the westerly sector, while it was vice versa for the south-south-westerly wind. For the westerly sector, there was a clear influence of fetch from the Prince Gustav Channel, which was often free of sea ice. For the south-south-westerly wind, it is likely that measured Q E would be higher than calculated Q E , as the ground is waterlogged in that direction (Figure 1).

Surface Energy Budget Impact on Ground Thermal Regime
The main heat source for ground warming on JRI was Q net and Q S↓ , as evidenced by the co-occurrence of the most intensive ground warming with the largest Q net (Figures 2 and 4). For the most part, Q net supplied heat to Q G , both turbulent fluxes and ground warming. However, in 53% of the time, the ground was also contributing via Q G to warming the air, which also occasionally happened in the dry Arctic (western Greenland) [45]. On 11 March, Q net and Q H were the main energy sources, as the ground was cooler than the air above ( Figure 2). On that day, the available energy was penetrating well via the frozen ground due to higher thermal conductivity of frozen compared to unfrozen ground, which led to slight ground warming (Figure 4). From the surface energy budget components, the closest relationship of ground temperature was found to be with Q G .
The energy transfer via the ground was similar on the eastern and western sides of the AP, as demonstrated, for instance, by the correlations between Q net and ground temperatures on JRI (this study) and KGI [17]. At the measurement site on northern JRI, there was bare ground, while Prosek et al. [17] carried out measurements over regularly flooded ground covered with a 3-cm-thick vegetation cover consisting of 30% moss and 70% grass. It is, therefore, clear that the insulating effect of vegetation cover acted to delay the energy transfer, which might have consequently operated at the level of days rather than hours.
The specific climate conditions and the absence of vegetation cover on JRI study site strongly influenced the ground thermal regime, as the ground temperature at a depth of 10 cm was about 4 • C lower than for the corresponding period in Arctic lichen tundra (Figure 7 in [49]). The ground without vegetation on JRI cooled more during the night, when the ground temperature at a depth of 10 cm fell to −0.4 • C, compared with a minimum of 2.3 • C in the Arctic.
The ground thawed and froze frequently both on JRI and KGI [17], but while Prosek et al. [17] observed refreezing of the whole active layer profile on KGI within 9 days, on JRI it was only 3 days.

Snow Cover and Cloudiness Effects
On JRI, the insulating effect of snow on Q G and the ground thermal regime was considered as rather negligible, both in this study and by Hrbacek et al. [50]. However, the maximum snow height reached only 26 cm during 23 January to 20 March, 2018 on JRI. Langer et al. [29] noticed that, at the Lena River Delta, the insulating effect was the most pronounced during winter, while, later in the year, the heat stored in the soil was depleted, meaning that snow insulation was of minor relevance. Furthermore, the high albedo of snow prevented the ground from receiving Q S↓ [44], which would have reduced ground cooling.
The radiation paradox [51] was not observed on JRI, as the effect of cloudiness was to cool the ground surface and soil layer, since the decrease in Q S↓ was not compensated for by an increase in Q L↓ . Such conditions also occurred in summer on Svalbard [52]. For individual days, the difference in Q net and Q G , due to lower cloudiness on JRI, could have been up to 120 and 17.5 W m −2 , respectively, while the difference in Q net was only 100 W m −2 across the tundra in Svalbard [43]. The cloud cover is usually less dense on the eastern side of the AP, as is documented, for instance, by the difference of 32 W m −2 in Q L↓ between this study and Hurd Peninsula glaciers on Livingston Island [14]. That could mean a significant increase in Q net and Q G on JRI with respect to the western side of the AP.
Neither the snow nor cloudiness affected Q H in summer significantly, both on the eastern and western sides of the AP [13]. One reason might have been that Q H responses to wind speed, as was shown earlier (Figure 3), and air-surface temperature differences, which do not always have to be connected to cloud cover development. Unlike on KGI [13], less energy was partitioned to Q E on the day with snow than on the day with extensive cloud cover.

Conclusions
In this study, the surface energy budget and its impact on the ground thermal regime were studied on JRI during January-March 2018. The main findings were: • Mean Q net reached 102.5 W m −2 , while the highest mean daily Q net was 253 W m −2 . Mean Q G was only 0.5 W m −2 and yet 46% of the time the ground was a heat source for the atmosphere.
Mean Q E was only by 0.4 Wm −2 higher than mean Q H . Mean Q E was up to 19.6 W m −2 higher than mean Q H when the wind blew from west-south-westerly to north-west-northerly sector (236.5 • to 348.5 • ), showing the influence of increased moisture availability from the sea.

•
The ground thermal regime was affected by surface energy budget components to the depth of 50 cm. The strongest relationship was found between Q G and the ground surface temperature, with a delay growing along with increasing depth. The active layer refroze at the end of February after a sequence of three days with continuously negative Q G .

•
The case studies have shown that an increase of cloud cover led to a decrease of both mean daily Q net and Q G which caused cooling in ground thermal profile. On a clear-sky day, the situation was vice versa.
• Insulation by snow cover was not observed, as Q G was 13.2 W m −2 below average and the ground was cooling in the whole profile during a sample day with snow cover.

•
By comparison with other studies, we concluded that Q net was higher on the eastern side of the AP due to less cloudiness affected by regional atmospheric patterns. Mean Q G reached similar values on the eastern and western side of the AP and comprised about 5% of Q net , approximately four-times less than was observed in the Arctic summer months. Mean Q E was, on both sides of the AP, approximately the same as Q H , even though the ratio of Q H to Q net varied by 15% among the studies.
Summer energy fluxes (Q H and Q E ) are crucial for the AP glaciers due to warming effect and efficient snow and ice melting (e.g., [9,14]), although they were mostly heat sources for the atmosphere in our study site. Ice-free areas might, accordingly, be an important source of energy in the AP region and should be considered in regional climate models. As all the surface energy budget studies from ice-free regions in Antarctica so far have only been seasonal, a deeper understanding of the surface energy budget would be gained from a full-year data set. The reported climate evolution of the AP region is currently slightly uncertain [8]. Should the warming observed since the 1950s continue, existing ice-free areas are more likely to be enlarged. Consequently, a study of the interannual and spatial variability in the surface energy budget in ice-free areas in Antarctica would provide valuable new information about how these new environments might influence the land-atmosphere interactions and the future climate of Antarctica.
The data that support the findings of this study are available from the corresponding author upon reasonable request.