Evaluation of Simplified Polarimetric Decomposition for Soil Moisture Retrieval over Vegetated Agricultural Fields

This paper investigates a simplified polarimetric decomposition for soil moisture retrieval over agricultural fields. In order to overcome the coherent superposition of the backscattering contributions from vegetation and underlying soils, a simplification of an existing polarimetric decomposition is proposed in this study. It aims to retrieve the soil moisture by using only the surface scattering component, once the volume scattering contribution is removed. Evaluation of the proposed simplified algorithm is performed using extensive ground measurements of soil and vegetation characteristics and the time series of UAVSAR (Uninhabited Aerial Vehicle Synthetic Aperture Radar) data collected in the framework of SMAP (Soil Moisture Active Passive) Validation Experiment 2012 (SMAPVEX12). The retrieval process is tested and analyzed in detail for a variety of crops during the phenological stages considered in this study. The results show that the performance of soil moisture retrieval depends on both the crop types and the crop phenological stage. Soybean and pasture fields present the higher inversion rate during the considered phenological stage, while over canola and wheat fields, the soil moisture can be retrieved only partially during the crop developing stage. RMSE of 0.06–0.12 m3/m3 and an inversion rate of 26%–38% are obtained for the soil moisture retrieval based on the simplified polarimetric decomposition.


Introduction
For agricultural fields, soil moisture is an important factor for several physical processes such as water conservation, soil erosion and surface runoff.In contrast to the conventional soil moisture acquisitions consisting of point samplings which are time and labor consuming, Synthetic Aperture Radar (SAR) presents the potentials to extract (sub-) surface information with high spatial and temporal resolution.Furthermore, polarimetric SAR extends the observation space, providing more promising approaches for the retrieval of the soil parameters under different crop conditions, types and phenological stages [1,2].Nevertheless, over agricultural vegetated fields, the effects of vegetation layers (mainly described by shape, structure, orientation and phenological developing stage) and the underlying soils (mainly described by soil moisture and surface roughness) are coherently superimposed in the measured SAR signature [1][2][3][4].As a consequence, the soil moisture retrieval over vegetated agricultural soils is a challenge addressed by several authors in the past decades [1][2][3][4][5][6][7][8][9].On one hand, the vegetation layer is modeled and removed using incoherent radiative transfer approach [6].Such a model considers the total backscattering coefficient as an incoherent sum of the scattering contributions from the vegetation volume scattering, the soil scattering and their interaction.For instance, the well-known water-cloud model [7], which is widely used for soil moisture retrieval, is a simple first-order radiative transfer solution which neglects the multiple scattering, and assumes that the vegetation and the underlying soil contributions are incoherently added up.On the other hand, based on the polarimetric SAR measurements, the target decomposition methodologies [10][11][12] are used to remove the vegetation component and estimate the underlying soil moisture from the ground scattering component.For example, the three-component model-based polarimetric decomposition [10] has been adapted [1][2][3] for soil moisture estimation under vegetation.In addition to the aforementioned soil moisture retrieval methods from radiative transfer models and polarimetric decomposition, empirical relationships [5,[13][14][15][16][17] were also developed between polarimetric SAR-derived parameters, vegetation and/or soil variables in order to improve the soil moisture retrieval in an easier way.
Considering the known sensitivity of SAR signature to soil and vegetation parameters, the SMAP (Soil Moisture Active Passive) Validation Experiment in 2012 (SMAPVEX12) was established for a launch assessment of SMAP products [18,19].Within the framework of SMAPVEX12, both airborne radar data and expansive ground measurements of soil and vegetation characteristics were collected in order to evaluate, enhance and develop soil moisture retrieval algorithms.
In this paper, we will take advantage of the polarimetric Uninhabited Aerial Vehicle Synthetic Aperture Radar (UAVSAR) data of SMAPVEX12 collected during most part of crop growth season to evaluate an existing polarimetric decomposition soil moisture retrieval algorithm [1].Indeed, both the surface and dihedral scattering components are used in [1] to retrieve soil moisture.Nevertheless, low accuracy and retrieval rate were obtained from the dihedral component which is a function of vegetation and soil characteristics.Furthermore, at a 45 ˝incidence angle, the dihedral component presents an ambiguity to separate the vegetation and soil moistures, hindering the robustness of soil moisture retrieval over vegetated fields.If the dihedral component is included in the retrieval process [1], the overall accuracy decreased.Thus, this paper proposes a simplification of [1] for soil moisture retrieval at L-band by using only the surface scattering component.This simplified algorithm for soil moisture retrieval is evaluated over a variety of crops, from the early development to the mature stage with maximum height (and biomass).The extensive field data and UAVSAR time series allow this study to focus on the detailed analysis of the model-based decomposition which is conceptually proposed and briefly analyzed in [1].For instance, in addition to the contribution of vegetation biomass and height, the effect of vegetation water content is analyzed on both the decomposed volume and surface scatterings.Thus, with the detailed data of SMAPVEX12, a thorough reasoning can be conducted in this study, and it can justify the inferences in [1] to confirm eventually the potential of model-based decomposition for soil moisture estimation under vegetation.In this paper, Section 2 describes the datasets acquired in the SMAPVEX12 campaign.Then, the retrieval method is proposed in Section 3. The main results are analyzed and discussed in Section 4 and the conclusion is presented in Section 5.

Study Site
The study site is that of the SMAPVEX12 experiment [18,19] shown in Figure 1.Since SMAPVEX12 is already documented, we only focus on information relevant for the understanding of this paper.For more details on the SMAPVEX12 site and dataset acquisition, the reader should refer to [18,19].The study site covers approximately 15 km ˆ70 km and is located within the larger Red River watershed in Manitoba, Canada [18].It is dominated by several annual crops, including cereals (around 32.2% of area), canola (13.2%), corn (7.0%) and soybean (6.7%).In addition, approximately 16.4% of the site is occupied by perennial cover (such as grassland and pasture), and some areas are also covered by wetlands and forests [18].Furthermore, the soil texture varies across the site, thus providing a large range of soil moisture conditions.These diversities in land cover (especially, different crop types) and soil types represent a good opportunity to develop and test the soil moisture retrieval algorithm over vegetation-covered soil.

Airborne SAR Data Acquisition and Processing
During SMAPVEX12, the L-band UAVSAR data were acquired in fine quad-polarization configuration, at 25 ˝to 65 ˝incidence angles range, and from 17 June to 17 July 2012.Since the same flight line (#31606) was used to conduct this study, a given field was observed with the same incidence angle from one acquisition date to another.Table 1 presents the temporal availability of UAVSAR measurements.In this study, the calibrated multi-looked cross products (MLC) with spatial resolution of 5.0 m in range by 7.2 m in azimuth are used.The MLC data were firstly extracted as a coherency matrix T3 using PolSARpro4.2 [20].Then, a boxcar filter with 7-by-7 windows was applied to reduce the speckle effect.As the study area has no topographic variation [19], no compensation of terrain slope was applied on the SAR images.

Ground Measurements
They consisted of both soil and vegetation characteristics measurements.They were conducted during six weeks from 6 June to 17 July 2012 in order to provide suitable ground measurements for interpreting the collected remote sensing data (including the UAVSAR measurements) with the changes in both soil and vegetation characteristics, developing and testing soil moisture retrieval algorithms under several vegetation types and crop growing stages including the peak season of vegetation cover [18].
(1) Soil Moisture: Almost coincident with each UAVSAR overpass, the volumetric soil moistures of 55 fields were measured at about 6 cm depth (probe length) using hand-held hydra probes (except for the acquisition on 19 June 2012).Furthermore, the meteorological observations near the test site were simultaneously recorded, providing the experimental conditions for UAVSAR data acquisition and related ground measurements.In Figure 2, the temporal variation of the measured volumetric soil moisture is in agreement with that of the rainfall amount.The peak of the volumetric soil moisture appeared around Julian day 168 as a consequence of successive rain events followed by a significant precipitation of about 30 mm on Julian 168.(2) Surface Roughness: Soil surface roughness parameters described by the RMS height and the auto correlation length were measured over the 55 fields using a 1 m long profilometer [18].Table 2 presents the range of these measurements for canola, corn, pasture, soybean and wheat.The observed field to field variability of roughness parameters results mainly from agricultural practices [19].The roughness measurements are subsequently used (Section 4.1) to compensate the attenuation of the microwave propagation in order to analyze the scattering mechanisms.(3) Vegetation parameters: Measurement of vegetation properties was also an important aspect of the SMAPVEX12 campaign with a main focus on representative crops (canola, corn, soybeans and wheat).The vegetation height and biomass were measured for analyzing the vegetation effects on the L-band microwave signal and subsequently on soil moisture retrieval.Besides, in a post campaign activity, maps of the vegetation water content over the SMAPVEX12 site were derived from SPOT-4 and RapidEye data with an accuracy ranging from 0.13 kg/m 2 to 0.80 kg/m 2 depending on the crop types [21].For the illustration, Figure 3 shows the measured vegetation height and net wet biomass [18] and the derived vegetation water content [21] of representative crops for the phenological stage considered in this study.

Method
Polarimteric decomposition is usually used to deconvolve the different scattering components over a vegetation-covered soil.For soil moisture retrieval using SMAPVEX12 datasets, the present study evaluates the applicability of a simplified version of an existing model-based polarimetric decomposition technique [1] which was developed from Freeman-Durden/Yamaguchi target decomposition methods [10,11].

Theoretical Aspects of Polarimetric Decomposition
Assuming the reflection symmetry for natural media, the coherency matrix [T3] is modeled as a sum of three scattering mechanisms [1,2], namely the surface, the dihedral and the volume scatterings.
rT3s " » --- (1) Surface scattering: In Equation (1), the surface scattering mechanism is modeled by the X-Bragg model.The surface roughness term δ (the empirical value of π/6 in [1] is also optimal in our study) is introduced into the X-Bragg model by rotating the coherency matrix around line of sight (LOS).Consequently, in contrast to the classical Bragg model, the X-Bragg model is applicable for a wider range of surface roughness conditions (ks < 1 with k and s, the wavenumber and the surface rms height, respectively).In addition, two other parameters are used to describe the surface scattering component: β is the normalized difference of the Bragg scattering in H and V polarizations and f s is the single scattering amplitude [1]: where m s accounts for the microwave attenuation due to the vegetation layer and surface roughness.The Bragg scattering coefficients (R h and R v ) depend on soil dielectric constant ε s and incidence angle θ as follows: R h " cosθ ´aε s ´sin 2 θ cosθ `aε s ´sin 2 θ and R v " pε s ´1qpsin 2 θ ´εs p1 `sin 2 θqq pε s cosθ `aε s ´sin 2 θq 2 .
(2) Dihedral scattering: For the dihedral scattering component, it is parameterized using α (the normalized difference of the combined ground-trunk Fresnel reflection in H and V polarizations) and f d which is the dihedral scattering amplitude [1]: where the Fresnel coefficients are given by: R ih " ε i cosθ i `bε i ´sin 2 θ i and the subscript i denotes ground (g) and trunk (t), as in Equations ( 4) and ( 5), with θ g = π/2 ´θs .In addition, the dihedral scattering amplitude f d is corrected from an attenuation factor (L s ) due to surface roughness [1]: where k is the wave number and s is the RMS height of soil roughness.
(3) Volume scattering: For the volume scattering component, three canonical orientation cases are assumed, namely: vertical, random and horizontal [11].They lead respectively to the following volume coherency matrices [Tv v ], [Tv r ] and [Tv h ] [1]: Finally, using Equation (1), the scattering powers of surface (single bounce) P S , dihedral (double bounce) P D and volume scattering P V are calculated as [1]:

Simplified Polarimetric Decomposition-Based Soil Moisture Retrieval
This paper only focused on the analysis of the surface scattering component (after removing the volume contribution appropriately) to retrieve soil moisture.The idea stemmed from the previous study [1], which indicated that few pixels could be inverted from the dihedral scattering component.Furthermore, the soil moisture retrieval from the dihedral scattering component is based on the two-dimensional minimization in the f d -α plane, since the dihedral scattering component is dependent on both the soil and vegetation dielectric constant.Nevertheless, the dihedral scattering intensity f d is disturbed by the microwave propagation through the vegetation canopy.Therefore, f d must be appropriately compensated for a better retrieval accuracy [1].Otherwise, the retrieved soil moisture from the dihedral scattering component is underestimated.In addition, an ambiguity is presented in the dihedral component around incidence angle 45 ˝to separate the contributions of vegetation and soil.Consequently, this paper aims to simplify the polarimetric-based soil moisture retrieval [1] by using only the surface scattering component, once the volume scattering contribution was removed.Figure 4 presents a flowchart of the algorithm described as follows:

‚
First, the study site is categorized into bare soils and vegetated soils.The X-Bragg model [22] which is an extended Bragg model in terms of roughness conditions is adopted for bare soil retrievals [1].

‚
For the vegetated soils, the orientation of the crop spatial distribution is determined for each pixel according to the Pr = 10¨log10(S vv S vv */S hh S hh *) value [11].

‚
Then, in order to restrict the negative powers of the different scattering components, the theory of non-negative eigenvalues [23] is considered for the ground component by setting eigen (Tg) = 0. ‚ For a given field, only the pixels corresponding to dominant surface scattering (satisfying the condition Re (S hh S vv *) > 0) are selected to estimate the soil moisture.An additional physical constraint ´1 < β < 0 must be satisfied.
‚ Finally, a minimization process is implemented for deriving the soil dielectric constant from β data .Then, the dielectric constant is converted into soil moisture using [24].Therefore, in this study, the unique feature lies in the pixel selection strategy for soil moisture retrieval using the surface scattering component.In this way, the potential of surface component inversion with only one dielectric constant, without any contribution of the dihedral component with two dielectric constants, is investigated.To summarize, for each individual agricultural field, the pixels are selected for soil moisture retrieval if they satisfy simultaneously these two conditions:

‚
The surface scattering is dominant in the remaining ground scattering matrix after removing the volume contribution.

‚
The values of β data extracted from UAVSAR data are physically correct

Assessement of the Volume Scattering Removal Using Entropy and α Angle
Furthermore, it is well known that an important aspect of the abovementioned soil moisture retrieval algorithm over vegetated soils relies on the modeling of the volume component.Thus, in this study, the success in modeling and removing the volume scattering is evaluated by using the Cloude-Pottier decomposition method [20,25].In this decomposition, the scattering entropy (H) and the mean scattering mechanism indicator (α angle) are derived from the eigenvalues λ i (i = 1, 2, 3) and the corresponding eigenvectors e i of the coherency matrix: where the normalized eigenvalue p i and the eigenvector e i are given by p i " λ i and e i " re i1 , e i2 , e i3 s.This eigen-based decomposition algorithm is applied to both the original coherency matrix [T3] and the ground coherency matrix [Tg].In order to evaluate the success in modeling the volume scattering, the 2-D distributions of the investigated fields before and after removing the volume component are compared using the entropy/α plane.

Statistical Index for the Retrieval Process Analysis
Finally, the performances of the retrieval procedure are assessed based on Pearson product moment correlation coefficient and RMSE [26]: where the index i describes the investigated day, O i is the mean value of the measured soil moisture for one crop type on a given day, P i is the mean value of the retrieved soil moisture from the surface scattering component for that crop type, O is the mean value of the measured soil moisture for one crop type during all of the investigated days, P is the mean value of the retrieved soil moisture for one crop type during all of the investigated days, and N is the number of acquisition dates.The R is used to quantify the level of agreement between model predictions and ground measurements, while the RMSE is used to quantify the magnitude of the modeling average error.The next section describes the results obtained from the simplified polarimetric-based soil moisture retrieval and the evaluation of the method.

Results and Discussion
This section presents the results obtained from the application of the simplified polarimetric soil moisture retrieval method to the temporal datasets of UAVSAR acquired during SMAPVEX12.To understand and evaluate the advantages and the limitations of this algorithm, the temporal variations of the scattering mechanisms (before and after the removal of the vegetation volume component) are first analyzed.Then, soil moisture retrieval results are presented and discussed.

Temporal Variation of Scattering Mechanisms
Figure 5 shows the false color composite of the normalized surface, dihedral, and volume scattering mechanisms (Equation ( 7)) for Julian days 169 (17 June 2012), 185 (3 July 2012), and 199 (17 July 2012) corresponding to three representative days (beginning, middle, and end of the campaign).A classification map [18] of the considered five crop types is provided to relate the scattering characteristics to crop types (on this map, the areas in white color are covered by other crops and the forested site).Generally, it can be noted that the surface scattering intensity tends to decrease from June to July 2012, since the fast crop development during the investigation period and the well-developed crops' structure (height and biomass are increased) in July reduce the penetration depth of the microwave signal.This decrease of the surface scattering component is especially pronounced in the near range, as the single scattering component is more significant at low incidence angle (corresponding to near range).A strong dihedral scattering component, like in [1], is observed at the beginning of the campaign (Figure 5a) for several fields, mainly those covered by corn.Indeed, according to the SMAPVEX12 ground measurements, the vegetation height and biomass are still low at this time, leading to a double bounce between vegetation stalk and the underlying soil.At mid-campaign, the dihedral scattering component becomes less significant due to the developed leaf layer and the increased complexity of vegetation structure which reduce the dihedral scattering mechanism.The volume scattering component is predominant over the fields covered by canola, wheat and corn.For different crop types, Figure 6 shows the temporal variation of each scattering mechanism (surface, dihedral, and volume scatterings) derived as the mean value of the normalized scattering power.The backscattering mechanisms of soybean are in accordance with the temporal variation of plant penology (Figure 3), which indicates the increase in volume scattering power and the decrease in surface scattering power.As the vegetation water content of soybean is low during the developing stage considered in this study (less than 1 kg/m 2 ), the temporal evolution of its volume scattering is mainly affected by the changes in soybean structure.Regarding the corn, Figure 6b shows that its dihedral scattering power is significant at the beginning of the campaign as a result of the already developed corn stalk supporting only a few leaves.From Julian day 175 (June 23), the dihedral scattering component decreases due to the well-developed corn leaves which attenuate the electromagnetic waves (followed by a decline of surface scattering component from Julian day 177 to 187).For the pasture, the temporal evolution of scattering mechanisms is quiet stable, which is probably due to the low and stable biomass indicated by the ground measurements.For wheat fields, the increase in surface scattering (Figure 6a) is in agreement with the quasi stable wheat's water content value (2.2 kg /m 2 ) and its increasing height [1].This corresponds to a vegetation drying out, which enhances the surface scattering and decreases the volume scattering (Figure 6c) due to a more pronounced signal penetration.
For the fields covered by canola, the surface scattering power increases while the volume scattering power decreases during the growing period (Figure 6).These temporal behaviors which seem unexpected are closely related to the temporal evolution of height, biomass and vegetation water content (Figure 3) as well as to the correlation between these parameters.Indeed, during the crop growth cycle, canola's height, biomass and vegetation water content increase until reaching the mature stage which is then followed by a drying out process.Therefore, the decreasing vegetation water content from the maturation stage leads to a decrease in the volume scattering and an increase in the surface scattering.Consequently, due to the combined effect of the biomass, height and the vegetation water content on the volume scattering, the latter does not increase during the growing season as expected; in turn, neither does the surface scattering during the growing season.In addition, other reasons arising from the polarimetric decomposition can explain these results [1,2], particularly those obtained before the canola's drying out.For instance, the high and sophisticated temporal changes in the shapes and structures of canola during the considered phenological stage may not be fully exploited by the only three fixed volume coherency matrices.The scattering characteristics are affected by several factors including the vegetation water content, height, biomass and density.The scatterplots in Figures 7 and 8 present the behaviors of the volume and surface scatterings with respect to vegetation water content.They show that as the vegetation water content increases, the volume scattering power increases while the surface scattering power decreases respectively.Different volume/surface scattering behaviors are obtained even for similar vegetation water content values due to the combined effects of the influencing factors.Furthermore, the dynamic of scattering mechanisms to vegetation water content is significantly dependent on the crop growing stage.The low and no correlations observed between the vegetation water content and volume (or surface scattering) during the middle and the end of the campaign (Figure 7b,c, Figure 8b,c), respectively, suggest the loss of explanatory power of the adopted linear regression.

Evaluation of the Modeling of the Volume Scattering and Its Removal from the Measured Coherency Matrix
In this paper, a main challenge lies in the modeling of the vegetation contribution and its removal in order to better retrieve the soil moisture from the remaining surface scattering component.To understand the vegetation scattering, the polarimetric parameters entropy and α scattering angle (Section 3.3) are used to build a 2-D plane, which is subsequently divided into nine zones for different scattering characteristics.For comparison purposes, for several dates, Figures 9 and 10 show the results of 2-D statistical distribution density of crops obtained before and after the removal of the volume component respectively, with an over plot of the mean values of entropy and α angle for each test field.The color represents the statistical distribution density (from low density in blue to high density in red), while the dashed lines denote the boundary of the nine divided scattering characteristics.Through the comparison, the effect of volume contribution can be understood.
Figure 9 shows that soybean fields present the highest dynamic, especially for its entropy values.It can be seen in Figure 9a, which corresponds to the beginning of campaign, that the soybean fields are dominated by the surface scattering (α angle is less than 45 ˝), and they are well discriminated from the other crop types.With the growth of the soybean, the entropy and the α angle increase due to changes in crop structure (Figure 9b).At the end of the campaign, when the crops are in the mature stage, Figure 9c shows that the field distribution is more clustered, thus hampering any discrimination between the crops.Comparing Figures 9 and 10 a decrease in entropy values clearly appears for all crops.This result is expected since the remaining ground coherency matrix [Tg] used to generate the Figure 10 represents the surface and dihedral scatterings, and consequently provides lower entropy than the results (Figure 9) obtained from the measured coherency matrix [T3].Nevertheless, in Figure 10, some of the scattering α angles which are still located within the dipole/volume scattering region are probably related to the residual volume scatterings.These can be explained not only from the limitation of this polarimetric decomposition algorithm using fixed volume orientation cases, but also from the strict narrow criterion on α angle.Indeed, α angle used to characterize the scattering mechanisms stems from a purely empirical segmentation of the entropy/α plane [20,25].
Furthermore, the performance of the volume removal depends on both crop phenological stage and crop types due to different plant geometry and density.Better results were obtained for moderate vegetation height and biomass condition.For example, during mid-campaign, the peak value of 2-D distribution density is shifted significantly into the surface scattering region in the entropy/α plane (Figure 10b).In contrast, at the beginning and the end of the campaign, the volume contribution can only be partially removed.Indeed, in this study, the vegetation structure is modeled as dipoles considering three fixed orientation cases.However, the low vegetation (lacking strict stem) and the high vegetation status (with abundant leaves and ears) seem to be represented by a spherical-like shape [1].As for the effect of crop types on the performance of the volume scattering modeling and removal, good results are obtained for soybean and pasture.We expect therefore that these crops match better with the predefined and fixed volume coherency matrix.Nevertheless, the volume contribution of wheat and canola can only be partially removed by the proposed volume model (Equation ( 6)).It is likely that their orientations are not sufficiently represented by the only three fixed volume components.In order to further analyze the effect of the three fixed vegetation orientation cases on the removal of the volume component, the exemplary results obtained over canola fields on day 185 (3 July 2012) are compared through the entropy/α plane presented in Figure 11.For this crop and at this phenological stage particularly, even though there is no clear tendency in the results shown in Figure 11, vertical orientation seems to provide the best one.Therefore, any failure in considering the appropriate crop orientation will negatively impact the modeling and the removal of the volume component.

Soil Moisture Retrieval Results
They are obtained from the fields characterized by a dominant surface scattering component, after the removal of the volume scattering component.
(1) Percentage of dominant surface scattering after removing the volume component: In order to ensure both the retrieval accuracy and spatial inversion rate, a given field was selected for soil moisture inversion, if more than 10% of its area was covered by pixels satisfying the dominant surface scattering criteria.Since for canola, corn and wheat's fields the dominant surface scattering cases were low at the beginning of the campaign (Figure 12), thresholding the dominant surface scattering cases at 10% allowed retrieval of the soil moisture for almost all the crops.However, as the soil moisture retrieval algorithm is a simplified one, which is only based on the dominant surface scattering, the accuracy is compromised in the conditions of the highly dominant dihedral scattering component.Consequently, for these fields characterized by more than 10% of dominant surface scattering cases, a spatial average of soil moisture was computed from the retrieved soil moisture obtained over the pixels with the dominant surface scattering.Considering Figure 12 which presents the percentage of the dominant surface scattering case for different crop types, the surface scattering component could be used as inputs in the retrieval process, except for canola fields before Julian day 173 (June 21), when the dominant surface scattering case is lower than 10%.Due to the complementarity between the dominant surface scattering and the dihedral scattering component in the remaining signal, Figure 12 also indicates an important dihedral scattering component for canola, wheat and corn from Julian day 169 to 178 (from June 16 to 26).Under such conditions, less performance of the simplified polarimetric soil moisture retrieval algorithm will be obtained, which is mainly due to the over-subtraction of volume scattering contribution.(2) Evaluation of the simulated β parameter: This step is important since the soil moisture is derived from the normalized difference of Bragg scattering in H and V polarization β (Figure 4). Figure 13 compares the β data values which are derived from UAVSAR data with corresponding simulation values obtained using the known incidence angles of UAVSAR and the field average of the measured soil moisture as inputs.In Figure 13, each β data value is a field average for a given sampling field (it is computed by averaging the values of all the pixels within the targeted field, except the pixels near the field boundary to remove boundary effects).Disagreements are found between the simulated and the derived β parameter for canola, wheat and corn crops, while for soybean and pasture the results of β data and β simu match better.This result is in line with the analysis of the scattering mechanisms conducted in Section 4.2, specifically the limitations observed in the volume scattering modeling for canola and wheat and its good performance for soybean and pasture.However, for soybean and pasture, the weak deviation observed between the simulated and the derived values of β confirms the ability of the surface and volume scattering formulations, as described in [1], to model the phenological development of short vegetation such as soybean and pasture.Therefore, it is important to highlight that the applicability of the simplified polarimetric soil moisture retrieval algorithm is dependent on both the crop type and the crop development stage.(3) Analysis of the retrieved soil moisture from the surface scattering component: Figure 14 presents the maps of the retrieved soil moisture for the beginning, middle and end of the agricultural campaign, while the invalid pixels are to white color.Two main reasons can explain these invalid pixels: First, some pixels are characterized by a dominant dihedral scattering after the removal of the volume component.Indeed, in this study, the dihedral component is disregarded due to its ambiguity and disturbance by the wave propagation (the dihedral intensity f d needs to be compensated properly [1,2]) as well as the low improvement in the soil moisture inversion rate that could be obtained.Second, for some pixels characterized by a dominant surface scattering, the corresponding β value is less than ´1, which is non-physically correct [1].However, the range (5% to 50%) of the retrieved soil moisture observed in Figure 14 is in agreement with that of the in situ soil moisture measurements [18].Furthermore, while some inversion results are provided as field-average values [5], it is interesting to underline the ability of this algorithm to inform about the spatial variability in soil moisture at the field scale (Figure 14).As for the retrieval rate in Figure 15, it varies from 26% to 38% before the mature stage of the crops and then decreases, thereby leading to high amount of white pixels in Figure 14c.A analysis is also conducted by comparing the mean values of the retrieved soil moisture on each field with the in situ soil moisture measured during the SMAPVEX12 experiment.However, in agreement with the scattering mechanism analysis (Section 4.2) and β simulation values (Figure 13), Figure 16 shows that better results are obtained over the fields covered by soybean and pasture than those covered by other crop types, especially at the beginning of the campaign.Besides, Figure 16 shows that the canola fields are not invertible for some dates.On Julian day 199, the inversion succeeded over three canola fields (Figure 16c) with an RMSE value around 0.06 m 3 /m 3 .For the corn fields, the accuracy of most of the retrieved results is less than ˘0.1 m 3 /m 3 .As for the wheat fields, it should be noted that the in situ measurements are missed on Julian day 185, thus the comparison with the retrieved results cannot be performed (Figure 16b).Figure 16a,c show that over the wheat fields the results are underestimated mainly due to the bias of β values (Figure 13), which stems from the partial removal of the volume scattering component (Figure 10).

Validation
In Figure 17, for each crop type, the time series of the mean value of the retrieved soil moisture are presented, along with that of in situ soil moisture measurements:

‚
For the canola (Figure 17a), soil moisture cannot be inverted for the most part of the phenological period (around 80%).Even if after removing the volume scattering contribution the surface scattering component is dominant over the canola fields (Figure 12), the second condition that its β data values should be confined between ´1 and 0 is generally not respected.Therefore, the obtained low values of β data (less than ´1) are considered as non-theoretically correct for soil moisture inversion.These invalid low values of the extracted β data are a consequence of the limitations observed in the modeling of the volume scattering of canola fields.However, although the soil moisture can be retrieved on some days (e.g.around the days of June 23 and July 14) with a RMSE of 0.076 m 3 /m 3 and a correlation of 0.62, the low inversion rate (20 %) limits the application of this algorithm to canola crops.

‚
For the corn (Figure 17b), the temporal evolution of the measured and the retrieved soil moisture agrees well up to Julian day 190 (July 8th).According to the temporal evolution of vegetation characteristics presented in Figure 3, high height (more than 150 cm), biomass (more than 3.5 kg/m 2 ), vegetation water content (more than 2.8 kg /m 2 ) and fast changes in corn structure can explain the higher difference observed between the measured and the retrieved soil moisture after Julian day 190.

‚
For the pasture (Figure 17c), RMSE of 0.098 m 3 /m 3 is obtained.However, there is no correlation between the temporal variation of the retrieved and the measured soil moisture.A possible reason may be the large std values of the measured soil moisture over pasture fields (Figure 16).Besides, compared to the agricultural fields, the pasture characteristics are more heterogeneous.

‚
For the soybean (Figure 17d), the retrieval is achieved with a correlation of 0.6 and RMSE of 0.078 m 3 /m 3 .The temporal variation of the three scattering mechanisms of soybean (Figure 6) is in agreement with its phenological development (Figure 3), and the simulations and the extracted values of β matched (Figure 13).Thus, the three predefined fixed vegetation orientation cases are considered suitable to model the vegetation contribution for soil moisture inversion from the full polarimetric signature of soybean fields.

‚
For the wheat (Figure 17e), the inversion is only achieved for 70% of the investigated period, with a correlation of 0.66 and RMSE of 0.12 m 3 /m 3 .Using single polarization ENVISAT ASAR data normalized at 20 ˝, higher retrieval accuracies (RMSE of 5.3% vol.and 6.4% vol., for respectively HH and VV) were obtained by [27] over wheat fields characterized by vegetation water content ranging from 0.08-1.15kg/m 2 .Compared to our dataset, both the lower incidence and the lower vegetation water content are very helpful for improving the retrieval accuracy.
With RMSE values varying from m 3 /m 3 , Figure 17 shows that the performance of the retrieval algorithm is crop type and phenological stage dependent.The soil moisture is underestimated over most of the crop fields, except for the corn.This result can be explained by the dipole modeling of crop shape [1] and the predefined fixed volume coherency matrices (only three types).Indeed, it is reported in [1] that assuming a dipole particle shape will lead to an underestimation of soil moisture.In [1], the soil moisture retrieval over corn fields is the best with RMSE of 0.08 m 3 /m 3 during the temporal comparison with ground measurements.For the fields covered by canola and wheat, the RMSEs are around 0.10-0.14m 3 /m 3 .In this study, by using only the surface scattering mechanism and neglecting the dihedral scattering mechanism which is affected by both the characteristics of vegetation and the underlying soils, the obtained RMSEs are 0.064 m 3 /m 3 for corn, 0.076 m 3 /m 3 for canola and 0.12 m 3 /m 3 for wheat.Nevertheless, in terms of performances, a comparison cannot be made between the soil moisture retrievals conducted in this study and in [1], due to the difference in the experimental conditions including the location of test areas, the crop growth stages and the soil and vegetation characteristics.

5.
This study investigates a simplification of the polarimetric decomposition algorithm of [1] for soil moisture retrieval over different agricultural fields of the SMAPVEX12 site by using only the surface scattering component.In comparison to [1], the volume scattering intensity is corrected by the eigen-based algorithm in order to restrict the negative power in the decomposed scattering mechanisms.Furthermore, to our knowledge, it is one of the first studies where a large dataset of multi-temporal SAR airborne measurements acquired over several vegetation covers (canola, corn, pasture, soybean and wheat) is used in order to discuss the potential of an L-band soil moisture decomposition algorithm for fields with different vegetation characteristics and a high dynamic of soil moisture.
Although a simplified soil moisture retrieval procedure is used, the obtained inversion rate of 26%-38% is higher than that of 5%-25% in [1], but still lower than that of 30%-49% in [2] for single acquisition.Difference in agricultural field conditions and in UAVSAR acquisitions could be an explanation.Furthermore, the simplification made by using only the dominant surface scattering component allows the soil moisture estimation, without any ambiguity that could be encountered by integrating the dihedral component in the retrieval process.
The main limitation of this simplification consists of the non-applicability of this simplified algorithm when the dihedral scattering mechanism is dominant.For canola fields, this occurs at the beginning of SMAPVEX12 campaign.At this stage, less than 10% of the pixels within the canola field are characterized by dominant surface scattering.Nevertheless, with the vegetation growth, the dihedral component decreases.However, the algorithm can retrieve the soil moisture only if the extracted β values are physically correct as a result of appropriate removal of the volume scattering component.This second condition allows the inversion from the proposed physically based algorithm only if the field conditions are in agreement with the scattering processes.The discontinuities observed in Figure 17a,e express the violation of this condition.
In addition to the crop phenology, it clearly appears that the performance of the soil moisture retrieval algorithm also depends on the crop types.The algorithms seem to be more suitable for crop types with a structure similar to soybean (short height and low biomass), as the volume scattering contribution can be properly modeled and removed.However, as shown in Figure 1, the investigated fields are distributed along the range.Consequently, the retrieval accuracy is affected by the incidence angle varying from the near range (25 ˝) to the far range (65 ˝).Even if a given field is observed with the same incidence angle from one acquisition date to another, the fields of a given crop type are observed with different incidence angles.These incidence angle variations are not taken into account in the evaluation of the retrieval accuracy, since the main objective of this study is to analyze the potential of polarimetric decomposition for soil moisture retrieval under various vegetation characteristics.Comparing our results with those from a preliminary inversion using the same UAVSAR dataset and two SMAP candidate radar-only algorithms which were applied to a few fields of the SMAPVEX12 campaign [19], the accuracies of the soil moisture retrieval are similar.
Due to the difficulty to account for the real complexity of crop structure during the vegetation growing period, this study also confirms that the modeling of volume scattering contribution is the major challenge for this kind of algorithm.A solution for removing a varying volume scattering is recently proposed and consists of modifying the particle shape and orientation for each resolution cell to individually adapt the vegetation structure [3].However, such modeling of the general volume scattering may require more parameters to describe the shape, orientation and phenological stage of crop.It might lead to an optimization for volume scattering subtraction instead of direct inversion of soil moisture [3].Ongoing work is specially focused on improving the vegetation volume scattering component for a better soil moisture retrieval over vegetated agricultural fields.

Figure 1 .
Figure 1.The SMAPVEX12 study site with delineation of UAVSAR swath (gray color image of T 11 in dB) and corresponding land cover.

Figure 2 .
Figure 2. Temporal evolution of the measured volumetric soil moisture along with daily precipitation amount and the availability of UAVSAR acquisitions (black arrows).

Figure 3 .
Figure 3. Temporal variation of vegetation (a) height, (b) wet biomass and (c) vegetation water content.

Figure 4 .
Figure 4.The simplified polarimetric-based soil moisture retrieval over agricultural fields.The dashed box indicates the process for removing the volume scattering component.

Figure 5 .
Figure 5. RGB color composition of the normalized three scattering mechanisms for Julian day (a) 169 (17 June 2012); (b) 185 (3 July 2012); (c) 199 (17 July 2012).Dihedral scattering power in red, volume scattering power in green and surface scattering power in blue.The incidence angle in the range direction varies from 25 ˝to 65 ˝.(d) Classification map of the considered five crop types (the areas in white color are covered by other crops and the forested site).

Figure 6 .
Figure 6.Temporal variation of normalized scattering powers of (a) surface scattering, (b) dihedral scattering and (c) volume scattering component for different crop types.The points (corresponding to right y-axis) shows the height of different crop types.

Figure 11 .
Figure 11.Distribution of canola fields on entropy/α plane after the removal of volume scattering component on Julian day 185 (3 July 2012) for three vegetation orientations (vertical, horizontal and random).

Figure 12 .
Figure 12.Percentage of dominant surface scattering case after removing volume scattering component, for different crop types.

Figure 15 .
Figure 15.The retrieval rate by using surface scattering component.

Figure 17 .
Figure 17.Comparison between the retrieved and measured soil moisture during the agricultural campaign for (a) Canola; (b) Corn; (c) Pasture; (d) Soybean; (e) Wheat.The retrieved and measured soil moisture are averaged for each crop type, and the discontinuities in the curves correspond to unsuccessful retrieval periods.
*: no corresponding ground truth measurement (use the field data of 17 June).

Crops RMS Height (cm) Correlation Length (cm) Min Mean Max Min Mean Max
[2] resulted three values of volume scattering intensity (fv 3 , fv 4 , and fv 5 ) are compared to the previous estimated values (fv 1 , fv 2 ).The appropriate estimation of volume scattering intensity is considered as the minimum value of these five calculated volume scattering intensities (fv 1 , fv 2 , fv 3 , fv 4 , and fv 5 )[2].
‚Ground scattering component [Tg] is then obtained by subtracting the appropriate volume component from the measured coherency matrix.