Despite the increased frequency of extreme climate events including their significance in Nepal's socio-economy, climate studies have seldom considered extremes, and even fewer have considered them in combination with temperature and precipitation. This study aimed at examining the trend of climate variables in Gandaki Province, Nepal. Daily temperature and precipitation data of five stations between 1990 and 2020 were analyzed. Modified Mann–Kendall and Sen's slope methods were used to detect trend and magnitude. The Mann–Whitney–Pettitt test was used to detect abrupt changes, and the Pearson correlation coefficient was used to find the correlation. The result showed an increasing trend and a significant abrupt change in the maximum temperature for all stations. A decreasing trend in the minimum temperature was observed in the Himalayas and the Hill region, whereas an increasing trend was seen in Siwalik and Terai regions. The Jomsom station, however, behaved differently by showing an increasing trend in precipitation and the number of rainy days. The majority of the temperature indices showed an increasing trend unlike precipitation indices, which showed a mixed result. The maximum five-day precipitation and consecutive dry days showed a significant positive correlation with altitude. The results indicate an increase in the frequency and intensity of extreme climate conditions in Gandaki Province.

  • The paper analyzes the trend in the mean and the extreme value of temperature and precipitation.

  • High altitude region could experience an extreme climate condition with frequent heavy rainfall and drought periods, while Siwalik and Terai regions could experience frequent drought periods.

  • The results will help the relevant stakeholders to understand the change happening in the Gandaki Province.

Climate change is a global issue that poses a significant risk to humanity. A long-term change in climate has been observed at a global scale caused by the increase in greenhouse gases (GHGs), particularly carbon dioxide (CO2) (Loo et al. 2015). This change has led to variations in precipitation patterns, temperature, wind patterns, and an increase in the frequency of extreme events (Loo et al. 2015). As per the Intergovernmental Panel on Climate Change (IPCC) report of 2022, the global average temperature rose by 1.1 °C between 2011 and 2020 compared with the average temperature in the late 19th century (Pörtner et al. 2022). The report highlights that this temperature rise is caused by the increase in GHGs (Pörtner et al. 2022). To understand the influence of climate on the frequency and intensity of extreme weather events, historical and future climate data have been widely studied (Felix et al. 2021). Variations in temperature and precipitation trends can easily alter the hydrological cycle, so these two variables play a crucial role in understanding the climate (Felix et al. 2021). The IPCC suggested that studies on historical and future climate trends are essential for developing strategies to cope with the impacts of climate change (Felix et al. 2021).

The University of East Anglia's Climatic Research Unit has found that there has been a significant rise in the global temperature following a decline from 1945 to 1973 (Loo et al. 2015). This increase is believed to be linked to anthropogenic causes such as carbon emissions and urbanization (Loo et al. 2015). This temperature change has altered the precipitation pattern in terms of its type, amount, intensity, and frequency (Azadani 2012). Eastern North America, southern South America, and northern Europe have seen an increase in precipitation, while the Mediterranean, most of Africa, and southern Asia have seen a decrease in precipitation (Azadani 2012). Additionally, the global climate model predicts that extreme precipitation events will become more intense during the 21st century worldwide (Li et al. 2019). Due to the influence of atmospheric circulation changes, extreme precipitation events are likely to intensify in the tropics and subtropics (Li et al. 2019).

Climate change has intensified many extreme events in South Asia as several South Asian countries have reported an increase in temperature and precipitation in recent years (Naveendrakumar et al. 2019). Naveendrakumar et al. (2019) has mentioned that India has experienced higher-than-usual temperature, with an average warming of 0.51 °C. Another study by Sheikh et al. (2015) found that there was a higher precipitation trend over the Himalayan region of South Asia.

The altitude of Nepal varies from 68 m in the south to 8,848 m in the north within a range of 193 km (Upadhayaya & Baral 2020). This variation in altitude within a short elevation has made climate change a very sensitive issue in Nepal (Upadhayaya & Baral 2020). As a result, Nepal is ranked the fourth most-vulnerable country to climate change (Joshi et al. 2019). Studies have shown that the mean annual temperature of Nepal has been rising by 0.06 °C/year, and the climate in the Himalayan region of Nepal is changing faster than the global average (Upadhayaya & Baral 2020). A study by Shrestha et al. (1999) on 49 stations in Nepal concluded that there was consistent and continuous warming after the mid-1970s. The study found that the warming was most pronounced in the higher altitude areas of Nepal and less significant in the Terai and Siwalik regions (Shrestha et al. 1999). Additionally, the warming was more pronounced in the winter season than in other seasons (Shrestha et al. 1999). However, the study did not find any significant trend in precipitation (Shrestha et al. 1999). A study on extreme climate in Nepal by Awasthi & Owen (2020) showed a significant positive and negative trend in the hot and cold extreme indices. The study showed a higher increase in the extreme temperature in the Himalayan and mountain regions. However, the study showed an insignificant increasing trend in the precipitation indices (Awasthi & Owen 2020).

Most research on climate change impact in Nepal has focused on the national level and these studies have often only focused on a single climate variable, i.e., temperature, precipitation, or extremes. While Nepal has seen an increase in extreme events, there have been few studies on climate extremes in the country. Also, few studies have been conducted with a primary focus on all three climate variables in Gandaki Province. Hence, this study was carried out to assess the climate change scenario in Gandaki Province by analyzing:

  • the trend in temperature, precipitation, and their extremes and

  • the correlation between the extreme climate indices.

The study will provide a better understanding of the trend of climate variables and their extremes in Gandaki Province. The public and private sectors can utilize this information to effectively formulate policies to adapt and mitigate the impacts of a changing climate.

Study area

The study area for this research is Gandaki Province (Figures 1 and 2). Gandaki Province is one of the seven federal provinces, which was established on 20 September 2015. The headquarters of the province are Pokhara. Geographically, the region is located between 27°20′N to 29°20′N and 82°52′E to 85°12′E. The highest point in the province is the Dhaulagiri Himal, which stands at 8,167 m, and the lowest point is the Triveni Susta of Nawalpur.

The total area of the province is 21,504 km2, which is approximately 14.57% of the entire country's area (Neupane et al. 2021). The Himalayan region accounts for 5,819 km2 (26.8%) of the total area, the hilly terrain covers 14,604 km2 (67.2%), and the Terai terrain makes up 1,310 km2 (6%) (Neupane et al. 2021).

The province stretches from the plains of Terai in the south to the cold Himalayan peaks in the north, resulting in a wide range of climatic conditions. This variation in climate, coupled with different agroecological regions, leads to a diverse range of vegetation, from subtropical/tropical to temperate/tundra. The Kaski district of Gandaki Province has the highest rainfall (2,711 mm) in Nepal, whereas the lowest rainfall is recorded in the Mustang district (258 mm) (Neupane et al. 2021).

Gandaki Province has a total population of 2,403,016, which is 9.06% of Nepal's total population (Neupane et al. 2021). The male population is 948,029, whereas the female population is 1,144,124 (Neupane et al. 2021). The region's urban population is 1,452,186 (60.5%), whereas the rural population is 943,652 (39.5%) (Neupane et al. 2021). The population of the Himalayan region is 19,990 (0.8%), of the hilly region is 2,072,162 (86.5%), and of the Terai region is 310,864 (12.7%). Around 14.19% of the total population of Gandaki Province is living in poverty (Neupane et al. 2021).

Data

Daily data for minimum temperature, maximum temperature, and precipitation for 30 years (1990–2020) were acquired from the Department of Hydrology and Meteorology (DHM), Government of Nepal. These data were collected for five stations within Gandaki Province, i.e., Chame, Jomsom, Lumle, Khudibazar, and Dumkauli of Manang, Mustang, Kaski, Lamjung, and Nawalparasi districts. The stations were selected in a way that represents all the climatic zones of the province as far as possible. Of the five stations, two stations are located in the Himalayas above 2,500 metres above sea level (masl), one station is located in the hilly region at an elevation of 1,738 masl, one station is located in the Siwalik region at an elevation of 838 masl, and one station is located in the Terai region at an elevation of 183 masl. The elevation of the station ranges from 183 masl (Dumkauli) to 2,741 masl (Jomsom). CMIP 6 data were acquired from the World Bank climate change knowledge portal, https://climateknowledgeportal.worldbank.org/.

All the stations have missing data (Table 1). On average, the stations have 5.4% missing data (0.3%–15.8%) for temperature and 1.9% (0.3%–4.5%) missing data for precipitation. The Chame station has the shortest temperature and precipitation record (1990–2012), whereas the other stations have data from 1990 to 2020.

Quality control and data fill

Ensuring the quality and completeness of the data was an essential step for this study. The temperature and precipitation data were inspected using RClimDex for outliers. RClimDex identifies outliers and erroneous data such as negative precipitation value and maximum temperature less than minimum temperature. We chose four standard deviations as a threshold for identifying outliers in RClimDex (Shrestha et al. 2017). The identified outliers were replaced by the average values from the previous day and the next day (DHM 2017).

Missing data for precipitation were filled with APHRODITE V1901 data (0.25° × 0.25°), which were obtained from https://www.chikyu.ac.jp/precip/. APHRODITE data are constructed from observed precipitation data along with additional pre-existing datasets (Yatagai et al. 2012). APHRODITE data were first bias-corrected using the linear scaling bias correction method (Shrestha 2015). Since APHRODITE data were not available after 2015, missing data after 2015 were filled by taking the average monthly values from previous years. Missing data for temperature were filled with CHELSA V2.1 data, which were obtained from https://envicloud.wsl.ch/. CHELSA is a global fine-scale climate dataset (Karger et al. 2017). CHELSA data were also first bias-corrected using the linear scaling bias correction method (Shrestha 2015).

Annual data and seasonal data for all the stations were calculated using the daily data. The seasons considered for this study are winter, i.e., December, January, and February (DJF); pre-monsoon, i.e., March, April, and May (MAM); Monsoon, i.e., June, July, and August (JJA); and post-monsoon, i.e., September, October, and November (SON).

Data analysis

Trend and step change analysis

Trend analysis is essential to determine and quantify the magnitude of the trend present in the climate time series data. Non-parametric tests, the Mann–Kendall test, and the Thiel–Sen slope test are widely used to detect the trend and quantify the magnitude of the trend, respectively (Felix et al. 2021). The Mann–Kendall test is used to assess trends and their significance in time series data. It does not require data to be normal, and it can tolerate missing data and outliers (Karki et al. 2017). The Thiel–Sen slope method is used to calculate the slope of the trend and is less sensitive to outliers and missing values (Karki et al. 2017). It has also been recommended by the World Meteorological Organization (WMO) to assess trends in environmental data time series (Rustum et al. 2017). In this study, we have used the modified Mann–Kendall test (Hamed & Rao 1998) to remove the serial correlation effect that is present in the time series data. To detect an abrupt change, the Mann–Whitney–Pettitt (MWP) test is used and to analyze the correlation, Pearson's correlation coefficient is used.
Figure 1

Map of Nepal showing Gandaki Province.

Figure 1

Map of Nepal showing Gandaki Province.

Close modal
Figure 2

Location of stations used in the study.

Figure 2

Location of stations used in the study.

Close modal
Figure 3

Annual temperature trend of Chame.

Figure 3

Annual temperature trend of Chame.

Close modal
Figure 4

Annual temperature trend of Khudibazar.

Figure 4

Annual temperature trend of Khudibazar.

Close modal
Figure 5

Annual temperature trend of Dumkauli.

Figure 5

Annual temperature trend of Dumkauli.

Close modal
Figure 6

Annual precipitation trend of Jomsom.

Figure 6

Annual precipitation trend of Jomsom.

Close modal
Figure 7

Maximum temperature trend of Chame by season.

Figure 7

Maximum temperature trend of Chame by season.

Close modal
Figure 8

Minimum temperature trend of Chame by season.

Figure 8

Minimum temperature trend of Chame by season.

Close modal
Figure 9

Number of FD trend.

Figure 9

Number of FD trend.

Close modal
Figure 10

Number of ID trend.

Figure 10

Number of ID trend.

Close modal
Figure 11

Minimum value of maximum temperature (TXn) trend.

Figure 11

Minimum value of maximum temperature (TXn) trend.

Close modal
Figure 12

DTR trend.

Figure 13

Maximum RX1 day trend.

Figure 13

Maximum RX1 day trend.

Close modal
Figure 14

CDD trend.

Figure 15

CWD trend.

Figure 16

Correlation between the trend magnitude of extreme climate indices.

Figure 16

Correlation between the trend magnitude of extreme climate indices.

Close modal

Extreme climate indices

The 17 extreme climate indices used in this study are summarized in Table 2.

RClimDex was used for the calculation of indices. It is a freely available R package with a graphical user interface and has been widely used (Shrestha et al. 2017). First, the upper and lower threshold values for maximum and minimum temperatures were defined. The upper threshold value for maximum and minimum temperatures was taken as the average temperature of MAM and DJF, respectively (Shrestha et al. 2017).

Table 1

Data availability and missing data

StationData availability (in year)Temperature missing data (%)Precipitation missing data (%)
Jomsom 1990–2020 5.5 3.5 
Chame 1990–2012 4.7 1.1 
Lumle 1990–2020 0.3 0.3 
Khudibazar 1990–2020 15.8 4.5 
Dumkauli 1990–2020 0.9 0.3 
StationData availability (in year)Temperature missing data (%)Precipitation missing data (%)
Jomsom 1990–2020 5.5 3.5 
Chame 1990–2012 4.7 1.1 
Lumle 1990–2020 0.3 0.3 
Khudibazar 1990–2020 15.8 4.5 
Dumkauli 1990–2020 0.9 0.3 
Table 2

Summary of 17 climate indices used in this study

Index classificationClimate indexIDDescriptionUnit
Cold duration Frost days FD Number of days when Tmin < 0 °C Days 
Icing days ID Number of days when Tmax < 0 °C Days 
Heat duration Summer days SU Number of days when Tmax > 25 °C Days 
Tropical nights TR Number of days when Tmin > 20 °C Days 
Intensity indices – temperature Coldest nighttime temperature TNn Monthly minimum value of daily minimum temperature °C 
Warmest nighttime temperature TNx Monthly maximum value of daily minimum temperature °C 
Coldest daytime temperature TXn Monthly minimum value of daily maximum temperature °C 
Warmest daytime temperature TXx Monthly maximum value of daily maximum temperature °C 
Diurnal temperature range DTR Monthly mean difference between Tmax and Tmin °C 
Intensity indices – precipitation Maximum 1-day precipitation RX1 day Maximum 1-day precipitation mm 
Maximum 5-day precipitation RX5 day Maximum consecutive 5-day precipitation mm 
Simple daily intensity index SDII Annual total precipitation divided by the number of wet days mm/day 
Total wet day precipitation PRCPTOT Total precipitation (PCP) in wet days (PCP ≥1 mm) mm 
Frequency indices – precipitation Heavy precipitation day R10 Number of days when precipitation ≥10 mm Days 
Very heavy precipitation day R20 Number of days when precipitation ≥20 mm Days 
Duration indices – precipitation Consecutive dry days CDD Maximum number of consecutive days with precipitation less than 1 mm Days 
Consecutive wet days CWD Maximum number of consecutive days with precipitation ≥1 mm Days 
Index classificationClimate indexIDDescriptionUnit
Cold duration Frost days FD Number of days when Tmin < 0 °C Days 
Icing days ID Number of days when Tmax < 0 °C Days 
Heat duration Summer days SU Number of days when Tmax > 25 °C Days 
Tropical nights TR Number of days when Tmin > 20 °C Days 
Intensity indices – temperature Coldest nighttime temperature TNn Monthly minimum value of daily minimum temperature °C 
Warmest nighttime temperature TNx Monthly maximum value of daily minimum temperature °C 
Coldest daytime temperature TXn Monthly minimum value of daily maximum temperature °C 
Warmest daytime temperature TXx Monthly maximum value of daily maximum temperature °C 
Diurnal temperature range DTR Monthly mean difference between Tmax and Tmin °C 
Intensity indices – precipitation Maximum 1-day precipitation RX1 day Maximum 1-day precipitation mm 
Maximum 5-day precipitation RX5 day Maximum consecutive 5-day precipitation mm 
Simple daily intensity index SDII Annual total precipitation divided by the number of wet days mm/day 
Total wet day precipitation PRCPTOT Total precipitation (PCP) in wet days (PCP ≥1 mm) mm 
Frequency indices – precipitation Heavy precipitation day R10 Number of days when precipitation ≥10 mm Days 
Very heavy precipitation day R20 Number of days when precipitation ≥20 mm Days 
Duration indices – precipitation Consecutive dry days CDD Maximum number of consecutive days with precipitation less than 1 mm Days 
Consecutive wet days CWD Maximum number of consecutive days with precipitation ≥1 mm Days 
Table 3

Summary of the annual temperature trend

StationWeather station data
CMIP6 model projection
Maximum temperature trend (°C/year)Minimum temperature trend (°C/year)Maximum temperature trend (°C/year)Minimum temperature trend (°C/year)
Jomsom 0.02 −0.01 0.03 0.03 
Chame 0.14 − 0.12 0.03 0.03 
Lumle 0.04 −0.003 0.03 0.03 
Khudibazar 0.03 0.01 0.03 0.03 
Dumkauli 0.03 0.02 0.03 0.03 
StationWeather station data
CMIP6 model projection
Maximum temperature trend (°C/year)Minimum temperature trend (°C/year)Maximum temperature trend (°C/year)Minimum temperature trend (°C/year)
Jomsom 0.02 −0.01 0.03 0.03 
Chame 0.14 − 0.12 0.03 0.03 
Lumle 0.04 −0.003 0.03 0.03 
Khudibazar 0.03 0.01 0.03 0.03 
Dumkauli 0.03 0.02 0.03 0.03 

Bold values are significant at α = 0.05.

Table 4

Summary of annual precipitation trend

StationPrecipitation trend (mm/year)
Weather station dataCMIP6 model projection
Jomsom 2.42 1.51 
Chame −1.8 1.69 
Lumle −15.68 1.51 
Khudibazar −7.16 1.7 
Dumkauli −5.38 0.6 
StationPrecipitation trend (mm/year)
Weather station dataCMIP6 model projection
Jomsom 2.42 1.51 
Chame −1.8 1.69 
Lumle −15.68 1.51 
Khudibazar −7.16 1.7 
Dumkauli −5.38 0.6 

Bold values are significant at α = 0.05.

Table 5

Summary of the number of rainy days trend

StationNumber of rainy days trend (days/year)
Jomsom 0.5 
Chame −0.43 
Lumle −0.25 
Khudibazar − 0.5 
Dumkauli −0.24 
StationNumber of rainy days trend (days/year)
Jomsom 0.5 
Chame −0.43 
Lumle −0.25 
Khudibazar − 0.5 
Dumkauli −0.24 

Bold values are significant at α = 0.05.

Table 6

Summary of the MWP test for change-point detection in annual trend

StationChange-point year
Maximum temperatureMinimum temperaturePrecipitation
Jomsom 2005 1994 2012 
Chame 2000 2004 2002 
Lumle 2000 2001 2004 
Khudibazar 2008 1997 2004 
Dumkauli 2014 1997 2007 
StationChange-point year
Maximum temperatureMinimum temperaturePrecipitation
Jomsom 2005 1994 2012 
Chame 2000 2004 2002 
Lumle 2000 2001 2004 
Khudibazar 2008 1997 2004 
Dumkauli 2014 1997 2007 

Bold values are significant at α = 0.05.

Table 7

Summary of the seasonal temperature trend

StationTemperatureDJF (°C/year)MAM (°C/year)JJA (°C/year)SON (°C/year)
Jomsom Max 0.022 0.021 0.003 0.048 
Min 0.002 0.016 −0.016 −0.034 
Chame Max 0.23 0.18 0.036 0.12 
Min −0.07 −0.024 − 0.31 −0.06 
Lumle Max 0.055 0.028 0.035 0.052 
Min −0.1 −0.013 0.013 −0.002 
Khudibazar Max 0.065 −0.002 −0.005 0.031 
Min 0.026 −0.004 −0.003 0.007 
Dumkauli Max 0.028 0.051 0.048 0.03 
Min 0.025 0.041 0.002 
StationTemperatureDJF (°C/year)MAM (°C/year)JJA (°C/year)SON (°C/year)
Jomsom Max 0.022 0.021 0.003 0.048 
Min 0.002 0.016 −0.016 −0.034 
Chame Max 0.23 0.18 0.036 0.12 
Min −0.07 −0.024 − 0.31 −0.06 
Lumle Max 0.055 0.028 0.035 0.052 
Min −0.1 −0.013 0.013 −0.002 
Khudibazar Max 0.065 −0.002 −0.005 0.031 
Min 0.026 −0.004 −0.003 0.007 
Dumkauli Max 0.028 0.051 0.048 0.03 
Min 0.025 0.041 0.002 

Bold values are significant at α = 0.05.

Table 8

Summary of the seasonal precipitation trend

StationDJF (mm/year)MAM (mm/year)JJA (mm/year)SON (mm/year)
Jomsom 0.67 0.69 1.5 −0.19 
Chame −2.96 −0.47 7.93 −5.32 
Lumle −0.71 0.12 −6.08 −5.35 
Khudibazar −0.81 0.14 −5.62 −0.53 
Dumkauli −0.2 4.53 −3.02 −2.92 
StationDJF (mm/year)MAM (mm/year)JJA (mm/year)SON (mm/year)
Jomsom 0.67 0.69 1.5 −0.19 
Chame −2.96 −0.47 7.93 −5.32 
Lumle −0.71 0.12 −6.08 −5.35 
Khudibazar −0.81 0.14 −5.62 −0.53 
Dumkauli −0.2 4.53 −3.02 −2.92 

Bold values are significant at α = 0.05.

Table 9

Summary of the trend in annual temperature indices

IndicesUnitJomsomChameLumleKhudibazarDumkauli
FD Days 0.63 
ID Days 
SU Days 0.24 0.66 1.02 0.15 
TR °C 0.06 
TNn °C 0.02 0.09 − 0.06 
TNx °C 0.03 − 0.43 
TXn °C 0.57 0.08 0.09 
TXx °C 0.06 −0.1 0.03 −0.01 
DTR °C 0.11 0.22 0.04 0.003 
IndicesUnitJomsomChameLumleKhudibazarDumkauli
FD Days 0.63 
ID Days 
SU Days 0.24 0.66 1.02 0.15 
TR °C 0.06 
TNn °C 0.02 0.09 − 0.06 
TNx °C 0.03 − 0.43 
TXn °C 0.57 0.08 0.09 
TXx °C 0.06 −0.1 0.03 −0.01 
DTR °C 0.11 0.22 0.04 0.003 

Bold values are significant at α = 0.05.

Table 10

Summary of the trend in annual precipitation indices

IndicesUnitJomsomChameLumleKhudibazarDumkauli
RX1 day mm 0.37 −0.71 −0.15 − 1.05 −0.91 
RX5 day mm 0.58 0.46 0.73 −1.79 −1.87 
SDII mm/day −0.05 0.08 0.02 0.12 
PRCPTOT mm 3.77 −3.41 −5 −6.91 −5.36 
R10 Days 0.27 −0.17 −0.33 0.26 
R20 Days −0.29 −0.19 −0.06 0.09 
CDD Days 0.82 0.9 0.43 0.19 −0.47 
CWD Days 0.08 0.24 0.47 −0.18 
IndicesUnitJomsomChameLumleKhudibazarDumkauli
RX1 day mm 0.37 −0.71 −0.15 − 1.05 −0.91 
RX5 day mm 0.58 0.46 0.73 −1.79 −1.87 
SDII mm/day −0.05 0.08 0.02 0.12 
PRCPTOT mm 3.77 −3.41 −5 −6.91 −5.36 
R10 Days 0.27 −0.17 −0.33 0.26 
R20 Days −0.29 −0.19 −0.06 0.09 
CDD Days 0.82 0.9 0.43 0.19 −0.47 
CWD Days 0.08 0.24 0.47 −0.18 

Bold values are significant at α = 0.05.

Table 11

Summary of the correlation between extreme climate indices and altitude

IndicesAltitude
FD 0.55 
ID −0.4 
SU −0.37 
TR −0.72 
TNn 0.84 
TNx −0.48 
TXn 0.49 
TXx −0.08 
DTR 0.83 
RX1 day 0.68 
RX5 day 0.89 
SDII −0.5 
PRCPTOT 0.71 
R10 0.13 
R20 −0.6 
CDD 0.97 
CWD 0.48 
IndicesAltitude
FD 0.55 
ID −0.4 
SU −0.37 
TR −0.72 
TNn 0.84 
TNx −0.48 
TXn 0.49 
TXx −0.08 
DTR 0.83 
RX1 day 0.68 
RX5 day 0.89 
SDII −0.5 
PRCPTOT 0.71 
R10 0.13 
R20 −0.6 
CDD 0.97 
CWD 0.48 

Bold values are significant at α = 0.05.

Annual trend

Temperature

Based on the result of the annual trend in temperature, we can see a significant trend in the maximum and minimum temperatures at Chame, Khudibazar, and Dumkauli (Table 3).

The maximum temperature at Chame showed a significant increasing trend at the rate of 0.14 °C/year (Table 3, Figure 3). This is the highest positive trend seen among the five stations and indicates a strong warming trend at high altitudes. Similar results were published by Khand et al. (2020), DHM (2017), Shrestha et al. (1999), and Bajracharya et al. (2011), where a higher warming trend was observed in the high-altitude region. A similar trend can also be observed in the whole Hindu Kush Himalayan region (Shrestha & Aryal 2011).

Snow albedo feedback influences the temperature at higher altitudes (Pepin et al. 2022). The reduction of snow cover will alter the albedo of the region, causing an increase in the temperature (Shrestha et al. 1999). A study by Gurung et al. (2017) showed a decreasing trend in the snow cover in the Gandaki basin. This reduction of snow cover may be contributing to the increasing temperature trend in the high-altitude region (Shrestha et al. 1999). Similarly, the maximum temperature at Khudibazar showed a significant increasing trend at the rate of 0.03 °C/year (Table 3, Figure 4).

The minimum temperature at Chame showed a significant decreasing trend at the rate of 0.12 °C/year (Table 3, Figure 3), which is the highest negative trend seen among the five stations selected in this study. Similar results were published by DHM (2017) and Bajracharya et al. (2011). However, the minimum temperature at Dumkauli showed a significant increasing trend at the rate of 0.02 °C/year (Table 3, Figure 5).

A decreasing trend in the minimum temperature can be observed in the Himalayan and Hill regions, whereas an increasing trend can be seen in the Siwalik and Terai regions (Table 3). Similar findings were published by DHM (2017).

Comparing our findings with the multi-model ensemble CMIP6 climate model projection for the year 2015–2100, the model result showed a significant increasing trend in the annual maximum and minimum temperatures (Table 3). Under the SSP2-4.5 scenario, the annual maximum and minimum temperatures over the selected stations are projected to increase at the rate of 0.03 °C/year. In the case of maximum temperature, both our findings and CMIP6 projection results showed a similar trend; however, in the case of minimum temperature, the CMIP6 projection showed a different trend in the Himalayan and Hill regions compared with our findings.

Precipitation

The precipitation at Jomsom showed a significant increasing trend at a rate of 2.42 mm/year (Table 4, Figure 6). This result is consistent with the findings published by Bhadra et al. (2021). A similar increasing trend in precipitation was also observed in the Tibetan Plateau (Zhan et al. 2017). A study by Sigdel et al. (2022) over the Gandaki basin also showed a significant increasing trend in precipitation over the Himalayan region.

Fort (2014) has mentioned increasing temperature trends at higher altitudes as a cause for increasing precipitation. The increasing temperature trend could influence permafrost distribution causing less snowfall and more rainfall.

The climate in the Trans Himalayan region is believed to be influenced by subtropical westerlies and the Indian summer monsoon (Zhu et al. 2015). So, the increase in precipitation could be because of the weather system developed by the collision of monsoons and westerlies supported by jet streams (The Kathmandu Post 2021).

Jomsom and Chame both lie in the Trans-Himalayan region, but the precipitation trend is different for these two stations (Table 4). In an area of complex topography, especially with drastic elevation changes, precipitation is locally determined (You et al. 2015). Factors such as mountain shape and relief and land cover affect precipitation in such regions (You et al. 2015). Also, these deep valleys establish their own wind system causing different precipitation patterns (Chen et al. 2018). Further research needs to be carried out to find out which factors play a key role in the distribution of precipitation in areas of complex topography.

The remaining stations showed a decreasing trend but were not statistically significant. The highest decreasing trend is seen in the Lumle station of Kaski district at the rate of 15.68 mm/year, which is similar to the findings published by DHM (2017).

Comparing our findings with the multi-model ensemble CMIP6 climate model projection for the years 2015–2100, the model result showed a significant increasing trend in the annual precipitation at Jomsom, Chame, Lumle, and Khudibazar and an insignificant increasing trend at Dumkauli (Table 4). Under the SSP2-4.5 scenario, the highest increasing trend is seen at the Khudibazar station. Our findings and the CMIP6 projection showed an increasing trend at the Jomsom station; however, the result was different for other stations.

The analysis of the number of rainy days, i.e., precipitation >2.5 mm in a day, showed a significant trend at Jomsom and Khudibazar (Table 5). The number of rainy days is significantly increasing and decreasing by 0.5 days/year at Jomsom and Khudibazar, respectively. All the remaining stations showed decreasing trends but these were not statistically significant.

The increase in the maximum and minimum temperature, decrease in precipitation, and the number of rainy days at Khudibazar and Dumkauli indicate that heatwaves and drought events could become more common in the Siwalik and Terai regions.

Change-point detection

The result from the MWP test to detect an abrupt change in the temperature and precipitation is presented in Table 6. There was no consistent year detected at the stations; however, all five stations showed a significant change-point year for maximum temperature. Due to the lack of metadata, further research needs to be carried out to understand the significance of the detected change-point years.

Seasonal trend

Temperature

The seasonal maximum and minimum temperatures showed a significant increasing trend at most of the stations (Table 7). The change in the maximum and minimum temperatures at Chame is generally higher compared with the other stations. During the pre-monsoon and post-monsoon seasons, Chame showed a significantly high increasing trend in the maximum temperature in comparison with the other stations. Similar results were published by Marahatta et al. (2009). The highest significant increase in the maximum temperature is seen during the winter season at Chame (Table 7, Figure 7), which is similar to the findings published by Konchar et al. (2015). Shrestha & Aryal (2011) also found that warming in the winter season in high-altitude regions is more than in the other seasons.

Chame showed a decreasing trend in the minimum temperature with a significant decreasing trend during the monsoon season (Figure 8). Marahatta et al. (2009) also found a high decreasing trend in the minimum temperature at Chame during the monsoon season.

The increase in the seasonal maximum temperature was generally higher than that of the seasonal minimum temperature, which is similar to the findings published by Khand et al. (2020) and Mishra et al. (2014).

Precipitation

The seasonal precipitation showed a significant increasing trend at Jomsom during the monsoon season and Dumkauli during the pre-monsoon season at the rate of 1.5 and 4.53 mm/year, respectively (Table 8).

Except for Chame, all the stations showed an increasing trend in the pre-monsoon season (Table 8). This trend is similar to the findings published by Karki et al. (2017) and Sigdel et al. (2022). This rise in precipitation in the pre-monsoon season indicates an increase in intense thunderstorms over the region (Karki et al. 2017).

There is an increasing trend in precipitation seen only in the Himalayan region during the monsoon season (Table 8). Similar findings were published by Karki et al. (2017) where a significant increasing trend was observed in the Himalayan region of central Nepal during the monsoon season.

The post-monsoon and winter seasons mostly showed a decreasing trend over the region (Table 8). Similar findings were published by Panthi et al. (2015) where a decreasing trend in post-monsoon and winter seasons was observed over the Gandaki River basin from 1981 to 2012. Winter crops may suffer due to a lack of soil moisture as a result of decreasing trends of precipitation in the post-monsoon and winter seasons (Panthi et al. 2015).

Trend of extreme climate indices

Temperature indices

Out of nine temperature indices used, five indices showed a significant trend (Table 9).

Jomsom showed an increasing trend (statistically insignificant) in frost days (FDs) at the rate of 0.63 days/year (Table 9, Figure 9). However, a study by Manandhar et al. (2012) over the Kali Gandaki River basin for the period of 1978–2007 found a significantly decreasing trend in the FD at Jomsom. The difference in the trend could be due to the data taken from the different timeframes. Icing days (IDs) were not recorded at any station (Figure 10), which is similar to the findings published by Manandhar et al. (2012).

The number of summer days (SUs) showed a significant increasing trend at Jomsom, Lumle, and Khudibazar stations at the rate of 0.24, 0.66, and 1.02 days/year, respectively. The minimum value of daily minimum temperature (TNn) and the maximum value of daily minimum temperature (TNx) showed a significantly decreasing trend at Dumkauli and Chame. The TNn and TNx values decreased at the rate of 0.06 and 0.43 °C/year at Dumkauli and Chame, respectively.

The minimum value of daily maximum temperature (TXn) showed a significant increasing trend at Chame, Lumle, and Khudibazar at the rate of 0.57, 0.08, and 0.09 °C/year, respectively (Table 9, Figure 11). This trend is similar to the findings published by Shrestha et al. (2017). This warming trend in Himalayan, Hill, and Siwalik regions could be due to GHG emission and changes in cloud cover (Shrestha et al. 2017).

The mean difference between maximum and minimum temperatures (diurnal temperature range, DTR) showed a significant increasing trend at Chame, Jomsom, and Lumle at the rate of 0.22, 0.11, and 0.04 °C/year, respectively (Table 9, Figure 12). This trend is similar to the findings published by Shrestha et al. (2017), where ten of the 12 hill/mountain stations of the Koshi River basin showed a significant increasing trend in the DTR. The positive trend in the DTR indicates that the maximum temperature is increasing faster than the minimum temperature.

Precipitation indices

Out of eight precipitation indices used, three indices showed a significant trend (Table 10).

The maximum one-day precipitation (RX1 day) showed a decreasing trend at all the stations except Jomsom (Table 10, Figure 13). However, a significant trend was seen only at the Khudibazar station where the RX1 day value decreased at the rate of 1.05 mm/year. This trend is similar to the findings obtained by Karki et al. (2017) where a decreasing trend was observed in the Hill and Siwalik regions of central Nepal from 1970 to 2012.

The consecutive dry days showed an increasing trend (statistically insignificant) at all stations except Dumkauli (Table 10, Figure 14). Similar findings were obtained by Shrestha et al. (2017) where 47 of the 50 stations showed an increasing trend in the consecutive dry days (CDD). The increase in CDD could lead to an increase in the frequency of drought.

The total precipitation in wet days (PRCPTOT) showed a significant increasing trend at Jomsom at the rate of 3.77 mm/year, and the remaining stations showed a decreasing trend (statistically insignificant) (Table 10).

The consecutive wet days (CWD) showed an increasing trend at Jomsom, Chame, and Lumle, a decreasing trend at Khudibazar, and no trend at Dumkauli (Table 10, Figure 15). Jomsom showed a significant increasing trend at the rate of 0.08 days/year. Similar results were obtained by Shrestha et al. (2017), where 33 stations showed an increasing trend and 17 stations showed a decreasing trend.

The finding indicates that extreme climate conditions with heavier rainfall and increasing drought period could become more common at higher altitude regions. It also indicates that the trend in temperature extremes is more significant than precipitation extremes, which is similar to the findings obtained by Shrestha et al. (2017).

Correlation between the trend magnitude of extreme climate indices

The correlation between the annual trend magnitude of extreme climate indices is shown in Figure 16.

Based on the result, the DTR showed a high positive correlation (statistically significant) with TNn. TNn showed a high positive correlation (statistically significant) with CDD. TNx showed a high negative and positive correlation (both statistically significant) with TXn and TXx, respectively. TXn showed a high negative correlation (statistically significant) with TXx. FD showed a high positive correlation (statistically significant) with PRCPTOT. SU showed a high negative correlation (statistically significant) with R10.

Correlation between extreme climate indices and altitude

The correlation between the trend magnitude of extreme climate indices with altitude is shown in Table 11.

RX5 day and CDD showed a significant positive correlation with altitude. This finding provides further evidence for the frequent occurrence of extreme climate conditions at the higher altitude region.

The study focused on the assessment of climate change in Gandaki Province, Nepal, by analyzing the trend in the average and extreme values of maximum temperature, minimum temperature, and precipitation as well as analyzing the correlation between the extreme climate indices from 1990 to 2020.

The result showed an increasing trend in the maximum temperature at all stations, where Chame and Khudibazar showed a significant increasing trend. For minimum temperature, the Himalayan and Hill regions showed a decreasing trend with a significant decreasing trend observed at Chame, and Siwalik and Terai regions showed an increasing trend with a significant increasing trend observed at Dumkauli. The seasonal maximum and minimum temperatures showed a significant increasing trend in most of the stations except for Chame, which showed a significant decreasing trend in minimum temperature in the monsoon season. The increase in the seasonal maximum temperature was found to be higher than the seasonal minimum temperature. Also, a significant change-point year was detected for the maximum temperature at all stations.

The result showed a significant increasing trend in precipitation at Jomsom and a decreasing trend in the remaining stations. The seasonal precipitation showed an increasing trend in the Himalayan region during the monsoon season. The number of rainy days showed an increasing trend at Jomsom and a decreasing trend in the remaining stations.

The majority of the temperature indices showed an increasing trend. The number of summer days showed an increasing trend, where significant changes were found at Jomsom, Lumle, and Khudibazar. The minimum value of maximum temperature also showed an increasing trend, where significant changes were found at Chame, Lumle, and Khudibazar. Similarly, the diurnal temperature range also showed an increasing trend, where significant changes were found at Jomsom, Chame, and Lumle.

The precipitation indices showed mixed results, which suggest complex processes in precipitation extremes. The maximum one-day precipitation showed a decreasing trend at all stations except Jomsom. However, maximum five-day precipitation showed an increasing trend in the Himalayan and Hill regions and a decreasing trend in Siwalik and Terai regions. The number of very heavy rainfall days showed a decreasing trend in Himalayan, Hill, and Siwalik regions and an increasing trend in the Terai region. The CDD and CWD showed an increasing trend in the Himalayan and Hill regions. RX5 day and CDD showed a significant positive correlation with altitude

The analysis presented here indicates that frequent and intense extreme climate conditions will increase in the future, causing significant loss and damage to lives and properties in the province. Hence, further local-level studies should be carried out to better assess and understand the changes in climate variables. Mitigation plans and strategies should be developed and implemented to tackle the changing climate.

All relevant data are included in the paper or its Supplementary Information.

The authors declare there is no conflict.

Awasthi
R. P.
&
Owen
J. S.
2020
Observed climate extreme in Nepal
.
The Geographic Base
7
,
1
14
.
https://doi.org/10.3126/TGB.V7I0.34262
.
Azadani
F. N.
2012
Modeling the Impact of Climate Change on Water Resources. Case Study: Arkansas River Basin in Colorado
.
Master's thesis
,
Colorado State University
,
Fort Collins, CO, USA
. .
Bajracharya
T. R.
,
Acharya
S.
&
Ale
B. B.
2011
Changing climatic parameters and its possible impacts in hydropower generation in Nepal (a case study on Gandaki River Basin)
.
Journal of the Institute of Engineering
8
(
1–2
),
160
173
.
https://doi.org/10.3126/JIE.V8I1-2.5108
.
Bhadra
R.
,
Neupane
B.
&
Khadka
U. R.
2021
Induced impacts of climate change on livelihood and migration in upper Himalayas: a case of Mustang, Nepal
. In:
Handbook of Climate Change Management
(Filho, W. L., Luetz, J. & Ayal, D., eds)
,
Springer, Cham
,
Switzerland
.
https://doi.org/10.1007/978-3-030-22759-3_303-2
.
Chen
R.
,
Han
C.
,
Liu
J.
,
Yang
Y.
,
Liu
Z.
,
Wang
L.
&
Kang
E.
2018
Maximum precipitation altitude on the northern flank of the Qilian Mountains, northwest China
.
Hydrology Research
49
(
5
),
1696
1710
.
https://doi.org/10.2166/NH.2018.121
.
Department of Hydrology and Meteorology (DHM)
2017
Observed Climate Trend Analysis of Nepal (1971–2014)
.
Department of Hydrology and Meteorology, Government of Nepal
,
Kathmandu, Nepal
. .
Felix
M. L.
,
Kim
Y. K.
,
Choi
M.
,
Kim
J. C.
,
Do
X. K.
,
Nguyen
T. H.
&
Jung
K.
2021
Detailed trend analysis of extreme climate indices in the upper Geum River Basin
.
Water
13
(
22
),
3171
.
https://doi.org/10.3390/W13223171
.
Gurung
D. R.
,
Maharjan
S. B.
,
Shrestha
A. B.
,
Shrestha
M. S.
,
Bajracharya
S. R.
&
Murthy
M. S. R.
2017
Climate and topographic controls on snow cover dynamics in the Hindu Kush Himalaya
.
International Journal of Climatology
37
(
10
),
3873
3882
.
https://doi.org/10.1002/JOC.4961
.
Hamed
K. H.
&
Rao
A. R.
1998
A modified Mann-Kendall trend test for autocorrelated data
.
Journal of Hydrology
204
(
1–4
),
182
196
.
https://doi.org/10.1016/S0022-1694(97)00125-X
.
Joshi
A.
,
Farquhar
S.
,
Assareh
N.
,
Dahlet
L.
&
Landahl
E.
2019
Climate change in Lamjung District, Nepal: meteorological evidence, community perceptions, and responses
.
Environmental Research Communications
1
(
3
),
031004
.
https://doi.org/10.1088/2515-7620/ab1762
.
Karger
D. N.
,
Conrad
O.
,
Böhner
J.
,
Kawohl
T.
,
Kreft
H.
,
Soria-Auza
R. W.
,
Zimmermann
N. E.
,
Linder
H. P.
&
Kessler
M.
2017
Climatologies at high resolution for the earth's land surface areas
.
Scientific Data
4
(
1
),
170122
.
https://doi.org/10.1038/sdata.2017.122
.
Karki
R.
,
ul Hasson
S.
,
Schickhoff
U.
,
Scholten
T.
&
Böhner
J.
2017
Rising precipitation extremes across Nepal
.
Climate
5
(
1
),
4
.
https://doi.org/10.3390/CLI5010004
.
Khand
P. B.
,
KC
V. K.
,
Dhungana
A. R.
&
Dhungana
S. M.
2020
Assessing the trend of climate change in the western Nepal
.
American Journal of Environmental Sciences
16
(
6
),
112
121
.
https://doi.org/10.3844/AJESSP.2020.112.121
.
Konchar
K. M.
,
Staver
B.
,
Salick
J.
,
Chapagain
A.
,
Joshi
L.
,
Karki
S.
,
Lo
S.
,
Paudel
A.
,
Subedi
P.
&
Ghimire
S. K.
2015
Adapting in the shadow of Annapurna: a climate tipping point
.
Journal of Ethnobiology
35
(
3
),
449
471
.
https://doi.org/10.2993/0278-0771-35.3.449
.
Li
C.
,
Zwiers
F.
,
Zhang
X.
,
Chen
G.
,
Lu
J.
,
Li
G.
,
Norris
J.
,
Tan
Y.
,
Sun
Y.
&
Liu
M.
2019
Larger increases in more extreme local precipitation events as climate warms
.
Geophysical Research Letters
46
(
12
),
6885
6891
.
https://doi.org/10.1029/2019GL082908
.
Loo
Y. Y.
,
Billa
L.
&
Singh
A.
2015
Effect of climate change on seasonal monsoon in Asia and its impact on the variability of monsoon rainfall in Southeast Asia
.
Geoscience Frontiers
6
(
6
),
817
823
.
https://doi.org/10.1016/J.GSF.2014.02.009
.
Manandhar
S.
,
Pandey
V. P.
&
Kazama
F.
2012
Hydro-climatic trends and people's perceptions: case of Kali Gandaki River Basin, Nepal
.
Climate Research
54
(
2
),
167
179
.
https://doi.org/10.3354/CR01108
.
Marahatta
S.
,
Dongol
B. S.
&
Gurung
G. B.
2009
Temporal and Spatial Variability of Climate Change over Nepal (1976–2005)
.
Practical Action
,
Kathmandu, Nepal
. .
Mishra
B.
,
Babel
M. S.
&
Tripathi
N. K.
2014
Analysis of climatic variability and snow cover in the Kaligandaki River Basin, Himalaya, Nepal
.
Theoretical and Applied Climatology
116
(
3–4
),
681
694
.
https://doi.org/10.1007/S00704-013-0966-1
.
Naveendrakumar
G.
,
Vithanage
M.
,
Kwon
H. H.
,
Chandrasekara
S. S. K.
,
Iqbal
M. C. M.
,
Pathmarajah
S.
,
Fernando
W. C. D. K.
&
Obeysekera
J.
2019
South Asian perspective on temperature and rainfall extremes: a review
.
Atmospheric Research
225
,
110
120
.
https://doi.org/10.1016/J.ATMOSRES.2019.03.021
.
Neupane
S.
,
Khatri
L.
,
Bhusal
A.
,
Neupane
G.
&
Shivakoti
S.
2021
Climate Investment Plan for the Agriculture Sector: A Decision Support Tool for Scaling Up Climate-Smart Agriculture Technologies and Practices in Gandaki Province, Nepal
.
Local Initiatives for Biodiversity, Research and Development (LI-BIRD), Pokhara and Ministry of Land Management, Agriculture, Cooperative and Poverty Alleviation (MoLMACPA), Gandaki Province, Nepal. Available from: https://hdl.handle.net/10568/116349 (accessed 13 June 2022)
.
Panthi
J.
,
Dahal
P.
,
Shrestha
M. L.
,
Aryal
S.
,
Krakauer
N. Y.
,
Pradhanang
S. M.
,
Lakhankar
T.
,
Jha
A. K.
,
Sharma
M.
&
Karki
R.
2015
Spatial and temporal variability of rainfall in the Gandaki River Basin of Nepal Himalaya
.
Climate
3
(
1
),
210
226
.
https://doi.org/10.3390/CLI3010210
.
Pepin
N. C.
,
Arnone
E.
,
Gobiet
A.
,
Haslinger
K.
,
Kotlarski
S.
,
Notarnicola
C.
,
Palazzi
E.
,
Seibert
P.
,
Serafin
S.
,
Schöner
W.
,
Terzago
S.
,
Thornton
J. M.
,
Vuille
M.
&
Adler
C.
2022
Climate changes and their elevational patterns in the mountains of the world
.
Reviews of Geophysics
60
(
1
),
e2020RG000730
.
https://doi.org/10.1029/2020RG000730
.
Pörtner
H.-O.
,
Roberts
D. C.
,
Poloczanska
E. S.
,
Mintenbeck
K.
,
Tignor
M.
,
Alegría
A.
,
Craig
M.
,
Langsdorf
S.
,
Löschke
S.
,
Möller
V.
&
Okem
A.
2022
Summary for policymakers
. In:
Climate Change 2022: Impacts, Adaptation and Vulnerability. Contribution of Working Group II to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change
(
Pörtner
H.-O.
,
Roberts
D. C.
,
Tignor
M.
,
Poloczanska
E. S.
,
Mintenbeck
K.
,
Alegría
A.
,
Craig
M.
,
Langsdorf
S.
,
Löschke
S.
,
Möller
V.
,
Okem
A.
&
Rama
B.
, eds).
Cambridge University Press
,
Cambridge
,
UK
and New York, NY, USA
, pp.
3
33
.
Rustum
R.
,
Adeloye
A. J.
&
Mwale
F.
2017
Spatial and temporal trend analysis of long term rainfall records in data-poor catchments with missing data, a case study of Lower Shire floodplain in Malawi for the period 1953–2010
.
Hydrology and Earth System Sciences Discussions
[preprint]. https://doi.org/10.5194/HESS-2017-601
.
Sheikh
M. M.
,
Manzoor
N.
,
Ashraf
J.
,
Adnan
M.
,
Collins
D.
,
Hameed
S.
,
Manton
M. J.
,
Ahmed
A. U.
,
Baidya
S. K.
,
Borgaonkar
H. P.
,
Islam
N.
,
Jayasinghearachchi
D.
,
Kothawale
D. R.
,
Premalal
K. H. M. S.
,
Revadekar
J. V.
&
Shrestha
M. L.
2015
Trends in extreme daily rainfall and temperature indices over South Asia
.
International Journal of Climatology
35
(
7
),
1625
1637
.
https://doi.org/10.1002/JOC.4081
.
Shrestha
M.
2015
Linear Scaling Bias Correction (Version 1.0) [Microsoft Excel file]. https://dx.doi.org/10.13140/RG.2.1.3365.8967
.
Shrestha
A. B.
&
Aryal
R.
2011
Climate change in Nepal and its impact on Himalayan glaciers
.
Regional Environmental Change
11
(
Suppl. 1
),
65
77
.
https://doi.org/10.1007/S10113-010-0174-9
.
Shrestha
A. B.
,
Wake
C. P.
,
Mayewski
P. A.
&
Dibb
J. E.
1999
Maximum temperature trends in the Himalaya and its vicinity: an analysis based on temperature records from Nepal for the period 1971–94
.
Journal of Climate
12
(
9
),
2775
2786
.
https://doi.org/10.1175/1520-0442(1999)012<2775:MTTITH>2.0.CO;2
.
Shrestha
A. B.
,
Bajracharya
S. R.
,
Sharma
A. R.
,
Duo
C.
&
Kulkarni
A.
2017
Observed trends and changes in daily temperature and precipitation extremes over the Koshi river basin 1975–2010
.
International Journal of Climatology
37
(
2
),
1066
1083
.
https://doi.org/10.1002/JOC.4761
.
The Kathmandu Post
2021
Mountain districts are receiving more rainfall this monsoon. The Kathmandu Post (19 June). Available from: https://kathmandupost.com/climate-environment/2021/06/19/mountain-districts-are-receiving-more-rainfall-this-monsoon (accessed: 29 November 2022)
.
Upadhayaya
R. P.
&
Baral
M. P.
2020
Trends of climate change in some selected districts of western Nepal
.
Janapriya Journal of Interdisciplinary Studies
9
(
1
),
148
158
.
https://doi.org/10.3126/JJIS.V9I1.35284
.
Yatagai
A.
,
Kamiguchi
K.
,
Arakawa
O.
,
Hamada
A.
,
Yasutomi
N.
&
Kitoh
A.
2012
APHRODITE: constructing a long-term daily gridded precipitation dataset for Asia based on a dense network of rain gauges
.
Bulletin of the American Meteorological Society
93
(
9
),
1401
1415
.
https://doi.org/10.1175/BAMS-D-11-00122.1
.
You
Q.
,
Min
J.
,
Zhang
W.
,
Pepin
N.
&
Kang
S.
2015
Comparison of multiple datasets with gridded precipitation observations over the Tibetan Plateau
.
Climate Dynamics
45
(
3–4
),
791
806
.
https://doi.org/10.1007/S00382-014-2310-6
.
Zhan
Y. J.
,
Ren
G. Y.
,
Shrestha
A. B.
,
Rajbhandari
R.
,
Ren
Y. Y.
,
Sanjay
J.
,
Xu
Y.
,
Sun
X. B.
,
You
Q. L.
&
Wang
S.
2017
Changes in extreme precipitation events over the Hindu Kush Himalayan region during 1961–2012
.
Advances in Climate Change Research
8
(
3
),
166
175
.
https://doi.org/10.1016/J.ACCRE.2017.08.002
.
Zhu
L.
,
X.
,
Wang
J.
,
Peng
P.
,
Kasper
T.
,
Daut
G.
,
Haberzettl
T.
,
Frenzel
P.
,
Li
Q.
,
Yang
R.
,
Schwalb
A.
&
Mäusbacher
R.
2015
Climate change on the Tibetan Plateau in response to shifting atmospheric circulation since the LGM
.
Scientific Reports
5
(
1
),
13318
.
http://dx.doi.org/10.1038/srep13318
.
This is an Open Access article distributed under the terms of the Creative Commons Attribution Licence (CC BY-NC-ND 4.0), which permits copying and redistribution for non-commercial purposes with no derivatives, provided the original work is properly cited (http://creativecommons.org/licenses/by-nc-nd/4.0/).