Climate change and spatio-temporal trend analysis of climate extremes in the homogeneous climatic zones of Pakistan during 1962-2019

Climate extremes, such as heat waves, droughts, extreme rainfall can lead to harvest failures, flooding and consequently threaten the food security worldwide. Improving our understanding about climate extremes can mitigate the worst impacts of climate change and extremes. The objective here is to investigate the changes in climate and climate extremes by considering two time slices (i.e., 1962–1990 and 1991–2019) in all climate zones of Pakistan by utilizing observed data from 54 meteorological stations. Different statistical methods and techniques were applied on observed station data to assess changes in temperature, precipitation and spatio-temporal trends of climatic extremes over Pakistan from 1962 to 2019. The Mann-Kendal test demonstrated increasing precipitation (DJF) and decreasing maximum and minimum temperatures (JJA) at the meteorological stations located in the Karakoram region during 1962–1990. The decadal analysis, on the other hand, showed a decrease in precipitation during 1991–2019 and an increase in temperature (maximum and minimum) during 2010–2019, which is consistent with the recently observed slight mass loss of glaciers related to the Karakoram Anomaly. These changes are highly significant at 5% level of significance at most of the stations. In case of temperature extremes, summer days (SU25) increased except in zone 4, TX10p (cold days) decreased across the country during 1962–1990, except for zones 1 and 2. TX90p (warm days) increased between 1991–2019, with the exception of zone 5, and decreased during 1962–1990, with the exception of zones 2 and 5. The spatio-temporal trend of consecutive dry days (CDD) indicated a rising tendency from 1991 to 2019, with the exception of zone 4, which showed a decreasing trend. PRCPTOT (annual total wet-day precipitation), R10 (number of heavy precipitation days), R20 (number of very heavy precipitation days), and R25mm (very heavy precipitation days) increased (decreased) considerably in the North Pakistan during 1962–1990 (1991–2019). The findings of this study can help to address some of the sustainable development goals related climate action, hunger and environment. In addition, the findings can help in developing sustainable adaptation and mitigation strategies against climate change and extremes. As the climate and extremes conditions are not the uniform in all climate zone, therefore, it is suggested to the formers and agriculture department to harvest crops resilient to the climatic condition of each zone. Temperature has increasing trend in the northern Pakistan, therefore, the concerned stakeholders need to make rational plans for higher river flow/flood situation due to snow and glacier melt.


Introduction
Climate change has been attributed to human-caused activities, especially increased greenhouse gases emissions after the mid-twentieth century [1][2][3][4]. Climate change is causing a number of threats, including an increase in the frequency of extreme events and other disasters [5,6]. Evidence of observed changes in extremes such as heavy precipitation, heatwaves, tropical cyclones, and droughts, and, in particular, the attribution of these extremes to human influence, has strengthened since AR5 (Fifth Assessment Report) [6]. The societal infrastructure is becoming more vulnerable to climate extremes, which are exacerbated by climate change [7][8][9]. Climate extremes such as droughts, floods, heatwaves, and typhoons have geographically different economic, health, social, ecological, and environmental impacts [10][11][12][13][14][15].
Many approaches have been developed to explore variability and trends of climate extremes using observed precipitation and temperature data around the world. For example, Grima et al. [16] investigated the variability of temperature and precipitation in the Upper Huai River Basin of China using yearly time series data from 1960 to 2016, concluded a decreasing trend in precipitation and a considerably increasing trend in temperature at all meteorological stations. Bougara et al. [17] investigated the variability of rainfall in the Tafna basin in northeastern Algeria using observed data from 1979 to 2011, resulted in heterogeneous results in different stations. Asfaw et al. [18] performed variability and trend analyses of temperature using data from 1901 to 2014 in the Woleka sub-basin in northcentral Ethiopia and found a significant increasing trend in mean and minimum average temperatures. Other studies that are based on observe station data on variability and trend analysis of temperature and precipitation in various regions of the globe include: [19][20][21][22][23][24].
Most of the studies conducted in Pakistan considered a sub-region of Pakistan for their analyses [25][26][27][28][29]. For instance, Alam et al. [25] considered three meteorological stations in the north-western part (Khyber Pakhtoonkhwa province) and used Mann Kendall (MK) test for trend analysis using daily observed data for the duration of 1960-2020. Safdar et al. [27] investigated the climate change indicators and spatio-temporal shift in Monsoon pattern by using observed data for the duration of 1980-2017. They concluded a warming tendency and a decrease in the rainfall in the Monsoon region of Pakistan. Hussain et al. [26] investigated variability and trends in precipitation and temperature over the northern part of Pakistan by using observed meteorological data (duration of 1955-2016 for some station and 1995-2012 for other stations). They concluded elevation-dependent conclusions in the target location, for instance, they noted increasing trend in maximum temperature and precipitation at station with elevation 1500-2800 meters in contrast to elevation less than 1500 meters during 1955-2016. Further studies which investigated variability and trends analysis in the sub-regions of the country include [28][29][30].
Studies related to precipitation and temperatures' extremes in different parts of the world revealed that these extremes are increased in terms of frequency, intensity and persistence [31][32][33]. The study of Sun et al. [33] concluded that the extreme precipitation has increased at two-zone and for each station in the HKH region. Moreover, the results of this study may help to address some of the Sustainable Development Goals (SDGs) of the United Nation Development Programme (UNDP) related to climate action, water, hunger and environment. The remaining paper is structured as: section 2 is for study area and data; section 3 presents the methodology; section 4 consists of results; section 5 is reserved for discussion and section 6 presents conclusions about this study.

Study area and data
The target location of this study is Pakistan which has 881,913 square kilometers area and approximately 222 million inhabitants [49][50][51]. Pakistan lies in the northern hemisphere with latitude of 23.4-37.3N and longitude of 60.5-78.5E. It has a diversity in climate as it has Arabian Sea in the south and the world's highest mountains (Himalaya, Hindukush and Karakoram) in the north. The country has four climate seasons (Winter = December-February (DJF), Spring = March-May (MAM), Summer = June-August (JJA), Autumn = September-November (SON)) and consequently significant spatial and temporal variation in the climate. In Fig 1, zone 1 comprises cold regions include northern and central-western parts of the country. Zone 2 comprised of Islamabad Capital Territory (ICT), some parts of Azad-Jammu and Kashmir and Khyber Pakhtoonkhwa which is part of the monsoon dominated region and receive maximum rainfall during the monsoon season. Zone 3 is the south-western part of the country and has least rainfall in a year as compared to the other regions of the country. Zone 4 includes the south-eastern and coastal areas in the south of the country. Zone 5 is the central-eastern part of the country and due to its fertile land, it is considering the food basket of the country. Observed daily climate data (maximum, minimum temperature and precipitation) for the period of 1962-2019 was acquired from Pakistan Meteorological Department (PMD) for fiftyfour (54) meteorological stations across the country. The spatial distribution of meteorological stations is shown in Fig 1. The data was divided into two independent climatic periods (1962-1990 and 1991-2019) for the assessment of climate change and extremes analysis over Pakistan.

Inconsistencies in the data
There were some missing values in the maximum, minimum temperature and precipitation particularly at the stations situated at high elevation. There are various statistical methods available to address the issue of missing values, for instance, multivariate regression, interpolation, imputation etc. In this study, imputation method (a regression-based approach) is used at multiple stations to preserve the spatial dependence structure. Multiple chains have been run and the convergence is assessed in each chain after a specified number of iterations. The number of iterations depend on researcher and the nature of the study [52]. However, the convergence of the process can be assessed at each iteration. The process can be stopped if there are no significant changes during the successive iterations. It is important to note that the inhomogeneous data need to be homogenized before the start of further analysis [53]. Inhomogeneities in data may occur due to breaks, changes or disorder in measuring instruments, observing practices of enumerators, location, and environment of the stations. The standard normal homogeneity tests (SNHT) are useful tools to assess the homogeneity, consistency and to detect the spikes and outliers in daily time series [54][55][56][57]. The commonly used SNHT for testing the homogeneity of data include Buishand Range test, Pettitt test, and Von Neumann Ratio test. In this study we have used Buishand Range test for testing the homogeneity of the data which can identify the month/year where the gap occurs. For further details about the above-mentioned homogeneity tests, we refer to [58].

Climate change assessment
For climate change assessment, the average differences between the two time slices are tested for possible statistical significance at a 5% level of significance. This is carried out under the hypothesis testing procedures of the student t-test. Under this procedures, the null hypothesis states that there is no difference between the average climatology of two-time slices where the alternative hypothesis states opposite to the null hypothesis. The hypothesis under the student t-test can be formulated as: Alternatively, it can be written as: Where H 0 and H A represent null and alternative hypothesis, respectively. The notation μ  and μ 1991-2019 represent the average climatology during 1962-1990 and 1991-2019, respectively. Here the average climatology means average temperature or precipitation for the mentioned time duration. The change between the two time slices is significant at 5% level of significance if the p-value of the test is less than or equal to 0.05. To assess the detailed behaviors of zone-wise climatology, this analysis was carried out on a seasonal as well as on annual basis.

Trend analysis
For trend analysis, a non-parametric test, the modified version of MK (mMK) test has been used initially proposed by [59] and widely uses in hydrology and climatology for trend analysis (e.g., [60,61]). The null hypothesis in MK test is "no monotonic trend is present in the data". Which means that if the p-value of the test is less than or equal to 0.05, then it is concluded that there is a trend in the data at 5% level of significance. For further details about the test statistic and calculation of the MK test, we refer to [62].

Climate extremes analysis
Climate extremes related to temperature and precipitation were calculated for both time slices. To assess the changes in climate extremes, probability density functions were drawn for both time slices and a comparison is made between them. To investigate spatial variability and changes, the changes in climate extremes were interpolated to the target location by using a geostatistical interpolation method called Kriging. The analysis provides information about the changes in climate extremes during the two-time slices over the full domain. Finally, the statistical significance of the climate extremes was performed by using a 5% level of significance. A comparison was made between climate extremes in both time slices for temperature and precipitation in each climate zone in terms of statistical significance. For definition and unit of measurement of climate extremes, we refer to [63,64]. A summary of the climate extremes used in this study is given in Annex-01.

Climate change assessment and extremes' analysis in HKH region
HKH comprises one of the world's largest non-polar glaciers [65][66][67], therefore, climate change and extremes can have significant impacts on snow and glacier melt which can further alter the hydrological regime of this region. Different studies conducted about trends in temperature, precipitation and extremes over Pakistan recently. For instance, the results of Bhatti et al. [68] suggested that the linear trend of all precipitation extreme indices showed an increasing tendency over Pakistan during 1980-2016. The study of Ullah et al. [69] illustrated that the spatial distributions of Tmin and Tmax showed a warming trend over the whole country, however, the seasonal and annual Tmin (Tmax) exhibited sharp increasing trends in the southern, southwestern and southeastern (northern and southwestern mountainous) regions of Pakistan. Elevation dependent significant fluctuations have been noted in long-term precipitation trend (during 1980-2016) in Pakistan with these categories: decreases in 500-1000m and 1000-1500m elevation zones; significant increasing trend in elevation zones of �500m and�1500 m [70]. Strong variability has been noted in precipitation during 1986-2005 over Pakistan, particularly in the northern parts of the country. For instance, a positive trend in the precipitation variability was noted during the winter season in Himalaya and Karakorum region while a negative trend was noted during pre-monsoon season in the Karakorum and Himalayas regions [71]. Studies related to extremes analysis suggested that the precipitation and temperature extremes like warm nights, summer days, warm days, warmest day, and coldest day increased during 1980-2016 over Pakistan [72,73]. Archer [74] investigated that with the rise of 1˚C in average temperature can cause an increase of 16-17% in river flow in the upper Indus Basin.
Few previous studies compared temperature and precipitation trends on local scale glaciers in the Karakoram and Himalayas in the UIB [75,76]. These findings suggest that slight temperature and precipitation supported the marginal retreating trends. The glaciers changes in correlation with the local climate improves the assessment of contemporary and future glaciers behavior under the changing climate. Our long-term climate assessment would improve the cryosphere changes in the glacier's rich mountainous region in the context of the Karakoram Anomaly.
There are seven meteorological stations with long-term observations in this region as shown in Table 1. The long-term meteorological changes were investigated by using observed data at each station. The significance of difference (average changes) between the two periods (1962-1990 and 1991-2019) was tested by utilizing a t-test to see whether it is significantly changing or not for each station separately on a seasonal basis. This can help to see the significance of the changes in average temperature and precipitation on a seasonal basis in different parts of the Karakoram region. In addition, the mMK test was used to identify a monotonic trend in temperature and precipitation. Finally, climate extremes were calculated for each station to see whether it is behaving differently than the zone wise climate extremes. This special part may disclose the latest status of climate change in different parts of the Karakoram and may update the present status of Karakoram anomaly which was first investigated by [43]. The details of meteorological station in Northern Pakistan are given in Table 1, with four stations are in the Karakoram and the remaining in the Himalaya and Hindukush region. The step-bystep methodology of this study is presented in Fig 2.

Results
Before the start of formal analysis, the data was investigated carefully for missing values and inhomogeneity. After filling the missing values using imputation, the results of SNHT test show that the data is homogeneous in time. Therefore, there was no need of any adjustment. The climate change extreme results of this study are discussed zone-wise and thereafter the new findings about the Karakoram anomaly. Considering the importance of the Karakoram region, climate extremes analysis is also presented for each station separately in this region.

Estimation
Zone 1. The difference between the two-time slices (1991-2019 and 1962-1990) was tested for possible significance in temperature and precipitation on a seasonal basis. A significant increase in temperature during DJF and MAM and a significant decrease during JJA and

PLOS ONE
Climate change and extremes analysis using observed data (1962-2019) over Pakistan SON was observed as shown in Table 2. Precipitation increased significantly during all seasons except JJA. The mMK tests' results for maximum, minimum temperature and precipitation

PLOS ONE
Climate change and extremes analysis using observed data (1962-2019) over Pakistan are shown in Table 3. For zone 1, the mMK test for precipitation's change shows significantly decreasing trend with value -0.0618. The values of mMK statistic for changes in maximum and minimum temperature are 0.0752 and 0.0762, respectively, indicating increasing trends. The mMK test's results support the Karakoram anomaly during the 1962-1990 given in Table 5, however, the decadal analysis during the 1990-2019 offers complex results (discussed in detail under the Climate Change assessment and extremes' analysis in HKH region). During the last decade, the trend analysis show that precipitation and temperature have decreasing and increasing trends, respectively, at the stations located in the Karakoram region. Results related to climates' extremes are given in  (Table 2) show a significant increase in average temperature during MAM and JJA, increase during DJF, and insignificant decrease during SON. The mMK results for changes in maximum and minimum temperature indicate significantly increasing trend with values 0.0293 and 0.0938, respectively. On the other hand, average precipitation increased significantly throughout in all seasons. The mMK tests' results (Table 3) show a significantly decreasing trend in precipitation with a magnitude of -0.0141.
Interestingly, this is the only region with the increasing number of frost days and decreasing GSL in both time slices. The SU25, TMAXMean, TN90p, TR20 and TX90p increased during The average number of CWD increased from 22 to 26 in the recent period. The standard deviation of CWD has increased from 8 (during 1962-1990) to 11 (during 1991-2019). The CDD (Fig 9) slightly reduced, with insignificant variability. The maximum number of FD0 recently increased from 34 to 47.
Zone 3. The average temperature significantly increased in this region, except a significant decreased for JJA with a value of -0.42. In contrast, the average precipitation increased significantly throughout except an insignificant increase during the DJF. The results of mMK test show decreasing trends in precipitation change and minimum temperature change while an

PLOS ONE
Climate change and extremes analysis using observed data (1962-2019) over Pakistan increasing trend for maximum temperature change with the values -0.0890, -0.0637 and 0.0645, respectively. Both GSL and SU25 increased in the recent period. TMAXMean has decreasing and significantly increasing trend in the first and second-time slice, respectively. The results of TMIN-Mean are in contrast with TMAXMean. TN10P, TN90P, TNN, TNX, TR20 decreased in the recent time period. A significantly decreasing and increasing trend is noted in TX10P and TX90P in the second time slice, respectively. In addition, most of the precipitation extremes indicate an increasing trend during the first time slice except CDD, CSDI and WSDI with decreasing trend. In contrast, most of these extremes decreased except CDD, R10mm, R25mm and WSDI which have increasing trend during the 1991-2021. R95P, R99P and PRCPTOT have decreasing trend in the second time slice.
The average CWD and their standard deviation reduced from 8 and 6 during 1962-1990 to 6 and 3 during 1991-2019. The CDD (Fig 9) drastically decreased from 330 to 170 during 1962-1990 to 1991-2019 with decreased variability in the recent period. The average number of FD0 (Fig 10) reduced from 55 to 34 during 1962-1990 to 1991-2019.
Zone 4. The average precipitation and temperature changes in zone 4 are given in Table 2. The average temperature changed significantly with an increase during DJF, MAM, SON and The average number of CWD (Fig 8) has increased from 7 to 8 annually. The standard deviation of CWD increased from 3 to 5 during the 1991-2019. The maximum number of CWD is 19 during 1991-2019. The maximum number of CDD reduced from 225 to 149 during 1962-1990 and 1991-2019, respectively. No significant change with the average number of zero is observed in the distribution of FD0 (Fig 10).
Zone 5. The average temperature increased significantly except JJA with an insignificant increase. Changes in average precipitation, except MAM are highly significant. There is a significantly decreasing trend in precipitation and maximum temperature while significantly increasing trend for minimum temperature in this region. The values of mMK test statistic for The average number of CWD slightly increased from 15 to 17 annually during 1962-1990 to 1991-2019 (Fig 8). In contrast, the standard deviation of CWD decreased from 8 to 6 during

Climate change and extremes' analysis in HKH region (Present climatic status of Karakoram Anomaly)
The widely known Karakoram anomaly was investigated to see its present status by utilizing the updated observed meteorological data from the North of Pakistan as shown in Table 1. The stations are located between an altitude of 1,250 m a.s.l. at Chilas and 2,317 m a.s.l. at Skardu station. The data was analyzed seasonally to assess the contemporary climate change. The average temperature has an interesting pattern, specifically in winter (DJF), a significant increase at high altitude station (Skardu) of 0.97 and a significant decrease at the lowest station (Chilas) of -0.56 given in Table 4. Almost similar pattern is observed during MAM with the maximum increase at Astore and decrease at Chilas station. During JJA, average temperature decreased at all meteorological stations. Maximum and minimum decrease in the average temperature is noted at Gopis (with the magnitude of -1.42˚C) and Darosh (with magnitude of -0.01˚C), respectively. These changes are highly significant except at Chilas meteorological station. In SON, the average temperature has an increasing trend at all stations other than Bunji, Chilas and Gopis with decreasing trend. Table 5 shows the seasonal trends in precipitation, maximum, and minimum temperature during 1962-1990. During DJF, precipitation has significantly increasing trend except Gopis where the trend is significantly decreasing. Maximum and minimum temperature results show significantly decreasing trend during JJA except Astore with an insignificant decrease in maximum temperature. These results are in line with Karakoram Anomaly during the 1962-1990.

PLOS ONE
Climate change and extremes analysis using observed data (1962-2019) over Pakistan trend. The comparison of these three decades for maximum temperature is shown in Table 7. The station-wise climate extremes in the northern Pakistan are shown in Fig 11. The results show that SU25 increased throughout except Skardu during 1991-2019. TN10P (cold nights) decreased significantly at some of the stations located in the Northern Pakistan during 1991-2019. TN90P (warm nights), TX10P (cold days) TR20 (tropical nights) have decreasing trend during 1991-2019 in most parts of this region. TX90P (warm days) has mixed trend, decreasing at Darosh, Gilgit and Skardu and increasing trend at the remaining stations.

PLOS ONE
Climate change and extremes analysis using observed data (1962-2019) over Pakistan

Discussion
Zone 1 is rich with water resources (snow and glaciers). Therefore, the changes in climate and climate extremes may have significant impacts on water availability in dry summers, the frequency of flooding, and Glacier Lake Outburst Floods (GLOFs) in the future. Zone 2 comprise Monsoon dominated region receiving more than 60% precipitation during Monsoon season [77][78][79]. The seasonal changes in precipitation and temperature may affect the demand and supply of water and ultimately the crop's production. In addition, Zone 3 is the south-west part of Pakistan receiving low amount of rainfall round the year. Precipitation extremes may contribute less to cause flooding in the region. Zone 4 faces shortage of rainfall and have severe rainfall decrease, which can further exacerbate the shortage of water for drinking, agriculture and other purposes. Zone 5 is mainly comprised of Punjab province dominating the agriculture production of the country. According to our results, R99P has increasing trend in this region with the likelihood of intense rainfall and may cause flash flooding in the future.
As UIB is melt water dominated, therefore, the latest information of the cryosphere is important. In addition, changes in climate variables can significantly affect downstream water's availability for millions of inhabitants [44,80]. Most of the HKHs' glaciers are retreating since the mid-19 th century [41,42,81]. In addition, some recent studies found that mass balance of the glaciers in the Karakoram are nearly in balance during 1971-2010 [41,82,83]. Muhammad et al. [76] derived the mass balance of glaciers in the sub-basins of UIB with negative mass balance except Shigar and Hunza in the central and western Karakoram, respectively. In another study, Muhammad et al. [84] investigated the mass balance of glaciers in Astore (a sub basin of UIB) and concluded that the glaciers in the basin (western Himalaya) are nearly in balance during 2000-2016 similar to the neighboring glaciers in the adjacent Karakoram. In addition to the direct glacier observations, climate data also show interesting signals and indirectly indicate glacier's health status. For instance, [44] noted that mean temperature decreased at the stations (Bunji, Gilgit and Skardu) lying in the Karakoram region during 1961-2000. Their analysis shows that a decrease of 1˚C in mean summer temperature can cause a reduction of~20% in summer runoff. Fowler and Archer [44] concluded an expansion in contrast to the widespread retreat in the remaining parts of the HKH region. Latif et al. [85] investigated trends in maximum, minimum temperature in the northern part of Pakistan. They observed that maximum temperature has mixed (upward and downward) trend while the cooling in minimum temperature become stronger in the recent past. Ali et al. [86] investigated temperatures' changes in the UIB by considering meteorological data for two stations. They noted that the mean temperature increased by 0.63˚C in Skardu while decreased by − 0.137˚C in Gilgit in the recent past. Shahzaman et al. [87] concluded in their study that drought has been observed over central and northern Pakistan during 1982-2019 by integrating various data types. Proportion of South Asian population exposer to heatwaves is increasing under the shared socioeconomic scenarios in the twenty first century, however, the northern areas of Pakistan are unaffected presently [88]. The mass balance of the glacier is heavily dependent on the winter accumulation and summer ablation strongly influenced by precipitation and temperature among other factors [89]. This study provides up to date changes in precipitation and temperature in some parts in the northern Pakistan (4 meteorological stations are in the Karakoram region). We found significant increase in winter average precipitation and decrease in summer average temperature particularly at the meteorological stations in the Karakoram region during the 1962-1990 (Table 1). However, the mMK test indicates a decreasing trend (during DJF) in precipitation and increasing trend (during JJA) in temperature in the recent period of 2010-2019 in the Karakoram which is not aligned with the Karakoram Anomaly. Our findings are an indication of glaciers mass loss in future in the Karakoram region. These stations are situated in the valley below the glacier terminating zones, however, an indirect clue can be inferred based on climate change and extreme events regarding future trends. Our results suggest further investigation of glaciers mass balance in the region in future for better understanding of glaciers' health and associated impacts. The downstream water availability in UIB is dominated by snow and glacier melt water particularly in summer, therefore, the climate change assessment can indicate indirect signals for water availability. In addition, climate extremes can provide guidelines to policy makers and relevant stakeholder for future climate trends and meltwater resources in the country.

Conclusions
During DJF and MAM, average temperature increased in all climate zones. A decreasing trend during JJA in zone 1, 3, 4 and during SON in zones 1 and 2 were noted in the average temperature. Precipitation indicated mixed trend in different zones. Mann-Kendal test indicated trend in maximum and minimum temperature except zone-4 where the maximum temperature has no trend. On the average, CWDs increased in all zones except in zone 3 where it decreased. The number of SU25 increased in all climate zone except zone 4. In the Karakoram region, temperature decreased significantly during JJA. Also, precipitation increased significantly in this region except a significant decrease in Gilgit during DJF. MK test suggested increasing trend in precipitation and decreasing trend in maximum and minimum temperature during 1962-1990 and closely in line with Karakoram Anomaly. However, precipitation (in DJF) and temperature (in JJA) have decreasing and increasing trend, respectively, during the most recent period (2010-2019), which do not support the notion of Karakoram Anomaly. These findings suggested that the glacier may lose mass in future rather than gaining or stable conditions. Climate extremes have mixed trends over North Pakistan. For instance, summer days increased, and cold nights have mixed trend. Warm nights decreased except in western Karakoram.
The key takeaways are: 1. A change is noted in the Climate (maximum, minimum temperature and precipitation) across the country during the analysis period. 2. TN10p, TX10p and TR20, CSDI and PRECPTOT have decreasing trend particularly during the 1991-2019 in most parts of the country.
3. Precipitation has increasing trend in DJF while temperature has decreasing trend in JJA during the 1962-1990 at the stations located in the Karakoram region.
4. During the last decade (2010-2019), precipitation and temperature have decreasing and increasing trend, respectively, at the stations located in the Karakoram region.
5. Based on the points in 3 and 4, the glaciers in the Karakoram region may loses mass balance slightly in the near future, which is in contrast to the famous Karakoram Anomaly.