Effects of Land Use Changes on Streamflow and Sediment Yield in Atibaia River Basin—SP, Brazil

The Soil and Water Assessment Tool (SWAT) is often used to evaluate the impacts of different land use scenarios on streamflow and sediment yield, but there is a need for some clear recommendations on how to select the parameter set that defines a given land use scenario and on what is the most appropriate methodology to change the selected parameters when describing possible future conditions. This paper reviews the SWAT formulation to identify the parameters that depend on the land use, performs a sensitivity analysis to determine the ones with larger impacts on the model results and discusses ways to consider future land use conditions. The case study is the Atibaia river basin, with 2838 km2 (Sao Paulo, Brazil). The parameters identified by sensitivity analysis with the largest impacts on streamflow and sediment yield were the initial curve number for moisture condition II (CN), maximum canopy storage for each land use (CANMX) and the cover and management factor (USLE_C). The identification and appropriate parameter change can provide real estimates of the magnitudes in the land use changes, which were verified in this study. Such information can be used as an instrument for proposing improvements in the basin’s environmental quality and management.


Introduction
Understanding the implications of land use change is critical to the river basin planning and remediation efforts that are occurring all over the globe. The impacts of land use modification on river hydrological behavior need to be evaluated when reviewing current and future water availability, assessing river basin degradation processes or identifying the most appropriate measures to control the impacts caused by the disorderly land use change [1][2][3].
The effects of land use changes are empirically known but it is often very difficult to explicitly quantify them [3]. Geographic information system (GIS) and distributed hydrological models may contribute to overcome this challenge [4].
By offering the possibility to assess the behavior of a watershed under different scenarios, hydrological models may be used to predict the consequences of land use changes on the simulated hydrological processes [5]. The Soil Water Assessment Tool (SWAT) [6,7] is the focus of the paper, given its popularity for land use impacts assessment and because its formulation is based on physical parameters, most of which can in theory be measured in the field, therefore facilitating the description of different land use situations. Other popular hydrology models that may be used to assess land use Table 1. Soil Water Assessment Tool (SWAT) components.

SWAT Components Description
Hydrology Four main processes are considered: runoff, evapotranspiration, soil water movement and groundwater. All of these processes are accounted for in the model's water balance equation.

Weather
Climate is the main process inducing agent of the terrestrial phase of the hydrological cycle. The model requires daily data and monthly data of various meteorological variables.

Erosion/Sedimentation
Sediment yield is calculated for each HRU using the Modified Universal Soil Loss Equation. Vegetation cover and crop residues are considered when estimating soil particles detachment and transport.

Land use and Plant growth
Plant growth is simulated using an Erosion Productivity Impact Calculator (EPIC) simplification and occurs only on days when the average daily temperature exceeds a plant-specific base temperature [22].

Nutrients and Pesticides
Nitrogen and phosphorus movement and transformation is traced within the basin. In addition, pesticide loading, and bacterial contamination can also be computed.

Management practices
Crop cultivation, growth and grazing are simulated, as well as irrigation and nutrient and pesticide applications. Soil protection offered by vegetation throughout the year and the deposition of crop remains on the soil after harvest is considered.

Main channel processes
The movement of water, sediment, nutrients, and pesticides through the channel is computed. The channel dimensions are usually assumed constant, but there is an option to assume them dependent on erosion and deposition.
Water bodies In addition to channels, water may flow through four types of water bodies: ponds, wetlands, depressions, and reservoirs.
This study used SWAT 2012, released on 2018, integrating an ArcView GIS interface.

SWAT Land Use Parameters
Most studies found in the literature cite CN as the most important parameter of the SWAT model related to land use, which is used to calculate the runoff volume and infiltration. CN values range between 0 and 100, with higher CN values associated with areas with higher runoff potential such as urban districts. Low values indicate large retention and soil infiltration capacity and low runoff potential, occurring for example in forest areas [7,23].
The CN parameter is variable for each HRU and is computed from the combination of the hydrologic soil group and land use [24]. The combination table proposed by the SCS is stored in the model database and when new land use maps are loaded to the model, the model automatically determines new CN values for each HRU. Land use maps may be obtained from satellite images [10][11][12][13][14]25] or from hypothetically designed maps to predict possible future changes [3,26,27].
Local CN values change from event to event due to antecedent moisture conditions, as SWAT adjusts the daily CN according to the retention parameter that varies with soil profile water content or with accumulated plant evapotranspiration. We have decided to adjust CN according to the plant evapotranspiration because the alternative often results in over prediction of runoff in shallow soils. By calculating daily CN as a function of plant evapotranspiration, the value is less dependent on soil storage and more dependent on antecedent climate [28,29].
During model calibration and validation, CN values can be changed by directly altering the SCS original table, but this is usually not done. Alternatively, the database values are multiplied by an HRU dependent percentage change value, therefore ensuring the relative physical meaning of CN for different soils or HRUs [30].
The maximum canopy storage (CANMX) parameter is also mentioned by several studies as important in simulating the impacts of land use changes in the streamflow [12,13,31]. The CANMX is the maximum amount of water that can be stored in the canopy and trunks of fully developed trees [31]. Its value controls the density of plant cover so that it significantly affects infiltration and evapotranspiration.
A CANMX value is defined for each crop and a first estimate is stored in the model database. The canopy storage in each day depends on the leaf area index of the specific crop [7] (Equation (1)): where CANMX day is the canopy storage at day t (mm), LAI is the leaf area index and LAI mx the maximum leaf area index. LAI increases until the maximum leaf area index (LAI mx ) is achieved, then remains constant until onset of senescence, after which it declines to zero at harvest [32]. The LAI represents the structural properties of the plant canopy and impacts the exchange of energy and mass fluxes between the surface and the atmospheric boundary layer [33][34][35].
SWAT estimates LAI from the canopy height, h c (cm), the fraction of the plant's maximum leaf area index, f r LAImx and the fraction of potential heat units accumulated for the plant on a given day in the growing season, f r PHU [36]. The fraction of potential heat units accumulated by a given day is calculated from the heat unit accumulated on day i (heat units), HU, and the total heat units required for plant maturity (heat units), PHU. Heat units are calculated from maximum and minimum air temperature and from the plant-specific base temperature. No growth occurs for average temperatures at or below the crop base temperature [7,19].
The values of the parameters describing the behavior of each crop, such as LAI mx , h c,mx and PHU, are stored in the model database. This database can be added or changed if different crops must be considered or if there is information particular to the study area, but usually the crop specific parameters are not changed during model calibration. The value of the maximum canopy storage (CANMX) can however be changed and defined as a function of the soil cover. The values have been obtained and adapted from literature recommending CANMX value for different land uses [37,38].
The erosion/sedimentation sub model of SWAT uses the Modified Universal Soil Loss Equation (MUSLE) method [39] to determine soil erosion and sediment yield from each HRU. The simulation of the impacts of land use changes on sediment yield can be studied by changing the USLE C parameter, which reflects the protection given to the soil by surface cover [40,41]. High values of USLE_C represent landscape management practices that need to be improved, while lower values demonstrate that the management practices used are environmentally favorable, as is the case in forest areas [7,42,43].
SWAT updates USLE_C daily (Equation (2)) as plant cover varies during the growth cycle of the plant [7]: USLE_C day is the value for the cover and management factor for the land cover updates daily computed by SWAT; USLE_C is the value for the cover and management factor for the land cover given by the user; and rsd sur f is the amount of residue on the soil surface (kg/ha) contained on the top 10mm of soil, which depends on weather conditions such as precipitation, temperature, solar radiation, humidity and wind speed [44].
SWAT assigns default values for the cover and management factor for each crop, however these values were not consistent with the reality of the study area located in Brazil. The values of the USLE_C were assigned manually, based on studies carried out in the study region [36,[45][46][47].
SWAT performance was assessed using the following indicators: Nash and Sutcliffe efficiency [55], the percent bias (Pbias), and the coefficient of determination (r2).
The Nash and Sutcliffe efficiency (NSE) evaluates a normalized difference between the model results and the observed results (Equation (3)): where Y obs t is the observed value, Y obs t is the average observed value, Y sim t is the computed value, Y sim t is the average computed value. The NSE values ranges between −∞ and the optimal value of 1, with values less than 0 indicating that it is better to use the observed average than the values predicted by the model (Moriasi, 2007).
The percent of Pbias (PBIAS) indicates the average tendency of the simulated flows to be higher or lower than the observed flow, with values close to zero indicating a good adjustment of the simulated results to the observed data (Equation (4)): The coefficient of determination (r2) represents the proportion of the observed data variance that is explained by the model. The coefficient ranges between 0 and 1, with higher values representing a better performance of the model and less variation in the error [56].

Land Use Scenario Simulation
Land use change impacts on streamflow were evaluated by importing each land use scenario, and new land use maps were imported into SWAT. As CN, CANMX and USLE_C were setup to be solely dependent on land use, their values were automatically updated by the model.

Watershed Data
The Atibaia river basin has an area of 2837.3 km 2 , covering several municipalities of São Paulo and Minas Gerais states ( Figure 1). According to [57], the total population of the basin is approximately 372,456 habitants, served mostly from surface water sources.
The watershed altitude ranges from 509 m, in the West, to 2029 m, in the East. A 12.5 m resolution DEM, obtained from Advanced Land Observing Satellite-ALOS, was used to describe topography.
Currently, the river basin rural areas are mainly covered by pastures, sugar cane and mixed forests. In recent decades, agriculture activity has been growing, with and intensification of land use and the substitution of perennial and semi-perennial agriculture (pasture, orange, corn, coffee, etc) mainly by sugarcane [15,16].
The Lower Atibaia sub-basin has conditions of geology, pedology and slope that favors infiltration, namely thick soil profiles (>20 m) of sandy texture and low compactness. The Intermediate/Lower Atibaia sub-basin local conditions result in thinner alteration profiles (10 to 20 m) in slabby terrain and shallow profiles (<5 m) in the steeper parts of the basin. The Intermediate/Upper Atibaia sub-basin has thin soils (up to 5m), of medium to fine texture and low consistency with less noticeable infiltration conditions. The Upper Atibaia sub-basin has shallow soil profiles (<2 m thick) of fine texture and medium consistency, with limited infiltration conditions.
The hydrometric data were obtained from the databases of two government agencies: ANA (Agência Nacional de Águas) and DAEE (Departamento de Águas e Energia do Estado de São Paulo). Four hydrographic gages (4009, 3003, 3006 and 3007) and 17 rain gages ( Figure 1) were selected due to their location and data completeness. Data on temperature, humidity, wind speed and solar radiation was obtained from a weather gage located near to the basin and operated by ESALQ (Escola Superior de Agricultura Luiz de Queiroz).
Water 2020, 12, x FOR PEER REVIEW 6 of 19 m) in slabby terrain and shallow profiles (<5 m) in the steeper parts of the basin. The Intermediate/Upper Atibaia sub-basin has thin soils (up to 5m), of medium to fine texture and low consistency with less noticeable infiltration conditions. The Upper Atibaia sub-basin has shallow soil profiles (<2 m thick) of fine texture and medium consistency, with limited infiltration conditions. The hydrometric data were obtained from the databases of two government agencies: ANA (Agência Nacional de Águas) and DAEE (Departamento de Águas e Energia do Estado de São Paulo). Four hydrographic gages (4009, 3003, 3006 and 3007) and 17 rain gages ( Figure 1) were selected due to their location and data completeness. Data on temperature, humidity, wind speed and solar radiation was obtained from a weather gage located near to the basin and operated by ESALQ (Escola Superior de Agricultura Luiz de Queiroz). The basin main land uses are pasture, natural vegetation (seasonal semi-deciduous forest), and agricultural crops, such as orange, coffee, corn, soy, and sugarcane, followed by urban areas, reforestation areas (eucalyptus) and water bodies [62,63]. From 1990 to 2016, various crops such as coffee and oranges have been replaced by sugar cane, due to the great economic incentive for ethanol production.

Description of Land Use Scenarios
The studied land use scenarios were defined based on a comprehensive study done by the [64] that describes the current land use in the studied area and projects possible future land use changes. Scenario I represents the current situation with a predominance of forest and pasture in the rural areas of the Atibaia river basin. Based on this scenario, two possible future scenarios were studied.
Scenario II represents a possible land use evolution if the current trends persist, in which the urban and agriculture areas increase and forest areas recede [5,[64][65][66]. The urban area increases by 20%, while 50% of the native forest is replaced by agricultural areas and many crops are substituted by sugar cane.
Scenario III represents a desirable and more beneficial projection of future land use, while still accepting the same urban growth of the trend scenario. To achieve this goal, pasture and all range-grasses areas are substituted by forest areas. Figure 2 and Table 2 describe the land use distribution of each scenario in the main sub-basins of the Atibaia river basin. Figure 3 presents the land use maps used as inputs of the SWAT model. Table 3 presents the average CN, CANMX and USLE_C for each sub-basin and each scenario. The current conditions at the Lower Atibaia sub-basin are represented by an USLE_C of 0.001, indicating the absence of conservation practices in the watershed, mainly due to pasture and sugar cane areas without proper management. The conditions at the other sub-basins with smaller areas of sugar cane are better with an USLE_C within the range of 0.004-0.005. Future scenarios II and III assume an improvement of management practices. Scenario II represents a possible land use evolution if the current trends persist, in which the urban and agriculture areas increase and forest areas recede [5,[64][65][66]. The urban area increases by 20%, while 50% of the native forest is replaced by agricultural areas and many crops are substituted by sugar cane.
Scenario III represents a desirable and more beneficial projection of future land use, while still accepting the same urban growth of the trend scenario. To achieve this goal, pasture and all rangegrasses areas are substituted by forest areas. Figure 2 and Table 2 describe the land use distribution of each scenario in the main sub-basins of the Atibaia river basin. Figure 3 presents the land use maps used as inputs of the SWAT model.      Table 3 presents the average CN, CANMX and USLE_C for each sub-basin and each scenario. The current conditions at the Lower Atibaia sub-basin are represented by an USLE_C of 0.001, indicating the absence of conservation practices in the watershed, mainly due to pasture and sugar cane areas without proper management. The conditions at the other sub-basins with smaller areas of sugar cane are better with an USLE_C within the range of 0.004-0.005. Future scenarios II and III assume an improvement of management practices.

SWAT Performance
The sensitivity analysis indicated the 10 parameters that stream flow is most sensitive to, as measured by p-value and t-stat values ( Table 4). The t-stat is used to identify the relative significance of each parameter, with a larger absolute value meaning greater sensitivity. The p-value determines the significance of the sensitivity, and values close to zero indicate the most significant parameters [67]. Among these 10 parameters, CN and CANMX are the two land-use related parameters that significantly affect stream flow. The cover and management factor (USLE_C) from the MUSLE equation does not affect streamflow but is the most important parameter regarding the sediment yield. The initial curve number values for moisture condition II (CN) were reduced by 30% homogeneously across the basin to support base flow and infiltration, and available water capacity of the topsoil layer (SOL_AWC) was reduced by 30% uniformly across the basin to decrease the soil holding capacity, to delay the flow reaching the river and to increase the base flows. The soil evaporation compensation factor (ESCO) was reduced to 0.6 to promote evaporation from the deepest soil layers.
The groundwater "revap" coefficient (GW_REVAP) controls the water flow from the shallow aquifer to the unsaturated zone. This parameter was set to 0.04, which means more water available for the base flow. The deep aquifer percolation fraction (RCHRG_DP) was set to 0.12, which indicates that 12% of soil percolation water is directed to a deep aquifer.
The maximum canopy storage (CANMX) was defined between 0 and 80 to describe different canopy storage capacity to intercept the precipitation for each crop. The base flow alpha factor (ALPHA_BF) is used by the model to calculate the base flow and was adjusted to 0.001.
The slope length for lateral subsurface flow (SLSOIL), associated with the interflow source, was set at 40 m and the surface runoff lag coefficient (SURLAG), related to the daily surface runoff amount that discharges into the main channel, was set 1. Weighting coefficient for calculating retention dependent of plant evapotranspiration (CNCOEF) was adjusted to 1. 6 To reproduce sediment yield, the cover and management factor (USLE_C) was defined between 0 and 1 for different land uses. The USLE equation support practice parameter (USLE_P) was set to 0.8, recognizing the presence of conservation practices in the watershed.
The channel erodibility factor (CH_COV1) and the channel coverage factor (CH_COV2) were adjusted to 0.1 representing a low vulnerability to channel erosion. The sediment concentration in lateral flow and groundwater flow (LAT_SED) was set at 3000 to represent the sediment yield in lateral and groundwater flow.
The model calibration and validation process is influenced by the choice of the objective function and affected by the equifinality problem [68]. To alleviate these issues, a manual approach was adopted which considered the model's ability to reproduce streamflow and sediment yield at several monitoring stations. By carefully selecting the parameter values and considering multiple sites, two objectives and the equifinality problem are mitigated. A complete description of the calibration process and the parameters selected is presented in [20,21]. Table 4. Sensitivity analysis of SWAT model parameters before the calibration and validation process and for the final results.

Parameters
Before SWAT Simulation  Figure 4 compares the simulated and observed streamflow at the most downstream gage, obtained during calibration and validation. The model can adequately simulate the overall variation of stream, as well as both its minimum and the maximum values, although some extreme high daily values are overestimated. The Nash-Sutcliffe coefficient is usually higher than 0.5 when computed from daily values and higher than 0.7 when computed from monthly values ( Table 6). The model results are slightly worse for the upstream sub-basin. Figure 5 compares observed sediment yield values with the model simulated results, showing that SWAT can replicate the scale and the variation pattern of the sediment yield, although not with precision. Table 7 presents the computed performance indicators for sediment yield. The low monitoring frequency of sediment transport (bimonthly) and the fact that most data has been collected during low flows hinders a detailed evaluation of the model performance.
Based on these streamflow and sediment yield results, for both the calibration and validations periods, we concluded that the model is adequate for estimating the impacts of different land use scenarios. a Multiplying factor to be applied to the parameter original value. Values adapted by (1) [46] (2) [47] (3) [45]. (4) [36].      The conversion of forests to agriculture, sugarcane and urban uses (scenario II) leads to an increase of streamflow, as evapotranspiration decreases and surface flow, percolation and

Streamflow
The conversion of forests to agriculture, sugarcane and urban uses (scenario II) leads to an increase of streamflow, as evapotranspiration decreases and surface flow, percolation and groundwater flow increases (Table 8). Conversely the replacement of pasture, range-grasses and agriculture by forest (scenario III) leads to an increase in evaporation and surface flow and to a decrease of streamflow, lateral flow, percolation and groundwater flow. This result is in line with findings that show that forest development promotes infiltrations and reduces surface streamflow [14,[69][70][71][72][73][74][75].
While land-cover change may have a moderate impact on average annual flow, it can significantly influence seasonal and monthly streamflow. Figure 6 presents the average monthly streamflow estimated at the four hydrometric stations for each land use scenario. The precipitation seasonal variability of precipitation leads to high flows from December to April and to low flows from July to October.
In general, scenario II shows an increase in the streamflow which is higher during wet months. The increase of agricultural and urban areas and the decrease of soil cover associated with this scenario (Table 4) results in higher average CN values and lower CANMX values. This trend is mainly associated with substitution of forest areas by agriculture. In turn, these changes lead to a faster water movement throughout the river basin, a lower retention time, a decrease of evapotranspiration and to an increase of surface and lateral flow (Table 9). This change of the runoff-precipitation ratio is more striking in the rainy periods and less in the dry periods.
The significant reduction in evapotranspiration in the upper sub-basin is mainly due to the replacement of forests by agriculture. The forest areas (Forest mixed-FRST) have a Leaf Area Index (LAI) value of 5 mm and the canopy height (h c,mx ) of 6 cm, while the agricultural areas (Agriculture Land Generic-AGRL) have a LAI of 3 mm and hc of 1 cm. The reduction in the value of these parameters influences the canopy storage and contributes to the decrease in the precipitation interceptions and evaporation in all sub-basins.
In scenario III, the increase of forest areas over the pasture and range grasses areas, results in lower CN values and higher CANMX values, mainly due to increases in the LAI and h c,mx values, respectively from 4 mm and 0.5 cm for pasture and 2.5 mm and 1 cm for range grasses to 5 mm and 6 cm for forests.
However, the simultaneous trend towards larger urban areas counterbalances the increase of evapotranspiration and decrease of lateral flow and streamflow due to forestation, and scenario III becomes very similar to scenario I. The small increase of streamflow in January and the small decrease in the dry season small change is due to the increase in temperature during the wet periods that has significant impact on the growing stage of the vegetation in function of the evapotranspiration, and consequently in this period the base flow becomes important to contribute to the streamflow.
These results are in line with the studies performed by [3,13,76].  Table 9 compares the annual average sediment yield estimated from each scenario for the 1990-2016 period.

Sediment Yield
The replacement of forest areas by perennial agriculture (scenario II) leads to an increase in annual average sediment yield values (+24%) because the original forest areas ensure a superior soil cover over a longer period of the year. Even when perennial agricultures are grown with good conservation practice, these areas produce high amounts of sediment due to their relatively low soil cover, as compared to forest areas. This result was also obtained by [77,78].
Scenario III shows a slight increase in annual average sediment yield values (+2%). Despite the increase of forest areas that ensures soil cover, the significant urban growth contributes to the increase in surface flow and, consequently, to a larger amount of sediment that is transported into rivers. According to [79], the increased runoff in urban areas is very effective at eroding the available sediment sources.  Figure 7 plots the average monthly sediment yield for each scenario. The sediment yield behavior is very similar to the streamflow because the SWAT model uses to the MUSLE equation, which computes the sediment yield from runoff volume and peak flow rate [80]. While the USLE and RUSLE models estimate average annual gross erosion as a function of rainfall energy, MUSLE uses a runoff factor that is recalculated every day, which is appropriate for the simulation of erosion and sediment yield within each HRU at a daily time step [7,[80][81][82]. Sediment yield prediction also is improved because runoff is a function of antecedent moisture condition as well as rainfall energy [38,83]. The need for delivery ratios required by USLE is eliminated because the runoff factor  Table 9 compares the annual average sediment yield estimated from each scenario for the 1990-2016 period.

Sediment Yield
The replacement of forest areas by perennial agriculture (scenario II) leads to an increase in annual average sediment yield values (+24%) because the original forest areas ensure a superior soil cover over a longer period of the year. Even when perennial agricultures are grown with good conservation practice, these areas produce high amounts of sediment due to their relatively low soil cover, as compared to forest areas. This result was also obtained by [77,78].
Scenario III shows a slight increase in annual average sediment yield values (+2%). Despite the increase of forest areas that ensures soil cover, the significant urban growth contributes to the increase in surface flow and, consequently, to a larger amount of sediment that is transported into rivers. According to [79], the increased runoff in urban areas is very effective at eroding the available sediment sources.  Figure 7 plots the average monthly sediment yield for each scenario. The sediment yield behavior is very similar to the streamflow because the SWAT model uses to the MUSLE equation, which computes the sediment yield from runoff volume and peak flow rate [80]. While the USLE and RUSLE models estimate average annual gross erosion as a function of rainfall energy, MUSLE uses a runoff factor that is recalculated every day, which is appropriate for the simulation of erosion and sediment yield within each HRU at a daily time step [7,[80][81][82]. Sediment yield prediction also is improved because runoff is a function of antecedent moisture condition as well as rainfall energy [38,83]. The need for delivery ratios required by USLE is eliminated because the runoff factor represents energy used in detaching and transporting sediment [7].
Ref. [81] show a comparison of the erosion prediction models from USLE, MUSLE and RUSLE in a Mediterranean watershed, in the case of Wadi Gazouana (NW of Algeria).

Conclusions
The effects of land use change on the streamflow and sediment yield of the Atibaia river basin were estimated using SWAT. The model parameters that have the largest impacts on streamflow and sediment yield are the initial curve number for moisture condition II (CN), the maximum canopy storage for each land use (CANMX) and the cover and management factor (USLE_C). Other parameters also related to land use are the maximum leaf area index (LAImx), the canopy height (ℎ , ), the total heat units required for plant maturity ( ), and the amount of residue on the soil surface ( ). The model was setup to automatically update these parameters from each land use map.
Two different land use change scenarios were applied to the study basin and the streamflow and sediment yield outputs were compared with the current situation. The expansion of perennial agriculture and urban areas at the expense of forest areas leads to an increase of streamflow, as percolation and groundwater flow increases and evapotranspiration decreases. The increase of surface runoff and streamflow leads to an increase of sediment yield. The expansion of forest areas over pasture, range-grasses and agriculture leads to a decrease in streamflow and to a slight increase in sediment yield as evaporation and surface flow increase and the lateral flow, the deep percolation, the groundwater flow and the percolation decreases. Both land use scenarios maintained the seasonal variation of streamflow and sediment yield, following the precipitation pattern in the wet and dry periods.
The identification and appropriate parameters change in the SWAT model can provide real estimates of the magnitudes in the land use changes, which were verified in this study. Such information can be used as an instrument for proposing improvements in the basin environmental quality and management.
It should also be noted that land use changes often lead to other significant alterations of the river basin, namely in the drainage network. A trend towards urbanization brings along changes in the river channels alignments, profiles, cross-sections, bed and bank materials and conveyance capacity. In general, the artificialization of the river network speeds the water velocity through the river network, leading to higher peak flows. A trend towards a renaturalization of the river basin may lead to changes in the opposite direction. The impacts of these alterations may be significant, and although they can also be simulated in SWAT by modifying the flow routing parameters, this was not performed in the current study.

Conclusions
The effects of land use change on the streamflow and sediment yield of the Atibaia river basin were estimated using SWAT. The model parameters that have the largest impacts on streamflow and sediment yield are the initial curve number for moisture condition II (CN), the maximum canopy storage for each land use (CANMX) and the cover and management factor (USLE_C). Other parameters also related to land use are the maximum leaf area index (LAI mx ), the canopy height (h c,mx ), the total heat units required for plant maturity (PHU), and the amount of residue on the soil surface (rsd sur f ). The model was setup to automatically update these parameters from each land use map.
Two different land use change scenarios were applied to the study basin and the streamflow and sediment yield outputs were compared with the current situation. The expansion of perennial agriculture and urban areas at the expense of forest areas leads to an increase of streamflow, as percolation and groundwater flow increases and evapotranspiration decreases. The increase of surface runoff and streamflow leads to an increase of sediment yield. The expansion of forest areas over pasture, range-grasses and agriculture leads to a decrease in streamflow and to a slight increase in sediment yield as evaporation and surface flow increase and the lateral flow, the deep percolation, the groundwater flow and the percolation decreases. Both land use scenarios maintained the seasonal variation of streamflow and sediment yield, following the precipitation pattern in the wet and dry periods.
The identification and appropriate parameters change in the SWAT model can provide real estimates of the magnitudes in the land use changes, which were verified in this study. Such information can be used as an instrument for proposing improvements in the basin environmental quality and management.
It should also be noted that land use changes often lead to other significant alterations of the river basin, namely in the drainage network. A trend towards urbanization brings along changes in the river channels alignments, profiles, cross-sections, bed and bank materials and conveyance capacity. In general, the artificialization of the river network speeds the water velocity through the river network, leading to higher peak flows. A trend towards a renaturalization of the river basin may lead to changes in the opposite direction. The impacts of these alterations may be significant, and although they can also be simulated in SWAT by modifying the flow routing parameters, this was not performed in the current study.