Modelling the Effects of Historical and Future Land Cover Changes on the Hydrology of an Amazonian Basin

: Land cover changes (LCC) affect the water balance (WB), changing surface runoff (SurfQ), evapotranspiration (ET), groundwater (GW) regimes, and streamﬂow (Q). The Tapaj ó s Basin (southeastern Amazon) has experienced LCC over the last 40 years, with increasing LCC rates projected for the near future. Several studies have addressed the effects of climate changes on the region’s hydrology, but few have explored the effects of LCC on its hydrological regime. In this study, the Soil and Water Assessment Tool (SWAT) was applied to model the LCC effects on the hydrology of the Upper Crepori River Basin (medium Tapaj ó s Basin), using historical and projected LCC based on conservation policies (GOV_2050) and on the “Business as Usual” trend (BAU_2050). LCC that occurred from 1973 to 2012, increased Q by 2.5%, without noticeably altering the average annual WB. The future GOV_2050 and BAU_2050 scenarios increased SurfQ by 238.87% and 300.90%, and Q by 2.53% and 2.97%, respectively, and reduced GW by 4.00% and 5.21%, and ET by 2.07% and 2.43%, respectively. Results suggest that the increase in deforestation will intensify ﬂoods and low-ﬂow events, and that the conservation policies considered in the GOV_2050 scenario may still compromise the region’s hydrology at a comparable level to that of the BAU_2050.


Introduction
Over the last four decades, the Brazilian Amazonian ecosystem has been impacted by logging, pasture ranching, mining, expansion of the road network, and agricultural exploration [1][2][3]. Although the deforestation rates of large areas have declined over recent years, small-scale clearings have increased by 34% in a recent 14 year period, representing new challenges for forest conservation [4]. In the Tapajós Basin (southeastern Amazon), these small-scale clearings are associated with small-scale gold mining activities, which also compromise the water quality in the river [5,6]. Furthermore, the water resources of important Amazonian basins, such as Madeira's, Xingu's, and Tapajós River Basins, are also threatened by both built and planned dams, and by increasing mining activities [7,8]. Specifically, in the Tapajós River Basin, existing and planned dams may affect the The main economic activity in the Crepori Basin is small-scale gold mining, with several gold mining sites spread all over the basin, and even outside the established Gold Mining Reserve ( Figure  1a). On the other hand, logging, agriculture, and pasture are still incipient, resulting in low largescale deforestation rates [6]. However, sketchy and illegal infrastructures used by miners, such as small roads and airstrips, are important drivers for increasing human settlements [5,6,37] that potentially lead to an increase in future deforestation [1,38], mainly under a less effective environmental legislation. Moreover, the construction of a series of dams in the Tapajós River are planned in the region near the Crepori River Basin [8]. In addition, gold prices and new mining legislation can possibly stimulate the opening of new gold mining sites in the region [6], fueling the expectations of population growth and intense LCC in the near future [38,39].

Method and Data
The Soil and Water Assessment Tool (SWAT) [40] was applied in this study to model the impacts of LCC on the hydrologic regime of the UCRB. Within the model framework, land cover scenarios of LCC that occurred over the last 40 years were applied to represent past and recent LCC, whereas The main economic activity in the Crepori Basin is small-scale gold mining, with several gold mining sites spread all over the basin, and even outside the established Gold Mining Reserve (Figure 1a). On the other hand, logging, agriculture, and pasture are still incipient, resulting in low large-scale deforestation rates [6]. However, sketchy and illegal infrastructures used by miners, such as small roads and airstrips, are important drivers for increasing human settlements [5,6,37] that potentially lead to an increase in future deforestation [1,38], mainly under a less effective environmental legislation. Moreover, the construction of a series of dams in the Tapajós River are planned in the region near the Crepori River Basin [8]. In addition, gold prices and new mining legislation can possibly stimulate the opening of new gold mining sites in the region [6], fueling the expectations of population growth and intense LCC in the near future [38,39].

Method and Data
The Soil and Water Assessment Tool (SWAT) [40] was applied in this study to model the impacts of LCC on the hydrologic regime of the UCRB. Within the model framework, land cover scenarios of LCC that occurred over the last 40 years were applied to represent past and recent LCC, whereas future LCC scenarios were applied considering two different policy conditions: "Business as Usual", and "Governance" [31].

Data Used for Model Set Up
Hydrological modelling with SWAT requires spatial data that represent key basin characteristics affecting its hydrology, such as topography, land cover, soil types, and climate (Table 1). While rainfall is the primary climatic driver, observed streamflow data were also required for model calibration and validation.

Climate data
Daily precipitation estimates were provided by the TMPA 3B42 v.7 product from the Tropical Rainfall Measuring Mission (TRMM) [41], whereas the daily maximum and minimum air temperatures were derived from the ERA Interim Daily product [42]. The temperature dataset was cross-referenced with in situ measurements from the weather station located near the Basin (at Itaituba city) [32]. It was found that the ERA Interim estimates were overestimating the minimum temperatures from 1998 to 2000 by 1.58 • C, whereas the maximum temperatures were underestimated by 7.58 • C for the entire series. Therefore, a bias correction was applied to the data series, by summing −1.58 • C to the minimum temperatures from 1998 to 2000, and 7.58 • C to the maximum temperatures of the entire series. Similarly, the precipitation data from the TRMM were also cross-referenced to the rain gauge measurements at Itaituba station, but no systematic biases were observed. Moreover, rainfall estimates from the 3B42 product were validated for the whole of the Tapajós Basin by Collischonn et al. [46]. The authors state that TRMM estimates are reliable and may be used as rainfall input for hydrological modelling in the Tapajós region where the rain gauge network is sparse, which is the case in the UCRB. Since no station data were available on wind speed, relative humidity, and solar radiation in the region, the Hargreaves method [47] was selected to calculate the basin's potential evapotranspiration.

Elevation and topography
The elevation data was derived from the Digital Elevation Model (DEM) (Figure 1c) provided by the Shuttle Radar Topography Mission (SRTM), with 1 arc second (~30 m) of spatial resolution [43]. The watershed boundary, its slopes, and the drainage network were retrieved from the DEM, and used to define the Hydrological Response Units (HRU). HRUs were, in turn, defined as areas with homogeneous land cover, soil type, and slope, for which the model simulated the hydrological variables, such as streamflow, evapotranspiration, and groundwater [40].

Soil data and parameterization
The soil map made available by the Brazilian Institute of Geography and Statistics (IBGE) [44] ( Figure 1a) was combined with the soil texture and depth information retrieved from the databases of the RADAMBRASIL project [48], the Brazilian Agricultural Research Corporation (EMBRAPA) [49,50], and the Brazilian National Agricultural Research Department of the Brazilian Ministry of Agriculture (DNPEA/MA) [51]. Soil hydrological groups and maximum root depths were also derived from the literature [52,53]. Pedotransfer functions, specifically derived for tropical and Amazonian soils, were applied to estimate key soil hydraulic parameters, such as the Saturated Hydraulic Conductivity [54][55][56][57], Soil Bulk Density [58], and Available Water Capacity [54,57,59].

Definition of Land cover scenarios
The impacts of LCC on the UCRB streamflow were assessed by running the hydrologic simulations using four land cover scenarios that represent three periods (past, recent, and future land cover conditions). The first scenario was the land cover map corresponding to the year 1973 [6] ( Figure 2), which represents a quasi-pristine land cover condition, since LCC in the UCRB has become noticeable only at the beginning of the 1980's as result of the gold rush [60]. This land cover scenario was defined by classifying Landsat5/TM scenes featuring the UCRB at 1973. The small gold mining sites present in the basin at 1973, as well as those present in the other scenarios, were classified as barren land, since the SWAT model database does not have a specific land cover class for gold mining sites. The second scenario represents more recent LCC and is represented by the UCRB's land cover maps of 1998, 2003, and 2010 ( Figure 2), which correspond to the years with the highest LCC for the period from 1998 to 2012 [61]. Landsat5/TM scenes of these years were processed using supervised classification techniques to map the land cover classes as: Forest, pasture, barren land, and open water. These maps were updated during the simulation and, therefore, are hereafter referred as one single scenario: 1998-2003-2010. The third and fourth scenarios represent potential future LCC, projected for 2050, under different policy conditions: The "Business as Usual" (BAU_2050) and the "Governance" (GOV_2050) scenarios ( Figure 2). These scenarios were developed by Soares-Filho et al. [31], using the SimAmazonia model and are based on policy-sensitive simulations of future patterns of deforestation at the Amazon Basin, from 2002 to 2050 [31,62]. The BAU_2050 scenario was built by projecting the deforestation rates, estimated from remote sensing monitoring images from 1997 to 2002, and adding the effect of paving major roads in the region. In addition to this, the GOV_2050 scenario also assumed a 50% limit imposed for deforested land within each basin's sub region, as well as assuming that existing and proposed Protected Areas play a decisive role in limiting deforestation [31,62].

Observed streamflow
Observed streamflow data was required for model calibration and validation. The only stream gauge located at the UCRB ( Figure 1a) recorded observed streamflow at the basin's outlet, mainly from 2003 to 2012. The data is made available by the Brazilian National Water Agency (ANA) [35] and was divided into two datasets: From 2003 to 2009, for model calibration, and from 2010 to 2012, for model validation.

Model Set Up
SWAT divides the UCRB into sub-basins based on topographical information, which are, in turn, divided into HRU, based on their land cover and soil type. In this study, the UCRB was divided into 33 sub-basins with a total of 274 HRUs. Equation (1) is the central equation in SWAT that governs the water balance in the basin, as represented by the different hydrological variables [40].
where SW t is the final soil water content (mm H 2 O); SW 0 is the initial soil water content (mm H 2 O); t is the time (days); R day i is the precipitation on the day i

Sensitivity Analysis, Model Calibration and Validation, and Scenarios Application
The period between 1998 and 2002 was used for the model 'spin-up' stage, when the model is run until physically sensible initial conditions are set, especially for the soil moisture. The period between 2003 and 2009 was used for calibration, whereas the period between 2010 and 2012 was used for model validation. Land cover maps of 1998, 2003, and 2010 were updated during the model spinup, calibration, and validation. Thus, this dynamic land cover map is referred to as a single scenario: The automatic calibration of a large number of parameters included in conceptual and semidistributed models, such as SWAT, can be ineffectively time-consuming and computationally demanding. Therefore, the sensitivity analysis was performed to select the set of parameters that have a larger influence on the model's results, optimizing the calibration procedure. Also, according to Yang et al. [63], model calibration should consider knowledge of hydrological processes in the basin. Thus, key hydrological processes, such as evapotranspiration and groundwater flow, should be adjusted first, using a manual calibration procedure (trial and error) guided by expert knowledge and initial values derived from the literature [64][65][66][67][68][69]. This procedure was applied until the average rate of annual evapotranspiration and the average rate of annual groundwater flow were comparable to rates reported for similar regions near the Crepori Basin [30,[70][71][72][73]. Then, using the SWAT-CUP tool [74], the parameters' global sensitivity analysis was performed. In this approach, the parameter

Sensitivity Analysis, Model Calibration and Validation, and Scenarios Application
The period between 1998 and 2002 was used for the model 'spin-up' stage, when the model is run until physically sensible initial conditions are set, especially for the soil moisture. The period between 2003 and 2009 was used for calibration, whereas the period between 2010 and 2012 was used for model validation. Land cover maps of 1998, 2003, and 2010 were updated during the model spin-up, calibration, and validation. Thus, this dynamic land cover map is referred to as a single scenario: 1998-2003-2010.
The automatic calibration of a large number of parameters included in conceptual and semi-distributed models, such as SWAT, can be ineffectively time-consuming and computationally demanding. Therefore, the sensitivity analysis was performed to select the set of parameters that have a larger influence on the model's results, optimizing the calibration procedure. Also, according to Yang et al. [63], model calibration should consider knowledge of hydrological processes in the basin. Thus, key hydrological processes, such as evapotranspiration and groundwater flow, should be adjusted first, using a manual calibration procedure (trial and error) guided by expert knowledge and initial values derived from the literature [64][65][66][67][68][69]. This procedure was applied until the average rate of annual evapotranspiration and the average rate of annual groundwater flow were comparable to rates reported for similar regions near the Crepori Basin [30,[70][71][72][73]. Then, using the SWAT-CUP tool [74], the parameters' global sensitivity analysis was performed. In this approach, the parameter sensitivity was assessed by performing a t-test on a regression relating the objective function values (regarding the model's output) against the parameters values, which were sampled using the Latin Hypercube Sampling (LHS) [74]. Next, the most sensitive parameters were ranked, from the most sensitive to the least sensitive, and included in the automatic calibration, which was performed using the Parameter Solution (ParaSol) optimization algorithm [75]. The model was then validated against an observed streamflow dataset that was not used in the calibration procedure. Finally, the model performance was assessed for the entire period of the rainfall data used for the simulations (2003 to 2012), using the Coefficient of determination (R 2 ) (Equation (2)), the Nash-Sutcliffe Efficiency (NSE) [76] (Equation (3)), Percent Bias (PBIAS) [77] (Equation (4)), and the Root Mean Squared Error Observations Standard Deviation Ratio (RSR) [78] (Equation (5)), as recommended by Moriasi et al. [79].
where n is the number of observed data, Q obs,i is the observed streamflow on the day or month i, Q obs is the average observed streamflow, Q sim,i is the streamflow simulated on the day or month i, and Q sim is the average simulated streamflow.
To evaluate the impacts of LCC on the hydrological regime of the UCRB, the calibrated model was run separately for each land cover scenario (past: 1973; recent: 1998-2003-2010; and future: 2050 'BAU_2050', and 'GOV_2050'), using the same rainfall dataset (corresponding to the period between 1998 and 2012, with 1998 to 2002 used for model spin-up-i.e., with ten years' simulation period) for each run.

Sensitivity Analysis, Model Calibration, Validation, and Performance Assessment
As result of the sensitivity analysis, 14 parameters were defined as significantly sensitive (p-value < 0.05) ( Table 2), including several groundwater parameters. All 14 most sensitive parameters were included in the automatic calibration, which resulted in the calibrated range and the best calibrated values presented in Table 2. Since data and detailed information about the subsurface water in the Crepori region are scarce, parameter values related to groundwater (ALPHA_BF, ALPHA_BNK, GW_DELAY, GWQMN, and REVAPMN_FRSE) ( Table 2) were only defined in the calibration procedure, respecting the limits of physically meaningful values indicated by Neitsch et al. [40]. Sensitive soil parameters were calibrated for the first layers and according to the land cover type (Table 2), since land cover substantially affects the soil's physical-chemical properties, particularly in the first soil layers [30,80,81].
The calibrated values for soil bulk density (SOL_BD), available water capacity (SOL_AWC), and saturated hydraulic conductivity (SOL_K) agrees with values indicated by Tomasella and Hodnett [36] for Amazonian soils (0.7-1.2 Mg/m 3 for bulk density, around 0.7 mm·H 2 O/mmsoil for available water capacity, and values up to 1000 mm/h for saturated hydraulic conductivity). High saturated hydraulic conductivity values are particularly calibrated for soils covered with forest, which is associated to high root development that enhances the soil permeability [53,82]. The curve number parameter (CN2) is related to the runoff production, and their relatively low values calibrated for Red-Yellow Ferralsols under forested areas (Table 2) confirms the high infiltration expected in forested areas [82,83]. On the other hand, high values of the curve number adjusted for Plinthosols under forested areas (Table 2) were associated to the soil type, since the Plinthosols are characterized by the plinthic horizon, which reduces the infiltration and, consequently, increases the runoff [84]. The adjusted value of the effective hydraulic conductivity in the main channel alluvium (CH_K2) ( Table 2) indicates that the Crepori river bed material is mainly characterized by sand and gravel mixture with a low silt-clay content [85], which is corroborated by ICMBio [33] and is also in accordance with the classification of Tapajós Basin streams as clear-water rivers with low concentrations of suspended materials [86]. No parameter related to the deep aquifer was significantly sensitive, indicating that aquifers eventually present in the UCRB have low productivity and low interaction with surface hydrology. This is also corroborated by the fact that the basin is located at the Hydrogeological Province of the Brazilian Shield [87], which typically presents fractured aquifers that are usually deep, sealed, and of low production [88]. The best calibrated parameter values ( Table 2) led to a mean evapotranspiration rate of around 54% of the mean annual rainfall, which corresponds to those rates reported in the literature for other Amazon regions like the UCRB [70][71][72][73]. The groundwater parameters were adjusted in such a way that the streamflow simulated for the dry months equaled that observed at the river gauge.
The simulated streamflow matched the observed streamflow and was in synchrony with the variability in the rainfall data (Figure 3), for both calibration and validation periods ( Figure 3, Table 3). The model's performance for the calibration, validation, and the entire period of rainfall used for model simulations (2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012) can be classified as 'good' and 'very good', according to Moriasi et al. [79] ( Table 3). Red-Yellow Ferralsols under forested areas (Table 2) confirms the high infiltration expected in forested areas [82,83]. On the other hand, high values of the curve number adjusted for Plinthosols under forested areas (Table 2) were associated to the soil type, since the Plinthosols are characterized by the plinthic horizon, which reduces the infiltration and, consequently, increases the runoff [84]. The adjusted value of the effective hydraulic conductivity in the main channel alluvium (CH_K2) ( Table 2) indicates that the Crepori river bed material is mainly characterized by sand and gravel mixture with a low silt-clay content [85], which is corroborated by ICMBio [33] and is also in accordance with the classification of Tapajós Basin streams as clear-water rivers with low concentrations of suspended materials [86]. No parameter related to the deep aquifer was significantly sensitive, indicating that aquifers eventually present in the UCRB have low productivity and low interaction with surface hydrology. This is also corroborated by the fact that the basin is located at the Hydrogeological Province of the Brazilian Shield [87], which typically presents fractured aquifers that are usually deep, sealed, and of low production [88]. The best calibrated parameter values ( Table 2) led to a mean evapotranspiration rate of around 54% of the mean annual rainfall, which corresponds to those rates reported in the literature for other Amazon regions like the UCRB [70][71][72][73]. The groundwater parameters were adjusted in such a way that the streamflow simulated for the dry months equaled that observed at the river gauge. The simulated streamflow matched the observed streamflow and was in synchrony with the variability in the rainfall data (Figure 3), for both calibration and validation periods ( Figure 3, Table  3). The model's performance for the calibration, validation, and the entire period of rainfall used for model simulations (2003-2012) can be classified as 'good' and 'very good', according to Moriasi et al. [79] (Table 3).   The Q curve simulated for all scenarios displays a slight delay in comparison to the rainfall curve (Figure 4a,b), indicating that the UCRB has a lag-time response to the rainfall events of approximately 1 month. The starting and ending month of both high-flow (January to June) and low-flow (June to December) remained unchanged amongst the Q curves simulated for all LCC scenarios (Figure 4b). This indicates that, on average, none of the LCC scenarios used in this study was sufficient to noticeably change the river's seasonal streamflow dynamics of the UCRB. The percent differences between Q values, simulated using the 1998-2003-2010 scenario and that simulated using 1973 scenario (Figure 4c), show that the deforestation that occurred between 1973 and 2012 was not sufficient to noticeably increase the Q at an annual scale, and resulted in an increase in streamflow of about only 1% throughout the seasons. On the other hand, the change in Q, resulting from simulations with the GOV_2050 and BAU_2050 scenarios, shows higher values during the high-flow season and at the beginning of the low-flow season, with maximum increase rates of 4% and 5%, respectively (Figure 4c). During most of the low-flow season (from mid-August to mid-November), Q, simulated using the GOV_2050 and BAU_2050 scenarios, showed a reduction of up to 10% and 12%, respectively, for the month of October, whereas the Q, simulated using the 1998-2003-2010 scenario, showed the smallest increase rate (0.2%) for the same month (Figure 4c). The reduction of streamflow during the low-flow season is even larger when analyzing the percent differences between the monthly streamflow simulated for each month of the time series (from 2003 to 2012) ( Figure 5). The curve simulated for all scenarios displays a slight delay in comparison to the rainfall curve (Figure 4a,b), indicating that the UCRB has a lag-time response to the rainfall events of approximately 1 month. The starting and ending month of both high-flow (January to June) and lowflow (June to December) remained unchanged amongst the curves simulated for all LCC scenarios (Figure 4b). This indicates that, on average, none of the LCC scenarios used in this study was sufficient to noticeably change the river's seasonal streamflow dynamics of the UCRB. The percent differences between values, simulated using the 1998-2003-2010 scenario and that simulated using 1973 scenario (Figure 4c), show that the deforestation that occurred between 1973 and 2012 was not sufficient to noticeably increase the at an annual scale, and resulted in an increase in streamflow of about only 1% throughout the seasons. On the other hand, the change in , resulting from simulations with the GOV_2050 and BAU_2050 scenarios, shows higher values during the highflow season and at the beginning of the low-flow season, with maximum increase rates of 4% and 5%, respectively (Figure 4c). During most of the low-flow season (from mid-August to mid-November), , simulated using the GOV_2050 and BAU_2050 scenarios, showed a reduction of up to 10% and 12%, respectively, for the month of October, whereas the , simulated using the 1998-2003-2010 scenario, showed the smallest increase rate (0.2%) for the same month (Figure 4c). The reduction of streamflow during the low-flow season is even larger when analyzing the percent differences between the monthly streamflow simulated for each month of the time series (from 2003 to 2012) ( Figure 5).    Since the monthly mean streamflow was not averaged over the entire simulation period, the percent differences showed in Figure 5 display higher absolute values than those presented in Figure 4c. The 1998-2003-2010 scenario still shows relatively small streamflow increase and reduction rates (up to 2.5% and less than 0.5%, respectively), compared to the baseline (simulation using 1973 land cover scenario). Nevertheless, the overall percent differences between the streamflow simulated for the 1998-2003-2010 scenario and the 1973 scenario are positive, indicating that the small LCC between both scenarios was sufficient to depict a trend of an increase in streamflow, at a monthly time-scale analysis ( Figure 5). The streamflow simulated using the GOV_2050 and BAU_2050 scenarios showed increases as high as 11% and 22%, respectively, throughout the series (during the high-flow season), and reduction rates as high as 19% and 32%, respectively, during the low-flow season ( Figure 5). Regarding the other water balance components (WBC) as an annual average, the increased UCRB deforestation depicts a trend of increasing Surface Runoff (SurfQ) and Average Annual Streamflow (Q AA ), with decreasing Groundwater (GW) and Evapotranspiration (ET) (Figure 6). The percent differences of the WBC between the 1973 and 1998-2003-2010 scenarios are small (Figure 6), since the deforested areas in these scenarios are relatively similar (Figure 2) and, therefore, did not affect the basin's hydrology in a noticeable way. On the other hand, the WBC simulated for the GOV_2050 and BAU_2050 scenarios varies noticeably amongst all of them, especially for the SurfQ and GW (Figure 6), resulting in higher percent differences ( Figure 6).

Streamflow Simulations Corresponding to Land Cover Scenarios
Water 2018, 10, x FOR PEER REVIEW 11 of 18 Since the monthly mean streamflow was not averaged over the entire simulation period, the percent differences showed in Figure 5 display higher absolute values than those presented in Figure  4c. The 1998-2003-2010 scenario still shows relatively small streamflow increase and reduction rates (up to 2.5% and less than 0.5%, respectively), compared to the baseline (simulation using 1973 land cover scenario). Nevertheless, the overall percent differences between the streamflow simulated for the 1998-2003-2010 scenario and the 1973 scenario are positive, indicating that the small LCC between both scenarios was sufficient to depict a trend of an increase in streamflow, at a monthly time-scale analysis ( Figure 5). The streamflow simulated using the GOV_2050 and BAU_2050 scenarios showed increases as high as 11% and 22%, respectively, throughout the series (during the high-flow season), and reduction rates as high as 19% and 32%, respectively, during the low-flow season ( Figure 5). Regarding the other water balance components (WBC) as an annual average, the increased UCRB deforestation depicts a trend of increasing Surface Runoff (SurfQ) and Average Annual Streamflow (QAA), with decreasing Groundwater (GW) and Evapotranspiration (ET) ( Figure  6). The percent differences of the WBC between the 1973 and 1998-2003-2010 scenarios are small (Figure 6), since the deforested areas in these scenarios are relatively similar ( Figure 2) and, therefore, did not affect the basin's hydrology in a noticeable way. On the other hand, the WBC simulated for the GOV_2050 and BAU_2050 scenarios varies noticeably amongst all of them, especially for the SurfQ and GW (Figure 6), resulting in higher percent differences ( Figure 6).     -2003-2010-1973GOV_2050-1973BAU_2050-1973. BAU_2050-1973and GOV_2050-1973 are plotted at the secondary scale. Positive percent values indicate an increase of Q following deforestation since 1973, whereas negative percent values indicate a decrease of the streamflow due to deforestation since 1973.

Discussion
A hydrological modelling approach was implemented to study the possible impacts of contemporary and future LCC on the hydrology of the UCRB. The possibility of updating the land cover of the SWAT model during the model run helped in capturing the LCC dynamic and achieving good overall performance results. The manual adjustment of parameters, prior to the automatic calibration, based on previous knowledge of the physiographic characteristics of the basin, helped in defining the optimal parameter set and improve the reliability of the simulations. Although small, the deforestation that occurred in the basin between 1973 and 2012 was sufficient to depict a trend of increases in the streamflow at a monthly time-scale analysis ( Figure 5), corresponding to the increase in deforested areas, which is in accordance with the findings from previous studies in Amazon basins and other tropical regions [29,30,89,90]. This trend was more clearly displayed when simulating discharge during the high-flow season corresponding to the future BAU_2050 and GOV_2050 scenarios (Figure 4c). At the same time, the future scenarios also show a dramatic percent decrease in Q during the low-flow season, with simulations corresponding to the BAU_2050 scenario resulting in the highest percent change in Q, in relation to simulations corresponding to the 1973 quasi-pristine scenario (Figures 4c and 5). This shows that the potential deforestation by 2050, based on both the BAU_2050 and GOV_2050 scenarios, can remarkably change Q and other components of the water balance, especially by increasing SurfQ and Q, and decreasing GW and ET ( Figure 6). SurfQ was the only WBC that showed a substantial increase in the annual average values, following the increase of deforestation ( Figure 6), indicating that the SurfQ was the main component responsible for the slight increase in Q AA throughout the scenarios. Forest removal in tropical Basins typically reduces interception and ET, and increases SurfQ and Q during the rainy season [29,90,91]. However, during the dry season, Q is reduced when deforested areas increase (Figure 4c). Even though, during the dry season, deforestation has similar effects on interception, ET, and SurfQ, the rainfall rates are lower in this season and, therefore, the increase of SurfQ does not seem to noticeably contribute to Q. Besides, the GW is reduced as result of deforestation, since most of the water that enters the basin during the rainy season (as rainfall) is lost as SurfQ and Q, and does not contribute to GW storage. Furthermore, it is reported that the amount of water that infiltrates into the soil decreases from 130 mm to 60 mm throughout the dry season, rarely reaching more than 1m of Amazonian soil depth [92], which contributes to the critical reduction of GW in that season. Finally, since Q during the dry season is more dependent on the GW, the reduction of the GW results in the reduction of Q during the dry season.
The findings of this study support the "Sponge Effect Hypothesis" [80,91], which identifies tropical forests as being responsible for reducing peak discharge during rainy seasons (due to low SurfQ and high ET), while acting as a 'sponge' during the dry season, allowing higher soil infiltration, increasing GW, and, therefore, resulting in higher Q than those from deforested areas [80,91]. However, while it is generally understood that deforestation in tropical regions usually leads to an increase in high flows, the impacts of deforestation on low flows is known to be more difficult to predict [93].
The estimated past and potential future impacts of LCC on the Q and water balance indicate that deforestation increases in the UCRB may result in higher streamflow and lower low-flows during high-flow and low-flow seasons, respectively, even when considering successful implementation of environmental regulations and policies, as in the GOV_2050 scenario. In addition, changes in the water balance may also impact the water quality of the Crepori River, especially during the low-flow season, which is when small-scale gold mining activity intensifies, resulting in high sediment loads to the river [5]. The reduction of Q, simulated for the low-flow season in the GOV_2050 and BAU_2050 scenarios, may increase sediment concentration, since there will be less discharge in the river to dilute and wash the sediments and other pollutants. For the various dam projects that are under consideration in the Tapajós River, downstream of the Crepori River [8], these impacts on the streamflow and, potentially, on the water quality can increase maintenance costs of hydropower plants and periods with lower water volume storage, potentially resulting in lower power generation rates [27,94]. On the other hand, the predicted increase in streamflow during the high-flow season represents an increase in flood risks in cities and indigenous lands along the river. It is also worth stating that the forest removal may reduce rainfall rates, potentially reducing Q and also affecting the climate in distant areas [19][20][21][22]. This climate feedback can be even more critical, since the reduction of rainfall rates hampers forest regeneration.
Since the focus of this study was to analyze the impacts of LCC alone on the streamflow, climate implications of deforestation were not taken into account. However, more accurate predictions of LCC impacts on the hydrology of the UCRB must consider the climate feedback, as well as finer time scale analysis, such as a daily scale. These should be the focus of other future studies.

Conclusions
This study investigated the impacts of historical and potential future deforestation on the hydrology of the Upper Crepori River Basin (UCRB) using a hydrological modelling approach. The main findings of the study are: Deforestation that occurred in the UCRB since human exploration in the early 1970 s until recent years (1998-2012) was small, and was not sufficient to change the streamflow more than 2.5%, in an annual basis, and not sufficient to noticeably alter the annual water balance of the basin.
Hydrologic model simulation, with land cover scenarios developed by Soares-Filho et al. [30] for 2050, in both conditions of 'Business as Usual' and 'Governance' led to an increase in streamflow during the high-flow season and a decrease in streamflow in the low-flow season, intensifying river high-flows and low-flows events.
The projected deforestation for 2050 is also found to alter the water balance of the basin, increasing surface runoff and streamflow, and reducing groundwater and evapotranspiration.
Since both future land cover scenarios, GOV_2050 and BAU_2050, led to very similar results of the hydrological variables, this study shows that, even if conservation policies are respected, the expected future deforestation can still compromise the water resources in the UCRB, directly affecting local human communities, the biota, and possible planned hydropower generation.
In general, this study presents important analysis and information for water resources management and land use planning in a portion of the Tapajós Basin, regarding future water availability with respect to land cover change. Further studies are still needed to better understand the impacts of deforestation on the hydrology, including possible climate feedbacks and finer time-scale analysis.