Spatiotemporal Variations in Drought and Vegetation Response in Inner Mongolia from 1982 to 2019

: Drought events cause ecological problems, including reduced water resources and degraded vegetation. Quantifying vegetation responses to drought is essential for ecological management. However, in existing research, the response relationships (correlations and lags) were typically determined based on Pearson correlation coefﬁcient and the resultant lag times were constrained by the spatial and temporal resolutions of the analyzed data. Inner Mongolia is an important ecological barrier in northern China. Ecological security is one of the most concerned issues of the region’s sustainable development. Herein, we combined Global Inventory Modeling and Mapping Studies (GIMMS) normalized difference vegetation index (NDVI3g) with Systeme Probatoire d’Observation de la Terra-vegetation (SPOT-VGT) NDVI data through spatial downscaling. The obtained 1 km-resolution NDVI dataset spanning Inner Mongolia from 1982 to 2019 was used to represent the reﬁned vegetation distribution. The standardized precipitation evapotranspiration index (SPEI) derived from gridded meteorological data was used to measure drought over the same period. We investigated the spatiotemporal characteristics of vegetation and drought in the region in the past 38 years. We then discussed changes in different vegetation responses to drought across eastern Inner Mongolia using cross wavelet transform (XWT) and wavelet coherence (WTC). The results reveal that in 82.4% of the study area, NDVI exhibited rising trends, and the SPEI values exhibited declining trends in 78.5% of the area. In eastern Inner Mongolia, the grassland NDVI was positively correlated with SPEI and signiﬁcantly affected by drought events, while NDVI in forestlands, including shrubs, broad-leaved forests, and coniferous forests, was negatively correlated with SPEI in the short term and weakly affected by drought. The NDVI lag times behind SPEI in grasslands, coniferous forests, and broad-leaved forests were 1–1.5, 4.5, and 7–7.5 months, respectively. These ﬁndings provide a scientiﬁc foundation for environmental preservation in the region.


Introduction
Vegetation maintains water and soil conditions, regulates the climate, and controls the energy balance of terrestrial ecosystems which rely heavily on it for their survival [1]. Climate variations can lead to changes in vegetation which in turn is an "indicator" of global climate change [2]. It is therefore of great significance to quantify the spatial patterns and temporal dynamics of the important relationship between them [3].
According to the sixth assessment report released by the United Nations Intergovernmental Panel on Climate Change (IPCC) in 2021, the global average surface temperature in 2011-2020 has increased by 1.09 • C compared with pre-industrialization (1850-1900) [4]. In the context of global warming, the frequent occurrence of drought has a huge impact on the vegetation ecosystem [5]. Meteorological drought is a disaster caused by water shortages

NDVI Data
The GIMMS NDVI3g dataset is the longest available NDVI time series in the world [11]. It contains 816 images with a temporal resolution of 15 days covering the period from 1982 to 2015. However, because its spatial resolution is approximately 8 km, capturing detailed changes in vegetation using this dataset is difficult [11]. We downloaded these images from the National Aeronautics and Space Administration (NASA) (https://ecocast.arc.nasa.gov/data/pub/gimms/, accessed on 11 December 2021). The maximum value composite (MVC) method was used to regenerate a new image by extracting the maximum pixel value from several images. This method can eliminate the influences of clouds and the atmosphere [36]. In this work, MVC was applied to synthesize monthly data, and a total of 408 images were obtained.
The monthly 1 km-resolution SPOT-VGT NDVI dataset contains 252 images covering the 1999-2019 period. We downloaded these images from the Resource and Environmental Science and Data Center of the Chinese Academy of Sciences (https://www.resdc.cn/, accessed on 11 December 2021). The MVC was used to process the GIMMS NDVI3g and

Thematic Data
The National Qinghai-Tibet Plateau Science Data Center (https://data.tpdc.ac.cn/ zh-hans/, accessed on 11 December 2021) provided gridded meteorological data on the monthly scale from 1982 to 2019, including 1 km-resolution monthly temperature and precipitation data. These data were generated using a spatial downscaling scheme based on long-time-series, low-spatial-resolution (approximately 55 km) meteorological data from the Climate Research Unit (CRU). The downscaled results obtained for China were evaluated based on data recorded at 496 meteorological stations, and the monthly temperature and precipitation deviations were 0.82-1.28 • C and 13.3 mm, respectively [37]. The total number of downloaded monthly temperature and precipitation images was 912.
The vegetation type data (CCI-LC) came from the global land cover dataset of the European Space Agency, which was created during the second climate change action phase (http://maps.elie.ucl.ac.be/CCI/viewer/, accessed on 11 December 2021). The data covered a time period from 1992 to 2019, and the pixel size was 300 × 300 m. Figure  herein have a long-time span  and that vegetation types inevitably change, this paper adopted vegetation type data at four different times: 1992, 2000, 2010, and 2019.

Methodology
Based on the three types of datasets listed above, the flowchart of the methodology is shown in Figure 2. First, the GIMMS NDVI3g and SPOT-VGT NDVI datasets were combined to obtain an annual NDVI dataset with a 1 km resolution from 1982 to 2019, and the accuracy was verified (Section 3.1). SPEI was calculated based on the gridded meteorological data at the same time, including monthly temperature and precipitation data (Section 3.2). Second, a trend analysis was conducted to examine the change dynamics of vegetation and drought over the past 38 years (Section 3.3). Then, by calculating the correlation coefficient (r) between NDVI and SPEI for each pixel in the whole Inner Mongolia from 1982 to 2019, the association between vegetation changes and simultaneous drought was understood (Section 3.4). Finally, the correlation and hysteresis characteristics between the vegetation and drought of four vegetation types, including grasslands, shrubs, broad-leaved forests, and coniferous forests in eastern Inner Mongolia were explored using XWT and WTC wavelet analysis techniques (Section 3.5).
Remote Sens. 2022, 14, x FOR PEER REVIEW 6 of 25 Figure 2. Flowchart of the methods used in this study. In the figure, TR represents temporal resolution and SR represents spatial resolution. Parallelograms represent datasets and rectangles represent research methods.

Spatial Downscaling
The regression model method is a way to establish a linear or nonlinear relationship between forecast variables and predictors, and it is the earliest and most widely used statistical downscaling method [38]. With the help of the IDL programming language, the linear regression equation was established based on the overlapping SPOT-VGT NDVI and resampled GIMMS NDVI3g in 2000 to 2015 [39]. It can be expressed as follows: where GN i represents the GIMMS NDVI3g data in the ith year, SN i represents the SPOT-VGT NDVI data in the ith year, a and b are constants, and c i is the residual of the linear equation. Based on Equation (1) the 8 km-resolution GIMMS NDVI3g data spanning from 1982 Flowchart of the methods used in this study. In the figure, TR represents temporal resolution and SR represents spatial resolution. Parallelograms represent datasets and rectangles represent research methods.

Spatial Downscaling
The regression model method is a way to establish a linear or nonlinear relationship between forecast variables and predictors, and it is the earliest and most widely used statistical downscaling method [38]. With the help of the IDL programming language, the linear regression equation was established based on the overlapping SPOT-VGT NDVI and resampled GIMMS NDVI3g in 2000 to 2015 [39]. It can be expressed as follows: where GN i represents the GIMMS NDVI3g data in the ith year, SN i represents the SPOT-VGT NDVI data in the ith year, a and b are constants, and c i is the residual of the linear equation.
Based on Equation (1) the 8 km-resolution GIMMS NDVI3g data spanning from 1982 to 1999 were downscaled to 1 km-resolution data. Then, the downscaled NDVI data from 1999 was taken as an example to illustrate the downscaling accuracy because the SPOT-VGT NDVI data in 1999 did not participate in the linear regression analysis described above.

SPEI Calculation
Vicente-Serrano et al. [40] presented SPEI, which takes into account the influence of temperature on potential evapotranspiration, apart from precipitation, and improves the performance of the precipitation-based SPI. The calculation process is as follows [40,41].
First, based on the gridded monthly mean temperature from 1982 to 2019, the potential evapotranspiration (PET) was calculated by using the Thornthwaite approach [42].
where PET is the potential evapotranspiration, K is the correction coefficient of the latitude and month function, T i is the monthly average temperature, and H is the annual heat index.
A is a constant, determined by H, where Second, the difference between monthly precipitation and evapotranspiration was calculated: where D i is the difference between precipitation and evapotranspiration, P i is the monthly precipitation, and PET i is the monthly evapotranspiration. Third, the log-logistic probability distribution of three parameters was used to normalize the D i data sequence. The calculation formula of the log-logistic probability distribution function is as follows: where the α, β, and γ parameters are fitted by the linear moment method. Finally, the cumulative probability density was standardized: where ω is the cumulative probability function value of the evapotranspiration precipitation derivation function. c 0 = 2.515517, c 1 = 0.802853, c 2 = 0.010328, d 1 = 1.432788, d 2 = 0.189269, and d 3 = 0.001308. SPEI can be calculated at monthly, seasonal, annual, and biennial scales. Here, to assess annual drought conditions, the 12-month SPEI (SPEI-12) was calculated with the help of the Python programming language. According to the National Meteorological Drought Grade standard released in 2017 (GB/T20481-2017), the SPEI outputs were categorized as shown in Table 1. To identify the vegetation and drought trends in time and space, Theil-Sen median trend analysis was utilized with the help of R language. Theil-Sen median trend analysis is a reliable nonparametric approach for calculating statistical trends that is superior to regression analysis in terms of preventing error [43,44]. The calculation formula is as follows: where x j and x i are the NDVI or SPEI values for the jth year and the ith year, respectively. If the trend value (β) is greater than 0, it is considered to indicate a rise; otherwise, it indicates a decline.

Mann-Kendall Trend Test
The Mann-Kendall trend test is a nonparametric rank test method that is unaffected by missing values or outliers and widely used to judge the significance of long-time-series data trends [45,46].
var(s) = n(n − 1)(2n + 5) 18 (13) sign Where Z (−∞,+∞) is the statistic. It is significant when Z is larger than 1.96 or less than −1.96; otherwise, it is not significant. var(s) is the variance, n is the length of the time series, and sign is a symbolic function.

Pearson Correlation Analysis
To clarify the correlation between vegetation change and drought, we used the Pearson correlation coefficient to analyze the pixel-by-pixel correlation between vegetation NDVI and SPEI in Inner Mongolia from 1982 to 2019 with the help of R language, and used a t test to test the significance of the correlation coefficient. The Pearson correlation between NDVI and drought index indicates the impact of drought on vegetation [6].

Wavelet Analysis
Wavelet analysis is an effective tool for obtaining the time-frequency characteristics, including the periodic, correlation, and hysteresis characteristics, of two time series in the fields of meteorology, hydrology, and astronomy [47,48]. With the help of Python programming language, XWT and WTC, components of the wavelet analysis method, were used herein to identify the correlation and lag characteristics between the NDVI and SPEI sequences.
XWT can display the common power of the two sequences in the time-frequency domain [28,49]. Supposing that the two time series are represented by X ={x 1 , x 2 , x 3 , . . . , x n } and Y = y 1 , y 2 , y 3 , . . . , y n , the corresponding continuous wavelet transforms with frequency scale s can be denoted as W X n (s) and W Y n (s). The cross wavelet spectrum is therefore defined as follows: where W Y * n (s) is the complex conjugate of W Y n (s). The cross wavelet power spectrum is W XY n (s) . WTC may be applied to detect remarkable coherence regardless of power level [28]. Its phase spectrum reflects the correlation and lag between the two time series. The WTC spectrum is defined as follows: where s is the scale and S is the smoothing operator. The phase angle in the phase spectrum can reflect the lag characteristics of two time series in different time domains, and the correlation between the two time series can be analyzed by considering the sign of the phase angle [28]. The average phase angle is defined as follows: where a i is a single phase angle in the wavelet transform time domain and i is the phase angle label.
To explore the response relationships between different types of vegetation changes and drought using XWT and WTC, based on the four-phase vegetation type data mentioned in Section 2.2.2, this paper extracted four vegetation types with relatively large distribution areas in eastern Inner Mongolia (Figure 1b), including grassland, shrub, broadleaved forest, and coniferous forest. At the same time, the NDVI and SPEI datasets were resampled to 300 m using the nearest neighbor method in the resampling tool in ArcGIS 10.5 (https://desktop.arcgis.com/zh-cn/arcmap/10.5/tools/data-management-toolbox/ resample.htm, accessed on 11 December 2021), with the same spatial resolution as the vegetation type data, so as to facilitate the extraction of NDVI and SPEI sequences at different vegetation types. Then, XWT and WTC were applied to the NDVI and SPEI sequences of the four vegetation types.

Accuracy Verification of the Downscaled NDVI Data
A total of 1000 pixels were randomly selected to compare the downscaled data and the SPOT-VGT NDVI in 1999 ( Figure 3). The correlation coefficient was 0.98, indicating that the downscaled NDVI data were reliable.

Distribution of Vegetation
The 38 year mean NDVI calculated using the annual NDVI data from 1982 to 2019 was low in western Inner Mongolia and high in the eastern region, exhibiting obvious regional heterogeneity (Figure 4a). The maximum value in eastern Inner Mongolia was 0.93, distributed in the Daxinganling Mountains. The minimum value was 0, corresponding to areas where lakes were located. The maximum value in the central region was 0.89, distributed in the eastern part of Xilin Gol League, and the minimum value was 0.03, distributed in the western part of Xilin Gol League. The maximum value in the western region was 0.84, distributed in the Hetao Plain, and the minimum value was 0, found in the desert regions.
The classified mean NDVI values fall into four categories, as shown in Table 2 and Figure 4b [5]. A total of 33.9% of the entire area was covered by areas with high vegetation and mostly found in eastern Inner Mongolia. Medium vegetation cover occupied 27.9% of the total area, which was mostly found in the grasslands of central and eastern Inner Mongolia. The areas with low vegetation coverage accounted for 22.4% of the total area and were mostly found in the desert steppe areas. Areas with no vegetation coverage accounted for 15.8% and were predominantly located in the desert regions of the west.

Vegetation Characteristics in Inner Mongolia from 1982 to 2019 4.2.1. Distribution of Vegetation
The 38 year mean NDVI calculated using the annual NDVI data from 1982 to 2019 was low in western Inner Mongolia and high in the eastern region, exhibiting obvious regional heterogeneity (Figure 4a). The maximum value in eastern Inner Mongolia was 0.93, distributed in the Daxinganling Mountains. The minimum value was 0, corresponding to areas where lakes were located. The maximum value in the central region was 0.89, distributed in the eastern part of Xilin Gol League, and the minimum value was 0.03, distributed in the western part of Xilin Gol League. The maximum value in the western region was 0.84, distributed in the Hetao Plain, and the minimum value was 0, found in the desert regions.
The classified mean NDVI values fall into four categories, as shown in Table 2 and Figure 4b [5]. A total of 33.9% of the entire area was covered by areas with high vegetation and mostly found in eastern Inner Mongolia. Medium vegetation cover occupied 27.9% of the total area, which was mostly found in the grasslands of central and eastern Inner Mongolia. The areas with low vegetation coverage accounted for 22.4% of the total area and were mostly found in the desert steppe areas. Areas with no vegetation coverage accounted for 15.8% and were predominantly located in the desert regions of the west.

Vegetation Changes in Time and Space
The annual NDVI showed a rising trend, with an average value of 0.45 from 1982 to 2019, indicating that the vegetation status improved throughout this period ( Figure 5

Vegetation Changes in Time and Space
The annual NDVI showed a rising trend, with an average value of 0. 45     the whole, the vegetation coverage was the highest in summer and the lowest in winter. NDVI in summer and autumn in Inner Mongolia showed a weak rising trend. In spring and winter, NDVI showed a weak declining trend. The results of the Theil-Sen median trend analysis and Mann-Kendall test were superimposed to study the spatial vegetation variations. The classification and statistics of annual NDVI spatial change trends are shown in Table 3. Areas with rising trends in NDVI changes accounted for 82.4%, while areas with declining trends accounted for only 17.6%. The NDVI changes exhibited rising trends in most areas, corresponding to an improved vegetation status, and 45.1% of the total areas were significantly improved. The spatial distribution of the annual NDVI change trends ( Figure 7) showed that areas with rising NDVI (indicating improved vegetation) were distributed in each prefecture-level city. Obvious vegetation improvement was observed in the Hetao Plain and the eastern study area, while the central region exhibited slight vegetation improvement. The abundant water resources in the Hetao Plain and eastern Inner Mongolia are the reasons for the vegetation improvements. In addition, the execution of ecological construction projects, for instance, returning farmlands to forestlands and grasslands and restoring vegetation in mining areas, has greatly promoted the improvement in vegetation. The northwestern region of Alxa League in western Inner Mongolia, a portion of Xilin Gol League in central Inner Mongolia, and the Daxinganling Mountains to the east all had decreased NDVI, indicating that vegetation degradation had occurred. The Alxa League The results of the Theil-Sen median trend analysis and Mann-Kendall test were superimposed to study the spatial vegetation variations. The classification and statistics of annual NDVI spatial change trends are shown in Table 3. Areas with rising trends in NDVI changes accounted for 82.4%, while areas with declining trends accounted for only 17.6%. The NDVI changes exhibited rising trends in most areas, corresponding to an improved vegetation status, and 45.1% of the total areas were significantly improved. The spatial distribution of the annual NDVI change trends (Figure 7) showed that areas with rising NDVI (indicating improved vegetation) were distributed in each prefecturelevel city. The reasons for the vegetation degradation in the Daxinganling Mountains are likely due to its location as a transition zone from a plateau region to a plain region where the ecological environment is fragile, and anthropic activities, such as unreasonable logging measures, have occurred.  The spatial variation in vegetation in each season ( Figure 8) and the area proportion of NDVI change trend classification (Table 4) showed that from 1982 to 2019, the spatial distribution of NDVI change trends in Inner Mongolia were quite different in each season. The areas with rising trends of NDVI in spring accounted for 25.1% of the total area, mainly in Hulunbuir and the eastern Ordos. The declining areas were widely distributed, accounting for 74.9% of the total area. The spatial distribution of the change trend of summer NDVI was similar to that of annual NDVI. In autumn, the areas of NDVI with rising trends accounted for 58.1%, and the areas with declining trends accounted for 41.9%. The declining areas were mainly distributed in central and western Inner Mongolia. In winter, NDVI in most regions showed declining trends (83.9%), while the rising areas accounted for only 16.1%, mainly distributed in eastern and southern Inner Mongolia. On the whole, the change trend of NDVI in spring and winter was mainly declining, and it was mainly rising in summer and autumn.  (Table 4) showed that from 1982 to 2019, the spatial distribution of NDVI change trends in Inner Mongolia were quite different in each season. The areas with rising trends of NDVI in spring accounted for 25.1% of the total area, mainly in Hulunbuir and the eastern Ordos. The declining areas were widely distributed, accounting for 74.9% of the total area. The spatial distribution of the change trend of summer NDVI was similar to that of annual NDVI. In autumn, the areas of NDVI with rising trends accounted for 58.1%, and the areas with declining trends accounted for 41.9%. The declining areas were mainly distributed in central and western Inner Mongolia. In winter, NDVI in most regions showed declining trends (83.9%), while the rising areas accounted for only 16.1%, mainly distributed in eastern and southern Inner Mongolia. On the whole, the change trend of NDVI in spring and winter was mainly declining, and it was mainly rising in summer and autumn. Remote Sens. 2022, 14, x FOR PEER REVIEW 14 of 25  Inner Mongolia is a leading region with regard to the construction of a national ecological civilization. In recent years, the region has effectively implemented a series of  Inner Mongolia is a leading region with regard to the construction of a national ecological civilization. In recent years, the region has effectively implemented a series of ecological projects, such as preventing and controlling sand deposition, returning grazing lands to grasslands, and returning farmlands to forestlands. These projects have promoted an increase in vegetation coverage. However, overgrazing, reclamation, deforestation, and massive mining phenomena are ongoing in some areas. Natural disasters such as diseases, swarms of pests, and forest fires caused either directly or indirectly by the anthropic destruction of forest ecosystems are commonplace, thereby inducing vegetation degradation. Thus, the impacts of human activities should not be underestimated. ecological projects, such as preventing and controlling sand deposition, returning grazing lands to grasslands, and returning farmlands to forestlands. These projects have promoted an increase in vegetation coverage. However, overgrazing, reclamation, deforestation, and massive mining phenomena are ongoing in some areas. Natural disasters such as diseases, swarms of pests, and forest fires caused either directly or indirectly by the anthropic destruction of forest ecosystems are commonplace, thereby inducing vegetation degradation. Thus, the impacts of human activities should not be underestimated.

Drought Monitoring Results
Drought occurred when the SPEI values were less than −0.5. The annual drought (SPEI-12) fluctuations ( Figure 9) showed that the study area had a humid environment before 1999, followed by drought conditions. The drought conditions in Inner Mongolia showed an aggravating trend, especially after 2000, and the drought frequency increased significantly. During the same time period, Wang et al. [50], An et al. [51], and Qin et al.

Drought Changes in Time and Space
The maximum SPEI was 1.3 in 1990, and the minimum was −1.45 in 2000 ( Figure 10). In 2000, the total precipitation was 239.4 mm, and the average temperature was 3.7. This was the driest year from 1982 to 2019. The annual SPEI showed a declining-fluctuating trend, and the rate was −0.034 per year, meaning that the study area became increasingly arid. A drought event occurred in 1999, followed by droughts from

Drought Changes in Time and Space
The maximum SPEI was 1.3 in 1990, and the minimum was −1.45 in 2000 ( Figure 10). In 2000, the total precipitation was 239.4 mm, and the average temperature was 3.7. This was the driest year from 1982 to 2019. The annual SPEI showed a declining-fluctuating trend, and the rate was −0.034 per year, meaning that the study area became increasingly arid. A drought event occurred in 1999, followed by droughts from 2000 to 2002, from 2004 to 2007, and in 2009, 2011, and 2017. The drought frequency was approximately 29%.
The annual SPEI change trend and the significance distribution from 1982 to 2019 ( Figure 11) showed that the SPEI values ranged from −0.071 to 0.027. A total of 78.5% of the entire area was covered by areas with declining SPEI trends, indicating that most areas tended to be arid, among which 53.1% showed no significant declining trend; that is, the aridity was not obvious in these areas, while the remaining areas showed obvious aridification, including the Inner Mongolia Plateau and the western desert areas. The eastern edge of the study area, Ordos and its surrounding areas, accounted for 21.5% of the overall area with rising SPEI trends, although no portion of this area reached a significant level; that is, humidification was not obvious in any of these areas. The increasing rainfall in the east, the Yellow River water supply project in Ordos, and the implementation of artificial rain in Inner Mongolia all played positive roles in alleviating drought conditions. The annual SPEI change trend and the significance distribution from 1982 to 2019 ( Figure 11) showed that the SPEI values ranged from −0.071 to 0.027. A total of 78.5% of the entire area was covered by areas with declining SPEI trends, indicating that most areas tended to be arid, among which 53.1% showed no significant declining trend; that is, the aridity was not obvious in these areas, while the remaining areas showed obvious aridification, including the Inner Mongolia Plateau and the western desert areas. The eastern edge of the study area, Ordos and its surrounding areas, accounted for 21.5% of the overall area with rising SPEI trends, although no portion of this area reached a significant level; that is, humidification was not obvious in any of these areas. The increasing rainfall in the east, the Yellow River water supply project in Ordos, and the implementation of artificial rain in Inner Mongolia all played positive roles in alleviating drought conditions.   The annual SPEI change trend and the significance distribution from 1982 to 2019 ( Figure 11) showed that the SPEI values ranged from −0.071 to 0.027. A total of 78.5% of the entire area was covered by areas with declining SPEI trends, indicating that most areas tended to be arid, among which 53.1% showed no significant declining trend; that is, the aridity was not obvious in these areas, while the remaining areas showed obvious aridification, including the Inner Mongolia Plateau and the western desert areas. The eastern edge of the study area, Ordos and its surrounding areas, accounted for 21.5% of the overall area with rising SPEI trends, although no portion of this area reached a significant level; that is, humidification was not obvious in any of these areas. The increasing rainfall in the east, the Yellow River water supply project in Ordos, and the implementation of artificial rain in Inner Mongolia all played positive roles in alleviating drought conditions.

. Correlation between Vegetation Changes and Simultaneous Drought in Inner Mongolia
By calculating the correlation coefficient (r) between NDVI and SPEI for each pixel from 1982 to 2019, the association between vegetation changes and simultaneous drought was discovered. The significance of these correlations was then tested using the p value obtained from the t test, and the significance was separated into four categories (Table 5). According to the spatial distributions of the correlation coefficient results between NDVI and SPEI (Figure 12a), the range of r is −0.74-0.79. Positive correlation regions covered 63.6% of the research area, whereas negative correlation areas covered 36.4%, indicating that the vegetation status was well correlated with drought. According to Figure 12b, the correlations in 76.1% of the areas with positive correlations were not significant. Most of the areas with a significant correlation were located in the center of the study area, and they appear sporadically in the eastern cities. Grasslands can be found in abundance in these regions. Considering that grasslands have a low tolerance for drought, the occurrence of drought can seriously impact grasslands, causing NDVI to decline. The negatively correlated areas were distributed mainly in the Alxa League and Hetao Plain in western Inner Mongolia and part of Hulunbuir in eastern Inner Mongolia. Hulunbuir featured a diverse range of broad-leaved and coniferous forests, both of which had the ability to strongly resist drought. In addition, these areas are relatively cold, and low temperatures might have a considerable impact on vegetation [53]. The impact of human activities is also greater, whether positive (shelter forest projects) or negative (excessive logging). The Alxa League is a desert area with relatively sparse vegetation and little precipitation. Drought has no obvious impact on vegetation in the area. The vegetation in the Hetao Plain has improved mainly due to anthropic irrigation with water sourced from the Yellow River and ecological restoration projects. The results of Pearson correlation analysis in Section 4.4.1 showed a correlation between NDVI and SPEI, but large differences between different vegetation types were Gu et al. [34] analyzed the impact of drought on vegetation in Inner Mongolia from 1982 to 2015, and pointed out that NDVI was negatively correlated with drought in northern forest areas and western desert areas of Inner Mongolia, and positively correlated with drought in central grassland. The findings are consistent with the results of this study. Based on the meteorological data of 41 stations in Inner Mongolia and GIMMS NDVI, Yang and Yang [33] discussed the response of vegetation to drought from 1982 to 2013. The correlation results of NDVI and SPEI obtained in this study are similar to the above research results in terms of spatial distribution, but the area of positive correlation is 18% lower than it. The difference may be caused by the difference in data spatial resolution and research period.

Relationships between Different Vegetation Types and Drought in Eastern Inner Mongolia
The results of Pearson correlation analysis in Section 4.4.1 showed a correlation between NDVI and SPEI, but large differences between different vegetation types were found. Therefore, XWT and WTC were used to further explore the response relationships between different types of vegetation changes and drought.
In the cross wavelet power spectrum and wavelet coherent condensation spectrum of SPEI and NDVI (Figures 13 and 14), the arrows pointing to the right indicate positive correlations between SPEI and NDVI, while the arrows pointing to the left indicate negative correlations. The downward (upward) arrows indicate that SPEI changed 1/4 (3/4) of a cycle ahead of NDVI [54]. The average phase angle can be transformed into a lag time between SPEI and NDVI.
According to Figure 13a, the grassland SPEI and NDVI time series had a significant resonance period of 2.2-3.1 years from 2004 to 2008. The period was mainly in phase during these years. For shrubs, broad-leaved forests, and coniferous forests, the cross wavelet power spectra of SPEI and NDVI were similar and exhibited two significant resonance periods (Figure 13b-d). Figure 14a showed that SPEI had significant coherence with the grassland NDVI in the 5.5-8.3-year period from 2005 to 2011, with a notable positive correlation. The phase differences between SPEI and NDVI from 1985 to 2002 during the 2.3-4.9-year period were 30-45 degrees, meaning that the grassland NDVI lagged behind SPEI by 1-1.5 months. Figure 14b,d shows that SPEI were anti-phase with NDVI of shrubs and coniferous forests in the 2-2.7-year period from 1994 to 2000; that is, these values were significantly negatively correlated. The phase difference between SPEI and NDVI of coniferous forests from 2005 to 2014 in the 3.3-4.9-year period was 135 degrees; that is, NDVI of coniferous forests lagged behind SPEI by 4.5 months. Figure 14c showed that the phase differences between SPEI and NDVI of broad-leaved forests from 1995 to 2004 during the 0-3.1-year period were 210-225 degrees. In other words, NDVI of broad-leaved forests lagged behind SPEI by 7-7.5 months. The resonance periods of SPEI and NDVI of the three types of forestland were all short-lived, exhibiting short-term negative correlations and indicating that drought conditions have relatively little impact on forestlands. However, the phase relationship between the grassland SPEI and NDVI varied little, forming a stable positive correlation and showing that drought had a great influence on grasslands.
Liu et al. [32] analyzed the correlation between vegetation growth and drought in Inner Mongolia from 1998 to 2013, and the results showed that drought had a large impact on vegetation changes in Inner Mongolia, but less impact on eastern forests. When drought occurs, vegetation can reduce water loss by reducing respiration, photosynthesis, and growth rate, so as to adapt to drought stress [55]. Due to the different structures and functions of different vegetation types, their adaptability to the climatic environment is different [14]. Usually, grassland vegetation has small plants, and a shallow root system and is sensitive to external environmental disturbances. When drought reduces the effective water content of the soil, grassland vegetation absorb water with difficulty, and leaves will wilt or even die [56]. Forest vegetation has a large leaf area and deep and developed roots, enabling it to adapt to drought and improve the water storage capacity under drought stress [56]. In addition, the forestland in Inner Mongolia is mainly distributed in the northeastern region. Studies have shown that the temperature in this area has a greater impact on vegetation, and drought events have a small impact on vegetation [53]. Our results of drought impact on vegetation based on wavelet analysis are consistent with the above conclusions. and growth rate, so as to adapt to drought stress [55]. Due to the different structures and functions of different vegetation types, their adaptability to the climatic environment is different [14]. Usually, grassland vegetation has small plants, and a shallow root system and is sensitive to external environmental disturbances. When drought reduces the effective water content of the soil, grassland vegetation absorb water with difficulty, and leaves will wilt or even die [56]. Forest vegetation has a large leaf area and deep and developed roots, enabling it to adapt to drought and improve the water storage capacity under drought stress [56]. In addition, the forestland in Inner Mongolia is mainly distributed in the northeastern region. Studies have shown that the temperature in this area has a greater impact on vegetation, and drought events have a small impact on vegetation [53]. Our results of drought impact on vegetation based on wavelet analysis are consistent with the above conclusions.

Anthropogenic Factors on Forestland in Eastern Inner Mongolia from 1982 to 2019
The results of wavelet analysis showed that in eastern Mongolia, NDVI and SPEI of forestland were negatively correlated. In order to further understand the reasons, we counted the average values of NDVI and SPEI of the pixels where the three types of forestland are located in the eastern part of Inner Mongolia, and made a line graph of changes (Figures 15 and 16).
As shown in the gray shadows in Figures 15 and 16, the forestland NDVI showed a declining trend from 1990 to 1999, when the areas were humidified. In the 1990s, excessive logging and land reclamation in the forested areas of the Daxinganling Mountains caused a considerable loss of forestland resources. Therefore, the main reason for declines in NDVI was the destruction of forestlands by humans. Considering the adverse consequences of overcutting and land reclamation, the first stage of the natural forest protection project was fully completed from 2000 to 2010 [57], as shown by the blue shadows in Figures 15 and 16. During this period, NDVI of forestlands showed a rising trend. Tong et al. [58] pointed out that since 2000, the project of converting farmland to forest in Inner Mongolia has been implemented on a large scale. However, our study found that most of the years were under dry conditions with SPEI values of less than -0.5. Such climatic conditions are not conducive to vegetation growth. However, the vegetation growth from 2000 to 2010 was relatively good, indicating that the national restoration project policy brought some improvements in vegetation coverage. The wavelet analysis results in Section 4.4.2

Anthropogenic Factors on Forestland in Eastern Inner Mongolia from 1982 to 2019
The results of wavelet analysis showed that in eastern Mongolia, NDVI and SPEI of forestland were negatively correlated. In order to further understand the reasons, we counted the average values of NDVI and SPEI of the pixels where the three types of forestland are located in the eastern part of Inner Mongolia, and made a line graph of changes (Figures 15 and 16).
As shown in the gray shadows in Figures 15 and 16, the forestland NDVI showed a declining trend from 1990 to 1999, when the areas were humidified. In the 1990s, excessive logging and land reclamation in the forested areas of the Daxinganling Mountains caused a considerable loss of forestland resources. Therefore, the main reason for declines in NDVI was the destruction of forestlands by humans. Considering the adverse consequences of overcutting and land reclamation, the first stage of the natural forest protection project was fully completed from 2000 to 2010 [57], as shown by the blue shadows in Figures 15  and 16. During this period, NDVI of forestlands showed a rising trend. Tong et al. [58] pointed out that since 2000, the project of converting farmland to forest in Inner Mongolia has been implemented on a large scale. However, our study found that most of the years were under dry conditions with SPEI values of less than -0.5. Such climatic conditions are not conducive to vegetation growth. However, the vegetation growth from 2000 to 2010 was relatively good, indicating that the national restoration project policy brought some improvements in vegetation coverage. The wavelet analysis results in Section 4.4.2 also revealed that NDVI and SPEI of forestlands were negatively correlated for a short period, mainly due to human activities (positive or negative). In the study of Sun et al. [59], it was pointed out that NDVI and SPEI in the eastern Inner Mongolia were negatively correlated, mainly due to human activities, such as deforestation, man-made fires, and ecological engineering, etc. It is consistent with the findings of this paper.

Study Limitations
Compared with previous studies, the research period of this paper is longer, the data resolution is higher, and the lag time obtained by time-frequency technology (wavelet analysis) is more accurate. These results can be used as a basis for formulating an ecosystem construction program in Inner Mongolia. However, there are also some limitations.
Firstly, the land surface temperature (LST) affects the energy distribution between the ground and vegetation and determines the surface air temperature [60]. LST is an im-

Study Limitations
Compared with previous studies, the research period of this paper is longer, the data resolution is higher, and the lag time obtained by time-frequency technology (wavelet analysis) is more accurate. These results can be used as a basis for formulating an ecosystem construction program in Inner Mongolia. However, there are also some limitations.
Firstly, the land surface temperature (LST) affects the energy distribution between the ground and vegetation and determines the surface air temperature [60]. LST is an important factor in global climate change and vegetation growth. Therefore, we should fur-

Study Limitations
Compared with previous studies, the research period of this paper is longer, the data resolution is higher, and the lag time obtained by time-frequency technology (wavelet analysis) is more accurate. These results can be used as a basis for formulating an ecosystem construction program in Inner Mongolia. However, there are also some limitations.
Firstly, the land surface temperature (LST) affects the energy distribution between the ground and vegetation and determines the surface air temperature [60]. LST is an important factor in global climate change and vegetation growth. Therefore, we should further study the effect of drought on vegetation in combination with surface temperature. Second, we only considered the effects of drought conditions on vegetation due to changes in precipitation and temperature. However, in addition to precipitation and temperature, soil, terrain, human activities, and other factors will also affect vegetation change [61]. According to the results of this study, human activities greatly impact the vegetation in Inner Mongolia. In follow-up research, we should further explore its impact on the vegetation change in combination with multiple influencing factors, especially human activities.

Conclusions
This study presents the spatiotemporal characteristics of vegetation and drought in Inner Mongolia from 1982 to 2019. The findings show that NDVI had a rising change trend, showing an improved vegetation status in 82.4% of the study area. The Hetao Plain and the eastern study area had the most obvious vegetation improvement, whereas the middle region had only a slight improvement. SPEI showed declining trends tending toward aridification in 78.5% of the study area. The Inner Mongolia Plateau and the western desert region were the primary areas that witnessed significant aridification.
Based on the relationships between different vegetation and drought in eastern Inner Mongolia, the NDVI and SPEI values of grasslands had a significant positive correlation, while the NDVI values of forestlands (shrubs, broad-leaved forests, and coniferous forests) were negatively correlated with the SPEI values during some short periods of time. These short-period fluctuations in forestland NDVI were related to human activities. XWT and WTC results show that grassland NDVI lagged behind SPEI by 1-1.5 months, coniferous forest NDVI lagged behind SPEI by 4.5 months, and broad-leaved forest NDVI lagged behind SPEI by 7-7.5 months.