Characterization of igneous terranes by zircon dating: implications for the UHP relicts occurrences and suture identification in the Central Rhodope, Northern Greece

One of the key for the understanding of an orogenic belt is the characterization of the terranes involved and the identification of the suture(s) separating crustal blocks: these are essential information for large-scale paleo-reconstructions. In addition, the structural relationships between the terranes involved in the collisional processes and the eventual UHP relicts may provide first order inputs to exhumation models of subducted rocks. The structure of the Rhodope Massif (northern Greece and southern Bulgaria) results from the stacking of high-grade nappes during a continental collision, which age is comprised between Latest-Jurassic and Early-Cenozoic. UHP and HP relicts, associated with oceanic and ultramafic material, suggest the presence of a dismembered suture zone within the massif. The location of this suture remains unclear; furthermore, up to now, the UHP and eclogitic localities represent isolated spots and no synthesis on their structural position within the massif has been proposed. The first aim of this work is to define the relationships between HP-UHP relicts, crustal blocks, shear zones and amphibolitic material. To achieve this objective, we characterized the accreted blocks in terms of protoliths ages of the orthogneisses mainly along two cross sections on the Greek part of the belt. Geochemical affinities of meta-igneous rocks served as a complementary tool for terrane characterization and geodynamic interpretation. Single-zircon Pb-Pb evaporation and zircon U-Pb SHRIMP dating of orthogneiss protoliths define two groups of intrusion-ages: Permo-Carboniferous and Late Jurassic-Early Cretaceous. Structurally, these two groups correspond to distinct units: the Late Jurassic gneissic complex overthrusts the one bearing the Permo-Carboniferous orthogneisses. Mylonites, eclogites, amphibolites of oceanic affinities, and UHP micaschists, mark a “melange” zone, intensively sheared towards the SW, which separates the two units. Thus, we interpret them as two distinct terranes, the Rhodope and Thracia terranes, separated by the Nestos suture. The correlation of our findings in northern Greece to the Bulgarian part of the Massif suggests a northern rooting of the Nestos Suture. This configuration results of the closure of a marginal oceanic basin of the Tethys system by a north-directed subduction. This interpretation is supported by the geochemical affinities of the orthogneisses: the Late-Jurassic igneous rocks formed by subduction-related magmatism, pprobably the same north-directed subduction that gave rise to the UHP metamorphism of the metasediments of the “melange” zone. It is noteworthy that the UHP-HP relicts seem to be restricted to the contact between the two terranes suggesting that the UHP relicts are exhumed only within the suture zone. Furthermore, the singularity of the suture suggests that the Late-Jurassic subduction explains the occurrence of UHP and eclogite relicts in the Central Rhodope despite the large age range previously attributed the UHP and/or HP stage.

The development of stylolites, from small-scale heterogeneities to multi-scale roughness

Stylolites are rough paired surfaces, indicative of localized stress-induced dissolution under a non-hydrostatic state of stress, separated by a clay parting which is believed to be the residuum of the dissolved rock. These structures are the most frequent deformation pattern in monomineralic rocks and thus provide important information about low temperature deformation and mass transfer. The intriguing roughness of stylolites can be used to assess amount of volume loss and paleo-stress directions, and to infer the destabilizing processes during pressure solution. But there is little agreement on how stylolites form and why these localized pressure solution patterns develop their characteristic roughness.rnNatural bedding parallel and vertical stylolites were studied in this work to obtain a quantitative description of the stylolite roughness and understand the governing processes during their formation. Adapting scaling approaches based on fractal principles it is demonstrated that stylolites show two self affine scaling regimes with roughness exponents of 1.1 and 0.5 for small and large length scales separated by a crossover length at the millimeter scale. Analysis of stylolites from various depths proved that this crossover length is a function of the stress field during formation, as analytically predicted. For bedding parallel stylolites the crossover length is a function of the normal stress on the interface, but vertical stylolites show a clear in-plane anisotropy of the crossover length owing to the fact that the in-plane stresses (σ2 and σ3) are dissimilar. Therefore stylolite roughness contains a signature of the stress field during formation.rnTo address the origin of stylolite roughness a combined microstructural (SEM/EBSD) and numerical approach is employed. Microstructural investigations of natural stylolites in limestones reveal that heterogeneities initially present in the host rock (clay particles, quartz grains) are responsible for the formation of the distinctive stylolite roughness. A two-dimensional numerical model, i.e. a discrete linear elastic lattice spring model, is used to investigate the roughness evolving from an initially flat fluid filled interface induced by heterogeneities in the matrix. This model generates rough interfaces with the same scaling properties as natural stylolites. Furthermore two coinciding crossover phenomena in space and in time exist that separate length and timescales for which the roughening is either balanced by surface or elastic energies. The roughness and growth exponents are independent of the size, amount and the dissolution rate of the heterogeneities. This allows to conclude that the location of asperities is determined by a polimict multi-scale quenched noise, while the roughening process is governed by inherent processes i.e. the transition from a surface to an elastic energy dominated regime.rn

Palaeo wind system reconstruction of the last glacial period over Europe, using high resolution proxy data and model-data-comparison

The atmosphere is a global influence on the movement of heat and humidity between the continents, and thus significantly affects climate variability. Information about atmospheric circulation are of major importance for the understanding of different climatic conditions. Dust deposits from maar lakes and dry maars from the Eifel Volcanic Field (Germany) are therefore used as proxy data for the reconstruction of past aeolian dynamics.rnrnIn this thesis past two sediment cores from the Eifel region are examined: the core SM3 from Lake Schalkenmehren and the core DE3 from the Dehner dry maar. Both cores contain the tephra of the Laacher See eruption, which is dated to 12,900 before present. Taken together the cores cover the last 60,000 years: SM3 the Holocene and DE3 the marine isotope stages MIS-3 and MIS-2, respectively. The frequencies of glacial dust storm events and their paleo wind direction are detected by high resolution grain size and provenance analysis of the lake sediments. Therefore two different methods are applied: geochemical measurements of the sediment using µXRF-scanning and the particle analysis method RADIUS (rapid particle analysis of digital images by ultra-high-resolution scanning of thin sections).rnIt is shown that single dust layers in the lake sediment are characterized by an increased content of aeolian transported carbonate particles. The limestone-bearing Eifel-North-South zone is the most likely source for the carbonate rich aeolian dust in the lake sediments of the Dehner dry maar. The dry maar is located on the western side of the Eifel-North-South zone. Thus, carbonate rich aeolian sediment is most likely to be transported towards the Dehner dry maar within easterly winds. A methodology is developed which limits the detection to the aeolian transported carbonate particles in the sediment, the RADIUS-carbonate module.rnrnIn summary, during the marine isotope stage MIS-3 the storm frequency and the east wind frequency are both increased in comparison to MIS-2. These results leads to the suggestion that atmospheric circulation was affected by more turbulent conditions during MIS-3 in comparison to the more stable atmospheric circulation during the full glacial conditions of MIS-2.rnThe results of the investigations of the dust records are finally evaluated in relation a study of atmospheric general circulation models for a comprehensive interpretation. Here, AGCM experiments (ECHAM3 and ECHAM4) with different prescribed SST patterns are used to develop a synoptic interpretation of long-persisting east wind conditions and of east wind storm events, which are suggested to lead to an enhanced accumulation of sediment being transported by easterly winds to the proxy site of the Dehner dry maar.rnrnThe basic observations made on the proxy record are also illustrated in the 10 m-wind vectors in the different model experiments under glacial conditions with different prescribed sea surface temperature patterns. Furthermore, the analysis of long-persisting east wind conditions in the AGCM data shows a stronger seasonality under glacial conditions: all the different experiments are characterized by an increase of the relative importance of the LEWIC during spring and summer. The different glacial experiments consistently show a shift from a long-lasting high over the Baltic Sea towards the NW, directly above the Scandinavian Ice Sheet, together with contemporary enhanced westerly circulation over the North Atlantic.rnrnThis thesis is a comprehensive analysis of atmospheric circulation patterns during the last glacial period. It has been possible to reconstruct important elements of the glacial paleo climate in Central Europe. While the proxy data from sediment cores lead to a binary signal of the wind direction changes (east versus west wind), a synoptic interpretation using atmospheric circulation models is successful. This shows a possible distribution of high and low pressure areas and thus the direction and strength of wind fields which have the capacity to transport dust. In conclusion, the combination of numerical models, to enhance understanding of processes in the climate system, with proxy data from the environmental record is the key to a comprehensive approach to paleo climatic reconstruction.rn