Retrieval of Soil Moisture by Integrating Sentinel-1A and MODIS Data over Agricultural Fields

Soil moisture is an essential variable in the land surface ecosystem, which plays an important role in agricultural drought monitoring, crop status monitoring, and crop yield prediction. High-resolution radar data can be combined with optical remote-sensing data to provide a new approach to estimate high-resolution soil moisture over vegetated areas. In this paper, the Sentinel-1A data and the Moderate Resolution Imaging Spectroradiometer (MODIS) data are combined to retrieve soil moisture over agricultural fields. The advanced integral equation model (AIEM) is utilized to calculate the scattering contribution of the bare soil surface. The water cloud model (WCM) is applied to model the backscattering coefficient of vegetated areas, which use two vegetation parameters to parameterize the scattering and attenuation properties of vegetation. Four different vegetation parameters extracted from MODIS products are combined to predict the scattering contribution of vegetation, including the leaf area index (LAI), the fraction of photosynthetically active radiation (FPAR), normalized difference vegetation index (NDVI), and the enhanced vegetation index (EVI). The effective roughness parameters are chosen to parameterize the AIEM. The Sentinel-1A and MODIS data in 2017 are used to calibrate the coupled model, and the datasets in 2018 are used for soil moisture estimation. The calibration results indicate that the Sentinel-1A backscattering coefficient can be accurately predicted by the coupled model with the Pearson correlation coefficient (R) ranging from 0.58 to 0.81 and a root mean square error (RMSE) ranging from 0.996 to 1.401 dB. The modeled results show that the retrieved soil moisture can capture the seasonal dynamics of soil moisture with R ranging from 0.74 to 0.81. With the different vegetation parameter combinations used for parameterizing the scattering contribution of the canopy, the importance of suitable vegetation parameters for describing the scattering and attenuation properties of vegetation is confirmed. The LAI is recommended to characterize the scattering properties. There is no obvious clue for selecting vegetation descriptors to characterize the attenuation properties of vegetation. These promising results confirm the feasibility and validity of the coupled model for soil moisture retrieval from the Sentinel-1A and MODIS data.


Introduction
Soil moisture plays an important role in the water and energy cycles of the land surface ecosystem and is considered as an important variable in Earth Science [1,2]. In agriculture, soil moisture is a critical component, which is essential for many applications, such as crop status monitoring, crop yield prediction, agricultural drought monitoring, and water resource management [3,4]. Due to the dynamic variability of soil moisture in time and space, dense spatial and temporal measurements at a regional and global scale remain difficult. Fortunately, remote sensing provides an operational method for soil moisture monitoring and estimation. Due to the high sensitivity to soil moisture and the ability of all-time and all-weather observation, microwave remote sensing has been widely applied Water 2020, 12 in soil moisture retrieval over the bare or vegetated land surface [5,6]. This includes the synthetic aperture radar (SAR), which, due to its high spatial resolution, can be potentially used for providing soil moisture products with high spatiotemporal resolution. However, the backscattering coefficient measured by SAR was affected by many factors, such as the sensor parameters (frequency, polarization, and incidence angle), soil parameters (surface roughness, soil dielectric constant, and soil texture), and the vegetation parameters (vegetation water content, vegetation coverage, vegetation structure, and others) [7,8]. For bare soil, some surface scattering models have been proposed to simulate the backscattering coefficient and retrieve soil moisture, such as the Oh model [9], Dubois model [10], Shi model [11], and the advanced integral equation model (AIEM) [12,13]. The common problem faced by these models in practice is how to parameterize the surface roughness parameters. Until now, many methods have been proposed for eliminating the effects of surface roughness parameters. Bai et al. [14] utilized the in-situ surface roughness parameters to parameterize the AIEM. Zribi and Dechambre [15] proposed a combined roughness parameter based on multi-angle Radarsat-1 data for simultaneously retrieving the soil moisture and roughness parameters. Bai and He [16] utilized the VV and HH co-polarization to remove the root mean square (RMS) height based on the Dubois model. Su et al. [17] proposed the effective roughness parameters to reduce the unknowns for soil moisture retrieval. With time-series Sentinel-1A data, Bai et al. [4] used the effective roughness parameters to parameterize the AIEM for soil moisture retrieval. Among these models and methods, the AIEM and effective roughness parameters have been widely used for soil moisture retrieval, which has achieved satisfactory results for bare soil surface.
The above-surface scattering models were proposed based on the interaction between the microwave signal and the bare soil surface, which have good performance for the case of a bare soil surface. In vegetated areas, the canopy will produce complex volume scattering and therefore reduce the sensitivity of the radar signal to the soil moisture [8,18]. It is necessary to eliminate the scattering effect of the vegetation canopy for soil moisture retrieval in vegetated areas. To date, some semi-empirical and theoretical models have been proposed to model the backscatter from the canopy. For example, the water cloud model (WCM) considers the vegetation scattering, interaction between the vegetation and the soil surface, and the direct soil backscatter [19]; the ratio method assumes that the ratio between the bare soil and measured backscattering coefficient is only dependent on the vegetation [20,21]; the Tor Vergata model considers the vegetation as an ensemble of discrete lossy scatterers, whose electromagnetic behaviors can be simulated by simple geometrical shapes [6,22]. The implementation of the aforementioned vegetation-scattering models requires parameterizing the scattering contribution of vegetation. From the perspective of operability and practicality, the WCM has been demonstrated to be suitable for soil moisture estimation over grassland and agricultural areas [5,[23][24][25][26][27][28][29].
The WCM needs two vegetation parameters for parameterizing the scattering contribution of vegetation. In practical applications, the descriptors are often parameterized by the optical vegetation parameters, such as the leaf area index (LAI) [5,27], vegetation water content (VWC) [23,24], normalized difference vegetation index (NDVI) [14,25], enhanced vegetation index (EVI) [14,25], fraction of photosynthetically active radiation (FPAR) [26], leaf water area index (LWAI) [27], and fuel moisture content [28]. Bai et al. [14,25] investigated some related works and found that the performance of these optical vegetation parameters for helping soil moisture retrieval is heavily related with the type of study area, such as vegetation type and vegetation growth state (dense or sparse). Qiu et al. [29] examined the impact of vegetation parameters from two different satellites on soil moisture retrieval. In the aforementioned works, the scattering and attenuation of vegetation were parameterized by the same vegetation parameters. He et al. [26] found that the scattering and attenuation properties should better be parameterized by different vegetation parameters. From the perspective of theoretical analysis, Park et al. [30] found that the particle moisture content and the VWC can be the optimal vegetation descriptors. These studies demonstrated that the optical vegetation parameters can accurately estimate the scattering contribution of vegetation. The integration of SAR data and optical remote-sensing Water 2020, 12, 1726 3 of 18 data may provide operational methods for soil moisture retrieval in vegetated areas. Multiple studies have focused on soil moisture retrieval by assuming that the scattering and attenuation properties are characterized by the same vegetation indicators [5,23,24]. However, relatively few studies have evaluated the performance of soil moisture retrieval while the vegetation-scattering and attenuation properties are expressed by different vegetation parameters, which may improve the accuracy of soil moisture estimation.
The objective of this paper was to assess the capability of different vegetation parameters for characterizing the scattering and attenuation properties during soil moisture retrieval in agricultural fields by the integration of Sentinel-1A data and the Moderate Resolution Imaging Spectroradiometer (MODIS) data. The physically based AIEM and semi-empirical WCM models are selected to obtain a coupled scattering model. Four vegetation parameters extracted from the MODIS products, namely LAI, FPAR, NDVI, and EVI were chosen as vegetation descriptors to characterize the scattering and attenuation properties of vegetation, respectively. Meanwhile, the whole study period was divided as a period with bare soil surface and a period with vegetation cover by refereeing to the LAI values. The structure of this paper is organized as follows. In Section 2, the study area and remote-sensing data, including Sentinel-1A and MODIS data, are introduced. In Section 3, the scattering models of AIEM and WCM are reviewed, and the calibration for the coupled mode and soil moisture retrieval is presented. The calibration results for the coupled model, soil moisture retrieval results, and error analysis are described in Section 4. Finally, the conclusions are drawn in Section 5.

Study Area and In-Situ Measurements
In this paper, the REMEDHUS soil moisture monitoring network (41 • 8 56"-41 • 27 21" N, 5 • 13 29"-5 • 35 31" W) in the central Douro Basin of Spain was selected as the test site ( Figure 1) [31], which covers an area of approximately 1300 km 2 with flat terrain. The altitude ranges from 700 to 900 m and can be classified as a continental semi-arid Mediterranean climate. The annual precipitation is about 385 mm, and the maximum monthly average precipitation occurs in May, with about 47 mm rainfall. and the minimum case occurs in August, with about 11 mm rainfall. The average temperature is 12 • C and the average potential evaporation is 908 mm. In this area, the main types of vegetation covers include rain-fed cereals in winter and spring (78%), irrigated crops (5%) and vineyards (3%) in summer, and some scattered forest and pastures (13%). The growth cycle of rain-fed cereals includes the sowing period in autumn, the development period in spring, and the maturity period in early summer.
From 1999 to 2009, the REMEDHUS ground monitoring network has set up 23 ground monitoring stations. The distribution of stations in this area is uneven based on the division of regional geomorphic units. The monitoring instrument of soil moisture was the Stevens Hydra probe, with the monitoring depth of 0-5 cm, and the data were automatically collected every hour. The soil moisture data collected by the REMEDHUS ground monitoring network were applied to many research fields, such as the calibration and validation of the Soil Moisture and Ocean Salinity (SMOS) soil moisture products [31], the validation of remote-sensing soil moisture products [32], and the scale transformation of soil moisture products [33]. The in-situ soil moisture data of REMEDHUS were downloaded from the international soil moisture network (station Las Bodegas). From the time series of the daily average soil moisture at Las Bodegas from 01/01/2014 to 28/12/2018 (Appendix A Figure A1), we found that the soil moisture presents periodic seasonal variation and that there is no obvious distinction between dry and wet years. It is deemed that one-year data covered the entire growth cycle of crops and were able to effectively validate the coupled model and retrieve soil moisture. Therefore, in-situ soil moisture data collected from 01/01/2017 to 31/12/2018 were used, which are presented in Figure 2. It is assumed that the measured soil moisture is representative of the fields studied. It can be seen from this figure that the soil moisture content in the coverage area of the station Las Bodegas is relatively low (below 0.25 m 3 /m 3 ) and there are obvious fluctuations with the seasons. The soil moisture is significantly increased in May due to abundant rainfall and decreased in August due to less rainfall. were used, which are presented in Figure 2. It is assumed that the measured soil moisture is representative of the fields studied. It can be seen from this figure that the soil moisture content in the coverage area of the station Las Bodegas is relatively low (below 0.25 m 3 /m 3 ) and there are obvious fluctuations with the seasons. The soil moisture is significantly increased in May due to abundant rainfall and decreased in August due to less rainfall.

Sentinel-1A Data and Preprocessing
The Sentinel-1 mission, the Earth Radar Observation for the Copernicus joint initiative of the European Commission and the European Space Agency, comprises a constellation of two satellites, operating day and night and performing C-band SAR imagining. The Sentinel-1A is the first satellite of the European Space Agency (ESA)'s Copernicus plans, which was launched in April 2014. In this paper, two years of Sentinel-1A data with descending orbit were downloaded in the study area from were used, which are presented in Figure 2. It is assumed that the measured soil moisture is representative of the fields studied. It can be seen from this figure that the soil moisture content in the coverage area of the station Las Bodegas is relatively low (below 0.25 m 3 /m 3 ) and there are obvious fluctuations with the seasons. The soil moisture is significantly increased in May due to abundant rainfall and decreased in August due to less rainfall.

Sentinel-1A Data and Preprocessing
The Sentinel-1 mission, the Earth Radar Observation for the Copernicus joint initiative of the European Commission and the European Space Agency, comprises a constellation of two satellites, operating day and night and performing C-band SAR imagining. The Sentinel-1A is the first satellite of the European Space Agency (ESA)'s Copernicus plans, which was launched in April 2014. In this paper, two years of Sentinel-1A data with descending orbit were downloaded in the study area from

Sentinel-1A Data and Preprocessing
The Sentinel-1 mission, the Earth Radar Observation for the Copernicus joint initiative of the European Commission and the European Space Agency, comprises a constellation of two satellites, operating day and night and performing C-band SAR imagining. The Sentinel-1A is the first satellite of the European Space Agency (ESA)'s Copernicus plans, which was launched in April 2014. In this paper, two years of Sentinel-1A data with descending orbit were downloaded in the study area from 1  Table A1) [34]. The specific parameters for the Sentinel-1A images are listed in Table 1. The batch pre-processing of the Sentinel-1A data was based on Python coding, which included applying orbit file, thermal noise removal, radiometric calibration, speckle filtering (Lee sigma with window sizes of X and Y were both 7 and the sigma value was 0.9), geometric corrections (Range-Doppler terrain correction), radiometric normalization, resample (1 km with bilinear sampling technique), and re-projection (WGS 1884 Universal Transverse Mercator coordinates).

MODIS Data and Preprocessing
To characterize the scattering and attenuation properties of vegetation, four different vegetation parameters were extracted from the MODIS products [35], namely the LAI, FPAR, NDVI, and the EVI, which were used to describe the growing status of vegetation.
The LAI is defined as the one-sided green leaf area per unit ground area in broadleaf canopies and as on-half the total needle surface area per unit ground area in coniferous canopies [36,37]. FPAR is defined as the fraction of incident photosynthetically active radiation absorbed by the green elements of a vegetation canopy [37]. The MCD15A3H Collection 6 (DOI: http://doi.org/10.5067/ MODIS/MOD15A3H.006) products included LAI and FPAR, which was composited every 4 days at 1 km spatial resolution [38]. A main look-up-table (LUT)-based procedure and a back-up algorithm were used to obtain the LAI/FPAR from the surface reflectance or vegetation indices [37]. Based on ground truth, the LAI products have been assessed over widely distributed fields, and the FPAR products have been estimated using relatively few independent measurements.
The NDVI/EVI is directly computed from surface reflectance, which is developed to provide vegetation conditions. Global MYD13A2 Collection 6 (DOI: http://doi.org/10.5067/MODIS/MYD13A2. 006) products included the NDVI and EVI, which was composited every 16 days at 1 km spatial resolution [39]. The vegetation indices products were validated over different study areas and time periods based on ground truth.
The preprocessing for the MODIS products was conducted with the following procedures. The original 1 km sinusoidal projection grid data were projected onto the WGS 1984 UTM coordinates. A third-order Savitzky-Golay filter was applied to smooth the original MODIS products and only the "good quality" data were used to reduce the influence of cloud on vegetation parameters [40]. The cubic spline interpolation technique was utilized to obtain the vegetation parameters on the Sentinel-1A acquisitions. Figure

AIEM
The AIEM was implemented to calculate the backscattering coefficient of the bare soil surface [12,13], which describes the backscattering behavior of a random rough surface with a very wide range of validity. This model has been applied and validated in many works [5,20,21,24,25,27,29]. The inputs and outputs of AIEM are conceptually described as follows: where f is the frequency, θ stands for the incidence angle, pp indicates the polarization mode (VV), ε is the soil dielectric constant, s is the root mean squares (RMS) height, l is the correlation length, ACF

AIEM
The AIEM was implemented to calculate the backscattering coefficient of the bare soil surface [12,13], which describes the backscattering behavior of a random rough surface with a very wide range of validity. This model has been applied and validated in many works [5,20,21,24,25,27,29]. The inputs and outputs of AIEM are conceptually described as follows: where f is the frequency, θ stands for the incidence angle, pp indicates the polarization mode (VV), ε is the soil dielectric constant, s is the root mean squares (RMS) height, l is the correlation length, ACF represents the auto-correlation function, and σ o soil (dB) is the backscattering coefficient of the bare soil surface. The Dobson dielectric mixing model was applied to determine the relationship between the soil moisture and the soil dielectric constant [41]. The texture parameters refer to the harmonized world soil database (HWSD), which was created by the Food and Agricultural Organization (FAO) and an International Institute for Applied System Analysis (IIASA). The topsoil composition for the station Las Bodegas is composed of about 36.0% sand, 21.0% clay, and 1.41 g/cm 3 bulk density.

WCM
The WCM, a semi-empirical backscatter model, was conducted to model the backscatter of vegetated areas [19]. This semi-empirical model has been widely applied for soil moisture retrieval due to its efficient performance. For a given polarization, it divides the whole backscatter into three parts, the vegetation-scattering contribution σ o veg (dB), double-bounce scattering components between the vegetation and the underlying soil surface σ o veg+soil (dB), and the direct soil backscatter attenuated by the vegetation layer σ o soil (dB), which can be written as follows: where τ 2 stands for the two-way attenuation through the vegetation, θ is the incidence angle, V 1 and V 2 are the vegetation descriptors, which represent the scattering and attenuation properties of vegetation, respectively, and A and B are the model coefficients which depend on the vegetation descriptor and sensor configuration. In this paper, the WCM was simplified by one common assumption that the interaction between the vegetation and the underlying soil can be neglected: The V 1 and V 2 are cross-parameterized by the LAI, FPAR, NDVI, and the EVI, respectively. The backscattering coefficient of bare soil in Equation (5) was calculated by the AIEM combined with the Dobson model.

Model Calibration and Soil Moisture Retrieval
In this paper, the WCM and AIEM were coupled to predict the backscattering coefficient and retrieve the soil moisture from the Sentinel-1A data and MODIS data. The WCM is used for simulating the vegetation backscatter and AIEM is used for simulating the soil backscatter. The coupled model was parameterized by the effective roughness parameters, which have been proved to be effective for parameterizing the surface roughness parameters [5,25]. The coupled model calibration and soil moisture retrieval were implemented by the following steps (Figure 4).
Step 1: The study period was distinguished as a vegetated cover period and bare soil surface period. The LAI was used to separate the whole study period as a vegetation-covered period when the LAI was larger than 0.4 m 2 /m 2 , and the bare soil surface period when the LAI is less than 0.4 m 2 /m 2 .
Step 2: The backscattering coefficient of the bare soil surface was estimated. The AIEM and Dobson model were implemented to calculate the backscattering coefficient of the bare soil surface with a series of input parameters. The frequency and incidence angle were set as 5.405 GHz and 40 • , which are the frequency and incidence angle of Sentinel-1A. The s-values ranged from 0.5 to 2.0 cm with an increment of 0.1 cm. The l-values ranged from 5.0 to 20 cm with an increment of 1 cm. The soil moisture ranged from 0.01 to 0.40 m 3 /m 3 with an interval of 0.01 m 3 /m 3 , which was set according to the maximum and minimum values of the soil moisture in the study area. Provided the input parameters were correct, the AIEM could predict the backscattering coefficient for the bare soil surface.
Step 3: The backscattering coefficient of the vegetation was simulated. The WCM was conducted to model the whole backscattering coefficient of the vegetated areas with a simulated backscattering coefficient of the bare soil surface. The WCM was parameterized by the LAI, FPAR, NDVI, and EVI, respectively. The model coefficients A and B were solved using the least square method.
Step 4: Establishment of the cost function S. The Sentinel-1 data and MODIS data acquired in 2017 were used for model calibration. During the period with the bare soil surface, the backscattering coefficient was calculated by the AIEM. During the period with vegetation cover, the backscattering coefficient was simulated by the combination of WCM and AIEM. The following cost function S was defined to minimize the Sentinel-1A observations and simulations with the coupled model, which was used for calibrating the model coefficient (A and B) and the effective roughness parameters (s and l): where σ o sim and σ o S1 indicate the simulated backscattering coefficients and the Sentinel-1A observations, and n is the number of the sampling period.
Step 5: The s varied from 0.5 to 2.0 cm, and the l varied from 5.0 to 20 cm. A 16 by 16 matrix was obtained, which was filled with the computed S-values. The index position of the minimum S was searched within this 16 by 16 matrix. The roughness parameters (s and l) corresponding to the searched index position were considered to be the effective roughness parameters. At this moment, the coupled model was calibrated with determined roughness parameters (s and l) and model coefficients (A and B).
Step 6: Soil moisture was retrieved. The soil moisture was retrieved using the calibrated coupled model from the Sentinel-1A data and the MODIS data acquired in 2018. The LUT method was used to retrieve the soil moisture based on the database produced in Step 2.  Step 3: The backscattering coefficient of the vegetation was simulated. The WCM was conducted to model the whole backscattering coefficient of the vegetated areas with a simulated backscattering coefficient of the bare soil surface. The WCM was parameterized by the LAI, FPAR, NDVI, and EVI, respectively. The model coefficients A and B were solved using the least square method.
Step 4: Establishment of the cost function S. The Sentinel-1 data and MODIS data acquired in Step 7: Accuracy evaluation: the accuracy of the coupled model calibration and soil moisture retrieval was evaluated using the following four statistical indicators, namely the mean bias (BIAS), mean absolute error (MAE), root mean square error (RMSE), and the Pearson correlation coefficient (R): where V obs and V est stand for the measured backscatter (or soil moisture) and the estimated backscatter (or soil moisture), and n is the number of sampling periods.

Calibration Results for the Coupled Model
The Sentinel-1A and MODIS datasets acquired in 2017 are used for model calibration, and the datasets collected in 2018 are used for soil moisture retrieval. The calibration of the coupled model (WCM and AIEM) is conducted by following the procedures in Section 3.3. The effective roughness parameters (s and l), WCM coefficients (A and B), and model calibration accuracy are given in Table 2. Figures 5 and 6 present the Sentinel-1A observations and the simulated backscattering coefficients using the coupled model based on the effective roughness parameters.
coupled model from the Sentinel-1A data and the MODIS data acquired in 2018. The LUT method was used to retrieve the soil moisture based on the database produced in Step 2.
Step 7: Accuracy evaluation: the accuracy of the coupled model calibration and soil moisture retrieval was evaluated using the following four statistical indicators, namely the mean bias (BIAS), mean absolute error (MAE), root mean square error (RMSE), and the Pearson correlation coefficient where obs V and est V stand for the measured backscatter (or soil moisture) and the estimated backscatter (or soil moisture), and n is the number of sampling periods.

Calibration Results for the Coupled Model
The Sentinel-1A and MODIS datasets acquired in 2017 are used for model calibration, and the datasets collected in 2018 are used for soil moisture retrieval. The calibration of the coupled model (WCM and AIEM) is conducted by following the procedures in Section 3.3. The effective roughness parameters (s and l), WCM coefficients (A and B), and model calibration accuracy are given in Table  2. Figures 5 and 6 present the Sentinel-1A observations and the simulated backscattering coefficients using the coupled model based on the effective roughness parameters.  From Table 2, it can be seen that all the values of BIAS are positive. From Figure 5, it is found that the location of the mean for the Sentinel-1A observations is higher than the values of the backscattering coefficient simulated by the coupled model with different vegetation parameters, which indicates that the simulated backscattering coefficients are slightly underestimated. This phenomenon is very obvious in the periods between 5 November 2017 and 11 December 2017. From Figure 5, it can be seen that the backscattering coefficient observed by Sentinel-1A in these two dates is relatively high compared with other dates, which is seen as outliers. The maximum and minimum values of MAE are 1.081 and 0.769 dB, respectively. The values of RMSE range from 0.996 to 1.401 dB. The Pearson correlation coefficient (R) has a minimum value of 0.58 and a maximum value of 0.81. From Figure 6, it is found that the simulated backscattering coefficient is consistent with Sentinel-1A observations, and the dynamics of backscattering coefficients are captured. From the perspective of statistical metrics, it can be said that the coupled model can predict the backscattering coefficient in the study area. percentiles (top and bottom of the blank box, respectively), and the outlier ranges (edges of the top and bottom whiskers). The "LAI_FPAR" indicates that the V1 and V2 are characterized by the LAI and FPAR, respectively.   It is noted that the calibrated results are similar between the different vegetation parameters. This is mainly because there is a certain correlation between the selected vegetation parameters (see Appendix A Table A2). The R between the NDVI and EVI is 0.97 and the value between the LAI and the FPAR is 0.91, which means there is a very strong correlation between them. There is a moderate correlation between the LAI and NDVI (0.59), the LAI and EVI (0.68), the FPAR and NDVI (0.63), and the FPAR and EVI (0.67). Therefore, the optimum parameters, statistical results, and backscattering coefficient dynamics are similar for the coupled model. When the LAI is chosen as V 1 , the statistical results are improved. The values of R are 0.75, 0.76, and 0.81 when V 2 is parameterized by the FPAR, NDVI, and the EVI, respectively. When the LAI is chosen as V 2 , the statistical results are relatively poor. The values of R are 0.61, 0.58, and 0.60 when V 1 is selected as the FPAR, NDVI, and the EVI, respectively. The calibration results for other combinations of the vegetation parameters are atypical. When V 1 and V 2 are chosen as LAI and EVI, the calibration result is optimized with R equal to 0.81. When V 1 and V 2 are parameterized by two different vegetation parameters, the statistical results are better than the results obtained from the case with the same vegetation parameters. Based on the above discussion, it can be concluded that the LAI is recommended to parameterize the scattering properties of vegetation, and the optimum combination of vegetation parameters should be considered for parameterizing the coupled model.
It is worth noting that most of the effective roughness parameters selected for the coupled model fall within the ranges of the defined ranges of surface roughness parameters, which also verified the rationality of the defined range. Simultaneously, it was found that the selected effective roughness parameters derived from each vegetation combination exhibited no discernible difference. This phenomenon may ascribe to the strong correlation between the vegetation parameters. The effective roughness parameters selected as inputs for the AIEM has accurately predicted the measured backscattering coefficients, which confirm the validity and effectiveness of this method.

Soil Moisture Retrieval Results
Based on the calibrated model, the soil moisture is estimated from the Sentinel-1A and MODIS data collected in 2018. The time series of the in-situ and retrieved soil moisture are displayed in Figures 7 and 8, and the statistical indicators are listed in Table 3.     Table 3. Statistical metrics between the measured and retrieved soil moisture.    Table 3, it is found that the accuracy of the soil moisture estimation is nearly the same when the different vegetation parameter combinations are used. The average values of the BIAS, MAE, and the RMSE are 0.018, 0.048, and 0.055 m 3 /m 3 . The positive BIAS indicates that the retrieved soil moisture is systematically underestimated compared with observations, which can be easily seen in Figure 7. After 15 July 2018, the retrieved soil moisture is generally underestimated. Before 15 July 2018, the retrieved soil moisture is generally overestimated. This phenomenon may be ascribed to the seasonal rainfall (more rainfall in May and less rainfall in August). On 9 February 2018, the retrieved soil moisture reached the critical value of the database, which is a bit abnormal. The backscattering coefficient on this date is very low compared with other days, which may have been caused by some unknown reason, other than the soil moisture and vegetation. The minimum and maximum values of R are 0.74 and 0.81, respectively. From Figure 8, it can be seen that the trend of retrieved soil moisture is consistent with the in-situ observations. The seasonal dynamics of the soil moisture have been captured and the retrieval results are promising.

V1 V2 BIAS (m 3 /m 3 ) MAE (m 3 /m 3 ) RMSE (m 3 /m 3 ) R (-)
Since there is no discernible difference between the effective roughness parameter selected for each combination of vegetation parameters, the soil moisture retrieval accuracies show no noticeable sensitivity to the choice of particular vegetation parameters. Nevertheless, when the LAI and FPAR are selected as V 1 , the results of the soil moisture retrieval are slightly better. When the NDVI and EVI are selected as V 1 , the results are slightly poorer. The other combinations of vegetation parameters for the soil moisture retrieval are atypical. From the results of the model calibration and the soil moisture retrieval, it can be concluded that the LAI is suitable for characterizing the scattering properties of vegetation. This finding is in agreement with the results of previous works [14,27]. It is difficult to give a clear answer in which vegetation parameters can effectively characterize the attenuation properties of vegetation. Based on the foregoing discussion, it is concluded that optimum vegetation parameters are needed for the high accuracy retrieval of soil moisture.

Error Analysis
From Figures 5 and 6, it can be seen that the Sentinel-1A observations on the dates of 5 November 2017 and 11 December 2017 are extremely high compared with other observations. At the same time, it is found that the retrieved soil moisture is generally underestimated after 15 July 2018. These deviations for the model calibration and soil moisture retrieval can be analyzed from the following three aspects.
First, the biggest uncertainty may come from the effective roughness parameters. It is well known that this method can be used for eliminating the effect of soil roughness and achieved good performance in numerous studies, such as in Bai et al. [5,25], who used this method to parameterize the roughness parameters and retrieve the soil moisture. However, the validity of the effective roughness parameter was based on the assumption that the roughness parameters of the soil surface were unchanged during the study period. The study areas in the works of Bai et al. [5,25] are located in the prairie areas, which is less subject to human interference. Moreover, our study area is different. It belongs to the farming land, which is inevitably subject to human interference. Field works, such as sowing, irrigation, cultivation, reward, and others, may change the surface roughness. This may be the biggest problem faced by soil moisture retrieval in crop-covered areas.
Second, partial uncertainty may come from satellite data. On the dates of 5 November 2017 and 11 December 2017, the backscattering coefficients observed by Sentinel-1 are extremely high compared to the values on other dates. These cases are very difficult to simulate based on the coupled model, which only takes into account the effect of soil moisture, vegetation, and surface roughness. It is guessed that the sharp increase in the backscattering coefficient may be caused by other factors, such as field work or precipitation. Furthermore, the MODIS products may also bring some uncertainty, which may cause the bias for the estimation of the vegetation-scattering contribution. This probably could be another explanation for the inconsistency between the model simulations and the satellite observations.
Finally, the possible uncertainty may be caused by the difference of sensing depth between the SAR and in-situ measurements. The SAR signal can perceive the soil information below the surface. However, the perceived depth may change with the growth of vegetation. It can be said that the soil depth observed by the Sentinel-1A SAR is not fixed but changes with vegetation growth. By contrast, the in-situ soil moisture measurement is fixed at the depth of 0-5 cm. The inconsistency of sensing depth may explain the underestimation of soil moisture.
In summary, satellite data, soil-scattering models, vegetation-scattering models, and retrieval methods have to be carefully selected to achieve a high-precision soil moisture retrieval, particularly for developing operational retrieval methods.

Conclusions
In this paper, the WCM and AIEM were coupled to retrieve the soil moisture from the time-series Sentinel-1A and MODIS data in agricultural fields. Four different vegetation parameters were combined to characterize the scattering and attenuation properties of vegetation, including the LAI, FPAR, NDVI, and the EVI. The roughness parameters needed by the AIEM are parameterized by the effective roughness parameters. The soil moisture is retrieved based on the calibrated model using a LUT method.
The maximum and minimum values of RMSE are 1.401 dB and 0.996 dB, and the R ranges from 0.58 to 0.81. The statistical metrics indicate that the simulated backscattering coefficient is in good agreement with the Sentinel-1A observations. The calibration results are relatively better when LAI is used to parameterize the scattering properties of vegetation. The RMSE between the retrieved soil moisture and the in-situ measurements ranges from 0.069 to 0.092 m 3 /m 3 , and the R ranges from 0.74 to 0.81. These results indicate that the trend of retrieved soil moisture is consistent with the in-situ soil moisture, and the seasonal dynamics of soil moisture was captured. The coupled model calibration and soil moisture retrieval are performed with different combinations of the LAI, FPAR, NDVI, and the EVI as the vegetation descriptors. The LAI is recommended to characterize the scattering properties of vegetation. The vegetation descriptor selected for describing the attenuation properties of vegetation has to be further investigated. It is reasonable to suggest that the scattering and attenuation properties of vegetation are parameterized by two different vegetation parameters.
The results in this study provide a preliminary assessment for estimating the soil moisture from the Sentinel-1A and the MODIS data in agricultural fields. The inconsistencies still exist between the Sentinel-1A observations and the model simulation, and between the in-situ measurement and retrieved soil moisture. The reasons may ascribe to three aspects, including the effective roughness parameters, the satellite data, and the sensing depth. Therefore, more work is needed to perfect the method used in this paper. For this purpose, the applicability of the effective roughness method has to be studied over different agricultural fields. Furthermore, the uncertainty due to the satellite data has to be addressed, and the sensing depth needs to be considered in the coupled model. Further research shall be conducted to address more details on the combination of Sentinel-1A and MODIS for soil moisture retrieval.

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