Assessing the Drought Variability in Northeast China over Multiple Temporal and Spatial Scales

: Long-term drought variation provides a scientiﬁc foundation for water resource planning and drought mitigation. However, the spatiotemporal variation characteristics of drought in northeast China (NEC) are unclear. We conducted a comprehensive assessment of drought status and trends based on the Standardized Precipitation Evapotranspiration Index (SPEI) in NEC from 1990 until 2018. The ﬁndings show that: (1) the drying trend peaked in 2001, and then exhibited a mitigation tendency before drying again after 2013. The implementation of ecological restoration projects is primarily responsible for drought mitigation. (2) The areas with wetting and drying trends in the future would cover 86% and 17% of NEC, respectively. (3) There is a time lag between improved vegetation and the trend shift from dry to wet. (4) Spring and winter revealed wet trends within 71% and 84% of NEC, respectively, showing high sensitivity and resilience to drought, while 92–93% of NEC displayed dry tendencies during the summer and autumn seasons. The drought-affected area was the highest in summer and lowest in autumn. (5) The interannual drought severity was highest in May and June. (6) The highest drought impacts and trends occur within shrub and grass and sparsely vegetated land, as well as middle-temperate semiarid regions (M-semiarid). (7) The warmer the temperature zone, the more sensitive it is towards drought under the same hydrological conditions, showing a high drought-affected area. The drier the land, the higher the drought-affected area within the same temperature zone, with pronounced drought trends during the spring and summer seasons. Our ﬁndings highlight the need for the government to more explicitly develop drought mitigation strategies in accordance with NEC’s spatiotemporal drought variations and speciﬁcally the need to concentrate on droughts in M-semiarid regions occurring in summer, particularly in May and June.


Introduction
Drought is a fundamental and recurring extreme climate event over terrestrial land that has severe direct and indirect impacts on the ecological, agricultural, and socioeconomic sectors [1][2][3]. Global drylands have expanded under a warming climate [4,5] and high greenhouse gas emission scenarios when the climate is expected to be warmer than it is now. Most land will become drier by 2100, with the exception of India and northern tropical Africa [4]. Observations from the Climatic Research Unit (CRU) have also revealed a significant positive drying trend in land areas under meteorological drought [6]. Meteorological drought is defined by precipitation deficits over a region for an extended period of time, possibly combined with increased potential evapotranspiration [1], can last for a wide range of time scales, and is frequently spatially variable due to landscape and climatic heterogeneity. The effects of different temporal scales on the efficacy of drought indices are well documented [7][8][9][10]. For example, severe droughts, even if only for a short period, would have a disastrous impact on agriculture during crop growth stages [11] and increase the risk of fire. In subsequent years after the severe drought, changes in hydro-meteorological variables are observable. The increase of water demand exacerbates the effects of future droughts [12,13]. Moreover, light and moderate droughts with long durations would have devastating impacts on ecosystems and the water supply [14]. As such, improved observational quality is required for the efficient monitoring and management of droughts in many regions of the world, due to potential high drought impacts, combined with complicated drought characteristics.
With regard to China, although the aridity trend has slowed over the past three decades, it has become more pronounced in the southwest and northeast of China [15]. Northeast China (NEC) is a typically vulnerable region to climate change and has been described as a major ecological barrier for northern China and even northeast Asia, playing a significant role in biodiversity conservation and regulating the regional climate of northeast Asia [16,17]. Previous research has shown that NEC is one of the hotspots of warming areas in East Asia [18,19], and has a high incidence of drought disasters [20]. Almost 45.4% of the ecosystem is at moderate or high risk of meteorological drought [21]. To date, there have been many studies on drought index trends and drought characteristics [22][23][24]. However, little attention has been paid to investigating the spatiotemporal mutation points of drought trends and the driving factors. The spatiotemporal characteristics of drought types that may have serious consequences for vegetation growth and ecosystems have also not been clarified. In addition, NEC has abundant vegetation resources and diverse climate regions characterized by different hydrothermal conditions. Revealing the drought status and trends of different land cover types and climate regions is beneficial for local governments to accurately monitor and mitigate the impact of drought by taking appropriate measures.
The atmosphere's evaporative demand is expected to rise as global temperatures rise [25]. Thus, univariate risk assessment methods based solely on rainfall (i.e., the Standardized Precipitation Index (SPI), China Z Index (CZI)) may understate the risk of extreme drought events [26,27]. The Standardized Precipitation Evapotranspiration Index (SPEI) is calculated by normalizing the difference between precipitation and potential evapotranspiration based on a specific distribution, and is more sensitive to the drought severity caused by rapid temperature rises [28]. Moreover, unlike the Palmer Drought Severity Index (PDSI), where the identification of rapidly emerging droughts is limited due to the fixed time scale of 9 months, SPEI provides multiple time scales (1-48 months) [28]. Therefore, it facilitates the quantitative comparison of drought incidence over time and space [29] and thus has been widely used to quantify drought onset, duration, magnitude, and spatial extent globally [30][31][32]. However, high spatial resolution SPEI datasets have not yet been developed in NEC.
Here, SPEI-aggregated multi-time scales with a spatial resolution of 1 km were calculated to explore the annual and quarterly trends and their mutation points from 1990 to 2018. The future dry-wet trend was also investigated. Additionally, we evaluated the land area dynamics of monthly severe drought, quarterly moderate drought, and annual light drought during the study period. The drought area and trend of different climate regions and land cover types were further clarified. We aimed to explore the reliable spatiotemporal drought variation characteristics in NEC, which has important practical implications for coping with future drought risk and developing more detailed drought mitigation strategies to protect the ecological environment.  (Figure 1a). It is an important grain-producing and ecological conservation area in China [17]. The topography is characterized by central plains surrounded by mountain ranges. The mountain ranges are mainly composed of the Greater Khingan Mountains in the west, the Lesser Khingan Mountains in the north, and the Changbai Mountains in the east (Figure 1b). NEC is dominated by a temperate continental monsoon climate with hot and rainy summers, as well as cold and dry winters. The annual average precipitation ranges from 300 to 1000 mm, and the annual average temperature varies between −3 and 10 • C, with an uneven spatial distribution. From south to north, NEC spans the warm, middle, and cold temperate zones. From east to west, NEC broadly transitions from humid to subhumid, and semiarid areas. These various hydrothermal conditions result in a variety of climatic regions in NEC. Using the classification criteria of temperature zones and arid/humid regions, NEC can be divided into six climatic regions under the scheme of climate regionalization for the period 1981-2010 [33], which updates the original climatic regionalization. The climate regions in NEC include the cold temperate humid (C-humid, I), middle temperate humid (M-humid, II), warm temperate humid (W-humid, III), middle temperate subhumid (M-subhumid, IV), warm temperate subhumid (W-subhumid, V), and middle temperate semiarid (M-semiarid, VI) regions. The main characteristics of climate regions in NEC are shown in Table 1.  The SPEI calculated over a specific timescale represents the cumulative water balance over the total months of the specified time period (e.g., SPEI12 represents the 12-month cumulative water balance) [8,28]. To reflect the monthly, quarterly, and annual drought variation, we calculated SPEI for 1 month (SPEI01), 3 months (SPEI03), and 12 months (SPEI12) using the monthly mean precipitation and potential evapotranspiration data from 1990 to 2018. SPEI decreases indicate drier conditions, and vice versa. A drought occurs when the SPEI value is less than minus one standard deviation. A comprehensive explanation on the SPEI calculation procedure can be found in Vicente-Serrano et al. [28]. The monthly mean precipitation and potential evapotranspiration data at 1 km spatial resolution from 1990 to 2018 were obtained from the Loess Plateau SubCenter, National Earth System Science Data Center, National Science and Technology Infrastructure of China (http://loess.geodata.cn, accessed on 20 March 2021). The monthly mean precipitation data that we downloaded were spatially downscaled to 1 km from the 30 Climatic Research Unit dataset with the climatology dataset of WorldClim using the delta spatial downscaling method [34]. They were evaluated using observations collected from 496 weather stations across China for the period 1951-2016, and the verification results are credible [34]. The monthly average potential evapotranspiration data that we downloaded were calculated from air temperature data using the Hargreaves method [35].

Materials and Methods
The drought categories of the SPEI are shown in Table 2. In this study, we extracted the land area of three different drought types that may have detrimental impacts on ecosystems and vegetation growth: monthly severe drought (SPEI01 < −1.5), quarterly moderate drought (SPEI03 < −1), and annual light drought (SPEI12 < −0.5). The drought area of different climate regions and land cover types were further investigated using annual light drought as an example.

Trend Analysis and Significance Testing
The Theil-Sen (TS) estimator [36] and Mann-Kendall (MK) test [37][38][39], which are commonly used to detect trends in hydro-meteorological time series, were performed in this study to examine the magnitude of raster-based annual and quarterly drought trends and their significance. These methods do not require the data to conform to a normal distribution, are insensitive to outliers, and can objectively characterize the overall trend of change. Z values that are positive or negative indicate an increasing or decreasing tendency, respectively. Furthermore, the absolute value of Z > 1.96 indicates a statistically significant trend at the 95% confidence level.
Linear regression was then applied to investigate the quarterly temporal trends based on SPEI03. Meanwhile, the 5-year moving average method was applied to smooth time series fluctuations and to highlight trends.

Mutation Point Detection
The mutation points of the raster-based annual dry-wet trend in NEC from 1990 to 2018 were extracted using piecewise linear regression. The time series is fitted before and after the trend turns assuming that the shortest period for the trend to turn is 5 years. The potential mutation point is determined by the point with the smallest sum of residual squares. It is then assessed whether the trend before and after the mutation point passes the significance test (p < 0.1). If both pass the significance test, the potential mutation point becomes the actual mutation point.
Furthermore, the sequential MK test method was used to determine abrupt changes in the temporal quarterly dry-wet trend. As a result, a forward (UF) and backward (UB) curve will be formed. The time series statistical curve is denoted by UF, and the inverse time series statistical curve is denoted by UB. If UF > 0, the data sequence exhibits an upward trend, and vice versus. The significance level of 0.05 corresponds to a threshold value of ±1.96 for the UF and UB curves. A UB crossing the threshold indicates that the series has a significant trend over this period. The mutation point is determined as the time corresponding to the intersection of the two curves within the threshold [40].

Hurst Exponent
Hurst exponent is an effective statistical method for quantitatively assessing the durability of changes in time series data over long periods and is widely applied in climatology and vegetation studies [41][42][43]. In this study, the Hurst exponent based on the rescaled range (R/S) analysis method was used to calculate the Hurst exponent of SPEI12 from 1990 to 2018, and then to study the spatial distribution and persistence of SPEI12 change trend in the future. The complete calculation procedure is shown in Hurst [41] as well as Mandelbrot and Wallis [42]. Hurst exponent values ranging from 0 to 1 can be classified into three categories: if 0 < H < 0.5, the time series has anti-persistence, and the sequence has an opposite trend in the future. When H equals 0.5, changes in the time series are stochastic, and the trend in the past has no effects on the future. If 0.5 < H < 1, the time series has persistence, and the future trend will be consistent with the changing trend observed during the study period. The closer H is to 0, the more inconsistent the time series; conversely, the closer H is to 1, the more consistent the sequence.

Spatiotemporal Dynamics of Dry-Wet Conditions in Northeast China
To show the spatial variation of annual dry-wet conditions, the study period was divided into six phases at 5-year intervals. The first phase was the wettest in NEC, with light wet dominating (Figure 2a). NEC transitioned from light wet to near-normal, and then to light and moderate drought from the first to the third phase (Figure 2a-c). The region of light and moderate drought accounted for 78.51% in the third phase, which was the highest compared to the other phases, with moderate drought primarily occurring in the Songliao Plain's midsection (Figure 2c, Table 3). Drought conditions eased in the fourth stage compared to the previous phase, with light and moderate drought covering 32.31% of the total area, primarily in the west and north (Figure 2d). The drought was clearly relieved from the fourth to the fifth phases, with NEC exhibiting near-normal and light wet conditions (Figure 2d,e). In the last phase, the drought worsened once more, with light and moderate drought covering 6.39% of the total area ( Figure 2f).    NEC's quarterly moderate drought area varied by season each year (Figure 4a). On the timeline of the study period, spring and winter exhibited the earliest drought events, followed by summer, and the last occurred in autumn. When the proportion of drought area for each season during the study period was superimposed, we discovered that summer had the largest proportion, followed by spring and winter, and autumn had the lowest. Moreover, there was generally only one season per year when the proportion of drought area exceeded 50%, but there were two (summer and autumn) and three (summer, autumn, and winter) seasons in 2001 and 2007, respectively.
Similarly, the monthly extreme drought area increased, declined, and then increased again over the last three decades (Figure 4b). Drought conditions with a severe drought percent greater than 40% were most prevalent throughout 2000-2009, notably in May and June (Figure 4b). We could further differentiate drought areas in percent in different climatic regions and land types using the annual light drought as an example ( Figure 5). The warmer the temperature zone was, the more sensitive it was towards drought under the same hydrological conditions. Within the same temperature zone, the drier the land was, the higher the proportion of land experiencing drought was. Forests contained the smallest share of drought area compared to other land cover types, followed by croplands and wetlands, shrub and grass, and sparsely vegetated lands showing the highest droughtaffected area.

Annual Dry-Wet Trends and Mutation Points
Spatial patterns of SPEI12 trends can roughly be characterized through their east-west contrasts (Figure 6a). Wetting in the eastern region and drying in the western region were found to be the most significant. A decreasing trend was present in 86.3% of the entire area. Only 13.7% of NEC, mostly concentrated in Sanjiang Plain and its southern mountains, exhibited an increasing tendency. Additionally, the Hurst exponent, which predicts the trend in future dry/wet conditions, was derived based on SPEI12 data from 1990 to 2018. In the entire NEC, the Hurst exponent was less than 0.5 (Figure 6b). The dry-wet trend in the future is therefore inconsistent with the current situation. Sen's slope analysis and the Hurst exponent results were superimposed to determine the spatial distribution of forthcoming dry-wet trends (Figure 6b). The areas with increased and decreased trends would cover 86.3% and 16.7% of NEC in the future, respectively.

Quarterly Dry-Wet Trends and Mutation Points
The SPEI of spring and winter increased at a linear rate of 0.1/10a and 0.3/10a from 1990 to 2018, respectively, showing a wetting trend (Figure 8a,d). The SPEI of summer and autumn decreased at a linear rate of −0.2/10a and −0.1/10a, respectively, showing a drying trend (Figure 8b,c). From the UF curve, it can also be seen that the SPEI in spring and winter showed an increasing trend, while that in autumn and winter showed a decreasing trend (Figure 9). Between the critical line at the 0.05 significance level, the UF and UB curves for spring and winter intersected in 2003-2004, which is the onset of the mutation. After 2004, the UF curve rose obviously, but did not exceed the critical line (±1.96) (Figure 9a,d). The UF and UB curves in summer and autumn intersected in 1991 and 1994-1995, respectively, and exceeded the critical line (−1.96) in 2000-2001, with a significant decreasing trend (Figure 9b,c).  The regions with increasing SPEI in spring and winter occupied 71.1% and 83.6% of the entire area, respectively (Figure 10a,d). Moreover, 1.7% of NEC experienced a significant dry trend in spring, mostly in the southwest of NEC, whereas the eastern part of NEC experienced a notable wet trend (Figure 10a). The areas tending to be dry in winter were primarily found in the northern part of NEC, whereas 22.6% of the total area became significantly wetter, primarily in the central and eastern parts of NEC (Figure 10d). In addition, 92.6% and 92.1% of the total area displayed a drying trend in summer and autumn, respectively (Figure 10b,c). NEC demonstrated a transition from a dry trend to a wet trend from west to east during the summer (Figure 10b). In contrast to the wet trend in summer observed in the Sanjiang Plain, the wetter regions in autumn were primarily found in the Songliao Plain's central region and the Greater Khingan Mountain's southern region (Figure 10c). For various climate regions and land cover types, we extracted the average annual and quarterly SPEI trends ( Figure 11). Overall, a strong drought trend was visible within shrub and grass and sparsely vegetated lands (mainly in M-semiarid), followed by wetlands, forests, and croplands. For the climate regions, M-semiarid, W-subhumid, and C-humid (primarily forest) showed the highest drought trend. Specifically, spring showed a wet trend except for W-subhumid and M-semiarid. In winter, all climate regions apart from C-humid showed a wet trend. All climate regions and land cover types revealed annual and seasonal drying trends in summer and autumn. In the same temperature zone, the drier the land, the stronger the drought trend in spring and summer.

Driving Factors of Drought Dynamics and Mutation Points
NEC has exhibited a general dry tendency over the past three decades, as evidenced by the percentage of annual drought areas increasing by 0.7% per year and 86.3% of the area exhibiting a decreasing trend of SPEI12 (Figures 3b and 6). It is noteworthy that this trend is not constant throughout time and space. Diagnosing mutation points of drought dynamics and exploring their drivers are important for developing strategies to address climate risks. The spatial distribution of SPEI12 at 5-year intervals showed a gradual drying of NEC that peaked in 2000-2004, and then decreased and increased again in 2014-2018. (Figure 2). Additionally, from 2002 to 2013, the annual, quarterly, and monthly drought areas all displayed a tendency for mitigation (Figures 3b and 4). In this study, variations in potential evapotranspiration and precipitation can be directly linked to variations in drought conditions, and temperature has a significant impact on potential evapotranspiration. Pearson correlation analysis ( Figure 12) revealed that temperature and potential evapotranspiration correlated negatively with SPEI12 showing values of −0.29 and −0.80, respectively, and precipitation correlated positively (0.93). More deeply, changes in climate factors are often closely linked to human activities. The enormous demand for food and economic growth since the 1950s has resulted in natural vegetation degradation, the loss of wetlands, and prominent desertification and salinization in NEC [44,45], which may have contributed significantly to the continued drought trend from 1990 to 2001. In response to the massive deforestation and associated ecological issues, NEC implemented the Natural Forest Conservation Project and the Grain to Green Project in 1998 and 1999, respectively. Commercial logging in NEC's natural forests came to an end with the implementation of these two projects. On the other hand, studies have reported trends in greening and productivity growth in NEC [46,47]. The vegetation impacted by human activities had a degradation trend from 1990 to 2001 and an increase after 2001, with 93% of this improved trend being linked to human activity [48]. Increasing vegetation cover not only increases soil organic carbon and total nitrogen, but also enhances ecosystem services, lessens soil erosion, and modifies regional climate [49,50]. This is supported by the improvement in climatic conditions from 2002 to 2013 [48] and the decreased trend of drought area at multiple time scales (Figures 3b and 4). However, most of the planted trees that are utilized for revegetation are still young, have underdeveloped root systems, and can only draw water from surface soil. The consequent large increase in vegetation greenness could lead to an excessive loss of water from the surface soil. As the trees mature, their deeper root systems may cause excessive evaporation of deep soil water, which would limit the vegetation growth and impact the regional climate [51][52][53]. This might be connected to the drought's resurgence after 2014. Moreover, greenhouse gas (mostly CO 2 ) emissions, nitrogen deposition, and other adverse human activities are also important causes of warming and drying of the climate [32,54,55].
We identified the moment and region in NEC where the trend of SPEI12 turned (p < 0.01). In 1999-2000, the western Liaohe Plain (dominated by grassland, cropland, and sparsely vegetated land) and the eastern Lesser Khingan Mountain (dominated by forest) underwent a dry to wet transition. Similarly, the western Sanjiang Plain (dominated by forest and cropland) experienced the same shift in 2004 (Figure 7). These regions contain rich natural resources and fertile black soil, and relief from the drought contributes to improved food production and ecosystem functions and enriches biodiversity. The turning point of vegetation trends in the western Liaohe Plain, eastern Lesser Khingan Mountain, and western Sanjiang Plain occurred in 1999, 1998, and 2000, respectively (Figure 13), which is the time when the ecological restoration project started to be implemented. According to Zhang et al. [56], the implementation of ecological restoration projects considerably promoted the growth of vegetation cover in the Horqin sand (mostly in the western Liaohe Plain) and prevented the spread of sandy areas. Mao et al. [46] also reported that the natural ecosystems in the western NEC were greatly improved in terms of sandstorm prevention, habitat provision, and food production between 2000 and 2015. The time lag between vegetation improvement and drought mitigation points to a delayed impact of vegetation on promoting the climate.

Quarterly and Monthly Drought Differentiation Characteristics
On the timeline of the study period, spring and winter revealed the earliest drought occurrences (Figure 4a), and wet trends were seen in 71.1% and 83.6% of NEC ( Figure 10), respectively, implying a high sensitivity and resilience to drought in spring and winter. Areas with a drought trend in summer and autumn represented 92.6% and 92.1% of NEC, respectively ( Figure 10). Nevertheless, in contrast to other seasons, the proportion of drought area was highest in summer and lowest in autumn during the study period (Figure 4a), suggesting that the summer drought is the worst and that, despite the drought trend being widely distributed in autumn, the drought area is the lowest. This can be explained by the fact that plants enter the peak growth period in summer and the evapotranspiration rate is at its highest. Even though there is a certain amount of rainfall (average precipitation within summer from 1990 until 2018 in NEC: 398-2190 mm; data not shown), soil moisture still decreases dramatically. Autumn is characterized by a slowdown in plant growth and a reduction in soil water usage, yet despite this, soil moisture still shows a decreasing trend, since precipitation was less than evapotranspiration. Zhang et al. [57] also reported the persistent summertime drought in northern China since 1990 under anthropogenic warming. In addition, monthly severe drought areas of a percentage greater than 40% were most prevalent in May and June (Figure 4b). This implies that a focus should be laid on the severe droughts in May and June, which are a critical period for vegetation growth, and a response is necessary to develop appropriate drought mitigation measures (e.g., irrigation, artificial rainfall, etc.) in a timely manner, adapted to different vegetation types.

Spatial Variation of Drought in Northeast China
We provided a spatial meteorological perspective on which regions in NEC have more severe drought conditions (area and trends) and raised awareness on different climate regions and land types facing drought risk. Overall, shrub and grass areas and sparsely vegetated land (mainly in M-semiarid) indicate the highest drought trend, followed by wetlands, croplands, and forests (mainly in humid land) ( Figure 11). The proportion of drought areas also indicated that shrub and grass and sparsely vegetated lands had the highest percentage, followed by cropland and wetland, and forests showed the lowest percentages (Figure 5b).
Drought profoundly affects ecohydrological processes along the soil-plant-atmosphere continuum on timescales ranging from hours to years [58], thus affecting both the functioning and structure of ecosystems. Plants in water-stressed areas can adapt to arid conditions by gaining access to groundwater [59]. Groundwater accessibility is determined by root systems (depth and density). As a result, deep-rooted forests have greater access to soil water pools than shrubs and grass. When there is a drought, forests can continue to take up the available water from deep soil [60], making them resistant to water imbalance. Meanwhile, this may cause a prolonged drought period for forests, experiencing a longer drought legacy than shrubs and grass (ranging between 1 and 4 years). This can lead to incomplete and delayed forest growth recovery and/or increased canopy disturbance, which is especially noticeable in regions with scarce water resources [61]. Therefore, drought occurring in forests in M-semiarid regions should be paid attention to and improved, such as by strengthening monitoring and irrigation or developing regional ecological environment restoration strategies to reduce the damage of the effects of drought on the forest ecosystem. However, longer drought legacy effects within forests do not necessarily result in weaker drought resilience. Instead, a stronger drought resilience in forests [62] was observed compared to shrubs and grass at timescales ranging from 1 to 3 years.
In contrast, shallow-rooted plants such as shrubs and grass typically take up water from the top-middle soil, where soil moisture quickly responds to precipitation pulses [63]. When drought occurs, the top-middle soil is quickly depleted. Grass organs store insufficient water and carbon to sustain growth during and after a severe drought, making them vulnerable to drought variations [64]. In the meantime, they would recover quickly from drought, and a drought legacy effect is only detectable in semi-arid regions with a maximum time lag of one year [64,65]. This varying pattern of water use can also account for the documented inverse relationship between plant size and severe drought-induced damage [66]. The different ecohydrological responses of various vegetation types to drought should be considered by governments in developing drought mitigation strategies.
Additionally, NEC has the biggest area of freshwater marshes in China [67,68] and plays an important role in conserving biodiversity and regulating the regional climate in northeast Asia [69,70]. Short-term water shortages may encourage the development of wetland vegetation and lead to the weakening of wetland ecological functions [71]. Long-term drought in wetland ecosystems may degrade them into other ecosystems such as meadows and barelands, completely losing their original ecosystem functions [72]. The relatively low drought trend and area of cropland may benefit from the regulation of human activities.
The C-humid (dominated by forests) regions also exhibited a strong trend of drought in both annual and other seasons, except spring. However, unlike M-semiarid areas, short-term precipitation deficits in C-humid areas may result in higher temperatures, thus promoting vegetation development [73,74], as temperature and solar radiation are limiting factors in humid areas or at high latitudes and altitudes [71]. The government should therefore concentrate on the annual drought trends in C-humid regions. The ratio of drought areas in different climatic regions is expected to reflect the trend of drought severity with warming or drying. Our findings indicate that the warmer the temperature zone, the more sensitive it is to drought under the same hydrological conditions, and the higher the ratio of the drought area. The drier the land, the higher the ratio of the drought area under the same temperature zone (Figure 4a), and the stronger the drought trend in spring and summer ( Figure 11). This pattern provides a basis for the development of spatially differentiated drought mitigation strategies. We suggest enhancing drought monitoring in warmer and drier lands, improving and using forecasts for multiple drought types at different time scales, and enhancing regional ecological restoration to improve regional climates.

Conclusions
Drought monitoring at multiple spatial and temporal scales is helpful for understanding the inherent mechanism of hydrological extremes and taking corresponding measures to mitigate drought impacts. In this study, we used trend analysis, the Hurst exponent, and mutation point detection methods to determine the spatiotemporal variation in drought status and trends in NEC from 1990 until 2018 based on the Standardized Precipitation Evapotranspiration Index (SPEI). We found that NEC has exhibited a general dry tendency over the past three decades. The percentage of annual drought-affected areas increased by 0.7%/a and 86.3% of the area exhibited a drying trend. Future dry-wet trends are inconsistent with the current situation. The areas with wetting and drying trends would cover 86.3% and 16.7% of NEC, respectively. From the timeline, the drying trend peaked in 2001, and then exhibited a mitigation tendency before drying again after 2013. The implementation of ecological restoration projects is primarily responsible for drought mitigation. Moreover, we identified the moment and region in NEC where the trend of SPEI12 turned. Western Liaohe Plain and eastern Lesser Khingan Mountain experienced a shift from a drought to a humid trend between 1999 and 2000, and western Sanjiang Plain experienced the same shift in 2004. We proved that there is a time lag between improved vegetation and the trend shift.
The multi-spatiotemporal meteorological perspective in this study raised awareness of different climate regions and land types facing drought risk. Our findings highlight that drought in shrub and grass and sparsely vegetated land (mainly in M-semiarid regions) occurring in summer, particularly in May and June, should be paid attention to and improved, such as by strengthening monitoring, irrigation, artificial rainfall, and developing regional ecological environment restoration strategies. Moreover, the government should concentrate on the annual drought trends in C-humid areas and drought occurring in forests in M-semiarid regions to reduce the damage of the effects of drought legacy to the forest ecosystem. In addition, we suggest enhancing drought monitoring in warmer and drier lands and enhancing regional ecological restoration to improve regional climates.