Trends of Runo ﬀ Variation and E ﬀ ects of Main Causal Factors in Mun River, Thailand During 1980–2018

: Mun River is the largest tributary of the Mekong River in Thailand and provides abundant water resources not only for an important agricultural area in Thailand but also for the lower Mekong River. To understand how the runo ﬀ of Mun River responds to climate change and human activities in recent decades, this study performed a detailed examination of the characteristics of runo ﬀ variation based on measurements at two hydrological gauging stations on the main stem of Mun River during 1980–2018. Using the Mann-Kendall test, Morlet wavelet transform and Double Cumulative Curve methods, this study identiﬁes that the variation of annual runo ﬀ of Mun River encountered an abruption in 1999 / 2000, with an increased trend taking place since then. Furthermore, a detailed assessment of the e ﬀ ects of the variations in rainfall, temperature, evaporation, and land use types extracted from remote sensing images at the basin scale reveals that a signiﬁcant reduction in forest area and slight reductions in evaporation and farmland area taking place since 1999 can lead to an increase in the runo ﬀ of Mun River, while the dramatic increase in garden area since 1999 tends to make the runo ﬀ decrease.


Introduction
Detecting trends of variations in hydrological variables and identifying main causes for the variations have been given emphasis in recent years, with an increasing scientific consensus on global and regional climate change due to human activities [1]. A detailed understanding of the variation of river runoff so as to reveal the changes in the spatio-temporal distribution and quantity of surface water resource is critical [2], since they deeply affect ecology, economy and people's livelihood along rivers [3]. The quantity and regime of river runoff depend on a variety of climate factors, including rainfall, temperature, evaporation, wind, etc. [4]. The impact of human activities is mainly reflected through land use change [5]. Current changes in river runoff vary on different scales and the report of IPCC (Intergovernmental Panel on Climate Change) in 2013 predicted that global climate change is likely to exert substantial impacts on surface water resource across the world [6]. In some regions water availability will increase but in many densely populated areas, particularly urbanizing regions with intensive conflicts among different water users, water scarcity will increase [7]. During 1936During -1998, the study of runoff variation in Yenisei River showed a significantly increasing trend, while the precipitation in the area essentially showed a decreasing trend; nevertheless, the runoff of Lena River shows an increasing trend, but the regional precipitation shows only a slight increasing trend [8]. As a result, there is no sufficient evidence to prove that the changes in runoff are determined solely by precipitation because the rainfall-runoff relationships of rivers vary from one basin to another. Some studies on the variation of annual runoff in the Yangtze River basin of China showed an insignificant increasing trend after 1960 mainly due to the rainfall brought by monsoon climate in the middle and lower reaches [9]. Yet the annual runoff in the Yellow River basin of China showed a continuously decreasing trend from 1950 to 2014, which was caused largely by human activities that took a large amount of water for irrigation and industrial production [10]. A statistically significant decreasing trend was detected based on runoff data from 925 seagoing rivers including the Mekong River from 1948 to 2004, and yet the impact of human activities on global river runoff was found far less than the effect of climate change [11]. In the prediction of future river runoff variation in Thailand and Malaysia, machine learning models and distributed hydrological models coupled with future climate scenarios are employed and the results show that the variation of forest area induced by land use change is a very influential factor [12][13][14].
As the largest tributary of the Mekong River in Thailand, Mun River inputs about 20 × 10 9 m 3 of freshwater per year to the Mekong River on average [15]. This is important for ecosystems and people living in the middle and lower reaches of the Mekong River. The Mun River basin is an important agricultural area in Thailand, where climate change and human activities exert considerable impacts on surface water resource [16]. Some researchers used statistical downscaling methods in combination with hydrological models to model and simulate future rainfall process in the Mun River and the adjacent Chi River basins, and yet found that the responsive relationships of runoff to rainfall in different climate models were different considerably [17]. The Standardized Precipitation and Evapo-transpiration Index (SPEI) was used to assess the impacts of recent climate variations on rice yields in the Mun River basin, and the results demonstrate that the drought index SPEI detects soil moisture deficiency and crop stress in rice better than precipitation or precipitation based indicators [18]. It is clear that a thorough assessment of runoff variability in the Mun River basin is still lacking. In addition, the entire Lancang-Mekong River basin has undergone rapid urbanization in recent decades while being dominated traditionally with agriculture and so it is urgent to know if Mun River [19], the largest tributary of the Mekong River in Thailand, has or is going to supply lesser amounts of freshwater to the lower Mekong River.
The main objectives of this study are twofold and the first is to gain knowledge on if the amount of freshwater input to the lower Mekong River from Mun River has reduced. The second objective is to identify main factors of climate change and human activities that play important roles in influencing the variation of runoff of Mun River over the last several decades. To achieve the objectives, this study presents a detailed evaluation of the characterization of runoff variation in Mun River based on long-term measurements at local weather and hydrological gauging stations. In addition, a detailed analysis is performed to determine the changes of local land use types in different time periods using remote sensing images and consequently the possible influences of the changes of local climate factors and different land use types on the variation of runoff of Mun River are evaluated.

Research Area
Mun River, one of the large tributaries of the Mekong River, is situated in the northeastern Thailand and lies between latitudes 14 • and 16 • and longitudes 101 • 30' and 105 • 30 , with a length of approximately 760 km and a basin area of 82,000 km 2 , respectively [20] (Figure 1). The river originates on the Khorat plateau near Nakhon Ratchasima in the northwest and flows towards east through southern flat area until it joins the Mekong River in Ubon Ratchathani province, Thailand [21]. The basin is an important agricultural area in Thailand, about 80%, 10%, and 5% of the total land area being farmland, forest and garden, respectively, in 2018 [22][23][24]. Water 2020, 12, x FOR PEER REVIEW 3 of 17 The Mun River basin has a tropical savanna climate which is most significantly affected by tropical monsoons in Asia [25]. The climate and hydrology within the basin show significant seasonal differences brought upon by the seasonal monsoons, forming distinct rainy (flooding) and dry seasons that last from June to October, and November to May of the following year, respectively [15]. Rainfall in the rainy season accounts for 92-98% of the total annual rainfall and only 2-8% in the dry season. The peak rainfall generally occurs during August to September. The spatial distribution of rainfall is gradually increasing from west to east, and the average annual rainfall is 1300-1500 mm ( Figure 2). The annual temperature in the basin is largely higher than 18°C, with an annual average temperature of about 25°C [26].

Data Collection
Hydrological monitoring in the Mun River basin has been conducted since 1980 at CHUMPHON (simply as CH) station in the middle reach and UBON (simply as UB) station in the lower reach, and Figure 1 shows the locations of the two hydrological gauging stations on the main stem of Mun River. Meteorological stations (about 159 stations) in the Mun River basin were set up since 1950. Due to data availability, daily rainfall data at 9 meteorological stations (including Lam The Mun River basin has a tropical savanna climate which is most significantly affected by tropical monsoons in Asia [25]. The climate and hydrology within the basin show significant seasonal differences brought upon by the seasonal monsoons, forming distinct rainy (flooding) and dry seasons that last from June to October, and November to May of the following year, respectively [15]. Rainfall in the rainy season accounts for 92-98% of the total annual rainfall and only 2-8% in the dry season. The peak rainfall generally occurs during August to September. The spatial distribution of rainfall is gradually increasing from west to east, and the average annual rainfall is 1300-1500 mm ( Figure 2). The annual temperature in the basin is largely higher than 18 • C, with an annual average temperature of about 25 • C [26]. The Mun River basin has a tropical savanna climate which is most significantly affected by tropical monsoons in Asia [25]. The climate and hydrology within the basin show significant seasonal differences brought upon by the seasonal monsoons, forming distinct rainy (flooding) and dry seasons that last from June to October, and November to May of the following year, respectively [15]. Rainfall in the rainy season accounts for 92-98% of the total annual rainfall and only 2-8% in the dry season. The peak rainfall generally occurs during August to September. The spatial distribution of rainfall is gradually increasing from west to east, and the average annual rainfall is 1300-1500 mm ( Figure 2). The annual temperature in the basin is largely higher than 18°C, with an annual average temperature of about 25°C [26].

Data Collection
Hydrological monitoring in the Mun River basin has been conducted since 1980 at CHUMPHON (simply as CH) station in the middle reach and UBON (simply as UB) station in the lower reach, and Figure 1 shows the locations of the two hydrological gauging stations on the main stem of Mun River. Meteorological stations (about 159 stations) in the Mun River basin were set up since 1950. Due to data availability, daily rainfall data at 9 meteorological stations (including Lam

Data Collection
Hydrological monitoring in the Mun River basin has been conducted since 1980 at CHUMPHON (simply as CH) station in the middle reach and UBON (simply as UB) station in the lower reach, and Figure 1 shows the locations of the two hydrological gauging stations on the main stem of Mun River. Meteorological stations (about 159 stations) in the Mun River basin were set up since 1950.
Due to data availability, daily rainfall data at 9 meteorological stations (including Lam Phra Phloeng, Makham Thao, Phi Mai, Mun River, Nong Tha Chok, Huai Lam Phok, Huai Wang Daeng, Ubon, and Regional R.I.D) and daily evaporation data at 13 stations measured respectively since 1960 and 1980 were selected. All of this information is provided by the Royal Irrigation Department (R.I.D) of, Thailand. The digital elevation model used in this study is obtained from the Geospatial Data Cloud (China), while temperature data at 3 stations and Land use data were obtained respectively from Globe Summary of the Day (GSOD) and Servir Mekong website (Table 1).

Methods
To investigate the temporal variation of runoff in Mun River, both parametric and non-parametric test methods were employed [27]. As a non-parametric method, the Mann-Kendall trend test is effective in dealing with long time-series hydrological and meteorological data that are not in normal distributions [28]. It can identify quantitatively whether a certain runoff change process has an abrupt point or has a stable change trend [29]. Comparatively, as one of the parametric test methods, the Double Cumulative Curve method is applicable to analyze the relationship between rainfall and runoff [30,31]. In addition, the spatial Kriging interpolation method is used to examine the spatial distribution of rainfall in terms of the calculated results using Thiessen polygon method, while the wavelet analysis is applied to study the periodicity of rainfall and runoff [32][33][34][35].

Mann-Kendall Non-Parametric Trend Test
For time series H 0 containing n samples (x 1 , x 2, . . . x n ), which are independent and random, the Mann-Kendall statistic S is given as The application of the trend test is conducted for a time series x i that is ranked from i = 1,2,...n−1 and x j , which is ranked from j = I + 1,2,...n. Each of the data point x i is taken as a reference point which is compared with the rest of the data points x j so that It has been documented that when n ≥ 8, the statistic S is approximately normally distributed with the mean [36,37]: Water 2020, 12, 831 5 of 16 The test statistics Z is computed as where Z follows a standard normal distribution. A positive (negative) value of Z signifies an upward (downward) trend. A significance level α is also utilized for testing either an upward or downward monotone trend (a two-tailed test). If |Z| > 1.28, 1.64 and 2.32, then the confidence test of 90%, 95%, and 99% is passed, respectively.

Mann-Kendall Non-Parametric Mutation Test
To perform a Mann-Kendall non-parametric mutation test, the rank column U m of the sequence needs to be constructed in the following manner: When the sample value i is greater than the value at the moment j, the number of values will be accumulated. When x 1 , x 2 , and x n are independently and continuously distributed, the following statistics can be determined: Var(U m ) = n * (n − 1) * (2n + 5) 72 (10) where UF m is a standard normal distribution and UF 1 = 0, and E(U m ) and Var(U m ) are the mean and variance of a cumulative count of U m , respectively. For a given significance level, the critical value U a can be obtained by querying the normal distribution table. If |UF i | > U a (U 0.05 = 1.96, U 0.01 = 2.58), it means that the sequence has an obvious trend of increase or decrease [38]. Repeat the calculation process above to obtain the inverted sequence statistical value UB m and when the UF m and UB m curves intersected between the critical boundary, the intersection point can be taken as the beginning of mutation [39].

Morlet Wavelet Transform Method
Morlet wavelet function is a kind of complex-valued and continuous wavelet and so continuous wavelet transform is more suitable for signal feature extraction and further analysis [40,41]. Complex-valued wavelet gives the information of phase and amplitude of time series and is able to eliminate the false oscillation caused by real wavelet transform coefficient, so that more accurate results can be obtained. The detailed expression of the method is: Water 2020, 12, 831 where f b is the bandwidth, f c is the central frequency, and i is the imaginary part. The expression of the continuous wavelet transform of f(t) takes a form of: (12) In order to further determine the main period of the time series, the wavelet squared difference can be obtained by integrating the square of all wavelet coefficients about alpha through formula (12) [42].
The wavelet squared diagram reflects the distribution of signal wave energy with the time scale. While the wavelet squared difference determines the signal peak value, the year corresponding to the peak value is regarded as the main cycle [43]:

Double Cumulative Curve Method
The method of Double Cumulative Curve has been used in the analysis of long-term trend of hydrological sequence [44]. Combined with Mann-Kendall test, this method can effectively identify the variation period of the slope of a double cumulative curve [45]. Theoretically, the slope of the double cumulative curve for runoff against rainfall will not change significantly when only rainfall takes effect on runoff, while the slope may change when other factors such as land use and other human activity factors exert an influence on runoff. According to the time and variation of the slope, the influencing degree of a certain factor on runoff can be evaluated.
Based on the mutation year of a double cumulative curve and choosing a period with no or little significant slope offset as a reference period (others are changing periods), the following linear regression equation between cumulative runoff and cumulative rainfall can be established: where Q is the runoff and P is the rainfall.
Assuming that the trend of climate change in a changing period is consistent with that in the reference period [46], the contributions of rainfall and human activities to runoff variation can then be calculated by using the following formulas, respectively: where ∆Q h is the quantity of runoff influenced by human activities, ∆Q is the contribution of rainfall, subscript 1 and 2 represent, respectively, the reference period and changing period, i stands for the value calculated, and j means the value measured.

Spatial Kriging Interpolation Method
Similar to Inverse Distance Weighting (IDW) method, Kriging method uses a weighting, which assigns more influence to the nearest data points in the interpolation of values for unknown locations [47]. Kriging method considers not only the relationship between the positions of the target points and the given points, but also the spatial correlation of the variables. The two-step process of Kriging method begins with semi variance estimations and then performs the interpolation [48].
The aim of Kriging method is to estimate the value of an unknown real-valued function Z at a point (m, n), with the given values of the function at some other points {(m 1 , n 1 ), (m 2 , n 2 )... (m 3 , n 3 )}.
Water 2020, 12, 831 7 of 16 A Kriging estimator is assumed to be linear because the predicted value Z(m, n) is a linear combination that can be written as: where γ i is a weight coefficient of the solutions of a system of linear equations which minimize the difference between the estimated value Z 0 and the true value Z 0 at the point (m 0 , n 0 ), while ε (m, n) is the error of prediction as: where random process F(m, n) also satisfies the condition of unbiased estimation:

Annual and Interannual Variation of Runoff
where is a weight coefficient of the solutions of a system of linear equations which minimize the difference between the estimated value and the true value at the point (m0, n0), while (m, n) is the error of prediction as: where random process F(m, n) also satisfies the condition of unbiased estimation: Figure 3 shows the distribution of monthly average flow discharge at CH and UB, two hydrological stations located in the main stem of Mun River, based on flow data measured during 1980 to 2018. The low and high flow discharges at the two stations in the study period all exhibit seasonal distribution characteristics. The rainy season is from August to November, during which the ratio of the maximum monthly flow discharge to the minimum can reach up to 600:1 and the runoff in the season accounts for about 77% of the total annual runoff. The monthly average flow discharge in the rainy and dry seasons in the Mun River basin fluctuated in wide ranges. The average flow discharge at CH station in the rainy season is 99.6 m³/s, 2.4 times of the counterpart in the dry season (Figure 3a). The annual average flow discharge at UB station in the rainy season is 1136.4 m³/s, four times of the counterpart in the dry season (Figure 3b). The average annual flow discharge at UB station is 640.7 m³/s, almost ten times of the counterpart at CH station in the upper reach. The peak flow discharge occurs in general from September to October, slightly behind the peak of the rainfall.  The average annual runoff at CH and UB stations in the period 1980-2018 is 2.18 × 10 9 m³ and 19.5 × 10 9 m³, with a runoff yield per square kilometer of 76.8 × 10 3 m³ and 183.17 × 10 3 m³, respectively ( Figure 4). Based on the Mann-Kendall trend test (Table 2), though no significant increasing or decreasing trend can be detected from the runoff variation during 1980-2018, a significantly decreasing trend and an increasing trend in the variation of the average annual runoff can be observed for the period before 1999 and after 2000, respectively. The Cv (coefficient of  (Table 2), though no significant increasing or decreasing trend can be detected from the runoff variation during 1980-2018, a significantly decreasing trend and an increasing trend in the variation of the average annual runoff can be observed for the period before 1999 and after 2000, respectively. The Cv (coefficient of variation) of the annual runoff at CH and UB stations over the period 1980-2018 is 0.4 and 0.6, respectively, while the annual runoff at variation of rainfall in the whole Mun River basin.  According to the results of an abruption analysis using Mann-Kendall test (Figure 5a) and the cumulative anomaly runoff change curve between CH and UB stations (Figure 5b   According to the results of an abruption analysis using Mann-Kendall test (Figure 5a) and the cumulative anomaly runoff change curve between CH and UB stations (Figure 5b variation) of the annual runoff at CH and UB stations over the period 1980-2018 is 0.4 and 0.6, respectively, while the annual runoff at both CH and UB stations reached maximum values in 2000. This general variation trend and high variability of annual runoff are highly consistent with the variation of rainfall in the whole Mun River basin.   According to the results of an abruption analysis using Mann-Kendall test (Figure 5a) and the cumulative anomaly runoff change curve between CH and UB stations (Figure 5b Figures 6 and 7 show the real part of runoff wavelet transform coefficients and wavelet square differences at CH and UB stations (blue polylines), and the real part of rainfall coefficients at three representative weather stations along the main stem of Mun River (Phi Mai Barrage station, Mun River station and Regional R.I.D shown in Figure 1). (Figure 7).

Temporal Variation of Runoff and Rainfall
The runoff and rainfall periodicities in the Mun River basin are summarized in Table 3, in which the period with a duration of 16-0 years is classified as long-term scale, the period of 10-15 years as intermediate scale, and the period of less than 10 years as short-term scale, corresponding to long, middle and short periods, respectively. It can be seen from Table 2 that in the Mun River basin the periodic changes of rainfall and runoff exhibit similar patterns.  corresponding to cycles of 5, 11 and 15 years, among which the 11-year cycle is the first main cycle (Figure 7). The runoff and rainfall periodicities in the Mun River basin are summarized in Table 3, in which the period with a duration of 16-0 years is classified as long-term scale, the period of 10-15 years as intermediate scale, and the period of less than 10 years as short-term scale, corresponding to long, middle and short periods, respectively. It can be seen from Table 2 that in the Mun River basin the periodic changes of rainfall and runoff exhibit similar patterns.  It can be noted from Figure 6 that the variation of annual runoff in the Mun River basin exhibits characteristic timescales of 6-8, 15-16 and 30 years, respectively, corresponding to the cycles of 6, 11 and 15 years, among which the 11-year cycle is the first main cycle. The annual rainfall in the Mun River basin also exhibits three characteristic timescales (8, 14-20 and 30 years), respectively corresponding to cycles of 5, 11 and 15 years, among which the 11-year cycle is the first main cycle (Figure 7).
The runoff and rainfall periodicities in the Mun River basin are summarized in Table 3, in which the period with a duration of 16-0 years is classified as long-term scale, the period of 10-15 years as intermediate scale, and the period of less than 10 years as short-term scale, corresponding to long, middle and short periods, respectively. It can be seen from Table 2 that in the Mun River basin the periodic changes of rainfall and runoff exhibit similar patterns. In comparison with Period A, Figure 8a shows that the value of slope in the double cumulative curve of runoff against rainfall is larger than in Period B, meaning an increase in runoff relative to the change in rainfall in Period B. Moreover, Table 4 shows that runoff and rainfall in period B increases 23.1% and 12.1%, respectively, at two stations, while the average annual rainfall increases slightly at one station (6.2%) and yet decreases insignificantly at the other station (−0.2%). These inconsistencies in the relationship of runoff and rainfall in the Mun River basin demonstrate clearly that the slight increase in annual runoff change in Period B in Mun River is caused not by the change in annual rainfall but by the other factors. Figure 8 presents the double cumulative curve of runoff against rainfall and the cumulative curves of runoff and rainfall in the Mun River basin in the period 1980-2015, respectively. According to the results obtained earlier and shown in Table 3 and Figures 2, 6, and 7, the runoff sequence in the Mun River basin can be divided into two periods: Period A (1980-1999) and Period B (2000)(2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015). In comparison with Period A, Figure 8a shows that the value of slope in the double cumulative curve of runoff against rainfall is larger than in Period B, meaning an increase in runoff relative to the change in rainfall in Period B. Moreover, Table 4 shows that runoff and rainfall in period B increases 23.1% and 12.1%, respectively, at two stations, while the average annual rainfall increases slightly at one station (6.2%) and yet decreases insignificantly at the other station (−0.2%). These inconsistencies in the relationship of runoff and rainfall in the Mun River basin demonstrate clearly that the slight increase in annual runoff change in Period B in Mun River is caused not by the change in annual rainfall but by the other factors.

Temperature and Evaporation
The average, maximum, and minimum temperature variations at three representative meteorological stations (i.e., Nakhon Ratchasima, Surin, and Ubon Rathathani stations) located respectively in the upper, middle, and lower reaches of Mun River are shown in Figure 9a. The annual evaporation change is shown in Figure 9b. Although slightly increasing trends are exhibited in the variations of all three variables, the corresponding variation ranges are very small (Figure 9a). In comparison with the average value of evaporation in Period A, the average value in Period B reduces about 100 mm, or −6% (Figure 9b). According to the water balance principle in a river basin, the reduction of evaporation normally exerts a positive effect on the increase of runoff [49,50]. In the Mun River basin, hence, the evaporation decrease can benefit the recovery of runoff, although not to a very significant degree because only a very minor reduction takes place in evaporation.
The average, maximum, and minimum temperature variations at three representative meteorological stations (i.e., Nakhon Ratchasima, Surin, and Ubon Rathathani stations) located respectively in the upper, middle, and lower reaches of Mun River are shown in Figure 9a. The annual evaporation change is shown in Figure 9b. Although slightly increasing trends are exhibited in the variations of all three variables, the corresponding variation ranges are very small ( Figure  9a). In comparison with the average value of evaporation in Period A, the average value in Period B reduces about 100 mm, or −6% (Figure 9b). According to the water balance principle in a river basin, the reduction of evaporation normally exerts a positive effect on the increase of runoff [49,50]. In the Mun River basin, hence, the evaporation decrease can benefit the recovery of runoff, although not to a very significant degree because only a very minor reduction takes place in evaporation.   Table 5 shows the areas of different land use types in typical years during 1987-2018 and the area changes during 1987-2000 and 2000-2018, respectively. It can be noticed that farmland and wetland changed only slightly during the two periods. In contrast, the area of forest decreased significantly, while the areas of grassland, garden and settlement all increased very significantly during the two periods. Because the areas of settlement, grassland and other land use types in total occupy a very small proportion of the whole land area, the most significant area changes occurred in forest and garden during 2000-2018, with the former decreasing significantly and the latter  Table 5 shows the areas of different land use types in typical years during 1987-2018 and the area changes during 1987-2000 and 2000-2018, respectively. It can be noticed that farmland and wetland changed only slightly during the two periods. In contrast, the area of forest decreased significantly, while the areas of grassland, garden and settlement all increased very significantly during the two periods. Because the areas of settlement, grassland and other land use types in total occupy a very small proportion of the whole land area, the most significant area changes occurred in forest and garden during 2000-2018, with the former decreasing significantly and the latter increasing dramatically. Because the farmland in the Mun River basin is largely in the forms of rice, sugarcane and cassava plantation, an increase in farmland area can induce higher water deprivation and vice versa [51]. While there are studies that demonstrate that a reduction in the area of forest or farmland can lead to an increase in runoff due to the reduction in the water storage capacity of the land [52,53], an increase in the area of garden commonly leads to the use of more water and so can result in a reduction in runoff [12,54]. Hence, it is clear that the significant reduction in forest area and the slight reduction in farmland area in the Mun River basin during 2000-2018 helped to increase the runoff of Mun River, while the dramatic increase in garden area tends to reduce the runoff.

Conclusions
As the largest tributary of the Mekong River in Thailand, the Mun River provides abundant water resources to the lower Mekong River. Nevertheless, the Mun River basin is an important agricultural area in Thailand and has been subject to climate change and human activities, in particular in recent decades. To understand how the runoff of Mun River responds to climate change and human activities, this study performed a detailed examination of the characteristics of runoff variation based on measurements at two hydrological gauging stations on the main stem of Mun River during 1980-2018. Importantly, this study also presents a detailed investigation of the effects of potential factors on the variation of runoff at the basin scale. These detailed analyses yield the following results:

1.
In the variation of annual runoff of Mun River during 1980-2018, an abruption occurred in 1999/2000 and since then the variation of annual runoff of Mun River takes an increased trend; 2.
The variations of rainfall and runoff across the Mun River basin exhibit similar periodic patterns, while the average annual runoff in Period B (during 2000-2018) is significantly larger than in Period A (during 1980-1999) and the average annual rainfall in the basin varies in a consistent trend in the two periods; 3.
The average, maximum and minimum temperature observed at representative meteorological stations in the upper, middle and lower drain basin vary all in slightly increasing trends, while the average evaporation in Period B reduces slightly; 4.
The land use types in the Mun River basin vary in very different forms during 1987-2015, with forest area showing a significant decrease, garden area increasing dramatically, farmland area reducing slightly, and the other types (settlement, wetland, grassland, etc.) changing insignificantly during 2000-2015.
5. The significant reduction in forest area and slight reductions in evaporation and farmland area since 1999 all help to increase the runoff of Mun River, although the dramatic increase in garden area since 1999 tends to decrease the runoff.
Under the effects of climate change and human activities, a very large number of rivers have encountered the situation of a significant reduction in runoff output. Hence, it is very interesting to know that the runoff of Mun River to the Mekong River has increased since 2000. Although this study identifies that the significant reduction in forest area and slight reductions in evaporation and farmland area since 1999 all helped to increase the runoff of Mun River, the findings are based on qualitative inferences and a more detailed quantitative assessment using basin-scale hydrologic models and other approaches is required in future studies.

Conflicts of Interest:
The authors declare no conflict of interest.