Next Article in Journal
Deep Learning-Based Framework for Soil Moisture Content Retrieval of Bare Soil from Satellite Data
Next Article in Special Issue
Recent Advances in Crop Disease Detection Using UAV and Deep Learning Techniques
Previous Article in Journal
Solar Photovoltaic Hotspot Inspection Using Unmanned Aerial Vehicle Thermal Images at a Solar Field in South India
Previous Article in Special Issue
Estimation of Anthocyanins in Leaves of Trees with Apple Mosaic Disease Based on Hyperspectral Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Modified Temperature Vegetation Dryness Index (mTVDI) for Agricultural Drought Assessment Based on MODIS Data: A Case Study in Northeast China

1
College of Geo-Exploration Science and Technology, Jilin University, Changchun 130026, China
2
Jilin Institute of GF Remote Sensing Application, Changchun 130012, China
3
Key Laboratory of Water Cycle and Related Land Surface Processes, Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing 100101, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2023, 15(7), 1915; https://doi.org/10.3390/rs15071915
Submission received: 9 March 2023 / Revised: 29 March 2023 / Accepted: 1 April 2023 / Published: 3 April 2023
(This article belongs to the Special Issue Crop Disease Detection Using Remote Sensing Image Analysis II)

Abstract

:
Satellite-based drought indices have been shown to be effective and convenient in detecting drought conditions. The temperature vegetation dryness index (TVDI) is one of the most frequently used drought indices; however, it is not suitable for areas with high fractional vegetation cover (FVC). In this study, a modified temperature vegetation dryness index (mTVDI) was constructed by using the multispectral vegetation dryness index (MVDI) proposed by a PROSAIL simulation and water stress experiments which was based on the theory of the TVDI and utilized MODIS data. Compared with the TVDI, the mTVDI presents a more triangular feature space, as well as obviously increased R2 values for dry and wet edges (from 0.37–0.90 to 0.53–0.91 for dry edges and from 0.00–0.77 to 0.24–0.80 for wet edges). The mTVDI was evaluated using standardized precipitation evapotranspiration indices (SPEIs), precipitation, potential evapotranspiration (PET), and the crop water deficit index (CWDI), and the results confirmed that the mTVDI can better reflect the actual spatial changes, compared to the TVDI, under high FVC, as well as presenting an increased Pearson correlation coefficient (by 0.06–0.10) when compared with SPEIs. Moreover, the good performance of the mTVDI in major drought events indicates its reliability and accuracy for drought monitoring. Overall, the mTVDI is a reliable and accurate satellite-based dryness index suitable for high FVC conditions.

1. Introduction

Drought is one of the most serious natural disasters, causing serious water and food security problems all over the world [1,2,3,4]. The mechanisms of drought are complex and poorly understood. From different perspectives, droughts can be divided into meteorological, agricultural, hydrological, and socioeconomic droughts [1,2,5]. Given the severe impacts of drought on agriculture, it is urgent and important to develop effective methods for agricultural drought assessment.
In recent years, remote sensing has been shown to be capable of providing fast, economical, non-destructive, and spatio-temporal measurements of agricultural drought. A number of drought indices have been proposed to monitor agricultural drought based on visible (VIS), near-infrared (NIR), shortwave infrared (SWIR), and thermal infrared (TIR) remote sensing data. The drought indices can be divided into vegetation, temperature, and temperature–vegetation indices [6,7].
A vegetation index can detect drought conditions by combining spectral information from various regions that are sensitive to the physiology and biophysical characteristics of vegetation. The normalized difference vegetation index (NDVI) is the earliest and most widely used vegetation index, which uses non-linear stretching to enhance the reflectance contrast in the NIR and red spectral regions [8]. When the fractional vegetation cover (FVC) is higher than 80%, the increase in the NDVI is delayed, showing a saturation state, which leads to decreased sensitivity of vegetation detection under high FVC [9]. To overcome the saturation phenomenon of the NDVI, further vegetation indices have been established, such as the enhanced vegetation index (EVI) [9,10,11]. These indices are often referred to as greenness indices because they reflect the greenness of the vegetation [12,13].
In addition to the change in vegetation index, stomatal closure on leaves under water stress can lead to an increase in canopy temperature, which occurs earlier and provides a more sensitive indicator than a change in the vegetation index. Based on this physiological characteristic of crops, many indices based on land surface temperature (LST) have been proposed for drought monitoring [14]. These indices mainly include the temperature condition index (TCI) [15,16] and normalized difference temperature index (NDTI) [17]. Temperature indices have better performance than vegetation indices in accurately monitoring soil moisture due to evapotranspiration and leaf temperature [18]; however, a temperature index cannot distinguish the differences in drought resistance among different vegetation types, leading to certain inaccuracies in drought assessment [19].
In recent years, various approaches to combine vegetation and temperature indices for the detection of drought have been explored [16,18,20,21]. The most widely used temperature–vegetation indices include the vegetation supply water index (VSWI), vegetation health index (VHI), and temperature vegetation dryness index (TVDI). The TVDI is constructed according to the NDVI–LST triangle space, relying on the relationship between LST and the NDVI (which are usually negatively correlated) [22,23,24]. Compared with the individual vegetation and temperature indices, the TVDI is more reliable for drought monitoring and impact assessment [25,26,27]. However, the TVDI has limitations regarding drought monitoring in high FVC areas, due to the saturation of the NDVI [28,29]. Therefore, it is necessary to modify the TVDI in order to improve the accuracy and reliability of agricultural drought monitoring.
The aim of this study is to improve the accuracy of drought monitoring in high FVC areas. We present a modified drought index, termed the modified temperature vegetation dryness index (mTVDI), based on the multispectral vegetation dryness index (MVDI) and LST.

2. Materials and Methods

2.1. Study Area

The study area covers an area of 1,245,202 km2 in Northeast China (38°40′—53°20′N, 111°37′—135°5′E). The terrain of the study area is mainly mountainous and plain, among which the Northeast plain is rich in black soil resources and fertile soil, making it an important grain base in China. The main climate type of the region is temperate monsoon climate, with an annual average temperature of 5.4 °C. The annual average precipitation is about 600 mm, which gradually decreases from east to west. The eastern part is humid, while the western part is semi-humid. The main local crop types in drylands are maize and soybean, and the growing period is from May to September. In this paper, the areas of corn and soybean were masked. In the extracted areas, the FVC was low in May due to seeding, and the FVC gradually increased with the growth of crops from June to September.
The most frequent natural disaster in Northeast China throughout the 20th century has been drought, which causes severe water shortages, environmental degradation, and negative socioeconomic effects. Extreme drought occurrences have become increasingly frequent, leading to significant agricultural losses [30,31].

2.2. Data

The Moderate Resolution Imaging Spectroradiometer (MODIS) data from May to September, 2007–2017, was chosen from the Google Earth Engine platform (https://code.earthengine.google.com, accessed on 1 June 2022). The MVDI and NDVI data were obtained from the MOD09A1 product corresponding to the Terra platform, with a spatial resolution of 500 m and a temporal resolution of 8 days. The LST data were generated, with a spatial resolution of 1 km and a temporal resolution of 8 days, using the MYD11A2 product corresponding to the Aqua platform. To match the spatial resolution of LST, the resolution of the MOD09A1 product was resampled to 1 km by cubic convolution interpolation. In general, the temporal resolution of the mTVDI was 8 days. When compared with other indices, such as the SPEI, the mTVDI was averaged monthly, in order to match the temporal scale.
The crop classification data were obtained from the 10 m crop type maps of the major crops (i.e., maize, soybean, and rice) in Northeast China in 2017 [32]. The dataset used the Google Earth Engine platform, Sentinel-2 remote sensing image, and random forest algorithm to classify major crops (i.e., maize, soybean, and rice) in Northeast China from 2017 to 2019. The resultant maps have overall accuracies spanning from 0.81 to 0.86.
The meteorological data were selected from the National Meteorological Centre of China Meteorological Administration (NMCCMA; https://data.cma.cn, accessed on 1 July 2022). In particular, daily precipitation, temperature, and sunshine duration were downloaded from 115 meteorological stations (Figure 1c) from May to September, 1978–2017.
The documented records of major drought events in Northeast China were obtained from the China Drought Dataset (CDD), which is released by China Engineering Science and Technology Knowledge Center (https://www.ckcest.cn, accessed on 15 July 2022).

2.3. Methodology

2.3.1. The Technical Route

The flowchart of this study is shown in Figure 2. We used the PROSAIL model to simulate the canopy spectral reflectance response of drought crops. Through water stress experiments of maize and soybean, it was shown that the spectral characteristic bands can detect the change of crops in response to drought. Using characteristic bands based on simulation and experiments, the MVDI was established. Then, we constructed the mTVDI by constructing the triangular feature space of the MVDI and LST based on MODIS data. We used the standardized precipitation evapotranspiration index (SPEI), the TVDI, precipitation, potential evapotranspiration (PET), the crop water eeficit index (CWDI), and the major drought events to evaluate the mTVDI. Finally, we used the mTVDI for monitoring drought in Northeast China.

2.3.2. Canopy Reflectance Simulation by PROSAIL

Model simulation can help to quantify the contributions of the canopy’s biophysical and biochemical properties to canopy reflectance. Since about thirty years ago, researchers have studied plant canopy spectral and directional reflectance in the solar domain using the PROSAIL radiative transfer model, which takes into account all biochemical parameters of leaves, vegetation canopy structure, and bidirectional scattering characteristics. PROSAIL-5B [33], which combines the 4SAIL and PROSPECT-5B [34,35] models, was employed in this study. The former defines how light is transmitted on the leaf’s surface and within, and is appropriate for both thick and sparse leaf vegetation. The latter model, known as the scattering by arbitrarily inclined leaves (SAIL) model, has been frequently used to investigate the vegetation canopy. The plant canopy in the model is thought to be a mixed medium, made up of isotropic leaves that are uniformly horizontal, endlessly long, and arbitrarily oriented. The transmission of light through the vegetation canopy is explained by radiative transfer theory, and canopy reflectance can be calculated.
For this study, 168 canopy spectral reflectance curves were simulated by successively changing the PROSAIL input parameters—equivalent water thickness (EWT), chlorophyll a+b content (Cab), and the leaf area index (LAI)—which are the most sensitive to water stress. The values and steps of the parameters are given in Table 1.
A sensitivity analysis was conducted to identify the characteristic bands of the spectral reflectance, which are strongly responsive to parameters. The sensitivity calculation formula is as follows:
S = i n 1 ( R i R i + 1 / R i )
where R i is the spectral reflectance of the corresponding parameter.
Figure 3 shows the canopy spectral reflectance and sensitivity analysis, as simulated by the PROSAIL model.
The EWT can be used for the characterization of crop water content and, further, to characterize crop drought [36]. It can be seen, in Figure 3a,b, that the range sensitive to crop water is in the NIR and SWIR, with the sensitive bands being 970, 1200, 1450, 1895, and 1980 nm. The function of the Cab is to convert light energy into chemical energy that plants can use. As can be seen in Figure 3c,d, the sensitive bands in this regard were 610 and 670 nm. It can be seen, in Figure 3e,f, that the spectral reflectance increased in the VIS and SWIR, decreased in the NIR, and changed changes in the LAI. The sensitive bands were 670, 780, 1450, and 1930 nm. This can be explained by interior scattering within the leaves at the air, cell, and water interfaces [37]. The PROSAIL simulation indicated that the spectral responses were generated by the different characteristic changes of crops subjected to water stress.

2.3.3. Water Stress Experiments

Experiments were conducted in the potted proving ground of Jilin University, Jilin, China (43°57′N, 125°15′E) during May–September of 2020–2021, using maize (Xianyu 335) and soybean (Jipin 609). To ensure stress environments, the maize and soybean were subjected to four water stress levels by adjusting the field capacity (FC) in the following ways: no water stress (N: 70–80% FC), light water stress (L: 50–60% FC), moderate water stress (M: 40–50% FC), and severe water stress (S: 30–40% FC). Nine replicates were set for each treatment, arranged in parallel (3 × 3). The samples were planted in flower pots, with one maize plant in each pot and four soybean plants in each pot. The quantity and quality of soil in the flowerpots were basically the same. The pots were 30 cm in diameter, 40 cm in height, and 35 cm in soil thickness. The soil was black in texture and the soil moisture retention was 23.8%. The soil moisture content of maize and soybean was controlled by an awning and irrigation.
The hyperspectral data were collected from the crop canopy by using an ASD Field Spec 3 portable spectrometer. This spectrometer had the following characteristics: a spectral region of 350–2500 nm, a spectral resolution of 3 nm in the VIS, 8.5 nm in the NIR, and 6.5 nm in the SWIR, a spectra sampling interval of 1.4 nm (350–1000 nm) and 2 nm (1000–2500 nm), and a 1 nm spectra resampling interval [38]. Reflectance measurements were made on sunny and windless days at 10:00–14:00 Beijing time. The field of view (FOV) was 25° and the distance between the optical head of the spectrometer and the top of the crop was kept at 1 m for all observations. During the measurements, standard whiteboard calibration was performed every 30 min. The spectral characteristics of the crops were measured during the critical period of crop growth (see Table 2).
The spectral reflectance of maize and soybean under water stress have similar characteristics in different growth stages. Figure 4 shows examples of the spectral reflectance of maize and soybean under different water stress levels. With an increase in water stress, the spectral reflectance of the maize and soybean canopy increased in the VIS and SWIR, while decreasing in the NIR. These spectral changes were related to the physiological and biochemical characteristics of crops under water stress. For example, with an increase in water stress, the decrease in chlorophyll content led to the weakening of spectral absorption which, in turn, led to the reflectance increase in VIS, while the decrease in the LAI led to the reflectance decrease in NIR, and the decrease in water content was reflected by the increased reflectance in the SWIR [6,39].
Figure 5 shows the sensitivity analysis of the canopy spectra of soybean and maize under different water stress levels. It can be seen that the sensitive bands of water stress for soybean were 480, 670, 780, 1200, 1450, and 2020 nm, while the sensitive bands of water stress for maize were 430, 670, 780, 1200, and 1450 nm. This is the result of a combination of crop water, structure, and chlorophyll. As can be seen in Figure 4 and Figure 5, the results of the water stress experiments for soybean and maize verified the results of the PROSAIL simulation, confirming that spectral information can be used to accurately detect the state of water stress in crops.

2.3.4. Multispectral Vegetation Dryness Index (MVDI)

Based on the results of the PROSAIL simulation and water stress experiments, we summarized the spectral response characteristic bands caused by different physiological and biochemical changes. The 430, 480, 610, and 670 nm wavelengths were chosen as indicators for chlorophyll content; the 970, 1200, 1450, 1895, 1980, and 2020 nm wavelengths were chosen as indicators for water content; and the 670, 780, 1450, and 1930 nm wavelengths were chosen as indicators for variability in crop structure.
Crops respond physiologically and biochemically to water stress in a variety of ways, such as by stunting growth, coiling or necrosis of photosynthetically active sections, wilting or, in extreme circumstances, a reduction in leaf area due to severe defoliation. In order to monitor crop drought, some bands may be helpful in spotting changes in crop water, structure, and chlorophyll. Determining an index based on a mixture of bands from spectral regions responsive to various crop attributes altered by water stress would, therefore, be valuable. To formulate a new spectral index, we combined the characteristic bands from different spectral responses. By combining different bands of chlorophyll content and crop structure, we found that the reflectance difference between 780 nm and 670 nm was most sensitive to water stress. On this basis, the difference between 780 nm and 670 nm and the water content bands were combined. Finally, the hyperspectral vegetation dryness index (HVDI), which is most relevant to water stress, was selected, as shown in Equation (2):
HVDI = ( R 780 R 670 ) / R 1200
where R x is the reflectance at the corresponding wavelength.
Considering the application of the HVDI, we transform it into a multispectral form that can be applied to the MODIS. The MODIS has been widely used for drought monitoring, due to its advantages of large imaging width, wide spectral range, and high updating frequency [1,9,40,41]. Therefore, MODIS imagery was selected in this study, such that the HVDI in Equation (2), by changing the narrow band channel to the wide band channel, as shown in Equation (3):
MVDI = ( R NIR 1 R Re d ) / R NIR 2
where MVDI is the multispectral HVDI based on MODIS, R Red denotes the MODIS red band (620–670 nm), and R NIR 1 and R NIR 2 are the MODIS near-infrared bands ( NIR 1 : 841–876 nm, NIR 2 : 1230–1250 nm). Although the center of the MODIS channel is off-HVDI, the features in both channels are comparable, as we discuss in Section 4.1.
In general, the value of the MVDI is greater than 0; the smaller the MVDI, the stronger the water stress.

2.3.5. Modified Temperature Vegetation Dryness Index (mTVDI)

On the basis of the TVDI theory [23,42,43], the mTVDI was developed by combining the MVDI and LST to construct a triangle feature space (Figure 6). In the MVDI–LST feature space, point P represents dry bare soil and point Q represents wet bare soil. With an increase in FVC, the LST decreases. Point R represents where the FVC is highest and soil moisture is sufficient, such that evaporation and transpiration are maximal. Therefore, the dry edge indicates low soil water availability and small surface evapotranspiration, while the wet edge denotes adequate soil moisture and normal vegetation growth.
The MVDI–LST feature space has a number of contours related to soil moisture (Figure 6). Based on this correlation, the mTVDI is defined as follows:
mTVDI = ( LST LST min ) / ( LST max LST min )
LST min = a 1 + b 1 × MVDI
LST max = a 2 + b 2 × MVDI
where LST is the land surface temperature; LST min and LST max are the lowest and highest LST corresponding to a certain MVDI, respectively, which can be determined by extracting the wet edge using Equation (5) and dry edge using Equation (6) through linear regression analysis; and a 1 , b 1 , a 2 , and b 2 are undetermined coefficients. The mTVDI value ranges from 0 to 1, with values closer to 1 indicating more severe drought.
In this study, LST was extracted by using 0.01 as the interval of the MVDI, and a 95% confidence interval was selected to leave out possible abnormal LST. Then, a series of maximum and minimum LSTs corresponding to the MVDI were calculated. After that, the dry and wet edges of the mTVDI were extracted using the linear regression method. The performance evaluation of the mTVDI is discussed in the Results and Discussion sections.

2.3.6. Other Drought Indices

In addition, nine vegetation indices that have been evaluated by other researchers using various vegetation data and have shown stable performance were calculated in order to comparatively assess the effectiveness of the MVDI. According to the number of bands and the selected algorithm, the spectral indices were classified into three types: structural, chlorophyll, and water indices (Table 3). We compared the MVDI with the nine conventional vegetation indices in Section 3.1.
In addition, the TVDI, SPEI, and crop water deficit index (CWDI) were calculated to assess the effectiveness of the mTVDI. The TVDI is a method for measuring soil moisture based on the triangle feature space of NDVI–LST, which is common and reliable for use in drought monitoring [22].
The SPEI is a drought index calculated according to precipitation and temperature, which can be calculated at various time scales and has been widely used in drought research [51,52]. In this study, the time scales of one month (SPEI1), three months (SPEI3), and six months (SPEI6) were calculated for comparison with the mTVDI.
The CWDI is an indicator used to characterize agricultural drought through the relationship between precipitation and PET [53,54]. The CWDI was used to verify the drought monitoring results of the mTVDI.

2.3.7. Statistical Metrics

In order to quantitatively evaluate the ability of the MVDI and the nine other vegetation indices to identify the degree of water stress, the normalized mean distance was introduced. According to the Fisher criterion, the relative distance between categories can be used as a measure of category separability; the greater the distance, the greater the separability [55]. The normalized mean distance principle formula is as follows [56]:
d norm = μ m μ n σ m + σ n
where μ m and μ n are the respective mean values of the vegetation indices under different water stress levels and σ m and σ n are the respective variances of the vegetation indices under different water stress levels.
To quantitatively evaluate the performance of the mTVDI with that of the TVDI and SPEI, the Pearson correlation coefficient was adopted as a statistical metric. Its calculation formula is as follows:
P X , Y = C o v ( X , Y ) V a r ( X ) V a r ( Y )
where X and Y are the variables for comparison, C o v is the covariance of X and Y, and V a r is the variance.
To verify the match between the HVDI bands and MODIS channels, a residual analysis was adopted as an evaluation metric for deviation. Its calculation formula is as follows:
ε = y y ^
where ε is the residual between the values of the HVDI and MVDI, and y and y ^ are the values of the HVDI and MVDI.

3. Results

3.1. Comparison of the MVDI with Other Vegetation Indices

The MVDI and nine other vegetation indices were selected to identify crop water stress levels using the hyperspectral data measured in the water stress experiments (Figure 7). As maize and soybean have similar spectral characteristics (Figure 4 and Figure 5), they are discussed together. As can be seen in Figure 7, the NDVI, EVI, PRI, SIPI, WI, NDII, and NDWI can identify water-stressed crops in the early, middle, and late stages, but could not accurately identify the water stress degree in the mature stage. The MCARI values of control and water-stressed crops presented no constant rule, except in the early stage, which cannot allow for accurate identification of water-stressed crops. The MVDI and RDVI could identify water-stressed crops throughout the whole growth period.
Moreover, to quantitatively evaluate the ability of the MVDI and the other nine vegetation indices to identify the degree of water stress, the normalized mean distance was calculated to judge the index recognition ability, as shown in Table 4. In the early, late, and mature stages, the mean distance values for the MVDI were 4.61, 2.68, and 1.46, respectively, which were all larger than those of the other indices. In the middle stage, the mean distance of the MVDI was 1.31, which was second only to that of the PRI (1.53). However, as the PRI could not distinguish the water stress levels in the mature stage; it could be considered that the MVDI presented the best comprehensive performance.
In conclusion, the MVDI had a stronger ability to identify the physiological and biochemical changes of crops than the NDVI, RDVI, EVI, PRI, SIPI, MCARI, WI, NDII, and NDWI, and it had strong stability throughout the whole growth period of the crops.

3.2. Evaluating the Feature Space of the MVDI–LST

Figure 8 shows the MVDI–LST space of the mTVDI and the NDVI–LST space of the TVDI. The results show that the shape of the MVDI–LST space was more triangular than the NDVI–LST space, especially from June to September when the FVC is high. Moreover, the linear fitting of dry and wet edges in the MVDI–LST space was better than that in the NDVI–LST. The dry edge R2 values for the NDVI–LST space were 0.90, 0.78, 0.72, 0.45, and 0.37, while those for the MVDI–LST space were 0.91, 0.88, 0.90, 0.82, and 0.53. The wet edge R2 values for the NDVI–LST space were 0.77, 0.01, 0.00, 0.00, and 0.49, while those for the MVDI–LST space were 0.80, 0.26, 0.52, 0.24, and 0.35.

3.3. Applying the mTVDI in Drought Monitoring in Northeast China

Given the excellent results of the MVDI–LST feature space for the proposed mTVDI, MODIS data were used to calculate the mTVDI in order to assess the spatio-temporal distribution of droughts in Northeast China in 2017 (Figure 9). Figure 9 shows the evolution process of drought over the crop growth period. The western part of Northeast China experienced widespread drought in May, which may have been caused by insufficient precipitation in the early period. Then, the drought spread in June, due to an increase in water demand as crops began to grow rapidly. As temperatures increased, the drought covered almost the whole of Northeast China in July, which may have been due to the increasing evapotranspiration and water demand. From August to September, the drought finally eased in Northeast China, except for the persistent drought in the southwest part.
Moreover, we used precipitation, PET, and the CWDI to explain the spatio-temporal distribution of droughts in Northeast China monitored by the mTVDI (Figure 10). Due to the spatial difference in drought occurrence, we divided Northeast China into four regions (Heilongjiang, Jilin, Liaoning, and Neimenggu) for comparison and analysis. The mTVDI showed an aggravating trend of drought from May to July, which was caused by the continuous low precipitation and the increase in crop PET due to an increase in temperature. From August to September, the mTVDI indicated that the drought eased significantly, which was caused by the significant increase in precipitation in August and the decrease in crop PET due to the lower temperatures, consistent with the results of the CWDI. Among the four regions in Northeast China, the mTVDI showed that Neimenggu suffered the worst drought, which was caused by the low precipitation and high PET during the whole crop-growing period. The next worst droughts were in Liaoning and Jilin, while the drought was lightest in Heilongjiang.
In addition, the distribution of drought in the four regions may be related to crop types. Studies have shown that soybean is more drought-resistant than maize; this is because, in the growth process of soybean, its small leaf area leads to lower plant transpiration and, thus, less water is needed [57,58,59]. The north and east parts of Heilongjiang are soybean-growing regions, as shown in Figure 1b, and the drought tolerance of soybean plays a regulating role in drought. Jilin and Liaoning are mostly maize-growing regions and, thus, are less resistant to drought. In Neimenggu, although the soybean-growing region is wide, the drought is more severe, which may be mainly caused by the intolerance of soybean with respect to high temperatures [59,60].

3.4. Comparison of the mTVDI with Other Drought Indices

To evaluate the performance of the mTVDI and its consistency with other drought indices, the mTVDI was compared with the TVDI and various SPEIs, and spatio-temporal distribution maps of drought were created, as shown in Figure 9.
In general, these indices behaved similarly, in terms of the spatial trends of drought. However, there were some differences in the duration, distribution, and degree of drought. For example, both the mTVDI and the SPEIs indicated that the drought in Northeast China gradually intensified from May to July and gradually eased from August to September; however, the TVDI results showed that the drought was most severe in May, relieved from June to July, and intensified from August to September. This may be caused by NDVI saturation in areas with high FVC (usually an LAI greater than 3 [61]). In addition, the results of the mTVDI, SPEI3, and SPEI6 indicated that the most severe drought occurred in the southwest of Northeast China in July, while SPEI1 and the TVDI presented no significant drought in the south. This may be due to the bias caused by SPEI1, taking into account only one month’s precipitation and PET.
In order to compare the performance of the mTVDI and TVDI quantitatively, Pearson correlation coefficients were calculated between the mTVDI, TVDI, and SPEIs, as shown in Table 5. It can be seen that the mTVDI had a stronger correlation with SPEIs than the TVDI. For example, from May to September, the average Pearson correlation coefficients of the mTVDI and SPEIs were −0.75, −0.63, −0.70, −0.64, and −0.75 (p < 0.01), which were all greater than the correlation coefficients with the TVDI. The results for the mTVDI, TVDI, and SPEIs showed that the mTVDI has a more accurate spatial distribution than the TVDI under high FVC.

3.5. Performance of the mTVDI in the Major Drought Events

To further verify the validity and applicability of the mTVDI, we looked for records of drought events that had occurred in Northeast China (Table 6). Due to the small number of drought records, we only obtained drought event data for the four considered regions (Heilongjiang, Jilin, Liaoning, and Neimenggu) during 2007–2010. The mTVDI values of the month when the drought events occurred were calculated, with mTVDI values closer to 1 indicating more drought. In addition, SPEI3 was calculated to assess the accuracy of the mTVDI due to its reliability in agricultural drought monitoring [3,42,62,63,64]. For ease of comparison, the SPEI results were normalized into the range from 0 to 1, with values closer to 0 indicating more drought. The results are shown in Figure 11, which demonstrates the drought records using red marks.
The results showed that the drought events identified by the mTVDI and SPEI3 were basically consistent with those recorded in the CDD. For example, in the four drought events in Liaoning, the mTVDI showed high values and SPEI3 showed low values, both of which represent drought. The performance of the mTVDI regarding the major drought events indicated its reliability and accuracy in drought monitoring.

4. Discussion

4.1. Matching between HVDI Bands and MODIS Channels

In this study, the MVDI was proposed for application to MODIS on the basis of the HVDI. It is based on the fact that the HVDI bands and the MVDI channels have the same meanings: R 670 of the HVDI and R Red of the MVDI represent the strong reflectance absorption in the red spectral region, R 780 of the HVDI and R NIR 1 of the MVDI represent the high reflectance platform in the NIR, and R 1200 of the HVDI and R NIR 2 of the MVDI represent the water absorption peaks in the NIR [10,36]. However, the bands of the HVDI do not exactly correspond to the central wavelength of the channels of MODIS; thus, it is necessary to verify the match between the HVDI bands and MODIS channels. The HVDI and MVDI were calculated using the hyperspectral data measured in the water stress experiments, and the MVDI channels were obtained by convolution of the spectral response function of MODIS. The results of the HVDI and MVDI, and the corresponding residual analysis, are shown in Figure 12.
As can be seen in Figure 12a–d, there is a good consistency between the HVDI and MVDI at different stages of crop growth. The residual plots provided a reasonable pattern in which residuals are distributed randomly and unpredictably (Figure 12e–h) [65]. The results show that the HVDI bands match well with MODIS channels; thus, the MVDI can express the spectral response characteristic represented by the HVDI.

4.2. Anti-Saturation of the mTVDI in High FVC

In the MVDI–LST feature space of the mTVDI, the MVDI reflects not only the red edge but also the water absorption characteristics of crops [36,61]. According to previous research, there are four significant absorption peaks—at about 970, 1200, 1450, and 1950 nm—which can be used for the diagnosis of the water status in plants [66,67]. However, only the water absorption bands at 1200 and 1450 nm were detected in the water stress experiments of maize and soybean (Figure 5), which may be caused by the difference in canopy geometry and the influence of the soil background [68,69]. In addition, although absorption at 1200 nm is not as strong as that at 1450 nm, it is less susceptible to atmospheric influence in satellite applications [49]. Therefore, the addition of 1200 nm to the MVDI can better reflect the physiological and biochemical state of vegetation than the NDVI. The MVDI incorporates water information based on vegetation information, which makes the mTVDI perform good for drought monitoring in high FVC.
On the other hand, the NDVI uses non-linear stretching to enhance the reflectance contrast in the NIR and red spectral regions, which causes saturation in high FVC [9,70]. Under high FVC conditions, the increase in the NDVI slows down with the increase in FVC. The form of the MVDI formula also avoids the saturation effect for the mTVDI in high FVC conditions.

4.3. Limitation and Future Research

In this study, we verified the ability of the mTVDI in drought monitoring of major crops in Northeast China and achieved good results. However, the mTVDI still has some uncertainties, such as those relating to the adaptability of other vegetation types and other areas. Therefore, the mTVDI should be further validated by utilizing multiple vegetation types and wider areas for drought monitoring.
More recently, hyperspectral remote sensing technology has played an important role in quantitative remote sensing because of its high spectral resolution. In terms of monitoring and application of drought, it also has advantages that traditional multispectral remote sensing does not have, which can capture and discriminate subtle drought response signals and monitor vegetation drought [36,55,71]. In future research, a hyperspectral-modified TVDI should be constructed by the triangular feature space of the HVDI and LST, and further verified and applied to hyperspectral images.

5. Conclusions

In this study, in order to improve the applicability of drought monitoring based on the TVDI under high FVC, we constructed the mTVDI from the MVDI based on simulation and experimental data. By comparing the feature spaces of the mTVDI and TVDI, the results indicated that the mTVDI presented a more triangular feature space than the TVDI, and had obviously increased the R2 values for the dry and wet edges (from 0.37–0.90 to 0.53–0.91 for dry edges and from 0.00–0.77 to 0.24–0.80 for wet edges). The Pearson correlation coefficient of the mTVDI was 0.06–0.10 higher than that of the TVDI; thus, it can better reflect the actual spatial changes, indicating its practical value in high FVC areas. As a result, the more trustworthy results generated using the mTVDI can aid in better field management, irrigation planning, and agricultural drought risk assessment.

Author Contributions

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

Funding

This work was supported by the National Key Research and Development Program of China (grant No. 2020YFA0714103) and the Scientific and Technological Development Scheme of Jilin Province (grant No. 20210201138GX).

Data Availability Statement

Data sharing is not application to this article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Du, L.; Tian, Q.; Yu, T.; Meng, Q.; Jancso, T.; Udvardy, P.; Huang, Y. A comprehensive drought monitoring method integrating MODIS and TRMM data. Int. J. Appl. Earth Obs. Geoinf. 2013, 23, 245–253. [Google Scholar] [CrossRef]
  2. Chang, K.Y.; Xu, L.; Starr, G.; Kyaw Tha Paw, U. A drought indice reflecting ecosystem responses to water availability: The Normalized Ecosystem Drought Index. Agric. For. Meteorol. 2018, 250–251, 102–117. [Google Scholar] [CrossRef]
  3. Kogan, F.N. Global drought watch from space. Bull. Am. Meteorol. Soc. 1997, 78, 621–636. [Google Scholar] [CrossRef]
  4. Yang, J.; Chang, J.; Wang, Y.; Li, Y.; Hu, H.; Chen, Y.; Huang, Q.; Yao, J. Comprehensive drought characteristics analysis based on a nonlinear multivariate drought index. J. Hydrol. 2018, 557, 651–667. [Google Scholar] [CrossRef]
  5. Hao, C.; Zhang, J.; Yao, F. Combination of multi-sensor remote sensing data for drought monitoring over Southwest China. Int. J. Appl. Earth Obs. Geoinf. 2015, 35, 270–283. [Google Scholar] [CrossRef]
  6. Gerhards, M.; Schleft, M.; Mallick, K.; Udelhoven, T. Challenges and future perspectives of multi-/hyperspectral thermal infrared remote sensing for crop water-stress detection: A review. Remote Sens. 2019, 11, 1240. [Google Scholar] [CrossRef] [Green Version]
  7. Chaerle, L.; Van Der Straeten, D. Imaging techniques and the early detection of plant stress. Trends Plant Sci. 2000, 5, 495–501. [Google Scholar] [CrossRef]
  8. Tucker, C.J.; Newcomb, W.W.; Los, S.O.; Prince, S.D. Mean and inter-year variation of growing season normalized difference vegetation index for the Sahel 1981–1989. Int. J. Remote Sens. 1991, 12, 1133–1135. [Google Scholar] [CrossRef]
  9. Qin, Q.; Wu, Z.; Zhang, T.; Sagan, V.; Zhang, Z.; Zhang, Y.; Zhang, C.; Ren, H.; Sun, Y.; Xu, W.; et al. Optical and Thermal Remote Sensing for Monitoring Agricultural Drought. Remote Sens. 2021, 13, 5092. [Google Scholar] [CrossRef]
  10. Huete, A.; Didan, K.; Miura, T.; Rodriguez, E.P.; Gao, X.; Ferreira, L.G. Overview of the radiometric and biophysical performance of the MODIS vegetation indices. Remote Sens. Environ. 2002, 83, 195–213. [Google Scholar] [CrossRef]
  11. Saleska, S.R.; Didan, K.; Huete, A.R.; Da Rocha, H.R. Amazon forests green-up during 2005 drought. Science 2007, 318, 612. [Google Scholar] [CrossRef] [Green Version]
  12. Jackson, T.J.; Chen, D.; Cosh, M.; Li, F.; Anderson, M.; Walthall, C.; Hunt, E.R. Vegetation water content mapping using Landsat data derived normalized difference water index for corn and soybeans. Remote Sens. Environ. 2004, 92, 475–482. [Google Scholar] [CrossRef]
  13. AghaKouchak, A.; Farahmand, A.; Melton, F.S.; Teixeira, J.; Anderson, M.C.; Wardlow, B.D.; Hain, C.R. Remote sensing of drought: Progress, challenges and opportunities. Rev. Geophys. 2015, 53, 452–480. [Google Scholar] [CrossRef] [Green Version]
  14. Gutman, G.G. Towards monitoring droughts from space. J. Clim. 1990, 3, 282–295. [Google Scholar] [CrossRef]
  15. Bhuiyan, C.; Singh, R.P.; Kogan, F.N. Monitoring drought dynamics in the Aravalli region (India) using different indices based on ground and remote sensing data. Int. J. Appl. Earth Obs. Geoinf. 2006, 8, 289–302. [Google Scholar] [CrossRef]
  16. Jain, S.K.; Keshri, R.; Goswami, A.; Sarkar, A.; Chaudhry, A. Identification of drought–vulnerable areas using NOAA AVHRR data. Int. J. Remote Sens. 2009, 30, 2653–2668. [Google Scholar] [CrossRef]
  17. McVicar, T.R.; Jupp, D.L. The current and potential operational uses of remote sensing to aid decisions on drought exceptional circumstances in Australia: A review. Agric. Syst. 1998, 57, 399–468. [Google Scholar] [CrossRef]
  18. Goetz, S.J. Multi-sensor analysis of NDVI, surface temperature and biophysical variables at a mixed grassland site. Int. J. Remote Sens. 1997, 18, 71–94. [Google Scholar] [CrossRef]
  19. Moran, M.S. Thermal infrared measurement as an indicator of planet ecosystem health. Therm. Remote Sens. Land Surf. Processes 2004, 1, 257–282. [Google Scholar]
  20. Price, J.C. On the analysis of thermal infrared imagery: The limited utility of apparent thermal inertia. Remote Sens. Environ. 1985, 18, 59–73. [Google Scholar] [CrossRef]
  21. Kogan, F.N. Application of vegetation index and brightness temperature for drought detection. Adv. Space Res. 1995, 15, 91–100. [Google Scholar] [CrossRef]
  22. Sandholt, I.; Rasmussen, K.; Andersen, J. A simple interpretation of the surface temperature/vegetation index space for assessment of surface moisture status. Remote Sens. Environ. 2002, 79, 213–224. [Google Scholar] [CrossRef]
  23. Liu, L.; Liao, J.; Chen, X.; Zhou, G.; Su, Y.; Xiang, Z.; Wang, Z.; Liu, X.; Li, Y.; Wu, J.; et al. The Microwave Temperature Vegetation Dryness Index (MTVDI) based on AMSR-E brightness temperatures for long-term drought assessment across China (2003–2010). Remote Sens. Environ. 2017, 199, 302–320. [Google Scholar] [CrossRef]
  24. Liu, Y.; Yue, H. The Temperature Vegetation Dryness Index (TVDI) Based on Bi-Parabolic NDVI-Ts Space and Gradient-Based Structural Similarity (GSSIM) for Long-Term Drought Assessment Across Shaanxi Province, China (2000–2016). Remote Sens. 2018, 10, 959. [Google Scholar] [CrossRef] [Green Version]
  25. Zhao, S.; Cong, D.; He, K.; Yang, H.; Qin, Z. Spatial-temporal variation of drought in China from 1982 to 2010 based on a modified Temperature Vegetation Dryness Index (mTVDI). Sci. Rep. 2017, 7, 1–12. [Google Scholar] [CrossRef] [Green Version]
  26. Qi, S.; Wang, C.; Niu, Z. Evaluating soil moisture status in China using the temperature/vegetation dryness index (TVDI). J. Remote Sens. 2003, 7, 420–427. [Google Scholar]
  27. Sun, L.; Sun, R.; Li, X.; Liang, S.; Zhang, R. Monitoring surface soil moisture status based on remotely sensed surface temperature and vegetation index information. Agric. For. Meteorol. 2012, 166, 175–187. [Google Scholar] [CrossRef]
  28. West, H.; Quinn, N.; Horswell, M. Remote sensing for drought monitoring & impact assessment: Progress, past challenges and future opportunities. Remote Sens. Environ. 2019, 232, 111291. [Google Scholar]
  29. Sha, S.; Guo, N.; Li, Y.; Han, T.; Zhao, Y. Introduction of Application of Temperature Vegetation Dryness Index in China. J. Arid Meteorol. 2014, 32, 128–134. [Google Scholar]
  30. Cong, D.M.; Zhao, S.H.; Chen, C.; Duan, Z. Characterization of droughts during 2001–2014 based on remote sensing: A case study of Northeast China. Ecol. Inform. 2017, 39, 56–67. [Google Scholar] [CrossRef]
  31. Guo, E.L.; Liu, X.P.; Zhang, J.Q.; Wang, Y.F.; Wang, C.L.; Wang, R.; Li, D.J. Assessing spatiotemporal variation of drought and its impact on maize yield in Northeast China. J. Hydrol. 2017, 553, 231–247. [Google Scholar] [CrossRef]
  32. You, N.; Dong, J.; Huang, J.; Du, G.; Zhang, G.; He, Y.; Yang, T.; Di, Y.; Xiao, X. The 10-m crop type maps in Northeast China during 2017–2019. Sci. Data 2021, 8, 41. [Google Scholar] [CrossRef] [PubMed]
  33. Feret, J.B.; Gitelson, A.A.; Noble, S.D.; Jacquemoud, S. PROSPECT-D: Towards modeling leaf optical properties through a complete lifecycle. Remote Sens. Environ. 2017, 193, 204–215. [Google Scholar] [CrossRef] [Green Version]
  34. Feret, J.B.; Franois, C.; Asner, G.P.; Gitelson, A.A.; Martin, R.E.; Bidel, L.P.; Ustin, S.L.; le Maire, G.; Jacquemoud, S.p. PROSPECT-4 and 5: Advances in the leaf optical properties model separating photosynthetic pigments. Remote Sens. Environ. 2008, 112, 3030–3043. [Google Scholar] [CrossRef]
  35. Jacquemoud, S.; Verhoef, W.; Baret, F.; Bacour, C.; Zarco-Tejada, P.J.; Asner, G.P.; Francois, C.; Ustin, S.L. PROSPECT + SAIL models: A review of use for vegetation characterization. Remote Sens. Environ. 2009, 113, 556–566. [Google Scholar] [CrossRef]
  36. Vinod, P.; Ratnadeep, R. Hyperspectal remote sensing for agriculture: A review. Int. J. Comput. Appl. 2017, 172, 30–34. [Google Scholar] [CrossRef]
  37. Hunt, J.; Ramond, E.; Rock, B.N. Detection in changes in leaf water content using near and mid-infrared reflectance. Remote Sens. Environ. 1989, 30, 45–54. [Google Scholar]
  38. Xu, Z.Y.; Chen, S.B.; Lu, P.; Wang, Z.B.; Li, A.Z.; Zeng, Q.H.; Chen, L.W. Optimizing a standard spectral measurement protocol to enhance the quality of soil spectra: Exploration of key variables in lab-based VNIR-SWIR spectral measurement. Remote Sens. 2022, 14, 1558. [Google Scholar] [CrossRef]
  39. Chaves, M.M.; Oliveira, M.M. Mechanisms underlying plant resilience to water deficits: Prospects for water-saving agriculture. J. Exp. Bot. 2004, 55, 2365–2384. [Google Scholar] [CrossRef] [Green Version]
  40. Li, C.; Adu, B.; Wu, J.; Qin, G.; Li, H.; Han, Y. Spatial and temporal variations of drought in Sichuan Province from 2001 to 2020 based on modified Temperature Vegetation Dryness Index (TVDI). Ecol. Indic. 2022, 139, 108883. [Google Scholar] [CrossRef]
  41. Keshavarz, M.R.; Vazifedoust, M.; Alizadeh, A. Drought monitoring using a soil wetness deficit index (SWDI) derived from MODIS satellite data. Agric. Water Manag. 2014, 132, 37–45. [Google Scholar] [CrossRef]
  42. Petropoulos, G.; Carlson, T.; Wooster, M.; Islam, S. A review of Ts/VI remote sensing based methods for the retrieval of land surface energy fluxes and soil surface moisture. Prog. Phys. Geogr. 2009, 33, 224–250. [Google Scholar] [CrossRef] [Green Version]
  43. Carlson, T.N.; Gillies, R.R.; Perry, E.M. A method to make use of thermal infrared temperature and NDVI measurements to infer surface soil water content and fractional vegetation cover. Remote Sens. Rev. 1994, 9, 161–173. [Google Scholar] [CrossRef]
  44. Tucker, C.J. Red and photographic infrared linear combinations for monitoring vegetation. Remote Sens. Environ. 1979, 8, 127–150. [Google Scholar] [CrossRef] [Green Version]
  45. Roujean, J.L.; Breon, F.M. Estimating PAR absorbed by vegetation from bidirectional reflectance measurements. Remote Sens. Environ. 1995, 51, 375–384. [Google Scholar] [CrossRef]
  46. Gamon, J.A.; Penuelas, J.; Field, C.B. A narrow-waveband spectral index that tracks diurnal changes in photosynthetic efficiency. Remote Sens. Environ. 1992, 41, 35–44. [Google Scholar] [CrossRef]
  47. Penelas, J.; Filella, I. Visible and near-infrared reflectance techniques for diagnosing plant physiological status. Trends Plant Sci. 1998, 3, 151–156. [Google Scholar] [CrossRef]
  48. Daughtry, C.S.T.; Walthall, C.L.; Kim, M.S.; Brown de Conlstoun, E.; McMurtrey, J.E., III. Estimating corn leaf chlorophyll concentration from leaf and canopy reflectance. Remote Sens. Environ. 2000, 74, 229–239. [Google Scholar] [CrossRef]
  49. Gao, B. NDWI—A normalized difference water index for remote sensing of vegetation liquid water from space. Remote Sens. Environ. 1996, 58, 257–266. [Google Scholar] [CrossRef]
  50. Penuelas, J.; Isla, R.; Filella, I.; Araus, J.L. Visible and near-infrared reflectance assessment of salinity effects on barley. Crop Sci. 1997, 37, 198–202. [Google Scholar] [CrossRef]
  51. Vicente-Serrano, S.M.; Beguería, S.; López-Moreno, J.I. A multi-scalar drought index sensitive to global warming: The standardized precipitation evapotranspiration index–SPEI. J. Clim. 2010, 23, 1696–1718. [Google Scholar] [CrossRef] [Green Version]
  52. Rouault, M.; Richard, Y. Intensity and spatial extension of droughts at different time scale since 1921 in South Africa. In Proceedings of the EGS—AGU—EUG Joint Assembly, Nice, France, 6–11 April 2003. [Google Scholar]
  53. Cao, X.J.; Yin, G.H.; Gu, J.; Ma, N.N.; Wang, Z.H. Application of WNN-PSO model in drought prediction at crop growth stages: A case study of spring maize in semi-arid regions of northern China. Comput. Electro. Agric. 2022, 199, 107155. [Google Scholar]
  54. Yang, X.J.; Liu, Y.; Bai, W.; Liu, B.C. Spatiotemporal assessment of drought related to soybean production and sensitivity analysis in Northeast China. J. Appl. Meteorol. Clim. 2017, 56, 937–952. [Google Scholar] [CrossRef]
  55. Tong, Q.X.; Zhang, B.; Zheng, L.F. Hyperspectral Remote Sensing; Higher Education Press: Beijing, China, 2006. [Google Scholar]
  56. Swain, P.H.; Davis, S.M. Remote Sensing: The Quantitative Approach; McGrow-Hill Inc: New York, NY, USA, 1978. [Google Scholar]
  57. Wang, C.; Linderholm, H.W.; Song, Y.; Wang, F.; Ren, G. Impacts of drought on maize and soybean production in northeast china during the past five decades. Int. J. Environ. Res. Public Health 2020, 17, 2459. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Huang, Y.; Li, J.; Wang, Y.; Zhang, X.Y.; Li, X.Y. Impact simulation of drought on maize growth and yield in different growth stages. J. Agric. Catastrophol. 2019, 9, 47–49. [Google Scholar]
  59. Dong, S.; Jiang, Y.; Dong, Y.; Wang, L.; Wang, W.; Ma, Z.; Yan, C.; Ma, C.; Liu, L. A study on soybean responses to drought stress and rehydration. Saudi J. Biol. Sci. 2019, 26, 2006–2017. [Google Scholar] [CrossRef]
  60. Buezo, J.; Sanz-Saez, Á.; Moran, J.F.; Soba, D.; Aranjuelo, I.; Esteban, R. Drought tolerance response of high-yielding soybean varieties to mild drought: Physiological and photochemical adjustments. Physiol. Plant. 2019, 166, 88–104. [Google Scholar] [CrossRef] [Green Version]
  61. Horler, D.N.H.; Dockray, M.; Barber, J. The red edge of plant leaf reflectance. Int. J. Remote Sens. 1983, 4, 273–288. [Google Scholar] [CrossRef]
  62. Zhao, H.; Xu, Z.; Zhao, J.; Huang, W. A drought rarity and evapotranspiration-based index as a suitable agricultural drought indice. Ecol. Indic. 2017, 82, 530–538. [Google Scholar] [CrossRef]
  63. Tian, L.Y.; Yuan, S.S.; Quiring, S.M. Evaluation of six indices for monitoring agricultural drought in the south-central United States. Agric. For. Meteorol. 2018, 249, 107–119. [Google Scholar] [CrossRef]
  64. Cao, Y.J.; Chen, S.B.; Wang, L.; Zhu, B.X.; Yu, Y. An agricultural drought index for assessing droughts using a water balance method: A case study in jilin province, northeast china. Remote Sens. 2019, 11, 1066. [Google Scholar] [CrossRef] [Green Version]
  65. Leal, L.H.C.; Violante, F.G.M.; Carvalho, L.G.; Garrido, B.C.; Rego, E.C.P.; Sarmanho, G.F.; Rocha, W.F.C. A new methodology for proficiency testing scheme interpretation based on residual analysis. Accredit. Qual. Assur. 2022, 27, 349–357. [Google Scholar] [CrossRef]
  66. Danson, F.M.; Steven, M.D.; Malthus, T.J.; Clark, J.A. High-spectral resolution data for determining leaf water content. Int. J. Remote Sens. 1992, 13, 461–470. [Google Scholar] [CrossRef]
  67. Patel, C.K.N.; Tam, A.C. Optical absorption coefficients of water. Nature 1979, 280, 302–304. [Google Scholar] [CrossRef]
  68. Holben, B.N.; Schutt, J.B.; Mcmurtrey, J. Leaf water stress detection utilizing thematic mapper bands 3, 4 and 5 in soybean plants. Int. J. Remote Sens. 1983, 4, 289–297. [Google Scholar] [CrossRef]
  69. Everitt, J.H.; Nixon, P.R. Canopy reflectance of two drought-stressed shrubs. Photogramm. Eng. Remote Sens. 1986, 52, 1189–1192. [Google Scholar]
  70. Liang, S.; Li, X.; Wang, J. Advanced Remote Sensing: Terrestrial Information Extraction and Applications; Academic Press: Cambridge, MA, USA, 2012. [Google Scholar]
  71. Houborg, R.; Fisher, J.B.; Skidmore, A.K. Advances in remote sensing of vegetation function and traits. Int. J. Appl. Earth Obs. Geoinf. 2015, 43, 1–6. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The study area. (a) Geographical position of Northeast China; (b) main upland crop types and range in Northeast China; and (c) digital elevation model data (DEM) and the meteorological stations in Northeast China.
Figure 1. The study area. (a) Geographical position of Northeast China; (b) main upland crop types and range in Northeast China; and (c) digital elevation model data (DEM) and the meteorological stations in Northeast China.
Remotesensing 15 01915 g001
Figure 2. Flowchart of the study.
Figure 2. Flowchart of the study.
Remotesensing 15 01915 g002
Figure 3. Spectral reflectance of plant canopy simulated by the PROSAIL model. (a) Simulated spectral reflectance under different EWT; (b) sensitivity analysis of EWT; (c) simulated spectral reflectance under different Cab; (d) sensitivity analysis of Cab; (e) simulated spectral reflectance under different LAIs; and (f) sensitivity analysis of the LAI.
Figure 3. Spectral reflectance of plant canopy simulated by the PROSAIL model. (a) Simulated spectral reflectance under different EWT; (b) sensitivity analysis of EWT; (c) simulated spectral reflectance under different Cab; (d) sensitivity analysis of Cab; (e) simulated spectral reflectance under different LAIs; and (f) sensitivity analysis of the LAI.
Remotesensing 15 01915 g003
Figure 4. Spectral reflectance of crop canopy at different water stress levels. (a) Canopy spectral reflectance of soybean at the branching stage and experimental photographs; and (b) canopy spectral reflectance of maize at the milking stage and experimental photographs (the missing parts of spectral reflectance are the atmospheric water vapor absorption channel during field measurement).
Figure 4. Spectral reflectance of crop canopy at different water stress levels. (a) Canopy spectral reflectance of soybean at the branching stage and experimental photographs; and (b) canopy spectral reflectance of maize at the milking stage and experimental photographs (the missing parts of spectral reflectance are the atmospheric water vapor absorption channel during field measurement).
Remotesensing 15 01915 g004
Figure 5. Sensitivity analysis of crop canopy spectral reflectance at different water stress levels. (a) Soybean; and (b) maize.
Figure 5. Sensitivity analysis of crop canopy spectral reflectance at different water stress levels. (a) Soybean; and (b) maize.
Remotesensing 15 01915 g005
Figure 6. The MVDI–LST triangle space and definition of the mTVDI.
Figure 6. The MVDI–LST triangle space and definition of the mTVDI.
Remotesensing 15 01915 g006
Figure 7. (a–j) Recognition ability of vegetation indices under different water stress levels during the main growth stage of the crop. (a) The NDVI; (b) RDVI; (c) EVI; (d) PRI; (e) SIPI; (f) MCARI; (g) WI; (h) NDII; (i) NDWI; and (j) MVDI. The four water stress levels are as follows: no water stress (N: 70–80% FC), light water stress (L: 50–60% FC), moderate water stress (M: 40–50% FC), and severe water stress (S: 30–40% FC).
Figure 7. (a–j) Recognition ability of vegetation indices under different water stress levels during the main growth stage of the crop. (a) The NDVI; (b) RDVI; (c) EVI; (d) PRI; (e) SIPI; (f) MCARI; (g) WI; (h) NDII; (i) NDWI; and (j) MVDI. The four water stress levels are as follows: no water stress (N: 70–80% FC), light water stress (L: 50–60% FC), moderate water stress (M: 40–50% FC), and severe water stress (S: 30–40% FC).
Remotesensing 15 01915 g007
Figure 8. The feature spaces of the TVDI and mTVDI in Northeast China during the growing season (May–September) of 2017. (a,c,e,g,i) Dry and wet edges of the NDVI–LST feature space; and (b,d,f,h,j) dry and wet edges of the MVDI–LST feature space.
Figure 8. The feature spaces of the TVDI and mTVDI in Northeast China during the growing season (May–September) of 2017. (a,c,e,g,i) Dry and wet edges of the NDVI–LST feature space; and (b,d,f,h,j) dry and wet edges of the MVDI–LST feature space.
Remotesensing 15 01915 g008
Figure 9. Spatio-temporal drought distribution in Northeast China, indicated by different drought indices, throughout the growth period (May–September) in 2017.
Figure 9. Spatio-temporal drought distribution in Northeast China, indicated by different drought indices, throughout the growth period (May–September) in 2017.
Remotesensing 15 01915 g009
Figure 10. Comparison of the mTVDI with precipitation, PET, and the CWDI during the growing season (May−September) of 2017.
Figure 10. Comparison of the mTVDI with precipitation, PET, and the CWDI during the growing season (May−September) of 2017.
Remotesensing 15 01915 g010
Figure 11. Validation of the mTVDI and SPEI3 with the CDD (red marks represent the recorded drought events).
Figure 11. Validation of the mTVDI and SPEI3 with the CDD (red marks represent the recorded drought events).
Remotesensing 15 01915 g011
Figure 12. Comparison and residual analysis of the HVDI and MVDI. (ad) The comparison between the HVDI and MVDI in the early, middle, late, and mature stages of crops; and (eh) the residual analysis between the HVDI and MVDI in the early, middle, late and mature stages of crops.
Figure 12. Comparison and residual analysis of the HVDI and MVDI. (ad) The comparison between the HVDI and MVDI in the early, middle, late, and mature stages of crops; and (eh) the residual analysis between the HVDI and MVDI in the early, middle, late and mature stages of crops.
Remotesensing 15 01915 g012
Table 1. Parameter settings used in the PROSAIL model.
Table 1. Parameter settings used in the PROSAIL model.
ModelParameterQuantityValueStepUnit
PROSPECTNLeaf structure parameter2.5--
CabChlorophyll a+b content10~1001μg/cm2
CarCarotenoid content8-μg/cm2
CWEquivalent water thickness0.005~0.030.001cm
CmDry matter content0.009-g/cm2
CbpBrown pigments content0--
SAILLAILeaf area index1~60.1-
LIDFLeaf inclination distribution functionBeta distribution--
SL.Hot spot parameter0.01--
sza of θSSolar zenith angle30-deg
vza of θVViewing zenith angle10-deg
raa of θSVRelative azimuth angle90-deg
Table 2. Dates of spectral reflectance measurements.
Table 2. Dates of spectral reflectance measurements.
StageMaizeSoybean
Growth PeriodAcquisition Time
(2021)
Growth PeriodAcquisition Time
(2020)
EarlyJointing24 JuneSeedling16 June, 28 June
MiddleTasseling14 JulyBranch9 July, 15 July
LateMilking9 AugustPod bearing21 July, 25 July, 6 August
MatureMaturing30 AugustPod filling12 August, 20 August, 26 August
Table 3. Vegetation indices used in this study.
Table 3. Vegetation indices used in this study.
AbbreviationIndex NameEquationReferences
Structural Indices
1. NDVINormalized difference vegetation index ( R 800 R 680 ) / ( R 800 + R 680 ) Tucker [44]
2. RDVIRenormalized difference vegetation index ( R 800 R 670 ) / ( R 800 + R 670 ) Roujean [45]
3. EVIEnhanced vegetation index 2.5 ( R 800 R 670 ) / ( R 800 + 6 R 670 7 R 450 + 1 ) Huete [9]
Chlorophyll Indices
4. PRIPhotochemical reflectance index ( R 570 R 531 ) / ( R 570 + R 531 ) Gamon [46]
5. SIPIStructure-intensive pigment index ( R 800 R 445 ) / ( R 800 + R 680 ) Penelas [47]
6. MCARIModified chlorophyll absorption ratio index [ ( R 700 R 670 ) 0.2 ( R 700 R 550 ) ] / ( R 700 / R 670 ) Daughtry [48]
Water Indices
7. NDIINormalized difference infrared index ( R 819 R 1649 ) / ( R 819 + R 1649 ) Jackson [12]
8. NDWINormalized difference water index ( R 860 R 1240 ) / ( R 860 + R 1240 ) Gao [49]
9. WIWater index R 900 / R 970 Penuelas [50]
Table 4. The normalized mean distance between different water stress levels during the main growth stages of the crop (empty cell indicates that the index could not distinguish different water stress levels).
Table 4. The normalized mean distance between different water stress levels during the main growth stages of the crop (empty cell indicates that the index could not distinguish different water stress levels).
TypeIndexEarly StageMiddle Stage
N–LL–MM–SAverageN–LL–MM–SAverage
Structural IndicesNDVI2.770.932.492.061.270.860.971.03
RDVI2.520.491.441.481.240.660.250.72
EVI2.560.581.671.601.260.750.430.81
Chlorophyll IndicesPRI3.401.222.372.331.221.112.251.53
SIPI2.140.731.701.521.290.830.760.96
MCARI0.830.961.821.20----
Water IndicesNDII3.354.184.113.880.870.951.191.00
NDWI2.152.932.832.640.900.621.511.01
WI2.732.281.602.200.790.431.160.79
This Study IndexMVDI4.324.475.034.611.590.891.451.31
TypeIndexLate stageMature stage
N–LL–MM–SAverageN–LL–MM–SAverage
Structural IndicesNDVI1.552.642.442.21----
RDVI2.571.221.551.781.221.091.431.25
EVI2.301.761.932.00----
Chlorophyll IndicesPRI1.000.441.961.13----
SIPI1.583.272.382.41----
MCARI--------
Water IndicesNDII0.410.050.270.24----
NDWI0.840.780.520.71----
WI0.820.860.380.69----
This Study IndexMVDI2.552.652.832.681.601.391.401.46
Table 5. Pearson correlation coefficients between the mTVDI, TVDI, and SPEIs.
Table 5. Pearson correlation coefficients between the mTVDI, TVDI, and SPEIs.
MonMethodHeilongjiangJilinLiaoningNeimengguMean
SPEI1SPEI3SPEI6SPEI1SPEI3SPEI6SPEI1SPEI3SPEI6SPEI1SPEI3SPEI6
MaymTVDI−0.83−0.86−0.84−0.39−0.75−0.76−0.62−0.78−0.88−0.61−0.80−0.82−0.75
TVDI−0.80−0.81−0.80−0.38−0.65−0.66−0.62−0.68−0.77−0.60−0.72−0.74−0.69
JunmTVDI−0.36−0.39−0.48−0.52−0.69−0.68−0.64−0.72−0.81−0.69−0.78−0.78−0.63
TVDI−0.37−0.32−0.40−0.48−0.57−0.58−0.66−0.66−0.46−0.50−0.52−0.48−0.55
JulmTVDI−0.68−0.65−0.68−0.86−0.76−0.81−0.72−0.69−0.51−0.75−0.70−0.66−0.71
TVDI−0.66−0.62−0.65−0.68−0.68−0.72−0.64−0.50−0.38−0.66−0.60−0.58−0.61
AugmTVDI−0.54−0.79−0.80−0.45−0.59−0.52−0.73−0.67−0.66−0.57−0.68−0.66−0.64
TVDI−0.47−0.79−0.770.29−0.35−0.61−0.68−0.66−0.69−0.29−0.60−0.69−0.56
SepmTVDI−0.77−0.74−0.81−0.62−0.69−0.89−0.65−0.76−0.83−0.68−0.73−0.84−0.75
TVDI−0.75−0.63−0.73−0.65−0.62−0.76−0.66−0.71−0.69−0.69−0.65−0.72−0.69
Table 6. The documented records of the major drought events from the China Drought Dataset (CDD).
Table 6. The documented records of the major drought events from the China Drought Dataset (CDD).
RegionDrought Month
HeilongjiangSeptember 2008, September 2010
JilinJuly 2009, August 2009, September 2009, September 2010
LiaoningSeptember 2008, July 2009, August 2009, September 2009
NeimengguJune 2007, July 2007, September 2009
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Dai, R.; Chen, S.; Cao, Y.; Zhang, Y.; Xu, X. A Modified Temperature Vegetation Dryness Index (mTVDI) for Agricultural Drought Assessment Based on MODIS Data: A Case Study in Northeast China. Remote Sens. 2023, 15, 1915. https://doi.org/10.3390/rs15071915

AMA Style

Dai R, Chen S, Cao Y, Zhang Y, Xu X. A Modified Temperature Vegetation Dryness Index (mTVDI) for Agricultural Drought Assessment Based on MODIS Data: A Case Study in Northeast China. Remote Sensing. 2023; 15(7):1915. https://doi.org/10.3390/rs15071915

Chicago/Turabian Style

Dai, Rui, Shengbo Chen, Yijing Cao, Yufeng Zhang, and Xitong Xu. 2023. "A Modified Temperature Vegetation Dryness Index (mTVDI) for Agricultural Drought Assessment Based on MODIS Data: A Case Study in Northeast China" Remote Sensing 15, no. 7: 1915. https://doi.org/10.3390/rs15071915

APA Style

Dai, R., Chen, S., Cao, Y., Zhang, Y., & Xu, X. (2023). A Modified Temperature Vegetation Dryness Index (mTVDI) for Agricultural Drought Assessment Based on MODIS Data: A Case Study in Northeast China. Remote Sensing, 15(7), 1915. https://doi.org/10.3390/rs15071915

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