Next Article in Journal
The Effect of Woody Encroachment on Evapotranspiration in a Semi-Arid Savanna
Previous Article in Journal
Soil Erosion Quantification using Machine Learning in Sub-Watersheds of Northern Portugal
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Trends and Variabilities in Rainfall and Streamflow: A Case Study of the Nilwala River Basin in Sri Lanka

by
Ravindu Panditharathne
1,2,
Miyuru B. Gunathilake
3,4,*,
Imiya M. Chathuranika
1,
Upaka Rathnayake
1,
Mukand S. Babel
5 and
Manoj K. Jha
6
1
Department of Civil Engineering, Faculty of Engineering, Sri Lanka Institute of Information Technology, New Kandy Road, Malabe 10115, Sri Lanka
2
Faculty of Technology, Wayamba University of Sri Lanka, Kuliyapitiya 27411, Sri Lanka
3
Hydrology and Aquatic Environment, Division of Environment and Natural Resources, Norwegian Institute of Bioeconomy and Research, 1430 Ås, Norway
4
Water, Energy and Environmental Engineering, Faculty of Technology, University of Oulu, 90570 Oulu, Finland
5
Water Engineering and Management, Asian Institute of Technology, Pathumthani 12120, Thailand
6
Civil, Architectural and Environmental Engineering, 456 McNair Hall, North Carolina Agricultural and Technical State University, Greensboro, NC 27411, USA
*
Author to whom correspondence should be addressed.
Hydrology 2023, 10(1), 8; https://doi.org/10.3390/hydrology10010008
Submission received: 6 November 2022 / Revised: 21 December 2022 / Accepted: 23 December 2022 / Published: 29 December 2022
(This article belongs to the Topic Hydrology and Water Resources in Agriculture and Ecology)

Abstract

:
Rainfall is one of the dominating climatic parameters that affect water availability. Trend analysis is of paramount significance to understand the behavior of hydrological and climatic variables over a long timescale. The main aim of the present study was to identify trends and analyze existing linkages between rainfall and streamflow in the Nilwala River Basin (NRB) of Southern Sri Lanka. An investigation of the trends, detection of change points and streamflow alteration, and linkage between rainfall and streamflow were carried out using the Mann–Kendall test, Sen’s slope test, Pettitt’s test, indicators of hydrological alteration (IHA), and Pearson’s correlation test. Selected rainfall-related extreme climatic indices, namely, CDD, CWD, PRCPTOT, R25, and Rx5, were calculated using the RClimdex software. Trend analysis of rainfall data and extreme rainfall indices demonstrated few statistically significant trends at the monthly, seasonal, and annual scales, while streamflow data showed non-significant trends, except for December. Pettitt’s test showed that Dampahala had a higher number of statistically significant change points among the six rainfall stations. The Pearson coefficient correlation showed a strong-to–very-strong positive relationship between rainfall and streamflow. Generally, both rainfall and streamflow showed non-significant trend patterns in the NRB, suggesting that rainfall had a higher impact on streamflow patterns in the basin. The historical trends of extreme climatic indices suggested that the NRB did not experience extreme climates. The results of the present study will provide valuable information for water resource planning, flood and disaster mitigation, agricultural operations planning, and hydropower generation in the NRB.

1. Introduction

Climate is a key factor that affects environmental systems, socioeconomic conditions, and water resource availability [1]. The changes in rainfall patterns will directly affect streamflow and thereby domestic, agricultural, and industrial water needs [2]. Moreover, streamflow will also be affected by anthropogenic activities [2,3], such as land-use change, operation of dams and reservoirs, and direct water extraction from surface water and groundwater systems [4]. Hence, identifying and analyzing the long-term trends of meteorological and hydrologic data will be useful for water resource planning and management [5], flood protection and disaster mitigation [3,6], and agricultural operations [2]. Trend analysis will be valuable to eliminate errors in approximations in designing hydraulic structures under assumed fixed hydrometeorological variables [2].
Many studies in different geographic regions of the world were directed toward identifying trends and variabilities in rainfall and streamflow and their associated linkages [2,3,6,7,8,9,10,11]. Mersin et al. [10] stated that the variations in the frequency and magnitude of rainfall caused biotic and abiotic disturbances in the environment. Kastridis et al. [11] investigated the relationship between climate and tree growth for a tree species called A. Borisii-regis in the Mediterranean. They found that rainfall was the key driving factor for tree growth during the study period. Ademe et al. [2] demonstrated that the change in the water flow of the Birr River in Ethiopia was not only influenced by the change in rainfall but was also due to changes in land cover and land use, as well as human interventions, such as upstream water abstraction. In another study by Chaluka et al. [3], it was found that changes in rainfall influenced the alterations in streamflow patterns. Bellabas et al. [8] used a climate elasticity model and a hydrologic model to examine the effects of anthropogenic activities and changes in climate on streamflow. The results revealed that anthropogenic reasons were the dominant causes for the alterations in streamflow. In contrast to the above results, several others, such as Hannaford [12] and Wang et al. [13], found that the variations in rainfall significantly influenced streamflow patterns. Moreover, studies such as those by Azari et al. [14], Dey and Mishra [15], and Xu et al. [16] found that climate change had impacts on streamflow changes to varying degrees. Most of the trend analysis studies [2,3,17] used the Mann–Kendall test and Sen’s slope estimator to study rainfall and streamflow trends. Pettitt’s test was used for the detection of changing points in a hydrometeorological time series [18,19,20,21]. Other trend analysis methods, such as Spearman’s rho and the linear regression test, were used by Fentaw et al. [9] and Coloiero [22]. However, some of these should be performed under certain assumptions, for instance, when the data is normally distributed and there are specific criteria on the length of the data series [22]. The indicators of hydrologic alteration (IHA) are commonly used to identify the hydrological impacts of human activities and to provide recommendations for environmental flow management [23,24,25].
Sri Lanka is an agrarian country that is highly dependent on rainfed and irrigation water. According to the Annual Report of the Central Bank of Sri Lanka, in 2021, the agriculture sector contributed 6.9% of the gross domestic product. Nearly 27.3% of Sri Lankans’ engage in the agricultural sector as their livelihood. Sri Lanka experiences two major monsoon periods, which are the northeast monsoon (NEM; December to February) and the southwest monsoon (SWM; May to September). The two inter-monsoon periods are the first inter-monsoon (FIM; March to April) and the second inter-monsoon (SIM; October to November) [26].
Several studies, including Abeysingha [27], Perera et al. [28], Alahacoon and Edirisinghe [29], and Ruwangika et al. [30], studied rainfall and streamflow trends in Sri Lanka. These studies identified an increasing rainfall trend over the country that was most prominent in the eastern, southeastern, north, and north–central areas. Jayasekara and Abeysingha [17] found that there was a significant association between streamflow and rainfall variations for 70% of gauging stations in the Kelani River Basin. Chathuranika et al. [26] found that the climate and streamflow conditions of the Nilwala River Basin are expected to change in the future relative to the current conditions. Even though rainfall and streamflow trends studies were carried out in Sri Lanka, a handful of them focused on extreme rainfall indices, while none of the document studies used IHA parameters to assess the shifts in hydrologic regimes. Rainfall and streamflow trends and variabilities and their existing linkages have not been assessed for the Nilwala River Basin (NRB), which is one of the major river basins in the southern part of the island. Therefore, this study aimed to analyze long-term rainfall and streamflow trends, detect change points, and identify hydrological variables and their linkage over the NRB. The findings of this study will be helpful for both public and private sectors that are involved in water resource planning and development, disaster management, agricultural development, etc.

2. Materials and Methods

2.1. Study Area

The Nilwala River Basin (NRB) is in the southern part of the country between latitudes 5°55′ N and 6°13′ N and between longitudes 80°25′ E and 80°38′ E [4]. The river originates from Panilkanda in Deniyaya at an altitude of 988 m above the MSL (mean sea level), flows about 72 km through agricultural, urban, and other land uses, and finally drains into the Indian Ocean in Matara [31]. The total basin area is about 1010 km2. The annual discharge of the river is more than 800 million cubic meters (MCM). The mean annual precipitation of the upper part of the NRB is about 3000 mm, while in the lower part, it is about 1900 mm [26]. Figure 1 demonstrates rainfall and hydrological stations in the NRB.

2.2. Rainfall and Streamflow Data

2.2.1. Observed Data

Daily observed rainfall data from 1991–2014 (24 years) for six stations, namely, Dampahala, Kamburupitiya, Kekenadura, Kirama, Goluwawatta, and Deniyaya, were collected from the Department of Meteorology of Sri Lanka. Daily observed discharge data for the Pitabeddara hydrological station was collected from the Department of Irrigation of Sri Lanka for the same period.

2.2.2. Gridded Data

Missing rainfall data were filled using two versions of the Asian Precipitation-Highly Resolved Observational Data Integration towards Evaluation of Water Resources (APHRODITE) products, namely, V1101_MA and V1901_MA, with the same resolutions of 0.25° × 0.25°. APHRODITE products are available in three main geographical domains: Monsoon Asia (MA), Middle East (ME), and Russia (RU) (http://aphrodite.st.hirosaki-u.ac.jp/, accessed on 2 August 2022) [26]. The high temporal and spatial resolution in APHRODITE when compared with other gridded-based products and well-developed quality control methods influenced researchers to use this gridded product. Yatagai et al. [32], Yatagai et al. [33], and many other studies used APHRODITE data products for climatological studies, to validate satellite data, and to downscale low-resolution model data. Table 1 provides the general information of the stations used in the present study.

2.3. Methodology

Initially, the data quality and consistency were checked, and missing rainfall data were filled in using the APHRODITE precipitation data. Both hydrological and meteorological data were categorized into monthly, seasonal, and annual timescales. Mann–Kendall, modified Mann–Kendall, Sen’s slope estimator, and Pettitt’s tests were performed to identify the trends in the rainfall data and streamflow data to compute their magnitudes and to detect change points in the time series data using XLSTAT software (available at https://www.xlstat.com/en/, accessed on 5 August 2022) [34]. In addition, 5 extreme rainfall indices were calculated using the RClimdex software (available at https://www.climdex.org/, accessed on 9 August 2022) [35] and computed trends and magnitudes were found using MK and Sen’s slope tests. The correlation between rainfall and streamflow was analyzed using Pearson’s correlation coefficient. Thereafter, indicators of hydrological alteration (IHA) software was used to analyze the variations of the 16 selected hydrological parameters. Figure 2 below demonstrates the overall methodology of the study.

2.4. Trend Analysis Methods

2.4.1. Mann–Kendall Test

The non-parametric Mann-Kendall test [36,37] was used to identify the rainfall trends. The Mann–Kendall test static S was calculated using the following Equations (1) and (2):
S = i = 1 n 1 j = i + 1 n s g n ( x j x i ) ,
s g n ( x j x i ) = { + 1       i f   x j > x i 0         i f   x j = x i 1       i f   x j < x i
where xj and xi are the data values at times j and i (j > i) and n is the length of the data series.
The Mann–Kendall Z statics and the variance Var(S) were calculated using the following Equations (3) and (4):
Z = { s 1 V a r ( S ) 0 s + 1 V a r ( S ) ,             i f [ S > 0 S = 0 S < 0 ]
The Mann–Kendall test accepts the null hypothesis if −ZZcrZ is the critical value of the normalized statics Z at a 5% confidence level. The negative and positive values of Z and S statistics indicate decreasing and increasing tendencies, respectively.
V a r ( S ) = 1 18 [ n ( n 1 ) ( 2 n + 5 ) ]

2.4.2. Modified Mann–Kendall Test

Importantly, the influence of the accuracy of values in time series data on each other is highly affected by the serial correlation. The serial correlation is present in hydrologic data sets, such as streamflow and water levels [38,39,40,41]. The presence of a positive serial correlation will increase the ability to show trends from a certain level other than the actual status, and sometimes it shows trends that have no significant trends in the actual scenario [42]. Therefore, a serial correlation check should be conducted before applying the MK test. If a serial correlation is present, the modified Mann–Kendall test proposed by Hamed and Rao [42] should be applied to the time series data to eliminate the effects of serial correlation. Trend analysis of the hydrological data was checked using the modified Mann–Kendall (MMK) test and the modified variance can be calculated using the following Equations (5) and (6):
V a r ( S ) = V a r ( S ) n n
where Var(S)* is the modified variance and the correction factor n/n*’ is given by the following Equation (6):
n n = 1 + 2 n ( n 1 ) ( n 2 ) j = 1 n 1 ( n k ) ( n k 1 ) ( n k 2 ) r k R
where n is the actual number of observations, n is the effective number of observations to count for the auto-correlation data, and r k R is the lag k autocorrelation coefficient of the rank of data.

2.5. Sen’s Slope Estimator

Sen’s slope estimator [43] is also a non-parametric test that is widely used to compute the magnitude of trend series. The slope of the data series is computed using the following Equations (7) and (8):
T i = X j X k j k       For i = 1 ,   2 ,   3 N
β = { T N + 1 2             N   i s   o d d 1 2 ( T N 2 + T N + 1 2 )     N   i s   e v e n
where xj and xi are data values at the time j and k (j < k). N is the number of data pairs (xj, xk), 1 ≤ j < k ≤ n. β defines Sen’s slope value. Positive β values signify an increasing trend and negative β values signify a decreasing trend

2.6. Change Point Detection

The nonparametric Pettitt’s test developed by Pettitt [44] was used to identify abrupt changes in time rainfall data within the study area. The method was derived using the Mann–Whitney statistic U t ,   n and two test samples from the same population, namely, x 1 , …, x t and x t + 1 , …, x n . The test statistic U t ,   n can be obtained using the following Equations (9) and (10):
U t ,   n = U t 1 ,   n + J = 1 n s g n ( x t x j )
where t = 2, 3, …, n; n is the length of the time series; and
If   ( x t x j ) > 0 ,   s g n ( x t x j ) = 1 If   ( x t x j ) = 0 ,   s g n ( x t x ) = 0 If   ( x t x j ) < 0 ,   s g n ( x t x j ) < 1
The test statistic quantifies the number of times when the first sample exceeds the second sample. The null hypothesis test indicates there are no changes in the data series, while the alternative hypothesis indicates the existence of changing points in the data series. The commonly used 5% significant level was used for the analysis. The test statistic K n and the associate probability P used in the test can be obtained using the following Equations (11) and (12), respectively:
K n = m a x 1 t n | U t ,   n |
P 2 e x p { 6 ( K n ) 2 ( n 3 + n 2 ) }

2.7. Indicators of Hydrological Alteration (IHA)

IHA software was developed by the US Nature Conservancy to measure the extent of hydrological changes due to climatic and human influences [21]. This tool has the ability to calculate 33 IHA parameters under five categories: magnitudes of monthly water conditions, magnitude and duration of annual extreme water conditions, the timing of annual extreme water conditions, frequency and duration of high/low pulses, and rate/ frequency of water condition changes. A 5% significance level p-value was used to evaluate and compare the parameters’ consistency.

2.8. RClimdex

RClimdex software was developed and maintained by Xuebin Zhang and Feng Yang at the Climate Research Branch of the Meteorological Service of Canada and designed to obtain 27 climate extreme indices that are recommended by the Expert Team for Climate Change Detection Monitoring and Indices (ETCCDMI) [4]. This study used five selected extreme precipitation indices: consecutive dry days (CDD), consecutive wet days (CWD), annual total wet day precipitation (PRCPTOT), number of days above 25 mm (R25), and maximum 5-day precipitation amount (Rx5) [45].

3. Results and Discussion

3.1. Correlation between Observed and APHRODITE Data

Since rainfall was missing from some of the rainfall stations, gridded precipitation data were used to fill them. To check the reliability, the correlation coefficient was calculated between the observed data and APHRODITE products V1901 (from January 2004 to March 2015) and V1101 (from February 1991 to September 1991). The longest and continuous data periods that were common for both the observed and APHRODITE data were selected on this basis. The results of the correlation are shown in Table 2 below. Most of the stations showed a moderate-to-very-strong relationship in the correlation analysis. Hence, we could justify the use of APHRODITE data for our study. The types of correlations were classified as per the rule of thumb for interpreting the correlation coefficient. The classification for the Pearson correlation coefficient was as follows: a positive very strong correlation was 0.8–1.0, a positive strong correlation was 0.6–0.8, a positive strong moderate correlation was 0.4–0.6, a positive moderate correlation was 0.2–0.4, and a positive weak or insignificant correlation was 0–0.2 [3].

3.2. Trend Analysis of the Rainfall

Trend analysis was carried out for the monthly, seasonal, and annual scales using Mann–Kendall and Sen’s slope tests. The significant trends are denoted in bold font in Table 3. Dampahala station showed significant increasing trends in March, September, November, and December, with magnitudes of 12.2 mm/yr, 8.96 mm/yr, 14.26 mm/yr, and 11.65 mm/yr, respectively. Kamburupitiya showed significant decreasing trends of 3.97 mm/yr in July. Kirama demonstrated significant decreasing trends in April, June, and July, with magnitudes of 7.23 mm/yr, 4.63 mm/yr, and 4.61 mm/yr, respectively. Deniyaya also revealed a significant 13.83 mm/yr decreasing trend in May, while Kekenadura and Goluwatta did not show any significant trends at the monthly scale. In general, most of the stations showed decreasing patterns at the monthly scale. Interestingly, Deniyaya showed only decreasing trends and Dampahala showed only increasing trends. Dampahala showed a significant increasing trend of 73.85 mm/yr, while Deniyaya showed a significant decreasing trend of 70.3 mm/yr at the annual scale. At the seasonal scale, Dampahala station revealed significant increasing trends during the NEM and FIM, with magnitudes of 23.43 mm/yr and 17.42 mm/yr, respectively. During the SWM, Kamburupitiya, Kirama, and Deniyaya showed significant 16.83 mm/yr, 11.68 mm/yr, and 33.19 mm/yr decreasing trends, respectively. During the SIM period, only the Deniyaya station revealed a significant decreasing trend, with a magnitude of 13.89 mm/yr. The highest magnitudes of increasing trends at the monthly (November) and annual scales and for the NEM and FIM were identified in the Dampahala station. Deniyaya exhibited the highest significant decreasing trends at the monthly (May) and annual scales and for the SWM and SIM. Other stations, except for Dampahala, experienced decreasing trends at the monthly and annual scales and for the SWM and FIM.

3.3. Trend Analysis of the Streamflow

Trend analysis of the streamflow at the Pitabeddara hydrologic station was computed for the monthly, seasonal, and annual scales. According to the results in Table 4, only one significant trend was observed, which was in December with a magnitude of 0.59 m3s−1/yr. During the other months, non-significant increasing and decreasing trends were observed. Significant trends were not observed on the annual or seasonal scale. However, non-significant increasing trends were observed for the annual scale and the NEM and FIM. Non-significant decreasing trends during the SWM and SIM were seen.

3.4. Trend Analysis of the Extreme Rainfall Indices

Five selected extreme rainfall indices, namely, consecutive dry days (CDD), consecutive wet days (CWD), annual total wet day precipitation (PRCPTOT), number of days above 25 mm (R25), and maximum 5-day precipitation amount (Rx5), were computed using the RClimdex software. According to the results shown in Table 5, Dampahala station revealed significant increasing trends in PRCPTOT and R25, with magnitudes of 74.55 mm/yr and 1.64 days/yr. In Deniyaya, significant decreasing trends in PRCPTOT and R25 were observed, with magnitudes of 71.75 mm/yr and 1.33 days/yr. Kamburupitiya showed a 0.7 days/yr significant decreasing trend for CWD. No significant trends were observed for the Kekenadura, Kirama, and Goluwatta stations for CDD, CWD, PRCPTOT, and R25.
Table 6 and Table 7 present the trend results of Rx5. According to the trend results of Rx5, significant increasing trends in March, September, and December were found for the Dampahala station, with values of 5.76 mm/yr, 4.92 mm/yr, and 5.25 mm/yr. Kamburupitiya revealed a 3.38 mm/yr significant increasing trend in December, while in the annual scale analysis, Kamburupitya revealed a significant decreasing trend, with a magnitude of 3.63 mm/yr. Kirama and Goluwatta stations also revealed decreasing trends in July and May, with magnitudes of 3.34 mm/yr and 4.0 mm/yr, respectively. Deniyaya revealed a comparatively high number of significant decreasing trend events in February, April, May, June, and November, and at the annual scale, with magnitudes of 4.24 mm/yr, 3.13 mm/yr, 5.24 mm/yr, 4,88 mm/yr, 4.81 mm/yr, and 6.45 mm/yr respectively. Generally, Rx5 showed more significant decreasing trends in both monthly and annual scales at most stations.

3.5. Change Point Detection in the Rainfall Data

Pettitt’s test results, which showed statistically significant changes at the annual and seasonal scales, are shown in Figure 3. Most stations did not show significant increasing or decreasing changes at both scales. The annual scale results of the Pettit’s test at Dampahala revealed a significant increasing shift in the year 1998 and Deniyaya revealed a significant decreasing shift in the year 2000. However, the Dampahala station showed significant increasing shifts in 1998 during the NEM, FIM, and SIM. Moreover, Kamburupitiya and Deniyaya showed significant decreasing changes in 1998 and 1999, respectively.
Table 8 shows a few significant changes in the monthly Pettitt’s test results for the Dampahala, Kamburupitiya, Kirama, and Deniyaya rainfall stations. The Dampahala station revealed significant increasing changes during March, September, and December in 1999, 1999, and 1997, respectively. The Kamburupitiya station revealed a significant decreasing shift in July 1998. Moreover, the Kirama station showed significant decreasing shifts in June 2000 and July 2005, and the Deniyaya station showed a significant decreasing change in October 1997.

3.6. Linkage between Rainfall and Streamflow

Pearson’s correlation was used to analyze the relationship between rainfall and streamflow during 1991–2014 at the monthly, annual, and seasonal timescales. Considering the contribution of rainfall to the Pitabeddara streamflow station, only the Dampahala, Kirama, and Deniyaya stations were chosen to find the existing linkages between rainfall and streamflow. The results of the correlation analysis are given in Table 9.
The Pearson correlation results indicated that February, March, May, and September showed positive very strong correlations, while January, June, November, and December showed positive strong correlations between rainfall and streamflow. Moreover, at the monthly scale, April, July, and October showed positive strong moderate correlations between rainfall and streamflow, while August had a positive moderate correlation. The annual scale results showed a positive strong moderate correlation between rainfall and streamflow. At the seasonal scales, a positive very strong correlation was observed for the NEM, and the FIM showed a positive strong moderate correlation. Both the SWM and SIM periods showed positive strong correlations between the rainfall and streamflow data.

3.7. Indicators of Hydrological Alteration (IHA)

According to the literature based on the NRB, no major obstruction or dam was built. However, we decided to investigate the changes in the flow regime before and after 2003. This year was selected purely arbitrarily. Two-period parametric analysis was used in this IHA method. The pre-impact period was chosen from 1991–2003 and the post-impact period was selected from 2004–2014. This section discusses 16 selected hydrological parameters out of the 33 that fall under the main IHA parameter groups of magnitude and duration of annual extreme water conditions, timing of annual extreme water conditions, and frequency and duration of high and low pulses. According to the IHA analysis, the mean annual flow during the post-impact period increased slightly from 15.87 m3s−1 to 16.23 m3s−1. Table 10 given below demonstrates the statistics for the pre-impact and post-impact periods for the selected IHA parameters.
Although notable changes were not observed for minimum flows in IHA group 2, maximum flows decreased significantly during the post-impact period. For instance, the 1-day maximum, 3-day maximum, and 7-day maximum flows decreased by 27.94%, 30.82%, and 23.87%, respectively. The only notable change in minimum flow was the 90-day minimum flow with an increase of 18.26%.
Results of the IHA group 3 parameters in Table 9 showed that the annual 1-day minimum flow was recorded on the 113th Julian date in the calendar for the pre-impact period and the 119th and 120th days for the post-impact period; this showed that the date shifted a little bit forward in the post-impact period. Julian’s date of each annual 1-day maximum for the pre-impact was recorded during the 211th and 212th days, and for the post-impact period, it was recorded during the 261st and 262nd days.
According to the results in Table 10 under the IHA group 3 category, 14.56 low pulses were found in the pre-impact period, which decreased to 13.36 per year during the post-impact period. The coefficient of variation for the low pulses was increased by 69.54% from the pre-impact to post-impact period. The number of high pulses for each year in the pre-impact period was found to be 12.54 per year, which increased in the post-impact period up to 13.73 per year. The deviation of the coefficient of variation also decreased by 27.29% from the pre-impact to post-impact period.

3.8. Discussion

According to the results found in the present study, most of the stations showed decreasing rainfall trend patterns at the monthly scale. Generally, the annual scale and the FIM and SWM exhibited decreasing rainfall trend patterns, NEM exhibited an increasing tendency and SIM exhibited average results. A previous study that was carried out during 1987–2017 by Nisansala et al. [46] reported similar results for the NEM, SWM, and SIM seasons. Wickramagamage [47] also reported an increasing rainfall trend during the NEM season and decreasing trends in the SWM for Sri Lanka during 1981–2010. However, Nisansala et al. [46] and Wickramagamage [47] also showed different results compared with the present study at the annual and seasonal scales for rainfall trends by showing an increasing tendency. In the present study, the Dampahala station showed increasing trends for rainfall, while all the other stations showed decreasing trend patterns. Similar to the present study, contrasting directions of the magnitude of trends in nearby stations of the same basin were reported by Khaniya et al. [48] and Pawar and Rathnayake [49]. These contrasting results might have been because of local rainfall events. Mehta and Yadav [50] demonstrated that the magnitude of climate variability change across spatial scales. According to the results of the present study, the streamflow trend analysis did not show significant trends, except in December. However, non-significant increasing trends were demonstrated at the annual scale and the NEM and FIM, but not for the SWM SIM. Dinethra and Basnayake [51] showed that the discharge of the Nilwala River increased during 2004–2013. When considering the extreme rainfall indices, they also showed few significant trends. However, extreme rainfall indices trend patterns also showed similar contrasting results, especially for the Dampahala and Deniyaya stations. These contrasting results might have been due to variations in elevations, as explained by Bizuneh [52]. The NRB comprises lowlands and mountains. Due to the windward and leeward sides of the mountains, these types of contrasting results in the direction of rainfall trends can happen. This is because the windward side normally receives higher rainfall, while the leeward side of the mountain gets lower rainfall. Considering rainfall, extreme rainfall, and streamflow trend results, we concluded that rainfall was not only the influencing factor for the changes in streamflow patterns for the NRB. Other factors, such as the density of physical features, watershed characteristics, and vegetation cover, can be influential as well [2].

4. Conclusions

The main objective of this study was to analyze the trends that were present in the rainfall and streamflow records and to identify the correlation between the rainfall and streamflow over the Nilwala River Basin. Six rainfall stations and one hydrological gauge station were chosen based on the data availability. There was a considerably good correlation present between the rainfall and streamflow in the upper Nilwala Basin, indicating that rainfall was the main driver for the changes in the streamflow. This study showed that the NRB did not face extreme climatic events from 1990 to 2014. Variations in the topographical features in the basin might lead to contrasting results in rainfall trends within the NRB. Five out of six stations showed decreasing trends in rainfall, suggesting that, in general, rainfall had been decreasing during the 25 years between 1990 and 2014. The insignificant trends in rainfall and streamflow suggested that the climate and hydrologic regimes were not altered during this period. Future research is advocated for in the NRB considering that other possible factors, such as the impacts of land-use change, could also be influential.

Author Contributions

Conceptualization, M.B.G.; methodology, R.P. and M.B.G.; software, R.P.; formal analysis, R.P.; writing—original draft preparation, R.P.; writing—review and editing, M.B.G., R.P., I.M.C., M.K.J., M.S.B. and U.R.; supervision, M.B.G., M.S.B., U.R. and M.K.J.; funding acquisition, M.B.G. and U.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research was carried out under the financial support of a SLIIT Research Grant: FGSR/RG/FE/2022/02.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The climatic data used in this research study are available upon request for research purposes.

Acknowledgments

The authors of this study would like to thank the Sri Lanka Institute of Information Technology (SLIIT), Wayamba University of Sri Lanka, and the Norwegian Institute of Bioeconomy Research (NIBIO) for providing the environment to conduct this research.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fentaw, F.; Hailu, D.; Nigussie, A. Trend and Variability Analysis of Rainfall & Stream Flow Series at Tekeze River Basin, Ethiopia. Int. J. Sci. Eng. Res. 2017, 8, 665–680. [Google Scholar]
  2. Ademe, D.; Alamirew, T.; Rotich, J.; Thi, N.; Linh, T.; Gebrie, T. Analysis of Rainfall and Streamflow Trend and Variability over Birr River Watershed, Abbay Basin, Ethiopia. Environ. Chall. 2022, 7, 100528. [Google Scholar] [CrossRef]
  3. Chauluka, F.; Singh, S.; Kumar, R. Rainfall and Streamflow Trends of Thuchila River, Southern Malawi. Mater. Today Proc. 2021, 34, 846–855. [Google Scholar] [CrossRef]
  4. Zhai, R.; Tao, F. Contributions of Climate Change and Human Activities to Runoff Change in Seven Typical Catchments across China. Sci. Total Environ. 2017, 605–606, 219–229. [Google Scholar] [CrossRef] [PubMed]
  5. Alemu, M.M.; Bawoke, G.T. Analysis of Spatial Variability and Temporal Trends of Rainfall in Amhara Region, Ethiopia. J. Water Clim. Chang. 2020, 11, 1505–1520. [Google Scholar] [CrossRef]
  6. Kliment, Z.; Matoušková, M. Long-Term Trends of Rainfall and Runoff Regime in Upper Otava River Basin. Soil Water Res. 2008, 3, 155–167. [Google Scholar] [CrossRef] [Green Version]
  7. Arrieta-Castro, M.; Donado-Rodríguez, A.; Acuña, G.J.; Canales, F.A.; Teegavarapu, R.S.V.; Kaźmierczak, B. Analysis of Streamflow Variability and Trends in the Meta River, Colombia. Water 2020, 12, 1451. [Google Scholar] [CrossRef]
  8. Charifi Bellabas, S.; Benmamar, S.; Dehni, A. Study and Analysis of the Streamflow Decline in North Algeria. J. Appl. Water Eng. Res. 2021, 9, 20–44. [Google Scholar] [CrossRef]
  9. Fentaw, F.; Melesse, A.M.; Hailu, D.; Nigussie, A. Precipitation and Streamflow Variability in Tekeze River Basin, Ethiopia. In Extreme Hydrology and Climate Variability; Elsevier: Amsterdam, The Netherlands, 2019; pp. 103–121. [Google Scholar] [CrossRef]
  10. Mersin, D.; Tayfur, G.; Vaheddoost, B.; Safari, M.J.S. Historical Trends Associated with Annual Temperature and Precipitation in Aegean Turkey, Where Are We Heading? Sustainability 2022, 14, 3380. [Google Scholar] [CrossRef]
  11. Kastridis, A.; Kamperidou, V.; Stathis, D. Dendroclimatological Analysis of Fir (A. Borisii-Regis) in Greece in the Frame of Climate Change Investigation. Forests 2022, 13, 879. [Google Scholar] [CrossRef]
  12. Hannaford, J. Climate-Driven Changes in UK River Flows: A Review of the Evidence. Prog. Phys. Geogr. 2015, 39, 29–48. [Google Scholar] [CrossRef] [Green Version]
  13. Wang, G.; Xia, J.; Chen, J. Quantification of Effects of Climate Variations and Human Activities on Runoff by a Monthly Water Balance Model: A Case Study of the Chaobai River Basin in Northern China. Water Resour. Res. 2009, 45, 1–12. [Google Scholar] [CrossRef] [Green Version]
  14. Azari, M.; Moradi, H.R.; Saghafian, B.; Faramarzi, M. Climate Change Impacts on Streamflow and Sediment Yield in the North of Iran. Hydrol. Sci. J. 2016, 61, 123–133. [Google Scholar] [CrossRef]
  15. Dey, P.; Mishra, A. Separating the Impacts of Climate Change and Human Activities on Streamflow: A Review of Methodologies and Critical Assumptions. J. Hydrol. 2017, 548, 278–290. [Google Scholar] [CrossRef]
  16. Xu, F.; Zhao, L.; Niu, C.; Qiu, Y. Effect of Climate Change and Anthropogenic Activities on Streamflow Indicators in a Tropical River Basin in Southern China. Water 2022, 14, 304. [Google Scholar] [CrossRef]
  17. Jayasekara, S.M.; Abeysingha, N.S.; Meegastenna, T. Streamflow Trends of Kelani River Basin in Sri Lanka (1983–2013). J. Natl. Sci. Found. Sri Lanka 2020, 48, 449. [Google Scholar] [CrossRef]
  18. Kale, S.; Ejder, T.; Hisar, O.; Mutlu, F. Climate Change Impacts on Streamflow of Karamenderes River (Çanakkale, Turkey). Mar. Sci. Technol. Bull. 2016, 5, 1–6. [Google Scholar] [CrossRef]
  19. Abungba, J.A.; Khare, D.; Pingale, S.M.; Adjei, K.A.; Gyamfi, C.; Odai, S.N. Assessment of Hydro-Climatic Trends and Variability over the Black Volta Basin in Ghana. Earth Syst. Environ. 2020, 4, 739–755. [Google Scholar] [CrossRef]
  20. Salarijazi, M. Trend and Change-Point Detection for the Annual Stream-Flow Series of the Karun River at the Ahvaz Hydrometric Station. Afr. J. Agric. Res. 2012, 7, 4540–4552. [Google Scholar] [CrossRef] [Green Version]
  21. Zhang, W.; Yan, Y.; Zheng, J.; Li, L.; Dong, X.; Cai, H. Temporal and Spatial Variability of Annual Extreme Water Level in the Pearl River Delta Region, China. Glob. Planet. Chang. 2009, 69, 35–47. [Google Scholar] [CrossRef]
  22. Caloiero, T.; Coscarelli, R.; Ferrari, E. Application of the Innovative Trend Analysis Method for the Trend Analysis of Rainfall Anomalies in Southern Italy. Water Resour. Manag. 2018, 32, 4971–4983. [Google Scholar] [CrossRef]
  23. Mathews, R.; Richter, B.D. Application of the Indicators of Hydrologic Alteration Software in Environmental Flow Setting1: Tools for Environmental Flow Setting. JAWRA J. Am. Water Resour. Assoc. 2007, 43, 1400–1413. [Google Scholar] [CrossRef]
  24. Barbalić, D.; Kuspilić, N. Trends of Indicators of Hydrological Alterations. GRADEVINAR 2014, 66, 613–624. [Google Scholar] [CrossRef] [Green Version]
  25. Gao, Y.; Vogel, R.M.; Kroll, C.N.; Poff, N.L.; Olden, J.D. Development of Representative Indicators of Hydrologic Alteration. J. Hydrol. 2009, 374, 136–147. [Google Scholar] [CrossRef]
  26. Chathuranika, I.M.; Gunathilake, M.B.; Azamathulla, H.M.; Rathnayake, U. Evaluation of Future Streamflow in the Upper Part of the Nilwala River Basin (Sri Lanka) under Climate Change. Hydrology 2022, 9, 48. [Google Scholar] [CrossRef]
  27. Abeysingha, N.S. A Review of Recent Changes in Rainfall Trend in Sri Lanka. Trop. Agric. Res. Ext. 2022, 25, 1–13. [Google Scholar] [CrossRef]
  28. Perera, A.; Ranasinghe, T.; Gunathilake, M.; Rathnayake, U. Comparison of Different Analyzing Techniques in Identifying Rainfall Trends for Colombo, Sri Lanka. Adv. Meteorol. 2020, 2020, 8844052. [Google Scholar] [CrossRef]
  29. Alahacoon, N.; Edirisinghe, M. Spatial Variability of Rainfall Trends in Sri Lanka from 1989 to 2019 as an Indication of Climate Change. Int. J. Geo-Inf. 2021, 10, 84. [Google Scholar] [CrossRef]
  30. Ruwangika, A.M.; Perera, A.; Rathnayake, U. Comparison of Statistical, Graphical, and Wavelet Transform Analyses for Rainfall Trends and Patterns in Badulu Oya Catchment, Sri Lanka. Complexity 2020, 2020, 7146593. [Google Scholar] [CrossRef]
  31. Dhanapala, L.; Gunarathna, M.H.J.P.; Kumari, M.K.N.; Ranagalage, M.; Sakai, K.; Meegastenna, T.J. Towards Coupling of 1D and 2D Models for Flood Simulation—A Case Study of Nilwala River Basin, Sri Lanka. Hydrology 2022, 9, 17. [Google Scholar] [CrossRef]
  32. Yatagai, A.; Arakawa, O.; Kamiguchi, K.; Kawamoto, H.; Nodzu, M.I.; Hamada, A. A 44-Year Daily Gridded Precipitation Dataset for Asia Based on a Dense Network of Rain Gauges. SOLA 2009, 5, 137–140. [Google Scholar] [CrossRef] [Green Version]
  33. Yatagai, A.; Kamiguchi, K.; Arakawa, O.; Hamada, A.; Yasutomi, N.; Kitoh, A. APHRODITE: Constructing a Long-Term Daily Gridded Precipitation Dataset for Asia Based on a Dense Network of Rain Gauges. Bull. Am. Meteorol. Soc. 2012, 93, 1401–1415. [Google Scholar] [CrossRef]
  34. Marie, M.; Yirga, F.; Haile, M.; Ehteshammajd, S.; Azadi, H.; Scheffran, J. Time-Series Trend Analysis and Farmer Perceptions of Rainfall and Temperature in Northwestern Ethiopia. Environ. Dev. Sustain. 2021, 23, 12904–12924. [Google Scholar] [CrossRef]
  35. Mutiibwa, D.; Vavrus, S.J.; McAfee, S.A.; Albright, T.P. Recent Spatiotemporal Patterns in Temperature Extremes across Conterminous United States. J. Geophys. Res. Atmos. 2015, 120, 7378–7392. [Google Scholar] [CrossRef]
  36. Mann, H.B. Nonparametric Tests Against Trend. Econometrica 1945, 13, 245–259. [Google Scholar] [CrossRef]
  37. Kendall, M.G. Rank Correlation Methods, 4th ed.; Griffin: London, UK, 1970. [Google Scholar]
  38. Yue, S.; Pilon, P.; Cavadias, G. Power of the Mann–Kendall and Spearman’s Rho Tests for Detecting Monotonic Trends in Hydrological Series. J. Hydrol. 2002, 259, 254–271. [Google Scholar] [CrossRef]
  39. Zhang, X.; Harvey, K.D.; Hogg, W.D.; Yuzyk, T.R. Trends in Canadian Streamflow. Water Resour. Res. 2001, 37, 987–998. [Google Scholar] [CrossRef]
  40. Birsan, M.-V.; Molnar, P.; Burlando, P.; Pfaundler, M. Streamflow Trends in Switzerland. J. Hydrol. 2005, 314, 312–329. [Google Scholar] [CrossRef]
  41. Kulkarni, A.; Storch, H.V. Monte Carlo Experiments on the Effect of Serial Correlation on the Mann-Kendall Test of Trend. Meteorol. Z. 1995, 4, 82–85. [Google Scholar] [CrossRef]
  42. Hamed, K.H.; Rao, A.R. A Modified Mann-Kendall Trend Test for Autocorrelated Data. J. Hydrol. 1998, 204, 182–196. [Google Scholar] [CrossRef]
  43. Sen, P.K. Estimates of the Regression Coefficient Based on Kendall’s Tau. J. Am. Stat. Assoc. 1968, 63, 1379–1389. [Google Scholar] [CrossRef]
  44. Pettitt, A.N. A Non-Parametric Approach to the Change-Point Problem. Appl. Stat. 1979, 28, 126. [Google Scholar] [CrossRef]
  45. Zhang, X.; Yang, F. RClimDex (1.0) User Manual; Climate Research Branch Environment Canada: Downsview, ON, Canada, 2004. [Google Scholar]
  46. Nisansala, W.D.S.; Abeysingha, N.S.; Islam, A.; Bandara, A.M.K.R. Recent Rainfall Trend over Sri Lanka (1987–2017). Int. J. Climatol. 2020, 40, 3417–3435. [Google Scholar] [CrossRef]
  47. Wickramagamage, P. Spatial and Temporal Variation of Rainfall Trends of Sri Lanka. Theor. Appl. Climatol. 2016, 125, 427–438. [Google Scholar] [CrossRef]
  48. Khaniya, B.; Jayanayaka, I.; Jayasanka, P.; Rathnayake, U. Rainfall Trend Analysis in Uma Oya Basin, Sri Lanka, and Future Water Scarcity Problems in Perspective of Climate Variability. Adv. Meteorol. 2019, 2019, 3636158. [Google Scholar] [CrossRef] [Green Version]
  49. Pawar, U.; Rathnayake, U. Spatiotemporal Rainfall Variability and Trend Analysis over Mahaweli Basin, Sri Lanka. Arab J. Geosci. 2022, 15, 370. [Google Scholar] [CrossRef]
  50. Mehta, D.; Yadav, S.M. An Analysis of Rainfall Variability and Drought over Barmer District of Rajasthan, Northwest India. Water Supply 2021, 21, 2505–2517. [Google Scholar] [CrossRef]
  51. Dinethra, A.G.T.; Basnayake, D.M.L.A. Contribution of climate change and human activities on the runoff changes in the upper catchment of Nilwala basin, Sri Lanka. In Proceedings of the 20th IAHR APD Congress, Colombo, Sri Lanka, 28–31 August 2016. [Google Scholar]
  52. Mekuriaw Bizuneh, A. Climate Variability and Change in the Rift Valley and Blue Nile Basin, Ethiopia: Local Knowledge, Impacts, and Adaptation; UAMR Studies on Development and Global Governance; Logos Verlag Berlin GmbH: Berlin, Germany, 2013. [Google Scholar]
Figure 1. Location map of Nilwala River Basin with rainfall and hydrological stations.
Figure 1. Location map of Nilwala River Basin with rainfall and hydrological stations.
Hydrology 10 00008 g001
Figure 2. Overall methodology.
Figure 2. Overall methodology.
Hydrology 10 00008 g002
Figure 3. Pettitt’s test results of the annual and seasonal rainfall data: (a) Deniyaya annual, (b) Dampahala annual, (c) Dampahala SIM, (d) Dampahala NEM, (e) Dampahala FIM, (f) Kamburupitiya SWM, and (g) Deniyaya SWM.
Figure 3. Pettitt’s test results of the annual and seasonal rainfall data: (a) Deniyaya annual, (b) Dampahala annual, (c) Dampahala SIM, (d) Dampahala NEM, (e) Dampahala FIM, (f) Kamburupitiya SWM, and (g) Deniyaya SWM.
Hydrology 10 00008 g003
Table 1. General information of the stations used.
Table 1. General information of the stations used.
StationLatitude (N)Longitude (E)Elevation (m.a.s.l)PeriodMissing %
Meteorological stations
Dampahala6.2780.64176

1991–2014





1991–2014
25.03
Kamburupitiya6.0880.562445.91
Kekenadura5.9780.57493.49
Kirama6.2280.671223.79
Goluwawatta6.1080.481617.10
Deniyaya6.3380.5539924.66
Hydrological station
Pitabaddara6.2080.48270.205
Table 2. Results of correlation check—observed data vs. APHRODITE data.
Table 2. Results of correlation check—observed data vs. APHRODITE data.
StationV1901 (January 2004–March 2005)V1101 (February 1991–September 1991)
CorrelationCorrelation TypeCorrelationCorrelation Type
Dampahala0.47Positive moderate0.12Positive weak
Kamburupitiya0.79Positive very strong0.69Positive strong
Kekenadura0.85Positive very strong0.61Positive strong
Kirama0.28Positive weak0.02Positive weak
Goluwawatta0.77Positive very strong0.38Positive moderate
Deniyaya0.45Positive moderate0.75Positive very strong
Table 3. Trend analysis results for the rainfall data.
Table 3. Trend analysis results for the rainfall data.
TimescaleDampahalaKamburupitiyaKekenaduraKiramaGoluwattaDeniyaya
p-ValueSen’s SlopeTrend Typep-ValueSen’s SlopeTrend Typep-ValueSen’s SlopeTrend Typep-ValueSen’s SlopeTrend Typep-ValueSen’s SlopeTrend Typep-ValueSen’s SlopeTrend Type
January0.1882.88NSIT0.9410.11NSIT0.7470.53NSIT0.5511.39NSIT0.9410.37NSIT0.637−1.33NSDT
February0.1087.08NSIT0.823−0.40NSDT0.7470.89NSIT0.823−0.19NSDT0.9800.49NSIT0.063−9.49NSDT
March0.00512.20SIT0.2481.72NSIT0.1361.57NSIT0.2052.55NSIT0.9410.50NSIT0.333−2.78NSDT
April0.1328.46NSIT0.248−3.93NSDT0.209−2.18NSDT0.039−7.23SDT0.713−1.25NSDT0.188−10.12NSDT
May0.3384.86NSIT0.268−3.57NSDT0.447−3.08NSDT0.980−0.20NSDT0.160−4.61NSDT0.039−13.83SDT
June0.6771.59NSIT0.063−4.45NSDT0.7880.91NSIT0.050−4.63SDT0.864−0.44NSDT0.118−8.26NSDT
July0.6771.24NSIT0.044−3.97SDT0.175−2.16NSDT0.002−4.61SDT0.192−3.38NSDT0.078−3.01NSDT
August0.3903.77NSIT0.7500.69NSIT0.2442.73NSIT0.785−0.42NSDT0.5401.51NSIT0.535−2.22NSDT
September0.0208.96SIT0.192−4.33NSDT0.573−1.31NSDT0.338−2.02NSDT0.192−5.17NSDT0.823−0.72NSDT
October0.3642.62NSIT0.9030.69NSIT0.6022.44NSIT0.9030.66NSIT0.607−2.09NSDT0.078−7.56NSDT
November0.01714.26SIT0.826−0.62NSDT0.677−2.28NSDT0.3383.55NSIT0.418−2.58NSDT0.188−4.32NSDT
December0.00311.65SIT0.1205.60NSIT0.1084.66NSIT0.3383.49NSIT0.5731.50NSIT0.980−0.06NSDT
Annual0.03473.85SIT0.314−13.94NSDT0.7135.03NSIT0.677−8.62NSDT0.228−17.06NSDT0.002−70.30SDT
NEM0.01323.43SIT0.2686.49NSIT0.1454.05NSIT0.7504.93NSIT0.750−1.81NSDT0.160−15.75NSDT
FIM0.01517.42SIT0.477−2.96NSDT0.826−0.74NSDT0.314−6.48NSDT0.747−2.10NSDT0.056−17.21NSDT
SWM0.26817.64NSIT0.007−16.83SDT0.338−8.45NSDT0.050−11.68SDT0.088−14.58NSDT0.003−33.19SDT
SIM0.06316.44NSIT0.941−0.40NSDT0.6373.22NSIT0.4773.47NSIT0.447−6.36NSDT0.031−13.98SDT
* NSIT refers to a non-significant increasing trend, NSDT refers to a non-significant decreasing trend, SIT refers to a significant increasing trend, SDT refers to a significant decreasing trend.
Table 4. Trend analysis results for the streamflow data.
Table 4. Trend analysis results for the streamflow data.
TimescaleKendall’s Taup-ValueSen’s SlopeTrend Type
January0.1520.2980.17NSIT
February0.1230.3990.13NSIT
March0.2170.1370.25NSIT
April0.1230.3990.15NSIT
May−0.0580.691−0.15NSDT
June−0.0800.567−0.15NSDT
July−0.1300.372−0.18NSDT
August−0.0510.728−0.02NSDT
September−0.0140.921−0.10NSDT
October−0.1960.067−0.30NSDT
November0.1010.4870.26NSIT
December0.3190.0290.59SIT
Annual0.0940.5190.10NSIT
NEM0.2750.0590.31NSIT
FIM0.2030.1650.20NSIT
SWM−0.1590.275−0.12NSDT
SIM−0.0360.804−0.03NSDT
* NSIT refers to a non-significant increasing trend, NSDT refers to a non-significant decreasing trend, SIT refers to a significant increasing trend, SDT refers to a significant decreasing trend.
Table 5. Trend analysis results for the extreme rainfall indices.
Table 5. Trend analysis results for the extreme rainfall indices.
StationExtremesKendall’s Taup-ValueSen’s SlopeTrend Type
DampahalaCDD−0.1150.440−0.14NSDT
CWD−0.2680.070−0.38NSDT
PRCPTOT0.3120.03474.55SIT
R250.3960.0071.64SIT
KamburupitiyaCDD0.1290.3840.20NSIT
CWD−0.4620.002−0.70SDT
PRCPTOT−0.1450.338−14.53NSDT
R250.1060.4710.23NSIT
KekenaduraCDD−0.1720.243−0.31NSDT
CWD−0.0610.6880.00No trend
PRCPTOT0.0650.6775.60NSIT
R250.0300.8420.00No trend
KiramaCDD0.1280.3840.57NSIT
CWD−0.2820.065−0.12NSDT
PRCPTOT−0.0580.713−8.45NSDT
R25−0.0220.8810.00No trend
GoluwattaCDD−0.0690.637−0.10NSDT
CWD−0.1340.369−0.15NSDT
PRCPTOT−0.1740.248−18.10NSDT
R25−0.2740.065−0.60NSDT
DeniyayaCDD0.1010.5000.09NSIT
CWD−0.0440.765−0.04NSDT
PRCPTOT−0.4540.002−71.75SDT
R25−0.4420.003−1.33SDT
* NSIT refers to a non-significant increasing trend, NSDT refers to a non-significant decreasing trend, SIT refers to a significant increasing trend, SDT refers to a significant decreasing trend.
Table 6. Trend analysis results for the extreme rainfall indices (Rx5).
Table 6. Trend analysis results for the extreme rainfall indices (Rx5).
TimescaleDampahalaKamburupitiyaKekenadura
p-ValueSen’s SlopeTrend Typep-ValueSen’s SlopeTrend Typep-ValueSen’s SlopeTrend Type
January0.0822.64NSIT0.6420.56NSIT0.6070.54NSIT
February0.3142.48NSIT0.9030.12NSIT0.2151.49NSIT
March0.0025.76SIT0.7100.39NSIT0.2281.10NSIT
April0.1324.51NSIT0.096−2.56NSDT0.508−0.81NDIT
May0.0963.05NSIT0.268−1.99NSDT0.145−2.45NDIT
June0.3381.88NSIT0.143−2.58NSDT0.8640.27NSIT
July0.6070.70NSIT0.447−0.44NSDT0.079−1.46NDIT
August0.4421.38NSIT0.6070.75NSIT0.1241.56NSIT
September0.0084.92SIT0.447−1.43NSDT0.359−1.78NDIT
October0.0702.56NSIT0.710−0.35NSDT0.785−0.22NDIT
November0.0783.87NSIT0.573−0.59NSDT0.9210.07NSIT
December0.0075.25SIT0.0303.38SIT0.1922.16NSIT
Annual0.2444.32NSIT0.030−3.63SDT0.228−1.11NDIT
* NSIT refers to a non-significant increasing trend, NSDT refers to a non-significant decreasing trend, SIT refers to a significant increasing trend, SDT refers to a significant decreasing trend.
Table 7. Trend analysis results for the extreme rainfall indices (Rx5).
Table 7. Trend analysis results for the extreme rainfall indices (Rx5).
TimescaleKiramaGoluwattaDeniyaya
p-ValueSen’s SlopeTrend Typep-ValueSen’s SlopeTrend Typep-ValueSen’s SlopeTrend Type
January0.3851.35NSIT0.788−0.25NSDT0.508−1.10NDIT
February0.9010.12NSIT0.750−0.98NSDT0.027−4.24SDT
March0.4130.80NSIT0.980−0.01NSDT0.070−3.19NDIT
April0.160−2.68NSDT0.9800.09NSIT0.050−3.13SDT
May0.9010.27NSIT0.014−4.00SDT0.024−5.24SDT
June0.130−1.62NSDT0.447−1.65NSDT0.050−4.88SDT
July0.003−3.34SDT0.607−0.95NSDT0.172−1.08NDIT
August0.7660.31NSIT0.941−0.05NSDT0.087−2.05NDIT
September0.290−1.57NSDT0.107−3.35NSDT0.286−1.59NDIT
October0.6770.78NSIT0.175−2.68NSDT0.9410.11NSIT
November0.8620.10NSIT0.314−2.32NSDT0.039−4.81SDT
December0.1601.84NSIT0.6730.66NSIT0.8620.62NSIT
Annual0.209−1.41NSDT0.070−5.52NSDT0.014−6.45SDT
* NSIT refers to a non-significant increasing trend, NSDT refers to a non-significant decreasing trend, SIT refers to a significant increasing trend, SDT refers to a significant decreasing trend.
Table 8. Pettitt’s test results of monthly rainfall data.
Table 8. Pettitt’s test results of monthly rainfall data.
StationMonthp-ValueChange PointShift
DampahalaMarch0.0141999Upward
September0.0401999Upward
December0.0201997Upward
KamburupitiyaJuly0.0311998Downward
KiramaJune0.0132000Downward
July0.0072005Downward
DeniyayaOctober0.0341997Downward
Table 9. Pearson correlation results.
Table 9. Pearson correlation results.
TimescaleCorrelationLinkage between Rainfall and Streamflow
January0.79PSC
February0.85PSVC
March0.85PSVC
April0.40PSMC
May0.85PSVC
June0.72PSC
July0.46PSMC
August0.33PMC
September0.80PSVC
October0.56PSMC
November0.61PSC
December0.71PSC
Annual0.55PSMC
NEM0.81PSVC
FIM0.58                     PSMC
SWM0.62                     PSC
SIM0.67                     PSC
PSC—positive strong correlation; PVSC—positive very strong correlation; PSMC—positive strong moderate correlation; PMC—positive moderate correlation.
Table 10. IHA results.
Table 10. IHA results.
IHA ParametersMeansCoefficient of Variation
Pre-ImpactPost-ImpactDeviation Factor (%)Pre-ImpactPost-ImpactDeviation of C.V (%)
Magnitude and Duration (IHA Group 2)
1-day minimum2.3752.322−2.2230.47690.3378−29.17
3-day minimum2.4662.424−1.710.45940.3298−28.22
7-day minimum2.7662.723−1.5640.4490.3073−31.56
30-day minimum4.1714.2241.2660.41710.3226−22.67
90-day minimum6.5647.76218.260.38650.2785−27.96
1-day maximum189.2136.3−27.940.58650.3681−37.24
3-day maximum129.989.85−30.820.73250.4233−42.21
7-day maximum82.7963.03−23.870.59190.334−43.56
30-day maximum43.2740.66−6.0360.31910.358912.48
90-day maximum28.4627−5.1380.24470.2198−10.2
No. of zero days00-00-
Baseflow index0.17220.17370.85780.4010.3362−16.14
Timing (IHA Group 3)
Date of each annual 1-day minimum113199.547.240.18570.222819.96
Date of each annual 1-day maximum211.5261.527.280.26930.2606−3.227
Frequency (IHA Group 4)
No. of low pulses in each year14.5413.36−8.0810.27310.462869.54
No. of high pulses in each year12.5413.739.4810.42530.3092−27.29
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Panditharathne, R.; Gunathilake, M.B.; Chathuranika, I.M.; Rathnayake, U.; Babel, M.S.; Jha, M.K. Trends and Variabilities in Rainfall and Streamflow: A Case Study of the Nilwala River Basin in Sri Lanka. Hydrology 2023, 10, 8. https://doi.org/10.3390/hydrology10010008

AMA Style

Panditharathne R, Gunathilake MB, Chathuranika IM, Rathnayake U, Babel MS, Jha MK. Trends and Variabilities in Rainfall and Streamflow: A Case Study of the Nilwala River Basin in Sri Lanka. Hydrology. 2023; 10(1):8. https://doi.org/10.3390/hydrology10010008

Chicago/Turabian Style

Panditharathne, Ravindu, Miyuru B. Gunathilake, Imiya M. Chathuranika, Upaka Rathnayake, Mukand S. Babel, and Manoj K. Jha. 2023. "Trends and Variabilities in Rainfall and Streamflow: A Case Study of the Nilwala River Basin in Sri Lanka" Hydrology 10, no. 1: 8. https://doi.org/10.3390/hydrology10010008

APA Style

Panditharathne, R., Gunathilake, M. B., Chathuranika, I. M., Rathnayake, U., Babel, M. S., & Jha, M. K. (2023). Trends and Variabilities in Rainfall and Streamflow: A Case Study of the Nilwala River Basin in Sri Lanka. Hydrology, 10(1), 8. https://doi.org/10.3390/hydrology10010008

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop