Acessibilidade / Reportar erro

Temporal analysis of drought coverage in a watershed area using remote sensing spectral indexes

Abstract

The development of several time series analysis programs using satellite images has provided many applications based on resources from geostatistics field. Currently, the use of statistical tests applied to vegetation indexes has enabled the analysis of different natural phenomena, such as drought events in watershed areas. The objective of this article is to provide a comparative analysis between NDVI and EVI vegetation index data made available by MOD13Q1 project of MODIS sensor for drought mapping using vegetation condition index (VCI) in the Serra Azul stream sub-basin, MG. The methodology adopted the Cox-Stuart statistical test for seasonality analysis and Pearson's linear correlation to verify the influence of different indexes on delimitation of drought in a watershed. The results indicated the NDVI vegetation index as more efficient than EVI in spatial characterization of studied watershed region, mainly in identification of seasonality. The VCI proved to be highly feasible for monitoring drought in study period between 2013 and 2018, allowing the effective delimitation of drought conditions in the Serra Azul stream sub-basin. In addition, the effectiveness of MODIS sensor data in characterizing drought events that affected the study area was proven.

Keywords:
Geostatistics; Remote Sensing Seasonality Drought Index; Vegetation Index

INTRODUCTION

The development of several Earth observation programs has provided, through spectral analyzes of the environment, increasingly adequate responses to different natural variations and anthropogenic actions that occur on planet´s surface. One of the important factors that helped in expansion of this knowledge is data recording in time series that make it possible to analyze the records of patterns related to different processes, such as biogeophysical, meteorological and other cycles (POTTER et al., 2003POTTER, C. et al. Major disturbance events in terrestrial ecosystems detected using global satellite data sets. Global Change Biology, v. 9, n. 7, p. 1005-1021, 2003. https://doi.org/10.1046/j.1365-2486.2003.00648.x
https://doi.org/10.1046/j.1365-2486.2003...
).

The concept of time series comes from the set of observations recorded over a given time interval. In the scope of environmental sciences, the availability of this type of data allows the development of several applications, such as rain or air temperature analysis using probabilistic or holistic approaches. In recent times, remote sensing technology has contributed significantly to these analyzes due to its characteristic ability to produce databases containing information from Earth's surface in the form of time series (DAVIES; CHATFIELD, 1990DAVIES, N.; CHATFIELD, C. The analysis of time series: an introduction. 6. ed. New York Washington: CHAPMAN & HALL/CRC, 1990. v. 74. https://doi.org/10.2307/3619403
https://doi.org/10.2307/3619403...
; SAUSEN; LACRUZ, 2015SAUSEN, T. M.; LACRUZ, M. S. P. Sensoriamento remoto para desastres. 1. ed. São Paulo: Oficina de Textos, 2015.).

In time series analysis, remote sensing data applications can be performed from the perspective of earth's surface analysis, through variables that are used directly in drought monitoring such as vegetation indexes (DECHANT; MORADKHANI, 2014DECHANT, C. M.; MORADKHANI, H. Toward a reliable prediction of seasonal forecast uncertainty: Addressing model and initial condition uncertainty with ensemble data assimilation and Sequential Bayesian Combination. Journal of Hydrology, v. 519, n. PD, p. 2967-2977, 2014. http://dx.doi.org/10.1016/j.jhydrol.2014.05.045
http://dx.doi.org/10.1016/j.jhydrol.2014...
, 2015). Vegetation indexes are mathematical formulations that use remote sensing spectral data to estimate the behavior of vegetation cover in a region. These formulations allow to analyze the vegetation activity as well as foliage variation in terms of seasonality (BONIFACIO; DUGDALE; MILFORD, 1993BONIFACIO, R.; DUGDALE, G.; MILFORD, J. R. Sahelian rangeland production in relation to rainfall estimates from Meteosat. International Journal of Remote Sensing, v. 14, n. 14, p. 2695-2711, 1993.; FORMAGGIO, SANCHES, 2017).

The analysis of spatial variability of data to obtain estimates in non-sampled locations is performed using geostatistics. In this process, exploratory techniques based on descriptive analysis of calculations from descriptive statistics are used. In time series research field, trend analysis is of great relevance in environmental matters, as they are able to identify the influence of seasonality on some parameters under study (SOARES, 2000SOARES, A. Geoestatística para as ciências da terra e do ambiente. Rio de Janeiro: Instituto Superior Técnico, 2000.; VIEIRA et al., 1983VIEIRA, S. et al. Geoestatistical theory and application to variability of some agronomical properties. Hilgardia, v. 51, n. 3, p. 1-75, 1983. https://doi.org/10.3733/hilg.v51n03p075
https://doi.org/10.3733/hilg.v51n03p075 ...
). In these types of studies, the MODIS sensor has been widely used in research involving sugarcane species, with use of drought rates (DUFT; PICOLI, 2018DUFT, D. G.; PICOLI, M. C. A. Uso de imagens do sensor MODIS para identificação da seca na cana-de-açúcar através de índices espectrais. Scientia Agraria, v. 19, n. 1, p. 52, 2018. https://doi.org/10.5380/rsa.v19i1.54005
https://doi.org/10.5380/rsa.v19i1.54005...
), as well as in studies on pasture areas (CUNHA et al., 2017CUNHA, A. P. M. A. et al. Avaliação de indicador para o monitoramento dos impacos da seca em áreas de pastagens no Semiárido do Brasil. Revista Brasileira de Cartografia, v. 69, n. 1, p. 89-106, 2017.).

Currently, several studies use time series analysis with vegetation index data to monitor natural disasters, such as drought and water scarcity in river basins. In studies related to drought, the vegetation indexes are applied to monitor the event through the Vegetation Condition Index (VCI), created by Kogan (1995KOGAN, F. N. Application of vegetation index and brightness temperature for drought detection. Advances in Space Research, v. 15, n. 11, p. 91-100, 1995. https://doi.org/10.1016/0273-1177(95)00079-T
https://doi.org/10.1016/0273-1177(95)000...
), and applied using MODIS sensor data in event analysis of drought in biological reserves (BRANCO, 2016BRANCO, E. R. F. Ocorrências de seca e tendências da vegetação na reserva biológica de sooretama e zona de amortecimento, no estado do Espírito Santo, Brasil. Dissertação de mestrado ed. Jerônimo Monteiro: Universidade Federal do Espirito Santo, 2016.), as well as in comparison between drought delimitation indexes (JIAO et al., 2019JIAO, W. et al. A new station-enabled multi-sensor integrated index for drought monitoring. Journal of Hydrology, v. 574, n. April, p. 169-180, 2019. https://doi.org/10.1016/j.jhydrol.2019.04.037
https://doi.org/10.1016/j.jhydrol.2019.0...
; ZHANG et al., 2017ZHANG, L. et al. Studying drought phenomena in the Continental United States in 2011 and 2012 using various drought indices. Remote Sensing of Environment, v. 190, p. 96-106, 2017. http://dx.doi.org/10.1016/j.rse.2016.12.010
http://dx.doi.org/10.1016/j.rse.2016.12....
). In addition, this technique is applied with the use of other remote sensing data, such as NOAA-NESDIS, in analyzes of vegetation tendency to drought variation (XU et al., 2020XU, Z. et al. Trends in Global Vegetative Drought from Long-Term Satellite Remote Sensing Data. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, v. 13, p. 815-826, 2020. https://doi.org/10.1109/JSTARS.2020.2972574
https://doi.org/10.1109/JSTARS.2020.2972...
) and variability of vegetation and temperature indexes in Brazilian regions (GOMES et al., 2019GOMES, A. R. S. et al. Estudo da Relação entre a Variabilidade dos índices de Vegetação e Temperatura da Região Nordeste do Brasil. Revista Brasileira de Meteorologia, v. 34, n. 3, p. 359-368, 2019. https://doi.org/10.1590/0102-7786343051
https://doi.org/10.1590/0102-7786343051...
).

In remote sensing researches, within geostatistics scope, the Cox-Stuart test can assess time series tendency to present seasonality through the application of a Ho hypothesis. The trend is calculated through the differences between pairs of time series variables, extracted from the original sample of time series. A positive or negative sign is associated with each pair, and equal values are eliminated. In the Ho hypothesis, the total number of negative and positive signs is expected to be similar, so that a series is considered to be trendless (COX; STUART, 1955COX, D. R.; STUART, A. Some quick sign tests for trend in location and dispersion. Biometrika, v. 42, n. 1/2, p. 80, 1955.; DETZEL et al., 2011DETZEL, D. H. M. et al. Estacionariedade das afluências às usinas hidrelétricas brasileiras. Revista Brasileira de Recursos Hídricos, v. 16, n. 3, p. 95-111, 2011. https://doi.org/10.21168/rbrh.v16n3.p95-111
https://doi.org/10.21168/rbrh.v16n3.p95-...
). The application of this technique has been developed in studies that use data from MODIS sensor (ARANTES et al., 2017ARANTES, T. B. et al. Effectiveness of BFAST algorithm to characterize time series of dense forest, agriculture and pasture in the amazon region. Theoretical and applied engineering, v. 1, n. 1, p. 10-19, 2017.) and in assessment of vegetation behavior for environmental planning and sustainable management with SPOT data (CHAVES; MATAVELI; JUSTINO, 2014CHAVES, M. E. D.; MATAVELI, G. A. V.; JUSTINO, R. C. Uso da modelagem estatística para monitoramento da vegetação no Parque Nacional da Serra da Canastra, Minas Gerais. Caderno de Geografia, v. 24, n. 1, p. 120-132, 2014. https://doi.org/10.5752/P.2318-2962.2014v24nespp120
https://doi.org/10.5752/P.2318-2962.2014...
).

The objective of this work is to provide a comparative analysis between data of normalized difference vegetation index (NDVI) and enhanced vegetation index (EVI) available in MOD13Q1 project of MODIS sensor for drought mapping by determining the vegetation condition index (VCI). The adopted approach involved the application of Cox-Stuart statistical test for seasonality analysis and the use of Pearson's linear correlation to verify the indexes influence in the Serra Azul stream basin drought delimitation.

METHODS

The drought variation was identified through data from the MOD13Q1 product, the methodology used was divided into five analyzes, namely: selection of the study area, obtaining data from the MODIS sensor, analysis of the vegetation indexes time series, calculation of the vegetation condition index and correlation of the results generated with meteorological data (Figure 1).

Figure 1
Flowchart of the applied methodology.

Drought analysis in the Serra Azul stream sub-basin

The study area is the Serra Azul stream sub-basin which covers 447.83 km² of drainage area along the municipalities of Mateus Leme, Igarapé, Juatuba and Itaúna (Figure 2). The region is located in state of Minas Gerais, Brazil between 20º 15' and 20º 00' south latitude parallels and 44º 15' and 44º 35' west longitude meridians. Inside this basin there is a supply reservoir owned by Companhia de Saneamento de Minas Gerais (COPASA), responsible for water supply of Belo Horizonte city and adjacent locations (DUTRA; BRIANEZI; COELHO, 2020DUTRA, D. J.; ELMIRO, M. A. T.; GARCIA, R. A. Comparative analysis of methods applied in vegetation cover delimitation using Landsat 8 images. Sociedade & Natureza, v. 32, n. July, p. 699-710, 9 out. 2020. https://doi.org/10.14393/SN-v32-2020-56139
https://doi.org/10.14393/SN-v32-2020-561...
; DUTRA; ELMIRO; GARCIA, 2020).

According to the Köeppen classification, the region has an Aw-type climate. The “A” refers to the tropical summer, that is, areas with an average monthly temperature above 18oC. The “w” means that the region has a period of four to six months of dry season (DUBREUIL et al., 2018DUBREUIL, V. et al. Os tipos de climas anuais no Brasil: Uma aplicação da classificação de Köppen de 1961 a 2015. Confins. Revue franco- brésilienne de géographie/Revista franco-brasilera de geografia, v. 37, 2018.). According to IBGE climate classification (IBGE, 2019), the region is inserted in Tropical Zone Central Brazil Sub-hot, semi-humid with four to five dry months that characterizes it as semi-humid.

Figure 2
Serra Azul stream sub-basin location

In relation to land use, the study area is located in the Alto São Francisco region, which is characterized by Rock vegetation and Atlantic forest in the neighborhood of water supply reservoir. In addition, the region is under a large intervention of human activities, due to the expansion of urban areas, agricultural activities and mining (Figure 3). These activities cause the suppression of vegetation, impermeability of the soil and several impacts in the region, such as increased surface temperature and water scarcity events (DUTRA; BRIANEZI; COELHO, 2020DUTRA, D. J.; ELMIRO, M. A. T.; GARCIA, R. A. Comparative analysis of methods applied in vegetation cover delimitation using Landsat 8 images. Sociedade & Natureza, v. 32, n. July, p. 699-710, 9 out. 2020. https://doi.org/10.14393/SN-v32-2020-56139
https://doi.org/10.14393/SN-v32-2020-561...
; DUTRA; ELMIRO; GARCIA, 2020; MINAS-GERAIS, 2015).

Figure 3
Variation of land use in the Serra Azul stream sub-basin, with the presence of the following points: 1 (Serra Azul reservoir); 2, 5 and 4 (forest area); 3 and 8 (urban area); 6 (agriculture); and 7 (mining area), where the green areas correspond to vegetation regions and the yellow areas are non-forest

MODIS sensor data

The data from MODIS sensor was obtained on United States Geological Survey (USGS) platform. Originally, the images are provided in Hierarchical Data Format (HDF) ,sinusoidal cartographic projection and 250m spatial resolution. In this way, QGIS 3.4.6 software tools were used to convert data to Tagged Image File Format (TIFF) and UTM (Universal Transverse Mercator) cartographic projection, 23 South Zone and SIRGAS 2000 horizontal datum. The project data are expressed in the range of -10000 to 10000, therefore it was necessary to reschedule data using a scale factor (0.0001) to obtain vegetation indexes values in conventional numerical range from -1 to 1.

Analysis of vegetation indexes time series: Cox-Stuart test

For application of the proposed statistical analyzes, a monthly based time series was created for the period between 01/01/2013 and 12/31/2018, corresponding to 72 months. Vegetative vigor information for each pixel was extracted from this monthly time series, referring to location of meteorological stations of the Instituto Nacional de Meteorologia (INMET) and Agência Nacional de Águas (ANA). Figure 4 shows the location of these stations.

In order to verify the existence of some tendency to seasonality in relationship between amplitude and mean of time series data, the Cox-Stuart test was applied, using RStudio software. This test developed by Cox and Stuart (1955COX, D. R.; STUART, A. Some quick sign tests for trend in location and dispersion. Biometrika, v. 42, n. 1/2, p. 80, 1955.) was applied in studies of Morettin and Toloi (2006MORETTIN, P. A.; TOLOI, C. M. C. Análise de séries temporais. 4. ed. New Jersey: Blucher, 2006.). According to Morettin and Toloi (1981), in applications of time series analysis referring to vegetation index data, the test aims to assess whether or not a time series presents trends in seasonality using the following hypotheses:

H0: no trend, the number of positive and negative signs are equal.

H1: trend, the number of positive and negative signs are different.

Vegetation condition index (VCI) calculation

In order to analyze the occurrence of drought using vegetation condition index (VCI), the data referring to vegetation index that showed the best response in statistical analyzes were organized according to seasons (Table 1). Dates were sorted in ascending order according Julian days (continuous scale of the days of the year).

Figure 4
Location of ANA and INMET stations

Table 1
Seasons, Julian days and Gregorian days referring to vegetation indexes images obtained from MODIS sensor.

Vegetation changes are not easily identified with direct use of vegetation indexes, so VCI is used, as a more effective index, in order to identify the productivity of an ecosystem. The VCI measures drought condition according to interference of climate in vegetative vigor of a region (DU et al., 2013DU, L. et al. A comprehensive drought monitoring method integrating MODIS and TRMM data. International Journal of Applied Earth Observation and Geoinformation, v. 23, n. 1, p. 245-253, 2013. http://dx.doi.org/10.1016/j.jag.2012.09.010
http://dx.doi.org/10.1016/j.jag.2012.09....
). For preparation of VCI, the averages of NDVI and EVI were calculated for each season throughout the analyzed period (2013 a 2018), according to Equation 1.

X ¯ = 1 N I i = 1 N I X i v (1)

where

NI = number of images for each season;

X¯ = average of a numerical data set; and

Xiv = vegetation index values.

In order to obtain the drought occurrences over the full studied period, the methodology proposed by Kogan (1995KOGAN, F. N. Application of vegetation index and brightness temperature for drought detection. Advances in Space Research, v. 15, n. 11, p. 91-100, 1995. https://doi.org/10.1016/0273-1177(95)00079-T
https://doi.org/10.1016/0273-1177(95)000...
a) was applied, according to Equation 2. First, the average images were grouped according to season, throughout the entire time series. Then, the average of all pixels referring to maximum and minimum values for each season was calculated for each type of vegetation index.

V C I = X ¯ e X ¯ min e X ¯ max e X ¯ min e * 100 (2)

where

VCI = Vegetation Condition Index (%);

X¯e = average per season of vegetation index for a given year;

X¯mine = overall average of minimum values of vegetation index for a given season; and

X¯maxe = overall average of maximum values of vegetation index for a given season.

The VCI, presenting values from zero to one hundred, allows analysis of vegetation index data (NDVI and EVI) in a short period of time, as well as the evaluation of changes in long term. According to Kogan (1995KOGAN, F. N. Application of vegetation index and brightness temperature for drought detection. Advances in Space Research, v. 15, n. 11, p. 91-100, 1995. https://doi.org/10.1016/0273-1177(95)00079-T
https://doi.org/10.1016/0273-1177(95)000...
a), the closer de VCI value to zero the greater the influence of drought phenomenon in a region. Following this approach, drought is classified according to Table 2.

Table 2
Class ranges of Vegetation Condition Index (VCI) values and the corresponding classification

Correlation between meteorological data and NDVI, EVI and VCI indexes

A correlation analysis between vegetation indexes and meteorological data was used to select the best vegetation index (NDVI or EVI) to be adopted in the VCI calculation. Three meteorological variables (air temperature, evapotranspiration and rainfall) were used in this analysis for selecting the best vegetation index. Meteorological data were obtained from the stations of INMET and ANA. Equation 3, proposed by Pearson (1982PEARSON, K. The grammar of science. Nature, v. 46, n. 1185, p. 247-247, 1982. https://doi.org/10.1038/046247b0
https://doi.org/10.1038/046247b0...
) and Pearson, Fischer and Inman (1994) was used to calculate the correlation.

r = i = 1 n ( x i x ¯ ) ( y i y ¯ ) i = 1 n ( x i x ¯ ) 2 i = 1 n ( y i y ¯ ) 2 (3)

where

r = Pearson’s correlation coefficient;

x¯ = x sample average; and

y¯= y sample average.

According to Figueiredo Filho and Silva Júnior (2009), Pearson's correlation coefficient (r) can present both negative and positive correlations with a range of values between -1 and 1. The closer to one, greater the linear association between variance of the indexes used and the closer to zero, lower this association. If the correlation is negative, the association will be inverse (Table 3).

Table 3
Classification of Pearson's correlation coefficient

RESULTS AND DISCUSSIONS

Best vegetation index selection for drought analysis using Cox-Stuart Test

The analysis of time series data showed the presence of seasonality patterns over the time. It was possible to observe peaks of greater vegetative vigor and moments of fall in values, caused by dry periods (Figure 5). According to Chaves, Mataveli and Justino (2014CHAVES, M. E. D.; MATAVELI, G. A. V.; JUSTINO, R. C. Uso da modelagem estatística para monitoramento da vegetação no Parque Nacional da Serra da Canastra, Minas Gerais. Caderno de Geografia, v. 24, n. 1, p. 120-132, 2014. https://doi.org/10.5752/P.2318-2962.2014v24nespp120
https://doi.org/10.5752/P.2318-2962.2014...
), this variation is directly associated with precipitation, since the beginning of rains contributes to the increase of vegetative vigor in a region.

Figure 5
Variation of NDVI and EVI in the Serra Azul stream sub-basin for period from 2013 to 2018

Figure 6
Difference of transformed series from Cox-Stuart Test, (a) EVI and (b) NDVI

Through the analysis of difference between pairs of variables obtained from stationary series transformed from original vegetation index data, it was identified that NDVI registered a greater tendency to present seasonality patterns in comparison with EVI data (Figure 6). Furthermore, it was observed that NDVI data showed less variation between their intervals compared to EVI data.

The hypothesis results (Table 4) demonstrated that when p-value is greater than 0.05, in a 95% confidence interval, there is no statistical evidence to reject the hypothesis. That is, the EVI values of time series did not show seasonality, whereas in NDVI data, the presence of seasonality was identified, as p-value presented values below 0.05. In this case, the null hypothesis was rejected. According to Chaves, Mataveli and Justino (2014CHAVES, M. E. D.; MATAVELI, G. A. V.; JUSTINO, R. C. Uso da modelagem estatística para monitoramento da vegetação no Parque Nacional da Serra da Canastra, Minas Gerais. Caderno de Geografia, v. 24, n. 1, p. 120-132, 2014. https://doi.org/10.5752/P.2318-2962.2014v24nespp120
https://doi.org/10.5752/P.2318-2962.2014...
) and Gow et al. (2016GOW, L. J. et al. A detection problem: Sensitivity and uncertainty analysis of a land surface temperature approach to detecting dynamics of water use by groundwater-dependent vegetation. Environmental Modelling and Software, v. 85, p. 342-355, 2016. https://doi.org/10.1016/j.envsoft.2016.09.003
https://doi.org/10.1016/j.envsoft.2016.0...
), rejecting Ho hypothesis means that variation in vegetation behavior does not change substantially, indicating that there is a linearity of vegetative vigor present in region.

Table 4
Analysis of hypotheses of Cox-Stuart test for vegetation indexes obtained from MODIS sensor

In regions of lower vegetation density, such as transitions from rural to urban space, NDVI is more sensitive to seasonality. This fact occurs because in areas with intense modification of land use, the vegetative vigor has a low saturation in response to the increase in biomass, thus NDVI tends to present a better spectral response in these types of areas when compared to EVI (HUETE et al., 2002HUETE, A. et al. Overview of the radiometric and biophysical performance of the MODIS vegetation indices. Remote Sensisng of Enviroment, v. 83, p. 195-213, 2002.). Studies by Nora and Santos (2010NORA, E. L. D.; SANTOS, J. E. Análise da dinâmica sazonal de duas formaçoes florestais do bioma Mata Atlântica com base em índices de vegetação. Perspectiva, v. 34, n. 125, p. 41-51, 2010.) applied to areas characterized by anthropic changes, demonstrate that NDVI tends to present greater variability in relation to sensitivity of canopy changes when compared to EVI. This relationship was demonstrated by Branco (2016BRANCO, E. R. F. Ocorrências de seca e tendências da vegetação na reserva biológica de sooretama e zona de amortecimento, no estado do Espírito Santo, Brasil. Dissertação de mestrado ed. Jerônimo Monteiro: Universidade Federal do Espirito Santo, 2016.) where, contrary to results of present study, the EVI from MODIS sensor data presented the best performance. The main reason for these differences in results is because the region analyzed by Branco (2016) presents a lesser anthropic intervention, mainly because the object of study area corresponds to a biological reserve.

Analysis of drought in study area from NDVI time series

In the researches carried out by Cunha et al. (2017CUNHA, A. P. M. A. et al. Avaliação de indicador para o monitoramento dos impacos da seca em áreas de pastagens no Semiárido do Brasil. Revista Brasileira de Cartografia, v. 69, n. 1, p. 89-106, 2017.) and Leivas et al. (2014LEIVAS, J. F. et al. Monitoramento da seca 2011/2012 no nordeste brasileiro a partir do satélite Spot-vegetation e TRMM. Revista Engenharia Na Agricultura - Reveng, v. 22, n. 3, p. 211-221, 2014. https://doi.org/10.13083/1414-3984.v22n03a03
https://doi.org/10.13083/1414-3984.v22n0...
), in Brazilian territorial areas, it was evident that VCI can be used as an important indicator for analyzing different aspects of drought in tropical regions. This index is capable of exploring different water stress conditions in several types of landscapes. The results of vegetation condition index in the Serra Azul stream region enabled to identify the presence of seasonality over the time, distributed over the analyzed VCI classes (Figure 7). In drought months (winter and autumn) this drought was more intense in anthropized areas and in areas adjacent to Serra Azul reservoir, due to lower vegetation cover present.

Figure 7
Variation of drought classes along time series based on VCI

The class of Severe Drought showed a percentage lower than 1% in all seasons analyzed. The Extreme Drought class had largest extensions in years 2014 and 2015 in study region. This class of drought was also found in other regions in years 2014 to 2016, as shown by studies of Uttaruk and Laosuwan (2017UTTARUK, Y.; LAOSUWAN, T. Drought detection by application of remote sensing technology and vegetation phenology. Journal of Ecological Engineering, v. 18, n. 6, p. 115-121, 2017. https://doi.org/10.12911/22998993/76326
https://doi.org/10.12911/22998993/76326...
) who related a reduction in VCI to an increase in extent of drought areas in a given region.

Most of study region showed VCI values in No Drought or Mild Drought classes (Figure 8). The areas in Extreme Drought class, with VCI values less than 20%, appeared in summer season, between 2014 and 2016, coinciding with period of water crisis that occurred in sub-basin, reported in Minas Gerais (2015), mainly in areas close to water supply reservoir. The spring season showed regions of Extreme Drought class only in years 2013 to 2014. In autumn season, the presence of drought areas was identified between 2013 and 2017 and the same situation occurred in winter season.

Figure 8
Area extensions in ha of drought classes provided by VCI for each season from year 2013 to 2018

With the application of VCI, it was possible to identify drought variation in the Serra Azul reservoir (Figure 9). According to Florenzano (2013FLORENZANO, T. G. Iniciação em sensoriamento remoto. 3o edição ed. São Paulo: Oficina de Textos, 2013.), Formaggio and Sanches (2017) and Ponzoni, Shimabukuro and Kuplich (2015PONZONI, F. J.; SHIMABUKURO, Y. E.; KUPLICH, T. M. Sensoriamento remoto da vegetação. São José dos Campos: Oficina de Textos, 2015.), as the spectral response of water is low in red and infrared bands, in VCI calculation the water regions tend to have a lower value, being associated with regions of Extreme Drought class. In this regard, the results indicated that as areas of Extreme Drought class became more evident, the tendency of water body areas to decrease also increased. As seen in Figure 8, the large concentration of Extreme Drought areas in 2013 favored a decrease in volume of water in reservoir, a fact that was documented for years 2014 and 2015 at the time of water scarcity episode suffered by region (MINAS GERAIS, 2015). In the same way, the increase in areas of Moderate and Mild Drought in reservoir region, between 2016 and 2017, contributed to increase in this area in 2018.

The research shows that changes in land use and natural changes have led to worsening periods of drought in the sub-basin, especially in years of water scarcity, as reported by Minas Gerais (2015). According to Dutra (2021DUTRA, D. J. Uso so sensoriamento remoto para análise de eventos de seca em bacias hidrográficas: estudo de caso na sub-bacia do ribeirão Serra Azul. MG. Dissertação de mestrado ed. Belo Horizonte: Universidade Federal de Minas Gerais (Dissertação). Programa de Pós-graduação em Análise e Modelagem de Sistemas Ambientais da Universidade Federal de Minas Gerais, 2021.), the influence of land use variations and natural events on the worsening of the drought phenomenon in the study region was identified based on water balance calculations in the sub-basin. In this regard, the changes associated to urban and agricultural expansion in region caused soil impermeability and a consequent decrease of basin maximum storage volume. Thus, in years when there is a decrease of water entry in system, due to movement of air masses that provide low precipitation, the region presented a decrease in water volume and increase of the worsening drought classes. This causes water scarcity episodes and difficulty in supplying water to the population due to the decrease in the Serra Azul reservoir volume.

According to Kamble et al. (2019KAMBLE, D. B. et al. Drought assessment for kharif rice using standardized precipitation index (SPI) and vegetation condition index (VCI). Journal of Agrometeorology, v. 21, n. 2, p. 182-187, 2019.) and Eyoh, Okeke and Ekpa (2019EYOH, A.; OKEKE, F.; EKPA, A. Assessment of the effectiveness of the Vegetation Condition Index (VCI) as an indicator for monitoring drought condition across the Niger Delta region of Nigeria using AVHRR / MODIS NDVI. European Journal of Earth and Environment, v. 6, n. 1, p. 12-18, 2019.), Yulistya, Wibowo and Kusratmoko (2019YULISTYA, V. D.; WIBOWO, A.; KUSRATMOKO, E. Assessment of agricultural drought in paddy field area using Vegetation Condition Index (VCI) in Sukaresmi District, Cianjur Regency. IOP Conference Series: Earth and Environmental Science, v. 311, n. 1, 2019. https://doi.org/10.1088/1755-1315/311/1/012020
https://doi.org/10.1088/1755-1315/311/1/...
) and Baniya et al. (2019BANIYA, B. et al. Spatial and temporal variation of drought based on satellite derived vegetation condition index in Nepal from 1982-2015. Sensors (Switzerland), v. 19, n. 2, 2019. https:/doi.org/10.3390/s19020430
https:/doi.org/10.3390/s19020430...
), the variation of VCI classes occurs due to the ability of this index to demonstrate the worst and the best drought condition over a time period due to the variation in precipitation and temperature in the region. The rains have a strong interaction with vegetation, being able to directly influence in vegetative vigor of a region. Thus, the greater the storage of water in vegetation structure, greater the VCI value and less the chance of this vegetation suffering from drought events.

Figure 9
Classification of VCI for the Serra Azul reservoir region in dry seasons (winter and autumn)

The results obtained allowed to identify drought events extent and influence of seasonality over seasons in the sub-basin through the application of VCI. It was observed that region presented periods of extreme drought, with a worsening between 2014 and 2016, allowing identifying the water scarcity process that occurred in study region since regions of extreme and severe drought were located close to the regions of reservoir and without vegetation presence in region.

Comparative analysis of NDVI, EVI and VCI indexes for drought monitoring in watershed regions

The results indicated that the VCI showed a better correlation with meteorological variables, between moderately and strongly positive, when compared to data of vegetation indexes (Table 5). This correlation has been corroborated in several international studies such as those of Zambrano et al. (2016ZAMBRANO, F. et al. Sixteen years of agricultural drought assessment of the biobío region in chile using a 250 m resolution vegetation condition index (VCI). Remote Sensing, v. 8, n. 6, p. 1-20, 2016. https://doi.org/10.3390/rs8060530
https://doi.org/10.3390/rs8060530...
) and Gomes et al. (2019GOMES, A. R. S. et al. Estudo da Relação entre a Variabilidade dos índices de Vegetação e Temperatura da Região Nordeste do Brasil. Revista Brasileira de Meteorologia, v. 34, n. 3, p. 359-368, 2019. https://doi.org/10.1590/0102-7786343051
https://doi.org/10.1590/0102-7786343051...
), that demonstrated the VCI has a strong correlation with meteorological data, being considered a very useful tool for analysis in regions where there is a worsening drought. According to Xu et al. (2020XU, Z. et al. Trends in Global Vegetative Drought from Long-Term Satellite Remote Sensing Data. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, v. 13, p. 815-826, 2020. https://doi.org/10.1109/JSTARS.2020.2972574
https://doi.org/10.1109/JSTARS.2020.2972...
) and Gu et al. (2019GU, L. et al. The contribution of internal climate variability to climate change impacts on droughts. Science of the Total Environment, v. 684, p. 229-246, 2019. https://doi.org/10.1016/j.scitotenv.2019.05.345
https://doi.org/10.1016/j.scitotenv.2019...
), the normalization of NDVI data in the VCI calculation process contributes to a moderate and strong correlation when compared with meteorological data. Normalization, allows the analysis of drought and seasonal trends presented in a region over the climatic variation.

Table 5
Pearson's correlation coefficient (r) between indexes and meteorological variables for period from 2013 to 2018.

The correlation of vegetation indexes, such as NDVI and EVI is not a straightforward function in relation to climatic variables, so arising discrepancies when correlated. The results of research identified that vegetation indexes presented a correlation between moderate and weak with climatological data. According to Cohen (1998COHEN, J. No Statistical Power Analysis for the Behavioral Sciences. Second ed. Mahwah: Lawrence Erlbaum, 1998.), Dancey and Reidy (2006DANCEY, C. P.; REIDY, J. Estatística sem Matemática para Psicologia: usando SPSS para Windows. Porto Alegre: Penso, 2006.), Paranhos et al. (2014PARANHOS, R. et al. Desvendando os Mistérios do Coeficiente de Correlação de Pearson: o Retorno. Leviathan (São Paulo), n. 8, p. 66, 2014. https://doi.org/10.11606/issn.2237-4485.lev.2014.132346
https://doi.org/10.11606/issn.2237-4485....
), Nora and Santos (2010NORA, E. L. D.; SANTOS, J. E. Análise da dinâmica sazonal de duas formaçoes florestais do bioma Mata Atlântica com base em índices de vegetação. Perspectiva, v. 34, n. 125, p. 41-51, 2010.) and Quesada et al. (2017QUESADA, H. B. et al. Análise da vegetação ripária em bacia hidrográfica utilizando Índice de Vegetação Normalizada (NDVI) no município de Maringá-PR. Geo UERJ, v. 0, n. 31, p. 439-455, 2017. https://doi.org/10.12957/geouerj.2017.26737
https://doi.org/10.12957/geouerj.2017.26...
), meteorological variables usually do not show strong correlation values (above 0.8) when compared to vegetation indexes. This is because plant areas only show physiological changes after effective occurrence of a certain environmental phenomenon. The changes are not manifested immediately during the occurrence of the event.

Results showed that EVI is more dependent on precipitation variable while the NDVI is more dependent on temperature. This relationship was corroborated in studies presented by Kafer and Rex (2020KÄFER, P. S.; REX, F. E. Avaliação espectral e temporal de remanescentes da mata atlântica com dados Spot-Vgt e variáveis meteorológicas. BIOFIX Scientific Journal, v. 5, n. 1, p. 13-22, 2020. https://doi.org/10.5380/biofix.v5i1.67235
https://doi.org/10.5380/biofix.v5i1.6723...
) and Nora and Santos (2010NORA, E. L. D.; SANTOS, J. E. Análise da dinâmica sazonal de duas formaçoes florestais do bioma Mata Atlântica com base em índices de vegetação. Perspectiva, v. 34, n. 125, p. 41-51, 2010.), who reported a better correlation between NDVI and air temperature than precipitation due to the fact that forest formations have resistance to short periods of water scarcity. According to Chaves, Mataveli and Justino (2014CHAVES, M. E. D.; MATAVELI, G. A. V.; JUSTINO, R. C. Uso da modelagem estatística para monitoramento da vegetação no Parque Nacional da Serra da Canastra, Minas Gerais. Caderno de Geografia, v. 24, n. 1, p. 120-132, 2014. https://doi.org/10.5752/P.2318-2962.2014v24nespp120
https://doi.org/10.5752/P.2318-2962.2014...
), as the vegetation develops, the standard deviation of vegetation index data in a time series increases. The lowest values of these variables tend to appear in periods of drought, especially in winter, between August and October, being related to lack of rain, which causes a decrease in values of vegetation indexes.

It was found in the results that VCI has a better interaction with climate data if compared with data on vegetation indexes, such as EVI and NDVI. The normalization of data performed in the process of VCI calculation allows a better identification of seasonal variations and influence of climatic variables on vegetation of a region.

FINAL CONSIDERATIONS

Through the statistical analysis adopted in research methodology, NDVI was identified as a more effective vegetation index than EVI in characterizing drought for transition regions from rural to urban space. The vegetation condition index (VCI) proved to be feasible and suitable for monitoring drought over study period between 2013 and 2018 having effectively delimited states of drought in study area, also making clear association of index with meteorological data. Finally, the widely available data from MODIS sensor that allowed reaching these results have been demonstrated as suitable and effective in characterizing these drought events investigated in study area.

ACKNOWLEGMENTS

This study was financed by the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES), a Brazilian government agency for improvement of higher education human resources - Finance Code 001.

REFERENCES

  • ARANTES, T. B. et al. Effectiveness of BFAST algorithm to characterize time series of dense forest, agriculture and pasture in the amazon region. Theoretical and applied engineering, v. 1, n. 1, p. 10-19, 2017.
  • BANIYA, B. et al. Spatial and temporal variation of drought based on satellite derived vegetation condition index in Nepal from 1982-2015. Sensors (Switzerland), v. 19, n. 2, 2019. https:/doi.org/10.3390/s19020430
    » https:/doi.org/10.3390/s19020430
  • BHUIYAN, C.; KOGAN, F. N. Monsoon variation and vegetative drought patterns in the Luni Basin in the rain-shadow zone. International Journal of Remote Sensing, v. 31, n. 12, p. 3223-3242, 2010. https://doi.org/10.1080/01431160903159332
    » https://doi.org/10.1080/01431160903159332
  • BONIFACIO, R.; DUGDALE, G.; MILFORD, J. R. Sahelian rangeland production in relation to rainfall estimates from Meteosat. International Journal of Remote Sensing, v. 14, n. 14, p. 2695-2711, 1993.
  • BRANCO, E. R. F. Ocorrências de seca e tendências da vegetação na reserva biológica de sooretama e zona de amortecimento, no estado do Espírito Santo, Brasil. Dissertação de mestrado ed. Jerônimo Monteiro: Universidade Federal do Espirito Santo, 2016.
  • CHAVES, M. E. D.; MATAVELI, G. A. V.; JUSTINO, R. C. Uso da modelagem estatística para monitoramento da vegetação no Parque Nacional da Serra da Canastra, Minas Gerais. Caderno de Geografia, v. 24, n. 1, p. 120-132, 2014. https://doi.org/10.5752/P.2318-2962.2014v24nespp120
    » https://doi.org/10.5752/P.2318-2962.2014v24nespp120
  • COHEN, J. No Statistical Power Analysis for the Behavioral Sciences. Second ed. Mahwah: Lawrence Erlbaum, 1998.
  • COVELE, P. A. Aplicação de índices das condições de vegetação no monitoramento em tempo quase real da seca em Moçambique usando NOAA-AVHRR. GEOUSP - Espaço e Tempo, n. 29, p. 85-95, 2011.
  • COX, D. R.; STUART, A. Some quick sign tests for trend in location and dispersion. Biometrika, v. 42, n. 1/2, p. 80, 1955.
  • CUNHA, A. P. M. A. et al. Avaliação de indicador para o monitoramento dos impacos da seca em áreas de pastagens no Semiárido do Brasil. Revista Brasileira de Cartografia, v. 69, n. 1, p. 89-106, 2017.
  • DANCEY, C. P.; REIDY, J. Estatística sem Matemática para Psicologia: usando SPSS para Windows. Porto Alegre: Penso, 2006.
  • DAVIES, N.; CHATFIELD, C. The analysis of time series: an introduction. 6. ed. New York Washington: CHAPMAN & HALL/CRC, 1990. v. 74. https://doi.org/10.2307/3619403
    » https://doi.org/10.2307/3619403
  • DECHANT, C. M.; MORADKHANI, H. Analyzing the sensitivity of drought recovery forecasts to land surface initial conditions. Journal of Hydrology, v. 526, p. 89-100, 2015. http://dx.doi.org/10.1016/j.jhydrol.2014.10.021
    » http://dx.doi.org/10.1016/j.jhydrol.2014.10.021
  • DECHANT, C. M.; MORADKHANI, H. Toward a reliable prediction of seasonal forecast uncertainty: Addressing model and initial condition uncertainty with ensemble data assimilation and Sequential Bayesian Combination. Journal of Hydrology, v. 519, n. PD, p. 2967-2977, 2014. http://dx.doi.org/10.1016/j.jhydrol.2014.05.045
    » http://dx.doi.org/10.1016/j.jhydrol.2014.05.045
  • DETZEL, D. H. M. et al. Estacionariedade das afluências às usinas hidrelétricas brasileiras. Revista Brasileira de Recursos Hídricos, v. 16, n. 3, p. 95-111, 2011. https://doi.org/10.21168/rbrh.v16n3.p95-111
    » https://doi.org/10.21168/rbrh.v16n3.p95-111
  • DU, L. et al. A comprehensive drought monitoring method integrating MODIS and TRMM data. International Journal of Applied Earth Observation and Geoinformation, v. 23, n. 1, p. 245-253, 2013. http://dx.doi.org/10.1016/j.jag.2012.09.010
    » http://dx.doi.org/10.1016/j.jag.2012.09.010
  • DUBREUIL, V. et al. Os tipos de climas anuais no Brasil: Uma aplicação da classificação de Köppen de 1961 a 2015. Confins. Revue franco- brésilienne de géographie/Revista franco-brasilera de geografia, v. 37, 2018.
  • DUFT, D. G.; PICOLI, M. C. A. Uso de imagens do sensor MODIS para identificação da seca na cana-de-açúcar através de índices espectrais. Scientia Agraria, v. 19, n. 1, p. 52, 2018. https://doi.org/10.5380/rsa.v19i1.54005
    » https://doi.org/10.5380/rsa.v19i1.54005
  • DUTRA, D. J. Uso so sensoriamento remoto para análise de eventos de seca em bacias hidrográficas: estudo de caso na sub-bacia do ribeirão Serra Azul. MG. Dissertação de mestrado ed. Belo Horizonte: Universidade Federal de Minas Gerais (Dissertação). Programa de Pós-graduação em Análise e Modelagem de Sistemas Ambientais da Universidade Federal de Minas Gerais, 2021.
  • DUTRA, D. J.; BRIANEZI, D.; COELHO, C. W. G. A. Uso de Geotecnologias para Análise da Dinâmica da Vegetação da Sub-bacia do Ribeirão Serra Azul, MG. Anuário do Instituto de Geociências - UFRJ, v. 43, n. 4, p. 283-292, 18 dez. 2020. http://dx.doi.org/10.11137/2020_4_283_292
    » http://dx.doi.org/10.11137/2020_4_283_292
  • DUTRA, D. J.; ELMIRO, M. A. T.; GARCIA, R. A. Comparative analysis of methods applied in vegetation cover delimitation using Landsat 8 images. Sociedade & Natureza, v. 32, n. July, p. 699-710, 9 out. 2020. https://doi.org/10.14393/SN-v32-2020-56139
    » https://doi.org/10.14393/SN-v32-2020-56139
  • EYOH, A.; OKEKE, F.; EKPA, A. Assessment of the effectiveness of the Vegetation Condition Index (VCI) as an indicator for monitoring drought condition across the Niger Delta region of Nigeria using AVHRR / MODIS NDVI. European Journal of Earth and Environment, v. 6, n. 1, p. 12-18, 2019.
  • FIGUEIREDO FILHO, D. B.; SILVA JÚNIOR, J. A. Desvendando os mistérios do coeficiente de correlação de Pearson (r). Revista Política Hoje, v. 18, n. 1, p. 115-146, 2009.
  • FLORENZANO, T. G. Iniciação em sensoriamento remoto. 3o edição ed. São Paulo: Oficina de Textos, 2013.
  • FORMAGGIO, A. R.; SANCHES, I. D. Índices espectrais de vegetação x agricultura. Sensoriamento remoto em agricultura. 1. ed. São Paulo: Oficina de Textos, 2017a. p. 95-113.
  • FORMAGGIO, A. R.; SANCHES, I. D. Sensoriamento remoto em agricultura. 1o ed. São José dos Campos: Oficina de Textos, 2017b.
  • GOMES, A. R. S. et al. Estudo da Relação entre a Variabilidade dos índices de Vegetação e Temperatura da Região Nordeste do Brasil. Revista Brasileira de Meteorologia, v. 34, n. 3, p. 359-368, 2019. https://doi.org/10.1590/0102-7786343051
    » https://doi.org/10.1590/0102-7786343051
  • GOW, L. J. et al. A detection problem: Sensitivity and uncertainty analysis of a land surface temperature approach to detecting dynamics of water use by groundwater-dependent vegetation. Environmental Modelling and Software, v. 85, p. 342-355, 2016. https://doi.org/10.1016/j.envsoft.2016.09.003
    » https://doi.org/10.1016/j.envsoft.2016.09.003
  • GU, L. et al. The contribution of internal climate variability to climate change impacts on droughts. Science of the Total Environment, v. 684, p. 229-246, 2019. https://doi.org/10.1016/j.scitotenv.2019.05.345
    » https://doi.org/10.1016/j.scitotenv.2019.05.345
  • HUETE, A. et al. Overview of the radiometric and biophysical performance of the MODIS vegetation indices. Remote Sensisng of Enviroment, v. 83, p. 195-213, 2002.
  • IBGE. Dados vetoriais. Avaible in: https://www.ibge.gov.br/geociencias-novoportal/informacoes-ambientais/climatolo gia /15817-clima.html?=&t=downloads Access in: jul, 2020.
    » https://www.ibge.gov.br/geociencias-novoportal/informacoes-ambientais/climatolo gia /15817-clima.html?=&t=downloads
  • JIAO, W. et al. A new station-enabled multi-sensor integrated index for drought monitoring. Journal of Hydrology, v. 574, n. April, p. 169-180, 2019. https://doi.org/10.1016/j.jhydrol.2019.04.037
    » https://doi.org/10.1016/j.jhydrol.2019.04.037
  • KÄFER, P. S.; REX, F. E. Avaliação espectral e temporal de remanescentes da mata atlântica com dados Spot-Vgt e variáveis meteorológicas. BIOFIX Scientific Journal, v. 5, n. 1, p. 13-22, 2020. https://doi.org/10.5380/biofix.v5i1.67235
    » https://doi.org/10.5380/biofix.v5i1.67235
  • KAMBLE, D. B. et al. Drought assessment for kharif rice using standardized precipitation index (SPI) and vegetation condition index (VCI). Journal of Agrometeorology, v. 21, n. 2, p. 182-187, 2019.
  • KOGAN, F. N. Application of vegetation index and brightness temperature for drought detection. Advances in Space Research, v. 15, n. 11, p. 91-100, 1995. https://doi.org/10.1016/0273-1177(95)00079-T
    » https://doi.org/10.1016/0273-1177(95)00079-T
  • LEIVAS, J. F. et al. Monitoramento da seca 2011/2012 no nordeste brasileiro a partir do satélite Spot-vegetation e TRMM. Revista Engenharia Na Agricultura - Reveng, v. 22, n. 3, p. 211-221, 2014. https://doi.org/10.13083/1414-3984.v22n03a03
    » https://doi.org/10.13083/1414-3984.v22n03a03
  • MINAS GERAIS. Portaria IGAM no 014, de 08 de abril de 2015. Declara situação crítica de escassez hídrica superficial na porção hidrográfica localizada no reservatório Serra Azul e a sua bacia de contribuição.: Belo Horizonte. , 2015
  • MORETTIN, P. A.; TOLOI, C. M. C. Análise de séries temporais. 4. ed. New Jersey: Blucher, 2006.
  • MORETTIN, P. A.; TOLOI, C. M. C. Modelos para previsão de séries temporais. Rio de Janeiro: Instituto de Matemática Pura e Aplicada, 1981.
  • NORA, E. L. D.; SANTOS, J. E. Análise da dinâmica sazonal de duas formaçoes florestais do bioma Mata Atlântica com base em índices de vegetação. Perspectiva, v. 34, n. 125, p. 41-51, 2010.
  • PARANHOS, R. et al. Desvendando os Mistérios do Coeficiente de Correlação de Pearson: o Retorno. Leviathan (São Paulo), n. 8, p. 66, 2014. https://doi.org/10.11606/issn.2237-4485.lev.2014.132346
    » https://doi.org/10.11606/issn.2237-4485.lev.2014.132346
  • PEARSON, K. The grammar of science. Nature, v. 46, n. 1185, p. 247-247, 1982. https://doi.org/10.1038/046247b0
    » https://doi.org/10.1038/046247b0
  • PEARSON, K.; FISCHER, R.; INMAN, H. F. Karl Pearson and R. A. Fischer on Statistical Test: A 1935 exchange from nature. The American Statistician, v. 48, n. 1, p. 2-11, 1994. https://doi.org/10.1080/00031305.1994.10476010
    » https://doi.org/10.1080/00031305.1994.10476010
  • PONZONI, F. J.; SHIMABUKURO, Y. E.; KUPLICH, T. M. Sensoriamento remoto da vegetação. São José dos Campos: Oficina de Textos, 2015.
  • POTTER, C. et al. Major disturbance events in terrestrial ecosystems detected using global satellite data sets. Global Change Biology, v. 9, n. 7, p. 1005-1021, 2003. https://doi.org/10.1046/j.1365-2486.2003.00648.x
    » https://doi.org/10.1046/j.1365-2486.2003.00648.x
  • QUESADA, H. B. et al. Análise da vegetação ripária em bacia hidrográfica utilizando Índice de Vegetação Normalizada (NDVI) no município de Maringá-PR. Geo UERJ, v. 0, n. 31, p. 439-455, 2017. https://doi.org/10.12957/geouerj.2017.26737
    » https://doi.org/10.12957/geouerj.2017.26737
  • SAUSEN, T. M.; LACRUZ, M. S. P. Sensoriamento remoto para desastres. 1. ed. São Paulo: Oficina de Textos, 2015.
  • SOARES, A. Geoestatística para as ciências da terra e do ambiente. Rio de Janeiro: Instituto Superior Técnico, 2000.
  • UTTARUK, Y.; LAOSUWAN, T. Drought detection by application of remote sensing technology and vegetation phenology. Journal of Ecological Engineering, v. 18, n. 6, p. 115-121, 2017. https://doi.org/10.12911/22998993/76326
    » https://doi.org/10.12911/22998993/76326
  • VIEIRA, S. et al. Geoestatistical theory and application to variability of some agronomical properties. Hilgardia, v. 51, n. 3, p. 1-75, 1983. https://doi.org/10.3733/hilg.v51n03p075
    » https://doi.org/10.3733/hilg.v51n03p075
  • XU, Z. et al. Trends in Global Vegetative Drought from Long-Term Satellite Remote Sensing Data. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, v. 13, p. 815-826, 2020. https://doi.org/10.1109/JSTARS.2020.2972574
    » https://doi.org/10.1109/JSTARS.2020.2972574
  • YULISTYA, V. D.; WIBOWO, A.; KUSRATMOKO, E. Assessment of agricultural drought in paddy field area using Vegetation Condition Index (VCI) in Sukaresmi District, Cianjur Regency. IOP Conference Series: Earth and Environmental Science, v. 311, n. 1, 2019. https://doi.org/10.1088/1755-1315/311/1/012020
    » https://doi.org/10.1088/1755-1315/311/1/012020
  • ZAMBRANO, F. et al. Sixteen years of agricultural drought assessment of the biobío region in chile using a 250 m resolution vegetation condition index (VCI). Remote Sensing, v. 8, n. 6, p. 1-20, 2016. https://doi.org/10.3390/rs8060530
    » https://doi.org/10.3390/rs8060530
  • ZHANG, L. et al. Studying drought phenomena in the Continental United States in 2011 and 2012 using various drought indices. Remote Sensing of Environment, v. 190, p. 96-106, 2017. http://dx.doi.org/10.1016/j.rse.2016.12.010
    » http://dx.doi.org/10.1016/j.rse.2016.12.010

Publication Dates

  • Publication in this collection
    01 Sept 2021
  • Date of issue
    2021

History

  • Received
    26 Feb 2021
  • Accepted
    01 June 2021
  • Published
    09 June 2021
Editora da Universidade Federal de Uberlândia - EDUFU Av. João Naves de Ávila, 2121 - Bloco 5M – Sala 302B, 38400902 - Uberlândia - Minas Gerais - Brasil, +55 (34) 3239- 4549 - Uberlândia - MG - Brazil
E-mail: sociedade.natureza@ig.ufu.br