This study investigated the seasonal variations of the normalized difference vegetation index (NDVI) and its relationships with climatic variables and topography in a small-scale (20 km×20 km) area (i.e., Tsogt-Ovoo village) within the desert steppe zone of Mongolia using in-situ observed climate data and satellite remote sensing data. We found that the topography is very important for vegetation growth in the desert steppe although the summer precipitation is the constraining factor. The unexpectedly high NDVI (upto 0.56), as well as the high aboveground biomass, in the valley bottom was primarily resulted from the topography-modulated redistribution of overland flow after relatively heavy precipitation events during the growing season. This makes the valley bottoms in desert steppes not only reliable feeding resources for livestock but also heavens for wild lives. But, the detected large standard deviation of annual maximum NDVI (NDVImax) from 2000 to 2013 in the valley bottom in response to rather variable precipitation implies that the valley bottoms under desert steppe climates are more vulnerable to climatic change.
The vegetation of an ecosystem not only provides essential ecological, economic, and social products and services but also plays an important role in land-air coupling processes (Ma et al., 2010; Ignatavič ius et al., 2013). Consequently, the coupling relationships between climate variables and vegetation parameters have been extensively investigated (Miyazaki et al., 2004; Li et al., 2007; Munkhtsetseg et al., 2007; Yamanaka et al., 2007). The relevant investigations have recently been further advanced and expanded through using satellite remote sensing data (Yu et al., 2003; Zhang et al., 2003; Pettorelli et al., 2005; Zhu et al., 2013).Among the relevant investigations, land cover analysis using satellite remote sensing data is a well-developed technique for calculating land cover fractions within a pixel of the remotely-sensed data (Ridd, 1995; Yue et al., 2006; Weng and Lu, 2009; Mayes et al., 2015). Calculating land cover fractions has been widely adopted to explore the vegetation responses to climate variations in hyper-arid and semi-arid areas, especially to precipitation variations (Fensholt et al., 2012; Yin et al., 2012; Zhao et al., 2015). The normalized difference vegetation index (NDVI) has been widely used for understanding the dynamic relationships between vegetation parametersand climatevariables(Guo et al., 2014; Jarlan et al., 2014).
The instrumental records show that the temperature has been steadily rising since 1940 in Mongolia and the precipitation has no rising trend at all (Dagvadorj et al., 2014). This simply implies that the effective soil moisture might have been declining due to rising temperature via evaporation enhancement at large spatial scales. However, the effective soil moisture and the associated vegetation growth at small spatial scales might have been somewhat independent from the large-scale patterns. In fact, it was reported that local climate and local topography are often very important in modulating the effective soil moisture and the associated vegetation growth (Wang et al., 2001, 2003; Ni, 2003; Zhou et al., 2003; Miyazaki et al., 2004; Zhang et al., 2005; Cheng et al., 2011; Li and Guo, 2012; Vandandorj et al., 2015).
The Mongolian desert steppeis presently facing severe environmental challenges including desertification and overgrazing (Chuluun and Ojima, 2002; Stumpp et al., 2005; Wang et al., 2006; Lu et al., 2009; Zhang et al., 2009). The desert steppe occupies about 23% of the total territory in Mongolia (Iwasaki, 2006; Johnson et al., 2006; Vandandorj et al., 2015). The vegetation productivity level is rather low in the desert steppe, but the economic significance and ecological importance are not negligible (Saizen et al., 2010). Existing studies have explored the small-scale relationships between the vegetation productivity levels and the controlling or modulating factors in other arid and semi-arid areas. For example, Li and Guo (2014) reported that vegetation growth in semi-arid areas is controlled by such factors as soil moisture, soil temperature, soil types and others. Their study (Li and Guo, 2014) stressed that topography can be one of the most important factors in modulating vegetation growth in semi-arid areas. Our study focused on a small-scale (20 km× 20 km) area (i.e., Tsogt-Ovoo village) within the desert steppe zone of Mongolia and attempted to explore the topographic influence on vegetation growth using in-situ observed climate data and satellite remote sensing data.
The observation site is situatedwithin Tsogt-Ovoo soum (44° 23′ 04′ ′ N, 105° 16′ 60′ ′ E; 1237 m a.s.l.), a village in Umnugovi aimag (a province) in the northern part of the Gobi Desert, Mongolia. The valley, where Tsogt-Ovoo soum (the village) is situated, runs approximately from northwest to southeast with a width of about 10 km (Fig. 1). The observationsite is on the northern slope of the valley. The natural vegetationaround the site was categorized as desert steppe (Yunatov, 1976) and the vegetation coverage ranged from 5% to 15% in July 2011 (Ishizuka et al., 2012). The plant communities were dominated by Anabasis brevifolia, Reaumuria soongarica, Salsola passerina, Nitraria sibirica, and Kalidium cuspidatum (Ishizuka et al., 2012).
Long-term climatic data were collected from the Tsogt-Ovoometeorological station, a surface synoptic observatory operated by the National Agency for Meteorology and Environmental Monitoring of Mongolia. The station is located approximately 5.3 km north of the center of the observation site. The precipitation was measured twice a day using a glass rain gauge and the snow depth was measured at three positions inside the meteorological station once every morning during the snow-covering season.Air temperature and humidity were measured using HMP-1550-10 (Vaisala Corp., Helsinki, Finland) at 1.98 m above the ground. The soil temperature and soil moisture were measured at the depths of 1.0, 2.5, 5.0, 10.0, 15.0, 20.0, 30.0, and 50.0 cm using platinum sensors (Pt100, Climatec Inc., Toshima, Tokyo, Japan) and amplitude domain reflectometry soil moisture sensors (ThetaProbe ML2x, Delta-T Devices Ltd., Burwell, Cambridge, UK), respectively. A data logger (CR1000-XT, Campbell Scientific Inc., North Logan, UT, USA) was used to record the data.
The Terra moderate resolution imaging spectroradiometer (Terra-MODIS) NDVI, which has 16-day composite data for the entire country of Mongolia between 2000 and 2013, was used in this study and the Terra-MODIS image has a fine resolution (250 m). The aboveground biomass in 2012 and 2013 were in-situmeasured by the Information and Research Institute of Meteorology, Hydrology and Environment, Mongolia. The measurements were conducted in a 50 m× 50 m fenced pasture on the 4th, 14th and 24th of each month during the growing season (June to September).
To investigate the topographic effects on vegetation coverage, we traced the variations in the vegetation coverageinthe valley bottom and in the uplandsfrom 2000 to 2013 using standard deviation (SD). To identify the topographic positions with the largest vegetation variability, wecalculated the distribution of the annual maximum NDVI’ s SD between 2000 and 2013 using Equation 1.

Where N is the number of years; xiis the maximum NDVI in the ith year; and µ is the mean annual maximum NDVI between 2000 and 2013.
The trendsin NDVI variations weredetected via the nonparametric Mann-Kendall test and the changing rates were evaluated based on Sen’ s slope. The Mann-Kendall test has been used to detect the trend of normally- or abnormally-distributed time series in meteorological variables (Sayemuzzaman and Jha, 2014) and also in ecological variable. The Mann-Kendall test statistic (S) is computed as,

Where n is the data length; NDVIj and NDVIk are the annual NDVI in year j and k (j> k), respectively.
The variance of S is defined as,

The test statistic Z is calculated as,

In this study, the statistical trend of meanNDVI was evaluated using the Zvalue. A negative Z value indicates a downward trend, a positive Z value shows an upward trend, and a zero Z value means neither detectable downward trend nor detectable upward trend. The slope of the trend (i.e., amount of change per year) was computed using Sen’ s nonparametric method. Time series data f(t)with a continuous monotonic increasing or decreasing trend can be expressed as,

Where Qt is a slope and B is a constant. Qt is estimated by Equation 7.

Where j and k are the year number. Sen’ s estimator of slope Qtis the median of these N values of Qt, calculated as follows:

Where N=n(n-1)/2 is the number of calculated slopes; and n is the length of data (14 in this study).
The annual precipitation at the meteorological station was 116.0 mm in 2012 (wet year) and 60.1 mm in 2013 (dry year). As shown in Figure 2a, 92% of the total precipitation (106.6 mm) fell during the rainy season from May to Septemberin 2012 and 89% of the total precipitation (53.4 mm) fell during the rainy season in 2013. The maximum monthly precipitation in 2012 was observed in July (49.4 mm), whereas the maximum in 2013 was observed in June (22.0 mm). Thetotal precipitation of April and May were small in both 2012 (3.4 mm) and 2013 (1.3 mm).
The study area was normally covered by snow from preceding mid-November to succeeding early March. The largest snow depth occurred in January and the number of snow-covered days was usually less than 50 days (Jambaajamts, 1989). A large amount of snow was observed in the 2012/2013 winter from November of 2012 to February of 2013 (Fig. 2b) and a total of 57 snow-covered days was recorded from 20 December of 2012 to 18 February of 2013. Because the soil temperature was consecutively below 0° C from proceeding early November to succeeding early March (Fig. 2b), the observed decrease in snow depth was speculated to be caused by sublimation rather than melting.
The soil moisture during the late July to September in 2012 was higher than that during the same period in 2013 at all the observed depths (Fig. 2b) and the variations in soil moisture seemed to have closely followed the variations in the precipitation (Fig. 2a). Because both the air and soil temperatures were below 0° C from preceding early November to succeeding early March, the soil moistures immediately after the freezing period (i.e., succeeding late March and early April)was similar to the soil moisture immediately before the freezing period (i.e., preceding October). This phenomenon is known as soil moisture memory (Seneviratne and Koster, 2012). Because of a higher amount of precipitation in the summer of 2012, the soil moisture in the succeeding spring (i.e., spring of 2013) was higher than that in the preceding spring (i.e., spring of 2012). It should be particularly noted that the winter snowfall made no contribution to the soil moisture at all. However, after the soil temperature exceeded 0° C in late March orearly April, snowfall, together with precipitation, can effectively increase the soil moisture. Overall, the variations in soil moisture at shallower depths (i.e., 1, 5, and 10 cm)sensitively responded to the rainfall (see Fig. 2a) and the soil moisture at deeper depth (i.e., 20 cm) was also respondent to the rainfall but with a lower sensitivity.
Vegetation was sparse around the observation site with a vegetation coverage less than 15% (Ishizuka et al., 2012). Because of the variations in topography and the associated variations in meso-scale hydrological conditions in the valley, the spatial heterogeneity of vegetation coveragecould be rather large. Here in this study, we used the NDVI to evaluatethe spatial heterogeneity of vegetation coverage at two different spatial scales: (1) a 20 km× 20 km area centered in the observation site; and (2) a 250 m× 250 m area corresponding to one pixel of the Terra-MODIS image including the observation site. As shown in Figure 3a, the average spatial NDVI in the 20 km× 20 km area had drasticallyincreased from late July and reached the maximum value of 0.17 in the early August of 2012. But, nosuch an increase was observed in 2013. The same patterns were observed for the aboveground biomass that was almost exactly corresponding with the variations in NDVI both in 2012 and 2013. However, the NDVI at the 250 m× 250 m area did not exhibit any drastic increasesboth in 2012 and 2013 (Fig. 3b). Keep in mind that the 250m× 250 m area corresponding to one pixel of the Terra-MODIS image included the observation site that was situated on the northern slope of the valley (i.e., not in the valley bottom).
Phenologically, vegetation starts to grow after the daily mean air temperature exceeds 5° C, normally around late March to early April. But, the NDVIdid not increase until late June, implying that the amount of soil moisture was insufficient for vegetation growth until late June. The best growing time was July-September when the soil moisture was sufficiently recharged by precipitation. For example, the precipitation were 49.4 mmin July and 31.6 mm in August of 2012 and the aboveground biomass drastically increased in late July and reached the maximum value of 210 kg/hm2 in early September (upper panel of Fig. 4). However, both the NDVI and the aboveground biomassdid not increase during the best growing time in 2013 (i.e., July-September) simply due to the fact the precipitation was not sufficient to recharge the soil moisture both in July (17.5 mm) and August (9.5 mm) (lower panel of Fig. 4).
Again, the average spatial NDVI in the 20 km× 20 km area had drasticallyincreased from late July and reached the maximum value of 0.17 in the early August of 2012. But, nosuch an increase was observed in 2013. It is quite notable that the NDVI exceeded 0.33 in the valley bottom from 28 July to 28 August of 2012 and the NDVI (> 0.33) can be considered to be unexpectedly high for a desert steppe. The unexpectedly high NDVI was primarily resulted from the topography-modulated redistribution of overland flow after relatively heavy precipitation events. Specifically, the slope overland flow can be readily concentrated into the valleybottom, forming a temporal pond (Fig. 5a) or a seasonal wetland (Fig. 5b). It should be stressed here that although the NDVI in July and August of 2012 (> 0.33) was unexpectedly high in the valley bottom, the NDVI of the same period was rather low in the uplands (see Fig. 3b), suggesting that the heterogeneity of the NDVI is rather large in the desert steppe and that the large heterogeneity was resulted from the variations in topography and the associated variations in meso-scale hydrological conditions.
The aforementioned differences of NDVI (see Figs. 3 and 4), between non-peak growing season (i.e., from 23 Aprilto 8 June) and peak growing season (i.e., from 28 Julyto 28 August), between the uplands and the valley bottom, and between 2012 and 2013, can be better expressed by plotting the pixel-based NDVI values against the pixel-corresponding altitudes. Figure 6showedthat the altitude did not play any detectable role in the NDVI distribution simply because the relief or altitude difference is not significantly large. In other words, neither orographic precipitation nor environmental lapse rate of temperature had a detectable imprint on the NDVI under the desert steppe climate. Instead, contrast of NDVI between the uplands and the valley bottom was rather striking only during the peak growing season (i.e., from 28July to 28 August) of 2012. Again, the unexpectedly high NDVI in the valley bottom within the context of a desert steppe climate was primarily resulted from the topography-modulated redistribution of overland flow after relatively heavy precipitation events during the peak growing season of 2012. Once more, the unexpectedly low NDVI in the valley bottom even during the peak growing season of 2013 was simply due to the fact the growing-season precipitation of 2013 was not sufficient to recharge the soil moisture.
To investigate the topographic effects on vegetation coverage, we compared the variations in the vegetation coverageinthe valley bottom and in the uplandsfrom 2000 to 2013 using standard deviation (SD). The SD of a14-year average NDVI maximum (NDVImax) (from 2000 to 2013) unequivocally demonstrated that the SD of NDVImaxwas by far larger in the valley bottom than in the uplands (Fig. 7). Specifically, the SD ranged from 0.05 to 0.18 in the valley bottom, whereas the SD was less than 0.05 in the uplands. This means that the vegetation coverage in the valley bottom has a much greater sensitivity to climate change than thatin the uplands and that the elevated sensitivity in the valley bottom was once again resulted from the topography-modulated redistribution of overland flow after relatively heavy precipitation events during the peak growing season. This also implies that the valley bottoms within desert steppe climates are more vulnerable to climatic change although they (i.e., valley bottoms) are more reliablefeeding resources for livestock.
The vegetation coverage in the desert steppe zone of Mongolia hasmuch higher variation in comparison with other vegetation zones in the country because of the large-amplitude variations of precipitation (Vandandorj et al., 2015). Our research demonstrated that the topography is very important for vegetation growth at small spatial scales. For example, the unexpectedly high NDVI and associated high aboveground biomass wereprimarily resulted from the topography-modulated redistribution of overland flow after relatively heavy precipitation events. Moreover, the overland flow from the surrounding high terrains was able to carry the eroded fine soils and organic materials (dead vegetation and livestock dung, etc.) down to the valley bottom, effectively increasing the soil fertility for vegetation growth and thus making the valley bottom a primary feeding source of pastoral livestock in Mongolia (Miller and Sheehy, 2008).
As aforementioned, the trendsin NDVI variations weredetected via the nonparametric Mann-Kendall test and the changing rates were evaluated based on Sen’ s slope.Figure 8presents the mean NDVI trends and the mean changing rates during the four observing periods from 2000 to 2013: (i) 23 April to 7 May; (ii) 25 May to 8 June; (iii) 28 July to 12 August; and (iv) 13 August to 28 August. The overall mean NDVI (2000-2013) in the uplands exhibited an increasing trend (i.e., Sen’ s slope> 0) from 23 Aprilto 12August and a decreasing trend from 12 August onward (i.e., Sen’ s slope< 0). However, the overall mean NDVI (2000-2013) in the valley bottom increased only during the peak precipitation period from 28 July to 12 August. Our interpretation for the longer-term difference of NDVI between the valley bottom and the uplands is as follow. For the uplands, the NDVI increasing trend was a sensitive response to the precipitation increasing trend from 23 April to 12 August and the decreasing trend after 12 August was most likely caused by persistently-high temperature-resulted soil moisture deficit. For the valley bottom, the drastic jump in NDVI from 28 July to 12 August was most likely an exaggerated response to heavy precipitation of the same time period. The decrease of NDVI in the valley bottom after August was also probably caused by persistently-high temperature-resulted soil moisture deficit.
This study found that the precipitation and the associated soil moisture in July and August are the most important factors affecting vegetation growth in the desert steppe of Mongolia. The winter snowfall made no contribution to the soil moisture simply due to the fact that the winter snow was mostly exhausted by sublimation. However, after the soil temperature exceeded 0° C in late March or early April, snowfall, together with precipitation, can effectively increase the soil moisture.
We have also explored the effects of the topography on vegetation growth and the associated aboveground biomass using in-situ observed climate data and Terra-MODIS NDVI data and several conclusions can be drawn. First, the topographic features at small spatial scales are very important for vegetation growth in the desert steppe although the summer precipitation is the constraining factor. Second, the unexpectedly high NDVI (up to 0.56), as well as the aboveground biomass, in the valley bottom under the desert steppe climate isprimarily resulted from the topography-modulated redistribution of overland flow after relatively heavy precipitation events during the growing season. Third, rather large standard deviation (SD) of the annual maximum NDVI (NDVImax) in the valley bottom in response to rather variable precipitation implies that the valley bottoms within desert steppe climates are more vulnerable to climatic change although they (i.e., valley bottoms) are more reliablefeeding resources for livestock.
This research was financially supported by the Japan Society for the Promotion of Science RONPAKU Program (MECS-11319), a Budget Request of Tottori University, and sponsored by special coordination funds from the Ministry of Education, Culture, Sports, Science and Technology of the Japan.
The authors have declared that no competing interests exist.
| 1 |
|
| 2 |
|
| 3 |
|
| 4 |
|
| 5 |
|
| 6 |
|
| 7 |
|
| 8 |
|
| 9 |
|
| 10 |
|
| 11 |
|
| 12 |
|
| 13 |
|
| 14 |
|
| 15 |
|
| 16 |
|
| 17 |
|
| 18 |
|
| 19 |
|
| 20 |
|
| 21 |
|
| 22 |
|
| 23 |
|
| 24 |
|
| 25 |
|
| 26 |
|
| 27 |
|
| 28 |
|
| 29 |
|
| 30 |
|
| 31 |
|
| 32 |
|
| 33 |
|
| 34 |
|
| 35 |
|
| 36 |
|
| 37 |
|
| 38 |
|
| 39 |
|
| 40 |
|
| 41 |
|
| 42 |
|
| 43 |
|

