Spatio-Temporal Dynamics and Physico-Hydrological Trends in Rainfall, Runoff and Land Use in Para í ba Watershed

: The detection and monitoring of changes in land use and land cover play a crucial role in understanding land degradation and are fundamental to preserving agroecosystems. Their association with hydrological information allows essential responses to changes in hydrological patterns to be identiﬁed, contributing to water security in watersheds. Therefore, this study aimed to assess spatio-temporal dynamics and physico-hydrological trends in rainfall, runoff and land use in the Para í ba watershed. The study was conducted in the Para í ba watershed, using land use data and information from pluviometric and ﬂuviometric stations with temporal series of more than 30 years. The Mann-Kendall statistical test was adopted to verify trends. Results indicate annual reduction trends for both native forest area and water bodies in the Para í ba watershed. On the other hand, the area designated for agriculture showed a signiﬁcant increase. The correlation analysis between water bodies and forests (R 2 = 0.63) highlights a strong association between the decrease in forest area and the reduction in water availability, inﬂuencing the decrease in annual ﬂow. These results serve as a warning to expand water resource management for the region, aiming to preserve and to enhance sustainable use. Therefore, the implementation of conservation measures, monitoring procedures, and adequate management is required to face the challenges imposed by climate change and land use and occupation, ensuring the water availability for the future.


Introduction
The Intergovernmental Panel on Climate Change (IPCC) has warned that the phenomenon of global warming will likely result in changes in the frequency of extreme weather events on a global scale [1].This process of climate change at a local level can be accelerated due to anthropogenic activities, such as the deforestation of native Caatinga vegetation, which can lead to a reduction of approximately 5% in precipitation at a local [2] and space scale [3], not only in Brazil, but also in several regions of the world, causing changes in hydrological dynamics, such as what happened in Jamaica, in the Caribbean region [4] in the arid district of Ananthapuramu, in the state of Andhra Pradesh, India [5], Peninsular Malaysia [6], Northeast Asia [7] and Nigeria [8].Indeed, modification in land utilization and land cover constitute primary drives for global land degradation.The identification and surveillance of these alterations have profoundly contributed to increase our comprehension of this practice and assisting in the management of natural resources in each region.Such knowledge is essential for the formulation of more efficient policies aimed at mitigating soil degradation and conserving the environment.By identifying and addressing the causes of land degradation, we can work towards sustainable land management practices and the conservation of natural resources [9].Furthermore, the identification of degradation factors is essential for the preservation of the agroecosystem and for making decisions aimed at mitigating actions, with a focus on environmental conservation and the enhancement of agricultural practices in the Brazilian Semiarid region [10].
The Caatinga biome has suffered a loss of approximately 37% of its native vegetation cover between 1985 and 2021 due to unregulated exploitation of natural resources for agricultural expansion and wood and charcoal production, according to data from the land use and land cover mapping project in Brazil-MapBiomas [11,12].
The spatio-temporal diagnosis of climatic variability and environmental changes in arid and semiarid regions is of paramount importance in mitigating pressures and impacts on natural ecosystems.In this context, the utilization of remote sensing techniques and satellite imagery has proven indispensable in detecting these changes and the processes of environmental degradation affecting vegetation cover and water resources through physical-hydrological parameters at the surface.These data have been instrumental in monitoring and decision-making concerning land use and natural resource management across Brazil [9].Droughts and climate variability can lead to increased poverty rates in a country due to changes in agroecological zones and land cover.To mitigate these effects on vegetation and ecosystem services, immediate measures are imperative.Such actions include implementing integrated watershed management, promoting climate-adapted agricultural practices, expanding forested areas, and conducting appropriate monitoring activities.Additionally, compliance with international environmental conventions and the use of spatio-temporal information to preserve a sustainable environment play a crucial role [13].
Precipitation is an essential component in hydrological processes, being directly associated to erosion, flooding, and aquifer recharge.In this sense, the spatio-temporal distribution of rainfall has a significant impact on the hydrological cycle, water security, agriculture, water resource management within the context of watersheds [3].In recent years, large alterations in both the rainfall regime have occurred across various areas within the Brazilian Semiarid region, which have implications for the management of natural resources and agricultural activities in the [9].
In the semiarid region of Paraíba State, the qualitative assessment of environmental conflicts and water resource utilization reveals that intensive agricultural and livestock activities have played a significant role in modifying the water dynamics, particularly regarding the intensive extraction of natural resources, resulting in a myriad of environmental problems and conflicts related to water use [14].
In this context, the temporal evolution of land use and water dynamics in the region becomes essential, especially in the context of the river basin [14].Some research addresses analysis of the spatio [3,6,15] and temporal dynamics of rainfall [3,5,13] and water flows [16], as well as land use and cover [11,17,18].However, studies that integrate these processes are still limited in the Semiarid region of Brazil.Like this, the aim of this research was to assess the spatio-temporal dynamics and physical-hydrological trends of rainfall, runoff, and land use in the Paraiba watershed, in the context of climate change and human activities.

Study Area
The study was conducted in the Paraiba watershed, located between latitudes 8 • 0 00 and 8 • 45 00 S and longitudes 37 • 45 00 and 38 • 45 00 W. This basin covers an area of 20,123.4km 2 , ranging from the Sertão region to the Zona da Mata, and discharging into the sea in the state of Paraíba.For the research, five pluviometric and fluviometric stations with temporal series of more than 30 years were selected (Figure 1).

Study Area
The study was conducted in the Paraiba watershed, located between latitudes 8°0′00″ and 8°45′00″ S and longitudes 37°45′00″ and 38°45′00″ W. This basin covers an area of 20,123.4km 2 , ranging from the Sertão region to the Zona da Mata, and discharging into the sea in the state of Paraíba.For the research, five pluviometric and fluviometric stations with temporal series of more than 30 years were selected (Figure 1).

Open Access Fluviometric and Pluviometric Data
The fluviometric and pluviometric data were collected through the Hidroweb website, a tool that is part of the National Information System on Water Resources, managed by the National Water Agency.This application allows access to various telemetric data collected by the National Hydrometeorological Network, such as information on rainfall and river flows, using the codes of the observed monitoring stations [19].To assist in the download and processing of the data, the hydrobr package was used [20].

Satellite-Derived Data Available (MapBiomas Brazil)
The land use data used in the MapBiomas Brazil project were obtained from the Landsat satellite series and processed in Google Earth Engine (GEE) environment and were used to create thematic maps of land use and land cover.This spatio-temporal monitoring of the natural environment of Brazilian biomes was driven by the implementation of a low-cost and open-access methodology.
MapBiomas provides monitoring of native vegetation areas, including mapping of loss and/or resilience of natural vegetation cover in biomes.Additionally, it also covers monitoring of pasture areas, assessing the quality and degradation of these areas, as well

Open Access Fluviometric and Pluviometric Data
The fluviometric and pluviometric data were collected through the Hidroweb website, a tool that is part of the National Information System on Water Resources, managed by the National Water Agency.This application allows access to various telemetric data collected by the National Hydrometeorological Network, such as information on rainfall and river flows, using the codes of the observed monitoring stations [19].To assist in the download and processing of the data, the hydrobr package was used [20].

Satellite-Derived Data Available (MapBiomas Brazil)
The land use data used in the MapBiomas Brazil project were obtained from the Landsat satellite series and processed in Google Earth Engine (GEE) environment and were used to create thematic maps of land use and land cover.This spatio-temporal monitoring of the natural environment of Brazilian biomes was driven by the implementation of a low-cost and open-access methodology.
MapBiomas provides monitoring of native vegetation areas, including mapping of loss and/or resilience of natural vegetation cover in biomes.Additionally, it also covers monitoring of pasture areas, assessing the quality and degradation of these areas, as well as monitoring the evolution of irrigation and mapping of different agricultural crops.Other monitored aspects include water resources, such as the water supply condition of rivers, streams, lakes, lagoons, and strategic reservoirs, as well as urban infrastructure [9] and application for detecting dry and rainy seasons and identifying vulnerable areas [21].

Mann-Kendall Test
The trend analysis was performed for the annual precipitation series, total annual runoff, and land use (forested areas, farming (agriculture and pasture), and water) using the non-parametric statistical test of Mann-Kendall [22,23].The test statistic (S) is described by Equation ( 1): where, "n" is the number of data points, "xj" and "xi" represent each of the measurements at different time steps i and j, with i = j; and "sgn(xj − xi)" is defined by Equation ( 2): If the dataset is uniformly and independently distributed, then the mean of S is zero, and the variance of S is determined by Equation (3): where, n-is the dataset number; t-is the number of data points with repeated values in a given group; q-is the number of groups containing repeated values.
In a long temporal series, the statistical value "S" can be represented by "Z" according to the following conditions of the Equation (4): When −1.96 ≤ Z ≤ 1.96, the null hypothesis (H0) is accepted, indicating that there is no statistically significant trend in the temporal series.The trend is considered significant at a 90% confidence level when |Z| > 1.64, at a 95% confidence level when |Z| > 1.96, and at a 99% confidence level when |Z| > 2.58.A positive value of Z indicates an increasing trend in the sequence, while a negative value of Z reflects a decreasing trend.

Statistical and Multivariate Data Analysis
An additional statistical procedure employed was Principal Component Analysis (PCA), used to assess the relationship between soil and rainfall, runoff and land use in the watershed aspect [24][25][26].PCA employs an orthogonal transformation to convert a set of potentially correlated observations into a set of linearly uncorrelated variables.This transformation is defined in such a way that the first principal component has the highest possible variance, and each subsequent component, in order, has the highest possible variance while being orthogonal to the previous components [27].The covariance matrix is assessed, allowing for the calculation of eigenvectors and eigenvalues associated with the matrix.Eigenvalues are computed to determine the contribution of each variable to the total variance of the principal components.A high contribution value indicates that the variable is well-represented in the principal component [15].PCA generates two factors, with the first one representing the primary component (Dimension 1 or Dim1).The second factor is orthogonal to the first and is denoted as dimension 2 (Dim2).Each dimension or component is the result of the individual contributions of each of the original variables [28].

Temporal Analysis of Land Use in the Paraíba Watershed
The analysis results indicate statistically significant trends (95% confidence level) in the annual reduction of native forest area (Z = −2.23)(Figure 2a) and water bodies (Z = −2.89)(Figure 2c) in the Paraíba watershed.On the other hand, the agricultural area showed a significant increase at a 90% confidence level (Z = 1.89) (Figure 2b), indicating an expansion of agricultural activity that may be impacting the hydrodynamic dynamics of the region.The correlation analysis between water bodies and forests reveals a determination coefficient of 0.63 (Figure 2d), justifying that the decrease in forest area is also associated with a reduction in water availability.
dimension or component is the result of the individual contributions of each of the original variables [28].

Temporal Analysis of Land Use in the Paraíba Watershed
The analysis results indicate statistically significant trends (95% confidence level) in the annual reduction of native forest area (Z = −2.23)(Figure 2a) and water bodies (Z = −2.89)(Figure 2c) in the Paraíba watershed.On the other hand, the agricultural area showed a significant increase at a 90% confidence level (Z = 1.89) (Figure 2b), indicating an expansion of agricultural activity that may be impacting the hydrodynamic dynamics of the region.The correlation analysis between water bodies and forests reveals a determination coefficient of 0.63 (Figure 2d), justifying that the decrease in forest area is also associated with a reduction in water availability.By adopting a moving average of the first 5 years of the series (1985-1989) for land use data and comparing it with the total area of the Paraíba watershed (20,123.4km 2 ), in addition to comparing it with the last 5 years of the series (2017-2021), it was observed that there was a 12% reduction in forest area compared to the initial average area of the first 5 years (11,428.2km 2 ).Simultaneously, there was a 14% increase in agricultural area By adopting a moving average of the first 5 years of the series (1985-1989) for land use data and comparing it with the total area of the Paraíba watershed (20,123.4km 2 ), in addition to comparing it with the last 5 years of the series (2017-2021), it was observed that there was a 12% reduction in forest area compared to the initial average area of the first 5 years (11,428.2km 2 ).Simultaneously, there was a 14% increase in agricultural area (8133.1 km 2 ), while the area of water bodies showed a 51% reduction (184.1 km 2 ) (Figure 2).These findings indicate that the Paraíba watershed has undergone significant changes in terms of vegetation cover and land use over the 36 years.
The reduction in the water bodies area may reflect environmental degradation in the basin, as clearly observed in the spatio-temporal land use maps (Figure 3).The reduction in the water bodies area may reflect environmental degradation in basin, as clearly observed in the spatio-temporal land use maps (Figure 3).

Temporal Analysis of Rainfall
It was observed that extensive vegetation suppression has occurred to the northw southwest, and east of the Epitácio Pessoa Reservoir.Despite the significant reduct (51%) (Figure 2) in water bodies observed in recent years, among the evaluated poi only point 736022 showed an increasing trend in rainfall (Z = 1.73) at a 90% significa

Temporal Analysis of Rainfall
It was observed that extensive vegetation suppression has occurred to the northwest, southwest, and east of the Epitácio Pessoa Reservoir.Despite the significant reduction (51%) (Figure 2) in water bodies observed in recent years, among the evaluated points, only point 736022 showed an increasing trend in rainfall (Z = 1.73) at a 90% significance level.In the remaining points, no significant change was identified in the studied hydrological series (Figure 4).
Geographies 2023, 3, FOR PEER REVIEW 7 level.In the remaining points, no significant change was identified in the studied hydrological series (Figure 4).

Temporal Analysis of Runoff Stations in the Paraíba Watershed
Only 38830000 station did not show a decreasing trend in runoff according to the Mann-Kendall test.The other stations (38850000, 38860000, 38880000, and 38895000) exhibited trends at a 99% confidence level, with values of (−2.67, −3.73, −4.85, −2.86 m 3 per year), respectively (Figure 5).

Temporal Analysis of Runoff Stations in the Paraíba Watershed
Only 38830000 station did not show a decreasing trend in runoff according to the Mann-Kendall test.The other stations (38850000, 38860000, 38880000, and 38895000) exhibited trends at a 99% confidence level, with values of (−2.67, −3.73, −4.85, −2.86 m 3 per year), respectively (Figure 5).
Upon examining the magnitude of the trend, it was observed that there was a reduction of 54%, 61%, 64%, and 59% in the annual flow at the runoff stations 38850000, 38860000, 38880000, and 38895000, respectively.Upon examining the magnitude of the trend, it was observed that there was a reduction of 54%, 61%, 64%, and 59% in the annual flow at the runoff stations 38850000, 38860000, 38880000, and 38895000, respectively.

Multivariate Analysis of Rainfall, Runoff and Land Use
The principal components, Dim1 and Dim2, were responsible for explaining 71.8% of the total variability, with a special emphasis on Dim1, which contributed 57.7% of the total data variability.In this component, we notice that variables related to runoff are positioned on the positive side of variables such as precipitation, surface water, and forest cover.On the other hand, farming is opposed to forest cover (Figure 6a) and negative correlation of farming with forest cover, rainfall, and runoff (Figure 6b).This suggests that the reduction in forest area is associated with the expansion of farming, especially in the southwestern and eastern regions of the basin (Figure 3), implying a trend of reduced runoff in four out of the five monitoring points evaluated in the Paraíba watershed (Figure 5).

Multivariate Analysis of Rainfall, Runoff and Land Use
The principal components, Dim1 and Dim2, were responsible for explaining 71.8% of the total variability, with a special emphasis on Dim1, which contributed 57.7% of the total data variability.In this component, we notice that variables related to runoff are positioned on the positive side of variables such as precipitation, surface water, and forest cover.On the other hand, farming is opposed to forest cover (Figure 6a) and negative correlation of farming with forest cover, rainfall, and runoff (Figure 6b).This suggests that the reduction in forest area is associated with the expansion of farming, especially in the southwestern and eastern regions of the basin (Figure 3), implying a trend of reduced runoff in four out of the five monitoring points evaluated in the Paraíba watershed (Figure 5).
Forest had a significant positive impact on runoff (9%), demonstrating the importance of Caatinga vegetation in conserving water resources and influencing the hydrological dynamics of the watershed (Figure 6a).
The Spearman correlation reveals that farming exerts a negative influence on both rainfall and runoff, strengthening the considerations obtained in the PCA analysis, with emphasis on the significant negative correlation with forest cover, reaching a coefficient of −0.98 (Figure 6b).
The pluviometric station 735015 (RainD) and the fluviometric station 38880000 (FlowD) showed the most solid relationship, with smaller angles between the vectors.It is relevant to observe that in this region, changes in land cover over time were less variable.On the other hand, the vectors of the pluviometric station 736020 (RainA) and the fluviomet-ric station 38830000 (FlowA) showed a greater separation (larger angles) between them, which is associated with the area that suffered the greatest negative impact in reducing Caatinga vegetation.Forest had a significant positive impact on runoff (9%), demonstrating the importance of Caatinga vegetation in conserving water resources and influencing the hydrological dynamics of the watershed (Figure 6a).
The Spearman correlation reveals that farming exerts a negative influence on both rainfall and runoff, strengthening the considerations obtained in the PCA analysis, with emphasis on the significant negative correlation with forest cover, reaching a coefficient of −0.98 (Figure 6b).
The pluviometric station 735015 (RainD) and the fluviometric station 38880000 (FlowD) showed the most solid relationship, with smaller angles between the vectors.It is relevant to observe that in this region, changes in land cover over time were less variable.On the other hand, the vectors of the pluviometric station 736020 (RainA) and the fluviometric station 38830000 (FlowA) showed a greater separation (larger angles) between them, which is associated with the area that suffered the greatest negative impact in reducing Caatinga vegetation.

Temporal Analysis of Land Use
The annual reduction in native forest, water bodies and a significant increase in agricultural area in the Paraíba river basin (Figure 2) are results that are crucial for understanding environmental changes in the region and can provide support for the implementation of conservation and sustainable natural resource management measures.
The reduction in forest area and the increase in agricultural area may be related to agricultural expansion in the region.
In the Semiarid region of Africa, vegetation degradation has been found to be correlated with increased drought conditions in the region.Therefore, it is recommended that measures be immediately implemented to mitigate the adverse impact of climate variability and worsening aridity on vegetation and ecosystem services [13].This study corroborates our own research findings as well as those presented by [30], who sought to discern

Temporal Analysis of Land Use
The annual reduction in native forest, water bodies and a significant increase in agricultural area in the Paraíba river basin (Figure 2) are results that are crucial for understanding environmental changes in the region and can provide support for the implementation of conservation and sustainable natural resource management measures.
The reduction in forest area and the increase in agricultural area may be related to agricultural expansion in the region.
In the Semiarid region of Africa, vegetation degradation has been found to be correlated with increased drought conditions in the region.Therefore, it is recommended that measures be immediately implemented to mitigate the adverse impact of climate variability and worsening aridity on vegetation and ecosystem services [13].This study corroborates our own research findings as well as those presented by [30], who sought to discern the historical patterns of drought in the Greater Horn of Africa and recorded a pervasive tendency towards escalating drought.
The spatio-temporal maps of land use show a reduction in water bodies (Figure 3).This pattern of substantial land use change has been highlighted by [9] for Campina Grande, PB municipality.Over time, there has been a reduction in forested areas, as indicated by the spectral condition of Caatinga vegetation cover.The significant reduction of natural vegetation in the semiarid region is evidenced by the expansion of agricultural areas (cropland and pasture) and urban infrastructure over the last three decades.
According to [31], the low percentage of Caatinga vegetation cover and the predominance of farming are correlated with animal density by municipalities, indicating the strong influence of this activity on vegetation cover reduction.[18] identified a similar pattern of land use change and expansion of farming areas, listing the main consequences of these land cover alterations as erosion and soil degradation, rural-urban migration, fragmentation of farming lands, climate change, and reduced crop yields.
Due to the considerable portion of the territory of the Paraíba River Basin located in a semi-arid environment, irrigation appears as an extremely important technology to increase the production and diversification of agricultural crops.The agricultural areas identified and mapped correspond to 21.42% of productive establishments [32].
The main agricultural crops grown in the region include tomatoes, peppers, beans, cabbage, lettuce and onions, both during and between harvests, as well as bananas and coconuts, using drip, furrow and micro-sprinkler irrigation systems.It is important to highlight the significant presence of these activities in the vicinity of the Epitácio Pessoa Dam [33].
According to MapBiomas [12], through land cover and land use monitoring, the data indicates that the São Francisco watershed experienced a significant loss of its natural water surface between 1985 and 2020.According to the results, the loss was 50%, but it's important to note that human actions, such as artificially increasing water surface in reservoirs, contributed to a 13% increase in water surface.This means that the actual reduction was 4%, with the largest losses occurring in the Upper and Lower São Francisco regions, with reductions of 19% and 21%, respectively.

Temporal Analysis of Rainfall
In the Brazilian Northeast, trends of both reduced and increased rainfall, as well as a decrease in rainy days, have been observed [3].In our study (Figure 4), although no significant reduction was observed at the evaluated points, a negative slope of the trend line was observed at 4 stations (736020, 735124, 735015, and 735036).It is crucial to continue monitoring in this regard, not only concerning rainfall but also regarding reservoirs.
Regarding the Epitácio Pessoa Reservoir, one of the most important reservoirs in the Paraíba watershed, an evaluation of the relationship between rainfall from 1972 to 2016 in the reservoir's contributing basin and its volume revealed that historically, its recharge is concentrated more in the first months of the year during the rainy season.Its volume exhibits high variability and is conditioned by multiple water uses from the reservoir, particularly in cities that demand substantial volumes [34].
In order to prevent the collapse of Epitácio Pessoa Reservoir due to the arrival of dead-level volume, as a result of drought and the high multiple use of the reservoir's water, which is responsible for supplying the important city of Campina Grande and the neighboring cities, the contribution of waters from the São Francisco River Integration Project with the Northern Northeast Basins (PISF) has improved the quantitative and qualitative water security in relation to physical and chemical parameters, particularly in relation to ammonia, EC, chlorides, hardness, sulfate and STD [16,34,35].
Correia et al. [36] identified eleven municipalities within the contributing basin of the Epitácio Pessoa Reservoir where land use is predominantly oriented towards agricultural and livestock activities, resulting in low coverage of native vegetation.Due to these agricultural and livestock practices, soil leaching, and erosion occur, directly impacting the river channel, leading to siltation and eutrophication of the water body.Additionally, the development of cyanobacteria is also observed because of these activities, contributing to the degradation of the aquatic ecosystem.These environmental issues underscore the need for the implementation of sustainable management practices and conservation measures to protect the natural resources of the region.

Temporal Analysis of Runoff Stations
The high reductions in annual runoff (Figure 5) indicate that the Paraíba watershed has undergone significant changes in terms of vegetation cover, land use, and rainfall variability.This scenario may have negative implications for water availability in the region, especially during extended drought periods.Similarly, Fernandes et al. [37] found that changes to land use coverage within the Sao Francisco River Basin had stark financial implications.This research revealed a loss of around US$3.2 billion from reduced carbon sequestration and stocks, as well as disruptions to water usage dynamics.
The modeling of future scenarios in the context of watersheds raises concerns about the intense land use over time, and such climate changes have a direct impact on the water resources scarcity [38].
The collected and analyzed information is of fundamental importance for hydrological modeling.In this sense, several models have the potential to provide solid answers and contribute to decision-making related to water resource management, considering changes in land use and incorporating climate trends [39].It is worth noting that precipitation, as the main source of water supply for the watershed and responsible for the balance of the water budget, must be thoroughly evaluated in future simulations to support decisionmaking [40].

Multivariate Analysis of Rainfall, Runoff and Land Use
The Caatinga forest has a significant impact on runoff regulation (Figure 6), in accordance with the study of [2], which indicated that deforested areas led to a reduction of approximately 22 mm of rainfall per year compared to areas covered by native Caatinga vegetation, representing a decrease of about 5%.
Principal Component Analysis (PCA) showed potential for hydrological assessment, exploring the associations between precipitation and runoff, in the Darong River watershed, situated in Guangxi province, China [26].Additionally, it was utilized to analyze the relationships among precipitation events, surface runoff, and soil erosion on the Pisha sandstone slopes of the Loess Plateau in China, where it has been demonstrated that forests play a pivotal role in erosion control [25].Therefore, PCA offers revelant insights for the modelling and management of stormwater within watersheds [24].
In the São Francisco watershed, the runoff water monitoring in the Verde Grande Sub-Basin over time shows a reduction in runoff of 80.5% from 1950 to 1976 and 69.2% from 1979 to 2000 [41].An analysis of potential droughts in prospective scenarios up to 2030 has also suggested that the reservoir volume would decrease, potentially reaching critical levels by that point in time [42].
The significant trend of decreasing vegetation cover and runoff rates at various monitoring points in the Paraíba River Watershed, along with the negative correlation between runoff and the increase in degraded areas, demands proactive attention from administrative and governmental institutions in the state.It is imperative to intensify initiatives for the restoration of degraded areas, with particular emphasis on preserving riparian forests.As highlighted by [43], riparian forests play a fundamental role as a natural barrier for water runoff, attenuating its direct impact on the soil and preventing erosion and the deposition of solid particles in riverbeds, thus avoiding sedimentation.They play a crucial role in maintaining environmental balance.
The significant changes in the landscape, especially due to deforestation, have consequences for hydrological behavior, triggering erosional processes, sedimentation of water bodies, and a reduction in their runoff rates, among other associated phenomena.These changes have a negative impact on the local ecosystem, adversely affecting the quality and availability of water resources and compromising environmental sustainability [17].

Conclusions
From 1985 to 2020, the Paraíba watershed underwent significant changes in land use, with a notable reduction in native vegetation cover and a decrease in water bodies.These alterations have resulted in substantial reductions in annual runoff along the main surface flow channel.This scenario serves as a warning about the importance of expanding water resource management in the region to preserve and sustainably utilize the available natural resources.
The spatio-temporal monitoring of land use has enabled the identification of the most vulnerable areas in the Paraíba watershed.Therefore, there is an urgent need to adopt appropriate conservation, monitoring, and management measures to address the challenges posed by changes in land and climate use, to ensure the sustainable availability of water and soil.

Figure 2 .
Figure 2. Mann-Kendall Trend for forest areas (a), agriculture (b), water bodies (c), and correlation between forest area and water body area (d) in the Paraíba watershed.

Figure 2 .
Figure 2. Mann-Kendall Trend for forest areas (a), agriculture (b), water bodies (c), and correlation between forest area and water body area (d) in the Paraíba watershed.

(
8133.1 km 2 ), while the area of water bodies showed a 51% reduction (184.1 km 2 ) (Fig2).These findings indicate that the Paraíba watershed has undergone significant chan in terms of vegetation cover and land use over the 36 years.

Figure 3 .
Figure 3. Land use for forested areas, farming (agriculture and pasture), and water bodies in Paraíba watershed.

Figure 3 .
Figure 3. Land use for forested areas, farming (agriculture and pasture), and water bodies in the Paraíba watershed.