Infiltration (Hydrogeologie)

Infiltration bezeichnet den Prozess des Eindringens von Niederschlägen in den Erdboden; er ist ein wichtiger Teilprozess des Wasserkreislaufes. Anschlussprozesse können Grundwasserneubildung und Abflussbildung sein. Die Infiltrationsmenge wird gemessen mit den Einheiten Millimeter pro Sekunde oder praktischer pro Minute.

Einleitung

Die Infiltration wird maßgeblich von zwei Kräften beeinflusst, der Gravitation und der Kapillarität bzw. dem Matrixpotential. Weitere entscheidende Einflussfaktoren sind die Temperatur im Boden, die Wassersättigung des Bodens (Vorfeuchte), der Bedeckungsgrad (Vegetation) sowie die Niederschlagsintensität. So haben sandige Böden (63µm-2mm Korngröße), mit relativ großen, gut verbundenen Poren und einem großen Porenvolumen (30-45%) die höchsten Infiltrationsraten. Da eine Vegetationsdecke das Auftreffen von Regentropfen abdämpft, verhindert sie effektiv das oberflächliche Verschlämmen von Poren und schafft durch Ihr Wurzelnetzwerk zusätzliche weitere, gut wasserleitende Grobporen. Damit hat sie einen deutlich positiven Einfluss auf die Infiltrationseigenschaften eines Bodens (höchste Infiltrationsraten lassen sich in Waldböden verzeichnen).

Messung

Die Messung der potentiellen Infiltration geschieht mittels sog. Infiltrometer. Dabei handelt es sich meist um einen Ring mit einem

Ringinfiltrometer

definierten Querschnitt, welcher in den Boden eingestochen wird und dann bis zu einer ebenfalls definierten höhe mit Wasser gefüllt wird. Da über die Einfüllhöhe und den Ringdurchmesser das eingefüllte Wasservolumen bekannt ist, kann die potentielle Infiltration (Infiltration bei maximalem Wasserdargebot) über die Zeit die vergeht bis alles Wasser aus dem Ring versickert ist bestimmt werden.

Doppelringinfiltrometer

Durch die Verwendung von Doppelringinfiltrometern wird versucht einen seitlichen Verlust von Sickerwasser zu minimieren, dieser würde bei der Messung zu höheren Infiltrationsraten führen. Beim Doppelringinfiltrometer, wird zur Messung nur das innere Volumen betrachtet.

Aus

$f(t)=- \frac{dh(t) }{dt}= \frac{Q(t) }{A}$

folgt nach dem Gesetz von Darcy:

$k_f= \frac{l }{t}= \frac{h(0) }{h(t)}$

mit:

k_f: Durchlässigkeitsbeiwert [cms ^{− 1}];

l: Länge der Fließstrecke [cm];

t: Zeit bis der Anfangswasserstand h(0) auf den Wasserstand zur Zeit t (h(t)) abgesunken ist;

h(0): Anfangswasserstand zur Zeit t=0;

h(t): Endwasserstand zur zeit t

Abschätzen der tatsächlichen Infiltration

Abschätzung über die Wasserbilanz

Sind alle übrigen Parameter bekannt, kann die Infiltration F als Restglied der Wasserhaushaltsgleichung berechnet werden.

$F=B_I+P-T-ET-S-R-I_A-B_O \!\,$

mit:

F :Infiltration

B_I: Randzufluss (unterirdischer oder oberirdischer Zufluss in das betrachtete System)

B_O: Randabfluss (unterirdischer oder oberirdischer Abfluss aus dem System)

P :Niederschlag

ET: Evapotranspiration

S: Speicher

I_A: kurzzeitiger oberflächlicher Rückhalt (z.B. Muldenrückhalt)

R: Oberflächenabfluss

Je nach Bewertung der örtlichen Verhältnisse kann diese Gleichung auch um einzelne Glieder vereinfacht werden.

Green-Ampt Infiltrationsmodell

Das Green-Ampt Infiltrationsmodell bietet eine Abschätzung der tatsächlichen Infiltration unter Berücksichtigung verschiedener Bodenparameter, als da wären: Saugspannung, Porosität, hydraulische Leitfähigkeit und dem bodenunabhängigen Parameter Zeit. Die halb-physikalische Formel approximiert den Infiltrationsprozess durch ein Stufenprofil mit vollständiger Wassersättigung und der sog. Transportzone. Nur der Wassergesättigte Teil des Bodens wird betrachtet (k_f konstant). Die Formel leitet sich aus der Darcy-Weisbach-Gleichung ab, wobei das Prinzip auf der sog. Gradientenmethode basiert. Zu Beginn wirkt eine hohe Saugspannung (Matrixpotential), die mit der Zeit schwächer wird.

$\int_0^{F(t)} {1-\psi\,\Delta\theta\over F+\psi\,\Delta\theta}\, dF = \int_0^t K\,dt$

mit:

ψ: Saugspannung an der Feuchtefront;

θ: Wassergehalt (Vorfeuchte);

K: hydraulische Leitfähigkeit;

F: Akkumuliertes Volumen des infiltrierten Wassers.

Durch Integration kann die Gleichung entweder nach dem Infiltrationsvolumen oder nach der Anfangsinfiltrationsrate aufgelöst werden.

$F(t)=Kt+\psi \, \Delta\theta \ln \left[1+{F(t)\over \psi \, \Delta\theta}\right]$

$f(t)=K\left[{\psi \, \Delta\theta\over F(t)}+1\right]$

Modell nach Horton

Der Anfangsinfiltrationswert f₀ gilt zu Beginn des Infiltrationsprozesses. Mit zunehmender Dauer verringert sich dessen Auswirkung und es stellt sich der End-Infiltrationswert $f_\infty$ ein. Das Model beschreibt eine exponentielle Abnahme der Infiltrationsrate bis zur End-Infiltrationsrate bei Sättigung des Bodens.

$f(t)=f_{\infty} + (f_0 - f_{\infty}) * e^{-k*t}\!\,$.

mit:

f₀: Anfangsinfiltrationsrate (t = 0);

f(t): max. Infiltrationsrate zum Zeitpunkt t [mm * h ^{− 1}]  ;

$f_\infty$: max Infiltrationsrate bei Sättigung (t → $\infty$);

k: Rückgangskonstante;

Ähnlich wie bei dem Modell nach Green und Ampt, besitzt auch das Horton-Modell eine Komponente die den Rückgang der Infiltrationsrate beschreibt und eine die den konstanten Teil berücksichtigt. Die Infiltrationsrate wird nie 0 sondern nähert sich einem vom Gravitationspotential abhängenden Endwert an. Dieser Endwert müsste theoretisch gleich der gesättigten hydraulischen Leitfähigkeit k_f sein welche bei Green-Ampt explizit dafür eingesetzt wird (für $f_\infty=k_f$). Die Modellparameter für Horton werden im allg. abgeschätzt, für unbewachsenen feinsandigen Ton bzw. für Grasboden gilt:

Parameter	Ton	Grasboden
f0	210 mm/h	900 mm/h
$f_\infty$	2 mm/h	290 mm/h
k	0.8/min	2/min

Weitere Modelle

Weitere gängige Modelle sind:

Modell nach Kostiakov

Empirisches Modell

$f(t) = akt^{a-1}\!$

mit:

a

k als empirische Parameter.

Literatur

• Siegfried Dyck, Gerd Peschke: Grundlagen der Hydrologie. 3. stark bearbeitete Auflage. Verlag für Bauwesen, Berlin 1995, ISBN 3-345-00586-7.