Analyzing Rainfall Trends Using Statistical Methods across Vaippar Basin, Tamil Nadu, India: A Comprehensive Study

: The Vaippar basin in southern India is economically important for rainfed and irrigated agriculture, mainly depending on the northeast monsoon (NEM) during October–December, and any changes in rainfall patterns directly affect crop ecosystems. This study aimed to analyze spatio-temporal rainfall changes using the monthly data from 13 scattered rain gauge stations in the Vaippar basin, India. They were converted into gridded rainfall data by creating 26 equally spaced grids with a spacing of 0.125 ◦ × 0.125 ◦ for the period between 1971 and 2019 through interpolation technique. Three methods, namely Simple Linear Regression (SLR), Mann–Kendell/modified Mann– Kendell (MK/MMK), and Sen’s Innovation trend analysis (ITA), were employed to detect trends and magnitudes for annual and seasonal gridded rainfall series. The results showed significant trends at 2.3%, 7.7%, and 44.6% of grid points using SLR, MK/MMK, and ITA methods, respectively. Notably, ITA analysis revealed significant trends in annual and NEM rainfall at 57.69% and 76.92% of the grid points, respectively, at a 5% significance level. The southwestern and central parts of the basin exhibited a higher number of significant upward trends in annual rainfall. Similarly for the NEM season, the south-eastern, central, and extreme southern parts experienced significant upward trend. The western part of the basin exhibited significantly upward trend with a slope value of 2.03 mm/year, while the central part showed non-significant downward trend with a slope value of − 1.89 mm/year for the NEM series. This study used the advantage of ITA method, allowing for exploration of monotonic/non-monotonic trends, as well as subtrends of low, medium, and high rainfall segments within the series. The key findings of this study serve as a scientific report from a policy perspective, aiding in the preparation and management of extreme climate effects on land and water resources in the Vaipaar basin.


Introduction
Whether it is rainfed agriculture or irrigated agriculture, rainfall remains the primary source of water for crop production.In countries like India, which heavily rely on agriculture, achieving grain self-sufficiency has been a significant accomplishment.However, the production is resource intensive, focused mainly on cereals, and biased towards specific regions, all while facing increasing stress on water resources [1].The state Tamil Nadu, one of the major contributors to food production in India, includes 17 major river basins, with approximately 2.4 million hectares irrigated by surface water through major, medium, and minor schemes [2].According to a recent study, the rainfed cropland area in Tamil Nadu was estimated to be 2.57 million hectares, accounting for 19.66% of the total geographical area of the state.The Virudhunagar (0.14 million hectares) and Thoothukudi district (0.17 million hectares), located within the Vaippar basin, contributed to higher rainfed cropland areas [3], and have been considered in this research work.
In the context of rainfall variability, uncertain distribution of rainfall poses a serious obstacle to agriculture [4].The availability of rainwater for both rainfed and irrigated agriculture is becoming scarce due to uncertain rainfall patterns.Moreover, extreme climate-induced hazards such as droughts and floods are becoming more common due to hydro-climatological variability [5].The trends of these events have been linked to changes in rainfall patterns [6].Therefore, understanding the historical pattern of rainfall variation and referring to scientific reports on rainfall trends are crucial for planning water conservation efforts and formulating mitigation measures for extreme climate events [7].Additionally, these details are essential for minimizing underestimation or overestimation of design parameters for water infrastructure [8].
In the recent years, trends in hydro-meteorological variables have gained considerable attention in different parts of the world [4].Different statistical methods are presently available for detecting trends of hydro-meteorological variables and used by many researchers.Each method has its own merits and demerits.Generally, the trend detection methods are divided into parametric and non-parametric methods.Many scientists use non-parametric methods for trend analysis because these methods are not sensitive to outliers [9].They can be particularly useful for analyzing non-normally distributed series with missing values [10], and they do not rely on any assumptions about the nature of the data [11].The non-parametric methods such as the Mann-Kendall (MK) test, modified MK test, Spearman rank order correlation (SRC), Kendall rank correlation (KRC), and innovative trend analysis (ITA) are commonly followed by many researchers for trend detection at different time scales and significance levels [12][13][14][15].Equally, parametric methods such as simple linear regression (SLR) method are extensively used for detecting the monotonic long-term trend in the rainfall time series [16,17].The main advantage of SLR method is that it measures the statistical significance for testing hypothesis on the estimated slope and also provides the magnitudes of the parameters considered for analysis [18].Most of the researchers used non-parametric tests for trend assessment [19] and some of the studies explored the comparison of non-parametric and parametric tests for trends [20].The magnitude of rainfall trend (mm/year) can be determined using Sen's Slope estimator (SSE) and Simple Linear Regression (SLR) tests [15].
Recently, the innovative trend analysis (ITA) method developed by Sen (2012) [8] has gained more attention around the world and it was test verified by many researchers [14,17,44,47,52,63].The main advantage of ITA method is analyzing trends by providing graphical forms of presentation and without any limitations, such as non-normality, serial correlation, and size of data in the time series.In addition, ITA method provides a robust and powerful result with minimum error.It is also possible to detect the monotonic and non-monotonic trends in a way that time series are divided into different subcategories of the time series such as high, medium, and low zones [69,70].
This study aims to explore the rainfall variations, spatial patterns of trends, and their magnitudes of annual and seasonal rainfall series of the Vaippar basin at micro-level using classical rainfall statistical methods.The uniqueness in this research work is the use of spatially interpolated gauge rainfall data applied at micro-level for trend detection and identification of subtrends within the rainfall series over space and time.

Study Area
The Vaippar River Basin, located in the southern part of Tamil Nadu, India is a significant river basin in the region.Situated between latitudes 8 • 97 ′ N and 9 • 78 ′ N and longitudes 77 • 24 ′ E and 78 • 37 ′ E, it encompasses a total catchment area of 5320 km 2 .The basin is bounded on the north by the Vaigai and Gundar basins, on the south by the Tamaraparanibasin, on the west by the Western Ghats, and on the east by the Gulf of Mannar (Bay of Bengal).The basin area spans across four districts, with Virudhunagar accounting for 68%, Thoothukudi for 20%, Madurai for 7%, and Tirunelveli for 5% of the total area.The Vaippar River originates from the Echamalai mottai, Neduntheri mottai, and Kiladiparai hill ranges of the Western Ghats, situated near Sivagiri in Tirunelveli district.It begins at an elevation of 165 m above mean sea level and flows predominantly in an easterly and southeasterly direction for a distance of 146 km before joining the Gulf of Mannar.The catchment area of the Vaippar basin encompasses hilly regions such as Kodaliparai mottai, Vasudevanallur reserve forest, Periyasudangi malai, and others.These mountain ranges lie in the rain shadow regions of the Western Ghats, resulting in relatively low rainfall.The entire catchment area of the Vaippar basin is within the boundaries of Tamil Nadu state.The basin has been further divided into 13 sub-basins, namely: (1) Nichabanadhi, (2) Kalingalar, (3) Deviar, (4) Nagariyar, (5) Sevalperiyar, (6) Kayalkudiar, (7) VallampattiOdai/Uppodai, (8) SindapalliUppodai, (9) Arjunanadhi, (10) Kousiganadhi, (11) Uppathurar, (12) Senkottaiyar, and (13) Vaippar.The location of the Vaippar basin is depicted in Figure 1.
The basin experienced frequent drought with a range from four to eight years per drought [71].Agricultural land covers approximately 74% of the total geographical area, while forested areas account for 10% of the area.Wasteland occupies 8% of the total geographical area, while settlements and water bodies together cover less than 8% of the basin's total geographical area.Out of the total agriculture area in the basin, cultivable land represents 43%.This cultivable land is primarily utilized for cultivating water-intensive crops such as paddy, sugar cane, and banana [72].Additionally, cotton, non-paddy, and dry crops are also grown in the basin.The basin's total irrigated area accounts for 24% of the cultivable land.The remaining 76% of the cultivable land relies mainly on rainwater for irrigation [71].The basin experienced frequent drought with a range from four to eight years per drought [71].Agricultural land covers approximately 74% of the total geographical area, while forested areas account for 10% of the area.Wasteland occupies 8% of the total geographical area, while settlements and water bodies together cover less than 8% of the basin's total geographical area.Out of the total agriculture area in the basin, cultivable land represents 43%.This cultivable land is primarily utilized for cultivating water-intensive crops such as paddy, sugar cane, and banana [72].Additionally, cotton, non-paddy, and dry crops are also grown in the basin.The basin's total irrigated area accounts for 24% of

Data Used
Daily rainfall data from 13 rain gauge stations distributed across the Vaippar basin were collected for this study.The data were obtained from the State Ground and Surface Water Resources Data Center, Public Works Department, and Water Resources Organization in Chennai.The location details of each rain gauge station and the period of data utilized in this study are presented in Table 1.Standard quality control measures were undertaken to identify outliers and errors in rainfall data for all 13 rain gauges.Any potential outliers, including missing data and instrument errors, were checked and corrected.To facilitate the analysis, the daily rainfall data for each station was processed to derive monthly data.Subsequently, the monthly rainfall series was further organized into seasonal rainfall series.The seasons considered for this analysis were as follows: Southwest Monsoon (SWM) from June to September, Northeast Monsoon (NEM) from October to December, winter season in January and February, summer season from March to May, and an annual rainfall series.The space-time correlation analysis of the mean monthly rainfall data with the rain gauge stations was conducted and is presented in Table A1.A strong correlation exists within the mean monthly rainfall series among the rain gauge stations.

Methodology
In this study, the monthly rainfall data collected from the 13 rain gauge stations were spatially interpolated to generate gridded rainfall data.A grid with a spacing of 0.125 • in both latitude and longitude was employed to create 26 individual grids that spanned the entire basin area.The rainfall data available for each location in each year during the period specified in Table 1 underwent spatial interpolation.Using spatially interpolated monthly rainfall data, gridded rainfall data was prepared.Using the gridded rainfall data, trends in the monthly, annual, and seasonal rainfall series were estimated.Three different methods were employed for trend analysis: Mann-Kendall (MK) test, simple linear regression (SLR), and Sen's innovative trend analysis (ITA) method.The trends were evaluated by estimating the slopes using Sen's slope, SLR, and ITA methods to quantify the magnitude and direction of the trends in the rainfall data, and the percentage change of magnitude of trends from mean rainfall.The overall methodology used in this study is depicted as a flowchart in Figure 2.

Spatial Interpolation of Rainfall
Spatial interpolation techniques play a crucial role in estimating values at locations where no observed data is available, using known data values.Several methods, such as inverse distance weighting (IDW), splines, and kriging, are commonly employed for spatial analysis of various variables [73][74][75].In the present study, the IDW approach was specifically utilized for spatially interpolating rainfall data across the Vaippar basin.This method considers the proximity of known data points to the location of interest and assigns weights accordingly, resulting in an interpolated surface that provides valuable in-

Raingauge data
Interpolated data

Changes in mean rainfall
Figure 2. The overall process involved for rainfall characterization.

Spatial Interpolation of Rainfall
Spatial interpolation techniques play a crucial role in estimating values at locations where no observed data is available, using known data values.Several methods, such as inverse distance weighting (IDW), splines, and kriging, are commonly employed for spatial analysis of various variables [73][74][75].In the present study, the IDW approach was specifically utilized for spatially interpolating rainfall data across the Vaippar basin.This method considers the proximity of known data points to the location of interest and assigns weights accordingly, resulting in an interpolated surface that provides valuable insights into the spatial distribution of rainfall within the basin.
The IDW interpolation technique assigns weights to each control point based on the inverse of their distances from the interpolated point.This means that closer control points have a greater influence on the interpolated value compared to those farther away.The technique assumes that each control point has a local influence that diminishes as distance increases.In IDW interpolation, the output value at a given location is determined by a specified number of nearest points or all points within a specified radius.The power parameter of the IDW method controls the significance of the surrounding points in determining the interpolated value.A higher power value reduces the influence of distant points, resulting in a stronger emphasis on the values of nearby points in the interpolation process [76].
The general form of IDW approach [77,78] is given in Equation (1): where zn i is the new value for grid i; z j is the value of m nearest neighbours; d j is the distance to m nearest neighbours; p is the exponent of distance.For this study, the exponent of distance was taken as two for spatial interpolation of the rainfall.This study employed the Inverse Distance Weighting (IDW) technique instead of the Thiessen polygon method to estimate the gridded rainfall across the basin.The Thiessen polygon method often results in a crude approximation of rainfall spatial variation due to attribute variations associated with each rain gauge station, ranging from 14 to 30% [79,80].To overcome this limitation, the Vaippar basin was divided into smaller square grids.Given the constraints of rainfall data availability and the uneven distribution of rain gauge stations within the basin, spatial interpolation of data at smaller grids was necessary to address these challenges and obtain a more accurate representation of rainfall patterns.
The total area of the Vaippar basin was divided into 26 grids, each with dimensions of 0.125 • (205 km 2 ).These grids accounted for approximately 3.85% of the total area, equivalent to 5320 km 2 .Monthly rainfall data recorded at 13 stations were spatially interpolated using the IDW method in QGIS 3.30.2resulting in gridded rainfall data for the period from 1971 to 2019.The gridded rainfall data was then utilized to analyze the spatial patterns of significant trends in rainfall over the study period.

Rainfall Variability
In this study, the variability of rainfall has been assessed using the coefficient of variation (CV).The CV is a statistical measure that represents the ratio of the standard deviation to the mean of a rainfall data series.It is often used to quantify the relative variability or dispersion of rainfall.A higher CV indicates greater variability, while a lower CV suggests more stability or consistency in the data [27,81].The CV for the rainfall events is expressed as a percentage (%).As mentioned by Hare (2003) [82], the CV can be used to classify the degree of variability of rainfall events into three categories: (i) when the CV is less than 20, rainfall events are considered to have less variability or consistent pattern of rainfall; (ii) when the CV is between 20 and 30, rainfall events are classified as having moderate variability or fluctuation; (iii) when the CV exceeds 30, rainfall events are categorized as having high variability indicating a less stable or more unpredictable rainfall.

Simple Linear Regression
Simple linear regression (SLR) analysis is a widely used parametric model for identifying monotonic long-term trends in monthly, seasonal, and annual rainfall time series [83].SLR analysis establishes a relationship between two variables, often employed to determine the slope of hydro-meteorological variables over time.A positive slope indicates an upward trend, while a negative slope suggests a downward trend.The primary advantage of SLR analysis is its ability to provide a measure of significance through hypothesis testing on the slope, along with quantifying the rate of change.To obtain the percentage change during the specified period, the slope values are multiplied by the duration of the study period in years [16].
In this study, the simple linear regression models for the annual and seasonal time series were developed between the gridded rainfall data and time for trend detection [40,84,85].The test statistic "t", which follows a Student's t-distribution with (n − 2) degrees of freedom (where "n" represents the length of the data points in the time series), was calculated to examine the significance of trend.The null hypothesis (Ho) of a zero slope is rejected when the calculated test statistic "t" value surpasses the critical value "t α/2 " at a given significance level.

Mann-Kendall (MK) Trend Analysis
The non-parametric Mann-Kendall (MK) test is a commonly used method for exploring trends in hydro-meteorological data.The trends identified by this test in the time series are monotonic, meaning they can be increasing or decreasing without assuming linearity.The test was originally proposed by Mann (1945) [86] and the test-statistic distribution was subsequently derived by Kendall (1975) [87].The null hypothesis for this test assumes that the data is independent and randomly ordered.One advantage of this test is that it does not require any assumption of normality.However, it only indicates the direction of significant trends, without providing information about their magnitude [86,87].The MK test statistic S indicates the number of positive differences minus the number of negative differences for all the considered differences.Using S value, its mean, and variance, the standardized test statistic Z was calculated for annul and seasonal rainfall series.A positive or negative Z value indicates an upward/downward trend [88].The Z statistic follows a normal distribution.To test for statistical significance, compare the calculated Z value with critical values obtained from the standard normal distribution tables at two significance levels (α = 5% and α = 10%).The null hypothesis (Ho) of no trend is rejected if the calculated value of Z is greater than the table value of Z 1−α/2 [13,30,89].

Modified Mann-Kendall Test
In the Modified Mann-Kendall (MMK) test proposed by Hamed and Rao (1998) [90], the autocorrelation coefficient (ACC) of the rainfall series is computed and then subjected to testing at various levels of significance.If the ACC value is determined to be significant, the MMK test is subsequently [13,23,89].
Initially, the modified variance was computed by integrating the Auto-Covariance Correction (ACC) at lag-i.Subsequently, this modified variance was utilized in the original Mann-Kendall (MK) test, leading to the calculation of the modified Mann-Kendall Z c statistic.To assess the significance of the trend, the modified Mann-Kendall Z c statistic was subjected to testing against threshold levels.For instance, significance levels of 10% and 5% were represented by threshold values of 1.645 and 1.96, respectively.

Autocorrelation Analysis of Time Series
The autocorrelation (serial correlation coefficient, r) analysis is a useful tool for evaluating the presence of randomness and periodicity within a time series at different lag periods [91].To calculate autocorrelation coefficients (ACC), the normalized anomaly of the rainfall series, obtained from gridded rainfall data, along with the long-term average and standard deviation of annual and seasonal rainfall, can be utilized.Positive or negative If the ACC is close to zero for different time lags, it suggests that the data points in the time series are randomly distributed, indicating no dependence.On the other hand, the presence of serial correlation, as examined through autocorrelation analysis, is an important preliminary test before performing the Mann-Kendall (MK) trend analysis.This is because the MK analysis requires the input data to be serially independent.Otherwise, the presence of positive serial correlation in the data can lead to an overestimation of trend significance.The lag-1 ACC specifically detects the serial correlation in the data series.The statistical significance of serial correlation can be tested by employing a normally distributed statistic at a significance level of "α", taking into account the lag period and the length of the series [94].
The ACC (r k ) of the rainfall time series at lag-k can be computed using the rainfall data of the time series and the total length of the time series, where "k" represents the time lag.The value of r k ranges between +1 and −1.A value of zero for r k signifies that the series is random for all lag-k values [15].To evaluate the statistical significance of the ACC in the time series data, testing should be conducted at both upper and lower confidence limits.If the ACC turns out to be statistically insignificant, the Mann-Kendall (MK) test can then be applied to the original time series [95,96].The null hypothesis (Ho) suggests the absence of serial correlation within the time series, while the alternative hypothesis (Ha) suggests the presence of some serial correlation rather than pure randomness.In this study, the hypothesis was examined using the lag-1 autocorrelation coefficient.The upper and lower confidence limits can be calculated [45,97] at a significance level denoted as α (e.g., for α = 10%, Z = ±1.645,and for α = 5%, Z = ±1.96)[45].When the r k value falls within the range of the critical values {(r k ) upper < r k < (r k ) lower}, the null hypothesis Ho: r k = 0 is rejected, while Ha: r k ̸ = 0 is accepted.This indicates that the time series is not random and exhibits some level of persistence or serial correlation [11].

Sen's Slope
The magnitude of trends in monthly, seasonal, and annual time series was determined using Sen's slope, a non-parametric method introduced by Sen (1968) [98].Sen's slope requires equally spaced data in a time series [11].One of the main advantages of Sen's slope over simple linear regression (SLR) is its resilience to data errors, outliers, and extreme observations [15].This method is particularly useful when assuming a linear trend line.
The slope estimate Q for each pair of data values is calculated using where x j and x k are the data values at time j and k, respectively (j > k).If there are n data values x j in the time series, there will be a total of N = n(n − 1)/2 slope estimates Q i .The N values of Q i are then ranked from the smallest to the largest.The Sen's slope is determined based on whether N is odd or even . A positive value of Q indicates an upward trend, while a negative value represents a downward trend.A Q value of zero indicates no trend in the time series [15,30].

Innovative Trend Analysis
The Innovative Trend Analysis (ITA) method has been successfully utilized for trend detection in hydro-meteorological variables [8,52].This method is simple, allowing for easy identification and visualization of trends in high, medium, and low data set on the trend line [70].Unlike non-parametric trend identification tests, the ITA method does not require restrictive assumptions such as data series independence, normality, or data length [8,47,99].It involves plotting all data points of the time series in a Cartesian coordinate system and comparing them with a diagonal straight 1:1 line [44,100].The construction procedure for ITA is provided below: (a) Divide the monthly, annual, or seasonal rainfall time series into two equal halves.
Arrange each half series in ascending order.
(b) Place the first half of the time series on the X-axis and the second half on the Y-axis.Plot the rainfall series as a scatter diagram.(c) Draw a straight 45 • line diagonally in the scatter graph representing a 1:1 relationship.
Divide the plot into upper and lower half triangles.(d) If the scattered points align perfectly on the 45 • line, the series does not exhibit a trend.(e) When the scattered points lie within the upper half triangle, the rainfall series is said to have an increasing trend.(f) When the scattered points lie within the lower half triangle, the rainfall series is said to have a decreasing trend.(g) A time series is considered to have a monotonic trend if all scattered points lie above or below the upper or lower half triangles.Non-monotonic trends occur when some scattered points are located in the upper half triangle and others in the lower half triangles [8].(h) The trend of low, medium, and high values can be observed in the scatter graph.(i) To test the significance of the trend, a null hypothesis (Ho) is considered: there is no significant trend if the calculated slope value (S) is below the critical value (S cri ).The alternative hypothesis (Ha) states that there is a significant trend when S > S cri .(j) The slope (S) of the trend is computed using the formula [101]: where, y 1 and y 2 represent the mean of the first and second halves of the rainfall time series, and n is the number of data points in the rainfall time series.(k) The standard deviation (σ s ) of the trend slope is calculated using the formula [101]: where σ represents the standard deviation and ρ y 1 y 2 is the cross-correlation coefficient between the means of the two-half series.(l) By utilizing the confidence limits (S cri ) for a standard normal probability density function, the confidence limits of the trend slope are calculated at a significance level of α using the formula [101]: If the slope value (S) falls beyond the lower or upper confidence limits, the null hypothesis of no significant trend is rejected at the α significance level.

Percentage Change in Magnitude of Trend from Mean Rainfall
The percentage change in magnitude of trend from mean rainfall of seasonal and annual time series can be calculated using the following expression [31,102,103]: where PC represents the percentage change of rainfall (%), n is the length of time series in years, S is the magnitude of the trend slope calculated using methods such as Simple Linear Regression (SLR) coefficient, Sen's slope, or Innovative Trend Analysis (ITA) slope, and x is the mean value of the time series.This equation allows for quantifying the percentage of change in rainfall over the specified time period, relative to the mean value of the time series.

Rainfall Variability
The mean rainfall and coefficients of variation (CV) for the monthly, seasonal, and annual rainfall of 26 grid points in Vaippar during the period 1971-2019 is presented in Tables A2 and A3.During the study period, the Vaippar basin experienced an annual average rainfall of 762.57mm.The highest rainfall of 891.93 mm was observed in the G08 grid point, while the lowest rainfall of 570.2 mm occurred in the G18 grid point.Out of 26 grid points, 12 grid points located mostly on the western side of the basin experienced annual rainfall that exceeded the average annual rainfall.The NEM is the major rainy season in the study area, contributing approximately 54.7% of the annual rainfall.Within the NEM season, the majority of the rainfall, around 85%, occurred in two months: October and November.The month of December contributed to approximately 15% of the summer rainfall.Among the grid points, the highest NEM rainfall of 503.54 mm was recorded in the G01 grid point, while the lowest NEM rainfall of 361.39 mm was recorded in the G18 grid point.Apart from the NEM season, the remaining seasons, namely SWM, winter, and summer, contributed 19.5%, 5.4%, and 20.4% of the annual rainfall, respectively.The maximum monthly rainfall was observed in October, accounting for 24.1% of the annual rainfall.The month of November followed closely with 22.3% contribution, while September contributed 9.9%, and December contributed 8.3% of the annual rainfall.
In the Nagariyar sub basin (G01 grid point), the maximum rainfall was recorded during the months of February, March, April, November, December, and the winter season.The G08 grid point (part of Arjunanadhi) experienced high rainfall in January, September, and October, annually, and during the summer season.On the other hand, the G20 grid point (part of Kousiganadhi) recorded higher rainfall in June, July, August, and during the Southwest Monsoon (SWM) season.The grid points located in the lower part of the Sinkottaiyar sub basin (G25) and Sindapalli Uppodai sub basin (G14 and G18) registered the minimum rainfall.The G25 grid point recorded the minimum monthly rainfall from January to September, November, SWM, winter, and summer seasons.The G18 grid point experienced the minimum rainfall during October, while the G14 grid point recorded the minimum rainfall in December.
The coefficient of variation (CV) of monthly, seasonal, and annual rainfall calculated for each grid point is shown in Table A3.Rainfall variability as classified by Hare (2003) [82], shows the CV for annual mean rainfall ranges from 22.9% to 31.5%, indicating moderate to high variability in rainfall distribution across all grids.The G09 grid point recorded the highest CV, while the G06 grid point recorded the lowest CV for annual rainfall.Most of the grid points exhibited moderate variability in annual rainfall, except for G08, G09, and G14, which are located at the centre of the basin.Rainfall variability was found to be greater in seasonal rainfall compared to annual rainfall.However, it is worth noting that the NEM season experienced relatively lower variability compared to the other seasons.Among the months, October (51.47%) and November (61.33%) exhibited lower coefficients of variation (CV) compared to the other months.

Trends of Annual and Seasonal Rainfall Series
The temporal trends were identified using the SLR, MK/MMK, and ITA methods at different grid points of Vaippar basin for annual and seasonal rainfall series.The calculated values of the SLR t test, MK/MMK test statistic (Z test), and ITA (slope values) were spatially mapped for each grid point by IDW interpolation method using QGIS 3.30.2software.The magnitude of the trends was identified and percentage changes in mean rainfall were also calculated.The comparison of number of grid points expressing the significant trends and correlation among the different trend methods were also attempted.The results are discussed in the subsequent sections.

Trends of Annual and Seasonal Rainfall Series by Simple Linear Regression
The spatial distributions of trends in annual and seasonal rainfall, detected by the SLR method at the 5% and 10% significance levels, are shown in Figure 3.The analysis of temporal trends using the simple linear regression for annual and seasonal rainfall series showed that approximately 73% of grid points exhibited non-significant upward trends in the annual rainfall series.Although the NEM rainfall is a major contributor to the annual rainfall, the trend pattern was not similar, with 50% of grid points showing a non-significant downward trend.The winter and summer rainfall series demonstrated that 57.7% and 76.9% of grid points, respectively, displayed non-significant upward trends.
trends and correlation among the different trend methods were also attempted.The results are discussed in the subsequent sections.

Trends of Annual and Seasonal Rainfall Series by Simple Linear Regression
The spatial distributions of trends in annual and seasonal rainfall, detected by the SLR method at the 5% and 10% significance levels, are shown in Figure 3.The analysis of temporal trends using the simple linear regression for annual and seasonal rainfall series showed that approximately 73% of grid points exhibited non-significant upward trends in the annual rainfall series.Although the NEM rainfall is a major contributor to the annual rainfall, the trend pattern was not similar, with 50% of grid points showing a nonsignificant downward trend.The winter and summer rainfall series demonstrated that 57.7% and 76.9% of grid points, respectively, displayed non-significant upward trends.Among the five-rainfall series (annual and seasonal), only two series (SWM and summer) exhibited 7.69% of significant trends at the extreme ends of the basin.In the SWM series, an upward trend was observed at the G03 grid point and a downward trend at the G23 grid point, both at a 10% significance level.The summer rainfall series showed an upward trend at the G01, G02, G14, and G26 grid points at a 10% significance level, and at the G04, G07, and G25 grid points at a 5% significance level.A significant downward trend was noticed at the G23 grid point for the summer series.Consequently, the expected decrease in rainfall will not have a significant impact on water availability, as these sea- Among the five-rainfall series (annual and seasonal), only two series (SWM and summer) exhibited 7.69% of significant trends at the extreme ends of the basin.In the SWM series, an upward trend was observed at the G03 grid point and a downward trend at the G23 grid point, both at a 10% significance level.The summer rainfall series showed an upward trend at the G01, G02, G14, and G26 grid points at a 10% significance level, and at the G04, G07, and G25 grid points at a 5% significance level.A significant downward trend was noticed at the G23 grid point for the summer series.Consequently, the expected decrease in rainfall will not have a significant impact on water availability, as these seasons contribute very little to the annual rainfall.

Trends of Annual and Seasonal Rainfall Series by MK/MMK Test
The Mann-Kendall (MK) method was applied to the annual and seasonal rainfall series at different grid points of the Vaippar basin to identify significant trends at the 5% and 10% significance levels using the Z-test.The MK test was conducted for the annual and seasonal series, taking into account the auto-correlated non-significant series at lag-1.The modified MK (Z c ) test was performed only for statistically significant auto-correlated series.
Autocorrelation analysis was conducted to select the appropriate trend analysis method and evaluate the performance of both the original and normalized rainfall series.The autocorrelation coefficient for 26 grid points at lag-1 period for annual and seasonal rainfall series were worked out and the correlogram is presented in Figure 4.The upper and lower bound were decided by the 95% confidence interval to test the limits of the autocorrelation coefficient.The autocorrelation was considered as significant if it is greater than or lower than ±0.28.Since autocorrelation was found to be significant for five rainfall series viz ANL-G10, ANL-G12, SWM-G25, WIN-G23 and WIN-G25, the modified MK test was performed for these five series.and seasonal series, taking into account the auto-correlated non-significant series at lag-1.
The modified MK (Zc) test was performed only for statistically significant auto-correlated series.
Autocorrelation analysis was conducted to select the appropriate trend analysis method and evaluate the performance of both the original and normalized rainfall series.The autocorrelation coefficient for 26 grid points at lag-1 period for annual and seasonal rainfall series were worked out and the correlogram is presented in Figure 4.The upper and lower bound were decided by the 95% confidence interval to test the limits of the autocorrelation coefficient.The autocorrelation was considered as significant if it is greater than or lower than ±0.28.Since autocorrelation was found to be significant for five rainfall series viz ANL-G10, ANL-G12, SWM-G25, WIN-G23 and WIN-G25, the modified MK test was performed for these five series.
The spatial pattern of trends in annual and seasonal rainfall, identified using the MK methods, is presented in Figure 5. Compared to the SLR test and ITA method, very few statistically significant trends were observed in the annual and seasonal rainfall series for both the MK and MMK tests.The temporal patterns of trend detected by MK test indicated that approximately 73% of the grid points for the annual series and 58% of the grid points for the NEM series showed non-significant upward trends.Furthermore, 54% of the grid points for the SWM series displayed a non-significant downward trend.
The MK test detected a downward trend at the 5% significance level in the G23 grid point for the summer series.A total of five rainfall series were tested using the modified MK method, which detected a significant downward trend at the G23 grid point for the SWM series, while the same grid point exhibited a significant upward trend for the winter season.The spatial pattern of trends in annual and seasonal rainfall, identified using the MK methods, is presented in Figure 5. Compared to the SLR test and ITA method, very few statistically significant trends were observed in the annual and seasonal rainfall series for both the MK and MMK tests.The temporal patterns of trend detected by MK test indicated that approximately 73% of the grid points for the annual series and 58% of the grid points for the NEM series showed non-significant upward trends.Furthermore, 54% of the grid points for the SWM series displayed a non-significant downward trend.

G1...G26 Grid Points
The MK test detected a downward trend at the 5% significance level in the G23 grid point for the summer series.A total of five rainfall series were tested using the modified MK method, which detected a significant downward trend at the G23 grid point for the SWM series, while the same grid point exhibited a significant upward trend for the winter season.

Trends of Annual and Seasonal Rainfall Series by ITA Method
The grid-wise trend parameters for the annual and seasonal rainfall series of the Vaippar basin, as detected by the ITA method, are presented in Tables 2-4.These tables provide valuable insights into the trends observed in the rainfall patterns across different grids within the basin.As observed from the tables, the slope values of annual and seasonal rainfall series during the period 1971-2019 fall outside the lower and upper confidence limits (CL) for the particular grid point, suggesting existence of a significant trend in the rainfall pattern.
The slope values of the annual rainfall series are presented in Table 2.Among the 26 grid points analyzed, it is noteworthy that 15 grid points (57.69%) exhibited significant trends at both the 5% and 10% significance levels.Out of these significant trends, 11 grid points (42.3%) displayed a significant upward trend, while four grid points showed a significant downward trend.The ITA analysis of the seasonal rainfall series displayed in Tables 2-4 revealed that significant trends were observed in the NEM and summer rainfall series in 20 and 21 grids i.e., 76.92 and 80.77%, respectively.In contrast, the SWM and winter rainfall series detected significant trends in nine and eight grids i.e., 34.62 and 30.77%, respectively.For the NEM season, 53.85% of the grids showed a significant upward trend, while 23% of the grid points showed a significant downward trend.In the case of summer rainfall, a significant upward trend was observed in 73.1% of the grid points.SWM exhibited a significant upward trend in 26.9% of the grid points.Winter, on the other hand, experienced a significant downward trend in 19.2% of the grid points.These findings suggest that the NEM and summer seasons experienced more widespread and pronounced changes in rainfall patterns compared to the SWM and winter seasons.The spatial distribution of slope values obtained through the ITA method, along with their significance for seasonal and annual rainfall series, is illustrated in Figure 6.Grid points located in the southwestern and central parts (G01, G02, G03, G09, G10) of the basin exhibited a higher number of significant trends in annual rainfall.In terms of the NEM season, significant trends were observed across almost all parts of the basin.The southeastern, central, and extreme southern parts experienced a significant upward trend, while the eastern parts displayed a significant downward trend.Significant upward trends were identified in the southwestern parts for the SWM season.For the summer season, significant upward trends were observed in the western, southern, and eastern parts of the basin.Notably, the western parts exhibited a significant upward trend during the winter season.Among the grid points, G10 stood out with a notably higher number of significant downward trends in monthly, annual, and seasonal rainfall series.This indicates a consistent and significant decrease in rainfall at that specific location.From the spatial analysis of the Vaippar basin, it is evident that the southeastern, central, and extreme southern parts have exhibited a positive increase in annual rainfall.Moreover, when considering the different seasons, it is notable that the NEM season has displayed widespread positive trends across various parts of the basin.It is important to highlight that these positive trends in rainfall can have significant implications for water availability, agricultural productivity, and overall ecological balance within the Vaippar basin.

Identification and Nature of Subtrends by ITA Method
One of the most important features of the ITA method is its ability to identify the subtrend of a rainfall series [8].Additionally, the scattered plot allows for the detection of both monotonic and non-monotonic upward or downward trends in a series.To conduct a comprehensive analysis for identifying subtrends within a rainfall series, scatter diagrams can be used, with data points plotted on 1:1 line graphs.These scatter diagrams are then divided into three segments based on rainfall depth: low, medium, and high rainfall segments.In the scenario where the annual series exhibits a combination of different trend patterns within the series, it is referred to as a non-monotonic trend.Conversely, if the series demonstrates a consistent upward or downward pattern, it is considered a monotonic trend.Additionally, it is also possible for the series to exhibit no trend if the scatter points align closely with the 1:1 line or if the slope values are zero.
The grid-wise subtrends and nature of trends observed in the annual and seasonal rainfall series for the Vaippar basin are displayed in Table 5.In the table, the subtrends of rainfall series, categorized as low, medium, and high rainfall segments, are indicated by upward or downward arrows.The nature of the trend is represented by whether it is monotonic or non-monotonic (upward/downward).Analyzing the annual and NEM series, it can be observed that five grid points (G01, G02, G03, G07 and G25) located at the southwestern parts of the basin exhibited monotonically upward trends.On the other hand, for the SWM, 10 grid points showed monotonically downward trends.In the case From the spatial analysis of the Vaippar basin, it is evident that the southeastern, central, and extreme southern parts have exhibited a positive increase in annual rainfall.Moreover, when considering the different seasons, it is notable that the NEM season has displayed widespread positive trends across various parts of the basin.It is important to highlight that these positive trends in rainfall can have significant implications for water availability, agricultural productivity, and overall ecological balance within the Vaippar basin.

Identification and Nature of Subtrends by ITA Method
One of the most important features of the ITA method is its ability to identify the subtrend of a rainfall series [8].Additionally, the scattered plot allows for the detection of both monotonic and non-monotonic upward or downward trends in a series.To conduct a comprehensive analysis for identifying subtrends within a rainfall series, scatter diagrams can be used, with data points plotted on 1:1 line graphs.These scatter diagrams are then divided into three segments based on rainfall depth: low, medium, and high rainfall segments.In the scenario where the annual series exhibits a combination of different trend patterns within the series, it is referred to as a non-monotonic trend.Conversely, if the series demonstrates a consistent upward or downward pattern, it is considered a monotonic trend.Additionally, it is also possible for the series to exhibit no trend if the scatter points align closely with the 1:1 line or if the slope values are zero.
The grid-wise subtrends and nature of trends observed in the annual and seasonal rainfall series for the Vaippar basin are displayed in Table 5.In the table, the subtrends of rainfall series, categorized as low, medium, and high rainfall segments, are indicated by upward or downward arrows.The nature of the trend is represented by whether it is monotonic or non-monotonic (upward/downward).Analyzing the annual and NEM series, it can be observed that five grid points (G01, G02, G03, G07 and G25) located at the southwestern parts of the basin exhibited monotonically upward trends.On the other hand, for the SWM, 10 grid points showed monotonically downward trends.In the case of the summer rainfall series, a monotonically upward trend was identified at 13 grid points.
For the annual rainfall series, upward trends were observed in the low rainfall segment for 21 grid points, in the medium rainfall segment for 14 grid points, and in the high rainfall segment for five grid points.As for the NEM rainfall, 24 grid points showed upward trends in the low rainfall segment, 15 grid points in the medium rainfall segment, and five grid points in the high rainfall segment.The high rainfall segments, as classified for annual and NEM rainfall, were found in the G01, G02, G03, G07, and G25 grid points, which are located in the southwestern parts of the basin.On the other hand, the SWM season exhibited the minimum number of grid points registering upward trends in the low, medium, and high rainfall segments.A higher number of grid points experienced upward trends in the low and medium rainfall segments for the winter series.Conversely, the summer rainfall series exhibited upward trends in the low, medium, and high rainfall segments for a larger number of grid points.

Magnitude of Trend in Rainfall Series
The magnitude of the trend (β), computed using the SLR method, and the Sen's slope (Q) are presented in Table 6.Considerable variability was observed in both the physical value and sign of the magnitude for different grid points.The SLR slopes indicate that the highest positive magnitude of 2.25 mm/year was recorded at the G02 grid point, while the lowest negative magnitude of −2.62 mm/year was observed at the G08 grid point for annual rainfall series.Furthermore, it is worth noting that these magnitudes indicate the absence of a statistically significant trend in the annual rainfall series.The upward trend was found to be more pronounced in the northern and southern parts of the basin.However, in the case of NEM rainfall, which contributes significantly to the annual rainfall, the trend magnitude did not follow a similar pattern as the annual rainfall.Grid point G25 (1.45 mm/year) and G08 (−1.52 mm/year) registered higher positive and negative magnitudes for the NEM series, respectively.Additionally, for the NEM series, the northern parts of the basin exhibited a negative magnitude, while the southern parts showed a positive magnitude.The magnitudes of Sen's slope (Q) calculated for annual and seasonal rainfall as presented in Table 6 showed higher positive magnitudes of the non-significant trends.For the annual rainfall series, the highest magnitude of 3.89 mm/year was observed at the G14 grid point, while the lowest magnitude of −1.6 mm/year was recorded at the G23 grid point for a non-significant trend.Positive magnitudes were noticed for more than 69.2% of the grids in the annual rainfall series.Higher positive magnitudes were observed in the southwestern and central parts of the basin for the annual series.In the case of NEM, the highest non-significant magnitude of 1.74 mm/year was observed at the G14 grid point, while the lowest magnitude of −1.66 mm/year was recorded at the G08 grid point.Positive magnitudes were noticed for more than 57.7% of the grids in the NEM rainfall series.Negative magnitudes were noticed in the western and central parts of the basin for the NEM series, while the remaining part of the basin recorded positive magnitudes.
The slope values estimated using the ITA method for the annual and seasonal series are presented in Tables 2-4.The slope values of the annual series indicate that 42.3% of the grids exhibited positive magnitudes.Grid point G02 displayed a significantly upward trend with a slope of 4.47 mm/year, while G09 showed a non-significant downward trend with a slope value of −4.24 mm/year for the annual series.Regarding the seasonal rainfall, the NEM and summer series exhibited predominantly positive magnitudes, while the SWM and winter series showed negative magnitudes.Grid point G02 exhibited a significantly upward trend with a slope of 2.03 mm/year, whereas G09 showed a non-significant downward trend with a slope value of −1.89 mm/year for the NEM series.A similar pattern of positive and negative magnitudes, similar to the annual rainfall, was observed for the NEM series.Positive magnitudes were observed in the western, central, and eastern parts of the basin for both the annual and NEM series.Among the three methods considered, the ITA method consistently estimated a higher magnitude of trend for all the rainfall series compared to the other methods.

Percentage Change in Magnitude of Trend
The percentage change in magnitude of trend from the mean values over the study period calculated for the seasonal and annual rainfall series using the SLR slope (β), Sen's slope (Q), and ITA slope (S) are presented in Table 7.Among three methods, the ITA slope exhibited a higher percentage change in the magnitude of the trend compared to SLR and Sen's slope.Additionally, it was observed that the percentage changes in magnitude from mean rainfall were lower in the annual and NEM rainfall series compared to the SWM, winter, and summer series.Positive values of percentage change in rainfall were more prominent in the annual, NEM, and summer series across all three methods.Specifically, in the annual series, the extreme western, central, and eastern parts of the basin showed positive magnitudes, while the northern side exhibited negative magnitudes for all three methods.Regarding the NEM series, the western parts displayed negative magnitudes, and the eastern parts showed positive magnitudes for SLR and Sen's slope methods.However, for the ITA slope method, this pattern was reversed, with positive magnitudes following a similar pattern to the annual series.

Comparison of Trend Methods
The number of grid points that exhibited a significant upward or downward trend in seasonal and annual rainfall series is displayed in Figure 7.A total of 130 rainfall series, including annual and four seasonal series, were analyzed for trends in 26 grid points using the SLR, MK/MMK, and ITA methods.The test results were compared based on the number of significant upward or downward trends identified.Among the three methods used, namely SLR, MK/MMK, and ITA, significant trends were observed in 2.3%, 7.7%, and 44.6% of the grid points, respectively.Notably, the SLR and MK/MMK methods detected significant trends only in the SWM, winter, and summer series, while no trends were identified in the annual and NEM series.It is noteworthy that SLR and MK/MMK methods only identified significant trends in the SWM, winter, and summer series, while they did not detect trends in the annual and NEM series.It is worth mentioning that all ten significant trends detected by the SLR test were also identified by the ITA method.However, the significant trend detected by the MK test did not align with the ITA method.One significant trend identified by the MMK method was also detected by the ITA method.
Table 8 presents the correlation between the test statistics of the MK and SLR methods with the slope values of the ITA method, regardless of the level of significance.It was observed that the MK test showed good correlation with the ITA method in the annual and summer series.Similarly, the correlation between the SLR test and the ITA method was found to be good for the annual, SWM, winter, and summer series.The study findings suggest that the ITA method outperforms traditional trend detection methods.The ITA method proved beneficial in detecting many significant trends that could not be identified by the traditional methods in the annual and seasonal rainfall series.It effectively revealed hidden trends in rainfall series across the grid points.Many researchers have validated the Sen's (2012) [8] ITA methodology for different hydro-meteorological variables in various parts of the world [4,14,17,50,89,[104][105][106][107][108][109][110].Unlike classical methods, the ITA method does not require prewhitening prior to its application [8].Trends of low, medium, and high data can be easily observed using this method [99,111].Among the three methods used, namely SLR, MK/MMK, and ITA, significant trends were observed in 2.3%, 7.7%, and 44.6% of the grid points, respectively.Notably, the SLR and MK/MMK methods detected significant trends only in the SWM, winter, and summer series, while no trends were identified in the annual and NEM series.It is noteworthy that SLR and MK/MMK methods only identified significant trends in the SWM, winter, and summer series, while they did not detect trends in the annual and NEM series.It is worth mentioning that all ten significant trends detected by the SLR test were also identified by the ITA method.However, the significant trend detected by the MK test did not align with the ITA method.One significant trend identified by the MMK method was also detected by the ITA method.
Table 8 presents the correlation between the test statistics of the MK and SLR methods with the slope values of the ITA method, regardless of the level of significance.It was observed that the MK test showed good correlation with the ITA method in the annual and summer series.Similarly, the correlation between the SLR test and the ITA method was found to be good for the annual, SWM, winter, and summer series.The study findings suggest that the ITA method outperforms traditional trend detection methods.The ITA method proved beneficial in detecting many significant trends that could not be identified by the traditional methods in the annual and seasonal rain-fall series.It effectively revealed hidden trends in rainfall series across the grid points.Many researchers have validated the Sen's (2012) [8] ITA methodology for different hydrometeorological variables in various parts of the world [4,14,17,50,89,[104][105][106][107][108][109][110].Unlike classical methods, the ITA method does not require prewhitening prior to its application [8].Trends of low, medium, and high data can be easily observed using this method [99,111].One disadvantage of this test is that it must be applied to each recorded series individually [112].The possibility of presenting results graphically in this method enables the easy observation of hidden subtrends and helps in identifying trends in extreme values [63,113].
The findings of the present study indicate that a majority of the grid points in the western and eastern parts of the Vaippar basin show a significant increasing trend in rainfall.This increasing trend is significant as it can contribute to an increase in runoff, which can be utilized for water management purposes, particularly for tapping and harnessing the runoff water in addition to the existing water conservation structures.
Conversely, it was also observed that some grid points exhibit a significant downward trend in rainfall.This decline in rainfall presents various challenges and implications for surface and groundwater management.One of the key challenges is the increased reliance on groundwater extraction for crop irrigation, leading to the depletion of groundwater reserves, drought-related issues, and diminished soil moisture [44].
Given these observations, it is crucial to carefully consider the changes in rainfall patterns and trends in long-term catchment-scale water management strategies.Adapting and planning for these changes can help mitigate the potential risks associated with declining precipitation, such as implementing measures to enhance water conservation, exploring alternative water sources, and promoting sustainable agricultural practices that optimize water usage.Long-term water management strategies should take into account the evolving rainfall patterns to ensure efficient and sustainable utilization of water resources in the Vaippar basin.

Conclusions
This study aimed to analyze the variation and trends of seasonal and annual rainfall series in the Vaippar basin using gridded rainfall data from 1971 to 2019.The SLR, MK/MMK, and ITA methods were employed to examine the trends, magnitudes, subtrends, and nature of the trend.In order to account for spatial variability, monthly rainfall data from 13 rain gauge stations within the basin were spatially interpolated using the inverse distance weighing (IDW) method under GIS environment.To further refine the spatial representation, the basin was divided into 26 grids, each covering approximately 205 km 2 , and gridded rainfall data was generated from the interpolated gauge data.The key findings and conclusions of this study are summarized below.
The basin experienced moderate variability in annual rainfall, lower variability during the NEM season, and higher variability in other seasons.Among the three methods (SLR, MK/MMK, and ITA), significant trends were detected in 2.3%, 7.7%, and 44.6% of the grid points, respectively.The SLR and MK/MMK methods detected significant trends only in the SWM, winter, and summer series.The significant trend detected by the MK test did not align with the ITA method, but the significant trends detected by the SLR test were consistent with the ITA method.The ITA method indicated that 57.69%, 76.92%, and 80.77% of the grid points exhibited significant trends at 5% and 10% significance levels in annual, NEM, and summer rainfall, respectively.However, the SWM and winter series showed less than 35% significant trends.For the NEM season, 53.85% of the grids displayed a significant upward trend, which is a positive sign for improving water management.Grid points located in the southwestern and central parts of the basin showed a higher number of significant trends in annual rainfall.In terms of the NEM season, the southeastern, central, and extreme southern parts experienced a significant upward trend.
In the annual and NEM series, grid points located in the southwestern parts of the basin exhibited monotonically upward trends.Approximately 19.3% of the grid points in the southwestern parts of the basin showed upward trends in high rainfall segments, as classified for annual and NEM rainfall.Grid points in the western part of the basin exhibited a significantly upward trend with a slope of 2.03 mm/year, while the central part showed a non-significant downward trend with a slope value of −1.89 mm/year for the NEM series.Among the three methods considered, the ITA method consistently estimated a higher magnitude of trend for all the rainfall series compared to the other methods.Compared with traditional methods, the ITA method, which represents rainfall series graphically without making any assumptions, detected trends that were not identified by traditional methods.It facilitated the identification of monotonic or non-monotonic upward/downward trends and trends in low, medium, and high rainfall segments.
From the spatial analysis of the Vaippar basin, it is evident that the southeastern, central, and extreme southern parts have exhibited a positive increase in annual rainfall.Moreover, when considering different seasons, it is notable that the NEM season has displayed widespread positive trends across various parts of the basin.These positive rainfall trends have significant implications for water availability, agricultural productivity, and overall ecological balance within the Vaippar basin.The significant findings of this study will serve as a crucial scientific reference for policymakers, assisting in the preparation and management of extreme climate effects on land and water resources within and around the Vaipaar basin.

Figure 1 .
Figure 1.Geographical location of study area along with (a) rain gauge stations and (b) grid points in subbasin.

Figure 1 .
Figure 1.Geographical location of study area along with (a) rain gauge stations and (b) grid points in subbasin.

Figure 7 .
Figure 7. Grid points exhibited significant trend in Monthly, Seasonal and Annual rainfall series.

Figure 7 .
Figure 7. Grid points exhibited significant trend in Monthly, Seasonal and Annual rainfall series.

Table 1 .
Details of rain gauge station geographical location and data used.

Table 2 .
Grid-wise trend parameters detected by the ITA method for the annual and SWM rainfall series of the Vaippar basin.

Table 3 .
Grid-wise trend parameters detected by the ITA method for the NEM and winter rainfall series of the Vaippar basin.

Table 4 .
Grid-wise trend parameters detected by the ITA method for the summer rainfall series of the Vaippar basin.
S-ITA slope; σ-standard deviation; ρ-correlation coefficient between the means of two half series; σs-standard deviation of trend slope; LCL/UCL-lower/upper confidence limits; Ho: no significant trend; Ha: there is a significant trend; **-significance level at 5 and 10%.

Table 5 .
Subtrends of annual and seasonal rainfall series for Vaippar basin.

Table 7 .
Percentage change in magnitude of trend from mean for SLR, Sen' Slope, and ITA slopes for Vaippar basin.

Table 8 .
Correlation between the test statistics of MK and SLR with the slope values.

Table 8 .
Correlation between the test statistics of MK and SLR with the slope values.