5.1 Überblick
5.2 Makroporenfluss
5.3 Matrixfluss
5.4 Bodenwasserentzug durch Verdunstung
5.5 Hypodermischer Abfluss
5.6 Drainabfluss
5.7 Kapillarer Aufstieg

Die Modellierung der Bodenwasserdynamik erfolgt mit einem Mehrschicht-Kapazitätsmodell nach Koitzsch (1977) und Glugla (1969), welches die Bodenwasserdynamik mittels abgeleiteter bodenspezifischer Kennwerte wie Feldkapazität und Permanenter Welkepunkt zur Charakterisierung der Wasserspeicherung unter bestimmten Spannungsverhältnissen beschreibt. Das Originalmodell wurde jedoch insbesondere hinsichtlich der Terme zur Abbildung der Verdunstungsintensität überarbeitet, da im Rahmen des PSCN-Moduls der aktuelle Transpirationsbedarf durch vegetationsspezifische Pflanzenmodelle berechnet wird. Außerdem wurde das Modell um einen Ansatz zur Beschreibung der hypodermischen Abflussbildung, um ein Makroporenflussmodell, einen Ansatz zur Abbildung des kapillaren Aufstiegs und um ein Drainagemodell erweitert.

Grundannahme des empirischen Makroporenflussmodells ist, dass die Bildung von Makroporen von den Bodeneigenschaften, der Landnutzung und der aktuellen Bodenfeuchte beeinflusst wird. Liegen keine Eingangsinformation zu den Makroporen der einzelnen Bodenhorizonte vor, so werden diese modellintern anhand der Trockenrohdichte, dem Skelett- und dem Tonanteil abgeschätzt. Dabei wird angenommen, dass Makroporen nur bei einer Lagerungsdichte von kleiner als 2 g/cm³ und erst ab einen Tonanteil von 20 % auftreten (k_t=0 für Ton(z) < 20 %).

Bezüglich der Landnutzung wird vorausgesetzt, dass Makroporen vor allem unter Wald, unter Wiesen und bei konservierender Bodenbearbeitung auftreten.

Für jede Bodenschicht z_j wird eine Makroporosität MP (0 ≤ MP ≤ 0.9) wie folgt berechnet:

mit

r_q Schrumpfterm bei Dürre
z_max maximale Tiefe der landnutzungsgebundenen Makroporen [mm]
z aktuelle Tiefe [mm]
t Zeit
q Wassergehalt [mm]
q_WP Wassergehalt am Permanentwelkepunkt [mm]
q_FK Wassergehalt bei Feldkapazität [mm]
k_s,k_t,k_v Parameter (k_v=0, wenn z>z_max)

In der obersten Bodenschicht ist MP der Anteil des infiltrierenden Niederschlags, der sofort in die darunter liegende Schicht weitergeleitet wird. Je nach Makroporosität der weiteren Schichten wird dieser weiter geleitet, bzw. kann anteilig in die Bodenmatrix infiltrieren.

Erreicht in den unteren Bodenhorizonten der Schichtwassergehalt die Sättigungsgrenze, so kann zusätzlich zur Perkolation durch die Bodenmatrix der Anteil MP(z_j,t) des mobilen Wassers in die darunter liegende Schicht transportiert werden. Dieser Makroporenfluss spielt insbesondere in Böden mit hohem Skelettgehalt im Unterboden eine Rolle.

Aufgrund der hohen Sensitivität der Bodenwasserhaushaltsmodellierung gegenüber der Makroporosität wird anstelle dieser internen Abschätzung empfohlen, die Makroporosität bei den Bodeneingangsdaten mit anzugeben.

Dem Mehrschicht-Kapazitätsmodell zur Beschreibung der Feuchtedynamik in der Bodenmatrix liegt die Annahme zugrunde, dass es erst nach Überschreiten der Feldkapazität zu einem Abfluss (vertikal und horizontal) aus der betrachteten Bodenschicht kommt. Anstelle der Kenntnis der hydraulischen Eigenschaften des Bodens wird nur noch der von der Bodenart abhängige Leitfähigkeitsparameter λ benötigt.

Die Änderung des Wassergehaltes q(z_j,t) der Schicht z_j in der Zeit t ergibt sich zu:

Nebenbedingungen:

 für 

Anfangswerte:

I Infiltrationsrate [mm/d]
P^d Abfluss über das Drainagesystem [mm/d]
P^g Perkolationsrate (Gravitationsabfluss) [mm/d]
P^h horizontaler (hypodermischer) Abfluss [mm/d]
P^k Kapillaraufstieg [mm/d]
P^m Makroporenfluss [mm/d]
E Wasserentzug durch Bodenevaporation und Transpiration [mm/d]  Θ Wassergehalt [mm]
Θ_WP Wassergehalt am Permanentwelkepunkt [mm]
Θ_PV Porenvolumen [mm]
Θ_FK Wassergehalt bei Feldkapazität [mm]

Die aktuelle Perkolationsrate P^g(z_j,t) aus der Schicht z_j ergibt sich zu:

BT(z_j,t) – Bodentemperatur der Schicht z_j [mm]

Die Wahl der Feldkapazität q_FK hat einen großen Einfluss auf das Simulationsergebnis. Der Leitfähigkeitsparameter λ lässt sich entweder aus dem Gehalt an abschlämmbaren Teilchen oder aus der gesättigten hydraulischen Leitfähigkeit ableiten (Glugla, 1969). Vereinfachend kann nach Koitzsch (1977) bodenartspezifisch ein Wert zwischen 0 und 1,3 mm^-1d^-1 bei Schichtmächtigkeiten von 10 cm angesetzt werden. Erfolgt keine Angabe zu diesem Parameter, so wird er modellintern aus der gesättigten hydraulischen Leitfähigkeit k [mm/h] berechnet.

λ – Leitfähigkeitsparameter nach Koitzsch
k – gesättigte hydraulische Leitfähigkeit [mm/h]

Der aktuelle Entzug E(t) setzt sich aus dem Bodenevaporationsanteil E_s und dem Transpirationsanteil E_TR zusammen.

Der Evaporationsentzug wird durch die potenzielle Bodenevaporation (Kapitel 2, Gl. 5) gesteuert und erfolgt bis zu einer bodenspezifischen Tiefe z_E. Standard für z_E in Nichtwaldböden ist 30 cm. Dieser Standardwert gilt nur für Böden ohne organische Auflage mit einer gesättigten Leitfähigkeit von weniger als 360 cm/d im Oberboden. Bei Auftreten von Schichten höherer Leitfähigkeit wird davon ausgegangen, dass diese die kapillare Nachlieferung aus den darunter liegenden Schichten unterbinden und somit als Sperrschicht für die Bodenevaporation wirken.

Es wird angenommen, dass aus der obersten Rechenschicht ein nur vom aktuellen Bodenwassergehalt limitierter Verdunstungsentzug stattfinden kann. Reicht dieser nicht zur Deckung des Bedarfs aus, so erfolgt ein Entzug aus tieferen Schichten. Bei Vorhandensein einer organischen Auflage (Waldböden, konservierende Bodenbearbeitung) kann nur aus dieser evaporiert werden.

Die evaporierte Wassermenge aus den einzelnen Schichten z_j oberhalb der maximalen Tiefe _ZE wird durch das Minimum aus dem verbleibenden Evaporationsbedarf und dem verfügbaren Bodenwasser oberhalb des Welkepunktes limitiert. Letzteres wird durch die Reduktionsfunktion R_SE(z) realisiert:

Es_p potenzielle Bodenevaporation [mm]
R_SE Reduktionsfunktion
ΘWassergehalt [mm]
Θ_WP Wassergehalt am Permanentwelkepunkt [mm]
Θ_FK Wassergehalt bei Feldkapazität [mm]

Der Transpirationsanteil E_TR(t) wird in seiner Menge und vertikalen Verteilung durch das Vegetationsmodell in Abhängigkeit vom pflanzenverfügbaren Wasser berechnet. Im generischen Vegetationsmodell erfolgt die Berechnung von E_TR(t) nach Gleichung 31 (Kapitel 4.4.3).

Liegt der aktuelle Bodenwassergehalt nach Realisierung von Entzug und Perkolation immer noch oberhalb der Feldkapazität, so kann auch ein horizontal gerichteter Abfluss P^h in Abhängigkeit von der Schichtneigung erfolgen, wenn die hydraulische Leitfähigkeit der unterliegenden Schicht kleiner als in der betrachteten Schicht ist, deren Matrix-Porenvolumen sehr klein ist (z.B. bei hohem Skelettanteil) bzw. deren gesamter Porenraum wassergesättigt ist.

α Neigung der Bodenschicht [°]
λ Leitfähigkeitsparameter nach Koitzsch [-]
BT Bodentemperatur [° C]
k^h Korrekturfaktor für den hypodermischen Abfluss [-]

Melioration hat einen wesentlichen Einfluss auf den Gebietswasserhaushalt, aber auch auf den Stoffaustrag (insbesondere Stickstoff und Phosphor) aus landwirtschaftlich genutzten Flächen. Genaue Informationen zu meliorierten Flächen, Art und Dauer der Meliorationsmaßnahmen liegen jedoch zumeist nicht flächendeckend vor. Alternativ zu einer genauen Kartierung der entwässerten Flächen können potenziell drainierte Flächen berücksichtigt werden. Voraussetzung ist eine Zuordnung der einzelnen Standorte zu einem Drainagetyp (Graben oder Rohrdrainage) und der Entwässerungswahrscheinlichkeit (Anschlussgrad). Entscheidend bei der Modellierung ist die Drainagetiefe. Entsprechend den Literaturangaben wird allgemein von einer Tiefe von 1 m bzw. 40 cm bei landwirtschaftlich genutzten Niedermoorstandorten ausgegangen (s. Möller et al., 2009).

Bei Sättigung (Feldkapazität) der Bodenschicht in dieser Tiefe wird jedes Zusatzwasser sofort dem Drainagesystem zugeführt und im gleichen Zeitschritt abflusswirksam. Liegt die Entwässerungswahrscheinlichkeit unter 100 %, so wird nur der entsprechende Anteil der Bodenlösung über Feldkapazität dem Vorfluter zugeführt. Der Rest steht weiterhin für die Verdunstung und Tiefenversickerung zur Verfügung.

Wichtig für die Abbildung der Entwässerung ist die Konsistenz der Eingangsdatenbasis zur Landnutzung, den Meliorationsinformationen und den mittleren Grundwasserflurabständen.

Der kapillare Aufstieg aus gesättigten tieferen Bodenschichten in den Wurzelraum kann für die Wasserversorgung der Vegetation von Bedeutung sein. Entsprechend des Kapazitätsansatzes zur Beschreibung der Bodenwasserdynamik kann zur Simulation dieses Phänomens nur ein empirischer Ansatz über Tabellenfunktionen, wie sie in der Bodenkundlichen Kartieranleitung (Ad-hoc-AG Boden, 2005) angegeben sind, genutzt werden. Kapillaraufstieg in die aktuelle Rechenschicht wird berechnet, wenn zwei Bedingungen erfüllt sind:

• Der aktuelle Bodenwassergehalt der Rechenschicht liegt unterhalb eines Schwellenwertes, der mit (Θ_WP -Θ_FK)/2 angenommen wird.
• Die Schicht liegt im Bereich der für dieses Bodenprofil nach der Bodenkundlichen Kartieranleitung berechneten Distanz, die das Grundwasser bei gegebener Wasserspannung im Boden gegen die Schwerkraft überwinden kann (kapillare Aufstiegshöhe).

Sind diese Bedingungen erfüllt, so wird eine Auffüllung des Schichtwassergehaltes maximal bis zu diesem Grenzwert entsprechend der kapillaren Aufstiegsrate simuliert. Die kapillaren Aufstiegsraten und -höhen werden in der Bodenkundlichen Kartieranleitung bodenartspezifisch angegeben. Für die Übertragung dieser Angaben auf geschichtete Böden wird für die kapillare Aufstiegshöhe die Bodenart der an den Grundwasserspiegel angrenzenden Bodenschicht berücksichtigt. Die kapillare Aufstiegsrate ergibt sich als das Minimum der jeweiligen Aufstiegsraten der Bodenschichten zwischen der Zehrschicht und dem Grundwasserspiegel.