Scattering Characteristics of X-, Cand L-Band PolSAR Data Examined for the Tundra Environment of the Tuktoyaktuk Peninsula , Canada

In this study, polarimetric Synthetic Aperture Radar (PolSAR) data at X-, Cand L-Bands, acquired by the satellites: TerraSAR-X (2011), Radarsat-2 (2011), ALOS (2010) and ALOS-2 (2016), were used to characterize the tundra land cover of a test site located close to the town of Tuktoyaktuk, NWT, Canada. Using available in situ ground data collected in 2010 and 2012, we investigate PolSAR scattering characteristics of common tundra land cover classes at X-, Cand L-Bands. Several decomposition features of quad-, co-, and cross-polarized data were compared, the correlation between them was investigated, and the class separability offered by their different feature spaces was analyzed. Certain PolSAR features at each wavelength were sensitive to the land cover and exhibited distinct scattering characteristics. Use of shorter wavelength imagery (X and C) was beneficial for the characterization of wetland and tundra vegetation, while L-Band data highlighted differences of the bare ground classes better. The Kennaugh Matrix decomposition applied in this study provided a unified framework to store, process, and analyze all data consistently, and the matrix offered a favorable feature space for class separation. Of all elements of the quad-polarized Kennaugh Matrix, the intensity based elements K0, K1, K2, K3 and K4 were found to be most valuable for class discrimination. These elements contributed to better class separation as indicated by an increase of the separability metrics squared Jefferys Matusita Distance and Transformed Divergence. The increase in separability was up to 57% for Radarsat-2 and up to 18% for ALOS-2 data.


Introduction
Polarimetric Synthetic Aperture Radar (PolSAR) data from an increasing number of different satellite systems has become available-or will become available in the near future-for up-to-date Earth observation and environmental monitoring.Microwave data, e.g., acquired by Sentinel-1, ALOS-2, or in the future, by the RADARSAT Constellation Mission, are capable of delivering remote sensing data at high spatial (<10 m) and temporal resolutions (<10 days); independent of weather and illumination conditions.Therefore they are well suited for characterizing and monitoring the dynamic nature of the land surface, especially in vast and remote regions like the Arctic.In light of the challenges associated with a changing climate and environment, such investigations are crucial for assessing and comprehending temporal and spatial changes.Several studies have therefore investigated the use of SAR, PolSAR, and SAR interferometry (InSAR) for characterizing Arctic land surfaces, particularly that of tundra environments.
Table 1 provides an overview of select studies that involved the use of active microwave imaging of Arctic tundra environments.As can be observed, shortwave C-and X-Band data were frequently used for the characterization of land cover and shoreline types, and most studies incorporated analysis of polarimetric information to relate observed values with certain types or states of the land surface [1][2][3][4][5][6][7][8][9][10][11][12][13][14][15][16].C-Band and L-Band data stacks have been used primarily to determine surface movements in permafrost regions using interferometry [17][18][19][20][21].This is likely driven by the opening of the ALOS archive and the capacity of the L-Band microwaves to penetrate the relatively small tundra vegetation, typically resulting in higher interferometric coherence, and leading to a more reliable estimate of surface movements.The seasonal thawing and freezing of the active layer is also a promising area for InSAR applications, and estimation of the active layer thickness/variations across the entirety of the Arctic is of particular importance considering recent warming trends [22][23][24].
In light of this previous research, and our preliminary investigations [6,7,11,13,14], we analyze X-, C-and L-Band PolSAR data in order to characterize scattering properties of select tundra land cover classes for a test site in the Arctic.In addition, this study incorporates quad-polarized data of ALOS-2, a novel compilation of in situ data for the test site, and a complete utilization of the Kennaugh Matrix approach, recently presented by [31].The Kennaugh Matrix approach offers a unified framework for processing polarimetric information of different polarization modes (quad-, dual-and compact-polarized data).It can be used to represent targets both incoherently and coherently, and can be converted into all of the well-established decomposition models, for all wavelengths [25].Information on the Kennaugh Matrix framework is provided in the subsequent Section 2.2.2 Polarimetric SAR Data and Decompositions of this manuscript and in Appendix A.
The objectives of this research are therefore: First, process and analyze decomposition features of quad-and dual-polarized data of different sensors at three different frequencies.Second, investigate the backscattering of generalized tundra land cover classes for quad-and dual-polarized data of X-, C-and L-Band data.Third, investigate the correlation among PolSAR features of quad-and dual-polarized decomposition techniques.Fourth, benchmark and rank all PolSAR (decomposition) features in terms of class separability, and identify feature spaces and parameters that are most meaningful for characterizing the tundra land cover.
This manuscript is structured as follows: The subsequent section provides details on the materials and methods, as well as information on the location and environment of the test site selected for this research.Further, this section describes and lists the available land cover reference information, and the PolSAR data, including: quad-and dual-polarized data of Radarsat-2 (R-2), TerraSAR-X (TSX), ALOS and ALOS-2.Subsequently, the data processing and all polarimetric decompositions applied to the data are described, as well as separability measures/metrics used to analyse them, including: Transformed Divergence and Jefferys Matusita Distance.
The third section presents the results of the correlation analyses of C-Band R-2 and L-Band ALOS-2 data.Backscatter characteristics and select decomposition features of the land cover classes are presented and analysed via box-and scatterplots for X-, C-and L-Band data.Afterward, the results of the separability analysis and the feature selection are presented.Section four discusses the main findings, while section five provides a summary of the study; major conclusions are drawn and an outlook on future work is given.

Test Site Description
The study area is located at the northern extent of the mainland of the Northwest Territories, Canada (see Figure 1a).The region is part of the Canadian Arctic and lies entirely north of the tree-line along the coastal tundra plains of the Southern Arctic Ecozone [32].The climate here is characterized cold winters, followed by short and cool summers.The mean annual air temperature at the climate station Tuktoyaktuk is −1 • C (1971-2000), and the mean air temperature between October and April is below −10 • C. The average precipitation is about 150 mm [33].The ground surface is characterized by the presence of continuous permafrost and its thickness is estimated to be up to 600 m [34].Therefore, the soils of the Tuktoyaktuk Peninsula are Cryosols, with an active layer thickness of several centimeters to decimeters.The soils developed on glacial deposits of Pleistocene to Holocene age [32], and current morphodynamics are dominated by periglacial processes.The land surface is therefore characterized by low-lying and flat coastal plains, rolling hills, thermokarst lakes and pingos, and extensive networks of high-and low-centered ice wedges [33].The ground surface is also characterized by tundra vegetation, with upland tundra usually composed of short herbaceous vegetation and shrubs (dwarf shrubs up to tall shrubs).The wetland vegetation (grasses, sedges or rushes) is frequently at or near water bodies, e.g., at drained lakes, or in the flat and low-lying intertidal zone.Depending on the coastal currents, the beach zone is characterized by fine sandy material, mixed sediments dominated by gravel, pebble or cobble and driftwood accumulations [33].
to Holocene age [32], and current morphodynamics are dominated by periglacial processes.The land surface is therefore characterized by low-lying and flat coastal plains, rolling hills, thermokarst lakes and pingos, and extensive networks of high-and low-centered ice wedges [33].The ground surface is also characterized by tundra vegetation, with upland tundra usually composed of short herbaceous vegetation and shrubs (dwarf shrubs up to tall shrubs).The wetland vegetation (grasses, sedges or rushes) is frequently at or near water bodies, e.g., at drained lakes, or in the flat and low-lying intertidal zone.Depending on the coastal currents, the beach zone is characterized by fine sandy material, mixed sediments dominated by gravel, pebble or cobble and driftwood accumulations [33].

Land Cover Reference
In situ data on the land cover of the test site was acquired during two ground truth campaigns in the summer months of 2010 and 2012.The field work was organized and conducted by Carleton University, (Ottawa, ON, Canada), the NWRC (National Wildlife Research Centre, Ottawa, ON, Canada), and the University of Wuerzburg (Institute of Geography and Geology, Wuerzburg, Germany).Combining the land cover information of our preliminary studies [6,7,11,14], the land cover reference was merged to a common database showing the eight land cover classes listed and described in Table 2.
During the field campaigns information on the land cover of representative homogenous locations was recorded, categorized and mapped.The in situ classification of the land cover was completed by field experts, and within the frameworks of [35,36].As specified in Table 2 the cut off criteria for the separation of Shrub (ST) and Herb Dominated Tundra (HT) was the presence of dwarf shrubs with a height greater or less than 0.25 m.The separation between Sand (BS) and Mixed Sediment (BM) was based on the presence of pebble, which had to cover >50% of the surface; the bare ground samples were considered homogeneous if 10% or less were "other" materials or vegetation.Locations were classified as Driftwood Accumulation (BW) if more than 80% of the ground were cover by driftwood.Wetland Vegetation Communities (WT) were dominated by grasses, sedges or rushes and frequently occurred inland at creeks and drained lakes.The locations classified as Inundated Low Lying Tundra (WI) exhibited tundra vegetation communities at or near a water body and were most prominent in the low-lying coastal supratidal north of the town of Tuktoyaktuk.Further, Figure 2 provides example photographs of select land cover classes.In total, information from more than fifty ground truth sites were available.Additionally, the number of samples was increased using high resolution airborne imagery with less than one meter spatial resolution provided by [34].The generation of the land cover reference database was completed and locations of homogenous land cover were digitized using the airborne imagery.The reference information was then available in polygon format.Afterward, random sampling was applied, in order to generate 200 samples for each of the land cover classes listed in Table 2.Note that each individual point was selected to represent homogenous information for a certain land cover class (i.e., areas of mixed land covers were avoided).Figure 1a shows the locations of some sites visited and indicates the centers of the manually digitized polygons that exhibited homogeneous land coverage in the airborne imagery.
of homogenous land cover were digitized using the airborne imagery.The reference information was then available in polygon format.Afterward, random sampling was applied, in order to generate 200 samples for each of the land cover classes listed in Table 2.Note that each individual point was selected to represent homogenous information for a certain land cover class (i.e., areas of mixed land covers were avoided).Figure 1a shows the locations of some sites visited and indicates the centers of the manually digitized polygons that exhibited homogeneous land coverage in the airborne imagery.Additionally, a second set of 50,000 land samples was randomly generated, representing approximately 10% of all land pixels inside the common coverage of TSX, ALOS, ALOS-2 and R-2.This second set was generated independent of the land cover and was used to estimate the correlations among the PolSAR features.It is assumed that this sample represents the natural distribution of the relevant land cover classes.

Polarimetric SAR Data and Decompositions
PolSAR data from R-2, TSX, ALOS and ALOS-2 were available for the test site.Table 3a lists the main acquisition parameters and shows that all data was acquired in the summer months, during the growing season of the tundra vegetation.Note that the ALOS-2 data was acquired at a steep incidence angle (28°) in 2016, while the data of the other sensors were acquired with incidence angles between 34° and 40° in 2010/2011; along with the in situ reference data.Changes of the land cover were considered to be of less relevance for the analysis, considering the recent studies on the decadal changes in composition of the tundra vegetation here [37,38] the spatial resolution of the data, and the rather broadly defined classes.Figure 1b shows a false-color composite of the HV intensities of TerraSAR-X (TSX), R-2 and ALOS.Additionally, a second set of 50,000 land samples was randomly generated, representing approximately 10% of all land pixels inside the common coverage of TSX, ALOS, ALOS-2 and R-2.This second set was generated independent of the land cover and was used to estimate the correlations among the PolSAR features.It is assumed that this sample represents the natural distribution of the relevant land cover classes.

Polarimetric SAR Data and Decompositions
PolSAR data from R-2, TSX, ALOS and ALOS-2 were available for the test site.Table 3a lists the main acquisition parameters and shows that all data was acquired in the summer months, during the growing season of the tundra vegetation.Note that the ALOS-2 data was acquired at a steep incidence angle (28 • ) in 2016, while the data of the other sensors were acquired with incidence angles between 34 • and 40 • in 2010/2011; along with the in situ reference data.Changes of the land cover were considered to be of less relevance for the analysis, considering the recent studies on the decadal changes in composition of the tundra vegetation here [37,38] the spatial resolution of the data, and the rather broadly defined classes.Figure 1b shows a false-color composite of the HV intensities of TerraSAR-X (TSX), R-2 and ALOS.The following processing steps were applied to the PolSAR data: First, synthetic dual-polarized data (HH/VV, VV/VH and HH/HV) were generated from the R-2 and the ALOS-2 quad-polarized data.These synthetic dual-polarized datasets were thus not affected by temporal variations and showed identical speckle characteristics (compared to the quad-polarized data from which the subsets were taken), as such, this allowed for direct comparison of class separability as a function of polarization diversity, as opposed to differences in moisture, and plant phenology.Second, the Sinclair scattering matrices of all dual-and quad-polarized data of all wavelengths were converted to the corresponding Kennaugh Matrices [31].Third, the data were multi-looked (minimum of four looks) in order to generate pixels with square ground range resolution.Forth, a simple boxcar filter with a window size of 3 × 3 pixels was applied.Fifth, the data were terrain corrected and geocoded using the Range-Doppler Approach [39].All data were transformed to UTM WGS1984 Zone 8 coordinate system with 12 m spatial resolution using the TanDEM-X intermediate digital elevation model (DEM) and the projected local incidence angle derived from this DEM [40].The data were processed as sigma nought intensities.
All of the preceding steps were completed in SNAP 5.0 (Sentinel Application Platform) released by the European Space Agency (ESA), Paris, France.The terrain corrected Kennaugh Matrices were then used to generate the polarimetric channels, the Yamaguchi Decomposition [27], the Eigen-decomposition with the features Entropy/Alpha/Anisotropy, the Two Component Ground-Volume Decomposition of [30] and the Two Component Surface-Diherdal Decomposition of [14] using IDL 8.5 and ENVI 5.3.All intensity features were scaled to decibels [dB].The above mentioned decompositions are explained in more detail in the subsequent paragraphs.
Kennaugh Matrix-For quad-polarized data the Kennaugh Matrix (Mueller Matrix, respectively [41]) describes the relation between the radiated and received wave as a symmetric 4 × 4 matrix using ten real elements (K0-K9).It is the linear transformation of the four-dimensional Stokes vector ( [42] p. 43 ff.and p. 83 ff.) in the backscatter-alignment coordinate system.Unlike the Covariance or the Coherency Matrix, the Kennaugh Matrix can describe both coherent and incoherent targets [42,43].The elements K0, K1 and K2 are intensity-based elements, while K3 and K4 are based on the cross-polarized intensity and the co-polarized phase information.The elements K5, K6, K7, K8 and K9 are phase-only elements that tend to provide unique information from natural targets.All elements of the full Kennaugh Matrix can be grouped as follows [31]: First, the total intensity (K0); second, the absorption elements that describe the loss of polarization during the scattering process (K1, K2, K3); third, diattenuation elements that describe the change of the relation between two amplitude values during reflection (K4, K5, K6); fourth, retardance elements that describe the phase delay during scattering in a particular direction (K7, K8, K9).The definition of the Kennaugh Matrix and its elements for quad-polarized (A1), HH/VV-polarized (A2) and HH/HV-or VV/VH-polarized data (A3) are shown in Appendix A in accordance to [31].The Kennaugh Matrix elements are linear combinations of the Coherency Matrix and combinations of K0, K1, K2 and K3 describe the diagonal elements of the Coherency Matrix (T 11 , T 22 , T 33 ), while combinations of K4, K5, K6, K7, K8 and K9 describe off-diagonal elements of the Coherency Matrix (T 12 , T 13 , T 21 , T 23 , T 31 , T 32 ) [42].The conversions of the Kennaugh Matrix to 3 × 3 Coherency Matrix (T) of quad-polarized (A4) and 2 × 2 T of HH/VV-polarized data (A5) are shown in Appendix A in accordance to [31].To generate all Kennaugh matrix elements requires quadrature polarized data, thus only a portion of can be generated using dual polarized data.For HH/VV-polarized data the Kennaugh Matrix consists of the elements K0, K3, K4 and K7.For HH/HV-and VV/VH-polarized data the Kennaugh Matrix consists of the elements K0, K1, K5 and K6 [31].Figure 1c shows as a false-color composite of the Kennaugh Matrix elements K1, K2 and K3, which were processed using the quad-polarized R-2 data.For the purpose of comparison Figure 1d shows a Landsat TM true-color RGB composite acquired in summer 2011, concurrent with the R-2 imagery.
Eigen-decomposition-The Eigen-decomposition approach is a frequently used to process PolSAR data [25,26].It decomposes the incoherent signal (usually stored in the Covariance or Coherency Matrix) using eigenvalues (λ) and eigenvectors (u x ) (( 1) and ( 2)).In the formula H denotes the conjugate transpose.Note that the eigenvalues of the Covariance or Coherency Matrix are the same, while the eigenvectors differ.For dual-polarized data two eigenvalues and eigenvectors are obtained (1); while for quad-polarized data three eigenvalues and eigenvectors are obtained (2) when reciprocity is anticipated due to a monostatic acquisition geometry.Consequently, identical scattering from HV and VH is assumed.
Entropy/Alpha/Anisotropy-The Eigen-decomposition was used to process additional features that describe scattering processes [25,26,42,43].The polarimetric Entropy and Alpha Scattering angle describes the scattering properties of incoherent (natural) scatterers.Entropy ( 3) and ( 4) can be understood as the degree of randomness of the scattered signal and is described by the logarithmic sum of the pseudo probabilities p of the eigenvalues, and ranges from zero to one.The polarimetric Alpha scattering angle is calculated as the sum of the inverse cosine of the absolute value of the first eigenvector element and is weighted by the pseudo probabilities p (5). Cloude and Pottier also showed a third feature for quad-polarized data that is calculated via the ratio between the normalized difference of the second and third eigenvalue: the Anisotropy ( 6), which indicates the relevance of secondary scattering processes.Anisotropy, understood in the quad-polarimetric sense, is unavailable for dual-polarized data.In the formulas n is equal to two for dual-polarized data and three for quad-polarized data.Note that Entropies of the Covariance or Coherency Matrix are the same, but the Alpha scattering angles are different due to the differences between the Eigenvectors.
Model-Based Decompositions-Besides the Kennaugh Matrix elements, the polarimetric intensities and the Eigen-decomposition features, three Model-based decompositions were applied to the data which apply simplified, pre-defined scattering models.For the quad-polarized data of R-2 and ALOS-2, the Three Component Yamaguchi Decomposition [27] was applied.This approach decomposes the total backscattered energy P Total into the intensities of surface scattering (P surface ), double bounce scattering (P double bounce ) and volume scattering (P volume ) (7).This frequently used approach is suitable for comprehending and characterizing predominant scattering processes in nature.
P Total = P surface + P double bounce + P volume (7) As shown by [14] the approach of Yamaguchi can be adopted for HH/VV-polarized data, by decomposing the total backscattered energy P Total into the intensities of surface scattering (P surface ) and double bounce scattering (P double bounce ) (8).The correlation between the corresponding features of this decomposition and the Yamaguchi Decomposition are then a function of the presence and power of volume scattering processes [14].Specifically, features are more highly correlated if volume scattering is negligible.
P Total = P surface + P double bounce (8) For HH/VV-polarized data the approach of [30] can be applied as an alternative dual-polarimetric decomposition technique.The approach involves a synthetized HV channel and the polarimetric H (see the preceding paragraph) (see also [44]).This technique decomposes the total backscattered energy P Total into contributions from scattering from ground (P ground ) and from vegetation (P volume ) (9).
The two component decompositions of [14,30] were applied to the X-Band HH/VV data of TSX, to the synthetic HH/VV data of R-2 and ALOS-2.Table 3c lists all the polarimetric data that were used in this study, and provides abbreviations that are used hereafter to refer to each decomposition element.The descriptions of Figures 3 and 4 list all features that were processed for a certain type of polarized data, e.g., for HH/HV or HH/VV data.
The two component decompositions of [14,30] were applied to the X-Band HH/VV data of TSX, to the synthetic HH/VV data of R-2 and ALOS-2.Table 3c lists all the polarimetric data that were used in this study, and provides abbreviations that are used hereafter to refer to each decomposition element.The descriptions of Figures 3 and 4 list all features that were processed for a certain type of polarized data, e.g., for HH/HV or HH/VV data. ) ranging from 0.0 (no correlation) to 1.0 (full linear correlation, determination respectively).R 2 was estimated using 50,000 randomly distributed samples on land (roughly 10% of all land pixels).Note that dual-polarimetric data of Radarsat-2 were derived as polarimetric subsets and thus are not affected by temporal variations.Feature abbreviations are as follows (see Table   [31]), VOL2 (volume scattering of [30]), GRD2 (ground scattering of [30]), DBL2 (double bounce of the [14]), ODD2 (surface scattering of the [14]).) ranging from 0.0 (no correlation) to 1.0 (full linear correlation, determination respectively).R 2 was estimated using 50,000 randomly distributed samples on land (roughly 10% of all land pixels).Note that dual-polarimetric data of Radarsat-2 were derived as polarimetric subsets and thus are not affected by temporal variations.Feature abbreviations are as follows (see Table   [31]), VOL2 (volume scattering of [30]), GRD2 (ground scattering of [30]), DBL2 (double bounce of the [14]), ODD2 (surface scattering of the [14]).) ranging from 0.0 (no correlation) to 1.0 (full linear correlation, determination respectively).R 2 was estimated using 50,000 randomly distributed samples on land (roughly 10% of all land pixels).Note that dual-polarimetric data of ALOS-2 were derived as polarimetric subsets and thus are not affected by temporal variations.Feature abbreviations are as follows (see Table   [31]), VOL2 (volume scattering of [30]), GRD2 (ground scattering of [30]), DBL2 (double bounce of [14]), ODD2 (surface scattering of [14]).

Correlation, Class Separability, and Feature Selection
The correlations between the above listed decomposition features were examined using the dB-scaled sigma nought intensity values of the calibrated data, where applicable.Correlations were estimated using a random sample of 50,000 points over land (see Section 2.2.1.Land Cover Reference); thus values and analyses were completed independent of the land cover classes of interest.The squared linear Pearson Correlation Coefficient (R 2 ) was used in all cases.The coefficient ² is defined as the squared ratio between the covariance (Cov) of two variables (i;j) and the product of the ) ranging from 0.0 (no correlation) to 1.0 (full linear correlation, determination respectively).R 2 was estimated using 50,000 randomly distributed samples on land (roughly 10% of all land pixels).Note that dual-polarimetric data of ALOS-2 were derived as polarimetric subsets and thus are not affected by temporal variations.Feature abbreviations are as follows (see Table   [31]), VOL2 (volume scattering of [30]), GRD2 (ground scattering of [30]), DBL2 (double bounce of [14]), ODD2 (surface scattering of [14]).

Correlation, Class Separability, and Feature Selection
The correlations between the above listed decomposition features were examined using the dB-scaled sigma nought intensity values of the calibrated data, where applicable.Correlations were estimated using a random sample of 50,000 points over land (see Section 2.2.1.Land Cover Reference); thus values and analyses were completed independent of the land cover classes of interest.The squared linear Pearson Correlation Coefficient (R 2 ) was used in all cases.The coefficient R 2 is defined as the squared ratio between the covariance (Cov) of two variables (i;j) and the product of the individual standard deviations (σ i σ j ) (10).R 2 is frequently used to quantify the degree of determination between two variables, though still can be interpreted as a coefficient that quantifies the correlation.R 2 ranges from zero to one; a value of one (zero) indicates perfect (no) linear correlation and a maximum (minimum) determination, 100% (0%) of the explained variance, respectively [45].
Cov(i, j) All PolSAR data were then used in separability analyses to quantify the ability of the polarimetric information to discriminate the land cover classes considered in this research.The Transformed Divergence (TD) (11) [45,46], Bhattacharyya Distance (BD) (12) [47], and Jefferys Matusita Distance (JD) (13) [45,46] were processed for each PolSAR feature space and each wavelength for all land cover classes [47].The features are processed for two classes c and d by assessing the classes' mean vectors M (14) and the classes' covariance matrix V (15) for a given set of features (as a minimum, two features are required).In the formula tr denotes the trace of a matrix, formula det denotes the determinant of a matrix, T refers to the matrix/vector transpose, and Cov denotes the covariance.The separability features TD and JD have been shown to act as meaningful predictors for classification potential, thus a high separability indicates greater potential for class discrimination [48,49].
The metrics TD and JD can further be used for feature selection in order to identify those that are most meaningful for class separation among a given set of features.This can be achieved by calculating the increase in separability (SI) (16): displayed as the amount of separability (SP) that is gained when a feature of interest (x), e.g., K0, is added to an existing feature space (K).The average increase in separability can be processed by averaging the SI values of each possible feature combination, e.g., the increase in separability when K0 is added to {K1, K2} or {K1, K3} or . . ., {K1, K2, K3} or {K1, K2, K4} or . . ., {K1, K2, K3, K4} or {K1, K2, K3, K5} or . . ., and so on.
The separability metrics were employed to demonstrate the differences between the PolSAR features, to gauge their use in classification, and to determine which land cover classes can be separated with the PolSAR features.All of the investigated separability distances require normally distributed data, or at least symmetrically distributed data.Such symmetric distribution properties can be assumed for most of the investigated features.

Corrleation
Correlations among the decomposition features of dual-and quad-polarized data were investigated prior to the assessment of the backscatter characteristics of the land cover classes and the separability of classes.The squared Pearson Correlation Coefficient (R 2 ) was derived using the 50,000 randomly distributed samples on land, and which represented 10% of all land pixels inside the common coverage of the TSX, ALOS, ALOS-2 and R-2 imagery.Results were drawn as correlation matrices for the features of C-Band data of R-2 in Figure 3, for features of L-Band data of ALOS-2 in Figure 4.
For C-Band data (Figure 3), it was observed that Kennaugh Matrix elements K4, K5, K6, K7, K8 and K9 of the quad-polarized data showed the lowest correlations among all other investigated decomposition features.An explanation for this is that the Kennaugh elements usually are uncorrelated, and that the elements K5 to K8 are phase-based elements, which are more or less stochastically fluctuating over natural targets [31].The R 2 values were less than 0.4 with the exception of K4 and K7.These features were highly correlated with K4 and K7 of HH/VV-polarized data (R 2 values of about 0.9).This high degree of correlation is because K4 holds the relation of HH to VV, which is not kept in other decompositions, and the same applies for K7.As well, the correlations between K5 and K5 of HH/HV and VV/VH were moderately high (R 2 values of about 0.6).Similarly, the R 2 values of K8 and K6 of HH/HV and VV/VH were around 0.5.The reason for these observations are most likely the similar polarimetric behavior (diagonal diattenuation [31]) only with different input channels (HH/HV and VV/VH, respectively).
The same observations were made for the L-Band data of ALOS-2 (Figure 4); however, correlations between the Kennaugh Matrix elements K4, K5, K6, K7, K8 and K9 of quad-, HH/VV-, HH/HV-and VV/VH-polarized data were generally higher.For example, correlation between K7 of quad-polarized Kennaugh Matrix and K3 of HH/VV-polarized Kennaugh Matrix showed R 2 values of about 0.5.Among the Kennaugh Matrix elements K0, K1, K2 and K3 of quad-, HH/VV-, HH/HV-and VV/VH-polarized data, the following distinct linear correlations were observed: K0 of quad-and HH/VV-polarized Kennaugh Matrix showed R 2 values greater than 0.9 in the C-and L-Band; K3 of quad-and HH/VV-polarized Kennaugh Matrix showed R 2 values greater than 0.8 in the C-and L-Band; K0 of quad-, HH/HV-and VV/VH-polarized Kennaugh Matrix showed R 2 values greater than 0.8 in the C-and L-Band.Again the correlations between the Kennaugh Matrix elements of quad-, HH/VV-, HH/HV-and VV/VH-polarized data were generally higher in the L-than in the C-Band.Thus, most likely the L-Band data is more "stable" in a polarimetric sense due to a longer wavelength.Further, the high correlation of K0 of quad-and HH/VV-polarized Kennaugh Matrix is present since HH and VV record the vast majority of backscatter, while the HV contribution is negligible.
Among the model-based (power) decomposition features of quad-and HH/VV-polarized data, good correspondence between the DBL3 and DBL2 (R 2 values of about 0.8), the ODD3 and ODD2 (R 2 values of about 0.7), the VOL3 and VOL2 (R 2 values of about 0.7 (L-Band) and 0.4 (C-Band)) and the ODD3 and GRD2 (R 2 values of about 0.8 (L-Band) and 0.6 (C-Band)) was observed for both C-and L-Band data.The R 2 values between any of the model-based (power) decomposition features and any other polarimetric feature were lower than these observations, with the exception of VOL3 and HV showing R 2 values of about 0.95 (C-and L-Band).The reason for this observation can be seen in the low proportion of volume scattering for the tundra environment, making the influence of the cross-polarization component negligible, and decomposition features of quad-and HH/VV-polarized data highly correlated.
With respect to the Eigen-decomposition features, ENT and polarimetric Alpha scattering angles (ALPT/ALPC) were highly correlated between the ENT and the ALPT of HH/VV-polarized data, and this was true for both C-and L-Band (R 2 values of about 0.7).Additionally, ENT and ALPT were highly correlated with each other, with R 2 values of about 0.7-0.8(quad-and HH/VV-polarized data) and ENT and ALPC were moderately correlated with R 2 values of about 0.4-0.5 (HH/HV-and VV/VH-polarized data).This is most likely because most reflection is recorded in HH and VV intensities, making their contributions higher than the intensities of HV or VH.This leads to a high correlation between the HH/VV-and quad-polarized decomposition features.
In summary, this assessment indicated that Kennaugh Matrix elements K0, K1, K3, K4 and K7 of quad-and HH/VV-polarized data of C-and L-Band were highly correlated and thus can be used interchangeably in some cases, e.g., for image classification.The Kennaugh Matrix elements K5 and K6 of quad-and cross-polarized data showed lower correlation coefficient values.The correlation was generally higher in the L-band, compared to the C-Band, which is likely a result of less interaction between the incident wave and the vegetation body of the long L-Band microwaves; less volume scattering occurs.

Backscatter Characterisics
Figure 5 provides boxplots-showing the minimum, lower quartile (25%), median (50%), upper quartile (75%), maximum-of the land cover classes (see Table 2) for select polarimetric features of X-Band (TSX), C-Band (R-2) and L-Band (ALOS and ALOS-2).Figure 5a-i display the backscatter characteristics of the land cover classes concerning the HH, HV and VV sigma nought intensities in decibels (dB).Figure 5j-o shows the information of the model-based (power) decomposition features of the Yamaguchi Decomposition of C-Band (R-2) and L-Band (ALOS-2) as DBLB3, VOL3 and ODD3; in dB. Figure 5p-x shows the boxplots of the land cover classes for the Kennaugh Matrix elements K0, K3 and K4 of HH/VV-polarized X-Band (TSX), quad-polarized C-Band (R-2) and quad-polarized L-Band (ALOS-2) data in dB.
The class OL showed a unique range of intensity values in the K0 and VV of X-and C-Band, the VOL3 and ODD3 of C-Band, and the K4 and HV of X-, C-and L-Band data.The scattering differences between water and land were clearly pronounced, as water was generally characterized by a low intensity value.This is because water surface was relatively calm, thus it was not observed as rough; the shallow angle incidence angles of the X-and C-Band data, and the longer wavelength of the L-Band.If the water surface were to become rough due to higher wind speeds, higher intensity values of K0 would be observed, thus complicating the separation of the classes.In such cases it is assumed that K3 and K4 will still be suitable to separate land from water, since both are indicators for double bounce scattering, typically minimal for water.Further, BS and BM showed increased HH, VV, HV, DBL3 and ODD3 scattering at C-and L-Band compared to other classes.The range of values observed for the BW class was more unique for L-Band features, than X-and C-Band features.The most distinct values were for L-Band VOL3; as the BW's median value exceeded +5 dB, compared to the median value of all other land cover classes.Scattering from BW at L-Band is therefore characterized by high intensity values for HH, VOL3 and ODD3, with medians of about −10 dB and −7 dB; whereas at C-Band the scattering from BW is characterized by high HH and ODD3 intensities.The median intensity of BW is comparably low in X-Band HH and differences of BW's statistics to BS and HT are less pronounced.Thus, independent of wavelength and decomposition technique, the data were sensitive to scattering differences between land and water (OL), and between sandy bare ground (BS) and mixed non-vegetated sediment (BM).BM had higher backscatter than BS, which can be attributed to the higher surface roughness of BM (grain sizes of 2.0 mm-256.0mm) compared to BS (grain sizes of 0.0625 mm-2.0 mm), which leads to higher backscatter intensities.The difference between HT and ST scattering is characterized by increased HV, VOL3, and K0 intensities at X-and C-Band; however, there is substantial overlap in their distributions and the differences between median values were small; ranging from +2 dB (X-Band HV) to +3 dB (C-Band VOL3) between HT and ST.The largest differences in HT and ST statistics were nevertheless found for VOL3 of the C-Band data, but the data ranges of HT and ST (lower quantile to upper quantile) also overlap the ranges of BM, BW, WI and WT.Therefore, the X-and C-Band showed higher volume scattering intensities from shrub dominated tundra (ST) compared to herb dominated tundra (HT).This is likely due to a higher proportion of volume scattering in the shrub plants, which is caused by the relatively short wavelength.Contrarily, the L-Band HV and volume scattering intensities (VOL3) were not sensitive to this difference.2).
The L-Band HV channel and the volume scattering intensity of the Yamaguchi decomposition showed the same median values and comparable data ranges (lower to upper quantile) for these two land cover classes.At L-Band, the difference between HT and ST was better expressed via the VV channel, and the surface scattering intensity of the Yamaguchi decomposition via lower intensities of the ST compared to HT.This indicates an absence of volume scattering processes and a full penetration of the vegetation by the L-Band microwaves.Assuming that this observation is not caused by temporal variations present in the ALOS-2 data-which were acquired in 2016, while C-, X-Band and the land cover reference data were acquired in 2010/2011, the signal can be interpreted to represent backscattering mostly from the ground, as it is assumed that major changes in land cover type present have not occurred in this time.
The differences between the land cover classes WI and WT were characterized by increased DBL3, K0, HH and VV scattering in X-and C-Band and by K3 in X-, C-and L-Band.The statistics of WT showed a clear separation from the other land cover classes in the DBL3 and K3 of C-Band and the HH and K3 of X-Band.The difference between the WT's median value and the median value of any other land cover class exceeded +5 dB in the X-Band HH and C-Band DBL3.The differences between wetland (WT) and inundated low-lying tundra (WI) was observed as higher HH and VV intensities-and the Kennaugh Matrix element K4 accordingly-in X-and C-Band and the double bounce intensity of the Yamaguchi decomposition.Further, both classes were characterized by comparably low values of the Kennaugh Matrix element K3, which points to distinct double bounce scattering, since S HH S * VV is a known discriminator for this type of scattering (compare [27,31]).In contrast, both classes showed low intensities in the HH, VV and double bounce of the Yamaguchi decomposition at L-Band.Accordingly, K3 and K4 were less distinct and no double bounce scattering was present, when using the L-Band, which again is most likely due to the relatively short statured vegetation, the high penetration depth, and the absence of interactions between incident microwaves and the water surface and vegetation canopy.
In addition to the boxplots, Figure 6 shows scatterplots of the Kennaugh Matrix elements K0, K3 and K4 of X-, C-and L-Band data in order to investigate the scattering characteristics of the land cover classes in a multivariate feature space.Figure 6a-i shows the position of the land cover reference in the K0/K3 (left column), K0/K4 (center column) and K3/K4 (right column) feature spaces of X-Band (a-c), C-Band (d-f) and L-Band (g-i).These results show that the feature space K0/K3 of X-and C-Band facilitates the differentiation of the classes: OL, BS and WT; however, the position of values for BM, HT, ST and WI were indiscriminant from others.The feature space K0/K4 (Figure 6b,e,h) shows increased distance between samples of HT/ST and WI-especially at X-Band; however, a substantial degree of overlap between the samples of BM and BW, and the samples of HT, ST and WI was present.K0 provided the best separation between land cover classes at X-Band, C-Band, and L-Band.This is unsurprising since K0 of X-and C-Band showed a high positive linear correlation (R 2 of 0.8) (Figure 6j).The combination of short-and longwave SAR facilitated the separation of the WI, BM and BW samples, and the features of X-/C-and L-Band showed no linear correlation (R 2 less than 0.1) (Figure 6k,l).In summary, the X-, C-and L-Band data exhibit distinct scattering characteristics for the different land cover classes.All PolSAR data were sensitive to the OL, BS and BM coverage; additionally, L-Band data were most sensitive to the BW.The X-and C-Band features were suited to pronounce differences in WI and WT, and HT and ST coverage via the features HV, VOL3, DBL3, K0 and K3.In summary, the X-, C-and L-Band data exhibit distinct scattering characteristics for the different land cover classes.All PolSAR data were sensitive to the OL, BS and BM coverage; additionally, L-Band data were most sensitive to the BW.The X-and C-Band features were suited to pronounce differences in WI and WT, and HT and ST coverage via the features HV, VOL3, DBL3, K0 and K3.

Class Separability and Feature Selection
Multivariate assessment was completed for all feature spaces of the decomposition elements of dual-and quad-polarized X-, C-and L-Band datasets.The separabilities between classes (for the feature spaces of interest) were quantified for all possible combinations of variables (41) using the Transformed Divergence (TD) (Table 4) and the squared Jefferys Matusita Distance (JD) (Table 5).The feature spaces were ranked in descending order based on the average separability (AV) by the feature space of interest.For the purpose of comparison the tables also list separabilities achieved with the multispectral Landsat TM data using the six spectral bands (thermal information was excluded).
Results showed that the quad-polarized Kennaugh Matrices of ALOS-2 and R-2 offered the best separation of all land cover classes, followed by the HH/VV-polarized Kennaugh Matrices of ALOS-2, R-2 and TSX.The separability distances TD and JD further indicated that AV of ALOS-2 was comparable to the AV offered by multispectral data.As well, among the different PolSAR decompositions, the use of all Kennaugh Matrix elements was more beneficial for class separation than using the features of the model-based (power) Decompositions, Eigen-decompositions, or the intensities of the polarimetric channels.For C-and L-Band it was further observed that TD and JD of the Kennaugh Matrix decreased from quad-, to HH/VV-, to VV/VH-to HH/HV-polarized data.
For X-Band it was observed that TD and JD of the Kennaugh Matrix decreased from HH/VV-to HH/HV-polarized data.The separability of the Eigen-decomposition features ENT, ALPT or ALPC and ANI was low, and these feature spaces, as indicated by JD, offered the lowest separability between classes among all investigated feature spaces; independent of the wavelength (X-, C-or L-Band).With HH/VV-polarized data, the high correlation of Entropy and the Alpha scattering angles was observed by others [50], though.Another reason for this might be the lack of diversity of scattering processes in this rather "bare" landscape, thus the Entropy/Alpha feature space remains "unfilled" to a certain degree.Specifically, the tundra landscape examined in this research, offers a minor depolarizing, and low entropy environment.
The separability distance JD further outlined that the average class separability decreased from ENT/ALPT/ANI (quad) to ENT/ALPT (HH/VV) to ENT/ALPC (VV/VH) to ENT/ALPC (HH/HV).This might simply be a function of intensity, which decreased from quad to HH/VV to VV/VH to HH/HV, since with lower intensities there is also lower information content.Among the land cover classes the classes OL, BS, BW and WT were shown to be the land cover classes with the highest average separability, thus the PolSAR data were especially suited to characterize these classes.The lowest average separability was observed for the land cover classes BM and WI, while separability of HT and ST was moderately high.
Table 6 draws the average increase of the separability features TD and JD for the ten elements of the quad-polarized Kennaugh Matrices (K0-K9) of C-Band (left column) and L-Band (right column) data.This assessment indicates that the average increase in separability when a feature of interest (K0-K9) is added to an existing feature space.This metric was used to identify the most important elements of the Kennaugh Matrix for class separation.For both C-and L-Band it was observed that K0, K1, K2, K3 and K4 were more important for class separation than K5, K6, K7, K8 and K9.An explanation for this observation is that K0, K1 and K2 are intensity-based elements and are thus value are generally stable.K3 and K4 also use the cross-polarized intensity and the co-polarized phase, and are therefore relatively stable as well.Contrarily, K5, K6, K7, K8 and K9 are phase-based elements and therefore generally unstable in natural environments.With respect to the Coherency Matrix, K0 to K4 explain the diagonal elements which can be associated with dominant scattering processes, including: surface, double bounce and volume scattering.Hence (K5-K9, representing the off-diagonal elements of the Coherency Matrix, have minor relevance for the class separability.Table 4. Average class separability measured as Transformed Divergence (TD) for classes OL (Permanent Water Bodies), BS (Sand), BM (Mixed Sediment), BW (Driftwood Accumulation), HT (Herb Dominated Tundra), ST (Shrub Dominated Tundra), WI (Inundated Low Lying Tundra) and WT (Wetland) (see Table 2).TD is ranging from 0 to 2000; higher values indicate better class separation.The feature spaces are ranked in descending order of the average separability (AV).The black bars are scaled linearly between the minimum and maximum of AV.The colors from red to yellow to green correspond to the 10%, 50% and 90% quantiles of the AV.Table 4. Average class separability measured as Transformed Divergence (TD) for classes OL (Permanent Water Bodies), BS (Sand), BM (Mixed Sediment), BW (Driftwood Accumulation), HT (Herb Dominated Tundra), ST (Shrub Dominated Tundra), WI (Inundated Low Lying Tundra) and WT (Wetland) (see Table 2).TD is ranging from 0 to 2000; higher values indicate better class separation.The feature spaces are ranked in descending order of the average separability (AV).The black bars are scaled linearly between the minimum and maximum of AV.The colors from red to yellow to green correspond to the 10%, 50% and 90% quantiles of the AV.   2).JD is ranging from 0 to 2; higher values indicate better class separation.The feature spaces are ranked in descending order of the average separability (AV).The black bars are scaled linearly between the minimum and maximum of AV.The colors from red to yellow to green correspond to the 10%, 50% and 90% quantiles of the AV.
Table 5.Average class separability measured as squared Jefferys Matusita Distance (JD) for classes OL (Permanent Water Bodies), BS (Sand), BM (Mixed Sediment), BW (Driftwood Accumulation), HT (Herb Dominated Tundra), ST (Shrub Dominated Tundra), WI (Inundated Low Lying Tundra) and WT (Wetland) (see Table 2).JD is ranging from 0 to 2; higher values indicate better class separation.The feature spaces are ranked in descending order of the average separability (AV).The black bars are scaled linearly between the minimum and maximum of AV.The colors from red to yellow to green correspond to the 10%, 50% and 90% quantiles of the AV.
Appl.Sci.2017, 7, 595 20 of 28 Table 5.Average class separability measured as squared Jefferys Matusita Distance (JD) for classes (Permanent Water Bodies), BS (Sand), BM (Mixed Sediment), BW (Driftwood Accumulation), HT (Herb Dominated Tundra), ST (Shrub Dominated Tundra), WI (Inundated Low Lying Tundra) and WT (Wetland) (see Table 2).JD is ranging from 0 to 2; higher values indicate better class separation.The feature spaces are ranked in descending order of the average separability (AV).The black bars are scaled linearly between the minimum and maximum of AV.The colors from red to yellow to green correspond to the 10%, 50% and 90% quantiles of the AV.For C-Band it was further observed that K0 and K1 offered the highest increase for the separation of ST, while for L-Band K1 and K3 were more important for the separation of this land cover class.At C-Band, the information of K2 and K3 was beneficial for the separation of BW, WI and WT.The land cover class HT was best characterized by the elements K0 and K2 at C-Band and K0, K1 and K3 at L-Band.This means that K0 (total backscattered intensity), K1 (absorption element showing the difference between co-and cross-polarized intensities), K2 and K3 (absorption elements that describe the loss of polarization during the scattering process) and K4 (diattenuation element showing the difference between HH and VV intensities) are good descriptors for the examined For C-Band it was further observed that K0 and K1 offered the highest increase for the separation of ST, while for L-Band K1 and K3 were more important for the separation of this land cover class.At C-Band, the information of K2 and K3 was beneficial for the separation of BW, WI and WT.The land cover class HT was best characterized by the elements K0 and K2 at C-Band and K0, K1 and K3 at L-Band.This means that K0 (total backscattered intensity), K1 (absorption element showing the difference between co-and cross-polarized intensities), K2 and K3 (absorption elements that describe the loss of polarization during the scattering process) and K4 (diattenuation element showing the difference between HH and VV intensities) are good descriptors for the examined tundra land cover classes.The elements K7, K8 and K9 (descriptors of the phase delay during the scattering in a certain For C-Band it was further observed that K0 and K1 offered the highest increase for the separation of ST, while for L-Band K1 and K3 were more important for the separation of this land cover class.At C-Band, the information of K2 and K3 was beneficial for the separation of BW, WI and WT.The land cover class HT was best characterized by the elements K0 and K2 at C-Band and K0, K1 and K3 at L-Band.This means that K0 (total backscattered intensity), K1 (absorption element showing the difference between co-and cross-polarized intensities), K2 and K3 (absorption elements that describe the loss of polarization during the scattering process) and K4 (diattenuation element showing the difference between HH and VV intensities) are good descriptors for the examined tundra land cover classes.The elements K7, K8 and K9 (descriptors of the phase delay during the scattering in a certain direction) play a minor role in the separation of classes, as phase delays happen during volume propagation.Since tundra vegetation has a relatively short stature (height), phase delays due to volume propagation are less likely.Table 6.Average increase of the separability features: (a-b) Transformed Divergence (TD) and (c-d) squared Jefferys Matusita Distance (JD) for Kennaugh-Matrix elements K0-K9 of quad-polarimetric C-Band Radarsat-2 (left column) and L-Band ALOS-2 (right column) data and for classes OL (Permanent Water Bodies), BS (Sand), BM (Mixed Sediment), BW (Driftwood Accumulation), HT (Herb Dominated Tundra), ST (Shrub Dominated Tundra), WI (Inundated Low Lying Tundra) and WT (Wetland) (see Table 2).The features are ranked in descending order of the average increase in separability (AV).The last column displays the AV in percent (%) The colors from red to yellow to green correspond to the 10%, 50% and 90% quantiles of the AV data range.The metric displays the average increase in separability when a feature of interest (K0-K9) is added to an existing feature space (see Section 2.3 Correlation, Class Separability and Feature Selection).
Appl.Sci.2017, 7, 595 21 of 28 tundra land cover classes.The elements K7, K8 and K9 (descriptors of the phase delay during the scattering in a certain direction) play a minor role in the separation of classes, as phase delays happen during volume propagation.Since tundra vegetation has a relatively short stature (height), phase delays due to volume propagation are less likely.for the tundra land cover investigated in this research.Contrary to this, the quad-polarized Kennaugh Matrix elements K5 and K6 were less correlated to the corresponding elements of the cross-polarized Kennaugh Matrix.Nevertheless, there are still benefits associated with the Kennaugh Matrix, since all kinds of PolSAR data can be stored, processed and analyzed in the same manner.It also provides a unified framework without any loss of information, and the capacity to interpret decomposed elements in a coherent and incoherent way since any other incoherent or coherent scattering matrix can be derived if necessary [31,42].The Stokes coordinate system used for the definition of the Kennaugh Matrix seems to offer an appropriate approach to characterize the environment investigated in this research.
For the examined tundra land cover of the Tuktoyaktuk Peninsula, it was further shown that the elements of the Two Component Decompositions of [14] and [30] were highly correlated with the corresponding elements of the Yamaguchi Decomposition and-with lower significance-volume scattering, and HV intensity, respectively.Thus, the HH/VV-polarized data provide crucial information for describing the land covers considered in this research.As pointed out by [14], the correspondence of these quad-and HH/VV-polarized decomposition features is a function of the presence and influence of volume scattering processes, relative to contributions from the ground.Thus, due to the relatively short stature (height) of tundra vegetation, the observed correlations were high due to a lack of a significant volume scattering component.Further, the correlations between the features were generally higher at L-Band, compared to features at C-Band.This can be attributed to the longer wavelength of the ALOS and ALOS-2 sensors, and the absence/weakness of random scattering processes as the penetration depth is higher and volume scattering is less likely (thus the volume component is small relative to surface scattering).
The backscatter characteristics of the tundra land cover classes were examined via box-and scatterplots of the individual PolSAR features.It was shown that X-, C-and L-Band data exhibit distinct scattering characteristics for the different land cover classes.Results indicate that the L-Band data were more sensitive to the bare ground classes; thus, it is better suited to investigate and monitor ground properties, e.g., soil moisture, or the surface heave and subsidence (via InSAR) caused by the freezing and thawing of the active layer (compare [17,20,21]); especially in sites dominated by shrubs.In contrast, use of short wavelengths (X-and C-Band) is beneficial for characterizing tundra and wetland vegetation.This observation is in accordance with other studies [9,12,15].
It is worth noting the clear distinction of the land cover class: driftwood accumulation (BW) in the L-Band data.The coverage of BW is characterized by non-vegetated, dead woody debris, and frequently such accumulations exhibit a very high surface roughness, since dead wood and stems pile up more than a meter high (compare Figure 2c).Even though this should be a clearly visible target, and distinct feature in the PolSAR data, the position of BW is less clear in the feature spaces of X-and C-Band compared to the position of BW in the L-Band feature space.For BW the highest HH and HV intensity values (derived from ALOS in 2010) were found among all land cover classes.As well, the scattering from this type of coverage was characterized by high volume intensities of the Yamaguchi decomposition (derived from ALOS-2 in 2016) at L-Band.The dielectric and geometric properties of the driftwood accumulations facilitate high intensity scattering at L-Band, thus this type of coverage is a "rough" target at L-band but not in C-and X-Band (i.e., because the logs are much larger than incident C-and X-Band microwaves).
Even though the L-Band data showed limited value for characterizing the land cover classes HT, ST and WT using a single feature, the ALOS-2 quad-polarized and HH/VV-polarized data offered the feature space with the highest class separability; as indicated by the Transformed Divergence (TD) and squared Jefferys Matusita Distance (JD).However, since the ALOS-2 data were acquired at a steeper incidence angle and with a delay of six years, a true comparability of these results cannot be guaranteed.These results are therefore surprising, since one would assume a change of the land cover over time and an increasing dissimilarity between the reference and the PolSAR measurement with increasing temporal difference.Still, the ALOS HH/HV-polarized data acquired in 2010 showed a fairly good separability (Rank 17 in JD, Rank 18 in TD), and the data were observed to be more valuable for class separation than the C-Band (Rank 23 in JD, Rank 26 in TD), or X-Band (Rank 20 in TD, Rank 16 in JD) HH/HV-polarized data.
All separability features indicated that the Kennaugh Matrix was the most favorable feature space among all examined decompositions, which is in accordance with the expectation that full PolSAR information is better suited for class separability than is available via Entropy/Alpha, or the Two/Three Component Decomposition models, for instance.Among the Model-based Decompositions, the Yamaguchi Decomposition of quad-polarized data exceeded the separability offered by the Two Component Decomposition models.Thus, cross-polarized information is important for class discrimination, even though volume scattering processes play a minor role for the tundra environment investigated.Thus, perhaps differences in roughness/geometry play a more important role.

Conclusions
Results from this analysis indicate that the quad-polarized Kennaugh Matrix elements K0, K1, K3, K4 and K7 were highly correlated with corresponding elements of the dual-polarized Kennaugh Matrices; therefore, to a certain extent, dual-polarized data provide a useful substitute for the full quad-polarized data.The Kennaugh Matrix offers a unified framework to store, process and analyze PolSAR data in the same manner, and the Kennaugh elements offer comparable information from dualor quad-polarized data.Thus, there is nearly no difference between the two acquisitions modes when using Kennaugh elements.
Among the investigated Model-based Decompositions and the Eigen-decompositions the features of the Two Component Decompositions models of [14] (based on HH/VV dual-polarized data) were highly correlated with the corresponding elements of the Yamaguchi Decomposition (based on quad-polarized data).Independent of the wavelength and polarization mode, the Eigen-decomposition features Entropy and the Alpha scattering angles were highly correlated and of less value for class separation.Therefore, this approach does not seem suitable for this low depolarizing as well as low entropy environment.
The X-, C-and L-Band data exhibit distinct scattering characteristics for the different land cover classes.The PolSAR data of all wavelengths are sensitive to the land cover classes: open water (OL), sand (BS) and mixed sediment (BM); L-Band data were most sensitive to the BW; X-and C-Band features were most sensitive to the inundated low-lying tundra (WI) and wetland WT, and herb dominated tundra (HT) and shrub dominated tundra (ST).The use of shorter wavelengths (X-and C-Band) is beneficial for characterizing wetland vegetation.The L-Band data exhibited the differences of the bare ground classes BS, BM and BW best.Thus, in accordance to previous studies L-Band data are favorable for InSAR applications in this region, due to the observed distinct surface scattering and the low volume scattering contribution.In contrast, C-and X-Band data are favorable for the characterization of the tundra land cover due to the observed sensitivity of the cross-and co-polarized information to tundra vegetation.
Nevertheless, the assessment of the class separability pointed out that PolSAR data of any wavelength-also of L-Band-were valuable for class separation and PolSAR information is beneficial for class discrimination.The results showed that quad-polarized data of ALOS-2 and R-2 offered the best separation of the land cover classes, followed by the HH/VV-, HH/HV-or VV/VH-polarized data of ALOS-2, R-2 and TSX.Further, full PolSAR information is better suited for class separation than less diverse polarimetric feature spaces, like all dual-polarimetric measurements (HH/VV, HH/HV or VV/VH).The Kennaugh Matrices offered the highest class separability among the investigated decompositions, and among the ten elements of the quad-polarized C-and L-Band Kennaugh Matrix the elements K0, K1, K2, K3 and K4 were found to be most valuable for class discrimination.This also indicates that the phase-relation between HH and VV (K3, K4) provides crucial information for separating the investigated tundra land cover classes, since it contains the distinction of surface from diplane scattering.Further, the intensity-based information of the elements K0, K1 and K2, which explain the diagonal elements of the Coherency Matrix, are favorable for class discrimination.
In light of the results presented in this manuscript, future work should focus on investigating the combined use of short-and long-wave PolSAR data, e.g., of C-/X-Band and L-Band.It is anticipated that such multi-frequency data will provide complementary information useful for accurate classification and the description of land surface parameters, as well as biophysical parameters of the tundra vegetation.In this context, the combination of PolSAR information via a multi-sensor approach seems very promising, since it will combine dielectric, and geometrical properties of the targets.
An interesting future question will be to also address the use of hybrid-polarimetric/ compact-polarimetric data that can be synthetically generated from quad-polarimetric data, also using the Kennaugh Matrix approach [31].The question will be how such data perform compared to quad-, or dual-polarized data.
As well, the potential for land cover classification should be addressed, e.g., via the Random-Forest approach that was shown to provide an interesting classification framework also for PolSAR data [16,51].In this context, upcoming studies should further acknowledge if the Random-Forest approach is appropriate and essential for a successful PolSAR classification.
The inclusion of Sentinel-1 C-Band PolSAR data is another option, as the Interferometric Wide-Swath mode provides large spatial coverage at high spatial resolutions, and the planned continuity of the Sentinel SAR systems will offer the capacity to support long term monitoring and consistent remote observations of Arctic land covers.However, as shown in this study the VV/VH polarization mode, employed by Sentinel-1 over most parts of the Canadian Arctic, seems less suited for characterizing of the tundra land cover classes; thus the use of a multi-frequency or multi-sensor approach is advisable.
In summary, the SAR data of all wavelengths-also of the L-Band-were shown to provide important information about the tundra environment and utilization of such remotely sensed information is strongly recommended.PolSAR data provide unique information on dielectric, and geometrical properties that can help to increase the information space.Whenever possible dual-or quad-polarized data should be used, as polarimetry was shown to be of high value and importance.

Figure 1 .
Figure 1.Location of the test site: Tuktoyaktuk (Northwest Territories, Canada) and RGB composites of remote sensing imagery: (a) elevation and slope of intermediate TanDEM-X DEM, coverage of TerraSAR-X, Radasat-2, ALOS and ALOS-2 imagery, locations of in situ field work in 2010 and 2012, locations of land cover reference samples derived from high resolution ortho-photos, extent indicator of the subsequent sub-figures (red rectangle); (b) RGB false-color composite of sigma nought HV intensities [dB] of X-Band (TerraSAR-X, 2011), C-Band (Radarsat-2, 2011) and L-Band (ALOS, 2010); (c) RGB false-color composite of Kennaugh Matrix Elements K2, K1 and K3 of quad-polarized Radarsat-2 data (see Section 2.2.2 Polarimetric SAR Data and Decompositions for more details on the Kennaugh Matrix Elements); (d) Landsat TM (2011) true-color composite of red, green and blue surface reflectance.RGB channels are stretched linearly between 1 and 99% of the data range.

Figure 1 .
Figure 1.Location of the test site: Tuktoyaktuk (Northwest Territories, Canada) and RGB composites of remote sensing imagery: (a) elevation and slope of intermediate TanDEM-X DEM, coverage of TerraSAR-X, Radasat-2, ALOS and ALOS-2 imagery, locations of in situ field work in 2010 and 2012, locations of land cover reference samples derived from high resolution ortho-photos, extent indicator of the subsequent sub-figures (red rectangle); (b) RGB false-color composite of sigma nought HV intensities [dB] of X-Band (TerraSAR-X, 2011), C-Band (Radarsat-2, 2011) and L-Band (ALOS, 2010); (c) RGB false-color composite of Kennaugh Matrix Elements K2, K1 and K3 of quad-polarized Radarsat-2 data (see Section 2.2.2 Polarimetric SAR Data and Decompositions for more details on the Kennaugh Matrix Elements); (d) Landsat TM (2011) true-color composite of red, green and blue surface reflectance.RGB channels are stretched linearly between 1 and 99% of the data range.

Table 2 .
Land cover classes considered in the analysis; description, class abbreviations (Abbr.), and class color coding.Bare Ground samples were considered homogeneous if 10% or less were "other" materials or vegetation.The letters "W", "B", "T" of the class abbreviations refer to Wetland, Bare Ground, and Tundra land cover classes.

Table 4 .
Average class separability measured as Transformed Divergence (TD) for classes OL (Permanent Water Bodies), BS (Sand), BM (Mixed Sediment), BW (Driftwood Accumulation), HT (Herb Dominated Tundra), ST (Shrub Dominated Tundra), WI (Inundated Low Lying Tundra) and WT (Wetland) (see Table2).TD is ranging from 0 to 2000; higher values indicate better class separation.The feature spaces are ranked in descending order of the average separability (AV).The black bars are scaled linearly between the minimum and maximum of AV.The colors from red to yellow to green correspond to the 10%, 50% and 90% quantiles of the AV.

Table 2 )
. The features are ranked in descending order of the average increase in separability (AV).The last column displays the AV in percent (%) The colors from red to yellow to green correspond to the 10%, 50% and 90% quantiles of the AV data range.The metric displays the average increase in separability when a feature of interest (K0-K9) is added to an existing feature space (see Section 2.3 Correlation, Class Separability and Feature Selection).