Extreme Climate Event and Its Impact on Landscape Resilience in Gobi Region of Mongolia

: The dzud , a speciﬁc type of climate disaster in Mongolia, is responsible for serious environmental and economic damage. It is characterized by heavy snowfall and severe winter conditions, causing mass livestock deaths that occur through the following spring. These events substantially limit socioeconomic development in Mongolia. In this research, we conducted an analysis of several dzud events (2000, 2001, 2002, and 2010) to understand the spatial and temporal variability of vegetation conditions in the Gobi region of Mongolia. The present paper also establishes how these extreme climatic events a ﬀ ect vegetation cover and local grazing conditions using the seasonal aridity index ( a AI Z ), time-series Moderate Resolution Imaging Spectroradiometer (MODIS) Normalized Di ﬀ erence Vegetation Index (NDVI), and livestock data. We also correlated a AI Z , NDVI, and seasonal precipitation in the varied ecosystems of the study area. The results illustrate that under certain dzud conditions, rapid regeneration of vegetation can occur. A thick snow layer acting as a water reservoir combined with high livestock losses can lead to an increase of the maximum August NDVI. The Gobi steppe areas showed the highest degree of vulnerability to climate, with a drastic decline of grassland in humid areas. Another result is that snowy winters can cause a 10 to 20-day early peak in NDVI and a following increase in vegetation growth. During a drought year with dry winter conditions, the vegetation growth phase begins later due to water deﬁciency, which leads to weaker vegetation growth. Livestock loss and the reduction of grazing pressure play a crucial role in vegetation recovery after extreme climatic events in Mongolia.


Introduction
Featuring a mix of monsoon, and harsh Siberian climates, the south Gobi of Mongolia is a unique desert steppe and nomadic living area, and one that is sensitive to global climate changes. Over the past 20 years, meteorological disasters in the Mongolian plateau have become increasingly severe, and its dynamic succession have become increasingly investigated partly through large-scale monitoring. The interaction between drought-dzud and vegetation has been an essential concern for pastureland use in Mongolia. Originally, dzud was the Mongolian term for a natural, winter-related climate disaster.
Dzuds are now defined, biophysically, as anomalous climatic or land surface conditions (i.e., snow and ice cover) that lead to reduced accessibility or availability of pastures [1]. It often causes very high Table 1. Types of dzuds in Mongolia and their local nomenclature [7,8].

Type of Dzud Weather Condition Effects
Tsagaan (white) Average thickness of snow layer on pasture land exceeds 21 cm in high mountains and forest regions, 16 cm in steppes, and 10 cm in the Gobi region. Snow density reaches 200 kg/m 3 or greater in any region.
Prevents access to grass. The most common form of dzud and the most disastrous when it affects large areas.

Khar (black)
No snow during winter and the monthly or ten-day average temperature is 5.0 • C lower.
Causes water supply shortages, often exacerbated by lack of winter grass.
Tumur (iron or ice) Snow cover melts and freezes due to rapid changes in temperature creating an ice cover that prevents livestock from grazing. Snow density reaches 0.30 g/cm 3 or greater.
Prevents access to grass.

Khuiten (cold)
Air temperature is 5 to 10 • C lower than the monthly mean temperature for several consecutive days.
Extreme cold and strong freezing winds prevent animals from grazing. Animals spend most of their energy maintaining their body heat.
Turen (hoof) Extremely dry weather. Causes complete depletion of grass due to drought and/or trampling and heavy grazing. Khavsarsan (combined) Two or more of the above occurring simultaneously.
The United Nations Convention to Combat Desertification (UNCCD) defines desertification as "land degradation in arid, semi-arid and dry sub-humid areas resulting from various factors, including climatic variations and human activities" [14]. Desertification, therefore, denotes a form of land degradation specific to drylands. The UNCCD defines drought as "the naturally occurring phenomenon that exists when precipitation has been significantly below normally recorded levels, causing serious hydrological imbalances that adversely affect land production systems" [14]. Droughts increase the frequency and severity of dzud [1,15]. A drought mechanisms analysis revealed that severe water scarcity results from a combination of insufficient precipitation, high evaporation, and over-exploitation of water resources [16]. It is important to note that reduced winter precipitation promotes drier summer conditions by reducing the soil water available for evapotranspiration in summer [17].
The livestock sector of Mongolia has struggled with the harsh climate, overgrazing, and an imbalance in livestock species. These factors have led to significant damage to regional vegetation [18]. The winter of 2010 was the "most severe winter in nearly five decades" [19], in which more than 10 million livestock died across the country. This represented a national disaster for a country where animal husbandry is the traditionally dominant form of livelihood, providing 16% of the Gross Domestic Product (GDP) for the country and employment to 366,000 herders, as well as indirect economic support to about one million people or about a third of the 3.1 million inhabitants [20]. The impact on many herders was particularly severe during the three consecutive dzud winters of 1999 to 2002, when 30% of the livestock in the nation perished [21,22]. More recently, the dzuds of 2009 to 2010 were even harsher. [21] analyzed the changing climate over the past 70 years (Appendix A), stating that the frequency and range of dzuds and summer droughts are expected to increase, particularly in central and northern Eurasia [23,24]. The effects can vary according to the type of dzud ( Table 1).
Analysis of extreme climate events has become an urgent task to understand dzud mechanisms in Mongolia. By studying the mechanisms of dzuds, we will be able to prepare herders and their livestock for dzud effects. Limited research has examined the response of vegetation conditions to dzuds and summer droughts for different ecosystems on the Mongolian plateau using long-term satellite observations [25,26]. One of the most cost-effective remote sensing approaches to detect vegetation changes across large scales is the Normalized Difference Vegetation Index (NDVI) [27]. NDVI features a nonlinear combination of the red and near-infrared (NIR) spectral irradiance, exhibiting a relationship with vegetation and green biomass. The NDVI index has been widely used to monitor global vegetation phenology from space across large scales [28,29]. This study focuses on the impacts of dzud-related factors (e.g., drought and snowfall), as well provides a more in-depth analysis of its spatial and temporal distribution and losses. The goal of this paper is to describe the effects of dzud on the vegetation growth cycle in the subsequent year. To comprehend the mechanical concept of dzud after droughts, and how vegetation responds to various dzud situations, we used MODIS-NDVI data (2000)(2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013) to estimate the spatiotemporal changes in NDVI. In this paper, we (1) determine the spatiotemporal variations in NDVI after dzud events and the length of the growing season, and (2) assess the impact of drought (using the climatic variables Aridity Index ( a AI z ), precipitation, and temperature) on vegetation throughout the dzud period.

Study Area
The study area comprises the Bayankhongor and Gobi-Altai provinces (aimags) in southwestern Mongolia (Figure 1), with a total size of 258,200 km 2 [30]. This is part of the central Asian desert belt [31] and has an arid ecosystem [32,33]. It has an elevation that ranges from the sea level to 4000 m.a.s.l. and the average temperature is 17 • C in July and −24 • C in January. The area has become drier and hotter in recent years, due to the impact of climate change [34]. The Siberian high-pressure cell in this area causes cold and dry conditions in winter [35]. The annual average rainfall varies from 200 to 300 mm, with an annual average wind speed of 3.1 m/s. The climate datasets and maps of vegetation zones were provided by the Institute of Meteorology and Hydrology of Mongolia. Six meteorological stations were examined for each province ( Table 2). The dominant plant species are Stipa grandis, Leymus chinensis, and multiple species of Artemisia spp. and Festuca spp. [36]. Annual germination occurs between April and July, depending on the antecedent soil moisture and rainfall [33,37,38].

Methods
In this study, we combined Moderate Resolution Imaging Spectroradiometer (MODIS) NDVI satellite data products, climate data from the Mongolian meteorological stations, spatial snow-cover data from global climate data sets, and livestock data from the National Statistics Office of Mongolia [39]. We examined extreme climate events NDVI trends, and the relationship between NDVI and aridity index ( a AI z ) to characterize vegetation growth cycles after dzud events. Figure 2 presents a framework of the study that comprises pre-processing and meteorological data processing to compute a AI z , remote sensing data processing of NDVI, precipitation, temperature, snowfall, and livestock mortality. We examined vegetation variation in the Gobi region by using MODIS NDVI (MOD13Q1) data [40], which provides global data records with a spatial resolution of 250 m and a time resolution of

MODIS NDVI Data (Vegetation Dataset)
We examined vegetation variation in the Gobi region by using MODIS NDVI (MOD13Q1) data [40], which provides global data records with a spatial resolution of 250 m and a time resolution of 16 days from the year 2000 to present. We used data to 2013 to extend over the most severe dzud years. The MOD13Q1 products are derived from atmospherically corrected bi-directional surface reflectance function (BDRF) that have been masked for clouds, shadows, water and aerosols [40]. Comprehensive information and characteristics of MODIS NDVI can be obtained from [41] in TIFF format [42]. In this study, we only processed NDVI data during dzud related years, specifically during and after dzud events to estimate the mean monthly August NDVI value. In addition, we ensured that the areas were at least 500 × 500 m in order to fit the spatial resolution of the MODIS sensor (250 × 250 m). NDVI is defined by the following equation: where red is the visible light of the red wavelength (from 400 to 700 nm) and NIR is the intensity of the near-infrared wavelength (from 700 to 1100 nm). The formula uses the specific reflective behavior of plant surfaces. NDVI values are normalized to a range from −1 to 1, where gross values over 0.1 represent vegetation. Higher value indicates more vegetation, while negative values indicate water, ice, or snow surfaces [43]. The MODIS NDVI data were then transformed into the same geographic coordinate system (UTM48N) to ensure coverage equality. We used the "buffer" tool in ArcGIS (with a radius of 10 km) around the climate stations, afterward, computed the average MODIS NDVI for each area.

Meteorological Data (Climate Dataset)
The National Agency of Meteorology and Environment Monitoring of Mongolia (NAMHEM) provided climate datasets (i.e., precipitation and temperature) from twelve different stations ( Figure 1). The station data were measured and documented for all study regions as average monthly temperature and total monthly rainfall. There were some limitations in data availability at (Khukhmorit, Bugat, and Tooroi) stations. The data from 1990 to 2013 were selected for aridity index analysis, while the data from 1985 to 2014 were selected for climate trend analysis in Bayankhongor and from 1990 to 2014 for the climate trend analysis in Gobi-Altai. The maximum and minimum temperatures were arithmetically derived from this data set. Rainfall and temperature records were averaged by season: summer (June, July, August), autumn (September, October, November), winter (December, January, February), and spring (March, April, May). We divided the time range (30 years) into three groups; 1985-2014, 1985-1999, 2000-2014, and trends were detected. The snowfall map was derived from the Climatic Research Unit (CRU) TS 3.22 global climate records for the period 1999 to 2010. The dataset was provided as an accumulation of the precipitation (i.e., snowfall) amount from November to March. The gridded CRU TS 3.22 data are monthly variations in weather from 1901 to 2013 with high resolution (0.5 × 0.5 degree) grids, created by the Climatic Research Unit (CRU) at the University of East Anglia [44]. The gridded data, along with the monthly station observations, are freely available at the British Atmospheric data center [45].

Socio-Economic Data (Statistical Dataset)
The livestock data used is from 2000 to 2013, and was provided by the National Administrative Department of Statistics [46]. The Ministry of Livestock Husbandry composed the mortality data of the livestock each December, and the National Statistical Office (NSO) of Mongolia (1990-2015) collected the total livestock numbers for all provinces. We combined the statistical results to check the spatial and temporal variability of dzud events and the relationship between climate variability and livestock loss. To examine livestock mortality, livestock data were first converted to equivalent 'sheep forage units' (SFU) to standardize feed requirements, as each type of livestock requires different amounts of feed [47]. The conversion rates to SFU were: 1 camel = 5 SFU, 1 horse = 7 SFU, 1 cow = 6 SFU, 1 sheep = 1 SFU and 1 goat = 0.9 SFU. The percent relative mortality for each province was calculated as the ratio of livestock deaths during winter and spring to the number of livestock at the beginning of the year. Consequently, the examined stock loss can be either a factor of the dzud of the prior year (in the case of drought temperatures) or of the current year (in the case of winter temperatures).

The Implementation of the Aridity Index
Drought indices often integrate precipitation, temperature, and other variables, but may emphasize different aspects of drought and must be carefully selected concerning the drought characteristics. For this purpose, we assessed changes in drought characteristics by focusing on the aridity index (AI z ). Due to the local arid climatic conditions, we applied an adapted approach of the aridity index based on the respective drought years [48]. This method was originally developed for the southern regions of Mongolia, to define the drought intensity during the summer months to estimate the relevant connection to pasture yields. In our study, it is intended to link spatially variable climatic circumstances and evaluate specific summer characteristics before dzud events. Hence, the aridity index was applied to each province [49] to develop a zonal adapted seasonal aridity index ( a AI z ). Subsequently, the a AI z was used to show the drought risk and periodic temporal patterns in NDVI. We then arranged the standardized a AI Z values into three groups, represented by values less than zero (dry years), values between zero and five (normal years), and values greater than five to show the wet years.
This method is based on daily temperature and precipitation values to assess the summer weather. The aridity index itself is based on a supply and demand concept of the water balance equation. The annual mean evapotranspiration and runoff rates are regulated by the amount of available energy (demand) and precipitation (supply). Equation (2) defines AI as a function of precipitation (P), the supply, and temperature (H i ), the demand.
where P7 is the cumulative precipitation in July, while H i 6 and H i 7 are the average maximum daily temperatures during June and July, respectively. The basic theory behind this approach lies in the disproportionate significance of the July precipitation on biomass growth and its maximum in mid to late August. Within the southern parts of the study area, rainfall is usually limited to the summer months of July and August, with rainfall in August having a lower impact on the maximum growth stated, due to a response period of 10 to 20 days [50]. Evapotranspiration is highest in June due to the long days, high levels of incoming solar radiation, and high temperatures. Evapotranspiration can also be high in July due to the hot conditions [48]. This aridity index approach demonstrates a stronger correlation to pasture yield in the southern part of the study area (i.e., arid climates) than in the central and northern parts. An expanded AI calculation as a zonal component (AI Z ) is beneficial for regional assessment in the north, where precipitation increases in July and August, and a growth period spans from May to September (P May-Sept ). The denominator was extended by an empirically derived number (10) to compensate for different locations [51,52]. Equation (3) represents the zonal aridity index.
where H i (Tmax June-July ) signifies the sum of the maximum daily temperatures from June and July, and P r is the total precipitation during the respective season (in this case, summer). The seasonal precipitation (P r ) was applied to the zonal aridity index formula because seasonal and annual precipitation is the main factor that controls the spatial distribution of plant activity [53]. Greater negative values indicate increased aridity (i.e., drought). With this approach, we have developed our own regionally adapted index. In our study, the daily temperature data from June and July were accessible to find the maximum monthly temperature. The adjusted aridity index ( a AI Z ) is: where 6m and 7m stand for the average daily temperature in June and July, respectively. This value allows us to estimate the drought intensity at individual study area stations. Lower values of the a AI Z generally represent drier conditions in summer and a lower than maximum NDVI in August. Nevertheless, the prospect of direct comparability of the study area stations is problematic, since they are subject to different mean values and standard deviations. The steppe vegetation area shows a particularly strong correlation in terms of the a AI z , whereas the more arid desert regions show a lower correlation. This suggests that the vegetation that is present is very sparse or only present in small areas, so that no significant vegetation signal can be detected. Thus, the calculated a AI Z values confirm the spatial variability of climatic conditions. Pearson's correlation coefficients were applied to determine the sensitivity of a AI Z to different seasonal conditions as well as the relationship between mid-August NDVI and seasonal precipitation.

Spatio Temporal MODIS NDVI Analysis
We examined the spatial variability of MODIS NDVI data for dzud events in 2009/2010. Figure 3 shows the NDVI spatial patterns in southwestern Mongolia (encompassing twelve different stations). The NDVI vegetation value of dry, shrub, desert regions were typically under 0.1. While NDVI values, associated with different vegetation types, varying environmental conditions, and station locations followed a similar trend. The maps designate that NDVI increased approximately (62.2%) in 2010, mostly in high mountain (steppe) regions in the northern part of Bayankhongor and the north-central part of Gobi-Altai. These results suggest that MODIS NDVI data showed a strong recovery after the 2010 dzud events (Figure 3b). The NDVI data from 2010 exhibited maximal vegetation growth (mid-August). This early growth phase depends on the weather conditions and the quantity of the snow cover. Through the exceptionally wet year of 2010 ( a AI z = 0.91), vegetation reached its maximum interim growth. A dzud may be a determinant for a more substantial vegetation growth phase (outside of extremely arid desert zones). NDVI spatial distribution maps after dzud events also allowed delineation of wet areas in the northwestern and southwestern dry regions of Mongolia.       Positive NDVI trends can only be expected in the presence of relatively high precipitation. In 2002, a lack of precipitation led to substantial NDVI losses compared to the previous years. This pattern was widespread between the different station locations. The northern region vegetation had a stronger negative reaction to extreme droughts, while the southern region vegetation was well adapted to short-term and long-term droughts. The examination of vegetation response indicates that, generally, drought and water supply due to snow accumulation are essential in terms of the impact of dzud on vegetation development.  Positive NDVI trends can only be expected in the presence of relatively high precipitation. In 2002, a lack of precipitation led to substantial NDVI losses compared to the previous years. This pattern was widespread between the different station locations. The northern region vegetation had a stronger negative reaction to extreme droughts, while the southern region vegetation was well adapted to short-term and long-term droughts. The examination of vegetation response indicates that, generally, drought and water supply due to snow accumulation are essential in terms of the impact of dzud on vegetation development. adapted to short-term and long-term droughts. The examination of vegetation response indicates that, generally, drought and water supply due to snow accumulation are essential in terms of the impact of dzud on vegetation development.

The Relationship between MODIS NDVI, Aridity Index (aAIz), and Seasonal Precipitation
In order to check how MODIS NDVI represents the spatial differences in vegetation, we examined the relationship between NDVI and the climate variables (adapted aridity index (aAIz) and seasonal precipitation; Figure 6). The results show that there is a significant relationship in the spatial patterns of NDVI; MODIS NDVI statistically significant (p < 0.001) correlated with aAIz and seasonal precipitation in the steppe regions with R 2 = 0.64 and R 2 = 0.59, respectively (Figure 6a,b). Furthermore, lower correlations were found in arid desert regions with R 2 = 0.51 and R 2 = 0.56, respectively ( Figure  6c,d). The correlation for different vegetation types is higher between NDVI and aAIz than with NDVI and seasonal precipitation.

The Relationship between MODIS NDVI, Aridity Index ( a AI z ), and Seasonal Precipitation
In order to check how MODIS NDVI represents the spatial differences in vegetation, we examined the relationship between NDVI and the climate variables (adapted aridity index (aAIz) and seasonal precipitation; Figure 6). The results show that there is a significant relationship in the spatial patterns of NDVI; MODIS NDVI statistically significant (p < 0.001) correlated with aAIz and seasonal precipitation in the steppe regions with R 2 = 0.64 and R 2 = 0.59, respectively (Figure 6a,b). Furthermore, lower correlations were found in arid desert regions with R 2 = 0.51 and R 2 = 0.56, respectively (Figure 6c,d). The correlation for different vegetation types is higher between NDVI and aAIz than with NDVI and seasonal precipitation. The aAIz performed as expected, highlighting drought years and dzud years. Droughts caused reduced biomass and led to a decrease in the amount of hay/forage [15]. We compared drought occurrences with aAIZ and created threshold criteria for the aAIZ values. From 2000 to 2002, in particular, the effects of drought are apparent in these study regions. The dry years also contain massive dzud disasters, confirming the connection between the summer drought and the intensity of dzud events. The results of developed aAIZ index (with their accompanying stations) is represented from 1990 to 2013 in Table 3. When there were high rainfall totals during the growth period, the mean The a AI z performed as expected, highlighting drought years and dzud years. Droughts caused reduced biomass and led to a decrease in the amount of hay/forage [15]. We compared drought occurrences with a AI Z and created threshold criteria for the a AI Z values. From 2000 to 2002, in particular, the effects of drought are apparent in these study regions. The dry years also contain massive dzud disasters, confirming the connection between the summer drought and the intensity of dzud events. The results of developed a AI Z index (with their accompanying stations) is represented from 1990 to 2013 in Table 3. When there were high rainfall totals during the growth period, the mean NDVI was 0.72 points above average. Therefore, a strong influence of July precipitation on the maximal growth phase in mid/late August should be evident.  (Figure 7a,b). However, NDVI increased in the steppe and mountain regions for most of the investigated area (Figure 7c). The contribution of water supply in winter, a decrease of grazing pressure (due to livestock loss), and recovery in spring-summer rainfall played a minor role in vegetation regeneration in the semi-arid and desert regions.  (Figure 7a,b). However, NDVI increased in the steppe and mountain regions for most of the investigated area (Figure 7c). The contribution of water supply in winter, a decrease of grazing pressure (due to livestock loss), and recovery in spring-summer rainfall played a minor role in vegetation regeneration in the semi-arid and desert regions.

Climate Condition Analysis
We examined the spatial distribution of CRU snow data during dzud events. Figure 8

Climate Condition Analysis
We examined the spatial distribution of CRU snow data during dzud events. Figure 8 shows  We also examined the trends of annual mean precipitation and temperature as well as monthly precipitation and temperature. The mean annual precipitation and temperature differ significantly over the observation period. Figure 9 shows the time series of mean annual precipitation and temperature from 1985 to 2014 in Bayankhongor province and from 1990 to 2014 in Gobi-Altai province. Annual mean precipitation was 202 mm from 1985 to 2014, 226 mm from 1985 to 1999, and 179 mm from 2000 to 2014 in Bayankhongor province (Figure 9a,b). During the combined droughtdzud years, less precipitation was observed, i.e., 124. As for temperature, a general trend towards a warmer climate was recognized in the period from 1985 to 2014. This result is consistent with [31]. We also examined the trends of annual mean precipitation and temperature as well as monthly precipitation and temperature. The mean annual precipitation and temperature differ significantly over the observation period. Figure 9 shows the time series of mean annual precipitation and temperature from 1985 to 2014 in Bayankhongor province and from 1990 to 2014 in Gobi-Altai province. Annual mean precipitation was 202 mm from 1985 to 2014, 226 mm from 1985 to 1999, and 179 mm from 2000 to 2014 in Bayankhongor province (Figure 9a,b). During the combined drought-dzud years, less precipitation was observed, i.e., 124. As for temperature, a general trend towards a warmer climate was recognized in the period from 1985 to 2014. This result is consistent with [31].   [1985][1986][1987][1988][1989][1990][1991][1992][1993][1994][1995][1996][1997][1998][1999]. We determined the mean yearly precipitation during 2000-2014 is between 15 and 20% lower than 1985-1999. Contrary to the expected strong warming tendency in the winter months, the last 10-15 years showed a cooling in winter, which could be related to the Siberian cold high. This circumstance, coupled with the high drought intensity, was a fundamental factor in the occurrence of the consecutive dzuds from 2000 to 2002. The average annual temperature of the last 15 years (2000-2014) is higher than the temperature from 1985 to 1999. We detected the largest increase in annual temperature in the drier regions (southern area). The mean winter temperature from 2000 to 2014 (−3 to 2 °C) is cooler than the mean winter temperature from 1985 to 2014 (−0 to 9 °C). The summer season shows the biggest gain in temperature, while the gain during autumn is quite low. Several studies have discovered an increase in Eurasian snow cover over the past decades [54,55] and increasing trends in anomalous fall melt events [56] through remote sensing.  [1985][1986][1987][1988][1989][1990][1991][1992][1993][1994][1995][1996][1997][1998][1999]. We determined the mean yearly precipitation during 2000-2014 is between 15 and 20% lower than 1985-1999. Contrary to the expected strong warming tendency in the winter months, the last 10-15 years showed a cooling in winter, which could be related to the Siberian cold high. This circumstance, coupled with the high drought intensity, was a fundamental factor in the occurrence of the consecutive dzuds from 2000 to 2002. The average annual temperature of the last 15 years (2000-2014) is higher than the temperature from 1985 to 1999. We detected the largest increase in annual temperature in the drier regions (southern area). The mean winter temperature from 2000 to 2014 (−3 to 2 • C) is cooler than the mean winter temperature from 1985 to 2014 (−0 to 9 • C). The summer season shows the biggest gain in temperature, while the gain during autumn is quite low.
Several studies have discovered an increase in Eurasian snow cover over the past decades [54,55] and increasing trends in anomalous fall melt events [56] through remote sensing.

Livestock Mortality Analysis
When considering the livestock mortality rate during consecutive dzud events, the most severe livestock loss was from 2000 to 2002 and from 2009 to 2010 ( Figure 10). The highest livestock loss rate was nearly 60% in Bogd, Shinejinst, Ekhiingol, Bayankhongor, and Tooroi regions during the 2000/2002 dzud events. Regions further west in Mongolia contained the most mortality hotspots, with 80% of the total mortality [39]. High mortalities indicate poor body conditions of the livestock associated with adverse vegetation conditions due to drought or pasture degradation. During the 2010 dzud, the herds could not access pasture due to deep snow [12]. The increase in pasture intensity reduces vegetation cover and it reinforced the effects of animal hooves [57].

Livestock Mortality Analysis
When considering the livestock mortality rate during consecutive dzud events, the most severe livestock loss was from 2000 to 2002 and from 2009 to 2010 ( Figure 10). The highest livestock loss rate was nearly 60% in Bogd, Shinejinst, Ekhiingol, Bayankhongor, and Tooroi regions during the 2000/2002 dzud events. Regions further west in Mongolia contained the most mortality hotspots, with 80% of the total mortality [39]. High mortalities indicate poor body conditions of the livestock associated with adverse vegetation conditions due to drought or pasture degradation. During the 2010 dzud, the herds could not access pasture due to deep snow [12]. The increase in pasture intensity reduces vegetation cover and it reinforced the effects of animal hooves [57].

Discussion
Our findings are consistent with other studies indicating that summer droughts and cold, snowy winters seem to be the major factors leading to dzud disasters [13,25,58]. As other studies documented that Mongolia is facing severe economic damage and huge livestock losses (Appendix B) (e.g., [59]) dzud might influence the vegetation phenology and livestock mortality. Our study finds that dzud, a combination of drought and harsh winter, contributes to the annual high livestock mortality in southwestern Mongolia, which is consistent with earlier studies [7,60]. The climate projections show that these trends will be aggravated in the medium term and those extreme events, such as drought and dzud, will become more frequent and more intense, with severe environmental, social, and economic impacts [13]. In addition, drought/dry years and lack of meltwater in the spring may cause high grazing pressure due to less vegetation resource for livestock. Moreover, reductions in vegetation cover due to high livestock grazing could negatively impact ecosystem function and increase the vulnerability [61,62]. We also noticed that livestock grazing is difficult after a cold dzud with a high snowfall that does not melt and covers the grassland during the entire winter. Vegetation is susceptible to the length of the rainy season. The years with low NDVI had very dry conditions during the following summer. Most investigated stations with dry conditions (minimum AI of −1) showed a low NDVI. We observed that the highest increase in annual temperature was found in the drier regions in the south of the study area, while stations in the northern parts of the study region exhibited a lower mean annual temperature due to lower temperatures during the winter. Years with drought risk and long dry periods heavily influence the local water budget by causing high evapotranspiration, increased capillary effect, and short-term soil salinization. This suggest that a

Discussion
Our findings are consistent with other studies indicating that summer droughts and cold, snowy winters seem to be the major factors leading to dzud disasters [13,25,58]. As other studies documented that Mongolia is facing severe economic damage and huge livestock losses (Appendix B) (e.g., [59]) dzud might influence the vegetation phenology and livestock mortality. Our study finds that dzud, a combination of drought and harsh winter, contributes to the annual high livestock mortality in southwestern Mongolia, which is consistent with earlier studies [7,60]. The climate projections show that these trends will be aggravated in the medium term and those extreme events, such as drought and dzud, will become more frequent and more intense, with severe environmental, social, and economic impacts [13]. In addition, drought/dry years and lack of meltwater in the spring may cause high grazing pressure due to less vegetation resource for livestock. Moreover, reductions in vegetation cover due to high livestock grazing could negatively impact ecosystem function and increase the vulnerability [61,62]. We also noticed that livestock grazing is difficult after a cold dzud with a high snowfall that does not melt and covers the grassland during the entire winter. Vegetation is susceptible to the length of the rainy season. The years with low NDVI had very dry conditions during the following summer. Most investigated stations with dry conditions (minimum AI of −1) showed a low NDVI. We observed that the highest increase in annual temperature was found in the drier regions in the south of the study area, while stations in the northern parts of the study region exhibited a lower mean annual temperature due to lower temperatures during the winter. Years with drought risk and long dry periods heavily influence the local water budget by causing high evapotranspiration, increased capillary effect, and short-term soil salinization. This suggest that a low NDVI and the low soil moisture conditions can induce drought, and these impacts tend to be a precondition for dzud. Our adapted aridity index showed a considerable regional difference in NDVI from the climatic effects. The NDVI was closely correlated with a AI z , while seasonal precipitation showed a lower correlation with NDVI ( Figure 6). In general, the maximum vegetation growth, irrespective of the annual precipitation conditions, remains stable. When there is a lack of snow cover, there is a lack of water available at the beginning of the growth phase. This is likely reducing the pasture growth due to water stress. For instance, a large amount of livestock died at a later time due to weakness or dehydration, as they were burdened during the winter due to weakened fodder supply. This is limited mainly by winter precipitation and the livestock loss rate, which has emerged as the most important condition of vegetation regeneration. Accumulated snow cover will increase the amount of soil moisture in spring and result in a positive effect on the vegetation cycle. We recognize that our study was limited to investigating recent climate related impacts on vegetation cover. Long-term climate trends were not assessed and may play a significant role, particularly if severe dzud events are becoming more frequent.

Conclusions
This study examined the spatiotemporal pattern in NDVI and the cycle of the growing season after different dzud events detected from MODIS data. Our analysis quantifies the understanding of different dzud mechanisms and their influence on vegetation and livestock grazing. Our results reveal spatiotemporal patterns of vegetation response to dzud and combined summer drought-dzud, and it can help assess the resilience of vegetation in different ecological zones. To clarify the role of dzud climate drivers, we assessed the relationship between the aridity index and summer NDVI (as a proxy for vegetation conditions) and found a significant correlation. The summer conditions represented by aridity index and drought risk have an impact on livestock mortality, which is heightened by dzud. The drought associated dzud years corresponded with lower summer NDVI. Specifically, we assume that during dry winter conditions, the growth phase begins later due to water deficiency and leads to a weaker and slightly later vegetation growth peak. However, the dzuds from 2009 to 2010 coincide with greater livestock losses corresponded and higher summer NDVI. Alleviating the impacts of climatic stresses (drought and dzud) on vegetation will be a crucial challenge in arid and semi-arid regions of Mongolia.
The consequences of changes in the frequency and intensity of environmental disasters have become a considerable issue for regional herders and the well-being of local communities. As a further study, the combined approach of herder's ecological knowledge and remote sensing is an opportunity to explore the speed of vegetation reaction to dzud and droughts.
Author Contributions: O.V. analyzed the data and performed the experiments; she also computed the statistical analysis in R-software. M.K. helped to conceive and design the statistical analysis. T.R. and S.R.F. contributed analysis tools and corrected the final manuscript. All authors together developed and discussed the manuscript and finally wrote the paper. All authors have read and agreed to the published version of the manuscript.  Sources: [59].