An observational study of oceanic eddy generation mechanisms by tall deep-water islands (Gran Canaria)

[EN]Oceanic eddy generation by tall deep-water islands is common phenomenon. It is recognized that these eddies may have a significant impact on the marine system and related biogeochemical fluxes. Hence, it is important to establish favourable conditions for their generation. With this objective, we present an observational study on eddy generation mechanisms by tall deep-water islands, using as a case study the island of Gran Canaria. Observations show that the main generation mechanism is topographic forcing, which leads to eddy generation when the incident oceanic flow is sufficiently intense. Wind shear at the island wake may acts only as an additional eddy-generation trigger mechanism when the impinging oceanic flow is not sufficiently intense. For the case of the island of Gran Canaria we have observed a mean of ten generated cyclonic eddies per year. Eddies are more frequently generated in summer coinciding with intense Trade winds and Canary Current.

A dynamic model of oceanic sulfur (DMOS) applied to the Sargasso Sea: Simulating the dimethylsulfide (DMS) summer paradox

[EN]A new one-dimensional model of DMSP/DMS dynamics (DMOS) is developed and applied to the Sargasso Sea in order to explain what drives the observed dimethylsulfide (DMS) summer paradox: a summer DMS concentration maximum concurrent with a minimum in the biomass of phytoplankton, the producers of the DMS precursor dimethylsulfoniopropionate (DMSP). Several mechanisms have been postulated to explain this mismatch: a succession in phytoplankton species composition towards higher relative abundances of DMSP producers in summer; inhibition of bacterial DMS consumption by ultraviolet radiation (UVR); and direct DMS production by phytoplankton due to UVR-induced oxidative stress. None of these hypothetical mechanisms, except for the first one, has been tested with a dynamic model. We have coupled a new sulfur cycle model that incorporates the latest knowledge on DMSP/DMS dynamics to a preexisting nitrogen/carbon-based ecological model that explicitly simulates the microbial-loop. This allows the role of bacteria in DMS production and consumption to be represented and quantified. The main improvements of DMOS with respect to previous DMSP/DMS models are the explicit inclusion of: solar-radiation inhibition of bacterial sulfur uptakes; DMS exudation by phytoplankton caused by solar-radiation-induced stress; and uptake of dissolved DMSP by phytoplankton. We have conducted a series of modeling experiments where some of the DMOS sulfur paths are turned “off” or “on,” and the results on chlorophyll-a, bacteria, DMS, and DMSP (particulate and dissolved) concentrations have been compared with climatological data of these same variables. The simulated rate of sulfur cycling processes are also compared with the scarce data available from previous works. All processes seem to play a role in driving DMS seasonality. Among them, however, solar-radiation-induced DMS exudation by phytoplankton stands out as the process without which the model is unable to produce realistic DMS simulations and reproduce the DMS summer paradox.

Velocity structure and seismic anisotropy in the crust and upper mantle from Receiver Function analysis: three case studies in Italy

The research for this PhD project consisted in the application of the RFs analysis technique to different data-sets of teleseismic events recorded at temporary and permanent stations located in three distinct study regions: Colli Albani area, Northern Apennines and Southern Apennines. We found some velocity models to interpret the structures in these regions, which possess very different geologic and tectonics characteristics and therefore offer interesting case study to face. In the Colli Albani some of the features evidenced in the RFs are shared by all the analyzed stations: the Moho is almost flat and is located at about 23 km depth, and the presence of a relatively shallow limestone layer is a stable feature; contrariwise there are features which vary from station to station, indicating local complexities. Three seismic stations, close to the central part of the former volcanic edifice, display relevant anisotropic signatures­­­ with symmetry axes consistent with the emplacement of the magmatic chamber. Two further anisotropic layers are present at greater depth, in the lower crust and the upper mantle, respectively, with symmetry axes directions related to the evolution of the volcano complex. In Northern Apennines we defined the isotropic structure of the area, finding the depth of the Tyrrhenian (almost 25 km and flat) and Adriatic (40 km and dipping underneath the Apennines crests) Mohos. We determined a zone in which the two Mohos overlap, and identified an anisotropic body in between, involved in the subduction and going down with the Adiratic Moho. We interpreted the downgoing anisotropic layer as generated by post-subduction delamination of the top-slab layer, probably made of metamorphosed crustal rocks caught in the subduction channel and buoyantly rising toward the surface. In the Southern Apennines, we found the Moho depth for 16 seismic stations, and highlighted the presence of an anisotropic layer underneath each station, at about 15-20 km below the whole study area. The moho displays a dome-like geometry, as it is shallow (29 km) in the central part of the study area, whereas it deepens peripherally (down to 45 km); the symmetry axes of anisotropic layer, interpreted as a layer separating the upper and the lower crust, show a moho-related pattern, indicated by the foliation of the layer which is parallel to the Moho trend. Moreover, due to the exceptional seismic event occurred on April 6th next to L’Aquila town, we determined the Vs model for two station located next to the epicenter. An extremely high velocity body is found underneath AQU station at 4-10 km depth, reaching Vs of about 4 km/s, while this body is lacking underneath FAGN station. We compared the presence of this body with other recent works and found an anti-correlation between the high Vs body, the max slip patches and earthquakes distribution. The nature of this body is speculative since such high velocities are consistent with deep crust or upper mantle, but can be interpreted as a as high strength barrier of which the high Vs is a typical connotation.

Constraints on mantle melting from major and trace element systematics in residual abyssal peridotites

ZusammenfassungDie Bildung von mittelozeanischen Rückenbasalten (MORB) ist einer der wichtigsten Stoffflüsse der Erde. Jährlich wird entlang der 75.000 km langen mittelozeanischen Rücken mehr als 20 km3 neue magmatische Kruste gebildet, das sind etwa 90 Prozent der globalen Magmenproduktion. Obwohl ozeanische Rücken und MORB zu den am meisten untersuchten geologischen Themenbereichen gehören, existieren weiterhin einige Streit-fragen. Zu den wichtigsten zählt die Rolle von geodynamischen Rahmenbedingungen, wie etwa Divergenzrate oder die Nähe zu Hotspots oder Transformstörungen, sowie der absolute Aufschmelzgrad, oder die Tiefe, in der die Aufschmelzung unter den Rücken beginnt. Diese Dissertation widmet sich diesen Themen auf der Basis von Haupt- und Spurenelementzusammensetzungen in Mineralen ozeanischer Mantelgesteine.Geochemische Charakteristika von MORB deuten darauf hin, dass der ozeanische Mantel im Stabilitätsfeld von Granatperidotit zu schmelzen beginnt. Neuere Experimente zeigen jedoch, dass die schweren Seltenerdelemente (SEE) kompatibel im Klinopyroxen (Cpx) sind. Aufgrund dieser granatähnlichen Eigenschaft von Cpx wird Granat nicht mehr zur Erklärung der MORB Daten benötigt, wodurch sich der Beginn der Aufschmelzung zu geringeren Drucken verschiebt. Aus diesem Grund ist es wichtig zu überprüfen, ob diese Hypothese mit Daten von abyssalen Peridotiten in Einklang zu bringen ist. Diese am Ozeanboden aufgeschlossenen Mantelfragmente stellen die Residuen des Aufschmelz-prozesses dar, und ihr Mineralchemismus enthält Information über die Bildungs-bedingungen der Magmen. Haupt- und Spurenelementzusammensetzungen von Peridotit-proben des Zentralindischen Rückens (CIR) wurden mit Mikrosonde und Ionensonde bestimmt, und mit veröffentlichten Daten verglichen. Cpx der CIR Peridotite weisen niedrige Verhältnisse von mittleren zu schweren SEE und hohe absolute Konzentrationen der schweren SEE auf. Aufschmelzmodelle eines Spinellperidotits unter Anwendung von üblichen, inkompatiblen Verteilungskoeffizienten (Kd's) können die gemessenen Fraktionierungen von mittleren zu schweren SEE nicht reproduzieren. Die Anwendung der neuen Kd's, die kompatibles Verhalten der schweren SEE im Cpx vorhersagen, ergibt zwar bessere Resultate, kann jedoch nicht die am stärksten fraktionierten Proben erklären. Darüber hinaus werden sehr hohe Aufschmelzgrade benötigt, was nicht mit Hauptelementdaten in Einklang zu bringen ist. Niedrige (~3-5%) Aufschmelzgrade im Stabilitätsfeld von Granatperidotit, gefolgt von weiterer Aufschmelzung von Spinellperidotit kann jedoch die Beobachtungen weitgehend erklären. Aus diesem Grund muss Granat weiterhin als wichtige Phase bei der Genese von MORB betrachtet werden (Kapitel 1).Eine weitere Hürde zum quantitativen Verständnis von Aufschmelzprozessen unter mittelozeanischen Rücken ist die fehlende Korrelation zwischen Haupt- und Spuren-elementen in residuellen abyssalen Peridotiten. Das Cr/(Cr+Al) Verhältnis (Cr#) in Spinell wird im Allgemeinen als guter qualitativer Indikator für den Aufschmelzgrad betrachtet. Die Mineralchemie der CIR Peridotite und publizierte Daten von anderen abyssalen Peridotiten zeigen, dass die schweren SEE sehr gut (r2 ~ 0.9) mit Cr# der koexistierenden Spinelle korreliert. Die Auswertung dieser Korrelation ergibt einen quantitativen Aufschmelz-indikator für Residuen, welcher auf dem Spinellchemismus basiert. Damit kann der Schmelzgrad als Funktion von Cr# in Spinell ausgedrückt werden: F = 0.10×ln(Cr#) + 0.24 (Hellebrand et al., Nature, in review; Kapitel 2). Die Anwendung dieses Indikators auf Mantelproben, für die keine Ionensondendaten verfügbar sind, ermöglicht es, geochemische und geophysikalischen Daten zu verbinden. Aus geodynamischer Perspektive ist der Gakkel Rücken im Arktischen Ozean von großer Bedeutung für das Verständnis von Aufschmelzprozessen, da er weltweit die niedrigste Divergenzrate aufweist und große Transformstörungen fehlen. Publizierte Basaltdaten deuten auf einen extrem niedrigen Aufschmelzgrad hin, was mit globalen Korrelationen im Einklang steht. Stark alterierte Mantelperidotite einer Lokalität entlang des kaum beprobten Gakkel Rückens wurden deshalb auf Primärminerale untersucht. Nur in einer Probe sind oxidierte Spinellpseudomorphosen mit Spuren primärer Spinelle erhalten geblieben. Ihre Cr# ist signifikant höher als die einiger Peridotite von schneller divergierenden Rücken und ihr Schmelzgrad ist damit höher als aufgrund der Basaltzusammensetzungen vermutet. Der unter Anwendung des oben erwähnten Indikators ermittelte Schmelzgrad ermöglicht die Berechnung der Krustenmächtigkeit am Gakkel Rücken. Diese ist wesentlich größer als die aus Schweredaten ermittelte Mächtigkeit, oder die aus der globalen Korrelation zwischen Divergenzrate und mittels Seismik erhaltene Krustendicke. Dieses unerwartete Ergebnis kann möglicherweise auf kompositionelle Heterogenitäten bei niedrigen Schmelzgraden, oder auf eine insgesamt größere Verarmung des Mantels unter dem Gakkel Rücken zurückgeführt werden (Hellebrand et al., Chem.Geol., in review; Kapitel 3).Zusätzliche Informationen zur Modellierung und Analytik sind im Anhang A-C aufgeführt

Geological and Structural evolution of the Eurasia Africa plate boundary in the Gulf of Cadiz Central Eastern Atlantic Sea.

Iberia Africa plate boundary, cross, roughly W-E, connecting the eastern Atlantic Ocean from Azores triple junction to the Continental margin of Morocco. Relative movement between the two plate change along the boundary, from transtensive near the Azores archipelago, through trascurrent movement in the middle at the Gloria Fracture Zone, to transpressive in the Gulf of Cadiz area. This study presents the results of geophysical and geological analysis on the plate boundary area offshore Gibraltar. The main topic is to clarify the geodynamic evolution of this area from Oligocene to Quaternary. Recent studies have shown that the new plate boundary is represented by a 600 km long set of aligned, dextral trascurrent faults (the SWIM lineaments) connecting the Gloria fault to the Riff orogene. The western termination of these lineaments crosscuts the Gibraltar accretionary prism and seems to reach the Moroccan continental shelf. In the past two years newly acquired bathymetric data collected in the Moroccan offshore permit to enlighten the present position of the eastern portion of the plate boundary, previously thought to be a diffuse plate boundary. The plate boundary evolution, from the onset of compression in the Oligocene to the Late Pliocene activation of trascurrent structures, is not yet well constrained. The review of available seismics lines, gravity and bathymetric data, together with the analysis of new acquired bathymetric and high resolution seismic data offshore Morocco, allows to understand how the deformation acted at lithospheric scale under the compressive regime. Lithospheric folding in the area is suggested, and a new conceptual model is proposed for the propagation of the deformation acting in the brittle crust during this process. Our results show that lithospheric folding, both in oceanic and thinned continental crust, produced large wavelength synclines bounded by short wavelength, top thrust, anticlines. Two of these anticlines are located in the Gulf of Cadiz, and are represented by the Gorringe Ridge and Coral Patch seamounts. Lithospheric folding probably interacted with the Monchique – Madeira hotspot during the 72 Ma to Recent, NNE – SSW transit. Plume related volcanism is for the first time described on top of the Coral Patch seamount, where nine volcanoes are found by means of bathymetric data. 40Ar-39Ar age of 31.4±1.98 Ma are measured from one rock sample of one of these volcanoes. Analysis on biogenic samples show how the Coral Patch act as a starved offshore seamount since the Chattian. We proposed that compression stress formed lithospheric scale structures playing as a reserved lane for the upwelling of mantle material during the hotspot transit. The interaction between lithospheric folding and the hotspot emplacement can be also responsible for the irregularly spacing, and anomalous alignments, of individual islands and seamounts belonging to the Monchique - Madeira hotspot.

Geochronology, geochemistry and isotopic compositions of the volcanic rocks from oceanic (Hawaii) and continental (Eifel) intra-plate environments

The Eifel volcanism is part of the Central European Volcanic Province (CEVP) and is located in the Rhenish Massif, close to the Rhine and Leine Grabens. The Quaternary Eifel volcanism appears to be related to a mantle plume activity. However, the causes of the Tertiary Hocheifel volcanism remain debated. We present geochronological, geochemical and isotope data to assess the geotectonic settings of the Tertiary Eifel volcanism. Based on 40Ar/39Ar dating, we were able to identify two periods in the Hocheifel activity: from 43.6 to 39.0 Ma and from 37.5 to 35.0 Ma. We also show that the pre-rifting volcanism in the northernmost Upper Rhine Graben (59 to 47 Ma) closely precede the Hocheifel volcanic activity. In addition, the volcanism propagates from south to north within the older phase of the Hocheifel activity. At the time of Hocheifel volcanism, the tectonic activity in the Hocheifel was controlled by stress field conditions identical to those of the Upper Rhine Graben. Therefore, magma generation in the Hocheifel appears to be caused by decompression due to Middle to Late Eocene extension. Our geochemical data indicate that the Hocheifel magmas were produced by partial melting of a garnet peridotite at 75-90 km depth. We also show that crustal contamination is minor although the magmas erupted through a relatively thick continental lithosphere. Sr, Nd and Pb isotopic compositions suggest that the source of the Hocheifel magmas is a mixing between depleted FOZO or HIMU-like material and enriched EM2-like material. The Tertiary Hocheifel and the Quaternary Eifel lavas appear to have a common enriched end-member. However, the other sources are likely to be distinct. In addition, the Hocheifel lavas share a depleted component with the other Tertiary CEVP lavas. Although the Tertiary Hocheifel and the Quaternary Eifel lavas appear to originate from different sources, the potential involvement of a FOZO-like component would indicate the contribution of deep mantle material. Thus, on the basis of the geochemical and isotope data, we cannot rule out the involvement of plume-type material in the Hocheifel magmas. The Ko’olau Scientific Drilling Project (KSDP) has been initiated in order to evaluate the long-term evolution of Ko’olau volcano and obtain information about the Hawaiian mantle plume. High precision Pb triple spike data, as well as Sr and Nd isotope data on KSDP lavas and Honolulu Volcanics (HVS) reveal compositional source variations during Ko’olau growth. Pb isotopic compositions indicate that, at least, three Pb end-members are present in Ko’olau lavas. Changes in the contributions of each component are recorded in the Pb, Sr and Nd isotopes stratigraphy. The radiogenic component is present, at variable proportion, in all three stages of Ko’olau growth. It shows affinities with the least radiogenic “Kea-lo8” lavas present in Mauna Kea. The first unradiogenic component was present in the main-shield stage of Ko’olau growth but its contribution decreased with time. It has EM1 type characteristics and corresponds to the “Ko’olau” component of Hawaiian mantle plume. The second unradiogenic end-member, so far only sampled by Honololu lavas, has isotopic characteristics similar to those of a depleted mantle. However, they are different from those of the recent Pacific lithosphere (EPR MORB) indicating that the HVS are not derived from MORB-related source. We suggest, instead, that the HVS result from melting of a plume material. Thus the evolution of a single Hawaiian volcano records the geochemical and isotopic changes within the Hawaiian plume.

Numerical models of trench migration for lateral heterogeneous subducting plates

The aim of this Thesis is to investigate the effect of heterogeneities within the subducting plate on the dynamics of subduction. In particular, I study the motion of the trench for oceanic and continental subduction, first, separately, and, then, together in the same system to understand how they interact. The understanding of these features is fundamental to reconstruct the evolution of complex subduction zones, such as the Central Mediterranean. For this purpose, I developed 2D and 3D numerical models of oceanic and continental subduction where the rheological, geometrical and compositional properties of the plates are varied. In these models, the trench and the overriding plate move self-consistently as a function of the dynamics of the system. The effect of continental subduction on trench migration is largely investigated. Results from a parametric study showed that despite different rheological properties of the plates, all models with a uniform continental crust share the same kinematic behaviour: the trench starts to advance once the continent arrives at the subduction zone. Hence, the advancing mode in continental collision scenarios is at least partly driven by an intrinsic feature of the system. Moreover, the presence of a weak lower crust within the continental plate can lead to the occurrence of delamination. Indeed, by changing the viscosity of the lower crust, both delamination and slab detachment can occur. Delamination is favoured by a low viscosity value of the lower crust, because this makes the mechanical decoupling easier between crust and lithospheric mantle. These features are observed both in 2D and 3D models, but the numerical results of the 3D models also showed that the rheology of the continental crust has a very strong effect on the dynamics of the whole system, since it influences not only the continental part of plate but also the oceanic sides.

Partial melting within the lower crust – Constraints from bulk rock geochemical data and trace element mineral analyses of granulites and migmatites from Central Finland

This PhD thesis concerns geochemical constraints on recycling and partial melting of Archean continental crust. A natural example of such processes was found in the Iisalmi area of Central Finland. The rocks from this area are Middle to Late Archean in age and experienced metamorphism and partial melting between 2.7-2.63 Ga. The work is based on extensive field work. It is furthermore founded on bulk rock geochemical data as well as in-situ analyses of minerals. All geochemical data were obtained at the Institute of Geosciences, University of Mainz using X-ray fluorescence, solution ICP-MS and laser ablation-ICP-MS for bulk rock geochemical analyses. Mineral analyses were accomplished by electron microprobe and laser ablation ICP-MS. Fluid inclusions were studied by microscope on a heating-freezing-stage at the Geoscience Center, University Göttingen. Part I focuses on the development of a new analytical method for bulk rock trace element determination by laser ablation-ICP-MS using homogeneous glasses fused from rock powder on an Iridium strip heater. This method is applicable for mafic rock samples whose melts have low viscosities and homogenize quickly at temperatures of ~1200°C. Highly viscous melts of felsic samples prevent melting and homogenization at comparable temperatures. Fusion of felsic samples can be enabled by addition of MgO to the rock powder and adjustment of melting temperature and melting duration to the rock composition. Advantages of the fusion method are low detection limits compared to XRF analyses and avoidance of wet-chemical processing and use of strong acids as in solution ICP-MS as well as smaller sample volumes compared to the other methods. Part II of the thesis uses bulk rock geochemical data and results from fluid inclusion studies for discrimination of melting processes observed in different rock types. Fluid inclusion studies demonstrate a major change in fluid composition from CO2-dominated fluids in granulites to aqueous fluids in TTG gneisses and amphibolites. Partial melts were generated in the dry, CO2-rich environment by dehydration melting reactions of amphibole which in addition to tonalitic melts produced the anhydrous mineral assemblages of granulites (grt + cpx + pl ± amph or opx + cpx + pl + amph). Trace element modeling showed that mafic granulites are residues of 10-30 % melt extraction from amphibolitic precursor rocks. The maximum degree of melting in intermediate granulites was ~10 % as inferred from modal abundances of amphibole, clinopyroxene and orthopyroxene. Carbonic inclusions are absent in upper-amphibolite facies migmatites whereas aqueous inclusion with up to 20 wt% NaCl are abundant. This suggests that melting within TTG gneisses and amphibolites took place in the presence of an aqueous fluid phase that enabled melting at the wet solidus at temperatures of 700-750°C. The strong disruption of pre-metamorphic structures in some outcrops suggests that the maximum amount of melt in TTG gneisses was ~25 vol%. The presence of leucosomes in all rock types is taken as the principle evidence for melt formation. However, mineralogical appearance as well as major and trace element composition of many leucosomes imply that leucosomes seldom represent frozen in-situ melts. They are better considered as remnants of the melt channel network, e.g. ways on which melts escaped from the system. Part III of the thesis describes how analyses of minerals from a specific rock type (granulite) can be used to determine partition coefficients between different minerals and between minerals and melt suitable for lower crustal conditions. The trace element analyses by laser ablation-ICP-MS show coherent distribution among the principal mineral phases independent of rock composition. REE contents in amphibole are about 3 times higher than REE contents in clinopyroxene from the same sample. This consistency has to be taken into consideration in models of lower crustal melting where amphibole is replaced by clinopyroxene in the course of melting. A lack of equilibrium is observed between matrix clinopyroxene / amphibole and garnet porphyroblasts which suggests a late stage growth of garnet and slow diffusion and equilibration of the REE during metamorphism. The data provide a first set of distribution coefficients of the transition metals (Sc, V, Cr, Ni) in the lower crust. In addition, analyses of ilmenite and apatite demonstrate the strong influence of accessory phases on trace element distribution. Apatite contains high amounts of REE and Sr while ilmenite incorporates about 20-30 times higher amounts of Nb and Ta than amphibole. Furthermore, trace element mineral analyses provide evidence for magmatic processes such as melt depletion, melt segregation, accumulation and fractionation as well as metasomatism having operated in this high-grade anatectic area.

Geochemistry, petrogenesis and tectonic setting of ophiolites and mafic-ultramafic complexes in the Northeastern Aegean Region

In this study two ophiolites and a mafic-ultramafic complexes of the northeastern Aegean Sea, Greece, have been investigated to re-evaluate their petrogenetic evolution and tectonic setting. These complexes are: the mafic-ultramafic complex of Lesvos Island and the ophiolites of Samothraki Island and the Evros area. In order to examine these complexes in detail whole-rock major- and trace-elements as well as Sr and Nd isotopes, and minerals were analysed and U-Pb SHRIMP ages on zircons were determined. The mafic-ultramafic complex of Lesvos Island consists of mantle peridotite thrusted over a tectonic mélange containing metasediments, metabasalts and a few metagabbros. This succession had previously been interpreted as an ophiolite of Late Jurassic age. The new field and geochemical data allow a reinterpretation of this complex as representing an incipient continental rift setting that led to the subsequent formation of the Meliata-Maliac-Vardar branches of Neotethys in Upper Permian times (253 ± 6 Ma) and the term “Lesvos ophiolite” should be abandoned. With proceeding subduction and closure of the Maliac Ocean in Late Jurassic times (155 Ma) the Lesvos mafic-ultramafic complex was obducted. Zircon ages of 777, 539 and 338 Ma from a gabbro strongly suggest inheritance from the intruded basement and correspond to ages of distinct terranes recently recognized in the Hellenides (e.g. Florina terrane). Geochemical similar complexes which contain rift associations with Permo-Triassic ages can be found elsewhere in Greece and Turkey, namely the Teke Dere Thrust Sheet below the Lycian Nappes (SW Turkey), the Pindos subophiolitic mélange (W Greece), the Volcanosedimentary Complex on Central Evia Island (Greece) and the Karakaya Complex (NW Turkey). This infers that the rift-related rocks from Lesvos belong to an important Permo-Triassic rifting episode in the eastern Mediterranean. The ‘in-situ’ ophiolite of Samothraki Island comprises gabbros, sparse dykes and basalt flows as well as pillows cut by late dolerite dykes and had conventionally been interpreted as having formed in an ensialic back-arc basin. The results of this study revealed that none of the basalts and dolerites resemble mid-ocean ridge or back-arc basin basalts thus suggesting that the Samothraki ophiolite cannot represent mature back-arc basin crust. The age of the complex is regarded to be 160 ± 5 Ma (i.e. Oxfordian; early Upper Jurassic), which precludes any correlation with the Lesvos mafic-ultramafic complex further south (253 ± 6 Ma; Upper Permian). Restoration of the block configuration in NE Greece, before extensional collapse of the Hellenic hinterland and exhumation of the Rhodope Metamorphic Core Complex (mid-Eocene to mid-Miocene), results in a continuous ophiolite belt from Guevgueli in the NW to Samothraki in the SE, thus assigning the latter to the Innermost Hellenic Ophiolite Belt. In view of the data of this study, the Samothraki ophiolite represents a rift propagation of the Sithonia ophiolite spreading ridge into the Chortiatis calc-alkaline arc. The ophiolite of the Evros area consists of a plutonic sequence comprising cumulate and non-cumulate gabbros with plagiogranite veins, and an extrusive sequence of basalt dykes, massive and pillow lavas as well as pyroclastic rocks. Furthermore, in the Rhodope Massif tectonic lenses of harzburgites and dunites can be found. All rocks are spatially separated. The analytical results of this study revealed an intra-oceanic island arc setting for the Evros ophiolitic rocks. During late Middle Jurassic times (169 ± 2 Ma) an intra-oceanic arc has developed above a northwards directed intra-oceanic subduction zone of the Vardar Ocean in front of the Rhodope Massif. The boninitic, island arc tholeiitic and calc-alkaline rocks reflect the evolution of the Evros island arc. The obduction of the ophiolitic rocks onto the Rhodope basement margin took place during closure of the Vardar ocean basins. The harzburgites and dunites of the Rhodope Massif are strongly depleted and resemble harzburgites from recent oceanic island arcs. After melt extraction they underwent enrichment processes by percolating melts and fluids from the subducted slab. The relationship of the peridotites and the Evros ophiolite is still ambiguous, but the stratigraphic positions of the peridotites and the ophiolitic rocks indicate separated origin. The harzburgites and dunites most probably represent remnants of the mantle wedge of the island arc of the Rhodope terrane formed above subducted slab of the Nestos Ocean in late Middle Jurassic times. During collision of the Thracia terrane with the Rhodope terrane thrusting of the Rhodope terrane onto the Thracia terrane took place, whereas the harzburgites and dunites were pushed between the two terranes now cropping out on top of the Thracia terrane of the Rhodope Massif.

Transform Tectonics and Non-Volcanic Oceanic Islands

Oceanic islands can be divided, according to their origin, in volcanic and tectonic. Volcanic islands are due to excess volcanism. Tectonic islands are mainly formed due to vertical tectonic motions of blocks of oceanic lithosphere along transverse ridges flanking transform faults at slow and ultraslow mid-ocean ridges. Vertical tectonic motions are due to a reorganization of the geometry of the transform plate boundary, with the transition from a transcurrent tectonics to a transtensive and/or transpressive tectonics, with the formation of the transverse ridges. Tectonic islands can be located also at the ridge–transform intersection: in this case the uplift is due by the movement of the long-lived detachment faults located along the flanks of the mid-ocean ridges. The "Vema" paleoisland (equatorial Atlantic) is at the summit of the southern transverse ridge of the Vema transform. It is now 450 m bsl and it is capped by a carbonate platform 500 m-thick, dated by 87Sr/86Sr at 10 Ma. Three tectonic paleoislands are on the summit of the transverse ridge flanking the Romanche megatrasform (equatorial Atlantic). They are now about 1,000 m bsl and they are formed by 300 m-thick carbonate platforms dated by 87Sr/86Sr, between 11 and 6 Ma. The tectonic paleoisland “Atlantis Bank" is located in the South-Western Indian Ridge, along the Atlantis II transform, and it is today 700 m bsl. The only modern example of oceanic tectonics island is the St. Paul Rocks (equatorial Atlantic), located along the St. Paul transform. This archipelago is the top of a peridotitic massif that it is now a left overstep undergoing transpression. Oceanic volcanic islands are characterized by rapid growth and subsequent thermal subsidence and drowning; in contrast, oceanic tectonic islands may have one or more stages of emersion related to vertical tectonic events along the large oceanic fracture zones.

Volatile organic compounds at the air-sea interface : gas exchange rates, oceanic emissions and the effect of ocean acidification

Oceans are key sources and sinks in the global budgets of significant atmospheric trace gases, termed Volatile Organic Compounds (VOCs). Despite their low concentrations, these species have an important role in the atmosphere, influencing ozone photochemistry and aerosol physics. Surprisingly, little work has been done on assessing their emissions or transport mechanisms and rates between ocean and atmosphere, all of which are important when modelling the atmosphere accurately.rnA new Needle Trap Device (NTD) - GC-MS method was developed for the effective sampling and analysis of VOCs in seawater. Good repeatability (RSDs <16 %), linearity (R2 = 0.96 - 0.99) and limits of detection in the range of pM were obtained for DMS, isoprene, benzene, toluene, p-xylene, (+)-α-pinene and (-)-α-pinene. Laboratory evaluation and subsequent field application indicated that the proposed method can be used successfully in place of the more usually applied extraction techniques (P&T, SPME) to extend the suite of species typically measured in the ocean and improve detection limits. rnDuring a mesocosm CO2 enrichment study, DMS, isoprene and α-pinene were identified and quantified in seawater samples, using the above mentioned method. Based on correlations with available biological datasets, the effects of ocean acidification as well as possible ocean biological sources were investigated for all examined compounds. Future ocean's acidity was shown to decrease oceanic DMS production, possibly impact isoprene emissions but not affect the production of α-pinene. rnIn a separate activity, ocean - atmosphere interactions were simulated in a large scale wind-wave canal facility, in order to investigate the gas exchange process and its controlling mechanisms. Air-water exchange rates of 14 chemical species (of which 11 VOCs) spanning a wide range of solubility (dimensionless solubility, α = 0:4 to 5470) and diffusivity (Schmidt number in water, Scw = 594 to 1194) were obtained under various turbulent (wind speed at ten meters height, u10 = 0:8 to 15ms-1) and surfactant modulated (two different sized Triton X-100 layers) surface conditions. Reliable and reproducible total gas transfer velocities were obtained and the derived values and trends were comparable to previous investigations. Through this study, a much better and more comprehensive understanding of the gas exchange process was accomplished. The role of friction velocity, uw* and mean square slope, σs2 in defining phenomena such as waves and wave breaking, near surface turbulence, bubbles and surface films was recognized as very significant. uw* was determined as the ideal turbulent parameter while σs2 described best the related surface conditions. A combination of both uw* and σs2 variables, was found to reproduce faithfully the air-water gas exchange process. rnA Total Transfer Velocity (TTV) model provided by a compilation of 14 tracers and a combination of both uw* and σs2 parameters, is proposed for the first time. Through the proposed TTV parameterization, a new physical perspective is presented which provides an accurate TTV for any tracer within the examined solubility range. rnThe development of such a comprehensive air-sea gas exchange parameterization represents a highly useful tool for regional and global models, providing accurate total transfer velocity estimations for any tracer and any sea-surface status, simplifying the calculation process and eliminating inevitable calculation uncertainty connected with the selection or combination of different parameterizations.rnrn