Ancient Greek harbours used as geo-archives for palaeotsunami research : case studies from Krane (Cefalonia), Lechaion (Gulf of Corinth) and Kyllini (Peloponnese)

Since historical times, coastal areas throughout the eastern Mediterranean are exposed to tsunami hazard. For many decades the knowledge about palaeotsunamis was solely based on historical accounts. However, results from timeline analyses reveal different characteristics affecting the quality of the dataset (i.e. distribution of data, temporal thinning backward of events, local periodization phenomena) that emphasize the fragmentary character of the historical data. As an increasing number of geo-scientific studies give convincing examples of well dated tsunami signatures not reported in catalogues, the non-existing record is a major problem to palaeotsunami research. While the compilation of historical data allows a first approach in the identification of areas vulnerable to tsunamis, it must not be regarded as reliable for hazard assessment. Considering the increasing economic significance of coastal regions (e.g. for mass tourism) and the constantly growing coastal population, our knowledge on the local, regional and supraregional tsunami hazard along Mediterranean coasts has to be improved. For setting up a reliable tsunami risk assessment and developing risk mitigation strategies, it is of major importance (i) to identify areas under risk and (ii) to estimate the intensity and frequency of potential events. This approach is most promising when based on the analysis of palaeotsunami research seeking to detect areas of high palaeotsunami hazard, to calculate recurrence intervals and to document palaeotsunami destructiveness in terms of wave run-up, inundation and long-term coastal change. Within the past few years, geo-scientific studies on palaeotsunami events provided convincing evidence that throughout the Mediterranean ancient harbours were subject to strong tsunami-related disturbance or destruction. Constructed to protect ships from storm and wave activity, harbours provide especially sheltered and quiescent environments and thus turned out to be valuable geo-archives for tsunamigenic high-energy impacts on coastal areas. Directly exposed to the Hellenic Trench and extensive local fault systems, coastal areas in the Ionian Sea and the Gulf of Corinth hold a considerably high risk for tsunami events, respectively.Geo-scientific and geoarcheaological studies carried out in the environs of the ancient harbours of Krane (Cefalonia Island), Lechaion (Corinth, Gulf of Corinth) and Kyllini (western Peloponnese) comprised on-shore and near-shore vibracoring and subsequent sedimentological, geochemical and microfossil analyses of the recovered sediments. Geophysical methods like electrical resistivity tomography and ground penetrating radar were applied in order to detect subsurface structures and to verify stratigraphical patterns derived from vibracores over long distances. The overall geochronological framework of each study area is based on radiocarbon dating of biogenic material and age determination of diagnostic ceramic fragments. Results presented within this study provide distinct evidence of multiple palaeotsunami landfalls for the investigated areas. Tsunami signatures encountered in the environs of Krane, Lechaion and Kyllini include (i) coarse-grained allochthonous marine sediments intersecting silt-dominated quiescent harbour deposits and/or shallow marine environments, (ii) disturbed microfaunal assemblages and/or (iii) distinct geochemical fingerprints as well as (iv) geo-archaeological destruction layers and (v) extensive units of beachrock-type calcarenitic tsunamites. For Krane, geochronological data yielded termini ad or post quem (maximum ages) for tsunami event generations dated to 4150 ± 60 cal BC, ~ 3200 ± 110 cal BC, ~ 650 ± 110 cal BC, and ~ 930 ± 40 cal AD, respectively. Results for Lechaion suggest that the harbour was hit by strong tsunami impacts in the 8th-6th century BC, the 1st-2nd century AD and in the 6th century AD. At Kyllini, the harbour site was affected by tsunami impact in between the late 7th and early 4th cent. BC and between the 4th and 6th cent. AD. In case of Lechaion and Kyllini, the final destruction of the harbour facilities also seems to be related to the tsunami impact. Comparing the tsunami signals obtained for each study areas with geo-scientific data from palaeotsunami events from other sites indicates that the investigated harbour sites represent excellent geo-archives for supra-regional mega-tsunamis.

Numerical investigations on atmosphere-biosphere interactions

Der Austausch von Spurengasen und Aerosolpartikeln zwischenAtmosphäre und Biosphäre spielt eine wichtige Rolle in derAtmosphärenphysik und -chemie. Wälder repräsentieren sowohleine signifikante Senke als auch Quelle für Spurengase undPartikel und tragen somit maßgeblich zu derenatmosphärischem Budget bei. Strahlungsnebel beeinflußt durchAufnahme, Entfernen und Prozessieren von Aerosolpartikelnund löslichen Spurengasen deren Konzentrationen in derGasphase. In dieser Arbeit wird erstmalig ein Modell präsentiert,welches die Simulation des Austausches zwischen Atmosphäreund Biosphäre unter Berücksichtigung der dynamischenWechselwirkung zwischen Strahlungsnebel, Blattflächenwasserund Mehrphasenchemie ermöglicht. Numerische Fallstudien mitfolgenden Schwerpunkten werden präsentiert: - Einfluß von Vegetation und Blattflächenwasser auf diezeitlichen und räumlichen Schwankungen derGrößenabhängigkeit der Flüssigphasenkonzentrationen inNebeltropfen, - Einfluß von Blattflächenwasser auf dieTrockendepositionsflüsse von Ammoniak im Wald - Simulationenwurden mit einem neuen dynamischen Depositionsmodelldurchgeführt und mit dem Widerstandsansatz verglichen -, - Einfluß von physikalischen und chemischen Prozessen aufdie Reduktion von NO- und Isoprenemissionen aus demWaldbestand verglichen mit den primären Emissionen.

Influences of deep convective cloud systems on tropospheric trace gases and photochemistry over the tropical West Pacific

A numerical model for studying the influences of deep convective cloud systems on photochemistry was developed based on a non-hydrostatic meteorological model and chemistry from a global chemistry transport model. The transport of trace gases, the scavenging of soluble trace gases, and the influences of lightning produced nitrogen oxides (NOx=NO+NO2) on the local ozone-related photochemistry were investigated in a multi-day case study for an oceanic region located in the tropical western Pacific. Model runs considering influences of large scale flows, previously neglected in multi-day cloud resolving and single column model studies of tracer transport, yielded that the influence of the mesoscale subsidence (between clouds) on trace gas transport was considerably overestimated in these studies. The simulated vertical transport and scavenging of highly soluble tracers were found to depend on the initial profiles, reconciling contrasting results from two previous studies. Influences of the modeled uptake of trace gases by hydrometeors in the liquid and the ice phase were studied in some detail for a small number of atmospheric trace gases and novel aspects concerning the role of the retention coefficient (i.e. the fraction of a dissolved trace gas that is retained in the ice phase upon freezing) on the vertical transport of highly soluble gases were illuminated. Including lightning NOx production inside a 500 km 2-D model domain was found to be important for the NOx budget and caused small to moderate changes in the domain averaged ozone concentrations. A number of sensitivity studies yielded that the fraction of lightning associated NOx which was lost through photochemical reactions in the vicinity of the lightning source was considerable, but strongly depended on assumptions about the magnitude and the altitude of the lightning NOx source. In contrast to a suggestion from an earlier study, it was argued that the near zero upper tropospheric ozone mixing ratios which were observed close to the study region were most probably not caused by the formation of NO associated with lightning. Instead, it was argued in agreement with suggestions from other studies that the deep convective transport of ozone-poor air masses from the relatively unpolluted marine boundary layer, which have most likely been advected horizontally over relatively large distances (both before and after encountering deep convection) probably played a role. In particular, it was suggested that the ozone profiles observed during CEPEX (Central Equatorial Pacific Experiment) were strongly influenced by the deep convection and the larger scale flow which are associated with the intra-seasonal oscillation.

Exploring physical processes related to past climate proxies: lake sediments and stable water isotopes

Proxy data are essential for the investigation of climate variability on time scales larger than the historical meteorological observation period. The potential value of a proxy depends on our ability to understand and quantify the physical processes that relate the corresponding climate parameter and the signal in the proxy archive. These processes can be explored under present-day conditions. In this thesis, both statistical and physical models are applied for their analysis, focusing on two specific types of proxies, lake sediment data and stable water isotopes.rnIn the first part of this work, the basis is established for statistically calibrating new proxies from lake sediments in western Germany. A comprehensive meteorological and hydrological data set is compiled and statistically analyzed. In this way, meteorological times series are identified that can be applied for the calibration of various climate proxies. A particular focus is laid on the investigation of extreme weather events, which have rarely been the objective of paleoclimate reconstructions so far. Subsequently, a concrete example of a proxy calibration is presented. Maxima in the quartz grain concentration from a lake sediment core are compared to recent windstorms. The latter are identified from the meteorological data with the help of a newly developed windstorm index, combining local measurements and reanalysis data. The statistical significance of the correlation between extreme windstorms and signals in the sediment is verified with the help of a Monte Carlo method. This correlation is fundamental for employing lake sediment data as a new proxy to reconstruct windstorm records of the geological past.rnThe second part of this thesis deals with the analysis and simulation of stable water isotopes in atmospheric vapor on daily time scales. In this way, a better understanding of the physical processes determining these isotope ratios can be obtained, which is an important prerequisite for the interpretation of isotope data from ice cores and the reconstruction of past temperature. In particular, the focus here is on the deuterium excess and its relation to the environmental conditions during evaporation of water from the ocean. As a basis for the diagnostic analysis and for evaluating the simulations, isotope measurements from Rehovot (Israel) are used, provided by the Weizmann Institute of Science. First, a Lagrangian moisture source diagnostic is employed in order to establish quantitative linkages between the measurements and the evaporation conditions of the vapor (and thus to calibrate the isotope signal). A strong negative correlation between relative humidity in the source regions and measured deuterium excess is found. On the contrary, sea surface temperature in the evaporation regions does not correlate well with deuterium excess. Although requiring confirmation by isotope data from different regions and longer time scales, this weak correlation might be of major importance for the reconstruction of moisture source temperatures from ice core data. Second, the Lagrangian source diagnostic is combined with a Craig-Gordon fractionation parameterization for the identified evaporation events in order to simulate the isotope ratios at Rehovot. In this way, the Craig-Gordon model can be directly evaluated with atmospheric isotope data, and better constraints for uncertain model parameters can be obtained. A comparison of the simulated deuterium excess with the measurements reveals that a much better agreement can be achieved using a wind speed independent formulation of the non-equilibrium fractionation factor instead of the classical parameterization introduced by Merlivat and Jouzel, which is widely applied in isotope GCMs. Finally, the first steps of the implementation of water isotope physics in the limited-area COSMO model are described, and an approach is outlined that allows to compare simulated isotope ratios to measurements in an event-based manner by using a water tagging technique. The good agreement between model results from several case studies and measurements at Rehovot demonstrates the applicability of the approach. Because the model can be run with high, potentially cloud-resolving spatial resolution, and because it contains sophisticated parameterizations of many atmospheric processes, a complete implementation of isotope physics will allow detailed, process-oriented studies of the complex variability of stable isotopes in atmospheric waters in future research.rn

Potential vorticity and moisture in extratropical cyclones : climatology and sensitivity experiments

The development of extratropical cyclones can be seen as an interplay of three positive potential vorticity (PV) anomalies: an upper-level stratospheric intrusion, low-tropospheric diabatically produced PV, and a warm anomaly at the surface acting as a surrogate PV anomaly. In the mature stage they become vertically aligned and form a “PV tower” associated with strong cyclonic circulation. This paradigm of extratropical cyclone development provides the basis of this thesis, which will use a climatological dataset and numerical model experiments to investigate the amplitude of the three anomalies and the processes leading in particular to the formation of the diabatically produced low-tropospheric PV anomaly.rnrnThe first part of this study, based on the interim ECMWF Re-Analysis (ERA-Interim) dataset, quantifies the amplitude of the three PV anomalies in mature extratropical cyclones in different regions in the Northern Hemisphere on a climatological basis. A tracking algorithm is applied to sea level pressure (SLP) fields to identify cyclone tracks. Surface potential temperature anomalies ∆θ and vertical profiles of PV anomalies ∆PV are calculated at the time of the cyclones’ minimum SLP and during the intensification phase 24 hours before in a vertical cylinder with a radius of 200 km around the surface cyclone center. To compare the characteristics of the cyclones, they are grouped according to their location (8 regions) and intensity, where the central SLP is used as a measure of intensity. Composites of ∆PV profiles and ∆θ are calculated for each region and intensity class at the time of minimum SLP and during the cyclone intensification phase.rnrnDuring the cyclones’ development stage the amplitudes of all three anomalies increase on average. In the mature stage all three anomalies are typically larger for intense than for weak winter cyclones [e.g., 0.6 versus 0.2 potential vorticity units (PVU) at lower levels, and 1.5 versus 0.5 PVU at upper levels].rnThe regional variability of the cyclones’ vertical structure and the profile evolution is prominent (cyclones in some regions are more sensitive to the amplitude of a particular anomaly than in other regions). Values of ∆θ and low-level ∆PV are on average larger in the western parts of the oceans than in the eastern parts. In addition, a large seasonal variability can be identified, with fewer and weaker cyclones especially in the summer, associated with higher low-tropospheric PV values, but also with a higher tropopause and much weaker surface potential temperature anomalies (compared to winter cyclones).rnrnIn the second part, we were interested in the diabatic low-level part of PV towers. Evaporative sources were identified of moisture that was involved in PV production through condensation. Lagrangian backward trajectories were calculated from the region with high PV values at low-levels in the cyclones. PV production regions were identified along these trajectories and from these regions a new set of backward trajectories was calculated and moisture uptakes were traced along them. The main contribution from surface evaporation to the specific humidity of the trajectories is collected 12-72 hours prior to therntime of PV production. The uptake region for weaker cyclones with less PV in the centre is typically more localized with reduced uptake values compared to intense cyclones. However, in a qualitative sense uptakes and other variables along single trajectories do not vary much between cyclones of different intensity in different regions.rnrnA sensitivity study with the COSMO model comprises the last part of this work. The study aims at investigating the influence of synthetic moisture modification in the cyclone environment in different stages of its development. Moisture was eliminated in three regions, which were identified as important moisture source regions for PV production. Moisture suppression affected the cyclone the most in its early phase. It led to cyclolysis shortly after its genesis. Nevertheles, a new cyclone formed on the other side of a dry box and developed relatively quickly. Also in other experiments, moisture elimination led to strong intensity reduction of the surface cyclone, limited upper-level development, and delayed or missing interaction between the two.rnrnIn summary, this thesis provides novel insight into the structure of different intensity categories of extratropical cyclones from a PV perspective, which corroborates the findings from a series of previous case studies. It reveals that all three PV anomalies are typically enhanced for more intense cyclones, with important regional differences concerning the relative amplitude of the three anomalies. The moisture source analysis is the first of this kind to study the evaporation-condensation cycle related to the intensification of extratropical cyclones. Interestingly, most of the evaporation occurs during the 3 days prior to the time of maximum cyclone intensity and typically extends over fairly large areas along the track of the cyclone. The numerical model case study complements this analysis by analyzing the impact of regionally confined moisture sources for the evolution of the cyclone.

Massenspektrometrische Wasserdampfmessung in der oberen Troposphäre und unteren Stratosphäre

Diese Arbeit beschreibt die Entwicklung des flugzeuggetragenen Atmosphärischen Ionisations-Massenspektrometers AIMS-H2O zur Messung von Wasserdampf in der oberen Troposphäre und unteren Stratosphäre (UTLS) und erste Flugzeugmessungen mit dem Instrument. Wasserdampf beeinflusst das Klima in der UTLS aufgrund seiner Strahlungseigenschaften und agiert als wichtiger Parameter bei der Bildung von Zirruswolken und Kondensstreifen. Deshalb sind genaue Wasserdampfmessungen für das Verständnis vieler atmosphärischer Prozesse unerlässlich. Instrumentenvergleiche wie sie im SPARC Report No. 2 und dem Bericht der AUQAVIT Kampagne zusammengefasst sind, haben gezeigt, dass große Abweichungen zwischen einzelnen Methoden und Instrumenten bestehen. Diese Unsicherheiten limitieren das Verständnis des Einflusses von Wasserdampf auf die Dynamik und die Strahlungseigenschaften in der UTLS. Die in dieser Arbeit vorgestellte Entwicklung einer neuen Messmethode für Wasserdampf mit dem Massenspektrometer AIMS-H2O ist deshalb auf die genaue Messung niedriger Wasserdampfkonzentrationen in der UTLS fokussiert. Mit AIMS H2O wird Umgebungsluft in einer neu entwickelten Gasentladungsquelle ionisiert. Durch eine Reihe von Ionen-Molekül-Reaktionen entstehen H3O+(H2O) und H3O+(H2O)2 Ionen. Diese Ionen werden genutzt, um die Wasserdampfkonzentration in der Atmosphäre zu bestimmen. Um die erforderliche hohe Genauigkeit zu erzielen, wird AIMS H2O im Flug kalibriert. In dem zu diesem Zweck aufgebauten Kalibrationsmodul wird die katalytische Reaktion von Wasserstoff und Sauerstoff auf einer Platinoberfläche genutzt, um definierte Wasserdampfkonzentrationen für die Kalibration im Flug zu erzeugen. Bei ersten Messungen auf der Falcon während der Kampagne CONCERT 2011 konnte dabei eine Genauigkeit von 8 bis 15% für die Messung der Wasserdampfkonzentration in einem Messbereich von 0,5 bis 250 ppmv erreicht werden. Die Messfrequenz betrug 4 Hz, was einer räumlichen Auflösung von etwa 50 m entspricht. Der Vergleich der Messung des Massenspektrometers mit dem Laserhygrometer Waran zeigt eine sehr gute Übereinstimmung im Rahmen der Unsicherheiten. Anhand zweier Fallstudien werden die Messungen von AIMS H2O während CONCERT 2011 detailliert analysiert. In der ersten Studie werden zwei Flüge in eine stratosphärische Intrusion über Nordeuropa untersucht. In dieser Situation wurde stratosphärische Luft bis hinunter auf 6 km Höhe transportiert und war dadurch mit der Falcon erreichbar. Es konnte gezeigt werden, dass AIMS-H2O sehr gut für die genaue Messung niedriger Wasserdampfkonzentrationen, in diesem Fall bis etwa 3,5 ppmv, geeignet ist. Der Vergleich der Messung mit Analysen des ECMWF Integrated Forecast Systems zeigt eine gute Übereinstimmung der gemessenen Wasserdampfstrukturen mit der dynamischen Tropopause. Unterschiede tauchen dagegen beim Vergleich der Wasserdampfkonzentrationen in der unteren Stratosphäre auf. Hier prognostiziert das Modell deutlich höhere Feuchten. Die zweite Fallstudie beschäftigt sich mit der Verteilung der relativen Feuchte in jungen Kondensstreifen im Vergleich zu ihrer direkten Umgebung. Dabei wurde für drei Messsequenzen im Abgasstrahl von Flugzeugen beobachtet, dass die relative Feuchte innerhalb des Kondensstreifens im Vergleich zur Umgebung sowohl bei unter- als auch übersättigten Umgebungsbedingungen in Richtung Sättigung verschoben ist. Die hohe Anzahl an Eispartikeln und die damit verbundene große Eisoberfläche in jungen Kondensstreifen führt also zu einer schnellen Relaxation von Gasphase und Eis in Richtung Gleichgewicht. In der Zukunft soll AIMS-H2O auch auf HALO für die genaue Messung von Wasserdampf bei ML-CIRRUS und weiteren Kampagnen eingesetzt werden.