Spatiotemporal Variability in the Hydrometeorological Time-Series over Upper Indus River Basin of Pakistan

This paper investigates the spatiotemporal variability in hydrometeorological time-series to evaluate the current and future scenarios of water resources availability from upper Indus basin (UIB). Mann–Kendall and Sen’s slope estimator tests were used to analyze the variability in the temperature, precipitation, and streamflow time-series data at 27 meteorological stations and 34 hydrological stations for the period of 1963 to 2014. The time-series data of entire study period were divided into two equal subseries of 26 years each (1963–1988 and 1989–2014) to assess the overlapping aspect of climate change acceleration over UIB. The results showed a warming pattern at low altitude stations, while a cooling tendency was detected at high-altitude stations. An increase in streamflow was detected during winter and spring seasons at all hydrological stations, whereas the streamflow in summer and autumn seasons exhibited decreasing trends. The annual precipitation showed a significant decreasing trend at ten stations, while a significant increasing trend was observed at Kohat station during second subseries of the study period. The most significant winter drying trends were observed at Gupis, Chitral, Garidopatta, and Naran stations of magnitude of 47%, 13%, 25%, and 18%, respectively, during the second subseries. The annual runoff exhibited significant deceasing trends over Jhelum subbasin at Azad Pattan, Chinari, Domel Kohala, Muzaffarabad, and Palote, while within Indus basin at Chahan, Gurriala, Khairabad, Karora, and Kalam in the second time-series. It is believed that the results of this study will be helpful for the decision-makers to develop strategies for planning and development of future water resources projects.


Introduction
Tibetan Plateau comprises three major mountainous areas of Asia, that is, Hindukush, Karakoram, and Himalaya (HKH), also known as the "third pole" or "roof of the world" because of the massive volumes of recurrent snow and glacial ice storage in its high-altitude basins [1,2]. e upper Indus basin (UIB) encompasses a huge constituency of hilly areas, and water resources originating from this region are of decisive significance to the interests of Pakistan. e upper Indus river system is of high importance to sustainable water supply for large populations located in the lower reaches of Indus river in Pakistan. Being an agricultural country with heavy population growth, there is a great stress on water resources to meet the food and fiber requirement for the people.
e satisfaction of irrigation system, domestic consumption, and hydropower needs are dependent on water resources from the UIB and its tributaries. e huge river basins of south and southeast Asia rely on summer monsoonal wet regimes downstream, but the UIB is dependent on melted water from glaciers and snow-fed catchments [3]. If anthropogenic interventions causing greenhouse gas emissions into the atmosphere persist, the average global surface temperature may rise from 0.2 to 0.5°C/decade in the next few decades over UIB [4]. is temperature upsurge in the Himalayan region is greater than the average global rise of 0.74°C over the last century [5]. Large-scale warming of the earth surface over the last ten decades or even more is also indicated by several researchers [6][7][8]. Not only did such high-scale warming affect the global circulation patterns but also direct affects occur in local climatic settings with changes in distribution and characteristics of precipitation and temperature [9]. Changes vary in space and time domains as affected by local climatic and topographic settings [10][11][12][13].
ese spatiotemporal variations in the climatic variables have motivated this study in which we aim to assess possible acceleration of climate changes and related hydrological impacts over UIB.
Numerous climate change studies have been carried out over the UIB. For instance, significant decreasing temperature trends were detected during the monsoon season and warming during the premonsoon season [14][15][16]. Moreover, warming winter and cooling summer trends were found over UIB; however, these trends lack a definite pattern of precipitation [17,18]. However, Bocchiola and Diolaiuti [16] argued that winter warming and summer cooling trends were overstated in earlier studies and were only restricted to Gilgit and Bunji stations, respectively. ey reported increasing (insignificant) precipitation trends over the Chitral and Northwest Karakoram regions and a drying pattern over the UIB. Some recent studies have reported general increasing temperature trends during the premonsoon instead during December-February. Latif et al. [19] stated that the decreasing precipitation trends dominated over the UIB spatially and temporally.
Archer [20] and Lutz et al. [21] observed that a rise of 1°C in mean summer temperature over UIB has directly affected the glacier melt and resulted in an increase of 17% in summer runoff for the river Shyok (basin area 65.025 km 2 ) and a 16% increase for the river Hunza (basin area 13.925 km 2 ), respectively. e mean temperature in the UIB has increased over the last century; however, long-term persistent trends (>100 years) have not been detected [15,16]. Also, the statistically significant (p < 0.05) increases in winter maximum temperature of 0.27, 0.55, and 0.51°C per decade were observed at Gilgit, Skardu, and Dir in the UIB. e climate change impacts are projected in a semiarid climate for the upper Indus basin rises of annual mean temperature up to 4.8°C and 18% increase in annual precipitation by the end of the 21 st century [22]. ese findings indicated that warming is more pronounced in the Astore (5.4°C) and Gilgit (5.4°C) river basins compared to the Skardu (4.9°C) river basin in the summer season. Projected changes in precipitation in the Hunza (+19%) and Gilgit (+21%) river basins are quite similar but changes in the Astore (113%) river basin are comparatively large. Moreover, long-term precipitation and temperature series in the UIB exhibited a complex season-dependent spatial correlation structure indicating large differences in climate change as affected by local climatic and topographic settings [13,23,24].
Glacial melt from the Karakoram region dominated the flows into the main Indus river system [25] and thus changes in precipitation and temperature directly affect streamflow discharges. Archer and Fowler [15,23] showed falling river flows from the highest central Karakoram watersheds over the past two decades and suggested that much winter precipitation is going into long-term storage by glacier surges. Khattak et al. [17] showed increasing trends in the winter flows with the increase of mean maximum winter temperature and a negative trend in mean summer flows. Sharif et al. [26] observed some effect of timing and magnitude of flows by climate change in the high elevated areas of the UIB. Flows are affected by initial snowpack conditions by winter precipitation at the beginning of the spring snow melt derived streamflow and trends in temperature that differ seasonally. It is noted that findings of these studies generally rely on time-series consisting of several decades, thus preventing assessment of aspects of acceleration of climate changes and related hydrological impacts. It is expected that the global climate change may have significant impacts on the regime of hydrologic extremes. Hence, the design and management of water resources systems should be adopted to changing hydrologic extremes. However, most of the studies conducted over UIB were restricted to limited number of stations and data lengths. erefore, the number of hydrometeorological stations was extended depending on the available data and coverage of entire Indus, Kabul, and Jhelum river basins was incorporated to investigate the variations in the hydrometeorological time-series.
is study investigates the spatiotemporal variability in the hydrometeorological time-series data and hydrological impacts over UIB by using Mann-Kendall and Sen's slope estimator tests. Moreover, the impacts of local topographical setting and altitude variations on runoff contributions from glacier melt, snow melt, and monsoon precipitation were also evaluated.

Description of Study Area
e Indus river basin is one of the world's largest transboundary river basins with a total drainage area of about 1.08 × 10 6 km 2 [27,28]. e UIB contributes to half of the surface water disposal in Pakistan depending upon the melted water resources from Hindukush-Karakoram-Himalayan (HKH) region. e URB is located within the terrestrial ranges of 33°40ʹ to 37°12ʹ N and 70°30ʹ to 77°, 30ʹ E, in the mountainous ranges of Hindu-Kush, Karakoram, Himalaya, and Tibetan Plateau [12,17]. ese ranges jointly host 11,000 glaciers [29], which make it one of the world's most glaciated areas, with roughly 22,000 km 2 of glacier surface area [30]. e altitude in the UIRB varies from 200 m to 8500 m.a.s.l with an average elevation of 3750 m a.s.l., covering a catchment area of 286,000 km 2 . e research area is considered as the prime source of fresh water for Pakistan and plays a vibrant role in the sustainable economic development of the country. e location of climatic and streamflow stations is presented in Figure 1. e major subbasins of the UIB are discussed in the following paragraphs.

Shyok and Shigar Subbasins of UIB.
e eastern and central part of the Karakoram is covered by the Shyok and Shigar basins, respectively. About 24% of the area of the Shyok river basin is covered with snow [31]. Similarly, onethird of the area of Shigar basin is covered by glaciers including the world's largest glaciers and ice masses and almost 25 to 90% of the area is covered with snow. e basic source of precipitation of these two basins is westerly disturbances during winter and spring seasons followed by the summer monsoon intruding during various intervals [3,19,32,33]. is study used the discharge data of Shyok river at Yugo and Indus river at Shigar for the trend analysis of flows within these two basins.

Astore and Hunza Subbasins of UIB.
e Astore basin is located in the western Himalayan and Hunza basin is in the western Karakoram ranges. e glacial coverage within these basins is less than snow coverage as compared to the Shyok and Shigar basins. e glaciers and permanent ice cover within Hunza basin is 28% and 14%, which exhibits almost 21% and 3% of the total UIB glacial coverage within Hunza and Astore basins, respectively [3,31]. ree high-altitude stations are installed within Hunza basin, that is, Khunjrab, Naltar, and Ziarat, while river discharge is measured at Danyior bridge at Hunza river. ere is only one climatic station (Astore) installed in Astore basin measuring temperature and precipitation data in this basin. Discharge data of Astore river at Doyian are used in this study for the trend analysis installed by Water and Power Development Authority (WAPDA).

Gilgit Subbasin of UIB.
e Gilgit subbasin ranges between 35.8 and 37 E and 72.5 to 74.4 N comprehends eastern part of Hindukush range and drains towards southeast to join Indus river. e discharge of Gilgit river is measured at Gilgit hydrometric station and at the confluence of Hunza and Gilgit river, which is called Alam bridge. e drainage area of this basin incorporates 12000 km 2 with an elevation range from 1481 to 7134 m m.a.s.l. Four climatic stations are installed in this area, that is, Gilgit, Gupis, Yasin, and Ushkore, by the Pakistan Meteorological Department (PMD) and WAPDA. In this study, data at two stations (Gilgit and Gupis) were used from 1960 to 2014. Hasson et al. [3] reported that the Gilgit basin receives maximum precipitation at Ushkore (3151 m) and minimum precipitation at Gilgit (1460 m) due to the westerly disturbances and summer monsoon.

Jehlum Subbasin of UIB.
e Mangla basin is located on the southern slope of the Himalayas with elevation ranging from 300 m to 6282 m m.a.s.l. and has basin area of around 33425 km 2 at Mangla dam.
is dam serves hydropower generation and regulates the flow from Mangla reservoir. About 55% of the area lies in Indian held Kashmir

Data and Methods.
In this study, hydrological timeseries data of 34 stream gauges and meteorological data of 27 stations for the period of 1963 to 2014 were collected from WAPDA and PMD. e information regarding the location of each stream gauge station, area of subbasins, and mean annual streamflow is presented in Table 1.
Moreover, information of mean annual maximum and minimum temperature (T max , T min ) and precipitation (P) of 27 climatic stations is shown in Table 2. Streamflow measurements in the UIB are carried out by WAPDA with the earliest records commencing from 1960. e stream gauges have a broad range of drainage area from 262 km 2 to 286,000 km 2 . e study area contained three major basins, namely, Jhelum, Indus, and Kabul. e locations of different hydrological and meteorological stations are presented in Figure 1.
e hydrometeorological time-series data of entire study period  were divided into two equal subseries, that is, 1963 to 1988 and 1989 to 2014, to analyze the aspects of acceleration of climate change. Mean monthly, seasonal, and annual values of Tmax, T min , P and Q st were dereived from the daily time-series data. To analyze the seasonal variations in the hydrometeorological time-series data, four seasons were defined as winter (December, January, and February (DJF)), spring (March, April, and May (MAM)), summer (June, July, and August (JJA)), and autumn (September, October, and November (SON)) seasons.

Change Detection.
To detect climate impacts and possible climate change acceleration over the past 52 years, Student's t-test was selected for sample means, whereas the "F-test" was selected to assess aspects of variability on timeseries. e nonparametric Mann-Whitney U test was used to assess aspects of the distribution of observations.

Student's t-Test and F-Test.
e application of Student's t-test aims to detect the change in the mean values of two 26-year periods for statistical significance. F-test was used to detect the variability in the time-series. Both tests were performed at the 90% confidence interval. e t-statistic is calculated by using equation (1), when the variances of both time-series have similar values. (1) x 1 and x 2 are mean values of two subseries; n 1 and n 2 indicate number of observations; s 1 and s 2 are the standard deviations; subscripts 1 and 2 indicate the periods 1963-1988 and 1989-2014; and S p is the standard deviation, which is given as If the variances for the two periods are different, then the t-statistic is used and given in the following equation:

Mann-Whitney U-Test.
e nonparametric Mann-Whitney U test [34] is selected to detect the shift in the temperature, precipitation, and streamflow time-series data.
e Mann-Whitney (MW) U test statistic [35,36] is given in the following equation: We have where U 1 and U 2 are the total count of samples 1 and 2; and R 1 and R 2 are the rank sums of sample 1 and sample 2, respectively. When the null hypothesis, H 0 , is true and n 1 and n 2 are both larger than 8, U is considered approximately normally distributed with mean of E(U) and variance of V(U) given as Advances in Meteorology 5

Relative
Changes. e relative change (%) in the annual and seasonal temperature, precipitation, and streamflow was assessed by using the following equation: relative change � mean of 2 nd period − mean of 1 st period mean of 1 st period .

Mann-Kendall Test. A nonparametric rank-based
Mann-Kendall (MK) trend analysis test was used to evaluate the variations in the hydrometeorological timeseries data over UIB [41,42]. e main advantage of the MK test is that there are no assumptions about the statistical distribution of the sample data. Since the method is rank-based, extreme data points in the hydrometeorological time-series will not largely affect the results. e MK test statistic (S) is given by where "n" denotes the length of a dataset and X j and X k are the sequential data values at times j and k. Positive value of S indicates an increasing (upward) trend, and negative value of S reveals a decreasing (downward) trend in the time-series data.
where t k is the number of tied values in the qth group and the sign " " represents the summation of all the tied groups. However, if there are no tied groups in the data, then this may be ignored. After calculating the variance Var(S) from equation (4), the standardized test statistic (Z mk ) value is calculated by using the following equation: A positive value of Z mk indicates an upward trend (i.e., increasing), whereas a negative value indicates a downward trend (i.e., decreasing). e test statistic (S) follows the standard normal distribution, where probability of observing a value higher than the test statistic Z mk is tested under the null hypothesis, H 0, that there is no trend for chosen α-level of significance. H 0 is rejected if the absolute value of Z mk > Z 1−a/2 at the α-level is significant.

Sen's Estimator of Slope.
Sen's nonparametric method [43] was used to estimate the magnitude of trends in the time-series data. e slope of "n" pairs of data can be first estimated by using the following equation: Sen's estimator is the median, Q med , of the N pairs of Q i . In the procedure, N values of Q i are ranked from smallest to largest and Sen's estimator is determined by using the following equation: if N was odd, Q med is tested by a two-sided test at the 100(1 − α)% confidence interval and the true slope may be obtained by the nonparametric test. Data were processed using an Excel macro named MAKESENS [44].

Variability in Temperature. Student's t-test, F-test, and
Mann-Whitney U test were used to detect the percentage change between two consecutive (26-year each) hydrometeorological time-series. Table 3 showed the results of maximum and minimum temperatures (T max and T min ) and precipitation variables at seasonal and annual scale. e results of testing are combined with assessment of relative changes as indicated by percent change. e results showed inconsistent variations in all variables at different stations and seasons. e annual Tmax during winter and spring season decreased (−22% and −27%) between two subseries at Bagh and Naran stations, respectively, whereas at Murree and Peshawar stations, it increased by 10% and 20%, respectively. It was noted that the changes in T max at seasonal scale were quite high in magnitude as compared to T min during the same seasons. For T max , most significant      in T min by using Student's t-test. Results of the F-test and U test also suggested that the climate for the 2 nd period was quite different from the 1 st period with most changes that are statistically significant at 90% confidence level. Significant changes in variance are indicated at most of the stations over UIB in T min (Figure 2). e summary of the trend analyses and the spatial variation in annual, winter, spring (premonsoon), summer (monsoon), and autumn (postmonsoon) maximum and minimum temperature are presented in Figures 3 and 4. Most of the stations exhibit increasing trends in annual maximum temperature for the 1st period. e increasing trends were found at 56% stations, out of which only 4% were significant. Similarly, decreasing trends were found at 44% stations (19% significant). Cherat, Gujar Khan, and Kakul have the highest decreasing rate (1.2, 0.8, and 0.6°C per decade at 99.9%, 99%, and 95% significant level, respectively). e stations show warming trends at a magnitude of 0.1 to 0.5°C·decade −1 . In the 2nd period, increasing trends were observed at 85% stations, out of which 44% stations exhibit significant trends. e decreasing trends were found at Bagh, Gujar Khan, Naran, and Parachinar at the rates of 7.3, 0.3, 0.3, and 0.7°C per decade, respectively, but only Bagh and Parachinar exhibited significant trends at 99.9% and 95% level of confidence interval, respectively. At annual scale, almost all of the stations exhibited warming trends. However, these findings needed further validation to analyze the clear scenarios of climate change acceleration in the study area. erefore, for more detailed trends of climate change, three-month seasonal analysis was carried out and similar warming patterns for winter, spring, and autumn were found during the 1 st period . MK test detected significant trends at 25, 24, and 22 stations at 90%-99.9% significant level in winter, spring, and autumn temperature time-series as shown in Figures 3 and 4. Inconsistent trends were detected during the summer season. e maximum temperature has decreased at 70% stations, out of which 33% stations exhibit significant trends in summer season.
e spring season showed the highest rate of warming as compared to other seasons. More increasing trends were observed during the 1 st period as compared to the 2 nd period. Increasing trends were also revealed at 81%, 89%, 67%, and 84% (15%, 70%, 19%, and 11% significant) during winter, spring, summer, and autumn seasons, respectively. Negative trends were observed at 19%, 11%, 33%, and 26% (11%, 4%, 15%, and 15% significant). Strong indications of climate change acceleration during second period were observed as compared to the first period. In the 1 st period, trend analysis proposed the existence of decreasing trends in annual minimum temperature at 22 stations. Most of the stations exhibited decreasing trends in the annual minimum temperature at 59% stations (41% significant). Only three stations showed significant increasing trends: Bunji, Chilas, and Peshawar. Bunji station showed highest warming trend of 1.5°C per decade. In the 2 nd period, trend analysis with the MK test displayed trend existence at 25 stations. More increasing trends were found as compared to the first period. ese increasing trends were found at 56% (19% significant) and the decreasing trends were found at 44% (7% significant). Trend analysis in seasonal minimum temperature during the 1 st period revealed that winter and spring seasons exhibited more increasing trends at 63% and 67% stations (30% and 22% significant), respectively. Meanwhile in summer and autumn there was a pattern of decreasing trends at 74% and 67% stations (41% and 41% significant), respectively. As compared to the 1 st period winter and summer seasons showed decreasing minimum temperature, while spring and autumn seasons indicated increasing minimum temperature in the 2 nd period . Figures 5 and 6 revealed that at 67% stations (26% significant) and 52% stations (11%) warming trends were found for the winter and autumn seasons. For summer season, 74% stations (67% significant) exhibited cooling trends. e cooling rates varied from 0.1°C·decade −1 to 1.9°C·decade −1 .

Variability in Precipitation.
Significant differences were observed at Naran and Gupis stations for all seasons, but a different pattern was revealed during winter and summer. It was observed that percent change values that are statistically significant are relatively large at few stations with values in the range of +25% to −25%. Highest increase in percent changes of precipitation was detected at Gupis and Naran stations during all seasons but these changes became negative and quite lower at low-altitude stations. e results of analysis by applying Mann-Kendall test and Sen's slope estimator methods in the annual precipitation time-series were summarized for two consecutive 26-year periods, that is, 1963-1988 and 1989-2014. e annual precipitation increased significantly at five stations, while it decreased at four stations during the first period. It was noted that the Gupis station exhibited significant increasing precipitation at the rate of 32% per year with 99% level of confidence. In the 2 nd period at two stations the annual precipitation has increased significantly but decreased at ten stations (Table 3). e highest increasing trend was observed with a magnitude of 47% per year with 99.9% level of significance at Kohat station, while the highest decreasing trend was revealed with a magnitude of 26% at 95% confidence interval at Risalpur station as shown in Figures 7 and 8. In the 1 st period the MK nonparametric test showed negative trends (Figures 7 and 8) in precipitation time-series during winter and autumn seasons at 59% staions (15% significant) and 59% stations (11% significant) and positive trends at 74% stations (11% significant) and 74% stations (30% significant) during spring and summer seasons, respectively. e most significant winter drying patterns were revealed at Gupis, Chitral, Garidopatta, and Naran stations at rates of 47%, 13%, 25%, and 18%, respectively, during the second period. Spring and summer seasons showed decreasing trends at 93% stations (48% significant) and 78% stations (22% significant), respectively, during the 2 nd period. It was observed that 63% stations (11% significant) exhibited increasing trends. In autumn seasons 63% insignificant stations showed decreasing trends as shown in Figures 7 and 8.

Variability in Streamflow.
e annual runoff in Kurram, Soan, and Indus subbasins decreased by 18%, 13%, and 12%, respectively; however, the runoff variations are found to be statistically significant in Indus subbasin. e winter season showed the largest variations compared to other seasons. Moreover, all subbasins showed positive variations during winter season except for Kurram river subbasin as shown in Table 4. e summer flows have been decreased in all rivers. Combined change detection results for Qst at annual base only show small relative changes (−29% to 11%) at most of the stations. ese results suggest that time-series have not notably changed over time. In Swat river a relative change of +22% was observed, which was found to be significant using Student's t-test, F-test, and U test. At seasonal scale, most of the changes are positive in winter season but negative in summer. In winter season, changes are largest (up to 69%) and statistically significant, whereas in summer season changes are negative and statistically insignificant. Chakdara station of Kabul river basin showed most substantial relative change of 69% in winter season, which was significant for all three tests. In summary, analyses on change detection in general indicate acceleration of climate change. e results of annual mean streamflow at 34 stations using MK test of two consecutive 26-year periods are presented in Figures 9 and 10. During the 1 st period, increasing trends were observed at 56% stations (11% significant) and decreasing trends at 44% stations (11% significant). However, only seven stations revealed significant decreasing trends. e highest decreasing trends were revealed at Jhansi Post at a rate of 43% during the 1 st period, that is, 1963-1988. e decreasing trends in annual mean streamflows were found at 77% stations (43% significant) and the increasing trends at 23% stations, which are statistically insignificant. All nine tributaries of Jhelum basin (Naran, Garhi Habibullah, Muzaffarabad, Chinari, Domail, Kohala, Azad Pathan, Kotli, and Palote) exhibited significant decreasing trends. e five tributaries of Indus basin (Karora, Gurriala, Khairabad, Chahan, and Massan) have also shown significant decreasing trends. At Kabul basin only Kalam showed significant decreasing trends. e highest decreasing significant trends were found at Palote station at magnitude of 43% during the whole study period. Winter mean flows have significantly increased at four stations and decreased at seven stations. e highest significant increasing trend was observed at Massan station of Indus river, while decreasing trend was revealed at Jhansi Post station of Kabul river at rates of 15% and 38% for the first and second periods, respectively. All three major rivers exhibited increasing trends at Azad Pattan in Jhelum, Besham in Indus, and Nowshera in Kabul; however, significant trends were detected in the Indus river. During spring season, significant streamflow trends were detected over 10 stations (5 increasing and five decreasing). e Brandu river at Daggar showed significant increasing trend at a rate of 18%, whereas the Bara river exhibited decreasing trend at Jhansi Post station of 39%. Most of decreasing trends were observed in summer and autumn seasons as shown in Figures 9 and 10. In summer and autumn seasons, 57% stations (9% significant) and 60% stations (31% significant) exhibited decreasing trends, respectively. In the second period for winter season mean flows have increased at rate of 54% (14% significant) and decreased at rate of 46% (11% significant) of the data period   basins revealed the significant deceasing trends at the rates of 15%, 23%, and 46%, respectively. Insignificant decreasing trends were found only at Poonch river at the rate of 18%.

Discussions and Conclusions
UIB is a region that is famous for conflicting signals of climate and contrasting hydrological regime [15]. e basic reason for this anomalous behavior is the difference between accumulation patterns of this region as reported by various researchers on the basis of geodetic mass balance and remote sensing data acquisition studies [12,24]. e results of this study predicted that the climate change and acceleration over UIB may seriously affect streamflow in Indus river. ere are primarily three sources of streamflow in UIB, i.e., glacier melt (Hunza, Shigar, and Shyok subbasins) followed by snowmelt (Astore and Gilgit subbasins) and precipitation. e hydrology of the Karakoram and Himalyan ranges is different as reported by various authors. e basic difference is between the accumulation patterns of these two regions, Karakoram mountain receives precipitation by the westerly disturbances, and Himalyan range is controlled by the summer monsoon [23]. Climatic characteristics vary across the Indus basin by large topographic variations from very high elevated to low elevated areas facing Himalaya mountainous slopes. Himalayans cover is one of the most dynamic and complex mountain ranges in the world and is also vulnerable to global warming and increasing human activities. Uncertainties in the rate and magnitude of climate change and potential impacts prevail, but there is no doubt that it is gradually and powerfully changing the ecological and socioeconomic landscape in the Himalayan region, particularly in streamflows. e impact of warmer climate over snow fed basins is opposite to the impact on the glacier fed basins: snow fed basins are more sensitive regarding reduction in the availability of water due to a compound effect of increase in evaporation and decrease in snow melt. e present study investigated the existence of trends and relative changes in the annual and seasonal maximum and minimum temperature, precipitation, and streamflow for two consecutive 26-year periods. e results of this study revealed that climate change is occurring remarkably with warming trends in the lower part of Mangla catchment, whereas cooling trends were observed at the higher elevation regions. e prevailing trends, caused by climate change, influencing the flows should be considered by the water managers for better water management in a water-scarce country like Pakistan.
Most of the river gauges during winter (DJF)) showed the significant increasing river flows during the first-time series. Mukhopadhyay and Khan [45] reported that precipitation is controlled by elevation; precipitation is almost negligible below 2500 m and there is little water yield between 1000 and 2500 m. ey defined "mid-altitude melt" as water that generates due to seasonal snows during intense winter between elevations of 2500 and 3500 m. ese increased river flows are more associated with increased precipitation than temperature during this period. ese results are partially consistent with the findings of Khattak et al. [17] who also reported increased winter flows due to increased precipitation and temperature causing early melting of snow. But we find a significant decrease in river flows during the second, third, and fourth data periods, which is consistent with decreased mean temperature during the same periods.
We observed significant changes in the second period as compared to the first period within UIB, which are consistent and in agreement with the global warming trends reported by Hasson et al. [32]. Our results indicated high variability of trends in annual and seasonal minimum temperature, but these trends appeared to be more significant and higher in magnitude, particularly during the second period. Similarly, we observed significant wetting pattern of trends in annual and seasonal precipitation at higher altitude region, but significant drying trends were revealed by lower-or mid-altitude stations. Most of the stations within UIB exhibited significant drying trends, which are consistent with the recent studies carried out in this region [19]. e annual runoff has been decreased significantly within Jhelum river basin at Azad Pattan, Chinari, Domel Kohala, Muzaffarabad, and Palote and within Indus at Chahan, Gurriala, Khairabad, Karora, and Kalam during the second data period. Similarly, we observed that seasonal runoff has been decreased significantly in all seasons except winter during second data series. e increasing trends of winter runoff are more associated with westerly precipitation as compared to concurrent melting temperature because hydrology of the UIB is dominated by winter precipitation (westerly disturbances) as compared to summer monsoon offshoots. Moreover, the decreasing trends of summer runoff are attributed to decreased melting rate consistent with summer cooling reported by various authors [46,47]. e decreasing trends in summer discharge show least melting rates in summer, resulting in stability of glaciers and consequently positive basin storage. Climate changes occur most noticeably in terms of temperature and precipitation over the UIB according to various authors. Moreover, this study found the spring season to be quite dry, supporting the idea of declining precipitation (reported by numerous studies carried out earlier in this area). Downstream areas in the lower portions of the drainage basin (where most of the population depends on the agriculture) are being affected by decreasing rainfall and its impacts on crop sowing and harvesting times. ere will be more stress on available water resources (which are already scarce) if precipitation does not show any significant upsurge: increased dryness could further stress agricultural production. To avoid this potentially distressing situation from getting worse, water resources management must play an important role to ensure the best utilization of available resources, for example, flood control, building dams and reservoirs, lining of canals and water courses, and conservative surface irrigation (trickle and sprinkler irrigation).
Data Availability e hydrometeorological time-series data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that they have no conflicts of interest.