Evolution of Backscattering Coefficients of Drifting Multi-Year Sea Ice during End of Melting and Onset of Freeze-up in the Western Beaufort Sea

Backscattering coefficients of Sentinel-1 synthetic aperture radar (SAR) data of drifting multi-year sea ice in the western Beaufort Sea during the transition period between the end of melting and onset of freeze-up are analyzed, in terms of the incidence angle dependence and temporal variation. The mobile sea ice surface is tracked down in a 1 km by 1 km region centered at a GPS tracker, which was installed during a field campaign in August 2019. A total of 24 Sentinel-1 images spanning 17 days are used and the incidence angle dependence in HH- and HV-polarization are −0.24 dB/deg and −0.10 dB/deg, respectively. Hummocks and recently frozen melt ponds seem to cause the mixture behavior of surface and volume scattering. The normalized backscattering coefficients in HH polarization gradually increased in time at a rate of 0.15 dB/day, whereas the HV-polarization was relatively flat. The air temperature from the ERA5 hourly reanalysis data has a strong negative relation with the increasing trend of the normalized backscattering coefficients in HH-polarization. The result of this study is expected to complement other previous studies which focused on winter or summer seasons in other regions of the Arctic Ocean.


Introduction
The Arctic Sea ice cover plays an important role in mediating the global climate through the ice-albedo feedback mechanism [1]. Over recent decades, the Arctic climate has been shifting rapidly and sea ice extent and thickness have substantially declined with more and more extreme cases [2][3][4][5]. In particular, the September extent has shown an exceptional decline trend with higher frequency of the lowest extents and prediction of the ice-free September Arctic Ocean after the year 2100 [6][7][8]. As sea ice extent and thickness decreases, thick multi-year sea ice (MYI) has been replaced by thinner first-year sea ice (FYI) [5,9].
Active remote sensing that operates in the microwave portion of the electromagnetic spectrum has been used to produce sea ice information. Synthetic aperture radar (SAR) for the most part has been an This study focuses on a large sea ice floe in the marginal ice zone of the western Beaufort Sea, north of Alaska, where a field campaign was carried on from August 12 to 15 in 2019 by ice breaker research vessel (IBRV) Araon (Figure 1a). This campaign was the first of two field campaigns that were carried out during the Arctic survey in 2019. Based on an airborne survey by a helicopter and field surveys, the surface of ice floe contained hummocks and wide melt ponds, of which surfaces were frozen at that time ( Figure 1b). The surface water drainage could occur well before the onset of freeze-up on FYI, as early as in June [27,28]; however, the perforation was not observed within 100 m radius from the IBRV Araon based on the field campaign. The average air temperature during the field campaign, measured from IBRV Araon, was -0.89 ± 0.70 ℃. The sea ice floe was determined as a multi-year ice, based on the salinity, temperature, and density measurements of 4 ice cores with a diameter of 9 cm, which were extracted at 3 different sites. Each ice core was sliced into 10-cm pieces, where small holes were drilled amid each piece for temperature measurement immediately after the ice was extracted. Ice temperatures were measured in those small holes, using a digital thermometer with an accuracy of ±0.1℃. The salinity was measured using a YSI Pro30 salinometer with an accuracy of 0.1 ppt. The upper part of the ice floe consisted of a fresh low-density layer (600 kg・m -3 for the lowest) lying over higher density layers, and salinity increases to 3 ppt with the thickest ice core of 3.5 m (Figure 2a   recorded every hour and remotely transmitted to the web-based server (www.spotmyglobalstar.com) via satellite communication. The operating temperature range of the buoy was between -30 and 60 ℃ and the expected operation period was 12 months. The tracker was firmly fixed in the ice using wood frames and weights, to prevent loss from strong wind or marine mammals. The time and coordinates of the tracker on the ice floe were used to select Sentinel 1A/B satellite images that simultaneously capture the sea ice floe within its swath.

Processing of Sentinel-1 C-Band SAR Data
Sentinel-1 extra wide (EW) Ground Range Detected Medium Resolution (GRDM) SAR images in HH-and HV-polarization covers over a 400 km swath at 40 m spatial resolution [29]. The incidence angle ranges from 18.9° to 47.0°, combining 5 sub-swaths. Sentinel-1 images are radiometrically and geometrically calibrated using European Space Agency's Sentinel-1 Toolbox software, so that the pixel values can be related to the radar backscatter of the sea ice surface. In the radiometric calibration, the backscattering coefficients of Sentinel-1 data are retrieved after applying the precise orbit information and removing the thermal noise. The backscattering coefficients are

Installation of a GPS Tracker on a Sea Ice Floe
A sea ice drifter buoy, developed by the Yotta Ocean in South Korea, with an integrated GNSS positioning system and Argos satellite-based data transmission system, was installed on the surface of the sea ice floe during the field campaign on 12 August, 2019. The coordinates of the tracker were recorded every hour and remotely transmitted to the web-based server (www.spotmyglobalstar.com) via satellite communication. The operating temperature range of the buoy was between -30 and 60°C and the expected operation period was 12 months. The tracker was firmly fixed in the ice using wood frames and weights, to prevent loss from strong wind or marine mammals. The time and coordinates of the tracker on the ice floe were used to select Sentinel 1A/B satellite images that simultaneously capture the sea ice floe within its swath.

Processing of Sentinel-1 C-Band SAR Data
Sentinel-1 extra wide (EW) Ground Range Detected Medium Resolution (GRDM) SAR images in HH-and HV-polarization covers over a 400 km swath at 40 m spatial resolution [29]. The incidence angle ranges from 18.9 • to 47.0 • , combining 5 sub-swaths. Sentinel-1 images are radiometrically and geometrically calibrated using European Space Agency's Sentinel-1 Toolbox software, so that the pixel values can be related to the radar backscatter of the sea ice surface. In the radiometric calibration, the backscattering coefficients of Sentinel-1 data are retrieved after applying the precise orbit information and removing the thermal noise. The backscattering coefficients are filtered to reduce speckle noise using a Lee filter with a window size of 11 by 11 pixels. The radiometrically calibrated images are then projected to the WGS84 NSIDC Sea Ice Polar Stereographic North (EPSG:3413), with a pixel size of 100 m. For every radiometrically and geometrically calibrated Sentinel-1 EW image in HH and HV-polarization, the backscattering coefficients and incidence angle at corresponding coordinates of the GPS tracker are extracted. The range of incidence angle is expected to be well distributed, benefitting from swath overlap at high latitudes [29]. Sentinel-1 data with the temporal difference from the GPS point less than 1-hour are selected, in order to minimize spatial error from continuously drifting sea ice. For example, 6-hour difference from the tracker creates over 10 km offset from the actual ice floe. Besides, interpolation of the recorded GPS coordinates is not carried out, owing to non-linear drift (speed and rotation) of the ice floe. Then, mean (shown as σ 0 pp in Tables 1 and 2) and standard deviation (shown as S(σ 0 pp ) in Tables 1 and 2) of the backscattering coefficients around extracted coordinates with a radius of 500 m are calculated, where pp denotes polarization, HH or HV. The sample window is decided for comparison with previous studies (i.e., Mahmud et al. [21]; Mäkynen and Karvonen [30]), and the spatial constraint is necessary to isolate from FYI and large open water, which could possibly be present around the ice floe of interest.
The incidence angle dependence is modelled as a first-order polynomial, with respect to the backscattering coefficients of HH-and HV-polarization [18]. The R 2 and root mean square error between the observed and modelled backscattering coefficients are calculated for the goodness of fit. The p-value for the coefficient is determined to test significance level. The backscattering coefficients are normalized (shown as N(σ 0 pp ) in Tables 1 and 2) for further analysis, using the derived regression model with respect to a particular reference angle, as in Mahmud et al. [21] In this study, the reference angle for the normalization of the backscattering coefficients is 33 • , the middle-angle of Sentinel-1 EW image [18]. Then, a fist-order polynomial regression model is used to determine a temporal trend of the normalized backscattering coefficients over the surface of multi-year sea ice during the observation period. The R 2 and root mean square error between the observed and predicted backscattering coefficients are calculated. The p-value for the coefficient is determined to test the significance level.

Atmospheric Effects on the Backscattering Coefficients
2-m temperature, instantaneous 10-m wind gust, and mean sea level pressure from the ERA5 hourly reanalysis dataset on single levels are analyzed, in an attempt to comprehend the weather-induced variation of backscattering coefficients. ERA5 is the fifth-generation atmospheric reanalysis of the global climate by the European Centre for Medium-range Weather Forecasts (ECMWF), which has a horizontal resolution of 0.25 • by 0.25 • and temporal resolution of 1 hour [31]. The hourly profile of each atmospheric variable at the coordinates of the installed GPS tracker is extracted by the linear interpolation between the grids using the recorded GPS data. Then, the parameter values of the atmospheric variables at corresponding times of Sentinel-1 data are estimated by linear interpolation between the profile points. Finally, 13-day records of air temperature and pressure from a buoy, installed in the inner part of the ice floe, are used to validate the reanalysis datasets. The recording period was much shorter than the GPS tracker, because the buoy was recovered on the way back owing to the frequent visit of a polar bear at the deployed site. The Spearman correlation coefficient with p-value is used to analyze the degree of association between the Sentinel-1 data and the atmospheric variables. The atmospheric variables from the reanalysis dataset are chosen based on the availability of in situ measurements for validation and the weather report from IBRV Araon during the field survey.

Results
The GPS tracker successfully transmitted its coordinates every hour from 12 August, 2019, and kept on recording its location until November of the same year. We visually checked every available Sentinel-1 EW which showed a continuous decrease in the ice concentration around the ice floe. A total of 24 images spanning 17 days, from 15 to 31 August, 2019, were selected for analysis. Figure 3 shows radiometrically calibrated Sentinel-1A/B images in HH-polarization which were used for the analysis of the incidence angle dependence, and later of temporal variation. A red square inside each image indicates a sample window of 1 km by 1 km which backscatter statistics were calculated. It is clear that the incidence angle induces substantial difference between the images that were acquired on the same date. It became unclear whether the ice floe completely melted, or the GPS tracker was detached from an ice floe after 1 September 2019 (Figure 4b). Besides, the images between 12 and 15 August 2019 were not used, because they were contaminated by a strong backscattering signal of IBRV Araon in the vicinity of the tracker.
Remote Sens. 2020, 12, x FOR PEER REVIEW 6 of 15 The GPS tracker successfully transmitted its coordinates every hour from 12 August, 2019, and kept on recording its location until November of the same year. We visually checked every available Sentinel-1 EW which showed a continuous decrease in the ice concentration around the ice floe. A total of 24 images spanning 17 days, from 15 to 31 August, 2019, were selected for analysis. Figure 3 shows radiometrically calibrated Sentinel-1A/B images in HH-polarization which were used for the analysis of the incidence angle dependence, and later of temporal variation. A red square inside each image indicates a sample window of 1 km by 1 km which backscatter statistics were calculated. It is clear that the incidence angle induces substantial difference between the images that were acquired on the same date. It became unclear whether the ice floe completely melted, or the GPS tracker was detached from an ice floe after 1 September 2019 (Figure 4b). Besides, the images between 12 and 15 August 2019 were not used, because they were contaminated by a strong backscattering signal of IBRV Araon in the vicinity of the tracker. a descending orbit, centered at the coordinates of the GPS tracker on the ice floe, shown as a red square. The time difference between the tracker and Sentinel-1 image is 6 minutes and a bright signal near the GPS coordinates is IBRV Araon, which stayed at the site until 15 August, 2019, for a field survey on the ice floe. For 17 days, the ice floe drifted about 100 km to the southeast, moving around the edge of the Chukchi Plateau, where the Beaufort Gyre usually moves counterclockwise ( Figure  1a). The total distance of drift was 245 km, with the highest instantaneous drifting speed being 2 km/h on 25 August in 2019. Notably,  0 HH and 0 HV of Sentinel-1 over the multi-year sea ice surface showed a strong linear relationship with the incidence angle.  0 are plotted against the well-distributed incidence angles ranging from 21° and 46°, with the standard deviation as error bars (Figure 5a). The minimum, maximum, and mean of  0 HH are −14.85 dB, −7.89dB, and −11.39 dB, with the standard deviation of 2.00 dB, respectively; and  0 HV are −24.54 dB, −21.31 dB, and −23.07 dB, with the standard deviation of 0.98 dB, respectively ( Table 1). The spatial standard deviations of HH-and HV-polarization ranged from 0.28 and 0.39 dB to 4.96 and 4.79 dB.
Then,  0 HH and  0 HV decreased as the incidence angle increased with a slope of −0.24 ± 0.02 dB/deg and −0.10 ± 0.02 dB/deg, respectively. The root mean square error from the linear least squares fitting line was 0.73 dB and 0.66 dB, with R 2 of 0.82 and 0.59 for HH and HV, respectively, with the confidence level of 99 percent (p-value < 0.01) ( Table 2).
Using the derived polynomial regression model as described in Section 2.2.3, radiometrically calibrated backscattering coefficients were normalized with respect to the reference angle of 33°, the middle-incidence angle of Sentinel-1 EW mode. After normalization, the difference between minimum and maximum and standard deviation of the mean backscattering coefficients in HH-polarization, N( 0 HH), decreased, as shown in Table 1. The mean backscattering coefficients in HV-polarization, N( 0 HV), exhibited subtle changes. The difference between the minimum and  Notably, σ 0 HH and 0 HV of Sentinel-1 over the multi-year sea ice surface showed a strong linear relationship with the incidence angle. σ 0 are plotted against the well-distributed incidence angles ranging from 21 • and 46 • , with the standard deviation as error bars (Figure 5a). The minimum, maximum, and mean of σ 0 HH are −14.85 dB, −7.89dB, and −11.39 dB, with the standard deviation of 2.00 dB, respectively; and σ 0 HV are −24.54 dB, −21.31 dB, and −23.07 dB, with the standard deviation of 0.98 dB, respectively ( Table 1). The spatial standard deviations of HH-and HV-polarization ranged from 0.28 and 0.39 dB to 4.96 and 4.79 dB.
Then, σ 0 HH and σ 0 HV decreased as the incidence angle increased with a slope of −0.24 ± 0.02 dB/deg and −0.10 ± 0.02 dB/deg, respectively. The root mean square error from the linear least squares fitting line was 0.73 dB and 0.66 dB, with R 2 of 0.82 and 0.59 for HH and HV, respectively, with the confidence level of 99 percent (p-value < 0.01) ( Table 2).
Using the derived polynomial regression model as described in Section 2.2.3, radiometrically calibrated backscattering coefficients were normalized with respect to the reference angle of 33 • , the middle-incidence angle of Sentinel-1 EW mode. After normalization, the difference between minimum and maximum and standard deviation of the mean backscattering coefficients in HH-polarization, N(σ 0 HH ), decreased, as shown in Table 1. The mean backscattering coefficients in HV-polarization, N(σ 0 HV ), exhibited subtle changes. The difference between the minimum and maximum normalized backscattering coefficients of the ice surface was around 3 dB both in HH-and HV-polarization. Table 1. Minimum, maximum, mean, and standard deviation of the un-normalized (σ 0 HH ), normalized (N(σ 0 HH )) backscattering coefficients and spatial standard deviations (S(σ 0 HH )) in a 1 km by 1 km region of interest in HH-and HV-polarization.

Min [dB]
Max  Table 2. The statistical relations, R 2 , the Spearman coefficient ρ, and p-value, for the backscattering coefficients (σ 0 ), the normalized backscattering coefficients (N(σ 0 )), and the spatial standard deviations (S(σ 0 )), in both polarizations with the atmospheric variables.   (Table 2). R 2 was 0.76 and the root mean square error was 0.44 dB, when fitted using linear least squares regression (Figure 5b). N(σ 0 HV ) also showed an increasing trend as time progressed, although the slope was relatively small, 0.06 ± 0.02 dB/day with R 2 of 0.24 and root mean square error of 0.53 dB with a confidence level of 95 percent (Table 2), when fitted using linear least squares regression. Figure 6 shows the time series of 2-m air temperature, maximum instantaneous wind gust at 10 m, and mean sea level pressure from the hourly ERA5 reanalysis dataset, with the in situ air temperature and pressure measured at a height of 2 m on a buoy (dashed red lines). The reanalysis datasets correlate well with the in situ measurements. The research period for which Sentinel-1 SAR data were gathered is denoted by a horizontal black bar at the bottom of each panel. Black vertical lines over the horizontal bar indicate the acquisition times of Sentinel-1A/B data. Some lines are indistinguishable due to close proximity between Sentinel-1 A and B acquisitions. The air temperature is smoothed by a Gaussian filter with a width of 24-hour to represent a general pattern. The air temperature fell below freezing on 13 Aug 2019 for the first and the diurnal variation stayed below freezing on 20 August 2019 and thereafter. It should be noted that the wind gust was considerably high after 20 August 2019, with strong low pressure. The highest wind gust was 21.72 m/s on 25 August 2019, when the temperature and pressure dropped largely ( Figure 6). The statistical relations, the Spearman coefficient  and p-value, between the normalized backscattering coefficients in both polarizations, as well as the spatial standard deviations and the atmospheric variables, are listed in Table 2. The normalized backscattering coefficients in both polarizations are strongly correlated with the air temperature. Especially, N( 0 HH) shows stronger negative trend with higher confidence level than N( 0 HV). The correlation coefficient between S( 0 HV) and the pressure at mean sea level show moderate, but negative, relation at the confidence level of 95 percent (p-value = 0.015). The rest show low correlation. The research period for which Sentinel-1 SAR data were gathered is denoted by a horizontal black bar at the bottom of each panel. Black vertical lines over the horizontal bar indicate the acquisition times of Sentinel-1A/B data. Some lines are indistinguishable due to close proximity between Sentinel-1 A and B acquisitions. The air temperature is smoothed by a Gaussian filter with a width of 24-hour to represent a general pattern. The air temperature fell below freezing on 13 Aug 2019 for the first and the diurnal variation stayed below freezing on 20 August 2019 and thereafter. It should be noted that the wind gust was considerably high after 20 August 2019, with strong low pressure. The highest wind gust was 21.72 m/s on 25 August 2019, when the temperature and pressure dropped largely ( Figure 6). The statistical relations, the Spearman coefficient ρ and p-value, between the normalized backscattering coefficients in both polarizations, as well as the spatial standard deviations and the atmospheric variables, are listed in Table 2. The normalized backscattering coefficients in both polarizations are strongly correlated with the air temperature. Especially, N(σ 0 HH ) shows stronger negative trend with higher confidence level than N(σ 0 HV ). The correlation coefficient between S(σ 0 HV ) and the pressure at mean sea level show moderate, but negative, relation at the confidence level of 95 percent (p-value = 0.015). The rest show low correlation.

Discussion
High  0 of the multi-year sea ice in the western Beaufort Sea in late August, both before and after normalization, represent the characteristics of multi-year sea ice well, mainly because of the

Discussion
High σ 0 of the multi-year sea ice in the western Beaufort Sea in late August, both before and after normalization, represent the characteristics of multi-year sea ice well, mainly because of the rough and tilted surfaces of the hummocks [32][33][34]. The incidence angle dependence of the co-polarized backscattering coefficients (σ 0 HH ) suggests that the coefficients during the study period are a combination of the surface and volume scattering. In general, the incidence angle dependence is weak (−0.16 dB/deg) when the volume scattering is dominant (i.e., multi-year sea ice in winter) [21] and becomes stronger as the surface scattering becomes more dominant in summer, attributed to the melting of the top layer [18,19,21,30,35]. In this study, the incidence angle dependence of the multi-year ice of interest still leans to that in summer, although the air temperature fell below freezing. The drop of the air temperature must have increased the volume scattering from bubbles in the upper part of the ice increases [24]. However, the contribution of volume scattering was not significant enough to dramatically weaken the incidence angle dependence.
The cross-polarized backscattering coefficients (σ 0 HV ) showed weaker incidence angle dependence than the co-polarized, also lower than that of Mäkynen and Karvonen [30] and Geldsetzer et al. [36]. However, the high magnitude of σ 0 HH , staying above the system noise equivalent for Sentinel-1 EW mode data throughout the research period, supports the relatively high confidence level of the regression model (Table 2) [30]. Previous results with σ 0 HV are scarce; thus, making a conclusion at this time is difficult. Besides, the analysis of the incidence angle dependence for σ 0 HV on neither FYI nor MYI has been done, whereas it is known that the incidence angle dependence for σ 0 HH is low when the volume scattering is dominant [21].
The increasing trend of σ 0 HH is likely linked to the atmospherically forced freezing events. Freezing of the snow and ice top layer would increase contributions of volume scattering from the coarse-grained snow-ice medium between the snow basal layer and ice top [19,37]. The air temperature fell below freezing at the acquisition of Sentinel-1 data (08~09 local time) and the statistical results show that the air temperature has strong negative relation with the increasing trend of N(σ 0 HH ) (r = −0.72, p-value = < 0.01). Previous research by Jeffries et al. [35] and Winebrenner et al. [24] showed a similar increasing trend in late August. They also measured the Arctic sea ice SAR signature during the melt-season in the western Beaufort Sea using ERS-1 SAR data (VV-pol), that also operated in a similar C-band (5.6 GHz). The central frequency of Sentinel-1 is 5.405 GHz.
Although the air temperature mostly stayed below freezing during the study period, it is worth checking whether the surface of melt pond was actually in a frozen state. The correlation between σ 0 and high wind gust during the study period was expected to be significantly high, because a strong surface wind can also substantially increase the C-Band backscattering [28,38]. However, the the surfaces of melt ponds were frozen at least during the field campaign, validated by visual observation, and the correlation was very low. In addition, there is about a 3 dB difference between 17 and 30 August 2019 when the wind gusts were similar. Therefore, it is unlikely that the surface of melt pond was covered with liquid water during the study period.
The spatial standard deviations of the co-polarized backscattering coefficients within the region of interest around the tracker infer homogeneity of the surface, expecting to be an indirect measurement of physical change of the surface. The correlation is expected to be high for the wind gust, because of the exceptionally high wind gust during the research period ( Figure 6); however, the correlation is weak, with the p-value showing some feasibility. Thus, wind gust has a subtle effect on the surface changes that affect the backscattering of the radar signal.
On the other hand, the high confidence level with a negative relation between the pressure and S(σ 0 HV ) is an intriguing result. Initially, the pressure was tested in an attempt to find any relationship between the SAR backscatter and the Arctic storm. Arctic storms have recently been of interest for sea ice owing to potential impact on the coupled Arctic sea-ice system [39,40]. In this context, the combined effect of high wind speed and possible precipitation at low pressure was expected to perturb the backscatter. However, an explanation for the high correlation between S(σ 0 HV ) and the pressure is still at large at the moment.
Finally, it should be noted that this ice floe, which was thicker than 3 m, likely disintegrated in about 2 weeks, despite strong evidence of the onset of freeze-up of the ice surface during the study period. The ocean heat in the Beaufort Sea can induce continuous bottom melting after surface melting has ended [41,42]. It has been reported that extreme events, such as the record low level of the sea ice extent in 2019, might have been caused by a large amount of ocean heat [43,44]. The average near-surface water temperature measured by a temperature sensor at a seawater intake near the bow of IBRV Araon between 13 and 16 August 2019 was −1.27°C ± 0.05°C. Considering that that the bottom temperature of sea ice is near −1.7°C, the ocean to ice-bottom heat flux is 71.54 Wm -2 [45]. Besides, the average sea surface temperature from the National Oceanic and Atmospheric Administration (NOAA) Optimum Interpolation Sea Surface Temperature version 2 (OISSTv2) product was 1~7°C warmer than the 1982-2019 August mean in the Beaufort Sea [46,47]. Although an in situ ice thickness measurement is not available, the thickness likely decreased over time and the large heat flux at the bottom might have affected the brine volume and dielectric constant in the ice layer, where the microwave interacts to increase the backscattering [20,41]. The effect of the warmer ocean on the backscattering requires additional investigation with regard to the geophysical properties of sea ice such as salinity and temperature profile in time series, along with in situ measurements of the environmental parameters and satellite images for a longer period of time.
Finally, there is feasibility for the incidence angle dependence as an indicator of which scattering is more dominant over sea ice. The incidence angle dependence is strong for the surface scattering, whereas the dependence is weak for volume scattering. We need to build up more datasets to analyze how fast the dominant scattering changes over time. In addition, the atmospheric forcing is probably the most significant driving factor, and yet, additional datasets are required. The melt pond fraction must be taken into consideration as well.

Conclusions
The substantial decline of the Arctic Sea ice extent and thickness in recent decades has replaced multi-year sea ice with thinner first-year sea ice, to a vast extent. The precise classification of multi-year and first-year sea ice can help better perceive the situation of on-going sea ice changes. High-resolution SAR data are known to distinguish such ice types; however, precise radiometric calibration is required regarding the incidence angle. This study tracked down the backscattering coefficients of a multi-year sea ice surface in motion, using a GPS tracker which was installed during a field campaign in August 2019 and 24 Sentinel-1 EW images spanning 17 days. The incidence angle dependence of the multi-year sea ice surface in HH-and HV-polarization were −0.24 dB/deg and −0.10 dB/deg, respectively. The surface with hummocks and recently frozen melt ponds had the incidence angle dependence of which were the mixture of surface and volume scattering. The normalized backscattering coefficients in HH polarization gradually increased in time at a rate of 0.15 dB/day as the shape of ice floe became indistinguishable. The temporal variation of the HV-polarization was relatively insignificant. Air temperature had strong negative relation with the normalized backscattering coefficients in HH-polarization, as the contribution of volume scattering increased. Other external atmospheric forces such as air pressure at mean sea level and wind gust unlikely had influences on the backscattering. Low correlation between the backscatter and wind gust suggests that the melt pond surfaces stayed frozen throughout the study period. Homogeneity of the surface had significant relation, only with the pressure at mean sea level, which was investigated using the standard deviations of the backscattering coefficients within the region of interest. Lastly, the most conspicuous finding is that this thick ice floe, which was thicker than 3 m, most likely disintegrated in about 2 weeks, despite strong evidence of the onset of freeze-up of the ice surface. The relatively warmer ocean temperature in the region and large heat flux might have induced continuous bottom melting. Any influence of the warm water on the dielectric constant in the layer where the microwave interacted will require additional investigation