The Protective Role of a Riparian Vegetation to the Sources of a Rural Watershed

A recent review of the Brazilian forest law has suggested the possibility of decrease the riparian buffer width along the banks of waterways. This vegetation can trap solutes which had been deposited in the soil and which would otherwise be carried to the waterways in runoff water. In the current study, we applied the AVSWAT model to the Rio Bonito micro-basin of city of Descalvado, State of Sao Paulo, Brazil, to verify the amount of nitrogen (N) and phosphorous (P) entering the water resource, and to evaluate the environmental protection provided by the riparian vegetation zones. Our results demonstrate this protection is effective in reducing the amount of nitrogen (N) and phosphorous (P) carried to the river. The scenario produced by the data could be used to diagnose the contamination of this particular watershed and to indicate the measures that should be adopted to ensure the restoration and preservation of riparian vegetation zones.


Introduction
One of the greatest challenges we face in the 21 st century is how to conciliate the growing world population and its demand for ever more comfortable lifestyles, with the scarcity of natural resources and the problem of pollution particularly as the growth in world population is followed by increased land occupation and more intense production.
Pollution of water resources due to rural economic activities is especially complex because it is often difficult to identify the source of the pollution. This type of pollution is essentially brought about by major hydrologic events and associated not only with the use and occupation of the soil, but also with the pedological characteristics of the particular watershed. Improper land use structure and management will lead to soil erosion and nutrient loss with surface runoff, resulting in a large area of non-point source pollution in the watershed causing loss of nitrogen and phosphorus and water eutrophication [1]. When erosion occurs in agricultural areas, the eroded soil can carry nutrients, but also fertilizers or pesticides [2]. Rainwater carries both organic and inorganic materials that were deposited in the soil and may be soluble in the water sources, significantly increasing its load of pollutants.
Amongst the various chemical elements carried to the waterways, nitrogen (N) and phosphorous (P) can be highlighted as being particularly important, in that they are directly involved in the metabolism of the aquatic life.
In agricultural areas, the riparian vegetation plays an important role in the structure and function of stream ecosystems [3] and can greatly improve the quality of the water of nearby agricultural streams by reducing nutrient leaching in groundwater and surface water runoff, by trapping sediments in surface runoff and thereby decreasing organic or inorganic material loading in surface water even when they comprise little of the total catchment area [4]. The load of nitrate can be effectively reduced by about 30% by establishing a watershed vegetation buffer [5].
In Brazil, changes have recently been made to the Brazilian federal forest law which sets standards for the protection and use of forests and other forms of vegetation. These changes endanger the water quality of the sources, because they prescribe much less maintenance or reforestation of riparian vegetation than was the case in previous legislation. Prescribing a consistent riparian buffer width has been difficult because there are considerable variations in slope, vegetation composition, soils, etc., determining the level of acceptable risk. Hence, the importance of maintaining critical riparian zone structural and functional attributes relies on the intensity of management activities and the resource(s) needing protection [6]. As a result of inadequate or incomplete knowledge of riparian zone width, farmers are often faced with the challenge of making decisions about appropriate riparian zone protection without a science-based understanding of the structural and functional characteristics that define the riparian area in order to protect stream water quality [7].
In this scenario of uncertainty arising from the recent change in Brazilian federal forest law, the hydrosedimentological model-AVSWAT (ArcView Soil and Water Assessment Tool) has emerged as an important tool for the analysis and evaluation of the environmental protection that riparian zones provide to our watersheds.
We are unaware of any studies that have examined the entry of chemical elements in the sources of the micro basin of Rio Bonito of city of Descalvado, State of São Paulo, Brazil. Thus, this study aims to evaluate the scenario simulated by the AVSWAT model, considering the contribution of nitrogen (N) and phos-  [8].
For modeling purposes, AVSWAT subdivides the basin based on its topography, soil type and land use, and thus preserves the spatially distributed parameters of the entire basin, as well as their homogeneous characteristics. The standard process used to divide a basin into sub-basins is to specify boundary areas, which are the minimum areas requiring drainage to a point, forming a channel [9]. Each sub-basin can be parameterized by AVSWAT using a series of Hydro-  It presents two options: either simulating the nutrients not considering the transformation processes that occur on the river or considering such transformations.

Monitoring Routine of Nitrogen (N) and Phosphorous (P) Done by the Model
Nitrogen, essential for plant growth, appears in three forms in the soil: organic nitrogen in the humus; mineral forms bound on colloids; and in solution. It can be added to soil by fertilization, fixation by bacteria and rain, and it can be removed by plants, lixiviation, volatilization, denitrification and erosion.
The AVSWAT monitors five forms of nitrogen in soil: two inorganic forms ( 4 NH + and 3 NO − ) and three organic forms (recent-associated with vegetables and biomass residues, stable and active-associated with humic substances). The runoff to the main channel carries organic nitrogen bound to soil particles, which is associated with the sediment load in the HRU. Changes in the sediment load will be reflected directly in the organic nitrogen load. The total of organic nitrogen carried along with sediments to the channel is calculated by the loading function modified by [11]: where: orgN surf is the total of organic nitrogen transported to the main channel by the runoff (Kg N/ha); conc orgN is the organic nitrogen concentration in the first 10 mm of topsoil (g N/ton de solo); sed is the sediment load on a given day (ton); area hru is the HRU area (ha); Σ N is the nitrogen enrichment ratio, The organic nitrogen concentration in the topsoil, conc orgN , is calculated by the function: Figure 1. SWAT soil nitrogen pools and process that move nitrogen in and out of pools Source: modified from [10]. The AVSWAT monitors six forms of phosphorous in the soil: three inorganic forms (stable, active and in solution) and three organic forms (recent-related to plant residues, stable and active-associated to humic substances). The AVSWAT simulates the processes of mineralization, decomposition, immobilization, absorption of inorganic phosphorous and lixiviation. As with nitrogen, the organic and mineral forms of phosphorous bound to soil particles are carried by the runoff to the main channel; maintaining the same relation to the sediment load generated by the HRU as described for the nitrogen, i.e., changes in the sediment load will reflect directly in phosphorous load.
where: minP act,surf is the total of phosphorous (kg P/ha) in the active mineral form in the first 10 mm of topsoil; minP sta,surf is the total of phosphorous (Kg P/ha) in the stable mineral form in the first 10 mm of topsoil; orgP hum,surf is the total of phosphorous (kg P/ha) in the humic organic form in the first 10 mm of topsoil; orgP frsh,surf is the total of phosphorous (kg P/ha) in the organic form (fresh) in the first 10 mm of topsoil; ρ b is the density (Mg/m 3 ) of the topsoil; depth surf is the depth of the topsoil (10mm). The micro basin is distinguished by its intense poultry industry including areas of high farm concentration, being one of the most important chicken an eggs producers in the Brazil [12] [13]. The farms are concentrated in the following eight sub-basins generated by the model: 4, 10, 20, 22, 24, 28, 30 e 39 ( Figure  5), basically following the path of the main channel of the micro basin ( Figure  4). This pattern suggests the continued practice of waste and sewage disposal into the aquifers, a fact that seriously damages the quality of water in many areas of agricultural production in Brazil, as shown in the case of the State of Santa Catarina, and more recently, the State of Goiás, [14] [15].  Each one of the forty one sub-basins of the micro basin, numbered as show in Figure 4 and Figure 5, has a different use and activity. However this can be related to economic exploration of soil, there is a predominance of sugar cane, orange orchards, and pastures. The micro basin soil types have a predominance of the Latosoils group, which is predominantly clay as well as others with some significant sand content, which are not the most susceptible regarding erosion (an environmental aspect that is associated with the nutrients' supply in water sources) that may occur due to the use of land already explored, and to the slope of the area, as shown in Figure 6.

Characterization of the Study Area
This information plan, containing the ranges of slopes was obtained by calculating the continuous values of these slopes through interpolation and data processing of the planialtimetric map, reclassifying them at intervals adapted to the classes of slope limiting agricultural work, allowing a preliminary qualitative evaluation of the susceptibility to erosion, which can result in a silting-up process. It can be seen in the Figure 6 that the south end of the basin is the most susceptible area and that some of the streams are in thalwegs, whose protection demands greater care.

Model Calibration and Statistical Analyzes
In order to calibrate the AVSWAT model, we considered several sources such as [16] and [17], and also conducted some trial and error tests to identify the variables which appeared to be more sensitive. As a result, the following variables   Universal Soil Lost Equation, CN, ALPHA_BF, SOL_Z, SOL_NO3, the later four being the most relevant to the study. [17] suggests several methods to evaluate the calibration of a model and later verify it: the mean absolute error; the determination of coefficient R 2 (that indicates what percentage of total variation is explained by linear regression); the Nash-Sutcliffe coefficient and the volumetric difference between the simulated and observed outflows. In the current paper, following [16], we used the Nash-Sutcliffe efficiency coefficients, COE (Equation (5)) in addition to the simulated data deviation related to the observed data, Dv (Equation (6)), which we used to compare some of the parameters obtained in the basin.   (5) in which E m is the observed event, E s is the simulated event created by the model; E is the mean observed event in the simulation period, and n is the number of events.
The deviation is presented by Equation (6): in which E is the observed event in the analyzed period, and E * is the simulated event in the period. These two parameters were used to evaluate the model performance.
Results showed that the model had a good estimative potential in relation to the observed data. In general a COE varying from 0.70 to 0.80 indicates a good fit [18]. The COEs for nitrogen and phosphorous were 0.76 and 0.74, respective- cates that the simulated amounts of nitrogen and phosphorous production were higher than the observed ones.

Nitrogen (N) and Phosphorous (P) Volumes Generated from Animal Dejects
We used the methodology proposed by [19] to calculate the volumes of nutrients (N and P) arising from poultry waste in different farms within the micro basin. The AVSWAT model considers the monthly fertilizer data, meaning that the final values of nitrogen (N) and phosphorous (P), shown in Table 1, were divided by 12, thus representing the soil fertilization in a homogeneous manner.
The use of this methodology for the nutrients data entry does not influence the model simulation results, since they were parameterized in annual values.
In order to input the fertilizer data, we used each HRU on each sub-basin where there was a farm, i.e., the total volume of animal waste generated in the sub-basins was proportionally distributed over each HRU in this same sub-basin.
To meet the entry unity required by the model (values of fertilizer entry expressed in kilograms per hectare-Kg/ha), the area of their respective HRUs was divided by the total volume of nutrients for each sub-basin; the results are expressed in Kg/ha.

Simulation of Nitrogen (N) and Phosphorous (P) Entry into the Source
For the purpose of simulation, we used only the results obtained for the nitrogen (N) and phosphorous (P) nutrients in the organic phase, since animal waste provids them "in natura".  In this research, we did not consider mineral forms resulting from the biochemical processes that the organic matter undergoes.
To analyze the effect of fertilization when using animal waste, we simulated the nutrients entries into the sources in two phases. First, we simulated the nutrient's entry without considering the waste entry, so that the model would consider the composition of each soil type and its characteristics nutrients. In the second phase, we went on simulating the nutrient's entry into the channel, considering soil fertilization using animal waste.
In this sense, the Table 2 and Table 3

Discussion
The relation between the nutrients carried to the channel and the percentage of the sub-basin covered by vegetation demonstrates the protective function of this vegetation. This indicated that, in some sub-basins, this area was responsible for the retention and accumulation of sediments, and therefore nutrients.
It can be inferred that the runoff slows when it reaches areas of dense vegetation. This is characteristic of the forests, savannas and reforestation areas that   It was noted, however, that for sub-basin 30, the vegetation was shown to be ineffective in the control of the supply of nutrients to the sub-basin's channel.
This fact is explained by two simultaneous conditions: the first condition, is due to the high concentration of nitrogen (N) and phosphorus (P) released from animal wastes, since it is the only sub-basin that has poultry and pork farms; the second one is the steep slopes found there, as noted in the item relating to the characterization of the study area (the southern end, where the sub-basin is a significant portion of its territory with slopes greater than 20%-ranges "20% -45%" and ">45%").
One likely reason for the heterogenicity of the protective function of riparian vegetation in all sub-basins is the structure of this barrier, which in some cases is little more than remnants or fragments reducing its effectiveness, besides the possibility of being distributed in areas of steep slopes.
In analyzing the results, the sub-basins 4, 20, 10 and 30 require special attention, since they are the most important providers of nitrogen (N) and phosphorous (P) in the watershed. In these areas, one must carefully evaluate fertilization practices, especially the predominant industries of pasture and sugar cane, since other solutes and pollutants could reach the sources of watershed.
Likewise, it is recommended that these critical areas should be inspected and monitored in order to confirm that farmers are fulfilling their legal obligation to ensure a barrier of riparian vegetation in the vicinity of sources most directly affected.
Finally, our hope is that our results will help to form the basis of discussion on the protective role of riparian vegetation, a role that has been underestimated by the recent reform of the Brazilian forest law.

Conclusions
A characteristic on the area surveyed is that the runoff, initiated by the erosion processes, slows when it reaches areas of vegetation along the banks of rivers, retaining the sediments and nutrients that otherwise would be carried by the water to the sources of the sub-basin.
The protective function of riparian vegetation may be compromised by the presence of steep slopes on the ground hence this represents a spatial configuration of the terrain that can increase runoff and nutrient load to the sources.
The existence of areas of riparian vegetation in the watershed has a direct correlation to the nutrient levels carried to the water sources.
This supports the literature on the subject and strengthens the role of this type of vegetation in protecting water resources, generating knowledge to decision makers that define the landscape involving riparian vegetation.