Soil Moisture Retrievals Using Multi ‐ Temporal Sentinel ‐ 1 Data over Nagqu Region of Tibetan Plateau

: This paper presents an approach for retrieval of soil moisture in Nagqu region of Tibetan Plateau using VV ‐ polarized Sentinel ‐ 1 SAR and MODIS optical data, by coupling the semi ‐ empirical Oh ‐ 2004 model and the Water Cloud Model (WCM). The Oh model is first used to estimate the surface roughness parameter based on the hypothesis that the roughness is invariant among SAR acquisitions. Afterward, the vegetation water content (VWC) in the WCM is calculated from the daily MODIS NDVI data obtained by temporal interpolation. To improve the performance of the model, the parameters A, B, and α of the WCM are analyzed and optimized using randomly selected half of the sampled dataset. Then, the soil moisture is retrieved by minimizing a cost function be ‐ tween the simulated and measured backscattering coefficients. The comparison of the retrieved soil moisture with the ground measurements shows the determination coefficient R 2 and the Root Mean Square Error (RMSE) are 0.46 and 0.08 m 3 /m 3 , respectively. These results demonstrate the capability and reliability of Sentinel ‐ 1 SAR data for estimating the soil moisture over the Tibetan Plateau.


Introduction
Soil moisture is an important variable in the earth ecosystem and a driving parameter for hydrologic, climate, and ecological models [1][2][3][4]. The spatial distribution and temporal evolution of soil moisture vary significantly [5,6], so it is of significance to study the spatiotemporal dynamic changes of soil moisture in a wide range. The traditional soil moisture information collected via the ground point measurements provides limited insight into the spatial distribution patterns of soil moisture at a short time interval.
With penetration capability through the vegetation canopy and atmospheric layers coupled with the advantage of all-weather observations [7], microwave remote sensing technology has achieved success in monitoring soil moisture over large-area [8]. Passive microwave remote sensing data on soil moisture retrieval has a long history and is relatively mature, but can be limited in terms of providing high-resolution detail. The synthetic aperture radar (SAR) in active microwave remote sensing has high spatial resolution and can capture the subtle spatial information features of the soil surface [9]. The Sentinel-1 mission composed of a constellation of two satellites provides high-resolution synthetic aperture radar (SAR) data. The double satellites increase their temporal resolution to six days [10]; thus, Sentinel-1 has received extensive attention in soil moisture retrieval [11]. Many studies have shown the appreciable capability of C-band Sentinel-1 SAR data in estimating soil moisture and vegetation parameters. For example, Vijay et al. [12] retrieved soil moisture through a modified water cloud model (WCM) and an improved inversion algorithm. Kumar et al. [9] estimated winter wheat crop growth parameters using time series of Sentinel-1A SAR data.
The SAR backscattering coefficient is sensitive to soil dielectric constant, which is directly related to soil moisture content, laying a solid physical foundation for microwave remote sensing retrieval of soil moisture [13]. However, in addition to being controlled by soil dielectric properties, the backscattered signal is strongly influenced by vegetation and surface roughness [14]. The quantification of the respective contributions from vegetation and surface roughness affects the accuracy of soil moisture retrieval. Nevertheless, the interaction between ground object and radar scattering signal is complex, and the scattering mechanisms are sophisticated due to the overlaps among surface, volume, and dihedral scatterings. This leads to relatively large uncertainty in the subsequent soil moisture retrieval [15,16]. The developments of the soil moisture retrieval methods based on SAR observation have been a hotspot and difficult area.
Vegetation scattering models describe the microwave scattering process within the vegetation canopy. So the coupling of the vegetation scattering model and bare soil scattering model may eliminate the influence of surface roughness and vegetation cover to improve the accuracy of retrieval. For bare soil, scattering models mainly include empirical/semi-empirical models, such as the Oh-2004 model [17] and Dubois model [18], and theoretical models such as the Integral equation model (IEM) and the Advanced IEM (AIEM) [19][20][21]. These models were evaluated on the TerraSAR-X data over the bare soils, and the Oh model provided the best results [22]. The Oh-2004 model simulates the radar backscattering coefficient in term of soil moisture to avoid the intermediate processing of converting the dielectric constant from soil moisture. In this model, only the root mean square height (RMSH) was used to characterize the surface roughness, which reduces the number of unknown parameters and the uncertainty in the inversion process. Compared to theoretical models such as IEM and AIEM, this model has a simple formula and a higher efficiency in retrieval process. For vegetation cover areas, the commonly used scattering models mainly include the WCM [21] and the MIMICS model [23]. The WCM supposes the backscattering coefficient is composed of the contributions from the vegetation canopy and soil surface. The MIMICS takes into account additional backscattering components, such as plant-ground and ground-plant interactions, as well as ground-plant-ground scattering contributions. However, for the low grass condition over the Tibetan Plateau, the WCM with a smaller number of unknown parameters is more suitable to retrieve the soil moisture than the MIMICS [22,[24][25][26]. Consequently, the WCM and Oh models are often coupled to estimate soil moisture in the vegetated area. The WCM describes the vegetation characteristics using the Vegetation Water Content (VWC) which can be estimated from optical remote sensing data [27].
In [26], different radiative transfer models were evaluated over the wheat fields using the Sentinel-1 data. The LAI along with the dynamic empirical parameters seems better to capture the variation of the backscattering coefficients than the static empirical parameters. Furthermore, the WCM was calibrated using C-band Sentinel-1 data over the agricultural and grass lands. The NDVI extracted from the Sentinel-2 optical data was used as a vegetation descriptor in the WCM, and the IEM was applied to interpret the surface scattering component. Based on the radiative transfer models, a variety of SAR soil moisture retrieval algorithms were formed, such as change detection, optimization cost function, linear/nonlinear interpolation, lookup table, Bayesian posterior estimation, and artificial neural network [28]. The change detection was used to monitor surface soil moisture over the heterogeneous areas [29]. The cost function algorithm was optimized to estimate soil moisture [30,31]. The statistical relationships between radar backscattering coefficient and surface soil moisture were approximated as linear [32], exponential [33], or polynomial [34]. The process of lookup table retrieval is similar to that of least squares interpolation retrieval [35,36]. The Bayesian posterior estimation and artificial neural network were also evaluated to retrieve the soil moisture from remote sensing data [37]. Among them, the cost function optimization method which estimates the surface parameters by iteratively minimizing system errors, is considered with high universality and reliability [38,39]. The optimal estimates are constrained by physical significance and global optimization status [10,40].
Considered as the "Third Pole of the Earth", the Tibetan Plateau is the highest one in the world. Soil moisture is one of the most sensitive factors affecting precipitation patterns and water cycles in this region due to its unique topographical features, so an adequate understanding of the soil moisture over the Tibetan Plateau is necessary. Considerable pioneer soil moisture retrieval works have been done over the Tibetan Plateau, either through radiometer [41], or radar signals [42]. However, it is still necessary to evaluate the impact of different vegetation and soil parameterization on the accuracy of soil moisture retrievals. Compared to [42] which used AIEM for surface component and LAI for characterization of the vegetation, we used the Oh model for surface component and applied an empirical relationship between the VWC and NDVI. The purpose of our study is to analyze the potential of the combination between Sentinel-1 SAR and MODIS imageries for soil moisture retrieval on the Tibetan Plateau. We couple the Oh-2004 [17] model into the WCM [15,43] to simulate SAR backscatter, then construct a cost function to minimize the difference between the model simulated and Sentinel-1 observed backscatters and to search for the optimal solution for soil moisture retrieval.
The novelty of this paper relies in the following aspects: 1) We utilized a modified volume scattering component in which a radar-shadow coefficient is introduced for the WCM model. 2) The Sentinel-1 SAR data were combined with MODIS imageries to retrieve soil moisture over the Tibetan Plateau. Since the soil moisture is one of the most sensitive factors affecting precipitation patterns and water cycles in this region, an adequate understanding of the soil moisture is necessary considering the unique topographical features. 3) A transfer function between VWC and NDVI [44] was evaluated to estimate the VWC for the study area. This may benefit the selection of transfer function to obtain the VWC parameter which is essential to parameterize the radar backscattering coefficient. Therefore, the current study deepens the understanding on the optical-microwave jointing estimation of soil moisture over the Tibetan Plateau. The outline of this paper was given as follows: Section 2 presents the details of the study area, in situ moisture, and Sentinel-1 SAR and Terra MODIS data. Section 3 describes the details of methodology, including vegetation and soil backscatter decomposition, model calibration, and soil moisture retrieval, followed by Section 4 for result analysis. Section 5 discusses the potential uncertainty in soil moisture retrieval, followed by the conclusions in Section 6.

Study Site
The study site is located in Nagqu, northern Tibet, in the hinterland of the Tibetan Plateau ( Figure 1). Nagqu plays a key role in terrestrial energy balance. The selected site has an average altitude of above 4500 m. Despite this altitude, the terrain is relatively flat, with rolling hills varying only a few dozen meters in elevation. The vegetation cover is mainly high-situated pasture. Its climate is characterized by extreme cold and hypoxia, dry and windy. The annual average temperature is −0.9 ℃ to −3.3 ℃, the annual relative humidity is 48-51% and the annual precipitation is 380 mm [45]. From November to the following year in March, it is a dry and windy period. From April to October, the climate is relatively warm, and it is the golden growing season of grassland. For this mountainous area, the soil moisture dynamics have a large spatial variability [46]. Therefore, this area is selected for our study to develop and validate a soil moisture retrieval algorithm using the Sentinel-1 and MODIS data.

In situ Soil Moisture
The in situ soil moisture data are obtained from the Nagqu soil moisture and temperature network in the central Tibetan Plateau. The network contains 56 soil moisture stations ( Figure 1). This network (ISMN: https://ismn.geo.tuwien.ac.at/en/) provides soil moisture and temperature with a time interval of 30 min at four different depths (0~5 cm, 10 cm, 20 cm, and 40 cm) from 2011 to 2016 [47,48]. The details of the Nagqu network can be found in Yang et al. [49]. Considering the potential penetration depth at C-band, this study selects the in situ soil moisture measurements at the top soil layer (0~5 cm) for the model evaluation and validation. However, there is a serious data loss, as part of the Sentinel-1 images can not completely cover the 56 stations. By appropriate screening, 43 stations are selected for further research. In addition, we have made quality control of the in situ soil moisture data according to ISMN Quality Flags. Figure 2 shows the temporal evolution of the mean values of the measured soil moisture and temperature over the 43 stations from April to November 2015, which indicate a seasonal variability.

Sentinel-1 SAR Data
Sentinel-1 is the European Space Agency's (ESA) Copernicus Earth Observation satellite. It consists of two satellites carrying C-band synthetic aperture radar with a central frequency of 5.405 GHz. Since the revisit period of a single satellite is 12 days, and that of A/B twin satellites reduces to 6 days. The SAR has a variety of imaging modes: Striped Map (SM), Interferometric Wide Swath (IW), and Extra Wide Swath (EM). We adopt the Ground Range Detected (GRD) product of IW mode with VV and VH dual-polarized modes and spatial resolution 5 m × 20 m. Many studies have shown that the cross-polarized VH channel is more sensitive to vegetation change, while co-polarized VV is sensitive to both vegetation and soil moisture, and thus the latter has greater potential in soil moisture retrieval [50,51]. Therefore, only VV polarized backscatter signal is selected in this study.
The Sentinel-1 Level 1 products including 38 scenes from April to November 2015, were downloaded from the Sentinel Data Hub website (https://scihub.copernicus.eu/). The full preprocessing chain comprises: 1) radiometric calibration to reduce radiation deviation; 2) refined Lee speckle filter with 5 × 5 windows to reduce the speckle noise; and 3) geometric corrections (Range-Doppler terrain correction) to correct the distortion of ground objects in radar images with 30 m SRTM DEM. The preprocessing is conducted using the SNAP software.
Given the spatial heterogeneity of soil moisture, buffer zones with a 200 m radius are established for each of the 43 ground stations in the study area. Then, the backscattering coefficients at VV polarized and the associated local incidence angle of each station are extracted and averaged for the model development.

MODIS Data
MODIS has two sensors that operate on the Terra and Aqua spacecraft, respectively. Both the Terra NDVI (MOD13Q1) and Aqua NDVI (MYD13Q1) are level-3 grid data in the Sinusoidal projection mode at 250 m spatial resolution and the 16-day temporal resolution derived from daily surface reflectance by selecting the high-quality pixels (https://search.earthdata.nasa.gov/search) [52]. Since the Terra NDVI data are used to characterize the vegetation dynamics in the current study, we perform a temporal interpolation on the Terra NDVI products in order to match them to the Sentinel-1 data acquisition time. Before the analysis, the two datasets need to be preprocessed via the MODIS Re-Projection Tool (MRT), including re-projection, file format conversion and spatial composition. Then the NDVI data for the soil moisture stations are extracted using the same buffer established for the Sentinel-1 SAR data. To verify the accuracy of the interpolation result, the daily mean NDVI derived from the interpolation of Terra NDVI data are compared with Aqua NDVI data in the corresponding period. Figure 3 shows a high consistency between the interpolated NDVI values and the Aqua NDVI products over a representative station L02. The study area is covered by grass, and the temporal variation of the interpolated NDVI can well capture the phenological rhythms of the grass growths ( Figure 3). Consequently, the NDVI dataset obtained by interpolation will be used to describe the vegetation influences on the radar backscattering coefficient.

Methodology
In this paper, the Oh-2004 model and the WCM are coupled to simulate the backscattering coefficients at VV polarization for Sentinel-1 SAR data. MODIS data are used to compute the vegetation contribution in the SAR backscattering coefficients. By optimizing a cost function constructed by the model simulations and SAR observations, the soil moisture was retrieved. The implementation of the coupled model calibration and soil moisture retrieval is illustrated in Figure 4, and detailed development of the coupled model is described in the following sections.

Vegetation and Soil Backscatter Decomposition
Since the vegetation disturbance is an important factor affecting the accuracy of soil moisture retrieval, we apply the WCM to correct the vegetation effect. This model is an empirical model proposed by Attema and Ulaby for agricultural crops, and has a wide range of applications [43]. The model makes the following hypothesis: (1) the vegetation is represented as a homogeneous horizontal cloud of identical water spheres, uniformly distributed throughout the space defined by the soil surface and the vegetation height. (2) the double bounce power is relatively low due to the shallower penetration depth at Cband, so it is theoretically appropriate to ignore the double-bounce component. According to [53], the dihedral scattering power at C-band is less than 10% during the entire phenological stages. In [54], the dihedral scattering power was also found low at C-band, and can be neglected in the WCM. Therefore, the total backscatter signal received by radar can be expressed as: (1) where is the total SAR backscattering coefficient, and are the backscatters from vegetation and underlying soil, respectively.

Surface Component
The backscattering from the bare soil surface is simulated by the semi-empirical Oh-2004 model. In this model, the surface roughness is described using surface Root Mean Square height (RMSH), without considering the autocorrelation length. From the aspect of soil moisture retrieval, reducing the number of unknowns reduces the uncertainty of the retrieval [55]. Consequently, the VV-polarized backscattering coefficient of the bare soil surface can be expressed as: where SM is the soil moisture; θ is incidence angle; k and s are the wavenumber and RMSH (cm), respectively. The validity of Oh model was tested for 0.04 < SM < 0.29 m 3 /m 3 , 0.13 < ks < 6.98 and 10°< θ < 70° [25,26].

Volume Component
In a radar beam, two or more canopies usually overlapped, leading to overestimation of the vegetation backscatter due to the canopy layover and radar shallow. Thus, we utilize a modified volume scattering component in which a radar-shadow coefficient is introduced to account for the effect of vegetation layover [15]: where θ is the incidence angle, VWC (kg/m 2 ) is the vegetation water content as computed by NDVI; is the two-way vegetation transmissivity, which can be obtained by parameterization of VWC. A and B are the empirical parameters of the model depending on the canopy type, α is the radar-shadow coefficient depending on vegetation type and land use.
The NDVI was used to estimate the VWC, which comprised of the foliage water content derived from the NDVI, and the stem water content derived from annual NDVI extremes along with a coefficient of stem factor [44]: where the NDVI is obtained from the MODIS-Terra MOD13Q1 product, NDVImax and NDVImin represent the annual maximum and minimum values of the NDVI for a given location. The stem factor as a product of average vegetation height and the ratio of sapwood area to leaf area, represents an estimate of the peak amount of water residing in the stems, and depends on vegetation types. The default value of stem factor is set at 1.5 when the vegetation type is grassland [56]. However, considering the low grass covering the study area, this default value may be too high. Therefore, in this study, we adjusted the factor to 0.3 [57]. By adjusting the stem factor, the Eq. (5) match well to the ground measurements, and maintain the reasonable temporal behavior of VWC in term of vegetation growth phenology.

Model Calibration and Soil Moisture Retrieval
The parameters (A、B, and α) in Eq. (3) were determined based on randomly selected 50% of the data. The remaining 50% data were used to retrieve the soil moisture, followed by the validation against the ground measurements. To obtain globally optimal solutions for different parameters, we minimized a cost function J (dB) was constructed using the simulated backscattering coefficients( )and radar observations( ):

Surface Roughness
The surface roughness was regarded stable over a short time [30], and optimized using Oh-2004 model applied to the SAR data acquired under bare and unfrozen soil conditions. We choose the Sentinel-1 image and the ground measured soil moisture collected on April 1, 2015 for which the soils are not covered by vegetation as input data to estimate the roughness [31]. In the process, roughness ks is set as a free parameter ranging between 0.1 and 3.0 to iteratively match the model simulations to the Sentinel-1 SAR backscattering coefficient.

Parameters A、B and α
To examine the influence of the parameters A、B and α of the WCM on soil moisture retrieval, we conducted a synthetic simulation experiment to observe the sensitivity of the SAR backscattering coefficient to the variations of the involved vegetation parameters. The experiment consists of three individual tests. For each test, only one parameter varies within its range, the other two parameters are fixed. The ranges of the parameters variation are set according to Table 1 from the findings in Bindlish and Barros [15]. In all three tests, the value of VWC ranges from 0 to 1.0 with an interval of 0.2 kg/m 2 . Based on the simulation, the cost function is used again to calibrate the parameters.

Retrieval of Soil Moisture
In the retrieval process, the calibrated parameters A、B and α and the estimated surface roughness are brought into Eq. (6), and the soil moisture is set as a free parameter to iteratively match the simulated backscattering coefficients to the Sentinel-1 SAR observations. Finally, the retrieved soil moisture is validated against the ground measurements.

Results
This section analyzes the estimated surface roughness, parameters of the WCM and soil moisture by applying the proposed approach to the Sentinel-1 SAR and MODIS data.

Sensitivity of Sentinel-1 Backscatters to Surface Parameters
To explore the response of SAR backscattering coefficients to the variations of soil moisture and NDVI, a simple regression analysis is first conducted. As shown in Figure  5, there is a relatively weak positive correlation between the backscattering coefficient and soil moisture (Figure 5a). The deviation of the points from the fitted line is attributed to the differences in vegetation condition, surface roughness and local incidence angle at different stations. A similar weak sensitivity of the C-band SAR backscattering coefficient to measured soil moisture is observed in many previous studies due to the joint effects of soil and vegetation parameters [29,32,58]. In Figure 5b, the influence of the vegetation on the backscattering power was observed. In addition, with an increasing NDVI, an increas-ing trend is observed in soil moisture (Figure 5c), confirming close coupling between vegetation and soil moisture. Thus, in order to retrieve soil moisture from C-band Sentinel-1 SAR data, it is necessary to better parameterize the vegetation contribution.

Estimated Surface Roughness
In this process, we applied Eq. (2) to the VV polarized backscattering and ground measurements, resulting in surface roughness values for 43 soil moisture ground stations (Figure 6a). In addition, the mean soil temperature (Figure 6b) and NDVI (Figure 6c Based on the optimized surface roughness parameters and the measured soil moisture, the backscattering coefficient is simulated. In Figure 7, a relatively high correlation between the simulated backscattering and Sentinel-1 SAR data at VV polarized channel was found. It indicates the feasibility of optimizing the assumed stable surface roughness parameters prior to retrieving the soil moisture.

Optimized Vegetation Parameters A, B, and α
Before the optimization, the sensitivity of simulated SAR backscattering coefficients to the variation of the vegetation parameters (A, B, and α) is analyzed for six different VWC levels. Figure 8 shows the response of simulated backscatters to the three vegetation parameters. As can be seen in Figure 8a,c, almost no changes are observed with increasing values of parameters A and α, but with different VWC levels, the backscatters show a certain difference. In contrast, in Figure 8b, the backscattering coefficients decrease with increasing values of parameter B and the variation amplitude increases with the increase of VWC. It indicates that parameters A and α are not sensitive to soil moisture retrieval, so parameter B should be mainly calibrated in the subsequent process. Similar observations for the vegetation parameters were found in previous studies [59]. Thus, the values of parameters A and α, are set to the "all land uses" values in Table 1. In Figure 9, the parameter B iterates between 0.03 to 0.14 with 100 intervals, and the lowest J value was obtained when B is equal to 0.05.

Retrieved Soil Moisture
After the calibration of the coupled Oh-2004 and WCM model, the remaining 50% Sentinel-1 data are used for soil moisture retrieval. As shown in Figure 10, the retrieved soil moisture shows consistency with observations over the ground stations. However, two outliers circled in red are found, with a large overestimation. The NDVI values for the two outliers are 0.25 and 0.33, which are lower than the averaged 0.53 for all the stations. On the other hand, their respective ks values are 1.0 and 0.96, which are significantly higher than the averaged 0.3 for all the stations, but still in the validity range of the Oh model. Nevertheless, the combined effects of higher roughness and lower vegetation condition result in a similar level of radar backscattering coefficient as the mean value across all the stations. Except for these two outliers, a determination correlation R 2 = 0.46 and root-mean-square error RMSE = 0.08 m 3 /m 3 were obtained. This indicated that the combination between Sentinel-1 SAR and MODIS optical remote sensing data is feasible to retrieve soil moisture over the Tibetan plateau. Furthermore, Figure 11 shows the temporal profiles of the retrieved and measured soil moisture over four stations for a more detailed analysis. To a certain extent that the soil moisture is closely related to the MODIS NDVI. This result is consistent with our previous sensitivity analysis ( Figure 5). The more water exists in the soil, the more vegetation grows which increases MODIS NDVI values. Vegetation growth and enhanced density also help to maintain the water particles in the soil layer. This indicates that the interaction between soil moisture and NDVI is a bidirectional process. The temporal variation of trends of the retrieved soil moisture is consistent with the ground measurement. How-ever, the deviation was observed over station M05 probably due to the variation of topography. Over the station M05, the incidence angle is above 50°, higher than the average value 38.6° across all the stations. As the incidence angle increases, the vegetation contribution dominates the backscattering power, while the soil contribution decreases leading thus to an underestimation of soil moisture.

Discussion
In our study, the Oh-2004 and the WCM are coupled to estimate the soil moisture from VV-polarized Sentinel-1 SAR and MODIS data in Nagqu. We propose to decompose the backscattering coefficients into vegetation canopy volume scattering and soil surface scattering mechanisms, and obtain soil moisture from the latter soil-related components. As can be seen from the retrieval results with a determination coefficient R 2 = 0.46 and the RMSE = 0.08 m 3 /m 3 , this accuracy is within a similar level as the literature. For instance, Bai et al. [42] estimated soil moisture over Tibetan Plateau by coupling the WCM and the AIEM. The retrieval results are with R 2 equal to 0.36 and 0.67, and RMSE of 0.073 m 3 /m 3 and 0.055 m 3 /m 3 for the ascending and descending data, respectively. Compared to [60] which applied the IEM model to multi-angular ASAR data over Tibetan Plateau, our method with the vegetation correction obtain better results. Chai et al. [61] used RADAR-SAT-2 imagery to model and map soil moisture of Plateau Pasture, with R 2 ranged from 0.59 to 0.7 and RMSE ranged from 0.043 to 0.048 m 3 /m 3 .
Nevertheless, the performance of the approach is subject to incident angle, vegetation and surface roughness parameters. The radar incidence angle determines the direction and transmission path of the radar beam. With the increase of incident angle, the distance of microwave propagation, and the amount of vegetation penetrated gradually increased [62]. Therefore, the soil contribution decreases while vegetation contribution predominates at a high incidence angle. There is an approximately exponential dependence between surface roughness and backscattering coefficient, and this dependence is more significant with the increase of incident angle [63]. Thus, the high incident angle makes it challenging to monitor soil moisture. This phenomenon was confirmed over the station M05 where an underestimation of soil moisture was found due to the high incidence angle.
The performance is also affected by assumptions made for the development of the algorithm. In our current study, the surface roughness which is derived by hypothetical bare and unfrozen soil conditions is assumed to be temporally stable during the Sentinel-1 acquisitions. However, the roughness parameter may be modified by rainfall or human factors [64]. Previous researches have indicated that stable roughness parameters are unable to characterize the changeability of the soil surface [65]. It is better to update the effective roughness parameters for each acquisition. In our case with a single polarization radar configuration, the elimination of the roughness effect is challenging, since the roughness significantly varies with local incidence angle, microwave frequency and polarization. Thus, multi-dimensional SAR configurations are necessary to refine the parameterization of the soil surface roughness on the soil moisture retrievals [66,67]. Furthermore, in the equation system, the estimation accuracy of effective roughness parameters has an important impact on the calibration of the WCM parameters A, B, and α, and thus affect the subsequent retrieval of soil moisture.
In addition to the impact of surface roughness, parameters A, B, and α in the WCM also have an impact on soil moisture retrieval. Since the coefficients of the WCM model are empirically fitted, the representability of the training data impacts the performance of the model. For test sites with different vegetation characteristics, the empirical coefficients need to be tuned. In our study, to reduce the number of known parameters, we first carried out a simulation analysis, and concluded that the VV-polarized backscattering coefficient had the highest sensitivity to the change of parameter B compared to the two other parameters. Based on this conclusion, the subsequent calibration was focused on parameter B, while A and α are directly taken from the previous studies. Furthermore, the Vegetation Water Content was derived from the NDVI. Nevertheless, with the increase of NDVI value, the sensitivity of radar signal to soil moisture decreased, increasing the uncertainty in extracting the soil moisture information [68]. In the following work, it is necessary to refine the parameterization of vegetation parameters in the WCM by extensive vegetation ground measurements.
Furthermore, in our current experiment, we randomly selected the training data to calibrate the parameters in WCM, which were then used in the further retrieval of soil moisture. We repeated the test for 10 times and found that the calibration with different training datasets may result in slightly different values of calibrating parameters. Although the results of soil moisture retrieval are affected by the selection of training datasets, this effect is relatively small over the same study area. The uncertainty caused by selecting different datasets is partially compensated by the associated parameter calibration.
Finally, the optimization algorithm used in the calibration step impacts the accuracy of soil moisture retrieval. Although we used a cost function to determine the vegetation parameters in WCM, other sophisticating optimization approaches such as the genetic algorithm may be considered as alternative approaches. In addition, the iterative loop setting and the convergence conditions also influence the optimization of the WCM parameters and the accuracy of subsequent soil moisture retrieval. In the current paper, as the number of iterations reaches around 100, the RMSE of the retrieval result tends to converge stably and the calculation time is relatively small.

Conclusion
In this paper, we combined the Oh-2004 model with the WCM to conduct soil moisture retrieval in Nagqu region of Tibetan Plateau using the sentinel-1 SAR and MODIS optical data. The main conclusions are summarized as follows: We carry out a linear interpolation on the MODIS NDVI 16-days composites to obtain the daily NDVI, then substitute it into an empirical model (Eq. (5)) to calculate the corresponding vegetation water content. The result demonstrates the temporal interpolation for limited measured vegetation growth parameters is feasible to capture the vegetation dynamics. The retrieved and measured soil moisture were close related to the NDVI, indicating the coupling mechanism of the soil-vegetation system.
As for roughness, we assume that it is stable over a relatively short time, and estimate it by inverting the backscattering coefficient of a SAR image acquired at bare soil conditions. Since the land cover in the study area is mainly grassland, and the variation of land surface attributes is small. Therefore, this assumption is reasonable to reduce the surface roughness effect on soil moisture retrieval.
The observed inconsistencies between the in situ measurement and retrieved soil moisture mainly attribute to the representativity of model calibration and the quality of vegetation parameters. Additional investigation of these factors is likely to provide a deeper insight into soil moisture retrieval. Although the Sentinel-1 SAR-based soil moisture retrieval remains challenging, the precise scale of SAR data is significantly important for agricultural and water management applications, as well as for the further development of land surface and hydrologic process models.