Theoretical Evaluation of Water Cloud Model Vegetation Parameters

: The advantage of implementing the Water Cloud Model (WCM) is in being able to express complex scattering characteristics in a vegetated area with simple bulk vegetation descriptors. However, there has been a lack of understanding or consensus about the optimal set of vegetation descriptors. In this paper, the original and improved expressions of WCM are evaluated and the optimal vegetation descriptors are presented by examining the relationship between WCM vegetation parameters and the theoretical scattering model predictions. In addition, the condition-speciﬁc regression relationship between bulk vegetation descriptors and theoretical scattering and attenuation coe ﬃ cients, expressed by the A and B parameters in the WCM, is analyzed in relation to the shape, size, and orientation distribution of the scatterer. Furthermore, the inﬂuence of radar observation conditions on the parameterization of the WCM is presented. The results show that the particle moisture content and the vegetation water content can be the optimal vegetation descriptors, denoted by the V 1 and V 2 variables in the WCM, respectively.


Introduction
Due to its all-weather imaging and vegetation penetration capabilities, remote sensing of vegetated areas with Synthetic Aperture Radar (SAR) has a great potential for retrieving bio-and geo-physical parameters related to vegetation and the underlying soil surface. The measured SAR signal backscattered from natural targets is given by the superposition of many elementary scatterers. Consequently, in order to estimate physical properties of scatterers from SAR data, it is necessary to investigate different scattering contributions among total backscattered signals by modeling the interactions of microwaves in the vegetated areas.
Theoretical scattering models, such as the discrete scattering model [1,2] and the radiative transfer model [3][4][5], have been used to assess the relative importance of microwave scattering mechanisms. In the theoretical models, the vegetation canopy is, in general, represented as an ensemble of randomly distributed dielectric particles, in a layer of specific height, overlying a dielectric ground layer. The total backscatter can be obtained by an incoherent sum of several scattering contributions, including direct backscattering from the vegetation layer, direct backscattering from the underlying rough surface, scattering interaction between the vegetation and the ground surface, and ground-vegetation-ground multiple bounce. The theoretical scattering models have been used successfully to interpret scattering characteristics of vegetated areas and to predict radar signals in relation to the biophysical properties of plants. However, they are usually complex and inconvenient because of a large number of input parameters which lead to complexity in resolving the inverse problem.

Water Cloud Model
The Water Cloud Model represents a vegetation layer as a collection of identical spherical particles, uniformly distributed throughout the volume layer. The total backscattering coefficient, σ 0 , can be expressed by the incoherent sum of the scattering contribution of the vegetation layer, σ 0 v , and the scattering contribution of the soil layer, σ 0 s , attenuated by the attenuation coefficient T 2 , given as follows: where σ v and κ e are the volume backscattering and extinction coefficients of the vegetation layer, h is the vegetation height, θ is the incidence angle, and σ 0 s is the soil backscattering coefficient. Neglecting multiple scattering effects, σ v and κ e , can be given by the sum of the backscattering cross section, S, and the extinction cross section, k, of a single particle, such as Remote Sens. 2019, 11, 894 3 of 17 where N is the number of particles per unit volume (m −3 ). The simplification of the vegetation layer permits the relation of the scattering characteristics to bulk vegetation parameters. In the case of uniformly distributed water particles, the extinction coefficient can be assumed to be proportional to the total water content in the unit volume W (kg/m 3 ), such as [6] κ e = B 1 W.
In a similar way, the ratio σ v /2κ e , which determines the vegetation backscattering coefficient, can be assumed to be independent of water content, such as However, the simple expressions for the scattering and extinction coefficients in Equations (5) and (6) are only valid for the cloud of water particles. In order to apply the simple vegetation model to practical SAR observations, many researches have used more general expression [7][8][9][10][11][12][13][14][15][16], such as The unknown coefficients A 2 and B 2 in the general WCM can be dependent on the canopy type. Since there have been no theoretical basis to determine A 2 and B 2 , they need to be empirically determined. The scattering and extinction coefficients are represented by bulk vegetation parameters V 1 and V 2 . Due to the heterogeneity of vegetation structure, several different sets of vegetation parameters have been proposed in the literature, e.g., VWC, LAI, canopy height (h), and particle moisture content (m g ).
In order to examine the best set of vegetation parameters and the corresponding A 2 and B 2 , let us firstly assume a vegetation layer consisting of cloud of spherical dielectric particles. The backscattering and extinction cross sections for a spherical Rayleigh scatterer of radius, r, and dielectric constant, ε v , are [18] Here, the dielectric constant of a particle, ε v , can be related to the gravimetric moisture content, m g , through the Ulaby and El-Rayes model [19], where W w and W d are the wet and dry mass of a particle. By comparing Equations (1)-(3) and (7), the A 2 V 1 and the B 2 V 2 terms can be expressed as follows: To understand the coefficients A 2 and B 2 and the vegetation parameters V 1 and V 2 , variations of A 2 V 1 and B 2 V 2 terms are examined in relation to several vegetation descriptors, such as the particle moisture content, m g , the areal density, N a (m −2 ), and the VWC (kg/m 2 ), as shown in Figures 1 and 2. The areal density and VWC are defined as where ρ w is the particle wet density (∼ 500 kg/m 3 , [20]) and particle volume v p = (4π/3)r 3 . Figures 1 and 2 show calculated A 2 V 1 and B 2 V 2 parameters plotted against different vegetation descriptors. The dielectric constant, ε v , for calculating WCM parameters is derived from the Ulaby and EL-Rayes model at 5 GHz. It is shown that the A 2 V 1 term is highly related to m g , while it is independent of particle density. Since it exhibits a quadratic-like dependence on m g , the A 2 V 1 term can be rewritten in relation to m g , such as Here, the coefficient A 2 varies with the particle size. According to the least square method, A 2 varies from 0.021 for r = λ/50 to 2.6 for r = λ/10, while E = 1.9 independently with the particle size. On the contrary, the B 2 V 2 term exhibits a linear relationship with VWC, as shown in Figure 2. In this case, we can rewrite the B 2 V 2 term as The coefficient B 2 is not a constant but determined by the moisture content of particles. It varies from 0.27 for dry particles (m g = 0.1) to 0.014 for wet particles (m g = 0.9).

Figures 1 and 2 show calculated
and parameters plotted against different vegetation descriptors. The dielectric constant, , for calculating WCM parameters is derived from the Ulaby and EL-Rayes model at 5 GHz. It is shown that the term is highly related to , while it is independent of particle density. Since it exhibits a quadratic-like dependence on , the term can be rewritten in relation to , such as Here, the coefficient varies with the particle size. According to the least square method, varies from 0.021 for = /50 to 2.6 for = /10, while = 1.9 independently with the particle size. On the contrary, the term exhibits a linear relationship with VWC, as shown in Figure 2. In this case, we can rewrite the term as The coefficient is not a constant but determined by the moisture content of particles. It varies from 0.27 for dry particles ( = 0.1) to 0.014 for wet particles ( = 0.9).   where is the particle wet density (~500 kg/m 3 , [20]) and particle volume = (4 /3) .

Figures 1 and 2 show calculated
and parameters plotted against different vegetation descriptors. The dielectric constant, , for calculating WCM parameters is derived from the Ulaby and EL-Rayes model at 5 GHz. It is shown that the term is highly related to , while it is independent of particle density. Since it exhibits a quadratic-like dependence on , the term can be rewritten in relation to , such as Here, the coefficient varies with the particle size. According to the least square method, varies from 0.021 for = /50 to 2.6 for = /10, while = 1.9 independently with the particle size. On the contrary, the term exhibits a linear relationship with VWC, as shown in Figure 2. In this case, we can rewrite the term as = VWC = 6 10 Im − 1 + 2 1 VWC.
The coefficient is not a constant but determined by the moisture content of particles. It varies from 0.27 for dry particles ( = 0.1) to 0.014 for wet particles ( = 0.9).

WCM Parameters for Non-Spherical Particles
The WCM tries to express complicated canopy backscatters in terms of a small number of vegetation parameters. However, the WCM with identical spherical particles predicts identical co-polarization and no cross-polarization backscatters, which are far from the real SAR observables. In order to figure out more realistic WCM model parameters, backscatter simulations with physical models considering non-spherical particles are presented in this section.
Based on the solution of the radiative transfer equation [3][4][5], the pq-polarization backscattering contribution of the vegetation layer can be written as In the case of a canopy layer consisting of randomly oriented dielectric non-spherical particles assuming identical size and shape, the volume scattering coefficient, σ v,pq , and the extinction coefficient, κ e,pq , can be obtained by statistical averaging over the scatterer's orientation distribution p(α, β), such as The orientation distribution of the azimuthal angle, α, and the inclination angle, β, are shown in Figure 3a. All scatterers are assumed to be distributed uniformly in the azimuth, such as p(α) = 1/2π. For a small non-spherical particle, the scattering amplitude, S pq (ŝ,î), on the incident,î, and scattered, s, directions and the extinction cross section, k p (î), can be calculated by using the Generalized Rayleigh-Gans (GRG) approximation [21]. It has been widely used (e.g., [22][23][24][25]) for modeling radar scattering properties from leaves. To evaluate the effects of particle shape, size, and orientation on the parameterization of WCM, two types of vegetation, composed of needle shaped ( Figure 3b) and circular disk shaped (Figure 3c) particles, are considered in this study.

WCM Parameters for Non-Spherical Particles
The WCM tries to express complicated canopy backscatters in terms of a small number of vegetation parameters. However, the WCM with identical spherical particles predicts identical copolarization and no cross-polarization backscatters, which are far from the real SAR observables. In order to figure out more realistic WCM model parameters, backscatter simulations with physical models considering non-spherical particles are presented in this section.
Based on the solution of the radiative transfer equation [3][4][5], the pq-polarization backscattering contribution of the vegetation layer can be written as In the case of a canopy layer consisting of randomly oriented dielectric non-spherical particles assuming identical size and shape, the volume scattering coefficient, , , and the extinction coefficient, , , can be obtained by statistical averaging over the scatterer's orientation distribution ( , ), such as The orientation distribution of the azimuthal angle, , and the inclination angle, , are shown in Figure 3a. All scatterers are assumed to be distributed uniformly in the azimuth, such as ( ) = 1/2 . For a small non-spherical particle, the scattering amplitude, (,), on the incident, ̂ , and scattered, ̂ ,directions and the extinction cross section, () , can be calculated by using the Generalized Rayleigh-Gans (GRG) approximation [21]. It has been widely used (e.g., [22][23][24][25]) for modeling radar scattering properties from leaves. To evaluate the effects of particle shape, size, and orientation on the parameterization of WCM, two types of vegetation, composed of needle shaped ( Figure 3b) and circular disk shaped ( Figure 3c) particles, are considered in this study.

Effect of Particle Shape
For non-spherical particles, the constants and will vary with polarization. In this study, polarization dependent model parameters will be examined for co-polarization cases, since the crosspolarization signal generally contains no scattering contribution from the soil surface, which has been of great interest in WCM applications. The polarization dependent model constants , and , can be estimated by fitting WCM parameters to the simulation results of the GRG-based physical

Effect of Particle Shape
For non-spherical particles, the constants A 2 and B 2 will vary with polarization. In this study, polarization dependent model parameters will be examined for co-polarization cases, since the cross-polarization signal generally contains no scattering contribution from the soil surface, which has been of great interest in WCM applications. The polarization dependent model constants A 2,pq and B 2,pq can be estimated by fitting WCM parameters to the simulation results of the GRG-based physical model. The WCM components are related to the vegetation attenuation and scattering terms of the physical model as The volume scattering and the extinction coefficients in the left-hand side of the equation are calculated by the GRG approximation. A number of scattering and extinction coefficients of different vegetation, composed of either needles or disks, are simulated using the parameters given in Table 1. In the simulation, the thickness (t) of needles and disks are set to be 0.1 cm and 0.01 cm, respectively. In addition to the uniform azimuth angle distribution, all scatterers are assumed to be distributed uniformly in the inclination angles, in the calculation of the GRG-based model.  Figure 4 illustrates variations of κ e,HH h and κ e,VV h terms by the GRG-based model with respect to different vegetation descriptors, such as m g , N a , and VWC. The simulation is performed at 5 GHz frequency. It shows a strong linear relationship between κ e,pq h terms and VWC in both HH-and VV-polarizations, regardless of the particle shape. Consequently, the VWC is selected as the optimal V 2 parameter and the corresponding B 2,pq coefficient can be determined by the slope of the least square line. Results show that the particle shape affects significantly on the B 2,pq coefficients of the WCM. A vegetation canopy composed of disk-shaped particles is described by a higher B 2,pq coefficient than if it is composed of needle shaped particles. In addition, despite the uniform distribution assumption, the estimated B 2,HH is slightly higher than B 2,VV in the case of disk-shaped particles. On the other hand, in the case of needle shaped particles, the estimated B 2,VV is slightly higher than B 2,HH .
To determine the V 1 parameter and to estimate corresponding A 2,pq coefficients, variations of the σ v,pq /2κ e,pq terms are examined against different vegetation descriptors, as shown in Figure 5. As with the spherical particle, the σ v,pq /2κ e,pq is strongly correlated with the m g and exhibits a quadratic-like relationship. Based on a non-linear expression, such as the A 2,pq and E parameters can be estimated by least square fit. The estimated E is about 1.9 independently with the particle shape, size, and density. The estimated model coefficient A 2,pq is affected by particle shape. The disk shaped particle exhibits higher A 2,pq than the needle shaped scatter.
The model coefficients at HH-and VV-polarizations, A 2,HH and A 2,VV , are not significantly different in both the disk and needle shaped particles.    at 5 GHz frequency. It shows a strong linear relationship between , ℎ terms and VWC in both HH-and VV-polarizations, regardless of the particle shape. Consequently, the VWC is selected as the optimal parameter and the corresponding , coefficient can be determined by the slope of the least square line. Results show that the particle shape affects significantly on the , coefficients of the WCM. A vegetation canopy composed of disk-shaped particles is described by a higher , coefficient than if it is composed of needle shaped particles. In addition, despite the uniform distribution assumption, the estimated , is slightly higher than , in the case of disk-shaped particles. On the other hand, in the case of needle shaped particles, the estimated , is slightly higher than , . To determine the parameter and to estimate corresponding , coefficients, variations of the , /2 , terms are examined against different vegetation descriptors, as shown in Figure 5. As with the spherical particle, the , /2 , is strongly correlated with the and exhibits a quadratic-like relationship. Based on a non-linear expression, such as

Effect of Particle Orientation
In the previous experiment, non-spherical particles are assumed to be uniformly distributed throughout vegetation layer. However, in many agricultural crops and forests, the scatterers in a volume layer may have a preferred orientation distribution. In order to consider more realistic vegetation structure in the parameterization of WCM, scattering and attenuation coefficients from particles with preferred orientation distribution are considered in this part. The distribution of the inclination angle p(β) is used to denote the orientation of particles. For simplicity, it is assumed to be a uniform distribution over the interval β 1 ≤ β ≤ β 2 , such as P(β) = 1/(β 2 − β 1 ). When β varies in a range around 0 • , the needle shaped particles are vertically oriented, while the disk-shaped particles are nearly horizontal. Figure 6 shows variations of the B 2,pq coefficients in the oriented volume. Based on previous analysis, VWC is selected for the vegetation descriptor V 2 and B 2,pq is obtained by fitting the linear function. For the inclination distribution, we consider the three following special cases: 0

Effect of Particle Orientation
In the previous experiment, non-spherical particles are assumed to be uniformly distributed throughout vegetation layer. However, in many agricultural crops and forests, the scatterers in a volume layer may have a preferred orientation distribution. In order to consider more realistic vegetation structure in the parameterization of WCM, scattering and attenuation coefficients from particles with preferred orientation distribution are considered in this part. The distribution of the inclination angle ( ) is used to denote the orientation of particles. For simplicity, it is assumed to be a uniform distribution over the interval ≤ ≤ , such as ( )=1/( − ). When varies in a range around 0°, the needle shaped particles are vertically oriented, while the disk-shaped particles are nearly horizontal. Figure 6 shows variations of the , coefficients in the oriented volume. Based on previous analysis, VWC is selected for the vegetation descriptor and , is obtained by fitting the linear function. For the inclination distribution, we consider the three following special cases: 0°≤ ≤90°, 0°≤ ≤60°, and 0°≤ ≤30°. The first case corresponds to the uniform inclination angle distribution, as with the previous experiment. The third case corresponds to either nearly vertically oriented needle or nearly horizontally oriented disk particles. It is seen from Figure 6 that the linear relationship between , ℎ terms and VWC is maintained throughout different orientation distributions. However, for horizontal disks and vertical It is seen from Figure 6 that the linear relationship between κ e,pq h terms and VWC is maintained throughout different orientation distributions. However, for horizontal disks and vertical needles, oriented volume effects in wave propagation through vegetation layer lead to the differences of the B 2,pq coefficients between HH-and VV-polarizations. Particularly, the B 2,HH coefficient, which is associated to the attenuation of HH-polarized signal, is more affected by the orientation of the scatterer than the B 2,VV coefficient. In comparing needle and disk-shaped particles, the orientation angle effect is more significant in the needle shaped particle. The B 2,HH coefficient decreases significantly as the needle shaped particle is oriented vertically.
Next, the effect of particle orientation on the A 2,pq coefficients is evaluated in Figure 7. As in the previous analysis, m E g is used for the vegetation descriptor. As with the B 2,pq coefficients, orientation distribution affects the A 2,pq coefficients of the needle shaped particle more, especially at HH-polarization. The estimated A 2,pq coefficients decrease significantly as the needle shaped particles are oriented vertically. In addition, it is seen that the estimated E parameter is influenced by the particle orientation distribution and becomes dependent on the radar polarization. It decreases from 1.9 for uniform distribution to 1.0 in HH-polarization and 1.6 in VV-polarization for the vertically oriented needle. On the other hand, the orientation distribution hardly affects the E parameter in the case of disk shaped particles. The estimated A 2,pq coefficients slightly decrease as the disk is oriented horizontally.
Remote Sens. 2018, 10, x FOR PEER REVIEW 9 of 17 needles, oriented volume effects in wave propagation through vegetation layer lead to the differences of the , coefficients between HH-and VV-polarizations. Particularly, the , coefficient, which is associated to the attenuation of HH-polarized signal, is more affected by the orientation of the scatterer than the , coefficient. In comparing needle and disk-shaped particles, the orientation angle effect is more significant in the needle shaped particle. The , coefficient decreases significantly as the needle shaped particle is oriented vertically.
Next, the effect of particle orientation on the , coefficients is evaluated in Figure 7. As in the previous analysis, is used for the vegetation descriptor. As with the , coefficients, orientation distribution affects the , coefficients of the needle shaped particle more, especially at HH-polarization. The estimated , coefficients decrease significantly as the needle shaped particles are oriented vertically. In addition, it is seen that the estimated parameter is influenced by the particle orientation distribution and becomes dependent on the radar polarization. It decreases from 1.9 for uniform distribution to 1.0 in HH-polarization and 1.6 in VV-polarization for the vertically oriented needle. On the other hand, the orientation distribution hardly affects the parameter in the case of disk shaped particles. The estimated , coefficients slightly decrease as the disk is oriented horizontally.

Effect of Radar Observation Condition
Two experiments in the previous section discussed parameterization of the WCM in relation to the shape and orientation of scatterers. In addition to the environmental condition of the scatterer, however, radar observation parameters affect the unknown model constants in the WCM. In particular, radar scattering mechanisms in vegetated areas can be sensitive to the incidence angle and the frequency of the transmitted signal. In order to further evaluate the effect of observation conditions, we analyzed the model parameters of WCM that were estimated under different incidence angles and frequencies of the radar signal. Figure 8 shows dependence of the estimated model parameters on the incidence angle and frequency. As with the estimation procedure discussed in the previous section, the model parameters were determined by fitting B 2,pq VWC and A 2,pq m E g terms in the WCM to the simulation results of the GRG-based physical model. In this figure, the shape and orientation distribution of scatterers were kept constant, such as the moderately vertical (0 • ≤ β ≤ 60 • ) needle shaped scatterers. In addition, according to the previous analysis, the E parameter in the WCM may also depend on the radar polarization. Therefore, the E parameter for the non-linear fitting was replaced by the polarization dependent parameter E pq .

Effect of Radar Observation Condition
Two experiments in the previous section discussed parameterization of the WCM in relation to the shape and orientation of scatterers. In addition to the environmental condition of the scatterer, however, radar observation parameters affect the unknown model constants in the WCM. In particular, radar scattering mechanisms in vegetated areas can be sensitive to the incidence angle and the frequency of the transmitted signal. In order to further evaluate the effect of observation conditions, we analyzed the model parameters of WCM that were estimated under different incidence angles and frequencies of the radar signal. Figure 8 shows dependence of the estimated model parameters on the incidence angle and frequency. As with the estimation procedure discussed in the previous section, the model parameters were determined by fitting , VWC and , terms in the WCM to the simulation results of the GRG-based physical model. In this figure, the shape and orientation distribution of scatterers were kept constant, such as the moderately vertical (0°≤ ≤60°) needle shaped scatterers. In addition, according to the previous analysis, the E parameter in the WCM may also depend on the radar polarization. Therefore, the parameter for the non-linear fitting was replaced by the polarization dependent parameter .    It is seen that the B 2,VV parameter increases as an increase of the incidence angle, while the B 2,HH parameter is not affected by the incidence angle variation. The A 2,pq parameters increase as an increase of the incidence angle. The estimated A 2,pq is consistently higher for the HH-polarization case than the VV-polarization case. In addition, the difference between A 2,HH and A 2,VV is larger at a high incidence angle region. The estimated E pq parameters are not very sensitive to the incidence angle as compared to other WCM parameters. The E HH parameter is consistently higher than the E VV parameters.
On the other hand, it is observed from Figure 8d-f that the WCM parameters vary significantly accordingly with the radar frequency. The B 2,pq parameters show quadratic like behaviors as a function of the frequency. The A 2,pq increase as an increase of the frequency and there is a significant difference between A 2,HH and A 2,VV especially in high frequency. In the case of the E pq parameters, the estimated E HH and E VV generally decrease as an increase of the frequency. Particularly, the relationship between the vegetation descriptor m g and the simulated σ v,pq /2κ e,pq term becomes nearly linear for the high frequency and HH-polarization case.

Comparison with Previous Studies
The above results on the evaluation of WCM parameters suggest the optimal set of vegetation descriptors. However, the model constants, which relate vegetation descriptors to microwave scattering and attenuation terms in the vegetation layer, can vary significantly for different environmental and radar observation conditions. To understand possible variability of the WCM parameters, additional experiments on the estimation of the B 2,pq , A 2,pq , and E pq parameters were carried out at different environmental and observation conditions. In this experiment, we set up 27 configurations, as summarized in Table 2. As with the simulations in Section 3.1, the scattering and attenuation terms were calculated using vegetation parameters, given in Table 1, under different radar frequency, particle orientation distribution, and particle thickness. Here, the simulation was performed for needle shaped particles and the incidence angle was set to be 30 • . The estimation results listed in Table 2 are also illustrated in Figures 9 and 10. In these figures, the red and blue circles are the estimated WCM parameters for HH-and VV-polarizations, respectively. Figure 9 shows variations of the B 2,pq parameter for different simulation configurations. In order to examine whether these WCM parameters derived from theoretical analysis can explain actual parameters obtained from field experimental datasets, the B 2,pq parameters reported in the literature are added in these figures. Among various studies presenting B 2,pq parameters, three studies [9, 11,13], which used the same vegetation descriptor, i.e., the VWC, were selected for the comparison. In [9], the B 2,HH,[9] and B 2,VV, [9] values were determined using L-band (θ = 35 • ) and C-band (θ = 23 • ) SAR data, respectively, acquired over agricultural areas. In [11], the B 2,HH, [11] value was obtained using multi-temporal X-band (θ= [28.3 • , 32.5 • ) SAR data acquired over an experimental farm area.
The L-band B 2,HH [9] shows good agreement with the estimated B 2,HH for the moderately vertical needle cases in Figure 9. However, the C-band B 2,VV [9] and X-band B 2,HH [11] are slightly lower than the estimated B 2,pq values of this study. In [13], several B 2,HH, [13] values for different land cover types were presented using L-band (θ = 38.49 • ) scatterometer data, acquired globally. It is shown from Figure 9 that the B 2,HH, [13] values, which range from 0.01 to 0.03, are in agreement with the estimated B 2,HH of this study for the vertically oriented needle case. Figure 10 illustrates variations of the A 2,pq and E pq parameters estimated for different simulation configurations. In this case, it was not able to compare those parameters with in situ information-based parameters since we could not find study cases where the m E g was used as the V 1 vegetation descriptor. Nonetheless, it is worth noting that both of the A 2,pq parameters obtained in this study and those reported in [13] reveal the relatively larger variability than other WCM parameters. The estimated A 2,pq parameters vary significantly, both with the particle size and the orientation distribution. The E pq parameters mainly vary with the particle orientation. In particular, if there is a preference for vertically oriented scatterers in the vegetation layer, the E HH becomes close to one, which indicates a linear relationship between GRG-based model calculation and the VWC of the vegetation layer. The estimation results listed in Table 2 are also illustrated in Figures 9 and 10. In these figures, the red and blue circles are the estimated WCM parameters for HH-and VV-polarizations, respectively. Figure 9 shows variations of the , parameter for different simulation configurations. In order to examine whether these WCM parameters derived from theoretical analysis can explain actual parameters obtained from field experimental datasets, the , parameters reported in the literature are added in these figures. Among various studies presenting , parameters, three studies [9, 11,13], which used the same vegetation descriptor, i.e. the VWC, were selected for the comparison. In [9], the , ,[ ] and , ,[ ] values were determined using Lband ( = 35°) and C-band ( = 23°) SAR data, respectively, acquired over agricultural areas. In [11], the values for different land cover types were presented using L-band ( = 38.49°) scatterometer data, acquired globally. It is shown from Figure 9 that the , ,[ ] values, which range from 0.01 to 0.03, are in agreement with the estimated , of this study for the vertically oriented needle case. Figure 10 illustrates variations of the , and parameters estimated for different simulation configurations. In this case, it was not able to compare those parameters with in situ information-based parameters since we could not find study cases where the was used as the vegetation descriptor. Nonetheless, it is worth noting that both of the , parameters obtained in this study and those reported in [13] reveal the relatively larger variability than other WCM parameters. The estimated , parameters vary significantly, both with the particle size and the orientation distribution. The parameters mainly vary with the particle orientation. In particular, if there is a preference for vertically oriented scatterers in the vegetation layer, the becomes close to one, which indicates a linear relationship between GRG-based model calculation and the VWC of the vegetation layer.

Validity of WCM
Both the WCM and the theoretical model for examining the WCM parameters assume that the total backscattered signals is composed of the direct vegetation backscattering and direct ground backscattering. Consequently, the WCM predictions are problematic when the observed signal contains other scattering contributions. In order to evaluate the validity of the WCM, we consider that the total backscattered signal can be expressed by the incoherent sum of four scattering mechanisms, as with the general first-order solution of the radiative transfer equation [3][4][5], as follows: The first two terms in the right-hand side correspond to the scattering mechanisms of the WCM, such as the direct backscattering from the vegetation ( σ ) and the underlying rough surface attenuated by vegetation ( σ ). The other two terms are the scattering interaction between the vegetation and the ground surface, including the vegetation-ground and ground-vegetation double bounces (σ ) and the ground-vegetation-ground multiple bounce (σ ). The vegetation scattering and extinction terms in the vegetation layer were calculated using the first order solution of the radiate transfer equation [3][4][5]. The Integral Equation Method (IEM) [25] was used for the backscattering from the ground surface. Figure 11 illustrates an example of the variations of vegetation scattering mechanisms as a function of the incidence angle at a C-band frequency. In this simulation, the vegetation layer was assumed to be composed of uniformly distributed needle shaped scatterers ( =5 cm, =0.2 cm; =0.5 gg −1 , =3000 m −3 , ℎ= 3 m). In addition, the root-mean-square height, the correlation length, and the volumetric moisture content of the soil surface were 0.5 cm, 5 cm, and 0.3 cm 3 cm −3 , respectively. It is seen that the relative contribution of elementary scattering mechanisms to the total observed signal can vary significantly with the radar observation conditions. The ground scattering mechanism is an important contributor at the low incidence angle, whereas the relative contribution of the volume scattering component increases as the incidence angle increases. We noticed that there

Validity of WCM
Both the WCM and the theoretical model for examining the WCM parameters assume that the total backscattered signals is composed of the direct vegetation backscattering and direct ground backscattering. Consequently, the WCM predictions are problematic when the observed signal contains other scattering contributions. In order to evaluate the validity of the WCM, we consider that the total backscattered signal can be expressed by the incoherent sum of four scattering mechanisms, as with the general first-order solution of the radiative transfer equation [3][4][5], as follows: The first two terms in the right-hand side correspond to the scattering mechanisms of the WCM, such as the direct backscattering from the vegetation (σ 0 V ) and the underlying rough surface attenuated by vegetation (σ 0 G ). The other two terms are the scattering interaction between the vegetation and the ground surface, including the vegetation-ground and ground-vegetation double bounces (σ 0 VG ) and the ground-vegetation-ground multiple bounce (σ 0 GVG ). The vegetation scattering and extinction terms in the vegetation layer were calculated using the first order solution of the radiate transfer equation [3][4][5]. The Integral Equation Method (IEM) [25] was used for the backscattering from the ground surface. Figure 11 illustrates an example of the variations of vegetation scattering mechanisms as a function of the incidence angle at a C-band frequency. In this simulation, the vegetation layer was assumed to be composed of uniformly distributed needle shaped scatterers (l = 5 cm, t = 0.2 cm; m g = 0.5 gg −1 , N = 3000 m −3 , h = 3 m). In addition, the root-mean-square height, the correlation length, and the volumetric moisture content of the soil surface were 0.5 cm, 5 cm, and 0.3 cm 3 cm −3 , respectively. It is seen that the relative contribution of elementary scattering mechanisms to the total observed signal can vary significantly with the radar observation conditions. The ground scattering mechanism is an important contributor at the low incidence angle, whereas the relative contribution of the volume scattering component increases as the incidence angle increases. We noticed that there can be a significant amount of contribution from the scattering interaction between vegetation and the ground, particularly at HH-polarization and high incidence angle regions.
In addition to the observation condition, each scattering mechanism in the observed signal can be affected by the environmental condition of the scatterer, such as vegetation structures and the amount of water content. To examine the validity of the WCM, we define the relative contributions of the direct vegetation, direct ground, and the vegetation-ground interaction terms as  Figure 12 shows variations of the relative contributions as a function of the VWC under different particle orientation distributions. These simulations were carried out at a C-band for different vegetation parameters, as listed in Table 1. The contribution P D indicates the implausibility of the WCM assumption. It remains low in the case of the uniform volume layer. However, the scattering contribution of vegetation-ground interaction increases significantly with an increase of the VWC in the case of the oriented volume. In comparing HH-and VV-polarizations, the HH-polarized signal can contain more vegetation-ground interaction components, as it can be anticipated from the analysis on the B 2,HH coefficient in Section 3.2.
The relative contribution of the vegetation-ground interaction component can lead to an erroneous result of the WCM based prediction of the total backscattered signal. To further examine possible problems in the backscatter prediction of the WCM model, we define the prediction error as the difference between σ 0 Total and σ 0 Figure 13 shows the prediction errors as a function of VWC for different orientation distributions. In the case of HH-polarization, the prediction error can be up to about 6dB for the oriented volume, whereas it is less than about 1.5 dB for the random volume. The WCM based prediction error in VV-polarization can reach up to about 3 dB in the case of oriented volume. It is worth noting that, in practice, estimation of the WCM vegetation parameters and soil moisture contents are carried out based on observed signals, which contain various scattering mechanisms. Consequently, prediction errors of the WCM model for the observed signals, particularly in the case of the vegetation with preferred orientation distribution and a high VWC level, may lead to overestimation of vegetation parameters.

Conclusions
In this study, a comprehensive theoretical analysis of the vegetation parameters in the WCM is presented by examining the relationship between WCM vegetation parameters and the theoretical scattering model predictions. The advantage of implementing the WCM is in being able to express complex scattering characteristics a in vegetated area with simple bulk vegetation descriptors. However, there has been a lack of understanding or consensus about the optimal set of vegetation descriptors denoted by and variables in the WCM. By comparing the theoretical scattering

Conclusions
In this study, a comprehensive theoretical analysis of the vegetation parameters in the WCM is presented by examining the relationship between WCM vegetation parameters and the theoretical scattering model predictions. The advantage of implementing the WCM is in being able to express complex scattering characteristics a in vegetated area with simple bulk vegetation descriptors. However, there has been a lack of understanding or consensus about the optimal set of vegetation descriptors denoted by V 1 and V 2 variables in the WCM. By comparing the theoretical scattering and attenuation coefficients with several vegetation descriptors, the m g and VWC are selected as the optimal vegetation descriptors for V 1 and V 2 in the WCM, respectively.
Results show the linear relationship between the VWC and the theoretical model predictions, regardless of the scatterers' and radar observation conditions. The regression coefficient denoted by B 2,pq in the WCM varies with observation conditions, as well as vegetation structure. A vegetation canopy composed of disk shaped particles is described by higher B 2,pq than the needle shaped particles. In the case of oriented volume, like horizontal disks or vertical needles, the difference of the B 2,pq coefficients between HH-and VV-polarization becomes significant.
On the other hand, there is the power law relationship between the m g and the theoretical model predictions. Two regression parameters defining the power law relationship, denoted by the A 2,pq and E pq , are also investigated in relation to the various environmental conditions. The A 2,pq parameter varies significantly by the particle shape, size, and orientation. Higher A 2,pq values can be expected in the volume composed of randomly distributed large particles. An interesting part of the A 2,pq m E pq g term in the WCM is that the form of the power law relationship is also affected by the environmental conditions. It can be a nearly quadratic form in the case of random volume, while the exponent tends to decrease in the oriented volume.
These results provide an insight in the microwave scattering and attenuation process in the vegetation and will be helpful to predict and to interpret the SAR signal in the vegetated areas. Nonetheless, due to the significant variabilities of model parameters, it is still difficult to use the WCM for inversion schemes without in situ experimental data. The relationship between the model constants A 2,pq , E pq , and B 2,pq of the WCM and the vegetation structures suggest the possibility of narrowing down variabilities with the use of prior knowledge of the land-cover where the vegetation layer may have a preferred shape, size, and orientation distribution. Consequently, further evaluation of the inversion of the WCM model will be carried out in future research with the aid of other independent data, such as the optical and multi-frequency SAR data. Funding: This work was supported by Hanwha Systems Corporation for the core technology research project based on the image formation using flight platform radar in remote sensing. It was carried out as part of cooperation contract number U-17-020, Korea.