Trend Analysis of Annual and Seasonal River Runoff by Using Innovative Trend Analysis with Signiﬁcant Test

: This study investigated the temporal patterns of annual and seasonal river runoff data at 13 hydrological stations in the Lake Issyk-Kul basin, Central Asia. The temporal trends were analyzed using the innovative trend analysis (ITA) method with signiﬁcance testing. The ITA method results were compared with the Mann-Kendall (MK) trend test at a 95% conﬁdence level. The comparison results revealed that the ITA method could effectively identify the trends detected by the MK trend test. Speciﬁcally, the MK test found that the time series percentage decreased from 46.15% in the north to 25.64% in the south, while the ITA method revealed a similar rate of decrease, from 39.2% to 29.4%. According to the temporal distribution of the MK test, signiﬁcantly increasing (decreasing) trends were observed in 5 (0), 6 (2), 4 (3), 8 (0), and 8 (1) time series in annual, spring, summer, autumn, and winter river runoff data. At the same time, the ITA method detected signiﬁcant trends in 7 (1), 9 (3), 6(3), 9 (3), and 8 (2) time series in the study area. As for the ITA method, the “peak” values of 24 time series (26.97%) exhibited increasing patterns, 25 time series (28.09%) displayed increasing patterns for “low” values, and 40 time series (44.94%) showed increasing patterns for “medium” values. According to the “low”, “medium”, and “peak” values, ﬁve time series (33.33%), seven time series (46.67%), and three time series (20%) manifested decreasing trends, respectively. These results detailed the patterns of annual and seasonal river runoff data series by evaluating “low”, “medium”, and “peak” values.


Introduction
Climate change affects regional water resources [1]. Water resources are an essential foundation for society's sustainable development, the economy, and the ecological environment [2][3][4]. As a unique geographic component, the mountain area is the birthplace of most rivers in arid regions. Almost all rivers originate in the mountainous regions [5]. Central Asia's water resources are mainly recharged by melting ice and snow in the mountainous areas, making them more sensitive to climate change [6][7][8]. Human activities have led to global environmental changes [9]. Due to population growth and disproportional spatial distribution, developing industries, agricultural expansion, and increased urban construction, there is increased pressure on water and other natural resources because the anthropogenic demand for water is growing [10]. The river runoff is one of the primary forms and components of water resources and the most critical water source for production and living [11,12]. Affected by climate variability and anthropogenic activities, the river used in trend calculations and is better than traditional methods, even when relying on time series.
As one of the crucial links in the water cycle, river runoff is affected by both climate (precipitation, temperature, evaporation) and human activities (projects involving the underlying surface, topography, and water conservancy, among others), and is highly complex, dynamic, and nonlinear [60,61]. Lake Issyk-Kul is an endorheic mountain lake [62]. Accordingly, climate change and human activities are affecting the lake's surface runoff [63]. Many researchers have pointed out that surface runoff changes are the main driving factor in the water level fluctuations of Lake Issyk-Kul [64][65][66]. Exploring river runoff variability of in this endorheic lake basin is essential for developing sustainable water resources in the study area [67]. Still, there have been few studies on the characteristics of river runoff changes. To the best of our knowledge, there is no published research on the use of ITA of river runoff in the Lake Issyk-Kul Basin. Moreover, there are some deficiencies in the existing studies: (1) Most previous research focused on annual river runoff changes, while only a few investigations performed detailed analyses on monthly runoff changes; (2) The collected river runoff data durations have not been long, and the number of hydrological stations used in the studies has been limited; and (3) The "low", "medium", and "peak" trends of river runoff were not explored. Keeping the above reviews in mind, this study was conducted with the following specific objectives: (i) to analyze the temporal trend in seasonal and annual river runoff time series data using the MK test, Theil-Sen Median (TSM), and new ITA method with significance test; (ii) to compare the results of MK method with the new ITA.

Study Area
Lake Issyk-Kul ( Figure 1) is an inland saltwater lake located in the basins of the Teskey Ala-Too Range and Kungey Ala-Too mountain ranges in the Western Tian Shan Mountains in northeastern Kyrgyzstan [68]. The center of the lake is 42.42 • N, 77.33 • E. Lake Issyk-Kul is 60.1 km from north to south and 178 km from east to west, with an area of 6236 km 2 . Its water storage capacity can reach 1730×10 8 m 3 , the perimeter of the lake shoreline is approximately 669 km, the average lake depth is 278.4 m, and the maximum lake depth is 702 m [64,69]. Among the world's mountain lakes, Lake Issyk-Kul is second only to Lake Titicaca in South America in terms of area and ranks first of lake depth and lake water volume. The climate in the lake area is mild and dry. The average temperature in July is 15-25 • C, and the average temperature in January is −6 • C. The Lake Issyk-Kul drainage basin occupies a surface area of 22,080 km 2 , and there are 118 large and small rivers with a surface runoff of about 36.7 × 10 8 m 3 [62]. The main river flows are Chong-Koi-Suu, Cholpon-Ata, Chong-Ak-Suu, Ak-Suu, Chong-Uryukty, Pzhergalan, Ak-sai, Tossor, Ton, Dzhuuku, Tamga, Karakol, and Chong-Kyzyl-Su [63].
Lake Issyk-Kul is a closed inland lake with no outlet. Many rivers are used for irrigation before flowing into the lake, but rivers are not the only source of lake water. The annual surface water supply is 13 × 10 8 m 3 , and the annual total of precipitation and groundwater is 330 mm. The annual evaporation of the lake is 820 mm.

Data Sources
This paper obtained the monthly and annual time series data (river runoff, covering the 1943-2012 period) selected from 13 hydrological stations surrounding Lake Issyk-Kul, as shown in Figure 1. This study compiled observational data provided by the National Science & Technology Infrastructure Center-Data Sharing Infrastructure of Earth System Science and the Department of Water Resources and Irrigation of the Ministry of Agriculture and Land Reclamation of the Kyrgyz Republic and the National Statistical Committee of the Kyrgyz Republic [63]. In this study, data from the four prominent seasons-spring (March-April-May), summer (June-July-August), autumn (September-October-November), and winter (December-January-February)-were analyzed. Also, the trends of seasonal and annual data series. Table 1 lists the arithmetical evaluation results, including mean, maximum, minimum, standard deviation, kurtosis, and skewness of annual and seasonal river runoff time series (m 3 /s). These results indicate that the river runoff in spring and autumn generally exhibits positive skewness (the river runoff has a highly skewed distribution, e.g., Ton annually, Chong-Urykty in winter), and the annual river runoff also displays positive kurtosis.

Data Sources
This paper obtained the monthly and annual time series data (river runoff, covering the 1943-2012 period) selected from 13 hydrological stations surrounding Lake Issyk-Kul, as shown in Figure 1. This study compiled observational data provided by the National Science & Technology Infrastructure Center-Data Sharing Infrastructure of Earth System Science and the Department of Water Resources and Irrigation of the Ministry of Agriculture and Land Reclamation of the Kyrgyz Republic and the National Statistical Committee of the Kyrgyz Republic [63]. In this study, data from the four prominent seasons-spring (March-April-May), summer (June-July-August), autumn (September-October-November), and winter (December-January-February)-were analyzed. Also, the trends of seasonal and annual data series. Table 1 lists the arithmetical evaluation results, including mean, maximum, minimum, standard deviation, kurtosis, and skewness of annual and seasonal river runoff time series (m 3 /s). These results indicate that the river runoff in spring and autumn generally exhibits positive skewness (the river runoff has a highly skewed distribution, e.g., Ton annually, Chong-Urykty in winter), and the annual river runoff also displays positive kurtosis.

Methodology
This study used the MK trend test [40,41] and selected theŞen's trend method [46,55] to evaluate the annual and seasonal, temporal trends and features of the river runoff. The data analysis is developed in the DPS (Data Processing System) and SPSS (Statistical Product and Service Solutions) software. These methods are described in the following sections.

Trend Test
The MK trend test is a rank-based nonparametric test employed to detect significant trends in meteorological and hydrological time series [40,41]. The MK trend test statistic S is defined as: where, x k is the rank of the kth observation, x i is the rank of the ith observations, n is the number of observations, and sgn (θ) = sign (sgn), recorded as sgn (x k − x i ) and measured as: Assuming that the series with sample size n ≥ 8 is the test statistic S, then if H 0 is valid, S is approximately normally distributed, and P is the number of tied groups, and t i denotes the number of data values in the ith tied group (i = 1, 2, 3 . . . n). A tied group is a set of sample data having the same value. With variance var[S] computed as: Water 2021, 13, 95 6 of 19 After calculating the var[S] of the data series, the standardized statistic Z c value is calculated using the following formula: The standardized statistical Z c value follows the standard normal distribution (Z), with a variance of one (σ 2 = 1) and a mean value of zero (µ = 0). The null hypothesis (H 0 ) assumes the river runoff observations have no trend, while the alternate hypothesis (H 1 ) assumes the river runoff data series have no trend. Therefore, under the two-sided condition, H 0 and H 1 are tested at a significant level (α): If |Z| > 2.58, the trend is significant at the confidence level of α = 99%; if |Z| > 1.96, the trend is significant at the confidence level of α = 95%. If ± Z c > ± Z α/2 , then H 1 is accepted, and H 0 is rejected. A negative (positive) value of Z c indicates a downward (upward) trend in the observed data.

Innovative Trend Analysis Method
In this approach, a recorded observation time series is separated into two halves from the initial data to the final data, and the two sub-series are then sorted in ascending order. Based on the two-dimensional Cartesian coordinate system, the first half-series (x i ) is located on the horizontal axis (x-axis), and the second half-series (x j ) is located on the vertical axis (y-axis). The range of the two axes should be equal, as shown in Figure 2a. The 1:1 (45 • ) line divides the diagram into two similar triangles. If the data series accumulates on the 1:1 line, it is concluded that there is no trend (trendless). If the data series falls above the 1:1 line (upper triangular area), it can be assumed that the data series exhibits an upward trend. If the data series accumulates below the 1:1 line (lower triangular area), there is a downward trend in the time-series data [46,55]. Therefore, the variable domain of each semi-sequence can be classified into different patterns ("low", "medium", and "peak"). The scope of each group can be determined either qualitatively from experts' opinions or quantitatively by dividing the range of change into subcategories. The closer the scatter points are to the 1:1 (45 • ) line, the weaker the trend slope [46,55]. The straight-line trend slope calculated by the following equation [46]: where s is the straight-line slope of the trend line represented by the ITA method; n is the number of data year; and are the arithmetic averages of the first half and the second half (y, sequence) of the dependent variable. Figure 2b is a graphical representation for understanding the concept of trend slope. As shown in Figure   The straight-line trend slope calculated by the following equation [46]: where s is the straight-line slope of the trend line represented by the ITA method; n is the number of data year; y 2 and y 1 are the arithmetic averages of the first half and the second half (y, sequence) of the dependent variable. Figure 2b is a graphical representation for understanding the concept of trend slope. As shown in Figure 2b, the substitution of the numerical values as n = 70 and the arithmetic averages from Figure 2b as y 2 = 36 and y 1 = 42 in to Equation (5), yields s = 2 × (36 − 42)/70 = −0.1714. The standard deviation of sampling slope is defined as: The standard deviation of the sampling slope is σ s ; the standard deviation of whole the time series is σ, and ρ y 1 y 2 is the correlation coefficient between the ascendingly sorted two halves arithmetic averages calculates as: Finally, the lower and upper confidence limits (CL) of a standard normal PDF with zero mean and standard deviation are s cri , then the confidence limits (CL) of the trend slope at the α level of significance are calculated by applying the following equation [70]: where CL represents the lower and upper confidence limits at α level of significance; s is the slope value, and σ s is the standard deviation of sampling slope. A negative (positive) value of s indicates a decreasing (increasing) trend in the data points.

Trend Variations of Annual River Runoff in the Lake Issyk-Kul Basin
This study used the annual river runoff data from 1943-2012 at the 13 hydrological stations in the Issyk-Kul Basin to evaluate the annual and seasonal trends by employing the MK test and TSM, as presented in Table 2. These results show that the magnitudes of the trends at most of the selected river stations were positive (with the exceptions of the Pzhergalan and Tossor stations) and varied from 0.001-0.031 m 3 /s. Notably, the stations in the northern portion of the study area (i.e., Chong-Koi-Suu, Cholpon-Ata, Ak-Suu, Chong-Urykty, and Chong-Ak-Suu) exhibited large slopes (the mean slope was 0.011 m 3 /s). In contrast, the stations in the basins' southeastern parts (i.e., Pzhergalan, Karakol, Chong-Kyzyl-Suu, and Dzhukuu) displayed higher slopes (the mean absolute slope was 0.028 m 3 /s).
After the MK test, the annual river runoff trends were determined by applying the ITA method, as listed in Table 3. It can be seen from this table that the slope (s) of yearly method, verifying its ability to detect the hidden patterns in data series. The ITA method's graphical results for all 13 stations from 1943-2012 are presented in Figure 3. The increasing tendency of annual river runoff at most stations is indicated by the data points falling above the 1:1 (45 • ) line in the Cartesian coordinate system.

Trend Variations of Spring River Runoff in the Lake Issyk-Kul Basin
The river runoff trends in the spring season were identified using the classical MK trend test and TSM, and the test statistics are summarized in Table 2. These findings demonstrate that the river runoff was dominated by positive trends, with seven stations, After applying the MK trend test, the ITA method was used for the spring season for river runoff. The test results are listed in Table 4. The ITA results indicated nine stations (Cholpon-Ata, Ak-Suu, Chong-Koi-Suu, Chong-Uryukty, Chong-Kyzyl-Su, Dzhuuku, Ak-sai, Tossor, and Ton) exhibited significant positive trends. In comparison, three stations (Chong-Ak-Suu, Pzhergalan, and Tamga) displayed significant negative trends. Still, employing the MK trend test, the values had both negative and positive directions. The graphical ITA results for all the 13 river stations during the spring seasons from 1943-2012 are presented in Figure 4. As can be seen from this figure, for Tossor and Karakol, most of the data points for the spring season runoff fell on the no-trend line.

Trend Variations of Summer River Runoff in the Lake Issyk-Kul Basin
The MK trend test and TSM statistics obtained for the summer season river runoff are listed in Table 2 Table 2 summarizes the MK test statistics obtained for the autumn river runoff. As illustrated by these results, the positive values were more dominant, and some were statistically significant at the 99% confidence level. Specifically, significant positive trends were found at Cholpon-Ata (Z c = 3.31), Chong-Ak-Suu (Z c = 2.81), Ak-Suu, (Z c = 4.99), Chong-Koi-Suu (Z c = 6.37), Chong-Urykty (Z c = 3.26), Chong-Kyzyl-Suu (Z c = 6.10), Dzhuuku (Z c = 4.18), and Ton (Z c = 5.48) (61.5%), mainly scattered on the northern rivers of the Lake Issyk-Kul Basin. The summer river runoff trends were similar to the annual river runoff results due to the concentrated river runoff during summer in this region. The value at Karakol (Z c = 1.95) was also very close to the positive significance limit (Z c = 1.96). Besides, the ITA values are provided in Table 4 and graphically illustrated in Figure 6. As shown in Table 4, the positive values of slope (s) were predominated, with negative trends only found at Pzhergalan, Karakol, Tossor, and Tamga.
Water 2021, 13, x FOR PEER REVIEW 13 of 21 Figure 6. Analysis of the ITA for autumn river runoff at 13.
These results indicate that the increasing annual river runoff trends at most stations were mainly due to river runoff in autumn and winter. It is worth noting that, with the exceptions of Dzhuuku, Karakol, and Tamga, winter river runoff at all stations displayed significant increasing trends. The ITA results for all 13 river stations are graphically illus-

Trend Variations of Winter River Runoff in the Lake Issyk-Kul Basin
The MK test statistics computed for winter river runoff are listed in Table 2. These results revealed that, significant positive trends were found at nine stations (69.23%), including Pzhergalan (Z c = 6.10), Ak-Suu (Z c = 5.39), Chong-Koi-Suu (Z c = 6.10), Chong-Kyzyl-Suu (Z c = 6.28), Tossor (Z c = 3.44), Ton (Z c = 8.14), Chong-Urykty (Z c = 3.34) at the 99% level, and the Chong-Ak-Suu (Z c = 2.43) and Ak-Sai (Z c = 6.28) stations at the 95% level. In contrast, a significant negative trend was detected at Dzhuuku (Z c = −2.54) at the 95% level. The value of Cholpon-Ata (Z c = 1.90) was close to the positive significance limit (Z c = 1.96). The other two stations (Karakol and Tamga) exhibited negative trends, but neither was significant.
These results indicate that the increasing annual river runoff trends at most stations were mainly due to river runoff in autumn and winter. It is worth noting that, with the exceptions of Dzhuuku, Karakol, and Tamga, winter river runoff at all stations displayed significant increasing trends. The ITA results for all 13 river stations are graphically illustrated in Figure 7 and summarized in Table 4. Negative trends were identified at Chong-Koi-Suu and Chong-Urykty, while all other stations showed positive tendencies.
These results indicate that the increasing annual river runoff trends at most stations were mainly due to river runoff in autumn and winter. It is worth noting that, with the exceptions of Dzhuuku, Karakol, and Tamga, winter river runoff at all stations displayed significant increasing trends. The ITA results for all 13 river stations are graphically illustrated in Figure 7 and summarized in Table 4. Negative trends were identified at Chong-Koi-Suu and Chong-Urykty, while all other stations showed positive tendencies.

Comparison of Trend Analysis Methods
Verify the reliability of the seasonal and annual river runoff trends in the study area, the ITA method (Tables 3 and 4), and two well-known methods of TSM and MK trend test (Table 2) were used. According to the temporal distribution of the MK test, significantly increasing (decreasing) trends were observed in 5 (0), 6 (2), 4 (3), 8 (0), and 8 (1) time series in annual, spring, summer, autumn, and winter river runoff data. At the same time, the ITA method detected significant trends in 7 (1), 9 (3), 6(3), 9 (3), and 8 (2) time series in the study area. The significant positive and negative trends of the seasonal and annual river runoff datasets were found to display similar decreases from north to south using either the ITA method (39.2% to 29.4%; Tables 3 and 4) or the MK trend test (46.15% to 25.64%; Table 2). Corresponding graphs of these results are displayed in Figures 3-7. For the summer and winter seasons, the data series (orange and yellow points) of most hydrological stations fell above the 1:1 (45 • ) line ( Figures 5 and 7, respectively).
Besides, in Figure 8, the analysis of the ITA results (Figure 8f-j,l) has compared with the MK trend test (Figure 8a-e,k). Among the 65 annual (13 gauging stations annual data series) and seasonal data series (13 gauging stations spring, summer, autumn, and winter data series) studied, significant trends of 39 data series (60%) were found through the MK trend test. However, using the ITA method, 51 data series (78.46%) showed a significant trend. All significant trends identified by the MK trend test can also be detected by the ITA method. It can be seen from Figure 8a-e,k and Figure 8f-j,l that the trend results of the MK trend test and ITA in the spring, summer, winter, and annual runoff series are similar. However, the ITA method can be used to effectively identify many significant trends that cannot be identified by the MK trend test (especially in spring and autumn). It should be said that the ITA method is beneficial for analyzing many hidden changing trends of river runoff data series in the study area.  Table 2). Corresponding graphs of these results are displayed in Figures 3-7. For the summer and winter seasons, the data series (orange and yellow points) of most hydrological stations fell above the 1:1 (45°) line ( Figure 5 and Figure 7, respectively). Besides, in Figure 8, the analysis of the ITA results (Figure 8f-j, and l) has compared with the MK trend test (Figure 8a-e, and k). Among the 65 annual (13 gauging stations annual data series) and seasonal data series (13 gauging stations spring, summer, autumn, and winter data series) studied, significant trends of 39 data series (60%) were found through the MK trend test. However, using the ITA method, 51 data series (78.46%) showed a significant trend. All significant trends identified by the MK trend test can also be detected by the ITA method. It can be seen from Figure 8a-e, and k and Figure 8f-j, and l that the trend results of the MK trend test and ITA in the spring, summer, winter, and annual runoff series are similar. However, the ITA method can be used to effectively identify many significant trends that cannot be identified by the MK trend test (especially in spring and autumn). It should be said that the ITA method is beneficial for analyzing many hidden changing trends of river runoff data series in the study area.  Table 5 summarizes the results of the ITA method for "low", "medium", and "peak." As presented in Table 5, the "low", "medium", and "peak" values of all 13 river stations were individually evaluated by the ITA method. As provided in Table 2, there is trendless (no trend) in annual and seasonal data series according to the MK test, but increasing and decreasing trends are seen with the ITA method ( Table 5). The ITA method results revealed that increasing trends occurred in 89 of the 195 data series at the 13 hydrological stations (45.64%), while decreasing trends were found in 15 data series (7.69%). In comparison, 24 data series (26.97%) exhibited increasing patterns for "peak" river runoff, 25 data series (28.09%) displayed increasing patterns for "low" river runoff, and 40 data series (44.94%) showed increasing patterns for "medium" river runoff. According to the "low", "medium", and "peak" river runoff, five data series (33.33%), seven data series (46.67%), and three data series (20%) manifested downward trends, respectively. These results provided detailed information concerning the annual and seasonal river runoff patterns by evaluating the "low", "medium", and "peak" values. Table 5. Results of the ITA method for annual and seasonal river runoff in the "low", "medium", and "peak" in the Lake Issyk-Kul baisin .  Table 5 summarizes the results of the ITA method for "low", "medium", and "peak." As presented in Table 5, the "low", "medium", and "peak" values of all 13 river stations were individually evaluated by the ITA method. As provided in Table 2, there is trendless (no trend) in annual and seasonal data series according to the MK test, but increasing and decreasing trends are seen with the ITA method ( Table 5). The ITA method results revealed that increasing trends occurred in 89 of the 195 data series at the 13 hydrological stations (45.64%), while decreasing trends were found in 15 data series (7.69%). In comparison, 24 data series (26.97%) exhibited increasing patterns for "peak" river runoff, 25 data series (28.09%) displayed increasing patterns for "low" river runoff, and 40 data series (44.94%) showed increasing patterns for "medium" river runoff. According to the "low", "medium", and "peak" river runoff, five data series (33.33%), seven data series (46.67%), and three data series (20%) manifested downward trends, respectively. These results provided detailed information concerning the annual and seasonal river runoff patterns by evaluating the "low", "medium", and "peak" values. Table 5. Results of the ITA method for annual and seasonal river runoff in the "low", "medium", and "peak" in the Lake Issyk-Kul baisin (1943-2012).

Discussion
The change of river runoff characteristics has always been the focus and hotspot of hydrological and water resource scientific research. In this study, the ITA method presented byŞen [46] was applied at annual and seasonal time series recorded from 13 hydrological stations (Figure 1) surrounding the lake Issyk-Kul. The river is the product of climate. Suppose the weather station data describes the climatic characteristics of different locations. In that case, the river flow describes the interaction between precipitation and temperature in the river catchment area on a larger scale [71].
River discharge is determined by rainfall, which falls in the catchment area, and by evaporation, which depends on air temperature [72,73]. As reported by Alifujiang, Abuduwaili [74], the ITA method's application at three sites in the Issyk-Kul Basin revealed that the precipitation exhibited significant increasing trends in summer and winter. Moreover, utilization of the MK trend test and ITA method demonstrated that the percentage of time series with significant positive trends of monthly precipitation was decreased from north to south by Alifujiang.et.al [74]. This study revealed that in the summer and winter seasons, the data series of most hydrological stations fell above the 1:1 (45 • ) line. The significant positive and negative trends of the seasonal and annual river runoff datasets were found to display similar decreases from north to south using either the ITA method or the MK trend test, both of which indicated that river runoff is affected by precipitation. Thus, external hydrological cycle processes mainly determine the precipitation in the Issyk-Kul Basin [66].
Rivers that are mainly supplied by glacial meltwater generally have a significant increase in runoff during the summer (from June to September), with more floods during the flood season and an increase in annual runoff [5]. According to the origin of Lake Issyk-Kul, there are four types of rivers, and most of the rivers are glacial-snow type [75]. Thus, some measures should take in the summer that is vulnerable to high floods in the future. For the Lake Issyk-Kul basin, the river stations on those 13 rivers are located above the water withdrawal to irrigation. Alymkulova, Abuduwaili [65] pointed out that agricultural water demand increased in summer (from May to July or August), while in October, it rises again. Our results showed most of the river stations have an increasing trend in summer. The results revealed that the time of the high volume of water intake is not under decreasing river runoff.

Conclusions
In this study, to verify the reliability of the seasonal and annual river runoff trends in the Lake Issyk-Kul basin, Central Asia, the ITA method and the well-known method of MK trend test were used. The ITA method can be used to effectively identify many significant trends that cannot be determined by the MK trend test. It should be said that the ITA method is beneficial for analyzing many hidden changing trends of river runoff data series in the study area. The conclusions as follows: (1) The MK trend test results show that in 39-time series, there were significant positive and negative trends (among them, northern river in 18-time series, south-eastern rivers in 11-time series, southern rivers in 10-time series) on seasonal and annual river runoff. (2) The ITA method results revealed that in 51-time series (in which the stations of the northern part of the lake basin in 20-time series, south-eastern parts in 16-time series, and southern regions in 15-time series), there were significant positive and negative trends on seasonal and annual river runoff. Specifically, the MK test found that the time series percentage decreased from 46.15% in the north to 25.64% in the south, while the ITA method revealed a similar rate of decrease, from 39.2% to 29.4%. (3) According to the temporal distribution of the MK test, significantly increasing (decreasing) trends were observed in 5 (0), 6 (2), 4 (3), 8 (0), and 8 (1) time series in annual, spring, summer, autumn, and winter river runoff data. At the same time, the ITA method detected significant trends in 7 (1), 9 (3), 6(3), 9 (3), and 8 (2) time series in the study area. The comparison results revealed that the ITA method could effectively identify the trends detected by the MK trend test. (4) According to the ITA method, the "peak" values of 24 time series (26.97%) exhibited increasing patterns, 25 time series (28.09%) displayed increasing patterns for "low" values, and 40 time series (44.94%) showed increasing patterns for "medium" values. According to the "low", "medium", and "peak" values, five time series (33.33%), seven time series (46.67%), and three time series (20%) manifested decreasing trends, respectively. These results detailed the patterns of annual and seasonal river runoff data series by evaluating "low", "medium", and "peak" values.

Author Contributions:
The research article is the joint work with five authors, J.A. designed the experiments, and Y.A. analyzed the trends of river runoff variability by employingŞen's method and wrote the manuscript. Y.G. helped to review and prepare this paper for publication. All authors have read and agreed to the published version of the manuscript.