Spatial distribution and estimation of rainfall trends and erosivity in the Epitácio Pessoa reservoir catchment, Paraíba, Brazil

The main goals of this study are to better understand the spatial and temporal variabilities in rainfall and to identify rainfall trends and erosivity for the period from 1963 to 1991 in the Epitácio Pessoa reservoir catchment, which is located in Paraíba, northeastern Brazil. This study analyzes annual rainfall trends on a regional scale by using monthly data from 13 rainfall stations. For this purpose, the nonparametric Mann–Kendall and Sen methods were used in the analysis. Descriptive statistics methods and interpolation techniques were also used for spatial–temporal analysis of the annual rainfall. A detailed statistical analysis applied to the time series of all the stations indicates that the rainfall presents substantial annual spatial–temporal variability and a negative trend (decrease) in the mean rainfall at most of the rainfall stations in the catchment during the study period. The results only showed a positive trend for the Soledade and Pocinhos stations. The distribution of positive and negative trends in the Epitácio Pessoa reservoir catchment is extremely irregular, and the changes in the study area are more significant compared to those identified in other studies.


Introduction
Global warming, and therefore changes in annual rainfall, has attracted the attention of researchers in different regions of the world (Mullick et al. 2019; Thomson et al. 2019;Morbidelli et al. 2018). One of the most significant consequences of global warming would be an increase or decrease in the magnitude and frequency of annual rainfall (Croitoru et al. 2013). Climate change is a very important topic and one of the main challenges at present. In addition, both the large regions of the world and the semiarid regions need more attention.
In the past decades, climate variations have been analyzed with special attention and the results of these studies are of interest in several areas of knowledge, especially for semiarid regions, where optimized management of water resources should be pursued and land degradation should be avoided, for example. The potential applications of this topic are numerous, including in the management of water resources and in the analysis of changes in ecosystems or in economic activities (Li et al. 2018;Moghim 2018;Guo et al. 2019). Changes in annual rainfall have been identified in several studies conducted at local, regional and global scales (Silva et al. 2015;Gocic and Trajkovic 2013;Hamlaoui-Moulai et al. 2013;Rutebuka et al. 2020).
In this study, rainfall variability is directly related to regional water management, which is of great socioeconomic and environmental importance, mainly in semiarid regions like Paraíba State. More than 85% of the area of Paraíba State is considered semiarid, and this state has historically suffered from constant periods of drought due to the position of the intertropical convergence zone, the El Niño-Southern oscillation and sea surface temperature, which are interrelated mechanisms that cause variability in rainfall and weather (Alves et al. 2017). The behavior of rainfall in time and space in this basin directly influences the water volume in the reservoirs within the region; furthermore, the semiarid northeastern Brazil is one of the world's most densely populated dryland regions (Pilz et al. 2018). In addition, severe drought events might deplete the drinking water supply, as has occurred in recent continuous dry years (2012)(2013)(2014)(2015)(2016)(2017).
The history of the Brazilian semiarid region, which could be considered one of the regions in the world most vulnerable to climate change (Marengo et al. 2018), is marked by droughts, water scarcity and consequently their effects, which manifest themselves in the most varied forms, such as environmental degradation, increased rural unemployment and hunger or poverty in or migration from the affected areas. Because of the irregular rainfall distribution in the region, which is below 600 mm per year in much of the region, the semiarid zone faces a chronic problem of water scarcity, which is an obstacle to the development of economic activities (Marengo 2008). Cycles of severe droughts and floods usually affect the region at intervals ranging from a few years to decades. Statistically, there are 18-20 years of drought every 100 years, which makes it difficult to analyze the variability of hydrological processes and consequently to obtain a more effective prediction of extreme events in the region (Galvão et al. 2005). In addition, estimations of future climate based on different atmospheric circulation models show that there will be global warming and an increase in climatic aridity in many zones globally. In the case of northeastern Brazil, projections predict an increase in droughts and decrease in water storage (Marengo and Bernasconi 2015;Fernandez et al. 2019).
More recently, much of Paraíba State, located in northeastern Brazil, has been experiencing one of the greatest droughts in its history (2012)(2013)(2014)(2015)(2016)(2017), which directly affects the local population (Gomes et al. 2017;Santos et al. 2019a, b;Dantas et al. 2020;Alvalá et al. 2017;Marengo et al. 2018). The risks arising from climate variability, whether natural or of anthropogenic origin, have raised great concern in the scientific and political milieus, in the media and in the general population, especially in regions that are more susceptible to natural risks, such as the Brazilian semiarid region ). There are some studies regarding the rainfall behavior trends, analysis of climate extremes and drought indices in the semiarid zone of Brazil (da Barbosa and Kumar 2016;Alves et al. 2017;Tomasella et al. 2018;Santos et al. 2019b;Marçal et al. 2019;Costa et al. 2020), mainly observing droughts in the last decades of the 20th century. However, none of the work provides a detailed study regarding rainfall trends, rainfall erosivity, land degradation and water management in the Epitácio Pessoa reservoir. The Epitácio Pessoa reservoir is a major reservoir in the semiarid region of Paraíba State that supplies 26 cities, including Campina Grande, one of the main cities in northeastern Brazil (De Medeiros et al. 2019). The annual cycle of rain and streamflow in the country varies between basins, and in fact, the interannual climate variability associated with El Niño or La Niña, or the variability in the tropical and southern Atlantic sea surface temperatures, may generate climatic anomalies that generate large droughts, such as those in the years 1877, 1983and 1998in northeastern Brazil (Andreoli and Kayano 2007. The droughts generated crises related to the public water supply in cities due to the constant reduction in the water stored behind dams after the beginning of the droughts (Dantas et al. 2020). As this region does not have groundwater recharge due to the crystalline basement and because it is a headwater basin, water storage depends on precipitation, which makes this basin extremely dependent on the region's drought events and the good management of water resources of the basin. In this sense, the analysis of the trend in the rainfall time series, coupled with geographic information systems, has been an important tool to identify and measure trends and to identify the existence of any deviations from trends. Thus, this study aims to identify possible trends in the rainfall time series and the erosivity in the Epitácio Pessoa reservoir catchment.

Characterization of the study area
The Epitácio Pessoa reservoir catchment drainage basin is located in the central portion of Paraíba State between 6° 50′ 00′′ S-8° 20′ 00′′ S and 36° 00′ 00′′ W-37° 25′ 00′′ W, constituting an area of 12,406 km 2 and a perimeter of 436 km (Fig. 1). This basin is located in the semiarid region of northeastern Brazil, more precisely on the Borborema Plateau. The Epitácio Pessoa reservoir catchment has smooth undulating relief and elevations ranging between 400 and 600 m above sea level; the mean monthly temperature is always above 26 °C . Its main rivers are the Upper Paraíba and the Taperoá rivers, which have intermittent flow regimes and originate in Serra do Teixeira County and Monteiro County, respectively, and terminate in the Epitácio Pessoa reservoir.
The predominant vegetation in the basin is caatinga (dry tropical forest) that presents a mosaic of different physiognomies with mostly small-sized vegetation together with cochineal cactus, cotton and agave crops (the latter two are of more commercial nature) and temporary crops such as maize and beans, which represent a small area of the basin (Souza et al. 2004;Santos et al. 2017). The Caatinga is a well-recognized ecological region that lies in the semiarid region hinterland of northeastern Brazil and that covers approximately 912,529 km 2 .
The climate in the Epitácio Pessoa reservoir catchment, according to De Medeiros et al. (2019) is equatorial tropical type (2d), hot semiarid, with a substantial temporal, spatial and interannual variability in rainfall; the dry season can last 8 months (June-January), and the wet season can last 4 months (February-May), as shown in Fig. 2. The historical rainfall in the basin indicates that the region presents an annual mean that ranges between 350 and 600 mm/year, with a higher concentration in a period of 2-4 months (January-April), corresponding to 65% of the total annual rainfall (Srinivasan and Paiva 2009). Regarding the climatic influence in the region, the uniform action of the Atlantic masses provides a uniform temperature and rainfall gradient in the region. The higher rainfall rates in this area occur due to the intertropical convergence zone, mesoscale convective complexes, upper tropospheric cyclonic vortexes and instability lines (Dantas et al. 2020).

Rainfall data and drought severity/anomalies
For this study, monthly rainfall data for 13 rainfall gauges were obtained from the Paraíba State Water Management Executive Agency (Agência Executiva de Gestão de Águas do Estado da Paraíba-AESA) for the period from 1963 to 1991 (Table 1). The choice of these rainfall gauges was based on the amount of missing data during the studied period. The total percentage of missing rainfall data was 2.8%, ranging from 0 to 4% among the stations. Any missing data were filled into form a complete time series for the available period. To fill the missing data, the regional weighting method was used (Eq. 1). If X is the station with missing data and a, b and c are the surrounding stations, the weighted average record at the three surrounding rainfall gauges were used to determine the rainfall of X gauge as: where P X is the variable that will represent the corrected data, M X is the arithmetic mean of the station with missing data, M a , M b and M c are the arithmetic means of neighboring stations, and P a , P b and P c are the data that were used to fill the missing values. These data are from the stations neighboring the station with missing data using linear distance (Table 2) from the same year.
The regional weighting method is generally used to fill in monthly to annual series, where the missing data of a station are filled using a weighting based on the data from

Fig. 2
Monthly mean rainfall hyetograph in the Epitácio Pessoa reservoir catchment  at least three neighboring stations, which should be from climatological regions similar to that of the station of interest and have a data series of at least 10 years ). The El Niño index was used to analyze the rainfall anomalies. Data from the Golden Gate Weather Services database (GGWS, 2019) were used to apply this index. The El Niño index indicates the intensity of the incidence of El Niño/La Niña phenomena from the mean sea surface temperature over three consecutive months in the El Niño 3.4 region in the Pacific Ocean. The range for each intensity class is as follows: (a) weak: 0.5-0.9, (b) moderate: 1.0-1.4, (c) strong: 1.5-1.9 and (d) very strong: ≥ 2.0. The El Niño index exhibits variation along its historical time series, which indicates high or low intensity of the reduction in rainfall in northeastern Brazil.

Analysis of rainfall trends
Statistical trend analysis methods assume the independence of the data to test the correlation between the data, which ensures the reliability of the results . Descriptive statistics techniques show the central position of the data (median) and whether there is symmetry or asymmetry in addition to showing the outliers, i.e., the values that diverge greatly from the other values in the series or are inconsistent with the series, which aid in understanding the data distribution. Among the various methods of trend analysis, two nonparametric methods were chosen: Mann-Kendall (Mann 1945;Kendall 1975), suggested by the World Meteorological Organization, and Sen's test (Sen 1968), which defines the magnitude of the trend. This integration of these two tests is used to detect significant trends and estimate the true slope of the trend. These methods have been widely used in hydrometeorological time series (Tabari et al. 2011;Gocic and Trajkovic 2013;Silva et al. 2015;Kisi et al. 2018). The Mann-Kendall method is applied in cases where the data values x i of a time series can obey the following model: where ƒ(t i ) is an increase or decrease in the monotonic function, and ε i are the residuals assumed to have the same distribution with zero mean. It is then assumed that the variance in the distribution is constant. In this method, the null hypothesis (H 0 ) is tested, i.e., the observations x i are ordered randomly against the hypothesis H 1 where there is a trend (Salmi et al. 2018). For the computation of this statistical test, both S statistics described in Santos et al. (2018) and the normal approximation (Z statistics) are exploited. Next, the variance S is calculated by the following equation: where n is the number of data values, p represents the number of groups with repeated data, and t i is the number of data values of the ith group. The values of S are used to calculate the test Z S : The presence of a trend is then evaluated by the value of Z s . A positive value indicates an upward trend, and a negative value corresponds to a downward trend. If the null hypothesis is rejected, the level of significance (α) is calculated. Significance levels of 0.05 (*), 0.10 (**) and 0.15 (+) were used in this work.
The other nonparametric method used was Sen's method. To estimate the slope of a trend using this method, the following equation was used: (2) where Q is the slope, and B is a constant. To estimate B, the n values of the difference x i − Q ti are calculated. The confidence intervals are then calculated at two different levels to assure the reliability of the test. Sen's method also makes it possible to identify whether or not there was a change in the trend and its magnitude, which is a technique widely used for this type of analysis Kisi et al. 2018).
To obtain the slope estimate given in the above equation, it is necessary to calculate the slope of all the data pairs: where x j and x k are the values of x in periods j and k, respectively. If there are n values x j in the time series, then there are N = n(n − 1)/2 estimates of slope Q i . Thus, the slope estimated by the Sen's method is the median of these N values of Q i . The N values of Q i are classified from lowest to highest and the Sen's estimate is expressed as follows: where Q med reflects data trend, while its value indicates the steepness of the trend. To determine whether the median slope is significantly different from zero, the confidence interval of Q med should be determined at a specific probability. The confidence interval about the time slope can be computed as follows: where Z 1−α/2 is obtained from the standard normal distribution table. In this paper, the confidence interval was computed at three significance levels (0.05, 0.10 and 0.15). M 1 = (n − C α )/2 and M 2 = (n − C α )/2 were then computed. The lower and upper limits of the confidence interval, Q min and Q max , are the M 1 th largest and the (M 2 + 1)th largest of the n-ordered slope estimates, respectively (Gilbert 1987). The slope Q med is significantly different from zero if the two limits (Q min and Q max ) have similar signs.
The Pettitt test uses the Mann-Whitney U test to find a single change point in the median of a series and returns a p value for the change point (Pettitt 1979). This test is robust to outliers and skewed distributions (Ryberg et al. 2020). In rare instances, the Pettitt test returns two change points because of ties in the data; in these instances, we used the first change point. The ranks of the y 1:n (r 1:n ) are used to calculate the Pettitt statistic, which is based on the Mann-Whitney two-sample test: The probable change point U is located at the k value at which the absolute value of P(k) is maximized.
The approximate probability for a two-sided test is calculated according to:

Spatial analysis of rainfall
Spatial analysis is the process of manipulating spatial information to extract new information and meanings from the original data (Golla et al. 2019). The geostatistical interpolation method was used to determine the spatial distribution of rainfall in the catchment. The interpolation of data is of paramount importance to assess the spatial distribution of the data. Thus, rainfall and erosivity interpolation maps were generated, as were interpolation maps for the trend test results for each rainfall station for the studied period. To obtain the distribution of the rainfall data using spatial interpolation, the inverse distance-weighted (IDW) interpolator was used, which is expressed by: where N is the number of data, v j is the value of point j, d j is the distance between the values of point j, and w(d) is the weighting function. To obtain the value of w(d), the following equation was used: where d min is the minimum distance, and d max is the maximum distance. The use of d min prevents the possibility of values with d = 0, and the use of d max avoids the use of very distant points.
Interpolation using the IDW interpolator is based on the weighted distance of sample points, i.e., the surface interpolated by this method is created based on the use of a weighting coefficient that controls the spatial estimation of values for the sample points. The main advantage of the method is simplicity, leading to variable results for a wide variety of data, with no issues for results that exceed the range of significant values (Caruso and Quarta 1998). This method is based on spatial dependence, i.e., it assumes that the closer an individual sample point is to the other, the greater the correlation of this individual sample point with its neighbors. Thus, greater weight is assigned to the closer individual sample points than to those furthest from the point to be interpolated.
The model consists of multiplying the values observed by the inverse of their respective distances to the reference point for interpolation of the values (Eq. 12). This statistical model considers the existence of an effect of distance and of any other factor, which is represented by the letter "p", i.e., the distance is raised to the power "p" so that different interpolated values can be obtained for the same distance. Rainfall spatialization using IDW was performed using ArcGIS 10.1 ® software.

Rainfall erosivity calculation
To calculate the rainfall erosivity, monthly rainfall data between January 1963 and December 1991 from 13 rainfall stations were used. The monthly rainfall erosivity of each station was determined using Eq. 14 (Da Silva 2004): where R m is the rainfall erosivity factor (MJ mm ha −1 h −1 y −1 ), P m is the average monthly precipitation depth (mm), and P is the average annual precipitation (mm). Figure 3 presents the annual variability in rainfall at the studied stations. High variability is observed between 1963 and 1991 in the Epitácio Pessoa reservoir catchment, with an annual mean of 495.98 mm and a standard deviation of 237.50 mm. The results show modest temporal trends but strong and coherent spatial trends in the rainfall data from the Epitácio Pessoa reservoir catchment. Essentially, the findings agree with studies of other semiarid regions (Santos et al. 2019b;Ramos and Durán 2014), which highlighted intraand interannual rainfall variability. The annual rainfall variations of approximately 350 mm are considered low, which caused an increase in the number of degraded areas during the study period, as reported by Silva et al. (2018). This expansion of the degradation area has been accelerated due to the severe droughts that recurred in this region in past decades, as shown by Tomasella et al. (2018), Awange et al. (2016) and Martins et al. (2018). Figure 4 represents the rainfall variability at the analyzed stations. The results show that several years exhibited large water deficits or above-average rainfall throughout the region. The Prata station showed the highest annual rainfall values when compared to the other stations. It can be observed that the number of outliers is low (not exceeding two per series) and that the Prata station stands out as having the greatest data dispersion.  Table 3 presents the mean values for each station and the variance (standard and mean deviations) of the data from the mean. The Congo station stands out due to its low standard deviations, which indicate that the values are closer to the mean, whereas the Camalaú and Prata stations presented the highest standard deviations. However, the Prata station indicated the highest standard deviation, unlike the other stations, which indicated relatively distant values. It should be noted that results for northeastern Brazil, especially the semiarid region, can be partly explained by variation in the general atmospheric circulation. Figure 5 shows the El Niño anomalies and the drought occurrence for northeastern Brazil between 1963 and 1991. Indeed, several studies have linked the decrease in rainfall to the El Niño-Southern Oscillation (ENSO), which influences the climate in South America and other parts of the world. In Brazil, this phenomenon affects parts of the northern, northeastern and extratropical areas of southern Brazil (Souza et al. 2018;da Silva et al. 2011). The climatic variability associated with the ENSO phases, especially rainfall anomalies, is quantifiable, and this information can be effectively used by crop managers to  decrease the associated risk or to make better use of forthcoming favorable climatic conditions (Cunha et al. 2001). According to de Souza et al. (2005), the rainfall in northeastern Brazil is strongly modulated by the combined effects of sea surface temperature anomalies in the Pacific and Atlantic Oceans. Therefore, the connection between the Pacific and Atlantic Oceans has received particular interest, and some researchers have suggested that the ENSO-related atmospheric teleconnection plays an important role in Atlantic climatic variability (Ropelewski and Halpert 1987;Saravanan and Chang 2000). Figure 6 shows the spatial distribution of rainfall in the Epitácio Pessoa reservoir catchment between 1963 and 1991. The highest and lowest rainfall values were identified in the southeastern and northeastern portions of the catchment, respectively. An analysis of the locations of the stations indicates that most of the closer stations tend to fall within the average rainfall range in the periods studied, except for the Cabaceiras (307.66 mm) and Prata (694.79 mm) stations, which are located far from the area of influence of the mean, and the Congo (446.07 mm) station, because it is located between the two stations.

Spatial-temporal analysis of annual rainfall trends
The Mann-Kendall and Sen's methods were applied to analyze the rainfall trends. Table 4 shows the results of trends and significance for the analyzed data. The values that reject the null hypothesis are highlighted in this table, i.e., values that indicate a trend. The Soledade, Coxixola, Caraúbas and Desterro stations showed positive trends (upward). Using the results for the Sen slope, indications of decreasing trends were observed in the annual rainfall series for the Pocinhos, Gurjão, Serra Branca, São José dos Cordeiros, Congo, Cabaceiras, Camalaú and Prata rainfall stations. The Camalaú, Barra de São Miguel and Prata Negative trends (downward) with higher significance levels were recorded at the Camalaú, Barra de São Miguel and Prata stations, with − 1.970, − 1.707 and − 1.444 for the Z S test and − 11.300, − 7.208 and − 9.808 for the Q test, respectively. At the other stations (Soledade, Serra Branca, São José dos Cordeiros, Pocinhos, Gurjão, Coxixola, Congo, Caraúbas, Cabaceiras and Desterro), the null hypothesis was accepted and there was no evidence of a trend. It can be observed that the negative slope of trends calculated using the Q test ranged from − 0.383 (São José dos Cordeiros station) to − 11.300 (Camalaú station). Climate change plays a key role in the rainfall variability in the environments of the region, especially including the Caatinga biome, where it compromises water and energy security and the subsistence agriculture in the region. According to Marengo et al. (2018), the Caatinga biome presents multiple stressors on natural and human systems that are derived in part from significant changes in climate and are influenced by land use and cover change. This variability also influences the extreme events  Figure 7a-m present the graphical representations of the Sen and Pettitt tests (magnitude and slope), which indicate the annual rainfall variability for each rainfall station. This absence of a significant trend in rainfall during the analyzed period is probably due to the substantial variability of rainfall that can be observed graphically (Fig. 7). The figures indicate that there is no change in the trend of the series studied, which are monotonic. In addition, they show the substantial variability in rainfall during the study period, especially for the Soledade (a), Coxixola (f), Caraúbas (h) and Desterro (m) stations, which also show evidence of trends in the Z S test. The results of the Pettitt test are presented in Table 5, which indicates that the null hypothesis H 0 (homogeneous data) was accepted for all time series. These results are close to those obtained by Sena andLucena (2013), Xavier et al. (2016) and Alves et al. (2017), who also analyzed rainfall trends in part of the studied basin.
This part of the paper is intended to describe regional variations in the spatial and temporal variability in the rainfall, i.e., to identify regions that exhibit similar rainfall variations. Figure 8a shows the spatial distribution of the results of the Z S test for the Epitácio Pessoa reservoir catchment. The lowest values occur in the Prata, Barra de São Miguel and Camalaú stations, which are located in the central and southeastern regions of the basin. The intermediate values cover a large portion of the basin, including the Congo, Serra Branca, São José dos Cordeiros, Cabaceiras, Gurjão and Pocinhos stations. However, the highest values calculated using the Z S test were identified more precisely in the central and northern regions of the catchment and refer to the Soledade, Coxixola, Caraúbas and Desterro stations. This spatial consistency can be explained by the substantial inhomogeneity of the rainfall regime in the catchment. Figure 8b shows the spatial distribution of the data resulting from Sen's method, which shows the spatial similarity of the data between the two methods. The highest values were observed at stations in the central, northern and eastern portions of the catchment. These results can be explained by the substantial variability of the rainfall index in the studied region.  Figure 9 shows the spatial distribution of the average rainfall erosivity for the wet and dry seasons, the average monthly erosivity and annual trends of rainfall erosivity. The results show that the erosivity values were within the range of 480-1200 MJ mm ha −1 h −1 yr −1 for the wet season (Fig. 9a) and 93-164 MJ mm ha −1 h −1 yr −1 for the dry season (Fig. 9b), i.e., a difference greater than 500% between the two seasons. In comparison, the erosivity values for the wet season increased from east to west and from west to east in the dry season. The results from this study show that erosivity in the region was low because of the intensity and amount of rainfall; however, many extreme events occur in the region that cause serious impacts to the environment such as land degradation. In the wet season, the average annual erosivity density was highest around the Prata rainfall station along the state border and at higher elevations. Erosivity density varied not only among the stations but also between months because these conditions change in the dry season and because the area around the Pocinhos rainfall station, located in the northeastern zone (as opposed to the Prata rainfall gauge, which is located in the southwestern zone), is the area with the highest erosivity values (164 MJ mm ha −1 h −1 yr −1 ). Figure 9c shows that the highest erosivity values occurred between February and April and that there is a greater than 70% concentration of erosivity in that period. These results suggest that there must be preventive measures during this period to protect the soil from  . 8 a Spatial distribution of the results obtained using the Mann-Kendall method for the rainfall stations used in this study and b spatialization of the results obtained by Sen's method for the rainfall stations used in this study erosion and avoid soil degradation (Fig. 9c). The results presented in this study regarding erosivity trends indicate that the trend is positive for four rainfall stations and is negative for nine stations, while the results showed no trends for the catchment study area in its entirety (Fig. 9d).

Conclusions
In this study, the Mann-Kendall trend and Sen's slope methods were used to investigate the temporal-spatial trends, rainfall variability and erosivity index for 13 rainfall stations in which each station provided data for over a 30-year period in the Epitácio Pessoa reservoir catchment. This work confirms the large spatial variability of trends in rainfall. Based on the results obtained, positive and negative trends were identified. The low spatial coherence of the trends found in the time series is consistent with the mixed pattern of positive Fig. 9 a Spatial distribution of mean rainfall erosivity for the wet season, b spatial distribution of mean rainfall erosivity for dry season, c mean monthly erosivity and d mean annual erosivity for the study area and negative changes at the global scale, especially those identified in northeastern Brazil. The Mann-Kendall and Sen's tests showed significant (downward) trends at the 15% significance level at the Camalaú, Barra de São Miguel and Prata stations. Exceptions were recorded in the other stations, whose trend tests were not statistically significant and in which changes in the time series occurred in different periods. The most noteworthy results observed were that (a) most of the calculated trend slopes are statistically nonsignificant and less than 25% of them are statistically significant at α = 0.15 and (b) the distribution of negative and positive trend slopes in the catchment is extremely irregular. The nonparametric Mann-Kendall test used in this study allowed a more comprehensive analysis of the climatic trends of rainfall time series, providing important information to help support more conclusive conclusions about the anomalies observed in some climatic parameters. Erosivity is the core driver of most soil erosion processes through its influence on particle detachment and overland flow in the region. There were no trends of erosivity and reduced rainfall in the eastern portion of the study area for the wet season, with a stronger core at the Camalaú station, while such trends occurred in the dry season in the western portion. These results need to be combined with information on land use in the catchment to identify possible causes of this variation. For example, if other methods were used to fill the gaps of missing data or rainfall data from other sources were used, the results may be affected.