Inﬂuence of Land Use Changes on the Longaví Catchment Hydrology in South-Center Chile

: During recent decades, the South-Central part of Chile has shown strong vulnerability due to the effects of land use change (LUC). The interaction of these changes with local hydrology has not been adequately investigated and is poorly understood, especially in mountainous areas under irrigated agriculture. We applied the SWAT + agrohydrological model to study the effects of LUC on hydrological ﬂuxes in the Longaví catchment, Maule region, South-Central Chile. Land use maps (LUMs) from 1997, 2009, and 2016 were used in conjunction with a 41-year (1979–2019) hydro-meteorological series of daily observations as forcing data. The dominant changes in land use during the study period relate to agriculture, shrublands, forestry of exotic species, and urban sprawl. First, the LUM of 1997 was used for model setup, sensitivity analysis, calibration, and validation. Second, the impact of LUC documented through LUMs 2009 and 2016 was analyzed. Our analysis clearly reveals that the overall water balance and internal moisture redistribution in the Longaví catchment have been considerably affected by decreases in precipitation, changes in land use and water use practices. Unless a comprehensive regulatory system is introduced that addresses current climatic conditions and territorial use, it is likely that the decrease in water resources will persist and worsen through climate changes.


Introduction
Human-induced processes and natural system changes can cause alterations in the landscape with corresponding environmental impacts on both natural and human systems [1,2]. In particular, land use changes (LUC) attributable to human activities can have potential impacts on catchment hydrology [3], and as a result on water resources availability [4]. At basin scale, particular land management practices and alterations in the use of the territory can affect the partitioning and redistribution of water between various flow pathways and system components [5]. Affected processes include precipitation interception by the canopy, where the type of vegetation can influence the interception rate (i.e., native or exotic species) [6,7]. In addition to evapotranspiration, infiltration, runoff generation, and surface flow pattern [8], as well as soil erosion [9].
In Chile, a country that is known for its important mining sector, agricultural crop production, and forestry [10], anthropogenic changes in land use are a very common practice and are justified at the institutional level by the need to support the socio-economic development of specific regions [11]. These types of land use changes are intertwined with the local climate, the productive capacity of particular areas, soil conditions, and the possibility of fostering exotic tree plantations (i.e., Pinus radiata and Eucalyptus globulus) [10]. In addition, some policy management allowed the transition from natural forest cover to other types of land use, including exotic tree plantations, intensive irrigated agriculture, and urban development, particularly from 1994 to 2014 [12][13][14]. For instance, 67% of Chilean temperate forests were lost by deforestation between 1975 and 2000 [15]. As a result, a large part of the country has undergone rapid land use change, with large areas that were originally covered by native vegetation becoming the subject of commercial plantation forestry [12,16] and intensive agriculture exploitation [11,15].
To satisfy agricultural demand, the irrigation water use is estimated to be around 77% to 85% of the total available freshwater in the country [17,18]. However, the overall lack of monitoring, among others due to scarcity of observing stations, hampers the quantification and estimation of available water, including surface flow and groundwater extractions [17,19]. In addition, an annual rainfall deficit from 25 to 45% has been reported since 2010, increasing water scarcity in the center of the country [18]. As consequence, the lack of an efficient water use governance leads to conflicts related to water overuse [20,21], resulting in unequal redistribution among water users affecting the basin hydrology.
To understand and evaluate the impacts of LUC on runoff and water availability, hydrological models constitute important and widely used investigation tools [22,23]. Particularly physical-based hydrological models allow one to describe the complex interaction between LUC and various components of the hydrological cycle in detail [24,25]. Continuous-time models, such as the Soil and Water Assessment Tool + (SWAT+) based on the original version of SWAT [26,27], support the evaluation and quantification of the impacts of different land management practices on water resources with varying soils and land use over long periods of time [28]. Thanks to the possibility to assess specific land management practices in agricultural systems and forest production, among many other processes, SWAT has been widely applied in different watersheds around the world [29][30][31].
The application of SWAT in Chile is not new. Some studies have been carried out using SWAT as a hydrological model with the aim of assessing the impacts of land use changes in Chile. These applications remain limited to mountain areas in proximity to the coast [32][33][34][35] or to the analysis of the impacts of climate change on snow accumulation [36,37]. Thus, the applicability of SWAT+ related to land use change in the rugged terrain of South-Central Chile has not yet been investigated, particularly in data-poor regions and/or irrigated agricultural areas.
Our chosen study area, the Longaví catchment, is a sub-entity of the Loncomilla river, located in the mountainous area of South-Central Chile. The main human activities in the area are related to agriculture and exotic tree plantations, established especially during the last few decades. Unfortunately, the impact of such land use management practices is not well understood, particularly with regard to the effects on specific hydrological processes, soil erosion, and overall basin water balance. Although it has become evident that the flow rate of rivers and precipitation has decreased, there is no information on the impacts of both precipitation and land use change effects on the catchment hydrology. The aim of this study is, therefore, a process-based analysis to evaluate how the combined effect of changes in precipitation and land use changes in the Longaví catchment has affected hydrology over a total period of 30 years. Given the similarity of territory and land management practices, our results can be extrapolated to similar systems in the region.

Materials and Methods
The Longaví catchment is located in the Maule region, VII Region of Chile, between latitude 35 • 49 S and longitude 71 • 47 W. The watershed encompasses an area of 1387 km 2 ( Figure 1). The geology of the basin is characterized by the presence of intrusive rocks, sedimentary, volcano, and volcano-sedimentary sequences [38], and elevations ranging from 102 m to 3174 m above sea level. The river basin is dominated by temperate-Mediterranean climate and annual total precipitations of around 1669 mm year −1 , with a dry summer of 74 mm year −1 (January, February, March), a rainy autumn of 746 mm year −1 (April, May, June) and winter of 682 mm year −1 (July, August, September), with regular snow accumulation in the mountain areas, and a variable spring of 746 mm year −1 (October, November, December) [39]. Concerning temperatures, maximum mean values of 24.6 • C for the summer season and minimum mean values of 0.5 • C for the winter season have been recorded in the study area over 41 years of continuous hydro-meteorological observation (temperatures and precipitation) [39]. With reference to existing local water management infrastructure, the Bullileo Dam operated by "Junta de Vigilancia del Río Longaví y sus Afluentes" (Longaví river surveillance board) stores water for irrigation mainly during drought periods. In addition, the Digua channel transfers water from the Longaví river to the Digua Dam, which is situated outside the watershed boundaries.

Swat+ Model Description
The Soil and Water Assessment Tool + (SWAT+) is a semi-distributed hydrological model that has been especially developed for agro-hydrological impact studies on management practices and climate-water interactions, as well as the transport of nutrients and pesticides on a basin scale [28]. SWAT + reproduces most of the physical processes of the hydrological cycle in different time steps on the basis of multiple data inputs, such as hydro-meteorological forcing, topography, soil parameters, land use, and land management information, including crop rotation [28,40]. With this input information, the hydrologic cycle can be simulated based on the basis of 1-dimensional water mass balance and steadystate momentum conservation across adjacent vertical soil columns (Equation (1)) [26,31]. The water balance equation for a soil column cell is stated as follows [26]: where SW t is the soil water content (mm) in time step t, SW o is the initial soil water content (mm), Pr is the amount of precipitation on day i (mm), Surq is the amount of surface runoff on day i (mm) and is estimated using the Soil Conservation Service (SCS) curve number equation, Latq is the amount of lateral subsurface flow to the channel on day i (mm), and it is originated from the saturated zone of soil layers and contributes to the stream flow calculated in each layer by a kinematic storage model, ET is the amount of actual evapotranspiration on day i (mm), Perc is the amount of percolation of soil water from the bottom of the soil profile on day i (mm), and is calculated as the sum of three terms: Flo (Return Flow), Rchrg (Recharge to deep aquifer) and Revap (Plant water uptake and evaporation). Additionally, the sum of Surq and Latq can be expressed as WY (Water yield).

Model Setup
In this study, SWAT+ (v. 2.0.4) was interfaced with QSWAT+ (v. 2.0. 6), an opensource graphical user interface [41], and the sub-basin scheme derived from the digital elevation model (DEM). To represent the river and natural flow paths, a topographic analysis is performed, leading to a drainage direction map with respective downhill cell interconnectivity. The headwaters streams and the river network were defined considering the point of union with the Loncomilla River and a threshold area ( Figure 1). The final network extent was verified against the network observable from remotely sensed images. To complete the model, Bullileo Dam was added as a storage entity to the network structure with its specific area and volume information.
Hydrological Response Units (HRUs) as basic spatial modeling entities were generated on the basis of raster images and topographic analysis by merging slope maps from the DEM, soil type, and land use maps. Meteorological forcing data at the daily time step for the hydrological processes have been pre-processed and readied for model use. For the calculation of the potential evapotranspiration (PET), the elaborate Penman-Monteith formulation was used. The input parameters are presented in Table 1. The DEM was obtained from the Shuttle Radar Topography Mission (STRM) dataset with 90 × 90 m spatial resolution [42]. The slopes of the land surface were subdivided into six categories ( Figure 2a) based on topographic information extracted from the DEM. The data on soil properties required for completion of the input database of the model were acquired during a field campaign in the catchment. The field data set was complemented with soil properties data from agricultural science studies of the VIIth administrative region that were conducted by the Natural Resources Information Center (CIREN) in 1997. In summary, a comprehensive soil type map for the study area was generated, which includes the 21 soil types described by CIREN [43] (Figure 2b).
In addition, soil samples under different land use types were taken in the field to corroborate and complement the soil physicochemical information provided by CIREN ( Table 2). Soil data like bulk density (BD), soil carbon content (CBN), saturated hydraulic conductivity (K), pH, and soil texture were determined through fieldwork in the framework of this study. The analysis of the soil samples was carried out using the Laboratory of Soils and Foliar Analysis of Pontificia Universidad Católica de Valparaíso (PUCV). According to CIREN [43], the predominant soil textures in the catchment are loamy sand soils (71.5%) developed from basic volcanic materials, loamy silt soils (10.5%) originated from volcanic ashes and sedimentary deposits, and loam soils (10.2%) with an alluvial origin and sediment deposits.  The hydro-meteorological input data for the years 1979 to 2019 were acquired from the Catchment Attributes and Meteorology for Large-sample Studies-Chile dataset (CAMELS-CL), taking as a reference 10 climate stations for precipitation and temperatures (daily minimum and maximum) located in the surrounding area of the study area ( Figure 2c, Table 1). Due to scarcity of in-situ data in Longaví, time series of wind velocity and relative air humidity from the nearby stations of the mountain area of Mataquito river basin, Maule region, were used. These data were supplied by "Dirección General de Agua" (DGA), the Central Water Directorate. Solar radiation data were obtained from "Explorador Solar" of the Chilean Ministry of Energy [45].
Three discharge stations are operated in the Longaví river network. Rio Longaví en Longitudinal (Longitudinal), located between 36 • 00 S latitude and 71 • 43 W longitude, has only a very short data record with few months of discharge measurements covering the period 1979 to 1985. Because of the limited usability of those records, two additional discharge stations were used instead with longer records: Quiriquina station, located between 36 • 23 S latitude and 71 • 46 W longitude, and El Castillo station, located between 36 • 15 S latitude and 71 • 20 W longitude [39]. Both the Quiriquina and El Castillo stations have discharge data available that cover the 1979 to 2019 period. Nevertheless, the time series are partially incomplete, with some gaps and unreliable discharge peak values.

Land Use Maps and Land Use Management in SWAT+
The land use (LU) maps of 1997, 2009, and 2016 were obtained from the National Forest Corporation (CONAF) [44], following the LU classifications into 12 categories ( Figure 3). The main land use types in the basin (Table 3) were agriculture with 416.8 km 2 (LU 1997), in addition to mixed deciduous and evergreen forests that covered an area of 414.7 km 2 (LU 2009) and 415.5 km 2 (LU 2016). For the studied periods, the most dominant land use types were agriculture (>28%), followed by forest mixed (>27%) and pine tree plantations (>9%). Between the years 1997 and 2016, the most important increase took place with mixed forests and pine tree plantations, with values that range between 2.67% (37.0 km 2 ) and 1.70% (23.7 km 2 ) respectively, followed by urban developments with 0.55% (7.6 km 2 ) of incremental extension. The highest areal reduction was recorded for shrublands with −2.50% (−34.7 km 2 ) and agriculture with a reduction of −1.42% (−19.7 km 2 ), followed by evergreen forests with −0.23% (−3.2 km 2 ).
The management schedules for agriculture, forest mixed and pine tree plantations were implemented in the model as pre-defined operations using conditions and actions inserted into the SWAT+ decision tables system [40]. Following the decision tables logic for automatic management schedules, the conditions of "plant" and "harvest" every 20 years were implemented for forest mixed. In the case of pine tree plantations, an automatic management schedule with "plant" and "harvest and kill" every 20 years were assigned in dependence on the pine tree maturity level. For agriculture, an automatic management schedule with "plant", "harvest and kill", "rotation and reset" for every year was established, taking as a reference local information from the traditional crops in the studied area. In addition, an auto-application of the furrow irrigation system was assigned under soil-water stress and plant maturity status, whereby the channel is used as a source to satisfy irrigation water demand.

Model Sensitivity Analysis, Calibration and Validation
A 41-year period (1979 to 2019) was selected in accordance with available forcing records for simulating runoff processes on a daily time step basis in Longaví, reserving one full year for model warm-up under the land use conditions of 1997. El Castillo and Quiriquina discharge stations were selected for model sensitivity analysis and multi-site model calibration from 1985 to 1986. In addition, El Castillo and Quiriquina stations were used jointly for validation purposes from 1996 to 1998 period, considering the availability and continuity of the data during the study period. A longitudinal station was used to evaluate the baseline simulation output in the lower part of the catchment, mainly to verify the flow after satisfying irrigation demand further upstream.
Global sensitivity analysis and calibration procedures were performed using the tool SWAT+ Toolbox v0.7.6, a sensitivity and automatic calibration module developed with basis on the existing calibration tool IPEAT+ [46]. Following indications from the literature [46,47], a total group of 16 hydrological parameters was selected with the goal of identifying the most sensitive ones for the studied system (Table 4). After identifying the most sensitive parameters, SWAT+ Toolbox was set up and run with 1000 iterations for the calibration period (1985)(1986). After successful calibration, the model was validated over a three-year period (1996-1998) by benchmarking the output against observations recorded at El Castillo and Quiriquina stations in the upper and middle part of the catchment. In this context, usual efficiency criteria were adopted as model performance metrics. For additional verification, cumulative mass analysis on different hydrological fluxes was used. To validate model performance for the entire catchment upstream of the Longaví basin outlet, a Longitudinal station was used to evaluate the irrigation system under the reference scenario. The calibration of the irrigation system was tested under different operating conditions referring to the longitudinal station discharges as a benchmark record. Afterward, a manual calibration of the irrigation system was performed by considering the water demand of the plants under different irrigation schedules.
To examine the effects of past land use changes on the catchment hydrology, the model was set up with land use maps LU 2009 and LU 2016 without any further modification under the same period of hydro-meteorological data (1979-2019). The output consistency of the model was verified for the periods of 2006 to 2009 and 2012 at El Castillo and Quiriquina stations using the different LU maps, respectively. Moreover, changes in components of the catchment hydrology were evaluated on a decadal basis for each LU map (1990-1999, 2000-2009, 2010-2019). The water balance was analyzed by examining the hydrological model components Surq, Latq, WY, Rchrg, and actual ET.

Efficiency Criteria for Model Evaluation
For model performance verification, the simulated and observed time series from El Castillo and Quiriquina stations were used. Hydrographs were compared first visually, and in a second step, quantitatively using statistical indices such as the coefficient of determination (R 2 ) and goodness of fit measures such as the Nash-Sutcliffe efficiency (NSE), relative Nash-Sutcliffe efficiency (rNSE), the standard deviation ratio of RMSE observations (RSR), and the percent of the model BIAS (PBIAS). R 2 is used to measure the consistency of the simulated model output against the observed data (Equation (2)). The values of R 2 vary between 0 and 1, less error variance is indicated by higher values [48,49]: where Q obs i and Q sim i represent the observed and simulated flow during each day. Whereas Q obs and Q sim are the observed and simulated means, respectively.
NSE is a standardized statistical method that determines the relative magnitude of the residual variance compared to the measured data variance (Equation (3)). NSE values vary from −∞ to 1, where 1 corresponds to a perfect match between the observed and simulated values, and NSE ≤ 0.5 indicate unsatisfactory model performance [49,50]: The rNSE is a modification of the NSE (Equation (4)), which calculates the efficiency between the simulated and observed data without gaps in the time series [51]: RSR is calculated as the ratio of Root Mean Square Error (RMSE) and standard deviation of the measured data (Equation (5)). An RSR of zero indicates the optimal value, while RSR > 0.7 represents unsatisfactory model performance [50].
PBIAS measures the estimation bias of the model (Equation (6)). PBIAS values can be positive or negative, indicating underestimation and overestimation, respectively, while the zero values represent the best simulation performance of the model [50].
In order to analyze the calibration and validation periods, recommended ranges of values for each efficiency criterion at daily time step were used [48][49][50].

Sensitivity Analysis
A simulation of 41 years period (1979 to 2019) was made using SWAT+ based on the land use map of 1997, considering in total 40 sub-basins, 255 channels, and 6119 HRUs. Even with field values, the K parameter turned out to be highly sensitive during the sensitivity analysis, possibly due to soil heterogeneity and the regionalization of soil parameters in the model. Therefore, to achieve a better performance of the model, seven parameters were calibrated using a multi-site calibration method considering a minimum and maximum range of sensitive values (Table 5).

SWAT+ Model Calibration and Validation
Satisfactory model performance, given poor data availability, was achieved for the multi-site calibration method with an NSE value of 0.52. Similarly, acceptable model performance was achieved for both the El Castillo and Quiriquina stations for calibration and validation periods ( Table 6).
R 2 values are considered good and satisfactory, respectively, albeit they change only marginally from 0.63 in the simulation to 0.64 for the calibration period at the El Castillo station and from 0.58 to 0.57 at the Quiriquina station. According to the classifications of Moriasi et al. [50], during the calibration period, the model showed a satisfactory performance with NSE of 0.53 at El Castillo station (Figure 4a) and 0.50 at the Quiriquina station. In addition, values of rNSE were considerably improved from 0.81 to 0.87 for El Castillo and from 0.75 to 0.83 for Quiriquina station, in both cases rNSE values can be classified as very good. In the case of RSR, with a value of 0.69 for El Castillo and 0.70 for Quiriquina, it can be considered also satisfactory. In the case of PBIAS, the resulting values have been considered unsatisfactory. However, the poor performance can be explained by bias in some unreliable peaks of daily discharge values.   According to the classifications of Moriasi et al. [50], for the validation period, the model showed a satisfactory performance with an NSE of 0.64 for both El Castillo and Quiriquina stations (Figure 4b). In the case of rNSE, values of 0.70 and 0.61 are considered good for the El Castillo station and satisfactory for the Quiriquina station. In the case of RSR, values of 0.58 and 0.59 were considered good for both El Castillo and Quiriquina stations, respectively. Similarly to the calibration period, PBIAS were considered unsatisfactory for both stations, mainly associated with the bias of unreliable discharge peaks.

Model Performance Under Land Use Change
For the evaluation of land use changes, each land use map was replaced in the original model without any further modification. This intervention modified the total number of HRUs created with LU 1997. Starting from 6119 originally, for LU 2009 and LU 2016, the number increased to 6638 and 6708 HRUs, respectively. El Castillo and Quiriquina stations were used to validate the performance of the model with both new land use maps, taking observations for each corresponding decade (Table 7)

Past Land Use Change Impacts on the Catchment Hydrology
The results based on each land use map were separated by each respective decade (1990-1999, 2000-2009, 2010-2019). The impacts of past land use changes on the catchment hydrology are presented on a monthly, seasonal, and yearly basis. Additionally, the changes between the periods 1990 and 2019 are presented at the basin level.
The changes on the monthly scale can be found in Figure 5 and Table 8. Differences in monthly precipitation (PRECIP) have been simulated during the period 1990-2019. In particular the lowest values were observed under LU 2016. The most determinant changes during the period 1990-1999 and 2010-2019 are related to a reduction of rainfall during the months between March to July, while an increase in precipitations has been estimated during the months of January and February. In accordance with precipitation changes, a reduction in surface flow (SURQ) with an accumulated decrease during the May to July period (116.2 mm) and an increase from November to March period (19.1 mm) was simulated. As a result, there was a reduction of 117.6 mm of annual accumulated SURQ values between LU 1997 and LU 2016. In the case of lateral flow (LATQ), some differences have been simulated for each month considering reductions and increases, in line with the changes in SURQ. Therefore, water yield (WY) decreased during the same months as SURQ and LATQ. In the case of groundwater recharge (RCHRG), and similarly with the behavior of other parameters, a reduction was achieved for the period 2010-2019 under LU 2016 configuration. For evapotranspiration (ET), a systematic increase has been simulated for the monthly study period between LU 1997 and LU 2016, and an accumulated ET value of 43 mm, with exception of the months of April and December with reduction of 5.6 mm.  For evaluation on a seasonal scale (Figure 6), the seasons were defined as Summer (January, February, March), Autumn (April, May, June), Winter (July, August, September), and Spring (October, November, December). In the case of annual changes (Figure 7, Table 9), variation in precipitation values has been simulated during the three decades of study, particularly a major reduction was estimated in the last decade. In the decade of 2000-2009, an increase in PRECIP has been obtained in comparison to LU 1997; specifically, an accumulated yearly average value of 2010.2 mm is the maximum value for all the studied periods. Similarly, maximum values of SURQ, LATQ, WY, and RCHRG were obtained during the period with LU 2009.  The minimum accumulated yearly mean value for precipitation was obtained during the period 2010-2019 with LU 2016. In addition, the coefficient of variation is shorter during the last period, this indicates that the variability in rainfall amount has been lower in the last decade. In addition, a notorious reduction in the hydrological fluxes SURQ, LATQ, WY, and RCHRG emerges from simulations performed during the last decade of this study. Conversely, the higher yearly mean value of the ET flux was obtained with LU 2016.
At the basin level (Figure 8), the differences between periods under LU 1997 and LU 2016 occurred synchronously with a drastic decrease in precipitation, especially in the mountainous part of the catchment. As to be expected, this leads to a reduction in the amount of SURQ in the middle part and at the outlet of the catchment. While a decrease in LATQ has been calculated in the middle part of the watershed in addition to the mountain part of the catchment. As consequence, reduction of water yield was simulated in different areas of the catchment. Similarly, RCHRG shows a tendency to decrease during the study period. Moreover, an increase in evapotranspiration has been calculated for some areas of the catchment, particularly in areas affected by land use changes related to shrubland reduction, increases in agricultural areas, and the expansion of pine tree plantations. The water extraction to supply the irrigation demand in the agricultural areas has increased during the last decades, contributing to a decrease in surface flow, water yield, and groundwater recharge, in addition to increases in evapotranspiration have been simulated in areas with crops of the catchment. Thus, for groundwater recharge and evapotranspiration, the strongest changes are to be expected mostly by remotion of shrublands and the establishment of pine trees and agricultural crops.

SWAT+ Model Results
SWAT+ is a new and restructured version of the original SWAT, a broadly accepted and used tool for the small watershed to river basin-scale modeling. In relation to the complexity of the input parameters, for our study case, the software showed some restrictions, specifically concerning the meteorological input. The CAMELS-CL data set [39] was selected as the most suitable regional-scale input data set. However, due to some gaps in the data records, especially discharge, flow rate values can be underestimated. Other authors have also described analogous situations of poor peak estimation using different hydro-meteorological data, for example the CHIRPS database, which they attribute to systematic overestimation of precipitation in the mountain areas of Chile [35,52]. Nevertheless, it was possible to obtain satisfactory results with CAMELS-CL in terms of calibration and validation with the LU 1997, and, therefore, with the validation of LU 2009 and LU 2016 land covers. This supported our decision to evaluate the impacts of land use change on the hydrology of the catchment using the CAMELS-CL dataset. However, some runoff peaks in the observed data for both El Castillo and Quiriquina stations were very difficult to model adequately, a fact that we attribute to an overestimation of discharge at these stations.

Model Efficiency Criteria
Model performance was evaluated under different efficiency criteria following recommended ranges of values at a daily time step. The model was successfully calibrated and validated under the land use map of 1997. In addition, the model output coherence was successfully validated under LU 2009 and LU 2016. Following the recommendations of Moriasi et al. [49], R 2 at daily time step can be considered satisfactory with values >0.5, at the contrary to monthly time step where the satisfactory qualification is achieved with values >0.7. Therefore, R 2 values were considered satisfactory and good for the evaluated periods even with unreliable discharge peaks and gaps in observed data. Similarly, NSE values were encouraging for all cases, even with the gaps in observed data. In the case of rNSE, an useful efficiency criteria under data scarcity conditions [51], even with the inclusion of some unreliable runoff peaks, rNSE values were classified as satisfactory, good, and very good. Unlike NSE, rNSE calculation includes observed means (Q obs ) and ignores the gaps in observed data, thus, leading to higher indicator values. In the case of PBIAS, unsatisfactory results were calculated for all cases, possibly due to the short-term period of available data, in addition to the unreliable peak runoff values and some gaps from observed discharge data. The underestimation of discharge peaks, on the other hand, could be explained by different factors, like lack of snow measurement stations in the catchment, measurement errors, or heterogeneity of the hydro-meteorological data, particularly due to the complex geomorphology of mountain areas. Additionally, unreliable peak events have been reported in observed discharges for other catchments in Chile [33,35]. However, we calculate the PBIAS parameter using values in daily time steps, in contrast to other studies that use seasonal, monthly, or yearly values [49,50,53]. Therefore, the classification can be overestimated and the impact of the time-step choice becomes more relevant, particularly in data-poor regions.

Impacts of the LUC on the Catchment Hydrology
In addition to the evident combined effect of precipitations and LUC on surface runoff during the study period, it is important to mention that the changes in the remaining hydrological fluxes can also be influenced by different factors. This was explained by Yin et al. [54], who emphasized that in addition to climate parameters, other controlling factors, such as the total number and the spatial arrangement of SWAT+ Hydrological Response Units (HRUs) as basic spatial modeling entities, may also impact certain basin-internal hydrological fluxes like surface runoff. For instance, the change of spatial resolution of the model by increasing the number of HRUs for each land use map, by simultaneously leaving other parameters unchanged, leads to a potential increase in surface runoff [55]. This indicates that such type of change in model structure can affect the simulated water yield and streamflow estimation by the model. Over our chosen study period, a maximum number of 589 HRUs was reached for the Longaví by QSWAT+ by changing the land use map from LU 1997 to LU 2016. As confirmed by other researchers [55,56], the model is clearly sensitive to the number and arrangement of HRUs for a particular catchment. However, contrary to the findings of Her et al. [55], runoff and other hydrological fluxes decreased in our application from 2010 to 2019, even with the new configuration of 589 HRUs. Therefore, to keep the number of HRUs steady and not affect the water fluxes, the use of decision tables implemented in SWAT+ for the development of land use scenarios seems to be an alternative approach to configure SWAT+ dynamically [40]. Based on our experience, additional studies on the development of procedures for the evaluation of land use scenarios with SWAT+ are recommended, especially in data-poor areas subject to continuous land use changes due to varying agricultural practices and forestry production.
Together with decreases in precipitation, our modeling results further suggest that the land use changes related to the expansion of mixed forest and pine tree areas at the expense of shrubland and agricultural areas contribute to an overall alteration of water redistribution in the catchment, by affecting groundwater recharge and evapotranspiration. Such redistribution shifts further in the future. Under past and current management practices related to deforestation, where the fragmentation of native forests and a decrease in patch density have been steadily increasing across Central Chile [57], the scenario is not encouraging. In addition, a decline in precipitation, associated with an artificial increase in evapotranspiration due to land use changes during the last 30 years, may have led to a net decrease in surface flow, lateral flow, and overall catchment water yield, as suggested by the SWAT+ simulations. A very similar trend has also been observed by Martínez-Retureta et al. [35], who observed annual increases in ET, and a decrease in water yield on Chilean coastal basins in conjunction with exotic tree plantations.
It is known that mountain ecosystems are characterized by particular natural regulatory functions [58], especially in areas with forested landscapes, by controlling the local climate [59], and water quantity and quality. The climate regulation and provision of clean water on the local scale can be affected through land use changes by disturbing the equilibrium of the local micro-climate by changing variables such as evapotranspiration, surface runoff, or groundwater recharge. In the Longaví catchment, a clear decrease in groundwater recharge has been reproduced by SWAT+ simulations during the last decade, in particular under the LU 2016 land use scenario. It is known that the impacts of aquifer exploitation can decrease water availability [60]. As a consequence, it is necessary to pay attention to unsustainable groundwater use to avoid problems related to (drinking) water scarcity in the future. Therefore, it is imperative to implement a comprehensive regulatory system that considers current climatic conditions and territorial needs. Additionally, policies and institutions for the protection of natural ecosystems need to be reinforced, with a clear emphasis on ecosystem services, such as the natural regulation of water flows and the protection of its quality, as well as environmental governance under conditions of an uncertain future climate.
Although unsustainable water use is not apparent, collateral long-term consequences are foreseeable and include progressive ground and surface water depletion. Moreover, the reduction of surface flow can affect the future transfer of water from the Longaví catchment to the Digua dam through inter-basin water transfer. Improvements in discharge measurements are necessary to avoid unreliable runoff peaks and gaps in the time series. In addition, there is a lack of snow measuring points in mountain areas; therefore, we highlight the importance of continuing to support the development and maintenance of hydro-meteorological and discharge stations in Chile. The actual quantification of precipitation decline due to climate change and the quantification of intrinsic uncertainties of the climate change projections for the region are a matter of ongoing investigation.

Conclusions
The paper presents the results of a study aimed at quantifying the combined effect of precipitation and the influence of three consecutive decades of land use changes in the hydrology of Longaví catchment, central Chile. Three land use maps of the basin were jointly used to track the progressive transition from predominantly native vegetation toward intensive agricultural and forestry exploitation. The catchment hydrology was analyzed by simulating hydrological fluxes with the agro-hydrological model SWAT+. The model was calibrated and validated using daily discharge records. Noticeable alterations of the catchment hydrology by combined effects of the decline in precipitation together with land use transitions have been changing the partitioning of surface flow, lateral flow, groundwater recharge, and evapotranspiration. As a result, the overall catchment water balance is altered after each new land use configuration. The principal changes in land use occurred through the expansion of areas devoted to urban sprawl, crop plantations, forestry, and pine tree production. These modifications combined with decreases in precipitation, caused negative trends in surface flow and hence a progressive decline in mean annual water yield between 2000 and 2019. The study confirms that land use transition has progressively affected internal water redistribution, with strongest impacts on groundwater recharge and evapotranspiration. In this context, it is important to consider the potential long-term effects of agriculture, forest, and pine tree production as important factors that have the potential to affect not only the hydrology in the Longaví catchment but also that of neighboring basins in the future due to inter-basin transfers. Further studies related to the potential effects of climate change on the region are a matter of ongoing research.

Conflicts of Interest:
The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Abbreviations
The following abbreviations are used in this manuscript: