Mapping Plant Functional Types in Floodplain Wetlands: An Analysis of C-Band Polarimetric SAR Data from RADARSAT-2

The inclusion of functional approaches on wetland characterizations and on biodiversity assessments improves our understanding of ecosystem functioning. In the Lower Paraná River floodplain, we assessed the ability of C-band polarimetric SAR data of contrasting incidence angles to discriminate wetland areas dominated by different plant functional types (PFTs). Unsupervised H/α and H/A/α Wishart classifications were implemented on two RADARSAT-2 images differing in their incidence angles (FQ24 and FQ08). Obtained classes were assigned to the information classes (open water, bare soil and PFTs) by a priori labeling criteria that involved the expected interaction mechanisms between SAR signal and PFTs as well as the relative values of H and α. The product obtained with the shallow incidence angle scene had a higher accuracy than the one obtained with the steep incidence angle product (61.5% vs. 46.2%). We show how a systematic analysis of the H/A/α space can be used to improve the knowledge about the radar polarimetric response of herbaceous vegetation. The map obtained provides novel ecologically relevant information about plant strategies dominating the floodplain. Since the obtained classes can be interpreted in terms of their functional features, the approach is a valuable tool for predicting vegetation response to floods, anthropic impacts and climate change.


Introduction
The largest wetlands of South America are associated with the floodplains of the Amazonas, Orinoco and Paraná-Paraguay Rivers [1]. These ecosystems feature high biodiversity and are important for society due to flood regulation, habitat provision for fish and wildlife, natural forage provision for cattle, and carbon sequestering, among other factors [2,3]. Floodplain wetlands have been mapped and monitored with optic data in several opportunities (e.g., [4][5][6]). However, optical data often fail in  Table 1) of 30 January 2011; and (d) Pauli representation of the steep incidence angle scene (beam FQ08, see Table 1) of 2 February 2011. In both (c) and (d), colors were assigned as follows: Red = |HH − VV|; Green = |HV + VH| = 2|HV|; Blue = |HH -VV|. The background image shows the near infrared band from a Landsat 5TM scene of 1 February 2011.

SAR Data Acquisition and Processing
Two polarimetric RADARSAT-2 scenes (C-band, frequency 5.405 GHz, wavelength 5.6 cm) were acquired on 30 January and 2 February, 2011 ( Table 1). The overlapping area of the scenes acquired covers 462 km 2 of freshwater wetlands, resulting in the region of interest for further analyses ( Figure 1). Incidence angles were contrasting: the first scene had a relatively shallow incidence angle (Beam FQ24) and the second one a steep incidence angle (FQ8) ( Table 1). Between acquisition dates,  Table 1) of 30 January 2011; and (d) Pauli representation of the steep incidence angle scene (beam FQ08, see Table 1) of 2 February 2011. In both (c) and (d), colors were assigned as follows: Red = |HH´VV|; Green = |HV + VH| = 2|HV|; Blue = |HH -VV|. The background image shows the near infrared band from a Landsat 5TM scene of 1 February 2011.

SAR Data Acquisition and Processing
Two polarimetric RADARSAT-2 scenes (C-band, frequency 5.405 GHz, wavelength 5.6 cm) were acquired on 30 January and 2 February, 2011 ( Table 1). The overlapping area of the scenes acquired covers 462 km 2 of freshwater wetlands, resulting in the region of interest for further analyses (Figure 1). Incidence angles were contrasting: the first scene had a relatively shallow incidence angle (Beam FQ24) and the second one a steep incidence angle (FQ8) ( Table 1). Between acquisition dates, the accumulated precipitation was negligible (2 mm, data from Servicio Meterológico Nacional (SMN), Argentina) and hydrometric level in the Paraná River increased only 13 cm (data from Prefectura Naval Argentina), remaining far below the level for which water begins to enter the floodplain [47]. Wind speed and wind direction (which can affect water roughness and the orientation of plant leaves) were similar in the two dates: gentle breezes between 12 and 15 km/h with east or northeast direction (data from SMN). Thus, backscattering differences between scenes were assigned only to the differing incidence angle.  Figure 2 shows a flow diagram for the pre-processing and processing of the two RADARSAT-2 scenes. PolSARpro [48] was chosen for most of the steps, since it can handle coherence and covariance matrices, filters, decompositions, segmentations, classifications and statistical analyses. Each SLC scene was imported and digital numbers were converted to sigma-nought values. Next, the T3 coherence matrix and the C3 covariance matrix were geocodified with MapReady [49]. Since the study area is relatively flat, no terrain correction with a digital elevation model was applied. The images were filtered for speckle reduction using the adaptive Lee filter [50] with a 7ˆ7 pixel window size. Geocoding before filtering diminished image artifacts. The intensities of the backscattering coefficients in HH, HV and VV were extracted and transformed to decibels (dB). Note that HV = VH if reciprocal targets are assumed [51].

Information Classes: Sampling and Characterization
Vegetation sampling was conducted during the summers of 2010-2011 and 2011-2012 in 39 georeferenced sampling sites distributed along the region of interest (Figure 1c,d). Dominant plant species were classified into plant functional types (PFTs) by means of a trait-based approach. The PFT classification was derived from a broader study, in which 16 functional traits were screened in 44 species that dominated 72 sites (39 sites in the study area of this work and 33 sites in surrounding areas of the Paraná River floodplain) (Morandeira and Kandus, in prep.). For some species, more than one population was sampled. Next, these populations were hierarchically classified into PFTs. Thus, we consider that the identified PFTs were representative of the Lower Paraná River floodplain. PFTs were characterized and assigned to a life strategy (competitor, ruderal, or stress-tolerant) [52,53] by interpreting their mean trait values and trait ranges. Five PFTs covered homogeneous areas of at least 1000 m 2 in the study area of this work, and are described in Table 2. The location of three sites with bare soil (cattle corrals or areas with high cattle transit) and of ten sites of open water (shallow lakes and rivers) was also recorded. This led to a total of 52 field observations, which were used for accuracy assessments in next steps. Site selection was limited by accessibility, but was independent of classification results (sampling was conducted before performing the classifications). Homogenous areas around these georeferenced sites were identified by visual interpretation of Pauli representations [51] derived from each scene. In order to assess whether the information classes (PFTs, bare soil and water) can be differenced by their backscattering coefficients in HH, HV and VV, permutation ANOVAs (999 permutations) followed by a posteriori Tukey's tests [54] were carried out. Permutation ANOVAs were run with the open-source software R [55] and the package lmPerm [54].

Information Classes: Sampling and Characterization
Vegetation sampling was conducted during the summers of 2010-2011 and 2011-2012 in 39 georeferenced sampling sites distributed along the region of interest (Figure 1c,d). Dominant plant species were classified into plant functional types (PFTs) by means of a trait-based approach. The PFT classification was derived from a broader study, in which 16 functional traits were screened in 44 species that dominated 72 sites (39 sites in the study area of this work and 33 sites in surrounding areas of the Paraná River floodplain) (Morandeira and Kandus, in prep.). For some species, more than one population was sampled. Next, these populations were hierarchically classified into PFTs. Thus, we consider that the identified PFTs were representative of the Lower Paraná River floodplain. PFTs were characterized and assigned to a life strategy (competitor, ruderal, or stress-tolerant) [52,53] by interpreting their mean trait values and trait ranges. Five PFTs covered homogeneous areas of at least 1000 m 2 in the study area of this work, and are described in Table 2. The location of three sites with bare soil (cattle corrals or areas with high cattle transit) and of ten sites of open water (shallow lakes and rivers) was also recorded. This led to a total of 52 field observations, which were used for accuracy assessments in next steps. Site selection was limited by accessibility, but was independent of classification results (sampling was conducted before performing the classifications). Homogenous areas around these georeferenced sites were identified by visual interpretation of Pauli representations [51] derived from each scene. In order to assess whether the information classes (PFTs, bare soil and water) can be differenced by their backscattering coefficients in HH, HV and VV, permutation ANOVAs (999 permutations) followed by a posteriori Tukey's tests [54] were carried out. Permutation ANOVAs were run with the open-source software R [55] and the package lmPerm [54].

H/A/α Segmentation
The Cloude-Pottier polarimetric decomposition expresses the 3ˆ3 complex coherence matrix as a sum of eigenvalues (λ i ) and eigenvectors (e i ) [56]. Three components can be obtained from the T3 matrix: entropy (H), angle α and anisotropy (A) [57]. The entropy ranges between 0 and 1: values near 0 indicate a single scattering mechanism, whereas values near 1 indicate a random combination of different scatterers [56]. The average α component ranges between 0˝and 90˝and depends on the dominant scattering mechanism: surface scattering for values near 0˝, dipole scattering for values near 45˝and double-bounce scattering for values near 90˝ [58]. Lastly, low A values are obtained for a main scattering mechanism (case H = 0 or H = 1) and two non-negligible secondary mechanisms, whereas high A values correspond to two main scattering mechanisms with the same probability of occurrence and a third mechanism of minor importance [59].
After extracting the three Cloude-Pottier components, the pixel density distribution of each scene was plotted in an H/α plane. An H/α segmentation was obtained by dividing the plane into nine zones (a combination of low/medium/high H and low/medium/high α) [56]. Thresholds between the zones are somewhat arbitrary. A curve delimits mathematically infeasible solutions (high H and low α), so that eight zones remain feasible. By introducing A, it is possible to obtain two H/α planes, one for A > 0.5 and the other for A < 0.5. Thus, an H/A/α segmentation leads to 16 zones with distinct scattering behavior.

Wishart Unsupervised Classifications on the Coherence Matrix
The Wishart H/α or H/A/α unsupervised classifiers are non-hierarchical statistical algorithms that require a Cloude-Pottier segmentation product [58]. In each scene, iteration optimization procedures were initialized from the centroids of the H/α or of the H/A/α segmentations (8 or 16 feasible classes, respectively). The classifications were computed on the T3 matrix by means of a complex algorithm based on the Maximum Likelihood classifier, with a convergence level of 99%, a maximum of 100 iterations and a 3ˆ3 pixel window size. Thus, two Wishart unsupervised classifications were obtained for each scene, one derived from the H/α plane and the other from the H/A/α plane.
Classes derived from the unsupervised H/α classification were assigned to the information classes (water, bare soil or a PFT) according to their H and α values and following a priori predictions. We considered that H increased with randomness and geometric anisotropy of the scatterers, whereas α varied with the dominant interaction mechanism between the signal and the target (low α for surface scattering, medium α for volume scattering and high α for double-bounce scattering). The water information class was expected to have low values of both H and α, and to be identified due to the known location of the main rivers and shallow lakes. Surface incoherent scattering was expected in bare soil targets, with low α and intermediate H values. Regarding PFTs, those presenting double-bounce scattering (mainly due to vertical stems) were expected to show the highest α (PFT A). Previous works have addressed the contribution of double-bounce scattering in Schoenoplectus californicus [60], which is a species included in PFT A and with a similar structure to Cyperus giganteus. Besides, the overall biomass and the degree of randomness in the canopy (distribution of leaves and stems) were assumed to be related to H (see Table 2 for predictions of the scattering mechanisms of PFTs). PFTs with broadleaf leaves and stems were expected to have higher randomness in biomass distribution (particularly PFT C, but also PFT B), and biomass was expected to be related to volume scattering. In this sense, we assigned Wishart classes to PFTs assuming an increasing order of H: PFT D < PFT E < PFT B < PFT C, from short and low biomass grasslands (PFT D) to tall and high biomass broadleaf marshes (PFT C). Dihedral scatterers such as ships or buildings were associated to very high α and very low H (these targets occurred in the segmented scenes but were outside the overlapping study area). This labeling procedure was a working hypothesis: the expected H and α values for different information classes could be inferred from the structural features of the targets. Accuracy assessment with field data (see Section 2.3.4) was used to discuss whether this hypothesis was accurate for our purpose. A second class labeling procedure, not associated with hypotheses on the vegetation-signal interaction, was conducted by maximizing the Kappa index.
The discrimination of targets with high entropy might be improved by including the anisotropy component (as suggested by [51]). Thus, a progressive classification procedure was carried out for those classes derived from the H/α Wishart classification with very high entropy values (H > 0.90). The H/α Wishart classes were further divided taking into account the classes from the unsupervised H/A/α Wishart classification. Next, whether PFT discrimination was improved using H/A/α classes to re-allocate pixels to the information classes was analyzed.

Accuracy Assessment
The accuracy of the Wishart unsupervised classifications was assessed in the 52 georeferenced sampling sites included in the region of interest (39 sites of PFTs, three of bare soil, ten of open water) (see Section 2.3.1 for details on site sampling). Overall accuracy was computed as the number of pixels correctly classified out of the total number of pixels. To evaluate how much the resulting accuracy differed from a random assignment of class labels, the Kappa index was computed [61]. The accuracy observed was also expressed as the proportion of the accuracy that may be obtained with an easier and simpler classification ("naïve" classification, sensu [62]). In this study, the accuracy of a straightforward Isodata unsupervised clustering of SAR data expressed in the canonical polarization base was considered for comparison purposes [63]. The Isodata classification was carried out on the pixels of the multipolarization stack (HH, HV and VV) (seven classes, 100 iterations or a convergence level of 0.99). A median filter with a 3ˆ3 window was applied on the product of the Isodata classification so that it could be compared to the Wishart product. Both Wishart and Isodata are non-hierarchical classifiers and compute the distances of the pixels to their respective clusters centroids in each of the iterations. Thus, the accuracy of the Wishart classification compared to the accuracy of this multipolarization stack might be understood as a measure of the improvement of the classification due to the use of the complex polarimetric information contained in the T3 matrix. This is particularly important, since polarimetric information extraction techniques are more complex and time-costly and therefore their use should be justified.

Description of the Scenes and Field Samples
In both RADARSAT-2 scenes, co-polarizations (HH and VV) showed similar backscattering coefficients for most of the field samples belonging to different classes (Figure 3), except for PFT E (tall grasslands, see Table 2) for which backscattering was higher and more variable in HH than in VV. The analysis of field samples representative of each class suggests that using only one polarization is quite inappropriate for classifying the region of interest, especially with regard to bare soil discrimination ( Figure 3). However, HV (or VH) seems to be the preferred polarization for discriminating wetland classes with RADARSAT-2 scenes. Differences between information classes were higher in the cross-polarization band (HV) than in co-polarization bands (HH or VV). This issue has been pointed out in [11]: the depolarizing character of the vegetation differs according to vegetation structure, biomass distribution and flooding state.
Characterized mainly by mirror reflection, water had low backscattering coefficients. However, water scattering was higher in the steep angle scene than in the shallow one, as expected for a slightly rough surface. The highest backscattering was observed for PFT A samples (bulrush marshes), probably due to the contribution of double-bounce scattering that was reported and modeled in [60]. The lowest scattering among vegetation classes occurred for PFT D, grouping short plants with low cover and biomass, which probably allowed large signal penetration up to rough soil or water. Bare soil backscattering was always similar to the backscatter of at least one PFT, thus suggesting difficulty to discriminate bare soil in further classifications (especially if polarimetric information is not included). Although the backscattering of some classes showed differences between the scenes (Figure 3), these differences were not evident in the Pauli representations (Figure 1c,d) and high cross-polarization response (high |HV|) predominated in both scenes.

H/α Segmentations
The interaction between PFT targets and the polarimetric SAR signal was heavily influenced by the incidence angle (Figure 4), assuming no substantial changes in vegetation, hydrometric or wind conditions during the two days between scene acquisitions (see Section 2.2). In line with the predominance of volume scattering at C-Band from vigorous herbaceous vegetation, most of the pixels of both scenes had high H (H > 0.8) and intermediate α values (30˝< α < 50˝). Thus, segmenting the H/α plane in eight zones was not accurate for discriminating PFTs and bare soil: the non-water areas were grouped in the same class in both segmented images (Figure 4a,b). The density of pixels in zones with low α and low or intermediate H, thus presenting simple interaction mechanisms, was higher in the steep angle scene than in the shallow angle scene (Figure 4c vs. Figure 4d). In the segmented image from the steep angle scene (Figure 4b), those pixels were located in wetland areas with high H values in the shallow angle scene. These differential signal-target interactions according to the incidence angle may indicate the presence of flooded vegetation. Indeed, with a shallow incidence angle, volume scattering with one or more vegetation layers occurred, so that the signal was scattered before reaching  Table 2 for a description of their features): (a) shallow incidence angle scene, HH polarization; (b) steep incidence angle scene, HH polarization; (c) shallow incidence angle scene, HV polarization; (d) steep incidence angle scene, HV polarization; (e) shallow incidence angle scene, VV polarization; and (f) steep incidence angle scene, VV polarization. In each figure, the water class is plotted in a different scale than that of the bare soil and Plant Functional Types (PFT) classes. In the boxplots, boxes show the medians (thick line) and the 1st and 3rd quartiles (limits of the boxes); whiskers show the maximum and minimum values; outliers were omitted for clarity. The letters onto the boxplots represent significant differences in the backscattering coefficients (ANOVA tests, followed by Tukey comparisons) at 0.01 significance level.
Characterized mainly by mirror reflection, water had low backscattering coefficients. However, water scattering was higher in the steep angle scene than in the shallow one, as expected for a slightly rough surface. The highest backscattering was observed for PFT A samples (bulrush marshes), probably due to the contribution of double-bounce scattering that was reported and modeled in [60]. The lowest scattering among vegetation classes occurred for PFT D, grouping short plants with low cover and biomass, which probably allowed large signal penetration up to rough soil or water. Bare soil backscattering was always similar to the backscatter of at least one PFT, thus suggesting difficulty to discriminate bare soil in further classifications (especially if polarimetric  Table 2 for a description of their features): (a) shallow incidence angle scene, HH polarization; (b) steep incidence angle scene, HH polarization; (c) shallow incidence angle scene, HV polarization; (d) steep incidence angle scene, HV polarization; (e) shallow incidence angle scene, VV polarization; and (f) steep incidence angle scene, VV polarization. In each figure, the water class is plotted in a different scale than that of the bare soil and Plant Functional Types (PFT) classes. In the boxplots, boxes show the medians (thick line) and the 1st and 3rd quartiles (limits of the boxes); whiskers show the maximum and minimum values; outliers were omitted for clarity. The letters onto the boxplots represent significant differences in the backscattering coefficients (ANOVA tests, followed by Tukey comparisons) at 0.01 significance level.
Few pixels with high α (indicating double-bounce scattering) occurred in both scenes. The ones with low H were located in the middle of the Paraná River or outside the region of interest in urban areas, pointing out dihedral reflectors (e.g., ships, buildings, metal roofs). Besides, the pixels with high α and intermediate or high H were probably related to double-bounce from vegetation targets. Their locations were coincident in both segmented images, but the covered area was slightly higher in the shallow incidence angle scene than in the steep one.
in the shallow angle scene. These differential signal-target interactions according to the incidence angle may indicate the presence of flooded vegetation. Indeed, with a shallow incidence angle, volume scattering with one or more vegetation layers occurred, so that the signal was scattered before reaching water below vegetation. On the contrary, the steep signal was less attenuated by the vegetation and mainly interacted with water below it. Few pixels with high α (indicating double-bounce scattering) occurred in both scenes. The ones with low H were located in the middle of the Paraná River or outside the region of interest in urban areas, pointing out dihedral reflectors (e.g., ships, buildings, metal roofs). Besides, the pixels with high α and intermediate or high H were probably related to double-bounce from vegetation targets. Their locations were coincident in both segmented images, but the covered area was slightly higher in the shallow incidence angle scene than in the steep one.
These analyses of the H/α space improved the previous knowledge about the response of herbaceous wetland vegetation to SAR signal. Most of the previous studies have described high H (d) steep incidence angle scene, density of pixels in the H/α plane; (e) shallow incidence angle scene, detail of the centroids assigned to bare soil and plant functional types (PFTs) classes (H > 0.80 and α > 0.36); (f) steep incidence angle scene, detail of the centroids assigned to bare soil and PFTs classes (H > 0.80 and α > 0.36). In both (a) and (b), the background image shows the near infrared band from a Landsat 5TM scene of 1 February 2011. In (c) and (d), the H/α planes were divided in nine zones and the centroids of the classes obtained with the Wishart classifier are shown; DR stands for "dihedral reflector" and W stands for "water". In (e) and (f), BS stands for "bare soil" and the letters A-D indicate the corresponding PFT. Class labeling followed a priori criteria (see Section 2.3.3) and was a working hypothesis assessed with field data (see Section 3.4).
These analyses of the H/α space improved the previous knowledge about the response of herbaceous wetland vegetation to SAR signal. Most of the previous studies have described high H areas as forests and high α areas as cities (e.g., [23,64,65]). However, high biomass vegetation assessed in this study featured very high H. According to the a priori labeling criteria, the location of the PFT A pixels in the H/α plane was similar to the location of pixels with floodable trees and shrubs in the Amazonas floodplain [32]. High α values in pixels dominated by the equisetoid herbs of PFT A (e.g., Schoenoplectus californicus in bulrush marshes) presented high α values such as predicted due to its high double-bounce scattering in presence of a water film [60].

Unsupervised Wishart H/α and H/A/α Classifications
The H/α Wishart classifiers converged to eight classes by shifting most of the centroids of the segmented plane to areas with high H and intermediate α (Figure 4c,d). In both classifications, the centroid of one of the classes had very high α (α > 60˝, dihedral reflectors, big ships). Thus, this class was excluded in further analyses. Next, the class with the lowest H and α was assigned to the open-water class. Among the remaining classes, the one with the lowest H (and low α) was assumed to belong to bare soil. The remaining classes were labeled as stated in Section 2.3.3 (Figure 4e,f). Both the bare soil and the PFT D (short grasslands, see Table 2) classes had lower α in the steep angle scene than in the shallow angle one.
In the shallow angle classification product, four classes had H > 0.90 and were provisionally assigned to PFT E, A, B and C (in increasing order of H) (Figure 4e). For these four classes with H > 0.90, a progressive classification procedure was followed by analyzing H/A/α subclasses, i.e., those classes obtained with the H/A/α classifier within each H/α class (see scheme in Figure 2). Within PFT B, five H/A/α subclasses occurred. The geographic location of one out of these H/A/α subclasses was coincident with field records of sites dominated by PFT A, whereas other H/A/α subclass was coincident with field records of sites dominated by PFT C. The remaining H/A/α classes were retained in PFT B class. Thus, PFT discrimination was slightly improved by reassigning these two subclasses. The number of reassigned pixels was a minority compared to the total pixels in PFT B class. None of the pixels of PFT A, E and C were reassigned according to the H/A/α subclasses. The second labeling procedure (i.e., maximizing the Kappa index when assigning classes obtained to information classes) was coincident with the labeling carried out by means of our a priori criteria.
In the steep angle classification product, only one class had a mean H higher than 0.90 (PFT C) (Figure 4f), so that the progressive classification on H/A/α planes was omitted. The two labeling procedures led to different results. Class assignment to information classes following the a priori criteria is shown in Figure 4d,f. Besides, Kappa index was maximized for this scene when pixels first assigned to PFT B were re-assigned to PFT E, and vice versa. This result suggests that our physical hypotheses on vegetation were contrary to the observed pattern for these two PFTs, at least for the steep incidence angle scene. Our a priori criteria could be misleading due to the difficulty to infer mean H and α values from biomass and non-quantitative vegetation structure data (distribution of biomass, degree of randomness in the canopy). In addition, the fact that our working hypothesis seems to be more accurate for the shallow incidence angle scene than for the steep one, is consistent with the general good performance of shallow incidence angle scenes for discriminating herbaceous wetland vegetation in the literature [11].

Comparison between Incidence Angles and Accuracy Assessment
The spatial pattern of PFT classes differed according to incidence angles (Figure 5a vs. Figure 5b). Both maps were quite homogeneous without much noise (i.e., zones of vegetation were identifiable). The zonation of vegetation is typical in wetland ecosystems, where abiotic constraints are important for plant community development, and has been described for the study area [18]. The classification products obtained with the shallow and the steep incidence angle scenes were coincident in 30.0% (a concordance of 17.8% in comparison with the randomly expected concordance, as estimated by the Kappa index). Water was almost equally classified in both scenes. Most of the differences between the resulting products arose for PFT E (tall grasslands), for which we identified a total area of 31 km 2 in the shallow angle product and an area of 149 km 2 in the steep angle one. Besides, 70.4% of the pixels assigned to PFT D in the shallow angle product were assigned to PFT B in the steep angle one, confounding short grasslands and short broadleaf marshes. In addition, much of the areas identified as PFT C in the shallow angle product were identified as PFT E in the steep angle one, and vice versa; confounding tall broadleaf marshes and tall grasslands. This mismatch between the maps suggests that broadleaf and graminoid herbs interact differently with the SAR signal depending on the incidence angle.
confounding short grasslands and short broadleaf marshes. In addition, much of the areas identified as PFT C in the shallow angle product were identified as PFT E in the steep angle one, and vice versa; confounding tall broadleaf marshes and tall grasslands. This mismatch between the maps suggests that broadleaf and graminoid herbs interact differently with the SAR signal depending on the incidence angle.  Table 2) of 30 January 2011; (b) Steep incidence angle scene (beam FQ08, see Table 2) of 2 February 2011. PFT stands for "plant functional type". For details on the classes, see Table 1. Maps correspond to the a priori labeling criteria.
Accuracy was higher for the shallow incidence angle scene than for the steep one. On the one hand, the product obtained with the shallow incidence angle scene had an overall accuracy of 61.5% or of 52.4% excluding the water class ( Table 3). The class with lowest commission error was PFT B, whereas the one with lowest omission error was PFT C (Table S1). PFT D and PFT E (short and tall grasslands, respectively) were confounded. Sixty percent of the pixels belonging to PFT E were omitted and assigned to other PFTs and to bare soil. The accuracy differed from the chance-adjusted expected accuracy and the Kappa index was 54.8% or 42.5% (including or not water class, respectively) ( Table 3). In comparison with the Isodata classification on HH, HV and VV, the Wishart classifier improved the accuracy by 20.7% (including the water class) or by 34.2% (without the water class), in line with previous comparisons of the performance of Wishart classifications derived from quad-pol imagery and of classifiers based on standard linear dual-pol imagery [66]. Regarding the omission errors (Table S1), some of the pixels classified as open water in marsh areas might actually correspond to flooded vegetation. In areas of the Amazonas floodplain dominated by macrophytes, mirror reflection has been reported for emergent biomass (above water level) up to 200 g/m 2 , from which volume scattering starts increasing [13]. Thus, some the flooded areas of PFT D and E could be underestimated and assigned to the water class. Table 3. Accuracy assessment of the classification products obtained with the RADARSAT-2 scenes. For the shallow incidence angle scene, the two class labeling procedures led to the same result, whereas for the steep incidence angle scene the results differed if the a priori labeling criteria were used or the Kappa index was maximized.   Table 2) of 30 January 2011; (b) Steep incidence angle scene (beam FQ08, see Table 2) of 2 February 2011. PFT stands for "plant functional type". For details on the classes, see Table 1. Maps correspond to the a priori labeling criteria.

Class labeling
Accuracy was higher for the shallow incidence angle scene than for the steep one. On the one hand, the product obtained with the shallow incidence angle scene had an overall accuracy of 61.5% or of 52.4% excluding the water class ( Table 3). The class with lowest commission error was PFT B, whereas the one with lowest omission error was PFT C (Table S1). PFT D and PFT E (short and tall grasslands, respectively) were confounded. Sixty percent of the pixels belonging to PFT E were omitted and assigned to other PFTs and to bare soil. The accuracy differed from the chance-adjusted expected accuracy and the Kappa index was 54.8% or 42.5% (including or not water class, respectively) ( Table 3). In comparison with the Isodata classification on HH, HV and VV, the Wishart classifier improved the accuracy by 20.7% (including the water class) or by 34.2% (without the water class), in line with previous comparisons of the performance of Wishart classifications derived from quad-pol imagery and of classifiers based on standard linear dual-pol imagery [66]. Regarding the omission errors (Table S1), some of the pixels classified as open water in marsh areas might actually correspond to flooded vegetation. In areas of the Amazonas floodplain dominated by macrophytes, mirror reflection has been reported for emergent biomass (above water level) up to 200 g/m 2 , from which volume scattering starts increasing [13]. Thus, some the flooded areas of PFT D and E could be underestimated and assigned to the water class. Table 3. Accuracy assessment of the classification products obtained with the RADARSAT-2 scenes. For the shallow incidence angle scene, the two class labeling procedures led to the same result, whereas for the steep incidence angle scene the results differed if the a priori labeling criteria were used or the Kappa index was maximized. On the other hand, the product obtained with the steep incident angle scene had an overall accuracy of 46.2% or of 35.0% excluding the water class (Table 3). Kappa indexes were low and, when excluding the water class, the 95% confidence interval included the zero. High omission errors were observed for all the PFTs and were the highest for PFT E (Table S1). Optimum class labeling in terms of Kappa maximization improved the results and led to an overall accuracy of 53.3% and a Kappa index of 42.9%. Even with the Kappa maximization criterion, the accuracy was higher for the shallow incidence angle scene than for the steep one, being the former the more reliable product.
The results suggest that our hypotheses between PFT and SAR signal interactions were accurate only for the case of shallow incidence angles. In addition, our results were consistent with the reported preference of shallow incidence angle SAR signals for discriminating vegetation types and of steep incidence angle signals for detecting water below marsh vegetation and for flood monitoring [11,12]. In our study, at least two observations support this fact: in comparison to the shallow incidence angle scene, the steep one showed a higher density of pixels with low H and low α (Figure 4d vs. Figure 4c); and the centroids of classes PFT B and C featured lower H and α. The observed pattern can be explained by considering that steep SAR signal is less attenuated by vegetation layers and can interact to soil or water below the vegetation.
Regarding double-bounce scattering, the interaction mechanism seems to be more common for targets observed with high incidence angle scenes: pixels with medium H and high α in the shallow incidence angle scene showed low H and low α in the steep one. In the case of flooded herbaceous vegetation with low emergent biomass, a vertical incidence angle may diminish the probability of a first bounce in vegetation and a second one in the water (or vice versa), increasing the probability of reaching water. This observation is in marked contrast to what happens in forest environments, where for a given polarization and frequency, double-bounce scattering is more important in steep than in shallow incidence angles scenes [12,25,67].
In the Lower Paraná River floodplain, previous vegetation maps performed with optic data have overestimated the forest coverage at the expense of herbaceous wetlands [68]. Thus, the map obtained with the shallow incidence angle scene is a better start point than optical-derived maps for monitoring wetland vegetation in the floodplain. In the future, the study of SAR/optical fusion schemes will be considered. Our results suggest that most of the area was dominated by a ruderal strategy or an intermediate ruderal-competitor strategy (PFTs B, C and E; 58.7% of the total area or 63.8% excluding water). Ruderal plants are characterized by fast growth rate, short life span and high seed production, thus are favored by chronically disturbed but potentially productive environments [52]. The hydrological regime periodically disturbs the marsh areas, along with cattle activity. Ruderal plants can quickly regenerate between two ordinary flooding events, thus exploiting and nutrients brought by the river and taking advantage of resource opportunities. In addition, plant species with intermediate ruderalness-competitiveness may play a crucial role in vegetation adaptation to wetland states or in disturbance resistance. The fact that most of the floodplain seems to be covered by ruderal or intermediate ruderal-competitor plants gives to the environment a high ability to recover after disturbances and ordinary floods [69,70].

Conclusions
In this study, the ability of C-Band polarimetric data to monitor plant functional types in the Lower Paraná River floodplain was assessed. Although it is easier to perform a non-polarimetric analysis of the data, our results showed that analyzing the complex polarimetric information improves the resulting accuracy in comparison with the use of a multipolarization stack (HH, HV and VV). Moreover, the use of polarimetric decompositions facilitated the interpretation of the interaction mechanisms between the SAR signal and the scatterers. The analysis of the H/α and H/A/α spaces improved the previous knowledge about the response of herbaceous wetland vegetation in the Lower Paraná River floodplain, particularly with regard to the interpretation of classifications results. The analysis of the two images with different incidence angles allowed confirming an important result widespread in the literature: the backscattering measured at shallow incidence angles is more sensitive to herbaceous vegetation structure. Moreover, a conclusion of this study is that a priori hypotheses on the physical interaction between SAR signal and vegetation types should be dependent on the incidence angle. We conclude that using a C-Band shallow incidence angle scene is more accurate for the study of herbaceous PFTs in the Paraná River floodplain than using a steep incidence angle scene.
Taking into account the overall complexity of the landscape in the area, where herbaceous plants dominate, the shallow incidence angle product adds significant information while compared with previous available maps. The map obtained provides novel ecologically relevant information, suggesting that ruderal and ruderal-competitor broadleaf herbs dominate the floodplain. Mapping PFTs is not a common application in remote sensing yet, especially in wetlands, but it is a promising tool for assessing wetland functional features and for predicting vegetation responses to floods, anthropic impacts and climate changes in large floodplains. In the future, this approach can be used to assess changes in functional diversity or to evaluate how plant strategies respond to environmental changes.
Supplementary Materials: The following are available online at www.mdpi.com/2072-4292/8/3/174, Table S1: Error matrices of the Wishart classification products obtained with the RADARSAT-2 scenes. Overall accuracy and Kappa indexes are provided in Table 3.