Hydrological modeling of a detention reservoir : flood control and aquifer recharge

In this paper, we proposed studying a small urban drainage catchment, named Pilot Basin of Mirassol (PBM) in Natal, RN Brazil whose outlet is a detention and infiltration reservoir (DIR). The rainfall-runoff transformation processes, water accumulation in DIR and the process of infiltration and percolation in the soil profile until the free aquifer were modeled; and, from rainfall event observations, water levels in DIR, free aquifer water level measurements, and also, parameter values determination, it was able to calibrate and modeling these combined processes. The mathematical modeling was carried out from the use of a distributed rainfall-runoff model, and besides, we developed a model to simulate the soil percolation in an unsaturated porous media. Continuous simulation was run over a period of eighteen months in time intervals of one minute. The generated hydrographs were transformed into inlet volumes to the DIR and then, it was carried out water balance in these time intervals, considering infiltration and percolation of water in the soil profile. As a result, we obtain an evaluation on the storage water process in DIR as well as the infiltration of water, redistribution into the soil and the groundwater aquifer recharge, in continuous temporal simulation. We found that the DIR has good performance in order to storage storm water (floods) and contributing to the local aquifer recharge process.


INTRODUCTION
Disorderly and intense occupation of urban areas without integrated planning of various infrastructures for sustainable development has increased drainage problems during heavy rain events.This fact is notorious especially in areas of periphery or suburban, mainly for lack of efficient planning usually observed in Brazilian urban catchments.
The occurrence of urban territorial expansion coupled with the lack of enforcement of following appropriate legislation and supervision to ensure proper disciplining of soil use and occupation is an inducer of flooding problems.
Drainage problem distributes along the natural lines of surface runoffs, which depends strongly on catchment topography and waterproofing degree (RIGHETTO et al., 2009a).
New conceptions for urban drainage regarding the use of hygienists systems and urban compensatory techniques to counter balance waterproofing by green areas and the use of lower topographic areas as local accumulation of drainage water are recently being discussed and applied.
Thus, the spread of the concept and application of low-impact techniques in providing urban drainage come gradually changing urban drainage systems, particularly when there is integration between urban infrastructure and city landscape adaptations.
In the city of Natal, several natural lakes such as New Lagoon (Lagoa Nova), Dry Lagoon (Lagoa Seca), Jacob Pond and Meadow Lagoon (Lagoa da Campina) disappeared.From the 1970s, there was a significant increase in occupancy areas in surrounding of these lakes, by construction of parks, schools, public and home buildings, inducing the transfer of the natural urban drainage system of the city to gallery networks.However, these practices intensified flooding episodes by increasing observed frequencies on road system.
In storm events, these occupied ponds resurface forming flooded areas because they are closed depressions where runoffs are directed as result of intensified waterproofing surface.The reduction of lake surfaces makes them unable to retain large runoff volumes even in more frequent lower rainy events; and thus, causing floods in their neighborhoods.
The Natal city drainage catchments practically are inserted in the municipal territory, because Natal receives small external surface water contribution to their urban area, which can be neglected in the global water balance of the drainage system.For water balance effect, there is a division between Natal drainage areas: open and closed catchments.
According to PDDMA (NATAL, 2009), closed drainage areas without runoff outlets are located in the western and southern areas of the city, including the catchments of Dunes Park, Jaguarari Lagoon, San Vale / Satellite City and Lagoinha, with total drainage area of 3,787.2ha, corresponding to 26.7% of the city territory.
Hydrogeological characteristics of most existing drainage areas in Natal -drainage basins closed with permeable soils, favoring the development of practical induction infiltration process.These drainage elements of preserving infiltration areas in lots and ponds of detention and infiltration reservoir (DIR) are very efficient to receive water from the public drainage system and promoting induced aquifer recharge.
The DIR are open structures that comprise the macro drainage urban system.The purpose of DIR is to minimize the impact of reducing natural storage capacity of the basin.In general, in closed drainage catchment, the aquifer recharge is enhanced by the drainage system.
With soil sealing, withholdings by detention, infiltration and decreased evapotranspiration, there is increased runoff and, consequently, the volumes that are directed to the DIR.Thus, the DIR also come to be an important instrument for aquifer recharge.
The infiltration is an important hydrological component responsible for the process of transferring surface water into the soil.The ability to infiltration depends on the soil characteristics and the state of moisture content of topsoil in the unsaturated zone (TUCCI, 2007).
By introducing LID in drainage systems, one has to consider the coupling between rainfall/runoff and localized storage/infiltration, taking advantage of depression and green areas to hold runoffs produced in relatively small plots that compose the drainage catchment.
With respect to DIR it is important to characterize the infiltration process that occurs in ponds of the macro drainage systems.This characterization is essential for reservoir dimension determination, its efficiency in flood control and in the estimation of local aquifer recharge.
The motivation for this work was to evaluate the operation of a DIR considering the ongoing process of rainfall-runoff events in the Pilot Basin of Mirassol (PBM), the accumulation of water in DIR, the infiltration, the percolation process and the local aquifer recharge in the Dunas-Barreiras groundwater, an important Natal's water resource.
This work is based on mathematical modeling of rainfall-runoff hydrological processes in the PBM and the infiltration, percolation and recharge in the DIR-M (The Detention and Infiltration Reservoir of Mirassol).
Data surveys were conducted, by monitoring rainfall events, water levels in the reservoir and the aquifer water table.Laboratory tests permitted estimating values for the main parameters of the soil used for modeling.
The integrated rainfall-runoff-storage-percolation-recharge modeling is the main contribution of this work for urban hydrology area.Because it was developed in an urban pilot watershed, the surveys and systematization of observations give solid foundation to study investigative work of the hydrological processes with possible inclusion of qualitative aspects and management in urban watersheds.
The developed mathematical model is sufficiently robust in terms of idealization -topology of the area, parameters considered and estimated from observations and the possibility to widespread use to other watersheds as it has been used among others, for example, by the SWMM application (JAMES; ROSSMAN; JAMES, 2010; LIMA; RIGHETTO, 2013).
Hydrological modeling of a detention reservoir: flood control and aquifer recharge Its idealization has targeted to achieve applicability in terms of ease of use and amount of information necessary to use the model.Enables introduce storage areas (LID) to reduce runoffs due to the urbanization process.It just need to define blocks or storage lots and allow runoffs street entrance to the interior of these areas.
As a contribution towards the process of storing, infiltration and recharge, the model requires parameters for the recharge process in one-dimensional domain despite a region typically three-dimensional; however, in this case, an adjusted parameter is proposed.
Besides the construction and testing of the integrated model, this study aims to understand the ongoing process of storing and drainage infiltration in a reservoir detention in order to analyze the process of percolation and timely recharge the groundwater aquifer as well as to evaluate DIR efficiency in control flood events.

RAINFALL-DISCHARGE MODEL
The distributed hydrologic model proposed by Righetto (2014) and applied to urban catchments is composed by water balance and resistance flow equations in a discretized domain formed by cells of the drainage area.The cells are blocks or street section units.
With the sequence of precipitation in time for each time interval (adopted the interval of 1 min) it is carried out water balance in each block, taking the precipitated volume, accumulated storage volume and the runoff to the street.Runoff coefficient values vary according to the accumulated water on the block.The larger the volume of water stored, the greater the runoff coefficient.Briefly, the following equations are used for the block modulus: where: q C : Block runoff coefficient; max C : Maximum runoff coefficient; min C : Minimum runoff coefficient; 1 a : Adjusted parameter; arm H : Water height in the block, m; q Q : Block runoff to street, m 3 /min; q A : Block area, m 2 ; q α and q β : Parameters of discharge equation (Manning equation).
The mathematical module for streets consists of taking the contributing flows from adjacent blocks and upstream streets directly linked to a particular street section.In this section, it is carried out water balance in time intervals of 1 minute, taking the direct precipitation, accumulation of water in the gutters and the flow discharged by this street section.The following equations summarize this module: where: m Q : Discharge contribution from adjacent upstream street, m 3 /min; r C : Street tunoff coefficient; P: Rainfall at time interval, mm; r A : Hydraulic area of street section; m 2 ; r Q : Outlet street discharge, m 3 /min; r H : Depth of water flow in the street, mm; : Discharge parameter of Manning equations, where I r is the street slope (m/m), n the Manning coefficient and .α'r geometric parameter of hydraulic street section.
As pointed out Righetto (1998), analytical solutions of Richards's equation are restricted to simple problems.On the other hand, several numerical techniques are available for obtaining numerical solution related to the equation of infiltration, percolation and redistribution of water in the soil.Among the numerical techniques, the present model incorporates finite differences precedure.
Figure 1 shows the discretization domain containing N cells and N 1 + nodes, so that each cell i has equal length, denoted z ∆ .
Figure 1.Discretization domain of soil profile.

Righetto and Freitas
The domain was limited at node N that was located at 3000 cm (30.0 meters) deep from the pond bottom.It was adopted a constant value for Δz, equal to 2.0 cm; thus, N = 1500.This cell size was adequate to ensure stability of the dynamic modeling of water percolation in unsaturated porous profile.
The discretization terms of the Richards equation using the finite difference approach, results in the following algebraic equation: where: (10) j : Time index and i: Space index.Substituting Equations 8 and 9 into 7 results in an expression that permits the determination of Ψ + through the solution of a tridiagonal linear system (RIGHETTO, 1998): For internal nodes, 1 i N 1 ≤ ≤ − , the parameters , A , e i B are expressed by: ( ) Boundary conditions may be of Dirichlet type (specified pressure or suction head) or of Neumann type (specified flux) or of Cauchy type, (flux is a function of the phreatic level position).
The following conditions were adopted in the simulations of this study: a) Lagoon basement (z 0 = ): It was considered that ( ) ) and depth of phreatic level above 23.00 meters: 1, 1.00 For this condition, the exploitation discharge value, q, is constant and equal to 0,001 cm/min, estimated by phreatic water level measurements.c) Base condition 2 ( máx z z = ) and depth of phreatic water level less than 23.00 m: where: β : Adjusted parameter for lateral percolation; s K : Saturated hydraulic conductivity, cm/min; ZLF: Phreatic water level depth, cm.
The base boundary condition 2 was implemented because model tests show that base boundary condition 1 (constant exploitation flow) after a certain volume of infiltrated water was not enough for decreasing water table and then, allowing soil profile to become completely saturated.
However, this situation is not observable.It is known that the local aquifer recharge forms a saturated cone-shaped profile and the horizontal lateral flow intensify with decreasing the depth of the regional water table.So it was introduced a law to estimate the discharge of exploitation or flow redistribution side, taking a linear growth as the groundwater level rise relative to regional water level, estimated as equal to 23,00 m.For this case, the flow redistribution was defined by the expression: where: 1 β is an adjusted parameter; s K is the saturated hydraulic conductivity (cm/min); ZLF phreatic water table depth (cm.) System of equations in finite differences together initial and boundary conditions is solved to determine the pressure or suction head distribution, which permits to determine the moisture distribution as a function of time and thus, allowing to simulate Hydrological modeling of a detention reservoir: flood control and aquifer recharge the movement of soil water.No special procedure is required in the saturated soil zone, because in this case, C (Ψ) becomes equal to zero and the Richards model reduces to a Laplacian.For each time step, the system of equations to be solved consists of a tridiagonal linear system of equations; in this work, solved by Thomas algorithm elimination method.

RESERVOIR WATER BALANCE
Matched to the unidimensional flow model in unsaturated porous media (percolation model) there is the model for water balance of water stored in the reservoir.Processing rainfall-discharge model one obtains water inlet volume to the IDR, t a V along time (eighteen months at time intervals of one minute).
Dividing this amount by the reservoir water surface area at time t, the increase in height, t h ∆ , due to this affluent volume is obtained.
For each time, t, it is used the value of the hydraulic conductivity, t 1 K and the heads t 0 Ψ , t 1 Ψ referring respectively to the surface of the discretization node (lake bottom) and the adjacent node of the first layer or soil profile cell.The infiltration is determined by Darcy's equation: where: t F is the infiltration at time t (cm.min -1) .
After calculation of infiltration, water balance in DIR is obtained by the following water balance equation:

STUDY AREA
The Mirassol pilot basin (PBM) is located near Federal University of Rio Grande do Norte and has been studied since 2003 (RIGHETTO et al., 2009b).
Figure 2 illustrates the drainage area of Mirassol pilot basin, with the location of the DIR, object of this study.The pilot basin is fully urbanized, occupied by residential of middle class and located in Capim Macio neighborhood, south of Natal, RN. Figure 3 shows the watershed delineation and it is topography.
PBM relief is representative of the region with predominant lacustrine formation associated with the presence of dunes.The layer of topsoil, frankly sandy, has high infiltration rates.The soil horizon B is characterized by a lower permeability profile called Barreiras and is an important aquifer water supply of the city.As a whole, the region has high water storage capacity due to the presence of Dunas-Barreiras aquifer.The basin area is 0.14 square kilometers with a maximum elevation of 61.10 meters and altitude at the downstream end of 34.90 meters.
Runoff is directed to a DIR (Mirassol lagoon), which is situated in the region with topographic elevation lower than other areas of the catchment.
There is only one small length drainage collectors, located a few dozen meters upstream this reservoir.The entrance of flow

Righetto and Freitas
is through open small channel connect to a 1.00 m diameter pipe (RIGHETTO et al., 2009b).
Above maximum height, surplus water is transferred through a gallery section connected to another reservoir located downstream of the studied lagoon.
As foregoing mentioned, the standard occupancy of MPR is single residential -family, class type B and C, with predominance of 200 m 2 lot area (LIMA, 2011).It is noteworthy that those single residential buildings densely occupy the entire area.

Blocks and streets runoff coefficients
Runoff coefficient values used in rainfall-discharge modeling were initially estimated as recommended by the Manual of Drainage Plan Director of Natal, as indicated in Table 1.
For street sections with asphalt, paving it adopted a runoff coefficient of 0.90 was adopted.For the other street sections of the basin, the pavement is covered with cobblestones and, then, a runoff coefficient of 0.50 was considered.For Block 1 (considered atypical) it was performed a weighted average of various surfaces that have been identified (roofs, green areas and permeable paving) resulting in a runoff coefficient of 0.27.
For the other blocks, it was preliminarily considered values of runoff coefficients suggested by PDDMA (NATAL, 2009) for high percentage waterproofed lots.
By default occupancy and aerial photos analysis, it was considered that the blocks are 85% waterproofed resulting in a runoff coefficient of 0.71.These values were adopted preliminary as runoff coefficient (C) and after altered within a range (Cmin and Cmax) based on calibration simulations of PBM rainfalldischarge modeling.

Precipitation data
The rainfall data were obtained from an automatic rainfall station located within a public area -the State School Floriano Cavalcanti (Floca), inserted in the drainage basin of this study.The observation period began in April 2014 and ended in September 2015, with storage resolution of one minute.

PBM geometry
PBM has been schematically discretized (spatially) as illustrated in Figure 4 in order to define the topology and numbering blocks and street sections.

Soil physical properties
Firstly, it was performed granulometric test with three soil samples taken from the reservoir bottom at depths soil profile of 2, 4 and 6 meters.The tests permit to conclude that the three samples have similar particle size distributions.Thus, the measurements indicated that the average sand, silt and clay are approximately 84%, 11.5% and 4.5%, respectively.With these average values and making use of the triangle textural classification, this soil is defined as free sand.
The soil textural classification was used in the RETC freeware, available on its developer site, the PC-PROGRESS ENGINEERING SOFTWARE DEVELOPER.According to the website of PC-PROGRESS, this software is used to analyze retention functions of water and hydraulic conductivity for unsaturated soils.This freeware makes use of Rosetta database (SCHAAP; LEIJ; VAN GENUCHTEN, 2001), which allows predictions of hydraulic soil parameters from textural classification.
It was observed during numerical tests of the percolation model that infiltration was above the values observed in the field.In order to bring infiltration near the observed field values, s θ e s K (provided by RETC) were adjusted to values of 0.40 cm 3 / cm 3 and 0.18 cm / min, respectively.
The equations that define the moisture content in the soil (θ (Ψ)), the hydraulic conductivity (K (θ)) and the specific hydraulic capacity (C (θ)), based on the theory of van Genuchten (1980) and Mualem (1976) model are as follows: where: Ψ is the suction or pressure head (cm); ( ) K θ : Hydraulic conductivity as function of soil moisture content (θ ), cm/min; The graphics that describe the behavior of these equations for the DIR-Mirassol are shown in Figure 5.
In this study, it was essential to consider in the surroundings of the study area existing groundwater exploitation carried out by the local urban water supply service (CAERN).
Groundwater level observations were carried out in the period from 01/10/2015 to 23/11/2015 (dry season) in order to determine weekly the aquifer exploitation flow (CARDOSO, 2015).On average, it was observed that the phreatic level varies approximately -17.0 cm per week.With this value, the exploitation flow, q, was estimated from the test flow percolation model.
The value of q resulting in a weekly decay of water table of -17.0 cm was 0.001 cm / min.This value was adopted as the basis of flow exploitation to the hydrogeological study of the DIR.

PROCEDURES
The study includes monitoring precipitation, conducted in the automatic station installed in PBM as well as water level measurements in the reservoir and in phreatic aquifer.
The data were stored every minute (compatible with the discretization in time used in modeling) and monthly collected during the period from April 2014 to September 2015.
Much of this study was devoted to computer modeling process.The modeling was performed using the Microsoft Developer Studio program using the Fortran 90 programming language.The rainfall-discharge model developed by Righetto (2014) underwent changes in its input data and calibration in Hydrological modeling of a detention reservoir: flood control and aquifer recharge relation to the expression that defines the coefficient a runoff variation over time of occurrence of rain event.

Modeling
With the rainfall data at intervals of 1 min the modelling was processed in order to determine the discharge flow (m 3 /min) at any street section and then, the temporal discharge to the DIR.
The inflows to DIR, in turn, are used as input data for the one-dimensional modeling the soil profile.For each flow generated by rainfall-runoff model (including zero), it is carried out water balance of the reservoir.Coupled to the water balance model, the infiltration model and percolation of water into the unsaturated porous medium is simulated through which is determined temporal position of phreatic water level at each time interval of 1 min.

Global simulation
The evaluation of the models presented here, as well as verification of the DIR performance and estimates of groundwater recharge were done by analysis of the overall simulation over a period of 18 months (April/2014 to September/2015), with observations of rainfall and water levels in the lagoon and in the surface of the phreatic groundwater.
To overall simulation, it is set an initial condition for percolation model in unsaturated soil; so that it was assumed that the initial soil subsurface just below the reservoir is in the field capacity and, therefore, with a small reservoir bottom flow to the water table.
As initial condition (t = 0), it was considered that the suction head, Ψ is in field capacity, Ψcc equal to -50 cm that happens in the reservoir bottom and the surface of the groundwater phreatic level defined by suction head equal zero.Then where: LF z is the depth of the water table.The z coordinate has value equal zero in the reservoir bottom with gravity direction.
This situation defines a constant flux at field capacity; then, for i from 0 to N LF-1 , ( ) For i varying from N LF +1 until NI (last domain node) the hydrostatic distribution was considered, i.e.: Over time, t≥0, it was considered and Lagoa H , obtained by simulating rainfall-discharge processing model in the watershed were stored to be used in the simulation of DIR water balance and soil percolation-groundwater recharge.
The last layer in the profile is considered waterproof, i.e., zero flow.The flow of exploitation was defined according to the phreatic surface level and adjusted according to values observed in a well installed beside DIR Figures 6 and 7 provide the main results of the global simulation for the period of April to September 2014, in which the largest annual rainfall event occurred.
Each graphic was constructed in daily scale with corresponding value variables results.Precipitation, P (mm); Volume affluent, V afluente (m 3 / day); maximum daily water level in the reservoir, H Lagoa Máx.(cm.); average level of the water table, ZLF Médio (x 10 -1 cm).
The lines corresponding to P, V afluente , H Lagoa Máx. and ZLF Médio vary with time due to rainy events that occurred during the simulated period and also due to the soil saturated degree between the reservoir bottom and the phreatic surface..
In Figure 6 it is important to note that there was a relative delay in recharging process due to the fact that virtually all the Righetto and Freitas soil column between the pond bottom and the water table was relatively dry (θ≅0,10), a fact that requires a longer time to forward saturation progressed and reached the phreatic water.Therefore, within the first three months of simulation, there is the beginning of the recharge period takes place in mid-June.
Also in this figure, one should note the rainfall events that occurred between the days 06/13/2014 and 15/06/2014.These events were marked in Natal memory due to the fact of the occurrence of several flooded areas in the city and have caused serious landslide occurred on Mãe Luiza neighborhood just at the time it was being held the World Cup Soccer.For this reason, this exceptional storm in Natal is called by rain 2014 Cup.
The rain gauge installed in BPM recorded three rainfall events during this period, which amounted to 343.7 mm with a cumulative time of 38 hours, estimated in total as 100 years return period event.
Finally, due to the fact of the occurrence of an event of that magnitude associated with the presence of previous rainfall events contributing to the increase in the degree of saturation of the soil, it is observed that recharging began to occur 5 days after the beginning of the rainfall event.Thus, this time period between infiltration and groundwater recharge occurrence can be considered small compared to the mean time of about 20 days.
In Figure 6, it is also observed that since beginning of July until early August has increased occurrence of rainfall events.
In this period, the presence of rain increases the saturation degree in soil profile, because in short time intervals occurred accumulation of water of different events, unlike the dry period; for this reason, the aquifer recharges are more frequent and intense.
The rest of the period (Figure 7), i.e., mid-August and September (except precipitation at 08 and 09 September) it is observed the beginning of the dry season, decreasing both the frequency and the magnitude of rainy events.
In the remaining months, the simulation confirms the permanence of the dry season with lower magnitudes of rainy events throughout the observation period.However, due to the fact of the occurrence of large rainfall events and recharges that preceded this period, coupled with the presence of a few rainfall events, although of low magnitude, recharges were sufficient to prevent excessive decay of level of the phreatic water table.During this period, it is assumed that the equilibrium of the groundwater surface was maintained more by moisture redistribution than the presence of rain itself.
The estimated level in the DIR was 1.84 m and occurred on April 11, 2014.Although it was an event with a total of 56.4 mm and duration of 150 minutes, it resulted in higher reservoir water level due to the highest peak inlet discharge observed along all the simulated period.
In all simulation, the phreatic water level as a function of time is strongly influenced by the exploitation flow rate.That is, the lateral flow happens from the aquifer soil profile as a whole.Through the parameter β, the model shows greater or lesser intensity of aquifer recharge, so that the groundwater may answer more or less depending on the recharge.
Because the study focus on unidimensional recharging process, the simulations showed that, the β parameter should assume values above 0.25.Values below this would lead the entire soil profile become saturated during intense rainfall events and long duration.It is possible, but not observed in the time period studied.
Moreover, β values greater than 0.75 would decrease water level excessively far superior to values of the water level observed through the observation well.

Local simulation
To better explain the behavior of the percolation of water in the soil profile, it was highlighted a period of strong recharge between day 11/04 and 21/04 2014.
Figure 8 shows the rainfall distribution and the corresponding hydrograph generated by the rainfall-discharge model.This event Figure 9 shows the evolution of the moisture profile into the soil, the front of wetness over a period of 10 consecutive days from the start of the main event (04.11.2014 -02h27min).
At time t = 0 h, it is observable the wetting front of the main event located practically in the reservoir bottom (Z = 0.0 m) as well as the existence of moisture front located 4.0 m below (from earlier rainfalls, which occurred on 8, 9 and 11 April).
Time t = 16 h corresponds to the instant at which the reservoir empties completely.It also shows the advance of the wetting front of the main event.It is noted that the saturated soil layer reaches a thickness of approximately 6.5 m.
At time t = 24 h, in the reservoir bottom soil moisture returns to field capacity (upper boundary condition when the lagoon is dry) and there is the advance of the main wetness with narrowing gap between the existing fronts.
At time t = 48 h, it is noted the appearance of a wetting front above the main front due to occurrence of 2.30 mm and 5 minutes duration rainfall, on April 12, 2015.
At time t = 120 h, the main front wetness merges with the front of initial time, i.e., t = 0h, besides highlighting the slow advance of the front from the day of the event (Abril 12, 2014).
In addition, in the time t = 240h, corresponds to the moment that begins the phreatic aquifer recharge

Rainfall design simulation
In order to assess the performance of the infiltration process that occurs in the reservoir, it was simulated an extreme storm in the PBM and thus, getting the maximum water level in the reservoir and its efficiency evaluated.It was adopted in this simulation a triangular rainfall distribution so that the maximum intensity positioned at half rain duration time.
The hyetograph refers to the intensity duration of 5 days and return period T = 50 years.The rainfall-discharge modeling was run for obtaining the hydrograph corresponding to 50 years  Righetto and Freitas return period along 120 h.For this return period and 5 rainy days, the total rainfall is equal to 67mm.With these values, the reservoir water balance was simulated, considering the process of infiltration and soil percolation.By means of coupled model, the models are processed concurrently within a single calculus environment.
Figure 10 shows the time dependent water levels in the DIR-Mirassol resulting from storm rainfall i(tc, T) or i(120h, 50anos).
During this period, the maximum level can lead to flooding around the reservoir if the water transfer mechanisms are inadequate, poor or damaged.
The proper characteristic of Natal ponds in flood control is the topographic elevation of about 30 m above the sea level.This elevation avoid most lagoons do not become contours of the potentiometric surface of the Dunas-Barreiras aquifer, preserving the high infiltration capacity of LDI.
V : Total volume infiltrated by field capacity condition; S : Water balance.
This error can be attributed to the temporal discretization process and by restricting water content to values larger than field capacity, which maintain residual flow even in the condition of dry reservoir.

CONCLUSION
The integrated model (rainfall-runoff-storage-infiltrationrecharge) developed for PBM proved to be effective and easy to use as tool for the analysis of drainage systems involving hydraulic infrastructure and LID installation.
The hydrological simulations over a sequence of 1 min time intervals and eighteen months duration showed adequate to evaluate flood control of Mirassol reservoir as well as infiltration and percolation of water to phreatic aquifer.
Two main aspects of the reservoir worth mentioning in relation to damping storage water despite its relatively small volume capacity: high hydraulic conductivity of the soil and position of the phreatic water table.
The infiltration capacity is high because the lagoon is located over dune formation and has adequate topography once the water table depth does not reach values that might compromise infiltration and reduction of water storaged in the reservoir.
However, in situations of high intensity and persistent rainy days, it is expected that all the region could contribute to the recharge and, in this situation, the lagoon can become ineffective and vulnerable to flooding.This fact occurs in several lagoons in the city of Natal, the ones located in lower elevation relative to sea level.
In the construction of percolation model, the parameter β was introduced to consider the lateral dispersion of the infiltrated water.This parameter is highly sensitive to lower or increase the variation of the water table and thus the infiltration process and recharging.
It is important to emphasize that the design and detention reservoir operation should consider the existence of groundwater and its influence on the infiltration process.
It is noteworthy that the pond emptying time for persistent rainfall events significantly alter its infiltration capacity.With the presented model, it is possible to evaluate the operation of DIR for the case of intense and successive rainy days.
Regarding the rainfall-runoff model, the natural complexity of the drainage basin was considered and relatively simplified.The division of drainage area into blocks and streets and, from

2a:
Adjusted geometric parameter; ( ) r f H : Convertion function of depth of water in the street to hydraulic depth taking into account pavement, gutter and sidewalk;
50 m, which corresponds to the volume of water stored equal to 4,450 m 3 .This level defines the reservoir overflow outlet, which transfers the surplus to another DIR.
Saturated and residual soil water content, values equal to 0.40 e 0.057, respectively; s K : Saturated hydraulic conductivity, value equal to 0.18 cm/min; , ne m α : Soil parameters, equals to 0.124, 2.28 and 0.56, respectively.

Figure 9 .
Figure 9. Evolution of soil moisture profile.
It was supposed that a very low values in relation to field capacity occurred in a very thin soil superficial layer and is not relevant for the determination of soil infiltration.
Ψ value can assume field capacity cc Ψ value (lagoon is dry) or t Lagoa H value (lagoon with storage water).máx z z = Incremental height at time interval between t e t+Δt, (cm).
H: Water level in reservoir at time t, (cm); t h ∆ :