Precipitation Variations under a Changing Climate from 1961–2015 in the Source Region of the Indus River

The source region of the Indus River (SRIR), which is located in the Hindukush, Karakoram and Himalayan (HKH) mountainous range and on the Third Pole (TP), is very sensitive to climate change, especially precipitation changes, because of its multifarious orography and fragile ecosystem. Climate changes in the SRIR also have important impacts on social and economic development, as well as on the ecosystems of the downstream irrigation areas in Pakistan. This paper investigates the changes in precipitation characteristics by dividing the daily precipitation rate into different classes, such as light (0–10 mm), moderate (10.1–25 mm) and heavy precipitation (>25 mm). Daily precipitation data from gauging and non-gauging stations from 1961–2015 are used. The results of the analysis of the annual precipitation and rainy day trends show significant (p < 0.05) increases and decreases, respectively, while light and heavy precipitation show significant decreasing and increasing trends, respectively. The analysis of the precipitation characteristics shows that light precipitation has the highest number of rainy days compared to moderate or heavy precipitation. The analysis of the seasonal precipitation trends shows that only 18 stations have significant increasing trends in winter precipitation, while 27 stations have significant increasing trends in summer precipitation. Both short and long droughts exhibit increasing trends, which indicates that the Indus Basin will suffer from water shortages for agriculture. The results of this study could help policymakers cope with floods and droughts and sustain eco-environmental resources in the study area.


Introduction
The Indus Basin is among the world's largest transboundary river basins. The southern part of this basin has subtropical arid and semiarid to temperate subhumid climates on the plains of Punjab and Sindh provinces, while the northern part has high-altitude mountainous highlands with humid climatic conditions. The annual precipitation ranges from approximately 2000 mm on the mountain slopes of the northern part to 100-500 mm in the lowlands of the southern part [1]. The upper part of this basin is the source region of the Indus River (SRIR). The SRIR is a distinctive area that has a as the trends of rainy days and short-and long-duration dry and wet spells. Daily precipitation data from 1961 to 2015 covering the entire SRIR, which includes 50 stations, were collected.

Study Area
The Indus Basin is a transboundary basin with a total area of approximately 1.12 million km 2 that is distributed among four countries (Pakistan, India, China, and Afghanistan). The largest part of this basin lies in Pakistan, while the smallest part lies in Afghanistan. According to Immerzeel [4] and Lutz [30], who studied the source regions of all Asian basins in the Himalayas, the Indus Basin heavily relies on high-altitude water resources. The semiarid climate downstream of the Indus Basin, along with the high demands for water for agriculture, industry and the quickly growing population, has made the Upper Indus Basin a very important region [31]. The upper part of the Indus Basin, which is the SRIR, is selected as the study area to investigate the spatiotemporal characteristics of precipitation. The SRIR is located within the terrestrial range of 67.71-82.4° E and 30.4-37.07° N in the HKH mountainous range and on the TP. Approximately 11,000 glaciers are distributed in these ranges [2], making it one of the world's most glaciated areas [32].
According to different studies on glaciers in the Western Himalayas, glacial mass balances are decreasing and glacier snouts are retreating [33]. The glaciated area of the Kabul River basin decreased by 40% from 1962-2002 [34]. A large number of glaciers in the north and west of Karakoram are decreasing, while some glaciers in the central part of Karakoram are showing an increase in mass balance [35,36]. The elevation in the study area, which was derived from a 90 m resolution digital elevation model (DEM) from the Shuttle Radar Topography Mission (SRTM) project governed by the US National Aeronautics and Space Administration (NASA), varies from 200 to 8069 m above sea level ( Figure 1). The study area was divided into four elevation zones with an altitudinal difference of 1500 m between two successive zones for the zone-wise analysis of precipitation ( Figure 2). Zone 4 has the highest area, while zone 1 has the lowest area. The SRIR is the main source of fresh water for agricultural, industrial and domestic use in Pakistan. The total catchment area of this source region is approximately 439,784 km 2 , and it plays a vibrant part in the development of the country.  The Indus River and its tributaries, i.e., the Jhelum, Chenab, Ravi and Satluj Rivers on the eastern side and the Kabul River on the western side, originate from the SRIR and flow through the Indus Basin down to the Arabian Sea. The Jhelum River originates from Verinag in the southeastern part of Kashmir. The source of the Chenab River is near Bara Lacha Pass in the Zaskar Range. The Ravi River originates from the Kullu hills in Himachal Pradesh. The Satluj River initiates in the Manasarovar Lakes in the western parts of Tibet. In this study, we discuss the SRIR along with its main sub-basins. The Indus River basin is the largest basin in this source region, followed by the Kabul River basin, Satluj River basin, Jhelum River basin and Chenab River basin. These river basins were delineated from the STRM 90 m resolution DEM at the corresponding hydrologic stations, i.e., Nowshera for the Kabul River, Tarbela for the Upper Indus River, Mangla for the Jhelum River, Marala for the Chenab River and Suni for the Satluj River ( Figure 2).

Data Collection
The daily precipitation data from 28 gauging stations for the period of 1961-2015 were collected from the Pakistan Meteorological Department (PMD) and Water and Power Development Authority (WAPDA). The elevations of these stations vary from 283 m at the Mangla gauging station to 2317 m at the Skardu gauging station. These stations cover most of the key area of the SRIR, which lies within the boundaries of Pakistan. These precipitation gauging stations cover a north-south range of approximately 100 km from the Gupis to Astore gauging stations and an east-west range of approximately 200 km from the Skardu to Gupis gauging stations [2]. The missing data accounted for less than 0.001% of the data, and these missing data were approximated by developing simple linear correlations with nearby stations. As the SRIR is a transboundary basin shared by Afghanistan, Pakistan, China and India, the PMD and WAPDA are Pakistani organizations that have installed precipitation gauging stations only within the boundaries of Pakistan. Due to the unavailability of precipitation data from Afghanistan, China and India in their respective parts of the SRIR, Asian Precipitation-Highly Resolved Observational Data Integration Towards Evaluation of Water Resources (APHRODITE) precipitation data from 22 non-gauging stations were used in this study. The Indus River and its tributaries, i.e., the Jhelum, Chenab, Ravi and Satluj Rivers on the eastern side and the Kabul River on the western side, originate from the SRIR and flow through the Indus Basin down to the Arabian Sea. The Jhelum River originates from Verinag in the southeastern part of Kashmir. The source of the Chenab River is near Bara Lacha Pass in the Zaskar Range. The Ravi River originates from the Kullu hills in Himachal Pradesh. The Satluj River initiates in the Manasarovar Lakes in the western parts of Tibet. In this study, we discuss the SRIR along with its main sub-basins. The Indus River basin is the largest basin in this source region, followed by the Kabul River basin, Satluj River basin, Jhelum River basin and Chenab River basin. These river basins were delineated from the STRM 90 m resolution DEM at the corresponding hydrologic stations, i.e., Nowshera for the Kabul River, Tarbela for the Upper Indus River, Mangla for the Jhelum River, Marala for the Chenab River and Suni for the Satluj River ( Figure 2).

Data Collection
The daily precipitation data from 28 gauging stations for the period of 1961-2015 were collected from the Pakistan Meteorological Department (PMD) and Water and Power Development Authority (WAPDA). The elevations of these stations vary from 283 m at the Mangla gauging station to 2317 m at the Skardu gauging station. These stations cover most of the key area of the SRIR, which lies within the boundaries of Pakistan. These precipitation gauging stations cover a north-south range of approximately 100 km from the Gupis to Astore gauging stations and an east-west range of approximately 200 km from the Skardu to Gupis gauging stations [2]. The missing data accounted for less than 0.001% of the data, and these missing data were approximated by developing simple linear correlations with nearby stations. As the SRIR is a transboundary basin shared by Afghanistan, Pakistan, China and India, the PMD and WAPDA are Pakistani organizations that have installed precipitation gauging stations only within the boundaries of Pakistan. Due to the unavailability of precipitation data from Afghanistan, China and India in their respective parts of the SRIR, Asian Precipitation-Highly Resolved Observational Data Integration Towards Evaluation of Water Resources (APHRODITE) precipitation data from 22 non-gauging stations were used in this study. These

APHRODITE Precipitation Data
APHRODITE is a state-of-the-art, high-resolution daily precipitation dataset that was developed by a consortium between the Research Institute for Humanity and Nature Japan and the Meteorological Research Institute of the Japan Meteorological Agency from a dense rain gauge observational network in Asia. APHRODITE precipitation data have been used by Tahir [37] along with the snowmelt runoff model in one sub-basin of the SRIR. In this study, we used the latest and improved versions of the daily dataset for Monsoon Asia (APHRO_MA_V1101 and APHRO_MA_V1101EX_R1) covering where M and AP are the measured and APHRODITE precipitation data, respectively; M and AP are the average values of the measured and APHRODITE precipitation data, respectively. APHRODITE precipitation data are considered perfect if the scores of the statistical indices i.e., BIAS, RMSE, CC, SE and IOA, are 0, 0, 1, 0 and 1, respectively. Data homogeneity was checked using a double mass curve by developing the annual amount from daily precipitation time series [38] at all stations in the study area. The double mass curve method disclosed a straight line at all stations, showing no breakpoints, which confirms the reliability of the precipitation time series. [39,40] also confirmed such results by using the double mass curve method in their studies conducted in the Yellow River basin, China.
For the analysis of different characteristics of precipitation in the SRIR, the daily precipitation rate was divided into different classes, e.g., 0-10 mm day −1 precipitation was described as light precipitation, 10.1-25 mm day −1 precipitation was described as moderate precipitation, while precipitation greater than 25 mm day −1 was described as heavy precipitation. A similar classification of daily precipitation has been used by different researchers for the analysis of precipitation in different regions [12,13,[40][41][42].

Selection of Climatic Variables
In this study, we used different precipitation classes and the number of rainy days of each class as prime indicators to determine the precipitation class that occurred most frequently. Light precipitation (0-10 mm) leads to drought-like conditions, while heavy precipitation (>25 mm) leads to flood-like conditions. We also evaluated the variability in the amount of precipitation during different seasons and the frequency of dry spells.

Trend Analysis
The objective of trend analysis is to determine the changes (increasing, decreasing, or trendless) in the values of climatic variables. Nonparametric methods are mostly used to determine the magnitudes of the trends in climatic data. The Mann-Kendall (MK) test and Sen's slope were used to determine the trends in precipitation in the study area. The trends in the time series of precipitation in the study area were also determined by the trend rate, which represents the slope of the trend. The trend rate m 1 was calculated through linear regression using Equation (6).
where x represents the different precipitation variables from 1961-2015. The significance of m 1 was verified by t-test. A p-value smaller than 0.05 of the associated t-value indicates that the associated m 1 is significant at the 95% confidence level [43]. A positive value of m 1 indicates an increasing trend, while a negative value of m 1 indicates a decreasing trend [44].

Mann-Kendall Trend Analysis Test
In this study, we have employed the Mann-Kendall (MK) test to check the trends in long-term time series of observed and APHRODITE precipitation, annual precipitation, seasonal precipitation and number of rainy days from 1961-2015 in the source region of Indus River. The MK test was also used to see trends in short and long duration dry spells and wet spells. Trends in different classes of precipitation in the study area and its main sub-basins are also determined using the MK test. This is an extensively used and accepted as a common effective technique to assess the existence of statistically significant trends in climatic and hydrologic data [45,46].
The MK test is a nonparametric test used for the identification of trends in climatic time series. The test has very low sensitivity to outliers, and a normal distribution of climatic time series data is not required for this test [47]. The MK test statistics (S) are computed by Equation (7).
where n is the length of the data set, X j and X k are data values at times j and k, and sgn is the sign function that takes on values of −1, 0, and +1. The resultant value of S indicates increasing or decreasing trends in climatic data sets.
where p = tied group, t k = number of observations in the kth group After the variance Var(S) is calculated from Equation (9), the standardized test statistic (Z S ) is then calculated using Equation (10).
The standardized test statistic (Z S ) indicates the significance of the trend. Z S is used to test the null hypothesis, H o if Z S > Z α/2 , where α represents the confidence level. α = 5% is used in this study.

Sen's Slope Test
The magnitudes of the trends in climatic time series data were calculated by Sen's slope test. The slope of the pair of data "n" was first calculated using Equation (11).
where X j and X k are data values at times j and k. The median of the "n" values of Q is Sen's estimator of slope. A positive value of Q indicates an increasing trend, while a negative value indicates a decreasing trend in the climatic time series data. The slopes (Q) of the "n" values were ranked from low to high, and Sen's estimator can be calculated using Equation (12).

Evaluation of the APHRODITE Precipitation Data
APHRODITE precipitation data from non-gauging stations in the eastern and western parts of the study area were used due to the unavailability of station data from these regions. The measured daily precipitation data from gauging stations from 1961-2015 were used for the evaluation of APHRODITE precipitation data. The results of the statistical indices that were calculated using Equation (1) to Equation (5) and used to evaluate the APHRODITE precipitation data are shown in Table 1, while the trends and magnitudes of the average monthly observed and APHRODITE precipitation data, which were calculated using the MK test, are presented in Figure 3.
The average values of the CC and IOA statistical indices were 0.7 and 0.8, which are close to their perfect values of 1. The average BIAS value was 2.3, showing an only 2.3% overestimation of the APHRODITE precipitation data ( Table 1). The average monthly measured and APHRODITE precipitation data have similar trends and magnitudes ( Figure 3) at all stations except for only five (Balakot, Gari Dupata, Gujjar Khan, Khandar and Plandri), where the measured precipitation exhibited a slightly larger decreasing magnitude than the APHRODITE precipitation. The average values of the statistical indices showed good agreement between the measured and APHRODITE precipitation data.

Spatial Distribution of Annual Precipitation Amounts
The SRIR has very distinctive climatic conditions, and the precipitation in this area exhibits great spatial and temporal variations. The distribution of precipitation during different months of a year in the study area and its main sub-basins from 1961 to 2015 are shown in Figure 4. The highest amount of precipitation in the study area (72.8 mm) was received in August, while the lowest amount of precipitation (13.4 mm) was received in November. According to a study conducted by Amin [48] in the northern areas of Pakistan which is part of the SRIR, about 21% of observed rainfall occurred in the month of March while about 37% of the observed rainfall occurred in the month of July and was almost the same in August.

Spatial Distribution of Annual Precipitation Amounts
The SRIR has very distinctive climatic conditions, and the precipitation in this area exhibits great spatial and temporal variations. The distribution of precipitation during different months of a year in the study area and its main sub-basins from 1961 to 2015 are shown in Figure 4. The highest amount of precipitation in the study area (72.8 mm) was received in August, while the lowest amount of precipitation (13.4 mm) was received in November. According to a study conducted by Amin [48] in the northern areas of Pakistan which is part of the SRIR, about 21% of observed rainfall occurred in the month of March while about 37% of the observed rainfall occurred in the month of July and was almost the same in August.

Spatial Distribution of Annual Precipitation Amounts
The SRIR has very distinctive climatic conditions, and the precipitation in this area exhibits great spatial and temporal variations. The distribution of precipitation during different months of a year in the study area and its main sub-basins from 1961 to 2015 are shown in Figure 4. The highest amount of precipitation in the study area (72.8 mm) was received in August, while the lowest amount of precipitation (13.4 mm) was received in November. According to a study conducted by Amin [48] in the northern areas of Pakistan which is part of the SRIR, about 21% of observed rainfall occurred in the month of March while about 37% of the observed rainfall occurred in the month of July and was almost the same in August.  The Jhelum River basin in the study area received the highest contribution of monthly precipitation, with maximum and minimum amounts of 165.3 mm (July) and 25.8 mm (November), respectively. The monthly amount of precipitation in the study area was lowest in the UIRB, with a maximum value in August (31.8 mm) and a minimum value in November (6.5 mm). Figure 4 clearly indicates that there are two clusters of peak monthly precipitation in the study area, i.e., February to April due to westerly air disturbances and July to August due to monsoonal effects.
The spatial distribution of the average annual precipitation in the study area is shown in Figure 5. The average annual precipitation was spatially interpolated through kriging interpolation, which uses semivariogram and data values from sample locations to create an unbiased estimation of the data values in unsampled locations [49]. According to Tabios III and Salas [50] and Phillips [51], kriging is a better interpolation technique than inverse distance weighting (IDW) for precipitation estimation. Different precipitation stations in the study area exhibited great variations in annual precipitation [2,52] due to diverse climatic circumstances, e.g., monsoon, westerly, and orographic effects. Many of the precipitation stations in the middle of the study area received the highest amount of annual precipitation. The average annual amount of precipitation decreased extensively from the center, i.e., Western Himalayas, to the eastern side, i.e., the TP, western side, i.e., Hindukush, and the northern side, i.e., Karakoram, in the study area with a range from 113 to 1963 mm ( Figure 5). The Jhelum River basin in the study area received the highest contribution of monthly precipitation, with maximum and minimum amounts of 165.3 mm (July) and 25.8 mm (November), respectively. The monthly amount of precipitation in the study area was lowest in the UIRB, with a maximum value in August (31.8 mm) and a minimum value in November (6.5 mm). Figure 4 clearly indicates that there are two clusters of peak monthly precipitation in the study area, i.e., February to April due to westerly air disturbances and July to August due to monsoonal effects.
The spatial distribution of the average annual precipitation in the study area is shown in Figure  5. The average annual precipitation was spatially interpolated through kriging interpolation, which uses semivariogram and data values from sample locations to create an unbiased estimation of the data values in unsampled locations [49]. According to Tabios III and Salas [50] and Phillips [51], kriging is a better interpolation technique than inverse distance weighting (IDW) for precipitation estimation. Different precipitation stations in the study area exhibited great variations in annual precipitation [2,52] due to diverse climatic circumstances, e.g., monsoon, westerly, and orographic effects. Many of the precipitation stations in the middle of the study area received the highest amount of annual precipitation. The average annual amount of precipitation decreased extensively from the center, i.e., Western Himalayas, to the eastern side, i.e., the TP, western side, i.e., Hindukush, and the northern side, i.e., Karakoram, in the study area with a range from 113 to 1963 mm ( Figure 5). The annual precipitation and number of stations at different elevations are shown in Figure 6a, while the trends and slope of the annual precipitation calculated using the MK test are presented in Figure 6b. Elevation zone 2 (1500-3000 m) has the highest value, while zone 4 (4500-6000 m) has the lowest value of annual precipitation. Most of zone 2 lies in the Western Himalayas, which provide a barrier against monsoon depressions and receives a large amount of precipitation, while zone 4 lies in the TP, which is a high-altitude arid grassland with an annual precipitation of only 100-300 mm. Annual precipitation exhibits increasing trends in all altitudinal zones, with significant increases in zone 3 and zone 4 with magnitudes of 4 mm/year and 0.84 mm/year, respectively, because of climate change. Yan and Liu [53] reported a significant increase in temperature over the TP at an elevation of 4000 m. The annual precipitation and number of stations at different elevations are shown in Figure 6a, while the trends and slope of the annual precipitation calculated using the MK test are presented in Figure 6b. Elevation zone 2 (1500-3000 m) has the highest value, while zone 4 (4500-6000 m) has the lowest value of annual precipitation. Most of zone 2 lies in the Western Himalayas, which provide a barrier against monsoon depressions and receives a large amount of precipitation, while zone 4 lies in the TP, which is a high-altitude arid grassland with an annual precipitation of only 100-300 mm. Annual precipitation exhibits increasing trends in all altitudinal zones, with significant increases in zone 3 and zone 4 with magnitudes of 4 mm/year and 0.84 mm/year, respectively, because of climate change. Yan and Liu [53] reported a significant increase in temperature over the TP at an elevation of 4000 m.

Trends in Annual Precipitation, Seasonal Precipitation, and Rainy Days
The precipitation observed at all stations was averaged to illustrate the temporal distribution of precipitation in the SRIR [54,55]. The temporal distribution of the annual amount of precipitation and the total number of rainy days in the study area and its main sub-basins are presented in Figure 7. The maximum and minimum amounts of annual precipitation in the study area were 753.33 mm in 2015 and 344.39 mm in 2001, respectively. This large difference between the maximum and minimum amount of annual precipitation in the study area was due to the drought, which started in 1999 and ended in 2002 [56]. The annual precipitation exhibited a significantly (P < 0.05) increasing trend in the study area and its four sub-basins, i.e., the Kabul, Chenab, Jhelum and Upper Indus River basins. The Satluj River basin exhibited a nonsignificant decreasing trend. Rainy days exhibited significantly decreasing trends in the study area and all of its sub-basins except the Kabul basin, which exhibited a nonsignificant decreasing trend (Figure 7).
Seasonal precipitation in the study area was calculated by taking into account the standard months of each season. The four seasons in this region were defined as December to February for winter, March to May for spring, June to September for summer and October to November for autumn. The trends in seasonal precipitation, annual precipitation and rainy days that were calculated using the MK test in the study area are presented in Figure 8. During the winter, 35 out of the 50 stations showed increasing trends in precipitation with significantly increasing trends (α = 5%) at 18 stations. Decreasing trends were observed at 15 stations with one station (India 5) showing a significant decrease. Significantly increasing trends were observed at 27 stations during the summer. The significantly increasing trends in precipitation at 54% of the stations during the summer and 36% of the stations during the winter will increase the chances of flooding in the Indus Basin. The eastern part of the study area mostly exhibited decreasing trends, while increasing trends in precipitation were observed during the winter, spring, summer, and autumn in the western part. (a) Elevation-dependent precipitation, (b) trends and slope of precipitation at different elevations.

Trends in Annual Precipitation, Seasonal Precipitation, and Rainy Days
The precipitation observed at all stations was averaged to illustrate the temporal distribution of precipitation in the SRIR [54,55]. The temporal distribution of the annual amount of precipitation and the total number of rainy days in the study area and its main sub-basins are presented in Figure 7. The maximum and minimum amounts of annual precipitation in the study area were 753.33 mm in 2015 and 344.39 mm in 2001, respectively. This large difference between the maximum and minimum amount of annual precipitation in the study area was due to the drought, which started in 1999 and ended in 2002 [56]. The annual precipitation exhibited a significantly (p < 0.05) increasing trend in the study area and its four sub-basins, i.e., the Kabul, Chenab, Jhelum and Upper Indus River basins. The Satluj River basin exhibited a nonsignificant decreasing trend. Rainy days exhibited significantly decreasing trends in the study area and all of its sub-basins except the Kabul basin, which exhibited a nonsignificant decreasing trend (Figure 7).
Seasonal precipitation in the study area was calculated by taking into account the standard months of each season. The four seasons in this region were defined as December to February for winter, March to May for spring, June to September for summer and October to November for autumn. The trends in seasonal precipitation, annual precipitation and rainy days that were calculated using the MK test in the study area are presented in Figure 8. During the winter, 35 out of the 50 stations showed increasing trends in precipitation with significantly increasing trends (α = 5%) at 18 stations. Decreasing trends were observed at 15 stations with one station (India 5) showing a significant decrease. Significantly increasing trends were observed at 27 stations during the summer. The significantly increasing trends in precipitation at 54% of the stations during the summer and 36% of the stations during the winter will increase the chances of flooding in the Indus Basin. The eastern part of the study area mostly exhibited decreasing trends, while increasing trends in precipitation were observed during the winter, spring, summer, and autumn in the western part.
Significantly increasing trends in annual precipitation were observed at 44% of the stations, while significantly decreasing trends were observed at only 4% of the stations, which were India1 and India6. On the other hand, for the number of rainy days, the majority of the stations (58%) showed significantly decreasing trends, whereas only 6% of the stations showed significantly increasing trends (Figure 8). Annual precipitation exhibited increasing trends while the number of rainy days exhibited decreasing trends, showing that a large amount of precipitation will be received in a small number of days, which will increase the chances of flash floods.

Characteristics of Different Precipitation Classes
Daily precipitation data from 50 stations from 1961 to 2015 were divided into light, moderate and heavy precipitation. This study investigated the different characteristics of these precipitation classes. Table 2 shows the trends in the different precipitation classes. Light precipitation exhibited significantly decreasing trends in the study area, moderate precipitation exhibited nonsignificant increasing trends, while heavy precipitation exhibited significantly increasing trends. The heavy precipitation trend increased as the number of rainy days decreased, which means that there will be heavy precipitation for a limited number of days, increasing the chances of floods. Ahmad [57] also reported more severe changes (increasing trend) in heavy precipitation during winter and summer in his study in the UIRB. The contribution of each precipitation class was calculated by taking the average accumulated precipitation of that class to the total seasonal precipitation. The contributions of light, moderate, and heavy precipitation in the study area were 48%, 28%, and 24%, respectively. More than 50% of seasonal precipitation in autumn was light precipitation. The highest contribution of heavy precipitation was in summer, while the lowest contribution was in autumn. The Satluj and Upper Indus sub-basins had the highest contributions from light precipitation, i.e., 76-89% and 59-70%, respectively, while the lowest contributions were from heavy precipitation, i.e., 2-5% and 10-19%, respectively (Table 3). Table 3. Contributions of different precipitation classes. Light  45  46  46  55  Moderate  30  30  25  27  Heavy  25  23  29  18   Kabul   Light  42  42  49  56  Moderate  34  33  24  27  Heavy  25  25  27  17   UIRB   Light  70  59  64  70  Moderate  18  23  19  20  Heavy  12  19  17  10   Jhelum   Light  32  36  32  43  Moderate  33  36  26  34  Heavy  35  29  42  23   Chenab   Light  39  48  67  63  Moderate  35  31  24  26  Heavy  26  21  10  12   Satluj   Light  76  83  79  89  Moderate  21  15  16  8  Heavy  3  2  5  3 The distributions of wet or rainy days for light, moderate, and heavy precipitation in the study area and its main sub-basins are shown in Figure 9. The average number of rainy days for light, moderate, and heavy precipitation are 121, 9 and 4, respectively. In the study area, the Jhelum River basin has the highest number of rainy days for heavy and moderate precipitation, i.e., 10 and 21, respectively, while the UIRB has the lowest number of rainy days for heavy (1 day) and moderate (3 days) precipitation. The distributions of wet or rainy days for light, moderate, and heavy precipitation in the study area and its main sub-basins are shown in Figure 9. The average number of rainy days for light, moderate, and heavy precipitation are 121, 9 and 4, respectively. In the study area, the Jhelum River basin has the highest number of rainy days for heavy and moderate precipitation, i.e., 10 and 21, respectively, while the UIRB has the lowest number of rainy days for heavy (1 day) and moderate (3 days) precipitation.

Trend Analysis of Wet and Dry Spells
Precipitation is very important for the restoration of ecosystems and for maintaining the atmospheric balance. According to De Groeve [58], heavy floods will hit Pakistan after a period of three years. Most of the precipitation in the study area took place as light and heavy precipitation. According to Iqbal [40], precipitation with a duration of 1-2 days was demarcated as short-duration wet spells, and precipitation with a duration of more than 2 days was defined as long-duration wet spells. Gong [41] used short-and long-duration precipitation to examine the continuation of precipitation events. The results of the trend analysis of short-and long-duration wet and dry spells are demonstrated in Figure 10. Short-duration wet spells exhibited increasing trends at most of the precipitation stations (80% of the stations), with significantly increasing trends at 48% of the stations. Long-duration wet spells exhibited decreasing trends at 82% of the stations, with 56% of the stations showing significantly decreasing trends ( Figure 10). Long periods with continuous no precipitation will turn the land into desert and stress the ecosystem. Jabeen [59] reported the existence of drought conditions after a period of six years in Pakistan. Periods with no precipitation for less than 10 days

Trend Analysis of Wet and Dry Spells
Precipitation is very important for the restoration of ecosystems and for maintaining the atmospheric balance. According to De Groeve [58], heavy floods will hit Pakistan after a period of three years. Most of the precipitation in the study area took place as light and heavy precipitation. According to Iqbal [40], precipitation with a duration of 1-2 days was demarcated as short-duration wet spells, and precipitation with a duration of more than 2 days was defined as long-duration wet spells. Gong [41] used short-and long-duration precipitation to examine the continuation of precipitation events. The results of the trend analysis of short-and long-duration wet and dry spells are demonstrated in Figure 10. Short-duration wet spells exhibited increasing trends at most of the precipitation stations (80% of the stations), with significantly increasing trends at 48% of the stations. Long-duration wet spells exhibited decreasing trends at 82% of the stations, with 56% of the stations showing significantly decreasing trends ( Figure 10). Long periods with continuous no precipitation will turn the land into desert and stress the ecosystem. Jabeen [59] reported the existence of drought conditions after a period of six years in Pakistan. Periods with no precipitation for less than 10 days were considered short-duration dry spells (short drought), while periods with no precipitation for more than 10 days were considered long-duration dry spells (long drought). were considered short-duration dry spells (short drought), while periods with no precipitation for more than 10 days were considered long-duration dry spells (long drought). It is very alarming that short droughts, i.e., short-duration dry spells and long droughts, i.e., long-duration dry spells both exhibited increasing trends. Approximately 58% of the stations showed increasing trends in short-duration dry spells, with 8% of the stations showing significantly increasing trends, whereas 78% of the stations showed increasing trends in long-duration dry spells with significant increasing trends at 48% of the stations (Figure 10).

Conclusions
This study focused on characteristics and trends in precipitation and was carried out in the SRIR using daily precipitation data from 1961 to 2015. The SRIR is a transboundary basin, the major part of which is in Pakistan, followed by India, China, and Afghanistan. The water for agriculture, industry, and human consumption in the Indus Basin comes from this source region. Data from 28 gauging stations were collected from the PMD and the WAPDA, which covers only the Pakistani part of this source region. Due to the unavailability of precipitation data from other countries, APHRODITE precipitation data from 22 non-gauging stations (daily data, 1961-2015) were used for the analysis of different characteristics of precipitation in the study area along with the data from gauging stations. The validity of the APHRODITE precipitation data was checked by comparing it with the observed precipitation data from the gauging stations. The results showed that the APHRODITE precipitation data strongly corresponded to the measured precipitation data, as the values of the statistical indices were close to their perfect values.
Annual precipitation exhibited a significantly increasing trend, suggesting that the study area will provide enough surface water for agriculture in the Indus Basin if properly managed. The annual maximum amount of precipitation in the study area was recorded during 2015 and was 753.33 mm, while the minimum amount of annual precipitation was 344.39 mm, which was recorded in 2001. The results of the trend analysis for annual precipitation showed increasing trends in the Kabul, Chenab, Jhelum and Upper Indus sub-basins, which were significant at the 95% confidence level. The rainy days exhibited a decreasing trend, meaning that there will be precipitation for a limited number of days of a year, which increases the chances of flash floods in the study area. Most of the stations in the study area showed a significantly increasing trend in the winter and summer seasons, while nonsignificant increasing and decreasing trends were observed during the autumn and spring seasons. It is very alarming that short droughts, i.e., short-duration dry spells and long droughts, i.e., long-duration dry spells both exhibited increasing trends. Approximately 58% of the stations showed increasing trends in short-duration dry spells, with 8% of the stations showing significantly increasing trends, whereas 78% of the stations showed increasing trends in long-duration dry spells with significant increasing trends at 48% of the stations (Figure 10).

Conclusions
This study focused on characteristics and trends in precipitation and was carried out in the SRIR using daily precipitation data from 1961 to 2015. The SRIR is a transboundary basin, the major part of which is in Pakistan, followed by India, China, and Afghanistan. The water for agriculture, industry, and human consumption in the Indus Basin comes from this source region. Data from 28 gauging stations were collected from the PMD and the WAPDA, which covers only the Pakistani part of this source region. Due to the unavailability of precipitation data from other countries, APHRODITE precipitation data from 22 non-gauging stations (daily data, 1961-2015) were used for the analysis of different characteristics of precipitation in the study area along with the data from gauging stations. The validity of the APHRODITE precipitation data was checked by comparing it with the observed precipitation data from the gauging stations. The results showed that the APHRODITE precipitation data strongly corresponded to the measured precipitation data, as the values of the statistical indices were close to their perfect values.
Annual precipitation exhibited a significantly increasing trend, suggesting that the study area will provide enough surface water for agriculture in the Indus Basin if properly managed. The annual maximum amount of precipitation in the study area was recorded during 2015 and was 753.33 mm, while the minimum amount of annual precipitation was 344.39 mm, which was recorded in 2001. The results of the trend analysis for annual precipitation showed increasing trends in the Kabul, Chenab, Jhelum and Upper Indus sub-basins, which were significant at the 95% confidence level. The rainy days exhibited a decreasing trend, meaning that there will be precipitation for a limited number of days of a year, which increases the chances of flash floods in the study area. Most of the stations in the study area showed a significantly increasing trend in the winter and summer seasons, while nonsignificant increasing and decreasing trends were observed during the autumn and spring seasons.
The precipitation data from all the stations were divided into different classes, i.e., 0-10 mm, 10.1-25 mm and >25 mm classes. The light precipitation class (0-10 mm) exhibited decreasing trends in the SRIR and all its sub-basins, with significant decreases in the SRIR, Jhelum and Satluj sub-basins. The heavy precipitation class (>25 mm) exhibited increasing trends in the SRIR and all its sub-basins, which was significant at α = 5%, which indicates increased flood risks in the study area. The precipitation in the study area was mainly contributed by light precipitation, followed by moderate and heavy precipitation.
Precipitation with a duration of one to two days exhibited increasing trends at 40 stations, of which 24 stations exhibited significant increasing trends. Precipitation with a duration of more than two days showed decreasing trends at 41 stations, of which 28 stations showed significantly decreasing trends, indicating that most days in a year will be dry, which will convert the land into desert. Analysis of droughts was also carried out in the study area. Periods of less than 10 days without precipitation were considered short droughts, and periods without precipitation for more than 10 days were considered long droughts. Short droughts exhibited increasing trends at 29 stations, while long droughts exhibited increasing trends at 39 stations out of the total 50 stations in the study area.
From the findings of this study, it is concluded that precipitation in the study area is mainly contributed by the light class followed by the moderate and heavy classes. The light precipitation class exhibited decreasing trends, whereas the heavy precipitation class exhibited increasing trends, which increased the chances of floods in the region. Periods with no precipitation also exhibited increasing trends, which indicates increased risks of droughts in the future. On the basis of these findings, it is recommended to construct more structures, i.e., water storage reservoirs, to cope with floods and droughts. The use of water-efficient technologies is also recommended in the Indus Basin to cope with shortages of irrigation water in the scenario of increasing droughts in the future in the source region of the Indus River.