Abstract

Climate change is the biggest environmental challenge that plays out through changing intensity, duration, and frequency of extreme events. To fulfill the research gap in understanding and quantifying the recent changes in precipitation extremes over the Narayani river basin of central Nepal, the long term daily precipitation data from 1980 to 2018 were run in ClimPACT2 an R software package to calculate ET-SCI extreme precipitation indices. In this study physically relevant 14indices obtained from 23 stations were examined for their spatial and temporal variation. Before the calculation of indices different scientific tools and programs were used to check the data quality and homogeneity. The results suggest that the variations of extreme indices throughout the study area are quite different from that of seasonal and annual patterns to some extent. The monsoonal precipitation was mostly concentrated in the central part of the basin within the Middle Mountain region (Lumle and its surroundings). Especially the lowlands (Terai and Siwaliks) and including some parts of middle mountains the precipitation intensity-based indices like as, percentile indices (R95p) and absolute indices (RX1day, RX3day, RX5day) were in the increasing trends, but the frequency of precipitation like threshold indices (R1mm, R10mm, R20mm) along with the duration of precipitation seemed to be decreasing. This implies that the lowlands regions bringing about rainfall related hazards like floods and soil degradation with inundation and may cause possible impact on agriculture and livelihood due to intense rainfall and prolongation of dry spells with the weakeing of rainfall duration (days/year). However, the light to moderate precipitation and associated days over the high altitude and that could be the possible cause of landslides. This study also highlights the suggestion that there may be a possible impact on agriculture facilities, food security, and water scarcity in the eastern part of the basin due to the significant decreasing trend of annual total wet days precipitation (PRCPTOT).

1. Introduction

Climate change is the biggest environmental challenge attracting global attention and takes the planet into ‘uncharted territory’ with the announcement of 2016 as the warmest year ever ¹. As a result, heat waves, cold waves, droughts, and floods which used to be occasional events are becoming more regular than in the past ². A changing climate leads to change in the frequency, intensity, spatial extent, duration and timing of weather and climate events and produces the unexpected climatic ³. Increases in extreme climate events such as heat waves, heavy rain, drought, and associated wildfires have significant impacts on human health, livelihoods, assets, and ecosystems ⁴ and therefore, the more extreme events, the more likely it is to cause societal or environmental damage ⁵.

The responses of precipitation extremes (heavy precipitation events) to climate change have been the subject of extensive study because of the potential impacts on human society and ecosystems ³. The study ^{6, 7} highlights the importance of advancing the understanding of the character of precipitation based on recent warming. The authors pointed out the increasing precipitation extremes like heavy precipitation being somehow associated with recent warming. Especially the heavy precipitation event depends on the available atmospheric moisture and is governed by the Clausius-Clapeyron equation.

In the meantime, several research papers, and reports published results on the analysis of precipitation extremes on global or continental or local scales to evaluate the changes in precipitation characteristics. Reference ⁸ brought the most up-to-date and comprehensive global picture of the observed trends in extreme in precipitation along with temperature based on the high-quality station data supplied by numerous scientists worldwide and found a widespread and significant increase in heavy precipitation indices. The findings were accepted by the Special Report on Extreme Events of Intergovernmental Panel on Climate Change (IPCC) ³ and also valid with the findings of Fifth Assessment Report fifth ⁴, which highlighted the general increasing trends of extreme precipitation events in most parts of the earth. Reference ⁹ also found that on a global scale the intensity, frequency, and duration of extreme precipitation is increasing on average. The Fifth assessment report of the Intergovernmental Panel on Climate Change ⁴ also stated that extended intervals of monsoon failures and dry spells have struck India and Southeastern Asia, in the last few years leading to prolonged and intense droughts, which are a recurring feature of Holocene paleoclimate, also, the frequency of heavy precipitation events are increasing while light rain events are decreasing. Tree ring based precipitation studies also showed the frequency of drought event is increasing in the recent decades in the western Himalaya ¹⁰ as well as the northern and southern slope of central the Himalayan region ^{11, 12, 13}.

The comprehensive statistical analysis of precipitation observations during the South Asia summer monsoon season showed the significant increases in the frequency of dry spells (low rainfall) and intensity of wet spells (heavy rainfall) and decreases in the intensity of dry spells in South Asia ¹⁴. In similar studies from the Hindu Kush Himalayan (HKH) showed heterogeneity in the amount of precipitation throughout the region, however regional average annual intense precipitation amount and consecutive wet-days, regional average of RX1DAY, RX3DAY, and RX5DAY is increasing ¹⁵. However, trends of consecutive dry-days were mostly opposite to those of consecutive wet-days.

Nepal, situated in the central Himalayan with steep to pography is more susceptible to the developments of heavy rainfall and drought events however as compared to the rest of South Asia studies on the observed precipitation extremes over Nepal are rare and sporadic. Before the last decade the Disaster Risk Reduction (DDR) strategy and policies were given very low priority in Nepal. Nepal’s first disaster report ¹⁶ shows that climate-related disasters accounted for almost one-fourth of the deaths and three-fourth of economic losses. Furthermore, Reference ¹⁷ also reported that Nepal is more vulnerable country due to extreme precipitation related hydro-meteorological disaster such as; flash floods and landslides. The research reports ^{18, 19} highlights the rising drought events as a consequence of changes in precipitation characteristics across the southern Himalaya (Nepal) and which could be the potential cause of the recent failure of the agriculture system in some parts of the country. For instance, the characteristics of the precipitation and its distribution throughout the country are unique. The annual precipitation distribution over Nepal possesseda similar pattern to the monsoonal distribution and delineated all the precipitation pocket zones (Lumle, Gumthang, and Num) as well as the low precipitation zones on the leeward side of high mountains (Mustang, Manang, Dolpa,and Everest) ²⁰. On average all the seasonal and annual precipitation trends of Nepal were insignificantly decreasing with the highest decreasing trend in the post-monsoon season than other seasons ^{21, 22}. Except for the Terai region of Nepal, the annual precipitation decreased at other physiographic zones with the highest negative trend in the High Mountain region. Reference ²³ used 26 stations to calculate the extreme precipitation across Nepal for the period of 1971-2006 and found that the increasing trend in the extreme precipitation events at most of the stations. Similarly, by the analysis of precipitation datasets from 210 stations for the period 1970-2012 across Nepal found the decreasing trend in the annual precipitation, wet days and high-intensity-related precipitation extremes (R95, RX1day), together with an increasing trend in the consecutive dry days over the central and eastern middle mountains ²⁰. The study ²⁴ analyzed the precipitation events and found that the number of extreme events per station is neither significantly increasing nor decreasing between 1971 and 2010 across Nepal. The analysis of spatio-temporal variability of extreme precipitation in Nepal by ²⁵ found that peak annual precipitation between elevation 2,000 and 3,500 m above sea level (masl), but 1-day extreme precipitation peaks were found at the lower elevation in the southern foothills with its highest intensity in a central region of the country. The research output of ²⁶ on the northern and southern slopes of the central Himalayas from 1975 to 2009 found that, the increased trend of all the indices except the maximum 1-day precipitation and simple daily precipitation intensity indices on the southern slopes. A recent study in the Roshi river basin found the increasing trend of consecutive dry days but heavy rainfall days (RX5day) and precipitation intensity on wet days were in a decreasing trend. The stream flow pattern was severely affected by these changes in extreme precipitation events ²⁷. Over the whole Koshi river basin, an increase in total annual rainfall and rainfall intensity was determined with no clear long-term linear trend while analyzing the recent 36 years of records ²⁸. Furthermore, the number of consecutive dry days increased at almost all the stations of the Koshi basin. The recent two decades of study in the same river basin ²⁹ found a solo increasing trend of consecutive dry days but, the intensity indices of extreme precipitation was declined across Eastern Nepal. A similar type of study in the West Rapti river basin shows that the number of consecutive dry days and consecutive wet days were in the upward trend in at most of the stations and the extremely wet days, the annual total precipitation and maximum1-day precipitation trends were declined ³⁰. While in the Kali-Gandaki river basin except for the total annual precipitation index all the indices were in increasing trend till 2007 ³¹.

Based on the above-mentioned research findings, there has not been a recent and specific assessment of observed changes in climate extremes over the Narayani River basin of central Nepal, where meteorological observing systems are among the best in Nepal. The study ³² analyzed the spatial and temporal trend of precipitation and leaves a gap in understanding and quantifying the recent precipitation extremes over the basin. Therefore, the present study will evaluate the details of precipitation extremes changing in terms of spatial and trend features including intensity, frequency, and duration by the nature of occurrence from the basin and fulfill the regional gap in the climatic extremes focusing on precipitation.

2. Study Area

The Narayani river basin is a transboundary river basin, also known as the Gandaki river basin. It covers the three countries China in the northern part, Nepal in the middle part and India in the southern part, having the total drainage area of 46300km², among which 72% lies in Nepal alone, remaining 18% and 10% lies in India and China respectively (Figure 1). The basin has a high topographic variation with elevation ranging from 60 masl in the south to higher than 8000 masl in the north where it passes through the high Himalayas, which contain the Dhaulagiri (8167 masl) and Annapurna (8091 masl) peaks. Narayani River includes seven tributaries namely, Kali Gandaki, Budhi Gandaki, Trishuli, Marsyandi, Daraudi, Seti, and Madi. It contains 1710 glaciers with an area of 2285km^{2 33}.

According to Nepal Electricity Authority ³⁴ the Narayani river system contributes about 50% of the total hydropower production of the country. Furthermore, it has a high potential for investment in hydropower development as well as irrigation in the downstream regions. The basin is also characterized by high poverty levels and low socio-economic indicators. Avalanches, landslides, glacial lake outburst floods (GLOFs) and other floods triggered by heavy rainfall, and drought are the major climatic issues in the upstream area (>2,500masl) due to the abrupt rise of topography, extremely rugged terrain, steep slopes, and deeply cut valleys. The Siwalik zone of the basin has dry and unconsolidated soil materials that are highly prone to erosion and landslides, while the Terai area is characterized by vast flood plains with deciduous forest and agriculture and subject in many places to annual inundation. A study carried out by the International Centre for Integrated Mountain Development (ICIMOD), Nepal ³⁵ also highlights the recent climate change and its toll on agriculture, water resources, tourism and biodiversity in many parts of the basin.

The climate in the basin varies from humid tropical in the low altitude of the southern part to polar tundra in the higher altitude northern belt ³⁶. The basin experienced four seasons, namely, winter (December-February), pre-monsoon (March-May), monsoon (June-September), and post-monsoon (October-November). Winter is a cold season with the precipitation over the basin is received by the western disturbances. Pre-monsoon is characterized by hot, dry and westerly windy weather with mostly localized precipitation in a narrow band. The monsoon is characterized by moist southeasterly monsoonal winds coming from the Bay of Bengal and occasionally from the Arabian Sea. Post-monsoon is a well-known dry season with the driest month, November. The spatial distribution of precipitation in the basin is highly heterogeneous (Figure 2). For instance, the annual precipitation varies from less than 270mm for the driest regions (Mustang, Manangand Dolpa, located at the leeward-side north of the Annapurna) to above 5401mm in and around the Lumle region. Precipitation patterns in the basin are directly associated with the summer monsoon, with about 78% of the annual precipitation falling between June and September.

3. Materials and Methods

The observed daily precipitation dataset of 74 meteorological stations for the period of 1980-2018 was collected from the Department Of Hydrology and Meteorology (DHM), Nepal. Within this study period the stations that feature data gaps for more than (a) three days in a month; (b) two weeks within a year, (c) three consecutive years or (d) six years in total were omitted from the analysis. All the selected stations within the study area are presented in Figure 1. Prior to the indices calculation, the errors were minimized through the quality control (QC) procedures using ClimPACT2 ³⁷ and RHtest ³⁸ was used to perform the homogeneity test. The monthly time-series data of each station were employed for RHtest. The stations used in this study were also tested by ²⁰ for homogeneity using the same toolkit. After the quality control and homogeneity test, the daily precipitation data from 23 stations (Table 1) were subjected to the indices calculation with the reference period of 1981-2010 under the ClimPACT2 toolkit.

ClimPACT2 gives both core and non-core Expert Team on Sector-Specific Climate Indices (ET-SCI). For our analysis 14, ET-SCI precipitation indices were selectedas presented in Table 2. The non-parametric Mann-Kendall trend test ^{39, 40} was used to assess the significance of the trend in the time series and Sen’s slope estimate ⁴¹was used to determine the true slope. The Excel template is generally known as MAKESENS ⁴² was used for Mann-Kendall test and Sen’s slope estimation. The trends were considered at 5% level of significance. For the spatial analysis of seasonal and extreme precipitation indices over the basin, we use the most frequently used interpolation method known as the ordinary kriging method in ArcGIS platform ⁴³ to interpolate the values in un-sampled location.

For better understanding all the indices (Table 2) were grouped under five broad categories:

i. Percentile-based indices (R95p, R95p TOT, R99p, R99p TOT)

ii. Absolute indices (RX1day, RX3day, RX5day)

iii. Threshold indices (R10mm, R20mm, R1mm)

iv. Duration indices (CDD, CWD)

v. Other indices (PRCPTOT, SDII).

4. Result and Discussion

The average seasonal distribution of precipitation over the basin with the precipitation trends for the long term 1980 to 2018 are presented in Figure 2 and Figure 3 respectively. The annual average precipitation distribution over the Narayani river basin was found to be in the range of 270 mm to 5401 mm. The Northwest part of the basin around Jomsom meteorological station on the leeward side of the Annapurna-Machhapuchare range received minimum precipitation but the Lumle area (central part of the basin) that lies on the windward side and thus received maximum precipitation. This result is also supported by some earlier studies ^{20, 21, 22, 32, 36}. The middle mountain receives the maximum amount of precipitation whereas the High Himalaya receives the minimum amount of rainfall due to the altitudinal dependency of precipitation ^{20, 21, 22}. The seasonal contribution of average precipitation over the basin from winter, pre-monsoon, monsoon and post-monsoon, are 3.125%, 14.26%, 78.9 and 3.68%, respectively. The winter precipitation shows distinct features with the gradual increase of precipitation up to the HM (<2500 masl), while for other seasons including annual patterns the highest value limited within the MM. The maximum value of annual average precipitation found as 2800 mm in Middle Mountain and amongst which 80% is contributed by the summer monsoon, 14% by pre-monsoon and the other two seasons with 3% each.

It is interesting to note that the winter precipitation at all the stations below 2300 m elevation was found to be in decreasing trends (Figure 3) with less than 0 to -1.68mm y^-1 with the average trend of -0.83mm y^-1 and above such elevation, it was increased by 0.54 mm y^-1. The winter precipitation over the basin was decreased by -0.71 mm y^-1 on average. The winter precipitation usually plays a crucial role in supporting the winter agriculture system of the basin; due to the decreasing trend of winter precipitation, there could be a decreased in winter production.

Jomsom and Late belong to the high elevation meteorological stations of this study area and where seasonal increasing trends of precipitation were observed except in post-monsoon. Except for pre-monsoon, the precipitation trends were decreased for the summer monsoon and post-monsoon season with the rate of -2.16 mm y^-1 and -0.18 mm y^-1 respectively. Even the pre-monsoon shows positive trend 0.93 mm y^-1, the average annual precipitation over the basin was in an insignificant decreasing trend of -2.28 mm y^-1. This annual declination rate of precipitation found to be greater than the annual average precipitation trend over Nepal ²¹. Amongst the 23 meteorological stations, 15 shows the decreasing trends on average of -5.53 mm y^-1 for the long term annual records with maximum negative and positive trends of -13.85 mm y^-1 and 15.89 mm y^-1 observed within a 30 km areal distance of Ridi Bazar in Ruru-chhetra a holly place and Musikot Gulmi respectively. This finding describes the orographic effect on monsoon.

The results of extreme precipitation indices including ever-recorded one-day extreme precipitation (evRX1day) over the Narayaniriver basin obtained for the period of 1980 to 2018 are presented through the spatial interpolation.

Percentile indices monitoring the intensity of precipitation extremes are widely used to assess the recent changes in precipitation extremes and are important for a number of potential impacts, including floods, soil erosion, water resources, agriculture, and overall ecosystems. The percentile-based extreme precipitation indices are: very wet days (R95p) and extremely wet days (R99p) and these indices represent the amount of precipitation falling above the 95th and 99th percentile of records respectively. The R95p and R99p indices were found to be in the range of 60-1181 mm and 23 to 337 mm each year with respective rainfall. Looking at the spatial patterns of both the R95p and R99p (Figure 4) the rainfall was concentrated in the central part of the basin or the middle mountain region and minimum in the northwestern part with slight variation in areal coverage of maximum R99p towards the southeast extend than that of R95p.

The absolute precipitation indices are maximum 1-day precipitation amount (RX1day), maximum 3 days precipitation amount (RX3day) and maximum 5 days precipitation amount (RX5day). The spatial distribution of RX1day, RX3day and RX5day observed in a similar pattern with a slight variation of intensity with the range of 35 to 209 mm, 51 to 359 mm and 56 to 486 mm respectively each year as presented in Figure 5. It is clear that the absolute indices (RX1day, RX3day, and RX5day) were intense around the central part of the basin area along with the south and southeast portion, while their intensity decreased throughout the HH region. The ever recorded 1-day maximum precipitation in Figure 5 (d) shows that maximum 1-day precipitation was found to be more intense toward the southern belt of the basin with a peak value of 516mm in the record. It indicates that the lowland region in the southern part of the basin is much more vulnerable to flood-related disasters like; loss of life and properties. This result is strongly supported by some previous research works ^{17, 20}.

The threshold indices here represent the number of days at which precipitation value falls below or above the fixed threshold such as number of wet days (WD or R1mm), number of heavy precipitation days (R10mm) and the number of very heavy precipitation days (R20). The spatial distributions of such indices are given in Figure 5. The annual frequency of heavy precipitation day (R10mm) and very heavy precipitation day (R20mm) indices have a range of 46-182 days and 2-80 days each year respectively. The central part of the basin (or Middle Mountain) around Lumle and Pokhara said to be the pocket area concerning the concentration of all these indices but its south and northwest parts seem to have a lower frequency of heavy rainfall. The increase in the frequency of threshold indices could possess some serious impact on the Middle Mountain region, landslide as a major and agriculture losses are secondary.

The precipitation duration indices are representative of the periods of excessive wetness or dryness of the basin. It includes the consecutive dry days (CDD) and consecutive wet days (CWD). As shown in Figure 6, the CDD for the Narayani river basin has a frequency range of 43-86 days each year. It is clear that in the southern part including Terai and Siwalik region of the basin numbers of consecutive dry days were highest and seem to be decreased northward by the Middle Mountain (MM). It indicates how the agricultural hub of the Narayani river basin gets dominated by maximum dry days. It is interesting to note down that some part of High Himalaya (HH) in western part around Jomsom possesses high frequency of CDD than the central and eastern belt. The study of the nearby river basin in eastern Nepal also shows that the HH area was comparatively drier than any other parts ²⁹. The range of consecutive wet days (CWD) index was 4-48 days each year (Figure 7). The frequency distribution shows the opposite in nature to the CDD index with the greatest frequency toward the center of the basin in the MM region and decreased north and southward.

In this category, the annual precipitation total (PRCPTOT) and simple daily intensity index (SDII) have been analyzed. The PRCPTOT index was found in the range of 265-5394mm each year. According to Figure 8 the maximum concentration of PRCPTOT found at MM as a center for the Northeast curve of low-value are asemerged from Northwest. The simple daily intensity index (SDII) has a range of 6 to 32 mm/day each year and its spatial distribution shows similar patterns as R95p and R99p, which proves that the high intensity of precipitation about the center part of the basin. From these results, it is obvious that water resources and hydrology of the southern belt of the basin are under stress.

Various rainfall related extreme indices were estimated for the 23 meteorological stations in the Narayani river basin for the period of the recent 4 decades; the trends of such indices are summarized in Table 2. The percentage of stations with different trend features of precipitation indices is presented in Figure 9, where 100% of station distribution covers all 23 selected stations.

The annual total precipitation on very wet days (R95p) for the recent 4 decades exhibits the trends of -0.05mm y^-1. This is consistent with the result reported by ^{29, 31}. R95p found to be in negative trend at 60.9% of total stations and including 4.3% of station represented for significant positive trend (4.11 mm y^-1) 39.1% have an insignificant negative trend. The R95p TOT index was non-significant with both positive (60.9% station) and negative trends (39.1% stations). Trend analysis for R99p was detected, but no trends were found for the majority of the stations due to large inter-annual variability that resulted in zero values for R99p. They were therefore omitted from the analysis. The precipitation fraction due to very wet days (R95p TOT) found to be increasing by 0.02% per year. Which may affect agriculture, food security, water resources and hydrology of the basin.

On average over the Narayani basin, absolute indices (RX1day, RX3day, and RX5day) were in increasing trend with the intensity of 0.11 mm y^-1, 0.16 mm y^-1 and 0.04 mm y^-1 respectively. The result is consistent with the study of ^{15, 23, 28}. RX1day shows an increasing trend at 47.8% of total stations but, for RX3days and RX5days more than 50% stations have positive trends. Among the stations having a positive trend of RX1day Pokhara airport (1.19 mm y^-1) was statistically significant. On the other hand, KhudiBazar station (-0.90 mm y^-1) belongs to a group of negative trend stations of RX1day. In the case of the maximum consecutive3-day and consecutive 5-day precipitation, Chapakot station possesses the highest negative trends of -1.83 mm per year and -2.43 mm y^-1 (statistically significant). Furthermore, Ridi bazar station is also has a significant negative trend (-1.6 mm y^-1) that belongs to the Siwalik range, other stations were statistically insignificant.

The average trends of absolute indices in Siwalik and High Himalaya range were decreasing, while in Terai and Middle Mountain those indices were in increasing trends, this is consistent with the results reported for whole Nepal ^{20, 21} and the adjoining Koshi river basin ^{28, 29}. As seen from Table 2, these absolute indices are mostly increased in low lands of the Terai region and thus bringing about rainfall related hazards like floods, soil degradation with inundation and may cause the possible impacts on agriculture and livelihood.

Out of 23 stations, 17 (73.9%). shows the decreasing trend of wet day (WD) index including 34.8% statistically significant with maximum negative trend of -1 days per year. This finding is partially supported by ^{20, 29}. In case of heavy precipitation days (R10mm) and very heavy precipitation days (R20mm) about 70% of stations exhibit maximum decreasing trends of -0.481 and -0.263 days per year respectively. This indicates that along with the normal precipitation days and heavy precipitation days within the basin are in decreasing trend at most of the stations. It should be noted that R1mm and R10mm increased more at high elevation (>2000 masl) with a maximum frequency of 0.86 and 0.03 days per year. While in the Middle mountain and Siwalik region they have the greatest declined rate of 0.5. In the case of the Terai region, R1mm seems to be decreased sharply with the maximum frequency of -1.0 days per year. These results suggest that in addition to Terai most of the Middle Mountain and Siwalik range moderate to heavy rainfall events are decreasing, while in High Mountain and Himalaya region the events seem to be increasing.

The consecutive wet days (CWD) index showed decreasing on an average of -0.02 days per year, but on the contrary, the consecutive dry days (CDD) were found to be increasing by 0.56 days per year for the last four decades within the study area. As presented in Figure 9, about 91.3% of stations represented for the increasing trend of consecutive dry days (CDD), out of which for 17.4% stations trends were only statistically significant. Remaining 8.7% of stations indicated for insignificant decreasing trend of CDD. In the case of consecutive wet days (CWD) 34.8% of stations were in increasing trend with only 8.7% of represented for significant trend thus, 65.2% of the stations have decreasing trends. These results for CDD and CWD are mostly satisfied by ^{15, 29}. It is worth noting that all the stations above the 2000m elevation showed anincreasing inthe CWD trend and all the stations in the low land area of Terai have a decreasing trend including many stations of Siwalik and MM. Furthermore, the frequency of consecutive dry days in HH station (Jomsom) reveals the highest negative trend (-0.25 days per year) and it was greater than 0.12 days per year below 2000m. These increasing trends of CDD and decreasing trends of CWD below 2000masl during the last 4 decades correlate much with the recent droughts events throughout the country identified by ^{18, 19}. The increasing trend of CDD and decreasing trend of CWD may affect the crop yield, cause wildfire by prolonged drought, and also the health-related problem by increasing of dust particles in the air within the basin.

The annual total wet day precipitation PRCPTOT decreased on average of -2.25 mm y^-1 with negative trends at 65% of total stations. PRCPTOT in High Himalaya exhibits a positive trend (1.94 mm y^-1) while in the Siwalik range it was decreased by -6.82 mm y^-1. One of the stations in Siwalik range (RidiBazar) showing a sharp significant decreasing trend of -13.86 mm y^-1. On the other hand, Musikot station in Middle Mountain displayed the highest significant increasing trend of 15.93 mm y^-1. In addition to the Jomsom as a high elevation station, the Late station of High Mountain also possesses a significant positive trend of 5.65 mm y^-1, this indicates the increasing pattern of precipitation in the northwestern part of the basin. There is a decreasing trend of PRCPTOT at all four stations situated in the eastern part of the basin within the range of -1.38 to -6.69 mm y^-1.

The simple daily intensity index (SDII) increased on average by 0.02 mm/day. SDII is increased at about 60% stations (13% statistically significant) and decreased at about 40% stations (4.3% statistically significant). The maximum increased trend was observed at Balewa Girwari station by 0.26 mm/day and a decreasing trend was observed to be -0.11 mm/day at Chapkot station. SDII was detected as a decreasing trend at all the stations above 2000 masl and increasing trends were observed at all stations of Terai regions (2 statistically significant). This result indicates that the precipitation intensity is rapidly increasing within the Terai than in other parts of the basin. The overall increasing pattern of SDII indicates that the precipitation is concentrated within a short period of time within the study area ^{20, 28}.

5. Conclusion

We have calculated and analyzed some physically relevant ET-SCI indices and average values based on the daily precipitation data of 23 stations from 1980 to 2018 over the Narayani river basin. The main findings can be summarized as:

• The spatial distribution of seasonal precipitations shows that the precipitation is highly concentrated at the central part of the basin (Lumle and its surroundings) and minimum precipitation was observed at the leeward side of the High Mountain (Jomsom and its surroundings). It received about 78.9 % of the total annual precipitation in monsoon season. Except for pre-monsoon the precipitation trends were in decreasing trends of -0.71 mm yr^-1, -2.28 mm yr^-1 and -0.18 mm yr^-1 respectively during winter, monsoon and post-monsoon season.

• The annual total precipitation on very wet days (R95p) decreased on average by -0.05 mm y^-1, but the precipitation fraction due to very wet days (R95pTOT) increased on average 0.02% per year. This indicates that the percentage of precipitation contribution of R95p throughout the basin was still positive.

• All the absolute indices RX1day, RX3day, and RX5day showed upward trends of 0.11 mm per year, 0.16 mm y^-1 and 0.04 mm y^-1 respectively. These indices were intense around the central part of the basin along with the south and southeast portion. It resembles the intense precipitation therein and thus bringing about rainfall related hazards like floods, soil erosion and more.

• The threshold indices over the basin exhibited decreasing trends, most notably, the normal precipitation days (R1mm) and heavy precipitation days (R10mm) were increased more at higher elevation (>2000 masl) while, most of the stations below 2000 masl were in decreasing trends. Which indicates that, middle mountain, Siwalik and Terai region leading towards drier conditions.

• More than 90% stations within a Narayani river basin represented for increasing trends (on average 0.56 days per year) of consecutive dry days (CDD) and on the contrary 65% of stations indicted for decreasing trends (on an average of -0.02 days per year) for the consecutive wet days(CWD). It is worth noting that all the stations above the 2000 m elevation showed increasing in the CWD trend and all the stations in the low land area of Terai have decreasing trends including many stations of Siwalik and MM. An increasing trend of CDD and the decreasing trend of CWD linked with the livelihood of the people in the area, which affects the agriculture activities through prolongation of the drought period.

• The simple daily intensity index (SDII) increased on average of 0.02 mm/day (at about 60% stations). SDII was detected as a decreasing trend at all the stations above 2000m and increasing trends were observed at all stations of Terai regions. But, mixed patterns were observed in Mountainous regions.

• PRCPTOT in High Himalaya exhibits a positive trend (1.94 mm y^-1) while in the Siwalik range it was decreased by -6.82 mm y^-1. There is a decreasing trend of PRCPTOT at all four stations (Hetauda, Chisapani-Gadi, Dhunibesi, Kakani) situated in the eastern part of the basin with the range of -1.38 to -6.69 mm y^-1. It suggests that there may be a possible impact on the overall agriculture facilities and food security including water scarcity in the eastern belt of the basin.

Acknowledgements

The authors like to thank Department of Hydrology and Meteorology (DHM), Government of Nepal forthe meteorological data. The corresponding author (BD) is thankful to Kathmandu Center for Research and Education, Chinese Academy of Sciences-Tribhuvan University, (KCRE) for the support.

Conflicts of Interests

The authors declare that there is no conflict of interests regarding this paper.

References