Spatial Variability of Rainfall Trends in Sri Lanka from 1989 to 2019 as an Indication of Climate Change

Analysis of long-term rainfall trends provides a wealth of information on effective crop planning and water resource management, and a better understanding of climate variability over time. This study reveals the spatial variability of rainfall trends in Sri Lanka from 1989 to 2019 as an indication of climate change. The exclusivity of the study is the use of rainfall data that provide spatial variability instead of the traditional location-based approach. Henceforth, daily rainfall data available at Climate Hazards Group InfraRed Precipitation corrected with stations (CHIRPS) data were used for this study. The geographic information system (GIS) is used to perform spatial data analysis on both vector and raster data. Sen’s slope estimator and the Mann–Kendall (M–K) test are used to investigate the trends in annual and seasonal rainfall throughout all districts and climatic zones of Sri Lanka. The most important thing reflected in this study is that there has been a significant increase in annual rainfall from 1989 to 2019 in all climatic zones (wet, dry, intermediate, and Semi-arid) of Sri Lanka. The maximum increase is recorded in the wet zone and the minimum increase is in the semi-arid zone. There could be an increased risk of floods in the southern and western provinces in the future, whereas areas in the eastern and southeastern districts may face severe droughts during the northeastern monsoon. It is advisable to introduce effective drought and flood management and preparedness measures to reduce the respective hazard risk levels.


Introduction
The global climate system has undergone unprecedented changes due to global warming since the 1950s, and continuous climate change will affect many countries, increasing the risk of flood and drought and other natural disasters, which will finally adversely affect the agricultural economy [1][2][3][4]. Climate change studies have become an essential part of understanding and predicting climate change as it exacerbates the effects of natural disasters such as floods, droughts, heavy rainfall, and temperature [5].
Numerous studies conducted in the recent past have shown that there is a growing rainfall variability globally, regionally, and sub-regionally [6][7][8][9][10][11]. North and Central Asia, the eastern regions of North and South America, and northern Europe are particularly prone to increased rainfall [12][13][14], while Northern and Southern Africa, the Mediterranean, and some other regions receive less rainfall. Although a significant increase in rainfall has been observed in Brazil's northern Amazon region, it has been shown to show a tendency to decrease rainfall as it extends to the southern regions [15]. Considering the context of India, it seems that for different states, the rainfall trend can be identified as both positive and negative [15][16][17][18][19]. In addition, some studies have shown that urban temperature is rising in Sri Lanka [20]. Analysis of rainfall trends has become a research field that has attracted great interest and enthusiasm in recent times [21]. However, as the global or regional scale declines to the national or sub-national level, further localized variations can specific area. In this study, the determination of the validity of CHRPS data as a primary requirement is also introduced.
The study used statistical methods such as the Mann-Kendall (M-K) test and Sen's slope estimator to analyze the rainfall trend at the district and climatic zone level in Sri Lanka [8][9][10][11][12][13][14]. Even though the validity of CHIRPS precipitation data has already been established by many researchers in the global context, which is also applied for the countries in the tropics [49][50][51], an additional cross validation process was perforemed using location-specific rainfall data through the linear regression method.

Study Area
Sri Lanka is an island (Figure 1) which is located below the Indian subcontinent between latitudes 5 • and 10 • N and longitudes 79 • and 82 • E. According to the Department of Census and Statistics, by 2020, Sri Lanka's total population increased to slightly over 21 million, and the total land area is about 65,000 square kilometers [52]. Furthermore, the main administrative units of Sri Lanka comprise nine provinces, 25 districts, and 331 divisional secretariat divisions (DSD). The coastline, coastal plains, inland plains, southcentral mountain ranges, and isolated mountain ranges are the main topographical features of the island. The altitude ranges from 0 to 2505 m, and the highest peak is Piduruthagala. There are 103 rivers in Sri Lanka, covering 103 major and minor river basins.

Data
CHIRPS data were used primarily for this study and are available daily. The daily data were then downloaded from the Climate Hazard Center and these data are available as daily-grided rainfall from 1989 to 2019 for this study [48]. In generating CHIRPS data, precipitation data estimated by satellites with a spatial resolution of 5 km were integrated with location-based rainfall.
Location-specific rainfall data (daily rainfall data) from nine (09) gauge points as a ground-measured source were downloaded from the https://www.irrigation.gov.lk/ (ac- There are four major climatic zones in Sri Lanka: the "wet zone", "intermediate zone", "dry zone" and "semi-arid zone". The country's "wet zone", the western face of the mountain range, and the southwest windy slopes receive an average of 2500 mm of rainfall each year. The "dry zone" of Sri Lanka receives between 1200 and 1900 mm of rainfall, mainly through the northeast monsoon, and the semi-arid zone (northwestern and southeastern coasts) receives the lowest rainfall, which is between 800 and 1200 mm.
The monsoon winds in the Indian Ocean and the depression and cyclones in the Bay of Bengal primarily affect the increase or decrease in rainfall in Sri Lanka [53]. Various studies have introduced four rainy seasons based on rainfall variations in Sri Lanka over a few decades [54,55]. Depending on the months of rain received, the rainy seasons are the first inter-monsoon (FIM) from March to April and the second inter-monsoon (SIM) is from October to November, while December to February constitute the northeastern monsoon (NEM), and the southwestern monsoon (SWM) occurs from May to September.

Data
CHIRPS data were used primarily for this study and are available daily. The daily data were then downloaded from the Climate Hazard Center and these data are available as daily-grided rainfall from 1989 to 2019 for this study [48]. In generating CHIRPS data, precipitation data estimated by satellites with a spatial resolution of 5 km were integrated with location-based rainfall.
Location-specific rainfall data (daily rainfall data) from nine (09) gauge points as a ground-measured source were downloaded from the https://www.irrigation.gov.lk/ (accessed on 31 December 2020) website and converted into monthly accumulation. Data from these nine rainfall stations out of 33 total stations and three years of data (2015 to 2017) were used for the aforementioned cross validation process. The precipitation gauge points used for the analysis were Hanwlla, Glencourse, Deraniyagala, Kitulgala, Holambuwa (Kelani river basin), Pitabaddara, Urawa (Nilwala river basin), Siyanmbalanduva (Heda Oya river), and Tantirimale (Malwatu Oya river), respectively ( Figure 1). The nine rainfall stations mentioned above are those with the highest accuracy, as the rest of the stations have significant data gaps for the considered three-year period. Moreover, it is important to state that these station data, freely available at the Irrigation Department of Sri Lanka, have been available since 2015. On the other hand, Sri Lanka received heavy rainfall as well as severe droughts during the period between 2015 and 2017 [56][57][58]. Therefore, this three-year data analysis represents the validity that CHIRPS data shows during drought as well as heavy rainfall. It is important to mention that through the cross-validation process, we attempt to further emphasize the relevance of CHIRPS data in the Sri Lankan context.

Methodology
Daily rainfall data downloaded from the aforementioned nine rainfall stations were used for the cross-validation process as an additional task which did not affect the findings of this study. Since those station-based rainfall data are available at a daily level, the monthly rainfall was calculated using the daily accumulation of a month. Then, the linear regression was used to observe the correlation between the monthly CHIRPS and the station data.
The downloaded CHIRPS data were in the NetCDF (Network Common Data Format) file format. R-Studio software was used to process these data, as an annual stack of 365 layers. Seasonal and annual rainfall data were generated using CHIRPS daily rainfall accumulation. Since these rainfall data are in raster format, the geo-spatial data analysis method such as cell statistics was used to complete the above process. Subsequently, the seasonal and annual average rainfall values for each year (1989 to 2019) for each district and climate zone used in this study were calculated with the zonal statistics method of spatial statistics. Therefore, the geographical information system (GIS) played a vital role during the spatial data processing. Then, the rainfall trends were analyzed separately at district and climatic zone level. Moreover, the precipitation trend direction (increasing-positive or decreasing-negative) and its magnitude were also calculated at the district and climate zone level using Mann-Kendall's test and Sen's slope estimator, respectively.

Mann-Kendall Test
The Mann-Kendall probability test can be abbreviated as the M-K test. The leastsquare linear regression method is more commonly used than other statistical methods to determine long-term trends for seasonal, annual, and daily maximum rainfall [59]. However, since it is a parametric test, it usually requires a normally distributed timeseries dataset to conduct the test. Among the various time-series trend detection tests, the M-K test is most commonly used to determine the trends of rainfall, temperature, river discharge, etc. However, the M-K test [60] is a non-parametric test, which means that it works for all types of distributions and it indicates whether the direction of inclination of the given parameter is positive or negative, but the magnitude of the inclination is often estimated using Sen's Slope estimate [61]. Mann-Kendall's statistics "S" can be calculated by Equations (1) and (2).
where x j and x i are sequential data values and n is the number of time-series data. The positive and negative values of S represent a "high tendency" and "low inclination", respectively. If the number of records used for the test is greater than 10 (N > 10), then the average distribution statistic is approximately equal to zero and the variance of S is calculated using Equation (3). Usually, more than ten records are used to analyze the trend of the time series.
where t i is the number of bonds specified in sample i. Then, Mann-Kendall's statistics Zc (standard average deviation) for N > 10 is calculated using Equation (4).
Zc follows the standard normal distribution. A positive Zc value shows an increasing trend, while the negative Zc data give a decreasing trend for the period. If Zc |> zα/2, the null hypothesis is rejected.

Sen's Slope Estimation
The magnitude of the direction generated by the M-K test can be calculated using Sen's slope estimator. Sen's slope estimator does not calculate the general straight and linear trend value of the time series data, where the true slope which changes in time is estimated using a procedure developed by Sen (1968). The test is widely used to estimate the magnitude of the trend in rainfall ranges over time. Slope pairs can be calculated for all data using Equation (5); where T i is the slope and x j and x k are the data values at time j and k, respectively. The mean of the n values of T i is encoded as Sen's slope estimator (Q i ) and is calculated using Equation (6); n is 0dd, , 1 2 T (n/2) + T (n+2)/2 , n is even

Results
As described in the sections above, a cross validation process was performed on CHIRPS data as an additional task. The CHIRPS data were analyzed according to the steps stated above and mapped the spatial distribution of rainfall trends at the district and climatic zone levels. This is beyond just presenting general statistical parameters.

Cross Validation of CHIRPS Data
Although nine rainfall locations were used to determine the correlation of rainfall, the monthly rainfall variations from 2015 to 2017 for the Hanwella, Pitabadra, Siyambalanduwa, and Tantirimale rainfall stations are shown in Figure 2. Hanwella and Pitabadra are located in the wet zone and Siyambalanduwa and Tantirimale can be observed in the north and south of the dry zone, respectively. In 2015, the SWM was average, and in 2016 as well as in 2017, it was clear that there was an increase in wetness. Nevertheless, the other two years in the dry zone apart from 2015 were marked as dry. The other important thing is that the second monsoon rains are better during these three years for wet zone locations compared to the dry zone. Since CHIPRS data are used entirely for this study, performing a cross validation process to observe the correlations with location-based data would give an initial positive impression on the study. For this purpose, scatter-plot charts were used to determine the correlation between ground-measured rainfall data and CHIRPS data, as shown in Figure  3, and the Pearson correlation coefficient was used to determine the magnitude of the correlation. Instead of performing this for the entire study period, only three years, 2015, 2016, and 2017, were used, which cover all climatic extremes, such as droughts, unusual rainfall, and severe floods, during that particular year range.
For rainfall locations which were selected for cross validation, the corresponding CHIRPS pixel's precipitation values were extracted to determine the accuracy between them. The maximum correlation on the coefficient is represented by Figure 3b and is approximately 0.91, Station-based rainfall measurements and CHIRPS data show an excellent correlation; therefore, cross validation performed on CHIRPS data provides a positive impression for the study. Since CHIPRS data are used entirely for this study, performing a cross validation process to observe the correlations with location-based data would give an initial positive impression on the study. For this purpose, scatter-plot charts were used to determine the correlation between ground-measured rainfall data and CHIRPS data, as shown in Figure 3, and the Pearson correlation coefficient was used to determine the magnitude of the correlation. Instead of performing this for the entire study period, only three years, 2015, 2016, and 2017, were used, which cover all climatic extremes, such as droughts, unusual rainfall, and severe floods, during that particular year range.
For rainfall locations which were selected for cross validation, the corresponding CHIRPS pixel's precipitation values were extracted to determine the accuracy between them. The maximum correlation on the coefficient is represented by Figure 3b and is approximately 0.91, Station-based rainfall measurements and CHIRPS data show an excellent correlation; therefore, cross validation performed on CHIRPS data provides a positive impression for the study. rainfall, and severe floods, during that particular year range.
For rainfall locations which were selected for cross validation, the corresponding CHIRPS pixel's precipitation values were extracted to determine the accuracy between them. The maximum correlation on the coefficient is represented by Figure 3b and is approximately 0.91, Station-based rainfall measurements and CHIRPS data show an excellent correlation; therefore, cross validation performed on CHIRPS data provides a positive impression for the study.   Figure 4 shows the changes in the annual rainfall of Sri Lanka, as well as wet, dry, intermediate, and semi-arid zones from 1989 to 2019. Furthermore, the spatial and temporal variability of rainfall can be identified, all of which usually show a variety of variance patterns, but overall, there is a tendency for rainfall to increase. Considering only the average annual rainfall of Sri Lanka, it can be seen that the maximum rainfall was received in 2010, 2011, and 2014. During these three years, the country experienced major floods [58].  Figure 4 shows the changes in the annual rainfall of Sri Lanka, as well as wet, dry, intermediate, and semi-arid zones from 1989 to 2019. Furthermore, the spatial and temporal variability of rainfall can be identified, all of which usually show a variety of variance patterns, but overall, there is a tendency for rainfall to increase. Considering only the average annual rainfall of Sri Lanka, it can be seen that the maximum rainfall was received in 2010, 2011, and 2014. During these three years, the country experienced major floods [58].  The other notable feature is the decrease in the total rainfall in Sri Lanka from 2015 to 2019. However, in 2016 and 2017, the wet zone received the highest recorded rainfall in its 30-year history. In contrast to the wet zone, the dry zone showed relatively low rainfall from 2015 to 2019, resulting in drought/dry conditions during the "Maha" season in 2016 to 2018. Therefore, the rainfall variation is very close to the terrestrial reality when viewed at the country or region level, as shown in Figure 4. On the other hand, this is one of the excellent examples to show the accuracy of CHIRPS data. The other notable feature is the decrease in the total rainfall in Sri Lanka from 2015 to 2019. However, in 2016 and 2017, the wet zone received the highest recorded rainfall in its 30-year history. In contrast to the wet zone, the dry zone showed relatively low rainfall from 2015 to 2019, resulting in drought/dry conditions during the "Maha" season in 2016 to 2018. Therefore, the rainfall variation is very close to the terrestrial reality when viewed at the country or region level, as shown in Figure 4. On the other hand, this is one of the excellent examples to show the accuracy of CHIRPS data. Table 1 shows the descriptive statistical parameters called mean, maximum, minimum, standard deviation, median, and coefficient of variance of the annual rainfall at the district level. The highest average annual rainfall between 1989 and 2019 was recorded in Kalutara district (3340.2 mm), and the lowest was recorded in Kilinochchi district (1109.7 mm). The observed specialty of the analysis of the maximum and minimum rainfall time series is that the maximum average rainfall value of 4489.4 mm was recorded in the year 2019 in the Kalutara district. Further study shows that the lowest maximum, 1619.8 mm in the Jaffna district, was observed in 1993, which is almost three times lower than the peak recorded in the Kalutara district in 2019. Considering the four climatic zones of Sri Lanka, the highest annual average rainfall in the wet zone is 2876.7 mm, followed by the intermediate zone (1979.3 mm), the dry zone (1520.6), and the semi-arid zone (1264.6 mm), and it is represented by Table 2. Another notable feature of this is that the average annual rainfall variability is higher in the wet zone than in the dry zone. The coefficient of variance (CV) is generally used to understand the inter-annual variability and variability level of rainfall. The CV value varies from 12.1 to 25.9 in all districts of Sri Lanka, with a maximum in Kilinochchi and a minimum in Nuwara Eliya. According to the [37] classification of CV, 6 out of 25 districts fall into the high variability class and the rest into the middle class. The peculiarity here is that all districts in the wet and intermediate zones belong to the middle CV class only, with about 50% of the districts in the dry and semi-arid zones showing high variability. The districts with the highest rainfall variability are Mannar, Trincomalee, Mullaitivu, Kilinochchi, Jaffna, and Batticaloa. Since rainfed agriculture is common in most of these districts, it can have a significant impact due to high rainfall variability.

Descriptive Statistics Generated Using CHIRPS Data
Sen's slope and Kendall's tau values were calculated at a 0.05 significance level using the M-K trend test for the annual rainfall (1987-2019) covering 25 districts in Sri Lanka. Figure 5 shows how the values of Sen's slope and Kendall's tau vary for each district. The two-tailed Mann-Kendall trend test used for this study is of 5% significance (Kendall's tau) and has 95% confidence in Sen's slope. The results of the M-K trend test on the annual rainfall show that there is a trend for an increase in rainfall in all districts and it is very clear that except for four districts, the other entire districts show a significant rainfall trend. The districts that do not show a significant trend are Ampara, Batticaloa, Jaffna, and Kilinochchi (Table A1).  (Figure 1), there is an increase in rainfall in both the wet and dry zones. Nevertheless, it is essential to study them separately to understand the precipitation trend in those climatic zones.

Discussion
This section mainly discusses the trend of seasonal rainfall over the past 30 years in Sri Lanka. Analysis of the seasonal rainfall trend helps farmers to manage their cultivation practices better and plan the crops. Furthermore, these analyses are of great help in managing long-term hazards such as floods and drought. The climate of Sri Lanka can be divided into two monsoons and two inter-monsoon seasons, namely, the southwestern  (Figure 1), there is an increase in rainfall in both the wet and dry zones. Nevertheless, it is essential to study them separately to understand the precipitation trend in those climatic zones.

Discussion
This section mainly discusses the trend of seasonal rainfall over the past 30 years in Sri Lanka. Analysis of the seasonal rainfall trend helps farmers to manage their cultivation practices better and plan the crops. Furthermore, these analyses are of great help in managing long-term hazards such as floods and drought. The climate of Sri Lanka can be divided into two monsoons and two inter-monsoon seasons, namely, the southwestern monsoon (SWM), the northeastern monsoon (NEM), and the first and second inter-monsoon (FIM, SIM), respectively.
The two-tailed Mann-Kendall trend test was also employed to identify seasonal rainfall trends with 5% significant and 95% confidence levels. The variation of the values determined by the M-K test for the four seasons of rainfall is graphically illustrated in Figure 6. For the Kalutara and Hambantota districts, the variation in annual, SWM, NEM, FIM, and SIM rainfall is illustrated in Figure 7, respectively, and the Sen's slope is shown together with the rainfall variation. This gives an obvious idea of the direction and magnitude of the rainfall trend. An important point that emerges from the analysis of M-K values in NEM is that no district shows a statistically significant trend in rainfall. The most important thing to note here is that about 40% of the districts represent a decreasing trend. The districts showing these negative values belong to all four climatic zones-wet, dry, intermediate, and semiarid-namely Ampara, Badulla, Batticaloa, Colombo, Galle, Hambantota, Kalutara, Matale, Matara, and Moneragala, respectively. The districts of Ampara, Moneragala, Batticaloa, and Badulla show a negative trend despite the fact that NEM is expected to receive more rainfall. Ultimately, this indicates that these districts are more prone to drought during the "Maha" season. Sri Lankan farmers had to experience this phenomenon in 2014, 2016, and 2017 [56].
When further interpreting the results of the M-K test, as shown in Figure 6 and Table  A3, there is an increased tendency for precipitation to increase for all districts during the FIM season. However, 17 districts (68%) out of the total districts show a statistically significant increasing trend (at 0.05 significant level). The highest trend of rainfall of 9.39 mm/year is recorded in the Kegalle district and the lowest, 1.27 mm/year, is in the Jaffna district. Batticaloa, Galle, Hambantota, Jaffna, Kilinochchi, Mannar, and Matara are the seven districts which do not show a statistically significant rainfall trend.
During the SIM, only Galle and Matara districts show a statistically significant rainfall trend of 6.28 and 6.09 mm/year, respectively. This shows that there has been no significant change in the intensity of rainfall in both the NEM and SIM seasons in the last 30 years. Table 2 shows the statistical parameters such as average, maximum, minimum, standard deviation, mean, and variance coefficients of long-term (1989-2019) annual rainfall based on different climatic zones. Average precipitation values in the wet, dry, and intermediate zones are still well within their classification ranges, but the average value in the semi-arid zone is higher than the maximum in its range. Furthermore, analysis of SWM NEM FIM SIM However, the results of the M-K trend test show that there is an increase in rainfall in all districts of Sri Lanka during the SWM, but only seven districts show a significant trend. All these districts fall in the country's wet zone, namely Colombo, Gampaha, Galle, Kalutara, Ratnapura, Matara, and Kegalle, respectively. Similar to the annual rainfall trend, the highest trend was recorded in SWM from Kalutara (21.943 mm/year) district and the lowest from Kilinochchi (0.978 mm/year) district (Table 3).
An important point that emerges from the analysis of M-K values in NEM is that no district shows a statistically significant trend in rainfall. The most important thing to note here is that about 40% of the districts represent a decreasing trend. The districts showing these negative values belong to all four climatic zones-wet, dry, intermediate, and semiarid-namely Ampara, Badulla, Batticaloa, Colombo, Galle, Hambantota, Kalutara, Matale, Matara, and Moneragala, respectively. The districts of Ampara, Moneragala, Batticaloa, and Badulla show a negative trend despite the fact that NEM is expected to receive more rainfall. Ultimately, this indicates that these districts are more prone to drought during the "Maha" season. Sri Lankan farmers had to experience this phenomenon in 2014, 2016, and 2017 [56].   (Table A2). From the above calculations, it can be concluded that there is a significant increase in rainfall from 1989 to 2019 in all climatic zones of Sri Lanka. Maps of annual rainfall trends of the district and climatic zones are shown in Figure 9. When further interpreting the results of the M-K test, as shown in Figure 6 and Table 3, there is an increased tendency for precipitation to increase for all districts during the FIM season. However, 17 districts (68%) out of the total districts show a statistically significant increasing trend (at 0.05 significant level). The highest trend of rainfall of 9.39 mm/year is recorded in the Kegalle district and the lowest, 1.27 mm/year, is in the Jaffna district. Batticaloa, Galle, Hambantota, Jaffna, Kilinochchi, Mannar, and Matara are the seven districts which do not show a statistically significant rainfall trend.
During the SIM, only Galle and Matara districts show a statistically significant rainfall trend of 6.28 and 6.09 mm/year, respectively. This shows that there has been no significant change in the intensity of rainfall in both the NEM and SIM seasons in the last 30 years. Table 2 shows the statistical parameters such as average, maximum, minimum, standard deviation, mean, and variance coefficients of long-term (1989-2019) annual rainfall based on different climatic zones. Average precipitation values in the wet, dry, and intermediate zones are still well within their classification ranges, but the average value in the semi-arid zone is higher than the maximum in its range. Furthermore, analysis of the coefficient of variance demonstrates that greater inter-annual variability cannot be detected in any climate zone. Figure 8 shows non-parametric M-K test results and the Sen's slope conducted with a 0.05 significance level for wet, dry, intermediate, and semi-arid climate zones of Sri Lanka. It is noteworthy that all of these climatic zones show a statistically significant trend. The highest trend (31.301 mm/year) is in the wet zone and the lowest (11.549 mm/year) is in the semi-arid zone. Increases in rainfall of 14.521 mm/year in the dry zone and 17.27 mm/year in the inter-zone are also shown (Table A2). From the above calculations, it can be concluded that there is a significant increase in rainfall from 1989 to 2019 in all climatic zones of Sri Lanka. Maps of annual rainfall trends of the district and climatic zones are shown in Figure 9.   Although the districts of Jaffna and Kilinochchi were included in the dry zone, the annual average precipitation of the two was found to be less than the rainfall margin (1250-1750 mm) in the dry zone. These two districts seem to be moving out of the dry zone and into the semi-arid zone. Further study of the intermediate zone and its distribution in the districts shows that although the Kurunegala District represents about 70% of the intermediate zone's land area, the average annual rainfall of the district is less as 1762 mm. A similar situation is observed in the Monaragala district as well (Table 1). That is, this value is best suited to represent the dry zone. However, in Badulla, Matale, and Kandy districts, the average rainfall is well within the intermediate-zone class margin. Considering all these factors, the main point to be concluded is that the margin separating the dry zone and the intermediate-zone has changed slightly.  Although the districts of Jaffna and Kilinochchi were included in the dry zone, the annual average precipitation of the two was found to be less than the rainfall margin (1250-1750 mm) in the dry zone. These two districts seem to be moving out of the dry zone and into the semi-arid zone. Further study of the intermediate zone and its distribution in the districts shows that although the Kurunegala District represents about 70% of the intermediate zone's land area, the average annual rainfall of the district is less as 1762 This study suggests that the methodology followed in this study could be adopted for developing countries in tropical and subtropical regions that suffer from a lack of local rainfall measurements to analyze rainfall variability over time. Moreover, rainfall station data do not show a good spatial variability of rainfall in terms of the distribution of rainfall stations. This study emphasizes that monthly CHIRPS data can be used as a suitable replacement for station rainfall for precipitation trend analysis to represent the effect of climate variability in Sri Lanka. The findings of this study have the potential to be used as an indication of climate change in Sri Lanka and to provide guidance for decision-makers on disaster risk management and mitigation processes.

Conclusions
For the past 31 years, from 1989 to 2019, each district's rainfall trends covering 25 districts of Sri Lanka have been studied using the M-K test and Sen's slope estimator. The results of the M-K non-parametric statistics test show an increase in annual rainfall for every district in the country. However, the annual rainfall in all districts except four in the dry zone increased significantly during the study period, and Jaffna, Batticaloa, Kilinochchi, and Ampara can be identified as the districts which did not show such an increase. Seasonal precipitation analysis shows that only the wet zone districts show a significant increase in SWM rainfall. However, 40% of the districts in the NEM show a negative trend and the other 60% show a positive trend, but no statistically significant trend in any district. Districts showing negative trends are spread over all climatological zones in Sri Lanka. Considering the trends in rainfall over the past 31 years in the climatic zones of Sri Lanka, it can be concluded that there is a significant tendency for precipitation to increase in all climatic zones (wet, dry, intermediate, and semi-arid) of Sri Lanka. Here, the maximum increase is recorded in the wet zone and the minimum increase is recorded in the semi-arid zone. A good example to prove this is the impact on the salt pans in the Hambantota district. That is to say, the salt production decreases as the district receives rainfall throughout the year. Further confirming these trends, it has been pointed out that the incidence of lightning in Sri Lanka has also increased from 1998 to 2014 [62]. Considering all these factors, the main point to be concluded is that the margin separating the dry zone and the intermediate-zone has changed slightly.
However, some studies [21,39,42] have shown that there is a decrease in annual rainfall trends in some districts. Nevertheless, this study seems to have generated a slightly opposite opinion. The main reason for this is that the study used raster data (CHIRPS data) generated by taking into account the spatial variability of rainfall, but previous studies used data provided by precipitation stations that did not consider spatial variability. Using one or two rainfall stations for one province or district, it is not very successful to apply those values to the whole district or province. Based on the results of this study, it can be concluded that there is a tendency for an increase in annual rainfall trends in all four climate zones.
With increasing rainfall at both the annual and seasonal level, it is likely that there will be an increased risk of floods in the southern and western provinces in the future. Finally, with the results of this study, the main point that should be brought to the attention of decision-makers is that the eastern and southeastern districts may face severe droughts in the future due to the declining rainfall trend in NEM. Therefore, it is advisable to introduce effective drought and flood management and preparedness measures with special attention to these areas. However, the findings from rainfall trend analysis through CHIRPS data can be used effectively to capture the effect of rainfall on climate change in Sri Lanka. Data Availability Statement: Daily grided rainfall data (CHIRPS data) are available at Climate Hazards Group InfraRed Precipitation corrected with stations data (https://data.chc.ucsb.edu/ products/CHIRPS-2.0/global_daily/netcdf/p05/ (accessed on 31 December 2020)). Location-specific ground-measured rainfall data are available at the Irrigation Department of Sri Lanka (https://www. irrigation.gov.lk (accessed on 31 December 2020)).

Acknowledgments:
The authors are thankful to the Climate Hazard Groups for providing the longterm observatory data (CHIRPS data), and the Irrigation Department of Sri Lanka for providing the relevant ground-measured rainfall data. Assistance and facilities provided by the International Water Management Institute and Department of Physics, University of Colombo are highly acknowledged.

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