Assessing Impacts of Climate Variability and Reforestation Activities on Water Resources in the Headwaters of the Segura River Basin (SE Spain)

Climate change and the land-use and land-cover changes (LULC) resulting from anthropic activity are important factors in the degradation of an ecosystem and in the availability of a basin’s water resources. To know how these activities affect the quantity of the water resources of basins, such as the Segura River Basin, is of vital importance. In this work, the Soil and Water Assessment Tool (SWAT) was used for the study of the abovementioned impacts. The model was validated by obtaining a Nash–Sutcliffe efficiency (NSE) of 0.88 and a percent bias (PBIAS) of 17.23%, indicating that SWAT accurately replicated monthly streamflow. Next, land-use maps for the years of 1956 and 2007 were used to establish a series of scenarios that allowed us to evaluate the effects of these activities on both joint and individual water resources. A reforestation plan applied in the basin during the 1970s caused that the forest area had almost doubled, whereas the agricultural areas and shrubland had been reduced by one-third. These modifications, together with the effect of climate change, have led to a decrease of 26.3% in the quantity of generated water resources, not only due to climate change but also due to the increase in forest area.


Introduction
Water resources have suffered an important decrease in recent decades due to the effects of human activities and climate change [1].This is especially the case in semi-arid regions, where water resources have become a critical element of socioeconomic development [2].Many recent studies have studied the driving factors of the hydrological cycle, reaching the conclusion that climate and land-use/land-cover changes (LULC) are two critical factors that affect the hydrological processes in the basin [3].On the one hand, climate change projections during this century in the Mediterranean region reveal a rise in temperature accompanied by a decrease in precipitation, resulting in a decrease of the water resources by over 20% [4].In fact, recent studies in the headwaters of the Segura River Basin (HWSRB) [5] predicted a decrease of over 50% by the end of the century.Furthermore, human activities such as agricultural irrigation expansion or urbanization affect the availability of water resources through their influence on runoff.This is especially true in Mediterranean Europe, where human activities have modified the landscape [6][7][8].Developing a better understanding of the influence of LULC on the hydrological processes of a watershed is essential for sustainable water resource management.This is particularly true in the Segura River Basin, which is located in a semi-arid region of southeast Spain, and where we can find the most water-stressed basin in Western Europe [9].
The HWSRB are relevant due to its 9% water resource contribution, in spite of the fact that it covers only 1.2% of the area over the total watershed [5].
The use of a physically based distributed hydrological model is a common approach to evaluating the impact of LULC changes on water resources [10].From the hydrological models available, we chose the Soil and Water Assessment Tool (SWAT) model due to the fact that this model has been used worldwide to simulate the hydrological cycle and has been used extensively to evaluate the impacts of LULC changes on watershed hydrology [11].Likewise, a multitude of studies have addressed separating the influences of human activities and climate variability on water resources using the SWAT model [3,[11][12][13].However, only two references have been found on the application of the SWAT model to assess the impacts of LULC changes in Spain.Molina-Navarro et al. [14] analyzed the potential effects of climate change and land-use management scenarios on water discharge and the water quality of the Pareja Reservoir located in the upper Tagus River Basin.Salmoral et al. [15] assessed the water-related impacts of land-use changes and agricultural practices, combining the use of the SWAT model and a water footprint assessment in the Genil River Basin (southern Spain).
Significant LULC have occurred in the southeastern part of Spain since the 1970s as a result of the progressive abandonment of dryland agricultural activities and the implementation of reforestation plans increasing forest cover [16].The strategy for reforesting the HWSRB was to introduce pine species; pinus nigra was used in upper parts of the basin, whereas pinus pinaster and pinus halepensis were used in the intermediate and lower parts of the basin [17].This increase of vegetation cover leads not only to a decrease in sediment yield [18] but also to a decrease in runoff generation.Several studies have analysed the potential water yield reduction due to forestation activities.Bosch and Hewlett [19] studied the effect of vegetation changes on the water yield on 94 different basins.They concluded that eucalypt and pine forest types cause a reduction of 40 mm in water yield per 10% change in cover, but in the case of shrubland, this reduction is four times lower (10 mm).Other authors added new basins in addition to those reviewed by Bosch and Hewlett; for example, Sahin and Hall [20] analysed 145 basins using a fuzzy linear regression analysis reaching similar conclusions, whereas Brown et al. [21] reviewed 166 basins focusing not only on water yield but also on low flows.Sun et al. [22] analysed the potential magnitude of annual water yield to forestation across China, suggesting that this reduction can vary from 30% in tropical regions to 50% in semiarid regions.Llorens et al. [23] studied a Mediterranean mountain basin showing a reduction in runoff up to 18% due to the reforestation activities using pine trees.
Due to the high water stress suffered in the Segura River Basin, the impact of forest restoration on water yield must be evaluated.Therefore, the objectives of this study were as follows: (1) to determine trends in the annual and monthly precipitation and temperature; (2) to calibrate and validate the SWAT model in a data-scarce river basin with a Mediterranean climate, such as the HWSRB; and (3) to analyze the impact of climate variability and reforestation activities on runoff and evapotranspiration (ET).

Study Area Description
The Segura River Basin represents one of the most arid zones of the Mediterranean area, presenting great heterogeneity in its flow regimes [24].It is located in the southeastern part of Spain (Figure 1).The SWAT model was applied to the 235 km 2 of the HWSRB that flow into the Anchuricas Reservoir, which has a capacity of 6 hm 3 and exists to generate electricity.As the terrain is mountainous, it presents steep slopes with an elevation range from 898 to 1912 m.According to data from 1951 to 2015, the average annual precipitation is above 800 mm, ranging from 381 mm to 1447 mm.As a typical Mediterranean climate, precipitation is received mainly during spring and autumn while it is insignificant during summer.The average annual temperature was 11.8 • C.
The mean monthly flow of this river leaving the mountains varies between 0.74 m 3 /s in August to 4.82 m 3 /s in February.Land cover in the watershed is mostly a forest-dominated area and features Mediterranean shrubland vegetation, which covers about 81% of the basin.The rest of the land is mainly used for range purposes.The main soil type is Rendzic Leptosol, which presents a variable depth but always less than 50 cm, abundant stoniness, high carbon levels and good drainage [5].The study area is characterized by the fundamental role of the groundwater in the surface hydrology due to the fact that a significant part of the streamflow comes from groundwater sources.The large volume of precipitations and the abundance of outcropping carbonate rocks have caused the formation of karstic systems [25].
Sustainability 2018, 10, x FOR PEER REVIEW 3 of 13 always less than 50 cm, abundant stoniness, high carbon levels and good drainage [5].The study area is characterized by the fundamental role of the groundwater in the surface hydrology due to the fact that a significant part of the streamflow comes from groundwater sources.The large volume of precipitations and the abundance of outcropping carbonate rocks have caused the formation of karstic systems [25].

The Precipitation and Temperature Trend Analysis
In this study, Microsoft Excel template application MAKESENS (Version 1.0) was used to calculate the trends and changing rates of temperature and precipitation.This application includes the non-parametric Mann-Kendall (MK) test for the trend and Sen's non-parametric method for the magnitude of the trend.More details regarding the MAKESENS model can be found in Salmi et al. [26].The MK test was used to detect the change in the annual and monthly temperature and precipitation.The advantage of the MK test, which is one of the most widely used nonparametric tests for detecting trends in hydroclimatic series [27], is that it does not require the data to be distributed normally and it has a low sensitivity to abrupt peaks due to inhomogeneous time series [28].The complete description and formulae of the MK test can be obtained from Partal and Kahya [29].The magnitudes of the estimated changes in the trend of temperature and precipitation were estimated by using a simple non-parametric procedure that Sen developed.The Sen [30] method alleviates the consequences of anomalous trends by using the median of the series of slopes as the judgmental foundation [13].The is given by the following: where and are data values at times j and k, respectively.Estimator is the median over all combinations of record pairs for the entire dataset.A positive slope magnitude indicates an increasing tendency, whereas a negative slope magnitude indicates a decreasing tendency in the time series.

Description of SWAT Model
Hydrological modeling was performed using the SWAT extension for the QGIS interface, called QSWAT [31].SWAT is a continuous-time, long-term, semi-distributed and physically based model.SWAT considers the heterogeneity of a watershed by dividing it into multiple sub-watersheds based

The Precipitation and Temperature Trend Analysis
In this study, Microsoft Excel template application MAKESENS (Version 1.0) was used to calculate the trends and changing rates of temperature and precipitation.This application includes the non-parametric Mann-Kendall (MK) test for the trend and Sen's non-parametric method for the magnitude of the trend.More details regarding the MAKESENS model can be found in Salmi et al. [26].The MK test was used to detect the change in the annual and monthly temperature and precipitation.The advantage of the MK test, which is one of the most widely used nonparametric tests for detecting trends in hydroclimatic series [27], is that it does not require the data to be distributed normally and it has a low sensitivity to abrupt peaks due to inhomogeneous time series [28].The complete description and formulae of the MK test can be obtained from Partal and Kahya [29].The magnitudes of the estimated changes in the trend of temperature and precipitation were estimated by using a simple non-parametric procedure that Sen developed.The Sen [30] method alleviates the consequences of anomalous trends by using the median of the series of slopes as the judgmental foundation [13].The Q i is given by the following: where x j and x k are data values at times j and k, respectively.Estimator Q i is the median over all combinations of record pairs for the entire dataset.A positive slope magnitude indicates an increasing tendency, whereas a negative slope magnitude indicates a decreasing tendency in the time series.

Description of SWAT Model
Hydrological modeling was performed using the SWAT extension for the QGIS interface, called QSWAT [31].SWAT is a continuous-time, long-term, semi-distributed and physically based model.SWAT considers the heterogeneity of a watershed by dividing it into multiple sub-watersheds based on the river network and topography; subsequently, sub-watersheds are divided into hydrologic response units (HRUs) consisting of homogeneous land use, land management, and soil characteristics.These HRUs represent percentages of the sub-watershed area and are not identified spatially within a SWAT simulation [32].SWAT was designed to predict the impact of land management practices on the hydrology of agricultural watersheds as well as on sediment and contaminant transport in these watersheds [33].This model can simulate surface and subsurface flows, sediment generation and deposition, and nutrient fate and movement through rivers and landscapes.The hydrologic routines within SWAT account for unsaturated zone processes (i.e., infiltration, evaporation, plant uptake, lateral flow, and percolation), groundwater flows and snow accumulation and melt [34].SWAT simulates the hydrologic cycle based on water balance, which is controlled by climate, moisture, and energy inputs, such as daily precipitation, the maximum/minimum air temperatures, solar radiation, the wind speed, and the relative humidity [35].The water balance equation that governs the hydrological components of the SWAT model [36] is as follows: where SW t is the final water-soil content (mm), SW 0 is the initial water-soil content (mm), t is the time in days, R day is the amount of precipitation on day i (mm), and Q surf is the amount of surface runoff on day i (mm); ET i is the actual evapotranspiration on day i (mm); W seep is the percolation on day i; and Q gw is the amount of baseflow on day i (mm).The SWAT model offers three options to calculate potential evapotranspiration: Priestley-Taylor, Hargreaves, and Penman-Monteith.In this study, potential evapotranspiration was simulated using the Hargreaves method due to the fact that it requires only maximum and minimum temperatures.Actual evapotranspiration is calculated as the sum of actual plant transpiration, actual soil evaporation, and interception.Firstly, SWAT evaporates any precipitation intercepted by the plant canopy and then calculates the transpiration, sublimation, and soil evaporation using an approach similar to that of Ritchie [37].

Input Data for Hydrological Modelling
To study the effect of climate variability and reforestation activities, land-use data in 1956 and 2007 at a scale of 1:25,000 were provided by the Andalusian Network of Environmental Information (REDIAM) belonging to the Regional Government of Andalusia [38].In addition to land cover, a 25 m resolution digital elevation model from the National Geographic Institute [39] was used.The soil data were obtained from the Harmonized World Soil Map [40] with a spatial resolution of about 1 km.The SWAT model was driven by meteorological data, including precipitation and temperature data, from various sources.The precipitation was obtained from the Spanish National Meteorological Agency (AEMET) grid, version 1.0, which provides daily rainfall in Spain for the period of 1951 to 2016 with a spatial resolution of 5 km.More information about this dataset can be found in Peral-García et al. [41].Temperature data were collected from the fifth version of the high-resolution (approximately 10 km) gridded dataset called SPAIN02, where data are available from 1951 to 2016.Detailed documentation of the development and analysis of the SPAIN02 dataset can be found in Herrera et al. [42].The discharge data at the catchment outlet were available on the Hydrographical Study Centre website [43].

Model Setup
The subdivision into hydrologic response units (HRUs, unique computational units of land coverage and soil types with homogeneous hydrologic responses) was performed with the land use, soil, and slope maps.Threshold levels of 10% were established to eliminate minor land uses, soils and slopes in each sub-basin and to facilitate model processing, and the SWAT interface identified 44 HRUs.These threshold values are used by the interface to eliminate minor land uses, soil types or slopes in each sub-basin.
To assess the impacts of reforestation activities, the parameters we calibrated were based on land-use data in 1956 that belonged to the natural period of 1951-1970.Discharge data from 1 January 1964 to 31 December 1970 and from 1 January 1954 to 31 December 1963 were selected to run the model validation and calibration respectively.Three first years (1951)(1952)(1953) were used to allow the model parameters to reach equilibrium as a warm-up period.The sensitivity analysis and automatic calibration was done in this study using the sequential uncertainty fitting program (SUFI-2) in the SWAT Calibration and Uncertainty Programs [44].First, a global sensitivity analysis was performed to identify the most influential parameters in governing the streamflow.A ranking of parameter sensitivities was obtained after 500 model runs to see the impact of each parameter on the objective function [45].The automatic calibration procedure was used to determine the best parameter values based on the observed discharge, and using the Nash-Sutcliffe coefficient (NSE) as the objective function.A total of 1000 simulations were run in two steps of 500 simulations, with the parameters adjusted after the first 500 simulations.

Performance Evaluation Criteria
Model performance was assessed quantitatively using the NSE, the root mean square error (RMSE), the observations standard deviation ratio (RSR), and the percent bias (PBIAS), and qualitatively using graphical time series plots.
The RSR standardizes the RMSE using the RSR (STDEV obs ).The RSR is calculated as follows: where n is the total number of observations, Q obs i and Q sim i are the observed and simulated discharges at ith observation, respectively, and Q is the mean of the observed data over the simulation period.The RSR varies from the optimal value of 0 to a large positive value.The lower the RSR, the lower the RMSE, and the better the model simulation performance.
The NSE [46] is a normalized statistic that determines the relative magnitude of the residual variance compared with the measured data variance.The NSE is dimensionless, and its values range from negative infinity to 1, with an optimal value of 1 [47].It indicates a 1:1 line fit between observed and simulated data and is computed as follows: The PBIAS [48] measures the average tendency of the simulated data to be larger or smaller than their observed counterparts [47] and is calculated as: For model evaluation, we used the criteria that Moriasi et al. [49] proposed about the performance rating for the described statistics for a monthly time step.According to these criteria, the calibration/validation performance can be considered satisfactory when the NSE is greater than 0.5, the RSR is less than 0.5, and the PBIAS range is ±25%.Based on these scenarios, we analyzed the impacts on water resources due to climate variability (the differences between B and A), reforestation (the difference between C and B), and all factors (the difference between C and A) respectively.This approach of "one factor at a time" has been successfully applied in many studies to evaluate the effects of land-use change and climate variability on hydrology [11,12,50].Based on these scenarios, we analyzed the impacts on water resources due to climate variability (the differences between B and A), reforestation (the difference between C and B), and all factors (the difference between C and A) respectively.This approach of "one factor at a time" has been successfully applied in many studies to evaluate the effects of land-use change and climate variability on hydrology [11,12,50].

Climate Variability in the HWSRB
Climate trend analysis requires a minimum of 15-30 years for accurate results [51].In this study, the annual and monthly precipitation, temperature, and streamflow trended for the period 1951-2015 were analyzed using the MK test.The values of the MK test statistic (Test Z) and Sen's slope for various months are provided in Table 1.In terms of precipitation, it can be seen that the greatest slope was in the months of December, March, and June.According to Paredes et al. [52], the strong decrease of rainfall in March is mostly associated with a positive trend of the North Atlantic oscillation (NAO) index.It is observed that eight of the 12 months revealed a decreasing trend.On a seasonal timescale, the most important decreases take place in the summer and winter.Del Río et al. [53] previously suggested these trends.The annual precipitation decreased across the study area with a slope of 2.22 mm/year.Streamflow and precipitation show a similar decreasing trend in the study area.Influenced by the interaction between reforestation and climate variability factors, a large decrease in streamflow is occurring between October and April.
Annual maximum air temperature indicated a significant upward trend at the 5% significant level.

Climate Variability in the HWSRB
Climate trend analysis requires a minimum of 15-30 years for accurate results [51].In this study, the annual and monthly precipitation, temperature, and streamflow trended for the period 1951-2015 were analyzed using the MK test.The values of the MK test statistic (Test Z) and Sen's slope for various months are provided in Table 1.In terms of precipitation, it can be seen that the greatest slope was in the months of December, March, and June.According to Paredes et al. [52], the strong decrease of rainfall in March is mostly associated with a positive trend of the North Atlantic oscillation (NAO) index.It is observed that eight of the 12 months revealed a decreasing trend.On a seasonal timescale, the most important decreases take place in the summer and winter.Del Río et al. [53] previously suggested these trends.The annual precipitation decreased across the study area with a slope of 2.22 mm/year.Streamflow and precipitation show a similar decreasing trend in the study area.Influenced by the interaction between reforestation and climate variability factors, a large decrease in streamflow is occurring between October and April.Annual maximum air temperature indicated a significant upward trend at the 5% significant level.In addition, the highest magnitude of the trend was found in June.This result is in accordance with those that Del Río et al. [54] proposed.With the exception of winter months, the maximum air temperature showed a rising trend during the remaining months.The decreasing magnitude of the minimum air temperature was smaller compared with the increasing magnitude of the maximum air temperature.Compared with their monthly changes, the decreasing magnitudes in the cold months were larger than they were in warm months.The decrease of the minimum air temperature was especially significant in March (0.01 level of significance).

Land-Use Change in the HWSRB
Table 2 shows the land uses in the HWSRB in 1956 and 2007.Two main trends of LULC existed from 1956 to 2007: the decrease of shrubland and agricultural land and the increase of forests.Compared with 1956, agricultural land, shrubland, and woodland-shrub transitions decreased in areal extent by respective amounts of 11.56, 14.79, and 15.60 km 2 .In contrast, during the period of 1956−2007, the forest area showed an increase of 34.26 km 2 .Urban areas, water bodies, grassland, and barren land changed little compared with their baselines of 1956.

Validation of the SWAT model
The global sensitivity analysis found that the most influential parameters were the effective hydraulic conductivity in tributary channel alluvium (CH_K1), the U.S. Soil Conservation Service curve number (CN2), the moist bulk density (SOL_BD), the saturated hydraulic conductivity (SOL_K), the groundwater delay (GW_DELAY), the maximum canopy storage (CANMX), the lateral flow travel time (LAT_TTIME), and the available water capacity of the soil layer (SOL_AWC).Similar parameters were chosen for calibration in basins with similar characteristics [55].Even the sensitivity analysis did not give the soil evaporation compensation factor (ESCO) as one of the most sensitive parameters, but it was considered to be a significant parameter because it is highlighted in the literature and in previous work made for this basin [5,56].Sensitivity analysis was followed by automatic calibration.The optimal values and range of parameters are presented in Table 3.According to Table 4, the SWAT model performed well in both the calibration and the validation period, accurately simulating the discharge after the sensitive parameters were optimized.The R, NSE, PBIAS, and RSR of the calibration period (1954)(1955)(1956)(1957)(1958)(1959)(1960)(1961)(1962)(1963) for monthly runoff were 0.87, 0.86, −14.11%, and 0.38 respectively.Overall, the model had good performance not only during the calibration period but also during the validation period.This is not usual because parameters are optimized during calibration, and better performance is expected during this period [57,58].The results suggest that the SWAT model performed well and can be used for further analysis in the HWSRB.To illustrate the model calibration and validation, observed and simulated monthly discharge data in the Anchuricas Reservoir for the calibration (1954)(1955)(1956)(1957)(1958)(1959)(1960)(1961)(1962)(1963) and validation (1964)(1965)(1966)(1967)(1968)(1969)(1970) periods are plotted in Figure 3.
As can be seen in Figure 4, we have also validated the model by comparing the simulated results using the 2007 land-use map (scenario C) with observed streamflows during the period 1996-2015.The values of model evaluation statistics such as NSE, R, PBIAS, and RSR were 0.72, 0.85, 21.93, and 0.53, respectively.These values indicate that the calibrated model satisfactorily simulates the monthly runoff.

Impacts of Land-Use Change and Climate Variability on Runoff and ET at HWSRB
Table 5 shows the annual mean runoff and ET simulated by SWAT under the three different scenarios analyzed in this study.Following the recommendation of Li et al. [59], the simulated data rather than observed data were used to compare the hydrological effects under the three hypothetical scenarios presented in Section 2.4.Compared with the baseline scenario (Scenario A), simulated runoff in Scenario C decreased by 95.8 mm.This reduction is due to the combined effects of climate variability and the reforestation plan.The contrast between Scenarios A and B showed the influence of the climate variability.Surface runoff decreased significantly by 61.7 mm due to climate variation, which accounted for 64.4% of the total runoff reduction.The results revealed that the reforestation plan and climate variability both decreased the runoff.It should be noted that the contribution of climate variability (61.7 mm) is almost two times greater than that of the reforestation plan (34.1 mm).In terms of ET, the combined effect of land-use change and climate variation increased ET by 31.8 mm.The climate variation caused a decrease of 13.2 mm in the ET mainly due to the significant decrease in the precipitation shown in the trend analysis, and in spite of the slight temperature increase.It should also be noted that the reforestation plan contributed to an increase in the ET values by about 45 mm, exceeding by far the impact of climate variability on ET.These results are consistent with those of previous studies [22,23].

Conclusions
Considering the HWSRB, as a case study, the aim of this study was to distinguish the influences of human activities and climate variability on water resources by using the SWAT model.From the results of this study, the following conclusions can be stated: The main trend in land cover change was the conversion of agricultural areas and shrubland to forest, both of which are a result of the reforestation plan carried out during the 1970s.

3.
In general, climate change and reforestation activities tend to decrease the streamflow of the HWSRB.The contributions of climate change (64.4%) were larger than that of the reforestation plan (35.6%).4.
The results also revealed that the reforestation plan had a higher impact than climate variability did on ET in the HWSRB.

5.
Although reforestation plans can result in decreased soil erosion rates, runoff is also reduced.
The results obtained from this study may have strong implications for a basin that is already suffering from high stress.The impacts of the reforestation plan should be incorporated into water resource management plans to develop sustainable strategies.6.
Future reforestation plans in this area should strengthen the native shrub type of vegetation instead of increasing the forested area in order to preserve the water yield.

Figure 1 .
Figure 1.Location of the headwaters of the Segura River Basin (HWSRB) and its digital elevation model (DEM).

Figure 1 .
Figure 1.Location of the headwaters of the Segura River Basin (HWSRB) and its digital elevation model (DEM).
is the deviation of data being evaluated, expressed as a percentage.The optimal value of the PBIAS is 0.0, with low-magnitude values indicating an accurate model simulation.Positive values indicate model underestimation bias, and negative values indicate model overestimation bias.

2. 4 .
Framework for Separating Effects of Climate Change and LULC As can be seen in Figure 2, to study climatic and land uses' impacts on water resources in the HWSRB, three simulation experiments were set up as follows: scenario A used land use and climate conditions around the 1950s and 1960s (land use in 1956 and weather/climate for 1951-1970); scenario B fixed land use in 1956 and actual weather conditions for 1996-2015; and scenario C used land use in 2007 and actual climate conditions for 1996-2015.

Sustainability 2018 , 13 2. 4 .
10, x FOR PEER REVIEW 6 of Framework for Separating Effects of Climate Change and LULC As can be seen in Figure 2, to study climatic and land uses' impacts on water resources in the HWSRB, three simulation experiments were set up as follows: scenario A used land use and climate conditions around the 1950s and 1960s (land use in 1956 and weather/climate for 1951-1970); scenario B fixed land use in 1956 and actual weather conditions for 1996-2015; and scenario C used land use in 2007 and actual climate conditions for 1996-2015.

Figure 2 .
Figure 2. Scenarios setup and research framework.

Figure 2 .
Figure 2. Scenarios setup and research framework.

Figure 4 .
Figure 4. Observed and simulated monthly discharge data for the period 1996−2015 (scenario C) in Anchuricas Reservoir.

Table 1 .
Trend analysis results.Test Z is the Mann-Kendall (MK) test statistic; Q i is the Sen's slope estimator; a indicates a significance level of 0.1; b indicates a significance level of 0.05; c indicates a significance level of 0.01; d indicates a significance level of 0.001.

Table 2 .
The conversion percentage (%) of each land-use type from 1956 to 2007.

Table 4 .
Performance of Soil and Water Assessment Tool (SWAT) model during calibration and validation periods in the HWSRB.

Table 5 .
Simulated average annual runoff and ET under various scenarios (mm).