Flujo no saturado

Zona de suelo

Zona no Zona
saturada o intermedia
vadosa

Zona
capilar

Zona
saturada
 Relaciones suelo - agua

1. Sólida
Fases: 2. Líquida
3. Gaseosa

Volumen de control: volumen sólidos + volumen agua +

volumen aire y vapor de agua
 Vg Aire Wg≈ 0
Vv
Vw Agua Ww
V W
Vs Suelo Ws

Contenido volumétrico de humedad: θ = Vw V

Volumen de vacíos: Vv = V − Vs = Vw + Vg
Porosidad del suelo: n = Vv V (n : 0.12 → 0.55)
Relación de vacíos: e = Vv V s
Grado de saturación: S = Vw Vv
 Vg Aire Wg≈ 0
Vv
Vw Agua Ww
V W
Vs Suelo Ws

θ = Vw V Vv = V − Vs = Vw + Vg n = Vv V S = Vw Vv e = Vv V s

→ θ = nS

Contenido de humedad en peso: w = Ww Ws

Peso unitario (densidad): Gravedades específicas:

Total: γT =W V Masa: Gm = γ T γ 0
Sólidos: γ s = Ws Vs Sólidos: Gw = γ w γ 0
Agua: γ w = Ww Vw Agua: G =γs γ0
Seco: γ seco = Ws V γ 0 = densidad de agua a 4°C
 FACTORES QUE AFECTAN EL MOVIMIENTO DEL
AGUA

1. Textura: distribución de tamaños de partículas del suelo

2. Estructura: ordenamiento y agregación de las partículas
definen el tamaño de poros
1. Poros capilares
2. Poros no capilares
3. Clase de minerales de arcilla
4. Cantidad de materia orgánica
Flujo en medios porosos
 Triángulo Arcilla: < 0.002 mm
Textural Limo: 0.002 – 0.05 mm
Arena: 0.05 – 2 mm

Arena muy fina: 0.05 – 0.1 mm

Lim
Arcilloso
Arena fina: 0.1 – 0.25 mm

o(%
)

Arena media: 0.25 – 0.5 mm

(%

)
Arcillo
la

Arcillo limoso Arena gruesa: 0.5 – 1 mm

cil

arenoso
Ar

Franco Franco
Arena muy gruesa: 1 – 2 mm
Franco arcillo- arcilloso arcillo-limoso
arenoso
Franco Franco
Franco limoso
Ar
fra enoso arenoso Limo
n
Arenoso co

Arena (%)
 Agua en el suelo
Fenómenos (adsorción, tensión superficial, ósmosis)
a) Adsorción: El agua puede ser retenida por superficies sólidas,
generalmente mediante fuerzas electrostáticas. Estas fuerzas
pueden ser fuertes, y en arcillas pueden ser responsables de la
retención significativa de agua contra fuerzas de succión altas.
b) Tensión superficial: La fuerza de tensión superficial es
responsable del ascenso capilar de agua en un tubo circular de
diámetro pequeño. La altura de la columna capilar de agua aumenta
en forma inversamente proporcional al diámetro del tubo.
Altura
capilar es paire
inversa al
radio
capilar
p B = p C = p aire = 0
p A = p B − ρ gh 2

p capilar = p aire − p A
 Agua en el suelo
El agua es almacenada en el espacio poroso del suelo. La
estructura de los poros es más compleja que el caso del
tubo simble, pero la capilaridad es un efecto igualmente
importante. La magnitud de las fuerzas de capilaridad
depende del tamaño de los poros, i.e. poros pequeños
pueden desarrollar fuerzas de capilaridad altas.
 Agua en el suelo
Fenómenos
Como en un tubo capilar, presiones de succión alta pueden
desarrollarse en el suelo - hasta del orden de 10000 cm H2O. Por lo
tanto es conveniente usar una escala logarítmica para describir la
succión en el suelo. Se utiliza la escala pF: pF = 1 (-10 cmH2O),
pF = 2 (-100 cmH2O), pF = 3 (-1000 cmH2O), etc.
b) Osmosis: Si una membrana es permeable para agua pura, pero
impermeable para una solución de agua, y si tal membrana semi-
permeable separa una solución de agua pura, entonces el agua se
moverá a través de la membrana (por difusión o osmosis),
diluyendo la solución. Membrana

Solución débil Solución fuerte

Agua
Dilución
Por lo tanto la diferencia en las soluciones através de la membrana
genera una presión efectiva. Para contrarrestar el movimiento debe
aplicarse una presión a la solución más fuerte (presión osmótica).
 Agua en el suelo
Fuerzas gobernantes
a) En un suelo saturado: El movimiento del agua es controlado por
fuerzas que provienen de la distribución hidráulica de presión
(positiva, +) y de la influencia de la gravedad.
b) En un suelo no saturado: Los fenómenos de adsorción y tensión
superficial ejercen fuerzas que se oponen al drenaje de agua del
suelo. (En la práctica no siempre es posible distinguir entre las dos
fuerzas, de efecto similar). La adsorción en general solo es
importante para suelos arcillosos secos. Los efectos de ósmosis en
el suelo surjen predominantemente debido a la acción de las raíces.
Asi las principales fuerzas gobernantes son las de capilaridad y
gravedad.
El efecto de la capilaridad es producir succión, o presión negativa
(con relación a la presión atmosférica). El flujo no saturado está
controlado por una distribución de presión negativa (-) y fuerzas
gravitacionales.
 1. Gravedad, G
Fuerzas:
2. Capilaridad, C

Infiltración y exfiltración
Zona de suelo
G C G C

Zona no Zona Percolación

saturada o intermedia
(G domina)
vadosa

Zona C nuevamente importante

capilar G C
Zona
saturada
 Agua en el suelo
Potencial gravitacional
Expresando potencial como energía potencial / unidad de volumen
Φ g = ρw ⋅ g ⋅ z
i.e. el trabajo realizado contra la gravedad, donde z es la distancia
vertical sobre un datum específico.
Potencial de presión
a) Si el agua está bajo una superficie libre, la presión hidrostática
es positiva con respecto a la presión atmosférica. La presión
hidrostática está dada por, Φ p = ρ w ⋅ g ⋅ψ donde ψ es la
profundidad bajo la superficie libre.
b) Si el agua está a una presión menor a la atmosférica, por ejemplo
debido a capilaridad, la presión se denomina presión de succión o
tensión. El potencial de presión hidrostático (negativo) está dado
por − ρ w ⋅ g ⋅ ψ donde ψ es la altura sobre la superficie libre.
 Agua en el suelo
Potencial de presión
-Ψ2

Ψ2 - Potencial presión negativo

Potencial presión (+) + P. atmosférica

Ψ1
Ψ1

Si existe una diferencia en potencial total Φ t , entonces existirá

flujo del sitio de potencial alto al de potencial bajo.
Bajo equilibrio hidrostático no hay movimiento y por lo tanto no
existe un gradiente de potencial total, i.e. Φ g (incrementando con z)
se balancea con Φ p (disminuyendo con z).
El potencial es comúnmente descrito en forma más simple por una
altura equivalente de agua (cabeza),
Φ t = ρ w ⋅ g ⋅ψ + ρ w ⋅ g ⋅ z
h (cab. total) = ψ (cab. presión) + z (cab. posición)
 Estado energético del agua intersticial

El agua en el suelo se mueve desde puntos donde tiene un nivel

energético mayor hacia puntos con nivel energético menor.
Φt = potencial total = estado energético del agua en la matriz de
suelo
Φ t = Φ p + Φ g + Φ osm + Φ gas
Φp = potencial matricial debido a las fuerzas interactivas locales
entre agua y suelo
Φg = potencial gravitacional por acción de la gravedad
Φosm = potencial osmótico por diferencias de concentración de sales
Φgas = potencial neumático por cambios de la presión externa del
aire incluido en la matriz de suelo
Definidos con respecto a un estado energético de referencia del
agua:
superficie de agua libre bajo condiciones atmosféricas
1. Presión atmosférica
2. 20 grados centígrados
3. Elevación relativa < DATUM cero

Usualmente se considera:
Φt = Φ p + Φ g
porque:
ψosm ∼ = 0 bajo suposición de agua homogénea
ψgas ∼ ψatm Æ ψgas = 0 por nivel energético de referencia
 Potencial en términos de cabeza hidráulica de agua:

ψg = distancia vertical desde el DATUM hasta el punto

= -z si se toma el eje vertical z positivo hacia arriba con cero
en el DATUM

ψ es negativo pues se requiere hacer trabajo en contra de las

fuerzas matriciales para drenar agua

H = Ψ =ψ − z
Cabeza gravitacional

Cabeza de presión de agua intersticial

Cabeza hidráulica
 Potenciales (cabezas) de presión y gravitacional

superficie
DATUM

ψg = -z Suelo no saturado
Ψ=H
ψ
Probeta porosa de
cerámica

Ψ = f ( textura del suelo, θ )

h = -ψ = succión
DATUM superficie

θ1

θ2 > θ1
Suelo no saturado

OJO: Se insinúa una relación inversa entre h y θ

 Agua en el suelo

Perfil de equilibrio de humedad del suelo

a) Para suelo uniforme:
Para que el perfil de humedad del suelo no saturado sobre la tabla de
agua esté en equilibrio, el potencial de presión debe disminuir
con el incremento de la altura sobre la tabla de agua de tal forma
que el potencial total (pot. de presión + gravitacional) se
mantenga constante.
i.e. si z es la altura sobre la tabla de agua.

Por lo tanto la distribución del contenido Ζ

de humedad es la misma característica de (=-Ψ )
humedad con -Ψ reemplazada por z
ψ = −z Tabla de agua

θ
 Agua en el suelo
Perfil de equilibrio de humedad del suelo
El suelo estará efectivamente saturado hasta una altura sobre la tabla
de agua al cual la succión es igual valor de succión de entrada de
aire. Esta altura de material saturado encima de la tabla de agua
(definida por Ψ = 0) se conoce como la franja de capilaridad.
a) Para suelo no-uniforme: Si el suelo contiene capas horizontales
con diferentes curvas características de humedad, la distribución de
humedad de equilibrio tendrá una forma más compleja. El potencial
de presión deberá satisfacer sin embargo la condición, ψ = − z de
tal forma que el gradiente de potencial total sea cero.
Pero para un valor dado de tensión, las
curvas características diferentes darán Ζ
Arena
contenidos de humedad diferentes, (=-Ψ )
ocurriendo discontinuidades en el Arcilla
perfil de humedad en las interfases. Tabla de agua

θ
 Contenido de humedad del suelo

Se utilizan comúnmente tres términos para describir condiciones

limitantes de la humedad del suelo:
a) Capacidad de campo
b) Punto de marchitamiento
c) Déficit de humedad del suelo DHS
 Contenidos de humedad del suelo
Si la humedad en el suelo en la zona radicular es muy baja la
presión de succión de las raíces (h) debe ser alta
⇒ h alta implica menor capacidad de bombeo de la raíz

El bombeo tiene un límite que corresponde al punto de

marchitamiento : contenido de humedad para el cual la raiz no
puede sacar más agua y la planta entra en estado de stress. El punto
de marchitamiento define el límite inferior de humedad del suelo
disponible para la planta.
El stress se define como el estado en el cual la planta se marchita y
no se recupera cuando las hojas se mantienen en una atmósfera
húmeda durante la noche. El punto de marchitamiento ocurre más o
menos a una presión de succión de 15 atmósferas (depende del tipo
de planta y condiciones atmosféricas).
 Contenido de humedad del suelo
Se han definido otros dos estados de contenido de humedad del
suelo, Capacidad de campo: “la cantidad de agua contenida en el
suelo después de que el exceso de agua ha drenado y la tasa de flujo
vertical ha decrecido materialmente”. Siguiendo la saturación del
perfil del suelo, después de dos o tres días de drenaje, el perfil
alcanzará un estado en el cual cambios adicionales de humedad
ocurren lentamente y el drenaje gravitacional es suficientemente
pequeño y puede ser despreciado. Si subsecuentemente más agua
entra al perfil de suelo (cuando está a la capacidad de campo) el
volumen adicional de agua pasará directamente como drenaje
gravitacional.
⇒ La capacidad de campo define el límite superior de contenido de
humedad del suelo al cual el agua está disponible para las plantas.
Déficit de humedad del suelo: la magnitud de humedad a la cual el
suelo está más seco que la capacidad de campo, e.g. Si DHS = 50
mm implica que se requieren 50 mm de agua para que el perfil de
suelo retorne a la capacidad de campo.
Medición del contenido de humedad y del potencial
Medición del contenido de humedad del suelo
a) Determinación gravimétrica
b) Bombardeo de neutrones
c) Atenuación de rayos Gamma.
Medición del potencial
a) Tensiómetro
b) Bloques de resistencia eléctrica
c) Psicrómetro de suelo.
 Agua en el suelo
Tensiómetro
 Agua en el suelo
Flujo no saturado
Los principios fundamentales de flujo saturado son aplicables a flujo
no saturado.
a) Potencial: Como en flujo saturado, ocurre flujo de una región de
potencial alto a una de potencial bajo, pero para flujo no saturado el
potencial de presión es negativo debido a tensión en el agua del
suelo, i.e., Φ = z + ψ donde ψ es (-)

ψ A =ψ B ψC > ψB ψ C << ψ D
z B > z A zC > z B z D > zC
Flujo B → A No Flujo Flujo C → D
(Φ B > Φ A ) (Φ C = Φ B ) (Φ C > Φ D )
 Agua en el suelo
Flujo no saturado
Ejemplo. Un tiempo después de una lluvia intensa sobre un
suelo bien drenado con el nivel freático a 1.3m por debajo
de la superficie del terreno se midió el siguiente registro de
potencial matricial (metros de cabeza de agua) mediante
tensiómetros ubicados a las profundidades indicadas:
z(m)-0.05-0.10-0.20-0.30-0.40-0.50-0.60-0.70-0.80-1.0-1.30
ψ(m)-1.85-1.62-1.34-1.10-0.89-071-0.59-0.49-0.40-0.23 0.00
Grafique la cabeza de succión y la cabeza piezométrica
con la profundidad y de éstas gráficas deduzca como está
drenando el perfil de suelo. Si la conductividad hidráulica
saturada es de 0.4 md-1 calcule el flujo de agua a través del
nivel freático.
 Agua en el suelo
Flujo no saturado
Ejemplo. El potencial de agua en el suelo es monitoreado en una
parcela de pasto mediante cinco tensiómetros que están conectados a
un manómetro agua-sobre-mercurio conectados a un tanque común
de mercurio 50 cm sobre el suelo. Las profundidades de los
tensiómetros son 20, 40, 60, 80 y 100 cm por debajo de la superficie
del suelo. En un instante de tiempo dado las lecturas (altura de
mercurio bajo el nivel del suelo) son 23.0, 11.9, 12.3, 13.5 y 15.1 cm
respectivamente.
a) Evalúe y esquematice la distribución de potencial total con la
profundidad y por lo tanto comente sobre la respuesta del perfil.
b) Si la conductividad hidráulica no saturada K (cms-1) está
representada por la función,
1
K= 3 para ψ > 10
10 + (ψ − 10) 3

donde ψ es la cabeza de presión (cm H2O) (ψ es la cabeza de

succión), estime el flujo en cm día-1 en la base del perfil medido.
 Curva de Retención de Humedad

Saturado Æ presión agua intersticial ≥ 0

Suelo
No saturado Æ presión agua intersticial < 0
(succión con h > 0)

Campana hermética con presión

inicial atmosférica

Muestra de suelo inicialmente saturado

Resultado: a mayor succión la muestra

libera más agua y por lo tanto menor es el
Reductor de presión
contenido de humedad remanente
 Agua en el suelo
Curva característica de retención de humedad del suelo

Considere una muestra de suelo (no expansivo), inicialmente

saturado, con agua. Entonces la presión de poros > p. atmosférica.
Si se remueve agua, se alcanza el estado para el cual el suelo está
apenas saturado. La presión de poros en la superficie en este caso
es cero (p. atmosférica).
Si se aplica una succión pequeña adicional, la respuesta del agua en
los poros alcanza el estado en el cual fuerzas de tensión superficial
actuarán para resistir la succión aplicada. Si la succión aplicada se
incrementa, las fuerzas de tensión en los poros más grandes se
vencerán, y drenará agua del suelo (el valor se conoce como
succión de entrada de aire).
 Agua en el suelo
Curva característica de retención de humedad del suelo
A medida que se incrementa la succión, agua adicional saldrá de
los poros más grandes, dejando solo los poros más finos (de menor
radio y mayores fuerzas de tensión superficial) llenos. Por lo tanto
a mayor succión, menor contenido de agua, y la forma de la
relación entre las dos variables reflejará la distribución del tamaño
de los poros.
Contenido de
humedad
volumétrica

0
Entrada
Presión de poros Ψ de aire
 Agua en el suelo
Curva característica de retención de humedad del suelo
Diferencias en los tipos de suelo se reflejan en diferencias en la
curva característica de humedad.
Contenido de
humedad
volumétrica

θ
Suelo
arcilloso
Suelo
arenoso
0

Presión de poros Ψ
La arcilla generalmente tiene mayor contenido de humedad debido
a su mayor porosidad, pero para tensiones moderadas de agua del
suelo, el suelo arenoso libera mucha más agua, debido
generalmente a los espacios mayores entre poros.
 107
arcilla
106
h
(cm) 105
limo
104

103
arena

102

101
0 θ 0.5
Experimento alternativo: mismo esquema pero con muestra
inicialmente seca la cual se satura a medida que se reduce la
succión
Resultado: gráficas similares pero
no iguales

Histéresis: contenido de
humedad diferente para igual
valor de succión

Explicación: se requiere más energía para drenar agua que

para incorporarla en la matriz de suelo
 Agua en el suelo
Curva característica de humedad del suelo
Si se considera el proceso inverso de re-humedad, o sorción, se
obtiene la siguiente curva característica:
Contenido de
humedad
volumétrica

0
Presión de poros Ψ

La curva característica de humedad exhibe histéresis, - una relación

no-única que depende de la historia previa de humedad y secado.
La histéresis puede explicarse por atrapamiento de aire y por el así
llamado efecto “botella” que surje debido a que muchos poros tienen
entradas relativamente estrechas. Asi, se alcanza un contenido de
humedad dado a menor succión durante el llenado de los poros.
 Agua en el suelo
Curva característica de humedad del suelo
La presencia de histéresis es una complicación mayor en el análisis
de la humedad del suelo. La curva de secado, seguida por un periodo
de tiempo considerable después de una respuesta rápida de
incremento de humedad por lluvia o irrigación, es (+) importante.
 107

106
h salida
(cm) 105

104
entrada
103

102

101
0 θ 0.5
 Descripción analítica de las curvas de retención de humedad

1. Feddes:

a1 (b1 −θ )
h=e para θ1 ≤ θ ≤ θ s
a2 (b2 −θ )
h=e para θ 21 ≤ θ ≤ θ1s
a3 (b3 −θ )
h=e para θ ≤ θ 2 s
 2. Brooks y Corey:
107
−1 m
h = hs s h
106

(cm) 105
m = índice de distribución de
104
poros 1
m
103

hs 102

101
1 10-1 s 10-2
 Curva de conductividad hidráulica
K es función del medio poroso y del fluido
Si condición saturada Æ vacíos están llenos de agua y Q
es máximo

Si condición no saturada Æ vacíos están llenos de agua y

aire. Éste va a representar una resistencia al flujo de agua y
además el área efectiva de flujo va a ser menor Æ Q es
menor. Entre menor sea el contenido de humedad Æ
menor es la conductividad:

K = f (θ ) θ = f (h )
pero

⇒ K = f (h )
 102
arena

101
K
(cm/día) 100
arcilla
10-1

limo
10-2

10-3

10-4 0
10 103 h (cm) 107
 Descripción analítica de las curvas de conductividad
hidráulica
1. Feddes:

K = Ks si h ≤ ha
− n ( h − hs )
K = K se si ha ≤ h ≤ hlim
K = ah −1.4 si h ≥ hlim
 2. Brooks y Corey:
10-4

K (θ ) = K s s c
K 10-5

(cm/s)
10-6
c = índice de conexión de poros
10-7
c

2 + 3m 10-8
c=
1

m 10-9

10-10
1 10-1 s 10-2
 Agua en el suelo
Flujo no saturado
b) Conductividad: Para flujo saturado Ksat = cte; todos los poros
están llenos de agua. Para flujo no saturado, K decrece a medida que
θ decrece, i.e. K = K(θ ). Los poros más grandes se desocupan
primero, dejando agua en los poros menores solamente.

log (K(θ))

log (-ψ )

Puede ocurrir un decrecimiento importante de K a succiones

relativamente bajas. Debido a la baja conductividad a succión alta, se
pueden requerir gradientes de potencial muy grandes para obtener
fujo significativo.
 Agua en el suelo
Flujo no saturado
Relaciones K versus θ para diferentes suelos pueden dar
aparentemente resultados sorprendentes, e.g. Ksat(arena) >> Ksat(arcilla)
(poros más grandes, tasas de transmisión mayores en arena).
Pero a medida que la succión aumenta, los poros más grandes se
desocupan primero y por lo tanto puede ocurrir que K(θ )(arena) <<
K(θ )(arcilla).
Arena a succión alta puede ser una capa impermeable.

Arena K(θ) Arena

log (K(θ))
Arcilla Limo

log (-ψ ) θ
 Agua en el suelo
Flujo no saturado
c) Ley de Darcy para material no saturado: Existe cierta discusión
de la aplicabilidad de la ley de Darcy a baja velocidad. Es un área
difícil en la cual obtener buenos datos experimentales.
Darcy para flujo saturado, si Kx = Ky = Kz = K:
∂Φ
Vz = − K etc
∂z
Extensión a flujo no saturado,
∂Φ
Vz = − K (θ ) etc., donde θ = contenido de humedad (1)
∂z
Potencial,
Φ =ψ + z donde ψ = succión (negativa) (2)
Continuidad
∂θ ⎛ ∂V ⎞
= −⎜ z ⎟ etc. (3)
∂t ⎝ ∂z ⎠
∂θ ∂ ⎛ ∂Φ ⎞
De (1) y (3), = ⎜ K (θ ) ⎟ etc.
∂t ∂z ⎝ ∂z ⎠
 Agua en el suelo
Flujo no saturado
De (2), ∂θ ∂ ⎛ ∂ (ψ + z ) ⎞
= ⎜ K (θ ) ⎟
∂t ∂z ⎝ ∂z ⎠
Y por lo tanto,
∂θ ∂ ⎛ ∂ψ ⎞ ∂K
= ⎜ K (θ ) ⎟+
∂t ∂z ⎝ ∂z ⎠ ∂z
Ecuación válida para infiltración vertical uni-dimensional.
Ahora, ∂ψ dψ ∂θ
=
∂z dθ ∂z
Gradiente de humedad
dθ
Pendiente curva caracaterística = C (θ )
dψ
∂ψ 1 ∂θ ∂ψ K (θ ) ∂θ
Por lo tanto = ó K (θ ) = Difusividad =
∂z C (θ ) ∂z ∂z C (θ ) ∂z D (θ )
 Agua en el suelo
Flujo no saturado
Ecuación de Richards, válida para infiltración vertical uni-
dimensional,
∂θ ∂ ⎛ ∂θ ⎞ ∂K (θ )
= ⎜ D(θ ) ⎟+
∂t ∂z ⎝ ∂z ⎠ ∂z