Geoprocessing Applied in the Estimate of Infiltration and Surface Runoff in HPP’s Contribution Watershed

The knowledge of the hydrological regime contributes to the prediction of periods of higher or lower water supply as a function of the frequency and volume of rainfall in relation to the pedological, land cover, and land use conditions. This work aims at evaluating the spatio-temporal variability of the infiltration and runoff potential of an area formed by 4 sub-basins tributary to the São Marcos River in Cristalina (GO). The SCS/NRCS-USDA (2004) method was used, using geoprocessing resources, considering soils in dry, semi-humid and humid conditions. The estimates indicate that rainfall with an average of 6.83 mm, in the September-October transition, in the face of soils in the dry condition does not provide surface runoff. Thus, the variability of pedological characteristics associated with the various conditions of coverage and land use indicates that even the areas with lower infiltration potential begin to contribute effectively with their respective channels in the transition from October to November, when rainfall events are more frequent and significant, and the soils are already in conditions of greater humidity.


Introduction
The knowledge about the hydrological behavior of hydrographic systems allows planning and management, as well as better equating the use of water resources, energy generation [2]. This is a complex but necessary equation; since observing the climatological and hydrological conditions of each area, it seeks to contemplate the most diverse activities, given the importance played by each of them.
The accentuated climatic seasonality prevailing in the Midwest region of Brazil results in a marked asymmetry regarding the distribution and availability of water in river systems, with a notable reflection on the water supply throughout the year [3]. This implies a greater concern about the volume available in the springs and, consequently, about the capacity of surface and subsurface reserves to meet the demands that tend to increase during the dry season. In this sense, it is necessary to evaluate the influence of climate seasonality on the rate of soil wetting in view of the increasing demands for water resources in order to reduce the effect of the asymmetry of rain supply in the behavior of hydrological parameters, especially the contribution of surface and subsurface flows to the flow of water sources [4].
The integration of climatological, geomorphological, pedological, land cover and land use data and variables using models operated in a geoprocessing environment represents one of the major advances in the field of modeling environmental processes and phenomena. In recent decades, the association of mathematics with the advent of computing has resulted in substantial advances in georeferenced spatial data processing technology [5]. By contemplating the spatial and temporal variability of environmental constraints at different scales or levels of detail, these technologies allow the elaboration of true hydrological radiography about the interaction of rainwater with soil and terrain properties. Consequently, they enable a dynamic approach to the behavior of infiltration and runoff [6]. Thus, based on the behavior of the spatial variability of the conditioning factors, it is possible to establish zoning proposals in which it is possible to ensure patterns of hydrological behavior and, consequently, work with greater predictability of water supply.
The relationships involving rainfall, soil infiltration and runoff on the ground are very complex because besides not being fully amenable to direct observations in the field, they vary according to the conditions of each environment [7] [8]. However, the analysis of the continuous distribution of the main variables and their correlation with the behavior of the processes directly monitored in the field and under specific environmental conditions, allows representative models to be developed and applied with reasonable levels of precision and accuracy [9] [10]. Among these models, the CN (Curve Number) method proposed by the SCS (Soil Conservation Service), currently NRCS-USDA (Natural Resources Conservation Service-United States Department of Agriculture), stands out, which since 1954 has undergone several improvements and with wide applica-tion in several countries [11]. With this, it has become a widely used model as it is based on the amount of rainfall resulting from a given rainfall event, the infiltration capacity of soils and losses resulting from local environmental factors such as microrelief, as well as other losses of lesser impact, having presented good results in applications in Brazil [12].
Thus, diagnoses resulting from hydrological studies constitute the starting point for the definition of mitigating or mitigating actions, insofar as they allow the identification of the factors that condition critical situations and, consequently, locate the environments in which such situations predominate. In environments where rainfall supply tends to be scarce or of markedly asymmetric distribution, a better understanding of the hydrological functions performed by soils in relation to the storage capacity and maintenance of the perenniality of springs and drainage channels is necessary, even in periods of drought [13] [14].
Thus, the integrated assessment of hydrological parameters in areas of direct influence of reservoirs, whatever their use, allows, for example, greater predictability of hydrological conditions, which provides more efficient management of surface water resources, observing the existing dependence between rainfall regime and behavior of water in contact with the soil [15].
In this sense, the present work aims to evaluate the spatial and temporal variability of the behavior of the hydrological parameters infiltration potential, effective precipitation and runoff under rainfall conditions representative of the pre-rainy season, beginning of the rainy season and full rainy season in the contribution areas of the right bank basins of the São Marcos River in the municipality of Cristalina, state of Goiás, Brazil. The knowledge about the spatial and temporal variability of the behavior of these parameters is of great importance for a better understanding of the relationship between rainfall occurrence and volume and the response time in reservoir level rise. This relationship operated by geoprocessing resources provides a greater quantity of calculations and, consequently, greater spatial and temporal detailing of the results. This provides a more detailed understanding of the interdependence of conditioning factors, especially rainfall intensity and duration, permeability and porosity of soils, the influence of relief and land cover and use on the hydrological balance.

Location and Characterization of the Study Area
The study area comprises 4 sub-basins on the right bank of the São Marcos River in the municipality of Cristalina, state of Goiás, Brazil, and that drain directly into the reservoir of HPP Batalha-FURNAS, as shown in Figure 1. It comprises an area with the predominance of agriculture, especially irrigated, with emphasis on the high number of irrigation pivots [16]. The higher and flatter parts are  Alluvial Deposits composed of sand and gravel [17]. Considering the World Reference Base for Soil Resources [18], in the higher and residual parts the Plintossols (Petric) predominates, followed by lower and flatter areas with the pre-

Methodology
The methodology consisted in evaluating the effects resulting from the occurrence of rainfall events in view of the conditions of relief, soil type, cover and E. D. Nunes et al. Journal of Geographic Information System use, as well as antecedent moisture conditions (AMC) in the spatial variability of the estimates of infiltration potential, effective precipitation and peak flow. For this purpose, the approach occurred at the slope scale contemplating the specific contribution area to its respective drainage channel. For the pre rainy season situation (AMC I), the soil was considered to be in the dry condition where the accumulated rainfall in the last 5 days is <13 mm. For the following period, beginning of the rainy period and of greater regularity of rainfall (AMC II) the soil was considered in the intermediate condition in which the accumulated rainfall in the last 5 days is ≥13 and ≤than 53 mm. For the full rainy period (AMC III) the soil was considered to be in a humid condition or close to saturation where the accumulated rainfall in the last 5 days is >53 mm.

Spatial Variability of Effective Precipitation Estimates
The effective precipitation was determined based on the proposal of the Soil Conservation Service (SCS), now Natural Resources Conservation Service (NRCS) of the U. S. Department of Agriculture, with its first version published in 1954 [23]. This proposal has undergone several revisions, one of them in 1993, and has been widely applied, including by the Federal Highway Administration (FHA-U.S. Department of Transportation) [24], consisting of the ratio between the rainfall resulting from the rainfall event and the infiltration capacity of soils, as presented in Equation (1).
where: Pe = effective precipitation, or resulting water sheet height, in mm; P = precipitation or rainfall height resulting from a given event, in mm; S = soil infiltration potential, in mm; and Ia = 0.2S abstraction or initial loss considered.
As observed in the first part of Equation (1), in the effective precipitation an initial loss corresponding to 20% of the infiltration potential is considered, due to interception by vegetation, retention in micro relief, as well as other forms of wetting of the surface environment. This implies that rainfall events with total volume less than 20% of the soil infiltration potential do not provide surface runoff.

Estimation of Intensity, Duration, and Accumulated Rainfall Volume
The values of total precipitation accumulated in the rainfall event in each of the three seasonal situations considered were determined by analyzing the average values of precipitation over the last 15 years from the meteorological station code A036 operated by the National Institute of Meteorology in the municipality of Cristalina-GO. For the condition preceding the beginning of the rainy season, the precipitation was determined based on the average of rainfall events that occurred in late September and early October, which are around 6.83 mm. For the beginning of the rainy season and more regular rainfall, the rainfall was determined based on the average of rainfall events in the second half of October, Journal of Geographic Information System which were around 22.83 mm. For the rainy season, the rainfall height was determined by means of the relationship Intensity-Duration-Frequency (IDF) [25], using Equation (2) where: i = average of maximum rainfall intensities, in mm/h; T = return time considered, in years; t = time of concentration, in min; and K, a, b and c adjustment coefficients specific to the climatological station.
For this purpose, the rainfall time was considered equal to the longest runoff concentration time found in the basins. The rainfall height resulted from the product of the precipitation intensity and the duration of the rain. Since this is not a study related to safety works, a return time of 25 years was considered in the present proposal.
For the concentration time we used the relationship between flow or basin length and average slope [27] using Equation (3).
where: Tc = runoff concentration time, in min; Lw = length of the basin main flow line, in km; and S = average basin slope in m/m.

Determination of Runoff Number (CN) and Infiltration Potential (S)
The S values were estimated based on the CN (runoff number) values, according to the Natural Resources Conservation Service (NRCS-USDA) [23], as presented in Equation (4).
where: S = infiltration potential, in mm; CN = Curve Number, dimensionless; and 25,400 and 254 are constants originating from the model. For this purpose, the CN values were determined based on the evaluation of soil types-especially texture and depth-and its classification into Hydrological Groups (HG) as presented in Table 1 [28] [29] the conditions of cover and use and the antecedent humidity. In this method, the CN values vary from 0 (low flow capacity) to 100 (high flow capacity). It is noteworthy that for this calculation the soil was considered in three conditions of antecedent moisture: 1) the situation preceding the beginning of the rainy period with the soil in the dry condition (5-day precipitation < 13 mm = AMC I); 2) the one corresponding to the beginning of the rainy period and greater rainfall regularity, with the soil in the intermediate condition (13 mm ≤ 5-day precipitation ≤ 53 mm = AMC II); and 3) referring to the full rainy period, with the soil in the wet condition (5-day precipitation > 53 mm = AMC III).
An example of the association of land cover classes and land use conditions with the Hydrologic Groups and their resulting runoff number values can be Journal of Geographic Information System

Estimation of Runoff Concentration Time and Peak Flows
The peak flow rate was determined as proposed in Equation (7)

Flow Length, Concentration Time and Contribution Area
Considering Relating the length of the strand with the length of the channels, it can be seen that the channels are up to 15 times more extensive than the flows lines developed along the strand.
A direct consequence of this behavior is the tendency of simultaneous contribution of the area of the lower order basins to their respective channels. However, considering the elongated character of the higher order basins, the influence of the hydrologic effects resulting from each sub-basin tends to occur in a gradual way [32]. These characteristics, associated with the low altimetric gradient along the main channels, confer to these basins low probability of occurrence of peak flows or overload in the main channels [33] [34]. Consequently, the high morphometric asymmetry between the contribution area at the slope level and the area of the higher-order basins suggests the need for long-lasting and higher-volume rainfall events so that rainfall can contribute, via hydrographic subsystems, to the reservoir.
The relief is distributed in three main altimetric levels, with slopes ranging from 0% to 8%, which transition from plateaus to pediments and from these to plains, by means of amphitheaters. The highest slopes occur in the transition zones between surfaces, which are small and rarely range between 20% and 43%.
Although in the transition between one plateau and another slopes that can reach 43% predominate, these are of little occurrence and the effects resulting from the high runoff coefficients developed upstream tend to be attenuated by the high infiltration potentials prevailing along the plains.

Rainfall Height, Hydrologic Groups, Cover and Use Conditions, and Behavior of Hydrologic Parameters
Taking as a reference the adopted methodology, the average precipitation of rainfall events that precede the rainy season, as well as the association of hydrological groups with the various conditions of land use, it can be seen that the predominant rainfall at the end of September and beginning of October is not sufficient to produce runoff. As observed in Figure 3,    In these conditions it is already possible to verify the occurrence of effective precipitation with a height of up to 7.4 mm, especially in the steeper areas and less developed soils, such as Cambisol and, mainly, Leptsol, both with low infiltration potential. As for the area corresponding to the fluvial plain, the effective precipitation is null due to the topographic, pedologic and cover and use conditions, unfavorable to the formation of surface runoff. Similar behavior can be verified for surface runoff estimates, which already reach values up to 0.3 in the steepest areas, with the predominance of litholic neosols. In these conditions of effective precipitation and runoff coefficients above zero, it is also possible to notice the occurrence of surface runoff, whose convergence in areas with high specific contribution, reaches flow peaks of up to 1.3 m 3 /s. Likewise, along Figure 6 it is possible to obtain a quantitative summary of the hydrological parameters previously mentioned. Also noteworthy are the occurrence of areas corresponding to the high runoff number classes; the tendency of low infiltration capacity of the soils; the consequent occurrence of effective precipitation; the occurrence of surface runoff above zero; and, mainly, the occurrence of peak flows of up to 1.3 m 3 /s, notably in the convergence segments of the slope flow, in the headwaters of the drainage channels.
In the rainy season and considering the soil in the humid condition, the maximum value of the runoff number went from 91 to 95.9, with the class from 80.1 to 90 becoming predominant in approximately 81.7% of the area. Under these conditions there is also a tendency for a generalized reduction in the infiltration potential, with a greater predominance of the class ranging from 25 to 60 mm. Even so, even under conditions of high antecedent humidity, a high infiltration potential is verified, especially along the plains of the main drainage channels. This finding corroborates the influence of the association of alluvial and hydromorphic soils, sometimes deep, developed along wide plains of low altimetric gradient and mostly covered by forests, such as gallery and riparian forests, with high infiltration potential.

Conclusions
The morphometric and pedological conditions related to soil cover, use, and management, suggest hydrographic systems with high hydrological potential, as well as considerable travel time, especially for the four higher-order channels. This potential is conditioned by the predominance of extensive areas of flat relief associated with soils with reasonable depth and texture varying from clayey to very clayey and with use predominantly by agriculture in level terrace systems.
These conditions act to keep the surface runoff at low velocities and, consequently, provide more time for concentration, which is essential to the infiltration process.
With reference to the evaluation of the set of hydrological variables, it is concluded that the area of the sub-basins on the right bank of the São Marcos River constitutes an area with high hydrological potential, even in humid soil conditions. This fact can be corroborated when one associates the elongated format of the higher-order basins with their long and wide fluvial plains that aggregate a high infiltration potential, even in the rainy season. If considered the response time of hydrological conditioning and its contribution to the increase in reservoir volume, this implies more time for the rainy season to contribute more effectively to the reservoir. However, this high infiltration capacity in most of the basins contributes to the postponement of the effects of droughts on the flow of drainage channels, as long as the existing dams along them release the minimum background flow.
This work is the beginning of a series of initiatives aimed at better understanding the relationship between the rainfall regime and the behavior of reservoir levels. In this sense, we highlight the need for advances in the methodological proposal, using field data for validation that are representative of the variability of hydrological conditions for the entire reservoir contribution watershed and in greater temporal detail. These adjustments are essential for the structuring of a Geographic Information System applied to the spatial-temporal monitoring of hydrological conditions throughout the basin and its relationship with the behavior of the reservoir level.