Sensitivity of Radar Altimeter Waveform to Changes in Sea Ice Type at Resolution of Synthetic Aperture Radar

Radar altimetry in the context of sea ice has mostly been exploited to retrieve basin-scale information about sea ice thickness. In this paper, we investigate the sensitivity of altimetric waveforms to small-scale changes (a few hundred meters to about 10 km) of the sea ice surface. Near-coincidental synthetic aperture radar (SAR) imagery and CryoSat-2 altimetric data in the Beaufort Sea are used to identify and study the spatial evolution of altimeter waveforms over these features. Open water and thin ice features are easily identified because of their high peak power waveforms. Thicker ice features such as ridges and multiyear ice floes of a few hundred meters cause a response in the waveform. However, these changes are not reflected in freeboard estimates. Retrieval of robust freeboard estimates requires homogeneous floes in the order of 10 km along-track and a few kilometers to both sides across-track. We conclude that the combination of SAR imagery and altimeter data could improve the local sea ice picture by extending spatially scarce freeboard estimates to regions of similar SAR signature.


Introduction
Active radar remote sensing plays a key role for high-resolution observations of the Arctic sea ice cover in the order of a hundred meters to a kilometer. Synthetic aperture radar (SAR) imagery carries information, for example, about sea ice concentration and extent, sea ice types and sea ice drift, e.g., [1]. Radar altimetry is the primary tool to estimate sea ice thickness with good accuracy on basin-scale [2,3]. The altimeter measures the sea ice freeboard, i.e., the distance between the ocean and the sea ice surface, and ice thickness can then be calculated assuming hydrostatic equilibrium and auxiliary information about sea ice density, i.e., sea ice type, and snow load [4]. To reduce the dependence on auxiliary information, the direct retrieval of sea ice type from altimeter waveform data has been investigated and showed promising results, but also a large variability of waveform parameters [5][6][7][8]. The introduction of the delay/Doppler or SAR altimeter has significantly increased the along-track resolution of spaceborne altimetric measurements [9,10]. This results in better sensitivity to small features on the surface [11] and enables the ability to investigate the response of altimeter waveforms to small-scale variations of sea ice types such as large ridges, refrozen leads or multiyear sea ice floes embedded in first-year sea ice on spatial scales from a few hundred meters to a few tens of kilometers.
Altimetric data alone is, however, not sufficient to unambiguously determine the origin of waveform changes over sea ice, because of temporal and spatial variability of waveforms and the impact of off-nadir specular reflectors [6,7,12]. Near-coincidental SAR imagery has proven a valuable source for verification of lead detection algorithms in altimeter data processing, e.g., [13,14]. In the late 1980s, two studies have been undertaken with data from Seasat and Geosat to relate SAR backscatter and altimetric data over sea ice, but both were limited to one SAR image as a reference [15,16]. Presently however, Copernicus' Sentinel-1 and Sentinel-3 two satellite constellations in addition to Radarsat-2/Radarsat-Constellation and CryoSat-2 provide an increasing and long term (Copernicus program) data availability especially at lower latitudes (below 81.5 • ). Furthermore, near real time availability of data [17] facilitates the use for short term observations of the sea ice cover. This allows us to perform, to the best of our knowledge, the first comprehensive study of the waveform evolution over different sea ice types with SAR imagery as high-resolution validation data.
SAR imagery gives a high-resolution, two-dimensional picture of the current sea ice situation, but sea ice thickness can only be inferred from sea ice signatures with ambiguities [18]. Meanwhile, altimeter data provides a means to obtain sea ice thickness, but is sparse in nature due to the one-dimensionality of the measurements. To overcome these spatial limitations, averaging in space and time is commonly applied to obtain standard products of hemispherical sea ice freeboard and thickness, e.g., [19], at the expense of small-scale features. Combining these two information sources therefore has several potential benefits. On the one hand, sea ice type information, necessary for freeboard to thickness conversion [4], can be provided by SAR imagery in high-resolution [20] to reduce uncertainties in this critical step. On the other hand, the scarce thickness retrievals could be locally expanded to neighboring areas with similar sea ice signature in the SAR image regarding backscatter intensity and texture.
In this paper, we use near-coincidental Sentinel-1 SAR imagery and CryoSat-2 altimeter data of three winter seasons in the Beaufort Sea to investigate the sensitivity of altimeter waveforms to small-scale (up to 10-20 km) variations of sea ice. We evaluate whether SAR imagery aids the interpretation of the sea ice situation and identification of sea ice types. Furthermore, we derive dimensions of such features to obtain meaningful freeboard estimates. The study area and data are presented in Section 2. The backscatter and altimeter parameters are presented in Section 3 and analyzed in Section 4, focusing first on the distinct waveform properties of the different ice types (4.1) and then on local features (4.2). We conclude with a discussion on how SAR imagery and altimetry can be used to augment each other's ability for sea ice monitoring. Figure 1 shows the study area within the Beaufort Sea, located in the Arctic. Blue areas represent the Sentinel-1 images and the red lines the CryoSat-2 altimeter ground tracks. The Beaufort Sea was chosen for this study because first-year and multiyear sea ice are present during wintertime in this area. Furthermore, a belt of larger multiyear ice floes embedded in first-year ice can usually be found in this area and is thus ideal for studying the ability to detect small-scale changes in sea ice. We collected near-coincidental Sentinel-1 SAR imagery and CryoSat-2 altimeter data of three winter seasons (Feb/Mar 2016, Nov/Dec/Jan 2016/17 and Jan/Feb/Mar 2018) with an average time lag from 34 min ranging from a few to about 80 min. For drift speeds up to about 6 cm/s the motion will be below the along-track resolution of the altimeter (about 300 m) for all time lags. Higher drift speeds are possible in this area [21] and therefore sea ice drift can exceed the altimeter footprint size for larger time lags.

Sentinel-1 A/B
We used 73 Copernicus Sentinel-1 A/B dual-polarization C-band (5.405 GHz = 5.6 cm) SAR images taken between February 2016 and March 2018. Table 1 summarizes the temporal distribution of the data where a sequence denotes up to three sequential images connected to the same near-coincidental altimeter track. All images are extra wide swath ground range detected with medium resolution, i.e., 40 m × 40 m pixel size and range and azimuth resolution of 93 m and 87 m, respectively. The polarization is horizontal transmit and receive (co-polarization, HH) and horizontal transmit and vertical receive (cross-polarization, HV). The data is preprocessed using the Sentinel Application Toolbox (SNAP) applying the following steps: precise orbit → noise removal → (optional) slice assembly → calibration to σ 0 → terrain correction to UTM with ASTER 1 sec GDEM.

CryoSat-2
We used 35 CryoSat-2 altimeter tracks that were near-coinciding with the Sentinel-1 images to analyze waveform response and behavior for different ice types and small-scale features. The SIRAL altimeter aboard CryoSat-2 is a K u -band (13.575 GHz = 2.2 cm) altimeter that operates in SAR (delay/Doppler) mode over sea ice. Sequences of radar bursts are emitted by the radar forming a stack of waveforms that view the same area on the ground from different viewing angles. The range corrected waveforms are averaged, and these multi-looked waveforms are provided alongside so-called stack parameters that describe the angular behavior of delay time integrated power [10]. The pulse limited footprint size in SAR mode is approximately 300 m along-track and 1.5 km across-track [22]. The distance between range bins of the multi-looked waveforms is 23.42 cm and a full waveform of 256 bins covers a range window of 60 m, whereas the range resolution is about 47 cm [23]. The high orbital inclination of the satellite (92 • ) allows observations up to 88 • N in the Arctic. CryoSat data provided by ESA is used in this study. Waveform parameters are obtained from level 1b Baseline C data and freeboard measurements are taken from level 2 data. For the analysis in this study no further processing of the data is performed.

SAR Image and Waveform Parameters
Comparison of SAR imagery and altimeter waveform data requires the spatial overlap of both data sets. Consecutive SAR images of the same orbit were assembled into one single image and corresponding altimeter tracks located using the geolocation information of the product. Rectangular polygons with a length of 300 m in along-track and 1.5 km in across-track direction were calculated around each altimeter track point. This size resembles the footprint size of the altimeter associated with the first peak of the return [22]. The mean value of backscatter intensity σ 0 at co-and cross-polarization of this polygon is calculated for comparison with parameters describing the shape of the altimeter waveform.
SAR instruments and altimeters perceive sea ice surface information differently because of their viewing geometry, e.g., side-looking for SAR imagery and nadir-looking for distance measuring altimeters [24,25]. Smooth surfaces, i.e., flat open water or thin ice, have a low backscatter intensity in SAR imagery while altimeter waveforms show a narrow peak with high power. With increasing roughness of the surface, SAR backscatter intensities increase, and altimeter waveforms become more diffuse with lower peak power and larger power values in the tail of the waveform. Figure 2 shows examples of typical waveforms for leads, first-year and multiyear sea ice. Note the differences of the scaling of the y-axis.  In this study, we used pulse peakiness (PP), scaled inverse mean power (IMP) and stack standard deviation (SSD) as parameters to characterize the shape of the waveform. The first two are derived from the waveform while SSD is a stack parameter that is provided as auxiliary information with the data set. The maximum peak power is a helpful parameter for visual waveform interpretation but contains similar information as PP. PP is used in the quantitative analysis whereas peak power is used for qualitative description of spatial waveform evolution.
PP is a measure for how much power is contained in the maximum power compared to the mean power of the waveform and is calculated according to [26]: where WF i is the waveform and 256 is the number of bins of the waveform. PP is high for returns from smooth surfaces and decreases with surface roughness. SSD describes the angular variation in power for the stacked waveforms [10] and is low for cases of specular surface scattering and increases for rougher surfaces where returns from larger incidence angles contribute significantly to the signal. We also introduced a new parameter, scaled inverse mean power IMP that reflects the total power contained in a waveform and is calculated as follows: This parameter is scaled by 2 · 10 −13 to avoid too small values and hence increase readability. IMP shows a similar behavior as SSD but enhances the contrast when waveforms become more similar in terms of peak power. Leads are expected to have a high PP and low values of IMP and SSD because of the specular nature of the nadir return. First-year sea ice has an increased surface roughness but is smoother than multiyear sea ice, which decreases PP but increases IMP and SSD. Multiyear sea ice or rough surfaces therefore have the lowest PP but the highest IMP and SSD values.

Freeboard
Sea ice freeboard estimates are used to assess the potential to retrieve sea ice thickness from altimeter data over small-scale sea ice features. The procedure to convert the measured altimeter waveform to surface elevations is called retracking and an overview of commonly applied methods can be found in [27]. Common to all these methods is to identify the point on the leading edge of the waveform that corresponds to the nadir return. Freeboard estimates used in this study are solely from the provided CryoSat-2 level 2 product that uses a threshold retracker referenced to the first peak of a smoothed waveform [28]. Distortions of the leading edge, e.g., by off-nadir specular reflectors or mixtures of ice types, can cause false elevation estimates [12,29]. Figure 3 shows an example of waveform evolution over a multiyear sea ice floe taken on 18 February 2018. Time separation between the SAR image and altimeter measurement is just a few minutes. The circular, darker SAR signature between samples 502 and 505 probably represents a smoother ice surface. The highly negative freeboard for sample 502 is related to the deformed, slightly kinked leading edge. The strong peak can no longer be distinguished from the first peak as in the waveforms before and after. While this is an exaggerated example with highly negative freeboard, this can also cause an underestimation of sea ice thickness in less affected cases. Freeboard measurements represent radar freeboard, i.e., the scattering horizon of the radar waves that usually falls somewhere between the air-snow and snow-ice interface [26]. Sea ice thickness is of primary interest but estimation of it requires auxiliary information about snow load and physical and electromagnetic properties of the sea ice, which is beyond the scope of this study.

Sample Selection
First-year and multiyear samples were manually distinguished in the SAR imagery based on spatial context, texture and backscatter intensity. Multiyear sea ice (MYI) is characterized by high backscatter values in co-and cross-polarization, distinct texture and large floe sizes. First-year sea ice (FYI) has a much larger spectrum of characteristic SAR appearances that are usually lower in backscatter than for MYI especially at cross-polarization. High backscatter FYI, e.g., brash ice or frost flower covered ice, is distinct by its more homogeneous texture compared to MYI. Ice charts were used as auxiliary information in the decision process. Only those waveform samples that show a homogeneous ice type within the footprint area of the altimeter were considered for further analysis. Furthermore, specular returns of the altimeter where rejected based on PP and peak power to avoid the influence of smooth features below the resolution of the SAR images. We also included larger leads, i.e., areas in between distinct floes, with lower SAR backscatter that are large enough to contain the footprint area of the altimeter. These leads could have a high peak return because it is characteristic for the specular return of smooth ice. Table 2 summarizes the selected samples for each ice type and season. Differences in sample numbers between the seasons reflect differences in the available near co-incidental data sets.

Waveform Parameter Statistics
The first step of our analysis focusses on the altimeter waveform and SAR image parameters of different ice types. Figure 4 shows the waveform parameter statistics (a-c) and HV backscatter coefficient (d) of FYI, MYI and large leads accumulated over the three winter seasons 2015/16, 2016/17 and 2017/18. The three ice types are clearly distinguished by the waveform parameters, but the distributions also show a large overlap. FYI waveforms have a large PP, i.e., most of the backscattered power is located in the maximum peak, and a low SSD value, i.e., the contributions to the multi-looked waveform are confined to a smaller incidence angle range. Furthermore, FYI waveforms have a larger mean power compared to MYI as shown by the low values of IMP. MYI waveforms contain energy from a larger range of incidence angles, i.e., high SSD value, and the low PP indicates a smoother distribution of backscattered energy over the entire waveform. Leads are characterized by the highest PP, low SSD and high mean power. IMP enhances the contrast for FYI and MYI but reduces it for FYI and leads. FYI and MYI are well separated by the mean backscatter coefficient of the SAR image since the age of sea ice has a significant impact on the SAR backscatter intensity [24]. Large leads have a bimodal distribution of HV backscatter intensities that largely overlaps with the FYI values. The lower peak is caused by leads of new ice whereas the second peak at higher backscatter values represents leads that have started to form thicker ice but still contain areas of thin ice that cause a highly peaked response in the altimeter.  The robustness of these characteristics can be checked by investigating their temporal variability for each ice type. Figures 5 and 6 show the variation of waveform parameters and HV backscatter for FYI and MYI in the three different winter seasons. FYI shows a larger variability between the different years for PP and SSD waveform parameters, while the cross-polarization backscatter coefficient is relatively stable. IMP of FYI is less variable then for MYI because of generally higher peak power in FYI waveforms. The larger tail of PP in the winter 2015/16 represents waveforms with higher peak power. Part of the observed variability is probably induced by the fact that FYI is much more variable in surface appearances that affect the scattering process. For MYI the PP varies very little and only the winter 2016/17 shows a deviation for SSD and IMP. The peak at lower IMP originates from waveforms over MYI in the beginning of November 2016, and therefore relatively early in the season. Slightly more texture is visible in the SAR image and refrozen melt ponds could have increased the specular component. It should be noted here that the limitation to near-coincidental altimeter data and manual sample selection also introduce variabilities to the waveform parameters. The complexity of altimeter interaction with the sea ice surface and the dependence of altimeter waveforms on small-scale specular scatterers in the footprint, e.g., [25], make it difficult to assess the reasons of certain waveform changes based on the available SAR imagery. A more in-depth analysis of potential causes in yearly variability is therefore not feasible in this study.
Apart from the presented findings for different ice types, the following observations have been made for the discarded samples and are based on the comparison of the SAR images and the altimeter waveforms during the sampling process. FYI contains a large proportion of lead-like waveforms, i.e., waveforms with high PP, low SSD and high peak power that do not correspond to a recognizable lead in the SAR image. Therefore, more samples have been discarded over FYI because they could not be unambiguously attributed to FYI. Waveforms of mixtures of ice types within the altimeter footprint, e.g., transition from FYI to MYI at floe boundaries or small refrozen leads within MYI, resemble more the waveforms of FYI. More details can be found in Section 4.2. This holds also true for the transition from FYI to MYI areas where floe sizes are smaller, and a mixture of thicker FYI and second year ice is present.  Figure 7 shows the freeboard distributions of different ice types and winter seasons based on our selected samples. The solid lines represent kernel density estimates of the histograms. Samples with zero freeboard, 59% of FYI and 24% of MYI samples, as well as with negative freeboard, 3% of FYI and 6% of MYI samples, have been excluded from the analysis. The high percentage of zero freeboard samples is caused by samples that are relatively lead-like but have been retained in our analysis because of their low peak power. Negative freeboards are possible if an ice floe gets submerged by a heavy snow load, but are unlikely to be measured by the radar as scattering probably will occur before the ice-snow interface. Errors in retracking the correct surface elevation because of distorted waveforms are more likely especially for high negative values. Expectedly, MYI ice shows a higher freeboard compared to FYI but also has a significant portion in the range of FYI freeboard values. These are probably caused by the influence of off-nadir smooth areas below SAR resolution as demonstrated in Section 3.2 and could therefore not be avoided in the selection process. Further examples of freeboard underestimation for mixed ice types and apparent homogeneous areas in SAR images are shown in Section 4.2. We see that FYI freeboard rarely exceeds 40 cm for the selected samples. Freeboard estimates for FYI and MYI are stable between the three winter seasons and are therefore not significantly affected by the variations in waveform parameters. No significant difference is recognizable between the early winter measurements of 2016/17 compared to the other two winters. Slightly higher MYI freeboard for the winter season 2016/17 could be caused by samples of older sea ice compared to the other winters.
In this section, we demonstrated that waveform parameters principally distinguish between FYI, MYI and large leads. FYI shows larger variations of these parameters for three winter seasons compared to MYI, but has more consistent HV backscatter. The principal difference in waveform parameters over different ice types forms the basis to study the waveform response to small-scale sea ice features such as ridges and MYI floes, which we now present in Section 4.2.

Local Waveform Variability
Of particular interest is the response of altimeter waveforms to spatial changes of the sea ice surface. In the following we present representative examples over a large ridge, and for different sized MYI floes that are embedded in FYI. Described findings are typical of these features and occurred for many samples in our dataset.
The first example in Figure 8 shows the local waveform variation and SAR imagery over a large ridge on 2016-03-15. The time lag between the SAR image and the altimeter track is below five minutes, i.e., for average drift speeds between 6 cm/s and 10 cm/s that means 18 m to 30 m. The ridge is clearly visible in the center of the cross-polarization SAR image by its higher backscatter caused by ice facets oriented towards the radar (Figure 8a-c). The altimeter on the other hand shows a lower peak intensity around the ridge because more energy is scattered away from the nadir-looking instrument (Figure 8d). PP is low over the ridge, while SSD and inverse mean power reach higher values (Figure 8e,f). This is characteristic for more diffuse scattering. Waveforms over FYI have a higher peak power with high PP and low SSD and inverse mean power. Notable is the large negative freeboard occurring for sample 25. Figure 9 shows the waveforms in the vicinity of the negative freeboard and closer inspection of waveform 25 shows a multiple kinked leading edge towards the first peak that significantly delays the retracking point. The smoother sea ice areas, indicated by the darker tone in the SAR image, below the ridge begin to have an impact on the waveform and subsequent waveforms have an increased peak power. Observations over other ridges show that those features need to be sufficiently large, i.e., minimum 3 altimeter footprints, to not be affected by the more specular scattering of the FYI. The freeboard, expected to be larger over the ridge area, does, however, not show any significant changes compared to the surrounding FYI. Figure 10 shows the SAR image and the spatial variation of altimeter waveforms of an old MYI floe embedded in thinner FYI. The data was taken on 22 February 2018 with a time separation of about 20 min, i.e., for average drift speeds between 6 cm/s and 10 cm/s that means 72 m to 120 m. The floe size along-track is approximately 24 km. In the SAR image, the MYI floe is characterized by high backscatter in co-and cross-polarization (Figure 10a) shows co-polarization) while the surrounding FYI has a lower backscatter. In the upper part of the MYI floe, some lead and low backscatter features are visible that are also reflected by reduced mean backscatter intensity and an increased standard deviation (light blue area in Figure 10b,c). The lower part of the floe looks relatively homogeneous. The waveforms in the interior of the floe show a low peak power with an exception around sample number 40 (Figure 10d). In this area the peak power is increased compared to before and after, but no feature is visible in the SAR image. The altimeter is very sensitive to smooth areas within the footprint that are too small to be resolved by the SAR but affect the measured waveform. Over the lower part of the MYI floe the freeboard estimate reflects the higher thickness of the MYI, but in the upper part of the floe, where smoother sea ice features affect the waveforms, freeboard is reduced. Therefore, in areas of mixed ice types sea ice freeboard more likely represents the smoother and thinner part of the ice. Negative freeboard is also observed in areas of changing sea ice surface at the beginning and end of the MYI floe. This example demonstrates the difficulties of obtaining robust freeboard even over relatively large and apparently homogeneous features, but careful observation of the changes in waveform shape facilitates explanation of the observations.  Figure 11 shows the SAR image and waveform spatial evolution over a small MYI floe that is split in two by a band of thinner ice with high SAR backscatter and smooth texture, i.e., brash ice or thin ice with a rough surface or frost flowers atop. The data was taken on 2018-02-24 with a time separation of about 10 min, i.e., for average drift speeds between 6 cm/s and 10 cm/s that means 36 m to 60 m. In the SAR image (Figure 11a) only the spatial context separates the MYI floe from the high backscatter FYI. A subtle difference between co-and cross-polarization channel mean backscatter can be observed. In the first one the backscatter around the feature is slightly increased, while for the latter one the backscatter is slightly decreased (Figure 11b,c around sample 20). The waveform peak power decreases with the start of the MYI floe but increases significantly over the thinner FYI area (Figure 11d)). That means that the FYI area must include areas of smooth ice, e.g., low density layer of frost flowers or open water between small flows that are responsible for the high power return. PP and SSD show only slight variations between the FYI before the floe and the MYI floe itself, while the difference in IMP is more pronounced (Figure 11e,f). The freeboard at the first part of the floe is mostly negative and only the smaller second part shows an increase of freeboard compared to the surrounding FYI.
The three presented test cases show that altimetric waveforms respond well to changes of the sea ice surface in the order of a few along-track altimeter footprints. Retrieval of a robust freeboard estimation of floes, however, requires significant larger features of a few kilometers to overcome uncertainties in estimating the surface elevation.

Discussion
Near-coincidental high-resolution SAR imagery significantly aids the interpretation of altimeter waveforms in particular for small-scale features in the sea ice cover. Waveform derived parameters such as PP, IMP and SSD can principally discern between FYI, MYI and leads. The sensitivity of altimeters waveforms to smooth areas within the footprint, however, complicates the classification task. We introduced the IMP parameter to improve the separability between FYI and MYI for cases where the peak power of the two ice types are similar (see Figure 4). Large temporal and spatial variability of these parameters, however, impedes the establishment of universally valid thresholds. Spatio-temporal changes can have a multitude of sources, e.g., differences in snow depth and composition [30,31], changes in surface roughness in the order of the radar wavelength [24] or temperature variations [32]. To circumvent this problem, adaptive thresholding has been proposed for larger scale assessments of ice type variability from altimeter data [7]. For small-scale investigations, the relative spatial evolution of waveforms corresponds well to changes of sea ice surface types as determined from SAR imagery (see . Therefore, usage of SAR imagery offers possibilities to assess the reasons of sea ice type misclassification from CryoSat-2 data [7,8] and inconsistencies in the spatial distribution of waveform parameters and sea ice type products [6] left unanswered in previous studies. From the comparison of the altimeter data with the SAR images we observed that in areas of predominantly FYI, waveforms with high PP and peak power are common that are not linked to visible features in the SAR data. The increased peak power can be attributed to sub-resolution specular reflectors, e.g., thin smooth ice and leads that dominate the waveform. These findings are in line with observations from conventional (non-SAR) altimetry over the marginal ice zone [33]. Similarly, we found that in mixtures of FYI and MYI, waveforms tend to be more similar to FYI waveforms. The return from the smoother components are still strong enough to dominate the waveform in those cases. Difficulties in waveform interpretation over non-homogeneous surfaces are in accordance with earlier observations from Seasat data [15,34].
The complexity of altimeter waveform interpretation of mixed surfaces also showed an impact on freeboard estimates. Close to boundaries of distinct MYI floes in FYI, mixtures of ice types within a floe or small floes of MYI in FYI, freeboard estimates are lower compared to the interior of a homogeneous floe and similar to FYI freeboard measurements. While freeboard may vary over such a floe, this behavior is observed for nearly all floes in our study. Thus, the measurements most likely reflect the range to the thinner component either at nadir or off-nadir positions [35]. The latter ones often cause highly negative freeboard values because the range estimate is no longer associated with the nadir return. Underestimation of freeboard for mixed ice cases has been demonstrated in a simulation study based on CryoSat-2 parameters [36]. To reduce biases in freeboard retrievals due to ambiguous waveforms, returns with increased PP and peak power over FYI or waveforms with widened leading edge are usually discarded for sea ice thickness retrievals in standard gridded products [3,37]. Filtering of waveforms and spatial averaging reduces the effect of outliers but especially in the transition from FYI and MYI more research is necessary to quantify those effects on sea ice thickness retrievals. This is of great importance to reduce uncertainties in sea ice thickness estimates. While ambiguous waveforms introduce errors to freeboard estimates, local changes of sea ice characteristics leading to such irregular waveforms contain information about small-scale features of the sea ice surface.
Leads of thin ice or open water even smaller than the resolution of the SAR imagery are easily discerned by the altimeter due to their high PP and peak power. Caution needs to be taken to separate real leads and lead-like waveforms that are characterized by lower peak powers [38]. Lead detection is a necessary step in all sea ice freeboard retrievals and crucial to determine an optimal local sea surface height as a reference for the freeboard [39]. FYI within MYI floes most likely originates from refrozen leads or cracks in the ice and even small features are easily detected by waveform changes towards peakier more powerful returns. Freeboard estimates can be biased by off-nadir impact of the smoother reflectors [29,35]. Small-scale features with an increase in backscatter intensity in the SAR imagery, i.e., large ridges and old ice floes, are characterized by a decrease in PP and peak power in the altimeter waveform sequence. Features a few hundred meters in size can usually be identified but can only be unambiguously assigned to the correct cause by the SAR imagery. Freeboard measurements are, however, not significantly impacted by these small features compared to the surrounding FYI because of the reasons explained above for mixed ice cases. Acquisition of robust freeboard estimates over such features therefore requires this feature on the one hand to be large enough that the altimetric footprint is not affected by surrounding areas and on the other hand to provide enough valid freeboard measurements to obtain a meaningful average value. Based on the analysis of many MYI floes in our data set, we suggest that a floe size in the order of 10 km in along-track and a few kilometers to both sides across-track with a minimum of 10 valid freeboard measurements should act as a guideline for robust estimates of floe freeboard and hence ice thickness.
Altimetric waveforms can reflect small-scale changes of the sea ice surface in size scales from a fraction of the altimeter footprint for leads to a few hundred meters for ridges and old ice floes. The complexity of the waveform and the leading edge in particular introduced by non-homogeneous altimeter footprints, however, impede reliable freeboard estimates on these small scales. SAR imagery as additional data aids the interpretation of the altimeter data and facilitates the understanding of waveform responses to surface changes.

Conclusions
Using near-coincidental SAR imagery and altimetric data over three different winters in the Beaufort Sea, we showed that small-scale changes of a few hundred meters in sea ice surface related to leads, large ridges or ice floes of different sea ice types are recognizable by concurrent changes of altimetric waveforms. SAR imagery is a great asset to unambiguously identify the cause of these changes and to better understand the waveform behavior. SAR images thereby also facilitate the identification of regions potentially biased by the influence of thin ice features in older and thicker ice. Reliable freeboard estimates are possible for MYI floes in the order of 10 km in along-track and a few kilometers to both sides of the altimeter track. Freeboard estimates could be locally expanded to areas of similar signatures in the SAR image. The combination of information from SAR imagery and radar altimetry therefore has the potential to complete the picture of the sea ice situation on a local scale. SAR imagery however cannot solve the challenges imposed on freeboard to thickness retrieval by the properties of the snow layer [26], but improve the discrimination of ice types necessary for the calculation.
Retracking methods and selection procedures of valid waveforms for freeboard estimation have been optimized for hemispherical and basin-scale sea ice thickness retrieval [3,26]. Conversion of the here presented findings to useful additions to sea ice information from SAR imagery would benefit from the development of methods for locally optimized retrackers and improved waveform interpolation to better identify kinks in the leading edge of waveforms.