Characterizing Meteorological Droughts in Nepal: A Comparative Analysis of Standardized Precipitation Index and Rainfall Anomaly Index

: Drought is an environmental disaster related to the extremes (on a drier side) in hydromete-orology. The precipitation amount modulates drought in Nepalese river basins. It is vital for efﬁcient water resources management to quantify and understand drought. This paper aims to characterize the droughts in Nepal based on standard precipitation index (SPI) and rainfall anomaly index (RAI) using daily precipitation data and assess their impacts on annual crop yields. We used 41 years (1975–2015) of daily precipitation data to compute monthly means and then the drought indices, namely, SPI and RAI, at 123 stations across Nepal. Results showed that the northwest and eastern regions experienced drought compared to other regions, although the severity and duration were shorter. For stations 101 and 308, we found extreme drought events after 2005 for SPI-1, SPI-3, and SPI-6. However, for SPI-6, extreme drought was also observed in 1989 and 1994 at both stations. The year 1992 was one of the severest drought years for the western and northwest regions, where the severity crossed more than − 2.0 for all SPI months. Similar to SPI, RAI also revealed a similar degree of drought in the country. RAI showed that the eastern region depicted a higher degree of severity of drought compared to other areas beyond 2004. The lesser severity is also seen in the far west part beyond 2005. The results showed that SPI and RAI could equally be used to analyze drought severity. More frequent drought incidents have been observed after 2010 at all the considered precipitation stations. With the increase in the drought severity, the crop yield (such as paddy, maize, barley, millet, and wheat) has been directly impacted. These results might be signiﬁcant for planning water resource and irrigation water management systems. extreme drought events shows the rising tendency of such events in the last few years. During the last ﬁve years (2011–2015), the total number of events has increased from 243 in 2006–2010 to 441. The total number of drought events are 183 (1985–1990), 231 (1991–1995), 175 (1996–2000), and 212 (2001–2005), which implies an increasing nature of the number of extreme events. The monthly and seasonal variations for different decades—1985–1995, 1996–2005, 2006–2015 are also analyzed. The analysis shows that the number of events is higher than in earlier decades. Many such events are observed from May to September, implying the shifting of rainy days exacerbated by the changing climate at regional and local levels.


Introduction
Drought is a recurrent natural hazard caused by precipitation deficiency resulting in damage to crops, water resources, economies, and human lives [1] over an extended space-time domain. Drought has a significant impact on various aspects of the global ecosystem [2][3][4][5]. The drought phenomena can be meteorological, hydrological, agricultural, or/and socio-economic [6][7][8]. Meteorological drought is defined as the lack of precipitation over a region for a more extended period, and hydrological drought is related to the low water supply in streams, reservoirs, and groundwater levels [9]. Agricultural drought usually refers to declining soil moisture and consequent crop failure without reference to surface water resources. In addition, socioeconomic drought is associated with the inability of water resources systems to meet water demands [10,11]. A significant relationship Earth 2022, 3 411 the current study also tests the suitability of RAI and SPI to characterize the drought events in the case of Nepal; it is the first of this kind of study in the country.
In this context, spatio-temporal characterization of drought using different techniques is essential for crop planning, irrigation scheduling, irrigation infrastructure development, and planning of irrigation facilities. That information helps the policymakers, farmers, and researchers for early planning of the crop calendar and associated irrigation scheduling in that area. Furthermore, proper planning could increase crop production, retain the rural area population, and motivate them to engage in agriculture. Hence, this study uses long-term rainfall data to characterize the drought using SPI and RAI techniques across Nepal.

Study Area
The study area lies between 26 • 12' N and 30 • 27' N latitudes and 80 • 04' E and 88 • 12' E longitudes. The country covers 147,516 km 2 , ranging from 60 m in the south to 8848 m above mean sea level in the north, thereby demonstrating high climate variability. The country is divided into seven provinces: (1) Province 1, (2) Madhesh, (3) Bagmati, (4) Gandaki, (5) Lumbini, (6) Karnali, and (7) Sudurpaschim, extending from east to west ( Figure 1). The country is rich in water resources reflecting the importance of water-intensive services such as agriculture production, inland navigation, and hydropower production. Almost two-thirds of the population is engaged in agriculture, although it is in a low development stage with low competitiveness and low productivity [47]. One of the significant enabling environments for enhancing agricultural productivity is irrigation facilities like improved seeds, fertilizers, innovative tools, and farm mechanization techniques.
In Nepal, approximately 40% of the total agricultural areas are irrigated, and the rest of the land's agricultural activities solely depends on the rainfall. Karnali, Narayani, and Koshi are the major rivers that flow within the country with high river flow and potential for navigation. The country also poses medium-sized rivers such as Bagmati, Babai, Rapti, Kamala, Kankai, etc., which act as a source of surface water for irrigation and hydropower generation. The country has six seasons: spring (mid-March-mid-May, Basanta), summer (mid-May-mid-July, Grishma), rainy (mid-July-mid-September, Barsha), autumn (mid-September-mid-November, Sharad), pre-winter (mid-November-mid-January, Hemanta), and winter (mid-January-mid-March, Shishir), with spatial and temporal variability on the amount of the rain received. The rainy (Barsha) season receives almost 80% of the rainfall, leaving the remaining seasons as moderate and lesser rainfall seasons. In the Tarai (southern Nepal), summer temperatures exceed 37 • C and are higher in some areas, and winter temperatures range from 7 • C to 23 • C in Tarai. In mountainous regions, hills, and valleys, summers are temperate, whereas winter temperatures can plummet below zero.

Data
We collected the required rainfall data from the Department of Hydrology and Meteorology (DHM), Kathmandu, Nepal, at a daily time scale for 41 years (1975-2015). We then filled in the missing daily rainfall data using an average station approach, one of the popular methods used in the Himalayan region. The daily data are then converted to a monthly time scale for calculating the SPI values in R-programming language.
Similarly, we obtained the annual crop data from the Ministry of Agricultural Development for 1980-2011. In Nepal, cereal crops such as maize, paddy, millet, barley, and others are prepotent. Among these, paddy is water intensive and depends on rainfall and irrigation infrastructures. Paddy is grown all across the country but is more produced in the Tarai (southern belt) and hilly regions (northern belt). The Tarai region is also considered the food basket of the country. The temporal change in the rainfall pattern alters rice production in Tarai. Millet and wheat require a temperate and tropical climate for cultivation. The crops like maize and millet are grown in Tarai and Hilly regions during dry seasons. We chose the crops of paddy, maize, wheat, and barley for the analysis as they

Data
We collected the required rainfall data from the Department of Hydrology and Meteorology (DHM), Kathmandu, Nepal, at a daily time scale for 41 years (1975-2015). We then filled in the missing daily rainfall data using an average station approach, one of the popular methods used in the Himalayan region. The daily data are then converted to a monthly time scale for calculating the SPI values in R-programming language.
Similarly, we obtained the annual crop data from the Ministry of Agricultural Development for 1980-2011. In Nepal, cereal crops such as maize, paddy, millet, barley, and others are prepotent. Among these, paddy is water intensive and depends on rainfall and irrigation infrastructures. Paddy is grown all across the country but is more produced in the Tarai (southern belt) and hilly regions (northern belt). The Tarai region is also considered the food basket of the country. The temporal change in the rainfall pattern alters rice production in Tarai. Millet and wheat require a temperate and tropical climate for cultivation. The crops like maize and millet are grown in Tarai and Hilly regions during dry seasons. We chose the crops of paddy, maize, wheat, and barley for the analysis as they are grown widely across the country throughout the year. The crop yield of any region is more dependent on the amount of rainfall, land use and management practices, and irrigation systems.

Methods
We show the overall methodological framework that characterizes the meteorological drought of Nepal in Figure 2. We analyzed the impact of meteorological droughts on different major crops such as paddy, maize, barley, and millet at various districts in Nepal. The detailed description and calculation method of the approaches used in the current study are described in subsequent sections.

Standardized Precipitation Index (SPI)
Standardized precipitation index (SPI) is a broadly used drought index to characterize meteorological drought and was developed by Mckee et al. [24]. The SPI application is not limited to drought only; itis also used for frequency analysis and climate impact studies. As per the World Meteorological Organization (WMO) recommendation, SPI is a meteorological drought index in which precipitation is the primary influencing climatic parameter [48]. The SPI quantifies observed precipitation as a standardized departure from a selected probability distribution function that models the raw precipitation data. The raw precipitation data are typically fitted to a gamma or a Pearson type III distribution and then Earth 2022, 3 413 transformed to a normal distribution. Therefore, the SPI values can be interpreted as the number of standard deviations the observed anomaly deviates from the long-term mean. It uses monthly precipitation aggregates at various time scales (1, 3, 6, 12, 18, and 24 months, etc.). The SPI value ranges from −2.00 (dry) to +2.00 (wet). As a result, the index shows a different degree of severity, as represented in Table 1. The calculation is based on the gamma distribution. We used the SPEI package in R-programming software to calculate SPI at various periods [49].

Methods
We show the overall methodological framework that characterizes the meteorological drought of Nepal in Figure 2. We analyzed the impact of meteorological droughts on different major crops such as paddy, maize, barley, and millet at various districts in Nepal. The detailed description and calculation method of the approaches used in the current study are described in subsequent sections.

Figure 2.
An overall methodological framework to access the meteorological drought characteristics in Nepal using two drought indices.

Standardized Precipitation Index (SPI)
Standardized precipitation index (SPI) is a broadly used drought index to characterize meteorological drought and was developed by Mckee et al. [24]. The SPI application is not limited to drought only; itis also used for frequency analysis and climate impact studies. As per the World Meteorological Organization (WMO) recommendation, SPI is a meteorological drought index in which precipitation is the primary influencing climatic parameter [48]. The SPI quantifies observed precipitation as a standardized departure from a selected probability distribution function that models the raw precipitation data. The raw precipitation data are typically fitted to a gamma or a Pearson type III distribution and then transformed to a normal distribution. Therefore, the SPI values can be interpreted as the number of standard deviations the observed anomaly deviates from the long-term mean. It uses monthly precipitation aggregates at various time scales (1, 3, 6, 12, 18, and 24 months, etc.). The SPI value ranges from −2.00 (dry) to +2.00 (wet). As a result, the index shows a different degree of severity, as represented in Table 1. The calculation is based on the gamma distribution. We used the SPEI package in R-programming software to calculate SPI at various periods [49].  Figure 2. An overall methodological framework to access the meteorological drought characteristics in Nepal using two drought indices.

SPI Calculation
We computed the SPIs for both space and time using the method developed by Mckee et al. [24]. The method uses monthly precipitation aggregates at various time scales (1, 3, 6, 12, 18, and 24 months, etc.). For an instance, to calculate SPI for a 3-month time scale, the precipitation accumulated from month j − 2 to month j is summed and attributed to month j. At this time scale, the first two months of the data time series are missing. Then, the normalization procedure is conducted, in which an appropriate probability density function is first fitted to the long-term time series of aggregated precipitation. Finally, the fitted function is used to calculate the cumulative distribution of the data points, which are finally transformed into standardized normal variates. This procedure is repeated for all needed time scales. The maximum-likelihood (ML) estimation method [50] was used to fit gamma probability distribution functions to each time series. The detail of the computation is given by Equations (1)-(5) as follows: The probability density function of the gamma distribution is defined as where x > 0 is the amount of precipitation, α > 0 and β > 0 are the shape and scale parameters, respectively, and Γ(α) is the gamma function. Detailed explanations on the gamma distribution can be found in Lloyd-Hughes and Saunders [51] and Guttman [26].
To fit the distribution parameters, α and β are estimated from the sample data using the approximation for ML defined by Thom [52] which are given by: where x is mean precipitation and A is given by A = ln(x) − n −1 ∑ ln(x). For a given month and time scale, the cumulative probability G(x) of an observed amount of precipitation is given by Let us consider t = x/β, the expression is reduced to the incomplete gamma function as follows: The gamma distribution is not defined for x = 0, and, the probability of zero precipitation q = P(x = 0) being positive, the actual probability of non-exceedance H(x) is calculated as follows: Finally, the cumulative probability distribution H(x) is transformed into the standard normal distribution to yield the SPI.

Rainfall Anomaly Index (RAI)
The RAI, developed by Rooy [25], is used to classify the positive and negative severities in rainfall anomalies. It is also considered an index of remarkable procedural simplicity (like SPI) because it requires only precipitation data [38,53]. Araújo et al. [54] employed the RAI concept to analyze the climatology at the Paraiba river basin. Similarly, Sanches et al. [55] also applied the RAI concept to analyze the precipitation variability for more than seven decades  in Alegrete. These studies concluded that RAI could be a vital tool in analyzing the precipitation data over the study area from the research results. Similarly, Santos et al. [56] utilized RAI to analyze the Mamanguape River basin's climate. They found that the river basin can be divided into three regions in terms of precipitation patterns. The research results further concluded that RAI is an alternative tool in observing the area's precipitation pattern. The degree of severity for the RAI index is shown in Table 2. The degree of severity varies from −4 (extremely dry) to +4 (extremely humid/wet). We used Microsoft Excel to compute the RAI values using Equations (6) and (7). We computed RAI for both annual and monthly periods using Equations (6) and (7): for the negative anomaly (7) where N is the current monthly/seasonal/annual precipitation; N is the average monthly/seasonal/annual precipitation of historical series; M is the average of ten highest monthly/seasonal/annual precipitations; X is the average of ten lowest monthly/seasonal/annual precipitations; N − N represents the positive anomaly and negative anomaly based on positive or negative values.
RAI is usually used to compare the deviation in precipitation in different regions.

Identification of Contributing Determinants for Agricultural Production
Irrigation facilities and soil moisture conditions are the critical determinants for agricultural production. Suppose there are no irrigation facilities at the field level; in that case, crop production mainly depends upon the soil moisture condition and available rainfall on that particular crop period, excluding other enabling environments like an improved seed, fertilizer, tools, and techniques. If the agricultural area is rainfed, the meteorological drought index could relate to agricultural production. Figure 1 shows the average annual precipitation across different provinces, and Figure 3 shows the temporal distribution of average monthly precipitation for 41 years from 1975 to 2015. Gandaki province's central region receives the maximum amount (>3500 mm) of the annual rainfall ( Figure 1). In contrast, the northwest of the same province receives the minimum (<200 mm) amount of yearly rainfall. On average, Karnali province receives rainfall in the range of 200-1600 mm, whereas Bagmati province and Province 1 receive the annual rainfall in a range of 600-3500 mm. The monthly rainfall is more than 200 mm from June to September, also a rainy/monsoon season ( Figure 3). Approximately 80% of the total rainfall is received in these four months only. The remaining eight months (October to May) get monthly rainfall of less than 200 mm. It clearly shows July experiences the maximum rainfall (602 mm) followed by August (585 mm), June (452 mm), and September (394 mm). High temporal heterogeneous distribution of the rainfall reflects the possibility of extreme events such as drought and flood, which need to be addressed for proper watershed management.

Temporal Analysis
On the one-month time scale, SPI-1 is plotted for 41 years (1975-2015) for two representative meteorological stations, station 101 and station 308, to represent the drought state at the point level. These two stations were chosen as representative because they Earth 2022, 3 416 receive relatively less annual and monsoon rainfall than other stations. Figure 4 shows SPI-1 values, temporal distribution, two months of January and August, and cumulative distribution function (CDF) at Station 101. The station is located in the far western region in Sudurpaschim province. Similarly, Figure 5 shows the same for station 308, located in the mid-western region at Karnali province. The figures show that 1984, 1997, and 2012 experienced severe drought at station 101, whereas station 308 experienced severe drought in 2008 and 2014. On average, both the stations experienced near-normal drought (−0.99 > SPI < 0.99) throughout the study period, however, both the intensity and frequency were intense in the later years.
Earth 2022, 3, rainfall in the range of 200-1600 mm, whereas Bagmati province and Province 1 r the annual rainfall in a range of 600-3500 mm. The monthly rainfall is more than 20 from June to September, also a rainy/monsoon season ( Figure 3). Approximately 8 the total rainfall is received in these four months only. The remaining eight month tober to May) get monthly rainfall of less than 200 mm. It clearly shows July exper the maximum rainfall (602 mm) followed by August (585 mm), June (452 mm), an tember (394 mm). High temporal heterogeneous distribution of the rainfall reflec possibility of extreme events such as drought and flood, which need to be address proper watershed management.

Temporal Analysis
On the one-month time scale, SPI-1 is plotted for 41 years (1975-2015) for two sentative meteorological stations, station 101 and station 308, to represent the dr state at the point level. These two stations were chosen as representative becaus receive relatively less annual and monsoon rainfall than other stations. Figure 4 SPI-1 values, temporal distribution, two months of January and August, and cumu distribution function (CDF) at Station 101. The station is located in the far western r in Sudurpaschim province. Similarly, Figure 5 shows the same for station 308, loca the mid-western region at Karnali province. The figures show that 1984, 1997, and experienced severe drought at station 101, whereas station 308 experienced drought in 2008 and 2014. On average, both the stations experienced near-normal dr (−0.99 > SPI < 0.99) throughout the study period, however, both the intensity an quency were intense in the later years. The monthly distribution of the one-month drought shows the lesser intensity of SPI wet values and higher intensity of SPI drought values in January. The number of wet days is higher during August. The high number of wet days in August might be because of the monsoon season (June-September). Monthly distribution of the drought events at station 101 shows high frequency in January. The station experienced continuous drought from 2006 to 2011 in January. The same station experienced drought from 2001 to 2006, even in August. Some years (1979,1997,2001,2009, and 2011), even in August, experienced moderate to severe drought (+1 < SPI < +2) at station 308, with the number of drought events getting more frequent beyond 2000. Simultaneously, we found that the number and intensity of moderate to severe droughts are high during January. As a result, an increase in the number and intensity of droughts directly impact crop yield and other hazards, such as increased heat strokes.
The CDF for stations 101 and 308 is plotted for the wet (August) and dry (January) months. The probability value ranged from 0.07 to 0.94 against the variation of SPI values (−1.45 to +1.54) in January at station 101. Similarly, the probabilistic variation of CDF varied from 0.03 to 0.99 with the variation of SPI values of −1.92 to +2.32 in August at the same station. In January, the probability ranged from 0.12 to 0.92 against the SPI −1.17 to +1.4. Similarly, the probability ranged from 0.02 to +1.00 for the SPI −1.97 +2.69 in August. The other stations also show a similar nature in probability distribution compared to the The monthly distribution of the one-month drought shows the lesser intensity of SPI wet values and higher intensity of SPI drought values in January. The number of wet days is higher during August. The high number of wet days in August might be because of the monsoon season (June-September). Monthly distribution of the drought events at station 101 shows high frequency in January. The station experienced continuous drought from 2006 to 2011 in January. The same station experienced drought from 2001 to 2006, even in August. Some years (1979,1997,2001,2009, and 2011), even in August, experienced moderate to severe drought (+1 < SPI < +2) at station 308, with the number of drought events getting more frequent beyond 2000. Simultaneously, we found that the number and intensity of moderate to severe droughts are high during January. As a result, an increase in the number and intensity of droughts directly impact crop yield and other hazards, such as increased heat strokes.
The CDF for stations 101 and 308 is plotted for the wet (August) and dry (January) months. The probability value ranged from 0.07 to 0.94 against the variation of SPI values (−1.45 to +1.54) in January at station 101. Similarly, the probabilistic variation of CDF varied from 0.03 to 0.99 with the variation of SPI values of −1.92 to +2.32 in August at the same station. In January, the probability ranged from 0.12 to 0.92 against the SPI −1.17 to +1.4. Similarly, the probability ranged from 0.02 to +1.00 for the SPI −1.97 +2.69 in August. The northwest region has low precipitation compared to the mid-western and eastern regions, which have moderate and high precipitation. The temporal distribution of various precipitation stations revealed the SPI values in the northwest region is in drought between 2000 and 2010, with increased severity of drought in 2002. Similarly, there has been no severe drought in the eastern region for the past 41 years, but the onset and termination of the drought are frequent, i.e., in cyclic order. The mid-western region has less drought time compared to the other areas. We found that the annual average precipitation in the western region is low, and it is moderately high in the eastern region. For most of the monsoon period, 1975-1985, the recorded precipitation is unavailable in the eastern region, which might be the reason for the drought for the first ten years. For other years, severe drought is not observed in the eastern and western regions. This fact reveals that these areas have the possibility of rain-fed activities such as agriculture and hydropower production (where the head is available). The northwest is under drought but not severe, which might be the consequence of arctic oscillation and persistent warming of the Indian Ocean [57]. The drought severity is high in the country after 2000, which is in agreement with the result Earth 2022, 3 418 obtained by Hamal et al. [58]. In a nutshell, drought in the country will significantly affect the environment, economy, and social livelihoods [59].
Earth 2022, 3,10 The other stations also show a similar nature in probability distribution compared to the SPI values. The CDF distribution of both dry and wet months shows that the SPI events were fitted entirely with the normal distribution. The northwest region has low precipitation compared to the mid-western and eastern regions, which have moderate and high precipitation. The temporal distribution of various precipitation stations revealed the SPI values in the northwest region is in drought between 2000 and 2010, with increased severity of drought in 2002. Similarly, there has been no severe drought in the eastern region for the past 41 years, but the onset and termination of the drought are frequent, i.e., in cyclic order. The mid-western region has less drought time compared to the other areas. We found that the annual average precipitation in the western region is low, and it is moderately high in the eastern region. For most of the monsoon period, 1975-1985, the recorded precipitation is unavailable in the eastern region, which might be the reason for the drought for the first ten years. For other years, severe drought is not observed in the eastern and western regions. This fact reveals that these areas have the possibility of rain-fed activities such as agriculture and hydropower production (where the head is available). The northwest is under drought but not severe, which might be the consequence of arctic oscillation and persistent warming of the Indian

Temporal Fluctuation of SPI Events
The temporal fluctuation of SPI values for stations 101 and 308 is shown in Figure 6. For most of the years, there is severe drought (i.e., SPI > −2) in all the SPI evolutions in the considered meteorological stations. However, both the stations experienced extreme drought (SPI < −2) after 2005 in SPI-1, SPI-3, and SPI-6 time evolutions. For SPI-6, extreme drought was observed in 1989 and 1994 at both stations. In contrast, results of the SPI-12 show no extreme drought at either station but a wet year after 2005. Figure 7 shows the spatial distribution of the one-month SPI in half-decadal time (5-year) scales. We observed from normal drought to moderate wetness in all periods. In particular, normal droughts were detected in more expansive areas in 2005, 2010, and 2015. We found that the eastern region (Province 1) and northwest region (Karnali and Sudurpaschim province) have normal drought, whereas the remaining area experienced normal wet in 2005. The spatial distribution of the normal drought was extended to Madhesh province and Bagmati province in the years 2010 and 2015. The temporal fluctuation of SPI values for stations 101 and 308 is shown in Figure 6. For most of the years, there is severe drought (i.e., SPI > −2) in all the SPI evolutions in the considered meteorological stations. However, both the stations experienced extreme drought (SPI < −2) after 2005 in SPI-1, SPI-3, and SPI-6 time evolutions. For SPI-6, extreme drought was observed in 1989 and 1994 at both stations. In contrast, results of the SPI-12 show no extreme drought at either station but a wet year after 2005.  year) scales. We observed from normal drought to moderate wetness in all periods. In particular, normal droughts were detected in more expansive areas in 2005, 2010, and 2015. We found that the eastern region (Province 1) and northwest region (Karnali and Sudurpaschim province) have normal drought, whereas the remaining area experienced normal wet in 2005. The spatial distribution of the normal drought was extended to Madhesh province and Bagmati province in the years 2010 and 2015.  Figures S1-S4). The figures depict that the severity of drought increases as the SPI time scale increases. The central (Bagmati and Gandaki provinces) and eastern (Province 1) regions were found to have a higher degree of drought than other regions in all the SPI time scales. However, the average monthly distribution of SPI values ranged from +0.52 (July) to +0.02 (January). SPI spatial distribution showed that the northern regions were less vulnerable to drought compared to the central and eastern regions, which is equally represented by Cavus and Aksoy [60]. The result further stated that the degree of drought severity reduced with elevation. The monsoon season's unpredictable change might have resulted in the drought vulnerability in some monsoon years [61]. The spatial distribution of SPI drought showed the increasing drought trend in  Figures S1-S4). The figures depict that the severity of drought increases as the SPI time scale increases. The central (Bagmati and Gandaki provinces) and eastern (Province 1) regions were found to have a higher degree of drought than other regions in all the SPI time scales. However, the average monthly distribution of SPI values ranged from +0.52 (July) to +0.02 (January). SPI spatial distribution showed that the northern regions were less vulnerable to drought compared to the central and eastern regions, which is equally represented by Cavus and Aksoy [60]. The result further stated that the degree of drought severity reduced with elevation. The monsoon season's unpredictable change might have resulted in the drought vulnerability in some monsoon years [61]. The spatial distribution of SPI drought showed the increasing drought trend in the country's far western region beyond 2005. A similar trend was observed in the research conducted by Khatiwada and Pandey [59] in the Karnali River Basin (Karnali Province) in the far west. The country's central region also received normal drought beyond 2005, as observed by Dahal et al. [46]. Figure 8 shows the monthly SPI-1 events for 1975-2015 on a decadal time scale. The figure reveals the variation in SPI-1 values in the different months at different years. The SPI-1 variation was found to vary from severe drought (−2) to severe wet (+2). The drought was more pronounced in September and October than the other months in the later years. By analyzing the intra-annual variation of the mean SPI-1 values, we find that 1975-2015 had near normal and moderate (1≥ SPI-1 ≤ 1.5) wet years. The years 1985The years , 1986The years , 1988The years , and 1997 were moderate wet years with average SPI-1 values varying from +1.31 to +1.91. The other years were near-normal years with SPI-1 values ranging from +0.08 to +0.91. SPI-1 variation was found to vary from severe drought (−2) to severe wet (+2). The drought was more pronounced in September and October than the other months in the later years. By analyzing the intra-annual variation of the mean SPI-1 values, we find that 1975-2015 had near normal and moderate (1≥ SPI-1 ≤ 1.5) wet years. The years 1985The years , 1986The years , 1988The years , and 1997 were moderate wet years with average SPI-1 values varying from +1.31 to +1.91. The other years were near-normal years with SPI-1 values ranging from +0.08 to +0.91. Most of the mean SPI-1 values for the annual times vary between 0 and 1. The minimum and maximum SPI-1 values were -0.86 (near normal) and +2.52 (extreme wet). The near-normal droughts have been observed periodically since the 1990s (1993,1995,2002,2004,2011,2013), implying the recurrence of drought events in the future. In contrast, the years with extreme wet events are fewer (1985, 1986, 1988, 1989, 1995, and 1997). The number of extreme wet events has not been observed since 1997, which implies reduced precipitation. The year 1997 was reported to be an extremely wet year, whereas, on average, none of the years was reported to have an extremely dry year. Therefore, most of the historical years were observed to have near-normal years. Most of the mean SPI-1 values for the annual times vary between 0 and 1. The minimum and maximum SPI-1 values were -0.86 (near normal) and +2.52 (extreme wet). The near-normal droughts have been observed periodically since the 1990s (1993,1995,2002,2004,2011,2013), implying the recurrence of drought events in the future. In contrast, the years with extreme wet events are fewer (1985, 1986, 1988, 1989, 1995, and 1997). The number of extreme wet events has not been observed since 1997, which implies reduced precipitation. The year 1997 was reported to be an extremely wet year, whereas, on average, none of the years was reported to have an extremely dry year. Therefore, most of the historical years were observed to have near-normal years.

Temporal Analysis
Temporal analysis of the rainfall anomaly index (RAI) over Nepal was carried out for 41 years (1975-2015), as shown in Figure 9, portraying the drought and rainy years. It can visualize the periods where drought and rainfall events occur throughout the temporal period. Most of the year, RAI is positive (>0), whereas some years experience a dry year before 2005. Most of the years after 2005 experienced dry years, whereas only two years (2008 and 2011) were humid. The positive values observed in Figure 9 represent rainy or wet years, and the negative values represent the dry years with different degrees of intensity. The occurrence of 23 years with a positive RAI, varying from extremely wet to humid, and 18 years with a negative RAI, ranging between very dry and dry, was observed. In other words, the drought years and rainy years are comparatively equal in numbers for the 41-year period. On analyzing the 41 years, drought years lasted for three years (1975)(1976)(1977)  wet years, and the negative values represent the dry years with different degrees of intensity. The occurrence of 23 years with a positive RAI, varying from extremely wet to humid, and 18 years with a negative RAI, ranging between very dry and dry, was observed. In other words, the drought years and rainy years are comparatively equal in numbers for the 41-year period. On analyzing the 41 years, drought years lasted for three years (1975)(1976)(1977) before 2005, whereas the drought year increased to a maximum of 4 years in the later period of analysis, i.e., 2012-2015. However, none of the years had an extremely dry year, i.e., RAI > −4. Further, the state of dryness increased with time, which stimulates water scarcity after 2005.  Figure 10. June to September is treated as rainy months, and the remaining months are exhibited as drought months for both drought and rainy years. The monthly distribution of RAI shows the month of June and September experienced drought months compared to the remaining months. For most of the drought month for both drought and rainy years, dryness intensity does not exceed the "very dry" condition. For example, the driest years in the historical series ( Figure 10) presented precipitation above average only in July and August for 1993 and 2015, whereas it continued for June and September for 1993. However, the drought year 2015 experienced dryness in June and September (intensity label as dry). In this light, it is possible to infer that there are months when RAI tends to be negative in very dry years, even in the rainy season.  Figure 10. June to September is treated as rainy months, and the remaining months are exhibited as drought months for both drought and rainy years. The monthly distribution of RAI shows the month of June and September experienced drought months compared to the remaining months. For most of the drought month for both drought and rainy years, dryness intensity does not exceed the "very dry" condition. For example, the driest years in the historical series ( Figure 10) presented precipitation above average only in July and August for 1993 and 2015, whereas it continued for June and September for 1993. However, the drought year 2015 experienced dryness in June and September (intensity label as dry). In this light, it is possible to infer that there are months when RAI tends to be negative in very dry years, even in the rainy season.

Spatial Analysis
The spatial distribution of RAI over the country in five years for the 41 years is shown in Figure 11. For 1985-2000, RAI intensity shows small traces of dryness in the central and western parts of the country. Beyond 2000, the intensity has shifted to dry compared to 1985-2000, with the eastern region more susceptible to intense dryness, i.e., extremely dry. From 1975 to 1980, the intensity of RAI had a mixed nature, partially dry and humid. This might be because some stations operated late or could not record the rainfall data correctly during the early period of installation. The temporal distribution of RAI showed a similar trend of increasing drought after 2005 as indicated by SPI. However, the spatial distribution showed that the magnitude of the drought severity posed by RAI is high compared to those offered by SPI. The result further revealed that the eastern region had more pronounced drought compared to the country's western region.

Drought Characteristics: Severity, Frequency, and Duration of Drought Events
The frequency of extreme drought events (SPI-1 < −1.5) is observed for different temporal scales ( Figure 12A -1985-1995, 1996-2005, 2006-2015 are also analyzed. The analysis shows that the number of events is higher than in earlier decades. Many such events are observed from May to September, implying the shifting of rainy days exacerbated by the changing climate at regional and local levels.

Spatial Analysis
The spatial distribution of RAI over the country in five years for the 41 years is shown in Figure 11. For 1985-2000, RAI intensity shows small traces of dryness in the central and western parts of the country. Beyond 2000, the intensity has shifted to dry compared to 1985-2000, with the eastern region more susceptible to intense dryness, i.e., extremely dry. From 1975 to 1980, the intensity of RAI had a mixed nature, partially dry and humid. This might be because some stations operated late or could not record the rainfall data correctly during the early period of installation. The temporal distribution of RAI showed a similar trend of increasing drought after 2005 as indicated by SPI. However, the spatial distribution showed that the magnitude of the drought severity posed by RAI is high compared to those offered by SPI. The result further revealed that the eastern region had more pronounced drought compared to the country's western region. In contrast, during April, the number of extreme events was higher during the 1985-1995 decade as some parts of the country exhibited less rainfall in 1989. The seasonal dispersion of the extreme drought events showed a similar tendency. The number of drought events increased by 2.28 times in 2006-2015 compared to 1996-2005 during JJAS months, whereas the annual events multiplied by 1.77 during the same period. The inclination of drought events in the later decades (or later years) signifies the risk of drier days.
The drought characteristics, frequency, and duration of SPI events were analyzed ( Figure 12). The number of drought events is found to increase in both annual and decadal periods. The number of decadal drought events shows the increasing rate of drought events during the 1990s, which started reducing during 1990-2000 and again began to increase in 2000-2010. The number of drought events in the first five years of the 2010-2020 decade is more than 50, close to the number of drought events in 2000-2010. The increase in the number of warm years and decrease in the number of wet days might be the potential cause for the rise in the number of drought events in later decades. The maximum number of annual and decadal SPI-1 drought events less than −1.5 in the study region for the period 1975-2015 are 150 and 500, respectively. Similarly, the annual and decadal mean SPI-1 events (<−1.5) are 42.17 and 6.14 events per year, respectively.   Figure 13 shows the paddy crop's temporal distribution at selected districts in each province and the average paddy yield in seven provinces. Wesee that the paddy yield started to drop after 2005, which might be due to various factors such as rainfall deficit, lack of fertilizers, lack of cultivable land for crop yield, and improper land management practices. In contrast, the paddy yield's temporal distributionincreased in selected districts of Gandaki province and Lumbini. In addition, the paddy yield was reduced to 1. the number of extreme events. The monthly and seasonal variations for different decades-1985-1995, 1996-2005, 2006-2015 are also analyzed. The analysis shows that the number of events is higher than in earlier decades. Many such events are observed from May to September, implying the shifting of rainy days exacerbated by the changing climate at regional and local levels. In contrast, during April, the number of extreme events was higher during the 1985-1995 decade as some parts of the country exhibited less rainfall in 1989. The seasonal dispersion of the extreme drought events showed a similar tendency. The number of drought events increased by 2.28 times in 2006-2015 compared to 1996-2005 during JJAS months, whereas the annual events multiplied by 1.77 during the same period. The inclination of drought events in the later decades (or later years) signifies the risk of drier days.
The drought characteristics, frequency, and duration of SPI events were analyzed ( Figure 12). The number of drought events is found to increase in both annual and decadal periods. The number of decadal drought events shows the increasing rate of drought events during the 1990s, which started reducing during 1990-2000 and again began to

Impact on Crop Production
The annual crop production for major crops in selected districts of the country is shown in Figure 14 and Supplementary Materials, Figure S5. The districts are so selected that they represent the entire country. The temporal distribution of crop yield shows a reduced trend Earth 2022, 3 426 after 2007 for rice, maize, millet, and barley. During 1990-2000, rice production was reduced in the Bhaktapur district, attenuated by reduced rainfall and increasing built-up area. A rising crop production (yield) trend for major crops is observed with some noticeable spikes. An elevated crop yield is due to multiple factors, including improved seeds, fertilizers, and agricultural equipment. The establishment of irrigation facilities (overland or sprinkler/drip) in recent times helps farmers combat the negative consequences of climate anomalies, particularly drought. A negative spike in rice production in Bhaktapur (Bagmati province) was observed in 2009 during the latter part of the analysis period. 2009 is characterized by a moderate drought from both the drought indices. Similarly, Surkhet and Dolpa (Karnali province) significantly reduced rice production in 2010. They are located in the western region, when 2010 was characterized as a drought year in the western parts of Nepal. In particular, Surkhet had reduced the production of barley, maize, and millet in 2010.
increase in 2000-2010. The number of drought events in the first five years of the 2010-2020 decade is more than 50, close to the number of drought events in 2000-2010. The increase in the number of warm years and decrease in the number of wet days might be the potential cause for the rise in the number of drought events in later decades. The maximum number of annual and decadal SPI-1 drought events less than −1.5 in the study region for the period 1975-2015 are 150 and 500, respectively. Similarly, the annual and decadal mean SPI-1 events (<−1.5) are 42.17 and 6.14 events per year, respectively. Figure 13 shows the paddy crop's temporal distribution at selected districts in each province and the average paddy yield in seven provinces. Wesee that the paddy yield started to drop after 2005, which might be due to various factors such as rainfall deficit, lack of fertilizers, lack of cultivable land for crop yield, and improper land management practices. In contrast, the paddy yield's temporal distributionincreased in selected districts of Gandaki province and Lumbini. In addition, the paddy yield was reduced to 1.  Figure 15 shows the temporal distribution of considered drought indices (SPI-1 and RAI) over the studied period in different districts of each province. The left top is Province 1, and the left bottom is Sudurpaschim province. An apparent variation of SPI-1 and the crop yield is observed for the considered districts, whereas a distinct variation with the RAI could not be observed. For example, Dhading district in Bagmati province increased the paddy yield during 1991-2006 (wetter side), whereas a reduction was observed during 2007-2011 (drier side). A similar trend is observed for other districts and crops across the country. Lower regions of the country are flatter in topography (like Sarlahi), where rice is cultivated two times a year and other crops have multiple harvests. They are the country's food granary, and irrigation facilities are readily available compared to other physiographical regions of Nepal [44]. Throughout the study period, a steady rise in crop production could be found in Jhapa. A considerable spike could be found in barley production in Jhapa in the year 1998. Two noticeable spikes are observed in 1995 and 2006. It is challenging to find the causes behind these spikes because these are affected by multiple factors, including natural and human interventions. However, a clear indication of a regular or routine break could be seen during a drought period. A precise analysis on an event basis at a different location demands digging/dissecting the underlying causes. Our study, though general, might be useful for sufficient water resources and irrigation facilities management. Considering the likely drought effects, a proper selection of crops and their rotations can be assessed for different locations at different times. Inter-basin water transfer is another viable option for managing water balance between deficit and surplus Interestingly, reduced production of rice was compensated for by an elevated maize production in Dolpa in 2010. The result suggests the switching of the crop from rice to maize in 2010. Short growing seasons and low temperatures limit crop growth in high mountains like Dolpa [43]. Maize, barley, and millet are grown at higher elevations, limiting the crop in a year. They are mainly rain-fed, and irrigation facilities are limited. Cold-tolerant rice is grown in such areas; for instance, Jumli Marshi (traditional rice varieties) is produced on the mountain slopes of Jumla (2830 m asl), the highest altitude rice-growing location in the world [62]. In mid-hills like Surkhet, crop cultivation is practised in both summer and winter. Thus, crop-shifting in mid-hills is not significantly apparent as it is in high mountain districts. Figure 15 shows the temporal distribution of considered drought indices (SPI-1 and RAI) over the studied period in different districts of each province. The left top is Province 1, and the left bottom is Sudurpaschim province. An apparent variation of SPI-1 and the crop yield is observed for the considered districts, whereas a distinct variation with the RAI could not be observed. For example, Dhading district in Bagmati province increased the paddy yield during 1991-2006 (wetter side), whereas a reduction was observed during 2007-2011 (drier side). A similar trend is observed for other districts and crops across the country. Lower regions of the country are flatter in topography (like Sarlahi), where rice is cultivated two times a year and other crops have multiple harvests. They are the country's food granary, and irrigation facilities are readily available compared to other physiographical regions of Nepal [44]. Throughout the study period, a steady rise in crop production could be found in Jhapa. A considerable spike could be found in barley It is challenging to find the causes behind these spikes because these are affected by multiple factors, including natural and human interventions. However, a clear indication of a regular or routine break could be seen during a drought period. A precise analysis on an event basis at a different location demands digging/dissecting the underlying causes. Our study, though general, might be useful for sufficient water resources and irrigation facilities management. Considering the likely drought effects, a proper selection of crops and their rotations can be assessed for different locations at different times. Inter-basin water transfer is another viable option for managing water balance between deficit and surplus areas. A naturally occurring drought, if well understood, can be managed, and the negative consequences can be mitigated.

Observed Crop Data
Earth 2022, 3,21 areas. A naturally occurring drought, if well understood, can be managed, and the negative consequences can be mitigated.

Conclusions
This study used two drought indices, namely SPI and RAI, to map Nepal's metereological drought condition both temporally and spatially. Although SPI has been extensively applied for drought mapping, RAI has been less explored in earlier studies. Thus, this study used 41 years (1975-2015) of 124 station precipitation data distributed heterogeneously over Nepal to analyze both SPI and RAI. Analyzing all the meteorological stations, we found that the average annual SPI-1 value ranged from −0.20 for most of the year to +2. 21 (1997) at station 101. Similarly, at station 308, the annual average SPI-1 value ranged from −0.20 to +2.31. The temporal distribution of SPI shows the increasing tendency of drought events beyond 2009 for SPI-1, 2006 for SPI-3, and 2004 for SPI-6 and SPI-12. The temporal SPI result concluded that the drought situation has risen in Nepal after 2004, which might be the consequence of anthropogenic activities. The temporal distribution showed an increasing trend of drought in recent years. Hence, research on drought impacted by climate change seems necessary to project the future temporal drought. The spatial distribution of SPI showed that eastern regions are subjected to drought under different SPIs. The result from the observed rainfall exhibits the near-normal SPI conditions in most of the regions in Nepal. The result concludes that the spatial distribution of the drought is not very severe across the country. However, if the trend continues and a proper mitigation approach is not adopted, the situation might lean towards severity and may cause an impact on water resource management and associated fields. The spatial and temporal results of RAI showed a similar trend to SPI in the drought level. The wet months (June-September) show the positive RAI anomaly, whereas the dry months (October-May) show the negative RAI anomaly. The RAI spatial distribution showed the drier condition towards the eastern region, whereas the western region showed wetter conditions similar to those portrayed by SPI. Both SPI and RAI can equally explain the meteorological dryness/wetness. The variation in the crop yield had been observed through the variation of SPI events. With the increase in the wetter years, the crop yield increased in the Dhading district, wheras the crop yield reduction was found in the drier years. We finally conclude that the meteorological variation directly modulates the crop yield.