Next Article in Journal
Mining and Restoration Monitoring of Rare Earth Element (REE) Exploitation by New Remote Sensing Indicators in Southern Jiangxi, China
Next Article in Special Issue
Intravaginal Devices and GNSS Collars with Satellite Communication to Detect Calving Events in Extensive Beef Production in Northern Australia
Previous Article in Journal
Indirect Measurement of Forest Canopy Temperature by Handheld Thermal Infrared Imager through Upward Observation
Previous Article in Special Issue
A Pilot Study to Estimate Forage Mass from Unmanned Aerial Vehicles in a Semi-Arid Rangeland
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Use of a Satellite-Based Aridity Index to Monitor Decreased Soil Water Content and Grass Growth in Grasslands of North-East Asia

1
Arid Land Research Center, Tottori University, Tottori 680-0001, Japan
2
Graduate School of Engineering, Nagasaki University, Nagasaki 852-8521, Japan
*
Author to whom correspondence should be addressed.
Remote Sens. 2020, 12(21), 3556; https://0-doi-org.brum.beds.ac.uk/10.3390/rs12213556
Submission received: 18 September 2020 / Revised: 27 October 2020 / Accepted: 29 October 2020 / Published: 30 October 2020

Abstract

:
Numerous simulation studies of the effect of global warming on arid regions have indicated that increases in temperature and decreases in precipitation will trigger water shortages, drought, and further aridification. In north-east Asia, especially China and Mongolia, the area of degraded land has increased since 2000. Land use in arid regions is mainly natural grasslands for grazing. Growth in this land use is limited by the precipitation amount and intensity. To develop sustainable management of grasslands, it is essential to examine the relationship between water consumption and the growth patterns of the grasses. This study examined the applicability of a satellite-based aridity index (SbAI) as a way to measure the water consumption and growth of grasslands in China and Mongolia. The effective cumulative reciprocal SbAI was strongly correlated with the cumulative decreased soil water content in the root zone and changes in the normalized difference vegetation index in Shenmu, China. Application of the effective cumulative reciprocal SbAI to grasslands in Mongolia and in north-east Asia revealed a high correlation between the effective cumulative reciprocal SbAI and changes in the normalized difference vegetation index (NDVI). The effective cumulative reciprocal SbAI might be suitable for the detection of water consumption and growth in grasslands from satellite data alone.

1. Introduction

Drylands are highly vulnerable to climate change. In particular, increasing global warming can increase temperatures and decrease precipitation in drylands at high latitudes [1,2,3], and their effects on ecosystems can increase [4,5]. Drought has been increasing in Mongolia since 2000 [6,7]. The main land use in the arid regions of north-east Asia, especially in China and Mongolia, is natural grasslands, and the main source of livelihood is stock grazing on grasslands. However, the growth of natural grasslands depends on unstable rainfall [8,9,10]. The coverage of natural grasslands has a profound effect on the occurrence of dust outbreaks [11,12,13], which bring damage to agriculture, stock grazing, and human and animal health, not only locally but also in neighboring regions such as Japan and Korea [14,15,16]. An early warning and monitoring system based on numerical models, remote sensing, and weather forecasts is urgently needed to guard human well-being and to control natural hazards in those regions. Such models must take into account the growth of vegetation (in this case, grasslands).
Numerical models of vegetation growth are different in climatic features [17]. In humid regions like Japan where precipitation is abundant (1718 mm average from 1971 to 2000: [18]), the quantity of growth depends on the temperature. In arid regions with an adequate sunshine duration and temperature, however, it depends strongly on precipitation and thus soil water content (SWC). Changes in the normalized difference vegetation index (∆NDVI) were closely related to accumulated rainfall and decreased SWC in natural grasslands in China [19]. Thus, observed values of rainfall and SWC as model inputs will affect the accuracy of numerical simulations of plant growth.
However, meteorological data are limited in arid regions, especially rainfall, leading to uncertainties in the simulated values of SWC. With a high resolution and frequency, satellite data offer advantages in monitoring vegetation growth in arid regions [20]. For example, the Moderate Resolution Imaging Spectroradiometer (MODIS) and Copernicus Missions (specifically Sentinel 1 or 2) has provided data since 2000 and 2014, respectively. Because lengthy cloudless periods are common in arid regions, much of the MODIS or Sentinel data are usable for global analyses [21]. However, implementation of remote sensing and imagery techniques to monitor and detect the effects of extreme weather events and climate change on vegetation growth is still underdeveloped [22].
In Japan, where precipitation is abundant, the effective cumulative temperature (ECT) is commonly used to simulate vegetation growth [23,24,25]. ECT is the effective cumulative daily mean temperature above a threshold (e.g., 10 °C in [25]) and is related to plant growth [24]. ECT has been developed into developmental indices, calculated by summing the developmental rate with respect to time [26,27,28].
Because vegetation growth in arid regions depends on precipitation or SWC, we got the idea from ECT and developmental index methods to replace the temperature index with a water index obtained from satellite data to monitor growth. Products using passive or active microwave sensors, such as the Advanced Microwave Scanning Radiometer–EOS (AMSR-E) and AMSR-2, are used to detect SWC [20,29,30,31,32]. However, since their resolution is low (e.g., 50 km in AMSR-2), the accuracy of the soil moisture estimation is low [31,32]. However, the satellite products usable for SWC estimation, from thermal inertia calculated from land surface temperature (LST), have a high resolution and high observation frequency, and give good agreement with observed SWC [33,34,35], thus offering advantages for monitoring drought, desertification, and dust occurrence [6,36,37].
A satellite-based aridity index (SbAI) was proposed [38], which is an indicator of the thermal inertia calculated from day–night LST differences, to monitor the extent of the aridity or wetness. The SbAI has been validated and used to identify degraded land and drought in north-east Asia and globally [6,36,39]. Here, we used the ECT method, replacing temperature with the SbAI calculated from satellite data, to monitor vegetation growth and consumed SWC in north-east Asia, focusing on the natural grasslands of China and Mongolia.

2. Methods

2.1. Study Areas and Periods

2.1.1. Shenmu, China

The observation site was an experimental station in the Liudaogou River Basin (7 km2) in Shenmu District, Shaanxi Province, China (38°47′ N, 110°21′ E, 1224 m above sea level; Figure 1) [40]. The average annual temperature from 1957 to 1989 was 8.4 °C, with an average annual rainfall of 437 mm [41]. The region is semi-arid [42].
Data of the SWC in Shenmu, China, were recorded from June 2004 to March 2010, but we used the SWC data from 2005 to 2007, a timespan when the data quality (data within the operating life (=5 years) of the SWC instruments) and a continuous data acquisition (without missing data due to some troubles) could be assured. This period fell during the Chinese Government’s Grain for Green Program, intended to reduce the farmland area and promote greening, so the effects of human-managed grazing were low. The station was covered by Stipa bungeana Trin., with the root zone within 34 cm of the surface.
The soil texture is a sandy loam (48% sand, 39% silt, and 13% clay). The volumetric SWC was 0.43 m3 m−3, field capacity was 0.19 m3 m−3, and permanent wilting point was 0.05 m3 m−3. The SWC at ten depths (Delta-T; ML2X; 6, 10, 18, 26, 34, 42, 50, 58, 66 and 100 cm) and rainfall (Davis Inst. Corp; Rain Collector 2) were sampled every minute with a datalogger (Campbell Scientific; CR10X), except for the SWC at 12:00 and 24:00 Beijing standard time (BST). A soil moisture sensor was installed vertically from the soil surface at a depth of 6 cm; therefore, the averaged SWC up to 6 cm was measured. Except for this installation at 6 cm, the sensors were set horizontally. The averaged SWC within the root zone (0–34 cm) was used for analysis. The total rainfall was 308 mm in 2005, 385 mm in 2006, and 481 mm in 2007.

2.1.2. Bayan-Unjuul, Mongolia

We selected a validation site at Bayan-Unjuul, northern Mongolia (47°02′38.7″ N, 105°56′56.1″ E, 1200 m above sea level; Figure 1), where a 300-m × 300-m fence was erected in June 2004 to exclude livestock [43]. The selected period to validate the proposed method was 2017 to 2019, the most recent three years; 2017 featured a prolonged period of severe dry weather, when an estimated 80% of the country was affected by drought [44], and 2018 featured heavy rainfall [45].
The average annual temperature from 1995 to 2005 was 0.1 °C and the average annual rainfall was 163 mm [43]. Most rainfall is concentrated from May to August (124 mm). The soils are classified as Kastanozems. The soil texture is sandy loam. The field capacity was 0.20 m3 m−3 and the permanent wilting point was 0.05 m3 m−3, close to those at Shenmu. The station was covered by short grasses, including Stipa krylovii and Agropyron cristatum, and the shrub species Caragana microphylla. The root zone was limited to 25 cm by a stiff CaCO3 layer (caliche) at that depth [12].

2.2. Data and Calculation of the NDVI and SbAI

The NDVI is defined as follows [46]:
NDVI = ρ NIR ρ RED ρ NIR + ρ RED   ,
where ρNIR is surface reflectance from MODIS band 2 (841 to 876 nm) and ρRED is that from MODIS band 1 (620 to 670 nm).
The 16-day NDVI was calculated with MOD09A1 8-day global surface-reflectance data in the near-infrared and visible (red) wavelength channels at a 500-m resolution. The 8-day NDVI with the highest score (i.e., little influence of cloud, sensor angle, and aerosol concentration) was selected as the 16-day NDVI. The spatial resolution was changed to 1 km by averaging 4 cells. Finally, the monthly NDVI was calculated by averaging the 16-day NDVI.
The SbAI is defined as follows [38]:
SbAI = T s R s ,
where
R s = 1 r S 0 cos θ c .
Here, ∆Ts is the difference in LST between day and night, and Rs is the absorbed solar radiation calculated from broadband albedo r, the solar constant S0 (1367 W m−2), and the solar zenith angle at the Sun’s apex, θc. ∆Ts is an indicator of both surface dryness and other factors, especially solar radiation. To compensate for the effect of solar radiation, ∆Ts is divided by Rs in Equation (2). For a dry surface, SbAI is large because ∆Ts is large, and ∆Ts is large because the low land surface wetness causes the thermal inertia to be low.
The concept and assumptions of the relationships among SWC, vegetation density, and SbAI are shown in Figure 2a. When the grass density is low, the little rainfall from March to June maintains the SWC in only the surface soil, not the deeper soil [40]. Therefore, SbAI indicates mainly the thermal inertia of the surface soil. When heavy rainfall from July to August infiltrates down to the root zone, the grass density becomes high [40]. Water uptake by the grass roots becomes active, resulting in an increased water content of the grass and thus increased thermal inertia, which is monitored by SbAI.
SbAI was calculated by the method of [38] using the MODIS data products MOD09CMG (daily global surface-reflectance) and MOD11C1 (daily global day/night LST). The spatial resolution of both products was 0.05°. For both the NDVI and SbAI, cells covering each site were used for analysis.
To plot the spatial distribution of the precipitation in north-east Asia, we used the CPC Global Unified Precipitation data provided by NOAA/OAR/ESRL PSL (Boulder, CO, USA; https://psl.noaa.gov/). The spatial resolution was 0.5° × 0.5°.

2.3. Effective Cumulative Reciprocal SbAI

ECT is calculated as follows:
E C T = T d T b ,
where Td is the daily averaged temperature (°C) and Tb is threshold temperature (°C) determined by the plant species.
Numerous studies report the relationships among rainfall, SWC, and vegetation growth (Figure 2b, red arrows) in natural grasslands of north-east Asia [9,19,43,47,48]. Starting with the ECT method, we replaced temperature with the effective cumulative reciprocal of SbAI (CRSbAI) obtained from satellite data (Figure 2b) to monitor the SWC consumption and vegetation growth as follows:
C R S b A I = k = i k = j S b A I 1 k 1 / 0.03 ,
where SbAI−1 is the reciprocal of SbAI, because SbAI decreases with increasing thermal inertia; i is the initial month of growth (March in Shenmu, April in Bayan-Unjuul); j is the month when growth reaches the maximum (August at both sites); k is month; and 0.03 is the threshold value of the SbAI in extreme dry conditions in which wind erosion is likely [6,38].
As the SbAI cannot always be derived owing to the effect of cloud cover, we used a 16-day moving average to obtain the continuous, daily SbAI.
This study asked the following two questions (Figure 2b, blue arrows):
  • Can the CRSbAI replace the cumulative decreased SWC (CdSWC)?
  • Is the CRSbAI related to vegetation growth?
We developed these concepts at Shenmu and verified them at Bayan-Unjuul. Indices which represent vegetation biomass should be based on observed above-ground dry matter. However, we used the NDVI, which is strongly correlated with the observed dry matter in grasslands of north-east Asia [48,49,50].

3. Results

3.1. Seasonal Change in the NDVI, SWC, and SbAI in Shenmu

In Shenmu, the NDVI increased gradually to a maximum of 0.35 to 0.4 in August (Figure 3). The total rainfall from March to August was 215 mm in 2005, 309 mm in 2006, and 303 mm in 2007. Although the total rainfall was higher in 2007 than in 2005, the NDVI was lowest in August 2007. This was caused by rainfall intensity and thus a higher SWC in June to July (DOY 152 to 212) (Figure 4). In June to July, the SWC approached field capacity (0.19 m3 m−3) and stayed high in 2005 and 2006, but was lower (<0.15 m3 m−3) in 2007. Concentrated rainfall in June to July strongly increased subsequent vegetation biomass in the grasslands of Mongolia and China [19,47].
Although the seasonal change in the daily SbAI varied, the trend of the 16-day moving-average SbAI mirrored that of the SWC (one was high when the other was low; Figure 4). The SbAI had a significant correlation with SWC (r = 0.405, p < 0.05). Here, r is the correlation coefficient, and p is the p-value.

3.2. Relationships among CRSbAI, CdSWC, and ∆NDVI in Shenmu

The relationships of CRSbAI (daily) with CdSWC within the root zone and ∆NDVI were examined in Shenmu from 2005 to 2007 (Figure 5). Here, “CRSbAI (daily)” means the CRSbAI from March to March, to April, to May, to June, to July, and to August, calculated from the daily SbAI with the 16-day moving average. The CdSWC covers March to March, to April, to May, to June, to July, and to August. ∆NDVI is the NDVI in March, April, May, June, July, and August minus the NDVI in February, the month with no grass growth.
The CRSbAI was significantly correlated with the CdSWC (r = 0.984, RMSE = 8.7 mm, p < 0.001) and ∆NDVI (r = 0.966, RMSE = 0.022, p < 0.001). Here, RMSE is the root-mean-squared error. Although the rainfall amount and intensity differed among the years, CRSbAI integrated them into a single linear relationship. An index of thermal inertia based on LST was showed to be most highly related to SWC within the root zone during grass growth in China [51]; they concluded that wet soil may cause both high transpiration and evaporation at the ground surface. These findings and our results in Figure 5 may validate the concept shown in Figure 2.
The final CdSWC values were 147 mm in 2005, 169 mm in 2006, and 142 mm in 2007 (Figure 5, blue arrows). The CdSWC is the evapotranspiration (ET), and its ratio to the total rainfall during study period was 68% in 2005, 55% in 2006, and 47% in 2007. These ratios are close to the results obtained in past studies in grasslands of California (54–59%; [52]), South Africa (40–76%; [53]), and Mongolia (66%; [54]).
The CRSbAI (monthly) was also significantly and highly correlated with CdSWC (r = 0.981, RMSE = 9.4 mm, p < 0.001) and ∆NDVI (r = 0.966, RMSE = 0.022, p < 0.001) (Figure 6). Here, the CRSbAI (monthly) was calculated from the monthly average observed SbAI. Therefore, the CRSbAI (monthly) can be calculated easily from the observed data without statistical processing (such as a 16-day moving average).

3.3. Validation of the CRSbAI Method in Bayan-Unjuul

In Bayan-Unjuul, the seasonal NDVI was steady from March to June and then increased rapidly to August (Figure 7). This is a general trend in grass growth in Mongolia, caused by the rainfall concentrated in June and July [19,47]. In Shenmu, however, the NDVI increased slowly from March to August (Figure 3) because of frequent rainfalls. This is a crucial difference in grass growth between the sites. The severe drought in 2017 [44] explains the lower NDVI in July and August 2017 when compared to 2018 and 2019.
Both the CRSbAI (daily) and CRSbAI (monthly) agreed well with ∆NDVI (daily: r = 0.982, RMSE = 0.024, p < 0.001) (monthly: r = 0.960, RMSE = 0.030, p < 0.001) (Figure 8), as at Shenmu (Figure 5 and Figure 6). The calculation began in April (March in Shenmu) because the soil was still frozen in March. From April to June, when the rainfall was low, the small CRSbAI corresponded to the small ∆NDVI. From July, however, during the rainy season, as the CRSbAI increased rapidly, so too did the ∆NDVI. This validation result in Bayan-Unjuul indicates that the CRSbAI responds well to the water condition of the land surface assumed in Figure 2a.
NDVI peaked at around 0.3 in Bayan-Unjuul (Figure 8); the corresponding CRSbAI values were 1300 to 1600 daily and 50 to 70 monthly. Although the ∆NDVI peaked at around the same value in Shenmu (Figure 5 and Figure 6), the corresponding CRSbAI values were 3000 to 3500 daily and 100 to 130 monthly, double those in Bayan-Unjuul. This means that the grass in Bayan-Unjuul has a high water-use efficiency by fitting its physiological activity to the late rainfall.

3.4. Simulation of CdSWC in Bayan-Unjuul

As the sandy loam and dominant grass species were comparable between Shenmu and Bayan-Unjuul, we assumed that the only serious difference between the sites was rainfall. The daily and monthly CdSWC values in Shenmu (Table 1) were calculated as follows:
C d S W C d a i l y = 8 · 10 6 C R S b A I d a i l y 2 + 0.0739 · C R S b A I d a i l y   0 < CRSbAI daily < 3500
C d S W C m o n t h l y = 0.0035 · C R S b A I m o n t h l y 2 + 1.8271 · C R S b A I m o n t h l y   0   < CRSbAI monthly < 140 .
The daily and monthly values were close, differing by only 10 mm among years.
The CdSWC was lowest in 2017 (54 to 66 mm), owing to drought, and that in 2018 was highest (95 to 108 mm), owing to heavy rainfall [45]. As the CdSWC (monthly) is the ET, the ET was 66 mm in 2017, 108 mm in 2018, and 93 mm in 2019. The Institute of Meteorology, Hydrology, and Environment records rainfall in Bayan-Unjuul. From April to August it recorded 122 mm in 2017, 154 mm in 2018, and 159 mm in 2019. On the Mongolian dry steppe near Bayan-Unjuul (Kherlenbayan-Ulaan: 47°12.838′ N, 105°44.240′ E, 1235 m above sea level), the ratio of ET to precipitation is 66% [54]. At that site, the soil is a Kastanozem and the predominant vegetation is S. krylovii, the same as at Byan-Unjuul. When this ratio was used, the ET was calculated as 81 mm in 2017, 102 mm in 2018, and 105 mm in 2019, close to the ET as CdSWC (root-mean-squared error of ±11 mm)

4. Discussion

Here, we will discuss the applicability evaluation of the CRSbAI over north-east Asia from 2017 to 2019, the three most recent years; that is, the spatial distribution of the CRSbAI and CdSWC will be examined using the calculation method obtained in Shenmu and Bayan-Unjuul, and validated using past published data.
Figure 9 maps the distributions of the CRSbAI (monthly) (cumulated from April to August) and NDVI in August from 2017 to 2019. Although the response of ∆NDVI to the CRSbAI differed among years (Figure 5, Figure 6 and Figure 8), the distribution of the CRSbAI corresponded well with that of the NDVI:
  • NDVI < 0.2 ≡ CRSbAI < 40;
  • 0.2 < NDVI < 0.4 ≡ 40 < CRSbAI < 80;
  • 0.4 < NDVI < 0.6 ≡ 80 < CRSbAI < 120.
In particular, the Gobi Desert, south central Mongolia, Inner Mongolia, the Loess Plateau of China, and irrigated farmland along the Yellow River were well delineated by the CRSbAI. The area where the CRSbAI < 60 (purple) was large in 2017, when severe drought prevailed, and the area with a CRSbAI ≈ 80 (blue) increased in 2018 owing to the heavy rainfall. Thus, the CRSbAI reflected differences in rainfall intensity.
Figure 10 shows the distribution of CdSWC calculated with Equation (7) and the total precipitation (Pr) from April to August 2017 to 2019, using the CPC Global Unified Precipitation data. Since the range of CRSbAI (monthly) was approximately 0 to 150 mm (Figure 6), areas where the CRSbAI were >150 mm were not available. However, areas where the CRSbAI were <150 mm corresponded closely to the Pr distribution:
  • 0 < CdSWC < 50 ≡ 0 < Pr < 80;
  • 50 < CdSWC < 70 ≡ 80 < Pr < 130;
  • 70 < CdSWC < 110 ≡ 130 < Pr < 210;
  • 110 < CdSWC ≡ 210 < Pr.
It is of note that the distribution of the CdSWC corresponded closely to the arid or semi-arid regions [42,55], mostly the Gobi Desert and natural grasslands where dust outbreaks arise [56,57,58].
As the CdSWC is the ET, the ET in the Gobi Desert was <30 mm during April to August. The ET was sensitive to rainfall in south-central Mongolia, Inner Mongolia (China), and the Loess Plateau of China (Figure 10). Drought in 2017 caused a low ET and NDVI, not only over a wide area in Mongolia, but also in Inner Mongolia. Thus, we can easily predict that such an effect of drought can trigger outbreaks of dust in these areas.
It is not clear whether Equation (7) can be applied to all areas in Figure 10. However, the ET observed by the eddy covariance method from May to August 2014 in Tsogt-Ovoo (desert steppe; 44°23′04″ N, 105°16′59″ E, 1232 m above sea level) is one reason for the validation [59]. The annual rainfall in 2014 was 72 mm and normal (the average from 1981 to 2010 was 75 mm), and the observed ET was 50 mm. This value is close to the ET values of the cell including this point (23 mm in 2017, 75 mm in 2018, and 38 mm in 2019). Another validation result was reported by [54], who calculated the ET by the eddy covariance method in Kherlenbayan-Ulaan, Mongolia, in 2003; the average annual rainfall from 1993 to 2002 was 181 mm, and the value in 2003 (248 mm) was larger than normal; the observed ET was 120 mm from April to August, only slightly higher than our CdSWC value of 103 mm in 2018, when the rainfall was heavier than normal.
The ratio of the CdSWC to Pr was 52% to 63% across their spatial distribution (Figure 10), close to the values in the dry steppe [54] and other grasslands [52,53].

5. Conclusions

We used a satellite-based aridity index (SbAI) calculated only from satellite data, the effective cumulative reciprocal SbAI (CRSbAI), to monitor the SWC and grass growth in north-east Asia, calibrating it in natural grasslands in Shenmu, China, and validating its applicability in Bayan-Unjuul, Mongolia. We asked the following questions:
  • Can the CRSbAI replace CdSWC?
  • Is the CRSbAI related to grass growth and represented by ∆NDVI?
The CRSbAI was significantly correlated with CdSWC in Shenmu. Although the rainfall differed among years, the CRSbAI integrated the rainfall values in a single linear relationship. As the CdSWC calculated from the CRSbAI is the ET, the ET in Bayan-Unjuul was 66 mm in 2017, 108 mm in 2018, and 93 mm in 2019. As the ET is 66% of the precipitation in the Mongolian steppe, the estimated ET was 81, 102, and 105 mm, respectively. These respective values are close (the root-mean-squared error is within ±11 mm). The ratio of CdSWC to Pr across north-east Asia was 52% to 63%, close to the ratios in other grasslands reported in past studies.
The CRSbAI was strongly related to ∆NDVI at both Shenmu and Bayan-Unjuul, although grass growth responded differently owing to differences in rainfall intensity; also, the spatial distribution of the CRSbAI corresponded well to that of the NDVI in August. This result suggests that the CRSbAI can be used to predict the distribution of the NDVI qualitatively. To predict the ∆NDVI quantitatively, their relationship should be examined in the respective target regions.
Verification of our assumptions (above two questions) is not yet complete. However, the CRSbAI could offer an index of water consumption and grass growth, and may be useful for development of early warning and monitoring systems for drought. We hope that the usefulness of our method will be confirmed by other researchers, and that it will serve as the basis for an improved system based on remote sensing techniques that will promote sustainable development in north-east Asia.

Author Contributions

Conceptualization, R.K.; methodology, R.K.; software, M.M.; validation, R.K. and M.M.; formal analysis, R.K.; resources, M.M.; data curation, R.K. and M.M.; writing—original draft preparation, R.K.; writing—review and editing, R.K.; visualization, R.K.; supervision, R.K.; project administration, R.K.; funding acquisition, R.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by a Grant-in-Aid for Scientific Research, grant number KAKENHI 19H04239.

Acknowledgments

We sincerely thank Amarsaikhan Davaadorj of the Institute of Meteorology, Hydrology, and Environment, Mongolia, and Batjargal Buyantogtoh of Tottori University for preparing the precipitation data from Bayan-Unjuul. We appreciate the invaluable comments from five reviewers and the academic editor of this paper.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. IPCC. Climate change 2013: The physical science basis. In Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; Stocker, T.F.D., Qin, G.-K., Plattner, M., Tignor, S.K., Allen, J., Boschung, A., Nauels, Y., Xia, V.B., Midgley, P.M., Eds.; Cambridge University Press: Cambridge, UK, 2013; p. 1535. [Google Scholar]
  2. Feng, S.; Fu, Q. Expansion of global drylands under a warming climate. Atmos. Chem. Phys. 2013, 13, 10081–10094. [Google Scholar] [CrossRef] [Green Version]
  3. Koutroulis, A.G. Dryland changes under different levels of global warming. Sci. Total Environ. 2019, 655, 482–511. [Google Scholar] [CrossRef] [PubMed]
  4. Bréon, F.M.; Luneau, G. Atlas du Climat face Aux Défis du Réchauffement; Éditions Autrement: Paris, France, 2015; p. 96. [Google Scholar]
  5. Cherlet, M.; Hutchinson, C.; Reynolds, J.; Hill, J.; Sommer, S.; Von Maltitz, G. (Eds.) World Atlas of Desertification; Publication Office of the European Union: Luxembourg, 2018; p. 248. [Google Scholar]
  6. Kimura, R. Validation and application of the monitoring method for degraded land area based on a dust erodibility in eastern Asia. Int. J. Remote Sens. 2017, 38, 4553–4564. [Google Scholar] [CrossRef]
  7. Hessl, A.E.; Anchukaitis, K.J.; Jelsema, C.; Cook, B.; Byambasuren, O.; Leland, C.; Nachin, B.; Pederson, N.; Tian, H.; Hayles, L.A. Past and future drought in Mongolia. Sci. Adv. 2018, 4, e1701832. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Natsagadorj, L.; Dulamsuren, J. Some aspects of assessment of the dzud phenomena. Pap. Meteorol. Hydrol. 2001, 23, 3–18. [Google Scholar]
  9. Iwasaki, H. NDVI prediction over Mongolian grassland using GSMaP precipitation data and JRA-25/JCDAS temperature data. J. Arid Environ. 2009, 73, 557–562. [Google Scholar] [CrossRef]
  10. Du, C.; Shinoda, M.; Tachiiri, K.; Nandintsetseg, B.; Komiyama, H.; Matsushita, S. Mongolian herders’ vulnerability to dzud: A study of record livestock mortality levels during the severe 2009/2010 winter. Nat. Hazards 2018, 92, S3–S17. [Google Scholar] [CrossRef]
  11. Kimura, R.; Bai, L.; Wang, J. Relationships among dust outbreaks, vegetation cover, and surface soil water content on the Loess Plateau of China, 1999–2000. Catena 2009, 77, 292–296. [Google Scholar] [CrossRef]
  12. Kimura, R.; Shinoda, M. Spatial distribution of threshold wind speeds for dust outbreaks in northeast Asia. Geomorphology 2010, 114, 319–325. [Google Scholar] [CrossRef]
  13. Kimura, R. Factors contributing to dust storms in source regions producing the yellow-sand phenomena observed in Japan from 1993 to 2002. J. Arid Environ. 2012, 80, 40–44. [Google Scholar] [CrossRef]
  14. Fang, Z.Y.; Zhu, F.K.; Jiang, J.X.; Qian, Z.A. Studies on Dust Storm in China; China Meteorological Press: Beijing, China, 1997; p. 158. [Google Scholar]
  15. Kwon, H.J.; Cho, S.H.; Chun, Y.; Lagarde, F.; Pershagen, G. Effects of the Asian dust events on daily mortality in Seoul, Korea. Environ. Res. Sect. 2002, A 90, 1–5. [Google Scholar] [CrossRef]
  16. Onishi, K.; Sekiyama, T.T.; Nojima, M.; Kurosaki, Y.; Fujitani, Y.; Otani, S.; Maki, T.; Shinoda, M.; Kurozawa, Y.; Yamagata, Z. Prediction of health effects of cross-border atmospheric pollutants using an aerosol forecast model. Environ. Int. 2018, 117, 48–56. [Google Scholar] [CrossRef]
  17. Jones, H.G. Plants and Microclimate: A Quantitative Approach to Environmental Plant Physiology, 3rd ed.; Cambridge University Press: Cambridge, UK, 2014; p. 423. [Google Scholar]
  18. MLIT. 2007. Available online: https://www.mlit.go.jp/river/pamphlet_jirei/bousai/saigai/kiroku/suigai/suigai_3-1-1.html/ (accessed on 24 August 2020).
  19. Kimura, R. Relationships between NDVI, rainfall, and soil water content of the natural grassland in the Liudaogou river basin of the Loess Plateau. J. Arid Land Stud. 2011, 21, 1–6. [Google Scholar]
  20. Jones, H.G.; Vaughan, R.A. Remote Sensing of Vegetation; Oxford University Press: New York, NY, USA, 2010; p. 353. [Google Scholar]
  21. Kimura, R.; Moriyama, M. Determination by MODIS satellite-based methods of recent global trends in land surface aridity and degradation. J. Agric. Meteorol. 2019, 75, 153–159. [Google Scholar] [CrossRef] [Green Version]
  22. Cogato, A.; Meggio, F.; De Antoni Migliorati, M.; Marinello, F. Extreme Weather Events in Agriculture: A Systematic Review. Sustainability 2019, 11, 2547. [Google Scholar] [CrossRef] [Green Version]
  23. Kanda, E.; Torigoe, Y.; Kobayashi, T. A simple model to predict the developmental stages of rice panicles using the effective accumulative temperature. Jpn. J. Crop. Sci. 2002, 71, 394–402. [Google Scholar] [CrossRef] [Green Version]
  24. Yasumoto, S.; Maki, N.; Kojima, M.; Ohshita, Y. Changes in developmental duration of direct-seeded rice in a well-drained paddy field in response to late planting. Plant Prod. Sci. 2017, 20, 279–287. [Google Scholar] [CrossRef] [Green Version]
  25. Li, B.; Ishii, Y.; Idota, S.; Tobisa, M.; Niimi, M.; Yang, Y.; Nishimura, K. Yield and quality of forages in a triple cropping system in southern Kyushu, Japan. Agronomy 2019, 9, 277. [Google Scholar] [CrossRef] [Green Version]
  26. Horie, T.; Nakagawa, H. Modelling and prediction of developmental process in rice. Jpn. J. Crop Sci. 1990, 59, 687–695. [Google Scholar] [CrossRef] [Green Version]
  27. Homma, K.; Maki, M.; Hirooka, Y. Development of a rice simulation model for remote-sensing (SIMRIW-RS). J. Agric. Meteorol. 2017, 73, 9–15. [Google Scholar] [CrossRef] [Green Version]
  28. Yoshida, R.; Fukui, S.; Fukui, S.; Yamazaki, T. Applicability of meteorological ensemble forecasting to predict summer cold damage in rice growth. J. Agric. Meteorol. 2020, 76, 128–139. [Google Scholar] [CrossRef]
  29. Jackson, T.J. Measuring surface soil moisture using passive microwave remote sensing. Hydrol. Process. 1993, 7, 139–152. [Google Scholar] [CrossRef]
  30. Koike, T.; Nakamura, Y.; Kaihotsu, I.; Davva, N.; Matsuura, K.; Tamagawa, K.; Fujii, H. Development of an advanced microwave scanning radiometer (AMSR-E) algorithm of soil moisture and vegetation water content. Annual Journal of Hydraulic Engineering. JSCE 2004, 48, 217–222. [Google Scholar]
  31. Fujii, H.; Koike, T.; Imaoka, K. Improvement of the AMSR-E algorithm for soil moisture estimation by introducing a fractional vegetation coverage dataset derived from MODIS data. J. Remote Sens. Soc. Jpn. 2009, 29, 282–292. [Google Scholar]
  32. Kaihotsu, I.; Asanuma, J.; Aida, K.; Oyunbaatar, D. Evaluation of the AMSR2 L2 soil moisture product of JAXA on the Mongolian Plateau over seven years (2012–2018). SN Appl. Sci. 2019, 1, 1477. [Google Scholar] [CrossRef] [Green Version]
  33. Verhoef, A. Remote estimation of thermal inertia and soil heat flux for bare soil. Agric. For. Meteorol. 2004, 123, 221–236. [Google Scholar] [CrossRef]
  34. Verstraeten, W.W.; Veroustraete, F.; Van der Sande, C.J.; Grootaers, I.; Feyen, J. Soil moisture retrieval using thermal inertia, determined with visible and thermal space borne data, validated for European forests. Remote Sens. Environ. 2006, 101, 299–314. [Google Scholar] [CrossRef]
  35. Matsushima, D.; Kimura, R.; Shinoda, M. Soil moisture estimation using thermal inertia: Potential and sensitivity to data conditions. J. Hydrometeorol. 2012, 13, 638–648. [Google Scholar] [CrossRef]
  36. Kimura, R. Global detection of aridification or increasing wetness in arid regions from 2001 to 2013. Nat. Hazards 2020, 103, 2261–2276. [Google Scholar] [CrossRef]
  37. Matsushima, D.; Kimura, R.; Kurosaki, Y.; Ganzorig, U.; Shinoda, M. A method for estimating the threshold wind speed for dust emissions as a function of soil moisture. Bound. Layer Meteorol. 2020, 175, 237–257. [Google Scholar] [CrossRef]
  38. Kimura, R.; Moriyama, M. Application of a satellite-based aridity index in dust source regions of northeast Asia. J. Arid Environ. 2014, 109, 31–38. [Google Scholar] [CrossRef]
  39. Kimura, R. Global distribution of degraded land area based on dust erodibility determined from satellite data. Int. J. Remote Sens. 2018, 39, 5859–5871. [Google Scholar] [CrossRef]
  40. Kimura, R.; Fan, J.; Zhang, X.; Takayama, N.; Kamichika, M.; Matsuoka, N. Evapotranspiration over the grassland field in the Liudaogou basin of the Loess Plateau, China. Acta Oecologica 2006, 29, 45–53. [Google Scholar] [CrossRef]
  41. Academia Sinica and Ministry of Water Resources. Memoir of Northwestern Institute of Soil and Water Conservation; Shaanxi Scientific and Technological Press: Xi’an, China, 1993; Volume 18, p. 144. (In Chinese)
  42. UNEP. World Atlas of Desertification; Arnold: London, UK, 1997; p. 182. [Google Scholar]
  43. Shinoda, M.; Nachinshonhor, G.U.; Nemoto, M. Impact of drought on vegetation dynamics of the Mongolian steppe: A field experiment. J. Arid Environ. 2010, 74, 63–69. [Google Scholar] [CrossRef]
  44. FAO. Special Report—FAO/WFP Crop and Livestock Assessment Mission to Mongolia; FAO/WFP: Rome, Italy, 2017; p. 43. [Google Scholar]
  45. JMA (Japan Meteorological Agency). 2018. Available online: https://www.jma.go.jp/jma/press/1812/21d/worldclim2018.html/ (accessed on 21 August 2020).
  46. Rouse, J.W.; Haas, R.H.; Schell, J.A.; Deering, D.W.; Harlan, J.C. Monitoring the Vernal Advancement and Retrogradation (Green Wave Effect) of Natural Vegetation; NASA/GSFC Final report; NASA Goddard Space Flight Center : Greenbelt, MD, USA, 1974. [Google Scholar]
  47. Miyazaki, S.; Yasunari, T.; Miyamoto, T.; Kaihotsu, I.; Davaa, G.; Oyunbaatar, D.; Natsagdorj, L.; Oki, T. Agrometeorological conditions of grassland vegetation in central Mongolia and their impact for leaf area growth. J. Geophys. Res. 2004, 109, D22106. [Google Scholar] [CrossRef]
  48. Nandintsetseg, B.; Shinoda, M.; Kimura, R.; Ibaraki, Y. Relationship between soil moisture and vegetation activity in the Mongolian steppe. SOLA 2010, 6, 029–032. [Google Scholar] [CrossRef] [Green Version]
  49. Kawamura, K.; Akiyama, T.; Yokota, H.; Tsutsumi, M.; Yasuda, T.; Watanabe, O.; Wang, S. Comparing MODIS vegetation indices with AVHRR NDVI for monitoring the forage quantity and quality in Inner Mongolia grassland, China. Jpn. Soc. Grassl. Sci. 2005, 51, 33–40. [Google Scholar] [CrossRef]
  50. Nakano, T.; Bavuudorj, G.; Nachinshonhor, G.U.; Shinoda, M. Monitoring aboveground biomass in semiarid grasslands using MODIS images. J. Agric. Meteorol. 2013, 69, 33–39. [Google Scholar] [CrossRef] [Green Version]
  51. Mohamed, A.A.; Kimura, R.; Shinoda, M.; Moriyama, M. Diurnal surface temperature difference index derived from ground-based meteorological measurements for assessment of moisture availability. J. Arid Environ. 2011, 75, 156–163. [Google Scholar] [CrossRef]
  52. Baldocchi, D.D.; Xu, L.; Kiang, N. How plant functional-type, weather, seasonal drought and soil physical properties alter water and energy fluxes of an oak-grass savanna and an annual grassland. Agric. Forest Meteorol. 2004, 123, 13–39. [Google Scholar] [CrossRef] [Green Version]
  53. Everson, C.S. The water balance of a first order catchment in the montane grasslands of South Africa. J. Hydrol. 2001, 241, 110–123. [Google Scholar] [CrossRef]
  54. Li, S.G.; Asanuma, J.; Kotani, A.; Davaa, G.; Oyunbaatar, D. Evapotranspiration from a Mongolian steppe under grazing and its environmental constraints. J. Hydrol. 2007, 333, 133–143. [Google Scholar] [CrossRef]
  55. Kimura, R.; Moriyama, M. Recent trends of annual aridity indices and classification of arid regions with satellite-based aridity indices. Remote Sens. Earth Syst. Sci. 2019, 2, 88–95. [Google Scholar] [CrossRef]
  56. Sun, J.; Zhang, M.; Liu, T. Spatial and temporal characteristics of dust storms in China and its surrounding regions, 1960–1999: Relations to source area and climate. J. Geophys. Res. 2001, 106, 10325–10333. [Google Scholar] [CrossRef]
  57. Miao, L.; Luan, Y.; Luo, X.; Liu, Q.; Moore, J.C.; Nath, R.; He, B.; Zhu, F.; Cui, X. Analysis of the phenology in the Mongolian plateau by inter-comparison of global vegetation datasets. Remote Sens. 2013, 5, 5193–5208. [Google Scholar] [CrossRef] [Green Version]
  58. Kimura, R. Satellite-based mapping of dust erodibility in northeast Asia. Nat. Hazards 2016, 92, 19–25. [Google Scholar] [CrossRef]
  59. Kimura, R.; Abulaiti, A.; Mano, M.; Matsushima, D. Surface heat flux analysis in Gobi Desert steppe, Mongolia—An observation study. SOLA 2016, 12, 175–180. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Locations of Shenmu, China, and Bayan-Unjuul, Mongolia, with the distribution of the NDVI in August 2019.
Figure 1. Locations of Shenmu, China, and Bayan-Unjuul, Mongolia, with the distribution of the NDVI in August 2019.
Remotesensing 12 03556 g001
Figure 2. Concepts of this study. (a) Relationships among the vegetation density, soil water content (SWC), and satellite-based aridity index (SbAI). Red upward arrow, thermal inertia; green dashed box, root zone; green solid box, grass canopy. (b) Relationships among the cumulative rainfall, cumulative decreased SWC (CdSWC), effective cumulative 1/SbAI (CRSbAI), and vegetation growth (∆NDVI). Red arrows, relationship among cumulative rainfall, CdSWC, and vegetation growth; blue arrows, research concept: (1) Can the CRSbAI replace CdSWC? (2) Is the CRSbAI related to vegetation growth?
Figure 2. Concepts of this study. (a) Relationships among the vegetation density, soil water content (SWC), and satellite-based aridity index (SbAI). Red upward arrow, thermal inertia; green dashed box, root zone; green solid box, grass canopy. (b) Relationships among the cumulative rainfall, cumulative decreased SWC (CdSWC), effective cumulative 1/SbAI (CRSbAI), and vegetation growth (∆NDVI). Red arrows, relationship among cumulative rainfall, CdSWC, and vegetation growth; blue arrows, research concept: (1) Can the CRSbAI replace CdSWC? (2) Is the CRSbAI related to vegetation growth?
Remotesensing 12 03556 g002
Figure 3. Monthly NDVI in Shenmu in 2005 to 2007.
Figure 3. Monthly NDVI in Shenmu in 2005 to 2007.
Remotesensing 12 03556 g003
Figure 4. Seasonal change in the SbAI and SWC averaged within the root zone in Shenmu in 2005 to 2007.
Figure 4. Seasonal change in the SbAI and SWC averaged within the root zone in Shenmu in 2005 to 2007.
Remotesensing 12 03556 g004
Figure 5. Relationships of CRSbAI (daily) with CdSWC and ∆NDVI in Shenmu from 2005 to 2007. Blue arrows show the final CdSWC values.
Figure 5. Relationships of CRSbAI (daily) with CdSWC and ∆NDVI in Shenmu from 2005 to 2007. Blue arrows show the final CdSWC values.
Remotesensing 12 03556 g005
Figure 6. Relationships of CRSbAI (monthly) with CdSWC and ∆NDVI in Shenmu from 2005 to 2007.
Figure 6. Relationships of CRSbAI (monthly) with CdSWC and ∆NDVI in Shenmu from 2005 to 2007.
Remotesensing 12 03556 g006
Figure 7. Monthly NDVI in Bayan-Unjuul from 2017 to 2019.
Figure 7. Monthly NDVI in Bayan-Unjuul from 2017 to 2019.
Remotesensing 12 03556 g007
Figure 8. Relationships of CRSbAI (daily) and CRSbAI (monthly) with ∆NDVI in Bayan-Unjuul from 2017 to 2019.
Figure 8. Relationships of CRSbAI (daily) and CRSbAI (monthly) with ∆NDVI in Bayan-Unjuul from 2017 to 2019.
Remotesensing 12 03556 g008
Figure 9. Spatial distributions of the CRSbAI (monthly) and NDVI in August in north-east Asia from 2017 to 2019. White areas indicate where cloud cover prevented calculation.
Figure 9. Spatial distributions of the CRSbAI (monthly) and NDVI in August in north-east Asia from 2017 to 2019. White areas indicate where cloud cover prevented calculation.
Remotesensing 12 03556 g009
Figure 10. Spatial distribution of the CdSWC (mm) calculated from the CRSbAI (monthly) and precipitation (Pr, mm) from April to August 2017 to 2019, using the CPC Global Unified Precipitation data. White areas indicate where the CRSbAI > 150 mm (therefore, data not available).
Figure 10. Spatial distribution of the CdSWC (mm) calculated from the CRSbAI (monthly) and precipitation (Pr, mm) from April to August 2017 to 2019, using the CPC Global Unified Precipitation data. White areas indicate where the CRSbAI > 150 mm (therefore, data not available).
Remotesensing 12 03556 g010
Table 1. The CdSWC to August calculated from the CRSbAI and Equations (6) and (7) with the ∆NDVI (August minus March) in Bayan-Unjuul from 2017 to 2019.
Table 1. The CdSWC to August calculated from the CRSbAI and Equations (6) and (7) with the ∆NDVI (August minus March) in Bayan-Unjuul from 2017 to 2019.
ΔNDVICRSbAI
(Daily)
CRSbAI
(Monthly)
CdSWC (mm)
(Daily)
CdSWC (mm)
(Monthly)
20170.219794395466
20180.27915496895108
20190.2811302578393
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kimura, R.; Moriyama, M. Use of a Satellite-Based Aridity Index to Monitor Decreased Soil Water Content and Grass Growth in Grasslands of North-East Asia. Remote Sens. 2020, 12, 3556. https://0-doi-org.brum.beds.ac.uk/10.3390/rs12213556

AMA Style

Kimura R, Moriyama M. Use of a Satellite-Based Aridity Index to Monitor Decreased Soil Water Content and Grass Growth in Grasslands of North-East Asia. Remote Sensing. 2020; 12(21):3556. https://0-doi-org.brum.beds.ac.uk/10.3390/rs12213556

Chicago/Turabian Style

Kimura, Reiji, and Masao Moriyama. 2020. "Use of a Satellite-Based Aridity Index to Monitor Decreased Soil Water Content and Grass Growth in Grasslands of North-East Asia" Remote Sensing 12, no. 21: 3556. https://0-doi-org.brum.beds.ac.uk/10.3390/rs12213556

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop