Aeolian geomorphodynamics in endorheic basins of the Mongolian Gobi Desert

Die vorliegende Arbeit ist im Zuge des DFG Projektes Spätpleistozäne, holozäne und aktuelle Geomorphodynamik in abflusslosen Becken der Mongolischen Gobi´´ entstanden. Das Arbeitsgebiet befindet sich in der südlichen Mongolei im nördlichen Teil der Wüste Gobi. Neben einigen Teilen der Sahara (Heintzenberg, 2009), beispielsweise das Bodélé Becken des nördlichen Tschads (z.B. Washington et al., 2006a; Todd et al., 2006; Warren et al., 2007) wird Zentralasien als ein Hauptliefergebiet für Partikel in die globale Zirkulation der Atmosphäre gesehen (Goudie, 2009). Hauptaugenmerk liegt hierbei besonders auf den abflusslosen Becken und deren Sedimentablagerungen. Die, der Deflation ausgesetzten Flächen der Seebecken, sind hauptsächliche Quelle für Partikel die sich in Form von Staub respektive Sand ausbreiten. Im Hinblick auf geomorphologische Landschaftsentwicklung wurde der Zusammenhang von Beckensedimenten zu Hangdepositionen numerisch simuliert. Ein von Grunert and Lehmkuhl (2004) publiziertes Model, angelehnt an Ideen von Pye (1995) wird damit in Betracht gezogen. Die vorliegenden Untersuchungen modellieren Verbreitungsmechanismen auf regionaler Ebene ausgehend von einer größeren Anzahl an einzelnen punktuellen Standorten. Diese sind repräsentativ für die einzelnen geomorphologischen Systemglieder mit möglicherweise einer Beteiligung am Budget aeolischer Geomorphodynamik. Die Bodenbedeckung durch das charakteristische Steinpflaster der Gobi - Region, sowie unter anderem Korngrößenverteilungen der Oberflächensedimente wurden untersucht. Des Weiteren diente eine zehnjährige Zeitreihe (Jan 1998 bis Dez 2007) meteorologischer Daten als Grundlage zur Analyse der Bedingungen für äolische Geomorphodynamik. Die Daten stammen von 32 staatlichen mongolischen Wetterstationen aus der Region und Teile davon wurden für die Simulationen verwendet. Zusätzlich wurden atmosphärische Messungen zur Untersuchung der atmosphärischen Stabilität und ihrer tageszeitlichen Variabilität mit Mess-Drachenaufstiegen vorgenommen. Die Feldbefunde und auch die Ergebnisse der Laboruntersuchungen sowie der Datensatz meteorologischer Parameter dienten als Eingangsparameter für die Modellierungen. Emissionsraten der einzelnen Standorte und die Partikelverteilung im 3D Windfeld wurden modelliert um die Konvektivität der Beckensedimente und Hangdepositionen zu simulieren. Im Falle hoher mechanischer Turbulenz der bodennahen Luftschicht (mit einhergehender hoher Wind Reibungsgeschwindigkeit), wurde generell eine neutrale Stabilität festgestellt und die Simulationen von Partikelemission sowie deren Ausbreitung und Deposition unter neutraler Stabilitätsbedingung berechnet. Die Berechnung der Partikelemission wurde auf der Grundlage eines sehr vereinfachten missionsmodells in Anlehnung an bestehende Untersuchungen (Laurent et al., 2006; Darmenova et al., 2009; Shao and Dong, 2006; Alfaro, 2008) durchgeführt. Sowohl 3D Windfeldkalkulationen als auch unterschiedliche Ausbreitungsszenarien äolischer Sedimente wurden mit dem kommerziellen Programm LASAT® (Lagrange-Simulation von Aerosol-Transport) realisiert. Diesem liegt ein Langargischer Algorithmus zugrunde, mittels dessen die Verbreitung einzelner Partikel im Windfeld mit statistischer Wahrscheinlichkeit berechnet wird. Über Sedimentationsparameter kann damit ein Ausbreitungsmodell der Beckensedimente in Hinblick auf die Gebirgsfußflächen und -hänge generiert werden. Ein weiterer Teil der Untersuchungen beschäftigt sich mit der geochemischen Zusammensetzung der Oberflächensedimente. Diese Proxy sollte dazu dienen die simulierten Ausbreitungsrichtungen der Partikel aus unterschiedlichen Quellregionen nach zu verfolgen. Im Falle der Mongolischen Gobi zeigte sich eine weitestgehende Homogenität der Minerale und chemischen Elemente in den Sedimenten. Laser Bebohrungen einzelner Sandkörner zeigten nur sehr leichte Unterschiede in Abhängigkeit der Quellregionen. Die Spektren der Minerale und untersuchten Elemente deuten auf graitische Zusammensetzungen hin. Die, im Untersuchungsgebiet weit verbreiteten Alkali-Granite (Jahn et al., 2009) zeigten sich als hauptverantwortlich für die Sedimentproduktion im Untersuchungsgebiet. Neben diesen Mineral- und Elementbestimmungen wurde die Leichtmineralfraktion auf die Charakteristik des Quarzes hin untersucht. Dazu wurden Quarzgehalt, Kristallisation und das Elektronen-Spin-Resonanz Signal des E’1 - Centers in Sauerstoff Fehlstellungen des SiO2 Gitters bestimmt. Die Untersuchungen sind mit dem Methodenvorschlag von Sun et al. (2007) durchgeführt worden und sind prinzipiell gut geeignet um Herkunftsanalysenrndurchzuführen. Eine signifikante Zuordnung der einzelnen Quellgebiete ist jedoch auch in dieser Proxy nicht zu finden gewesen.

Methoden zur Analyse, Interpretation und Modellierung des elektrischen Widerstandes von archäologischen Steinstrukturen

In der Archäologie werden elektrische Widerstandsmessungen routinemäßig zur Prospektion von Fundstellen eingesetzt. Die Methode ist kostengünstig, leicht anwendbar und liefert in den meisten Fällen zuverlässige und leicht zu interpretierende Ergebnisse. Dennoch kann die Methode die archäologischen Strukturen in manchen Fällen nur teilweise oder gar nicht abbilden, wenn die bodenphysikalischen und bodenchemischen Eigenschaften des Bodens und der archäologischen Strukturen dies nicht zulassen. Der spezifische elektrische Widerstand wird durch Parameter wie Wassergehalt, Bodenstruktur, Bodenskelett, Bodentextur, Salinität und Bodentemperatur beeinflusst. Manche dieser Parameter, wie z.B. der Wassergehalt und die Bodentemperatur, unterliegen einer saisonalen Veränderung. Die vorliegende Arbeit untersucht den spezifischen elektrischen Widerstand von archäologischen Steinstrukturen und evaluiert die Möglichkeit, auf Grundlage von Geländemessungen und Laboranalysen archäologische Strukturen und Böden als numerische Modelle darzustellen. Dazu wurde eine Kombination von verschiedenen bodenkundlichen, geoarchäologischen und geophysikalischen Methoden verwendet. Um archäologische Strukturen und Bodenprofile als numerische Widerstandsmodelle darstellen zu können, werden Informationen zur Geometrie der Strukturen und ihren elektrischen Widerstandswerten benötigt. Dabei ist die Qualität der Hintergrundinformationen entscheidend für die Genauigkeit des Widerstandsmodells. Die Geometrie der Widerstandsmodelle basiert auf den Ergebnissen von Rammkernsondierungen und archäologische Ausgrabungen. Die an der Ausbildung des elektrischen Widerstands beteiligten Parameter wurden durch die Analyse von Bodenproben gemessen und ermöglichen durch Pedotransfer-Funktion, wie die Rhoades-Formel, die Abschätzung des spezifischen elektrischen Widerstandes des Feinbodens. Um den Einfluss des Bodenskeletts auf den spezifischen elektrischen Widerstand von Bodenprofilen und archäologischen Strukturen zu berechnen, kamen die Perkolationstheorie und die Effective Medium Theory zum Einsatz. Die Genauigkeit und eventuelle Limitierungen der Methoden wurden im Labor durch experimentelle Widerstandsmessungen an ungestörten Bodenproben und synthetischen Materialien überprüft. Die saisonale Veränderung des Wassergehalts im Boden wurde durch numerische Modelle mit der Software HYDRUS simuliert. Die hydraulischen Modelle wurden auf Grundlage der ermittelten bodenkundlichen und archäologischen Stratigraphie erstellt und verwenden die Daten von lokalen Wetterstationen als Eingangsparameter. Durch die Kombination der HYDRUS-Ergebnisse mit den Pedotransfer-Funktionen konnte der Einfluss dieser saisonalen Veränderung auf die Prospektionsergebnisse von elektrischen Widerstandsmethoden berechnet werden. Die Ergebnisse der Modellierungsprozesse wurden mit den Geländemessungen verglichen. Die beste Übereinstimmung zwischen Modellergebnissen und den Prospektionsergebnissen konnte für die Fallstudie bei Katzenbach festgestellt werden. Bei dieser wurden die Modelle auf Grundlage von archäologischen Grabungsergebnissen und detaillierten bodenkundlichen Analysen erstellt. Weitere Fallstudien zeigen, dass elektrische Widerstandsmodelle eingesetzt werden können, um den Einfluss von ungünstigen Prospektionsbedingungen auf die Ergebnisse der elektrischen Widerstandsmessungen abzuschätzen. Diese Informationen unterstützen die Planung und Anwendung der Methoden im Gelände und ermöglichen eine effektivere Interpretation der Prospektionsergebnisse. Die präsentierten Modellierungsansätze benötigen eine weitere Verifizierung durch den Vergleich der Modellierungsergebnisse mit detailliertem geophysikalischem Gelände-Monitoring von archäologischen Fundstellen. Zusätzlich könnten elektrische Widerstandsmessungen an künstlichen Mauerstrukturen unter kontrollierten Bedingungen zur Überprüfung der Modellierungsprozesse genutzt werden.

Die Rolle der Globalstrahlung im Klimasystem Südwestdeutschlands

Die Untersuchungen umfassen die Periode 1981 – 2000 und basieren hauptsächlich auf Daten des Deutschen Wetterdienstes (DWD). Relativwerte der Globalstrahlung beziehen sich auf die Rayleigh-Atmosphäre. Das Regressionsmodell nach Angström ermöglicht die Erweiterung des Meßnetzes. In linearer und nichtlinearer Regression und Korrelation ist die Globalstrahlung entweder abhängige (Sonnenscheindauer, Bewölkung) oder unabhängige Variable (Lufttemperatur, Bodentemperatur). Ihre Intensität in Abhängigkeit von Großwetterlagen, Großwettertypen und Luftmassen wird diskutiert. Diesbezüglich werden mit der Linearen Diskriminanzanalyse ähnliche Großwetterlagen und Stationen in signifikant unterschiedenen Gruppen zusammengefaßt, getrennt nach Sommer- und Winterhalbjahr. Abhängig von der Zeit betrachtet, enthalten Globalstrahlung, direkte und diffuse Sonnenstrahlung, Lufttemperatur, Bewölkung und Niederschlag signifikante zyklische Variationen, die gegebenenfalls klimatologisch relevant sind. Weiteren Aufschluß ergeben deshalb die Zeitreihenanalysen. Autokorrelation-Spektralanalysen (ASA) der genannten Variablen werden in integrierten Spektren dargestellt. Hinweise auf die zeitliche Konstanz signifikanter Varianzmaxima enthalten die Spektren der dynamischen (gleitenden) ASA.

Landwirtschaftliche Schaderregerprognose mit Hilfe von geographischen Informationssystemen

In der hier vorliegenden Arbeit wurde am Beispiel der Kraut- und Knollenfäule an Kartoffeln Phytophthora infestans und des Kartoffelkäfers Leptinotarsa decemlineata untersucht, ob durch den Einsatz von Geographischen Informationssystemen (GIS) landwirtschaftliche Schader¬reger¬prognosen für jeden beliebigen Kartoffelschlag in Deutschland erstellt werden können. Um dieses Ziel zu erreichen, wurden die Eingangsparameter (Temperatur und relative Luftfeuchte) der Prognosemodelle für die beiden Schaderreger (SIMLEP1, SIMPHYT1, SIMPHYT3 and SIMBLIGHT1) so aufbereitet, dass Wetterdaten flächendeckend für Deutschland zur Verfügung standen. Bevor jedoch interpoliert werden konnte, wurde eine Regionalisierung von Deutschland in Interpolationszonen durchgeführt und somit Naturräume geschaffen, die einen Vergleich und eine Bewertung der in ihnen liegenden Wetterstationen zulassen. Hierzu wurden die Boden-Klima-Regionen von SCHULZKE und KAULE (2000) modifiziert, an das Wetterstationsnetz angepasst und mit 5 bis 10 km breiten Pufferzonen an der Grenze der Interpolationszonen versehen, um die Wetterstationen so häufig wie möglich verwenden zu können. Für die Interpolation der Wetterdaten wurde das Verfahren der multiplen Regression gewählt, weil dieses im Vergleich zu anderen Verfahren die geringsten Abweichungen zwischen interpolierten und gemessenen Daten aufwies und den technischen Anforderungen am besten entsprach. Für 99 % aller Werte konnten bei der Temperaturberechnung Abweichungen in einem Bereich zwischen -2,5 und 2,5 °C erzielt werden. Bei der Berechnung der relativen Luftfeuchte wurden Abweichungen zwischen -12 und 10 % relativer Luftfeuchte erreicht. Die Mittelwerte der Abweichungen lagen bei der Temperatur bei 0,1 °C und bei der relativen Luftfeuchte bei -1,8 %. Zur Überprüfung der Trefferquoten der Modelle beim Betrieb mit interpolierten Wetterdaten wurden Felderhebungsdaten aus den Jahren 2000 bis 2007 zum Erstauftreten der Kraut- und Knollenfäule sowie des Kartoffelkäfers verwendet. Dabei konnten mit interpolierten Wetterdaten die gleichen und auch höhere Trefferquoten erreicht werden, als mit der bisherigen Berechnungsmethode. Beispielsweise erzielte die Berechnung des Erstauftretens von P. infestans durch das Modell SIMBLIGHT1 mit interpolierten Wetterdaten im Schnitt drei Tage geringere Abweichungen im Vergleich zu den Berechnungen ohne GIS. Um die Auswirkungen interpretieren zu können, die durch Abweichungen der Temperatur und der relativen Luftfeuchte entstanden wurde zusätzlich eine Sensitivitätsanalyse zur Temperatur und relativen Luftfeuchte der verwendeten Prognosemodelle durchgeführt. Die Temperatur hatte bei allen Modellen nur einen geringen Einfluss auf das Prognoseergebnis. Veränderungen der relativen Luftfeuchte haben sich dagegen deutlich stärker ausgewirkt. So lag bei SIMBLIGHT1 die Abweichung durch eine stündliche Veränderung der relativen Luftfeuchte (± 6 %) bei maximal 27 Tagen, wogegen stündliche Veränderungen der Temperatur (± 2 °C) eine Abweichung von maximal 10 Tagen ausmachten. Die Ergebnisse dieser Arbeit zeigen, dass durch die Verwendung von GIS mindestens die gleichen und auch höhere Trefferquoten bei Schaderregerprognosen erzielt werden als mit der bisherigen Verwendung von Daten einer nahegelegenen Wetterstation. Die Ergebnisse stellen einen wesentlichen Fortschritt für die landwirtschaftlichen Schaderregerprognosen dar. Erstmals ist es möglich, bundesweite Prognosen für jeden beliebigen Kartoffelschlag zur Bekämpfung von Schädlingen in der Landwirtschaft bereit zu stellen.