Aerosol-Wolken-Wechselwirkungen stellen einen der wesentlichen Unsicherheitsfaktoren bei Verständnis und Quantifizierung der geographischen Verteilung von Wolken- und Niederschlagseigenschaften, aber auch des Strahlungsantriebs des globalen Klimawandels dar. Die grundlegende Idee des Projekts ist es, regional unterschiedliche Trends in anthropogenenen Emissionen von Aerosolen zu nutzen, um deren Einfluss auf Trends in Wolken-, Niederschlags- und Strahlungsgrößen zu bestimmen. Hierzu sollen verschiedene Szenarien in Multi-Klimamodell-Ensembles ('historische' Simulationen mit allen Strahlungsantrieben und 'Aerosol'-Simulationen mit allen Antrieben außer anthropogenem Aerosol) analysiert werden und mit Beobachtungsdaten verglichen werden. Konkret werden vier Fragen untersucht:(i) Welche Beziehung besteht zwischen regionalen Trends in Aerosolemissionen und Wolken-Strahlungs-Effekten? - Diese Studien analysieren Simulationen aus dem Multi-Modell-Ensemble.(ii) Wie erfolgreich reproduzieren die Modelle beobachtete Trends? Hier werden die Klimamodelle mit Beobachtungsdaten verglichen.(iii) Welchen Einfluss haben Emissionstrends für Aerosole und resultierende Strahlungsantriebe auf die atmosphärische Zirkulation? Simulationen mit dem Aerosol-Klima-Modell ECHAM6-HAM2 sollen für drei Zeitscheiben durchgeführt und analysiert werden.(iv) Welche Rolle spielen Emissionstrends für Änderungen in Extremniederschlägen in Südost-Asien? - Mit speziellen Simulationen sollen die verschiedenen Hypothesen getestet werden.
Der weit nach Süden vordringende Keil Südamerikas ist weltweit die einzige nennenswerte Landmasse zwischen ca. 45° und 60°Süd. Das senkrecht zur Hauptwindrichtung verlaufende Andengebirge stellt eine wirksame Barriere für die Westwinddrift dar und hat einen bestimmenden Einfluss auf die hemisphärische Zirkulation sowie das lokale Wettergeschehen. Das Gebirge zwingt die maritimen Luftmassen zum Aufsteigen, was häufig mit intensiven Steigungsregen auf der Luvseite der Anden einhergeht. Durch die Überströmung des Gebirges kommt es zur Ausbildung von speziellen Prozessgefügen in der atmosphärischen Strömung sowohl auf der Meso- als auch regionaler Skala. Der damit einhergehende Transport und Austausch von Energie und Masse beeinflusst maßgeblich die Entstehung und den Ausfall von Hydrometeoren. Trotz der starken Wechselwirkung zwischen Strömung, Topographie und Niederschlag wurde in Patagonien darüber bisher nur wenig geforscht. Das vorgeschlagene Forschungsvorhaben leistet daher einen Beitrag zum Verständnis der Wechselwirkung zwischen dynamischen Prozessen und der räumlichen und zeitlichen Variabilität von Niederschlag in dieser Region. Ziel des Projektes ist die Quantifizierung wichtiger Prozesse die neue Aufschlüsse über die relevanten Mechanismen liefern soll. Anhand von hochauflösenden numerischen Simulationen werden an Einzelfallstudien die dynamischen und thermodynamischen Eigenschaften der atmosphärischen Strömung im Detail analysiert. Begleitende Sensitivitätsstudien mit vereinfachten analytischen Modelle werden zudem Aussagen zu den Auswirkungen der atmosphärischen Variabilität auf die Niederschlagsverteilung liefern. Das aus der Studie gewonnene Prozessverständnis ist eine wichtige Grundlage für weiterführende Forschungsarbeiten im Bereich der Hydrologie, Glaziologie und Ökologie.
Das Thermosphären/Ionosphären (T/I) System wird sowohl von oben (solar, geomagnetisch), als auch von unten stark beeinflusst. Einer der wichtigsten Einflüsse von unten sind Wellen (z.B. planetare Wellen, Gezeiten, oder Schwerewellen), die größtenteils in der Troposphäre bzw. an der Tropopause angeregt werden. Die vertikale Ausbreitung der Wellen bewirkt hierbei eine vertikale Kopplung der T/I mit der unteren und mittleren Atmosphäre. Vor allem der Einfluss von Schwerewellen (GW) ist hierbei weitestgehend unverstanden. Einer der Gründe hierfür ist, dass GW sehr kleinskalig sind (einige zehn bis zu wenigen tausend km) - eine Herausforderung, sowohl für Beobachtungen, als auch für Modelle. Wir werden GW Verteilungen in der T/I aus verschiedenen in situ Satelliten-Datensätzen ableiten (z.B., sowohl in Neutral-, als auch in Elektronendichten). Hierfür werden Datensätze der Satelliten(-konstellationen) SWARM, CHAMP, GOCE und GRACE verwendet werden. Es sollen charakteristische globale Verteilungen bestimmt, und die wichtigsten zeitlichen Variationen (z.B. Jahresgang, Halbjahresgang und solarer Zyklus) untersucht werden. Diese GW Verteilungen werden dann mit von den Satelliteninstrumenten HIRDLS und SABER gemessenen Datensätzen (GW Varianzen, GW Impulsflüssen und Windbeschleunigungen durch GW) in der Stratosphäre und Mesosphäre verglichen. Einige Datensätze (CHAMP, GRACE, SABER) sind mehr als 10 Jahre lang. Räumliche und zeitliche Korrelationen zwischen den GW Verteilungen in der T/I (250-500km Höhe) und den GW Verteilungen in der mittleren Atmosphäre (Stratosphäre und Mesosphäre) für den gesamten Höhenbereich 20-100km werden untersucht werden. Diese Korrelationen sollen Aufschluss darüber geben, welche Höhenbereiche und Regionen in der mittleren Atmosphäre den stärksten Einfluss auf die GW Verteilung in der T/I haben. Insbesondere Windbeschleunigungen durch GW, beobachtet von HIRDLS und SABER, können zusätzliche Hinweise darauf geben, ob Sekundär-GW, die mutmaßlich in Gebieten starker GW Dissipation angeregt werden, in entscheidendem Maße zur globalen GW Verteilung in der T/I beitragen. Zusätzlich wird der Versuch unternommen, sowohl GW Impulsfluss, als auch Windbeschleunigungen durch GW aus den Messungen in der T/I abzuleiten. Solche Datensätze sind von besonderem Interesse für einen direkten Vergleich mit von globalen Zirkulationsmodellen simulierten GW Verteilungen in der T/I. Diese werden für eine konsistente Simulation der T/I in Zirkulationsmodellen (GCM) benötigt, stellen dort aber auch eine Hauptunsicherheit dar, da eine Validierung der modellierten GW durch Messungen fehlt.
Der Klimawandel ist eine der Hauptherausforderungen für die Menschheit im 21. Jahrhundert. Seine Auswirkungen sind vielschichtig wobei der anwachsende Massenverlust von Gletschern außerhalb der großen Eisschilde sowie deren bedeutender Beitrag zum Meeresspiegelanstieg zu den am stärksten hervorstechenden zählt. Diesbezüglich sind die Gletscher und Eiskappen der Arktis aufgrund ihres großen Volumens und ihrer großen Oberfläche, die als Kontaktfläche zum Klima- und Ozeanantrieb und damit zum Klimawandel selber fungiert, von besonderer Bedeutung. Da die Arktis darüber hinaus diejenige Region der Erde mit dem höchsten, prognostizierten, zukünftigen Temperaturanstieg ist, wird erwartet, daß sich die Bedeutung der arktischen Eismassen für den Meeresspiegelanstieg auch in Zukunft fortsetzt oder sogar noch steigern wird.Die großen Gletscher der Nordpolarregion umgeben den arktischen Ozean in ähnlichen Breitenlagen, weisen aber in jüngster Zeit ein inhomogenes Verhalten auf. Diese Tatsache legt eine räumliche Variabilität der klimatischen und ozeanischen Antriebsmechanismen der Gletschermassenbilanz innerhalb der zirkumarktischen Regionen nahe und offenbart damit die Diversität der Einflüsse des Klimawandels. Bezüglich der Variabilität der Antriebsmechanismen weist Svalbard in der Arktis eine einzigartige Lage auf. Es liegt an der Grenze zwischen kalten, polaren Luftmassen und Ozeanwassern und den Einflüssen des Westspitzbergenstroms, welcher der hauptsächliche Warmwasserlieferant für das arktische Umweltsystem ist. Darum verspricht das Erforschen der Reaktionen der Gletscher auf Svalbard auf die Veränderlichkeit des Klima- und Ozeanantriebs bedeutende Einblicke in die komplexe Kausalkette zwischen Klimawandel, der Variabilität der Klima- und Ozeanbedingungen in der Arktis und der Reaktion der arktischen Landeismassen. Das Ziel des Projektes ist es eine zuverlässige Abschätzung der räumlichen und zeitlichen Variabilität der klimatischen Massenbilanz aller Gletscher und Eiskappen auf Svalbard zu erreichen und diese mit dem Klima- und Ozeanantrieb in Verbindung zu setzen. Dazu wird ein räumlich verteiltes, von statistisch downgescalten Klimadaten angetriebenes Model zur Berechnung der klimatischen Massenbilanz aufgesetzt. Die Massenbilanz aller Gletscherflächen auf Svalbard wird für den Zeitraum 1948-2013 modelliert und die zeitlich variablen Felder von Ablation, Akkumulation, wiedergefrorenem Schmelzwasser und klimatischer Massenbilanz für anschließende geostatistische Studien genutzt. Diese Studien werden potentielle Einflüsse der raumzeitlichen Variabilität von großräumigen Mustern des Luftdrucks, der Meereisbedeckung und der Meeresoberflächentemperatur auf die Variabilität der Gletschermassenbilanz auf Svalbard identifizieren und analysieren. Auch Telekonnektionen zu fernen Modi der atmosphärischen Zirkulation werden durch Studien bezüglich der potentiellen Einflüsse verschiedener atmosphärischer Zirkulationsindizes in die Betrachtungen einbezogen.
Interne Schwerewellen (SW) verbinden verschiedene Schichten der Atmosphäre von der Troposphäre bis zur Thermosphäre und treiben die großskalige Zirkulation der mittleren Atmosphäre an. Viele der für SW relevanten Prozesse, von ihrer Entstehung über die Ausbreitung bis zur Dissipation sind jedoch unvollständig verstanden und, wegen der geringen typischen Wellenlänge, meist schlecht in numerischen Wettervorhersage- und Klimamodellen repräsentiert. GWING ist eines der Projekte der Forschergruppe MS-GWaves, die darauf abzielt, unser Verständnis der oben angesprochenen multi-skalaren dynamischen Schwerewellenprozesse zu verbessern, um letztendlich eine einheitliche Parametrisierung der in Atmosphärenmodellen nicht auflösbaren Schwerewellen (und ihrer Effekte) von der Entstehung bis zur Dissipation zu entwickeln. Um hierzu beizutragen, ist das zentrale Ziel von GWING die Entwicklung und Anwendung des atmosphärischen Zirkulationsmodells UA-ICON. Mit diesem Modell integriert GWING das in der Forschergruppe MS-GWaves entwickelte Wissen. In der zweiten Phase von GWING stehen zwei übergeordnete wissenschaftliche Fragen im Fokus: a) Welche Bedeutung haben Eigenschaften von Schwerewellen, die in klassischen Parametrisierungen nicht berücksichtigt werden, also insbesondere horizontale und nicht-inständige Propagation sowie die Wechselwirkung transienter Wellen mit dem Grundstrom? b) Welche Rolle spielen Schwerewellen für die globale Zirkulation und ihre Variabilität? Um diese Fragen zu beantworten, werden wir UA-ICON global sowohl mit einer Maschenweite von etwa 20 km (d.h. mit Auflösung von SW bis etwa 100 km Wellenlänge) als auch mit grober Auflösung, dafür aber mit der State-of-the-art Parametrisierung MS-GWaM nutzen. Weiterhin werden spezielle Beobachtungsepisoden mit sehr hoch (ca. 1,5 km) aufgelösten Nestern simuliert. Zur Evaluation und Analyse werden diese Modellsimulationen mit Beobachtungen der Partnerprojekte zusammengeführt. Die wesentlichen Entwicklungsziele für UA-ICON in Phase 2 des Projekts sind dementsprechend die Implementierung von MS-GWaM (entwickelt im Partnerprojekt 3DMSD), die Einführung physik-basierter Schwerewellenquellen (zusammen mit 3DMSD und SV) und eine verbesserte Behandlung von SW bei sehr hoher Modellauflösung. Die Nutzung der verschiedenen UA-ICON-Konfigurationen wird schließlich erlauben, die Bedeutung bisher vernachlässigter Eigenschaften von SW zu untersuchen, d.h. die erste der oben genannten Fragestellungen zu beantworten. Ein spezielles Ziel im Rahmen von GWING ist diese Untersuchung für Episoden plötzlicher Stratosphärenerwärmungen, die durch sich schnell ändernde und zonal nicht symmetrische Bedingungen des Grundstroms gekennzeichnet sind. Im Hinblick auf die zweite übergeordnete Fragestellung, wird sich GWING auf a) die Rolle der SW und einer hohen Modellausdehnung für die Simulation von Zirkulationsänderungen bei globaler Erwärmung und b) die Rolle für die Güte von Wettervorhersagen konzentrieren.
Kopplungsprozesse zwischen der Ionosphäre und der neutralen Atmosphäre spielen eine wichtige Rolle für die dynamischen Prozesse in der oberen Atmosphäre. Neue Fortschritte im Verständnis dieser Prozesse wurden erreicht seitdem Satelliten im erdnahen Orbit kontinuierlich hochgenaue Daten der thermosphärischen und ionosphärischen Parameter (z.B. Massendichten, zonale Winde und Elektronendichteprofile) bereitstellen. Mit diesem Projekt planen wir die Beobachtung der Auftretenshäufigkeit und Eigenschaften sporadischer E Schichten auf globaler Skala. Die Untersuchungen basieren auf GPS Radiookkultationen der Satelliten CHAMP, GRACE, TerraSAR-X, TanDEM-X und FORMOSAT-3/COSMIC. Seit dem Start des Satelliten CHAMP im Jahre 2001 wurden mehr als 5 Millionen der Radiookkultationsprofile aufgezeichnet, was ermöglicht, dass das Auftreten und die Eigenschaften der sporadischen E Schichten in hoher räumlicher Auflösung analysiert werden können. Weiterhin ermöglicht die Zeitreihe erste statistische Trendanalysen der genannten Parameter. Während der Durchführung des Projektes soll der momentan genutzt numerischer Algorithmus zur Detektion von sporadischen E Schichten um ein Modul erweitert werden, der ermöglichen wird auch Rückschlüsse auf die Eigenschaften der Schichten zu ziehen. Globale Beobachtungen der Intensitäten sporadischer E Schichten existieren aktuell nicht und werden von uns zum erstmalig bereitgestellt werden. Diese Datenbasis kann genutzt werden, um statistische Änderungen im Verhalten der sporadischen E Schichten zu Untersuchen. Ebenfalls werden wir untersuchen, ob Abhängigkeit der sporadische Eigenschaften von anderen geophysikalischen Parametern, wie beispielsweise die Abnahme des Erdmagnetfeldes, der Solarzyklus, atmosphärische Gezeiten, Meteoreinfall oder Plamadichteabnahmen in der Ionosphäre zu finden sind.
Schwerewellen (GWs) sind zu kleinskalig, um in den heutigen Wetter- und Klimamodellen aufgelöst zu werden. Sie müssen daher parametrisiert werden, da sie einen starken Einfluss auf die Dynamik der großen Skalen haben. Parametrisierungen existieren für orographisch und konvektiv erzeugte GWs, während für die GW-Quellen entlang großskaliger Jets noch keine etablierte Parametrisierung vorliegt. Die Quellen resultieren aus einer spontanen Imbalance (SI) der großskaligen quasi-geostrophischen Strömung. Die Untersuchung von Schwerewellenabstrahlung durch SI ist schwierig, da die GWs in ein sehr komplexes zeitabhängiges Strömungsfeld eingebettet sind, mit einer großen Zahl von interagierenden Prozessen. Auch die Validierung von Parametrisierungen wird dadurch erschwert. Daher kombinieren wir Theorie und numerische Modellierung mit ergänzenden Laborexperimenten. Laborexperimente garantieren eine Reproduzierbarkeit der betrachteten großskaligen Strömungssituation. Die direkte Korrespondenz zwischen den experimentellen Daten und den Modelldaten und die erwähnte Reproduzierbarkeit machen das Laborexperiment zu einem idealen Prüfstand für Parametrisierungen und für die Untersuchung klimarelevante Prozesse. Das differenziell beheizte rotierende Zylinderspalt-Experiment, welches an der BTU (Brandenburg Technische Universität Cottbus-Senftenberg) aufgebaut und betrieben wird, stellt die Referenzdaten für Benchmark-Simulationen an der GU-F (Goethe Universität Frankfurt) und dem IAP (Leibniz Institut für Atmosphärische Physik, Kühlungsborn) bereit. Dabei stehen Experimente im Vordergrund, die zeigen sollen, welche baroklinen Strömungen eine besonders ausgeprägte GW-Abstrahlung aufweisen. Ergänzend dazu werden idealisierte numerische Simulationen an der GU-F und dem IAP durchgeführt, um die Variabilität der GWs und den Abstrahlungsprozess zu untersuchen. Wichtig ist dabei, einen Zusammenhang zwischen verschiedenen großskaligen Strömungen und der mesoskaligen GW-Quelle herzustellen und diesen Zusammenhang mittels grob aufgelöster Wellenstrahlenmodelle zu validieren. Ziel ist es, eine skalenabhängige SI-Parametrisierung zu konstruieren. Diese Parametrisierung soll mit Hilfe der Labor-Referenzdaten validiert werden. Begleitet wird dies von einer Analyse grob- und feinaufgelöster Daten aus UA-ICON Simulationen. Schließlich soll die Parametrisierung an das Wellenstrahlenmodell MS-GWaM angekoppelt werden, welches in UA-ICON implementiert ist.
Die Erforschung der Land-Atmosphären-Rückkopplung für heterogene Landoberflächen benötigt die Erfassung biophysikalisch bedeutsamer Wärme- und Impulsbilanzen nahe der Unterlage, die das Projekt motivieren. Dieser Energie- und Massenaustausch ist eine essenzielle Komponente in der Prozesskette der Rückkopplungen. Trotz beachtlicher Fortschritte im Verständnis und Quantifizierung der Auswirkungen einfacher Heterogenitätssprünge bei orthogonaler und paralleler Anströmung besteht eine Wissenslücke für die am häufigsten auftretenden Queranströmungen. Weder der Prozess selbst, noch die Auswirkungen der Queranströmung auf den turbulenten Transport sind momentan verstanden oder quantifiziert. Das Zentralziel ist, diese Queranströmungen zu detektieren und ein konzeptionelles Verständnis durch Messungen über den landwirtschaftlichen LAFO-Flächen zu erlangen. Die Zentralhypothese ist, dass Queranströmung zu einer Zone systematisch erhöhten turbulenten Austausches innerhalb einer Übergangsgrenzschicht führt, in der die Land-Oberflächen-Rückkopplung verstärkt ist. Die Forschung wird durch 4 Unterhypothesen und -ziele geleitet. Mithilfe der hochauflösenden (Sekunden und Dezimeter) Fiber-Optic Distributed Sensing (FODS) Technik erfolgen räumlich kontinuierliche Messungen der Lufttemperatur in der Bodenschicht und im Pflanzenbestand, sowie der turbulenten Windgeschwindigkeit und Turbulenter Kinetischer Energie über Hunderte Meter hinweg. FODS misst ebenfalls (10er-Sekunden, Dezimeter) Bodentemperaturen und -feuchtigkeit nahe der Oberfläche, sowie die solare Lichtverteilung kombiniert mit Netzwerken klassischer Punktsensoren über der Winterweizen- und Maisparzelle von LAFO. Unter der Leitung von P5 baut das LAFI-Team die FODS-Glasfaserharfen auf und P1 betreibt sie. Das wissenschaftliche Programm besteht aus 4 eng verzahnten Arbeitspaketen (WP): WP1 führt eine Post-field Kalibration mithilfe unserer eigens entwickelten Kalibriersoftware durch. WP2 erforscht die in den Raum-Zeit-Daten enthaltenen Strukturen und Flüsse abhängig vom Anströmwinkel. Ergänzt durch die verbesserte Energiebilanz (P6), die fernerkundlichen Daten zur Landoberfläche (P2) und Verdunstung (P4) wird die Zentralhypothese überprüft. Die Lidarmessungen aus P1 schaffen einzigartige Beobachtungssynergien. In WP3 werden die Erkenntnisse u.a. mit der Hilfe von Deep Learning zu einem konzeptionellem Gesamtrahmen zusammengefügt. Für die zweite Phase sehen wir Erweiterungen hin zu anderen Landoberflächenheterogenitäten, der Kohlenstoffbilanz und messtechnische Weiterentwicklungen des FODS vor. Dieses Projekt befasst sich mit den zentralen LAFI-Zielen O1 bis O3, OS und OE im Falle eines klimatischen Extremums. Zusätzlich zu den oben genannten projektübergreifenden Kooperationen stehen wir im engen Austausch mit P11 wegen der Licht- und Temperaturvertikalprofile.
Die Einbeziehung von physikalischen Erhaltungsgesetzen ist bereits seit längerem als wichtiger Punkt bei der Entwicklung von numerischen Wettervorhersagemodellen bekannt. In Datenassimilation im weiteren Sinne ist die Berücksichtigung von Erhaltungsgesetzen und die Analyse deren Bedeutung erst seit kurzem ein zentraler Gegenstand der Forschung. Numerische Atmosphärenmodelle sind heute in der Lage, kleinräumige und deshalb hochgradig nichtlineare Dynamik und Physik aufzulösen. Vorhersagen hängen dabei sehr sensible von den Anfangs- und Randbedingungen ab. Die Datenassimilation für Modelle, die viele Skalen darstellen und in die sowohl zeitlich als auch räumlich hochaufgelöste Beobachtungen eingehen verlangt nach einer Überprüfung und der Weiterentwicklung der zur Zeit in numerischen Wettervorhersagemodelle genutzten, auf weniger nicht lineare Anwendungen ausgelegte Methoden. Das primäre Ziel dieses Projektes ist die Entwicklung von Ensemble-basierten Datenassimilations-Algorithmen welche Eigenschaften von nichtlinearen dynamischen Systemen, zum Beispiel die Erhaltung von Masse, Drehimpuls, Energie und Enstrophie berücksichtigen. Insbesondere werden die folgenden zwei Problemstellungen bearbeitet. Gemeinsam mit Kollegen haben wir kürzlich zeigen können, dass Erhaltung von Masse und Positivität wichtige Nebenbedingungen für Datenassimilations-Algorithmen sind. Beide Bedingungen sind in einem neuen Algorithmus, dem quadratic programming ensemble Kalman filter, umgesetzt worden. Tests an linearer Dynamik wurden ebenfalls bereits durchgeführt. Im Rahmen des vorgeschlagenen Projektes soll dieser Algorithmus in einer idealisierten Konfiguration für die Assimilation von Radarreflektivität mit dem nichthydrostatischen, konvektionsauflösenden COSMO-DE-Modell erweitert, implementiert und evaluiert werden. Zweitens, wird untersucht wie sich Datenassimilations-Algorithmen wie der Ensemble Kalman Filter und der quadratic programming ensemble Kalman filter auf die zu erhaltenden Größen im Fall von Experimenten mit einem idealisierten, nichtlinearen zweidimensionalen Flachwassermodell auswirken, sowie ob und wie diese Ensemble-basierten Algorithmen modifiziert werden können um Lösungen mit den vorgeschriebenen Eigenschaften zu erhalten. Es ist zu erwarten, dass Erhalt von Masse und Entropie die nicht lineare Energie Kaskade im System verbessert. Der mögliche Einfluss auf die Exaktheit von Vorhersagen wird auch untersucht.
| Organisation | Count |
|---|---|
| Bund | 162 |
| Europa | 3 |
| Kommune | 1 |
| Land | 4 |
| Wissenschaft | 117 |
| Type | Count |
|---|---|
| Daten und Messstellen | 9 |
| Förderprogramm | 159 |
| unbekannt | 5 |
| License | Count |
|---|---|
| Offen | 169 |
| Unbekannt | 4 |
| Language | Count |
|---|---|
| Deutsch | 110 |
| Englisch | 105 |
| Resource type | Count |
|---|---|
| Archiv | 6 |
| Datei | 3 |
| Dokument | 1 |
| Keine | 111 |
| Webseite | 53 |
| Topic | Count |
|---|---|
| Boden | 137 |
| Lebewesen und Lebensräume | 156 |
| Luft | 173 |
| Mensch und Umwelt | 173 |
| Wasser | 139 |
| Weitere | 173 |