## Servicios Personalizados

## Revista

## Articulo

## Indicadores

- Citado por SciELO
- Accesos

## Links relacionados

- Citado por Google
- Similares en SciELO
- Similares en Google

## Compartir

## Obras y proyectos

##
*versión On-line* ISSN 0718-2813

### Obras y Proyectos no.13 Concepción 2013

#### http://dx.doi.org/10.4067/S0718-28132013000100003

**On the proceeding of a saturation front under ponded conditions**

**Sobre el avance de un frente de saturación bajo condiciones de apozamiento**

**Benjamin Sinaba ^{1}, Bernhard Becker^{2}, Wiebke Klauder^{3}, Iván Salazar^{4} and Holger Schüttrumpf^{1}**

^{1} Lehrstuhl und Institut fur Wasserbau und Wasserwirtschaft, RWTH Aachen University, Mies-van-der-Rohe-Strasse 152056 Aachen, Germany, sinaba@iww.rwth-aachen.de, schuettrumpf@iww.rwth-aachen.de

^{2} Deltares, Rotterdamseweg 185 2629 HD Delft, P.O. Box 177, 2600 MH Delft, The Netherlands, bernhard.becker@deltares.nl

^{3} RWE Power AG, Zum Gut Bohlendorf, 50935 Bergheim, Germany, wiebke.klauder@rwe.com

^{4} Departamento de Ingeniería Civil, Universidad Católica del Norte, Avda. Angamos 0610, Antofagasta, Chile, isalazar@ucn.cl

*The vertical infiltration of water is a variable saturated flow process until the soil under an inundated area is completely saturated. However, in several modelling studies that regard the river-aquifer interactions, the vertical infiltration under ponded conditions is modelled as a saturated soil flow process. In order to find out whether the assumption of saturated conditions under an inundation area is suitable or not, the saturated time of the soil below an inundated area is determined. This is realised by evaluating an analytical infiltration model with two example soils under different initial conditions and for different ponding depths. The model results show, that saturation of several meters of sandy soil under ponded conditions takes place within a period of hours. For the clay example, the saturation process happens within the order of days, and if the clay is initially dry, it can take even weeks. This leads to the conclusion that saturated model approaches are suitable for sand and gravel aquifers if they are not covered with a clay layer. While the ponding depth influences the infiltration rate and the progression of the saturation front in the sandy soil, for initially wet clay the infiltration rates were found to change only marginally with the ponding depth.*

* Keywords: infiltration, ponded conditions, inundation, vadose zone, Richards's equation, Haverkamp infiltration model*.

*La infiltración vertical de agua es un proceso de flujo parcialmente saturado hasta que el suelo ubicado bajo un área inundada esté completamente saturado. Sin embargo, en muchos estudios de modelación sobre interacción de agua superficial-subterránea, el proceso de infiltración vertical bajo condiciones de apozamiento es modelado como un proceso de flujo en suelo saturado. Con la finalidad de estudiar si esta simplificación es adecuada, se determina el tiempo de saturación de un suelo ubicado bajo un área inundada. Esto es realizado mediante la evaluación de un modelo de análisis de infiltración a través de dos ejemplos de suelo bajo diferentes condiciones iniciales y para diferentes profundidades de apozamiento. Los resultados del modelo indican que para un suelo arenoso bajo condiciones de apozamiento, la saturación de varios metros se lleva a cabo en horas. Para el caso de suelo arcilloso el proceso de saturación ocurre en días, y si la arcilla se encuentra inicialmente seca puede durar incluso semanas. Esto lleva a concluir que la simplificación de un modelo saturado es adecuado para suelos arenosos y gravosos, siempre y cuando no sean subyacentes a un estrato de arcilla. En suelos arenosos la profundidad de apozamiento influye en la velocidad de infiltración y en el avance del frente de saturación, mientras que en suelos arcillosos inicialmente húmedos la profundidad de apozamiento influye marginalmente en la velocidad de infiltración.*

* Palabras clave: infiltración, condiciones de apozamiento, inundación, zona vadosa, ecuación de Richards, modelo de infiltración de Haverkamp*.

**Introduction**

In recent years, floods have been causing significant damages in urban areas of Germany. One of the most disastrous flood events occurred in August 2002 in the Elbe river catchment. Rising groundwater due to high water levels of neighbouring rivers (subsurface flood) became a significant problem for the City of Dresden (see SLUG, 2003; Huber *et al,* 2003; Sommer and Ullrich, 2004; Kreibich and Thieken, 2008). In this case, not only bank storage effects (Freeze, 1972) were responsible for water inflow from the river into the adjacent aquifers, but also vertical infiltrations from inundated areas contributed considerably to the groundwater head rise (SLUG, 2003).

Predictive numerical simulations are an appropriate way to support regional planning or risk assessment with respect to groundwater head rising in alluvial aquifers due to high water in neighbouring rivers (Becker, 2010; Becker *et al.,* 2012). Until the complete saturation of the soil below the inundated area the vertical infiltration under ponded conditions is a variable saturated flow process that can be represented by the Richards' equation (see Freeze and Cherry, 1979). However, in several studies the infiltration from an inundated area into the aquifer is assumed to be saturated by modelling the infiltration with a first order head boundary condition or third order leakage boundary condition in a saturated groundwater flow model *(e.g.* Homann, 2006; Monninkhoff and Kernbach, 2006; van Linn, 2006; Bauer *et al.,* 2006; Koch *et al.,* 2007; Peetz and Steckel, 2007; Sommer and Eulitz, 2007). For practical applications such a saturated approach is often chosen for the following reasons: it is computationally less expensive to assume saturated flow only, there is a lack of program features to handle infiltration processes in the used groundwater flow simulation program and in practice the additional parameters necessary for variable saturated flow are often not available for the modelling area.

This article addresses the question of how fast the saturation front under an inundated area proceeds and how long it takes until the soil below is completely saturated. A one-dimensional analytical infiltration model is evaluated for two soil examples under different initial conditions and different values of ponding depth. For hydrological studies and studies within the context of irrigation usually the infiltration rate and the cumulative infiltration are the parameters of major interest. We present diagrams with the progression of the saturation front. Within the frame of modelling subsurface floods these diagrams help to decide whether the simplification of the process of infiltration under ponded condition as saturated flow is feasible or a more advanced approach like the Richards model is necessary.

Initially, the mathematical background of the infiltration model used is introduced. Afterwards, the parameters of the example soil for which the model is evaluated are given. Subsequently, the model results are discussed. The article closes with a summary and conclusions.

**Mathematical infiltration model**

Haverkamp *etal.* (1990) present a mathematical infiltration model for infiltration under ponded conditions. The fact that it is based on physical soil parameters makes this model appropriate for this study, because parameters of example soils from the literature can be applied. Haverkamp *et al.* (1990) derive the cumulative infiltration I as a function of time t based on the one-dimensional Richards' equation:

where h is the ponding depth, *i.e.* the water level above the top ground surface, being assumed as constant over time. θ_{i} is the initial moisture content and θ_{s} denotes the soil water content at saturation. The symbol Ki stands for the initial hydraulic conductivity under unsaturated conditions and K_{s} is the hydraulic conductivity at saturation. The sorptivity S describes the impact of the matrix potential ψ on the infiltration of water into the soil. It embodies the influence of the capillarity capacity of dry soil in one single parameter, whereby S takes the SI-unit ms^{-0.5} (see Philip, 1957). Strictly, the sorptivity has meaning only in relation to an initial state of the porous medium and an imposed boundary condition (Philip, 1969). According to Freeze and Cherry (1979), ψ _{a} is denoted as the air entry pressure head, which is a constant soil parameter. For fine grained soils, the relationship between saturated moisture content and negative pressure head holds over a slight range ψ > ψ _{a} where the soil is practically saturated, but the pressure head is still negative.

Inserting the formulation of the time dependent infiltration rate q(t)

into equation (1) leads to (Barry *et al.,* 1995):

For this formulation, limits of application are given by the assumption of homogeneous isotropic soil. So layered aquifers can not be modeled and the model infiltration formula does not take into account macropore flow (see Beven and Germann, 1982). Furthermore, the horizontal movement of the soil water is not considered due to the one-dimensional approach. Finally, the model does not allow the groundwater table capillary fringe at the lower end of the model soil column as it assumes homogeneous initial conditions θ_{i} (Gillham, 1984).

A formula for the position of a saturated front z(t) is derived from the Laplace differential equation by Schmitz *et al.* (1989):

The saturation front is preceded by a partly saturated transition zone (Haverkamp *et al.,* 1990; Schmitz *et al.,* 1989). So the position of the saturation front is defined as the intersection between the nearly saturated zone and the partially saturated zone (Schmitz *et al.,* 1989). The space coordinate z is oriented in vertical direction of infiltration.

**Soil parameters of the model**

Chong and Green (1983) introduce an equation for sorptivity S(B) based on a combination of the Philip's equation (Philip, 1957) and the Green-Ampt's equation (Green and Ampt, 1911).

For its estimation the saturated hydraulic conductivity K_{s}, the moisture content at saturation θ_{s }, the initial water content Bi and the corresponding suction head ψ(θ_{i}) have to be known. For K* in equation (5) Chong and Green (1983) use the Green-Ampt-infiltration model parameter K_{s}. But the Green-Ampt-model assumes a piston shaped proceeding of the infiltration water, so the parameter does not account for a transition zone. Therefore, Rawls and Brakensiek (1983) suggest to use

to account for the conductivity conditions in the partly saturated transition zone.

The relation between water content B and the suction head ψ is given by the van Genuchten approach (van Genuchten, 1978),

where θ_{r} is the residual water content, α and n are fitting parameters. The relation between unsaturated hydraulic conductivity K to the corresponding water content θ is given by Schaap and van Genuchten (2006):

The empirical parameter L regards the tortuosity of the soil. S_{e} is the relative saturation to a corresponding water content θ(ψ). In van Genuchten (1978), the relative saturation S_{e} is noted as:

To identify the air entry pressure ψ _{a} , a water content has to be selected, whose value is near the water content θ_{s} at saturation. This is necessary, because the van Genuchten model disregards a negative pressure head at complete saturation. According to the van Genuchten model the total soil saturation matches to zero suction head.

With the equations (4) to (7), the evaluation of the infiltration model with different initial soil moisture content θ_{i} is possible. The implementation of the model into a computer code is facilitated by the dimensionless representation of the model as a function of the infiltration parameter γ (Haverkamp *et al.,* 1990):

The relationship between infiltration time and dimensionless infiltration time is given by

Solving equation (11) allows the computation of the dimensionless cumulative infiltration I* by

where q_{s}* is the dimensionless infiltration rate. The dimensionless cumulative infiltration and the cumulative infiltration are related with

and the implicit non-dimensional formulation of the infiltration rate is

Finally, the infiltration rate can be obtained from the dimensionless infiltration rate q_{s}* with

To model the infiltration process, the main working step is the solution of equation (11) in order to obtain the infiltration rate for given discrete time steps. Equation (11) is a combination of fractured rational functions and an elementary transcendental function and mathematically defined in a domain of

However, an unambiguous and physically feasible solution is given by q_{s}* > 1. The bisection method has been chosen to solve these equations, because a solution interval can easily be set (see Sinaba, 2007).

The implementation of the equations into a computer code written in Matlab is called Haverkamp-infiltration-model in this paper. Input values are soil properties, points in time, and initial conditions. Model results are infiltration rate, cumulative infiltration, and the progression of the saturation front.

**Example soils**

The two example soils for this study correspond to the textural classes of the USDA (2009). A sand example soil is given with a composition of 95 % sand, 2 % silt, and 3 % clay. The second example is a clay soil, which consists of 53.2 % clay, 31.7 % silt, and 15.1 % sand. Table 1 summarizes the soil properties and initial conditions for the example soils. Values for the parameters θ_{i}, θ_{s}, α, n, K_{s} and L were obtained from Schaap (2009). All other parameters were obtained from the equations listed above. Initial moisture contents were chosen out of the range of field moisture contents given by Kolymbas (1998), 0.05 < B < 0.1 for sand and 0.2 < B < 0.3 for clay.

Table1: Soil parameters of the example sand and clay for initially dry and wet conditions, respectively (Sinaba, 2007)

**Results**

**Infiltration rates**

Figure 1 shows infiltration rates over time, which have been computed with the Haverkamp infiltration model for the sand example. Curves for wet and dry initial moisture contents are given for ponding depths between 0 and 2 m. As expected, infiltration rates for initially dry sand are generally higher than for wet sand, because in the beginning of the infiltration process wet sands have a lower soil pressure head ψ. This leads to smaller sorptivity values which reduces the infiltration rate.

Figure 1: Infiltration rates q with ponding depths of h = 0.0, 0.5,

1.0, 1.5, and 2.0 m and saturated hydraulic conductivity K_{s} for

the example sand initially a) dry and b) wet

Increasing the ponding depth increases the infiltration rates, because high ponding depths cause high pressure gradients at the soil surface, although the infiltrated water reduces the matrix potential. However, the higher the ponding depth, the lower is the influence of increasing the ponding depths on the infiltration rate for both wet and drv sand.

The infiltration rate with a ponding depth of 0.0 m represents the case of a precipitation event where the precipitation intensity is equal to the infiltration capacity q. In this case, the infiltration rate converges to a value which is equal to the saturated conductivity K_{s}.

Different to sand, the infiltration rates for initially wet example clay only change marginally with varying ponding depths and hense they can barely be distinguished (Sinaba, 2007), for this reason the infiltration curves are not shown here. For this reason the suction head ψ of the clay is much higher than in the sand, computed infiltration rates for the clay is lower, because a higher initial moisture content Bi and a lower value of hydraulic conductivity K_{s} compensate the suction head value, resulting in sorptivity values almost equal to the sand example. Furthermore, clay has a lower air entry pressure ψ_{a}.

**Progression of the saturation front**

Figure 2 shows the progression of the saturation front for initially dry and wet sand examples under different ponding depths. The dashed line denotes a saturated depth of one meter.

Figure 2: Position of the saturation front z with ponding depths

of 0.0, 0.5, 1.0, 1.5, and 2.0 m for the a) dry and b) wet example sand

Without a water layer at the soil surface, the saturation curve is nearly linear and the infiltration period to saturate one meter of the soil from above the surface, accounts about 11000 s (3 h) if the soil is initially dry and 6000 s (1.7 h) if it is already wet. Increasing the ponding depth increases the velocity of the saturation front for both the wet and the dry sand. However, the higher the ponding depth, the lower is the impact of the initial moisture content on the progression of the saturation front. Also the impact of increasing the ponding depth on the progression of the saturation front reduces for large ponding depths.

While for the example clay the influence of varying ponding depths on the infiltration rate is only small, the ponding depth changes heavily the velocity of the saturation front progression for the wet example clay (Figure 3). Since the hydraulic conductivity depends on the moisture content of the soil, the wet clay has a higher hydraulic conductivity than the dry clay. The saturation front proceeds faster in the wet clay and the impact of the hydraulic gradient on the infiltration process and the proceeding of the saturation front is larger. With high ponding depths, saturation of an initially wet soil column of one meter takes time spans in the order of days, while the saturation of a one meter column of the initially dry example clay needs about one week.

Figure 3: Position of the saturation front z with ponding depths

of 0.0, 0.5, 1.0, 1.5, and 2.0 m for the a) dry and b) wet example clay

**Transition zone**

Neglecting the transition zone that precedes the saturation front, a water amount of

would saturate a soil column with a length of z. Under this assumption, the infiltrating water fills the available pore space completely with a sharp, piston-shaped wetting front like assumed in the Green-Ampt-model. The comparison of the infiltration time after a cumulative infiltration of I = (θ_{s} -θ_{i} ). 1 m has been achieved with the time after the saturation front has reached a depth of one meter allows conclusions on the shape of the transition front. Table 2 shows the results of this comparison: for high ponding depths the dry and wet example sand values only differ slightly (Table 2 only shows the data for wet sand). So a nearly piston-shaped wetting front of the saturation front can be concluded for the sand soil. Because the two time values differ significantly for dry clay, a large fraction of the infiltrated water is stored in the transition front that precedes the saturation front in the clay soil. However, displaying this transition front is not possible with the used infiltration formulas, because the equations do not provide any information on the moisture distribution.

Table 2: fnfiltration time t in seconds after the cumulative

infiltration has reached the value of I = (0_{S} - Θ, )-Im and the

time after the saturation front has reached i m for two example soils

**Summary and conclusions**

In several studies the infiltration of water from an inundated area into an aquifer below the inundation area is modelled by applying saturated flow boundary conditions. A prerequisite of this assumption is that the pore space under the inundated area is initially saturated or at least the saturation front under the inundation area proceeds fast enough to justify neglecting the variably saturated flow processes. In order to estimate how long it takes to saturate the pore space under an inundated area, an analytical infiltration model has been used for two different example soils under different initial conditions. While many infiltration-related studies address hydrological questions or irrigation problems where the infiltration rate or the cumulative infiltration is of major interest, this study presents diagrams with the position of the saturation front over time. These diagrams can help practitioners dealing with infiltration to decide whether the simplification of the process of infiltration under ponded condition as a saturated flow is feasible or a more advanced approach like the Richards's model is necessary. The motivation for this work arises from the predictive modelling of subsurface floods, but the results of this study might also be of interest for other applications where infiltrations may play a main role, for example the infiltration of contaminants from copper mining waste dumps (see *e.g.* Salazar, 2011).

The results show that the saturation of several meters of sandy soil below an inundation takes place within the order of hours. For many groundwater models this will cover only a few simulation time steps, so the assumption of saturation under an inundated area will apply for groundwater models of sand or gravel aquifers. The higher the ponding depth, the faster does the saturation front proceed, however, the impact of increasing the ponding depth on the progression of the saturation front reduces for large ponding depths.

For the clay soil example, orders of days were found to be the time period necessary to saturate one meter of soil. If the clay is initially dry, it takes even weeks to saturate a clay soil of one meter. Thus, the assumption of saturation below an inundated area will not be applicable for aquifers that are covered with a dry clay layer. However, if the clay layer is initially wet, the saturation front moves through a clay layer of one meter thickness comparatively fast, but it still takes several hours.

For the analyzed sand example, the ponding depth does not only influence the velocity of the progression of the saturation front, it also has a remarkable influence on the infiltration rate, while for initially wet clay, the infiltration rates were found to change only marginally with the ponding depth. The infiltration into sand occurs with a piston-shaped saturation front nearly without a transition zone. During infiltration into clay a larger fraction of the infiltrated water is distributed in the transition zone that precedes the saturation front. Thus, the arrival time of infiltrating water at a specific position, *e.g.* the groundwater level, can be estimated roughly by the progression of the saturation front for sandy soil. This is not possible for clay, because a distinct transition zone precedes the saturation front here.

As the assumption of saturated conditions under an inundated area does not apply in general, applicable approaches to handle the infiltration from an inundated area for groundwater models are desirable. With such modelling features the opportunities to perform predictive simulations of subsurface flood events will improve.

**References**

Barry, D.A., Parlange, J.Y., Haverkamp, R. and Ross, P.J. (1995). Infiltration under ponded conditions: 4. An explicit predictive infiltration formula. *Soil Science* **160**(1): 8-7. [ Links ]

Bauer, P., Gumbricht, T. and Kinzelbach, W. (2006). A regional coupled surface water/groundwater model of the Okavango Delta, Botswana. *Water Resources Research* **42,** No. W04403 [ Links ]

Becker, B.P.J. (2010). *Zur gekoppelten numerischen Modellierung von unterirdischem Hochwasser.* Dissertation RWTH Aachen University [ Links ]

Becker, B.P.J., Forberig, S., Flogel, R., Schüttrumpf, H. and Kongeter, J. (2012). On the determination of groundwater levels for hazard maps of groundwater head rise induced by high water. *Wasser Wirtschaft* **12,** 10-16 [ Links ]

Beven, K. and Germann, P. (1982). Macropore and water flow in soils. *Water Resources Research* **18** (5): 1311-1325. [ Links ]

Chong, S.K. and Green, R.E. (1983). Sorptivity measurement and its application. In: *Proceedings of the National Conference on Advances in Infiltration,* American Society of Agricultural Engineers, 12-13 Dec. 1983, Chicago, Illinois, 82-91. [ Links ]

Freeze, R. A. (1972). Role of Subsurface Flow in Generating Surface Runoff. 1. Base Flow Contributions to Channel Flow. *Water Resources Research* **8,** No. 3, 609-623. [ Links ]

Freeze, R.A. and Cherry, J.A. (1979). *Groundwater.* Prentice Hall, Englewood Cliffs, New Jersey. [ Links ]

Gillham, R.W. (1984). The capillary fringe and its effect on water table response. *Journal of Hydrology* **67,** 307-324. [ Links ]

Green, W.H. and Ampt, G.A. (1911). Studies on Soil Physics: 1. Flow of Air and Water through Soils. *The Journal of Agricultural Science* 4, 1-4. [ Links ]

Haverkamp, R., Parlange, J.Y., Starr, J.L., Schmitz, G. and Fuentes, C. (1990). Infiltration under ponded conditions: 3. A predictive equation based on physical parameters. *Soil Science* 149, No. 5, 292-300. [ Links ]

Homann, C. (2006). *Automatische Optimierung von Grundwasserhaltungen in von Hochwasser beeinflussten Aquiferen.* Dissertation RWTH Aachen University. [ Links ]

Huber, G., Hiller, G. and Braune, A. (2003). Konzepte des Hochwasserschutzes für die Bauten des Freistaates Sachsen im Historischen Stadtkern von Dresden. In: Umweltamt Landeshauptstadt Dresden and Dresdner Grundwasserforschungszentrum e.V. (eds.), *Hochwassernachsorge Grundwasser Dresden,* Wissenschaftliche Tagung zum BMBF-Forschungsprojekt/8.Oktober 2003/ Dresden, Rathaus/Tagungsband. Dresden Druckerei und Verlag Christoph Hille, 57-62. [ Links ]

Koch, A., Meon, G. and Leibundgut, C. (2007). Die gekoppelte Simulation der 2D-Oberflachenwasser- und 3D-Grundwasserstr0mung im Rahmen eines Expertensystems zur Modellierung der hydrodynamischen Wechselwirkungen einer Flusslandschaft. In: Meon, G., Sauter, Führbôter and Wolf (eds.): *Niedersachsisches Grundwasserkolloquium* 21.-22. Februar 2007 in Braunschweig. Grundwasser, Hochwasser und Salzwasser/Konzepte zur Risikominimierung. Stuttgart: E. Schweizerbart'sche Verlagsbuchhandlung (Nägele u. Obermiller), Zentralblatt für Geologie und Paläontologie/Teil I Allgemeine, Angewandte, Regionale und Historische Geologie Vol. 2007, No. 1/2, 99-114. [ Links ]

Kolymbas, D. (1998). *Geotechnik: Bodenmechanik und Grundbau.* Springer, Berlin. [ Links ]

Kreibich, H. and Thieken, A.H. (2008). Assessment of damage caused by high groundwater inundation. *Water Resources Research* 44, No. 9. http://dx.doi.org/10.1029/2007WR006621. [ Links ]

Monninkhoff, B. and Kernbach, K. (2006). Coupled surface water - groundwater modeling for planning of flood retention in the Lower Havel area. In: Proceedings of the: *International* FEFLOW *User Conference,* September 10-15, 2006 Berlin, Germany, 115-124. [ Links ]

Peetz, J.-V. and Steckel, B. (2007). 3-code coupling in flood simulation. In: *MpCCI 8th user forum* proceedings, february 13 and 14, 2007, 108-116. http://www.mpcci.de/fileadmin/mpcci/Userforum/MpCCI_8th_UserForum.pdf [ Links ]

Philip, J.R. (1957). Theory of infiltration: 1. The infiltration equation and its solution. *Soil Science* 83(5), 345-358. [ Links ]

Philip, J.R. (1969). *Theory of infiltration.* Advances in Hydroscience 5, Academic Press, 215-296. [ Links ]

Rawls, J.W. and Brakensiek, D.L. (1983). A procedure to predict Green and Ampt infiltration parameters. In: *Proceedings of the National Conference on Advances in Infiltration,* American Society of Agricultural Engineers: 12.-13. Dec. 1983, Chicago, Illinois, 102-112. [ Links ]

Salazar, I. (2011). *Beitrag zur Verringerung der Schadstoffkonzentration in heterogenen porosen Medien.* Dissertation RWTH Aachen University. [ Links ]

Schaap, G. M. (2009). Rosetta: A computer program for estimating soil hydraulic parameters with hierarchical pedotransfer functions. url: http://ag.arizona.edu/research/rosetta/download/rosetta.pdf, checked 2009-12-10. [ Links ]

Schaap, M.G. and van Genuchten, M.T. (2006). A modified Mualem-van Genuchten formulation for improved description of the hydraulic conductivity near saturation. *Vadose Zone Journal* 5, No. 1, 27-34. [ Links ]

Schmitz, G., Haverkamp, R., Parlange, J.Y. and Liedl, R. (1989). Infiltration bei eingestauter Bodenoberflaeche. 1. Analytische Lôsung der Infiltrationsgleichung und praktische Berechnung der Versickerung. *Zeitschriftfür Kulturtechnik und Landentwicklung* 30, No. 2, 121-131. [ Links ]

Sinaba, B. (2007). *Zur Infiltration bei eingestauter Bodenoberflache.* Diplomarbeit RWTH Aachen University SLUG Sächsisches Landesamt für Umwelt und Geologie (2003). *Einfluss des August Hochwassers 2002 auf das Grundwasser.* Materialien zur Wasserwirtschaft. [ Links ]

Sommer, T. and Ullrich, K. (2004). Das unsichtbare Hochwasser - Auswirkungen des August-Hochwassers 2002 auf das Grundwasser im Stadtgebiet von Dresden. In: 34. IWASA *Internationales Wasserbau-Symposium Aachen* 2004 : Hochwasserschutz - eindeichen oder ausweichen. Mitteilungen des Inst. f. Wasserbau und Wasserwirtschaft, RWTH Aachen. Aachen: Shaker, 203-224. [ Links ]

Sommer, T. and Eulitz, K. (2007). Auswirkungen des Hochwassers auf das Grundwasser/Phänomene und Modellansätze. In: Meon, Sauter, Führbôter and Wolf (eds.): *Niedersachsisches Grundwasserkolloquium* 21.-22. Februar 2007 in Braunschweig: Grundwasser, Hochwasser und Salzwasser/Konzepte zur Risikominimierung. Stuttgart: E. Schweizerbart'sche Verlagsbuchhandlung (Nägele u. Obermiller), Zentralblatt für Geologie und Paläontologie/Teil I Allgemeine, Angewandte, Regionale und Historische Geologie Vol. 2007, No. 1/2, 7-15. [ Links ]

van Linn, A. (2006). *Automatische Optimierung zur Bewertung* USDA United States Department of Agriculture (2009). USDA Risikoanalyse einer Hochwasserschutzma fi nahme. Dissertation RWTH Aachen University. Textural Classification Chart. url: http://www.gintsoftware.com/Jdownloads/reports/usda_textural_class.pdf, checked 2009-12-10. [ Links ]

van Genuchten, R. (1978). *Calculating the unsaturated hydraulic conductivity with a new closed-form analytical model.* Water Resources Program: Department of Civil Engineering, Princeton University. [ Links ]

Fecha de entrega: 29 de junio 2012 Fecha de aceptación: 29 de enero 2013.