Soil Moisture Retrieval Using SAR Backscattering Ratio Method during the Crop Growing Season

: Soil moisture content (SMC) is an indispensable basic element for crop growth and development in agricultural production. Obtaining accurate information on SMC in real time over large agricultural areas has important guiding signiﬁcance for crop yield estimation and production management. In this study, the paper reports on the retrieval of SMC from RADARSAT-2 polarimetric SAR data. The proposed SMC retrieval algorithm includes vegetation correction based on a ratio method and roughness correction based on the optimal roughness method. Three vegetation description parameters (i.e., RVI, LAI, and NDVI) serve as vegetation descriptors in the parameterization of the algorithm. To testify the vegetation correction result of the algorithm, the water cloud model (WCM) was compared with the algorithm. The calibrated integrated equation model (CIEM) was employed to describe the backscattering from the underlying soil. To improve the accuracy of SMC retrieval, the CIEM model was optimized by using the optimal roughness parameter and the normalization method of reference incidence angle. Validation against ground measurements showed a high correlation between the measured and estimated SMC when the NDVI serves as vegetation descriptor (R 2 = 0.68, RMSE = 4.15 vol.%, p < 0.01). The overall estimation performance of the proposed SMC retrieval algorithm is better than that of the WCM. It demonstrates that the proposed algorithm has an operational potential to estimate SMC over wheat ﬁelds during the growing season.


Introduction
Soil moisture content (SMC) is an important parameter in the change process of land surface systems, playing a critical role in agricultural management, surface water cycling, and energy exchange in near-earth space [1][2][3][4][5]. SMC is an indispensable basic element for crop growth and development in agricultural production. The SMC can be used as an effective index of farmland health. The content of soil moisture is directly related to the health status of crops. Drought can lead to soil drying and cracking, hinder crop growth, reduce crop yield, and threaten food security. Therefore, access to accurate SMC information on farmland is of great significance for crop yield estimation, drought monitoring, and agricultural production management. Traditionally, SMC monitoring and investigation methods were ground-based point measurements or interpolated lines or grids. The implementation is not only time-consuming and labor-intensive, but also difficult to obtain spatially continuous soil moisture information in time, which has great limitations for its operational use [6]. Satellite remote sensing has the advantages of being fast, macro, and dynamic, providing an effective means to obtain soil moisture information.
For visible-infrared bands, the SMC estimation is from the relationships between soil spectral reflectance and soil moisture or vegetation traits occur under water stress, such as the single spectral analysis method [7][8][9] and the vegetation index method [10][11][12]. For thermal-infrared bands, the SMC is derived from the parameters related to the soil thermal properties, such as the thermal inertia method [13][14][15] and the temperature index method [16,17]. However, the optical/thermal methods are associated with some inherent limitations, such as vegetation interference, night effects, poor temporal resolution, and cloudy conditions [18]. Comparing with optical/thermal remote sensing, Synthetic Aperture Radar (SAR) has all-weather day-and-night penetration capabilities, and can provide effective and cost-efficient means of monitoring and assessing SMC over vegetation cover conditions [19][20][21]. SAR remote sensing retrieves SMC by establishing relationships between the radar backscattering coefficient and soil moisture, because the amount of water held in the soil is directly proportional to the dielectric constant of soil [22]. The basis of SMC retrieval is the soil dielectric constant. However, besides the soil dielectric constant, SAR signal is also affected by other factors, such as surface roughness, vegetation cover, and radar system parameters. Therefore, if we only pay attention to the relationship between soil moisture and the SAR backscattering coefficient and ignore the influence of other parameters, we can only obtain a trend of overall change, which is far from meeting the requirements of SMC retrieval. There are still some challenges in retrieving SMC using SAR data over vegetated areas. The observed SAR signals include the combined information of vegetation moisture and SMC because the vegetation canopy contains water [23,24]. The vegetation canopy volume scattering, the interaction between soil and vegetation and the observed backscattering are highly nonlinear [25][26][27], which increases the difficulty of SMC retrieval. When using SAR to retrieve SMC under vegetation cover, the SMC is underestimated if vegetation effect is not considered [28,29]. Therefore, removing the influence of vegetation is the key in SMC retrieval using SAR.
To remove the vegetation effects, a number of methods has been established by scholars. These methods can be generally classified into two categories: methods based on polarization decomposition and methods based on vegetation scattering models [30]. Hajnsek et al. [31] used surface, dihedral, and vegetation scattering to interpret the scattering processes and to retrieve SMC. The potential of using polarimetric decompose for estimating volumetric SMC under agricultural vegetation was proved. Combining polarimetric decomposition techniques with the concept of multi-angularity, Jagdhuber et al. proposed a retrieval technique for SMC under low vegetation cover [32]. He et al. estimated SMC in agricultural crop fields from fully polarimetric SAR data through the polarimetric decomposition of a SAR coherency matrix. Martino et al. [33] used a polarimetric two-scale two-component model (PTSTCM) to describe the surface scattering component, and a randomly (uniformly, vertically, or horizontally) oriented thin-dipole model to describe the volume scattering contribution from the vegetation layer that covers the scattering surface. The problem of negative power in SMC retrieval is significantly mitigated. It is noteworthy that polarization decomposition method for soil moisture retrieval requires fully polarization data. Vegetation scattering models describing the scattering processes of microwave signals have been established (such as Michigan microwave canopy scattering (MIMICS) [34]; Water cloud model (WCM) [35]; the distorted Born approximation (DBA) [36]; the Tor Vergata model (TVM) [37]), which can also be used to remove the influence of vegetation cover on soil moisture retrieval. Among these models, WCM has been widely used for SMC retrieval over vegetated areas because it only requires a few parameters and has high accuracy [30].
SAR backscattering includes mixed signals of vegetation and soil from agricultural fields, while vegetation contributions vary at different growth stages [38]. In addition, soil surface roughness also affects radar backscattering [39]. Over vegetated areas, soil surface roughness cannot be easily measured directly, which has caused some problems for soil moisture retrieval. Therefore, the key issue in the quantitative estimation of SMC is to separate the effect of vegetation and soil roughness. In this paper, a method was developed to retrieve SMC over agriculture fields during the crop growing season. The proposed method combined the optimal roughness method and ratio method. The ratio method was used to separate the scattering contribution of vegetation from total backscattering. The vegetation correction procedure is embedded within a SMC retrieval algorithm, in which the surface scattering component is formulated by the Calibrated Integrated Equation Model (CIEM) [40]. For the problem that roughness cannot be measured, the optimal roughness was defined by a global search method. Finally, the developed method is applied to estimate SMC over agriculture fields. In addition, the proposed model was compared with the WCM to evaluate its effectiveness and accuracy.

Ground Measurements
The study area is a relatively flat agricultural region located at London (Longitude: 81 • 5 W, Latitude: 42 • 7 N), southwestern Ontario, Canada ( Figure 1). The selected area is a rain-fed agricultural area; therefore, SMC will not be affected by irrigation.
Remote Sens. 2022, 14, x FOR PEER REVIEW 3 of 22 moisture retrieval. Therefore, the key issue in the quantitative estimation of SMC is to separate the effect of vegetation and soil roughness. In this paper, a method was developed to retrieve SMC over agriculture fields during the crop growing season. The proposed method combined the optimal roughness method and ratio method. The ratio method was used to separate the scattering contribution of vegetation from total backscattering. The vegetation correction procedure is embedded within a SMC retrieval algorithm, in which the surface scattering component is formulated by the Calibrated Integrated Equation Model (CIEM) [40]. For the problem that roughness cannot be measured, the optimal roughness was defined by a global search method. Finally, the developed method is applied to estimate SMC over agriculture fields. In addition, the proposed model was compared with the WCM to evaluate its effectiveness and accuracy.

Ground Measurements
The study area is a relatively flat agricultural region located at London (Longitude: 81°5′W, Latitude: 42°7′N), southwestern Ontario, Canada ( Figure 1). The selected area is a rain-fed agricultural area; therefore, SMC will not be affected by irrigation. Field campaigns were conducted in the study area from May to July in 2019 simultaneously with each RADARSAT-2 SAR data acquisition. A total of eight field visits was carried out, and 240 sample data were collected throughout the wheat growing season. Selected photos taken during in the field campaigns are shown in Figure 2. The spatial distance between any two adjacent sampling points is no less than 30 m. A time domain reflectometry (TDR) instrument was used to measure the SMC (vol.%) at 0-5 cm of the soil surface. Each sampling point was measured six times, and the average value of the six measurements was taken as the soil moisture reading of the sampling location. Realtime kinematic (RTK) carrier phase division technology was used to record the longitude and latitude information of each sampling point, so that the image parameter information corresponding to the sampling point can be extracted from the remote sensing image. The Field campaigns were conducted in the study area from May to July in 2019 simultaneously with each RADARSAT-2 SAR data acquisition. A total of eight field visits was carried out, and 240 sample data were collected throughout the wheat growing season. Selected photos taken during in the field campaigns are shown in Figure 2. The spatial distance between any two adjacent sampling points is no less than 30 m. A time domain reflectometry (TDR) instrument was used to measure the SMC (vol.%) at 0-5 cm of the soil surface. Each sampling point was measured six times, and the average value of the six measurements was taken as the soil moisture reading of the sampling location. Real-time kinematic (RTK) carrier phase division technology was used to record the longitude and latitude information of each sampling point, so that the image parameter information corresponding to the sampling point can be extracted from the remote sensing image. The leaf area index (LAI) was extracted by using an optical fisheye camera (Nikon d300s) with a focal length of 10.5 mm and can-eye software. Unfortunately, due to instrument failure, the LAI could not be obtained on June 16. In addition, the normal difference vegetation index (NDVI) was extracted from Sentinel-2. However, in the two days before and after the SAR observation, only six cloudless optical images were available (Table 1)

RADARSAT-2 Data
During the wheat growing season (May-July 2019), eight RADARSAT-2 SAR acquired over the study area were provided by the Canadian Space Agency thro SOAR (Science and Operational Applications Research) program. The nominal spa olution of RADARSAT-2 images was approximately 8 m. The satellite revisit peri 24 days, and the orbit period was 100.7 min. For details of the acquired RADA SAR images, please see reference [20]. The RADARSAT-2 SAR preprocessing ma cludes radiometric calibration and speckle filtering.
The local incidence angle of SAR usually has an effect on the surface backsca coefficient, which should not be ignored in SMC estimation. In order to limit the in of different local incident angles, the theoretical method proposed by Ulaby et which is based on Lambert's law for optics [42], was used to normalize the backsca coefficients in this study:

RADARSAT-2 Data
During the wheat growing season (May-July 2019), eight RADARSAT-2 SAR images acquired over the study area were provided by the Canadian Space Agency through the SOAR (Science and Operational Applications Research) program. The nominal spatial resolution of RADARSAT-2 images was approximately 8 m. The satellite revisit period was 24 days, and the orbit period was 100.7 min. For details of the acquired RADARSAT-2 SAR images, please see reference [20]. The RADARSAT-2 SAR preprocessing mainly includes radiometric calibration and speckle filtering.
The local incidence angle of SAR usually has an effect on the surface backscattering coefficient, which should not be ignored in SMC estimation. In order to limit the influence of different local incident angles, the theoretical method proposed by Ulaby et al. [41], which is based on Lambert's law for optics [42], was used to normalize the backscattering coefficients in this study: where θ is the local incidence angle, and θ re f is the reference incidence angle.

Methods
In this study, the scattering model CIEM was used to simulate the backscattering coefficients of soil surface covered by winter wheat. After the vegetation correction, the soil moisture estimation and accuracy verification were carried out using a look-up table approach. The workflow of SMC retrieval is shown in Figure 3. The specific experiment included the following steps:  Figure 3. Workflow of soil moisture retrieval.

Ratio Method
The ratio of bare soil backscattering contribution to the observed backscattering can be affected by vegetation cover and remote sensing sensor configuration [43]. Therefore, for areas with vegetation cover, when the sensor configuration parameters were fixed, the ratio of the surface backscattering to the total backscattering coefficients was only related to some vegetation description parameters, which can be expressed using the following (1) The model was constructed based on the observed total coefficients, the vegetation parameters (e.g., RVI, LAI and NDVI), and the bare soil backscattering coefficients simu-Remote Sens. 2022, 14, 3210 6 of 21 lated by the CIEM model. Then, the undetermined coefficients of the model were fitted using the least-squares method.
(2) For each sampling date, 70% of the data were randomly selected as the training set and the rest as the validation set. Basic statistics (e.g., mean, standard deviation, etc.) of the training and validation sets were calculated and compared, and they were not statically different.
(3) After fitting the undetermined coefficients of the model, the estimated bare soil backscattering could be obtained. Using a look-up table approach, the SMC of the training set was estimated based on the minimum cost function strategy (Equation (2)) and the simulated bare soil backscatter coefficients. Then the RMSE between the measured and estimated SMC was calculated as the estimation accuracy of SMC on the training set.
(4) All the RMS heights (s) within the given roughness parameter interval were traversed and steps (1) and (2) repeated until s reached the maximum value in the given range. The estimation accuracy of SMC on the training set was then compared with different s, and the s with the highest estimation accuracy was taken as the optimal roughness parameter of the study area on the sampling date. The optimal roughness was taken as the effective roughness of the model.
(5) For the validation set, the estimated backscattering coefficients of bare soil were obtained using the constructed model. Then the SMC of the validation set was estimated using the minimum cost function strategy and the backscattering coefficients of bare soil simulated by the CIEM model with the optimal roughness parameter. Finally, the accuracy evaluation was conducted using the measured SMC on the validation set. (6) The SMC mapping of the study area was conducted using the constructed model.

Ratio Method
The ratio of bare soil backscattering contribution to the observed backscattering can be affected by vegetation cover and remote sensing sensor configuration [43]. Therefore, for areas with vegetation cover, when the sensor configuration parameters were fixed, the ratio of the surface backscattering to the total backscattering coefficients was only related to some vegetation description parameters, which can be expressed using the following equation: where pp represents the HH or VV polarization mode, σ • ppsoil is the backscattering coefficient of bare soil, and σ • pp represents the total backscattering coefficient. The ratio function F on the right side of Equation (3) does not have a fixed expression form. Scholars have used a variety of ratio functions in order to better characterize the backscattering contribution of vegetation in the ratio equation [3,43,44]. A combination of the power exponential function and linear function could better characterize the vegetation scattering [45].
where V is the vegetation description parameter, a, b and c are the undetermined coefficients of the ratio function.
For the selection of vegetation description parameters in the ratio function, the current research has not reached a consistent conclusion. Some studies believed that the vegetation water content (VWC) had a strong attenuation effect on the scattered radiation of the microwave signal in the vegetation layer, so VWC was suitable for vegetation parameter description. There is a strong correlation between radar vegetation index (RVI) and VWC [46]. The RVI can be obtained directly from SAR images without being affected by weather conditions. Therefore, RVI was selected as the vegetation parameter to parameterize the ratio function in this study. The definition of RVI is shown as follows: where σ • HV , σ • HH and σ • VV are the backscattering coefficient of different polarization modes. In addition to RVI, leaves have been considered to be able to control the scattering and attenuation of microwave signals in the canopy to a certain extent [39], especially for herbaceous vegetation. Therefore, LAI was used to describe the scattering characteristics of vegetation. NDVI can well represent vegetation greenness information, reflecting a growth state of vegetation, and does not depend on field measurements which can be obtained from optical remote sensing data. Therefore, NDVI was also used to parameterize the ratio function. These three vegetation parameters were compared with for the SMC estimation performance. LAI was derived from field measurements, and NDVI was extracted from Sentinel-2. Since the data of LAI and NDVI did not cover all sampling dates, only the SMC of the sampling dates with data was estimated when using these two parameters.

WCM
The WCM is a semi-empirical backscattering model proposed by Attema and Ulaby [35]. The WCM divides the observed total backscattering into the three parts: the vegetation canopy scattering component, the multiple scattering component between vegetation and the soil surface, and the surface scattering component attenuated by the vegetation, which can be expressed as follows: where σ • can is the observed total backscattering coefficient, σ • veg is the scattering component from the vegetation canopy, σ • veg+soil is the multiple scattering component formed between vegetation and soil, and σ • soil is the surface scattering component attenuated by vegetation layer. τ 2 represents the two-way attenuation through the vegetation canopy. A and B are the undetermined coefficients of the model, which are generally obtained by leastsquared fitting according to the measured data. V 1 and V 2 represent the direct scattering and the attenuation properties of vegetation, respectively. In this study, V 1 and V 2 were parameterized by the three vegetation parameters (RVI, LAI and NDVI) that were used to construct the ratio function. Among the scattering components in the WCM model, σ • veg+soil accounts for a very small proportion compared to the remaining two components, which can usually be ignored.
Using the McLaughlin series to expand the parameter τ 2 , the WCM can be further simplified [47].
Ignoring the multiple scattering components between soil and vegetation, and retaining the first two terms of the McLaughlin series expansion, the WCM can be expressed as follows: Taking σ • can as the dependent variable and σ • soil as the independent variable, let a = 2AB and b = −2B/cosθ, Equation (8) can be simplified to the following expression: where V represents the vegetation description parameter. According to the simplified WCM model (Equation (9)), the functional relationship between the backscattering coefficient of bare soil and the total backscattering coefficient can be established after obtaining the undetermined coefficients.

CIEM Model
The integrated equation model (IEM) is a highly applicable bare soil scattering model [48], which constructs the functional relationship between co-polarized backscattering coefficients and surface roughness parameters (e.g., RMS height (s) and correlation length (L)), incident angle and wavelength, and soil dielectric constant. The field measurement of the surface RMS height (s) and correlation length (L) are greatly affected by human factors. Baghdadi et al. [40] constructed the functional expression between s and L based on a large number of roughness data from field surveys and C-band SAR data combined with simulated value of IEM model, and proposed a calibrated integrated equation model (CIEM), which reduced the dependence on s or L.
The empirical relationship between correlation length and RMS height in the CIEM model establish for C-band SAR can be expressed as follows: where L is the surface correlation length, s is RMS height, pp is the HH or VV polarization mode, and a, b, c and d are the undetermined coefficients. According to the different polarization modes, the values are as follows: a HH = 4.026, d HH = 1.551, a VV = 3.289, d VV = 1.222, b HH = b VV = −1.774, c HH = bc VV = −0.0025,

Optimal Roughness Parameters
The correlation length (L) in the CIEM is expressed as an empirical function of the RMS height (s), reducing the number of unknown roughness parameters. Therefore, only the s was used as the roughness parameter in this study. The roughness parameter is an important parameter in microwave scattering model of bare soil. However, the measurement process of roughness parameters is complex and easily affected by human factors and the measuring instruments used, especially in densely vegetated areas [49]. In this study, the optimal roughness parameter based on the global search method [50] was used to parameterize the bare soil scattering model. The idea of this method is to assume that the roughness parameters in the study area have constant values during a certain period of time. The solution process of the optimal roughness parameters adopts the following method: given a roughness parameter range, for this study, the range of the RMS height was set to 1 mm-30 mm and the step size was 1 mm. Then the bare soil scattering model was constructed using the roughness parameter within the given range to estimate SMC. The SMC estimation accuracy with different roughness conditions was obtained, and finally the RMS height corresponding to the highest estimation accuracy was taken as the optimal roughness parameter of the study area. For the estimation of the overall SMC in the study area, the optimal roughness parameter was essentially used as an empirical parameter of the bare soil scattering model to improve the estimation performance of the overall SMC, and there would be deviations from the real roughness parameters at different Remote Sens. 2022, 14, 3210 9 of 21 locations in the study area. Therefore, it could not represent the real roughness conditions of all sampling points.

Accuracy Assessment
The accuracies of different models were compared and evaluated using statistical indicators, including the coefficient of determination (R 2 , Equation (12)) and the root-meansquare error (RMSE, Equation (13)).
where n is the number of sample points; y i andŷ i represent the measured SMC value and estimated SMC value of the i-th sample site respectively; and y is the mean value of the measured SMC values.

Optimal Roughness Parameter Selection
The SMC estimation accuracies on the training set with different RMS heights and different vegetation parameters were obtained (Figures 4 and 5). Figures 4 and 5 show the RMSE of the estimation results on the training set of different sampling dates when the model was the ratio method and WCM, respectively. When the RMSE reached the minimum, the RMS height corresponding to the horizontal axis was considered to be the optimal roughness parameter in the current phase. Table 2 shows the fitting coefficients of the two models (Ratio method and WCM) with the condition of optimal roughness parameters. Tables 3 and 4 show the optimal roughness parameters of the study area for different sampling dates when using the ratio method and WCM.
As shown in Figures 4 and 5, the SMC estimation accuracy was jointly affected by the vegetation parameters and the RMS height. For the ratio method, no matter which vegetation parameters was used as the input parameter of the model, the RMSE curve on the training set is generally parabolic, and it shows a trend of rapid decline at first and then a much gentler trend (See Figure 4). It shows that the estimation accuracy of SMC was more sensitive to the change of the RMS height when the RMS height was small.
As shown in Table 3, when different vegetation parameters were used to parameterize the ratio method, different values of the optimal roughness parameter were obtained. It is noteworthy that the values of optimal roughness parameters obtained on most sampling dates were large. Figure 4 shows that RMSE decreases with an increase in optimal roughness parameters. The variation of backscattering coefficients simulated by the CIEM with the RMS heights is not linear, the RMSE on the training set estimated by the ratio method was not strictly monotonically decreasing. The RMSE of estimated SMC would reach the global minimum with a smaller roughness parameter on some dates. Remote Sens. 2022, 14, x FOR PEER REVIEW 10 of 22  4, comparing with the ratio method, the optimal roughness of WCM was smaller. It suggested that the optimal roughness parameter not only changed with the vegetation parameters and the acquisition time of the SAR data, but also depended on the models used in the estimation process to describe the scattering characteristics of vegetation. In general, the optimal roughness parameter, as an empirical coefficient in the scattering model, can effectively improve the estimation performance of the SMC.       As shown in Figure 5, For WCM, the RMSE variation curves estimated on the training set are significantly different from that obtained from the ratio method. When RVI or LAI was used the vegetation parameter, the estimated RSME curves shows an upward trend on most sampling dates. With the gradual increase of the RMS height, the change of the RMSE tended to be gentle, the sensitivity of estimation accuracy to the change in RMS height decreased, and the estimation accuracy reached the highest when the roughness was low. When NDVI was used as vegetation description parameter, the RMSE curves of the training set on different dates almost tended to be horizontal; that is, the estimation accuracy did not change significantly with the increase in RMS height. As shown in Figure 5h, for the SMC estimation of the training set on July 10, the estimation performance with three vegetation parameters has little change with different RMS heights. This is because the value of measured SMC on July 10 was low, and the backscattering of the bare soil simulated by CIEM is not sensitive to the change of the RMS height, so that the estimated SMC values with different roughness conditions were almost the same. As shown in Table 4, comparing with the ratio method, the optimal roughness of WCM was smaller. It suggested that the optimal roughness parameter not only changed with the vegetation parameters and the acquisition time of the SAR data, but also depended on the models used in the estimation process to describe the scattering characteristics of vegetation. In general, the optimal roughness parameter, as an empirical coefficient in the scattering model, can effectively improve the estimation performance of the SMC. Figures 6-11 show the scatter plots between the measured and the estimated SMC on the validation set from ratio method or WCM using different vegetation parameters. Table 5 shows the estimation accuracy of SMC on the validation set. Since the data of LAI and NDVI did not cover all sampling dates, when comparing the estimation accuracy under different vegetation parameters, only the sampling dates that contained vegetation parameters data were compared.            As shown in Figures 6-11 and Table 5, the accuracy of SMC retrieval was the lowest when RVI was used as vegetation descriptor (for ratio method, R 2 = 0.48, RMSE = 5.90 vol.%; for WCM R 2 = 0.34, RMSE = 7.17 vol.%). The accuracy of SMC retrieval was the highest when NDVI was used as vegetation descriptor (for ratio method, R 2 = 0.65, RMSE = 4.35 vol.%; for WCM R 2 = 0.64, RMSE = 4.33 vol.%). In order to compare the estimation performances of models parameterized by RVI, LAI and NDVI, the SMC estimation accuracies were studied. Since NDVI data were not available for 9 May and 20 May and LAI data were not available for 6 June, these dates were excluded from the comparison. As shown in Figures 6-11 and Table 5, the performances of using LAI and NDVI were better than using RVI. As shown in Table 5, the estimation accuracy based on NDVI was higher than that based on LAI.         The estimation performances based on NDVI were better than other vegetation parameters, which may be because healthy vegetation has unique absorption and reflection characteristics for the red and near-infrared (NIR) band. The NDVI calculated by optical data could better describe the vegetation parameter in the model. As shown in Table 5, comparing the results of ratio method and WCM, it suggests that there is no significant difference in estimation accuracy of the two models when NDVI was used to parameterize the models. When using RVI or LAI, the estimation accuracy of the ratio method is higher than that of WCM. In general, the performance of the ratio method is better than WCM.

Soil Moisture Estimation Results
In addition, for a specific study area, the SMC estimation accuracy would be affected by vegetation parameters. Selecting appropriate vegetation parameters would help to improve the estimation performance of SMC.

Reference Incidence Angle
The difference in distance between the ground object and the sensor would cause a difference in the incident angle, and the local incident angle would also be affected by the terrain. Since polarization signature can be affected by incidence angle, a normalization was performed at a reference incidence angle to reduce the effect of local incidence angle effect on the backscatter coefficients. In this study, the normalization method of incidence angle proposed by Ulaby et al. [41] was used to reduce the effect of local incidence angle on backscatter coefficients. The reference incidence angle was set to 30 • . However, there was no unified conclusion on the selection of the reference incidence angle. In order to analyze the influence of the reference incidence angle on the estimation accuracy of the SMC estimation proposed in this study, the range of the reference incidence angle was 20 • to 40 • , and the step size was set to 1 • . The SMC on the validation set was retrieved by using the technical route constructed in Figure 3. When using different vegetation indices with different reference incidence angles, R 2 and RMSE between measured and estimated SMC values on the validation are shown in Figures 12 and 13. The incidence angle corresponding to the minimum RMSE was taken as the optimal reference incidence angle. Table 6 shows the SMC estimation accuracy of the validation set when the optimal reference incidence angle was used. with different reference incidence angles, R 2 and RMSE between measured and estimated SMC values on the validation are shown in Figures 12 and 13. The incidence angle corresponding to the minimum RMSE was taken as the optimal reference incidence angle. Table 6 shows the SMC estimation accuracy of the validation set when the optimal reference incidence angle was used.  As shown in Figure 12, when estimating SMC based on the ratio method, the reference incidence angle had a significant impact on the SMC estimation results. The estimation accuracy on the validation set generally increased first and then decreased when the RVI, LAI, and NDVI served as vegetation descriptors. Meanwhile, it could be seen that, when RVI was used as vegetation descriptor, the variations in R 2 and RMSE curves are larger than that when the other two vegetation parameters were used. This may be due to the change of the simulated backscattering coefficients of the bare soil when the reference incidence angle changed. The RVI was also changed due to the influence of the normalization of the incidence angle. The parameters of dependent and independent variable changed when constructing the model. When LAI or NDVI serves as the descriptor, only the simulated backscattering coefficients of the bare soil changed, and the independent variable did not change. Therefore, when RVI was used as a vegetation description, RVI was more affected by the reference incidence angle than LAI or NDVI.   As shown in Tables 5 and 6, the estimation accuracies of SMC were improved when using the optimal reference incident angle. Taking the ratio method combined with RVI as an example, it is known in the previous analysis that when the reference incident angle was 30 • , the estimation performance based on RVI was the worst (R 2 = 0.48, RMSE = 5.90 vol.%). When using the optimal reference angle, the estimation accuracy of SMC was improved to R 2 = 0.56, RMSE = 5.10 vol.%. Therefore, when the models combined the optimal roughness parameter and the vegetation parameter for SMC retrieval, the reference incidence angle can also be used as another adjustment parameter to help improve the estimation accuracy.

Regional Soil Moisture Mapping
In this study, LAI data were only available for a limited number of sampling points. Although NDVI could be extracted from optical images, it could not cover all sampling dates due to cloud interference in the study area. Therefore, only RVI was used as vegetation description parameter to construct the models for SMC mapping in this study. The SMC maps are shown in Figure 14. In general, the spatial patterns of the SMC maps obtained by the two models are consistent, and they conform to the trend of the measured SMC throughout the growing season. Comparing Figure 14a,b, the soil moisture values retrieved by WCM are higher in the area with high SMC, reaching the critical value of 50 vol.% in most areas, which shows an apparent overestimation compared with the measured soil moisture. In contrast, there is no obvious overestimation area in the SMC map based on the ratio method. On the last sampling date, both models showed an obvious underestimation. This is mainly because the SMC measured on July 10 was low, and the CIEM model is not sensitive to soil changes when the SMC is low.

Limitations and Potential Improvements
Due to various constrains, soil moisture measurements did not cover the entire whe growing season in this study. Further experiments can be conducted in the future to co

Limitations and Potential Improvements
Due to various constrains, soil moisture measurements did not cover the entire wheat growing season in this study. Further experiments can be conducted in the future to collect and analyze soil moisture information for the whole growing season of wheat to evaluate the effectiveness and applicability of the method proposed in this study. The effectiveness of the proposed methods for other crop types is also worth further verification in future studies.
From a spatial coverage point of view, the study area in this study was limited to a single wheat field, so the optimal roughness parameter was used to characterize the roughness condition of the whole study area. When conducting soil moisture inversion over a large agricultural area, the root-mean-square height or correlation length of the surface can be obtained by inverting the surface roughness, which and then be used to construct an improved surface scattering model.

Conclusions
In this study, based on RADARSAT-2 SAR data and the measured SMC information, we investigated the potential of a simple method combined with three different vegetation parameters to estimate SMC under vegetation cover. In addition, the WCM was also used to retrieve SMC for comparison. Based on the study, the following conclusions can be drawn: (1) As an empirical coefficient in the scattering model, the optimal roughness parameter could effectively improve the estimation accuracy of SMC. (2) A reference incident angle could be used as an effective adjustment parameter to improve the estimation accuracy of SMC. (3) Comparing with LAI and RVI, NDVI is more suitable as the vegetation description parameter of the model. (4) In general, the estimation performance of the ratio method is better than that of WCM.
Despite the encouraging results of this study, there are certain limitations that need to be addressed in future studies. The measured data in this study only covered part of the wheat growing period. In the future, it will be necessary to evaluate the effectiveness of the proposed method over longer wheat growth cycles. In addition, the same method should be tested on other crop types to further verify its effectiveness.

Conflicts of Interest:
The authors declare no conflict of interest.