MODIS and in Situ Data Analysis for Mongolia Drought

Drought is a period of abnormally dry weather with a serious shortage of water supply. Drought indices can be an advantageous indicator to assess drought for taking further response actions. However, drought indices based on ground meteorological measurements could not completely reveal the land use effects over a regional scale. On the other hand, the satellite-derived products provide consistent, spatial and temporal comparisons of global signatures for the regional-scale drought events. This research is to investigate the drought signatures over Mongolia by using satellite remote sensing imagery. The evapotranspiration (ET), potential evapotranspiration (PET) and two-band Enhanced Vegetation Index (EVI2) were extracted from MODIS data. Based on the standardized ratio of ET to PET (ET/PET) and EVI2, the Modified Drought Severity Index (MDSI) anomaly during the growing season from May–August for the years 2000–2013 was acquired. Fourteen-year summer monthly data for air temperature, precipitation and soil moisture content of in situ measurements from sixteen meteorological stations for four various land use areas were analyzed. We also calculated the percentage deviation of climatological variables at the sixteen stations to compare to the MDSI anomaly. Both comparisons of satellite-derived and observed anomalies and variations were analyzed by using the existing common statistical methods. The results demonstrated that the air temperature anomaly (T anomaly) and the precipitation anomaly (P anomaly) were negatively (correlation coefficient r =  ́0.66) and positively (r = 0.81) correlated with the MDSI anomaly, respectively. The MDSI anomaly distributions revealed that the wettest area occupied 57% of the study area in 2003, while the driest (drought) area occurred over 54% of the total area in 2007. The results also showed very similar variations between the MDSI and T anomalies. The highest (wettest) MDSI anomaly indicated the lowest T anomaly, such as in the year 2003, while the lowest (driest) MDSI anomaly had the highest T anomaly in 2007. By comparing the MDSI anomaly and soil moisture content at a 10-cm depth during the study period, it is found that their correlation coefficient is 0.74.


Introduction
Mongolia has a severe climate with high seasonal variations and a low annual average in rainfall [1].Based on climatic characteristics, the meteorological drought impacts on agricultural products are precipitation shortages, differences between actual and potential evapotranspiration (PET) and other physiologic factors of the soil and biomass.Hence, poor biomass production would cause huge economic losses for the livestock industry due to the drought impacts in summer [2].For example, tens of millions of livestock were lost after the drought and harsh winter in the last decade, and thus, the lives of hundreds of thousands of people were forced to change.Monitoring and even forecasting of drought becomes crucial to Mongolia.Drought indices that consist of a large amount of data, such as precipitation, snowpack, streamflow and other water supply indicators, can be applied to monitor drought severities in a comprehensive framework [3].
It is known that the major limitation of the ground in situ data is the observational data that cannot be effectively used to monitor and evaluate the drought events over a large area.By contrast, remote sensing techniques provide a better solution for monitoring drought conditions over a large scale at early stages because of its wide coverage, high spatial resolution, easy access and relatively lower expense [4].In the field of remote sensing, reflectance and thermal datasets of satellite images are used for developing drought indices.One of the most reliable sources of remote sensing data is the Moderate Resolution Imaging Spectroradiometer (MODIS) sensor on board the Terra and Aqua satellites that provides information on vegetation and surface energy conditions [5].Furthermore, the National Oceanic Atmospheric Administration (NOAA) Advanced Very High Resolution Radiometer (AVHRR)-derived Normalized Difference Vegetation Index (NDVI) and other related indices (e.g., NDVI anomaly, integrated and standardized NDVI, Vegetation Condition Index, etc.) have been successfully used to identify and monitor areas at the regional and local scales [6][7][8][9][10][11][12][13].
Satellite vegetation indices (VIs), especially NDVI and Enhanced Vegetation Index (EVI), have been commonly used for monitoring global vegetation photosynthetic activities [14][15][16] and estimating gross primary production [17].For instance, VIs from MODIS sensors have great improvements in spatial, spectral and radiometric measurements of surface vegetation conditions [18], and NDVI could be used to detect the changes of surface vegetation conditions [19][20][21][22], to monitor the maize green leaf area index [23], to estimate evapotranspiration [24,25], to conduct flood mapping and disaster loss assessment [26,27] and to perform environmental vulnerability assessment [28].In addition, many studies showed that NDVI is a useful index for studying vegetation and ecosystems in semi-arid environments where vegetation cover is less than 30% [29,30].To improve upon NDVI, the EVI is optimized to enhance the vegetation signal and to reduce the soil background influences [31,32].Furthermore, Kim et al., 2010 [33], investigated the multi-sensor spectral bandpass dependencies of EVI, EVI2 and NDVI using Hyperion satellite images acquired over a range of vegetation conditions.
Yao et al. [34] developed the Evaporative Drought Index (EDI), which is based on evapotranspiration (ET) and PET of MODIS data, to estimate the deviation of the ratio between ET and PET (ET/PET) from unity.The EDI utilizes reanalysis of meteorological data and remotely-sensed data, and the MODIS data MOD16 ET/PET products were used to quantify water fluxes between terrestrial ecosystems and the atmosphere [35][36][37][38].
Mu et al. [39] derived the Drought Severity Index (DSI) to integrate operational MOD16 ET/PET and MOD13 NDVI products to calculate a remotely-sensed global drought index and then evaluated drought events at regional scale.The advantage of DSI is in that it could combine a wide range of historic climate data to quantify how ecosystems respond to drought events.
Using both satellite-derived and observed data, we examined drought occurrence over Mongolia.In our method, we used the two-band Enhanced Vegetation Index (EVI2) to replace the NDVI that was traditionally used.A new drought index, the Modified Drought Severity Index (MDSI) anomaly, which was derived from MOD16 ET/PET and MOD13 EVI2 based on the MODIS dataset and compared to ground-based observation data, is thus proposed.EVI2 is adopted to represent the vegetation status because EVI2 exhibits a nearly identical behavior for well-preprocessed MODIS data, like EVI [16], and EVI has a higher sensitivity to moderate-to-high vegetation biomass [31].The observation data were obtained from the sixteen meteorological stations over Mongolia.The time series of both satellite-derived and observed data were acquired for the period from May-August of 2000-2013.In order to compare the evaluation of the drought occurrence of the study area, the land use area was divided into four various land use areas, namely the forest steppe, steppe, high mountains and desert steppe.

Study Area
The land use of Mongolia is classified into five subtypes: the high mountainous region (4.5%), forest steppe (23%), steppe (28%), desert steppe (28%) and desert (16%) [40].The northern part of Mongolia is covered by forested mountain ranges with a dry sub-humid climate, whereas the southern part encompasses the Gobi desert with lower elevation and a drier climate [41], as illustrated in Figure 1.
Remote Sens. 2016, 8, 509 3 of 16 series of both satellite-derived and observed data were acquired for the period from May-August of 2000-2013.In order to compare the evaluation of the drought occurrence of the study area, the land use area was divided into four various land use areas, namely the forest steppe, steppe, high mountains and desert steppe.

Study Area
The land use of Mongolia is classified into five subtypes: the high mountainous region (4.5%), forest steppe (23%), steppe (28%), desert steppe (28%) and desert (16%) [40].The northern part of Mongolia is covered by forested mountain ranges with a dry sub-humid climate, whereas the southern part encompasses the Gobi desert with lower elevation and a drier climate [41], as illustrated in Figure 1.In Mongolia, precipitation is significantly different in different regions, and temperature increases from north to south.The southern parts of Mongolia have higher average temperatures, which range from 20-25 °C.By contrast, the northern parts have lower temperatures of 15-20 °C.The annual total precipitation is around 400 mm in the mountainous areas, 250-300 mm in the forest steppe regions, 150-250 mm in the steppe regions and less than 100 mm in the desert area.Most of rainfall, 85%-90%, would happen in the summer season.It is the main source of soil moisture, but the amount of rainfall is insufficient for the growing season period [42,43].
Clearly, Mongolia's climatic feature is a function of latitude characterized by the natural vegetation cover in the geographical zones.These zones are divided into natural vegetation belts, which are located at different altitudes (from mountains to plain steppe) and latitudes (from the north to south).The natural vegetation has a very short growing season (April-September), and pasture would start to grow from the middle of May to the middle of August [44].

MODIS Products
We selected a fourteen-year period to observe the dynamic range of the drought index during the terrestrial growing season with different climatological conditions.Accordingly, a monthly composite of the 1-km spatial resolution MOD16A3 [45] and MOD13A3 [46] products from MODIS for free from the University of Montana's Numerical Terradynamic Simulation group [45] and the National Aeronautics and Space Administration (NASA) Earth Observing System (EOS) [47] were downloaded.Both MOD16A3 and MOD13A3 products were acquired from May-August of 2000-2013 with 224 images in total.In Mongolia, precipitation is significantly different in different regions, and temperature increases from north to south.The southern parts of Mongolia have higher average temperatures, which range from 20-25 ˝C.By contrast, the northern parts have lower temperatures of 15-20 ˝C.The annual total precipitation is around 400 mm in the mountainous areas, 250-300 mm in the forest steppe regions, 150-250 mm in the steppe regions and less than 100 mm in the desert area.Most of rainfall, 85%-90%, would happen in the summer season.It is the main source of soil moisture, but the amount of rainfall is insufficient for the growing season period [42,43].
Clearly, Mongolia's climatic feature is a function of latitude characterized by the natural vegetation cover in the geographical zones.These zones are divided into natural vegetation belts, which are located at different altitudes (from mountains to plain steppe) and latitudes (from the north to south).The natural vegetation has a very short growing season (April-September), and pasture would start to grow from the middle of May to the middle of August [44].

MODIS Products
We selected a fourteen-year period to observe the dynamic range of the drought index during the terrestrial growing season with different climatological conditions.Accordingly, a monthly composite of the 1-km spatial resolution MOD16A3 [45] and MOD13A3 [46] products from MODIS for free from the University of Montana's Numerical Terradynamic Simulation group [45] and the National Aeronautics and Space Administration (NASA) Earth Observing System (EOS) [47] were downloaded.Both MOD16A3 and MOD13A3 products were acquired from May-August of 2000-2013 with 224 images in total.
MOD16A3 produces a monthly composite of MODIS Terrestrial ET and PET products.The ratio of ET to PET is commonly used as an indicator of terrestrial water availability, and associated wetness or drought was modified by Mu et al. [39] and can be described as: The MOD13A3 product is generated using the 16-day composite period MODIS Terra vegetation index (VI) output.It uses the time-weighted average of the reflectance fields in the 16-day composited VI product that falls within a particular month and the ones that overlap in the beginning and end of each calendar month as determined by Huete et al. [48].EVI can really reflect the vegetation conditions, which can reduce the effect on vegetation information by the soil background through adopting the principle of the soil adjusted vegetation index as presented by Wang et al. [49].The modification of EVI, EVI2 [16], is calculated from the red and near-infrared (NIR) bands as follows: where ρ red and ρ NIR are the surface reflectance of the red and NIR bands, respectively; 2.5 is a gain factor, and 1 is a canopy background adjustment as described by Huete et al. [48,50]; 2.4 is the coefficient of red reflectance in the denominator as determined by Jiang et al. [16].
Both ET/PET and EVI2 of the MODIS images were merged and projected using the HDF-EOS to GeoTIFF Conversion Tool (HEG).This tool allows the projection from the MODIS Sinusoidal (SIN) grid to the Universal Transverse Mercator (UTM) Zone 48 N projection in our analysis.

In Situ Measurement Data
In situ measurement data based on air temperature, precipitation and soil moisture observations were obtained at the sixteen meteorological observations stations from 2000-2013 for the growing season.The vegetation growing season data were acquired between May and August.Air temperature ( ˝C) and precipitation (mm) data were used to acquire their corresponding monthly mean, which is based on daily observation data.Soil moisture content (by gravimetric method) was obtained at a 10-cm depth at a monthly interval (27th and 28th day of each month) for four consecutive months (May-August).These in situ observation data were provided by the Information and Research Institute of Meteorology, Hydrology and Environment (IRIMHE) of Mongolia (http://www.icc.mn/)[51].There are no missing air temperature and precipitation data in the time period of investigation.Note that, in the years 2001 and 2002, 3-5-times of the soil moisture data were missing.
The sixteen stations are widely distributed across the different climatic regions of Mongolia, consisting of various land use areas, forest steppe, steppe, high mountains and desert steppe, as illustrated in Figure 1.The geographical locations, elevations and various land use area classifications of the sixteen stations are represented in Table 1.
Air temperature, precipitation and soil moisture content data were averaged over the monthly intervals from May-August to match the MODIS monthly composite period.For interpolation, continuous grids of the MDSI anomaly value with the same geographical location and resolution (0.1 ˝ˆ0.1 ˝) as the in situ measurement data were created.The sixteen stations are listed in order from north to south and west to east.

Methodology
In this study, we used ET/PET and EVI2 derived from MODIS and climatological variables measured from ground stations to investigate the spatiotemporal drought occurrences of satellite-derived and observed anomalies during the growing season of 2000-2013.The overall description of the study consists of three main parts: (1) the spatial variations of the MDSI anomaly; (2) the temporal variations of the climatological variables' anomaly; and (3) the relationship between the MDSI anomaly with meteorological variables and soil moisture content.

Anomaly of MDSI Description
A new monthly composite MDSI anomaly for ET/PET and EVI2 variables was calculated on a pixel by pixel basis.We investigated the spatial variations of the monthly MDSI anomaly (Equation ( 6)) for the entire dataset (Z Ratio , Z EVI2 ) used to evaluate changes in the variation during study years.The MDSI anomaly for a variable is calculated as follows.The standardized ratio (Z Ratio ) is calculated from Equation (1): The standardized EVI2 (Z EVI2 ) is calculated from Equation (2): The Z Ratio and Z EVI2 terms are then summed as: A new remotely-sensed MDSI anomaly is finally calculated as the standardized Z terms: where Z represents the variable value for the sum of the Z Ratio and Z EVI2 terms.Z and σ Z are the long-term mean (2000-2013) and standard deviation for the Z Ratio and Z EVI2 terms, respectively.Thus, the MDSI anomaly is a dimensionless index theoretically ranging from unlimitedly positive values to unlimitedly negative values for wet to dry period deviations, respectively, from prevailing conditions.

Anomaly of Climatological Variables Description
To understand the temporal variations of climatological variables for the region of interest, we generated a time series of meteorological data of air temperature and precipitation averaged over the sixteen stations' scale.As described by Dhar et al. [52], monthly anomalies for the entire air temperature and precipitation datasets were used to evaluate changes in the variation during the study period.The percentage deviations of variables from their long-term mean (1981-2010) and the air temperature anomaly (T anomaly) and the precipitation anomaly (P anomaly) are given by the following equations: P anomaly " ˆPm ´Pm where T anomaly and P anomaly are the percentage deviations from the long-term mean of air temperature and precipitation, respectively.T m and P m are the monthly air temperature and precipitation values, respectively.T m , and P m are the long-term means of the monthly air temperature and precipitation, respectively.
Monthly anomalies based on the percentage deviation (T anomaly, P anomaly) were then used to categorize the drought condition.The category of the drought condition would have a T anomaly value of 0.0 or above, while the P anomaly would be ´0.0 or less.We then merged four consecutive (May-August) monthly data for the growing season of the study years to obtain the mean values of MDSI, T and P anomalies and soil moisture content.

Relationship between MDSI Anomaly and Climatological Variables
Proposed approaches were followed to explore the relation between the MDSI anomaly and climatological variables.For the category of drought condition, remotely-sensed and ground observation data were utilized.We obtained the monthly MDSI anomaly from MODIS data and compared to the monthly air temperature and precipitation, as well as their anomalies of the percentage deviation and soil moisture content from ground observation data.In addition, to quantify the comparison results, the Pearson's correlation coefficients (r) for the MDSI anomaly vs. the T anomaly, the MDSI anomaly vs. the P anomaly and the MDSI anomaly vs. the soil moisture content were computed for the sixteen stations.The formula of the coefficient of Pearson (r) is given in Equation ( 9) as: r " where X i and Y i are the individual derivations and measurements of variables X and Y, respectively.X and Y are the means of X and Y, respectively.The linear Pearson's correlation (r) was investigated at the monthly timescale (May-August) for satellite-derived and observed anomalies, as well as soil moisture content values for the period of 2000-2013, and statistical significance level (p) was chosen at the 90%-99% in this research.

Climatological Variables among Stations
Table 2 presents the temporal variations of climatological variables, including air temperature, precipitation and soil moisture content during the growing season for the study area.The results showed that the average air temperature of sixteen stations in Mongolia was approximately 17.0 ˝C; the mean precipitation was about 149 mm; and the average soil moisture content was 6.7% during the study period (2000)(2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013).In four various land uses, temperatures were increased from the high mountains to desert steppe, but precipitation was decreased from the forest steppe to desert steppe.The results demonstrated that the forest steppe had lower temperature, but higher precipitation and soil moisture content, while the desert steppe had the opposite trends.Furthermore, the high mountains region had the lowest temperature and soil moisture content values due to its unique terrain features and dependence on rainfall.Table 3 presents the correlation matrix among different climate parameters, and the results showed that the sixteen stations had strong and significant correlations.The inter-station correlation matrix (r) between climatological variables showed that the correlation coefficient of Pearson (r) (Equation ( 9)) ranged from negative 0.77 to positive 0.78 with significance p < 0.001.Basically, the air temperature tends to be in hot and dry conditions, while the precipitation tends to occur in cool and wet periods.This indicates that the air temperature and precipitation were inversely correlated with each other.Therefore, surface soil moisture was related to the amount of precipitation.In addition, the temporal variations of the T and P anomalies' results will be shown in Section 4.3.Note that the negative correlation between T and soil moisture may indicate stronger plant-available water stress as temperature rises [53].Table 3 presents the correlation matrix among different climate parameters, and the results showed that the sixteen stations had strong and significant correlations.The inter-station correlation matrix (r) between climatological variables showed that the correlation coefficient of Pearson (r) (Equation ( 9)) ranged from negative 0.77 to positive 0.78 with significance p < 0.001.Basically, the air temperature tends to be in hot and dry conditions, while the precipitation tends to occur in cool and wet periods.This indicates that the air temperature and precipitation were inversely correlated with each other.Therefore, surface soil moisture was related to the amount of precipitation.In addition, the temporal variations of the T and P anomalies' results will be shown in Section 4.3.Note that the negative correlation between T and soil moisture may indicate stronger plant-available water stress as temperature rises [53].

Anomaly of MDSI Variations
The averaged MDSI anomaly was calculated as the monthly standardized ET/PET and EVI2 throughout the growing season, four consecutive months from May-August, for the years 2000-2013.The growing season of the MDSI anomaly was estimated for both wet and dry periods in the study area.The MDSI anomaly indicates the categories of wet and dry conditions with negative values referring to drier conditions as compared to normal conditions and positive values implying wetter conditions.To classify drought intensity, the anomalies for ET/PET and EVI2 were standardized and compared.In this study, we used the MDSI classification ranges corresponding to the Palmer Drought Severity Index (PDSI) [54] drought severity categories, where the W1-2 (W, wet) and D1-2 (D, dry) categories indicate progressively wetter and drier conditions, respectively.Thus, the MDSI anomaly value ranges between unlimitedly positive values (>2) and unlimitedly negative values (<´2).Then, the MDSI anomaly value would be around ´0.5 or less for a drought condition and 0.5 or greater for a wet condition.For normal conditions, it would be between 0.49 and ´0.49[39,54], as shown in Table 4.   3h) and the southwest region in 2009 (Figure 3j) were less than ´0.5.In addition, the summers in 2003 and 2011-2013 (Figure 3d,l-n) had more favorable hydrologic conditions than the other years (Figure 3a,c,e,i,k) with fairly positive MDSI anomaly values of 0.5 or greater.
The values of the MDSI anomaly were attributed to five categories (Table 4), namely severely wet (>2), moderately wet (0.5 to 1.99), normal (0.49 to ´0.49), moderate drought (´0.5 to ´1.99) and severe drought (<´2), in this study.As illustrated in Figure 4, the frequency distribution of the MDSI anomaly showed a wet period in the years 2003 and 2011-2013 with 85%-92% of the total area identified as normal and moderately wet, while 8%-15% of it was moderate drought conditions.In contrast, dry periods appeared in the years 2001, 2007 and 2009, and 45%-54% of the total area was identified as moderate drought, but 46%-55% of it was observed as normal and moderately wet conditions.In the other years, 18%-32% of the total area occurred as moderate drought, while more areas (68%-82%) were identified as normal and moderately wet.Furthermore, a maximum of 9% occupied high intensity wetness (severely wet >2) in the year 2003.In the other years, wetness regions do not exceed 0.4% of the total area.Furthermore, the high intensity dryness of MDSI anomaly (severe drought <´2) does not exceed 0.2% of the total area during the study period.Finally, the performance of the MDSI anomaly distributions revealed that the wettest area occupied 57% of the total area in 2003, while the driest area occurred 54% of the total area in 2007.The MDSI anomaly was computed as the standardized Z terms represented in Equation (6).

Anomaly of MDSI and Climatological Variable Variations
The climatological variables' anomaly was used to identify wet and dry conditions based on Equations ( 7) and ( 8) of the percentage deviations and compared to the MDSI anomaly, as shown in Figure 5. Figure 5a showed that the T anomaly was negatively correlated with the MDSI anomaly (r = ´0.66 with p < 0.01) during the growing season of 2000-2013.Besides, by comparing MDSI with T anomalies, the highest MDSI anomaly value (+0.8) indicated the lowest T anomaly (2%), as observed in the year 2003, while the lowest MDSI anomaly value (´0.8) had the highest T anomaly (20%), as seen in the year 2007.These results are also consistent with the temporal variations, as shown in Figure 4. Figure 5b illustrates that there are variations in the MDSI with P anomalies in different regions.Both the distributions and ranges of MDSI with P anomalies varied with year.A comparison between them indicates that they are positively correlated with a coefficient of 0.81 and significance p < 0.001.

Anomaly of MDSI and Climatological Variable Variations
The climatological variables' anomaly was used to identify wet and dry conditions based on Equations ( 7) and ( 8) of the percentage deviations and compared to the MDSI anomaly, as shown in Figure 5. Figure 5a showed that the T anomaly was negatively correlated with the MDSI anomaly (r = −0.66 with p < 0.01) during the growing season of 2000-2013.Besides, by comparing MDSI with T anomalies, the highest MDSI anomaly value (+0.8) indicated the lowest T anomaly (2%), as observed in the year 2003, while the lowest MDSI anomaly value (−0.8) had the highest T anomaly (20%), as seen in the year 2007.These results are also consistent with the temporal variations, as shown in Figure 4. Figure 5b illustrates that there are variations in the MDSI with P anomalies in different regions.Both the distributions and ranges of MDSI with P anomalies varied with year.A comparison between them indicates that they are positively correlated with a coefficient of 0.81 and significance p < 0.001.There are sixteen stations distributed over four various land use areas with different soil moisture conditions during the study period.As an example, the soil moisture content at 10 cm was selected to compare to the MDSI anomaly on a station basis, as shown in Figure 6.The temporal variations showed a better relationship between the MDSI anomaly and soil moisture content, except for 2001 and 2002.It can be found that the MDSI anomaly had a strong significant correlation (r = 0.74 with p < 0.01) with surface soil moisture.It is feasible to monitor the drought status by using the MDSI anomaly, which is calculated from the ratio of ET to PET and the EVI2 feature space.In There are sixteen stations distributed over four various land use areas with different soil moisture conditions during the study period.As an example, the soil moisture content at 10 cm was selected to compare to the MDSI anomaly on a station as shown in Figure 6.The temporal variations showed a better relationship between the MDSI anomaly and soil moisture content, except for 2001 and 2002.It can be found that the MDSI anomaly had a strong significant correlation (r = 0.74 with p < 0.01) with surface soil moisture.It is feasible to monitor the drought status by using the MDSI anomaly, which is calculated from the ratio of ET to PET and the EVI2 feature space.In addition, the temporal and spatial MDSI anomaly clearly shows that the soil moisture content of the drier seasons is linked to rainfall, as well as the characteristics and topography of various land use areas.

Anomalies among Different Land Use Areas
Table 5 presents the correlations of the MDSI anomalies with different climatological variables and soil moisture values over all of Mongolia for the study period.It indicates that the values of Pearson's correlation coefficient (r) were −0.66 between the MDSI and T anomalies and 0.81 between the MDSI and P anomalies with significance p < 0.001 over all of Mongolia.In four various land use areas, the correlation coefficients varied from −0.67 to −0.52 between the MDSI and T anomalies, and 0.67-0.83between the MDSI and P anomalies.However, T and P with MDSI anomalies had a low value of −0.45 (p > 0.05) and 0.34 (p > 0.10) in the desert steppe and forest steppe, respectively.Soil moisture content can reflect precipitation, evaporation, infiltration and runoff.In turn, it acts as a strong control on the partitioning of water between the atmosphere and the surface.As shown in Table 5, the MDSI anomaly and soil moisture content appears to have a high correlation of 0.74 as p < 0.01 over Mongolia.As for individual land use area, the correlations of MDSI anomaly with soil moisture content were rather high 0.63 for forest steppe, 0.64 for high mountains and 0.72 for desert steppe with significance p < 0.01, except with a low value of 0.39 (p > 0.10) for the steppe.
Figure 7 illustrates the temporal variations of the MDSI anomaly for four various land use areas and related to the climate conditions at the station scale.It represents that the MDSI anomaly exhibited high and low temporal variations in response to wet and dry conditions with climatic parameters.In Figure 7a

Anomalies among Different Land Use Areas
Table 5 presents the correlations of the MDSI anomalies with different climatological variables and soil moisture values over all of Mongolia for the study period.It indicates that the values of Pearson's correlation coefficient (r) were ´0.66 between the MDSI and T anomalies and 0.81 between the MDSI and P anomalies with significance p < 0.001 over all of Mongolia.In four various land use areas, the correlation coefficients varied from ´0.67 to ´0.52 between the MDSI and T anomalies, and 0.67-0.83between the MDSI and P anomalies.However, T and P with MDSI anomalies had a low value of ´0.45 (p > 0.05) and 0.34 (p > 0.10) in the desert steppe and forest steppe, respectively.Soil moisture content can reflect precipitation, evaporation, infiltration and runoff.In turn, it acts as a strong control on the partitioning of water between the atmosphere and the surface.As shown in Table 5, the MDSI anomaly and soil moisture content appears to have a high correlation of 0.74 as p < 0.01 over Mongolia.As for individual land use area, the correlations of MDSI anomaly with soil moisture content were rather high 0.63 for forest steppe, 0.64 for high mountains and 0.72 for desert steppe with significance p < 0.01, except with a low value of 0.39 (p > 0.10) for the steppe.

Conclusions
The MDSI anomaly was developed for Mongolia based on the standardized ET/PET and EVI2 index extracted from the MODIS dataset during the growing season from May-August of 2000-2013.Furthermore, the fourteen-year summer monthly data for temperature, precipitation and soil

Conclusions
The MDSI anomaly was developed for Mongolia based on the standardized ET/PET and EVI2 index extracted from the MODIS dataset during the growing season from May-August of 2000-2013.Furthermore, the fourteen-year summer monthly data for temperature, precipitation and moisture content of in situ measurements from sixteen meteorological stations were analyzed for four various land use areas (forest steppe, steppe, high mountains and desert steppe).The percentage deviations of climatological variables on the sixteen stations of four various land use areas were computed and compared to the MDSI anomaly during the study years.
The analysis of the MDSI anomalies revealed different temporal occurrences of dry (2001,2007,2009) and wet (2003,(2011)(2012)(2013) periods that dominated 35%-57% of the total area; while a maximum of 57% exhibited being the wettest in the year 2003, and 54% occurred as the driest in year 2007.Note that the year 2004 with higher precipitation and lower temperature was categorized as a normal year possibly because of "the memory effect" of higher soil moisture preserved in the ground in the year 2003.The exploration of the relationship between the MDSI anomaly performance with climatological anomalies indicated that the T anomaly had a negative correlation of ´0.66 (p < 0.01), whereas the P anomaly had a more strongly positive correlation of 0.81 (p < 0.001).Overall, the results showed that MDSI and T anomalies were negatively correlated, while MDSI and P anomalies were positively correlated, during the study period.As the MDSI anomaly distributions demonstrated, the highest (wettest) MDSI anomaly and lowest T anomaly values were seen in 2003, while the lowest (driest) MDSI anomaly and the highest T anomaly were observed in 2007.Besides, the correlations for the MDSI anomaly with T and P anomalies were analyzed among the four land use areas.MDSI and T anomalies indicated negative correlations, which were generally decreasing from forest steppe to desert steppe (r = ´0.67 to r = ´0.45).By contrast, MDSI and P anomalies had positive correlations (r = 0.67 to r = 0.83) with significance from p < 0.01-0.001except for the forest steppe.
Rainfall would increase with the surface soil moisture content, which in turn would result in a high MDSI anomaly.Furthermore, the MDSI anomaly and soil moisture content correlations ranged 0.63-0.74with significance p < 0.01 for various land use areas, except with a low value of 0.39 (p > 0.10) in the steppe during the study period.In addition, the responses of the MDSI anomaly to in situ measurements of precipitation and soil moisture strongly depend on the natural vegetation in each land use area.Comparisons across land use areas indicated that various features of natural plants are utilized in response to growing season hydrologic condition variations.Therefore, these results further proved that the MDSI anomaly can be used as an indicator through its correlation with soil moisture content to evaluate drought event.
The MDSI anomaly analysis and interpretations conducted here show a vital importance for both satellite-derived and observed anomalies.A comparison between the performance and limitations of MDSI and climatological variables' anomalies over Mongolia was presented.It demonstrated that the proposed algorithm can be applied for drought evaluation in the harsh climate regions.While demonstrating the usefulness of the MDSI for drought investigation, its performance against the existing drought indices, such as the Drought Severity Index [39], the Effective Drought Index [55] and the Bowen ratio [56], is worthy of further investigation.In additional, the use of the proposed index in the prediction of drought events, such as the four-dimensional water index [57], is practical and critical for improving the welfare of human beings.

Figure 2
Figure 2 shows the time series plots of the climatic variables during the growing season.The average air temperature was between 15.6 ˝C (year 2003) and 18.3 ˝C (2007); mean total precipitation ranged from 114 mm (2001)-204 mm (2013); and soil moisture content was from 5.3% (2000, 2002)-8.6%(2012) during the study period.Furthermore, the figure shows, based on meteorological long-term data, that the wettest years were 2003 and 2011-2013, while the driest years were 2001, 2007 and 2009; the rest of the years were normal years.Generally speaking, dry weather implies lower rainfall and higher temperature, so that drought occurred in the years 2001 and 2007.By contrast, wet weather accompanies higher precipitation and lower temperature, such as the years 2003, 2012 and 2013.One interesting question is raised though: why is the year 2004 categorized as a normal year?It is observed that the precipitation of 2004 was relatively low, but the soil moisture content and MDSI anomaly of 2004 was not low.This fact may indicate "the memory effect"; the high precipitation in 2003 possibly contributed to the higher soil moisture in 2004.Such a memory effect may prevail relatively easily since the ground is frozen in winter and preserves the soil moisture well till the next year.

Figure 2 .
Figure 2. Time series of the climate data for air temperature (°C: red open squares), precipitation (mm: blue bars) and soil moisture content (%: light red open triangles) averages over the sixteen stations during the growing season of 2000-2013.

Figure 2 .
Figure 2. Time series of the climate data for air temperature ( ˝C: red open squares), precipitation (mm: blue bars) and soil moisture content (%: light red open triangles) averages over the sixteen stations during the growing season of 2000-2013.

Figure 3
Figure 3 illustrates the growing season variations of the MDSI anomaly in the study area.The results represent the variations of the MDSI anomalies, following the total rainfall amount averaged over Mongolia, for the growing season from 2000-2013, as shown in Figure 2. The negative MDSI anomaly values that appeared in the south part of Mongolia in 2001 (Figure 3b), the southeast part in 2005 and 2006 (Figure 3f,g), the central area in 2007 (Figure3h) and the southwest region in 2009 (Figure3j) were less than ´0.5.In addition, the summers in 2003 and 2011-2013 (Figure3d,l-n) had more favorable hydrologic conditions than the other years (Figure3a,c,e,i,k) with fairly positive MDSI anomaly values of 0.5 or greater.The values of the MDSI anomaly were attributed to five categories (Table4), namely severely wet (>2), moderately wet (0.5 to 1.99), normal (0.49 to ´0.49), moderate drought (´0.5 to ´1.99) and severe drought (<´2), in this study.As illustrated in Figure4, the frequency distribution of the MDSI anomaly showed a wet period in the years 2003 and 2011-2013 with 85%-92% of the total area identified as normal and moderately wet, while 8%-15% of it was moderate drought conditions.In contrast, dry periods appeared in the years 2001, 2007 and 2009, and 45%-54% of the total area was identified as moderate drought, but 46%-55% of it was observed as normal and moderately wet conditions.In the other years, 18%-32% of the total area occurred as moderate drought, while more areas (68%-82%)

Figure 3 .
Figure 3. (a-n) MDSI anomaly distributions during the growing season for the years 2000-2013: The MDSI anomaly was computed as the standardized Z terms represented in Equation (6).

Figure 3 .
Figure 3. (a-n) MDSI anomaly distributions during the growing season for the years 2000-2013:The MDSI anomaly was computed as the standardized Z terms represented in Equation (6).

Figure 3 .
Figure 3. (a-n) MDSI anomaly distributions during the growing season for years 2000-2013: The MDSI anomaly was computed as the standardized Z terms represented in Equation (6).

Figure 5 .
Figure 5.Time series for the MDSI anomaly with different climatological variables over the sixteen stations during 2000-2013: (a) anomalies of MDSI (green open circles) and T (red opened squares); and (b) anomalies of MDSI (green open circles) and P (blue bars).The horizontal lines show that the MDSI anomaly value is −0.5 or less for drought conditions, whereas horizontal dashed lines show the percentage deviation values of T (0.0 or above) and P (0.0 or less) for drought condition.

Figure 5 .
Figure 5.Time series for the MDSI anomaly with different climatological variables over the sixteen stations during 2000-2013: (a) anomalies of MDSI (green open circles) and T (red opened squares); and (b) anomalies of MDSI (green open circles) and P (blue bars).The horizontal lines show that the MDSI anomaly value is ´0.5 or less for drought conditions, whereas horizontal dashed lines show the percentage deviation values of T (0.0 or above) and P (0.0 or less) for drought condition.
Remote Sens. 2016, 8, 509 11 of 16drier seasons is linked to rainfall, as well as the characteristics and topography of various land use areas.

Figure 6 .
Figure 6.Time series for the MDSI anomaly (green open circles) with soil moisture content (%: light red open triangles) over the sixteen stations in 2000-2013: the horizontal line shows that the MDSI anomaly value is −0.5 or less for drought condition.
, the forest steppe had the highest rainfall compared to the other land use areas.Soil moisture content reached roughly 11%-14%, while the amount of precipitation was about 274-324 mm in 2003, 2006, 2011 and 2012.In addition, the MDSI anomaly attains its maximum at about two (>1.0) in the years 2003 and 2011.The corresponding minimum MDSI anomaly was at −1.4,

Figure 6 .
Figure 6.Time series for the MDSI anomaly (green open circles) with soil moisture content (%: light red open triangles) over the sixteen stations in 2000-2013: the horizontal line shows that the MDSI anomaly value is ´0.5 or less for drought condition.

Figure 7
Figure7illustrates the temporal variations of the MDSI anomaly for four various land use areas and related to the climate conditions at the station scale.It represents that the MDSI anomaly exhibited high and low temporal variations in response to wet and dry with climatic parameters.In Figure7a, the forest steppe had the highest rainfall compared to the other land use areas.Soil moisture content reached roughly 11%-14%, while the amount of precipitation was about 274-324 mm in 2003, 2006, 2011 and 2012.In addition, the MDSI anomaly attains its maximum at about two (>1.0) in the years 2003 and 2011.The corresponding minimum MDSI anomaly was at ´1.4, while the minimum soil moisture content was about 4.6% in 2002.By contrast, the steppe (Figure7b) had a lower rainfall amount during the growing season.Both precipitation and soil moisture variables demonstrated a very similar trend with year in response to the MDSI anomaly.As an example, the minimum MDSI anomaly falls below ´1.5, while the minimum soil moisture content was about 3.6% in the year 2007.The high mountains and the desert steppe (Figure7c,d) had the lowest rainfall and stable variations as compared to the three parameters.In the high mountains, maximum rainfall ranged approximately 205-206 mm, while maximum the MDSI anomaly was between +0.7 and +1.9, and soil moisture content was about 7%-8% in 2003 and 2005.Particularly, the lowest rainfall (62 mm), the MDSI anomaly (´1.5) and soil moisture content (3.3%) occurred in 2009.However, in the desert steppe, the lowest amplitude of the MDSI anomaly (´0.9) and rainfall (25 mm and 62 mm) trends look very similar, such as in 2005 and 2006.The lowest rainfall reached around 25 mm in 2005, while the lowest soil moisture content was about 2.4% during 2004 and 2006.In addition, it is observed that, in 2002, the precipitations of the forest step and high mountain shown in Figure7a,c were the lowest, while the precipitation of steppe in Figure7bwas not low and even higher than those in2000,  2001, 2003 and 2004.This outstanding contrast possibly indicates that the regional dependence of precipitation prevails over the large-scale or climatic impacts.Remote Sens. 2016, 8, 509 12 of 16

Figure 7 .
Figure 7. (a-d) Temporal variations of the MDSI anomaly among different land use areas for 2000-2013: (a) forest steppe; (b) steppe; (c) high mountains; and (d) desert steppe.The MDSI anomaly symbols correspond to the average of four consecutive months' (the growing season) value calculated in the pixel region at each station.The vertical axes depict the standardized MDSI anomaly (green open circles) of the pixel region.Precipitation (mm: blue bars) is accumulated during the growing season, while soil moisture content (%: light red open triangles) is acquired during the growing season periods to match the MDSI anomaly value in the same periods.The horizontal line shows that the MDSI anomaly value −0.5 or less for drought condition.

Figure 7 .
Figure 7. (a-d) Temporal variations of the MDSI anomaly among different land use areas for 2000-2013: (a) forest steppe; (b) steppe; (c) high mountains; and (d) desert steppe.The MDSI anomaly symbols correspond to the average of four consecutive months' (the growing season) value calculated in the pixel region at each station.The vertical axes depict the standardized MDSI anomaly (green open circles) of the pixel region.Precipitation (mm: blue bars) is accumulated during the growing season, while soil moisture content (%: light red open triangles) is acquired during the growing season periods to match the MDSI anomaly value in the same periods.The horizontal line shows that the MDSI anomaly value ´0.5 or less for drought condition.

Table 1 .
Meteorological observation stations with information about the geographical location, elevation and various land use areas of Mongolia.

Table 2 .
Average and standard deviations (SD) of air temperature ( ˝C), precipitation (mm) and soil moisture content (%) over the sixteen stations divided into four various land use areas (forest steppe, steppe, high mountains and desert steppe) for the 2000-2013 growing seasons.
) during the study period.Furthermore, the figure shows, based on meteorological long-term data, that the wettest years were 2003 and 2011-2013, while the driest years were 2001, 2007 and 2009; the rest of the years were normal years.Generally speaking, dry weather implies lower rainfall and higher temperature, so that drought occurred in the years 2001 and 2007.By contrast, wet weather accompanies higher precipitation and lower temperature, such as the years 2003, 2012 and 2013.One interesting question is raised though: why is the year 2004 categorized as a normal year?It is observed that the precipitation of 2004 was relatively low, but the soil moisture content and MDSI anomaly of 2004 was not low.This fact may indicate "the memory effect"; the high precipitation in 2003 possibly contributed to the higher soil moisture in 2004.Such a memory effect may prevail relatively easily since the ground is frozen in winter and preserves the soil moisture well till the next year.

Table 2 .
Average and standard deviations (SD) of air temperature (°C), precipitation (mm) and soil moisture content (%) over the sixteen stations divided into four various land use areas (forest steppe, steppe, high mountains and desert steppe) for the 2000-2013 growing seasons.Land Use

Table 3 .
Correlation (r) matrix among the air temperature, precipitation and soil moisture content over the sixteen stations during the growing season of 2000-2013.Correlation was between two climatic parameters with a significance level of p < 0.001.

Table 4 .
Categories of wet (W), dry (D) and normal (WD) conditions of the MDSI anomaly.

Table 5 .
Correlations among the MDSI anomaly with different climatological variables over all of Mongolia during the growing season of 2000-2013.

Table 5 .
Correlations among the MDSI anomaly with different climatological variables over all of Mongolia during the growing season of 2000-2013.