Characterization and Evaluation of MODIS-Derived Crop Water Stress Index (CWSI) for Monitoring Drought from 2001 to 2017 over Inner Mongolia

: Inner Mongolia is one of the main green production bases of agricultural and animal husbandry products. Due to factors such as natural geographical location, drought occurs frequently in Inner Mongolia. Based on the MOD16 product and the method of crop water stress index (CWSI) combined with multi-year precipitation and temperature data, the spatial and temporal distribution characteristics and major inﬂuencing factors of drought in Inner Mongolia from 2001 to 2017 were analyzed. In order to provide effective scientiﬁc basis for drought control and drought resistance in Inner Mongolia for decision. The results showed that: (1) during 2001–2017, the average annual CWSI in Inner Mongolia had a strong spatial heterogeneity, which showed a trend of gradual increase from northeast to southwest. The annual average CWSI was 0.7787 and showed a ﬂuctuating downward trend for Inner Mongolia. (2) The CWSI of every 8d during one year in Inner Mongolia showed the double-peak trend, reaching its maximum of 0.9043 in the 121st day. In addition, the average CWSI of every 8d was 0.6749. (3) In Inner Mongolia, the average CWSI of different land-use types showed little difference and ranged from small to large: woodland (0.5954) < cropland (0.7733) < built-up land (0.8126) < grassland (0.8147) < unused land (0.8392). (4) The average correlation coefﬁcients between CWSI and precipitation, temperature respectively were − 0.53 and 0.18, which indicated that CWSI was highly correlated with precipitation in Inner Mongolia. and the annual average CWSI values are below 0.7725 and there is no drought. The ETs of Hunshandake sandy land at the southern end of Xilinguole grassland in central Inner Mongolia and Yinshan Mountain range in Hetao plain of the Yellow River are relatively low, while the PET and CWSI are relatively high and the multi-year average CWSI is above 0.8985, which is characterized by severe drought and extreme drought. The results of this study may provide some practical reference for the aim to prevent drought and maintain grassland productivity in a warming climate.


Introduction
Drought is a recurring natural disaster that frequently occurs and lasts for a long time and has an important impact on economic development and ecological environment construction, especially on agricultural production [1,2]. Drought is generally divided into four types: meteorological drought, agricultural drought, hydrological drought and socio-economic drought [3]. At present, the main research fields are meteorological drought and agricultural drought but agricultural drought is more complex. Daily, monthly and annual weather and climate disturbances in agricultural region have a significant impact on local and regional scales [4].
With the development of remote sensing technology, many scholars had used satellite land observation data to monitor the dynamic changes of various land surface [5][6][7]. Compared with traditional observation methods, remote sensing monitoring provides continuous temporal and spatial monitoring of vegetation cover, soil moisture level and drought severity [8]. Therefore, drought indices derived from remote sensing-based measurements have been widely applied to monitor agricultural drought conditions, including Sustainability 2021, 13, 916 2 of 17 temperature vegetation drought Index (TVDI), anomaly vegetation index (AVI), Vegetation condition index, VCI), crop water stress index (CWSI), perpendicular drought index (PDI), et al. [9][10][11]. The crop water stress index (CWSI) is based on the surface energy balance theory, which has high measurement accuracy and clear physical definition [12]. It is an important factor that is superior to other drought indices. In addition, the CWSI also considers the effects of wind speed, water vapor pressure, sunshine hours and temperature [11]. However, its calculation is complex and requires the participation of conventional meteorological data and ground observation data, which limits the application of CWSI [13]. In 2011, the National Aeronautics and Space Administration (NASA) research team made significant advance in evapotranspiration (ET) retrieval from moderate MODIS data based on the Penman-Monteith (P-M) equation [14,15] and the global MOD16 data were released with high temporal-spatial resolution and free access [16]. It provides a new idea for the calculation of CWSI. At present, many scholars have validated MOD16 products, which shows that the data has a wide range of applications and high reliability. In addition, it had been widely applied in some field such as water use efficiency, monitoring of wet and dry patterns, estimating actual evapotranspiration and water deficit footprint and some scholars made significant progress [17][18][19]. Hence, MOD16 data has more extensive application value in drought monitoring.
The Inner Mongolia is one of the major producing areas of agricultural and animal husbandry products in China. However, the frequent occurrences of droughts pose a significant threat to agricultural and animal husbandry production in the region. Therefore, in order to perform drought monitoring and research for the Inner Mongolia, the development of the appropriate drought monitoring model to analyze the spatial-temporal evolution and severity of droughts has great importance to agricultural, animal husbandry production and economic development. In this study, CWSI was used to analyze the spatial-temporal changes of drought in Inner Mongolia from 2001 to 2017 and the major influencing factors of drought in Inner Mongolia were discussed. The purpose of this study was to evaluate CWSI in the Inner Mongolia area of drought monitoring effectiveness, provide scientific basis for drought control and drought resistance decision.

Study Area
The Inner Mongolia (37 • 24 N-53 • 23 N, 97 • 12 E-126 • 04 E) is located in the northern frontier of China. With 1,183,000 km 2 in size, it is about 13% of the total area of the China. It spans the northeast, north and northwest regions of the China, adjacent to eight provinces in China and adjacent to Mongolia and Russia ( Figure 1). The climate in the area, which gradually transits from humid and semi-humid areas to semi-arid and arid areas from east to west, is mainly temperate continental monsoon climate, with annual precipitation between 46 mm to 489 mm. However, annual average evaporation is between 942 and 4138 mm and the degree of drought is more serious [20]. Geomorphological types are complex and diverse, mainly in plateau, with an average elevation of more than 1000 m. Due to the obvious spatial differences in climate, the vegetation types from east to West are forest, grassland and desert in turn. Because of the special geographical location, topography and climatic conditions, the ecological environment in this area is extremely fragile, grassland degradation, serious soil erosion and unbalanced distribution of water resources [21][22][23][24][25]. Drought is one of the most serious natural disasters in the study area.

Data Sources and Processing
The data used in this paper include land use data, meteorological data, remote sensing data and relative soil moisture data (Table 1).

Data Sources and Processing
The data used in this paper include land use data, meteorological data, remote sensing data and relative soil moisture data (Table 1). (1) Land use data The land use data used in this study are the standard grid data products with 1 km resolution provided by Resource and Environment Science Data Center of Chinese Academy Using ENVI and ArcGIS, the land cover data in China were cut and reclassified and classified into six first-class land types: cropland, woodland, grassland, water body, built-up land and unused land ( Figure 1). The spatial resolution of the data is consistent with that of MOD16.
(2) Meteorological data Meteorological data (temperature, precipitation) of 44 weather stations in Inner Mongolia from 2001 to 2017 were provided by the China meteorological data sharing network. The basic information and the specific geographical location of 44 meteorological stations is shown in Table 2 and Figure 1. Based on the longitude and latitude data of weather stations in this study area, the spatial distribution of meteorological elements was realized by Kriging interpolation under the ArcGIS and the grid data of meteorological elements were obtained and the spatial resolution was consistent with MOD16 data. The time series of annual average precipitation and temperature in this study area were calculated. (3) Remote sensing data MOD16 products synthesized in 8 days from 2001 to 2017 were provided by NASA. The satellite orbital number is h25v03/h25v04/h25v05/h26v03/h26v04/h26v05/h27v04, including actual evapotranspiration (MOD16ET) and potential evapotranspiration (MOD-16PET) data. With the help of MRT software provided by NASA, the original HDF format was transformed into GeoTiff format and the SIN projection was transformed into the WGS~1984/Geographic longitude and latitude coordinate system, which was mosaic and clipping. According to the data usage instructions provided by the website, the invalid values in the data were eliminated and the real values were restored. The ET/PET data of 8D were weighted averagely to obtain the ET and PET values of each year, month in this study area.
(4) Relative soil moisture data The data set of crop growth and development in China provided by China Meteorological Data Network is based on the ten-day monthly Agrometeorological data reported by China Agrometeorological Stations since 1991 (http://data.cma.cn/data/). There were 31 observatories in Inner Mongolia with sampling depths of 10, 20, 50, 70 and 100 cm and the observation time was every ten days. Due to the influence of weather, man-made, instruments and other factors, some stations lack data integrity. Therefore, based on the comprehensive consideration of data integrity, month continuity and 20 cm relative soil moisture as the drought classification standard of the National Meteorological Administration, 16 stations with complete data were finally screened out, including 2003-2009 and April-September 2011 (After thawing and before freezing) 10-day 20 cm relative soil moisture data, according to the climatic conditions in this study area, the abnormal values were eliminated and the 10-day data were weighted averagely to synthesize monthly relative soil moisture data. The basic information and the specific 16 agrometeorological stations is shown in Table 3 and Figure 1. Crop water stress index (CWSI), based on the principle of surface energy balance, has high measurement accuracy and clear physical significance, which is one reason why it is superior to other drought indexes. In addition, CWSI does not need to consider that the study area has all land use types, moreover, factors such as the vegetation cover of the underlying surface, wind speed, water vapor pressure, sunshine duration and air temperature are fully considered [25,26]. Jackson et al. studied CWSI and defined CWSI as [27,28]: where ET is the actual evapotranspiration, PET is the potential evapotranspiration, CWSI is between 0 and 1, the larger the value is, the more arid and water-deficient the region is. Conversely, the smaller the value, the wetter it is.

Linear Tendency Estimation
The temporal linear propensity ratio of CWSI for each pixel from 2001 to 2017 is calculated by using linear tendency estimation. The calculation formula is as follows: where slope is linear tendency rate; t i = 2001, 2002, . . . . . . , 2017; n is the total length of the annual sequence (n = 17); CWSI i is the cumulative value of annual CWSI of the ith year; When slope > 0, with the increase of time, CWSI tends to increase, whereas CWSI tends to decrease. According to the significance test and linear tendency estimation of correlation coefficient, the change trend of CWSI was divided into seven grades (Table 4).

Correlation Coefficient Method
The correlation between CWSI of each pixel in Inner Mongolia and precipitation, temperature is calculated. The simple correlation coefficient is used to express the correlation. The formula is as follows: where r is a simple correlation coefficient between the two variables; CWSI is the average value of CWSI from 2001 to 2017; x i is the precipitation or temperature value of the ith year; x is the average value of precipitation or temperature from 2001 to 2017; other quantities have the same meaning as before. Moreover, the proportion of area, correlation between CWSI and precipitation, temperature positively and negatively, was calculated by ENVI statistical analysis function.

Drought Suitability Evaluation and Classification Standard
China Meteorological Administration takes the relative soil moisture at 20 cm as the grading standard for drought [29]. Therefore, scholars tested the accuracy of drought monitoring by analyzing the correlation between CWSI and relative soil moisture at 20 cm. While the degree of drought was divided by the linear regression equation between CWSI and relative soil moisture at 20 cm [8,11].

Drought Monitoring Results Test
In order to ensure the accuracy of CWSI drought monitoring results in this study area, the CWSI calculated based on MOD16 products was tested. The results showed that relative soil moisture at 20 cm can effectively reflect the accuracy of CWSI drought detection. Therefore, CWSI is used as a drought monitoring index. Figure 2 shows the correlation between monthly CWSI and relative soil moisture at 20 cm in 16 stations from 2003 to 2009 and from April to September 2011. The results showed that the correlation coefficient between CWSI and relative soil moisture at 20 cm is −0.86, which indicated that the CWSI calculated by MOD16 has good validation accuracy and applicability in Inner Mongolia and can be used for drought monitoring in this area. In addition, the results show that the linear regression equation between CWSI and relative soil moisture (RSM) at 20 cm is y =−0.0063x + 1.1505, where y is CWSI and X is soil relative humidity. Finally, the classification standard for drought is determined. The classification is as shown in Table 5.

Interannual Spatio-Temporal Distribution Characteristics of CWSI
According to the interannual variation of evapotranspiration and crop water shortage index in Inner Mongolia (Figure 3), ET and PET fluctuated little. Among them, the range of ET fluctuation is 211-282 mm, the average ET for many years is 246 mm. In addition, the range of PET fluctuation is 1024-1210 mm, the average PET for many years is 1110 mm. The range of the difference between PET and ET is 742-981 mm and the average difference for many years is 864 mm. ET is increasing slowly, PET is decreasing slowly, DET is decreasing obviously, CWSI is decreasing obviously, indicating that the difference between PET and ET is increasing slowly. The trend of drought in Mongolia has slowed down. The average CWSI for many years is 0.7787. The maximum value appears in 2001 (0.8189) and the minimum value in 2016 (0.7243). As a whole, Inner Mongolia is in a mild state of drought. In addition, the results show that the linear regression equation between CWSI and relative soil moisture (RSM) at 20 cm is y = −0.0063x + 1.1505, where y is CWSI and X is soil relative humidity. Finally, the classification standard for drought is determined. The classification is as shown in Table 5.

Interannual Spatio-Temporal Distribution Characteristics of CWSI
According to the interannual variation of evapotranspiration and crop water shortage index in Inner Mongolia (Figure 3), ET and PET fluctuated little. Among them, the range of ET fluctuation is 211-282 mm, the average ET for many years is 246 mm. In addition, the range of PET fluctuation is 1024-1210 mm, the average PET for many years is 1110 mm. The range of the difference between PET and ET is 742-981 mm and the average difference for many years is 864 mm. ET is increasing slowly, PET is decreasing slowly, DET is decreasing obviously, CWSI is decreasing obviously, indicating that the difference between PET and ET is increasing slowly. The trend of drought in Mongolia has slowed down. The average CWSI for many years is 0.7787. The maximum value appears in 2001 (0.8189) and the minimum value in 2016 (0.7243). As a whole, Inner Mongolia is in a mild state of drought.    (Figure 5a-c). ET shows a strip-like change in the east, the middle and the West on the whole, while PET and CWSI shows a trend of change in the west, the middle and the East on the whole. The trend of CWSI and ET is opposite, which is consistent with that of PET. The annual average ET fluctuation ranges from 39 mm to 677 mm, PET fluctuation ranges from 375 mm to 1874 mm and CWSI fluctuation ranges from 0.3312 to 0.9991. The ET in the Great Hinggan Mountains and the foreststeppe ecotone in eastern Inner Mongolia are significantly higher than those in other regions, while the PET and CWSI values are significantly lower than those in other regions     The variation trend of CWSI is opposite with that of ET but it is same as PET can be explained by related theory analysis that evapotranspiration is complementa the definition of CWSI [30]. When the underlying surface is wet enough, ET is eq PET and CWSI = 0. When the water content is insufficient, the ET decreases and the action between the land surface and the atmosphere results in the increase of PET a increase of CWSI; when there is no water supply, there is basically no evapotranspi and chamaephyte will wither or die and CWSI equals to 1. In the eastern part of Mongolia, the vegetation coverage and precipitation are higher in the Da Hinggan M tains and the interlaced areas of woodland and grassland. The lower average annua perature makes the ET value and PET value larger, so the CWSI is smaller in this are Hunshandak sandy land at the southern end of the Xilin Gol Prairie in the centra and the Yinshan Mountains in the Hetao plain of the Yellow River has higher solar tion, longer sunshine hours, less precipitation, higher temperature and insufficient supply on the underlying surface, resulting in a lower ET value and a larger PET so CWSI is larger here.
The change trends of ET, PET and CWSI from 2001 to 2017 were calculated b mula (2), the results are shown in Figure 6a-c. According to Table.4, the change tr CWSI is divided and the results are shown in Figure.6d. The Figure 6a shows tha 2001 to 2017, the annual spatial change rate of ET ranged from −31.39 to 37.92 mm/a an average annual change rate of 3.75 mm/a. Among the annual rate of change, 55.8 the area shows a significant change. Among them, the area with a negative change only accounts for 1.62% and the area with a positive change trend is as high as 98.38 the area larger than 8 mm/a accounts for 7.23%, mainly distributed in the woodlan The variation trend of CWSI is opposite with that of ET but it is same as PET. This can be explained by related theory analysis that evapotranspiration is complementary and the definition of CWSI [30]. When the underlying surface is wet enough, ET is equal to PET and CWSI = 0. When the water content is insufficient, the ET decreases and the interaction between the land surface and the atmosphere results in the increase of PET and the increase of CWSI; when there is no water supply, there is basically no evapotranspiration and chamaephyte will wither or die and CWSI equals to 1. In the eastern part of Inner Mongolia, the vegetation coverage and precipitation are higher in the Da Hinggan Mountains and the interlaced areas of woodland and grassland. The lower average annual temperature makes the ET value and PET value larger, so the CWSI is smaller in this area. The Hunshandak sandy land at the southern end of the Xilin Gol Prairie in the central part and the Yinshan Mountains in the Hetao plain of the Yellow River has higher solar radiation, longer sunshine hours, less precipitation, higher temperature and insufficient water supply on the underlying surface, resulting in a lower ET value and a larger PET value, so CWSI is larger here.
The change trends of ET, PET and CWSI from 2001 to 2017 were calculated by Formula (2), the results are shown in Figure 6a-c. According to Table 4, the change trend of CWSI is divided and the results are shown in Figure 6d. The Figure 6a shows that from 2001 to 2017, the annual spatial change rate of ET ranged from −31.39 to 37.92 mm/a, with an average annual change rate of 3.75 mm/a. Among the annual rate of change, 55.84% of the area shows a significant change. Among them, the area with a negative change trend only accounts for 1.62% and the area with a positive change trend is as high as 98.38% and the area larger than 8 mm/a accounts for 7.23%, mainly distributed in the woodland and grassland ecotone and the area with relatively dense farmland (Figure 1). Most of the areas in Hulunbeir failed to pass the significance test, which indicates that the trend of change is obscure. The annual change rate of PET space is between −104.08 and 86.11 mm/a, the a change rate is −3.95 mm/a and there is a significant change in the area of 33.08% annual rate of change. The area with positive change trend only accounts for 1.99% area with a negative change trend is as high as 98.01% and the area smaller than −10 accounts for 85.57%. It is mainly distributed in the low hilly area in the southeast o Gol League and the meadow steppe area in the east, while most other areas do no the significance test and the change trend is obscure (Figure 6b).
The annual change rate of CWSI space is between −0.0296 and 0.0195, the a change rate is −0.0047 and there is a significant change in the area of 79.81% of the a rate of change. The area with positive change trend only accounts for 0.08%. The area a negative change trend is as high as 99.92% and the area smaller than 0.0089 accoun 11.65%. The range of variation gradually decreases from east to west and the dr mitigation trend is obvious in the meadow steppe in the east of Xilin Gol League b in the Yinshan Mountain range and Hunshandak sandy land (Figure 6c-d). Figure 7 shows the annual average change of ET, PET and CWSI in 8 d Mon from 2001 to 2017 and as shown in the graph, we can see that the yearly variation of Inner Mongolia area is a single peak curve and the maximum occurs on day 217(at mer)(14.77 mm) ; the yearly change of PET also shows a single peak pattern but d day 145 to 217 (late spring to summer), the PET fluctuates slightly and the large appeared on the 145th day (late spring) (49.75 mm). For specific explanation is that shows no drought on day 1 to72 (winter, early spring), day 169 to 248 (summer, autumn) and day 289 to 365 (366 days a year) (late autumn, winter), which may be c The annual change rate of PET space is between −104.08 and 86.11 mm/a, the annual change rate is −3.95 mm/a and there is a significant change in the area of 33.08% of the annual rate of change. The area with positive change trend only accounts for 1.99%. The area with a negative change trend is as high as 98.01% and the area smaller than −10 mm/a accounts for 85.57%. It is mainly distributed in the low hilly area in the southeast of Xilin Gol League and the meadow steppe area in the east, while most other areas do not pass the significance test and the change trend is obscure (Figure 6b).

Spatial-Temporal Distribution Characteristics of CWSI in the Year
The annual change rate of CWSI space is between −0.0296 and 0.0195, the annual change rate is −0.0047 and there is a significant change in the area of 79.81% of the annual rate of change. The area with positive change trend only accounts for 0.08%. The area with a negative change trend is as high as 99.92% and the area smaller than 0.0089 accounts for 11.65%. The range of variation gradually decreases from east to west and the drought mitigation trend is obvious in the meadow steppe in the east of Xilin Gol League but not in the Yinshan Mountain range and Hunshandak sandy land (Figure 6c-d). Figure 7 shows the annual average change of ET, PET and CWSI in 8 d Mongolia from 2001 to 2017 and as shown in the graph, we can see that the yearly variation of ET in Inner Mongolia area is a single peak curve and the maximum occurs on day 217 (at summer) (14.77 mm); the yearly change of PET also shows a single peak pattern but during day 145 to 217 (late spring to summer), the PET fluctuates slightly and the large value appeared on the 145th day (late spring) (49.75 mm). For specific explanation is that CWSI shows no drought on day 1 to 72 (winter, early spring), day 169 to 248 (summer, early autumn) and day 289 to 365 (366 days a year) (late autumn, winter), which may be caused by the lowest temperature and the lesser potential evapotranspiration on day 1 to 72 (winter, early spring) and day 289 to 365 (366 days a year) (late autumn, winter). But in the period of day 169 to 248, the temperature is higher, the evapotranspiration is stronger and the precipitation is abundant, so the surface water can be replenished in time. CWSI had different degrees of drought on 73 to 168 and day 249 to 288 and on day 113 to 144, it appears to be severe drought. The main reason may be that it is in the spring snow melt period, the temperature rises and the appropriate temperature and moisture are used for forage returning to green, so there is no enough temperature and moisture for evapotranspiration [31][32][33]. While during day 249 to 288, it is in the autumn vegetation withering period, the temperature begins to decrease and the evapotranspiration of vegetation weakens [34]. Therefore, in this period, the overall situation of Inner Mongolia is drought. forage returning to green, so there is no enough temperature and moisture for evapotranspiration [31][32][33]. While during day 249 to 288, it is in the autumn vegetation withering period, the temperature begins to decrease and the evapotranspiration of vegetation weakens [34]. Therefore, in this period, the overall situation of Inner Mongolia is drought. In Inner Mongolia, CWSI fluctuates obviously in the year (Figure 8). In the eastern and central regions (except Baotou city), the changes to the Bayannur city and the whole region are basically the same. However, the CWSI in Baotou city and western region (except Bayannur city) shows trapezoidal curves with no obvious bottom value, which increase rapidly in day 1 to 112, stabilize between 0.8648 and 0.9801 in day 113-280 and decrease rapidly in day 281-365 (366 days a year). This may be due to the lower precipitation, higher temperature, weaker actual evapotranspiration and stronger potential evapotranspiration in these areas during the rainy season, in turn, the value of CWSI does not fluctuate significantly at 113-280 days.

Characteristics of CWSI Changes in Different Land Use Types
The overall difference in CWSI of different land use types is small (Figure 9). The annual average of CWSI is from small to large: woodland (0.5954) < cropland (0.7733) < built-up land (0.8126) < grassland (0.8147) < unused land (0.8392); the annual mean values of ET for different land use types were: woodland (352.03 mm) > cropland (266.84 mm), > In Inner Mongolia, CWSI fluctuates obviously in the year (Figure 8). In the eastern and central regions (except Baotou city), the changes to the Bayannur city and the whole region are basically the same. However, the CWSI in Baotou city and western region (except Bayannur city) shows trapezoidal curves with no obvious bottom value, which increase rapidly in day 1 to 112, stabilize between 0.8648 and 0.9801 in day 113-280 and decrease rapidly in day 281-365 (366 days a year). This may be due to the lower precipitation, higher temperature, weaker actual evapotranspiration and stronger potential evapotranspiration in these areas during the rainy season, in turn, the value of CWSI does not fluctuate significantly at 113-280 days.
Sustainability 2021, 13, x FOR PEER REVIEW forage returning to green, so there is no enough temperature and moisture for ev spiration [31][32][33]. While during day 249 to 288, it is in the autumn vegetation w period, the temperature begins to decrease and the evapotranspiration of ve weakens [34]. Therefore, in this period, the overall situation of Inner Mongolia is In Inner Mongolia, CWSI fluctuates obviously in the year (Figure 8). In the and central regions (except Baotou city), the changes to the Bayannur city and th region are basically the same. However, the CWSI in Baotou city and western re cept Bayannur city) shows trapezoidal curves with no obvious bottom value, w crease rapidly in day 1 to 112, stabilize between 0.8648 and 0.9801 in day 113decrease rapidly in day 281-365 (366 days a year). This may be due to the lower tation, higher temperature, weaker actual evapotranspiration and stronger potent otranspiration in these areas during the rainy season, in turn, the value of CWSI fluctuate significantly at 113-280 days.

Characteristics of CWSI Changes in Different Land Use Types
The overall difference in CWSI of different land use types is small (Figure annual average of CWSI is from small to large: woodland (0.5954) < cropland ( built-up land (0.8126) < grassland (0.8147) < unused land (0.8392); the annual mea of ET for different land use types were: woodland (352.03 mm) > cropland (266.8

Characteristics of CWSI Changes in Different Land Use Types
The overall difference in CWSI of different land use types is small (Figure 9). The annual average of CWSI is from small to large: woodland (0. Surface drought is mainly related to land use type, geographical location and climate. For Inner Mongolia, most of the woodland are located in the Da Hinggan Mountains, with abundant rainfall, lower temperature and better water conservation ability. The ET value is higher and PET value is lower here, thus the CWSI value is lower, which shows that this area has the strongest drought resistance capability and the lowest risk of drought. Cropland plants artificial vegetation and has irrigation water supply, plus the ET value is higher, so the CWSI value is low. The vegetation cover of built-up land is less and because of the heat island effect and canyon effect, the temperature here is higher than that of surrounding area and strong wind is easier to happen in partial urban town and the largest PET value indicates that the risk of drought in built-up land is high; because the grassland are small, short and dense, evapotranspiration makes it consume more water, while unused land is mostly desert and sandy land, vegetation coverage is low, precipitation is rare and evapotranspiration results in faster water loss, thus the risk of drought is higher in grassland and unused land.

Correlation between CWSI and Climate Factors
The spatial distribution of meteorological elements in Inner Mongolia is as follows: The precipitation gradually decreases from northeast to southwest and the high value region appears in the Da Hinggan Mountains. The precipitation gradually decreased from northeast to southwest and the high value region appeared in the Da Hinggan Mountains (Figure 10). Correlation coefficient analysis is used to analyze the spatial distribution of the correlation between the annual average CWSI value and the annual average precipitation and the annual mean temperature on the pixel scale from 2001 to 2017, as shown in Figure 11 Through statistical analysis, the average correlation coefficient r of CWSI with precipitation and temperature is −0.53 and 0.18, respectively. Comparatively, the correlation between CWSI and precipitation is stronger. Surface drought is mainly related to land use type, geographical location and climate. For Inner Mongolia, most of the woodland are located in the Da Hinggan Mountains, with abundant rainfall, lower temperature and better water conservation ability. The ET value is higher and PET value is lower here, thus the CWSI value is lower, which shows that this area has the strongest drought resistance capability and the lowest risk of drought. Cropland plants artificial vegetation and has irrigation water supply, plus the ET value is higher, so the CWSI value is low. The vegetation cover of built-up land is less and because of the heat island effect and canyon effect, the temperature here is higher than that of surrounding area and strong wind is easier to happen in partial urban town and the largest PET value indicates that the risk of drought in built-up land is high; because the grassland are small, short and dense, evapotranspiration makes it consume more water, while unused land is mostly desert and sandy land, vegetation coverage is low, precipitation is rare and evapotranspiration results in faster water loss, thus the risk of drought is higher in grassland and unused land.

Correlation between CWSI and Climate Factors
The spatial distribution of meteorological elements in Inner Mongolia is as follows: The precipitation gradually decreases from northeast to southwest and the high value region appears in the Da Hinggan Mountains. The precipitation gradually decreased from northeast to southwest and the high value region appeared in the Da Hinggan Mountains (Figure 10). Correlation coefficient analysis is used to analyze the spatial distribution of the correlation between the annual average CWSI value and the annual average precipitation and the annual mean temperature on the pixel scale from 2001 to 2017, as shown in Figure 11 Through statistical analysis, the average correlation coefficient r of CWSI with precipitation and temperature is −0.53 and 0.18, respectively. Comparatively, the correlation between CWSI and precipitation is stronger. (Figure 10). Correlation coefficient analysis is used to analyze the spatial distributi the correlation between the annual average CWSI value and the annual average pr tation and the annual mean temperature on the pixel scale from 2001 to 2017, as sho Figure 11 Through statistical analysis, the average correlation coefficient r of CWSI precipitation and temperature is −0.53 and 0.18, respectively. Comparatively, the co tion between CWSI and precipitation is stronger.  The Correlation between CWSI and precipitation, temperature from 2001 to 201 calculated by Formula (3). The results are shown in Figure 11. Statistical results show there is a negative correlation between CWSI and precipitation as a whole, which in ing that CWSI decreases with the increase of precipitation; the statistical results show the area with negative correlation is 96.58% and the area of 65.95% is significantly neg correlation at the significant level of 0.05, which is concentrated in the steppe area only 0.02% of the area is positively correlated, that is mainly because CWSI is affect vegetation coverage, temperature, sunshine hours, solar radiation, slope, elevation, speed and other factors. There is a positive correlation between CWSI and temper which indicating that CWSI increases with the increase of temperature. Among them area with positive correlation accounted for 82.04% of the total area, while at the si cant level of 0.05 it only accounted for 5.65% of the total area. It is concentrated in Barhu Banner, Chifeng City, Ulanchab City, Hohhot City and Baotou City. It show these regions are greatly affected by temperature. According to the analysis above, is closely related to rainfall and temperature, the increase of precipitation makes the tive soil moisture increase and the water supply is sufficient. In addition, the lower perature reduces evapotranspiration and decreases the rate of water loss, which is ficial to the growth of vegetation. Therefore, in the regions with more precipitatio lower temperature, the vegetation cover degree is higher, the CWSI value is smalle the drought degree is less.

Discussion
The existing research on the temporal and spatial variation of drought in Inner golia is mainly focused on the long-term meteorological observation data, such a analysis of the characteristics and causes of temporal and spatial changes of SPEI an Based on the Penmen-Monteith formula and meteorological data, Zhang Xuting and ers analyzed the temporal and spatial variation characteristics of drought in Inner golia from 1960 to 2015 by calculating SPEI in different time scales combined with va statistical methods, and the relationship between drought characteristics and Pacif terdecadal Oscillation Index is discussed [34]. An Qiang and others analyzed the tem and spatial variation characteristics of drought in Inner Mongolia from 1958 to 20 calculating SPEI in different time scales and combined with the influence of contin drought, the comprehensive evaluation index of drought intensity was put forward Liu Jiyao and others used monthly observational meteorology data, combined with trend analysis and SPI analysis to study the spatial and temporal variations of drou Inner Mongolia from 1961 to 2013 [36]. These studies illustrate in detail the variati drought in Inner Mongolia on a certain scale and have important practical significan drought monitoring and decision making but they may be confined to the "point" or the lack of large-scale and long-term continuous observation data, therefore they a general understanding of the drought situation in the whole research area. In this p based on the analysis of the temporal and spatial variation characteristics of ET, The Correlation between CWSI and precipitation, temperature from 2001 to 2017 are calculated by Formula (3). The results are shown in Figure 11. Statistical results shows that there is a negative correlation between CWSI and precipitation as a whole, which indicating that CWSI decreases with the increase of precipitation; the statistical results show that the area with negative correlation is 96.58% and the area of 65.95% is significantly negative correlation at the significant level of 0.05, which is concentrated in the steppe area; And only 0.02% of the area is positively correlated, that is mainly because CWSI is affected by vegetation coverage, temperature, sunshine hours, solar radiation, slope, elevation, wind speed and other factors. There is a positive correlation between CWSI and temperature, which indicating that CWSI increases with the increase of temperature. Among them, the area with positive correlation accounted for 82.04% of the total area, while at the significant level of 0.05 it only accounted for 5.65% of the total area. It is concentrated in Chen Barhu Banner, Chifeng City, Ulanchab City, Hohhot City and Baotou City. It shows that these regions are greatly affected by temperature. According to the analysis above, CWSI is closely related to rainfall and temperature, the increase of precipitation makes the relative soil moisture increase and the water supply is sufficient. In addition, the lower temperature reduces evapotranspiration and decreases the rate of water loss, which is beneficial to the growth of vegetation. Therefore, in the regions with more precipitation and lower temperature, the vegetation cover degree is higher, the CWSI value is smaller and the drought degree is less.

Discussion
The existing research on the temporal and spatial variation of drought in Inner Mongolia is mainly focused on the long-term meteorological observation data, such as the analysis of the characteristics and causes of temporal and spatial changes of SPEI and SPI.
Based on the Penmen-Monteith formula and meteorological data, Zhang Xuting and others analyzed the temporal and spatial variation characteristics of drought in Inner Mongolia from 1960 to 2015 by calculating SPEI in different time scales combined with various statistical methods, and the relationship between drought characteristics and Pacific Interdecadal Oscillation Index is discussed [34]. An Qiang and others analyzed the temporal and spatial variation characteristics of drought in Inner Mongolia from 1958 to 2019 by calculating SPEI in different time scales and combined with the influence of continuous drought, the comprehensive evaluation index of drought intensity was put forward [35]. Liu Jiyao and others used monthly observational meteorology data, combined with linear trend analysis and SPI analysis to study the spatial and temporal variations of drought in Inner Mongolia from 1961 to 2013 [36]. These studies illustrate in detail the variation of drought in Inner Mongolia on a certain scale and have important practical significance in drought monitoring and decision making but they may be confined to the "point" scale or the lack of large-scale and long-term continuous observation data, therefore they lack a general understanding of the drought situation in the whole research area. In this paper, based on the analysis of the temporal and spatial variation characteristics of ET, PET, CWSI, we use high-resolution remote sensing products and meteorological observation data of continuous time period to further discuss the variation trend of CWSI over the years, which reflects the spatial and temporal distribution of drought in the region.
The results can reflect more information of surface drought characteristics. Moreover, it can provide a new method for large scale inhomogeneous drought monitor.
Compared to meteorological data, MOD16 data has a shorter time scale, which is only after 2000. What's more, due to the influence of cloud, shadow and vegetation cover on MOD16 data, some pixels have no value phenomenon, especially in winter, which is difficulty to study the drought in winter.
The relationship between ET, PET and CWSI in different land use types is discussed in this paper: forest land ET is the largest, PET is the smallest and CWSI is the smallest, which shows that the ability of forest land to resist drought is the strongest and it is of great significance to regional water conservation. The author has reclassified the land use data of 1 km resolution provided by the Chinese Academy of Sciences and the 30 m resolution land use data provided by Tsinghua University according to the unified standard and the CWSI extracted from different land use types were compared separately ( Figure 12). It can be seen from the figure that the results obtained from the two data sources are basically the same and the average absolute error is only −0.0034, which indicates that the reliability of CWSI extracted from 1 km land use data provided by the Chinese Academy of Sciences is higher; this paper may underestimate the drought resistance of woodland and construction land, the possibility of drought risk on bare land and overestimate the drought resistance of cropland and grassland. This paper discusses the correlation between CWSI and precipitation and temperature and CWSI in Inner Mongolia has a significant correlation with precipitation, indicating that precipitation has a greater impact on CWSI, which is consistent with Park and others view that precipitation is the main cause of drought [37]. However, the process of CWSI change is very complex, which is influenced by many meteorological factors, such as altitude, slope, slope direction and vegetation coverage and the influence mechanism of these factors on different regions is different [38]. Therefore, according to the actual situation of the surface of the study area, in order to strengthen the ability of drought prevention and drought resistance and improve the regional climate environment, the focus of further research is the contribution of different factors to CWSI. lution land use data provided by Tsinghua University according to the unified standard and the CWSI extracted from different land use types were compared separately ( Figure 12). It can be seen from the figure that the results obtained from the two data sources are basically the same and the average absolute error is only −0.0034, which indicates that the reliability of CWSI extracted from 1 km land use data provided by the Chinese Academy of Sciences is higher; this paper may underestimate the drought resistance of woodland and construction land, the possibility of drought risk on bare land and overestimate the drought resistance of cropland and grassland. This paper discusses the correlation between CWSI and precipitation and temperature and CWSI in Inner Mongolia has a significant correlation with precipitation, indicating that precipitation has a greater impact on CWSI, which is consistent with Park and others view that precipitation is the main cause of drought [37]. However, the process of CWSI change is very complex, which is influenced by many meteorological factors, such as altitude, slope, slope direction and vegetation coverage and the influence mechanism of these factors on different regions is different [38]. Therefore, according to the actual situation of the surface of the study area, in order to strengthen the ability of drought prevention and drought resistance and improve the regional climate environment, the focus of further research is the contribution of different factors to CWSI.

Conclusions
In this paper, MOD16 data are used to calculate crop water shortage index (CWSI) and the spatial and temporal variation characteristics of meteorological factors are analyzed, possible reasons of droughts were discussed and explored. The correlation between CWSI and the measured relative soil moisture data (p < 0. 05, r = −0.86) was negative. It is

Conclusions
In this paper, MOD16 data are used to calculate crop water shortage index (CWSI) and the spatial and temporal variation characteristics of meteorological factors are analyzed, possible reasons of droughts were discussed and explored. The correlation between CWSI and the measured relative soil moisture data (p < 0. 05, r = −0.86) was negative. It is shown that the CWSI calculated by MOD16 products can be used for drought monitoring in Inner Mongolia. The conclusions are as follows: (1) Droughts in 2001-2017 showed a significant decreasing trend (p < 0.05) at the interannual scale in terms of CWSI monitoring and the study area showed mild drought. The eight days variation characteristics of CWSI indicates that severe droughts occurred mainly during 113 to 144 days, this time are also the key period of natural forage regreening period. Therefore, frequent droughts will affect grassland productivity. It is suggested that the local government should pay more attention to drought resistance during the regreening period. (2) The annual change of droughts indicates that droughts occurred mainly in Ulanchab City, Ordos City, Baotou City, Bayan Nur, Alxa League and Wuhai City. While the drought intensity of Hulunbeier, Hinggan League, Chifeng City and Tongliao City were relatively weak. Spatially, more crops can be observed mainly in northern Ordos, southern Baotou and southern Ulanchab. These areas are also the key regions of severe drought and moderate drought. Frequent droughts will definitely affect crop yields. It is suggested that the local government should strengthen the use of irrigation water in crop growth period to prevent the economic losses caused by drought. (3) The annual average of droughts in different land use/cover types indicates that woodland (0.5954) < cropland (0.7733) < built-up land (0.8126) < grassland (0.8147) < unused land (0.8392). Hence, the ability of woodland to resist drought is the strongest and the possibility of drought in unused land is the highest. Therefore, the possibility of drought in Inner Mongolia can be reduced through afforestation, closing hillsides for afforestation. (4) The average correlation coefficients between CWSI and precipitation, temperature are −0.53 and 0.18. This indicates that CWSI is more sensitive to precipitation. Spatially, droughts of grassland and cropland will definitely be more sensitive to precipitation, while droughts of woodland is little or no affected by precipitation. Droughts of Chen Barhu Banner, Chifeng City, Ulanchab City, Hohhot City and Baotou City will definitely be more sensitive to temperature, while droughts of other areas is little or no affected by temperature. Therefore, precipitation has a certain indicative effect on drought monitoring of grassland and crops in the study area, while temperature also has a certain indicative effect on drought monitoring of Chen Barhu banner, Chifeng City, Ulanchab City, Hohhot City and Baotou city.
The results of this study may provide some practical reference for the aim to prevent drought and maintain grassland productivity in a warming climate.