Wie wird die NatUrWB-Referenz bestimmt?

Um zu verstehen, was hinter der bestimmten naturnahen urbanen Wasser-Bilanz (NatUrWB) Referenz steht, ist im Folgendedien Vorgehensweise erläutert. Für ein Planungsgebiet in Siedlungen wird als Referenzzustand des Wasserhaushaltes der Wasserhaushalt des gleichen Gebietes mit nicht urbaner Landnutzung definiert. Hierbei wird die umliegende heutige Landnutzung als Kulturland ohne Siedlungsanteile auf gleichen Böden in ähnlicher Lage angenommen. Diese Herangehensweise wurde von einer Arbeitsgruppe im Rahmen des Projektes WaSiG entwickelt1.
Grundsätzlich basieren die Werte auf Simulationen mit dem RoGeR_WB_1D Wasserbilanz-Modell der Uni Freiburg. Hierfür musste ein Set aus Wetterdaten und Bodenprofilen erstellt werden. Im Folgenden werden die eingehenden Daten näher erläutert.

Wetterdaten

Um einzelne Flächen zu simulieren, wurde zuerst aus den Wetter-Stationen des Deutschen Wetterdienstes (DWD) die Stationen mit ausreichend und qualitativ guten Daten für die Wasserbilanz-Simulationen ausgewählt. Diese Selektion ergab 754 mögliche Stationen um den Zeitraum vom 01.01.2009 bis zum 31.12.2019 mit 10-minütigen Niederschlagsdaten simulieren zu können.
Bestehende Datenlücken einzelner Stationen wurden anschließend durch eine Regionalisierung von benachbarten Stationen aufgefüllt. Hierfür wurden die regionalisierten, langjährigen Rasterdaten vom DWD, die den Zeitraum von 1991 bis 2020 berücksichtigen, herangezogen. Daraus wurden die jährlichen Wichtungsfaktoren für die Temperatur- und die Evapotranspirationsdaten zwischen den Stationen bestimmt. Für den Niederschlag wurde hierfür ein halbjährlicher Wichtungsfaktor für Sommer (April-September) und Winter (Oktober-März) bestimmt.
Die gefüllten Niederschlag-Zeitreihen werden anschließend noch nach Richter (1995) korrigiert.

Anschließend wurden die Wetterdaten auf die langjährigen Stationsmittelwerte normiert, um die jährlichen Schwankungen auszugleichen. Dies ist in den folgenden zwei Diagrammen näher erläutert. Hier wurde jeweils für die Temperaturen bzw. den Niederschlag, der Unterschied zum langjährigen Mittelwert dargestellt. Dabei wird der Niederschlag als Quotient des gemessenen jährlichen Mittelwertes zum langjährigen (1991-2020) Mittelwert angezeigt. Liegt der Quotient bei 1, entspricht das gewählte Jahr einem durchschnittlichen Klimajahr. Ist der Quotient über eins, handelt es sich um ein niederschlagsreicheres Jahr als der langjährige Durchschnitt. Bei der Temperatur wird die absolute Differenz zum langjährigen Mittelwert angezeigt.
Des Weiteren ist der Mittelwert dieses Quotienten, bzw. Differenz, über die 11 Jahre Simulationsdauer als gestrichelte Linie dargestellt. Durch die orangenen Linien lässt sich erkennen, dass der Zeitraum von 2009 bis 2019 trockener und wärmer als der langjährige Mittelwert ist. Die blauen Linien zeigen den Verlauf der korrigierten Wetterdaten an. Dieser liegt für den Niederschlag im Mittel bei eins und entspricht daher dem langjährigen Mittel.


Ein Diagramm des Verhältnisses des jährlichen Mittelwertes zum langjährigen Mittelwert des Niederschlags
Ein Diagramm des Verhältnisses des jährlichen Mittelwertes zum langjährigen Mittelwert der Temperatur

Die Wetterdaten weisen also weiterhin eine jährliche Schwankung auf, entsprechen aber über elf Jahre gemittelt dem langjährigen Klimamittel. Folglich sind die erhaltenen NatUrWB-Referenzwerte auch als langjährige Mittelwerte anzusehen.

Bodenprofile

Um die Böden für die Simulationen der Wasserbilanz zu parametrisieren, wurde die Bodenübersichtskarte (BÜK) der Bundesanstalt für Geologie und Rohstoffe herangezogen. In dieser sind für ganz Deutschland die Böden in Bodengesellschaften gruppiert und ihr geografisches Auftreten definiert. Jede Bodengesellschaft besteht aus mehreren Bodenprofilen mit den jeweiligen Bodenhorizonten.
Die Ableitung der RoGeR Parameter aus der BÜK ist in Anlehnung an Steinbrich et al. (2018)51/a> geschehen. Hierbei wurde für jeden Bodenhorizont mithilfe der Pedotransferfunktionen von Wessolek et al. (2009)2 "die gesättigte hydraulische Leitfähigkeit (Ksat), die nutzbare Feldkapazität (nFK), die Luftkapazität (LK), sowie die Funktion des Horizontes als Basis des Feinbodens (Gestein, Stauhorizont, Grundwasserbeeinflussung) ermittelt."
Diese Horizontparameter wurden anschließend zu einem einzigen fiktiven Horizont zusammengefasst, da RoGeR_1D_WBM nur mit einem Horizont rechnet. Hierfür wurde die Bodenmächtigkeit bis zum ersten Gesteins-Horizont oder dem ersten grundwasserbeeinflussten Horizont gewählt. Für die hydraulische Leitfähigkeit der Tiefenperkolation wurde die minimal vorkommende hydraulische Leitfähigkeit der Horizonte gewählt, da der Horizont mit der geringsten Leitfähigkeit ausschlaggebend ist, wieviel Wasser versickern kann. Ansonsten wurde das gewichtete arithmetische Mittel der Parameter bestimmt.

Da die BÜK nicht lückenlos ist und nicht alle Böden parametrisiert sind, wurde diese mit der geologischen Übersichtskarte verschnitten. Anschließend wurden die definierten Bodengesellschaften ohne Hauptbodenprofil (nur das Profil mit dem größten Flächenanteil) durch Nachbar-Polygone ersetzt. Hierfür wurden zuerst versucht die Bodengesellschaft eines Nachbar-Polygon gleicher Geologie zu nutzen. Bei mehreren zutreffenden Polygonen wurde der Nachbar mit der größten gemeinsamen Grenze herangezogen. Wurde kein Nachbar mit gleicher Geologie ausfindig gemacht, wurde das Nachbar-Polygon mit der größten gemeinsamen Grenze genutzt. Dieser Prozess verschlechtert das Ergebnis, war aber notwendig ist, um in manchen Gebieten eine naturnahe Wasserbilanz zu berechnen. Daher wird der Anteil dieses Schrittes als Warnmeldung dem Ergebnis hinzugefügt.

Modellgebiete

Um die Polygone zu definieren, für die eine Wasserbilanz-Simulation gemacht wurde, wurde zunächst die im Abschnitt davor gefüllte BÜK herangezogen (1.). Zunächst wurden Wasser-Polygone (Seen & Flüsse) ausgenommen, da diese sich durch die Urbanisierung nicht ändern außer sie würden überbaut werden, was der Datengrundlage aber nicht zu entnehmen wäre. In diesem Fall wäre das NatUrWB-Modell nicht anwendbar, wobei die Bestimmung der Wasserbilanzgrößen recht einfach manuell zu bestimmen ist, da es zu 100 % der potenziellen Evapotranspiration kommen würde und der Rest zum Oberflächenabfluss zu zählen wäre.

Der Niederschlag ist sehr stark von der Topographie abhängig und die Evapotranspiration stark von der Exposition. Weiterhin hat das Gefälle einen großen Einfluss auf den modellierten Oberflächen- und Zwischenabfluss. Da für die Modellierung jeweils nur ein Set an Klima-Wichtungsparametern und ein Gefälle angesetzt wird, sollten die definierten Modellgebiete daher möglichst wenig Variabilität dieser Parameter aufweisen. Aus diesem Grund wurden die Polygone der BÜK250 weiter aufgeteilt um möglichst homogene Gebiete zu beschreiben.

Hierfür wurde zuerst aus einem Digitalen Geländemodel5, das auf eine Auflösung von 20 m skaliert wurde, das Gefälle und die gesamt Strahlungsbilanz im Verhältnis zu einem flachen Gelände (Expositionsfaktor) bestimmt. Je Raster wurde folgender Arbeitsablauf durchlaufen:

  1. Die Standardabweichung des Parameters je Polygon wurden mit zonaler Statistik bestimmt.
  2. Polygone die über einem gewissen Grenzwert (siehe Tabelle) liegen wurden ausgewählt.
  3. Je Polygon wurden Klassengrenzen bestimmt, die das jeweilige darunterliegende Raster in gleich große Gruppen von maximal X Einheiten (siehe Tabelle) einteilt.
  4. Aus dem jeweiligen Raster wurden die Konturlinien dieser Klassengrenzen ermittelt
  5. Das Polygon wurde mittels dieser Konturlinien aufgeteilt.
  6. Damit nicht zu viele kleine Polygone entstehen, wurden Polygone die kleiner als 80 ha sind vereinfacht indem Sie dem Nachbarpolygon zugeordnet wurden.

Je nach Raster wurden hierfür andere Grenzwerte genutzt die in der folgenden Tabelle aufgelistet sind.

Liste der genutzten Grenzen der Standardabweichungen und Gruppeneinteilung je Raster.
# Raster Grenze der Standardabweichung Einteilung in Gruppen von maximal
1 Topografie 50 m 100 m
2 Gefälle 10 ° 15 °
3 Expositionsfaktor 15 % 20 %

Die folgende Abbildung stellt diesen Prozess Beispielhaft für ein Polygon dar:

beispielhafte Darstellung der Vorgehensweise zur Aufteilung der Polygone in Abhängigkeit des Rasters

Weiterhin werden Polygone die größer als 1 km2 sind mit den Thiessen-Polygonen der Wetterstationen verschnitten, damit eine gute Zuordnung der Polygone zu den Wetterstationen möglich ist. Anschließend wird jedem Polygon die jeweils nächst gelegene Wetterstation zugeordnet.

Die sich so ergebenden Polygone werden im folgenden Modellgebiete genannt.

Regionalisierung und Gefälle je Modellgebiet

Die Wichtungsfaktoren für den Niederschlag, die Temperatur und die Evapotranspiration wurden für jedes der erstellten Modellgebiete bestimmt. Für Temperatur und Evapotranspiration wurde eine jährliche Wichtung angewendet und für den Niederschlag wurde eine Sommer- (April-September) und ein Winter-Wichtungsparameter (Oktober-März) bestimmt. Hierfür wurden die langjährigen regionalisierten 1 x 1 km Raster des DWDs der Jahre 1991-2020 genutzt. Mittels zonaler Statistik wurde für den Rastermittelpunkt jedes Polygons das jeweilige langjährige Mittel bestimmt und mit dem langjährigen Stationsmittelwert verglichen.

Für die Temperatur und die Evapotranspiration wurde zusätzlich zu diesem Höhenfaktor noch ein Expositionsfaktor bestimmt. Hierfür wurde die Sonneneinstrahlung aus dem DGM55 und einem fiktiven flachen DGM bestimmt. Aus diesen

Das Gefälle wurde anhand des digitalen Geländemodells mit in einer Auflösung von 5 x 5 m5, das auf 20 m hochskaliert wurde, bestimmt. Anschließend wurde ebenfalls mit einer zonalen Statistik der Mittelwert der sich im Modellgebiet befindenden Rastermittelpunkte bestimmt.

Landnutzung

Für jedes dieser Modellgebiete wurden 13 verschiedenen Landnutzungstypen definiert. Hierüber wurden die Makroporendichte in Abhängigkeit der Landnutzung, die Verdunstungstiefe und der jahreszeitlich variierende Interzeptionsspeicher definiert.
Die definierte Landnutzungstype sind:

  1. Felsen ohne Vegetation
  2. Flächen mit spärlicher Vegetation
  3. Strände- Dünen und Sandflächen
  4. Heiden und Moorheiden / Wald-Strauch-Übergangsstadien
  5. Nicht bewässertes Ackerland
  6. Weinbauflächen
  7. Obst- und Beerenobstbestände
  8. Wiesen und Weiden / Natürliches Grünland
  9. Landwirtschaftlich genutztes Land mit Flächen natürlicher Bodenbedeckung von signifikanter Größe
  10. Laubwälder
  11. Mischwälder
  12. Nadelwälder
  13. Sümpfe / Torfmoore

Abfrage

Bei der Abfrage der Ergebnisse wird eine als naturnah angenommene Landnutzungsverteilung ermittelt. Hierfür wurden die Modellgebiete mit den Naturraumeinheiten des Hydrologischen Atlasses Deutschlands3 und dem Corine Land Cover Layer44 verschnitten. Für jedes Verschnitt-Polygon wird anschließend in der gleichen Naturraumeinheit nach den nicht urbanen Landnutzungen gesucht. Deren Verteilung wird als die naturnahe Verteilung für dieses Polygon angesehen und die verschiedenen Modellergebnisse werden hiermit gewichtet gemittelt.

Um auf die drei Hauptkomponenten des Wasserhaushaltes (Evapotranspiration, Grundwasserneubildung und Abfluss) zu kommen, musste der durch das RoGeR_WB_1D Modell bestimmte Zwischenabfluss aufgeteilt werden. Dieser Anteil beschreibt das Wasser, das in der ungesättigten Bodenzone horizontal abfließt. Je nach Gegebenheit fließt dieses Wasser anschließend dem Grundwasser oder einem nahe gelegenen Oberflächengewässer zu. Bei einem hohen Grundwasserspiegel ist davon auszugehen, dass der Zwischenabfluss zu einer schnellen Abflussreaktion in einem nahegelegenen Gewässer führt. Daher wurde das Kriterium so gewählt, dass der Zwischenabfluss dem Direktabfluss zugezählt wird, wenn der GW-Spiegel im Bodenprofil ansteht. Liegt das Grundwasser tiefer als die definierte Bodensohle, wird der Zwischenabfluss der Grundwasserneubildung hinzugezählt.

Abschließend werden die verschiedenen Ergebnisse aller Polygone je nach Anteil im Untersuchungsgebiet gewichtet gemittelt, um auf eine NatUrWB-Referenzwert zu kommen.


  1. Steinbrich, Andreas; Henrichs, Malte; Leistert, Hannes; Scherer, Isabel; Schuetz, Tobias; Uhl, Mathias; Weiler, Markus (2018): Ermittlung eines naturnahen Wasserhaushalts als Planungsziel für Siedlungen. In: Hydrologie und Wasserbewirtschaftung 62 (6), S. 400–409. DOI: 10.5675/HYWA_2018.6_3. Online Verfügbar unter http://www.hywa-online.de/ermittlung-eines-naturnahen-wasserhaushalts-als-planungsziel-fuer-siedlungen-2/
  2. Wessolek, Gerd; Kaupenjohann, Martin; Renger, Manfred (2009): Bodenphysikalische Kennwerte und Berechnungsverfahren für die Praxis. In: Bodenökologie und Bodengenese (40). Online verfügbar unter https://www.researchgate.net/profile/Gerd_Wessolek/publication/294427537_Bodenphysikalische_Kennwerte_und_Berechnungsverfahren_fur_die_Praxis/links/56c0c55c08ae2f498ef99662.pdf, zuletzt geprüft am 10.02.2021.
  3. Bundesministerium für Umwelt Naturschutz und Reaktorsicherheit (Hg.) (2003): Hydrologischer Atlas von Deutschland (HAD). 3. Lieferung. Bonn.
  4. Bundesamt für Kartographie und Geodäsie (2018): Corine Land Cover 5 ha. CLC5.
  5. GeoBasis-DE: Bundesamt für Kartographie und Geodäsie. 2016. Digitales Geländemodell Gitterweite 5 m. DGM5. Produktstand: 2016. Online verfügbar unter https://gdz.bkg.bund.de/index.php/default/digitale-geodaten/digitale-gelandemodelle/digitales-gelandemodell-gitterweite-5-m-dgm5.html, zuletzt geprüft am 08.11.2021.
© Max Schmit (2021)