## 1.

## Introduction

Wheat productivity is related to food security, which plays an important role in social stability and economic development. Thus, the monitoring of growth and output in wheat areas is critical. Due to the pervasive cloud presence in some locations where wheat is planted, continuous monitoring using optical remote sensing is impossible. With the advantage of acquiring data at all-weather conditions, microwave remote sensing is helpful. Moreover, with the development of the microwave backscattering model, one can not only model the backscatter under different growing stages, but also invert the model to estimate the wheat yield.

It is well known that the backscattering coefficient is not only affected by the radar system parameters such as frequency, polarization, and incident angles, but also the surface parameters such as soil roughness and moisture, and presence and structure of vegetation. In order to monitor the surface vegetation, researchers have applied many different theories such as the vector radiative transfer (VRT) theory, analytic wave theory of random media, and the discrete scatter theory to understand the scattering mechanism from vegetated surface or wheat field in this case.

The Michigan Microwave Canopy Scattering (MIMICS) model based on the VRT theory has been developed to simulate the radar backscattering from a forested environment for the forest.^{1} The model is successful and is widely used.^{2}3.4.^{–}^{5} Toure et al.^{6} modified the MIMICS model for the wheat at the L- and C-bands in 1994. Their model explained horizontal–horizontal polarization (HH) backscattering mechanisms well at the L- and C-bands. A coherent polarimetric microwave scattering model for grassland was developed.^{7}^{,}^{8} Marliani et al.^{9} used the coherent electromagnetic model to compute the backscattering coefficient of wheat and sunflowers during the growing season, which demonstrates the potential of the interferometric observation on the crop classification algorithms based on scattering mechanisms. Champion et al.^{10} built the semiempirical model to simulate the backscattering coefficient of a wheat crop and confirmed that the radar configuration had a similar role to the canopy structure in the control of the radar backscattering coefficient. Cookmartin et al.^{11} and Brown et al.^{12} presented a comprehensive multilayer second-order radiative transfer model, which showed that the second-order terms contributed no more than 0.5 dB to the backscattering coefficient at all polarizations from wheat. Picard et al.^{13}14.15.^{–}^{16} also presented a second-order backscatter model for wheat canopies based on numerical solution of the multiple scattering Foldy–Lax equation and omitted the leaves and ears. Del Frate et al.^{17} developed an algorithm based on the radiative transfer theory to monitor the soil moisture and growth cycle of wheat.

The above models focused on the microwave backscattering characters of wheat in various aspects using a scatter mechanism and experimental data, and promoted the development of the microwave scattering theory. To a certain degree, their models simulated the components as ideal scatters and had to omit the influence of other components, thus complex formation of backscatter presents a tough task for modeling research. Closely examining the above models, one observes some deficiency in the modeling of backscattering especially the underestimation of the cross-polarized term. One possible cause is the exclusion of wheat ears in the modeling. Therefore, a microwave backscattering model of winter wheat-based on the VRT theory is presented. A wheat ear is one of the scattering components in modeling. Descriptions of the model development and verification as well as model prediction of wheat yield are detailed.

## 2.

## Model

In modeling, winter wheat is divided into two types of scattering components (wheat stem and wheat leaves) before the heading period. After the heading, model component for wheat ears is added as the third component (Fig. 1). The soil under the wheat crop is assumed to be a random rough surface. The stem is modeled as a vertical finite length cylinder and the leaf is modeled as an ellipsoid with a limited leaf inclination angle. The wheat stem and leaf are distributed evenly in the vertical range from $z=-d$ to $z=0$ before the heading period. The wheat ears modeled as elliptical scatterers within a layer between $z=0$ and $z=-h$ are inserted ($d$: aggregate thickness; $h$: ears thickness) (Fig. 1). M0 could be taken as the soil contribution; M1 represents the individual contribution from wheat leaves, while M2 indicates stems; M3 indicates wheat ears’ response, and M4 can be explained as mixed contributions from wheat parameters and other corresponding parameters. As demonstrated in the previous studies as well as this study, the one-layer scattering model is suitable to simulate the backscatter from wheat, especially before the heading stage. The addition of the second layer to the model should accommodate for the scattering contribution from wheat ears after the heading. The new vertical variable for the extinction matrix and phase matrix should account for full polarimetric scattering of the nonuniform distributed and elliptical scatter,^{18} which is the cause for the improvement in the simulation of cross-polarized backscattering.

The VRT is given in Eq. (1):^{18}

## (1)

$$\frac{\mathrm{d}\mathsf{I}(\mathsf{r},\widehat{s})}{\mathrm{d}s}=-{\overline{\mathsf{k}}}_{e}(\widehat{s})\xb7\mathsf{I}(\mathsf{r},\widehat{s})+\int \mathrm{d}{s}^{\prime}\overline{\mathsf{P}}(s,{s}^{\prime})\xb7\mathsf{I}(r,{s}^{\prime}).$$The VRT equation describes the scattering, absorbed, multiple scattering, and transferred processes of electromagnetic wave intensities [Stokes intensities $\mathsf{I}(\mathsf{r},\widehat{S})$] along the direction. The extinction matrix ${\overline{\mathsf{k}}}_{e}$ describes the attenuation of $\mathsf{I}(\mathsf{r},\widehat{S})$ caused by scatter and absorption. The phase matrix describes the multiple scattering progresses from direction $({\theta}^{\prime},{\varphi}^{\prime})$ to direction $(\theta ,\varphi )$.

By means of the forward-scattering theorem, the extinction matrix can be written as^{18}

## (2)

$${\overline{\mathsf{k}}}_{e}(\theta ,\phi )=\frac{4\pi}{k}{n}_{0}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{Im}[\widehat{p}\xb7\overline{\mathsf{S}}(\theta ,\varphi ;\theta ,\varphi )\xb7\widehat{p}]\phantom{\rule{0ex}{0ex}}=\frac{2\pi}{k}{n}_{0}\left[\begin{array}{cccc}2\text{\hspace{0.17em}}\mathrm{Im}\langle {S}_{vv}^{0}\rangle & 0& \mathrm{Im}\langle {S}_{vh}^{0}\rangle & -\mathrm{Re}\langle {S}_{vh}^{0}\rangle \\ 0& 2\text{\hspace{0.17em}}\mathrm{Im}\langle {S}_{hh}^{0}\rangle & \mathrm{Im}\langle {S}_{hv}^{0}\rangle & \mathrm{Re}\langle {S}_{hv}^{0}\rangle \\ 2\text{\hspace{0.17em}}\mathrm{Im}\langle {S}_{hv}^{0}\rangle & 2\text{\hspace{0.17em}}\mathrm{Im}\langle {S}_{vh}^{0}\rangle & \mathrm{Im}\langle {S}_{vv}^{0}+{S}_{hh}^{0}\rangle & \mathrm{Re}\langle {S}_{vv}^{0}-{S}_{hh}^{0}\rangle \\ 2\text{\hspace{0.17em}}\mathrm{Re}\langle {S}_{hv}^{0}\rangle & -2\text{\hspace{0.17em}}\mathrm{Re}\langle {S}_{vh}^{0}\rangle & \mathrm{Re}\langle {S}_{hh}^{0}-{S}_{vv}^{0}\rangle & \mathrm{Im}\langle {S}_{vv}^{0}+{S}_{hh}^{0}\rangle \end{array}\right],$$The phase matrix is a transfer matrix which describes the scattering energy transfer from direction $({\theta}^{\prime},{\varphi}^{\prime})$ to direction $(\theta ,\varphi )$, which has a similar definition as that of the Mueller matrix. The matrix is given by^{18}

## (3)

$$\overline{P}(\theta ,\varphi ;{\theta}^{\prime},{\varphi}^{\prime})\phantom{\rule{0ex}{0ex}}={n}_{0}\left[\begin{array}{cccc}\langle {|{S}_{vv}|}^{2}\rangle & \langle {|{S}_{vh}|}^{2}\rangle & \mathrm{Re}\langle {S}_{vv}{S}_{vh}^{*}\rangle & -\mathrm{Im}\langle {S}_{vv}{S}_{vh}^{*}\rangle \\ \langle {|{S}_{hv}|}^{2}\rangle & \langle {|{S}_{hh}|}^{2}\rangle & \mathrm{Re}\langle {S}_{hv}{S}_{hh}^{*}\rangle & -\mathrm{Im}\langle {S}_{hv}{S}_{hh}^{*}\rangle \\ 2\text{\hspace{0.17em}}\mathrm{Re}\langle {S}_{vv}{S}_{hv}^{*}\rangle & 2\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{Re}\langle {S}_{vh}{S}_{hh}^{*}\rangle & \mathrm{Re}\langle {S}_{vh}{S}_{hh}^{*}+{S}_{vh}{S}_{hv}^{*}\rangle & -\mathrm{Im}\langle {S}_{vv}{S}_{hh}^{*}-{S}_{vh}{S}_{hv}^{*}\rangle \\ 2\text{\hspace{0.17em}}\mathrm{Im}\langle {S}_{vv}{S}_{hv}^{*}\rangle & 2\text{\hspace{0.17em}}\mathrm{Im}\langle {S}_{vh}{S}_{hh}^{*}\rangle & \mathrm{Im}\langle {S}_{vv}{S}_{hh}^{*}+{S}_{vh}{S}_{hv}^{*}\rangle & \mathrm{Re}\langle {S}_{vv}{S}_{hh}^{*}-{S}_{vh}{S}_{hv}^{*}\rangle \end{array}\right].$$According to the VRT theory, the extinction matrix mentioned in Eq. (4) is given by Eq. (2). Thus, we can get the VRT equations of this model^{18}^{,}^{19}

## (4)

$$\theta \frac{d}{dz}{I}_{r}^{+}(\theta ,\varphi ,z)=-{K}_{r}^{+}(\theta ,\varphi )\xb7{I}_{r}^{+}(\theta ,\varphi ,z)+{F}_{r}^{+}(\theta ,\varphi ,z),$$## (5)

$$-\theta \frac{d}{dz}{I}_{r}^{-}(-\theta ,\varphi ,z)=-{K}_{r}^{-}(-\theta ,\varphi )\xb7{I}_{r}^{-}(-\theta ,\varphi ,z)+{F}_{r}^{-}(-\theta ,\varphi ,z).$$In Eqs. (4) and (5), ${I}_{r}^{+}(\theta ,\varphi ,z)$ and ${I}_{r}^{-}(-\theta ,\varphi ,z)$ are the upward and downward Stokes intensities, ${K}_{r}^{\pm}$ is the extinction matrix, and ${F}_{r}^{+}(\theta ,\varphi ,z)$ and ${F}_{r}^{-}(-\theta ,\varphi ,z)$ are the scattering source functions which represent the scattering contribution by the wheat layer. The phase matrix mentioned in Eq. (6) is given by Eq. (3).

These scattering source functions are^{20}^{,}^{21}

## (6)

$${F}_{r}^{+}(\theta ,\varphi ,z)=\frac{1}{\theta}[{\int}_{0}^{2\pi}{\int}_{0}^{1}{P}_{r}(\theta ,\varphi ;{\theta}^{\prime},{\varphi}^{\prime}){I}_{r}^{+}({\theta}^{\prime},{\varphi}^{\prime},z)\mathrm{d}{\theta}^{\prime}\mathrm{d}{\varphi}^{\prime}\phantom{\rule{0ex}{0ex}}+{\int}_{0}^{2\pi}{\int}_{0}^{1}{P}_{r}(\theta ,\varphi ;-{\theta}^{\prime},{\varphi}^{\prime}){I}_{r}^{-}(-{\theta}^{\prime},{\varphi}^{\prime},z)\mathrm{d}{\theta}^{\prime}\mathrm{d}{\varphi}^{\prime}],$$## (7)

$${F}_{r}^{-}(-\theta ,\varphi ,z)=\frac{1}{\theta}[{\int}_{0}^{2\pi}{\int}_{0}^{1}{P}_{r}(-\theta ,\varphi ;{\theta}^{\prime},{\varphi}^{\prime}){I}_{r}^{+}({\theta}^{\prime},{\varphi}^{\prime},z)\mathrm{d}{\theta}^{\prime}\mathrm{d}{\varphi}^{\prime}\phantom{\rule{0ex}{0ex}}+{\int}_{0}^{2\pi}{\int}_{0}^{1}{P}_{r}(-\theta ,\varphi ;-{\theta}^{\prime},{\varphi}^{\prime}){I}_{r}^{-}(-{\theta}^{\prime},{\varphi}^{\prime},z)\mathrm{d}{\theta}^{\prime}\mathrm{d}{\varphi}^{\prime}],$$In order to solve Eqs. (4) and (5), one needs to have the boundary conditions given as^{20}^{,}^{21}

## (8)

$${I}_{r}^{-}(-\theta ,\varphi ,z=0)={I}_{0}\delta (\theta -{\theta}_{0})\xb7\delta (\varphi -{\varphi}_{0}),$$Solving Eqs. (4)–(9) iteratively as a group, one can obtain the first-order solution^{20}

## (10)

$${I}^{bs}({\theta}_{0},{\varphi}_{0})={I}_{1}^{+}({\theta}_{0},{\varphi}_{0}\pm \pi ,z=0)=T({\theta}_{0},{\varphi}_{0}){I}_{0},$$## (11)

$$T({\theta}_{0},{\varphi}_{0})=\frac{1}{{\theta}_{0}}Q{D}_{r}({\theta}_{0},{\varphi}_{0}\pm \pi ;-h/{\theta}_{0})\xb7{Q}^{-1}{R}_{w}({\theta}_{0},{\varphi}_{0}\pm \pi )Q{A}_{1}{Q}^{-1}\phantom{\rule{0ex}{0ex}}\xb7{R}_{w}({\theta}_{0},{\varphi}_{0}){D}_{r}({\theta}_{0},{\varphi}_{0};-h/{\theta}_{0}){Q}^{-1}\phantom{\rule{0ex}{0ex}}+\frac{1}{{\theta}_{0}}Q{D}_{r}({\theta}_{0},{\varphi}_{0}\pm \pi ;-h/{\theta}_{0})\xb7{Q}^{-1}{R}_{w}({\theta}_{0},{\varphi}_{0}\pm \pi )Q{A}_{2}{Q}^{-1}\phantom{\rule{0ex}{0ex}}+\frac{1}{{\theta}_{0}}Q{A}_{3}{Q}^{-1}{R}_{w}({\theta}_{0},{\varphi}_{0})\xb7Q{D}_{r}({\theta}_{0},{\varphi}_{0};-h/{\theta}_{0})\xb7{Q}^{-1}\phantom{\rule{0ex}{0ex}}+\frac{1}{{\theta}_{0}}Q{A}_{4}{Q}^{-1}+Q{A}_{5}{Q}^{-1},$$## (12)

$${D}_{r}{(\theta ,\varphi ;-z/\theta )}_{(i,i)}=\mathrm{exp}[-{\lambda}_{i}(\theta ,\varphi )z/\theta ],$$${A}_{1}$ to ${A}_{5}$ (Ref. 20) in Eq. (11) are given as

## (13)

$${A}_{1}={\int}_{-h}^{0}{D}_{r}(-\theta ,\varphi ;-({z}^{\prime}+h)/\theta ){Q}^{-1}{P}_{r}(-\theta ,\varphi ;{\theta}_{0},{\varphi}_{0})Q{D}_{r}({\theta}_{0},{\varphi}_{0};-({z}^{\prime}+h)/{\theta}_{0})\mathrm{d}{z}^{\prime}$$## (14)

$${A}_{2}={\int}_{-h}^{0}{D}_{r}(-\theta ,\varphi ;-({z}^{\prime}+h)/\theta ){Q}^{-1}{P}_{r}(-\theta ,\varphi ;{\theta}_{0},{\varphi}_{0})Q{D}_{r}(-{\theta}_{0},{\varphi}_{0};{z}^{\prime}/{\theta}_{0})\mathrm{d}{z}^{\prime}$$## (15)

$${A}_{3}={\int}_{-h}^{0}{D}_{r}(\theta ,\varphi ;{z}^{\prime}/\theta ){Q}^{-1}{P}_{r}(\theta ,\varphi ;{\theta}_{0},{\varphi}_{0})Q{D}_{r}({\theta}_{0},{\varphi}_{0};-({z}^{\prime}+h){\theta}_{0})\mathrm{d}{z}^{\prime}$$## (16)

$${A}_{4}={\int}_{-h}^{0}{D}_{r}(\theta ,\varphi ;{z}^{\prime}/\theta ){Q}^{-1}{P}_{r}(\theta ,\varphi ;-{\theta}_{0},{\varphi}_{0})Q{D}_{r}(-{\theta}_{0},{\varphi}_{0};{z}^{\prime}/{\theta}_{0})\mathrm{d}{z}^{\prime}$$## (17)

$${A}_{5}={D}_{r}(\theta ,\varphi ;-h/\theta ){Q}^{-1}{R}_{w}(\theta ,\varphi ;-{\theta}_{0},{\varphi}_{0})Q{D}_{r}(-{\theta}_{0},{\varphi}_{0};-h/{\theta}_{0})\mathrm{.}$$Therefore, the backscattering coefficient of wheat is

## 3.

## Experiments

The study area is at the crop research and demonstration site of Qianjin Town, Qionglai County, Sichuan Province, China. The location is 30°24′22.29″N and 103°32′15.97″E. The topography is gentle with an elevation at 487 m above mean sea level. The site is jointly managed and owned by the University of Electronic Science and Technology of China and Sichuan Agriculture University. Winter wheat is one of the research and demonstration crops, and its growth period is about 200 days from November to May. The studied wheat area is 3 ha.

A scatterometer (Table 1) was used to collect the backscattering coefficient during the period of wheat growth. Measurements at 11 time intervals (Table 2) were completed. At each time interval, full polarized C-band data were collected with elevation angles ranging from 10 deg to 70 deg with an increment of 3 deg. The C-band is chosen because its wavelength is similar to that of the wheat and it is sensitive to wheat parameters such as leaves. The measured backscattering coefficients’ value is the average of 6 to 9 measurements (in 6 to 9 different measurements, the azimuths are different at the same incident angle). It should be noted that the booting and filling stages are the most critical growth phases to determine the wheat productivity. The booting stage is the final stage of vegetative growth, which lays a solid foundation for the reproductive growth if the growth is flourishing [Fig. 2(a)]. In the filling stage [Fig. 2(b)], the wheat stores starch, protein, and other organic materials produced by photosynthesis in the grain, which is the crucial moment for the harvest.

## Table 1

System parameters of C-band ground-based scatterometer.

Center frequency | 5.3 GHz |
---|---|

Polarization | HH/VV/VH/HV |

Bandwidth | $\ge 0.6\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{GH}$ |

Work system | frequency-modulated continuous-wave |

Operation slant range | $\ge 100\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}$ |

Transmit power | $\ge 50\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{mw}$ |

Dynamic range | $\ge 45\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{dB}$ |

Measurement accuracy | 0.5 dB |

Range of elevation angle | 0 deg to 90 deg |

Range of azimuth angle | 0 deg to 359 deg |

Receiving antenna beam width | 6 deg |

Transmitting antenna width | 8 deg |

## Table 2

Measurement schedule at 11 time intervals.

Growth phase | Date | Growing days since seeding |
---|---|---|

Infantile state | 13 Nov 2010 | 15 |

Two leaves | 3 Dec 2010 | 35 |

The early stage of tillering | 25 Dec 2010 | 57 |

The late stage of tillering | 14 Jan 2011 | 77 |

Jointing stage | 25 Feb 2011 | 119 |

Flag leaf stage | 4 Mar 2011 | 126 |

Booting stage | 19 Mar 2011 | 141 |

Heading stage | 1 Apr 2011 | 153 |

Flowering stage | 16 Apr 2011 | 168 |

Filling stage | 29 Apr 2011 | 181 |

Mature stage | 18 May 2011 | 200 |

After the processing of the initial field data, values at the booting and filling stages were obtained as shown in Table 3. All data were mean values. The stalk moisture and ear moisture were derived from the fresh weight and dry weight of the stalk and ears, respectively. The correlation length and roughness of the ground surface were obtained by digitizing the surface (vertical) profile photographed against white boards at the sample site.

## Table 3

The parameter list of the wheat and ground surface data at the booting and filling stages.

Model parameters | At booting stage | At filling stage |
---|---|---|

Leaf thickness (m) | 0.0002 | 0.00026 |

Leaf minor axis (m) | 0.0175 | 0.0168 |

Leaf major axis (m) | 0.229 | 0.206 |

Leaf moisture (%) | 88.3 | 62.3 |

Leaf density ($/{\mathrm{m}}^{3}$) | 1490 | 1392 |

Stalk length (m) | 0.455 | 0.67 |

Stalk radius (m) | 0.0028 | 0.0026 |

Stalk moisture (%) | 89.7 | 74.2 |

Stalk density ($/{\mathrm{m}}^{3}$) | 298 | 348 |

Height of wheat crops (m) | 0.621 | 0.86 |

Ear length (m) | 0.14 | |

Ear radius (m) | 0.0052 | |

Ear moisture (%) | 63.1 | |

Ear density ($/{\mathrm{m}}^{3}$) | 2137 | |

Surface correlation length (m) | 0.08072 | 0.08072 |

Surface roughness (m) | 0.00636 | 0.00636 |

Soil moisture (%) | 27.4 | 22.9 |

The parameterization of the structural parameters of wheat for extinction and phase matrices was on the basis of Jin and Xu.^{20}

## 4.

## Results and Discussions

Even though scatterometer and field measurements were conducted at 11 time intervals, modeling at the booting and filling stages was given next because of the stated importance in wheat productivity. The growth and output situations are highly related to the booting stage and filling stage. Therefore, the simulation results in these two stages are compared with the measured values and simulation results of the modified MIMICS model. When putting MIMICS into the wheat scattering characteristic research, the common methods simplify its tall tree structure. We take the wheat stem as the main branch, wheat leaves as tree leaves, ears as secondary branches. The trunk has to be omitted since it has a large body.

## 4.1.

### Comparison Between Simulation Values Modeled with Measured Values

Observed and measured C-band HH, vertical–vertical polarization (VV), and horizontal–vertical polarization (HV) backscatter coefficients versus incidence angles at the booting and filling stages were shown in Fig. 3. Figures 3(a), 3(c), and 3(e) were at the booting stage, and Figs. 3(b), 3(d), and 3(f) were at the filling stage. Results in HH-, VV-, and HV-polarizations were the first, second, and third rows, respectively. Figure 3(a) is the HH-polarization, while the second row is the VV-polarization, and the third row is the VH-polarization. As shown in Fig. 3, both the observed and simulated backscatter coefficients vary with the change of incidence angles. At each incidence angle and each polarization, there might be a general agreement between the observed and modeled values since the maximum difference was less than 3 dB overall (Fig. 3). Multipolarized simulation results and measured values at the C-band vary with incidence angles. Because the HV-polarization curve is similar to that of the VH-polarization, it has not been shown here. The maximum difference between the measured values and simulation results is no more than 3 dB. Thus, it can be concluded that the simulation results from the wheat microwave backscattering model agree well with the values measured in a real-world situation. In this experiment, numerous variables affect the measurement accuracy, thus the range of 3 dB error could be taken as the normal error in the microwave scattering measurement. Around the previously mentioned two stages, since backscatter is sensitive to the flag leaves, the flag leaves might cause the oscillation of backscattering with the variation of incidence angles. The backscattering coefficients decreased and produced the wave valley in these curves. The previous researchers show that the backscattering coefficient decreases with an increase in the incidence angles. Otherwise, the curve was different from the actual situation, the reason of which could be the simplification of the extinction matrix and phase matrix. Figure 3 shows the simulation and the measured values. The main reason for this fit that is the measured data selected from the measured data base can fits the simulation well. Meanwhile, the measured data had quite a bit of difference compared to the simulation results from modified MIMICS. We consider it suitable to compare the measured data to the simulation results by our model.

## 4.2.

### Backscattering Coefficients Comparison among Values by Model, Values from MIMICS and Values Measured by Scatterometer

The MIMICS model was developed for radar backscattering in a forested environment. The usefulness of this model for vegetation such as wheat was explored. When the model is applied to wheat, there is no obvious stratification between the stalk and canopy. The tree trunk does not exist. In this study, a wheat stem was treated as a main branch in MIMICS model, a wheat leaf as a “tree” leaf, and a wheat ear as a secondary branch. The tree trunk is ignored. With these modifications, modeled results were obtained. As compared to the observed values, there was an overall agreement in the copolarized backscatter coefficients. However, modeled HV backscattering coefficients were significantly lower than those observed (Fig. 4). A possible cause was attributed to the cross-polarized backscattering simulated in modified MIMICS and underestimated in the modification. The simulation value of the model built in the paper is flat because the simplification of the extinction matrix and phase matrix make the model less effectiveness than the modified MIMICS model. The cross-polarized backscattering is primarily derived from multiscattering phenomena. However, when the MIMICS model is applied to simulate the wheat scattering value, several M4 values contributing to cross-polarized backscattering are omitted. The removal of trunks leads to no second-order scattering item for the trunk–ground interactions. The second-order term is another important contributor to the cross-polarized backscattering. Second, the secondary branch and the leaf have the same ranges in MIMICS, but the ears and the leaves of wheat have different ranges. If the wheat ear is taken as the secondary branch, errors result.

## 4.3.

### HV Backscattering Coefficients After Booting Stage

After the booting stage or 141 days since seeding in this study, the model component for the wheat ear was added to the model. The assessment was carried out with the incidence angles ranging from 20 deg to 70 deg. The component contributed greatly to cross-polarized backscatter, whereas the component of the copolarized scattering could be limited. The cross-polarized backscatter was elevated by the inclusion of this component (Fig. 5). At an incident angle of 30 deg, the elevated value increased as the growing season continued after booting. The different values between the simulation values without ears and those actually measured reached a maximum at 200 growth days, the mature stage in the wheat growth cycle.

## 4.4.

### Discussion

The model proposed in this paper simulated the backscatter at two stages, the booting and filling stages, respectively. Focusing on the booting stage, the research studied the microwave characteristics of wheat and prepared a comparison of the backscattering values between wheat with and without ears. Because protein, starch, and organic matter are assimilated in the grain by photosynthesis in the filling stage, the research chose the filling stage as the study stage. The research aim is to propose a model that can better simulate the backscattering coefficients for wheat with ears. The stages with and without ears could be compared in the model. We could find the result without ears had quite a difference from the measured ones, while the results with ears agreed well with the measured results (Fig. 5), so the simulation results of the model showed the ears play a key role in the backscatter contribution after booting stage, and could verify the assumption we proposed the exclusion of wheat ears in modeling could cause the underestimation of cross polarization. Taking wheat ears as a necessary modeling component is acceptable. However, the model has some disadvantages compared with the modified MIMICS model. The copolarization simulation results were not applicable and had worse effects than the modified MIMICS model of wheat. The model should be modified to fit the full-polarization data and promoted to fit normal situations such as moisture and roughness of soil, topography, and wheat parameters. The next step for developing the model is adding the microwave characteristics of wheat with ears to the yields’ impact factors (soil moisture, temperature, soil moisture, variety and growth parameters wheat, and so on) and build the wheat yields’ estimation model.

## 5.

## Conclusion

A backscattering model for wheat has been developed using the VRT theory. The layer of wheat crop is comprised of two kinds of scattering particles before the heading period, the wheat stem and the leaf of wheat, and three kinds of scattering particles after the heading period, the wheat stem, the wheat leaf, and the wheat ear. To account for the influence on wheat ear, two submodels were designed. Before the booting stage, there was no wheat ear. The first submodel consisted of wheat stem and leaf layer over ground surface. After booting, the model component for the wheat ear was added. Thus, the second submodel included the stem, leaf, and ear over the ground surface.

Compared to the measured values, the accuracy of the proposed model prove that the proposed model has a better cross-polarization than MIMICS by comparing the simulation results of the proposed model to the measured values and simulation results of the MIMICS model.

On the basis of this model, the paper analyzed the cross-polarized backscattering coefficients and found that the wheat ears have considerable impact on the cross-polarization. The results may be used to interpret remote sensing images, monitor the crop growth, and make predictions of the crop production.

## Acknowledgments

This work is supported by the National Natural Science Foundation of China (Nos. 41371340 and 41071222). The Prior Research Program of the 12th Five-year Civil Aerospace Plan (D040201-04).

## References

## Biography

**Bo Huang** is studying in the School of Automation Engineering at the University of Electronic Science and Technology of China. His research fields include microwave measurement and modeling, and crop growth monitoring.