Satellite Soil Moisture for Agricultural Drought Monitoring : Assessment of SMAP-Derived Soil Water Deficit Index in Xiang River Basin , China

Agricultural drought can have long-lasting and harmful impacts on both the ecosystem and economy. Therefore, it is important to monitor and predict agricultural drought accurately. Soil moisture is the key variable to define the agricultural drought index. However, in situ soil moisture observations are inaccessible in many areas of the world. Remote sensing techniques enrich the surface soil moisture observations at different tempo-spatial resolutions. In this study, the Level 2 L-band radiometer soil moisture dataset was used to estimate the Soil Water Deficit Index (SWDI). The Soil Moisture Active Passive (SMAP) dataset was evaluated with the soil moisture dataset obtained from the China Land Soil Moisture Data Assimilation System (CLSMDAS). The SMAP-derived SWDI (SMAP_SWDI) was compared with the atmospheric water deficit (AWD) calculated with precipitation and evapotranspiration from meteorological stations. Drought monitoring and comparison were accomplished at a weekly scale for the growing season (April to November) from 2015 to 2017. The results were as follows: (1) in terms of Pearson correlation coefficients (R-value) between SMAP and CLSMDAS, around 70% performed well and only 10% performed poorly at the grid scale, and the R-value was 0.62 for the whole basin; (2) severe droughts mainly occurred from mid-June to the end of September from 2015 to 2017; (3) severe droughts were detected in the southern and northeastern Xiang River Basin in mid-May of 2015, and in the northern basin in early August of 2016 and end of November 2017; (4) the values of percentage of drought weeks gradually decreased from 2015 to 2017, and increased from the northeast to the southwest of the basin in 2015 and 2016; and (5) the average value of R and probability of detection between SMAP_SWDI and AWD were 0.6 and 0.79, respectively. These results show SMAP has acceptable accuracy and good performance for drought monitoring in the Xiang River Basin.


Introduction
Drought is a natural disaster that occurs with high frequency and which has long-lasting impacts on agriculture production, the ecological environment, and the economy [1].Agricultural drought can occur in many parts of the world but usually develops slowly and causes widespread devastation and economic loss [2].During the growing seasons from 1988 to 2001, Canada suffered more than 5 billion dollars in economic losses per year due to agricultural drought [3].Nebraska, USA, experienced moisture data to establish a support vector regression model with climate index to predict agricultural drought in the Xiang River Basin.The results showed that the second half of the years studied (from 1963 to 2012) suffered more drought events than the first half for all categories of droughts.By means of mutation analysis, wavelet analysis, variance index, empirical orthogonal functions, and rotated empirical orthogonal functions, Zhang et al. [43] analyzed the changes in drought and flood events in the Xiang River Basin over 48 years from 1961 to 2008.The results showed that the frequency of droughts and floods had declined since the 21st century, and this mutation had begun during the 1980s.Ma et al. [45] studied the past precipitation and temperature trends of the Xiang River Basin, and made future predictions for 2021 to 2050.Their study pointed out that more droughts will possibly occur in the future due to an insignificant increase in precipitation and significant increase in temperature.Therefore, it is of great significance to identify the temporal-spatial distribution characteristics of dry and wet areas and to monitor agricultural drought in the Xiang River Basin for disaster prevention and economic development of this basin as well as the Yangtze River Basin.
Remote Sens. 2019, 11, x FOR PEER REVIEW 3 of 21 droughts.By means of mutation analysis, wavelet analysis, variance index, empirical orthogonal functions, and rotated empirical orthogonal functions, Zhang et al. [43] analyzed the changes in drought and flood events in the Xiang River Basin over 48 years from 1961 to 2008.The results showed that the frequency of droughts and floods had declined since the 21st century, and this mutation had begun during the 1980s.Ma et al. [45] studied the past precipitation and temperature trends of the Xiang River Basin, and made future predictions for 2021 to 2050.Their study pointed out that more droughts will possibly occur in the future due to an insignificant increase in precipitation and significant increase in temperature.Therefore, it is of great significance to identify the temporal-spatial distribution characteristics of dry and wet areas and to monitor agricultural drought in the Xiang River Basin for disaster prevention and economic development of this basin as well as the Yangtze River Basin.Although there are many studies on drought monitoring in the Xiang River Basin, few studies have evaluated the accuracy of SMAP soil moisture using assimilated soil moisture data from the China Land Soil Moisture Data Assimilation System.Moreover, few studies have used the soil water deficit index obtained from SMAP to predict agricultural drought and to reveal the temporal-spatial characteristics of agricultural drought in this basin in the past three years.In addition, although some studies have been conducted to evaluate the efficiency of SMAP on drought monitoring around world [12,37,46], due to the high spatial variation and uncertainty of satellite-based soil moisture data, the potential capability of applying SMAP for drought monitoring on this specific region should be evaluated.Therefore, this study is aimed at investigating the accuracy of SMAP soil moisture data and the efficiency of the soil water deficit index obtained from SMAP used for drought monitoring in the Xiang River Basin.The remaining structure of the paper is organized as follows: the descriptions of in situ observations and SMAP soil moisture data are presented in Section 2; the methodology and evaluation indices are introduced in Section 3; Section 4 provides the results and discussion; and conclusions are drawn in Section 5.

Study Area and In Situ Data
The Xiang River is one of the largest tributaries of the Yangtze River.This river is located in Southern China (24.5°-28.5°N;111.0°-114.5°W),which is the largest river in Hunan Province and the Dongting Lake water system.The location of the Xiang River Basin and the distribution of the Although there are many studies on drought monitoring in the Xiang River Basin, few studies have evaluated the accuracy of SMAP soil moisture using assimilated soil moisture data from the China Land Soil Moisture Data Assimilation System.Moreover, few studies have used the soil water deficit index obtained from SMAP to predict agricultural drought and to reveal the temporal-spatial characteristics of agricultural drought in this basin in the past three years.In addition, although some studies have been conducted to evaluate the efficiency of SMAP on drought monitoring around world [12,37,46], due to the high spatial variation and uncertainty of satellite-based soil moisture data, the potential capability of applying SMAP for drought monitoring on this specific region should be evaluated.Therefore, this study is aimed at investigating the accuracy of SMAP soil moisture data and the efficiency of the soil water deficit index obtained from SMAP used for drought monitoring in the Xiang River Basin.The remaining structure of the paper is organized as follows: the descriptions of in situ observations and SMAP soil moisture data are presented in Section 2; the methodology and evaluation indices are introduced in Section 3; Section 4 provides the results and discussion; and conclusions are drawn in Section 5.

Study Area and In Situ Data
The Xiang River is one of the largest tributaries of the Yangtze River.This river is located in Southern China (24.5 • -28.5 • N; 111.0 • -114.5 • W), which is the largest river in Hunan Province and the Dongting Lake water system.The location of the Xiang River Basin and the distribution of the meteorological stations are shown in Figure 1.The Xiangtan discharge station is the control station of the Xiang River Basin.The surface area of the basin above the Xiangtan Station is about 82,375 km 2 .The Xiang River Basin is dominated by the Pacific monsoon climate.The average rainfall is 1400-1700 mm and the average temperature is 17 • C [1].The uneven distribution of rainfall has led to frequent floods and droughts in the Xiang River Basin [1].The area of cultivated land in the basin is 1.735 million hm 2 , which is a relatively well-developed agricultural production area in Hunan Province and an important commodity grain base in China [47].
Fourteen meteorological stations in the Xiang River Basin were used and the data were obtained from the National Meteorological Information Center (https://data.cma.cn/site/index.html).Daily precipitation, daily average temperature, max and min daily temperature, average pressure, average relative humidity, and daylight hours from 2015 to 2017 were provided by these stations.The Penman-Monteith equation was employed to calculate the potential evapotranspiration.

Reference Soil Moisture from China Land Soil Moisture Data Assimilation System
The China Land Soil Moisture Data Assimilation System (CLSMDAS) simulates and calculates soil moisture at different depths (0-5 cm, 0-10 cm, 10-40 cm, 40-100 cm, 100-200 cm) by using the Ensemble Kalman Filtering (EnKF) assimilation method and a land-surface model.The inputs of the land-surface model were comprised mainly of forcing data (i.e., FY-2C precipitation product and FY2C surface incident solar radiation product) and observation data (i.e., Advanced Microwave Scanning Radiometer for EOS (AMSR-E) soil moisture product and station observation soil moisture data).The output of the land-surface model was soil moisture data gridded over China with a spatial resolution of 10 km and temporal resolution of one day.The CLSMDAS soil moisture dataset was validated and the results showed that the assimilated soil moisture and in-situ observations were quite consistent in most parts of China's southern, eastern, and central regions [48][49][50].Moreover, the simulated high-quality soil moisture grid-point data can serve as important basic information for monitoring climate changes, including droughts.The CLSMDAS system presents reasonable temporal-spatial distribution of soil moisture and matches well with severe drought events in the southwest part of China [48][49][50].The soil moisture data derived from CLSMDAS has been widely applied for drought monitoring [48][49][50].The soil moisture dataset from CLSMDAS has been available since 19 January 2017.Therefore, the 0-5 cm soil moisture data from CLSMDAS in the Xiang River Basin during 1 April to 30 November 2017 were used as a reference, as there were no in situ soil moisture data available for this study.The 0-5 cm soil moisture dataset is provided online by the China Meteorological Data Sharing Service System: http://data.cma.cn/.

SMAP L2 Soil Moisture
The Soil Moisture Active Passive (SMAP) mission was launched by the National Aeronautics Space Agency (NASA) in January 2015.It provides a global mapping of soil and landscape freeze/thaw states, as well as the net ecosystem exchange of carbon, by coarser resolution (40 km) radiometer observations [51].The SMAP can provide the top 0-5 cm soil moisture data using the L-band radiometer.
The SMAP satellite mission was used to retrieve the global surface soil moisture information at different levels with a high-target accuracy (0.04 m 3 m −3 ) and a revisit of 1-3 days [35].There are four levels of the SMAP satellite mission.Level 1B and 1C have geolocated instrument measurements within 12 hours of acquisition; Level 2 products are geophysical retrievals of soil moisture on a fixed-Earth grid within 24 hours of acquisition; Level 3 products are actively used for soil moisture data within 50 hours of acquisition; and Level 4 are the products of surface and root zone soil moisture and carbon net ecosystem exchange within 7 and 14 days of acquisition, respectively [13].Therefore, Level 2 products were used in this study because of there daily temporal resolution.Specifically, the passive SMAP Level 2 L-band radiometer descending soil moisture product (SMAP L2 Radiometer Global Daily 36-km EASE-Grid soil moisture version 1, SMAP_L2_SM_P) was used to analyze drought in the Xiang River Basin.The dataset can be downloaded freely online: https://earthdata.nasa.gov/.

Methodology
The 0-5 cm soil moisture data from CLSMDAS was used as a reference dataset to evaluate the performance of 0-5 cm SMAP soil moisture.The atmospheric water deficits (AWD) calculated from meteorological data were used as reference datasets to evaluate the performance of SWDI.After the evaluation, the SMAP-derived soil water deficit index (SMAP_SWDI) was calculated to estimate the drought conditions in the Xiang River Basin.For comparison, AWD was also used to estimate the drought conditions in the basin.
The analysis was accomplished on both temporal and spatial scales.The steps are as follows: 1.
Evaluate the accuracy of SMAP soil moisture and SMAP_SWDI; 2.
Analyze the temporary and spatial drought conditions in the Xiang River Basin based on SMAP_SWDI; 3.
Calculate the AWD values with in situ meteorological data; 4.
Compare the performance of AWD and SMAP_SWDI on both temporal and spatial scales.

Soil Water Deficit Index (SWDI)
The agriculture drought is quantified by SWDI using the top 0-5 cm surface soil moisture.As described in previous studies, SWDI has shown a good performance in defining drought levels and severity; for example, the SWDI calculated using satellite soil moisture data showed very good agreement with the AWD, and the values of SWDI can adequately capture the drought dynamics, especially the intensity [12,52].The SWDI is calculated as follows: where θ is the time series of SMAP soil moisture (m 3 /m 3 ) and 0-5 cm soil moisture data from CLSMDAS, θ FC , θ WP and, θ AWC represents the field capacity, wilting point, and available water capacity, respectively.There are several ways to define θ FC and θ WP : (1) the 5th and 95th soil moisture data of the time series denote θ WP and θ FC ; (2) the soil moisture at a soil water potential of −33 kPa and −1500 kPa are considered equal to θ WP and θ FC ; and (3) θ WP and θ FC are calculated by basic soil physical characteristics such as the proportion of clay and sand via pedo-transfer functions [53].These three methods all show good performance in defining θ WP and θ FC , and the first method is the simplest way for calculation [53][54][55].Agricultural drought is linked with soil moisture availability during the growing season [53].Therefore, the 5th and 95th soil moisture data of the selected time series of the growing season (April to November) were used to obtain annual θ WP and θ FC .The daily SMAP_SWDI was computed based on the time series of every grid of this dataset and then the daily SWDI of grids were averaged for the whole Xiang River Basin.The SWDI was analyzed on a weekly timescale for it is the timescale used in the irrigation schedule [56], therefore, the average daily SWDI would transfer to a weekly SWDI.If the values of the SWDI are negative, it indicates the occurrence of drought in the Xiang River Basin.According to Martinez-Fernandez et al. [53], the values of the SWDI and the corresponding drought categories are shown in Table 1.

Percentage of Drought Weeks (PDW)
A drought week is defined when severe drought occurs with the values of weekly SMAP_SWDI lower than −5 [12].The percentage of drought weeks (PDWs) can be represented as the drought conditions of the Xiang River Basin.Percentage of drought weeks can be calculated as: where D represents the number of drought weeks and W represents the total study weeks, which were 35 in this study.

Atmospheric Water Deficit (AWD)
Atmospheric water deficit is the difference between precipitation (P) and potential evapotranspiration (ET 0 ).Atmospheric water deficit is claimed to be a suitable index to reflect the drought condition related to meteorological parameters [56].The AWD value is calculated by P minus ET 0 on a weekly scale: where i represents the week of the study period, P i and ET i denotes the sum of precipitation and sum of evapotranspiration (ET) of week i, respectively.The daily ET was computed using the Penman-Monteith equation.The Penman-Monteith equation is useful for computing ET based on temperature, humidity, wind speed, and solar radiation.It is rather stable over large areas, especially uniform surfaces [57].This index indicates drought when AWD is lower than 0. Extreme drought occurs when AWD is lower than −50 mm [53].

Evaluation Indices
The Pearson correlation coefficient (R-value), BIAS, ubRMSE, and the drought weeks' probability of detection (POD) were calculated to evaluate the accuracy and the reliability of the SMAP soil moisture and the performance of SMAP on drought monitoring.The three indices are defined as follows: where X i and Y i represent the in situ dataset and remote sensing dataset, respectively, and X and Y represent the mean value of these two datasets, respectively.In Equation ( 6), A represents the weeks that the droughts are both indicated by AWD and SMAP_SWDI, and C represents the weeks that AWD indicates drought while SMAP_SWDI indicates no drought.As long as in situ soil moisture data and precipitation or alternative datasets are available, the methodology in this study can be widely applied to evaluate the capability of SMAP on drought monitoring on other regions.

Evaluation of SMAP Soil Moisture
The grid scale and regional scale were both adopted to evaluate the accuracy of SMAP soil moisture.At the grid scale, 0-5 cm soil moisture data from CLSMDAS and 0-5 cm SMAP soil moisture were used for evaluation.Daily soil moisture on each grid obtained from SMAP and CLSMDAS between 1 April to 30 November 2017 were used for calculation of R-value.As for the regional scale, the average value of SMAP soil moisture was compared with the averaged CLSMDAS soil moisture by averaging all grids of these two datasets over the Xiang river basin.
At the grid scale, the daily CLSMDAS soil moisture and SMAP soil moisture were compared with each other on every grid via R-value (Figure 2a).In terms of R-value, around 70 percent of grids performed well (R-value > 0.6), 20 percent performed fairly (0.4 < R-value < 0.6), and just 10 percent performed poorly (R-value < 0.4).Higher consistency between CLSMDAS and SMAP soil moisture was mostly observed in the southwest part of the Xiang River Basin, which is located in a region with relatively low altitude compared to surrounding regions.Poor performance was mainly illustrated in the eastern basin with the highest altitude.These results suggest that topographic factors play an important role in the performance of SWDI for agricultural drought, and the results are consistent with Paredes-Trejo and Barbosa [9], who showed that agricultural drought indices obtained from satellite soil moisture tends to show poor performance in mountain regions.A possible reason is that the surface water in mountainous regions is limited due to the high potential for infiltration because of complex topography [34, 58,59].In terms of BIAS values (Figure 2b), about 85 percent of grids were between −20 % and 20%, 51 percent were between −10 % and 10%, and only 3.5 percent were over 30%.In addition, the SMAP soil moisture demonstrated good performance with ubRMSE values ranging from 0.03 m 3 •m −3 to 0.08 m 3 •m −3 .At the regional scale, the time series of averaged SMAP soil moisture along with the averaged CLSMDAS soil moisture on a daily scale are shown in Figure 3.The R-value, BIAS, and ubRMSE between these two averaged datasets (R-value = 0.62, BIAS = 6.04 %, ubRMSE = 0.033 m 3 •m −3 ) for this basin indicates that the SMAP soil moisture has a relatively good performance on capturing the features of CLSMDAS soil moisture.During May, the values of CLSMDAS soil moisture are slightly lower than SMAP soil moisture.Soil moisture from these two datasets show similar variations and features in middle June and July.However, significant difference occurs after August.The values of CLSMDAS soil moisture are higher than that of SMAP soil moisture, even though there are similar variations between these two datasets.A possible reason for the big gaps between these two datasets is that the L-band of SMAP cannot penetrate massive vegetation canopies [60].
Above all, SMAP soil moisture is satisfactory in terms of R-value, BIAS, and ubRMSE for both point scale and regional scale, with acceptable accuracy and relatively good performance in most of the Xiang River Basin.
with Paredes-Trejo and Barbosa [9], who showed that agricultural drought indices obtained from satellite soil moisture tends to show poor performance in mountain regions.A possible reason is that the surface water in mountainous regions is limited due to the high potential for infiltration because of complex topography [34, [58][59].In terms of BIAS values (Figure 2b), about 85 percent of grids were between −20 % and 20%, 51 percent were between −10 % and 10%, and only 3.5 percent were over 30%.In addition, the SMAP soil moisture demonstrated good performance with ubRMSE values ranging from 0.03 m 3 m −3 to 0.08 m 3 m −3 .
(a) (b) At the regional scale, the time series of averaged SMAP soil moisture along with the averaged CLSMDAS soil moisture on a daily scale are shown in Figure 3.The R-value, BIAS, and ubRMSE between these two averaged datasets (R-value = 0.62, BIAS = 6.04 %, ubRMSE = 0.033 m 3 m −3 ) for this basin indicates that the SMAP soil moisture has a relatively good performance on capturing the features of CLSMDAS soil moisture.During May, the values of CLSMDAS soil moisture are slightly lower than SMAP soil moisture.Soil moisture from these two datasets show similar variations and features in middle June and July.However, significant difference occurs after August.The values of CLSMDAS soil moisture are higher than that of SMAP soil moisture, even though there are similar variations between these two datasets.A possible reason for the big gaps between these two datasets is that the L-band of SMAP cannot penetrate massive vegetation canopies [60].Above all, SMAP soil moisture is satisfactory in terms of R-value, BIAS, and ubRMSE for both point scale and regional scale, with acceptable accuracy and relatively good performance in most of the Xiang River Basin.

Temporal and Spatial Analysis of SMAP_SWDI
At the temporal sale, the weekly time series of SMAP_SWDI from 2015 to 2017 are presented in Figure 4. Figure 4 illustrates that the weekly time series of SMAP_SWDI of three different years show a similar trend.The values of SMAP_SWDI were within −6 to 0 between April and mid-June, and then rapidly declined to −10 during July and October.After October, these values apparently increased to the range of −4 to −7. Figure 4 presents that the values of SMAP_SWDI were focused on

Temporal and Spatial Analysis of SMAP_SWDI
At the temporal sale, the weekly time series of SMAP_SWDI from 2015 to 2017 are presented in Figure 4. Figure 4 illustrates that the weekly time series of SMAP_SWDI of three different years show a similar trend.The values of SMAP_SWDI were within −6 to 0 between April and mid-June, and then rapidly declined to −10 during July and October.After October, these values apparently increased to the range of −4 to −7. Figure 4 presents that the values of SMAP_SWDI were focused on −4 to 0 during April and June of 2016 and 2017.According to Table 1, it indicates that there were no severe droughts in the Xiang River Basin at this stage.It coincides with Reference [61], which suggests that the Xiang River Basin may encounter less risk of drought due to heavy rainfalls during April to June.However, the values of SMAP_SWDI in April of 2015 were lower than that of 2016 and 2017, which suggests that the severity and magnitude of the agricultural drought in 2015 was more serious than that in 2016 and 2017.
There is an apparently dynamic trend where the values of SMAP_SWDI decreased from −1 to −10 during mid-June and August for all three years.Based on Table 1, the corresponding drought category changed from mild drought to severe drought in this period, and the condition of severe drought existed until the end of September.The possible reason for the mutation may be the high air temperature in the basin between June to September and high evaporation from the soil surface as well as high evapotranspiration of crops in the development stage and mid-season stage (Figure 5).Xiao et al. [62] also confirmed that the temperature gradually increased from April to September in the Xiang River Basin, and that there was no significant increase in rainfall in most meteorological stations from June to July.This demonstrates that the Xiang River Basin was threatened by the high risk of agricultural drought during April and September.In terms of interannual changes, the values of SMAP_SWDI in 2015 varied largely from week to week, and the mutation was half a month ahead compared to 2016 and 2017.Such findings indicate that the severity of agricultural drought in 2015 was stronger than 2016 and 2017.The values of SMAP_SWDI during the three years studied were stable between mid-May and June, August to October, and at the end of November.Therefore, the conditions of agricultural drought in mid-May, early August, and at the end of November in the Xiang River Basin were more representative than other stages.Thus, the spatial distributions of SMAP_SWDI in mid-May, early August, and at the end of November from 2015 to 2017 in this basin were selected to be compared.
As shown in Figure 6, during mid-May, the severity of droughts gradually decreased from 2015 to 2017.Severe droughts (−10 < SWDI < −5) were only detected in the southern and northeastern parts of the basin in 2015, and less moderate droughts occurred in 2016 and 2017.
Early August basically represents the driest period in the Xiang River Basin.Early August of 2016 suffered the most severe and wide range of droughts, with the driest area mainly in the northern part of the basin.However, the same category of drought as 2016 mainly occurred in the southern part of the basin in 2015.During early August of 2017, the severity of agricultural drought in the central region was relieved compared to 2016, but the severe drought was shifted from the north to the eastern part of the basin.The severity of agricultural drought was reduced in both the northeast and southwestern part of the basin in early August of 2017, but there remains a small area with severe drought in the middle basin.As for the period at the end of November during 2015 and 2017, severe drought was mainly concentrated in the northern part of the basin.Drought conditions in 2015 were apparently less severe than that in 2016 and 2017.The distribution of different categories of droughts in 2016 is similar to that of 2017, with severe droughts occurring in the northwest and moderate droughts occurring in the eastern part of the basin.
and southwestern part of the basin in early August of 2017, but there remains a small area with severe drought in the middle basin.As for the period at the end of November during 2015 and 2017, severe drought was mainly concentrated in the northern part of the basin.Drought conditions in 2015 were apparently less severe than that in 2016 and 2017.The distribution of different categories of droughts in 2016 is similar to that of 2017, with severe droughts occurring in the northwest and moderate droughts occurring in the eastern part of the basin.

Temporal and Spatial Analysis of PDWs
The annual PDWs of the Xiang River Basin showed a gradual downward trend from 2015 to 2017 (from 65% to 57%) in Figure S1.This downward trend is consistent with Zou and Ren [63] who showed that most river basins in Southern China demonstrated a slightly downward trend in drought in recent years.
After calculating the PDWs of each grid, the spatial change in PDWs can be analyzed from 2015 to 2017 (Figure 7). Figure 7 illustrates that PDWs followed an increasing trend from the northeast to the southwest in 2015 and 2016.Figure 7a shows that the severe conditions of drought weeks mainly occurred in the northwest and southwest of the basin.As for 2016 (Figure 7b), the severity of drought weeks was alleviated in the western region compared to that of 2015, while the condition of drought weeks in the east was more serious than that in 2015.In 2017, the PDW values were reduced over the entire basin, especially the western and southern parts, but there remained some parts of the basin that had high PDW values (Figure 7c).The severe condition of drought weeks in 2017 was mainly in small parts of the northwest and northeast, but the extent of the drought weeks was significantly reduced compared to the previous two years.
From the perspective of the entire basin, the northeastern part of the basin was in a state of fewer drought weeks over the past three years.This result is consistent with the findings of Zhang et al. [7], who claims that the values of Pan evaporation in the Xiang River Basin declined from the southwest to northeast, indicating that the northeast part may not suffer high PDWs compared with the other parts of the basin.Significant changes can be seen in the western and southern parts of the basin, which show a process of progressively weakening PDWs upstream of the basin.
From the perspective of the entire basin, the northeastern part of the basin was in a state of fewer drought weeks over the past three years.This result is consistent with the findings of Zhang et al. [7], who claims that the values of Pan evaporation in the Xiang River Basin declined from the southwest to northeast, indicating that the northeast part may not suffer high PDWs compared with the other parts of the basin.Significant changes can be seen in the western and southern parts of the basin, which show a process of progressively weakening PDWs upstream of the basin.The PDW has four categories: (1) 0~25% (i.e., it represents less than nine weeks of drought events); (2) 25~50%; (3) 50~75%; (4) 75~100% (i.e., it represents more than 26 weeks of drought events).A higher percentage of grids that are located in the intervals of 50%~75% and 75%~100% indicates a higher percentage of drought events [12].The PDW in four intervals were quantified in the Xiang The PDW has four categories: (1) 0~25% (i.e., it represents less than nine weeks of drought events); (2) 25~50%; (3) 50~75%; (4) 75~100% (i.e., it represents more than 26 weeks of drought events).A higher percentage of grids that are located in the intervals of 50%~75% and 75%~100% indicates a higher percentage of drought events [12].The PDW in four intervals were quantified in the Xiang River Basin.It can be seen from Figure 8 that the percentage of grids in the four intervals present different values of PDWs in each year.Most of the PDW values were between 50% and 75%.In 2015, except for the interval of 50-75%, the percentage of grids in other intervals stayed in a smaller range (Figure 8a).In 2016 and 2017, there are few cases where the proportion of PDWs was more than 75% or less than 50% (Figure 8b,c).These results suggest that the severity of drought weeks have eased in the past three years.However, since most of the PDWs were still in the interval of 50% to 75%, attention should be paid to drought monitoring in the Xiang River Basin in the coming years.
(c) The PDW has four categories: (1) 0~25% (i.e., it represents less than nine weeks of drought events); (2) 25~50%; (3) 50~75%; (4) 75~100% (i.e., it represents more than 26 weeks of drought events).A higher percentage of grids that are located in the intervals of 50%~75% and 75%~100% indicates a higher percentage of drought events [12].The PDW in four intervals were quantified in the Xiang River Basin.It can be seen from Figure 8 that the percentage of grids in the four intervals present different values of PDWs in each year.Most of the PDW values were between 50% and 75%.In 2015, except for the interval of 50-75%, the percentage of grids in other intervals stayed in a smaller range (Figure 8a).In 2016 and 2017, there are few cases where the proportion of PDWs was more than 75% or less than 50% (Figure 8b and 8c).These results suggest that the severity of drought weeks have eased in the past three years.However, since most of the PDWs were still in the interval of 50% to 75%, attention should be paid to drought monitoring in the Xiang River Basin in the coming years.

Comparison between SMAP_SWDI and AWD
To verify the reliability of SMAP_SWDI, a comparison was conducted with the agricultural drought index AWD.On the one hand, the AWD is a widely applied drought index using atmospheric data.On the other hand, the AWD is a suitable tool to capture drought dynamics related to soil water storage [58].The AWD is used to evaluate the efficiency of soil moisture from SMAP on agricultural drought monitoring.To estimate the fitting relationship between SMAP_SWDI and AWD at the point scale, 13 stations were chosen from the basin.For each station, the values of R,

Comparison between SMAP_SWDI and AWD
To verify the reliability of SMAP_SWDI, a comparison was conducted with the agricultural drought index AWD.On the one hand, the AWD is a widely applied drought index using atmospheric data.On the other hand, the AWD is a suitable tool to capture drought dynamics related to soil water storage [58].The AWD is used to evaluate the efficiency of soil moisture from SMAP on agricultural drought monitoring.To estimate the fitting relationship between SMAP_SWDI and AWD at the point scale, 13 stations were chosen from the basin.For each station, the values of R, POD, and PDWs of AWD and SMAP_SWDI are shown in Table 2.The values of R between AWD and SMAP_SWDI varied from 0.42 to 0.68 (average R-value = 0.6).Around six stations performed well (R-value > 0.6) and seven stations performed fairly (0.4 < R-value < 0.6), indicating that there is a relatively high agreement between these two indices.The relatively good correlation between SMAP_SWDI and AWD is similar to the results obtained in several previous studies [12,57].The highest value of R was observed at Yongzhou Station (R-value = 0.68), and its POD was 0.79.The lowest value of R belonged to Yuanjiang (R-value = 0.43) with a very low value of POD (0.23).Almost every station performed reasonably on POD, varying from 0.7 to 0.93, except Chenzhou and Yuanjiang.The best performance was observed at Nanyue Station with the value of 0.93.The results indicate that SMAP_SWDI shows a good capture of AWD in the detection of drought.
In addition, a drought week is defined when SMAP_SWDI is lower than −5 or AWD is lower than 0 mm [12].It can be observed that the PDWs obtained by SWDI were higher than that obtained by AWD at most stations.This is similar to the conclusions obtained from a previous study [55] From Table 2, about 11 stations show high PDWs with SMAP_SWDI.Figure 9 shows the distribution of the numbers of PDWs of AWD and SMAP_SWDI in four intervals.We can see that the proportion of PDWs obtained by SMAP_SWDI was mainly between 50% and 75%, while that of AWD was between 25% and 50%.These indicate that SMAP may overestimate PDWs compared to in situ AWD.
The weekly evolution of AWD and SMAP_SWDI at randomly selected stations is shown in Figure 10.Stations Malingpo (Figure 10a), Pingjiang (Figure 10b), Changsha (Figure 10c), and Changning (Figure 10d) are selected as examples, with the values of R varying from 0.45 to 0.6.Figure 10 illustrates that SMAP_SWDI and AWD displayed similar seasonal cycles.This result is due to the top 5 cm surface soil moisture being directly related to the atmospheric water content [64].
Among these four sites, Changning (Figure 10d) showed relatively worse performance compared with the other three stations.From Figure 10a-c, the values of weekly SMAP_SWDI and AWD were much higher when the basin experienced weekly heavy rainfall.The values of AWD were more sensitive to heavy rainfall according to Figure 10b,c.The changes of values of SMAP_SWDI and AWD were consistent with rainfall dynamics, especially among the weeks of precipitation peaks and droughts.POD, and PDWs of AWD and SMAP_SWDI are shown in Table 2.The values of R between AWD and SMAP_SWDI varied from 0.42 to 0.68 (average R-value = 0.6).Around six stations performed well (R-value > 0.6) and seven stations performed fairly (0.4 < R-value < 0.6), indicating that there is a relatively high agreement between these two indices.The relatively good correlation between SMAP_SWDI and AWD is similar to the results obtained in several previous studies [12,57] In addition, a drought week is defined when SMAP_SWDI is lower than −5 or AWD is lower than 0 mm [12].It can be observed that the PDWs obtained by SWDI were higher than that obtained by AWD at most stations.This is similar to the conclusions obtained from a previous study [55] From Table 2, about 11 stations show high PDWs with SMAP_SWDI.Figure 9 shows the distribution of the numbers of PDWs of AWD and SMAP_SWDI in four intervals.We can see that the proportion of PDWs obtained by SMAP_SWDI was mainly between 50% and 75%, while that of AWD was between 25% and 50%.These indicate that SMAP may overestimate PDWs compared to in situ AWD.
The weekly evolution of AWD and SMAP_SWDI at randomly selected stations is shown in Figure 10.Stations Malingpo (Figure 10a), Pingjiang (Figure 10b), Changsha (Figure 10c), and Changning (Figure 10d) are selected as examples, with the values of R varying from 0.45 to 0.6.Figure 10 illustrates that SMAP_SWDI and AWD displayed similar seasonal cycles.This result is due to the top 5 cm surface soil moisture being directly related to the atmospheric water content [64].
Among these four sites, Changning (Figure 10d) showed relatively worse performance compared with the other three stations.From Figure 10a, 10b, and 10c, the values of weekly SMAP_SWDI and AWD were much higher when the basin experienced weekly heavy rainfall.The values of AWD were more sensitive to heavy rainfall according to Figure 10b and 10c.The changes of values of SMAP_SWDI and AWD were consistent with rainfall dynamics, especially among the weeks of precipitation peaks and droughts.

Conclusions
In this study, the accuracy of satellite-based soil moisture from SMAP was evaluated with the soil moisture from CLSMDAS used as the reference dataset.In addition, the potential applications of SMAP for drought monitoring in the Xiang River Basin was investigated with a specific agricultural drought index (SWDI).A temporal and spatial analysis was conducted by comparing the SMAP_SWDI and the AWD based on meteorological variables.These two indices were estimated weekly and compared in the growing season (April to November) from 2015 to 2017.The main conclusions are as follows: 1.
The SMAP soil moisture had acceptable accuracy and relatively good performance in most of Xiang River Basin.In terms of R-values between the soil moisture datasets obtained from CLSMDAS and SMAP, around 70% performed well and only 10% performs poorly at the grid sale.At the regional scale, the SMAP soil moisture captured the features of CLSMDAS soil moisture with the R-value 0.62.

2.
Severe droughts mainly occurred from mid-June to the end of September.The severity of droughts gradually decreased in mid-May from 2015 to 2016.Severe droughts were only detected in the south and northeast of the Xiang River Basin in 2015.Early August of 2016 suffered the most severe and wide range of droughts, with the driest area mainly concentrated in the northern part of the basin.Severe drought was mainly concentrated in the northern part of the Xiang River Basin at the end of November during 2015 and 2017.

3.
The SMAP showed relatively good performance in drought monitoring in the Xiang River Basin with a high Pearson correlation coefficient (mean value equals to 0.6) and high drought weeks probability of detection (vary from 0.7 to 0.9) between SWDI and AWD.

Figure 1 .
Figure 1.Xiang River Basin and distribution of the meteorological stations.

Figure 1 .
Figure 1.Xiang River Basin and distribution of the meteorological stations.

Figure 2 .
Figure 2. R-values (a) and BIAS (b) values between Soil Moisture Active Passive (SMAP) soil moisture and China Land Soil Moisture Data Assimilation System (CLSMDAS) soil moisture datasets.

Figure 2 .
Figure 2. R-values (a) and BIAS (b) values between Soil Moisture Active Passive (SMAP) soil moisture and China Land Soil Moisture Data Assimilation System (CLSMDAS) soil moisture datasets.

Figure 3 .
Figure 3. Regional comparison between SMAP soil moisture and CLSMDAS soil moisture.

Figure 3 .
Figure 3. Regional comparison between SMAP soil moisture and CLSMDAS soil moisture.
Remote Sens. 2019, 11, x FOR PEER REVIEW 9 of 21 compared to 2016 and 2017.Such findings indicate that the severity of agricultural drought in 2015 was stronger than 2016 and 2017.

Figure 4 .
Figure 4. Weekly time series of SMAP_SWDI in the Xiang River Basin.Figure 4. Weekly time series of SMAP_SWDI in the Xiang River Basin.

Figure 4 .
Figure 4. Weekly time series of SMAP_SWDI in the Xiang River Basin.Figure 4. Weekly time series of SMAP_SWDI in the Xiang River Basin.

Figure 4 .
Figure 4. Weekly time series of SMAP_SWDI in the Xiang River Basin.

Figure 5 .
Figure 5. Average monthly temperature and potential evapotranspiration from April to November between 2015 and 2017.

Figure 6 .
Figure 6.Spatial distribution of SMAP_SWDI during mid-May, early August, and at the end of November in the Xiang River Basin from 2015 to 2017.

Figure 6 .
Figure 6.Spatial distribution of SMAP_SWDI during mid-May, early August, and at the end of November in the Xiang River Basin from 2015 to 2017.
. The highest value of R was observed at Yongzhou Station (R-value = 0.68), and its POD was 0.79.The lowest value of R belonged to Yuanjiang (R-value = 0.43) with a very low value of POD (0.23).Almost every station performed reasonably on POD, varying from 0.7 to 0.93, except Chenzhou and Yuanjiang.The best performance was observed at Nanyue Station with the value of 0.93.The results indicate that SMAP_SWDI shows a good capture of AWD in the detection of drought.

Figure 9 .
Figure 9. Distribution of the PDWs of atmospheric water deficit (AWD) and SMAP_SWDI in four intervals.

Figure 9 .
Figure 9. Distribution of the PDWs of atmospheric water deficit (AWD) and SMAP_SWDI in four intervals.

Table 2 .
Values of CC, POD, and PDWs of AWD and SMAP_SWDI in different stations.