Long-Term Hydro–Climatic Trends in the Mountainous Kofarnihon River Basin in Central Asia

: Hydro–climatic variables play an essential role in assessing the long-term changes in streamﬂow in the snow-fed and glacier-fed rivers that are extremely vulnerable to climatic variations in the alpine mountainous regions. The trend and magnitudinal changes of hydro–climatic variables, such as temperature, precipitation, and streamﬂow, were determined by applying the non-parametric Mann–Kendall, modiﬁed Mann–Kendall, and Sen’s slope tests in the Kofarnihon River Basin in Central Asia. We also used Pettitt’s test to analyze the changes during the 1951–2012 and 1979–2012 time periods. This study revealed that the variations of climate variables have their signiﬁcant spatial patterns and are strongly regulated by the altitude. From mountainous regions down to plain regions, the decadal temperature trends varied from − 0.18 to 0.36 ◦ C / decade and the variation of precipitation from − 4.76 to − 14.63 mm yr − 1 per decade. Considering the temporal variation, the temperature trends decreased in winter and signiﬁcantly increased in spring, and the precipitation trends signiﬁcantly decreased in spring but signiﬁcantly increased in winter in the high-altitude areas. As consequence, total streamﬂow in headwater regions shows the obvious increase and clear seasonal variations. The mean monthly streamﬂow decreased in fall and winter and signiﬁcantly increased in the spring and summer seasons which can be attributed to the inﬂuence of global warming on the rapid melting of snow and ice. Although the abrupt change points in air temperature and precipitation occurred around the 1970s and 1990s in the low-altitude areas and 2000s in the high-altitude areas during the 1951–2012 and 1979–2012 periods, the general trends of hydro–climatic variables keep consistent. This study beneﬁts water resource management, socio–economic development, and sustainable agricultural planning in Tajikistan and its downstream countries. in fall and The modiﬁed MK test result exhibited a signiﬁcant trend in the annual streamﬂow in the upstream of the KRB. The streamﬂow and in the and glacier-fed mountainous KRB of the and We detected change points in the upstream ﬂow in 1991 during both the 1951–2012 and 1979–2012 This study analyzed only the hydro-climatic spatiotemporal variability and did not consider the impact of anthropogenic activities. In the future, relevant indicators should be selected for further study, and the impact of anthropogenic activities should be analyzed. We revealed that the hydro–climatic variability in the upstream areas of the KRB impacts downstream water ﬂuctuations. Our ﬁndings in this study will provide a scientiﬁc theoretical basis for water issues under ongoing climate change and help us to improve the water resources management in the KRB of Central Asia.


Introduction
The mountainous Kofarnihon River Basin (KRB) in Central Asia has a strong local contrast due to its mountainous topography and is considered the most vulnerable territory in Central Asia to

Study area
The Kofarnihon River Basin is located in Central Asia and lies between 37 • and 39 • N and 68 • and 70 • E. It is one of the northwest tributaries of the Amu Darya River (ADR) and the transboundary river between Tajikistan and Uzbekistan. Before confluence with the ADR, the total length of the Kofarnihon River (KR) is 387 km. The total area of the KRB is 11,590 km 2 . The city of Vahdat and the capital city of Tajikistan, Dushanbe, are located on the KRB. Water is used for drinking in these well-populated cities, for irrigation of the lands of the valley, and for fishing. Kofarnihon River irrigates 44.958 hectares of land in Tajikistan [3] and 70% of the KRB area is situated in the mountainous region. The main tributaries of this river are Varzob and Khonaqo in the northwest and Elok in the northeast. Elevation in the basin ranges from 304 m to 4830 m above sea level. The climate of Tajikistan is determined by its geographical position within the continent of Eurasia on the verge of subtropical and temperate zones. Its characteristic features are a high intensity of solar radiation, aridity, low cloud cover, a long duration of sunshine, sharp fluctuations in daily and seasonal temperatures, and significant dustiness of the air. However, these features manifest themselves differently in individual regions due to the peculiarities in the orographic structures and diversity of the heights of the Earth's surface. In the cold season, Water 2020, 12, 2140 4 of 25 a polar front is generally located over Tajikistan and Central Asia. Weather conditions form under the influence of cold dry air arriving in the form of cyclones from the Atlantic Ocean. The climate of the KRB is continental under the domination of the westerlies, with very high local contrast due to the geographical terrain. The KRB is occupied by a mountainous climate in the most of its territory, which is characterized by moderate winters in hilly dissected terrain areas, cold winters in mountainous areas, and summer seasons with relatively large annual temperature variations. With an increase of elevation, the temperature decreases, and at different altitudes and features, the precipitation has different patterns. The location of the KRB is shown in Figure 1.
precipitation has different patterns. The location of the KRB is shown in Figure 1.
The average water flow in the KR is 164 m 3 /s. The minor management of water flow in the KR begins in late February and continues until July, when the flow is at its maximum. The river flow starts to decline in July and continues until January when the river flow is in its lowest condition. The flooding period mostly starts in the first 10 days of March and ends in the third 10 days of September. The river flow from July to September is relatively insignificant. In the fall (October-November) and spring (March-May) seasons, the water flow in the KR is stable. Fall and spring are the two growing seasons in Tajikistan, when precipitation is more intensive. The winter months in the mountains of the KRB are a period of precipitation accumulation in the form of snow mass. The low water flow period occurs in the August and September months when the melting of snow concludes. The conditions for ice melting depend on the climatic features in the relevant period of time. Ice melting can mostly be observed in the areas of the Anzob and Faizobod meteorological stations. In these territories, the climate is not sufficiently humid and has moderately warm summers and mild winters. Isambay station is located in the middle reaches of the KR. In the territory of the Isambay meteorological station in the middle stream of the KR, the climate is not humid enough; the area has very warm summers and mild winters. The average water flow in the KR is 164 m 3 /s. The minor management of water flow in the KR begins in late February and continues until July, when the flow is at its maximum. The river flow starts to decline in July and continues until January when the river flow is in its lowest condition. The flooding period mostly starts in the first 10 days of March and ends in the third 10 days of September. The river flow from July to September is relatively insignificant. In the fall (October-November) and spring (March-May) seasons, the water flow in the KR is stable. Fall and spring are the two growing seasons in Tajikistan, when precipitation is more intensive. The winter months in the mountains of the KRB are a period of precipitation accumulation in the form of snow mass. The low water flow period occurs in the August and September months when the melting of snow concludes. The conditions for ice melting depend on the climatic features in the relevant period of time. Ice melting can mostly be observed in the areas of the Anzob and Faizobod meteorological stations. In these territories, the climate is not sufficiently humid and has moderately warm summers and mild winters. Isambay station is located in the middle reaches of the KR. In the territory of the Isambay meteorological station in the middle stream of the KR, the climate is not humid enough; the area has very warm summers and mild winters.

Data and Materials
In this study, we utilized the monthly temperature and precipitation from six meteorological stations and monthly river flow data from one hydrological station in the KRB of Central Asia ( Table 1). The monthly data of the temperature and precipitation were supplied from the Agency of Hydrometeorology of the Committee for Environmental Protection under the Government of the Republic of Tajikistan. The hydrological data for the Dahana station were derived from the Ministry of Energy and Water Resources of the Republic of Tajikistan. A 30 m spatial resolution of the Digital Elevation Model was used from the website of the Shuttle Radar Topographic Mission to demarcate the catchment boundary and river network in the KRB [55]. The land-use type of the KRB at a 300 × 300 m grid-scale was computed from the Envisat Medium Resolution Imaging Spectrometer [56]. The land cover types for the basin are water body (0.01%), forest (0.82%), urban area (1.55%), snow and ice (1.94%), herbaceous cover (4.75%), bare land (10.27%), shrubland (11.29%), grassland (33.76%), and agriculture land (35.62%). For the KRB, the soil type was determined from version 1.2 of the Harmonized World Soil Database (HWSD) with a scale of 1:5,000,000. The HWSD is freely available on the website of the Food and Agriculture Organization of the United Nations [57]. In this study, the most prominent soil types were cambisols (48.69%), leptosols (26.60%), calcisols (9.79%), anthrosols (8.37%), shifting sand dunes (3.53%), and arenosols (3.02%). Mountain gorges, with parents based on soil, include terraces of alluvial deposits, while valleys mainly feature loess rocks, which reach depths of 100 m and below. Throughout the Kofarnihon River, soil ranges change (from mountain brown to light gray soils and gray-brown soils) with a decrease in the elevation.
In Central Asia, the complexity of the geological structure of the land of the KRB and the diversity of its minerals are due to its location at the junction of the two largest geosynclinal-folded mobile zones, the Ural-Mongolian and Mediterranean, which respectively include the Caledonian-Hercynian North and Mid-Tien Shan and the Herzino-Cimmerian North and South Pamir geosynclinal fold belts. The KRB is based on the geological developments related to Central (Southern Tien Shan and Gissaro-Alai) and South-Western Tajikistan. Central Tajikistan (Gissaro-Alai) is part of the Hercynian geosynclinals, the folded region of South Tien Shan. The geological structure of the latter part mainly includes Middle and Upper Paleozoic and, to a lesser extent, Precambrian, Lower Paleozoic, Mesozoic, and Cenozoic deposits. South-Western Tajikistan is represented by the Tajik depression, the Epigercin Alpine sub-platform, and the eastern flank of the Turan Plate [58].
Tajikistan is included at the center of the origin of cultivated plants, in the Southwest Asia and Central Asian foci. Its flora now includes 5000 species. It is here that the homelands of the most important grains and leguminous crops, such as soft wheat, peas, lentils, and chickpeas, were established. It is recognized that the foothills and mountain zones of Tajikistan are centers of the origin of fruit crops, such as Juglans Regia, Zizyphus Sativa, Pistacia vera L., Pyrus communis L., Amygdalus communis L., Elaeagnus Angustifolia, and Vitis vinifera L. Tajikistan is also notable as the secondary center of origin for a variety of fruit crops: Armeniaca Vulgaris, Malus sieversii, Vitis vinifera, Morus alba L., and Punica granatum. Significant in this region are the diversity of wild species and forms of Berberis L., Prunus Sogdiana, Crataegus L., Eleagnus L., and Hippophae rhamnoides L [59].
The inter-annual variation of temperature and precipitation for all six meteorological stations during the 1951-2012 time period is shown in Figure 2. The seasonal temperature and precipitation variability were characterized in winter (December, January, and February), spring (March, April, and May), summer (June, July, and August), and fall (September, October, and November). The coldest and hottest months for all meteorological stations were January and July, with a mean minimum and mean maximum temperature from −12.1 • C to +9.7 • C at Anzob, from −17 • C to +24.2 • C at Hushyori, from +0.4 • C to +25.5 • C at Faizobod, from +2.3 • C to +27.1 • C at Dushanbe, from +1.9 • C to +30.4 • C at Isambay, and from +1.8 • C to +31.4 • C at Shaartuz. The absolute maximum and minimum temperature at all climate stations ranged from +24 • C to −10 • C at Anzob, from +39 • C to +3 • C at Hushyori, from +41 • C to −26 • C at Faizobod, from −27 • C to +46 • C at Dushanbe, from +47 • C to −27 • C at Isambay, and from +47 • C to −22 •  The inter-annual variation of temperature and precipitation for all six meteorological stations during the 1951-2012 time period is shown in Figure 2. The seasonal temperature and precipitation variability were characterized in winter (December, January, and February), spring (March, April, and May), summer (June, July, and August), and fall (September, October, and November). The coldest and hottest months for all meteorological stations were January and July, with a mean minimum and mean maximum temperature from −12.1 °C to +9. 7   The global temperature trends according to the IPCC report present the temperature per decade; the CRUTEM4.1.1.0 dataset was explored by Jones et al. [9], GHCNv3.2.0 by Lawrimore et al. [10], and GISS by Hansen et al. [11] and Rohde et al. [12] during the five historical time periods, as shown in Figure 3 and Table 2. At a global scale (after IPCC, 2013) [60], the precipitation trends per The global temperature trends according to the IPCC report present the temperature per decade; the CRUTEM4.1.1.0 dataset was explored by Jones et al. [9], GHCNv3.2.0 by Lawrimore et al. [10], and GISS by Hansen et al. [11] and Rohde et al. [12] during the five historical time periods, as shown in Figure 3 and Table 2. At a global scale (after IPCC, 2013) [60], the precipitation trends per decade were investigated by applying various datasets; the CRU TS 3.10.01 was used by Mitchell and Jones [21], GHCN V2 by Vose [22], GPCC V6 by Becker [23] and Smith et al. [24] in two past time horizons, as shown in Figure 4. In addition, according to the IPCC 2013 [60], Dai et al. [31] studied the record data for 925 gauging station across the globe on the largest rivers. Dai et al. [31] reported that the discharges in the 200 rivers showed statistically significant trends in the period of 1948-2004, with a downward trend in 45 rivers and upward trends in 19 rivers. Stahl et al. [61] analyzed the streamflow datasets of more than 400 small near-natural catchments in Europe to investigate the river discharge trends for the period of 1962-2004. The analyses of Stahl et al. showed negative streamflow trends in the southern and eastern regions and positive trends elsewhere. The temporal trends of natural river discharge varied with the period of study. The comprehensive analyses during the 20th century IPCC Report (2013) led to low confidence in an increasing trend of global river discharge [60]. In addition, in this study, we present streamflow trend analyses to improve our understanding of the hydrological regime in the upstream area of the ungauged Kofarnihon River Basin against the background of the variability of the river regime of other rivers in other regions across the globe. that the discharges in the 200 rivers showed statistically significant trends in the period of 1948-2004, with a downward trend in 45 rivers and upward trends in 19 rivers. Stahl et al. [60] analyzed the streamflow datasets of more than 400 small near-natural catchments in Europe to investigate the river discharge trends for the period of 1962-2004. The analyses of Stahl et al. showed negative streamflow trends in the southern and eastern regions and positive trends elsewhere. The temporal trends of natural river discharge varied with the period of study. The comprehensive analyses during the 20th century IPCC Report (2013) led to low confidence in an increasing trend of global river discharge [60]. In addition, in this study, we present streamflow trend analyses to improve our understanding of the hydrological regime in the upstream area of the ungauged Kofarnihon River Basin against the background of the variability of the river regime of other rivers in other regions across the globe.    The significance levels tested are *** for p < 0.001, ** for p < 0.01, * for p < 0.05, NS for p ≥ 0.1.

Methodology
Globally, statistical and gridding methods have been widely used to estimate the long-term temperature, precipitation, and streamflow trends [62][63][64]. In this study, based on measured station data, the statistical methods of Mann-Kendall (MK) test, modified Mann-Kendall (MMK) test, and Sen's slope were applied to determine the trends in a long-term time series of temperature, precipitation, and river flow. To determine abrupt changes, the Pettitt's test was also used.

Mann-Kendall (MK) Test
The MK test is a non-parametric test for identifying trends in time series data [65]. Using linear regression, the trend rate m 1 is determined applying Equation (1): where x indicates the temperature and precipitation for the periods of 1951-2012 and 1951-2008. Using a t-test, the significance of m 1 was verified. The positive and negative values of m 1 show increasing and decreasing trends of temperature, precipitation, and flow in particular time series [66]. Yue and Wang [67] reported that the strength of the trend depends on the magnitude of the trend, the sample size, the adjusted significance level, and the number of alterations over a time series. The MK test statistics (S) are determined by Equation (2): where n is the length of the dataset, X j and X k are the data values at times j and k, and sgn is the sign function that takes on values of −1, 0, and +1. The resulting t value of S shows upward or downward trends in the climatic datasets: where p is the tied group, and t k is the number of observations in the k th group. As determined by Equation (5), the (STS) standardized test statistic (Z S ) can be determined by where Z S shows the significance of the trend. The STS is then applied to test the null hypothesis, H 0 if Z s > Z α/2 , while α shows the confidence level. In this study, if a trend in the temperature, precipitation, or flow time series is statistically significant at a significance level of α = 0.01 (or 99% confidence intervals), α = 0.05 (or 95% confidence intervals), and α = 0.

Sen's Slope Test
Using Sen's slope test, the magnitude of the trends over the time series was calculated. The slope of the pair of data "n" was first calculated using Equation (6): where X j and X k are the data values at times j and k. The median of the "n" values of Q is Sen's estimator of slope. The positive and negative values of Q present an increasing and a decreasing trend in the temperature, precipitation, and flow datasets. The slopes (Q) of the "n" values were ranked from low to high. Sen's estimator can be calculated using Equation (7): if n is even.

Pettitt's Test
To detect a single change point in temperature, precipitation and flow time series with continuous data the Pettitt [52] approach is applied. Usually, when an adjacent change point exists in a series, the maximum value K T will be determined as the change point: If U t,T increases with time t, then the sequence does not have a change point; conversely, if U t,T shows a decreasing trend against time t, then a change point occurred in the series. The change point of the series is located at K T , provided that the statistic is significant. The significance probability of K T is approximated for p ≤ 0.05 with

Results
Based on the long-term measured station data in the ungauged Kofarnihon River Basin, the trends and magnitudinal changes of monthly and annual temperature, precipitation, and streamflow were determined by applying the non-parametric Mann-Kendall, modified Mann-Kendall, and Sen's slope tests. We also used Pettitt's test to analyze the changes during the 1951-2008, 1951-2012, and 1979-2012 time periods. In this study, the availability of observational hydro-climatic measured data in the KRB and its utilization in regional analysis were essential to derive comprehensive outcomes of the hydro-climatic variability in recent decades. In addition, we presented the temperature, precipitation, and streamflow trend analyses from the IPCC report 2013 solely for demonstration purposes, as the period of our study in the KRB was the same as that in the IPCC report 2013 (Tables 3-6). In continental and global scale trend analyses most studies were performed based on gridded data, which has more uncertainty than observational data. The global trend analysis according to the IPCC report 2013 was described based on the datasets of CRUTEM, GHCNv3, GISS, and Berkeley for the temperature trend; CRUTS, GHCN V2, GPCC V6, and Smith et al. for the precipitation trends; and Dai et al. [31] and Stahl et al. [61] for the streamflow trends. Goddard Institute of Space Studies (GISS) [11] 0.188 ± 0.032 0.267 ± 0.054 Berkeley [12] 0.175 ± 0.029 0.254 ± 0.049 The significance levels tested are *** for p < 0.001, ** for p < 0.01, * for p < 0.05, NS for p ≥ 0.1.

Trend Analysis of the Seasonal Temperature in the Kofarnihon River Basin
The seasonal temperature trends at the six climate stations during the 1951-2012 and 1979-2012 periods in the Kofarnihon River Basin are shown in Table 2. The results of the original MK test showed a significant increasing trend for the mean monthly temperature during the 1951-2012 period in spring (March, April, and May), summer (June, July, and August), and fall (September, October, and November) at the five climate stations, except for Anzob station, where a decreasing trend was observed in winter (December, January, and February). In the low-altitude areas, the change in the mean monthly temperature is somewhat complicated, particularly at the Dushanbe station. At this station, the mean monthly temperature showed an increasing trend during the 1951-2012 period in all months. The 1979-2012 period presented a decreasing trend over nine months, except for February, March, and April. The original MK test of the seasonal mean temperature for the periods of 1951-2012 and 1979-2012 demonstrated that in the Kofarnihon River Basin, the summer temperature trend significantly increased at high-altitudes and that the winter temperature trend significantly decreased during the 1951-2012 period. In the low-altitude areas, the spring and fall temperature trends increased, and the winter temperature decreased at all six stations during the 1979-2012 period. The findings of the winter cooling in our basin are in agreement with Li et al. [68] and Wang et al. [69]. The Eurasian winter cooling trend during the 1998-2012 period is caused internal atmospheric variability and constitutes an extreme climate event. Furthermore, the observed reduction in Artic sea ice enhances the probability of Eurasian winter and thus increases the probability of an extreme Eurasian winter cooling trend [68]. East Asia experienced more cold winters and significantly negative surface air temperature anomalies during the recent strong East Asian winter monsoon epoch spanning the period 2004-2012 [69]. Cooling has generally coincided with negative winter temperature trends in eastern Eurasia over the last two decades, possibly contributing to the observed regional cooling. Enhanced wintertime blocking activities around the Ural mountain region and the diminished Arctic sea ice concentrations are suggested to be the internal atmospheric processes responsible for this phenomenon and the external drivers for the recent re-amplification of the East Asian winter monsoon [69]. The significance levels tested are ** for p < 0.01, * for p < 0.05.

Trend Analysis of the Annual and Decadal Temperature
We computed the temperature trends per decade during the 1951-2012 and 1979-2012 time periods, which correspond to the periods of the global temperature trends based on the trend per decade (T/D), as shown in Table 3.
The results of the annual temperature trend per decade, the Sen's slope, the abrupt change detection, and the original and modified Mann-Kendall trend tests are shown in Table 4. The mean annual temperature trend showed an increasing trend during the 1951-2012 period. The highest temperature trends of 0.25 • C/decade and 0.23 • C/decade were observed at the Dushanbe and Shaartuz stations; the lowest of 0.01 • C/decade was found at Anzob, with 0.07 • C/decade at Hushyori. According to the original and modified MK tests, the mean annual temperature showed an increasing trend during the 1951-2012 period, which is in line with the linear trend of the global land surface air temperature (Table 3). Since the 1930s, the mean annual temperatures have increased, and the average rate of variation has remained around 0.1 • C per decade in Tajikistan [20]. For the five climate stations under study, the mean annual temperature showed a significant increasing trend during the 1951-2012 period, except for the Anzob station, where a non-significant trend was found. Anzob is located at the highest altitude compared to the other reference stations in the catchment.
The temperature trend per decade during the 1979-2012 period at all stations ranged from 0.36 • C/decade to −0.18 • C/decade ( Figure 5). Our results for the temperature trend per decade mostly correspond with the linear trends in global air surface temperature. There is a very small mismatch in the value ranges between our results in Table 4 and the datasets in Table 2. This small inconsistency in the ranges of the values can be ascribed to the use of observational data in the current study, which provide more accurate outcomes at the basin scale. Furthermore, the result of the original and modified MK tests showed a significant increase in the mean annual temperature trends at the Hushyori, Isambay, and Shaartuz stations, while a significant decreasing trend was found at the Anzob station during the 1979-2012 period. This discrepancy in the temperature trend could be attributed to the topography of the basin. The results of the Pettitt's test showed that a change in temperature from low to high-altitudes occurred in various periods. The Pettit test results showed that an abrupt change in the high-altitude areas occurred later than that in the middle-and low-altitude areas of the Kofarnihon River Basin.

Trend Analysis of the Seasonal Precipitation in the Kofarnihon River Basin
The original MK test result of the seasonal precipitation trend during the 1951-2008 period in the Kofarnihon River Basin in Central Asia is shown in Table 5. The seasonal precipitation showed decreasing trends at all six climate stations in spring (March, April, and May). In winter (December, January, and February); however, we found an increasing trend for seasonal precipitation during the 1951-2008 period at most climate stations. Our results for seasonal precipitation reveal that in the high-and low-altitude areas, the spring precipitation significantly decreased; the winter precipitation significantly increased only in the high-altitude areas in the Kofarnihon River Basin in Central Asia. Table 6 presents a global analysis of the linear trend in the annual precipitation of the northern hemispheric and annual global land precipitation during the 1951-2008 period. In this study, the period of the annual and decadal precipitation trends was also chosen 1951-2008.

Trend Analysis of the Annual and Decadal Precipitation
The results of the modified and original MK tests showed a decreasing trend in the annual and decadal precipitation during the 1951-2008 period. Among the six climate stations, the decreasing trend of the decadal precipitation (P/D) ranged from −14.63 mm year −1 per decade to −4.76 mm year −1 per decade ( Figure 6). The negative P/D trend in the Kofarnihon River Basin was stronger (Table 7) compared to the indications of the global analyses of the northern hemisphere and the global land areas ( Table 6). The modified MK test presented a significant decreasing trend in the annual precipitation solely in the low-altitude areas, such as the Isambay and Shaartuz stations, while a non-significant trend was found in the high-altitude areas. The results of the Pettitt's test showed the change point to occur in 1994 at the Anzob station during the 1951-2008 time period. During this period, the precipitation changes at the Faizobod and Shaartuz stations occurred in 1999, while at the Hushtori, Dushanbe, and Isambay stations, the precipitation change occurred in 1969. The Pettitt's test results demonstrated that in this mountainous area, the earlier changes occurred detected in the low-altitude areas. This change could be attributable to the earlier effects of climate change on the low-altitude areas of the mountainous Kofarnihon River Basin in Central Asia.
correspond with the linear trends in global air surface temperature. There is a very small mismatch in the value ranges between our results in Table 4 and the datasets in Table 2. This small inconsistency in the ranges of the values can be ascribed to the use of observational data in the current study, which provide more accurate outcomes at the basin scale. Furthermore, the result of the original and modified MK tests showed a significant increase in the mean annual temperature trends at the Hushyori, Isambay, and Shaartuz stations, while a significant decreasing trend was found at the Anzob station during the 1979-2012 period. This discrepancy in the temperature trend could be attributed to the topography of the basin. The results of the Pettitt's test showed that a change in temperature from low to high-altitudes occurred in various periods. During the 1951-2012 time period, the change point at most stations appeared after 1976 and during the 1979-2012 period (after the 1990s), except for the Anzob station, where the change point occurred in 2007 during the 1951-2012 and 1979-2012 periods. The Pettit test results showed that an abrupt change in the high-altitude areas occurred later than that in the middle-and low-altitude areas of the Kofarnihon River Basin.   The significance levels tested are *** for p < 0.001, ** for p < 0.01, NS for p ≥ 0.1.  Based on the data from the six measured stations, the annual mean precipitation and annual mean air temperature were analyzed using Sen's slope (Figure 7). The mean air temperature increased from the upstream areas to the downstream areas of the KRB. The increased rate of air temperatures around the Dushanbe station was the highest in the KRB. This station is located in the capital city of Tajikistan, where urban built-up land is the main land use type near the station. The increased rate of annual precipitation in the KRB increased gradually from upstream areas to downstream areas. The annual precipitation decreased, and the annual mean air temperature increased, throughout the KRB.

The Trend Analysis of the Seasonal Streamflow in the Kofarnihon River Basin
The seasonal streamflow trend analysis was performed in the Varzob River based on the Dahana discharge station data for the 1951-2012 and 1979-2012 periods, as shown in Table 8. The Varzob River is an upstream and northwest tributary of the Kofarnihon River Basin. The results of the original MK test showed an increasing trend in the mean monthly streamflow in winter, spring, and summer, while the summer streamflow presented a significant increasing trend during the 1951-2012 period. The summer and spring streamflow showed an increasing trend, and the fall and winter streamflow indicated a decreasing trend during the 1979-2012 period. Our results show that in spring, the mean monthly streamflow significantly increased, and the fall streamflow significantly decreased, during the 1951-2012 and 1979-2012 periods in the upstream area of the KRB.
temperatures around the Dushanbe station was the highest in the KRB. This station is located in the capital city of Tajikistan, where urban built-up land is the main land use type near the station. The increased rate of annual precipitation in the KRB increased gradually from upstream areas to downstream areas. The annual precipitation decreased, and the annual mean air temperature increased, throughout the KRB.

The Trend Analysis of the Seasonal Streamflow in the Kofarnihon River Basin
The seasonal streamflow trend analysis was performed in the Varzob River based on the Dahana discharge station data for the 1951-2012 and 1979-2012 periods, as shown in Table 8. The Varzob River is an upstream and northwest tributary of the Kofarnihon River Basin. The results of the original MK test showed an increasing trend in the mean monthly streamflow in winter, spring, and summer, while the summer streamflow presented a significant increasing trend during the 1951-2012 period. The summer and spring streamflow showed an increasing trend, and the fall and winter streamflow indicated a decreasing trend during the 1979-2012 period. Our results show that in spring, the mean monthly streamflow significantly increased, and the fall streamflow significantly decreased, during the 1951-2012 and 1979-2012 periods in the upstream area of the KRB.  The significance levels tested are ** for p < 0.01, * for p < 0.05.

Trend Analysis of the Annual Streamflow in the Kofarnihon River Basin
The modified MK test showed a significant increasing trend for the annual streamflow in the upstream areas of the KRB during the 1951-2012 and 1979-2012 time periods. Increased flow upstream can be ascribed to the rapidly melting snow and ice caused by the continuously increasing air temperature. The large mass of water formed in the Varzob River experiences little demand during the growing season because the mountainous terrain in which the river flows does not have fertile valley plots of land. The water formed from the Varzob River almost completely flows downstream of the Kofarnihon River. In the downstream area of the KRB, the annual streamflow showed a decreasing trend according to Lobanova et al. This discrepancy indicates that, over time and space, anthropogenic activities have strictly disturbed the water distribution in the middle and downstream areas of the KRB. The decreasing trend of the annual streamflow in the downstream areas of the KRB is attributed to the increased water use with rapid population growth and the degradation of land under irrigation. Most of the populated cities are located in the middle and downstream areas, such as the capital city of Tajikistan-Dushanbe and the city of Vahdat, along with several well-populated districts in the territory of the KRB. The large-scale water intake and water consumption from the river occurs in the middle and downstream portions of the river. In the upstream areas of the KRB, the streamflow change point according to the Pettitt's test result occurred in 1991 during the 1951-2012 and 1979-2012 periods, as shown in Table 9 and Figure 8. The long-term observation data in the KR showed that the annual runoff of the river is 30-50 l/s per km 2 . The flooding period in this river starts at the end of March and the beginning of April, while the duration of this period depends on the snow cover and the water reserves in the snow. Peak floods occur in the May-June months. The reserves of the seasonal snow cover deplete in the month of July and begin to melt many small glaciers, depending on the warm air temperature in the The long-term observation data in the KR showed that the annual runoff of the river is 30-50 l/s per km 2 . The flooding period in this river starts at the end of March and the beginning of April, while the duration of this period depends on the snow cover and the water reserves in the snow. Peak floods occur in the May-June months. The reserves of the seasonal snow cover deplete in the month of July and begin to melt many small glaciers, depending on the warm air temperature in the summer season. The duration of the flooding is 160-180 days, and the increase in flooding waves lasts 60-80 days.
The glaciers and snow cover significantly affect the water cycle and flow process of the KR; these processes can be significantly changed by climatic variation. Using field observations, we observed that the territory of the Anzob climate station is covered by snow for a period of eight months each year. The long-term data analysis of the snow-fed and glacier-fed river showed a constantly increasing trend of the streamflow in the upstream tributaries of the KRB during the 1951-2012 and 1979-2012 time periods. The result of the Pettitt's test showed that abrupt changes in streamflow occurred in 1991 in the upstream areas of the KRB during the 1951-2012 and 1979-2012 periods. We conclude that the increasing streamflow in the upstream regions of the KRB may be attributable to an increase in the snow and glaciers melting in the basin, which is influenced by the continuously increasing air temperature in the region. Studying the recent changes in the water cycle and flow regime of the Kofarnihon River is essential for Tajikistan, which uses water from the KR for drinking, irrigation, and generating electricity. It is also important for downstream countries, such as Uzbekistan and Turkmenistan, to adapt their national agricultural strategies to anomalous hydro-climatic conditions.

Discussion
The Kofarnihon River is the one of the northwest tributaries of the Amu Darya River. The climate of the basin is continental under the influence of westerly wind, leading to significant seasonal changes in temperature and precipitation [20]. These changes are due to the mountainous topography of the basin, which has a very high local contrast [1,70]. The mountainous area of this basin is the most vulnerable territory in Central Asia to hydro-climatic changes [2]. The tributaries of the Amu Darya River, such as the Panj River, Vakhsh River, and Kofarnihon River, are characterized by a seasonal cycle of river flow; the streamflow increases rapidly with the beginning of the melting period in March and April, the maximum river flow occurring in the summer season [71]. The present study examined the recent trends of the decadal, annual, and seasonal temperature, precipitation, and streamflow during the 1951-2012 and 1979-2012 time periods in the Kofarnihon River Basin in Central Asia.
Our results for the MK test revealed that the summer temperature trend significantly increased, and the winter temperature trend significantly decreased, during the 1951-2012 period in the high-altitude regions of the KRB. In the low-altitude areas, the spring and fall temperature trends increased, and the winter temperature decreased, at all six stations during the 1979-2012 period. Our results for the seasonal temperature trends also support the findings of Mannig et al. [17], who reported that the climate change signal of temperature is the largest in winter in the mountainous terrains of Central Asia. Piao et al. [33] reported that a reduction in the warming trend mostly occurred in the winter in the northern hemisphere. Similar results were presented by Pachauri et al. [72] for the global winter temperature trends during recent decades. Further, Li et al. [68] and Wang et al. [69] confirmed that in the Eastern Asia, the winter temperature decreased during the 1998-2012 and 2004-2012 periods. The decreasing probability of an extreme Eurasian winter temperature trend could be attributed to a reduction in Arctic sea ice enhancing the probability of Eurasian winters [68]. In this study, the mean annual temperature trend mostly showed an increasing trend during the 1951-2012 and 1979-2012 periods. These findings are in line with the results Chevallier et al. [19], who confirmed that the overall homogeneous trend of temperature increased in the mountainous region of Tajikistan. Aalto et al. [20] noted that since the 1930s, the mean annual temperatures have increased in Tajikistan. Giese et al. [14] and Christensen et al. [15] demonstrated that Central Asia has faced above-average global warming since the 1950s; the mean temperature increased by 1.2-2.1 • C, and during the 20th century, it increased by 1-2 • C. Our findings also agree with the global land surface air temperature trends shown in Tables 3 and 4 [9][10][11][12]. The climate changes in Tajikistan are mostly limited by a variety of topographical situations due to adiabatic processes [73]. The temperature changes are likely influenced by other factors, such as solar radiation and cold-air drainage in mountainous terrains [74]. Another possible cause of the rising temperatures in the Kofarnihon River Basin could be the increasing concentration of aerosols and greenhouse gases in the regional atmosphere. Xin et al. [75] reported that over China and Central Asia, a rising trend in the concentration of aerosols in the atmosphere could cause a substantial increase in temperature. Bollasina et al. [76] confirmed that in Asia, the concentration of atmospheric aerosols has increased steadily. Similarly, Folini et al. [77] reported that aerosol emissions in the 20th century increased in association with the enormous growth in population and industrialization.
Seasonal precipitation showed a decreasing tend in the spring season in both high-and low-altitude regions, while the winter precipitation presented a significant increasing trend in the high-altitude regions during the 1951-2008 period in the Kofarnihon River Basin in Central Asia. Similarly, a significant increasing trend of the winter precipitation was found by Chen et al. [28] and Song et al. [30] for the Central Asia region, where they defined an increasing trend during the 1930-2009, 1960-1991, and 1960-2013 time periods. Pendergrass et al. [78] reported that the global winter precipitation increased during the second half of the 20th century and attributed this increase to the role of increasing moisture counteracted by weakening circulation. Our results for the annual precipitation trend showed a decreasing trend over the catchment at all stations under study. This result for the annual precipitation follows the same statistical trends observed over Central Asia [27]. Li et al. [27] noted that in Central Asia, at the end of 20th century, there was a persistent decreasing trend in annual precipitation. Chen et al. [28] studied annual precipitation using the monthly gridded method of precipitation from the Climatic Research Unit (CRU) during the 1930-2009 period and concluded that, generally, precipitation increased over the past 80 years in the westerly circulation-dominated regions of arid Central Asia. According to Ahmed et al. [29], the modified Mann-Kendall trend test results showed a significant increasing trend in annual precipitation during the 1961-2015 period in the Yangtze River of the Qinghai-Tibetan Plateau. In the Asian region, these contradictions in our findings may be associated with the increasing concentration of anthropogenic absorbing aerosols and the westerly system [76,77,79]. Donat et al. reported that in the dry regions of the globe, precipitation strongly increased, while in wet regions, only small changes were observed over the period of 1951-2010 [25]. We revealed that the decreasing trend of the annual precipitation in the Kofarnihon River Basin in Central Asia was more significant during the 1951-2008 period compared to the same period in previous studies on the northern hemisphere and at global scales, where different datasets were applied by Mitchell et al. [21], Vose et al. [22], Becker et al. [23], and Smith et al. [24] (Tables 6 and 7).
The Varzob River is one of the main tributaries and upstream of the Kofarnihon River. In this study, a trend analysis was performed based on the long-term streamflow data from the Dahana discharge station of the Varzob River. The seasonal streamflow trend in the upstream area of the KRB showed an increasing trend in spring and summer but a decreasing trend in fall and winter. In spring, the streamflow showed a significant increasing trend during the 1951-2012 period. Our results for the streamflow trend upstream of the KR contradicts the trend of the streamflow in the downstream areas of the KR, which was studied by Lobanova et al., who reported that the flow decreases in winter and fall; however, in spring, the flow increases and does not change during the summer season in the downstream portion of the KR [38]. Highly-populated cities of Tajikistan are located in the middle stream and downstream areas of the KRB, where the drinking and irrigating land use of water from the KR decreased the water availability downstream. Previous studies on the Yangtze River in Asia by Piao et al. [33] and Pachauri et al. [72] described an increasing trend of the streamflow in summer and a decreasing trend in fall. Jaweso et al. [80] reported that the monthly MK test exhibited a significant decreasing trend of the streamflow in the dry season and a significant increasing trend in the wet season during the 1980-2008 period in the Upper Omo-Ghibe River Basin in Ethiopia. The results of the modified MK test showed a significant increasing trend of the annual streamflow during the 1951-2012 and 1979-2012 periods in the upstream areas of the Kofarnihon River. Our results of the annual streamflow trends agree with the results of previous studies that addressed the hydrological regime variations of the main rivers in Central Asia [34,38,81]. Shen et al. [35] noted that the MK test trend showed a significant increasing trend in the annual streamflow in four glacierized basins (Toxkon, Kunmalik, Kaidu and Huangshuigou) in the southern Tian Shan Mountains in Central Asia, especially during the spring and winter months. Zhang et al. [82] confirmed that the MK test results of the annual streamflow between the 1970s and the 1990s indicate a decreasing trend in the Yarlung Zangbo River in the Qinghai-Tibetan Plateau in Central Asia. Further, the MK trend test exhibited a significant decreasing trend of the annual streamflow in the rain-fed Yihe River in the central part of China during 1960-2008 [61]. Piao et al. reported a decreasing trend of the streamflow in the Yellow River in Asia since 1960s [33]. The rain-fed river basins showed that the streamflow increased in the Lanyang River Basin, Youluo River Basin, and Shangping River Basin. However, the streamflow decreased in the Keelung, Dahan, and Fengshan river basins during the 1980-2017 period in the East Asian regions [83]. Dai et al. [31] investigated the increasing and decreasing trends of the streamflow in different rivers across the globe and reported that in 45 rivers, the streamflow trend increased, while in 19 rivers, the streamflow decreased. The IPCC report 2013 noted that confidence is low for an increasing trend in global river discharge during the 20th century [60]. However, Rashid et al. [37] confirmed that the glacier-fed and snow-fed Lidder and Sindh rivers in the Himalayan Mountains exhibited a statistically significant depletion in their streamflow during the 1962-2018 period. According to Hooke [6] and Schaefli et al. [4], an increasing flow in mountainous rivers can be mostly attributed to global warming, which causes severe snowmelt and the retreat of glaciers. The Kofarnihon River in this study is a snow-fed and glacier-fed river; thus, these increases in summer water flow can be ascribed to the influence of global warming on the rapid melting of snow and ice.
The trends of hydro-climatic variables, such as temperature, precipitation, and streamflow, were studied in the mountainous KRB of Central Asia. We conclude that the mean annual temperature showed an increasing trend, while the annual precipitation showed a decreasing trend. The annual mean streamflow trends upstream showed an increasing trend during the 1951-2012 and 1979-2012 periods. In addition, this study revealed that in the mountainous region of the KR, the streamflow mostly depends on snowmelt and glacial melt, not rainfall. This indicates that climate warming in this mountainous area accelerated the melting of snow and glaciers, which caused an increase in river discharge. The use of hydro-climatic trend analysis based on observational data at the basin scale is essential for local agronomists, hydropower engineers, and decision makers to develop new strategies for anomalous hydro-climatic conditions. Further investigations could focus on combining the trend analysis of observational data with global climate model (GCM) data coupled with a hydrological model to predict the impacts of climate change and anthropogenic activities on water resources in the KRB.

Conclusions
In this study, the trends, abrupt changes, and periods of annual and seasonal temperature, precipitation, and streamflow series were analyzed. We applied non-parametric Mann-Kendall, modified Mann-Kendall and Sen's slope tests to monthly and annual hydro-climatic variables to determine the trends and magnitudinal changes in the Kofarnihon River Basin in Central Asia. We also used Pettitt's test to analyze the changes during the 1951-2012 and 1979-2012 time periods. The following conclusions can be summarized from the above analysis results.
(1) The mean monthly temperature showed a significant increasing trend in spring, summer, and fall at the five climate stations, except for Anzob station, where we found a decreasing trend in winter during the 1951-2012 period. The seasonal mean temperature indicated that in the high-altitude areas of the KRB, the summer temperature trend significantly increased, and the winter temperature trend decreased. The mean annual temperature showed an increasing trend, with the highest trends of 0.25 • C/decade and 0.23 • C/decade at the Dushanbe and Shaartuz stations, respectively; the lowest of 0.01 • C/decade and 0.07 • C/decade was found at the Anzob and Hushyori stations. The change points mainly occurred in the 1970s and 1990s for all stations, except for Anzob station, where the change point appeared in 2007. The Pettitt's test results showed that the abrupt changes in the high-altitude areas occurred later than those in the middleand low-altitude areas in the KRB. This study analyzed only the hydro-climatic spatiotemporal variability and did not consider the impact of anthropogenic activities. In the future, relevant indicators should be selected for further study, and the impact of anthropogenic activities should be analyzed. We revealed that the hydro-climatic variability in the upstream areas of the KRB impacts downstream water fluctuations. Our findings in this study will provide a scientific theoretical basis for water issues under ongoing climate change and help us to improve the water resources management in the KRB of Central Asia.