Trends in the Diurnal Temperature Range Over the Southern Slope of Central Himalaya: Retrospective and Prospective Evaluation

: The Diurnal Temperature Range (DTR) profoundly affects human health, agriculture, eco ‐ system, and socioeconomic systems. In this study, we analyzed past and future changes in DTR using gridded Climate Research Unit (CRU) datasets for the years 1950–2020 and an ensemble means of thirteen bias ‐ corrected Coupled Model Intercomparison Project Phase 6 (CMIP6) models under different Shared Socioeconomic Pathways (SSP1 ‐ 2.6, SSP2 ‐ 4.5, and SSP5 ‐ 8.5) scenarios for the rest of the 21st century over the southern slope of Central Himalaya, Nepal. Furthermore, the po ‐ tential drivers (precipitation and cloud cover) of seasonal and annual DTR were studied using cor ‐ relation analysis. This study found that the DTR trends generally declined; the highest decrease was observed in the pre ‐ monsoon and winter at a rate of 0.09 °C/decade ( p ≤ 0.01). As expected, DTR demonstrated a significant negative correlation with cloudiness and precipitation in all four sea ‐ sons. Further, the decreased DTR was weakly related to the Sea Surface Temperature variation (SST) in the tropical Pacific and Indian Oceans. We found that the projected DTR changes in the future varied from a marginal increase under the SSP1 ‐ 2.6 (only pre ‐ monsoon) scenario to continued sig ‐ nificant decreases under SSP2 ‐ 4.5 and SSP5 ‐ 8.5. Insights based on retrospective and prospective evaluation help to understand the long ‐ term evolution of diurnal temperature variations. Trends in the Diurnal Temperature Range Over Southern Slope of and


Introduction
The surface air temperature is an essential climatic variable mediating the diversity of biological, physical, and chemical processes guaranteeing the sustainability of life on earth [1,2]. Due to ongoing climate change, the global land surface temperature increased by 0.99 °C from 1850-1900 to 2001-2020; recent decades have been successively warmer than previous decades [3]. The observed variation in the mean temperature (Tmean) trend corresponds to either the maximum temperature (Tmax) or the minimum temperature (Tmin) trend; otherwise, it is relative to both temperature trends. The complicated climate and leading to extreme weather conditions [36]. Different phases of the Coupled Model Intercomparison Project (CMIP) models have been developed to understand future climatic changes arising from natural, unforced variability or in response to changes in radiative forcing in a multi-model context. For instance, 20 different Coupled Model Inter-Comparison Project phase 5 (CMIP5) models have been used to project DTR globally, revealing a reduction in the globally averaged DTR in most of these models [37]. Further, a decreasing DTR trend is observed in CMIP phase 6 (CMIP6) models; however, CMIP6 models can simulate DTR better than CMIP5 [38]. Furthermore, CMIP6 models can unfold future climatic conditions due to improved emissions scenarios, land use scenarios, and computational advances [39,40]. Moreover, the better-simulating capacity of CMIP6 for temperatures over the South Asian region was already verified by [41]. However, none of these studies has estimated the future DTR trends using CMIP6 models; therefore, this study aims to analyze the projected (2021-2100) annual and seasonal DTR changes using 13 bias-corrected CMIP6 models over Nepal.
Overall, this study attempts to answer three research questions: (a) how did annual and seasonal DTR trends vary over the last 71 years (1950-2020)? (b) What are the possible factors responsible for the annual and seasonal DTR change? And (c) how will climate models project future DTR trends (2021-2100) over Nepal?

Study Area
Nepal (26°22′ to 30°27′ N 80°04′ to 88°12′ E), a south Asian country located between China and India, is a geographically diverse country (Figure 1a). The altitude varies from below 100 m to over 8000 m above sea level within a short latitudinal distance. It extends from ~885 km east to west and 130 km to 260 km from north to south, covering 147,516 km 2 . The country's climate varies from subtropical to alpine [42]. The seasons in Nepal are categorized as pre-monsoon (March-May), monsoon (June-August), post-monsoon (September-November), and winter (December-February) [43]. Pre-monsoon is hot and dry, with prevailing westerly windy weather and generally localized precipitation. Monsoon is characterized by humid, southeasterly winds blowing from the Bay of Bengal and by widespread precipitation [42,44]. During the post-monsoon season, rainfall is significantly reduced, with November usually the driest month [45][46][47]. Winter (December-February) is generally dry, and some rains occur due to westerly disturbances, primarily in the western highlands [17]. (c) anomalies of annual Tmean in CRU (trend; 0.12 °C/decade, p < 0.01) and CMIP6 (trend; 0.11 °C/decade, p < 0.01) historical data over the study region.

Gridded Data
Daily gridded measurements of precipitation, Tmin, Tmax, Tmean, and cloud cover data over Nepal from the Climate Research Unit CRU TS v.4.05 (released on 17 March 2021) were used in this study [48]. CRU gridded datasets are derived by the interpolation of monthly climate anomalies from extensive weather station observations around the globe and provide climatological data at a 0.5° resolution over all land domains except for Antarctica. These datasets are updated regularly from several sources, i.e., CLIMAT messages exchanged internationally between World Meteorological Organization (WMO) countries, obtained as quality-controlled files via the UK Met Office; Monthly Climatic Data for the World (MCDW) summaries, obtained from the US National Oceanographic and Atmospheric Administration (NOAA), via its National Climate Data Centre (NCDC); and updates on the minimum and maximum temperatures for Australia, obtained from the Bureau of Meteorology (BoM). In addition, ad hoc collections of stations are incorporated (after quality control checks, including location, correspondence to existing holdings, and outlier checking). The CRU dataset is spatially and temporally complete over land, but where station data are unavailable or very limited, the interpolated gridded values are relaxed towards the mean for the baseline period 1961-1990 [48,49]. This is the case for DTR before 1950 in Nepal. Thus, we used CRU data only from 1950 to 2020 for the current study. More details on the datasets can be found in [48] and at https://crudata.uea.ac.uk/cru/data/hrg/ (accessed on 1 September 2021).

Climatic Indices
This study used the monthly Indo-Pacific climatic indices of NINO3.4 and DMI between 1950 and 2020 developed by the National Centers for Environmental Information (NOAA). We used the NINO3.4 and DMI from the preceding winter and autumn, as they appeared to be stronger. These datasets are freely available on the NOAA's Physical Sciences Laboratory: https://psl.noaa.gov/gcos_wgsp/Timeseries/ (accessed on 1 September 2021).

Model Simulations
We used 13 bias-corrected CMIP6 models (ACCESS-CM2, ACCESS-ESM1-5, BCC-CSM2-MR, CanESM5, EC-Earth3, EC-Earth3-Veg, INM-CM4-8, INM-CM5-INM-CM5-0, MPI-ESM1-2-HR, MPI-ESM1-2-LR, MRI-ESM2-0, NorESM2-LM, NorESM2-MM) for Nepal, with a 0.25° × 0.25° spatial resolution. These are the bias-corrected models developed for the South Asian countries by Mishra, Bhatia, and Tiwari [39]. Indian Meteorological Department (IMD) gridded datasets were used for the Indian region, and data generated by Sheffield, et al. [50]  The spatial distribution of the mean monthly temperatures over Nepal between 1950 and 2020 is presented in Figure 1a. Temperature variation over the country is directly associated with altitude, local winds, and season [29,52,53]. The Tmean varies considerably from south to north; the maximum (minimum) temperature is observed in the flat southern strip (high mountains or the Himalayas in the north) ( Figure 1a). The seasonal cycle of the CRU and CMIP6 historical for temperatures across Nepal is presented in Figure 1b. CMIP6 slightly underestimated the monthly variation of the CRU temperature; however, both CRU and CMIP6 historical mean annual temperature anomalies displayed very similar magnitudes (~0.12 °C/decade, p < 0.01) of increasing trends during 1951-2014, indicating that CMIP6 can represent the temperature variation and trend over Nepal (Figure 1c).

Methodology
Daily DTR data were obtained by subtracting the daily minimum temperature from the daily maximum temperature. Further, the monthly, seasonal, and annual DTR were calculated by averaging daily DTR for each month. For this purpose, the daily DTR was averaged for each month; subsequently, the average DTR for each month was arranged according to season (pre-monsoon, monsoon, post-monsoon, and winter). For the historical (1951-2014) and projected (2021-2100) CMIP6 datasets, 13 bias-corrected models were averaged to generate the Tmax and Tmin daily time series. These CMIP6 data were also validated across Nepal using CRU datasets (Figure 1b,c). Next, spatially averaged annual and seasonal mean DTRs were obtained through a yearly and seasonal average of daily DTR. Next, the temporal averages were conducted for annual and seasonal (pre-monsoon, monsoon, post-monsoon, and winter) analysis to obtain the DTR time series for Nepal.
The Mann-Kendall (MK) statistical test [54,55] was used to assess the monotonic upward or downward trend of the climatic variables for seasonal and annual timescales. The test was proposed by Mann [54], further studied by Kendall [55], and approved by Hirsch and Slack [56]. The MK statistical test compares the relative magnitudes of sample data rather than the values themselves and uses them to identify the monotonic trends present in time series data. It tests either to reject the null hypothesis (Ho) or to accept the alternative hypothesis (Ha). The Ho of the test is that there is no trend in the data, and the Ha represents a monotonic increasing or decreasing trend [54]. Further, the test statistic (S) is computed using Equations (1) and (2), and the variance is calculated using Equation (3). Many studies have used the MK test to detect a trend [4,5,[57][58][59]. (1) where xi and xj are the values of sequence i, j; n is the number of observations for a particular climate data. If n is higher than 8, statistic S approximates to normal distribution. The mean of S is 0, and the variance of S can be acquired as follows: Next, the test statistic Z is generated using Equation (4).
Positive and negative Z values in Equation (4) indicate increasing and decreasing trends, respectively. For a two-tailed test, the null hypothesis (Ho) is rejected at a given significance level α for an absolute value of |Z| ≥ Z1 − α/2. In this study, the α is taken at p < 0.1 (90%), p < 0.05 (95%), and p < 0.01 (99%). Furthermore, the magnitude of a time series trend was estimated by Sen's non-parametric estimator (Equation (4)); the trend is calculated by Equation (5); where is Sen's slope [60], x denotes the climatic variable, and I and j are indices. The values > 0 and < 0 indicate upward and downward trends in a time series, respectively.
Additionally, a Pearson correlation coefficient (R 2 ; Equation (6)) was performed using the spatially averaged time series of cloud cover, precipitation, and sea surface temperature for their possible effects on or association with DTR.
where R 2 = Pearson correlation coefficient and Cov(x,y) represents the covariance of the climatic variables x and y. The values σx and σy denote the standard deviation of x and y, respectively.

Seasonal and Annual DTR Trends over Nepal
The MK statistical test and Sen's Slope estimator were performed using the CRU datasets to present the annual and seasonal DTR trends and their magnitudes across Nepal from 1950-2020. An increasing annual and seasonal Tmax and Tmin trends were observed across Nepal. The Tmin increment rate was faster, ~2 times higher than the Tmax ( Figure  2 and Table 1). Meanwhile, [29] also mentioned the warming trend in all seasons, except the cooling in the winter season, especially over the southern lowland of Nepal. Decreasing seasonal DTR trends were observed during the study period, with a significant decrease (Figures 2a,d,g,j) in the pre-monsoon (−0.09 °C/decade) and winter (−0.09 °C/decade) compared to the monsoon (−0.07 °C/decade), and post-monsoon season (−0.06 °C/decade) ( Table 1). It was found that, in the pre-monsoon and winter season, the Tmax increased at a rate of 0.06 and 0.11 °C/decade, respectively, which is smaller than Tmin (0.15 and 0.20 °C/decade, respectively) (Table 1). Similarly, a significant decreasing DTR trend at a rate of −0.08 °C/decade was observed in the annual timescale ( Figure 2m, Table  1). Furthermore, a recent study noted the decreasing trend of annual DTR at a rate of 0.02 °C/decade during 1991-2016 in India [57]. The decreasing trend and higher DTR variability were observed mainly after the 1950s. This variation is similar to the global DTR trend (−0.036 °C/decade) from 1901-2014, mainly due to the significant decrease in DTR from 1951 to 2014 [6]. Overall, the results support the previous conclusion that the decreasing DTR is mainly due to the Tmin increasing faster than the Tmax [4,61]. A low DTR can exert adverse health impacts on human beings [10]; it has been observed that the increased occurrence of cardiovascular diseases is more pronounced in females and the elderly than males and adults [62].   (4) and (5)).

Factors Affecting DTR
The MK and Sen's slope estimator was used to reveal the annual and seasonal cloud cover and precipitation trends and their magnitudes, presented in Table 2. It was found that cloud cover increased across Nepal, with a significant increasing trend in the premonsoon (0.56%/decade) and winter (0.42%/decade) at 99% and 95% confidence levels, respectively. However, an insignificant increasing trend was observed for the monsoon and post-monsoon seasons. The annual cloud cover displayed an increasing rate of 0.36%/decade (p < 0.05) during the study period. Similarly, increased cloud cover was observed over the Indo-Gangetic plains and northeast India during 1961-2010 [63]. The increase in the cloud cover over the country represents a cooling effect during those months [29]. Moreover, the precipitation trend over Nepal exhibited a decreasing trend, except during the pre-monsoon season. The pre-monsoon precipitation increased at a 3.90 mm/decade rate, significant at a 99% confidence level. By contrast, a larger decreasing rate was observed in the monsoon season (−4.36 mm/decade, p < 0.01) than in the post-monsoon (−0.79 mm/decade, p = 0.50) and winter (−0.42 mm/decade, p = 0.44). Moreover, the annual precipitation was observed to decrease by −0.48 mm/decade and p = 0.52. Similarly, a study that used 143 stations over Nepal also demonstrated decreasing annual precipitation, except for the western high-hills from 2001-2016 [45]. Table 2. Trends of seasonal and annual total cloud cover (%/decade) and precipitation (mm/decade) in Nepal from 1950-2020. Trend values are calculated using Sen's slope estimator, and the statistical significance of the trends is assessed by MK test (Equations (4) and (5)). Previous studies have already highlighted that the considerable DTR variation is due to the total cloud cover and precipitation [61,64,65]; therefore, their correlation (Equation (6)) with the annual and seasonal DTR is presented in Figure 3. The DTR correlation with the cloud cover and precipitation was performed at different confidence levels. The weaker correlation in the annual timescale (−0.13, p = 0.28) indicates that precipitation may not influence DTR on an interannual timescale; however, this variation may change on the decadal and multidecadal timescales [64]. Conversely, a significant negative correlation was observed on a seasonal timescale for winter (−0.42, p < 0.01), pre-monsoon (−0.54, p < 0.01), monsoon (0.25, p < 0.05), and post-monsoon (−0.50, p < 0.01). The winter, premonsoon, and post-monsoon seasons are comparatively drier than the monsoon; therefore, the result indicates low DTR during high precipitation and vice-versa. Further, it was suggested that the precipitation deficit led to an increase in DTR extremes [66].

Seasons Cloud Precipitation
We further analyzed the possible relationship between DTR and the cloud cover over Nepal. A significant negative correlation was observed for all the seasons and annual time scales, with the cloud cover at a 99% confidence level. The correlation was <−0.55 for the monsoon, post-monsoon, winter, and annual, except for the pre-monsoon (−0.41); this indicates that cloud cover can substantially affect DTR variation over the country ( Figure  3). Similarly, the inverse relationship between DTR and precipitation and cloud cover is reported in northern Pakistan [67] and Thailand [68]. Further, the cloud cover demonstrated increasing trends during all the seasons ( Table 2), indicating that increased cloud cover results in a decrease in Tmax and an increase in Tmin. Moreover, air pollution and aerosols are highest during the winter in Nepal [69]; the increase of pollution in the air combined with temperature inversion creates smog and fog (https://airlief.com/air-pollution-during-winter/ (accessed on 1 September 2021)). The increasing cloud cover/fog over the country may have reduced the DTR by dampening the diurnal cycle of the radiation balance, as suggested by [70,71].
Additionally, to analyze the cloud and precipitation's relationship with the DTR trend, we developed a linear regression model between the interannual variations of DTR and the climatic variables (DTR = a + b*cloud + c*precip). Comparing the regression-based trends in Appendix Table A1 with the observed DTR trends in Table 1, the regression captures the sign of the DTR trends but not their full magnitude (particularly in the monsoon and post-monsoon seasons). The results further revealed that the significant decrease in DTR during pre-monsoon can be linked to the significantly increased precipitation over Nepal, as high precipitation dampens the Tmax through evaporative cooling [65]. Further, an increasing trend of cloud cover during pre-monsoon and winter matches well with the long-term decreasing trend of DTR. The increased cloudiness and precipitation reduced the DTR, mainly for the pre-monsoon and winter seasons, but additional factors such as attenuating solar radiation by increasing aerosols might also have played a role during other seasons.

Relation with Climatic Indices
The seasonal and annual DTR over Nepal is correlated with the NINO3.4 in winter and DMI in autumn, as presented in Table 3. The summer DTR corresponds positively with the ENSO but with a weak relationship (0.27, p < 0.05). A previous study has shown that the warmer surface air temperature of India is associated with the positive sea surface temperature (SST) anomalies over the equatorial eastern-central Pacific, representing El Niño effects [72]. La Nina Phase can be associated with decreased DTR, which induces stronger low-level winds and greater moisture availability, and reduces incoming radiation; the opposite is the case for warmer temperatures during the El Nino phase [73]. On the other hand, weak-to-moderate influences of DMI were observed in the spring, autumn, winter, and annual phases, in which the ENSO influence was not observed. The results indicate that the positive DMI phase (warming in the western Indian Ocean > cooling in the eastern Indian Ocean) can reduce DTR. Similarly, the temperature during spring, autumn, and winter can be influenced by the Indian Ocean SST [73]. Overall, the decrease in DTRs is weakly related to the Sea Surface Temperature (SST) variation in the tropical Pacific and the Indian Ocean.

Future Projection of DTR
A significant decreasing trend in annual DTR was observed in both CMIP6 historical (−0.09 °C/decade, p < 0.01) and CRU data sets (−0.07 °C/decade, p < 0.01). The results indicate that CMIP6 can simulate the temporal trend of DTR across Nepal (Appendix - Figure  A1). The ensemble mean of the 13 bias-corrected CMIP6 datasets was used to project (2021-2100) the annual and seasonal trends of DTR, Tmax, and Tmin under three different scenarios. Under the SSP1-2.6, a significantly increasing trend in projected Tmax compared with Tmin for pre-monsoon leads to an increase in DTR at a rate of 0.02 °C/decade (p < 0.05) (Figures 4a-c and Table 4). The projected Tmax and Tmin feature the same increasing trends, presenting no change in DTR trend during the monsoon, post-monsoon, and winter seasons (Figures 4d-l). Annually, DTR is likely to increase at a rate of 0.01 °C/decade (p < 0.05) due to a slightly higher increasing trend in Tmax than in Tmin ( Figure  4m-o, Table 4). A previous study projected an annual and seasonal temperature increase across Nepal for all three scenarios, in which the mean temperature is likely to increase by 1.3-4.5 °C [41]. The projected DTR, Tmax, and Tmin under the SSP2-4.5 scenario are presented in Figure 5. The result highlights that Tmin (monsoon, post-monsoon, winter, and annual) will feature a higher increasing trend than Tmax except for pre-monsoon. On a seasonal scale, the DTR is projected to decrease significantly at a rate of −0.08 °C, −0.05 °C, −0.02 °C per decade for monsoon, post-monsoon, and winter, respectively (Figures 5d-l, Table 4). However, Tmax and Tmin will feature a similar warming trend at a rate of 0.24 °C/decade, presenting no significant change in DTR during pre-monsoon (Figures 5a-c). In addition, a recent study projected an increase in summer temperature over south Asia in all scenarios [41]. The results further indicate that a higher decreasing trend in DTR is likely to be observed in the monsoon season compared to other seasons. Annually, the DTR will also exhibit a significantly decreasing trend, at a rate of −0.04 °C/decade until 2100 (Figures  5m-o, Table 4). These findings are generally consistent with [37], according to which the globally averaged reduction in DTR is projected to increase due to greenhouse gases. In the SSP5-8.5 scenario, Tmin will increase significantly more than Tmax, exhibiting a significantly decreasing DTR seasonal and annual trend ( Figure 6). Both Tmin and Tmax are projected to increase sharply in high-emission scenarios; this may lead to harsh weather conditions in the future. The increase in temperature was estimated over Central Nepal; moreover, Tmin and Tmax will be more pronounced in higher altitudinal stations than lower altitude stations [74]. Further, a recent study also estimated worsening extreme heat conditions in the south and east Asia by the end of this century [75]. The decreasing rate of DTR trend will be higher in monsoon (−0.21 °C/decade) than post-monsoon (−0.17 °C/decade), winter (−0.11 °C/decade), and pre-monsoon (−0.03 °C/decade) (Figures 6a-l ,  Table 4). Similarly, the DTR will decrease significantly on an annual timescale at a rate of −0.14 °C/decade, which is significant at a 99% confidence level (Figures 6m-o, Table 4). Overall, the decreasing rate of DTR is most likely to be higher in high (SSP5-8.5) than medium (SSP2-4.5) and low (SSP1-2.6) greenhouse emission scenarios. Further, unabated human-induced global warming is expected to enhance the decreasing rate of DTR. Alarmingly, the decreasing DTR trend is already causing an increase in mortality rates [10,76]; in addition to the DTR-related deaths, which are projected to increase by 1.4-10.3% in 2090-99, depending upon the variability, this pattern may vary between countries and regions [1].  Table 4. Trends of projected (2021-2100) seasonal and annual average diurnal temperature range (DTR), maximum temperature (Tmax), and minimum temperature (Tmin) (unit: o C/decade) under SSP1-2.6, SSP2-4.5, and SSP5-8.5 over Nepal. Trend values are calculated using Sen's slope estimator, and the statistical significance of the trends is assessed by MK test (Equations (4) and (5)).

Conclusions and Policy Suggestions
This study investigated historical (1950-2020) and projected (2021-2100) DTR trends on annual and seasonal timescales over the southern slope of Central Himalaya, Nepal, using CRU and 13 bias-corrected CMIP6 models. An increasing trend in Tmax and Tmin on seasonal and annual timescales was observed from 1950-2020 in Nepal; the Tmin's increasing trend was ~2 times higher than that of the Tmax. This variational trend further decreased the DTR (0.08 °C/decade, p ≤ 0.01) on an annual timescale. Similarly, the magnitude of seasonal DTR has declined over Nepal, with a more significant decrease in premonsoon (−0.09 °C/decade, p ≤ 0.01) and winter (−0.09 °C/decade, p ≤ 0.01) compared to the monsoon (−0.07 °C/decade, p ≤ 0.01) and post-monsoon seasons (−0.06 °C/decade, p ≤ 0.05). A stronger relationship between DTR and precipitation was observed for dry seasons (winter, pre-monsoon, and post-monsoon) than for the wet season (monsoon). Meanwhile, a significant negative correlation was observed between DTR and cloud cover for all seasons, indicating that high cloud cover is responsible for decreasing DTR across Nepal. Further, weak-to-moderate correlations were observed between the DTR and the sea surface temperature variations in the tropical Pacific and Indian Oceans.
The CMIP6 historical and CRU exhibited a similar magnitude of temperature and DTR, indicating that the CMIP6 datasets can simulate temperature trends across Nepal. The ensemble of climate models projects an increase in Tmax and Tmin on seasonal and annual timescales in SSP1-2.6, SSP2-4.5, and SSP5-8.5 scenarios. In the SSP1-2.6 scenario, the DTR trend is projected to increase only in pre-monsoon. A decreasing DTR trend is anticipated in SSP2-4.5, with a higher increasing rate in the monsoon season followed by post-monsoon, winter, and no trend in the pre-monsoon. A similar DTR variation is projected in SSP5-8.5, with pre-monsoon featuring a smaller DTR decreasing trend (−0.03 °C/decade, p < 0.01). An increasing annual DTR trend is expected to be observed in the SSP1-2.6 scenario. By contrast, the annual DTR might decrease in the SSP2-4.5 and SSP5-8.5 scenarios. Our study suggests that the DTR is likely to decrease faster in SSP5-8.5 than in SSP2-4.5 and SSP1-2.6, indicating that high greenhouse emissions could lead to a higher increase in Tmin than in Tmax. In the context of ongoing climate change, this study could thus be helpful for the preparation and implementation of preparedness and adaptation strategies to counter imminent DTR-related threats. Moreover, this study only analyzed the temporal variation of DTR over Nepal; we recommend that further studies focus on the spatial pattern of DTR changes over the Central Himalayan region using fine-resolution satellite/reanalysis datasets.  Data Availability Statement: The CRU gridded precipitation, temperature, and cloud cover datasets used in this study can be freely accessed from https://crudata.uea.ac.uk/cru/data/hrg/ (accessed on 1 September 2021). The gridded bias-corrected CMIP6 datasets for Nepal were generated by Mishra et al. (2020) and are freely available at https://zenodo.org/record/3873998#.YKOjiqgzY2w (accessed on 1 February 2021).