Next Article in Journal
Mapping the Distribution of Shallow Groundwater Occurrences Using Remote Sensing-Based Statistical Modeling over Southwest Saudi Arabia
Next Article in Special Issue
Using Apparent Electrical Conductivity as Indicator for Investigating Potential Spatial Variation of Soil Salinity across Seven Oases along Tarim River in Southern Xinjiang, China
Previous Article in Journal
Crop Monitoring Using Satellite/UAV Data Fusion and Machine Learning
Previous Article in Special Issue
Limitations in Validating Derived Soil Water Content from Thermal/Optical Measurements Using the Simplified Triangle Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Capability of Integrating Optical and Microwave Data for Detecting Soil Moisture in an Oasis Region

1
College of Resources and Environment Sciences, Xinjiang University, Urumqi 830046, China
2
Key Laboratory of Oasis Ecology under Ministry of Education, Xinjiang University, Urumqi 830046, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2020, 12(9), 1358; https://doi.org/10.3390/rs12091358
Submission received: 19 February 2020 / Revised: 22 April 2020 / Accepted: 24 April 2020 / Published: 25 April 2020
(This article belongs to the Special Issue Monitoring Soil Degradation by Remote Sensing)

Abstract

:
In the earth ecosystem, surface soil moisture is an important factor in the process of energy exchange between land and atmosphere, which has a strong control effect on land surface evapotranspiration, water migration, and carbon cycle. Soil moisture is particularly important in an oasis region because of its fragile ecological environment. Accordingly, a soil moisture retrieval model was conducted based on Dubois model and ratio model. Based on the Dubois model, the in situ soil roughness was used to simulate the backscattering coefficient of bare soil, and the empirical relationship was established with the measured soil moisture. The ratio model was used to eliminate the backscattering contribution of vegetation, in which three vegetation indices were used to characterize vegetation growth. The results were as follows: (1) the Dubois model was used to calibrate the unknown parameters of the ratio model and verified the feasibility of the ratio model to simulate the backscattering coefficient. (2) All three vegetation indices (Normalized Difference Vegetation Index (NDVI), Vegetation Water Content (VWC), and Enhanced Vegetation Index (EVI)) can represent the scattering characteristics of vegetation in an oasis region, but the VWC vegetation index is more suitable than the others. (3) Based on the Dubois model and ratio model, the soil moisture retrieval model was conducted, and the in situ soil moisture was used to analyze the accuracy of the simulated soil moisture, which found that the soil moisture retrieval accuracy is the highest under VWC vegetation index, and the coefficient of determination is 0.76. The results show that the soil moisture retrieval model conducted on the Dubois model and ratio model is feasible.

1. Introduction

Oases are important parts of the ecosystem that constitute arid and semi-arid regions; they are the “source of life” throughout the arid and semi-arid regions. China’s oases are mainly distributed in the five northwestern provinces, especially in the Xinjiang Uygur Autonomous Region (hereinafter referred to as “Xinjiang”). The Xinjiang oases covers 8% of the total area, but bears more than 90% of the cultivated land, population, and Gross Domestic Product (GDP) of Xinjiang. Meanwhile, the inner ecological environments of oases are fragile; once they are affected or destroyed by human or natural causes, it is difficult for them to restore to normal.
Soil moisture is an important parameter in meteorological, agro-environmental, and hydrological studies [1]. In meteorological research, accurate detection of surface soil moisture is important and can be applied in a variety of hydrological and climate models (e.g., rainfall mapping, drought monitoring, vegetation water requirement analysis, etc.). In agricultural research, the most important thing is to estimate the crop yields; in addition, the surface soil moisture can be introduced into the basin hydrological model as a priori knowledge (which can predict the deep soil moisture content in order to study vegetation water stress, etc.). Finally, soil moisture can help predict natural disasters, such as debris flows, landslides, etc.; meanwhile, surface soil moisture mapping plays a crucial role in improving local and even global climate predictions [2]. Moreover, soil moisture is an important indicator for monitoring land degradation and drought [3], an indispensable key factor in the land water cycle [4,5,6], and the basic guarantee for vegetation growth and development. It has guiding significance in estimating crop yield and drought monitoring [7,8].
In arid and semi-arid regions, the spatial and temporal soil moisture retrieval can be used to better understand the changes of soil moisture migration, to further understand the ecological effects of vegetation and, simultaneously, improve soil moisture holding capacity and vegetation water use efficiency (and, accordingly, play an important role in the sustainable development of an oasis ecology). Therefore, the importance of soil moisture for the sustainable development of an oasis is self-evident, and real-time monitoring of large-scale soil moisture is of great significance [9].
At present, high-precision retrieval of spatial and temporal distribution of soil moisture is an urgent problem to be solved [10,11], to which the hotspots and difficulties lies in the soil moisture retrieval under vegetation coverage. At present, the soil moisture monitoring method is based on the soil moisture monitoring site, which is time-consuming and labor-intensive, and the observation site is limited and unevenly distributed. On the other hand, with the rapid development of remote sensing technology, more and more people use remote sensing methods to monitor soil moisture. After several decades development of remote sensing technology, the data sources of soil moisture retrieval mainly include optical, microwave and thermal infrared [12]. Optical remote sensing has developed rapidly and has become a high time resolution, high spectral resolution, and high spatial resolution data [13,14], but it is susceptible to the external environment, such as atmosphere, clouds, fog, etc. It is susceptible to leaves and stems under vegetation-covered areas; with its strong penetrating ability, microwave remote sensing can realize all-day observation under various meteorological conditions, and most importantly, it is sensitive to soil moisture changes, making microwave remote sensing widely used in large-scale mapping of soil moisture in arid and semi-arid regions [15].
The microwave remote sensing has two kinds of observing systems, active and passive missions. The passive mission provides global soil moisture products with a high temporal resolution (2–3 days) and high retrieval accuracy [16], but with a poor spatial resolution. The active mission provide a high spatial resolution data, but can be influenced by the soil and radar’s system. Taking the study area into account, it was not fit for the passive remote sensing soil moisture derive, so the active mission was selected as the data source. Over the past few decades, active microwave remote sensing data have been successfully used for estimating soil moisture due to a finer spatial resolution [17]. Several soil moisture retrieval algorithms have been developed and tested for multiple Synthetic Aperture Radar (SAR) satellites operated at the L/C/X-bands [18], such as Advanced Land Observing Satellite-2 (ALOS-2) [16,19], Radarsat-2 [20,21], Advanced Synthetic Aperture Radar (ASAR) [22], Sentinel-1 [23,24,25,26,27], and Multi-temporal X-band SAR (TerraSAR-X) [28,29]; studies that have used these data to estimate bare soil moisture have achieved promising results. Based on this, several soil moisture models have been proposed to the bare land, including the statistical models (Oh model [30] and Dubois model [31]) and physical models (the integral equation model (IEM) [32] and advanced IEM (AIEM) [33,34,35,36,37]). Combined with the vegetation scattering model [38,39], which can be used to quantify vegetation attenuation of radar signals in radiative transfer function models [40], the bare soil scattering model can adapt to the vegetation areas. To eliminate the effect of vegetation, some studies [39,41,42,43,44] have attempted to use additional vegetation information provided by optical remote sensing, which has been widely used to derive information of vegetation properties. Furthermore, other studies [43,45,46,47] have suggested that the accuracy of soil moisture estimates was significantly improved when considering synergy between radar and optical data as compared to estimates from SAR data only.
The aim of this paper is to verify the feasibility of the soil moisture retrieval model based on bare surface scattering models and vegetation scattering models in an oasis region, in order to accurately derive soil moisture in the oasis of the arid region, which can better understand the response of the surface parameters and the canopy to backscatter. Based on this purpose, the Ugan-Kuqa River Delta Oasis in the Aksu region of Xinjiang was selected as our target area, with microwave data as the main and supplemented by optical data. We analyzed the relationship between backscattering coefficient and soil moisture based on microwave and optical data; therefore, a surface soil moisture retrieval model was conducted. The basic of this work is, coupled, the Dubois model and ratio model, using the surface roughness, soil moisture, and radar system parameters to simulate the backscattering coefficient of the bare soil to remove the backscatter contribution of vegetation and conduct a surface soil moisture retrieval model. This method can provide a theoretical support for the sustainable development of oases.

2. Materials and Methods

2.1. Overview of the Study Region

Currently, the ecological degradation problem in the lower reaches of the Tarim River Basin is serious; the upper reaches of the Tarim River Basin (Ugan-Kuqa River Delta Oasis) was selected as the study area. The main reason we selected this area is that the use of water resources is closely related to the ecological problems and this area an important food production base in Xinjiang. In this oasis region, vegetation coverage will affect the soil moisture retrieval accuracy based on optical remote sensing, and microwave remote sensing can provide a promising method of soil moisture derive.
The Ugan-Kuqa River Delta Oasis(latitude 41°00′ N–41°50′ N, longitude 82°10′ E–83°30′ E (hereinafter referred to as “this area”) is located on the northern edge of the Tarim Basin in Xinjiang. Luntai County lies to the east, Baicheng County lies to the west, the north bank of the Tarim River lies to the south, and the Qiulitage Mountain in the southern foothills of the Tianshan Mountains lies to the north (Figure 1). It is a typical fan-shaped plain oasis in China. This oasis has a temperate continental arid climate with abundant light and heat resources. The climate is very dry with a deficiency of rainfall. To be specific, the annual total radiation is 6.11 × 105 J cm2, the annual sunshine hours are approximately 2.85 × 103 h, the sunshine rate is 65%, the annual average temperature is 11.3 °C, the potential evaporation is 2356 mm, the average annual precipitation is 55.45 mm, and the drying index is 42 [27]. This oasis agriculture is mainly based on planting. Artificial vegetation mainly includes cotton, corn, and winter wheat, while natural vegetation is halophyte, which mainly includes Kalidium foliatum, Halocnemum strabilaceum, Halostachys capsica, Alhagi sparsifolia, Populus euphratica, Tamarix spp, Phragmites communis.
The land cover use of the study area is shown in Figure 2. As shown in Figure 2, there are six types of land cover, which include built-up, agriculture, natural vegetation, bare land, water, salinized. The artificial vegetation is mainly planted in the center of this oasis; meanwhile, the natural vegetation is mainly in the oasis-desert ecotone region where the soil salinization is serious.

2.2. Observation Data

2.2.1. Soil Moisture, Texture, and Bulk

The soil samples were collected at a depth of 10 cm in July 2017, which used a random clustered sampling technique to allocate 69 samples in the Ugan-Kuqa River Delta Oasis. Each sample should be representative of the soil surrounding a larger area (which had a minimum distance of 2 km to the next site to avoid autocorrelation). The samples were selected away from the trails and consisted of homogenous land use types (Table 1). Throughout the in situ soil sampling, the weather is clear and cloudless with no rainfall. The landscape in this region includes agriculture, natural vegetation, bare land, water, salinized. Soil data process were as follows:
  • The soil sample of each site was uniformly mixed and placed into an aluminum box, which was weighed and recorded, then transported back to the laboratory for drying until the soil was completely dehydrated, cooled to room temperature, and weighed in a cool place. The soil moisture was calculated. The range of measured soil moisture value was 0.00 to 47.59%, and the average was 24.17%.
  • The collected soil samples were transported to the laboratory, the plant and gravel impurities were removed after air drying, the soil was ground and sieved with a 0.5 mm aperture sieve; physical and chemical analyses were then conducted to measure soil texture and soil bulk density.

2.2.2. Soil Roughness

At present, there are two methods to obtain soil roughness: the probe method and laser scanning method. In this study, the probe method was used to obtain the surface roughness, which is a section panel with a centimeter grid (Figure 3a). The specific steps of soil roughness data process were divided into four steps:
  • Using a camera, the soil roughness was captured, ensuring that the samples board was perpendicular to the ground and that the camera was parallel to the soil surface.
  • The captured photos of soil roughness were vectorized in MATLAB (Figure 3b), and then a manual correction processing was performed to get a vector picture of the soil roughness (Figure 3c).
  • The discrete height value of each point was extracted by every 1 cm of the vectorized data, and data were saved in Excel for the next step.
  • The data of the discrete height value were analyzed according to the equations of root mean square height and related length; therefore, the correlation length and the root mean square height was obtained. The range of the correlation length value was 13.73 cm to 42.16 cm, and the average was 24.09 cm. The range of the root mean square height value was 0.32 cm to 3.69 cm, and the average was 1.68 cm.

2.3. Satellite Data

2.3.1. Sentinel-1 Data

Sentinel-1 is an imaging radar mission providing continuous all weather, day-and-night imagery at C-band, which launched on April 3, 2014. This SAR instrument carries two satellites (Sentinel-1A and Sentinel-1B), and may acquire data in four exclusive modes: Stripmap (SM), Interferometric Wide swath (IW), Extra Wide swath (EW), Wave (WV).
The microwave data source used in this study was single polarization, which included VV (vertical transmit and receive) and HH (horizontal transmit and receive) polarization modes. The data were acquired on July 15, 2017, coinciding with the soil sample acquisition times that were used. The data used in this study were the so-called interferometric wide swath mode at a spatial resolution of 20 by 20 m, and the azimuth width was 250 km. These data were downloaded from https://scihub.copernicus.eu/dhus/#/home.
The Sentinel-1 image was an ellipsoid geocoded Ground Range Detected (GRDH), preprocessed by the SARscape 5.2.1® module developed by SARmap in ENVI 5.3® software. Following these procedures, the processed SAR data with a 30 m spatial resolution were projected onto World Geodetic System (WGS) 1984 Universal Transverse Mercator (UTM) coordinates. The preprocessing of Sentinel-1 data are as follows:
  • Filtering and denoising processing (refined-Lee filtering, 3 pixels by 3 pixels window) was performed to eliminate speckle noise.
  • Geocoding was performed, using digital elevation maps for geometric fine correction.
  • Radiation calibration was conducted to obtain the backscattering coefficient in the multipolarization mode of the target region.

2.3.2. Landsat-8 OLI Data

Landsat-8, which was launched in February 2013 by National Aeronautics and Space Administration (NASA), carries two sensors, the OLI Land Imager and the Thermal Infrared Sensor (TIRS). For this study, OLI data from July 14, 2017, were downloaded from https://earthexplorer.usgs.gov/. The preprocessing of Landsat-8 OLI data includes radiation correction and atmospheric correction, which were conducted using the ENVI 5.3® software (Exelis Visual Information Solutions, Broomfield, CO, United States). Following these procedures, the processed OLI data with a 30 m spatial resolution were also projected onto WGS 1984 UTM coordinates.
The three indices ((Normalized Difference Vegetation Index (NDVI), Vegetation Water Content (VWC), and Enhanced Vegetation Index (EVI)) were selected as the vegetation index, and the equation is as follows:
NDVI = (b5 − b4) / (b5 + b4)
VWC = (b5 − b6) / (b5 + b6) × 2.15 + 0.32
EVI = 2.5 × (b5 − b4) / (b5 + b4 − 7.5 × b2 + 1)
In the equation, b2 represents the value of band 2; b4 represents the value of band 4; b5 represents the value of band 5; b6 represents the value of band 6.

2.4. Methodology

SAR is the main way to detect the surface soil moisture by active microwave remote sensing. The basic principle is that the backscattering coefficient obtained by microwave remote sensing has a good correlation with the soil dielectric constant property, and the dielectric constant of the soil is mainly controlled by soil moisture. Therefore, based on this, the surface soil moisture can be effectively derived [48].
In this paper, the Dubois Model and the Ratio Model were used to derive soil moisture in the Ugan-Kuqa River Delta Oasis. The Ratio Model is used to simulate the vegetation backscattering, and the Dubois is used to simulate the bare soil backscattering. The implementation of the two models’ soil moisture retrieval can be seen in Figure 4, and detailed descriptions of the two models are introduced in the following sections.
Step 1: based on the Dubois model, the in situ surface roughness, and the other parameters, such as the in situ soil moisture, related length, and system configuration of the Sentinel-1 data were used to simulate the bare soil-backscattering coefficient under two polarizations.
Step 2: the backscattering coefficient of bare soil simulated by the Dubois model, the total backscattering coefficient of the Sentinel-1 image, and three optical vegetation indices (VWC, NDVI, and EVI) were introduced into the ratio model, and the unknown coefficient in the ratio model was calculated using the least squares algorithm.
Step 3: based on the in situ soil moisture and the backscattering coefficient of the bare soil simulated by the ration model, the empirical equations for the two polarizations were established (Equation (4)).
d σ p p s o i l o + e = m v
In the equation, d and e represents the unknown coefficient, m v represents the soil moisture, and p p represents the polarization mode.
Step 4: based on the ratio model, because the radar-backscattering coefficient is known, the soil-backscattering coefficient can be obtained. Combined with Equation (4), the soil moisture was derived under two polarizations. Using Equation (5), the soil moisture was derived.
m v = ( m v v v + m v h h ) / 2

2.4.1. Dubois Model

In 1995, an empirical model was conducted based on microwave scatterometry data (including multiple frequencies, multi-polarization, multi-angle) under different rough conditions [31].
This model can be used in the following range:
30 < θ < 65
0.3   cm < S < 3.0   cm
1.5   GHz < f < 11   GHz
In the equation, θ represents the radar incident angle, S represents the root mean square height, f represents the radar frequency. Under this condition, the simulated scatter coefficient has good accuracy, but under a tough soil roughness, it is not suitable to use this model to simulate.
The equation of the Dubois model is:
σ h h 0 = 10 2.75 cos 1.5 θ sin 5 θ 10 0.028 ε tan θ ( k s sin θ ) 1.4 λ 0.7
σ v v 0 = 10 2.35 cos 3 θ sin 3 θ 10 0.046 ε tan θ ( k s sin θ ) 1.1 λ 0.7
In the equation, σ h h 0 represents the backscattering coefficient of the HH polarization mode, σ v v 0 represents the backscattering coefficient of the VV polarization mode, k represents the radar wave number, λ represents the wavelength, s represents the root mean square height, and θ represents the radar incident angle.
Based on a complex electromagnetic scattering model, the theoretical model was established, which has a certain physical basis. It can be widely applied to different SAR sensors and has good adaptability. Compared with the theoretical model, the empirical model has advantages, which can simulate the backscattering coefficient of the soil beyond the threshold range of the theoretical model, and the model has a better precision under a low vegetation coverage. Despite the above advantages, the model still has the condition that the soil moisture is underestimated [49,50], and the model ignores the influence of the correlation length on the backscattering coefficient [51], so the universality of the constructed model needs to be improved.

2.4.2. Topp Model

In 1980, an empirical relationship between soil water content and soil dielectric constant was conducted [52], based on the in situ data, which does not consider the effects of soil texture and temperature on soil moisture, and also ignores the influence of soil imaginary part on soil moisture. The empirical model relational equation is as follows:
ε = 3.03 + 9.3 m v + 146.0 m v 2 76.7 m v 3
m v = 5.3 × 10 2 + 2.92 × 10 2 ε 5.5 × 10 4 ( ε ) 2 + 4.3 × 10 6 ( ε ) 3
In the equation, m v represents the soil water content, and ε represents the real part of the soil dielectric constant.
According to the above equation, it can be seen that the soil moisture can be obtained only when the real part of the soil dielectric constant is required. The model has simple expression and simple calculation, which is suitable for soil under lighter soil roughness, and in an oasis region, the imaginary part of the dielectric constant is negligible.

2.4.3. Ratio Model

In 2008, Joseph et al. used a microwave scatterometer to monitor the growth of corn throughout its growth cycle; the ratio model was conducted [53]. The model assumes that the ratio of the backscattering coefficient of the soil to the backscattering coefficient of the radar is only related to the vegetation indices under certain conditions of the microwave sensor. Therefore, the ratio model can be used to eliminate the influence of vegetation canopy on radar backscattering coefficient. Based on this, under vegetation coverage, this model can be used to establish the empirical relationship between the backscattering coefficient and soil moisture, and then derive the soil moisture. The ratio model replaces the surface backscattering, secondary scattering, and multiple scattering of vegetation by the total backscattering of the vegetation canopy. The ratio model can be simplified as follow:
σ p p s o i l o σ p p o = f ( v e g e t a t i o n , a n t e n n a   c o n f i g u r a t i o n )
In the equation, σ p p s o i l o represents the backscattering coefficient of bare soil, σ p p o represents the total backscattering coefficient of the radar, and p p represents the polarization mode. When the radar system parameters are certain, the ratio model is simplified as follows:
σ p p s o i l o σ p p o = f ( v e g e t a t i o n )
In order to better separate the scattering contribution of vegetation, many researchers have proposed different optimal models (such as linear models, nonlinear models, etc.) to characterize the scattering contribution of vegetation. Joseph use quadratic and exponential equations to characterize vegetation coverage, in which the ratio model was modified [53].
a V 2 + exp ( b V ) = f ( v e g e t a t i o n )
In 2010, a constant term was added to the model by Joseph, which have a better usage under sparse vegetation coverage or the primary stage of vegetation growth [54].
a V 2 + exp ( b V ) + c = f ( v e g e t a t i o n )
In 2012, Prakash et al. [55] used a polynomial function to characterize the scattering component of vegetation.
a V 2 + b V + c = f ( v e g e t a t i o n )
Combined with the specific situation of the study area, where the inner oasis is mostly covered by low-lying crops and the oasis-desert ecotone region is under sparse vegetation coverage, we intend to use a combination of linear and power exponential functions to characterize vegetation scattering contributions. The equation is as follows:
a V 2 + b V c = f ( v e g e t a t i o n )
In the equation, a, b, and c all represent unknown coefficients, and V represents optical vegetation indices. There is no uniform standard for the best choice of optical vegetation indices.
In 2011, Lievens et al. [38] selected three vegetation indexes, including LAI, VWC, and NDVI, which were used to initialize the water cloud model (WCM) model, and the soil moisture was derived. After precision analysis, it was found that the accuracy of using LAI to derive the soil moisture was much higher than the other two vegetation indices, which proved that the leaves of vegetation affected vegetation scattering, and affected the interaction between radar signals and vegetation canopy. In 2015, Kweon et al. [56] modified the water-cloud model by adding new parameters (the average and standard deviation of leaf angle distribution) to accurately estimate the backscattering coefficients with the angular effect of scattering particles in a vegetation canopy. It proved that vegetation water content affects the penetration depth of radar, and further affects the interaction between radar signal and vegetation canopy. In 2012, Mattar et al. [42] chose the NDVI index, obtained from optical measurements, for removing the vegetation effect, and concluded with positive results.
Therefore, three different vegetation indices are used in the ratio model, respectively, VWC, NDVI, and EVI. The three vegetation indices are used to estimate the scattering contribution of vegetation to better derive soil moisture.

2.4.4. Minimum-Maximum (Min-Max) Normalization Approach

The Min-Max normalization is a linear transformation of the original data, so that the result falls in the interval [0,1], the equation used is:
S n o r m a l i z a t i o n   = S M i n M a x M i n
where S n o r m a l i z a t i o n   is normalized value of soil moisture, S is derived value of soil moisture, M a x is maximum value of soil moisture data,   M i n is minimum value of soil moisture data.

2.4.5. Accuracy Assessment

The correlation coefficient (R) is as follows:
R ( X , Y ) = C o v ( X , Y ) V a r [ X ] V a r [ Y ]
In the equation, C o v ( X , Y ) represents the covariance of X and Y, V a r [ X ] represents the variance of X, and V a r [ Y ] represents the variance of Y.
The root mean square error (RMSE) is as follows:
RMSE = 1 n = 1 n ( y i y i ^ ) 2
In the equation, n represents the number of samples,   y i represents the value of the i-th sample, y i ^ represents the mean value of the sample.

3. Results

3.1. Removing the Vegetation Effect

In this study, the ratio model was used to eliminate the backscatter contribution of vegetation, in which the growth of vegetation is characterized by three vegetation indices of VWC, NDVI, and EVI, and the three vegetation indices were derived by optical data. The linear and power function were used to characterize the ratio model, which was the ratio of the backscattering coefficient of the soil and of the imagery, and the function also use the least square method to solve the unknown parameters in the ratio model.
Figure 5, Figure 6 and Figure 7 show the process of removing the vegetation effect in July. Figure 5 shows the retrieval of the three vegetation indices. Figure 5a shows the VWC vegetation index, Figure 5b shows the NDVI vegetation index, and Figure 5c shows the EVI vegetation index. Figure 6 shows the ratio scatter plots of σ h h s o i l o / σ h h o to vegetation indices. Figure 6a shows the ratio scatter plot under VWC vegetation index, Figure 6b shows the ratio scatter plot under NDVI vegetation index, and Figure 6c shows the ratio scatter plot under EVI vegetation index. Figure 7 shows the ratio scatter plots of σ v v s o i l o / σ v v o to vegetation indices. Figure 7a shows the ratio scatter plot under VWC vegetation index, Figure 7b shows the ratio scatter plot under NDVI vegetation index, and Figure 7c shows the ratio scatter plot under EVI vegetation index.
As shown in Figure 5, Figure 6 and Figure 7, the ratio of the backscattering coefficient has a negative correlation with the vegetation indices, and the comprehensive function can well characterize the relationship between the two models. When the vegetation indices are small, the ratio of backscattering coefficient decreases rapidly with the increase of vegetation indices. When it reaches the critical point, the ratio of backscattering coefficient decreases slowly and the saturation tends to be gentle. This rule is most pronounced when the vegetation index is NDVI.
After analysis, we can conclude that the ratio of the backscattering coefficient shows a rapid decrease when the vegetation indices are small (the main reason being the influence of vegetation and soil roughness on the backscattering coefficient). The main reason for the saturation is that the attenuation of vegetation and soil roughness to backscatter coefficient tends to be saturated, which is the same as Joseph et al.’s research [53,54].
The unknown coefficients in the ratio equations of three different vegetation indices were obtained by the least square method (Table 2). So the regression equation of two polarizations under three vegetation indices were obtained. As shown in the regression equation, it can be found that coefficients a and c are negative values; coefficient b is positive value.

3.2. Constructing the Soil Moisture Retrieval Model

In order to evaluate the effectiveness of the ratio model of the soil moisture retrieval in this oasis, the relationship between the soil-backscattering coefficient of the ratio model and bare soil- backscattering coefficient of the Dubois model was established. Figure 8 shows the function of the backscattering coefficient of three vegetation indices under HH polarization mode. Figure 9 shows the function of the backscattering coefficient of three vegetation indices under VV polarization mode. Table 3 shows R and RMSE between the simulated values of the backscattering coefficient.
As shown in Figure 8 and Figure 9, it can be found that the backscattering coefficient simulated by the Dubois model and the ratio model have a fine linear relationship, and the scatter plots are basically clustered near the 1:1 line, which shows that the backscattering coefficient simulated by ratio model is feasible. As shown in Table 3, the quantitative analysis of the fitting results of the backscattering coefficient shows that when the vegetation index is VWC, the backscattering coefficient has the highest correlation, the root mean square error has the smallest, and the modeling effect is the best.
There are two functions of the Dubois model in this research. One is to simulate the backscattering coefficient of bare soil and solve the unknown coefficient in the ratio model. The other is to simulate the backscattering coefficient of bare soil and establish the empirical relationship of soil moisture. Therefore, according to the Dubois model, it can be concluded that when the soil roughness is certain, the backscattering coefficient of bare soil is controlled by soil moisture. Based on this, the Dubois model is used to establish the relationship between soil backscattering coefficient and soil moisture. Moreover, 80% of the in situ soil data (55 soil samples) were used to build the empirical relationship between the in situ soil moisture and backscattering coefficient. Figure 10 shows the scatter plot of the backscattering coefficient of the Dubois model and in situ soil moisture under HH and VV polarization modes.
Figure 10 shows that there is a linear relationship between the backscattering coefficient simulated by the Dubois model and the in situ soil moisture. Therefore, we used a linear regression model to represent the functional relationship between the backscattering coefficient and soil moisture, which as shown in following Equation (4). The unknown coefficient of the equation was solved and the results are shown in Table 4. The coefficient d represents the slope and the coefficient e represents the intercept. From Table 4, the regression equation can be concluded. The regression equation is m v = 0.96 σ h h s o i l o + 30.18 under HH polarization, the regression equation is m v = 0.92 σ v v s o i l o + 23.72 under VV polarization. As shown in Table 4, it can be found that the empirical relationship between the two polarizations and soil moisture are basically the same.

3.3. Soil Moisture Retrieval

The backscattering coefficient of the soil was obtained by using the ratio model; meanwhile, the empirical model (Equation (4)) was used to derive soil moisture based on the Sentinel-1 microwave data. In order to distinguish the changes of soil moisture, the results of surface soil moisture were normalized to [0,1] using a Min-Max normalization approach (Figure 11).
As shown in Figure 11, the soil moisture in the oasis-desert ecotone region (southeast) is low (below 0.2). In the oasis region (middle and north), it is moderate (about 0.3); this part has the largest proportion. In the oasis region (west), it is high (above 0.3). The soil moisture in the inner oasis is much higher than that in the oasis-desert ecotone region.
In order to see the difference between the soil moisture retrieval distribution, the pixel number of each interval was computed, and the percentage was computed. The results are listed in Table 5. In Table 5, it can be seen that the percentage remains highest when the soil moisture is between 0.2 and 0.4. When the soil moisture is between 0 and 0.2, the percentage is about 25%. When the soil moisture is higher than 0.4, the pixel number has a small percentage.
There were differences in the soil moisture retrieval under the three vegetation indices, which proved that the model conducted was suitable for deriving soil moisture in this oasis region. All of these findings can be seen in Figure 11.
Moreover, 20% of the in situ soil data (14 soil samples) were used to verify the accuracy of the soil moisture model constructed in this paper. The soil moisture simulated by the Dubois model were correlated with the in situ soil moisture under three vegetation indices. The results are shown in Figure 12, and their error metrics, i.e., bias, RMSE, are listed in Table 6.
Figure 12 shows that the correlation coefficient between the simulation soil moisture and in situ soil moisture under three different vegetation indices. The soil moisture retrieval accuracy is the highest under the VWC vegetation index, and the coefficient of determination reached 0.76, which was significantly correlated. When NDVI and EVI vegetation indices are used to eliminate vegetation scattering contribution, the soil moisture retrieval accuracy is not much different, and the coefficient of determination is about 0.73, which is slightly weaker than the result under the VWC vegetation index.
As shown in Table 6, the error metrics were calculated under three vegetation indices. A positive bias indicates that the retrieved soil moisture is systematically overestimated compared with the in situ measurements. The RMSE represents the absolute errors between the observations and simulations. In general, the accuracy of the soil moisture estimation is promising; the model built using VWC vegetation index works best. The inner oasis is the most impacted by the choice of the vegetation indices, and the land cover is agriculture.

4. Discussion

According to the analyses of the soil moisture retrieval accuracy, which is not high, there are uncertain factors in the retrieval process. The main aspects are as follows: (1) it is difficult to maintain the same time for field survey data, microwave data, and the optical data. During this period, it is assumed that the vegetation index remain unchanged, but July is a period of vigorous vegetation growth, which will lead to a miscalculation of vegetation scattering. (2) The parameter error of soil moisture is caused by the nonstandard sampling and the uncertainty of data processing in the laboratory. (3) There are systematic errors in the model itself, such as the backscattering coefficients simulated by the Dubois model and the ratio model, which will affect the accuracy of the soil moisture retrieval. Nevertheless, this method is still feasible if the purpose is to obtain the soil moisture in an oasis region.
Note that two models have equal importance to the soil moisture derived. The Dubois model can establish an empirical relationship and the ratio model can eliminate the vegetation scattering contribution. Among them, the ratio model is suitable for low shrubs, crops, and is not fit for arbor. However, when it comes to the soil moisture derived under forest coverage, the continuous Michigan Microwave Canopy Scattering (MIMICS) model is commonly used to simulate the scattering characteristics of forests, which is suitable for arbor. In this situation, different vegetation coverage could influence the detection results. Therefore, this method may not be suitable for temperate forests or tropical forests.

5. Conclusions

The soil moisture retrieval model was conducted based on the Topp dielectric constant model, the Dubois model and ratio method; then, the Sentinel-1 data were used to derive the soil moisture in an oasis region, and the feasibility of this method was verified by using the in situ soil moisture. Based on the findings of this study, the following conclusions can be drawn:
  • Based on the Dubois model, in situ soil roughness, the soil moisture, and radar system configuration parameters were used to simulate the backscattering coefficient of bare soil, which solved the unknown coefficient of the ratio model, on one hand, and verified the feasibility of the ratio model on the other hand.
  • The VWC, NDVI, and EVI were the three vegetation parameters that were selected to characterize the scattering contribution of vegetation; based on the ratio model, the least squares method was used to calibrate the unknown coefficient. Based on this model, the unknown coefficient was certain; the backscattering coefficient of bare soil was simulated and fitted with the backscattering coefficient, simulated by the Dubois model, which found that the ratio model could effectively simulate the scattering characteristics of bare soil.
  • The Dubois model was used to simulate the backscattering coefficient of the soil; an empirical relationship between it and the in situ soil moisture was established. Based on the model we conducted, the spatial soil moisture was successfully derived.
  • The in situ soil moisture was used to analyze the accuracy of the soil moisture derived in this oasis under three vegetation indices. The results show that the soil moisture retrieval accuracy was the highest under the VWC vegetation index, and the soil moisture retrieval accuracy was not much different under NDVI and EVI vegetation indices. The results show that the soil moisture model conducted on the Dubois model and ratio model is feasible.

Author Contributions

Conceptualization, S.H. and J.D.; funding acquisition, J.D.; software, B.L. and J.W.; methodology, B.L. and X.G.; formal analysis, S.H.; investigation, X.G. and J.Z. (Jie Zou); data curation, J.W. and J.Z. (Junyong Zhang); writing—original draft preparation, S.H.; project administration, S.H. and J.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China [No. 41961059 and No.41771470] and the Xinjiang University Excellent Doctoral Innovation Project (XJUBSCX-2016013).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kornelsen, K.C.; Coulibaly, P. Advances in soil moisture retrieval from synthetic aperture radar and hydrological applications. J. Hydrol. 2013, 476, 460–489. [Google Scholar] [CrossRef]
  2. Gao, F.; Wang, J.; Sun, C.; Wen, J. Advances in Study on Microwave Remote Sensing of Soil Moisture. Remote Sens. Technol. Appl. 2001, 16, 97–102. [Google Scholar]
  3. Procházka, P.; Hönig, V.; Maitah, M.; Pljučarská, I.; Kleindienst, J. Evaluation of Water Scarcity in Selected Countries of the Middle East. Water 2018, 10, 1482. [Google Scholar] [CrossRef] [Green Version]
  4. Zhe, Y.; Lei, T.; Yan, Q.; Chen, Z.; Dong, Y.J.C. A near-infrared reflectance sensor for soil surface moisture measurement. Comput. Electron. Agric. 2013, 99, 101–107. [Google Scholar]
  5. Lo, M.H.; Famiglietti, J.S. Irrigation in California’s Central Valley strengthens the southwestern US water cycle. Geophys. Res. Lett. 2013, 40, 301–306. [Google Scholar] [CrossRef] [Green Version]
  6. Hong, Z.; Zhang, W.; Yu, C.; Zhang, D.; Li, L.; Meng, L. SWCTI: Surface Water Content Temperature Index for Assessment of Surface Soil Moisture Status. Sensors 2018, 18, 2875. [Google Scholar] [CrossRef] [Green Version]
  7. Machiwal, D.; Gupta, A.; Jha, M.K.; Kamble, T. Analysis of trend in temperature and rainfall time series of an Indian arid region: Comparative evaluation of salient techniques. Theor. Appl. Climatol. 2019, 136, 301–320. [Google Scholar] [CrossRef]
  8. Jin, H.; Zhu, Q.; Zhao, X.; Zhang, Y. Simulation and Prediction of Climate Variability and Assessment of the Response of Water Resources in a Typical Watershed in China. Water 2016, 8, 490. [Google Scholar] [CrossRef] [Green Version]
  9. Sorooshian, S.; Li, J.; Hsu, K.L.; Gao, X. How significant is the impact of irrigation on the local hydroclimate in California’s Central Valley? Comparison of model results with ground and remote-sensing data. J. Geophys. Res. Atmos. 2011, 116. [Google Scholar] [CrossRef] [Green Version]
  10. Fisher, J.B.; Melton, F.; Middleton, E.; Hain, C.; Anderson, M.; Allen, R.; Mccabe, M.; Hook, S.; Baldocchi, D.; Townsend, P.A. The Future of Evapotranspiration: Global requirements for ecosystem functioning, carbon and climate feedbacks, agricultural management, and water resources. Water Resour. Res. 2017, 53, 2618–2626. [Google Scholar] [CrossRef]
  11. Qiaozhen, M.U.; Zhao, M.; Steven, W. Improvements to a MODIS global terrestrial evapotranspiration algorithm. Remote Sens. Environ. 2011, 115, 1781–1800. [Google Scholar]
  12. Lei, F.; Crow, W.T.; Kustas, W.P.; Dong, J.; Yang, Y.; Knipper, K.R.; Anderson, M.C.; Gao, F.; Notarnicola, C.; Greifeneder, F.; et al. Data assimilation of high-resolution thermal and radar remote sensing retrievals for soil moisture monitoring in a drip-irrigated vineyard. Remote Sens. Environ. 2020, 239, 111622. [Google Scholar] [CrossRef] [PubMed]
  13. Zhao, J.; Zhong, Y.; Hu, X.; Wei, L.; Zhang, L. A robust spectral-spatial approach to identifying heterogeneous crops using remote sensing imagery with high spectral and spatial resolutions. Remote Sens. Environ. 2020, 239, 111605. [Google Scholar] [CrossRef]
  14. Tong, X.-Y.; Xia, G.-S.; Lu, Q.; Shen, H.; Li, S.; You, S.; Zhang, L. Land-cover classification with high-resolution remote sensing images using transferable deep models. Remote Sens. Environ. 2020, 237, 111322. [Google Scholar] [CrossRef] [Green Version]
  15. Zhang, K.; Chao, L.-J.; Wang, Q.-Q.; Huang, Y.-C.; Liu, R.-H.; Hong, Y.; Tu, Y.; Qu, W.; Ye, J.-Y. Using multi-satellite microwave remote sensing observations for retrieval of daily surface soil moisture across China. Water Sci. Eng. 2019, 12, 85–97. [Google Scholar] [CrossRef]
  16. Meyer, T.; Weihermüller, L.; Vereecken, H.; Jonard, F. Vegetation Optical Depth and Soil Moisture Retrieved from L-Band Radiometry over the Growth Cycle of a Winter Wheat. Remote Sens. 2018, 10, 1637. [Google Scholar] [CrossRef] [Green Version]
  17. Bai, X.; He, B.; Li, X.; Zeng, J.; Wang, X.; Wang, Z.; Zeng, Y.; Su, Z. First Assessment of Sentinel-1A Data for Surface Soil Moisture Estimations Using a Coupled Water Cloud Model and Advanced Integral Equation Model over the Tibetan Plateau. Remote Sens. 2017, 9, 714. [Google Scholar] [CrossRef] [Green Version]
  18. Zhang, L.; Meng, Q.; Yao, S.; Wang, Q.; Zeng, J.; Zhao, S.; Ma, J. Soil Moisture Retrieval from the Chinese GF-3 Satellite and Optical Data over Agricultural Fields. Sensors 2018, 18, 2675. [Google Scholar] [CrossRef] [Green Version]
  19. Rahman, M.S.; Di, L.; Yu, E.; Lin, L.; Zhang, C.; Tang, J. Rapid Flood Progress Monitoring in Cropland with NASA SMAP. Remote Sens. 2019, 11, 191. [Google Scholar] [CrossRef] [Green Version]
  20. Gherboudj, I.; Magagi, R.; Berg, A.A.; Toth, B. Soil moisture retrieval over agricultural fields from multi-polarized and multi-angular RADARSAT-2 SAR data. Remote Sens. Environ. 2011, 115, 33–43. [Google Scholar] [CrossRef]
  21. Chai, X.; Zhang, T.; Shao, Y.; Gong, H.; Liu, L.; Xie, K. Modeling and Mapping Soil Moisture of Plateau Pasture Using RADARSAT-2 Imagery. Remote Sens. 2015, 7, 1279–1299. [Google Scholar] [CrossRef] [Green Version]
  22. Pratola, C.; Barrett, B.; Gruber, A.; Dwyer, E. Quality Assessment of the CCI ECV Soil Moisture Product Using ENVISAT ASAR Wide Swath Data over Spain, Ireland and Finland. Remote Sens. 2015, 7, 15388–15423. [Google Scholar] [CrossRef] [Green Version]
  23. Dabrowska-Zielinska, K.; Budzynska, M.; Tomaszewska, M.; Malinska, A.; Gatkowska, M.; Bartold, M.; Malek, I. Assessment of Carbon Flux and Soil Moisture in Wetlands Applying Sentinel-1 Data. Remote Sens. 2016, 8, 756. [Google Scholar] [CrossRef] [Green Version]
  24. Alexakis, D.D.; Mexis, F.-D.K.; Vozinaki, A.-E.K.; Daliakopoulos, I.N.; Tsanis, I.K. Soil Moisture Content Estimation Based on Sentinel-1 and Auxiliary Earth Observation Products. A Hydrological Approach. Sensors 2017, 17, 1455. [Google Scholar] [CrossRef] [Green Version]
  25. Gao, Q.; Zribi, M.; Escorihuela, M.J.; Baghdadi, N. Synergetic Use of Sentinel-1 and Sentinel-2 Data for Soil Moisture Mapping at 100 m Resolution. Sensors 2017, 17, 1966. [Google Scholar] [CrossRef] [Green Version]
  26. Liu, Z.; Li, P.; Yang, J. Soil Moisture Retrieval and Spatiotemporal Pattern Analysis Using Sentinel-1 Data of Dahra, Senegal. Remote Sens. 2017, 9, 1197. [Google Scholar] [CrossRef] [Green Version]
  27. Huang, S.; Ding, J.; Zou, J.; Liu, B.; Chen, W. Soil Moisture Retrival Based on Sentinel-1 Imagery under Sparse Vegetation Coverage. Sensors 2019, 19, 589. [Google Scholar] [CrossRef] [Green Version]
  28. Azza, G.; Mehrez, Z.; Nicolas, B.; Bernard, M.; Pascal, F.; Zohra, C.J.R.S. Retrieval of Both Soil Moisture and Texture Using TerraSAR-X Images. Remote Sens. 2015, 7, 10098–10116. [Google Scholar]
  29. Wang, C.; Zhang, Z.; Paloscia, S.; Zhang, H.; Wu, F.; Wu, Q. Permafrost Soil Moisture Monitoring Using Multi-Temporal TerraSAR-X Data in Beiluhe of Northern Tibet, China. Remote Sens. 2018, 10, 1577. [Google Scholar] [CrossRef] [Green Version]
  30. Oh, Y.; Sarab, K.; Ulaby, F.T. Semi-Empirical Model of the Ensemble-Averaged Differential Mueller Matrix for Microwave Backscattering from Bare Soil Surfaces. IEEE Trans. Geosci. Remote Sens. 2002, 40, 1348–1355. [Google Scholar] [CrossRef] [Green Version]
  31. Dubois, P.C.; Van Zyl, J.; Engman, T. Measuring soil moisture with imaging radars. IEEE Trans. Geosci. Remote Sens. 1995, 33, 915–926. [Google Scholar] [CrossRef] [Green Version]
  32. Fung, A.K.; Li, Z.; Chen, K.S. Backscattering from a randomly rough dielectric surface. IEEE Trans. Geosci. Remote Sens. 1992, 30, 356–369. [Google Scholar] [CrossRef]
  33. Pettinato, S.; Santi, E.; Paloscia, S.; Pampaloni, P.; Fontanelli, G. The Intercomparison of X-Band SAR Images from COSMO? SkyMed and TerraSAR-X Satellites: Case Studies. Remote Sens. 2013, 5, 2928–2942. [Google Scholar] [CrossRef] [Green Version]
  34. Choker, M.; Baghdadi, N.; Zribi, M.; El Hajj, M.; Paloscia, S.; Verhoest, N.E.C.; Lievens, H.; Mattia, F. Evaluation of the Oh, Dubois and IEM Backscatter Models Using a Large Dataset of SAR Data and Experimental Soil Measurements. Water 2017, 9, 38. [Google Scholar] [CrossRef]
  35. Wu, T.D.; Chen, K.S.; Shi, J.; Fung, A.K. A transition model for the reflection coefficient in surface scattering. IEEE Trans. Geosci. Remote Sens. 1998, 39, 2040–2050. [Google Scholar]
  36. He, L.; Jing, M.C.; Chen, K.S. Simulation and SMAP Observation of Sun-Glint over the Land Surface at the L-Band. IEEE Trans. Geosci. Remote Sens. 2017, 55, 2589–2604. [Google Scholar] [CrossRef]
  37. Chen, K.S.; Wu, T.D.; Tsang, L.; Li, Q.; Shi, J.; Fung, A.K. Emission of rough surfaces calculated by the integral equation method with comparison to three-dimensional moment method simulations. IEEE Trans. Geosci. Remote Sens. 2003, 41, 90–101. [Google Scholar] [CrossRef]
  38. Lievens, H.; Verhoest, N.E. On the Retrieval of Soil Moisture in Wheat Fields from L-Band SAR Based on Water Cloud Modeling, the IEM, and Effective Roughness Parameters. IEEE Geosci. Remote Sens. Lett. 2011, 8, 740–744. [Google Scholar] [CrossRef]
  39. He, B.; Xing, M.; Bai, X. A Synergistic Methodology for Soil Moisture Estimation in an Alpine Prairie Using Radar and Optical Satellite Data. Remote Sens. 2014, 6, 10966–10985. [Google Scholar] [CrossRef] [Green Version]
  40. Wang, C.; Qi, J. Biophysical estimation in tropical forests using JERS-1 SAR and VNIR imagery. II. Aboveground woody biomass. Int. J. Remote Sens. 2008, 29, 6827–6849. [Google Scholar] [CrossRef]
  41. Notarnicola, C.; Angiulli, M.; Posa, F. Use of radar and optical remotely sensed data for soil moisture retrieval over vegetated areas. IEEE Trans. Geosci. Remote Sens. 2006, 44, 925–935. [Google Scholar] [CrossRef]
  42. Mattar, C.; Wigneron, J.-P.; Sobrino, J.A.; Novello, N.; Calvet, J.C.; Albergel, C.; Richaume, P.; Mialon, A.; Guyon, D.; Jimenez-Munoz, J.C.; et al. A combined optical–microwave method to retrieve soil moisture over vegetated areas. IEEE Trans. Geosci. Remote Sens. 2012, 50, 1404–1413. [Google Scholar] [CrossRef]
  43. Xu, C.; Qu, J.J.; Hao, X.; Wu, D. Monitoring Surface Soil Moisture Content over the Vegetated Area by Integrating Optical and SAR Satellite Observations in the Permafrost Region of Tibetan Plateau. Remote Sens. 2020, 12, 183. [Google Scholar] [CrossRef] [Green Version]
  44. Santi, E.; Paloscia, S.; Pettinato, S.; Notarnicola, C.; Pasolli, L.; Pistocchi, A. Comparison between SAR Soil Moisture Estimates and Hydrological Model Simulations over the Scrivia Test Site. Remote Sens. 2013, 5, 4961–4976. [Google Scholar] [CrossRef] [Green Version]
  45. Li, J.; Wang, S. Using SAR-Derived Vegetation Descriptors in a Water Cloud Model to Improve Soil Moisture Retrieval. Remote Sens. 2018, 10, 1370. [Google Scholar] [CrossRef] [Green Version]
  46. Attarzadeh, R.; Amini, J.; Notarnicola, C.; Greifeneder, F. Synergetic Use of Sentinel-1 and Sentinel-2 Data for Soil Moisture Mapping at Plot Scale. Remote Sens. 2018, 10, 1285. [Google Scholar] [CrossRef] [Green Version]
  47. Bousbih, S.; Zribi, M.; El Hajj, M.; Baghdadi, N.; Lili-Chabaane, Z.; Gao, Q.; Fanise, P. Soil Moisture and Irrigation Mapping in A Semi-Arid Region, Based on the Synergetic Use of Sentinel-1 and Sentinel-2 Data. Remote Sens. 2018, 10, 1953. [Google Scholar] [CrossRef] [Green Version]
  48. Ulaby, F.T.; Batlivala, P.P.; Dobson, M.C. Microwave Backscatter Dependence on Surface Roughness, Soil Moisture, and Soil Texture: Part I-Bare Soil. IEEE Trans. Geosci. Electron. 1978, 16, 286–295. [Google Scholar] [CrossRef]
  49. Panciera, R.; Tanase, M.A.; Lowell, K.; Walker, J.P. Evaluation of IEM, Dubois, and Oh radar backscatter models using airborne L-band SAR. IEEE Trans. Geosci. Remote Sens. 2013, 52, 4966–4979. [Google Scholar] [CrossRef]
  50. Merzouki, A.; McNairn, H.; Pacheco, A. Evaluation of the Dubois, Oh, and IEM radar backscatter models over agricultural fields using C-band RADARSAT-2 SAR image data. Can. J. Remote Sens. 2010, 36, S274–S286. [Google Scholar] [CrossRef]
  51. Oh, Y.; Sarabandi, K.; Ulaby, F.T. An empirical model and an inversion technique for radar scattering from bare soil surfaces. IEEE Trans. Geosci. Remote Sens. 1992, 30, 370–381. [Google Scholar] [CrossRef]
  52. Topp, G.C.; Davis, J.; Annan, A.P. Electromagnetic determination of soil water content: Measurements in coaxial transmission lines. Water Resour. Res. 1980, 16, 574–582. [Google Scholar] [CrossRef] [Green Version]
  53. Joseph, A.T.; van der Velde, R.; O’Neill, P.E.; Lang, R.H.; Gish, T. Soil moisture retrieval during a corn growth cycle using L-band (1.6 GHz) radar observations. IEEE Trans. Geosci. Remote Sens. 2008, 46, 2365–2374. [Google Scholar] [CrossRef] [Green Version]
  54. Joseph, A.; van der Velde, R.; O’neill, P.; Lang, R.; Gish, T. Effects of corn on C-and L-band radar backscatter: A correction method for soil moisture retrieval. Remote Sens. Environ. 2010, 114, 2417–2430. [Google Scholar] [CrossRef]
  55. Prakash, R.; Singh, D.; Pathak, N.P. A fusion approach to retrieve soil moisture with SAR and optical data. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2011, 5, 196–206. [Google Scholar] [CrossRef]
  56. Kweon, S.K.; Oh, Y. A Modified Water-Cloud Model with Leaf Angle Parameters for Microwave Backscattering from Agricultural Fields. IEEE Trans. Geosci. Remote Sens. 2015, 53, 2802–2809. [Google Scholar] [CrossRef]
Figure 1. Geographical location of the study area: (a) the Xinjiang Uygur Autonomous Region; (b) the Ugan-Kuqa River Delta Oasis with standard false-color composite (data source: Landsat-8 Operational Land Imager (OLI), July 14, 2017).
Figure 1. Geographical location of the study area: (a) the Xinjiang Uygur Autonomous Region; (b) the Ugan-Kuqa River Delta Oasis with standard false-color composite (data source: Landsat-8 Operational Land Imager (OLI), July 14, 2017).
Remotesensing 12 01358 g001
Figure 2. Land cover of study area (data source: Landsat-8 OLI, July 14, 2017).
Figure 2. Land cover of study area (data source: Landsat-8 OLI, July 14, 2017).
Remotesensing 12 01358 g002
Figure 3. Soil roughness data processing: (a) sketch of roughness measurement; (b) image vectorization; (c) manual correction processing.
Figure 3. Soil roughness data processing: (a) sketch of roughness measurement; (b) image vectorization; (c) manual correction processing.
Remotesensing 12 01358 g003
Figure 4. Flowchart for the soil moisture retrieval.
Figure 4. Flowchart for the soil moisture retrieval.
Remotesensing 12 01358 g004
Figure 5. Three vegetation indices retrieval in July. (a) Vegetation Water Content (VWC); (b) Normalized Difference Vegetation Index (NDVI); (c) Enhanced Vegetation Index (EVI).
Figure 5. Three vegetation indices retrieval in July. (a) Vegetation Water Content (VWC); (b) Normalized Difference Vegetation Index (NDVI); (c) Enhanced Vegetation Index (EVI).
Remotesensing 12 01358 g005
Figure 6. The ratio scatter plots of σ h h s o i l o / σ h h o to vegetation indices (a) VWC; (b) NDVI; (c) EVI.
Figure 6. The ratio scatter plots of σ h h s o i l o / σ h h o to vegetation indices (a) VWC; (b) NDVI; (c) EVI.
Remotesensing 12 01358 g006
Figure 7. The ratio scatter plots of σ v v s o i l o / σ v v o to vegetation indices (a) VWC; (b) NDVI; (c) EVI.
Figure 7. The ratio scatter plots of σ v v s o i l o / σ v v o to vegetation indices (a) VWC; (b) NDVI; (c) EVI.
Remotesensing 12 01358 g007
Figure 8. Scatter plot of backscattering coefficient of HH polarization of Dubois model and ratio model (a) VWC; (b) NDVI; (c) EVI.
Figure 8. Scatter plot of backscattering coefficient of HH polarization of Dubois model and ratio model (a) VWC; (b) NDVI; (c) EVI.
Remotesensing 12 01358 g008
Figure 9. Scatter plot of backscattering coefficient of VV polarization of Dubois model and ratio model (a) VWC; (b) NDVI; (c) EVI.
Figure 9. Scatter plot of backscattering coefficient of VV polarization of Dubois model and ratio model (a) VWC; (b) NDVI; (c) EVI.
Remotesensing 12 01358 g009
Figure 10. The scatter plot of the backscattering coefficient of the Dubois model and in situ soil moisture (a) HH polarization (b) VV polarization.
Figure 10. The scatter plot of the backscattering coefficient of the Dubois model and in situ soil moisture (a) HH polarization (b) VV polarization.
Remotesensing 12 01358 g010
Figure 11. The spatial of soil moisture (a) VWC; (b) NDVI; (c) EVI.
Figure 11. The spatial of soil moisture (a) VWC; (b) NDVI; (c) EVI.
Remotesensing 12 01358 g011aRemotesensing 12 01358 g011b
Figure 12. The correlation of the soil moisture between simulation and measured values (a) VWC; (b) NDVI; (c) EVI.
Figure 12. The correlation of the soil moisture between simulation and measured values (a) VWC; (b) NDVI; (c) EVI.
Remotesensing 12 01358 g012aRemotesensing 12 01358 g012b
Table 1. Features of the land cover classes.
Table 1. Features of the land cover classes.
ClassSitesDescription
Agriculture21Planting crops
Natural vegetation 19Almost all vegetation is shrubs with a high vegetation coverage
Bare land10No vegetation cover
Water5Rivers, reservoirs
Salinized14No vegetation cover, but there are salt shells on the surface
Table 2. The fitting result of the unknown coefficient of the ratio equation.
Table 2. The fitting result of the unknown coefficient of the ratio equation.
Vegetation IndicesHH PolarizationVV Polarization
abcabc
VWC−0.231.15−0.38−0.261.13−0.40
f = 0.23 V 2 + 1.15 V 0.38 f = 0.26 V 2 + 1.13 V 0.40
NDVI−0.161.02−0.02−0.151.01−0.24
f = 0.16 V 2 + 1.02 V 0.02 f = 0.15 V 2 + 1.01 V 0.24
EVI−0.090.99−0.23−0.100.98−0.24
f = 0.09 V 2 + 0.99 V 0.23 f = 0.10 V 2 + 0.98 V 0.24
Table 3. Fitting result of backscattering coefficient.
Table 3. Fitting result of backscattering coefficient.
SARVegetation IndexHH PolarizationVV Polarization
RRMSERRMSE
Sentinel-1VWC0.833.020.733.65
NDVI0.813.260.723.79
EVI0.813.300.723.84
Table 4. The unknown coefficient of empirical model.
Table 4. The unknown coefficient of empirical model.
deRRMSE
HH polarization0.9630.180.754.56
VV polarization0.9223.720.723.27
Table 5. Pixel number of each interval of soil moisture map.
Table 5. Pixel number of each interval of soil moisture map.
VWCNDVIEVI
0–0.21,509,041 (24.36%)1,577,622 (25.47%)1,662,995 (26.85%)
0.2–0.43,853,753 (62.22%)3,749,477 (60.53%)3,703,183 (59.79%)
0.4–0.6818,423 (13.21%)841,098 (13.58%)799,751 (12.91%)
0.6–0.812,282 (0.20%)25,158 (0.41%)27,313 (0.44%)
0.8–1151 (0.01%)295 (0.01%)408 (0.01%)
Count6,193,6506,193,6506,193,650
Table 6. Statistical metrics between in situ and simulation soil moisture.
Table 6. Statistical metrics between in situ and simulation soil moisture.
ModelCoefficient of Determination (R2)RMSE (vol. %)Bias
VWC0.762.040.86
NDVI0.732.431.09
EVI0.742.391.14

Share and Cite

MDPI and ACS Style

Huang, S.; Ding, J.; Liu, B.; Ge, X.; Wang, J.; Zou, J.; Zhang, J. The Capability of Integrating Optical and Microwave Data for Detecting Soil Moisture in an Oasis Region. Remote Sens. 2020, 12, 1358. https://doi.org/10.3390/rs12091358

AMA Style

Huang S, Ding J, Liu B, Ge X, Wang J, Zou J, Zhang J. The Capability of Integrating Optical and Microwave Data for Detecting Soil Moisture in an Oasis Region. Remote Sensing. 2020; 12(9):1358. https://doi.org/10.3390/rs12091358

Chicago/Turabian Style

Huang, Shuai, Jianli Ding, Bohua Liu, Xiangyu Ge, Jinjie Wang, Jie Zou, and Junyong Zhang. 2020. "The Capability of Integrating Optical and Microwave Data for Detecting Soil Moisture in an Oasis Region" Remote Sensing 12, no. 9: 1358. https://doi.org/10.3390/rs12091358

APA Style

Huang, S., Ding, J., Liu, B., Ge, X., Wang, J., Zou, J., & Zhang, J. (2020). The Capability of Integrating Optical and Microwave Data for Detecting Soil Moisture in an Oasis Region. Remote Sensing, 12(9), 1358. https://doi.org/10.3390/rs12091358

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop