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.
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 m
3s
−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 m
3s
−1 to 16.23 m
3s
−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].