RELACION

AGUA - SUELO - PLANTA

[Link]. Luis Rázuri Ramírez

Profesor Titular
Luis Rázuri Ramírez
CIDIAT-ULA 1
 CONCEPTOS BÁSICOS SOBRE SUELO-AGUA

El suelo como sistema de tres fases

El término suelo se refiere a la capa más superficial de la

corteza terrestre formado por material fragmentado y
meteorizado.

En la formación de un suelo intervienen procesos físicos de

degradación de las rocas y minerales, así como procesos
biológicos que aportan materia orgánica y contribuyen a la
evolución continua del sistema.

Desde un punto de vista físico el suelo es el resultado del

equilibrio entre tres fases: sólida, líquida y gaseosa.
Luis Rázuri Ramírez
CIDIAT-ULA 2
 CONCEPTOS BÁSICOS SOBRE SUELO-AGUA

En la mayoría de los casos los sólidos, que incluyen tanto

partículas orgánicas como inorgánicas, constituyen entre el
35-75% del total del volumen del suelo.

El resto del volumen es espacio hueco en forma de capilares,

poros o fisuras, que puede ser llenado por aire o agua.

Los sólidos, aire y el agua son fundamentales para el óptimo

desarrollo de las plantas y la proporción de cada componente
en cada momento depende de factores edáficos (del suelo
como tal), climáticos y humanos.

Luis Rázuri Ramírez

CIDIAT-ULA 3
 CONCEPTOS BÁSICOS SOBRE SUELO-AGUA

La fracción sólida del suelo se compone de una fracción

mineral y otra orgánica. Cerca de la superficie una parte
substancial del volumen de suelo puede estar ocupado por
material vivo (raíces de plantas, gusanos, bacteria,
microflora, microfauna, etc).

La proporción y actividad de estos elementos depende

obviamente del suelo, tipo de cubierta, época del año y
puede tener un efecto importante sobre la distribución y
movimiento del agua del suelo.

Luis Rázuri Ramírez

CIDIAT-ULA 4
 El sistema suelo-agua

 El suelo es un sistema constituido de tres

fases distintas: fase sólida, líquida y
gaseosa;
 Las fases líquidas y gaseosas son
complementarias, ocupando el mismo espacio.
 CONCEPTOS BÁSICOS SOBRE SUELO-AGUA

Los suelos, en virtud de su contenido en materia orgánica

(MO) y profundidad, se clasifican en:

suelos orgánicos (>30% MO en peso, >40cm espesor)

suelos inorgánicos (el resto)

Aunque miles de hectáreas de suelos orgánicos se dedican

actualmente al cultivo, la mayoría de los suelos agrícolas
son suelos inorgánicos o minerales.

El rango normal de contenido de materia orgánica de los

suelos agrícolas varía de < 1% a 5%.

Los suelos de las regiones áridas tienen contenidos bajos de

MO (<1%).
Luis Rázuri Ramírez
CIDIAT-ULA 6
Los suelos inorgánicos contienen también materia
orgánica en forma de humus
El humus es la sustancia compuesta por productos
orgánicos, de naturaleza coloidal, que proviene de la
descomposición de los restos orgánicos (hongos y
bacterias).
Se caracteriza por su color negruzco debido a la gran
cantidad de carbono que contiene. Se encuentra
principalmente en las partes altas de los suelos con
actividad orgánica

Luis Rázuri Ramírez

CIDIAT-ULA 7
 El perfil del suelo

La parte visible del suelo, muy importante, es la

superficie. El examen de esta zona revelará procesos que
ocurren en dicha superficie, pero no describen
necesariamente el carácter del suelo en su conjunto.
 El perfil del suelo

Para ello hemos de examinar el suelo en profundidad (por

ejemplo abriendo una calicata). La sección vertical de un
suelo se denomina perfil del suelo.

El perfil del suelo es raramente uniforme en profundidad y

típicamente se compone de una sucesión de capas o
estratos más o menos uniformes.

Cuando estas zonas diferenciadas han surgido de los

procesos de evolución interna del suelo (pedogénesis) se
denominan horizontes.
 El perfil del suelo
El horizonte superior de un suelo,
horizonte A, al ser el de mayor actividad
biológica se enriquece en materia
orgánica y suele presentar colores más
oscuros que los horizontes inferiores.
Este horizonte superior recibe agua
directamente en su cara exterior. El
agua altera, solubiliza y desplaza a capas
más inferiores algunos de sus materiales
(como arcillas o carbonatos) en el
proceso llamado eluviación.
A continuación se encuentra el horizonte B, donde se acumulan
algunos de los materiales que se lavan del horizonte A en el
proceso llamado iluviación.
Luis Rázuri Ramírez
CIDIAT-ULA 10
 El perfil del suelo
Bajo este horizonte se puede encontrarse el horizonte C
compuesto por el material parental del suelo alterado
(fragmentos meteorizados de la roca subyacente o sedimentos
antiguos).

Por último se encontraría la roca madre o algún material

sedimentario antiguo sobre los que se desarrollo el suelo
moderno.

La secuencia A-B-C se reconoce claramente en algunos casos,

siendo poco acusada e incluso indistinguible en otros.

Cada uno de los horizontes anteriores puede estar

caracterizado no sólo por cambios de color sino por cambios en
sus propiedades físicas, en especial las relacionadas con el
contenido y redistribución del agua en el suelo, de interés en
nuestro caso Luis Rázuri Ramírez
11
CIDIAT-ULA
Un caso particular en el perfil de un suelo lo constituye la
presencia de un horizonte duro y/o impermeable, duripan,
(hardpan) de distintos orígenes que a veces se forman en un
suelo y que cambian radicalmente su comportamiento
(encharcamientos, mala aireación de las raíces, etc)
Pueden formarse por acumulación de caliza en climas
áridos, por efecto del laboreo, o incluso por degradación
(pérdida de estructura) de las arcillas por salinidad.
La información obtenida de la cuidadosa descripción del
perfil del suelo será inestimable para realizar un buen
diseño y manejo del riego y del drenaje.

Luis Rázuri Ramírez

CIDIAT-ULA 12
 Fracción mineral del suelo. Textura

La fracción mineral del suelo se clasifica según el tamaño de

sus partículas aisladas.

La distribución del tamaño de las partículas de un suelo se

denomina textura.

Los principales tipos de partículas de un suelo son arcilla (A),

limo (L) y arena (a) que abarcan los tamaños desde 0 a 2 mm.

Partículas mayores de 2 mm no suelen considerarse en las

determinaciones analíticas y se denominan grava.

Existen varias clasificaciones texturales aunque la más

difundida actualmente es la USDA (United States Department
of Agriculture, 1951), dada en la Tabla siguiente
Luis Rázuri Ramírez
CIDIAT-ULA 13
 Fracción mineral del suelo. Textura
Rango de tamaños
Textura (mm)

Arcilla 0 – 0,002

Limo 0,002 – 0,05

Arena 0,05 – 2,0

Arena muy fina 0,05 – 0,1

Arena fina 0,1 – 0,2

Arena media 0,2 – 0,5

Arena gruesa 0,5 – 1,0

Arena muy gruesa 1,0 – 2,0

Grava > 2,0

Desde el punto de vista cualitativo la textura representa la

sensación al tacto (húmedo o seco) de un suelo determinado.
Desde un punto de vista cuantitativo, interpretado por una
clasificación textural, la textura representa un atributo
natural y permanente del suelo y uno de los más usados para
describir la composición física de un suelo.
Luis Rázuri Ramírez
CIDIAT-ULA 14
 Clasificación de los suelos

0,002 0,05 0,1 0,25 0,5 1,0 2,0 mm

muy fina fina mediana gruesa muy gruesa
arcilla limo grava
arena
Clasificación según USDA

0,002 0,02 0,2 2,0 mm

fina gruesa
arcilla limo grava
arena
Clasificación según ISSC
 Fracción mineral del suelo. Textura

La Figura, ilustra la relación entre los tamaños de las partículas del suelo

Figura. Partículas del suelo (Hillel,1994)

La fracción arcilla resulta mucho más determinante en el

comportamiento de un horizonte y en consecuencia del suelo y
de la repuesta de las plantas que las restantes fracciones, y
así se pone de manifiesto en la denominación de la clase
textural.

Debe haber unas tres veces más de arena gruesa que arcilla
para que el suelo tenga unas propiedades condicionadas por la
arena (Gaucher, 1971).
Luis Rázuri Ramírez
CIDIAT-ULA 16
 Fracción mineral del suelo. Textura

La distribución del tamaño de partículas o granulometría,

expresado como porcentaje en peso de las partículas de un
determinado tamaño, se determina mediante un análisis
granulométrico y se representa en un diagrama semilogarítmico.

Luis Rázuri Ramírez

CIDIAT-ULA 17
 Fracción mineral del suelo. Textura

El "porcentaje menor que d" en la

Figura, representa el porcentaje
que pasaría por un tamiz o luz de
malla de ese diámetro d.
Suelos con una concentración de
partículas dentro de un rango
pequeño se denominan suelos mal
gradados (por ejemplo el suelo
gravoso de la Figura).

Los suelos agrícolas se clasifican normalmente como, arenoso (a),

arcilloso (A), limoso (L), franco (F) o mezclas de éstos.
Estas designaciones texturales se determinan en base a los porcentajes
(en masa) de cada uno de los principales tipos de partículas leídos en el
triángulo textural , USDA
Luis Rázuri Ramírez
CIDIAT-ULA 18
Fracción mineral del suelo. Textura

Luis Rázuri Ramírez

CIDIAT-ULA 19
 Fracción mineral del suelo. Textura

Por ejemplo un suelo con un 30% de arena, 35% de limo y 35%

de arcilla se clasificaría como Franco Arcilloso (FA).
El término franco se refiere a un suelo con una composición
equilibrada de cada uno de los 3 tipos principales de partículas
y que presenta condiciones físicas óptimas para el cultivo.
A veces se emplean los términos "suelo pesado" y "suelo ligero"
para referirse a suelos arcillosos y arenosos respectivamente,
aunque son términos equívocos y se prefieren los términos
"suelo fuerte" y "suelo suelto" respectivamente.

Luis Rázuri Ramírez

CIDIAT-ULA 20
1 Triángulo textural

90 10
80 20
70 Muy 30
)

Lim
arcilloso
(%

60 40

o(
la

Arcilloso
cil

50 Arcillo- 50

%)
Ar

40 Arcillo- limoso 60
arenoso Franco-arcilloso Franco-arcillo-
30 Franco-arcillo- limoso 70
arenoso
20 80
Franco
10 Franco-arenoso Franco-limoso
90
Arenoso A. franco Limoso
0
0 90 80 70 60 50 40 30 20 10 1
Arena (%)
1
Ejemplo de puntos ubicados en el
triángulo textural

90 10
80 20
70 30
)
(%

Lim
60 40
la
cil

o(
50
Ar

50

%)
40 60
30 70
20 80
10 90
0
90 80 70 60 50 40 30 20 10
0 1
Arena (%)
 Coordenadas para localización de los
1 X=limo+arcilla*(cos 60)
puntos en el triángulo textural
Y=arcilla * (sen 60)

90 10
80 20
) 70 30
(%

Lim
60 40
la
cil

o(
Ar

50 50

%
)
40 60
30 70
20 Y 80
10 90
60º
0
X
(0,0) 80 70 60 50 40 30 20 10
90
0 0 1
Arena (%)
 Coordenadas de un punto dado en
el triángulo textural

Como se trata de un triángulo equilátero las

coordenadas cartesianas pueden ser fácilmente
determinadas utilizando relaciones trigonométricas

X  Limo  Arcilla  cos (60 ) 0

Y  Arcilla  sen (60 ) 0

X = distancia horizontal del punto en relación al origen

en el sistema cartesiano de coordenadas;
Y = ordenada vertical del punto en relación al origen en
el sistema cartesiano de coordenadas.
 Ejemplo de cálculo para ubicación de
un punto en el triangulo textural

Ejemplo.
Calcule las coordenadas cartesianas del punto
en el triángulo textural, correspondiente a 20%
de arcilla, 40% de limo y 40% de arena.

Solución.

X  Limo  Arcilla  cos (60o ) Y  Arcilla  sen (60o )

X  40  20  0,5 Y  20  3 / 2
X  50,0 Y  17,3
 1
Ubicación del punto en el
triangulo textural

90 10
80 20
70 Muy 30

Lim
arcilloso

(%
60 40

o(
la Arcilloso
cil
50 Arcillo- 50

%)
Ar

40 Arcillo- limoso 60
arenoso Franco-arcilloso Franco-arcillo-
30 Franco-arcillo- 70
limoso
arenoso
20 80
Franco
10 Franco-arenoso Franco-limoso
Y =17,3 90
Arenoso A. franco Limoso
0 0
90 80 70 60 50 40 30 20 10
0 1
Origen Arena (%)
(0,0) X = 50,0
0 50
 Ejemplos de clasificaciones en el
triangulo textural
N° Arcilla (%) Limo (%) Arena (%) Clasificación textural
1 4,0 84,0 12,0 Limoso
2 72,0 17,0 11,0 Muy arcilloso
3 18,0 64,0 18,0 Franco limoso
4 39,0 18,0 43,0 Franco arcilloso
5 31,0 52,0 17,0 Franco arcillo limoso
6 26,0 16,0 58,0 Franco arcillo arenoso
7 13,0 10,0 77,0 Franco arenoso
8 23,0 28,0 49,0 Franco
9 45,0 19,0 36,0 Arcilloso
10 41,0 42,0 17,0 Arcillo limoso
11 37,0 9,0 54,0 Arcillo arenoso
12 4,0 3,0 93,0 Arenoso
13 12,0 3,0 85,0 Areno francoso
Es importante señalar que la textura es una propiedad
fundamental de los suelos ya que permite inferir otras
propiedades como:
•la capacidad de retención de agua disponible para las plantas y
su suministro;
•facilidad de circulación del agua;
•facilidad para el laboreo;
•riesgo de formación de costra superficial que condiciona la
infiltración y la nascencia de cultivos;
•riesgo de erosión eólica e hídrica;
•capacidad para almacenar nutrientes y sales;
•capacidad para admitir aguas residuales y otros residuos
líquidos;
•orden de magnitud de la superficie específica.
Luis Rázuri Ramírez
CIDIAT-ULA 28
 Superficie específica de las partículas del suelo

Se denomina superficie específica ,s, a la superficie total

intersticial de los poros del suelo por unidad de masa del
mismo

s = As/Mt (L2M-1)

Suelos compuestos por partículas dispersas de pequeño

tamaño presentan un gran valor de s.

La superficie específica ha de entenderse como la superficie

expuesta a los procesos de intercambio físico y químico que
tienen lugar en el suelo, por lo que se relaciona con el carácter
reactivo del mismo.

La relación con el tamaño de partícula depende de la forma de

la misma.
Luis Rázuri Ramírez
CIDIAT-ULA 29
Para un medio compuesto de N partículas esféricas de radio R y
densidad ρ, se tiene:
s = a / m = 3 / ρ R

lo que muestra que la superficie específica de una cantidad arbitraria

de esferas idénticas es inversamente proporcional al radio de la
esfera.

Para el caso de partículas planas en forma de disco, expresadas como

cilindros de grosor T y radio R, y asumiendo que R>>T, se tiene:

s = a / m ≈ 2 / ρ T

luego en este caso el área específica es inversamente proporcional al

grosor de las partículas

Las ecuaciones anteriores se utilizan para estimar la superficie

específica de las partículas del suelo asociadas a las distintas clases
texturales.
Luis Rázuri Ramírez
CIDIAT-ULA 30
Las arcilla es la clase de partículas del suelo que exhibe la
mayor superficie del mismo, como se comprueba en la
segunda ecuación (los "discos" de arcilla pueden tener
espesores de 10-6 mm)
La fracción orgánica del suelo a menudo presenta superficies
altamente reactivas que exhiben las características de
materiales con una enorme superficie específica (hasta 1000
m2/g).
Por lo tanto la superficie específica real de un suelo está
muy influenciada por su contenido en materia orgánica y
arcilla.

Luis Rázuri Ramírez

CIDIAT-ULA 31
 Estructura del suelo

Las partículas del suelo, en especial las de mayor superficie

específica, con ayuda de la materia orgánica del suelo pueden
asociarse para formar grupos que presentan propiedades que
difieren de las de las partículas individualmente consideradas.

Esta propiedad del suelo se denomina estructura y es fundamental

en la interpretación de las propiedades macroscópicas que exhiben
los suelos .

Luis Rázuri Ramírez

CIDIAT-ULA 32
 Estructura del suelo

Cuando un suelo está encharcado, el tamaño de los espacios

entre las partículas del suelo (poros) está relacionados con
su capacidad de transmisión del agua.
Esto significaría que partículas pequeñas (tipo arcilla)
tendrían poros pequeños y por lo tanto la transmisión del
agua sería lenta.
Sin embargo en ocasiones se encuentran suelos arcillosos
que presentan una aceptable permeabilidad, lo que
normalmente está relacionado con el buen estado de
agregación o estructura del suelo bien desarrollada.

Luis Rázuri Ramírez

CIDIAT-ULA 33
 Relaciones
masa-volumen
 Relaciones masa-volumen del suelo

Volúmenes (cm3):
M
MA
VA = volumen de aire AR
AIRE ar

VH = volumen de agua
VS = volumen de sólidos V
VA MH Ma bH
g
VT = volumen total del suelo ar MMT
t Z
z

Masas (gr): V
VaH SOLUÇÃO
AGUA
g
VVT M
MS
MA = peso de aire t
s

MH = peso de agua
VVS SÓLIDOS
SÓLIDOS
MS = peso de sólidos
s

MT = peso total del suelo

Z = altura del estrato del suelo (mm);

H = altura de la lámina de agua en el suelo (mm).
 Masa específica de las partículas sólidas

Masa específica, también conocida como densidad real o

densidad de las partículas del suelo (rr), expresado en
gr/cm3, representa la relación entre la masa de las
partículas del suelo seco (MS), en relación a su volumen
(VS).

MS
rr 
VS

Valores típicos en suelos

minerales:

rr = 2,60 a 2,68 gr/cm3

 Masa específica de las partículas sólidas

Minerales Densidad real (gr/cm3)

Caolinita 2,60 – 2,68

Hematita 4,90 – 5,30

Magnetita 5,60 - 6,00

Mica 2,70 - 3,30

Cuarzo 2,65 - 2,66

 Masa específica de las partículas
sólidas

En laboratorio, la determinación de la densidad

real puede ser hecha, evaluando el volumen de
alcohol necesario para completar todos los
espacios vacíos de una muestra de suelo seco
al aire y que ha sido pasado por un tamiz con
una malla menor o igual a 2 mm.
 Masa específica del suelo o densidad
aparente
También conocida por densidad aparente o densidad del suelo
(ra), expreso en gr/cm3, representa la relación entre la masa de
las partículas del suelo seco (MS) más la masa del aire (MA), en
relación al volumen ocupado por las partículas sólidas (VS) más
el volumen ocupado por el espacio poroso (VA).

MS  MA MS
ra  
VS  VA VT

Densidad
Tipos de suelo
suelo (gr/cm3)
Arenosos 1,25 – 1,50
Arcillosos y francos 1,00 – 1,40
Húmicos 0,75 – 1,00
Turbosos 0,20 – 0,45
 Determinación de la densidad del suelo
o densidad aparente

Utilice un cilindro
metálico con diámetro
conocido;
En la foto del lado, el
cilindro metálico tiene
una altura de 5,1 cm y un
volumen de 100 cm3.
 Determinación de la densidad del suelo o
densidad aparente

Utilice un barreno
propio para sacar
muestras de suelo
con estructura no
disturbada.
 Determinación de la densidad del suelo
o densidad aparente

Para sacar la muestra de

suelo con la estructura no
disturbada, puede insertar
el propio cilindro
directamente.
En caso de que el suelo se
encuentre seco, es
necesario humedecerlo
para facilitar la penetración
del cilindro afectando lo
menos posible la
estructura natural. cilindro metálico
 Determinación de la densidad del suelo
o densidad aparente

Cuando se desee hacer el muestreo en profundidad,

deberán hacerse unos hoyos con un diámetro más
ancho que el muestreador para sacar la muestra.
 Determinación de la densidad del suelo
o densidad aparente

El exceso del suelo debe ser cuidadosamente

eliminado, utilizando un instrumento cortante,
de modo que el volumen final colectado sea
exactamente igual al volumen del cilindro.
 Determinación de la densidad del
suelo o densidad aparente

Las muestras de suelo

deberán ser colocadas
para secar en una estufa a
una temperatura
controlada de 105 ºC.
Déjela secando por 48
horas para garantizar que
toda el agua contenida en
el espacio poroso sea
complemente vaporizada.
 Determinación de la densidad del
suelo o densidad aparente
 Después que las muestras de suelo se encuentran
completamente secas, entonces hay que tomar su
peso final para efectuar el cálculo de la densidad
aparente.
 Ejemplo de cálculo de la densidad del
suelo o densidad aparente

Ejemplo:
Suponga que una muestra de suelo seco, arenoso,
después de permanecer durante 48 horas en la estufa
con temperatura controlada de 105 ºC, pesa 233,2 gr y
que el cilindro vacío pesa 81,0 gr.
¿Cuál es la densidad aparente de este suelo?
Solución:
MS
ra 
VS
233,2  81,0
ra   1,52 gr / cm 3
100
 Ejemplo de cálculo de densidad del
suelo o densidad aparente

Ejemplo:
Suponga que una muestra no disturbada de un suelo
arcilloso contiene solamente 3% de materia orgánica;
después de permanecer durante 48 horas en la estufa con
temperatura controlada de 105 ºC, tiene un peso de 180,2
gr. El peso del cilindro vacío es de 81,0 gr.
¿Cuál es la densidad aparente de este suelo?

Solución:
MS
ra 
VT
180,2  81,0
ra   0,99 gr / cm 3
100
 Porosidad total del suelo

Es todo el espacio vacío disponible para almacenamiento

del agua o del aire. Este parámetro puede ser estimado o
calculado, utilizando las informaciones sobre la densidad
real y la densidad aparente.
VA VT  VS VS
p  1
VT VT VT
VS
1
MS rr
p 1  1
VT 1
ra
MS
ra
p 1
rr
 Ejemplo de porosidad total del suelo

Supóngase que un suelo tiene una densidad aparente de

1,40 gr/cm3 y una densidad real de 2,65 gr/cm3. Calcule su
espacio poroso total.
r
p  1  a
rr
1,40
 1
2,65
 0,472 cm3 / cm3
Esto significa que hay un volumen potencialmente
disponible para el almacenamiento de agua y aire de
47,2%, o sea, cada 100 cm3 de suelo puede almacenar
47,2 cm3 de agua y aire;
En la práctica, ese es el espacio poroso disponible
Medidas del contenido de
agua en el suelo
 Medidas del contenido de agua en el suelo

La cantidad de agua
presente en el suelo
(MT o VT) puede ser
expresa como una
fracción del volumen
total (VT) o de la masa
total (MS) del suelo
seco.
 Humedad gravimétrica
o en base a peso seco

Es definido como la cantidad

total del agua (masa o volumen)
presente en el suelo y se
expresa en base al peso total
del suelo seco.

MH El componente agua:
Pw  VH  MH / rH
MS rH  1,0 g / cm3
VH  MH
 Muestreo del suelo para determinar la
humedad gravimétrica

El muestreo del suelo

para la determinación
de la humedad
gravimétrica puede
ser hecho en
estructura disturbada,
utilizando un barreno.
 Pasos para la evaluación de la
humedad gravimétrica

 La muestra colectada en el campo debe ser

almacenada en un recipiente o en una bolsa plástica y
mantenida herméticamente cerrada para prevenir
cualquier pérdida de agua, hasta la toma del peso
húmedo.
 Normalmente, la toma del peso húmedo es realizada
conjuntamente con el recipiente utilizado en el
transporte del suelo hasta el laboratorio, de modo que
también es necesario conocer el peso del recipiente.
 En seguida, la muestra es puesta a secar hasta que
toda el agua se haya vaporizado y entonces se toma el
peso seco.
 Ejemplo de cálculo de la humedad
gravimétrica
Ejemplo:
Suponga que una muestra de suelo fue retirada para determinación
de su humedad gravimétrica.
Los siguientes datos fueran obtenidos:
a) Peso del suelo húmedo: 87,21 gr.
b) Peso del recipiente utilizado: 27,11 gr.
c) Peso del suelo seco: 79,82 gr.
¿Cuál es la humedad gravimétrica en esta muestra de suelo?

Solución:
MH
Pw  ; M H  87,21  79,82  7,39 gr; M S  79,82  27,11  52,71 gr
MS
7,39
Pw   0,1402 gr / gr  14,02%
52,71
 Humedad volumétrica

Es definido como la cantidad

total de agua (volumen)
presente en el suelo y se
expresa en base al volumen
total del suelo.

VH
V 
VT El componente agua:
(cm3 de agua/cm3 de suelo) VH  M H / r H
r H  1,0 gr / cm 3
VH  M H
 Relación entre humedad volumétrica
y humedad gravimétrica

VH M H VH
V  Pw  
VT MS MS

MS
ra  MS  ra VT
VT

VH VH
Pw   Pw 
MS r a VT
VH
 r a Pw  V  r a Pw
VT
 Ejemplo de cálculo de la humedad
volumétrica
Ejemplo:
Suponga que una muestra de suelo fue retirada para determinación
de su humedad volumétrica.
Los siguientes datos fueran obtenidos:
a) Humedad gravimétrica: 14,02%
b) Densidad aparente del suelo: 1,41 gr/cm3

¿Cuál es la humedad volumétrica en esta muestra de suelo?

Solución:
 V  r a g
 V  1,41 14,02 
 V  19,77%
60
61
 Humedad volumétrica en la saturación

En principio, la humedad del

suelo en el punto de saturación
es siempre menor que su
porosidad total, porque por
mejor que sea el proceso de
saturación habrá algún espacio
del vacío del suelo que se
queda ocupado con el aire
S  p
entrapado y que no fue
ra
desalojado por el agua. S  1
rr
 Procedimientos para obtención de s

 Ponga las muestras del suelo para saturar en

una bandeja, durante por lo menos 24 horas,
manteniendo una capa de agua suficiente para
cubrir hasta dos tercios de la altura del cilindro;
 Quite, cuidadosamente, con una espátula curva,
la muestra de suelo para ser pesada, retirando
también el papel filtrante;
 Usando una espátula recta, retire todo suelo que
ha sido retenido en el papel filtrante y
recolóquelo de vuelta en el cilindro antes que
sea efectuado el pesaje;
 Las muestras de suelo deberán ser colocadas
para secar en una estufa con temperatura
controlada en 105º C, en cerca de 24-48 horas
para vaporización completa del agua contenida
en el suelo.
 Ejemplo de cálculo de S y p

Ejemplo:
Suponga que una determinada muestra de suelo fue
retirada para determinación de su humedad saturada y
de su porosidad total.
Los siguientes datos fueran obtenidos:
a) Peso del cilindro metálico: 82,9 gr
b) Peso del suelo húmedo más cilindro: 239,38 gr
c) Peso del suelo seco más cilindro: 183,04 gr
d) Densidad real: 2,76 gr/cm3

Calcule la humedad saturada, gravimétrica y

volumétrica, la porosidad total y la densidad aparente
 Ejemplo de cálculo de S y p

Datos: Pc  82,9 gr
M S  M H  Pc  239,38 gr
M S  Pc  183,04 gr
r r  2,76 gr / cm 3
Solución:
M
Pw  H ra 
MS  s  r a Pw
MS VS  s  1,000,5626
239,38  183,04 183,04  82,9
Pw  ra   s  0,5626 cm 3 / cm 3
183,04  82,9 100
Pw  0,5626 gr / gr r a  1,00 gr / cm 3

ra 1,00
p 1  1  0,6377 cm / cm
3 3

rr 2,76
 Relaciones de masa y volumen en un suelo

Las distintas relaciones de masa y volumen entre las tres

fases del suelo dan lugar a una serie de definiciones que han
sido consideradas útiles para caracterizar el estado físico de
un suelo que se resumen en la Tabla.

Luis Rázuri Ramírez

CIDIAT-ULA 66
Parámetro Símbolo Definición

Luis Rázuri Ramírez

CIDIAT-ULA 67
1) Densidad específica (ρs) Expresa la relación entre la masa de la
fracción sólida del suelo y el volumen que ésta ocupa en el sistema.

Para la mayoría de los suelos minerales ρs varia entre 2,6 a 2,7 g/cm3.
La presencia de materia orgánica disminuye el valor. La densidad de los
suelos se expresa con frecuencia en relación a la densidad del agua a
4°C (1 g/cm3) con lo que se eliminan las unidades (valor adimensional).

2) Densidad aparente (ρa ) Expresa la relación entre la masa de la

fracción sólida del suelo y el volumen total del suelo. Para la mayoría
de los suelos minerales ρa varia de 1,1 g/cm3 (suelos arcillosos) a 1,6
g/cm3 (suelos arenosos).

La ρa varia además en función de la estructura del suelo (grado de

compactación del suelo, presencia de fisuras, o grado de expansión o
contracción la arcillas), pero siempre tiene un valor inferior al de la ρs
debido a la presencia de poros en el suelo.
Luis Rázuri Ramírez
CIDIAT-ULA 68
3) Porosidad (p). El volumen de poros (o vacío) de un suelo está
constituido por el volumen de agua (Va) más el volumen de aire (Vg) en
un momento determinado (Figura).

La porosidad se define como el volumen de poros en una unidad de

volumen de suelo inalterado. La porosidad de los suelos minerales
varia desde 25% para suelos compactados hasta 70% para suelos no
consolidados. La Tabla siguiente resume los rangos de porosidad
normales en suelos de diferente textura.

4) Relación de vacío (e) se define como la relación entre el volumen

de poros (o vacío) y el volumen de sólidos en una unidad de volumen de
suelo. Esta relación se emplea normalmente en mecánica de suelos,
como medida de la compactación de un suelo (el denominador de la
expresión permanece constante al compactar una unidad de volumen
de suelo). Su valor varia desde 0,7 a 1,5 en arcillas bien
estructuradas.
Luis Rázuri Ramírez
CIDIAT-ULA 69
5) Contenido de agua (Ө y Pw). Puede ser expresado como la
relación entre el volumen de agua y volumen total del suelo
(contenido volumétrico de humedad, (Ө) o la relación entre la masa
de agua del suelo y la masa seca del suelo (contenido gravimétrico
de humedad, Pw).
Este concepto de contenido de agua es muy importante en el
estudio del agua en el suelo que se tratará posteriormente.

Luis Rázuri Ramírez

CIDIAT-ULA 70
Valores representativos de las propiedades físicas del suelo según texturas
(Schwab, 1996)

Cap. campo [Link]

Textura p(%) ρa (gr/cm )
3
θ (333) θ (15000)
(cm3/cm3) (cm3/cm3)
A 32-42 1.65 0.15 0.08
F-A 40-47 1.50 0.21 0.12
F 43-49 1.40 0.31 0.17
a-F 47-51 1.35 0.36 0.18
L-a 49-53 1.30 0.40 0.20
a 51-55 1.25 0.44 0.21

Luis Rázuri Ramírez

CIDIAT-ULA 71
 POTENCIAL DEL AGUA EN EL SUELO

Concepto de succión de agua en el suelo

Bajo la superficie del suelo las raíces de las plantas

comparten el espacio con las tres fases del suelo (partículas
sólidas, agua y aire).

En este medio el agua no se mueve libremente como lo hace

en la superficie sino que está más o menos retenida en la
matriz del suelo. Esta matriz condiciona como el agua se
mueve y como las raíces acceden a ella.

El agua es naturalmente atraída por las partículas del suelo

y se fija a la superficie de éstas y dentro de los espacios
capilares o poros de diverso tamaño entre las partículas
del suelo.
Luis Rázuri Ramírez
CIDIAT-ULA 72
 Concepto de succión de agua en el suelo

Cuando el suelo está muy húmedo, la mayoría de los poros

grandes están llenos de agua, el agua se mueve con bastante
libertad en el suelo y las raíces pueden absorberla fácilmente.
Cuando el suelo se seca, el agua que queda en el suelo es
atrapada con más fuerza en los capilares de menor tamaño del
suelo.

Luis Rázuri Ramírez

CIDIAT-ULA 73
 Concepto de succión de agua en el suelo

La Figura ilustra cómo la fuerza que se requiere para extraer el agua

de los poros de menor tamaño es cada vez mayor a medida que el suelo
se seca.

A causa de ésto las plantas encuentran cada vez más difícil obtener
suficiente agua cuando el suelo se seca.
Luis Rázuri Ramírez
CIDIAT-ULA 74
 La tensión del agua en el suelo

La matriz del suelo, por intermedio de

su intricado sistema de poros, ejerce
una fuerza de retención del agua, que
limita, en mayor o menor grado el
movimiento del agua en el suelo
 Fuerzas de retención del agua

ADHESIÓN: Fuerzas de atracción entre las moléculas de agua y

partículas de suelo; electrostáticas
COHESIÓN: Fuerzas de atracción entre moléculas de agua

Expresión de la energía de retención

Las plantas ejercen cierta fuerza por unidad de área de suelo para
absorber agua (Presión= fuerza por unidad de área). La unidad para
expresar presión es el bar o cb:
1 bar= 106 dinas/cm2
Pascal en el SIU:
1Pa = 1 Newton/m2 (Newton = kg/m/s2)
1Pa = 10-5 bares
1Mpa = 106 Pa = 10 bares
76
Fuerzas que actúan en la retención del
agua en el suelo
 Fuerza capilar en el suelo

2
T
R
En que,
T = tensión (Pa);
 = tensión superficial entre el
aire y agua [0,072 N/m];
R = radio del capilar [m]

Ejemplo: Calcule la tensión ejercida en el agua por

un capilar de radio igual a 10 micrones.

2 2 0,072 
T   14.400 Pa  14,4 kPa
R 10 / 1000000
 Concepto de succión de agua en el suelo

Cuando el agua se aloja sólo en capilares extremadamente

pequeños las plantas no pueden extraerla, se marchitan y
mueren.
Esto ocurre incluso en casos en los que en términos absolutos
todavía hay una cantidad considerable de agua en el suelo, pero
las plantas no pueden extraerla (típico de suelos muy
arcillosos).
Se define como succión del suelo, ψ, a la energía que hay que
emplear para desplazar el agua del suelo una distancia dada y
que a su vez determina en que dirección se mueve.

Luis Rázuri Ramírez

CIDIAT-ULA 79
 Concepto de succión de agua en el suelo

Este concepto es empleado cuando el suelo está por debajo

de saturación (todos los poros llenos de agua) con lo que la
succión tiene valores de presión negativa o tensión.

El concepto de succión corresponde al de potencial de

presión o mátrico del suelo (con el signo cambiado) como se
definirá a continuación.

Luis Rázuri Ramírez

CIDIAT-ULA 80
 Concepto de potencial de agua

El potencial del agua se refiere a su estado energético.

Hay varias fuerzas que actúan sobre el agua en el suelo

y que por lo tanto afectan su potencial. Entre ellas las
más importantes son: a) las fuerzas derivadas de la
interacción del agua con las partículas sólidas, b) las
fuerzas de interacción agua-moléculas en solución, c)
fuerzas del campo gravitatorio.

Excepto en suelos muy húmedos o salinos, las mayores

cuantitativamente son las de interacción con los
sólidos del suelo.

Luis Rázuri Ramírez

CIDIAT-ULA 81
 Concepto de potencial de agua

El potencial del agua tiene una definición termodinámica

que no emplearemos explícitamente para no complicar
innecesariamente el tema. Podemos definir el potencial del
agua como «aquella capacidad de hacer trabajo por
unidad de masa para mover esa masa de agua sin
solutos desde su estado hasta el de agua libre».

Luis Rázuri Ramírez

CIDIAT-ULA 82
 Concepto de potencial de agua

donde ψ se refiere a potencial y los subíndices t, g, p y o

indican: total, gravitatorio, presión y osmótico,
respectivamente.

Para simplificar, se incluye en el potencial de presión los dos

casos: el potencial mátrico (que se da en condiciones de no
saturación y que corresponde a una tensión o presión
negativa) y el hidrostático o de presión propiamente (que se
da en condiciones de saturación).
Luis Rázuri Ramírez
CIDIAT-ULA 83
 Unidades del potencial del agua

Según la definición de potencial las unidades son de

energía/masa, es decir : Joule/Kg. (J=N.m)

Sin embargo, debido a que la densidad del agua es la unidad, es

más corriente emplear las unidades de energía/volumen, que
numéricamente son las mismas que las anteriores (excepto por
el factor mil, ya que 1 m3= 103 Kg).

La ventaja de operar con joule/m3 = Newton/m2 = Pascal, que

son unidades de presión, conceptualmente más cómodas de
emplear.

Otras veces se emplean las unidades energía/peso que

equivalen a una longitud (altura de una columna de agua). En
estas unidades se tiene que 1 bar = 105 Pa = 10 m de
columna de agua (1000 cm H2O= 75 cm Hg).
Luis Rázuri Ramírez
CIDIAT-ULA 84
 Componentes del potencial total de agua

Potencial gravitatorio: es la que tiene el agua en virtud de su

posición sobre un plano de referencia dado, derivada de la
fuerza que la gravedad ejerce sobre el agua. Por tanto, la
altura a que se encuentre el agua afecta a su energía. La
energía potencial gravitatoria de una masa m a una altura z
(sobre el nivel de referencia) es:
E m g z  rV g z
V es el volumen y g es la aceleración debida a la gravedad.
Así el potencial será:
ψg = E /m = gz, expresado como energía por unidad de masa
ψg = E /V =ρ gz, expresado como energía por unidad de volumen
ψg = E /mg = z, expresado como energía por unidad de peso

Luis Rázuri Ramírez

CIDIAT-ULA 85
 Componentes del potencial total de agua
2) Potencial de presión: en el caso más frecuente de suelo no
saturado (poros parcialmente llenos por agua y aire) este
potencial es debido a las fuerzas de adhesión de las moléculas
de agua a las partículas sólidas del suelo y a las fuerzas de
cohesión entre las moléculas de agua lo que se conoce como
capilaridad.

Cuando el suelo se encuentre insaturado, el agua está bajo

tensión (o presión negativa) y se denomina también potencial
mátrico.

Si el suelo está saturado, entonces el agua en cada punto

soporta la presión hidrostática que depende de la altura de la
columna de agua que gravita sobre el punto. En ambos casos, la
manera más corriente de expresar este potencial es en
unidades de presión o altura Luis
deRázuri
[Link]írez
CIDIAT-ULA 86
 Componentes del potencial total de agua

La tensión mátrica en el suelo no saturado es debida al

principio de capilaridad ilustrado en la Figura. La capilaridad se
manifiesta como la elevación del agua en el interior de un
capilar (tubo de poco diámetro) o la retención de parte de agua
en su interior cuando se drena el tubo.

Luis Rázuri Ramírez

CIDIAT-ULA 87
 Componentes del potencial total de agua

Cuando se introduce un capilar de radio r en el agua se origina

un tensión en el interior del poro cilíndrico como consecuencia
de la formación de una superficie cóncava de radio de
curvatura R = r / cos α cuyo valor es:

P  2 / R  2 (cos  ) / r

donde ∆P significa diferencia de presión entre el aire en

contacto con la superficie cóncava y el agua, σ es la tensión
superficial del agua (0,0728 N/m a 20 ºC) y α es el ángulo de
contacto del agua con la superficie del sólido; si el agua moja
perfectamente el sólido α = 0 (p.e. agua pura-vidrio).

Luis Rázuri Ramírez

CIDIAT-ULA 88
 Componentes del potencial total de agua
Dado que la presión hidrostática es proporcional a la
profundidad bajo la superficie libre del líquido (P = ρghc),
podemos inferir que la tensión hidrostática (presión negativa)
en un punto en el interior de un tubo capilar al nivel de la
superficie libre del líquido en el exterior (punto A en Figura)
es proporcional a la altura hc sobre la superficie libre del
agua,
hc  2 (cos ) / rgr)

Para un capilar de vidrio con

agua pura a 20°C, la ecuación
anterior se puede simplificar a
hc 1.48xl0 -5 /r si se utiliza
unidades del SI.
Luis Rázuri Ramírez
CIDIAT-ULA 89
Componentes del potencial total de agua

Luis Rázuri Ramírez

CIDIAT-ULA 90
Componentes del potencial total de agua

Luis Rázuri Ramírez

CIDIAT-ULA 91
 El agua del suelo en condiciones hidrostáticas
En condiciones hidrostáticas (en
el equilibrio y sin movimiento de
agua) todos los puntos del perfil
de suelo tienen el mismo valor
de H, y los potenciales se
distribuyen según la Figura.
En la figura se representa la
superficie libre del líquido (SL)
en la superficie del suelo, por lo
tanto los poros de la columna
están totalmente llenos de agua
y se dice que el suelo está
saturado.
El potencial de presión en el punto A es igual a la distancia entre el
punto y el nivel freático con signo positivo en saturación (por debajo
del SL) y negativo en insaturación (por encima del SL), hA en la Figura
Luis Rázuri Ramírez
CIDIAT-ULA 92
También es igual a la altura a la que el agua ascenderá
(descenderá en caso de suelo por debajo de saturación) sobre
ese punto en un tubo abierto con su extremo insertado en ese
punto. Se llama piezómetro a un tubo de dichas características.
El potencial gravitatorio puede definirse como la distancia de
ese punto al plano de referencia con el signo correspondiente
(ZA en la Figura).
Se denomina nivel freático a la línea imaginaria que uniría los
puntos con potencial de presión igual a 0.
En condiciones hidrostáticas la elevación del nivel freático
coincidirá con la altura de agua en el piezómetro.

Luis Rázuri Ramírez

CIDIAT-ULA 93
Consideremos ahora una columna
de suelo como el de la Figura, en
la que el nivel freático se
encuentra por debajo de la
superficie, a una altura sobre el
plano de referencia z=L. Los
poros por encima del nivel
freático contienen agua y aire.
Se dice que el suelo está
insaturado.
El potencial de presión en un punto por encima del nivel
freático es negativo (potencial matricial). Dicho potencial
expresado en términos de energía por unidad de peso es
precisamente la succión ψ,
  h
Luis Rázuri Ramírez
CIDIAT-ULA 94
El punto B por encima del nivel
freático (Figura) se encuentra
en la zona insaturada y su
potencial de presión es
negativo (tensión). Si el suelo
está en condiciones estáticas, el
potencial hidráulico es, como en
el caso anterior, constante e
igual a la distancia desde el plano
de referencia al nivel freático,
L.
H  z  h  z   cte
En el suelo insaturado la medida del potencial de presión se
efectúa mediante tensiómetros
Luis Rázuri Ramírez
CIDIAT-ULA 95
 ¿Cómo medir la tensión del agua en el
suelo?
Baja tensión

Cápsula porosa

Suelo T = - 10 cm

Tubo lleno de agua

 ¿Cómo medir la tensión del agua en el
suelo?

Muestra de
suelo

Plato
poroso

Columna de agua
pendiente

Tubo lleno de
agua
Baja tensión
 ¿Cómo medir la tensión del agua en el
suelo?

Media tensión 13,65L  L  h  P  T

T  12,65L  h  P (cm H 2 0)
 0,098 12,65L  h  P  (kPa)

L (altura de mercurio)

h (altura del recipiente de mercurio)

Nivel de
referencia

P (profundidad del tensiómetro)

T (tensión del agua en el suelo)

 ¿Cómo medir la tensión del agua en el
suelo?

El agua es colectada
Presión del aire aplicado en la bureta hasta que
el equilibrio es
Fuente de alcanzado
presión

Media tensión

Trampa para burbujas de aire

 Curva característica de humedad del suelo
(o de succión)

Se denominan curva característica del suelo a la curva que

relacionan el contenido de agua con la succión, ψ. En general,
estas curvas se obtienen en el laboratorio utilizando
muestras inalteradas de suelo y un equipo de placas porosas
de presión.

El contenido de agua de un suelo a bajas succiones es función

principalmente de su estructura, mientras que a altas
succiones el factor determinante es la textura (esto es
debido a que la estructura determina en gran medida la
proporción de poros grandes, mientras que la textura
determina la proporción de poros pequeños).

Luis Rázuri Ramírez

CIDIAT-ULA 100
 Curva característica de humedad del suelo
(o de succión)
En la Figura se presentan curvas de este tipo para dos casos
obtenidos en Tenerife. Una de las curvas representa un suelo
arcilloso donde el contenido de humedad permanece alto incluso
a altas succiones. Por el contrario, el arenoso-gravoso, es un
material suelto de origen volcánico, cuyo contenido de humedad
desciende muy rápidamente a valores muy bajos de succión (<10
cm H2O).

Luis Rázuri Ramírez

CIDIAT-ULA 101
 Capacidad de campo

Según la definición de Veihmeyer y Hendrickson (1949) es " la

cantidad de agua retenida en el suelo después que el agua en
exceso ha drenado y la velocidad de flujo a capas profundas
ha cesado prácticamente, lo cual ocurre en los 2-3 días
después de una lluvia o riego en suelos permeables de
estructura y textura uniforme"

El principal problema de esta definición es que, en general, la

redistribución no cesa totalmente con el tiempo. Tradicionalmente
se ha considerado que es constante para cada suelo. Se ha usado
este concepto para determinar la dosis de riego. Este concepto es
de mayor utilidad en suelos arenosos con buen drenaje.

Algunos métodos de laboratorio para determinar la humedad a

capacidad de campo son:
a) centrifugación a l000g.
b) succión de 1/3 bar mediante una placa de Richards.
Luis Rázuri Ramírez
CIDIAT-ULA 102
Los principales factores que afectan a la capacidad de campo
son:

• Textura del suelo

• Cantidad y tipo de arcilla: montmorillonita aumenta el agua
retenida.
• Profundidad de frente húmedo y humedad inicial.
• Presencia de estratos de diferente textura
• Evapotranspiración.

Otro factor que afecta a la capacidad de campo (y al punto de

marchitez) es el contenido de materia orgánica del suelo.

Aunque en general las diferencias en contenido de humedad del

suelo en función de su materia orgánica no son muy grandes,
estas se hacen más notables cuanto menor es el contenido de
arcilla del suelo, lo que afectaLuisaRázuri
su Ramírez
comportamiento hidrológico.
CIDIAT-ULA 103
 Punto de marchitez permanente

El punto de marchitamiento permanente es el punto de

humedad mínima en el cual una planta no puede seguir
extrayendo agua del suelo y no puede recuperarse de la
pérdida hídrica aunque la humedad ambiental sea
saturada

Es la tensión máxima que puede realizar un cultivo para

extraer el agua del suelo. A partir de allí, esa planta en
esas condiciones de humedad no tendrá posibilidades de
abastecerse de agua

Luis Rázuri Ramírez

CIDIAT-ULA 104
En capítulo anterior se incluyó un estudio del contenido de humedad
para el caso de un suelo drenado al equilibrio que interesa contrastar
con esta definición de capacidad de campo. ¿Es la capacidad de campo
una propiedad de todo el perfil del suelo cuando éste está drenado al
equilibrio? (ver figura).

Luis Rázuri Ramírez

CIDIAT-ULA 105
• Constantes de humedad del suelo:
• Capacidad de Campo:
-30 KPa = -30 cb = -0.03 Mpa

• Punto de Marchitez permanente:

-1500KPa = -1500 cb = -1.5 MPa

• Agua “disponible” para la planta:

(CC – PMP)

• Agua “fácilmente disponible” para la planta:

[(CC-PMP)/2]

-30 a -70 cb = -30 a -70 KPa = -0.03 a -0.05 MPa

Luis Rázuri Ramírez

CIDIAT-ULA 106
 El agua en el suelo
SUELO ARCILLOSO

53 % SATURACIÓN

CAPACIDAD
35 %
DE CAMPO

17 % P. M. P.

DRENAJE

AGUA = GRAVITACIONAL DISPONIBLE

Luis Rázuri Ramírez NO DISPONIBLE
CIDIAT-ULA 107
Luis Rázuri Ramírez
CIDIAT-ULA 108
Luis Rázuri Ramírez
CIDIAT-ULA 109
 Agua
Aire Aire Aire

Partícula
de suelo

Porcentaje de Capacidad de Porcentaje de

Saturación campo marchitamiento

1/10 atmósfera 1/3 atmósfera 15 atmósferas

Luis Rázuri Ramírez

CIDIAT-ULA 110
Luis Rázuri Ramírez
CIDIAT-ULA 111
 -1500
CURVA RETENCION
HUMEDAD

Ψm
(kPa)
-100

-
33

0 sat
50 75 100
HumedadLuisaprovechable
Rázuri Ramírez
(%)
CIDIAT-ULA 112
Curvas de retención/desabsorción de agua
Muestran la relación entre el contenido de humedad y la tensión del
agua del suelo. La curva es característica de cada suelo, pues influyen
propiedades como textura (sup. específica), estructura, MO,
configuración del espacio poroso.

Es relevante la sup. específica de las arcillas, ya que la adsorción del

agua es un fenómeno superficial, así el área expuesta, la densidad de
carga y los cationes saturantes son de importancia. Ej.

Luis Rázuri Ramírez

CIDIAT-ULA 113
 70
Retención del agua del suelo

60

50 Agua
disponible
40
% Agua

Arcilla 30%
30

20
Arena 7%
10

0
0,033 1,5 MPa
Tensión
Luis Rázuri Ramírez
CIDIAT-ULA 114
Retención de agua en suelos de Costa Rica

Prof. % agua ap H2O

Suelo
cm 0,34 15,1 g/cm3 Útil
Alajuela plano 0-18 70,5 30,0 0,60 40,5
Alajuela ondul. 0-18 63,5 30,0 0,83 33,5
Grecia 0-18 50,0 29,5 0,90 20,5
Ciruelas 0-18 56,0 35,0 0,79 21,0
Paraíso 0-30 46,0 32,0 0,75 14,0
Cervantes 0-30 84,0 62,0 0,56 22,0
Birrisito 0-50 88,0 51,5 0,49 36,5
Instituto 0-30 42,0 27,0 1,82 15,0

Luis Rázuri Ramírez

CIDIAT-ULA 115
 Humedad Aprovechable en función de la textura

Capacidad de Marchitez Humedad

campo permanente aprovechable
Textura del suelo
(CC) (PMP) (CC-PMP)
Arenoso 9 4 5
Franco arenoso 14 6 8
Franco 22 10 12
Franco arcilloso 27 13 14
Arcillo arenoso 31 15 16
Arcilloso 35 17 18

Luis Rázuri Ramírez

CIDIAT-ULA 116
 Agua disponible en el
suelo
30
40

24

CC 30 Contenido
Agua agua
18

suelo Agua disponible cm/m suelo

20
(%) 12

PMP
6
10
Agua no disponible
0

Arena franco franco franco franco arcilloso

arenoso limoso arcilloso
Luis Rázuri Ramírez
CIDIAT-ULA 117
 MOVIMIENTO DEL AGUA EN EL SUELO

PMP CC PE

Agua Agua Agua

Higroscopica Capilar Gravitacional

MPa
3.1 -1.5 -0.03 0
No disponible Disponible
Luis Rázuri Ramírez
Exceso
CIDIAT-ULA 118
AGUA UTIL O HUMEDAD APROVECHABLE

D
D
Luis Rázuri Ramírez
CIDIAT-ULA
 MEDIDA DE LA HUMEDAD DEL SUELO

Existen en la actualidad unos 15 métodos para la medida del

agua del suelo que pueden ser clasificados en métodos
volumétricos (aquellos que dan estimaciones de la humedad
volumétrica del suelo) o métodos tensiómétricos (aquellos que
lo hacen de la tensión o potencial del suelo).

Se describirán los siguientes:

•Tensiómetros
•Bloques de resistencia eléctrica (Bouyoucos)
•Sonda de neutrones
•TDR (Time Domain Reflectometry)
Si se desea un buen estudio reciente de comparación de 23
sensores utilizados con frecuencia en agricultura, se
encontrará en el libro de Yoder et al (1998).
Luis Rázuri Ramírez
CIDIAT-ULA 120
 Tensiómetros (medida de h + z)

Principio: El tensiómetro es un tubo lleno de agua

desgasificada que insertado verticalmente en el suelo une
hidráulicamente el vacuómetro de su extremo superficial con el
agua del suelo a través de la cápsula de porcelana porosa
situada en su extremo inferior.

El agua del suelo transmite su tensión a la columna de agua del

tensiómetro, produciendo un desplazamiento del agua del
interior hacia el suelo.

El agua de la columna a su vez tira del mecanismo del

vacuómetro produciendo del desplazamiento mecánico de la
aguja del dial (vacuómetro de Boume), una depresión en la
columna del manómetro de mercurio, o un impulso eléctrico a
través de un transductor de presión (Figura).
Luis Rázuri Ramírez
CIDIAT-ULA 121
El Tensiómetro

Luis Rázuri Ramírez

CIDIAT-ULA 122
 Tensiómetro con medidor del tipo
vacuómetro

Depósito de agua

Vacuómetro

Cuerpo del
tensiómetro

Cápsula de
cerámica
 Tensiómetro con medidor con escala
de mercurio

Escala de
mercurio

Soporte de
sustentación

Cuerpo del
tensiómetro
Tensiómetro con medidor con escala
de mercurio
 Tensiómetro con medidor electrónico

Transductor
de presión

Cuerpo del
tensiómetro

Medidor
electrónico
portátil
 Tensiómetro con medidor electrónico
con aguja de perforación

Tapa de
goma con
silicona
Medidor
electrónico
portátil
Cuerpo del
tensiómetro
Tensiómetros (medida de h + z)

Luis Rázuri Ramírez

CIDIAT-ULA 128
 Tensiómetros (medida de h + z)

El vacuómetro del tensiómetro da una

lectura Hm, que no es igual al potencial de
presión en la cápsula (h) directamente al
haber una columna de líquido entre ambos
de longitud |z|+|ho|.
El líquido en el tubo entre la cápsula y el
vacuómetro está en equilibrio hidrostático
y por lo tanto la presión en el tubo
aumenta linealmente con la profundidad y
el potencial de presión en la cápsula del
tensiómetro es:

h  H m  z abs  ho abs

Luis Rázuri Ramírez

CIDIAT-ULA 129
 Tensiómetros (medida de h + z)
Cuando se añade agua al suelo la succión decrece, parte del agua entra
en la cápsula produciendo una subida del nivel del agua en el tubo que
se transmite al manómetro como una bajada de succión.

Nótese que como el líquido del tensiómetro está en equilibrio estático,

si se considera el plano de referencia en la superficie del suelo (z>O
hacia arriba) el potencial hidráulico en la cápsula sería,

H  h  z  H m  ho
Nótese que si se pudiera despreciar ho, se llega a H ≈Hm+O=Hm (el
potencial medido en el vacuómetro es igual al potencial hidráulico).
Nótese también que si tuviéramos dos tensiómetros A y B con los
vacuómetros a la misma altura sobre el suelo (hoA=hoB), la diferencia
de potencial hidráulico entre ellos seria la diferencia entre las
lecturas de sus vacuómetros:
H  H A  H B  H mA  H mB
Luis Rázuri Ramírez
CIDIAT-ULA 130
 Tensiómetros (medida de h + z)

Ejemplo: Consideremos un tensiómetro enterrado a una

profundidad de 0,6 m tiene el vacuómetro a 20 cm sobre el suelo
y da una lectura de succión = 9 cbar.

La lectura manométrica en términos de potencial es Hm=-0.9 m, el

potencial de presión en la capsula resulta entonces h = -0,9 + |-0,6|
+ |0,2| = -0,1 m, y el potencial hidráulico (considerando el plano de
referencia en la superficie y positivo hacia arriba) seria H = h + z =
-0.1 -0.6 = -0,7m.

Descripción: El tensiómetro consta de un tubo de plexiglás

normalmente transparente, con un extremo en el que se encuentra la
cápsula de porcelana porosa y en el otro un tapón y el mecanismo del
manómetro. El interior se llena con agua desgasificada hasta por
encima del punto de inserción del vacuómetro.
Luis Rázuri Ramírez
CIDIAT-ULA 131
 Tensiómetros (medida de h + z)
Uso: Previo a su instalación hay que proceder a saturar la cápsula y
eliminar el aire que pudiera quedar almacenado en el vacuómetro y las
paredes del tubo y cápsula. Esto se realiza llenando el tensiómetro con
agua desgasificada y sumergiendo la cápsula en agua haciendo vacío a la
vez varias veces con una bomba de succión.

Durante la instalación hay que asegurarse de que se produzca un

contacto íntimo entre la cápsula y el suelo circundante. Esto se suele
hacer perforando en el suelo con una barrena un agujero vertical de un
diámetro próximo al del tubo (5 cm suele ser suficiente). La tierra
extraída se tamiza y se mezcla con agua preparando un barro fino que
se coloca en el fondo del agujero donde se introduce la cápsula. El
resto del espacio se rellena con suelo normal.

En el cuello del tensiómetro, en la superficie se puede colocar un anillo

de bentonita u otro material impermeable para evitar que el agua de la
superficie se introduzca hasta la cápsula por las paredes exteriores de
tubo. Luis Rázuri Ramírez
CIDIAT-ULA 132
 Tensiómetros (medida de h + z)

Desventajas:
• Necesitan revisión y rellenado frecuente. Si el nivel del líquido cae
por debajo del punto de inserción del manómetro (lo que puede pasar
en periodos secos si no se revisan con frecuencia) hay que sacarlos
del suelo y empezar el proceso de instalación desde el principio.

• Los vacuómetros de Bourne no son muy precisos, pero este

inconveniente se puede superar con transductores de presión.

• La lectura es de potencial del suelo, no del contenido de humedad.

Para conocerlo hay que disponer de una curva de succión. El fenómeno
de histéresis en la curva de succión del suelo puede introducir
errores en las lecturas según se haga la lectura durante el secado o
mojado del suelo.

• Trabaja sólo en el rango útil de 0 a - 85 KPa.

Luis Rázuri Ramírez
CIDIAT-ULA 133
 Tensiómetros (medida de h + z)

Ventajas:

• Precio reducido.

• Es una medida directa y muy intuitiva del potencial del suelo

• Su lectura se pueden automatizar con transductores de presión y

un equipo registrador.

•Su instalación no resulta especialmente difícil.

Luis Rázuri Ramírez

CIDIAT-ULA 134
 Bloques de Resistencia Eléctrica
(medida de θp + θo)
Principio: Un bloque de yeso se equilibra (enterrándolo) con la solución
del suelo en condiciones en las cuales tanto agua como solutos son
intercambiados.
El potencial mátrico del bloque (igual al del suelo circundante en el
equilibrio) se calcula interpolando en una relación previamente calibrada
entre la impedancia eléctrica del bloque (función de su contenido de
agua) y el potencial hídrico del bloque.
La conductividad eléctrica del sensor está determinada principalmente
por su contenido de humedad hasta un límite máximo de salinidad de la
solución del suelo (2 mmho/cm). A partir de ese límite ambos factores
(contenido de humedad y contenido en sales) afectan las lecturas de la
impedancia del bloque.
Descripción: El aparato consta de las siguientes partes: i) sensor de
medida: formado por dos electrodos colocados en el interior de un
bloque de yeso a una separación determinada; ii) cable conector de
longitud variable; iii) terminal de lectura, en el que se distinguen la
escala de lectura, y botones de puesta en marcha y calibrado a cero.
Luis Rázuri Ramírez
CIDIAT-ULA 135
 Equipos para medidas en el suelo
BLOQUES DE YESO CON SENSORES PARA MEDICIONES
DE LA RESISTENCIA ELÉCTRICA EN EL SUELO
 Bloques de Resistencia Eléctrica
(medida de θp + θo)
Calibración: La utilización del aparato requiere la obtención de una
curva de calibración que recoge la relación de impedancia eléctrica
medida frente a una succión del suelo o su contenido de humedad.

Para su obtención existen dos procedimientos alternativos: i)

determinación gravimétrica (en laboratorio o campo); ii) calibración
frente a curvas de succión (en laboratorio según lo explicado)

Estos dos tipos de calibración arrojarán dos tipos de curvas. En general

se considera más precisa la calibración frente a succión en el suelo
usando placas de porcelana porosas.

Uso: Se humedece el bloque a saturación antes de su instalación en

campo. Se entierra el bloque a la profundidad de lectura deseada
procurando que el suelo quede bien compactado alrededor del bloque y
que se produce la mínima perturbación en el suelo. Pasado un tiempo en
el que se deja equilibrar con la solución del suelo, se pueden hacer
lecturas. La lectura se efectúa directamente sobre la escala y se
interpola el valor en la gráfica de calibración.
Luis Rázuri Ramírez
CIDIAT-ULA 137
 Bloques de Resistencia Eléctrica
(medida de θp + θo)
Desventajas o Críticas:
• La salinidad del suelo (en exceso de 2 mmhos/cm) introduce errores
en las lecturas.
• La calibración del aparato es tediosa.
• Los bloques de yeso envejecen en el suelo, especialmente en suelos
salinos y muy húmedos, produciéndose una degradación de las lecturas
hasta su desintegración.
• El fenómeno de histéresis en la curva de succión del suelo puede
introducir errores en las lecturas según se haga la lectura durante el
secado o mojado del suelo
• Es difícil especificar la exactitud de las medidas debido a las
muchas fuentes de error a la que están sometidas (calibración,
materiales, instalación, envejecimiento, histéresis), en general se
considera que pueden producirse errores de 2-100% en las lecturas.
• Temperaturas inferiores a 0°C tienen efecto sobre las lecturas
(incremento).
Luis Rázuri Ramírez
CIDIAT-ULA 138
 Bloques de Resistencia Eléctrica
(medida de θp + θo)

Ventajas:

• Precio reducido

• Medida in situ no destructiva (sin muestreo)

• No precisa operaciones de mantenimiento

Luis Rázuri Ramírez

CIDIAT-ULA 139
 Sonda de neutrones (medida de θ)
Principio: Se basa en el principio de termalización de los
neutrones de alta energía ("rápidos") emitidos por una fuente
radiactiva (como el americío-berilio) al golpear con las
moléculas de agua del medio.

Los neutrones rápidos chocan con los átomos de hidrógeno del

agua y forman lo que se denomina neutrones termalizados o
"lentos", de menor energía.

Las múltiples colisiones en un punto del suelo forman una nube

termalizada de neutrones de tamaño fijo pero de densidad
dependiente del número de colisiones, o lo que es igual del
número de átomos de hidrógeno (moléculas de agua) presentes.

Un contador de electrones "lentos", instalado en las

inmediaciones del emisor mide la nube. La medida se presenta
como la relación entre electrones rápidos (emitidos) y lentos
"recibidos". Luis Rázuri Ramírez
140
CIDIAT-ULA
 Sonda de neutrones (medida de θ)
Descripción: El aparato consta de un contador externo y de
una sonda emisor. Receptor de forma cilíndrica que se
introduce en un tubo metálico vertical instalado de forma
permanente en el suelo.

Calibración: Se realiza midiendo la relación de conteo a varias

profundidades en un suelo dado en diferentes momentos y
tomando simultáneamente muestras de suelo para determinar
su humedad gravimétrica mediante secado en estufa.

La relación obtenida entre conteo y contenido de humedad es

normalmente lineal y específica para cada suelo.

Uso: La sonda se desplaza a varias profundidades dentro del

tubo donde se realizan lecturas de conteo que se traducen a
humedad mediante la recta de calibración obtenida
anteriormente. Luis Rázuri Ramírez
141
CIDIAT-ULA
 Sonda de neutrones (medida de θ)
Crítica:
• Contiene una fuente radiactiva por lo que su uso implica un riesgo y
sólo está permitido a personal autorizado cualificado. Está sujeto a
inspección periódica por la Junta de Energía Nuclear con su
consiguiente protocolo.
• Precio elevado (en tomo a los $8.500-10.000).

• Por su naturaleza radiactiva no puede ser automatizado al no poder

dejarse la sonda sin atención en el campo.

• Requiere calibración.
Ventajas:
• Es muy resistente y se han descrito casos de 30 años de uso
continuado.
• La calibración sólo es necesaria una vez para cada suelo ya que no se
degrada el sistema de medida.
• Muy utilizado. Aunque en retroceso, aporta resultados muy
reproducibles.
Luis Rázuri Ramírez
CIDIAT-ULA 142
 Time Domain Refectometry (TDR)
(medida de θ)

Principio: El TDR es un método no destructivo e instantáneo

de medida de humedad en el suelo. El método tiene sus
orígenes en el desarrollo del RADAR en 1941 y en los años 80
fue puesto a punto para los suelos por Topp et al. (1984).

El método se basa en que existe una relación única entre el

contenido volumétrico de humedad (θ) y la constante
dieléctrica de los suelos (Ka).

La determinación del contenido de humedad por este método

está basado en el alto valor de Ka para el agua (Ka=80)
comparado con el pequeño valor para los sólidos del suelo
(Ka=5).

Luis Rázuri Ramírez

CIDIAT-ULA 143
 Time Domain Refectometry (TDR)
(medida de θ)

El valor de la constante dieléctrica global del suelo está por

tanto condicionada por el contenido de agua del mismo:
pequeñas variaciones en el contenido de humedad provocan
cambios importantes en el valor de esta propiedad. La
constante dieléctrica del medio puede determinarse por el
tiempo de desplazamiento, t (del orden de 100 ps a 20 ns),
de una onda electromagnética a través del material.

Luis Rázuri Ramírez

CIDIAT-ULA 144
 Time Domain Refectometry (TDR)
(medida de θ)

La onda electromagnética se genera en forma de un plano

mediante un pulso rápido en los extremos una sonda
compuesta por dos (ó 3) barras paralelas de metal de longitud
l que actúan como una línea transmisora, y la ecuación
resultante es:

K a  (c.t / 2l ) 2  (1.5x10 8 t / l ) 2

para l en metros, t en segundos y c velocidad de la luz (3x108

m/s)
El pulso electromagnético es reflejado al llegar al final de la
sonda y vuelve al extremo donde se generó inicialmente. El
tiempo de desplazamiento es el retraso entre la generación
del pulso y la lectura de vuelta en el receptor, afectado por la
constante dieléctrica y por consiguiente por el contenido de
humedad del suelo. Luis Rázuri Ramírez
CIDIAT-ULA 145
 Time Domain Refectometry (TDR)
(medida de θ)

Esto se ilustra en la Figura, donde puede verse como el pulso se

atenúa durante el transporte de la onda hasta su vuelta
después de recorrer la sonda de TDR (subida de la señal en
torno a t=6ns), y el tiempo de retraso sería 2t = t1-t2=2,38 ns.

Luis Rázuri Ramírez

CIDIAT-ULA 146
 Time Domain Refectometry (TDR)
(medida de θ)

Los trabajos de Topp et al. (1984) mostraron que para la

mayoría de los suelos minerales puede emplearse la ecuación
universal expresada como:

   5.3x102  2.92 x102 Ka  5.5x104 Ka 2  4.3x106 Ka3

Descripción: Existen varios tipos de aparatos basados en la

misma tecnología. Los primeros aparatos estaban basados en
los comprobadores de cable Tektronix 1502C a los que se
añadían la sondas de fabricación casera (dos barras de metal
sobre un soporte, cable y un transformador de señal (balun).

Luis Rázuri Ramírez

CIDIAT-ULA 147
 Time Domain Refectometry (TDR)
(medida de θ)

El Tektronix consta de una pantalla sobre la que se

representa una grafica del pulso enviado así como del
recibido. El tiempo de desplazamiento se obtiene
directamente de la gráfica y el contenido de humedad se
obtiene aplicando la fórmula universal.
Los nuevos sistemas (SoilMoisture, IMKO, etc) tienen sus
unidades de lectura especialmente desarrolladas para este
propósito y proporcionan medidas del contenido de humedad
directamente presentando además la capacidad de almacenar
lecturas, modo automático de lectura, etc.

Luis Rázuri Ramírez

CIDIAT-ULA 148
 Time Domain Refectometry (TDR)
(medida de θ)
Calibración: Para la mayoría de los suelos minerales se emplea la
ecuación universal por lo que no requiere calibración. Para valores
precisos se recomienda una calibración con muestra inalterada o in-
situ.

Uso: Las sondas se entierran a las distintas profundidades de

estudio. Normalmente se perfora el suelo usando un muestreador en
tubo, similar al empleado para la colocación de tensiómetros o
piezómetros, y se desciende la sonda hasta su fondo.

Se clavan las guías asegurando un buen contacto con el suelo y se

recubre la sonda con el suelo extraído dejando los cables sobre la
superficie protegidos por un cajetín de la acción del agua, animales,
etc.

En el momento de hacer las lecturas se conectan los cables a la

unidad de lectura y se efectúa la lectura sobre la pantalla de forma
instantánea tras la generación de la onda (según modelos).
Luis Rázuri Ramírez
CIDIAT-ULA 149
 Time Domain Refectometry (TDR)
(medida de θ)
Críticas:
• Precio, varía de $7.500 a $ 10.000 + según sistema y número de
sondas.
• Requiere instalación cuidadosa
Ventajas:
• Gran exactitud de medida (1 %)
• Medida in situ no destructiva (sin muestreo)
• Permite la automatización de las medidas
• Muy rápida
• Gran sencillez de uso tras su instalación
• No necesita calibración para la mayoría de los suelos minerales.
• El sistema está basado en componentes electrónicos sin piezas
móviles lo que elimina una fuente de desgaste y mantenimiento,
permitiendo medidas durante largos periodos de tiempo.
• La lectura de humedad es independiente de la temperatura,
densidad del suelo, salinidad y composición mineral del suelo.
Luis Rázuri Ramírez
CIDIAT-ULA 150
 Equipos para medidas en el suelo

MODELOS DE TDR DE LA COMPAÑIA DELTA-T

Luis Rázuri Ramírez
CIDIAT-ULA 152
Equipos para medidas en el suelo

OTROS MODELOS DE TDR

154
 Colocación de los aparatos de medida de
humedad

La colocación de los aparatos varía con el objetivo de la

medida que se pretende realizar, así como con el espesor de
suelo que se pretende controlar y con la distribución de
humedad en el perfil del suelo.

Como ejemplo, en riego localizado conviene disponer de más

de un punto de control de humedad, situados a diferentes
profundidades (Figura). Ello permitirá controlar la extensión
y la profundidad del bulbo para las dosis de riego aplicado a
lo largo del ciclo del cultivo y ajustar la dosis en periodos de
mínima/máxima demanda, de manera que se mantenga el
tamaño del bulbo deseado.

Luis Rázuri Ramírez

CIDIAT-ULA 155
Luis Rázuri Ramírez
CIDIAT-ULA 156
Localización de los tensiómetros
 FLUJO DEL AGUA EN EL SUELO

• EL CONTENIDO DE AGUA Y LA VELOCIDAD CON QUE ÉSTA

SE MUEVE EN EL INTERIOR DEL PERFIL DEL SUELO DEPENDE
EN GRAN MEDIDA DE LAS CARACTERISTICAS FÍSICAS DEL
SUELO.

Luis Rázuri Ramírez

CIDIAT-ULA 158
EL FLUJO DE AGUA EN EL SUELO
SE DESARROLLA
INCESANTEMENTE, CON
VELOCIDADES Y DIRECCIONES
VARIABLES

Luis Rázuri Ramírez

CIDIAT-ULA 159
PRINCIPALES PROCESOS DE FLUJO DEL
AGUA EN EL SUELO

* Infiltración
* Redistribución
* Drenaje o percolación
* Evaporación
* Movimiento hacia las raíces

Luis Rázuri Ramírez

CIDIAT-ULA 160
 MOVIMIENTO DE AGUA EN EL SUELO

• CLASIFICACIÓN DEL AGUA PRESENTE

EN EL SUELO:

– AGUA DE HIDRACIÓN (NO DISPONIBLE)

– AGUA HIGROSCÓPICA (NO DISPONIBLE)
– AGUA CAPILAR (DISPONIBLE)
– AGUA GRAVITACIONAL (DISPONIBLE POR POCO
TIEMPO)

Luis Rázuri Ramírez

CIDIAT-ULA 161
Luis Rázuri Ramírez
CIDIAT-ULA 162
 REDISTRIBUCIÓN

• Desde otra parte del suelo con potencial

mayor (menos negativo).
• Liberación de agua desde el sistema de
raíces.
• Entre sistemas de raíces en asociación de
cultivos.
• Alza hidráulica (absorción profunda con
liberación superficial).

Luis Rázuri Ramírez

CIDIAT-ULA 163
 MOVIMIENTO DE AGUA EN EL SUELO

Flujo saturado: Ley de Darcy

Experimento de Darcy

A mediados del siglo XIX Henri Darcy, un ingeniero francés,

fué contratado para limpiar las aguas públicas en su municipio
de Dijón (Francia), muy contaminadas por la pujante industria
de fabricación de mostaza y otros usos públicos.

Su idea fue utilizar filtros de arena para limpiar las aguas

antes de introducirlas en el sistema de distribución. Sin
embargo pronto encontró que con los filtros de por medio la
cantidad de agua disponible por unidad de tiempo (caudal) se
reducía enormemente.

Luis Rázuri Ramírez

CIDIAT-ULA 164
Tras una serie de experimentos con columnas de arena (Figura) estableció en
1856 que el flujo del agua a través de un medio poroso saturado, como en la
columna de arena, es análogo al flujo en una tubería.

La relación encontrada en sus experimentos fue que el flujo unitario (q = Q/A)

[LT-1] a través de la columna de arena saturada (todos sus poros llenos de agua)
de sección A, era directamente proporcional a la diferencia entre las alturas de
los depósitos de alimentación y recogida del agua (∆H) e inversamente
proporcional a la longitud, L, de la columna.
Luis Rázuri Ramírez
CIDIAT-ULA 165
La relación entre las variables observadas se conoce como la Ley de
Darcy y es la ecuación fundamental del flujo en medio saturado:

V H1  H 2 H H 2  H1 dH
q   Ks  Ks   Ks   Ks   K s H
tA L L x2  x1 ds

El flujo unitario es llamado también flujo de velocidad de Darcy

porque asume que la descarga ocurre a través de todo el material sin
atender al hecho de que los sólidos ocupan la mayor parte del volumen
(tiene unidades de velocidad, LT-1).
La ley es una relación para cuantificar la energía (diferencias de
potenciales entre depósitos) que se requiere para mover agua a
través del medio. La energía se disipa a causa de la fricción entre el
agua que se mueve y las paredes de los poros. La ecuación, establece
que esta pérdida de energía es proporcional al flujo.

Luis Rázuri Ramírez

CIDIAT-ULA 166
 Flujo no saturado: Ecuación de Darçy-Buckingham

Este es el de más interés en la mayoría de los suelos en que no hay

problemas de drenaje y normalmente se encuentran en estado
insaturado. Si se construyera un aparato experimental como el de la
Figura, veríamos que la columna de suelo está insaturada a ambos
lados y seguimos observando flujo.

Luis Rázuri Ramírez

CIDIAT-ULA 167
Para calcular el flujo insaturado en la
Figura, comenzaríamos con el suelo
saturado y los depósitos a la altura de la
columna. Dichos depósitos, unidos con
mangueras flexibles a la columna, se
bajarían por debajo de la columna a
ambos lados provocando una succión y
salida del líquido de la columna hacia
dichos depósitos.
El líquido quedaría desplazado por el aire
que entraría por los conductos en la
parte superior de la columna, con lo que
tendríamos condiciones de insaturación.
Si mantuviéramos los depósitos en una posición fija, con el depósito de la izquierda
más elevado que el de la derecha como en la Figura, observaríamos que se produce
flujo de izquierda a derecha que mediríamos a la salida del depósito derecho.
Igualmente observaríamos que al ser la succión distinta a ambos lados de la
columna, también lo será su contenido de humedad, que variará a lo largo de la
columna como se ve en el diagrama inc1uído en el centro de dicha Figura.
Luis Rázuri Ramírez
CIDIAT-ULA 168
Siguiendo una metodología como la propuesta por Darcy (variando
longitudes y sección de columna), manteniendo la posición de los
depósitos fija, veríamos que el flujo sigue una ecuación de la misma
forma que la ecuación de Darcy.
Al variar la posición relativa de los depósitos, en definitiva cambiando
las condiciones de succión de la columna o lo que es lo mismo su
contenido de humedad, obtendríamos otra relación distinta a la
anterior, caracterizada por un nuevo valor de la conductividad
hidráulica (K).
Si se compilara todos los resultados para distintos niveles de
insaturación en la columna (distintos gradientes hidráulicos y valores
de K resultantes) se tendría la ecuación de flujo in saturado o
ecuación de Darcy-Buckingham,
H H
q   K ( )   K ( h)
x x
Luis Rázuri Ramírez
CIDIAT-ULA 169
 H H
q   K ( )   K ( h)
x x
En este caso K es una función del contenido de humedad y no una
constante como en el caso de flujo saturado. La gran dependencia
que K tiene de θ (y por tanto de h) se manifiesta en que al pasar la
tensión de 0,1 a 1 bar, K puede disminuir hasta 100.000 veces, como
se observa en la Figura.

Luis Rázuri Ramírez

CIDIAT-ULA 170
Nótese que en el flujo saturado (ψ < 0) los suelos de mayor
conductividad (K=Ks) son los arenosos (picón en Figura b),
mientras que en el insaturado, excepto a valores bajos de
tensión, los suelos arcillosos suelen tener una K mayor que
los arenosos para un mismo potencial mátrico (ψ >10 en Fig
b.).

Algunas de las expresiones matemáticas utilizadas en la

descripción de la conductividad insaturada en los suelos se
recogen en la Tabla siguiente:

Luis Rázuri Ramírez

CIDIAT-ULA 171
 Tabla : Ecuaciones empleadas para describir la conductividad insaturada

Ecuación Conductividad hidráulica insaturadaa

αψ
Gardner Kse

Campbell Ks (θ/P)n
   r  2
Brooks y Corey K s  n ;  3 
  s  r  
2
1/ 2   1/ m m 

   r  
1  1     
van Genuchten K s   
r
  
  s  r      s  r   
   

a) donde α, λ y n son constantes empíricas, Ks es la conductividad hidráulica

saturada y P es la porosidad; ψ la succión (-h); θs y θr son los contenidos de
humedad a saturación y residual respectivamente.
Luis Rázuri Ramírez
CIDIAT-ULA 172
Es importante recordar que en el flujo insaturado no todos
los poros están llenos de agua.

Esto hace que la proporción de área transversal que conduce

el agua sea menor que en el flujo saturado y la tortuosidad
del flujo aumente.

Los conductos que más agua conducen en flujo saturado (los

de mayor diámetro) son los primeros en vaciarse a
insaturación y por lo tanto los que no conducen agua. Es
importante recordar que en el flujo insaturado el potencial
de presión (mátrico) del suelo es negativo.

La distribución de valores de conductividad hidráulica en un

campo suele seguir una función de distribución de
probabilidad log-normal, y su media se aproxima por la media
geométrica. Luis Rázuri Ramírez
CIDIAT-ULA 173
 Ecuación general del flujo no saturado.
Ecuación de Richards

La ecuación unidimensional del balance de masa o de continuidad

Entradas (l) - Salidas (O) = Cambio de almacenamiento (ΔS)
aplicada a un volumen elemental (∆x en 1 dimensión) de suelo insaturado
resulta:
I=q
 q/  x
q q+∆q
O = q+∆q = q + ∆x
∆x

x ∆S = ∆t (Volumen de agua)≈∆t(θ ∆x)=∆x ∆t(θ)=∆x  θ/  t

donde vol. de agua = θ V suelo ≈ θ∆x (en una dimensión)

Sustituyendo en la ecuación de continuidad se tiene:

 q

t x
donde t es tiempo, q es flujo y x es la dimensión considerada
Luis Rázuri Ramírez
CIDIAT-ULA 174
La ecuación de continuidad expresa que la variación en el contenido de
humedad en dicho volumen elemental en el tiempo ( / t ) es igual a la
variación de flujo (q / x) entre la entrada y salida del mismo.

Combinando la ecuación de Darcy-Buckingham para flujo no saturado y la

ecuación de continuidad anterior se obtiene para un sistema
unidimensional con flujo horizontal (H = h+z →dH = dh+dz = dh+0 = dh),

  H  h
 ( K ( h) )  ( K ( h) )
t x x x x
En el caso de flujo vertical (dirección de eje z), el gradiente es
(h+z)/z, y por tanto la ecuación anterior pasa a ser:
  H  h
 ( K ( h) )  ( K ( h) )
t z z z z
Lo que se conoce como ecuación de Richards unidimensional. Aquí
consideramos z creciente en el sentido hacia arriba.
Luis Rázuri Ramírez
CIDIAT-ULA 175
Para simplificar el tratamiento matemático de la ecuación de Darcy o
la última y convertirlas en ecuaciones similares a las de la difusión
de calor, para las cuales ya se han desarrollado soluciones analíticas,
se hace el siguiente cambio:
 h 

x  x
Ahora la ecuación de Darcy se convierte en:
h 
q   K ( )
 x

Si ahora se define, la Difusividad hidráulica como:

h
D ( )  K ( )

Luis Rázuri Ramírez
CIDIAT-ULA 176
Las ecuaciones de Darcy-Buckingham se convierten en:
  
 ( D( ) )
t x x
   K ( )
 ( D( ) )
t z z z
A D(θ) se le denomina difusividad hidráulica y a la ecuación anterior
forma difusiva o forma-θ de la ecuación de Richards.
Además de la simplificación matemática, el empleo de D(θ) es
conveniente porque su rango de variación es menor que el de K(θ) al
variar el contenido de humedad.
Conviene recordar que aunque se emplea el término difusividad
(tomado por analogía con el análisis de flujo de calor) en realidad el
movimiento de agua en el suelo no es por difusión sino por flujo de
masa o convectivo.
Luis Rázuri Ramírez
CIDIAT-ULA 177
Otra formulación posible de la ecuación de Richards, válida para todo
el rango de h (suelo saturado e insaturado), se obtiene haciendo el
siguiente cambio:
  h h
  C ( h)
t h t t
donde C(h) se denomina capacidad, y la ecuación de Darcy-
Buckingham-2, se puede expresar como:

h  h K (h)
C ( h)  ( K ( h)  )
t z z z

Esta ecuación se denomina forma capacitiva o forma-h de la

Ecuación de Richards. Esta ecuación sin embargo obliga a usar dos
funciones de h que varían muy rápidamente al variar ésta, lo que
resulta en problemas de cara a su solución numérica.

Luis Rázuri Ramírez

CIDIAT-ULA 178
Infiltración del agua
en el suelo

Luis Rázuri Ramírez

CIDIAT-ULA 179
La conductividad hidráulica está afectada por la
estructura del suelo; es mayor en un suelo poroso,
fracturado y agregado, que en un suelo compacto y
denso.

La conductividad es función no sólo de la porosidad total,

sino principalmente del tamaño de los poros conducentes.

En muchos suelos, la conductividad hidráulica no es

constante en el tiempo, debido a procesos químicos,
físicos y biológicos

Luis Rázuri Ramírez

CIDIAT-ULA 180
 INFILTRACIÓN
Descripción
Se denomina infiltración al proceso de entrada de agua en el suelo.
Durante una tormenta la relación que hay entre el ritmo de aporte
de agua a la superficie del suelo y el de infiltración determina que
proporción de agua entra en el suelo y cual queda en superficie
disponible para la escorrentía superficial.
La infiltración instantánea real (f) se define como el flujo de agua
que entra en el suelo por unidad de superficie. Es por lo tanto el
valor que q, flujo en medio poroso, alcanza en la superficie del suelo
(f = q) y sus dimensiones son de velocidad [LT-1].
La infiltración acumulada real (F) tiene unidades de longitud [L].
Se llama infiltración potencial, infiltrabilidad o capacidad de
infiltración o (fp) a la velocidad de infiltración que se produce
cuando el suelo está encharcado con una lámina de agua pequeña en
superficie. La capacidad de infiltración acumulada bajo
encharcamiento se denota porLuis
Fp.Rázuri Ramírez
CIDIAT-ULA 181
Si la superficie insaturada del
suelo es bruscamente encharcada
(por lluvia intensa o inundación),
entonces el gradiente hidráulico
(ecuación de Darcy) en la capa
superficial es muy grande (∆H>>L
para L pequeño).
Sin embargo, a medida que se va
infiltrando el agua el gradiente se
reduce y por lo tanto la velocidad
de infiltración disminuye en el
tiempo hasta aproximarse a un valor final que teóricamente es
igual al valor de la conductividad hidráulica del suelo saturado
(∆H≈L; ∆H/L≈1; f = Ks)
Luis Rázuri Ramírez
CIDIAT-ULA 182
En la práctica el valor final puede ser menor que Ks debido
al aire atrapado en los poros de suelo a medida que el
frente húmedo va descendiendo en el suelo desde la
superficie.

Luis Rázuri Ramírez

CIDIAT-ULA 183
 EN EL PRIMER MINUTO

ΔH/L es grande:
ΔH = Saturado – seco
L = pequeño

PROFUNDIDAD
DE LAS RAICES

Luis Rázuri Ramírez

CIDIAT-ULA 184
 EN LA MITAD DEL PROCESO

L
ΔH/L es intermedio:
ΔH = Saturado – seco
L = es mas grande
PROFUNDIDA
D DE LAS
RAICES

Luis Rázuri Ramírez

CIDIAT-ULA 185
 AL FINALIZAR EL PROCESO

L ΔH/L es 1:
ΔH = Saturado
L = es máximo

PROFUNDIDAD
DE LAS RAICES EL SUELO ESTA
SATURADO
VI =f= Ks

Luis Rázuri Ramírez

CIDIAT-ULA 186
Factores que afectan a la infiltración
Según Skaggs y Khaleel (1982) son los siguientes:
•Tipo de suelo (textura): A mayor conductividad hidráulica
saturada (Ks) mayor será la velocidad de infiltración f.
•Contenido inicial de humedad (θi): Cuanto más seco esté
inicialmente el perfil mayor será la velocidad de infiltración
inicial.

Luis Rázuri Ramírez

CIDIAT-ULA 187
• Intensidad de lluvia: Valores superiores a la capacidad de
infiltración potencial fp (definida por todos los demás
factores) provocan exceso de lluvia en la superficie
encharcamiento o escorrentía, pero no afectan a la
infiltración real, f .

•Sellado y costra superficial: La estructura de la superficie

del suelo influye mucho en la conductividad hidráulica y esta
determina en gran medida la velocidad de infiltración. Una
superficie compactada limita mucho esta velocidad.

•Aire atrapado: Reduce sensiblemente la capacidad de

infiltración, sobre todo en fases iniciales del proceso de
infiltración y resulta en valores finales de f inferiores al valor
de la conductividad hidráulica del suelo saturado.
Luis Rázuri Ramírez
CIDIAT-ULA 188
Presencia de horizontes en el suelo: Tanto estratos
arenosos como arcillosos pueden ser un obstáculo para la
infiltración, aunque por diferentes motivos: el estrato
arcilloso retarda el flujo debido a su baja conductividad en
saturación, mientras que el estrato arenoso dificulta el avance
del frente húmedo debido a que en él hay condiciones de no
saturación y en estas condiciones el suelo arenoso suele tener
menor conductividad hidráulica que el arcilloso.

Longitud del periodo de infiltración: La velocidad de

infiltración disminuye en el tiempo hasta un valor mínimo, con
los que el valor instantáneo de f depende del tiempo.

Luis Rázuri Ramírez

CIDIAT-ULA 189
 Consideraciones de importancia

 Infiltración del agua en el suelo y su aplicación en el

diseño y operación del riego superficial

 El movimiento cíclico del agua en el suelo comienza con

la entrada de agua al perfil a través del proceso de
infiltración, continua con su almacenamiento temporario
en la zona de raíces...

 Termina con su remoción del suelo por el drenaje, la

evaporación o la absorción por las plantas.

Luis Rázuri Ramírez

CIDIAT-ULA 190
 El agua puede entrar a través de toda la superficie
uniformemente (inundación o lluvia). también puede
penetrar a través del perímetro mojado de surcos o
acequias e incluso infiltrarse hacia arriba en el perfil
desde una napa freática.

 La velocidad de infiltración determina la cantidad de

agua de escurrimiento superficial y con ello, el peligro
de erosión hídrica.

Luis Rázuri Ramírez

CIDIAT-ULA 191
 La velocidad de infiltración también determina los
tiempos de riego (duración de cada evento) y el diseño
de los sistemas de riego (tamaño de la unidad de riego
y caudales óptimos).

 La velocidad de infiltracion del suelo disminuye con el

tiempo, durante el proceso de infiltración.

 La disminucion en la velocidad de infiltración es

consecuencia de una disminución inevitable del
gradiente de succión matrica, que ocurre a medida que
tiene lugar el proceso.
Luis Rázuri Ramírez
CIDIAT-ULA 192
Curva de Velocidad de
Infiltración

Luis Rázuri Ramírez

CIDIAT-ULA 193
Luis Rázuri Ramírez
CIDIAT-ULA 194
Luis Rázuri Ramírez
CIDIAT-ULA 195
LA VELOCIDAD DE INFILTRACION DEL SUELO ES ALTA EN LOS
PRIMEROS ESTADOS DEL PROCESO, TIENDE A DECRECER EN
FORMA MONÓTONA, ACERCANDOSE A UNA VELOCIDAD
CONSTANTE. (Capacidad de infiltración final), cuyo valor es = a la
conductividad hidráulica saturada del suelo.

LA CURVA DE VELOCIDAD DE INFILTRACION EN EL TIEMPO ES

FUNCIÓN DEL CONTENIDO DE AGUA INICIAL DEL SUELO y
DE LA CAPACIDAD DE ALMACENAMIENTO DE ÉSTE (TEXTURA,
ESTRUCTURA Y ESTRATIFICACIÓN DEL PERFIL.

LA INFILTRACION ACUMULADA ES LA INTEGRACION EN EL

TIEMPO DE LA VELOCIDAD DE INFILTRACIÓN. PRESENTA UNA
DEPENDENCIA EXPONENCIAL EN EL TIEMPO, CON UNA
PENDIENTE QUE DECRECE GRADUALMENTE.

Luis Rázuri Ramírez

CIDIAT-ULA 196
DISTRIBUCION DE LA HUMEDAD DURANTE LA INFILTRACION

LA SUPERFICIE DEL SUELO ESTA SATURADA SOLAMENTE

HASTA UNA PROFUNDIDAD DE ALGUNOS mm HASTA UN PAR
DE cm.

MAS ABAJO DE LA ZONA SATURADA EXISTE UNA ZONA DE

ELONGACION, APARENTEMENTE UNIFORME EN CONTENIDO
DE AGUA Y BASTANTE CERCANA A LA SATURACION; SE
CONOCE COMO ZONA DE TRANSICION.

Luis Rázuri Ramírez

CIDIAT-ULA 197
Luis Rázuri Ramírez
CIDIAT-ULA 198
 DESPUES DE LA ZONA DE TRANSICION, ESTÁ LA ZONA DE
HUMEDECIMIENTO, EN LA CUAL EL CONTENIDO DE AGUA EN EL
SUELO DISMINUYE CON LA PROFUNDIDAD, HASTA EL FRENTE DE
HUMEDECIMIENTO

Luis Rázuri Ramírez

CIDIAT-ULA 199
Luis Rázuri Ramírez
CIDIAT-ULA 200
 MEDICION DE LA INFILTRACION

TODOS LOS METODOS PARA MEDIR LA INFILTRACION EN

CONDICIONES DE CAMPO SE BASAN EN LA INUNDACION Y
ESTANCAMIENTO DEL AGUA EN LA SUPERFICIE y EN
MEDICIONES DE CAUDAL DE ENTRADA y SALIDA EN SURCOS
y MELGAS

METODOS PARA MEDIR LA VELOCIDAD DE INFILTRACION

* CILINDROS INFILTROMETROS
* SURCOS INFILTROMETROS
* METODO DE REPRESA

Luis Rázuri Ramírez

CIDIAT-ULA 201
Luis Rázuri Ramírez
CIDIAT-ULA 202
Luis Rázuri Ramírez
CIDIAT-ULA 203
Luis Rázuri Ramírez
CIDIAT-ULA 204
 Modelos aproximados
Ecuación de Philip: Infiltración horizontal.

Partiendo de la ecuación general del flujo horizontal:(1)

  
 ( D ( ) )
t x x
Esta ecuación se puede resolver bajo unas condiciones iniciales y de
contorno como las que existen en la infiltración horizontal: (2)
 i x  0 t  0
 o x  0 t  0
donde θi y θ0 son la humedad volumétrica inicial y final,
respectivamente.

Según Philip, con el cambio de variable (transformación de

Boltzmann): (3)
 ( )  x t
1/ 2

Luis Rázuri Ramírez

CIDIAT-ULA 205
La ecuación inicial (1) se convierte en: (4)
d d d
  ( D( ) )
2d d d
con las condiciones de contorno: (5)

 i   
 o   0

La ecuaciones (4 y 5) requiere técnicas especiales para su resolución. La

infiltración acumulada F es:
0
F   xd
i

Sustituyendo (3) en (6) se tiene:

Luis Rázuri Ramírez

CIDIAT-ULA 206
 0
F (t )  t   ( ) d  s t1/ 2
1/ 2
i

donde s es el valor de la integral y denominado sorptividad. Para el

cálculo de la velocidad de infiltración se obtiene:

dF 1 1/ 2
f (t )   st
dt 2

Luis Rázuri Ramírez

CIDIAT-ULA 207
Ecuación de Philip: Infiltración vertical.

Aplicando la ecuación general de flujo vertical (forma difusiva de la

ecuación de Richards),
   K ( )
 ( D ( ) 
t  z z z
y las condiciones iniciales y de contorno similares al caso horizontal
(ec. 5), Philip llegó a unas soluciones aproximadas por una serie
infinita en la que se despreciaban los términos de orden 2 en
adelante:
F (t )  st
1/ 2
 K st
f (t ) 1 / 2st 1/ 2  K s

Para los valores pequeños de t el primer término de f(t) predomina

sobre el segundo. Lo contrario ocurre para los valores grandes de t.
El valor de Ks es precisamente la conductividad hidráulica del suelo
a saturación (Tabla).
Luis Rázuri Ramírez
CIDIAT-ULA 208
 Tabla : Valores de Ks para la ecuación de Philip.

Tipo de suelo Ks(mm/h)

Arenoso >20

Franco 5-10

Arcilloso 1-5

Obsérvese que a medida que t →∞ entonces, f(t) →Ks.

Los dos términos de la ecuación de Philip representan los efectos
del potencial de succión y el gravitacional respectivamente.
Esta es la razón por la que en el caso anterior de infiltración
horizontal aparece sólo el primer término (ecuaciones 7 y 8).

Luis Rázuri Ramírez

CIDIAT-ULA 209
 Ecuación de Horton
Horton presentó uno de los primeros modelos empíricos (basado en
comportamiento experimental) de la infiltración de los suelos en una
ecuación dependiente de tres parámetros,
dF
fp   f c  ( f 0  f c ) e kt
dt
donde fc es el valor final de la infiltración instantánea, f0 es el valor
inicial de infiltración instantánea (t=0), y k es un parámetro del suelo
que controla el ritmo en el que la infiltración desciende con el tiempo
[T-1] (Figura).

Skaggs y Khaleel (1982) señalan que los parámetros f0 y k deben

depender del contenido de humedad inicial del suelo, así como de la
intensidad de lluvia, y para suelos homogéneos fc será ligeramente
menor que la conductividad hidráulica saturada.
Es importante señalar que la bondad de ésta o cualquier otra de las
ecuaciones presentadas dependerá de la calidad de los parámetros
usados en las mismas, que deben ser evaluados experimentalmente.
Luis Rázuri Ramírez
CIDIAT-ULA 210
Luis Rázuri Ramírez
CIDIAT-ULA 211
 Modelo de Green-Ampt

Para la deducción de la ecuación se

consideran las siguientes hipótesis
para describir la infiltración en un
suelo:
1. Suelo encharcado con una lámina
de agua ho desde el principio (t=0)

2. Frente de avance de humedad plano

(flujo en forma de pistón).

3. Suelo profundo y homogéneo, con

un contenido inicial de humedad
uniforme en profundidad.

Aplicando la ecuación de Darcy, en el perfil de suelo anterior (desde la

superficie al punto de profundidad L), haciendo q =fp, se tiene:

Luis Rázuri Ramírez

CIDIAT-ULA 212
 H1  H 2 (h  0)  (S av  L h  S av  L
f p  Ks ( )  Ks ( 0 )  Ks ( 0 ) (1)
L L L

y considerando que ho≈O y teniendo en cuenta que para el perfil

anterior la infiltración acumulada en el tiempo es Fp = (θs – θi)L, se
obtiene una expresión para la infiltración instantánea bajo
encharcamiento o capacidad de infiltración,

K s S av M B
f p  Ks   A (2)
Fp Fp
donde, Ks es la conductividad hidráulica saturada de Darçy,
M = θs - θi, es la porosidad llenable de agua, Sav es la succión
efectiva en el frente de avance, y θs, θi son los contenidos de
humedad a saturación e inicial, respectivamente.

Luis Rázuri Ramírez

CIDIAT-ULA 213
La ecuación anterior (2) no contiene la variable tiempo (t), necesaria
para obtener perfiles de infiltración en el tiempo

Se puede obtenerla ya que fp = dFp/dt, e integrando la última

ecuación con la condición Fp(t=0)=0 , tenemos:

 Fp 
K s t  Fp  S avm ln 1   (3)
 S a vM 

que es una ecuación implícita de la infiltración acumulada F, es

decir no se puede despejar F=F(t). Esto obliga a resolver la ecuación
iterativamente mediante aproximaciones sucesivas o a tantear
valores de F para obtener la t deseada. Una vez obtenida la F para
la t deseada, la infiltración instantánea se obtiene sustituyendo en
la ecuación (2)
Luis Rázuri Ramírez
CIDIAT-ULA 214
Los valores de los parámetros de la ecuación de Green-Ampt deben
medirse experimentalmente para obtener los mejores resultados
(Skaggs, 1969). En ausencia de datos experimentales pueden
estimarse partiendo de la textura del suelo según proponen diversos
autores (Rawls and Brakensiek, 1982; McCuen et al, 1981). En la
Tabla se presentan algunos de los valores estimados.

Valores de los parámetros de Green-Ampt para distintas texturas (adpt. Rawls and
Brakensiek, 1982)

Textura(USDA) Ks(m/s)x10-6 Sav(m) θe?θs(m3/m3)

Arcillosa 0.167 0.0639-1.565 0.427-0.523
(0.3163) (0.475)
ArciUo-Arenosa 0.333 0.0408-1.402 0.370-0.490
(0.2390) (0.430)
Franco-Arcillosa 0.556 0.0479-0.9110 0.409-0.519
(0.2088) (0.464)

Luis Rázuri Ramírez

CIDIAT-ULA 215
Textura(USDA) Ks(m/s)x10-6 Sav(m) θe? θs(m3/m3)
Arcillo-Limosa 0.278 0.0613-1.394 0.425-0.533
(0.2922) (0.479)
Franco-Arcillo- 0.556 0.0567-1.315 0.418-0.524
Limosa (0.2730) (0.471)
Franco-Arcillo- 0.0442-1.080 0.332-0.464
0.833
Arenosa (0.2185) (0.398)
Franca 3.67 0.0133-0.5938 0.375-0.551
(0.0889) (0.463)
Franco-Limosa 1.89 0.0292-0.9539 0.420-0.582
(0.1668) (0.501)
Franco-Arenoso 6.06 0.0267-0.4547 0.351-0.555
(0.1101) (0.453)
Arenosa-Franca 16.6 0.0135-0.2794 0.363-0.506
(0.0613) (0.437)
Arenosa 65.4 0.0097-0.2536 0.374-0.500
(0.0495) (0.437)

Luis Rázuri Ramírez

CIDIAT-ULA 216
Ejemplo 1

Calcular la infiltración instantánea (fp) y la acumulada (Fp) durante

una hora en un suelo franco limoso encharcado con θs = 0,501,
θi = 0,161, Sav = 16,7 cm, Ks = 0,65 cm/h

Como sabemos, M = θs - θi,= 0,340; SavM= 5,68 cm.

Para calcular la F al final de la hora (t=1) se reordena la ecuación 3,

 Fp 
Fp  K s t  S av M ln 1   (4)
 S av M 

Luis Rázuri Ramírez

CIDIAT-ULA 217
y se prueba un valor de Fp en la parte derecha de Ec. 4 y se obtiene
un nuevo valor de Fp que posteriormente se vuelve a introducir en la
misma ecuación, y así, sucesivamente hasta que el valor que sale de Fp
es parecido al que entra hasta un nivel de precisión elegido (2
decimales suficiente).

Tabla. Procedimiento de iteración para el método de Green-Ampt (Ejemplo 1).

iteración F iteración F iteración F

O 0.65 9 3.089 18 3.165
1 1.265 10 3.116 19 3.166
2 1.792 11 3.134 20 3.1
3 2.208 12 3.146 21 3.167
4 2.515 13 3.153 22 3.167
5 2.732 14 3.158 23 3.167
6 2.881 15 3.161 24 3.167
7 2.980 16 3.163 25 3.167
8 3.046 17 3.165 26 3.167

Luis Rázuri Ramírez

CIDIAT-ULA 218
Con el valor obtenido de Fp(t=1)=3,17 cm se entra en la Ec. 2 y se
obtiene fp= 1,81 cm/h. (hoja de cálculo)

Siguiendo el procedimiento anterior para valores t=0; 0,1; 0,2; 0,3;

0;4; así hasta 1,0 horas se representan los valores en la Figura.
Obsérvese cómo el frente húmedo (L) va avanzando en profundidad
(eje Y de la derecha).

Luis Rázuri Ramírez

CIDIAT-ULA 219
 Infiltración de lluvia constante sin
encharcamiento inicial

Uno de los factores que afectan al comportamiento de la infiltración

en el tiempo es la intensidad de lluvia, i.

Sin embargo la relación entre infiltración instantánea y acumulada no

depende de la i. La dificultad con que se encuentra es que en
condiciones naturales el suelo no se encharca al principio de un suceso
de lluvia, en contra de la hipótesis (1) de la ecuación de G-A.

Mein y Larson (1973) aplicaron la ecuación de G-A a éste caso

calculando primero la infiltración acumulada (F = Fp) que hace falta
para llegar al tiempo de encharcamiento (t = tp)

El subíndice "p" en las infiltraciones (Fp,fp) denota las

correspondientes al suelo encharcado según la teoría de G-A, su
ausencia (F, f) indica las infiltraciones reales para un suelo esté o no
encharcado.

Luis Rázuri Ramírez

CIDIAT-ULA 220
Existen tres principios invo1ucrados en el proceso: i) antes de que
llegue el tiempo de encharcamiento toda la lluvia se infiltra; ii) la tasa
de infiltración f es función de la acumulada F; iii) el encharcamiento
ocurre cuando la tasa de infiltración, f, es menor o igual que la
intensidad de lluvia, i

Luis Rázuri Ramírez

CIDIAT-ULA 221
Así, para una precipitación constante tendremos:
f  i.............................0  t  t p
f  f p  i.....................t  t p
f  f p ( Ec.1..........t  t p (5)
en el momento de encharcamiento (t=tp), al sustituir f = fp = i,
Fp = [Link] en la Ec. 2, se tendrá un tiempo de encharcamiento,
K s .S av .M
tp  ; Fp  i.t p (6)
i(i  K s )

Para obtener la tasa de infiltración f después del encharcamiento,

se construye una curva de infiltración potencial empezando en un
tiempo to (tiempo equivalente) de manera que f y F, en tp , sean
iguales a las que se observan bajo una lluvia que empezara en el
tiempo 0 (líneas gruesas punteadas en la Figura).

Luis Rázuri Ramírez

CIDIAT-ULA 222
Al introducir este concepto, Mein y Larson (1971) propusieron una
ecuación análoga a la Ec. 4,
 Fp 
K s (ti  t p  to )  Fp  MS av ln 1   (7)
 MSav 

de donde se puede calcular el t0 tras conocer tp y Fp con la Ec. 6 al

hacer t=tp,
Fp M .S av  Fp 
to   
ln 1   (8)
Ks Ks  M .S av 

Nótese que si el tiempo al encharcamiento calculado al principio (tp)

es mayor que el tiempo total de la tormenta, entonces no hay
encharcamiento y la ecuación de Green-Ampt no es aplicable, en su
lugar F=i.∆t y f=i.

Luis Rázuri Ramírez

CIDIAT-ULA 223
 Ejemplo 2

Para los datos del ejemplo anterior (pero sin encharcamiento

superficial), calcular la infiltración para una intensidad de lluvia
constante de 2,5 cm/h

tp= 0,65x 5,68/[2,5(2,5-0,65)]=0,798 horas

Fp= tp.i = 0,80 x 2,5 = 1,996 cm

to= 1,996/0,65 – 5,68/0,65 ln(1+1,996/5,86) = 0,439 horas

De la Ec. 7 se puede calcular F=F(t) a partir de t>tp siguiendo el

procedimiento del apartado anterior (nótese que cada valor de F para
cada t se obtuvo de un proceso de iteración como el de la Tabla
anterior y la Figura)

Luis Rázuri Ramírez

CIDIAT-ULA 224
Tabla. Cálculo de la infiltración para suelo sin encharcamiento con lluvia constante.

t F f t F f
0.000 0.000 1.100 2.670 2.033
0.100 0.250 2.500 1.200 2.869 1.937
0.200 0.500 2.500 1.300 3.058 1.857
0.300 0.750 2.500 1.400 3.240 1.789
0.400 1.000 2.500 1.500 3.416 1.731
0.500 1.250 2.500 1.600 3.587 1.679
0.600 1.500 2.500 1.700 3.752 1.634
0.700 1.750 2.500 1.800 3.914 1.593
0.798 1.996 2.500 1.900 4.071 1.557
0.900 2.239 2.299 2.000 4.225 1.524
1.000 2.461 2.150

Luis Rázuri Ramírez

CIDIAT-ULA 225
 Redistribución del agua después de la
infiltración
Al acabar la infiltración, el movimiento del agua en el suelo no cesa
sino que sigue en función de los gradientes de potencial hidráulico.

Las capas superficiales, inicialmente más húmedas, pierden agua

mientras que otras más profundas aumentan su θ. La redistribución
es más o menos rápida dependiendo del tipo de suelo.

La redistribución determina cuanta agua pasará la zona radicular en

un periodo determinado después del riego.

1. Redistribución en perfiles húmedos

En este caso se puede suponer que el gradiente del potencial de

presión con la profundidades es muy pequeño al estar el perfil
húmedo en profundidad comparado con el gradiente de potencial de
gravedad, entonces resulta que:

Luis Rázuri Ramírez

CIDIAT-ULA 226
 H  z
q  ( K ( ) )   K ( ) (h  z )   K ( )   K ( )
z z z
(El signo negativo indica flujo hacia abajo ya que se definió z
creciente hacia arriba)

El caso más general (unidimensional) sólo se puede resolver con la

ecuación de flujo unidimensional.
  H
 ( K ( ) )
t z z
Si suponemos que inicialmente el perfil está saturado y que
durante la redistribución no hay evaporación ni infiltración,
entonces las condiciones iniciales y de contorno serían:
t  0, z  0,   s ( perfil inicialmen te saturado)
t  0, z  0, q  0 (evap. ó inf iltración  0)
t  0, z  zb , q  K ( b )
Luis Rázuri Ramírez
CIDIAT-ULA 227
donde la última condición de contorno supone que en el fondo del
perfil considerado H=z (se mantiene saturado en el punto z=zb, h=0
en el nivel freático).

Los resultados de este modelo (Hillel, 1994) son:

• El suelo arenoso contiene menos agua en saturación pero la pierde

más rápidamente al principio; al cabo de dos días ha perdido casi el
doble que el franco y unas cinco veces lo del arcilloso.

• Con el tiempo la velocidad de drenaje es: arcilloso>franco>arenoso.

• La redistribución es afectada por la presencia de estratos menos

permeables

Luis Rázuri Ramírez

CIDIAT-ULA 228
2. Redistribución en perfiles no húmedos.

En este caso la rapidez de redistribución depende además de las

propiedades hidráulicas del suelo, de la profundidad inicial del
suelo mojado y de la humedad de las capas más profundas.

Si el mojado inicial es poco profundo y el suelo profundo está

relativamente seco tenemos:

• La redistribución es rápida.

• El frente húmedo avanza rápido al principio y está bien

delimitado; poco a poco se hace más difuso y lento.

• En la redistribución se dan fenómenos de histéresis que

complican el análisis

Luis Rázuri Ramírez

CIDIAT-ULA 229
 Ecuación de Kostiakov

Se ha propuesto varias ecuaciones empíricas para expresar la

velocidad de infiltración como una función del tiempo, que pueden ser
representadas por una curva de forma hiperbólica. La ecuación de
Kostiakov (1932) expresa la velocidad de infiltración en un punto:

I at b

donde:
I es la velocidad de infiltración (L T-1) en mm min -1 o mm h-1

t es el tiempo de infiltración (1) en min o en h

a es un coeficiente, el cual representa la velocidad de

infiltración a t = 1,0; expresado en mm min-(1+b)
ó mm h-(1+b)
b es un exponente sin dimensión, siempre es negativo
con valores que van de 0 a - 1,0.
Luis Rázuri Ramírez
CIDIAT-ULA 230
Al integrar la ecuación entre límites t = 0 y t, se obtiene la
infiltración acumulada:
a b1
I cum  t
b 1
La ecuación también puede ser expresada así

I cum  A t B

Generalmente, las ecuaciones anteriores se ajustan muy bien a

la mayoría de las condiciones de la práctica del riego por
superficie. Sin embargo, hay algunos casos en los cuales la
velocidad de entrada alcanza un valor constante dentro del
período de infiltración; si es así, la ecuación primera se
convierte en:
I at c
b

Luis Rázuri Ramírez

CIDIAT-ULA 231
donde: c es la velocidad de infiltración, constante para t = ∞
La infiltración acumulada, lcum será por lo tanto:
a b1
I cum  t ct
b 1
con b = - 0,50 esta ecuación es equivalente a la ecuación de Philip

La ecuación de Kostiakov ha sido usada en suelos y riego,

principalmente por resultar práctica.

En la actualidad, distintas ecuaciones usadas para proyectar los

métodos de riego por superficie involucran los parámetros de la
ecuación de Kostiakov, especialmente el exponente b ó B.

Luis Rázuri Ramírez

CIDIAT-ULA 232
 Integrando

Infiltración Básica ∆ ≤ 10%

Derivando

Despejando t
Tiempo básico

Luis Rázuri Ramírez Infiltración básica en 233

Infiltración básica en mm/hr CIDIAT-ULA mm/min
 Ecuación del Servicio de Conservación de Suelos

El Servicio de Conservación de Suelos de los Estados Unidos ha

obtenido en pruebas de campo una gran cantidad de curvas
representativas de muy diversas condiciones edáficas, que han
sido agrupadas en "familias de curvas", las que responden a la
siguiente ecuación:
I cum  A´t  C´
B´

La Tabla incluye los parámetros de la ecuación para diferentes

familias de infiltración y la Figura las curvas respectivas. Los
datos de una prueba de campo se representan en la Figura, lo
que permite identificar la curva que más se aproxime a los
puntos respectivos y, por lo tanto, la ecuación para Icum = f(t)
en la Tabla.
Luis Rázuri Ramírez
CIDIAT-ULA 234
 Ecuación del Servicio de Conservación de Suelos

Figura. Familias de curvas de infiltración acumulada, según el procedimiento del

Servicio de Conservación de Suelos (1974).
Luis Rázuri Ramírez
CIDIAT-ULA 235
Luis Rázuri Ramírez
CIDIAT-ULA 236
Infiltración promedio
La velocidad de infiltración promedio, I p, es la relación entre
la infiltración acumulada y el tiempo de entrada.
I cum
Ip 
t
Sustituyendo la ecuación de Icum en la ecuación anterior se obtiene
la ecuación de infiltración promedio en un punto determinado.

a b
Ip  t
b 1
Entonces, Ip, es la velocidad promedio de entrada de agua en el
suelo en un período t.

Luis Rázuri Ramírez

CIDIAT-ULA 237
Infiltración básica

La infiltración básica, Ib es otro parámetro que merece

consideración debido a su importancia al proyectar el riego a nivel de
predio.

De acuerdo con el Departamento de Agricultura de los Estados

Unidos, Servicio de Conservación de Suelos, la velocidad de entrada
básica es el valor instantáneo cuando la velocidad de cambio de
entrada, para un período estándar, es 10% o menos de su valor.

El tiempo en el cual I = Ib se encuentra igualando la primera

derivada de la ecuación de velocidad de infiltración con la ecuación
multiplicada por 0,1:

dI
  0,1 I
dt
Luis Rázuri Ramírez
CIDIAT-ULA 238
 Eficiencia de Riego
 Definición: Es la relación entre el agua
consumida por los cultivos y el agua derivada
desde la fuente
Vc
Ef 
Vd
 Componentes de la eficiencia
1. Eficiencia de conducción
2. Eficiencia de distribución
3. Eficiencia de aplicación
Ef  EfcxEfdxEf a
Luis Rázuri Ramírez
CIDIAT-ULA
Luis Rázuri Ramírez
CIDIAT-ULA
Luis Rázuri Ramírez
CIDIAT-ULA
Luis Rázuri Ramírez
CIDIAT-ULA
Luis Rázuri Ramírez
CIDIAT-ULA
Luis Rázuri Ramírez
CIDIAT-ULA
 Eficiencia de aplicación, lámina bruta de riego
En general, cuando se aplica un riego, no toda el agua
queda almacenada en la zona radicular, sino parte se
pierde por evaporación, escorrentía y percolación
profunda, siendo muy diferente la cuantía de cada tipo
de pérdida según el sistema de riego.

Conceptualmente la idoneidad del riego depende de:

 Cantidad de agua almacenada en la zona radicular
producido por el riego.
 Pérdidas por percolación profunda y escorrentía
superficial.
 Uniformidad de la lámina infiltrada.
 Déficit de humedad delLuissuelo después del riego.
Rázuri Ramírez
CIDIAT-ULA
 ÁREA QUE RECIBE MÁS QUE LA LAMINA NECESARIA (%)

Pev
0 50 100

Hb = db

ÁREA ADECUADAMENTE REGADA

LAMINA REQUERIDA
Hba = cantidad de agua infiltrada

(LÁMINA NETA)
Hr = dn
Hr + R

i
ÁREA
DISTRIBUCIÓN DEL AGUA
INFRARREGADA
INDICADA CON UN Cu

LAMINA DE AGUA

Diagrama típico de la distribución del agua en un riego a presión

(aspersión y localizado)
Luis Rázuri Ramírez
CIDIAT-ULA
Esta distribución de agua se ajusta, en la generalidad
de los casos del riego por aspersión y del riego
localizado, a una distribución normal (Gausiana), no
siendo así en el caso del riego por superficie.

Es importante mencionar que la falta de uniformidad

en los sistemas de riego a presión es consecuencia de :

 Variación de fabricación de los emisores.

Diferencia de presiones en las unidades de riego.
Envejecimiento de las partes.
Obstrucciones que puedan presentarse.
Luis Rázuri Ramírez
CIDIAT-ULA
Si Hn es la lámina media de agua que ha quedado
almacenada en la zona radicular (que será siempre
menor que Hr, lámina requerida, cuando exista una
zona de déficit)
Hp, la lámina media percolada por debajo de la zona
radicular y Hd, la lámina media de agua que representa
el déficit en la zona infra regada.

La calidad del riego, para el deseado porcentaje de

área bien regada o sobre regada puede definirse,
basándose en una serie de parámetros como:

Luis Rázuri Ramírez

CIDIAT-ULA
 Eficiencia de aplicación: Efa= Hn/Hba
Eficiencia de distribución: ED = Hr/Hba
Coeficiente de déficit: Cd = Hd/Hr
Factor de disponibilidad: Fd = Hn/Hr
Porcentaje de percolación: Cp = Hp/Hb
Porcentaje de escorrentía: Ce = He/Hb

Como: Hn + Hd = Hr ; se tendrá que Hn/Hr + Hd/Hr = 1

Por lo tanto: Cd + Fd = 1
También se cumplirá que: Cd = 1 - Fd

La producción máxima se logrará sólo si Cd = 0

Luis Rázuri Ramírez
CIDIAT-ULA
250
251
252
253