Observed Trends and Variability of Temperature and Precipitation and Their Global Teleconnections in the Upper Indus Basin, Hindukush-Karakoram-Himalaya

: Having an extreme topography and heterogeneous climate, the Upper Indus Basin (UIB) is more likely to be affected by climate change and it is a crucial area for climatological studies. Based on the monthly minimum temperature (T min ), maximum temperature (T max ) and precipitation from nine meteorological stations, the spatiotemporal variability of temperature and precipitation were analyzed on monthly, seasonal, and annual scales. Results show a widespread signiﬁcant increasing trend of 0.14 ◦ C/decade for T max , but a signiﬁcant decreasing trend of − 0.08 ◦ C/decade for T min annually, during 1955–2016 for the UIB. Seasonally, warming in T max is stronger in winter and spring, while the cooling in T min is greater in summer and autumn. Results of seasonal T max indicate increasing trends in winter, spring and autumn at rates of 0.38, 0.35 and 0.05 ◦ C/decade, respectively, while decreasing in summer with − 0.14 ◦ C/decade. Moreover, seasonal T min results indicate increasing trends in winter and spring at rates of 0.09 and 0.08 ◦ C/decade, respectively, while decreasing signiﬁcantly in summer and autumn at rates of − 0.21 and − 0.22 ◦ C/decade respectively for the whole the UIB. Precipitation exhibits an increasing trend of 2.74 mm/decade annually, while, increasing in winter, summer and autumn at rates of 1.18, 2.06 and 0.62 mm/decade respectively. The warming in T max and an increase in precipitation have been more distinct since the mid-1990s, while the cooling in T min is observed in the UIB since the mid-1980s. Warming in the middle and higher altitude (1500–2800 m and >2800 m) are much stronger, and the increase is more obvious in regions with elevation >2800 m. The wavelet analysis illustrated sporadic inter-annual covariance of seasonal T max , T min and precipitation with ENSO, NAO, IOD and PDO in the UIB. The periodicities were usually constant over short timescales and discontinuous over longer timescales. This study offers a better understanding of the local climate characteristics and provides a scientiﬁc basis for government policymakers.


Introduction
Global climate change receives substantial attention from scholars and practitioners due to its effects in temperature and precipitation [1]. Precipitation is one of the most essential variables due to its critical significance in global climate systems and energy Precipitation variations at different time periods are previously highlighted in the UIB [23,26,27]. The authors of Ref. [26] detected significant increasing trends in winter, summer and annual precipitation during 1961-1999 while contrasting results were detected by [24] during 1980-2006. Furthermore, Ref. [17] observed a precipitation increase in western UIB during 1995-2012. Moreover, Ref. [28] observed a more robust decreasing precipitation at higher elevations. These contrasting responses of precipitation observed by different researchers observed using different time-period data, need to be assessed using the available historical records in the UIB.
On the other hand, The North Atlantic Oscillation (NAO) and El Niño Southern Oscillation (ENSO), Indian Ocean Dipole (IOD) and Pacific Decadal Oscillations (PDO) have broad implications in the regional climate of the UIB. These four climate indices need a better understanding of local Tmax, Tmin and precipitation in this region. Here, NAO is an indicator of oscillations from the Atlantic Ocean (mid-latitude westerlies). At the same time, ENSO and PDO are the influence index from the Pacific Oscillations and IOD as an index from the Indian Ocean, respectively. The 2010 floods in Pakistan were associated to a strong La Niña [29]. The authors of Ref. [29] identified the positive and negative correlation of NAO with winter and summer precipitation, respectively, in the UIB [26]. The authors of Ref. [30] explained the winter precipitation and ENSO connections where they found the stable disparity is between higher winter precipitation in El Niño years and condensed early winter precipitation in La Niña years. Furthermore, Ref. [31] explained the effects of NAO on winter precipitation and Southern Oscillation Index (SOI) on the monsoon in the adjoining Himalayas. He further explained the correlation between temperature and precipitation until the 1960s which weakened due to quasi-biennial oscillations, ENSO and sunspot cycles, in this period.
According to the aforementioned studies, many researchers studied the trends and variability in temperature and precipitation. Yet, there are different contrasting results at different timescales. The identified gaps in previous studies are as follows: First, most of the previous studies focused on the trends either the inter-annual or the seasonal variations. Second, most of the studies were only limited to trends. In this study, the available climate data are employed to assess the trends and magnitudes in monthly, seasonal, and annual scales. The motivation and main reason this study is novel from previously published research are: First, the time series anomalies were assessed for temperature and precipitation in the UIB. Second, start of trend, abrupt change, and significant trends were On the other hand, The North Atlantic Oscillation (NAO) and El Niño Southern Oscillation (ENSO), Indian Ocean Dipole (IOD) and Pacific Decadal Oscillations (PDO) have broad implications in the regional climate of the UIB. These four climate indices need a better understanding of local T max , T min and precipitation in this region. Here, NAO is an indicator of oscillations from the Atlantic Ocean (mid-latitude westerlies). At the same time, ENSO and PDO are the influence index from the Pacific Oscillations and IOD as an index from the Indian Ocean, respectively. The 2010 floods in Pakistan were associated to a strong La Niña [29]. The authors of Ref. [29] identified the positive and negative correlation of NAO with winter and summer precipitation, respectively, in the UIB [26]. The authors of Ref. [30] explained the winter precipitation and ENSO connections where they found the stable disparity is between higher winter precipitation in El Niño years and condensed early winter precipitation in La Niña years. Furthermore, Ref. [31] explained the effects of NAO on winter precipitation and Southern Oscillation Index (SOI) on the monsoon in the adjoining Himalayas. He further explained the correlation between temperature and precipitation until the 1960s which weakened due to quasi-biennial oscillations, ENSO and sunspot cycles, in this period.
According to the aforementioned studies, many researchers studied the trends and variability in temperature and precipitation. Yet, there are different contrasting results at different timescales. The identified gaps in previous studies are as follows: First, most of the previous studies focused on the trends either the inter-annual or the seasonal variations. Second, most of the studies were only limited to trends. In this study, the available climate data are employed to assess the trends and magnitudes in monthly, seasonal, and annual scales. The motivation and main reason this study is novel from previously published research are: First, the time series anomalies were assessed for temperature and precipitation in the UIB. Second, start of trend, abrupt change, and significant trends were assessed for annual and seasonal variation in the whole UIB. Third, accessing the homogeneity and single change-point detection for all stations on annual and seasonal scale. Fourth, the wavelet analysis is used to comprehensively assess the significant periodicity modes and their connections with global climate indices, i.e., ENSO, NAO, IOD and PDO's influence in the regional temperature and precipitation of the UIB. We intend to offer and assess the detailed information on temperature and precipitation variability and change point analysis in both time and space in the UIB during 1955-2016, and its association with atmospheric circulation. Such information can help the government policymakers better understand the local and global climate variation characteristics and provide a theoretical basis for effectively adapting and mitigating climate change in the UIB.

Study Area
The UIB spans from latitude 31 to 37 • N and longitude 72 to 82 • E in the HKH, with the Hindukush region in the west and the TP in the east, the Himalayan and Karakoram ranges are in the south and north respectively ( Figure 1). The UIB is home to the most extensive perennial ice (22,000 km 2 ) outside the polar region [30]. The elevation range within the UIB is extreme, spreading from the second highest point in the world, 8611 m at the K-2 Peak, to just over 569 m in the valleys [32]. Topographic altitude typically greatly affects climatic variables. The northern valley floors of the UIB are thus arid, with an annual precipitation of just 100-200 mm, which rises to 600 mm at an altitude of 4400 m, and glaciological studies indicate annual accumulation peaks of 1500 to 2000 mm at 5500 m for higher altitudes [32,33]. With its minimum in the summer snow melt season (June-September) and its high in the winter snow accumulation period (December-February), the average annual snow cover area in the UIB varies from 10% to 70% [34]. The extreme climate events, e.g., flash floods, GLOFs and droughts, are becoming more common in this region from the last decade making the survival of mountain livelihood very difficult and endangered [35].

Data
Monthly T max , T min and precipitation data of 9 meteorological stations were collected from Pakistan Meteorological Department (PMD) and Water and Power Development Authority (WAPDA) ( Table 1). The spatial distribution of meteorological stations is displayed in Figure 1. Out of the nine stations, three locates were in high-altitude (>2800 m), these stations were installed by WAPDA in 1995, and the data are available from the period 1995-2012. These stations installed by WAPDA are automated weather stations, known as data collection platform (DCPs), that measure snow in millimeters of the water equivalent [36]. The remaining six stations were in the elevation of 1250-2210 m and the time period in these stations was from 1955 to 2016. Climate data before 1955 was missing values in the six stations; therefore, we used the data from 1955 onwards which had no missing values. The authors of Ref. [37] and Ref. [17] also used data from 1960 onwards due to these inconsistencies. However, there were no missing values in the data of three high-altitude climate stations from 1995 to 2012 [17]. The seasons were classified as follows: winter (December, January, February); spring (March, April, May); summer (June, July, August) and autumn (September, October, November). The regional mean T max , T min and precipitation of annual and seasonal time series were calculated by taking the average values of all stations. For regional means, only six stations with data from 1960 to 2016 were considered, because the remaining three stations with data from 1995 to 2012 were affecting the time series.

Autocorrelation Test
Autocorrelation method is a very valuable tool in the study of a climatological time series [38]. Prior to running a pattern evaluation, the time series of T max , T min and precipitation data must be checked for autocorrelation [39]. The findings may be greatly changed by the use of nonparametric tests to identify patterns. Both T max , T min and precipitation time series data were then first checked for the presence of an autocorrelation coefficient [40]. The autocorrelation coefficient value was evaluated at a 95% confidence interval against the null hypothesis, using a two-tailed test. T max , T min and precipitation time series data were then first checked for the presence of an autocorrelation coefficient (r1) at a 5% significance level using a two-tailed test [40].
If r1 falls between the upper and lower confidence level limits (lag-1 ≥ 0.254 or lag-1 ≤ −0.254), the Mann-Kendall test can be applied directly prior to any pre-processing for removal of autocorrelation. However, if there is any autocorrelation persist in the time series, it can be removed using (X 2 − r 1 X 1 , X 3 − r 2 X 2 . . . . X n − r n X n ) [27]. The significant autocorrelations were removed using the pre-whitening method [27].

Mann-Kendall (MK) Test and Sen's Slope Estimator
The non-parametric MK test [41,42] is used to detect meteorological variables statistically significant trends. It requires independent data, whereas this technique is robust against outliers and missing values [27]. This test was used earlier to detect trends in climate related studies [2,4,15,43,44]. This test's null hypothesis suggests no trend in the time series, while the alternate hypothesis indicates an increasing or decreasing trend. We used 95% (±1.96) significance level for statistically significant results. In contrast, the negative values indicate a decreasing trend in the time series [45]. Sen's slope estimator [46] is used to find the magnitude of slope in the meteorological variables [43,44,47].

Sequential Mann-Kendall (SQMK) Test
The SQMK test developed by [48] is used to define the start of a trend, abrupt change, and significant trend time in a time series [49]. The forward (UFK) series and backward (UBK) lines estimate the significant trend and abrupt change. UFK in this study is a standardized variable with a standard deviation of zero mean and function. Hence, its sequential activity fluctuates around the point of zero. UFK is similar to the Z values that start from the first to the last data point. Likewise, UBK's value is computed backward, starting from the endpoint to the first time-series point. This test considers the relative value of all terms in the time series (x 1 , x 2 , . . . , x n ). UFK and UBK's declining values specify a negative trend, whereas the increasing values characterize a positive trend in the time series. The point where both UFK and UBK intersect each other is a possible significant change in that year in the time series. The threshold values in this study are ±1.96 (p = 0.05) significance level [50].

Pettitt's Test
The Pettitt's test [51] is widely used to detect a single change point in continuous hydroclimatological time-series data. This test is most commonly used to detect a significant shift in a climatic time series [3]. In this test, the null hypothesis denotes homogenous data, whereas the alternate hypothesis denotes the occurrence of failure in the time series under consideration [52]. This change point is deemed statistically significant if the computed p value is less than 0.05.

Continuous Wavelet (CWT) and Cross Wavelet Transform (XWT)
Wavelet analysis provides excellent knowledge about different frequency intervals; thus, it is an appropriate tool to study the low-frequency events and the nature of a time series [53,54]. It is a suitable tool for assessing the local variations in power in a time series [55]. The CWT [55] emphasizes the one-dimensional variable (e.g., precipitation). Therefore, we used it in this study to assess the climatic time series in the time-frequency domain [56]. Meteorological studies typically apply the Morlet function as a mother wavelet because it provides a good compromise between the time and frequency resolution [56,57]. The Cross wavelet transform (XWT) is another method to assess geophysical time series data on different levels. It can effectively evaluate the correlation between two signals and display the phase structure and particulars of the signal in the time domain and frequency [58]. More details about the theoretical background of XWT power spectra are explained by [55,56].

Autocorrelation Exclusion
Autocorrelation was applied before using the MK test to ensure any significant serial correlation in temperature and precipitation; results of annual and seasonal autocorrelation at different time lags are presented in Figure 2. The results obtained from autocorrelation for temperature and precipitation showed no significant autocorrelation at the first lag, however, a few significant correlations persist at the later lags. Therefore, to advert the possible effects of autocorrelation on an MK test, all the significant autocorrelations were removed.

Basic Information of T max , T min and Precipitation
Mean annual T max , T min and precipitation in the UIB is displayed in Figure 3. T max in region of elevation below 1500 m varies between 20 and 30 • C, while Gupis, Astore and Skardu in the Hindukush and Himalayan region reaches 15-20 • C. However, stations with a high-altitude >2800 m in Karakoram region show big regional differences, with T max in Khunjerab below 0 • C, while Naltar reach up to 10-15 • C. The changing pattern of T min is obvious owing to the altitude effects. T min in stations of elevations above 2800 m is below 0 • C except Naltar. However, it is 10-15 • C in the stations with low elevations. All the stations with elevation more than 1500 m in the Karakoram and Himalayan region receive precipitation in the form of snow during late autumn, winter and sometimes in the advent of spring as well. Naltar receives a higher precipitation of 800 mm due to both snow and rainfall. However, the valley-based stations with elevation below 1500 m, etc. Gilgit, Bunji and Chilas receive precipitation less than 200 mm.

Spatial Distribution of the Trends in T max and T min 4.3.1. Trends in Monthly T max and T min
Results of Sen's slope and the Mann-Kendall analysis of monthly T max and T min are presented in Tables 2 and 3, respectively. T max is cooling on a regional scale in July-September, with a significant trend of −0.36 • C/decade for the whole UIB in August (Table 2). However, the remaining months mostly have significant warming trends. T max variations also show largely spatiotemporal disparities. The greater significant warming trends usually occur in January and February with a trend of 0.44 and 0.35 • C/decade for Skardu and Gilgit, respectively. While the larger cooling magnitudes are detected in September with significant decreases of −1.63 and −1.70 • C/decade for Naltar and Ziarat, respectively. Compared to T max , the discrepancy patterns of T min in the UIB demonstrate better consistency spatially. The prominent trend for most stations is increasing from January to May, while a downward trend is evident in most stations during the remaining months (Table 3). Geographically, stations in higher elevations generally exhibit larger magnitudes, indicating more rapid warming in this region.

Basic Information of Tmax, Tmin and Precipitation
Mean annual Tmax, Tmin and precipitation in the UIB is displayed in Figure 3. Tmax in region of elevation below 1500 m varies between 20 and 30 °C, while Gupis, Astore and Skardu in the Hindukush and Himalayan region reaches 15-20 °C. However, stations with a high-altitude >2800 m in Karakoram region show big regional differences, with Tmax in Khunjerab below 0 °C, while Naltar reach up to 10-15 °C. The changing pattern of Tmin is obvious owing to the altitude effects. Tmin in stations of elevations above 2800 m is below 0 °C except Naltar. However, it is 10-15 °C in the stations with low elevations. All the stations with elevation more than 1500 m in the Karakoram and Himalayan region receive precipitation in the form of snow during late autumn, winter and sometimes in the advent of spring as well. Naltar receives a higher precipitation of 800 mm due to both snow and rainfall. However, the valley-based stations with elevation below 1500 m, etc. Gilgit, Bunji and Chilas receive precipitation less than 200 mm.  Tables 2 and 3, respectively. Tmax is cooling on a regional scale in July-September, with a significant trend of −0.36 °C/decade for the whole UIB in August (Table 2).  For more information and profound insights of temperature changes at different timescales, trends in annual and seasonal temperature were also calculated. Figure 4 depicts the spatial distribution of the Mann-Kendall test for annual and seasonal T max at 9 meteorological stations in the UIB. Observed results of annual T max indicate clear warming with a significant increasing trend of 0.14 • C/decade in the UIB. However, the same increasing significant trend for whole UIB is evident in winter and spring at the rates of 0.38 and 0.35 • C/decade (Table 4). On seasonal scale, all the stations display widespread warming in both winter and spring ( Figure 4). This implies that warming in the two seasons contributes the most for the annual warming. More significant warming appears in Chilas during winter at a rate of 0.16 • C/decade and in Bunji during spring at a rate of 0.20 • C/decade, all these stations are located in low elevations. The cooling trends are dominated by summer, with the most significant trend of −0.22 • C/decade occurring in Chilas. For autumn, the variation trends of T max were less spatially consistent throughout UIB, with mixed downward and upward trends being detected.    Observed results of annual T min indicate significant decreasing trend in annual at rate of −0.08 • C/decade in the UIB. The same negative significant trend for whole UIB is evident in summer and autumn at the rates of −0.21 and −0.22 • C/decade, except in winter and spring with non-significant positive trends for the whole UIB (Table 4). For T min , winter and spring are in a synchronous condition with T max (Figure 4). However, there are less stations showing increasing trends for T min compared to T max . Summer cooling is also common in T max and T min , but the trend magnitudes are stronger in T min at most stations. No dominated trend is detected for autumn ( Figure 4). T min shows large spatial heterogeneity with a mixed pattern of increasing and decreasing trends. Naltar, Ziarat and Khunjerab demonstrate warming trends of 0.61, 0.58 and 0.85 • C/decade, respectively; indicating warming is more prominent in higher elevations. As mentioned earlier, the UIB has the most snow and ice outside of the Polar Regions. Warming in higher elevation region may cause far-reaching eco-environmental problems such as glacier melting, permafrost degradation and vegetation change [4]. Overall, warming in T max is more pronounced than T min for winter and spring, while the cooling in T min is more pronounced than T max for summer.

Trends in Monthly Precipitation
Precipitation is not changing much with significant trends at a few stations (Table 5). Regional trends of precipitation for the UIB are increasing in all months with significant positive trend of 1.18 and 1.30 mm/decade is observed in June and September respectively. Precipitation at Khunjerab increased significantly in February and December with trends of 5.67 and 8 mm/decade. Such significant trends are also found in Ziarat in July and September at rates of 14.13 and 10.75 mm/decade, respectively. Further, no significant increasing/decreasing trends are observed in the low-and mid-altitude (<2800 m) stations during October and February.

Trends in Annual and Seasonal Precipitation
Results of Mann-Kendall test and Sen's estimator for annual and seasonal precipitation in the UIB are displayed in Figure 4. The annual mean precipitation in the UIB is increasing in all stations, except in Astore, and greater magnitudes occur at high-altitude (>2800 m), such as Naltar and Khunjerab. However, winter has the largest number of stations with significant increasing trends among all seasons. Yet, a precipitation decrease is noticeable only in spring in the region, and a trend of −0.15 mm/decade is observed as a whole. In summer and autumn, most stations have consistently upward trends with only a significant trend of 4.39 mm/decade being observed in Gupis. In spring, an insignificant decrease is observed in Naltar and Astore during spring months, while the same scenario is evident in Khunjerab during summer.

Temporal Variations in Anomaly Time Series for Temperature and Precipitation
Anomalies for annual and seasonal T max and T min are presented in Figure 5. Overall, the anomaly series of T max exhibited increasing trend during 1955-2016, except in summer. A prominent increasing trend occurred during 1996−2016 with positive values, except in 2012 ( Figure 5). The changing patterns in winter, spring and autumn, to a large extent, are very similar to those in a year. The only difference is that warming in autumn is a little earlier since the mid-1980s. The anomaly series in summer indicates a well-balanced fluctuant pattern until 1985, but a cooling tendency thereafter. Observed anomalies of T min are presented in Figure 5. The UIB annually experienced a warming trend from 1955 to the mid-1980s, but a variable cooling trend afterwards. Such variations also occurred in summer and autumn. Here, winter and spring exhibited neither a stable nor a gradual pattern during the entire time period. the anomaly series of Tmax exhibited increasing trend during 1955-2016, except in summer. A prominent increasing trend occurred during 1996−2016 with positive values, except in 2012 ( Figure 5). The changing patterns in winter, spring and autumn, to a large extent, are very similar to those in a year. The only difference is that warming in autumn is a little earlier since the mid-1980s. The anomaly series in summer indicates a well-balanced fluctuant pattern until 1985, but a cooling tendency thereafter. Observed anomalies of Tmin are presented in Figure 5. The UIB annually experienced a warming trend from 1955 to the mid-1980s, but a variable cooling trend afterwards. Such variations also occurred in summer and autumn. Here, winter and spring exhibited neither a stable nor a gradual pattern during the entire time period.  Anomalies for annual and seasonal precipitation are shown in Figure 5. No clear patterns are observed in the UIB for the annual precipitation until the mid-1980s, but an increasing trend is apparent during 1985-2016. Such increasing patterns also exists in winter and summer with the anomalies being dominated by positive values from mid-1980s to 2016. For spring and autumn, the anomaly values keep changing, without any clear trend.

Abrupt Changes in Annual and Seasonal T max , T min and Precipitation in the UIB
The SQMK analysis provides a detailed picture for start of trend, abrupt changes and significant trends as presented in Figure 6. The intersection of a forward (UFK) and backward (UBK) line represents a significant change point. The major significant change points are marked with an arrow. In general, 1955-1995 is the period where the trends in the annual and seasonal T max are decreasing, however, increasing afterwards except a constant decrease in summer during 1955-2016. and summer with the anomalies being dominated by positive values from mid-1980s to 2016. For spring and autumn, the anomaly values keep changing, without any clear trend.

Abrupt Changes in Annual and Seasonal Tmax, Tmin and Precipitation in the UIB
The SQMK analysis provides a detailed picture for start of trend, abrupt changes and significant trends as presented in Figure 6. The intersection of a forward (UFK) and backward (UBK) line represents a significant change point. The major significant change points are marked with an arrow. In general, 1955-1995 is the period where the trends in the annual and seasonal Tmax are decreasing, however, increasing afterwards except a constant decrease in summer during 1955-2016.  The abrupt change points for annual, winter, spring, summer and autumn T max occurred in 1992, 1996, 1991, 2014 and 1992 respectively. The abrupt change and significant trends for annual, winter, spring and autumn occurred during 1990-2000, except for summer. However, the abrupt change and significant trends for T min began in the mid-1970s. Overall, T min showed a decreasing trend earlier for annual, summer and autumn starting from the 1980s to mid-1995. The abrupt change points for annual, winter, spring, summer and autumn T min occurred in 1980, 2015, 2014, 1983 and 1974 respectively. In the UIB, temporal variation in annual and seasonal precipitation were more complex in behavior than T max and T min . Precipitation in the UIB in general displayed a decreasing trend during 1955-1990 for annual, winter and summer. However, there is a general rise in summer from the mid-1990s. Here, the abrupt change points for annual, spring, summer and autumn occurred in 1985,2009,1983 and 2015 respectively. However, winter has two abrupt change points, 1989 and 2015.

Change Point in Annual and Seasonal T max , T min and Precipitation at Various Stations
Change point results for annual and seasonal T max , T min and precipitation are shown in Table 6. For T max , there is an upward trend after the mid-1990s in all stations, except Skardu, during annual, winter and spring. Chilas and Bunji displayed a negative trend after 1985. However, Skardu revealed positive trends in T max annually and seasonally during 1975-1985. For T min , the negative trends for seasonal and annual trends were prominent during 1970-1980 in all stations, except in Astore and Chilas. These both stations displayed positive trends after their respective change point years. Furthermore, precipitation unveiled steady trend patterns in seasonal and annually in various stations. Here, Gilgit and Astore displayed negative trends after change points from 1973 and 1975 respectively. However, Gupis displayed positive trends annually and in summer after 1991 and 1992 respectively.

Relationship between Temperature and Precipitation with Elevation
The elevation dependency of climatic warming was found in many previous studies [4,[59][60][61]. To further investigate the impacts of elevation variation on climate, the trends of annual T max , T min and precipitation in different elevations (<1500 m, 1500-2800 m, >2800 m) are computed (Figure 7). Warming in T max is more pronounced since the mid-1990s, the variations during three time periods (1955-1994, 1995-2012 and 1955-2012) are also analyzed for more details. Compared to regions with altitudes <1500 m, the regional trends for T max , T min and precipitation were greater in regions with elevation of 1500-2800 m at the trends of 0.25 • C/decade, −0.12 • C/decade and 6.84 mm/decade during 1955-2012 (Figure 7). During 1955-1994, T max increased, while precipitation decreased at <2800 m (Figure 7). In addition, the impacts of elevation variation on precipitation were more distinct during 1995-2012. Precipitation decreased in regions with altitudes <2800 m, and increased in higher elevation (>2800 m). Overall, the effects of elevation variation on climate are obvious. Temperature and precipitation in the UIB exhibit obvious changes in different phases. As shown in Figure 7, the regional trends for T min were negative during 1955-2012 in <1500 m. However, clear warming trends of 0.20 and 0.48 • C/decade are observed for T min during 1995-2012 in <1500 m and >2800 m.

Monthly T max , T min and Precipitation Variability Modes
The periodic results of regionally averaged monthly T max , T min and precipitation were analyzed by using CWT spectra in the UIB displayed in Figure 8. Observed wavelet spectra results reveal periodic inter-annual significant oscillations for various periods during 1955-2016 in the UIB. Sporadic significant inter-annual oscillations for T max and T min for the entire period ranged from 8 to 16 months (i.e., 0.8-1.4 years). However, the significant inter-

Monthly Tmax, Tmin and Precipitation Variability Modes
The periodic results of regionally averaged monthly Tmax, Tmin and precipitation were analyzed by using CWT spectra in the UIB displayed in Figure 8. Observed wavelet spectra results reveal periodic inter-annual significant oscillations for various periods during 1955-2016 in the UIB. Sporadic significant inter-annual oscillations for Tmax and Tmin for the entire period ranged from 8 to 16 months (i.e., 0.8-1.4 years). However, the significant inter-annual oscillations for precipitation occurred on a time scale of 8-16 months (0.6-1.3 years)., and one during 2005 at 64 months, (5.3 years). Precipitation showed sporadic but significant oscillations during a <1-year time span.

Monthly Tmax, Tmin and Precipitation Variability Modes
The periodic results of regionally averaged monthly Tmax, Tmin and precipitation were analyzed by using CWT spectra in the UIB displayed in Figure 8. Observed wavelet spectra results reveal periodic inter-annual significant oscillations for various periods during 1955-2016 in the UIB. Sporadic significant inter-annual oscillations for Tmax and Tmin for the entire period ranged from 8 to 16 months (i.e., 0.8-1.4 years). However, the significant inter-annual oscillations for precipitation occurred on a time scale of 8-16 months (0.6-1.3 years)., and one during 2005 at 64 months, (5.3 years). Precipitation showed sporadic but significant oscillations during a <1-year time span.   Figure 9 displays the XWT between T max , T min and precipitation with ENSO, IOD, NAO and PDO during 1955-2016 in the UIB. Both T max and T min share almost the same inter-annual coherence of 8-16 months (i.e., 0.8-1.4 years) with all climate indices. However, precipitation in the UIB shared multiple sporadic significant coherence periods ranging from inter-annual to inter-decadal time scale with ENSO, IOD, NAO and PDO. Among all climate indices, ENSO shared the most intermittent significant oscillation followed by IOD and NAO. However, PDO shared the least significant periods among the listed indices. ENSO had numerous periodic oscillations during 1955-2016 ranging from 8 to 64 months, while, IOD had the same significant coherence from 8 to 40 months. However, NAO shared a significant inter-decadal coherence from 8 to 128 months (i.e., 1-10 years). There is inter-decadal significant coherence in PDO and monthly precipitation in the range of~128 months; however, these significant coherence periods are smaller compared to those in NAO in the corresponding period.

Association between Climate Indices with Monthly T max , T min and Precipitation
NAO and PDO during 1955-2016 in the UIB. Both Tmax and Tmin share almost the same inter-annual coherence of 8-16 months (i.e., 0.8-1.4 years) with all climate indices. However, precipitation in the UIB shared multiple sporadic significant coherence periods ranging from inter-annual to inter-decadal time scale with ENSO, IOD, NAO and PDO. Among all climate indices, ENSO shared the most intermittent significant oscillation followed by IOD and NAO. However, PDO shared the least significant periods among the listed indices. ENSO had numerous periodic oscillations during 1955-2016 ranging from 8 to 64 months, while, IOD had the same significant coherence from 8 to 40 months. However, NAO shared a significant inter-decadal coherence from 8 to 128 months (i.e., 1-10 years). There is inter-decadal significant coherence in PDO and monthly precipitation in the range of ~128 months; however, these significant coherence periods are smaller compared to those in NAO in the corresponding period.

Discussion
With the background of global warming, both Tmax and Tmin show consistent warming at a regional and global scale in recent decades, with the latter warming faster than the former [4,[62][63][64][65]. However, climate change in UIB exhibits incomplete synchronization. This study reveals widespread significant warming for Tmax but cooling for Tmin during 1955-2016 (Figures 4 and 5). These warming and cooling results in Tmax and Tmin are consistent with the results of [22] and [21]. The warming in Tmax is more obvious during 1995-2016 compared to the whole time period, especially in high-altitude stations. Further, warming in Tmax is particularly evident in winter and spring, while cooling in Tmin is much stronger in summer and autumn which supports the findings of [17]. The authors of Ref. Right-pointing arrows indicate that the two signals are in phase, while left is for anti-phase signals. Thick contours denote 5% significance levels against red noise.

Discussion
With the background of global warming, both T max and T min show consistent warming at a regional and global scale in recent decades, with the latter warming faster than the former [4,[62][63][64][65]. However, climate change in UIB exhibits incomplete synchronization. This study reveals widespread significant warming for T max but cooling for T min during 1955-2016 (Figures 4 and 5). These warming and cooling results in T max and T min are consistent with the results of [21,22]. The warming in T max is more obvious during 1995-2016 compared to the whole time period, especially in high-altitude stations. Further, warming in T max is particularly evident in winter and spring, while cooling in T min is much stronger in summer and autumn which supports the findings of [17]. The authors of Ref. [5] Observed winter warming in the neighboring region of Gansu, China. Where, such warming is ascribed to greenhouse gas emission [7]. Rate of deforestation increased during the recent decade, which leads to an increase in surface temperature [66]. The cloud condition, particularly the longwave radiation changes are also causing the temperature increase in the UIB, which has a considerable influence in the near ground surface temperature [67].
Summer cooling was identified for both T max and T min in this paper, which was in line with the findings of [17] in UIB. The authors of Ref. [68] observed summer cooling in surface temperature due to rapid increase in irrigation. Results of monthly T max indicate warming from November to May, both in low-altitude and high-altitude stations. This summer cooling also associates the positive mass balance of glaciers, especially in the Karakoram [9,13,69]. On the other hand, Ref. [14] found no significant mass loss in the glaciers of the Astore basin in the North-West Himalayan area, which is adjacent to the Karakoram region. This might be attributed to changes in the local climate of this location. Furthermore, T max and T min also suggests summer cooling in Astore with increasing precipitation. The authors of Ref. [70] observed increased streamflow in winter and spring, and decreased during the summer and fall seasons in the UIB. This further indicates that the influence of prevailing climatic trends and their influential status in this region, which is mentioned above in associations to the glaciers mass balance.
In the UIB, the increasing trends are observed in annual and seasonal precipitation except spring (Figure 4). These results are in good agreement with previous studies conducted in the UIB [23,27], and in neighboring regions [5,71]. The authors of Ref. [27] observed a mixed pattern of increasing and decreasing trends in the winter precipitation in the UIB, but our results indicate a significant positive trend of 2.3 mm/decade. Winter precipitation in the UIB is transported through the westerly winds. The HKH and TP mountain ranges in the UIB act as a barrier to the westerly; therefore, downpours substantial precipitation that accounts for concentrated snowfall during winters in the UIB [36]. The increasing precipitation indicates that monsoons in the region is becoming stronger. In some regions of South Asia, the monsoons are becoming less frequent but more intense in summer, and this increases the risk of climate extreme events in the region [72]. Precipitation increase (and decrease) throughout winter (spring) is owing to changes in westerly winds; spring drying is primarily related with the strengthening and northward movement of the mid-latitude storm, as well as an increase in the number of dry days [73]. It demonstrates that changes in the time periods may lead to inconsistencies in findings. The ongoing climate changes exert adverse effects on river flows, ecosystem services, biodiversity and glacial melting.
Our results displayed a strong correlation between regional climate with ENSO, NAO, IOD and PDO. T max and T min share almost a same inter-annual coherence of~1-1.4 years with all climate indices. Moreover, precipitation shared significant coherence ranging from inter-annual to inter-decadal time scale. Among them, ENSO and IOD had inter-annual periodic coherences, while NAO and PDO also shared significant inter-decadal periods ranging from~1 to 10 years. A recent study conducted by [74] reveals the impact of ENSO on temperature extremes in Pakistan. He also demonstrates that the frequency of daylight (and nighttime) warm extremes is favored during La Niña events, when the anomalous Pacific Ocean SSTs are warmer than normal at 20 • C. Furthermore, [75] found a significant influence of ENSO in the year-round precipitation in Pakistan. They further observed the influence of AMO and PDO in the UIB. Furthermore, [76] found major correlations between the mean temperature and global teleconnections during spring with NAO and ENSO. They further explained that the NAO might have implications on the certain monsoon months, particularly in August. In addition, another study on the temperature patterns over Pakistan observed a strong correlation of NAO in the spring in central and northern Pakistan; particularly with Tmax in the UIB, while ENSO affects from mid to late spring [77]. The NAO and ENSO play a significant role in winter precipitation in the southern parts of Central Asia. At the same time, the warm phase of the AMO can significantly affect the Indian monsoon rainfall, which affects the southern part of Central Asia [78,79]. Another study explains the NAO and ENSO precipitation signal in the Central Southwest Asia are mostly linked with an intensification through westerlies originating in the Eastern Mediterranean and Middle East regions and moving eastward across the Central Southwest Asia region (including northern Pakistan) during the positive NAO and the warm ENSO phases [80]. Both ENSO and NAO have effects on the winter precipitation in the Central Southwest Asia [81]. The increase in northward moisture transport across the Arabian Sea is believed to be the cause of Pakistan's strong upward trend in rainfall in the central monsoon area [82]. The authors of Ref. [26] observed a strong positive correlation between winter NAO and the winter precipitation in Karakoram and a negative correlation between the monthly NAO index and the summer rainfall at some stations in the UIB. Although we concluded our findings regarding the selected climate indices with the monthly T max , T min and precipitation in the UIB, the general behavior is much more complicated than observed. Therefore, further research is required to investigate the characteristics of regional climatic patterns contemporary with the global atmospheric patterns. This region needs more attention to assess the complex climatic patterns and their teleconnections with large-scale climate oscillations.

Conclusions
Based on the monthly temperature and precipitation datasets from nine meteorological stations in the UIB, the non-parametric Mann-Kendall and Sen's slope estimator are applied to analyze the spatial-temporal variations in T max , T min , and precipitation during 1955-2016. Observed results present different characteristics in both time and space as well as under different altitudes.

•
A widespread warming for T max and a cooling with large spatial heterogeneity for T min are detected during 1955-2016, which is not in a parallel direction to that of global findings on global warming, as T max (T min ) in the UIB increased (and decreased) significantly between 1955 and 2016. • For T min , warming trends are more prominent in stations located at above 2800 m.

•
For T max , warming in winter and spring contributes the most for the warming yeararound, while for T min , cooling in summer and autumn especially in November is much stronger. Summer cooling is identified for both T max and T min .

•
The total precipitation amount does not change significantly in UIB during 1960-2012 against its large inter-annual and inter-decadal variability.

•
On the temporal scale, the warming in T max is more obvious during 1995-2016. The decreasing trends in T min become clear from mid-1980s.

•
Increasing trends for T max and precipitation occur in regions with elevations of 1500-2800 m compared to elevations <1500 m during 1955-2012.

•
The results of wavelet transform illustrated strong correlations between ENSO, NAO, IOD and PDO with monthly T max , T min and precipitation in the UIB. • Significant inter-annual oscillations for T max and T min ranged from~1 to 1.4 years, and~1-6 years for precipitation. In general, the periodicities in the CWT observations were spatially continuous.

•
The XWT results for T max and T min share nearly the same inter-annual coherence of 8-16 months (i.e., 0.8-1.4 years) with all climate indices. However, precipitation in the UIB shared inter-annual to inter-decadal significant coherences. • ENSO shared most intermittent significant periods followed by IOD and NAO. However, PDO shared the least significant periods among the listed indices. • ENSO and IOD had inter-annual periodic coherences, while NAO and PDO shared significant inter-decadal periods ranging from~1 to 10 years.
This study is only limited to the climate data of nine meteorological stations in the region. More studies are suggested to better understand the variability, trends and change points using updated in situ data (1995-2021) from 20 high meteorological stations (DCPs) installed over various altitudes of the UIB. A detailed research on the influence of large-scale climate variability controls is also mandatory to better explore the prevailing influences of climate indices in this region.