Next Article in Journal
Ambient-Visible-Light-Mediated Enhanced Degradation of UV Stabilizer Bis(4-hydroxyphenyl)methanone by Nanosheet-Assembled Cobalt Titanium Oxide: A Comparative and DFT-Assisted Investigation
Next Article in Special Issue
The Sensitivity of Meteorological Dynamics to the Variability in Catchment Characteristics
Previous Article in Journal
Budyko-Type Models and the Proportionality Hypothesis in Long-Term Water and Energy Balances
Previous Article in Special Issue
Hydrological Response to Meteorological Droughts in the Guadalquivir River Basin, Southern Iberian Peninsula
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Developing a Combined Drought Index to Monitor Agricultural Drought in Sri Lanka

1
Department of Ecology and Conservation Biology, Texas A&M University, College Station, TX 77843, USA
2
The World Bank, Washington, DC 20433, USA
3
Department of Agriculture, Peradeniya 20400, Sri Lanka
4
Department of Civil Engineering, University of Peradeniya, Natural Resource Management Centre, Peradeniya 20400, Sri Lanka
*
Author to whom correspondence should be addressed.
Water 2022, 14(20), 3317; https://doi.org/10.3390/w14203317
Submission received: 18 September 2022 / Revised: 13 October 2022 / Accepted: 18 October 2022 / Published: 20 October 2022

Abstract

:
Developing an agricultural drought monitoring index through integrating multiple input variables into a single index is vital to facilitate the decision-making process. This study aims to develop an agricultural drought index (agCDI) to monitor and characterize the spatial and temporal patterns of drought in Sri Lanka. Long-term (1982 to 2020) remote sensing and model-based agroclimatic input parameters—normalized difference vegetation index (NDVI), land surface temperature (LST), 3-month precipitation z-score (stdPCP), and evaporative demand drought index (EDDI)—were used to develop agCDI. The principal component analysis (PCA) approach was employed to qualitatively determine the grid-based percentage contribution of each input parameter. The agCDI was apparently evaluated using an independent dataset, including the crop yield for the major crop growing districts and observed streamflow-based surface runoff index (SRI) for the two main crop growing seasons locally, called Yala (April to September) and Maha (October to March), using 20-years of data (from 2000 to 2020). The results illustrate the good performance of agCDI, in terms of predominantly capturing and characterizing the historic drought conditions in the main agricultural producing districts both during the Yala and Maha seasons. There is a relatively higher chance of the occurrence of moderate to extreme droughts in the Yala season, compared to the Maha season. The result further depicts that relatively good correlation coefficient values (> 0.6) were obtained when agCDI was evaluated using a rice crop yield in the selected districts. Although the agCDI correlated well with SRI in some of the stations (>0.6), its performance was somehow underestimated in some of the stations, perhaps due to the time lag of the streamflow response to drought. In general, agCDI showed its good performance in capturing the spatial and temporal patterns of the historic drought and, hence, the model can be used to develop agricultural drought monitoring and an early warning system to mitigate the adverse impacts of drought in Sri Lanka.

1. Introduction

Drought is a recurring natural hazard in Sri Lanka that triggers water shortage that consequently has a diversified impact on agricultural production and other sectors [1]. When severe to extreme drought events occur, Sri Lanka suffers economic losses and damages, including crop failure, lack of clean and adequate quantity of drinking water, inadequate hydropower generation, and disease outbreaks [2]. The primary cause of drought is a persistent precipitation deficit or failure of monsoon rainfall that potentially prompts the subsequent decline in soil moisture, groundwater, and other surface water storage and flows accessible for agriculture, domestic, and other uses [3,4]. Previous studies have illustrated the frequent occurrence of drought in Sri Lanka [5,6,7]. For example, Nianthi [5] reported that the frequency of occurrence of extreme drought in Sri Lanka is once every ten years, even though local to regional droughts occur once every 3 to 4 years. Although there are several contributing factors to the occurrences and severity of drought, the failure of the monsoons and effects of El Niño–Southern Oscillation are the primary causes that prompt droughts in Sri Lanka. For example, the failure of monsoons led to the 2007 drought in parts of Sri Lanka, and El Niño–Southern Oscillation events were connected to the 1982-1983 drought event, where more than one million people were affected [8]. Rainfed agriculture is the first sector vulnerable to drought, and its failure often significantly influences the agrarian economy, the livelihood of the community, and the household level income [9].
According to a World Bank report, the agricultural sector employed approximately 24% of the total labor force in Sri Lanka in 2020, even though its contribution to the gross domestic product (GDP) has been declining since 1974 [10]. Water resources, including the drinking water supply, are also vulnerable to drought and climate change and sensitive to the seasonal variability of the rainfall [11]. There were significant agricultural losses in 1988, 2001, and 2004 largely attributed to the occurrence of severe drought events [12]. Drought-related impacts surpass other natural disasters (e.g., floods, landslides, etc.) and have accounted for the highest percentage (~52%) of the agricultural losses from 1974 to 2013 [2]. In terms of the total number of populations affected by the historical drought events, more than one million farmers were affected by drought in 2013 [2] and 2016-2017 [1,13]. This highlights the need for the development of drought-resilient farming systems to mitigate its adverse impacts on socioeconomic sectors through developing effective drought monitoring and early warning systems.
Wilhite and Glantz [14] categorized drought into four types, including meteorological, agricultural, hydrological, and socioeconomic droughts, although their definitions vary per discipline. Meteorological drought measures the degree of dryness due to persistent deficits of rainfall, in comparison to the threshold value. Agricultural drought, on the other hand, measures the deficits of soil moisture and evapotranspiration, largely due to the shortage of rainfall, while hydrological drought characterizes the effects of precipitation shortfalls on the groundwater, streamflow, reservoir, and lake levels. Socioeconomic drought is linked to drought impacts on the socioeconomic sectors. These droughts can be characterized using the drought monitoring indices. There have been significant efforts made to develop drought monitoring systems at the local to global spatial scales using drought indices [15,16,17,18]. The effectiveness of these drought indices, however, depends largely on the availability of input data at reasonable spatial and temporal scales, beyond the difference in the conceptual frameworks. Most of the drought indices use a single input variable, traditionally obtained from local gauging stations, to characterize the different aspects of drought. Despite the good quality and accuracy of the station’s data, the lack of capturing the continuous drought condition across space, due to their sparse locations, limits the wide application to characterize drought at a larger spatial scale, particularly in developing countries [19]. Therefore, remote sensing products are extensively applied to bridge these gaps and to monitor drought at a finer spatial resolution covering larger extents [20]. Some of the remote sensing products even blend station-based observations to enhance the data quality and accuracy. In addition to meteorological observations, the remote sensing-based approaches make other key parameters paramount to monitoring droughts available, such as the vegetation health, soil moisture, and water levels [21]. The potential application of remote sensing products has been highlighted in water resource management, particularly in data-scarce regions in developing countries [22,23].
In Sri Lanka, there are several drought monitoring studies that have been conducted to characterize meteorological and agricultural droughts using the most common and widely used drought indices. For example, Alahacoon et al. [1] attempted to monitor meteorological and agricultural droughts using the standardized precipitation index (SPI), rainfall anomaly index (RAI), and vegetation health index (VHI) drought indices. There were also other studies conducted to monitor meteorological drought [4,13,24,25,26,27,28], agricultural drought [1,29,30,31,32], hydrological drought [25], and socioeconomic drought [33]. However, there is still limited effort to evaluate the impacts of drought on the agricultural sector using the integrated drought approaches by combining several hydroclimatic and agricultural input variables.
Since droughts are evaluated based on the local impacts, the broad spectrum of drought in a given sector cannot be fully captured using a single index that uses largely single input variables [34]. This prompted the great attention of researchers towards developing the combined drought index using several biophysical input variables under a “convergence of evidence” framework [35]. Thus far, few attempts have been made to monitor drought by developing the combined drought index, and the effectiveness of the index has been reported to monitor drought in real-time across the US and Canada [36,37,38], East Africa [39,40], the Middle East and North Africa regions [34], and India [41]. Different approaches, such as regression tree, PCA, and expert opinion and judgment, were used to quantify the percentage contributions of the input variables. Such drought monitoring approaches have never been used to develop a combined drought index in Sri Lanka, as per the knowledge of the authors. Therefore, this study aimed to develop a combined drought index to monitor agricultural drought in Sri Lanka using the remote sensing-based input variables, including CHIRPS rainfall, NDVI, LST, and EDDI. An objective approach was followed to estimate the grid-based percentage contributions of each input variable in the process of developing the combined drought indicators. The method followed in this study can be applied to other study regions to develop a drought monitoring tool that can be used to characterize the severity of the historic drought conditions.

2. Materials and Methods

2.1. Description of the Study Area

Sri Lanka is an island country located at the southern tip of the Indian sub-continent at a geographic coordinate of 7°00 N and 81°00E. It covers a boundary area of 65,610 sq km, and the elevation ranges between 0 m (Indian Ocean) to 2524 m (Pidurutalagala peak in the central massif). According to the Department of Meteorology of Sri Lanka, there are four rainfall seasons, including northeast monsoon (NEM, December–February), first inter-monsoon (FIM; March–April), southwest monsoon (SWM; May–September), and the second inter-monsoon (SIM; October–November). The rainfall originates mainly from monsoonal, conventional, and weather systems formed in the Bay of Bengal (i.e., low-level atmospheric disturbances and depressions) and is spatially variable, due to orographic influences [42]. Figure 1 shows the three climatic zones in Sri Lanka (wet, intermediate, and dry zones), delineated based on hydroclimatic variables, mainly rainfall and other biophysical characteristics. The wet zone receives relatively high mean annual precipitation (>2500 mm) during the SWM and SIM seasons, as compared to the dry zone (less than 1750 mm during the SIM and NEM and the intermediate zone between 1750 and 2500 mm). The temperature reaches a maximum during summer months, with the highest values ranging between 32 °C and 35 °C in April and May, while the lowest temperature ranges between 18 °C and 24 °C in January [43]. The temperature further declines to the range of 9 °C and 21 °C in the central highlands throughout the year. Maha (October to March) and Yala (April to the end of September) are the two main crop-growing seasons [1,44,45]. Maha growing season covers the two consecutive rainfall seasons, SMI and NEM, while the Yala season combines FIM and SWM seasons [46]. The dry zone receives only FIM rainfall from mid-March to early May, since SWM rainfall is not onset in the dry zone. In addition, the dry zone receives only a small amount of rain for two months during the Yala season; thus, there is no significant agricultural production during this season, and supplementary irrigation is required to boost crop cultivation. Rice is the primary food crop grown in Sri Lanka. The per capita arable land is decreasing, due to the increase in population, despite agricultural land being the dominant land use (50%) in Sri Lanka [45].

2.2. Data Used

Table 1 shows the summary of input variables used in this study. These input variables represent the different components of the water budget and are viable for drought monitoring and evaluation. A detailed description of each input variable is presented in the subsequent subsections.

2.2.1. Precipitation

CHIRPS is a blended product developed by the U.S. Geological Survey (USGS) and the Climate Hazards Group at the University of California, Santa Barbara (UCSB). It combines precipitation climatology and quasi-global geostationary TIR satellite observations from the Climate Prediction Center (CPC) and the National Climate Forecast System, version 2 (CFSv2) [47]. Good performance and skillful application of this product have been validated in Sri Lanka [48,49], as well as in other parts of the world [50,51,52,53,54]. The long-term gridded data were acquired from the Climate Hazards Center data portal from 1982 to 2020, at 5 km × 5 km pixel size.

2.2.2. Normalized Difference Vegetation Index (NDVI) and Land Surface Temperature (LST)

Gridded smoothed weekly time series of NDVI and LST were acquired from National Oceanic and Atmospheric Administration (NOAA) Satellite Applications and Research (STAR) from 1982 to 2020. The data are a blended product of Visible Infrared Imaging Radiometer Suite (VIIRS, 2013–present) and Advanced Very High-Resolution Radiometer (AVHRR, 1982–2012), freely accessible from the NOAA data portal. The spatial resolution of the data is 4 km × 4 km, and the time series data were extracted for Sri Lanka from the global product. These products are widely applied to monitor the vegetation condition in response to changing weather and climate [55].

2.2.3. Evaporative Demand Drought Index (EDDI)

The EDDI is a drought index that characterizes the anomalous atmospheric evaporative demand across space and time [56]. EDDI can be applied to monitor agricultural and hydrologic droughts. It provides near-real-time information on the persistence of water stress and excess relevance for the preparedness of natural hazards. Penman–Monteith FAO56 was used to estimate the reference evapotranspiration during EDDI computation. Its input parameters, including temperature, humidity, wind speed, and incoming solar radiation, were extracted from North American Land Data Assimilation System (NLDAS-2). The spatial and temporal resolutions of EDDI are presented in Table 1.

2.2.4. Agricultural Land Use

The land use shapefile, which shows the distributions and extents of the different land use classes in Sri Lanka, was acquired from the Natural Resource Management Center, Department of Agriculture, Peradeniya. Figure 2 shows mainly the agricultural land use classes. Forest, grassland, bare area, urban, and rocks were aggregated into the “other land use” class, while open water and wetlands were grouped into water bodies. The percentage area coverage of each land use type is shown in Table 2. Woody perennial crops are the dominant agricultural land use (17%), widely distributed in different districts. Paddy crop is cultivated mainly in the dry zone, while perennial crops dominate in the wet and intermediate zones. The land use map masks the non-agricultural pixels, while developing the drought model and presenting the result.

2.2.5. Crop Yield

Crop yield data were acquired from the Department of Census and Statistics of Sri Lanka and used to evaluate the agricultural drought index. The crop yields are at level 2 administration districts spatial scale and are accessible for cereal crops, vegetable, and paddy crops from 2000 to 2020. However, the rice crop yield is used in this study for selected five districts responsible for more than half of the country’s annual rice production. The crop yield data are at seasonal temporal scale for each district; hence, the data for Yala and Maha seasons were used in this study.

2.2.6. Streamflow Data

The streamflow data were acquired from the Irrigation Department of Sri Lanka for stations mainly located in the wet and intermediate zones. The data are at the daily time step and for a total of 20 gauging stations. However, the data of ten stations are used in this study from 2000 to 2020, based on the availability of continuous historical time series of records.

2.3. Method

This method is designed to develop a combined drought index to monitor agricultural drought across Sri Lanka. Four remote sensing-based input variables, including NDVI, LST, precipitation, and EDDI, were used in this study. These input variables were selected due to their usefulness and broader application to monitoring agricultural drought and due to the availability of long-term data, as per the recommendation of previous studies [57,58,59]. The inconsistency in the spatial resolution of the input variables was treated using inverse distance weight (IDW) methods, and coarser spatial resolution input variables (e.g., EDDI) were resampled to 4km. Except for EDDI, all other input variables are available at a weekly time step on a global scale. Thus, the monthly EDDI, currently available in the public domain for the study country, was used to represent each week within the month. Next, each input variable, except EDDI, was standardized to avoid any impacts on the units of measurement and to enhance the spatial comparison of each parameter [38]. Three-month aggregate period of precipitation, recommended for agricultural drought monitoring, was considered to derive the time series z-score values (stdPCP). Then, we applied the principal component analysis (PCA) to compute the percentage contribution of each input variable.
Although PCA is widely applied in the atmospheric and hydrologic sciences to reduce the dimensionality of observational data, it has also become popular in drought studies to integrate several hydro-climatic variables into a single drought index [38]. It follows an objective approach to quantify the percentage contribution of each input variable for each week and grid point. The number of PCs corresponds to the number of input variables, and PCs cannot be combined to accommodate the total variability in the data [38]. However, the first principal axis accommodated the immense variability of the data (>75%) and was used to develop the percentage contributions of each input variable [39,41]. The eigenvalues of the first PC were estimated using the PCA approach and then squared to quantify the percentage contribution of each input variable [40]. Figure 3 illustrates the steps to develop PCA-based agriculture drought monitoring (agCDI) for Sri Lanka. There are six main steps, including:
  • Standardize all the input variables based on long-term weekly mean and standardization values.
  • Estimate the eigenvectors/eigenvalues used to transform each variable into separate orthogonal principal components (PCs).
  • Determine the percentage contributions (weight) of each input variable in PC1 (a total of 52 grid-based maps of the percentage contribution of each input variable were developed and used as a weight to combine the input variables into a single combined drought index).
  • Compute the weighted sum of the input variables using Equation (1) and generate the time series CDI maps, then normalize to minimize the higher degree of variability if it exists in some weeks.
  • Assess the spatial and temporal patterns of drought.
  • Evaluate the agricultural drought maps using independent datasets, including crop yield and streamflow-based SRI.
Equation (1) shows the mathematical representations of the agricultural drought indices, respectively.
a g C D I Z , I = ω p c p , i s t d P C P z , i + ω B T , i s t d L S T z , i + ω N D V I , i s t d N D V I z , i +   ω E D D I , i E D D I z , i
where z and I represent the data length (1982–2020) and weeks (52 weeks), respectively; a g C D I Z , I is weekly agricultural combined drought index value for a particular year ω p c p , i ,   ω B T , i ,   ω N D V I , i ,   a n d   ω E D D I , i are the percentage contribution of the anomaly values of precipitation, LST, NDVI, and EDDI, respectively.
stdPCP is standardized precipitation; stdLST is standardized land surface temperature; stdNDVI is standardized NDVI; and EDDI is the anomaly of the atmospheric evaporation demand.

2.4. Spatial and Temporal Assessment of Drought

The time series of agCDI were extracted at level 2 administration boundaries of Sri Lanka to analyze the temporal patterns of drought at the district level and to compare with district level detrended crop yield. McKee et. al. [60] drought classification category (Table 3) was used to identify the drought severity classes of the agCDI. Bayissa et al. [39] developed the threshold of the combined drought index and reported similar severity classes as McKee et al. [60]. Thus, in this study, we used the drought severity classes in Table 3 to categorize agCDI into different severity classes. The spatial patterns of drought were also assessed in agricultural-dominated land use in Sri Lanka. The frequency of occurrence of mild, moderate, severe, and extreme drought was computed using the time series of agCDI at each grid location. The percentage of the number of occurrences of each drought category (i.e., mild, moderate, severe, and extreme), with respect to the total data length, was used to quantify the frequency of the occurrence of each drought category. These drought characteristics measure the severity and persistence of the agricultural and water supply droughts.

2.5. Evaluation of agCDI

The agCDI was evaluated using two independent datasets (i.e., crop yield and streamflow-based SRI) for selected districts and hydrological boundaries. According to the Department of Census and Statistics (DCS) of Sri Lanka, Ampara, Polonnaruwa, Kurunegala, Anuradhapura, and Hambantota (Figure 3) districts are the central rice-producing districts in the country. More than 50% of the annual rice production occurs in these five districts. Hence, agCDI was evaluated using these five districts’ rice crop yield data for the Yala and Maha seasons. The crop yield data were detrended first to eliminate the impacts of improvement in agricultural management practices on the crop yield. Detrending removes the upward linear increase in the crop yield, due to improvements in farming practices, and accounts only for weather-related changes in crop yield [40,61]. The average seasonal values of agCDI were extracted at these five districts and compared with the corresponding detrended crop yield to explore any temporal association between these variables.
Similarly, we use the streamflow data of selected stations located mainly in the wet and intermediate zones to generate the time series of SRI. These stations were selected based on the availability of continuous flow data from 2000 to 2020. Then, the hydrological boundary of each flow gauging station was delineated and used to extract the seasonal time series of agCDI. The locations and extents of the hydrological boundaries are shown in Figure 4. The name of each hydrological boundary corresponds to the name of the corresponding gauging station. The Pearson correlation coefficient was used to measure the strength of the statistical relationship between agCDI and detrended crop yield and SRI.

3. Results and Discussions

3.1. Identifying the Historical Drought Years

The grid-based percentage contributions of each of the input variables were computed for each week, and a total of 208 maps were generated. However, the resulting figures for the selected weeks in the Maha and Yala cultivation seasons are shown in Appendix A (Figure A1 and Figure A2). The results illustrate relatively higher contributions of stdNDVI and stdLST, as compared to stdPCP and EDDI. The latter variable contributions are higher in some pocket areas across the study area. In addition, the stdPCP contribution was higher, as compared to other variables in some weeks. Overall, the contribution of EDDI was relatively lower than the other input variables, perhaps due to its coarser spatiotemporal resolution, which impeded the capture of the variability of the evaporative demand across space. The objective-based approach of determining the percentage contribution of the input variables has an advantage over the expert judgment-based approach, which is susceptible to subjectivity and lacks capturing the spatial variability of the weights.
The weekly percentage contribution of each input variable has been utilized to generate the time series of the agricultural drought index (agCDI). Then, agCDI was used to identify the historical drought events in the past 38 years. The number of weeks with severe to extreme drought intensities each year was identified and compared to stipulate the historic drought years. Therefore, 1983, 1984, 1986, 1992, 1994, 2000, 2011, and 2016 were identified as the historical drought years, since a relatively large number of weeks were identified with severe to extreme drought conditions, as compared to the other years (Figure 5). The highest number of weeks of intense drought intensities were observed in 1994, and it covered a vast agricultural area across Sri Lanka. Similarly, a higher number of weeks of severe drought events were observed in 1983 and 1984, mainly in the southwest and central parts, where perennial agriculture (i.e., rubber, coconut, and tea) and other crops were dominantly growing and located relatively in the higher rainfall zones (wet and intermediate zones). Similar spatial patterns of the drought events were observed in 1986 and 1992, as of 1983 and 1984. Conversely, a smaller number of drought events were detected relatively in 2000, 2006, and 2009, except in pocket areas in the central, eastern, and northern parts, where there were more than 20 weeks of severe drought intensities indicated. The prevalence of severe to extreme drought was observed mainly in the dry zone in 2016. Predominantly, rainfed and irrigation-based crop production is taking place in the dry zone, particularly in the northern, north-central, eastern, and south-eastern parts of it. The historic drought years identified in this study are in line with the drought years indicated in previous studies. For example, Abeysingha and Rajapaksha [13] detected some of the drought years, such as 1982–1983, 1986–1987, 2000–2001, and 2016–2017, on their weather stations-based meteorological (SPI) drought assessment in the study period from 1970 to 2017. Ekanayake, and Perera [62] also identified the occurrence of severe drought events in 1986, 1992, 1995, and other years using SPI in the Anuradhapura districts in Sri Lanka. The 2012 and 2016 drought years were spotted in the study by Alahacoon and Ranagalage [62]. Therefore, agCDI has shown its capability to identify the historic drought years in Sri Lanka.
A similar analysis of the frequency of occurrence of drought or the number of weeks under severe and extreme drought conditions was carried out for the Yala and Maha crop growing seasons in Sri Lanka. According to the department of agriculture of Sri Lanka, rice is one of the popular crops that is grown widely, with average cultivation areas of ~748,027 ha in the Maha and 368,906 ha in Yala seasons, with an annual total of 1,116,933 ha, accounting for 34% of the total cultivated area of the island. Out of 25 administrative districts of the country, Ampara, Anuradhapura, Polonnaruwa, Kurunegala, and Hambanthota are the top rice crop growing districts in Sri Lanka. The historical drought years were also identified in both the Yala and Maha seasons (Figure 6). Relatively more severe drought events were captured in both seasons, except 2016, where more weeks under drought conditions were observed in the Maha season.

3.2. Frequency of Occurrence of Drought

Figure 7 illustrates the frequency of the occurrence of mild, moderate, severe, and extreme droughts in Yala and Maha, as well as the annual time scales in the agricultural-dominated areas in Sri Lanka. The occurrence of these drought categories in the study period (1982–2020) was taken into consideration to generate the maps. The frequency of the occurrence of drought in this section refers to the ratio of the number of weeks under each drought category to the total number of weeks of the corresponding seasons within the analysis period. Different legend categories were used for each drought type for better visualization of the spatial distribution of the frequency of the occurrence of drought. The results, in general, illustrate the frequent occurrence of mild droughts, as compared to the other drought categories, in the majority of the agricultural lands. There is a more than 30% chance of the occurrence of mild drought in the annual time scale, while the seasonal analysis demonstrated relatively more prevalent mild drought conditions in the Maha season, as compared to the Yala season. On the other hand, there is somewhat more of a chance of the occurrence of moderate drought (>10%) in the Yala season, compared to the Maha season. The result on an annual time scale further depicts relatively more frequent moderate drought (>8%) in some of the major rice-growing districts, such as Anuradhapura and Hambantota. Moderate drought occurred more frequently next to the mild droughts across majorities of the districts. Severe droughts occurred (2–4%) across the different crop growing districts at an annual time scale, although more than 4% frequent severe droughts have been observed in some pocket areas in central crop growing districts. Relatively more frequent severe droughts have occurred in the Yala season, as compared to the Maha season, in the majority of the districts. Similarly, extreme droughts occurred across a majority of the regions, with a frequency of occurrence of 2%. A higher frequency (>2%) of extreme droughts occurred mainly in the southern region (wet and intermediate zones) and some pocket areas in the dry zone, where agricultural practices are prevalent (eastern part). The frequency of an extreme drought above 3% is often considered significant, since the probability of the occurrence of extreme events is once in a while in a given analysis period [63]. The comparison of the seasonal frequency of the occurrence of severe to extreme drought intensities indicated relatively more frequent severe droughts in the Yala season, as compared to the Maha season, across the agricultural lands. Abeysingha and Rajapaksha [13] indicated more drought events during the Yala season, compared to the Maha season, on their SPI-based spatiotemporal drought assessment across Sri Lanka—the finding somehow supports the results of this study.
Figure 8 presents the frequency of the occurrence of moderate to extreme drought conditions in the dry, intermediate, and wet zones in Sri Lanka. The average areal values of the frequency of occurrences of moderate to extreme drought were generated for each zone. This analysis somehow infers the drought-prone zone in Sri Lanka. The result shows relatively more frequent droughts occur in the dry zone during the Yala season (16%), while a bit less frequent droughts occur in the wet zone. In addition, relatively more severe droughts occur during the Yala seasons in all three zones.

3.3. Spatial Patterns of Drought

The spatial patterns and severity of the selected historical drought events are presented in Figure 9 for the Yala and Maha seasons. The result illustrated prevalent severe drought conditions in most parts during the Yala season, except for the normal conditions in some pocket areas. The vast area had experienced abnormally drier conditions in 1983, 1994, and 2011. The drought in these years lingered in the Maha season, with relatively less severity, except in 2011, which indicated severe drought in November and December. On the other hand, the 2016 drought progressively developed in the Yala season and intensified in the Maha season. Most of the crop-producing districts were struck by moderate to extreme drought conditions in these years. Rajendram [6] reported the occurrence of extreme drought conditions in 1983 and 2016 on the drought analysis, based on SPI, using long-term (1881–2020) observed rainfall data. At the same time, the occurrence of severe drought in the other years (2011–2012 and 1994) was demonstrated by the Disaster Management Center of Sri Lanka [64].
Moreover, Chithranayana and Punyawardena [44] attempted to identify the drought-prone regions using the moisture availability index (MAI) and presented the relatively higher drought vulnerability of districts ascribed in the dry and intermediate zones—which corresponds with the findings of this study. The wet zone is relatively less at risk of drought, although it had experienced moderate to extreme drought conditions in some pocket areas in most of the drought events. Some of the severe drought events were linked to the occurrence of El Niño years. For example, El Niño events in 1982–1983 and 2015–2016 might have triggered the droughts in 1983–1984 and 2016. However, a detailed analysis is required to explore the connection between drought and El Niño occurrences in Sri Lanka. According to Zubair et al. [65], the rise in the Indian Ocean Sea temperature instigated warmer temperatures than usually expected during El Niño event, and this situation exacerbated the drier condition in 2016.

3.4. Evaluation of Drought with Crop Yield and SRI

An attempt has been made to evaluate the agCDI using the district-level crop yield data. The results obtained at selected major crop growing districts have been presented in this section for the Yala and Maha seasons. Figure 10 presents the resulting time series plots in the Yala and Maha seasons. In general, there is a good agreement between the agCDI and detrended crop yield in the five selected districts, with correlation coefficient values ranging between 0.61 to 0.77 and 0.68 to 0.81 in the Yala and Maha seasons, respectively. The agCDI was able to capture the patterns of the crop yield anomalies, except for some discrepancies in some of the drought events. Relatively higher correlation coefficient values were observed in the Ampara (0.77) and Hambantota (0.81) districts in the Yala and Maha seasons, respectively. Conversely, relatively lower correlation values have been noted in Kurunegala (0.61) and Anuradhapura (0.68) in the Yala and Maha seasons. An anticipated significant crop yield reduction was largely observed in the Ampara and Hambantota districts in the Yala season, due to the 2016 drought, while the drought impact intensified largely in these districts during the Maha season, as substantiated by the crop yield reduction. The failure of the rainfall, due to El Niño and the non-formation of usual low-level atmospheric disturbances and/or depressions in the Bay of Bengal, was the principal cause that triggered the 2016 drought. Paddy cultivation was severely affected because of the dry spell in this drought event. According to WFP (2016), 63% of the paddy cultivation was affected in the 2016 Maha season, with a persistently worsening impact on the main agricultural harvest in the next Yala season in 2017 [66]. The result may further improve with the availability of maps that show the specific location of the rice cultivation farmlands. The statistical significance of the correlation coefficient values was tested at a 95% confidence interval. The result indicated there is a positive relationship between the agCDI and detrended crop yield, with alpha values less than 0.05.
The correlation coefficient values between agCDI and paddy crop yield may further improve with the availability of locations and extents of the paddy cultivation land in each district. In the wet zone, perennial agricultural crops, such as tea, coconut, and rubbers, are the predominant agricultural land use, and the drought tolerance of these crops might also be contributed to largely less correlation in the southwestern part, where southwest monsoon rains are highly effective, as per the general rainfall climatology of the island (result not shown). In addition, the paddy cultivation in districts in the wet zone covers relatively less area, as compared to Dry Intermediate zones.
Additional effort has been made to validate the agCDI using another independent dataset, which is streamflow-based SRI. The number of stations used was selected based on the availability of the observed data with reasonable missing values. The hydrological boundary draining to each gauging station was delineated and used to extract the areal average time series values of agCDI for the Yala and Maha seasons. The spider plot in Figure 11 illustrates the correlation coefficient values between agCDI and SRI of each sample station. The result, in general, showed relatively lower correlation coefficient values, as compared to the crop yield, except for the Kude Oya station, where a relatively higher value of correlation coefficient (0.77) was observed. Since the streamflow is a response to the catchment process, there might be some lag time to respond to drought. Future studies that consider the time lag relationship and different aggregate periods of some of the input variables may improve the results in this section.

4. Conclusions

The agricultural drought index (agCDI) developed in this study is generally able to capture and characterize the spatial patterns of the historical drought conditions in Sri Lanka. The study presented the occurrences of moderate to extreme droughts in most agricultural lands across Sri Lanka. In addition, 1983, 1994, 2011, and 2016 were identified as some of the recent historic drought years, with severe drought conditions prevalent in the majority parts of the island. Some of the drought years followed the occurrence of El Niño and the non-formation of usual low-level atmospheric disturbances and depressions in the Bay of Bengal. In addition, the comparison of the seasonal frequency of the occurrence of severe to extreme drought intensities indicated relatively more frequent severe drought in the Yala season, as compared to the Maha season, across the agricultural lands. Therefore, the agCDI can be used to characterize and identify the spatial pattern of the agricultural drought in Sir Lanka.
The evaluation of the seasonal agCDI with the detrended crop yield data illustrated a strong linkage in the major rice growing districts during the Yala and Maha seasons. A correlation coefficient value of greater than 0.65 was observed in the majority of the selected districts, providing fair evidence for the utilization of agCDI in developing agricultural yield prediction models and decision-making processes by local governmental and non-governmental organizations. In addition, agCDI can be used to develop a drought forecasting framework to enhance the preparedness for future droughts and to mitigate its adverse impacts in the agricultural and other important socioeconomic sectors of Sri Lanka.
The potential usefulness of agCDI can be further improved using the additional input variables representing water balance components, such as groundwater, soil moisture, etc., although the input variables considered in this study somehow explain the different water balance components. In addition, comparing the outcome of this study using other data mining approaches, such as random forest, regression tree, etc., may provide better insights towards developing a robust drought monitoring system in future studies.

Author Contributions

The first author Y.B. analyzed the results and prepared the first draft of the manuscript. Y.B. and R.S. framed the research. All the authors including Y.B., R.S., G.J., A.B., A.S., S.A., R.P. and K.D.W.N. contributed to the interpretation of the results and review and provided insightful inputs at several stages of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research is financially supported by the Japan-Bank Program for Mainstreaming DRM in Developing Countries.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data used in this study can be shared without violating any data sharing policy of the corresponding organizations.

Acknowledgments

This paper was made possible with the financial support from the Japan-Bank Program for Mainstreaming DRM in Developing Countries, which is financed by the Government of Japan and receives technical support from the World Bank Tokyo Disaster Risk Management Hub and from the Global Water Security and Sanitation Partnership (GWSP), which is a multi-donor fund within the World Bank that produces cutting-edge research and analytics to create and deliver urgent, practical, and innovative solutions. We also acknowledge the Department of Census and Statistics, Department of Irrigation, Department of Agriculture, and Department of Meteorology of Sri Lanka for providing the hydroclimatic, crop yield, and land use data.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Percentage contribution (weight) of each input variable for the selected weeks in the Maha and Yala seasons.
Figure A1. Maha season.
Figure A1. Maha season.
Water 14 03317 g0a1
Figure A2. Yala season.
Figure A2. Yala season.
Water 14 03317 g0a2

References

  1. Alahacoon, N.; Edirisinghe, M.; Ranagalage, M. Satellite-Based Meteorological and Agricultural Drought Monitoring for Agricultural Sustainability in Sri Lanka. Sustainability 2021, 13, 3427. [Google Scholar] [CrossRef]
  2. Prasanna, R.P.I.R. Economic costs of drought and farmers’ adaptation strategies: Evidence from Sri Lanka. Sri Lanka J. Econ. Res. 2018, 5, 61. [Google Scholar] [CrossRef]
  3. WFP. Sri Lanka-Joint Assessment of Drought Impacts on Food Security and Livelihoods, March 2017. World Food Programme (WFP). 2017. Available online: https://www.wfp.org/publications/Sri_Lanka_Joint_Drought_Assessment (accessed on 5 February 2021).
  4. Manesha, S.; Vimukthini, S.; Premalal, K.H.M.S. Develop Drought Monitoring in Sri Lanka Using Standard Precipitation Index (SPI). Sri Lanka J. Meteorol 2015, 1, 64–71. [Google Scholar]
  5. Nianthi, R. Drought Risk Reduction in the Dry Zone of Sri Lanka. In Droughts in Asian Monsoon Region; Emerald Group Publishing: Bingley, UK, 2011; Volume 8, pp. 97–120. [Google Scholar] [CrossRef]
  6. Rajendram, K. Rainfall variability and drought in the dry and wet zones of Sri Lanka. World Sci. News 2021, 160, 172–189. [Google Scholar]
  7. Alahacoon, N.; Edirisinghe, M. Spatial Variability of Rainfall Trends in Sri Lanka from 1989 to 2019 as an Indication of Climate Change. ISPRS Int. J. Geo-Inf. 2021, 10, 84. [Google Scholar] [CrossRef]
  8. Zubair, L. El Niño-southern oscillation influences on rice production in Sri Lanka. Int. J. Clim. 2002, 22, 249–260. [Google Scholar] [CrossRef] [Green Version]
  9. Rao, C.S.; Gopinath, K.A.; Raju, B.M.K.; Rejani, R.; Venkatesh, G.; Kumari, V.V. Dryland Agriculture in South Asia: Experiences, Challenges and Opportunities. In Innovations in Dryland Agriculture; Springer: Cham, Switzerland, 2016; pp. 345–392. [Google Scholar] [CrossRef]
  10. World Bank. Agriculture, Forestry, and Fishing, Value Added (% of GDP)–Sri Lanka. 2020. Available online: https://data.worldbank.org/indicator/NV.AGR.TOTL.ZS?locations=LK (accessed on 24 February 2021).
  11. Ganapuram, S.; Nagarajan, N.; Balaji, V. Village-level drought vulnerability assessment using geographic information system (GIS). Int. J. Adv. Res. Comput. Sci. Softw. Eng. 2013, 3, 1–10. [Google Scholar]
  12. Ministry of Disaster Management. Sri Lanka National Report on Disaster Risk, Poverty and Human Development Relationship; Ministry of Disaster Management, Vidya Mawatha: Colombo, Sri Lanka, 2009.
  13. Abeysingha, N.S.; Rajapaksha, U.R.L.N. SPI-Based Spatiotemporal Drought over Sri Lanka. Adv. Meteorol. 2020, 2020, 1–10. [Google Scholar] [CrossRef] [Green Version]
  14. Wilhite, D.A.; Glantz, M.H. Understanding: The Drought Phenomenon: The Role of Definitions. Water Int. 1985, 10, 111–120. [Google Scholar] [CrossRef] [Green Version]
  15. Barichivich, J.; Osborn, T.; Harris, I.; van der Schrier, G.; Jones, P. Drought: Monitoring global drought using the self-calibrating Palmer Drought Severity Index. Bull. Am. Meteorol. Soc. 2019, 100, S39–S40. [Google Scholar]
  16. Sanchez, N.; González-Zamora, Á.; Martínez-Fernández, J.; Piles, M.; Pablos, M. Integrated remote sensing approach to global agricultural drought monitoring. Agric. For. Meteorol. 2018, 259, 141–153. [Google Scholar] [CrossRef]
  17. Wilhite, D.A.; Hayes, M.J.; Svoboda, M.D. Drought Monitoring and Assessment: Status and Trends in the United States. In Drought and Drought Mitigation in Europe; Springer: Berlin, Germany, 2000; pp. 149–160. [Google Scholar] [CrossRef]
  18. Ntale, H.K.; Gan, T.Y. Drought indices and their application to East Africa. Int. J. Climatol. A J. R. Meteorol. Soc. 2003, 23, 1335–1357. [Google Scholar] [CrossRef]
  19. Zhao, Q.; Chen, Q.; Jiao, M.; Wu, P.; Gao, X.; Ma, M.; Hong, Y. The Temporal-Spatial Characteristics of Drought in the Loess Plateau Using the Remote-Sensed TRMM Precipitation Data from 1998 to 2014. Remote Sens. 2018, 10, 838. [Google Scholar] [CrossRef] [Green Version]
  20. Rhee, J.; Im, J.; Carbone, G. Monitoring agricultural drought for arid and humid regions using multi-sensor remote sensing data. Remote Sens. Environ. 2010, 114, 2875–2887. [Google Scholar] [CrossRef]
  21. 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] [CrossRef]
  22. Sheffield, J.; Wood, E.F.; Pan, M.; Beck, H.; Coccia, G.; Serrat-Capdevila, A.; Verbist, K. Satellite remote sensing for water resources management: Potential for supporting sustainable development in data-poor regions. Water Resour. Res. 2018, 54, 9724–9758. [Google Scholar] [CrossRef] [Green Version]
  23. Thenkabail, P.S.; Gamage, M.S.D.N. The Use of Remote Sensing Data for Drought Assessment and Monitoring in Southwest Asia; IWMI: Colombo, Sri Lanka, 2004; Volume 85. [Google Scholar]
  24. Lyon, B.; Zubair, L.; Ralapanawe, V.; Yahiya, Z. Finescale Evaluation of Drought in a Tropical Setting: Case Study in Sri Lanka. J. Appl. Meteorol. Clim. 2009, 48, 77–88. [Google Scholar] [CrossRef]
  25. Abeysingha, N.S.; Wickramasuriya, M.G.; Meegastenna, T.J. Assessment of meteorological and hydrological drought; a case study in Kirindi Oya river basin in Sri Lanka. Int. J. Hydrol. Sci. Technol. 2020, 10, 429–447. [Google Scholar] [CrossRef]
  26. Pani, P.; Alahacoon, N.; Amarnath, G.; Bharani, G.; Mondal, S.; Jeganathan, C. Comparison of SPI and IDSI Applicability for Agriculture Drought Monitoring in Sri Lanka. In Proceedings of the 37th Asian Conference on Remote Sensing, Colombo, Sri Lanka, 17–21 October 2016; pp. 17–21. [Google Scholar]
  27. Sumanta, D.A.S.; Choudhury, M.R.; Gandhi, S.; Joshi, V. Application of earth observation data and standardized precipitation index-based approach for meteorological drought monitoring, assessment and prediction over Kutch, Gujarat, India. Int. J. Environ. Geoinform. 2016, 3, 24–37. [Google Scholar]
  28. Herath, H.M.R.C.; Premalal, K.H.M.S.; Kaumudie, A.L.I.; Sanjeewani, D.M.N. Analysis of Standard Precipitation Indices to Identify for Drought Condition in 2015. Sri Lanka J. Meteorol 2015, 1, 20–31. [Google Scholar]
  29. Gunda, T.; Hornberger, G.M.; Gilligan, J.M. Spatiotemporal Patterns of Agricultural Drought in Sri Lanka: 1881–2010. Int. J. Clim. 2015, 36, 563–575. [Google Scholar] [CrossRef]
  30. Sai, M.V.R.S.; Murthy, C.S.; Chandrasekar, K.; Jeyaseelan, A.T.; Diwakar, P.G.; Dadhwal, V.K. Agricultural drought: Assessment & monitoring. MAUSAM 2016, 67, 131–142. [Google Scholar] [CrossRef]
  31. Chathuranga, K.A.M.; Sandamali, K.U.J. Vegetation Condition Index based Agricultural Drought mapping over the past decade of Sri Lanka by utilizing the Satellite Remote Sensing. Built Environ. Spat. Sci. 2020. Available online: http://ir.kdu.ac.lk/handle/345/3251 (accessed on 17 September 2022).
  32. Jayawardhana, W.G.N.N.; Chathurange, V.M.I. Investigate the sensitivity of the satellite-based agricultural drought indices to monitor the drought condition of paddy and introduction to enhanced multi-temporal drought indices. J. Remote Sens. GIS 2020, 9, 272. [Google Scholar]
  33. De Silva, M.; Kawasaki, A. Socioeconomic Vulnerability to Disaster Risk: A Case Study of Flood and Drought Impact in a Rural Sri Lankan Community. Ecol. Econ. 2018, 152, 131–140. [Google Scholar] [CrossRef]
  34. Fragaszy, S.R.; Jedd, T.; Wall, N.; Knutson, C.; Fraj, M.B.; Bergaoui, K.; Svoboda, M.; Hayes, M.; McDonnell, R. Drought Monitoring in the Middle East and North Africa (MENA) Region: Participatory Engagement to Inform Early Warning Systems. Bull. Am. Meteorol. Soc. 2020, 101, E1148–E1173. [Google Scholar] [CrossRef] [Green Version]
  35. Svoboda, M.; LeComte, D.; Hayes, M.; Heim, R.; Gleason, K.; Angel, J.; Rippey, B.; Tinker, R.; Palecki, M.; Stooksbury, D.; et al. The drought monitor. Bull. Am. Meteorol. Soc. 2002, 83, 1181–1190. [Google Scholar] [CrossRef] [Green Version]
  36. Brown, J.F.; Wardlow, B.D.; Tadesse, T.; Hayes, M.J.; Reed, B.C. The Vegetation Drought Response Index (VegDRI): A New Integrated Approach for Monitoring Drought Stress in Vegetation. GIScience Remote Sens. 2008, 45, 16–46. [Google Scholar] [CrossRef]
  37. Tadesse, T.; Champagne, C.; Wardlow, B.D.; Hadwen, T.A.; Brown, J.F.; Demisse, G.B.; Bayissa, Y.A.; Davidson, A.M. Building the vegetation drought response index for Canada (VegDRI-Canada) to monitor agricultural drought: First results. GIScience Remote Sens. 2017, 54, 230–257. [Google Scholar] [CrossRef]
  38. Keyantash, J.A.; Dracup, J.A. An aggregate drought index: Assessing drought severity based on fluctuations in the hydrologic cycle and surface water storage. Water Resour. Res. 2004, 40, W09304. [Google Scholar] [CrossRef] [Green Version]
  39. Bayissa, Y.A.; Tadesse, T.; Svoboda, M.; Wardlow, B.; Poulsen, C.; Swigart, J.; Van Andel, S.J. Developing a satellite-based combined drought indicator to monitor agricultural drought: A case study for Ethiopia. GISci. Remote Sens. 2019, 56, 718–748. [Google Scholar] [CrossRef]
  40. Bayissa, Y.; Tadesse, T.; Demisse, G. Building A High-Resolution Vegetation Outlook Model to Monitor Agricultural Drought for the Upper Blue Nile Basin, Ethiopia. Remote Sens. 2019, 11, 371. [Google Scholar] [CrossRef] [Green Version]
  41. Kulkarni, S.; Wardlow, B.; Bayissa, Y.; Tadesse, T.; Svoboda, M.; Gedam, S. Developing a Remote Sensing-Based Combined Drought Indicator Approach for Agricultural Drought Monitoring over Marathwada, India. Remote Sens. 2020, 12, 2091. [Google Scholar] [CrossRef]
  42. Zubair, L.; Siriwardhana, M.; Chandimala, J.; Yahiya, Z. Predictability of Sri Lankan rainfall based on ENSO. Int. J. Climatol. A J. R. Meteorol. Soc. 2008, 28, 91–101. [Google Scholar] [CrossRef] [Green Version]
  43. Naveendrakumar, G.; Vithanage, M.; Kwon, H.-H.; Iqbal, M.C.M.; Pathmarajah, S.; Obeysekera, J. Five Decadal Trends in Averages and Extremes of Rainfall and Temperature in Sri Lanka. Adv. Meteorol. 2018, 2018, 1–13. [Google Scholar] [CrossRef]
  44. Chithranayana, R.; Punyawardena, B. Identification of drought prone agro-ecological regions in Sri Lanka. J. Natl. Sci. Found. Sri Lanka 2008, 36, 117. [Google Scholar] [CrossRef] [Green Version]
  45. Mapa, R.B.; Kumaragamage, D.; Gunarathne, W.D.L.; Dassanayake, A.R. Land use in Sri Lanka: Past, present and the future. In Proceedings of the 17th World Congress of Social Science (WCSS), Bangkok, Thailand, 14–21 August 2002; pp. 14–21. [Google Scholar]
  46. Punyawardena, B.V.R. Climate of the Dry Zone of Sri Lanka. In Soils of the Dry Zone of Sri Lanka; Special Publication No.7. Soil Science Society of Sri Lanka; Springer: Berlin/Heidelberg, Germany, 2010; pp. 9–26. [Google Scholar]
  47. Funk, C.C.; Peterson, P.J.; Landsfeld, M.F.; Pedreros, D.H.; Verdin, J.P.; Rowland, J.D.; Romero, B.E.; Husak, G.J.; Michaelsen, J.C.; Verdin, A.P. A quasi-global precipitation time series for drought monitoring. US Geol. Surv. Data Ser. 2014, 832, 1–12. [Google Scholar] [CrossRef] [Green Version]
  48. Alahacoon, N.; Edirisinghe, M. A comprehensive assessment of remote sensing and traditional based drought monitoring indices at global and regional scale. Geomat. Nat. Hazards Risk 2022, 13, 762–799. [Google Scholar] [CrossRef]
  49. Bandurathna, L.B.; Wang, L.; Zhou, X.; Cheng, Y.; Chen, L. Intraseasonal oscillation of the southwest monsoon over Sri Lanka and evaluation of its subseasonal forecast skill. Atmos. Ocean. Sci. Lett. 2021, 14, 100062. [Google Scholar] [CrossRef]
  50. Bayissa, Y.; Tadesse, T.; Demisse, G.; Shiferaw, A. Evaluation of Satellite-Based Rainfall Estimates and Application to Monitor Meteorological Drought for the Upper Blue Nile Basin, Ethiopia. Remote Sens. 2017, 9, 669. [Google Scholar] [CrossRef] [Green Version]
  51. Gao, F.; Zhang, Y.; Ren, X.; Yao, Y.; Hao, Z.; Cai, W. Evaluation of CHIRPS and its application for drought monitoring over the Haihe River Basin, China. Nat. Hazards 2018, 92, 155–172. [Google Scholar] [CrossRef]
  52. Saeidizand, R.; Sabetghadam, S.; Tarnavsky, E.; Pierleoni, A. Evaluation of CHIRPS rainfall estimates over Iran. Q. J. R. Meteorol. Soc. 2018, 144, 282–291. [Google Scholar] [CrossRef] [Green Version]
  53. Liu, J.; Shangguan, D.; Liu, S.; Ding, Y.; Wang, S.; Wang, X. Evaluation and comparison of CHIRPS and MSWEP daily-precipitation products in the Qinghai-Tibet Plateau during the period of 1981–2015. Atmos. Res. 2019, 230, 104634. [Google Scholar] [CrossRef]
  54. Abdelmoneim, H.; Soliman, M.R.; Moghazy, H.M. Evaluation of TRMM 3B42V7 and CHIRPS Satellite Precipitation Products as an Input for Hydrological Model over Eastern Nile Basin. Earth Syst. Environ. 2020, 4, 685–698. [Google Scholar] [CrossRef]
  55. Kogan, F. Remote Sensing for Food Security. In Remote Sensing for Food Security; Springer: Cham, Switzerland, 2019; pp. 51–73. [Google Scholar] [CrossRef]
  56. Hobbins, M.T.; Wood, A.; McEvoy, D.J.; Huntington, J.L.; Morton, C.; Anderson, M.; Hain, C. The Evaporative Demand Drought Index. Part I: Linking Drought Evolution to Variations in Evaporative Demand. J. Hydrometeorol. 2016, 17, 1745–1761. [Google Scholar] [CrossRef]
  57. Peters, A.J.; Walter-Shea, E.A.; Ji, L.; Vina, A.; Hayes, M.; Svoboda, M.D. Drought monitoring with NDVI-based standardized vegetation index. Photogramm. Eng. Remote Sens. 2002, 68, 71–75. [Google Scholar]
  58. Lu, J.; Carbone, G.J.; Gao, P. Mapping the agricultural drought based on the long-term AVHRR NDVI and North American Regional Reanalysis (NARR) in the United States, 1981–2013. Appl. Geogr. 2019, 104, 10–20. [Google Scholar] [CrossRef]
  59. Nanzad, L.; Zhang, J.; Tuvdendorj, B.; Nabil, M.; Zhang, S.; Bai, Y. NDVI anomaly for drought monitoring and its correlation with climate factors over Mongolia from 2000 to 2016. J. Arid Environ. 2019, 164, 69–77. [Google Scholar] [CrossRef]
  60. McKee, T.B.; Doesken, N.J.; Kleist, J. The Relationship of Drought Frequency and Duration to Time Scales. In Proceedings of the 8th Conference on Applied Climatology, Anaheim, CA, USA, 17–22 January 1993; No. 22. pp. 179–183. [Google Scholar]
  61. FAO. Production Variability and Losses. In Special: Agroclimatic Concepts. Sustainable Development Department (SD), R. Gommes, Ed.; Food Agriculture Organization of the United Nations (FAO): Rome, Italy, 1999; Available online: http://www.fao.org/nr/climpag/agroclim/losses_en.asp (accessed on 6 March 2022).
  62. Ekanayake, E.M.R.S.B.; Perera, K. Analysis of Drought Severity and Duration Using Copulas in Anuradhapura, Sri Lanka. Br. J. Environ. Clim. Chang. 2014, 4, 312–327. [Google Scholar] [CrossRef]
  63. Modarres, R.; Sarhadi, A.; Burn, D.H. Changes of extreme drought and flood events in Iran. Glob. Planet. Chang. 2016, 144, 67–81. [Google Scholar] [CrossRef]
  64. Chandrasekara, S.S.; Kwon, H.-H.; Vithanage, M.; Obeysekera, J.; Kim, T.-W. Drought in South Asia: A Review of Drought Assessment and Prediction in South Asian Countries. Atmosphere 2021, 12, 369. [Google Scholar] [CrossRef]
  65. Zubair, L.; Yahiya, Z.; Agalawatte, P.; Lokuhetti, R. The El Niño event of 2015/16 in Sri Lanka predictions, preparedness, communication, and impacts. Neela Har. Clim. Chang. Mag. Sri Lanka 2016, II, 40–46. [Google Scholar]
  66. WFP. Sri Lanka-Initial Rapid Assessment on Drought 2016. World Food Programme (WFP). 2016. Available online: https://www.wfp.org/publications/Sri_Lanka_Drought_Assessment (accessed on 10 April 2022).
Figure 1. Topography and climatic zones of Sri Lanka.
Figure 1. Topography and climatic zones of Sri Lanka.
Water 14 03317 g001
Figure 2. The spatial distribution of the agricultural land use in Sri Lanka.
Figure 2. The spatial distribution of the agricultural land use in Sri Lanka.
Water 14 03317 g002
Figure 3. Summary of the method used in this study.
Figure 3. Summary of the method used in this study.
Water 14 03317 g003
Figure 4. Location of the selected five districts (red boundaries) used for crop yield-based evaluation of the agCDI. In addition, the locations and extents of the hydrological boundaries (yellow background) are shown and labeled using the names of the river gauging stations.
Figure 4. Location of the selected five districts (red boundaries) used for crop yield-based evaluation of the agCDI. In addition, the locations and extents of the hydrological boundaries (yellow background) are shown and labeled using the names of the river gauging stations.
Water 14 03317 g004
Figure 5. Historical drought years identified based on the number of severe to extreme drought weeks in each year.
Figure 5. Historical drought years identified based on the number of severe to extreme drought weeks in each year.
Water 14 03317 g005
Figure 6. Historical drought years identified based on the number of weeks with severe to extreme drought conditions in Yala and Maha crop growing seasons.
Figure 6. Historical drought years identified based on the number of weeks with severe to extreme drought conditions in Yala and Maha crop growing seasons.
Water 14 03317 g006
Figure 7. Frequency of occurrence of mild, moderate, severe, and extreme droughts in Yala, Maha and annual.
Figure 7. Frequency of occurrence of mild, moderate, severe, and extreme droughts in Yala, Maha and annual.
Water 14 03317 g007
Figure 8. Frequency of occurrence of moderate to extreme drought conditions in dry, intermediate, and wet zones Yala and Maha seasons and annual time scale.
Figure 8. Frequency of occurrence of moderate to extreme drought conditions in dry, intermediate, and wet zones Yala and Maha seasons and annual time scale.
Water 14 03317 g008
Figure 9. Spatial patterns and severity of selected historic drought years (i.e., 1983, 1994, 2011, and 2016) in Yala and Maha seasons.
Figure 9. Spatial patterns and severity of selected historic drought years (i.e., 1983, 1994, 2011, and 2016) in Yala and Maha seasons.
Water 14 03317 g009
Figure 10. Evaluation of agCDI with detrended crop yield for the selected crop growing districts for the Yala and Maha seasons.
Figure 10. Evaluation of agCDI with detrended crop yield for the selected crop growing districts for the Yala and Maha seasons.
Water 14 03317 g010
Figure 11. Evaluation of agCDI with SRI largely in the southwest parts of the districts for Yala and Maha seasons.
Figure 11. Evaluation of agCDI with SRI largely in the southwest parts of the districts for Yala and Maha seasons.
Water 14 03317 g011
Table 1. Input variables used to monitor agricultural and water supply droughts and their spatial and temporal resolutions.
Table 1. Input variables used to monitor agricultural and water supply droughts and their spatial and temporal resolutions.
SnoInput VariablesResolutionData Length Source
Spatial Temporal
1CHIRPS or other satellite or observed rainfall (SPI- at various time scales) 5 kmdaily1982-presentUniversity of California, Santa Barbara
2Normalized difference vegetation index (NDVI)4 kmweekly1982-presentNOAA STAR
3Land surface temperature (LST)4 kmweekly1982-presentNOAA STAR
4Evaporative demand drought Index (EDDI)12.5 kmMonthly1982-presentNOAA
Table 2. Area coverage of the land use classes shown in Figure 2.
Table 2. Area coverage of the land use classes shown in Figure 2.
Land Use TypeArea Coverage (%)
Paddy13
Crop lands5
Vegetated area16
Perennial agriculture—rubber4
Perennial agriculture—tea3
Perennial agriculture—coconut6
Woody perennial crops17
Other land uses 28
Water bodies7
Table 3. Drought category based on McKee et al. [60].
Table 3. Drought category based on McKee et al. [60].
SPI ValuesDrought Category
−2.00 and lessExtreme drought
−1.50 to −1.99Severe drought
−1.00 to −1.49Moderate drought
0 to −0.99Near normal or mild drought
Above 0No drought
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bayissa, Y.; Srinivasan, R.; Joseph, G.; Bahuguna, A.; Shrestha, A.; Ayling, S.; Punyawardena, R.; Nandalal, K.D.W. Developing a Combined Drought Index to Monitor Agricultural Drought in Sri Lanka. Water 2022, 14, 3317. https://doi.org/10.3390/w14203317

AMA Style

Bayissa Y, Srinivasan R, Joseph G, Bahuguna A, Shrestha A, Ayling S, Punyawardena R, Nandalal KDW. Developing a Combined Drought Index to Monitor Agricultural Drought in Sri Lanka. Water. 2022; 14(20):3317. https://doi.org/10.3390/w14203317

Chicago/Turabian Style

Bayissa, Yared, Raghavan Srinivasan, George Joseph, Aroha Bahuguna, Anne Shrestha, Sophie Ayling, Ranjith Punyawardena, and K. D. W. Nandalal. 2022. "Developing a Combined Drought Index to Monitor Agricultural Drought in Sri Lanka" Water 14, no. 20: 3317. https://doi.org/10.3390/w14203317

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